Accurate and efficient integration for molecular dynamics simulations at constant temperature and pressure

Abstract

In molecular dynamics simulations, control over temperature and pressure is typically achieved by augmenting the original system with additional dynamical variables to create a thermostat and a barostat, respectively. These variables generally evolve on timescales much longer than those of particle motion, but typical integrator implementations update the additional variables along with the particle positions and momenta at each time step. We present a framework that replaces the traditional integration procedure with separate barostat, thermostat, and Newtonian particle motion updates, allowing thermostat and barostat updates to be applied infrequently. Such infrequent updates provide a particularly substantial performance advantage for simulations parallelized across many computer processors, because thermostat and barostat updates typically require communication among all processors. Infrequent updates can also improve accuracy by alleviating certain sources of error associated with limited-precision arithmetic. In addition, separating the barostat, thermostat, and particle motion update steps reduces certain truncation errors, bringing the time-average pressure closer to its target value. Finally, this framework, which we have implemented on both general-purpose and special-purpose hardware, reduces software complexity and improves software modularity.

EGFR oligomerization organizes kinase-active dimers into competent signalling platforms

“…The force constant was linearly scaled down to zero over at least 50 ns. Simulations were performed in the constant number (N), pressure (P), and temperature (T) (NPT) ensemble with T =310 K and P =1 bar using a variant70 of the Nosé–Hoover71 and the Martyna–Tobias–Klein algorithm72. Water molecules and all bond lengths to hydrogen atoms were constrained using M-SHAKE73.…”

Transmembrane Pore Structures of β-Hairpin Antimicrobial Peptides by All-Atom Simulations

“…The production simulations were carried out on the Anton 2 supercomputer 59 at the Pittsburgh Supercomputing Center using the Multigrator integration framework 60 . We used VMD to generate Anton 2 input files from the equilibrated NAMD output and parameter files, and then we used the viparr program to add the force field information to the Anton input files.…”

Atomistic Modeling of Ion Conduction through the Voltage-Sensing Domain of the Shaker K+ Ion Channel

“…The simulations were performed at constant temperature (300 K) and pressure (1 atm), using Nose-Hoover chains 41 and the Martyna-Tobias-Klein barostat 36 . The RESPA algorithm and the temperature and pressure controls were implemented using the multigrator scheme 42 .…”