Supernova Simulations Research Articles

Influence of Effective Nucleon Mass on Equation of State for Supernova Simulations and Neutron Stars

Abstract We investigate the influence of the effective nucleon mass on the equation of state (EOS), which is constructed for simulations of core-collapse supernovae and binary neutron star mergers, within the relativistic mean-field (RMF) framework. The study introduces a new RMF parameter set, TM1m, which is a modification of the TM1e model with an adjusted effective mass, maintaining the saturation properties of nuclear matter. The TM1m model, with a larger effective mass ratio (M */M ~ 0.8) compared to the TM1e model (M */M ~ 0.63), is employed to construct a new EOS table, EOS5. This EOS table is designed to offer insights into the influence of the effective nucleon mass on the EOS within a relativistic framework, particularly above the saturation density. The results of EOS5 are compared with those obtained from other models, including both relativistic and nonrelativistic approaches. The properties of cold neutron stars, calculated using the TM1m model, are compatible with the existence of a 2 M ⊙ pulsar and the latest constraints on the tidal deformability and radii of a canonical 1.4 M ⊙ neutron star, derived from astrophysical observations.

Hyperons during proto-neutron star deleptonization and the emission of dark flavoured particles

Complementary to high-energy experimental efforts, indirect astrophysical searches of particles beyond the standard model have long been pursued. The present article follows the latter approach and considers, for the first time, the self-consistent treatment of the energy losses from dark flavoured particles produced in the decay of hyperons during a core-collapse supernova (CCSN). To this end, general relativistic supernova simulations in spherical symmetry are performed, featuring six-species Boltzmann neutrino transport, and covering the long-term evolution of the nascent remnant proto-neutron star (PNS) deleptonization for several tens of seconds. A well-calibrated hyperon equation of state (EOS) is therefore implemented into the supernova simulations and tested against the corresponding nucleonic model. It is found that supernova observables, such as the neutrino signal, are robustly insensitive to the appearance of hyperons for the simulation times considered in the present study. The presence of hyperons enables an additional channel for the appearance of dark sector particles, which is considered at the level of the Λ hyperon decay. Assuming massless particles that escape the PNS after being produced, these channels expedite the deleptonizing PNS and the cooling behaviour. This, in turn, shortens the neutrino emission timescale. The present study confirms the previously estimated upper limits on the corresponding branching ratios for low and high mass PNS, by effectively reducing the neutrino emission timescale by a factor of two. This is consistent with the classical argument deduced from the neutrino detection associated with SN1987A.

To Explode or Not Explode? An Analysis of the Explodability of Binary and Single Stars

Abstract Most massive stars in our universe are born in binary systems, and yet their fates are quite uncertain, due to a poor understanding of their evolution. It is well known, however, that the density profile of a collapsing core, particularly at the interface between silicon and oxygen shells, is important when determining the explodability of a star. While this has been studied for single stars, it has been rarely studied for binary stars. In this note, we explore the explodability of single and binary stars using self-consistent supernova simulations. We study the differences in the silicon–oxygen interface and compactness of the progenitor stars and how they affect the explodability. We find that binary-stripped stars with initial masses above 37 M ⊙ are not as explodable as their single-star counterparts.

Neutrino-driven core-collapse supernova yields in galactic chemical evolution

Abstract We provide yields from 189 neutrino-driven core-collapse supernova (CCSN) simulations covering zero-age main sequence masses between 11 and 75 M⊙ and three different metallicities. Our CCSN simulations have two main advantages compared to previous methods used for applications in Galactic chemical evolution (GCE). Firstly, the mass cut between remnant and ejecta evolves naturally. Secondly, the neutrino luminosities and thus the electron fraction are not modified. Both is key to obtain an accurate nucleosynthesis. We follow the composition with an in-situ nuclear reaction network including the 16 most abundant isotopes and use the yields as input in a GCE model of the Milky Way. We adopt a GCE which takes into account infall of gas as well as nucleosynthesis from a large variety of stellar sources. The GCE model is calibrated to reproduce the main features of the solar vicinity. For the CCSN models, we use different calibrations and propagate the uncertainty. We find a big impact of the CCSN yields on our GCE predictions. We compare the abundance ratios of C, O, Ne, Mg, Si, S, Ar, Ca, Ti, and Cr with respect to Fe to an observational data set as homogeneous as possible. From this, we conclude that at least half of the massive stars have to explode to match the observed abundance ratios. If the explosions are too energetic, the high amount of iron will suppress the abundance ratios. With this, we demonstrate how GCE models can be used to constrain the evolution and death of massive stars.

