Local fermion density in inhomogeneous free-fermion chains: a discrete WKB approach

本論文は、滑らかに変化するホッピングや磁場を持つ自由フェルミオン鎖の局所フェルミオン密度を、離散 WKB 近似を用いて任意の充填率やパラメータに対して解析的に導出する新たな手法を提案し、エンタングルメントエントロピーの抑制メカニズムの理解や従来の場の理論的手法を超えた解析的アプローチへの第一歩を提供するものである。

要約

量子の「住み分け」を解き明かす新しい地図

『不規則な自由フェルミオン鎖における局所フェルミオン密度：離散 WKB アプローチ』の解説

この論文は、「量子の世界で、粒子がどこにどれだけ住んでいるか」を、複雑な地形でも正確に予測する新しい地図（数式）を作ったという研究です。

少し難しい専門用語を、日常の風景に例えてお話ししましょう。

1. 舞台設定：量子の「アパート」と「住人」

まず、この研究の舞台は**「一次元の量子チェーン（鎖）」です。
これを「長いアパート」**と想像してください。

• 部屋（サイト）： 1 つ1 つの部屋が原子や電子の位置です。
• 住人（フェルミオン）： 部屋に住んでいるのは「フェルミオン」という粒子たちです。
• ルール： この粒子たちは「同じ部屋に 2 人は入れない（パウリの排他原理）」という厳しいルールを持っています。

通常、このアパートは「均一な高級マンション」で、廊下（ hopping: ホッピング）の広さや、部屋の磁気（磁場）はどこも同じです。しかし、近年の研究では、「廊下が狭くなったり広くなったりする」「部屋の磁気が場所によって違う」という、ごちゃごちゃした（不均一な）アパートに注目されています。

2. 問題：住人の「密度」がどうなるか？

このごちゃごちゃしたアパートに、ある決まった人数の住人（粒子）を入れたとき、**「どの部屋にどれくらい住人がいるか（密度）」**を計算するのは非常に難しいです。

• 均一な場合： 住人は均等に散らばります。
• 不均一な場合： 不思議なことが起きます。
  • 枯渇（Depletion）： 特定のエリアに**「誰も住んでいない（密度がゼロ）」**場所ができてしまう。
  • 飽和（Saturation）： 逆に、**「部屋が満員（密度が 1）」**になって、それ以上入れない場所ができてしまう。

以前の研究では、この現象を「低密度で磁場がない場合」だけしか説明できませんでした。「なぜ満員になるのか？」「どこで枯渇するのか？」を、どんな条件（どんな人数、どんな磁場）でも説明できる「万能な地図」は存在しませんでした。

3. 解決策：新しい「地形図」の描き方

著者たちは、従来の「フィールド理論（物理学の高度な地図）」ではなく、**「離散 WKB 法」**という新しいアプローチを使いました。

これを**「山登りのガイド」**に例えてみましょう。

• 従来の方法： 山全体を巨大な雲（場）として捉え、大まかな傾向を推測する。しかし、細かい地形（磁場やホッピングの変化）が複雑だと、雲の形が崩れて正確な地図が描けなくなる。
• この論文の方法（WKB 法）： 粒子が「波」のように振る舞うことに注目し、「粒子が登れる山（エネルギーが許される場所）」と「登れない崖（エネルギーが足りない場所）」を、一つ一つの階段（格子点）を数えながら、滑らかに繋いでいく方法です。

彼らは、粒子の「波」が、廊下の広さや部屋の磁気という「地形」に合わせて、どう曲がり、どこで止まり、どこで跳ね返るかを、微細なレベルで計算しました。

4. 発見：驚くほどシンプルな「住み分けの法則」

この新しいアプローチから導き出された結論は、驚くほどシンプルで美しい公式でした。

ある部屋の「住人の密度」は、その場所の**「地形の険しさ（ホッピング）」と「壁の高さ（磁場）」、そして「全体の水位（フェルミエネルギー）」**だけで決まります。

• 水位が低い場所： 水が引いて、**「干潟（枯渇）」**になります（密度 0）。
• 水位が高い場所： 水が溢れて、**「満ちた港（飽和）」**になります（密度 1）。
• 水位がちょうどいい場所： 波が穏やかに揺れて、**「適度な水深」**になります。

