Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Ejecta”

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 1,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
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 OSTI2025

Hydrodynamic Predictions for the Next Outburst of T Coronae Borealis: It Will Be the Brightest Classical or Recurrent Nova Ever Observed in X-Rays*

T Coronae Borealis (TCrB) is a recurrent nova with recorded outbursts in 1866 and 1946 and possible outbursts in 1217 and 1787. It is predicted to explode again in 2025 or 2026 based on multiple observational studies. The system consists of a massive ( M wd ≳ 1.35 M ⊙ ) white dwarf (WD) and a red giant (M3–M4 III). We have performed 1D hydrodynamic simulations with NOVA to predict the behavior of the next outburst. These simulations consist of a range of mass accretion rates onto ∼1.35 M ⊙ WDs, designed to bound the conditions necessary to achieve ignition of an explosion after an ≈80 yr interoutburst period. We have used both carbon–oxygen and oxygen–neon initial compositions, in order to include the possible ejecta abundances to be measured in the observations of the next outburst. As the WD in the TCrB system is observed to be massive, theoretical predictions reported here imply that the WD is growing in mass as a consequence of the thermonuclear runaway. Therefore, the secular evolution of the WD may allow it to approach the Chandrasekhar limit and either explode as a Type Ia supernova or undergo accretion-induced collapse, depending on its underlying composition. We have followed the evolution of just the WD, after removing the ejected matter from the surface layers. Our intent is to illuminate the mystery of the unique, second maximum in the two well-observed outbursts and we have found conditions that bracket the predictions.

Cataclysmic variable stars↗
DOE OSTI2025

Low-metallicity Nova Explosions: A Site for Weak rp -process Nucleosynthesis

Classical novae are common cataclysmic events involving a binary system of a white dwarf and a main-sequence or red giant companion star. In metal-poor environments, these explosions produce ejecta differently from their solar counterparts due to the accretion of subsolar-metallicity material onto the white dwarf. In particular, it has been suggested that the nucleosynthesis flow in such low-metallicity nova explosions extends up to the Cu–Zn region, much beyond the expected endpoint, around Ca, predicted for solar-metallicity classical novae. This behavior resembles a weak rp-process, and such nuclear activity has never been observed in accreting white dwarf binaries with typical accretion flows. In this work, we study the characteristics of the weak rp-process for four nova models with metallicities Z = 2 × 10 −9 , 10 −7 , 2 × 10 −6 , and 2 × 10 −5 , and explore the impact of the nuclear physics uncertainties via a Monte Carlo sensitivity study. We identify nuclear reactions whose uncertainties affect the production of intermediate-mass nuclei under these conditions. These reactions and relevant nuclear quantities are targets for measurements at stable or radioactive beam facilities to reduce their rate uncertainties.

Classical novae↗
DOE OSTI2025

Production of Radioactive 22 Na in Core-collapse Supernovae: The Ne-E(L) Component in Presolar Grains and Its Possible Consequences on Supernova Observations

Presolar graphite grains carry the isotopic signatures of their parent stars. A significant fraction of presolar graphites show isotopic abundance anomalies relative to solar for elements such as O, Si, Mg, and Ca, which are compatible with nucleosynthesis in core-collapse supernovae (CCSNe). Therefore, they must have condensed from CCSN ejecta before the formation of the Sun. Their most puzzling abundance signature is the 22 Ne-enriched component Ne-E(L), interpreted as the effect of the radioactive decay of 22 Na (T 1/2 = 2.6 yr). Previous works have shown that if H is ingested into the He shell and not fully destroyed before the explosion, the CCSN shock in the He-shell material produces large amounts of 22 Na. Here we focus on such CCSN models, showing a radioactive 26 Al production compatible with grain measurements, and analyze the conditions of 22Na nucleosynthesis. In these models, 22 Na is mostly made in the He shell, with a total ejected mass varying between 2.6 × 10 −3 M ⊙ and 1.9 × 10 −6 M ⊙ . We show that such 22 Na may already impact the CCSN light curve 500 days after the explosion, and at later stages it can be the main source powering the CCSN light curve for up to a few years before 44 Ti decay becomes dominant. Based on the CCSN yields above, the 1274.53 keV γ-ray flux due to 22 Na decay could be observable for years after the first CCSN light is detected, depending on the distance. This makes CCSNe possible sites to detect a 22 Na γ-ray signature consistently with the Ne-E(L) component found in presolar graphites. Finally, we discuss the potential contribution from 22 Na decay to the Galactic positron annihilation rate.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