Detecting Gravitational Wave Memory in the Next Galactic Core-Collapse Supernova.

We present an approach to detecting (linear) gravitational wave memory in a Galactic core-collapse supernova using current interferometers. Gravitational wave memory is an important prediction of general relativity that has yet to be confirmed. Our approach uses a combination of linear prediction filtering and matched filtering. We present the results of our approach on data from core-collapse supernova simulations that span a range of progenitor mass and metallicity. We are able to detect gravitational wave memory out to 10kpc. We also present the false alarm probabilities assuming an on-source window compatible with the presence of a neutrino detection.

Insights into the Production of 44Ti and Nickel Isotopes in Core-collapse Supernovae

We report nucleosynthetic results for both 44Ti and nickel isotopes for 18 three-dimensional (3D) core-collapse supernova (CCSN) simulations extended to ∼20 s after bounce. We find that many of our long-term models are able to achieve 44Ti/56Ni ratios similar to that observed in Cassiopeia A, and modern supernova models can synthesize up to 2 × 10−4 M ⊙ of 44Ti. Neutrino-driven winds and the fact that there can be simultaneous accretion and explosion in 3D models of CCSNe play central roles in its production. We conclude that the 44Ti underproduction problem in previous CCSN models is no longer an issue. In addition, we discuss the production of both 57Ni and stable nickel/iron ratios and compare our results to observations of SN 1987A and the Crab.

Estimate for the Bulk Viscosity of Strongly Coupled Quark Matter Using Perturbative QCD and Holography.

Modern hydrodynamic simulations of core-collapse supernovae and neutron-star mergers require knowledge not only of the equilibrium properties of strongly interacting matter, but also of the system's response to perturbations, encoded in various transport coefficients. Using perturbative and holographic tools, we derive here an improved weak-coupling and a new strong-coupling result for the most important transport coefficient of unpaired quark matter, its bulk viscosity. These results are combined in a simple analytic pocket formula for the quantity that is rooted in perturbative quantum chromodynamics at high densities but takes into account nonperturbative holographic input at neutron-star densities, where the system is strongly coupled. This expression can be used in the modeling of unpaired quark matter at astrophysically relevant temperatures and densities.

Nucleosynthesis in the Innermost Ejecta of Magnetorotational Supernova Explosions in Three Dimensions

Core-collapse supernova (CCSN) explosions powered by rotation and magnetic fields present an interesting astrophysical site for nucleosynthesis that potentially contributes to the production of r-process elements. Here we present yields of the innermost ejecta in 3D magnetorotational CCSN models simulated using the CoCoNuT-FMT code. Strong magnetic fields tap the rotational energy of the proto−neutron star and lead to earlier and more energetic (∼3 × 1051 erg) explosions than typical neutrino-driven CCSNe. Compared to a reference nonmagnetic model, the ejecta in the magnetorotational models have much more neutron-rich components with Y e down to ∼0.25. Our post-processing calculations with the reaction network SkyNet show significant production of weak r-process elements up to mass number ∼130. We find negligible differences in the synthesis of heavy elements between two magnetorotational models with different initial field strengths of 1010 and 1012 G, in accord with their similar explosion dynamics. The magnetorotational models produce about ∼0.19 and 0.14 M ☉ of radioactive 56Ni, on the low end of inferred hypernova nickel masses. The yields are publicly available at Zenodo (doi: 10.5281/zenodo.10578981) for comparison with stellar abundance patterns, inclusion in modeling galactic chemical evolution, and comparison with other yield calculations. Our results add to the yet-restricted corpus of nucleosynthesis yields from 3D magnetorotational supernova simulations and will help quantify yield uncertainties.

Supernova Explosions of the Lowest-mass Massive Star Progenitors

