AnisotropicsofteningofmagneticexcitationsinlightlyelectrondopedSr IrO 2 4 X. Liu,1,2,3,∗ M. P. M. Dean,3 Z. Y. Meng,1 M. H. Upton,4 T. Qi,5 T. Gog,4 Y. Cao,3 J. Q. Lin,1 D. Meyers,3 H. Ding,1,2 G. Cao,5 and J. P. Hill3 1BeijingNationalLaboratoryforCondensedMatterPhysicsandInstituteofPhysics,ChineseAcademyofSciences,Beijing100190,China 2Collaborative Innovation Center of Quantum Matter, Beijing, China 3Condensed Matter Physics and Materials Science Department, Brookhaven National Laboratory, Upton, New York 11973, USA 4Advanced Photon Source, Argonne National Laboratory, Argonne, Illinois 60439, USA 5Center for Advanced Materials, Department of Physics and Astronomy, University of Kentucky, Lexington, Kentucky 40506, USA 6 (Dated:May30,2016) 1 Themagneticexcitationsinelectrondoped(Sr La ) IrO withx = 0.03weremeasuredusingresonant 0 1−x x 2 4 inelasticX-rayscatteringattheIrL -edge. Althoughmuchbroadened,welldefineddispersivemagneticex- 2 3 citationswereobserved. Comparingwiththemagneticdispersionfromtheundopedcompound,theevolution y ofthemagneticexcitationsupondopingishighlyanisotropic. Alongtheanti-nodaldirection,thedispersionis a almostintact. Ontheotherhand,themagneticexcitationsalongthenodaldirectionshowsignificantsoftening. M Theseresultsestablishthepresenceofstrongmagneticcorrelationsinelectrondoped(Sr La ) IrO with 1−x x 2 4 closeanalogiestotheholedopedcuprates,furthermotivatingthesearchforhightemperaturesuperconductivity 7 inthissystem. 2 PACSnumbers:71.27.+a,74.25.Ha,78.70.Dm ] l e - Together with the tremendous research activity on the su- energy electron evolution as observed in the cuprates. An r t perconductingcuprates[1,2],effortstocomparethecuprates anisotropic pseudo gap opens on the Fermi surface, with the s . with other related systems have also been on-going for samesymmetryasthatofthecuprates. However, therelated t a decades. Such comparison serves as a natural approach to question of whether the magnetic correlations that are often m clarifytherolesofmultipleemergentphenomenainthephase implicated in HTS [13, 14, 18] are also analogous remains - diagramofthecuprates, includingmagneticfluctuations, su- largelyunexplored. d perconductivity, pseudo gap and charge density waves etc. Here we measure the magnetic excitations in electron n The5doxideSr IrO isanexcellentcandidateforsuchstudy. doped(Sr La ) IrO withx = 0.03usingIrL edgeres- o 2 4 1−x x 2 4 3 c This so called spin-orbit-coupling driven Mott insulator [3] onant inelastic X-ray scattering (RIXS). The samples at this [ is in close proximity to the single layered cuprate La CuO , dopingshowweaklymetallicbehaviorandastronglyreduced 2 4 both structure-wise [4] and electronically [5–8]. Sr IrO magnetic moment induced by field [15, 19]. In the ARPES 2 2 4 v hosts a single hole in the t2g manifold where a Mott gap measurements,asizableFermisurface[7,17]wasobserved. 2 is opened, assisted by strong spin-orbit coupling [3, 9, 10], The phase diagram summarized by X. Chen et al. [16] sug- 7 and its magnetic excitation spectrum can be well described geststhatthisdopingsitsattheborderofthemetal-insulator 1 using a Heisenberg model of effective spin-1 moments on a transition. Weobservethatdopinginducesdampedmagnetic 2 2 square lattice [11, 12]. With a minimum single band model, excitationswithastronglyanisotropicsofteningcomparedto 0 . Sr2IrO4 and La2CuO4 are strikingly similar [13], leading to the undoped compound. Along the anti-nodal direction, the 1 the proposal that this compound could also host unconven- dispersionisalmostintactwithdoping.Alongthenodaldirec- 0 6 tional high temperature superconductivity (HTS) upon dop- tion, however, themagneticexcitationsarestronglysoftened 1 ing[13,14]. Moreover,duetotheoppositesignsofthenext- byabout 27%. Thisphenomenologyiscloselyanaologousto : nearest-neighborhoppingintegralinthesetwosystems,theo- theholedopedcupratesandsupportstheoreticalnotionsthat v i reticalworkfurthersuggeststhattheelectrondopedSr2IrO4 electron-dopingiridatesmayinduceHTS[13,14]. X mightbemorecloselyanalogoustothoseofhole(ratherthan Thehighqualitysinglecrystalsusedweresynthesizedus- r electron)dopedcuprates[13,14]. ingaself-fluxtechnique[15].IrL -edgeRIXSmeasurements a 3 AlthoughthephasediagramofdopedSr IrO hasnotbeen were carried out at 27ID and 30ID of the Advanced Photon 2 4 fully explored, a large amount of experimental work sup- Source. The incident X-ray energy was set to 11.215 keV ports the hypothesis that the Fermiology of the doped iri- basedonaresonantconditionsurvey. Theoverallexperimen- datesiscloselyanalogoustothecuprates. Upondopingwith tal resolution was about 37 meV (FWHM) estimated from a La up to 6%, Sr IrO evolves from an antiferromagnetically Lorentzianfittingofthequasi-elasticline. Alldatapresented 2 4 ordered insulator to a paramagnetic [15] or percolative [16] werecollectedatlowtemperatureof20K.Forconvenience, metal.AT-linearresistivitywasobservedwithpotassiumsub- the reciprocal space here is indexed using I4/mmm nota- stitution [15]. Further, angle-resolved photoemission spec- tion,wherethe[1,0,0]directionisparalleltothenearestIr-Ir troscopy (ARPES) data from several groups [5–7, 17] has bond directions and, in analogy to the notation used in elec- shown convincingly that doping indeed drives a similar low tronic structure measurements [5–7, 20] we refer to this as 2 theanti-nodaldirection, whereastheBrillouinzonediagonal whose width was pre-determined as 37 meV FWHM from [1,1,0] is referred to as the nodal direction. The insensitiv- an off-resonant spectrum and kept fixed during the fitting. ityoftheRIXSspectratotheout-of-planedirectionwascon- The upturn on the high energy loss side was accounted firmed at multiple Q points with different out-of-plane mo- for with a Gaussian tail. For the magnetic excitation, we mentumtransfervalues,consistentwiththe2Dnatureofthis model the imaginary part of the magnetic susceptibility with material. a Lorentzian and note that the measured RIXS intensity, Figure 1 summarizes the RIXS spectra taken at different I(Q,ω), is proportional to this multiplied by the Bose- Q points along several cuts in reciprocal space. In addition Einsteinfactor, tothequasi-elasticlinescenteredatzeroenergyloss,twodis- Γ 1 persivefeaturescanbeobservedbelow1eV.Suchanobserva- I(Q,ω)∝ Q · , (1) tionissimilartothatontheundopedcompound[11,21]. The (ω−ωQ)2+Γ2Q 1−e−ω/kBT higherenergydispersivefeature(centeredaround0.7eV)was whereω andΓ aretheenergyandFWHMofthemagnetic Q Q suggested to originate from the so called spin-orbit exciton, excitation at Q respectively, and k is the Boltzmann con- B andthelowerenergydispersionwasassignedasmagneticin stant.ThisformwasnumericallyconvolvedwithaLorentzian nature. The magnetic dispersion maximum reaches approx- in order to account for the experimental resolution. Simi- imately 200 meV at the (π,0) point, similar to the magnon larformalismwassuccessfullyappliedtothedopedcuprates bandwidth obtained on the undoped sample. The relative [22,23,25]. strength of the magnetic excitations to the spin-orbit exciton excitationsobservedhereissimilartothatobservedintheun- dopedcompoundwiththesameexperimentalgeometry.Since thecreationoftheexcitonisexpectedtoscalewiththeunoc- cupiedt states,thedopingcouldresultinnomorethan6% 2g reduction of its spectral weight. Based on this argument, we concludethatthespectralweightofthemagneticexcitations arenotsignificantlysuppressedupondoping. (π,0) ) t ni (0,0) u . b ar (π/2,π/2) ( y t si n (π,π) e t n I (π,0) (π/2,π/2) FIG.2:Acomparisonofthedispersionofthemagneticexcitationsin (Sr La ) IrO withx = 0.03withtheundopedcompound. Top 1−x x 2 4 −0.2 0 0.2 0.4 0.6 0.8 1 pannel:thesolidsquaresanddiamondsarefromthecurrentworkon Energy loss (eV) thedoped andundoped samples, and theopen circlesare fromthe undoped compound extracted from Ref. [11]. Bottom pannel: the energy-momentumintensitymapofthemeasuredRIXSspectrawith FIG.1: RIXSenergylossspectrarecordedalongseveralcutsinre- thequasi-elasticlinesubtracted; thewhitedots(Ref.[11])andblack ciprocal space. A damped magnon excitation is observed dispers- diamonds(oursample)arefortheundopedcompound. ingfromlowenergyat(0,0)and(π,π)tohigherenergiesaround thezoneboundaryat(π,0)and(π/2,π/2). Thespin-orbitexciton modeisseenaround0.7eV. The obtained dispersion relation is shown in the top panel of Fig .2 (solid red squares), together with that from the un- The dispersion of the observed magnetic excitations was dopedcompound.Theopencirclesarethedataextractedfrom extracted by fitting the RIXS spectra up to 0.45 eV energy Ref. [11], and the black diamonds are the magnon energies loss. The quasi-elastic line was fitted with a Lorentzian, measuredonourownundopedsample[19].Upon3%doping, 3 the magnetic excitations in (Sr La ) IrO are modified in thecenterofthemagneticexcitationat(π,0)agreeswellwith 1−x x 2 4 astronglyanisotropicmanner.Alongthe(0,0)→(π,0)anti- thatoftheundopedcompound.For(π/2,π/2),thelargesoft- nodaldirection,themagneticdispersioninthedopedsample ening is clearly observable from the raw experimental data. followsthatoftheundopedcompoundclosely.Themaximum Our fitting gave a peak center of 77(±6) meV, compared to at (π,0) zone boundary determines a dispersion bandwidth 110meVfromRef. [11]and105meVfromourownundoped of ∼ 200 meV. In contrast, there is a significant “softening” sample. Suchanisotropicimpactfromdopingalsoappearsin along the (0,0) → (π,π) nodal direction with doping. At thebroadeningofthemagneticexcitations. Thefittedwidths the (π/2,π/2) zone boundary point, the magnetic excitation (FWHM)are155and214meVfor(π,0)and(π/2,π/2)re- is softened by about 27%. Such anisotropic response can be spectively. Themorebroadenedmagneticexcitationssuggest seenmoredirectlyintheenergy-momentumintensitymapin strongerscatteringalongthenodaldirection. the bottom panel of Fig .2 where the measured RIXS spec- The anisotropic response to doping implies that our ob- trawithonlythequasi-elasticlinesubtractedareshown. The served magnetic excitations are not from macroscopic phase comparisonpresentedinFig.2leadstoourtwomajorobser- separation, which might be a general concern in doped sys- vations. Firstly, the magnetic excitations from the isospin-1 tems. Rather,theintroducedcarriersmusteitherco-existwith 2 square lattice in the undoped compound evolve, but clearly orresideinnano-scaleproximitytotheobservedmagnetism. survive upon 3% doping where a sizable Fermi-surface has On the other hand, we do not rule out microscopic inhomo- well developed [7, 17]. Secondly, the response of the mag- geneity, whichhasbeencommonlyobservedinelectroncor- neticexcitationstodopingishighlyanisotropic. relatedsystemsupondoping[26]. Indeed,withscanningtun- Moredetailsofsuchanisotropicimpactfromdopingcanbe nelingmicroscopymeasurementsontheir5%dopedsample, seen in Fig. 3 where the RIXS spectra for two characteristic X.Chenetal. [16]observednano-scaleinhomogeneousdis- zoneboundarypoints,namely(π/2,π/2)and(π,0),areplot- tribution of the local density of states, which is further sup- ted.Thedatawasfittedwiththeschemedescribedearlier.The portedbytheworkbyY.J.Yanetal. [27]. Similarly, nano- fasterdownturnacrosszerotowardsthenegativeenergyloss scale inhomogeneity has also been observed in hole-doped sideofthefittedmagneticexcitationpeaks(redsolidlines)is cuprates[26],andmightbeanintrinsicfeatureofthesedoped duetotheBose-Einsteinfactor. Forcomparison,themagnon Mottinsulators. energies of our undoped sample are indicated with vertical Thedopingevolutionofthemagneticexcitationswereport dashedlines. here shares much in common with that has been observed in the hole doped cuprates. It has been shown [22–25] that 400 the magnetic excitations survive into the heavily overdoped ( π, 0 ) region in many cuprate families with both dispersions and 350 (π/2, π/2) spectralweightssimilartothoseofundopedcompounds. For 300 Sr IrO , theARPESmeasurements[7,17]showthatdoping ) 2 4 counts250 msigenVifircaanngtely. rAetfotrhmesdothpeinJgeflfev=el s21tubdaienddhinerae,ftehwe Jheufnfd=red21s y (200 bandhasalreadyextendedacrosstheFermienergy,resulting sit inasizableFermisurface[17].Theindirectchargegapseems n150 e tobeclosedwiththeappearanceoftheelectron-likepockets nt I100 around (π/2,π/2) and hole-like pockets around (π,0). To explainsuchadrasticresponse,A.delaTorreetal. [7]sug- 50 geststhatdopingstronglyweakenstheonsiteCoulombrepul- sion,leaving(Sr La ) IrO tobeintheweaklyinteracting 0 1−x x 2 4 −0.1 0 0.1 0.2 0.3 0.4 region. Such strong doping dependence of the Coulomb U Energy loss(eV) was also discussed for cuprates [28], which was challenged bymanyothers[29]. Ourresultsshowthat,althoughthedop- FIG. 3: RIXS spectra for (π/2,π/2) (open diamonds) and (π,0) ing significantly renormalizes the electron band structure in (open circles) respectively. Dashed lines through the data points (Sr La ) IrO , thesystemstillsupportsmagneticcorrela- showourfittingresults. Thetwopeaksdepictedbysolidredlines 1−x x 2 4 arethefittedmagneticexcitationpeaks,withhorizontalbarsshowing tionswithsimilardispersionandspectralintensity,indicating theirFWHMoftheLorentzianmodelinEq.1. Theverticaldashed thatithasasimilardegreeofelectron-electroncorrelationas linesindicatethemagnonenergiesatthesetwoQpointsinourun- theundopedcompound. dopedsample[19]. The observation of anisotropic modification to the mag- netic excitations is particularly interesting. We notice that At both Q points, the magnetic excitations are much such anisotropic softening of magnetic excitations along the broader than the instrumental resolution. This is in contrast nodal direction has been also observed with RIXS in hole- to the case of the undoped compound [11, 12, 19, 21], for doped superconducting cuprates as well [30, 31], where in whichthemagnonsareresolutionlimited. Thus,themagnetic doped Bi-based cuprates, the magnons collapse along the excitation is highly damped upon doping. From our fitting, nodal direction but persist along the anti-nodal direction in 4 the momentum space. Such behavior is captured in random phase approximation (RPA) calculations of the magnetic re- sponse in the cuprates based on itinerant quasiparticles [30– ∗ Electronicaddress:[email protected] 34]. Hence, such anisotropic softening of the magnetic ex- [1] P.A.Lee,N.NagaosaandX.-G.Wen,Rev.ofMod.Phys.78, citationsindopedcuprateswasassociatedwiththeemergent 17(2006). correlated itinerant, or partially itinerant, nature of the elec- [2] B. Keimer, S. A. Kivelson, M. R. Norman, S. Uchida and J. tronsinthenodaldirectioncloseto(π/2,π/2)(i.e.,theemer- Zaanen,Nature518,179(2015). gent Fermi arc). The spin fluctuations in the nodal direction [3] B. J. Kim, Hosub Jin, S. J. Moon, J.-Y. Kim, B.-G. Park, C. stronglydecayintotheemergentelectron-holecontinuumand S. Leem, Jaejun Yu, T.W. Noh, C. Kim, S.-J. Oh, J.-H. Park, V. Durairaj, G. Cao and E. Rotenberg, Phys. Rev. Lett. 101, arehencesoftenedanddamped. 076402(2008). Interestingly, such a picture is somehow more appropri- [4] M. K. Crawford, M. A. Subramanian, R. L. Harlow, J. A. ate in the context of doped Sr2IrO4 than that in the doped Fernandez-Baca, Z. R. Wang and D. C. Johnston, Phys. Rev. cuprates. TheCoulombrepulsionintheIr5dshellisweaker B49,9198(1994). thanthatintheCu3dshell,andARPESexperiments[7,17] [5] Y.K.Kim,O.Krupin,J.D.Denlinger,A.Bostwick,E.Roten- havereportedthattheSr IrO ismuchmoresensitivetowards berg,Q.Zhao,J.F.Mitchell,J.W.AllenandB.J.Kim,Science 2 4 345,187(2014). doping than the cuprates. In fact, a low doping already gen- [6] Y. K. Kim, N. H. Sung, J. D. Denlinger and B. J. Kim, Nat. eratesaFermisurfacewithcoherentnodalexcitationsandan Phys.doi:10.1038/nphys3503(2015). anti-nodalpseudogap. Hence, ifitisthecouplingtothecor- [7] A.delaTorre,S.McKeownWalker,F.Y.Bruno,S.Ricco´,Z. related itinerant quasiparticles (electron-hole continuum as- Wang,I.GutierrezLezama,G.Scheerer,G.Giriat,D.Jaccard, sociated with the emergent Fermi surface) that give rise to C.Berthod,T.K.Kim,M.Hoesch,E.C.Hunter,R.S.Perry,A. the anisotropic softening of magnetic excitations along the TamaiandF.Baumberger,Phys.Rev.Lett.115,176402(2015). nodaldirectioninunderdopedcurprates,thiscouplingshould [8] Y.Cao,Q.Wang,J.A.Waugh,T.J.Reber,H.Li,X.Zhou,S. Parham, N.C.Plumb, E.Rotenberg, A.Bostwick, J.D.Den- be stronger and more likely to happen in doped Sr IrO , as 2 4 linger,T.Qi,M.A.Hermele,G.CaoandD.S.Dessau,Nature the itinerant nature of the lightly doped Ir 5d electrons are Communications7,11367(2016). more pronounced, and indeed, this is what has been clearly [9] B.J.Kim,H.Ohsumi,T.Komesu,S.Sakai,T.Morita,H.Tak- observedintheRIXSdatapresentedinthiswork. agiandT.Arima,Science323,1329(2009). In conclusion, we have shown that 3% doped [10] H.Jin,H.Jeong,T.OzakiandJ.Yu,Phys.Rev.B80,075112 (Sr La ) IrO still hosts persistent magnetic excitations (2009). 1−x x 2 4 with similar spectral intensity to the undoped compound, [11] JunghoKim,D.Casa,M.H.Upton,T.Gog,Young-JuneKim, J.F. Mitchell, M. van Veenendaal, M. Daghofer, J. van den despitethefactthatasizableFermiSurfacehasformedatthis Brink, G. Khaliullin and B. J. Kim, Phys. Rev. Lett. 108, dopinglevel[7,17]. Atthesametime,theobservedmagnetic 177003(2012). excitations respond to doping anisotropically. Along the [12] LuukJ.P.Ament,G.KhaliullinandJ.vandenBrink,Phys.Rev. anti-nodal direction, the dispersion is almost intact with B84,020403(R)(2011). doping. Ontheotherhand,alongthenodalnextnearestIr-Ir [13] F.WangandT.SenthilPhys.Rev.Lett.106,136402(2011). direction, the magnetic excitations are strongly softened and [14] Z. Y. Meng, Y. B. Kim and H. Y. Kee, Phys. Rev. Lett. 113, more damped. These observations unambiguously point to 177003(2014). [15] M.Ge,T.F.Qi,O.B.Korneta,D.E.DeLong,P.Schlottmann, a strong coupling between the doped electrons and the mag- W.P.CrummettandG.Cao,Phys.Rev.B84,100402(2011). netism in this system. Such a behavior is closely analogous [16] X. Chen, T. Hogan, D. Walkup, W. Zhou, M. Pokharel, M. tohole-dopedcuprates,furthermotivatingthesearchforhigh Yao,W.Tian,T.Z.Ward,Y.Zhao,D.Parshall,C.Opeil,J.W. temperaturesuperconductivityinthissystem. Lynn,V.MadhavanandS.D.Wilson,Phys.Rev.B92,075125 Note added: During our submission, we noticed a similar (2015). observationreportedbyH.Gretarssonandcollaborators[35]. [17] YueCaoetal.,privatecommunication;inpreparation. [18] D.J.Scalapino,Rev.Mod.Phys.84,1383(2012). WethankT.Xiangforfruitfuldiscussions. X.L.acknowl- [19] seesupplementaryinformationfordetails. edges financial support from MOST (No. 2015CB921302) [20] A.Damascelli,Z.HussainandZ.-X.ShenRev.Mod.Phys.75, and CAS (Grant No: XDB07020200) of China. Z. Y. M. 473(2003). is supported by the National Natural Science Foundation of [21] Jungho Kim, M. Daghofer, A.H. Said, T.Gog, J. van den China (Grant Nos. 11421092 and 11574359). Both X. L. Brink, G. Khaliullin and B.J. Kim, Nat. Comms. DOI: andZ.Y.M.aresupportedbytheNationalThousandYoung- 10.1038/ncomms5453(2014). TalentsProgramofChina. TheworkatBrookhavenwassup- [22] M.LeTacon, G.Ghiringhelli, J.Chaloupka, M.MorettiSala, V.Hinkov,M.W.Haverkort,M.Minola,M.Bakr,K.J.Zhou, portedbytheU.S.DepartmentofEnergy,DivisionofMateri- S. Blanco-Canosa, C. Monney, Y. T. Song, G. L. Sun, C. T. alsScience,underContractNo. DE-AC02-98CH10886. The Lin,G.M.DeLuca,M.Salluzzo,G.Khaliullin,T.Schmitt,L. work at the University of Kentucky was supported by NSF BraicovichandB.Keimer,Nat.Phys.7725(2011). through Grant DMR-1265162. Use of the Advanced Photon [23] M.P.M.Dean,G.Dellea,R.S.Springell,F.Yakhou-Harris,K. SourcewassupportedbytheU.S.DepartmentofEnergy,Of- Kummer,N.B.Brookes,X.Liu,Y-J.Sun,J.Strle,T.Schmitt, ficeofScience,OfficeofBasicEnergySciences,underCon- L. Braicovich, G. Ghiringhelli, I. Boz˘ovic´ and J. P. Hill, Nat. tractNo. DE-AC02-06CH11357. Mat.12,1019(2013). 5 [24] M.P.M.Dean,A.J.A.James,R.S.Springell,X.Liu,C.Mon- [30] M.Guarise,B.DallaPiazza,H.Berger,E.Giannini,T.Schmitt, ney, K.J.Zhou, R.M.Konik, J.S.Wen, Z.J.Xu, G.D.Gu, H.M.Rønnow,G.A.Sawatzky,J.vandenBrink,D.Altenfeld, V. N. Strocov, T. Schmitt and J. P. Hill, Phys. Rev. Lett. 110, I.EreminandM.Grioni,Nat.Commun.55760(2014). 147001(2013). [31] M. P. M. Dean, A. J. A. James, A. C.Walters, V. Bisogni, I. [25] M.P.M.Dean,JournalofMagnetismandMagneticMaterials Jarrige, M.Hu¨cker, E.Giannini, M.Fujita, J.Pelliciari, Y.B. 376,3(2015). Huang,R.M.Konik,T.SchmittandJ.P.Hill,Phys.RevB90, [26] K.McElroy,JinhoLee,J.A.Slezak,D.-H.Lee,H.Eisaki,S. 220506(R)(2014). UchidaandJ.C.Davis,Science309,1048(2005). [32] A.J.A.James,R.M.KonikandT.M.Rice,Phys.Rev.B86, [27] Y.J.Yan,M.Q.Ren,H.C.Xu,B.P.Xie,R.Tao,H.Y.Choi, 100508(2012). N.Lee, Y.J.Choi, T.ZhangandD.L.Feng, Phys.Rev.X5, [33] R.ZeyherandA.Greco,Phys.Rev.B87,224511(2013). 041018(2015). [34] M. V. Eremin, I. M. Shigapov and H. T. D. Thuy, Journal of [28] C.Kusko,R.S.Markiewicz,M.LindroosandA.Bansil1,Phys. Physics:CondensedMatter25,345701(2013). Rev.B66140513(2002). [35] H. Gretarsson, N. H. Sung, J. Porras, J. Bertinshaw, C. Dietl, [29] A.J.Millis,A.Zimmers,R.P.S.M.Lobo,N.Bontempsand JanA.N.Bruin,A.F.Bangura,R.Dinnebier,JunghoKim,A. C.C.Homes,Phys.Rev.B72,224517(2005);S.S.Kancharla, Al-Zein,M.MorettiSala,M.Krisch,M.LeTacon,B.Keimer B.Kyung,D.Se´ne´chal,M.Civelli,M.Capone,G.Kotliarand andB.J.Kim,etal.,arXiv:1603.07547(2016). A.-M.S.Tremblay,Phys.Rev.B77,184516(2008);andrefer- encestherein.