Zhen Huang, Denis Golež, Hugo U. R. Strand, Jason Kaye
SciPost Phys. 19, 121 (2025) ·
published 7 November 2025
|
· pdf
We present an efficient separation of variables algorithm for the evaluation of imaginary time Feynman diagrams appearing in the bold pseudo-particle strong coupling expansion of the Anderson impurity model. The algorithm uses a fitting method based on AAA rational approximation and numerical optimization to obtain a sum-of-exponentials expansion of the hybridization function, which is then used to decompose the diagrams. A diagrammatic formulation of the algorithm leads to an automated procedure for diagrams of arbitrary order and topology. We also present methods of stabilizing the self-consistent solution of the pseudo-particle Dyson equation. The result is a low-cost and high-order accurate impurity solver for quantum embedding methods using general multi-orbital hybridization functions at low temperatures, appropriate for low-to-intermediate expansion orders. In addition to other benchmark examples, we use our solver to perform a dynamical mean-field theory study of a minimal model of the strongly correlated compound Ca$_2$RuO$_4$, describing the anti-ferromagnetic transition and the in- and out-of-plane anisotropy induced by spin-orbit coupling.
SciPost Phys. 10, 091 (2021) ·
published 26 April 2021
|
· pdf
We propose a method to improve the computational and memory efficiency of numerical solvers for the nonequilibrium Dyson equation in the Keldysh formalism. It is based on the empirical observation that the nonequilibrium Green's functions and self energies arising in many problems of physical interest, discretized as matrices, have low rank off-diagonal blocks, and can therefore be compressed using a hierarchical low rank data structure. We describe an efficient algorithm to build this compressed representation on the fly during the course of time stepping, and use the representation to reduce the cost of computing history integrals, which is the main computational bottleneck. For systems with the hierarchical low rank property, our method reduces the computational complexity of solving the nonequilibrium Dyson equation from cubic to near quadratic, and the memory complexity from quadratic to near linear. We demonstrate the full solver for the Falicov-Kimball model exposed to a rapid ramp and Floquet driving of system parameters, and are able to increase feasible propagation times substantially. We present examples with 262144 time steps, which would require approximately five months of computing time and 2.2 TB of memory using the direct time stepping method, but can be completed in just over a day on a laptop with less than 4 GB of memory using our method. We also confirm the hierarchical low rank property for the driven Hubbard model in the weak coupling regime within the GW approximation, and in the strong coupling regime within dynamical mean-field theory.
