possible electronic phase separation ( eps ) in the cuprates has usually been found in terms of stripe phases .
thus , neutron diffraction measurements find evidence for fluctuating stripe order in la@xmath0sr@xmath1cuo@xmath2 ( lsco ) associated with incommensurate inelastic neutron scattering@xcite , which can be transformed into long range charged stripe order@xcite by co - doping with nd or eu .
similar incommensurate peaks are found in other cuprates , and long range charge ordered states are also found@xcite in strongly underdoped samples of yba@xmath3cu@xmath4o@xmath5 ( ybco ) , with stripes parallel to the chains , while at higher doping short range stripe order is found at virtually all temperatures up to the pseudogap t@xmath6 .
however , eps can also manifest itself in the form of _ larger domains _ , particularly if the dopant ions are somewhat mobile and can follow the hole motion .
these domains , long known in la@xmath3cuo@xmath7 may recently have been found in scanning tunneling microscope ( stm)@xcite and microwave@xcite studies of bi@xmath3sr@xmath3cacu@xmath3o@xmath8 ( bi2212 ) .
stripes are typically interpreted either in terms of charged domain walls in antiferromagnets or of restricted eps . while the models are virtually indistinguishible in the underdoped regime , they diverge as doping increases .
a phase separation model implies the existence of a uniform end phase : at some doping antiferromagnetic ( afm ) stripes must disappear , leaving a uniform state similar to that found on the charged stripes . for slightly lower doping the afm stripes form a minority
 a situation not possible in a domain wall picture .
