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:cond-mat/0612172 · Physical Review B

Quantum Monte Carlo impurity solver for cluster dynamical mean-field theory and electronic structure calculations with adjustable cluster base

Kristjan Haule

Published Apr 23, 2007• 613 citations
arXivHTMLPDF

Abstract

We generalized the recently introduced impurity solver [P. Werner et al., Phys. Rev. Lett. 97, 076405 (2006)] based on the diagrammatic expansion around the atomic limit and quantum Monte Carlo summation of the diagrams. We present generalization to the cluster of impurities, which is at the heart of the cluster dynamical mean-field methods, and to realistic multiplet structure of a correlated atom, which will allow a high-precision study of actinide and lanthanide based compounds with the combination of the dynamical mean-field theory and band-structure methods. The approach is applied to both the two-dimensional Hubbard and $t\text{\ensuremath{-}}J$ models within cellular dynamical mean-field method. The efficient implementation of the algorithm, which we describe in detail, allows us to study coherence of the system at low temperature from the underdoped to overdoped regime. We show that the point of maximal superconducting transition temperature coincides with the point of maximum scattering rate, although this optimal doped point appears at different electron densities in the two models. The power of the method is further demonstrated in the example of the Kondo volume collapse transition in cerium. The valence histogram of the dynamical mean-field theory solution is presented, showing the importance of the multiplet splitting of the atomic states.

DOI
10.1103/physrevb.75.155113

Related papers

Dynamical variational Monte Carlo as a quantum impurity solver: Application to cluster dynamical mean field theory

Peter Rosenberg, David Sénéchal, A.–M. S. Tremblay +1 more

Two of the primary sources of error in the Cluster dynamical mean field theory (CDMFT) technique arise from the use of finite size clusters and finite size baths, which makes the development of impurity solvers that can treat larger systems an essential goal. In this work we introduce an impurity solver based on the recently developed dynamical variational Monte Carlo (dVMC) method. Variational Monte Carlo possesses a favorable scaling as a function of system size, which enables the treatment of systems beyond the reach of current exact diagonalization solvers. To benchmark the technique, we perform a systematic set of CDMFT calculations on the single-band one-dimensional Hubbard model. We compare to results obtained with an exact diagonalization solver for small clusters, and against the exact solution in the thermodynamic limit obtained by Lieb and Wu for larger clusters. The development of improved impurity solvers will help extend the reach of quantum cluster methods, which can be applied to a wide range of strongly correlated electron systems, promising new insights into their emergent behavior.

arXiv:2307.15738•Physical review. B./Physical review. B•Dec 8, 2023•4 citations

Quantum Monte Carlo impurity solvers for multi-orbital problems and frequency-dependent interactions

Hiroshi Shinaoka, Fakher F. Assaad, Nils Blümer +1 more

The solution of an auxiliary quantum impurity system is the computationally expensive step in dynamical mean field theory simulations of lattice models and materials. In this review, we discuss Monte Carlo based impurity solvers, which are suitable for a wide range of applications. In particular, we present an efficient implementation of the hybridization expansion approach, which enables the simulation of multiorbital impurity problems with off-diagonal and complex hybridizations, and dynamically screened (retarded) density-density interactions. As a complementary approach, we discuss an impurity solver based on the determinant Monte Carlo method, which scales favorably with inverse temperature and hence provides access to the very low temperature regime. The usefulness of these state-of-the-art impurity solvers is demonstrated with applications to the downfolding problem, i.e., the systematic derivation of dynamically screened interactions for low-energy effective models, and to pyrochlore iridates, where the spin-orbit coupling leads to complex hybridization functions in a multi-orbital system. As a benchmark for cluster extensions of dynamical mean field theory, we also present results from lattice Monte Carlo simulations for the momentum dependence of the pseudo-gap in the half-filled two-dimensional Hubbard model.

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

Continuous-time impurity solvers

James E. Gubernatis, Naoki Kawashima, Philipp Werner

