[論文レビュー] Convergence rates of supercell calculations in the reduced Hartree-Fock model

David Gontier, Salma Lahbabi|arXiv (Cornell University)|Jul 1, 2015

Spectral Theory in Mathematical Physics参考文献 13被引用数 11

ひとこと要約

本稿は、完全な結晶における縮小ハートリー・フォック（rHF）モデルのスーパセル計算について、周期的rHFモデルへの指数的収束速度を確立する。ブリユアンゾーン上での積分とそのリーマン和の差を分析することにより、結晶が絶縁体または半導体である場合、単位格子当たりのエネルギーおよび電子密度が、超胞サイズLに関して指数的かつ高速に収束することを証明する。誤差の上限は、正の定数Cとαに対してCe−αLの形をとる。

ABSTRACT

This article is concerned with the numerical simulations of perfect crystals. We study the rate of convergence of the reduced Hartree-Fock (rHF) model in a supercell towards the periodic rHF model in the whole space. We prove that, whenever the crystal is an insulator or a semi-conductor, the supercell energy per unit cell converges exponentially fast towards the periodic rHF energy per unit cell, with respect to the size of the supercell.

研究の動機と目的

• スーパセルrHF計算が周期的rHFモデルに指数的収束するという数値的観察を、数学的に厳密に裏付けること。
• 理論的解析で知られていたO(L−1)収束率と、数値的に観測された指数的収束の間のギャップを埋めること。
• 単位格子当たりのエネルギーおよび電子密度の両方について、L∞(R³)ノルムにおける定量的指数的収束速度を確立すること。
• 計算材料科学分野で広く用いられているが、未だ証明のない仮定、すなわち均一なブリユアンゾーンサンプリングが絶縁体および半導体系において小さな誤差をもたらすという仮定の数学的基盤を提供すること。

提案手法

• スーパセルrHFエネルギーを、ブリユアンゾーン上での周期的rHFエネルギー積分のリーマン和近似に再定式化する。
• エネルギー差の被積分関数が複素数ストリップ上で解析的かつ周期的であることを証明し、複素解析の道具を用いることを可能にする。
• 密度差の計算におけるリゾルベント作用素のトレースクラスノルムを制御するために、Kato–Seiler–Simonの不等式を適用する。
• Bloch–Floquet理論を用いて、ハミルトニアンおよび密度作用素をブリユアンゾーン上のファイバー束に分解する。
• スペクトルギャップと解析性に基づき、リゾルベント作用素（Bper₂およびBL₂）に対して一様なバインディングを確立し、指数的減衰推定に至る。
• CoulombカーネルG₁の滑らかさの性質を活用して、L²ノルムからL∞(R³)ノルムへの収束をブートストラップする。

実験結果

リサーチクエスチョン

• RQ1 なぜ、従来の理論的境界がO(L−1)にとどまっているにもかかわらず、スーパセルrHF計算の数値誤差は超胞サイズLに関して指数的に減少するのか？
• RQ2 絶縁体および半導体系において、スーパセルrHFエネルギーの単位格子当たりの値が周期的rHFエネルギーに指数的に収束することは、数学的に厳密に証明可能か？
• RQ3 スーパセルモデルにおける電子密度は、L∞(R³)ノルムにおいて、周期的rHF密度に指数的に高速に収束するか？
• RQ4 収束速度は、系のスペクトルギャップおよびブローチバンドの解析性にどのように依存するか？
• RQ5 平均場ハミルトニアンの収束が、エネルギーおよび密度と同様に指数的レートに従うことを示せるか？

主な発見

• スーパセルrHFモデルにおける単位格子当たりのエネルギーは、あるC > 0およびα > 0に対して、Ce−αLの形の誤差上限をもって周期的rHFエネルギーに指数的かつ高速に収束する。この収束は系のサイズに依存しない。
• スーパセルモデルにおける電子密度は、L∞(R³)ノルムにおいて、周期的rHF密度に同様にCe−αLの誤差上限をもって指数的かつ高速に収束する。
• 収束速度は、系のスペクトルギャップおよび複素平面におけるブローチバンドの解析性によって支配され、これによりフーリエ係数の減衰が保証される。
• 平均場ハミルトニアンHLは、作用素ノルムにおいて、周期的ハミルトニアンH0に指数的レートCe−αLで収束する。
• ミニマックス原理により、スーパセルハミルトニアンの固有値が、周期的ハミルトニアンの固有値に同じ指数的レートで収束することが示唆される。
• 線形およびrHFモデルにおける結晶シリコンの数値シミュレーションにより、エネルギーおよび密度誤差が超胞サイズLの増加に伴い、予測された指数的減衰を確認した。