Direct-Space Corrections Enable Fast and Accurate Lorentz-Berthelot Combination Rule Lennard-Jones Lattice Summation.

Long-range lattice summation techniques such as the particle-mesh Ewald (PME) algorithm for electrostatics have been revolutionary to the precision and accuracy of molecular simulations in general. Despite the performance penalty associated with lattice summation electrostatics, few biomolecular simulations today are performed without it. There are increasingly strong arguments for moving in the same direction for Lennard-Jones (LJ) interactions, and by using geometric approximations of the combination rules in reciprocal space, we have been able to make a very high-performance implementation available in GROMACS. Here, we present a new way to correct for these approximations to achieve exact treatment of Lorentz-Berthelot combination rules within the cutoff, and only a very small approximation error remains outside the cutoff (a part that would be completely ignored without LJ-PME). This not only improves accuracy by almost an order of magnitude but also achieves absolute biomolecular simulation performance that is an order of magnitude faster than any other available lattice summation technique for LJ interactions. The implementation includes both CPU and GPU acceleration, and its combination with improved scaling LJ-PME simulations now provides performance close to the truncated potential methods in GROMACS but with much higher accuracy.

Performance enhancements for GROMACS nonbonded interactions on BlueGene.

Several improvements to the previously optimized GROMACS BlueGene inner loops that evaluate nonbonded interactions in molecular dynamics simulations are presented. The new improvements yielded an 11% decrease in running time for both PME and other kinds of GROMACS simulations that use nonbonded table look-ups. Some other GROMACS simulations will show a small gain.

Optimization of parameters for molecular dynamics simulation using smooth particle-mesh Ewald in GROMACS 4.5.

Based on our critique of requirements for performing an efficient molecular dynamics simulation with the particle-mesh Ewald (PME) implementation in GROMACS 4.5, we present a computational tool to enable the discovery of parameters that produce a given accuracy in the PME approximation of the full electrostatics. Calculations on two parallel computers with different processor and communication structures showed that a given accuracy can be attained over a range of parameter space, and that the attributes of the hardware and simulation system control which parameter sets are optimal. This information can be used to find the fastest available PME parameter sets that achieve a given accuracy. We hope that this tool will stimulate future work to assess the impact of the quality of the PME approximation on simulation outcomes, particularly with regard to the trade-off between cost and scientific reliability in biomolecular applications.

Ensuring Mixing Efficiency of Replica-Exchange Molecular Dynamics Simulations.

We address the question of constructing a protocol for replica-exchange molecular dynamics (REMD) simulations that make efficient use of the replica space, assess whether published applications are achieving such "mixing" efficiency, and provide a how-to guide to assist users to plan efficient REMD simulations. To address our first question, we introduce and discuss three metrics for assessing the number of replica-exchange attempts required to justify the use of a replica scheme and define a "transit number" as the lower bound for the length of an efficient simulation. Our literature survey of applications of REMD simulations of peptides in explicit solvent indicated that authors are not routinely reporting sufficient details of their simulation protocols to allow readers to make independent assessments of the impact of the method on their results, particularly whether mixing efficiency has been achieved. Necessary details include the expected or observed replica-exchange probability, together with the total number of exchange attempts, the exchange period, and estimates of the autocorrelation time of the potential energy. Our analysis of cases where the necessary information was reported suggests that in many of these simulations there are insufficient exchanges attempted or an insufficiently long period between them to provide confidence that the simulation length justifies the size of the replica scheme. We suggest guidelines for designing REMD simulation protocols to ensure mixing efficiency. Two key recommendations are that the exchange period should in general be larger than 1 ps and the number of exchange attempts should be chosen to significantly exceed the transit number for the replica scheme.