SciPost Submission Page
Single and twoparticle observables in the Emery model: a dynamical meanfield perspective
by YiTing Tseng, Mário O. Malcolms, Henri Menke, Marcel Klett, Thomas Schäfer, Philipp Hansmann
Submission summary
Authors (as registered SciPost users):  Philipp Hansmann · Thomas Schäfer 
Submission information  

Preprint Link:  https://arxiv.org/abs/2311.09023v2 (pdf) 
Date submitted:  20240911 14:10 
Submitted by:  Hansmann, Philipp 
Submitted to:  SciPost Physics 
Ontological classification  

Academic field:  Physics 
Specialties: 

Approaches:  Theoretical, Computational 
Abstract
We investigate dynamical meanfield calculations of the threeband Emery model at the one and twoparticle level for materialrealistic parameters of high$T_c$ superconductors. Our study shows that even within dynamical meanfield theory, which accounts solely for temporal fluctuations, the intrinsic multiorbital nature of the Emery model introduces effective nonlocal correlations. These correlations lead to a nonCurielike temperature dependence of the magnetic susceptibility, consistent with nuclear magnetic resonance experiments in the pseudogap regime. By analyzing the temperature dependence of the static dynamical meanfield theory spin susceptibility, we find indications of emerging oxygencopper singlet fluctuations, explicitly captured by the model. Despite correctly describing the hallmark of the pseudogap at the twoparticle level, such as the drop in the Knight shift of nuclear magnetic resonance, dynamical meanfield theory fails to capture the spectral properties of the pseudogap.
Author indications on fulfilling journal expectations
 Provide a novel and synergetic link between different research areas.
 Open a new pathway in an existing or a new research direction, with clear potential for multipronged followup work
 Detail a groundbreaking theoretical/experimental/computational discovery
 Present a breakthrough on a previouslyidentified and longstanding research stumbling block
