Strong coupling impurity solver based on quantics tensor cross interpolation

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.