朗之万方程

在统计物理中， 朗之万公式（保罗·朗之万，1908年) 是一个描述自由度的子集的时间演化的随机微分方程。 这些自由度，通常是那些在与系统的其他（微观的）变量相比，变化较缓慢的集体（宏观的）变量。 快速变化（微观）的变量导致了朗之万公式的随机性。

布朗运动为原型
原朗之万公式 描述了布朗运动，因受到流体分子的碰撞，粒子在流体中做无规则运动，
: m\frac{d^{2}\mathbf{x}}{dt^{2}}=-\lambda \frac{d\mathbf{x}}{dt}+\boldsymbol{\eta}\left( t\right).
这里，自由度是粒子的位置 \mathbf{x} ，m 表示粒子的质量。作用在粒子上的力表达成正比于粒子速度（斯托克斯定律）的粘滞力，和一个表示流体分子碰撞影响的噪声项\boldsymbol{\eta}\left( t\right) （随机微分方程中表示随机过程的术语在物理背景中的命名）的和。这个力（涨落力）\boldsymbol{\eta}\left( t\right) 具有高斯分布，其相关函数
: \left\langle \eta_{i}\left( t\right)\eta_{j}\left( t^{\prime}\right) \right\rangle =2\lambda k_{B}T\delta _{i,j}\delta \left(t-t^{\prime }\right) ,
其中 k_B 是玻尔兹曼常数,  T是温度， \eta_i\left( t\right) 是矢量 \boldsymbol{\eta}\left( t\right)的第 i 分量， δ-函数形式的时间相关性，表示假设该力在时刻 t， 与其他任何时刻完全不相关。这是一个近似，实际上随机力有一个与分子碰撞时长相对应的非零的相关时间。但是，朗之万方程是用来描述“宏观”微粒在很长时间尺度下的运动，并且在这种极限情况下 \delta-相关 和朗之万方程是精确的。

朗之万方程的另一个典型特征是在随机力的相关函数中导致了阻尼系数\lambda 出现，这一现象也被称为爱因斯坦关系。

数学方面
一个严格的 \delta-关系的涨落力 \boldsymbol{\eta}\left( t\right) 不是通常数学意义上的可微函数，即使它的一阶导数 d\mathbf{x}/dt 在这种极限下也没有定义。 要求朗之万方程在这种情况下的解释，可参见条目伊藤积分。

一般朗之万方程
经典力学有一个对一般朗之万方程的形式推导 。这个一般方程在临界动力学 和非平衡统计力学的其他领域扮演了核心角色。上述描述布朗运动的方程是一般朗之万方程的特殊情况。

一个推导一般朗之万方程的必要条件是对自由度不同快慢类型的标准划分（熵理论认为影响系统的变量可以分为快变量和慢变量）。例如，在液体中可以在几次碰撞时间内达到局部热力学平衡，但对于守恒量的密度，比如质量和能量，却需要长得多的时间去达到平衡。因此守恒量的密度，尤其是它们的长波分量，是慢变量的候选者。技术上来说这种划分是以Zwanzig投影算子 来实现的，它是推导中的必要工具。 推导不完全严格，因为它依赖于(貌似可信的)假设，类似于其他基本的统计力学中的假设。

令 A=\{A_i\} 表示慢变量。 则一般朗之万式表示为
: \frac{dA_{i}}{dt}=k_{B}T\sum\limits_{j}{\left[ {A_{i},A_{j}}\right] \frac{{d}\mathcal{H}}{{dA_{j}}}}-\sum\limits_{j}{\lambda _{i,j}\left( A\right) \frac{d\mathcal{H}}{{dA_{j}}}+}\sum\limits_{j}{\frac{d{\lambda _{i,j}\left(A\right) }}{{dA_{j}}}}+\eta _{i}\left( t\right).
涨落力 \eta_i\left( t\right) 服从高斯分布，其相关函数为
: \left\langle {\eta _{i}\left( t\right) \eta _{j}\left( t^{\prime }\right) }\right\rangle =2\lambda _{i,j}\left( A\right) \delta \left( t-t^{\prime}\right).
这暗示了阻尼系数 \lambda 具有昂萨格倒易关系 \lambda_{i,j}=\lambda_{j,i} 。 d\lambda_{i,j}/dA_j 对 A 的依赖性在大多数情况下可以忽略不计。符号 \mathcal{H}=-ln\left(p_0\right)表示了系统的哈密顿量，其中p_0\left(A\right) 是变量 A 的平衡概率分布。最后，[A_i,A_j]是慢变量 A_i 和 A_j 的泊松括号在慢变量空间投影。

