【计算材料学-从算法原理到代码实现】视频教程 | 3.4.3.2_Ewald求和方法
Ewald核心想法:给每个点电荷“加并减”一团高斯屏蔽电荷,把库仑相互作用拆成实空间的短程部分和倒空间的平滑长程部分。
计算步骤:
在实空间只求和短程项(衰减很快,截断即可)。
在倒空间对平滑长程项做傅里叶求和(少量波矢就收敛)。
加上自能校正和必要的背景/零频项处理。
可调参数:分裂参数控制“实空间 vs 倒空间”的工作量分配;同时选择合适的实空间截断和倒空间截断来达成精度。
适用条件:体系整体必须电中性(或加入均匀背景);天然适配周期性边界条件。
优点:两部分都快速收敛、精度可控;能量与力一致。
扩展
立即观看