Hybrid quantum-classical matrix-product state and Lanczos methods for electron-phonon systems with strong electronic correlations: Application to disordered systems coupled to Einstein phonons

この論文は、電子相関を数値的に厳密に扱い光学フォノンを古典的に近似するハイブリッド量子古典手法（時間依存ランチョス法と行列積状態法）を開発し、強い電子-フォノン結合が乱雑な系における多体局在の不安定化と非局在化を引き起こすことを示した。

要約

1. 舞台設定：電子と振動する床

まず、物質の中を想像してください。

• 電子：部屋の中を飛び回る「元気な子供たち」。
• 格子（原子）：子供たちが立つ「床」や「壁」。
• フォノン（格子振動）：床が揺れたり、壁が振動したりすること。

通常、電子は床の上を滑らかに走りますが、この研究では**「床が激しく揺れている（強い振動）」かつ「子供同士がぶつかり合っている（強い相互作用）」**という、とてもカオスな状況を扱っています。

さらに、この部屋には**「ランダムに置かれた障害物（不純物）」**もあります。

• 障害物がある状態：子供たちが自由に動けず、特定の場所に**「閉じ込められてしまう（局在化）」**現象が起きます。これを「アインシュタインの局在」や「多体局在（MBL）」と呼びます。

2. 問題点：計算が難しすぎる！

このカオスな状況を正確に計算しようとするとき、従来の方法には大きな壁がありました。

• 量子力学の壁：電子は「波」のような性質を持つため、すべての可能性を同時に計算する必要があります。子供が 100 人いれば、その組み合わせは天文学的な数になり、スーパーコンピューターでも計算しきれません。
• 振動の壁：床の振動も量子力学で扱うと計算量が爆発します。

3. 解決策：新しい「ハイブリッド・シミュレーション」

そこで、この論文の著者たちは**「量子と古典を混ぜ合わせた新しい方法」**を開発しました。

① 電子は「量子」、床は「古典」

• 電子（子供たち）：まだ「量子力学」のルール（確率や波）に従って厳密に計算します。
• 床（振動）：「古典力学」のルール（ボールが転がるような単純な動き）として扱います。
  • 例え：子供たちの動きは複雑な魔法（量子）で計算しますが、床の揺れはただの「揺れる板」として扱います。

② 「多経路エレンフェスト法」というテクニック

床の揺れを正確に表現するために、**「何百回も同じ実験を繰り返す」**という手法を使います。

• イメージ：
  1. 床の揺れ方を「パターン A」「パターン B」「パターン C」……と何百通りもランダムに決めます。
  2. それぞれのパターンに対して、電子（子供たち）がどう動くかを計算します。
  3. 最後に、すべての結果を**「平均」**して、本当の答えを導き出します。

これを**「ランチョス法」と「MPS（行列積状態）」**という 2 つの異なる計算テクニックで行い、両方が同じ結果を出すことを確認しました。これは、新しい計算方法が正しいことを証明する「ベンチマーク（基準）」です。

4. 発見：「揺れ」が閉じ込めを解き放つ

この新しい方法を使って、実際に「閉じ込められた電子」がどうなるかを調べました。

• 実験：

  • 障害物（不純物）だらけの部屋で、電子が動けなくなっている状態（局在化）を作ります。
  • そこに、床を**「揺らす（電子と振動を結合させる）」**操作を加えます。

• 結果：

  • 床が揺れると、「閉じ込められていた電子が、少しずつ動き出し始めました！」
  • 完全に自由になったわけではありませんが、「じっとしている」状態から「ゆっくりと広がる（亜拡散）」状態に変わりました。

• なぜ？（メカニズムの解説）

  • 床が揺れると、一時的に「障害物」の位置や強さが変わります。
  • 電子にとっては、「あ、今なら通り抜けられる！」という瞬間が生まれるのです。
  • 電子は電子同士でぶつかり合いながら、この「揺らぎ」に乗って、少しずつ逃げ出していくのです。

5. この研究の意義

