Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “supernova”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Explosion Mechanism of Core-collapse Supernovae: Role of the Si/Si–O Interface

Abstract We present a simple criterion to predict the explodability of massive stars based on the density and entropy profiles before collapse. If a pronounced density jump is present near the Si/Si–O interface, the star will likely explode. We develop a quantitative criterion by using ∼1300 1D simulations where ν -driven turbulence is included via time-dependent mixing-length theory. This criterion correctly identifies the outcome of the supernova more than 90% of the time. We also find no difference in how this criterion performs on two different sets of progenitors, evolved using two different stellar evolution codes: FRANEC and KEPLER. The explodability as a function of mass of the two sets of progenitors is very different, showing: (i) that uncertainties in the stellar evolution prescriptions influence the predictions of supernova explosions; (ii) the most important properties of the pre-collapse progenitor that influence the explodability are its density and entropy profiles. We highlight the importance that ν -driven turbulence plays in the explosion by comparing our results to previous works.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Tracer Particles for Core-collapse Supernova Nucleosynthesis: The Advantages of Moving Backward

Abstract After decades, the theoretical study of core-collapse supernova explosions is moving from parameterized, spherically symmetric models to increasingly realistic multidimensional simulations. However, obtaining nucleosynthesis yields based on such multidimensional core-collapse supernova simulations is not straightforward. Frequently, tracer particles are employed. Tracer particles may be tracked in situ during the simulation, but often they are reconstructed in a post-processing step based on the information saved during the hydrodynamic simulation. Reconstruction can be done in a number of ways, and here we compare the approaches of backward and forward integration of the equations of motion to the results based on inline particle trajectories. We find that both methods agree reasonably well with the inline results for isotopes for which a large number of particles contribute. However, for rarer isotopes that are produced only by a small number of particle trajectories, deviations can be large. For our setup, we find that backward integration leads to better agreement with the inline particles by more accurately reproducing the conditions following freeze-out from nuclear statistical equilibrium, because the establishment of nuclear statistical equilibrium erases the need for detailed trajectories at earlier times. Based on our results, if inline tracers are unavailable, we recommend backward reconstruction to the point when nuclear statistical equilibrium was last applied, with an interval between simulation snapshots of at most 1 ms for nucleosynthesis post-processing.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Multidimensional Radiation Hydrodynamics Simulations of Pulsational Pair-instability Supernovae

Stars with masses of 80–130 M ⊙ can encounter pulsational pair-instability at the end of their lives, which triggers consecutive episodes of explosive burning that eject multiple massive shells. Collisions between these shells produce bright transients known as pulsational pair-instability supernovae (PPI SNe) that may explain some extreme supernovae. In this paper, we present the first 2D and 3D radiation hydrodynamics simulations of PPI SNe with the CASTRO code. Radiative cooling causes the collided shells to evolve into thin, dense structures with hot spots that can enhance the peak luminosity of the SN by factors of 2–3. The light curve peaks at 1.9–2.1 × 10 43 erg s -1 for 50 days and then plateaus at 2–3 × 10 42 erg s -1 for 200 days, depending on the viewing angle. The presence of 12 C and 16 O and the absence of 28 Si and 56 Fe in its spectra can uniquely identify this transient as a PPI SN in follow-up observations. Our models suggest that multidimensional radiation hydrodynamics is required to model the evolution and light curves of all shell-collision SNe, such as Type IIne, not just PPI SNe.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Fermi-LAT Detection of the Supernova Remnant G312.4-0.4 in the Vicinity of 4FGL J1409.1-6121e

