skip to main content


Title: Comparison of Electron Capture Rates in the N = 50 Region using 1D Simulations of Core-collapse Supernovae
Abstract

Recent studies have highlighted the sensitivity of core-collapse supernovae (CCSNe) models to electron-capture (EC) rates on neutron-rich nuclei near theN= 50 closed-shell region. In this work, we perform a large suite of one-dimensional CCSN simulations for 200 stellar progenitors using recently updated EC rates in this region. For comparison, we repeat the simulations using two previous implementations of EC rates: a microphysical library with parametrizedN= 50 rates (LMP), and an older independent-particle approximation (IPA). We follow the simulations through shock revival up to several seconds post-bounce, and show that the EC rates produce a consistent imprint on CCSN properties, often surpassing the role of the progenitor itself. Notable impacts include the timescale of core collapse, the electron fraction and mass of the inner core at bounce, the accretion rate through the shock, the success or failure of revival, and the properties of the central compact remnant. We also compare the observable neutrino signal of the neutronization burst in a DUNE-like detector, and find consistent impacts on the counts and mean energies. Overall, the updated rates result in properties that are intermediate between LMP and IPA, and yet slightly more favorable to explosion than both.

 
more » « less
Award ID(s):
2209429 1927130
NSF-PAR ID:
10377275
Author(s) / Creator(s):
; ; ; ; ; ;
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
939
Issue:
1
ISSN:
0004-637X
Format(s):
Medium: X Size: Article No. 15
Size(s):
["Article No. 15"]
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT

    The explosion outcome and diagnostics of core-collapse supernovae depend sensitively on the nature of the stellar progenitor, but most studies to date have focused exclusively on one-dimensional, spherically symmetric massive star progenitors. We present some of the first core-collapse supernovae simulations of three-dimensional massive star supernovae progenitors, a 12.5- and a 15-M⊙ model, evolved in three dimensions from collapse to bounce through explosion with the radiation-hydrodynamic code fornax. We compare the results using those starting from three-dimensional progenitors to three-dimensional simulations of spherically symmetric, one-dimensional progenitors of the same mass. We find that the models evolved in three dimensions during the final stages of massive star evolution are more prone to explosion. The turbulence arising in these multidimensional initial models serves as seed turbulence that promotes shock revival. Detection of gravitational waves and neutrinos signals could reveal signatures of pre-bounce turbulence.

     
    more » « less
  2. ABSTRACT

    We investigate the impact of strong initial magnetic fields in core-collapse supernovae of non-rotating progenitors by simulating the collapse and explosion of a $16.9\, \mathrm{M}_\odot$ star for a strong- and weak-field case assuming a twisted-torus field with initial central field strengths of ${\approx }10^{12}$ and ${\approx }10^{6}\, \mathrm{G}$. The strong-field model has been set up with a view to the fossil-field scenario for magnetar formation and emulates a pre-collapse field configuration that may occur in massive stars formed by a merger. This model undergoes shock revival already $100\, \mathrm{ms}$ after bounce and reaches an explosion energy of $9.3\times 10^{50}\, \mathrm{erg}$ at $310\, \mathrm{ms}$, in contrast to a more delayed and less energetic explosion in the weak-field model. The strong magnetic fields help trigger a neutrino-driven explosion early on, which results in a rapid rise and saturation of the explosion energy. Dynamically, the strong initial field leads to a fast build-up of magnetic fields in the gain region to 40 per cent of kinetic equipartition and also creates sizable pre-shock ram pressure perturbations that are known to be conducive to asymmetric shock expansion. For the strong-field model, we find an extrapolated neutron star kick of ${\approx }350\, \mathrm{km}\, \mathrm{s}^{-1}$, a spin period of ${\approx }70\, \mathrm{ms}$, and no spin-kick alignment. The dipole field strength of the proto-neutron star is $2\times 10^{14}\, \mathrm{G}$ by the end of the simulation with a declining trend. Surprisingly, the surface dipole field in the weak-field model is stronger, which argues against a straightforward connection between pre-collapse fields and the birth magnetic fields of neutron stars.

     
    more » « less
  3. Abstract

    We compare the core-collapse evolution of a pair of 15.8Mstars with significantly different internal structures, a consequence of the bimodal variability exhibited by massive stars during their late evolutionary stages. The 15.78 and 15.79Mprogenitors have core masses (masses interior to an entropy of 4kBbaryon−1) of 1.47 and 1.78Mand compactness parametersξ1.75of 0.302 and 0.604, respectively. The core-collapse simulations are carried out in 2D to nearly 3 s postbounce and show substantial differences in the times of shock revival and explosion energies. The 15.78Mmodel begins exploding promptly at 120 ms postbounce when a strong density decrement at the Si–Si/O shell interface, not present in the 15.79Mprogenitor, encounters the stalled shock. The 15.79Mmodel takes 100 ms longer to explode but ultimately produces a more powerful explosion. Both the larger mass accretion rate and the more massive core of the 15.79Mmodel during the first 0.8 s postbounce time result in largerνe/ν¯eluminosities and RMS energies along with a flatter and higher-density heating region. The more-energetic explosion of the 15.79Mmodel resulted in the ejection of twice as much56Ni. Most of the ejecta in both models are moderately proton rich, though counterintuitively the highest electron fraction (Ye= 0.61) ejecta in either model are in the less-energetic 15.78Mmodel, while the lowest electron fraction (Ye= 0.45) ejecta in either model are in the 15.79Mmodel.

     
    more » « less
  4. Abstract

    We study in detail the ejecta conditions and theoretical nucleosynthetic results for 18 three-dimensional core-collapse supernova (CCSN) simulations done by Fornax. Most of the simulations are carried out to at least 3 s after bounce, which allows us to follow their longer-term behaviors. We find that multidimensional effects introduce many complexities into the ejecta conditions. We see a stochastic electron fraction evolution, complex peak temperature distributions and histories, and long-tail distributions of the time spent within nucleosynthetic temperature ranges. These all lead to substantial variation in CCSN nucleosynthetic yields and differences from 1D results. We discuss the production of lighterα-nuclei, radioactive isotopes, heavier elements, and a few isotopes of special interest. Comparing pre-CCSN and CCSN contributions, we find that a significant fraction of elements between roughly Si and Ge are generally produced in CCSNe. We find that44Ti exhibits an extended production timescale as compared to56Ni, which may explain its different distribution and higher than previously predicted abundances in supernova remnants such as Cas A and SN1987A. We also discuss the morphology of the ejected elements. This study highlights the high-level diversity of ejecta conditions and nucleosynthetic results in 3D CCSN simulations and emphasizes the need for additional long-term (∼10 s) 3D simulations to properly address such complexities.

     
    more » « less
  5. Abstract

    We present numerical results from a parameter study of the standing accretion shock instability (SASI), investigating the impact of general relativity (GR) on the dynamics. Using GR hydrodynamics with GR gravity, and nonrelativistic (NR) hydrodynamics with Newtonian gravity, in an idealized model setting, we vary the initial radius of the shock, and by varying its mass and radius in concert, the proto-neutron star compactness. We investigate four compactnesses expected in a post-bounce core-collapse supernova (CCSN). We find that GR leads to a longer SASI oscillation period, with ratios between the GR and NR cases as large as 1.29 for the highest-compactness suite. We also find that GR leads to a slower SASI growth rate, with ratios between the GR and NR cases as low as 0.47 for the highest-compactness suite. We discuss implications of our results for CCSN simulations.

     
    more » « less