SciPost Submission Page
Multiloop functional renormalization group for the twodimensional Hubbard model: Loop convergence of the response functions
by Agnese Tagliavini, Cornelia Hille, Fabian B. Kugler, Sabine Andergassen, Alessandro Toschi, Carsten Honerkamp
This Submission thread is now published as
Submission summary
Authors (as registered SciPost users):  Sabine Andergassen · Cornelia Hille 
Submission information  

Preprint Link:  https://arxiv.org/abs/1807.02697v3 (pdf) 
Date accepted:  20181219 
Date submitted:  20181128 01:00 
Submitted by:  Hille, Cornelia 
Submitted to:  SciPost Physics 
Ontological classification  

Academic field:  Physics 
Specialties: 

Approaches:  Theoretical, Computational 
Abstract
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 selfenergy; (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 oneloop 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 "postprocessing" contraction of the interaction vertex at the end of the flow. The presented substantial refinement of fRGbased computation schemes paves a promising route towards future quantitative fRG analyses of more challenging systems and/or parameter regimes.
Author comments upon resubmission
\vskip 5mm
\noindent
We thank the Referee for carefully reading our manuscript and for posing constructive questions and suggestions to improve on the readibility of our work. We are glad that she/he finds {\sl ``[...] it quite impressive how the work improves in the fRG on so many fronts.''} and sees {``\sl [...] significant potential of this method for future applications to other (possibly more complicated) Hubbard models''}. Below we provide the clarification of all specific questions/remarks.
\vskip 5mm
\noindent
{\it 1. ``The Section 2.2 about the derivation of the multiloop fRG equations for the susceptibility is rather technical and I had problems understanding how this multiloop extension works. The section would profit much from a diagrammatic illustration of the scheme. In Fig. 1 an illustration is already given, but I didn't find it very enlightening."}\vspace{2mm}\\
\noindent
To this purpose, following the useful suggestion of the Referee, we added Fig.~1 in the revised version of the manuscript. The latter aims at giving a diagrammatic illustration how to construct the multiloop corrections upon the standard oneloop truncated flow equations of the susceptibility and the fermionboson vertex $\gamma_3$ in all physical channels. Similarly to the illustration of the multiloop flow equations for the channeldecomposed vertices in [Kugler Phys. Rev. B {\bf 97} (2018), Fig.~5], we proceed by showing the twoloop correctionterms (Fig.~1(b) in our revised version of the manuscript) and by generalizing the corrections to an arbitrary $2+\ell \geq 3$loop order (Fig.~1.(c)). We hope that Fig.~2 (Fig.~1 in the old version of the manuscript) can now be more clearly interpreted as an example of the lowest multiloop corrections, up to $\ell \!=\! 3$ and $\ell \!=\!2$ for the flow equation of the susceptibility and fermionboson vertex, respectively, at the leading order in the (bare) interaction.
\vskip 5mm
\noindent
{\it 2. ``As shown in Fig. 6 the real part of the self energy vanishes at two points at the Fermi surface; hence at these points the Fermi surface remains unchanged by the self energy feedback. Does this apply to the whole Fermi surface and for higher loop orders as well? In other words, does the perfect nesting effect always remain intact?''}\vspace{2mm}\\
The vanishing of the Fermi surface shift is caused by the particlehole symmetry of the model at half filling. In this situation, one can easily show that the selfenergy fulfill the symmetry property which is reported in the revised version of the manuscript as Eq.~53. In the text we added additional comments regarding the consequence of the particlehole symmetry on the structure of the selfenergy, and hence on the spectral properties of the system at the highly symmetric points of the Brillouine Zone, which have been shown in Fig.~7 (Fig.~6 in the old version of the manuscript). As specified in the text {\it `` This symmetry is not violated by any of the perturbative corrections and also not by the numerical implementation (e.g. the choice of kpoints in the BZ).."}. The confirmation that the particlehole symmetry of the selfenergy is preserved at higher looporder has been numerically displayed in Fig.~7 (Fig.~6 in the old version of the manuscript) where we added the $\ell=8$ mfRG results for the real part (left panel) and imaginary part (right panel) of the selfenergy.\\ \noindent Interestingly, the higher looporder results do not show any appreciable difference in the bandwidth renormalization with respect to the corresponding $\ell=1$ ones, and the multiloop corrections provide only a slight reduction of the scattering rate (more enhanced at the FS) for all ${\bf k}$points considered.
\vskip 5mm
\noindent
{\it 3. ``The authors argue that in the limit of large loop numbers, the MerminWagner theorem should be fulfilled. On the other hand, their numerical results in Fig. 10 still show a significant violation of MerminWagner even at l=8. The authors provide arguments why it is very challenging to suppress the pseudocritical temperature in the absence of extremely longrange fluctuations. However, given the fact that the fRG is often criticized because of the violation of MerminWagner, it would be desirable to see a more detailed discussion here. What I mean is, if a method that is known to fulfill MerminWagner still finds a finite ordering temperature, then this temperature should have a physical meaning. For example, it should be related to the extent of correlations taken into account and it should maybe even be possible to calculate/estimate it without solving the fRG equations. Can the authors comment on this?''}\vspace{2mm}\\
\noindent
We thank the Referee for the comment and we fully agree that this point deserves a more precise clarification (which we have also included, more shortly, in the revised version of our manuscript).
\noindent
The fulfillment of the MerminWagner (MW) Theorem in mfRG is formally guaranteed by the fact that for $\ell \rightarrow \infty$, one reproduces the PA per construction (which was demonstrated, in turn, to guarantee the MW theorem). Nonetheless, the extraction of the pseudocritical temperature from a set of a numerical mfRG data (e.g, our $\ell=8$ data of Fig.~11) is highly non trivial, because of other issues.
\noindent
One is physical, as the fulfillment of the MW theorem implies the emergence, at a very lowT scale (let's call it here : $T^{*}$), of antiferromagnetic (AF) fluctuations, growing exponentially (and no longer with a power law) for $T < T^{*}$. This is also reflected into the associated AF (inverse) susceptibility which will display an highT linear behavior for $T \gg T^{*}$ (where the spatial fluctuations are short range and an effective meanfield might work reasonably), which is eventually reverted in an exponential behavior for $T \ll T^{*}$.
\noindent The second reason is computational, as the convergence of the mfRG calculations w.r.t. to the momentum and frequencygrid becomes particularly though in the aforementioned lowT domain of quasilong range fluctuations. The ability to correctly capture the latter is essential to obtain the suppression of $T^{*}$ down to exactly $T\!=\!0$.
\noindent
On the basis of the above considerations, the interpretation of our multiloop mfRG data becomes rather clear. In the region where our higher loop calculations could be converged w.r.t. the underlying momentum/frequencygrids, we have also obtained loopconverged results (namely, $\ell=8$) to the PA. This indicates a systematic suppression of $\chi_{\text{AF}}$ w.r.t. the corresponding $\ell$loop results, due to the additional fluctuations included, which is a trend in perfect accordance with the theoretical expectations based on the MW theorem.
At the same time, we note that the temperatures considered are certainly larger than $T^{*}$, as our (loopconverged, PA) data for $\chi_{\text{AF}}^{1}(T)$ display a linear behavior.
Within the present implementation, we could not get, instead, momentumconverged $8$loop results in the (quasilong range fluctuations regime) of $T < T^{*}$ of the $U=2t$ halffilled case. Hence, it is simply not possible to extrapolate a corresponding $T_{\text{pc}}$ from our mfRG converged data, because they are all lying in the ``effective meanfield" regime. Such an extrapolation could be rather regarded as the instability of a renormalized meanfield theory (such as, e.g., the DMFT), which yields typically much lower instability temperatures than the standard (static) meanfield. This interpretation is fully supported by quantitative estimates in our case. In fact, by linearly extrapolating our $8$loop results, we would get a temperature $\sim 0.07 t$, which is much lower than the corresponding RPA estimates of $\sim 0.21 t$ [C. Eckhardt et al., PRB 2018)] but quantitatively in agreement with the DMFT results [see Sec. IV in RMP 90 025003 (2018)]. We note that, physically, this Tscale would be roughly consistent with the T interval where a sharp crossover from a paramagnetic metallic to a paramagnetic insulating state is found, and, hence, at weakcoupling, also with the aforementioned $T^{*}$scale.
The latter interpretation finds further support in the numerical results of diagrammatic extensions of DMFT applied to the lowT regime of the 2d Hubbard model [see: T. Sch\"afer et al, PRB 2015, and Sec. IV in RMP 90 025003 (2018)].
\vskip 5mm
\noindent
{\it 4. ``Besides quantitative changes in the results when applying their extended fRG scheme to the twodimensional Hubbard model, do the authors also find qualitatively new properties of this system? If yes, I think it should be better highlighted. ''}\vspace{2mm}\\
The 2D Hubbard model at halffilling represents an highly investigated system whose underlying features have been studied by means of a huge variety of methods (see, e.g., the extensive review by J. P. F. LeBlanc et. al. [Phys. Rev. X (2015)]).
At the present stage, the description of physical properties computed by means of the different approaches appears to match qualitatively, and in significant parameter region also quantitatively, if calculations are performed in the particlehole symmetric sector of this model, where the physics is controlled essentially by antiferromagnetic correlations.
The 2D Hubbard model at halffilling represents thus a rather solid and, at the same time, physically nottrivial situation to verify the applicability of our multiloopextended fRG scheme.
\noindent The fRG method itself is, however, not restricted to the particlehole symmetric sector, or to systems dominated by a specific type of correlations. On the contrary, one of the strength of the fRG, proven by the several 1loop studies in literature [cf. W. Metzner et al, RMP (2012)], is the possibility to analyze situations where the landscape of instabilities and emergent energy scales is not clearly settled a priori. To be more specific, we recall that, over the last decades, fRG has been extensively applied in more general situations, with a considerable number of studies of the 2D Hubbard model outofhalf filling. These studies were performed, however, with a lower quantitative control w.r.t. our presented algorithm (see Sec.~3), which motivates the effort of the extensions presented in our work.
\noindent Hence, summarizing the considerations above, because of the the strong leading antiferromagnetic fluctuations dominating the halffilling/particlehole symmetric case, the highloop corrections can not (and do not) lead to any qualitative change of the results found for $\ell=1$, {\sl but} for the fulfillment of the MerminWagner theorem, i.e., for the disappearance of any pseudocritical (magnetic) divergence wherever the mfRG calculation can be fully converged. The discussion of this point has been refined in the revised manuscript (see reply to the previous point). The highly studied sector of halffilling represents, in any case, the most natural testbed to benchmark and present the fundamental performance of our improved approach.
\noindent On a longer time perspective, nonparticlehole symmetric situations offer a very promising ``playground" for our newly developed mfRG algorithm. The situation in which fluctuations of different nature become simultaneously large represents an evident case, where multiloop corrections could play an crucial role in modifying, qualitatively, the phase diagram of the system. In fact, as discussed in the manuscript, one of the most important effect of the multiloop corrections is to change the balance between the strength of the different fluctuating channels. Hence, in the proximity of competing instabilities, such changes might lead to radically different outcomes w.r.t. some previous fRG studies [C. J. Halboth and W. Metzner, Phys. Rev. B {\bf 61} (2000), D. Zanchi and H. J. Schulz Phys. Rev. B {\bf 61},C. J. Halboth and W. Metzner, Phys. Rev. Lett. {\bf 24} (2000), C. Husemann and M. Salmhofer Phys. Rev. B {\bf 79} (2009), K.W. Giering and M. Salmhofer Phys. Rev. B {\bf 86} (2012)].
\noindent In the resubmitted version, we have amended the conclusion section with a similar statement to make this point clearer.
\vskip 5mm
\noindent
{\it 5. `It would be very interesting to know about the numerical efforts of these calculations. How long did they take and on how many CPUs did they run? Do these calculations require massive parallelization and supercomputing facilities?''}\vspace{2mm}\\
\noindent
To this purpose we have added an additional subsection, namely Sec.3.3 {\it "Performance of the code"}, which should provide the reader information regarding the multithreading parallelization of the code and a rough estimate of the computing time needed for the longest calulation perfomed in the paper. The scaling with respect to a different parameter choice can be deduced by looking at Fig.~3 (Fig.~2 in the old version of the manuscript).
\vskip 5mm
\noindent
{\it 6. `The numbering of the subsections in Sec. 4 is very strange as it starts with a zeroth subsection.''}\vspace{2mm}\\
\noindent
We thank the Referee for noticing this error which has occured in the formatting process. In the revised version of our manuscript, the enumeration of Sec.4 has been corrected.
List of changes
{\bf LIST OF CHANGES}
\vskip 5mm
\begin{enumerate}[label=\roman{*}), ref=(\roman{*})]
\item We have added the new Fig.~1 which illustrates the diagrammatic representation of the multiloop flow equation for the susceptibility and the fermionboson vertex.
\item We have added to Fig.~7 (left and right) mfRG results for the selfenergy corresponding to $\ell=8$ .
\item We have added a short discussion regarding the properties of the selfenergy at paticlehole symmetry referring to the $1\ell$ results (see Sec.~4.4) and to the mfRG results (see Sec. 4.5).
\item We have added a short sentence in Sec.~4.5 p.~2627 \textcolor{red}{@Ale:todo?}.
\item We have corrected the numbering of the subsections in Sec.~4 according to the remark (6) from the Referee.
\item We have added some sentences in the Conclusion, to answer question (5) from the Referee.
\item Funding information? Cluster at MPI.
\item We added Figure 14 (Appendix A) for a better understanding of the frequency and
momentum notations used.
\item We have corrected a typo in Eq.~53a)c) in Appendix A and in Eq.~92e) in Appendix F.
\item We have added Appendix G {\it ``Performance of the code"}.
\item We have corrected typos throughout the manuscript.
\end{enumerate}
Published as SciPost Phys. 6, 009 (2019)
Reports on this Submission
Report #1 by Anonymous (Referee 2) on 20181210 (Invited Report)
 Cite as: Anonymous, Report on arXiv:1807.02697v3, delivered 20181210, doi: 10.21468/SciPost.Report.730
Strengths
(1) Manuscript develops and discusses a very advanced manybody method
(2) This method is applied to a complicated model (the 2D Hubbard model)
(3) Many convergence and benchmark checks are performed
(4) The developed method has much potential for future applications
Weaknesses
(1) The results do not fulfill the MerminWagner theorem but the reasons for this failure are well discussed
Report
The resubmitted manuscript adequately addresses all suggestions and questions raised in my previous report. It features a more expanded discussion of various important aspects of this work. I remain to be impressed by the fact that this paper improves the fRG in so many different directions. I am, therefore, happy to recommend the manuscript for publication.
Sabine Andergassen on 20191025 [id 633]
A few minor corrections on this paper, including some typos in formulas, are included in un updated arXiv:1807.02697 version.