An Integral-based Technique to Accelerate the Monte Carlo Radiative Transfer Computation for Supernovae

We present an integral-based technique (IBT) algorithm to accelerate supernova (SN) radiative transfer calculations. The algorithm utilizes “integral packets,” which are calculated by the path integral of the Monte Carlo (MC) energy packets, to synthesize the observed spectropolarimetric signal at a given viewing direction in a 3D time-dependent radiative transfer program. Compared to the event-based technique (EBT) proposed by M. Bulla et al., our algorithm significantly reduces the computation time and increases the MC signal-to-noise ratio (S/N). Using a 1D spherical symmetric Type Ia SN ejecta model DDC10 and its derived 3D model, the IBT algorithm has successfully passed the verification of spherical symmetry and cross comparison on a 3D SN model with the direct-counting technique and EBT. Notably, with our algorithm implemented in the 3D MC radiative transfer code SEDONA, the computation time is faster than EBT by a factor of 10−30, and the S/N is better by a factor of 1.5−3, with the same number of MC quanta.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Elemental and Isotopic Yields from T Coronae Borealis: Predictions and Uncertainties

T Coronae Borealis (T CrB) is a symbiotic recurrent nova system expected to undergo its next outburst within the next 2 yr. Recent hydrodynamic simulations have predicted the nucleosynthetic yields for both carbon–oxygen (CO) and oxygen–neon (ONe) white-dwarf models, but without accounting for thermonuclear reaction rate uncertainties. We perform detailed Monte Carlo postprocessing nucleosynthesis calculations based on updated thermonuclear reaction rates and uncertainties from the 2025 evaluation. We quantify the resulting abundance uncertainties and identify the key nuclear reactions that dominate them. Our results show that both the CO and ONe nova models robustly produce characteristic CNO isotopes. More pronounced abundance differences emerge for elements with A ≥ 20. Sulfur is the most robust observational discriminator between the CO and ONe nova models, with a model-to-model difference of a factor of ≈30 and minimal sensitivity to reaction rate uncertainties. Neon, silicon, and phosphorus exhibit even larger abundance differences (factors of ≈150–250), providing strong diagnostic potential. While their predicted yields are subject to larger uncertainties, these remain smaller than the model-to-model differences, allowing these elements to serve as useful, though less precise, tracers of white-dwarf composition. Chlorine, argon, and potassium also differ between models, but the 1σ-abundance ranges for the CO and ONe models overlap, reducing their present usefulness as composition tracers. We find that only nine nuclear reactions dominate the abundance uncertainties of the most diagnostically important isotopes, and their influence is largely independent of the underlying white-dwarf composition. These results provide guidance for future experimental efforts and for interpreting ejecta compositions in the next eruption of T CrB.

Chemical Abundances↗
DOE OSTI2026

Revisiting the Supernova Engines in the 3C 397 and W49B Supernova Remnants

