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:2201.03099 · Physical review. B./Physical review. B

Interaction-expansion inchworm Monte Carlo solver for lattice and impurity models

Jia Li, Yang Yu, Emanuel Gull, Guy Cohen

Published Apr 18, 2022• 12 citations
arXivHTMLPDF

Abstract

Multiorbital quantum impurity models with general interaction and hybridization terms appear in a wide range of applications, including embedding, quantum transport, and nanoscience. However, most quantum impurity solvers are restricted to a few impurity orbitals, discretized baths, diagonal hybridizations, or density-density interactions. Here, we generalize the inchworm quantum Monte Carlo method to the interaction expansion, and we explore its application to typical single- and multiorbital problems encountered in investigations of impurity and lattice models. Our implementation generically outperforms bare and bold-line quantum Monte Carlo algorithms in the interaction expansion. For the systems studied here, our implementation remains inferior to the more specialized hybridization expansion and auxiliary field algorithms. The problem of convergence to unphysical fixed points, which hampers so-called bold-line methods, is not encountered in inchworm Monte Carlo.

DOI
10.1103/physrevb.105.165133

Related papers

Multiorbital Quantum Impurity Solver for General Interactions and Hybridizations

Eitan Eidelstein, Emanuel Gull, Guy Cohen

We present a numerically exact inchworm Monte Carlo method for equilibrium multiorbital quantum impurity problems with general interactions and hybridizations. We show that the method, originally developed to overcome the dynamical sign problem in certain real-time propagation problems, can also overcome the sign problem as a function of temperature for equilibrium quantum impurity models. This is shown in several cases where the current method of choice, the continuous-time hybridization expansion, fails due to the sign problem. Our method therefore enables simulations of impurity problems as they appear in embedding theories without further approximations, such as the truncation of the hybridization or interaction structure or a discretization of the impurity bath with a set of discrete energy levels, and eliminates a crucial bottleneck in the simulation of ab initio embedding problems.

arXiv:1907.08570•Physical Review Letters•May 22, 2020•37 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

Inchworm quasi Monte Carlo for quantum impurities

Hugo U. R. Strand, Joseph Kleinhenz, Igor S. Krivenko

The inchworm expansion is a promising approach to solving strongly correlated quantum impurity models due to its reduction of the sign problem in real and imaginary time. However, inchworm Monte Carlo is computationally expensive, converging as 1 / N where N is the number of samples. We show that the imaginary-time integration is amenable to quasi Monte Carlo, with parametrically better 1 / N convergence, by mapping the Sobol low-discrepancy sequence from the hypercube to the simplex with the so-called Root transform. This extends the applicability of the inchworm method to, e.g., multiorbital Anderson impurity models with off-diagonal hybridization, relevant for materials simulation, where continuous-time hybridization expansion Monte Carlo has a severe sign problem. Published by the American Physical Society 2024

arXiv:2310.16957•Physical review. B./Physical review. B•Sep 24, 2024•6 citations

Inchworm tensor train hybridization expansion quantum impurity solver

Yang Yu, André Erpenbeck, Dominika Zgid +3 more

The investigation of quantum impurity models plays a crucial role in condensed matter physics because of their wide-ranging applications, such as embedding theories and transport problems. Traditional methods often fall short in producing accurate results for multi-orbital systems with complex interactions and off-diagonal hybridizations. Recently, tensor-train-based integration and summation techniques have shown promise as effective alternatives. In this study, we use tensor train methods to tackle quantum impurity problems formulated within the imaginary-time inchworm hybridization expansion framework. We identify key challenges in the inchworm expansion itself and its interplay with tensor-train-based methods. We demonstrate the accuracy and versatility of our approach by solving general quantum impurity problems. Our results suggest that tensor-train decomposition schemes offer a viable path toward accurate and efficient multi-orbital impurity solvers.

arXiv:2505.16117•Physical review. B./Physical review. B•Jul 28, 2025•1 citation

Double-expansion impurity solver for multiorbital models with dynamically screened U and J

Karim Steiner, Yusuke Nomura, Philipp Werner

We present a continuous-time Monte Carlo impurity solver for multiorbital impurity models which combines a strong-coupling hybridization expansion and a weak-coupling expansion in the Hund's coupling parameter $J$. This double-expansion approach allows to treat the dominant density-density interactions $U$ within the efficient segment representation. We test the approach for a two-orbital model with static interactions, and then explain how the double expansion allows to simulate models with frequency dependent $U(\ensuremath{\omega})$ and $J(\ensuremath{\omega})$. The method is used to investigate spin-state transitions in a toy model for fullerides, with repulsive bare $J$ but attractive screened $J$.

arXiv:1506.01173•Physical Review B•Sep 10, 2015•19 citations

Dynamical mean field theory of the bilayer Hubbard model with inchworm Monte Carlo

Dolev Goldberger, Yehonatan Fridman, Emanuel Gull +2 more

