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.22975 · Physical Review Research

Two-particle calculations with quantics tensor trains: Solving the parquet equations

Stefan Rohshap, Marc K. Ritter, Hiroshi Shinaoka, Jan von Delft, Markus Wallerberger, Anna Kauch

Published Apr 24, 2025• 17 citations
arXivHTMLPDF

Abstract

We present an application of quantics tensor trains (QTTs) and tensor cross interpolation (TCI) to the solution of a full set of self-consistent equations for multivariate functions, the so-called parquet equations. We show that the steps needed to evaluate the equations (Bethe-Salpeter equations, parquet equation, and Schwinger-Dyson equation) can be decomposed into basic operations on the QTT-TCI compressed objects. The repeated application of these operations does not lead to a loss of accuracy beyond a specified tolerance and the iterative scheme converges even for numerically demanding parameters. As examples, we take the Hubbard model in the atomic limit and the single impurity Anderson model, where the basic objects in parquet equations, the two-particle vertices, depend on three frequencies, but not on momenta. The results show that this approach is able to overcome major computational bottlenecks of standard numerical methods. The applied methods allow for an exponential increase of the number of grid points included in the calculations, and a corresponding exponential reduction of the computational error, for a linear increase in computational cost.

DOI
10.1103/physrevresearch.7.023087

Related papers

Memory-Efficient Nonequilibrium Green’s Function Framework Built On Quantics Tensor Trains

M. Środa, Ken Inayoshi, Hiroshi Shinaoka +1 more

One of the challenges in diagrammatic simulations of nonequilibrium phenomena in lattice models is the large memory demand for storing momentum-dependent two-time correlation functions. This problem can be overcome with the recently introduced quantics tensor train (QTT) representation of multivariable functions. Here, we demonstrate nonequilibrium Green's function simulations within the GW and Migdal approximations with high momentum resolution, up to times which exceed the capabilities of standard implementations and are long enough to study, e.g., transient Floquet physics during multicycle electric field pulses and thermalization dynamics. The self-consistent calculation on the three-leg Kadanoff-Baym contour is fully self-contained, employing only QTT-compressed functions and input functions which are either generated directly in QTT form or obtained via quantics tensor cross interpolation.

arXiv:2412.14032•Physical Review Letters•Oct 14, 2025•8 citations

Adaptive sampling-based optimization of quantics tensor trains for noisy functions: Applications to quantum simulations

Kohtaroh Sakaue, Hiroshi Shinaoka, Rihito Sakurai

Tensor cross interpolation (TCI) is a powerful technique for learning a tensor train (TT) by adaptively sampling a target tensor based on an interpolation formula. However, when the tensor evaluations contain random noise, optimizing the TT is more advantageous than interpolating the noise. Here, we propose a new method that starts with an initial guess of TT and optimizes it using non-linear least-squares by fitting it to measured points obtained from TCI. We use quantics TCI (QTCI) in this method and demonstrate its effectiveness on sine and two-time correlation functions, with each evaluated with random noise. The resulting QTT exhibits increased robustness against noise compared to the QTCI method. Furthermore, we employ this optimized QTT of the correlation function in quantum simulation based on pseudo-imaginary-time evolution, resulting in ground-state energy with higher accuracy than the QTCI or Monte Carlo methods.

arXiv:2405.12730•SciPost Physics•Aug 14, 2025•3 citations

Multiscale Space-Time Ansatz for Correlation Functions of Quantum Systems Based on Quantics Tensor Trains

Hiroshi Shinaoka, Markus Wallerberger, Yuta Murakami +4 more

The correlation functions of quantum systems-central objects in quantum field theories-are defined in high-dimensional space-time domains.Their numerical treatment thus suffers from the curse of dimensionality, which hinders the application of sophisticated many-body theories to interesting problems.Here, we propose a multiscale space-time ansatz for correlation functions of quantum systems based on quantics tensor trains (QTTs), "qubits" describing exponentially different length scales.The ansatz then assumes a separation of length scales by decomposing the resulting high-dimensional tensors into tensor trains (also known as matrix product states).We numerically verify the ansatz for various equilibrium and nonequilibrium systems and demonstrate compression ratios of several orders of magnitude for challenging cases.Essential building blocks of diagrammatic equations, such as convolutions or Fourier transforms, are formulated in the compressed form.We numerically demonstrate the stability and efficiency of the proposed methods for the Dyson and Bethe-Salpeter equations.The QTT representation provides a unified framework for implementing efficient computations of quantum field theories.

arXiv:2210.12984•Physical Review X•Apr 27, 2023•53 citations

Nonequilibrium diagrammatic many-body simulations with quantics tensor trains

Matthias Murray, Hiroshi Shinaoka, Philipp Werner

