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:2410.19707

Solving quantum impurity models in the non-equilibrium steady state with tensor trains

Martin Eckstein

Published Oct 25, 2024• 1 citation
arXivHTMLPDF

Abstract

We discuss the evaluation of the integrals for intermediate-order diagrams in the self-consistent strong-coupling expansion on the Keldysh contour using Tensor Cross Interpolation (TCI). TCI is used to factorize the nested parts of the integrand, allowing the integral to be computed as a recursion of convolution integrals, which are efficiently evaluated using the Fast Fourier Transform. The evaluation of diagrams where all vertices lie on one branch of the Keldysh contour resembles the structure of the imaginary-time formalism. For diagrams with time arguments on both contour branches, we find that it can be advantageous to parametrize the integrals in terms of physical time arguments with an additional sum over Keldysh indices. We benchmark the solution in relevant test cases, including the single impurity Anderson model and an exactly solvable electron-boson model. While the bond dimension increases with diagram order, the TCI-based integration efficiently handles low-order diagrams, making it a promising approach to go beyond the non-crossing approximation in steady-state non-equilibrium dynamical mean-field theory simulations.

DOI
10.48550/arxiv.2410.19707

Related papers

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

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

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

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

Strong coupling impurity solver based on quantics tensor cross interpolation

Aaram J. Kim, Philipp Werner

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.

arXiv:2411.19026•Physical review. B./Physical review. B•Mar 10, 2025•15 citations

Learning Feynman Diagrams with Tensor Trains

Y. Núñez-Fernández, Matthieu Jeannin, Philipp T. Dumitrescu +4 more

We use tensor network techniques to obtain high-order perturbative diagrammatic expansions for the quantum many-body problem at very high precision.The approach is based on a tensor train parsimonious representation of the sum of all Feynman diagrams, obtained in a controlled and accurate way with the tensor cross interpolation algorithm.It yields the full time evolution of physical quantities in the presence of any arbitrary time-dependent interaction.Our benchmarks on the Anderson quantum impurity problem, within the real-time nonequilibrium Schwinger-Keldysh formalism, demonstrate that this technique supersedes diagrammatic quantum Monte Carlo by orders of magnitude in precision and speed, with convergence rates 1=N 2 or faster, where N is the number of function evaluations.The method also works in parameter regimes characterized by strongly oscillatory integrals in high dimension, which suffer from a catastrophic sign problem in quantum Monte Carlo calculations.Finally, we also present two exploratory studies showing that the technique generalizes to more complex situations: a double quantum dot and a single impurity embedded in a two-dimensional lattice.

arXiv:2207.06135•Physical Review X•Nov 16, 2022•87 citations
We discuss the evaluation of the integrals for intermediate-order diagrams in the self-consistent strong-coupling expansion on the Keldysh contour using Tensor Cross Interpolation (TCI). TCI is used to factorize the nested parts of the integrand, allowing the integral to be computed as a recursion of convolution integrals, which are efficiently evaluated using the Fast Fourier Transform. The evaluation of diagrams where all vertices lie on one branch of the Keldysh contour resembles the structure of the imaginary-time formalism. For diagrams with time arguments on both contour branches, we find that it can be advantageous to parametrize the integrals in terms of physical time arguments with an additional sum over Keldysh indices. We benchmark the solution in relevant test cases, including the single impurity Anderson model and an exactly solvable electron-boson model. While the bond dimension increases with diagram order, the TCI-based integration efficiently handles low-order diagrams, making it a promising approach to go beyond the non-crossing approximation in steady-state non-equilibrium dynamical mean-field theory simulations.
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.
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.
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.
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.
We use tensor network techniques to obtain high-order perturbative diagrammatic expansions for the quantum many-body problem at very high precision.The approach is based on a tensor train parsimonious representation of the sum of all Feynman diagrams, obtained in a controlled and accurate way with the tensor cross interpolation algorithm.It yields the full time evolution of physical quantities in the presence of any arbitrary time-dependent interaction.Our benchmarks on the Anderson quantum impurity problem, within the real-time nonequilibrium Schwinger-Keldysh formalism, demonstrate that this technique supersedes diagrammatic quantum Monte Carlo by orders of magnitude in precision and speed, with convergence rates 1=N 2 or faster, where N is the number of function evaluations.The method also works in parameter regimes characterized by strongly oscillatory integrals in high dimension, which suffer from a catastrophic sign problem in quantum Monte Carlo calculations.Finally, we also present two exploratory studies showing that the technique generalizes to more complex situations: a double quantum dot and a single impurity embedded in a two-dimensional lattice.