arXiv:0709.0902v1 [cond-mat.str-el] 6 Sep 2007 atfrtesigppaei e tt ob called be to state new a rele- is state phase RVB spingap the the that for insulating shows vant Mott work undoped present and The doped states. of family a for name under- anomaly. to areal attempts[14] theoretical this some stand been has measurements There other and ! counting hole on based expected ufc okt,wt nae htis Fermi that and area particles an with quasi pockets, for surface evidence present- remarkable 12], a oscillation cuprates[11, ing quantum underdoped recent statistics in very An- exclusion measurements iii) of i) notion and by early Haldane[15] the an inspired by in ii) is locking[7] phase, theory gap -charge spin of Our proposal characterise recent phase. to derson’s gap state reference spin new the a liquid temper- non-Fermi propose high We doped the over phase. as and such metal over incoherent phases crosses ature normal phase exotic this We how other superconductivity and to experimen- Tc phase explanation. this high new from how demands and emerges understand that to precise 13] need are also 12, there results[11, as tal time, the of its and phase spingap a lacking. of 10], is physics system, 9, neighbourhood, rich reference 8, the like’ [7, describes gas decades that ‘ideal 2 simple and last and attempts satisfactory the valiant in several of insights spite significant in pro- However, has superconducting Tc the state. high 6] of of 5, properties mechanism and the 4, superconductivity of 3, understanding deep insulators[2, state a Mott vided bond valence doped resonating in of (RVB) theory be- The has states clear. reference come new fur- for and need cuprates[1] a developments, in ther superconductivity in ad- Tc the gas high With Fermi of ideal interaction. vent attractive an weak of of instability presence superconduc- an the and as state viewed reference is is liquid tivity Pb, Fermi or a Hg like on exam- superconductors built For real for solids. theory complex BCS and ple, real shaping gas, of in Bose understanding role fundamental our ideal a played solid, have harmonic gas Fermi a ideal as such systems erence em liquid Fermi h unu pnlqi rRBsaei generic a is state RVB or liquid spin quantum The etrudrtnigo h pna hs saneed a is phase spingap the of understanding better A ref- ‘ideal’ solids, of theory quantum of days early From em e raby area sea Fermi nalwdniysi-hrelqi.Sisadcagsinter charges and Spins a liquid. holon spin-charge real density one low W adds a dopant in puzzles. hole with each statistics riddled where exclusion is phase, state, this superconducting study Tc high to spin- a ihpoete htaedffrn rmFermi from different are that properties with 3 2 Friliquid -Fermi lcrcladmgei rpriso nedpdcuprates underdoped of properties magnetic and Electrical 2 1 oetadnvlecuinstatistics exclusion novel and moment nttt fMteaia cecs ...Cmu,Chenna Campus, C.I.T. Sciences, Mathematical of Institute oo ar ihasio;ti oelk emo composit like hole this ; a with pairs holon A . 2 3 Frilqi:tesce fhg-ccuprates high-Tc of secret the liquid: -Fermi 3 2 g eni eetqatmoclaineprmnsa w diffe two at experiments oscillation quantum recent in seen , h and 1 = 2 3 ie agrthan larger times g s = 2 1 oafil nr ped ae)si iudvcu,result vacuum, liquid spin gaped) (pseudo inert fairly a to , .Baskaran G. 2 3 - g hole eni eetqatmoclainmeasurements. oscillation area quantum sea Fermi recent anomalous in an the- explains seen our it how formulating show After we ory violated. strongly system is theorem interacting non reference the em e.Ttlnme fhln scnevdbecause conserved is holons of number novel Total a condensa- into sea. and composite Fermi spinon (fermion) spinon-holon doped this one of pairing and tion as holon it view doped we one is theory, of our there In side, enters temperature locking. one spin-charge high a as from that regime, fact gap spin phenomenological the a as it we take Anderson[7], Following will paper. present the in a teraction equal as and call we holons which fluid doped , interacting doped of of fluid What number a liquid. is one spin inert have adds fairly we hole this to of doped spinon Each range one and the us. holon small in to interest negligibly gap) of (pseudo and temperature spinons holons thermal no of has spin density insulating which Mott phase, reference liquid a exists there hypothesis: c rssol rmtehl oin h - oe plus metal’. model ‘projected describes t-J repulsion this The coulomb range motion. short hole unscreened gap. the some from energy only finite arises ics dynam- a charge and have survives superexchange to Consequently continue excitations band iud n utne iud nsm udmna way. fundamental some in Further liquids Luttinger and liquids dynamically from arise to thought compos- fields. be spinon-holon gauge may RVB paired generated pairing of The liquid density low ite. A 1: FIG. uao.We oigfato sx is fraction doping When sulator. = ews osuyalgtyhl oe spin- doped hole lightly a study to wish We u hoyof theory Our 2 3 ed o oit h aueadoii fti in- this of origin and nature the into go not do We . eepana nmlu xaso of expansion anomalous an explain We . c n omanvlcletv tt,a state, collective novel a form and act 2 3 Frilqi sntaibtclycnetdto connected adiabatically not is liquid -Fermi rps oe eeec tt to state reference novel a propose e ntesi a hs,aprecursor a phase, gap spin the in doera pnn ihHaldane with spinon, real one nd 2 3 Frilqi trswt h following the with starts liquid -Fermi 0 1,India 113, 600 i are hre+ and +e charge carries e etdopings. rent ∗ Spinon << Holon ,uprHubbard upper 1, ing Luttinger . spin-charge 1 2 otin- Mott 2

π of charge conservation. However, doped spinon number ( 0 , ) ( π , π ) is not a conserved quantity in general; in fact, we ex- k σ pect them to form spin singlets and disappear into the vacuum. Our assumption is that spin-charge locking and k' σ' ( π , 0 ) coulomb repulsion between the composite keeps the spins Filled Fermi Sea of the doped spinons nearly free and uncorrelated (figure 1). Spinon and holon of RVB theory are topological excitations[5, 6]; this property manifests an underlying topological order[5]. Their addition or removal leads to FIG. 2: Two spinon states created by a Gutzwiller projected a shakeup of the many body states they live, in a global particle-hole excitation in a Fermi sea of half filled band of free fashion; for example, in 1D repulsive Hubbard model in a 2D square lattice. By construction, momentum these shakeups appear as self consistent change of all the of a given spinon is limited one half of the BZ (either inside pseudo momenta of spins and charges. We suggest that or outside the Fermi sea, which get identified after GW pro- jection). exclusion statistics, introduced by Haldane[15], is a key and minimal notion that needs to be incorporated, in or- der to solve our complex spin gap problem, where we do analysis and find number of allowed single spinon states not have the advantage of exact solutions or Bethe type in the Brillouin zone (BZ) of a spin liquid state. One di- ansatz. To characterise the spin-charge fluid we need mensional Haldane-Shastry model[16] is known to exhibit to find the exclusion statistics of the constituent spinon a remarkable property that the number of single spinon (g ) and holon (g ) separately first and then that of the s h states in k-space is only half the number of k-points in spinon-holon composite (g ). We find that exclusion hole the BZ. This implies an exclusion statistics of g = 1 for a statistics of the spinon-holon pair is a sum of their indi- s 2 spinon. We will show, using variational (Fermionic) RVB vidual exclusion statistics, g = g + g . hole s h wave functions for spinon states that allowed k-values of In Haldane’s theory of exclusion statistics, a reduction a single spinon is restricted to only half the number of in Hilbert space dimension of a given single quasi parti- points in BZ for RVB states in any dimension, imply- cle, due to presence of other quasi particles, determines ing that the mutual exclusion statistics of a spinon is 1 , g r , ..r , ..r 2 exclusion statistics parameters . Let Ψ( 1 i M ) be independent of their spin states. the wave function of a M identical quasi particle system in In RVB theory[2, 4], Gutzwiller projection and mean coordinate representation. For our purpose we consider field states (Slater determinants) play a fundamental a lattice system containing N lattice sites. We leave free role. Let us start with a simple mean field solution, a half the i-th quasi particle coordinate and freeze coordinates filled band of free Fermi sea (zero-flux solution), |FSi in 2 of the rest M-1 quasi particles at a typical configuration dimensions. A two spinon state[4], with momentum and C denoted by and expand Ψ in a set of single particle spin projections (k, σ) and (k , σ ) is defined as (figure basis wave function of the i-th particle {φC (r )} : ′ ′ ν i 2):

C C Ψ(r1, ..ri, ..rM )= X aν φν (ri) (1) ′ ′ |kσ, k′σ′i≡ PGck† σck σ¯ |FSi (2) ν=1,..Q It is easy to see that the above construction creates two The set of single particle wave functions {φC (r )} span a ν i unpaired spins with spin projections σ and σ′, in the one particle Hilbert space with dimension Q. For exclu- background of resonating valence bonds and maintaining sion statistics to become meaningful, Q should be inde- single electron occupancy at every site. As two spinons pendent of the configuration C. In general Q will change are created from a half filled band Fermi sea contain- as particles are added. Haldane defines the exclusion ing a particle-hole pair, the allowed momenta of the two statistics parameter through a suggestive differential re- spinons k and k′ are the two halves of the BZ (figure 1). lation ∆Q ≡−g ∆M. It is remarkable that in correlated Further, it is easy to show using equation (2) that (non slater determinantal) states like Laughlin quasi par- ticle states in fractional quantum Hall effect or spinon |kσ, k′σ′i ≡ |k + Q σ, k′ + Q σ′i (3) states of RVB theory g can be fractional; this is due to an intrinsic non-orthogonality of the ‘topological excita- where Q ≡ (π, π). The above completes the proof that tions’, which is transparent in real space basis. How ex- the spinon band covers only half the BZ. That is, there is clusion statistics, gets translated into k-space is a difficult one independent spinon state per two points in k-space. issue, as we will see below. This proof goes through, as long as the starting mean For holon it is easily shown that gh = 1. It follows field RVB state is a Slater determinant that fills any half from the hard core and band filling nature. To find the of mean field single particle states with up and down spin exclusion statistics of spinons, first we will offer a k-space electrons. 3

One can get a feel for the above state reduction, using a) b)             Haldane’s argument, which uses non-orthogonality of site     Spinon  Spinon π π π π localised basis for spinons. For an insulating RVB state _ , _ _ , _ Holon () 3 ( ) 22 _ −Fermi Liquid 2 2 containing N spins, Nsp spinons (unpaired spins), and 2 N Ns unbroken singlet bonds on a lattice, the Hilbert (0,0) (0,0) −  2        N Ns         space dimension of a single spinon is 1 + −2 , inde-   pendent of the spin projection of the spinon. This is Spinon Spinon understood as follows. A given spinon can occupy the site that it is initially on or it can be moved to any site that is part of a pair. But if |1{23}i represents a FIG. 3: a) spinon and holon spectrum in the π-flux (s+id) three-site wave function where site 1 is unpaired and RVB mean field solution. In a slave boson mean field theory, the doped spinons occupy and form a spinon Fermi sea and the sites 2 and 3 are paired, nonorthogonality means that 1 doped holons Bose condense at the Γ point. b) Our proposal: |1{23}i = √2 |2{13}i− |3{12}i. Hence only the combina- the spectrum of the spinon-holon paired composite and the 3 tion |2{13}i+|3{12}i is independent of |1{23}i. There is 2 -Fermi liquid they form in k-space. thus one extra independent spinon state per bond (half the number of sites the bonds occupy). This makes 1 gms = 2 , independent of their spin states. point and spinons build four identical Fermi pockets at What is the exclusion statistics of a spinon-holon com- π π the (± 2 , ± 2 ) points (Figure 3a). In the RVB gauge posite ? We answer this question by a heuristic argu- theory[10, 18, 19], which takes into account fluctuations ment. Fractional exclusion property is a manifestation of arising from double occupancy constraints, dynamically real space non-orthogonality of the quasi particles. The generated RVB gauge fields appear. In our picture, we spinon component of the composite continues to have the have matter (spinons and holons), interacting through same non-orthogonality property for every spinon-holon residual dynamically generated gauge fields. We view configuration in a composite. Thus it follows that mu- the locking phenomenon as formation of a weak, spinon- tual exclusion statistics of a composite is additive and holon , induced by dynamically generated independent of the spin states of the two : i.e., gauge fluxes (figure 1). We take the spinon-holon pairing 1 3 ghole = gs + gh = 2 +1 = 2 . That is, as we add a as a phenomenological fact and do not go into details in spinon-holon composite to the system, its holon compo- the present paper. As spinons and holons have their en- nent excludes one state and spinon component excludes π π ergy minima at (± 2 , ± 2 ) and (0, 0) in k-space, we find half a state for the next composite that enters. that the paired spinon-holon composites have their en- π π So far our discussion has been general. Now we use spe- ergy minima at the finite momenta, (± 2 , ± 2 ). cific spin liquid state relevant for 2D cuprates and get in- 3 Since there are four minima for our 2 fermions in formation about the location and volume Fermi sea pock- k-space, spinon-holon composite created by doping will 3 ets, that a collection of fermions with a 2 exclusion statis- build four Fermi pockets around these equivalent points tics will condense into. We start with the π-flux (s+id) in k-space (figure 3b). As we mentioned earlier, the mu- RVB mean field solution[17] of the spin- 1 Heisenberg 3 2 tual exclusion statistics of 2 is independent of the spin S S antiferromagnetic Hamiltonian H = J P ij i j , on a state of the holes. It means that our more exclusive square lattice with a nearest neighbour exchangeh i J. This 3 fermion occupies an average of 2 k-points, below the variational solution, after Gutzwiller projection, is known Fermi level. We incorporate this in the following fash- to give a good ground state energy. It has nodal spinon ion: each low energy k-state, below a Fermi momentum π π excitations at (± 2 , ± 2 ); that is, after linearisation, a is occupied by our fermion with an average occupancy of Dirac cone like spectrum is obtained: ǫsp(q) ≈ ¯hvspq. 2 3 , independent of the spin state. This leads to an average Here vsp ∼ Ja is the spinon Fermi velocity, a is the lat- 3 number of 2 k-state per fermion. Further, for simplicity, tice constant and q is measured from the four nodes. The we assume that the spins are dangling and form a param- linearly vanishing density of states of spinons as its en- agnetic gas in k-space. Knowing the technical difficulties ergy tends to zero ensures that we have a pseudo gap and in constructing the Bethe ansatz in 1D systems, it is clear small number thermal spinons at low temperatures. The that it will be rather hard to give a wave function rep- massive holon excitation occurs at the (0, 0) point in the 3 2 2 resentation of our 2 Fermi sea, in terms of underlying h¯ q BZ (fig 1). It has a quadratic dispersion, ǫh ≈ 2m∗ . Here electron coordinates. However, a realization that such a h¯2 fractional occupancy and novel Fermi sea is possible, is a m∗ ∼ ta2 is the holon effective mass. As discussed earlier, each added hole contributes one key step for a better understanding of the spin gap phase. π π real spinon and one real holon in their respective low en- Each Fermi sea around the (± 2 , ± 2 ) points will in- 3 xN ergy regions in k-space. In the slave boson mean field clude 2 4 k-points below its Fermi surface. In other 3 xN theory, without local constraints, holons and spinons are words the Fermi surface area of each hole pocket is 2 4 . non interacting: doped holons Bose condense at the Γ As each k-point below the Fermi sea has an average oc- 4

