エバルトの方法のソースを表示

{{出典の明記|date=2015年9月}}
'''エバルトの方法'''（エバルトのほうほう、{{lang-en-short|Ewald method}}）は、[[分子動力学法]]、[[量子化学的手法]]、[[バンド計算]]手法などで、[[単位胞]]（または[[スーパーセル法|スーパーセル]]）内の[[原子核]]（または[[イオン芯]]）同士の[[クーロン相互作用]]を効率良く計算する手法である。

[[実空間]]と[[逆格子空間]]のどちらでも発散してしまうクーロン相互作用を、実空間での計算が都合のよい部分と、逆格子空間での計算が都合のよい部分との2つに分けて別々に計算を行い、これら2つの計算結果の和が求めるべき答となる。

==具体例==
格子点<math>\boldsymbol{l}</math>におかれたイオンの作る静電ポテンシャル<math>\phi(\boldsymbol{r})</math>を例にとって説明する。<math>\phi(\boldsymbol{r})</math>は[[係数]]を省略すれば次の式で書ける。

:<math>\phi(\boldsymbol{r}) = \sum_{\boldsymbol{l}} \frac{1}{\left|\boldsymbol{r}-\boldsymbol{l}\right|}</math>

静電ポテンシャル<math>\phi(\boldsymbol{r})</math>を、関数<math>f(\boldsymbol{r})</math>を用いて、次のように分割する。

:<math>\phi(\boldsymbol{r}) = \sum_{\boldsymbol{l}} \frac{f(\boldsymbol{r})}{\left|\boldsymbol{r}-\boldsymbol{l}\right|} + \sum_{\boldsymbol{l}} \frac{1-f(\boldsymbol{r})}{\left|\boldsymbol{r}-\boldsymbol{l}\right|}</math>

ここで、<math>f(\boldsymbol{r})</math>として短距離で0に[[収束]]する関数をうまく選び、第1項の和が実空間で短距離で0に収束し、第２項の和が逆格子空間で短距離で0に収束すると計算上都合がよい。
実際には、<math>f(\boldsymbol{r})</math>として[[相補誤差関数]] erfc がよく使われる。''G''を任意の定数として<math>f(\boldsymbol{r})=\operatorname{erfc}(G\left|\boldsymbol{r}-\boldsymbol{l}\right|)</math>を、<math>1-f(\boldsymbol{r})</math>は[[誤差関数]]の定義式を代入すると

:<math>\phi(\boldsymbol{r}) = \sum_{\boldsymbol{l}} \frac{1}{\left|\boldsymbol{r}-\boldsymbol{l}\right|}\operatorname{erfc}(G\left|\boldsymbol{r}-\boldsymbol{l}\right|) + \sum_{\boldsymbol{l}} \frac{1}{\left|\boldsymbol{r}-\boldsymbol{l}\right|}\frac{2}{\sqrt{\pi}}\int_0^{G\left|\boldsymbol{r}-\boldsymbol{l}\right|} e^{-t^2}\,\mathrm dt</math>

第2項で<math>t=\left|\boldsymbol{r}-\boldsymbol{l}\right|\rho</math>とおいて変数変換する。

:<math>\phi(\boldsymbol{r}) = \sum_{\boldsymbol{l}} \frac{1}{\left|\boldsymbol{r}-\boldsymbol{l}\right|}\operatorname{erfc}(G\left|\boldsymbol{r}-\boldsymbol{l}\right|) + \int_0^G \sum_{\boldsymbol{l}} \frac{2}{\sqrt{\pi}} \exp(-|\boldsymbol{r}-\boldsymbol{l}|^2\rho^2) \,\mathrm d\rho </math>

第２項の被積分関数は、格子の周期を持つ<math>\boldsymbol{r}</math>の関数であるから、[[フーリエ級数]]に展開できる。単位格子の体積を<math>v_{\mathrm{cell}}</math>とおくと、次のように第２項を展開できる。

:<math>\begin{align}
 \phi(\boldsymbol{r}) & = \sum_{\boldsymbol{l}} \frac{1}{\left|\boldsymbol{r}-\boldsymbol{l}\right|}\operatorname{erfc}(G\left|\boldsymbol{r}-\boldsymbol{l}\right|) + \int_0^G \sum_{\boldsymbol{g}} \frac{2\pi}{v_{\mathrm{cell}}} \frac{1}{\rho^3} \exp(-|\boldsymbol{g}|^2/4\rho^2) \exp(i\boldsymbol{g}\cdot\boldsymbol{r}) \,\mathrm d\rho \\
& = \sum_{\boldsymbol{l}} \frac{1}{\left|\boldsymbol{r}-\boldsymbol{l}\right|}\operatorname{erfc}(G\left|\boldsymbol{r}-\boldsymbol{l}\right|) + \frac{4\pi}{v_{\mathrm{cell}}} \sum_{\boldsymbol{g}} \frac{\exp(-\left|\boldsymbol{g}\right|^2/4G^2)}{|\boldsymbol{g}|^2} \exp(i\boldsymbol{g}\cdot\boldsymbol{r})
\end{align} </math>

