Molecular Dynamics Simulations of Solutions at Constant Chemical Potential

Abstract

Molecular Dynamics studies of chemical processes in solution are of great value in a wide spectrum of applications, which range from nano-technology to pharmaceutical chemistry. However, these calculations are affected by severe finite-size effects, such as the solution being depleted as the chemical process proceeds, which influence the outcome of the simulations. To overcome these limitations, one must allow the system to exchange molecules with a macroscopic reservoir, thus sampling a Grand-Canonical ensemble. Despite the fact that different remedies have been proposed, this still represents a key challenge in molecular simulations. In the present work we propose the Constant Chemical Potential Molecular Dynamics (C$\mu$MD) method, which introduces an external force that controls the environment of the chemical process of interest. This external force, drawing molecules from a finite reservoir, maintains the chemical potential constant in the region where the process takes place. We have applied the C$\mu$MD method to the paradigmatic case of urea crystallization in aqueous solution. As a result, we have been able to study crystal growth dynamics under constant supersaturation conditions, and to extract growth rates and free-energy barriers.

8 pages, 8 figures (Supplementary Information: 6 pages, 7 figures). Typos and labelling corrected Ver. 3: Minor comments added in Sec. 3. References 13,36,38 added. Minor text changes and typos corrected

References in corpus (6)

• Canonical sampling through velocity-rescaling
• Well-Tempered Metadynamics: A Smoothly Converging and Tunable Free-Energy Method
• PLUMED: a portable plugin for free-energy calculations with molecular dynamics
• Hamiltonian adaptive resolution simulation for molecular liquids
• Molecular Dynamics Simulations of Solutions at Constant Chemical Potential
• Grand-canonical-like molecular-dynamics simulations by using an adaptive-resolution technique

Cited by in corpus (4)

• Crystal Nucleation in Liquids: Open Questions and Future Challenges in Molecular Dynamics Simulations
• Molecular Dynamics Simulations of Solutions at Constant Chemical Potential
• Molecular dynamics of open systems: construction of a mean-field particle reservoir
• Simulations of Crystal Nucleation from Solution at Constant Chemical Potential