Enhanced configurational sampling with hybrid non-equilibrium molecular dynamics–Monte Carlo propagator
[1];
[2];
- Univ. of Chicago, IL (United States)
- Argonne National Lab. (ANL), Argonne, IL (United States). Leadership Computing Facility
- Univ. de Lorraine (France); Univ. of Illinois at Urbana-Champaign, IL (United States)
- Univ. of Chicago, IL (United States); Argonne National Lab. (ANL), Argonne, IL (United States). Center for Nanoscale Materials
Molecular dynamics (MD) trajectories based on classical equations of motion can be used to sample the configurational space of complex molecular systems. However, brute-force MD often converges slowly due to the ruggedness of the underlying potential energy surface. Several schemes have been proposed to address this problem by effectively smoothing the potential energy surface. However, in order to recover the proper Boltzmann equilibrium probability distribution, these approaches must then rely on statistical reweighting techniques or generate the simulations within a Hamiltonian tempering replica-exchange scheme. The present work puts forth a novel hybrid sampling propagator combining Metropolis-Hastings Monte Carlo (MC) with proposed moves generated by non-equilibrium MD (neMD). This hybrid neMD-MC propagator comprises three elementary elements: (i) an atomic system is dynamically propagated for some period of time using standard equilibrium MD on the correct potential energy surface; (ii) the system is then propagated for a brief period of time during what is referred to as a “boosting phase,” via a time-dependent Hamiltonian that is evolved toward the perturbed potential energy surface and then back to the correct potential energy surface; (iii) the resulting configuration at the end of the neMD trajectory is then accepted or rejected according to a Metropolis criterion before returning to step 1. A symmetric two-end momentum reversal prescription is used at the end of the neMD trajectories to guarantee that the hybrid neMD-MC sampling propagator obeys microscopic detailed balance and rigorously yields the equilibrium Boltzmann distribution. The hybrid neMD-MC sampling propagator is designed and implemented to enhance the sampling by relying on the accelerated MD and solute tempering schemes. It is also combined with the adaptive biased force sampling algorithm to examine. Illustrative tests with specific biomolecular systems indicate that the method can yield a significant speedup.
- Research Organization:
- Argonne National Lab. (ANL), Argonne, IL (United States)
- Sponsoring Organization:
- National Science Foundation (NSF); USDOE Office of Science (SC), Basic Energy Sciences (BES)
- Grant/Contract Number:
- AC02-06CH11357; MCB-1517221
- OSTI ID:
- 1490226
- Alternate ID(s):
- OSTI ID: 1415336
- Journal Information:
- Journal of Chemical Physics, Vol. 148, Issue 1; ISSN 0021-9606
- Publisher:
- American Institute of Physics (AIP)Copyright Statement
- Country of Publication:
- United States
- Language:
- English
Web of Science
A Computational Predictive Approach for Controlling the Morphology of Functional Molecular Aggregates on Substrates
|
journal | September 2019 |
The Mechanism of Cholesterol Modification of Hedgehog Ligand
|
journal | November 2019 |
Thermodynamics of helix formation in small peptides of varying length in vacuo, in implicit solvent, and in explicit solvent
|
journal | December 2018 |
Similar Records
SU-E-T-58: A Novel Monte Carlo Photon Transport Simulation Scheme and Its Application in Cone Beam CT Projection Simulation
Dimension-free path-integral molecular dynamics without preconditioning