Ferroelectric quantum critical point in superconducting hydrides: The case of H$_3$S

機械学習ポテンシャルを用いた経路積分分子動力学を用いることで、本研究は、H$_3$Sにおける高温超伝導が、約134 GPaにある4次元イジング普遍類に属する強誘電量子臨界点の上方において、大きな核量子ゆらぎに支配された常誘電領域で発生することを明らかにしている。

要約

超高圧下で、硫黄と水素の原子からなる極めて高密度な小さな球体があると想像してみてください。この球体は、極限の圧力下で超伝導体、つまり電気抵抗ゼロで電気を流す物質へと変化します。長い間、科学者たちはなぜ「H3S」と呼ばれる特定の物質でこのような現象が起きるのか、その理由に頭を悩ませてきました。彼らは、特定の圧力（約155 GPa）で最も効果的に機能することは分かっていましたが、原子がどのように振る舞うかという「地図」が欠けていたのです。

この論文は、その欠けていた地図を描くようなものです。研究者たちは、原子を単なる「固体の球」としてではなく、「確率のぼやけた雲（量子効果）」として追跡するために、強力なコンピュータ・シミュレーションを用いました。彼らが発見したことを、分かりやすく説明します。

1. 「ぼやけた」原子と魔法の圧力

微小な原子の世界では、物事は固定的ではなく、ゆらぎ、震えています。研究者たちは、特定の圧力（約134 GPa）において、H3S中の水素原子がある「転換点」に達することを発見しました。

• 比喩： ボウルの中に置かれたボールを想像してください。ボウルが深ければ、ボールは中央に留まります。もしボウルを揺らしたり（熱）、押しつぶしたり（圧力）すれば、ボールは動き回るかもしれません。
• 発見： この「量子臨界点（QCP）」と呼ばれる転換点において、原子は最大の混乱状態にあります。原子は一箇所に落ち着いているわけでもなく、かといって完全にランダムなわけでもありません。彼らは、どちらに曲がるべきか決めることができない群衆のように、激しく「揺らいで（fluctuating）」いるのです。

2. 相転移：「対称的」から「非対称」へ

この物質には、主に2つの形状（相）が存在します。

• 「完璧にバランスの取れた」相（常誘電相）： 水素原子が硫黄原子の間の中央に位置しています。これはシーソーが完璧に均衡している状態のように、対称的です。
• 「偏った」相（強誘電相）： 水素原子が片側に押しやられます。シーソーが傾いてしまうのです。

論文によれば、「バランスの取れた」状態から「偏った」状態への転移は、超伝導が最も強く現れる場所と正確には一致しません。むしろ、超伝導のピークは「バランスの取れた」領域にありますが、原子が最も激しく揺れ動いている転換点のすぐ隣に位置しています。

3. 「超伝導のスイートスポット」

ここに大きな驚きがあります。

• 旧来の考え： 科学者たちは、超伝導のピークは、物質が「バランスの取れた」状態から「偏った」状態へと切り替わることで起こると考えていました。
• 新しい発見： 論文は、ピークが実際には「バランスの取れた」領域、つまり混沌のすぐ隣で起きていることを示しています。
• 比喩： サーファーを想像してみてください。最高の波は、穏やかで平坦な水面でもなく、混沌とした荒波でもありません。最高の波は、海が荒れ始めているまさにその場所にあります。この「荒々しさ（量子ゆらぎ）」が、電子をペアにし、抵抗なく流れるのを助けるのです。論文は、転換点付近での水素原子の激しい揺らぎが、超伝導へのブーストとして機能していることを示唆しています。

4. 「4Dイジング」のルールブック

研究者たちは、この転換点の背後にある数学を分析し、それが「4Dイジング普遍性クラス」として知られる非常に特定のルールブックに従っていることを見出しました。

• 比喩： さまざまなゲーム（チェス、チェッカー、囲碁など）を想像してください。それらは見た目は違っても、駒の動き方に関する共通の基礎的な論理に従っています。研究者たちは、これらの原子の振る舞いが、4次元における状態変化を記述するために用いられる、特定の複雑な数学モデルと同じ「論理」に従っていることを発見しました。これにより、彼らの発見が単なる偶然ではなく、物理学の根本的な法則であることを裏付けました。

5. なぜこれが「地図」にとって重要なのか

この研究の前まで、H3Sの地図はぼやけていました。研究者たちは、従来のメソッドでは計算コストが高すぎるシミュレーションを実行するために、新しいタイプの「AIの脳（機械学習ポテンシャル）」を使用しました。

• 彼らは、もし量子的な「ぼやけ（fuzziness）」を無視して（原子を固いビリヤードの球として扱うと）、正しい地図が得られないことを突き止めました。量子的なゆらぎによって、転移圧力は約50 GPaも大きくシフトしてしまうのです。
• これらのゆらぎを含めることで、彼らはついに「転換点（QCP）」を特定し、超伝導のピークが、転換点の上方にある「強い量子ゆらぎ」の領域に位置していることを示しました。

まとめ

この論文は、H3Sにおける超伝導のマジックが、単に物質の形が変わることによって引き起こされるのではないことを明らかにしています。そうではなく、物質が原子が激しく振動する「量子の転換点」のすぐそばに漂っているために起こるのです。これらの激しい振動が触媒として作用し、電気を完璧に流れるように助けています。研究者たちは、これがどこで起きているかを正確にマッピングし、それが特定の、普遍的な数学的ルールに従っていることを証明しました。

