-Quantum Effects or Theoretical Artifacts? A Computational Reanalysis of Hydrogen's High-Pressure Phase Stability and Properties

問題提起
高圧下における水素の理論的特性評価は、その元素特有の量子的な性質と標準的な密度汎関数理論（DFT）の限界により、依然として大きな課題となっている。低温における水素の金属化に関する実験的確認は、350から440 GPaの間で半金属または金属状態への転移を示唆する相反するデータが存在することから、困難を極めている。理論的には、標準的な一般勾配近似（GGA）汎関数、特にPBEは自己相互作用誤差に苦しみ、人工的な電子の非局在化を引き起こす。これは、水素の結合エネルギー、バンドギャップ、および分子相と原子相の相対的安定性の記述において不正確な結果をもたらす。先行研究では、分子相における動力学的不安定性や非調和の兆候を量子核効果（QNE）に帰しているが、これらが交換相関汎関数のアーティファクト（偽像）であるかどうかの範囲は明らかになっていない。

手法
著者らは、400から700 GPaの範囲における水素の低温（0 K）相図について、以下の3つの異なる理論レベルを用いて比較計算再解析を行った：

1. GGA-PBE: 一般的に用いられる汎関数。
2. Meta-GGA (R2SCANおよびSCAN0): 相関と勾配展開をより良く記述するために、運動エネルギー密度を取り入れた汎関数。
3. Meta-Hybrid-GGA (SCAN0): 完全な交換項を含む高次レベルの汎関数。周期系における完全交換を用いた構造最適化の膨大なコストを軽減するため、R2SCANで最適化された幾何構造に対して一点エネルギー計算として適用した。

本研究では、$Cmca$-12、$Cmca$-4、および$C2/c$（Phase III）構造を含む候補相の相対エンタルピーを評価した。さらに、調和フォノンスペクトルによる格子力学の解析を行い、量子原子論（QTAIM）および周期エネルギー分解解析（PEDA）を用いて結合特性を調査した。

主な結果

• 相の安定性と金属化:

  • PBE: 原子状の$I4_1/amd$相が約507 GPaで基底状態になると予測しており、これは拡散モンテカルロ（DMC）法による予測よりも大幅に早い。すべての相が金属的であると予測されている。
  • R2SCAN (Meta-GGA): 原子相よりも高い圧力（~675 GPa）まで分子相（$Cmca$-4）を安定化させる。$C2/c$相は基底状態とほぼ縮退しているが、400 GPaにおいて0.41 eVのバンドギャップを有する。
  • SCAN0 (Meta-Hybrid): R2SCANの幾何構造に適用した際、SCAN0は400 GPaにおいて$C2/c$相を実質的なバンドギャップ（1.75 eV）を持つ基底状態として特定した。このギャップは、$Cmca$-12への転移過程である約470 GPaで閉じる。この絶縁体から金属への転移は、425–440 GPa以上での金属化を示す実験的な赤外およびラマンデータとよく一致している。原子状の$I4_1/amd$相の安定化は、850 GPa以上で起こると外挿される。

• 格子力学と非調和性:

  • PBE計算は、700 GPaにおける$Cmca$-4相の動力学的不安定性（虚のフォノンモード）を予測しており、非調和補正の必要性を示唆している。
  • 対照的に、R2SCAN計算では700 GPaまで$Cmca$-4に虚のフォノンモードは示されない。GGAによって予測された動力学的不安定性とソフトフォノンモードはメタGGAレベルで消失しており、これらの特徴が量子核効果ではなく、汎関数のアーティファクトである可能性を示唆している。
  • R2SCANは、QNEによって引き起こされるものと同様の分子相の体積膨張を予測しているが、これは結合の伸長ではなく、分子間のパウリ排斥の増大に起因するとしている。

• 結合解析:

  • 分子内結合: PBEは分子内のH–H結合を人工的に弱め、R2SCANの結果と比較して約0.1 Å長くしている。PBEは電子密度を過度に非局在化させ、HOMO-LUMOギャップを減少させ、結晶内における分子の共有結合性を弱めている。
  • 分子間相互作用: PBEは電荷の非局在化を通じて分子間相互作用を強化するが、R2SCANは、より強い分子内結合とより弱い分子間結合を伴う、より局在化した分子特性を保持している。
  • QTAIMおよびPEDA: 解析により、PBEは結合臨界点においてより低い電子密度（$\rho$）およびラプラシアン（$\nabla^2\rho$）値を与え、拡散した密度を示すことが確認された。PEDAは、R2SCAN構造がより高いパウリ排斥と、分子間のより低い軌道相互作用を示すことを明らかにしており、これが系の体積膨張を駆動している。

意義と主張
本論文は、水素の高温相図に対する交換相関汎関数の選択への敏感さが、しばしば見落とされる重要な要因であることを主張している。著者らは以下の通り述べている：

1. 汎関数の欠陥 vs 量子効果: 分子状水素相において、これまで量子核効果に帰せられてきた多くの動力学的不安定性や非調和の兆候は、実際には標準的なGGA汎関数の欠陥に由来している可能性がある。
2. 実験との一致の改善: メタGGAおよびメタハイブリッド汎関数を用いることで、予測される金属化圧力および相転移の順序（$C2/c \to Cmca$-12）は、PBEベースのモデルよりも実験的観察と著しく良く一致する。
3. 手法としての必要性: ポテンシャルエネルギー曲面の正確な記述、具体的には分子内相互作用と分子間相互作用のバランスの記述は、高密度水素の相安定性と結合強度を信頼性高く評価するための前提条件である。本研究は、このような汎関数の敏感性が、超イオン性や超伝導水素ケージのような、水素を含む他の量子状態にも同様に存在することを示唆している。