埃瓦尔德求和

✍ dations ◷ 2025-05-17 10:39:16 #电磁学,位势论,傅里叶分析,计算化学,理论化学

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

埃瓦尔德求和将相互作用势(英语:Interatomic potential)表示为两部分之和:

其中, φ s r ( r ) {\displaystyle \varphi _{sr}(\mathbf {r} )} )的长程部分是有限的,但可能有简易的数学形式,如高斯分布。该方法假设短程势容易求和,因此需要重点考虑的是长程势。由于使用了傅里叶级数,该方法将周期性边界条件作为假设,此周期性系统的重复单元称为原胞,选择一个原胞作为中央原胞作为参考,其余单元称为镜像。

长程力的能量是中央原胞的电荷与晶格所有电荷间相互作用能(英语:Interaction energy)之和,因此可以表示为原胞和晶格的电荷密度的双重积分:

其中原胞的电荷密度 ρ u c ( r ) {\displaystyle \rho _{uc}(\mathbf {r} )} 是中央原胞中位置 r k {\displaystyle \mathbf {r} _{k}} 上的电量 q k {\displaystyle q_{k}} 之和:

总电荷密度 ρ TOT ( r ) {\displaystyle \rho _{\text{TOT}}(\mathbf {r} )} 是原胞及其镜像电量 q k {\displaystyle q_{k}} 之和:

这里, δ ( x ) {\displaystyle \delta (\mathbf {x} )} 表示狄拉克δ函数, a 1 {\displaystyle \mathbf {a} _{1}} a 2 {\displaystyle \mathbf {a} _{2}} a 3 {\displaystyle \mathbf {a} _{3}} 表示晶格矢量, n 1 {\displaystyle n_{1}} n 2 {\displaystyle n_{2}} n 3 {\displaystyle n_{3}} 的范围为所有整数。总电荷密度 ρ TOT ( r ) {\displaystyle \rho _{\text{TOT}}(\mathbf {r} )} 可以表示为 ρ u c ( r ) {\displaystyle \rho _{uc}(\mathbf {r} )} 与晶格函数 L ( r ) {\displaystyle L(\mathbf {r} )} 的卷积:

由于 ρ TOT ( r ) {\displaystyle \rho _{\text{TOT}}(\mathbf {r} )} 为卷积,其傅里叶变换为一个积:

其中晶格函数的傅里叶变换是狄拉克δ函数的另一个和:

其中定义倒空间矢量为 b 1   = d e f   a 2 × a 3 Ω {\displaystyle \mathbf {b} _{1}\ {\stackrel {\mathrm {def} }{=}}\ {\frac {\mathbf {a} _{2}\times \mathbf {a} _{3}}{\Omega }}} (周期性排列),其中 Ω   = d e f   a 1 ( a 2 × a 3 ) {\displaystyle \Omega \ {\stackrel {\mathrm {def} }{=}}\ \mathbf {a} _{1}\cdot \left(\mathbf {a} _{2}\times \mathbf {a} _{3}\right)} 为中心原胞的体积(几何形状通常为平行六面体), L ( r ) {\displaystyle L(\mathbf {r} )} L ~ ( k ) {\displaystyle {\tilde {L}}(\mathbf {k} )} 为实函数和偶函数。

为了简洁起见,定义有效单粒子势能:

因为其亦为卷积,其傅里叶变换是一个积:

其中定义了傅里叶变换:

现在,长程力的能量可以表示为单个电荷密度的积分:

使用帕塞瓦尔定理,能量亦可于倒空间中求和:

其中 k = m 1 b 1 + m 2 b 2 + m 3 b 3 {\displaystyle \mathbf {k} =m_{1}\mathbf {b} _{1}+m_{2}\mathbf {b} _{2}+m_{3}\mathbf {b} _{3}} 是最终的和值。

计算出 ρ ~ u c ( k ) {\displaystyle {\tilde {\rho }}_{uc}(\mathbf {k} )} 后, k {\displaystyle \mathbf {k} } 的和值或积分是显然的,可以很快地收敛。不能收敛的最常见原因是原胞不太明确,其必须为电中性,以避免无穷大的和。

