Reaching the thermodynamic limit of periodic CCSD cohesive energies and band gaps with denser Brillouin zone sampling

本論文は、熱力学的極限への凝集エネルギーおよびバンドギャップの信頼性の高い外挿を可能にするために、稠密なブリルアンゾーンサンプリング（最大216 k点）を実現する周期的なCCSD理論のスケーラブルな分散メモリ実装を提示しており、これにより8種の半導体および絶縁体に対して、実験データと比較してそれぞれ約0.1–0.2 eVおよび0.4 eVの誤差で決定的なベンチマーク値を提供している。

要約

あなたは、巨大なダイヤモンドや塩の塊のような、巨大で完全な結晶の挙動を理解しようとしていると想像してください。現実の世界では、これらの結晶は数兆個もの原子を含んでおり、非常に巨大です。しかし、コンピュータで研究する場合、科学者たちは通常、問題を小さく扱いやすい「ミニ結晶（スーパーセル）」へと縮小し、それが無限に繰り返されていると仮定しなければなりません。

このミニ結晶というアプローチには問題があります。それは、大陸全体の天気を理解しようとして、たった一つの裏庭だけを見ているようなものです。これでは全体像を見逃してしまい、「有限サイズ誤差」という、サンプルが小さすぎるために生じる間違いを引き起こしてしまいます。

課題：「ピクセル化された」結晶
量子化学の世界では、サンプルの「サイズ」は、ブリルアンゾーンと呼ばれる数学的な空間において、どれだけの点をサンプリングするかによって決まります。これらの点は、デジタル画像のピクセルのようなものです。

• 低解像度（少ないピクセル）： 結晶の性質（原子がどれほど強く結合しているか、あるいは電子がギャップを飛び越えるのにどれだけのエネルギーが必要かなど）について、ぼやけた不正確な画像が得られます。
• 高解像度（多いピクセル）： 明瞭で正確な結晶の画像が得られます。

問題は、CCSD（単一および二重励起を用いた結合クラスター法）と呼ばれる黄金標準の手法を用いてこれらの性質を計算することが、非常にコストがかかるということです。それは、電卓で4K映画をレンダリングしようとするようなものです。CCSDは計算負荷が非常に高いため、これまでの研究では「低解像度」の画像（小さなピクセルグリッド）を用いることしかできず、高解像度の画像がどのようになるかを推測するしかありませんでした。これらの推測はしばしば大きな誤差を生じさせてきました。

解決策：超強力なコンピュータ・チーム
著者らは、非常に効率的な新しいソフトウェアプログラムを構築しました。これは、大量のワーカー（最大12個のコンピュータノード、各9体96コアで動作）のように機能します。このチームは、これまでよりもはるかに大きな「画像」を扱うことができます。

小さなピクセルグリッドを見る代わりに、彼らはブリルアンゾーンにおいて最大216点（$6 \times 6 \times 6$のグリッド）をサンプリングすることができました。これは、ぼやけたサムネイルから高精細な4K画像へとアップグレードすることに相当します。この新しい明瞭さによって、彼らは推測に頼ることなく、真の「熱力学的極限」、すなわち完璧で無限の結晶の挙動をようやく捉えることができたのです。

判明したこと： 「黄金標準」のベンチマーク
この高精細アプローチを用いて、チームは8つの一般的な材料（酸化マグネシウム、シリコン、ダイヤモンドなど）について、2つの主要な性質を計算しました。

1. 結合エネルギー： 結晶を個々の原子へと引き離すために必要なエネルギー。
2. バンドギャップ： 電子が電気を導くために飛び越えなければならないエネルギー「ギャップ」（本質的に、その材料が絶縁体か半導体かを示すもの）。

彼らは、自らの高精細な結果を実世界の実験と比較しました。

• 結合エネルギー： 彼らの予測は現実に非常に近く、通常、誤差はわずか 0.1〜0.2 eV でした。彼らの結果は、原子が結合する強さをわずかに過小評価する傾向がありました。
• バンドギャップ： 彼らの予測もかなり良好で、誤差は約 0.4 eV でしたが、彼らはギャップをわずかに過大評価する（その材料が実際よりも少し優れた絶縁体であると予測する）傾向がありました。

「間接的」なパズル
一部の材料は「間接的」なバンドギャップを持っており、これらは計算がより困難です。それは、直接は見ることができない2点間の距離を測定しようとするようなものです。著者らは、標準的な推測方法（外挿法）がここでは失敗し、ギャップを過小評価してしまうことを発見しました。彼らは巧妙な「複合的」戦略を開発しました。つまり、まず直接の経路を測定し、次に間接的な経路に対する補正を加えることで、より正確な結果を得る方法です。

