the importance of neutrinos in the explosions of massive stars has long been recognized .
this significance is coupled with their ability to carry information to us about the processes and conditions in the core of the supernova so that the neutrino signal from the next galactic supernova will provide us with an opportunity to test the core collapse paradigm .
the explosion begins when the runaway process of electron capture in the core a massive star leads to a rapid compression and collapse that is only halted when the degenerate pressure of the neutrons and thermal pressure of non - degenerate particles kicks in at super - nuclear densities .
the rapid neutronization during the collapse leads to a large burst of @xmath0 neutrinos that identifies the beginning of the event . as nuclear densities
are reached the mean free path of the neutrinos becomes shorter than the size of the proto - neutron star and the neutrinos become trapped .
thermal processes within the core create a thermal bath of neutrino - antineutrino pairs of all flavors .
the neutrinos slowly diffuse from the core over a period of order @xmath1 seconds , carrying 99% of the gravitational binding energy of the core with them .
with such a large neutrino luminosity even a small number of neutrino interactions above the core can create an important impact on the explosion physics .
detection of the neutrinos from the next nearby supernova will offer us the opportunity to examine the internal evolution of the explosion since the neutrinos allow us to see down all the way to the proto - neutron star .
for example , if the proto - neutron star collapses to a black hole , the resulting neutrino signal will be altered @xcite .
in addition , although these neutrinos are only weakly interacting they can have significant effects upon the nucleosynthesis that occurs in the supernova : in the outer layers of the star neutrinos can cause a transformation of the elements synthesized during the preceding thermonuclear burning period in the `` neutrino process '' @xcite , it is thought that a neutrino driven wind occurs at late time in the supernova which may create sufficiently neutron rich material to produce the r - process elements @xcite , and long duration gamma ray bursts ( which are thought to be a rare type of supernova ) produce elements in conditions governed by strong neutrino fluxes @xcite .
the results of all these nucleosynthesis processes are quite sensitive to the neutrino spectra @xcite .    for reasons pertaining to both hydrodynamics and to element synthesis
it is important to understand the details of neutrino propagation in the supernova and the mixing between neutrino flavors .
neutrino oscillations also alter what we observe but , simultaneously , also allow the possibly of elucidating yet unknown neutrino oscillation physics , e.g. @xcite .
supernova neutrino flavor transformation is a rapidly developing field with the neutrino background terms the subject of intense study @xcite .
realistic possibilities include a complete mixing of all flavors @xcite or a partial oscillation between flavors @xcite .
an important part of any future observation of supernova neutrinos will be to understand this physics . in parallel , phase effects due to multiple resonances
have recently been recognized @xcite distorting the simple picture of neutrino propagation through the outer layers of the supernova that has been used in the past @xcite .
the explosion of the star leads to an evolution of the flavor mixing so that to appreciate the effects upon neutrino propagation it is important to understand the density profile .
observations of high pulsar velocities and polarized supernova light suggest an inherent asymmetry in the explosion mechanism .
though many possible mechanisms for generating the asymmetry have been proposed , recent work by blondin , mezzacappa & demarino @xcite has identified an instability of the standing accretion shock that leads to large dipole and quadrupole moments .
more recent work by blondin & mezzacappa @xcite and ohnishi , kotake , and yamada @xcite as well as scheck _
@xcite has both confirmed and furthered the understanding of this instability .
an alternative mechanism , acoustic heating , has also been discussed by burrows et al . @xcite .
the temporal evolution of the neutrino flavor mixing and the subsequent variation of the neutrino signal means that it may be possible to to detect prominent features of the supernova density profile and to use this information to learn important information about the explosion .
this possibility was demonstrated first by schirato and fuller @xcite who used a time - dependent , one - dimensional @xmath2 density profile to demonstrate that the supernova s forward shock wave reaches and disrupts the ` h ' resonance transformation layer a few seconds after the core bounce leading to a detectable change in the charged current neutrino signal .
this was followed and elaborated upon by dighe & smirnov @xcite , takahashi _ et al .
_ @xcite and fogli _ et al .
more recently tom ` as _
@xcite used a more sophisticated numerical model of a progenitor star to investigate the effect of both forward and reverse shocks on the neutrino signal .
tom ` as _ et al .
_ also present one second of a two - dimensional simulation result to show how the density profile can be greatly complicated by a deformed forward shock followed by strong convection currents . but
due to the intensive computational burden of supernova models it is numerically very taxing to to watch these the explosion propagate outward until they affect neutrino transformations .    in this paper
we present calculations which link hydrodynamical simulations of time dependent density profiles and phase - retaining neutrino oscillations .
our presentation begins in section [ section : simulations ] with the results from simulations of supernovae in both one and two dimensions where we artificially heat a density profile constructed so as to mimic the state of the supernova at the point where the accretion shock has stalled .
our two - dimensional explosions are heated aspherically so as to create an aspherical supernova .
though this 2d model lacks some structural features present in more sophisticated models that generate asphericity via standing accretion shock instabilities , we still obtain results that give a good idea of how density isotropy affects neutrino flavor transformations . in section iii
we pass a spectrum of neutrinos through the simulation results and show how the crossing probability is affected by the evolution of the profiles .
we finish in section iv by presenting our calculation of two neutrino detector signals - the positron spectra within a water cerenkov detector , and the ratio of charged current to neutral current event rates for a heavy water detector - that demonstrate our ability to extract information about the supernova .
our first task in determining the neutrino signal is to simulate the supernova so as to obtain the time evolution of the density profile . a supernova simulation is a complex and computationally intensive problem . at the present time
there is not yet a robust , self - consistent model for core - collapse supernovae , and the most sophisticated multi - dimensional simulations take a very long time to compute .
but we do not need to simulate the actual core collapse and formation of the proto - neutron star , these events occur deep within the core at a radius of @xmath3 or so , whereas the first neutrino resonance occurs at around @xmath4 for neutrino energies of order @xmath5 , and the second is even further out
. what is occurring in the core is irrelevant for our purposes , all we need from the simulation of the supernova is the propagation of the shock through the star .
for this reason we employ a simplistic and artificially energetic numerical simulation to create explosions with a range of features that might be present in reality .
we do not intend this model to be realistic itself , but rather a tool to help with the difficult problem of reconstructing density features of an exploding star based on changing neutrino signals detected at earth .
we use vh-1 @xcite , a hydrodynamic code based on woodward and colella s @xcite piecewise parabolic method . for both the one - dimensional and two - dimensional simulations we employ an exponentially coarsening radial grid that ranges from @xmath6 to @xmath7 where @xmath8 is the radius from the center of the star , and for the two - dimensional simulations an angular grid that ranges from @xmath9 to @xmath10 where @xmath11 is the polar angle .
the innermost radial grid increment is set to @xmath12 and increases by 1.3% thereafter .
by @xmath13 the radial increment has grown to @xmath14 and by @xmath15 we reach @xmath16 .
we map into the code a spherically symmetric progenitor profile intended to represent the state of the supernova about @xmath17 after the collapse .
beyond @xmath18 the profile is that of a low metallicity , @xmath19 progenitor model developed by heger et al @xcite that we have allowed to collapse further within vh-1 so that the infall velocities approached free fall .
interior to @xmath20 we splice in a slow , outward - moving standing accretion shock profile provided by heywood @xcite that is similar to that found in @xcite .
finally , inside @xmath21 we insert a dense , @xmath22 , polytropic core containing approximately @xmath23 .
the inner boundary condition inside the dense core at @xmath24 is reflecting .
figure ( [ fig:1da ] ) shows the full initial density profile .
the initial , @xmath25 , density profile used in our sn simulations .
the dense core ( a ) inside @xmath26 contains approximately @xmath23 .
the slow , outward - moving standing accretion shock ( b ) is located at @xmath27 . above that we have a collapsing , @xmath19 , progenitor ( c).,width=480 ]    to drive the explosion , we mimic neutrino heating by inserting energy into the region above a @xmath28 gain radius .
the energy deposition rate per unit volume , @xmath29 , is proportional to the density of material , falls as @xmath30 and decreases exponentially with time @xmath31 over a timescale @xmath32 set to @xmath33 .
for the 1d simulations then , @xmath34 the additional factors of @xmath35 and @xmath31 are inserted so as to ramp up the energy deposition over both distance and time thereby avoiding the violent disruption that occurs if we deposit energy either too suddenly or within a small volume .
for the 2d models we introduce an angular dependence into the energy deposition prescription so as to match the observation by blondin _
_ @xcite that small perturbations in standing accretion shock models can ultimately lead to aspherical shock modes .
we break the spherical symmetry of our 2d simulations by heating the initial profile with a combination of 75% spherical mode heating and 25% @xmath36 mode heating .
thus for the 2d cases @xmath29 becomes @xmath37 in either case the total energy , @xmath38 , input during the simulation was recorded .
below the gain radius we wish to maintain the spherical density / gravity conditions and to prevent mass from escaping the core so that we do not disturb the shock heating and convective flow in the outer star . in a real supernova the core is stabilized by neutrino emission but in lieu of full implementation of this cooling mechanism our simulations achieve the same result by forcing all radial and angular velocities below @xmath39 to zero . in this way we separate the evolution of the outer layers of the star from the complex behavior near the core .
the runtime of the 1d simulations is sufficiently short that we can vary the total energy deposition thereby obtaining a variety of results ranging from weak explosions through to the very powerful .
the 2d simulations take much longer to run and so we have just one data set with an explosion energy set to @xmath40 .    before presenting our simulation results
we mention that in general , numerical schemes tend to spread the shock front over several zones and they therefore become artificially softened .
schirato and fuller @xcite account for this artifact by steepening by hand the density profile at the shocks .
the shocks in our simulations are sufficiently steep to demonstrate the primary effect on neutrino mixing , therefore we do not implement a similar correction .
the density profiles of the sn we obtain from the simulations possess noticeable differences as the energy deposition changes . for weak explosions the profiles are monotonically decreasing functions of the radius , @xmath8 , with a single , forward facing , shock front that moves relatively slowly outwards . as the deposition energy increases a lower - density cavity forms behind the shock .
further increases in @xmath38 eventually lead to the formation of a reverse shock behind the lower - density zone .
all three features are present in the results of the 2d simulations but the higher dimensionality , coupled with the aspherical heating means , that this sn possesses a much more turbulent / chaotic profile .     the density as a function of the radius and angle in a 2d sn model at @xmath41 .
the forward shock is located to the left of ` a ' , the reverse shock is the step - up in density found to the right of ` b ' , and one of the many local cavities in the profile between the shocks is to the right of ` c'.,width=480 ]    in figure ( [ fig : rho.snv2d.t=2.5 ] ) we show a time slice of the 2d sn simulation .
the forward shock is seen at ` a ' , the ` reverse shock ' is at ` b ' and one of the many bubbles in the density profile is located at ` c ' .
we note that although our applied neutrino heating is symmetric about the equator ( @xmath42 ) the resulting density profile is not .
we attribute this to the inherent instability of multi - dimensional accretion shocks , potentially excited by simulation numerics . in what follows
we discuss the behavior and properties of the various features of the profile .
the forward shock is a generic feature of supernovae simulations .
a stalled forward shock was present in the initial profile inserted into the hydrodynamical code and the heating we introduced was meant to revive it s outward motion . after revival the shock propagates out through the star and is the explosion feature furthest from the proto - neutron star .
the forward shock is visible in figure ( [ fig : rho.snv2d.t=2.5 ] ) as the large jump in density at larger radii and we note that the forward shock in the 2d simulation is slightly oblate due to the aspherical heating of the material .
we study the behavior of the forward shock with various explosion strengths using the 1d simulations .     the density as a function of the radius in a 1d sn model with @xmath43 at @xmath44 ( dotted ) , @xmath45 ( solid ) , @xmath46 ( long dashed ) and @xmath47 ( dash dot).,width=480 ]    in fig .
( [ fig : rho.snv1d.1e35.t=0.9,1.8,3.6,7.2 ] ) we show four snapshots of the density profile in a one - dimensional simulation in which the total energy deposition was @xmath43 .
actual sn are thought to be more energetic than this so we regard the results from this simulation as being at the lower end of realistic possibilities .
this energy is also significantly less than that used in the simulation shown in fig . (
[ fig : rho.snv2d.t=2.5 ] ) and , consequently , the profile is much simpler .
the forward shock is clearly visible and we note that the fractional jump in density across it does not vary . for normal shocks
the density jump across the shock is given by @xmath48 where @xmath49 the mach number and @xmath50 is determined by the ratio of specific heats .
in our simulations we adopted a fixed equation of state and the mach number of the shock is always sufficiently large that we find that the density jump does not vary to any great extent with time and/or with @xmath38 .
but the position of the forward shock as a function of time is greatly influenced by @xmath38 with larger values leading to more rapidly moving forward shocks .
the density as a function of the radius in a 1d sn model with @xmath51 at @xmath44 ( dotted ) , @xmath45 ( solid ) , @xmath46 ( long dashed ) and @xmath47 ( dash dot).,width=480 ]    this can be seen in fig .
( [ fig : rho.snv1d.6e35.t=0.9,1.8,3.6,7.2 ] ) which are snapshots at the same moments as those in fig .
( [ fig : rho.snv1d.1e35.t=0.9,1.8,3.6,7.2 ] ) for a 1d simulation with @xmath51 .
this figure also displays the low - density cavity that can form behind the forward shock as @xmath38 increases .
this profile is very similar to the profiles used by fuller & schirato @xcite and fogli _ et al .
_ @xcite .
the heating that led to the regeneration of the forward shock continues to accelerate the material above the proto - neutron star even after the shock has been revived and is moving outwards .
a wind is created with a velocity that increases with radius . when the velocity of the material becomes larger than the local sound speed a reverse shock
is formed .
the reverse shock feature was not present in the ` initial ' profile : it develops only later . if a reverse shock is to form enough energy must be deposited to create a sufficiently strong wind .
a reverse shock can be seen in figure ( [ fig : rho.snv2d.t=2.5 ] ) where it is the step - up in density at the back of the turbulent zone behind the forward shock .
the reverse shock is oblate due to the aspherical heating of the material in the simulation .
the density as a function of the radius in a 1d sn model where @xmath52 at @xmath53 ( dotted ) , @xmath54 ( solid ) , @xmath55 ( long dashed ) , @xmath41 ( dash - dot ) and @xmath56 ( dash double dot).,width=480 ]    once again we can use 1d simulations to study this feature more easily .
no reverse shock was seen in the profiles shown figure ( [ fig : rho.snv1d.1e35.t=0.9,1.8,3.6,7.2 ] ) or ( [ fig : rho.snv1d.6e35.t=0.9,1.8,3.6,7.2 ] ) but when we increase @xmath38 to @xmath52 we obtain the results presented in figure ( [ fig : rho.snv1d.8e35.t=1,1.5,2,2.5,3 ] ) .
we notice that the forward shock is , again , moving more rapidly compared to the results shown in figures ( [ fig : rho.snv1d.1e35.t=0.9,1.8,3.6,7.2 ] ) and ( [ fig : rho.snv1d.6e35.t=0.9,1.8,3.6,7.2 ] ) due to the larger energy deposition @xmath38 .
the figure also indicates that the reverse shock is smaller than the forward shock . as @xmath38 increases further the density jump across the reverse shock increases and it moves closer to the forward shock .
we find that it is also possible for the reverse shock to penetrate to densities lower than the forward shock .
the density as a function of the radius in a 1d sn model with @xmath57 at @xmath58 ( dashed ) , @xmath54 ( solid ) and @xmath59 ( long dashed).,width=480 ]    both of these behaviors are shown in figure ( [ fig : rho.snv1d.2e36.t=1,1.5,2 ] ) where we show snapshots of the profile for a 1d simulation with @xmath57 .
the density jump across the reverse shock is still smaller than across the forward but they are almost equivalent , and the density immediately behind the reverse shock is lower than the density immediately in front of the forward shock .
it was mentioned by toms _
_ @xcite that this feature was also seen at early times in their simulations . but realistically sn may struggle to achieve explosions with comparable energies so we regard the results of this simulation as being at the upper end of possibilities .    in both figures ( [ fig : rho.snv1d.8e35.t=1,1.5,2,2.5,3 ] ) and ( [ fig : rho.snv1d.2e36.t=1,1.5,2 ] )
the reverse shock was driven radially outward by a wind generated by the material heated above the proto - neutron star .
our energy deposition decreased exponentially with time which led to a gradual decrease in wind strength . in turn , as the wind abates , the outward motion of the reverse shock slows and , eventually , we find that its motion can be completely halted .
this stalling of the reverse shock also occurred in the two - dimensional simulation with the shock stalling at different times for different polar angles .
we also find that in both 1d and 2d simulations the density jump across the reverse shock decreases as the shock is about to turn around . after the reverse shock stalls both the 1d and
2d simulations indicate that the reverse shock feature then moves back towards the core . in figure ( [ fig : rho.snv1d.8e35.t=1,1.5,2,2.5,3 ] ) we see this stalling of the reverse shock and the backwards motion can be seen after comparing the profiles at @xmath41 and @xmath56 .
the reverse shock in this simulation actually reached it s furthest radial position at about @xmath60 .
similar backwards motion for the reverse shock may be seen in figure ( 1 ) of toms _ et al .
_ @xcite .
eventually the reverse shock reaches the core whereupon our simulations indicate that it is reflected and subsequently becomes a weak forward shock .
the density as a function of the radius in the same 1d sn model shown in figure ( [ fig : rho.snv1d.8e35.t=1,1.5,2,2.5,3 ] ) i.e. with @xmath52 , at @xmath61 ( dotted ) , @xmath62 ( solid ) and @xmath63 ( long dashed).,width=480 ]    further snapshots of the density profile taken from the simulation with @xmath52 are shown in figure ( [ fig : rho.snv1d.8e35.t=2.8,3.1,4.3 ] ) where we see more clearly the backwards motion of the reverse shock and its later reflection .
the radial position at which the reverse shock stalls depends upon the energy deposition .
additional snapshots from the simulation with slightly larger energy deposition , @xmath57 , are show in figure ( [ fig : rho.snv1d.2e36.t=4,4.5,4.9 ] ) .
though the reverse shock in this simulation also attained it s maximum radial position at @xmath60 the figure shows that it was located further out before it was turned around .     the density as a function of the radius in the same 1d sn model shown in figure ( [ fig : rho.snv1d.2e36.t=1,1.5,2 ] ) , that is @xmath57 , at @xmath64 ( dotted ) , @xmath65 ( solid ) and @xmath66 ( long dashed).,width=480 ]
the reverse shock is an interesting feature of the sn and , in contrast with the forward shock , can move both outwards and inwards .
for very energetic explosion it may penetrate to lower densities than that immediately in front of the forward shock .
the density jump across the reverse shock can vary with time and does so most noticeably when the shock is about to stall .
the difference in the hydrodynamics between the one and two dimensional sn simulations can be traced back to the aspherical heating of the material above the proto - neutron star in the two dimensional case .
non - spherical heating can lead to turbulent fluid flow creating eddies and bubbles as shown in figure ( [ fig : rho.snv2d.t=2.5 ] ) .     the density at a polar angle of @xmath67 as a function of the radius in a 2d sn model at @xmath68 ( solid ) , @xmath69 ( dotted ) and @xmath70 ( dashed).,width=480 ]     the density as a function of the radius in a 2d sn model at @xmath41 .
the angular slices are at @xmath71 ( solid ) , @xmath72 ( dotted ) and @xmath73 ( dash dot).,width=480 ]    the neutrinos released by the proto - neutron star propagate along radial slices of the profile . in fig .
( [ fig : rho.snv2d.t=3.9,4.8,5.7.slice=25 ] ) we show the density profile from the 2-d model at the polar angle of @xmath74 for snapshots at @xmath68 , @xmath69 and @xmath70 while in figure ( [ fig : rho.snv2d.t=2.5.slice=45,105,165 ] ) we present the density along three different polar angles at @xmath41 , i.e. taken from figure ( [ fig : rho.snv2d.t=2.5 ] ) . in
both we see the forward and reverse shocks and , as mentioned earlier , figure ( [ fig : rho.snv2d.t=2.5.slice=45,105,165 ] ) also shows that the radial position of both the forward and reverse shocks varies with the polar angle .
figure ( [ fig : rho.snv2d.t=3.9,4.8,5.7.slice=25 ] ) also indicates that between @xmath69 and @xmath70 the reverse shock ( along this radial slice ) stalled and began to move back to the core .     the density at a polar angle of @xmath75 as a function of the radius in a 2d sn model at @xmath65 ( dash dot ) , @xmath76 ( solid ) , @xmath77 ( dotted ) , and @xmath47 ( dashed).,width=480 ]    for other radial slices the reverse shock turnaround time will be different : in figure ( [ fig : rho.snv2d.t=4.5,5.4,6.3,7.2.slice=125 ] ) we plot radial slices along the @xmath75 line of sight at various times .
the reverse shock in the figure , located just beyond @xmath78 , turns around between @xmath79 and @xmath80 .    in all our 2d figures the general shape of the density profile between the two shocks is reminiscent of that shown in figure ( [ fig : rho.snv1d.8e35.t=1,1.5,2,2.5,3 ] ) for the 1d simulation with a similar value for @xmath38 - i.e. a lower density region in front of the reverse shock , a wall of higher density material behind the forward shock - but
clearly there are large ` fluctuations ' upon this basic trend both between the shocks and behind the reverse shock .      the forward shock , the reverse shocks , contact discontinuities and the local bubbles / cavities are the features in the density profile with the greatest potential to alter the state of any neutrinos propagating through the sn since they represent the locations within the sn where the density gradient is largest .
large density gradients lead to non - adiabatic evolution of the neutrinos and significant differences compared to the neutrino propagation through the undisturbed profile .
the vast majority of the neutrinos emitted by the supernova last interact with matter at the neutrinosphere located at the surface of the proto - neutron star . as they propagate outwards a small percentage
will be absorbed by nucleons thereby transferring energy to the supernova and , it is thought , reviving the stalled shock .
but even though the remainder of the neutrinos survive the passage through the material overlying the proto - neutron star what emerges is not the same as what was emitted .
the change in the neutrinos is due to neutrino oscillations and the presence of matter modulates this mixing .
neutrino oscillations arise due to a distinction between the interaction eigenstates ( otherwise known as the flavor eigenstates @xmath81 ) and the eigenstates of the free hamiltonian ( known as the mass eigenstates ) with masses @xmath82 .
since there are three basis states a general neutrino wavefunction is described by three complex coefficients and evolves according to the schrodinger equation . in the vacuum
the hamiltonian is diagonal in the mass basis but possesses off - diagonal terms in the flavor basis that are the cause of flavor oscillations . in the presence of matter a potential , @xmath83 , that takes into account coherent forward scattering of the neutrinos , is included in the hamiltonian . since we are only concerned with mixing between active neutrino flavors ( i.e. all the flavors that have ordinary weak interactions ) we may subtract off the common neutral current contribution ( which will contribute only an overall phase ) leaving just the charged current contribution to the @xmath84 component of @xmath83 .
this contribution is the well - known @xmath85 where @xmath86 is fermi s constant , @xmath87 is the mass density , @xmath88 is the electron fraction and @xmath89 is the nucleon mass .
the effect of the matter upon the antineutrinos differs from that of the neutrinos - the potential @xmath90 that appears in the antineutrino hamiltonian has the same magnitude but the opposite sign . due to the inclusion of @xmath91
the hamiltonian is a function of position and is neither diagonal in the mass basis nor the flavor basis .
one may try and diagonalize @xmath92 but the unitary transformation that relates the flavor basis to the new ` matter ' basis also varies with the position .
consequently the gradient of the unitary transformation is non - zero and one finds that the schrodinger equation in this new matter basis - that was meant to diagonalize @xmath93 - picks up off - diagonal terms .
thus it is not possible , in general , to find a basis that diagonalizes the hamiltonian and therefore oscillations of the amplitudes that describe the wavefunction occur in every basis .    in general , the three complex components of the wavefunction oscillate simultaneously .
exactly what occurs depends upon the energy of the neutrino @xmath94 , the differences between the squares of the masses @xmath82 , the mixing angles that describe the relationship between the flavor and mass basis , and , of course , @xmath91 .
since there are three mass eigenstates there are three separate mass splittings @xmath95 ( though only two are independent ) and the relationship between the neutrino flavor and mass bases is described by the matrix @xmath96 which , in turn , is parameterized by three mixing angles @xmath97 and @xmath98 plus a cp - phase @xmath99 .
the structure of @xmath96 is @xmath100 where @xmath101 , @xmath102 .
mixing has been observed in the neutrinos emitted by the sun and the neutrinos produced by cosmic rays striking the atmosphere . both have been confirmed with terrestrial experiments .
each observation of neutrino mixing can be described by a single @xmath103 - @xmath11 pair of parameters and experimentally the ` solar ' and ` atmospheric ' mass splittings differ by around a factor of @xmath104 with the solar mass splitting , @xmath105 , being the smaller of the two .
this observation permits us to consider the evolution of the general , three component , neutrino wavefunction as being factored into spatially distinct , two - neutrino mixes .
the factorization simplifies matters greatly and from it one can demonstrate that there are two resonances in the supernova density profile : the so - called ` l ' resonance and the ` h ' resonance .
the l resonance , at lower density and thus further from the proto - neutron star , involves mixing between matter states @xmath106 and @xmath107 , the relevant mass splitting is @xmath108 - which is approximately the solar mass splitting @xmath105 - and the mixing angle is @xmath109 - which is approximately the mixing angle determined by the solar neutrino experiments @xmath110 .
the relevant mass eigenstates and mass splitting for the h resonance , at higher density and closer to the core , depend upon the sign of @xmath111 ( or , equivalently , @xmath112 ) and this is not currently known .
if @xmath111 is positive ( a normal hierarchy ) then the h resonance involves mixing between states @xmath107 and @xmath113 and the mass splitting @xmath111 .
if @xmath111 is negative ( an inverted hierarchy ) then at the h resonance it is the antineutrinos states @xmath114 and @xmath115 that mix and the relevant mass splitting is @xmath116 . in either case the mixing angle is the unknown , but small , @xmath117 .
quite generally the coherent matter basis wavefunction that arrives at a distance , @xmath118 , from the core of the sn is related to the initial matter basis wavefunction at the proto - neutron star via the equation @xmath119 while for the anti - matter states @xmath120 .
the matrices @xmath121 and @xmath122 are the s - matrices for the neutrinos and antineutrinos respectively , for a discussion of this approach to neutrino oscillations , see e.g. @xcite . to determine @xmath123 and @xmath124 we need to know both the initial states @xmath125 and @xmath126 and the two matrices @xmath121 and @xmath122 .
the matrices @xmath127 and @xmath128 can be factored as @xmath129 and @xmath130 which are sufficiently general to accommodate our lack of knowledge of the hierarchy .
the two matrices @xmath131 and @xmath132 represent the neutrino propagation from the surface of the sn through the vacuum to earth and both are diagonal in the mass basis .
the matrices @xmath133 , @xmath134 and @xmath135 represent the change to the initial wavefunction due to the neutrino s or antineutrino s passage through the h or l resonance have multiple h resonances therefore the matrices @xmath134 and @xmath135 represent their combined effect of all the h resonances and @xmath133 all the l resonances .
we assume that the neutrinos encounter all l - type resonances after all h - type resonances .
this is the case for all profiles we have generated ] .
there is no matrix @xmath136 because antineutrinos do not experience an l resonance in matter .
the two matrices @xmath134 and @xmath135 are evaluated somewhere between the h and l resonances while @xmath137 is evaluated at the surface of the sn .
these evaluation positions will not affect our result if they are sufficiently far from the resonances .
the structures of @xmath133 , @xmath134 and @xmath135 follow from our knowledge of the matter states that mix at either the h or l resonance and are thus @xmath138 after omitting irrelevant phases . from these matrices
we define @xmath139 which are the crossing probabilities for neutrinos or antineutrinos at the two resonances .
the resonances are said to be ` adiabatic ' or ` non - adiabatic ' depending upon whether the crossing probability is close to zero or closer to unity .
indeed these are the two natural values since the crossing probability is determined by the ratio of the resonance width to the local oscillation lengthscale and typically one is significantly larger than the other .
our interest now turns to the initial states .
the density at the proto - neutron star is so large that the matter eigenstates and the flavor eigenstates are strongly aligned there . from a full 3-neutrino mixing calculation
we find that the initial matter basis spectra for a normal hierarchy ( nh ) are @xmath140 , @xmath141 while all other states , @xmath142 and @xmath143 are equal to the @xmath144 spectrum which we call @xmath145 . with an inverted hierarchy ( ih )
the initial states are @xmath146 , @xmath147 and this time @xmath148 and @xmath149 are all equal to @xmath145 .    putting the initial spectra and definitions of the s - matrices together one finds that the flux of matter state @xmath150 a distance d from the supernova is given by @xmath151 where @xmath152 and @xmath153 are the initial spectra of the matter states .
but what are detected at earth are , of course , the flavor states . during their flight from the supernova
to earth any coherence between the matter eigenstates is lost so the flux of flavor @xmath154 that arrives at earth is the incoherent sum @xmath155 these flavor fluxes may be rewritten in terms of the fluxes emitted by the neutrinosphere by introducing the survival probabilities @xmath156 and @xmath157 for the electron neutrinos and antineutrinos respectively since these are the two flavors with distinct initial spectra . written this way the detectable flavor fluxes are @xmath158 , \\   & & f_{\bar{\nu}_e}(d , e ) = \frac{1}{4\,\pi\,d^{2}}\;\left[\,\bar{p}(e)\,\phi_{\bar{\nu}_e}(0,e ) + ( 1-\bar{p}(e))\,\phi_{\nu_x}(0,e)\right ] , \\ & & 4f_{\nu_x}(d , e ) = \frac{1}{4\,\pi\,d^{2}}\;\left [ ( 1-p(e))\,\phi_{\nu_e}(0,e ) + ( 1-\bar{p}(e))\,\phi_{\bar{\nu}_e}(0,e ) + ( 2+p(e)+\bar{p}(e))\,\phi_{\nu_x}(0,e)\right]\nonumber\\ & & \end{aligned}\ ] ] where we have adopted the notation of dighe & smirnov @xcite by denoting @xmath159 by @xmath160 . from equations ( [ eq : fi ] ) and ( [ eq : fbari ] ) , and the initial spectra , we find that @xmath156 and @xmath157 are related to the elements of @xmath161 and @xmath162 via @xmath163 where @xmath164 , @xmath165 for a normal hierarchy , and @xmath166 , @xmath167 for an inverted hierarchy . after inserting the definition of the crossing probabilities @xmath156 and @xmath157 become @xmath168 these formulae are quite general but from our knowledge of the various neutrino oscillation parameters we can be more specific about what exactly happens to the neutrinos and antineutrinos .
the solar mixing parameters are such that the neutrino l resonance is almost completely adiabatic until , perhaps , the very end of the supernova neutrino signal .
for this reason @xmath169 .
thus at the end we find that determining the neutrino fluxes at earth comes down to computing @xmath170 or @xmath171 .
the schrodinger equation forms a starting point by which @xmath170 or @xmath171 can be determined .
but if one navely applies , for example , a runge - kutta integrator to this equation one quickly discovers that this is a difficult problem from a numerical standpoint because such algorithms are simply not suited to differential equations where the solution is a highly oscillatory function .
a number of alternate methods have been developed for calculating @xmath170 or @xmath171 , such as application of the the landau - zener result or the semi - analytic method by balantekin & beacom @xcite , but for one reason or another these alternate approaches can break down for complex profiles .
some authors @xcite have recognized that the evolution of the neutrino wave can be recast as a scattering problem and we adopt in this paper the algorithm outlined in kneller & mclaughlin @xcite which computes the evolution of the neutrino wavefunction by a monte carlo integration . for this paper we selected @xmath172 and for @xmath117 we use @xmath173 as a representative value . with this choice of @xmath117
the unperturbed profile is just sufficiently adiabatic to suppress the transformation effects of all but the most significant density features .
the current experimental limit on @xmath117 is set by the chooz experiment @xcite , and for our choice of @xmath103 that limit is @xmath174 .
finally , in addition to the effects caused by the propagation through the supernova there are potentially earth matter effects that can arise .
these are straightforward to compute and we do not consider them here since their presence ( or absence ) in the neutrino signal is a function of the position of the supernova with respect to the detector when the event occurs .
with the sn simulations complete and a method for calculating the crossing probabilities in hand we turn to study the effects of various features in the profiles upon @xmath175 both as a function of time and energy .
the forward shock is a generic feature of all our sn simulations .
initially the shock is located close to the core and at high density and then , after it is regenerated , begins to move outwards and to lower densities as shown in the figures from section  [ section : simulations ] .
the density as a function of the radius in a 1d sn model with @xmath43 at @xmath44 ( dotted ) , @xmath45 ( solid ) , @xmath46 ( long dashed ) and @xmath47 ( dash dot ) .
the horizontal dashed lines are ( from top to bottom ) the resonance densities for @xmath176 , @xmath1 , @xmath177 , @xmath178 and @xmath179 neutrinos.,width=480 ]    in fig .
( [ fig : rho.snv1d.1e35.t=0.9,1.8,3.6,7.2.+rhores ] ) we reproduce the results from fig .
( [ fig : rho.snv1d.1e35.t=0.9,1.8,3.6,7.2 ] ) and also superimpose the resonance densities for @xmath176 , @xmath1 , @xmath177 , @xmath178 and @xmath179 neutrinos .
the radii where the horizontal dashed lines intercept the profile are the locations of the resonances .
the h resonance crossing probability @xmath175 as a function of neutrino energy for the 1d sn simulation with @xmath43 . in the top panel
the snapshot times are @xmath55 ( solid ) and @xmath41 ( dashed ) , in the middle panel @xmath56 ( solid ) and @xmath180 ( dashed ) and in the bottom panel @xmath181 ( solid ) and @xmath65 ( dashed).,width=480 ]    and in fig .
( [ fig : pc.snv1d.1e35.t=2,2.5,3,3.5,4,4.5 ] ) we show @xmath175 as a function of neutrino energy at six snapshots of this same simulation . the crossing probability as a function of the energy clearly evolves with time .
if @xmath117 is not too small then the evolution of the neutrinos through the undisturbed progenitor profile is almost adiabatic i.e. @xmath182 . as the forward shock arrives at the h - resonance for @xmath183 neutrinos - fig .
( [ fig : rho.snv1d.1e35.t=0.9,1.8,3.6,7.2.+rhores ] ) indicates this occurs at shortly after @xmath184 - the evolution becomes non - adiabatic i.e. @xmath185 because the density jump across the shock straddles the resonance densities of this neutrino energy . as time progresses and the shock moves outwards to lower densities we see from fig .
( [ fig : rho.snv1d.1e35.t=0.9,1.8,3.6,7.2.+rhores ] ) that the shock will begin to affect @xmath186 neutrinos at @xmath46 . the non - adiabaticity sweeps up through the neutrino spectrum from low energy to high . as time progress further eventually the shock ceases to affect the neutrinos of a particular energy and so their propagation returns to being adiabatic . for this simulation fig
( [ fig : rho.snv1d.1e35.t=0.9,1.8,3.6,7.2 ] ) indicates that the evolution of @xmath183 neutrinos returns to being adiabatic at around @xmath46 and that for @xmath186 the return occurs roughly at @xmath47 . note that the @xmath183 neutrinos were affected for a much briefer period ( @xmath187 ) than the @xmath186 neutrinos ( @xmath188 )
this is due to the shape of the density profile .
the h resonance crossing probability @xmath175 as a function of time for selected neutrino energies through the 1d sn model with @xmath43 .
the curves are : @xmath189 ( solid ) , @xmath190 ( long dashed ) , @xmath191 ( dash - dot ) , @xmath192 ( short dashed ) , and @xmath193 ( dash double - dot).,width=480 ]    a complimentary perspective is to look at how particular neutrino energies evolve with time such as those in fig .
( [ fig : pc.snv1d.1e35.e=10,15,20,25,30 ] ) . the temporary transition to non - adiabaticity for each neutrino
energy is clearly visible in the figure .
this figure makes it most obvious that the lower energies are affected before the higher and also that the duration of the non - adiabatic period increases with the neutrino energy .    the extent of the shock feature in neutrino energy seen in fig .
( [ fig : pc.snv1d.1e35.t=2,2.5,3,3.5,4,4.5 ] ) is related to the density jump , @xmath194 , across the shock . at any given time
, if @xmath195 is the highest neutrino energy affected by the shock and @xmath196 is the range of neutrino energies then the density jump across the shock is @xmath197 which is independent of the mixing parameters . from equation ( [ eq : delta rho ] ) we saw that , for normal shocks
, @xmath194 is a function of the mach number and ratio of specific heats @xmath50 .
if @xmath49 is large then @xmath194 , and consequently @xmath196 , are essentially just a function of @xmath50 but if the mach number @xmath49 is not too large - which may be the case - then the width of this feature could be used to infer @xmath49 if @xmath50 is known .     the shock position as a function of time , @xmath198 , for the 1d sn simulation where @xmath43 .
an artificial time delay of @xmath199 has been added to mimic the stalling of the shock at @xmath200.,width=480 ]    as the sn proceeds to explode @xmath195 moves up through the neutrino spectrum .
the correspondence between energy and resonance density means that if we know , or assume , a progenitor profile shape then @xmath195 measures the shock position @xmath201 . from the locations of @xmath195 shown in fig .
( [ fig : pc.snv1d.1e35.t=2,2.5,3,3.5,4,4.5 ] ) , and the initial profile , we show in fig .
( [ fig : rs.snv1d.1e35 ] ) @xmath201 as function of time for this particular simulation .
the change of @xmath201 with time of course measures the shock speed but we can also use @xmath202 to extrapolate the shock position both forward in time , if necessary , but , more interestingly , also backwards toward the proto - neutron star . if this extrapolation can be accomplished successfully then we can infer the moment , @xmath203 , when @xmath200 . the neutronization burst (
if detected ) supplies a zero of time so @xmath204 is the brief period the shock was stalled .
thus the neutrino signal , and in particular the location of the forward shock , contains evidence of whether a key component of the core - collapse supernova paradigm is correct and can provide quantitative data with which to compare with more sophisticated sn simulations than ours .
we illustrate this idea in fig .
( [ fig : rs.snv1d.1e35 ] ) where we have artificially included an offset in time .
the reverse shock was a feature seen in the more powerful explosions from section  [ section : simulations ] and the 2d result .
this feature was not present in the initial profile and appeared later on when when the velocity of the wind superseded the local sound speed .
after its formation the reverse shock initially moved outwards into the star but , as the energy deposition faded with time and the strength of the wind abated , eventually the reverse shock stalled and then headed back to the core . like
the forward shock , the reverse shock affects the adiabaticity of the evolution of the neutrino wavefunction through the profile so that neutrinos with resonance densities straddled by the density jump across the reverse shock evolve non - adiabatically as they pass through it .
note that any neutrino with a resonance density straddled by the density jump across the reverse shock will also experience two other resonances : one before the reverse shock and one after .
the density as a function of the radius in a 1d sn model with @xmath57 at @xmath58 ( dashed ) , @xmath54 ( solid ) and @xmath59 ( long dashed ) .
the horizontal dashed lines are ( from top to bottom ) the resonance densities for @xmath176 , @xmath177 and @xmath179 neutrinos.,width=480 ]    this is shown in fig .
( [ fig : rho.snv1d.2e36.t=1,1.5,2.+rhores ] ) where we have superimposed the resonance densities for @xmath176 , @xmath177 , and @xmath179 neutrinos upon the profiles shown in fig .
( [ fig : rho.snv1d.2e36.t=1,1.5,2 ] ) .
we see from the figure that it is possible for some neutrinos to be affected by both shocks .
if this occurs , and if we neglect the effects from other resonances in the profile , then we might expect that the net effect upon such doubly shock - affected neutrinos to be zero @xcite .
the h resonance crossing probability @xmath175 as a function of neutrino energy for three snapshots taken from the 1d sn model where @xmath52 . in the top panel
the time is @xmath45 , in the middle @xmath59 , and in the bottom panel @xmath60.,width=480 ]    however this expected cancellation is not seen in our results for @xmath175 as a function of the neutrino energy shown in fig .
( [ fig : pc.snv1d.8e35.t=1.8,2,2.4 ] ) for the model where @xmath52 .
the crossing probabilities plotted in the figure show some similarities to those plotted in fig .
( [ fig : pc.snv1d.1e35.t=2,2.5,3,3.5,4,4.5 ] ) at the higher neutrino energies affected by the forward shock where there is change from adiabatic to non - adiabatic evolution .
but for those energies where cancellation is naively expected we see instead that @xmath175 oscillates wildly
. these rapid oscillations are phase effects due to the interference between the two shocks .
similar rapid oscillations in the crossing probability have been seen previously in sn neutrino calculations by fogli _
@xcite in a profile with a forward shock and then a bubble cavity behind it , in the sn test case used by kneller & mclaughlin @xcite , and in the results of dasgupta & dighe @xcite where their significance was emphasized .
the h resonance crossing probability @xmath175 as a function of neutrino energy for three snapshots taken from the 1d sn model where @xmath57 . in the top panel
the time is @xmath205 , in the middle @xmath206 , and in the bottom panel @xmath207.,width=480 ]    fig .
( [ fig : pc.snv1d.2e36.t=1.1,1.4,3.0 ] ) we display the crossing probability for the model where @xmath57 at @xmath205 , @xmath206 and @xmath207 . again phase effects are seen .
compared to the results shown in fig .
( [ fig : pc.snv1d.8e35.t=1.8,2,2.4 ] ) for the weaker explosion where @xmath52 and in fig .
( [ fig : rs.snv1d.1e35 ] ) for @xmath43 , in this case there is no indication of the characteristic transition from adiabatic to non - adiabatic propagation associated with the forward shock .
this occurs because , as noted earlier about this model , the reverse shock penetrates to lower densities than the forward shock at these early times . for this particular model
the reverse shock affects the adiabaticity of a particular neutrino energy before the forward shock .
the presence of phase effects in @xmath175 are caused by interference between forward and reverse shocks .
but it would be nice to find in the neutrino signal a clean signature of the reverse shock that is not contaminated by the forward shock or other features of the profile .
a potential signature could arise from the behavior noted earlier which is that the reverse shock stalled as the energy deposition into the material above the proto - neutron star faded and then headed back towards the core .
as the reverse shock moves backwards its effects move down through the neutrino spectrum and might eventually cease to overlap with the forward shock .
since the density jump across the reverse shock becomes quite small by the time this behavior occurs the reverse shock will produce a narrow spectral feature .
for the simulation where @xmath52 , shown in fig .
( [ fig : rho.snv1d.8e35.t=1,1.5,2,2.5,3 ] ) , this expected pattern is difficult to see in the neutrino signal because the forward shock has not swept through the neutrino spectrum by the time the reverse shock turns around .
but in the simulation with slightly larger energy deposition , @xmath57 , the forward shock moves much more quickly and has largely swept through the h resonance region by the time the reverse shock starts to make its way back to the core . for this simulation
the backwards moving reverse shock is more visible in the neutrino signal . the crossing probability for this simulation at the snapshots shown in fig .
( [ fig : rho.snv1d.2e36.t=4,4.5,4.9 ] ) are presented in fig .
( [ fig : pc.snv1d.2e36.t=4,4.5,4.9 ] ) .
the backwards moving reverse shock is responsible for the phase effects at @xmath208 in the middle panel . and
then once the reverse shock has reached the core our boundary conditions there led to its reflection and so the reverse shock became a outward moving , weak forward shock and re - appeared in the neutrino spectrum as a spectral feature moving from low to high neutrino energies .
this is the cause of the bump in @xmath175 seen at @xmath209 in the bottom panel .
the h resonance crossing probability @xmath175 as a function of neutrino energy for the 1d sn model where @xmath57 . from top
to bottom the snapshot times are @xmath181 , @xmath65 and @xmath66.,width=480 ]    in summary , though it is apparent that the neutrino signal can vary considerably depending upon exactly how the reverse shock behaves the presence of strong phase effects - rapid oscillations with large amplitude in both time and energy of the crossing probability - are a notable signature of the presence of multiple shocks .
therefore , there exists tremendous potential in a future detection of the a supernova neutrino signal for probing the inner hydrodynamics of the explosion .
the notable difference between the one - dimensional and the two - dimensional models was the local density fluctuations that appeared in the latter .
these local density enhancements / cavities affect only neutrinos along particular lines of sight so if one could determine their presence in the signal we would obtain an indication of the asphericity of the explosion . the most obvious method would be to compare different lines of sight to the same sn but obviously this is not practical . instead we must find another signature of asphericity that one might use .
the h resonance crossing probability @xmath175 as a function of neutrino energy for a radial slice at @xmath74 through the 2d sn model . from top
to bottom the snapshot times are @xmath60 , @xmath76 and @xmath210.,width=480 ]    in figure ( [ fig : pc.snv2d.t=2.4,5.4,6.4 ] ) we show the crossing probability @xmath175 as a function of neutrino energy for the snapshot @xmath60 , @xmath76 and @xmath210 and a line of sight @xmath74 .
features in the figure resemble those seen in the 1d results : the forward shock leads to a range of energies with @xmath211 as seen in the top panel , phase effects due to presence of the reverse shock lead to the high ` frequency ' oscillations of @xmath175 seen in the middle panel , and the turn around of the reverse shock can be seen in the bottom panel at @xmath212 where it appears , as before , as a narrow range of high frequency phase effects moving down through the spectrum . while a detailed analysis of these results might indicate a difference that is due to asphericity there is no striking feature that one can point to .
the h resonance crossing probability @xmath175 as a function of neutrino energy for a radial slice at @xmath74 through the 2d sn model . from top
to bottom the snapshot times are @xmath213 , @xmath214 and @xmath215.,width=480 ]    however , at late times we found a change in the qualitative appearance of @xmath175 versus neutrino energy .
the high frequency phase effects no longer dominate the entire spectrum and now lower frequency portions show up .
we did not see anything comparable at late times in any of our 1d simulation results .
the changing quality of the phase effects are seen in fig .
( [ fig : pc.snv2d.t=7.4,8,9 ] ) which shows the crossing probability at the h resonance as a function of neutrino energy for the same line of sight at @xmath74 through the 2d model .
high frequencies can be seen above @xmath216 in the top panel and above @xmath217 in the middle panel but they have disappeared by @xmath215 . in all three panels
we see that at lower neutrino energies the curve is smoother with only narrow patches of high frequency phase effects that moved up through the spectrum .
the shocks have largely ceased affecting the lower neutrino energies by this time so it is the ` fluctuations ' upon the profile - both between the shocks and behind the reverse shock - which produce mildly adiabatic resonances , and/or resonances that are located very close to one another ( and possibly even overlap ) that give rise to these features in the figure .
the effect upon neutrino flavor transport of supernovae density fluctuations was considered by loreti et al @xcite and later by fogli _
et al _ @xcite and friedland & gruzinov @xcite .
these authors suggest that @xmath175 should tend to 1/2 in the presence of many small scale fluctuations .
we note , however , that due to the exponential growth of the radial grid spacing the scale of features we can resolve by @xmath218 is of order @xmath219 and this is somewhat larger than the scale considered by fogli _
et al _ and friedland & gruzinov . in our formulation of neutrino
mixing using s - matrices we can motivate this expected result of @xmath220 by dividing the density profile into domains such that within each there is a single neutrino resonance .
the passage of a neutrino ( or antineutrino ) through each interval is described by an s - matrix which has exactly the same structure as equations ( [ eq : shnu ] ) or ( [ eq : shnubar ] ) as appropriate .
the magnitudes of the elements @xmath154 and @xmath221 lie between zero and unity depending upon the adiabaticity of the resonance in that region .
the multiplication of all these s matrices together ( in the correct order ) to achieve the net effect of passing through all the resonances will tend to produce a net result where either @xmath222 and @xmath223 or @xmath224 and @xmath225 are of equal magnitude . in our simulations
we had cases where we found up to 13 resonances but have not been able to reproduce the anticipated limiting behavior . nevertheless , as shown in figs .
( [ fig : pc.snv2d.t=2.4,5.4,6.4 ] ) and ( [ fig : pc.snv2d.t=7.4,8,9 ] ) , we see a clear distinctive qualitative change in the crossing probabilities at late time in aspherical models .
since the detection of neutrinos from sn 1987a in imb @xcite and kamiokande @xcite the number and size of the neutrino detectors have both increased in scale so that many thousands or even tens of thousands of neutrinos will be detected from the next galactic supernova .
though many neutrino detectors are water cerenkov detectors that are best suited to detecting the electron antineutrinos from supernovae via the inverse-@xmath221 reaction on protons , other neutrino detection technologies have been demonstrated or explored that allow one to see the electron neutrinos and the @xmath226/@xmath32 flavors via nucleus - neutrino reactions such as @xmath227 , or neutral current reactions such as @xmath228 .
the detection of these other neutrino flavors would provide a great deal of valuable information we could then use to learn about the sn explosion .
for the purposes of demonstrating the temporal evolution of supernovae neutrino signals we shall consider two detector types : the water cerenkov detector because the majority of neutrino detectors are of this type , and a heavy water detector ( such as was used in the sudbury neutrino observatory ) since this detector was also considered by schirato & fuller @xcite .
the largest neutrino detectors are water cerenkov detectors .
these detectors are sensitive to both electron neutrinos and electron anti neutrinos via the reactions @xmath229 and @xmath227 @xcite with the inverse-@xmath221 reaction dominating the signal .
the cross sections for the neutrino - nucleon reactions are well known but a good deal of uncertainty exists in the important neutrino - nucleus cross sections .
a heavy water detector is sensitive to both charged - current and neutral current events so could be used to detect the @xmath226 and @xmath32 flavors .
however obtaining energy resolution for the neutral current events is difficult so the observable quantity considered by schirato & fuller @xcite was the ratio of total charged - current to neutral current events .
though temporal evolution of this ratio is somewhat smothered due to the fact that the event rates are integrated over the neutrino spectrum schirato & fuller demonstrated that variations of the ratio would occur at the level of @xmath230 .
lead detectors are also capable of detecting the @xmath226/@xmath32 flux and have been studied by fuller , haxton & mclaughlin @xcite engel , mclaughlin & volpe @xcite and kolbe et al @xcite .    to determine what any given detector will see we must fold in three effects : the initial neutrino spectra emitted by the proto - neutron star , the energy dependence of the cross sections in the detector and , thirdly , the energy resolution of the detector .
there is a considerable degree of uncertainty about the emitted neutrino spectra .
furthermore , neutrino self interactions - described in samuel @xcite , pastor , raffelt & semikoz @xcite and friedland & lunardini @xcite - in the region immediately above the proto - neutron star can alter the spectra after it s emission from the proto - neutron star but before the flux reaches the h resonance layer .
there is a possibility that in some cases the spectra may be ` equilibrated ' by the time they get to the h resonance @xcite . in such circumstances the changes in the density profile and the neutrino mixing described in previous sections
would not lead to any observable changes in a detector except for the gradual fading of the signal .
we shall consider the case where such equilibration does not occur and experiment with an ` alpha fit ' power law for the spectrum exiting the neutrinosphere motivated both by its analytic simplicity and by recent supernova core monte carlo simulations @xcite . for the time dependence of the luminosity we use an exponential decay with a time constant of @xmath231 .
the emitted differential spectra are then @xmath232\,\exp[-t/\tau ] .
\label{eq : powlaw}\ ] ] where @xmath233 is the luminosity , @xmath234 the mean energy and @xmath154 is the ` pinch ' parameter .
we adopt a hierarchy of mean energies : @xmath235 , @xmath236 and @xmath237 ( where , again , @xmath238 stands for either @xmath239 or @xmath240 ) and the luminosities are @xmath241 , @xmath242 .
fluxes of electron antineutrinos ( top panel ) and the sum of @xmath243 neutrino fluxes ( bottom panel ) for completely adiabatic , i.e. @xmath244 ( solid ) and non - adiabatic , i.e. @xmath245 ( dashed ) evolution .
the hierarchy is inverted , the pinch parameter was set to @xmath246 and the vacuum mixing angles are @xmath247 , @xmath248 .
the distance to the sn was set at @xmath249 and the time @xmath31 to @xmath250.,width=480 ]     the fluxes of electron neutrinos ( top panel ) and the sum of @xmath243 neutrino fluxes ( bottom panel ) for completely adiabatic , i.e. @xmath244 ( solid line ) and non - adiabatic , i.e. @xmath245 ( dashed ) evolution .
here the hierarchy is normal but otherwise all parameters are the same as those used in fig .
( [ fig : fnu.alpha=3.ih.phbar=0,1]).,width=480 ]    in figs .
( [ fig : fnu.alpha=3.ih.phbar=0,1 ] ) and ( [ fig : fnu.alpha=3.nh.ph=0,1 ] ) we present the neutrino spectra that would be observed at earth from a sn at a distance of @xmath251 if the evolution were either completely adiabatic or completely non - adiabatic for example through a single h resonance , i.e. @xmath252 or @xmath245 in either the inverted or normal hierarchy . the two curves in each panel of the figures therefore represent extremes in the sense that the actual flux at a given energy would fall somewhere in between the two curves .
we do not plot the electron antineutrino flux for the normal hierarchy or the electron neutrino flux for the inverted hierarchy because these fluxes are always adiabatic and would only fade exponentially with time .
in each panel of the two figures the pinch parameter was set to @xmath246 , other than a change in shape the general behavior of the curves at other values of @xmath154 is unchanged .
there are many points to notice in each figure .
the first is that for any given neutrino flux and hierarchy there is a particular energy where the two curves cross .
this will occur as long as the average energies of @xmath253 flavor neutrinos are greater than that of the @xmath0 and @xmath254 .
the crossing can be seen best in the @xmath255 spectrum in the normal hierarchy .
we define the energies where the curves cross as the critical energies .
the critical energies vary with both flavor and the hierachy . at the critical energy for a particular hiearachy and flavor the adiabaticity of the neutrinos
as they propagate through the profile is irrelevant and the flux of that flavor and at that energy through a detector is independent of @xmath175 .
the other significant feature to notice is the large difference between the two curves for the electron neutrinos and electron antineutrinos for energies below @xmath256 .
this wide separation of the extremal fluxes at low to medium neutrino energies , @xmath257 , suggests that we focus our attention upon this portion of the spectrum because it is here that changes in adiabaticity will lead to the most significant changes in the flux .
the flux at higher energies will also change with the evolving adiabaticity but to a lesser degree because the fluxes are so small .
thus if the evolution of the neutrinos evolves from being adiabatic to non - adiabatic then , in the case of a normal hierarchy , the flux of electron neutrinos will increase for those energies below the critical energy and decrease above .
similar behavior will occur for antineutrinos in the inverted hierarchy .
the decrement above the critical energy is difficult to see in the figures because the critical energies are all far into the tail of the spectra . for the @xmath258 neutrino
fluxes the opposite change occurs : a transition to non - adiabaticity leads to lower fluxes of these flavors at energies below the critical energy and higher fluxes above .
having determined the effects of the evolving density profile upon the neutrino mixing we are now in a position to examine the actual neutrino signal within a detector .
as discussed in section [ sec : pfcp ] , the actual neutrino flux of a particular flavor through a detector differs from that emitted from the neutrinosphere . in the simplest scenario ,
if the neutrinos propagate through the outer layers of the sn adiabatically then the survival probabilities @xmath259 and @xmath260 are found by setting @xmath261 in equations ( [ eq : ppbarnh ] ) and ( [ eq : ppbarih ] ) . for the normal hierarchy @xmath262 and @xmath263 : for the inverse hierarchy @xmath264 and @xmath265 .
using the value of @xmath266 we used to calculate @xmath175 and @xmath267 based upon the solar mixing parameters , we see that for the normal hierarchy the amount of the emitted @xmath268 flux present in the electron neutrino flux at a detector here on earth is minuscule . similarly we see that for an inverted hierarchy the amount of the emitted electron antineutrino flux present in any detected electron antineutrino flux is equally small .
thus if the mass hierarchy is normal , @xmath269 is not too small , and the neutrino propagation adiabatic then the electron neutrino detectors on earth will essentially be detecting the emitted @xmath270 and @xmath271 flux . or
if the hierarchy is inverted , again @xmath269 is not too small , and the neutrino propagation adiabatic then the electron antineutrino flux at earth is the emitted @xmath272 and @xmath273 . in order to see the emitted electron neutrino or antineutrino flux we would need to detect the @xmath270 , @xmath272 , @xmath273 and/or @xmath271 flavors or seek the small amount of the emitted @xmath268/@xmath254 present in the detected @xmath0/@xmath255 flux . in either case
, it will be very difficult to reconstruct the the @xmath268 or @xmath254 spectrum that was originally emitted from the proto - neutron star .
the evolution of the density profile changes this conclusion . depending upon the scenario ,
either @xmath259 ( for a normal hierarchy ) or @xmath260 ( for an inverted hierarchy ) will change with time and so the proportions of the emitted fluxes present in the fluxes at earth will also change .
note that @xmath259 and @xmath260 do not both evolve with time within in the same scenario ; one will stay constant .
this means there is a potential to separate the time evolution of the emitted flux from the neutrinosphere from the evolution of the density profile through which the neutrinos propagate .
as the star explodes and the forward shock eventually reaches the h resonances the passage of portions of the neutrino spectrum through the density profile becomes non - adiabatic .
whether it is the neutrinos or antineutrinos that are affected is determined by the hierarchy .
when the passage of the neutrinos is completely non - adiabatic we see that @xmath259 or @xmath260 change to @xmath264 or @xmath263 . again using @xmath274 and @xmath267 we see that @xmath275 , @xmath276 .
in effect the forward shock of such explosions acts like a window sweeping across either the neutrino or antineutrino spectra through which we observe an altered mixture of the emitted components . for neutrino or antineutrino energies below the critical energies
this leads to an enhanced @xmath0 or @xmath255 flux as shown in figs .
( [ fig : fnu.alpha=3.ih.phbar=0,1 ] ) and ( [ fig : fnu.alpha=3.nh.ph=0,1 ] ) .
the ratio of charged current event rates to neutral current event rates in a heavy water detector as a function of time for neutrino propagation through the @xmath67 radial slice of the the 2d simulation .
the top panel is for an inverted hierarchy , the lower panel for a normal hierarchy . in each
the lines represent different spectral parameters : @xmath277 ( solid ) , @xmath278 ( dotted ) , @xmath246 ( long dashed ) , @xmath279 ( dash dot ) and @xmath280 ( long dash - double dot).,width=480 ]     the ratio of charged current event rates to neutral current event rates in a heavy water detector as a function of time for neutrino propagation through the 1d simulation with @xmath43 .
the top panel is for an inverted hierarchy , the lower panel for a normal hierarchy . in each
the lines are : @xmath277 ( solid ) , @xmath278 ( dotted ) , @xmath246 ( long dashed ) , @xmath279 ( dash dot ) and @xmath280 ( long dash - double dot).,width=480 ]     the ratio of charged current event rates to neutral current event rates in a heavy water detector as a function of time for neutrino propagation through the 1d simulation with @xmath52 .
the top panel is for an inverted hierarchy , the lower panel for a normal hierarchy . in each
the lines are : @xmath277 ( solid ) , @xmath278 ( dotted ) , @xmath246 ( long dashed ) , @xmath279 ( dash dot ) and @xmath280 ( long dash - double dot).,width=480 ]     the ratio of charged current event rates to neutral current events rates in a heavy water detector as a function of time for neutrino propagation through the 1d simulation with @xmath57 .
the top panel is for an inverted hierarchy , the lower panel for a normal hierarchy . in each
the lines are : @xmath277 ( solid ) , @xmath278 ( dotted ) , @xmath246 ( long dashed ) , @xmath279 ( dash dot ) and @xmath280 ( long dash - double dot).,width=480 ]    the first detector type we consider is a heavy water , sno like , detector and the signal we focus upon is the ratio of charged current ( cc ) - both from electron neutrinos and electron antineutrinos - to neutral current ( nc ) event rates .
this ratio of event rates is the same signal considered by schirato & fuller @xcite even though the @xmath281 counters and then the nacl inserted into sno gave the detector some capability of distinguishing between the @xmath0 and @xmath255 events . to some extent the uncertainty of the initial spectra can be removed by considering this ratio of event rates though by no means is it eliminated . for the results that follow we calculate the event rates using the nakamura _
@xcite neutrino - deuteron cross sections .
the basic understanding of what happens to the @xmath0 or @xmath255 flux helps us understand the temporal evolution of the cc / nc ratio shown in fig .
( [ fig : ccnc.snv2d ] ) for the 2d simulation and for the 1d simulations in figs .
( [ fig : ccnc.snv1d.1e35 ] ) , ( [ fig : ccnc.snv1d.8e35 ] ) and ( [ fig : ccnc.snv1d.2e36 ] ) . for all combinations of hierarchy ,
pinch parameter and explosion energy the ratio begins to vary at around @xmath282 into the simulation .
the initial increase in the ratio is due to larger charged current rates ( the neural current rate is flavor blind ) caused by greater electron neutrino flavor content ( for the nh ) or electron antineutrino flavor ( for the ih ) .
the change is a clear signature of a transition from adiabaticity to non - adiabaticity that occurs when the forward shock reaches the h resonances for neutrinos with @xmath283 .
as expected , the ratio varies more strongly for an inverted hierarchy due to the greater change in the electron antineutrino spectrum when the profile becomes non - adiabatic than for the electron neutrinos as indicated by figs .
( [ fig : fnu.alpha=3.ih.phbar=0,1 ] ) and ( [ fig : fnu.alpha=3.nh.ph=0,1 ] ) .
more detailed examination of the figures reveals that for the particular case of the weak explosion - fig .
( [ fig : ccnc.snv1d.1e35 ] ) - there is evidence of the shock effects moving up through the neutrino spectra because at late times we see the ratio dips beneath it s @xmath284 value .    for the other simulations the cc / nc ratio exhibits a great deal of variability at the level of a few percent that continues for several seconds .
we also notice that the variability of the @xmath285 and @xmath286 1d simulations ceased after @xmath287 when the forward shock had swept through the spectrum and the reverse shock had turned around while the cc / nc ratio for the 2d simulation was sustained for the entire @xmath288 . while observing a period of fluctuations in the cc / nc ratio
would indicate that the explosion and the density profile was dynamic it is difficult to be more quantitative about the presence and motion of the various features of the profile from the use of this ratio of event rates alone .
the second example of a detector we examine is a water cerenkov type since this is the dominant supernova neutrino detection technology .
we focus upon the positron spectrum created by inverse @xmath221 reactions of antineutrinos upon protons .
the temporal variations of the neutrino signal we have discussed in the paper would therefore require an inverted hierarchy in order to be visible in this signal .
the positron spectrum seen in a detector , @xmath289 , is given by @xmath290 where @xmath291 is the number of protons in the detector , @xmath292 is the electron antineutrino flux and @xmath293 is the differential cross section . for our analysis we used the cross sections of strumia and vissanti @xcite .
the positron spectrum allows us to see the effect of the density profile upon the neutrino flux as a function of neutrino energy but , before we show our results , we firstly consider the relationship between the neutrino energy and the positron energy and our ability to resolve the features of the neutrino spectrum .
the ` mean ' positron energy @xmath294 as a function of the neutrino energy ( solid line ) given by equation ( 16 ) in strumia and vissanti @xcite .
the dashed line is @xmath295.,width=480 ]    the positron energy and neutrino energy are closely correlated as shown in fig .
( [ fig : meaneae.eps ] ) .
features in the neutrino spectrum at a particular energy will appear in the positron spectrum at more or less the same energy .
the super - k energy resolution error , @xmath296 , ( solid line ) and the range of positron energies , @xmath297 , ( dashed ) as a function of the positron kinetic energy.,width=480 ]    but the positron spectrum is smoothed relative to the neutrino spectrum by two effects : the distribution of positron energies for a given neutrino energy , and the energy resolution of the detector . for a given neutrino energy
the positrons emerge with a range of energies which we may invert so that for a given positron energy , @xmath298 , there is a range of neutrino energies , @xmath297 , from which it could have arisen .
this range . ] is given by @xmath299 where @xmath300 is the momentum of the positron and @xmath301 , @xmath302 and @xmath303 are the masses of the neutron , proton and electron respectively .
@xmath297 is shown in fig .
( [ fig : skerror , spread ] ) as a function of the positron kinetic energy .
and then for the energy resolution of the detector we used a filter so that the spectrum of positrons at energy @xmath298 is given by @xmath304 and we take @xmath305 to be a normalized gaussian with a mean @xmath306 and a variance @xmath307 given by @xmath308}\right)^{0.3467}.\ ] ] this formula is our fit to the energy resolution of the super - k detector @xcite .
the energy resolution of the detector , @xmath296 , is also shown in figure ( [ fig : skerror , spread ] ) as a function of the positron kinetic energy .
both smearing effects , the range , @xmath297 , of neutrino energies that can produce a positron with energy @xmath298 , and the detector resolution , @xmath296 , increase with @xmath298 and we notice that beyond @xmath309 it is the range @xmath297 that determines that determines the scale of the features we can resolve . below @xmath309 the detector resolution limits out ability to see features of the neutrino spectrum .
future detectors will improve upon the energy resolution used here and so - since this is also the energy range where we should focus our attention because it is where the electron antineutrino flux will vary the largest due to the changing adiabaticity - in all our figures we plot both the spectrum as it would be measured with the super - k energy resolution and with a detector that has perfect resolution .
of course there is a third effect one needs to consider : the poisson error due to the finite number of detected events .
this does not smooth out features of the spectrum but rather make them difficult to detect . in what follows
we shall ignore the poisson statistics .    .
in the top panel the times are @xmath250 ( solid ) and @xmath41 ( long dashed ) , in the middle panel @xmath310 ( solid ) and @xmath64 ( long dashed ) , and in the bottom panel @xmath311 ( solid ) and @xmath312 ( long dashed ) . in all panels
the short dashed lines indicate the spectrum convolved with our adopted detector energy resolution .
[ fig : faevse.1e35],width=480 ]    in fig .
( [ fig : faevse.1e35 ] ) we plot the positron spectrum as a function of positron energy in the case of the inverted hierarchy at various snapshot times using the neutrino flux we calculated for the weakest 1d simulation , @xmath313 . with this spectral information
, we can begin to see the features of the explosion .
the shock is clearly seen at @xmath314 at @xmath41 in the perfectly resolved spectrum .
by @xmath310 the sharp change in spectrum has moved up to @xmath315 and by @xmath64 we begin to see the backside of the shock at @xmath314 and again at @xmath311 at @xmath316 and @xmath312 at @xmath317 .
the shock is a less dramatic feature in the spectra convolved with our adopted detector energy resolution but is visible particularly at @xmath41 . given enough events even present detector technology would allow us to see the shock move up through the positron spectrum from which we could begin to glean information about the explosion .    .
the curves are @xmath189 ( long dashed ) , @xmath191 ( dash dot ) , and @xmath193 ( solid ) .
the short dashed lines indicate the spectrum convolved with our adopted detector energy resolution .
[ fig : faevst.1e35],width=480 ]    the explosion features also appear in the evolution of the positron spectrum at particular energies as a function of time .
this is shown in figure ( [ fig : faevst.1e35 ] ) where we plot the positron spectrum for @xmath189 , @xmath191 and @xmath193 .
the general trend for each energy is @xmath318 since this was the time dependence we used for the emitted spectra but on top of this exponential decay the shock is obviously visible for all three chosen energies in both filtered and unfiltered curves .
the width of the forward shock feature grows with @xmath94 as we noted earlier .    .
in the top panel the times are @xmath54 ( solid ) and @xmath319 ( long dashed ) , in the middle panel @xmath60 ( solid ) and @xmath310 ( long dashed ) , and in the bottom panel @xmath46 ( solid ) and @xmath66 ( long dashed ) . in all panels
the short dashed lines indicate the spectrum convolved with our adopted detector energy resolution.[fig : faevse.8e35],width=480 ]    in fig .
( [ fig : faevse.8e35 ] ) we plot the positron spectrum for the more powerful explosion with @xmath285 at various snapshot times .
the phase effects we saw in the crossing probability for this model appear also in the spectrum for perfect detector resolution .
they are largest in the range of @xmath320 where the difference between adiabatic and non - adiabatic electron antineutrino flux is largest and where the distribution of positron energies for a given neutrino energy is narrow . at higher energies
the phase effects largely disappear because of the smoothing effect of the positron distribution .    , and an inverted hierarchy .
the curves are @xmath191 ( dash dot ) , and @xmath193 ( solid ) .
the short dashed lines indicate the spectrum convolved with our adopted detector energy resolution.[fig : faevst.8e35],width=480 ]    this smoothing is also seen in the plot of the spectrum at various energies as a function of time shown in fig .
( [ fig : faevst.8e35 ] ) since the @xmath321 positrons has much greater variability than then @xmath322 positrons even in the filtered spectrum .    .
in the top panel the time is @xmath64 , in the middle panel @xmath65 , and in the bottom panel @xmath66 . in all panels the short dashed lines indicate the spectrum convolved with our adopted detector energy resolution.[fig : faevse.2e36],width=480 ]    we can even find the signature of a reverse shock that stalled and headed back to the core . in fig .
( [ fig : faevse.2e36 ] ) we plot the positron spectrum for the more energetic @xmath286 1d explosion .
the snapshot times are the same as those shown in fig .
( [ fig : rho.snv1d.2e36.t=4,4.5,4.9 ] ) for the density and fig .
( [ fig : pc.snv1d.2e36.t=4,4.5,4.9 ] ) for the crossing probability for this model .
the returning reverse shock appears at @xmath323 at @xmath65 which matches what we expect from the crossing probability seen in the middle panel of fig .
( [ fig : pc.snv1d.2e36.t=4,4.5,4.9 ] ) .
( solid ) and @xmath60 ( long dashed ) , in the middle panel @xmath210 ( solid ) and @xmath213 ( long dashed ) , and in the bottom panel @xmath214 ( solid ) and @xmath215 ( long dashed ) . in all panels
the short dashed lines indicate the spectrum convolved with our adopted detector energy resolution.[fig : faevse.snv2d],width=480 ]    finally in fig .
( [ fig : faevse.snv2d ] ) we plot the spectrum for the 2d simulation at various snapshot times selected so that they match as closely as possible the snapshots we have shown of the profile and/or the crossing probability .
the forward shock and phase effects can be seen here too in the unfiltered spectrum and the inward moving reverse shock is responsible for the phase effects seen at @xmath316 at @xmath210 and in fig .
( [ fig : pc.snv2d.t=2.4,5.4,6.4 ] ) .
other features upon the positron spectrum can also be related to various features shown in the crossing probability for this model .
the next galactic supernova has the potential to provide a great deal of information about both neutrino mixing and about the supernova .
the temporal and spectral evolution of the signal is altered as the star explodes due to the changes in the density profile of the star .
the forward shock should be easily detected even with present detector technology as either a feature sweeping up through the positron spectrum in a water cerenkov detector or as a sudden increase in event rates for a detector with no spectral information such as a heavy water detector .
once identified , and if the progenitor were known , the information can be used to measure the shock speed , the strength of the shock and potentially also the stalled shock time delay .
the reverse shock can be spotted qualitatively as the cause of fluctuations in either event rates or upon the positron spectrum .
there is also the possibility for the reverse shock to be seen if its motion is reversed and it heads back to the core .
the neutrino signal may also contain hints that the supernova was aspherical as a qualitative change in the phase effects .
the authors are very grateful to john blondin for his helpful advice with the hydrodynamic simulations .
this work was supported by the u.s .
department of energy at umn under grant de - fg02 - 87er40328 and at ncsu under grant de - fg02 - 02er41216 ..      j.  f.  beacom , r.  n.  boyd and a.  mezzacappa , , * 63 * , 073011 ( 2001 ) [ arxiv : astro - ph/0010398 ] . g.  c.  mclaughlin and r.  surman , , * 75 * , 023005 ( 2007 ) [ arxiv : astro - ph/0605281 ] .
d.  h.  hartmann , r.  d.  hoffman , s.  e.  woosley and w.  c.  haxton , nucl .
a * 527 * , 663c ( 1991 ) .
a.  heger , e.  kolbe , w.  c.  haxton , k.  langanke , g.  martinez - pinedo and s.  e.  woosley , phys .
b * 606 * , 258 ( 2005 ) [ arxiv : astro - ph/0307546 ] .
s.  e.  woosley , j.  r.  wilson , g.  j.  mathews , r.  d.  hoffman and b.  s.  meyer , , * 433 * , 229 ( 1994 ) .
r.  surman , g.  c.  mclaughlin and w.  r.  hix , , * 643 * , 1057 ( 2006 ) [ arxiv : astro - ph/0509365 ] .
j.  pruet , r.  surman and g.  c.  mclaughlin , , * 602 * , l101 ( 2004 ) [ arxiv : astro - ph/0309673 ] .
g.  c.  mclaughlin , g.  m.  fuller and j.  r.  wilson , , * 472 * , 440 ( 1996 ) [ arxiv : astro - ph/9701114 ] .
g.  c.  mclaughlin and g.  m.  fuller , , * 464 * , l143 ( 1996 ) [ arxiv : astro - ph/9603082 ] .
b.  s.  meyer , g.  c.  mclaughlin and g.  m.  fuller , , * 58 * , 3696 ( 1998 ) [ arxiv : astro - ph/9809242 ] .
t.  yoshida , t.  kajino , h.  yokomakura , k.  kimura , a.  takamura and d.  h.  hartmann , arxiv : astro - ph/0606042 . g.  m.  fuller , w.  c.  haxton and g.  c.  mclaughlin , , * 59 * , 085005 ( 1999 ) [ arxiv : astro - ph/9809164 ] .
a.  s.  dighe and a.  y.  smirnov , , * 62 * , 033007 ( 2000 ) [ arxiv : hep - ph/9907423 ]
. j.  engel , g.  c.  mclaughlin and c.  volpe , , * 67 * , 013005 ( 2003 ) [ arxiv : hep - ph/0209267 ] . c.  lunardini and a.  y.  smirnov , jcap
* 0306 * , 009 ( 2003 ) [ arxiv : hep - ph/0302033 ] .
a.  friedland and c.  lunardini , , * 68 * , 013007 ( 2003 ) [ arxiv : hep - ph/0304055 ] .
r.  f.  sawyer , , * 72 * , 045003 ( 2005 ) [ arxiv : hep - ph/0503013 ] .
h.  duan , g.  m.  fuller , j.  carlson and y.  z.  qian , , * 97 * , 241101 ( 2006 ) [ arxiv : astro - ph/0608050 ] .
s.  hannestad , g.  g.  raffelt , g.  sigl and y.  y.  y.  wong , , * 74 * , 105010 ( 2006 ) [ arxiv : astro - ph/0608695 ] .
j.  p.  kneller and g.  c.  mclaughlin , , * 73 * , 056003 ( 2006 ) [ arxiv : hep - ph/0509356 ] .
b.  dasgupta and a.  dighe , [ arxiv : hep - ph/0510219 ] .
blondin , j.  m. , mezzacappa , a. , & demarino , c.  2003 , , * 584 * , 971 blondin , j.m . and mezzacappa , a. , , * 642 * , 401 ( 2006 ) ohnishi , n. , kotake , k. and yamada , s. , , * 641 * , 1018 ( 2006 ) scheck , l. , kifonidis , k. , janka , h .- th . and mueller , e. , astro - ph/0601302    a.  burrows , e.  livne , l.  dessart , c.  d.  ott and j.  murphy , [ arxiv : astro - ph/0610175 ] .
schirato , r.c . and fuller , g.,m . ,
astro - ph/0205390 takahashi , k. , sato , k. , dalhed , h.e . and
wilson , j. r. , astropart .
phys , 20 , 1989 ( 2003 ) fogli , g. l. , lisi , e. , and mirizzi , a. and montanino , d. , , * d68 * , 033005 ( 2003 )        a.b .
balantekin , and j.f .
beacom , , * 54 * , 6323 ( 1996 ) a.n .
ioannisian and a.yu .
smirnov , , * 93 * , 241801 ( 2004 ) a.n .
_ et al . _
_ , , * 71 * , 033006 ( 2005 ) e.kh .
akhmedov , m. maltoni and a.yu .
smirnov , , * 95 * , 211801 ( 2005 )      f.  n.  loreti , y.  z.  qian , g.  m.  fuller and a.  b.  balantekin , , * 52 * , 6664 ( 1995 ) [ arxiv : astro - ph/9508106 ] .
fogli , g. , lisi , e. , mirizzi , a. and montanino , d. , , jcap , * 6 * , 12 ( 2006 ) friedland , a. and gruzinov , a. , arxiv : astro - ph/0607244 ( 2006 )