Suppressed correlation-spreading in a one-dimensional Bose-Hubbard model with strong interactions

強い相互作用を持つ一次元ボース・ハバード模型において、初期状態が二重占有密度波である場合、ドメインウォール励起による相関の伝播が抑制され、特に放物型トラップが存在するとその伝播がさらに凍結されることを示し、この現象が強相互作用極限における反強磁性トランスバース場イジング模型のスピン波励起の群速度によって記述可能であることを明らかにした。

要約

🧊 物語の舞台：「ボース・ハバード模型」というお菓子屋

まず、実験の舞台を想像してください。
長い通路（一次元）に、**「お菓子（原子）」**が置かれています。
この通路には、お菓子が隣同士に飛び移れる「トンネル（ジャンプ）」と、同じ場所に 2 つのお菓子が集まると激しく反発し合う「壁（強い相互作用）」があります。

研究者たちは、最初にお菓子を**「2 つ、0 個、2 つ、0 個……」**という規則正しい並び（密度波）に整えました。

• 2 つ ＝ 2 個のお菓子
• 0 個 ＝ 空っぽ

この状態から、お菓子たちがどう動き出すか（時間経過）を観察したのがこの研究です。

---

🚦 発見その 1：「ドッペルゲンガー」の交換ゲーム

通常、お菓子（原子）は 1 つずつ隣に飛び移ろうとします。しかし、「壁（相互作用）」が非常に強い場合、1 つずつ動くのはエネルギー的に大変で、ほとんど動きません。

代わりに、面白いことが起きます。
「2 つのお菓子（ドゥーブロ）」と「空っぽの場所（ホロン）」が、隣同士でペアになって入れ替わるのです。

• 元の状態： [2 個][0 個]
• 入れ替わり： [0 個][2 個]

これを**「ドゥーブロ・ホロン交換」と呼びますが、イメージとしては「雪だるまの列」で考えるとわかりやすいです。
雪だるまが 2 人組で並んでいる場所と、誰もいない場所が隣り合っていると、「2 人組」と「空き地」が、まるでダンスのように隣同士で位置を交換していく**のです。

この「交換の波」が、通路の端から端へと伝わっていきます。これを論文では**「ドメインウォール（境界壁）の移動」**と呼んでいます。

---

🧱 発見その 2：「壁」があると動きが止まる（非エルゴード性）

ここが今回の最大の発見です。
もし、この通路の両端に**「重力（トラップポテンシャル）」**がかかっていたらどうなるでしょうか？

• 重力がない場合： 「2 人組と空き地の交換」が端から端まで行き渡り、最終的にはお菓子が均一に散らばってしまいます（熱平衡）。
• 重力がある場合： 通路の端（エッジ）にいるお菓子たちが、まるで凍りついたように動けなくなります。

面白いことに、この重力は**「トンネルを渡る力」よりもずっと弱いにもかかわらず、端の動きを完全に止めてしまいます。
その結果、お菓子たちの「2 人組と空き地の交換」が、通路の中央付近でしか起こらなくなり、「情報の伝播（広がり）」が抑えられてしまいます。**

これを**「非エルゴード性（非平衡状態）」と呼びます。
つまり、「時間がいくら経っても、お菓子たちは元の整然とした並びを完全に忘れられず、動きが凍りついたままになる」**という現象です。

---

🔍 発見その 3：「相関」という目に見えない糸

研究者たちは、お菓子同士がどう「連動」しているか（相関）も調べました。

1. 単独のお菓子（単一粒子相関）： ほとんど動かない。隣のお菓子ともほとんど繋がっていない。
2. ペアのお菓子（ペア相関）： 少し動くが、やはり隣までしか広がらない。
3. 密度の波（密度 - 密度相関）： これが重要！「2 人組と空き地の交換」が進むと、「2 人組→空き地→2 人組」という波が、反対側の端まで伝わっていく様子が見えました。

しかし、重力（トラップ）をかけると、この「波」の広がりも抑えられ、**「遠くまで情報が届かない」**状態になります。

---

🧩 発見その 4：「魔法の鏡」で単純化

この複雑な原子の動きを理解するために、研究者たちは**「鏡」を使いました。
原子の「2 個と 0 個」の並びを、「上向き（↑）」と「下向き（↓）」の磁石（スピン）**に置き換えるのです。

• [2 個][0 個] → ↑
• [0 個][2 個] → ↓

すると、この複雑な原子の世界は、**「横から風が吹く磁石の模型（横磁場イジング模型）」**という、物理学でよく知られた単純なモデルに変換できました。

