SciPost Submission Page
Nonperturbative intertwining between spin and charge correlations: A "smoking gun" singlebosonexchange result
by Severino Adler, Friedrich Krien, Patrick ChalupaGantner, Giorgio Sangiovanni, Alessandro Toschi
This is not the latest submitted version.
This Submission thread is now published as
Submission summary
Authors (as registered SciPost users):  Severino Adler · Giorgio Sangiovanni 
Submission information  

Preprint Link:  https://arxiv.org/abs/2212.09693v4 (pdf) 
Date submitted:  20231117 22:03 
Submitted by:  Adler, Severino 
Submitted to:  SciPost Physics 
Ontological classification  

Academic field:  Physics 
Specialties: 

Approaches:  Theoretical, Computational 
Abstract
We study the microscopic mechanism controlling the interplay between local charge and local spin fluctuations in correlated electron systems via a thorough investigation of the generalized onsite charge susceptibility of several fundamental manyelectron models, such as the Hubbard atom, the Anderson impurity model, and the Hubbard model. By decomposing the numerically determined generalized susceptibility in terms of physically transparent singleboson exchange processes, we unveil the microscopic mechanisms responsible for the breakdown of the selfconsistent manyelectron perturbation expansion. In particular, we unambiguously identify the origin of the significant suppression of its diagonal entries in (Matsubara) frequency space and the slight increase of the offdiagonal ones which cause the breakdown. The suppression effect on the diagonal elements originates directly from the electronic scattering on local magnetic moments, reflecting their increasingly longer lifetime as well as their enhanced effective coupling with the electrons. Instead, the slight and diffuse enhancement of the offdiagonal terms can be mostly ascribed to multiboson scattering processes. The strong intertwining between spin and charge sectors is partly weakened at the Kondo temperature due to a progressive reduction of the effective spinfermion coupling of local magnetic fluctuations in the low frequency regime. Our analysis, thus, clarifies the precise mechanism through which the physical information is transferred between different scattering channels of interacting electron problems and highlights the pivotal role played by such an intertwining in the physics of correlated electrons beyond the perturbative regime.
List of changes
1. ) Slightly refined the wording of the abstract, for better readability.
2.) Added a footnote in section 2.3 referring to different ways of defining the crossover border of the local moment regime. A more extended explanation has been included at the end of Appendix D.
3. ) Added a second footnote in section 2.3 referring to the first order jump of the "fingerprint" at the Mott transition.
4. ) Added reference to PRB 105, 155151 (2022) in various points.
5.) Fixed the typo on page 8 (now on page 9) Fig. ??. It now reads Fig. 14
6.) Dropped the sentence "a relatively flat (Curie) behavior of the quantity $T\chi^{\text{sp}}_{\omega=0}(T)$," at the beginning of Sec.~3.2.
7.) Dropped ", yielding a much less suppressed local moment formation" at the end of Sec.~3.2.
8.) Slightly refined the wording in Sec.~4
9.) Changed $\chi^{\text{sp}}_{0}$ to $\chi^{\text{sp}}_{\omega=0}$ in Appendix B and slightly modified Eq.~(25) and the sentence above.
10.) Changed the scaling of the plots in Fig.~14 to fit the normalization of the rest of the manuscript. Clarified the meaning of the blue shaded area in the corresponding figure caption.
11.) Update/slight extension of the manuscript's bibliography.
12.) Corrected small typos in the manuscript.
Current status:
Reports on this Submission
Anonymous Report 2 on 202416 (Invited Report)
 Cite as: Anonymous, Report on arXiv:2212.09693v4, delivered 20240106, doi: 10.21468/SciPost.Report.8377
