[論文レビュー] Taming the sign problem of the finite density Hubbard model via Lefschetz thimbles

Maksim Ulybyshev, Christopher Winterowd+1 more|arXiv (Cornell University)|Jun 6, 2019

Complex Systems and Time Series Analysis参考文献 1被引用数 12

ひとこと要約

本稿では、ヘキサゴナル格子上の有限密度 Hubbard モデルにおける深刻な符号問題を軽減するために、Lefschetz のテイブルと正則勾配フローを組み合わせる手法を提案する。複素化された場の空間における統合経路を、臨界点に接続されたテイブル多様体に沿って変形することにより、符号のフラクチュエーションを低減し、安定な量子モンテカルロシミュレーションを可能にした。この手法により、従来の手法が失敗する低温および大きな化学ポテンシャル領域においても、正確対角化との整合性を達成した。

ABSTRACT

We study the sign problem in the Hubbard model on the hexagonal lattice away from half-filling using the Lefschetz thimbles method. We identify the saddle points, reduce their amount, and perform quantum Monte Carlo (QMC) simulations using the holomorphic gradient flow to deform the integration contour in complex space. Finally, the results are compared against exact diagonalization (ED). We show that the sign problem can be substantially weakened, even in the regime with low temperature and large chemical potential, where standard QMC techniques exhibit an exponential decay of the average sign.

研究の動機と目的

• 有限密度 Hubbard モデルにおける深刻な符号問題に、特に半フィルインでない領域、低温および大きな化学ポテンシャルの下で、標準的 QMC が失敗する状況に対処すること。
• Lefschetz のテイブルが、複素化された場空間への統合経路の変形を通じて、複雑な作用積分を単純化できるかどうかを調査すること。
• 正則勾配フロー手法が、符号フラクチュエーションを低減し、QMC シミュレーションにおけるエルゴード性を回復する効果をテストすること。
• 小スケール格子における正確対角化との比較を通じて、結果の妥当性を検証し、手法の正確性を保証すること。

提案手法

• 本手法は Lefschetz のテイブル分解を採用し、作用の複素臨界点に接続されたテイブル多様体上の積分の和として分配関数を表現する。
• 臨界点は複素場空間における ∂S/∂Φ = 0 を解くことで特定され、そのテイブル構造の関連性が分析される。
• 正則勾配フローを用いて、元の実数の統合多様体を複素的なテイブルに類似した多様体へと変形し、作用の虚数部を低減する。
• ハミルトニアン・モンテカルロと共役勾配ソルバー（HMC-CG）を用いて量子モンテカルロシミュレーションを実行し、フロー変形された配置の平均を通じて物理量を評価する。
• 符号問題の度合いは、⟨cos(Im S)⟩、⟨cos(Im ln det J))⟩、および結合符号 ΣG といった指標で定量化され、ジャコビアンのフラクチュエーションは分散 D_J でモニタリングされる。
• シミュレーションは U=2κ、μ=κ、β=20 および β=30 の 2×2×2×384 格子で実施され、正確対角化との比較が行われる。

実験結果

リサーチクエスチョン

• RQ1 Lefschetz のテイブルは、半フィルインでない領域において、有限密度 Hubbard モデルにおける符号問題の深刻さを効果的に軽減できるか？
• RQ2 正則勾配フロー手法は、統合多様体の構造をどのように変化させ、符号フラクチュエーションを抑制するか？
• RQ3 どのフロー時間および α パrameter において、テイブル境界への衝突またはジャコビアンのゼロ点への到達によってエルゴード性が崩壊するか？
• RQ4 低温（大きな β）および真空でない臨界点が存在する状況においても、この手法は安定的かつ正確に機能するか？
• RQ5 ジャコビアンのフラクチュエーションおよびその分散 D_J が、非常に低温においてシミュレーションの信頼性をどの程度損なうか？

主な発見

• α=0.8 の場合、HMC-GF シミュレーションによる運動エネルギーおよびスピン相関関数の結果は、誤差範囲内で正確対角化と一致し、正しさと安定性が確認された。
• α=1.0 の場合、複数の関連するテイブルとジャコビアンのゼロ点への衝突によりエルゴード性が失敗し、適切なサンプリングが不可能になったため、ED からのずれが生じた。
• ジャコビアンの分散 D_J は、α=0.8 のとき 1.011 から β=30 のとき 1.68 に増加しており、低温におけるフラクチュエーションの増大が示され、極めて低温では手法の限界が顕在化する可能性がある。
• α=0.8 では結合符号 ΣG が安定的かつ非ゼロを維持しており、符号問題の効果的な抑制が示されたが、α=1.0 ではその性能が劣化した。
• 本手法により符号問題が著しく弱体化した：従来の BSS-QMC ではこれらのパrameter で指数的符号減少が見られたが、HMC-GF では使用可能な統計が維持された。
• β=20 のとき 𝒯=0.05 から β=30 のとき 𝒯=0.08 にフロー時間を見積もる必要があり、同じ程度の Im S フラクチュエーションを維持するには、適応的フロー調整が必要であることが示された。