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

Strong coupling impurity solver based on quantics tensor cross interpolation

Aaram J. Kim, Philipp Werner

Published Mar 10, 2025• 15 citations
arXivHTMLPDF

Abstract

Numerical methods capable of handling nonequilibrium impurity models are essential for the study of transport problems and the solution of the nonequilibrium dynamical mean-field theory (DMFT) equations. In the strong correlation regime, the self-consistently resummed hybridization expansion is an appealing strategy, which however has been employed so far mainly in the lowest-order noncrossing approximation. At higher orders, standard implementations become numerically costly, but a significant speed-up can be achieved by evaluating multidimensional integrals in an approximate factorized form. Here we develop a one-crossing approximation solver based on the recently introduced quantics tensor cross interpolation, and demonstrate its accuracy and efficiency with applications to the Anderson impurity model and nonequilibrium steady-state DMFT calculations for the Hubbard model.

DOI
10.1103/physrevb.111.125120

Related papers

Quantics tensor cross interpolation for high-order strong-coupling expansions

Kanghyeon Kim, Lei Geng, Philipp Werner +1 more

Real-time impurity solvers enable the study of transport phenomena and the description of nonequilibrium lattice systems within the framework of dynamical mean-field theory (DMFT). They also provide direct access to the spectral functions of both equilibrium and nonequilibrium systems. A widely used approach is the self-consistent strong-coupling expansion, whose lowest-order implementation corresponds to the non-crossing approximation. Higher-order implementations, however, are computationally demanding because the number of diagram topologies grows factorially with expansion order, while the evaluation of self-energies and Green's functions requires increasingly high-dimensional integrations. Here, we demonstrate that the latter challenge can be mitigated by employing quantics tensor cross interpolation in a variable-separated framework. Compared with the previously used scale-separated approach, the new scheme yields substantially lower bond dimensions and capacitates self-consistent steady-state DMFT calculations up to fourth order. We illustrate its performance with representative results for both equilibrium and photo-doped systems. In addition, we analyze the convergence of the strong-coupling expansion in the challenging noninteracting limit by computing diagrams up to sixth order. At this order, the onset of the asymptotic regime of the strong-coupling expansion becomes apparent, which allows the application of extrapolation techniques.

arXiv:2608.23160•Aug 24, 2026•0 citations

Weak-coupling tensor cross interpolation impurity solver for nonequilibrium dynamical mean-field theory

Shuta Matsuura, Hiroshi Shinaoka, Philipp Werner +1 more

Simulating nonequilibrium quantum many-body systems remains a major challenge due to the exponential growth of the computational complexity with real time. Here we implement a nonequilibrium impurity solver based on the weak-coupling expansion and the tensor cross interpolation (TCI), and apply it to nonequilibrium dynamical mean-field theory (DMFT). The method approximates the integrands of the high-dimensional integrals arising in the weak-coupling expansion in a tensor-train form, enabling efficient evaluations without stochastic sampling and thereby mitigating the sign problem affecting continuous-time quantum Monte Carlo (CT-QMC) methods. Benchmark calculations for an exactly solvable nonequilibrium impurity model agree well with the exact results and reveal a low-rank structure of the integrands. When applied to interaction-quench problems in the half-filled Hubbard model, the method reproduces fast thermalization at a critical interaction strength with accuracy comparable to CT-QMC. Away from half filling, where the sign problem becomes even more severe, the present approach remains well controlled, revealing a crossover instead of a sharply defined fast thermalization point in the 3/4-filled case. The solver can also be applied to steady-state DMFT problems, yielding accurate spectral functions in the metallic regime without analytic continuation.

arXiv:2607.00702•Physical review. B./Physical review. B•Oct 14, 2026•0 citations

Accelerating a Strong-Coupling Non-Equilibrium Steady-State Impurity Solver using (Quantics) Tensor Trains

Bastian Schindler, Martin Eckstein

Including higher order diagrammatic corrections to the strong-coupling expansion is mainly limited by the evaluation of high-dimensional, time-ordered integrals. In this work we present and compare four different parametrizations of the integrands in order to obtain a low-rank (quantics) tensor-train representation using tensor cross interpolation. Particular emphasis is placed on a quantics time-difference formulation in which the required retarded convolutions are performed directly in quantics tensor-train form. Using controlled Gaussian benchmarks, we analyze the accuracy, bond dimensions, and computational scaling of the different approaches. We then validate the most promising formulations in self-consistent equilibrium and nonequilibrium DMFT calculations and demonstrate calculations up to the third order in the strong-coupling expansion. Finally, we extend the solver to impurity models with retarded density-density interactions and apply it within nonequilibrium extended DMFT. Our results show that tensor cross interpolation substantially reduces the cost of evaluating higher-order diagrams and provides a controlled, systematically improvable framework for nonequilibrium quantum impurity calculations.

