埃瓦尔德求和

title: 推导.
content: 埃瓦尔德求和将相互作用势表示为两部分之和:  其中,formula_2表示实空间中和值快速收敛的短程势,formula_3表示倒空间中和值快速收敛的长程势。所有量(如"r")的长程部分是有限的,但可能有简易的数学形式,如高斯分布。该方法假设短程势容易求和,因此需要重点考虑的是长程势。由于使用了傅里叶级数,该方法将周期性边界条件作为假设,此周期性系统的重复单元称为原胞,选择一个原胞作为中央原胞作为参考,其余单元称为镜像。  长程力的能量是中央原胞的电荷与晶格所有电荷间相互作用能之和,因此可以表示为原胞和晶格的电荷密度的双重积分:  其中原胞的电荷密度formula_5是中央原胞中位置formula_6上的电量formula_7之和:  总电荷密度formula_9是原胞及其镜像电量formula_10之和:  这里,formula_12表示狄拉克δ函数,formula_13、formula_14、formula_15表示晶格矢量,formula_16、formula_17、formula_18的范围为所有整数。总电荷密度formula_9可以表示为formula_5与晶格函数formula_21的卷积:  由于formula_9为卷积,其傅里叶变换为一个积:  其中晶格函数的傅里叶变换是狄拉克δ函数的另一个和:  其中定义倒空间向量为formula_26(周期性排列),其中formula_27为中心原胞的体积(几何形状通常为平行六面体),formula_21和formula_29为实函数和偶函数。  为了简洁起见,定义有效单粒子势能:  因为其亦为卷积,其傅里叶变换是一个积:  其中定义了傅里叶变换:  现在,长程力的能量可以表示为单个电荷密度的积分:  使用帕塞瓦尔定理,能量亦可于倒空间中求和:  其中formula_35是最终的和值。  计算出formula_36后,formula_37的和值或积分是显然的,可以很快地收敛。不能收敛的最常见原因是原胞不太明确,其必须为电中性,以避免无穷大的和。
title: 粒子网格埃瓦尔德(PME)方法.
content: 在计算机普及前,埃瓦尔德求和是理论物理的理论。然而,自20世纪70年代以来,埃瓦尔德求和在粒子系统的计算机模拟中被广泛使用,尤其是遵守平方反比定律的粒子相互作用,如重力和静电力。最近,粒子网格埃瓦尔德方法也用于计算兰纳-琼斯势的formula_38部分,以消除截断产生的伪影。其应用包括等离子体、星系及分子的模拟。  在粒子网格埃瓦尔德方法中,和标准埃瓦尔德求和相同,相互作用势被分为两部分formula_1,其基本思想是用实空间中短程力的直接求和formula_40(粒子部分),及倒空间中长程力的求和(埃瓦尔德部分),代替点粒子间相互作用的能量的直接求和:  其中formula_44和formula_45表示力和电荷密度的傅里叶变换。由于两个求和分别在实空间和倒空间中迅速收敛,它们可能被精确截断,且所需计算时间大幅减少。计算电荷密度的傅里叶变换formula_45可使用快速傅里叶变换,需在空间中的离散格子上(即网格部分)估计电荷密度。  由于埃瓦尔德方法隐含的周期性假设,粒子网格埃瓦尔德方法于物理系统中的应用需施加周期性。因此,该方法最适合用于空间范围内可以模拟为无限的系统。在分子动力学模拟中,常构造可以无限平铺形成镜像的电中性原胞;然而,为了正确解释这种近似效应,这些镜像被重新并入原始模拟原胞中,这种整体效应被称为周期性边界条件。 想象一个单位立方体,上表面与下表面有效接触,右侧面与左侧面有效接触,前表面与后表面有效接触。因此,原胞的尺寸必须足够大,以避免两个接触面间不正确的运动相关性,但仍需足够小以便计算。短程力与长程力间截断的定义也可以引入伪影。  电荷密度对网格的限制,使得粒子网格埃瓦尔德方法对电荷密度或势函数平滑变化的系统更有效。利用快速多极子方法可以更有效地处理局部系统或电荷密度波动较大的系统。
title: 偶极子.
content: 极性晶体(即原胞中具有净偶极子formula_47的晶体)的静电能为条件收敛,即取决于求和顺序。例如,若中央原胞的偶极与不断增加的立方体上的原胞偶极相互作用,则其能量收敛值并不会与考虑不断增大的球面时相等。大致来说,这种条件收敛是因为在半径为formula_48的壳上的偶极子数约为formula_49;偶极-偶极相互作用的强度约为formula_50;而两者相乘的结果是发散的调和级数formula_51。  这看似令人惊讶的结果并不与现实晶体能量有限的事实相违背,因为现实晶体并非无限,具有特定边界。具体而言,极性晶体的边界的有效表面电荷密度为formula_52,其中formula_53为表面法向量,formula_54为单位体积的净偶极矩。则中央原胞之偶极子与表面电荷密度formula_55的相互作用能formula_56可写为:  其中,formula_47和formula_59分别为原胞的净偶极矩和体积,formula_60为晶面上的无穷小区域,formula_61为中央原胞到无穷小区域的向量。此公式来自于对能量formula_62积分,其中formula_63表示无穷小电场,由无穷小的表面电荷formula_64产生(库仑定律):  负号来自于formula_61的定义,其指向电荷方向为正方向。
title: 历史.
content: 埃瓦尔德求和由德国物理学家保罗·彼得·埃瓦尔德于1921年发表,用于确定离子晶体的静电能及马德隆常数。
title: 复杂度.
content: 不同的埃瓦尔德求和具有不同的时间复杂度。直接求和的时间复杂度为formula_67,其中formula_68为系统中原子数。粒子网格埃瓦尔德方法的时间复杂度为formula_69。

Summary: 埃瓦尔德求和(英语:Ewald summation),是一种计算分子的周期性系统中长程力(如静电力)的方法,以德国物理学家保罗·彼得·埃瓦尔德命名。埃瓦尔德求和最初用于计算离子晶体的电势能,现在用于计算化学中计算长程力。埃瓦尔德求和是泊松求和公式的特殊形式,用倒空间中的等效求和代替实空间中相互作用能的总和。埃瓦尔德求和将相互作用势分为短程力和无奇点的长程力两部分,短程力在实空间中计算,长程力用傅里叶变换计算。与直接求和相比,此方法的优势为能量能够快速收敛,这意味着此方法在计算长程力时具有较高的精度和合理的速度,是计算分子的周期性系统中长程力的标准方法。此方法需要分子系统的电中性,以准确计算总库仑力。

Sommation d'Ewald: La sommation d'Ewald (ou parfois somme d'Ewald) est une méthode de calcul des énergies d'interaction de systèmes périodiques (et particulier des cristaux), et tout particulièrement les énergies électrostatiques. La sommation d'Ewald est un cas particulier de la formule sommatoire de Poisson, avec le remplacement de la sommation des énergies d'interaction dans l'espace réel par une sommation équivalente dans un espace de Fourier. L'avantage de cette approche est la convergence rapide de la sommation dans l'espace de Fourier comparée à son équivalent dans l'espace réel lorsque les interactions se font à longue distance. Les énergies électrostatiques comprenant à la fois des termes d'interactions de courtes et de longues portées, il est très intéressant de décomposer le potentiel d'interaction en termes de courte portée - dont la sommation se fait dans l'espace réel - et de longue portée - dont la sommation se fait dans l'espace de Fourier. La méthode fut développée par Paul Peter Ewald en 1921 (voir références) afin de déterminer l'énergie électrostatique (et, par là, la constante de Madelung) des cristaux ioniques.