The nature of the supernova remnants (SNRs) 3C 397 and W49B has long been a subject of debate, with prior studies offering conflicting interpretations between thermonuclear and core-collapse scenarios. To help settle this debate, we present a systematic, spatially resolved, spectroscopic analysis of both remnants using XMM-Newton. By applying multicomponent thermal models, we derive key physical properties including elemental abundances, ejecta temperatures, ambient densities, and explosion energetics. We compare the inferred metal abundance ratios to a wide range of core-collapse and thermonuclear nucleosynthesis models, including new models whose explosion energies differ from the canonical value of 10 51 erg. We find that the observed Fe/Si and Ca/Si ratios in both SNRs are best matched by certain thermonuclear models. However, no model fully reproduces the complete set of observed abundance patterns. In 3C 397, high Fe enrichment and spatial abundance variations suggest interaction with a dense progenitor environment, and W49B’s composition is overall consistent with a thermonuclear origin; however, both require a low-energy (∼10 50 erg) supernova explosion. We additionally map the Fe Kα line centroid energies and find a spread, with W49B falling within the core-collapse region—highlighting both environmental complexity and the limitations of this diagnostic for supernova classification. Our results highlight the need for caution in relying on any single diagnostic or nucleosynthesis model for supernova typing, underscore the need for improved nucleosynthesis models, and motivate future high-resolution, high-throughput observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

The Emission and Suppression of Line Features in Luminous Transients

Featureless optical and ultraviolet (UV) spectra are a puzzling signature to emerge from recent observations of luminous fast blue optical transients (LFBOTs) and some tidal disruption events (TDEs). We describe the landscape of source and gas properties that are expected to form H, He I , and He II emission lines, and map spectral types to the parameter space of luminosity and system radius. Using one-dimensional radiative transfer calculations, we show that high source luminosities (L > 10 44 erg s −1 ) and compact ejecta radii (r < 10 14 cm) produce featureless spectra, due to the high temperature and ionization state of the emitting medium. Intermediate luminosities and moderately compact systems can generate He II -dominated spectra, while lower luminosities and more extended atmospheres result in conspicuous H and He I emission. Large expansion velocities (v ≥ 0.1c) can further broaden lines such that they blend into the continuum. Featureless UV spectra may require even more extreme ionization environments or velocities in order to suppress the many intrinsically strong metal lines at those wavelengths. Applying this framework to understand the absence of features observed in LFBOTs and featureless TDEs, we find that, under the optically thick, quasi-thermal conditions considered here, nonhomologous, compact outflows are likely necessary for featurelessness to persist in optical and UV spectra.