arXiv:2608.13146•Aug 13, 2026•0 citations

Third-order strong-coupling impurity solver for real-frequency dynamical mean field theory: Accurate spectral functions for antiferromagnetic and photodoped states

Lei Geng, Aaram J. Kim, Philipp Werner

Strong-coupling expansions for impurity models up to third order are implemented here on the real-time axis using a quantics tensor cross interpolation framework. This approach provides an accurate impurity solver that can describe both equilibrium systems and nonequilibrium steady states in the intermediate to strong-interaction regime.

arXiv:2507.20385•Physical review. B./Physical review. B•Nov 8, 2025•3 citations

Tensor cross interpolation approach for quantum impurity problems based on the weak-coupling expansion

Shuta Matsuura, Hiroshi Shinaoka, Philipp Werner +1 more

We apply the tensor cross interpolation (TCI) algorithm to solve equilibrium quantum impurity problems with high precision based on the weak-coupling expansion. The TCI algorithm, a kind of active learning method, factorizes high-dimensional integrals that appear in the perturbative expansion into a product of low-dimensional ones, enabling us to evaluate higher-order terms efficiently. This method is free from the sign problem which quantum Monte Carlo methods sometimes suffer from, and allows one to directly calculate the free energy. We benchmark the TCI impurity solver on an exactly solvable impurity model and find good agreement with the exact solutions. We also incorporate the TCI impurity solver into the dynamical mean-field theory to solve the Hubbard model and show that the metal-to-Mott insulator transition is correctly described with comparable accuracy to the Monte Carlo methods. Behind the effectiveness of the TCI approach for quantum impurity problems lies the fact that the integrands in the weak-coupling expansion naturally have a low-rank structure in the tensor-train representation.

arXiv:2501.12643•Physical review. B./Physical review. B•Apr 28, 2025•3 citations

Nonequilibrium dynamical mean-field calculations based on the noncrossing approximation and its generalizations

Martin Eckstein, Philipp Werner

We solve the impurity problem which arises within nonequilibrium dynamical mean-field theory for the Hubbard model by means of a self-consistent perturbation expansion around the atomic limit. While the lowest order, known as the noncrossing approximation (NCA), is reliable only when the interaction $U$ is much larger than the bandwidth, low-order corrections to the NCA turn out to be sufficient to reproduce numerically exact Monte Carlo results in a wide parameter range that covers the insulating phase and the metal-insulator crossover regime at not too low temperatures. As an application of the perturbative strong-coupling impurity solver we investigate the response of the double occupancy in the Mott insulating phase of the Hubbard model to a dynamical change in the interaction or the hopping, a technique which has been used as a probe of the Mott insulating state in ultracold fermionic gases.

