the charge and magnetic instabilities of one dimensional materials are well understood in terms of fermi surface ( fs ) nesting .
here it is shown that for two - dimensional materials , features in maps of the bare susceptibility can be understood in terms of a series of fs ` nesting curves ' , and the dominant instabilities are generally related to _
double nesting _ features .
the analysis shows how to locate these nesting features in momentum ( @xmath2 ) space , often providing analytical expressions . as an application ,
the full evolution with doping of the leading magnetic instabilities for several families of cuprates is presented , both in the conventional hartree - fock ( hf ) plus rpa ( hf+rpa ) and in a gutzwiller approximation ( ga+rpa ) calculation .
the analysis provides a pseudogap candidate and makes the surprising prediction that the ` checkerboards ' seen in scanning tunneling microscopy ( stm ) studies are _ not _ the same phase as the ` stripes ' in la@xmath3sr@xmath4cuo@xmath5 ( lsco ) .
the present results should find extensive utilization .
first , for weakly correlated two - dimensional systems the hf+rpa results provide an essentially complete solution to the nesting problem , as long as the interaction @xmath0 is @xmath2-independent .
the ga+rpa extends the results into the intermediate coupling regime . in both cases ,
the magnetic instabilities are determined by zeroes of the stoner denominator , @xmath6 here for cuprates @xmath7 is the hubbard @xmath0 in hf+rpa , and a more complicated @xmath2-dependent object @xmath8 in the ga+rpa calculation .
thus , the leading hf+rpa instability is simply associated with the maximum of the bare susceptibility @xmath9 , while the leading gutzwiller instability can be shifted by the @xmath2-dependence of @xmath8 .
it will be clear that the same analysis can be extended to any two - dimensional material .
the cuprates appear to be in an intermediate coupling regime where the gutzwiller results can be expected to provide a good approximation to the phase diagram at @xmath10 .
thus , in the cuprates , recent quantum monte carlo ( qmc)@xcite and ` quasiparticle - gw ' ( qp - gw)@xcite calculations have been able to reproduce experimental arpes spectra of optimally and overdoped cuprates , starting essentially from lda bands and calculating the self - energy self consistently . in the qp - gw approach
, the self energy is calculated as a convolution over a renormalized rpa susceptibility . not only is the low - energy dispersion renormalization reproduced , but also the ` waterfall ' effect@xcite , which represents the dressing of low - energy quasiparticles by ( mainly ) spin fluctuations .
extension of these results to the hole - underdoped regime will require identification of the phase or phases responsible for the pseudogap .
since these are most likely to be incommensurate density wave or ` stripe ' phases , the qmc calculations have a severe problem of limited @xmath2-resolution , while the qp - gw calculations are ideally suited to handle this .
such an approach has had considerable success with electron - doped cuprates , where the leading instability is commensurate @xmath11 .
( we set the lattice constant @xmath12 . )
the paper is organized around the stoner criterion as follows .
sections 2 - 3 describe the calculations of the zero frequency bare susceptibility @xmath13 , and determining its maximum @xmath14 . in section 2
we introduce the concept of nesting curve , associated with the nesting criterion @xmath15 , where @xmath16 is the ( anisotropic ) fermi wave vector , and we demonstrate that @xmath14 is generally associated with a double nesting .
we apply this concept to the cuprates , taking the hopping parameters from tight - binding one - band fits to the lda dispersions , and show that over most of the doping range there is a unique @xmath14 , with a surprising electron - hole symmetry . in section 3 we extend this analysis to the more complicated hole doping regime for doping @xmath17 between half filled and the doping of the van hove singularity ( vhs ) . here
competing magnetic phases are found , and the dominant phase is sensitive to material parameters , being different for different cuprates .
section 4 describes the ga+rpa technique and introduces the corresponding @xmath18 .
section 5 presents the resulting gutzwiller magnetic phase diagrams .
it is found that the phase diagram of la@xmath3sr@xmath4cuo@xmath5 ( lsco ) is distinct from that of most other cuprates , but that for all cuprates , using a bare hubbard @xmath19 , the paramagnetic state is unstable at the ga level for all relevant dopings , including the overdoped regime .
a discussion is presented in section 6 , and conclusions in section 7 .
the present calculations are based on tight - binding parametrizations of typical dispersions for single - layer cuprates , including models of la@xmath3sr@xmath4cuo@xmath5 ( lsco ) , nd@xmath3ce@xmath4cuo@xmath5 ( ncco ) , and bi@xmath20sr@xmath20cu@xmath21o@xmath22 ( bi2201 ) , table i in appendix  [ sec : band - parameters ] . for ncco , bi-2201(2 ) , and lsco(2 ) ,
the parameters are based on a tight - binding fit to the lda bands@xcite , while sets bi-2201(1 ) and lsco(1 ) are direct fits to the experimental bands . in all cases ,
@xmath23 dispersion is neglected , approximating the cuprates as two - dimensional .
0.5 cm    the real part of the susceptibility @xmath24 is dominated by a series of plateaus , with the largest susceptibility systematically shifting from one plateau to another as a function of doping .
figure  [ fig:1 ] shows that these and related features dominate @xmath25 over the full doping range , and further reveals a striking quasi - electron - hole symmetry of the evolution . however , instead of being symmetric about @xmath26 , the evolution is symmetric about the doping of the van hove singularity ( vhs ) , @xmath27 , where the fermi energy coincides with the vhs , @xmath28 , and the density of states ( dos ) diverges logarithmically .
figure  [ fig:1](a ) shows @xmath25 calculated along the high symmetry axes of lsco , for the full electron - doping range @xmath29 from 1 to 2 , with an extension to the hole doping of the vhs ( blue curves ) . near half filling the susceptibility
is dominated by the well - known plateau@xcite near @xmath30 ( red and blue curves ) , which peaks at the vhs . as electron - doping increases , the peak shifts to a second plateau near @xmath31 ( green curves ) , then to a third near @xmath32 for a nearly full band ( brown curves ) .
the _ same _ sequence is followed for hole doping , fig .  [ fig:1](b ) .
in addition to these plateaus , there is an additional @xmath32-centered feature , which is prominent near the doping of the vhs , but has largely disappeared by half filling .
this feature will be referred to as the antinodal nesting ( ann ) plateau .
bi2201 displays the same plateaus , with the same quasi - electron - hole symmetry , fig .
[ fig:1](c , d ) .
a key difference is that the @xmath30-plateau of bi2201 is convex over a considerably wider doping range than in lsco , before turning concave near the vhs .
this has a consequence that the ann peak dominates over a wide doping range in bi2201 , but only in the immediate vicinity of the vhs in lsco . at the vhs ,
the bare susceptibility diverges logarithmicly both at @xmath32 ( where it is equal to the dos ) and at @xmath30 , but the latter divergence is quite weak and not apparent in the numerical calculations of fig .
1 . despite the apparent complexity of these evolving susceptibility patterns ,
the plateau evolution can be understood in detail , with analytic formulas for the positions of all dominant peaks .     and folded back into the first bz ( dashed ) .
( b ) electron - like fermi surface ( solid ) shifted by the double nesting vector @xmath33 and folded back into the first bz ( dashed ) .
the ( double ) nesting points are marked by a solid dot .
( c , d ) construction of nesting curves ( dashed ) for hole - like ( c ) and electron - like ( d ) fs s ( solid lines ) from condition eq .
[ eq : g2 ] .
the fs is shown in the first bz whereas the nesting curves are folded into the momentum space @xmath34 as defined by the dashed square .
the arrows indicate scattering processes which lead to ` double nesting ' as explained in the text .
the dotted lines corresponds to the boundary of the magnetic bz in the first and enlarged zone .
( a , c ) parameters for ncco , @xmath35 ; ( b , d ) parameters for lsco(1 ) , @xmath36.,title="fig:",width=340 ]   and folded back into the first bz ( dashed ) .
( b ) electron - like fermi surface ( solid ) shifted by the double nesting vector @xmath33 and folded back into the first bz ( dashed ) .
the ( double ) nesting points are marked by a solid dot .
( c , d ) construction of nesting curves ( dashed ) for hole - like ( c ) and electron - like ( d ) fs s ( solid lines ) from condition eq .
[ eq : g2 ] .
the fs is shown in the first bz whereas the nesting curves are folded into the momentum space @xmath34 as defined by the dashed square .
the arrows indicate scattering processes which lead to ` double nesting ' as explained in the text .
the dotted lines corresponds to the boundary of the magnetic bz in the first and enlarged zone .
( a , c ) parameters for ncco , @xmath35 ; ( b , d ) parameters for lsco(1 ) , @xmath36.,title="fig:",width=340 ] 0.5 cm    the generic evolution of the plateaus from @xmath37 can be understood with reference to the ` nesting curves ' , fig .  [
fig : new ] . for the generic case of two fermi surface segments
, a nesting curve can be defined as the locus of all points @xmath38 , where @xmath39 is a point on the @xmath40th fs , fs@xmath41 , with the restriction that when fs@xmath21 is shifted by @xmath2 it is tangent to fs@xmath20 . for the cuprates , there is usually just a single fs section [ an exception is given below in section iii , fig .  6 ] , and the nesting curves simplify to plots of @xmath42 for any fermi momentum @xmath16 , fig .
[ fig : new ] . since the point @xmath30 lies on horizontal and vertical planes of reflection symmetry for the susceptibility ,
the nesting curves must be supplemented by their reflections about the lines @xmath43 , @xmath44 .
[ fig:2 ] shows three sets of nesting curves corresponding to three hole dopings for lsco ( dispersion 1 ) , @xmath17 = 0.41 ( red curves ) , 0.62 ( blue ) , and 0.79 ( green ) . frames ( b - d ) show the corresponding susceptibility maps @xmath45 , and it can be seen that the ridges in @xmath25 are exactly given by the nesting curves . [ for convenience , the @xmath36 data is replotted in fig .
note that the susceptibility maps are plotted over the range of @xmath46 between 0 and @xmath47 , whereas the nesting curves are plotted over the wider range 0 to @xmath48 . ] furthermore , the dominant peaks in @xmath25 correspond to the intersection of two nesting curves . by drawing the original and @xmath2-shifted fermi surfaces ( fss ) , figs .
[ fig : new](a ) ,  ( b ) , it can be seen that this overlap corresponds to the simultaneous nesting of two different sections of fs .
hence the term ` double nesting ' .    as shown in figs .
[ fig : new](c ) ,  ( d ) this kind of ` double nesting ' can originate from either the scattering between points on opposite segments ( cf .
example in fig .
[ fig : new](c ) ) or adjacent segments ( cf .
example in fig .
[ fig : new](d ) ) of the fs . in the example shown in fig .
[ fig : new](c ) we denote the scattered states on the fs as @xmath49 and @xmath50 which yield the points on the nesting curves @xmath51 and @xmath52 . thus in this case
` double nesting ' ( @xmath53 ) occurs when @xmath54 which generally can only be fulfilled for points near the antinodes of hole - like fs s .
these are thus referred to as antinodal nesting ( ann ) features .
more common is the situation sketched in fig .
[ fig : new](d ) where we consider the scattered states @xmath49 and @xmath55 . in this case
the ` double nesting ' condition @xmath53 can only be fulfilled for @xmath56 ( which is trivial since initial and final states of the two scattering processes are identical ) or @xmath57 .
the latter condition implies that this kind of ` double nesting ' generally affects the nodal states and leads to scattering vectors @xmath2 close to @xmath58 , generating the @xmath30-plateau .
a slightly different point of view might help clarify the role of double nesting .
tangency of two fermi surface segments [ ` nesting ' ] means a stability of fs overlap , in that the nesting persists if one displaces the nesting vector in a particular direction .
then double nesting means tangency along two surface segments , so nesting persists if one displaces the nesting vector in two ( rather than one ) directions .
0.5 cm    the origin of the plateau transitions is now apparent : as hole doping increases , the fs shrinks to a small pocket near @xmath32 before disappearing .
the nesting curves shrink in a similar fashion , but with a doubled radius , since @xmath15 .
thus the dominant overlap shifts from near @xmath59 at the vhs ( fig .  2(b ) ) towards @xmath31 ( fig .
2(c ) ) , and finally towards @xmath32 , in fig .
2(d ) , thereby explaining the plateau evolution . for electron doping ,
the fermi surface ultimately closes at @xmath30 , leading to the same sequence of plateau transitions , as illustrated in fig .
[ fig:3 ] , where the nesting maps are directly superposed on the susceptibility curves . in the doping range relevant to the cuprates , the physics of ncco is dominated by the @xmath30 plateau , fig .
[ fig:3]a and red curves in fig .
1(a ) , which shrinks to a point at the end of the ` hot - spot ' regime .
0.5 cm    in all cases , the dominant peak lies along a high symmetry axis , and the doping dependence of its position can readily be found from the dispersion @xmath60 - 4t'c_x(a)c_y(a ) \nonumber \\
-2t''[c_x(2a)+c_y(2a ) ] \nonumber \\
-4t'''[c_x(2a)c_y(a)+c_y(2a)c_x(a ) ] \ > ,
\label{eq:2}\end{aligned}\ ] ] where @xmath61 and @xmath62 is an integer ( or half - integer ) . for this dispersion ,
the vhs is generally at @xmath31 , or @xmath63 .
the positions of the peaks in fig .
[ fig:2](b , c ) and fig .
[ fig:3](b , c ) lie along the zone boundary ( @xmath64 ) with @xmath65 , \label{eq:1x}\ ] ] with @xmath66 $ ] and @xmath67/[4t'']$ ] . for electron [ hole ] doping
the peak is exactly at @xmath31 when @xmath682(t-2t''')$ ] . beyond this point ,
the peak is at @xmath69 , @xmath70 given by eq .
[ eq:1x ] with @xmath712t'+2t''')/[4(t''-[+]2t''')]$ ] and @xmath722(t-2t''')]/[4(t''-[+]2t''')]$ ] .
since a peak at @xmath32 corresponds to ferromagnetism , the above reproduces the common finding that a nearly empty or nearly full band tends to be ferromagnetic .
the ann peak @xmath53 in fig .
[ fig : new]a satisfies @xmath73/[2(t'-2t'')]$ ] .
thus the peak susceptibility is generally associated with double nesting .
the only exception we have found to this is a tendency to remain commensurate for a finite doping range about high symmetry points such as @xmath11 or @xmath74,@xcite  a form of van hove nesting@xcite .
0.5 cm    while the above sequence is completely generic , holding for all cuprates and being electron - hole symmetric , the additional features associated with the vhs are considerably more variable .
[ fig:4 ] illustrates the low - hole doping regime for three dispersions proposed for the cuprates . at half filling ,
all three have a peak susceptibility on the @xmath30 plateau , associated with conventional staggered antiferromagnetism ( saf ) . at low doping the peak is either commensurate or @xmath75 incommensurate , with @xmath76 increasing with doping . for
finite hole doping a new feature emerges , a roughly @xmath77-shaped peak along the zone diagonal , at @xmath78 . from the nesting curves
, it can be seen that the peak is associated with nesting of the flat sections of the bands near @xmath31  hence the name antinodal nesting ( ann )  but the largest susceptibility lies along the zone diagonal , where both @xmath31 and @xmath79 nesting occur simultaneously , fig .
[ fig : new]a .
this ann feature has an interesting relation with ` hot spot ' physics@xcite .
a ` hot spot ' is a point of the fermi surface which is simultaneously on the antiferromagnetic ( af ) zone boundary [ diagonal of the paramagnetic brillouin zone ] . from fig .
[ fig : new](c ) it becomes apparent ( cf . large dot ) that the image of the af zone boundary in the extended @xmath80-map intersects the ` nesting curves ' exactly at the points of ` double nesting ' .
but the image of the af zone boundary gets folded in the @xmath2-map onto the zone diagonal , @xmath81 , thereby generating the usual ann feature .
thus all the diagonal ann peaks in @xmath2 arise from hot spots in @xmath82 . as the doping of the vhs
is approached the hot spots move toward @xmath83 and the nesting curve moves to @xmath32 . since the susceptibility at @xmath32 is equal to the dos , it diverges at the vhs , thereby controlling the magnetic instability .
the appearence of the ann peak leads to a competition between two different kinds of magnetic order , and the doping evolution of the susceptibility maps diverges . for most dispersions
studied , including the left and central columns of fig .
[ fig:4 ] , the ann intensity grows and the dominant peak changes discontinuously from the saf plateau to ann nesting . for the bi-2201 dispersion ( left column )
this ann peak becomes dominant at about @xmath17=0.2 , and evolves smoothly to @xmath32 at @xmath17=0.4 .
if one unfolds the nesting curve , one sees that this happens exactly at the doping of the vhs , when the fs passes through @xmath31 .
the central column , corresponding to an extreme dispersion proposed for lsco to enhance antinodal nesting , displays a much more complicated nesting map over a limited doping range close to the vhs , but also displays a dominant peak at @xmath32 exactly at the vhs .
this map will be explained below .
finally , for the right column , corresponding to a more conventional dispersion for lsco , the ann peak is weaker than the saf peak except in the immediate vicinity of the vhs . in this case , the dominant susceptibility peak remains commensurate at @xmath30 until the doping of the vhs , jumps to @xmath32 at the vhs , then jumps back and smoothly evolves to @xmath84 incommensurate [ this sequence can also be followed in fig .
1(a),(b ) ] .
note further , in fig .  5(i , j ) , that there is a wide doping range where the susceptibility peak remains commensurate , at @xmath30 .
[ since @xmath30 lies along several mirror planes , the corresponding susceptibility is in general a local maximum or minimum . ]    0.5 cm    the complicated nesting map of the middle column near the vhs is explained in fig .
[ fig:5 ] .
the dispersion is such as to produce an ` extended vhs',@xcite which first intersects the @xmath85 axis at a point @xmath86 . for larger hole doping , the fs has two sheets , one a squarish barrel centered at @xmath32 , the other a pocket centered at @xmath31 , fig .
[ fig:5](b ) . the intrasheet
@xmath15 nesting maps are shown by the red curves in fig .
[ fig:5](a ) .
in addition to the one associated with the barrel fs , there are two overlapping segments associated with the pockets , but translated by the @xmath2 folding to the vicinity of @xmath32 . as in the other frames of figs .  2 - 4
, these nesting curves exactly match some of the ridges seen in the susceptibility . however , there are additional ridges , associated with nesting between two different fs segments .
thus , the blue curves in fig .
[ fig:5](a ) represent barrel - pocket nesting  that is , the locus of the @xmath2-vectors needed to shift the pocket until it is tangent to some point on the barrel fs . for instance , the points labelled 1 ,  2 ,  3 represent translations equivalent to those shown in fig .
[ fig:5](b ) .
the new nesting curves can easily be found numerically , by requiring that the two fs sections have a common tangent at the point of osculation .
similarly , the green curve near @xmath30 represents inter - pocket nesting , with point 4 illustrating one particular nesting vector .
it can be seen that the full collection of nesting curves explains all of the ridge - like features seen in the susceptibility map , and in particular allow the determination of the points of maximal susceptibility , except for the above - noted commensurability effects .
note however , that despite these complications , the susceptibility peak moves to @xmath32 in a finite doping range about the point where the antinodal electron pocket shrinks to zero width .
the above discussion can be summarized : in the doping between half filling and the vhs , a new susceptibility peak arises , associated with nesting of the antinodal parts of the large fs .
two kinds of behavior are found : when the ann is dominant , the peak evolves to @xmath32 at the vhs . on the other hand , for lsco(2 ) with small @xmath87 ,
the ann peak is inherently weaker than the @xmath30 peak , in which case the susceptibility remains commensurate at @xmath30 from half filling almost to the vhs , then smoothly develops a @xmath88 incommensurability . yet
even here , in the immediate vicinity of the vhs , the ann peak at @xmath32 becomes dominant in a very limited doping range .
hence , for most dispersions the susceptibility will have strong fm fluctuations near the vhs , while for other dispersions the fluctuations remain mostly afm .
in the ga+rpa calculation@xcite , the electronic paramagnetic ground state energy is calculated in the ga@xcite and then expanded to second order in the on - site and intersite magnetic fluctuations in the spirit of vollhardt s fermi liquid approach@xcite .
response functions are computed using linear response in the presence of small external field.@xcite one obtains rpa like susceptibilities but with strong vertex corrections .
as shown in appendix  [ sec : spin - rotat - invar ] longitudinal and transverse susceptibilities are trivially related and lead to the same stoner criteria as required by spin rotational invariance . in terms of a tensor bare transverse susceptibility ( cf . appendix  [ sec : spin - rotat - invar ] and appendix c ) @xmath89 with @xmath90 , the ga+rpa dressed susceptibility @xmath91 is given by @xmath92 here the ratio of dressed to bare dispersion is given by @xmath93 , with the gutzwiller renormalization factor @xmath94 ^ 2}{1-x^2}\ ] ] which depends on the ga double occupancy variational parameter @xmath95 and the doping @xmath17 .
the interaction matrix is@xcite @xmath96 which is defined in appendix  [ sec : spin - rotat - invar ] .    while eq .
[ eq : g4 ] is a tensor equation , it can be expanded into the form of eq .
[ eq : g1 ] with @xmath97 and @xmath98 , @xmath99 with @xmath100 , @xmath101 .
the @xmath102 factors appear now in @xmath103 since we want to use the bare @xmath104 in eq .
[ eq : g1 ] .
all correlation effects are incorporated in the definition of @xmath103 which therefore can be viewed as a vertex corrected interaction term in the magnetic @xmath105 channel .
note that @xmath103 is not strictly a pure interaction term , but is weighted by kinetic terms which enhance its @xmath2-dependence .
0.5 cm    figure  [ fig:6 ] compares @xmath106 , @xmath107 , @xmath108 , and the derived @xmath109 , @xmath110  the results are fairly insensitive to doping or band parameters .
the sharp structures associated with nesting show up only in @xmath106 , while the other @xmath111 s are smooth , and the @xmath112 s contain only a weak structure imposed from @xmath106 .
hence @xmath103 remains a smooth , weakly varying function of @xmath2 , and the instabilities are controlled by the peaks in @xmath104 .
figure  [ fig:7](a ) shows how @xmath103 varies as the bare @xmath0 is increased . at small @xmath0 , @xmath113 and the conventional hf+rpa
is recovered .
again the results do not depend strongly on doping ( fig .
8(b ) ) or dispersion .
the general trend is that at large @xmath0 @xmath103 tends to saturate , and the large-@xmath2 components of @xmath103 are reduced most strongly , tending to favor instabilities nearer @xmath114 , suggestive of ferromagnetic domains .
[ fig:7 ] should be compared with fig .  2 of vilk and tremblay@xcite , who find a similar saturation .
since these authors employ a very different perscription for including vertex corrections , the similarity of our results gives considerable further confidence to the trends we find .
in addition since the gutzwiller approximation has a simple interpretation in terms of kinetic energy suppression due to correlation our results shed further light on the physical meaning of the vertex corrections
.    0.5 cm
0.5 cm    figure  [ fig:8 ] shows our main result , the magnetic phase diagrams of lsco and bi2201 calculated in the gutzwiller ga+rpa approach , based on the stoner criterion , eq .  1 .
the phase diagram for bi2201 is fairly generic , and we have found similar results for ncco , bi2212 [ neglecting bilayer splitting ] and scoc , while the lsco phase diagram is limited to the dispersion parameters of lsco(1 ) .
thus , for electron doping the phase diagram is very simple , dominated by simple af order with @xmath2 very close to @xmath11 , while for hole doping there is a competition between the saf and ann order , except in lsco , where saf order dominates .
it is possible that ann order is relevant to the ` checkerboard ' phase seen in stm studies of several cuprates , as discussed briefly below .
results for hf+rpa are similar , but since @xmath7 in eq .  1 is not bounded ( @xmath115 ) , there would be magnetic instabilities for all dopings , with transitions to the @xmath31 and @xmath32 plateaus , as in fig .
[ fig:1 ] .
in contrast , in the ga+rpa calculation @xmath103 saturates , and magnetic order exists only in a limited doping range in the cuprates , primarily in the range where the @xmath30 plateau exists .
however , this doping range encompasses the full doping range of relevance to cuprate physics  including the overdoped regime  if a doping independent @xmath116 is assumed .
thus , to have a magnetic quantum critical point near optimal doping , one must go beyond the one - band ga+rpa approach .
the limited doping range for the magnetic phases has also been seen in a @xmath117 study.@xcite analogously to hf+rpa , ga+rpa determines the transition line at which the mean field state becomes unstable , but the energies are evaluated at mean field level .
fluctuation contributions to the energy will most likely lower the energy of the paramagnetic ( disordered ) phase further reducing the stability range of the magnetic phases .
there have by now been countless calculations concerning nanoscale phase separation or stripe phases in cuprates , and this is becoming a paradigm in many other correlated electronic systems as well .
we believe , however , that the present results are unique in providing a systematic phase diagram , covering the full doping dependence@xcite and all possible @xmath1-vectors for realistic band dispersions , and indeed finding that the dominant instabilities may be different in different cuprates .
while the above has provided a thorough analysis of possible magnetic instabilities , nothing has been said about competing instabilities in the charge channel  charge density waves ( cdws ) or ` stripes ' .
of course , one possibility is that charge order arises as a secondary effect following the spin order , with @xmath118 .
the case where the charge instability is primary will be somewhat harder to analyze .
that is because the corresponding stoner factor analysis is likely to require an extension beyond the hubbard model , and hence involve greater uncertainty in the choice of @xmath119 . for charge instabilities ,
the corresponding stoner denominator is nothing but the ( zero frequency ) dielectric constant @xmath120 , where @xmath121 is a background dielectric constant and @xmath122 is a coulomb potential .
when @xmath122 is taken as either the hubbard @xmath0 or a long - range coulomb interaction,@xcite @xmath123 is always positive and there is no instability .
moreover , from the theory of dielectric stability , it is known that @xmath123 can not fall in the range between 0 and 1 [ equivalently @xmath124 , and an instability @xmath125 must be approached through negative values of @xmath123.@xcite while a purely electronic instability could still arise via inclusion of local field effects , the most natural situation arises when electron - phonon coupling is included .
a simple @xmath126 wave instability will be suppressed by the on - site coulomb repulsion as shown in ref .
but one can still have instabilities in other channels .
to the extent that the matrix elements which determine the channel have a smooth behavior one will have a stoner denominator similar to eq .
( [ eq : g1 ] ) with @xmath127 , where @xmath128 is the electron - phonon coupling parameter , @xmath129 the corresponding bare phonon propagator , @xmath130 a bare phonon frequency , and @xmath131 an electron - electron dielectric constant.@xcite    thus , since the instability is also controlled by a stoner factor , the results of the present paper will also apply for charge instabilities .
this is consistent with the common expectation that cdws and spin density waves ( sdws ) are controlled by the same nesting instabilities
. stated differently , the stoner criterion is a formal expression of the idea that @xmath132 for an instability , where @xmath133 , the interaction energy , is represented by @xmath7 while the kinetic energy @xmath134 is measured via @xmath135 .
the kinetic energy is closely tied to the band structure , and carries important material - specific information .
in contrast , the potential energy is fairly featureless , leading to smoothly varying @xmath119 .
thus , the locations of the instabilities are controlled by peaks in @xmath104 , and a @xmath2-dependent @xmath0 can only shift the dominant instability between two peaks of comparable height . in the remainder of this section ,
when we discuss comparison to experiment , it will be seen that the present model has correctly determined the dominant @xmath2 values , but that in several cases experiment points more toward cdws than sdws .
there have been a number of recent hints that the ` stripe ' order in underdoped la@xmath3a@xmath4cuo@xmath136 , a = sr ( lsco ) or ba ( lbco ) is not the same phenomenon as the ` checkerboard ' order found in sr@xmath20cuo@xmath20cl@xmath20 ( scoc ) and bi@xmath20sr@xmath20cacu@xmath20o@xmath137 ( bi-2212 ) .
thus , resonant soft x - ray scattering experiments find evidence of charge order in lbco ( stripes)@xcite but not in ca@xmath20cuo@xmath20cl@xmath20 ( checkerboards)@xcite , while evidence for time reversal symmetry breaking has been found in yba@xmath20cu@xmath138o@xmath139 ( ybco ) , but not in lsco@xcite .
the present results suggest a connection between the ann and checkerboard phases , in that the periodicity of the latter also scales with the antinodal nesting vector@xcite , and a further connection between the @xmath30 plateau and stripes .
thus it is quite interesting to observe that the ann phase is virtually absent for the lsco(2 ) dispersion , where evidence for conventional stripes is strongest .      in 1d systems
the susceptibility diverges at @xmath15 at all dopings , due to perfect ( flat band ) nesting .
the present results generalize this to 2d systems : the dominant instability is generally associated with double nesting , but since @xmath104 in general does not diverge , a finite coupling @xmath0 is required to drive an instability .
an exception is vhs nesting@xcite , for which @xmath104 has a logarithmic divergence . at the vhs , there are two competing instabilities , af at @xmath30 and fm at @xmath32 , corresponding to the magnetic branch of the @xmath140 phase diagram of the cuprates , and these are responsible for the commensurate pinning near @xmath30 and @xmath32 , respectively.@xcite doped away from the vhs , these instabilities evolve into the two dominant peaks of the susceptibility , and can be considered as ` generalized vhs nesting ' . as
such , they dominate the magnetic physics of the cuprates over the full doping range , coextensive with the limits of the @xmath141-plateau in the susceptibility .
recently , storey _ et al._@xcite proposed that the generic behavior of high - t@xmath142 cuprates could be understood if the pairing interaction ( or pairing energy cutoff ) falls off rapidly near the vhs .
a recent experiment@xcite does indeed confirm that @xmath143 scales with @xmath144 , but @xmath145 .
the present results suggest that most cuprates will have strong fm fluctuations near the vhs , which are incompatible with simple d - wave superconductivity .
empirical evidence for fm fluctuations has been noted previously@xcite .
while some previous calculations@xcite have found evidence of ferromagnetism near a vhs , others@xcite have suggested that fm instabilities were unlikely in competition with incommensurate susceptibility .
the present calculations confirm that a dominant fm susceptibility should be present near the vhs .
the advantage of two - dimensional materials is that it is straightforward to display the susceptibility maps and nesting curves . in 3d ,
the curves become _ nesting surfaces _ and the susceptibility maps are 4-dimensional .
one can speculate that the dominant susceptibility peaks correspond to triple - nesting points .
for quasi-2d materials , it should be possible to analyze a series of cuts perpendicular to the ( weakly - dispersing ) @xmath146-axis .      in principle
, the present results could provide information on nanoscale phase separation ( nps ) as well .
one model of nps is that there are different instabilities associated with particular dopings [ e.g. , half filling and the vhs ] , and that these two phases are more stable than uniform phases at intermediate dopings . in this case , the two end phases could still be described by nesting maps , only at particular dopings .
for instance , in the hf+rpa analysis of lsco(2 ) , the fm phase is stable only very close to the vhs , and could lead to nps with a second phase at the undoped insulator .
the stoner criterion determines which @xmath2 value is most unstable , and the minimum @xmath0 needed to drive that instability .
however , as @xmath0 increases above threshold , the @xmath2 of the ordered phase may shift .
thus , at half filling , when a full gap can be formed , @xmath2 will no longer be determined by best nesting , but by the largest gap .
this tends to favor more commensurate @xmath2 values , leading to a pinning of @xmath2 at these commensurate values over a wide range of parameters .
we find that as @xmath0 increases , there is a first order transition from an incommensurate phase with fermi surface pockets to a commensurate [ @xmath30 or @xmath31 ] phase that is fully gapped.@xcite away from half filling , luttinger s theorem ensures the persistence of a fermi surface , so nesting instabilities should persist over a wider range of @xmath0 s .
furthermore , as we have seen above , the @xmath2-value corresponding to the largest @xmath148 can shift with @xmath0 .
typically , for @xmath103 the shift is to a smaller @xmath2-value , associated with an instability towards fm order ( in the magnetic channel ) or nps ( in the charge channel ) .
the present results provide a constructive scheme for identifying the dominant nesting instabilities for any two dimensional material .
clearly , for fermi surfaces with multiple sections , a large number of nesting curves are possible , leading to extremely complicated susceptibility maps .
nevertheless , the present scheme will automatically sort out the possible double nesting peaks and follow their evolution with doping .
this should allow a much more detailed understanding of 2d phase diagrams , particularly for magnetic phases , where the interaction @xmath0 has negligible @xmath2-dependence .
the good agreement of the ga+rpa calculations with more detailed variational calculations at half filling@xcite and with exact and numerical results in @xmath117,@xcite and of the af+sc model in electron doped cuprates with experiments gives us confidence in the model .
accordingly , we note the following three points . first , within ga+rpa
, the paramagnetic state is unstable over the full hole - doping range in the cuprates , including overdoped . to avoid this conclusion , and restore a fermi liquid phase in the overdoped regime
, it may be necessary to include non gaussian fluctuations@xcite or an additional doping dependence for @xmath0 .
the origin of any such doping dependence lies outside the ga treatment of the hubbard model .
second , the @xmath11 phase , or its incommensurate extension , is unstable against a competing ann order in [ most ] hole - doped cuprates .
third , there is a material dependence to the phase diagram , and lsco may have a very different doping dependence than other cuprates .
the results have been applied to a number of model dispersions for the cuprates . a possible pseudogap candidate has been identified and a distinction between stripes and checkerboards proposed .
these findings will be discussed in greater detail in ensuing publications .
this work is supported by the us department of energy , office of science , basic energy sciences contract de - fg02 - 07er46352 , and benefited from the allocation of supercomputer time at nersc , northeastern university s advanced scientific computation center ( ascc ) .
rsm s work has been partially funded by the marie curie grant piif - ga-2008 - 220790 soqcs , while gs work is supported by the vigoni program 2007 - 2008 of the ateneo italo - tedesco deutsch - italienisches hochschulzentrum .
.i . band parameter sets [ cols="^,^,^,^,^,^",options="header " , ]     in table i , all hopping parameters are given in mev .
the @xmath149 term for lsco(1 ) is the coefficient of a term @xmath150 in eq .
ncco , lsco(2 ) , and bi2201(2 ) data sets are fit to lda band dispersions of the near - fermi level antibonding cuo@xmath20 band , and are appropriate for the gutzwiller analysis .
the data sets lsco(1 ) and bi2201(1 ) are taken from fits to experimental arpes data . to the extent that the bands are renormalized by a @xmath2-independent factor @xmath151 ,
the effective bare susceptibilities @xmath152 can be used in the gutzwiller analysis , while @xmath14 occurs at the same @xmath2 .
recently it has been suggested that the enhanced nesting , such as observed in lsco(1 ) , may be associated with stripe formation@xcite .
here we sketch the ga+rpa formalism and define the relevant quantities . for a review see ref .  .
the susceptibility can be computed in the longitudinal@xcite or in the transverse channel .
@xcite as shown bellow both results are equivalent as dictated by spin rotational invariance on a paramagnetic ( singlet ) ground state .    for later use
we define the density matrix associated with the unprojected slater determinant , @xmath153 , as @xmath154 and as a shorthand we define @xmath155 .
we start from the spin - rotational invariant gutzwiller energy functional for the hubbard model defined in terms of @xmath156 and the double occupancy in the gutzwiller variational state @xmath157 .
as derived e.g. in ref .
the functional reads , @xmath158 where we have defined the spinor operators @xmath159 and the @xmath160-matrix    @xmath161\cos\varphi_i \\ \frac{s_i^+}{2s_i^z}[z_{i\uparrow}-z_{i\downarrow}]\cos\varphi_i & z_{i\uparrow}\sin^2\frac{\varphi_i}{2}+z_{i\downarrow}\cos^2\frac{\varphi_i}{2 } \end{array } \right ) \\ \tan^2\varphi_i&=&\frac{s_i^+s_i^-}{(s_i^z)^2}.\end{aligned}\ ] ]    with the @xmath146 factors given by @xmath162    and for clarity spin expectation values are denoted by @xmath163 , @xmath164 , @xmath165 , and @xmath166 .    in the limit of a vanishing rotation angle
@xmath167 the @xmath160-matrix becomes diagonal and the renormalization factors reduce to those of the standard gutzwiller approximation .    for a homogeneous , paramagnetic system ( @xmath168 , @xmath169 ) the expansion of the gutzwiller energy functional up to second order in the particle - hole excitations reads , @xmath170
the first contribution is the free fermion part , @xmath171 @xmath172 denotes the dispersion of the gutzwiller quasiparticles whereas @xmath173 corresponds to the unrenormalized dispersions , i.e. @xmath174 .    the remaining part in eq .
( [ eq : ega2 ] ) is due to quasiparticle interactions and separates into contributions from the charge ( @xmath175 ) , the longitudinal ( @xmath176 ) and transverse ( @xmath177 ) spin channel . @xmath178 and @xmath179
have been derived in detail in ref .  .
for our purposes we report in the following the expansion in the longitudinal channel : @xmath180 and in the transverse magnetic channel which has been derived in ref .
@xmath181 the relevant fluctuations have been defined as @xmath182    the 11 element of the interaction kernel is given by @xmath183 with @xmath184 here @xmath185 is the brinkman - rice energy .
@xmath186 is similar , but requires separate averages for different components of the energy : @xmath187\bigl<{c_{kx}+c_{ky}}\bigr > -4t'c_{qx}c_{qy}\bigl < c_{kx}c_{ky}\bigr > \nonumber \\ & -&t''[c_{2qx}+c_{2qy}]\bigl<{c_{2kx}+c_{2ky}}\bigr > \nonumber \\ & -&2t'''[c_{qx}c_{2qy}+c_{qy}c_{2qx}]\bigl<{c_{kx}c_{2ky}+c_{ky } c_{2kx}}\bigr>,\nonumber \label{eq : g17}\end{aligned}\ ] ] in an obvious notation .
furtheron we define the susceptibility matrices for the bare time - ordered correlation functions both in the longitudinal @xmath188 and in the transverse channel @xmath189 where a hat has been added to distinguish fluctuations ( @xmath190 ) form operators ( @xmath191 ) .
the longitudinal susceptibility describes the @xmath192 , singlet to triplet excitations of the paramagnetic state while the transverse describes @xmath193 spin excitations .
spin rotational invariance dictates that these excitations should be degenerate .
one obtains for these correlation functions @xmath194 and @xmath195    for the non - interacting gutzwiller quasiparticles spin - rotational invariance is thus guaranteed from the relation @xmath196    this identity is preserved within the ga+rpa when we compute the interacting susceptibilities from the rpa series @xmath197 since from eqs .
( [ eq : vzz ] , [ eq : vpm ] ) the interaction kernels are related by @xmath198 .
clearly the energies of particle hole excitations described by these equations are degenerate .
in particular they lead to the same stoner criteria in both channels .
the expressions involve the following derivatives of @xmath199 in the longitudinal channel @xmath200 @xmath201 @xmath202 @xmath203 @xmath204 the derivatives in the transverse channel can be related to those in the longitudinal channel as follows , @xmath205 on the right we have given explicit expressions in terms of @xmath76 , the doping measured with respect to half - filling and the double occupancy @xmath95 in the paramagnetic state
.    99 a. macridin , m. jarrell , t. maier , and d.j .
scalapino , phys .
lett . * 99 * , 237001 ( 2007 ) .
markiewicz , s. sahrakorpi , and a. bansil , phys .
rev . b*76 * , 174514 ( 2007 ) .
f. ronning , k.m .
shen , n.p .
armitage , a. damascelli , d.h .
lu , z .- x .
shen , l.l .
miller , and c. kim , phys .
b*71 * , 094518 ( 2005 ) ; j. graf , g .- h .
gweon , k. mcelroy , s.y .
zhou , c. jozwiak , e. rotenberg , a. bill , t. sasagawa , h. eisaki , s. uchida , h. takagi , d .- h .
lee , and a. lanzara , phys .
lett . * 98 * , 067004 ( 2007 ) .
markiewicz , s. sahrakorpi , m. lindroos , hsin lin , and a. bansil phys .
rev . b*72 * , 054519 ( 2005 ) .
aligia , phys .
b*39 * , 6700 ( 1989 ) ; q. si , y. zha , k. levin , and j.p .
lu , phys .
rev . b*47 * , 9055 ( 1993 ) ; p. bnard , l. chen , and a .- m.s .
tremblay , phys .
rev . b*47 * , 15217 ( 1993 ) ; m. lavagna and g. stemmann , phys .
rev . b*49 * , 4235 ( 1994 ) ; f. onufrieva and p. pfeuty , phys .
b*61 * , 799 ( 2000 ) ; m.r .
norman , phys .
* , 14751 ( 2000 ) .
for an example of this extended commensurate nesting , see frames j ,  k of fig .
rice and g.k .
scott , phys .
lett . * 35 * , 120 ( 1975 ) .
r. hlubina and t.m .
rice , phys .
rev . b*51 * , 9253 ( 1995 ) .
andersen , o. jepsen , a.i .
liechtenstein , and i.i .
mazin , phys .
rev . b*49 * , 4145 ( 1994 ) .
a. diciolo , j. lorenzana , m. grilli , and g. seibold , phys .
b * 79 * , 085101 ( 2009 ) .
f. gnther , g. seibold , and j. lorenzana , phys .
98 * , 176404 ( 2007 ) . some previous determinations of stripe incommensurability vs doping , but for a limited range of @xmath1-values , are : j. lorenzana and g. seibold , phys
lett . * 89 * , 136401 ( 2002 ) ; g. seibold and j. lorenzana , phys .
rev . b*69 * , 134513 ( 2004 ) .
f. gebhard , phys .
b*41 * , 9452 ( 1990 )
. g. kotliar and a.e .
ruckenstein , phys .
lett . * 57 * , 1362 ( 1986 ) .
d. vollhardt , rev .
56 * , 99 ( 1984 ) .
y.m . vilk and a .- m.s .
tremblay , j. phys .
i * 7 * , 1309 ( 1997 ) .
markiewicz and a. bansil , phys .
rev . b*75 * , 020508(r ) ( 2007 ) .
dolgov , d.a .
kirzhnits , and e.g. maksimov , rev .
* 53 * , 81 ( 1981 ) .
p. abbamonte , a. rusydi , s. smadici , g.d .
gu , g.a . sawatzky , and d.l .
feng , nature phys . * 1 * , 155 ( 2005 ) .
s. smadici , p. abbamonte , m. taguchi , y. kohsaka , t. sasagawa , m. azuma , m. takano , and h. takagi , phys .
b*75 * , 075104 ( 2007 ) .
macdougall , a.a .
aczel , j.p .
carlo , t. ito , j. rodriguez , p.l .
russo , y.j .
uemura , s. wakimoto , g.m .
luke , phys .
* 101 * , 017001 ( 2008 ) .
shen , f. ronning , d.h .
lu , f. baumberger , n.j.c .
ingle , w.s .
lee , w. meevasana , y.kohsaka , m. azuma , m. takano , h. takagi , and z .- x .
shen , science * 307 * , 901 ( 2005 ) ; w.d .
wise , m.c .
boyer , k. chatterjee , t. kondo , t. takeuchi , h. ikuta , y. wang , and e.w .
hudson , nature physics * 4 * , 696 ( 2008 ) .
markiewicz and m.t .
vaughn , phys .
rev . b*57 * , r14052 ( 1998 ) .
storey , j.l .
tallon , and g.v.m .
williams , phys .
b * 76 * , 174522 ( 2007 ) .
palczewski , t. kondo , r. khasanov , n.n .
kolesnikov , a.v .
timonina , e. rotenberg , t. ohta , a. bendounan , y. sassa , a.v .
fedorov , s. pailhs , a.f .
santander - syro , j. chang , m. shi , j. mesot , h.m .
fretwell , and a. kaminski , phys .
b*78 * , 054523 ( 2008 ) .
a. kopp , a. ghosal , and s. chakravarty .
usa * 104 * , 6123 ( 2007 ) .
r. hlubina , s. sorella , and f. guinea , phys .
lett . * 78 * , 1343 ( 1997 ) ; r. hlubina , phys . rev . b*59
* , 9600 ( 1999 ) ; d. vollhardt , n. blmer , k. held , m. kollar , j. schlipf , and m. ulmke , z. phys .
b*103 * , 283 ( 1997 ) ; b. valenzuela , m.a.h .
vozmediano , and f. guinea , phys .
b*62 * , 11312 ( 2000 ) .
igoshev , a.a .
katanin , and v.yu .
irkhin , j.e.t.p .
, to be published , ( extended version : arxiv:0709.3219 ) .
m. hashimoto , t. yoshida , h. yagi , m. takizawa , a. fujimori , m. kubota , k. ono , k. tanaka , d.h .
lu , z .- x .
shen , s. ono , and yoichi ando , phys .
b*77 * , 094516 ( 2008 ) .
we thank tanmoy das for providing us with these values , based on lda .
norman , phys .
rev . b*75 * , 184514 ( 2007 ) .
zabolotnyy , a.a .
kordyuk , d.s .
inosov , d.v .
evtushinsky , r. schuster , b. buechner , n. wizent , g. behr , sunseng pyon , h. takagi , r. follath , s.v .
borisenko , europhys .
lett . * 86 * , 47005 ( 2009 ) .
markiewicz , phys .
* , 174518 ( 2004 ) .
markiewicz , j. lorenzana , and g. seibold , unpublished .
s. andergassen , s. caprara , c. di castro , and m. grilli , phys .
lett . * 87 * , 056401 ( 2001 ) .
j. lorenzana and g. seibold , low .
physics * 32 * , 320 ( 2006 ) .
g. seibold , f. becca , p. rubin , and j. lorenzana , phys . rev .
b * 69 * , 155113 ( 2004 ) .


wz sge - type dwarf novae are an unusual subtype of su uma - type dwarf novae ( @xcite ; @xcite ; @xcite ; @xcite ; @xcite ) . wz sge - type dwarf have properties common to su uma - type dwarf novae ( @xcite ; @xcite ; @xcite ) in that they show superoutbursts and superhumps .
the most unusual properties of wz sge - type dwarf novae are characterized by the extremely low frequency of outbursts ( usually once in years to decades ) , the extremely large outburst amplitude ( exceeding 6 mag ) and long ( in some cases reaching @xmath0100 d before the quiescence is reached ) duration , and the lack or extreme deficiency of normal ( short ) outbursts , which , in ordinary su uma - type dwarf novae , occur more frequently than superoutbursts ( cf .
@xcite ) .
the origin of such unique characters has been long sought , both from the observational and theoretical sides .
two basic ideas have been historically proposed , just like a descendant from historical discussions on the mechanism of dwarf nova outbursts .
one is the mass - transfer instability model , which assumes the outer atmosphere of the extreme low - mass secondary in wz sge - type dwarf novae become unstable to the irradiation from the outbursting accretion disk , thus giving rise to an enhanced mass - transfer from the secondary , which leads to a long - lasting superoutburst ( the best historical example of this application being @xcite ) .
the other is an application of the disk - instability model , which , with an assumption of very low quiescent viscosity of the accretion disk , can produce the main characteristics of wz sge - type dwarf novae , namely the long outburst intervals and the absence of normal outbursts ( @xcite ; @xcite ) .
the required quiescent viscosity ( @xmath1 ) to reproduce such long outburst intervals has been shown to be extremely small ( @xmath2 : @xcite ; @xmath3 : @xcite ) .
these values of the required quiescent viscosity has usually been considered too small for a disk with a developed magneto - rotational instability @xcite , although this problem is recently becoming more positively resolved by considering a suppression of this instability in a cold disk with a low electric conductivity @xcite .
the apparent difficulty with the low quiescent viscosity has historically invoked modifications of the disk - instability model exemplified by the truncated inner disk @xcite .
this model supposed a cold steady state - solution during quiescence , and a slight variation in mass - transfer rate was supposed to be the cause of infrequent outbursts . @xcite
proposed another possibility of the magnetically truncated disk .
although these truncated disks were successful in reproducing one of the characteristics of wz sge - type dwarf novae ( the long recurrence time ) , these models with a high @xmath1 ( @xmath0 0.01 ) have been shown by @xcite to meet difficulties in reproducing long - lasting superoutbursts .
although there have been attempts ( e.g. @xcite ) to reproduce long outbursts by considering an enhanced mass - transfer during superoutbursts , the resultant light curve is very different from observation .
furthermore , the supposed observational evidence for an enhanced mass - transfer is excluded by later careful analysis @xcite .
in addition to these outburst properties , wz sge - type dwarf novae are shown to have other properties common to all members , but are absent in ordinary su uma - type dwarf novae .
one is the presence of _ early superhumps _ during the earliest stage of superoutbursts of wz sge - type dwarf novae ( cf .
the early superhumps have a stable period almost identical with the orbital period [ during the best observed 2001 superoutburst of wz sge , the period of the early superhumps was reported to be very slightly shorter than the orbital period @xcite ] .
this feature was historically first detected and described in wz sge itself ( @xcite ; @xcite ) , although the origin of this variation was not properly discussed or identified at this time .
the common existence of early superhumps among wz sge - type dwarf novae has been progressively confirmed : al com ( @xcite ; @xcite ; @xcite ; @xcite ; @xcite ) , hv vir ( @xcite ; @xcite ) , rz leo ( @xcite ) , and the main focus of this paper , eg cnc ( @xcite ; @xcite ) .
the other is occasional  dips " ( transient fading episodes lasting for one to several days ) during the later stage of superoutburst , which was first clearly described in the 1978 outburst of wz sge ( e.g. @xcite ) , and was repeated in the 1995 outburst of al com ( @xcite ; @xcite ; @xcite ; @xcite ) .
@xcite repoeted that the stage after the dip phenomenon shares the characteristics common to superoutbursts of ordinary su uma - type dwarf novae : the presence of a likely precursor and the subsequent growth of usual superhumps . from these findings ,
the 1977 - 1978 outburst of wz sge and the 1995 outburst of al com were interpreted as
double superoutbursts " @xcite.http://www.kusastro.kyoto - u.ac.jp / vsnet / mail / alert1000/msg00195.html@xmath4 . ]    an argument naturally rises why such a sequence of superoutburst can be only realized in wz sge - type dwarf novae . a naive explanation may be a consequence of enhanced - mass transfer , but is rather difficult to reconcile considering the low proper angular momentum of the  fresh " transferred matter , accumulation of which will more easily produce normal outbursts than a second superoutburst .
@xcite , @xcite , and @xcite instead suggested a possibility of the cooling and heating waves originating before the termination of the entire superoutburst .
@xcite presented observational evidence in al com that the second superoutburst was triggered by a short outburst , which , from the absence of superhumps , was attributed to a normal outburst as seen in ordinary su uma - type dwarf novae .
independent observations during the second ( super)outburst are critically needed to confirm the universality of the claimed observational features , and to explain the unique outburst pattern of wz sge - type dwarf novae .
on 1996 november 30 , an alert ( vsnet - alert 599 ) through the vsnet alert system @xcite notifying a bright outburst of the suspected wz sge - type star , http://www.kusastro.kyoto - u.ac.jp / vsnet / mail / vsnet - alert/ + msg00599.html@xmath4 and @xmath5http://www.kusastro.kyoto - u.ac.jp / vsnet / dne / egcnc.html@xmath4 . ]
eg cnc @xcite announced the opening of a new story .
eg cnc was originally discovered by @xcite from his photographic archive on the occasion of an outburst in 1977 .
subsequent searches for other photographic materials between 1977 and 1983 , and searches for further outbursts by huruhata himself and a number of amateur astronomers have yielded a null result until 1995 .
@xcite correctly identified the quiescent counterpart on palomar sky survey prints , and suggested from its marked blue color and the apparent low frequency of outbursts that this object is similar to the dwarf nova wz sge .    upon the alert by schmeer , we started time - resolved ccd photometry using telescopes at osaka kyoiku university ( oku ) and ouda station @xcite , kyoto university .
the initial results proving the existence of early superhumps ( @xmath6 d ) and the subsequent development of usual superhumps ( @xmath7 d ) were already discussed in our earlier paper ( @xcite ) .
we discuss mainly on its peculiar late - stage phenomenon , particularly focusing on repetitive rebrightenings .
the methods of observation and reduction were the same as in @xcite .
the observations at osaka kyoiku university were carried out using the 50-cm telescope at the cassegrain focus ( focal length @xmath8 ) equipped with an astromed eev  88200 ccd camera ( @xmath9 pixels ) and a @xmath10-band filter .
the observations at ouda station were carried out using the 60-cm telescope at the cassegrain focus ( focal length @xmath11 ) equipped with a thomson th  7882 ccd camera ( @xmath12 pixels ) .
we used an interference filter designed to reproduce the johnson @xmath10 band .
the @xmath13 on - chip binning mode was used .
the frames obtained at osaka kyoiku university were reduced using the iraf .
the frames obtained at ouda station were reduced with a microcomputer - based automatic aperture and psf photometry packages developed by one of the authors ( tk ) .
we used gsc 1948.1452 ( @xmath10 = 12.82 ) as the local comparison star , whose constancy during the observation was confirmed to 0.01 mag using the check star gsc 1948.1193 ( @xmath10 = 13.86 ) .
heliocentric correction were applied before the following analysis .
table [ tab : log ] lists the log of ccd photometry .
ccccccc & start  end ( ut ) & @xmath14 & exp(s ) & mean @xmath10 mag & obs + 1996 & december & 2.625  2.791 & 195 & 30 - 90 & 12.38 & oku + & & 3.625  3.708 & 66 & 60 - 80 & 12.51 & oku + & & 7.679  7.811 & 208 & 45 & 12.42 & ouda + & & 8.715  8.870 & 176 & 60 & 12.48 & ouda + & & 9.629  9.755 & 195 & 45 & 12.54 & ouda + & & 11.635  11.870 & 345 & 45 & 12.76 & ouda + & & 13.746  13.877 & 166 & 60 & 13.18 & ouda + & & 14.683  14.882 & 322 & 45 & 14.32 & ouda + & & 20.666  20.860 & 112 & 120 & 16.13 & ouda + & & 21.550  21.871 & 478 & 45 - 60 & 13.29 & ouda + & & 22.770  22.824 & 49 & 80 - 90 & 14.10 & ouda + & & 23.555  23.856 & 264 & 90 & 15.39 & ouda + & & 24.555  24.870 & 277 & 90 & 15.92 & ouda + & & 25.586  25.697 & 102 & 60 & 15.59 & ouda + & & 26.700  26.883 & 140 & 45 - 60 & 16.23 & ouda + & & 27.619  27.883 & 322 & 60 & 16.22 & ouda + & & 30.600  30.871 & 320 & 45 - 90 & 15.84 & ouda + & & 31.519  31.745 & 195 & 90 & 16.30 & ouda + 1997 & january & 2.543  2.753 & 144 & 90 & 16.25 & ouda + & & 3.862  3.863 & 2 & 120 & 16.35 & ouda + & & 4.473  4.888 & 721 & 40 - 50 & 13.37 & ouda + & & 6.616  6.730 & 72 & 90 & 15.98 & ouda + & & 10.725  10.795 & 82 & 60 & 13.92 & ouda + & & 11.640  11.804 & 97 & 90 & 15.29 & ouda + & & 12.614  12.755 & 78 & 90 & 16.32 & ouda + & & 13.605  13.751 & 122 & 90 & 16.23 & ouda + & & 15.647  15.811 & 137 & 90 & 16.27 & ouda + & & 19.756  19.745 & 75 & 90 & 16.35 & ouda + & & 20.501  20.843 & 284 & 90 & 16.22 & ouda + & february & 1.664  1.669 & 5 & 90 & 17.17 & ouda + & & 3.433  3.439 & 5 & 120 & 17.38 & ouda + & & 5.443  5.454 & 6 & 180 & 17.29 & ouda + & & 27.512  27.697 & 87 & 90 - 120 & 17.82 & ouda + & & 28.490  28.617 & 61 & 150 & 17.78 & ouda + & march & 1.619  1.727 & 57 & 150 & 17.83 & ouda + & & 3.589  3.726 & 63 & 150 & 17.95 & ouda + & & 4.634  4.727 & 49 & 150 & 17.94 & ouda + & & 5.522  5.680 & 82 & 150 & 18.02 & ouda + & & 7.631  7.639 & 4 & 150 & 17.92 & ouda + & & 8.493  8.675 & 96 & 150 & 18.01 & ouda + & & 11.486  11.616 & 58 & 180 & 18.09 & ouda + & & 28.431  28.516 & 39 & 150 & 17.94 & ouda + & & 30.503  30.600 & 50 & 150 & 18.25 & ouda + & & 31.469  31.604 & 61 & 180 & 18.26 & ouda + & april & 1.435  1.602 & 58 & 240 & 18.15 & ouda + & & 13.455  13.463 & 5 & 150 & 17.97 & ouda + & & 14.467  14.472 & 4 & 150 & 17.73 & ouda + & & 24.442  24.450 & 5 & 170 & 18.45 & ouda + & & 25.459  25.463 & 3 & 150 & 18.55 & ouda +   +      ( 160mm,80mm)fig1.eps    figure [ fig : vis ] represents the overall course of the outburst from ccd and visual observations reported to the vsnet , including the available @xmath10-band ccd observations , averaged to one to a few points per night .
the object was first detected in outburst on 1996 november 30.917 ut at @xmath1512 @xcite .
the last available negative observation was done on november 22.685 ut , giving an upper limit of 13.3 ( t. watanabe , private communication ) .
the early rise was unfortunately not observed due to the interference by the moon .
the outburst thus started some time between november 22 and 30 .
the existence of early superhumps @xcite , which are known to rapidly decay within several days of the maximum light of a wz sge - type superoutburst , and the recorded maximum brightness ( @xmath16 = 11.8 ) close to that of the 1977 outburst , suggest that the detection by schmeer was made within a few days of the start of the outburst .    the main superoutburst lasted until december 15 ( hjd 2450433 ) , followed by a precipitous decline .
the object then faded close to or below @xmath10 = 16 . on december 20 (
hjd 2450438 ) , the first hint of rebrightening was recorded in the ouda data . on the next night , the object reached the peak brightness of @xmath10 = 13.29 , and soon started fading rapidly .
this rebrightening was first considered to be of the rather common phenomenon among short orbital - period su uma - type dwarf novae ( e.g. v1028 cyg : @xcite , see also @xcite ) , not necessarily confined to wz sge - type dwarf novae .
the unique property of eg cnc , however , emerged after this .
the star remained about 3 mag brighter than its quiescence , and again showed a rebrightening 7 days after its first one .
similar rebrightenings repetitively occurred six times in total ( counting those exceeding @xmath10 = 14 at maximum ) , separated by five to eleven days .
table [ tab : reblist ] summarizes the observed post - superoutburst rebrightenings .
although there was a suggestion of two rebrightenings after the 1994 superoutburst of another wz sge - type candidate , uz boo @xcite , , such a large number of repetitive brightenings had been never observed in any dwarf novae until this outburst of eg cnc .
ccc date of maximum & magnitude & @xmath17 ( d ) + 50439 & 13.3 & @xmath18 + 50446 & 13.4 & 7 + 50453 & 13.3 & 7 + 50458 & 13.2 & 5 + 50465 & 13.3 & 7 + 50475 & 13.2 & 10 +   +   +   +    after its sixth rebrightening on 1997 january 26 , eg cnc then entered a stage of slow , steady decline , keeping the rate of slow fading between rebrightenings .
this stage lasted at least until our last observation on 1997 april 25 ( 146 d after schmeer s detection of the outburst ) , when eg cnc had a magnitude of @xmath10 = 18.6 , which was still 0.6 mag above the quiescence ( cf .
eg cnc thus reproduced the long - fading tail of wz sge observed in 1978 ( cf .
@xcite ) and al com in 1995 ( @xcite ) , accompanied by numerous short rebrightenings instead of a superoutburst - looking second outburst in wz sge and al com .
figure [ fig : wzcomp ] shows the comparison of outburst light curves of three best - observed wz sge - type dwarf novae in the recent years .
as already described in @xcite , doubly humped early superhumps were observed during the earliest two nights of our observations . since the raw light curve was not presented in @xcite , we present the light curve in figure [ fig : earlylc ] . then  textbook " superhumps ( often referred to as _ common superhumps _ or _ ordinary superhump _ when discussing on wz sge - type stars ) emerged ( e.g. figure [ fig : shlc ] ) .
( 88mm,120mm)fig3.eps    ( 88mm,60mm)fig4.eps    ( 88mm,120mm)fig5.eps    ( 88mm,120mm)fig6.eps    figures [ fig : earlysh ] and [ fig : sh ] show the comparison of the profiles of these two kinds of superhumps between the recent three well - observed wz sge - type dwarf novae ( the data of al com from @xcite , hv vir from @xcite ) .
the periods ( 0.05877 d and 0.06038 d ) used in folding the data of eg cnc are our revised estimates , using the pdm ( phase dispersion minimization : @xcite ) technique to the datasets for 1996 december 23 ( the stage showing early superhumps ) and 1996 december 713 ( common superhumps ) , after subtracting the trend of linear decline in each dataset .
the new best periods were determined using a minimum determination algorithm developed by one of the authors ( tk ) .
the 1-@xmath19 errors for these periods were determined to be 0.00017 and 0.00001 d , respectively , using fernie s application to lafler - kinman class of statistics @xcite .
statistical f - tests yielded significance levels of the signals 87% ( @xmath14 = 261 ) and @xmath2099.9% ( @xmath14 = 805 ) , respectively .
as is evident from these figures , all of these three wz sge - type dwarf novae showed doubly peaked early superhumps , sometimes with a rather complex structure ( cf .
@xcite ) .
the smallness of the amplitude in eg cnc probably suggests that this early modulation had already started to decay at the time of observations , or that this system has a low orbital inclination ( cf .
@xcite ) . the existence of these early superhumps , whose periods are considered to be very close to the orbital periods ( wz sge : @xcite ; @xcite ; @xcite ; @xcite ) , al com : @xcite ) , is now regarded as an independent criterion for discriminating wz sge - type dwarf novae @xcite .
the common superhumps appearing later are very similar between different objects , and have almost identical profiles to those of ordinary su uma - type dwarf novae ( cf .
@xcite ; @xcite ) .
we then analyzed short - term variation during the post - superoutburst rebrightening stage .
we first divided the observations during these stages into the following three representative stages : ( 1 ) _ between rebrightenings _ ( faint state between repetitive rebrightenings ; data for 1996 december 26 - 31 and 1997 january 2 were used ) , ( 2 ) _ rebrightening maximum _ ( around peak brightness of rebrightenings ; the 1997 january 4 data were used ; note that the unique period identification is improbable because of the short duration of the data segment ) , and ( 3 ) _ fading tail _ ( 1997 february 27 and after ) .
we then removed the long - term trend of each dataset by subtracting a linear fit , and applied the pdm technique to the datasets .
the resultant theta diagrams , together with those of the main superoutburst , are shown in figure [ fig : pdm ] ( the weaker signals arounds 0.0569 and 0.0643 d are one - day aliases , which are not true signals ) .
( 88mm,120mm)fig7.eps    although the result of the second set is less conclusive in terms of period , the common presence of signals ( minima in theta ) around 0.0603 d suggest that the period of usual ( later ) superhumps during the main superoutburst persisted throughout the rebrightening phase and the fading tail .
statistical f - tests yielded significance levels of the signals 59% ( @xmath14 = 1041 ) , 87% ( @xmath14 = 716 ) , and 84% ( @xmath14 = 494 ) for the above three stages , respectively .
( 88mm,120mm)fig8.eps    figure [ fig : fold ] shows the phase - averaged light curves of individual data sets .
the hump profiles , though amplitudes are very different , resemble each other and that of usual superhumps .
this observation of strong superhump modulation after the main superoutburst safely preclude the possibility of the cooling white dwarf as a main source of the fading tail of wz sge - type dwarf novae , as suggested by @xcite .
( 88mm,60mm)fig9.eps    we further examined the evolution of these persistent superhumps with respect to the rebrightening stages .
we used the data between 1996 december 20 and 23 , which best covered the rebrightening ( the second rebrightening ) .
the resultant phase - averaged light curves are shown in figure [ fig : minifold ] , which again show the stability of the 0.0603-d persistent superhumps , despite that these persistent superhumps become weaker when the system brightens .
figure [ fig : miniamp ] shows the time - variation of the pulsed ( superhump ) flux and the total flux of the object .
both fluxes are shown in magnitude scale corresponding to apparent @xmath10 magnitudes ( i.e. 10% flux modulations at mean magnitude @xmath10 = 14.0 corresponds to a pulsed flux corresponding to @xmath10 = 16.5 ) .
despite the @xmath0 3 mag magnitude brightening of the system during the rebrightening , the pulsed ( superhump ) flux remained rather constant within a factor of two .
( 88mm,60mm)fig10.eps    cccccccc @xmath21 & hjd@xmath222450000 & source & @xmath23 & @xmath21 & hjd@xmath222450000 & source & @xmath23 + -39 & 422.882 & 2 & 0.007 & 128 & 432.959 & 2 & -0.008 + -38 & 422.943 & 2 & 0.008 & 129 & 433.020 & 2 & -0.007 + -23 & 423.853 & 2 & 0.011 & 159.5 & 434.862 & 2 & -0.008 + -7 & 424.809 & 2 & 0.000 & 209.5 & 437.890 & 2 & -0.002 + -6 & 424.872 & 2 & 0.003 & 214.5 & 438.188 & 1 & -0.006 + -5 & 424.935 & 2 & 0.005 & 215.5 & 438.252 & 1 & -0.003 + -4 & 424.993 & 2 & 0.003 & 216.5 & 438.311 & 1 & -0.004 + -3 & 425.055 & 2 & 0.004 & 223.5 & 438.728 & 2 & -0.010 + 0 & 425.236 & 1 & 0.004 & 230.5 & 439.158 & 1 & -0.003 + 1 & 425.299 & 1 & 0.007 & 258.5 & 440.856 & 2 & 0.003 + 10 & 425.839 & 2 & 0.003 & 259.5 & 440.913 & 2 & -0.000 + 11 & 425.900 & 2 & 0.003 & 260.5 & 440.974 & 2 & 0.000 + 12 & 425.961 & 2 & 0.004 & 261.5 & 441.039 & 2 & 0.005 + 13 & 426.022 & 2 & 0.005 & 262.5 & 441.091 & 1 & -0.004 + 17 & 426.263 & 1 & 0.004 & 263.5 & 441.163 & 1 & 0.008 + 18 & 426.324 & 1 & 0.004 & 264.5 & 441.214 & 1 & -0.002 + 25 & 426.743 & 2 & 0.000 & 265.5 & 441.288 & 1 & 0.012 + 26 & 426.805 & 2 & 0.002 & 291.5 & 442.855 & 2 & 0.008 + 27 & 426.861 & 2 & -0.002 & 307.5 & 443.820 & 2 & 0.006 + 28 & 426.923 & 2 & -0.001 & 387.5 & 448.654 & 2 & 0.005 + 32 & 427.165 & 1 & -0.001 & 487.5 & 454.700 & 2 & 0.008 + 33 & 427.225 & 1 & -0.001 & 503.5 & 455.675 & 2 & 0.017 + 41 & 427.707 & 2 & -0.002 & 521.5 & 456.760 & 2 & 0.014 + 42 & 427.768 & 2 & -0.002 & 538.5 & 457.787 & 2 & 0.014 + 43 & 427.828 & 2 & -0.002 & 586.5 & 460.690 & 2 & 0.016 + 44 & 427.892 & 2 & 0.001 & 668.5 & 465.638 & 2 & 0.009 + 45 & 427.950 & 2 & -0.001 & 684.5 & 466.599 & 2 & 0.003 + 65 & 429.162 & 1 & 0.002 & 720.5 & 468.783 & 2 & 0.011 + 67 & 429.282 & 1 & 0.001 & 738.5 & 469.858 & 2 & -0.002 + 68 & 429.343 & 1 & 0.002 & 752.5 & 470.704 & 2 & -0.002 + 93 & 430.846 & 2 & -0.006 & 853.5 & 476.812 & 2 & 0.003 + 94 & 430.909 & 2 & -0.003 & 868.5 & 477.716 & 2 & 0.001 + 95 & 430.970 & 2 & -0.003 & 885.5 & 478.740 & 2 & -0.003 + 96 & 431.031 & 2 & -0.002 & 934.5 & 481.700 & 2 & -0.004 + 100 & 431.272 & 1 & -0.003 & 1064.5 & 489.555 & 2 & -0.005 + 101 & 431.330 & 1 & -0.005 & 1083.5 & 490.699 & 2 & -0.009 + 104 & 431.510 & 2 & -0.007 & 1148.5 & 494.627 & 2 & -0.009 + 106 & 431.630 & 2 & -0.007 & 1513.5 & 516.678 & 2 & -0.015 + 107 & 431.691 & 2 & -0.007 & 1544.5 & 518.561 & 2 & -0.005 + 108 & 431.753 & 2 & -0.005 & 1613.5 & 522.717 & 2 & -0.019 + 110 & 431.870 & 2 & -0.009 & 1762.5 & 531.733 & 2 & -0.007 + 111 & 431.932 & 2 & -0.008 & 1794.5 & 533.669 & 2 & -0.005 + 112 & 431.992 & 2 &
-0.008 & 1811.5 & 534.705 & 2 & 0.004 + 116 & 432.233 & 1 & -0.009 & 1845.5 & 536.755 & 2 & -0.001 + 117 & 432.295 & 1 & -0.007 & 1853.5 & 537.252 & 2 & 0.013 + 118 & 432.355 & 1 & -0.008 & 1861.5 & 537.727 & 2 & 0.004 + 125 & 432.780 & 2 & -0.006 & 1877.5 & 538.701 & 2 & 0.011 + 126 & 432.842 & 2 & -0.004 & 1886.5 & 539.240 & 2 & 0.007 + 127 & 432.900 & 2 & -0.006 & & & & +   +   +   +    we then estimated the times of superhump maxima by eye examination of the light curve ( except at the peak of rebrightening ) and from the phase - averaged light curve ( around the rebrightening maximum ) .
the values , together with the literature values , are listed in table [ tab : shmax ] , with cycle counts ( @xmath21 ) starting with @xmath21 = 0 at hjd 2450425.236 .
the accuracy of maximum determination is typically 0.001 d during the main superoutburst , and 0.002 d during the rebrightening stage .    in order to see whether the observed maximum times can be expressed by a continuous period change of a coherent variation , or we need to introduce a discontinuous phase jump ( cf .
@xcite ) , we calculated residuals to linear fits to the observed times in both cases assuming the superhumps having no phase jump and a 0.5 phase jump .
the latter ( 0.5 phase jump , or phase reversal ) corresponds to the phenomenon what is called late superhumps ( @xcite ; @xcite ; @xcite ; @xcite ) .
( 88mm,60mm)fig11.eps    ( 88mm,120mm)fig12.eps    in the former case , the superhump period was confirmed to show a discontinuous change and larger @xmath23 values ( figure [ fig : oc1 ] ) , while the latter case showed a more reasonable , smooth period change when we assume @xmath24159 maxima are late superhumps .
this identification looks reasonable that @xmath21 = 159 corresponds to the termination of the main superoutburst , at which a 0.5 phase jump is known to be most frequently observed ( e.g. @xcite ; @xcite ) . in table
[ tab : shmax ] , we provide @xmath21 numbers assuming the 0.5 phase jump for @xmath24159 .
a linear regression to the observed superhump times ( table [ tab : shmax ] ) yielded the following ephemeris :    @xmath25    the @xmath23 s given in table [ tab : shmax ] were calculated against this linear regression ( see figure [ fig : oc2 ] ) .
the superhump period initially increased during the main superoutburst and the first rebrightening , and the period started to decrease following the main superoutburst . a similar sequence of the superhump period change is frequently observed in other wz sge - type dwarf novae or infrequently outbursting su uma - type dwarf novae ( e.g. @xcite ; @xcite ; @xcite ) .
a quadratic fit to the observed @xmath23 s for @xmath26210 yielded a definitely positive period derivative of @xmath27 .    to summarize , the 0.0603-d superhump signal persisted throughout the rebrightening stage and the fading tail
rebrightenings do not seem to affect the profile , superhump phase and pulsed flux .
these findings make a clear contrast to what are observed during the initial stages of an ordinary su uma - type superoutburst @xcite .
the repetitive brightenings less likely represent the start of another superoutburst in spite of accompanying superhumps .
instead , they bear more characteristics of normal outbursts . the persistent superhumps should be rather considered to be the consequence of a persistent eccentricity in the accretion disk , and rebrightening must have occurred independent of these superhumps .
as briefly discussed in section [ sec : intro ] , no confirmed outburst was observed since the 1977 ( discovery ) outburst until 1995 . in order to estimate the possibility of missed outbursts ( owing to the observational gaps ) we applied monte - carlo simulations on actual observations by the vsolj members .
the total number of observations was 238 , distributing between jd 2445344 and 2450933 , spanning @xmath015 yr .
the results were that 52% of simulated superoutbursts ( the light curve was assumed to identical with the present one ) and 86% of simulated normal outbursts ( maximum magnitude of 13.3 and decay rate of 1 mag d@xmath28 assumed ) were missed .
these probabilities can be considered as upper limits , since the vsolj database contains no observations by huruhata himself , and the object had been equally intensively monitored by a number of skilled visual observers , including p. schmeer .
though it is still premature to conclude on the previous occurrence of outbursts , the 1-@xmath19 lower limit of recurrence times of superoutbursts ( @xmath29 ) and normal outbursts ( @xmath30 ) can be set 4 and 1 yr , respectively , assuming the poisson statistic of the outburst occurrence .
the low number of outbursts has been also confirmed by the absence of a confirmed outburst in the vsnet observations after the 19961997 one .
these lower limits of recurrence times are located at the longest end of the distribution of @xmath29 among all su uma - type dwarf novae ( e.g. @xcite ; @xcite ) .      according to @xcite ,
the 1977 outburst started between 1977 november 9.8 and 12.8 ut . the brightest measured magnitude was 11.9 on 1977 november 12.8 .
since the photometric sequence of the comparison stars by huruhata closely matches the present scale ( m. huruhata , private communication ) , we can safely compare huruhata s values with our modern @xmath10-magnitudes .
the object declined to 12.4 mag on november 22.8 ut , 10.9 d after the initial detection .
the duration of the outburst indicates that the 1977 outburst was also a superoutburst . after an observational gap of 11 d , @xcite
further positively detected the variable on three films taken in two separate occasions in early december , giving estimates @xmath014.0 mag 2122 d and 25 d after the initial detection .
considering that the duration of the main superoutburst in 1996 did not exceed 23 d , even taking the observational gap into account , some of these december positive detections may be interpreted as a hint of rebrightenings , or a part of a second ( super)outburst . the main difference from the 19961997 outburst is the faintness ( 14.0 mag ) during this stage .
it is , however , clear no evidence was found regarding repeated rebrightenings as in 19961997 , and the 1977 outburst seems to more resemble the 1913 and 1946 outbursts of wz sge .
it would be important to note the same star at times exhibit different patterns of activity .      based on our own analysis
, we propose the period of early superhumps ( @xmath6 d ) to be the candidate orbital period .
@xcite alternatively proposed a different ( 0.05997 d ) period , which is extremely close to the superhump period .
we suspect that this period identification suffered from ambiguity from two reasons .
the spectroscopic period by @xcite was _ not _ directly from radial velocity variations but from a periodic variation of the asymmetry of the profile .
we must note that this method of period determination could have suffered from the remaining eccentricity in the disk , as would be naturally expected from the long persistence of the superhumps .
the photometric method was even more ambiguous ; our own analysis of the early - stage variations gave a different result @xcite , which was also favored by quiescent photometry @xcite .
we absolutely need an independent determination from _ true _ radial velocity observation before concluding on the true identification of the orbital period . in this paper
, we adopt our own identification ( @xmath6 d ) in the following discussion .
lcccc & eg cnc & al com & hv vir & wz sge + year of outbursts & 1977 , 1996 & 1892 , 1941 , 1961 & 1929 , 1939 , 1970 , & 1913 , 1946 , 1978 , + & & 1965 , 1974 , 1975 , & 1981 , 1992 & 2001 + & & 1995 , 2001 & & + maximum magnitude & 11.9 & 12.0 & 11.5 & 7.0 + minimum magnitude & 19.1 & 20.5 & 19.1 & 15.5 + outburst amplitude ( mag ) & 7.2 & 8.5 & 7.6 & 8.5 + supercycle ( yr ) & 19 ? & 1020 & 10 : & 2333 + normal outbursts ?
& no ? & yes & yes ?
& no + superoutburst duration ( d ) & @xmath20100 & @xmath2070 & @xmath2050 & 130 + long fading tail & yes & probably yes & yes & yes + type(s ) of rebrightening & repeated short & superoutburst - like & probably yes & long ( similar + & outbursts & & ( short ? ) & to al com ) + early superhumps & yes & yes & yes & yes + period of early superhumps ( d ) & 0.05877 & 0.05666 & 0.05708 & 0.05667 + period of common superhumps ( d ) & 0.06038 & 0.05722 & 0.05820 & 0.05721 + superhump excess ( % ) & 2.7 & 1.0 & 2.0 & 1.0 + quiescent humps & ? & double & ? & double+eclipses +   +   +   +   +    as @xcite pointed out , the optical behavior of the present outburst bears all the characteristics common to those of known wz sge - type dwarf novae .
table [ tab : wzcomp ] summarizes the comparison of the observed properties between these objects .
the late - stage behavior of superoutbursts has a wide range of diversity between objects , and even between different outbursts of the same object , in spite of the marked similarity of the general light curve and the time evolution of superhumps during the early stage of outbursts .
phenomenologically , the superoutbursts of wz sge - type dwarf novae seem to make a continuous sequence ranging from superoutburst without noticeable rebrightenings ( wz sge in 1913 and 1946 ) , ones accompanied by a single rebrightening ( eg cnc in 1977 ; possibly hv vir in 1992 , cf .
@xcite ) , double to multiple rebrightenings ( eg cnc in 19961997 ; uz boo in 1994 : @xcite ) , and to two distinct outbursts separated by a  dip " ( wz sge in 19781979 , 2001 ; al com in 1961 , 1995 and 2001 ) , some of which may be double superoutbursts .    a similar idea of a continuum of types of outbursts of large - amplitude , infrequently outbursting dwarf novae ( tremendous outburst amplitude dwaef novae or toads : @xcite ) was discussed in @xcite , who mainly treated the diversity of duration and brightness of superoutbursts in these systems .
an ideal model of wz sge - type dwarf novae should reproduce this wide and continuous variety of outburst activity by changing some input parameters .      in observationally interpreting the rebrightening phenomenon , discrimination of _ inside - out_-type ( heating wave originating from the inner accretion disk ) , or _
outside - in_-type outbursts ( @xcite ; the terminology corresponding to type b and type a outbursts in @xcite ) , is potentially important : different parameters of the disk instability model predict different origins of the instability ( e.g. @xcite ; @xcite ) , which may provide a potential observational diagnostic to the underlying mechanisms .
( 88mm,120mm)fig13.eps    a close examination of the light curve of each rebrightening revealed the presence of five rebrightenings which showed a rapid rise and one which showed a much slower rise .
figure [ fig : minicomp ] demonstrates enlarged light curves of each rebrightening .
the first five rebrightenings showed rather common outburst behavior : an abrupt rise and a 3-mag fade in 3 d. the time needed for the rise is well confined to less than 0.4 d from the observations of the fourth and fifth rebrightenings .
the sixth rebrightening is different , accompanied with a slow rise lasting more than a day before reaching maximum .
( 88mm,60mm)fig14.eps    ( 88mm,60mm)fig15.eps    a similar phenomenon , but with a lesser amplitude , was also observed in time - resolved ccd photometry of the first rebrightening ( figures [ fig : reb1 ] and [ fig : miniprof ] ) .
a fast - rising type rebrightening occurred at around hjd 2450439 , while a more slowly rising brightening followed at around hjd 2450443 ( not counted as one of six major rebrightenings because of its faintness ) shows much a slower rise and faint peak brightness ( the faintness of the peak was also confirmed from the subsequent ccd observations reported to the vsnet ) . 1 d ) of the phenomenon , as well as the continued slow rising for 7 hr , was much longer than the superhump period . ]
although whether or not this small brightening can be treated in the same context of other ( major ) rebrightenings is an open question , the slow rising stage of this small brightening was more analogous to the sixth rebrightening , possibly implying the common underlying mechanism .
the existence of two types ( fast - rising and slow - rising ) of rebrightenings may represent the coexistence of _ inside - out_-type and _ outside - in_-type outbursts , more usual ones being the faster _ outside - in_-type .
the sixth rebrightening may be related to the former ,  slow " rebrightening .
since the disk instability theory predicts the cycle length of normal outbursts , for a certain range of parameters , is roughly inversely proportional to the mass - transfer rate ( e.g. @xcite ) , it would be natural to attribute the increased outburst activity after the main superoutburst to an enhanced mass - transfer somehow caused by the superoutburst .
since the measured short outburst intervals ( 510 d ) are one of the shortest even among frequently outbursting su uma - type dwarf novae ( cf .
@xcite ; @xcite ) , this interpretation would require ( assuming no change in other parameters before and after the superoutburst ) a mass - transfer rate something between ordinary su uma - type dwarf novae and er uma stars ( @xcite ; @xcite ; @xcite ; @xcite ; @xcite ) .
the er uma stars are peculiar su uma - type stars which are proposed to have mass - transfer rates several times higher than in ordinary su uma - type dwarf novae ( @xcite ; @xcite ) .
although the initial pattern and amplitudes of rebrightenings of eg cnc closely resembled those of er uma stars , the enhanced mass - transfer and the enhanced hot spot brightness would not solely by themselves explain the later behavior , since the brightness of the system was shown to decrease gradually ( see figure [ fig : vis ] ) , even after the sixth rebrightening .
this smooth decrease in brightness should require _ smoothly decreasing _ frequency of rebrightenings , which is contrary to the observation if the system brightness actually reflects the brightness of the hot spot or the mass - transfer rate .
@xcite instead proposed a model to reproduce these rebrightenings , by assuming increased quiescent viscosity after the main superoutburst .
this model , by adjusting the time dependence of the quiescent viscosity and its radial dependence through the disk , could at least reproduce the observed rebrightenings without an assumption of an enhanced mass - transfer .
@xcite originally suggested that such increased quiescent viscosity is somehow related to the turbulent motion in the disk resulting from the tidal instability .
they also proposed the extremely low mass - transfer rate in wz sge - type dwarf novae would be responsible for the slow circularization of the eccentric disk , leading to the slow decrease in the quiescent viscosity after the outburst and multiple rebrightenings never seen in usual dwarf novae .
the underlying physical mechanism for the decay of the disk viscosity has been recently proposed by @xcite by considering a decay of mhd turbulence under the condition of the low magnetic reynolds numbers in the cold accretion disk @xcite .
a potential difficulty of this model @xcite is that the disk after the main superoutburst would not expand enough to trigger a second superoutburst as reported in al com @xcite .
[ note that the outburst of al com after the  dip " was not composed of recurring small rebrightenings , as observed in the 2001 superoutburst of wz sge . ]
@xcite and @xcite proposed the idea that some material can be piled up just behind the cooling wave , causing it to be reflected as a heating wave , which they considered to cause an additional brightening .    regarding the reason why such a condition is predominantly met in wz sge - type dwarf novae , or toads
, @xcite suggested that the stored disk material survives being depleted by accretion , owing to the low occurrence of normal outbursts in these systems .
however , the degree of depletion of the disk mass during one supercycle is rather expected to affect the initial condition of superoutburst , but not necessarily that of the late - stage condition of superoutburst .
as @xcite has shown , the late stage of wz sge - type superoutbursts is expected to follow the same time - evolution as those of usual dwarf novae .
the initial excessive matter is effectively accreted on a viscous time - scale ; it is not still convincingly clear how the low occurrence of normal outbursts ( or its precondition ) naturally explains the unique late - stage phenomenon in wz sge - type superoutbursts .
though the discussion was more limited to sxts , @xcite could reproduce a similar pattern of rebrightenings by modifying the radial dependence of the @xmath31 viscosity .
this model , however , would suffer , as in the model proposed by @xcite , from the depletion of high - angular momentum matter necessary to invoke an immediate second superoutburst .
we alternately consider a possibility that the matter in the outer accretion disk may be left from being accreted during the main superoutburst , by way of the premature quenching of the hot state in the outer disk .
we propose the weaker turbulence and resultant heating caused by the weaker tidal instability in the outer disk than in ordinary su uma - type dwarf novae can be responsible for this possibility .
this weaker tidal torque can be reasonably achieved under conditions of extreme mass ratios in wz sge - type dwarf novae , which can hold the 3:1 resonance radius well inside the roche lobe of the primary , or the maximum dimension of the expanded disk .
this possibility was originally proposed in 1997 by @xcite and @xcite , which was further explored by @xcite and @xcite
. the unusual presence of na d absorption line @xcite could be attributed to the cool reservoir in the disk required by this interpretation .
superoutbursts in ordinary su uma - type dwarf novae are believed to be triggered when the disk radius , upon the ignition of the outburst , first reaches the radius of the 3:1 resonance during the normally outbursting cycle @xcite .
the eccentricity wave @xcite caused by the tidal instability should naturally start in the outermost region of the disk , and should propagate inward , accounting for the general decrease in the superhump period during superoutbursts of ordinary su uma - type dwarf novae @xcite .    in wz sge - type dwarf novae ,
however , the absence of normal outbursts may lead to an accumulation of a large amount of matter , leading to more violent expansion of the disk during superoutburst enabling the disk expanding beyond the 3:1 resonance .
the above picture then predicts the eccentricity wave can propagate outward .
the presence of positive period derivatives of superhump periods in some wz sge - type dwarf novae and infrequently outbursting , large - amplitude su uma - type dwarf novae ( e.g. al com : @xcite ; v1028 cyg : @xcite ; sw uma : @xcite ; @xcite ; hv vir : @xcite ; see also @xcite for a recent summary of period derivatives in su uma - type dwarf novae ) would strengthen the idea of the relative difference in the location of generation and propagation of tidal instability between wz sge - type and ordinary su uma - type dwarf novae .
the extremely low mass - ratios in wz sge - type dwarf novae could play the key role in generating tidal instability well inside the entire disk radius .
we obtained extensive ccd photometry during the 19961997 superoutburst of eg cnc .
this outburst showed an unprecedented series of six major rebrightenings spaced by 510 d. during this rebrightening stage and the following slow fading tail , photometric modulations with a period of 0.0603 d persisted , which is equal to the superhump period observed in the latter course of the main superoutburst .
we suggest that the superhumps observed during the rebrightening stage and the fading tail are a  remnant " of usual superhumps , and are not a result of rebrightenings . by comparison with the 1977 outburst of this object and outbursts of other wz sge - type dwarf novae
, we propose an activity sequence of wz sge - type superoutbursts , in which the current outburst of eg cnc is placed between a single - rebrightening event and distinct outbursts separated by a dip .
we propose a new potential explanation of post - superoutburst behavior of the wz sge - type dwarf novae , by considering the effect of extreme mass ratios in these binaries on the generation and propagation of the tidal instability .
we are grateful to many amateur observers for supplying their vital visual and ccd estimates via the vsnet alert network , and especially to patrick schmeer for his detection and early notification of the long - awaited outburst .
we are grateful to the vsolj staffs and members providing earlier observations .
we are also grateful to jonathan kemp for providing up - to - date information , and people involved in extensive on - line discussions in _ vsnet - chat _ ; part of this log can be electronically accessed at http://www.kusastro.kyoto-u.ac.jp/vsnet/dne/egcnc.html and its related links .
we express our sincere gratitude to yoji osaki and shin mineshige for valuable discussions made at the annual meeting of the astronomical society of japan .
part of this work was supported by a research fellowship of the japan society for the promotion of science for young scientists ( dn ) .