The particles in the lipid bilayer undergo fluctuations from their equilibrium positions. These motions are best described in terms of correlated monolayer motions (undulations), anticorrelated monolayer motions (thickness fluctuations) and number density correlations, both within and between the monolayers. Particle-based computer simulations, atomistic or coarse-grained, can be used to calculate the spectra of these equilibrium fluctuations. This has conventionally been done by interpolating the particle structure onto a grid to facilitate for the discrete Fourier transform. However, the gridding introduces artifacts and an upper cutoff on the available wave vectors. Here, we instead analyze the simulations directly in Fourier space, circumventing the need for a grid approximation. The fluctuation spectra can then be calculated at high precision up to any wave vector desired, and the low wave vectors are only bounded by the size of the lateral simulation box. From the analysis, a picture emerges which bridges a continuum regime at low wave vectors, and a domain that can be attributed to the in-plane bilayer molecular structure at high wave vectors. We calculate spectra for undulations, thickness fluctuations, and the spectrum of the number density along the undulating surface, and we discuss how the shape of the spectra in the low-wave vector and high-wave vector regimes are related to membrane material constants, and give some predictions from an analytical theory based on equilibrium correlation functions.
Read full abstract