We here focus on the behavior of supernovae that technically explode in 1D (spherical symmetry). When simulated in 3D, however, the outcomes of representative progenitors of this class are quite different in almost all relevant quantities. In 3D, the explosion energies can be 2 to 10 times higher, and there are correspondingly large differences in the 56Ni yields. These differences between the 3D and 1D simulations reflect in part the relative delay to explosion of the latter and in the former the presence of protoneutron star convection that boosts the driving neutrino luminosities by as much as ∼50% at later times. In addition, we find that the ejecta in 3D models are more neutron-rich, resulting in significant weak r-process and 48Ca yields. Furthermore, we find that in 3D the core is an interesting, though subdominant, source of acoustic power. In summary, we find that though a model might be found theoretically to explode in 1D, one must perform supernova simulations in 3D to capture most of the associated observables. The differences between 1D and 3D models are just too large to ignore.

Diagnostics of 3D explosion asymmetries of stripped-envelope supernovae by nebular line profiles

ABSTRACT Understanding the explosion mechanism and hydrodynamic evolution of core-collapse supernovae (SNe) is a long-standing quest in astronomy. The asymmetries caused by the explosion are encoded into the line profiles which appear in the nebular phase of the SN evolution – with particularly clean imprints in He star explosions. Here, we carry out nine different supernova simulations of He-core progenitors, exploding them in 3D with parametrically varied neutrino luminosities using the prometheus-hotb code, hydrodynamically evolving the models to the homologous phase. We then compute nebular phase spectra with the 3D Non-Local Thermodynamic Equilibrium spectral synthesis code extrass (EXplosive TRAnsient Spectral Simulator). We study how line widths and shifts depend on progenitor mass, explosion energy, and viewing angle. We compare the predicted line profile properties against a large set of Type Ib observations, and discuss the degree to which current neutrino-driven explosions can match observationally inferred asymmetries. With self-consistent 3D modelling – circumventing the difficulties of representing $^{56}$Ni mixing and clumping accurately in 1D models – we find that neither low-mass He cores exploding with high energies nor high-mass cores exploding with low energies contribute to the Type Ib SN population. Models which have line profile widths in agreement with this population give sufficiently large centroid shifts for calcium emission lines. Calcium is more strongly affected by explosion asymmetries connected to the neutron star kicks than oxygen and magnesium. Lastly, we turn to the near-infrared spectra from our models to investigate the potential of using this regime to look for the presence of He in the nebular phase.

Structure Factors for Hot Neutron Matter from AbInitio Lattice Simulations with High-Fidelity Chiral Interactions.

We present the first abinitio lattice calculations of spin and density correlations in hot neutron matter using high-fidelity interactions at next-to-next-to-next-to-leading order in chiral effective field theory. These correlations have a large impact on neutrino heating and shock revival in core-collapse supernovae and are encapsulated in functions called structure factors. Unfortunately, calculations of structure factors using high-fidelity chiral interactions were well out of reach using existing computational methods. In this Letter, we solve the problem using a computational approach called the rank-one operator (RO) method. The RO method is a general technique with broad applications to simulations of fermionic many-body systems. It solves the problem of exponential scaling of computational effort when using perturbation theory for higher-body operators and higher-order corrections. Using the RO method, we compute the vector and axial static structure factors for hot neutron matter as a function of temperature and density. The abinitio lattice results are in good agreement with virial expansion calculations at low densities but are more reliable at higher densities. Random phase approximation codes used to estimate neutrino opacity in core-collapse supernovae simulations can now be calibrated with abinitio lattice calculations.

Chaos in inhomogeneous neutrino fast flavor instability

In dense neutrino gases, the neutrino-neutrino coherent forward scattering gives rise to a complex flavor oscillation phenomenon not fully incorporated in simulations of neutron star mergers (NSM) and core collapse supernovae (CCSNe). Moreover, it has been proposed to be chaotic, potentially limiting our ability to predict neutrino flavor transformations in simulations. To address this issue, we explore how small flavor perturbations evolve in the nonlinear regime of the neutrino quantum kinetic equation within a narrow centimeter-scale region inside a NSM and a toy neutrino distribution. Our findings reveal that paths in the flavor state space of solutions with similar initial conditions diverge exponentially, exhibiting chaos. This inherent chaos makes the microscopic scales of neutrino flavor transformations unpredictable. However, the domain-averaged neutrino density matrix remains relatively stable, with chaos minimally affecting it. This particular property suggests that domain-averaged quantities remain reliable despite the exponential amplification of errors. Published by the American Physical Society 2024