arXiv:1005.1872•Physical Review B•Sep 17, 2010•219 citations
Numerical methods capable of handling nonequilibrium impurity models are essential for the study of transport problems and the solution of the nonequilibrium dynamical mean-field theory (DMFT) equations. In the strong correlation regime, the self-consistently resummed hybridization expansion is an appealing strategy, which however has been employed so far mainly in the lowest-order noncrossing approximation. At higher orders, standard implementations become numerically costly, but a significant speed-up can be achieved by evaluating multidimensional integrals in an approximate factorized form. Here we develop a one-crossing approximation solver based on the recently introduced quantics tensor cross interpolation, and demonstrate its accuracy and efficiency with applications to the Anderson impurity model and nonequilibrium steady-state DMFT calculations for the Hubbard model.
Real-time impurity solvers enable the study of transport phenomena and the description of nonequilibrium lattice systems within the framework of dynamical mean-field theory (DMFT). They also provide direct access to the spectral functions of both equilibrium and nonequilibrium systems. A widely used approach is the self-consistent strong-coupling expansion, whose lowest-order implementation corresponds to the non-crossing approximation. Higher-order implementations, however, are computationally demanding because the number of diagram topologies grows factorially with expansion order, while the evaluation of self-energies and Green's functions requires increasingly high-dimensional integrations. Here, we demonstrate that the latter challenge can be mitigated by employing quantics tensor cross interpolation in a variable-separated framework. Compared with the previously used scale-separated approach, the new scheme yields substantially lower bond dimensions and capacitates self-consistent steady-state DMFT calculations up to fourth order. We illustrate its performance with representative results for both equilibrium and photo-doped systems. In addition, we analyze the convergence of the strong-coupling expansion in the challenging noninteracting limit by computing diagrams up to sixth order. At this order, the onset of the asymptotic regime of the strong-coupling expansion becomes apparent, which allows the application of extrapolation techniques.
Simulating nonequilibrium quantum many-body systems remains a major challenge due to the exponential growth of the computational complexity with real time. Here we implement a nonequilibrium impurity solver based on the weak-coupling expansion and the tensor cross interpolation (TCI), and apply it to nonequilibrium dynamical mean-field theory (DMFT). The method approximates the integrands of the high-dimensional integrals arising in the weak-coupling expansion in a tensor-train form, enabling efficient evaluations without stochastic sampling and thereby mitigating the sign problem affecting continuous-time quantum Monte Carlo (CT-QMC) methods. Benchmark calculations for an exactly solvable nonequilibrium impurity model agree well with the exact results and reveal a low-rank structure of the integrands. When applied to interaction-quench problems in the half-filled Hubbard model, the method reproduces fast thermalization at a critical interaction strength with accuracy comparable to CT-QMC. Away from half filling, where the sign problem becomes even more severe, the present approach remains well controlled, revealing a crossover instead of a sharply defined fast thermalization point in the 3/4-filled case. The solver can also be applied to steady-state DMFT problems, yielding accurate spectral functions in the metallic regime without analytic continuation.
Including higher order diagrammatic corrections to the strong-coupling expansion is mainly limited by the evaluation of high-dimensional, time-ordered integrals. In this work we present and compare four different parametrizations of the integrands in order to obtain a low-rank (quantics) tensor-train representation using tensor cross interpolation. Particular emphasis is placed on a quantics time-difference formulation in which the required retarded convolutions are performed directly in quantics tensor-train form. Using controlled Gaussian benchmarks, we analyze the accuracy, bond dimensions, and computational scaling of the different approaches. We then validate the most promising formulations in self-consistent equilibrium and nonequilibrium DMFT calculations and demonstrate calculations up to the third order in the strong-coupling expansion. Finally, we extend the solver to impurity models with retarded density-density interactions and apply it within nonequilibrium extended DMFT. Our results show that tensor cross interpolation substantially reduces the cost of evaluating higher-order diagrams and provides a controlled, systematically improvable framework for nonequilibrium quantum impurity calculations.
We apply the tensor cross interpolation (TCI) algorithm to solve equilibrium quantum impurity problems with high precision based on the weak-coupling expansion. The TCI algorithm, a kind of active learning method, factorizes high-dimensional integrals that appear in the perturbative expansion into a product of low-dimensional ones, enabling us to evaluate higher-order terms efficiently. This method is free from the sign problem which quantum Monte Carlo methods sometimes suffer from, and allows one to directly calculate the free energy. We benchmark the TCI impurity solver on an exactly solvable impurity model and find good agreement with the exact solutions. We also incorporate the TCI impurity solver into the dynamical mean-field theory to solve the Hubbard model and show that the metal-to-Mott insulator transition is correctly described with comparable accuracy to the Monte Carlo methods. Behind the effectiveness of the TCI approach for quantum impurity problems lies the fact that the integrands in the weak-coupling expansion naturally have a low-rank structure in the tensor-train representation.
We solve the impurity problem which arises within nonequilibrium dynamical mean-field theory for the Hubbard model by means of a self-consistent perturbation expansion around the atomic limit. While the lowest order, known as the noncrossing approximation (NCA), is reliable only when the interaction UUU is much larger than the bandwidth, low-order corrections to the NCA turn out to be sufficient to reproduce numerically exact Monte Carlo results in a wide parameter range that covers the insulating phase and the metal-insulator crossover regime at not too low temperatures. As an application of the perturbative strong-coupling impurity solver we investigate the response of the double occupancy in the Mott insulating phase of the Hubbard model to a dynamical change in the interaction or the hopping, a technique which has been used as a probe of the Mott insulating state in ultracold fermionic gases.