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
    you can find a particular algorithm in a 1998 j chem phys paper by chandler and dellago.

Know someone interested in this topic? Share this thread via Reddit, Google+, Twitter, or Facebook

Have something to add?

Similar Discussions: Langevin dynamics random force term generation algorithm
  1. Langevin curve fitting (Replies: 3)

  2. Metropolis Algorithm (Replies: 4)