Sensitivity of Simulations of Double-detonation Type Ia Supernovae to Integration Methodology

We study the coupling of hydrodynamics and reactions in simulations of the double-detonation model for Type Ia supernovae. When assessing the convergence of simulations, the focus is usually on spatial resolution; however, the method of coupling the physics together as well as the tolerances used in integrating a reaction network also play an important role. In this paper, we explore how the choices made in both coupling and integrating the reaction portion of a simulation (operator/Strang splitting versus the simplified spectral deferred corrections method we introduced previously) influences the accuracy, efficiency, and nucleosynthesis of simulations of double detonations. We find no need to limit reaction rates or reduce the simulation time step to the reaction timescale. The entire simulation methodology used here is GPU-accelerated and made freely available as part of the Castro simulation code.

Contribution of Neutrino-dominated Accretion Flows to the Cosmic MeV Neutrino Background

Neutrino-dominated accretion flows (NDAFs) are one of the important MeV neutrino sources and significantly contribute to the cosmic diffuse neutrino background. In this paper, we investigate the spectrum of the diffuse NDAF neutrino background (DNNB) by fully considering the effects of the progenitor properties and initial explosion energies based on core-collapse supernova (CCSN) simulations, and estimate the detectable event rate by the Super-Kamiokande detector. We find that the predicted background neutrino flux is mainly determined by the typical CCSN initial explosion energy and progenitor metallicity. For the optimistic cases, in which the typical initial explosion energy is low, the diffuse flux of the DNNB is comparable to the diffuse supernova neutrino background, which might be detected by upcoming larger neutrino detectors, such as Hyper-Kamiokande, the Jiangmen Underground Neutrino Observatory, and the Deep Underground Neutrino Experiment. Moreover, the strong outflows from NDAFs could dramatically decrease their contribution to the neutrino background.

Gray two-moment neutrino transport: Comprehensive tests and improvements for supernova simulations

In this work we extended an energy-integrated neutrino transport method to facilitate efficient, yet precise, modeling of compact astrophysical objects. We particularly focus on core-collapse supernovae. We implemented a gray neutrino-transport framework from the literature into FLASH and performed a detailed evaluation of its accuracy in core-collapse supernova simulations. Based on comparisons with results from simulations using energy-dependent neutrino transport, we incorporated several improvements to the original scheme. Our analysis shows that our gray neutrino transport method successfully reproduces key aspects from more complex energy-dependent transport across a variety of progenitors and equations of state. We find both qualitative and reasonable quantitative agreement with multi-group M1 transport simulations. However, the gray scheme tends to slightly favor shock revival. In terms of gravitational wave and neutrino signals, there is a good alignment with the energy-dependent transport, although we find 15-30<!PCT!> discrepancies in the average energy and luminosity of heavy-lepton neutrinos. Simulations using the gray transport are around four times faster than those using energy-dependent transport.

General-relativistic Radiation Transport Scheme in Gmunu. II. Implementation of Novel Microphysical Library for Neutrino Radiation—Weakhub

We introduce Weakhub, a novel neutrino microphysics library that provides opacities and kernels beyond conventional interactions used in the literature. This library includes neutrino–matter, neutrino–neutrino interactions and plasma process, along with corresponding weak and strong corrections. A full kinematics approach is adopted for the calculations of β-processes, incorporating various weak corrections and medium modifications due to the nuclear equation of state. Calculations of plasma processes, electron neutrino–antineutrino annihilation, and nuclear de-excitation are also included. We also present the detailed derivations of weak interactions and the coupling to the two-moment based general-relativistic multigroup radiation transport in the general-relativistic multigrid numerical (Gmunu) code. We compare the neutrino opacity spectra for all interactions and estimate their contributions at hydrodynamical points in core-collapse supernovae and binary neutron star (BNS) postmerger remnants, and predict the effects of improved opacities in comparison to conventional ones for a BNS postmerger at a specific hydrodynamical point. We test the implementation of the conventional set of interactions by comparing it to an open-source neutrino library NuLib in a core-collapse supernova simulation. We demonstrate good agreement with discrepancies of less than ∼10% in luminosity for all neutrino species, while also highlighting the reasons contributing to the differences. To compare the advanced interactions to the conventional set in core-collapse supernova modeling, we perform simulations to analyze their impacts on neutrino signatures, hydrodynamical behaviors, and shock dynamics, showing significant deviations.

