Electronic correlations and fluctuating lattice distortions in vanadium dioxide

技術要約：二酸化バナジウムにおける電子相関と格子歪みの揺らぎ

問題提起
二酸化バナジウム（VO$_2$）のような強相関物質における金属－絶縁体転移（MIT）は、しばしば構造変化を伴う。従来の理論的記述は、「コヒーレント格子近似」に依存することが多い。これは、格子をエネルギー最小化によって固定された、単一で決定論的な空間的に均質な歪み座標として表現するものである。この枠組みでは、構造の自由度は本質的なエントロピーを持たないが、転移エントロピーの電子系と格子系への分配が転移温度を決定づける決定的な要因であるという証拠が存在する。例えばVO$_2$において、高温相であるルチル相の熱力学的安定化は大きな振動エントロピーに起因するとされる一方で、他の推定では転移エントロピーを電子系に帰する議論もある。さらに、ナノスケールの共存や超高速な無秩序化の実験的観察は、静的な平均場格子記述では捉えきれない格子揺らぎや非断熱効果が重要な役割を果たしていることを示唆している。課題は、相関電子と揺らぐ格子歪みが、特に複数の結合した格子モードや非線形な電子－フォノン相互作用を伴う系において、等価な立場で同時に進化する枠組みを開発することにある。

手法
著者らは、相関電子と揺らぐ格子歪みの結合ダイナミクスに対処するため、動力学的平均場理論（DMFT）の確率的半古典的拡張を開発した。この手法は以下の主要な要素で構成される：

1. モデル系： VO$_2$の金属－絶縁体転移の本質的な物理を捉えるために、最小限の二軌道モデルを構築した。このモデルには以下が含まれる：

   • 電子の局所的なハバード反発（$U$）を持つ2つの電子バンド。
   • モノクリニックM1相の対称性の異なる歪み（二量化 $X_1$ および横方向の傾き $X_2$）を表す2つの局所格子座標（$X_1, X_2$）。
   • 結合： 二量化モードは電子密度と線形に結合し（電荷密度波不安定性を駆動）、傾きモードは軌道分極と二次的に結合する（結晶場として機能）。
   • 格子ポテンシャル： 互いに置換的な二量化を促進するための、サイト間フォノン結合を含む不適切強誘電体の現象論的ポテンシャルを使用。

2. 理論的枠組み：

   • 不均一DMFT： 系を、無限の配位数（$Z \to \infty$）を持つベテ格子上のインパリティ問題のアンサンブルにマッピングする。これにより、完全に不均一な格子問題を明示的に解くことなく、空間的に揺らぐ格子構成を扱うことが可能になる。
   • 確率的半古典近似（SCA）： 格子座標は確率的な古典的軌跡（$X_{cl}(t)$）として扱われ、電子の自由度は量子力学的に扱われる。
   • 非保存型ランジュバン方程式： 電子を積分消去し、量子的な格子場の影響関数を量子成分の2次まで展開することにより、格子の運動に関するランジュバン方程式を導出する。これらの方程式には以下が含まれる：
     • 決定論的な力（電子系からのヘルマン－ファインマン力）。
     • 電子による減衰（摩擦）。
     • 相関のある確率的ノイズ。
   • 自己整合性： 摩擦およびノイズ行列は、電子のグリーン関数および応答関数（遅延およびKeldysh成分）によって自己整合的に決定される。この定式化は、単一モード・線形結合の先行研究を、多モード・非線形結合へと拡張するものである。
   • 非平衡ダイナミクス： 電子系は量子ボルツマン方程式（QBE）を用いて発展させ、非熱的な電子分布を扱うことが可能であるが、本研究では平衡状態に焦点を当てる。

主な貢献

• 定式化の一般化： 本論文は、単一モード・線形電子－フォノン結合から、多モード・非線形結合を持つ系へと、確率的半古典的定式化を拡張した。これにより、電子の摩擦とノイズが瞬時の格子歪みに依存する、構成依存型の頂点（バーテックス）の取り扱いが可能となった。
• 統一された枠組み： 絶縁体、中間的な金属相、および高対称金属相を明示的に記述でき、電子相関と格子エントロピーの相互作用を考慮できる単一の理論的枠組みを提供する。
• 転移の微視的起源： このアプローチは、構造転移を駆動する電子エントロピーと格子エントロピーの役割を分離し、それらが異なる構造自由度に対して異なる温度スケールで作用することを示している。

結果
最小限のVO$_2$モデルにこの手法を適用した結果、温度上昇に伴い3つの明確な相を持つ平衡相図が得られた：

1. 絶縁体相（低温）： 二量化（$X_1$）と傾き（$X_2$）の両方の歪みが有限である。互置的な二量化が下部バンドに絶縁ギャップを開き、傾きが軌道分極を強化する。系は常磁性絶縁体である。
2. 中間低対称金属相（中間温度）： 二量化が特定の温度（モデル単位で $T \approx 0.4$）で消失し、下部バンドにおける金属的挙動を回復させる。しかし、傾き歪みは有意に高い温度（$T \approx 1.45$）まで有限に留まる。この相は、残留する軌道分裂を伴う金属状態であり、実験的に議論されている「モノクリニック金属」に類似している。
3. 高対称金属相（高温）： 傾き歪みが最終的に消失し、歪みのない高対称なルチル相が回復する。

転移メカニズムに関する主要な知見：

• 二量化の融解： 絶縁体から中間金属相への転移は、主に相関電子への結合によって駆動される。この転移は、格子揺らぎを無視するコヒーレント格子近似においても定性的に捉えられるが、確率的アプローチによる揺らぎの導入は、転移温度を定量的に再正規化（約0.7倍に低下）させる。
• 傾きの回復： 中間金属相から高対称金属への転移は、電子による平均場自由エネルギーのみでは説明できない。なぜなら、コヒーレント格子近似はすべての温度において有限の傾き歪みを予測するためである。代わりに、この第2の転移は揺らぐ格子のエントロピーによって駆動される。高温では、格子揺らぎが十分に重要となり、残存する傾いた状態を不安定化させる。
• 電子によるフィードバック： 熱平衡状態であっても、格子に作用する電子の減衰とノイズは、瞬時の格子構成に応じて強い軌跡間の変動を示す。これは、軌跡依存の減衰とノイズを保持することが、単一の現象論的な値を使用することよりも不可欠であることを示唆している。

意義と主張
本論文は、非線形結合を持つ相関電子－格子系に対する確率的半古典的枠組みを確立したと主張している。その主要な意義は、電子相関と格子エントロピーが異なる構造自由度に対して作用し、相図を定性的に作り変えることができることを示した点にある。著者らは、コヒーレント格子近似が中間金属相の出現を捉えることはできるものの、高対称相の完全な回復には、格子揺らぎとその随伴するエントロピーの導入が必要であることを示している。

著者らは、現在の研究は平衡状態に焦点を当てているものの、非平衡DMFTとQBEを使用しているため、本定式化が非熱的な電子分布を持つ駆動・光励起系に直接適用可能であることを控えめに述べている。彼らは、VO$_2$における超高速な光誘起絶縁体－金属転移の微視的なシミュレーションを自然な次の応用として挙げており、そこでは本枠組みを用いて、急速な格子無秩序化や過渡的な不均一無秩序状態を調査できるとしている。本論文は、実在のVO$_2$の全微視的詳細を解決することを目的としたものではなく、むしろ最小限のモデルにおける競合する要素間の本質的な相互作用を捉えることを目的としている。