Hubbard-$U$-corrected electron-phonon interactions in strongly correlated materials via the finite-displacement method

要約

活気あふれる都市を想像してください。そこでは常に二種類の交通渋滞が発生しています。一つは車同士が衝突すること（電子が他の電子を反発すること）によるもので、もう一つは車が道路の穴（ポットホール）に突っ込むこと（電子が振動する地面、すなわち「フォノン」に衝突すること）によるものです。超伝導体、つまりゼロ抵抗で電気を伝導する物質の世界において、科学者たちはこの二つの渋滞がどのように相互作用するかを知りたがっています。もしこれらが完璧に調和すれば、都市は交通が完全に流れる「スーパーハイウェイ」を実現できるのです。

長年、科学者たちはこれらの物質を研究するためにDFT（密度汎関数理論）と呼ばれる標準的な地図作成ツールを用いてきました。しかし、「強相関」物質（車同士が非常に攻撃的で絶えず衝突し合っている状況）において、この標準的な地図はしばしば不正確です。これを修正するため、科学者たちは地図にハバード Uと呼ばれる「補正因子」を追加しました。

問題は、科学者たちがこの補正を車（電子）に適用する方法は知っていたものの、ポットホール（フォノン）やそれらの間の衝突（電子 - フォノン結合）にどう適用すべきか分からなかったことです。彼らは車の地図を修正しながらも、攻撃的な運転によってポットホール自体の形状が変わるという事実を無視していました。

新しいアルゴリズム：完全なリモデル
この論文は、「ハバード U」補正をすべて、すなわち車、ポットホール、そしてそれらの間の衝突に適用する新しい手法（アルゴリズム）を導入します。これは単に車を修理するだけでなく、道路を再舗装し、交通規則を再設計して、すべてが整合性を取るよう同時に実行する建設チームのようなものです。

研究者たちはこの新しい「完全なリモデル」を、2 つの特定の物質でテストしました。

1. ニッケレート市（LaNiO₂）

• 謎: この物質は非常に低温で超伝導体となります。最近の研究では、異なる非常に高価な手法（GW 法）を用いて、「車とポットホールの間の衝突」が標準的な地図が予測したものの 5 倍も巨大であることが示唆されました。これは、衝突が超伝導性の主な原因であることを意味していました。
• 論文の発見: 著者たちが新しい「完全なリモデル」（DFT+U）を使用すると、衝突は依然として小さいことが分かりました。
• 比喩: GW 法が「車はポットホールに激しく衝突して、道路全体を揺らしている！」と言ったとすると、新しい手法は「実際には、車は普通に走行しているだけだ」と言います。
• なぜ違いが生じたのか? 二つの手法は都市の配置（フェルミ面）を異なって描きました。GW 法は、車が狭い角に追いやられ、大規模な衝突を引き起こすような配置を描きました。一方、新しい手法は、車に十分な移動の余地があり、衝突が軽微なままになるような配置を描きました。著者たちは、この物質において「衝突」だけでは超伝導性を説明するには弱すぎると結論付け、現象を駆動しているのは別の要因であると示唆しています。

2. ルテニウム二酸化物市（RuO₂）

• 謎: この物質は特定の基板上に成長させた薄膜です。実験では超伝導体となることが示されていますが、非常に低い温度（1.5 ケルビン）に限られます。しかし、標準的な地図（通常の DFT）は破滅を予言しました。道路が不安定すぎて崩壊する（虚数のフォノン・モード）と言われ、衝突があまりに暴力的であるため、都市ははるかに高い温度（30 ケルビン）で超伝導体になるべきだとしました。
• 論文の発見: 「完全なリモデル」（フォノンと衝突の両方にハバード U を追加）を適用すると、二つのことが起こりました。
  1. 道路が安定した: 「崩壊する道路」（虚数モード）は消えました。道路は現実世界で観察されるように、堅固で安定したものになりました。
  2. 衝突が鎮静化した: 暴力的な衝突は穏やかな衝突に変わりました。「衝突エネルギー」の総量は大幅に減少しました。
