Scaling of entanglement with distance between subsystems: detection of phase transitions in the $t'$-$\pi$-flux model

技術要約：相転移検出のための距離に伴う絡み合いのスケーリング

問題提起
強相関フェルミオン系における量子相転移の検出および分類は、標準的なランダウ＝ギンツブルグ秩序パラメータが存在しない、あるいは定義が困難な場合、依然として大きな課題となっている。エンタングルメント・エントロピーのシステムサイズ・スケーリングのような、エンタングルメントに基づくアプローチは理論モデルに対しては有効であることが証明されているが、相互作用するフェルミオン系においては深刻な制限に直面している。具体的には、厳密解は稀であり、量子モンテカルロ（QMC）法のような数値的手法は、しばしば小さなシステムサイズに制限される。さらに、広範な部分系の簡約密度行列（RDM）を構築し、保存することは計算量的に極めて困難である。近接する小さな部分系間のエンタングルメントは、特定のケース（例：スピン系、キタエフ・モデル）において転移の兆候を示すことが示されているが、これらの兆候が、特に強い相互作用や有限温度が存在するあらゆるシナリオにおいて、必ずしも持続したり検出可能であったりするとは限らない。

手法
著者らは、全システムサイズ $L$ を変化させる代わりに、2つの小さな（単一サイトまたはダイマーである）非連結な部分系間の距離 $d$ を変化させることで相転移を検出する「距離によるスケーリング」アプローチを提案している。

1. モデル系: 本研究では、正方格子上の $t'$-$\pi$-flux ハバードモデルに焦点を当てる。このモデルは、以下の3つの異なる転移を示す豊かな相図を持つことが選定理由である：

   • 半金属（SM）からバンド絶縁体（BI）（ガウス普遍性クラス）。
   • 半金属（SM）から反強磁性（AFM）モット絶縁体（MI）（グロス・ニューウェン・普遍性クラス）。
   • 磁気秩序・無秩序転移（MI から BI）（O(3) ウィルソン・フィッシャー普遍性クラス）。

2. エンタングルメント尺度:

   • フェルミオン的ネガティビティ ($N^F$): 二部グラフのエンタングルメントの主要な尺度として使用される。これは、部分時間反転操作を介して正の部分転置（PPT）基準をフェルミオン・モードへと拡張したものであり、局所的なパリティ破れが存在する場合でもエンタングルメントを確実に捉える。著者らは、対数ネガティビティ $\mathcal{N}^F = \log_2(2N^F + 1)$ を用いる。
   • スピン・エンタングルメント ($N^S$): 局所的な電荷を保存しないRDMの非対角要素を除去した、フィルタリングされたバージョンのネガティビティである。これにより、純粋にスピン揺らぎに由来するエンタングルメントを孤立させる。
   • エンタングルメント・エントロピー: 非相互作用極限（$U=0$）において、SM-BI転移の参照用として計算される。

3. 計算手法:

   • 補助場QMC法: シミュレーションは、有限温度（$T = 1/16t$）においてALFコード・フレームワークを用いて行われる。
   • RDM再構成: 著者らは、QMCデータから任意の等時間および時間変位相関関数を計算する一般的な手法を導入している。これには、フォック空間ベクトルに作用する線形演算子のためのコンピュータ代数システムが含まれ、観測量をフェルミオン・グリーン関数の積へと変換する。
   • ノイズ緩和: 大きな距離におけるQMCノイズに対処するため、Vidalらによる公式を拡張することで、ノイズを含むデータからネガティビティを抽出する方法を導出している。また、各観測量に対して揺らぎの小さいハバード・ストラトノビッチ・チャネル（スピンまたは電荷）を選択する戦略も採用している。
   • 解析: ネガティビティのチューニング・パラメータ（$t'$ または $U$）に対する一次および二次の対数微分が、主要な診断ツールとなる。

主な貢献

