关于Brownian motion在势场奇点聚集及HMC模拟与Gibbs测度偏差的技术咨询

嘿，这个问题戳中了蒙特卡洛模拟里一个挺容易踩坑的点——奇异势场下的遍历性偏差，我来帮你梳理清楚：

核心矛盾：理论Gibbs测度 vs 模拟的伪平衡 #

首先得明确：你预期的Gibbs测度，本质上是和$\exp(-\beta V(x))$成正比的（$\beta$是逆温度，$V(x)$是你的势场）。如果势场在原点有吸引性奇点（比如$V(x) \to -\infty$当$x \to 0$），只要$\exp(-\beta V(x))$在整个环面上的积分是有限的，这个Gibbs测度就是良定义的——它会在奇点附近有更高的密度，但绝对不会是“粒子几乎全时间待在原点”这种极端情况。你看到的模拟结果，大概率是HMC数值离散化引入的伪平衡，而非真实的连续Brownian运动的平衡分布。

连续情形的数学本质 #

对于连续Brownian运动在势场$V(x)$下的平衡态，对应的是Fokker-Planck方程的稳态解：
$$\rho(x) \propto \exp(-\beta V(x))$$
这个结论是严格的——只要$\exp(-\beta V(x))$可积，真实的平衡分布就完全由这个密度决定。举个具体例子：假设在2D环面上，原点附近势场是$V(r) \sim -\alpha \log r$（$\alpha>0$），那Gibbs测度密度就是$\propto r^{\beta \alpha}$，积分在环面上有限的条件是$\beta \alpha > -1$，这时候原点附近的密度确实更高，但粒子仍然会在整个环面上游走，只是在奇点附近停留时间更长而已。

为什么HMC会“卡”在奇点？ #

HMC的核心依赖是哈密顿动力学的辛守恒性和遍历性，但当势场有强奇点时，数值离散化的问题会被放大：

• 当粒子靠近原点时，势能急剧下降，根据总能量守恒，动能会急剧上升，粒子的运动速度会变得极大；
• 你的离散化步长如果固定，要么步长太大，直接跳过奇点区域（导致采样不到高密区），要么步长太小，数值积分的误差会让总能量不再严格守恒——粒子会逐渐损失动能，最终被困在势能最低的奇点附近，形成“伪平衡”。
  简单说：HMC的离散流在奇异势场下破坏了遍历性，没法遍历整个状态空间，自然得不到真实的Gibbs测度。

相关研究与解决方案 #

这种奇异势场下的蒙特卡洛偏差早就被研究过了，属于遍历性破缺的一类典型场景：

• 针对HMC的改进：可以用自适应步长HMC，在奇点附近自动缩小步长，或者使用更高阶的辛积分器（比如Symplectic Euler的高阶版本），减少能量守恒的误差；
• 势能正则化：给势场加一个小的正则项，比如$V_\epsilon(x) = V(x) + \epsilon |x|^2$，让原点处的势能不再是负无穷，先模拟$V_\epsilon$下的分布，再取$\epsilon \to 0$的极限；
• 对照验证：用Metropolis-Hastings（MH）算法做对照——MH对奇异势场的鲁棒性更好（虽然效率低），如果MH能得到符合Gibbs测度的结果，就能确认是HMC的数值问题；
• 理论层面：关于离散哈密顿系统在奇异势场下的遍历性，有不少研究证明，当势能奇点足够强时，离散化的流会失去遍历性，导致采样偏差。

怎么预测真实的平衡分布？ #

1. 严格计算Gibbs测度：先验证$\exp(-\beta V(x))$在环面上的积分是否有限，如果有限，归一化后的密度就是真实的平衡分布；
2. 修正模拟方案：调整HMC的参数（步长、积分器阶数），或者换用鲁棒性更强的采样算法，直到模拟结果逼近理论密度；
3. 小范围验证：可以先在奇点附近的小区域内计算粒子的停留概率，和理论密度的积分结果对比，确认是否一致。

备注：内容来源于stack exchange，提问作者Jesse