Abstract Gamma-ray emission provides constraints on the nonthermal radiation processes at play in astrophysical particle accelerators. This allows both the nature of accelerated particles and the maximum energy that they can reach to be determined. Notably, it remains an open question to what extent supernova remnants contribute to the sea of Galactic cosmic rays. In the Galactic plane, at around 312° of Galactic longitude, Fermi-LAT observations show an extended source (4FGL J1409.1−6121e) around five powerful pulsars. This source is described by one large disk of 0.°7 radius with a high significance of 45σin the 4FGL-DR3 catalog. Using 14 yr of Fermi-LAT observations, we revisited this region with a detailed spectro-morphological analysis in order to disentangle its underlying structure. Three sources have been distinguished, including the supernova remnant G312.4−0.4 whose gamma-ray emission correlates well with the shell observed at radio energies. The hard spectrum detected by the LAT, extending up to 100 GeV without any sign of a cutoff, is well reproduced by a purely hadronic model.

Astronomy & Astrophysics↗
DOE OSTI2024

Artificial Intelligence Assisted Inversion (AIAI): Quantifying the Spectral Features of 56 Ni of Type Ia Supernovae

Abstract Following our previous study of Artificial Intelligence Assisted Inversion (AIAI) of supernova analyses, we train a set of deep neural networks based on the 1D radiative transfer code TARDIS to simulate the optical spectra of Type Ia supernovae (SNe Ia) between 10 and 40 days after the explosion. The neural networks are applied to derive the mass of 56 Ni in velocity ranges above the photosphere for a sample of 124 well-observed SNe Ia in the TARDIS model context. A subset of the SNe have multi-epoch observations for which the decay of the radioactive 56 Ni can be used to test the AIAI quantitatively. The 56 Ni mass derived from AIAI using the observed spectra as inputs for this subset agrees with the radioactive decay rate of 56 Ni. AIAI reveals that a spectral signature near 3890 Å is related to the Niii4067Å line, and the 56 Ni mass deduced from AIAI is found to be correlated with the light-curve shapes of SNe Ia, with SNe Ia with broader light curves showing larger 56 Ni mass in the envelope above the photosphere. AIAI enables spectral data of SNe to be quantitatively analyzed under theoretical frameworks based on well-defined physical assumptions.

Astronomy & Astrophysics↗
DOE OSTI2024

Nucleosynthetic Analysis of Three-dimensional Core-collapse Supernova Simulations

Abstract We study in detail the ejecta conditions and theoretical nucleosynthetic results for 18 three-dimensional core-collapse supernova (CCSN) simulations done by F ornax . 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 that 44 Ti exhibits an extended production timescale as compared to 56 Ni, 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.

Astronomy & Astrophysics↗
DOE OSTI2024

Evaluating the Consistency of Cosmological Distances Using Supernova Siblings in the Near-infrared

Abstract The study of supernova (SN) siblings, supernovae with the same host galaxy, is an important avenue for understanding and measuring the properties of Type Ia SN Ia light curves (LCs). Thus far, sibling analyses have mainly focused on optical LC data. Considering that LCs in the near-infrared (NIR) are expected to be better standard candles than those in the optical, we carry out the first analysis compiling SN siblings with only NIR data. We perform an extensive literature search of all SN siblings and find six sets of siblings with published NIR photometry. We calibrate each set of siblings ensuring they are on homogeneous photometric systems, fit the LCs with the SALT3-NIR and SNooPy models, and find median absolute differences inμvalues between siblings of 0.248 and 0.186 mag, respectively. To evaluate the significance of these differences beyond measurement noise, we run simulations that mimic these LCs and provide an estimate for uncertainty on these median absolute differences of ∼0.052 mag, and we find that, statistically, our analysis rules out the nonexistence of intrinsic scatter in the NIR at the 99% level. When comparing the same sets of SN siblings, we observe a median absolute difference inμvalues between siblings of 0.177 mag when using optical data alone as compared to 0.186 mag when using NIR data alone. It is unclear if these results may be due to limited statistics or poor quality NIR data, all of which will be improved with the Nancy Grace Roman Space Telescope.

Astronomy & Astrophysics↗
DOE OSTI2024

How Population III Supernovae Determined the Properties of the First Galaxies

