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

Inchworm quasi Monte Carlo for quantum impurities

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

Published Sep 24, 2024• 6 citations
arXivHTMLPDF

Abstract

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

DOI
10.1103/physrevb.110.l121120

Related papers

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

Jia Li, Yang Yu, Emanuel Gull +1 more

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.

arXiv:2201.03099•Physical review. B./Physical review. B•Apr 18, 2022•12 citations

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

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

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

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

Taming the Dynamical Sign Problem in Real-Time Evolution of Quantum Many-Body Problems

Guy Cohen, Emanuel Gull, David R. Reichman +1 more

Current nonequilibrium Monte Carlo methods suffer from a dynamical sign problem that makes simulating real-time dynamics for long times exponentially hard. We propose a new "inchworm algorithm," based on iteratively reusing information obtained in previous steps to extend the propagation to longer times. The algorithm largely overcomes the dynamical sign problem, changing the scaling from exponential to quadratic. We use the method to solve the Anderson impurity model in the Kondo and mixed valence regimes, obtaining results both for quenches and for spin dynamics in the presence of an oscillatory magnetic field.

arXiv:1510.03534•Physical Review Letters•Dec 23, 2015•227 citations
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
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 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.
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.
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.
Current nonequilibrium Monte Carlo methods suffer from a dynamical sign problem that makes simulating real-time dynamics for long times exponentially hard. We propose a new "inchworm algorithm," based on iteratively reusing information obtained in previous steps to extend the propagation to longer times. The algorithm largely overcomes the dynamical sign problem, changing the scaling from exponential to quadratic. We use the method to solve the Anderson impurity model in the Kondo and mixed valence regimes, obtaining results both for quenches and for spin dynamics in the presence of an oscillatory magnetic field.