この「鏡」を通して見ると、「2 人組と空き地の交換」は、単に「磁石の向きが反転する（スピン・フラップ）」ことであることがわかりました。
さらに、この単純なモデルを使うと、**「波が伝わる速さ」**を正確に計算できることが証明されました。

---

📝 まとめ：何がすごいのか？

この論文は、**「強い力で結ばれた原子の列」において、「弱い重力さえあれば、情報の伝播が凍りつく」**という現象を解明しました。

• 日常の例え：
  想像してください。長い廊下に並んだ人々が、隣の人と手を取り合って「右・左・右・左」とリズムよく移動しようとしています（これが原子の動き）。
  しかし、廊下の端に**「軽い風（弱い重力）」が吹いているだけで、端にいる人々が「足が止まって動けなくなる」**。
  その結果、リズムの波が廊下の中央で止まり、廊下全体が「動かない状態」のまま固定されてしまいます。

この研究は、**「量子コンピュータ」や「新しい物質の設計」**において、情報をどう制御し、どう保存するか（あるいは壊さないようにするか）を考える上で、非常に重要なヒントを与えてくれます。

「強い相互作用」と「弱い重力」の組み合わせが、原子の世界に「凍りついた秩序」を作り出すという、一見矛盾しているような不思議な現象を、数式とシミュレーションで見事に描き出したのです。

技術概要

以下は、提示された論文「Suppressed correlation-spreading in a one-dimensional Bose-Hubbard model with strong interactions（強い相互作用を持つ 1 次元ボース・ハバードモデルにおける相関の伝播の抑制）」の技術的な要約です。

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

• 背景: 孤立した多体系量子系において、熱化（thermalization）とエルゴード性の破れは統計力学の基礎的な課題です。通常、固有状態熱化仮説（ETH）が満たされれば系は熱化しますが、積分可能系や多体局在（MBL）系、ヒルベルト空間の断片化（Hilbert Space Fragmentation: HSF）が生じる系などでは、ETH が破れ、非エルゴード的な振る舞いが観測されます。
• 問題設定: 従来の研究では、1 次元ボース・ハバードモデルにおいて「単一占有密度波状態」から始まるダイナミクスが速やかに熱化することが示されていました。一方、「二重占有密度波状態（$|2, 0, 2, 0, \dots\rangle$）」から始まる場合、強い相互作用と弱いトラップポテンシャルの組み合わせにより非エルゴード的になることが知られていましたが、その相関関数の時間発展における具体的な振る舞い、特に「相関の伝播（correlation spreading）」がどのように抑制されるかについては未解明な部分がありました。
• 目的: 二重占有密度波状態からのクエench（急激なパラメータ変化）ダイナミクスにおいて、強い相互作用領域での非エルゴード性の兆候、特に相関関数の空間的広がり（伝播）がどのように抑制されるかを解明すること。

2. 手法

• モデル: 1 次元ボース・ハバードモデル（開境界条件）に、放物線型のトラップポテンシャルを加えた系を考慮。
  • ハミルトニアン：ホッピング項（$J$）、オンサイト相互作用（$U$）、トラップポテンシャル（$\Omega$）。
  • 初期状態：二重占有密度波状態 $|\psi_0\rangle = |2, 0, 2, 0, \dots\rangle$。
• 数値計算:
  • 行列積状態（MPS）法を用いた厳密な時間発展シミュレーション。
  • アルゴリズム：時間発展ブロック消去法（TEBD）。
  • パラメータ設定：相互作用 $U/J = 40$（強相互作用領域）、格子サイズ $M=16$（または $12$）、最大結合次元$\chi_{max}=1600$、サイトあたりの最大占有数$n_{max}=4$。
  • トラップ強度：$\Omega/J = 0$（なし）および $0.02$（弱いトラップ）。
• 解析対象:
  • サイト占有数 $n_i(t)$ の時間発展。
  • 相関関数：単一粒子相関 $C_{sp}$、対相関 $C_{pair}$、密度 - 密度相関 $C_{dd}$。
  • 相関の広がり幅 $\sigma(t)$ の評価。
• 理論的アプローチ: 強相互作用極限（$U \gg J$）において、有効ハミルトニアンを導出し、系を横磁場イジングモデル（Transverse-Field Ising Model: TFI）へ写像することで、ダイナミクスを直感的に理解。

3. 主要な結果

A. サイト占有数のダイナミクス