• 技術面：「量子」と「古典」を賢く混ぜることで、以前は計算できなかった複雑な物質の動きを、現実的な時間でシミュレーションできるようになりました。
• 物理面：「電子が完全に止まっている（絶縁体）」状態でも、「振動（熱や音）」が加わると、電子が動き出す（導電性を持つ）可能性があることを示しました。
  • これは、「多体局在（MBL）」という、非常に安定した状態が、振動によって不安定になることを意味します。

まとめ

この論文は、**「複雑な電子の動きを、床の揺れを『古典的な揺れ』として扱いながら、何百回もシミュレーションを繰り返すことで解き明かす」**という新しいアプローチを紹介しました。

その結果、**「電子が閉じ込められていた場所でも、床が揺れることで、電子が少しずつ逃げ出すことができる」**という、物質の新しい性質を見つけ出しました。これは、将来の新しい電子デバイスやエネルギー変換材料の開発につながる重要な発見です。

技術概要

この論文は、強い電子相関と電子 - 格子相互作用（フォノン）を同時に扱うためのハイブリッド量子・古典シミュレーション手法を提案し、乱れ（不純物）が存在する系における電荷密度波（CDW）秩序の崩壊と多体局在（MBL）の安定性について研究したものです。

以下に、問題設定、手法、主要な貢献、結果、および意義について詳細な技術的サマリーを記します。

1. 問題設定 (Problem)

固体物理学や量子化学において、電子 - 格子結合は超伝導やポラロン形成など多くの現象で重要ですが、電子間の強い相関（電子 - 電子相互作用）と電子 - 格子結合、さらに乱れ（disorder）が共存する系の非平衡ダイナミクスを正確に記述することは理論的に極めて困難です。

• 既存手法の限界: 従来の行列積状態（MPS）やダイナミカル・メアンフィールド理論（DMFT）などの手法は、フォノンの量子性を厳密に扱う際に、特にフォノンエネルギーが帯域幅より小さい「断熱的（adiabatic）」な領域において計算コストが爆発的に増大する、あるいは収束が困難になるという課題を抱えています。
• 研究の目的: 電子相関を数値的に厳密に扱いながら、光学フォノンの自由度を古典的に近似することで、長時間・大規模な系の時間発展をシミュレーション可能な手法を開発し、乱れのある電子 - フォノン系における局在の安定性を検証すること。

2. 手法 (Methodology)

著者らは、電子系を量子力学、フォノン系を古典力学として扱うマルチ・トラジェクトリ・エレンフェスト（MTE: Multi-Trajectory Ehrenfest）法を基盤とし、これに 2 つの異なる量子多体手法を組み合わせるハイブリッド手法を提案しました。

• 基本アプローチ (MTE):
  • フォノンの運動を古典的な軌道（座標 $x_\ell$、運動量 $p_\ell$）として扱い、電子系はこれらの軌道に依存する時間依存ハミルトニアン下で量子力学に従って進化させます。
  • 量子的不確実性を考慮するため、初期状態のフォノン座標・運動量をウィグナー関数からサンプリングし、多数の独立した軌道（トラジェクトリ）をシミュレーションして物理量を平均化します。
• 提案された 2 つのハイブリッド手法:
  1. Lanczos-MTE 法: 時間依存ランチョス法（Krylov 部分空間法）を用いて電子系の時間発展を計算。有限サイズ系（1 次元）において電子相関を厳密に扱えます。長時間シミュレーションに優れます。
  2. TEBD-MTE 法: 時間発展ブロックデシメーション（TEBD）と行列積状態（MPS）を用いた手法。より大きな系サイズを扱えますが、エンタングルメントの増大により長時間シミュレーションには限界があります。
• モデル:
  • 1 次元格子におけるスピンレス・フェルミオン（t-V モデル、Heisenberg 鎖と等価）。
  • ホルシュタイン型電子 - フォノン結合（局所フォノン）。
  • 乱れ項（ランダムなオンサイトポテンシャル）。
  • 初期状態：電荷密度波（CDW）状態または単一粒子の局在状態。

