Thermalisation as Diffusion in Hilbert Space

この論文は、標準的なマルコフ近似やフェルミ黄金律を超えて、相互作用行列要素を確率変数としてモデル化することで、ヒルベルト空間内での拡散過程として熱化を記述する微視的理論を構築し、非マルコフ的な緩和時スケールと平衡条件を導出したことを報告しています。

要約

この論文は、**「量子の世界で、なぜ物が温まったり冷めたりするのか（熱化）」**という、物理学の大きな謎を解き明かそうとする挑戦的な研究です。

従来の物理学では、「熱くなるプロセス」は単純な「お風呂にお湯が混ざっていく」ようなイメージ（マルコフ過程）で説明されることが多かったのですが、この論文は**「もっと複雑で、予測不能な世界」**を扱っています。

以下に、専門用語を排し、日常の例え話を使ってこの研究の核心を解説します。

---

1. 物語の舞台：小さな温度計と巨大な海

まず、この研究のモデルを想像してください。

• 温度計（サーモメーター）: 小さな量子システム（例えば、たった 1 つの電子スピン）。
• 海（バス）: 温度計を取り囲む、巨大で複雑な量子の海（多くの粒子が絡み合っている状態）。

私たちが普段見ている「温度」は、この小さな温度計が、巨大な海とエネルギーをやり取りして、最終的に海と同じ温度になる現象です。

2. 従来の考え方 vs 新しい発見

従来の考え方（お風呂の例え）

昔の理論は、温度計がお風呂（海）に入ると、お湯の分子が次々とぶつかり、**「一定の速さ」**で温度計の温度が変わると考えていました。これは「フェルミの黄金律」と呼ばれるルールで、お湯が混ざり合う様子は非常にスムーズで予測しやすいものでした。

この論文の発見（迷路の例え）

しかし、この論文は**「お湯が混ざらない、あるいは非常に複雑に絡み合う世界」を扱っています。
ここでは、温度計が海とエネルギーをやり取りする様子は、「巨大な迷路を歩く」**ことに似ています。

• 迷路の壁（エネルギー準位）: 温度計が移動できる場所。
• 迷路の道（相互作用）: 温度計が次の場所へ移動するルート。

この研究では、迷路の道が**「ランダム」で、しかも「太さがバラバラ」**（一部は極太、一部は極細）であることが重要だと指摘しています。

3. 核心：ヒルベルト空間での「拡散」

この論文の最大の特徴は、**「熱化とは、ヒルベルト空間（量子の状態のすべて）における『拡散』である」**と捉え直した点です。

• ヒルベルト空間とは？: 温度計が取りうる「すべての可能性（状態）」の地図です。
• 拡散（Diffusion）: 温度計が、この地図の上を、ランダムに、しかしある法則に従って「ふらふらと歩き回る」こと。

「なぜ熱くなるのか？」
温度計が迷路を歩き回り、最終的に「最も歩きやすい場所（エネルギー的に安定な場所）」に落ち着くまで、時間がかかります。この**「歩き回る速さ」を決めるのが、「レベルの広がり（ブロードニング）」**という概念です。

• レベルの広がり ＝ 「迷路の道がどれだけ広く開けているか」
  • 道が広い（相互作用が強い） → すぐに次の部屋に行ける → 熱化が速い。
  • 道が狭い、あるいは分岐が複雑 → 行き詰まる → 熱化が遅い、あるいは止まる。

4. 重要な発見：「平均」では測れない世界

従来の理論は「平均的な道幅」を見て計算していましたが、この論文は**「道幅のバラつき（分布）」**そのものが重要だと説いています。

• レヴィ分布（Lévy distribution）の例え:
  道幅が「ほとんどが細い道」なのに、**「たまに超巨大なトンネル」**が現れるような世界を想定しています。
  • 従来の理論：「平均的な道幅」で計算する → 実際にはトンネルを通れるので、計算と現実にズレが生じる。
  • この論文の理論：「トンネルの存在」を含めた「道幅の全体的な分布」を考慮する → 実際の熱化の速さを正確に予測できる。

つまり、「稀に起こる巨大な出来事（トンネル）」が、全体の熱化のスピードを支配しているというのです。