二酸化チタンのテスト
この手法がより複雑な材料に対しても有効であることを証明するために、彼らはルチル型二酸化チタン（一般的な白い顔料および光触媒）に適用しました。彼らの計算は 4.17 eV のバンドギャップを予測しました。これは実験値（補正後で約3.9 eV）よりもわずかに高い値ですが、著者らは、この小さな誤差は手法の既知の限界と一致しており、完璧な精度を得るためにはさらに複雑な物理学（三重励起など）が必要である可能性を示唆しています。

結論
この論文は単に新しい数値を提供しているだけではありません。これは決定的なベンチマークを提供しています。大量のコンピュータ・チームを使用して「高精細」なデータを生成することで、著者らはCCSD法がどの程度機能するかについての、新しく信頼できる基準を確立しました。彼らは、この手法が非常に優れている一方で、現実の世界と比較した際に、約0.3〜0.4 eVの予測可能な「ぼやけ」（誤差）が依然として存在することを示しました。これにより、他の科学者たちが将来、同様の計算をどの程度信頼できるかを正確に知ることができるようになります。

技術概要

技術要約：周期的なCCSDにおける凝集エネルギーおよびバンドギャップの熱力学的極限への到達

問題提起
周期的な単一および二重励起を含む結合クラスター理論（CCSD）は、平均場理論である密度汎関数理論（DFT）に対して系統的な改善可能性を提供する、固体電子構造計算のための有望な非摂動的手法である。しかし、その熱力学的極限（TDL：有限サイズ誤差が排除される物理的に意味のある極限）への適用は、計算コストによって厳しく制限されてきた。周期的なCCSDのストレージおよび実行時間は、$O(n^4 N_k^3)$ および $O(n^6 N_k^4)$ （ここで $n$ は単位胞あたりの原子数、$N_k$ は$k$点数）としてスケールする。その結果、これまでの研究は、三次元固体に対して最大でも $4^3$ までの$k$点メッシュに限定されていた。この制限は、熱力学的極限への信頼できる外挿を妨げている。特にバンドギャップは、基底状態エネルギー（$N_k^{-1}$）よりも収束が遅い（$N_k^{-1/3}$）ためである。十分な密度の$k$点サンプリングにアクセスできないことは、凝集エネルギーおよびバンドギャップの予測において未定量化の有限サイズ誤差を生じさせ、CCSDレベルの理論における決定的なベンチマークの確立を阻害している。

手法
これらの制限に対処するため、著者らはPySCFフレームワークに基づいた、新しい最適化された分散メモリ実装の周期的なCCSDおよび方程式運動論的CCSD（EOM-CCSD）を開発し、適用した。主な手法の特徴は以下の通りである：

• 並列化と対称性： 本コードは分散メモリおよび共有メモリの両方の並列化を利用し、空間群対称性を活用して既約ブリルアンゾーンで計算を行う。これにより、ノードあたり1.5 TBのメモリを持つ最大12ノード（96コア）での計算を可能にしている。
• サンプリング密度： この実装により、ダブルゼータ（DZ）基底系では最大 $N_k = 6^3 = 216$、$TZ$基底系では$5^3$、$QZ$基底系では$4^3$ のブリルアンゾーン・サンプリングを近似なしで行うことができる。
• 外挿スキーム： 研究ではTDLに到達するために特定の外挿モデルを用いている：
  • 凝集エネルギーについては、ハートリー・フォック（HF）エネルギーを3パラメータモデル（$N_k^{-1} + N_k^{-2}$）で適合させ、相関エネルギー（MP2およびCCSD）には2パラメータモデル（$N_k^{-1}$）を用いる。
  • バンドギャップについては、3つのモデルを評価している：モデルA（主要項 $N_k^{-1/3}$）、モデルAB（第2次まで）、およびモデルABC（第3次までの $N_k^{-1}$ 項を含む）。
  • 「GW-EOM」スキームもテストされている。これは、EOM-CCSDと $G_0W_0@HF$ 計算の間の有限サイズ誤差係数の比例関係を仮定し、より小さな$k$メッシュを用いてTDLのバンドギャップを予測するものである。
