First-Principles Evidence for Strongly Correlated Superconductivity Driven by Structural Variations in La$_3$Ni$_2$O$_7$

第一原理計算により、La$_3$Ni$_2$O$_7$における高温超伝導が、圧力誘起の構造変化に伴う軌道局在とスクリーニング効果の競合によって増大する強い電子相関に駆動されていることを明らかにし、Ac$_3$Ni$_2$O$_7$との比較を通じて A サイト陽イオンの重要性を確証しました。

要約

🏗️ 1. 舞台設定：「二階建てのビル」と「変形する部屋」

まず、この物質の構造を想像してください。
La3Ni2O7 は、ニッケル（Ni）と酸素（O）が組み合わさった**「二階建てのビル」**のような構造をしています。

• 通常の状態（低圧力）： ビルの柱（ニッケルと酸素の結合）が少し傾いていて、部屋（結晶格子）が長方形（ひし形に近い）になっています。これは**「歪んだ状態」**です。
• 高圧力をかけると： 外部から強い圧力をかけると、ビルは圧縮され、柱がまっすぐに立ち、部屋は**「正方形」に整います。これを「構造相転移」**と呼びます。

実験では、この「まっすぐになった正方形の状態」で超伝導が起きることが分かっています。しかし、なぜ超伝導が起きるのか、その「スイッチ」がどこにあるのかは長年謎でした。

🔋 2. 核心の発見：「電子の喧嘩」と「壁の厚さ」

この研究の最大の見どころは、**「電子同士の関係性（相関）」**が、圧力によってどう変わるかを突き止めた点です。

電子は、同じ部屋（原子）に集まると、互いに反発し合います（これを「クーロン斥力」と言います）。

• 反発が弱いと： 電子は自由に動き回り、ただの金属になります。
• 反発が強いと： 電子は「喧嘩」を始め、奇妙な動き（超伝導など）をします。

研究者たちは、この物質をシミュレーションして、圧力をかけながら「電子の喧嘩の強さ」を測りました。すると、面白いことが分かりました。

📈 「右三角形」の謎を解く鍵

実験では、圧力と温度のグラフ上に、「右三角形」の形をした超伝導領域が現れます。

• 左端（圧力 14 GPa 付近）： 突然超伝導が始まる。
• 頂点（圧力 18 GPa 付近）： 最も高い温度（80 K）で超伝導になる。
• 右端（圧力 50 GPa 超）： 再び超伝導が弱まり、消える。

この研究は、この三角形の形が、**「電子の喧嘩の強さ（U/t 比）」**の変化と完全に一致することを発見しました。

🎭 3. 2 つの力が戦う「綱引き」

なぜ、圧力をかけると「喧嘩の強さ」が一度強くなり、その後弱まるのでしょうか？ ここには2 つの相反する力が働いています。

1. 力 A：「部屋を狭くする力」

   • 圧力をかけると、原子が押し合いへし合いになり、電子が閉じ込められます。
   • 結果： 電子同士の距離が縮まり、「喧嘩（反発）」が激しくなります。
   • これにより、超伝導が起きやすくなります。

2. 力 B：「壁を厚くする力（遮蔽効果）」

   • 圧力が高くなりすぎると、ビルを囲んでいる「壁（ラジウムイオン）」が電子の周りに集まり、電子同士の喧嘩を**「遮る（シールドする）」**ようになります。
   • 結果： 電子同士が互いの存在を認識しにくくなり、「喧嘩（反発）が弱まります。
   • これにより、超伝導が弱まります。

🌟 結論：

• 超伝導が始まる点： 「力 A（部屋を狭く）」が優勢になり、喧嘩が激しくなった瞬間。
• 超伝導が最強の点（頂点）： 「力 A」と「力 B」が絶妙なバランスになった瞬間（実験の 18 GPa と一致）。
• 超伝導が消える点： 「力 B（壁の遮蔽）」が勝ってしまい、喧嘩が弱まりすぎた瞬間。

この**「綱引きのバランス」**が、実験で見られた「右三角形」の超伝導領域を完璧に説明しているのです。

🌡️ 4. 温度の影響：「揺れるビル」

さらに、この研究は**「温度」**の影響もシミュレーションしました。