在计算机普及前,埃瓦尔德求和是理论物理的理论。然而,自20世纪70年代以来,埃瓦尔德求和在粒子系统的计算机模拟中被广泛使用,尤其是遵守平方反比定律的粒子相互作用,如重力和静电力。最近,粒子网格埃瓦尔德方法也用于计算兰纳-琼斯势的 r 6 {\displaystyle r^{-6}} 部分,以消除截断产生的伪影(英语:Artifact (error))。其应用包括等离子体、星系及分子的模拟。

在粒子网格埃瓦尔德方法中,和标准埃瓦尔德求和相同,相互作用势(英语:Interatomic potential)被分为两部分 φ ( r )   = d e f   φ s r ( r ) + φ r ( r ) {\displaystyle \varphi (\mathbf {r} )\ {\stackrel {\mathrm {def} }{=}}\ \varphi _{sr}(\mathbf {r} )+\varphi _{\ell r}(\mathbf {r} )} ,其基本思想是用实空间中短程力的直接求和 E s r {\displaystyle E_{sr}} (粒子部分),及倒空间中长程力的求和(埃瓦尔德部分),代替点粒子间相互作用的能量的直接求和:

其中 Φ ~ r {\displaystyle {\tilde {\Phi }}_{\ell r}} ρ ~ ( k ) {\displaystyle {\tilde {\rho }}(\mathbf {k} )} 表示力和电荷密度的傅里叶变换。由于两个求和分别在实空间和倒空间中迅速收敛,它们可能被精确截断,且所需计算时间大幅减少。计算电荷密度的傅里叶变换 ρ ~ ( k ) {\displaystyle {\tilde {\rho }}(\mathbf {k} )} 可使用快速傅里叶变换,需在空间中的离散格子(英语:Lattice (discrete subgroup))上(即网格部分)估计电荷密度。

由于埃瓦尔德方法隐含的周期性假设,粒子网格埃瓦尔德方法于物理系统中的应用需施加周期性。因此,该方法最适合用于空间范围内可以模拟为无限的系统。在分子动力学模拟中,常构造可以无限平铺形成镜像的电中性原胞;然而,为了正确解释这种近似效应,这些镜像被重新并入原始模拟原胞中,这种整体效应被称为周期性边界条件。 想象一个单位立方体,上表面与下表面有效接触,右侧面与左侧面有效接触,前表面与后表面有效接触。因此,原胞的尺寸必须足够大,以避免两个接触面间不正确的运动相关性,但仍需足够小以便计算。短程力与长程力间截断的定义也可以引入伪影(英语:Artifact (error))。

电荷密度对网格的限制,使得粒子网格埃瓦尔德方法对电荷密度或势函数平滑变化的系统更有效。利用快速多极子方法(英语:Fast multipole method)可以更有效地处理局部系统或电荷密度波动较大的系统。

极性晶体(即原胞中具有净偶极子 p u c {\displaystyle \mathbf {p} _{uc}} 的晶体)的静电能为条件收敛,即取决于求和顺序。例如,若中央原胞的偶极与不断增加的立方体上的原胞偶极相互作用,则其能量收敛值并不会与考虑不断增大的球面时相等。大致来说,这种条件收敛是因为在半径为 R {\displaystyle R} 的壳上的偶极子数约为 R 2 {\displaystyle R^{2}} ;偶极-偶极相互作用的强度约为 1 R 3 {\displaystyle {\frac {1}{R^{3}}}} ;而两者相乘的结果是发散的调和级数 n = 1 1 n {\displaystyle \sum _{n=1}^{\infty }{\frac {1}{n}}}

这看似令人惊讶的结果并不与现实晶体能量有限的事实相违背,因为现实晶体并非无限,具有特定边界。具体而言,极性晶体的边界的有效表面电荷密度为 σ = P n {\displaystyle \sigma =\mathbf {P} \cdot \mathbf {n} } ,其中 n {\displaystyle \mathbf {n} } 为表面法矢量, P {\displaystyle \mathbf {P} } 为单位体积的净偶极矩。则中央原胞之偶极子与表面电荷密度 σ {\displaystyle \sigma } 的相互作用能 U {\displaystyle U} 可写为:

