We have calculated a grid of photospheric phase atmospheres of Type Ia supernovae (SNe Ia) with metallicities from ten times to one thirtieth the solar metallicity in the C+O layer of the deflagration model, W7. We have modeled the spectra using the multi-purpose NLTE model-atmosphere and spectrum-synthesis code, PHOENIX. We show models for the epochs 7, 10, 15, 20, and 35 days after explosion. When compared to observed spectra obtained at the approximately corresponding epochs these synthetic spectra fit reasonably well. The spectra show variation in the overall level of the UV continuum with lower fluxes for models with higher metallicity in the unburned C+O layer. This is consistent with the classical surface cooling and line blocking effect due to metals in the outer layers of C+O. The UV features also move consistently to the blue with higher metallicity, demonstrating that they are forming at shallower and faster layers in the atmosphere. The potentially most useful effect is the blueward movement of the Si II feature at 6150Å with increasing C+O layer metallicity. We also demonstrate the more complex effects of metallicity variations by modifying the 54 Fe content of the incomplete burning zone in W7 at maximum light. We briefly address some shortcomings of the W7 model when compared to observations. Finally, we identify that the split in the Ca H+K feature produced in W7 and observed in some SNe Ia is due to a blending effect of Ca II and Si II and does not necessarily represent a complex abundance or ionization effect in Ca II.

We have performed ab initio neutrino radiation hydrodynamics simulations in three and two spatial dimensions (3D and 2D) of core-collapse supernovae from the same 15 M progenitor through 440 ms after core bounce. Both 3D and 2D models achieve explosions, however, the onset of explosion (shock revival) is delayed by ∼100 ms in 3D relative to the 2D counterpart and the growth of the diagnostic explosion energy is slower. This is consistent with previously reported 3D simulations utilizing iron-core progenitors with dense mantles. In the ∼100 ms before the onset of explosion, diagnostics of neutrino heating and turbulent kinetic energy favor earlier explosion in 2D. During the delay, the angular scale of convective plumes reaching the shock surface grows and explosion in 3D is ultimately lead by a single, large-angle plume, giving the expanding shock a directional orientation not dissimilar from those imposed by axial symmetry in 2D simulations. We posit that shock revival and explosion in the 3D simulation may be delayed until sufficiently large plumes form, whereas such plumes form more rapidly in 2D, permitting earlier explosions.

We present four ab initio axisymmetric core-collapse supernova simulations initiated from 12, 15, 20, and 25 M  zero-age main sequence progenitors. All of the simulations yield explosions and havebeen evolved for at least 1.2 s after core bounce and 1 s after material first becomes unbound. These simulations were computed with our CHIMERA code employing RbR spectral neutrino transport, special and general relativistic transport effects, and state-of-the-art neutrino interactions. Continuing the evolution beyond 1 s after core bounce allows the explosions to develop more fully and the processes involved in powering the explosions to become more clearly evident. We compute explosion energy estimates, including the negative gravitational binding energy of the stellar envelope outside the expanding shock, of 0.34, 0.88, 0.38, and 0.70 Bethe (B≡10 51 erg) and increasing at 0.03, 0.15, 0.19, and 0.52 B s 1 -, respectively, for the 12, 15, 20, and 25 M  models at the endpoint of this report. We examine the growth of the explosion energy in our models through detailed analyses of the energy sources and flows. We discuss how the explosion energies may be subject to stochastic variations as exemplfied by the effect of the explosion geometry of the 20 M  model in reducing its explosion energy. We compute the proto-neutron star masses and kick velocities. We compare our results for the explosion energies and ejected Ni 56 masses against some observational standards despite the large error bars in both models and observations.

We present an overview of four ab initio axisymmetric core-collapse supernova simulations employing detailed spectral neutrino transport computed with our CHIMERA code and initiated from Woosley & Heger (2007) progenitors of mass 12, 15, 20, and 25 M ⊙ . All four models exhibit shock revival over ∼ 200 ms (leading to the possibility of explosion), driven by neutrino energy deposition. Hydrodynamic instabilities that impart substantial asymmetries to the shock aid these revivals, with convection appearing first in the 12 M ⊙ model and the standing accretion shock instability (SASI) appearing first in the 25 M ⊙ model. Three of the models have developed pronounced prolate morphologies (the 20 M ⊙ model has remained approximately spherical). By 500 ms after bounce the mean shock radii in all four models exceed 3,000 km and the diagnostic explosion energies are 0.33, 0.66, 0.65, and 0.70 Bethe (B = 10 51 ergs) for the 12, 15, 20, and 25 M ⊙ models, respectively, and are increasing. The three least massive of our models are already sufficiently energetic to completely unbind the envelopes of their progenitors (i.e., to explode), as evidenced by our best estimate of their explosion energies, which first become positive at 320, 380, and 440 ms after bounce. By 850 ms the 12 M ⊙ diagnostic explosion energy has saturated at 0.38 B, and our estimate for the final kinetic energy of the ejecta is ∼ 0.3 B, which is comparable to observations for lower-mass progenitors.