• 温度が上がると、原子は熱でブルブルと震えます（熱揺らぎ）。
• 実験では、高温になるほど超伝導が起きにくくなる（三角形の右側が斜めに下がる）ことが知られていますが、これは**「熱で原子が揺れると、電子の『喧嘩』が乱され、超伝導の秩序が崩れるから」**だと分かりました。
• 特に、圧力が高い領域では、熱揺らぎに対して構造が安定しているため、超伝導が長く続きますが、ある限界を超えると急激に弱まります。

🧪 5. 応用：「ラジウム」を「アクチニウム」に置き換える

最後に、研究者たちは「もし、ラジウム（La）という元素を、少し大きい元素（アクチニウム：Ac）に置き換えていたらどうなるか？」をシミュレーションしました。

• 予想： 大きい元素を入れると、圧力をかけなくても「まっすぐな正方形」の状態になりやすいはず。
• 結果： 確かに、「自然な圧力（化学的圧力）」だけで超伝導に必要な状態に近づきました。
• しかし： 電子の「喧嘩」が少し弱まりすぎたため、超伝導の温度（Tc）は、ラジウム版よりも少し低くなる（80 K → 70 K 程度）と予測されました。

これは、「常温超伝導」を目指す新しい材料の設計図として非常に重要です。「どの元素を混ぜれば、圧力なしで超伝導が起きるか」のヒントが得られたのです。

---

📝 まとめ

この論文は、La3Ni2O7 という物質が、**「圧力というスイッチ」によって、電子同士の「喧嘩の強さ」**を微妙に調整していることを発見しました。

• 圧力 0 → 14 GPa： 部屋が狭くなり、喧嘩が激しくなり、超伝導がスタート。
• 圧力 18 GPa： 喧嘩が最高潮に達し、超伝導が最強になる。
• 圧力 50 GPa 以上： 壁が厚くなりすぎて喧嘩が止まり、超伝導が消える。

この「綱引き」のバランスを数値で証明したことで、なぜ超伝導が三角形の領域に現れるのか、そしてどうすればより高い温度で超伝導を起こせるのかという**「設計図」**が描けたのです。これは、未来の「常温超伝導」を実現するための重要な一歩と言えます。

技術概要

以下は、提示された論文「First-Principles Evidence for Strongly Correlated Superconductivity Driven by Structural Variations in La3Ni2O7（構造変化に起因する La3Ni2O7 の強相関超伝導に対する第一原理的証拠）」の技術的サマリーです。

1. 研究の背景と課題 (Problem)

• 背景: 近年、層状ニッケル酸化物 La3Ni2O7 において、14 GPa 以上の静水圧下で 80 K までの高温超伝導が観測され、注目を集めています。
• 課題:
  • 超伝導発現領域（圧力 - 温度相図）は「右三角形」のドーム型を示し、臨界温度（Tc）は 18 GPa で最大となりますが、その背後にある電子相関と構造変化のメカニズムは未解明です。
  • 実験的には、室温から低温にかけての構造相転移（Cmcm 相、Fmmm 相、I4/mmm 相）の境界や、有限温度における構造の揺らぎが明確にされていません。
  • 従来の密度汎関数理論（DFT）計算だけでは、圧力変化に伴う電子相関の増大と Tc の相関、および「右三角形」ドームの起源を定量的に説明できていませんでした。

2. 手法 (Methodology)

本研究では、以下の第一原理計算手法を組み合わせ、多角的な解析を行いました。

• エンタルピー最小化による構造最適化: 0〜100 GPa の静水圧下で、La3Ni2O7 の結晶構造を最適化し、0 K における相図を構築しました。
• 多体摂動論と有効ハミルトニアンの導出:
  • 最適化された電子構造を基に、制約付きランダム位相近似（cRPA）を用いて、低エネルギー領域（Ni の eg 軌道）の有効モデル（二重層ドイマーモデル）を構築しました。
  • 最大局在ワニエ関数（MLWFs）を定義し、オンサイト反発力（U）、ホッピング積分（t）、およびそれらの比（U/t）を圧力依存性として算出しました。
• 第一原理分子動力学（AIMD）シミュレーション:
  • 10〜100 K の温度範囲で AIMD を実施し、有限温度効果（格子定数の揺らぎ、結合角の揺らぎ）を考慮した構造相図を再構築しました。
