Synchronization of trajectories in canonical molecular-dynamics simulations: observation, explanation, and exploitation

For two methods commonly used to achieve canonical-ensemble sampling in a molecular-dynamics simulation, the Langevin thermostat and the Andersen [H. C. Andersen, J. Chem. Phys. 72, 2384 (1980)] thermostat, we observe, as have others, synchronization of initially independent trajectories in the same...

Full description

Saved in:
Bibliographic Details
Published in:The Journal of chemical physics 2004-04, Vol.120 (14), p.6363-6374
Main Authors: Uberuaga, Blas P, Anghel, Marian, Voter, Arthur F
Format: Article
Language:English
Subjects:
Algorithms
Computer Simulation
Molecular Conformation
Thermodynamics
Citations: Items that this one cites
Items that cite this one
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:For two methods commonly used to achieve canonical-ensemble sampling in a molecular-dynamics simulation, the Langevin thermostat and the Andersen [H. C. Andersen, J. Chem. Phys. 72, 2384 (1980)] thermostat, we observe, as have others, synchronization of initially independent trajectories in the same potential basin when the same random number sequence is employed. For the first time, we derive the time dependence of this synchronization for a harmonic well and show that the rate of synchronization is proportional to the thermostat coupling strength at weak coupling and inversely proportional at strong coupling with a peak in between. Explanations for the synchronization and the coupling dependence are given for both thermostats. Observation of the effect for a realistic 97-atom system indicates that this phenomenon is quite general. We discuss some of the implications of this effect and propose that it can be exploited to develop new simulation techniques. We give three examples: efficient thermalization (a concept which was also noted by Fahy and Hamann [S. Fahy and D. R. Hamann, Phys. Rev. Lett. 69, 761 (1992)]), time-parallelization of a trajectory in an infrequent-event system, and detecting transitions in an infrequent-event system.
ISSN:0021-9606
1089-7690
DOI:10.1063/1.1667473