Hartmut Hafermann

TRIQS: A Toolbox for Research on Interacting Quantum Systems

We present the TRIQS library, a Toolbox for Research on Interacting Quantum Systems. It is an open-source, computational physics library providing a framework for the quick development of applications in the field of many-body quantum physics, and in particular, strongly-correlated electronic systems. It supplies components to develop codes in a modern, concise and efficient way: e.g. Greens function containers, a generic Monte Carlo class, and simple interfaces to HDF5. TRIQS is a C++/Python library that can be used from either language. It is distributed under the GNU General Public License (GPLv3). State-of-the-art applications based on the library, such as modern quantum many-body solvers and interfaces between density-functional-theory codes and dynamical mean-field theory (DMFT) codes are distributed along with it.

Orthogonal polynomial representation of imaginary-time Green's functions

We study the expansion of single-particle and two-particle imaginary-time Matsubara Greens functions of quantum impurity models in the basis of Legendre orthogonal polynomials. We discuss various applications within the dynamical mean-field theory (DMFT) framework. The method provides a more compact representation of the Greens functions than standard Matsubara frequencies, and therefore significantly reduces the memory-storage size of these quantities. Moreover, it can be used as an efficient noise filter for various physical quantities within the continuous-time quantum Monte Carlo impurity solvers recently developed for DMFT and its extensions. In particular, we show how to use it for the computation of energies in the context of realistic DMFT calculations in combination with the local density approximation to the density functional theory (LDA+DMFT) and for the calculation of lattice susceptibilities from the local irreducible vertex function.

Fermi Condensation Near van Hove Singularities Within the Hubbard Model on the Triangular Lattice

The proximity of the Fermi surface to van Hove singularities drastically enhances interaction effects and leads to essentially new physics. In this work we address the formation of flat bands (Fermi condensation) within the Hubbard model on the triangular lattice and provide a detailed analysis from an analytical and numerical perspective. To describe the effect we consider both weak-coupling and strong-coupling approaches, namely the renormalization group and dual fermion methods. It is shown that the band flattening is driven by correlations and is well pronounced even at sufficiently high temperatures, of the order of 0.1-0.2 of the hopping parameter. The effect can therefore be probed in experiments with ultracold fermions in optical lattices.

Efficient Perturbation Theory for Quantum Lattice Models

We present a novel approach to long-range correlations beyond dynamical mean-field theory, through a ladder approximation to dual fermions. The new technique is applied to the two-dimensional Hubbard model. We demonstrate that the transformed perturbation series for the nonlocal dual fermions has superior convergence properties over standard diagrammatic techniques. The critical Néel temperature of the mean-field solution is suppressed in the ladder approximation, in accordance with quantum Monte Carlo results. An illustration of how the approach captures and allows us to distinguish short- and long-range correlations is given.

Superconductivity, antiferromagnetism and phase separation in the two-dimensional Hubbard model: A dual-fermion approach

The dual-fermion approach offers a way to perform diagrammatic expansion around the dynamical mean-field theory. Using this formalism, the influence of antiferromagnetic fluctuations on the self-energy is taken into account through ladder-type diagrams in the particle-hole channel. The resulting phase diagram for the (quasi-)two-dimensional Hubbard model exhibits antiferromagnetism and d-wave superconductivity. Furthermore, a uniform charge instability, i.e., phase separation, is obtained in the low doping regime around the Mott insulator. We also examine spin/charge density wave fluctuations including d-wave symmetry. The model exhibits a tendency towards an unconventional charge density-wave, but no divergence of the susceptibility is found.

Superperturbation solver for quantum impurity models

We present a very efficient solver for the general Anderson impurity problem. It is based on the perturbation around a solution obtained from exact diagonalization, using a small number of bath sites. We formulate a perturbation theory which is valid for both weak and strong coupling and interpolates between these limits. Good agreement with numerically exact quantum Monte Carlo results is found for a single bath site over a wide range of parameters. In particular, the Kondo resonance in the intermediate-coupling regime is well reproduced for a single bath site and the lowest-order correction. The method is particularly suited for low temperatures, alleviates analytical continuation of imaginary time data due to the absence of statistical noise and can be generalized to obtain dynamical quantities directly on the real axis.

