J. Phys. Chem. B, 108 (1), 363 -368, 2004. 10.1021/jp030458y S1089-5647(03)00458-9
Web Release Date: December 4, 2003

Not subject to U.S. Copyright. Published 2004 American Chemical Society

Pressure-Based Long-Range Correction for Lennard-Jones Interactions in Molecular Dynamics Simulations: Application to Alkanes and Interfaces

Patrick Lagüe and Richard W. Pastor

Laboratory of Biophysics, Center for Biologics Evaluation and Research, FDA, 1401 Rockville Pike, Rockville, Maryland 20852-1448

Bernard R. Brooks*

Laboratory of Structural Biology, Division of Computer Research & Technology, National Institutes of Health, Bethesda, Maryland 20892

Received: April 17, 2003

In Final Form: October 3, 2003

Abstract:

A straightforward method that accounts for the long-range Lennard-Jones (LJ) terms in constant pressure molecular dynamics simulations is presented. This long-range correction (LRC) consists of an additional applied pressure tensor which is periodically calculated from the difference of instantaneous pressures at the selected cutoff and a very long cutoff. It provides results that are nearly independent of the LJ cutoff distance at negligible additional calculation costs, and is particularly suited for anisotropic systems such as liquid/liquid interfaces or heterogeneous macromolecules where approximations based on spherically symmetric radial distribution functions are expected to fail. The utility of the method is demonstrated for a series of alkanes and water, and for interfaces including a lipid bilayer. The LRC increases densities and decreases isothermal compressibilities, with the changes larger for alkanes than for water (where the long-range interactions are dominated by electrostatic interactions). While implementation of the LRC will not necessarily improve agreement with a particular experiment, it will provide a baseline for improvements in a parameter set that are consistent with the long-range Lennard-Jones interactions.


Download the full text: PDF | HTML