Molecular dynamics simulations of phospholipid bilayers:: Influence of artificial periodicity, system size, and simulation time

This article investigates the convergence of structural and dynamical properties with system size and with time in molecular dynamics simulations of solvated phospholipid bilayers performed at constant volume under periodic boundary conditions using lattice-sum electrostatics. The electron density profile across the bilayer, the carbon-deuterium order parameters, and the surface tension are shown to be converged for a bilayer containing 36 lipids per leaflet and simulated over a period of 3-4 ns. Reasonable estimates for these properties can already be obtained from a system containing 16 lipids per leaflet. The convergence limit of 36 lipids per leaflet and the investigation of the correlation between lipid headgroup dipoles suggest a correlation length of about 3-5 nm in the lateral directions for a hydrated DPPC bilayer in the liquid-crystalline phase. Although these (relatively small) system sizes and (relatively short) time scales appear sufficient to obtain converged collective structural properties at constant volume, two restrictions should be kept in mind: (i) the relaxation times associated with the motion of individual lipids may be much longer and (ii) simulated properties converge significantly faster under constant volume conditions as compared to constant pressure conditions. Therefore, an accurate assessment of the dynamical properties of the system or of the relaxation of the bilayer under constant pressure conditions may require longer simulation time scales.