Current status:
Reports on this Submission
Report
The authors present the solution (oneparticle spectral function, selfenergy and uniform susceptibility) of the threeband Emery model within the singlesite DMFT approximation. In this model the electrons in the oxygen porbitals are treated as noninteracting and only the dd onsite repulsion is taken into account. The DMFT approximates the solution by neglecting dynamical nonlocal correlations, which in this case means correlations beyond the unit cell. The parameters of the Emery model are adapted to describe the cuprate superconductors as in earlier cluster DMFT studies of Refs. 34 and 35. More precisely, one set of parameters is used. The oneparticle spectra are discussed and compared directly with the oneband Hubbard model (also computed by the authors and placed in the Appendix) and indirectly with earlier cluster studies. The lack of momentum dependence in the selfenergy and hence absence of the pseudogap is addressed and at the oneparticle level the deficiencies (and also strengths) of DMFT carefully addressed. The authors then present twoparticle calculations on the example of static uniform susceptibility, computed by performing computations at finite magnetic field. The numerical results are then compared to the experimental measurements of the Knight shift and (in the Appendix) to the oneband Hubbard model in DMFT as well as to the asymmetric dimer toy model. They find that the asymmetric dimer represents the qualitative temperature dependence of the susceptibility much better than the DMFT calculation for the oneband Hubbard model, which the authors assign to the singlet formation captured well by the dimer or by the Emery model (even in DMFT), but apparently not in the DMFT for the oneband Hubbard model. The singlesite DMFT calculations are thought of as a starting point for a future D$\Gamma$A study of the Emery model.
I find already the DMFT results interesting and certainly deserving publication. What I miss in the current version of the manuscript is a thorough discussion of twoparticle properties as captured or not captured by DMFT. It is not clear to me why the uniform susceptibility is the correct quantity to represent the Knight shift in place of the local susceptibility. The authors should explain if this is related to the fact that the computational method used gives better results for uniform susceptibility, because it can be computed directly as a derivative and not from the correlation function (which I suspect to be the case), or to the actual way the experiments are performed. In case the authors address this point in detail, I recommend the paper for publication.
Minor points:
(1) I would appreciate a justification of the choice of the value of the onsite repulsion. Do the authors expect the results to depend on the precise value?
(2) Optional: A third plot in Fig. 6 with the donly spin susceptibility (if data are available for the relevant temperatures) could make the argumentation stronger.
(3) The Clogston–Jaccarino plot should be explained.
Recommendation
Ask for minor revision
Report
The authors address the physics of the threeband Emery model within the singlesite DMFT approximation. The calculations are performed for a set of parameters, which has been recently suggested as material realistic for the cuprate superconductors. As the authors state themselves, singlesite DMFT comes with limitations, foremost the failure of momentum resolving spectral functions crucially required to unravel pseudogap physics. However, this work can lay a ground on which other more advanced techniques can be rooted and should be seen as a first step towards more elaborate studies. T
Within the limitations of the singlesite DMFT approximation, the calculations appear to be competently carried out, and the results are presented concisely. I would nevertheless suggest several improvements to the readability, as in the current form, the manuscript only seems accessible to people at the close intersection between the DMFT community and the cuprate community.
1) In Fig. 1, the color code of the band structure is referred to as "orbital character". As far as I can see, no definition in the text is provided and should be added.
2) The authors repeatetly mention the ZhangRice singlet state without an explanation. The readability of the manuscript can be improved by incorporating a concise explanation of the ZhangRice singlet. Moreover, in Fig. 2 it should be explained why the zero energy peak corresponds to the Zhang Rice singlet. This is currently not explained well in the manuscript.
3) The DMFT method is explained only very briefly in section 2. This should be expanded and a discussion of strengths and limitations should be included.
4) The shown spectral functions show rather fine features, which is remarkable as they apparently have been obtained from analytical continuation. The manuscript should explain in further detail, to which accuracy the Matsubara Greens functions are calculated and how this leads to such precise predictions on the realfrequency Greens functions.
5) In general, all energies are reported in electronvolt. However, to broaden the scope for other computational studies it would not hurt to also include dimensionless units, e.g. in units of some hopping in the three band model.
6) In Fig 3., the authors mention that the selfenergy shows Fermi liquid behavior. The authors should specify explicitly what is meant by this and how this can be seen.
7) Similarly, the authors repeatedly mention nonCurie like behavior of the magnetic susceptibility. This should be made more specific, e.g. by saying that no divergence but a maximum is observed.
8) The authors show a ClogstonJaccarino plot but do not explain what this means. This plot should not be considered widespread knowledge and a short explanation would be helpful.
9) In the conclusion, the authors mention that a single band Hubbard model seems to give precise predictions for the nickelates. This statement should be weakened, as much of the community seems to agree that multiorbital effects in the nickelates are of crucial importance.
I consider the present manuscript a valuable contribution to the literature and am in favor of publication. However, I am sceptical about whether the paper meets the bar of being among the top 20 percent of PRB publications required for SciPost Physics. It is unclear which of the findings will survive a treatment with more accurate approximations beyond singlesite DMFT. Even if there is better agreement with cuprates than in the single band case, it is not excluded that this is a coincidence, even though it might be an excellent hint already. I would consider a study using cluster DMFT revealing the pseudogap behavior to meet this bar, but a singlesite DMFT calcution might just be too crude to induce confidence the the observed behavior is a trustworthy feature of the genuine twodimensional Emery model. I would, therefore, rather suggest publishing the manuscript in the SciPost Core series, upon my requested changes are implemented.
Recommendation
Accept in alternative Journal (see Report)
Strengths
1  Clear Presentation
2 The methodology is well established with its pros and cons. The manuscript leverages on this point
3 Solid information and simple message about a very important topic
Weaknesses
1 The connection with similar models (periodic Anderson model) is not discussed
Report
The manuscript presents a theoretical analysis of the Emery model for hightemperature superconductors using singlesite Dynamical MeanField Theory (DMFT). DMFT only introduces local selfenergy corrections, hence it is usually discarded as a method to study models for hightemperature superconductors, in which nonlocal correlations are universally considered fundamental. Yet the present work shows that, in contrast with the singleband Hubbard model, the threeorbital Emery model shows important nonlocal correlations owing to the multiorbital character. This point is well discussed in the manuscript in general terms and its consequences are explored. In particular, the authors find that the q=0 uniform magnetic susceptibility $\chi_{dd}$ is remarkably different from the singleband Hubbard model (shown in the appendix A) and, remarkably, it reproduces some trends of the experimental data in cuprates, reproducing in particular the qualitative behavior in the pseudo gap region. This is particularly surprising and intriguing, as the DMFT singleparticle spectra fail, as expected, to reproduce the pseudo gap. The authors discuss some arguments behind this seeming discrepancy (lack of feedback on singleparticle properties of the behavior of the twoparticle response given by the spin susceptibility) and the perspective to overcome it by using D$\Gamma$A.
I find this result very intriguing and I support the publication of the manuscript in SciPost Physics. The use of the very well established singlesite DMFT and the straightforward calculation of the uniform susceptibility computing the response to a small static magnetic field make the conclusions very solid and they identify an non trivial effect. I am also intrigued by the simple twosite model described in the appendix B. Indeed this oversimplified model captures the main trend of the DMFT calculation. This success also suggests that the results discussed here can be related to those for the periodic Anderson model (PAM). Indeed in Phys. Rev. B 85, 235110 (2012) the PAM has been solved within DMFT, finding the formation of ZhangRice singlets similar to those reported here. In that paper the authors compute the local (impurity) spin susceptibility which however has a Curielike behavior. So, I am wondering if this is the case also for the Emery model, while the nonCurie behavior is found in the q=0 response, which is notoriously different from the local counterpart in DMFT.
As I detail below, this comparison with the local susceptibility (and possibly with the PAM) could help to strengthen the manuscript.
Another suggestion I have is to consider a title/presentation which makes it more explicit that the twoparticle observables mentioned in the title coincide with the uniform susceptibility. The current version of the title gives the incorrect impression that the manuscript would show dynamical twoparticle observables which have been computed in a number of papers from the Vienna group and coworkers.
I detail below also other points that I would like the authors to address in order to meet the criteria for publication. A revised version meeting all my points should be published in SciPost Physics.
Requested changes
1 I would eliminate the sentence about static meanfield reproducing the undoped system. While an AFM is typically found, the meanfield description is quite far from the actual situation and from DMFT.
2 At line 1011 of page 2, I would stress more that the drop of the Knight shift is generally interpreted as resulting from the reduction of spectral weight, while the present calculation finds it without a spectral pseudo gap.
3 The authors should discuss what is the choice of the double counting used here and whether also the interactions are extracted from abinitio calculations. I would also mention the value of the dp repulsion even if this is neglected, since this is probably the most important simplification introduced in the model
4 The authors write that the susceptibility is computed for a value of h $\simeq$ 5 meV. Why $\simeq$ and not just equal? Do I understand correctly that the authors simply use a linear approximation connecting zero field with this value of h? Did they test if this value is small enough for the present calculations? (observing that no change is introduced by reducing h)
5 minor: I think that the ClogstoneJaccarino plot should be explained to the general readership
6 As mentioned above, I would suggest a title/presentation making more explicit that the twoparticle observables mentioned in the title coincide with the uniform susceptibility.
7 I would add if possible one example of the local spin susceptibility to compare with the q=0 response. This can also be connected to a comparison with the periodic Anderson model.
8 In Fig. 2 the x and y range of the panels (a) showing the selfenergy are cut in a rather extreme way. I would extend them to highlight better the evolution of the real and imaginary part of $\Sigma_{dd}$.
Recommendation
Publish (surpasses expectations and criteria for this Journal; among top 10%)