Eigenstate Thermalization Hypothesis (ETH) for off-diagonal matrix elements in integrable spin chains

本研究では、代数ベア Ansatz を用いて XXX 鎖の固有状態間の局所演算子の非対角行列要素を解析し、熱的マクロ状態内の要素が指数関数的に減衰し、異なるマクロ状態間の要素はより急速に減衰することを示し、両者の確率分布がグンベル分布で記述されることを明らかにしました。

要約

1. 舞台設定：カオスなパーティー vs 整然とした行進

まず、この研究が扱っている「量子システム」を二つのパーティーに例えてみましょう。

• カオスなパーティー（非積分可能系）：
  参加者が自由に動き回り、誰とでも話しかけ、予測不能な騒ぎになるようなパーティーです。この世界では、**「固有状態熱化仮説（ETH）」**というルールが成り立ちます。これは、「どんな部屋（状態）にいても、最終的には全員が同じような平均的な振る舞いをする（熱化する）」というルールです。
• 整然とした行進（積分可能系・この論文のテーマ）：
  ここは、全員が厳格なルールに従って、整然と行進しているような世界です（例：スピン1/2のXXX鎖）。通常、このような「整然とした」世界では、カオスな世界のような「熱化」は起きないはずだと考えられていました。

この論文の目的：
「整然とした行進（積分可能系）の中でも、実は『熱化』に近い現象が起きているのか？そして、その仕組みはカオスな世界とどう違うのか？」を、**「行列要素（Matrix Elements）」**という数値を調べることで突き止めようとしています。

2. 鍵となる道具：「行列要素」とは何か？

論文で扱っている「行列要素」とは、簡単に言うと**「ある状態から、別の状態へ『ジャンプ』する時の確率の大きさ」**です。

• 対角成分（同じ状態から同じ状態へ）：
  「その状態が、どれくらいエネルギーを持っているか」を表す値。
• 非対角成分（異なる状態間でのジャンプ）：
  「ある状態から、別の状態へ飛び移る時の『しやすさ』」を表す値。

この「ジャンプのしやすさ」が、システムがどのように熱平衡に達するかを決定づけます。

3. 発見された驚きの事実

研究者たちは、この「整然とした行進（XXX鎖）」の中で、**「ジャンプのしやすさ」**を巨大な計算機を使ってシミュレーションしました。その結果、以下のような面白いことがわかりました。

A. 同じグループ内でのジャンプ（同じ「マクロ状態」同士）

同じ温度やエネルギー密度を持つ状態同士でのジャンプを調べました。

• 発見： ジャンプの確率は、システムが大きくなる（人数が増える）につれて、**「指数関数的に急激に小さくなる」**ことがわかりました。
• アナロジー： 大きな広場で、同じグループの人同士が手をつなごうとしても、人数が増えれば増えるほど、偶然手が触れる確率は極端に低くなります。
• 重要点： この「小さくなる速さ」は、カオスな世界（ETH）とほぼ同じでした。つまり、「整然とした行進」でも、局所的には「熱化」の兆候が見られるのです。

B. 異なるグループ間でのジャンプ（異なる「マクロ状態」同士）

例えば、「絶対零度（氷のような状態）」と「無限高温（沸騰した状態）」のような、全く異なる性質を持つ状態同士でのジャンプを調べました。

• 発見： こちらは、同じグループ内よりもさらに激しく、爆発的に小さくなりました。
• アナロジー： 氷の塊と沸騰したお湯が、お互いの分子を交換しようとしても、それはほぼ不可能に近いほど確率が低いです。

4. 最大の驚き：分布の「形」が違う！

ここがこの論文の最も面白い部分です。ジャンプの確率の「大きさ」をグラフにすると、その**「形（分布）」**に大きな違いが見つかりました。

• カオスな世界（ETH）の予想：
  確率の分布は「ガウス分布（ベル型の鐘の曲線）」になるはず。これは、平均的な値の周りに、均等にばらつくイメージです。
• この論文の結果（整然とした世界）：
  分布の形は**「ギューベル分布（Gumbel distribution）」という、「右に長い尾を持つ、偏った形」**でした。
  • アナロジー：
    • ガウス分布（カオス）： 多くの人が平均的な身長で、背が高い人・低い人が均等に散らばっている。
    • ギューベル分布（積分可能）： 多くの人が平均的な身長だが、**「とんでもなく背が高い巨人」**がごく少数、ひょっこりと現れる形。

