Core-collapse supernovae are, despite their spectacular visual display,
neutrino events. Virtually all of the 10^53 ergs of gravitational binding
energy released in the formation of the nascent neutron star is carried away in
the form of neutrinos and antineutrinos of all three flavors, and these
neutrinos are primarily responsible for powering the explosion. This mechanism
depends sensitively on the neutrino transport between the neutrinospheres and
the shock. In light of this, we have performed a comparison of multigroup
Boltzmann neutrino transport (MGBT) and multigroup flux-limited diffusion
(MGFLD) in post-core bounce environments. Differences in the mean inverse flux
factors, luminosities, and RMS energies translate to heating rates that are up
to 2 times larger for Boltzmann transport, with net cooling rates below the
gain radius that are typically 0.8 times the MGFLD rates. These differences are
greatest at earlier postbounce times for a given progenitor mass, and for a
given postbounce time, greater for greater progenitor mass. The increased
differences with increased progenitor mass suggest that the net heating
enhancement from MGBT is potentially robust and self-regulated.Comment: 7 pages, 2 figures, 1 table; LaTex using iopconf.sty; To appear in:
  Proceedings of The Second Oak Ridge Symposium on Atomic & Nuclear
  Astrophysic

Bruenn, S. W.

Guidry, M. W.

Messer, O. E. B.

Mezzacappa, A.

English

arXiv

