Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:1402.0807 · Physical Review B

Efficient real-frequency solver for dynamical mean-field theory

Yi Lu, M. Höppner, O. Gunnarsson, M. W. Haverkort

Published Aug 5, 2014• 221 citations
arXivHTMLPDF

Abstract

We here present how a self-consistent solution of the dynamical mean field theory equations can be obtained using exact diagonalization of an Anderson impurity model with accuracies comparable to those found using renormalization group or quantum Monte Carlo methods. We show how one can solve a correlated quantum impurity coupled to several hundred uncorrelated bath sites, using a restricted active basis set. The number of bath sites determines the resolution of the obtained spectral function, which consists of peaks roughly spaced by the band width divided by the number of bath sites. The self-consistency cycles are fully performed on the real frequency axis and expressed as numerical stable matrix operations. The same impurity solver has been used on Ligand Field and finite size cluster calculations and is capable of treating involved Hamiltonians including the full rotational invariant Coulomb interaction, spin-orbit coupling and low symmetry crystal-fields. The proposed method allows for the calculation of a variety of correlation functions at little extra cost.

DOI
10.1103/physrevb.90.085102

Related papers

Exact diagonalization as an impurity solver in dynamical mean field theory

Yi Lu, M. W. Haverkort

The dynamical mean-field theory (DMFT) maps a correlated lattice problem onto an impurity problem of a single correlated site coupled to an uncorrelated bath. Most implementations solve the DMFT equations using quantum Monte-Carlo sampling on the imaginary time and frequency (Matsubara) axis. We will here review alternative methods using exact diagonalization, i.e., representing the many-body ground state of the impurity as a sum over Slater determinants and calculating Green’s functions using iterative Lanczos procedures. The advantage being that these methods have no sign problem, can handle involved multi-orbital Hamiltonians (low crystal symmetry, spin-orbit coupling) and – when working completely on the real axis – do not need a mathematically ill-posed analytical continuation. The disadvantage of traditional implementations of exact diagonalization has been the exponential scaling of the calculation problem as a function of number of bath discretization points. In the last part we will review how recent advances in exact diagonalization can evade the exponential barrier thereby increasing the number of bath discretization points to reach the thermodynamic limit.

The European Physical Journal Special Topics•Jul 1, 2017•17 citations

Solving the Multi-site and Multi-orbital Dynamical Mean Field Theory Using Density Matrix Renormalization

Y. Núñez-Fernández, K. Hallberg

We implement an efficient numerical method to calculate response functions of complex impurities based on the Density Matrix Renormalization Group (DMRG) and use it as the impurity-solver of the Dynamical Mean Field Theory (DMFT). This method uses the correction vector to obtain precise Green's functions on the real frequency axis at zero temperature. By using a self-consistent bath configuration with very low entanglement, we take full advantage of the DMRG to calculate dynamical response functions paving the way to treat large effective impurities such as those corresponding to multi-orbital interacting models and multi-site or multi-momenta clusters. This method leads to reliable calculations of non-local self energies at arbitrary dopings and interactions and at any energy scale.

Frontiers in Physics•Feb 26, 2018•9 citations

Solving the multi-site and multi-orbital Dynamical Mean Field Theory using Density Matrix Renormalization

Y. Núñez-Fernández, K. Hallberg

We implement an efficient numerical method to calculate response functions of complex impurities based on the Density Matrix Renormalization Group (DMRG) and use it as the impurity-solver of the Dynamical Mean Field Theory (DMFT). This method uses the correction vector to obtain precise Green's functions on the real frequency axis at zero temperature. By using a self-consistent bath configuration with very low entanglement, we take full advantage of the DMRG to calculate dynamical response functions paving the way to treat large effective impurities such as those corresponding to multi-orbital interacting models and multi-site or multi-momenta clusters. This method leads to reliable calculations of non-local self energies at arbitrary dopings and interactions and at any energy scale.

arXiv:1711.08745•Nov 23, 2017•1 citation

Natural-orbital impurity solver and projection approach for Green's functions

Yi Lu, Xiaodong Cao, P. Hansmann +1 more

