Undulation contributions to the area compressibility in lipid bilayer simulations.

It is here shown that there is a considerable system size-dependence in the area compressibility calculated from area fluctuations in lipid bilayers. This is caused by the contributions to the area fluctuations from undulations. This is also the case in experiments. At present, such a contribution, in most cases, is subtracted from the experimental values to obtain a true area compressibility. This should also be done with the simulation values. Here, this is done by extrapolating area compressibility versus system size, down to very small (zero) system size, where undulations no longer exist. The area compressibility moduli obtained from such simulations do not agree with experimental true area compressibility moduli as well as the uncorrected ones from contemporary or earlier simulations, but tend, instead, to be approximately 50% too large. As a byproduct, the bending modulus can be calculated from the slope of the compressibility modulus versus system-size. The values obtained in this way for the bending modulus are then in good agreement with experiment.