其中, p u c {\displaystyle \mathbf {p} _{uc}} V u c {\displaystyle V_{uc}} 分别为原胞的净偶极矩和体积, d S {\displaystyle dS} 为晶面上的无穷小区域, r {\displaystyle \mathbf {r} } 为中央原胞到无穷小区域的矢量。此公式来自于对能量 d U = p u c d E {\displaystyle dU=-\mathbf {p} _{uc}\cdot \mathbf {dE} } 积分,其中 d E {\displaystyle d\mathbf {E} } 表示无穷小电场,由无穷小的表面电荷 d q = d e f σ d S {\displaystyle dq{\stackrel {\mathrm {def} }{=}}\sigma dS} 产生(库仑定律):

负号来自于 r {\displaystyle \mathbf {r} } 的定义,其指向电荷方向为正方向。

埃瓦尔德求和由德国物理学家保罗·彼得·埃瓦尔德于1921年发表,用于确定离子晶体的静电能及马德隆常数。

不同的埃瓦尔德求和具有不同的时间复杂度。直接求和的时间复杂度为 O ( N 2 ) {\displaystyle O(N^{2})} ,其中 N {\displaystyle N} 为系统中原子数。粒子网格埃瓦尔德方法的时间复杂度为 O ( N log N ) {\displaystyle O(N\,\log N)}

相关

  • 舌形动物类舌形动物类(学名:Pentastomida),又称五口动物类,是动物界的一个分类。属于无脊椎动物,介乎于环节动物和节肢动物之间,均为寄生性动物。其特征是体软、扁而长,无色,透明,无足。身长从数
  • 布陀布陀(英语:Buto),又称勒托,是古埃及的城市,在今埃及谢赫村省境内,位于迪斯沃克附近的Ibtu村或者Abtu村,在亚历山大港以东95公里处,是旧石器时代到进入第一王朝之前下埃及的一个首要文
  • 珍妮特·纳波利塔诺珍妮特·纳波利塔诺(英语:Janet Napolitano,1957年11月29日-),美国政治家和律师,现是第20任加州大学校长。她曾是第21任亚利桑那州州长,是该州的第三位女性州长。2006年至2007年期间
  • 圣何塞地震圣何塞地震(英语:San Jose Earthquakes)是美国职业足球大联盟足球队,位于加利福尼亚州圣何塞。他们曾于1995年-2005年运作、2007年重新经营,并于2008年重新加盟大联盟。球队的主
  • 提康德罗加级提康德罗加级巡洋舰(英语:Ticonderoga-class cruiser)是美国海军所属第一种配备宙斯盾(AEGIS)系统的战舰,其特色为配备以AN/SPY-1舰用相控阵雷达为核心的整合式水面作战系统(宙斯盾
  • 彼得·毕克尔彼得·约翰·毕克尔 (Peter John Bickel,1940年-)是在加利福尼亚大学伯克利分校工作的著名美国统计学家 。毕克尔在横跨统计学和概率论的多个领域做出了重大贡献,这些领域包括:
  • 亚尔斯兰战记角色列表以下内容,记述著田中芳树所著作小说《亚尔斯兰战记》故事中的人物介绍列表。声为剧场版和OVA版/电视动画版。是指奉亚尔斯兰为主、被人称之为“解放王亚尔斯兰的十六翼将”的
  • 灯塔 (电影)《灯塔》(英语:)是一部2019年美国和加拿大合拍的黑白奇幻恐怖片,由罗柏·艾格斯执导并与弟弟麦斯·艾格斯(Max Eggers)共同撰写剧本。该片由威廉·达福和罗伯特·帕丁森主演,故事讲
  • 爱之初《爱之初》是2020年首播的中国大陆电视剧。该剧根据汪一洋小说《洋嫁》改编,讲述了俞飞鸿饰演的谢桥,在海外经历辛酸打拼、爱情浮沉后最终选择回国成就梦想与幸福的暖心故事
  • 裙带菜酒 (性文化)裙带菜酒(わかめ酒)是日本艺妓座敷艺的一种,和女体盛、花电车等一样是日本传统性文化的一部分。原本是烟花巷里像在正月这样的庆典中艺妓和风韵之徒之间的一种游戏。女性裸露身