We extend a previously proposed rotation and truncation scheme to optimize quantum Anderson impurity calculations with exact diagonalization [Y. Lu, M. H\"oppner, O. Gunnarsson, and M. W. Haverkort, Phys. Rev. B 90, 085102 (2014)] to density-matrix renormalization group (DMRG) calculations. The method reduces the solution of a full impurity problem with virtually unlimited bath sites to that of a small subsystem based on a natural impurity orbital basis set. The later is solved by DMRG in combination with a restricted-active-space truncation scheme. The method allows one to compute Green's functions directly on the real frequency or time axis. We critically test the convergence of the truncation scheme using a one-band Hubbard model solved in the dynamical mean-field theory. The projection is exact in the limit of both infinitely large and small Coulomb interactions. For all parameter ranges, the accuracy of the projected solution converges exponentially to the exact solution with increasing subsystem size.

arXiv:1909.02757•Physical review. B./Physical review. B•Sep 16, 2019•36 citations

Distributional exact diagonalization formalism for quantum impurity models

Mats Granath, Hugo U. R. Strand

We develop a method for calculating the self-energy of a quantum impurity coupled to a continuous bath by stochastically generating a distribution of finite Anderson models that are solved by exact diagonalization, using the noninteracting local spectral function as a probability distribution for the sampling. The method enables calculation of the full analytic self-energy and single-particle Green's function in the complex frequency plane, without analytic continuation, and can be used for both finite and zero temperature at arbitrary fillings. Results are in good agreement with imaginary frequency data from continuous-time quantum Monte Carlo calculations for the single-impurity Anderson model and the two-orbital Hubbard model within dynamical mean-field theory (DMFT) as well as real frequency data for self-energy of the single-band Hubbard model within DMFT using the numerical renormalization group. The method should be applicable to a wide range of quantum impurity models and particularly useful when high-precision real frequency results are sought.

arXiv:1201.6160•Physical Review B•Sep 10, 2012•36 citations

Coupled-cluster impurity solvers for dynamical mean-field theory

Tianyu Zhu, Carlos A. Jiménez-Hoyos, James D. McClain +2 more

We describe the use of coupled-cluster theory as an impurity solver in dynamical mean-field theory (DMFT) and its cluster extensions. We present numerical results at the level of coupled-cluster theory with single and double excitations (CCSD) for the density of states and self-energies of cluster impurity problems in the one- and two-dimensional Hubbard models. Comparison to exact diagonalization shows that CCSD produces accurate density of states and self-energies at a variety of values of $U/t$ and filling fractions. However, the low cost allows for the use of many bath sites, which we define by a discretization of the hybridization directly on the real frequency axis. We observe convergence of dynamical quantities using approximately 30 bath sites per impurity site, with our largest 4-site cluster DMFT calculation using 120 bath sites. We suggest that coupled-cluster impurity solvers will be attractive in ab initio formulations of dynamical mean-field theory.

arXiv:1905.12050•Physical review. B./Physical review. B•Sep 26, 2019•57 citations
We here present how a self-consistent solution of the dynamical mean field theory equations can be obtained using exact diagonalization of an Anderson impurity model with accuracies comparable to those found using renormalization group or quantum Monte Carlo methods. We show how one can solve a correlated quantum impurity coupled to several hundred uncorrelated bath sites, using a restricted active basis set. The number of bath sites determines the resolution of the obtained spectral function, which consists of peaks roughly spaced by the band width divided by the number of bath sites. The self-consistency cycles are fully performed on the real frequency axis and expressed as numerical stable matrix operations. The same impurity solver has been used on Ligand Field and finite size cluster calculations and is capable of treating involved Hamiltonians including the full rotational invariant Coulomb interaction, spin-orbit coupling and low symmetry crystal-fields. The proposed method allows for the calculation of a variety of correlation functions at little extra cost.
The dynamical mean-field theory (DMFT) maps a correlated lattice problem onto an impurity problem of a single correlated site coupled to an uncorrelated bath. Most implementations solve the DMFT equations using quantum Monte-Carlo sampling on the imaginary time and frequency (Matsubara) axis. We will here review alternative methods using exact diagonalization, i.e., representing the many-body ground state of the impurity as a sum over Slater determinants and calculating Green’s functions using iterative Lanczos procedures. The advantage being that these methods have no sign problem, can handle involved multi-orbital Hamiltonians (low crystal symmetry, spin-orbit coupling) and – when working completely on the real axis – do not need a mathematically ill-posed analytical continuation. The disadvantage of traditional implementations of exact diagonalization has been the exponential scaling of the calculation problem as a function of number of bath discretization points. In the last part we will review how recent advances in exact diagonalization can evade the exponential barrier thereby increasing the number of bath discretization points to reach the thermodynamic limit.
We implement an efficient numerical method to calculate response functions of complex impurities based on the Density Matrix Renormalization Group (DMRG) and use it as the impurity-solver of the Dynamical Mean Field Theory (DMFT). This method uses the correction vector to obtain precise Green's functions on the real frequency axis at zero temperature. By using a self-consistent bath configuration with very low entanglement, we take full advantage of the DMRG to calculate dynamical response functions paving the way to treat large effective impurities such as those corresponding to multi-orbital interacting models and multi-site or multi-momenta clusters. This method leads to reliable calculations of non-local self energies at arbitrary dopings and interactions and at any energy scale.
We implement an efficient numerical method to calculate response functions of complex impurities based on the Density Matrix Renormalization Group (DMRG) and use it as the impurity-solver of the Dynamical Mean Field Theory (DMFT). This method uses the correction vector to obtain precise Green's functions on the real frequency axis at zero temperature. By using a self-consistent bath configuration with very low entanglement, we take full advantage of the DMRG to calculate dynamical response functions paving the way to treat large effective impurities such as those corresponding to multi-orbital interacting models and multi-site or multi-momenta clusters. This method leads to reliable calculations of non-local self energies at arbitrary dopings and interactions and at any energy scale.
We extend a previously proposed rotation and truncation scheme to optimize quantum Anderson impurity calculations with exact diagonalization [Y. Lu, M. H\"oppner, O. Gunnarsson, and M. W. Haverkort, Phys. Rev. B 90, 085102 (2014)] to density-matrix renormalization group (DMRG) calculations. The method reduces the solution of a full impurity problem with virtually unlimited bath sites to that of a small subsystem based on a natural impurity orbital basis set. The later is solved by DMRG in combination with a restricted-active-space truncation scheme. The method allows one to compute Green's functions directly on the real frequency or time axis. We critically test the convergence of the truncation scheme using a one-band Hubbard model solved in the dynamical mean-field theory. The projection is exact in the limit of both infinitely large and small Coulomb interactions. For all parameter ranges, the accuracy of the projected solution converges exponentially to the exact solution with increasing subsystem size.
We develop a method for calculating the self-energy of a quantum impurity coupled to a continuous bath by stochastically generating a distribution of finite Anderson models that are solved by exact diagonalization, using the noninteracting local spectral function as a probability distribution for the sampling. The method enables calculation of the full analytic self-energy and single-particle Green's function in the complex frequency plane, without analytic continuation, and can be used for both finite and zero temperature at arbitrary fillings. Results are in good agreement with imaginary frequency data from continuous-time quantum Monte Carlo calculations for the single-impurity Anderson model and the two-orbital Hubbard model within dynamical mean-field theory (DMFT) as well as real frequency data for self-energy of the single-band Hubbard model within DMFT using the numerical renormalization group. The method should be applicable to a wide range of quantum impurity models and particularly useful when high-precision real frequency results are sought.
We describe the use of coupled-cluster theory as an impurity solver in dynamical mean-field theory (DMFT) and its cluster extensions. We present numerical results at the level of coupled-cluster theory with single and double excitations (CCSD) for the density of states and self-energies of cluster impurity problems in the one- and two-dimensional Hubbard models. Comparison to exact diagonalization shows that CCSD produces accurate density of states and self-energies at a variety of values of U/tU/tU/t and filling fractions. However, the low cost allows for the use of many bath sites, which we define by a discretization of the hybridization directly on the real frequency axis. We observe convergence of dynamical quantities using approximately 30 bath sites per impurity site, with our largest 4-site cluster DMFT calculation using 120 bath sites. We suggest that coupled-cluster impurity solvers will be attractive in ab initio formulations of dynamical mean-field theory.