TY - JOUR
T1 - Multiple timescale molecular dynamics with very large time steps
T2 - avoidance of resonances
AU - Abreu, C. R.A.
AU - Tuckerman, M. E.
N1 - Funding Information:
C.R.A.A. thanks the Federal University of Rio de Janeiro for a sabbatical leave and the New York University for a visiting scholar program membership. M.E.T. acknowledges funding from the National Science Foundation through grant no. CHE-1955381.
Funding Information:
C.R.A.A. thanks the Federal University of Rio de Janeiro for a sabbatical leave and the New York University for a visiting scholar program membership. M.E.T. acknowledges funding from the National Science Foundation through grant no. CHE-1955381.
Publisher Copyright:
© 2021, The Author(s), under exclusive licence to EDP Sciences, SIF and Springer-Verlag GmbH Germany, part of Springer Nature.
PY - 2021/11
Y1 - 2021/11
N2 - Reversible multiple timescale (MTS) integration algorithms have long been recognized as a straightforward way to increase efficiency and extend accessible timescales in molecular dynamics simulations without altering the ensemble distribution sampled. MTS methods are based on the idea that interatomic forces in a system drive motion on numerous timescales, and by decomposing force components according to these timescales and assigning an individual time step to each one, fast, computationally cheaper forces are evaluated more frequently than the slow, expensive forces. As it happens, the largest time step that can be employed in standard MTS methods is fundamentally limited by so-called resonance artifacts that originate in the fastest timescales. Thus, while it should be possible to assign the slowest timescales very large time steps approaching 100 fs in, for example, fully atomistic simulations, resonances impose a practical limit on this step size to around 5–10 fs, which allows for useful but only modest savings in computational overhead. This article will review the basic MTS approach and the origin of resonances and then will provide a perspective on how to solve the resonance problem for molecular dynamics simulations in different ensembles, showing how both statistical and dynamical properties can be generated with very large time steps.
AB - Reversible multiple timescale (MTS) integration algorithms have long been recognized as a straightforward way to increase efficiency and extend accessible timescales in molecular dynamics simulations without altering the ensemble distribution sampled. MTS methods are based on the idea that interatomic forces in a system drive motion on numerous timescales, and by decomposing force components according to these timescales and assigning an individual time step to each one, fast, computationally cheaper forces are evaluated more frequently than the slow, expensive forces. As it happens, the largest time step that can be employed in standard MTS methods is fundamentally limited by so-called resonance artifacts that originate in the fastest timescales. Thus, while it should be possible to assign the slowest timescales very large time steps approaching 100 fs in, for example, fully atomistic simulations, resonances impose a practical limit on this step size to around 5–10 fs, which allows for useful but only modest savings in computational overhead. This article will review the basic MTS approach and the origin of resonances and then will provide a perspective on how to solve the resonance problem for molecular dynamics simulations in different ensembles, showing how both statistical and dynamical properties can be generated with very large time steps.
UR - http://www.scopus.com/inward/record.url?scp=85120311265&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=85120311265&partnerID=8YFLogxK
U2 - 10.1140/epjb/s10051-021-00226-4
DO - 10.1140/epjb/s10051-021-00226-4
M3 - Review article
AN - SCOPUS:85120311265
SN - 1434-6028
VL - 94
JO - European Physical Journal B
JF - European Physical Journal B
IS - 11
M1 - 231
ER -