despite decades of study , type ia supernovae ( snia ) continue to be an active topic for astrophysical research . the accepted physical mechanism for these events , that they are the result of thermonuclear explosions of degenerate material
( hoyle & fowler 1960 ) in white dwarf stars , is very well established but owing to the considerable complexity of this process many aspects of these spectacular events remain poorly understood .
observationally , it is clear that all snia are not the same  there is significant diversity in both brightness and decay timescale ( phillips 1993 ) . aside from its direct relevance to the study of snia themselves , understanding this diversity has important implications to other branches of astrophysics  in particular to cosmology where inferences about the properties of distant snia play an important role in probing the rate of expansion of the universe .
some correlations between supernova properties are already well established from observations of nearby snia ( e.g. phillips 1993 ) .
however , it is becoming clear that the observed diversity is not adequately described by a single parameter ( benetti et al .
2004 , 2005 ) and that unraveling this diversity relies on a combination of careful observational study and sophisticated theoretical modelling .    in recent years
, there has been rapid development in the sophistication of numerical modelling of the hydrodynamics of snia explosions . in particular , while in earlier work supernovae explosions were modelled using first one - dimensional ( 1d ) computer codes ( e.g. nomoto , thielemann & yokoi 1984 ; hflich , wheeler & thielemann 1998 ) and later two - dimensions ( e.g. mller & arnett 1986 ; niemeyer , hillebrandt , woosley 1996 ) , the most up - to - date simulations are fully three - dimensional ( 3d , e.g. reinecke , hillebrandt & niemeyer 2002 ; gamezo et al .
2003 ; rpke 2005 ; rpke et al .
these multi - dimensional models are crucial for the understanding of realistic flame propagation and hence nucleosynthesis in snia explosions and have clearly demonstrated that the earliest 1d models greatly underestimate the likely complexity of real snia .
this naturally raises the question of whether and to what extent this multi - dimensional complexity may be responsible for the observed diversity of explosions .    to connect hydrodynamical explosion models and observations of real snia light curves or spectra requires modelling of radiation transport in the supernova .
the majority of the light escaping from a supernova explosion at around optical maximum originates from energy deposited in the ejecta by the absorption and compton scattering of @xmath0-rays emitted by radioactive isotopes .
given that all snia explosion models are significantly optically thick at the earliest times ( @xmath1  day ) , the emission , transport and deposition of radiation need only be followed for times after most of the complex dynamics have ceased and the ejecta is in near - homologous expansion ; this greatly simplifies the radiation transport calculation and allows it to be decoupled from the hydrodynamical simulation .
the level of sophistication of modern radiation transport computations for supernovae has gradually developed to keep pace with the development of the hydrodynamical explosion models : for earlier applications ( e.g. branch et al .
1981 , lucy 1987 , mazzali & lucy 1993 ) simple one - dimensional calculations were sufficient but there is now growing interest in fully three - dimensional simulations .
in particular , lucy ( 2005 ) has recently described and tested a sophisticated monte carlo approach to 3d , time - dependent radiation transport .
this method holds much promise for a wide range of applications and is already being adopted and developed in several contexts  for example , kasen , thomas & nugent ( 2006 ) have developed a code based on that described by lucy ( 2005 ) which incorporates a treatment of polarisation and lifts the assumption of grey opacity used in the earlier code .
independently , maeda , mazzali & nomoto ( 2006 ) have also adopted lucy s techniques and utilized them in multi - dimensional simulations of light curves for core collapse supernovae .    in this investigation ,
monte carlo calculations using methods similar to those of lucy ( 2005 ) are used to investigate aspects of the influence of the 3d structure of snia on their observable properties . in section  2 ,
the computer code used for the calculations is described .
although similar to the codes described by lucy ( 2005 ) and kasen et al .
( 2006 ) , there are key differences , specifically in the means by which the observable light curves are extracted from the monte carlo simulation .
this code is then used to investigate the effect of departures from spherical symmetry for two physically distinct cases .
first , in section  3 toy models are used to investigate the influence of large scale ( low - mode ) asymmetries on observable light curves .
secondly , a representation of a modern 3d explosion model is used to investigate the implications of the complex , relatively small - scale inhomogeneity predicted ( section  4 ) .
the emphasis here is on identifying and interpreting differential effects between 3d and 1d radiative transfer and thus several simplifying assumptions are made in the treatment of the micro - physics .
conclusions are drawn in section  5 .
the calculations presented in this paper have been performed using a monte carlo code which follows the propagation of radiative energy in 3d as a function of time .
this code is closely based on that described by lucy ( 2005 ) and thus only a brief summary of the operation of the code is given here , with particular emphasis on the the departures from the approach described by lucy ( 2005 ) . the nomenclature used by lucy for the monte carlo quanta ( radioactive pellets , @xmath0-packets , @xmath2-packets )
is adopted throughout .      to obtain light curves for a pre - specified supernova model ,
the code undertakes the following steps .
a computation domain is assigned which is large enough to encompass the physical extent of the model .
this is divided into a number of grid cells each of which is assigned an initial mass - density based on the input supernova model .
pellets of radioactive material are placed in these cells , again in accordance with the chosen model .
the radiative decays of these pellets by emission of @xmath0-rays and the subsequent propagation and thermalisation of the @xmath0-packets is followed via a monte carlo calculation .
a grey - opacity treatment is adopted for the propagation of radiation ( @xmath2-packets ) in all other wave - bands .    during the monte carlo calculation ,
estimators are determined for various properties of the radiation field in every grid cell ( see below ) .
no information is retained regarding the variations of these properties within individual grid cells .
after the monte carlo simulations are complete , the estimators are used to determine observer frame emissivities . a formal solution of the radiative transfer equation is then performed using these emissivities to determine the emergent radiation field .
the code currently produces two classes of output : `` bolometric '' ( ultraviolet - optical - infrared , uvoir)-light curves which are obtained from the behaviour of the @xmath2-packets , and @xmath0-ray spectra and light curves which are derived from the @xmath0-packets . for the applications described in sections  3 and 4 , only the @xmath2-packet light curves are needed , however the analogous means by which @xmath0-ray properties are obtained are also described here for future reference .      following lucy ( 2005 ) , a regular cartesian grid which expands with time is used .
however , we allow the grid to expand continuously rather than in discrete jumps at the end of a each time - step .
since the supernova ejecta is assumed to be in homologous expansion , this adds little complexity to computing the propagation of the monte carlo quanta and removes the need to check whether they skip across boundaries due to the modification of the grid at the start of each time - step . in the limit of small time - steps , this modification has no effect on the results of the calculation .
the code described by lucy ( 2005 ) obtained light curves by directly counting the number of quanta that escaped the computational domain during each time - step .
however , as discussed by lucy ( 2005 ) , that approach is not the most efficient with regard to minimizing monte carlo noise .
significantly higher quality spectra and light curves can be obtained by using the paths of the monte carlo quanta to compute estimators for the emissivities in the model which can then be used , _ post hoc _ , to obtain the intensity via a formal solution of the radiative transfer equation .
methods using such an approach have previously been used and shown to be highly successful ( see e.g. lucy 1999 ) .    for spherical models  such as the test case used by lucy ( 2005 )  the method of counting packets can always be used reliably , its inefficiency being countered by the use of a large number of quanta . for complex 3d structures , however , this approach becomes unsatisfactory since the angular dependence of the escaping radiation field can only be addressed by angular binning of the quanta .
this has further negative impact on the signal - to - noise of the computed spectra and rapidly becomes prohibitive if more than a few angle - bins are to be considered .
in contrast , the use of emissivity estimators and a formal solution of the radiative transfer equation allow the spectra and light curves to be computed correctly for individual lines - of - sight to the supernova without introducing extra monte carlo noise
. thus such methods are strongly favoured for the study of any models which depart from spherical symmetry .
the implementation of these methods for the calculation of uvoir light - curves and @xmath0-ray spectra are described in the next two sub - sections .
the light - curve for uvoir radiation as viewed by a distant observer in direction @xmath3 is given by :    @xmath4    here , the integral is performed over the plane perpendicular to @xmath3 and @xmath5 is the emergent intensity of a ray which is destined to reach the observer at time @xmath6 .
the intensity depends on @xmath7 , the impact - vector of the ray given by @xmath8 where @xmath2 is any position on the ray trajectory .
@xmath9 is the luminosity the observer would imply were they to assume that the emission was isotropic .
the true ( unobservable ) luminosity of the supernovae , @xmath10 , is given by    @xmath11    this @xmath10 is the quantity obtained by directly summing all emerging packets in a monte carlo simulation .    to evaluate equation  1 for a particular viewing direction ( @xmath3 ) at a particular time
( @xmath6 ) , a large sample of rays ( typically @xmath12 ) are chosen such that they cover the complete range of @xmath13 required by the size of the supernova ejecta .
these rays are launched simultaneously from a plane perpendicular to ( @xmath3 ) behind the supernovae having initially zero intensity .
the time of launch is chosen such that the ray with @xmath14 crosses the coordinate origin at time @xmath6 . for each ray , the emergent intensity is determined by numerically solving the radiative transfer equation along its trajectory :    @xmath15    where @xmath16 is the element of pathlength along the ray while @xmath17 and @xmath18 are the observer frame , direction- , time- and position - dependent emissivity and opacity , respectively .
since the trajectory is a light ray , the position and the time @xmath19 are related by @xmath20 .
the emergent intensity , @xmath21 , is the value of @xmath22 at the point where the ray trajectory finally leaves the supernova ejecta .    in all the calculations presented here ,
the uvoir opacity per unit density in a grid cell is assumed to be constant and isotropic in the co - moving frame .
thus the observer frame opacity , @xmath23 can be readily calculated .
in contrast , the emissivity @xmath24 is not known _ a priori _ and is obtained from a monte carlo simulation .    in the current version of the code , there are two distinct uvoir emissivity source terms : one due to thermalisation of @xmath0-ray packets ( @xmath25 ) and one due to scattering of uvoir photon packets ( @xmath26 ) .
following lucy ( 2005 ) , estimators for these emissivities in the co - moving frame for a particular grid cell during a particular time step are obtained from the monte carlo energy packet trajectories using :    @xmath27    and    @xmath28    in equation  4 , the summation runs over all the trajectories of @xmath0-ray packets within the cell ( which has volume @xmath29 ) during the time - step ( which has duration @xmath30 ) ; @xmath16 is the trajectory length and @xmath31 is the packet energy determined in the observer frame .
the packet is travelling in direction @xmath32 and the velocity of the ejecta at the mid - point of the trajectory is @xmath33 .
the co - moving frame @xmath0-ray thermalisation opacity ( @xmath34 ) is frequency - dependent and includes contributions from both compton scattering and photoelectric absorption .
note that the compton term accounts only for the rate at which @xmath0-rays transfer energy to compton electrons ( which are assumed to thermalise in situ ) . in equation  5 ,
the summation is over uvoir packet trajectories and @xmath35 is obtained from the adopted uvoir opacity coefficient .
it has been assumed here that the emissivity is isotropic in the co - moving frame and that terms @xmath36 and smaller can be neglected .
note that although the co - moving frame emissivity ( @xmath37 ) is isotropic , the observer frame emissivity ( @xmath17 ) is not , owing to the doppler terms in the transformation between frames .
the scheme described thus far is valid provided that the grid cells are individually optically thin . in practice , this condition is violated at early times when the ejecta is compact and dense . at such times , it becomes unacceptable to assign a uniform emissivity to each grid cell .
this can be overcome , however , by weighting each contribution to the emissivities with a factor which accounts for the probability of energy absorbed and re - emitted during the related trajectory escaping to infinity .
the weighting factor used here , which must be applied individually to each contributing pathlength included in the sums in equations 4 and 5 is    @xmath38    where @xmath39 and @xmath40 is the total @xmath2-packet optical depth from the mid point of the trajectory @xmath16 to the edge of the supernova in the direction @xmath3 .
this form of the weighting factor is valid provided that _ either _ the total optical depth across a grid cell is small _ or _ that the contributing pathlengths ( @xmath16 ) are all small compared to the typical length scale on which the physical properties ( e.g. mass density ) of the model vary .
therefore , in calculations which make use of these weighted estimators , the pathlengths that energy packets can travel in a single step are not permitted to exceed a predetermined maximum step size , @xmath41 . for each time step in the calculation ,
@xmath41 is set to _ either _ one - tenth the width of a grid cell _ or _ the distance corresponding to @xmath42 in the densest grid cell ; the larger of these two values is chosen .
furthermore , when the weighted emissivities are used , the opacity term in equation  3 must be neglected since the probability of the emitted energy escaping has already been addressed .
this weighting overcomes the problem of having optically thick cells at early times but it has two significant drawbacks .
first , the required computations of @xmath43 are time demanding and lead to a substantial increase in code execution time . secondly ,
since the weighting factor is angle dependent , separate sets of estimators are required for each viewing angle that is to be investigated ( in contrast , the unweighted estimators are independent of @xmath3 ) .
therefore , this method should only be used for computing the light - curve at early times  when the density in the grid cells is high  and be replaced with the unweighted scheme once the opacity of individual cells has dropped sufficiently .      to test the implementation of the method described above
, the code has been applied to the test model used by lucy ( 2005 ) .
this is a spherical model , based on that used by pinto & eastman ( 2000a ) , having total mass @xmath44 , @xmath45ni mass @xmath46 and a maximum velocity of @xmath47  km  s@xmath48 .
the distribution of ni is centrally peaked .
following lucy ( 2005 ) , the grey absorption cross - section of 0.1  @xmath49  g@xmath48 is adopted for uvoir radiation , and the photoelectric absorption coefficient is taken from ambwani & sutherland ( 1988 ) , adopting a mean value of @xmath50 .
a 100@xmath51 cartesian grid is adopted for this test and @xmath52 monte carlo packets were used in the calculations .    the light curve obtained by lucy ( 2005 ) for this model
is shown in figure 1 ( solid line ) .
this agrees very well with the light curve computed here by directly counting the number of monte carlo packets escaping this model ( diamond symbols in figure  1 ) .
it also matches well with the light - curve computed for a particular viewing angle using the weighted estimators described in section 2.4 ( the particular viewing direction was chosen randomly here ) .
however , when the estimators are not weighted ( i.e. if equations 4 and 5 are used directly ) , the computed light curve is significantly overestimated at early times ( @xmath53 days , in this case ) . as discussed in section 2.4 , this is because the grid cells of the model are optically thick at early times ; the maximum optical depth across one grid cell is plotted , as a function of time , in figure  2 .
this shows that significant errors arise if the unweighted estimators are used when this optical depth is @xmath54 .
the @xmath0-ray spectrum is obtained following the same principles as that used for the uvoir light curve described above .
the important difference being that the treatment of @xmath0-rays is fully frequency dependent .    to deal with the frequency dependence ,
a grid of frequency points is used to divide the spectrum into small frequency intervals .
one point in this grid is set to the rest frequency of each of the radioactive @xmath0-ray emission lines in the range of interest and the remaining points are spaced logarithmically between . to compute the spectrum , the same scheme of solving the radiative transfer equation along a set of rays trajectories through the model is used . here , however , the radiative transfer along each trajectory is computed multiple times , once for each frequency point in the frequency grid , thereby determining the emergent radiation field as a function of both frequency and time .
as for the grey - computations described above , the opacity term in the radiative transfer equation is known ( the sum of compton and photoelectric terms ) .
there are again two emissivity terms which need to be considered .
the first , direct emission of @xmath0-rays by radioactive decay can also be expressed analytically for every grid cell in terms of the half - lives of the radioactive isotopes and their initial concentrations in the cell .
the second emissivity term is due to compton down - scattering .
the treatment of this term requires that both the angular- and frequency - dependence of the compton process be considered .
it is determined via a set monte carlo estimators ( one per frequency interval per grid cell per time step ) ; the estimator for the frequency interval @xmath55 in the grid cell @xmath56 for timestep @xmath57 is given ( in the observer frame ) by    @xmath58    where the sum runs over all @xmath0-ray trajectories which lie in cell @xmath56 ( having volume @xmath59 ) and have frequency appropriate for scattering into the frequency interval @xmath55 ( of width @xmath60 ) during the timestep @xmath57 ( which has duration @xmath61 ) .
@xmath62 is the number density of target electrons in the cell , @xmath31 is the observer - frame energy of the @xmath0-ray packet and @xmath63 is the compton differential cross - section for scattering into the direction of the line - of - sight ( @xmath3 ) in the observer frame .
this cross - section depends upon the angle between the trajectory and @xmath3 and is determined by applying the klein - nishina formula for the cross - section in the co - moving frame .
in this section , two toy models of elliptical supernovae are used to investigate possible observational consequences of large scale asphericity in supernova explosions .
such an investigation is motivated by observational evidence for global asphericity in snia obtained via polarimetry ( see e.g. howell et al .
2001 ; wang et al .
the origin of this asphericity is not well known : mostly likely it is determined by the details of the explosion process itself but may have its roots in the properties of a rapidly rotating progenitor . here
, however , the objective is not to gain insight to the physical origin of such a geometry but rather to study how it might effect both the amplitude and shape of observed light curves in comparison with spherical explosions .
earlier calculations of radiation transport in elliptical supernovae have been discussed by hflich ( 1991 ) . in that work ,
a multi - dimensional monte carlo code was also used . however , the treatment of energy packet generation and emission was simplified via the use of a parameterised photosphere in contrast to the full treatment of @xmath0-transport and deposition employed here .
a simple elliptical model has been constructed , closely related to the spherical model used as a test case in section  2 .
the adopted model has the same total mass and same @xmath45ni mass as the spherical model .
it is also assumed to be in homologous expansion and to have uniform density .
however , the maximum velocity ( and hence spatial extent ) is taken to be smaller in the @xmath64-direction than in the @xmath65- and @xmath66-directions ( symmetry is still assumed under rotation about the @xmath64-axis ) .
such a model is intended as a simple description for cases in which either the explosion mechanism or the properties of the progenitor lead to a large - scale ( low angular mode ) departure from sphericity .
two particular realisations of the model will be considered here .
for both the maximum velocity in the @xmath65- and @xmath66-directions was kept at 10@xmath67  cm  s@xmath48 .
the models differ in the chosen maximum velocity in the @xmath64-direction : for one , this velocity was fixed at @xmath68  cm  s@xmath48 , thereby giving an axis ratio of 2:1 ; while for the other @xmath69  cm  s@xmath48 was adopted to give an axis ratio of 5:4 .
the axis ratio of 2:1 may be regarded as extreme
it is comparable to the axis ratio that might be present in a very rapidly rotating progenitor ( see e.g. yoon & langer 2005 ) but there is little evidence to suggest that this would be preserved during an explosion . as an extreme case , however , this model is useful for providing a clear indication of the sense in which asphericity can affect the light curve .
the second ratio adopted ( 5:4 ) is comparable to the @xmath70 per cent asphericity found by howell et al .
( 2001 ) for sn1999by and is therefore more likely to be indicative of the propertied of real snia explosions .    as discussed in section  1 , the interest here is in probing the effect of a 3d treatment of the radiative transfer compared to spherically symmetric 1d calculations and so the simplifying assumption of a constant , grey - uvoir opacity is retained .
for all the calculations discussed in the section , this grey - uvoir absorption cross - section remains fixed at @xmath71  @xmath49g@xmath48 .
light - curves have been computed for observers viewing the ellipsoidal supernovae from infinity along both the major- and minor - axes .
these were computed using the weighted - estimators method described in section  2.4 .
the two light curves for the model with axis ratio 2:1 are plotted in figure  3 , along with the angle - averaged light curve for the same model ( i.e. the light curve obtained from the arithmetic mean of the light curves seen by a large number of observers from random viewing angles ) .
when viewed down the short - axis , the light curve is considerably brighter than as observed down the long - axis . around maximum light ,
the difference in brightness is approximately a factor of 3.5 . at later times ,
as the supernova becomes less optically thick , the difference becomes smaller  approximately a factor of 2.5 at around 50 days .
the light curve peaks slightly earlier when viewed down the long axis ( at @xmath72 days compared to @xmath73 days if viewed down the short axis ) .
also , the light curve decays more slowly if viewed down the long axis ; this is characterised by the bolometric @xmath74-values is the change in magnitude between maximum light and 15 days after maximum light .
] of 0.55 ( long - axis ) and 0.64 ( short - axis ) .    qualitatively similar , but quantitatively smaller differences are seen in the light curves computed for different viewing angles using the model with the axis ratio 5:4 . for this model ,
representative light curves are show in figure  4 .
the scale of the angular variation is determined by the interplay of several effects .
consider viewing an opaque ellipsoid with axis ratio 2:1 and uniform surface brightness as characterised by uniform surface temperature .
one would expect to find the flux to be twice as large when viewing along the short axis compared to the long axis , simply due to the increase in apparent surface area .
the computed ratio exceeds this for two closely related reasons .
first , the choice of uniform density in the models means that the column density to any particular 2:1 ellipsoidal surface within the model is always less when viewing down the short axis than the long .
thus the optical depth to the surface is smallest along the short axis making the light curve brighter when viewed from this direction .
secondly , contours of constant radiation density ( or equivalently emissivity ) do not exactly follow the 2:1 ellipsoidal geometry of the model ; in the outer regions there can be significant departures , always in the sense that the radiation energy density is highest at the points of intersection by the short axis .
this is illustrated in figure  5 where the variation of the @xmath2-packet emissivity ( @xmath75 ) is shown along the ellipsoidal axes .
a higher radiation energy density on the short axis makes the light curve brighter when viewed down that axis and so enhances the variation in brightness with viewing angle .
the origin of this effect also lies in the lower column densities along the short axis . at late times
, a quasi - static description of the radiation field becomes valid since packets are not trapped in the model for a significant number of timesteps . under such circumstances , the lower opacity along the short axis means that photons preferentially diffuse in that direction , making the energy density on an 2:1 ellipsoidal surface peak at the points of intersection with the short axis . at earlier times , this effect is enhanced by the time - dependent nature of the calculations ; fewer packets manage to reach outer ellipsoidal shells than would be predicted in a quasi - static description because they have had insufficient time to diffuse far enough .
the difference in diffusion time means that this affects the energy density along the long axis more significantly than along the short axis and thus acts to enhance the effect expected from the quasi - static case .
as time passes , the shapes of the contours of radiation energy density evolve slightly and the outward decline becomes somewhat less steep ( this can be seen by comparing the @xmath76  days and @xmath77  days results in figure  5 ) .
however , throughout the time range considered here , a departure from 2:1 ellipticity remains .    since
both the effects described above are the result of angular variations in the optical depth to 2:1 ellipsoidal surfaces , both persist while the ejecta remains optically thick .
they slowly decrease in strength during the decay phase as expansion causes the optical depths to drop ; at very late times the entire ejecta will become optically thin to uvoir radiation such that the light curve will become independent of viewing angle .
for the model adopted here , however , this nebular phase will not begin until several hundred days later than the times considered here .
the viewing - angle dependency of the light curves computed from these simple models may have some interesting ramifications for understanding the observed properties of snia light curves . as pointed out by wang et al .
( 2003 ) in the context of sn2001el , directional dependence of the luminosity as predicted for elliptical models of supernovae would lead to dispersion in the observed peak magnitudes ( based on the earlier work by hflich 1991 and their implied asphericity of @xmath78  per cent for sn2001el , they speculate that this dispersion would be around 0.1  mag ) . the results obtained here support this argument and indicate that if the degree of asphericity were larger in some cases ( e.g. sn 1999by ; howell et al .
2001 ) the spread in the peak magnitude could be greater , @xmath79  mag .
furthermore , the full light curves computed here allow this dispersion relative to the known relationship between light curve shape and peak luminosity to be examined .
this trend , the so - called `` phillips relation '' following phillips ( 1993 ) , expresses the negative correlation between peak brightness and the @xmath74-parameter measured from observed light curves of snia .
for both the models considered here , the variation of @xmath74 with viewing angle is in the opposite sense to the standard relationship .
this is illustrated in figure  6 where the six light curves shown in figures  3 and 4 are represented as points in the @xmath74-@xmath80 plane .
the gradient of the standard phillips relation ( describing the mean observed relationship between @xmath74 and @xmath80 ) in the b - band is plotted for comparison in the figure .
this effect would lead to a detectable scatter about the phillips relation and thus may have a major role to play in understanding the diversity of supernova observations : the results plotted in figure  6 would suggest that if snia explosions were moderately elliptical ( such that the model with axis ration of 5:4 were approximately applicable ) , viewing angle effects could explain a scatter of several tenths of a magnitude about the mean relationship .
significant caution must be applied in interpreting this result since the grey - treatment adopted here does not allow band - limited light curves to be studied for direct comparison with observations  quantitative differences may occur if the frequency - dependence of the opacity were taken into account . furthermore , the models used here have predicted an angular variation of the radiation energy density and thus , by implication , the temperature of the ejecta .
such a variation further contradicts the used of a uniform opacity and highlights the need for the consideration of more detailed micro - physics .
also , the particularly simple model chosen ( uniform density with centrally concentrated @xmath45ni and time - independent mean opacity ) produces significantly smaller absolute values of @xmath74 than are typically observed  thus further work using more realistic models of aspherical supernovae are needed .
modern snia explosion models predict complex , three - dimensional sub - structure ( e.g. reinecke et al . 2002 ; gamezo et al . 2003 ; rpke 2005 ; rpke et al .
2006 ) within the explosion .
in contrast , most models used to compute light curves for comparison with observation have a smooth , one - dimensional density-/composition - profile .    in this section , the effect of the predicted inhomogeneity on model light curves will be investigated . in principle , there are two classes of effect which are of interest here : first , inhomogeneity can change the viewing - angle averaged light curve compared to that of a spherically symmetric model . secondly , inhomogeneity could also lead to angular dependence of the light curve  the degree of angular dependence may be reasonably expected to depend on a combination of the dynamic range of the inhomogeneity in density and the physical length scale of the variations .
the length scale is relevant since if only small scale inhomogeneity is present it will tend to be averaged out when integrating over solid angle and so not introduce significant angular variation of the light curve .    given that the code used here adopts several simplifying assumption  most importantly , perhaps , that of a grey absorption coefficient for uvoir radiation
the emphasis here is on understanding and assessing the differential effect of introducing 3d structure by comparing with an equivalent 1d model using a fixed set of well - understood approximations .
such a calculation is an important first step in understanding the role of 3d structure and is a useful starting point for further work where the micro - physics is improved ( see section 5 for further discussion ) .
the model used for this test calculation is based on the 3d explosion model computed by rpke ( 2005 ) . for this model ,
rpke ( 2005 ) followed the hydrodynamics of an exploding white dwarf star ( with total mass 1.4  @xmath81 ) for 10 seconds on a schematic grid of 512@xmath82 cartesian cells .
only one spatial octant was simulated and symmetry under reflection was assumed to describe the remaining octants ( thus numerically only 256@xmath82 grid cells were used ) .
the distribution of mass - density and mass - fraction of iron - group elements in each grid cell at the end of their simulation was made available for this work .
light curves have previously been simulated from 1d representations of this model by blinnikov et al .
( 2006 ) .    in order to make the calculations here tractable  in terms of both computer memory required and photon statistics in each grid cell  the model adopted here uses only 170@xmath51 grid cells .
it was obtained from the 512@xmath51 grid by first removing the outermost cell from both ends of the grid in each of the three cartesian directions resulting in a 510@xmath51 grid .
the mass in the cells removed in this process was negligible .
the resolution was then reduced by a factor of 3 by subdividing the 510@xmath51 grid into 3@xmath51 blocks and replacing each block with a single cell whose density was equal to the mean density of the original 27 cells .
given the grey treatment of uvoir radiation currently adopted in the code , it is not necessary to specify the detailed composition of the material in each grid cell .
however , it is necessary to specify the initial distribution of @xmath45ni which provides the source of radiative energy .
the hydrodynamics code used by rpke ( 2005 ) provides as estimate of the fractional mass of iron - group elements in each grid cell but does not give a reliable estimate of the breakdown of this material into specific isotopes and elements .
the @xmath45ni mass - fractions used here were obtained by adopting a constant ratio for the mass of @xmath45ni to the total mass of iron - group elements in all grid cells .
this ratio was fixed to yield a total @xmath45ni - mass of 0.28  @xmath81 , as derived by kozma et al .
( 2005 ) for this model ; note that the nucleosynthesis calculations described by kozma et al .
( 2005 ) do not produce compositional information in sufficient detail to reconstruct the full 3d - distribution of @xmath45ni - mass owing to the modest ( 27@xmath51 ) number of tracer particles they used in comparison to the number of grid cells in the 3d model .    in order that the differential effect of the 3d structure could be assessed , a 1d comparison model was made by averaging the 3d model over spherical shells , taking care to conserve the radial distributions of total mass and @xmath45ni - mass . in order to illustrate the degree of inhomogeneity in the 3d model ,
figure  7 shows the velocity distribution of density in the spherically averaged model with points indicating the actual densities of individual grid cells in the 3d model .
this shows that across most of the velocity range of the model , there is a spread in density of at least a factor of three between different cells with similar velocities .
this is comparable to the dynamic range of density for a given velocity in the original hydrodynamical model ; thus one may be confident that the model used here contains a fairly reliable representation of the inhomogeneity implied by the hydrodynamics .
we continue to adopt a grey uvoir - absorption cross - section .
however , to correctly evaluate the influence of 3d structure on the light curve , it is necessary to consider that compositional inhomogeneity would cause the cross - section per gram to be a function of position .
it goes beyond the scope of this paper to undertake full calculations of the composition dependence of the opacity ; instead , a simple one - parameter description of the opacity is adopted , following mazzali & podsiadolwski ( 2006 ) .
they consider the opacity to be determined by the mass fraction of iron - group elements on the assumption that the opacity per gram is a factor of ten higher for the iron - group than for lighter elements .
adopting the same assumption , the uvoir - absorption cross - section per gram used in this section is given by    @xmath83\ ] ]    where @xmath84 is the mass fraction of iron - group elements , which varies from cell to cell , and the normalisation factor @xmath85 is chosen such that the mean value of @xmath86 in the ejecta is fixed to @xmath87  @xmath49 g@xmath48 .
although crude , this parameterisation captures the essential physics that the heavy elements dominate the mean opacity and only requires the compositional information which is directly available from the explosion model ( namely , the total iron - group mass ) .
note that , for ease of comparison with the constant-@xmath86 calculations used in the earlier sections of this paper , the time - dependence imposed on @xmath86 by mazzali & podsiadolwski ( 2006 ) is not included here .    in view of the simplifications used  in terms of both numerical resolution and
particularly the very simple treatment of opacity  the results obtained here should not be regarded as a definitive prediction of the radiation properties of the rpke et al .
( 2005 ) explosion model .
rather , the emphasis is on using a model which has a reliable representation of the degree of inhomogeneity in a real explosion model to understand the role played by the complex structure in the radiation transport .
the light curve obtained with the 3d model described above is shown as the solid histogram in figure  8 ; this is the light curve averaged over viewing angles . at peak , this light curve is dimmer than those plotted in figures  1 , 3 and 4 ; a consequence of the lower @xmath45ni - mass in the present model .
the light curve peaks earlier ( at @xmath88  days ) as a result of differences in the distribution of @xmath45ni with velocity : in the model used here , the @xmath45ni - distribution is less centrally concentrated than that adopted in sections  2.5 and 3 .    to establish the influence of the inhomogeneity on the light curve , the 1d model obtained by spherically averaging the 3d model ( see above ) was also used to compute a light curve .
this light curve is shown by the dashed curve in figure  8 .    from the figure
, it is apparent that the inhomogeneity causes the light curve to be a little brighter at early times and very slightly fainter during the decay phase .
this occurs because in the complex 3d medium a fraction of the @xmath2-packets find trajectories that preferentially pass through lower density cells and therefore encounter lower opacities than is possible in a spherically symmetric model with the same total mass .
this allows some @xmath2-packets to escape more quickly making the light curve brighter during the pre - maximum rise phase .
it also means that slightly fewer packets remain trapped in the ejecta which leads to the marginal dimming in the post - maximum phase .
the scale of the effect in this model is rather modest : at most it amounts to a difference of @xmath89  mag and this occurs only at very early times , @xmath78  days before maximum light . at around maximum light ,
the effect is limited to a change in magnitude of around @xmath90  mag ( 5 per cent in luminosity ) .    in accordance with `` arnett s rule ''
( arnett 1982 ; arnett , branch & wheeler 1985 )  which states that the emitted bolometric luminosity is roughly equal to the instantaneous rate of generation of radioactive luminosity at maximum light
the increase in peak luminosity resulting from the 3d structure is accompanied by a faster rise to maximum ; the peak luminosity occurs approximately one day earlier in the 3d calculations than the 1d case .
the 3d model was also used to look for viewing - angle dependency of the light curve .
light curves were extracted for a set of eight randomly chosen viewing angles and compared to the angle - averaged light curve .
no significant ( @xmath91 per cent ) departures from the average light curve were found ; this can be understood since large scale departures from sphericity are not present in this model , partially as a result of the assumed reflection symmetries ( see above ) .
it is interesting to note that the results from the model considered here show that inhomogeneity would act in the sense of pushing model results closer to observation ; for example , a recent comparison of light curves computed from 1d models with observations ( blinnikov et al .
2006 ) does suggest that the 1d models produce light curves which rise too slowly and underestimate peak luminosities ( see their figures 8  13 ) .
direct comparisons with the calculations of blinnikov et al .
( 2006 ) are not possible here owing to the considerably more complex treatment of opacity that they adopt , but it may reasonably be speculated that their rise times would also be shorter if multi - dimensional effects were incorporated .
the scale of the effect here is rather too small to directly have major implications for the confrontation of models and observations .
however , the effect is not completely negligible and it is plausible that the calculations presented here may underestimate the true scale of this phenomenon ; in particular , both the reduction in resolution from that used in the explosion model and the grossly simplified one - parameter treatment of the opacity may both suppress 3d density and composition effects that would be present in a complete treatment .
furthermore , only one explosion model has been considered here and this model is constructed from a simulation describing only one octant ( see above ) . in the future , it will be necessary to examine a range of models including fully 3d models where the greater complexity may enhance the effects .
we have described and tested a new 3d , time - dependent monte carlo code for modelling radiation transport in snia .
the code adopts the methods presented by lucy ( 2005 ) and also incorporates a scheme which uses monte carlo radiation field estimators to allow observables to be extracted for specific viewing angles ; this approach helps to suppress monte carlo noise and eliminates the need to obtain intensities by angular binning of emergent monte carlo packets .
this code has been used to investigate two classes of three dimensional effects in snia models .
first , two elliptical snia toy models were used to investigate how large scale asphericity might alter observable light curves .
as expected from previous simplified treatments ( e.g. hflich 1991 ) , it was found that light curves were brighter when viewed down the minor - axis than the major - axis .
the brightness enhancement is largest at early times and becomes smaller during the decline phase as the supernova ejecta becomes less optically thick . for a model with axis ratio comparable to that suggested by polarisation data for real snia ( howell et al . 2001 ; wang et al . 2003 ) , the differences in peak brightness and light curve shapes between viewing angles is detectable and in principle could play a role in the observed scatter about the mean relationship between light curve brightness and width .
secondly , a model with structure based on the results of recent three - dimensional hydrodynamical simulations was used to study the effect of inhomogeneities in both density and composition on light curves .
it was found that 3d structure could lead to light curves which are both brighter at early times and which peaked sooner after the explosion that is found from 1d models .
for the particular model considered , the effect was rather modest ( @xmath92 5 per cent around maximum light ) but study of other explosion models is required to quantify the possible diversity in this effect .    considerable further work is required to fully understand the role of 3d effects on radiation transport in snia .
the use of the grey approximation in the treatment of uvoir radiation is currently the greatest limiting factor to the practical applicability of these results . in calculations involving a uniform grey absorption coefficient
, the line - of - sight opacity depends solely on the total column density . under realistic conditions for snia ejecta ,
the opacity is dominated by spectral lines ( see e.g. pinto & eastman 2000b ) and is thus a strong function of frequency , velocity gradient , composition and the ionization / excitation state of the plasma . indeed , when strong lines dominate , the opacity has little direct sensitivity to the column density and is instead mostly a function of the density of spectral lines in frequency space .
furthermore , photon escape from highly opaque material is facilitated by frequency redistribution to those regions of the spectrum where there are relatively few spectral lines ( see e.g. pinto & eastman 2000b for a detailed discussion ) ; in the grey calculations this means of escape is not available and all the energy packets must burrow through the imposed opacity .
given this , it is likely that the grey approximation overestimates the role of geometry in determining photon propagation .
however , further calculations are required to determine and understand this in detail . in the context of the inhomogeneous model discussion in section  4 ,
a more realistic treatment of the role played by compositional inhomogeneity would be of particular interest ; only a crude attempt at describing the variation of opacity with composition has been used here ( equation  8) and a sophisticated treatment involving the very real differences between the atomic properties of different elements may lead to interesting effects .    fortunately , as discussed by lucy ( 2005 ) , the monte carlo method is well suited to incorporating realistic physics in a complex geometry which will allow the interplay of more detailed micro - physics and 3d structure to be studied in the near future ; this would involve a more sophisticated treatment of the opacity , perhaps similar to the methods employed by blinnikov et al .
( 2006 ) or kasen et al .
from such a calculations it would be possible to extract not only a bolometric light curve , but also light curves for specific photometric bands , an important step in the progression towards confrontation of theory and observation .
light curve calculations also need to be performed for a wider range of explosion models than the exploratory set discussed here .
this is necessary in order to quantify how the various effects might differ in scale depending on various properties of the explosion and so help to understand what role they might play in establishing both the average properties and diversity of sn observations .
eventually , it would also be valuable to undertake a similar assessment of the role of 3d effects in full nlte modelling of uvoir - spectra ; much of the greatest diagnostic power lies in spectral modelling  including both intensity and polarimetry  and thus it is important that commonly used modelling assumptions , such as that of spherical symmetry , are studied in the context of modern explosion models .
i am very grateful to f. rpke for making available the data used for the trial 3d model described in section  4 and for advice relating to hydrodynamical models and their interpretation ; to w. hillebrandt for useful comments and suggestions relating to various aspects of this work ; to l. lucy for making available the numerical results of his test light curve computations ( used for figure  1 ) ; and to d. sauer and r. kotak for helpful discussions .
i also thank the referee , p. pinto , for his careful reading of this manuscript and constructive comments and suggestions .
arnett w. d. , 1982 , apj , 253 , 785 + arnett w. d. , branch d. , wheeler j. c. , 1985 , nature , 314 , + 337 + ambwani k. , sutherland p. , 1988
, apj , 325 , 820 + benetti s. , et al . , 2004 , mnras , 348 , 261 + benetti s. , et al . , 2005 , apj , 623 , 1011 + blinnikov s. i. , rpke f. k. , sorokina e. i. , gieseler m. , + reinecke m. , travaglio c. , hillebrandt w. , stritzinger + m. , 2006 , a&a , 453 , 229 + branch d. , falk s. w. , uomoto a. k. , wills b. j. , mccall + m. l. , rybski p. , 1981
, apj , 244 , 780 + branch d , 1992 , apj , 392 , 35 + gamezo v. n. , khokhlov a. m. , oran e. s. , chtchelkanova + a. y. , rosenberg r. o. , 2003 , science , 299 , 77 + hflich p. , 1991 ,
a&a , 246 , 481 + hflich p. , wheeler j. c. , thielemann f. k. , 1998 , apj , + 495 , 617 + howell d. a. , hflich p. , wang l. , wheeler j. c. , 2001 , + apj , 556 , 302 + hoyle f. , fowler w. a. , 1960 , apj , 132 , 565 + kasen d. , thomas r. c. , nugent p. , 2006
, apj , 651 , 366 + kozma c. , fransson c. , hillebrandt w. , travaglio c. , + sollerman j. , reinecke m. , rpke f. k. , spyromilio j. , + 2005 , a&a , 437 , 983 + lucy l. b. , 1999 , a&a , 345 , 211 + lucy l. b. , 1987 , in eso workshop on the sn 1987a , ed .
+ i. j. danziger , 417 , european southern observatory , + garching + lucy l. b. , 2005 , a&a , 429 , 19 + maeda k. , mazzali p. a. , nomoto k. , 2006 , apj , 645 , 1331 + mazzali p. a. , lucy l. b. , 1993 , a&a , 279 , 447 + mazzali p. a. , podsiadlowski ph .
, 2006 , mnras , 369 , l19 + niemeyer j. c. , hillebrandt w. , woosley s. e. , 1996 , apj , + 471 , 903 + nomoto k. , thielemann f .- k . , yokoi k. , 1984 , apj , 286 , + 644 + mller e. , arnett w. d. , 1986 , apj , 307 , 619 + phillips m. m. , 1993 , apj , 413 , l105 + pinto p. a. , eastman r. g. , 2000a , apj , 530 , 744 + pinto p. a. , eastman r. g. , 2000b , apj , 530 , 757 + reinecke m. , hillebrandt w. , niemeyer j. c. , 2002 , a&a , + 391 , 1167 + rpke f. k. , 2005 , a&a , 432 , 969 + rpke f. k. , gieseler m. , reinecke m. , travaglio c. , + hillebrandt w. , 2006 , a&a , 453 , 203 + wang l. , et al .
2003 , apj , 591 , 1110 + yoon s .- c .
, langer n. , 2005 , a&a , 435 , 967 + [ lastpage ]