• 結果: これが、超伝導性がなぜそれほど弱い（低温である）のかを説明します。「補正」は道路を硬化（フォノンの硬化）させ、車が衝突しにくくしました。これは実験的な現実と完璧に一致します。

大きな教訓
この論文は、「車」（電子）を修正するならば、「道路」（フォノン）と「交通規則」（結合）も同時に修正しなければならないと主張しています。

• もし車だけを修正する（「部分的」な補正）と、間違った答えを得る可能性があります。ルテニウム二酸化物の場合、部分的な修正は現実には存在しない超強力な超伝導体を予言することになります。
• 著者たちは、ニッケレートのような物質では、補正は配置をわずかに変えるものの、結果にはあまり影響を与えないことを示しています。一方、ルテニウム二酸化物のような物質では、道路の崩壊を防ぎ、超伝導性がなぜそれほど弱いのかを説明するために、この補正が不可欠です。

要約すれば、この論文は複雑な物質における電子と振動の相互作用をより一貫性を持って描く新しい方法を提供し、振動自体に対する「相関」効果を無視することは誤った予測につながることを示しています。

技術概要

技術的概要：有限変位法による強相関物質におけるハバード U 補正を施した電子 - 格子相互作用

問題提起
電子 - 電子相互作用と電子 - 格子相互作用の相互作用は、強相関物質の理解において中心的な役割を果たす。これらの系における電子状態の計算には、密度汎関数理論にハバード U を加えた手法（DFT+U）が標準的に用いられているが、電子 - 格子結合（g 行列）および格子振動スペクトルへの拡張は限定的であった。従来のアプローチでは、電子 - 電子相互作用と電子 - 格子相互作用を別々に扱い、格子振動には通常の DFT を、相関には経験的パラメータを用いることが多かった。銅酸化物やニッケル酸化物などの物質において電子 - 格子計算への完全な GW 補正が適用された例もあるが、計算コストは複雑な実在物質に対しては依然として許容しがたいほど高い。さらに、これらの系における相関効果の重要性は既知であるにもかかわらず、ハバード U 補正が格子振動スペクトルおよび電子 - 格子 g 行列に与える具体的な影響は体系的に調査されていない。

手法
著者らは、DFT+U 法と有限変位法を統合し、格子振動スペクトルおよび電子 - 格子 g 行列を計算するアルゴリズムを実装した。このアプローチは、特に大きな U 値が適用される場合に頻繁に遭遇する、ハバード U を加えた密度汎関数摂動論（DFPT+U）計算における収束問題を回避する。

主要な手法的特徴は以下の通りである：

• 有限変位アプローチ: この手法は、実空間における格子振動摂動ポテンシャルを得るために超格子 DFT+U 計算を利用し、これをフーリエ変換することで運動量空間における電子 - 格子 g 行列を導出する。
• 基底関数: 計算は、平面波基底に比べて大規模な超格子シミュレーションにおいて計算効率が高い数値原子軌道（NAO）基底関数を用いた OpenMX パッケージで行われる。
• 完全 U 補正フレームワーク: 決定的な点として、ハバード U 補正は、電子状態（バンド構造）、格子振動スペクトル（力定数）、および電子 - 格子 g 行列の 3 つの構成要素に対して一貫して適用される。これは、相関を電子状態のみに適用し、格子振動や結合行列は通常の DFT レベルのままにする「部分的 U 補正」アプローチとは区別される。
• 対象物質: このアルゴリズムは、2 つの代表的な相関物質、すなわち無限層ニッケル酸化物（特に 20% ホールドープされた LaNiO$_2$）およびルテニウム酸化物（RuO$_2$、TiO$_2$ 基板上の歪んだ RuO$_2$ 薄膜を含む）に適用された。

