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

TL;DR

This paper introduces a tensor train compression method for four-point vertices obtained from multipoint NRG, enabling high-resolution frequency evaluations crucial for advanced electronic structure calculations.

Contribution

It presents a novel tensor train compression technique using quantics tensor cross interpolation for four-point vertices from mpNRG, significantly improving frequency grid resolution capabilities.

Findings

Achieved high-accuracy vertex representations with low bond dimensions.

Enabled evaluations on frequency grids far beyond previous methods.

Demonstrated effectiveness on the single-impurity Anderson model.

Abstract

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 non-perturbatively across a wide range of interactions and temperatures. It gives access to a local, non-perturbative 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 its asymptotic decomposition. For imaginary frequencies, the full vertex can be represented to an accuracy on the order of $2\cdot 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 $\lesssim 2\%$. Our work marks another step toward tensor-train-based diagrammatic calculations for correlated electronic lattice models starting from a local, non-perturbative mpNRG vertex.