Grand-canonical molecular dynamics simulations powered by a hybrid 4D nonequilibrium MD/MC method: Implementation in LAMMPS and applications to electrolyte solutions - Sorbonne Université Accéder directement au contenu
Article Dans Une Revue The Journal of Chemical Physics Année : 2023

Grand-canonical molecular dynamics simulations powered by a hybrid 4D nonequilibrium MD/MC method: Implementation in LAMMPS and applications to electrolyte solutions

Résumé

Molecular simulations in an open environment, involving ion exchange, are necessary to study various systems, from biosystems to confined electrolytes. However, grand-canonical simulations are often computationally demanding in condensed phases. A promising method [L. Belloni, J. Chem. Phys. 151, 021101 (2019)], one of the hybrid nonequilibrium molecular dynamics/Monte Carlo algorithms, was recently developed, which enables efficient computation of fluctuating number or charge density in dense fluids or ionic solutions. This method facilitates the exchange through an auxiliary dimension, orthogonal to all physical dimensions, by reducing initial steric and electrostatic clashes in three-dimensional systems. Here, we report the implementation of the method in LAMMPS with a Python interface, allowing facile access to grand-canonical molecular dynamics simulations with massively parallelized computation. We validate our implementation with two electrolytes, including a model Lennard-Jones electrolyte similar to a restricted primitive model and aqueous solutions. We find that electrostatic interactions play a crucial role in the overall efficiency due to their long-range nature, particularly for water or ion-pair exchange in aqueous solutions. With properly screened electrostatic interactions and bias-based methods, our approach enhances the efficiency of salt-pair exchange in Lennard-Jones electrolytes by approximately four orders of magnitude, compared to conventional grand-canonical Monte Carlo. Furthermore, the acceptance rate of NaCl-pair exchange in aqueous solutions at moderate concentrations reaches about 3% at the maximum efficiency.
Fichier principal
Vignette du fichier
JCP23-AR-02525.pdf (659.2 Ko) Télécharger le fichier
Origine Fichiers produits par l'(les) auteur(s)

Dates et versions

hal-04237391 , version 1 (11-10-2023)

Identifiants

Citer

Jeongmin Kim, Luc Belloni, Benjamin Rotenberg. Grand-canonical molecular dynamics simulations powered by a hybrid 4D nonequilibrium MD/MC method: Implementation in LAMMPS and applications to electrolyte solutions. The Journal of Chemical Physics, 2023, 159 (14), pp.144802. ⟨10.1063/5.0168878⟩. ⟨hal-04237391⟩
42 Consultations
25 Téléchargements

Altmetric

Partager

Gmail Mastodon Facebook X LinkedIn More