( a note of caution : more recent domain wall models find a lower doping on the charged stripes , and as the doping decreases , the distinction between the two models gradually blurs . )
evidence for the termination of the afm stripes at a fixed doping@xcite and recent observations of larger domain structures@xcite taken together provide extremely strong support for a phase separation scenario .
at this point it is essential to better characterize the nature of the stripes , in particular the charged stripes , and to understand how their properties affect physical phenomena , in particular photoemission spectra.@xcite while fluctuations@xcite play an important role in the real cuprates , we have constructed ordered stripe arrays , for which detailed tight - binding calculations are possible@xcite .
such calculations can aid in elucidating the structure of the charged stripes , both for wide stripes  what stabilizes the preferred hole density  and for narrow stripes
how do quantum size effects ( qse ) modify properties of the stripes .
the present paper provides an extended analysis of these issues ; some of these results were summarized recently@xcite .
the paper is organized as follows .
section ii enumerates key issues which must be addressed in any eps model of stripes , including determination of the hole density on charged stripes .
it is found that the model can simultaneously account for a large number of experimental results .
the charged stripes are probably stabilized by some competing order , either magnetic or paramagnetic ( including charge density waves ) .
section iii shows that magnetic charged stripes can arise in a mean field hubbard model , and can be either ferromagnetic or a linear antiferromagnetic ( laf ) phase similar to the white - scalapino stripes@xcite .
calculations on stripe arrays find that the charged stripes lead to midgap states near the fermi level ; section iv shows how these data can be interpreted in terms of quantum size effect ( qse ) on single stripes .
section v presents the results of single stripe calculations : competing charge density wave ( cdw ) and superconducting orders can exist on paramagnetic stripes , but they are strongly modified by the qse .
on laf stripes , d - wave superconductivity and an unusual form of cdw are both found to persist down to the narrowest ( 2 cells wide ) stripes .
section vi includes results on stripe arrays : the laf stripes produce photoemission constant - energy maps in significantly better agreement with experiment . in the discussion , section vii ,
we summarize additional recent evidence which favors a phase separation model  in particular , evidence that afm stripes terminate slightly beyond optimal doping , where superconducting properties remain strong , and that a second regime of phase separation exists above optimal doping  and we show that experiments are consistent with our prediction that the superconducting gap grows in the underdoped regime .
conclusions are given in section viii .
a phase separation model of stripes is characterized by the two well - defined stable end phases between which phase separation takes place .
the insulating stripes are generally understood to be antiferromagnetic ( afm )  essentially the same as the mott insulator found in undoped cuprates .
the hole - doped stripes are assumed to have a finite doping , @xmath9 .
this simple idea has three experimentally verifiable characteristic features : ( i ) _ the stripe phase must terminate _
when @xmath10 ; ( ii ) there will be a _ crossover _ at a lower doping , @xmath11 from magnetic - dominated ( @xmath12 ) to charge - dominated ( @xmath13 ) stripe arrays ; ( iii ) some _ interaction _ on the charged stripes stabilizes the end phase at @xmath9 .
here we discuss current evidence for ( a ) the doping on the charged stripes , ( b ) evidence for a crossover , and ( c ) the nature of the dominant interaction stabilizing the charged stripes .
in addition , we ask what constraints the yamada plot puts on the model , and how a model based on eps compares to a domain wall model .
recent evidence suggests that the stripes and pseudogap terminate at the same doping @xmath9 while superconductivity persists to higher doping@xcite .
however , the proper choice of @xmath9 requires some discussion .
the neutron diffraction measurements of tranquada@xcite and yamada@xcite have established that charged stripes in la@xmath0sr@xmath1cuo@xmath14 ( lsco ) have an invariant topology over the doping range @xmath15 , acting as antiphase boundaries ( apbs ) for the afm stripes and having a net doping of 0.5 holes per stripe . however , there are two models for how this charge is distributed : either in one row with average hole density 0.5 hole per copper site or in two rows with 0.25 hole per copper .
these two alternatives are often somewhat simplistically referred to as site order vs bond order .
the strongest evidence distinguishing between the alternatives comes from x - ray data on the charge order@xcite at @xmath16=1/8 : non - observation of diffraction harmonics suggests a sinusoidal distribution of charge . for a 4-cu repeat distance , two insulating and two charged rows would be exactly sinusoidal , whereas one charged and three insulating rows should have signifigant harmonic content .
a similar conclusion was reached by @xmath17sr lineshape analysis@xcite . however , the charge ordering peaks are weak , and it remains possible that fluctuations or disorder could wash out the harmonics .
direct evidence for the density on the charged stripes is found from low temperature nqr measurements@xcite , which find values @xmath18 . while this is close to the lower value
, the small difference can also be understood : this is a _ local _ measurement , and it is expected that some holes will be pushed off onto the magnetic stripes .
indirect evidence favoring the lower hole density includes the fact that it is easier to understand the properties of afm stripes in terms of even - leg ladders ( e.g. , the afm stripes would be two coppers wide at @xmath16=1/8)@xcite , and that the stripe phase appears to terminate when the average doping approaches @xmath16=0.25 ( section vii.a.1 ) .
the lower doping is also more consistent with the tj model simulations of white and scalapino@xcite .
tallon@xcite finds an optimal doping at @xmath19 for all cuprates , with respect to which the stripes terminate at a doping @xmath20 .
however , it is hard to reconcile a common optimal doping with muon spin resonance data@xcite , which find @xmath21 is optimized at very different values of @xmath22 ( @xmath23 is the pair density , which seems to scale with carrier concentration@xcite , and @xmath24 effective mass ) for lsco and ybco .
we assume instead that @xmath25 scales with @xmath26 , so if @xmath19 for lsco , it is 0.21 for ybco , in good agreement with several estimates@xcite .
this also resolves a problem with the thermopower . while the thermopower appears to be universal for most cuprates , and the best means of estimating the doping is from the room temperature thermopower , lsco is anomalous in that ` overdoped ' samples still have high thermopower@xcite .
if the doping for ybco is rescaled as above , however , the thermopower data of lsco fall on the universal curve .
hence , the anomaly for lsco is not in the thermopower , but in a too low value of @xmath21 , which is accompanied by a too low value of @xmath25 .
it is likely that these features are associated with a competing ltt phase , which is most prominent in lsco , and which also leads to the most nearly static stripe correlations .
taking @xmath25=0.21 for ybco gives @xmath27 , which we believe holds for _ all cuprates _ , including lsco , section vii.a.1 .
this would lead to very wide charged stripes near optimal doping : the width of the charged stripes @xmath28 satisfies @xmath29 , or @xmath30 cu wide .
hence , models of isolated quasi - one - dimensional charged stripes are likely to be valid only in the far underdoped regime , while for the good superconductors a better model would be a metal with intrinsic weak links@xcite .      for @xmath31
, the crossover @xmath32 can be identified with the 1/8 anomaly , where both charged and afm stripes have their minimal width ( 2 cu atoms ) .
there is considerable evidence that the doping 1/8 , in addition to its special stability , acts as a crossover in the properties of the stripes .
thus , uchida , et al.@xcite , studying the hall coefficient @xmath33 , find a crossover from one - dimensional behavior ( @xmath34 as @xmath35 ) for @xmath36 to two dimensional behavior ( coupled charged stripes ) for @xmath13 . in ybco , the spin gap grows slowly with doping for @xmath12 , then more rapidly for @xmath13 ; this behavior can be understood in terms of coupled spin ladders , as the coupling changes with the width of the charged stripes@xcite . in eu substituted lsco@xcite , the meissner fraction is negligibly small for @xmath12 , then grows roughly linearly with doping until @xmath37 , staying large up to at least @xmath38 .
finally , the two - magnon raman peak in lsco has a splitting at low temperatures which has been associated with stripes@xcite , on analogy with similar observations in la@xmath0sr@xmath1nio@xmath2 ( lsno)@xcite .
for @xmath12 the ratio of the two peak frequencies is constant and consistent with a simple stripe model ; for @xmath39 the lower frequency starts decreasing with doping .
moreover , the higher frequency loses intensity with doping ; near @xmath40 , the _ intensity _ of one mode approximately disappears , while the _ frequency _ of the other mode extrapolates to zero .
= 0.25 0.5 cm = 0.48 0.5 cm    the yamada@xcite plot , fig .
[ fig : y1]a , provides a severe constraint on any model of stripes : in lsco the incommensurability @xmath41 is found to grow linearly with doping for @xmath12 but to saturate for @xmath42 .
furthermore , the saturation value is just the incommensurability expected for @xmath43-doped stripes , @xmath44 .
a similar saturation has been reported in ybco , but different groups find different values for @xmath45 : @xmath46 [ ] or @xmath47 [ ] . if one of these values proves correct , it would suggest some nonuniversality in @xmath9 , perhaps associated with bilayer splitting . in the domain wall model ,
the saturation is interpreted as evidence that the stripes stop changing in width , and additional holes leak into the antiferromagnetic background , washing out the stripes .    for @xmath48
, the eps model agrees with the yamada plot , fig .
[ fig : y1]b : increase in doping causes the afm stripes to narrow , with no changes in the charged stripes .
note that , for concreteness we have assumed that the charged stripes have laf order ( section iii ) ; these stripes naturally act as antiphase boundaries ( apbs ) , consistent with the neutron evidence . however , the naive prediction of the eps model is in disagreement with the yamada plot for @xmath49 , fig .
[ fig : y1]c : as the charge stripes get wider , the incommensurability should decrease , while the neutron peaks may broaden if the wider charge stripes do not act as apbs . this behavior is not observed .
however , the model can be simply modified to explain the observed saturation , fig .
[ fig : y1]d , e .
this would be a commensurability effect , with part of the sample pinned at 1/8 doping while the rest forms a different phase ( e.g. , at 1/4 doping ) where no stripes are present .
such behavior is well known in nickelates , where coexistence of 1/3 and 1/2 stripes is common .
evidence for such commensurability effects can be found by comparing@xcite the chemical potential @xmath17 in lsco@xcite and lsno@xcite . in both materials , @xmath17 is approximately independent of doping between half filling ( @xmath50 ) and @xmath51 , with @xmath52 in lsco , 1/3 in lsno , which has long - range charge order .
remarkably , the simple prediction of a macroscopic phase separation model ( @xmath17 constant throughout the phase separation regime ) is violated in lsno : @xmath17 is _ not _ constant for @xmath53 , although the stripe phase persists over the full doping range .
this is presumably a commensurability effect , reflecting the coexistence of 1/3 and 1/2 stripe phases for intermediate dopings .
the same doping dependence of @xmath17 is found in lsco : constant for @xmath54 , variable for @xmath49 .
indeed , the lsno data can be scaled to that of the cuprates , giving a common doping dependence @xmath55 in both compounds@xcite .
recently local charged domains have been directly visualized by stm studies in bi2212@xcite with gaps ranging from slightly overdoped to strongly underdoped .
such behavior is extremely difficult to explain in the domain wall picture , but is easy to understand in terms of eps , and indeed is consistent with the commensurability effects discussed above .
the gap distribution is found to be broad rather than bimodal , but that is expected since eps is very sensitive to charged impurities , and the local gap will correlate with the local density of dopant , presumably interstitial oxygen in bi2212 .
this sensitivity to impurities leads to the question of which came first : is a preexisting eps sensitive to local impurities , or does oxygen clustering@xcite provide the driving mechanism for domain formation ?
since the electronic inhomogeneity seems characteristic of most cuprates while there is considerable variety in the doping counterions , the simpler interpretation would appear to be that the eps is primary .
thus , in la@xmath3cuo@xmath7 ( lco ) , the interstitial oxygens are highly mobile , allowing the domains to grow to macroscopic size .
similar clusters form in ybco ( here associated with chain oxygens ) , but can be suppressed in fully oxygenated samples by quenching@xcite . on the other hand ,
well formed stripes appear when the doping counterions are least mobile , in lsco .
the domains in bi2212 would seem to be consistent , but a number of questions remain about the role of annealing .
while these domains are regularly found in stm studies when the samples are cleaved at low - temperature , they can be annealed out at high temperatures@xcite .
this could be caused by a ` melting ' of the eps , as is found in lco near room temperature@xcite , and can be tested by careful annealing studies .
alternatively , it may be a question of pinning the eps .
certainly , it is known that a weak domain disorder is necessary to pin vortices , to observe the vortex lattice via stm.@xcite also , the role of the interstitial oxygen in the superlattice modulations in bi2212 needs to be better understood .
the strength of this modulation suggests that it is associated with an ordering of the interstitial oxygen , but this should lead to a strong correlation between the domains and the superlattice , which seems not to be the case .
further evidence for electronic inhomogeneity comes from microwave measurements : anomalies in bi2212 have been interpreted in terms of similar electronic domains@xcite , suggesting that they are representative of the bulk , while measurements on other cuprates@xcite find similar anomalous behavior which was interpreted in terms of pinned cdw s , possibly stripe related .
a domain picture would also explain the persistence of nodal quasiparticles in the underdoped regime , at least down to 1/8 doping@xcite .
table i summarizes the discussion of the previous subsections , and compares the predictions of eps models     model & @xmath56&apb?&saturation & domains?&crossover & termination + & & & at @xmath57?&&_at _
@xmath58?&at @xmath59 ? + domain wall&fixed pattern&@xmath60&@xmath60&@xmath61&@xmath62&@xmath62 + & holes spread out & & & & & + eps ( simple)&charge stripes grow&@xmath63&@xmath61&@xmath61&@xmath60&@xmath60 + & in width & & & & & + eps ( commensurate & domain phase:&@xmath63&@xmath60&@xmath60&@xmath60&@xmath60 + pinning)&@xmath58 : stripes & & & & & + & @xmath64 : uniform & & & & & +    with domain wall models of stripes .
domain wall models arose in early unrestricted hartree - fock ( uhf ) calculations@xcite for the doped tj model , as heterogeneous ground states in which the holes are confined to domain walls between afm domains which act as apbs .
these domain walls are not driven by phase separation@xcite ; the phase is realized as a long - period modulated afm , with holes doping the rows of spins where the moment changes sign@xcite . when @xmath65 , there is one hole per row , a large hole doping which is inconsistent both with experiments on stripes in the cuprates , and with density matrix renormalization group ( dmrg ) calculations@xcite .
more recent calculations reduce the doping , in ways that bring the model close to the phase separation picture : ( 1 ) assuming cdw order along the charged stripes@xcite lowers @xmath9 to 0.5 , while ( 2 ) letting @xmath66 results in @xmath67.@xcite since the cdw order ( case 1 ) lowers the free energy of the charged stripes , it should be possible to dope the system all the way to a pure cdw at @xmath68 . in case 2 ,
it is not explained why a low density appears ; a good possibility is that this doping corresponds to the van hove singularity ( vhs ) , as in the phase separated stripes .
( this could be checked by varying @xmath69 and calculating @xmath70 . )
if such low - density domain walls do exist , the striped phase must somehow terminate when the average doping approached 0.2 , but how this happens is not explained .
however , this calculation@xcite is not fully uhf , being restricted to periodic arrays , and uhf calculations@xcite find instead ferromagnetic charged stripes which can be understood within a phase separation model ( see below ) .
since the stability of these @xmath71 domain wall stripes is unclear , in table i a comparison is made with the earlier , higher - hole - density domain wall models .
for such high doping , stripe termination at @xmath72 , or a major crossover near 1/8 , are both difficult to interpret .
stripes as apbs are more natural in domain wall models , but we shall show that they can also arise in eps models . note that the distinction between the models may be becoming moot , since in the domain wall models there is now a search for a ` hidden ' order parameter@xcite , and it is common for secondary order parameters to be expressed predominantly on domain walls of the primary order@xcite
.    in conclusion , assumption of a charged stripe doping @xmath73 reconciles the neutron diffraction data , evidence for a termination of the stripe phase near @xmath9 , and the 1/8 anomaly as a crossover effect near @xmath74 , while    commensurability effects can explain the saturation in the yamada plot and the stm observation of charge domains .
further supportive evidence will be presented in section vii.a . in the remainder of this paper
, we will assume that @xmath75
.      in any phase separation model , a key issue is understanding the nature of the charged stripes . indeed , since superconductivity seems to arise predominantly on these stripes , such understanding is likely to play an important role in elucidating the origin of the high superconducting transition temperatures . for the stripe phase to exist
, the doping @xmath9 must be particularly _ stable_. this can arise via an electronic _ instability _ , which opens up a gap over much of the fermi surface , making the electronic phase nearly incompressible .
this ` stability from instability ' is a fairly general feature , underlying , e.g. , hume - rothery alloys@xcite .
[ this is a modification of an argument due to anderson@xcite . ] here , we explore a number of candidates for the predominant electronic instability .    in a related paper@xcite , we will provide strong evidence that this ` hidden order ' is a form of cdw , which could include the flux phase .
however , here we will explore a wider variety of possibilities .
one issue is that in the low doping limit the charge stripes act as apbs for the afm stripes .
such an effect arises naturally if the charged stripes have some residual magnetic interaction , and we will explore this possibility . however , in nickelates charged stripes coupled to a cdw are found to act as apbs@xcite .
the large hubbard on - site repulsion @xmath76 plays an important role .
strong correlation effects lead to two classes of charged stripe phases : either the constraint against double occupancy leads to magnetic order ( magnetic charged stripes ) or kinetic energy dominates , leading to a magnetically disordered phase ( paramagnetic charged stripes ) , with correlations leading to reduced hopping , as in tj@xcite and slave boson@xcite models .
section iii will provide examples of both classes , denoted as class b and class a stripes , respectively .
class a stripes could be simply correlated metals ( as in tj or slave boson calculations ) or could have additional , e.g. , cdw , order .
a crossover from magnetic to correlated paramagnetic groundstate arises as a function of doping in models of itinerant ferromagnets@xcite .    for completeness
, class c stripes are defined as those arising not from phase separation but from long - range modulated afm order ( domain wall stripes ) .
such stripes have been described in detail@xcite , and will not be further considered here .    in the next section ,
we show that class b ( magnetic ) charge stripes can arise in a mean - field hubbard model@xcite , with the charged stripes displaying either ferromagnetic ( fm ) or linear antiferromagnetic ( laf ) ( ordering vector @xmath77 ) order .
the laf stripes are very similar to white - scalapino stripes@xcite . the fm phase is stabilized by vhs nesting
. this fm phase may be present in ruthenates@xcite , but is unlikely to be relevant for the cuprates ( for one thing , the fm stripes are likely to be diagonal , and do not form apbs@xcite , contrary to experiment ) .
we have suggested that other vhs - stabilized phases are more likely@xcite ( see also ref .  ) , and here we explore the properties of a class a cdw phase@xcite . at a doping @xmath78 ,
the effects of correlations are relatively weak , renormalizing the bandwidth by a factor of @xmath632 .
thus for the present calculations on paramagnetic stripes , it will be assumed that renormalized parameters are used , and other effects of strong correlations will be neglected .    while laf stripes are most stable when @xmath65 , we explore the possibility that they can be stabilized even when @xmath71 by on - stripe cdw or superconducting order .
ordinarily , the cdw phase is believed to _ compete _ with strong correlation effects , but we find ( section v.b ) that an unusual form of cdw phase can _ coexist _ with laf order : the charge density varies between zero and one ( not two ) holes per atom .      in a stripe array ,
the alternating stripes are associated with the two stable thermodynamic phases . here ,
we summarize the different ways these stripes are denoted in this paper .
the stripe with lower hole doping is variously referred to as ` insulating ' or ` antiferromagnetic ' ( afm ) .
( these are also known as magnetic stripes , since the charged stripes have considerably weaker magnetic order , but we will avoid that notation here . ) the stripe with higher hole doping is generically referred to as ` charged ' or ` hole - doped ' . at low temperatures ,
these stripes are also ` superconducting ' , but at high temperatures , they are stabilized by some ` hidden order ' , and one purpose of this paper is to explore a number of possible orders .
the orders fall into two groups : ` magnetic charged ' ( class b ) stripes could have fm or laf order ( the latter are white - scalapino - like stripes ) , while ` paramagnetic ' ( class a ) stripes could have cdw or flux - phase order .
strong coupling models would seem to be natural for producing phase separated or striped ground states .
any magnetic ordering avoids double occupancy , while changing from one form of magnetic order to another , via , e.g. , doping , requires highly collective spin rotations , as competing orders are orthogonal .
while superexchange leads to antiferromagnetic ( afm ) insulators at half filling , doping tends to favor textures with parallel spins ( e.g. , ferromagnetic ( fm ) ) to maximize hole hopping .
such ferron phases were introduced long before high @xmath21@xcite , but it remains controversial whether such states are ground states of the hubbard model@xcite . while the tj model does have phase separation for large @xmath79 , it is unclear whether such phases extend to the values @xmath80 expected for the cuprates@xcite .
we have found phase separated solutions of the hubbard model at mean - field level@xcite . while these solutions are metastable in uhf calculations , they closely resemble the ws stripes , and provide an interesting example of phase - separation mediated stripe phases .
we find a well defined surface tension for wide , isolated stripes , which decreases and changes sign as the stripes become narrower .
when the surface tension becomes negative , the stripes no longer remain straight , and spontaneously meandering solutions are found .
these solutions are found by considering only low - order commensurate phases , with wave vector @xmath81 0 or @xmath82 only .
the bare dispersion is @xmath83 with @xmath84 .
the hubbard interaction @xmath85 leads to magnetic order with a mean - field magnetization @xmath86 at wave vector @xmath87 , and the quasiparticle dispersion becomes @xmath88 where @xmath89 and @xmath90 for the cuprates , we expect@xcite @xmath91 , @xmath92 and @xmath93 .
the magnetization is found self - consistently from @xmath94 with fermi function @xmath95 .
the resulting free energy is @xmath96 with @xmath97 @xmath98    = 0.40 0.5 cm    the competing phases include afm for @xmath99 , fm with @xmath100 , and a linear antiferromagnet ( laf ) with @xmath101 . when the laf stripes are two cells wide , this laf phase closely resembles the white - scalapino stripes , fig .  [ fig : y1]e .
the afm state has lowest free energy at half filling , but ( for @xmath65 ) is unstable for finite hole doping . for @xmath65 ,
there is phase separation between the afm and laf phases , while for finite @xmath69 the phase separation is between afm and fm phases , fig .
[ fig:10 ] .
( when electron - phonon coupling is included , it is found that the fm phase is unstable with respect to a cdw phase@xcite . )    = 0.40 0.5 cm    the laf stripes for @xmath65 are discussed in ref .  .
for @xmath71 electron - hole symmetry is absent ; for hole doping @xmath102 , there is phase separation to a fm phase , consistent with recent simulations by vozmediano , et al.@xcite , who find a uniform fm phase at @xmath103 for @xmath104 , @xmath105 .
however , _ on the electron - doped side a uniform afm phase is stable over a large doping range _ , suggestive of the asymmetry found in the cuprates .
the dot - dashed curve shows that in the phase separation regime , the low energy physics can be approximated by the form of free energy assumed in ref .
, @xmath106 ( neglecting a term linear in @xmath16 ) with @xmath9 the hole doping of the uniform charged phase . the fm phase is stabilized by vhs nesting , inset to fig .  [ fig:10 ] , as found previously@xcite .
the tangent construction tends to select the fm phase at dopings somewhat away from optimal nesting ( fermi energy above midgap ) .
both regimes of phase separation seem to be driven by hole delocalization : one - dimensional ( along the laf rows ) when @xmath65 , two - dimensional for finite @xmath107 .
the resulting phase diagrams @xmath16 vs @xmath76 , fig .
[ fig:6 ] , are strikingly different . for @xmath65 , fig .
[ fig:6]a , the phase separation is between the afm and paramagnetic ( laf ) phase for @xmath108 ( @xmath109 ) , while for finite @xmath69 , fig .
[ fig:6]b , there is generally a vhs - stabilized fm phase . for small @xmath76 and @xmath71
, there is a regime where simple spin - density wave theory works and a uniform afm phase is stable , but when @xmath65 phase separation persists for all finite @xmath76 .
note that @xmath110 marks a crossover between class a ( paramagnetic ) and class b ( magnetic ) charged stripes .
the value @xmath110 is close to the @xmath111 expected in the cuprates , although for finite @xmath69 @xmath110 decreases , @xmath112 for @xmath113 , and the range of @xmath76 for which paramagnetic stripes are stable becomes very small . while these stripes are metastable in uhf calculations , we will show below that the stability of phase separating stripes can be enhanced by _
additional interactions _ beyond the pure hubbard model .
in ordered stripe arrays it is found@xcite that the afm stripes have a mott - hubbard gap , and the features near the fermi level are associated with the charged stripes .
this charge stripe dispersion shows a series of quasi - one - dimensional features which qualitatively resemble the bands of an isolated stripe produced by qse . in this section
, we make a quantitative comparison with isolated stripes , and explore the mechanism of qse - induced van hove splitting .    for a single stripe @xmath28 cu atoms wide ,
the dispersion is still given by eq .
[ eq:1 ] , but the allowed @xmath114 values are quantized , with @xmath115 , @xmath116 .
these are in fact the usual quantized bloch bands , but for large @xmath28 the quantization is not noticeable . for small @xmath28 , the dispersion appears as a series of @xmath28 overlapping one - dimensional ( 1d ) dispersions .
equation  [ eq:1 ] can be rewritten as @xmath28 1d dispersions @xmath117 these are the qse , with corresponding density of states ( dos ) shown in fig .
[ fig:1 ] .
notice that for @xmath118 , the vhs is readily apparent in the dispersion .
even down to @xmath119 , the vhs is clearly defined ( albeit only within a finite interval ) as the _ locus of energies where all subbands overlap_. in fact , _ the qse opens a gap at the vhs _ , effectively lowering the kinetic energy of the electrons just like a conventional ( e.g. , cdw ) gap .
the vhs splitting can be found from eq .
[ eq:1e ] : @xmath120 this splitting has two consequences : first , the vhs splitting enhances the stabilization of the striped phase ; but second , the qse gap competes with other gaps , such as cdw s and superconductivity . however , while the qse splits the vhs peak , substantial dos remains ungapped , so additional instabilities remain possible . this competition will be discussed further in the next section .
= 0.40 0.5 cm    it should be noted that this is the first direct demonstration that the vhs can be defined on a stripe only two atoms wide .
the definition is quite analogous to the standard definition in two dimensions : the point at which the bands cross over from electron - like to hole - like .
figure  [ fig:1a ] compares the gaps of a single stripe with those found in the ordered stripe array@xcite ; the array is labelled @xmath121 when the magnetic stripes are @xmath24 coppers wide and the charged stripes are @xmath122 coppers wide . in the array calculation ,
no competing order was introduced on the charged stripes , so the qse provides the only gap .
figure  [ fig:1a ] shows that there is a very good match for both 2 cu wide and 6 cu wide stripes , although for the 2 cu stripe , the vhs gap is somewhat larger for the single stripe than in the array . from eq .
[ eq:1e ] , the 1d bands have dos peaks at band bottom and top ; the band bottom corresponds to @xmath123  i.e. , the dispersion from @xmath124 is flat , at the energy corresponding to the lower dos peak .
( the intensity along @xmath125 is given by a structure factor , which does not directly come into a single stripe calculation . ) along @xmath126 one should see the 1d dispersion extrapolating to the band top at @xmath127 .
given the good agreement , it should be possible to analyze competing orders on a single stripe , for which the calculations are simpler ( no need to self - consistently adjust doping on each row to account for charging effects ) .
= 0.40 0.5 cm
in this section , we will develop two closely related models of the competition of cdw order and superconductivity _ on a single charged stripe_. the first is a class a paramagnetic stripe , stabilized by electron - phonon coupling@xcite , while the second has dominant electron - electron coupling , with ( class b ) magnetic charged stripe order .
the earlier calculation@xcite used a van hove stabilized cdw model to describe the doping dependence of the pseudogap . here
we reapply the model for a single stripe , by introducing the following modifications .
( 1 ) the calculation is carried out on a single , finite width stripe .
( 2 ) pinning to the vhs arises naturally , since all stripes are at the same doping .
( 3 ) for closer approximation to experiment , d - wave superconductivity is assumed .
( 4 ) correlation effects due to hubbard @xmath76 are neglected : previous slave boson calculations suggest that the main effect is a bandwidth renormalization of a factor @xmath632.@xcite    we briefly recall the energy dispersion and the gap equations of the model@xcite , generalized to d - wave . in terms of a function
@xmath128 the gap functions are @xmath129 for superconductivity , and @xmath130 for the cdw .
the energy eigenvalues are @xmath131 and their negatives , with @xmath132 @xmath133 , @xmath134 , @xmath135 , @xmath136 , and the nesting vector @xmath137 .
the gap equations are @xmath138\over 2 e_{+,k}}\tanh{{e_{+,k}\over 2k_bt } } \nonumber \\ -{e_{-,k}^2-\epsilon _ { \vec k+\vec q}^2-\theta_{\vec k+\vec q}[\delta_{\vec k}^2+g_{\vec k}^2]\over 2 e_{-,k}}\tanh{{e_{-,k}\over 2k_bt } } \bigr ) , \label{eq:23}\end{aligned}\ ] ] @xmath139 with interaction energies @xmath140 and @xmath141 , and @xmath142 , @xmath143 .
the cdw - superconducting competition was studied in bulk in ref .  .
the previous results are recovered for a sufficiently wide stripe ( @xmath63100 cu wide ) . for narrower stripes , it is found that the quantum confinement gap severely interferes with alternative gap formation .
figure  [ fig:19 ] illustrates how the various gaps vary with stripe width , near @xmath75 .
the data are plotted vs average doping , assuming a regular stripe array with 2-cu wide afm stripes ( which do not contribute to the gaps ) and n - cu wide charged stripes of doping @xmath9 , giving an average hole doping @xmath144 .
the cdw gap is most sensitive to stripe width , but in the narrowest stripes superconductivity is also suppressed ( the suppression is stronger for a d - wave gap ) .
strong instabilities are possible on the stripes , but they are shifted in doping away from @xmath75 .
thus , the cdw instability requires both @xmath145 and @xmath146 to be near the fermi level ; for a two - cell wide stripe , this is only possible near @xmath50 . on the other hand ,
superconductivity is possible anywhere , if the coupling is strong enough .
for a two - cell wide stripe , the optimal superconductivity arises when the fermi level is at the _ one - dimensional vhs _ at the edge of one of the stripe subbands .
this depends on @xmath69 , and for @xmath113 falls at @xmath147 ( a larger gap is found on the electron - doped side , @xmath148 ) .
it should be noted that , even though the superconductivity is assumed to be d - wave , in general a finite minimum gap is found on the stripe , even when the cdw gap is zero .
this is because the vanishing d - wave gap can be sampled only when the point @xmath149 is sufficiently close to the fermi surface , which in general requires @xmath28 to be odd ( recall that the allowed values of @xmath114 are integer multiples of @xmath150 ) or very large .
= 0.40 0.5 cm      in the above calculations , the @xmath151 s arose from electron - phonon coupling@xcite .
similar contributions follow from electron - electron coupling , in an extended hubbard model .
for instance , the near - neighbor coulomb repulsion has the following mean - field expansion :    @xmath152    where @xmath153 and @xmath154 .
the first two terms in eq .
[ eq:24b ] renormalize the chemical potential and the hopping @xmath155 respectively , and can be neglected .
the terms in @xmath156 comprise a pseudospin triplet of cdw - like distortions , with @xmath157 representing a cdw similar to the one discussed above , @xmath158 being related to the low - temperature tetragonal distortion , and @xmath159 an orbital antiferromagnet , closely related to the flux phase .
the remaining terms are a d - wave superconductor .
the coefficients of the terms must be found self - consistently by solving the gap equations :    @xmath160    @xmath161    @xmath162    @xmath163    the terms @xmath164 and @xmath165 have recently been discussed by valenzuela and vozmediano@xcite . a detailed discussion of the competition between the three cdw - like modes is given in ref . .
retaining only the superconducting term in eq .
[ eq:24b ] , the interaction can be derived from a quartic term @xmath166 with @xmath167 . assuming @xmath168 , the gap equations can be written in the form @xmath169 ( @xmath170 ) , with @xmath171 and @xmath172 .    for the uniform charged state ( infinitely wide stripe ) @xmath173 , and
the gap symmetry can be simply analyzed .
the symmetry can be either d - wave ( @xmath174 ) or extended s - wave ( @xmath175 ) , with the choice @xmath176 giving the largest gap .
the @xmath177 term is always bcs - like , having the opposite sign from @xmath178 , while the @xmath179 term has the opposite sign from @xmath177 , the integral being dominated by the regions near the vhs s .
hence , there are two possibilities : ( i ) attractive ( @xmath180 ) d - wave superconductivity or ( ii ) repulsive ( @xmath181 ) extended s - wave .
however , the latter would require @xmath182 , which does not arise in the present model , so only case ( i ) is possible .
these considerations readily generalize to a finite - width stripe , for which @xmath183 .
agterberg , et al.@xcite recently introduced a model for ` exotic ' superconductivity in multiband superconductors .
if the fermi surface consists of several inequivalent but degenerate pockets , the order parameter can consist of symmetry allowed superpositions of the order parameters of the individual pockets .
equation  [ eq:37 ] can be thought of as a form of exotic superconductivity , with the degenerate vhs s playing the role of hole pockets .
= 0.40 0.5 cm    in the above calculations , it was implicitly assumed that the doping is high enough that the only role of the on - site repulsion @xmath76 is to renormalize the band parameters .
however baskaran@xcite recently estimated that near optimal doping correlation effects remain stronger than the kinetic energy associated with hopping .
hence , it is important to look for stripe ground states which minimize this on - site repulsion ( class b stripes ) .
the linear antiferromagnet ( laf ) stripes discussed in section iii are a good candidate for the cuprate charged stripes : they closely resemble the white - scalapino stripes@xcite , have an appropriate doping , close to @xmath75 , include strong correlations , and have the additional advantage that a two - cell wide laf charge stripe acts as a natural apb for the af stripes . in this section
, we will explore these stripes , and show that they can be further stabilized by additional interactions .    a special form of strongly correlated cdw is found to exist on a laf .
the charge and spin distribution is shown in the insert to fig .
[ fig:20b ] , with the corresponding dispersion in figure  [ fig:20a]c .
there is a strong antiferromagnetic ordering on one sublattice , while most of the holes are confined on the other , nonmagnetic sublattice .
whereas in a conventional cdw the charge density is zero on one sublattice and two on the other , in this strong coupling case the hole density varies from @xmath184 to 1 , and there is no double occupancy , fig .
[ fig:20b ] . whereas the paramagnetic stripes were extremely sensitive to quantum confinement ,
these magnetic charged stripes are much less so : this cdw is stable almost independently of the stripe width . from fig .
[ fig:20a]c , it can be seen that the gapped fermi surface still has hole pockets near @xmath149 , which would lead to conducting stripes , consistent with optical properties@xcite .
however , it is only found near a hole doping @xmath68 , and so does not appear to be relevant for stripe physics in the cuprates .
= 0.40 0.5 cm    away from this doping , cdw instabilities are relatively weak and it is possible to stabilize d - wave superconductivity , figs .  [ fig:20a]b,[fig:20 ] .
while the overall dispersion varies with stripe width , the superconducting gap is also relatively insensitive to the width , and actually increases for the narrowest stripes , fig .
[ fig:20c ] .
note that the order parameter is not a pure d - wave , the gap along the stripe being larger .
such a large anisotropy is not consistent with tunneling measurements of the gap ; it is possible that the anisotropy is reduced by strong interstripe coupling . on the other hand ,
a large gap anisotropy has been found in ybco@xcite , where the stripes are aligned along the chain direction@xcite .
= 0.40 0.5 cm    = 0.40 0.5 cm    note in figs .
[ fig:20a]b ,  [ fig:20 ] that the combination of laf and d - wave order leads to a finite minimum gap over the full fermi surface .
while the pure laf phase is not stabilized by the vhs , the d - wave superconductivity is optimized when the fermi level is at the @xmath77 vhs  at essentially the _ same doping _ , @xmath185 , as the vhs on a paramagnetic stripe !    = 0.40 0.5 cm
= 0.40 0.5 cm    in ref .
, the fermi surface was calculated for a series of ordered stripe arrays .
these results can now be compared to experimental photoemission data@xcite . for this purpose
we replot the data as integrated spectral weight over a finite energy cut within energy @xmath186 of the fermi surface .
figure  [ fig:21]a shows a cut with @xmath186 = 200mev , for the model of a 1/8 doped stripe array ( i.e. , @xmath187)@xcite . the pattern is readily understood : the stripe superlattice leads to a number of quasi - one - dimensional bands ; however , due to structure factor effects they have significant intensity only near the original fermi surface , solid line in fig .
[ fig:21]a . for comparison with experiment ,
the calculated spectral weight is symmetrized @xmath188 in fig .
[ fig:21]b to represent a sample with regions of stripes running along both @xmath189 and @xmath190 directions .
finally , an empirical matrix element is included , fig .
[ fig:21]c , which extinguishes spectral weight along the zone diagonal , @xmath191 , similar to the matrix element assumed in analyzing bi2212 , refs .  .
the resulting fermi surface maps for several values of @xmath186 are illustrated in ref .
for paramagnetic charged stripes , and in figs .
[ fig:22 ] ( @xmath43 ) and [ fig:23 ] ( @xmath192 ) for laf stripes .
= 0.40 0.5 cm    for both models , the stripe band nearest @xmath77 is in good agreement with experiment : there is little dispersion perpendicular to the stripe , while the intensity falls off toward @xmath193 due to the structure factor effect . in general , the laf stripes are in better agreement with experiment , since the additional subbands predicted for paramagnetic stripes ( moving from @xmath77 toward @xmath194 ) are not seen in the experiment .
while the matrix element improves the agreement , theory suggests that this effect is present only for certain photon polarizations@xcite .
one disagreement with experiment for both models is that for shallow energy cuts ( 30 , 100 mev ) the experiment still finds a smeared dispersion rather than a sharp fermi surface .
this is presumably an effect of stripe fluctuations .
it should be noted that all the spectral weight in figs .
[ fig:21]-[fig:23 ] is associated with the charged stripes ; the lower hubbard band of the afm stripes lies below 0.5ev in lsco .
it is somewhat surprising that the spectral weight nearest the fermi level is near @xmath77 , since this is where the pseudogap arises .
nevertheless , our calculation reproduces both the ( quantum confinement ) pseudogap , fig .
[ fig:1a ] , and the spectral weight distribution .
the two most common interpretations for the pseudogap are in terms of either precursor pairing or a competing order parameter , which may be magnetic@xcite or charge - density wave ( cdw)@xcite .
we will argue in the next subsection that the latter possibility is more likely .
such competing instabilities arise naturally in a phase separation model ( e.g. , refs . ) , and we have long argued that the pseudogap is a manifestation of local phase separation@xcite . it should be noted that a pseudogap arises in the nickelates@xcite in conjunction with stripe fluctuations , and turns into a true gap at the charge ordering temperature .
therefore , the fact that the pseudogap closes in the overdoped regime in the cuprates strongly suggests that the stripe phase terminates at the same doping .
direct evidence for this has recently appeared@xcite .
figure [ fig:3]a compares the intensity @xmath195 of the inelastic neutron peaks near @xmath196 as a function of doping in ybco@xcite and lsco@xcite . in a stripe picture
, @xmath195 should be a measure of the fraction of material in afm stripes . remarkably , the intensity extrapolates to zero at nearly the same doping in both materials , even though @xmath197 peaks at substantially different dopings . for both materials ,
this is the doping at which the pseudogap closes .
while in ybco , these peaks have been interpreted as a fermi surface nesting effect@xcite , this can not explain the charge order fluctuations , and in lsco the tranquada data@xcite is associated with _
elastic _ magnetic peaks .
furthermore , a study of the neutron diffraction pair distribution function@xcite for lsco finds evidence for charge fluctuations , presumably associated with stripes .
the excess fluctuations are maximal near @xmath103 , and terminate near @xmath198 .
strong reductions of thermal conductivity associated with stripe scattering also terminate at a comparable doping@xcite , while a recent optical study@xcite finds evidence for a quantum critical point at a similar doping , @xmath199 .
moreover , in bi2212 , tokunaga , et al.@xcite have introduced a new crossover temperature @xmath200 based on cu nmr , below which afm correlations develop ; they find @xmath201 near @xmath40 .
a recent nqr study of the slowing of spin fluctuations in re substituted lsco@xcite finds that the effective spin stiffness @xmath202 ( or equivalently the effective exchange constant ) scales to zero at a comparable doping ; the inverted triangles in fig .
[ fig:3 ] show @xmath203 .
as might have been anticipated from sections ii , iv , the doping dependence of @xmath204 changes radically below @xmath205 . note that while the integrated neutron intensity scales approximately with the area fraction of charged stripes , @xmath206 scales to @xmath63460k as @xmath207 .
this is only 1/4 of the actual spin stiffness , @xmath208 in the undoped afm .
the change by nearly a factor of four is suggestive of a dimensional reduction ( lower coordination ) , but for an isolated straight spin ladder , a factor of two might have been expected .
= 0.33 0.5 cm    it is important to note the proximity of this termination of the phase separation regime to the vhs : the arrow in fig .
[ fig:3]a shows the doping at which the pseudogap in the heat capacity@xcite closes , leaving an approximately logarithmic peak@xcite , while the @xmath61 indicates the point at which photoemission finds the vhs crossing the fermi level@xcite .
termination of the stripe phase close to a vhs is an important prediction of our eps model of stripes .
figure  [ fig:3 ] has a bearing on the discussion given earlier of optimal doping in ybco , or equivalently of the doping @xmath16 at which the stripe phase terminates .
setting optimal doping at @xmath209 for all cuprates puts the termination at doping @xmath210 .
however , this value is clearly too low for lsco : when the stripes are gone , magnetic correlations should be weak , but there is evidence for long range magnetic order at @xmath211@xcite and @xmath212 , the latter coupled with a suppression of @xmath21@xcite .
on the other hand , both facts are compatible with @xmath213 , fig .
[ fig:3 ] .
the termination of the pseudogap has also been interpreted in terms of a quantum critical point ( qcp)@xcite . at first glance , it looks straightforward to distinguish a model of eps , which is by nature first order , from a qcp , which is typically second order . however , the phases are actually quite similar .
due to strong coulomb effects , the eps is restricted to nanoscopic scale , so that , e.g. , the superconducting critical temperature @xmath21 changes smoothly with doping . on the other hand , near a qcp with quenched disorder , long - time correlations can lead to effective spatial inhomogeneity@xcite .
we suggested earlier@xcite that the peak and hump features seen in photoemission from bi2212 were associated with the charged and the afm stripes respectively .
as such , the intensity of the peak should have the doping dependence predicted for charged stripes , with the intensity increasing from zero at half filling , approximately linearly with doping @xmath16 .
this has now been verified experimentally@xcite . moreover
, the maximum intensity of the spectral weight occurs at the same doping@xcite @xmath9 discussed above , where the stripe phase terminates .
remarkably , the peak spectral weight closely tracks @xmath21 , suggesting that _ the superconducting pairs live on the charged stripes _ , as predicted by several models@xcite .
consistent with this , a number of measures of the strength of superconductivity ( condensation energy , critical current ) are optimized at this same point@xcite where the charge stripe intensity is maximum and afm stripes vanish .
the fact that @xmath21 itself is actually optimized at a slightly lower doping may be a hint that stripes can enhance the superconducting gap , as found above , fig .
[ fig:20c ] .
the sharp falloff of spectral weight at higher doping may be an indication for a _ second regime of phase separation _ in overdoped samples , where the new phase is not intrinsically superconducting .
some additional evidence for such a phase separation is discussed below .
it is important to note that as soon as the pseudogap closes , experiments find a peak in the density of states consistent with that expected for a vhs .
this is discussed in ref . on p. 1203
21 ) ; more recent evidence is in ref . .
for even higher dopings there are hints of a new regime of phase separation .
this is an improtant prediction of the model : if stripes are stabilized by lowering the electronic energy at a vhs , then there could be a _ second range of phase separation _ in the overdoped regime , with the non - superconducting phase totally unrelated to the afm state at half filling .
earlier experimental evidence for this second regime of phase separation in the cuprates has been presented in references ( section 11.6 ) and ( see also wen , et al@xcite ) .
recently , loram@xcite has presented a detailed analysis of the heat capacity of la@xmath0sr@xmath1cuo@xmath2 over an extended doping range . in the range @xmath214 ,
the temperature dependence of the somerfeld constant @xmath215 is consistent with that expected for a uniform phase close to a vhs ( ref . ,
21 ) . for lower doping
, @xmath215 falls off too rapidly to be associated with rigid - band filling away from a vhs , but can be easily explained in terms of the opening of a pseudogap . for higher hole doping , @xmath215 again falls off@xcite too rapidly to be due to rigid band filling , although without an obvious pseudogap forming .
the near symmetry of @xmath215 between the overdoped and underdoped sides is suggestive of two regimes of phase separation .
further possible evidence for this second regime comes from tunneling studies@xcite which find a subdominant order parameter in optimally and over - doped ybco ( i.e. , in the same doping range where a split superconducting transition is found@xcite ) .
this could arise from an allender - bray - bardeen mechanism@xcite , with the d - wave superconducting stripes inducing s - wave ( or @xmath216 ) superconductivity on adjacent normal - fermi - metal stripes .
also , microwave measurements@xcite find a crossover near optimal doping from superconductor - insulator to superconductor - metal domains .
uemura@xcite has also recently proposed phase separation in the overdoped cuprates . in extending the uemura plot to the overdoped regime
, a ` boomerang ' effect is found : both @xmath21 and @xmath217 decrease in parallel@xcite .
here @xmath218 is the superfluid density and @xmath24 is the effective mass , and their ratio is extracted from @xmath17sr measurements of the penetration depth .
uemura@xcite suggests that this proportionality between @xmath21 and @xmath217 can most easily be accounted for in terms of ( nanoscale ) phase separation between the superconductor and a non - superconducting phase .
thus for both overdoped and underdoped cuprates@xcite , the scaling between @xmath21 and @xmath23 can be understood in terms of phase separation .
this has important implications for superconductivity : _ superconductivity arises only at a unique composition near optimal doping . _ the present results suggest that even at this doping , the superconductivity must compete with a second order parameter .
above , we have presented evidence that the optimum @xmath21 occurs at a different doping in lsco than for the other cuprates , while the stripes terminate at approximately the same doping for all cuprates .
further evidence for this can be found by comparing the superconducting fluctuations@xcite , fig .
[ fig:38a ] . in all three materials , lsco , ybco , and bi2212
, some superconducting fluctuations appear to persist in the low doping regime , to temperatures in excess of optimal @xmath21 ( although considerably less than the pseudogap temperature ) .
the fluctuating superconductivity in lsco is very similar to that found in the other cuprates . on the other hand ,
if one assumed a universal curve of superconductivity , with @xmath219 vs @xmath16 the same for all cuprates , then the scaled superconducting fluctuations in lsco would be at anomalously high temperatures .
it seems simpler to assume that in lsco long range superconducting order is anomalously suppressed ( as also suggested by baskaran@xcite ) , in which case there is no reason for optimal @xmath21 to fall at the same doping .    in bi2212 and ybco ,
the onset of superconducting fluctuations falls close to the ` strong pseudogap ' regime@xcite , and hence at temperatures substantially below the weak pseudogap temperature@xcite , making it unlikely that the weak pseudogap is associated with preformed pairs . in lsco
, there is some confusion , since the pseudogap temperatures reported by batlogg , et al.@xcite are considerably higher than those found by ido , et al.@xcite .
the ido @xmath220 s fall close to the fluctuation onset temperatures in fig .
[ fig:38a ] .
= 0.33 0.5 cm    the fluctuations observed in fig .
[ fig:38a ] are very suggestive of the theoretical results of fig .
[ fig:20c ] : thus , quantum size effects enhance the average gap on a narrow stripe , which will lead to an enhanced mean field transition temperature in the underdoped regime . however , the coupling between stripes also weakens , leading to a loss of interstripe coherence , so the pairing is only evident as enhanced fluctuations.@xcite
recent experiments have provided considerable evidence for the presence of stripes and eps in the cuprates , but there remain many questions of how universal these are , how they arise and vary with doping , and how they interact with superconductivity .
we have here elaborated our earlier@xcite model of stripes driven by frustrated phase separation , in particular adducing evidence that the doping on the charged stripes is close to @xmath221 , and that when the average doping approaches this value eps terminates . moreover , near @xmath32 there is a crossover in stripe properties : for @xmath12 the charged stripes are quantum confined , for @xmath13 the afm stripes are so confined .
this model can explain the 1/8 anomaly ( @xmath222 ) , the anomalous hall effect ( @xmath34)@xcite for @xmath54 ( charged stripes confined , hence one - dimensional ) , and the growing spin gap in ybco for @xmath49@xcite .    on the important issue of the _ structure _ of a charged stripe , we have explored a number of possibilities without coming to any final conclusions . while there is evidence that superconductivity lives on the charged stripes , there also appears to be a second instability on these stripes , which stabilizes the stripe phase while competing with superconductivity .
we have shown that a semiquantitative understanding can be achieved by looking at the properties of a single doped ladder , and we have discussed how a number of instabilities ( both cdw and cuperconducting ) vary with ladder width .
we showed that strong correlation effects could lead to charged stripes with a residual magnetic order , introduced a simple model for white - scalapino - like stripes , and found novel superconducting and cdw instabilities associated with such stripes .
we illustrated how stripe order would affect arpes spectra , both dispersions and fermi surface maps .
future studies will apply the model to describing other properties of the cuprates .
certain anomalous features of strong - correlation calculations may find an explanation in underlying phase separation .
thus , the vanishing of the renormalized hopping parameter @xmath155 near half filling in slave boson calculations may reflect the vanishing of the charged stripes at half filling@xcite , while the frequently - observed pinning of the vhs near the fermi level@xcite is consistent with vhs - stabilized charged stripes .
finally , the idea of a commensurability effect near 1/8 doping , leading to a coexistence of _ domains _ for @xmath49 , provides a simple explanation for a large variety of experimental findings , including the saturation of the yamada plot , the direct observation of domains in stm studies , and a variety of microwave anomalies .
this may also lead to a resolution of the combined puzzle of magnetic neutron scattering incommensurability and the neutron resonance peak .
a stripe model provides a natural explanation of the incommensurability for @xmath48 , including a stripe reorientation transition at the metal - superconductor transition near @xmath223 .
however , a band picture@xcite ( with @xmath224 close to a vhs ) provides a superior model for the combined , frequency - dependent incommensurability _ cum _ resonance peak found near optimal doping in ybco and bi2212 .
an eps crossover from stripes to domains near @xmath225 would provide a natural explanation of these phenomena .
interestingly , the only calculation@xcite to approximately describe the doping dependence is based on a slave boson model , which as noted above tends to mimic eps behavior as @xmath207 .
* acknowledgment : * these computations were carried out using the facilities of the advanced scientific computation center at northesatern university ( nu - ascc ) .
their support is gratefully acknowledged .
we had stimulating conversations with s. pan and j.c .
ck s research was supported in part by nsf grant nsf-9711910 .
publication 776 of the barnett institute .
k. yamada , c.h .
lee , k. kurahashi , j. wada , s. wakimoto , s. ueki , h. kimura , y. endoh , s. hosoya , g. shirane , r.j .
birgeneau , m. greven , m.a .
kastner , and y.j .
kim , phys .
b*57 * , 6165 ( 1998 ) .
tranquada , b.j .
sternlieb , j.d .
axe , y. nakamura , and s. uchida , nature * 375 * , 561 ( 1995 ) ; j.m .
tranquada , j.d .
axe , n. ichikawa , a.r .
moodenbaugh , y. nakamura , and s. uchida , phys . rev .
lett * 78 * , 338 ( 1997 ) .
mook , presented at the sixth international conference on the materials and mechanisms of superconductivity  high - temperature superconductivity ( m@xmath226s - htsc - vi ) , houston , tx , feb .
20 - 25 , 2000 ; h.a .
mook and f. doan , nature * 401 * , 145 ( 1999 ) , and cond - mat/0103037 .
t. cren , d. roditchev , w. sacks , j. klein , j .- b .
moussy , c. deville - cavellin , and m. lagus , phys .
* 84 * , 147 ( 2000 ) ; c. howald , p. fournier , a. kapitulnik , cond - mat/0101251 ; s.h . pan , j.p .
oneal , r.l .
badzey , c. chamon , h. ding , j.r .
engelbrecht , z. wang , h. eisaki , s. uchida , a.k .
gupta , k .- w .
hudson , k.m .
lang , and j.c .
davis , unpublished ; and j.c .
davis , presented at the sixth internat .
conf . on spectroscopies of novel superconductors ( sns01 ) ,
chicago , may 13 - 17 ( 2001 ) , unpublished .
j. orenstein , presented at sns01 , unpublished . j.l .
tallon , cond - mat/9911422 , j.l .
tallon , j.w .
loram , and g.v.m .
williams , cond - mat/9911423 , j.l .
tallon , j.w .
loram , g.v.m .
williams , j.r .
cooper , i.r .
fisher , j.d .
johnson , m.p .
staines , and c. bernhard , phys .
b*215 * , 531 ( 1999 ) .
salkola , v.j .
emery , and s.a .
kivelson , phys .
* 77 * , 155 ( 1996 ) .
s. caprara , c. di castro , m. grilli , a. perali , and m. sulpizi , physica c*317 - 8 * , 230 ( 1999 ) ; g. seibold , f. becca , f. bucci , c. castellani , c. di castro , and m. grilli , cond - mat/9906108 .
markiewicz , phys .
rev . b*62 * , 1252 ( 2000 ) .
zacher , r. eder , e. arrigoni , and w. hanke , phys . rev . lett . * 85 * , 2585 ( 2000 ) .
j. eroles , g. ortiz , a.v .
balatsky , and a.r .
bishop , cond - mat/0008341 .
r.s . markiewicz and c. kusko , int . j. mod
b*14 * , 3561 ( 2000 ) .
white and d.j .
scalapino , phys .
lett . * 80 * , 1272 ( 1998 ) , _ ibid . _ * 81 * , 3227 ( 1998 ) .
zimmermann , a. vigliante , t. niemller , n. ichikawa , t. frello , j. madsen , p. wochner , s. uchida , n.h .
andersen , j.m .
tranquada , d. gibbs , and j.r .
schneider , europhys .
lett . * 41 * , 629 ( 1998 ) ; t. niemoeller , n. ichikawa , t. frello , h. huennefeld , n.h .
andersen , s. uchida , j.r .
schneider , and j.m .
tranquada , cond - mat/9904383 ; t. niemller , h. hnnefeld , t. frello , n.h .
andersen , n. ichikawa , s. uchida , and j.r .
schneider , j. low - t . phys . * 117 * , 455 ( 1999 ) .
kojima , h. eisaki , s. uchida , y. fudamoto , i.m .
gat , a. kinkhabwala , m.i .
larkin , g.m .
luke , and y.j .
uemura , physica b*289 - 90 * , 343 ( 2000 ) .
teitelbaum , b. bchner , and h. de gronkel , phys .
* 84 * , 2949 ( 2000 ) .
j. tworzydlo , c.n.a .
van duin , and j. zaanen , proc . of 2nd intern .
conf . on stripes and high tc superconductivity , rome 98 , to be published , j. supercond .
( cond - mat/9808034 ) .
tallon , physica c*168 * , 85 ( 1990 ) .
uemura , g.m .
luke , b.j .
sternlieb , j.h .
brewer , j.f .
carolan , w.n .
hardy , r. kadano , j.r .
kempton , r.f .
kiefl , s.r .
kreitzman , p. mulhern , t.m .
riseman , d.ll .
williams , b.x .
yang , s. uchida , h. takagi , j. gopalakrishnan , a.w .
sleight , m.a .
subramanian , c.l .
chien , m.z .
cieplak , g. xiao , v.y .
lee , b.w .
statt , c.e .
stronach , w.j .
kossler , and x.h .
yu , phys .
lett . * 62 * , 2317 ( 1989 ) .
optical evidence that @xmath16 and the superelectron density @xmath23 both scale with @xmath21 is found in h.l .
liu , m.a .
quijada , a. zibold , y .- d .
yoon , d.b .
tanner , g. cao , j.e .
crow , h. berger , g. margaritondo , l. forr , beom - hoan o , j.t .
markert , r.j .
kelly , and m. onellion , j. phys .
11 * , 239 ( 1999 ) .
consistent with @xmath25 = 0.21 [ y. tokura , j.b .
torrance , t.c .
huang , and a.i .
nazzal , phys . rev . * b38 * , 7156 ( 1988 ) ] , 0.24 [ u. tutsch , p. schweiss , r. hauff , b. obst , th . wolf , and h. whl , j. low - t . phys . * 117 * , 951 ( 1999 ) ] , 0.2 [ m. merz , n. ncker , p. schweiss , s. schuppler , c.t .
chen , v. chakarian , j. freeland , y.u .
idzerta , m. klser , g. mller - vogt , and th . wolf , phys .
lett . * 80 * , 5192 ( 1998 ) ] , 0.21 [ s.l .
drechsler , j. malek , and h. eschrig , phys .
rev . b*55 * , 606 ( 1997 ) ] .
a recent critique of tallon s bond valence sum calculation@xcite also finds evidence for a larger @xmath25 , but concludes that the technique can not accurately determine the mobile hole fraction on the cuo@xmath3 planes [ o. chmaissem , y. eckstein , and c.g .
kuper , phys .
rev . b*63 * , 174510 ( 2001 ) .
tallon , j.r .
cooper , p.s.i.p.n .
de silva , g.v.m .
williams , and j.w .
loram , phys .
* 75 * , 4114 ( 1995 ) .
jones , d.c .
christen , j.r .
thompson , r. feenstra , s. zhu , d.h .
lowndes , j.m .
phillips , m.p .
siegal , and j.d .
budai , phys .
rev . b*47 * , 8986 ( 1993 ) .
t. noda , h. eisaki , and s. uchida , science * 286 * , 265 ( 1999 ) .
v. kataev , b. rameev , a. validov , b. bchner , m. hcker , and r. borowski , phys .
rev . b*58 * , 11876 ( 1998 ) .
s. sugai and n. hayamizu , j. phys .
solids * 62 * , 177 ( 2001 ) .
g. blumberg , m.v .
klein , and s .- w .
cheong , phys .
lett . * 80 * , 564 ( 1998 ) ; k. yamamoto , t. katsufuji , t. tanabe , and y. tokura , phys .
lett . * 80 * , 1493 ( 1998 ) ; and s. sugai , n. kitamori , s. hosoya , and k. yamada , j. phys .
67 * , 2992 ( 1998 ) .
m. fujita , k. yamada , h. hiraka , p.m. gehring , s.h .
lee , s. wakimoto , and g. shirane , cond - mat/0101320 .
m. arai , t. nishijima , y. endoh , a.w .
garrett , s. tajima , k. tomimoto , y. shiohara , c.d .
frost , and s.m .
bennington , cond - mat/9912233 ; m. arai , y. endoh , s. tajima , and s.m .
bennington , int .
b*14 * , 3312 ( 2000 ) .
p. dai , h.a .
mook , r.d .
hunt , and f. dogan , phys .
rev . b*63 * , 54525 ( 2001 ) .
a. ino , t. mizokawa , a. fujimori , k. tamasaku , h. eisaki , s. uchida , t. kimura , t. sasagawa , and k. kishio , phys .
lett . * 79 * , 2101 ( 1997 ) .
m. satake , k. kobayashi , t. mizokawa , a. fujimori , t. tanabe , t. katsufuji , and y. tokura , phys .
* , 15515 ( 2000 ) .
balatsky , presented at sns01 , unpublished .
a. erb , a.a .
manuel , m. dhalle , f. marti , j .- y .
genoud , b. revaz , a. junod , d. vasumathi , s. ishibashi , a. shukla , e. walker , o. fischer , r. fluekiger , r. pozzi , m. mali , and d. brinkmann , sol .
* 112 * , 245 ( 1999 ) . o. fischer , presented at sns01 , unpublished .
hammel , a.p .
reyes , s - w .
cheong , z. fisk , and j.e .
schirber , phys .
lett . * 71 * , 440 ( 1993 ) . c. kusko , z. zhai , n. hakim , r.s .
markiewicz , s. sridhar , d. colson , viallet - guillet , a. forget , yu.a .
nefyodov , m.r .
trunin , n.n .
kolesnikov , a. maignan , a. daignere , and a. erb , unpublished ; r.s .
markiewicz , unpublished . x.j .
zhou , t. yoshida , s.a .
kellar , p.v .
bogdanov , e.d .
lu , a. lanzara , m. nakamura , t. noda , t. takeshita , h. eisaki , s. uchida , a. fujimori , z. hussain , and z .- x .
shen , phys .
* 86 * , 5578 ( 2001 ) .
j. zaanen and a.m. oles , ann .
physik * 5 * , 224 ( 1996 ) ; c. nayak and f. wilczek , phys . rev . lett . * 78 * , 2465 ( 1997 ) ; v.j .
emery , s.a .
kivelson , and o. zachar , phys .
rev . b*56 * , 6120 ( 1997 ) ; s.a .
kivelson , e. fradkin , and v.j .
emery , nature * 393 * , 550 ( 1998 ) .
m. ichioka and k. machida , j. phys .
. jpn . * 68 * , 4020 ( 1999 ) .
b. valenzuela , m. vozmediano , and f. guinea , phys .
rev . b*62 * , 11312 ( 2000 ) , and m. vozmediano , personal communication . s. chakravarty , r.b .
laughlin , d.k .
morr , and c. nayak , phys .
rev . b*63 * , 94503 ( 2001 ) ; j. zaanen , cond - mat/0103255 .
kaplan and b.g .
vekhter ,  cooperative phenomena in jahn - teller crystals " ( plenum , n.y . , 1995 ) , p. 216 .
j. hafner , `` from hamiltonians to phase diagrams '' ( springer , berlin , 1987 ) .
anderson , _ `` van hove singularities as a source of anomalies in the cuprates '' _ , pup preprint 414 ( 1994 ) , and `` the theory of high - t@xmath227 superconductivity '' , ( princeton , university press , 1997 ) .
yi , z .- g .
bishop , and j.t .
gammel , phys .
rev . b*58 * , 503 ( 1998 ) .
t. tohyama , s. nagai , y. shibata , and s. maekawa , phys .
lett . * 82 * , 4910 ( 1999 ) .
markiewicz , phys .
rev . b*56 * , 9091 ( 1997 ) .
p. fazekas , foundations of physics * 30 * , 1999 ( 2000 ) . c. kusko and r.s .
markiewicz , unpublished , cond - mat/0102440 .
r. hlubina , s. sorella , and f. guinea , phys .
lett . * 78 * , 1343 ( 1997 ) ; r. hlubina , phys . rev . b*59 * , 9600 ( 1999 ) .
see also d. vollhardt , n. blmer , k. held , m. kollar , j. schlipf , and m. ulmke , z. phys .
b*103 * , 283 ( 1997 ) .
markiewicz , j. phys .
sol . * 58 * , 1179 ( 1997 ) .
j. fritzenktter and k. dichtel , in `` stripes and related phenomena '' , ed . by a.
bianconi and n. saini , ( kluwer academic / plenum , n.y . , 2000 ) , p. 369 .
visscher , phys .
rev . b*10 * , 943 ( 1974 ) ; e.l .
nagaev ,  physics of magnetic semiconductors " ( moscow , mir , 1983 ) .
g. su , phys .
rev . b*54 * , 8281 ( 1996 ) ; a.c .
cosentini , m. capone , l. guidoni , and g.b .
bachelet , phys .
rev . b*58 * , 14685 ( 1998 ) .
putikka , r.l .
glenister , and r.r.p .
singh , phys .
lett . * 73 * , 170 ( 1994 ) ; c.t .
shih , y.c .
chen , and t.k .
lee , phys .
rev . b*57 * , 627 ( 1998 ) ; and m. calandra , f. becca , and s. sorella , phys .
lett . * 81 * , 5185 ( 1998 ) .
c.s . hellberg and e. manousakis , ( a )
lett . * 78 * , 4609 ( 1997 ) , ( b ) phys
. rev . b*61 * , 11787 ( 2000 ) .
l. p. pryadko , s. kivelson and d. w. hone , phys .
lett , * 80 * , 5651 ( 1998 ) .
h. schulz , j. phys .
( paris ) * 50 * , 2833 ( 1989 ) ; d. poilblanc and t.m .
rice , phys .
rev . b*39 * , 9749 ( 1989 ) ; j. zaanen and o. gunnarsson , _ ibid .
_ b*40 * , 7391 ( 1989 ) ; m. inui and p.b .
littlewood , _ ibid .
_ b * 44 * , 4415 ( 1991 ) ; a. singh and z. teanovi , _ ibid .
_ b*41 * , 614 ( 1990 ) ; k. yonemitsu , a.r .
bishop , and j. lorenzana , phys .
rev . b*47 * , 8065 , 12059 ( 1993 ) ; and d. gra , k. rociszewski , and a.m. ole , phys .
rev . b*60 * , 7429 ( 1999 ) .
vergs , e. louis , m.p .
lpez - sancho , f. guinea , and a.r .
bishop , phys .
rev . b*43 * , 6099 ( 1991 ) .
m. kato , k. machida , h. nakanishi , and m. fujita , j. phys .
. jpn . * 59 * , 1047 ( 1990 ) .
laughlin , cond - mat/9709195 ; j. zaanen , j. phys .
solids * 59 * , 1769 ( 1998 ) .
markiewicz , c. kusko and v. kidambi , phys .
b.*60 * , 627 ( 1999 ) . c. balseiro and l. falicov , phys .
rev . b*20 * , 4457 ( 1979 ) . b. valenzuela and m.a.h .
vozmediano , phys .
b*63 * , 153103 ( 2000 ) .
markiewicz and c. kusko , unpublished , cond - mat/0102452 .
agterberg , v. barzykin , and l.p .
gorkov , europhys . lett .
* 48 * , 449 ( 1999 ) and phys .
rev . b*60 * , 14868 ( 2000 ) .
g. baskaran , cond - mat/0007137 .
markiewicz , unpublished .
feng , n.p .
armitage , k.m .
shen , a. damascelli , c. kim , f. ronning , z .- x .
shen , d.a .
bonn , r. liang , w.n .
hardy , a.i .
rykov , and s. tajima , phys .
lett . * 86 * , 4370 ( 2001 ) .
zhou , p. bogdanov , s.a .
kellar , t. noda , h. eisaki , s. uchida , z. hussain , and z .- x .
shen , science * 286 * , 268 ( 1999 ) .
d. orgad , s.a .
kivelson , e.w .
carlson , v.j .
emery , x.j .
zhou , and z .- x .
shen , phys .
lett . * 86 * , 4362 ( 2001 ) .
m. lindroos and a. bansil , personal communication .
a. kampf and j. schrieffer , phys .
rev . b*41 * , 6399 ( 1990 ) . s .-
c . zhang , science * 275 * , 1089 ( 1997 ) .
klemm , unpublished , and physica c*341 - 8 * , 839 ( 2000 ) .
c. castellani , c. di castro , and m. grilli , phys .
lett . * 75 * , 4650 ( 1995 ) .
markiewicz , physica c*169 * , 63 ( 1990 ) .
t. katsufuji , t. tanabe , t. ishikawa , y. fukuda , t. arima , and y. tokura , phys .
rev . b*54 * , 14230 ( 1996 ) .
p. bourges , cond - mat/0009373 .
p. bourges , l.p .
regnault , j.y .
henry , c. vettier , y. sidis , and p. burlet , physica b*215 * , 30 ( 1995 ) .
j. brinckmann and p.a .
lee , phys .
* 82 * , 2915 ( 1999 ) .
norman , cond - mat/9912203 ; y .- j .
kao , q. si , and k. levin , phys .
rev . b*61 * , 11898 ( 2000 ) .
boin , g.h .
kwei , h. takagi , and s.j.l .
billinge , phys .
lett . * 84 * , 5856 ( 2000 ) .
o. baberski , a. lang , o. maldonado , m. hcker , b. bchner , and a. freimuth , europhys . lett . * 44 * , 335 ( 1998 ) .
t. timusk , a.v .
puchkov , s. doyle , a.m. hermann , and n.n .
kolesnikov , unpublished .
y. tokunaga , k. ishida , k. yoshida , t. mito , y. kitaoka , y. nakayama , j. shimoyama , k. kishio , o. narikiyo , and k. miyake , physica b * 284 - 88 * , 663 ( 2000 ) .
a.w . hunt , p.m. singer , a.f .
cederstrm , and t. imai , cond - mat/0011380 .
loram , j.l .
luo , j.r .
cooper , w.y .
liang , and j.l .
tallon , physica c*341 - 8 * , 831 ( 2000 ) .
shen , presented at m@xmath226s - htsc - vi , houston , tx , feb .
20 - 25 , 2000 ; d.l .
feng , d.h .
shen , c. kim , h. eisaki , a. damascelli , r. yoshizaki , j .-
shimoyama , k. kishio , g. gu , s. oh , a. andrus , j. odonnell , j.n .
eckstein , and z .- x .
shen , science * 280 * , 277 ( 2000 ) .
y. koike , m. akoshima , m. aoyama , k. nishimaki , t. kawamata , t. adachi , t. noji , i. watanabe , s. ohira , w. higemoto , and k. nagamine , int . j. mod
b*14 * , 3520 ( 2000 ) .
c. sire , c.m .
varma , a.e .
ruckenstein , and t. giamarchi , phys .
* 72 * , 2478 ( 1994 ) .
d. belitz and t.r .
kirkpatrick , cond - mat/0106279 .
h. ding , j.r .
engelbrecht , z. wang , j.c .
campuzano , s .- c .
wang , h .- b .
yang , r. rogan , t. takahashi , k. kadowaki , and d.g .
hinks , cond - mat/0006143 .
i. martin , g. ortiz , a.v .
balatsky , and a.r .
bishop , int .
phys , b*14 * , 3567 ( 2000 ) . c. bernhard , j.l .
tallon , th .
blasius , a. golnik , and ch .
niedermayer , phys .
* 86 * , 1614 ( 2001 ) .
tallon , j.w .
loram , and g.v.m .
williams , cond - mat/0008295 ; however , note j. bobroff , w.a .
macfarlane , h. alloul , and p. mendels , cond - mat/0008396 .
j.b . goodenough and j .- s .
zhou , phys .
rev . b*49 * , 4251 ( 1994 ) .
wen , w.l .
yang , and z.x .
zhao , physica c*341 - 8 * , 1735 ( 2000 ) , and h.h .
wen , x.h .
chen , w.l .
yang , and z.x .
zhao , phys .
lett . * 85 * , 2805 ( 2000 ) .
cooper and j.w .
loram , j. phys .
i france * 6 * 2237 ( 1996 ) .
r. krupke and g. deutscher , phys .
lett . * 83 * , 4634 ( 1999 ) ; a. sharoni , g. koren , and o. millo , cond - mat/0103581 .
s. rusiecki , e. kaldis , e. jilek , and c. rossel , j. less comm .
metals * 164 * , 31 ( 1990 ) .
d. allender , j. bray , and j. bardeen , phys .
rev . b*7 * , 1020 ( 1973 ) .
uemura , physica c*341 - 8 * , 2117 ( 2000 ) .
uemura , a. keren , l.p .
luke , w.d .
wu , y. kubo , t. manako , y. shimakawa , m. subramanian , j.l .
cobb , and j.t .
markert , nature * 364 * , 605 ( 1993 ) .
markiewicz and b.c .
giessen , physica c*160 * , 497 ( 1989 ) .
j. corson , r. mallozzi , j. orenstein , j.n .
eckstein , and i. bozovic , nature * 398 * , 221 ( 1999 ) .
ong , y. wang , t. kakeshita , and s. uchida , nature * 406 * , 486 ( 2000 ) . c. bernhard , d. munzar , a. golnik , c.t .
lin , a. wittlin , j. humlek , and m. cardona , phys .
rev . b*61 * , 618 ( 2000 ) . c. meingast , v. pasler , p. nagel , a. rykov , s. tajima , and p. olsson , phys .
lett . * 86 * , 1606 ( 2001 ) .
g. baskaran , mod .
b*14 * , 377 ( 2000 ) .
b. batlogg and v.j .
enery , nature * 382 * , 20 ( 1996 ) .
v. barzykin and d. pines , phys .
b*52*,13585 ( 1995 ) ; j. schmalian , d. pines , and b. stojkovic , phys .
lett . * 80 * , 3839 ( 1998 ) .
b. batlogg , h.y .
hwang , h. takagi , r.j .
cava , h.l .
kao , and j. kwo , physica c*235 - 240 * , 130 ( 1994 ) ; j.l .
tallon , j.r .
cooper , p.s.i.p.n .
de silva , g.v.m .
williams , and j.w .
loram , phys .
lett . * 75 * , 4114 ( 1995 ) .
m. ido , n. momono , and m. oda , j. low .
* 117 * , 329 ( 1999 ) .
v.j . emery and s.a .
kivelson , nature * 374 * , 6521 ( 1995 ) .
markiewicz , j. phys .
* 2 * , 665 ( 1990 ) ; a. himeda and m. ogata , phys .
lett . * 85 * , 4345 ( 2000 ) .