[論文レビュー] Solutions of the Two-Dimensional Hubbard Model: Benchmarks and Results from a Wide Range of Numerical Algorithms

J. P. F. LeBlanc, Andrey E. Antipov+23 more|arXiv (Cornell University)|Jan 1, 2015

Physics of Superconductivity and Magnetism参考文献 137被引用数 581

ひとこと要約

本稿は、量子モンテカルロ法、図式モンテカルロ法、DMRG、動的平均場的手法を含む広範な数値的手法を用いて、2次元 Hubbard モデルの包括的なベンチマーク研究を提示している。弱い結合、半フィルイン、および半フィルインから大きく離れた状態では信頼できる結果が得られたが、半フィルイン付近および中程度の相互作用強さでは、競合する量子相の影響により、依然として一貫性のない不確実性が残存している。

ABSTRACT

Numerical results for ground-state and excited-state properties (energies, double occupancies, and Matsubara-axis self-energies) of the single-orbital Hubbard model on a two-dimensional square lattice are presented, in order to provide an assessment of our ability to compute accurate results in the thermodynamic limit. Many methods are employed, including auxiliary-field quantum Monte Carlo, bare and bold-line diagrammatic Monte Carlo, method of dual fermions, density matrix embedding theory, density matrix renormalization group, dynamical cluster approximation, diffusion Monte Carlo within a fixed-node approximation, unrestricted coupled cluster theory, and multireference projected Hartree-Fock methods. Comparison of results obtained by different methods allows for the identification of uncertainties and systematic errors. The importance of extrapolation to converged thermodynamic-limit values is emphasized. Cases where agreement between different methods is obtained establish benchmark results that may be useful in the validation of new approaches and the improvement of existing methods.

研究の動機と目的

• 2次元 Hubbard モデルを解くための多様な数値的手法の正確性と一貫性を評価すること。
• 異なる手法が一致するパラメータ領域を特定し、今後の理論的・アルゴリズム的開発のための信頼できるベンチマークを確立すること。
• 有限サイズ効果、図の切断、統計的サンプリングの不確実性に起因する不確実性を定量化すること。
• 特に半フィルイン付近および中程度の U における、系統的誤差と相の競合による手法間の不一致が顕在化するパラメータ空間の領域を強調すること。

提案手法

• 熱力学的極限における数値的に正確な結果を得るため、補助場量子モンテカルロ法（AFQMC）を用いた。
• Feynman 図を制御された近似で体系的に含めるために、ボトムラインおよびボルドライン図式モンテカルロ法（DiagMC）を適用した。
• 動的平均場理論を超えた長波長相関にアクセスするために、デュアルフェルミオン（DF）法を用いた。
• 有限クラスタにおける強い相関を、埋め込み自己エネルギー補正を用いて検討するために、密度行列埋め込み理論（DMET）および密度行列縮小群（DMRG）を用いた。
• 変分的手法と摂動的手法の比較のために、固定ノード拡散モンテカルロ法（FN-DMC）および非制限付きカップルクラスター（UCCSD）を用いた。
• 対称性の破れや不均一相を調査するために、多電子状態予測ハートリー・フォック法（MRPHF）とクラスタ動的平均場理論（DCA）を組み合わせた。

実験結果

リサーチクエスチョン

• RQ1 2次元 Hubbard モデルの基底状態エネルギーおよび二重占有数について、異なる数値的手法が熱力学的極限でどの程度一致するか。
• RQ2 さまざまな数値的手法における系統的誤差の主な要因は何か。また、それらは相互作用強さや電子密度に依存してどのように変化するか。
• RQ3 どのパラメータ領域（U, n）で複数の手法が誤差範囲内で一致するか。信頼できるベンチマークとしての妥当性が示される。
• RQ4 なぜ半フィルイン付近および中程度の U では不一致が持続するのか。これは数値的アーチファクトか、物理的相の競合によるものか。
• RQ5 数値的に正確な手法（例：AFQMC、DiagMC）の一致を用いて、DMRG や UCCSD などの近似手法を検証できるか。

主な発見

• 半フィルインにおいて、AFQMC、DiagMC および他の数値的に正確な手法の間で優れた一致が確認され、モット絶縁体相の信頼できるベンチマークが確立された。
• 弱い結合領域および半フィルインから大きく離れたキャリア濃度では、複数の手法が同じ基底状態エネルギーおよび二重占有数を、小さな不確実性で一致して得た。
• 半フィルイン付近および中程度の相互作用強さ（U ≈ 8–12）では、手法間で顕著な不一致が生じ、大きな不確実性と相の競合の可能性を示唆した。
• 弱い結合および半フィルインにおいて、Matsubara 軸上の自己エネルギーは AFQMC と DiagMC の両方で一貫した挙動を示し、動的性質の評価にこれらの手法を用いる妥当性が裏付けられた。
• 本研究では、DMRG や UCCSD の系統的誤差が中程度の U の領域で顕著に顕在化しており、特に半フィルイン付近では競合する秩序状態を捉えられていない可能性があると判明した。
• 著者らは、中程度の領域における残存する不確実性は、主に量子相転移に近接する物理的要因に起因しており、単なる数値的アーチファクトではないと結論づけた。