The nonequilibrium Green's function formalism provides a versatile and powerful framework for numerical studies of nonequilibrium phenomena in correlated many-body systems. For calculations starting from an equilibrium initial state, a standard approach consists of discretizing the Kadanoff-Baym contour and implementing a causal time-stepping scheme in which the self-energy of the system plays the role of a memory kernel. This approach becomes computationally expensive at long times, because of the convolution integrals and the large amount of computer memory needed to store the Green's functions. A recent idea for the compression of nonequilibrium Green's functions is the quantics tensor train representation. Here, we explore this approach by implementing equilibrium and nonequilibrium simulations of the two-dimensional Hubbard model with a second-order weak-coupling approximation to the self-energy. We show that calculations with compressed two-time functions are possible without any loss of accuracy, and that the quantics tensor train implementation shows a much improved scaling of the computational effort and memory demand with the length of the time contour.

arXiv:2312.03809•Physical review. B./Physical review. B•Apr 18, 2024•15 citations

Compressing local vertex functions from the multipoint numerical renormalization group using quantics tensor cross interpolation

Markus Frankenbach, Marc K. Ritter, Mathias Pelz +3 more

The multipoint numerical renormalization group (mpNRG) is a powerful impurity solver that provides accurate spectral data useful for computing local, dynamic correlation functions in imaginary or real frequencies nonperturbatively across a wide range of interactions and temperatures. It gives access to a local, nonperturbative four-point vertex in imaginary and real frequencies, which can be used as input for subsequent computations such as diagrammatic extensions of dynamical mean-field theory. However, computing and manipulating the real-frequency four-point vertex on large, dense grids quickly becomes numerically challenging when the density and/or the extent of the frequency grid is increased. In this paper, we compute four-point vertices in a strongly compressed quantics tensor train format using quantics tensor cross interpolation, starting from discrete partial spectral functions obtained from mpNRG. This enables evaluations of the vertex on frequency grids with resolutions far beyond the reach of previous implementations. We benchmark this approach on the four-point vertex of the single-impurity Anderson model across a wide range of physical parameters, both in its full form and in its asymptotic decomposition. For imaginary frequencies, the full vertex can be represented to an accuracy on the order of 2 × 10 − 3 with maximum bond dimensions not exceeding 120. The more complex full real-frequency vertex requires maximum bond dimensions not exceeding 170 for an accuracy of ≲ 2 % . Our work marks another step toward tensor-train-based diagrammatic calculations for correlated electronic lattice models starting from a local, nonperturbative mpNRG vertex.

arXiv:2506.13359•Physical Review Research•Aug 27, 2025•1 citation

Solving the Gross-Pitaevskii equation with quantic tensor trains: Ground states and nonlinear dynamics

Qian Can Chen, I-Kang Liu, Jheng Wei Li +1 more

© Q.-C. Chen et al. We develop a tensor network framework based on the quantic tensor train (QTT) format to efficiently solve the Gross-Pitaevskii equation (GPE), which governs Bose-Einstein condensates under mean-field theory. For sufficiently regular functions, the QTT approach enables exponentially fine grids with computational cost that grows only linearly. By adapting the time-dependent variational principle (TDVP) and gradient descent methods, we accurately handle the GPE’s nonlinearities within the QTT structure. Our approach enables high-resolution simulations with drastically reduced computational cost on highly refined grids. We benchmark the ground states and dynamics of BECs, including vortex lattice formation and breathing modes, demonstrating superior performance over regular-grid-based methods and stable long-time evolution due to saturating bond dimensions. This establishes QTT as a powerful tool for nonlinear quantum simulations.

