[論文レビュー] Hydrogen phase-IV characterization by full account of quantum anharmonicity

Tommaso Morresi, Rodolphe Vuilleumier+1 more|arXiv (Cornell University)|Mar 12, 2022

Quantum, superfluid, helium dynamics参考文献 50被引用数 8

ひとこと要約

本稿では、強い非調和性を示す分子結晶におけるフォノン周波数を、量子核効果を完全に考慮して高精度に計算するための、制約付き変位演算子を備えた新しい経路積分分子動力学（PIMD）フレームワークを提示する。マツバラグリーン関数の静的極限を一般化固有値問題を用いて計算することで、誤差が最小限に抑えられ、不適切な解析的継続を回避する。この手法により、固体水素の相IIIおよび相IVにおける振動モードを成功裏に予測し、実験的ラマンおよび赤外分光法データ（同位体シフトを含む）との直接比較によって、相IVの結晶対称性がC2/c-24であると特定した。

ABSTRACT

We devise a framework to compute accurate phonons in molecular crystals even in case of strong quantum anharmonicity. Our approach is based on the calculation of the static limit of the phononic Matsubara Green's function from path integral molecular dynamics simulations. Our method enjoys a remarkably low variance, which allows one to compute accurate phonon frequencies after a few picoseconds of nuclear dynamics, and it is further stabilized by the use of appropriate constrained displacement operators. We applied it to solid hydrogen at high pressure. For phase III, our predicted infrared (IR) and Raman active vibrons agree very well with experiments. We then characterize the crystalline symmetry of phase IV by direct comparison with vibrational data and we determine the character of its Raman and IR vibron peaks.

研究の動機と目的

• 大きな量子核揺らぎが生じる強非調和的分子結晶の振動的性質を予測するための調和的および摂動的手法の限界を克服すること。
• 固体水素のような分子固体において、核の量子効果（NQE）と非調和性を完全に考慮する第一原理的手法を開発すること。
• 長年の曖昧さが残る水素相IVの結晶対称性を、振動分光法を用いて特定すること。
• 制約付き一般化固有値定式化により、PIMDに基づくフォノン計算におけるサンプリング誤差と分散を低減すること。

提案手法

• 本手法は、第一原理力に基づく経路積分分子動力学（PIMD）シミュレーションを用いて、フォノンのマツバラグリーン関数の静的極限を計算する。
• 解析的継続を回避するため、逆マツバラグリーン関数へのアクセスにKubo変換された相関関数を用いる。
• 混合モードの干渉を抑えるために、制約付き変位演算子を用いた一般化固有値問題（GEV）を解く。
• 制約付き変位演算子により、 librational モードからのバイアスが低減され、分子結晶内の関連するバイブロンモードが分離される。
• 正規モードの最大局在化基準を用いることで、PIMD軌道における収束性が向上し、統計的分散が低減される。
• フォノン周波数は、逆Kubo変換グリーン関数の固有値から直接抽出され、強い非調和性下でも高精度が保証される。

実験結果

リサーチクエスチョン

• RQ1 理論的予測と実験的振動分光データに矛盾する中、水素相IVの正しい結晶対称性は何か？
• RQ2 強い量子非調和性と大きな核の量子揺らぎを示す分子結晶において、フォノン周波数をどのように高精度に計算できるか？
• RQ3 PIMDに基づく手法が、解析的継続に起因する数値的不安定性と高い分散を回避しながら、NQE補正されたフォノン周波数を捉えられるか？
• RQ4 分子固体における標準的PIMDフォノン計算において、分子内librationalモードがどの程度バイアスを引き起こすか？
• RQ5 ラマンおよび赤外ピークの同位体シフトは、相IVにおける競合する結晶構造を区別するためにどの程度有効か？

主な発見

• 制約付きPIMD（cPIMD）フレームワークは、数ピコ秒のシミュレーション時間で高精度なフォノン周波数を達成し、サンプリング分散が顕著に低減された。
• 固体水素相IIIでは、予測された赤外およびラマン活性バイブロン周波数が、実験測定値と定量的に一致した。
• 本手法により、実験的振動分光データとの直接比較に基づき、C2/c-24構造が相IVの正しい結晶対称性であると特定された。
• 相IVにおけるラマンおよび赤外バイブロンピークの振るまいと同位体シフトが、正しく予測され、実験観測と一致した。
• Cmca-12およびPbcn-48構造は、280 GPaおよび20 KにおけるPIMDシミュレーションで不安定であり、混合層構造への相転移が観察された。
• 調和近似は、相IVにおけるバイブロンモードを最大で1000 cm⁻¹低く見積もっていたが、cPIMD手法は強い非調和的補正とΓ点における軟化を捉えられた。