Report
In this theoretical work the authors investigate the Matsubara structure of the local charge susceptibility for some paradigmatic correlated systems.
The goal of this analysis is to identify the microscopic mechanisms controlling, on the twoparticle level, the physics of charge localization and its relation with the formation of local moments. To address this issue they exploit the numerical decomposition of the generalized susceptibility in terms of processes involving the exchange of single welldefined bosonic modes.
The main result of the study is to unveil the microscopic intertwined between different fluctuating channel. In particular, the study demonstrates that, in the local moment regime, the suppression of local charge fluctuations directly originates from the exchange of static spin fluctuations. The type of analysis performed nicely show that the effect is strongly concentrated for ν = ν ′ as a precise consequence of the long lifetime of the local moment in this parameter regime.
The topic of the research is of high interest for the community working in the wide field of manybody methods and correlated systems. Although highly technical, the paper addresses very intuitive and physical questions. The work performed here, not only provides valuable insight in the physics of correlated material, but also a new strategy to analyze correlation effects at a microscopic level. Such method could be definitely exploited in the future to study more realistic many body models (e.g. including multiorbital physics, non local interactions and more …).
The quality of the research is excellent and the presentation of the results is very good.
The organization of the work, the choice of models and parameters and the technical aspects of the computation are well discussed.
I also revised previous referees’ reports and corresponding answers, and I find the reply and relative revision to the paper satisfactory.
After consideration of all these aspects, I strongly recommend the paper for publication. In what follows I mention some minor suggestions that I invite the authors to consider:
1. Some physical aspects, although quite basic could be not obvious for not experts. For example what is the physical meaning of diagonal and notdiagonal component in the Matsubara space in terms of static and dynamic fluctuations, the physical meaning of the matric eigenvalues, as well as the need to invert the generalized susceptibility. The authors could spend a few words in the introductory sections to make the reader able to better appreciate the physics behind the quantities considered.
2. To help the reader I would use the extensive names of the models in Section 2.2 and define here the acronym used in the following part of the paper.
3. The paper is already very long and dense. I suggest the authors to have a last revision and try to rephrase/shorten a bit the text whenever is possible.
4. While the results of the approximations can be briefly mentioned within the discussion to further support the findings, I would suggest to remove “3.4 Insights from Approximations” from the main text and create an appendix.
Requested changes
see report
Anonymous Report 1 on 20231218 (Invited Report)
 Cite as: Anonymous, Report on arXiv:2212.09693v4, delivered 20231218, doi: 10.21468/SciPost.Report.8304
