Clemens Watzenböck, Martina Fellinger, Karsten Held, Alessandro Toschi
SciPost Phys. 12, 184 (2022) ·
published 7 June 2022
|
· pdf
We investigate the onset of a not-decaying asymptotic behavior of temporal magnetic correlations in the Hubbard model in infinite dimensions. This long-term memory feature of dynamical spin correlations can be precisely quantified by computing the difference between the zero-frequency limit of the Kubo susceptibility and the corresponding static isothermal one. Here, we present a procedure for reliably evaluating this difference starting from imaginary time-axis data, and apply it to the testbed case of the Mott-Hubbard metal-insulator transition (MIT). At low temperatures, we find long-term memory effects in the entire Mott regime, abruptly ending at the first order MIT. This directly reflects the underlying local moment physics and the associated degeneracy in the many-electron spectrum. At higher temperatures, a more gradual onset of an infinitely-long time-decay of magnetic correlations occurs in the crossover regime, not too far from the Widom line emerging from the critical point of the MIT. Our work has relevant algorithmic implications for the analytical continuation of dynamical susceptibilities in strongly correlated regimes and offers a new perspective for unveiling fundamental properties of the many-particle spectrum of the problem under scrutiny.
SciPost Phys. 6, 009 (2019) ·
published 18 January 2019
|
· pdf
We present a functional renormalization group (fRG) study of the two dimensional Hubbard model, performed with an algorithmic implementation which lifts some of the common approximations made in fRG calculations. In particular, in our fRG flow; (i) we take explicitly into account the momentum and the frequency dependence of the vertex functions; (ii) we include the feedback effect of the self-energy; (iii) we implement the recently introduced multiloop extension which allows us to sum up {\emph{all}} the diagrams of the parquet approximation with their exact weight. Due to its iterative structure based on successive one-loop computations, the loop convergence of the fRG results can be obtained with an affordable numerical effort. In particular, focusing on the analysis of the physical response functions, we show that the results become {\emph{independent}} from the chosen cutoff scheme and from the way the fRG susceptibilities are computed, i.e., either through flowing couplings to external fields, or through a "post-processing" contraction of the interaction vertex at the end of the flow. The presented substantial refinement of fRG-based computation schemes paves a promising route towards future quantitative fRG analyses of more challenging systems and/or parameter regimes.