Dynamical mean field theory allows access to the physics of strongly correlated materials with nontrivial orbital structure, but relies on the ability to solve auxiliary multiorbital impurity problems. The most successful approaches to date for solving these impurity problems are the various continuous time quantum Monte Carlo algorithms. Here, we consider perhaps the simplest realization of multiorbital physics: the bilayer Hubbard model on an infinite-coordination Bethe lattice. Despite its simplicity, the majority of this model's phase diagram cannot be predicted by using traditional Monte Carlo methods. We show that these limitations can be largely circumvented by recently introduced inchworm Monte Carlo techniques. We then explore the model's phase diagram at a variety of interaction strengths, temperatures, and filling ratios.

arXiv:2311.16887•Physical review. B./Physical review. B•Feb 23, 2024•8 citations
Multiorbital quantum impurity models with general interaction and hybridization terms appear in a wide range of applications, including embedding, quantum transport, and nanoscience. However, most quantum impurity solvers are restricted to a few impurity orbitals, discretized baths, diagonal hybridizations, or density-density interactions. Here, we generalize the inchworm quantum Monte Carlo method to the interaction expansion, and we explore its application to typical single- and multiorbital problems encountered in investigations of impurity and lattice models. Our implementation generically outperforms bare and bold-line quantum Monte Carlo algorithms in the interaction expansion. For the systems studied here, our implementation remains inferior to the more specialized hybridization expansion and auxiliary field algorithms. The problem of convergence to unphysical fixed points, which hampers so-called bold-line methods, is not encountered in inchworm Monte Carlo.
We present a numerically exact inchworm Monte Carlo method for equilibrium multiorbital quantum impurity problems with general interactions and hybridizations. We show that the method, originally developed to overcome the dynamical sign problem in certain real-time propagation problems, can also overcome the sign problem as a function of temperature for equilibrium quantum impurity models. This is shown in several cases where the current method of choice, the continuous-time hybridization expansion, fails due to the sign problem. Our method therefore enables simulations of impurity problems as they appear in embedding theories without further approximations, such as the truncation of the hybridization or interaction structure or a discretization of the impurity bath with a set of discrete energy levels, and eliminates a crucial bottleneck in the simulation of ab initio embedding problems.
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 inchworm expansion is a promising approach to solving strongly correlated quantum impurity models due to its reduction of the sign problem in real and imaginary time. However, inchworm Monte Carlo is computationally expensive, converging as 1 / N where N is the number of samples. We show that the imaginary-time integration is amenable to quasi Monte Carlo, with parametrically better 1 / N convergence, by mapping the Sobol low-discrepancy sequence from the hypercube to the simplex with the so-called Root transform. This extends the applicability of the inchworm method to, e.g., multiorbital Anderson impurity models with off-diagonal hybridization, relevant for materials simulation, where continuous-time hybridization expansion Monte Carlo has a severe sign problem. Published by the American Physical Society 2024
The investigation of quantum impurity models plays a crucial role in condensed matter physics because of their wide-ranging applications, such as embedding theories and transport problems. Traditional methods often fall short in producing accurate results for multi-orbital systems with complex interactions and off-diagonal hybridizations. Recently, tensor-train-based integration and summation techniques have shown promise as effective alternatives. In this study, we use tensor train methods to tackle quantum impurity problems formulated within the imaginary-time inchworm hybridization expansion framework. We identify key challenges in the inchworm expansion itself and its interplay with tensor-train-based methods. We demonstrate the accuracy and versatility of our approach by solving general quantum impurity problems. Our results suggest that tensor-train decomposition schemes offer a viable path toward accurate and efficient multi-orbital impurity solvers.
We present a continuous-time Monte Carlo impurity solver for multiorbital impurity models which combines a strong-coupling hybridization expansion and a weak-coupling expansion in the Hund's coupling parameter JJJ. This double-expansion approach allows to treat the dominant density-density interactions UUU within the efficient segment representation. We test the approach for a two-orbital model with static interactions, and then explain how the double expansion allows to simulate models with frequency dependent $U(\ensuremath{\omega})$ and $J(\ensuremath{\omega})$. The method is used to investigate spin-state transitions in a toy model for fullerides, with repulsive bare JJJ but attractive screened JJJ.
Dynamical mean field theory allows access to the physics of strongly correlated materials with nontrivial orbital structure, but relies on the ability to solve auxiliary multiorbital impurity problems. The most successful approaches to date for solving these impurity problems are the various continuous time quantum Monte Carlo algorithms. Here, we consider perhaps the simplest realization of multiorbital physics: the bilayer Hubbard model on an infinite-coordination Bethe lattice. Despite its simplicity, the majority of this model's phase diagram cannot be predicted by using traditional Monte Carlo methods. We show that these limitations can be largely circumvented by recently introduced inchworm Monte Carlo techniques. We then explore the model's phase diagram at a variety of interaction strengths, temperatures, and filling ratios.