Aspegren, Olivia [University of California, Berkel↗
DOE OSTI2026

Probing the γ -Ray Emission Origin of Two Star-forming Galaxies NGC 2403 and NGC 3424 with the Fermi-LAT

Star-forming galaxies (SFGs) are a subclass of γ-ray emitters, and a correlation between their γ-ray luminosity (L γ ) and the total infrared (IR) luminosity (L IR ) has been established based on the Fermi Large Area Telescope (LAT) data. NGC 2403 and NGC 3424 have been reported as outliers in the L γ –L IR correlation with light curves showing significant variability, which contrasts with the temporally stable γ-ray emission in other SFGs, originating primarily from cosmic rays interacting with interstellar medium. In this study, we reanalyze the γ-ray emission in the directions of NGC 2403 and NGC 3424 using more than 16.5 yr Fermi-LAT data. NGC 3424 is found to be spatially coincident with the detected γ-ray source, while NGC 2403 is significantly offset from the nearest γ-ray source, suggesting an implausible association. We confirm the previously reported variability of both γ-ray sources and the significant deviation from the L γ –L IR correlation when assuming an association of both γ-ray sources with the two galaxies. Our findings lend further support to the interpretation that their γ-ray emission is driven primarily by alternative radiative processes—rather than by star formation activity—such as the ejecta of the Type IIP supernova SN 2004dj in NGC 2403 interacting with a surrounding high-density shell and an obscured active galactic nucleus in NGC 3424.

Liu, Linjie [Chinese Academy of Sciences (CAS), Ku↗
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

Late-time Hubble Space Telescope Ultraviolet Spectra of SN 2023ixf and SN 2024ggi Show Ongoing Interaction with Circumstellar Material

We present far- and near-ultraviolet (UV) spectra of the type II supernovae (SNe) SN 2023ixf from days 199 to 722 and SN 2024ggi at days 41 and 232. Both SNe show broad, blueshifted, and asymmetric UV emission lines with an initial maximum velocity of ∼9000 km s −1 and narrow unresolved emission in C ivλλ1548.9, 1550.8. We compare the optical and UV emission-line profiles, showing that they evolve from two distinct velocity profiles to a single profile tracing the UV emission. We interpret this as shock power from interaction with circumstellar material coming to dominate over the radioactive-decay power from the inner ejecta. Comparing our observations to radiative transfer models with injected shock power, we find SN 2024ggi is best matched by P shock,abs = 1 × 10 41 erg s −1 at day 40; SN 2023ixf at day 300 and SN 2024ggi at day 200 are best matched by P shock,abs = 1 × 10 40 erg s −1 ; and SN 2023ixf at day 600 is best matched by Pshock,abs = 5 × 10 39 erg s −1 . From these models, we find that the mass-loss rate of both SNe increased just before the explosion. For SN 2023ixf, our mass-loss rates go from 4 × 10 −5 M ⊙ yr −1 at 600 yr before explosion to 2 × 10 −2 M ⊙ yr − 1 at 15 yr prior to explosion. For SN 2024ggi, we find a mass-loss rate of 9 × 10 −5 M ⊙ yr −1 at 150 yr before explosion and 1 × 10 −3 M ⊙ yr −1 at 30 yr before explosion.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Gamma-Ray Spectra of r -process Nuclei

The radioactive decay of unstable nuclei created in the rapid neutron capture process releases a large amount of γ-rays. When the ejecta are optically thick, these γ-rays may contribute to an associated kilonova. Once transparent, prominent spectral features will be directly observable in current and future γ-ray detectors. In this work, we study and compare the γ-ray spectra of different representative r-process trajectories across a broad range of timescales, from hours to 50,000 yr, identifying the nuclei that significantly contribute. We discuss these findings in the context of observability, noting that there are several practical challenges to connecting observed signatures to specific nuclei. However, if these challenges can be overcome, direct observation of γ-rays from r-process sites can provide insight into the fundamental physics underpinning the r-process.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
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 OSTI2023

The Nuclear Reaction Network WinNet

Abstract We present the state-of-the-art single-zone nuclear reaction network WinNet , which is capable of calculating the nucleosynthetic yields of a large variety of astrophysical environments and conditions. This ranges from the calculation of the primordial nucleosynthesis, where only a few nuclei are considered, to the ejecta of neutron star mergers with several thousands of involved nuclei. Here we describe the underlying physics and implementation details of the reaction network. We additionally present the numerical implementation of two different integration methods, the implicit Euler method and Gears method, along with their advantages and disadvantages. We furthermore describe basic example cases of thermodynamic conditions that we provide together with the network and demonstrate the reliability of the code by using simple test cases. With this publication, WinNet will be publicly available and open source at GitHub and Zenodo.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Radio Afterglow Detection and AI-driven Response (RADAR): A Federated Framework for Gravitational-wave Event Follow-up

The landmark detection of both gravitational waves (GWs) and electromagnetic (EM) radiation from the binary neutron star merger GW170817 has spurred efforts to streamline the follow-up of GW alerts in current and future observing runs of ground-based GW detectors. Within this context, the radio band of the EM spectrum presents unique challenges. Sensitive radio facilities capable of detecting the faint radio afterglow seen in GW170817, and with sufficient angular resolution, have small fields of view compared to typical GW localization areas. Additionally, theoretical models predict that the radio emission from binary neutron star mergers can evolve over weeks to years, necessitating long-term monitoring to probe the physics of the various postmerger ejecta components. These constraints, combined with limited radio observing resources, make the development of more coordinated follow-up strategies essential—especially as the next generation of GW detectors promises a dramatic increase in detection rates. Here, we present RADAR, a framework designed to address these challenges by promoting community-driven information sharing, federated data analysis, and system resilience, while integrating AI methods for both GW signal identification and radio data aggregation. We show that it is possible to preserve data rights while sharing models that can help design and/or update follow-up strategies. We demonstrate our approach through a case study of GW170817, and discuss future directions for refinement and broader application.

Gravitational waves↗
DOE OSTI2020

Monte-Carlo Neutrino Transport in Neutron Star Merger Simulations

Gravitational waves and electromagnetic signals from merging neutron star binaries provide valuable information about the the properties of dense matter, the formation of heavy elements, and high-energy astrophysics. To fully leverage observations of these systems, we need numerical simulations that provide reliable predictions for the properties of the matter unbound in these mergers. An important limitation of current simulations is the use of approximate methods for neutrino transport that do not converge to a solution of the transport equations as numerical resolution increases, and thus have errors that are impossible to quantify. Here, we report on a first simulation of a binary neutron star merger that directly solves the transport equations using Monte-Carlo techniques. We demonstrate that a computationally cheap Monte-Carlo transport can be implemented without ill effects from sampling noise. We apply Monte-Carlo transport to an unequal mass neutron star binary merger up to 5ms past merger, and report on the properties of the matter and neutrino outflows. Finally, we compare our results to the output of our best approximate `M1' transport scheme, demonstrating that an M1 scheme that carefully approximates the neutrino energy spectrum only leads to ~10% uncertainty in the composition and velocity of the ejecta, and ~20% uncertainty in the electron (anti)neutrino luminosities and energies. Here, the most significant disagreement found between M1 and Monte-Carlo results is a factor of ~2 difference in the luminosity of heavy-lepton neutrinos.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

The Relative Contribution to Heavy Metals Production from Binary Neutron Star Mergers and Neutron Star–Black Hole Mergers

The origin of the heavy elements in the universe is not fully determined. Neutron star–black hole (NSBH) and binary neutron star (BNS) mergers may both produce heavy elements via rapid neutron-capture (r-process). We use the recent detection of gravitational waves from NSBHs, improved measurements of the neutron star equation of state (EoS), and the most modern numerical simulations of ejected material from binary collisions to measure the relative contribution of NSBHs and BNSs to the production of heavy elements. As the amount of r-process ejecta depends on the mass and spin distribution of the compact objects, as well as on the EoS of the neutron stars, we consider various models for these quantities, informed by gravitational-wave and pulsar data. We find that in most scenarios, BNSs have produced more r-process elements than NSBHs over the past 2.5 billion years. If black holes have preferentially small spins, BNSs can produce at least twice the amount of r-process elements than NSBHs. If black hole spins are small and there is a dearth of low-mass (<5M ⊙ ) black holes within NSBH binaries, BNSs can account for the near totality of the r-process elements from binaries. For NSBH to produce a large fraction of r-process elements, black holes in NSBHs must have small masses and large aligned spins, which is disfavored by current data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Constraints on the Electromagnetic Counterpart of the Neutron-star–Black-hole Merger GW200115

We report the results of our follow-up campaign for the neutron-star—black-hole (NSBH) merger GW200115 detected during the O3 run of the Advanced LIGO and Advanced Virgo detectors. We obtained wide-field observations with the Deca-Degree Optical Transient Imager covering ~20% of the total probability area down to a limiting magnitude of w = 20.5 AB at ~23 hr after the merger. Our search for counterparts returns a single candidate (AT2020aeo), likely not associated with the merger. In total, only 25 sources of interest were identified by the community and later discarded as unrelated to the GW event. We compare our upper limits with the emission predicted by state-of-the-art kilonova simulations and disfavor high-mass ejecta (>0.1 M⊙), indicating that the spin of the system is not particularly high. By combining our optical limits with gamma-ray constraints from Swift and Fermi, we disfavor the presence of a standard short-duration burst for viewing angles ≲15° from the jet axis. Our conclusions are, however, limited by the large localization region of this GW event, and accurate prompt positions remain crucial to improving the efficiency of follow-up efforts.

79 ASTRONOMY AND ASTROPHYSICS↗