Capturing Nuclear Quantum Effects in Hydrogen Diffusion through MoS2 via Machine-Learning-Enhanced Path-Integral Simulations

技術要約：機械学習強化パス積分シミュレーションによるMoS₂を通じた水素拡散における核量子効果の捕捉

問題提起
層状二次元（2D）材料（二硫化モリブデン（MoS₂）など）を通じた水素輸送は、水素貯蔵、燃料電池、および同位体分離などの用途において極めて重要である。実験的研究では、MoS₂やh-BNのような材料における層間拡散や同位体選択的輸送が示されているが、正確な理論的記述は依然として困難である。従来の理論的研究の多くは原子核の古典的な扱い方に依存しており、核量子効果（NQE）を無視してきた。しかし、水素のような軽い原子は、室温であっても零点振動（ZPE）、トンネル効果、および非局在化といった顕著な量子挙動を示す。さらに、これらの系を正確にモデリングするには、van der Waals（vdW）相互作用と水素結合の両方を捉える必要があり、密度汎関数理論（DFT）における高度な交換相関汎関数（例：r2SCAN+rVV10）が必要となる。このような高精度DFTと、パス積分分子動力学（PIMD）および拡張サンプリングを、ねじれ二層構造のような構造的に複雑で大規模な系に対して組み合わせる計算コストは、極めて高い。

手法
これらの計算上の課題に対処するため、著者らはr2SCAN+rVV10 DFTを用いて生成されたコンパクトなデータセットに基づき、MACEフレームワークに基づく機械学習原子間ポテンシャル（MLIP）を開発した。このMLIPは、エネルギー安定性とフォノンバンド構造に関してDFT参照計算に対して検証され、良好な一致を示した。

本研究では、以下のハイブリッドシミュレーションフレームワークを用いた：

1. ウェル・テンパード・メタダイナミクス（WTMetaD）： PLUMEDを介して実装され、層間拡散の自由エネルギー障壁（$\Delta F$）を算出する。
2. パス積分分子動力学（PIMD）： 核量子効果（NQE）を明示的に組み込むため、i-PIパッケージを介して実装された。核は、300 Kにおいて、水素（および同等の収束を保証するための重水素）に対して$P=16$個のビードを持つリングポリマーとして表現された。
3. 系構成： 高対称バルクMoS₂ポリタイプ（$H_h^h$, $H_X^h$, $R_M^h$, および3R相）および、ねじれ角 $\theta = 3.89^\circ, 21.79^\circ, 38.21^\circ, 56.11^\circ$ のねじれ二層MoS₂（TBL-MoS₂）を調査した。
4. 同位体分析： 運動学的同位体効果（KIE）を定量化するために、プロチウム（H）と重水素（D）の両方の拡散をシミュレートした。

主要な貢献

• MLIP-PIMDフレームワークの開発： 著者らは、複雑な2D系における量子核動力学に必要な長さ/時間スケールと、DFTの精度との間の溝を埋める堅牢なワークフローを確立した。
• MoS₂におけるNQEの定量化： 本研究は、量子ゆらぎがどのようにMoS₂の異なるポリタイプにおける拡散景観を変化させるかを系統的に評価している。
• 同位体選択的メカニズム： 本研究は、同位体分離を駆動するメカニズムについて原子論的な洞察を提供し、特にHとDの間の自由エネルギー障壁の差を定量化している。
• モアレ駆動の拡散景観： ねじれ二層構造における水素輸送をマッピングし、局所的なスタッキング（スタッガード対エクリプス）およびソリトンネットワークによって支配される、空間的に変化する拡散環境を明らかにした。

主要な結果

• 拡散障壁へのNQEの影響： NQEの導入により、古典的なシミュレーションと比較して、水素拡散の自由エネルギー障壁が大幅に低下した。減少幅は、スタッキング構成に応じて数十meVからNearly 60 meVに及ぶ。
  • 最も安定な$H_h^h$相では、障壁は122 meV（古典）から64 meV（量子）へと低下した。
  • $H_X^h$相では、障壁は45 meVから18 meVへと減少した。
  • この自由エネルギー景観の軟化は、自己拡散係数（$D$）の著しい増加をもたらし、量子補正因子（QCF）は$H_h^h$相において最大9.43に達した。
• 運動学的同位体効果（KIE）： 顕著なKIEが観察された。古典的なシミュレーションでは重水素の方が水素よりも高い障壁を示すが、NQEの導入はこの差を増幅させる。プロチウムの量子自由エネルギー障壁は、重水素（$H_h^h$において95 meV）よりも約35 meV低くなっており（60 meV）、強い質量依存的な量子拡散を示している。
• ねじれ二層ダイナミクス：
  • 高角度ねじれ（$21.79^\circ, 38.21^\circ$）： これらの系は拡張された高対称ドメインを欠いており、その結果、拡散景観は中間的な構成によって支配される。したがって、拡散障壁は高く（量子状態で約120 meV）、拡散係数は（孤立した高対称相よりも少なくとも一桁）低くなる。
  • 低角度ねじれ（$3.89^\circ, 56.11^\circ$）： これらの系は、明確なモアレドメイン（スタッガード）とノード/ソリトン（エクリプス）を形成する。水素は、層間距離が小さく配位が有利なスタッガードドメインに優先的に存在する。
  • 空間的変動： 拡散障壁はモアレ超格子全体で大きく変化し、有利なスタッガードドメインと不利なソリトン/エクリプス領域の間で最大80 meVの差が生じる。$56.11^\circ$の系は、メタステーブルなノードと有利な$H_X^h$ドメインの存在により、$3.89^\circ$の系（約197 meV）と比較して低い平均障壁（約150 meV）を示す。

意義
本論文は、周囲温度においても、層状MoS₂を通じた水素輸送においてNQEが決定的な役割を果たすと結論付けている。これらの効果を無視することは、拡散障壁の過大評価および拡散係数の過小評価につながる。本研究は、量子的な増強の大きさが高精度なスタッキング環境に強く依存していることを示している。さらに、同位体選択的な輸送メカニズムを正確に記述するためには、量子核動力学の明示的な導入が不可欠である。最後に、本研究は、ねじれエンジニアリングを用いることで拡散景観を調整し、水素移動のための空間的に選択的なチャネルを作成できることを示している。提示されたMLIP-PIMDフレームワークは、エネルギー技術に関連する複雑なファンデルワールスヘテロ構造における量子輸送を調査するための、転用可能なアプローチとして位置付けられている。