この「偏った分布（ギューベル分布）」が見つかったことは、**「積分可能系では、カオスな世界とは異なる、独特な『熱化』のメカニズムが働いている」**ことを示しています。

5. なぜこれが重要なのか？

• 「ストリング（String）」の存在：
  積分可能系では、粒子が束になって「ストリング（紐）」のような状態を作ることがあります。これまでは、このストリングがあるせいで計算が難しすぎて、正確な答えが出せませんでした。
• 技術の進歩：
  この論文では、最新の数学的手法（代数的ベア Ansatz）を駆使して、この「ストリング」を含む複雑な計算を、従来の方法（対角化）では不可能だった巨大なサイズ（L=240 程度）まで実行することに成功しました。
• 結論：
  「ストリング」があっても、熱化の基本的なスケール（指数関数的な減衰）は変わらないが、**「確率の揺らぎ（分布の形）」**は、カオスな世界とは全く異なる独特のものになることが証明されました。

まとめ

この論文は、「整然とした量子世界（積分可能系）」においても、「熱平衡」という現象が、「カオスな世界」とは異なる、しかし驚くほど似た形で現れていることを明らかにしました。

• 同じ： 状態が変わる確率は、系が大きくなると急激に消え去る。
• 違う： その消え方の「揺らぎの形」は、カオスな世界の「平均的な鐘の曲線」ではなく、**「偏ったギューベル分布」**という、より劇的な形をしている。

これは、量子力学の「熱化」という現象が、単一のルールではなく、システムの性質（カオスか整然か）によって、**「同じような結果（熱化）に至るが、その道筋（統計的な振る舞い）は多様である」**ことを示唆しています。まるで、同じ目的地（熱平衡）に向かうのに、カオスな世界は「大勢で均一に進む」のに対し、整然とした世界は「少数の巨人が先導する」ような、異なる旅のスタイルを持っているようなものです。

技術概要

論文タイトル

積分可能スピン鎖における非対角行列要素の固有状態熱化仮説 (ETH)
（原題：Eigenstate Thermalization Hypothesis (ETH) for off-diagonal matrix elements in integrable spin chains）

1. 研究の背景と問題設定

孤立した量子多体系における非平衡ダイナミクスと熱化の理解は、現代物理学の重要な課題です。

• ETH の概要: 一般の（カオス的な）非積分系では、固有状態熱化仮説（ETH）が成立し、局所演算子の行列要素は特定の統計的構造（対角成分はエネルギーの滑らかな関数、非対角成分は指数関数的に抑制され、ランダム行列理論に基づくガウス分布に従う）を持つことが知られています。
• 積分可能系における課題: 積分可能系（例：XXX 鎖）では、無数の保存量が存在するため、標準的な ETH は破綻します。しかし、対角行列要素の振る舞いについては、代数ベテ Ansatz (ABA) を用いた研究が進んでいますが、非対角行列要素の統計的性質、特に積分可能系特有の「束縛状態（ストリング）」がその振る舞いにどう影響するかは未解明でした。
• 先行研究との対比: 最近、Lieb-Liniger モデル（連続体モデル）において、非対角行列要素の有限サイズスケーリングと確率分布が研究されました [26]。そこでは、同じ熱力学的マクロ状態内の状態間では $L \ln L$ の補正を伴う指数減衰、異なるマクロ状態間では $L^2$ の減衰が示されました。しかし、この結果が格子モデル（スピン鎖）、特に束縛状態（ストリング）を含む系に適用可能かどうかは不明でした。

2. 研究方法

本研究では、**等方性スピン 1/2 ハイゼンベルグ鎖（XXX 鎖）**を対象とし、以下の手法を用いて非対角行列要素を高精度に計算しました。

• 代数ベテ Ansatz (ABA) の活用:
  • 厳密対角化法では到達不可能な大規模系（$L \sim 500$）を扱うため、ABA に基づく局所演算子の行列要素の厳密な公式 [30, 31] を採用しました。
  • 計算コストは演算子のサポートサイズ $\ell$ に対して指数関数的ですが、粒子数 $M$ に対して多項式的（$O(M^\ell)$）にスケールするため、大規模計算が可能です。
• 対象とする演算子:
  • 1 スピン演算子: $E^{(11)}_i$（スピン 1 個に作用）。
  • 2 スピン演算子: $E^{(11)}_i E^{(22)}_{i+1}$（隣接する 2 スピンに作用）。
