Continuous-time quantum Monte Carlo impurity solvers

Abstract Continuous-time quantum Monte Carlo impurity solvers are algorithms that sample the partition function of an impurity model using diagrammatic Monte Carlo techniques. The present paper describes codes that implement the interaction expansion algorithm originally developed by Rubtsov, Savkin, and Lichtenstein, as well as the hybridization expansion method developed by Werner, Millis, Troyer, et al. These impurity solvers are part of the ALPS-DMFT application package and are accompanied by an implementation of dynamical mean-field self-consistency equations for (single orbital single site) dynamical mean-field problems with arbitrary densities of states. Program summary Program title: dmft Catalogue identifier: AEIL_v1_0 Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AEIL_v1_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: ALPS LIBRARY LICENSE version 1.1 No. of lines in distributed program, including test data, etc.: 899 806 No. of bytes in distributed program, including test data, etc.: 32 153 916 Distribution format: tar.gz Programming language: C++ Operating system: The ALPS libraries have been tested on the following platforms and compilers: • Linux with GNU Compiler Collection (g++ version 3.1 and higher), and Intel C++ Compiler (icc version 7.0 and higher) • MacOS X with GNU Compiler (g++ Apple-version 3.1, 3.3 and 4.0) • IBM AIX with Visual Age C++ (xlC version 6.0) and GNU (g++ version 3.1 and higher) compilers • Compaq Tru64 UNIX with Compq C++ Compiler (cxx) • SGI IRIX with MIPSpro C++ Compiler (CC) • HP-UX with HP C++ Compiler (aCC) • Windows with Cygwin or coLinux platforms and GNU Compiler Collection (g++ version 3.1 and higher) RAM: 10 MB–1 GB Classification: 7.3 External routines: ALPS [1], BLAS/LAPACK, HDF5 Nature of problem: (See [2].) Quantum impurity models describe an atom or molecule embedded in a host material with which it can exchange electrons. They are basic to nanoscience as representations of quantum dots and molecular conductors and play an increasingly important role in the theory of “correlated electron” materials as auxiliary problems whose solution gives the “dynamical mean field” approximation to the self-energy and local correlation functions. Solution method: Quantum impurity models require a method of solution which provides access to both high and low energy scales and is effective for wide classes of physically realistic models. The continuous-time quantum Monte Carlo algorithms for which we present implementations here meet this challenge. Continuous-time quantum impurity methods are based on partition function expansions of quantum impurity models that are stochastically sampled to all orders using diagrammatic quantum Monte Carlo techniques. For a review of quantum impurity models and their applications and of continuous-time quantum Monte Carlo methods for impurity models we refer the reader to [2]. Additional comments: Use of dmft requires citation of this paper. Use of any ALPS program requires citation of the ALPS [1] paper. Running time: 60 s–8 h per iteration. References: [1] A. Albuquerque, F. Alet, P. Corboz, et al., J. Magn. Magn. Mater. 310 (2007) 1187. [2] http://arxiv.org/abs/1012.4474 , Rev. Mod. Phys., in press.

[1]  W. Krauth,et al.  Dynamical mean-field theory of strongly correlated fermion systems and the limit of infinite dimensions , 1996 .

[2]  E. Müller-Hartmann,et al.  Correlated fermions on a lattice in high dimensions , 1989 .

[3]  Matthias Troyer,et al.  Continuous-time solver for quantum impurity models. , 2005, Physical review letters.

[4]  S.Rombouts,et al.  QUANTUM MONTE CARLO METHOD FOR FERMIONS, FREE OF DISCRETIZATION ERRORS , 1999 .

[5]  D. Vollhardt,et al.  Correlated Lattice Fermions in High Dimensions , 1989 .

[6]  Georges,et al.  Numerical solution of the d= , 1992, Physical review letters.

[7]  I. S. Tupitsyn,et al.  Exact, complete, and universal continuous-time worldline Monte Carlo approach to the statistics of discrete quantum systems , 1997, cond-mat/9703200.