5. 実験的な検証：3 つのシミュレーション

この理論が正しいか確認するために、著者は 3 つの異なる「迷路」でシミュレーションを行いました。

1. レヴィモデル: 道幅が極端に偏っている（太いトンネルがたまにある）迷路。
2. TFIM（トランスファード・イジングモデル）: すべてがつながった、カオスな迷路。
3. Imbrie モデル: 1 次元の、少し乱れた迷路（不純物がある）。

結果、**「迷路を歩き回る（拡散する）理論」と、「実際に迷路を歩いた（数値計算した）結果」**が、驚くほど一致しました。特に、従来の理論が破綻する「複雑で不規則な世界」でも、この新しい「拡散の理論」は正しく機能しました。

6. 結論：何がすごいのか？

この論文は、**「量子の世界での熱化は、単純なお湯の混ざり合いではなく、複雑な迷路を歩くような『拡散』のプロセスである」**ことを示しました。

• 従来の限界: 「平均的な速さ」でしか説明できなかった。
• この論文の貢献: 「道幅のバラつき」や「稀な巨大トンネル」を含めた「分布」を考慮することで、**「マルコフ的（単純な）ではない、より複雑な熱化」**を説明できるようになった。

これは、「乱れた物質（不純物が多い物質）」や「量子コンピュータの誤り」、あるいは**「なぜ一部の物質は熱くなりにくいのか（多体局在）」**といった、現代物理学の難問を解くための新しい「地図」となると期待されています。

---

一言でまとめると：
「熱くなるプロセスは、単純なお湯の混ざり合いではなく、**『道幅がバラバラな巨大迷路』を、その迷路の特性（広がり）に合わせてふらふらと歩き回る『拡散』**だと考えれば、複雑な量子現象も説明できるよ！」という新しい視点を提供した研究です。

技術概要

以下は、提示された論文「Thermalisation as Diffusion in Hilbert Space（ヒルベルト空間における拡散としての熱化）」の技術的な要約です。

1. 研究の背景と問題意識

閉じた量子系における熱化（thermalisation）の理解は現代物理学の中心的な課題の一つです。

• 従来のアプローチ: 固有状態熱化仮説（ETH）は、長時間極限での局所観測量の熱的値への収束を説明しますが、熱化の時間スケールを制御された形で予測することは困難です。また、有限サイズの数値計算から分布パラメータを抽出するのも難しい場合があります。
• 課題: 強い乱れを持つ相互作用系（多体局在：MBL の近傍など）では、ヒルベルト空間内の局在効果が重要となり、マルコフ近似やフェルミ黄金則（FGR）に基づくマスター方程式（リンドブラッド方程式）が破綻します。特に、相互作用行列要素の統計分布が広範で非ガウス的（例：レヴィ分布）な場合、緩和率の分布が広がり、従来の記述では実験的に観測される実空間の緩和とヒルベルト空間内の遷移を定量的に結びつけることができていません。
• 本研究の目的: マルコフ近似や FGR に依存しない、より一般的な熱化の微視的理論を構築し、ヒルベルト空間内のレベルの広がり（level broadening）の分布が熱化のダイナミクスをどのように制御するかを明らかにすること。

2. 手法と理論的枠組み

本研究では、小さな「温度計（thermometer）」と大きな「浴（bath）」からなる複合系をモデル化しました。

• モデル設定:

  • 全ハミルトニアン：$\hat{H} = \hat{H}_T + \hat{H}_B + \hat{V}_{int}$
  • 相互作用行列要素 $\hat{V}_{int}$ を、非相互作用基底において独立した確率変数としてモデル化します（有限の分散を仮定しない）。これは ETH を相互作用演算子に適用する考え方です。
  • 初期状態は定義されたエネルギーを持ち、密度行列は $\hat{\rho} = f(\hat{H})$ と書けます。

