numerous observations , from galaxy rotation curves to gravitational lensing , to cosmic microwave background radiation , all point to the existence of dark matter which is not made of ordinary atoms , but , rather , of some new , yet undiscovered particles .
in particular , the abundance of dark matter has been precisely determined recently by the wmap  @xcite .
however , the nature of dark - matter particles remains unknown .
one attractive candidate is a sterile neutrino with mass of several kev and a small mixing with the ordinary neutrinos .
if such a particle exists , it could be produced in the early universe with the right abundance to be dark matter either from neutrino oscillations  @xcite , or from some other processes , for example from inflaton decay  @xcite .
the same particle can explain the observed velocities of pulsars because its emission from a cooling neutron star would be anisotropic , hence providing the neutron star with a recoil velocity of a sufficient magnitude @xcite .
in addition , the neutrino kick can enhance the convection during the first second of the supernova explosion , which can increase the energy of the supernova shock and bring the supernova calculations in better agreement with the observations @xcite .
sterile neutrinos can also help the formation of supermassive black holes in the early universe @xcite .
although the dark - matter sterile neutrinos have interactions that are too weak to be discovered in the laboratory , some other sterile neutrinos may exist .
most models of neutrino masses introduce sterile ( or right - handed ) neutrinos to generate masses of the ordinary neutrinos via the seesaw mechanism  @xcite .
the name _ sterile _ was coined by  @xcite .
many seesaw models assume that sterile neutrinos have very large masses , which makes them unobservable .
however , a number of theoretical models predict the existence of sterile neutrinos at or below the electroweak scale  @xcite . there is an indication that some light sterile mass eigenstates may cause the neutrino flavor transformation seen in the los alamos liquid scintillator neutrino detector ( lsnd ) experiment @xcite .
the results of lsnd are being tested by the miniboone experiment .
one can expect that there are at least three light sterile neutrinos @xcite , in which case the neutrino oscillations can explain the baryon asymmetry of the universe @xcite .
although dark - matter sterile neutrinos are stable on cosmological time scales , they nevertheless decay . the mass matrix is assumed to have the seesaw form with the majorana mass term of the order of several kev and with a much smaller dirac mass term .
the diagonalization of such a mass matrix yields the light mostly active majorana mass eigenstates and the heavy sterile majorana eigenstates , which can decay into the light ones .
the dominant decay mode , into three light neutrinos , is `` invisible '' because the daughter neutrinos are beyond the detection capabilities of today s experiments .
the most prominent `` visible '' mode is decay into one active neutrino and one photon , @xmath0 . assuming two - neutrino mixing for simplicity
, one can express the inverse width of such a decay as @xmath1 where @xmath2 is the mass and @xmath3 is the mixing angle .
for mass units we set @xmath4 .    since this is a two - body decay , the photon energy is half the mass of the sterile neutrino .
the monochromatic line from dark matter decays can , in principle , be observed by x - ray telescopes .
no such observation has been reported , and some important limits have been derived on the allowed masses and mixing angles  @xcite .
these constraints are based on different astrophysical objects , from virgo and coma clusters , to large magellanic clouds , to the milky way halo and its components , etc .
there are different uncertainties in modeling the dark matter populations in these objects .
different groups have also used very different methodologies in deriving these bounds : from a conservative assumption that the dark - matter line should not exceed the signal , to more ambitious approaches that involved modeling the signal or merely fitting it with a smooth curve and requiring that the line - shaped addition not affect the quality of the fit .
the x - ray limits are usually shown based on the assumption that the sterile neutrinos constitute the entire dark matter ( @xmath5 , where @xmath6 and @xmath7 refer to the sterile neutrinos and the dark matter , respectively ) .
however , the x - ray observations exclude certain masses and mixing angles even if the sterile neutrinos make up only a fraction of dark matter .
the production via mixing @xcite provides the lowest possible abundance of sterile neutrinos . in fig .
[ fig : models ] we show both the bounds based on the assumption @xmath5 and the model - independent exclusion region  @xcite based on the production only via the @xcite mechanism .
@xcite used the analytical fit to the numerical calculation of @xcite , and re - scaled the x - ray flux in proportion to the amount of relic neutrinos produced for a given mass and mixing angle .
a different set of bounds comes from the observations of the lyman - alpha forest  @xcite , which limit the sterile neutrino mass from below .
based on the high - redshift data from sdss and some modeling of gas dynamics , one can set a limit as strong as 14  kev  @xcite .
however , the high - redshift data have systematic errors that are not well understood  @xcite , and more conservative approaches , based on the relatively low - redshift data , have led to some less stringent bounds  @xcite .
more recently , @xcite have reanalyzed the high - redshift data and arrived at the bound @xmath8  kev .
the mass bound as quoted depends on the production mechanism in the early universe .
the lyman - alpha observations constrain the free - streaming lengths of dark matter particles , not their masses . for each cosmological production mechanism ,
the relation between the free - streaming length and the mass is different  @xcite .
for example , the bound @xmath9  kev  @xcite applies to the production model due to  @xcite . even within a given cosmological scenario ,
there are uncertainties in the production rates of neutrinos for any given mass and mixing angle  @xcite .
these uncertainties may further affect the interpretation of the lyman - alpha bounds in terms of the sterile neutrino mass .
the x - ray bounds depend on a combination of mass and mixing angle as in eq .
( [ meanlife ] ) , regardless of one s cosmological assumptions , as long as all the dark matter is comprised of sterile neutrinos .
however , the x - ray bounds are subject to uncertainties in modeling the dark matter populations of different objects .
if the sterile neutrinos make up only a part of dark matter , none of the above bounds apply .
also , if inflation ended with a low reheating temperature , the bounds are significantly weaker  @xcite .
it should also be mentioned that the lyman - alpha bounds appear to contradict the observations of dwarf spheroidal galaxies  @xcite , which suggest that dark matter is warm and which would favor the 15  kev mass range for sterile neutrinos .
there are several inconsistencies between the predictions of n - body simulations of cold dark matter ( cdm ) and the observations @xcite .
each of these problems may find a separate independent solution .
perhaps , a better understanding of cdm on small scales will resolve these discrepancies .
it is true , however , that warm dark matter in the form of sterile neutrinos is free from all these small - scale problems altogether , while on large scales wdm fits the data as well as cdm .    several independent considerations , from the pulsar kicks to small - scale behavior of dark matter , favor sterile neutrinos as the dark matter candidate .
it is important , therefore , to examine the effects of sterile - neutrino dark matter on star formation and reionization of the universe .
@xcite have attempted to study such effects for some kind of warm dark matter , but they have used a linear matter power spectrum and assumed that the dark matter particles have a thermal distribution , which is not the case for sterile neutrinos @xcite .
even more importantly , they have not considered the effects of sterile neutrino decays and the resulting x - ray background .
it was shown by @xcite that the x - rays produced by the decays of relic sterile neutrinos can speed up the formation of molecular hydrogen , which can help the cooling of the gas and the star formation , which can lead , in turn , to a reionization of the universe at the redshift consistent with the wmap results .
more specifically , the wmap ( three years ) measurement of the reionization redshift @xmath10 @xcite has posed a challenge to theories of star formation . on the one hand
, stars have to form early enough to reionize gas at redshift  11 . on the other hand ,
the spectra of bright distant quasars imply that reionization must be completed by redshift  6 .
stars form in clouds of hydrogen which collapse at different times , depending on their sizes : the small clouds collapse first , while the large ones collapse last .
if the big clouds must collapse by redshift  6 , then the small halos must undergo the collapse much earlier .
it appears that the star formation in these small halos would have occurred at high redshift , when the gas density was very high , and it would have resulted in an excessive thompson optical depth @xcite . to be consistent with wmap , the efficiency for the production of ionizing photons in minihalos must have been at least an order of magnitude lower than expected  @xcite .
one solution is to suppress the star formation rate in small halos by some dynamical feedback mechanism .
the suppression required is by at least an order of magnitude .
alternatively , one can consider warm dark matter , in which case the number of small clouds is strongly suppressed .
however , for some wdm candidates ( for example , the gravitino ) , the lack of power on small scales in the spectrum of density perturbations leads to an unacceptable delay in the onset of star formation @xcite .
this is not the case for sterile neutrinos , which decay slowly during the `` dark ages '' preceding the birth of stars .
the decay photons could ionize gas , and , although this ionization is too weak to affect the wmap measurements directly  @xcite , the increase in the ionization fraction could cause a rapid production of molecular hydrogen at redshifts as high as 100  @xcite .
molecular hydrogen is the essential cooling agent that allows the gas to cool and collapse to form the first stars .
@xcite have not considered the effects of gas heating associated with the decays , which may , in some cases , thwart the effects of ionization .
the increase in the fraction of molecular hydrogen can help cooling the gas , but if the heating from sterile neutrino decays is significant , it can overcome the radiative cooling and stymie the collapse and star formation .
in this paper we will examine thermal evolution of gas clouds in detail , taking into account both effects of the sterile neutrino decays , namely the ionization and the heating of gas .
the paper is organized as follows . in section [ sec_h2 ]
we will discuss the role of sterile neutrino decays and the production of molecular hydrogen in the cooling of the primordial gas clouds .
next , in sections 3 - 6 we will briefly describe the top - hat overdensity simulation .
the results are shown in section  [ sec_results ] .
the existence of the first objects is a direct consequence of the growth of the primordial density fluctuations . at the beginning , there are linear density perturbations which expand with the overall hubble flow .
subsequently , these perturbations can grow and form primordial clouds .
clouds with enough density contrast decouple from this flow and start to collapse .
the kinetic energy of the infalling gas is dissipated through shocks and the cloud becomes pressure supported .
the subsequent evolution of the cloud is determined by its ability to cool sufficiently fast .
clouds that do not cool fast enough stay in a pressure - supported state and do not form stars .
the existence of the efficient cooling mechanism is necessary to continue the collapse of the cloud , its subsequent fragmentation and star formation .    in the absence of metals
, the most important cooling mechanism is cooling by hydrogen atoms and molecules .
this mechanism involves the collisional excitation of h or h@xmath11 , subsequent spontaneous de - excitation , and photon emission .
the emitted photons can escape and take away energy , hence reducing the kinetic energy of the collapsing cloud .
cooling by h dominates for temperatures above @xmath12  k , while the h@xmath11 cooling is more important at lower temperatures .
the temperature of the primordial cloud can be higher than @xmath12 k only during a short period of virialization of the baryonic gas .
therefore , h@xmath11 cooling is the most important cooling mechanism during almost entire evolution of the halos , and it is crucial for the star formation .    in primordial gas clouds , h@xmath11 molecules can form mainly through the coupled reactions @xmath13 or @xmath14 in which @xmath15 and @xmath16 act only as a catalyst .
x - ray radiation can increase the production of the h@xmath11 by enhancement of the ionization fraction .    as discussed above ,
sterile neutrinos with the mass of several kev and a small mixing are stable on cosmological time scales .
nevertheless some of them can decay .
the decay channel important for us is @xmath17 its inverse width is given in equation ( [ meanlife ] ) .
the energy of the photon is equal to @xmath18 where @xmath19 is the mass of the active neutrino .
x - ray photons from the sterile neutrinos decays can increase the ionization fraction and the h@xmath11 production .
this can subsequently lead to a speed - up of the gas cooling and star formation .
we will follow the evolution of the baryonic top - hat overdensity , the gas temperature and its h@xmath11 and @xmath15 fraction .
our goal is to juxtapose the evolution of the gas temperature in the primordial clouds in the @xmath20cdm model and the wdm model with kev sterile neutrinos .
we assume that the dark matter is smoothly distributed in the region we are considering .
we will use the parameters from the best fit to the three years wmap data @xcite , namely @xmath21 km @xmath22 @xmath23 , @xmath24 , @xmath25 , @xmath26 , @xmath27 , @xmath28 , @xmath29 .
let us take some spherical region with the top - hat overdensity @xmath30 .
it is well known that such an overdensity will evolve according to the following equation @xcite @xmath31 where the parameter @xmath32 is related to the redshift @xmath33 and redshift of virialization @xmath34 through @xmath35    further evolution depends on the type of the matter in the overdense region .
if it is the dark matter then after virialization its density remains constant forever .
the situation is different in the case of baryons .
if there exists an efficient cooling mechanism , then the density gradually increases . otherwise , the density remains constant , and there is no star formation in the halo .    in reality
an overdense region consists of two components namely , baryonic gas and dark matter .
baryonic density fluctuations are suppressed before recombination .
in contrast , the dark matter density fluctuations can start to grow well before this period , so its initial overdensity is larger than that of the baryons .
dark matter forms the potential wells into which the baryonic matter falls .
overall , the baryon density follows the dark matter density in the early universe .
therefore , one can consider the primordial fluctuations as the halos consisting of baryonic and dark matter , whose overall overdensity follow the top - hat solution ( eqs .
[ th1 ] and [ th2 ] ) .
of course , this assumption breaks down at the later stages of the halo evolution , e.g. after the virialization of the dark matter .
following @xcite , let us assume that the evolution of the baryonic and dark matter density is described by equations ( [ th1 ] ) and ( [ th2 ] ) until @xmath32 reaches @xmath36 .
this occurs at the redshift @xmath37 after this redshift , we assume that the density of the halo @xmath38 stays constant e.g. it is equal to the value @xmath39 where @xmath40 is the critical density of the universe .
we note that the overdensity @xmath30 in equation ( [ th1 ] ) is given by @xmath41 where @xmath42 is the combined average density of dark and baryonic matter .
we have taken the ratio of the baryonic to the dark matter mass in the halo to be the standard cosmological ratio of @xmath43 .
the density of the baryonic matter can be expressed as @xmath44
the gas temperature evolution is governed by the equation of energy conservation @xmath45 where @xmath46 , @xmath47 and @xmath20 are the pressure , the internal energy per unit mass and the cooling / heating function respectively .
we use the equation of state of perfect gas @xmath48 where @xmath49 is the adiabatic index .
the internal energy per unit mass can be described as @xmath50 where @xmath51 is the boltzmann constant , @xmath52 is the gas temperature , @xmath53 denotes the molecular weight , and @xmath54 is the mass of the hydrogen atom . introducing the number density of non dark matter particles @xmath55 we can write @xmath56 as @xmath57    if we use eqs .
( [ p ] ) , ( [ u ] ) and ( [ p2 ] ) than we can cast eq .
( [ u1 ] ) in the form @xmath58    the time @xmath59 can be translated to the redshift @xmath33 in the following way @xmath60 using eqs .
( [ tt ] ) and ( [ hubble ] ) , we obtain @xmath61 to calculate the gas temperature of the primordial cloud we have to integrate this equation .
one can assume that the gas temperature evolution is governed by eq .
( [ tz ] ) except for the redshift range between @xmath62 and @xmath34 . if we simply integrate eq.([tz ] ) to the redshift @xmath34 , than the gas temperature will be much lower than the virial temperature @xmath63 , at least for clouds with mass greater than some critical value  @xcite @xmath64 where @xmath65 , @xmath66 and @xmath67 are respectively the mean molecular weight during virialization , the halo mass and the overdensity of the halo at virialization given by @xcite , which is very well approximated by @xmath68 for redshifts @xmath69 .
therefore we have to take into account shocks and increase the gas temperature to the virial value .
it is reasonable to assume that the gas temperature evolution is linear between the redshift @xmath62 and @xmath34 and is described by the equation @xmath70 where @xmath71 is the temperature at the redshift @xmath62 .
let us consider the following five species : h , h@xmath72 , @xmath73 , @xmath15 , and @xmath74 with the mass fraction @xmath75 @xcite .
let us denote by @xmath76 ( @xmath77h , h@xmath72 , h@xmath11 , @xmath15 and he ) the number density of baryonic component @xmath78 .
we assume that @xmath79 e.g. there is no he or h@xmath11 ionization .
the number density of hydrogen can be written as @xmath80 if we introduce the number fraction of component `` @xmath78 '' @xmath81 then the particle number density @xmath55 and the baryon density @xmath56 are equal to @xmath82 and @xmath83 where @xmath84 is the hydrogen mass fraction .
using eqs .
( [ p2 ] ) , ( [ nbb ] ) and ( [ rob ] ) we can get @xmath85 at the virialization we have typically @xmath86 and @xmath87 , so one can take the mean molecular weight at this time to be @xmath88 .
we can express @xmath49 as @xmath89 where @xmath90 and @xmath91 .
to calculate the cooling / heating function @xmath20 in eq .
( [ tz ] ) we have to take into account all of the relevant chemical and thermal processes and follow the number density evolution of different baryonic components . all chemical and thermal processes included in our calculation
are summarized in tables [ tab : reaction ] and [ tab : cooling ] .
the evolution of @xmath92 and @xmath93 is described by equations @xmath94 @xmath95 where @xmath96 is the production rate of free electrons by the sterile neutrino decays which is calculated in the appendix . in eq .
( [ xe ] ) only the ionization and recombination of hydrogen are included .
all the rates @xmath97 are given in table [ tab : reaction ] .
we note in passing that neglecting he and h@xmath11 ionization renders the system with fewer free electrons , so the production of h@xmath11 is slower .
in addition , cooling by h@xmath98 is stronger than by h@xmath11 ( see last figure in @xcite ) .
therefore , the cooling we get in this approximation is somewhat less efficient than it should be in reality .
the molecular hydrogen number fraction can be calculated by integration of the following equation @xmath99 where @xmath100 is the molecular hydrogen formation rate via @xmath101 and h@xmath98 channel @xcite . to get eq .
( [ km ] ) we have to make the assumption that after formation of @xmath101 and @xmath102 they are used almost instantaneously in the next reactions ( reactions number 2 and 3 or 5 and 6 respectively , see table [ tab : reaction ] ) .
this is a very good approximation , because in reality the @xmath101 and @xmath102 fractions always remain small .
the cooling / heating function @xmath20 is given in the appendix ( eq . [ coolfunc ] ) .
it includes contributions from the following processes :    1 .
excitation of h by collisions with @xmath15 2 .
recombination to h 3 .
ionization of h by collisions with @xmath15 4 .
photo - ionization of h by the cmb radiation 5 .
excitation of h@xmath11 by collisions with h and h@xmath11 6 .
h@xmath11 formation cooling 7 .
dissociation of h@xmath11 by collisions with h , h@xmath11 and @xmath15 8 .
photo - dissociation of h@xmath11 by cmb photons 9 .
compton cooling 10 .
bremsstrahlung 11 .
photo - ionization of h by the photons from the sterile neutrino decays and h ionization by the secondary electrons produced due to these photons    we have not included the effect of the photo - dissociation of h@xmath11 by the photons from the sterile neutrino decays .
this is because the flux of these photons is small in comparison with the cmb radiation in the 10 - 30 ev range , which is the energy range relevant to the photo - dissociation of h@xmath11 .
the formulae for the contributions to the @xmath20 function from all of these processes are given in table [ tab : cooling ] .
we have assumed the ortho - h@xmath11 to para - h@xmath11 ratio to be 3 : 1 .
we start our calculations at the redshift of recombination @xmath103 , taking as initial conditions @xmath104 , @xmath105 and assuming that the gas temperature is equal to the cmb temperature .
evolution of the halo is insensitive for the initial value of @xmath106 .
we follow the evolution of the density , eqs .
( [ th1 ] ) and ( [ th2 ] ) , and the temperature , eq .
( [ tz ] ) with the cooling / heating function given in eq .
( [ coolfunc ] ) untill the redshift is equal to @xmath62 , _ cf . _
( [ z3pi ] ) . from that time
we assume that the density of the halo stays constant and is equal to its virial value ( [ densvir ] ) .
we assume that the gas starts to be shock heated in this period and reaches its virial temperature at @xmath34 .
the temperature of the gas is assumed to change linearly with the redshift between @xmath62 and @xmath34 , eq .
( [ telin ] ) . for the redshift lower than @xmath34
the gas temperature is calculated once again by integration of eq .
( [ tz ] ) .
following @xcite , we use the following criterion for the successful collapse : @xmath107 where @xmath108 @xcite .
the results of our analysis depend , of course , on the parameters we choose for the dark - matter sterile neutrinos . the x - ray flux , and the resulting increase in the ionization fraction depend on the combination ( @xmath109 ) .
the power spectrum on small scales , which determines the presence or absence of small clouds , depends on the free - streaming length , which is set by the mass ( in every given cosmological scenario ) .    based in part on the existing limits , we have chosen four illustrative benchmark cases to illustrate the effects of ionization and heating due to sterile neutrino decays on the cloud collapse .
we have compared the following four models :    1 .
cold dark matter with no additional radiation besides cmb ( cdm ) .
warm dark matter ( wdm ) with the sterile neutrino mass @xmath110 kev and mixing angle @xmath111 ( wdm1 ) .
3 .   wdm with @xmath112 kev and @xmath111 ( wdm2 ) .
4 .   wdm with @xmath113 kev and @xmath114 ( wdm3 ) .
the values of parameters chosen for these models are shown in fig.[fig : models ] .
some comments are in order .
the choice of wdm1 is obvious because it satisfies all the constraints , including the x - ray bounds and the lyman - alpha bounds .
the masses as large as @xmath110 kev can not be constrained by either chandra or xmm , and the best constraint comes from heao , which is weaker because of its lower energy resolution @xcite .
also , wdm1 corresponds to a very small free - streaming length and , on large scales , it creates structure indistinguishable from cdm .
the only difference between wdm1 and cdm is the additional x - ray flux .
therefore , the comparison of cdm and wdm1 cases provides a clean test of the effects of the increased ionization fraction @xcite .
also , wdm1 provides a realistic set of parameters consistent with all the bounds , as well as the pulsar kick mechanism due to resonant active - sterile neutrino conversions @xcite .
wdm2 is also consistent with all the existing limits , including the most ambitious lyman - alpha bounds  @xcite .
( we remind the reader that @xcite have obtained a 50% weaker limit using the same data . ) the x - ray background is lower in this case than in the case of wdm1 .
finally , wdm3 corresponds to masses that would be ruled out by the `` strong '' lyman - alpha bounds obtained by  @xcite and by  @xcite .
however , one has to remember that the relation between the masses and the free - streaming length depends on the production mechanism of sterile neutrinos .
for example , the free - streaming length in the  @xcite production mechanism is different from the  @xcite mechanism . even within the latter framework ,
additional entropy production can redshift the sterile neutrino momenta and eliminate the lyman - alpha bound  @xcite .
if dark matter is produced from a coupling between the sterile neutrinos and the inflaton  @xcite , the free - streaming length is , again , different .
wdm3 requires redshifting of the neutrino momenta by a factor @xmath115 , as compared to the thermal distribution .
this can be achieved , for example , by entropy production of the order of @xmath116 , which is comparable to the effect of change in the number of degrees of freedom from the electroweak scale to zero temperature  @xcite .
hence , wdm3 should not be considered an unrealistic scenario .
in addition , wdm3 is of interest because it is consistent with the pulsar kick mechanism via off - resonant active - sterile neutrino oscillations  @xcite , even if this sterile neutrino is not the dominant component of dark matter .
the results of the calculation are presented in figures [ fig : t ] , [ fig : xe ] and [ fig : xh2 ] .
several conclusions can be drawn from these results .
we find that large clouds , with masses greater than @xmath117 , collapse regardless of the details of cooling in all the models defined above .
one could expect this behavior because the virial temperature of such a cloud is higher than @xmath12  k , so the h cooling mechanism is efficient enough and the additional h@xmath11 does not play an important role .
there is , however , a dramatic difference when it comes to smaller clouds . for a virialization redshift @xmath118 ,
a cloud of @xmath119 collapses in model wdm1 , but not in the cdm or other wdm models .
this illustrates the importance of the increased h@xmath120 fraction on the cooling of gas .    in wdm2 and wdm3 model ,
because of the relatively low associated x - ray flux , the temperature evolution is practically unaffected by the photons from the dark matter decay .    from fig .
[ fig : xh2 ] one can draw another conclusion . in the large clouds ,
all of the h@xmath11 molecules are destroyed during the virialization process .
this is because of the very high virial temperature .
in contrast , there is a drastic increase in h@xmath11 fraction during the virialization of the clouds with the mass of @xmath121 . in both cases ,
most of molecular hydrogen is produced after virialization .
for clouds with @xmath122 and @xmath123 the critical h@xmath11 fraction ( see eq .
( 11 ) of @xcite ) is reached at virialization . on the other hand ,
in the clouds having equivalent mass with @xmath118 , the molecular hydrogen faction @xmath106 attains its virial critical value only for wdm1 model .
this happens at @xmath124 .    in all the wdm models the ionization and the molecular hydrogen fraction
are enhanced in comparison with the cdm model , as suggested by @xcite .
in addition , we can see that in wdm models the temperature of the collapsing clouds soon after virialization is lower then that in the cdm case , which speeds up the collapse , at least for the low - mass clouds .
we note in passing that the density of the cloud was kept constant after virialization , so the evolution of @xmath52 , @xmath92 and @xmath106 is reliable for the redshifts close to virialization , but not for much smaller redshifts .
it appears that wdm3 model is disfavored .
the increase in the cooling rate is almost negligible in comparison with the cdm case , while the suppression on the small scales is very strong .
therefore , this leads to an unacceptable delay in the onset of star formation , as discussed by @xcite .    in the case of wdm1 and wdm2 models ,
the cooling is enhanced by the additional fraction of molecular hydrogen .
however , we can not definitively prove that star formation occurs at the redshift consistent with wmap because of the simplifying assumptions and limitations of our analysis .
the main limitation is that the virialization redshift of a particular cloud is a fixed parameter in the simulation . in reality , it can be changed by the enhancement of the cooling , so the cloud can reach its virial stage much earlier .
moreover , the density of the cloud is kept constant after virialization , which is a simplifying assumption .
we can not calculate the real redshift of the cloud collapse and onset of the star formation . as a result
, we can not calculate how many clouds collapse at a given redshift and we can not calculate the redshift of reionization in wdm model reliably .
also , we have not included the additional x - rays from the decay of sterile neutrinos inside the collapsing clump .
this radiation can cause an additional increase in ionization and heating , which could change our results .
also the effects of opacity can be important in large clouds , which we have not taken into account .
these limitations will be addressed in the future work .
we have performed a detailed analysis of the cooling and collapse of primordial gas in the model with warm dark matter , taking into account both the increase in the fraction of molecular hydrogen  @xcite and the heating due to the sterile neutrino decay . as expected , the effect on the largest gas clouds is negligible , but the smaller clouds are , in fact affected .
we have performed the analysis for some benchmark cases which arise in realistic scenarios . for the largest clouds , the additional molecular hydrogen made no difference . for smaller clouds , especially for those with masses
@xmath125 , the increase in the x - ray background made the successful collapse possible in cases where it could not occur in the absence of sterile neutrino decays .
we thank g.  gilmore , m.  mapelli , m.  shaposhnikov , and s.  stachniewicz for very helpful discussions and comments . the work of p.l.b . and
j.s . was supported by the pierre auger grant 05  cu  5pd1/2 via desy / bmbf .
the work of a.k . was supported in part by the doe grant de - fg03 - 91er40662 and by the nasa atp grants nag  5 - 10842 and nag  5 - 13399 .
thanks the cern theory unit and epfl for hospitality during his visit .
thanks the organizers of marcel grossmann meeting for support .
in this section we shall calculate the photon spectrum produced by the radiative decays of the sterile neutrinos ( our derivation is similar to the one in @xcite ) .
the dominant channel of radiative decays is through the following reaction @xmath126 where @xmath127 is the sterile neutrino and @xmath128 is an ordinary active neutrino .
the photon energy is given by @xmath18 where @xmath129 and @xmath19 are respectively the mass of the sterile and an active neutrino .
the cosmic expansion redshifts the photon energy .
a photon that at redshift @xmath33 has energy @xmath130 was produced at redshift @xmath131 given by @xmath132    let @xmath133 be the energy flux ( in units of erg @xmath134 @xmath135 ) at redshift z of photons with energy e produced by the sterile neutrino decays .
the flux per unit energy and solid angle is given by @xmath136 where @xmath137 is the photon flux at redshift z. it is related to the photon flux at emission redshift @xmath131 through the following equation @xmath138 where we have included the factor of @xmath139 produced by the expansion of the universe .
the photon density emitted at @xmath131 is given by the usual decay law @xmath140 where @xmath141 is the mean sterile neutrino lifetime for the decay into an active neutrino and a photon ( see eq .
[ meanlife ] ) , @xmath142 is the number density of the sterile neutrino at redshift @xmath131 and @xmath143 is given by @xmath144 where @xmath145 is the hubble expansion rate at photon emission redshift @xmath131 . writing everything together we obtain @xmath146    choosing the time @xmath147 otherwise arbitrary we can write @xmath148 where @xmath149 is the age of the universe at redshift @xmath131 , @xmath150 and @xmath151 would be the number density of sterile neutrinos at redshift @xmath33 if they did not decay .
assuming that all of the dark matter consists of sterile neutrinos , @xmath151 can be written as @xmath152 where @xmath153    in the case of a flat universe filled with matter and a nonvanishing cosmological constant the age of the universe at redshift @xmath33 is given by @xcite @xmath154    taking @xmath155 and substituting eq .
( [ neu ] ) into eq .
( [ ff1 ] ) we can get the specific intensity ( in units of erg @xmath134 @xmath135 @xmath156 @xmath157 ) of the photons with frequency @xmath158 from the sterile neutrino decays at redshift @xmath33 @xmath159    we have to multiply this value by @xmath160 in order to correct specific intensity for the absorption by the background neutral hydrogen atoms .
the optical depth is given by @xmath161 where @xmath162 and @xmath163 are the cross section for the h ionization at energy @xmath164 and the hydrogen number density of the universe at redshift @xmath165 respectively ; @xmath164 is the energy of the photon at redshift @xmath165 , which was emitted in sterile neutrino decay at redshift @xmath131 and is given by the equation @xmath166    the mean hydrogen number density of the universe at redshift @xmath165 can be written as @xmath167 \left(1+z'\right)^3 \textrm { , } \label{ap2}\ ] ] where @xmath168 is its present value and is equal to @xmath169    the cross section for the h ionization is given by @xmath170 where @xmath171 @xmath172 and @xmath173 ev .    neglecting the @xmath174 term in eq .
( [ hubble ] ) , which is a very good approximation at high redshifts , and using eqs .
( [ energy ] ) , ( [ ap1 ] ) , ( [ ap2 ] ) and ( [ ap4 ] ) we can cast eq .
( [ abs ] ) into the form @xmath175\left(1+z'\right)^{-5/2 }   dz ' \textrm{.}\ ] ] as @xmath176 , let us take @xmath177 where @xmath178 is the recombination redshift .
it gives us @xmath179 \textrm{.}\ ] ]    photons emitted before recombination are compton scattered and thermalized .
for simplicity we assume that these photons are not visible .
the final specific intensity is equal to @xmath180 where @xmath181 @xmath182    we have assumed that the photons are absorbed only by the ly@xmath32 wedge , so the above equations are valid only for @xmath183 ev . to derive the photon flux at energy lower than 13.6 ev we have to take into account another hydrogen absorption lines .
the radiation with energy below 13.6 ev start to be important only at lower redshifts ( e.g. for @xmath184 kev at @xmath185 ) so it is not relevant in our calculations .
the photons from the decays of the sterile neutrino are mainly absorbed by neutral hydrogen atoms leading to their ionization .
the ionization rate due to these photons is enhance almost 100 times due to additional ionization by the secondary electrons which deposit almost 1/3 of its energy into ionization ( depending on the ionization fraction of the medium ) .
the energy of the absorbed photons partially goes into ionization and partially into heating and excitations .
we have adopted the approximation by @xcite , for which ionization rate and heating due to the photons from the decays of the sterile neutrino are respectively equal to @xmath186 \nonumber \\ & \times & c_{i}\left(1-x_{e}^{a_{i}}\right)^b_{i } + \int_{\nu_{th}^{h}}^{\infty}4\pi\sigma_{h}\left(\nu\right)\frac{i_{\nu } \left(z\right)}{h\nu } d\nu \label{ion}\end{aligned}\ ] ] @xmath187 \nonumber \\ & \times & c_{h}\left [ 1-\left(1-x_{e}^{a_{h}}\right)^b_{h } \right ] n_{h } \label{heat } \textrm{,}\end{aligned}\ ] ] where @xmath188\ ] ] is the cross section for he ionization , @xmath189  ev is its energy threshold ; @xmath190 , @xmath191 , @xmath192 , @xmath193 , @xmath194 , @xmath195 @xcite .
the h@xmath11 cooling function ( @xmath196 , @xmath197 and @xmath198 ) given in table [ tab : cooling ] has been computed assuming that the cmb temperature is much smaller than the gas temperature .
this approximation is not valid in our case since in the wdm model with kev sterile neutrinos the gas temperature can be very close or even lower than the cmb temperature because of relatively high ionization fraction .
the situation is very similar in cdm model , where the gas temperature strictly follows the cmb temperature at least at high redshifts .    in such conditions ,
the level population of h@xmath11 is strongly affected by stimulated emission and absorption .
molecules become an effective heating source for the gas , because the rate of collisional de - excitation of the vibrational levels is faster than their radiative decay . in order to correct the cooling / heating function for this effect
, we have to take the net cooling rate @xmath199 where @xmath200 represents either @xmath196 , @xmath197 , or @xmath198 .
we also note that @xmath201 for @xmath202 .
one should do the same for the @xmath203 ( see table [ tab : cooling ] ) .
however , the cooling from the collisional excitation of hydrogen atoms h and its subsequent spontaneous de - excitation is only important for the gas temperatures @xmath204 . in our case
the cmb temperature is well below @xmath205 k , so we can neglect the heating term for this cooling / heating mechanism .
abazajian , k.  2006a , , 73 , 063506 abazajian , k.  2006b , , 73 , 063513 abazajian , k. , fuller , g. m. , & patel , m.  2001 , , 64 , 023501 abazajian , k. , fuller , g. m. , & tucker , w. h.  2001 , , 562 , 593 abazajian , k. , koushiappas , s. m. , & zentner , a. r.  2006 , private communication    abel , t. , anninos , p. , zhang , yu . , & norman , m. l.  1997 , , 2 , 181 akhmedov , e. k. , rubakov , v. a. , & smirnov , a. y.  1998 , , 81 , 1359 allgood , b. , flores , r.  a. , primack , j.  r. , kravtsov , a.  v. , wechsler , r.  h. , faltenbacher , a. , & bullock , j.  s.  2006 , , 367 , 1781    appelquist , t. & shrock , r.  2002 ,   phys .
b 548 , 204 asaka , t. , kusenko , a. , & shaposhnikov , m.  2006 , arxiv : hep - ph/0602150 asaka , t. , laine , m. , & shaposhnikov , m.  2006 , arxiv : hep - ph/0605209 asaka , t. , blanchet , s. , & shaposhnikov , m.  2005 , phys .
b , 631 , 151 asaka , t. , & shaposhnikov , m.  2005 , phys .
b , 620 , 17 athanassopoulos , c. , et al .
( lsnd )  1998 , , 81 , 1774 barkovich , m. , dolivo , j. c. , & montemayor , r.  2004 , , 70 , 043005 barkovich , m. , dolivo , j.  c. , & montemayor , r.  2005 , in _ new research _
, nova science publishers [ arxiv : hep - ph/0503113 ] belokurov , v. , et al .  2006 , , 642 , l137    biermann , p. l. , & kusenko , a.  2006 , , 96 , 091301 bode , p. , ostriker , j. p. , & turok , n.  2001 , , 556 , 93 van den bosch , f. c. , & swaters , r. a.  2001 , , 325 , 1017 boyarsky , a. , neronov , a. , ruchayskiy , o. , & shaposhnikov , m.  2005 , arxiv : astro - ph/0512509      boyarsky , a. , neronov , a. , ruchayskiy , o. , shaposhnikov , m. , & tkachev , i.  2006c , arxiv : astro - ph/0603660 bryan , g. l. , & norman , m. l.  1998 , , 495 , 80 dalcanton , j. j. , & hogan , c. j.  2001 , , 561 , 35 de gouva , a. , phys .  rev .
d 72 , 033005 ( 2005 ) .
dodelson , s. , & widrow , l. m.  1994 , , 72 , 17        fryer , c. l. , & kusenko , a.  2006 , .
, 163 , 335 fuller , g. m. , kusenko , a. , mocioiu , i. , & pascoli , s.  2003 , , 68 , 103002 galli , d. , & palla , f.  1998 , , 335 , 403 gell - mann , m. , ramond , p. , & slansky , r.  1979 , in _ supergravity _ , ed . by freedman , d. ,
et al . , north holland    gelmini , g. , palomares - ruiz , s. , & pascoli , s.  2004 , , 93 , 081302 gilmore , g. , in proceedings of _ 7th ucla symposium on sources and detection of dark matter and dark energy in the universe _ , 22 - 24 feb 2006 , marina de rey , california , to appear .
governato , f. , et al .
2004 , , 607 , 688 haiman , z. , thoul , a. a. , & loeb , a.  1996 , , 464 , 523 haiman , z. , & bryan ,  g.  l.  2006 , arxiv : astro - ph/0603541 hansen , s.  h. , lesgourgues , j. , pastor , s. , & silk , j.  2002 , , 333 , 544        kitayama , t. , susa , h. , umemura , m. , & ikeuchi , s.  2001 , , 326 , 1353 kormendy , j. , & fisher , d. b.  2005 , in revista mexicana de astronomia y astrofisica conference series , pp .
101 - 108 , arxiv : astro - ph/0507525 kusenko , a. , & segr , g.  1997 , phys .
b396 , 197 kusenko , a. , & segr , g.  1999 , , 59 , 061302 kusenko , a. , segr , g. , & vilenkin , a.  1998 , phys .
b437 , 359 kusenko , a.  2004 , ijmp d 13 , 2065 kusenko , a.  2006 , arxiv : hep - ph/0609081 jena , t. , et al .
2005 , , 361 , 70                        sommer - larsen , j. , & dolgov , a. d.  2001 , , 551 , 608 sorel , m. , conrad , j. m. , & shaevitz , m.  2004 , , 70 , 073004 spergel , d.n .
( wmap )  2003 , , 148 , 175 spergel , d.n . , et al .
( wmap )  2006 , arxiv : astro - ph/0603449 strigari , l.  e. , bullock , j.  s. , kaplinghat , m. , kravtsov , a.  v. , gnedin , o.  y. , abazajian , k. , & klypin , a.  a.  2006 , arxiv : astro - ph/0603775      viel , m. , lesgourgues , j. , haehnelt , m.  g. , matarrese , s. , & riotto , a.  2005 , , 71 , 063534 viel , m. , lesgourgues , j. , haehnelt , m.  g. , matarrese , s. , & riotto , a.  2006 , arxiv : astro - ph/0605706 watson , c.  r. , beacom , j.  f. , yuksel , h. & walker , t.  p.  2006
, arxiv : astro - ph/0605424 wilkinson , m.  i. , kleyna , j.  t. , wyn evans , n. , gilmore , g.  f. , read , j.  i. , koch , a. , grebel , e.  k. , & irwin , m.  j.  2006 , eas publications series , 20 , 105 willman , b. , governato , f. , wadsley , j. , & quinn , t.  2004 , , 353 , 639 yanagida , t.  1980 , prog .
phys . , 64 ,
1103      lll +   + @xmath209 & @xmath210 & @xcite + @xmath211 & @xmath212 & @xcite + @xmath213 & @xmath214 & @xcite +   +   +   + @xmath215 & @xmath216 & + & @xmath217 & @xcite + h@xmath218 & @xmath219 & @xcite + h@xmath220 & @xmath221 & @xcite +   +   +   + h@xmath222 & @xmath223 + & @xmath224 & @xcite + & @xmath225 $ ] & @xcite + & @xmath226 + & @xmath227 & @xcite + h@xmath228 & @xmath229 & @xcite + & @xmath230 + & @xmath231 + & @xmath232 + & @xmath233 + h@xmath234 & @xmath235 & @xcite +   +   +   + h@xmath236 & @xmath237 & @xcite + & @xmath238 is the cmb radiation flux in erg @xmath134 @xmath135 + & at @xmath239 & + h@xmath240 & @xmath241 & @xcite + & @xmath242 ev , + & @xmath243 is the cmb intensity & + & @xmath244 + & @xmath245 and @xmath246 are the crossections for para- and ortho - h@xmath11 + & respectively , @xmath247 is the ortho- , to para - h@xmath11 ratio +   +   +   + @xmath248 & @xmath249 & @xcite + & @xmath250 ev + & @xmath251 + @xmath252 & @xmath253 & @xcite + @xmath254 & @xmath255 & @xcite +   +   +   +   + @xmath248 & see text eq .
( [ ion ] ) & + [ tab : reaction ]    ll +   +   +   + a ) excitation of h by collisions with @xmath15 + @xmath256 & @xcite +   + b ) recombination to h + @xmath257 & @xcite +   + c ) ionization of h by collisions with @xmath15 + @xmath258 & @xcite
+   + d ) photo - ionization of h + @xmath259 + @xmath260 ev +   +   +   +   +   + a ) collisional excitation of h@xmath11 +   +   + @xmath261 & @xcite +   +   +   + @xmath262 & @xcite +   + @xmath263 & @xcite + @xmath264 & @xcite + @xmath265 +    3\times10^{-24}\exp\left(-\frac{0.51}{t_{3}}\right)\right\}$ ] & @xcite + @xmath266 $ ] & @xcite + @xmath267 + 0.75\left[\frac{7}{3 } \gamma_{3}^{h_{2}}\exp\left[-\left(\frac{e^{r}_{31}}{kt}\right)\right]e^{r}_{31 } \right]$ ] & @xcite + @xmath268 & @xcite + @xmath269\right\}$ ] & @xcite + @xmath270 k , @xmath271 k & @xcite +   + b ) h@xmath11 formation cooling @xmath272 + @xmath273n_{h}n_{e}$ ] & @xcite + @xmath274 + 1.4x_{h_{2}}\exp\left[-\frac{12000}{t+1200}\right]\right\}$ ] & @xcite +   + c ) collisional dissociation of h@xmath11 + @xmath275 & @xcite +   + d ) photo - dissociation of h@xmath11 +   +   + @xmath276 & @xcite +   +   + @xmath277 $ ] & @xcite + @xmath278 ev , @xmath279 ev , @xmath280 is the ortho - h@xmath11 to para - h@xmath11 ratio +   +   + @xmath281 & @xcite +   +   +   +   +   + @xmath282 & @xcite + @xmath283 $ ] +   +   +   +   +   + @xmath284see text eq .
( [ heat ] ) & + [ tab : cooling ]