since the initial work of kobayashi and maskawa , the standard model ( sm ) of electroweak interactions is known to be able to accommodate the experimentally observed cp
violation through a unique phase , @xmath1 , in the cabibbo  kobayashi
maskawa mixing matrix ( ckm ) . however , the available experimental information , namely @xmath2 and @xmath3 , is not enough to establish this phase as the only source of cp  violation .
most of the extensions of the sm include new observable phases that may significantly modify the pattern of cp violation .
supersymmetry is , without a doubt , one of the most popular extensions of the sm .
indeed , in the minimal supersymmetric extension of the sm ( mssm ) , there are additional phases which can cause deviations from the predictions of the sm .
after all possible rephasings of the parameters and fields , there remain at least two new physical phases in the mssm lagrangian .
these phases can be chosen to be the phases of the higgsino dirac mass parameter ( @xmath4 $ ] ) and the trilinear sfermion coupling to the higgs , ( @xmath5 $ ] ) @xcite . in fact , in the so  called constrained minimal supersymmetric standard model ( cmssm ) , with strict universality at the grand unification scale , these are the only new phases present .    it was soon realized that for most of the cmssm parameter space , the experimental bounds on the electric dipole moments of the electron and neutron constrained @xmath6 to be at most @xmath7 .
consequently these new supersymmetric phases have been taken to vanish exactly in most studies of cmssm .
however , in the last few years , the possibility of having non - zero susy phases has again attracted a great deal of attention .
several new mechanisms have been proposed to suppress edms below the experimental bounds while allowing susy phases @xmath8 .
methods of suppressing the edms consist of cancellation of various susy contributions among themselves @xcite , non - universality of the soft breaking parameters at the unification scale @xcite and approximately degenerate heavy sfermions for the first two generations @xcite . in a recent work @xcite
, we showed that , in a model with heavy sfermions of the first two generations and in the large @xmath9 regime , @xmath2 and @xmath10 could receive very sizeable contributions from these new susy phases .
similar studies @xcite including a larger set of experimental constraints have reported the impossibility of such large supersymmetric contributions . and @xmath11 @xcite . ]    in this work , we are going to complete our previous analysis with the inclusion of all the relevant constraints in a cmssm scenario .
in doing so we adopt a different perspective .
we will assume from the very beginning that both supersymmetric phases are @xmath8 , ignoring for the moment edm bounds .
in these conditions , and taking into account other cp  conserving constraints , we will analyze the effects on the low energy cp
violation observables , especially @xmath2 and @xmath10 .
it should be noticed that the model used in @xcite can be easily obtained as a limit of the cmssm by decoupling the first two generations of squarks and neglecting the intergeneration mixing in the sfermion mass matrices .
hence , our results in the more general cmssm will include this model as a limiting case .
in the next section we study the new sources of flavor mixing present at the electroweak scale in any supersymmetric model . in section 3
we are going to analyze neutral meson mixing , i.e. @xmath12@xmath13 and @xmath14@xmath15 mixings , with large supersymmetric phases .
section 4 will be devoted to the study of the branching ratio of the decay @xmath16 . in section 5
we will show the impact of the measured @xmath16 branching ratio on the supersymmetric contributions to @xmath2 and @xmath10 .
section 6 will present our conclusions and in appendices a and b we collect , respectively , the formulas for the integration of relevant rge s and the different loop functions appearing in the text .
the issue of flavor changing neutral current ( fcnc ) interactions in the cmssm has been widely investigated in the literature . for the completeness of the discussion
, we briefly recall those properties which will be relevant for our analysis .
the cmssm is the simplest supersymmetric structure we can build from the sm particle content .
this model is completely defined once we specify the soft - supersymmetry breaking terms .
these are taken to be strictly universal at some ultra - high energy scale , which we take to be the grand unification scale ( @xmath17 ) .
@xmath18    that is , there is a common mass for all the scalars , @xmath19 , a single gaugino mass , @xmath20 , and all the trilinear soft - breaking terms are directly proportional to the corresponding yukawa couplings in the superpotential with a proportionality constant @xmath21 .
now , with the use of the renormalization group equations ( rge ) of the mssm , as explained in appendix a , we can obtain the whole supersymmetric spectrum at the electroweak scale .
all the supersymmetric masses and mixings are then a function of @xmath19 , @xmath20 , @xmath22 , @xmath23 , @xmath24 and @xmath9 .
we require radiative symmetry breaking to fix @xmath25 and @xmath26 @xcite with tree
level higgs potential . and
@xmath27 are switched off .
for this reason we prefer to give the relevant rge s with non  vanishing susy phases in appendix a. ]    it is important to notice that , even in a model with universal soft breaking terms at some high energy scale as the cmssm , some off  diagonality in the squark mass matrices appears at the electroweak scale . working on the basis where the squarks are rotated parallel to the quarks , the so  called super ckm basis ( sckm ) , the squark mass matrix is not flavor diagonal at @xmath28 .
this is due to the fact that at @xmath17 there exist two non - trivial flavor structures , namely the two yukawa matrices for the up and down quarks , which are not simultaneously diagonalizable .
this implies that through rge evolution some flavor mixing leaks into the sfermion mass matrices . in a general supersymmetric model
, the presence of new flavor structures in the soft breaking terms would generate large flavor mixing in the sfermion mass matrices .
however , in the cmssm , the two yukawa matrices are the only source of flavor change . as always in the sckm basis ,
any off - diagonal entry in the sfermion mass matrices at @xmath28 will be necessarily proportional to a product of yukawa couplings .
the rge s for the soft breaking terms are sets of linear equations , and thus , to match the correct quirality of the coupling , yukawa couplings or tri - linear soft terms must enter the rge in pairs , as we can see in eqs .
( a.1-a.3 ) in appendix a.    in fact , in the up ( down ) squark mass matrix the up ( down ) yukawas will also be diagonalized and so will mainly contribute to diagonal entries while off  diagonal entries will be due to the down ( up ) yukawa matrix .
this means , for instance , that in this model the off - diagonality in the @xmath29 matrix will roughly be @xmath30 . with @xmath31 a proportionality factor that typically is , @xmath32 as expected from the loop factor and the running from @xmath17 to @xmath28
nevertheless , we have to keep in mind that this is simply a typical estimate and the final value of @xmath31 can suffer a sizeable variation depending on many other factors not present in this simple estimate .    on the other hand , this has clear implications on the @xmath9 dependence of these off  diagonal entries of the sfermion mass matrices . in the basis where the down yukawa matrix is diagonal , we can write the up and down yukawas as , @xmath33 with @xmath34 the cabibbo - kobayashi - maskawa mixing matrix and @xmath35 the diagonalized mass matrices for the quarks
. we can see in this equation that , for @xmath36 , the up yukawa matrix will maintain similar values when going to large @xmath9 .
hence , the off  diagonal entries in the down squarks mass matrix will be roughly stable with @xmath9 . in the up squark mass matrix
we have the opposite situation and the @xmath9 dependence is very strong . in this case
the off  diagonal entries depend on the down yukawa matrix that grows linearly with @xmath9 for large @xmath9 .
this means that we can expect the flavor change in the up and down squark mass matrix to be similar when @xmath37 . while for @xmath38 the flavor change in the up mass matrix will be approximately @xmath39 times smaller ( see appendix a for details ) .
these points also apply to the left
right sub
matrices where again flavor changing entries will be due to the opposite isospin yukawa matrix .
in fact , this left
right sfermion mixing only appears after electroweak symmetry breaking .
the expression for these matrices in the sckm basis is , @xmath40 ) .
this is always true for all the generation changing entries that are produced by the @xmath41 matrices .
however , in the down mass matrix , this suppression can be partially compensated by a large value of @xmath9 in the diagonal terms proportional to @xmath42 .
these are all well  known facts in the different studies of fcnc processes in the framework of the cmssm @xcite and imply that flavor mixing is still dominantly given by the usual ckm mixing matrix in w - boson , charged higgs and chargino vertices .
in this work , we are especially interested in cp violating observables
. then we must also consider the presence of observable phases in the sfermion mass matrices . in the following
we will take the ckm matrix exactly real to isolate pure effects of the new supersymmetric phases @xcite .
the sfermion mass matrices contain several physical phases that give rise to cp violation phenomena .
in particular , before rge evolution , these phases ( @xmath23 , @xmath24 ) are confined to the left - right part of the sfermion mass matrix while both the left  left , @xmath43 , and right  right , @xmath44 , matrices are real and diagonal .
however this is not true anymore at @xmath28 : @xmath23 leaks into the off - diagonal elements of these hermitian matrices through rge evolution . from the explicit rge in the mssm , eq.([mqrge ] ) , it is clear that this phase only enters the @xmath45 evolution through the combinations @xmath46 or @xmath47 . at @xmath17
these matrices have a common phase , and so the combination @xmath48 is exactly real .
so , to the extent that the @xmath41 matrices keep a uniform phase during rge evolution , no phase will leak into the @xmath49 matrices .
however , we can easily see from eqs.(a.2 ) and ( a.3 ) that this is not the case , and different elements of the @xmath41 matrices are renormalized differently . in this equation
, we can see that only the terms involving two yukawa and one @xmath41 matrix can produce a mismatch in the phases .
moreover , these terms will only be important when there are no small yukawas involved .
then , we can expect a mismatch only on the off
diagonal elements involving the third generation . keeping this in mind , the general form of the @xmath49 matrix at @xmath28 in terms of the initial conditions is , @xmath50 where the coefficients @xmath51 are @xmath52 matrices with real numerical entries . in this expression
we can see that the presence of imaginary parts will be linked to the non - symmetric part of the @xmath53 matrices . as is clear from the mass matrices in appendix a ( eqs .
( a.5-a.7 ) and ( a.12-a.14 ) ) , these non  symmetric parts of @xmath54 are always more that three orders of magnitude smaller than the corresponding symmetric parts .
this means that , in the sckm basis , the imaginary parts of any mass insertion are present only in one part per @xmath55 , and are always associated with @xmath56 mi , eqs .
( a.8-a.11 ) and ( a.15-a.18 ) .
a very similar situation was also found by bertolini and vissani , in the cmssm with vanishing susy phases for the leakage of @xmath1 @xcite .
so , we conclude that in the processes we will consider , we can take both @xmath57 and @xmath58 as real to a very good approximation .    in the following
we will analyze the new effects of this model on indirect cp violation in @xmath59 and @xmath60 systems .
in doing so , we will use both the exact vertex mixing method and the mass insertion ( mi ) approximation @xcite .
notice that the mi approximation is extremely good in the case of the cmssm where all the off - diagonal entries are sufficiently small .
the size of these off  diagonal entries directly gives in the mi approximation the amount of flavor changing induced by the sfermion mass matrices .
a possible exception may arise in the stop and sbottom sectors that , in any case , could be diagonalized to ensure the validity of the mi approximation @xcite . as we will see in the next sections ,
this is frequently useful to understand the exact results obtained in the vertex mixing method .
in the sm neutral meson mixing arises at one loop through the well  known @xmath61box . however , in the cmssm , there are new contributions to @xmath62 processes coming from boxes mediated by supersymmetric particles .
these are charged higgs boxes ( @xmath63 ) , chargino boxes ( @xmath64 ) and gluino - neutralino boxes ( @xmath65 , @xmath66 ) .
the amount of the indirect cp violation in the neutral meson @xmath67 system is measured by the well  known @xmath68 parameter : @xmath69 where @xmath70 is the @xmath67@xmath71 mass splitting .
@xmath68 depends on the matrix elements of the @xmath62 hamiltonian , @xmath72 , which can be decomposed as @xmath73 where the relevant four  fermion operators are given by @xmath74 with @xmath75 for the @xmath59 and @xmath60systems respectively and @xmath76 as color indices . in the cmssm , these are the only three operators present in the limit of vanishing @xmath77 .    at this point , we are going to divide our discussion into two parts .
we analyze separately the effective operator @xmath78 that preserves chirality along the fermionic line , and the operators @xmath79 and @xmath80 that change chirality along the fermionic line . as we will see below the flavor mixing in the sfermion mass matrix and the experimental constraints on both kinds of operators are very different .      in eq.([df=2 ] ) , @xmath78 is the only operator present that does not involve a chirality change in the fermionic line .
with respect to the associated sfermion , no chirality change in the sfermion propagator will be needed , and so , the suppression associated with left
right sfermion mixing can be avoided . in general , @xmath81 can be decomposed as follows @xmath82    the usual sm @xmath61box , where all the couplings are purely left
handed can only contribute to this effective operator .
however , with @xmath83 , @xmath84 does not contain any complex phase and hence can not contribute to the imaginary part in @xmath85 . in any case
it will always be , in the cmssm , the dominant contribution to @xmath70 .
similarly , the charged higgs contribution , @xmath86 depends on the same combination of ckm elements with no other cp  violating phase @xcite .
so it will not contribute to our cp violating observable .
gluino and neutralino contributions to @xmath87 are specifically supersymmetric . they involve the superpartners of quarks and gauge bosons . here , the source of flavor mixing is not directly the usual ckm matrix .
it is the presence of off  diagonal elements in the sfermion mass matrices , as discussed in the previous section . from the point of view of cp violation
, we will always need a complex wilson coefficient . in the sckm basis
all gluino vertices are flavor diagonal and real .
this means that in the mi approximation , we need a complex mass insertion in one of the sfermion lines .
as explained in the previous section , these mi are proportional to yukawa couplings and real up to 1 part in @xmath88 .
the complete expressions for the gluino contributions to @xmath62 processes in the mi approximation can be found in @xcite .
the bounds obtained there for the real and imaginary parts of the mass insertions required to saturate @xmath89 and @xmath2 are , @xmath90 where @xmath91 is an average squark mass .    in the cmssm , as we can see in appendix a , these mass insertions are much smaller . in particular , the fact that the bound on @xmath89 , the real part of the mi , is satisfied implies that the imaginary parts are at least two orders of magnitude below the required value to saturate @xmath2 .
hence , no sizeable contributions to @xmath2 from gluino boxes are possible .
the situation in @xmath14@xmath15 mixing is completely analogous : assuming that the minimum phase from the mixing observable in the b  factories is around @xmath92 radians , we would need an imaginary contribution not more than one order of magnitude below the real one . with the arguments given above , this is clearly out of reach for gluino boxes in the cmssm .
neutralino contributions are generally smaller than gluino due to smaller couplings with the same source of flavor mixing . in fact , although neutralino vertices in the sckm basis also involve the complex neutralino mixings , any imaginary part on this operator will only be due to a complex mass insertion .
this can be seen in the explicit expressions in @xcite where all neutralino mixings in this operator appear in pairs with its complex conjugate counterpart .
finally , the charginos also contribute to @xmath93 . in this case ,
flavor mixing comes explicitly from the ckm mixing matrix , although off
diagonality in the sfermion mass matrix introduces a small additional source of flavor mixing .
@xmath94 \end{aligned}\ ] ]    where @xmath95 , @xmath96 , and @xmath97 represent the coupling of chargino , @xmath98 , and squark @xmath99 to the left  handed down quark , @xmath100 .
finally @xmath101 is an intermediate up
quark index associated with the factorization of the ckm mixing matrix .
the expression for this coupling is then , @xmath102 where @xmath103 and @xmath104 are @xmath105 matrices such that the @xmath106 unitary matrix @xmath107 diagonalizes the up  squark mass matrix , @xmath108 .
@xmath109 is one of the matrices that diagonalize the chargino mass matrix through a bi  unitary transformation @xmath110 , with , @xmath111 from these equations it is clear that @xmath112 will in general be complex , as both @xmath24 and @xmath23 are present in the different mixing matrices .
the loop function @xmath113 is given in eq .
( b.1 ) of appendix b.    the main part of @xmath114 in eq([chwc ] ) will be given by pure ckm flavor mixing , neglecting the additional flavor mixing in the squark mass matrix @xcite .
in this case , @xmath115 and @xmath116 , we have , @xmath117 \end{aligned}\ ] ] but , taking into account that @xmath118 is symmetric under the exchange of any pair of arguments we have , @xmath119 and so @xmath120 , is exactly real @xcite .
this is not exactly true in the cmssm , where there is additional flavor change in the sfermion mass matrices . here ,
some imaginary parts appear in the @xmath114 in eq([chwc ] ) .
being associated to the size of intergenerational sfermion mixings , these imaginary parts will be maximal for large @xmath9 . in fig([imk40 ] ) we show in a scatter plot the size of imaginary and real parts of @xmath114 in the k system for a fixed value of @xmath121 .
the region of susy parameters explored in this and all of the following scatter plots is @xmath122 and @xmath123 . with these initial conditions we impose that all squarks are heavier than @xmath124 with the exception of the stops that , as the charginos , are only required to be above @xmath125 .
furthermore we impose the constraint from the @xmath16 decay .
notice that , as we will see later , this is a conservative attitude , in the sense that other constraints that we do not impose could only make our conclusions stronger . under these conditions , we can see here that , in the cmssm , this wilson coefficient is always real up to a part in @xmath126 .
fig([imb40 ] ) is the equivalent plot for the case of @xmath14@xmath15 mixing . here , imaginary parts are relatively larger but , in any case , out of reach for the foreseen b  factories .    taking this into account , from the point of view of experimental interest , we will always neglect imaginary parts in the wilson coefficient @xmath127 within the cmssm .
notice that this would not apply in a general model with non - universality at the gut scale @xcite , and each particular model should be considered separately .      from the point of view of flavor change and cp violation ,
operators @xmath79 and @xmath80 are different from @xmath78 .
these two operators always involve a change in the chirality of the external quarks , and consequently also a change of the chirality of the associated squarks or gauginos .
in particular , this implies the direct involvement of the supersymmetric phases . on the other hand ,
these operators are suppressed by the presence of down quark yukawa couplings , and so can only be relevant in the region of large @xmath9 @xcite .
we can write the different contributions to @xmath128 and @xmath129 as , @xmath130 in first place , charged higgs contributes only to @xmath128 , but parallel to the discussion for @xmath131 , the absence of phases prevents it from contributing to @xmath132 .
gluino and neutralino boxes contribute both to @xmath79 and @xmath80 .
however flavor change will be given in this case by an off  diagonal left
right mass insertion . in the cmssm
these mi are always proportional to the mass and are never enhanced by large @xmath9 values ( see eq.([lr ] ) ) of the right handed squark .
this implies that these left
right flavor transitions from gluino will always be smaller in the cmssm than the corresponding chargino contributions , where flavor change is directly given by the ckm matrix .
in fact , this is already well
known for the case of @xmath16 decay @xcite , which is completely equivalent from the point of view of flavor change .
hence , the most important contribution , especially for light stop and chargino , will be the chargino box . before the inclusion of qcd effects
, it contributes solely to the coefficient @xmath129 ,    @xmath133    where @xmath134 is the coupling of chargino , @xmath98 , and squark , @xmath99 , to the right  handed down quark @xmath100 , with , @xmath135 and @xmath136 is given in eq .
( b.2 ) . unlike the @xmath114 wilson coefficient , @xmath137 is complex even in the absence of intergenerational mixing in the sfermion mass matrices @xcite .
in fact , the presence of flavor violating entries in the up  squark mass matrix hardly modifies the results obtained in their absence @xcite .
so , in these conditions we have , @xmath138\\ f_s(\alpha , k,\gamma , l , i , j)= \frac{m_q^2}{2 m_w^2 \cos^2 \beta }   h^{(\alpha , k)i}g^{(\alpha , k)j * }   g^{(\gamma , l)i * } h^{(\gamma , l)j } y_2(z_{k } , z_{l } , s_i , s_j )   \nonumber\end{aligned}\ ] ] where we have used ckm unitarity and degeneracy of the first two generations of squarks . due to the differences between @xmath139 and @xmath140 couplings , this contribution is always complex in the presence of susy phases .
the most relevant feature of eqs.([chwcr ] ) and ( [ chwcr0 ] ) is the explicit presence of the external quark yukawa coupling squared , @xmath141 .
this is the reason why this contribution is usually neglected in the literature @xcite .
however , as we showed in @xcite , this contribution could be relevant in the large @xmath9 regime . for instance
, in @xmath14@xmath15 mixing we have @xmath142 that for @xmath143 is larger than 1 and so , it is not suppressed at all when compared with the @xmath114 wilson coefficient .
this means that this contribution can be very important in the large @xmath9 regime @xcite and could have observable effects in cp violation experiments in the new b  factories .
however , in our previous work @xcite , we did not include the additional constraints coming from @xmath16 decay . in the next sections we will analyze the relation of @xmath132 with this decay , and the constraints imposed by its experimental measure .
the decay @xmath16 has already been extensively studied in the context of the cmssm with vanishing susy phases @xcite . because the branching ratio is a cp conserving observable , the presence of new phases will not modify the main features found in @xcite concerning the relative importance of the different contributions .
however , in the presence of the new susy phases , these contributions will have different phases and will be observable through the interference . as we will see next , the experimental constraints will also have a large impact on the imaginary parts of the decay amplitudes .
this decay is described by the following @xmath144 effective hamiltonian @xmath145 where the relevant operators are given by @xmath146 here @xmath147 , and the wilson coefficients @xmath148 can be decomposed according with the particles in the loop , @xmath149    among these contributions , the @xmath61 penguin is exactly the same as in the sm and it does not depend on any supersymmetric parameters , it is simply a function of sm couplings and masses .
this contribution is @xcite , @xmath150 with @xmath151 and @xmath152 the charge of the up quarks .
similarly , in the charged higgs penguins all the variables are known with the exception of @xmath153 .
again this contribution is unchanged by the inclusion of the new susy phases , @xmath154 where @xmath155 .
this contribution gives a sizeable correction to the @xmath16 decay that constrains the mass of the charged higgs in two higgs doublet models or in the mssm with low @xmath9 .
however , in the case of moderate  large @xmath9 , chargino contributions may partially compensate this charged higgs contribution relaxing the constraints @xcite .
in addition to the @xmath156 and charged higgs contributions analyzed above , there are three specifically supersymmetric contributions mediated by gluino , neutralino and chargino . in gluino or neutralino diagrams
flavor change is due to the off  diagonality in the sdown mass matrix .
being left
right flavor off  diagonal transitions , they are suppressed by the mass of the @xmath157 quark . indeed , smallness of gluino and neutralino contributions has already been established in @xcite where it was shown that , in the cmssm , such contributions are roughly one order of magnitude smaller than the chargino contribution .    together with the @xmath158 and charged higgs , the most important supersymmetric contribution will be , especially in the large
moderate @xmath9 regime , the chargino contribution . in the @xmath61 and
charged higgs contributions , the necessary chirality flip for the dipole amplitude is always proportional to @xmath159 . however , in the chargino penguin the chirality flip can be made either through a chargino mass insertion in the loop or through an external leg mass insertion proportional to @xmath159 .
in fact , as pointed out in @xcite , this enhancement due to @xmath160 is partially compensated by the presence of the @xmath157 yukawa coupling .
nevertheless , this compensation is only effective for low values of @xmath161 . in terms of the chargino  quark
squark couplings used in the previous section , these contributions are , @xmath162 with the loop functions defined in appendix b. similarly to the situation for the wilson coefficient @xmath129 , we can , to a very good approximation , neglect the presence of intergenerational mixing in the up  squark mass matrix @xcite , then , @xmath163 where , once more , we use ckm unitarity and degeneracy of the first two generations of squarks .    the second term in @xmath164 in eq . ( [ char ] ) , which corresponds to the chargino mass insertion in the loop , is dominant in the large @xmath9 regime .
notice that both @xmath165 and @xmath166 are products of the squark and chargino mixing matrices that can be @xmath167 ( in the case of flavor  diagonal stop mixings ) .
then , for stop and chargino masses around the electroweak scale , this term has an extra enhancement of @xmath168 .
this means that , for large @xmath9 , we can approximate these wilson coefficients as , @xmath169
as we have seen in section iii , chargino contribution to the @xmath129 wilson coefficient , eq.([chwcr0 ] ) , is the main contribution to indirect cp violation of the new supersymmetric phases for large values of @xmath9 .
however , if we compare this wilson coefficient with the chargino contribution to the decay @xmath16 , eqs .
( [ chwcr0 ] ) and ( [ char ] ) , we can see that both chargino contributions are deeply related .
in fact , if we make a rough approximation and assume that the two different loop functions involved are of the same order , i.e. , @xmath170 we would obtain , @xmath171 of course , this can not be considered as a good approximation . as we can see from their explicit expressions in appendix b ,
the loop functions are clearly different .
anyway , they can be expected to give results of the same order of magnitude .
so , the order of magnitude of @xmath129 is determined by the allowed values of @xmath172 , as we will explicitly show below .    to reach this goal , we will follow @xcite , where they constrain in a model  independent way new physics contributions to the wilson coefficients involved in the @xmath16 decay . in terms of these wilson coefficients , the branching ratio , @xmath0 is , @xmath173 + 0.161 re[\xi_{8 } ] + 0.083 re[\xi_{7}\xi_{8}^ { * } ] \nonumber\end{aligned}\ ] ] where @xmath174 . the different coefficients appearing in eq.([constrain ] ) are the sm renormalization group evolved contributions , that must be recovered in the limit @xmath175 .
the numerical values are taken from @xcite .
we have not taken into account the errors associated with the choice of the scale and the restrictions on the photon energy that do not modify our conclusions .
now , using the experimental measure , @xmath176 , we can constrain the allowed values of the complex variables @xmath177 and @xmath178 .
in fact , we can already see from eq .
( [ constrain ] ) , that in the approximation @xmath179 this is simply the equation of an ellipse in the @xmath180$]@xmath181 $ ] plane . in the case of supersymmetry with large @xmath9 ,
the new physics contribution to @xmath177 and @xmath178 will be mainly due to the chargino . the allowed values of @xmath182 directly constrain then the chargino contributions to @xmath183 and indirectly constrain the values of @xmath184 .    in figure [ figc7 ] ,
we show a scatter plot of the allowed values of @xmath185 versus @xmath186 in the cmssm for a fixed value of @xmath9 with the constraints from eq.([constrain ] ) .
notice , that a relatively large value of @xmath9 , for example @xmath187 , is needed to compensate the @xmath61 and charged higgs contributions and cover the whole allowed area with positive and negative values .
however , the shape of the plot is clearly independent of @xmath9 , only the number of allowed points and its location in the allowed area depend on the value considered . in this figure
we take @xmath188 because only a large value could give rise to observable cp violation @xcite .
the values of @xmath172 and @xmath189 used here are the values obtained in the cmssm for a given set of initial conditions . although we do not use the approximation @xmath179 this does not modify the elliptic shape of the plot .
figure [ figc3 ] shows the allowed values for a re  scaled wilson coefficient @xmath190 corresponding to the same allowed points of the susy parameter space in figure[figc7 ] .
as we anticipated previously , the allowed values for @xmath191 are close to the square of the values of @xmath172 in figure [ figc7 ] slightly scaled by different values of the loop functions .
this is the proof of the importance of the @xmath16 constraint on the chargino contributions to indirect cp violation .
we can immediately translate this result to a constraint on the size of the chargino contributions to @xmath132 .
@xmath192\end{aligned}\ ] ] in this expression @xmath193 , @xmath194 and @xmath195 denote the mass , mass difference and decay constant of the neutral meson @xmath196 .
the coefficient @xmath197 @xcite includes the rge effects from @xmath28 to the meson mass scale , @xmath42 , and @xmath198 is the b
parameter associated with the matrix element of the @xmath80 operator @xcite .    then , for the @xmath59 system , using the experimentally measured value of @xmath199
we obtain , @xmath200   \approx 0.4 \times 10^{-7 } im[\bar{c}_3]\end{aligned}\ ] ] given the allowed values of @xmath191 in fig.[figc3 ] , this means that in the cmssm , even with large susy phases , chargino can not produce a sizeable contribution to @xmath2 .
we have seen in section iii that also gluino and neutralino give negligible contributions in the cmssm or in a model without off  diagonal soft  breaking terms at the gut scale .
hence indirect cp violation in the kaon system will be mainly given by the usual sm box and the presence of a cp violating phase in the ckm matrix , @xmath1 is still needed .
the case of @xmath14@xmath15 mixing has a particular interest due to the arrival of new data from the b  factories .
in fact , as explained at the end of section iii and in @xcite , in the large @xmath9 regime chargino contributions to indirect cp violation can be very important .
however , for any value of @xmath9 we must satisfy the bounds from the @xmath16 decay .
so , if we apply these constraints to the @xmath14@xmath15 mixing , @xmath201   \approx 0.5 \times 10^{-3 } im[\bar{c}_3]\end{aligned}\ ] ] where once again , with the allowed values of fig.[figc3 ] we get a very small contribution to cp violation in the mixing . we must take into account that the mixing  induced cp phase , @xmath202 , measurable in @xmath14 cp asymmetries , is related to @xmath10 by @xmath203 .
the expected sensitivities on the cp phases at the b factories are around @xmath204 radians , so this supersymmetric chargino contribution will be completely out of reach .
gluino and neutralino contributions to indirect cp violation can also be discarded in the cmssm .
once again we have to conclude that no new contributions to indirect cp violation from the new susy phases will be observable in @xmath14 cp asymmetries in the framework of the cmssm .
recently , cdf @xcite has provided preliminary indications that @xmath205 is in agreement with the sm predictions . clearly , from the above result , eq .
( [ epsb ] ) , it appears that the cmssm contribution is too small by itself to account for this result .
in this work , the effects of non - vanishing supersymmetric phases on indirect cp violation in @xmath59 and @xmath60 systems have been analyzed within the cmssm .
we have found that operators involving only left  handed external quarks are not sensitive to these new phases at an observable level .
this is due to the absence of intergenerational mixings beyond those originated from the ckm matrix . on the contrary
, operators involving both right and left
handed quarks are in general complex , even in the absence of @xmath1 , and could be relevant in the large @xmath9 regime .
however , we have shown that these contributions are deeply related with the @xmath0 decay .
so , taking into account the constraints coming from this decay these contributions also turn out to be too small to be measured experimentally .
although these conclusions are specific for indirect cp violation , they could also be implemented for chargino mediated direct cp violation in the decays .
again , in these decays the same chargino  quark
squark couplings are involved and we can also expect a big impact of the @xmath16 constrain .
in fact , the conclusions reached in this paper are far more general .
the correlation between @xmath16 and susy induced indirect cp violation exists in any supersymmetric model with sufficiently small intergenerational mixings in the sfermion mass matrices .
this would include specifically all the models without new flavor structures beyond the usual ckm matrix at the gut scale and simplified models as the one the authors used in @xcite .    in summary , concerning the simpler supersymmetric models , like cmssm , the constraints coming from @xmath0 decay are sufficient to rule out pure supersymmetric indirect cp violation in @xmath59 and @xmath60 systems , even in the absence of any electric dipole moment constraints .
this has very important consequences for the supergravity induced models where a cancellation between different supersymmetric contributions allows large supersymmetric phases while respecting edm bounds @xcite . in these models , even in the regions of parameter space where this cancellation occurs , no observable effect of the large susy phases will appear on indirect cp violation experiments
. however , as pointed out by baek and ko @xcite , these phases would still be observable in cp asymmetries in the @xmath16 decay .
all this means that the presence of large susy phases is not sufficient to produce observable effects at the low energy experiments .
in particular , new sources of flavor change beyond the usual ckm matrix are needed . and
so , any deviation from the sm expectations at indirect cp violation experiments due to supersymmetry should be taken as a sign of non  universality of the soft breaking terms . in this context one
recalls the recent studies on superstring compactifications with non
universal gaugino masses @xcite .
in this appendix we describe the new features of the integration of rge s in the cmssm with non - vanishing susy phases relevant to our analysis .
the complete matrix form of the rg equations can be found in @xcite . using their notation and conventions , with the only change of @xmath206
, we will mainly concentrate on the left  left scalar  quark mass matrix and the tri  linear soft breaking coupling evolution , @xmath207 \nonumber + \\ & \frac{1}{2 } [ \tilde{y}_d\tilde{y}_d^\dagger m_q^2 + m_q^2 \tilde{y}_d \tilde{y}_d^\dagger + 2 ( \tilde{y}_d m_d^2 \tilde{y}_d^\dagger + \bar{\mu}^2_1   \tilde{y}_d\tilde{y}_d^\dagger +   \tilde{a}_d \tilde{a}_d^\dagger ) ] \end{aligned}\ ] ]    @xmath208 \end{aligned}\ ] ]    @xmath209 \end{aligned}\ ] ]    except for the yukawa coupling matrices , the rge s of all other quantities are linear @xcite .
this means , in particular , that rge s of all soft masses , though coupled , can be solved as a linear combination of the gut  scale parameters @xmath210 , @xmath211 and @xmath212 at any scale below @xmath213
. however , one notices that the initial conditions on the trilinear couplings require the knowledge of the particular yukawa texture at the unification scale . to do this
, we numerically integrate the yukawa rge with a given value of @xmath9 and in terms of the fermion masses and the ckm matrix . specifying the gut
scale initial conditions in this way , it is straightforward to compute all soft masses at @xmath214 for arbitrary values of @xmath210 , @xmath211 and @xmath212 .
thanks to the linearity of the corresponding rge s the soft masses at @xmath214 admit the following expansion , @xmath215 where the coefficients @xmath101 and @xmath51 are @xmath52 matrices with real numerical entries .
one notices that the matrices @xmath216 would be completely real were it not for the non  symmetric terms in the matrix @xmath217 .
however , it will be seen from the specific examples that this matrix remains nearly symmetric , and thus , cp  violating entries @xmath216 are extremely suppressed .
moreover , one notices that @xmath218 carries , in general , large cp violating phases ; however , these terms are effective only for intragenerational @xmath219-type mixings .
hence , this particular observation shows the importance of chargino contributions for cp violation in fcnc processes , as explained in sec .
ii .    as mentioned before , due to the non  linearity of the rge s for yukawa matrices , it is not possible to give a fully analytic solution for the soft mass parameters .
nevertheless , once we fix @xmath220 , we can numerically integrate the yukawa rge .
therefore , below we give semi  analytic solutions of rge s for @xmath221 and @xmath222 to illustrate the small and large @xmath220 regimes .
fixing @xmath223 , we get for the relevant @xmath51 matrices in eq.([etas ] ) ,    @xmath224    where the vanishing off  diagonal entries in the last matrix mean values smaller than @xmath225 in absolute magnitude . among the matrices involved in eq.([etas ] ) , @xmath226 is always the largest one for similar values of @xmath212 and @xmath227 .
so , it sets the scale of the matrix element while @xmath228 is the only one that can produce an imaginary part .
hence , we do not specify the other @xmath51matrices , which are not important for our discussion .
once we obtain the @xmath229 matrix with the help of eq.([etas ] ) we can get the values of the @xmath230 and @xmath231 in the sckm basis , that give the size of flavor change in the squark mass matrices compared with the diagonal elements . for @xmath221 ,
those elements of the squark mass  squared matrix causing @xmath232 transitions between first and second , as well as second and third generations , are given by    @xmath233    now , we repeat the same quantities above for @xmath222 :    @xmath234    @xmath235    a comparison of the corresponding quantities in @xmath221 and @xmath222 cases reveals the sensitivity of the results on @xmath220 . as explained in sec .
ii , @xmath236 remains nearly unchanged while @xmath237 assumes an order of magnitude enhancement as @xmath220 varies from 2 to 40 .
this change in @xmath237 affects various quantities as dictated by the differential equations ( a.1  a.3 ) .
in this appendix we collect the different loop function used in the text . the functions @xmath238 and @xmath239 entering @xmath60@xmath240 and @xmath59@xmath241 mixings are given by @xmath242 and @xmath243\,.\end{aligned}\ ] ]      @xmath245 , \\
f_{r}^{7}(x , y ) & = &          \frac{1}{x }            \left[q_{u}f_{4}(y / x)+f_{3}(y / x)\right ] ,                        \\    f_{l}^{8}(x , y ) & = &          \frac{1}{x}f_{2}(y / x ) ,                        \\    f_{r}^{8}(x , y ) & = &          \frac{1}{x } f_{4}(y / x)\,.\end{aligned}\ ] ]    m. dugan , b. grinstein and l.j .
hall , _ nucl .
phys . _ * b 255 * , 413 ( 1985 ) ; + s. dimopoulos and s. thomas , _ nucl
* b 465 * , 23 ( 1996 ) . t. ibrahim and p. nath , _ phys .
rev . _ * d57 * , 478 ( 1998 ) ; + t. ibrahim and p. nath , _ phys .
rev . _ * d58 * , 111301 ( 1998 ) ; + m. brhlik , g.j .
good and g.l .
kane , _ phys .
* d59 * , 115004 ( 1999 ) hep - ph/9810457 .
abel and j.m .
rev . _ * d 55 * , 1623 ( 1997 ) ; + s.  khalil , t.  kobayashi and a.  masiero , _ phys .
rev .  _ * d60 * , 075003 ( 1999 ) , hep - ph/9903544 .
s. dimopoulos and g.f .
giudice , _ phys .
* b 357 * , 573 ( 1995 ) ; + a. cohen , d.b .
kaplan and a.e .
nelson , _ phys .
* b 388 * , 599 ( 1996 ) ; + a. pomarol and d. tommasini , _ nucl
* b 466 * , 3 ( 1996 ) .
demir , a.  masiero and o.  vives , phys .
* 82 * ( 1999 ) 2447 , hep - ph/9812337 .
s.  baek and p.  ko , korea inst .
sci report no .
kaist - th-99 - 1 , april 1999 hep - ph/9904283 ; + s.  baek and p.  ko , phys . rev
* 83 * ( 1999 ) 488 , hep - ph/9812229 .
s.  baek , j. h.  jang , p.  ko and j. h.  park , korea inst .
sci report no .
kaist-99 - 04 , july 1999 , hep - ph/9907572 . n.k .
falck , z. phys . *
c30 * ( 1986 ) 247 .
s. bertolini , f. borzumati , a. masiero , and g. ridolfi , nucl . phys .
* b353 * ( 1991 ) 591 .
p.  cho , m.  misiak and d.  wyler , phys . rev .
* d54 * ( 1996 ) 3329 , hep - ph/9601360 g. c. branco , f. cagarrinho , and f. kruger , phys . lett .
* b459 * ( 1999 ) 224 , hep - ph/9904379 .
s.  bertolini and f.  vissani , phys . lett .
* b324 * ( 1994 ) 164 , hep - ph/9311293 .
.  hall , v.a .
kostelecky and s.  raby , nucl . phys .
* b267 * ( 1986 ) 415 .
a. j.  buras , a.  romanino , l.  silvestrini , nucl .
phys . * b520 * ( 1998 ) 3 .
f.  gabbiani , e.  gabrielli , a.  masiero and l.  silvestrini , nucl . phys .
* b477 * ( 1996 ) 321 hep - ph/9604387 .
branco , g.c .
cho , y. kizukuri and n. oshimo , _ phys .
_ * b 337 * , 316 ( 1994 ) ; + g.c .
branco , g.c .
cho , y. kizukuri , and n. oshimo , nucl .
* b449 * ( 1995 ) 483 .
f. borzumati , z. phys .
* c63 * , 291 ( 1994 ) ; + s. bertolini , f. vissani , z. phys . *
c67 * , 513 ( 1995 ) ; + t. goto , y.y .
keum , t. nihei , y. okada and y. shimizu , phys . lett . *
b460 * , ( 1999 ) 333 , hep - ph/9812369 .
a. l. kagan and m. neubert , eur .
j. * c7 * , 5 ( 1999 ) .
m. ciuchini _
_ , _ jhep _ * 10 * ( 1998 ) 008 , hep - ph/9808328 ; + r. contino and i. scimemi , _ eur .
phys .  j.  _ * c10 * , 347 ( 1999 ) , hep - ph/9809437 .
t. affolder _
_ , fermilab report no .
fermlab - pub-99/225-e , september 1999 , hep - ex/9909003 .
m. brhlik , l. everett , g. l. kane , and j. lykken,_phys .
lett . _ * 83 * , 2124 ( 1999 ) , hep - ph/9905215 .