この公式は、「低密度・無磁場」だけでなく、「高密度・強磁場」の状況でも、なぜ満員になるのか（飽和）を正確に予測しました。まるで、複雑な地形の川でも、水位さえ分かれば「どこが干上がって、どこが水没するか」が一目でわかるようなものです。

5. 応用：量子もつれ（Entanglement）の謎を解く鍵

この研究の最大の意義は、**「量子もつれ（エンタングルメント）」**という、量子力学の最も不思議な現象の理解に役立つ点です。

• 量子もつれ： 離れた 2 つの粒子が、まるで心電図のようにリンクしている状態。
• 発見： 住人が「枯渇（誰もいない）」や「飽和（満員）」しているエリアでは、量子もつれが抑えられてしまうことが知られています。

この論文で得られた「密度の地図」を使えば、**「どの部屋で、どのくらい粒子がリンクしているか」**を、従来の複雑な計算なしに、シンプルに予測できるようになります。特に、均一な世界ではない（臨界点から離れた）複雑な系でも適用できるため、新しい量子コンピュータやシミュレーターの設計に役立つはずです。

まとめ

この論文は、**「複雑で入り組んだ量子の世界でも、粒子の住み分けは、地形と水位のシンプルな関係で説明できる」**ことを示しました。

• 従来の地図： 限られた条件（低密度・無磁場）しか描けなかった。
• 新しい地図（この論文）： どんな条件（高密度・強磁場・不均一な地形）でも、「枯渇」と「飽和」の境界を正確に描ける。

これは、量子物理学の「地形図」を描く技術において、大きな一歩を踏み出したと言えます。まるで、複雑な街並みでも、天気と人口密度さえ分かれば、どこが混雑し、どこが空いているかを予測できるようになったようなものです。

技術概要

以下は、提示された論文「Local fermion density in inhomogeneous free-fermion chains: a discrete WKB approach（不均一自由フェルミオン鎖における局所フェルミオン密度：離散 WKB 法によるアプローチ）」の技術的な要約です。

1. 研究の背景と課題

• 対象系: 一次元量子スピン系、特に site 依存のホッピング強度 $J_n$ と外部磁場 $B_n$ を持つ不均一な自由フェルミオン（XX モデル）鎖。
• 背景: 不均一な XX 鎖は、光格子中の冷たい原子系や量子シミュレータにおいて重要であり、エンタングルメントエントロピーの振る舞いに関しても注目されている。
• 課題: 以前の研究（Mula ら）では、低充填率かつ磁場ゼロの極限において、局所フェルミオン密度の「枯渇（depletion）」現象を説明する経験的な式が導出された。しかし、任意の充填率、任意のホッピング強度、および外部磁場が存在する一般的な場合において、この現象（枯渇および「飽和（saturation）」）を解析的に記述する手法は確立されていなかった。また、従来の場の理論的手法は、臨界点近傍や低充填率などの特定の極限に依存しており、一般性に限界があった。

2. 提案手法：離散 WKB 近似

著者らは、場の理論に頼らず、離散モデルの固有関数が満たす 3 項漸化式に直接 WKB 型の近似を適用するという新規アプローチを提案した。

• 基本方針:
  • 格子間隔 $a \to 0$ の熱力学的極限（$N \to \infty$）を仮定し、離散的な変数を連続変数 $x$ として扱う。
  • 単粒子固有値問題（行列 $H$ の対角化）を、単粒子波動関数 $\phi_n$ に対する漸化式として記述する。
  • この漸化式に対して、WKB 近似（$\psi(x) = e^{iS(x)/a}$ のような Ansatz）を適用し、連続極限での波動関数の漸近形を導出する。
• 導出プロセス:
  1. 漸化式を連続極限で展開し、位相関数 $S(x)$ の微分方程式を導く。
  2. 波動関数の振幅と位相を決定し、転回点（turning points）$B(x) \pm 2J(x) = \epsilon$ における挙動を解析する。
  3. 波動関数が振動領域（$|B(x) - \epsilon| < 2J(x)$）と指数減衰領域（$|B(x) - \epsilon| > 2J(x)$）でどのように振る舞うかを特定し、接続公式を用いて波動関数の包絡線を求める。
  4. 境界条件（ディリクレ境界条件）から、状態密度 $D(\epsilon)$ とフェルミエネルギー $\epsilon_F$ の関係を導く。

