2017
DOI: 10.1021/acs.jctc.7b00178
|View full text |Cite
|
Sign up to set email alerts
|

Real Cost of Speed: The Effect of a Time-Saving Multiple-Time-Stepping Algorithm on the Accuracy of Molecular Dynamics Simulations

Abstract: To enhance efficiency in molecular dynamics simulations, forces that vary slowly are often evaluated less often than those that vary rapidly. We show that the multiple-time-step algorithm implemented in recent versions of GROMACS results in significant differences in the collective properties of a system under conditions where the system was previously stable. The implications of changing the simulation algorithm without assessment of potential artifacts on the parametrization and transferability of effective … Show more

Help me understand this report

Search citation statements

Order By: Relevance

Paper Sections

Select...
2
2
1

Citation Types

4
74
2

Year Published

2018
2018
2021
2021

Publication Types

Select...
7
1
1

Relationship

2
7

Authors

Journals

citations
Cited by 48 publications
(80 citation statements)
references
References 47 publications
4
74
2
Order By: Relevance
“…Previous work with the GROMOS force field used a twin-range cutoff scheme but its association with the current implementation of the integrator in GROMACS leads to artefacts in the collective properties of a lipid bilayer. 74 Despite some loss in computational efficiency, using a single-range cutoff scheme was found to counterbalance those artefacts. 74 To correct for the truncation of electrostatic interactions beyond 1.4 nm a reactionfield correction 75 with a relative dielectric constant (e r ) of 62 was applied.…”
Section: Simulation Parametersmentioning
confidence: 99%
See 1 more Smart Citation
“…Previous work with the GROMOS force field used a twin-range cutoff scheme but its association with the current implementation of the integrator in GROMACS leads to artefacts in the collective properties of a lipid bilayer. 74 Despite some loss in computational efficiency, using a single-range cutoff scheme was found to counterbalance those artefacts. 74 To correct for the truncation of electrostatic interactions beyond 1.4 nm a reactionfield correction 75 with a relative dielectric constant (e r ) of 62 was applied.…”
Section: Simulation Parametersmentioning
confidence: 99%
“…74 Despite some loss in computational efficiency, using a single-range cutoff scheme was found to counterbalance those artefacts. 74 To correct for the truncation of electrostatic interactions beyond 1.4 nm a reactionfield correction 75 with a relative dielectric constant (e r ) of 62 was applied. The systems with 0.04 and 0.4 M [H 3 O + ] were simulated in triplicate and the reference system without H 3 O + ions was simulated in duplicate, each time using randomly assigned starting velocities.…”
Section: Simulation Parametersmentioning
confidence: 99%
“…charge group, reaction-field and twin-cutoff schemes). 55,56 The latest version of GROMACS (v.2016) 57 replaced the leapfrog algorithm by a variation of the Verlet-I/r-RESPA as the multiple-step alternative method to separately integrate fast-from slow-evolving interactions. The advantage of the new Verlet method 57 is its time-reversibility and symplecticity, but at the cost of great dependency of the update frequency on the nature of the system.…”
Section: Introductionmentioning
confidence: 99%
“…There are examples of changes assumed to be without effect on the results of MD simulations turning out to be important in the end (e.g. Reißer et al (2017)). The hypothesis that AMBER's dependence on unphysical details of the system description has no practical importance would have to be validated for all possible applications of the force field.…”
Section: Motivationmentioning
confidence: 99%