Langevin dynamics random force term generation algorithm

iibewegung
Messages
14
Reaction score
0
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)
 
Last edited:
Physics news on Phys.org
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.
 
Back
Top