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.