在布朗运动的例子中，一个系统的状态可以有 \mathcal{H}=\mathbf{p}^2/\left(2mk_BT\right)，A=\{\mathbf{p}\} 或 A=\{\mathbf{x}, \mathbf{p}\} ， [x_i, p_j]=\delta_{i,j} 。 对 \mathbf{x} 的运动方程 d\mathbf{x}/dt=\mathbf{p}/m 是精确的， 其中没有涨落力 \eta_x 和阻尼力 \lambda_{x,p}.

例子
流体中的谐振子
一个不理想的谐振子会受到某些阻尼影响，由于涨落耗散定理，系统中一定会有一些波动。右图展示的是动量 p=mv 以及谐振子的位置 r 随时间演化的相图。 确定性的运动会沿着这条椭圆轨迹演化，并且不能与其他任何一条轨道交叉而不改变其能量。某些形式的阻尼的存在，例如分子流体环境（由扩散项和阻尼项为代表），会不断地从系统中得到或失去动能，导致一个谐振子的初始系综（图中虚线圈）会逐渐发散开，并最终成为正则系综（热平衡）。

电阻中的热噪声
上述的典型布朗颗粒，与约翰逊－奈奎斯特噪声，即由每个电阻中的热力学涨落引起的电压，有一个相似的类比。右图展示了包含一个电阻R和电容C的电路。这个电路的慢变量是电阻两端的的电压。其哈密顿量表示为\mathcal{H}=E/k_BT=CU^2/(2k_BT)，朗之万方程则表示为
: \frac{dU}{dt} =-\frac{U}{RC}+\eta \left( t\right),\;\;
\left\langle \eta \left( t\right) \eta \left( t^{\prime }\right)\right\rangle = \frac{2k_{B}T}{RC^{2}}\delta \left(t-t^{\prime }\right).
这个方程可以用来确定相关函数
: \left\langle U\left(t\right) U\left(t^{\prime }\right) \right\rangle
=\left( k_{B}T/C\right) \exp \left( -\left\vert t-t^{\prime }\right\vert
/RC\right) \approx 2Rk_{B}T\delta \left( t-t^{\prime }\right),
当电容C小可以忽略不计时成为白噪声（约翰逊噪声）。

临界动力学
二级相变的序参量的动力学在接近临界点时变慢，并且可以用朗之万方程描述推导得到
:
\frac{\partial P\left(A,t\right)}{\partial t}=\sum_{i,j}\frac{\partial}{\partial A_{i}}\left(-k_{B}T\left[A_{i},A_{j}\right]\frac{\partial\mathcal{H}}{\partial A_{j}}+\lambda_{i,j}\frac{\partial\mathcal{H}}{\partial A_{j}}+\lambda_{i,j}\frac{\partial}{\partial A_{j}}\right)P\left(A,t\right).
平衡分布 P(A) = p_0(A) = const\times\exp(-\mathcal{H}) 是一个平稳解。

路径积分
一个等价于朗之万方程的路径积分表述可以从相应的福克-普朗克方程得到，或通过将涨落力 \eta 的高斯分布 P^{(\eta)}(\eta)d\eta 转变成慢变量的概率分布，示意为 P(A)dA = P^{(\eta)}(\eta(A))\det(d\eta/dA)dA。如果以自然的（因果的）方式离散化朗之万方程，式中的函数行列式和相关的数学细节不证自明，A(t+\Delta t)-A(t)取决于 A(t) 而不是 A(t+\Delta t)。引入辅助反应变量 \tilde A 是方便的。等价于一般朗之万方程的路径积分表述为
: \int P(A,\tilde{A})\,dA\,d\tilde{A} = N\int \exp \left( L(A,\tilde{A})\right) dA\,d\tilde{A},
N 是归一化因子。路径积分表述没有引进任何新的东西，但它能够使用量子场论的工具，比如微扰论（摄动论）和重整化群方法（如果它们有意义的话）。

参见

• 朗之万动力学
• 泛函积分

阅读

• David Tong. Kinetic Theory Ch. 3.
• Applied Stochastic processes. M. Scott.

参考文献
延伸阅读

• W. T. Coffey (Trinity College, Dublin, Ireland) and Yu P. Kalmykov (Université de Perpignan, France, The Langevin Equation: With Applications to Stochastic Problems in Physics, Chemistry and Electrical Engineering (Third edition), World Scientific Series in Contemporary Chemical Physics - Vol 27.
• Reif, F. Fundamentals of Statistical and Thermal Physics, McGraw Hill New York, 1965. See section 15.5 Langevin Equation
• R. Friedrich, J. Peinke and Ch. Renner. How to Quantify Deterministic and Random Influences on the Statistics of the Foreign Exchange Market, Phys. Rev. Lett. 84, 5224 - 5227 (2000)
• L.C.G. Rogers and D. Williams. Diffusions, Markov Processes, and Martingales, Cambridge Mathematical Library, Cambridge University Press, Cambridge, reprint of 2nd (1994) edition, 2000.