• A サイトカチオンの置換: 化学的圧力を評価するため、より大きなイオン半径を持つアクチニウム（Ac）でランタン（La）を置換した Ac3Ni2O7 についても同様の計算を行い、A サイトカチオンの役割を解明しました。

3. 主要な結果 (Key Results)

A. 構造相転移と電子相関の相関

• 構造変化: 約 10 GPa で、斜方晶（Cmcm 相）から正方晶（I4/mmm 相）への第一種相転移が発生します。これに伴い、NiO6 八面体の傾きが消失し、c 軸方向への整列（θ1 = 180°）と、ab 面内の対称性回復（a=b）が起こります。
• 電子相関の増大: この構造転移点付近で、有効オンサイト反発力（U）が約 30% 増大し、ホッピング積分（t）に対する比（U/t）が急激に上昇します。
• Tc との対応: 計算された U/t の圧力依存性は、実験的に観測された超伝導ドーム（右三角形）と驚くほど一致します。
  • U/t の最大値は 18 GPa 付近に現れ、これは実験的な最大 Tc（80 K）の圧力と一致します。
  • 低圧側と高圧側で U/t が低下する傾向は、超伝導が抑制される領域に対応します。

B. 相関増大のメカニズム（軌道局在とスクリーニングの競合）

U/t の圧力依存性は、以下の 2 つの競合する効果によって説明されます。

1. 相関増強要因: 八面体歪みの減少、バンド分裂、および MLWF の局在化により、有効相互作用（U）が増大します。
2. 相関弱体化要因（高圧側）: 高圧域（60 GPa 以上）では、スペーサー層である La 5d 軌道（5dx2-y2）がフェルミ準位付近に降りてきて、Ni eg 軌道を効果的にスクリーニングします。これにより、実効的な相互作用が減少し、超伝導が抑制されます。

• この「局在化による相互作用増大」と「スペーサー帯によるスクリーニング増大」の競合が、右三角形の超伝導ドームの形状を決定づけています。

C. 有限温度効果と「右三角形」ドームの起源

• AIMD 結果により、Cmcm 相から I4/mmm 相への転移境界は、温度に対してほぼ垂直（負の傾き）であることが示されました。
• 高温側では、熱揺らぎにより八面体の整列（θ1）や格子定数の差（a-b）がわずかに乱れますが、超伝導領域内では I4/mmm 相が安定に存在します。
• 高圧側での Tc の低下は、圧力と温度の両方によって intralayer 結合角（θ2）の「バッキング（たわみ）」が増大し、スペーサー La によるスクリーニングがさらに強化されることに起因すると結論付けられました。

D. Ac3Ni2O7 による化学的圧力の検証

• Ac3Ni2O7（Ac はより大きなイオン半径）では、化学的圧力により I4/mmm 相が常圧付近で安定化し、臨界圧力が低下します。
• しかし、Ac 置換により高圧域でのスクリーニング効果がさらに強まり、La3Ni2O7 に比べて最大 U/t 値が低下します。
• これは、Ac3Ni2O7 が常圧近傍で超伝導を示す可能性はあるものの、その最大 Tc は La3Ni2O7 よりも低くなる（約 70 K 程度）ことを予測しており、実験的・機械学習的予測と整合します。

4. 貢献と意義 (Significance)

• メカニズムの解明: La3Ni2O7 の超伝導が、単なるバンド構造の変化ではなく、「構造相転移に起因する電子相関の劇的な増大」と「スペーサー帯によるスクリーニングの競合」によって駆動されていることを、第一原理計算から定量的に証明しました。
• 相図の再構築: 有限温度効果を組み込んだ AIMD により、実験的に議論されていた構造相境界の位置を明確にし、超伝導ドームの「右三角形」形状の物理的起源を解明しました。
• 材料設計への指針: A サイトカチオンのイオン半径が、化学的圧力を通じて電子相関とスクリーニングのバランスを制御し、超伝導特性（臨界圧力と Tc）を決定づけることを示しました。これは、常圧超伝導ニッケル酸化物の探索における重要な設計指針となります。

この研究は、強相関電子系における構造と電子状態の密接な関係を明らかにし、高温超伝導のメカニズム理解と新材料開発に重要な貢献を果たしています。