Efficient implementation of the continuous-time hybridization expansion quantum impurity solver

Abstract Strongly correlated quantum impurity problems appear in a wide variety of contexts ranging from nanoscience and surface physics to material science and the theory of strongly correlated lattice models, where they appear as auxiliary systems within dynamical mean-field theory. Accurate and unbiased solutions must usually be obtained numerically, and continuous-time quantum Monte Carlo algorithms, a family of algorithms based on the stochastic sampling of partition function expansions, perform well for such systems. With the present paper we provide an efficient and generic implementation of the hybridization expansion quantum impurity solver, based on the segment representation. We provide a complete implementation featuring most of the recently developed extensions and optimizations. Our implementation allows one to treat retarded interactions and provides generalized measurement routines based on improved estimators for the self-energy and for vertex functions. The solver is embedded in the ALPS-DMFT application package. Program summary Program title: ct-hyb Catalogue identifier: AEOL_v1_0 Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AEOL_v1_0.html Program obtainable from: CPC Program Library, Queen’s University, Belfast, N. Ireland Licensing provisions: Use of the hybridization expansion impurity solvers requires citation of this paper. Use of any ALPS program requires citation of the ALPS [1] paper. No. of lines in distributed program, including test data, etc.: 650044 No. of bytes in distributed program, including test data, etc.: 20553265 Distribution format: tar.gz Programming language: C++/Python. Computer: Desktop PC, high-performance computers. Operating system: Unix, Linux, OSX, Windows. Has the code been vectorized or parallelized?: Yes, MPI parallelized. RAM: 1 GB Classification: 7.3. External routines: ALPS [1, 2, 3], BLAS [4, 5], LAPACK [6], HDF5 [7] Nature of problem: Quantum impurity models were originally introduced to describe a magnetic transition metal ion in a non-magnetic host metal. They are widely used today. In nanoscience they serve as representations of quantum dots and molecular conductors. In condensed matter physics, they are playing an increasingly important role in the description of strongly correlated electron materials, where the complicated many-body problem is mapped onto an auxiliary quantum impurity model in the context of dynamical mean-field theory, and its cluster and diagrammatic extensions. They still constitutes a non-trivial many-body problem, which takes into account the (possibly retarded) interaction between electrons occupying the impurity site. Electrons are allowed to dynamically hop on and off the impurity site, which is described by a time-dependent hybridization function. Solution method: The quantum impurity model is solved using a continuous-time quantum Monte Carlo algorithm which is based on a perturbation expansion of the partition function in the impurity–bath hybridization. Monte Carlo configurations are represented as segments on the imaginary time interval and individual terms correspond to Feynman diagrams which are stochastically sampled to all orders using a Metropolis algorithm. For a detailed review on the method, we refer the reader to [8]. Running time: 1–8xa0h. References: [1] B. Bauer, L. D. Carr, H. G. Evertz, A. Feiguin, J. Freire, S. Fuchs, L. Gamper, J. Gukelberger, E. Gull, S.xa0Guertler, A. Hehn, R. Igarashi, S. V. Isakov, D. Koop, P. N. Ma, P. Mates, H. Matsuo, O. Parcollet, G.xa0Pawlowski, J. D. Picon, L. Pollet, E. Santos, V. W. Scarola, U. Schollwock, C. Silva, B. Surer, S. Todo, S. Trebst, M. Troyer, M. L. Wall, P. Werner and S. Wessel, Journal of Statistical Mechanics: Theory and Experiment 2011, P05001 (2011). [2] F. Alet, P. Dayal, A. Grzesik, A. Honecker, M. Korner, A. Lauchli, S. R. Manmana, I. P. McCulloch, F.xa0Michel, R. M. Noack, G. Schmid, U. Schollwock, F. Stockli, S. Todo, S. Trebst, M. Troyer, P. Werner, S. Wessel, J. Phys. Soc. Japan 74S (2005) 30. [3] A. Albuquerque, F. Alet, P. Corboz, P. Dayal, A. Feiguin, S. Fuchs, L. Gamper, E. Gull, S. Gurtler, A.xa0Honecker, R. Igarashi, M. Korner, A. Kozhevnikov, A. Lauchli, S. Manmana, M. Matsumoto, I.xa0McCulloch, F. Michel, R. Noack, G. Pawlowski, L. Pollet, T. Pruschke, U. Schollwock, S. Todo, S. Trebst, M. Troyer, P.xa0Werner and S. Wessel, J. Magn. Magn. Mater. 310, 1187 (2007), proceedings of the 17th International Conference on Magnetism The International Conference on Magnetism. [4] C. L. Lawson, R. J. Hanson, D. R. Kincaid, and F. T. Krogh, ACM Transactions on Mathematical Software 5, 324 (1979). [5] L. S. Blackford, J. Demmel, I. Du, G. Henry, M. Heroux, L. Kaufman, A. Lumsdaine, A. Petitet, and R.xa0C.xa0Whaley, ACM Trans. Math. Softw. 28, 135 (2002). [6] E. Anderson, Z. Bai, C. Bischof, S. Blackford, J. Demmel, J. Dongarra, J. Du Croz, A. Greenbaum, S.xa0Hammarling, A. McKenney, and D. Sorensen, LAPACK Users’ Guide, 3rd ed. (Society for Industrial and Applied Mathematics, Philadelphia, PA, 1999). [7] The HDF Group, Hierarchical data format version 5, http://www.hdfgroup.org/HDF5 (2000–2010). [8] E. Gull, A. J. Millis, A. I. Lichtenstein, A. N. Rubtsov, M. Troyer and P. Werner, Rev. Mod. Phys. 83, 349 (2011).