Abstract Massive Population III stars can die as energetic supernovae that enrich the early Universe with metals and determine the properties of the first galaxies. With masses of about 10 9 M ⊙ atz≳ 10, these galaxies are believed to be the ancestors of the Milky Way. This paper investigates the impact of Population III supernova remnants (SNRs) from both Salpeter-like and top-heavy initial mass functions (IMFs) on the formation of first galaxies with high-resolution radiation-hydrodynamical simulations with theENZOcode. Our findings indicate that SNRs from a top-heavy Population III IMF produce more metals, leading to more efficient gas cooling and earlier Population II star formation in the first galaxies. From a few hundred to a few thousand Population II stars can form in the central regions of these galaxies. These stars have metallicities of 10 −3 –10 −2 ,Z ⊙ , greater than those of extremely metal-poor (EMP) stars. Their mass function follows a power-law distribution with dN ( M * ) / dM * ∝ M * α , whereM * is stellar mass, andα= 2.66–5.83 and is steeper for a top-heavy IMF. We thus find that EMP stars were not typical of most primitive galaxies.

Astronomy & Astrophysics↗
DOE OSTI2024

Extrapolation of Type Ia Supernova Spectra into the Near-infrared Using Principal Component Analysis

Abstract We present a method of extrapolating the spectroscopic behavior of Type Ia supernovae (SNe Ia) in the near-infrared (NIR) wavelength regime up to 2.30 μ m using optical spectroscopy. Such a process is useful for accurately estimating K-corrections and other photometric quantities of SNe Ia in the NIR. A principal component analysis is performed on data consisting of Carnegie Supernova Project I & II optical and NIR FIRE spectra to produce models capable of making these extrapolations. This method differs from previous spectral template methods by not parameterizing models strictly by photometric light-curve properties of SNe Ia, allowing for more flexibility of the resulting extrapolated NIR flux. A difference of around −3.1% to −2.7% in the total integrated NIR flux between these extrapolations and the observations is seen here for most test cases including Branch core-normal and shallow-silicon subtypes. However, larger deviations from the observation are found for other tests, likely due to the limited high-velocity and broad-line SNe Ia in the training sample. Maximum-light principal components are shown to allow for spectroscopic predictions of the color-stretch light-curve parameter, s BV , within approximately ±0.1 units of the value measured with photometry. We also show these results compare well with NIR templates, although in most cases the templates are marginally more fitting to observations, illustrating a need for more concurrent optical+NIR spectroscopic observations to truly understand the diversity of SNe Ia in the NIR.

Astronomy & Astrophysics↗
DOE OSTI2024

Supernova Explosions of the Lowest-mass Massive Star Progenitors

Abstract 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 56 Ni 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 48 Ca 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.

Astronomy & Astrophysics↗
DOE OSTI2024

Modeling the Progenitor Stars of Observed Type IIP Supernovae

Abstract Type IIP supernovae (SNe IIP) are thought to originate from the explosion of massive stars >10M ⊙ . Their luminosity is primarily powered by the explosion energy and the radioactive decay energy of 56 Co, with the photosphere location regulated by hydrogen recombination. However, the physical connections between SNe IIP and their progenitor stars remain unclear. This paper presents a comprehensive study of SNe IIP and their progenitor stars by using the one-dimensional stellar evolution code,MESA. Our model grids consider the effects of stellar metallicity, mass, and rotation in the evolution of massive stars, as well as the explosion energy and 56 Ni production in modeling supernovae. To elucidate the observed SNe IIP and their origins, we compare their light curves (LCs) with our models. Furthermore, we investigate the impact of stellar parameters on LCs by considering stellar mass, metallicity, rotation, explosion energy, and 56 Ni production. We find that more massive stars exhibit longer plateaus due to increased photon diffusion time caused by massive ejecta. Higher metallicity leads to increased opacity and mass loss of progenitor stars. Rapid rotation affects internal stellar structures, enhancing convective mixing and mass loss, potentially affecting the plateau’s brightness and duration. Higher explosion energy results in brighter but shorter plateaus due to faster-moving ejecta. 56 Ni mass affects late-time luminosity and plateau duration, with larger masses leading to slower declines.

