ALPS 2.1 Examples:Paramagnetic Metal/ja

From ALPS
Jump to: navigation, search

常磁性金属と外挿誤差

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

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

dmft parm_hyb

and

dmft parm_int

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