[論文レビュー] Strong-coupling formula for momentum-dependent susceptibilities in dynamical mean-field theory

Junya Otsuki, Kazuyoshi Yoshimi+2 more|arXiv (Cornell University)|Dec 28, 2018

Physics of Superconductivity and Magnetism参考文献 38被引用数 28

ひとこと要約

本稿では、動的平均場理論（DMFT）における運動量依存する磁化率を計算するための強結合近似（SCL式）を提案する。この手法は、局所的相関を用いてバーティカル・サルツァー方程式を簡略化することで、計算コストを著しく削減する。物理的現象として運動エネルギー交換やRKKY相互作用を弱い結合から強い結合まで正確に捉えることができ、多軌道系におけるスピンおよび軌道磁化率の効率的計算を可能にする。本手法は3段階の実用的近似レベルを提供する。

ABSTRACT

Computing momentum-dependent susceptibilities in the dynamical mean-field theory (DMFT) requires solving the Bethe-Salpeter equation, which demands large computational cost. Exploiting the strong-coupling feature of local fluctuations, we derive a simplified formula that can be solved at a considerably lower cost. The validity and the physical meaning of the formula are confirmed by deriving the effective intersite interactions in the strong-coupling limit, such as the kinetic exchange and RKKY interactions. Furthermore, numerical calculations for single-orbital and multiorbital models demonstrate surprisingly wider applicability including weak-coupling region. Based on this formula, we propose three levels of practical approximations that can be chosen depending on complexity of problems. Simpler evaluations of spin and orbital susceptibilities in multiorbital systems thus become possible within DMFT.

研究の動機と目的

• バーティカル・サルツァー方程式に起因する計算コストの高さにより、通常は計算が非現実的となる、DMFTにおける運動量依存磁化率を計算するための計算効率の高い手法を開発すること。
• 局所的フラクチュエーションの強結合物理的特徴を活用し、運動エネルギー交換やRKKYといった主要な物理的相互作用を保持する簡略化された公式を導出すること。
• 多軌道系への拡張と、実際の材料に応用可能なスケーラブルな近似スキームの提供。
• 弱結合領域を含む広範な結合範囲において、本手法の妥当性を検証すること。特に、予想に反して弱結合領域でも高い精度を示す点に注目。

提案手法

• 局所的頂点関数の強結合的特徴を活用して、バーティカル・サルツァー方程式を簡略化し、RPAに類似した形に還元するが、入力として完全に相関された局所的磁化率を用いる。
• SCL式を導入：χ(q) ≈ χloc + χloc [X0(q) - X0,loc] χloc ここで、χloc は局所的磁化率、X0(q) は非相互作用の2粒子グリーン関数。
• 既知の強結合物理的現象（ハーバード模型における運動エネルギー交換、周期的アンダーソン模型におけるRKKY）を正確に再現できることを示して、公式の妥当性を検証。
• 3段階の実用的近似レベルを提案：SCL1（局所的 χ のみ）、SCL2（運動量依存の X0）、SCL3（u0(iω) の2極近似を用いて χloc を効率的に推定）。
• 連続時間量子モンテカルロ法とハッバー-I近似をイムプリンティ・ソルバーとして用い、周波数カットオフとフィッティング手順を導入して数値的安定性を確保。
• 多軌道系における χloc の効率的評価を可能にするために、主要モード関数 u0(iω) の2極近似を採用。

実験結果

リサーチクエスチョン

• RQ1 DMFTにおけるバーティカル・サルツァー方程式から、計算コストを著しく削減するが、主要な物理的性質を保持する簡略化された公式を導出可能か？
• RQ2 提案されたSCL式は、ハーバード模型と周期的アンダーソン模型における既知の強結合相互作用（運動エネルギー交換、RKKY）を正確に記述できるか？
• RQ3 SCL式は強結合領域を超えて、特に弱結合領域および多軌道系において、どの程度の性能を示すか？
• RQ4 本手法は、精度と複雑さを段階的に向上させる複数の近似レベルへ体系的に拡張可能か、実用的な材料シミュレーションへの応用が可能か？

主な発見

• SCL式は、弱い結合から強い結合まで広範なU/W範囲において、正確な運動量依存磁化率を提供する。これは、このような近似が強結合限界でのみ有効であると予想されるのとは対照的である。
• 本手法は、ハーバード模型における運動エネルギー交換相互作用および周期的アンダーソン模型におけるRKKY相互作用を正確に再現しており、物理的妥当性が裏付けられている。
• 2軌道模型の数値結果から、u0(iω) の2極フィットを用いたSCL3近似が、完全なバーティカル・サルツァー解と良好な一致を示す一方で、計算コストを大幅に削減できることを示している。
• u0(iω) の2極近似は、主要モード構造を効果的に捉えており、E± ≈ E±,min と近似されるフィッティングパラメータが、近似における最小励起エネルギーの妥当性を裏付けている。
• SCL3スキームでは、s0 の直接計算を回避し、χloc ∝ s0|u0|2 の関係を用いることで、完全なイムプリンティ問題を解かずに効率的な評価が可能になっている。
• 公式は2粒子和則を近似的に満たし、弱結合領域であっても χloc を信頼できる入力として使用可能であることが示唆されている。