BCS-BEC crossover in trapped one-dimensional Fermi-Hubbard chains: entanglement and correlation signatures from DMRG and effective-pairing theory

本論文は、DMRG シミュレーションと有効対形成理論を組み合わせることで、調和ポテンシャルに閉じ込められた一次元フェルミ・ハバード鎖における BCS-BEC クロスオーバーを特徴づけ、空間閉じ込めが相関パターンをどのように再編成し、条件付き相関関数を通じて識別可能な絶縁領域と超流動領域の共存をもたらすかを明らかにする。

要約

長い細い廊下を想像してください。そこには、小さな目に見えないダンサー（電子）が動き回ろうとしています。完璧で無限に続く廊下では、これらのダンサーは厳格なルールに従います。時には一人で踊り、時にはペアを組んでワルツを踊ります。物理学者はこれを「BCS-BEC クロスオーバー」と呼びます。これは、ダンサーが緩く結びついたパートナー（BCS）から、単一の単位として密着して接着されたもの（BEC）へと変化するスペクトルです。

しかし、現実世界では廊下は無限ではなく、壁を持っています。この論文では、研究者たちは、中央で狭くなり、両端で広がる「曲がった廊下」（調和トラップ）に閉じ込められたこれらのダンサーが何を起こすかを研究しました。この閉じ込めはすべてを変えます。

以下に、彼らの発見を簡単に説明します。

1. 設定：混雑した曲がったダンスフロア

研究者たちは、これらの電子を観察するために超高性能のコンピュータシミュレーション（DMRG と呼ばれる）を使用しました。また、数学に迷い込まないで物理学を理解するために、より単純な「おもちゃのモデル」（有効理論）も構築しました。

• トラップ: 廊下がボウルの形をしていると想像してください。ダンサーたちは自然と最も深い部分（中心）に座ろうとします。
• 相互作用: ダンサーたちは互いを無視したり、互いに押し返したり、互いに強く引き合ったりすることができます。研究者たちはペアがどのように形成されるかを見るために「引力」を最大限に高めました。

2. 2 つの極端なダンス

この論文は、電子が振る舞う 2 つの主要な方法を探索します。

• 「緩いワルツ」（BCS 領域）: 引力が弱いとき、電子はペアを形成しますが、それは部屋を横切って手をつなぐ遠距離のパートナーのようです。彼らは広がっており、ある程度独立して動きます。
• 「接着された双子」（BEC 領域）: 引力が非常に強いとき、電子は非常に強くくっつき、単一の重い物体のように振る舞います。彼らは同じ場所に接着されています。

3. 驚き：「絶縁性のコア」と「超流動の翼」

通常の無限の廊下では、フロア全体が同じように振る舞うはずです。しかし、曲がったトラップのため、この論文はシステムに奇妙な二面性があることを発見しました。

• 中心（絶縁体）: 廊下が混雑するにつれ、真ん中のダンサーたちはあまりにも密集して完全に動きを止めます。彼らは固体のブロックに凍りつきます。研究者たちはこれを絶縁領域と呼びます。誰も動けない交通渋滞のようなものです。
• 端（超流動）: ここが魔法です。中心が凍結しているにもかかわらず、廊下の端にいるダンサーたちは自由に踊り続けます。彼らは「超流動」（摩擦のない流れ）を形成します。
• 結果: あなたはサンドイッチを手に入れます。凍結して動かないコアが、流れ踊る殻に囲まれているのです。この論文はこれを複合 INS+SF 相と呼んでいます。

4. 違いをどう見分けたか

ダンサーたちが「緩いワルツ」をしているのか、「接着された双子」として振る舞っているのか、どうやって見分けるのでしょうか？研究者たちはデータを眺める新しい方法を考案しました。

• 「RMS 距離」（ペアの大きさ）: 彼らは、ペアを組む 2 人のダンサーが通常どれくらい離れているかを測定しました。

  • BCSモードでは、ペアは巨大です（部屋を横切って手をつなぐようなもの）。
  • BECモードでは、ペアは微小です（同じ場所に接着されている）。
  • この距離が引力を強めるにつれてどのように縮小するかを観察することで、彼らは一つのダンススタイルからもう一方への移行を明確に見ることができました。