Core-collapse supernovae are, despite their spectacular visual display, neutrino events. Virtually all ({approximately}99%) of the 10{sup 53} power ergs of gravitational binding energy released in the formation of the nascent neutron star is carried away in the form of neutrinos and antineutrinos of all three flavors, and these neutrinos are primarily responsible for powering the explosion. This mechanism depends sensitively on the neutrino transport between the neutrinospheres and the shock. In light of this, the authors have performed a comparison of multigroup Boltzmann neutrino transport (MGBT) and (Bruenn`s) multigroup flux-limited diffusion (MGFLD) in post-core bounce environments. Their analysis concentrates on those quantities central to the postshock matter heating stemming from electron neutrino and antineutrino absorption, namely the neutrino luminosities, RMS energies, and mean inverse flux factors. The authors show that MGBT yields mean inverse flux factors in the gain region that are {approximately}25% larger and luminosities that are {approximately}10% larger than those computed by MGFLD. Differences in the mean inverse flux factors, luminosities, and RMS energies translate to heating rates that are up to 2 times larger for Boltzmann transport, with net cooling rates below the gain radius that are typically {approximately}0.8 times the MGFLD rates. These differences are greatest at earlier postbounce times for a given progenitor mass, and for a given postbounce time, greater for greater progenitor mass. The increased differences with increased progenitor mass suggest that the net heating enhancement from MGBT is potentially robust and self-regulated

Messer, O. E.

UNT Digital Library

ORNL/CP-98110 NEUTRINO-DRIVEN SUPERNOVAE: BOLTZMANN NEUTRINO TRANSPORT AND THE EXPLOSION MECHANISM O.E.B. Messert,t, A. Mezzacappat,S, S. W. Bruenn$, andM. W. Guidryt,$ ?Physics Division, Oak Ridge National Laboratory, Oak Ridge, TN 3783 1-6354 $Department of Physics and Astronomy, University of Tennessee, Knoxville, TN 37996 $Department of Physics, Florida Atlantic University, Boca Raton, FL 33431 to be published in Proceedings of Second Oak Ridge Symposium on Atomic and Nuclear As6ophysks Oak Ridge, Tennessee December 2-6, 1997 RECEIVED *Oak Ridge National Laboratory is managed by Lockheed Martin Energy Research Corp. for the U.S. Department of Energy under Contract No. DE-AC05-960R22464. “The submitted manuscript has been authored by a contractor of the U.S. Government under contract No. DE-AC05-960R22464. Accordingly, the U.S. Government retains a nonexclusive, royalty-free license to publish or reproduce the published form of this contribution, or allow others to do so, for U.S. Government purposes.” This repon was prepared as an account of work sponsored by an agency of the United States Government. Neither the United States Government nor any agency thereof, nor any of their employees, makes any warranty, express or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or use- fulness of any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any spe- cific commercial product, proms, or service by trade name, trademark, manufac- turer, or otherwise does not necessarily constitute or imply its endorsement, m m -  mendation, or favoring by the United States Government or any agency thereof. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States Government or any agency thereof. Neutrino-driven supernovae: Boltzmann neutrino transport and the explosion mechanism 0. E. B. MesserT$ A. Memacappat$ S. W. Bruenns, and M. W. Guidryt$ t Physics Division, Oak Ridge National Laboratory, Oak Ridge, TN 37831 $ Department of Physics and Astronomy, University of Tennessee, Knoxville, TN 37996 5 Department of Physics, Florida Atlantic University, Boca Raton, FL 33431 Abstract. Core-collapse supernovae are, despite their spectacular visual dis- play, neutrino events. Virtually all (-99%) of the ergs of gravitational binding energy released in the formation of the nascent neutron star is carried away in the form of neutrinos and antineutrinos of all three flavors, and these neutrinos are primarily responsible for powering the explosion. This mechanism depends sensi- tively on the neutrino transport between the neutrinospheres and the shock. In light of this, we have performed a comparison of multigroup Boltzmann neutrino transport (MGBT) and (Bruenn’s) multigroup flux-limited diffusion (MGFLD) in post-core bounce environments. Our analysis concentrates on those quantities central to the postshock matter heating stemming from electron neutrino and an- tineutrino absorption, namely the neutrino luminosities, RMS energies, and mean inverse flux factors. We show that MGBT yields mean inverse flux factors in the gain region that are -25% larger and luminosities that are -10% larger than those computed by MGFLD. Differences in the mean inverse flux factors, luminosities, and RMS energies translate to heating rates that are up to 2 times larger for Boltz- mann transport, with net cooling rates below the gain radius that are typically -0.8 times the MGFLD rates. These differences are greatest at earlier postbounce times for a given progenitor mass, and for a given postbounce time, greater for greater progenitor mass. The increased differences with increased progenitor mass suggest that the net heating enhancement from MGBT is potentially robust and self-regulated. . 1. Introduction Ascertaining the core collapse supernova mechanism is a long-standing problem in astro- physics. The current paradigm begins with the collapse of a massive star’s iron core and the generation of an outwardly propagating shock wave that results from core rebound. Because of nuclear dissociation and neutrino losses, the shock stagnates. This sets the 1 stage for a shock reheating mechanism whereby neutrino energy deposition via electron neutrino and antineutrino absorption on nucleons behind the shock reenergizes it (Bethe & Wilson 1985; Wilson 1985). The shock reheating phase is essential to the supernova’s success, but it is pre- cisely this phase that is difficult to simulate realistically. During shock reheating, core electron neutrinos and antineutrinos are radiated from their respective neutrinospheres, and a small fraction of this radiated energy is absorbed in the exterior shocked mantle. The shock reheating depends sensitively on the electron neutrino and antineutrino lu- minosities, spectra (best characterized by the RMS energies), and angular distributions in the region behind the shock (e.g., see Burrows & Goshy 1993, Janka & Muller 1996, Mezzacappa et al 1998). These, in turn, depend on the neutrino transport in the semitransparent region encompassing the neutrinospheres, necessitating a neutrino transport treatment that is able to transit accurately and seamlessly between neutrino-thick and neutrino-thin regions. Various neutrino transport approximations have been implemented in simulating core collapse supernovae. The most sophisticated approximation, which naturally has been used in realistic one-dimensional simulations, is multigroup flux-limited diffusion (MGFLD; e.g., Bowers & Wilson 1982, Bruenn 1985, Myra et al 1987). MGFLD closes the neutrino radiation hydrodynamics hierarchy of equations at the level of the first moment (the neutrino flux) by imposing a relationship between the flux and the gradient of the neutrino energy density (the zeroth moment). For example, CA aU, 3 ar + ..., F, = --- where X is the neutrino mean free path, and U, and F, are the neutrino energy density and flux (Bruenn 1985). [Other forms for the flux-limiter A can be found in Bowers & Wilson (1982), Levermore & Pomraning (1981), and Myra et al (1987).] Whereas the limits X + 0 and X + 00 produce the correct diffusion and free streaming fluxes, it is in the critical intermediate regime where the MGFLD approximation is of unknown accuracy. Unfortunately, the quantities central to the postshock neutrino heating, i.e., the neutrino RMS energies, luminosities, and mean inverse flux factors, are determined in this regime, and given the sensitivity of the neutrino heating to these quantities, it becomes necessary to consider more accurate transport schemes. Moreover, in de- tailed one-dimensional simulations that have implemented elaborate MGFLD neutrino transport (e.g., see Bruenn 1993, Wilson & Mayle 1993, and Swesty & Lattimer 1994), explosions were not obtained unless the neutrino heating was boosted by additional phe- nomena, such as convection. This leaves us with at least two possibilities to consider: (1) Failures to produce explosions in the absence of additional phenomena, such as convec- tion, have resulted from neutrino transport approximation. (2) Additional phenomena may be essential in obtaining explosions. 2. Initial Models, Codes, and Methodology We begin with 15 Ma and 25 Ma precollapse models S15s7b and S25s7b provided by Woosley (1995). The initial models were evolved through core collapse and bounce using one-dimensional Lagrangian hydrodynamics and MGFLD neutrino transport coupled to the Lattimer-Swesty equation of state (Lattimer & Swesty 1991). The data at 106 ms and 233 ms after bounce for S15s7b and 156 ms after bounce for S25s7b were thermally and hydrodynamically frozen. Stationary-state neutrino distributions were computed for these profiles using both MGBT and MGFLD. The MGBT simulations were performed using BOLTZTRAN: a Newtonian gravity, O(v/c), three-flavor, Boltzmann neutrino transport code developed for the supernova problem and used thus far for studies of stellar core collapse (Mezzacappa and Matzner 1989, Mezzacappa & Bruenn 1993abc). The MGFLD simulations were performed using MGFLD-TRANS: a Newtonian gravity, O(v/c), three-flavor, MGFLD neutrino transport code, which has been used for both core collapse and postbounce evolution (e.g., Bruenn 1985, 1993). The MGBT simulations used 110 nonuniform radial spatial zones and 12 neutrino energy zones spanning a range between 5 and 225 MeV. The MGFLD used the same spatial and energy grids. Simulations with 20 energy zones spanning the same energy range were performed with BOLTZTRAN; no changes in the results presented here were seen. For the MGBT simulations there is an added dimension: neutrino direction cosine. Because MGBT computes the neutrino distributions as a function of direction cosine and energy for each spatial zone, the isotropy of the neutrino radiation field as a function of radius and neutrino energy is computed from first principles. This is one of the key features distinguishing MGBT and MGFLD. Because the isotropy of the neutrino radiation field is critical to the shock reheating/revival, four Gaussian quadrature sets (2-, 4-, 6-, and 8-point) were implemented in the MGBT simulations to ensure numerical convergence of the results. 3. Results For electron neutrino and antineutrino absorption on neutrons and protons, the neutrino heating rate (in MeV/nucleon) in the region between the neutrinospheres and the shock can be written as 1 F < EEe >< = > < E:e >< - > +-- 1 x p  Le Xn Lue E =  --A; 4 m 2  F A: 47rr2 (3) where: A t  = ii = G;p(g$ + 3g : ) /~ (hc )~m~;  G F / ( ~ c ) ~  = 1.166 x GeVF2 is the Fermi coupling constant; p is the matter density; gv = 1.0, g A  = 1.23; mB is the baryon 'mass; Xn,p are the free neutron and proton mass fractions; and L,,,, < E:e,De >, and F, F are the electron neutrino and antineutrino luminosities, RMS energies, and mean inverse flux factors, defined by 19.0 F Q, 17.0 E. $> 15.0 35.0 'u) 30.0 F c! V - r Q, 0 g 25.0 7 3 20.0 2.0 A %> -1.0 30.0 80.0 130.0 180.0 230.0 radius [km] Figure 1. RMS energies, luminosities, and mean inverse flux factors for model S15s7b at 233 ms after core bounce. In equations (4)-(6), f is the electron neutrino distribution function, which is a function of the electron neutrino direction cosine, pve, and energy, Eve. In equation (6), Vue and Fye are the electron neutrino energy density and flux. Corresponding quantities can be defined for the electron antineutrinos. Success or failure to generate explosions via neutrino reheating must ultimately rest on the three quantities defined in equations In Figure 1, at 233 ms after bounce for model S15s7b7 we plot the electron neutrino and antineutrino RMS energies, luminosities, and mean inverse flux factors as a function of radius for our ($-point Gaussian quadrature) MGBT and MGFLD runs. The energy- averaged neutrinospheres (at 57 km and 48 km, for electron neutrinos and antineutrinos, respectively), and the location of the shock (at 191 km), are indicated by arrows. The gain radius (neutrino-energy integrated), located at 98 km, is also marked by an arrow. For the electron neutrinos, the differences in RMS energies between MGBT and MGFLD are at  most 2% throughout most of the region plotted, although MGBT consistently gives higher energies. The differences between MGBT and MGFLD antineutrino RMS energies are smaller, and neither transport scheme yields consistently higher values. It should be noted that we expect larger differences when a fully hydrodynamic simulation is carried out, with MGBT giving harder spectra (Mezzacappa and Bruenn 1993a,c; see (4)-(6) * A . also Burrows 1998). In a static matter configuration, differences that result from different treatments of the neutrino energy shift measured by comoving observers do not occur. Significant differences between MGBT and MGFLD are evident when comparing the neutrino and antineutrino luminosities and mean inverse flux factors. Both transport methods compute similar electron neutrino luminosities until the- neutrinospheres are approached from below. The antineutrino luminosities for each transport method also coincide up to  this point. Just below the neutrinospheres, the MGBT luminosities diverge upward from the MGFLD luminosities, as MGFLD underestimates the neutrino flux, differing by 7% (4% for antineutrinos) at the neutrinospheres. After a decline from this maximum difference, the fractional difference grows from approximately 3% at the base of the gain region to a constant difference of 6% beyond about 170 km. Similar behavior is exhibited by the antineutrino luminosities, with the same fractional differences, 3% and 6%) obtained at the base of the gain region and near the shock, respectively. For the electron neutrinos, the fractional difference between < 1/F >MGFLD and < 1/F >MGBT is 2%, 8%, and 12% at the neutrinosphere, gain radius, and shock, respec- tively. Just above the shock, the difference converges to 10% and is maintained to the edge of the core. Focusing on the semitransparent region, < 1/F >MGFLD is greater below the gain radius until the gain radius is approached; i.e., the MGFLD neutrino radiation field is more isotropic than the MGBT radiation field in this region. At 80 km, as the gain radius is approached, MGFLD computes a sharp decrease in < 1/F >. This sharp decrease marks the radius at which the electron neutrino source is enclosed. The dip at 106 km and the sharp spike at 163 km in < 1/F >MGFLD are caused by local density perturbations. For the electron antineutrinos, the same features are seen in < 1/F >MGFLD. The fractional difference is 0%, 11%, and 11% at the neutrinosphere, gain radius, and shock, respectively. The initial sharp decrease in < 1/F >MGFLD occurs at a smaller radius; the point at which the electron antineutrino source is enclosed is at a smaller radius. These results are typical of all three time slices. Because each of the quantities plotted in Figure 1 is consistently greater for MGBT (while this is not strictly true for the antineutrino RMS energies in our stationary state comparisons, in a fully dynamical simulation these energies will be consistently higher for MGBT [Mezzacappa and Bruenn 1993a,c; see also Burrows 1998]), and because the neutrino heating rate is proportional to  each of them, MGBT yields a significantly higher heating rate. As an example, just above the gain radius for model S15s7b, at tpb = 233 ms and the net-heating peak, MGBT yields a heating rate from neutrino absorption that is (102%)2 x 110% x 112% of the MGFLD rate. In Figure 2, for MGBT and MGFLD, we plot the net neutrino heating rates as a function of radius for model S15s7b at tpb = 233 ms. (As discussed in Section 2, the results from four Gaussian quadrature sets are plotted to demonstrate numerical convergence.) These rates include the contributions from both the electron neutrinos and antineutrinos, and were computed using the following formulae: where E is the internal energy per gram; E,, $!, A!’, and j i  are the electron neutrino or antineutrino energy, zeroth angular moment, absorption mean free path, and emissivity, 5.0e+20 0) C m 0) .- e BOLTZTRAN (2 point GQ) ' BOLTZTRAN (4 point GQ) _. BOLTZTRAN (6 point GQ) --- v, sphere ---- - i BOLTZTRAN (8 point GQ) , - - MGFLD 50.0 100.0 150.0 200.0 radius [km] Figure 2. Net heating rates for model S15s7b at 233 ms after core bounce. respectively; i = 1 corresponds to electron neutrinos, and i = 2 corresponds to  elec- tron antineutrinos. Only the contributions from neutrino emission and absorption were included. The MGBT simulation yields a net heating rate just above the gain radius that is -1.3 times the MGFLD rate, and a net cooling rate below the gain radius that is consistently -0.8 times the MGFLD rate. The differences in net heating rate are even greater for the tpb = 106 ms time slice in our 15 Ma model and for the tpb = 156 ms time slice in our 25 Ma model (cf. Table 1 and Messer et al 1998). 4. Summary, Discussion, and Conclusions Comparing three-flavor MGBT and t hree-flavor MGFLD in post bounce supernova envi- ronments, we find that MGBT leads to a significant increase/decrease in the net heat- ing/cooling rate, particularly just above/below the gain radius. The MGBT net heating rate can be as much as 2 times the MGFLD net heating rate above the gain radius, with net cooling rates that are typically 0.8 times the MGFLD rate below. These differences stem primarily from differences in the neutrino luminosities and mean inverse flux fac- tors; the heating rate is linearly proportional to both these quantities, and differences in both add to produce a significant difference in the net heating rate. We also observe that the differences in the net heating rate are greatest at earlier postbounce times for a given progenitor mass, and at a given postbounce time, greater for greater progenitor mass. This is illustrated in Table 1. The enhancement in heating with increased progenitor mass suggests that the net heating enhancement from MGBT is potentially robust and self-regulated. In closing, our results are promising, and their ramifications for core collapse super- novae, and in particular, for the postbounce neutrino-heating, shock-revival mechanism, await one- and two-dimensional dynamical simulations with MGBT coupled to the core hydrodynamics. One-dimensional simulations are currently underway, and we plan to  report on them soon. Table 1. Maximum Net Heating/Cooling Rates ~~ Progenitor Mass [Ma] t p b  [ms] Maximum Net Heating Ratio Maximum Net Cooling Ratio 15 106 2.0 0.8 233 1.3 0.8 25 156 2.0 0.8 5. Acknowledgements BM was supported at the University of Tennessee under NASA grant NRA-96-04-GSFC-073. AM and MWG were supported at the Oak Ridge National Laboratory, which is managed by Lockheed Martin Energy Research Corporation under DOE contract DE-AC05-960R22464, and at the University of Tennessee, under DOE contract DE-FG05-93ER40770. SWB was supported at Florida Atlantic University under NASA grant NRA-96-04-GSFC-073 and NSF grant AST-9618423. The simulations presented in this paper were carried out on the Cray C90 at the National Energy Research Supercomputer Center, the Cray Y/MP at the North Carolina Supercomputer Center, and the Cray Y/MP and Silicon Graphics Power Challenge at the Florida Supercomputer Center. AM and SWB gratefully acknowledge the hospitality of the Institute for Theoretical Physics, Santa Barbara, which is supported in part by the National Science Foundation under grant number PHY94-07194. References Bethe H & Wilson J R 1985 ApJ295 14 Bowers R L & Wilson J R 1982 ApJS 50 115 Bruenn S W 1985 ApJS 58  771 Bruenn S W 1993 Nuclear Physics in the Universe (Bristol: IOP Publishing) p 31 Burrows A 1998 Proceedings of the 18th Texas Symposium on Relativistic Astrophysics Burrows A & Goshy J 1993 ApJ416 L75 Janka H-Th, & Muller, E. 1996 A&A 306 167 Lattimer J M & Swesty F D 1991 Nucl. Phys. A 535 331 Levermore C D & Pomraning G C 1981 ApJ 248 321 Messer 0 E B, Mezzacappa A, Bruenn S W, & Guidry M W 1998 ApJ submitted Mezzacappa A & Bruenn S W 1993a A p J 4 0 5  637 (Singapore: World Scientific Press) in press -1993b ApJ 405 669 -1993~ ApJ 410 710 Mezzacappa A, Calder A C, Bruenn S W, Blondin J M, Guidry M W, Strayer M R, & Umar A S ApJ495 911 Mezzacappa A & Matzner R A 1989 ApJ 343 853 Myra E, Bludman S, Hoffman Y, Lichtenstadt I, Sack N & Van Riper K 1987 ApJ 318 744 Swesty F D & Lattimer J M 1994 ApJ425 195 Wilson J R 1985 Numerical Astrophysics (Boston: Jones & Bartlett) p 422 Wilson J R & Mayle R W 1993 Phys. Rep. 227 97 Woosley S E 1995 private communication 7 M98005656 I11111111 Il 111ll lll 11111l llll 111111 I111 Publ. Date (11) Sponsor Code (1 8) U C Category (1 9) 19980702 046 DOE 

Neutrino-driven supernovae: Boltzmann neutrino transport and the explosion mechanism

Abstract

Similar works

Full text

Available Versions

UNT Digital Library