Ewald求和基础：库仑相互作用、Ewald分解、收敛参数
聊到周期性体系的静电相互作用，咱们得从最头疼的问题开始。
你想想看，分子动力学模拟里，每个原子都带着电荷。正负电荷之间有长程库仑作用，一算就牵扯到无限多个周期镜像。真把整个空间全算一遍？不可能。截断掉？误差又大得离谱。
Ewald求和就是来解决这个矛盾的。说白了，它把慢收敛的库仑求和拆成两个快收敛的部分：一个在实空间算，一个在倒易空间（也就是傅里叶空间）算。再补上自能修正和偶极修正，结果就精确了。
我早期做NAMD调参时，第一次看到Ewald参数的配置项其实也有点懵。后来摸清了门道，才发现它的逻辑非常清晰。
库仑相互作用的“远”与“近”
先看一个最简单的场景：两个电荷qᵢ和qⱼ，距离rᵢⱼ。库仑势能写成：
U_coul = (1 / 4πε₀) * Σ (qᵢ qⱼ) / rᵢⱼ求和要对所有周期镜像的原子对进行。问题是，1/r 衰减太慢了！你截断到10 Å，远端的贡献依然不可忽略。如果硬截断，体系能量就不守恒，温度都会飘。我见过有人调截断值随手写个8 Å，结果模拟跑着跑着就炸了。
Ewald分解：核心思路
Ewald 求和把1/r拆成两项：
1/r = erfc(β r) / r + erf(β r) / r其中β是Ewald分裂参数（也叫收敛参数）。等号右边第一项在实空间衰减很快，第二项在倒易空间衰减很快。这样我们就能：
- 实空间部分：只算近邻原子。用截断值rc截断，配合erfc函数快速收敛。
- 倒易空间部分：用傅里叶变换算，靠波矢求和。这部分打包了所有长程相互作用。
- 自能修正：把倒易空间里多算的“自己跟自己”的相互作用减掉。
总静电势能就等于这三项之和。嗯，这里要注意，第三项很小但它必须算，否则能量不对。
核心组件：一张图说清楚
下面这张图可以帮你理解Ewald求和的整体流程：
收敛参数β：一切的关键
β这个参数决定了实空间和倒易空间的分工。β越大，erfc(β r)衰减越快，实空间截断可以设得更小；但倒易空间里需要更多的波矢才能保证精度。反之亦然。
我调试过一个RNA水溶液体系，一开始β选了0.35 Å⁻¹，实空间截断设到12 Å，算起来很慢。后来把β提到0.45 Å⁻¹，截断降到9 Å，速度翻了近一倍，精度完全没掉。
具体来说：
| 参数 | 物理含义 | 调参影响 |
|---|---|---|
| β (Ewald分裂参数) | 控制实空间和倒易空间的分工 | β↑ → 实空间截断可减小，但倒易空间波矢数需增加 |
| rc (实空间截断) | 实空间计算的最大原子对距离 | rc↑ → 计算量↑，但实空间误差↓ |
| Kmax (最大波矢) | 倒易空间求和的波矢截断 | Kmax↑ → 精度↑，但FFT网格密度↑ → 计算量↑ |
PME 为何代替了原始 Ewald
原始的 Ewald 求和要用 O(N²) 的时间复杂度来算倒易空间波矢求和，N 一上去就扛不住了。PME 把电荷格点化到 FFT 网格上，用快速傅里叶变换来做，复杂度降到 O(N log N)。
我在项目中遇到过用原始 Ewald 算一个 10 万原子体系的情况，一个时间步要好几秒。换成 PME 后，同样的精度，速度提升了几十倍。说白了，没有 PME，现在的大规模模拟根本跑不动。
好，你现在应该对 Ewald 求和有了一个整体认识：为什么需要它、它怎么拆的、β 参数如何影响精度和效率。下一套组合拳——PME 的具体网格配置和 FFT 参数，我们紧接着就要攻克它。