Difference between revisions of "ALPS 2 Examples:Paramagnetic Metal"

From ALPS
Jump to: navigation, search
m (Text replace - "/tutorials2.0.0/" to "/tutorials2.1.0/")
Line 12: Line 12:
 
  dmft parm_int
 
  dmft parm_int
  
on the command line or by running the python scripts  [http://alps.comp-phys.org/static/tutorials2.0.0/dmft-06-paramagnet/tutorial6a.py tutorial6a.py] and [http://alps.comp-phys.org/static/tutorials2.0.0/dmft-06-paramagnet/tutorial6b.py tutorial6b.py] with the vispython interpreter. 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.''].
+
on the command line or by running the python scripts  [http://alps.comp-phys.org/static/tutorials2.1.0/dmft-06-paramagnet/tutorial6a.py tutorial6a.py] and [http://alps.comp-phys.org/static/tutorials2.1.0/dmft-06-paramagnet/tutorial6b.py tutorial6b.py] with the vispython interpreter. 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.''].

Revision as of 22:08, 10 May 2012


Paramagnetic metal and extrapolation errors

In this example we simulate the Hubbard model on the Bethe lattice with interaction U=3D/\sqrt{2} at a temperature \beta =32 \sqrt{2}/D 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 Diagonalizationr 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 tutorials/dmft-06-paramagnet in your ALPS install directory. You can run the simulations by executing

dmft parm_hyb

and

dmft parm_int

on the command line or by running the python scripts tutorial6a.py and tutorial6b.py with the vispython interpreter. At each DMFT iteration i 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..