Langevin dynamics random force term generation algorithm

[iibewegung]

Hi,

Can anyone tell me an algorithm to generate the stationary Gaussian distribution R(t) with

$$\langle R(t) \rangle = 0$$ (zero mean)
$$\langle R(t) R(t')^{T} \rangle = A \delta(t-t')$$, $$A = 2 \gamma k_B T m$$ (autocorrelation)

?

What I just wrote is from the Wikipedia article "Langevin dynamics"
and R(t) belongs to the simple Langevin equation
$$F(x) = m \ddot{x} = -\nabla U(x) - \gamma m \dot{x} + R(t)$$

 

[jwallace]

I have been looking for the same code. It is not exactly trivial. I found some code in the following book
"The Molecular Dynamics of Liquid Crystals,by G. R. Luckhurst, C. A. Veracini"
I have been looking for a DJVU copy of this book but haven't found one.

 

[nkpreketes]

you can find a particular algorithm in a 1998 j chem phys paper by chandler and dellago.

http://gold.cchem.berkeley.edu/Pubs/DC150.pdf