Translational diffusion coefficients are primary reporters on the dynamics in lipid membranes. However, diffusion coefficients of lipids and membrane proteins calculated in molecular dynamics (MD) simulations suffer from severe finite-size effects. We show that these finite-size effects can be quantified in terms of hydrodynamic theory and develop accurate correction formulas. Using large-scale coarse-grained MD simulations of up to 132 million particles, we assess the effects of box size and shape on translational diffusion in lipid membrane simulations. For the flat periodic boxes commonly used for membrane simulations, we find that the diffusion coefficients of lipids and membrane proteins increase logarithmically (and therefore without bound) with the box width. Only if the box size is increased also normal to the membrane surface, the diffusion coefficient converges. This dependence on both the size and the shape of the simulation box is captured quantitatively by a hydrodynamic theory. The corresponding correction formulas make it possible to perform meaningful comparisons of diffusion coefficients from simulations with different box sizes, and facilitate the comparison to experiment.
Read full abstract