Efficiency in nonequilibrium molecular dynamics Monte Carlo simulations
- Argonne National Lab. (ANL), Argonne, IL (United States)
- Univ. of Chicago, IL (United States). Dept. of Biochemistry and Molecular Biology; Argonne National Lab. (ANL), Argonne, IL (United States). Center for Nanoscale Materials
Hybrid algorithms combining nonequilibrium molecular dynamics and Monte Carlo (neMD/MC) offer a powerful avenue for improving the sampling efficiency of computer simulations of complex systems. These neMD/MC algorithms are also increasingly finding use in applications where conventional approaches are impractical, such as constant-pH simulations with explicit solvent. However, selecting an optimal nonequilibrium protocol for maximum efficiency often represents a non-trivial challenge. This work evaluates the efficiency of a broad class of neMD/MC algorithms and protocols within the theoretical framework of linear response theory. The approximations are validated against constant pH-MD simulations and shown to provide accurate predictions of neMD/MC performance. An assessment of a large set of protocols confirms (both theoretically and empirically) that a linear work protocol gives the best neMD/MC performance. Lastly, a well-defined criterion for optimizing the time parameters of the protocol is proposed and demonstrated with an adaptive algorithm that improves the performance on-the-fly with minimal cost.
- Research Organization:
- Argonne National Laboratory (ANL), Argonne, IL (United States)
- Sponsoring Organization:
- USDOE Office of Science (SC), Basic Energy Sciences (BES); National Institutes of Health (NIH)
- Grant/Contract Number:
- AC02-06CH11357; U54-GM087519
- OSTI ID:
- 1374573
- Journal Information:
- Journal of Chemical Physics, Vol. 145, Issue 13; ISSN 0021-9606
- Publisher:
- American Institute of Physics (AIP)Copyright Statement
- Country of Publication:
- United States
- Language:
- English
Web of Science
Similar Records
Enhanced configurational sampling with hybrid non-equilibrium molecular dynamics–Monte Carlo propagator
Generalized Metropolis acceptance criterion for hybrid non-equilibrium molecular dynamics—Monte Carlo simulations