Hi,(adsbygoogle = window.adsbygoogle || []).push({});

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]

**Physics Forums - The Fusion of Science and Community**

Dismiss Notice

Join Physics Forums Today!

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

# Langevin dynamics random force term generation algorithm

Loading...

Similar Threads for Langevin dynamics random | Date |
---|---|

I Boundary Regions in LAMMPS MD Simulation | Mar 14, 2017 |

I Scaling and the renormalization (sub)group | Dec 31, 2016 |

A Can a molecular dynamics simulation enter a limit cycle? | Oct 20, 2016 |

A Free energy density in liquid crystals | May 14, 2016 |

Langevin curve fitting | Jun 24, 2008 |

**Physics Forums - The Fusion of Science and Community**