Reaching the thermodynamic limit of periodic CCSD cohesive energies and band gaps

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

問題提起
周期的な結合クラスター理論（CCSD）は、単一および二重励起を伴う非摂動的でサイズ一貫性のある手法であり、固体電子構造の高精度な予測を可能にする。しかし、その適用は計算コストによって厳しく制限されてきた。ストレージ要件は $O(n^4 N_k^3)$、実行時間は $O(n^6 N_k^4)$ でスケールする（ここで $n$ は単位胞あたりの原子数、$N_k$ は $k$ 点数）。その結果、これまでの研究では通常、$k$ 点メッシュは $4^3$ を超えない範囲に限定されていた。この制限は、基底状態エネルギー（$N_k^{-1}$）と比較して収束が遅いバンドギャップ（$N_k^{-1/3}$）における有限サイズ誤差の確実な除去を妨げている。十分な密度の $k$ 点サンプリングにアクセスできないことは、間接バンドギャップや、小さなメッシュでの予測が漸近収束領域から外れてしまう材料において、熱力学的極限（TDL）への外挿を信頼できないものにする。

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

• 並列化: コードは空間群対称性を活用し、分散メモリおよび共有メモリの並列化を利用しており、最大12ノード（各ノード96コア、1.5 TBメモリ）で効率的に動作する。
• サンプリング密度: この実装により、ダブルゼータ（DZ）基底では最大 $6^3 = 216$ 個、$3$ 重ゼータ（TZ）では $5^3$ 個、$4$ 重ゼータ（QZ）では $4^3$ 個の $k$ 点を用いた計算が可能となった。
• 外挿スキーム: 本研究では、TDL収束のための複数の外挿モデルを評価している：
  • モデル A: 単項式（$N_k^{-1/3}$）。
  • モデル AB: 二項式（$N_k^{-1/3} + N_k^{-2/3}$）。
  • モデル ABC: 三項式（$N_k^{-1/3} + N_k^{-2/-2/3} + N_k^{-1}$）。
  • GW-EOM スキーム: 先行文献（Ref. 19）で提案された手法であり、EOM-CCSDの展開係数と $G_0W_0@HF$ バンドギャップの間の比例関係を仮定することで収束を加速させる。
• 基底関数およびポテンシャル: 計算には、相関一貫基底（cc-pVXZ）およびハートリー・フォック理論に最適化されたガウス型軌道（GTH）擬ポテンシャルを使用している。全電子（AE）による凍結コア（FC）近似と、擬ポテンシャル（PP）計算の両方を実施した。
• 間接ギャップ戦略: 間接バンドギャップに対して、著者らは複合的なアプローチを提案している。間接ギャップを直接外挿するのではなく（これは $k$ メッシュのシフトに敏感であるため）、価電子帯最大（VBM）における直接ギャップと、VBMと伝導帯最小（CBM）の間のオフセットを個別に外挿し、それらを合算する。

主な貢献および結果

1. ソフトウェア実装: 著者らは、標準的な（非局在化）軌道に基づいた、現時点で最大の周期的なCCSD計算を報告しており、8つの単純な半導体および絶縁体（MgO, LiCl, LiF, LiH, BN, BP, Si, C）に対して信頼できるTDL外挿を実現した。
2. 凝集エネルギー:
   • 凝集エネルギーは0.1 eV以内で収束した。
   • CCSDの結果は、零点エネルギー補正後の実験値と比較して平均0.1–0.2 eVの誤差を示しており、通常、凝集エネルギーを過小評価する傾向がある。
   • 擬ポテンシャルの誤差は小さく（ほとんどの固体で <0.05 eV）、主にハートリー・フォックレベルで決定されることが判明した。
   • CCSDの精度は、一般的な密度汎関数（PBE, PBEsol, SCAN）と同等またはそれ以上である。
3. バンドギャップおよび外挿解析:
   • 収束性: バンドギャップは凝集エネルギーよりも収束が困難である。本研究は、$N_k=2^3$ では外挿モデルにおける高次項を無視するには不十分であり、$N_k=3^3$ が一般に十分であることを確立した。
   • モデルの性能: $3^3$ から $6^3$ の $k$ メッシュを用いたABC外挿モデルが、収束したリファレンスとして機能する。単項（モデルA）の外挿は、小さなメッシュに限定されている場合、バンドギャップを一貫して0.2–0.5 eV過小評価する。
   • GW-EOM スキーム: $2^3$–$4^3$ メッシュから係数をフィッティングするGW-EOM-234スキームは、不安定な結果をもたらす。これは一部の材料（MgO, LiCl, LiF）では良好に機能するが、GWとCCSDの間の比例関係の仮定が崩れる場合、他の材料（BN, BP, Si, C）に対してバンドギャップを著しく過大評価する。
   • 間接ギャップ: 複合的アプローチ（直接ギャップ + オフセット）は、間接ギャップの $k$ メッシュ・シフトに対する敏感さを効果的に緩和し、約0.1 eVの不確かさで収束した値を提供した。
4. ベンチマーク値: 本論文は、CCSDの凝集エネルギーおよびバンドギャップに関する決定的なベンチマーク値を提供している。
   • バンドギャップの精度: EOM-CCSDは、零点再正規化（ZPR）補正後の実験値に対して、平均絶対誤差（MAE）約0.4 eVでバンドギャップを予測する。理論は通常、ギャップを過大評価する。
   • 外れ値: 岩塩型構造（MgO, LiH, LiF, LiCl）において顕著な偏差（0.8–1.0 eV）が見られる。単一励起の性質（$n_{1}^{IP}$）の分析は、これらの材料が他の材料（>95%）と比較して低い単一励起の重み（93–94%）を持っていることを示唆しており、これらの特定の系においては三重励起（例：CCSD(T)）の導入が必要であることを示している。
5. TiO$_2$ への適用: 本手法はルチル型TiO$_2$に適用され、TDLバンドギャップとして4.17 eVを予測した。これはZPR補正後の実験値を約0.3 eV上回っており、他の固体で見られる一般的な傾向と一致しているが、著者らは遷移金属の存在と異方的なメッシュが追加の不確実性を導入していると指摘している。

意義
本論文は、単純な固体の凝集エネルギーおよびバンドギャップに対して、初めて信頼できる収束したCCSDベンチマークを提供し、有限サイズ誤差および基底関数不完全性の誤差を効果的に排除したと主張している。標準的な外挿スキーム（モデルAなど）がバンドギャップを系統的に過小評価すること、およびGW-EOMスキームが普遍的に堅牢ではないことを示すことで、本研究は多くの体理論の精度を評価するための新しい基準を確立した。結果は、CCSDは凝集エネルギーに対しては非常に高精度であるが、バンドギャップに対しては約0.4 eVの系統的な過大評価を保持しており、これはイオン性物質における電子相関（特に三重励起）の不十分な取り扱い、または電子・フォノン相互作用に起因する可能性が高いことを示唆している。開発されたソフトウェアとデータは、将来の高精度な固体計算の基礎となるものである。