[8]  David J. Luitz,et al.  Weak-coupling continuous-time quantum Monte Carlo study of the single impurity and periodic Anderson models with s -wave superconducting baths , 2009, 0909.2656.

[9]  Emanuel Gull,et al.  Momentum-sector-selective metal-insulator transition in the eight-site dynamical mean-field approximation to the Hubbard model in two dimensions , 2009, 0909.1795.

[10]  J. C. Kimball,et al.  Perturbation Technique for the Anderson Hamiltonian , 1970 .

[11]  Emanuel Gull,et al.  Continuous-time auxiliary-field Monte Carlo for quantum impurity models , 2008, 0802.3222.

[12]  C. Marianetti,et al.  Electronic structure calculations with dynamical mean-field theory , 2005, cond-mat/0511085.

[13]  A. I. Lichtenstein,et al.  Continuous-time quantum Monte Carlo method for fermions: Beyond auxiliary field framework , 2004 .

[14]  Fye,et al.  Monte Carlo method for magnetic impurities in metals. , 1986, Physical review letters.

[15]  M. Troyer,et al.  Diagrammatic Monte Carlo for correlated fermions , 2009, 0907.0863.

[16]  Caffarel,et al.  Exact diagonalization approach to correlated fermions in infinite dimensions: Mott transition and superconductivity. , 1994, Physical review letters.

[17]  T. Pruschke,et al.  Numerical renormalization group method for quantum impurity systems , 2007, cond-mat/0701105.

[18]  F. Assaad,et al.  Diagrammatic determinantal quantum Monte Carlo methods: Projective schemes and applications to the Hubbard-Holstein model , 2007, cond-mat/0702455.

[19]  Emanuel Gull,et al.  Momentum-selective metal-insulator transition in the two-dimensional Hubbard model: An 8-site dynamical cluster approximation study , 2009 .

[20]  Anomalous normal-state properties of high-Tc superconductors: intrinsic properties of strongly correlated electron systems? , 1995, supr-con/9502001.

[21]  A. I. Lichtenstein,et al.  Continuous-time quantum Monte Carlo method for fermions , 2005 .

[22]  H. Bethe.  Statistical Theory of Superlattices , 1935 .

[23]  A. Honecker,et al.  The ALPS Project: Open Source Software for Strongly Correlated Systems , 2004, cond-mat/0410407.

[24]  Kristjan Haule.  Quantum Monte Carlo impurity solver for cluster dynamical mean-field theory and electronic structure calculations with adjustable cluster base , 2007 .

[25]  T. Pruschke,et al.  The Anderson model with finite Coulomb repulsion , 1989 .

[26]  Emanuel Gull,et al.  Local order and the gapped phase of the Hubbard model: A plaquette dynamical mean-field investigation , 2008, 0805.3778.

[27]  A. I. Lichtenstein,et al.  Dual fermion approach to nonlocal correlations in the Hubbard model , 2006, cond-mat/0612196.

[28]  P ? ? ? ? ? ? ? % ? ? ? ? , 1991 .

[29]  Emanuel Gull,et al.  Bold-line diagrammatic Monte Carlo method: General formulation and application to expansion around the noncrossing approximation , 2010, 1004.0724.

[30]  Philipp Werner,et al.  Efficient dynamical mean field simulation of the Holstein-Hubbard model. , 2007, Physical review letters.

[31]  Emanuel Gull,et al.  Performance analysis of continuous-time solvers for quantum impurity models , 2007 .

[32]  T. Pruschke,et al.  Quantum cluster theories , 2004, cond-mat/0404055.

[33]  A. Honecker,et al.  The ALPS project release 1.3: Open-source software for strongly correlated systems , 2007 .

[34]  Boris V. Svistunov,et al.  POLARON PROBLEM BY DIAGRAMMATIC QUANTUM MONTE CARLO , 1998 .

[35]  Philipp Werner,et al.  Hybridization expansion impurity solver: General formulation and application to Kondo lattice and two-orbital models , 2006 .