Astronomy & Astrophysics↗
DOE OSTI2024

Gravitational-wave and Gravitational-wave Memory Signatures of Core-collapse Supernovae

Abstract In this paper, we calculate the energy, signal-to-noise ratio (SNR), detection range, and angular anisotropy of the matter, matter memory, and neutrino memory gravitational-wave (GW) signatures of 21 three-dimensional initially nonrotating core-collapse supernova (CCSN) models carried to late times. We find that inferred energy, SNR, and detection range are angle-dependent quantities, and that the spread of possible energy, signal to noise, and detection ranges across all viewing angles generally increases with progenitor mass. When examining the low-frequency matter memory and neutrino memory components of the signal, we find that the neutrino memory is the most detectable component of a CCSN GW signal, and that DECIGO is best equipped to detect both matter memory and neutrino memory. Moreover, we find that the polarization angle between the h + and h × strains serves as a unique identifier of matter and neutrino memory. Finally, we develop a Galactic density- and stellar mass-weighted formalism to calculate the rate at which we can expect to detect CCSN GW signals with the Advanced Laser Interferometer Gravitational-Wave Observatory (aLIGO). When considering only the matter component of the signal, the aLIGO detection rate is around 65% of the total Galactic supernova rate, but increases to 90% when incorporating the neutrino memory component. We find that all future detectors (Einstein Telescope, Cosmic Explorer, DECIGO) will be able to detect CCSN GW signals from the entire Galaxy, and for the higher-mass progenitors even into the Local Group of galaxies.

Choi, Lyla (ORCID:0000000153860133)↗
DOE OSTI2025

The Dark Energy Bedrock All-sky Supernova Program: Cross Calibration, Simulations, and Cosmology Forecasts

Type Ia supernovae (SNe Ia) have been essential for probing the nature of dark energy; however, most SN analyses rely on the same low-redshift sample, which may lead to shared systematics. In a companion paper, we introduce the Dark Energy Bedrock All-Sky Supernova (DEBASS) program, which has already collected more than 500 low-redshift SNe Ia on the Dark Energy Camera, and present an initial release of 77 SNe Ia within the Dark Energy Survey (DES) footprint observed between 2021 and 2024. Here, we examine the systematics, including photometric calibration and selection effects. We find agreement at the 10 mmag level among the tertiary standard stars of DEBASS, DES, and Pan-STARRS1. Our simulations reproduce the observed distributions of DEBASS SN light-curve properties, and we measure a bias-corrected Hubble residual scatter of 0.08 mag, which, while small, is found in 10% of our simulations. We compare the DEBASS SN distances to the Foundation sample and find consistency with a median residual offset of 0.016 ± 0.019 mag. Selection effects have negligible impacts on distances, but a different photometric calibration solution shifts the median residual −0.015 ± 0.019 mag, highlighting calibration sensitivity. Using conservative simulations, we forecast that replacing historical low-redshift samples with the full DEBASS sample will improve the statistical uncertainties on dark energy parameters w 0 and w a by 30% and 24%, respectively, enhance the dark energy Figure of Merit by up to 60%, and enable a measurement of fσ 8 at the 25% level.

