Two-particle calculations with quantics tensor trains - solving the parquet equations

We present the first application of quantics tensor trains (QTTs) [1] and tensor cross interpolation (TCI) [2] 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 (QTCI) [3] compressed objects. The repeated application of these operations does not lead to the 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 leading to an exponentially improving computational error for a linear increase in computational cost.

[1] H. Shinaoka, et al., Physical Review X 13, 021015 (2023)

[2] Y. Nunez Fernandez, et al., Phys. Rev. X 12, 041018 (2022)

[3] M. Ritter, et al., Phys. Rev. Lett. 132, 056501 (2024)