ALPS 2.1 Examples:Paramagnetic Metal/ja

常磁性金属と外挿誤差

このチュートリアルでは、常磁性自己無頓着法を用いて相互作用$U=3D/\sqrt{2}$、温度$\beta =32 \sqrt{2}/D$の条件でBethe格子でのHubbard模型の計算をおこないます。自己エネルギーを求め、DMFTで説明したGeorges it et al.のFig.15と比較し、Hirsch-Fyeと厳密解との関係を調査します。Hirsch-Fyeアルゴリズムの２つの連続時間モンテカルロアルゴリズムCT-HYBとCT-INTは、離散化による誤差に悩むことはないでしょう。厳密解を再現できるといえます。

パラメータファイル、Pythonスクリプトはtutorials/dmft-06-paramagnetにあります。

dmft parm_hyb

and

dmft parm_int

コマンドラインでの実行、Pythonスクリプトでの実行はtutorial6a.py、tutorial6b.pyを参照してください。各DMFT反復$i$での自己エネルギーはselfenergy_iで表されます。計算結果をGeorges it et al.と比較してみてください。