Strengths
1) Accurate and stateofthe art numerics using various advanced methods
2) Informative and well thoughtthrough plots
3) Interesting question with convincing answer
Weaknesses
1) unnecessarily hard to understand for nonexperts since sentences are often too long and convoluted
Report
I have worked through the previous referee reports and the respective answers and in my opinion you, the authors, have addressed the points raised by the previous referees in a satisfactory manner. In particular, the debated formation of a local moment is clearly observed since the local spin operator becomes an (approximate) constant of motion. This is convincingly demonstrated by the dominance of the diagonal contribution in the respective SBE channel.
a) My main concern is about the convoluted and often hardtounderstand way to discuss the results. I urge you to go over the text again, try to sharpen the message and shorten sentences whereever possible. Here is an example from the bottom of page 12:
"At the same time, the lowfrequency bubble contribution is, as expected, almost vanishing in the HA, reflecting the groundstate, while it is sizable (though, generally smaller than in the perturbative regime) in the AIM or the DMFT solution of the Hubbard model, consistent with the Fermiliquid nature of their ground states."
Here is another example from page 5:
"The strong suppression of the diagonal entries in \chi, even down to
negative values, drives the breakdown of the selfconsistent perturbative description, as it is responsible for several sign flips (from positive to negative) of the eigenvalues of the generalized susceptibility and, hence, for corresponding divergences of irreducible vertex functions in the corresponding channel [1]. "
Similar constructions appear throughout the manuscript. You need to avoid this type of convoluted sentences if you want to give nonexperts a chance to understand your paper.
In the following I list a number of minor points that should be addressed to further improve the readability of the manuscript.
b) On page 14, please explain why you qualify a Green's function with a suppressed low frequency part as "insulating".
c) Please make it very clear where in the manuscript you talk about ONEline and TWOline irreducible vertices.
d) You talk about the sign change of the eigenvalues of \chi many times. Why not showing these eigenvalues as a function of T for the models considered if it is so important. Also the need to invert the generalized susceptibility could be explained for completeness. Yes, you may use formulas here!
e) You put yourself in the comfortable position that you only study models for which you can obtain exact correlation functions. But what are the implications of your findings for method development aiming at other more realistic models?
f) There is a typo in Fig. 11 (label "perTurbative") and in the end of the caption of Fig. 14.
Requested changes
see report
Author: Severino Adler on 20240117 [id 4249]
(in reply to Report 1 on 20231218)
We thank the Referee for the careful review of our manuscript, for her/his overall positive evaluation, as well as for her/his constructive observations. Below, we detail our Reply to all specific points raised in her/his report:
a) My main concern is about the convoluted and often hardtounderstand way to discuss the results. I urge you to go over the text again, try to sharpen the message and shorten sentences whereever possible. Here is an example from the bottom of page 12: "At the same time, the lowfrequency bubble contribution is, as expected, almost vanishing in the HA, reflecting the groundstate, while it is sizable (though, generally smaller than in the perturbative regime) in the AIM or the DMFT solution of the Hubbard model, consistent with the Fermiliquid nature of their ground states." Here is another example from page 5: "The strong suppression of the diagonal entries in $\chi$, even down to negative values, drives the breakdown of the selfconsistent perturbative description, as it is responsible for several sign flips (from positive to negative) of the eigenvalues of the generalized susceptibility and, hence, for corresponding divergences of irreducible vertex functions in the corresponding channel [1]. " Similar constructions appear throughout the manuscript. You need to avoid this type of convoluted sentences if you want to give nonexperts a chance to understand your paper. In the following I list a number of minor points that should be addressed to further improve the readability of the manuscript.
We thank the Referee for this comment. Following her/his suggestion we have now simplified long and convoluted sentences throughout the whole manuscript, improving the readability of our discussions. (s. enclosed .pdf file, with the highlighted changes in the text)
b) On page 14, please explain why you qualify a Green's function with a suppressed low frequency part as "insulating".
For $\text{i}\nu\rightarrow0$ the imaginary Green's function is directly connected to the spectral function at the Fermi level (namely: $\frac{1}{\pi} G(i\nu \rightarrow 0) = A(\omega=0)$). Therefore, if its value tends to zero, the spectral function at the Fermi level is correspondingly suppressed, featuring an insulating behavior.
c) Please make it very clear where in the manuscript you talk about ONEline and TWOline irreducible vertices.
We thank the Referee for his observation: The adjective ``irreducible" was indeed used in two different contexts, though always at the level of twoparticle quantities: (i) twoparticle vertex functions irreducible w.r.t. a cut of two fermionic lines (mostly in a given channel), which define the typical kernel of parquet/BetheSalpeter equations, and (ii) twoparticle vertex functions irreducible w.r.t. a cut of a single interaction line ($U$irreducibility) which naturally appear in the SBE decomposition. Following the observation of the Referee, in the revised manuscript we now always explicitly specify whether we are referring to twoparticle irreducible quantities (i.e. to twoparticle diagrams that do not fall apart if two fermionic lines are cut) or to Uirreducible quantities (i.e. to twoparticle diagrams that do not fall apart when a single interaction line is cut). At the same time, note that we have not used in the paper the word oneparticle irreducible for the oneparticle quantities, such as the selfenergy, as this diagrammatic specification was not needed for our discussions.
d) You talk about the sign change of the eigenvalues of $\chi$ many times. Why not showing these eigenvalues as a function of T for the models considered if it is so important. Also the need to invert the generalized susceptibility could be explained for completeness. Yes, you may use formulas here!
We did not show the results for the signchanges of the eigenvalues of the generalized charge susceptibility as extensive studies of this specific issue has already been performed in the literature for these models, e.g. in Ref. [1,3,7,9,20,21,2325,28] mentioned in the Introduction and in Sec. II of the revised manuscript. On the other hand, we agree with the Referee that showing the mathematical expression linking the vanishing eigenvalues of the generalized susceptibility and the divergences of the twoparticle vertex functions could definitely help the clarity of our discussions. Hence, in the revised manuscript, namely at p. 14, we have now reported the explicit definition of the twoparticle irreducible vertex function, given by the inversion of the corresponding BetheSalpeter equation. From that expression, it becomes fully transparent why a zero eigenvalue of a generalized susceptibility triggers a divergence in the twoparticle irreducible vertex of the corresponding sector.
e) You put yourself in the comfortable position that you only study models for which you can obtain exact correlation functions. But what are the implications of your findings for method development aiming at other more realistic models?
While the numerical/analytical results of our manuscript rigorously hold for the specific models considered, some general considerations can be plausibly drawn for more general cases, where no (numerically) exact solution is available, and may be used as a guidance for future studies beyond the framework of our work. Indeed, the SBE decomposition is quite general and strong spin fluctuations could manifest themselves in a nonlocal generalized charge susceptibility similarly as in the (purely local) cases studied here. Notice, however, that in the most general case, not only the physical susceptibilities, but also the Hedin vertices are nonlocal quantities. Therefore, it is not easy to predict how local or nonlocal spin fluctuations may be related to corresponding local and/or nonlocal divergences of the twoparticle irreducible vertices. These questions are certainly relevant for future investigations, as similar twoparticle irreducible vertex divergences in the charge channel have been reported to occur in the cluster DMFT (approximated) solution of the twodimensional Hubbard model even at lower Uvalues than those found in the purely local DMFT solution, cf. Refs. [9,16]. Interestingly, these divergences, beyond affecting the lowest frequency structure of the vertex functions, are additionally associated to the specific momentum structure of (commensurate) AF fluctuations (Ref. [16]). One could speculate, thus, that the highly nonperturbative physics of the twodimensional Hubbard model could be somewhat encoded in a ``generalization" of the strong link between different scattering sectors, which we discuss in our work, beyond the purely onsite physics. In that case, the role of the static local response might be replaced by strong antiferromagnetic and/or RVB fluctuations and could affect, depending on the parameters, the CDW and/or $d$wave superconduncting response of the systems. Future studies could shed light on these (at the moment purely speculative) ideas.
f) There is a typo in Fig. 11 (label "perTurbative") and in the end of the caption of Fig. 14
We thank the referee for noticing it and have fixed the typo in the revised manuscript.
Author: Severino Adler on 20240117 [id 4250]
(in reply to Report 2 on 20240106)We thank the Referee for reviewing our manuscript as well as for her/his overall positive assessment on our work, consistent to the positive evaluation of Referee 1. The Referee has asked us to consider specific points to be addressed prior to publications. We have thoroughly considered all of them (see the detailed reply below) and included the corresponding changes into the manuscript text and the appendix.
In the view of our study, some formal aspects, such as the distinction between diagonal and offdiagonal entries of the generalized susceptibilities, are important mostly for their effects in driving the signchange of the eigenvalues and the associated irreducible vertex divergences (whose link with the structure of the generalized susceptibility have been now more precisely highlighted in the last paragraph of p. 12 of the revised manuscript). Further, as the Referee correctly notes our analysis, based on the SBE decomposition of the generalized susceptibility, also allows to ascribe specific contributions to its diagonal, e.g. those associated to the exchange of a spin collective mode, to the static spin response of the system, which becomes the dominating damping factor of the onsite charge fluctuation in the local moment regime (cf. detailed discussion in Sec. 3.2 and in previous replies). At the same time, one cannot associate, in general, the whole diagonal/offdiagonal entries to one specific static/dynamic response of the systems, because more than one SBE/multiboson term, as well as the bubbleterm for the diagonal entries, might contribute to the different matrix diagonal/offdiagonal elements. In this context, in order to highlight the general properties, which characterize the frequency diagonal/off diagonal entries of the generalized charge susceptibility and, in particular, their relation with the bubble and the vertex correction terms, we have slightly extended the discussion right after Eq. (8) in Sec. 2 of the revised manuscript. As for the physical interpretation of the eigenvalues, apart from the suppression effects of onsite fluctuations we mention in our work, which has been extensively discussed, at the level of a significant ``observation'' in the recent literature, see Refs. [3, 25] in the revised manuscript, one should also recall the crucial effect that negative eigenvalues can have in lattice model, when lifting the condition of perfect particlehole symmetry: The lowest (negative) eigenvalue can trigger the onset of a phaseseparation instability, such that observed in the proximity of the Mott metalinsulator transition in the DMFT solution of the Hubbard model. Motivated by the Referee's question, we have inserted an extra sentence at p. 6 before Eq. (8), mentioning this important eigenvalue property and including the specific reference to the papers discussing this aspect (e.g. Ref. [21] and [24]). Eventually, we agree with the Referee that it is also important to explain to the readership, why one might actually need to invert the BetheSalpeter equation (and hence the generalized susceptibility matrix). In fact, as we now emphasize in Sec. 2.3 at the end of p. 5 of the revised manuscript, computing such inversions is necessary, e.g., to provide the input for diagrammatic approaches based on parquet equations, such as D$\Gamma$A and QUADRILEX. In fact, at intermediatetostrong coupling, these approaches might be directly affected by the divergences of the twoparticle irreducible vertex functions, triggered by vanishing eigenvalues. In this respect, we also mention in the revised manuscript, that for related reasons, other algorithmic schemes, such as bold diagrammatic Monte Carlo and the nested cluster approaches can face problems in the same parameter regime, as it is discussed in Refs. [2,9].
In the light of this Referee suggestion, we now use the full names also in Sec. 2.2 of the revised manuscript and define the acronyms there for a second time.
We thank to Referee for this comment, which has also been mentioned by Referee 1. We have worked throughout the whole manuscript and, whenever possible, we have simplified the text, by shortening and/or rewording too long and complex sentences, without compromising the completeness of information (s. enclosed .pdf file, with the highlighted changes in the text)
We agree with the Referee that the manuscript is long, and can thus benefit from the proposed rearrangement. Hence, we have followed her/his advise by moving Sec. 3.4 into a newly created appendix. In order to keep reference of most important information in the main part of the manuscript we have added a concise paragraph at the end of Sec. 3.3, which summarizes the results obtained via different approximations, with reference to the appendix.
Attachment:
manuscript_changes_5BFxpca.pdf