Accurate and efficient protocols for high-throughput first-principles materials simulations

技術要約：高スループット第一原理材料シミュレーションのための正確かつ効率的なプロトコル

問題提起
高スループットな材料探索の進展は、自動化された第一原理シミュレーション、特に密度汎関数理論（DFT）に依存している。これらのワークフローを数万種類の化合物へとスケールアップさせる際の決定的なボトルネックは、数値的な精度と計算効率のバランスを取るための、信頼できるシミュレーションパラメータの選択である。具体的には、k点サンプリング（ブリルアンゾーンの離散化）とスメアリング手法（特に金属におけるフェルミ準位での電子占有の不連続性を処理するために使用される）の相互作用が複雑である。

• サンプリング誤差： 金属系においては、不適切なk点サンプリングは「統計的」誤差を引き起こし、メッシュ密度に応じて結果がランダムに変動するため、収束には極めて高密度なメッシュが必要となる。
• 系統的誤差： スメアリングは、電子の占有関数を滑らかにするために仮想的な電子温度（$\sigma$）を導入し、k点の収束を加速させる。しかし、大きな$\sigma$値は、エントロピーの寄与によって全自由エネルギーが物理的な $\sigma \to 0$ の極限から逸脱するため、「系統的」誤差を導入する。
  現在のワークフローでは、固定されたパラメータを使用するか、スメアリングの最適化を行わないことが多く、サンプリング誤差を抑制することと系統的誤差を最小化することの間のトレードオフを制御できていない。

手法
著者らは、最適化されたk点およびスメアリングパラメータを選択するための、厳密なベンチマークフレームワークである**標準固体状態プロトコル（SSSP）**を開発した。

1. データセット： Materials Cloudデータベースからサンプリングされた269個の結晶構造（3Dおよび2D）を使用した。対称性に起因する力の消失を避けるため、各ユニットセル内の原子を結合長の5%分だけランダムに変位させた。
2. ワークフロー： Quantum ESPRESSO (v7.0) と SSSP PBEsol Efficiency 1.2 擬ポテンシャルライブラリを用いて、自動化されたDFT計算を実施した。
   • 変数： k点サンプリング密度は、逆格子空間における隣接するk点間の最大距離（$\mathring{A}^{-1}$単位）である $\lambda$ によってパラメータ化された。スメアリング温度（$\sigma$）は、Marzari-Vanderbilt cold smearing を用い、0.002から0.03 Ryの範囲とした。
   • 指標： 全エネルギー（$E$）、変位した原子にかかる力（$F$）、および擬似応力（$\Sigma$）を算出した。
3. 誤差解析： 著者らは誤差を2つの異なるカテゴリーに分離した。
   • サンプリング誤差 ($\varepsilon^k$)： 特定の $(\lambda, \sigma)$ で計算された特性の、同じ $\sigma$ における最も密なk点メッシュ ($\lambda_{ref}$) で得られた値からの偏差。
   • 系統的誤差 ($\varepsilon^s$)： 特定の $\sigma$（$\lambda_{ref}$ において）で計算された特性の、最小の $\sigma$ ($\sigma_{ref}$) における参照値からの偏差。
4. 最適化戦略： 著者らは、与えられた平均誤差に対して計算の壁時間（wall time）を最小化する「パレートフロント」のパラメータを探索した。決定的なのは、サンプリング誤差が系統的誤差の5%未満になるよう制約を課したこと（$\varepsilon^k < 0.05\varepsilon^s$）である。これにより、制御が困難なサンプリング誤差が全誤差予算を支配しないようにし、精度とコストの調整のための主要な変数として系統的誤差（$\sigma$ によって制御されるもの）を扱えるようにした。

主な貢献

1. 分離された誤差解析： 本論文は、統計的（サンプリング）誤差と系統的（スメアリング）誤差を個別に定量化および制御する手法を確立しており、これらが偶然に打ち消し合って収束の問題を隠蔽してしまうことを防いでいる。
2. SSSPプロトコルの定義： パレート最適化に基づき、精度と効率のトレードオフに合わせて調整された3つの異なるプロトコルを定義した。
   • Fast（高速）： テストや予備的な緩和に最適化されている。粗いメッシュと高いスメアリングを使用する。
   • Balanced（バランス型）： 高スループット・スクリーニングおよび一般的な用途に推奨される。
   • Stringent（厳格）： 計算コストが急激に上昇する前の、最高精度の選択肢。
3. 材料特異的なガイドライン：
   • 絶縁体： サンプリング誤差の影響を受けにくいため、Balanced プロトコルの恩恵を大きく受け、より粗いメッシュを利用できる。
   • f電子元素（ランタノイド/アクチノイド）を含む金属： フェルミ準位付近の電子状態密度の鋭いピークにより、高い系統的誤差を示す。精度を維持するためには、Stringent プロトコル（より低い $\sigma$ と密なメッシュ）が必要となる。
4. ソフトウェア統合： これらのプロトコルは AiiDA ワークフローフレームワークおよび Materials Cloud 上の Quantum ESPRESSO 入力生成器に統合されており、ユーザーは材料のクラスと目的のプロトコルに基づいて入力を自動生成できる。

結果

• 収束挙動： スメアリング温度を上げるとサンプリング誤差が効果的に抑制され、より粗いk点メッシュの使用が可能になることを本研究は確認している。しかし、サンプリング誤差と系統的誤差の和が、特定の計算コストに対して最小となる「スイートスポット」が存在する。
• ランタノイドの感度： ランタノイドを含む材料は、標準的なスメアリング温度において、他の金属と比較して著しく高い系統的誤差を示す。これらのシステムにおいて許容可能な誤差閾値を維持するには、Stringent プロトコルが必要である。
• プロトコルの性能：
  • Balanced プロトコルは、金属に対して約20.9 meV/Å の平均力誤差をもたらし、Stringent プロトコルよりも2.7倍高速である（Stringent は Balanced より2.7倍遅い）。
  • Stringent プロトコルは、金属に対して平均力誤差を約11.5 meV/Å まで低減させる。
  • 絶縁体の場合、Balanced プロトコルで十分であり、Stringent プロトコルを実行することは、誤差のわずかな改善に対して1.82倍高い計算コストを強いることになる。
• 転移性： 派生した特性（格子定数およびフォノン周波数）に関するベンチマークは、これらの「高次」特性における誤差が力（force）の誤差と直接相関していることを示しており、力を主要な収束指標として用いることの妥当性を検証している。

意義および主張
著者らは、SSSPプロトコルが、DFTパラメータの設定における堅牢で自動化されたソリューションを提供し、ケースバイケースの収束テストの必要性を排除すると主張している。これらのプロトコルによって、複雑なk点とスメアリングの選択を抽象化することで、以下のことが可能になる：

• 再現可能かつ一貫した高スループット・ワークフローの実現。
• 過剰に収束した計算への計算リソースの浪費、あるいは収束不足のデータによる信頼性の低いデータの生成の防止。
• 原理的にはコードに依存しない設計（ただし、Quantum ESPRESSO向けに実装されている）であり、その数値は他の平面波コード（VASP, Abinitなど）や、スメアリングを利用する全電子コードに対する第一近似として機能する。
• 材料設計のためのモジュール式構成要素としてワークフローを扱える、相互運用可能なソフトウェアエコシステムへの基礎的なステップである。

本論文は、これらのプロトコルは全エネルギー、力、および応力の自己整合的DFT計算に対して最適化されているものの、特定の線形応答特性（qグリッドを用いたフォノン計算など）については、追加の収束テストが必要になる可能性があることを強調している。