5.3.1.1 朗之万方程

文档摘要

5.3.1.1 朗之万方程 5.3.1.1 朗之万方程：当随机力不“随机”——一个被忽略的数值采样陷阱与它的三重修复术你有没有试过用朗之万动力学模拟一个纳米颗粒在微流控通道中的布朗运动，结果发现它在零外力下却持续向右漂移？有没有调过 LangevinIntegrator 的 gamma（摩擦系数）从 1 ps⁻¹ 调到 100 ps⁻¹，却发现扩散系数 D 不但没按 $D = \frac{kB T}{\gamma}$ 衰减，反而先降后升、甚至发散？有没有在 GPU 上跑完百万步轨迹后，用 MSD（均方位移）拟合出的 $D{\text{eff}}$ 比理论值低 37% ——而同事在 CPU 上复现同一参数却只差 2%？这不是你的模型错了。不是物理理解有偏差。甚至不是代码有 bug。