Difference between revisions of "ALPS 2 Examples:Paramagnetic Metal"
(5 intermediate revisions by 2 users not shown) | |||
Line 4: | Line 4: | ||
=Paramagnetic metal and extrapolation errors= | =Paramagnetic metal and extrapolation errors= | ||
− | In this example we simulate the Hubbard model on the Bethe lattice with interaction <math>U=3D/\sqrt{2}</math> at a temperature <math>\beta =32 \sqrt{2}/D</math> using a paramagnetic self-consistency. We will calculate the self-energy and compare it to Fig. 15 in the DMFT review by [http://dx.doi.org/10.1103/RevModPhys.68.13 Georges ''it et al.''], where Hirsch-Fye and Exact | + | In this example we simulate the Hubbard model on the Bethe lattice with interaction <math>U=3D/\sqrt{2}</math> at a temperature <math>\beta =32 \sqrt{2}/D</math> using a paramagnetic self-consistency. We will calculate the self-energy and compare it to Fig. 15 in the DMFT review by [http://dx.doi.org/10.1103/RevModPhys.68.13 Georges ''it et al.''], where Hirsch-Fye and Exact Diagonalization results are shown for the same system. In contrast to the Hirsch-Fye algorithm the two Continuous time Monte Carlo algorithms CT-HYB and CT-INT do not suffer from discretization errors and reproduce the ED-results. |
− | The parameter files and python scripts are located in the directory <tt>tutorials/dmft-06-paramagnet</tt> in your ALPS install directory. You can run the simulations by executing | + | The parameter files and python scripts are located in the subdirectories <tt>hyb</tt> and <tt>int</tt> of the directory <tt>tutorials/dmft-06-paramagnet</tt> in your ALPS install directory. You can run the simulations by executing (for the hybridization expansion version) |
− | + | alpspython tutorial6a.py | |
− | and | + | (vispython on Mac) and (for the interaction expansion version) |
− | + | alpspython tutorial6b.py | |
− | + | At each DMFT iteration <math>i</math> the self-energy is written to the file <tt>selfenergy_i</tt>. Plot the converged self-energy and compare your results to Fig. 15 in [http://dx.doi.org/10.1103/RevModPhys.68.13 Georges ''it et al.'']. Alternatively you may use the python script for this task as it is done in the tutorial [[ALPS 2 Tutorials:DMFT-02 Hybridization]]. |
Latest revision as of 16:12, 17 September 2013
Languages: |
English • 日本語 (ja) • 繁體中文 (zh-tw) • 简体中文 (zh) |
Paramagnetic metal and extrapolation errors
In this example we simulate the Hubbard model on the Bethe lattice with interaction at a temperature using a paramagnetic self-consistency. We will calculate the self-energy and compare it to Fig. 15 in the DMFT review by Georges it et al., where Hirsch-Fye and Exact Diagonalization results are shown for the same system. In contrast to the Hirsch-Fye algorithm the two Continuous time Monte Carlo algorithms CT-HYB and CT-INT do not suffer from discretization errors and reproduce the ED-results.
The parameter files and python scripts are located in the subdirectories hyb and int of the directory tutorials/dmft-06-paramagnet in your ALPS install directory. You can run the simulations by executing (for the hybridization expansion version)
alpspython tutorial6a.py
(vispython on Mac) and (for the interaction expansion version)
alpspython tutorial6b.py
At each DMFT iteration the self-energy is written to the file selfenergy_i. Plot the converged self-energy and compare your results to Fig. 15 in Georges it et al.. Alternatively you may use the python script for this task as it is done in the tutorial ALPS 2 Tutorials:DMFT-02 Hybridization.