Langevin方程数值积分与蒙特卡洛模拟技术.pdfVIP

Langevin方程的数值积分

方程：蒙特卡洛模拟

・L・厄和海伦・巴克

环境科学部，加州大学利弗莫尔，加利福尼亚州利弗莫尔94550

收稿日期：；修订日期：

本文推导了在外部力作用下布朗粒子普通Langevin运动方程的蒙特卡洛模拟技术。

这些方法允许在选择时间步长时有较大的自由度，时间步长仅受外部力变化率的

限制。这种方法被扩展到使用摩擦力项中的记忆函数的广义Langevin方程。推导

出了一般性的模拟技术，这些技术独立于记忆函数的形式。对于指数记忆函数的

情况，提出了一种需要较少空间的特殊方法。

1.引言

布朗运动的朗之万模型已被广泛用于描述在热浴中处于热力学平衡条件下的粒子

动态行为。该模型以统计方式描述布朗粒子与周围热浴之间的相互作用，而不是

尝试处理每个单独的布朗粒子‑热浴粒子相互作用。假设布朗粒子由于热浴而受到

的力由两部分组成：一部分是与布朗粒子速度成正比的摩擦力，另一部分是只能

通过其统计特性来描述的随机波动力。通过使用这种力在牛顿运动方程中，可以

得到朗之万方程。

求解朗之万方程与求解普通运动方程不同，后者在相空间中的解是一条唯一的轨

迹。朗之万方程的随机性质允许对于任何给定的