Perego, C; Salvalaglio, M; Parrinello, M; (2015) Molecular dynamics simulations of solutions at constant chemical potential. The Journal of Chemical Physics , 142 (14) , Article 144113. 10.1063/1.4917200.

Preview

Text Salvalaglio_Molecular_dynamics_JCP_CmuMD_2015.pdf - Published Version Download (3MB) | Preview

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μ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μ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.

Download activity - last month

Export as JSONXMLCSV

Download activity - last 12 months

Export as JSONXMLCSV

Downloads by country - last 12 months

Export as CSVXMLJSON

Archive Staff Only

View Item