朗之万方程

✍ dations ◷ 2025-04-03 12:45:37 #随机微分方程,方程,统计力学

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

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

这里,自由度是粒子的位置  x {\displaystyle \mathbf {x} }  , m {\displaystyle m} 表示粒子的质量。作用在粒子上的力表达成正比于粒子速度(斯托克斯定律)的粘滞力,和一个表示流体分子碰撞影响的噪声项 η ( t ) {\displaystyle {\boldsymbol {\eta }}\left(t\right)} (随机微分方程中表示随机过程的术语在物理背景中的命名)的和。这个力(涨落力) η ( t ) {\displaystyle {\boldsymbol {\eta }}\left(t\right)} 具有高斯分布,其相关函数

其中 k B {\displaystyle k_{B}} 是玻尔兹曼常数,   T {\displaystyle T} 是温度, η i ( t ) {\displaystyle \eta _{i}\left(t\right)} 是矢量 η ( t ) {\displaystyle {\boldsymbol {\eta }}\left(t\right)} 的第 i 分量, δ-函数形式的时间相关性,表示假设该力在时刻 t, 与其他任何时刻完全不相关。这是一个近似,实际上随机力有一个与分子碰撞时长相对应的非零的相关时间。但是,朗之万方程是用来描述“宏观”微粒在很长时间尺度下的运动,并且在这种极限情况下  δ {\displaystyle \delta } -相关 和朗之万方程是精确的。

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

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

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

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

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

涨落力 η i ( t ) {\displaystyle \eta _{i}\left(t\right)} 服从高斯分布,其相关函数为

这暗示了阻尼系数 λ {\displaystyle \lambda } 具有昂萨格倒易关系 λ i , j = λ j , i {\displaystyle \lambda _{i,j}=\lambda _{j,i}} d λ i , j / d A j {\displaystyle d\lambda _{i,j}/dA_{j}} A {\displaystyle A} 的依赖性在大多数情况下可以忽略不计。符号 H = l n ( p 0 ) {\displaystyle {\mathcal {H}}=-ln\left(p_{0}\right)} 表示了系统的哈密顿量,其中 p 0 ( A ) {\displaystyle p_{0}\left(A\right)} 是变量 A {\displaystyle A} 的平衡概率分布。最后, {\displaystyle } 是慢变量 A i {\displaystyle A_{i}} A j {\displaystyle A_{j}} 的泊松括号在慢变量空间投影。

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

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

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

这个方程可以用来确定相关函数

当电容C小可以忽略不计时成为白噪声(约翰逊噪声)。

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

其他的普适类(命名方法是像 “model A",..., "model J")包含一个扩散的序参量,有几个分量的序参量,其他临界变量和(或)来自泊松括号的贡献 。

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

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

如果 x ( t ) {\displaystyle x(t)} 保持有限值那么这个量是没有意义的。此外,使用斯特拉托诺维奇诠释,我们就可以摆脱第二项中的 η {\displaystyle \eta } 从而我们最终得到

这里我们利用了概率密度函数 p ( x ) {\displaystyle p(x)} 。通过显式计算平均值

第二项为分部积分(因此有负号)。因为需要对任意的函数 f {\displaystyle f} 都成立,所以我们必须有

因此恢复为玻尔兹曼函数

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

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

平衡分布 P ( A ) = p 0 ( A ) = c o n s t × exp ( H ) {\displaystyle P(A)=p_{0}(A)=const\times \exp(-{\mathcal {H}})} 是一个平稳解。

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

N {\displaystyle N} 是归一化因子。路径积分表述没有引进任何新的东西,但它能够使用量子场论的工具,比如微扰论(摄动论)和重整化群方法(如果它们有意义的话)。

延伸阅读

相关

  • 淡水河淡水河,流经中华人民共和国广东省深圳市和惠州市,是西枝江左岸支流,上游也称龙岗河,发源于深圳市梧桐山北麓,支流有坪山河、横岭水向北流经深圳市龙岗区和惠州市惠阳区,最后于惠州
  • 七松之役七松之役(Battle of Seven Pines)于1862年5月31日至6月1日爆发于美国维吉尼亚州,为南北战争的半岛会战(Peninsula Campaign)中的一场战事。南军发现北军有两个兵团孤立于奇克哈默
  • 慕尼黑大学诺贝尔奖得主列表慕尼黑大学(英语:;德语:)是德国著名大学之一。大学成立时间可追溯到1471年。慕尼黑大学历史上名人辈出,在自然科学领域有杰出表现。历史上共有34位校友荣获诺贝尔奖。
  • 埃塞俄比亚正统台瓦西多教会埃塞俄比亚正统台瓦西多教会(阿姆哈拉语:የኢትዮጵያ ኦርቶዶክስ ተዋሕዶ ቤተ ክርስቲያን,罗马化:Yäityop'ya ortodoks täwahedo bétäkrestyan)是一个基督教东方正统
  • 信息空间网络空间(英文:, Cyberspace),又译信息空间、赛博空间、电子世界、网络世界等,是哲学和电脑领域中的一个抽象概念,指在电脑以及数字网络里的虚拟现实。“网络空间”一词是控制论(cyb
  • 京都府第3区京都府第3区是日本众议院的选区,始于1994年。北海道 13 | 山形县 4 | 静冈县 9 | 岛根县 3 | 大分县 4福井县 3 | 山梨县 3 | 德岛县 3 | 高知县 3 | 佐贺县 3青森县 4 | 岩
  • 陈田陈田(1855年-?),字松珊,贵州省贵阳府贵筑县(今贵州省贵阳市)人,清朝政治人物、进士出身。光绪十二年,登进士;同年五月,改翰林院庶吉士。光绪十五年四月,散馆,授翰林院编修。次年,任国史馆协
  • 精神益生菌精神益生菌(英语:Psychobiotics)是爱尔兰的提摩太·G·迪南(Tithomy G. Dinan)等人在2013年所定义的名词,一般是指在适量食用后,透过与肠道菌的相互作用促进精神心理健康的活性微生
  • 博茨瓦纳经济博茨瓦纳经济在1966至1999年间每年平均增长9%左右,私人机构就业人数每年平均增长约10%。21世纪初,经济陷入停滞,直到2010年代初国内生产总值增长率才恢复到6-7%。博茨瓦纳是世界
  • 王永志王永志(1932年11月17日-),出生于辽宁省昌图县,航天技术专家,中国载人航天工程的开创者和学术技术带头人。40多年来在中国战略火箭、地地战术火箭以及运载火箭的研制工作中做出了突