• 距離スケーリング・プロトコル: 有限サイズ・スケーリング解析におけるシステムサイズ $L$ と同様に、2つの小さな部分系間の距離 $d$ が役割を果たすプロトコルを確立した。これにより、大きなシステムサイズのスケーリングが不可能な領域においても転移の検出が可能となる。
• アルゴリズムの進展: グリーン関数からマルチサイトの簡約密度行列を計算するための、QMCコード内での自動化されたウィック分解の開発。これにより、任意の部分サイズに対するレニー・エントロピーやネガティビティを超えたエンタングルメント尺度の計算が可能となった。
• ノイズ処理: ネガティビティを計算する際のQMCノイズを軽減するための具体的な技術を提供し、距離スケーリング解析の信頼性を確保した。

結果
本研究は、フェルミオン的ネガティビティの距離に関する対数微分に、相転移の鋭い兆候が現れることを示している：

1. SM-BI 転移（非相互作用、$U=0$）:

   • この転移は、相関の減衰が冪乗則（$d^{-2}$）から指数関数的減衰（$e^{-d/\xi}$）への変化によって特徴付けられる。
   • ネガティビティの対数微分 $\partial \log N^F / \partial t'$ は、臨界点（$t'=3$）において距離 $d$ に対して線形に増大するピークを示す。二階微分は、相関長指数 $\nu=1$ と一致するように、高さが $d$ に比例してスケールするピークを示す。
   • 大規模システムにおけるエンタングルメント・エントロピーの計算は、エリア則およびサブリーディングなスケーリング挙動を裏付けている。

2. MI-BI 転移（相互作用、$U=7$）:

   • この転移は O(3) ウィルソン・フィッシャー普遍性クラス（$\nu \approx 0.711$）に属する。
   • スピン・エンタングルメント ($N^S$) は最近接を超えると「突然死（sudden death）」を起こすが、フェルミオン的ネガティビティ ($N^F$) はゼロのままではない。
   • $N^F$ の対数微分は、転移において距離とともに増大するピークを示す。著者らは、これを、臨界的なスピン揺らぎが電荷セクターへ間接的に影響を与えているためであると考えている。これにより、電荷セクター自体はギャップがあっても、電荷相関の相関長の発散が駆動される。

3. SM-MI 転移（相互作用、$U$ を変化させる）:

   • この転移はグロス・ニューウェン・ハイゼンベルク普遍性クラスに属する。近年の大規模シミュレーションは $\nu \approx 1.11$ を示唆している。
   • $1 < \nu < 2$ という理論的期待と一致して、一次対数微分は発散するピークを示さないが、二次の対数微分は距離とともにスケールする明確なピークを示す。
   • スピン・エンタングルメントはここでも大きな距離で突然死を示し、この領域における距離スケーリングへの有用性が限定的であることを示している。

意義と主張
本論文は、小さな部分系間のフェルミオン的ネガティビティのスケーリングが、強相関フェルミオン系における量子相転移を検出するための、堅牢かつ普遍的な診断ツールとして機能することを主張している。

• 普遍性: 本手法は、対数微分の距離によるスケーリングを分析することにより、異なる普遍性クラス（ガウス、グロス・ニューウェン、O(3)）にわたる転移を検出することに成功している。これらの微分の挙動は、転移の相関長指数 $\nu$ と直接結びついている。
• 相互作用系への適用可能性: 大きなシステムサイズを必要とする従来の手法や厳密対角化とは異なり、本アプローチは、大きなシステムサイズへのアクセスが困難な相互作用モデルのQMCシミュレーションの制約条件下でも動作する。
• 秩序パラメータを超えて: 本手法は、相関長が発散する限り、従来の秩序パラメータを欠く転移に対しても有効である。
• 堅牢性: フェルミオン的ネガティビティは、スピン・エンタングルメントよりも優れていることが示された。なぜなら、それは大きな距離での突然死に苦しまず、電荷揺らぎを含むモデルの全物理を捉えることができるからである。

著者らは、このエンタングルメント尺度の「距離スケーリング」が、従来のスケーリング解析が計算量的に実行不可能な複雑な多体系において、臨界点を特定するための実用的かつ理論的根拠に基づいた手法を提供すると結論づけている。