Hybrid simulations of lateral diffusion in fluctuating membranes

In this paper we introduce a method to simulate lateral diffusion of inclusions in a fluctuating membrane. The regarded systems are governed by two dynamic processes: the height fluctuations of the membrane and the diffusion of the inclusion along the membrane. While membrane fluctuations can be expressed in terms of a dynamic equation which follows from the Helfrich Hamiltonian, the dynamics of the diffusing particle is described by a Langevin or Smoluchowski equation. In the latter equations, the curvature of the surface needs to be accounted for, which makes particle diffusion a function of membrane fluctuations. In our scheme these coupled dynamic equations, the membrane equation and the Langevin equation for the particle, are numerically integrated to simulate diffusion in a membrane. The simulations are used to study the ratio of the diffusion coefficient projected on a flat plane and the intramembrane diffusion coefficient for the case of free diffusion. We compare our results with recent analytical results that employ a preaveraging approximation and analyze the validity of this approximation. A detailed simulation study of the relevant correlation functions reveals a surprisingly large range where the approximation is applicable.