よって、<math>\operatorname{erfc}(x)</math>や<math>e^{-x}</math>といった、<math>x</math>について速やかに0に収束する関数が現れる形に変形することができた。
後は、各項がそれぞれ<math>\boldsymbol{l}</math>と<math>\boldsymbol{g}</math>に対して速く収束するように適当なGの値を選べば効率よく計算できる。

==粒子・メッシュ・エバルト (PME) 法==
エバルト法は、[[計算機]]の出現よりずっと以前に、[[理論物理学]]における手法として開発された。しかしながら、エバルト法は1970年代以降、粒子系の{{仮リンク|コンピュータシミュレーション|en|Computer simulation}}、特に[[重力]]や[[静電気学]]といった[[逆2乗の法則|逆2乗]]力を介して相互作用する粒子系において広範に使用されている。最近、PME法は打ち切りによるアーティファクトを除去するために[[レナード-ジョーンズ・ポテンシャル]]の<math>r^{-6}</math>部分の計算にも使用されている<ref>{{Citation |last=Di Pierro |first=M. |last2=Elber |first2=R. | last3=Leimkuhler |first3=B. |     title=A Stochastic Algorithm for the Isobaric-Isothermal Ensemble with Ewald Summations for all Long Range Forces. | journal=[[Journal of Chemical Theory and Computation]] | year=2015| doi=10.1021/acs.jctc.5b00648}}.</ref>。PME法は[[プラズマ]]、[[銀河]]、[[分子]]のシミュレーションに応用されている。

[[粒子・メッシュ法]]では、標準のエバルト和と同じく、包括的相互作用ポテンシャルが2つの項へと分離される。
:<math>\varphi(\boldsymbol{r}) \ \stackrel{\mathrm{def}}{=}\  \varphi_{\mathrm{sr}}(\boldsymbol{r}) + \varphi_{\mathrm{lr}}(\boldsymbol{r})</math>.
粒子・メッシュ・エバルト和の基本的考えは、点粒子間の相互作用エネルギーの直接和
:<math>
E_\text{TOT} = \sum_{i,j} \varphi(\boldsymbol{r}_{j} - \boldsymbol{r}_i) = E_{\mathrm{sr}} + E_{\mathrm{lr}}
</math>
を実空間における短距離ポテンシャルの直接和<math>E_{\mathrm{sr}}</math>（すなわち'''粒子・メッシュ・エバルト'''の'''粒子'''部分）
:<math>
E_{\mathrm{sr}} = \sum_{i,j} \varphi_{\mathrm{sr}}(\boldsymbol{r}_j - \boldsymbol{r}_i)
</math>
と長距離部分のフーリエ空間における和
:<math>
E_{\mathrm{lr}} = \sum_{\boldsymbol{k}} \tilde{\Phi}_{\mathrm{lr}}(\boldsymbol{k}) \left| \tilde{\rho}(\boldsymbol{k}) \right|^2
</math>
へと置き換えることである。<math>\tilde{\Phi}_{\ell r}</math>および<math>\tilde{\rho}(\boldsymbol{k})</math>は[[ポテンシャル]]および[[電荷密度]]のフーリエ変換を表す（すなわち'''エバルト'''部分）。どちらの和もそれぞれの空間（実空間ならびにフーリエ空間）において素早く収束するため、精度の損失がほとんどなく打ち切ることができ、必要な計算時間を大きく改善することができる。電荷密度場のフーリエ変換<math>\tilde{\rho}(\boldsymbol{k})</math>を効率的に評価するため、[[高速フーリエ変換]]が用いられる。高速フーリエ変換では、密度場は空間中の離散格子（すなわち'''メッシュ'''部分）上で評価される必要がある。

エバルト和はポテンシャルの周期性を仮定する。PME法の物理系への適用にあたり，ポテンシャルの周期的な対称性が必要となる。そのため，この手法は空間的に無限に広がる系(バルク固体など)のシミュレーションに適している。[[分子動力学]]シミュレーションでは、これは電荷が中性の単位セルを無限に並べることによって通常達成される。しかしながら、この近似の効果を適切に説明するために、無限に続くセルは元のシミュレーションセルへと再取り込みされる。これは[[周期的境界条件]]と呼ばれる。

密度場のメッシュへの制限は、密度の変動が「滑らか」な系（連続ポテンシャル関数を持つ系）についてPME法をより効率的にしている。局在した系、または密度の揺らぎが大きな系については、GreengardとRokhlinの{{仮リンク|高速多重極法|en|Fast multipole method}}を用いてより効率的に扱うことができる。

== 脚注 ==
{{reflist}}

==関連項目==
*[[エバルト項]]
*[[マーデルングエネルギー]]
*[[計算物理]]
*[[第一原理バンド計算]]

{{DEFAULTSORT:えはるとのほうほう}}
[[Category:固体物理学]]
[[Category:計算物理学]]
[[Category:理論化学]]