Acevedo, Maria [Duke Univ., Durham, NC (United Sta↗
DOE OSTI2026

The Chicago Carnegie Hubble Program: Improving the Calibration of Type Ia Supernovae with JWST Measurements of the Tip of the Red Giant Branch

We present distances to 10 supernova (SN) host galaxies determined via the tip of the red giant branch using JWST/NIRCam and the F115W, F356W, and F444W bandpasses. The majority of the analysis was conducted on photometric catalogs that had their absolute zero-points randomized to mask information on distance. The new F115W TRGB distances, anchored by the geometric maser distance to NGC 4258, agree well with our previously derived Hubble Space Telescope (HST) TRGB distances, differing by only 1% on average and 4% on a per-galaxy basis. The color-corrected F115W TRGB is therefore equally precise a method of distance measurement as, and offers unique advantages over, its color-insensitive, I-band counterpart. We use these distances to update four published H 0 calibrations and evaluate how different SN analyses, both within and across independent groups, yield different H 0 values. For our JWST sample of 11 SNe, we find consistent values of H 0 ≃ 69 km s −1 Mpc −1 , with a dispersion of just 0.6 km s −1 Mpc −1 across the updated calibrations. When we expand the sample to 24 by combining with HST TRGB measurements, the results from different SN analyses begin to diverge, with the H 0 based on Pantheon+ and the Carnegie Supernova Project II (CSP-II), respectively, increasing by +2.0 km s −1 Mpc −1 (3.1σ significance) and +0.8 km s −1 Mpc −1 (1.4σ significance). More independent analyses of low-redshift SNe and JWST observations of the TRGB are needed to improve our understanding of systematics in distance ladder determinations of H 0 .

Hoyt, Taylor J. [Lawrence Berkeley National Labora↗
DOE OSTI2026

Neutrino and Electromagnetic Signatures from Superluminous Supernovae: A Case Study for SN 2017egm

Superluminous supernovae (SLSNe) are rare transients that are ∼10–100 times more luminous than ordinary stellar explosions, reaching peak optical luminosities ∼10 44 –10 45 erg s −1 . The energy source powering SLSNe remains uncertain. In this work, we explore the multiwavelength and multimessenger signatures of the scenario in which SLSNe are powered by a newly born millisecond magnetar. We model the dynamical evolution and emission from the coupled system comprised of the magnetar, wind, nebula, and supernova ejecta, consistently evaluating the pair multiplicity of the wind and nebula regions, and the bulk wind Lorentz factor governing the e + − e − injection spectra in the nebula. We compute the thermal and nonthermal electromagnetic signatures, neutrino signatures, and investigate their detection prospects. For SN 2017egm, the nearest observed SLSNe, our prediction for high-energy gamma rays matches the recent detection by the Fermi Large Area Telescope. For neutrinos, using SN 2017egm a canonical SLSNe, we find that in the era of the Vera C. Rubin Observatory, a stacking analysis with upcoming neutrino observatories can lead to 3σ detection significance of neutrino events from a population of SLSNe within a decade of operation.

Mukhopadhyay, Mainak [Fermi National Accelerator L↗
DOE OSTI2026

Simulation to a Newborn Supernova Remnant from a Low-mass Iron Core Star

Supernova remnant observations show a high degree of asymmetry, mixing, and inhomogeneity. These asymmetries are seeded during the early seconds of the explosion and are further enhanced and modified as the shock and ejecta move through the stellar progenitor and into the circumstellar medium. We present simulations of a 9.6 M⊙ zero-metallicity progenitor initialized after shock revival and evolved for several years when the ejecta is in the circumstellar medium. A suite of 1D and 2D simulations examines the effects of neutron-star wind and radioactive decay heating. In 1D, decay heating forms a low-density bubble that suppresses the reverse shock. While in 2D, the heating is localized to metal-rich pockets, inflating them and compressing the surrounding material into dense shells. In 3D, the neutron-star wind and decay heating modify the plume morphology, producing more large-scale structures. The extended plume morphology leads to an asymmetrical shock breakout. After breakout, the leading plumes cannot keep up with the shock front, resulting in deceleration and fragmentation by the reverse shock while retaining the large-scale asymmetry. The projected ejecta morphology and velocities are strongly viewing angle dependent. The relatively uniform metal-rich distribution does not resemble the strongly inhomogeneous ejecta structure of Cas A. The 160-isotope decay network shows that 24.4% of the radioactive heating comes from decay chains other than the canonical 56Ni chain. The low explosion energy, low 56Ni yield, and Ni/Fe ratio greater than unity suggest an observational signature similar to an electron capture supernova.

Neopane, Sudarshan [University of Tennessee (UT)]↗
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↗