3. 主要な成果と結果

A. 局所フェルミオン密度の閉じた式

本研究の最大の成果は、任意の充填率 $\nu$、任意の site 依存ホッピング $J_n$、および磁場 $B_n$ に対して有効な、局所フェルミオン密度 $\langle c^\dagger_n c_n \rangle$ の解析的な漸近式を導出したことである。

連続極限における密度 $\rho(x, \epsilon_F)$ は以下の式で与えられる（離散サイト $n$ における期待値は $a\rho(na, \epsilon_F)$）：

$$\langle c^\dagger_n c_n \rangle \simeq \begin{cases} 0, & \epsilon_F \le B_n - 2J_n \quad (\text{枯渇領域}) \\ \frac{1}{\pi} \arccos\left( \frac{B_n - \epsilon_F}{2J_n} \right), & B_n - 2J_n < \epsilon_F < B_n + 2J_n \\ 1, & \epsilon_F \ge B_n + 2J_n \quad (\text{飽和領域}) \end{cases}$$

• 物理的意味:
  • 枯渇（Depletion）: フェルミエネルギーが局所的なバンド下端 ($B_n - 2J_n$) よりも低い場合、そのサイトの占有確率は 0 になる。
  • 飽和（Saturation）: フェルミエネルギーが局所的なバンド上限 ($B_n + 2J_n$) よりも高い場合、占有確率は 1 になる。
  • 中間領域: 上記の範囲内では、$\arccos$ 関数に従って滑らかに変化する。

B. 既存研究との整合性と拡張

• 磁場ゼロかつ低充填率の極限では、既存の経験式（Mula ら）と一致する。
• 新規発見: 高充填率や強い磁場下では、枯渇の「鏡像」とも言える**飽和（saturation）**現象が解析的に記述可能であることを示した。これは従来の近似式では再現できなかった。

C. 数値検証

提案された式は、以下の多様な不均一鎖モデルに対して、数値シミュレーションと極めて高い精度で一致することが確認された。

1. 均一鎖（Homogeneous chain）: 定常的な密度分布（Friedel 振動を除く）を正しく再現。
2. Krawtchouk 鎖: 特定の充填率で両端に枯渇または飽和領域が現れる複雑な構造を正確に予測。
3. Rainbow 鎖: 中心から指数関数的に減衰するホッピングを持つ系。枯渇領域の位置と幅を正確に記述。
4. Cosine 鎖および非対称 Cosine 鎖: ホッピングと磁場が周期的・非対称に変化する系。複数の転回点を持つ場合でも、密度プロファイルを高精度で再現。

4. 意義と将来展望

• 理論的枠組みの確立: 場の理論（特に共形場理論）が直接適用できないスケール不変性の欠如する系においても、エンタングルメントや相関関数を解析的に扱うための強力な枠組みを提供した。
• エンタングルメントエントロピーへの応用: 自由フェルミオン系では、相関行列の固有値からエンタングルメントエントロピーが計算できる。本研究で得られた波動関数の包絡線と密度分布を用いることで、任意の充填率・磁場におけるエンタングルメントエントロピーの漸近式を導出する道が開かれた。特に、臨界点から離れた領域でのエンタングルメントの「抑制（suppression）」メカニズムを、局所密度の枯渇・飽和を通じて理解できる。
• 一般性: この離散 WKB 法は、より一般的なフェルミオン鎖や、非局所効果を含む問題への拡張も可能である。

結論

本論文は、不均一な自由フェルミオン鎖の局所密度プロファイルを、任意の充填率と外部場に対して解析的に記述する初めての一般論を提示した。離散 WKB 近似を用いることで、場の理論的な近似に依存せず、枯渇と飽和の両現象を統一的に理解する枠組みを確立し、今後のエンタングルメント研究の基礎となる重要な成果である。