• 「エンタングルメント」（接続性）: 彼らはまた、廊下の左半分と右半分がどの程度「接続」されているかも調べました。

  • 中心が凍結（絶縁体になる）すると、左右の側面間の接続が突然断ち切られます。橋を切断したようなもので、両側はもう互いに「会話」できなくなります。この突然の断ち切りが、絶縁性のコアが形成される正確なタイミングを教えてくれます。

5. なぜ中心が凍結するのか

なぜ真ん中が詰まってしまうのでしょうか？

• 「実効的」トラップ: 電子が接着されたとき（BEC）、それらは重いボソンのように振る舞います。研究者たちは、この接着されたペアにとってトラップが実効的に強く感じられることを発見しました。まるでボウルが単一のダンサーよりもペアにとって深く、急勾配になるかのようです。
• 反発: ペア同士は互いに引き合っていますが、BEC ペアの「接着」された性質が、隣接するものをわずかに押し返します。これにより、凍結したコアの端近くで密度が上下する奇妙な振動が生じます。

発見のまとめ

この論文は、これらの量子ダンサーを曲がった空間に閉じ込めたときに何が起こるかを示しています。

1. 強い引力が彼らを接着させます（BEC）。
2. 混雑が中心を固体のブロックに凍結させます（絶縁体）。
3. 端は流動的で踊り続けます（超流動）。
4. 「緩いペア」と「接着されたペア」の間の移行は単なる滑らかなスライドではなく、ペアの大きさやシステムの接続性のあり方に明確な指紋を残します。

研究者たちは、これらの異なる振る舞いがどこで起こるかを正確にマッピングし、「これだけの混雑とこれだけの引力があれば、凍った中心に踊る翼が生まれる」と教えてくれる「地図」（相図）を作成しました。彼らは、単純な「おもちゃのモデル」が複雑なコンピュータシミュレーションと完全に一致することを証明し、トラップに押し込められた量子物質がどのように振る舞うかについての統合された図を描き出しました。

技術概要

技術的サマリー：閉じ込められた一次元フェルミ・ハバード鎖における BCS-BEC クロスオーバー

問題提起
バーディーン・クーパー・シュリーファー（BCS）からボース・アインシュタイン凝縮（BEC）へのクロスオーバーは、一様フェルミ・ハバード系において確立された現象であり、相互作用強度の増加に伴い、弱く束縛されたクーパー対から強く束縛された局在対へと進化することを特徴とする。しかし、光格子中の超低温原子を用いた実験的実現では、しばしば調和トラップなどの空間的閉じ込めが伴い、並進対称性が破れる。この不均一性により、対相関長、運動量分布、またはレゲット基準に基づく標準的なプローブが曖昧となり、BCS-BEC クロスオーバーの同定が複雑化する。閉じ込められた一次元（1D）系では、積分可能性の喪失により厳密な解析解が得られず、既存の数値的または平均場アプローチは、強く不均一な多体系における対称領域を分類するための物理的に直感的な基準の提供に苦慮している。

手法
著者らは、ゼロ温度における 1D 調和閉じ込めフェルミ・ハバード鎖の BCS-BEC クロスオーバーを、二重のアプローチを用いて調査する：

1. 数値シミュレーション： 行列積状態（MPS）形式の密度行列繰り込み群（DMRG）法を採用する。スピンバランスの取れた系に対して開放境界条件を用いたシミュレーションを ITensor ライブラリを用いて実施する。精度は結合次元によって制御され、初期値は $\chi_{max}=1500$ に設定され、段階的に増加させる。
2. 有効モデル： 物理的直観と解析的ベンチマークを提供するため、著者らは異なる相互作用領域に対して有効モデルを導出する：
   • 強結合領域（$|U| \gg t$）： 第二摂動論（ロウディン摂動論）を用いて、系を（強く束縛されたオンサイト対を表す）ハードコア・ボソンによる閉じ込め tight-binding モデルに写像する。このモデルには、再正規化されたホッピング（$\tau \approx 2t^2/|U|$）、増強された有効閉じ込め強度（$\tilde{k} = k|U|/2t$）、および最近接反発が含まれる。
   • 弱から中程度の結合領域（$|U| \lesssim 2t$）： 相互作用項を結合解除してボゴリューボフ準粒子を自己無撞着に解く平均場（MF）近似を用いる。
   • 非相互作用極限（$U=0$）： 閉じ込め tight-binding モデルの解析解を用いて、スペクトルおよび密度プロファイルの参照点を確立する。