Dual Fermion Approach to Susceptibility of Correlated Lattice Fermions

In this paper, we show how the two-particle Green function (2PGF) can be obtained within the framework of the Dual Fermion approach. This facilitates the calculation of the susceptibility in strongly correlated systems where long-ranged non-local correlations cannot be neglected. We formulate the Bethe-Salpeter equations for the full vertex in the particle-particle and particle-hole channels and introduce an approximation for practical calculations. The scheme is applied to the two-dimensional Hubbard model at half filling. The spin-spin susceptibility is found to strongly increase for the wavevector q = (�,�), indicating the antiferromagnetic instability. We find a suppression of the critical temperature compared to the mean-field result due to the incorporation of the non-local spin-fluctuations.

Improved Estimators for the Self-Energy and Vertex Function in Hybridization Expansion Continuous-Time Quantum Monte Carlo Simulations

measurement procedure for the self-energy and vertex function of the Anderson impurity model within the hybridization expansion continuous-time quantum Monte Carlo algorithm. The method is based on the measurement of higher-order correlation functions related to the quantities being sought through the equation of motion. For the case of interactions of densitydensity type, the additional correlators can be obtained at essentially no additional computational cost. In combination with a recently introduced method for ltering the Monte Carlo noise using a representation in terms of orthogonal polynomials, we obtain data with unprecedented accuracy. This leads to an enhanced stability in analytical continuations of the self-energy or in two-particle based theories such as the dual fermion approach. As an illustration of the method we reexamine the previously reported spin-freezing and high-spin to low-spin transitions in a two-orbital model with density-density interactions. In both cases, the vertex function undergoes signicant changes, which suggests signicant corrections to the dynamical mean-eld solutions in dual fermion calculations.

Mechanisms of finite-temperature magnetism in the three-dimensional Hubbard model

We examine the nature of the transition to the antiferromagnetically ordered state in the half-filled three-dimensional Hubbard model using the dual-fermion multiscale approach. Consistent with analytics, in the weak-coupling regime we find that spin-flip excitations across the Fermi surface are important, and that the strong coupling regime is described by Heisenberg physics. In the intermediate interaction, strong correlation regime we find aspects of both local and non-local correlations. We analyze the critical exponents of the transition in the strong coupling regime and find them to be consistent with Heisenberg physics down to an interaction of