技術概要

問題の背景
硫化水素（H$_3$S）は、約155 GPaにおいて203 Kという高い転移温度（$T_c$）を示す、著名な高温超伝導体である。その超伝導ドームは、銅酸化物やドープされたSrTiO$_3$と同様のピークを示すが、構造相転移と超伝導極大値との関係性を含め、相図の正確な性質については未解決のままである。先行研究では、水素の変位によって駆動される、体心立方構造の$Im\bar{3}m$相から三方晶の$R3m$相への強誘電相転移が示唆されていた。しかし、これらの理論では、この転移圧力が実験的な$T_c$のピークよりも大幅に低い（40 GPa以上）値に設定されていた。さらに、広範な温度および圧力範囲におけるH$_3$Sの包括的な相図の理解、特に熱的および核的量子効果（NQE）の相互作用に関する理解も不足していた。

手法
これらの課題に対処するため、著者らは核的量子効果を厳密に扱うためにパス積分分子動力学法（PIMD）を用いた。必要なシステムサイズと温度に対してPIMDを密度汎関数理論（DFT）などの第一原理電子構造計算と組み合わせることは、計算コストが極めて高いため、著者らは機械学習原子間ポテンシャル（MLIP）を利用した。具体的には、40から220 GPaの圧力範囲、100 Kおよび200 Kの条件におけるDFT-BLYP構成のデータセットを用いて、MACE（高次等変メッセージパッシングニューラルネットワーク）ポテンシャルを訓練した。

シミュレーションは、有限サイズスケーリング解析を可能にするために、$L \times L \times L$の基本単位格子（$L \in \{2, 3, 4\}$）を用いたi-PIパッケージによるNVTアンサンブルを用いて実施された。研究対象とした温度範囲は50 Kから300 Kである。主要な観測量には以下が含まれる：

• 秩序変数： 常誘電相と強誘電相を区別するための、全プロトン変位（$\Delta$）およびその絶対値（$\Delta_{abs}$）。
• 局所的観測量： 局所的なプロトン変位の分布および局所的な双極子モーメント。
• 動的特性： 非調和フォノン周波数を抽出し、遅延効果を研究するための虚時間フォノングリーン関数（$G_{ij}(\tau)$）。
• スケーリング解析： 臨界点付近における秩序変数の有限サイズスケーリングによる普遍性クラスの決定。

主な結果

1. 相図と量子臨界点（QCP）： 本研究は、温度が低下するにつれてより低い圧力へとシフトする強誘電転移線を示している。熱力学的極限および$T=0$ Kへの外挿により、著者らは$p_{QCP} \approx 134 \pm 2$ GPaにおける強誘電量子臨界点（QCP）を特定した。
2. 普遍性クラス： QCP付近の秩序変数の有限サイズスケーリングにより、臨界指数 $\nu \approx 0.483 \pm 0.014$ が得られた。この値は、ローレンツ不変性（$z=1$）を仮定した場合、4次元イジング普遍性クラス（$\nu = 1/2$）と一致する。
3. 実験的な$T_c$との乖離： 計算された強誘電転移線（例：200 Kにおいて124 GPa）は、実験的な超伝導ピークである155 GPaよりも著しく低い。したがって、実験的な$T_c$の極大値は、強誘電相ではなく、常誘電の$Im\bar{3}m$相内に位置している。
4. 量子ゆらぎの役割： QCPの周辺および常誘電相へと広がる領域は、大きな核的量子ゆらぎによって特徴付けられる。著者らはこれを、虚時間 $\tau = \beta/2$ における局所グリーン関数の分散と $\tau = 0$ における分散の比（$\sigma^2_g(\beta/2)/\sigma^2_g(0)$）を用いて定量化した。この比は、QCP付近の常誘電相でピークに達し、強い遅延効果と双極子モーメントのゆらぎを示している。これらは、プロトンが凍結する強誘電相では抑制される。
5. 古典的動力学との比較： 古典的分子動力学（MD）シミュレーションでは、強誘電転移は量子的な場合と比較して（200 Kにおいて約50 GPaシフトして）はるかに高い圧力で発生しており、NQEが常誘電相を安定化させ、転移圧力を低下させる上で決定的な役割を果たしていることが浮き彫りになった。
6. フォノン軟化： 転移を駆動する光学フォノンモードは、常誘電側から転移に近づくにつれて軟化を示し、その後、強誘電相に入ると急激な周波数ジャンプを示す。これは二次の相転移と一致している。

意義および主張
本論文は、H$_3$Sにおける高温超伝導が高強誘電相転移そのものではなく、QCP付近の強い量子ゆらぎと遅延効果に支配された常誘電領域で発生していることを確立した。著者らは、これらのゆらぎが、おそらく強誘電QCPへの近接性によって媒介されており、超伝導を強化している可能性を示唆している。このシナリオは、完全な周波数依存性や頂点補正をしばしば無視する標準的なミダル・エリアジアス理論が、H$_3$Sを記述するには不十分である可能性を意味している。強誘電QCPと強い光学フォノン軟化の存在は、希薄金属であるSrTiO$_3$に対して提案されているような、線形次を超えた電子対形成メカニズムの可能性を支持している。本研究は、核的量子効果と非摂動的な電子・フォノン結合評価の必要性を強調しつつ、超伝導性水和物の基礎となる構造的不安定性と超伝導ドームの位置関係を整合させた、詳細かつ量子力学的に正確な相図を提供している。