[論文レビュー] The plain and simple parquet approximation: single- and multi-boson exchange in the two-dimensional Hubbard model

Friedrich Krien, Anna Kauch|arXiv (Cornell University)|Feb 1, 2022

Physics of Superconductivity and Magnetism参考文献 51被引用数 22

ひとこと要約

本稿では、弱い結合における2次元 Hubbard モデルに対して、数値的に効率的なボソン化されたパーケット近似の実装を提示する。16×16格子上で運動量および周波数の完全な分解能を達成可能であり、パーケット方程式をボソン的交換図（単一ボソンおよび多ボソン交換）の観点から再定式化することで、収束が速くなり、計算コストが低減される。これにより、長距離にわたる反強磁性相関が存在する状況下でも、頂点関数の偏りのない評価と、截断ユニティ近似のベンチマークが可能になる。

ABSTRACT

The parquet approach to vertex corrections is unbiased but computationally demanding. Most applications are therefore restricted to small cluster sizes or rely on various simplifying approximations. We have recently shown that the bosonization of the parquet diagrams provides interpretative and algorithmic advantages over the original purely fermionic formulation. Here we present first results of the numerical implementation of this method by applying it to the half-filled Hubbard model on the square lattice at weak coupling. The improved algorithmic performance allows us to evaluate the parquet approximation for a $16 imes16$ lattice, retaining the full momentum and frequency structure of the various vertex functions. We discuss their symmetries and consider parametrizations of their momentum dependence using the truncated unity approximation.

研究の動機と目的

• 本稿の目的は、従来のフェルミオン的定式化における計算的ボトルネックを克服することで、Hubbard モデルに対する偏りのない大規模パーケット計算を可能にすることである。
• ボソン化されたパーケット形式が、標準的なパーケット実装の実用的で効率的な代替手段であることを検証することである。
• 截断ユニティ近似をボソン化された頂点関数に適用した際の収束性と精度を調査することである。
• 特にスピンチャンネルにおいて、長距離にわたる反強磁性相関が存在する状況下での手法の性能をベンチマークすることである。
• 将来の関連研究における頂点補正、Ward 恒等式、和則の検証のための基準フレームワークを提供することである。

提案手法

• 著者らは、U-不変頂点関数（˜Λ = Λ − U）を用いたパーケット方程式の再定式化により、頂点補正のボソン化された記述を可能にしている。
• 頂点関数を単一ボソン交換（SBE）と多ボソン交換（M）の寄与に分離し、SBE は太字の γ および W 量で表される。
• SBE 図式は正確に取り扱い、その漸近的減衰を用いて和計算におけるマツブア周波数の必要数を削減している。
• 頂点関数の完全な運動量および周波数依存性を保持しており、パラメータ化には截断ユニティ（TU）近似を M および ∆ 関数に適用している。
• アルゴリズムは、半充填状態における16×16正方格子に実装されており、運動量および周波数構造の完全な分解能が得られる。
• 截断ユニティの収束は、形式因子の数（Nℓ）を変化させることで評価され、異なるパラメータ化スキーム（Φsp、Msp、Msp + ∆spR）との比較がなされている。

実験結果

リサーチクエスチョン

• RQ1 ボソン化されたパーケット形式は、大規模 Hubbard モデル計算において、標準的なフェルミオン的定式化と比較して、どの程度計算効率が向上するか？
• RQ2 長距離にわたる反強磁性相関が存在する状況下でも、截断ユニティ近似がボソン化された頂点関数に適用された際の収束度合いはどの程度か？
• RQ3 弱い結合におけるスピンチャンネルにおいて、頂点関数（Φsp、Msp、∆sp）の運動量および周波数依存性はどのように振る舞うか？
• RQ4 截断ユニティの収束度合いから測定した場合、Msp の相対的重要性が相関長 ξ に依存するかどうか？
• RQ5 ボソン化されたフレームワークにおいて、截断ユニティ近似を用いて頂点の漸近的挙動を低周波数領域に効果的に拡張できるか？

主な発見

• ボソン化されたパーケット手法により、16×16格子上で完全なフェルミオン的パーケット計算が初めて実現され、運動量および周波数の完全な分解能が保持された。
• 相関長 ξ が増加しても、Msp と ∆sp の相対的重要性は概ね一定のままであるため、Msp に対する截断ユニティの収束が安定していることが示された。
• Msp + ∆spR に截断ユニティを適用した場合、U/t が大きくなる（例：U=4t）と収束が著しく遅くなるため、これらの項を組み合わせることでパラメータ化の効率が低下することが示された。
• 截断ユニティを Msp のみに適用した場合、収束が最も速く、これは Msp が ∆sp よりも弱い運動量依存性を示すという事実と整合している。
• 頂点の漸近的挙動は、高周波数に限らず、低周波数領域の頂点関数のパラメータ化に信頼性を持って用いることができる。
• 本手法により、頂点関数の偏りのない評価が可能となり、将来の関連研究における Ward 恒等式や和則の検証が促進される。