• 主要な導出プロセス:

  1. 拡散伝播関数（Diffusion Propagator）の導出:
     温度計の縮約密度行列の対角要素のダイナミクスを、グリーン関数の積を用いて記述します。長時間極限において、主要な寄与は異なる半平面にあるグリーン関数の極から来ます。
  2. キャビティ方程式（Cavity Equation）の適用:
     相互作用をランダム行列とみなし、大行列サイズ極限における「キャビティ近似」を用います。これは、対応するホッピングハミルトニアンが木グラフ上の運動を記述する場合に厳密であり、密な行列に対しても適用可能です。これにより、自己エネルギー（self-energy）の完全な分布を導出できます。
  3. 拡散近似の確立:
     グリーン関数の非対角要素を、木グラフ上の経路の和として評価し、ランダム行列要素の独立性を仮定して干渉項を無視します（梯子ダイアグラムの総和に相当）。これにより、実空間の相関関数 $D_{\alpha\beta}$ と、相互作用によって引き起こされるレベルの広がり（level broadening）$\Gamma$ の統計分布を結びつける式が得られます。

• 核心的な結果式:
  温度計の伝播関数は、遷移行列 $\hat{\Pi}$ を用いて以下のように表されます（Eq. 17）：
  $$D_{\alpha\beta} = \left( \hat{1} - \hat{\Pi} \right)^{-1}_{\alpha\beta} B_\beta$$
  ここで、$B_\beta$ は局所状態密度とレベルの広がりに依存する項であり、$\hat{\Pi}$ はレベルの広がり $\Gamma$ の統計に依存する遷移行列です。

3. 主要な貢献と発見

1. 非マルコフ的な熱化の一般化:
   熱化の時間スケールは、典型的なレベルの広がり（typical broadening）の逆数によって決定されます。この理論は、マルコフ近似が成立しない場合（例えば、行列要素の分散が発散するレヴィ分布の場合）でも有効であり、非マルコフ的なグローバル・バランス（global balance）の一般化を提供します。
2. 熱分布の出現:
   十分な広がりを持つ分布条件下では、温度計は熱的平衡状態に達することが示されました。特に、$\kappa \to 0$（長時間極限）における伝播関数の極の存在が確率保存を反映し、その極の位置と性質が最終的な熱分布（ボルツマン分布）を決定します。
3. 数値的検証:
   以下の 3 つの異なるモデルに対して、厳密対角化（Exact Diagonalization）による数値計算と理論予測を比較し、良好な一致を確認しました。
   • レヴィモデル（Levy model）: 相互作用行列要素がパレート分布（べき乗則）に従うモデル。分散が存在しない場合でも理論が成立することを示しました。
   • 全結合横場イジングモデル（TFIM）: カオス的なスピン系。ETH が成り立つ領域での熱化を記述します。
   • Imbrie モデル: 1 次元の相互作用乱れ系。弱い乱れ領域での拡散的な振る舞いを再現しました（強い乱れ領域では局在化により理論の精度が低下することも示唆されています）。

4. 結果の解釈と意義

• ヒルベルト空間拡散としての熱化:
  本研究は、熱化を単なる確率遷移ではなく、ヒルベルト空間内における「拡散過程」として捉え直す枠組みを提供しました。実空間での緩和が遅く見える場合でも（MBL 的振る舞い）、ヒルベルト空間内での拡散が進行し、それが広範な統計分布を通じて熱化を導くメカニズムを解明しました。
• 実験的・理論的インパクト:
  量子プロセッサ実験などで観測される「実空間では凍結しているように見えるが、ヒルベルト空間内では緩和が続く」という現象を、スピンガラス現象論や非マルコフダイナミクスと結びつけて説明できます。
• 今後の展望:
  • 縮約密度行列の非対角要素（位相消滅）への拡張。
  • 輸送観測量への適用による、乱れ相互作用系における遅いモード（slow modes）の特性評価。
  • カオス的な多体ダイナミクスの数値シミュレーションにおける拡散近似の実用的な簡略化。

結論

A.V. Lunkin によるこの論文は、マルコフ近似やフェルミ黄金則に依存しない、ヒルベルト空間内のレベル広がり分布に基づく熱化の微視的理論を確立しました。拡散伝播関数の導出と、レヴィ分布を含む多様なモデルでの数値的検証を通じて、非マルコフ的な系における熱化のメカニズムと時間スケールを統一的に記述する強力な枠組みを提供しています。これは、多体局在（MBL）やスピンガラス系を含む、複雑な量子系のダイナミクス理解における重要な進展です。