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

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

Shuta Matsuura, Hiroshi Shinaoka, Philipp Werner, Naoto Tsuji

Published Oct 14, 2026• 0 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1103/bvds-1gxf

Related papers

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

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

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

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

Continuous-time quantum Monte Carlo solver for dynamical mean field theory in the compact Legendre representation

Evan Sheridan, Cédric Weber, Evgeny A. Plekhanov +1 more

Dynamical mean-field theory (DMFT) is one of the most widely-used methods to treat accurately electron correlation effects in ab-initio real material calculations. Many modern large-scale implementations of DMFT in electronic structure codes involve solving a quantum impurity model with a continuous-time quantum Monte Carlo (CT-QMC) solver [Rubtsov et al., Phys. Rev. B 72, 035122 (2005); Werner et al., Phys. Rev. Lett. 97, 076405 (2006); Werner and Millis, Phys. Rev. B 74, 155107 (2006); Gull et al., Rev. Mod. Phys. 83, 349 (2011)]. The main advantage of CT-QMC is that, unlike standard quantum Monte Carlo approaches, it is able to generate the local Green's functions $G(\ensuremath{\tau})$ of the correlated system on an arbitrarily fine imaginary time $\ensuremath{\tau}$ grid, and is free of any systematic errors. In this work, we extend a hybrid QMC solver proposed by Khatami et al. [Phys. Rev. E 81, 056703 (2010)] and Rost et al. [Phys. Rev. E 87, 053305 (2013)] to a multiorbital context. This has the advantage of enabling impurity solver QMC calculations to scale linearly with inverse temperature $\ensuremath{\beta}$, and permit its application to $d$- and $f$-band materials. In addition, we present a Green's-function processing scheme which generates accurate quasicontinuous imaginary time solutions of the impurity problem which overcome errors inherent to standard QMC approaches. This solver and processing scheme are incorporated into a full DFT+DMFT calculation using the castep DFT code [Clark et al., Z. Kristallogr. 220, 567 (2005)]. Benchmark calculations for ${\mathrm{SrVO}}_{3}$ properties are presented. The computational efficiency of this method is also demonstrated.

arXiv:1811.12869•Physical review. B./Physical review. B•May 30, 2019•14 citations

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

Martin Eckstein

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.

arXiv:2410.19707•Oct 25, 2024•1 citation
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.
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.
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.
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.
Dynamical mean-field theory (DMFT) is one of the most widely-used methods to treat accurately electron correlation effects in ab-initio real material calculations. Many modern large-scale implementations of DMFT in electronic structure codes involve solving a quantum impurity model with a continuous-time quantum Monte Carlo (CT-QMC) solver [Rubtsov et al., Phys. Rev. B 72, 035122 (2005); Werner et al., Phys. Rev. Lett. 97, 076405 (2006); Werner and Millis, Phys. Rev. B 74, 155107 (2006); Gull et al., Rev. Mod. Phys. 83, 349 (2011)]. The main advantage of CT-QMC is that, unlike standard quantum Monte Carlo approaches, it is able to generate the local Green's functions $G(\ensuremath{\tau})$ of the correlated system on an arbitrarily fine imaginary time $\ensuremath{\tau}$ grid, and is free of any systematic errors. In this work, we extend a hybrid QMC solver proposed by Khatami et al. [Phys. Rev. E 81, 056703 (2010)] and Rost et al. [Phys. Rev. E 87, 053305 (2013)] to a multiorbital context. This has the advantage of enabling impurity solver QMC calculations to scale linearly with inverse temperature $\ensuremath{\beta}$, and permit its application to ddd- and fff-band materials. In addition, we present a Green's-function processing scheme which generates accurate quasicontinuous imaginary time solutions of the impurity problem which overcome errors inherent to standard QMC approaches. This solver and processing scheme are incorporated into a full DFT+DMFT calculation using the castep DFT code [Clark et al., Z. Kristallogr. 220, 567 (2005)]. Benchmark calculations for SrVO3{\mathrm{SrVO}}_{3}SrVO3​ properties are presented. The computational efficiency of this method is also demonstrated.
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.