Ookami: An A64FX Computing Resource

We present a look at Ookami, a project providing community access to a testbed supercomputer with the ARM-based A64FX processors developed by a collaboration between RIKEN and Fujitsu and deployed in the Japanese supercomputer Fugaku. We provide an overview of the project and details of the hardware, and describe the user base and education/training program. We present highlights from previous performance studies of two astrophysical simulation codes and present a strong scaling study of a full 3D supernova simulation as an example of the the machine’s capability.

Constraining the Onset Density for the QCD Phase Transition with the Neutrino Signal from Core-collapse Supernovae

The occurrence of a first-order hadron–quark matter phase transition at high baryon densities is investigated in astrophysical simulations of core-collapse supernovae, to decipher yet incompletely understood properties of the dense matter equation of state (EOS) using neutrinos from such cosmic events. It is found that the emission of a nonstandard second neutrino burst, dominated by electron antineutrinos, is not only a measurable signal for the appearance of deconfined quark matter but also reveals information about the state of matter at extreme conditions encountered at the supernova (SN) interior. To this end, a large set of spherically symmetric SN models is investigated, studying the dependence on the EOS and the stellar progenitor. General relativistic neutrino-radiation hydrodynamics is employed featuring three-flavor Boltzmann neutrino transport and a microscopic hadron-quark hybrid matter EOS class. Therefore, the DD2 relativistic mean-field hadronic model is employed, and several variations of it, and the string-flip model for the description of deconfined quark matter. The resulting hybrid model covers a representative range of onset densities for the phase transition and latent heats. This facilitates the direct connection between intrinsic signatures of the neutrino signal and properties of the EOS. In particular, a set of linear relations has been found empirically. These potentially provide a constraint for the onset density of a possible QCD phase transition from the future neutrino observation of the next galactic core-collapse SN, if a millisecond electron anti-neutrino burst is present around or less than 1 s.

The Physics of Core-Collapse Supernovae: Explosion Mechanism and Explosive Nucleosynthesis

Recent developments in multi-dimensional simulations of core-collapse supernovae have considerably improved our understanding of this complex phenomenon. In addition to that, one-dimensional (1D) studies have been employed to study the explosion mechanism and its causal connection to the pre-collapse structure of the star, as well as to explore the vast parameter space of supernovae. Nonetheless, many uncertainties still affect the late stages of the evolution of massive stars, their collapse, and the subsequent shock propagation. In this review, we will briefly summarize the state-of-the-art of both 1D and 3D simulations and how they can be employed to study the evolution of massive stars, supernova explosions, and shock propagation, focusing on the uncertainties that affect each of these phases. Finally, we will illustrate the typical nucleosynthesis products that emerge from the explosion.

Physical Correlations and Predictions Emerging from Modern Core-collapse Supernova Theory

In this paper, we derive correlations between core-collapse supernova observables and progenitor core structures that emerge from our suite of 20 state-of-the-art 3D core-collapse supernova simulations carried to late times. This is the largest such collection of 3D supernova models ever generated and allows one to witness and derive testable patterns that might otherwise be obscured when studying one or a few models in isolation. From this panoramic perspective, we have discovered correlations between explosion energy, neutron star gravitational birth masses, 56Ni and α-rich freezeout yields, and pulsar kicks and theoretically important correlations with the compactness parameter of progenitor structure. We find a correlation between explosion energy and progenitor mantle binding energy, suggesting that such explosions are self-regulating. We also find a testable correlation between explosion energy and measures of explosion asymmetry, such as the ejecta energy and mass dipoles. While the correlations between two observables are roughly independent of the progenitor zero-age main-sequence (ZAMS) mass, the many correlations we derive with compactness cannot unambiguously be tied to a particular progenitor ZAMS mass. This relationship depends on the compactness/ZAMS mass mapping associated with the massive star progenitor models employed. Therefore, our derived correlations between compactness and observables may be more robust than with ZAMS mass but can nevertheless be used in the future once massive star modeling has converged.