扩散矩阵平方根是否唯一？朗之万方程中A的定义与选择问题

先直接给结论：式(1)不能唯一确定矩阵$A$。下面我来拆解原因，以及说说怎么选符合物理需求的$A$：

为什么$A$不唯一？ #

迁移率矩阵$M_{ij}^{ab}$作为扩散矩阵的正比形式，本质是对称正定矩阵——毕竟扩散过程的系数必须正定，否则会出现物理上不可能的负扩散行为。对于任意对称正定矩阵$M$，满足$M = A_{ik}^{ac}A_{jk}{bc}$（对应你的式(1)的指标收缩形式）的矩阵$A$不止一个：

• 假设$A$是一个满足条件的平方根，那么对于任意正交矩阵$O$（满足$O^T O = I$，也就是转置等于逆矩阵），矩阵$AO$也会满足：
  $$(AO)(AO)^T = A O O^T A^T = A I A^T = A A^T = M$$
  这就意味着存在无穷多个符合要求的$A$。
• 举个简单的例子：如果$M$是2x2单位矩阵，那么$A$可以是单位矩阵，也可以是任意旋转矩阵、反射矩阵，这些都满足$A A^T = I$。

如何选择“正确的”$A$？ #

这里的“正确”其实是指符合物理约束且便于计算/分析的选择，核心要满足朗之万方程的涨落耗散定理——也就是随机力的关联函数必须严格对应迁移率矩阵。常见的选择方式有两种：

1. 对称正定平方根 #

因为$M$是对称正定的，它存在唯一的对称正定平方根$A$，满足$A = A^T$且$A2 = M$。这个平方根可以通过$M$的特征分解直接得到：
假设我们对$M$做特征分解：$M = V \Lambda V^T$，其中$V$是正交的特征向量矩阵，$\Lambda$是对角矩阵（对角元是$M$的正特征值），那么对称平方根就是：
$$A = V \Lambda^{1/2} V^T$$
这种选择的优势是形式对称，物理上直观，而且数值计算稳定——大多数科学计算软件都有直接计算对称正定矩阵平方根的工具，非常适合数值模拟场景。

2. 基于特征模式的平方根 #

如果我们希望随机力的分量对应于扩散过程的本征模式（每个模式的涨落相互独立），可以选择$A = V \Lambda^{1/2}$。此时，随机力$\sqrt{2k_B T} A \xi(t)$的各个分量就对应$M$的特征模式，每个模式的涨落幅度由$\sqrt{\lambda_i}$（$\lambda_i$是$M$的特征值）决定。这种分解在分析多粒子扩散的主导行为、模式耦合时特别有用。

核心原则 #

无论选哪种$A$，只要保证随机力的关联函数满足：
$$\langle f_i^a(t) f_k^c(t') \rangle = 2k_B T M_{ik}^{ac} \delta(t-t')$$
就完全符合物理要求。实际研究中，你可以根据具体问题场景选最方便的形式——比如做大规模数值模拟时对称平方根更易用，做模式分析时特征模式分解更合适。

内容的提问来源于stack exchange，提问作者Quantum spaghettification