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:0708.0718 · Journal of the Physical Society of Japan

Continuous-Time Quantum Monte Carlo Method for the Coqblin–Schrieffer Model

Junya Otsuki, Hiroaki Kusunose, Philipp Werner, Yoshio Kuramoto

Published Oct 24, 2007• 68 citations
arXivHTMLPDF

Abstract

An impurity solver based on a continuous-time quantum Monte Carlo method is developed for the Coqblin-Schrieffer model. The Monte Carlo simulation does not encounter a sign problem for antiferromagnetic interactions, and accurately reproduces the Kondo effect. Our algorithm can deal with an arbitrary number N of local degrees of freedom, becomes more efficient for larger values of N, and is hence suitable for models with orbital degeneracy. The dynamical susceptibility and the impurity t-matrix are derived with the aid of the Padé approximation for various values of N, and good agreement is found with other methods and available exact results. We point out that the Korringa-Shiba relation needs correction for a finite value of the exchange interaction.

DOI
10.1143/jpsj.76.114707

Related papers

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

Hybridization expansion impurity solver: General formulation and application to Kondo lattice and two-orbital models

Philipp Werner, Andrew J. Millis

A recently developed continuous time solver based on an expansion in hybridization about an exactly solved local limit is reformulated in a manner appropriate for general classes of quantum impurity models including spin exchange and pair hopping terms. The utility of the approach is demonstrated via applications to the dynamical mean field theory of the Kondo lattice and two-orbital models. The algorithm can handle low temperatures and strong couplings without encountering a sign problem.

arXiv:cond-mat/0607136•Physical Review B•Oct 11, 2006•438 citations

Superperturbation solver for quantum impurity models

Hartmut Hafermann, C. Jung, S. Brener +3 more

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.

arXiv:0809.1051•Europhysics Letters (EPL)•Jan 1, 2009•51 citations

Superperturbation solver for quantum impurity models

Hartmut Hafermann, Jung, C., S. Brener +3 more

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.

•Jan 1, 2009•36 citations

Multilevel algorithm for quantum-impurity models

Jaebeom Yoo, Shailesh Chandrasekharan A, Harold U. Baranger

A continuous-time path integral quantum Monte Carlo method using the directed-loop algorithm is developed to simulate the Anderson single-impurity model in the occupation number basis. Although the method suffers from a sign problem at low temperatures, the new algorithm has many advantages over conventional algorithms. For example, the model can be easily simulated in the Kondo limit without time discretization errors. Furthermore, many observables including the impurity susceptibility and a variety of fermionic observables can be calculated efficiently. Finally the new approach allows us to explore a general technique, called the multilevel algorithm, to solve the sign problem. We find that the multilevel algorithm is able to generate an exponentially large number of configurations with an effort that grows as a polynomial in inverse temperature such that configurations with a positive sign dominate over those with negative signs. Our algorithm can be easily generalized to other multi-impurity problems.

arXiv:cond-mat/0408123•Physical Review E•Mar 31, 2005•3 citations

Steady-state properties of multi-orbital systems using quantum Monte Carlo

André Erpenbeck, Thomas Blommel, Lei Zhang +3 more

A precise dynamical characterization of quantum impurity models with multiple interacting orbitals is challenging. In quantum Monte Carlo methods, this is embodied by sign problems. A dynamical sign problem makes it exponentially difficult to simulate long times. A multi-orbital sign problem generally results in a prohibitive computational cost for systems with multiple impurity degrees of freedom even in static equilibrium calculations. Here, we present a numerically exact inchworm method that simultaneously alleviates both sign problems, enabling simulation of multi-orbital systems directly in the equilibrium or nonequilibrium steady-state. The method combines ideas from the recently developed steady-state inchworm Monte Carlo framework [Erpenbeck et al., Phys. Rev. Lett. 130, 186301 (2023)] with other ideas from the equilibrium multi-orbital inchworm algorithm [Eidelstein et al., Phys. Rev. Lett. 124, 206405 (2020)]. We verify our method by comparison with analytical limits and numerical results from previous methods.

arXiv:2407.00771•The Journal of Chemical Physics•Sep 4, 2024•12 citations
An impurity solver based on a continuous-time quantum Monte Carlo method is developed for the Coqblin-Schrieffer model. The Monte Carlo simulation does not encounter a sign problem for antiferromagnetic interactions, and accurately reproduces the Kondo effect. Our algorithm can deal with an arbitrary number N of local degrees of freedom, becomes more efficient for larger values of N, and is hence suitable for models with orbital degeneracy. The dynamical susceptibility and the impurity t-matrix are derived with the aid of the Padé approximation for various values of N, and good agreement is found with other methods and available exact results. We point out that the Korringa-Shiba relation needs correction for a finite value of the exchange interaction.
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.
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.
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.
A continuous-time path integral quantum Monte Carlo method using the directed-loop algorithm is developed to simulate the Anderson single-impurity model in the occupation number basis. Although the method suffers from a sign problem at low temperatures, the new algorithm has many advantages over conventional algorithms. For example, the model can be easily simulated in the Kondo limit without time discretization errors. Furthermore, many observables including the impurity susceptibility and a variety of fermionic observables can be calculated efficiently. Finally the new approach allows us to explore a general technique, called the multilevel algorithm, to solve the sign problem. We find that the multilevel algorithm is able to generate an exponentially large number of configurations with an effort that grows as a polynomial in inverse temperature such that configurations with a positive sign dominate over those with negative signs. Our algorithm can be easily generalized to other multi-impurity problems.
A precise dynamical characterization of quantum impurity models with multiple interacting orbitals is challenging. In quantum Monte Carlo methods, this is embodied by sign problems. A dynamical sign problem makes it exponentially difficult to simulate long times. A multi-orbital sign problem generally results in a prohibitive computational cost for systems with multiple impurity degrees of freedom even in static equilibrium calculations. Here, we present a numerically exact inchworm method that simultaneously alleviates both sign problems, enabling simulation of multi-orbital systems directly in the equilibrium or nonequilibrium steady-state. The method combines ideas from the recently developed steady-state inchworm Monte Carlo framework [Erpenbeck et al., Phys. Rev. Lett. 130, 186301 (2023)] with other ideas from the equilibrium multi-orbital inchworm algorithm [Eidelstein et al., Phys. Rev. Lett. 124, 206405 (2020)]. We verify our method by comparison with analytical limits and numerical results from previous methods.