the radial instability of the electron - degenerate , chandrasekhar - mass core marks the beginning of the final episode in the life of a massive star with zero - age main - sequence ( zams ) mass in the range @xmath3 .
collapse ensues and , once fully dynamical , separates the core into a subsonically homologously contracting inner core and a supersonically collapsing outer core . at nuclear density , the repulsive component of the nuclear force leads to a stiffening of the equation of state ( eos ) .
this stabilizes the inner core , which overshoots its new equilibrium , then rebounds into the outer core , launching a strong hydrodynamic shock wave from its edge .
this instant in time is referred to as _
core bounce_. the inner core has a mass of @xmath4 at bounce and this material becomes the unshocked core of the protoneutron star .
the shock formed at core bounce propagates into the outer core , but the dissociation of nuclei in the accreted material into neutrons and protons and the electron capture on free protons in the region behind the shock ( the _ postshock _ region ) soon sap its might , driving it into submission to the ram pressure of accretion .
the shock stalls and turns into a standing accretion shock that must be revived to unbind the stellar envelope and drive a core - collapse supernova explosion .
neutrinos play a pivotal and dominant role in stellar collapse and core - collapse supernovae .
neutrinos and antineutrinos of all flavors carry away the @xmath5@xmath6 ( @xmath7 ) of gravitational binding energy of the remnant neutron star over tens of seconds after core bounce .
aided by multi - dimensional fluid instabilities , they probably deposit , within a few hundred milliseconds after bounce , sufficient energy in the region behind the stalled shock to revive the shock and make typical @xmath5@xmath8 core - collapse supernova explosion ( @xcite and references therein ) . only hyper - energetic ( i.e. , @xmath9 ) explosions may require a different mechanism @xcite , e.g. , rapid rotation combined with strong magnetic fields , which may lead to energetic jet - driven explosions ( e.g. , @xcite ) .    for a galactic or near - extragalactic core - collapse supernova ,
neutrinos offer the unique possibility of directly observing the dynamics and thermodynamic conditions prevalent in the supernova core .
together with gravitational waves ( see , e.g. , @xcite ) they will herald the next nearby supernova possibly hours before any telescope sensitive to electromagnetic waves will notice the event . in the probable case that the next galactic supernova occurs in a dust - enshrouded region and/or close or behind the galactic center ,
the supernova may be impossible to observe in broad bands of the electromagnetic spectrum , making neutrino and gravitational - wave observations even more important .
the observation of neutrinos from sn 1987a in the large magellanic cloud @xcite confirmed the basic picture of core collapse and early protoneutron star evolution ( e.g. , @xcite and references therein ) , but the small number and poor timing of the observed interactions did not allow far - reaching and robust conclusions on core - collapse supernova dynamics and the involved neutrino physics , nuclear physics and astrophysics
.    the situation will be completely different when the neutrino burst from a galactic core - collapse supernova reaches current and near - future neutrino detectors on earth .
super - kamiokande @xcite , icecube @xcite , lvd @xcite , borexino @xcite , kamland @xcite , sno+ @xcite , no@xmath10a @xcite and others will together see many thousands of neutrinos from a core collapse event at @xmath11 @xcite .
distance estimates based on the observed neutrino flux @xcite , sky localization @xcite , and triggering of gravitational - wave searches by the reconstruction of the time of core bounce @xcite will likely all be possible .
well - timed high - statistics coincident neutrino observations will allow to probe in detail a broad range of supernova astrophysics , nuclear physics , and neutrino physics ( see @xcite for overviews ) .
fast characteristic temporal variations in the preexplosion neutrino fluxes would be tell - tale signs of multi - dimensional fluid instabilities in the postshock region @xcite and/or early postbounce ring - down oscillations of a rapidly spinning protoneutron star @xcite .
while artificially - driven , spherically - symmetric core collapse simulations suggest that a sudden deep drop of the accretion - driven component of the neutrino luminosity ( primarily in @xmath12 and @xmath13 ) within a few hundred milliseconds would indicate the onset of explosion ( e.g. , @xcite ) , the situation is likely different in nature as multidimensional simulations show a much shallower drop due to simultaneous accretion and explosion of material .
this leads to a persistent accretion luminosity even after the explosion has begun @xcite .
the spectral characteristics and long - term spectral evolution of the neutrino flux could provide important constraints on the nuclear eos @xcite and/or the spin of the progenitor core @xcite .
the ability of some detectors to distinguish interactions of different neutrino flavors would lead to constraints on the neutron - to - proton ratio in the neutrino - driven wind phase , allowing an observational test of core - collapse supernovae as potential sites for @xmath14-process nucleosynthesis @xcite .    the neutrino signature of core collapse , of the subsequent core - collapse supernova evolution and of the protoneutron star cooling phase ,
is invariably intertwined with neutrino oscillation physics .
the robustness of all of the above mentioned observational conclusions will depend on our understanding of the impact of neutrino flavor oscillations .
neutrinos propagating from their emission site to detectors on earth may experience ( _ i _ ) so - called vacuum oscillations driven by neutrino mass differences @xcite , ( _ ii _ ) oscillations mediated by a resonance in @xmath10@xmath15 scattering ( the mikheyev - smirnov - wolfenstein [ msw ] effect ; @xcite ) , and ( _ iii _ ) oscillations due to @xmath10@xmath10 scattering ( @xcite ; see @xcite for a review ) .
vacuum and msw oscillations are well understood and their outcomes depend essentially only on neutrino mixing parameters , in particular the neutrino mass hierarchy and the mixing angles .
the @xmath10@xmath10-scattering driven oscillations , on the other hand , have a non - linear hamiltonian that may lead to so - called collective oscillations with very complex spatial and temporal outcome that remains to be fully understood ( see , e.g. , @xcite and references therein ) .
however , a number of recent studies suggest that collective oscillations may be completely or at least partially suppressed in the preexplosion accretion phase of ordinary core - collapse supernovae @xcite , but see @xcite and @xcite for discrepant results .
provided that collective oscillations can be ignored in the preexplosion phase and that the @xmath16 mixing angle indeed has the large value suggested by recent measurements @xcite , the neutrino mass hierarchy may be inferred from the qualitative shape of the early postbounce neutrino signal @xcite .
a preexplosion accretion phase with suppressed collective oscillations would also offer the opportunity to probe the structure of the progenitor star on the basis of the observed neutrino signal .
the details of the preexplosion neutrino emission have been discussed carefully , e.g. , by @xcite and @xcite and we shall not repeat them here .
it is , however , necessary to outline its most salient features . for simplicity
, we neglect neutrino oscillations in the following
.    in core collapse and in the subsequent postbounce evolution , emission of @xmath12 and @xmath13 occurs via charged and neutral currents , while heavy - lepton neutrinos @xmath17 are created exclusively via thermal neutral - current pair processes . before core
bounce , only @xmath12 are emitted from electron capture in the collapsing core .
milliseconds after core bounce , the shock breaks out of the @xmath12 neutrinosphere ( where the optical depth is @xmath18 ) and a strong burst of @xmath12 is emitted for @xmath5@xmath19 from rapid electron capture on the freshly abundant free protons behind the shock .
@xmath20 are copiously created in the hot interior of the protoneutron star after bounce and begin to diffuse out , leading to a steep rise , quick leveling and subsequent slow decay of the @xmath20 luminosity ( @xmath21 ) .
@xmath13 production via charged - current positron capture is initially suppressed due to the high degeneracy of the electrons .
the latter is partially lifted after bounce at the moderate - density , hot edge of the protoneutron star and @xmath22 rises , reaching or surpassing the value at which @xmath23 levels off after the neutronization burst decays .
the subsequent preexplosion luminosity can roughly be split into a diffusive component from the core and accretion luminosity ( @xmath24 \dot{m } / r_\nu$ ] , where @xmath25 is an approximate neutrinosphere radius ) from or from above the neutrinosphere @xcite .
@xmath21 is primarily diffusive , while @xmath23 and @xmath22 are dominated by accretion . in general , @xmath26 , but @xmath27 .
@xmath20 have the lowest opacity , since they interact only via neutral currents .
they decouple from matter at the smallest radii and highest temperatures and thus have the highest average energies @xmath28 .
@xmath13 have a slightly lower opacity than @xmath12 , leading to the well established neutrino energy hierarchy in the preexplosion phase @xmath29 .
the mean energy of all species grows with increasing postbounce time , reflecting the recession of the neutrinospheres due to the contraction of the protoneutron star .
considering that the accretion luminosity will scale with the postbounce accretion rate @xmath30 , one would naturally expect an increase of the detected neutrino interactions with increasing mass of the stellar core .
since higher accretion rates correspond to more material compressing and settling more rapidly on the protoneutron star , the latter s outer regions will be hotter .
thus the thermal neutral - current emission will be enhanced , leading to higher luminosities and higher mean neutrino energies .
the variation of the preexplosion neutrino signal with progenitor star zams mass was first discussed by @xcite based on the pioneering simulations of @xcite and @xcite .
these authors provided total emission characteristics and spectra that show a systematic increase of total energy emitted in neutrinos and mean @xmath13 energy with zams mass in the range from @xmath31 .
@xcite , before sn 1987a , carried out simulations of a range of progenitor stars with zams mass in the range @xmath32 .
they found that the @xmath12 neutronization burst shows little dependence on the progenitor , due to the rather universal homologous collapse and bounce dynamics .
furthermore , they mentioned , though did not discuss in detail , that the luminosities and mean neutrino energies increase as a function of iron core mass ( and not zams mass ) .
a more detailed and clear physical discussion was provided by @xcite , who contrasted the predicted neutrino signal from the early postbounce phase in two different progenitor core models with neutrino observations of sn 1987a .
he noted that there are significant uncertainties in connecting a given zams mass to precollapse structure . instead of a progenitor with an associated zams mass , he considered a massive ( and high - entropy ) @xmath33-@xmath34 iron core model and a lower - mass ( and lower - entropy ) @xmath35-@xmath34 iron core model in his spherically - symmetric ( 1d ) neutrino radiation - hydrodynamic simulations .
he showed that the more massive core leads to a consistently higher @xmath13 luminosity in both the accretion and diffusion sectors .
the water cherenkov detectors that observed neutrinos from sn 1987a are most sensitive to the inverse beta decay ( ibd ) reaction @xmath36 .
@xcite predicted a factor of two difference in the integrated number of early ibd interactions between the massive and the low - mass core in these detectors .
he concluded that the neutrino signal observed by these detectors from sn 1987a was most consistent with the low - mass core .
@xcite , who carried out a parameter study of quasi - hydrostatic protoneutron star cooling , considering various initial masses , ad - hoc accretion rates , and different nuclear eos , found a similar trend .
he showed that more massive cores , higher accretion rates , and softer eos lead to stronger , higher - energy neutrino emission .
some of his strongest emitters were cases in which eventually a black hole was formed .
liebendrfer and collaborators carried out a sequence of studies of the progenitor dependence of the neutrino signal using modern general relativistic 1d radiation - hydrodynamics simulations @xcite .
they showed that the @xmath12 neutronization burst is indeed almost independent of progenitor structure ( as first suggested by @xcite ) .
they also qualitatively and quantitatively connected the evolution of the postbounce preexplosion luminosity to the postbounce accretion rate , but did not discuss observational implications .
their results were corroborated by similarly sophisticated subsequent studies of @xcite .
of these , @xcite presented the most comprehensive analysis and also compared between 1d and axisymmetric ( 2d ) results .
they found that in 2d , convection in the protoneutron star alters the structure of the latter , affecting the neutrino emission starting @xmath5@xmath37 after bounce , but preserving the overall systematics with accretion rate .
@xcite also were the only authors to suggest that the accretion - rate dependence of luminosity and total emitted energy in the preexplosion phase could be used to infer the structure of the progenitor .
the other studies , being focused on aspects such as neutrino oscillations , black hole formation , or the late - time post - explosion evolution , did not consider observational consequences .
@xcite , using a limited set of three progenitor models ( \{@xmath38 , @xmath39 , @xmath40}@xmath41 at zams ) , found similar systematics as the aforementioned studies , but also carried out an analysis of the expected signal in various neutrino detectors in the first @xmath42 after bounce .
they computed ibd interaction rates for their @xmath40-@xmath34 and @xmath38-@xmath34 and found a factor of two more ibd interactions for the former , which would allow a high - confidence distinction between these progenitors for a galactic core collapse event .
however , their @xmath39-@xmath34 model yielded a postbounce neutrino signal very similar to that of their @xmath38-@xmath34 model and would be indistinguishable by neutrino observations alone .
this suggests that zams mass is not a good parameter to describe presupernova stellar structure ( cf . , @xcite ) .    in this article
, we present a fresh look at the progenitor dependence of the neutrino signature in the preexplosion accretion phase of core - collapse supernovae .
we perform 1d general relativistic radiation - hydrodynamics core collapse simulations of 32 progenitor models from the single - star solar - metallicity presupernova model suite of @xcite and follow the postbounce preexplosion evolution for @xmath43 . in zams mass , these models range from @xmath0 to @xmath1 , but guided by the previous results discussed in the above , we choose not to parameterize our simulations by zams mass .
instead we employ the compactness parameter @xmath44 ( for a relevant mass scale @xmath45 , measured at the time of bounce ) .
as shown in @xcite , and further explored in @xcite , @xmath46 is a quantitative stellar structure parameter that describes the postbounce accretion evolution to a remnant mass scale @xmath45 .
we demonstrate that the preexplosion neutrino emission is very well parameterized by the compactness .
the preexplosion luminosities and mean energies of all neutrino species increase essentially monotonically with increasing @xmath46 .
we compute predicted integrated ibd interactions for a galactic core - collapse supernova in the super - kamiokande detector and show that the clear systematics governed by @xmath46 carries over to observation , even when standard msw neutrino oscillations are taken into account .
our results thus indicate that  in the absence of complicated collective neutrino oscillations  a high - statistics detection of neutrinos from the preexplosion phase will allow , in principle , a tight constraint of the compactness of the progenitor star s core .
this , however , will require knowledge of the nuclear eos and of the rotation rate of the collapsed core , since , as we show , both can dilute the otherwise clear compactness - dependent neutrino emission systematics .
this paper is structured as follows . in
[ sec : methods ] , we discuss our general relativistic hydrodynamics code gr1d and introduce its extension to neutrino radiation - hydrodynamics in the two - moment approximation , nugr1d . the initial models and the employed eos are discussed in  [ sec : models ] . in
[ sec : boltzmanncompare ] , we present results from a benchmark collapse and postbounce simulation that allows us to compare with the previously published code comparison of @xcite to assess nugr1d s ability to reproduce results of full boltzmann neutrino transport .
we present the results of our simulations in  [ sec : results ] , analyze the dependence of the neutrino signal on progenitor compactness , discuss predicted ibd signals from a galactic core collapse event in the super - kamiokande detector , and explore potential degeneracies introduced by eos and rotation .
finally , in  [ sec : discussion ] , we critically summarize our work and conclude by contrasting our results with the early neutrino signal observed from sn 1987a .
we make use of the open - source 1d general relativistic hydrodynamics code ` gr1d ` ( @xcite ; available at http://www.stellarcollapse.org ) outfitted with an energy - dependent multi - species m1 neutrino transport scheme in which the zeroth and first moments of the neutrino distribution function are evolved .
we refer the reader to @xcite for details on ` gr1d ` and describe in the following our current implementation of the transport scheme . for this first application
, we neglect the computationally - expensive energy - coupling neutrino interactions and transport terms  these terms are undoubtedly important for making highly accurate predictions of the neutrino signature ( see , e.g. , @xcite ) , but are unlikely to affect the general trends we observe .
we will address them in future work , but provide a discussion on the consequences of neglecting these terms via a comparison to full boltzmann neutrino transport simulations in .
our m1 scheme closely follows @xcite , who formulate the m1 evolution equations in a closed covariant form .
the scheme is simplified greatly by neglecting the energy - coupling terms .
this further requires that the velocity dependent terms are also ignored . in this limit , and using the schwarzschild - like metric and radial - polar slicing of ` gr1d ` and setting @xmath47 , the coordinate frame evolution equations for the neutrino energy density , @xmath48 , and the neutrino flux vector , @xmath49 , simplify from equations  3.37 and 3.38 of @xcite to @xmath50 and @xmath51 where @xmath52 is the neutrino interaction source term ( see below ) , @xmath53 is the lapse function and @xmath54 .
@xmath55 is the neutrino pressure tensor and is taken to be an interpolation between the two limiting cases of free streaming and diffusion .
we follow @xcite , who express @xmath56 as @xmath57 where @xmath58 is the eddington factor , taken here to be the maximum entropy closure in a closed , analytic form @xcite , @xmath59 in the no - velocity limit for ` gr1d ` , @xmath60 .
the free streaming and diffusion limits of the neutrino pressure tensor are @xmath61 and @xmath62 , respectively .
we set out to solve the system of equations via standard hyperbolic methods borrowed from conservative hydrodynamic schemes @xcite .
transport variables live at cell centers and we employ piece - wise linear reconstruction to cell interfaces with van leer s limiter @xcite and the hlle approximate riemann solver @xcite for calculating the intercell fluxes .
a complication arises when solving the neutrino moment equations in the high - opacity limit . in this case
, the standard fluxes returned from the riemann solver are dominated by a numerical diffusion term .
we follow @xcite and modify the fluxes to correct for this .
essentially , the modified fluxes return the diffusion - limit flux in the high opacity limit .
the modified form of the neutrino energy density flux through the @xmath63 interface is given by @xmath64 the corresponding modified interface flux for the neutrino flux evolution equation is , @xmath65 with @xmath66    in these equations , @xmath67 controls the modification to the fluxes to account for the high opacity .
following @xcite , we take @xmath68 where @xmath69 is the sum of the scattering and absorptive opacities .
these opacities are strong functions of energy and are also species dependent .
we note that when @xmath67 is 1 , the intercell fluxes reduce to the standard hlle approximation .
the characteristic speeds needed for the hlle scheme are calculated in the same spirit as the neutrino pressure tensor @xcite , @xmath70 where in the zero velocity limit , @xmath71 and @xmath72 .
@xmath73 and @xmath74 are the maximum and minimum values , respectively , of these characteristic speeds evaluated from both the right and left reconstructed variables .    finally , the source terms in and [ eq : fluxevolution ] are taken from @xcite . in the zero velocity limit , @xmath75 where @xmath76 , @xmath77 , and @xmath78 are the neutrino emissivity , neutrino absorption opacity , and the neutrino scattering opacity , respectively
we precompute the neutrino interaction terms for each neutrinos species ( we treat @xmath12 , @xmath13 and @xmath79 ) and neutrino energy group in dense tabular form as a function of density @xmath80 , temperature @xmath81 , and electron fraction @xmath82 .
we then use linear interpolation for efficient on - the - fly interpolation .
we include all standard iso - energetic scattering processes , charged - current absorption and emission , and thermal pair - production processes @xcite in the calculation of the neutrino interaction terms . since the neutrino  matter interactions for heavy - lepton neutrinos and antineutrinos are slightly different , ` nulib ` averages the two values of the emissivities and opacities .
our library of neutrino interaction routines , which we call ` nulib ` , is open source and available as a github repository at http://www.nulib.org . `
nulib ` requires an eos for the evaluation of the emissivities and opacities .
our treatment of thermal pair - processes in gr1d warrants some comments . since we do not currently consider energy ( or species )
coupling for thermal emission processes such as electron
positron annihilation to a neutrino
antineutrino pair , we compute an emissivity based on the thermal content of the matter ignoring any final state neutrino blocking . to limit the neutrino energy density to the equilibrium value ( where neutrino - antineutrino annihilation rates are in equilibrium with the thermal pair production rates ) , we use kirchhoff s law to derive an _ effective _ absorption opacity for neutrino  antineutrino annihilation from the thermal emissivity , @xmath83 where @xmath84 is the thermal energy density of neutrinos with energy @xmath48 and @xmath85}+1)$ ] is the equilibrium neutrino distribution function with chemical potential @xmath86 . as we shall see
, this method performs well at predicting the thermal neutrino flux of the heavy - lepton neutrinos during the preexplosion phase .    in ` nugr1d ` , we first update the hydrodynamic variables to the @xmath87-th timestep .
we then compute the neutrino opacities and emissivities associated with the updated hydrodynamic variables .
we update the radiation field operator - split .
the flux term is solved explicitly , using the radiation moments of the @xmath88-th timestep .
we calculate the neutrino
matter interaction terms using the @xmath87 radiation moments via a local implicit update . with the @xmath87 radiation energy density source term
, we then update the energy density and electron fraction of the matter .
we use 24 energy groups , with lowest - energy group centers at 0.5mev and 1.5mev , and then spaced logarithmically up to 200mev for @xmath12 , @xmath13 and @xmath20 .
we note that for the highest energy bins it occasionally occurs that the evolved neutrino flux vector exceeds the evolved neutrino energy density .
this tends to occur in the most dynamic phases of our simulations and where the opacities vary significantly from one zone to the next .
when this is the case we limit the neutrino flux to the neutrino energy density .
we extract the radiation quantities in the coordinate frame at a radius of 500 km .
we employ the most recent non - rotating solar - metallicity single - star model set from the stellar evolution code ` kepler ` @xcite .
this model set contains the presupernova configuration of 32 stars ranging in zams mass from 12@xmath34 to 120@xmath34 .
we denote individual models by @xmath89wh07 , where @xmath90 corresponds to the integer zams mass of the model , e.g. , @xmath91wh07 is the @xmath92-@xmath34 model of this model set . in @xcite
, we investigated this and other model sets in the context of black hole formation . under the assumption of a failed core - collapse supernova
, we found a strong empirical relation between the properties of the presupernova structure and the evolution of the failing supernova , e.g. , the time to black hole formation .
this led to a clear prediction : if we observe black hole formation in a failed core - collapse supernova via neutrinos , the lifetime of the protoneutron star ( and thus of the neutrino signal ) relays direct information about the presupernova structure .
however , such a prediction , ( _ i _ ) requires a failed supernova , which may not be the norm , and ( _ ii _ ) has a strong dependence on the nuclear equation of state . the empirical parameter introduced in @xcite is the compactness of the progenitor , measured at the time of core bounce .
it is an inverse measure of the radial extent of a given mass coordinate at the time of bounce , @xmath93 where @xmath94 is the radial coordinate that encloses a baryonic mass of @xmath45 at the time of core bounce . in @xcite ,
we chose @xmath95 , since this is the relevant mass scale for black hole formation , i.e. , a typical maximum baryonic mass at which a range of eos can no longer support a neutron star against gravity . in this study , we primarily use @xmath96 .
the motivation for this is that during the postbounce preexplosion phase , the relevant mass scale , especially for models with relatively small compactness , is much less than 2.5@xmath34 . in this study , we choose 1.75@xmath34 because this is close to the average baryonic mass inside the shock at 200300ms after bounce for all models : in the two extreme models that span the space in compactness parameter ( model s12wh07 , [ @xmath97 and @xmath98 , on the lower end ; model s40wh07 [ @xmath99 and @xmath100 on the upper end ) , the baryonic mass accreted through the shock at 250 ms after bounce is 1.45 m and 2.05 m , respectively .
we further justify our motivation of using @xmath96 over @xmath101 in . in , we plot both @xmath96 and @xmath101 versus zams mass for all 32 considered models .
@xmath96 is provided in for all models .    for
, one notes that while @xmath96 and @xmath101 differ quantitatively , there is no significant qualitative difference between them .
the overall trends transcending individual models remain , including the two regions of high compactness near 2225@xmath34 and 3545@xmath34 .
@xmath96 simply provides a more fine - grained parameterization at the lower mass scale relevant in the first few hundred milliseconds after bounce .
note , however , that there are a few models that have similar @xmath101 , but rather different density structure at small enclosed masses and radii and , hence , a different @xmath96 .
models @xmath10214wh07 and @xmath10216wh07 are examples .     and @xmath101 .
the mapping between zams mass and precollapse structure is highly non - monotonic , making the former an ill - suited parameter for describing progenitor structure in core collapse simulations . ]    in this study we perform core collapse simulations with each progenitor and two eos .
we use the eos of @xcite with a nuclear incompressibility of 220mev .
the ls220 eos is based on a compressible liquid - drop model of the nucleus . of the publicly available nuclear eos , the ls220 eos
best matches the constraints from nuclear theory and astrophysical observations ( see fig .  1 of @xcite and @xcite ) .
we also employ the relativistic mean field eos of @xcite that is based on the tm1 parameter set .
it is very different from the ls220 eos .
the maximum neutrino - less @xmath103-equilibrium cold neutron star gravitational masses are 2.04@xmath34 and 2.24@xmath34 for the ls220 and hshen eos , respectively .
the radius of a neutrino - less @xmath103-equilibrium cold neutron star with a gravitational mass of 1.4@xmath34 using the ls220 eos is 12.7 km . for the hshen eos ,
the corresponding radius is 14.6 km . for details on our particular implementation and the treatment of the low - density eos
, we refer the reader to @xcite .
the eos tables , reader and interpolation routines are available from http://www.stellarcollapse.org .
since our implementation of neutrino transport is new and approximate , a comparison with published results of full boltzmann neutrino transport is warranted .
this will allow us to test the ability of our code to reproduce the neutrino luminosities and spectral properties in the preexplosion phase .
we compare ` nugr1d ` with the results of @xcite , a comparison study between two boltzmann neutrino transport codes ] .
the two codes , @xcite , and ` vertex ` @xcite , approach the neutrino transport problem in very different ways .
their results compare very well in the newtonian limit , but show significant quantitative differences in the general relativistic case .
subsequent modifications to the approximate general relativistic potential used in ` vertex ` @xcite have since removed many of the quantitative differences between the codes .
the general relativistic test case of @xcite was the collapse and early postbounce evolution of a @xmath39-@xmath34 ( at zams ) solar - metallicity progenitor of @xcite , referred to as model @xmath10215ww95 in the following .
they employed the ls180 eos @xcite and a baseline set of neutrino - matter interactions , including coupling of energy groups via inelastic scattering processes @xcite .
we repeat their test here , using the same initial conditions and eos , with our current approximations and compare the neutrino observables .
we emphasize again that the current version of ` nugr1d ` lacks inelastic neutrino
electron scattering and velocity dependent transport terms . both
are included in the simulations of @xcite .
our transport scheme evolves only the zeroth and first moment of the neutrino distribution function , using an analytic closure to truncate the series of moment equations , whereas @xcite solve the full boltzmann equation for neutrino transport .
15ww95 progenitor .
these luminosities ( top panels ) and energies ( bottom panels ) correspond to the comparison study of @xcite .
the left panels contain results for @xmath12 , the right panels show @xmath13 ( thick lines ) and @xmath20 ( thin lines ) results .
the inset plot in the upper left panel shows the @xmath12 luminosity around core bounce .
shown in solid black lines are luminosities and root mean squared energies obtained with ` nugr1d ` .
the blue dashed lines and red dashed - dotted lines are the results from @xcite using the ` vertex ` code @xcite and code @xcite , respectively .
the thin dashed black lines ( labeled as ` nugr1d ` * ) are curves from a test simulation using profile data from the ` vertex ` simulation 5ms before bounce as starting data .
a detailed discussion of the differences is provided in the text . ]    in , we show the luminosities ( top panels ) and root mean squared energies ( bottom panels ) of three neutrino species : @xmath12 ( left panels ) , and @xmath13 and @xmath20 ( right panels ) .
both the luminosity and root mean squared energies are defined in @xcite , section  4 .
the black solid lines are the results obtained with ` nugr1d ` , the red dashed - dotted lines are the predictions of , and the blue dashed lines are the ` vertex ` results .
overall , the agreement is good , however there are are several systematic differences :    ( _ i _ ) the magnitudes of the @xmath12 and @xmath13 luminosities in the early postbounce phase predicted by ` nugr1d ` agree well with the results but they are systematically lower than the ` vertex ` results .
this discrepancy , which also exists between ` vertex ` and the , has been further investigated in @xcite and has since been resolved .
the updated ` vertex ` code employs an improved general relativistic potential and gives comparable amplitudes to , and hence ` nugr1d ` .
( _ ii _ ) the time when the silicon  oxygen interface ( located at a baryonic mass coordinate of 1.43@xmath34 in model @xmath10215ww95 ) accretes through the shock , which is marked by a sharp drop in the @xmath12 and @xmath13 luminosities , is earlier ( at @xmath104ms ) in our simulations than in the simulations of @xcite ( @xmath105ms in the ` vertex ` simulations and @xmath13 luminosities begins around the same time as in the ` vertex ` simulations . ] ) .
we have carried out a number of tests to try to understand this difference .
we have varied the zero point of the internal energy , which affects the relativistic enthalpy entering the momentum and energy equations and alters the collapse time .
we have attempted different mappings of the presupernova stellar structure to our eulerian grid , we have tested variations in precollapse electron capture by parameterizing @xmath82 as a function of density in the collapse phase @xcite , and we have replaced the low density compositions with a pure silicon gas ( which replaces the nse equivalent of mainly iron and therefore better represents the low density matter ) .
none of these tests led to a change of the postbounce time to silicon
oxygen interface accretion by more than @xmath510ms .    without currently having the modeling technology to test it
, we suspect that the apparent difference may be due to the only other obvious difference between ` nugr1d ` and ` vertex ` : our lack of detailed neutrino physics during the collapse phase .
two facts lead us to this conclusion .
the collapse time discrepancy between the ` vertex ` ( @xmath106ms ) and ( @xmath107ms ) and ` nugr1d ` ( @xmath108ms ) can be explained by the lack of inelastic scattering .
@xcite gives a @xmath109ms in when neutrino - electron inelastic scattering is neglected .
this collapse time difference , @xmath540ms , is the difference in the silicon - oxygen interface accretion times .
also , when we map in the ` vertex ` profiles from @xcite at 5ms before bounce to ` nugr1d ` and continue the evolution we reproduce the postbounce time that the silicon - oxygen interface reaches the shock .
we show the result of this simulation with a thin dashed curve , marked as nugr1d * , in .
it is worth noting that there may be additional differences due to the simplistic treatment of the eos at low densities .
below a density of @xmath110g@xmath111 , ` vertex ` replaces the nuclear statistical equilibrium ( nse ) eos with an eos that specifically depends on the composition of the matter . assumes all matter below this density is silicon .
both have some treatment of nuclear burning . in ` gr1d ` , as discussed in @xcite , we assume nse compositions from the nuclear eos . at densities below the validity regime of the nuclear eos , we take the compositions at the lowest density point from the nuclear eos and use the timmes eos @xcite . for reference ,
the initial density of the silicon
oxygen interface is 0.41.0@xmath112g@xmath111 , where the range represents the extent in density space .
( _ iii _ ) the @xmath12 and @xmath13 root mean squared energies , predicted by ` nugr1d ` agree very well with the boltzmann transport results during the postbounce phase .
the difference seen in the @xmath20 root mean squared energy is similar to that observed by @xcite and @xcite when investigating the effects of inelastic neutrino  electron scattering . in the postbounce evolution , this interaction is expected to predominately affect the @xmath20 neutrino .
we currently ignore this process in ` nugr1d ` and note that the @xmath20 luminosity predicted by ` nugr1d ` still agrees well with the full boltzmann results .
( _ iv _ ) another difference between the evolution in ` nugr1d ` and the full boltzmann transport results arises in the collapse phase .
the lack of velocity terms and inelastic @xmath113 scattering significantly effects the composition of the inner core . in simulations with inelastic @xmath113 scattering ,
neutrinos from electron capture on free protons down - scatter off of electrons to lower energies .
since the optical depth is lower , these neutrinos can then escape , deleptonizing the core .
in our simulations , these high - energy neutrinos can not down - scatter and therefore can not escape .
deleptonization is suppressed until later phases , when central density and temperature are higher .
the lack of velocity - dependent terms delays full trapping by neutrino advection that would normally begin to occur at @xmath114g@xmath111 until nuclear densities , allowing for further deleptonization . at
bounce the central value of @xmath82 are @xmath50.22 compared to @xmath50.29 in ` vertex ` and .
the root mean squared @xmath12 energies predicted by ` nugr1d ` are higher during the prebounce phase , because the neutrinos do not experience the down - scattering via inelastic neutrino  electron scattering .
overall , we find that differences in the inner core structure and composition at bounce do not strongly present themselves in the neutrino signal after the collapse phase .
therefore , for this study , we find the current version of ` nugr1d ` to be acceptable , since our primary focus is the neutrino signal of the preexplosion accretion phase .
we do note , however , that the lack of energy - coupling terms in our transport can cause qualitative differences near black hole formation .
when energy - coupling terms , such as gravitational redshift , are included , then the @xmath12 and @xmath13 luminosities drop off at times very close to black hole formation , as seen in @xcite .
this effect can be captured with ` gr1d ` s leakage scheme in which it is trivial to include redshift terms @xcite .
capturing redshift in an energy - dependent transport scheme requires energy - group coupling .
we have also compared our m1 scheme to the results of @xcite for the 40@xmath34 model from @xcite using the ls180 eos . in this model a black hole forms within 500ms of bounce .
we find differences in the @xmath12 and @xmath13 luminosities of @xmath51020% in the last @xmath550ms , but good agreement in the early postbounce phase .
rcccccccccccc @xmath10212 & 0.235 & 19.24 / 14.19 / 7.73 & 1.02 / 0.92 & 2.13 / 1.78 + @xmath10213 & 0.383 & 22.18 / 16.58 / 8.76 & 1.25 / 1.09 & 2.53 / 2.07 + @xmath10214 & 0.537 & 25.19 / 19.35 / 9.24 & 1.36 / 1.20 & 3.06 / 2.49 + @xmath10215 & 0.580 & 25.51 / 19.59 / 9.17 & 1.30 / 1.16 & 3.13 / 2.59 + @xmath10216 & 0.338 & 18.91 / 13.72 / 8.20 & 1.11 / 0.95 & 2.00 / 1.68 + @xmath10217 & 0.383 & 19.93 / 14.54 / 8.57 & 1.20 / 1.02 & 2.13 / 1.78 + @xmath10218 & 0.738 & 28.66 / 22.26 / 10.30 & 1.55 / 1.36 & 3.62 / 2.92 + @xmath10219 & 0.544 & 23.36 / 17.49 / 9.37 & 1.43 / 1.23 & 2.67 / 2.18 + @xmath10220 & 0.944 & 29.39 / 22.70 / 11.08 & 1.79 / 1.55 & 3.64 / 2.91 + @xmath10221 & 0.325 & 18.48 / 13.41 / 8.09 & 1.04 / 0.91 & 1.95 / 1.65 + @xmath10222 & 0.972 & 30.06 / 23.29 / 11.26 & 1.82 / 1.58 & 3.76 / 3.00 + @xmath10223 & 1.256 & 42.00 / 33.88 / 14.68 & 2.29 / 2.00 & 5.99 / 4.67 + @xmath10224 & 1.167 & 39.14 / 31.35 / 13.66 & 2.16 / 1.89 & 5.44 / 4.24 + @xmath10225 & 1.040 & 32.30 / 25.28 / 11.87 & 1.93 / 1.68 & 4.15 / 3.28 + @xmath10226 & 0.727 & 24.97 / 18.84 / 10.10 & 1.59 / 1.38 & 2.88 / 2.33 + @xmath10227 & 0.783 & 25.67 / 19.45 / 10.21 & 1.61 / 1.40 & 3.01 / 2.42 + @xmath10228 & 0.640 & 25.86 / 19.83 / 9.54 & 1.36 / 1.21 & 3.16 / 2.62 + @xmath10229 & 0.556 & 23.24 / 17.38 / 9.43 & 1.45 / 1.22 & 2.64 / 2.17 + @xmath10230 & 0.760 & 26.61 / 20.30 / 10.19 & 1.60 / 1.40 & 3.20 / 2.58 + @xmath10231 & 0.687 & 25.23 / 19.09 / 9.95 & 1.56 / 1.35 & 2.96 / 2.40 + @xmath10232 & 0.883 & 28.10 / 21.58 / 10.70 & 1.71 / 1.49 & 3.43 / 2.75 + @xmath10233 & 0.965 & 30.36 / 23.57 / 11.28 & 1.82 / 1.58 & 3.82 / 3.04 + @xmath10235 & 1.129 & 35.93 / 28.50 / 12.77 & 2.06 / 1.79 & 4.83 / 3.79 + @xmath10240 & 1.328 & 48.32 / 39.45 / 17.51 & 2.49 / 2.21 & 7.23 / 5.75 + @xmath10245 & 1.300 & 46.42 / 37.81 / 16.56 & 2.42 / 2.14 & 6.86 / 5.44 + @xmath10250 & 0.701 & 25.60 / 19.42 / 10.00 & 1.57 / 1.32 & 3.02 / 2.45 + @xmath10255 & 0.577 & 23.56 / 17.65 / 9.52 & 1.48 / 1.25 & 2.69 / 2.20 + @xmath10260 & 0.461 & 21.52 / 15.88 / 9.07 & 1.33 / 1.12 & 2.36 / 1.95 + @xmath10270 & 0.755 & 25.83 / 19.58 / 10.17 & 1.60 / 1.39 & 3.04 / 2.46 + @xmath10280 & 0.591 & 23.75 / 17.80 / 9.63 & 1.48 / 1.26 & 2.71 / 2.22 + @xmath102100 & 0.792 & 29.90 / 23.33 / 10.59 & 1.61 / 1.42 & 3.82 / 3.07 + @xmath102120 & 0.474 & 22.31
/ 16.64 / 9.02 & 1.33 / 1.15 & 2.52 / 2.07 +    [ tab : results ]      we perform core collapse and early postbounce evolutions of all 32 models introduced in using both the ls220 and the hshen eos . in we
present the three neutrino luminosities and average energies for each model and eos as a function of postbounce time .
we do not expect a clear trend in the neutrino observables with zams mass .
however , we do expect trends based on the presupernova structure of the star , which is well encapsulated by the compactness parameter introduced in  [ sec : models ] . in the top set of panels
we show simulations run using the ls220 eos . in the bottom panels ,
simulations performed with the hshen eos are shown . to highlight that there is indeed a trend with presupernova structure , we color - code individual models according to their compactness parameter .
the mapping between line color and @xmath96 is provided on the right . to more directly highlight the eos dependence
, we include the luminosity and average energies of two models run with the ls220 eos in the hshen eos panels with thick dashed lines .
these models , @xmath10212wh07 and @xmath10240wh07 , have the lowest and highest compactness parameter in our model set , respectively .
we find that there is little variation in the peak luminosity of the @xmath12 neutronization burst signal .
for all 32 models simulated using the ls220 ( hshen ) eos , the peak amplitude varies by less than 3% ( 5% ) from the average .
this reflects the universal nature of the collapse of the inner core @xcite .
after the neutronization burst , the postbounce luminosities of all species increase systematically with increasing compactness parameter .
models with higher @xmath96 have higher temperatures throughout the protoneutron star @xcite .
this increases the diffusive neutrino luminosity and is best seen in the @xmath20 luminosities .
the postbounce accretion rate also increases with the compactness parameter .
higher accretion rates , and the deeper gravitational potential due to the higher protoneutron star mass , increase the accretion luminosity , which is most directly reflected in the @xmath12 and @xmath13 signals . after roughly 100ms , the average energy of the emitted neutrinos also shows an increasing trend with the compactness parameter .
the matter temperature at the neutrinosphere is higher in models with larger compactness , therefore a higher average neutrino energy is observed at infinity .
the neutrino luminosities and average energies from simulations using the hshen eos are systematically lower than the luminosities and average energies from simulations of the same model run with the ls220 eos .
this is clearly seen in the bottom set of panels in : models @xmath10212wh07-ls220 and @xmath10240wh07-ls220 have luminosities and average energies that are comparable to or larger than in the corresponding hshen models . for a fixed accretion rate ( or fixed progenitor model ) , the location of the neutrinosphere of each species influences the emitted luminosities and spectra . in models evolved with the stiff hshen eos ,
the neutrinospheres are located systematically at larger radii and lower matter temperatures than in models run with the softer ls220 eos .
for example , in the @xmath10212wh07 simulations , the rosseland - mean @xmath12 neutrinospheres at 200ms after bounce have radii and temperatures of @xmath535.3 km and @xmath54.86mev ; and @xmath539.5 km and @xmath54.46mev for the ls220 and the hshen eos , respectively .
the larger neutrinosphere radii are responsible for the lower accretion luminosity since the latter is set essentially by the product of the mass accretion rate and the gravitational potential at the protoneutron star surface @xcite .
the latter is located at larger radii in simulations using the hshen eos .
the differences in the neutrinosphere radii and temperatures between the ls220 and hshen eos also give an explanation for the systematically lower average neutrino energies seen in the hshen simulations .
matter at larger radii has been compressed less and therefore is cooler .
this leads to average neutrino energies that can be up to 5mev lower for the hshen eos than for the ls220 eos for the same progenitor model ( see ) .
the difference in the neutrino luminosities and average energies between the two eos is largest for models with @xmath115 , in which the high accretion rates lead to the accumulation of @xmath52@xmath34 of material inside the shock within @xmath5200300ms of bounce . in the case of the ls220 eos
, this leads to very high temperatures throughout the protoneutron star as it becomes more and more compact and closer to gravitational collapse to a black hole . in our simulations , the most compact model @xmath10240wh07 ( @xmath99 , @xmath116 ) forms a black hole 503ms after bounce .
the slightly less compact model @xmath10245wh07 ( @xmath117 , @xmath118 ) forms a black hole 563ms after bounce .
the high temperatures present in the ls220 simulations at these times will not be obtained until postbounce times @xmath119 in models using the hshen eos40wh07 and @xmath10245wh07 models are @xmath51.3s and @xmath51.4s , respectively , when using the hshen eos . ] .
+     most neutrino detectors are most sensitive to the electron antineutrino luminosity through the dominant ibd interaction , @xmath120 . in the left panel of
, we consider the cumulative emitted @xmath13 energy for each model using the ls220 eos .
we color code the models based on their compactness parameter and include two reference models that use the hshen eos , model @xmath10212wh07 and model @xmath10240wh07 .
the graphs shown in this panel are the integral of the graphs shown in the top center panel of .
it is obvious that the cumulative amount of @xmath13 energy emitted during the preexplosion phase strongly correlates with the compactness of the progenitor model .
for example , the amount of emitted @xmath13 energy from model @xmath10240wh07 ( @xmath99 ) is _ always _ between two and three times of that of model @xmath10212wh07 ( @xmath97 ) .
we make this point more quantitative in the center and right panels of .
in the center ( right ) panel we plot the cumulative emitted @xmath13 energy at 100 , 200 , 300 , and 400ms after bounce for both eos as a function of @xmath96 ( @xmath101 ) . for reference
, we present a subset of these numbers in .
we see a very clear correlation that depends only weakly on the chosen eos .
note , however , that for models with small compactness parameter ( @xmath121 ) , the correlation between the total emitted @xmath13 energy and the compactness parameter is not as strong after 400ms of postbounce evolution .
comparing the center and right panels justifies our choice of @xmath96 over @xmath101 as explained in .
the onset of an explosion will break the correlation observed in .
once it is launched , the accretion luminosity effectively turns off and only the diffusion luminosity remains .
one also expects this diffusion luminosity to show a correlation with the compactness of the progenitor , since the remnant protoneutron star s thermodynamic conditions , such as the central entropy and its mass , are essentially set by the presupernova structure .
however , it is currently unclear whether one should obtain a correlation between explosion time and the compactness parameter .
clarification will require a more complete understanding of the core - collapse supernova explosion mechanism and may require extensive parameter studies with fully self - consistent three - dimensional simulations .       and is provided in .
the dashed lines are results for models @xmath10212wh07 and @xmath10240wh07 run with the hshen eos . in the inset
we show the cumulative ibd interactions as a function of @xmath96 for each model and eos at four postbounce times : 100 , 200 , 300 , and 400ms . ]
we use the publicly available software ` snowglobes ` ] to predict the neutrino signal observed in earth - based neutrino detectors . `
snowglobes ` @xcite , which in turn relies on ` globes ` @xcite , is a set of routines that compute the interaction rates of supernova neutrinos in user - specified detector configurations .
variables include detector material ( e.g. , water , scintillator , argon , and lead ) , detector volume , detector response functions , and a host of relevant neutrino interactions .
for this investigation , we consider only ibd interactions in a water cherenkov detector .
we choose a detector mass of 32kt , the mass of water in super - kamiokande that is sensitive to core - collapse supernova neutrinos @xcite . for reference
, we use the ` wc100kt30prct ` smearing rates and efficiencies provided with ` snowglobes ` .
we construct ` snowglobes ` initial fluence data from our simulations binned in 5ms intervals .
we provide these energy - dependent fluences at http://www.stellarcollapse.org/m1prog for all models , neutrino species , and both eos in 5ms intervals up to 450ms after bounce .
we assume a fiducial galactic supernova distance of 10kpc . in , we show the cumulative number of interactions for each model run with the ls220 eos and , for reference , two models run with the hshen eos .
the lines are color - coded according to compactness parameter .
note that the vertical scale in this figure is in thousands of interactions .
we note that our answers agree with the total number of expected detected interactions in super - kamiokande from a galactic core - collapse supernova at 10kpc , which is estimated to be @xmath57000 @xcite .
to arrive at this number from our results , consider the lowest zams mass progenitor in our model set , model @xmath10212wh07 .
after 450ms of evolution , 15b of @xmath13 energy has been radiated ( ) , which corresponds to 2000 ibd detected interactions ( ) . for 50b of released energy ( @xmath51/6 of 300b , the fiducial energy released in neutrinos over the entire cooling phase ) , one would then expect @xmath57000 detected interactions .
however , as is clear from , the number of detected interactions from the next galactic supernova may be higher than this fiducial number .
more importantly , the rate of interactions in the preexplosion phase will give us detailed information on the progenitor core structure .    in order to more directly quantify the differences between variations in progenitor compactness and variations in the nuclear eos , we plot in the inset of the number of expected ibd detected interactions in a super - kamiokande - like water cherenkov detector at various postbounce times versus @xmath96 .
there is a well defined trend : the number of ibd interactions detected in the first 100 , 200 , 300 , and 400ms increases with the compactness parameter of the models . for reference , we include the expected number of interactions at 200 and 400ms for both eos in .
we find that the eos dependence of the expected number of interactions is similar to the eos dependence of the total emitted @xmath13 energy : the hshen eos leads to a lower number of interactions ( compare the inset of to the center panel of ) .
the dependence on eos is somewhat stronger here , since the the lower average @xmath13 energy predicted from stiffer eos translates into a reduced cross section in earth - based detectors .
in addition to the total number of interactions , a water cherenkov detector measures individual energies , and thus , allows for the reconstruction of the cumulative emitted @xmath13 energy over time .
this reconstruction will depend on the detector s response function and efficiency .
an additional independent path to experimentally probing the inner structure of the progenitor is via the total neutrino energy emitted in all species over the first 10s of seconds after the initial collapse .
this method requires a measurement of the total fluence of neutrinos of all species , not just electron antineutrinos .
examples of neutrino interactions capable of relaying such information are the mono - energetic de - excitation of a neutral - current neutrino - excitation of @xmath122c @xcite or neutrino - proton elastic scattering interactions @xcite .
such measurements would require good energy resolution , a significant source of carbon and/or a low energy threshold , for example , a liquid scintillator neutrino detector .
we note that even with a liquid scintillator detector , the dominant neutrino interaction is still ibd @xcite .
if such a measurement was made , and there is not a significant amount of rotation ( see the discussion on rotation in ) , one can immediately infer the gravitational binding energy of the remnant , since neutrinos carry away the vast majority ( @xmath1231% of the energy is predominantly shared among the kinetic energy of the explosion , the original binding energy of the unbound stellar mantle , and the binding energy of the iron core at the onset of collapse . ] ) of the gravitational binding energy . for typical nuclear eos like the ones considered here ,
this results in a one - to - one mapping of the released gravitational binding energy to the baryonic mass of the remnant , and , hence , the gravitational mass of the remnant .
this is most easily seen by fitting the gravitational binding energy of a cold ( t=0.1mev ) , neutrino - less @xmath103-equilibrium , non - rotating neutron star to its baryonic mass . from cold neutron star tov solutions using the ls220 eos
one can obtain an empirical fit to better than 3% above a baryonic mass of 1.15@xmath34 , @xmath124 a similar fit for the hshen eos gives , @xmath125 and is accurate to 5% above baryonic masses of 1.15@xmath34 .
below @xmath126 , the empirical quadratic fit is not as accurate .
however , all models considered here reach a baryonic protoneutron star mass of 1.15@xmath34 within @xmath510ms of bounce . hence , we believe that the above fits are acceptable for the iron - core core collapse events considered here .
we now make the assumption that an explosion launched at a particular postbounce time will result in a neutron star remnant with a baryonic mass equal to the baryonic mass that has accreted through the shock up until the time of the explosion .
this neglects any late - time fallback of material onto the protoneutron star which would lead to additional neutrino emission . @xcite and
@xcite predict fallback masses @xmath1275 - 10% of the initial protoneutron star remnant mass .
we also neglect any asymmetric mass accretion which may occur in the early explosion phase . in , we convert the baryonic mass enclosed by the shock to the _ total _ emitted neutrino energy using and [ eq : ebind_hshen ] .
we plot this for all progenitor models ( run with the ls220 eos ) as a function of the hypothetical time of explosion . as a concrete example , consider the situation where 300b ( shown as the dashed line in ) of total neutrino energy was observational inferred . this could correspond ( @xmath128 ) to a progenitor with a high compactness parameter that exploded at an early time , e.g. , model @xmath10240wh07 at 70ms or ( @xmath129 ) to a low progenitor with low compactness parameter that exploded at late times , e.g. , model @xmath10212wh07 at 400ms .
if we have an estimate of the explosion time , e.g. , via characteristic features in the neutrino observables , then we can use the combined measurement to probe the progenitor core structure .
this is further quantified in , where we choose three total emitted neutrino energies , 250 , 300 , and 350b , and determine the time at which the explosion must have been launched for a given compactness and the respective total emitted energy .
we plot this explosion time versus @xmath96 for all models and both eos . in general , for a fixed total emitted energy , as the compactness parameter of the progenitor increases , the explosion time must decrease . for progenitors with high @xmath96 ( @xmath130 )
, there is a clear mapping between the explosion time and the compactness parameter , given a specific total emitted neutrino energy and an eos . as was the case for the total ibd rates ,
there is some scatter at low @xmath96 ( @xmath131 ) , but there is still an overall trend .    , and the baryonic mass enclosed in the shock at any given time .
this defines the explosion time to be , in a lagrangian sense , the time at which the outermost final neutron star mass element accretes through the shock .
the color coding corresponds to @xmath96 , the color coding is provided in . ]
shows that there is a very strong eos dependence . while the baryonic mass inside the shock as a function of postbounce time does not vary strongly with eos
, the gravitational binding energy released does .
based on the above empirical fits to the gravitational binding energy , the ls220 eos leads to @xmath132 more energy release than the hshen eos for the same baryonic mass .
this method of determining the compactness parameter ( by combining an estimate of the explosion time and the total emitted neutrino energy ) can be used together with the methods described above ( using the postbounce , preexplosion ibd rates ) as a consistency check , or to break degeneracies , which we will discuss next in .
14% energy for a given baryonic mass . ]
there are a number of degeneracies and uncertainties that may prevent fully conclusive statements regarding the mapping from detected signal to progenitor core structure .
these include , nuclear eos , rotation , viewing angle , distance , and neutrino oscillations ( including collective oscillations ) .
we will discuss each one of these consecutively and independently , although all may be relevant in a generic situation .
_ nuclear eos _
: we have already briefly explored the dependence of the neutrino observables on the nuclear eos by taking two very different eos and comparing the emitted neutrino signal .
the total number of interactions predicted to be detected from a given progenitor model varies with eos ( cf .
, ) in such a way that a high compactness model paired with the hshen eos produces the same number of interactions as a model with a slightly lower compactness paired with a softer eos , like the ls220 . hence , there is a clear degeneracy between the progenitor compactness and the equation of state .
however , a water cherenkov detector will not only detect a given number of ibd interactions , but also their energy distribution . from this distribution , and the detector response , one can work backwards to reconstruct both the emitted average @xmath13 energy and the luminosity .
we show in how knowledge of the emitted energy spectrum can break the degeneracy with cumulative emitted @xmath13 energy and eos , provided the distance to the core collapse event is known ( see the discussion of distance uncertainties later in this section ) . for each model and eos ,
we plot the average @xmath13 energy as a function of cumulative emitted @xmath13 energy at three select postbounce times : 100 , 200 , and 300ms .
the data from each of the 32 progenitors fall on a unique line that is parameterized by the compactness . in principle , provided very large statistics and ultimate predictions of neutrino observables , this allows one to break the eos - progenitor degeneracy between eos . in reality , achieving an instantaneous measurement of the neutrino average energy to within a few % will be difficult , and there will likely be a large suite of nuclear eos that could reproduce the observables within error .
energy versus cumulative emitted @xmath13 energy each eos and progenitor .
this is shown for three postbounce times , 100ms ( left panel ) ; 200ms ( center panel ) ; and 300ms ( right panel ) .
solid ( open ) points correspond to the ls220 ( hshen ) eos .
the degeneracy between eos and cumulative @xmath13 energy is broken when the average energy is taken into account .
we find that is a robust prediction across all progenitors and postbounce preexplosion times . ]    _ rotation and viewing angle _ : the effect of rotation on the neutrino signal may be more difficult to disentangle from the effects of progenitor structure . to explore the effect of rotation ,
we perform 1.5d simulations of rotating core collapse .
these simulations treat rotation in a spherically symmetric way , using ` shellular rotation ' @xcite .
this approximation only captures the spherically - averaged centrifugal effect of rotation on the matter . as the amount of rotation
is increased , the protoneutron star becomes more and more centrifugally supported .
increased rotation leads to lower densities and temperatures throughout the protoneutron star .
this in turn effects the neutrino signal . in , we plot the @xmath13 luminosities and average @xmath13 energies determined from ` nugr1d ` for 12 rotating core collapse simulations .
we use model @xmath10215wh07 paired with the ls220 eos .
the initial rotation rate is assigned via @xcite @xmath133^{-1}\times 10^{16}\ \mathrm{cm}^2\ \mathrm{s}^{-1}\,,\label{eq : rotationlaw}\ ] ] where @xmath134 is the radius that encloses 1@xmath34 , which for the @xmath10215wh07 progenitor is 703 km , and @xmath135 is the specific angular momentum at infinity , in units of @xmath136 .
we vary @xmath135 from 0 to 3 in increments of 0.25 . for the @xmath10215wh07 progenitor ,
the initial central angular velocity is @xmath137rads@xmath138 , giving a range of @xmath139 from 0 to 6.08rads@xmath138 . at bounce ,
the @xmath140rads@xmath138 model has a rotation period of @xmath52.4ms , which decreases nearly linearly with increasing initial @xmath135 .
increased rotation leads to the overall reduction of both the @xmath13 luminosity and average energy .
for example , with @xmath139 = 0 , 2.03 , 4.05 , and 6.08rads@xmath138 the @xmath13 luminosity at 100ms is @xmath555 , @xmath552 , @xmath540 , and @xmath525bs@xmath138 , respectively .
while we do not show the @xmath12 and @xmath20 luminosities and average energies , they follow the @xmath13 trends . the only exception to this is the @xmath12 neutronization burst where rotation in our 1.5d approach does not significantly alter the @xmath12 luminosity .
our 1.5d approach to rotation does not capture the anisotropy of the neutrino emission , rather just the angle - averaged value .
@xcite considered a model with an initial rotation rate of @xmath141rads@xmath138 , roughly corresponding to our @xmath142 model .
they found that the ratio of polar to equatorial luminosity can be as large as 34 and that average neutrino energies in polar regions can be harder by @xmath14312mev than on the equator .
we find that it is not possible with the cumulative neutrino signals alone to break the degeneracy between rotation and the progenitor star
 eos combinations . the angle - averaged @xmath13 luminosity and average energy from a rotating system
