[论文解读] Fundamental solution of Fokker - Planck equation

Igor A. Tanski|ArXiv.org|Jul 4, 2004

Statistical Mechanics and Entropy参考文献 7被引用 10

一句话总结

本文推导了在速度相关阻尼及位置与速度空间中联合扩散作用下，非相互作用粒子的三维福克-普朗克方程的精确基本解。通过傅里叶变换与特征线法，构造出具有显式时间依赖系数的高斯解，并通过对称性分析与直接代入法验证。该解描述了在联合扩散与阻尼力作用下，初始为狄拉克函数分布的演化过程。

ABSTRACT

Fundamental solution of Fokker - Planck equation is built by means of the Fourier transform method. The result is checked by direct calculation. Changes: missed factor in (29), (30), corrected (31), (35), removed former (36), added explanation on prolonged operators (45), (46), (47).

研究动机与目标

• 推导三维相空间中位置与速度坐标联合扩散的福克-普朗克方程基本解的精确解析表达式。
• 建模在速度成比例阻尼与随机扩散作用下，非相互作用粒子系统的时变演化。
• 为位置与速度空间中初始狄拉克函数分布提供闭式解。
• 通过对称性分析与直接代入法验证解的正确性。

提出的方法

• 对福克-普朗克方程应用傅里叶变换，将其转化为动量空间中的一个一阶线性偏微分方程。
• 使用特征线法求解变换后的偏微分方程，识别出三个守恒量（动量 px, py, pz）以及涉及速度动量与初始条件的三个附加不变量。
• 通过积分特征方程，在傅里叶空间中推导解，并将傅里叶变换的对数表示为二次型。
• 利用逆傅里叶变换，结合卷积性质与高斯函数的已知变换对，完成求解。
• 将基本解构造为一个平移后的狄拉克函数（代表初始粒子位置与速度）与一个时间依赖的高斯核的卷积。
• 通过半逆向方法验证解：假设解具有高斯形式，代入原偏微分方程，确认所推导的系数满足由此导出的常微分方程组。

实验结果

研究问题

• RQ1 在速度相关阻尼及位置与速度空间中联合扩散作用下，三维福克-普朗克方程的精确基本解是什么？
• RQ2 在联合扩散与阻尼作用下，相空间中初始为狄拉克函数分布的解如何随时间演化？
• RQ3 能否通过傅里叶变换与特征线法推导并验证该解？
• RQ4 对应于空间与速度位移的微分算子是否为福克-普朗克方程的对称性？
• RQ5 所推导的解在直接代入原偏微分方程后是否满足？

主要发现

• 基本解为相空间中的高斯函数，其系数显式依赖于阻尼系数 α 与扩散系数 k，且随时间变化。
• 解的形式为一个平移狄拉克函数（用于描述初始位置与速度）与一个时变高斯核的卷积。
• 空间与速度分量的演化包含一个与 (1 − e^−αt)/α 成正比的漂移项，反映了速度向零的弛豫以及由初始速度引起的空间位移。
• 指数中二次型矩阵的行列式为 D(t) = αt(1 − e^−2αt) − 2(1 − e^−αt)²/(2α⁴)，该量控制了解的宽度随时间的演化。
• 通过验证所推导的系数 A(t)、B(t)、C(t) 与 Q(t) 满足将高斯假设代入福克-普朗克方程后得到的常微分方程组，完成了解的验证。
• 证明了微分算子 ∂/∂x、∂/∂y、∂/∂z 以及 (1/α)(1 − e^−αt)∂/∂x + e^−αt ∂/∂vx（以及 y、z 方向的类似形式）为该偏微分方程的对称性，从而可将解从零初始条件推广至任意初始位置与速度。