  In this paper we introduce a novel 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.