core collapse supernovae are dominated by neutrinos , therfore many supernova properties may depend on the nature of neutrino - nucleon interactions .
this provides an opportunity to study characteristic symmetries and features of the standard model weak interactions .
supernovae may provide macroscopic manifestations of charge conjugation @xcite or parity violation @xcite .    in the laboratory at high energies ,
antineutrino - nucleon cross sections are systematically smaller than neutrino - nucleon cross sections .
this is related to charge conjugation and parity violation in the standard model .
nevertheless , most core collapse supernova simulations still use the _ same _ lowest order cross section @xmath6 for both neutrinos and antineutrinos , @xmath7 .
\label{1}\ ] ] here @xmath8 is the ( incoming ) neutrino energy , @xmath9 the scattering angle and the vector @xmath10 and axial vector @xmath11 coupling constants are listed in table i.    in this paper , we discuss free space corrections to eq .
( [ 1 ] ) from nucleon recoil , weak magnetism , strange quarks and single nucleon form factors .
these corrections may be important for supernova simulations because they are present at all densities .
in contrast , the density dependent nucleon correlations discussed by @xcite and @xcite are important only at very high densities , well inside the neutrino sphere .
electron antineutrinos were detected from supernova 1987a via their capture on protons .
weak magnetism reduces this cross section and impacts the deduced neutron star binding energy and antineutrino temperature .
weak magnetism also changes the rate of antineutrino capture in the neutrino driven wind above a protoneutron star .
this may change the electron fraction @xmath12 ( ratio of electrons or protons to baryons ) and nucleosynthesis yields in the wind @xcite .    however , to confirm this change in @xmath12 , one should perform supernova simulations including weak magnetism because this may change the neutrino luminosities and emitted spectra which will also impact @xmath12 . in this paper , we present simple formulae so that these corrections can be incorporated in simulations .
we also give estimates of the change in emitted electron antineutrino spectrum and the differences expected between the spectrum of mu or tau antineutrinos and neutrinos .
the reduction in opacities from recoil and weak magnetism should increase luminosities of both neutrinos and antineutrinos as more energy is transported to the neutrino sphere .
this could increase heating behind the shock and more then compensate for the smaller antineutrino absorption cross section .
weak magnetism should also reduce the cooling from positron capture on neutrons further increasing the net heating behind the shock .
these changes in heating could make a simulation more likely to explode .
finally , weak magnetism allows mu and tau antineutrinos to escape the star faster than mu and tau neutrinos .
this should lead to a large muon number ( number of mu neutrinos minus antineutrinos ) or tau number of up to @xmath13 for the hot protoneutron star @xcite .
this is accompanied by a nonzero chemical potential for mu and tau neutrinos .
accurate treatments of neutrino transport are now feasible based on the boltzmann equation @xcite .
therefore , it is important to systematically improve neutrino opacities . in section
ii , we catalog a number of possible opacity corrections .    in section iii ,
we present exact neutrino - nucleon cross sections for both charged and neutral currents .
these are accurate to all orders in @xmath14 where @xmath8 is the neutrino energy and @xmath15 the nucleon mass .
we also present corrections from single nucleon form factors and strange quark contributions in the nucleon . in section
iv we discuss recoil and weak magnetism corrections to mu and tau neutrino , energy and lepton number fluxes and spectra .
section v discusses corrections for electron antineutrinos .
we also discuss possible changes in electron fraction @xmath12 and heating rates .
finally , we conclude in section vi that these weak magnetism and recoil corrections should be incorporated in future supernova simulations .
recently , accurate numerical algorithms have been developed to solve the boltzmann equation for neutrino transport @xcite .
the high accuracy of these simulations now warrants systematic improvements in the neutrino opacities .
therefore , we discuss a number of opacity corrections .
the corrections listed in table ii can be divided into two main groups .
numbers 1 through 5 are classified as model independent because they can be calculated exactly and are independent of the model used to describe dense matter .
in contrast , corrections 6 through 13 are model dependent .
not only do these corrections depend on the model , but it is important that these corrections be consistent with the equation of state .
corrections 1 , phase space , and 2 , matrix element are perhaps the simplest .
these are the only ones that are present in the limit of very low densities .
they represent more accurate calculations of free space neutrino - nucleon scattering . by phase space we mean corrections from the @xmath16 value of the reaction such as the neutron - proton mass difference or because the outgoing neutrino momentum is slightly different from the incoming momentum .
in this paper we focus on number 2 , corrections to the matrix element .
recoil corrections , 2a , arise because the nucleon is not infinity heavy and recoils slightly from the neutrino .
weak magnetism , 2b , arise from the parity violating interference between the weak magnetic moment of a nucleon and its axial current .
weak magnetism is important because it has opposite sign for neutrinos and antineutrinos .
thus weak magnetism increases the opacity for neutrinos and decreases the opacity for antineutrinos .
the first effects of the dense medium are 3 , pauli blocking of some outgoing nucleon or electron states , and 4 , the fermi or thermal motion of the initial nucleon or electron .
these can be calculated exactly .
the final model independent correction in table ii is 5 , coulomb corrections .
if heavy nuclei are present , the opacity may be dominated by coherent elastic scattering from nuclei . however , these nuclei have strong coulomb interactions and form a strongly correlated classical liquid
. the correlations between ions can greatly reduce the opacity for low energy neutrinos that have wavelengths comparable to the distance between ions .
this correction is model independent because the coulomb interaction is known .
for example , ref .
@xcite calculated the exact static structure factor of the ions with a simple monte carlo simulation . in principle , coulomb correlations are also present for neutrino electron scattering @xcite .
corrections 613 depend on the model used for the strong interactions between nucleons .
the mean field approximation , 6 , is perhaps the simplest way to treat the interactions .
nonrelativistic mean field models incorporate mean field effects in the density of states with an effective mass .
relativistic mean field models have both lorentz scalar and vector mean fields and the scalar field changes the effective dirac mass .
it is important that the opacities be consistent with the model used for the equation of state ( eos ) .
this can be achieved by using linear response theory @xcite to calculate the response of the medium to a neutrino probe .
the random phase approximation ( rpa ) or sum of ring diagrams provides the linear response of a mean field theory ground state .
therefore , if one uses a mean field eos and calculates the opacities in rpa , where the rpa effective interaction is the same as that used for the eos , then the opacities and eos will be consistent .
for example , the relativistic rpa opacities of ref .
@xcite are consistent with their mean field eos .
we note that the nonrelativistic calculations of burrows and sawyer @xcite include corrections 1 , 37 however they do not contain any matrix element correction , 2 .
likewise the relativistic rpa calculations of @xcite contain corrections 1 , 37 and 2a recoil but do not contain most of 2b weak magnetism . in the present paper
we discuss weak magnetism corrections .
in later work we will incorporate these into full nonrelativistic and relativistic rpa calculations .
there are many @xmath17 correlations beyond a simple rpa approximation , 8 , for example those responsible for @xmath18 @xcite .
in addition there are meson exchange current corrections , 9 .
the stony brook group has done some work calculating corrections from hyperons or other strange hadron components , 10 , @xcite .
more work could be done calculating opacities for meson condensates or other exotic phases of dense matter , 11 .
if superfluid or superconducting phases are present there will be pairing corrections to the opacity , 12 .
these can be important for the late time cooling of a neutron star @xcite or perhaps if there is a color superconducting phase @xcite .
finally there could be important corrections to opacities in nonuniform matter , 13 .
for example , in the low density pasta phases present in a neutron star inner crust or for a nonuniform meson condensate@xcite .
all of the above assumed magnetic fields do not play important roles .
if magnetic fields are important , they could substantially complicate the opacity , 14 , see for example @xcite .
magnetic fields might be important in jets or to break the spatial symmetry to explain neutron star kicks @xcite .    to conclude this section , now that accurate neutrino transport is possible
, it is important to systematically improve the neutrino opacities .
one should include weak magnetism and other free space corrections along with the many density dependent effects .
furthermore , these free space corrections can be included exactly and without any model dependence .
in this section we discuss recoil , weak magnetism , nucleon form factor , and strange quark corrections to the zeroth order cross section , eq .
. then we present differential , transport and total cross sections that are exact to all orders in the neutrino energy over the nucleon mass @xmath14 .
perhaps the simplest correction comes from the phase space of the scattered lepton because its energy @xmath19 is less than the incident energy @xmath8 .
if we use full relativistic kinematics for the nucleon and we assume it is originally at rest , then , @xmath20 \equiv   \eta k ,   \label{2}\ ] ] @xmath21^{-1 } , \label{3}\ ] ] with @xmath22 .    we expect corrections to depend on the small parameter , @xmath23 to the first order in @xmath24 , @xmath25
. \label{5}\ ] ] if one uses the phase space with @xmath19 but still uses the matrix element for @xmath26 then eq . (
[ 1 ] ) becomes , @xmath27 ,   \label{6}\ ] ] @xmath28 @xmath29 .
\label{8}\\\ ] ] we note that burrows and sawyer s work @xcite in the limit of low densities reduces to @xmath30 , see appendix i.      if one evaluates the matrix element to first order in e and includes the phase space corrections one has , @xmath31 .
\label{9}\ ] ] we note that the full recoil correction in eq . ( [ 9 ] ) is 50% larger than that for @xmath30 in eq .
( [ 8 ] ) .
equation ( 9 ) still neglects the parity violating interference between axial and vector currents .
this interference is dominated by weak magnetism .
the full weak vector current @xmath32 has dirac @xmath33 and pauli or tensor @xmath34 contributions , @xmath35 where @xmath36 is the momentum transfered to the nucleon .
the weak current @xmath32 is related by cvc ( conserved vector current ) to the electromagnetic current .
therefore the large anomalous magnetic moments of the proton and neutron give rise to a large @xmath37 .
the couplings @xmath10 and @xmath37 are collected in table i.    the cross section to the first order in @xmath24 is , @xmath38 , \label{11}\ ] ] with plus sign for @xmath39 and minus sign for @xmath40 .
weak magnetism increases @xmath39 and decreases @xmath40 , while recoil decreases both @xmath39 and @xmath40 cross sections .
therefore weak magnetism and recoil corrections approximately cancel for @xmath39 but add for @xmath40 .
note , burrows et al .
@xcite , see also for example @xcite , include weak magnetism for the charged currents but neglect it for neutral currents .
equation ( 11 ) is a good approximation for @xmath39
. however it can become negative for large @xmath40 energies . to cure this problem
we include corrections to all orders in @xmath24 .
the exact cross section can be calculated , see also @xcite , @xmath41 + c_a^2   \left [ 1+\eta^2 + \eta ( 1-x ) \right ]   \pm 2c_a ( c_v+f_2 ) ( 1-\eta^2)\nonumber \\ \left .
+ \eta^2 { e^2\over 2 } ( 1-x ) \left[f^2_2 ( 3-x ) + 4c_v f_2 ( 1-x ) \right ] \right\ } .
\label{12}\end{aligned}\ ] ]    in general the @xmath42 and @xmath43 terms are small because they only enter at order @xmath44 . however , their inclusion is important to ensure a positive cross section at high energies .
note , schinde considers the @xmath45 term but neglected all @xmath37 contributions @xcite .
the finite size of the nucleon can be included in eq .
( [ 12 ] ) by using appropriate form factors , @xmath46 that are functions of the four momentum transfer squared , @xmath47 for example , for @xmath48 , @xmath49 with @xmath50 1.26 and the axial mass @xmath51 1 gev .
form factors only enter at order @xmath44 .
they reduce the cross section by about 10% at @xmath52 , and even less at lower energies .
however , for completeness we include simple parameterizations of the form factors in appendix ii .
because the corrections are so small , we will ignore the form factors in the rest of this paper .
the opacity for @xmath53 from charge current reactions depends on the total cross section for @xmath54 , @xmath55 the zeroth order result is , @xmath56 while the correction factor from just phase space is , @xmath57 @xmath58+c_a^2   \left [ { 2e-1 \over e(1 + 2e ) } + { 1\over 2e^2 } ln ( 1 + 2e)\right]\right\ } / ( c_v^2 + 3c_a^2 ) , \label{19}\ ] ] @xmath59 the exact result is , @xmath60 @xmath61 , \label{22}\end{aligned}\ ] ] @xmath62 for @xmath3 using table i this is , @xmath63 while for @xmath53 we have , @xmath64 we note the main effect of weak magnetism and recoil is to reduce the @xmath65 opacity by the large amount in eq . ( [ 24 ] ) . for a typical neutrino energy @xmath66 20 mev this is a 15% reduction .
the full correction factor @xmath67 , is plotted in fig .
[ fig1 ] along with the lowest order forms in eqs .
( [ 24],[25 ] ) .
we see that the lowest order is fine for neutrinos but fails for antineutrinos above about 50 mev .      for mu and tau
neutrinos the scattering opacity depends on the transport cross section , @xmath68 note , strictly speaking this is true only in the limit @xmath69 nevertheless eq . (
[ 26 ] ) provides a good estimate when @xmath70 .
the zeroth order transport cross section is , @xmath71 the transport cross section including only phase space corrections is , @xmath72 @xmath73 + c^2_a \left
[ { 2e + \left(2e^2 -e-1\right ) ln \left(1 + 2e\right ) \over e^3 \left(1 + 2e\right ) } \right ] \right\ } \nonumber \\ / \left [ { 2\over 3 } \left ( c_v^2 + 5c_a^2 \right)\right ] , \label{29}\end{aligned}\ ] ] @xmath74 finally the exact transport cross section is @xmath75 @xmath76 \right .
\nonumber \\   + c_a^2 \left [ { 1+e\over 2e^3 } ln ( 1 + 2e ) -   { 10e^3 + 27e^2 + 18e + 3 \over 3e^2 ( 1 + 2e)^3}\right ] \nonumber \\
\pm   \left(c_v + f_2\right ) c_a \left [ { 1\over e^2 } ln \left ( 1 + 2e\right ) -   { 2 + 10e+ { 28\over 3}e^2 \over e ( 1 + 2e)^3 } \right ] \nonumber \\ + c_v f_2 \left[{1\over e^2}ln ( 1 + 2e ) -{2\over 3 }   \left ( { 3 + 15e + 22e^2\over e ( 1 + 2e)^3}\right)\right]\nonumber \\ \left
. + f^2_2 \left[{1\over 4e^2 } ln ( 1 + 2e ) +   { 8e^3 -22e^2 - 15e-3 \over 6e ( 1 + 2e)^3 } \right ] \right\ } / \nonumber \\ \left [ { 2\over 3 } \left(c^2_v + 5c_a^2 \right ) \right ] ,   \label{32}\end{aligned}\ ] ] @xmath77 the full correction , eq .
( [ 32 ] ) , is plotted in fig .
( [ fig2 ] ) .
corrections for @xmath78 are very close to those for @xmath79 elastic scattering .
we also plot the lowest order result , eq .
( [ 33 ] ) .
this fails for antineutrinos at high energies .
neutrino nucleon elastic scattering is sensitive to possible strange quark contributions in the nucleon .
strange quarks do not contribute to @xmath80 because the nucleon has no net strangeness
. however both @xmath81 and @xmath82 can be modified .
parity violating electron scattering constrains @xmath37 @xcite . in this section
we concentrate on @xmath11 because @xmath37 only contributes at order @xmath14 .
strange quark contributions are assumed to be isoscalar while the dominant contribution to @xmath11 is isovector .
therefore , strange quarks are expected to increase the cross section for @xmath79 scattering and decrease that for @xmath78 .
we write , @xmath83 with @xmath84 and the plus sign for @xmath79 and the minus sign for @xmath78 scattering .
the strange quark contribution @xmath85 is expected to be negative or zero .
the best limit comes from a brookhaven @xmath79 scattering experiment @xcite .
we assume , @xmath86 a better measurement of @xmath85 from a laboratory experiment would be very useful .
we make a simple estimate of strange quark effects starting from the zeroth order transport cross section in eq .
( [ 27 ] ) .
we consider the neutral current opacity for a mixture of neutrons and protons of electron fraction @xmath12 . for simplicity we set @xmath87 for @xmath79 scattering .
the ratio of the opacity with @xmath85 to that with @xmath88 is , @xmath89 + 10 g_a g_a^s(1 - 2y_e)+1-y_e\over 5g_a^2 + 1-y_e}.\ ] ] for @xmath90 , @xmath91 represents a 21% reduction in the opacity at @xmath92 or a 15% reduction at @xmath93 .
thus strange quarks could reduce opacities by 10 to 20% .
note , strange quarks do not contributed to charged current interactions .
in this section , we use the weak magnetism and recoil corrections for transport cross sections , eq .
( [ 32 ] ) to discuss muon and tau neutrino properties in supernovae .
first , we examine muon and tau lepton number currents .
because neutrinos now have shorter mean free paths than antineutrinos , there will be net @xmath94 and @xmath95 densities and nonzero @xmath94 and @xmath95 chemical potentials .
we then calculate the increase in the @xmath94 and @xmath95 energy flux because of weak magnetism and recoil .
this should increase the neutrino luminosities .
finally we examine corrections to the emitted spectra of @xmath94 , @xmath95 and @xmath96 @xmath97 .
we show that the @xmath98 are expected to be hotter than @xmath99 because of weak magnetism .    for simplicity , we neglect charged muons .
these could play some role at very high densities and or temperatures . without muons ,
our results are identical for muon and tau neutrinos .
therefore , we write @xmath5 in this section where @xmath100 could be @xmath101 or @xmath2 .
we start by considering high densities , well inside the neutrino sphere , where a simple diffusion approximation is valid .
earlier work @xcite calculated the lepton number current to first order in @xmath102 .
we extend this work to all orders .
the transport mean free path for a @xmath5 of energy @xmath8 is , @xmath103 with the zeroth order mean free path @xmath104 evaluated at an arbitrary reference energy @xmath105 , @xmath106^{-1 } , \label{35}\ ] ] and @xmath107 is the neutron number density .
note , for simplicity we assume pure neutron matter , @xmath108 . because the correction factor @xmath109 , eq .
( [ 32 ] ) , is so similar for @xmath110 and @xmath111 scattering our results are not expected to change much for @xmath112 .
the lepton number current for @xmath113 in a simple diffusion approximation is , @xmath114}. \label{36}\ ] ] here @xmath115 is the neutrino chemical potential .
for @xmath4 we have , @xmath116 this is now longer then @xmath117 because of the different weak magnetism correction factors @xmath118 .
the @xmath4 lepton number current has the same form as eq .
( [ 36 ] ) with @xmath115 replaced by @xmath119 , @xmath120 } .
\label{38}\ ] ] the total lepton number current is ,    @xmath121 } }   \right .
\left . - \int^\infty_0 { dk \over r^t_{\bar\nu } }   { 1 \over 1 + { \rm exp } \left [ ( k+\mu_{\nu_{x } } ) /t \right ] } \right ] .
\label{39}\end{aligned}\ ] ]    if @xmath122 the lepton number of the star ( number of @xmath5 - number of @xmath4 ) will rapidly change .
this buildup in lepton number gives rise to a nonzero chemical potential @xmath115 .
reference @xcite argued that the system rapidly reaches steady state equilibrium where @xmath123 .
we numerically solve for the chemical potential @xmath115 so that , @xmath124 this is shown in fig .
[ fig3 ] as a function of @xmath125 .
this chemical potential insures @xmath126 .    if we expand eq .
( [ 40 ] ) to first order in @xmath127 and use @xmath128 and @xmath129 expanded to the first order in @xmath14 one reproduces the lowest order chemical pot @xmath130 of ref .
@xcite , @xmath131 with @xmath132 for @xmath111 scattering .
the lowest order @xmath130 is seen in fig .
[ fig3 ] to be an excellent approximation to the exact result even at @xmath133 mev .
note , this nonzero @xmath115 can lead to muon number ( number of @xmath94 minus number of @xmath96 ) or tau number for the protoneutron star as large as @xmath13 @xcite .
we calculate the energy flux carried by @xmath5 by multiplying the integrand in eq .
( [ 36 ] ) by @xmath8 , @xmath134 we use eq .
( [ 41 ] ) for @xmath135 with @xmath136 and add the @xmath4 energy flux .
we assume the temperature @xmath125 has only a radial dependence so that , @xmath137
. \label{43 }   \end{aligned}\ ] ] in the absence of weak magnetism and recoil corrections , @xmath138 and the energy flux becomes , @xmath139 in fig .
[ fig4 ] we plot the ratio of the full energy current with weak magnetism and recoil to the zeroth order result @xmath140 .
_ we see that weak magnetism and recoil substantially enhances the energy flux .
_ this should raise the neutrino luminosity of the protoneutron star . at a temperature of 30 mev
the enhancement is over 50% .
many supernova simulations combine @xmath5 and @xmath4 into one effective species . in this approximation @xmath141 .
therefore , we evaluate @xmath142 in eq .
( [ 43 ] ) while setting @xmath143 .
this would correspond to using a ( geometric ) average for the correction factor , @xmath144 , \label{45}\ ] ] for both @xmath5 and @xmath4 .
we see that @xmath145 is slightly larger then @xmath140 .
this is because the chemical potential somewhat reduces the number of @xmath146 and hence their energy current .
even so , @xmath147 is a much better approximation then @xmath148 .
the above lepton number and energy fluxes are based on a diffusion approximation .
this is valid inside the neutrino sphere where @xmath149 is much less than the size of the system .
we now discuss weak magnetism and recoil corrections to the emitted @xmath5 and @xmath4 spectra .
we will use a simple monte carlo model of raffelt @xcite to estimate the emitted spectra .
we emphasize that this simple model should be checked against full simulations that include weak magnetism and recoil .
nevertheless , we expect the model to provide a first orientation .
raffelt @xcite discusses the formation of @xmath5 spectra . at high densities , reactions such as @xmath150 , @xmath151 , or @xmath152
keep @xmath5 in chemical and thermal equilibrium with the matter .
we define the energy sphere as the approximate location where these reactions , which all have small cross sections , become too slow to maintain thermal equilibrium .
next the @xmath5 propagate through a scattering atmosphere where neutrinos diffuse because of @xmath153 which has a much larger cross section .
however the small energy transfer in nucleon scattering is assumed to be too small to maintain thermal equilibrium . finally
, the neutrinos escape from the neutrino sphere .
note , the average spectrum for neutrinos and antineutrinos may be modified somewhat by nucleon - nucleon bremstrahlung in the scattering atmosphere @xcite .
however , we expect weak magnetism in nucleon scattering to provide an estimate for the _ difference _ between the spectrum of neutrinos and antineutrinos .
this should be explicitly checked in future work that includes weak magnetism for nucleon - nucleon bremstrahlung .
we model the energy sphere as a black body with temperature @xmath154 12 mev .
this high temperature is chosen so that the final temperature of the emitted @xmath5 agrees with detailed simulations .
alternatively , @xmath155 is the hot matter temperature well inside the neutrino sphere at densities high enough ( @xmath156 so that the @xmath5 are thermalized . we assume the same @xmath155 for both @xmath5 and @xmath4 because they are produced in pairs .
we next assume the @xmath5 propagate through a scattering atmosphere of optical depth , @xmath157 here @xmath158 is the thermally averaged optical depth in the absences of recoil and weak magnetism corrections .
raffelt assumes a boltzmann distribution for the initial thermalized neutrinos for which @xmath159 .
note , the optical depth is the thickness of the scattering atmosphere in units of the energy dependent mean free path .
a value of @xmath160 30 corresponds to one typical simulation @xcite .
raffelt calculates the survival probability @xmath161 for a @xmath5 of energy @xmath8 to make it through the atmosphere and escape .
otherwise it is assumed to be scattered back to the energy sphere and absorbed .
a fit to ref .
@xcite monte carlo results is , @xmath162^{-1 } \left [ 1- { 0.033 \over   1 + 1.5 ( \ell_\tau + .17)^2 + 0.5 ( \ell_\tau + .32)^6 } \right ] , \label{47}\ ] ] with @xmath163 .
therefore the final emitted spectrum is , @xmath164 for antineutrinos the optical depth is , @xmath165 and the @xmath4 spectrum is , @xmath166 the normalization @xmath167 is chosen so that @xmath168 .
this corresponds to steady state equilibrium and no net change in the lepton number .
note , @xmath167 is related to the chemical potential @xmath115 .
the @xmath5 and the @xmath4 spectra are shown in fig .
5 . both of these curves differ from a boltzmann distribution .
nevertheless the average energy can be characterized by a spectral temperature , @xmath169 table ii collects t values .
we find that t for @xmath4 is about 7% larger than that for @xmath5 almost independent of the choice of @xmath158 .
thus weak magnetism insures that the @xmath4 spectrum is hotter than that for @xmath5 . at high energies
@xmath170 is small .
therefore @xmath171 becomes significantly larger than @xmath172 .
this is shown in fig . 5 as an increasing ratio of @xmath173 to @xmath172 .
thus the high energy tail in the spectrum is expected to be antineutrino rich .
the difference between the @xmath5 and the @xmath4 spectra in fig .
[ fig5 ] is a macroscopic manifestation of charge conjugation , c , violation in the standard model . if the weak interactions had conserved c then @xmath174 . in principle , this difference is directly observable .
however , it may be very difficult to distinguish a detected @xmath5 from a @xmath4 . instead , a superposition of the @xmath5 and the @xmath4 spectra may be more easily observable .
weak magnetism , by separating the @xmath5 and the @xmath4 spectra and by increasing the high energy tail of the @xmath4 may lead to an observable broadening of the combined @xmath5 and @xmath4 spectrum .
the reduction in @xmath3 charged current interactions from weak magnetism and recoil can change the emitted @xmath3 spectrum , decrease the electron number current and increase the @xmath175 energy flux .
these changes , in turn , impact the @xmath12 in the neutrino driven wind which is a possible site for @xmath176process nucleosynthesis @xcite .
finally these changes could be important for analyzing the detected @xmath177 events from sn1987a or a future galactic supernova .
the increased @xmath3 mean free path from recoil and weak magnetism will increase the @xmath177 number current or decrease the total electron lepton number current @xmath178 .
this can increase the @xmath53 density in the star . for simplicity
, consider a situation where @xmath179 is near zero in a simulation without recoil or weak magnetism . then in steady state equilibrium , we expect recoil and weak magnetism to lead to a @xmath179 of order the @xmath115 from section iii eq .
. this chemical potential will lead to an increase in the electron fraction because , in @xmath180 equilibrium , @xmath181 .
this will lead to a small change in the electron fraction of order , @xmath182 .
\label{52}\ ] ] here @xmath183 is the electron chemical potential and @xmath184 the electron fraction without weak magnetism and @xmath12 is the new electron fraction . using eq .
( [ 41 ] ) for @xmath179 at t = 10 mev , @xmath185 = 15 mev one has , @xmath186 an increase of 10% .
note , eqs .
( [ 52],[53 ] ) slightly overestimate the change in @xmath12 because the charged current opacity for @xmath3 increases with @xmath12 .
nevertheless , weak magnetism could lead to a modest increase in @xmath12 .    in principle , weak magnetism will increase the lepton number diffusion time scale .
however this time scale is very sensitive to temperature .
we also expect weak magnetism to increase the @xmath3 energy flux in a similar way to fig .
4 . with weak magnetism ,
the protoneutron star should cool faster and the lower temperatures will decrease the lepton number diffusion time .
these contrary effects should be investigated in a full simulation .
we consider the following simple model @xcite to estimate the change in the @xmath3 spectrum with weak magnetism .
let the @xmath3 neutrino sphere be at a temperature @xmath187 without weak magnetism . for simplicity
we assume a boltzmann spectrum @xmath188 .
now we include weak magnetism and recoil described by @xmath189 .
this reduction in the opacity will shift the neutrino sphere inwards to higher densities and temperatures .
we define a new energy dependent temperature @xmath190 to include this shift , @xmath191^u t_0 . \label{54}\ ] ] here @xmath192 is the ratio of temperature to density gradients , @xmath193 , \label{55}\ ] ] for the protoneutron star at the neutrino sphere .
the extra -1 in the denominator follows because the optical depth involves a path integral from the neutrino sphere to infinity @xcite .
realistic values of @xmath192 could be @xmath194 0.25 to 0.35 @xcite .
the energy dependent temperature @xmath190 defines the emitted spectrum , @xmath195 if we expand @xmath196 , with @xmath197 and expand @xmath198 to first order in @xmath14 , @xmath199 e^{-k / t_0 } .
\label{57}\ ] ] we define a spectral temperature @xmath200 as , @xmath201 , \label{58}\ ] ] for @xmath192 = 0.3 and @xmath202 5 mev one has , @xmath203 thus weak magnetism and recoil can increase the emitted @xmath3 average energy by of order 10% .
it is important to check our simple model with full simulations .
figure 6 shows the full spectrum , eq .
( [ 57 ] ) .
we see that weak magnetism and recoil shift the strength to higher energies and increase the high energy tail significantly .
this is shown by a large ratio of the spectrum with weak magnetism to @xmath204 at high energies .
we now discuss some possible implications of weak magnetism on neutrino heating , electron fraction and neutrino cooling in the low density region outside the gain radius .
the gain radius is the point were cooling from neutrino emission balances heating from neutrino absorption and scattering
. the net amount of heating outside the gain radius could be important for the success of the explosion .
the absorption of @xmath3 outside the gain radius will be reduced by the smaller cross section from weak magnetism .
however the cross sections grows with the square of the energy so the new absorption rate may be proportional to , @xmath205 the first factor could increase the rate by about 20% while the second factor @xmath206 0.8 could decrease the rate by a similar amount .
therefore , the net rate of @xmath3 absorption might be little changed .
thus the electron fraction , which depends on the absorption rate , should not change greatly .
however , the @xmath3s absorbed have a higher energy so the net heating should increase by about 10% because of the higher average energy in eq .
( [ 59 ] ) .
this assumes the luminosity of the @xmath3 changed only because of the change in @xmath200 . if the luminosities of @xmath53 and @xmath3 increased further because of the increased energy fluxes in fig .
4 then the total heating could be even larger .
finally we discuss cooling from positron capture outside the gain radius , @xmath207 the rate for eq .
( [ 61 ] ) will be reduced by weak magnetism , @xmath208 at @xmath125 = 2 mev the factor in parenthesis reduces the cooling rate by 9% .
thus weak magnetism can change the @xmath3 luminosity to increase the heating while , at the same time , reducing the cross section to decrease the cooling .
these results should be checked with full simulations .
in this paper we examine recoil and weak magnetism corrections to @xmath39-nucleon interactions .
these are important because they are present at all densities , even at the relatively low densities near the neutrino sphere .
furthermore , the corrections are model independent .
we calculate simple , exact correction factors eqs . ( [ 12],[22],[32 ] )
to include recoil and weak magnetism in supernova simulations .
perhaps the most important effect of weak magnetism and recoil is to increase energy fluxes of both @xmath4 and @xmath3 antineutrinos .
we calculate , in a diffusion approximation , an increase of order 15% in the total energy flux for temperatures near 10 mev .
this should raise the neutrino luminosity .
weak magnetism and recoil will also change the emitted spectrum of @xmath0 and @xmath3 .
we estimate that @xmath4 will be emitted about 7% hotter than @xmath5 because @xmath4 have longer mean free paths .
likewise weak magnetism may increase the @xmath3 temperature by of order 10% .
this increase in temperature coupled with the increase in neutrino luminosity should overcome the reduced absorption cross section and increase the heating in the low density region outside of the neutrino sphere .
this , in turn , could be important for the success of an explosion .
we find large corrections .
however , supernova simulations are very complicated with many degrees of feedback . therefore , it is important to check our results with a full simulation that includes boltzmann neutrino transport and weak magnetism corrections .
we thank john beacom , hans - thomas janka and georg raffelt for extensive discussions of the physics .
we thank the institute for nuclear theory in seattle for its hospitality when this work was started and acknowledge financial support from doe grant de - fg02 - 87er40365 .
llll + reaction & @xmath10 & @xmath11 & @xmath37 + @xmath209 & @xmath210 0.035 & @xmath211 0.630 & @xmath212 1.019 + @xmath213 & @xmath214 & @xmath215 -0.630 & @xmath216 -.963 + @xmath217 & 1 & @xmath50 1.260 & @xmath218 3.706 + [ table1 ]    here @xmath50 1.260 , @xmath219 0.2325 , @xmath220 1.793 and @xmath221 1.913 .
lll + & & correction + 1 . & & phase space + 2 . & & matrix element + & a. & recoil + & b. & weak magnetism + & c. & form factors + & d. & strange quarks + 3 . & & pauli blocking + 4 . & & fermi / thermal motion of initial nucleons + 5 . & & coulomb interactions + 6 . & & mean field effects + 7 . & & nn correlations in rpa + 8 . & & nn correlations beyond rpa + 9 . & & meson exchange currents + 10 . & & other components such as hyperons + 11 . & & other phases such as meson condensates or quark matter + 12 . & & corrections from superfluid/ superconductor pairing + 13 . & & nonuniform matter + 14 . & & magnetic field effects + [ table2 ]    [ cols= " < , < , < , < " , ]     c.j .
horowitz and gang li , phys .
b443 * ( 1998 ) 58 .
c.j . horowitz and
gang li , phys .
d * 61 * ( 2000 ) 063002 .
a. vilenkin , 1979 , unpublished ; ap j * 451 * ( 1995 ) 700 . c. j. horowitz and j. piekarewicz , nuc
a640 * ( 1998 ) 281 .
phil arras , dong lai , phys rev .
d * 60 * ( 1999 ) 043001 .
sanjay reddy , madappa prakash , james m. lattimer and jose a. pons , phys .
c * 59 * , ( 1999 ) 2888 .
a. burrows and r.f sawyer , phys .
c * 58 * ( 1998 ) 554 ; phys .
c * 59 * ( 1999 ) 510 .
horowitz and gang li , phys .
lett . * 82 * ( 1999 ) 5198 .
horowitz , astro - ph/0108113 .
see for example , m. rampp and h .-
janka , ap jl * 539 * ( 2000 ) 33 , or s.w .
bruenn , k.r .
denisco , and a. mezzacappa , astro - ph/0101400 . c. j. horowitz , phys .
d * 55 * ( 1997 ) 4577 . c. j. horowitz , phys .
* 66 * ( 1991 ) 272 .
a. l. fetter and j. d. walecka , `` quantum theory of many particle systems '' , mcgraw - hill , ny 1971 , p 171 .
s. reddy , g. f. bertsch , and m. prakash , phys .
475b * ( 2000 ) 1 .
steen hannestad and georg raffelt , apj .
* 507 * ( 1998 ) 339 .
e. g. flowers , m. ruderman , and p. g. sutherland , apj .
* 205 * ( 1976 ) 541 ; a. d. kaminker , p. haensel and d. g. yakovlev , ast . and astrophys .
* 373 * ( 2001 ) l17 .
gregory w. carter and sanjay reddy , phys .
d * 62 * ( 2000 ) 103002 .
sanjay reddy , madappa prakash and james m. lattimer , phys .
d * 58 * ( 1998 ) 013009 .
stephen j. hardy , markus h. thoma , phys .
d * 63 * ( 2001 ) 025014 .
adam burrows , timothy young , philip pinto , ron eastman and todd a. tompson , apj . * 539 * ( 2000 ) 865 .
p. vogel , and j. f. beacom , phys .
d * 60 * ( 1999 ) 053003 .
ahrens et al . , phys .
d * 35 * ( 1987 ) 785 .
schinde , ap j suppl .
* 74*(1990 ) 249 .
b. mueller et al . , phys .
lett . * 78 * ( 1997 ) ; k. a. aniol et al . ,
* 82 * ( 1999 ) 1096 .
raffelt , astro - ph/0105250 .
m.j . musolf and t.w .
donnelly , nucl .
a546 * ( 1992 ) 509 .
in this appendix we discuss the low density limit of burrows and sawyer s calculations @xcite and show that it reduces to our phase space result , @xmath222 , in eq .
( [ 6 ] ) .
burrows and sawyer start by considering a pure vector interaction .
( 2 ) for the differential rate of neutrino scattering @xmath223 is , @xmath224\lambda^{00}s(q,\omega ) .
\label{a1.1}\ ] ] here the momentum transfered to the neutrons is @xmath225 and the energy transfered is @xmath226 .
the neutrino distribution is @xmath227 and @xmath228 is the neutrino tensor , eq .
( 3 ) of ref .
the dynamic structure function @xmath229 for the nuetrons is , @xmath230 2\pi \delta(\omega+\epsilon_p - \epsilon_{p+q } ) , \label{a1.2}\ ] ] where @xmath231 is the neutron distribution function and @xmath232 the energy of a neutron of momentum @xmath233 .
one can recover our eq .
( [ 6 ] ) by , ( a ) neglecting the final state neutrino pauli blocking , ( b ) integrating eq .
( [ a1.1 ] ) over @xmath234 , ( c ) assuming that the density is low enough so that @xmath235 in eq . ( [ a1.2 ] ) , @xmath236 @xmath237 here @xmath107 is the neutron density . finally ( d ) we assume , in the limit of low density , that the @xmath234 depdendence of @xmath229 is sharply peaked near , @xmath238 , with @xmath239 given by our eq .
( [ 2 ] ) . adding a simillar expression for the axial interactions gives , @xmath240 with @xmath222 given by our eq . ( [ 6 ] ) .
in this appendix we collect simple parameterizations of the single nucleon form factors for use in eq .
( [ 13 ] ) .
we write the weak form factors in terms of the electromagnetic dirac , @xmath241 , and pauli , @xmath242 , form factors for @xmath243 or @xmath244 .                      here , see eq .
( [ 14 ] ) , @xmath254 @xmath255 g / \left ( 1 + \tau \right ) ,   \label{a2.11}\ ] ] @xmath256 @xmath257 @xmath258 with @xmath259 @xmath260 and finally , @xmath261 these parameterizations are from ref .