•Aug 10, 2026•0 citations
We present an application of quantics tensor trains (QTTs) and tensor cross interpolation (TCI) to the solution of a full set of self-consistent equations for multivariate functions, the so-called parquet equations. We show that the steps needed to evaluate the equations (Bethe-Salpeter equations, parquet equation, and Schwinger-Dyson equation) can be decomposed into basic operations on the QTT-TCI compressed objects. The repeated application of these operations does not lead to a loss of accuracy beyond a specified tolerance and the iterative scheme converges even for numerically demanding parameters. As examples, we take the Hubbard model in the atomic limit and the single impurity Anderson model, where the basic objects in parquet equations, the two-particle vertices, depend on three frequencies, but not on momenta. The results show that this approach is able to overcome major computational bottlenecks of standard numerical methods. The applied methods allow for an exponential increase of the number of grid points included in the calculations, and a corresponding exponential reduction of the computational error, for a linear increase in computational cost.
One of the challenges in diagrammatic simulations of nonequilibrium phenomena in lattice models is the large memory demand for storing momentum-dependent two-time correlation functions. This problem can be overcome with the recently introduced quantics tensor train (QTT) representation of multivariable functions. Here, we demonstrate nonequilibrium Green's function simulations within the GW and Migdal approximations with high momentum resolution, up to times which exceed the capabilities of standard implementations and are long enough to study, e.g., transient Floquet physics during multicycle electric field pulses and thermalization dynamics. The self-consistent calculation on the three-leg Kadanoff-Baym contour is fully self-contained, employing only QTT-compressed functions and input functions which are either generated directly in QTT form or obtained via quantics tensor cross interpolation.
Tensor cross interpolation (TCI) is a powerful technique for learning a tensor train (TT) by adaptively sampling a target tensor based on an interpolation formula. However, when the tensor evaluations contain random noise, optimizing the TT is more advantageous than interpolating the noise. Here, we propose a new method that starts with an initial guess of TT and optimizes it using non-linear least-squares by fitting it to measured points obtained from TCI. We use quantics TCI (QTCI) in this method and demonstrate its effectiveness on sine and two-time correlation functions, with each evaluated with random noise. The resulting QTT exhibits increased robustness against noise compared to the QTCI method. Furthermore, we employ this optimized QTT of the correlation function in quantum simulation based on pseudo-imaginary-time evolution, resulting in ground-state energy with higher accuracy than the QTCI or Monte Carlo methods.
The correlation functions of quantum systems-central objects in quantum field theories-are defined in high-dimensional space-time domains.Their numerical treatment thus suffers from the curse of dimensionality, which hinders the application of sophisticated many-body theories to interesting problems.Here, we propose a multiscale space-time ansatz for correlation functions of quantum systems based on quantics tensor trains (QTTs), "qubits" describing exponentially different length scales.The ansatz then assumes a separation of length scales by decomposing the resulting high-dimensional tensors into tensor trains (also known as matrix product states).We numerically verify the ansatz for various equilibrium and nonequilibrium systems and demonstrate compression ratios of several orders of magnitude for challenging cases.Essential building blocks of diagrammatic equations, such as convolutions or Fourier transforms, are formulated in the compressed form.We numerically demonstrate the stability and efficiency of the proposed methods for the Dyson and Bethe-Salpeter equations.The QTT representation provides a unified framework for implementing efficient computations of quantum field theories.
The nonequilibrium Green's function formalism provides a versatile and powerful framework for numerical studies of nonequilibrium phenomena in correlated many-body systems. For calculations starting from an equilibrium initial state, a standard approach consists of discretizing the Kadanoff-Baym contour and implementing a causal time-stepping scheme in which the self-energy of the system plays the role of a memory kernel. This approach becomes computationally expensive at long times, because of the convolution integrals and the large amount of computer memory needed to store the Green's functions. A recent idea for the compression of nonequilibrium Green's functions is the quantics tensor train representation. Here, we explore this approach by implementing equilibrium and nonequilibrium simulations of the two-dimensional Hubbard model with a second-order weak-coupling approximation to the self-energy. We show that calculations with compressed two-time functions are possible without any loss of accuracy, and that the quantics tensor train implementation shows a much improved scaling of the computational effort and memory demand with the length of the time contour.
The multipoint numerical renormalization group (mpNRG) is a powerful impurity solver that provides accurate spectral data useful for computing local, dynamic correlation functions in imaginary or real frequencies nonperturbatively across a wide range of interactions and temperatures. It gives access to a local, nonperturbative four-point vertex in imaginary and real frequencies, which can be used as input for subsequent computations such as diagrammatic extensions of dynamical mean-field theory. However, computing and manipulating the real-frequency four-point vertex on large, dense grids quickly becomes numerically challenging when the density and/or the extent of the frequency grid is increased. In this paper, we compute four-point vertices in a strongly compressed quantics tensor train format using quantics tensor cross interpolation, starting from discrete partial spectral functions obtained from mpNRG. This enables evaluations of the vertex on frequency grids with resolutions far beyond the reach of previous implementations. We benchmark this approach on the four-point vertex of the single-impurity Anderson model across a wide range of physical parameters, both in its full form and in its asymptotic decomposition. For imaginary frequencies, the full vertex can be represented to an accuracy on the order of 2 × 10 − 3 with maximum bond dimensions not exceeding 120. The more complex full real-frequency vertex requires maximum bond dimensions not exceeding 170 for an accuracy of ≲ 2 % . Our work marks another step toward tensor-train-based diagrammatic calculations for correlated electronic lattice models starting from a local, nonperturbative mpNRG vertex.
© Q.-C. Chen et al. We develop a tensor network framework based on the quantic tensor train (QTT) format to efficiently solve the Gross-Pitaevskii equation (GPE), which governs Bose-Einstein condensates under mean-field theory. For sufficiently regular functions, the QTT approach enables exponentially fine grids with computational cost that grows only linearly. By adapting the time-dependent variational principle (TDVP) and gradient descent methods, we accurately handle the GPE’s nonlinearities within the QTT structure. Our approach enables high-resolution simulations with drastically reduced computational cost on highly refined grids. We benchmark the ground states and dynamics of BECs, including vortex lattice formation and breathing modes, demonstrating superior performance over regular-grid-based methods and stable long-time evolution due to saturating bond dimensions. This establishes QTT as a powerful tool for nonlinear quantum simulations.