3. 主要な貢献 (Key Contributions)

• 手法の確立とベンチマーク:
  • 電子相関を厳密に扱う Lanczos-MTE および TEBD-MTE 手法を実装し、非相互作用系（単一電子）において既存の MTE 手法や厳密対角化の結果と一致することを確認しました。
  • 相互作用系においても、ランチョス誤差や MPS の切断誤差（discarded weight）を制御することで、統計的誤差（トラジェクトリ平均）が支配的になるまで精度を向上させることを示しました。
• 断熱的極限における有効性の確認:
  • フォノン周波数 $\omega_0$ が電子ホッピング $t_0$ よりも十分に小さい（$\omega_0 \ll t_0$）断熱的領域において、これらのハイブリッド手法が信頼できる結果を与えることを示しました。

4. 結果 (Results)

シミュレーションは、乱れ（$W$）と電子 - フォノン結合強度（$\gamma$）を変化させて行われました。

• 単一電子の拡散（Anderson 局在の不安定性）:
  • 電子間相互作用がない場合（$V=0$）、$\gamma=0$ では Anderson 局在が観測されますが、$\gamma > 0$ とすると粒子は局在から解放され、時間とともに広がります。
  • 拡散のダイナミクスは**サブ拡散（sub-diffusive）**であり、平均二乗変位 $\sigma^2(t) \propto t^z$ の指数 $z$ は $0 < z < 1$（約 0.5 以下）となることが示されました。これは、古典的なフォノン浴との結合が局在を破壊し、遅い拡散運動を引き起こすことを意味します。
• 多体局在（MBL）の不安定性:
  • 電子間相互作用がある場合（$V>0$）、$\gamma=0$ では有限サイズ系において MBL が観測されます。
  • しかし、$\gamma > 0$ とすると、CDW 秩序パラメータ（不均衡 $I_{CDW}$）が時間とともに減衰し、MBL が不安定化することが示されました。
  • 減衰の特性: 秩序の崩壊はべき乗則に従います。
    • 弱い結合領域（$\gamma \ll t_0$）では、電子間相互作用が強いほど減衰が速くなります（相互作用が緩和を促進）。
    • 強い結合領域（$\gamma \gtrsim t_0$）では、ポラロン形成により粒子の移動度が低下し、逆に緩和が遅くなる傾向が見られました。
• 物理的メカニズム:
  • 古典フォノンの揺らぎが動的な乱れポテンシャルを生成し、隣接サイト間のエネルギー差を瞬間的にゼロにする（共鳴条件を満たす）ことで、電子が局在から脱出する「局所的な緩和プロセス」が駆動されると解釈されました。

5. 意義と結論 (Significance and Conclusion)

• 理論的意義: 電子相関、電子 - フォノン結合、乱れが共存する非平衡系のダイナミクスを研究するための、計算コストと精度のバランスの取れた新しい数値枠組みを提供しました。特に、断熱的極限における古典フォノン近似の有効性を、電子相関を厳密に扱う手法と組み合わせることで実証しました。
• 物理的知見:
  • 古典的なフォノン浴（熱浴）との結合は、1 次元の乱れ系における Anderson 局在だけでなく、多体局在（MBL）さえも不安定化し、サブ拡散的な輸送を引き起こすことを示しました。
  • これは、MBL が「完全な孤立系」でのみ安定であり、環境との結合（特に古典的な振動子との結合）によって容易に壊れる可能性を示唆しています。
• 将来展望:
  • TEBD-MTE 法は、平衡状態のスペクトル関数や有限温度での動的性質の計算への応用が期待されます。
  • 将来的には、より高度な軌道法（多配置エレンフェスト法など）と MPS を組み合わせることで、MTE の限界（過剰なコヒーレンスなど）を克服し、より正確な非断熱的ダイナミクスの記述が可能になると考えられています。

この研究は、強相関電子系における非平衡現象を理解するための重要なステップであり、実験的に観測される光励起後の緩和過程や、乱れのある物質中の輸送現象の解釈に寄与するものです。