Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “core-collapse supernovae”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Impacts of Multidimensional Progenitor Perturbations on Core-collapse Supernova Explosions

Numerical studies of core-collapse supernovae have demonstrated the importance of nonradial motions in precollapse progenitors on the explosion outcome. We use the Chimera neutrino radiation hydrodynamics code running seven two-dimensional simulations of 15 M⊙ progenitors with different progenitor structures introduced by different one- and two-dimensional precollapse stellar evolution environments to examine the impacts of stellar structure and nonspherical motion in the precollapse progenitor on the development of explosions. We compare the explosion evolution of these models in terms of shock dynamics, diagnostic energy, neutrino heating, accretion, explosion geometry, nuclear abundances, and turbulent convection. We also analyze how stochastic variation impacts our simulations. Contrary to results reported in prior studies examining the impacts of multidimensional progenitors, we observe similar shock revival times and explosion development in our simulations despite differences in initial compositions and structures. We find no discernible impact from the accretion of nonradial perturbations from a multi-D progenitor onto the stalled shock in the revival and strength of explosion, as fully developed neutrino-driven convection behind the stalled shock is similar for all our models. For models with physically sourced noise in the iron core, a strong oscillation of the shock occurs after bounce and deflects infall laterally, and accelerates the saturation of the lateral turbulent kinetic energy. An examination of model stochasticity shows that any prior expected impacts on explosive outcome due to convection-related perturbations lie below the detectable threshold of numerical variation.

Chen, Chien-Hui [North Carolina State University]↗
DOE OSTI2026

Selected 3D Flash-X Checkpoints for the Long-Time Evolution of a 9.6 Solar-Mass Core-Collapse Supernova Model

This dataset contains selected 3D Flash-X checkpoint files from the long-time evolution of a low-energy core-collapse supernova explosion of a 9.6 Msun zero-metallicity, low-mass iron-core progenitor. The checkpoints span the shock-breakout phase through the young-remnant phase, ending at approximately 3 yr after core bounce. The dataset is intended for follow-up analysis and post-processing, especially radiation-transport calculations using the hydrodynamic and compositional structure of the ejecta. For a full description of the numerical setup, physical assumptions, limitations, and interpretation of the simulation, users should refer to the associated paper.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Production of Heavy α -Elements and 44 Ti in Cas A: Comparison to Abundances from 1D Core-collapse Supernova Models and Evidence for Carbon–Oxygen Shell Mergers

The merger between the carbon (C) and oxygen (O) shells hours to days before the collapse of a massive star significantly changes its nucleosynthesis, which is reflected in the elemental ratios observed in supernova remnants (SNRs). We present a nucleosynthesis study of 44 Ti production in core-collapse supernovae (CCSNe), highlighting large silicon (Si), sulfur (S), calcium (Ca), and, most importantly, argon (Ar) to neon (Ne) ratios as diagnostics for carbon–oxygen (C–O) shell mergers. We compare yields from eight different sets of CCSNe models to observations of Cassiopeia A (Cas A), and show that C–O shell mergers are consistently the models that best match X-ray and infrared observations. These models produce high Ar/Ne ratios (≳0.1), due to 20 Ne depletion and production of 36 Ar and 38 Ar, while lower ratios are obtained from nonmerger cases. Based on the Ar/Ne diagnostic, we compare the range of expected 44 Ti produced by C–O shell mergers, which is up to ∼20%–30% of the overall 44 Ti, but expected to be located outside the reverse shock. Based on the sets of models considered, the photon flux expected from the 44 Ti synthesized in the C–O shell merger in Cas A is below the NuSTAR and COSI detection limits, compatible with current limits locating most of the 44 Ti interior to the reverse shock, but might be detectable from proposed missions like ASCENT. Finally, for the SNR of 1987A, a dominant C–O merger origin of the observed 44 Ti is unlikely based on the observed redshift in its 44 Ti line.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

An Exploration of the Equation of State Dependence of Core-collapse Supernova Explosion Outcomes and Signatures