• 特異点の除去とストリング仮説:
  • 一様極限（homogeneous limit）やストリング偏差（string deviations）の極限を取る際に生じる「虚特異点」を慎重に除去する手法を適用しました。
  • 1 スピン演算子については、ストリングを含む状態でも解析的に極限を処理可能ですが、2 スピン演算子については計算の複雑さから、**ストリングを含まない固有状態（束縛状態を持たない状態）**に限定して数値計算を行いました。
• サンプリング手法:
  • 熱力学的マクロ状態（ギブスアンサンブル）から固有状態をサンプリングするために、メトロポリス法を用いたヒルベルト空間モンテカルロ法を採用しました。
  • 対象としたマクロ状態：無限温度（$\beta=0$）、有限温度（$\beta=0.5$）、およびゼロ温度（基底状態）。

3. 主要な結果

A. 同じマクロ状態内の固有状態間の非対角行列要素

• スケーリング則:
  • 行列要素の対数 $M_{ij} = \ln |\langle E_i | O | E_j \rangle|^2$ は、系サイズ $L$ に対して $M_{ij} \sim -L$ のように指数関数的に減衰します。
  • 詳細な解析により、$M_{ij} \approx -c L \ln L - c_0 L$ というスケーリングが確認されました。これは Lieb-Liniger モデルの結果 [26] と定性的に一致しており、束縛状態（ストリング）の存在がスケーリングの定性を変化させないことを示しました。
• 確率分布関数 (PDF):
  • 対数行列要素 $M_{ij}$ の分布は、グンベル分布 (Gumbel distribution) で非常に良く記述されます。
  • これは、非積分系における ETH が予測する「ガウス分布」とは明確に異なります。
  • 異なるマクロ状態（異なる温度や粒子密度）や異なる演算子（1 スピン vs 2 スピン）においても、グンベル分布が成立することが確認されました。

B. 異なるマクロ状態間の固有状態間の非対角行列要素

• スケーリング則:
  • 異なるマクロ状態（例：無限温度と基底状態）に属する固有状態間の行列要素は、同じマクロ状態内の場合よりも急速に減衰します。
  • 具体的には、$|\langle E_i | O | E_j \rangle| \propto \exp(-d L^2)$ のように、$L^2$ に比例する指数減衰を示しました。
• 分布:
  • この場合も、対数行列要素の分布はグンベル分布で記述されることが確認されました。

C. ガウス性からの逸脱の検証

• 比率 $\Gamma$ の解析:
  • 行列要素がガウス分布に従う場合、比率 $\Gamma = \overline{|O_{ij}|^2} / (\overline{|O_{ij}|})^2$ は $\pi/2$ に収束するはずです。
  • 本研究の数値計算では、$\Gamma$ が $\pi/2$ から大きく逸脱し、系サイズ $L$ の増加とともに発散する傾向（$\Gamma \sim 10^5$ 程度）が観測されました。
  • これは、積分可能系における非対角行列要素の統計が非ガウス的であることを強く示唆しており、その分布がグンベル分布であることと整合します。

4. 貢献と意義

1. 積分可能格子モデルへの ETH スケーリングの適用:
   • 連続体モデル（Lieb-Liniger）で示された非対角行列要素のスケーリング則（$L \ln L$ 補正や $L^2$ 減衰）が、束縛状態を持つ**格子モデル（XXX 鎖）**でも成立することを初めて実証しました。
2. 統計分布の明確化:
   • 積分可能系における非対角行列要素の統計が「ガウス分布」ではなく「グンベル分布」に従うことを、大規模数値計算によって確立しました。これは、積分可能系と非積分系の決定的な違いの一つとして位置づけられます。
3. 技術的ブレイクスルー:
   • ABA 公式を用いて、ストリングを含む状態（1 スピン演算子）やストリングを含まない状態（2 スピン演算子）において、$L \sim 500$ までの大規模系での行列要素計算を可能にしました。これは従来の厳密対角化法（通常 $L \lesssim 20$）を遥かに凌駕する規模です。
4. 将来の展望:
   • この結果は、自由確率論（Free Probability Theory）との関係性や、量子クエンチ後の緩和ダイナミクス、および初期状態と固有状態の重なり（overlaps）の統計的理解への新たな道を開いています。

結論

本論文は、積分可能スピン鎖における非対角行列要素の振る舞いを、大規模な代数ベテ Ansatz 計算によって詳細に解明しました。その結果、積分可能系においても ETH 的な指数減衰スケーリングは保たれるものの、その統計的性質（グンベル分布への従属）は非積分系（ガウス分布）とは根本的に異なることが示されました。これは、積分可能系の熱化メカニズムと非積分系との本質的な違いを浮き彫りにする重要な成果です。