We present the gravitational wave signatures for a suite of axisymmetric core collapse supernova models with progenitor masses between 12 and 25 M⊙. These models are distinguished by the fact that they explode and contain essential physics (in particular, multi-frequency neutrino transport and general relativity) needed for a more realistic description. Thus, we are able to compute complete waveforms (i.e. through explosion) based on non-parameterized, first-principles models. This is essential if the waveform amplitudes and time scales are to be computed more precisely. Fourier decomposition shows that the gravitational wave signals we predict should be observable by AdvLIGO across the range of progenitors considered here. The fundamental limitation of these models is in their imposition of axisymmetry. Further progress will require counterpart three-dimensional models.

We present a comparison between several simulation codes designed to study the core-collapse supernova mechanism. We pay close attention to controlling the initial conditions and input physics in order to ensure a meaningful and informative comparison. Our goal is three-fold. First, we aim to demonstrate the current level of arXiv:1806.04175v1 [astro-ph.HE] 11 Jun 2018Global Core-Collapse Supernova Comparison 2 agreement between various groups studying the core-collapse supernova central engine. Second, we desire to form a strong basis for future simulation codes and methods to compare to. Lastly, we want this work to be a stepping stone for future work exploring more complex simulations of core-collapse supernovae, i.e., simulations in multiple dimensions and simulations with modern neutrino and nuclear physics. We compare the early (first ∼500 ms after core bounce) spherically-symmetric evolution of a 20M progenitor star from six different core-collapse supernovae codes: 3DnSNe-IDSA, AGILE-BOLTZTRAN, FLASH, Fornax, GR1D, and PROMETHEUS-VERTEX. Given the diversity of neutrino transport and hydrodynamic methods employed, we find excellent agreement in many critical quantities, including the shock radius evolution and the amount of neutrino heating. Our results provide an excellent starting point from which to extend this comparison to higher dimensions and compare the development of hydrodynamic instabilities that are crucial to the supernova explosion mechanism, such as turbulence and convection.

We present the gravitational waveforms computed in ab initio two-dimensional core collapse supernova models evolved with the Chimera code for progenitor masses between 12 and 25 M . All models employ multi-frequency neutrino transport in the ray-by-ray approximation, state-of-the-art weak interaction physics, relativistic transport corrections such as the gravitational redshift of neutrinos, two-dimensional hydrodynamics with the commensurate relativistic corrections, Newtonian self-gravity with a general relativistic monopole correction, and the Lattimer-Swesty equation of state with 220 MeV compressibility, and begin with the most recent Woosley-Heger nonrotating progenitors in this mass range. All of our models exhibit robust explosions. Therefore, our waveforms capture all stages of supernova development: 1) a relatively short and weak prompt signal, 2) a quiescent stage, 3) a strong signal due to convection and SASI activity, 4) termination of active accretion onto the proto-neutron star, and 5) a slowly increasing tail that reaches a saturation value. Fourier decomposition shows that the gravitational wave signals we predict should be observable by AdvLIGO for Galactic events across the range of progenitors considered here. The fundamental limitation of these models is in their imposition of axisymmetry. Further progress will require counterpart three-dimensional models, which are underway.

We have conducted a series of numerical experiments with the spherically symmetric, general relativistic, neutrino radiation hydrodynamics code Agile-BOLTZTRAN to examine the effects of several approximations used in multidimensional core-collapse supernova simulations. Our code permits us to examine the effects of these approximations quantitatively by removing, or substituting for, the pieces of supernova physics of interest. These approximations include: (1) using Newtonian versus general relativistic gravity, hydrodynamics, and transport; (2) using a reduced set of weak interactions, including the omission of non-isoenergetic neutrino scattering, versus the current state-of-the-art; and (3) omitting the velocity-dependent terms, or observer corrections, from the neutrino Boltzmann kinetic equation. We demonstrate that each of these changes has noticeable effects on the outcomes of our simulations. Of these, we find that the omission of observer corrections is particularly detrimental to the potential for neutrino-driven explosions and exhibits a failure to conserve lepton number. Finally, we discuss the impact of these results on our understanding of current, and the requirements for future, multidimensional models.