Anharmonicity and Nonadiabaticity in Hydride Superconductors

技術要約：水素化物超伝導体における非調和性と非断熱性

問題提起
高圧下における水素に富む化合物は、H$_3$SやYH$_6$のように転移温度（$T_c$）が200 Kを超え、LaH$_{10}$のように室温に近い値を示す、高温超伝導体の最も顕著なクラスである。これらの材料の発見と合理化はミダル・エリアシャス（ME）理論によって導かれてきたが、計算された$T_c$値と実験値との定量的一致は一様ではない。いくつかのベンチマークにおいて、最先端の調和近似および断熱近似に基づくME計算は、$T_c$を過大評価している。

従来のME理論における2つの主要な限界が、これらの不一致の潜在的な原因として特定されている：

1. 調和近似： 従来の理論は、軽い水素原子が大きな零点運動を行う際に不可欠となる、量子的な非調和格子動力学を無視している。この運動は、フォノンスペクトルと電子－フォノン（e-ph）結合を強く繰り込む。
2. 断熱近似（ミダルの定理）： 従来の理論は、フォノンエネルギーが電子のエネルギースケールに対して十分に小さい場合にのみ正当化される、e-ph頂点補正を無視している。水素に由来するフォノンは100–250 meVに達することがあり、フェルミ準位（$\varepsilon_F$）付近の電子構造の変化のエネルギースケール、特にvan Hove特異点が存在する場合において、これと同程度になる。

本研究が取り組む中心的な問いは、特定の水素化物において、どの近似（調和 vs 非調和、断熱 vs 非断熱）がME予測の精度を制限しているのか、そして両方の補正が同時に必要であるかどうかを判断することである。

手法
著者らは、いくつかの高度な計算手法を組み合わせた第一原理フレームワークを採用している：

• 非調和フォノン： 機械学習原子間ポテンシャル（ML/SSCHA）によって加速された確率的自己整合調和近似（SSCHA）を用いて計算される。この手法は、量子核効果を考慮し、フォノンの周波数と固有ベクトルを繰り込む。
• 電子－フォノン結合： 行列要素は、EPWコードを用いてワニエ表現で補間される。
• 超伝導特性： 等方的なエリアシャス方程式は、フェルミ準位での近似ではなく、電子状態密度（DOS）のエネルギー依存性を保持するフルバンド幅（FBW）理論を用いて解かれる。
• 頂点補正： 非断熱効果を考慮するために、第一次のe-ph頂点補正が含まれ、ペアリングカーネルを修正する。
• 対象物質： 本フレームワークは、異なる結合環境を持つ代表的な水素化物のセットに適用される：H$_3$S（共有結合的なS–Hネットワーク）、YH$_6$およびYH$_9$（クラスレート状のHネットワーク）、LaBeH$_8$（BeH$_8$ユニットを持つ三元系水素化物）、および常圧のPdH/PdD/PdTシリーズ。

主要な寄与と記述子
本論文は、標準的なME理論を超える補正の必要性に基づいて水素化物を分類するための、2つの材料固有の無次元記述子を導入している：

1. 非調和繰り込み（$A_\lambda$）： $|\lambda_{anh} - \lambda_{har}| / \lambda_{har}$ と定義され、非調和格子動力学によって誘起されるe-ph結合の相対的な変化を定量化する。
2. 頂点比（$R_V$）： $|\lambda_V^{anh}| / \lambda_{anh}$ と定義され、非調和結合に対する一次の頂点寄与の大きさを定量化する。

これらの記述子により、調和近似および断熱近似の破綻に関する異なる領域へと材料を分類することが可能になる。

結果
本研究は、調査した水素化物の中に3つの明確な領域があることを明らかにしている：

1. 強い非調和性と非断熱性（H$_3$S, YH$_6$, YH$_9$）：

   • これらの高圧水素化物は、非調和および非断熱の両方の補正を必要とする。
   • 非調和性は水素由来の光学枝を硬化させ、e-ph結合（$\lambda$）を25–40%抑制し、$T_c$を低下させる。
   • 頂点補正はさらに$T_c$を減少させる。例えば、200 GPaにおけるH$_3$Sでは、調和近似による$T_c$（244 K）は、非調和性によって196 Kに減少し、さらに頂点補正によって184 Kとなり、実験値の184 Kと一致する。
   • YH$_6$およびYH$_9$では、これら両方の補正を組み合わせることで、計算された$T_c$値（それぞれ226 Kおよび231 K）が、実験による測定値（224 Kおよび239 K）と密接に一致する。

2. 強い非調和性だが断熱的（PdH, PdD, PdT）：

   • 常圧のPd水素化物は強い非調和性を示しており、これは調和的な不安定性を除去し、光学的なラトリングモードを硬化させる。
   • しかし、これらは低いフォノンエネルギーと$\varepsilon_F$付近の滑らかなDOSを持つため、断熱的であり、頂点補正は無視できる（$R_V \approx 0.08$）。
   • 非調和性の取り扱いは、逆同位体効果（質量が増加すると$T_c$が増加する現象：PdH < PdD < PdT）を再現し、$T_c$の値として8.0、10.5、11 Kをそれぞれ得ており、実験と極めて良く一致している。

3. 弱い非調和性と有限の頂点補正（LaBeH$_8$）：

   • LaBeH$_8$は弱い非調和性（$A_\lambda \approx 0.02$）を示し、$\lambda$のわずかな減少のみが見られる。
   • 頂点補正は有限であり（$R_V \approx 0.11$）、$T_c$を約4 K（163 Kから159 Kへ）減少させる。
   • これらの補正にもかかわらず、予測される$T_c$は実験値である約110 K（80 GPaにおいて）よりも大幅に高いままである。著者らは、非線形e-ph補正や異方的な取り扱いが明らかな解決策にはならない（それらはむしろ$T_c$を増加させる可能性があるため）と指摘している。

意義
本論文は、水素化物超伝導体の$T_c$を定量的に予測するために必要な補正の階層を特定するための、第一原理フレームワークを確立している。記述子$A_\lambda$と$R_V$を導入することで、著者らは調和近似と断熱近似がいつ破綻するかを判断するための実用的な診断ツールを提供している。本フレームワークは、H$_3$S、YH$_6$、YH$_9$、およびPd水素化物の実験的傾向を正常に再現しており、鋭い電子的特徴を持つ高圧水素化物においては、非調和性と非断熱的頂点補正の同時的な包含が極めて重要であることを示している。LaBeH$_8$における持続的な不一致は、特定の材料においては、これら高度な補正ですら不十分である可能性を示唆しており、高次の効果や代替メカニズムへのさらなる調査の必要性を示している。