This chapter, like the previous one, describes quantum Monte Carlo methods for the simulation of a system of interacting electrons at nonzero temperature. Here, we target a particular class of problems, characterized by a small number of correlated sites or orbitals, that we call impurity problems. While we could apply the methods of the previous chapter to such impurity models, the algorithms discussed here are more advantageous not only because they work in continuous imaginary time and hence lack the Trotter approximation error inherent to the previous methods, but also because they manipulate determinants of matrices of smaller size, which can be handled more efficiently. At the same time, the continuous-time approach is applicable to broader classes of impurity problems. Driving the development of these continuous-time impurity solvers was the desire to perform more efficient simulations of correlated lattice models within the framework of dynamical mean-field theory. This formalism maps the lattice problem onto an impurity problem, whose parameters are determined by a self-consistency condition. In some cases, the self-consistent solution gives a good approximation of the properties of the original lattice problem. When we combine this dynamical mean-field approximation with the local-density approximation for electronic structure calculations, we obtain a powerful scheme for electronic structure calculations of strongly correlated materials. In this chapter we discuss the most important classes of impurity models, sketch the basics of the dynamical mean-field approximation, and detail the different variants of continuous-time impurity solvers. Quantum impurity models A quantum impurity model describes an atom or molecule embedded in some host with which it exchanges electrons, spin, and energy. This exchange allows the impurity to make transitions between different quantum states, and in the presence of interactions within or between the impurity orbitals, these transitions lead to nontrivial dynamical properties. Quantum impurity models play a prominent role, for example, in the theoretical description of the magnetic and electric properties of dilute metal alloys and in theoretical studies of quantum dots and molecular conductors. These models also appear as an auxiliary problem whose solution yields the dynamical mean-field description of correlated lattice models.

Cambridge University Press eBooks•Jan 31, 2016•1 citation

Continuous-time quantum Monte Carlo solver for dynamical mean field theory in the compact Legendre representation

Evan Sheridan, Cédric Weber, Evgeny A. Plekhanov +1 more

Dynamical mean-field theory (DMFT) is one of the most widely-used methods to treat accurately electron correlation effects in ab-initio real material calculations. Many modern large-scale implementations of DMFT in electronic structure codes involve solving a quantum impurity model with a continuous-time quantum Monte Carlo (CT-QMC) solver [Rubtsov et al., Phys. Rev. B 72, 035122 (2005); Werner et al., Phys. Rev. Lett. 97, 076405 (2006); Werner and Millis, Phys. Rev. B 74, 155107 (2006); Gull et al., Rev. Mod. Phys. 83, 349 (2011)]. The main advantage of CT-QMC is that, unlike standard quantum Monte Carlo approaches, it is able to generate the local Green's functions $G(\ensuremath{\tau})$ of the correlated system on an arbitrarily fine imaginary time $\ensuremath{\tau}$ grid, and is free of any systematic errors. In this work, we extend a hybrid QMC solver proposed by Khatami et al. [Phys. Rev. E 81, 056703 (2010)] and Rost et al. [Phys. Rev. E 87, 053305 (2013)] to a multiorbital context. This has the advantage of enabling impurity solver QMC calculations to scale linearly with inverse temperature $\ensuremath{\beta}$, and permit its application to $d$- and $f$-band materials. In addition, we present a Green's-function processing scheme which generates accurate quasicontinuous imaginary time solutions of the impurity problem which overcome errors inherent to standard QMC approaches. This solver and processing scheme are incorporated into a full DFT+DMFT calculation using the castep DFT code [Clark et al., Z. Kristallogr. 220, 567 (2005)]. Benchmark calculations for ${\mathrm{SrVO}}_{3}$ properties are presented. The computational efficiency of this method is also demonstrated.

arXiv:1811.12869•Physical review. B./Physical review. B•May 30, 2019•14 citations

Strong-coupling solver for the quantum impurity model

Xi Dai, Kristjan Haule, Gabriel Kotliar

