the collapse of a massive star down to a neutron star releases over @xmath1erg of potential energy .
just one percent of this energy is required to power the observed type ib / c and type ii supernovae .
most of the energy is emitted in the form of neutrinos .
how this small fraction of the energy is converted into explosion energy is still a matter of debate , but it is believed by most that neutrinos play a role in depositing energy above the collapsed core to drive this explosion ( see fryer 2003 for a review ) .
whether or not neutrinos are important for the explosion , they do provide astronomers a window into the explosion mechanism behind supernovae .
neutrinos also provide a window into the material conditions at the core of a supernova explosion . with temperatures above 10mev and densities above nuclear densities ,
core - collapse supernovae make ideal laboratories for nuclear physics . to extract information about nuclear physics from such laboratories
, we must understand the systematics in our supernova experiment .
when the explosion engine is active , convective motions in the engine ( herant et al . 1994 ) make it very difficult to interpret the neutrino signal .
an observed elevated neutrino luminosity could be a change in the neutrino opacity or it could be a difference in the convective engine .
one approach to eliminate issues with convection would be to wait until all convective activity ceases .
a number of studies have now been presented following the evolution of a collapsing star over 0.5 - 1.5s after bounce ( fryer & heger 2000 ; burrows et al .
2006 ; scheck et al . 2006 ; and , in 3-dimensions , fryer & young 2007 ) . in all these calculations ,
hydrodynamic motions near the proto - neutron star continue to dramatically affect the neutrino emission through the end of the simulations . to achieve a clean neutrino signal ,
we must wait until after the launch of the explosion .
we must also wait a few seconds after the launch of the explosion because the proto - neutron star may experience deep convection ( keil et al .
1996 ) . after this time
( roughly a few seconds to 10 - 20s ) , the supernova finally seems to have achieved an ideal condition as a physics laboratory ( reddy et al .
after 10 - 20s , the neutrino signal will be too weak to do much experimental science , so we truly are limited to this narrow time window .
unfortunately , even at these late times , nature does not allow completely pristine conditions .
material falling back from the supernova explosion ( `` fallback '' ) may well produce a new round of convection and confusion to our neutrino signal .
this fallback has been studied both in its important role in calculating the initial mass of the neutron star formed in a supernova explosion ( e.g. fryer & kalogera 2001 ) and in estimating the r - rpocess yields in supernovae fryer et al .
( 2006 ) . in this paper , we study its role in determining the neutrino luminosity arising after the first few seconds of a supernova . in  2 , we review the history of supernova fallback and present new calculations of supernova explosions estimating the fallback for a range of explosion energies and stellar masses .
 3 describes the 2-dimensional code used to model the neutrino emission from the this fallback and shows results for the suite of simulations run for this paper .
different than the multi - dimensional simulations modeling stellar collapse , these simulations do not start at the onset of collapse ( or bounce ) and end at the launch of the explosion .
instead , they start a 2 - 10s ( depending on the explosion energy , etc . from  2 ) after the launch of the explosion when material begins to fall back onto the neutron star .
we conclude with a brief discussion on how neutrino signals can be used to help better understand the supernova explosion , neutron star birth masses and neutrino cross - sections .
the idea of fallback was first discussed by colgate ( 1971 ) to overcome nucleosynthesis issues arising from the supernova ejection of neutron rich material produced in stellar cores ( arnett 1971 , young et al .
colgate argued that the inner layers of the ejected material would deposit its energy to the stellar material above it , ultimately reducing its energy below that needed to escape the neutron star , and it would fall back onto the neutron star . in such a scenario , one would expect the inner material to fall back quickly ( within the first few to ten seconds ) .
it was argued that this material ( the neutron rich material from the initial explosion ) would accrete onto the neutron star , alleviating any nucleosynthesis issues .    since this work by colgate
, supernova explosion calculations have confirmed that fallback does occur ( bisnovatyi - kogan & lamzin 1984 ; woosley 1989 ; fryer et al .
1999 ; macfadyen et al . 2001 ) and new arguments for the cause of this fallback were suggested .
for example , woosley ( 1989 ) argued that when the supernova shock decelerates in the hydrogen layers of the star , it sends a reverse shock that drives fallback .
such a model argues that the fallback will happen at late times , long after it can affect the neutrino luminosity from the cooling proto - neutron star .
it also suggested that close binary systems ( where a star s hydrogen envelope was removed prior its collapse ) might experience a very different amount of fallback than the amount of fallback in the collapse of a single star .    which is the true cause of the fallback ? and more to the point for neutrino emission , when does fallback occur ?
this confusion has mostly exists because , without a quantitatively reliable explosion mechanism , scientists have artificially driven supernova explosions to be able to study the results of these explosions ( e.g. supernova light curves , nucleosynthetic yields , and fallback ) .
much of the past work ( e.g. woosley 1989 ; fryer et al .
1999 ; macfadyen et al .
2001 ) used piston driven explosions . by moving the piston out , scientists artificially lowered the amount of fallback and delayed this fallback considerably . these calculations
all predicted that fallback would occur more than 100s after the launch of the explosion .
but the piston - driven explosion mechanism may not accurately model the nature of the fallback .
young & fryer ( 2007 ) found that piston - driven explosions produced both different fallback rates and different nucleosynthetic yields than energy - driven explosions of the same final energy .
the errors in the yields or light - curves are on par with slight changes in the explosion energy .
such small errors seemed unimportant in matching the observations .
for fallback , the differences are much more dramatic .
not only does the amount of fallback change , but the timescale at which the fallback occurs can change by more than an order of magnitude , moving a falback time of a few hundred seconds down to just 3 - 15s .
this changes the fallback accretion rate by over an order of magnitude and ultimately determines whether or not fallback is important in estimating the observed neutrino signal .
it also means that the the amount of fallback will not change between binary and single stars ( unless binary interactions change the internal structure of the star ) .    in this paper
, we focus on fallback calculations from energy - injected explosions .
energy - injection is much better at mimicking the currently - favored supernova mechnanisms .
let s take the convection - enhanced , neutrino - driven supernova engine as an example ( see fryer 2003 for a review ) . in this mechanism ,
the basic energy source is neutrinos ( either diffusing out of the proto - neutron star core or from newly accreting material ) that heat the atmosphere above the neutron star .
this atmosphere is topped by the accretion shock of the infalling star .
when enough energy is deposited in this region , the accretion shock will be pushed outward and an explosion occurs .
convection aids this mechanism by both allowing heated material to rise ( cooling by adiabatic expansion instead of forcing it to continue to heat until it can cool by neutrino emission ) and allowing shocked material at the top of the convective region to flow down to the neutron star surface to accrete onto the neutron star ( emitting neutrinos ) or be heating to be part of the rising bubble .    an energy injection method ,
although not mimicking the effects of convection directly , can mimick the basic tenets of this model - heating just above the neutron star surface to blow off the accretion shock .
ideally , once the region above the neutron star becomes more rarefied , the energy injection will essentially halt ( aside from a weak , by supernova standards , neutrino - driven wind ) .
some groups have gone so far as to only inject this energy through a neutrino flux ( frlich et al .
2006 ) . in this manner
, once the region becomes rarefied , the energy injection drops naturally .
a piston - driven explosion can not mimic this effect well unless very specific attention is paid to the input of the piston .
the very different fallback results from piston explosions demonstrate just how far off such explosions are from the energy drive of the standard neutrino model .
this does not mean that by using energy - injection , we can produce a definitive fallback estimate for a supernova explosion .
the explosion energy is one of our primary uncertainties in estimating the fallback rate .
supernova scientists have yet to agree on the exact mechanism behind the supernovae and we are far from achieving quantitative predictions of these explosions . with accurate light - curve and spectra calculations , we may be able to estimate the explosion energy for a particular supernova based on its observations . by studying a range of explosion energies for a given progenitor , we can provide a template for the neutrino flux ( from fallback ) arising from these systems , allowing us to possibly extract this effect and once more focus a study on the physics of dense , hot matter .
our fallback calculations will all be based on fallback estimates from energy - driven explosions .
we use the same multi - step technique employed in young et al .
( 2006 ) and young & fryer ( 2007 ) .
we start with progenitor stars modeled to collapse ( either from heger et al .
2000 or young et al .
these progenitors are mapped into the 1-dimensional core - collapse code from herant et al .
this code includes equations of state valid from stellar densities up to nuclear densities , a 3-flavor flux - limited diffusion neutrino transport scheme , and a simple nuclear network . with this code
, we follow the collapse and formation of the proto - neutron star and the propogation of the shock produced when the collapse halts due to nuclear forces .    after this shock stalls ( as it loses its energy via neutrino losses ) , we have a structure defined by a shocked `` atmosphere '' produced by the now - stalled bounce shock above a dense proto - neutron star .
the proto - neutron star typically has a baryonic mass between 1.1 and 1.3m@xmath2 .
the edge of the proto - neutron star is determined by the mass where the density drops below @xmath3 . at these densities
, there is typically a well - defined edge where the density drops from @xmath4 to @xmath5 over a very narrow mass cut .
the exact location of this edge depends upon the progenitor mass and probably the exact code used to model this collapse phase ( e.g. 1-dimensional versus multi - dimensional results ) .
in general , 1-dimensional calculations have not produced supernova explosions . with
the proto - neutron star removed , we have also removed the energy source for any explosion . to induce an explosion in 1-dimension
, we must source in energy .
the simulations here source the energy directly into the inner 15 cells ( roughly 0.1m@xmath2 ) of the star .
we keep this energy source on for a limited time ( between 50 and 300ms ) , varying the energy injection rate and time to produce a range of explosion energies .
young & fryer ( 2007 ) found that such energy sourcing was more flexible than a simple neutrino enhancement to modeling the full range of proposed explosion mechanisms .
in our calculations , we use the shorter ( 50ms ) injection duration for the lowest energy explosions and the lowest mass stars and the longer ( 300ms ) injection duration for the more energetic explosions .
we follow this explosion for 400s .
this allows us to follow the shock as it moves well into the star ( and , for binary systems , out of the star surface ) .
our simulation space includes enough matter to ensure that for this 400s duration , the shock is well within the simulation space ( there are at least 100 zones between the final shock postion and the outer zone of our star ) . in the case of the binary systems , we have included a mass loss estimated from the stellar models .
typically , we model this star with roughly 2000 zones .
young et al . (
2008 ) have done a resolution study , comparing 1000,2000 and 4000 zones .
they do not find an appreciable difference in the ejected ( and hence fallback ) mass based on this resolution .
it will make a large difference on how the fallback mass tapers off .
however , note that with the lagrangian code in this paper , we can not accurately calculate low fallback rates .
but , as we shall see in our 2-dimensional models , it may well be that fallback powered explosions might cut off this low - rate fallback .
we use the same 1-dimensional code to model the fallback .
as material falls back onto the edge of our proto - neutron star and its density rises above @xmath6 , we remove the particle and add its mass to our proto - neutron star .
although we could follow the accretion and neutrino cooling of this material to higher densities with our code , as the density rises , the sound speed increases and the cell size decreases , both of which cause the time step to decrease . to make this problem tractable , we choose a low enough density that the material accretes before dropping the timestep below 1 microsecond .
even so , these simulations typically take between 1 million to 10 million timesteps . as we shall find in  3 , the modeling the true behavior of this matter requires modeling the accretion in multi - dimensions .
since these simulations are focused on calculating the infall rate ( not true accretion rate ) , our assumptions do not introduce large errors in our analysis .    in this paper
, we present the results from a suite of 1-dimensional explosion models models using 3 different progenitor masses and explosion energies ( for a summary , see table  [ tab : fallback ] ) .
table  [ tab : fallback ] also shows the peak accretion rates and total mass accreted for these models . note that we predict a range of neutron star masses based on both progenitor mass and explosion energy , in agreement with fryer & kalogera ( 2001 ) . but remember that the neutron star masses assume that all of the fallback remains on the neutron star . as we shall see in  3 ,
some of this matter is re - ejected . for normal explosion energies ,
it is likely that 12 - 15m@xmath2 stars produce neutron stars with gravitational masses in the 1.3 - 1.5m@xmath2 range .
if the explosions are stronger , the gravitational masses may well be as low as 1.2m@xmath2 .
for weak explosions , or more - massive stars , the remnant masses may be so large that the neutron star collapses to a black hole .
figure  [ fig : macc ] shows the mass accreted in the first 15s for our models .
note that in all casses , fallback occurs almost immediately ( with some delays of 2 - 7 seconds ) .
unless large amounts of fallback occur ( above 1m@xmath2 ) , the fallback is likely to be mostly over after 10 - 15s . with over 0.1m@xmath2 falling back in 10 - 15s , accretion rates above 0.01m@xmath2s@xmath7
are expected , corresponding to neutrino luminosities in excess of @xmath8 .
the bottom line is that , for normal supernovae , fallback occurs , at least for stars with initial masses at or above 12m@xmath2 .
this fallback occurs early , so it will definitely play some role in the neutrino emission at the 3 - 15s timescale .
an additional feature of the progenitor that affects the fallback and the neutrino luminosity is the angular momentum in the progenitor .
figure  [ fig:1drot ] shows the angular momentum of the inner 4m@xmath2 of a star for a variety of stars with and without magnetic braking ( heger et al .
2000,2005 ) . for neutron star accretion ,
the relevant angular momentum is that within the inner 1.4 - 2.0m@xmath2 .
such mass zones have low - angular momenta : a few@xmath9 for stars with magnetic braking , a few@xmath10 for stars without .
our typical simulations use a value of @xmath11 .
these low angular momenta had little effect on our results , and we ran some simulations with twice that amount . for our black hole systems , we ran even higher angular momenta . for neutron stars ,
the angular momentum is not enough for this material to form a true disk in the star , but it can alter the downflow and it is this effect that we would like to study in this paper .    note
that although there is quite a bit of structure in the angular momenta ( caused by incomplete angular momentum transport across elemental boundaries ) , in general , the angular momenta of the stars increases as one moves to higher and higher mass shell .
this is one reason why black - hole forming systems are more likely to produce asymmetric explosions than typical neutron - star forming systems .
the increased angular momentum means that angular momentum plays a bigger role in shaping the explosion , possibly producing larger asymmetries .
our 1-dimensional calculations provide us with a fallback rate .
if we assume that any fallback material emits all of the potential energy released from its downfall at the moment it hits the proto - neutron star , we can estimate the neutrino luminosity from the fallback . the accretion of 0.1m@xmath2 masses onto a 10 km
, 1.4m@xmath2 neutron star over 10s would correspond to a neutrino luminosity of @xmath12 . for our more massive stars
, this fallback can be ten times higher , corresponding to a neutrino luminosity of @xmath13 for over 10s after the launch of the explosion .
assuming pair - annihilation dominates the neutrino emission , this neutrino luminosity would be nearly 50% electron and 50% anti - electron neutrinos .
typical neutron star luminosities after 1s are below @xmath14 and can be as low as a few @xmath15 ( e.g. bruenn 1987 , keil & janka 1996 ) . even at 1s , fallback can dominate the neutrino emission if the fallback is heavy .
after @xmath165s our fallback estimates argue that fallback neutrinos dominate the neutrino emission .
but we have made several assumptions in this estimate for the neutrino luminosity .
first , we have assumed that all of the potential energy released is immediately emitted in neutrinos .
this energy can go into heating the proto - neutron star which will then cool on a longer timescale .
this energy also may go into ejecting other infalling material . as the material falls onto the proto - neutron star , it is shock heated and can rise .
fryer et al . ( 1996,2006 ) found that these rising shocked bubbles can accelerate above the escape velocity and actually be re - ejected .
so it is quite possible that a fraction of the `` fallback '' material does not end up on the neutron star if we account for these multi - dimensional effects .
this re - ejected material both does not contribute to the total energy available to produce neutrinos , but it takes some fraction of the potential energy released to power its explosion .
finally , @xmath17 and @xmath18 neutrinos may also be released and we must include these neutrinos in our energy budget . in this section ,
we study the effects of the above assumptions to get a more accurate neutrino luminosity from fallback .
our primary concern with our simple estimate of fallback neutrino emission is the fact that we ignore multi - dimensional effects .
we have known for some time that if a compact object is accreting mass with considerable angular momentum ( and inefficient cooling ) , outflows occur ( e.g. blandford & begelman 1999 ) . but even if the angular momentum is minimal , if the compact object is a neutron star , a sizable fraction of the infalling material can be re - ejected ( fryer et al .
2006 ) . to study the fate of supernova fallback ,
we want to understand a number of effects . for example , how do the results vary with accretion rate ? but we would like to also understand the role of angular momentum and the differences between hot or cold neutron stars . finally , numerical effects , such as boundary conditions , are bound to play a role .
before we discuss the results of our simulations , let s discuss the code used for these calculations and our tests of the physical and numerical effects .
our code must model the physics of downflows and allow us to answer the questions in the preceding paragraph .
first and foremost , we must follow the evolution in a multi - dimensional manner . as a first step , we use the two - dimensional smooth particle hydrodynamics code described in fryer et al .
( 1996,2006 ) .
we model the region from 10,000 km above the proto - neutron star down to the proto - neutron star surface .
typical runs range from an initial set of @xmath19 particles moving up to 50 - 70,000 particles by the end of the simulation .
the code includes an equation of state valid from densities below 1g@xmath20 up to nuclear densities ( including an estimate of nuclear statistical equilibrium ) .
neutrino transport is followed using flux - limited diffusion neutrino scheme for three neutrino species ( herant et al .
the neutrino emission and cross - sections are also outlined in herant et al .
( 1994 ) .
our simulations are set to mimic the conditions 2 - 10s after the launch of the explosion when material begins falling onto the proto - neutron star .
the initial condition begins with a set of particles ranging from 1000 km up to 10,000 km . the material is given velocities set to the free - fall velocity ( in our 1-dimensional simulations , the infall velocity is within 1 - 10velocity ) at their radial position with densities set to give the desired accretion rate ( recall , @xmath21 .
with @xmath22 within 1 - 10% of the free - fall velocity , we can accurately set up our initial conditions : our 3 separate rates correspond to 0.001,0.01,0.1m@xmath2s@xmath7 ( table  [ tab : neut ] ) .
for these simulations , we model only constant fallback rates .
note that in nature , the fallback rate varies quite a bit .
but we are focusing our study on the neutrino emission at peak fallback rates , and this suite of simulations will definitely bracket the range of results .
the entropy of the infalling material is generally set to a few .
the results are fairly insensitive to this initial entropy as the shock resets the entropy .
boundary conditions are probably the biggest uncertainty in our calculations .
particles are fed in through the outer boundary and are allowed to accrete through the inner boundary and be ejected out of the outer boundary .
the infall through the outer boundary is determined by a fixed infall rate .
this infall is only altered if material is flowing out of this outer boundary . at any point where an outflow occurs
, the inflow is temporarily halted .
this models the effect of outflows choking off the accretion .
the inner boundary is more difficult .
ideally , we would model the matter until its density reaches neutron star densities ( @xmath23 ) and it is mostly deleptonized . at such time , the matter has lost most of its energy and we can be sure that we have accounted for the total neutrino luminosity
. however , the sound speeds at such densities and the size of our lagrangian particles near the proto - neutron star would decrease the timestep to fractions of a microsecond .
such small timesteps prohibit us from following the evolution of the fallback for more than a fraction of a second .
instead we opt to use slightly less demanding criteria for the removal of particles on the inner boundary .
our standard set of models accretes particles whose density rises above @xmath24 with electron fractions below 0.3 .
this means that we are assuming any further energy released by the matter as it accretes goes into heating the neutron star which will cool on longer timescales . in our suite of models ,
we include a test of this boundary condition and find that the total neutrino luminosity is not too sensitive to our assumption for accretion .
we assume the axis of rotation lies along our axis of symmetry in our 2-dimensional caclulation .
each particle is given an angular momentum and it retains that angular momentum for the duration of the calculation . the angular momentum for each particle is determined to fit the angular velocities of our progenitor stars ( fig .
2 ) , so material along the axis has a low angular momentum whereas material in the equator has the highest angular momentum . the specific angular momentum @xmath25 is given by @xmath26 where @xmath27 is the angular velocity taken from the stellar models .
the inner material starts below the angular momentum given by @xmath27 , but by the end of our simulations , most of the material near the neutron star surface has the full angular momentum set by this @xmath27 value .
there is no angular momentum transport and hence , technically , no heating from this transport .
however , the angular momentum will slow the materials inflow , breaking the symmetry in the downflow .
fryer & heger ( 2000 ) also found that the angular momentum alters the instability criterion , preventing convection in the equatorial region where the angular momentum gradient is highest .
to study the effect of a hot neutron star , we introduced a non - zero neutrino flux arising from our inner , neutron - star , boundary .
in particular , we are interested in how the neutrinos from a hot neutron star affect the hydrodynamics , so we have chosen rather large luminosities ( @xmath28 electron neutrino luminosities with energies of 10mev with corresponding @xmath29 anti - electron neutrino luminosities with energies of 15mev ) .
these neutrinos are a boundary source for our flux - limited diffusion transport scheme and transported out of the system , heating the inflowing material .
however , in most of our simulations , the neutrino optical depth is fairly low , and the total energy deposited is also minimal ( see  3.2 ) .    finally , there is always concern that the artificial viscosity used in sph is introducing spurious effects into our calculations . for the most part
, our studies have found this numerical artifact to play a small role in results studying core - collapse supernovae , and we do not expect it to play a large role in these calcualtions .
nonetheless , we have included a simulation where we have increased both the bulk and von neumann - richtmyer viscosities by a factor of 2 ( 3.0,6.0 respectively versus our standard values of 1.5,3.0 ) .    with both numerical and physical effects to study
, we have run a suite of simulations to test the dependence of the neutrino luminosity on the initial conditions and on the numerics .
a summary of this suite of models , along with their basic results , is summarized in table  [ tab : neut ] .
our base model has an accretion rate of 0.01m@xmath2 s@xmath7 with an angular momentum equal to that shown in the circle in figure  [ fig:1drot ] .
we also include higher rotating runs with angular momenta set by equating the angular momenta to the value denoted by the square in figure  [ fig:1drot ] . in both cases ,
the angular momentum is low , so we do nt expect the formation of a full accretion disk .
we study the inner boundary in a number of ways .
we include a simulation where the criteria for accretion is more strict : @xmath30 with electron fractions below 0.1 .
we also include a set of simulations with an absorbing boundary at 100 km ( a black hole boundary condition ) .
we study accretion rates 10 times higher and ten times lower than our canonical rate .
finally , since the accretion occurs at early times , we have also included a few simulations where the neutron star itself is still emitting neutrinos . in   [ sec : dynamics ] , we compare the results on the dynamics for this suite of calculations . in   [ sec : neutrinos ] we compare the resulting neutrino luminosities .      before we discuss the neutrino fluxes from this suite of calculations , let s compare the dynamics in the accretion .
we expect material to shock as it hits the proto - neutron star , and some of this shocked material will begin to rise , driving convection .
these simulations are 2-dimensional . even without angular momentum ,
perturbations along our symmetry axis would allow the instability to develop stronger along this axis .
but we have studied this in some detail in core - collapse calculations and have minimized this effect ( see , for example , fryer & heger 2000 ) . however , the fact that the angular momentum axis also lies along the symmetry axis drives convection along our axis of symmetry .
the growth of convective instabilities is stabilized by angular momentum gradients and the deceleration of material ( along with the fact that our initial condition has an angular momentum gradient which is strongest along the equator ) , the convection initially grows strongest in the poles .
figure  [ fig : comp - nsvbh ] shows the results of our standard calculation showing this outflow .
this figure shows the evolution at 0.15s .
the right panel shows the evolution of the corresponding `` black hole '' simulation with the absorptive boundary . the true innermost stable circular orbit for a slowly rotatign 3m@xmath2 black hole is closer to 20 km , so our 100 km absorption radius understimates the activity from a real black hole .
however , we can see that although the angular momentum alters the inflow , it is insufficient to stop it , and the material continues to accrete directly onto the black hole .    as our neutron star model evolves the outflow expands , constraining the accretion to a funnel roughly @xmath31 from the plane .
figure  [ fig : comp - late ] shows the evolution of our standard model to 0.3 and 0.45s . over half of the inflowing material
is ultimately ejected .
not only does the ejected material not contribute to the energy available for neutrino emission , some of the energy of that accreted material goes toward accelerating this ejecta .
such outflows have been studied for over a decade in supermassive black hole systems such as active galactic nuclei ( see blandford & begelman 1999 ) for a review .
the amount of potential energy released in the accretion of material onto a compact object is enormous . if this energy is tapped ( either by viscous forces or neutrino emission ) , it can drive an explosion .
rockefeller et al . ( 2007 ) and fryer et al . (
2006 ) have recently applied this physics to stellar - massed compact objects to better understand their simulations of accreting systems . in rotating black hole systems
, material hangs up in a disk .
viscous forces that transport out angular momentum also transport energy , driving outflows .
what prevents this ejection is cooling ( either via photon radiation in the case of supermassive black holes or neutrinos in the case of most collapsing systems ) .
of course , for black hole systems that do not have a hard surface preventing the inflow of material , high angular momentum is required to prevent the energy from all flowing directly into the black hole . in neutron star systems , outflows ( or at least vigorous convection ) have been expected for more than a decade ( chevalier 1989 ; fryer et al .
1996 ) .
if the angular momentum were high enough , the infalling material would hang up in a disk .
this is most important in the black hole systems .
we ran a series of models increasing the specific angular momentum ( @xmath25 ) a factor of 3 and , for our black hole models , a factor of 10 .
macfadyen & woosley ( 1999 ) found that the specific angular momentum must be at least @xmath32 ( a factor of 1000 higher than our standard @xmath33 value ) .
as we can see from figure  [ fig : comp - rot ] , our factor of 10 increase in @xmath25 is not enough to change the fate of material falling back on our `` black hole '' simulation with its large absorptive boundary . in a true black hole system , the angular momentum available would increase as we move beyond 3m@xmath2 ( well above the angular momenta shown in fig .
[ fig:1drot ] ) .
this is one reason why the collapsar grb model argues for systems where the black hole mass exceeds 3m@xmath2 when the angular momentum in the star is sufficient to produce an accretion disk .
neutrinos from these systems have been considered in detail elsewhere ( macfadyen & woosley 1999 , hungerford et al .
2006 , rockefeller et al . 2007 ) and we do not study them in this paper .    if the neutron star is still hot and emitting neutrinos , it can also alter the inflow of fallback .
figure  [ fig : comp - hot ] shows two simulations with varying amounts of neutrino energy ( 18 and 36@xmath34 ) arising off the neutron star surface .
the effect on the dynamics is minimal , although a slight increase in the velocity ( and hence position at a given time ) can be seen . as we shall see below
, the neutrinos emanating from the proto - neutron star surface will dominate the total neutrino flux in such cases .
but we do nt expect these high luminosities at 10s and the accretion luminosity will dominate at these later times .
but at early times , the modification of the downflow on the neutrino opacities is the dominant effect .
finally , we have varied the infall rate from 0.001 to 0.1m@xmath2s@xmath7 .
the effect this has on the dynamics is shown in figure  [ fig : comp - mdot ] . as with many of our models ,
the nature of the dynamics is not altered significantly by these changes .
but we shall see that the neutrino flux very much depends on the accretion rate .
we found in the last section that the dynamical behavior of the fallback was relatively insensitive to both uncertainties in the numerics as well as uncertainties initial conditions : e.g. fallback rate and rotation .
but what about the neutrino luminosity ?
first let s study the uncertainties in the numerics .
fig  [ fig : nunum ] shows the neutrino luminosity and mean energy for 3 different models studying the numerical effects ( both the artificial viscosity and the treatment of the inner boundary ) on the calculation .
the variations in the viscosity and the inner boundary lead to differences that are less than a factor of 2 in the neutrino luminosity and 5% variations in the neutrino mean energy .
many of these errors could be dominated by the explicit transport scheme used in these calculations and it is possible that these errors can be significantly diminished with implicit schemes that are more stable .
figure  [ fig : nunum ] also shows the results for a fast rotating ns model ( ns - rot2 ) . with the low angular momenta in our calculations based on the rotation velocities in the inner core material of massive stars ( fig .
[ fig:1drot ] ) , rotation does not alter the neutrino luminosity noticeably .
black hole systems depend more sensitively on the rotation because it is the angular momentum that prevents the material from accreting directly into our black hole . but as we expected from the fact that our angular momentum is too low to produce accretion disks , the neutrino emission from our black hole systems is negligible ( fig .  [
fig : nubh ] ) . in such low - angular momentum systems
, we expect essentially no emission after the collapse of the neutron star down to a black hole .
contrast this to typical collapsar conditions , where the material falling back onto the black hole has enough angular momentum to hang up in a disk and its falback accretion rate is high enough to produce high - density , high - temperature structures . in this case
, the neutrino emission from a black hole can be quite large - with luminosities on par with the neutrino burst of the original collapse ( macfadyen & woosley 1999 ; hungerford et al .
2006 ; rockefeller et al . 2007 ) .    in our models ,
the @xmath17 and @xmath18 neutrinos tend to be an order of magnitude lower than the corresponding electron neutrino fluxes .
so our assumption that most of the neutrinos are emitted as electron or anti - electron neutrinos is reasonably valid .
hungerford et al . ( 2006 ) and rockefeller et al .
( 2007 ) found that for collapsar models , the @xmath17 and @xmath18 neutrinos make up a sizable fraction of the total emission . at higher accretion rates ( and higher angular momentum in the infalling material ) , the fraction of the luminosity emitted in @xmath17 and @xmath18 neutrinos will likely increase .    the strongest dependency on our initial conditions in the mass accretion rate .
figure  [ fig : numdot ] shows the neutrino emission from 3 different accretion rates onto our neutron star surface .
our basic potential energy released estimate would argue that the neutrino luminosity scales with the accretion rate .
and this trend is basically true for our simulations .
however , note that the electron and anti - electron neutrinos for our highest accretion rate ( 0.1m@xmath2s@xmath7 ) case are not quite an order of magnitude higher than our standard accretion rate ( 0.01m@xmath2s@xmath7 ) .
this is because the electron neutrinos become trapped in this highest accretion rate case , lowering the escaping luminosity .
the @xmath17 and @xmath18 remain nearly an order of magnitude higher as they are not trapped in any of our simulations .
figure  [ fig : nuhot ] shows the neutrino emission for 3 different `` hot '' neutron star models .
fallback contributes an additional 10 - 20% of the luminosity on average for the ns2-hot2 model , 20 - 40% to the ns2-hot1 model , and over 50% to the ns1-hot model .
the neutrino energy is also altered by an amount comparable to the change in the luminosity .
although a hot neutron star may dominate the neutrino luminosity , fallback clearly can contribute a sizable fraction of the observed neutrino flux .    finally , note that in general , our predicted fallback neutrino energies are high .
the most notable exception is that the mean neutrino energy is 5 - 10mev cooler for our low accretion - rate simulations . as the accretion rate decreases , so too will the mean energy of the emitted neutrinos .
fallback ( of at least 0.1m@xmath2 ) is likely to occur in normal ( @xmath35erg ) supernova explosions from stars with initial masses of 12m@xmath2 or greater . with energy - injected ( more realistic than piston - driven ) explosion calculations ,
this fallback occurs quickly , generally in the first 15s and peak accretion rates in the 0.01 - 0.1m@xmath2s@xmath7 range should be expected .
if the engine stays active long after the launch of the supernova shock , the total fallback will be lower .
but neutrino engines weaken significantly quickly after the shock is launched and the material above the neutrinosphere ( absorbing the energy ) is ejected .
magnetic - driven explosions could be very different .
the mechanism for this fallback is essentially that proposed by colgate ( 1971 ) whereby the outflowing material decelerates as it pushes against the material above it , ultimately decreasing its velocity below the escape velocity and causing it to fall back onto the proto - neutron star .
this means that the bulk of the fallback is fairly insensitive to the structure on the outer envelope of the star ( so fallback is roughly the same whether or not the star is in a binary ) .    our multi - dimensional simulations of this fallback suggest that not all of this fallback material is actually incorporated into the proto - neutron star .
some of it flows out of the system .
the corresponding neutrino luminosity is also lower , both because less material is accreted and some of the energy released goes toward driving outflows .
but the neutrinos from fallback could still dominate the neutrino emission from a proto - neutron star 10s into the explosion .    at this time , observations of neutrinos from supernovae are limited to those of supernova 1987a ( hirata et al .
1987 , bionta et al .
1987 ) . with only 20 neutrinos over @xmath1615s ,
it is difficult to place many constraints on the fallback .
some authors have claimed that claimed that the late - time neutrinos could not be easily explained by a cooling neutron star ( e.g. suzuki & sato 1987 ) , but others have argued that the observed neutrino signal is consistent with neutrinos diffusing out of a hot proto - neutron star ( burrows & lattimer 1987 ; bruenn 1987 ) . in our neutron - star forming models ,
fallback ends within the first 10 - 15s , comparable to the duration of the observed neutrino burst and the flux is consistent with our more standard accretion rates .
the observations are consistent with fallback , but could easily be explained by a neutron star without any fallback .
given that the progenitor star is believed to be greater than 15m@xmath2 , fallback is likely to have occured .
but with the current errors in the time - dependent luminosity , it is difficult to determine whether or not accretion is taking place .
if we assume fallback accretion did occur , the fact that the mean neutrino energy appears to be dropping with time in the observations suggests that the accretion rate is dropping at 10s .
but neutrinos are still observed at 10s . unless the fallback material has considerable angular momentum ,
the compact remnant is not a black hole at 10s .
in addition , because the accretion rate is dropping dramatically at this time , it is unlikely that the remnant will accrete much additional mass .
it will remain a neutron star . in our fallback calculations , systems that accrete enough material to form a black hole are either already a black hole at 10s or still accreting rapidly at 10s , neither of which is supported by the observations .
but we really need a better neutrino signal to say anything definitive .    with such high accretion rates , fallback can easily dominate the neutrino luminosity after a few seconds , more than doubling the emission from the neutron star at thearly times .
if we wish to study neutrino opacities in the supernova explosion , we will have to be able to calculate this fallback .
it may be possible to estimate the fallback rate from detailed study of supernova light - curves .
fortunately , any system that is detected in neutrinos will have a wealth of data in photons of all wavelengths .
if the fallback rate is high , the infalling material will alter the position of the neutrinosphere , and the problem becomes nearly as complex as the supernova explosion modeling itself .
we have only touched the surface on the difficulties in modeling such systems . in this paper , we have not addressed the production or feedback of magnetic fields , which may definitely alter the fate of the fallback .
we have also not discussed the nuclear yields from the ejecta ( a first paper on this is by fryer et al .
2007 and we plan future projects studying this nucleosynthesis ) .
however , note that the old `` wind '' r - process picture will not work in any system with fallback ( stars more massive than @xmath36m@xmath2 ) .
the matter trajectories just ca nt be explained by the wind solution in the cases where fallback dominates the motion near the neutron star .
we also have not studied the actual accretion onto the neutron star in detail .
note that the final neutron star masses in table  [ tab : fallback ] assumed no mass ejecta .
but as we have found in this study , over half of the fallback material may be ejected , leading to smaller neutron star masses and a narrower mass range for these neutron stars .
it is a pleasure to thank patrick young , frank timmes and aimee hungerford for useful conversations on this project .
this project was funded in part under the auspices of the u.s .
dept . of energy , and supported by its contract w-7405-eng-36 to los alamos national laboratory , and by a nasa grant swif03 - 0047 .
12 & & & & + & 0.55 & 0.33 & 1.64 & @xmath370.02 & 120 + & 0.92 & 0.23 & 1.54 & @xmath370.01 & 86 + & 2.4 & 0.030 & 1.34 & @xmath370.003 & 11 + & 3.6 & 0.015 & 1.33 & @xmath370.001 & 5.6 + 15 & & & & + & 0.9 & 0.34 & 1.75 & @xmath370.05 & 130 + & 1.65 & 0.25 & 1.66 & @xmath370.05 & 93 + & 7.0 & 0.11 & 1.52 & @xmath370.02 & 41 + & 17 & 0.08 & 1.49 & @xmath370.01 & 30 + 23 & & & & + & 1.25 & 2.2 & 3.9 & @xmath370.15 & 820 + & 2.0 & 1.75 & 3.45 & @xmath370.1 & 650 + & 2.5 & 0.86 & 2.56 & @xmath370.1 & 320 + & 6.6 & 0.02 & 1.72 & @xmath370.01 & 7.5 +      ns3 & @xmath38 & @xmath39 & @xmath40 & 0 . &
+ bh3 & @xmath38 & @xmath41 km & @xmath40 & 0 . & + ns2 & @xmath42 & @xmath39 & @xmath40 & 0 . & + ns2-rot2 & @xmath42 & @xmath39 & @xmath43 & 0 . & + ns2-hot1 & @xmath42 & @xmath39 & @xmath40 & 18 . & + ns2-hot2 & @xmath42 & @xmath39 & @xmath40 & 36 . & + ns2alpha & @xmath42 & @xmath39 & @xmath40 & 0 .
& @xmath44 + ns2bound & @xmath42 & @xmath45 & @xmath40 & 0 . &
@xmath44 + bh2 & @xmath42 & @xmath46 & @xmath40 & 0 . &
+ bh2-rot0 & @xmath42 & @xmath46 & 0 . & 0 . & + bh2-rot2 & @xmath42 & @xmath46 & @xmath43 & 0 .
& + bh2-rot10 & @xmath42 & @xmath46 & @xmath47 & 0 . & + ns1 & @xmath48 & @xmath39 & @xmath40 & 0 . & + ns1hot & @xmath48 & @xmath39 & @xmath40 & 36 . & + bh1 & @xmath48 & @xmath46 & @xmath40 & 0 . &