there are indications both from observations and theoretical studies that stars in many astrophysical environments may have the relative abundances of elements synthesized by nuclear @xmath4-capture reactions ( e.g.  c , o , ne , mg , si , ti , etc ) that are enhanced with respect to the solar ones . in brief , king ( 1994 ) found the mean [ o / fe ] ratio for a handful of halo stars to be @xmath8  dex .
similar results are found for globular cluster stars ( e.g.  [ @xmath9=0.18 $ ]  dex for 47 tuc , carney 1996 ) and galactic bulge k giants ( mcwilliam & rich 1994 , rich & mcwilliam 2000 ) . according to models of chemical evolution of the galaxy ,
bulge stars should be @xmath4-enhanced to reproduce the metallicity distribution ( matteucci & brocato 1990 , matteucci et al .
indications that the abundance ratios of elements in elliptical galaxies are not solar ( e.g. @xmath10 = 0.3 - 0.7 $ ]  dex ) come from worthey et al .
also synthetic metal line strengths ( weiss , peletier & matteucci 1995 ) from high - metallicity models show that solar - scaled models do not fit the observed mg2 ( 5180   ) and fe ( 5270 and 5335   ) line indices .    given these premises ,
it is of paramount importance to include the enhancement of @xmath4-elements in the population synthesis models designed to describe those galaxies . in the past
this has been done but using models that suffered from many limitations .
firstly , till recently there were no extended sets of opacities for @xmath4-enhanced mixtures covering the complete relevant range of stellar temperatures .
such opacity calculations are nowadays available .
secondly , it has often been assumed that the behaviour of @xmath4-enhanced isochrones can be easily mimicked by suitably scaling the metallicity  ( salaris et al .
1993 ) . however , the validity of such a re - scaling has only been confirmed for low metallicities , and it holds only for some particular cases like all @xmath4-elements being enhanced by the same factor .
weiss et al .  ( 1995 ) already pointed out the differences in the evolution due to @xmath4-enhancement at solar and super - solar metallicities .
recent work by salaris & weiss ( 1998 ) and vandenberg et al .
( 2000 ) also indicates that the simple re - scaling of the  content becomes risky at high metallicities .
moreover , we should keep in mind that there is no _ a priori _ reason , from the point of view of both nucleosynthetic and chemical evolution theories , to expect that all @xmath4-elements are enhanced by the same factor . from observations of halo and bulge stars ,
the degree of enhancement of different @xmath4-elements is , in general , not constant with differences amounting to @xmath11  dex .
therefore in this paper we intend to re - analyze the impact of the enhancement of @xmath4-elements on the properties of metal - rich stellar populations ( whose metallicity ranges from half - solar to three times - solar ) .
the analysis is made for a particular choice of abundance ratios , for which self - consistent opacities are available .
we start calculating extended sets of stellar evolutionary tracks ( sect .
[ sec_tracks ] ) , with total metal contents @xmath12 , 0.019 , 0.040 and 0.070 both for solar - scaled ( same abundance ratios as in the solar mix ) and @xmath4-enhanced patterns of abundances .
particular attention is paid to the description of the adopted initial element abundances and sources of opacity .
based on those stellar models we derive large grids of isochrones and single stellar populations for which we present magnitudes , colours and mass - to - light ratios in many pass - bands of the johnson and wfpc2 photometric systems ( sect .
[ sec_isochrones ] ) .
extensive tabulations of stellar models , isochrones and photometric properties are provided in electronic format for the purpose of general use .
finally , a short summary is given in sect .
[ sec_conclusion ] .
vandenberg et al .
( 2000 ) recently presented extensive sets of evolutionary tracks for mixtures , limited to the range of low - mass stars ( @xmath13 ) .
a constant enhancement of all @xmath4-elements has been assumed ( also in their opacity tables ) . in the present work , instead
, the choice of initial abundances has been guided by observational results .
it is worth remarking that , within the present errors , none of both choices seems to be strongly preferred .
we calculate four grids of stellar models with initial masses from 0.15 to 20 @xmath14 and initial chemical compositions @xmath0 $ ] , @xmath1 $ ] , @xmath2 $ ] , @xmath3 $ ] .
the metallicities and the helium - to - metal enrichment ratio @xmath15 are chosen in such a way to secure consistency with the girardi et al .
( 2000 ) models .
the helium - to - metal enrichment law is @xmath16 .    for the above values of @xmath17
$ ] , we compute both solar - scaled and @xmath4-enhanced sets of stellar tracks .
the evolutionary phases covered by the grids extend from the zero age main sequence ( zams ) up to either the start of the thermally pulsing asymptotic giant branch ( tp - agb ) phase or carbon ignition .
the input physics of the present stellar models is the same as in girardi et al .
( 2000 ) , apart from differences in ( i ) the opacities , and ( ii ) the rates of energy loss by plasma neutrinos . a short description of the updated ingredients and the way we deal with the enhancement of @xmath4-elements is presented in the sections below .
for a fixed total metallicity @xmath18 , the solar - scaled abundance ratios of metals are taken from grevesse & noels ( 1993 ) . in the @xmath4-enhanced case
, we adopt a mixture in which only the elements resulting from nuclear @xmath4-capture reactions ( o , ne , mg , si , s , ca , ti , ni ) are enhanced with respect to the solar abundance ratios .
this is the same mixture as used by salaris et al .
( 1997 ) and salaris & weiss ( 1998 ) to study the age of the oldest metal - poor and disk metal - rich globular clusters . for the sake of clarity , table [ elements ] lists the abundances for the solar - scaled and @xmath4-enhanced mixtures .
columns ( 2 ) and ( 4 ) show the abundance @xmath19 of the elements in logarithmic scale , @xmath20 , where @xmath21 is the abundance by number . columns ( 3 ) and ( 5 ) display the abundance by mass fraction @xmath22 , relative to the metallicity @xmath18 . columns ( 6 ) and ( 7 ) list the enhancement ratio ( relative either to fe or the total metal fraction ) in spectroscopic notation .
the enhancement factors come from the determinations of chemical abundances in metal - poor field stars by ryan et al .
( 1991 , cf .  also salaris & weiss 1998 ) . to derive the initial abundances of different isotopes we use the isotopic ratios ( by mass ) of anders & grevesse ( 1989 ) .
c|cc|cccc element & & + & @xmath19 & @xmath23 & @xmath19 & @xmath23 & @xmath24 $ ] & @xmath25 $ ] + ( 1 ) & ( 2 ) & ( 3 ) & ( 4 ) & ( 5 ) & ( 6 ) & ( 7 ) + o & 8.870 & 0.482273 & 9.370 & 0.672836 & 0.50 & @xmath260.1442 + ne & 8.080 & 0.098668 & 8.370 & 0.084869 & 0.29 & @xmath270.0658 + mg & 7.580 & 0.037573 & 7.980 & 0.041639 & 0.40 & @xmath260.0441 + si & 7.550 & 0.040520 & 7.850 & 0.035669 & 0.30 & @xmath270.0558 + s & 7.210 & 0.021142 & 7.540 & 0.019942 & 0.33 & @xmath270.0258 + ca & 6.360 & 0.003734 & 6.860 & 0.005209 & 0.50 & @xmath260.1441 + ti & 5.020 & 0.000211 & 5.650 & 0.000387 & 0.63 & @xmath270.1875 + ni & 6.250 & 0.004459 & 6.270 & 0.002056 & 0.02 & @xmath270.3371 + c & 8.550 & 0.173285 & 8.550 & 0.076451 & 0.00 & @xmath270.3553 + n & 7.970 & 0.053152 & 7.970 & 0.023450 & 0.00 & @xmath270.3553 + na & 6.330 & 0.001999 & 6.330 & 0.000882 & 0.00 & @xmath270.3553 + cr & 5.670 & 0.001005 & 5.670 & 0.000443 & 0.00 &
@xmath270.3557 + fe & 7.500 & 0.071794 & 7.500 & 0.031675 & 0.00 & @xmath270.3557 +      the opacities we have adopted are the same as in salaris , deglinnocenti & weiss ( 1997 ) .
they are available at : ` http://www-phys.llnl.gov/v_div/opal/existing.html ` .
details are given below .
for high temperatures ( @xmath28 ) we use the radiative opal opacities by rogers & iglesias ( 1995 ) in the case of solar - scaled mixtures , and by iglesias & rogers ( 1996 ) for the ones . at low temperatures these
are combined with the molecular opacities provided by alexander & ferguson ( 1994 ) for the solar mixtures , and by alexander ( priv .
communication ) for the ones . in both cases ,
low- and high - temperature opacity grids are merged by interpolation in the temperature range from @xmath29 to @xmath30 .
for very high temperatures , @xmath31 , the weiss et al .
( 1990 ) opacities are used .
these rely on the los alamos opacity database .
the radiative opacities for c , o mixtures , needed in the post - main sequence phases , are from iglesias & rogers ( 1993 ) . finally , conductive opacities of electron - degenerate matter are from itoh et al .
( 1983 ) .
it is very important to emphasize that the present models are based on _ complete and consistent opacities _ ,
i.e.   the opacities ( both for solar and mixtures ) are generated for the same pattern of abundances adopted in the stellar models , both at low and high temperatures .
the importance of this issue is illustrated by fig .
[ compare_opacity_hr ] , which shows the effect in the hertzsprung - russell ( hr ) diagram of adopting different ( inconsistent ) opacities and/or mixtures and/or enhancement factors as input to a 1  @xmath14 main sequence model with @xmath32 $ ] .
some of the illustrated cases are often found in literature .
the tracks are evolved up to the solar age of 4.6 gyr .
the meaning of the labels is the following : the first label refers to high temperature opacity ( @xmath33 ) and can vary from opal - s , i.e.  solar - scaled mixture , to opal-@xmath4 , i.e.  enhanced according to salaris & weiss ( 1998 ; w95 ) or rogers ( priv .
communication ; r95 ) .
this latter case is a special @xmath4-enhanced mixture characterized , if compared with salaris & weiss ( 1998 ) , by more c , less o , more ne , more na ( which is not an alpha - element ) , less al , more s , and more fe .
the second label refers to the molecular opacity , with the same notation as for solar and @xmath4-enhanced mixtures .
af stands for alexander & ferguson ( 1994 ) .
the prescription enclosed in a box represents the one adopted by us .
another point to emphasize is the difference between the present opacities and those used by girardi et al .
two are the main sources of differences .
first the adoption of the itoh et al .
( 1983 ) conductive opacities instead of the older ones by hubbard & lampe ( 1969 ) .
interestingly , catelan et al .  ( 1996 ) show that itoh et al .
( 1983 ) conductive opacities should not be used in rgb models , because they are valid only for a liquid phase , and not for the physical conditions existing in the interiors of rgb stars .
however , we verified that the effect of using the itoh et al .
( 1983 ) opacities instead of the hubbard & lampe ( 1969 ) is to increase the core mass at the helium flash by 0.006  ( see also castellani et al .  2000 and cassisi et al .
1998 ) , which is a negligible effect .
the second difference lies in the numerical technique used to interpolate within the grids of the opacity tables . in this paper
we use the two - dimensional bi - rational cubic damped - spline algorithm ( see & weiss 1998 ; and spth 1973 ) , whereas girardi et al .
( 2000 ) adopted the smooth bi - parabolic spline interpolation ( ssi ) already introduced in our code by bressan et al .
we verified that the differences amount to about a few percent .    in general ,
evolutionary models of low mass differ very little by changing the opacity interpolation scheme
. however , massive stars ( say @xmath34 ) are much more sensitive to the interpolation algorithm .
this is shown in fig .
[ compare_tracks_6msun ] , where we plot two 4  and 6  models with @xmath35 $ ] calculated both with the ssi ( continuous lines ) and the damped - spline ( dashed lines ) scheme .
figure  [ compare_opacity_inner ] compares the effects of the different interpolations on the opacity at the physical conditions holding in the inner regions of a 6  star when the central hydrogen mass fraction is @xmath36 , i.e.  at the stage marked with a dot in fig .
[ compare_tracks_6msun ] . at this stage
the temperature at the schwarzschild border of the convective core is @xmath37 and the density is @xmath38 ( corresponding to the arrow in the fig .
[ compare_opacity_inner ] ) . for the temperatures indicated along the curves ,
these represent the opacity as a function of the quantity @xmath39 ( see rogers & iglesias 1992 for more details ) .
the continuous curve is based on the ssi interpolation , whereas the dashed curve refers to the method adopted in the present study .
the small difference in opacity at the core border is typical of the values found throughout the evolutionary sequences , and is responsible for the different hr patterns .
we notice that both interpolations are acceptable , because the percentage error associated to both of them is roughly 3 % , which is lower than the opacity uncertainty in the original opal tables ( the latter can be as high as 10 % at individual @xmath40 points ) .
this shows that even the fine opal grid for some applications is not sufficiently dense and more grid points might be needed .
this is also illustrated by the solar model by schlattl & weiss ( 1998 ) , where the inverse situation , i.e.  spline producing a nicer opacity run , is found .
the reason for the different luminosities in the high - mass models , is that in these stars the radiative gradient across the border of the convective core is shallow .
therefore , even slightly larger opacities at the schwarzschild border may cause an increase of the convective core and lead to higher luminosities .
this trend already begins during the core h - burning phase and remains during the subsequent evolution .
the effects are detectable starting from stars with mass larger than @xmath41   ( see fig .
[ compare_tracks_6msun ] ) and becomes dramatic over 60  .
it is worth recalling that the maximum mass of our stellar models is 20  .
energy losses by neutrinos are from haft et al .
compared with the previous ones of munakata et al .
( 1985 ) used by girardi et al .
( 2000 ) , neutrino cooling during the rgb is more efficient .
this causes an increase of 0.005  in the core mass at helium ignition .
we have checked this on a 1  model with @xmath42 .
differences of the same order have been found by other authors ( e.g.  haft et al .
1994 ; cassisi et al .  1998 ) .
the remaining physical input is the same as in girardi et al .
( 2000 ) , to whom we refer for further details . in summary , we follow the evolution of h , @xmath43he , @xmath44he , @xmath45c , @xmath46c , @xmath47n , @xmath48n , @xmath49o , @xmath50o , @xmath51o , @xmath52ne , @xmath53ne , mg , according to the nuclear reaction rates of caughlan & fowler ( 1988 ) and landr et al .
for mg we refer to the total content of @xmath54 mg , @xmath55 mg and @xmath56 mg .
the equation of state ( eos ) is that of a fully - ionized gas for temperatures higher that @xmath57  k. for lower temperatures we adopt the mihalas et al .
( 1990 and references therein ) eos .
calibration of the solar model fixes the value of the mixing length parameter @xmath4 of the bhm - vitense ( 1958 ) theory at 1.68 . overshooting both from the convective core and the convective envelope
are taken into account according to the ballistic algorithm by bressan et al .
table [ overshooting ] summarizes the values of the overshooting parameters for the core , @xmath58 , and envelope , @xmath59 we have adopted as a function of the initial mass .
as far as the mass - loss is concerned , this is taken into account and its effect on the stellar structure is calculated for masses higher than 6  .
we used the same prescription adopted by fagotto et al .
( 1994 ) , where more details can be found .
c|c|c @xmath60 & @xmath58 & @xmath59 + 0.5 - 1.0 & 0 & 0.25 + 1.1 - 1.4 & @xmath61 & 0.25 + 1.5 - 2.0 & 0.5 & 0.25 + 2.0 - 2.5 & 0.5 & @xmath62 + 2.5 - 20 & 0.5 & 0.7 +      for the sake of illustration , in figs .
[ hrtracksz008 ] , [ hrtracksz008_low ] and [ hrtracksz008_low_hb ] we present the complete set of evolutionary tracks for @xmath0 $ ] and @xmath4-enhanced mixtures .
[ hrtracksz008 ] shows the models of intermediate - mass stars from the zero age main sequence ( zams ) up to the tp - agb phase and those of massive stars from the zams up to carbon ignition .
[ hrtracksz008_low ] presents the low - mass tracks from the zams up to the rgb - tip .
finally , fig .
[ hrtracksz008_low_hb ] shows the corresponding he - burning phase from the zero age horizontal branch ( zahb ) up to the tp - agb phase .
some relevant physical quantities are given in table  [ timez008 ] , which lists , as a function of the initial mass , the lifetimes in years of the central h - burning and he - burning phases , @xmath63 and @xmath64 respectively , the core mass at the helium ignition , @xmath65 , and the core mass at the first thermal pulse on the agb phase , @xmath66 .
for intermediate - mass and massive stars , @xmath65 refers to the core mass at the he - ignition . for the same chemical composition table  [ dredgeup ]
reports the changes in the surface chemical abundances of those elements that suffered from nuclear processing via the pp - chain and cno - cycle in the deep interiors , after the completion of the first and the second dredge - up by envelope mixing .
the surface abundances of heavier elements are not altered by nuclear processing .
the last two columns give the ratios of @xmath47n and @xmath49o with respect to their initial values .
these ratios are indicative of the efficiency of the dredge - up episodes .
table  [ critical_masses ] displays the transition masses @xmath67 , @xmath68 and @xmath69 for models calculated with solar and initial composition . in brief
, @xmath67 is the maximum mass below which stars become fully convective during the central h - burning phase , @xmath68 is the maximum mass for the core he - flash to occur , and finally @xmath69 is the maximum mass for the central c - burning to occur in highly electron degenerate gas ( deflagration or detonation followed by disruption of the star ) .
c|cccc|cccc & & + @xmath60 & @xmath70yr & @xmath71yr & @xmath65 & @xmath72 & @xmath70yr & @xmath71yr & @xmath65 & @xmath72 + 0.6 & 5.45  @xmath73 & 1.05  @xmath74 & 0.4814 & 0.5108 & 4.48  @xmath73 & 1.08  @xmath75 & 0.4445 & 0.5398 + 0.7 & 3.14  @xmath73 & 1.04  @xmath74 & 0.4806 & 0.5159 & 2.45  @xmath73 & 1.05  @xmath75 & 0.4434 & 0.5614 + 0.8 & 1.87  @xmath73 & 1.03  @xmath74 & 0.4793 & 0.5192 & 1.42  @xmath73 & 1.02  @xmath75 & 0.4434 & 0.5756 + 0.9 & 1.18  @xmath73 & 1.02  @xmath74 & 0.4781 & 0.5258 & 8.28  @xmath76 & 1.01  @xmath75 & 0.4430 & 0.5926 + 1.0 & 7.64  @xmath77 & 1.02  @xmath74 & 0.4782 & 0.5272 & 5.11  @xmath76 & 1.01  @xmath75 & 0.4426 & 0.5769 + 1.1 & 5.26  @xmath77 & 1.01  @xmath74 & 0.4776 & 0.5317 & 3.64  @xmath76 & 9.92  @xmath57 & 0.4424 & 0.6049 + 1.2 & 3.90  @xmath77 & 1.00  @xmath74 & 0.4778 & 0.5312 & 2.89  @xmath76 & 9.87  @xmath57 & 0.4418 & 0.6079 + 1.3 & 3.19  @xmath77 & 9.83  @xmath78 & 0.4778 & 0.5291 & 2.39  @xmath76 & 9.71  @xmath57 &
0.4421 & 0.6049 + 1.4 & 2.68  @xmath77 & 9.79  @xmath78 & 0.4763 & 0.5314 & 2.01  @xmath76 & 9.89  @xmath57 & 0.4375 & 0.6064 + 1.5 & 2.28  @xmath77 & 1.07  @xmath74 & 0.4642 & 0.5273 & 1.71  @xmath76 & 1.11  @xmath75 & 0.4169 & 0.6022 + 1.6 & 1.87  @xmath77 & 1.14  @xmath74 & 0.4539 & 0.5245 & 1.40  @xmath76 & 1.23  @xmath75 & 0.3999 & 0.5947 + 1.7 & 1.57  @xmath77 & 1.26  @xmath74 & 0.4391 & 0.5194 & 1.16  @xmath76 & 1.36  @xmath75 & 0.3845 & 0.5875 + 1.8 & 1.33  @xmath77 & 1.44  @xmath74 & 0.4192 & 0.5065 & 9.76  @xmath75 & 1.46  @xmath75 & 0.3680 & 0.5913 + 1.9 & 1.15  @xmath77 & 3.10  @xmath74 & 0.3337 & 0.4893 & 8.28  @xmath75 & 1.70  @xmath75 & 0.3406 & 0.5849 + 2.0 & 9.90  @xmath74 & 2.71  @xmath74 & 0.3451 & 0.4909 & 7.11  @xmath75 & 2.15  @xmath75 & 0.3281 & 0.5730 + 2.2 & 7.65  @xmath74 & 2.10  @xmath74 & 0.3757 & 0.5104 & 5.33  @xmath75 & 1.71  @xmath75 & 0.3598 & 0.6114 + 2.5 & 5.44  @xmath74 & 1.44  @xmath74 & 0.4246 & 0.5452 & 3.66  @xmath75 & 1.21  @xmath75 & 0.4104 & 0.6117 + 3.0 & 3.39  @xmath74 & 7.36  @xmath78 & 0.5367 & 0.6483 & 2.16  @xmath75 & 6.60  @xmath57 & 0.5178 & 0.6844 + 3.5 & 2.31  @xmath74 & 3.98  @xmath78 & 0.6718 & 0.7806 & 1.40  @xmath75 & 3.77  @xmath57 & 0.6408 & 0.6844 + 4.0 & 1.67  @xmath74 & 2.48  @xmath78 & 0.8213 & 0.9375 & 9.71  @xmath57 & 2.34  @xmath57 & 0.7713 & 0.6844 + 4.5 & 1.27  @xmath74 & 1.67  @xmath78 & 0.9744 & 1.0910 & 7.12  @xmath57 & 1.56  @xmath57 & 0.9092 & 0.6844 + 5.0 & 1.00  @xmath74 & 1.21  @xmath78 & 1.1333 & 1.2497 & 5.45  @xmath57 & 1.10  @xmath57 & 1.0593 & 1.2062 + 6.0 & 6.72  @xmath78 & 6.73  @xmath79 & 1.4456 & - & 3.53  @xmath57 & 6.06  @xmath80 & 1.3904 & 1.5452 + 7.0 & 4.87  @xmath78 & 4.55  @xmath79 & 1.7673 & - & 2.51  @xmath57 & 3.89  @xmath80 & 1.7609 & - + 8.0 & 3.76  @xmath78 & 3.30  @xmath79 & 2.1265 & - & 1.91  @xmath57 & 2.74  @xmath80 & 2.1599 & - + 9.0 & 3.03  @xmath78 & 2.52  @xmath79 & 2.4985 & - & 1.53  @xmath57 & 2.06  @xmath80 & 2.5988 & - + 10.0 & 2.51  @xmath78 & 2.02  @xmath79 & 2.9138 & - & 1.28  @xmath57 & 1.62  @xmath80 & 3.0935 & - + 12.0 & 1.86  @xmath78 & 1.43  @xmath79 & 3.7804 & - & 9.59  @xmath80 & 1.14  @xmath80 & 3.9767 & - + 15.0 & 1.34  @xmath78 & 1.03  @xmath79 & 4.6706 & - & 7.10  @xmath80 & 8.26  @xmath81 & 5.3587 & - + 20.0 & 9.39  @xmath79 & 7.12  @xmath82 & 6.7187 & - & 5.17  @xmath80 & 6.01  @xmath81 & 7.7582 & - +    @xmath68 is here defined as the initial mass for which the core mass at he ignition has its minimum value ( see table  [ timez008 ] and fig .
[ mflash ] ) .
[ mflash ] shows that similar values for @xmath68 are found both in the and in the solar - scaled case .
compared with girardi et al .
( 2000 ) , we find higher @xmath69 values .
this is once more caused by the different interpolation schemes of the opacity tables .
the one currently in use yields smaller convective cores for the same initial mass .
it follows that with the present overshooting prescription the maximum mass for a star to develop an electron - degenerate core after the core he - burning phase is slightly higher than previously estimated .
ccccccccccccc @xmath83 & h & @xmath43he & @xmath44he & @xmath45c & @xmath46c & @xmath47n & @xmath48n & @xmath49o & @xmath50o & @xmath51o & @xmath84 & @xmath85 +   + all & 0.742 & 3.80@xmath86 & 0.250 & 1.37@xmath87 & 1.65@xmath86 & 4.24@xmath88 & 1.67@xmath89 & 3.85@xmath87 & 1.56@xmath89 & 8.68@xmath89 & 1.00 & 1.00 +   + 0.60 & 0.730 & 3.95@xmath87 & 0.258 & 1.37@xmath87 & 1.72@xmath86 & 4.24@xmath88 & 1.66@xmath89 & 3.85@xmath87 & 1.56@xmath89 & 8.68@xmath89 & 1.00 & 1.00 + 0.70 & 0.727 & 2.73@xmath87 & 0.263 & 1.36@xmath87 & 2.50@xmath86 & 4.26@xmath88 & 1.55@xmath89 & 3.85@xmath87 & 1.56@xmath89 & 8.68@xmath89 & 0.993 & 1.01 + 0.80 & 0.725 & 1.97@xmath87 & 0.265 & 1.32@xmath87 & 3.71@xmath86 & 4.59@xmath88 & 1.44@xmath89 & 3.85@xmath87 & 1.56@xmath89 & 8.64@xmath89 & 0.964 & 1.08 + 0.90 & 0.724 & 1.55@xmath87 & 0.267 & 1.27@xmath87 & 4.02@xmath86 & 5.18@xmath88 & 1.36@xmath89 & 3.85@xmath87 & 1.57@xmath89 & 8.52@xmath89 & 0.925 & 1.22 + 1.00 & 0.723 & 1.20@xmath87 & 0.267 & 1.22@xmath87 & 4.13@xmath86 & 5.72@xmath88 & 1.30@xmath89 & 3.85@xmath87 & 1.60@xmath89 & 8.30@xmath89 & 0.891 & 1.35 + 1.10 & 0.723 & 1.01@xmath87 & 0.268 & 1.17@xmath87 & 4.24@xmath86 & 6.35@xmath88 & 1.22@xmath89 & 3.85@xmath87 & 1.70@xmath89 & 8.04@xmath89 & 0.851 & 1.50 + 1.20 & 0.725 & 8.67@xmath88 & 0.267 & 1.14@xmath87 & 4.38@xmath86 & 6.70@xmath88 & 1.18@xmath89 & 3.85@xmath87 & 1.83@xmath89 & 7.89@xmath89 & 0.829 & 1.58 + 1.30 & 0.726 & 8.00@xmath88 & 0.265 & 1.11@xmath87 & 4.18@xmath86 & 7.00@xmath88 & 1.15@xmath89 & 3.85@xmath87 & 2.01@xmath89 & 7.72@xmath89 & 0.811 & 1.65 + 1.40 & 0.724 & 6.89@xmath88 & 0.268 & 1.02@xmath87 & 4.59@xmath86 & 8.06@xmath88 & 1.02@xmath89 & 3.85@xmath87 & 2.56@xmath89 & 7.28@xmath89 & 0.742 & 1.90 + 1.50 & 0.723 & 6.11@xmath88 & 0.268 & 9.81@xmath88 & 4.49@xmath86 & 8.53@xmath88 & 9.72@xmath90 & 3.84@xmath87 & 1.02@xmath86 & 7.02@xmath89 & 0.716 & 2.01 + 1.60 & 0.724 & 5.35@xmath88 & 0.268 & 9.56@xmath88 & 4.47@xmath86 & 9.02@xmath88 & 9.44@xmath90 & 3.82@xmath87 & 1.02@xmath86 & 6.86@xmath89 & 0.697 & 2.13 + 1.70 & 0.722 & 4.76@xmath88 & 0.269 & 9.34@xmath88 & 4.50@xmath86 & 9.66@xmath88 & 9.15@xmath90 & 3.77@xmath87 & 1.70@xmath86 & 6.76@xmath89 & 0.681 & 2.28 + 1.80 & 0.720 & 4.21@xmath88 & 0.272 & 9.04@xmath88 & 4.39@xmath86 & 1.04@xmath87 & 8.88@xmath90 & 3.72@xmath87 & 1.93@xmath86 & 6.60@xmath89 & 0.660 & 2.46 + 1.90 & 0.718 & 3.78@xmath88 & 0.274 & 9.01@xmath88 & 4.45@xmath86 & 1.09@xmath87 & 8.76@xmath90 & 3.67@xmath87 & 1.63@xmath86 & 6.55@xmath89 & 0.658 & 2.57 + 2.00 & 0.718 & 3.38@xmath88 & 0.274 & 8.95@xmath88 & 4.41@xmath86 & 1.12@xmath87 & 8.72@xmath90 & 3.64@xmath87 & 1.63@xmath86 & 6.52@xmath89 & 0.653 & 2.65 + 2.20 & 0.714 & 2.76@xmath88 & 0.278 & 8.88@xmath88 & 4.51@xmath86 & 1.20@xmath87 & 8.57@xmath90 & 3.57@xmath87 & 1.53@xmath86 & 6.45@xmath89 & 0.648 & 2.82 + 2.50 & 0.711 & 2.13@xmath88 & 0.280 & 8.72@xmath88 & 4.53@xmath86 & 1.27@xmath87 & 8.37@xmath90 & 3.51@xmath87 & 1.16@xmath86 & 6.36@xmath89 & 0.636 & 3.00 + 3.00 & 0.713 & 1.50@xmath88 & 0.279 & 8.65@xmath88 & 4.47@xmath86 & 1.30@xmath87 & 8.27@xmath90 & 3.48@xmath87 & 9.83@xmath89 & 6.30@xmath89 & 0.631 & 3.08 + 3.50 & 0.717 & 1.13@xmath88 & 0.275 & 8.78@xmath88 & 4.66@xmath86 & 1.27@xmath87 & 8.23@xmath90 & 3.50@xmath87 & 6.74@xmath89 & 6.42@xmath89 & 0.641 & 2.99 + 4.00 & 0.718 & 9.02@xmath86 & 0.274 & 8.80@xmath88 & 4.83@xmath86 & 1.27@xmath87 & 8.10@xmath90 & 3.50@xmath87 & 5.56@xmath89 & 6.43@xmath89 & 0.642 & 2.99 + 4.50 & 0.720 & 7.69@xmath86 & 0.272 & 9.13@xmath88 & 4.77@xmath86 & 1.22@xmath87 & 8.49@xmath90 & 3.51@xmath87 & 3.89@xmath89 & 6.64@xmath89 & 0.666 & 2.89 + 5.00 & 0.721 & 6.49@xmath86 & 0.270 & 9.49@xmath88 & 4.64@xmath86 & 1.18@xmath87 & 8.90@xmath90 & 3.51@xmath87 & 3.70@xmath89 & 6.73@xmath89 & 0.693 & 2.79 + 6.00 & 0.722 & 4.75@xmath86 & 0.270 & 8.75@xmath88 & 4.96@xmath86 & 1.28@xmath87 & 7.79@xmath90 & 3.50@xmath87 & 4.04@xmath89 & 6.40@xmath89 & 0.638 & 3.01 + 7.00 & 0.723 & 3.85@xmath86 & 0.269 & 8.55@xmath88 & 5.00@xmath86 & 1.32@xmath87 & 7.51@xmath90 & 3.47@xmath87 & 3.84@xmath89 & 6.28@xmath89 & 0.624 & 3.12 + 8.00 & 0.721 & 3.26@xmath86 & 0.271 & 8.41@xmath88 & 4.98@xmath86 & 1.37@xmath87 & 7.31@xmath90 & 3.43@xmath87 & 3.61@xmath89 & 6.15@xmath89 & 0.614 & 3.24 + 9.00 & 0.722 & 3.06@xmath86 & 0.270 & 8.83@xmath88 & 4.99@xmath86 & 1.31@xmath87 & 7.83@xmath90 & 3.45@xmath87 & 2.77@xmath89 & 6.42@xmath89 & 0.644 & 3.08 + 10.00 & 0.721 & 2.78@xmath86 & 0.271 & 8.79@xmath88 & 4.92@xmath86 & 1.33@xmath87 & 7.79@xmath90 & 3.43@xmath87 & 2.83@xmath89 & 6.37@xmath89 & 0.641 & 3.13 + 12.00 & 0.713 & 2.36@xmath86 & 0.279 & 8.62@xmath88 & 5.03@xmath86 & 1.44@xmath87 & 7.52@xmath90 & 3.33@xmath87 & 2.49@xmath89 & 6.21@xmath89 & 0.629 & 3.40 + 15.00 & 0.687 & 1.90@xmath86 & 0.305 & 8.23@xmath88 & 4.94@xmath86 & 1.67@xmath87 & 7.15@xmath90 & 3.12@xmath87 & 2.24@xmath89 & 5.81@xmath89 & 0.601 & 3.94 +   + 4.00 & 0.694 & 8.64@xmath86 & 0.298 & 8.46@xmath88 & 4.71@xmath86 & 1.42@xmath87 & 7.84@xmath90 & 3.37@xmath87 & 5.91@xmath89 & 6.17@xmath89 & 0.618 & 3.35 + 4.50 & 0.676 & 6.95@xmath86 & 0.316 & 8.35@xmath88 & 4.65@xmath86 & 1.51@xmath87 & 7.72@xmath90 & 3.29@xmath87 & 4.59@xmath89 & 6.07@xmath89 & 0.609 & 3.56 + 5.00 & 0.662 & 5.71@xmath86 & 0.330 & 8.48@xmath88 & 4.51@xmath86 & 1.56@xmath87 & 7.90@xmath90 & 3.22@xmath87 & 4.24@xmath89 & 6.01@xmath89 & 0.619 & 3.68 + 6.00 & 0.667 & 4.35@xmath86 & 0.325 & 8.08@xmath88 & 4.76@xmath86 & 1.60@xmath87 & 7.29@xmath90 & 3.22@xmath87 & 3.82@xmath89 & 5.86@xmath89 & 0.590 & 3.77 + 7.00 & 0.714 & 3.75@xmath86 & 0.278 & 8.35@xmath88 & 4.89@xmath86 & 1.41@xmath87 & 7.38@xmath90 & 3.40@xmath87 & 3.89@xmath89 & 6.11@xmath89 & 0.609 & 3.32 +    c|c|c|c|c|c @xmath18 & @xmath91 & mixture & @xmath67 & @xmath68 & @xmath69 + 0.008 & 0.250 & solar - scaled & 0.36 & 1.90 & @xmath92 + 0.019 & 0.273 & solar - scaled & 0.37 & 2.10 & @xmath93 + 0.040 & 0.320 & solar - scaled & 0.36 & 2.20 & @xmath94 + 0.070 & 0.390 & solar - scaled & 0.33 & 1.90 & @xmath94 + 0.008 & 0.250 & & 0.37 & 1.90 & @xmath92 + 0.019 & 0.273 & & 0.37 & 2.10 & @xmath92 + 0.040 & 0.320 & & 0.36 & 2.10 & @xmath94 + 0.070 & 0.390 & & 0.33 & 2.00 & @xmath93 +
the four sets of stellar models are used to calculate isochrones , integrated magnitudes and colours of single stellar populations with ages from @xmath78 to @xmath95  yr both in the johnson - cousins @xmath6 and hst - wfpc2 photometric systems .
the age range is large enough to describe young clusters and associations as well as old globular clusters .
the dense grid of stellar tracks allows us to construct detailed isochrones at small age steps ( @xmath96 ) .
this also makes the data - base suitable for the simulation of synthetic colour - magnitude diagrams ( e.g.  girardi 1999 ) .
before theoretical isochrones are constructed , the tp - agb phase is included in all tracks of @xmath97 .
therefore , this phase is present in the isochrones older than about @xmath75  yr . the reader is referred to girardi & bertelli ( 1998 ) and girardi et al .
( 2000 ) for all details about the synthetic algorithm used to follow the tp - agb evolution .      the hr diagram of fig .
[ compare_isoc ] shows a comparison between isochrones with chemical composition @xmath98 $ ] both for solar - scaled ( continuous lines ) and ( dotted lines ) stellar models .
for the sake of clarity we plot only the isochrones with ages between @xmath99 and @xmath100 , in steps of @xmath101 . as already noted by salaris & weiss ( 1998 ) , for solar and super - solar total metallicity
, the element abundance ratios among metals affect the evolution . from the hr diagram of fig .
[ compare_isoc ] , two main features can be singled out .
\(i ) @xmath4-enhanced isochrones have fainter and hotter turn - offs ( to ) .
this is shown by the entries of table  [ turnoff ] , which lists the turn - off luminosity and effective temperature of isochrones with @xmath102 but different mixtures of @xmath4-elements .
the reason is that at given central hydrogen content the stellar model with solar - scaled composition has a mean opacity higher than the one .
the point is illustrated in fig .
[ compare_inner ] where we compare the opacity profile across the 1  models with @xmath103 $ ] and both   ( dashed line ) and solar - scaled ( solid line ) mixtures .
higher opacities induce steeper radiative temperature gradients @xmath104 and lower surface temperatures @xmath105 for the same central conditions .
furthermore , steeper @xmath104s imply smaller burning regions and lower luminosities in turn .
this effect overwhelms the one induced by the slightly larger convective core due to the higher @xmath104 ( when a convective core is present ) . combining all those effects together , at the same evolutionary stage
, the solar - scaled model is fainter , cooler and older than the one .
the same trend is recovered in the isochrones as well .
the implications of the different turn - off temperatures and luminosities on the ages of globular clusters have already been extensively investigated ( see e.g.  salaris & weiss 1998 ) . for isochrones younger than @xmath76  yr , the differences in the turn - off properties are much less marked .
this is due both to ( 1 ) the weak metallicity dependence of the opacities in massive stars ( where electron - scattering becomes the main source of opacity ) and ( 2 ) the larger extension of central convection caused by the flatter profile of @xmath104 in stellar interiors .
\(ii ) @xmath4-enhanced isochrones are hotter than the solar - scaled ones throughout all evolutionary phases .
the effect increases with the metallicity .
it is mainly due to the higher opacities of the solar - scaled mixtures .
looking at the 10  gyr isochrone with chemical composition @xmath103 $ ] as an example , the temperature difference @xmath106 is about @xmath270.015 at the turn - off and increases to @xmath270.020 at the bottom of the rgb .
these differences in effective temperatures result into colour differences @xmath107 of 0.05  mag and 0.08  mag , at the turn - off and bottom of the rgb , respectively .
c|c|c|c|c & & + age & @xmath108 & @xmath109 & @xmath108 & @xmath109 + 7.00 & 4.491 & 4.424 & 4.473 & 4.428 + 7.20 & 4.109 & 4.367 & 4.084 & 4.370 + 7.40 & 3.749 & 4.309 & 3.722 & 4.312 + 7.60 & 3.418 & 4.253 & 3.380 & 4.255 + 7.80 & 3.095 & 4.197 & 3.046 & 4.198 + 8.00 & 2.762 & 4.143 & 2.728 & 4.143 + 8.20 & 2.390 & 4.087 & 2.367 & 4.088 + 8.40 & 2.090 & 4.030 & 2.028 & 4.031 + 8.60 & 1.791 & 3.976 & 1.733 & 3.978 + 8.80 & 1.492 & 3.924 & 1.436 & 3.926 + 9.00 & 1.207 & 3.874 & 1.154 & 3.877 + 9.20 & 0.853 & 3.834 & 0.839 & 3.837 + 9.40 & 0.587 & 3.808 & 0.528 & 3.812 + 9.60 & 0.383 & 3.787 & 0.377 & 3.795 + 9.80 & 0.379 & 3.776 & 0.406 & 3.785 + 10.00 & 0.240 & 3.759 & 0.232 & 3.770 + 10.20 & 0.083 & 3.740 & 0.082 & 3.753 +    it is worth remarking that in weiss et al .
( 1995 ) , the isochrones are found to have hotter and brighter turn - offs , and cooler rgbs ( at least in large portions of it ) than in the solar - scaled case ( see their fig .  4 ; notice that the line - types there are wrong : solid lines refer to the case ) .
the cooler rgbs are no longer found here or in vandenberg et al .
( 2000 ) and have been a consequence of the lack of low - temperature opacity tables for compositions in opacities in weiss et al .
this confirms how important is to use self - consistent opacities .
finally , one can ask whether an isochrone may be approximated by a solar - scaled one with a slightly lower metallicity . in order to answer this question
let us examine a set of 10 gyr old isochrones . in fig .
[ fig_isoc_10gyr ] we plot all the isochrones with this age , and with metallicity going from @xmath12 to @xmath110 for both the ( continuous lines ) and solar - scaled ( short dashed lines ) mixtures .
one can see that both sets of isochrones have different shapes , the differences becoming larger at higher metallicities .
more specifically , one can notice that the difference in effective temperatures at the turn - off point between an isochrone and a solar - scaled one , @xmath111 , increases with the metallicity , passing from @xmath112 to @xmath113 as @xmath18 increases from 0.008 to 0.07 .
on the contrary , the temperature difference at the base of the rgb , @xmath114 , is almost constant and of about @xmath115 .
for the same age we calculate , by interpolating among the grids of different metallicity , several solar - scaled isochrones with @xmath18 spanning from @xmath116 to @xmath117 , and selected among them the one best fitting the main sequence of the @xmath102 isochrone .
this turns out to be the interpolated isochrone with @xmath118 , which is shown by the long - dashed line in fig .
[ fig_isoc_10gyr ] .
clearly , this solar - scaled @xmath118 isochrone does not reproduce well the rgb sequence of the @xmath102 one . more specifically , although @xmath111 is very small for this isochrone pair , @xmath114 is not neglegible ( @xmath119 ) .
from this simple test we then conclude that for a fixed age and relatively high metallicities , an isochrone can not be reproduced by tuning the metallicity of a solar - scaled one .
the above temperature differences for rgb models ( isochrones ) bear very much on integrated spectral line indices of old metal - rich populations , which are commonly used to infer ages and metallicities of elliptical galaxies .
this will be the subject of a forthcoming paper .
theoretical luminosities and effective temperatures are translated to magnitudes and colours by means of extensive tabulations of bolometric corrections ( bc ) and colours obtained from properly convolving spectral energy distributions ( seds ) as a function of @xmath105 , @xmath120 and [ fe / h ] .
the procedure is amply described in bertelli et al .
( 1994 ) , bressan et al .
( 1994 ) , and tantalo et al .
( 1996 ) to whom the reader should refer for more details .
suffice it to recall here the following basic steps :    \(i ) the main body of the spectral library is from kurucz ( 1992 ) , however extended to the high and low temperature ranges . for stars with @xmath121
k pure black - body spectra are assigned , whereas for stars with @xmath122 k the catalogue of stellar fluxes by fluks et al .
( 1994 ) is adopted .
this catalogue includes 97 observed spectra for all m - spectral subtypes in the wavelength range @xmath123 , and synthetic photospheric spectra in the range @xmath124 .    the scale of  in fluks et al .
( 1994 ) is similar to that of ridgway et al .
( 1980 ) for spectral types earlier than m4 but deviates from it for later spectral types . since ridgway s et al .
( 1980 ) scale does not go beyond the spectral type m6 , no comparison for more advanced spectral types is possible .
the problem is further complicated by possible effects of metallicity .
the ridgway scale of  is based on stars with solar metallicity ( @xmath125 ) and empirical calibrations of the -scale for @xmath126 are not available .    to cope with this difficulty
, we have introduced the metallicity-  relation of bessell et al.(1989 , 1991 ) using the ( @xmath127 ) colour as a temperature indicator .
an interpolation is made between the  of bessell et al .
( 1989 ) and the ( @xmath127 ) colours given by fluks et al .
( 1994 ) for the spectral types from m0 to m10 .
\(ii ) at a given metallicity , surface gravity @xmath128 , and , bcs are determined by convolving the corresponding stellar sed with the response functions of the various pass - bands .
\(iii ) we recall that kurucz does not provide atmospheric models for both solar - scaled and mixtures .
it is not therefore possible to estimate the impact ( if any ) of atmospheres on the isochrone colours and magnitudes .
[ [ johnson - cousins - photometry . ] ] johnson - cousins photometry .
+ + + + + + + + + + + + + + + + + + + + + + + + + + +    the response functions for the various pass - bands in which magnitudes and colours are generated are from the following sources : buser & kurucz ( 1978 ) for the @xmath129 pass - bands , bessell ( 1990 ) for the @xmath130 and @xmath131 cousins pass - bands , and finally bessell & brett ( 1988 ) for the @xmath132 pass - bands .
re - normalization of the colours obtained by convolving the seds with pass - band , has been made by convolving the sed of vega and imposing that the computed colours strictly match the observed ones ( kurucz 1992 ) .    finally , the zero point of the bcs is fixed by imposing that the bc for the kurucz model of the sun is @xmath133 .    [
[ hst - wfpc2-photometry . ] ] hst - wfpc2 photometry .
+ + + + + + + + + + + + + + + + + + + + +    the transformation to the wfpc2 photometric system ( filters f170w , f218w , f255w , f300w , f336w , f439w , f450w , f555w , f606w , f702w , f814w , and f850lp ) requires some additional explanations .
the main parameters defining this photometric system are given in table  [ hstfilter ] .
columns ( 1 ) to ( 3 ) report the filter name , the mean wavelength @xmath134 and the r.m.s .
band width ( bw ) , respectively .
c|c|c|c|c|c filter & @xmath134 (  ) & bw (  ) & @xmath135 & @xmath136 & note + ( 1 ) & ( 2 ) & ( 3 ) & ( 4 ) & ( 5 ) & ( 6 ) + f170w & 1747 & 290 & @xmath2720.6926 & @xmath270.4070 & @xmath27 + f218w & 2189 & 171 & @xmath2720.8737 & @xmath270.2263 & @xmath27 + f255w & 2587 & 170 & @xmath2721.0940 & @xmath270.0060 & @xmath27 + f300w & 2942 & 325 & @xmath2721.1336 & + 0.0336 & wide u + f336w & 3341 & 204 & @xmath2721.1876 & + 0.0876 & u - john .
+ f439w & 4300 & 202 & @xmath2720.4325 & @xmath270.6675 & b - john .
+ f450w & 4519 & 404 & @xmath2720.6339 & @xmath270.4661 & wide b + f555w & 5397 & 522 & @xmath2721.0798 & @xmath270.0202 & v - john .
+ f606w & 5934 & 637 & @xmath2721.4250 & + 0.3250 & wide v + f702w & 6862 & 587 & @xmath2721.8594 & + 0.7594 & wide r + f814w & 7924 & 647 & @xmath2722.3309 & + 1.2309 & i - john .
+ f850lp & 9070 & 434 & @xmath2722.7578 & + 1.6578 & @xmath27 +    with the wfpc2 detector at least three kinds of magnitudes are commonly in use : @xmath137 , @xmath138 and @xmath139 ( see the _ synphot user s guide _ by white et al .
1998 , distributed by the stsdas group ) . in the following we
limit ourselves to consider only the theoretical counterparts of @xmath137 and @xmath138 .    for radiation with
flat spectrum and specific flux @xmath140 ( in @xmath141 ) impinging on the telescope , the @xmath137 magnitude is defined by @xmath142 which means that @xmath143 for @xmath144 @xmath145 .
let us then define the pass - band @xmath146 as the product of the filter transmission by the response function of the telescope assembly and detector in use . in this case , @xmath140 is related to the _ counts _ rate @xmath147
i.e.  the number of photons per second registered by the detector  through @xmath148 where @xmath149 is the effective collecting surface of the telescope , and @xmath150 and @xmath151 are the planck constant and speed of light , respectively .
@xmath147 , the quantity actually measured at the telescope , can be easily converted to a specific flux by means of eq .
( [ countrate ] ) .
therefore , specific fluxes @xmath140 and @xmath152 magnitudes can be defined for incoming spectra of any shape . in the case of a stellar source of radius @xmath130 , located at a distance @xmath153 , and emitting a specific flux @xmath154 , eq .
( [ stmagdef ] ) can be replaced by @xmath155    \label{stmagdefstar}\ ] ] which defines the apparent @xmath152 magnitude of a star .    in our case
, we derive the absolute @xmath156 magnitudes by simply locating our synthetic stars at a distance of @xmath157  pc in eq .
( [ stmagdefstar ] ) .
the specific fluxes @xmath158 are taken from the library of synthetic spectra in use as a function of @xmath105 , surface gravity @xmath159 , and chemical composition .
similarly , we can define the @xmath160 magnitudes @xmath161    \label{vegamagdefstar}\ ] ] by imposing that the zero points @xmath135 are such that the synthetic magnitudes of vega match the ground - based @xmath162 johnson apparent magnitudes for this star , for the filters which are closest in wavelength to each other ( see column 6 in table  [ hstfilter ] ) .
we adopt the values 0.02 , 0.02 , 0.03 , 0.039 , 0.035  mag for the apparent magnitudes of vega in the @xmath162 johnson system ( see holtzman et al .
1995 ) . for the uv filters
, vega is assumed to have an apparent magnitude equal to zero .
we remind the reader that if we assume the apparent vega magnitude to be equal to zero in _ all _ the pass - bands , we would obtain @xmath163 and eq .  ( [ vegamagdefstar ] ) would become @xmath164 this latter is the definition of _ vegamag _ in the synphot package distributed with the stsdas software . in this paper
we do not make use of this calibration , but prefer to follow the slightly different one defined by holtzman et al .
( 1995 ) and described above .    to get the final calibration
, we insert in eq .
[ vegamagdefstar ] the specific flux of vega calculated by castelli & kurucz ( 1994 ) with atmosphere models , and converted it to an absolute flux on the earth assuming for vega the angular diameter of @xmath165  mas ( code et al .  1976 ) .
the absolute flux we have adopted differs ( a few percent ) from the one by hayes ( 1985 ) used by holtzman et al .  (
once more , the absolute magnitude @xmath166 is obtained locating the synthetic star at a 10  pc distance .
column ( 4 ) of table [ hstfilter ] lists the zero points we have obtained , whereas column ( 5 ) gives the difference @xmath167 between the zero point of the @xmath137 and @xmath138 magnitudes .
the conversion from one magnitude scale to another is thus possible .
finally , we calculate the bolometric corrections bc which allow one to convert theoretical bolometric magnitudes both to @xmath168 and @xmath169 : @xmath170 \nonumber \\    & + & 21.1   + 2.5 \log \frac{\int p(\lambda ) \,i_{\rm k}(\lambda)\ , \lambda   \,{\mbox{d}}\lambda}{\int p(\lambda)\ , \lambda \,{\mbox{d}}\lambda } \\ { \rm bc}(p)^{m_{\rm vega}}_{t_{\rm eff},g } & = & { \rm bc}(p)^{m_{\rm st}}_{t_{\rm eff},g } -21.1 - z_p \label{bc_stmag}\end{aligned}\ ] ] where @xmath171 is the stefan - boltzmann constant . for the solar luminosity and absolute bolometric magnitude
we adopt @xmath172 ( bahcall et al .
1995 ) and @xmath173 .
owing to the presence of contaminants inside the wfpc2 ( see baggett & gonzaga 1998 and holtzman et al .  1995 ) , @xmath146 changes slowly with time ; in particular the uv throughput degrades , changing the photometric performances . to cope with this drawback
, small corrections are usually added to the definition of the instrumental magnitudes in order to bring the magnitudes back to the optimal conditions ( see holtzman et al .
1995 ) . because of this , in the calculation of the bolometric corrections we do not use the present - day pass - bands @xmath146 provided by synphot but insert the pre - launch pass - bands , as in holtzman et al .
( 1995 ) .
in this section we present the integrated magnitudes , colours and mass - to - light ratios for single stellar populations ( ssp ) of the same age and chemical compositions of the isochrones above , both in the johnson - cousins and wfpc2 photometric systems .
the integrated magnitudes and colours are computed assuming an initial mass function @xmath174 and adopting the normalization @xmath175 where @xmath176 and @xmath177 are the lower and upper limits of zero age main sequence stars . by doing so
, the integrated magnitudes of these ideal ssps can be easily scaled to produce integrated magnitudes for stellar populations of any arbitrary total _ initial _
mass @xmath178 , by simply adding the factor @xmath179 .
the theoretical mass - to - light ratios @xmath180 are calculated following the same procedure as in alongi & chiosi ( 1990 ) and chandar et al .
( 1999 ) .
the fate of a single star depends on its initial mass ( binary stars are neglected here ) . in a brief and over - simplified picture of stellar evolution
the following mass limits and groupings can be identified : ( i ) stars more massive than @xmath181 have a lifetime shorter than the current estimate of the hubble age of the universe , say @xmath182 gyr .
stars lighter than this limit are not of interest here because once formed they live for ever .
( ii ) stars more massive than @xmath69 explode as supernovae leaving a neutron star remnant of @xmath183 .
the possibility that massive stars ( @xmath184 ) may end up as black holes is not considered here .
( iii ) stars less massive than @xmath185 and more massive than @xmath186 terminate their evolution as white dwarfs of suitable masses @xmath187 that depend on the initial mass and efficiency of mass loss during the rgb and agb phases . with the current estimates of the mass loss efficiency @xmath188 , i.e.   @xmath189 .
the possibility , however , exists that @xmath190 .
in such a case stars in the mass interval @xmath185 to @xmath69 reach the c - ignition stage and deflagrate as supernovae leaving no remnant .
for the purposes of the present study , this latter case is neglected and @xmath191 is always assumed .
finally , let @xmath192 be the initial mass of the most massive star still alive in a ssp of a certain age .
it is worth recalling that all the above mass limits depend on the initial chemical composition and that this dependence can be properly taken into account .
therefore , at any given age @xmath193 the total mass @xmath194 of a ssp is made of two contributions : @xmath195 where the first right - hand side term is the total mass of the stars still alive at the age @xmath193 ( i.e.  burning a nuclear fuel ) and the second one is the total mass in stellar remnants , i.e.   white dwarfs or neutron stars depending on the initial mass of the progenitor ( i.e.  @xmath196 for @xmath197 , and @xmath198 for @xmath199 ) .
we adopt an upper mass limit @xmath200 .
the choice for the lower mass limit , @xmath176 , is explained below .
it is worth recalling here that , due to eq .
[ eq_imf ] , @xmath201 is the current mass of an ideal ssp with initial mass equal to @xmath202 .
@xmath201 decreases as a function of time from its initial value , due to the mass lost by both stellar winds ( for all masses ) and supernova explosions ( from massive stars ) .
finally , the total luminosity @xmath203 of the ssp in any pass - band @xmath204 is obtained by integrating , along the isochrone , the luminosity of individual stars whose number is given by the imf @xmath205    to proceed further , one has to specify the initial mass function and @xmath176 . here
, we adopt the salpeter ( 1955 ) law with slope @xmath206 over the whole mass range @xmath176 to @xmath177 . concerning @xmath176 , the choice @xmath207
is made by imposing that the observational value @xmath208 of the lmc cluster ngc  1866 is matched ( see girardi & bica 1993 ) .
the age and metallicity of ngc  1866 are @xmath74 years and @xmath12 , respectively , whereas its @xmath209 is close to the mean value of 0.20 for lmc clusters of the same age ( see battinelli & capuzzo - dolcetta 1989 , and references therein ) .
we have computed extended sets of evolutionary tracks and isochrones at relatively high metallicities . a major improvement , with respect to previous calculations of tracks , is that we made use of self - consistent opacities , i.e.  computed with the same chemical compositions over the complete range of temperatures .
the main result is that in general all evolutionary phases ( isochrones ) have higher effective temperature with respect to solar - scaled models ( isochrones ) of same metal content @xmath18 .
the temperature shift is caused by the lower opacities of the mixtures .
moreover , we find that at relatively high metallicities and old ages , an isochrone can not be mimicked by simply using a solar - scaled isochrone of lower metallicity .
obvious limitation of the present results is that they refer to a particular choice for the chemical compositions derived from observations of a particular sample of metal - poor field stars .
future observational data may suggest a different partition of metals in stars , and hence different results for the stellar models and corresponding isochrones .    testing a large number of possible combinations of ratios , though feasible ,
would not be of practical use owing to the present uncertainty on the abundance ratios .
therefore , the results of this study ought to be taken as indicative of the _ potential _ effect of element ratios .
moreover , it is likely that variations of single @xmath4-elements are less important than the whole problem of including or not the enhancement of @xmath4-elements in the initial chemical composition . finally , it is still not clear whether differences in enhancements are significant compared to the errors in the determinations .
[ [ retrieval - of - the - data - sets . ] ] retrieval of the data sets .
+ + + + + + + + + + + + + + + + + + + + + + + + + + +    complete tabulations of the relevant data for all evolutionary sequences , isochrones , integrated magnitudes and colours and mass - to - light rations together with useful summary tables can be obtained either upon request to the authors or downloaded from the address ` http://pleiadi.pd.astro.it ` .
the layout of the tables of stellar models and isochrones is the same as in girardi et al .
( 2000 ) as far as the johnson - cousins photometric system is concerned . for the wfpc2 system
the layout is similar , with table headers allowing to easily identify the pass - band .
we like to thank a.  bressan for many suggestions and useful conversations , h.  schlattl for his help with the opacity interpolation code , h.  aussel and m.  zoccali for their useful clarifications about the hst / wfpc2 photometric system , j.  holtzman for providing us with the pre - launch wfpc2 pass - bands , and the anonymous referee for useful suggestions .
l.g .  acknowledges support from the alexander von humboldt - stiftung during his stay at mpa .
thanks mpa for the warm hospitality and support .
this study has been funded by the italian ministry of university , scientific research and technology ( murst ) under contract `` formation and evolution of galaxies '' n.   9802192401 .
alexander d.r . , ferguson j.w . , 1994 , apj 437 , 879                                                                                schlattl h. , weiss a. , 1998 , in `` proc .
neutrino astrophysics '' , ringberg castle , tegernsee , germany , 20 - 24 oct 1997 , ed .
m.  altmann , w.  hillebrandt , h .- t.janka and g.raffelt ( sfb astroteilchenphysik , technical university munich )