[論文レビュー] Universal term of Entanglement Entropy in the $π$-flux Hubbard model

Yuan Da Liao|arXiv (Cornell University)|Jul 20, 2023

Machine Learning in Materials Science被引用数 11

ひとこと要約

本論文は、行列式量子モンテカルロ法（DQMC）内に制御可能な増分アルゴリズムを導入し、相互作用フェルミオンの第2 Rényi エンタングルメントエントロピー（EE）を正確に計算する。指数的な分散問題に対処し、追加の構成空間を必要とせず、並列化可能でスケーラブルなEE計算を可能にする。

ABSTRACT

Researchers in physical science aim to uncover universal features in strongly interacting many-body systems, often hidden in complicated observables like entanglement entropy (EE). The non-local nature of EE makes it challenging to compute numerically, necessitating the development of an unbiased and convenient algorithm. In this paper, we use quantum Monte Carlo to reveal that the coefficient of variation in direct EE calculations increases exponentially with system size, leading to inaccuracies. To address this issue, we develop a power incremental algorithm and a technique for straightforwardly calculating the universal term of EE, successfully evaluating the EE of a 2D Hubbard model. Our numerical results demonstrate the consistency of the universal coefficient of EE from sharp corners at the Gross-Neveu quantum critical point and for free Dirac fermions. Our method can also be applied to other unstable observables, such as partition functions, entanglement spectra, and negativity, thereby fostering computational and theoretical progress.

研究の動機と目的

• 強く相関した系におけるエンタングルメントエントロピーの普遍的特徴を探る。
• 相互作用フェルミオンの第2 Rényi EEをDQMCで統計的に安定かつ実用的に計算する方法を開発する。
• 追加の構成空間に依存せず、増分ステップの数を定量的に扱える手段を提供する。
• グリーン関数の行列式に関連する他の量（例：エンタングルメントスペクトル、ネゲティビティ）にもアプローチを拡張する。

提案手法

• GroverのEE定義における det g_A^{s1,s2} の不安定性を分析し、系サイズとともに指数的分散が生じることを同定する。
• W(λ_k, det g_A^{s1,s2}) = (det g_A^{s1,s2})^{λ_k} および λ_k = k/N_λ を用いた制御可能な増分スキームを導入する。
• log det g_A^{s1,s2} が正規分布に従うことを示し、その平均 μ(L) と分散 σ^2(L) をべき法則として適合させ、分散を制御するために N_λ ≈ |μ| を設定する。
• (det g_A^{s1,s2})^{λ_k} および (det g_A^{s1,s2})^{1/N_λ} を用いて分散を適切に抑える並列化可能な増分比 Z(λ_{k+1})/Z(λ_k) を導出する。
• μ(L) のスケーリングから N_λ ∼ 0.5 L^{1.35} を定量的に決定し、追加の構成空間を増やすことなくプロジェクターDQMC内でこのアプローチを実装する。
• N_λ が増えるとEEの収束を示し、従来法と比較して計算効率が向上することを示す。

実験結果

リサーチクエスチョン

• RQ1 相互作用フェルミオンのエンタングルメントエントロピーを、追加の構成空間を導入せずにDQMCで正確に計算できるか。
• RQ2 det g_A^{s1,s2} の分布はどうなるか、log det g の正規分布特性を分散制御に利用できるか。
• RQ3 最適な増分ステップ数 N_λ の系サイズ L に対するスケーリングはどうなるか、計算コストにどう影響するか。
• RQ4 グリーン関数の行列式に依存する他の量（エンタングルメントスペクトル、ネゲティビティなど）へ方法を拡張できるか。

主な発見

• det g_A^{s1,s2} の直接評価は、結合が大きく系サイズが大きくなると希なスパイクと不安定な統計を示す。
• log det g_A^{s1,s2} は正規分布に従い、μ(L) ~ -0.50 L^{1.35} および σ(L) ~ 0.67 L^{0.77}。
• det g の変動係数は指数的に増加するが、(det g)^{1/N_λ} または log det g のCVは N_λ が L^{1.35} にスケールする場合に系サイズとともに減衰する。
• (det g)^{λ_k} および (det g)^{1/N_λ} を用いた増分スキームは、追加の構成空間を必要とせず、増分間で並列化可能な正確なEEの計算を可能にする。
• この手法は標準的なDQMCと同程度のCPU時間で信頼性の高い第2 Rényi EE推定を達成し、N_λ を選ぶ定量的な方法（約 0.5 L^{1.35}）を提供する。
• このアプローチは、エンタングルメントスペクトルやネゲティビティなど、グリーン関数の行列式に結びつく他の観測量へ拡張可能である。