主要な結果

• 無限層ニッケル酸化物（LaNiO$_2$）:

  • 電子状態: Ni-$d$ 軌道に対するハバード U 補正は、主に他の Ni-$d$ バンドをフェルミ準位から押しやるが、主要な Ni-$d_{x^2-y^2}$ バンドはほぼ固定されたまま残る。その結果、フェルミ面は中程度に変化し、小さなゾーン中心の電子ポケットを誘起するが、支配的な Ni-$d_{x^2-y^2}$ シートは維持される。
  • 格子振動と結合: 格子振動スペクトルは力学的に安定しており、U 補正に対してほとんど感応しない。したがって、全電子 - 格子結合（$\lambda$）はわずかに増加するのみである（U=0 eV で 0.18 から U=7 eV で 0.28 へ）。
  • 超伝導: 完全な U 補正を施しても、推定される超伝導転移温度（$T_c$）は 3 K 未満のままであり、実験値である 10–30 K よりもはるかに低い。
  • GW 計算との比較: これらの結果は、電子 - 格子結合が 5 倍に増加すると予測する最近の完全 GW 研究とは対照的である。著者らは、この不一致の原因を、DFT+U が予測するフェルミ面トポロジー（Ni-$d_{x^2-y^2}$ が支配的）と、La-$d$ 由来のシートを含む再構成を予測する GW との間の違いに帰着させる。

• ルテニウム酸化物（RuO$_2$）:

  • 安定性: 通常の DFT 計算は、実験的に観測された歪んだ RuO$_2$ の直方体 $Pnnm$ 構造に対して、虚数の格子振動モード（格子不安定性）を予測する。Ru-$d$ 軌道に対するハバード U 補正の導入は、これらの虚数モードを除去し、構造を力動的に安定化させる。
  • 格子振動の硬化: U の増加は、特に低周波モードにおいて顕著な格子振動の硬化を引き起こす。
  • 結合の減少: この格子振動の硬化は、全電子 - 格子結合を著しく減少させ、$\lambda$ を（U=1 eV で 0.8 から U=5 eV で 0.2 へ）低下させる。
  • 超伝導: 結合の減少により、推定される $T_c$ は（低 U で）約 30 K から数ケルビンへと低下し、実験的に観測される約 1.5 K と一致する。これは、通常の DFT が予測する大きな結合と、低い実験的 $T_c$ の間の不一致を解決する。
  • 部分補正と完全補正: 「部分的 U 補正」アプローチ（DFT 格子振動を使用）は RuO$_2$ に対して失敗し、虚数モードの存在により物理的に不適切な大きな結合をもたらす。これは、格子振動スペクトル自体に相関補正を適用する必要性を浮き彫りにする。

意義と主張
本論文は、電子状態だけでなく、格子振動スペクトルおよび電子 - 格子結合行列においてもハバード U 補正を取り入れた、相関物質における電子 - 格子相互作用を計算するための完全な自己無撞着アルゴリズムを提供すると主張している。

著者らは、相関が電子 - 格子結合に影響を与える 2 つの異なるメカニズムを強調していると述べている：

1. 電子チャネル: フェルミ面トポロジーと散乱位相空間を変更する（LaNiO$_2$ では minor な効果として観測）。
2. 格子振動チャネル: 格子ダイナミクスを直接再規格化し、格子振動の硬化または軟化をもたらす（RuO$_2$ では支配的な効果として観測）。

本研究は、相関物質における電子 - 格子相互作用の信頼性の高い記述には、電子および格子のすべての自由度に一貫して相関効果を含める必要があると結論付けている。格子振動セクターにおけるこれらの効果を無視することは、偽の格子不安定性や超伝導転移温度の過大な過大評価といった、質的に誤った予測につながる可能性がある。この研究は、複雑な実在物質におけるこれらの性質を調査するための、完全な GW 法に代わる計算効率の高い代替手段を提供する。