2 3 cupancy of 3 , we get the total number of holes as xN. the 2 Fermi sea is likely to loose part of the doped spinons Thus we see that the new ‘hole’ quasi particle, the as singlet pairs and build a strong superconducting am- weakly bound holon and spinon has no adiabatic con- plitude. Some finite temperature effects can be captured nection to ‘hole’ of the reference non interacting system by an ‘ideal gas’ like Haldane-Wu[22] distribution func- 3 having a total of N(1-x) electrons, even though they both tion of our 2 Fermi sea. We hope to address the above 1 carry charge +e and spin- 2 moment.. issues in future. To understand nature of spingap phase, devoid of complications arising from disorder, high quality single As this preprint was getting completed, a preprint by crystals of underdoped YBa Cu O [11]and two chain ∗ 2 3 6.5 P.W. Anderson [arXiv:0709.0656] has appeared, with a compound YBa Cu O [12]have been studied recently. 2 4 8 spirit similar to mine. Title of Anderson’s paper inspired These two underdoped compounds have different doping, me to replace a old title by the present one. x ≈ 0.1 and 0.125 respectively. In these studies a super- conducting state is destroyed by a large external mag- netic field, thereby accessing a spin gap normal state at sufficiently low temperatures, where inelastic mean free path of the quasi particles become sufficiently large, as [1] A. Bednorz and A. Muller, Z. Phys. B 64 189 (1986). to be able to perform complete cyclotron motions. The [2] P.W. Anderson, Science 235 1196 (1987) quantum oscillation measurements in underdoped YBCO [3] The Theory of Superconductivity in the High Tc Cuprate revealed presence of 4 closed Fermi surface pockets in the Superconductors, P.W. Anderson, (Princeton University Press, Princeton) BZ. A theoretical fit, using Fermi liquid formula, gives xN [4] G. Baskaran, Z. Zou and P.W. Anderson, Sol. State each hole pocket an area that is ≈ 1.5 4 . As observed by Commm. 63 973 (1987) the authors of the experiments, for both the compounds [5] S. Kivelson, J. Sethna and D. Rokshar, Phys. Rev. B 38 this area is nearly 1.5 times the value that one will get 8865 (1987) assuming that each added hole occupies one k-point in- [6] P. W. Anderson, G. Baskaran, T. Hsu and Z. Zou, Phys. side the Fermi sea. However, it is in excellent agreement Rev. B 58 2790 (1987) 3 xN [7] P.W. Anderson, Phys. Rev. Lett., 96 017001 (2006) with our formula Fermi pocket volume of 2 4 . C 24 Thus we believe that we have given a microscopic ex- [8] P.W. Anderson et al., J. Phys. R355 (2004) [9] G Baskaran, Iranian J. Phys. Res., 6 163 (2006) planation for the observed anomalous Fermi sea volume. [10] P. A. Lee, N. Nagaosa and X.-G. Wen, Rev. Mod. It will be very important to confirm our theory and mea- Phys. 78 17 (2006); X. -G. Wen, Quantum Field Theory sure the quasi particle Fermi surface shape and volume at of Many-Body Systems (Oxford University Press, Ox- low temperature regime, after destroying superconduc- ford)2004 tivity, using similar systems where quantum oscillation [11] Nicolas Doiron-Leyraud et al., Nature, 447 565 (2007); measurements have been performed. Existing ARPES H.D. Drew, arXiv:0708.3117 [12] E.A. Yelland et al., Cond-mat:07070057; A. F. Bangura results[20, 21] in the spin gap phase are above the super- et al., arXiv:0707.4461 conducting temperatures of the under doped regime. In [13] Yayu Wang, Lu Li and N. P. Ong, Phys. Rev. B 73 this region, a strong buildup of superconducting ampli- 024510 (2006) tude and associated vortex liquid phenomena are likely to [14] Wei-Qiang Chen et al., arXiv:0706.3556; A. Carrington 3 destroy the Fermi surface of our reference 2 Fermi liquid. and E.A. Yelland, arXiv:0706.4418 Having given a novel reference state to describe the low [15] F.D.M. Haldane, Phys. Rev. Lett., 67 937 (1991) 60 temperature spin gap phase, important tasks remain: i) [16] F.D.M. Haldane,Phys. Rev. Lett., 635 (1988); B.S. Shastry, Phys. Rev. Lett., 60 639 (1988) how superconducting instability arises in the 3 Fermi liq- 2 [17] I. Affleck and B. Marston, Phys. Rev. B 37 3774 (1988); uid at zero temperatures, as we increase doping ii) nature G. Kotliar, Phys. Rev., 37 3664 (1988); F.C. Zhang, C. of a suggested vortex liquid above superconducting Tc in Gros, T.M. Rice and H. Shiba, Supercond. Sc. and Tech., the spin gap phase. iii) nature of the quantum phase 1 36 (1988) transition to the over doped regime, where ARPES sees [18] G. Baskaran and P.W. Anderson, Phys. Rev. B 37 580 a large Fermi surface, as given by band structure results. (1988); G. Baskaran, Physica Scripta, T 27 53 (1989); 3 [19] P.B. Wiegman, Phys. Rev. Lett. 60 821 (1988) Constituents of 2 -fermion have conflicts, including ten- dency for spinon pair annihilation, implying different in- [20] J.C. Campuzano, M.R. Norman and M. Randeria, in Physics of Superconductors, vol. II, Ed. K.H. Bennemann stabilities. It is likely that at T=0, as we increase the and K.H. Ketterson (Springer, Berlin) 2004, page. 167; doping, superconducting instability arises beyond a crit- [21] A. Damascelli, Z.X. Shen and Z. Hussain, Rev. Mod. ical doping, by annihilation of doped spinons, and a con- Phys., 75473 (2003) sequent collapse of our spin-charge fluid into a charged [22] Y.S. Wu, Phys. Rev. Lett., 73 922 (1994); M.V.N. fluid. In the same way, at finite T in the spin gap phase, Murthy and R. Shankar, Phys. Rev. Lett., 72 3629 (1994)