the recent measurements of the wilkinson microwave anisotropy probe ( wmap ) satellite @xcite on the cosmic microwave background radiation ( cmbr ) have shown that the matter abundance in the universe is @xmath2 , where @xmath3 with @xmath4 being the energy density of the @xmath5-th species and @xmath6 the critical energy density of the universe and @xmath7 is the present value of the hubble parameter in units of @xmath8 .
the baryon abundance is also found by these measurements to be @xmath9 .
consequently , the cold dark matter ( cdm ) abundance in the universe is @xmath10 . the @xmath11 confidence level ( c.l . )
range of this quantity is then , roughly , @xmath12 .
taking @xmath13 , which is its best - fit value from the hubble space telescope ( hst ) @xcite , and assuming that the total energy density of the universe is very close to its critical energy density ( i.e. @xmath14 ) , as implied by wmap , we conclude that about @xmath15 of the energy density of the present universe consists of cdm .
the question then is , what the nature , origin , and composition of this important component of our universe is .
particle physics provides us with a number of candidate particles out of which cdm can be made .
these particles appear naturally in various particle physics frameworks for reasons completely independent from cdm considerations and are , certainly , not invented for the sole purpose of explaining the presence of cdm in the universe .
the basic properties that such a candidate particle must satisfy are the following : ( @xmath5 ) it must be stable or very long - lived , which can be achieved by an appropriate symmetry , ( @xmath16 ) it should be electrically and color neutral , as implied by astrophysical constraints on exotic relics ( like anomalous nuclei ) , but can be interacting weakly , and ( @xmath17 ) it has to be non - relativistic , which is usually guaranteed by assuming that it is adequately massive , although even very light particles such as axions can be non - relativistic for different reasons .
so , what we need as constituent of cdm is a weakly interacting massive particle .
there are several possibilities , but we will concentrate here on the major particle physics candidates which are the axion , the lightest neutralino , the axino , and the gravitino ( for other candidates , see e.g. ref .
note that the last three particles exist only in supersymmetric ( susy ) theories .    in sec .
[ sec : axion ] , we examine the possibility that the axions are constituents of cdm .
[ sec : mssm ] is devoted to outlining the salient features of the minimal supersymmetric standard model ( mssm ) , which will be used as a basic frame for discussing susy cdm . in sec .
[ sec : relic ] , we summarize the calculation of the relic abundance of the lightest neutralino , which is normally the lightest supersymmetric particle ( lsp ) , and investigate the circumstances under which it can account for the cdm in the universe . in secs .
[ sec : axino ] and [ sec : gravitino ] , we discuss , respectively , axinos and gravitinos as constituents of cdm . in sec .
[ sec : quasi ] , we present a susy grand unified theory ( gut ) model which solves the bottom - quark mass problem by naturally and modestly violating the exact unification of the third generation yukawa couplings .
we study the parameter space of the model which is allowed by neutralino dark matter considerations as well as some other phenomenological constraints . finally , in sec .
[ sec : concl ] , we summarize our conclusions .
the most natural solution to the strong @xmath18 problem ( i.e. the apparent absence of @xmath18 violation in strong interactions implied by the experimental bound on the electric dipole moment of the neutron ) is the one provided by a peccei - quinn ( pq ) symmetry @xcite . this is a global @xmath19 symmetry , @xmath20 , which carries qcd anomalies and is spontaneously broken at a scale @xmath21 , the so - called axion decay constant or simply pq scale .
astrophysical @xcite and cosmological constraints imply that @xmath22 .
the upper bound originates @xcite from the requirement that the relic energy density of axions does not overclose the universe
. it should be noted , however , that this upper bound can be considerably relaxed if the axions are diluted @xcite by entropy generation after their production at the qcd phase transition ( for more recent applications of this dilution mechanism , see e.g. ref .
@xcite ) .
the axion is a pseudo nambu - goldstone boson corresponding to the phase of the complex pq field , which breaks @xmath20 by its vacuum expectation value ( vev ) . after the end of inflation @xcite
, this phase appears homogenized over the universe ( supposing that the pq field is non - zero ) with a value @xmath23 , which is known as the initial misalignment angle .
naturalness suggests that @xmath23 is of order unity .
this angle remains frozen until the qcd phase transition , where the qcd instantons come into play .
they break explicitly the pq symmetry to a discrete subgroup @xcite since this symmetry carries qcd anomalies .
so , a sinusoidal potential for the phase of the pq field is generated and this phase starts oscillating coherently about a minimum of the potential .
the resulting state resembles pressureless matter consisting of static axions with mass @xmath24 , where @xmath25 is the qcd scale . for @xmath26 , the mass of the axion @xmath27 .
note that axions , although very light , are good candidates for being constituents of the cdm in the universe since they are produced at rest . also ,
they are very weakly interacting since their interactions are suppressed by the axion decay constant @xmath21 .
the relic abundance of axions can be calculated by using the formulae of ref .
@xcite , where we take the qcd scale @xmath28 and ignore the uncertainties for simplicity .
we find @xmath29 ( note that a primordial magnetic helicity , may @xcite influence this abundance ) .
so , for natural values of @xmath30 and @xmath26 , axions can contribute significantly to cdm , which can even consist solely of axions .
the main disadvantage of axionic dark matter is that it leads to isocurvature perturbations if the pq field emerges with non - zero ( homogeneous ) value at the end of inflation .
indeed , during inflation , the angle @xmath23 acquires a superhorizon spectrum of perturbations as all the almost massless degrees of freedom . at the qcd phase transition
, these perturbations turn into isocurvature perturbations in the axion energy density , which means that the partial curvature perturbation in axions is different than the one in photons .
the recent results of wmap @xcite put stringent bounds @xcite on the possible isocurvature perturbation .
so , a large axion contribution to cdm is disfavored in models where the inflationary scale is superheavy ( i.e. of the order of the susy gut scale ) and the pq field is non - zero at the end of inflation .
we now wish to turn to the discussion of the main susy candidates for dark matter : the lightest neutralino @xmath31 , the axino @xmath32 and the gravitino @xmath33 .
we will consider them mainly within the simplest susy framework , which is the mssm .
it is , thus , important to first outline the basics of mssm .
we consider the mssm embedded in some general susy gut model .
we further assume that the gut gauge group breaking down to the standard model ( sm ) gauge group @xmath34 occurs in one step at a scale @xmath35 , where the gauge coupling constants of strong , weak , and electromagnetic interactions unify .
ignoring the yukawa couplings of the first and second generation , the effective superpotential below @xmath36 is @xmath37 where @xmath38 and @xmath39 are the quark and lepton @xmath40 doublet left handed superfields of the third generation and @xmath41 , @xmath42 , and @xmath43 the corresponding @xmath40 singlets .
also , @xmath44 , @xmath45 are the electroweak higgs superfields and @xmath46 the @xmath47 antisymmetric matrix with @xmath48 .
the gravity - mediated soft susy - breaking terms in the scalar potential are given by @xmath49 @xmath50 , \label{vsoft}\ ] ] where the sum is taken over all the complex scalar fields @xmath51 and tildes denote superpartners .
the soft gaugino mass terms in the lagrangian are @xmath52 where @xmath53 , @xmath54 and @xmath55 are the bino , winos and gluinos respectively .
the lagrangian of mssm is invariant under a discrete @xmath56 matter parity symmetry under which all `` matter '' ( i.e. quark and lepton ) superfields change sign . combining this symmetry with the @xmath56 fermion number symmetry under which all fermions change sign
, we obtain the discrete @xmath56 r - parity symmetry under which all sm particles are even , while all sparticles are odd . by virtue of r - parity conservation , the lsp is stable and , thus , can contribute to the cdm in the universe .
it is important to note that matter parity is vital for mssm to avoid baryon- and lepton - number - violating renormalizable couplings in the superpotential , which would lead to highly undesirable phenomena such as very fast proton decay .
so , the possibility of having the lsp as cdm candidate is not put in by hand , but arises naturally from the very structure of mssm .
the susy - breaking parameters @xmath57 , and @xmath58 ( @xmath59 ) are all of the order of the soft susy - breaking scale @xmath60 , but are otherwise unrelated in the general case . however ,
if we assume that soft susy breaking is mediated by minimal supergravity ( msugra ) , i.e. supergravity with minimal khler potential , we obtain soft terms which are universal `` asymptotically '' ( i.e. at @xmath36 ) .
in particular , we obtain a common scalar mass @xmath61 , a common trilinear scalar coupling @xmath62 , and a common gaugino mass @xmath63 . the mssm supplemented by universal boundary conditions
is known as constrained mssm ( cmssm ) @xcite .
it is true that msugra implies two more asymptotic relations : @xmath64 and @xmath65 , where @xmath66 and @xmath67 is the ( asymptotic ) gravitino mass .
these extra conditions are usually not included in the cmssm . imposing them , we get the so - called very cmssm @xcite , which is a very restrictive version of mssm and
will not be considered in these lectures .
the cmssm can be further restricted by imposing asymptotic yukawa unification ( yu ) @xcite , i.e. the equality of all three yukawa coupling constants of the third family at @xmath36 : @xmath68 exact yu , which makes the cmssm considerably more predictive , can be obtained in guts based on a gauge group such as @xmath69 or @xmath70 under which all the particles of one family belong to a single representation with the additional requirement that the masses of the third family fermions arise primarily from their unique yukawa coupling to a single superfield representation which predominantly contains the electroweak higgs superfields . it should be noted that exact yu in the cmssm leads to unacceptable values for the bottom - quark mass @xmath71 and , thus , must be corrected in order to become consistent with observations .
we will ignore this problem for the moment , but we will return to it in sec .
[ sec : quasi ] .
now , we assume that our effective theory below @xmath36 is the cmssm with yu .
this theory depends on the following parameters ( @xmath72 ) : @xmath73 where @xmath74 with @xmath75 being the gut gauge coupling constant and @xmath76 is the ratio of the two electroweak vevs .
the parameters @xmath77 and @xmath36 are evaluated consistently with the experimental values of the electromagnetic and strong fine - structure constants @xmath78 and @xmath79 , and the sine - squared of the weinberg angle @xmath80 at @xmath81 . to this end
, we integrate @xcite numerically the renormalization group equations ( rges ) for the mssm at two loops in the gauge and yukawa coupling constants from @xmath36 down to a common but variable @xcite susy threshold @xmath82 ( @xmath83 are the stop - quark mass eigenstates ) . from @xmath84 to @xmath81 ,
the rges of the non - susy sm are used .
the set of rges needed for our computation can be found in many references ( see e.g. ref .
we take @xmath85 which , as it turns out , leads to gauge coupling unification at @xmath36 with an accuracy better than 0.1@xmath86 .
so , we can assume exact unification once the appropriate susy particle thresholds are taken into account .
the unified third generation yukawa coupling constant @xmath87 at @xmath36 and the value of @xmath88 at @xmath84 are estimated using the experimental inputs for the top - quark mass @xmath89 and the @xmath90-lepton mass @xmath91 .
our integration procedure of the rges relies @xcite on iterative runs of these equations from @xmath36 to low energies and back for every set of values of the input parameters until agreement with the experimental data is achieved .
the susy corrections to @xmath92 are taken from ref .
@xcite and incorporated at @xmath84 .
assuming radiative electroweak symmetry breaking , we can express the values of the parameters @xmath93 ( up to its sign ) and @xmath94 ( or , equivalently , the mass @xmath95 of the @xmath18-odd neutral higgs boson @xmath96 ) at @xmath84 in terms of the other input parameters by minimizing the tree - level renormalization group ( rg ) improved potential @xcite at @xmath84 .
the resulting conditions are @xmath97 where @xmath98 , @xmath99 are the soft susy - breaking scalar higgs masses .
we can improve the accuracy of these conditions by including the full one - loop radiative corrections to the potential from ref .
@xcite at @xmath84 .
we find that the corrections to @xmath93 and @xmath95 from the full one - loop effective potential are minimized @xcite by our choice of @xmath84 .
so , a much better accuracy is achieved by using this variable susy threshold rather than a fixed one .
furthermore , we include in our calculation the two - loop radiative corrections to the masses @xmath100 and @xmath101 of the @xmath18-even neutral higgs bosons @xmath7 and @xmath102 .
these corrections are particularly important for the mass of the lightest @xmath18-even neutral higgs boson @xmath7 . finally , the susy corrections to @xmath71 are also included at @xmath84 using the relevant formulae of ref .
@xcite . as
already mentioned , the predicted value of the bottom - quark mass is not compatible with experiment .
however , we will ignore this problem for the moment .
the sign of @xmath93 is taken to be positive , since the @xmath103 case is excluded because it leads @xcite to a neutralino relic abundance which is well above unity , thereby overclosing the universe , for all @xmath95 s permitted by @xmath104 .
we are left with @xmath105 and @xmath62 as free input parameters .
the lsp is the lightest neutralino @xmath106 .
the mass matrix for the four neutralinos is @xmath107 in the @xmath108 basis . here , @xmath109 , @xmath110 , and @xmath111 , @xmath112 are the mass parameters of @xmath53 , @xmath113 in eq .
( [ gaugino ] ) . in cmssm
, the lightest neutralino turns out to be an almost pure bino @xmath53 .
the lsps are stable due to the presence of the unbroken r - parity , but can annihilate in pairs since this symmetry is a discrete @xmath56 symmetry .
this reduces their relic abundance in the universe .
if there exist sparticles with masses close to the mass of the lsp , their coannihilation @xcite with the lsp leads to a further reduction of the lsp relic abundance .
it should be noted that the number density of these sparticles is not boltzmann suppressed relative to the lsp number density .
they eventually decay yielding an equal number of lsps and , thus , contributing to the relic abundance of the lsps .
of particular importance is the next - to - lsp ( nlsp ) , which , in cmssm , turns out to be the lightest stau mass eigenstate @xmath114 .
its mass is obtained by diagonalizing the stau mass - squared matrix @xmath115 in the gauge basis ( @xmath116 ) . here ,
@xmath117}}$ ] is the soft susy - breaking mass of the left [ right ] handed stau @xmath118}$ ] and @xmath119 the tau - lepton mass .
the stau mass eigenstates are @xmath120 where @xmath121 is the @xmath122 mixing angle .
the large values of @xmath0 and @xmath90 yukawa coupling constants implied by yu cause soft susy - breaking masses of the third generation squarks and sleptons to run ( at low energies ) to lower physical values than the corresponding masses of the first and second generation .
furthermore , the large values of @xmath88 implied by yu lead to large off - diagonal mixings in the sbottom and stau mass - squared matrices .
these effects reduce further the physical mass of the lightest stau , which is the nlsp .
another effect of the large values of the @xmath0 and @xmath90 yukawa coupling constants is the reduction of the mass @xmath95 of the @xmath18-odd neutral higgs boson @xmath96 and , consequently , the other higgs boson masses to smaller values .
we now turn to the calculation of the cosmological relic abundance of the lightest neutralino @xmath31 ( almost pure @xmath53 ) in the cmssm with yu according to the standard cosmological scenario ( for non - standard scenaria , see e.g. ref .
@xcite ) . in general
, all sparticles contribute to @xmath123 , since they eventually turn into lsps , and all the ( co)annihilation processes must be considered .
the most important contributions , however , come from the lsp and the nlsp .
so , in the case of the cmssm , we should concentrate on @xmath31 ( lsp ) and @xmath124 ( nlsp ) and consider the coannihilation of @xmath31 with @xmath114 and @xmath125 .
the important role of the coannihilation of the lsp with sparticles carrying masses close to its mass in the calculation of the lsp relic abundance has been pointed out by many authors ( see e.g. refs .
@xcite ) . here , we will use the method of ref .
@xcite , which was also used in ref .
note that our analysis can be readily applied to any mssm scheme where the lsp and nlsp are the bino and stau respectively .
in particular , it applies to the cmssm without yu , where we have @xmath88 as an extra free input parameter .
the relevant quantity , in our case , is the total number density @xmath126 since the @xmath114 s and @xmath125 s decay into @xmath31 s after freeze - out . at cosmic temperatures relevant for freeze - out ,
the scattering rates of these ( non - relativistic ) sparticles off particles in the thermal bath are much faster than their annihilation rates since the ( relativistic ) particles in the bath are considerably more abundant .
consequently , the number densities @xmath127 ( @xmath128 , @xmath114 , @xmath125 ) are proportional to their equilibrium values @xmath129 to a good approximation , i.e. @xmath130 . the boltzmann equation ( see e.g. ref .
@xcite ) is then written as @xmath131 where @xmath102 is the hubble parameter , @xmath132 is the `` relative velocity '' of the annihilating particles , @xmath133 denotes thermal averaging and @xmath134 is the effective cross section defined by @xmath135 with @xmath136 being the total cross section for particle @xmath5 to annihilate with particle @xmath137 averaged over initial spin states . in our case , @xmath134 takes the following form @xmath138 for @xmath139 , we use the non - relativistic approximation @xmath140 @xmath141 here @xmath142 , 1 , 1 ( @xmath128 , @xmath114 , @xmath125 ) is the number of degrees of freedom of the @xmath5-th particle with mass @xmath143 and @xmath144 with @xmath145 being the photon temperature .    using boltzmann equation ( which is depicted in eq .
( [ boltzmann ] ) ) , we can calculate the relic abundance of the lsp at the present cosmic time .
it has been found @xcite to be given by @xmath146 with @xmath147 here @xmath148 is the planck scale , @xmath149 is the effective number of massless degrees of freedom at freeze - out @xcite and @xmath150 with @xmath151 being the freeze - out photon temperature calculated by solving iteratively the equation @xcite @xmath152 the constant @xmath153 is chosen to be equal to @xmath154 @xcite .
the freeze - out temperatures which we obtain here are of the order of @xmath155 and , thus , our non - relativistic approximation ( see eq .
( [ ri ] ) ) is _ a posteriori _ justified .
away from s - channel poles and final - state thresholds , the quantities @xmath156 are well approximated by applying the non - relativistic taylor expansion up to second order in the relative velocity @xmath132 : @xmath157 actually , this corresponds @xcite to an expansion in s and p waves .
the thermally averaged cross sections are then easily calculated @xmath158 here , we approximated the masses of the incoming particles by the neutralino mass , i.e. @xmath159 . the reduced mass of the incoming particles is then equal to @xmath160 .
we also thermally averaged over the relative velocity rather than the separate velocities of the incoming particles , which would be more accurate .
using eqs .
( [ sigmaeff1 ] ) , ( [ sigmaeff2 ] ) , ( [ sigmaeff3 ] ) , and ( [ average ] ) , one obtains @xmath161 where we sum over @xmath162 , @xmath163 , and @xmath164 with @xmath165 , @xmath166 , and @xmath167 , @xmath168 given by @xmath169 here @xmath170 , 4 , 2 for @xmath171 , @xmath163 , and @xmath172 respectively .
.feynman diagrams [ cols="^,^,^ " , ]     [ tab : charges ]    a moderate violation of exact yu can be naturally accommodated in this model by adding a new higgs superfield @xmath173 with yukawa couplings @xmath174 .
actually , ( * 15,2,2 * ) is the only representation of @xmath175 , besides ( * 1,2,2 * ) , which possesses such couplings to the matter superfields . in order to give superheavy masses to the color non - singlet components of @xmath173
, we need to include one more higgs superfield @xmath176 with the superpotential coupling @xmath177 , whose coefficient is of the order of @xmath36 .
after the breaking of @xmath175 to @xmath34 , the two color singlet @xmath178 doublets @xmath179 , @xmath180 contained in @xmath173 can mix with the corresponding doublets @xmath181 , @xmath182 in @xmath7 .
this is mainly due to the terms @xmath177 and @xmath183 .
actually , since @xmath184 there are two independent couplings of the type @xmath183 ( both suppressed by the string scale @xmath185 , as they are non - renormalizable ) .
one of these couplings is between the @xmath186 singlets in @xmath187 and @xmath188 and the other between the @xmath186 triplets in these combinations .
so , we obtain two bilinear terms @xmath189 and @xmath190 with different coefficients , which are suppressed by @xmath191 . these terms together with the terms @xmath192 and @xmath193 from @xmath177 , which have equal coefficients , generate different mixings between @xmath181 , @xmath179 and @xmath182 , @xmath180 .
consequently , the resulting electroweak doublets @xmath44 , @xmath45 contain @xmath194 violating components suppressed by @xmath191 and fail to form a @xmath186 doublet by an equally suppressed amount .
so , yu is naturally and moderately violated .
unfortunately , as it turns out , this violation is not adequately large for correcting the bottom - quark mass within the framework of the cmssm with @xmath1 .    in order to allow for a more sizable violation of yu
, we further extend the model by including the superfield @xmath195 with the coupling @xmath196 . to give superheavy masses to the color non - singlets in @xmath195
, we introduce one more superfield @xmath197 with the coupling @xmath198 , whose coefficient is of order @xmath36 .
the superpotential terms @xmath198 and @xmath199 imply that , after the breaking of @xmath175 to @xmath34 , @xmath195 acquires a vev of order @xmath36 .
the coupling @xmath196 then generates @xmath186 violating unsuppressed bilinear terms between the doublets in @xmath176 and @xmath7 .
these terms can overshadow the corresponding ones from the non - renormalizable term @xmath183 .
the resulting @xmath186 violating mixing of the doublets in @xmath7 and @xmath173 is then unsuppressed and we can obtain stronger violation of yu .      to further analyze the mixing of the doublets in @xmath7 and @xmath173 , observe that the part of the superpotential corresponding to the symbolic couplings @xmath177 , @xmath196 is properly written as @xmath200 where @xmath201 is a mass parameter of order @xmath36 , @xmath202 is a dimensionless parameter of order unity , @xmath203 denotes trace taken with respect to the @xmath204 and @xmath178 indices , and the superscript t denotes the transpose of a matrix .
after the breaking of @xmath175 to @xmath34 , @xmath195 acquires a vev @xmath205 . substituting it by this vev in the above couplings , we obtain @xmath206 where the ellipsis in eq .
( [ mass ] ) contains the colored components of @xmath176 , @xmath173 and @xmath207 . inserting eqs .
( [ mass ] ) and ( [ triplet ] ) into eq .
( [ expmix ] ) , we obtain @xmath208 so , we get two pairs of superheavy doublets with mass @xmath201 . they are predominantly given by @xmath209 the orthogonal combinations of @xmath181 , @xmath210 and @xmath182 , @xmath211 constitute the electroweak doublets @xmath212 the superheavy doublets in eq .
( [ superdoublets ] ) must have vanishing vevs , which readily implies that @xmath213 and @xmath214 .
equation  ( [ ew ] ) then gives @xmath215 , @xmath216 . from the third generation
yukawa couplings @xmath217 , @xmath218 , we obtain @xmath219 where @xmath220 . from eqs .
( [ top ] ) and ( [ bottomtau ] ) , we see that yu is now replaced by the yquc @xmath221 for simplicity , we restricted ourselves here to real values of @xmath153 only which lie between zero and unity , although @xmath153 is , in general , an arbitrary complex quantity with @xmath222 .      below the gut scale @xmath36 , the particle content of our model reduces to this of mssm ( modulo sm singlets ) . we assume universal soft susy breaking scalar masses @xmath61 , gaugino masses @xmath63 , and trilinear scalar couplings @xmath62 at @xmath36 . therefore , the resulting mssm is the so - called cmssm @xcite with @xmath1 supplemented by the yquc in eq .
( [ minimal ] ) . with these initial conditions ,
we run the mssm rges @xcite between @xmath36 and a common variable susy threshold @xmath84 ( see sec .  [
sec : mssm ] ) determined in consistency with the susy spectrum of the model . at @xmath84 ,
we impose radiative electroweak symmetry breaking , evaluate the susy spectrum and incorporate the susy corrections @xcite to the @xmath0-quark and @xmath90-lepton masses .
note that the corrections to the @xmath90-lepton mass ( almost 4@xmath86 ) lead @xcite to a small reduction of @xmath88 . from @xmath84 to @xmath81 , the running of gauge and yukawa coupling constants is continued using the sm rges .    for presentation purposes , @xmath63 and @xmath61
can be replaced @xcite by the lsp mass @xmath223 and the relative mass splitting between this particle and the lightest stau @xmath224 ( recall that @xmath124 is the nlsp in this case ) . for simplicity , we restrict this presentation to the @xmath225 case ( for @xmath226 see refs .
so , our input parameters are @xmath223 , @xmath227 , @xmath153 , and @xmath88 .    for
any given @xmath228 in the range in eq .
( [ mbrg ] ) and with fixed @xmath89 and @xmath229 , we can determine the parameters @xmath153 and @xmath88 at @xmath84 so that the yquc in eq .
( [ minimal ] ) is satisfied .
we are , thus , left with @xmath223 and @xmath227 as free parameters .
restrictions on the parameters of our model can be derived by imposing a number of cosmological and phenomenological requirements ( for similar recent analyses , see refs .
these constraints result from    @xmath230 _ cdm considerations_. as discussed in sec .
[ sec : mssm ] , in the context of the cmssm , the lsp can be the lightest neutralino which is an almost pure bino .
it naturally arises @xcite as a cdm candidate .
we require its relic abundance , @xmath231 , not to exceed the @xmath11 c.l .
upper bound on the cdm abundance derived @xcite by wmap : @xmath232 we calculate @xmath231 using micromegas @xcite , which is certainly one of the most complete publicly available codes . among other things ,
it includes all possible coannihilation processes @xcite and one - loop qcd corrections @xcite to the higgs decay widths and couplings to fermions .
@xmath230 _ branching ratio of @xmath104_. taking into account the experimental results of ref .
@xcite on this ratio , @xmath233 , and combining @xcite appropriately the experimental and theoretical errors involved , we obtain the @xmath11 c.l .
range @xmath234 although there exist more recent experimental data @xcite on the branching ratio of @xmath104 , we do not use them here .
the reason is that these data do not separate the theoretical errors from the experimental ones and , thus , the derivation of the @xmath11 c.l .
range is quite ambiguous . in any case , the @xmath11 c.l .
limits obtained in ref .
@xcite on the basis of these latest measurements are not terribly different from the ones quoted in eq .
( [ bsgb ] ) . in view of this and
the fact that , in our case , the restrictions from @xmath233 are overshadowed by other constraints ( see sec .  [
sec : parameters ] ) , we limit ourselves to the older data .
we compute @xmath233 by using an updated version of the relevant calculation contained in the micromegas package @xcite . in this code , the sm contribution is calculated following ref .
the charged higgs ( @xmath235 ) contribution is evaluated by including the next - to - leading order ( nlo ) qcd corrections @xcite and @xmath88 enhanced contributions @xcite .
the dominant susy contribution includes resummed
nlo susy qcd corrections @xcite , which hold for large @xmath88 .
@xmath230 _ muon anomalous magnetic moment_. the deviation , @xmath236 , of the measured value of @xmath237 from its predicted value in the sm , @xmath238 , can be attributed to susy contributions , which are calculated by using the micromegas routine @xcite .
the calculation of @xmath238 is not yet stabilized mainly because of the instability of the hadronic vacuum polarization contribution . according to recent calculations ( see e.g. refs .
@xcite ) , there is still a considerable discrepancy between the findings based on the @xmath239 annihilation data and the ones based on the @xmath90-decay data .
taking into account the results of ref .
@xcite and the experimental measurement of @xmath237 reported in ref .
@xcite , we get the following @xmath11 c.l .
ranges : @xmath240 following the common practice @xcite , we adopt the restrictions to parameters induced by eq .
( [ g2e ] ) , since eq .
( [ g2 t ] ) is considered as quite oracular , due to poor @xmath90-decay data .
it is true that there exist more recent experimental data @xcite on @xmath237 than the ones we considered and more updated estimates of @xmath236 than the one in ref .
@xcite ( see e.g. ref .
@xcite ) . however , only the @xmath11 c.l .
upper limit on @xmath236 enters into our analysis here and its new values are not very different from the one in eq .
( [ g2e ] ) .
@xmath230 _ collider bounds_. here , as it turns out , the only relevant collider bound is the @xmath11 c.l .
lep lower bound @xcite on the mass of the lightest @xmath18-even neutral higgs boson @xmath7 : @xmath241 the susy corrections to the lightest @xmath18-even higgs boson mass @xmath100 are calculated at two loops by using the feynhiggsfast program @xcite included in the micromegas code @xcite .
we will now try to delineate the parameter space of our model with @xmath1 which is consistent with the constraints in sec .
[ sec : pheno ] .
the restrictions on the @xmath242 plane , for @xmath225 and the central values of @xmath243 and @xmath244 , are indicated in fig .
[ figa ] by solid lines , while the upper bound on @xmath223 from eq .
( [ cdmb ] ) , for @xmath245~{\rm gev}$ ] , is depicted by a dashed [ dotted ] line .
we observe the following :    * the lower bounds on @xmath223 are not so sensitive to the variations of @xmath228 . *
the lower bound on @xmath223 from eq .
( [ mhb ] ) overshadows all the other lower bounds on this mass . * the upper bound on @xmath223 from eq .
( [ cdmb ] ) is very sensitive to the variations of @xmath228 . in particular
, one notices the extreme sensitivity of the almost vertical part of the corresponding line , where the lsp annihilation via an @xmath96-boson exchange in the s - channel is @xcite by far the dominant process , since @xmath95 , which is smaller than @xmath246 , is always very close to it as seen from fig .
this sensitivity can be understood from fig .
[ figc ] , where @xmath95 is depicted versus @xmath223 for various @xmath228 s .
we see that , as @xmath228 decreases , @xmath95 increases and approaches @xmath246 .
the @xmath96-pole annihilation is then enhanced and @xmath231 is drastically reduced causing an increase of the upper bound on @xmath223 . * for low @xmath227 s , bino - stau coannihilations @xcite
take over leading to a very pronounced reduction of the lsp relic abundance @xmath231 , thereby enhancing the upper limit on @xmath223 .
so , we obtain the almost horizontal tail of the allowed region in fig .  [ figa ] .    for @xmath1 , @xmath225 , @xmath243 and @xmath244
, we find the following allowed ranges of parameters : @xmath247 the splitting between the bottom ( or tau ) and top yukawa coupling constants @xmath248 ranges between 0.25 and 0.29 .
it is interesting to note that our susy gut model gives rise @xcite naturally to a modified version of shifted hybrid inflation @xcite .
hybrid inflation @xcite , which is certainly one of the most promising inflationary scenarios , uses two real scalars : one which provides the vacuum energy density for driving inflation and a second which is the slowly varying field during inflation .
this scheme , which is naturally incorporated @xcite in susy guts ( for an updated review , see ref .
@xcite ) , in its standard realization has the following property @xcite : if the gut gauge symmetry breaking predicts topological defects such as magnetic monopoles @xcite , cosmic strings @xcite , or domain walls @xcite , these defects are copiously produced at the end of inflation . in the case of monopoles or walls , this leads to a cosmological catastrophe @xcite .
the breaking of the @xmath175 symmetry predicts the existence of doubly charged monopoles @xcite .
so , any ps susy gut model incorporating the standard realization of susy hybrid inflation would be unacceptable . one way to remedy
this is to invoke @xcite thermal inflation @xcite to dilute the primordial monopoles well after their production .
alternatively , we can construct variants of the standard susy hybrid inflationary scenario such as smooth @xcite or shifted @xcite hybrid inflation which do not suffer from the monopole overproduction problem . in the latter scenario ,
we generate @xcite a shifted inflationary trajectory so that @xmath175 is already broken during inflation .
this could be achieved @xcite in our susy gut model even before the introduction of the extra higgs superfields , but only by utilizing non - renormalizable terms .
the inclusion of @xmath249 and @xmath250 does not change this situation .
the inclusion of @xmath195 and @xmath197 , however , very naturally gives rise @xcite to a shifted path , but now with renormalizable interactions alone .
we showed that particle physics provides us with a number of candidate particles out of which the cdm of the universe can be made .
these particles are not invented solely for explaining the cdm , but they are naturally there in various particle physics models .
we discussed in some detail the major candidates which are the axion , the lightest neutralino , the axino , and the gravitino .
the last three particles exist only in susy theories and can be stable provided that they are the lsp .
the axion is a pseudo nambu - goldstone boson associated with the spontaneous breaking of a pq symmetry .
this is a global anomalous @xmath19 symmetry invoked to solve the strong @xmath18 problem .
it is , actually , the most natural solution to this problem which is available at present . the axions are extremely light particles and are generated at the qcd phase transition carrying zero momentum .
we argued that these particles can easily provide the cdm in the universe .
however , if the pq field emerges with non - zero value at the end of inflation , they lead to isocurvature perturbations , which , for superheavy inflationary scales , are too strong to be compatible with the recent results of the wmap satellite on the cmbr anisotropies
.    the most popular cdm candidate is , certainly , the lightest neutralino which is present in all susy models and can be the lsp for a wide range of parameters .
we considered it within the simplest susy framework which is the mssm whose salient properties were summarized .
we used exclusively the constrained version of mssm which is known as cmssm and is based on universal boundary conditions . in this case , the lightest neutralino is an almost pure bino , whereas the nlsp is the lightest stau .
we sketched the calculation of the neutralino relic abundance in the universe paying particular attention not only to the neutralino pair annihilations , but to the neutralino - stau coannihilations too .
it is very important for the accuracy of the calculation to treat poles and final - state thresholds properly and include the one - loop qcd corrections to the higgs boson decay widths and the fermion masses .
we find that two effects help us reduce the neutralino relic abundance and satisfy the wmap constraint on cdm : the resonantly enhanced neutralino pair annihilation via an @xmath96-pole exchange in the s - channel , which appears in the large @xmath88 regime , and the strong neutralino - stau coannihilation , which is achieved when these particles are almost degenerate in mass .
the axino , which is the susy partner of the axion , can also be the lsp in many cases since its mass is a strongly model - dependent parameter in the cmssm .
it is produced thermally by 2-body scattering or decay processes in the thermal bath , or non - thermally by the decay of sparticles which are already frozen out of thermal equilibrium . for small axino masses ,
tp is more important yielding a very narrow favored region in the parameter space . for larger axino masses ,
however , ntp is more efficient and the favored region in the parameter space becomes considerably wider .
one finds that , in the case of the cmssm , almost any point on the @xmath251 plane can be allowed by axino cdm considerations .
the required reheat temperatures though are quite small ( @xmath252 ) .
the mass of the gravitino is a practically free parameter in the cmssm .
so , the gravitino can easily be the lsp and , in principle , contribute to the cdm of the universe .
it is produced thermally by 2-body scattering processes in the thermal bath as well as non - thermally by the decay of the nlsp , which can be either the neutralino or the stau .
in contrast to the axino case , however , the nlsp can now have quite a long lifetime .
the electromagnetic showers resulting from the nlsp decay can destroy the successful predictions of bbn .
so , we obtain strong constraints which allow only very limited regions of the parameter space of the cmssm .
as it turns out , ntp in these regions is not efficient enough to account for cdm .
we can , however , make these regions cosmologically favored by raising @xmath253 to enhance tp of gravitinos .
we studied the cmssm with @xmath1 and @xmath225 applying a yquc which originates from a ps susy gut model .
this condition yields an adequate deviation from yu which allows an acceptable @xmath228 .
we , also , imposed the constraints from the cdm in the universe , @xmath104 , @xmath254 and @xmath100 .
we found that there exists a wide and natural range of cmssm parameters which is consistent with all the above constraints .
the parameter @xmath88 ranges between about 58 and 59 and the asymptotic splitting between the bottom ( or tau ) and the top yukawa coupling constants varies in the range @xmath255 for central values of @xmath228 and @xmath256
. the predicted lsp mass can be as low as about @xmath257 .
moreover , the model resolves the @xmath93 problem of mssm , predicts stable proton , generates the baryon asymmetry of the universe via primordial leptogenesis , and gives rise to a new version of shifted hybrid inflation which is based solely on renormalizable interactions .
we thank l. roszkowski , p. sikivie , and f.d .
steffen for useful suggestions . this work was supported by european union under the contract mrtn - ct-2004 - 503369 as well as the greek ministry of education and religion and the epeak program pythagoras .    # 1#2#3int .
mod . phys . * # 1 * ,  # 3  ( # 2 ) # 1#2#3phys .
b * # 1 * ,  # 3  ( # 2 ) # 1#2#3z .
c * # 1 * ,  # 3  ( # 2 ) # 1#2#3phys .
# 1 * ,  # 3  ( # 2 ) # 1#2#3rev .
# 1 * ,  # 3  ( # 2 ) # 1#2#3phys . rep . *
# 1 * ,  # 3  ( # 2 ) # 1#2#3phys . rev .
d * # 1 * ,  # 3  ( # 2 ) # 1#2#3nucl . phys . * b#1 * ,  # 3  ( # 2 ) # 1#2#3nucl .
b ( proc . sup . ) * # 1 * ,  # 3  ( # 2 ) # 1#2#3mod .
# 1 * ,  # 3  ( # 2 ) # 1#2#3mod .
lett . a * # 1 * ,  # 3  ( # 2 ) # 1#2#3annu .
nucl . part .
sci . * # 1 * ,  # 3  ( # 2 ) # 1#2#3sov .
# 1 * ,  # 3  ( # 2 ) # 1#2#3jetp lett . *
# 1 * ,  # 3  ( # 2 ) # 1#2#3acta phys .
# 1 * ,  # 3  ( # 2 ) # 1#2#3riv .
nuovo cim . *
# 1 * ,  # 3  ( # 2 ) # 1#2#3ann . phys . * # 1 * ,  # 3  ( # 2 ) # 1#2#3prog .
# 1 * ,  # 3  ( # 2 ) # 1#2#3astrophys .
j. lett . *
# 1 * ,  # 3  ( # 2 ) # 1#2#3nature * # 1 * ,  # 3  ( # 2 ) # 1#2#3astrophys .
j. * # 1 * ,  # 3  ( # 2 ) # 1#2#3astron .
j. * # 1 * ,  # 3  ( # 2 ) # 1#2#3mnras * # 1 * ,  # 3  ( # 2 ) # 1#2#3gen .
# 1 * ,  # 3  ( # 2 ) # 1#2#3science * # 1 * ,  # 3  ( # 2 ) # 1#2#3bull .
astron . soc . *
# 1 * ,  # 3  ( # 2 ) # 1#2#3_ibid .
_ * # 1 * ,  # 3  ( # 2 ) # 1#2#3comput .
. commun . *
# 1 * ,  # 3  ( # 2 ) # 1#2#3astropart . phys . *
# 1 * ,  # 3  ( # 2 ) # 1#2#3eur .
j. c * # 1 * ,  # 3  ( # 2 ) # 1#2#3nucl .
meth . a * # 1 * ,  # 3  ( # 2 ) # 1#2#3j
. high energy phys . * # 1 * ,  # 3  ( # 2 ) # 1#2#3j .
# 1 * ,  # 3  ( # 2 ) # 1#2#3astrophys .
* # 1 * ,  # 3  ( # 2 ) # 1#2#3lect .
notes phys . *
# 1 * ,  # 3  ( # 2 ) # 1#2#3j .
phys . a * # 1 * ,  # 3  ( # 2 ) # 1#2#3prog . part .
# 1 * ,  # 3  ( # 2 ) # 1#2#3jetp ( sov . phys . ) * # 1 * ,  # 2  ( # 3 )      c.l .
bennett et al : ; d.n .
spergel et al : w.l .
freedman et al : d.v .
ahluwalia - khalilova , d. grumiller : ; ; m. cirelli , n. fornengo , a. strumia : hep - ph/0512090 r. peccei , h. quinn : ; s. weinberg : ; f. wilczek : d.s.p .
dearborn , d.n .
schramm , g. steigman : j. preskill , m.b .
wise , f. wilczek : ; l.f .
abbott , p. sikivie : m. dine , w. fischler : p.j .
steinhardt , m.s .
turner : ; k. yamamoto : g. lazarides , c. panagiotakopoulos , q. shafi : ; g. lazarides , r.k .
schaefer , d. seckel , q. shafi : k. choi , e.j .
chun , j.e .
kim : ; k. dimopoulos , g. lazarides : a.h .
for a recent review see g. lazarides : p. sikivie : ; g. lazarides , q. shafi : ; h. georgi , m.b .
wise : m.s .
turner : l. campanelli , m. giannotti : astro - ph/0512458 h.v .
peiris et al : c. gordon , a. lewis : ; p. crotty , j. garca - bellido , j. lesgourgues , a. riazuelo : ; c. gordon , k.a .
malik : ; m. bucher , j. dunkley , p.g .
ferreira , k. moodley , c. skordis : ; k. moodley , m. bucher , j. dunkley , p.g .
ferreira , c. skordis : ; m. beltrn , j. garca - bellido , j. lesgourgues , a. riazuelo : k. dimopoulos , g. lazarides , d. lyth , r. ruiz de austri : ; g. lazarides , r. ruiz de austri , r. trotta : ; g. lazarides : g.l .
kane , c. kolda , l. roszkowski , j.d . wells : j.r .
ellis , k.a .
olive , y. santoso , v.c .
spanos : b. ananthanarayan , g. lazarides , q. shafi : ; . for
a more recent update see u. sarid : hep - ph/9610341 ( in : _ snowmass 1996 , new directions for high - energy physics _ )
gmez , g. lazarides , c. pallis : ; m. drees , m.m .
nojiri : v. barger , m. s. berger , p. ohmann : ; m. drees , m.m .
nojiri : ; m. olechowski , s. pokorski : d. pierce , j. bagger , k. matchev , r. zhang : g. gamberini , g. ridolfi , f. zwirner : h. baer , c. chen , m. drees , f. paige , x. tata : f.m .
borzumati , m. olechowski , s. pokorski : h. baer , m. brhlik , d. castao , x. tata : k. griest , d. seckel : c. pallis : ; ; hep - ph/0510234 m. drees , m.m .
nojiri : ; m. drees : hep - ph/9703260 j. ellis , t. falk , k.a .
olive : ; j. ellis , t. falk , g. ganis , k.a .
olive , m. schmitt : j. ellis , t. falk , k.a .
olive , m. srednicki : , ( e ) e.w .
kolb , m.s .
turner : _ the early universe _ ( addison - wesley , redwood city ca 1990 )    k. griest , m. kamionkowski , m.s .
turner : p. gondolo , g. gelmini : t. falk , k.a .
olive , m. srednicki : ; j. edsj , p. gondolo : t. nihei , l. roszkowski , r. ruiz de austri : m.e .
gmez , g. lazarides , c. pallis : t. nihei , l. roszkowski , r. ruiz de austri : v.d .
barger , c. kao : j.r .
ellis , t. falk , g. ganis , k.a .
olive , m. srednicki : phys .
b * 510 * , 236 ( 2001 ) ; l. roszkowski , r. ruiz de austri , t. nihei : ; a.b .
lahanas , d.v .
nanopoulos , v.c .
spanos : hep - ph/0211286 ( in : _ oulu 2002 , beyond the desert _ ) a. djouadi , j. kalinowski , m. spira : g. blanger , f. boudjema , a. pukhov , a. semenov : p. gondolo , j. edsj , l. bergstrm , p. ullio , e.a .
baltz : astro - ph/0012234 ( in : _ york 2000 , the identification of dark matter _ ) l. covi , j.e .
kim , l. roszkowski : ; l. covi , h .- b .
kim , j.e .
kim , l. roszkowski : ; a. brandenburg , f.d .
steffen : l. roszkowski : hep - ph/0102325 ( in : _ york 2000 , the identification of dark matter _ ) ; j.e .
kim : astro - ph/0205146 ( in : _ cape town 2002 , dark matter in astro- and particle physics _ ) k. rajagopal , m.s .
turner , f. wilczek : ; e.j .
chun , j.e .
kim , h.p .
nilles : l. covi , l. roszkowski , r. ruiz de austri , m. small : ; a. brandenburg , l. covi , k. hamaguchi , l. roszkowski , f.d .
steffen : l. covi , l. roszkowski , m. small : j. ellis , j.e .
kim , d. nanopoulos : ; t. moroi , h. murayama , m. yamaguchi : m. bolz , w. buchmller , m. plmacher : e. holtmann , m. kawasaki , k. kohri , t. moroi : phys .
d * 60 * , 023506 ( 1999 ) ; t. gherghetta , g.f .
giudice , a. riotto : ; t. asaka , k. hamaguchi , k. suzuki : ; j.l .
feng , a. rajaraman , f. takayama : ; ; j.l .
feng , s. su , f. takayama : ; l. roszkowski , r. ruiz de austri , k .- y .
choi : j. high energy phys . * 08 * , 080 ( 2005 ) ; d.g .
cerdeno , k .- y .
choi , k. jedamzik , l. roszkowski , r. ruiz de austri : hep - ph/0509275 m. bolz , a. brandenburg , w. buchmller : f.d .
steffen : hep - ph/0507003 g. lazarides , c. pallis : hep - ph/0404266 ( in : _ vrnjacka banja 2003 , mathematical , theoretical and phenomenological challenges beyond the standard model _ ) g. lazarides , c. pallis : hep - ph/0406081 g. lazarides , c. panagiotakopoulos : ; s. khalil , g. lazarides , c. pallis : l. hall , r. rattazzi , u. sarid : ; m. carena , m. olechowski , s. pokorski , c.e.m .
wagner : m. carena , d. garcia , u. nierste , c.e.m .
wagner : s.f .
king , m. oliveira : s. abel et al .
( sugra working group collaboration ) : hep - ph/0003154 h. baer , j. ferrandis , k. melnikov , x. tata : t. blaek , r. dermek , s. raby : ; d. auto et al : u. chattopadhyay , a. corsetti , p. nath : ; c. pallis : g. lazarides , q. shafi , c. wetterich : ; g. lazarides , q. shafi : m.e .
gmez , g. lazarides , c. pallis : ; c. pallis , m.e .
gmez : hep - ph/0303098 r. barate et al .
( aleph collaboration ) : ; k. abe et al .
( belle collaboration ) : ; s. chen et al .
( cleo collaboration ) : r. jeannerot , s. khalil , g. lazarides : r. jeannerot , s. khalil , g. lazarides , q. shafi : g. lazarides : hep - ph/0011130 ( in : _ cascais 2000 , recent developments in particle physics and cosmology _ ) ; r. jeannerot , s. khalil , g. lazarides : hep - ph/0106035 ( in : _
cairo 2001 , high energy physics _ ) i. antoniadis , g.k .
leontaris : g. lazarides , q. shafi : g.r .
dvali , g. lazarides , q. shafi : g. lazarides , c. panagiotakopoulos , q. shafi : n. ganoulis , g. lazarides , q. shafi : ; g. lazarides , q. shafi : g. lazarides , q. shafi : m. fukugita , t. yanagida : g. lazarides : ; g. lazarides , n.d .
vlachos : ; t. dent , g. lazarides , r. ruiz de austri : ; m.e .
gmez , c. pallis : hep - ph/0303094 ( in : _
hamburg 2002 , supersymmetry and unification of fundamental interactions _ ) j. ellis , k.a .
olive , y. santoso , v.c .
spanos : ; a.b .
lahanas , d.v .
nanopoulos : ; h. baer , c. balzs : ; u. chattopadhyay , a. corsetti , p. nath : h. goldberg : ; j.r .
ellis , j.s .
hagelin , d.v .
nanopoulos , k.a .
olive , m. srednicki : b. aubert et al .
( babar collaboration ) : hep - ex/0207074 ; hep - ex/0207076 j.r .
ellis , s. heinemeyer , k.a .
olive , g. weiglein : a.l .
kagan , m. neubert : ; p. gambino , m. misiak : m. ciuchini , g. degrassi , p. gambino , g. giudice : ; g. degrassi , p. gambino , g.f .
giudice : s. martin , j. wells : m. davier , hep - ex/0312065 ( in : _
pisa 2003 , sighad 03 _ )
a. hcker : hep - ph/0410081 ( in : _ beijing 2004 , ichep 2004 _ ) ; a. vainshtein : g.w . bennett et al .
( muon @xmath258 - 2 collaboration ) : phys . rev . lett . * 89 * , 101804 ( 2002 ) , ( e ) g.w .
bennett et al .
( muon @xmath258 - 2 collaboration ) : phys . rev . lett . * 92 * , 161802 ( 2004 ) ; a. aloisio et al .
( kloe collaboration ) : aleph , delphi , l3 and opal collaborations , the lep higgs working group for higgs boson searches : hep - ex/0107029 ( in : _
budapest 2001 , high energy physics _ ) ; lhwg - note/2002 - 01 , + http://lephiggs.web.cern.ch/lephiggs/papers/july2002_sm/index.html s. heinemeyer , w. hollik , g. weiglein : hep - ph/0002213 a.b .
lahanas , d.v .
nanopoulos , v.c .
spanos : a.d .
linde : e.j .
copeland , a.r .
liddle , d.h .
lyth , e.d .
stewart , d. wands : ; g.r .
dvali , q. shafi , r.k .
schaefer : ; g. lazarides , r.k .
schaefer , q. shafi : v.n .
senoguz , q. shafi : ; g. lazarides , c. panagiotakopoulos : g. t hooft : ; a. polyakov : jetp lett . * 20 * , 194 ( 1974 ) t.w.b .
kibble , g. lazarides , q. shafi : ya.b .
zeldovich , i.yu .
kobzarev , l.b .
okun : t.w.b .
kibble : g. lazarides , m. magg , q. shafi : g. lazarides , q. shafi : g. lazarides , c. panagiotakopoulos , q. shafi : ; d.h .
lyth , e.d .


the use of ultraviolet emission to shed light on the evolutionary history of galaxies is not straightforward .
the presence of dust in galaxies represents one of the major obstacles complicating a direct quantification of the star formation activity in local and high redshift galaxies .
absorption by dust grains reddens the spectra at short wavelengths completely modifying the spectral energy distribution of galaxies .
since the uv radiation is preferentially emitted by young stars that are generally more affected by attenuation from surrounding dust clouds than older stellar populations , rest - frame uv observations can lead to incomplete and/or biased reconstructions of the star formation activity and star formation history of galaxies affected by dust absorption , unless proper corrections are applied .
+ radiative transfer models suggest that the total - ir ( tir ) to uv luminosity ratio method ( i.e. @xcite ) is the most reliable estimator of the dust attenuation in star - forming galaxies because it is almost completely independent of the extinction mechanisms ( i.e. dust / star geometry , extinction law , see @xcite ) .
this method is based on the assumption that a fraction of photons emitted by newly formed young stars are absorbed by the dust .
the dust heats up and subsequently re - emits the energy in the mid- and far - infrared .
the amount of uv attenuation can thus be quantified by means of an energy balance . the reliability of this method has also made the @xmath1 ratio ideal to calibrate empirical methods to correct for dust attenuation when far - infrared observations are not available , like the well known @xmath2 relation ( @xcite , where @xmath3 is the slope of the uv continuum spectrum in the range 1300@xmath42600  ) . + although almost independent from dust - geometry and extinction law , the @xmath1 ratio method unfortunately depends on the age of the underlying stellar populations . in systems with low specific star formation rate ( ssfr )
the dust heating by old stars becomes important , and only a fraction of the dust emission is related to the ultraviolet absorption @xcite .
+ in the past this _ age effect _ has been generally considered negligible since uv observations were available only for active star - forming systems for which the @xmath1 vs. @xmath0 can be assumed to be independent from the star formation history ( sfh , @xcite ) .
however after the launch of the galaxy evolution explorer ( galex , @xcite ) this may not be the case any more .
galex is delivering to the community an unprecedented amount of uv data on local and high redshift galaxies covering the whole range of morphologies and star formation histories : from starbursts , to passive spirals and elliptical galaxies ( e.g. @xcite ) .
can we extend the use of the _ standard method _ ( calibrated on active star forming galaxies and starbursts ) to correct for dust attenuation to the thousands of galaxies detected by galex ? in recent years various studies have been undertaken to address these issues ( i.e. @xcite ) .
however more attention has been always given to the @xmath1 vs. @xmath3 relation than to the @xmath1 vs. @xmath0 relation .
@xcite suggested that quiescent galaxies deviate from the @xmath2 relation of starburst galaxies , because they tend to have redder ultraviolet spectra at fixed @xmath1 ratio .
they interpreted the different behavior of starbursts and normal galaxies as due to a difference in the star formation histories , proposing that the offset from the starburst @xmath2 relation can be predicted using the birthrate parameter @xmath5 ( @xcite ; i.e. the ratio of the current to the mean past star formation activity ) .    in this work
we use a different approach , investigating the dependence of the @xmath1 vs. @xmath0 relation on the galaxy star formation history , to quantify the influence of the age - independent correction on the interpretation of uv observations .
we will adopt a very simple dust geometry model and extinction law to show that the use of standard ( i.e. age independent ) methods can strongly bias our interpretation of uv properties of local and high redshift galaxies .
the main goal of this paper is to provide the community with some new empirical relations based on observable quantities suitable for deriving dust attenuation corrections taking into account the _ age effect_. more detailed geometries and dust models have been developed in the last decade ( e.g. @xcite ) .
however , since our results do not depend on the geometry or extinction law adopted ( as shown in appendix a ) , we decided to adopt a very simple approach in order to make our recipes useful for the widest possible range of astrophysical applications .
in section 2 we present our model and discuss the variation of the @xmath6 vs. @xmath7 relation with the age of the underlying stellar population . in section 3
we apply our model to a sample of nearby field and cluster galaxies and in section 4 we use this sample to test new empirical recipes to determine @xmath7 .
additional applications of these techniques are discussed in section 5 .
in order to investigate the correlation between the uv dust attenuation @xmath0 and the @xmath1 ratio for different stellar populations we used a library of spectral energy distributions ( seds ) obtained using the bruzual & charlot population synthesis models @xcite .
we adopted a @xcite initial mass function ( imf ) and a star formation history ( sfh ) a la sandage in the formalism used by @xcite : @xmath8 where @xmath9 is the star formation rate per unit mass , @xmath10 is the age of the galaxy ( we assumed @xmath10=13 gyr at the present epoch ) and @xmath11 is the time at which the star formation rate reaches the highest value over the whole galaxy history : short @xmath11 correspond to galaxies dominated by old stellar populations while long @xmath11 correspond to young ( i.e. star forming ) galaxies .
we investigated a range of @xmath11 between 0.1 and 25 gyr , with steps of 0.2 gyr for 0.1@xmath1210 gyr and 0.5 gyr for @xmath1310 gyr and metallicities in the range 0.02@xmath142.5 z@xmath15 in five steps : 0.02 , 0.2 , 0.4 , 1 , and 2.5 z@xmath15 .
the library obtained is able to reproduce all the seds typically observed in local galaxies @xcite and it is only a function of @xmath11 will be only used to indicate the shape of the sed and it must not be adopted to quantify the real age of the stellar populations . ] .
similar libraries can be obtained assuming different sfh ( e.g. exponential sfh , see appendix a ) .
+ each synthetic sed in our library has been reddened using different values of @xmath7 in the range 0@xmath16 4 mag , the typical range observed in normal star - forming galaxies ( e.g. @xcite ) .
higher uv attenuations are normally associated with strong starbursts and highly obscured objects which are outside the goal of the present work . for each @xmath7 ,
@xmath17 ( the attenuation at each @xmath18 ) has been derived assuming a large magellanic cloud ( lmc ) extinction curve @xcite : @xmath19 and a simple sandwich model for dust geometry @xcite , where a thin layer of dust of thickness @xmath20 is embedded in a thick layer of stars : @xmath21\left(1+{\rm e}^{-\tau_{dust}(\lambda ) \cdot { \rm sec}(i)}\right ) + & & { } \nonumber \\ { } + \left[\frac{\zeta(\lambda)}{\tau_{dust}(\lambda ) \cdot { \rm sec}(i)}\right ] \cdot \left(1-{\rm e}^{-\tau_{dust}(\lambda ) \cdot { \rm sec}(i)}\right ) \bigg )    ~~[{\rm mag}]\end{aligned}\ ] ] where @xmath22 is the galaxy inclination , @xmath23 is the optical depth and the dust to stars scale height ratio @xmath24 depends on @xmath18 ( in units of  ) as @xcite : @xmath25 in the case of the fuv band ( @xmath26  ) @xmath27 . in this case
@xmath28(fuv ) can be derived by inverting eq .
[ alambda ] : @xmath29 \cdot \big(0.0259 + 1.2002 \times a({\rm fuv } )   + 1.5543 \times   & & { } \nonumber \\        { }   \times a({\rm fuv})^2 - 0.7409 \times a({\rm fuv})^3 + 0.2246 \times
a({\rm fuv})^4 \big ) & & \end{aligned}\ ] ] using the lmc extinction law @xmath30 , we then derive : @xmath31 and compute the complete set of @xmath17 using eq.[alambda ] . in order to simplify the calculations in the following we assume @xmath22=0 .
we adopted the lmc extinction law as a compromise between a milky way extinction law , with a strong bump at 2175  , and the small magellanic cloud , with almost no bump at @xmath322000  .
recent analysis have also pointed out that a lmc extinction law can reproduce galex observations better than a mw or smc extinction law @xcite    the reddened sed is then subtracted from the original _ dust - free _ sed , providing ( once integrated over all wavelengths ) the amount of energy absorbed by dust which , assuming an energetic balance , corresponds to the total infrared radiation ( tir ) emitted by the galaxy . finally , by convolving the reddened sed with the galex - fuv filter we estimated the fuv flux and the @xmath6 ratio for each model .
the relations between @xmath7 and @xmath6 so obtained for different values of @xmath11 and stellar metallicity are shown in fig.[tauafuv ] .
+ the mean age of the stellar populations plays a crucial role in the @xmath6 vs @xmath7 relation , as already pointed out by several theoretical ( e.g. @xcite ) and observational ( e.g. @xcite ) studies : for the same @xmath6 ratio active star forming galaxies ( @xmath337 gyr ) are more attenuated than objects with low ssfr ( @xmath34 4 - 5 gyr ) or quiescent galaxies ( @xmath35 2 gyr ) . on comparison
, we find a very weak dependence ( @xmath360.1 mag ) of the @xmath6 vs. @xmath7 on stellar metallicity ( see fig.[tauafuv ] ) . in table
[ convfuv ] we provide the best polynomial fit ( averaged over the whole metallicity range considered ) to the @xmath7 and @xmath6 for different values of @xmath11 .    the origin of the _ age effect
_ is clearly visible in fig.[enbudg ] where we compare the amount of the energy absorbed by dust at long ( @xmath37 ) and short ( @xmath38 ) wavelengths for different values of @xmath11 .
for @xmath39 gyr the high - energy photons ( mainly uv ) contribute to less than the 50% of the total energy absorbed by dust and the far infrared emission is mainly due to the re - emission of the stellar radiation emitted by intermediate age stars in the optical . only for
@xmath40 gyr the uv radiation dominates the dust heating , contributing more than the 75% of the whole energy absorbed and re - emitted in the far infrared .
this also implies that the far - infrared radiation is not always a good star formation indicator . + we remark that our results do not strongly depend on the parameters adopted in our model ( e.g. dust geometry , attenuation law ) .
the variation in the @xmath6 vs. @xmath7 relations when different attenuation laws are considered is @xmath410.2 mag , as discussed in appendix a. this is confirmed by the good agreement between our model and the one of @xcite ( black dotted line in fig.[tauafuv ] ) in the case of strong star forming systems and starburst galaxies ( @xmath338 gyr ) .
our results are also independent of the shape of the adopted sfh and galaxy age .
in fact the @xmath6 vs. @xmath7 relation mainly depends on the specific star - formation and , for the same specific star formation activity , it is not strongly affected by the real shape of the past sfh ( see appendix a ) .
moreover , in fig.[tauafuv ] we show the @xmath6 vs. @xmath7 relation proposed by @xcite ( empty pentagons ) for a birthrate parameter @xmath420.06 ( the lowest birthrate parameter considered in their model , corresponding to @xmath344.1 gyr in our formalism ) .
the level of agreement ( @xmath320.2 mag ) between @xcite and our model is quite comforting , considering that the @xcite model is based on different sfh ( i.e. exponential+bursty model ) and geometry ( i.e. the time - dependent scenario proposed by @xcite ) . +
as already discussed by @xcite , the flux ratio method would not be applicable in the case of an embedded starburst in a galaxy with a second older , less embedded stellar population . at uv and ir wavelengths , the starburst would dominate , but at optical and near - ir wavelengths the older population would dominate
. however , at least for spiral galaxies in the local universe this seems not to be the case @xcite and smooth sfh like an a la sandage and an exponential well reproduce the observations ( e.g. @xcite ) .
we therefore conclude that our results are model - independent within 0.2 mag ( see appendix a ) .
in order to quantify the impact of the _ age effect _ on the estimate of the uv attenuation , we computed @xmath7 for a sample of spiral galaxies and compared our estimate with the one obtained using the standard relation calibrated on active star - forming galaxies ( e.g. @xmath338 gyr , @xmath434 mag , @xmath440.5 ) .
the sample here adopted is an extension of the optically selected sample described in @xcite , composed by late - type galaxies ( later than s0a ) including giant and dwarf systems extracted from the virgo cluster catalogue ( vcc , @xcite ) and from the cgcg catalogue @xcite .
the data include galaxies in the virgo , abell1367 and abell262 clusters and part of the coma - a1367 supercluster ( including cluster and field galaxies ) observed as part of the all - sky imaging survey ( ais ) and of the nearby galaxy survey ( ngs ) carried out by galex in two uv bands : fuv ( @xmath45 ) and nuv ( @xmath46 ) .
we include in our analysis only late - type galaxies detected in both nuv and fuv galex bands and in both 60 @xmath47 and 100 @xmath47 iras bands : 191 galaxies in total ( @xmath3270% in high density environments ) .
uv and far - infrared data have been combined with multifrequency data available . these are optical and near - ir h imaging @xcite , most of which are available from the goldmine galaxy database @xcite ( http://goldmine.mib.infn.it ) . data from uv to near - ir
have been corrected for galactic extinction according to @xcite . + we assume a distance of 17 mpc for the members of virgo cluster a , 22 mpc for virgo cluster b , and 32 mpc for objects in the m and w clouds @xcite .
members of the cancer , a1367 , and coma clusters are assumed to lie at distances of 65.2 , 91.3 , and 96 mpc , respectively .
isolated galaxies in the coma supercluster are assumed at their redshift distance , adopting @xmath48 = 75 @xmath49 .
the fuv attenuation of each galaxy in our sample has been computed using the following sed fitting procedure , implemented into the galaxy observed simulated sed interactive program ( gossip , @xcite ) . in order to take into account the _ age effect
_ we reddened each synthetic sed and fit them to the observed fuv - to - near - infrared sed obtained using all the available photometric magnitudes .
+ in details , for each galaxy we compute the observed @xmath6 ratio using iras and galex observations .
the tir flux emitted in the range 1 - 1000 @xmath47 , is obtained following @xcite : @xmath50 where @xmath51 is the far - infrared flux , defined as the flux between 42 and 122 @xmath47 @xcite : @xmath52\ ] ] and @xmath53 and @xmath54 are the iras flux densities measured at 60 and 100 @xmath47 ( in jansky ) . using the relations presented in the previous section , we then convert the observed @xmath6 ratio into @xmath7 for each value of @xmath11 and @xmath55 considered in our model and determine @xmath17 as described in sec.[model ] . each @xcite sed is then reddened with the @xmath17 obtained following this procedure and fitted to the observed ( i.e. not corrected for internal extinction ) sed , using a @xmath56 fitting technique similar to the one described in @xcite and assuming a conservative photometric error of 0.15 mag for each band .    for each galaxy , the value of @xmath56 determines the weight of a given model ( as @xmath57 ) and a probability distribution function ( pdf ) for @xmath11 can be build by combining weights for each model .
we then normalize the final pdfs and use the peak value ( i.e. the most probable one ) as the best estimate of @xmath11 and the range of @xmath11 containing 68.2% of the pdf s area ( equivalent to the use of constant @xmath56 contours ) as estimate of the 1 @xmath58 error .
the sed fitting provides not only the best value of @xmath11 but also the right estimate of @xmath7 and its uncertainty .
the combined use of the observed @xmath6 ratio for each galaxy and of the relations presented in the previous section makes this method ideal to estimate the uv attenuation for a wide range of morphological types and star formation histories .
the difference between the fuv attenuation obtained with the sed fitting technique ( @xmath59 ) and the one obtained by blindly applying a constant conversion calibrated on star forming galaxies and starbursts ( @xmath338 gyr , @xmath60 consistent with the one presented by @xcite ) is show in fig.[diffafuv_tau ] as a function of @xmath11 ( for clarity only the errors on a(fuv ) are shown ) . as already shown in fig.[tauafuv ] , for low values of @xmath11 the standard conversion overestimate the uv dust attenuation : in particular for @xmath32 30% ( 59/191 galaxies ) of our sample the use of a constant @xmath6 vs. @xmath7 relation leads to an overestimate of more than 0.5 mag of the uv attenuation and for @xmath3216% this systematic error exceeds 1 mag . in order to investigate the properties of the most deviating objects in fig.[diffafuv_tau ] we divide our sample according to their hi - deficiency and optical linear diameter @xmath61 : hi - def = @xmath62 @xcite ] : a value of @xmath63 = 0.4 has been used to separate healthy ( blue circles ) from hi - deficient star forming spirals ( red triangles ) .
hi - deficient spirals are among the most deviating objects in fig.[diffafuv_tau ] , consistent with the fact that gas deficient objects have normally a ssfr significantly lower than that expected from their luminosity @xcite .
however also some healthy spirals show a significant offset from the _ age independent _
@xmath6 vs. a(fuv ) relation .
these are massive early type spirals with low ssfr , similar to m31 in the local group ( i.e. @xmath647 mag , @xcite ) . + the @xmath6 ratio can not thus be considered as a good proxy of the uv attenuation for sample of galaxies spanning a wide range of ssfr , in particular in clusters , and relations like the @xmath2 relation @xcite can not be blindly used to determine @xmath7 .
this is shown in fig .
[ afuv_beta ] where we compare the @xmath7 vs. @xmath3 relation obtained from our sed fitting technique ( filled symbols ) with the one obtained when @xmath7 is computed using the @xmath7 vs. @xmath6 relation for strong star forming galaxies ( empty symbols ) . the @xmath3 parameter has been computed from the fuv - nuv colour following @xcite : @xmath3=2.201@xmath65 - 1.804 .
galaxies with @xmath66 ( @xmath670.75 ) are mainly hi - deficient galaxies and have a fuv dust attenuation @xmath322 mag lower than the one obtained using the standard @xmath68 relation : a difference significantly larger than the typical error on the estimate of a(fuv ) ( @xmath410.5 mag ) by our fitting technique .
the systematic errors introduced in the data is similar when the @xmath6 ratio is not available and @xmath7 is determined using empirical methods based on the @xmath69 colour like the one proposed by @xcite , and once again calibrated on strong star forming galaxies ( solid and dashed line in fig.[afuv_beta ] ) : for @xmath70 - 0.2 , a uv attenuation in the range 2.99@xmath163.32 mag is predicted , whereas our method gives an average value of @xmath71 1 mag .
this is mainly due to the fact that in low star forming systems the uv spectral slope is strongly contaminated by the old stellar populations , whose contribution increases the value of @xmath3 even if the fuv attenuation is low .
this result is consistent with the recent analysis of a sample of @xmath321000 galaxies selected from sdss presented by @xcite , who found that a large fraction of galaxies having red @xmath69 colours is also characterized by large 4000   break ( @xmath72 ) , suggesting that part of the dust heating comes from old stellar populations and not from extremely obscured star forming regions .
we can therefore conclude that the @xmath2 relation can be blindly used to estimate @xmath7 only for @xmath73 ( @xmath740.75 ) .
it clearly appears that a correct estimate of the uv attenuation requires information about the shape of the galaxy sed . a sed fitting technique like the one here described represents therefore the best method available to quantify @xmath7 and properly correct uv observations .
unfortunately sed fitting is only possible when large multiwavelength data sets are available , which is not always the case . for this reason we investigated the possibility to derive different techniques which can be used to estimate @xmath7 when sed fitting is not possible .
ideally , these techniques will have a small rms and not display the systematic shift relative to the galaxy age observed in fig.[diffafuv_tau ] .
the crucial step is to find a good proxy for @xmath11 ( and therefore the shape of the sed ) .
this is particularly important for @xmath75 gyr where the @xmath6 vs. @xmath7 relation strongly depends on the age of the stellar populations .
the first natural choice is to look for a colour with the largest possible dynamical range , sensitive to small changes in the shape of the sed @xcite .
therefore in fig.[fuvh_tau ] we plot the distribution of @xmath11 for our sample as a function of the observed ( i.e. not corrected for internal extinction ) @xmath76 colour . for @xmath137
gyr the two variables are not correlated and galaxies show approximately the same @xmath76 colour independently from the value of @xmath11 , reflecting the large error on the estimate of @xmath11 .
this is as expected since in young stellar populations ( @xmath137 gyr ) the variations in the observed @xmath76 are mainly due to dust attenuation and not to age .
however for @xmath777 gyr ( i.e. the range in which we are interested ) the @xmath76 is a very good proxy for @xmath11 ( pearson correlation coefficient @xmath78 ) . a simple
lest - square fitting gives us : @xmath79    which can be used to estimate @xmath11 from the @xmath76 colour .
the dispersion on this relation varies quite remarkably with the colour .
we therefore determine the typical dispersion within 0.5 mag wide fuv - h bins by combining the @xmath11 s pdf for each galaxy in the bin and estimating the 1@xmath58 error as described in the previous section .
the dispersion in the relation ( indicated by the shaded area in fig.[fuvh_tau ] ) is @xmath80 for 6@xmath819 mag , increasing significantly ( up to @xmath821 ) for redder or bluer colours .
similar relations are found when we consider different far - ultraviolet - optical colours and are presented in table [ taucol ] .
we therefore propose two different ways to determine @xmath7 depending on whether or not far infrared observations are available . in appendix
b we also provide similar recipes in order to determine @xmath83 in the case fuv observations are not available .
+ i ) * far - infrared observations are available . *
determine the observed @xmath6 ratio and use one of the relations presented in table [ taucol ] to determine @xmath11 .
finally use the value of @xmath11 so obtained to choose the conversion between the @xmath6 ratio and @xmath7 from the relations provided in table [ convfuv ] . as discussed
above we suggest using the colour covering the widest possible dynamical range , i.e. first choice @xmath76 , last choice @xmath84 .
the error on the estimate of the fuv dust attenuation ( @xmath85 ) obtained with this method depends on the observational errors on the fuv - optical / near - infrared colour and on the intrinsic dispersion of the colour-@xmath11 relation adopted . in order to estimate @xmath85 ,
for each galaxy in our sample we generated 1000 random galaxies having @xmath76 colour following a gaussian distribution centered on the observed colour with @xmath58=0.2 mag ( more accurate estimate of the @xmath76 colour would correspond to a lower uncertainty on @xmath7 ) .
for each @xmath76 colour the @xmath11 s pdf are used to random generate the correspondent value of @xmath11 vs @xmath11 relation in the estimate of @xmath11 . ] and then the final value of @xmath7 by applying the relations in table [ taucol ] . in fig.[test ] ( upper panel ) we compare this method with the sed fitting technique .
the shaded area indicates the 1 @xmath58 error on the estimate of @xmath7 .
as expected , the error on the estimate of @xmath7 varies significantly with @xmath11 : from @xmath32 + 0.1/-0.25 mag for @xmath336 gyr ( @xmath865 mag ) to a maximum of + 0.4/-0.8 mag for @xmath87 gyr ( 7@xmath818.5 mag ) .
even if the uncertainty for low @xmath11 is quite large ( due to the large variation of the @xmath6 vs. @xmath7 relation with @xmath11 ) , it is clear that our empirical method is able to remove the systematic overestimate of @xmath7 when the age independent @xmath6 vs. @xmath7 conversion is used ( filled symbols in fig.[test ] ) .
+ ii ) * far - infrared observations are not available . * in this case the first step is to find a way to estimate the @xmath6 ratio from the available observations .
recently various methods have been proposed based on the use of the @xmath69 colour ( or @xmath3 , @xcite ) , h - band luminosity @xcite , gas metallicity @xcite , effective surface brightness @xcite , @xmath72 or optical and ultraviolet colours @xcite .
once the @xmath6 ratio has been determined it is possible to proceed as described above : i.e. estimate @xmath11 from an ultraviolet - optical / near infrared colour ( see table [ taucol ] ) and than convert @xmath6 in @xmath7 using the proper relation in table [ convfuv ] .    to test this second method we estimate the @xmath6 ratio from the @xmath69 colour @xcite : @xmath88 and then determine @xmath11 from the @xmath76 colour .
the typical uncertainty in this method has been computed following the same procedure described in the previous point .
also in this case our recipe is able to remove the systematic overestimate on @xmath7 .
this is clearly evident in fig.[test ] ( lower panel ) where our method is compared with the age independent method based on the @xmath69 colour proposed by @xcite . unfortunately the uncertainty on the determination of @xmath7 considerably increases to @xmath32 + 0.5/-0.7 mag for @xmath336 gyr ( @xmath865 mag ) reaching a maximum of @xmath32 + 1/-1.2 mag for @xmath87 gyr ( 7@xmath818.5 mag ) when the total infrared luminosity is estimated from empirical relations involving colours or luminosities . as discussed by @xcite
this is mainly due to the intrinsic scatter in the empirical relation used to determine the @xmath6 ratio .
+ in the very unlikely scenario in which there is no sufficient data to estimate the @xmath6 ratio and/or @xmath11 as described above , the only possibility is to use the morphological type to obtain a rough estimate of @xmath7 . in table
[ afuv_type ] we give the average value and standard deviation of @xmath7 obtained for our sample in bins of morphological type . for our sample this method provides an estimate of @xmath7 with an average error of @xmath320.7 mag .
however we recommend to apply this technique only when accurate an morphological classification is available ( i.e. the local universe ) , otherwise the error on @xmath7 will be considerably larger than the one obtained for our sample .
the simple recipes presented here have the widest scientific application possible , requiring only an uv and uv - optical colour . in order to show their real power , in the following we will discuss two different applications of these methods to the study of the uv properties of nearby galaxies and compare them with the age independent techniques usually adopted .      the _ age effect _ on the @xmath6 vs. @xmath7 relation
is particularly important in panchromatic studies of resolved galaxies @xcite since in different regions the dust can be heated up by different stellar populations @xcite .
+ we test our method on the nearby spiral galaxy ngc 4569 ( m90 ) , the brightest spiral galaxy in the virgo cluster recently studied by @xcite .
ngc 4569 can be considered as the prototype of hi deficient galaxy having only about one - tenth of the atomic gas of a comparable field galaxy of similar type and size .
this galaxy has a truncated h@xmath89 and hi radial profile ( at a radius of @xmath32 5 kpc , @xcite ) and shows significant colour gradients with star formation activity only in the central 5 kpc .
therefore , given its significant age gradients , the contribution of uv photons to the dust heating probably varies with the distance from the galaxy center , making this object ideal to compare with our recipes with age independent techniques .
as described in @xcite we determined the @xmath6 ratio profile combining galex - fuv image and spitzer 24 , 70 , and 160 @xmath90 radial profiles .
then , we use the observed @xmath76 colour profile to estimate @xmath11 at each radius , determining which @xmath6 vs. @xmath7 relation to use at each radius . the @xmath7 profile so obtained is shown in fig.[n4569 ] ( left panel , empty circles ) and compared to the one estimated using the age independent @xmath6 vs. @xmath7 conversion @xcite .
as expected , our method gives a fuv attenuation @xmath321.5 mag lower than the standard conversion implying a factor @xmath324 difference in the sfr surface density profile obtained from the fuv profile ( fig.[n4569 ] , right panel ) .
the results obtained with our technique are supported by the very good agreement with the sfr density profile independently obtained from the h@xmath89 line ( corrected for extinction using balmer decrement ; fig .
[ n4569 ] right panel , empty triangles ) . in fact only a recent ( @xmath91 yr ) starburst @xmath32100 times stronger than the normal star formation rate in ngc  4569 could reconcile the difference between the sfr obtained from the h@xmath89 and from the fuv corrected with the standard methods @xcite : an extremely unlikely scenario as discussed by @xcite and @xcite .
+ this exercise shows the crucial importance of the age dependent @xmath6 vs @xmath7 relations in the panchromatic study of resolved galaxies .
the need of empirical methods to correct for dust attenuation is particular important for the study of uv properties of large samples of galaxies lacking far - infrared rest - frame data .
of particular importance for our understanding of galaxy evolution is the correct interpretation of the uv - optical colour magnitude ( cm ) relation , as recently discussed by @xcite , @xcite and @xcite .
it is very difficult to use a single correction technique able to deal simultaneously with starbursts , low ssfr objects and elliptical galaxies .
of particular importance is the correct estimate of @xmath7 for transition objects : galaxies with low ssfr for which the standard corrections calibrated on active star forming systems are likely to be not valid . in order to quantify the influence of the _ age effect _ on our interpretation of the uv - optical colour magnitude relation
we compare our method with the age independent recipes using a sample of galaxies extracted from galex observations of the coma cluster region ( gi - cycle 1 , cortese et al .
, in preparation ) .
this sample includes all galaxies detected in both fuv and nuv bands as well as with sdss - dr6 @xmath92 photometry and with spectroscopic redshift data available ( 833 galaxies ) . since for this large sample
far infrared observations are not available , we use the second method described in sec.[recipe ] , determining the @xmath6 ratio from the @xmath69 colour and @xmath11 from the @xmath93 colour .
the @xmath94 colour-@xmath95 magnitude relation obtained is compared with the observed one ( i.e. not corrected for internal extinction ) in fig.[cm ] ( left panel ) .
whereas galaxies in the blue sequence ( @xmath96 4 mag ) have an average fuv dust attenuation @xmath711.5 - 2 mag , galaxies in the red sequence ( @xmath976 mag ) show ( as expected ) a very low amount of attenuation ( @xmath71 0.5 mag ) and occupy almost the same parameter space as uncorrected data .
the separation between the red and blue sequences is therefore more evident after the correction for dust attenuation .
we remark that our technique should not be blindly applied to elliptical galaxies , since in these objects uv emission comes from old stellar populations and not young stars ( e.g. @xcite ) .
however if an accurate morphological classification is not available the use of the recipes presented here does not introduce a strong systematic bias in the data , as shown in fig.[cm ] ( left panel ) .
this result indicates that our procedures are reliable for old as well as young stellar populations , contrary to previous empirical methods , calibrate and valid only for active star - forming galaxies ( e.g. @xcite ) .
this is clearly evident in fig.[cm ] ( right panel ) where the cm relation obtained with our recipe is compared with the one determined by using the age - independent method to convert the @xmath69 colour into @xmath7 ( e.g. @xcite ) .
the difference between the two cm relations is quite impressive : while no significative difference is observed in blue - sequence galaxies , the red sequence shifts @xmath322 mag towards bluer colours and the gap between the blue and red sequence ( i.e. the so - called `` green valley '' , @xcite ) is considerably reduced .
a similar result has been recently shown by @xcite when comparing the cm relation determined using the balmer decrement as indicator of dust attenuation with the one estimated using the recipes proposed by @xcite .
this result shows how strong the bias can be when blindly using the recipes to estimate @xmath7 without taking into account the _ age effect _ , leading to an incorrect interpretation of the data and reconstruction of galaxy s evolution history .
we have investigated the dependence of the @xmath6 vs. @xmath7 relation on the average age of galaxy stellar populations . our simple method has shown that even for spiral galaxies the use of a standard ( i.e. not age dependent ) conversion of the @xmath6 into @xmath7 leads to a systematic overestimate ( i.e. @xmath981 mag ) of the dust attenuation in galaxies with low specific star - formation , mainly anemic cluster spirals .
this systematic bias strongly affects our interpretation of uv observations and can produce a significant overestimate ( up to a factor 10 ) of the star formation rate , particularly in hi deficient galaxies .
therefore we have developed different methods for determining the uv dust attenuation taking into account the age dependence of the @xmath6 vs. @xmath7 relation .
these recipes require only an uv colour and an uv - optical / near infrared colour providing an estimate of the uv attenuation with an average uncertainty varying from @xmath32 + 0.1/-0.25 mag for @xmath336 gyr ( @xmath865 mag ) to a maximum of + 0.4/-0.8 mag for @xmath87 gyr ( 7@xmath818.5 mag ) , when far - infrared observations are available , and from @xmath32 + 0.5/-0.7 mag for @xmath336 gyr ( @xmath865 mag ) to a maximum of @xmath32 + 1/-1.2 mag for @xmath87 gyr ( 7@xmath818.5 mag ) when far - infrared data are not available .
the small amount of multiwavelength data necessary for their application makes these methods extremely useful for the widest possible range of application eventually providing us with a less biased view of the uv properties of galaxies in the local universe and at higher redshift .
we want to thank the anonymous referee , whose comments and suggestions were extremely useful for improving the present manuscript .
we thank jonathan davies for useful discussions and comments on this manuscript .
lc is supported by the uk particle physics and astronomy research council .
galex ( galaxy evolution explorer ) is a nasa small explorer , launched in april 2003 .
we gratefully acknowledge nasa s support for construction , operation , and science analysis for the galex mission , developed in cooperation with the centre national detudes spatiales of france and the korean ministry of science and technology . this research has made use of the nasa / ipac extragalactic database , which is operated by the jet propulsion laboratory , california institute of technology , under contract to nasa and of the goldmine database
are the results obtained in this work and the recipes proposed to estimate the uv dust attenuation model - dependent ? in order to answer to this question ,
in this section we investigate the dependence of our results on the three free parameters entering our model : the shape of the extinction law , the dust geometry and the star formation history .    in fig.[modgeo ] ( left ) we show the @xmath6 vs @xmath7 relation obtained for different values @xmath11 , a sandwich geometry and three different extinction laws : lmc ( solid line ) , milky way ( dotted line ) and small magellanic cloud ( dashed line ) . as already pointed out by several authors ( e.g. @xcite )
, the @xmath6 vs @xmath7 relation is almost independent of the extinction law with a typical variation limited to @xmath410.2 mag , considerably lower than the effect of @xmath11 on the @xmath6 vs @xmath7 relation .
similar results are obtained if we investigate the effect of different dust geometries . in fig.[modgeo ] ( right ) we compare the results obtained for the sandwich model ( solid line ) with a simple slab geometry ( @xcite , dashed line ) .
we also added the results for a @xcite attenuation law ( dotted line ) . even in this case , we observed no dependence of the @xmath6 vs. @xmath7 relation on the dust geometry within @xmath410.2 mag .
similar results are obtained if we use the homogeneous and clumpy dust models proposed by @xcite ( not shown ) .    finally , we fitted the 191 galaxies in our sample with various sed libraries in order to test the dependence of the results presented in sec.[data ] and sec.[recipe ] from sfh and age .
the sed library used in sec.[model ] has been obtained by fixing the shape of the sfh ( assumed to be a la sandage ) , the age of the galaxy ( assumed equal to 13 gyr ) and the stellar initial mass function ( imf , assumed to be a salpeter imf ) .
we produced different set of sed libraries varying every time one of these three free parameters and we fitted them to the data .
in particular we considered an exponential sfh , a fixed ( @xmath10=13 gyr ) or free ( in the range 0@xmath9915 gyr ) age and a @xcite imf . in fig.[modsfh ] is presented the difference in the estimate of @xmath7 between the library adopted in sec.[model ] and some of the various combination adopted . for each combination
the standard deviation in the estimate of @xmath7 is lower than 0.1 - 0.2 mag .
this result was somehow expected since the @xmath6 vs. @xmath7 relation is only affected by the actual shape of the sed and not by the shape of the past ( older than 1 - 2 gyr ) sfh .
we remark that our results are not applicable only in the case of an embedded starburst in a galaxy with a second older , less embedded stellar population @xcite . at uv and ir wavelengths ,
the starburst would dominate , but at optical and near - ir wavelengths the older population would dominate .
however , at least for spiral galaxies in the local universe this seems not to be the case @xcite and smooth sfh like a la sandage and exponential well reproduce the observations ( e.g. @xcite )    we can therefore conclude that the results obtained from our model and the recipes presented in this work are model - independent at least within @xmath320.2 mag .
this uncertainty is still lower or equal to the uncertainty in the estimate of @xmath7 from our recipes .
in this appendix we discuss the possibility to extend our recipes to the nuv band in the case that galex - fuv observations are available .
the galex - nuv filter ( @xmath46 ) lies in a region where the sed of galaxies with @xmath100 4 gyr is strongly dependent on stellar metallicity for @xmath101 ( e.g. see fig .
6 in @xcite ) .
this strongly affects the @xmath102 vs @xmath83 relation which , for low values of @xmath11 , can vary of even 0.4 mag from @xmath103 to @xmath104 , apparently complicating the use of our model directly on the nuv band .
however for none of the galaxies with @xmath100 4 gyr in our sample the best fitting model has metallicity lower than @xmath103 , consistent with the fact that in the local universe evolved galaxies tend to be metal rich ( e.g. @xcite ) .
therefore for @xmath100 4 gyr we compute the average @xmath102 vs @xmath83 relation by combining only the relations obtained in the metallicity range @xmath105 .
this makes our recipes not valid for @xmath100 4 gyr and @xmath106 .
the average @xmath102 vs. @xmath83 relations for different values of @xmath11 so obtained are presented in table [ convnuv ] . in table
[ taunuv ] are given the relations to use to estimate @xmath11 from the nuv - optical / near - infrared colours .
their typical dispersion is consistent with the one observed in the fuv - optical / near - infrared colours vs. @xmath11 relations . as discussed in sec.[recipe ] the suggested recipes to determine @xmath83 result : + i ) * far - infrared observations are available . *
determine the observed @xmath102 ratio and use one of the relations presented in table [ taunuv ] to determine @xmath11 .
finally use the value of @xmath11 obtained with this procedure to choose the conversion between the @xmath6 ratio and @xmath7 from the relations provided in table [ convnuv ] .
the typical error on the estimate of @xmath83 ( computed through a montecarlo simulation similar to the one described in sec .
4 ) varies between @xmath32 + 0.1/-0.2 mag for @xmath338 gyr , to + 0.3/-0.7 mag for @xmath87 gyr .
+ ii ) * far - infrared observations are not available . * in this case the first step is to find a way to determine the @xmath102 from the observations available . once the @xmath102 ratio has been determined it is sufficient to proceed as described in the previous point : i.e. estimate @xmath11 from an ultraviolet - optical / near infrared colour and than convert @xmath102 in @xmath83 using the proper relation in table [ convnuv ] . the typical error on the estimate of @xmath83 ( computed through a montecarlo simulation similar to the one described in sec .
4 ) varies between @xmath32 + 0.4/-0.6 mag for @xmath338 gyr , and to @xmath32 + 0.9/-1 mag for @xmath87 gyr .
+ in table [ anuv_type ] are presented the average value of @xmath83 obtained for the different morphological types in our sample .
the morphological type can be used to estimate @xmath83 only if the to previous methods can not be applied and an accurate morphological classification is available .
the relationship between the @xmath6 ratio and the fuv attenuation @xmath7 obtained from our model for different values of @xmath11 .
solid and dashed lines show the relations for stellar metallicity @xmath55=2.5 z@xmath15 and @xmath55=0.02 z@xmath15 respectively .
the dotted black line indicates the age independent relation given by buat et al .
empty pentagons show the relation proposed by kong et al .
( 2004 ) for birthrate parameter @xmath420.06.,width=321 ]    the importance of the uv emission in the dust heating .
the ratio of the energy absorbed by dust at @xmath107 and at @xmath108 is shown as a function of @xmath11 ( i.e. the age of the universe at which the galaxy sfr peaks).,width=321 ]     the difference between the fuv attenuation for active star forming galaxies ( @xmath33 8 gyr ; consistent with buat et al .
( 2005 ) ) and the one obtained from our sed fitting technique as a function of @xmath11 .
circles indicate healthy spirals ( hi - def@xmath360.4 ) and triangles hi - deficient star forming objects ( hi - def@xmath980.4 ) .
the distributions of the difference for healthy ( solid histogram ) and deficient ( dashed histogram ) spirals are shown in the upper - right panel .
the unshaded area indicates the range of @xmath11 on which the standard @xmath6 vs. @xmath7 relation for star - forming galaxies is usually calibrated.,width=321 ]     the @xmath7 vs. @xmath3 relation for our optically selected sample .
a(fuv ) is computed from our sed fitting technique ( filled symbols ) and using the standard recipe for active star forming galaxies ( empty symbols ) respectively .
the solid and dashed lines indicate the relations between @xmath7 and @xmath69 proposed by salim et al .
symbols are as in fig.[diffafuv_tau].,width=302 ]     relation between the observed ( i.e. not corrected for internal extinction ) @xmath76 colour and @xmath11 for our sample .
the dotted line indicates the best linear fit for galaxies having @xmath777 gyr .
the shaded area indicate the typical dispersion of the data computed in fuv - h bins 0.5 mag wide .
symbols are as in fig.[diffafuv_tau].,width=302 ]     residuals in the estimate of @xmath7 obtained with the different methods described in sec .
the @xmath6 ratio is estimated using iras observations ( upper panel ) and the @xmath69 colour ( lower panel ) respectively .
the shaded areas show the typical uncertainties in the new recipes presented in this paper , estimated by generating 1000 random galaxies having @xmath76 colour following a gaussian distribution centered on the observed colour with @xmath58=0.2 mag ( see sec . 4 for details ) .
for comparison , the residuals obtained when an age independent method is used are indicated by the filled symbols .
symbols are as in fig.[diffafuv_tau].,width=302 ]           difference in the estimate of @xmath7 with the sed technique described in sec.[data ] ( @xmath10913 gyr , a la sandage sfh and salpeter imf ) and various combination of sfh ( sandage or exponential ) , age ( variable or fixed to 13 gyr ) and initial mass function ( imf , salpeter or chabrier ) .
the dashed lines indicates the 1@xmath58 difference in the estimate of @xmath7 .
, width=320 ]