Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Langevin dynamics random force term generation algorithm

  1. Aug 11, 2008 #1

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

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


    What I just wrote is from the Wikipedia article "Langevin dynamics"
    and R(t) belongs to the simple Langevin equation
    [tex]F(x) = m \ddot{x} = -\nabla U(x) - \gamma m \dot{x} + R(t)[/tex]
    Last edited: Aug 11, 2008
  2. jcsd
  3. Feb 23, 2009 #2
    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 havent found one.
  4. Dec 4, 2009 #3
Share this great discussion with others via Reddit, Google+, Twitter, or Facebook