An Investigation of Neutrino-driven Convection and the Core Collapse Supernova Mechanism Using Multigroup Neutrino Transport
- Theoretical and Computational Physics Group, Oak Ridge National Laboratory, Oak Ridge, TN 37831-6373 (United States)
- Department of Physics and Astronomy, University of Tennessee, Knoxville, TN 37996-1200 (United States)
- Department of Physics and Astronomy, Vanderbilt University, Nashville, TN 37235 (United States)
- Department of Physics, Florida Atlantic University, Boca Raton, FL 33431-0991 (United States)
- Department of Physics, North Carolina State University, Raleigh, NC 27695-8202 (United States)
We investigate neutrino-driven convection in core collapse supernovae and its ramifications for the explosion mechanism. We begin with a postbounce model that is optimistic in two important respects: (1) we begin with a 15 {ital M}{sub {circle_dot}} precollapse model, which is representative of the class of stars with compact iron cores; (2) we implement Newtonian gravity. Our precollapse model is evolved through core collapse and bounce in one dimension using multigroup (neutrino energy{endash}dependent) flux-limited diffusion (MGFLD) neutrino transport and Newtonian Lagrangian hydrodynamics, providing realistic initial conditions for the postbounce convection and evolution. Our two-dimensional simulation begins at 12 ms after bounce and proceeds for 500 ms. We couple two-dimensional piecewise parabolic method (PPM) hydrodynamics to precalculated one-dimensional MGFLD neutrino transport. (The neutrino distributions used for matter heating and deleptonization in our two-dimensional run are obtained from an accompanying one-dimensional simulation. The accuracy of this approximation is assessed.) For the moment, we sacrifice dimensionality for realism in other aspects of our neutrino transport. MGFLD is an implementation of neutrino transport that simultaneously (1) is multigroup and (2) simulates with sufficient realism the transport of neutrinos in opaque, semitransparent, and transparent regions. Both are crucial to the accurate determination of postshock neutrino heating, which sensitively depends on the luminosities, spectra, and flux factors of the electron neutrinos and antineutrinos emerging from their respective neutrinospheres. By 137 ms after bounce, we see neutrino-driven convection rapidly developing beneath the shock. By 212 ms after bounce, this convection becomes large scale, characterized by higher entropy, expanding upflows and lower entropy, denser, finger-like downflows. The upflows reach the shock and distort it from sphericity. The radial convection velocities at this time become supersonic just below the shock, reaching magnitudes in excess of 10{sup 9} cm s{sup {minus}1}. Eventually, however, the shock recedes to smaller radii, and at {approximately}500 ms after bounce there is no evidence in our simulation of an explosion or of a developing explosion. Our angle-averaged density, entropy, electron fraction, and radial velocity profiles in our two-dimensional model agree well with their counterparts in our accompanying one-dimensional MGFLD run above and below the neutrino-driven convection region. In the convection region, the one-dimensional and angle-averaged profiles differ somewhat because (1) convection tends to flatten the density, entropy, and electron fraction profiles, and (2) the shock radius is boosted somewhat by convection. However, the differences are not significant, indicating that, while vigorous, neutrino-driven convection in our model does not have a significant impact on the overall shock dynamics. The differences between our results and those of other groups are considered. These most likely result from differences in (1) numerical hydrodynamics methods; (2) initial postbounce models, and, most important; (3) neutrino transport approximations. We have compared our neutrino luminosities, rms energies, and inverse flux factors with those from the exploding models of other groups. Above all, we find that the neutrino rms energies computed by our multigroup (MGFLD) transport are significantly lower than the values obtained by Burrows and coworkers, who specified their neutrino spectra by tying the neutrino temperature to the matter temperature at the neutrinosphere and by choosing the neutrino degeneracy parameter arbitrarily, and by Herant and coworkers in their transport scheme, which (1) is gray and (2) patches together optically thick and thin regions. The most dramatic difference between our results and those of Janka and M{umlt u}ller is exhibited by the difference in the net cooling rate below the gain radii: Our rate is 2{endash}3 times greater during the critical 50{endash}100 ms after bounce. We have computed the mass and internal energy in the gain region as a function of time. Up to {approximately}150 ms after bounce, we find that both increase as a result of the increasing gain region volume, as the gain and shock radii diverge. However, at all subsequent times, we find that the mass and internal energy in the gain region decrease with time in accordance with the density falloff in the preshock region and with the flow of matter into the gain region at the shock and out of the gain region at the gain radius. Therefore, we see no evidence in the simulations presented here that neutrino-driven convection leads to mass and energy accumulation in the gain region. (Abstract Truncated)
- OSTI ID:
- 642457
- Journal Information:
- Astrophysical Journal, Vol. 495, Issue 2; Other Information: PBD: Mar 1998
- Country of Publication:
- United States
- Language:
- English
Similar Records
The interplay between proto{endash}neutron star convection and neutrino transport in core-collapse supernovae
General Relativistic Effects in the Core Collapse Supernova Mechanism