主要な貢献と結果

• 相図と転移： 本研究は、相互作用強度 $|U|$ と平均密度 $n$ の関数として相図を構築する。3 つの主要な領域を同定する：

  • BEC 領域： 強く束縛された対が支配的であり、低密度または非常に強い相互作用で生じる。
  • BCS 領域： 広がったクーパー対が特徴であり、中程度の密度と弱い相互作用で生じる。
  • 複合 INS+SF 相： 高密度において、トラップ中心は局在状態の形成により絶縁体（INS）となる一方、超流動（SF）状態は端部に存続する。この「超流動の翼」と呼ばれる現象は、相互作用強度にわたって頑健である。

• 新しい診断ツール：

  • 条件付き相関関数： 著者らは、反対スピン間の局所的な束縛を調べる 4 粒子相関関数 $P_j(r) = \langle c^\dagger_{j+r,\uparrow} c_{j,\uparrow} n_{j,\downarrow} \rangle$ を導入する。この分布の二乗平均平方根距離（$r_{rms}$）が、BCS 領域では緩やかな減衰（$\nu_1 \sim 0.1$）、BEC 領域では急激な減衰（$\nu_2 > 1$）という 2 つの異なるべき乗則スケーリングを示すことを実証する。クロスオーバーは、これらのスケーリング挙動間の遷移によって同定される。
  • エンタングルメントに基づくプローブ： 半鎖エンタングルメントエントロピー（$S_{L/2}$）は二重の役割を果たすことが示される。$S_{L/2}$ の急激な低下は、超流動から絶縁体への転移（中央絶縁障壁の形成）をシグナルする。さらに、$S_{L/2}$ の大きさを参照値（$2 \ln 2$）と比較することで、BCS 様（弱相関、$S_{L/2} \ge 2 \ln 2$）と BEC 様（強く束縛された対、$S_{L/2} < 2 \ln 2$）の領域を区別できる。

• 物理的メカニズム：

  • 端部超流動： 有効モデルは、強結合極限において、トラップによる並進対称性の破れに起因する BEC 様準粒子間の有効反発が、対をトラップ中心から押しやることを明らかにする。このメカニズムは、中心が絶縁体であっても端部に超流動状態が存続する理由を説明する。
  • 閉じ込め増強： 本研究は、引力相互作用が実効的にトラップポテンシャルを強化（$\tilde{k} > k$）し、中心密度を高め、$|U|$ の増加に伴って超流動 - 絶縁体転移をより低い密度へシフトさせることを確認する。
  • 特異な BCS 対： 中程度の領域において、平均場解析は、BCS 様対がトラップ中心を避け、空サイトと中央領域の界面に局在することを明らかにする。これは DMRG 結果で観測された密度振動と整合的である。

意義
本論文は、堅牢な数値結果と有効理論的記述を組み合わせることで、調和閉じ込め幾何学における BCS-BEC クロスオーバーの統一された図像を提供すると主張する。その主な意義は以下の点にある：

1. 同定課題の解決： 標準的なプローブが機能しない不均一系において、BCS 領域と BEC 領域を区別するための、相関関数のべき乗則スケーリングおよびエンタングルメントエントロピーに基づく、新しい物理的に透明な基準を提供する。
2. 端部状態の理解： 閉じ込め鎖の端部における超流動状態の頑健性に対する微視的説明を提供し、それらを対間の有効反発および閉じ込めポテンシャルにおける対の特定の空間構造に帰着させる。
3. 整合性： DMRG シミュレーションと有効モデル（平均場およびハードコア・ボソン記述を含む）との間に、定性的および定量的な整合性が存在することを示し、強結合極限における本質的な物理を捉えるためにこれらの簡略化されたモデルを使用することを検証する。
4. 実験的関連性： 発見事項は、準 1D 調和トラップにおける局所占有確率の既存の実験的観測、および最近の 2D 実験、特に非対原子の抑制と電荷密度波様変調の出現に関して、定性的に整合すると提示される。

著者らは、同定された特徴は主に調和閉じ込めの特定の形式ではなく、並進対称性の破れに由来することを結論付け、定性的な発見はより一般的なトラップポテンシャルにおいても存続する可能性を示唆している。