• 間接ギャップのための複合的アプローチ： 間接バンドギャップ（イオン化ポテンシャルと電子親和力の和）は、伝導帯最小（CBM）に対する$k$メッシュのシフト要求があるため収束が困難であることを認識し、複合的アプローチを提案している。これは、計算を価電子帯最大（VBM）における直接ギャップと、CBM-V辺オフセットに分離することで、より安定した外挿を可能にするものである。
• 基底関数系とポテンシャル： 計算には相関一貫基底系（cc-pVXZ, $X=D, T, Q$）およびGTH擬ポテンシャルを使用している。基底関数の不完全性による誤差は、外挿（例：$X^{-3}$ 形式）またはDZ、TZ、QZの結果の比較を通じて補正される。

主要な貢献と結果
著者らは、8つの単純な半導体および絶縁体（C, Si, BN, BP, MgO, LiH, LiF, LiCl）および新しく適用されたルチル型 $\text{TiO}_2$ について、収束したCCSD凝集エネルギーおよびバンドギャップを報告している。

• 凝集エネルギー： 本研究は、約0.1 eVに収束したCCSD凝集エネルギーの決定的なベンチマーク数値を提供している。

  • CCSDは、実験値（零点振動エネルギー補正後）と比較して、平均して0.1–0.2 eV、凝集エネルギーを過小評価する傾向がある一方、MP2は過大評価する傾向がある。
  • CCSDの精度は、標準的なDFT汎関数（PBE, PBEsol, SCAN）と同等か、それ以上に高いことが判明した。
  • 擬ポテンシャルの誤差は、ほとんどの固体において一般に小さい（<0.05 eV）が、BNおよびCではやや大きい。

• バンドギャップ： 本論文は、信頼できるTDL外挿には、従来アクセス可能であったものよりも大幅に大きな$k$メッシュ（$N_k \ge 3^3$）が必要であることを確立している。

  • 外挿モデル： 単純な単一項外挿（モデルA）を小さなメッシュ（$N_k \le 4^3$）に適用すると、バンドギャップを一貫して0.3–0.5 eV過小評価する。$N_k=3^3$–$6^3$で適合させた3パラメータのモデルABCが、最も堅牢なリファレンスとして特定された。
  • GW-EOMスキーム： $G_0W_0$ を用いてEOM-CCSDの有限サイズ誤差を補正するGW-EOMアプローチは、混合的な信頼性を示している。一部の材料（MgO, LiCl, LiF）ではうまく機能するが、他の材料（BN, BP, Si, C）ではギャップを0.4–0.5 eV過大評価しており、誤差係数の比例関係という仮定が普遍的に成立しないことを示している。
  • 間接ギャップ： 複合的アプローチは、間接ギャップ（例：BN）の収束問題を解決することに成功した。具体的には、VBMにおける直接ギャップの外挿は安定しており、CBMオフセットは急速に収束することが明らかになった。
  • 精度： 最終的なTDL EOM-CCSDバンドギャップは、ゼロ点振動補正後の実験値に対して、平均絶対誤差（MAE）が0.39 eV（擬ポテンシャル使用時）および0.42 eV（全電子計算）である。この手法は通常、バンドギャップを過大評価する。
  • $\text{TiO}_2$： ルチル型 $\text{TiO}_2$ について、TDL EOM-CCSDによる予測は4.17 eVであり、実験値（ZPR補正後、~3.9 eV）を約0.3 eV過大評価している。

• 単一励起の性質： 著者らは、実験からの偏差が大きい材料（具体的にはMgOやLiHのような岩塩構造）は、偏差が小さい材料（>95%）と比較して、単一励起の性質が低い（$n_{1}^{IP} \approx 93-94\%$）傾向にあることを指摘しており、これはこれらの系において三重励起が必要である可能性を示唆している。

意義
本研究は、標準的な軌道に基づく周期的なCCSD計算の中で最大規模のものであり、有限サイズ誤差および基底関数系の誤差を解決することに成功し、熱力学的極限におけるCCSDの凝集エネルギーおよびバンドギャップに関する最初の信頼できるベンチマークを提供した。本研究は、CCSDは凝集エネルギーに対しては高い精度を提供する一方で、バンドギャップへの適用は外挿スキームと$k$点密度に対して敏感であることを示している。著者らは、実用的な計算が $N_k \le 4^3$ に制限されている場合、利用可能な最も正確なアプローチはAB外挿モデルであるが、それでも完全な収束したTDLに対しては約0.25 eVの残留誤差を伴うと結論付けている。得られた結果は、EOM-CCSDと実験との間の残りの不一致（0.3–0.4 eV）が、有限サイズ誤差ではなく、電子相関の不十分な取り扱い（具体的には三重励起の必要性）または電子・フォノン相互作用に起因することを示唆している。これは、より安価な手法のベンチマークのための新しい基準を設定し、固体電子構造におけるCCSDの現在の限界を明確にするものである。