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

布朗运动为原型
-------

原朗之万公式 描述了 布朗运动，因受到流体分子的碰撞，粒子在流体中做无规则运动，

m\frac{d^{2}\mathbf{x}}{dt^{2}}=-\lambda \frac{d\mathbf{x}}{dt}+\boldsymbol{\eta}\left( t\right).
这里，自由度是粒子的位置 \mathbf {x}  ，m 表示粒子的质量。作用在粒子上的力表达成正比于粒子速度（斯托克斯定律）的粘滞力，和一个表示流体分子碰撞影响的噪声项{\displaystyle {\boldsymbol {\eta }}\left(t\right)} （随机微分方程中表示随机过程的术语在物理背景中的命名）的和。这个力（涨落力）{\displaystyle {\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是温度， {\displaystyle \eta \_{i}\left(t\right)} 是矢量 {\displaystyle {\boldsymbol {\eta }}\left(t\right)}的第 i 分量， δ-函数形式的时间相关性，表示假设该力在时刻 t， 与其他任何时刻完全不相关。这是一个近似，实际上随机力有一个与分子碰撞时长相对应的非零的相关时间。但是，朗之万方程是用来描述"宏观"微粒在很长时间尺度下的运动，并且在这种极限情况下 \delta -相关 和朗之万方程是精确的。

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

数学方面
----

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

一般朗之万方程
-------

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

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

令 {\displaystyle A=\{A\_{i}\}} 表示慢变量。 则一般朗之万式表示为

{\displaystyle {\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).}
涨落力 {\displaystyle \eta \_{i}\left(t\right)} 服从高斯分布，其相关函数为

{\displaystyle \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  具有昂萨格倒易关系 {\displaystyle \lambda \_{i,j}=\lambda \_{j,i}} 。 {\displaystyle d\lambda \_{i,j}/dA\_{j}} 对 A 的依赖性在大多数情况下可以忽略不计。符号 {\displaystyle {\mathcal {H}}=-ln\left(p\_{0}\right)}表示了系统的哈密顿量，其中{\displaystyle p\_{0}\left(A\right)} 是变量 A 的平衡概率分布。最后，{\displaystyle [A\_{i},A\_{j}]}是慢变量 A\_{i} 和 A\_{j} 的泊松括号在慢变量空间投影。

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

例子
--

一个谐振子相图展示了由郎之万方程决定的相空间传播
### 流体中的谐振子

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

包含一个电阻和一个电容的电路
### 电阻中的热噪声

上述的典型布朗颗粒，与约翰逊－奈奎斯特噪声，即由每个电阻中的热力学涨落引起的电压，有一个相似的类比。右图展示了包含一个电阻R和电容C的电路。这个电路的慢变量是电阻两端的的电压。其哈密顿量表示为{\displaystyle {\mathcal {H}}=E/k\_{B}T=CU^{2}/(2k\_{B}T)}，朗之万方程则表示为

{\displaystyle {\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).}
这个方程可以用来确定相关函数

{\displaystyle \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小可以忽略不计时成为白噪声（约翰逊噪声）。

### 临界动力学

二级相变的序参量的动力学在接近临界点时变慢，并且可以用朗之万方程描述。最简单的例子是具有非保守标量阶参量的普适类 "model A"，在轴向铁磁体中的实现

{\displaystyle {\begin{aligned}{\frac {\partial \varphi \left(\mathbf {x} ,t\right)}{\partial t}}&=-\lambda {\frac {\delta {\mathcal {H}}}{\delta \varphi }}+\eta \left(\mathbf {x} ,t\right),\\{\mathcal {H}}&=\int d^{d}x\left\{{\frac {1}{2}}\varphi \left[r\_{0}-\nabla ^{2}\right]\varphi +u\varphi ^{4}\right\},\\\left\langle \eta \left(\mathbf {x} ,t\right)\eta \left(\mathbf {x} ',t'\right)\right\rangle &=2\lambda \delta \left(\mathbf {x} -\mathbf {x} '\right)\delta \left(t-t'\right).\end{aligned}}}
其他的普适类（命名方法是像 "model A",... "model J"）包含一个扩散的序参量，有几个分量的序参量，其他临界变量和（或）来自泊松括号的贡献 。

### 重现玻尔兹曼分布

朗之万方程必须能重现玻尔兹曼分布。一维过阻尼布朗运动是一个有启发性的例子。在颗粒的惯性相对于阻尼力来说可以忽略不计时，就实现了过阻尼条件。在势能 V(x)下粒子的轨迹 x(t) 用朗之万方程描述

{\displaystyle \lambda {\frac {dx}{dt}}=-{\frac {\partial V(x)}{\partial x}}+\eta (t),}
噪音项的特征由 {\displaystyle \left\langle \eta (t)\eta (t')\right\rangle =2k\_{B}T\lambda \delta (t-t')} 决定，其中 \lambda  是阻尼系数。我们想要计算经过一段时间粒子位置的分布 p(x)。确定这个分布一个直接的方法是引进一个测试函数f，然后看这个函数在所有实现中的平均（统计均值）

{\displaystyle \lambda {\frac {d\left\langle f(x(t))\right\rangle }{dt}}=\left\langle f'(x(t))\lambda {\frac {dx}{dt}}\right\rangle =\left\langle -f'(x(t)){\frac {\partial V}{\partial x}}+f'(x(t))\eta (t)\right\rangle .}
如果 x(t) 保持有限值那么这个量是没有意义的。此外，使用斯特拉托诺维奇诠释，我们就可以摆脱第二项中的\eta 从而我们最终得到

{\displaystyle \left\langle -f'(x){\frac {\partial V}{\partial x}}+k\_{B}Tf''(x)\right\rangle =0,}
这里我们利用了概率密度函数 p(x)。通过显式计算平均值

{\displaystyle \int \left(-f'(x){\frac {\partial V}{\partial x}}p(x)+{k\_{B}T}f''(x)p(x)\right)dx=\int \left(-f'(x){\frac {\partial V}{\partial x}}p(x)-{k\_{B}T}f'(x)p'(x)\right)dx=0,}
第二项为分部积分（因此有负号）。因为需要对任意的函数 f 都成立，所以我们必须有

{\displaystyle {\frac {\partial V}{\partial x}}p(x)+{k\_{B}T}p'(x)=0,}
因此恢复为玻尔兹曼函数

{\displaystyle p(x)\propto \exp \left({-{\frac {V(x)}{k\_{B}T}}}\right).}

等价的技巧
-----

在涨落力的具体实现中，朗之万方程并不关心其本身的解，它关心的是在对涨落力取平均后慢变量的相关函数。这样的相关函数也可以用其他（等价的）技巧确定。

### 福克-普朗克方程

福克-普朗克方程是关于随机变量A的含时概率密度{\displaystyle P\left(A,t\right)}的一个确定性方程。对应上面一般朗之万方程的福克-普朗克方程可以由标准技巧推导得到

{\displaystyle {\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).}
平衡分布 {\displaystyle P(A)=p\_{0}(A)=const\times \exp(-{\mathcal {H}})} 是一个平稳解。

### 路径积分

一个等价于朗之万方程的路径积分表述可以从相应的福克-普朗克方程得到，或通过将涨落力 \eta  的高斯分布 {\displaystyle P^{(\eta )}(\eta )d\eta } 转变成慢变量的概率分布，示意为 {\displaystyle P(A)dA=P^{(\eta )}(\eta (A))\det(d\eta /dA)dA}。如果以自然的（因果的）方式离散化朗之万方程，式中的函数行列式和相关的数学细节不证自明，{\displaystyle A(t+\Delta t)-A(t)}取决于 A(t) 而不是 {\displaystyle A(t+\Delta t)}。引入辅助反应变量 {\displaystyle {\tilde {A}}} 是方便的。等价于一般朗之万方程的路径积分表述为

{\displaystyle \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.

参考文献
----

1. ↑  Langevin, P.. C. R. Acad. Sci. (Paris). 1908, **146**: 530–533.
2. ↑  Kawasaki, K.. J. Phys. A: Math. Nucl. Gen. 1973, **6**: 1289. Bibcode:1973JPhA....6.1289K. doi:10.1088/0305-4470/6/9/004.
3. ↑  Dengler, R.. 2015. arXiv:1506.02650v2 可免费查阅.
4. 1 2 3  Hohenberg, P. C.; Halperin, B. I.. Reviews of Modern Physics. 1977, **49** (3): 435–479. Bibcode:1977RvMP...49..435H. doi:10.1103/RevModPhys.49.435.
5. ↑  Zwanzig, R.. Phys. Rev. 1961, **124** (4): 983–992. Bibcode:1961PhRv..124..983Z. doi:10.1103/PhysRev.124.983.
6. ↑  J. Johnson, "Thermal Agitation of Electricity in Conductors", Phys.
7. ↑  Ichimaru, S.  1st. USA: Benjamin: 231, 1973, ISBN 0805387536
8. ↑  Janssen, H. K.. Z. Phys. B. 1976, **23**: 377. Bibcode:1976ZPhyB..23..377J. doi:10.1007/BF01316547.
**延伸阅读**

* 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.

|  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |
| --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- |
| 

| * 
* 论
* 编
概率论：随机过程 |
| --- |
|  |
| 离散时间（英语：） | 
* 伯努利过程
* 分支过程
* 中餐馆过程
* 高尔顿-沃特森过程（英语：）
* 独立同分布
* 马尔可夫链
* 莫兰过程（英语：）
* 隨機漫步
	+ 循环擦除随机游走（英语：）
	+ 自避行走

 |
|  |
| 连续时间（英语：） | 
* 贝塞尔过程
* 出生-死亡過程
* 维纳过程/布朗运动
	+ 布朗桥
	+ Excursion（英语：）
	+ 分数布朗运动（英语：）
	+ 几何布朗运动
	+ Meander（英语：）
* 柯西过程（英语：）
* Contact process（英语：）
* Cox process（英语：）
* Diffusion process（英语：）
* Empirical process（英语：）
* 费勒过程（英语：）
* 弗莱明-维奥过程（英语：）
* 伽马过程（英语：）
* 亨特过程（英语：）
* Interacting particle system（英语：）s
* 伊藤积分
* 伊藤過程
* 跳跃扩散（英语：）
* 跳跃过程
* 萊維過程
* Local time（英语：）
* 马尔可夫加过程（英语：）
* 麦基恩-弗拉索夫过程（英语：）
* 奥恩斯坦-乌伦贝克过程
* 泊松过程
	+ 复合泊松过程（英语：）
	+ 非齐次泊松过程
	+ 泊松点过程
* 施拉姆-勒夫纳演进
* 半鞅
* Sigma-martingale（英语：）
* Stable process（英语：）
* Superprocess（英语：）
* Telegraph process（英语：）
* Variance gamma process（英语：）
* 维纳过程
* Wiener sausage（英语：）

 |
|  |
| 离散时间与连续时间 | 
* 分支過程
* 高斯过程
* 隐马尔可夫模型（HMM）
* 馬可夫過程
* 鞅
	+ 鞅差序列（英语：）
	+ 局部鞅（英语：）
	+ Sub-
	+ Super-（英语：）
* Random dynamical system（英语：）
* Regenerative process（英语：）
* Renewal process（英语：）
* 白雜訊

 |
|  |
| 场及其它 | 
* 狄利克雷过程（英语：）
* 高斯隨機場（英语：）
* 吉布斯测度（英语：）
* 霍普菲尔德神经网络
* 易辛模型
* 马尔可夫网络
* 渗流理论
* 皮特曼-约尔过程（英语：）
* 点过程（英语：）
	+ Cox（英语：）
	+ 泊松过程
* 玻茨模型
* 随机场
* 随机图

 |
|  |
| 时间序列模型 | 
* ARCH模型
* ARIMA模型
* 自我迴歸模型
* ARMA模型
* 广义ARCH模型
* 移动平均模型

 |
|  |
| 金融模型 | 
* 布莱克-德尔曼-托伊模型（英语：）
* 布莱克-卡拉辛斯基模型（英语：）
* 布莱克-舒尔斯模型
* 陈模型
* Constant elasticity of variance (CEV)（英语：）
* 科克斯-英格索尔-罗斯模型 (CIR)（英语：）
* Garman–Kohlhagen（英语：）
* HJM框架
* 赫斯顿模型（英语：）
* Ho–Lee（英语：）
* 赫爾-懷特模型
* LIBOR市场模型（英语：）
* SABR volatility（英语：）
* 瓦西塞克模型（英语：）

 |
|  |
| 精算學 | 
* Bühlmann（英语：）
* Cramér–Lundberg（英语：）
* Risk process（英语：）
* Sparre–Anderson（英语：）

 |
|  |
| 等候理論 | 
* Bulk（英语：）
* Fluid（英语：）
* Generalized queueing network（英语：）
* M/G/1（英语：）
* M/M/1
* M/M/c（英语：）

 |
|  |
| 性质 | 
* 右连左极函数
* Continuous（英语：）
* Continuous paths（英语：）
* 遍历性
* Exchangeable（英语：）
* Feller-continuous（英语：）
* Gauss–Markov（英语：）
* 马尔可夫性质
* Mixing（英语：）
* Piecewise deterministic（英语：）
* 可预测过程
* 循序可测过程
* Self-similar（英语：）
* 平稳过程
* Time-reversible（英语：）

 |
|  |
| 极限定理 | 
* 中心极限定理
* Donsker's theorem（英语：）
* Doob's martingale convergence theorems（英语：）
* 遍历理论
* Fisher–Tippett–Gnedenko theorem（英语：）
* Large deviation principle（英语：）
* 大數法則
* 重对数律
* Maximal ergodic theorem（英语：）
* Sanov's theorem（英语：）

 |
|  |
| 不等式 | 
* Burkholder–Davis–Gundy（英语：）
* Doob's martingale（英语：）
* Kunita–Watanabe（英语：）

 |
|  |
| 工具 | 
* Cameron–Martin formula（英语：）
* 随机变量的收敛
* Doléans-Dade exponential（英语：）
* Doob decomposition theorem（英语：）
* Doob–Meyer decomposition theorem（英语：）
* Doob's optional stopping theorem（英语：）
* Dynkin's formula（英语：）
* 费曼-卡茨公式
* 右连左极函数
* Girsanov theorem（英语：）
* Infinitesimal generator（英语：）
* 伊藤积分
* 伊藤引理
* Kolmogorov continuity theorem（英语：）
* Kolmogorov extension theorem（英语：）
* Lévy–Prokhorov metric（英语：）
* Malliavin calculus（英语：）
* Martingale representation theorem（英语：）
* Optional stopping theorem（英语：）
* Prohorov theorem（英语：）
* 二次變差
* Reflection principle（英语：）
* Skorokhod integral（英语：）
* Skorokhod's representation theorem（英语：）
* 右连左极函数
* Snell envelope（英语：）
* 隨機微分方程
	+ Tanaka（英语：）
* 停时
* 隨機积分
* Uniform integrability（英语：）
* Usual hypotheses（英语：）
* 维纳空间
	+ Classical（英语：）
	+ Abstract
* 漂移项

 |
|  |
| 相关领域 | 
* 精算學
* 计量经济学
* 遍历理论
* 极值理论（EVT）
* Large deviations theory（英语：）
* 數理金融學
* 数理统计学
* 概率论
* 等候理論
* Renewal theory（英语：）
* Ruin theory（英语：）
* 统计学
* 随机分析
* 时间序列分析
* 机器学习

 |
|  |
| 
* 分类**分类**

 |

 |