We propose a fast impurity solver for the general quantum impurity model based on the perturbation theory around the atomic limit, which can be used in combination with the local density approximation (LDA) and the dynamical mean-field theory (DMFT). We benchmark the solver in the two-band Hubbard model within DMFT against quantum Monte Carlo (QMC) and numerical renormalization-group (NRG) results. We find that the solver works very well in the paramagnetic Mott insulator phase. We also apply this impurity solver to the DMFT study of the antiferromagnetic phase transition in the unfrustrated Bethe lattice. The Neel temperature obtained by the fast impurity solver agrees very well with the QMC results in the large Hubbard U limit. The method is a promising tool to be used in combination with the $\mathrm{LDA}+\mathrm{DMFT}$ to study Mott insulators starting from first principles.

arXiv:cond-mat/0412505•Physical Review B•Jul 8, 2005•26 citations

Cellular Dynamical Mean-field Theory in large Impurity Clusters

Marcel Klett

Even after decades of intense research, the single band Hubbard model representing the fundamental model for interacting quantum systems and electronic correlations remains an unsolved cornerstone paradigm in theoretical solid state physics. Within the non-perturbative dynamical mean-field theory (DMFT), the lattice problem is mapped to a self-consistent auxiliary quantum impurity model leading to a local approximation of the self-energy. Nonlocal correlations can be included by the real-space cluster extension of the DMFT (CDMFT). Enlarging the unit cell makes the calculation of a numerically exact solution in CDMFT very challenging. We here use a quantum Monte-Carlo approach in the imaginary-time space in order to solve the self-consistency equations for different cluster sizes. In order to restore the broken translational symmetry re-periodization schemes for the Green function, the self-energy, or its cumulants have been introduced. However, these suffer from ambiguity and may even lead to convergence problems when attempted inside the self-consistency loop. The comparison to numerically exact diagrammatic quantum Monte Carlo calculations shows that introducing a so-called centerfocused extrapolation (CFE) to approximate the lattice self-energy yields very accurate results. Moreover, the CFE converges faster with the cluster size than previous periodization schemes. We here perform a detailed CDMFT analysis for the single-band Hubbard model on the 2D square lattice, reaching real-space cluster sizes of up to 9x9 sites. In addition to spectral properties, we also compute two-particle correlation functions which are not accessible in DMFT. Using benchmarks against diagrammatic Monte Carlo at high temperature, we show that the cluster spin susceptibility can be extrapolated to the exact result at large cluster size. In particular, the exponential decay of spin-spin correlations is very well captured by CDMFT calculations, even when correlations extend beyond the size of the cluster. We further present results at lower temperature T and larger $U$ than the range currently accessible with diagrammatic Monte Carlo methods techniques. The CDMFT+CFE represents therefore a powerful computational tool to access the physics of non-local correlations beyond dynamical mean-field theory.