mimic the @xmath13 luminosity and average energy from a non - rotating system with a progenitor with lower compactness and/or a stiffer eos .
in addition , the rotational energy stored in the neutron star remnant can lead to an underestimate of the total emitted neutrino energy , again mimicking a lower - compactness progenitor and/or a stiffer eos .
however , all hope is not lost : stellar evolution and current constraints on pulsar birth spins ( e.g. , @xcite ) suggest that most massive stars will be spinning too slowly for rotation to have a strong effect on dynamics and neutrino signal .
@xcite , for example , predict a presupernova rotation rate of 0.2rads@xmath138 for a 15@xmath34 star , which is smaller than all rotation rates considered here . if , on the other hand , the progenitor is rapidly spinning , the degeneracy in the neutrino signal may be broken by coincident observations in gravitational waves that are able to constrain the rotation rate of the collapse core @xcite .
15wh07 for various initial central angular velocities .
we assign angular momentum to the @xmath10215wh07 progenitor via with values of @xmath135 ranging from 0 to 3 in 0.25 increments .
this gives initial central rotation rates of 0 to @xmath56.08 rad s@xmath138 in increments of @xmath50.51 rad s@xmath138 . in our simulations , increased rotation lowers both the neutrino luminosity and the average energy . ]
_ distance _ : if the electromagnetic signal is blocked , by , for example , the galactic center , we may not be able to obtain a reliable distance to the supernova . the @xmath13 flux at earth follows an inverse square law , while the spectral distribution will not change over the distances of relevance .
this means that degeneracy with distance is also hard to break with the ibd neutrino signal alone for the following reason : a high - compactness progenitor at a large distance can produce the same energy - integrated flux at earth as a low - compactness progenitor .
less neutrinos are emitted in the latter but the flux has not been diluted as much due to the closer distance .
if the nuclear eos is not known , the observed neutrino flux and energy spectra could be associated with a range of nuclear eos  distance combinations .
a progenitor at a large distance could produce the same neutrino flux and energy distribution as a lower - compactness progenitor with a softer eos . here
, the softer eos , which will give rise to an increased average energy , compensates for the decrease in the average energy produced by the smaller compactness .    future constraints on the nuclear eos could break this degeneracy since one could more reliably associate an observed energy distribution with a particular progenitor compactness ( ) . this may even provide a distance estimate , as would a detection of the neutronization burst signal @xcite .
it does not show a strong dependence on the progenitor model , but is difficult to detect since the @xmath12 cross sections in water cherenkov detectors are much lower than the ibd cross section , resulting in few interactions .
some detector materials have significantly larger @xmath12 cross sections and would be better suited to detect the neutronization burst , such as liquid argon @xcite and lead @xcite .    _
neutrino oscillations : _ first , considering only matter - induced neutrino oscillations , in the normal mass hierarchy , the electron antineutrino signal at earth is a composite spectrum of @xmath14470% of the original @xmath13 neutrinos and @xmath14530% of the original @xmath146 spectrum , where @xmath147 is the mixing angle between the mass eigenstates @xmath148 and @xmath149 . in the inverted neutrino mass hierarchy , the entire @xmath13 signal is replaced with the original @xmath150 signal @xcite . in either case
, we still expect the total number of ibd interactions to increase with the compactness parameter of the progenitor since both the @xmath13 and the @xmath20 luminosity increase .
if the hierarchy of the neutrino mass eigenstates is not determined by neutrino experiments before the next nearby core - collapse supernova , the early postbounce , preexplosion neutrino luminosities may provide an answer @xcite .
this relies on the systematically different rise times of the @xmath20 and @xmath13 signals .
we also point out a degeneracy in disentangling the spectral properties of the emitted neutrino spectra from the matter - induced oscillated neutrino spectra observed in earth detectors @xcite .
much more troublesome are the collective neutrino oscillations that arise from coherent neutrino - neutrino forward scattering .
collective oscillations are very sensitive to the energy spectra ( both the distribution and magnitude ) of all neutrino flavors and the background matter density .
since the governing equations are highly non - linear , there are currently no simple analytic expressions predicting the neutrino signal at earth based on the output of core - collapse simulations .
recent studies suggest that during the early postbounce , preexplosion phase , collective neutrino oscillations may be suppressed @xcite , however the community has not yet reached consensus , see , e.g. , @xcite and @xcite .
the next nearby core - collapse supernova will be extremely well observed in neutrinos .
super - kamiokande alone will observe @xmath5@xmath151 electron antineutrinos from a typical core - collapse supernova at a fiducial galactic distance of @xmath11 .
future detectors of the scale of the proposed hyper - kamiokande may see in excess of @xmath152 interactions .
such high - statistics observations will provide rich information on the neutrino signal .
comparison with theoretical model predictions will allow to falsify or constrain a broad range of hypotheses in core - collapse supernova astrophysics and nuclear / neutrino physics .
unexpected signal features may lead to the discovery of new physics .    in this study , our focus has been on the imprint of the progenitor star s structure on the neutrino signal in the postbounce preexplosion phase of core - collapse supernovae .
we have carried out a large set of spherically - symmetric radiation - hydrodynamics simulations of core collapse and the early postbounce phase with the goal of studying trends in the neutrino signal with variations in progenitor structure .
our results show , in agreement with previous work ( e.g. , @xcite ) , that the @xmath12 signal from the neutronization burst emerging shortly after bounce has very little progenitor dependence , due to the universal nature of homologous inner core collapse .
the neutrino signal in the postbounce preexplosion phase is determined primarily by the accretion luminosity of outer iron core , silicon shell , and oxygen shell material .
the postbounce accretion rate depends on the inner structure of the progenitor at the presupernova stage .
our results show that the preexplosion neutrino signal has an essentially monotonic dependence on progenitor structure described by a single parameter , the compactness @xmath153 ( where @xmath45 is a typical baryonic mass reaching the center over the timescale of interest ) . the greater a progenitor s @xmath46 , the higher are the emitted luminosities and average energies of all neutrino species . scaling in the same way
is the total emitted energy in neutrinos over the entire protoneutron star cooling phase for a given explosion time at which accretion is shut off .
these trends are robust and independent of the nuclear eos .
they are also rather insensitive to the particular choice of the reference mass @xmath45 in @xmath46 as long as it is in the range of typical neutron star baryonic masses ( @xmath5@xmath154 ) and we find @xmath96 to be a good choice .     and @xmath20 spectra .
the color coded lines denote different progenitors and are generated with ` snowglobes ` based on our simulations with the ls220 eos .
more details are provided in the text .
we overlay the cumulative detected interactions observed in kamiokande
ii from sn 1987a , assuming the first interaction denotes the time of core - bounce .
we caution the reader by noting that the very low number statistics make it hard to draw any conclusions regarding the progenitor of sn1987a . ]
the monotonic dependence of the preexplosion neutrino emission on progenitor compactness translates directly to the neutrino signal observed by detectors , provided collective neutrino oscillations do not lead to complicated swaps of flavor spectra that brake the dependence of the observed signal on progenitor structure .
neutrino observations of the next nearby core collapse event thus may , in principle , allow quantitative constraints on the inner structure of the progenitor star .
as an example with real neutrino data , we consider the early postbounce neutrino signal observed from sn 1987a by the kamiokande
ii experiment @xcite .
of the eleven interactions that were observed , the first four occurred within 323ms of each other . all interactions observed by kamiokande
ii are consistent with being ibd interactions @xcite .
we assume that the first interaction occurs at the onset of the postbounce phase , although the actual bounce time is likely somewhat earlier . in , we plot the cumulative number of detected interactions observed from sn 1987a in the first 500ms along with the ` snowglobes ` prediction from our simulations with the ls220 eos using a 2.14kt water cherenkov detector at 51.4kpc .
we use the efficiencies quoted in @xcite and the smearing matrices from the ` snowglobes ` detector configuration ` wc100kt30prct ` . to show the full range of the possible effects of msw neutrino oscillations ,
we show the expected number of interactions assuming no oscillations and assuming a complete switch of the @xmath13 and @xmath20 spectra ( as would be the case in the inverted mass hierarchy ; the normal hierarchy is a combination of these two signals ) . during the accretion
phase the expected number of detected interactions from an oscillated @xmath20 spectrum can be significantly smaller than the electron - type neutrino luminosity due to the smaller neutrino luminosity in the accretion phase . while the quantitative results obviously depend on neutrino oscillation details ,
the qualitative trend with @xmath46 is unbroken .
comparing our predictions with the interactions observed from sn 1987a one notes ( but must keep the very small - number statistics in mind ) that either the explosion must have occurred early in the postbounce phase and/or the progenitor must have had a relatively low @xmath96 . stated another way , that data strongly disfavors a late - time explosion in a high compactness model .
both of these statements are broadly consistent with the previous work of @xcite and @xcite .
it is also not inconsistent with the proposed zams mass of @xmath5@xmath155 for sanduleak  @xmath156  @xmath157 ( e.g. , @xcite ) , the blue supergiant progenitor star of sn 1987a .
however , the mapping between zams mass and stellar structure ( i.e. , compactness ) at the presupernova stage appears to be highly non - monotonic ( cf .
[ fig : xis ] ; @xcite ) .
this makes it very difficult to link the observed neutrino signal to zams mass without additional constraints from classical astronomical observations in the electromagnetic spectrum .    in the discussion of sn 1987a , we have ignored uncertainties regarding nuclear eos , rotation of the progenitor core , distance , and collective neutrino oscillations .
all affect the preexplosion neutrino signature in ways that may be degenerate with variations in the compactness parameter .
however , as we have shown , large variations in the stiffness of the nuclear eos in the density and temperature regime relevant in the preexplosion phase can be disentangled from @xmath46 via observations of the emitted neutrino spectrum .
rapid rotation , which decreases both ( angle - averaged ) luminosities and average neutrino energies , can be constrained by gravitational wave observations of a galactic event ( e.g. , @xcite ) .
distance uncertainties from electromagnetic observations , which translate into uncertainties in the absolute luminosities , can be reduced by exploiting the generic neutronization burst as a standard candle @xcite .
collective neutrino oscillations are not yet fully understood and have not been directly incorporated in neutrino radiation - hydrodynamics simulations .
they may lead to one or multiple energy - dependent swaps of spectra between flavors and could thus complicate the mapping between neutrino signal and compactness parameter .
however , the recent understanding suggests that collective oscillations may not be significant in the preexplosion phase ( @xcite ; but : @xcite and @xcite ) .    our goal with this study was to highlight overall trends of the preexplosion neutrino signal with progenitor structure in the limit of spherical symmetry .
we expect these overall trends to be robust and to carry over to the multi - dimensional case .
we did not aim at making precise and robust quantitative predictions for any individual model .
these are not possible with our current spherically - symmetric radiation - hydrodynamics treatment , which neglects inelastic scattering , redshift and velocity dependent terms ( e.g. , @xcite ) .
rotation is included only in an angle - averaged manner and other important multi - dimensional dynamics , in particular convection and the standing - accretion shock instability @xcite , can not be captured .
future work to remove these limitations will be necessary to produce the reliable signal predictions necessary for drawing quantitative conclusions from neutrino observations of the next nearby core - collapse supernova .
however , even with fully realistic and complete simulation codes , large parameter studies will be necessary to account for prevailing uncertainties in the nuclear eos and/or neutrino interaction physics .
we acknowledge helpful discussions with and input from john beacom , adam burrows , luc dessart , ken nomoto , ryan patterson , kate scholberg , mark vagins , and stan woosley .
cdo thanks the kavli institute for the physics and mathematics of the universe for hospitality while work on a draft of this article was carried out .
the computations were performed at caltech s center for advanced computing research on the cluster `` zwicky '' funded through nsf grant no .
phy-0960291 and the sherman fairchild foundation .
furthermore , computations were performed on louisiana optical network infrastructure computer systems under allocation loni_numrel06 , on the nsf xsede network under allocation tg - phy100033 , and on resources of the national energy research scientific computing center , which is supported by the office of science of the u.s .
department of energy under contract no .
de - ac02 - 05ch11231 .
eoc is supported in part by a post - graduate fellowship from the natural sciences and engineering research council of canada ( nserc ) .
this research is supported in part by the national science foundation under grant nos .
ast-0855535 and oci-0905046 , by the alfred p. sloan foundation , and by the sherman fairchild foundation .                      ,
a. , beroz , f. , carr , r. , duan , h. , friedland , a. , kaiser , n. , kneller , j. , moss , a. , reitzner , d. , scholberg , k. , webber , d. , & wendell , r. 2011 , snowglobes , revision 1373 , http://www.phy.duke.edu/~schol/snowglobes/