径向分布函数理论计算

我第一次做 径向分布函数理论计算是写水模型的 method 部分，审稿人问”你的 g_OO(r) 第一峰对应几配位”，我脱口而出 4，结果积分到第一谷算出来是 4.2，才意识到 g(r) 的配位数要靠积分而不是看峰高。那次让我把 径向分布函数理论计算 从”画条曲线”升级成”读得懂曲线在说什么”。

一、为什么 径向分布函数理论计算先懂定义

g(r) = ρ(r)/ρ_bulk，是局部数密度相对平均密度的比值。r 很小处 g(r)=0（粒子不能重叠），第一峰对应最近邻壳层，峰谷之间积分得到配位数 n(r) = 4πρ∫₀ʳ g(r’)r’²dr’。我做 径向分布函数理论计算时，默认先把”看峰”改成”积分配位数”，因为峰高受温度和密度影响大，配位数才稳健。

二、核心原理：g(r) 从哪来

平衡系综下 g(r) 由总关联函数推导，理论求解靠 Ornstein-Zernike 方程配 closure（Percus-Yevick 或 HNC 超网链）。PY 适合硬球、HNC 适合带电体系，但都是近似。关键认知：g(r) 的振荡衰减行为刻画了短程有序，衰减越快有序越短程；液体只在第一二壳层有结构，长程回归 1。我用它判断”这个体系有没有明显配位壳层”。

三、关键技术要点：截断与归一化

径向分布函数理论计算里最容易被忽略的是截断球壳的体积归一化和边界修正。盒子太小、粒子数少时，远距壳层统计涨落大，g(r) 在长程处乱抖。我默认要求体系粒子数 > 500、盒子 > 3 倍相关长度。另一个坑是只用单帧算，统计不够，必须跨多帧平均。曾经粒子数 200 算出来的第三壳层峰是噪声，误判成”长程有序”。

四、实操流程：从理论到可解释量

我的标准动作：明确研究对象（同类/异类 g_AB）→ 选理论或模拟路线 → 多帧统计平均 → 积分求配位数与峰位 → 和实验 X 射线/中子散射比对。最容易被砍的是边界修正，但 径向分布函数理论计算 的长程可信度全靠它，我默认保留。

补一个实例：NaCl 水溶液 g_NaO(r)，第一峰位 2.4 Å、积分配位数 5.8，和实验水合数 6 一致；若只看峰高会误以为配位很强。另一个例子是合金短程有序参数，从 g_AB 偏离随机混合的程度能反推有序度。如果只画不积，这些定量结论全丢了。这之后我做 g(r)，配位数积分是标配输出，不再只贴峰。

五、常见踩坑：壳层交叠与多峰

径向分布函数理论计算常见错是把相邻峰当独立壳层，实际是同一壳层不同取向的贡献交叠；或异类 g_AB 第一峰和同类 g_AA 峰位重叠时没分开算，配位数被双计。另一个坑是温度效应，高温 g(r) 峰变矮变宽，容易被误读成”结构消失”。

六、复盘：g(r) 是液体结构的指纹

回过头看，径向分布函数理论计算最值钱的是用一根曲线定量描述”谁和谁挨多近、配几位”，它是液体和合金结构的指纹。我接项目默认先分清同类/异类、算配位数再解读。把定义、归一化、积分三件事做对，g(r) 才读得懂。被证明有用的，是愿意为审稿人的一个问题回去翻统计力学的踏实。下次有人给我一张 g(r)，我会先问——你配位数积到哪个 r 了。