• トラップなし（$\Omega=0$）の場合:
  • 初期状態は高エネルギー状態であり、単一粒子のホッピングはエネルギー的に抑制される（$U$ に比例するエネルギー減少が必要）。
  • 代わりに、ドゥブロン（二重占有）とホロン（空サイト）の交換（ doublon-holon exchange）、すなわち対トンネリング（pair hopping）が支配的となる。
  • このプロセスは、空の端から「ドメインウォール（$|2,0\rangle$ と $|0,2\rangle$ の境界）」を励起し、反対側へ伝播させる現象として記述される。
  • 伝播速度は有効ホッピング強度 $\tilde{J} = 2J^2/U$ に依存し、時間スケールは $U/J$ の逆数に比例して遅くなる。
• トラップあり（$\Omega > 0$）の場合:
  • 端のサイトがエネルギー的に分離され、局在化する。
  • 特に左端ではドゥブロン - ホロン交換のエネルギーコストが最大となり、占有数が固定される。
  • 結果として、ドメインウォールの伝播が抑制され、系の中心付近でのみ動的変化が生じる。

B. 相関関数の振る舞いと伝播抑制

• 単一粒子・対相関関数（$C_{sp}, C_{pair}$）:
  • どちらのケースでも、相関は隣接サイト付近に強く局在しており、バルクへは急速に減衰する。これは対トンネリングが支配的であることと整合する。
  • トラップが存在すると、相関の全体的な強度が低下し、伝播距離がさらに制限される。
• 密度 - 密度相関関数（$C_{dd}$）:
  • 単一粒子相関とは異なり、長距離でも無視できない相関が維持される。
  • 時間発展に伴い、初期の密度波パターンとは逆転した「負の階段状パターン（negated staggered pattern）」が現れる。これはドメインウォールの励起と伝播の直接的な証拠である。
  • 相関の広がり幅 $\sigma(t)$ の解析:
    • トラップなし：$U/J$ が増加すると、$\sigma(t)$ の最初のピークがより遅い時間へシフトする（熱化の遅延）。しかし、十分長い時間では緩和（熱化）が見られる。
    • トラップあり：$U/J$ が増加しても $\sigma(t)$ は減少し続け、長時間平均しても飽和しない。相関の伝播が強く抑制され、非エルゴード的な状態が維持される。
    • この抑制は、トラップ強度がトンネリング率 $J$ よりも弱い場合（$\Omega/J = 0.02 \ll 1$）であっても発生する。

C. 有効モデル（横磁場イジングモデル）による記述

• 強相互作用極限において、ボース・ハバードモデルは有効ハミルトニアン（Eq. 4）へ、さらにスピン写像（$|2,0\rangle \to |\uparrow\rangle$, $|0,2\rangle \to |\downarrow\rangle$）により横磁場イジングモデル（TFI）（Eq. 5）へ写像できる。
• この写像により、ドメインウォールの励起は「スピン反転励起（spin-flip excitation）」として解釈できる。
• TFI モデルの励起スペクトルから計算される群速度 $v_g$ が、元のボース・ハバードモデルにおける相関伝播速度（$2v_g$）を非常に良く再現していることを確認した。
• この有効モデルは、格子サイズを半分に縮小できるため、より大きな系での計算や直感的な理解を可能にする。

4. 結論と意義

• 結論:
  1. 二重占有密度波状態からのダイナミクスにおいて、強い相互作用領域では対トンネリング（ドゥブロン - ホロン交換）が支配的であり、これがドメインウォールの伝播を引き起こす。
  2. 単一粒子・対相関は強く局在するが、密度 - 密度相関はドメインウォールの伝播を示す。
  3. 最も重要な発見として、非常に弱いトラップポテンシャル（$\Omega \ll J$）が存在するだけで、相関の伝播が著しく抑制され、非エルゴード的な振る舞いが長時間維持されることを示した。
  4. 強相互作用極限では、系が横磁場イジングモデルへ写像可能であり、そのスピン波の群速度が元の系の相関伝播速度を記述できる。
• 意義:
  • 従来の MBL や積分可能系とは異なるメカニズム（ヒルベルト空間の断片化とドメインウォールダイナミクス）による非エルゴード性の明確な実証。
  • 弱いトラップであっても、強い相互作用と組み合わせることで相関伝播を制御・抑制できる可能性を示唆し、量子情報伝播の制御や非平衡量子多体系の理解に寄与する。
  • 複雑なボース・ハバードモデルのダイナミクスを、より単純なスピンモデル（TFI）で捉える有効なアプローチを提供した。

この研究は、冷原子実験におけるパラメータ制御（特に相互作用強度とトラップポテンシャルの調整）を通じて、非エルゴード相を安定的に実現・観測するための指針となるものです。