ALPS 2 Tutorials:DMFT-07 Hirsch-Fye

Revision as of 16:04, 17 September 2013 by Jimriska (talk | contribs)

(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to: navigation, search

Tutorial 07: Hirsch Fye Code

We start by running a discrete time Monte Carlo code: the Hirsch Fye code. As an example we reproduce Fig. 11 in the DMFT review by Georges it et al.. The series of six curves shows how the system, a Hubbard model on the Bethe lattice with interaction U=3D/\sqrt{2} at half filling, enters an antiferromagnetic phase upon cooling.

The Hirsch Fye algorithm is described in here, and this review also provides an open source implementation for the codes. More information can also be found in Blümer's PhD. While many improvements have been developed (see e.g. Alvarez08 or Nukala09), the algorithm has been replaced by continuous-time algorithms eliminating systematic discretization errors.

The Hirsch Fye simulation will run for about 20 seconds per iteration. The python scripts for running this simulation can be found in the directory tutorials/dmft-07-hirschfye, a short version running simulations at 2 temperatures only (takes 5 minutes) or a version reproducing all 6 temperatures For evaluation you may use the same script as it is described in the tutorial ALPS_2_Tutorials:DMFT-02_Hybridization.

If you want to use Vistrails to run your DMFT-simulations, you can use dmft-07-hirschfye.vt.

The main parameters are same as those described in the tutorial ALPS_2_Tutorials:DMFT-02_Hybridization.

As a discrete time method, HF suffers from \Delta\tau - errors. Pick a set of parameters and run it for sucessively larger N.

Tutorial by Emanuel - Please don't hesitate to ask!