•Dec 29, 2020•0 citations
We generalized the recently introduced impurity solver [P. Werner et al., Phys. Rev. Lett. 97, 076405 (2006)] based on the diagrammatic expansion around the atomic limit and quantum Monte Carlo summation of the diagrams. We present generalization to the cluster of impurities, which is at the heart of the cluster dynamical mean-field methods, and to realistic multiplet structure of a correlated atom, which will allow a high-precision study of actinide and lanthanide based compounds with the combination of the dynamical mean-field theory and band-structure methods. The approach is applied to both the two-dimensional Hubbard and t\ensuremath-Jt\text{\ensuremath{-}}Jt\ensuremath-J models within cellular dynamical mean-field method. The efficient implementation of the algorithm, which we describe in detail, allows us to study coherence of the system at low temperature from the underdoped to overdoped regime. We show that the point of maximal superconducting transition temperature coincides with the point of maximum scattering rate, although this optimal doped point appears at different electron densities in the two models. The power of the method is further demonstrated in the example of the Kondo volume collapse transition in cerium. The valence histogram of the dynamical mean-field theory solution is presented, showing the importance of the multiplet splitting of the atomic states.
Two of the primary sources of error in the Cluster dynamical mean field theory (CDMFT) technique arise from the use of finite size clusters and finite size baths, which makes the development of impurity solvers that can treat larger systems an essential goal. In this work we introduce an impurity solver based on the recently developed dynamical variational Monte Carlo (dVMC) method. Variational Monte Carlo possesses a favorable scaling as a function of system size, which enables the treatment of systems beyond the reach of current exact diagonalization solvers. To benchmark the technique, we perform a systematic set of CDMFT calculations on the single-band one-dimensional Hubbard model. We compare to results obtained with an exact diagonalization solver for small clusters, and against the exact solution in the thermodynamic limit obtained by Lieb and Wu for larger clusters. The development of improved impurity solvers will help extend the reach of quantum cluster methods, which can be applied to a wide range of strongly correlated electron systems, promising new insights into their emergent behavior.
The solution of an auxiliary quantum impurity system is the computationally expensive step in dynamical mean field theory simulations of lattice models and materials. In this review, we discuss Monte Carlo based impurity solvers, which are suitable for a wide range of applications. In particular, we present an efficient implementation of the hybridization expansion approach, which enables the simulation of multiorbital impurity problems with off-diagonal and complex hybridizations, and dynamically screened (retarded) density-density interactions. As a complementary approach, we discuss an impurity solver based on the determinant Monte Carlo method, which scales favorably with inverse temperature and hence provides access to the very low temperature regime. The usefulness of these state-of-the-art impurity solvers is demonstrated with applications to the downfolding problem, i.e., the systematic derivation of dynamically screened interactions for low-energy effective models, and to pyrochlore iridates, where the spin-orbit coupling leads to complex hybridization functions in a multi-orbital system. As a benchmark for cluster extensions of dynamical mean field theory, we also present results from lattice Monte Carlo simulations for the momentum dependence of the pseudo-gap in the half-filled two-dimensional Hubbard model.
This chapter, like the previous one, describes quantum Monte Carlo methods for the simulation of a system of interacting electrons at nonzero temperature. Here, we target a particular class of problems, characterized by a small number of correlated sites or orbitals, that we call impurity problems. While we could apply the methods of the previous chapter to such impurity models, the algorithms discussed here are more advantageous not only because they work in continuous imaginary time and hence lack the Trotter approximation error inherent to the previous methods, but also because they manipulate determinants of matrices of smaller size, which can be handled more efficiently. At the same time, the continuous-time approach is applicable to broader classes of impurity problems. Driving the development of these continuous-time impurity solvers was the desire to perform more efficient simulations of correlated lattice models within the framework of dynamical mean-field theory. This formalism maps the lattice problem onto an impurity problem, whose parameters are determined by a self-consistency condition. In some cases, the self-consistent solution gives a good approximation of the properties of the original lattice problem. When we combine this dynamical mean-field approximation with the local-density approximation for electronic structure calculations, we obtain a powerful scheme for electronic structure calculations of strongly correlated materials. In this chapter we discuss the most important classes of impurity models, sketch the basics of the dynamical mean-field approximation, and detail the different variants of continuous-time impurity solvers. Quantum impurity models A quantum impurity model describes an atom or molecule embedded in some host with which it exchanges electrons, spin, and energy. This exchange allows the impurity to make transitions between different quantum states, and in the presence of interactions within or between the impurity orbitals, these transitions lead to nontrivial dynamical properties. Quantum impurity models play a prominent role, for example, in the theoretical description of the magnetic and electric properties of dilute metal alloys and in theoretical studies of quantum dots and molecular conductors. These models also appear as an auxiliary problem whose solution yields the dynamical mean-field description of correlated lattice models.
Dynamical mean-field theory (DMFT) is one of the most widely-used methods to treat accurately electron correlation effects in ab-initio real material calculations. Many modern large-scale implementations of DMFT in electronic structure codes involve solving a quantum impurity model with a continuous-time quantum Monte Carlo (CT-QMC) solver [Rubtsov et al., Phys. Rev. B 72, 035122 (2005); Werner et al., Phys. Rev. Lett. 97, 076405 (2006); Werner and Millis, Phys. Rev. B 74, 155107 (2006); Gull et al., Rev. Mod. Phys. 83, 349 (2011)]. The main advantage of CT-QMC is that, unlike standard quantum Monte Carlo approaches, it is able to generate the local Green's functions $G(\ensuremath{\tau})$ of the correlated system on an arbitrarily fine imaginary time $\ensuremath{\tau}$ grid, and is free of any systematic errors. In this work, we extend a hybrid QMC solver proposed by Khatami et al. [Phys. Rev. E 81, 056703 (2010)] and Rost et al. [Phys. Rev. E 87, 053305 (2013)] to a multiorbital context. This has the advantage of enabling impurity solver QMC calculations to scale linearly with inverse temperature $\ensuremath{\beta}$, and permit its application to ddd- and fff-band materials. In addition, we present a Green's-function processing scheme which generates accurate quasicontinuous imaginary time solutions of the impurity problem which overcome errors inherent to standard QMC approaches. This solver and processing scheme are incorporated into a full DFT+DMFT calculation using the castep DFT code [Clark et al., Z. Kristallogr. 220, 567 (2005)]. Benchmark calculations for SrVO3{\mathrm{SrVO}}_{3}SrVO3​ properties are presented. The computational efficiency of this method is also demonstrated.
We propose a fast impurity solver for the general quantum impurity model based on the perturbation theory around the atomic limit, which can be used in combination with the local density approximation (LDA) and the dynamical mean-field theory (DMFT). We benchmark the solver in the two-band Hubbard model within DMFT against quantum Monte Carlo (QMC) and numerical renormalization-group (NRG) results. We find that the solver works very well in the paramagnetic Mott insulator phase. We also apply this impurity solver to the DMFT study of the antiferromagnetic phase transition in the unfrustrated Bethe lattice. The Neel temperature obtained by the fast impurity solver agrees very well with the QMC results in the large Hubbard U limit. The method is a promising tool to be used in combination with the LDA+DMFT\mathrm{LDA}+\mathrm{DMFT}LDA+DMFT to study Mott insulators starting from first principles.
Even after decades of intense research, the single band Hubbard model representing the fundamental model for interacting quantum systems and electronic correlations remains an unsolved cornerstone paradigm in theoretical solid state physics. Within the non-perturbative dynamical mean-field theory (DMFT), the lattice problem is mapped to a self-consistent auxiliary quantum impurity model leading to a local approximation of the self-energy. Nonlocal correlations can be included by the real-space cluster extension of the DMFT (CDMFT). Enlarging the unit cell makes the calculation of a numerically exact solution in CDMFT very challenging. We here use a quantum Monte-Carlo approach in the imaginary-time space in order to solve the self-consistency equations for different cluster sizes. In order to restore the broken translational symmetry re-periodization schemes for the Green function, the self-energy, or its cumulants have been introduced. However, these suffer from ambiguity and may even lead to convergence problems when attempted inside the self-consistency loop. The comparison to numerically exact diagrammatic quantum Monte Carlo calculations shows that introducing a so-called centerfocused extrapolation (CFE) to approximate the lattice self-energy yields very accurate results. Moreover, the CFE converges faster with the cluster size than previous periodization schemes. We here perform a detailed CDMFT analysis for the single-band Hubbard model on the 2D square lattice, reaching real-space cluster sizes of up to 9x9 sites. In addition to spectral properties, we also compute two-particle correlation functions which are not accessible in DMFT. Using benchmarks against diagrammatic Monte Carlo at high temperature, we show that the cluster spin susceptibility can be extrapolated to the exact result at large cluster size. In particular, the exponential decay of spin-spin correlations is very well captured by CDMFT calculations, even when correlations extend beyond the size of the cluster. We further present results at lower temperature T and larger UUU than the range currently accessible with diagrammatic Monte Carlo methods techniques. The CDMFT+CFE represents therefore a powerful computational tool to access the physics of non-local correlations beyond dynamical mean-field theory.