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.
 
From the BCS theory of superconductivity is well known that the superfluid density smoothly decreases with increasing temperature. Annihilated superfluid carriers become normal and lose their momenta on lattice atoms. So if we induce a persistent supercurrent in a ring below Tc and after that slowly increase the temperature, we must observe a decrease in the actual supercurrent, because the density of electron pairs and total supercurrent momentum decrease. However, this supercurrent...
Back
Top