We explore, using a state-of-the-art simulation code in 3D and late-enough times to witness final observables, the dependence of core-collapse supernova explosions on the nuclear equation of state (EOS). Going beyond questions of explodability, we compare final explosion energies, nucleosynthetic yields, recoil kicks, and gravitational-wave and neutrino signatures using the SFHo and DD2 nuclear EOSs for a 9 M ⊙ /solar-metallicity progenitor star. The DD2 EOS is stiffer and has a lower effective nucleon mass. The result is a more extended protoneutron star (PNS) and lower central densities. As a consequence, the mean neutrino energies, final explosion energy, and recoil kick speed are lower. Moreover, the evolution of PNS convection differs between the two EOS models in significant ways. This translates in part into interestingly altered neutrino “light” curves and noticeably altered gravitational-wave signal strengths and frequency characteristics that may be diagnostic. The faster exploding model (SFHo) yields slightly more neutron-rich ejecta and more species with atomic weights between 60 and 90 and a weak r-process. However, this is merely a preliminary study. The next step is a more comprehensive and multiprogenitor set of 3D supernova simulations for various EOSs to late times when the observables have asymptoted. Such a future investigation will have a direct bearing on the neutron star and black hole birth mass functions and the quest toward a fully quantitative theory of supernova observables.

Rusakov, Aleksandr [Princeton University, NJ (Unit↗
DOE OSTI2026

thornado+FLASH-X: A Hybrid Discontinuous Galerkin–Implicit-explicit and Finite-volume Framework for Neutrino-radiation Hydrodynamics in Core-collapse Supernovae

We present neutrino-transport algorithms implemented in the toolkit for high-order neutrino-radiation hydrodynamics (thornado) and their coupling to self-gravitating hydrodynamics within the adaptive mesh refinement–based multiphysics simulation framework FLASH-X. thornado, developed primarily for simulations of core-collapse supernovae (CCSNe), employs a spectral, six-species two-moment formulation with algebraic closure and special-relativistic observer corrections accurate to $\mathcal{O}(v/c)$, and uses discontinuous Galerkin (DG) methods for phase-space discretization combined with implicit-explicit time stepping. A key development is a nonlinear neutrino–matter coupling algorithm based on nested fixed-point iteration with Anderson acceleration, enabling fully implicit treatment of collisional processes, including energy-coupling interactions such as neutrino–electron scattering and pair production. Coupling to finite-volume (FV) hydrodynamics is achieved through a hybrid DG-FV representation of the fluid variables and operator-split evolution within FLASH-X. The implementation is verified using basic transport tests with idealized opacities and relaxation and deleptonization problems with tabulated microphysics. Spherically symmetric CCSN simulations demonstrate accuracy and robustness of the coupled scheme, including close agreement with the CCSN simulation code Chimera. An axisymmetric CCSN simulation further demonstrates the viability of DG-based neutrino transport for multidimensional supernova modeling within FLASH-X. thornado’s neutrino-transport solver is GPU-enabled using OpenMP offloading or OpenACC, and all CCSN applications included in this work use the GPU implementation. Together, these results establish a foundation for future enhancements in physics fidelity, numerical algorithms, and computational performance, for increasingly realistic large-scale CCSN simulations.

Endeve, Eirik [Oak Ridge National Laboratory (ORNL↗
DOE OSTI2026

Enhanced antineutrino emission from 𝛽 decay in core-collapse supernovae with self-consistent weak decay rates

Nuclear weak-interaction rates are known to exert a prominent effect in the late stages of stellar collapse. Despite their importance, most studies to date on core-collapse supernovae (CCSNe) have focused primarily on the effects of electron captures, neglecting 𝛽 decay contributions. In this work, we present the first CCSNe simulation incorporating global 𝛽 decay rates from a microscopic theory. These are enabled by a large-scale evaluation of both electron capture and 𝛽 decay rates, obtained self-consistently utilizing the relativistic energy density functional theory and finite-temperature quasiparticle random-phase approximation. Including 𝛽 decay leads to a dramatic enhancement of the pre-bounce antineutrino signal as the antineutrino emissivity increases by more than two orders of magnitude and the luminosity by a factor of 50 relative to thermal emission alone, while the average antineutrino energy increases by over 1 MeV. Furthermore, it is expected that these new rates could help us constrain the model uncertainties related to weak-interaction processes, improving the prediction of antineutrino signal during the final stages of stellar death.

Beta decay↗
DOE OSTI2026

The Kinetic-Energy–Momentum–Mass 5-Flux of a Baryon Fluid in Bargmann Spacetimes

A Bargmann spacetime is a constrained five-dimensional setting that, while introducing no new physical degrees of freedom beyond those of ordinary four-dimensional spacetime, permits Galilei physics to be expressed with a tensor formalism that respects the distinction between mass and energy while affording the conceptual and technical advantages of a spacetime metric. This framework offers a route to a strong-field ‘Galilei general relativity’ approximating the usual Poincaré general relativity introduced by Einstein. In preparation for modeling core-collapse supernovae, where such an approximation would be useful, this work generalizes the kinetic-energy–momentum–mass 5-flux 𝒯 and its associated spacetime tensor law from a simple fluid of constant particle mass to a baryon fluid whose multiple nuclear species can interconvert rest mass and internal energy. The spacetime tensor law on Bargmann–Galilei spacetime 𝐵𝒢 and its decompositions relative to comoving (‘Lagrangian’) and fiducial (‘Eulerian’) observers are derived in detail. The formalism is rendered more suitable for core-collapse supernova modeling by an extension from strict 𝐵𝒢 to a regime that might be denoted as 𝐵𝒢+: microscopically Poincaré yet macroscopically Galilei. This extension accommodates energy generation by nuclear composition changes and allows comoving energy density and pressure to contribute relative to mass density, while preserving the simplifications of Galilei bulk fluid flow and the streamlined geometry governed by the Bargmann–Galilei spacetime metric.

Cardall, Christian [ORNL] (ORCID:000000020086105X)↗
DOE OSTI2026

Study of Key 57 Ni ⁢(𝑝,𝛾)⁢ 58 Cu Resonances to Understand 44 Ti Nucleosynthesis in Supernovae

An important validation of nucleosynthesis models of core-collapse supernovae is the comparison of radioisotope predictions to abundances inferred from observations of 𝛾 rays emitted in remnants. One such isotope, 44 Ti, is especially sensitive to the 57 Ni ⁢(𝑝,𝛾)⁢ 58 Cu reaction rate. Despite this importance, no experimentally constrained rate exists for this reaction. It is thus crucial to determine this rate. GRETINA ORRUBA: Dual Detectors for Experimental Structure Studies (GODDESS) was used along with the first use of the Enge split-pole spectrograph at Notre Dame to measure structure properties of 58 Cu via the 58 Ni ⁢( 3 He,𝑡)⁢ 58 Cu reaction. Here, the combined analysis of these complementary data sets allows for precise determination of 58 Cu level energies, and spins are constrained from analysis of the 𝛾-ray decay. 44 Ti yields were found to change by over 25% compared to previous estimates in model calculations using the new experimentally constrained reaction rate.

energy levels↗
DOE OSTI2026

A Search for Successful and Choked Jets in Nearby Broad-lined Type Ic Supernovae

The observational link between long gamma-ray bursts (GRBs) and broad-lined stripped-envelope core-collapse supernovae (SNe Ic-BL) is well established. Significant progress has been made in constraining what fraction of SNe Ic-BL may power high- or low-luminosity GRBs when viewed at small off-axis angles. However, the GRB–SN connection still lacks a complete understanding in the broader context of massive-star evolution and explosion physics. Models predict a continuum of outcomes for the fastest ejecta, from choked to ultrarelativistic jets, and observations from radio to X-rays are key to probing these scenarios across a range of viewing angles and velocities. Here, we present results from a coordinated radio-to-X-ray campaign targeting nearby (z ≲ 0.1) SNe Ic-BL designed to explore this diversity. With eight new radio-monitored events and updated data for one previously observed SN, we further tighten constraints on the fraction of SNe Ic-BL as relativistic as SN 1998bw (GRB 980425). We identify SN 2024rjw as a new radio-loud event likely powered by strong interaction with circumstellar material, and add evidence supporting a similar interpretation for SN 2020jqm. We also establish new limits on the properties of radio-emitting ejecta with velocities consistent with cocoons from choked jets, highlighting SN 2022xxf as a promising cocoon-dominated candidate. These results refine our understanding of the continuum linking ordinary SNe Ic-BL, engine-driven explosions, and GRBs, and contribute to building a sample that will inform future multimessenger searches for electromagnetic counterparts to high-energy neutrinos.

O’Dwyer, Tanner [Johns Hopkins Univ., Baltimore, M↗
DOE OSTI2026

Standard Candles for Supernova Neutrino Detection at DUNE

The Deep Underground Neutrino Experiment (DUNE) far detector is sensitive to $\mathcal{O}(10)$ MeV electron neutrinos through $ν_e$ charged-current reaction with argon. This capability is a unique window into the $ν_e$ component of a Galactic core-collapse supernova flux. Extracting the properties of the supernova spectrum is, however, limited by the poorly-known $ν_e$-Ar cross section. We propose a data-driven strategy that leverages $^8$B solar neutrinos and muon-decay-at-rest neutrinos as standard candles for this process. These calibration samples constrain both the low-energy and high-energy components of the cross section relevant for supernova detection. Our method reduces the reliance on nuclear models, which can bias the extraction of the spectral parameters by as much as 300$\%$.

Cheng, Ting [Fermilab]↗
DOE OSTI2026

Cosmic Overture Echoed in a Stellar Final Act: Implications of Nanohertz Gravitational Waves for Core Collapse Supernova Neutrinos

The growing evidence for nanohertz gravitational waves, from NANOGrav and other observations, may be pointing to a cosmological first-order phase transition at temperatures of 𝒪⁢(10−100)⁢MeV. Such an interpretation requires beyond the Standard Model dynamics in this energy range. If so, it may well be the case that galactic core-collapse supernova explosions would produce the key components related to the first-order phase transition, leaving detectable imprints on the spectrum of neutrinos emitted in the initial few seconds of the collapse. This might provide further evidence in support of the early universe interpretation of the nanohertz gravitational wave signal. The scenario proposed here is also suggestive of a low-mass seesaw mechanism to explain neutrino masses. We outline the prospects for future observations of Galactic supernovae to uncover the signals of this scenario, with further confirmation from future pulsar timing array measurements of primordial nanohertz gravitational waves.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

AT2025ulz and S250818k: Leveraging DESI Spectroscopy in the Hunt for a Kilonova Associated with a Subsolar-mass Gravitational-wave Candidate

On 2025 August 18, the LIGO–Virgo–KAGRA collaboration reported a subthreshold gravitational-wave candidate detection consistent with a subsolar-mass neutron star merger, denoted S250818k. An optical transient, AT2025ulz, was discovered within the localization region. AT2025ulz initially appeared to meet the expected behavior of kilonova emission, the telltale signature of a binary neutron star merger. The transient subsequently rebrightened after ∼5 days and was classified as a Type IIb supernova. In this work, we analyze the observations of its host galaxy obtained by the Dark Energy Spectroscopic Instrument (DESI). From the DESI spectrum, we obtain a secure redshift of z = 0.084840 ± 0.000006. If S250818k has an astrophysical origin, this places the transient within 2σ of the gravitational-wave distance and results in an integral overlap between the gravitational-wave alert and the transient location of log 10 $\mathcal{I} ≈ 3.9 - 4.2$. Our analysis of the host galaxy’s spectral energy distribution reveals a star-forming, dusty galaxy with stellar mass ∼10 10 M ⊙ , broadly consistent with the population of both short gamma-ray bursts and core-collapse supernova host galaxies. We also present our follow-up of DESI-selected candidate host galaxies using the Fraunhofer Telescope at the Wendelstein Observatory, and show the promise of DESI for associating or rejecting candidate electromagnetic counterparts to gravitational-wave alerts. These results emphasize the value of DESI’s extensive spectroscopic dataset in rapidly characterizing host galaxies, enabling spectroscopic host subtraction, and guiding targeted follow-up.

Hall, Xander J. [Carnegie Mellon University, Pitts↗
DOE OSTI2026

Continuum contribution to charged-current absorption of low-energy $ν_e$ on $^{40}$Ar

Accurate modeling of the absorption of tens-of-MeV $ν_e$ on $^{40}$Ar is needed to enable measurements of astrophysical neutrinos using large liquid argon time projection chamber (LArTPC) detectors, such as those planned for the Deep Underground Neutrino Experiment (DUNE). We revisit the MARLEY neutrino interaction model used in present estimates of DUNE sensitivity to supernova and solar neutrino signals. Multiple theoretical refinements are pursued, especially in the unbound continuum region of nuclear excitation energy. Inclusive charged-current neutrino-argon cross sections are calculated using a hybrid strategy. Nuclear transitions to unbound states are treated using a Hartree-Fock Continuum Random Phase Approximation (HF-CRPA) model, including forbidden contributions. Allowed transitions to low-lying discrete levels are also included using indirect measurements and approximate corrections for the momentum transfer dependence. Exclusive predictions are obtained by coupling these calculations with a statistical nuclear de-excitation model. The impact on observables of interest for DUNE and similar experiments is examined in terms of both total and differential cross sections. Our refined calculations predict a lower allowed portion of the cross section relative to the prior MARLEY model. At neutrino energies appreciably below 100 MeV, the inclusion of forbidden transitions does not fully compensate for the loss of allowed strength. For a representative neutrino burst from a galactic core-collapse supernova, our results suggest that MARLEY 1.2.0 overestimates the event yield in a DUNE-like detector by approximately 20%. However, because this overestimation is more severe at backwards angles, use of the charged-current $ν_e$-$^{40}$Ar reaction for supernova pointing may be more feasible than previously expected.

Gardiner, Steven [Fermilab]↗
DOE OSTI2026

Distributed Resilience in High-Energy Physics Data Acquisition

Historical experience in the High-Performance Computing community teaches us that as computing systems grow, the instance of failures goes from rare to a regular occurrence. A survey of the growth in the size and complexity of Data AcQuisition (DAQ) networks in High-Energy Physics (HEP) experiments reveals that these networks are scaling exponentially, trending to a point where automated fault handling should be considered over the current manual practice, especially given the rarity of data such as in DUNE's mission to observe core-collapse supernovae. We propose a general system, DiDAQt, which is designed to provide fault detection and handling in HEP DAQs specifically, through MPI-like primitives that allow it to be added easily to existing systems. We evaluate the scalability and response time of a prototype on the FABRIC national testbed, with results indicating sufficient scalability for current and near-future DAQs as well as practical response times (under 1 microsecond decision time).

Wolosewicz, A. [IIT, Chicago]↗
DOE OSTI2026

GenASiS: General Astrophysical Simulation System. II. Self-gravitating Baryonic Matter*

GenASiS (General Astrophysical Simulation System) is a code being developed initially and primarily, though not exclusively, for the simulation of core-collapse supernovae on the world's leading capability supercomputers. This paper---the second in a series---documents capabilities for Newtonian self-gravitating fluid dynamics, including tabulated microphysical equations of state treating nuclei and nuclear matter (`baryonic matter'). Computation of the gravitational potential of a spheroid, and simulation of the gravitational collapse of dust and of an ideal fluid, provide tests of self-gravitation against known solutions. In multidimensional computations of the adiabatic collapse, bounce, and explosion of spherically symmetric pre-supernova progenitors---which we propose become a standard benchmark for code comparisons---we find that the explosions are prompt and remain spherically symmetric (as expected), with an average shock expansion speed and total kinetic energy that are inversely correlated with the progenitor mass at the onset of collapse and the compactness parameter.

Cardall, Christian [ORNL] (ORCID:000000020086105X)↗
DOE OSTI2026

Pre-Equilibrium De-Excitations for Neutrino-Nucleus Interactions

The Deep Underground Neutrino Experiment (DUNE) is sensitive to MeV-scale energy depositions from low-energy astrophysical neutrinos, including those from core-collapse supernovae. Interpreting these detector signals requires accurate modeling of the nuclear de-excitation that follows from the neutrino-nucleus interaction. The MARLEY (Model of Argon Reaction Low-Energy Yields) event generator specializes in the low-energy regime. MARLEY currently assumes the residual nucleus equilibrates immediately after the primary interaction. This omits the intermediate pre-equilibrium stage in which energy redistributes among nucleons until statistical equilibrium is reached. While pre-equilibrium effects are well established for nucleon-induced reactions, they have not previously been studied for neutrino-nucleus interactions. This work addresses that gap by implementing a two-component exciton model, which is the first dedicated treatment of pre-equilibrium de-excitation for neutrino-nucleus interactions, restructured around MARLEY's existing class hierarchy to prepare for direct integration, including particle-hole state densities, internal transition rates, and pre-equilibrium particle emission. The calculations show encouraging agreement with the TALYS-2.2 nuclear reaction code for neutron-nucleus interactions. We further propose a concrete integration path into the full MARLEY event generator, including derived class structure and an extended event record for pre-equilibrium vertices in support of future reweighting. Remaining work focuses on refining the emission width calculation, adding $\gamma$-ray emission, and completing this integration to quantify the impact of pre-equilibrium effects on the expected low-energy neutrino signals in DUNE and similar experiments.

Visser, Erin [Michigan State U., East Lansing (mai↗