Search NASA⌕ Search

SEARCH · Search NASA

Results for “kilonova”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

A Unified Model of Kilonovae and Gamma-Ray Bursts in Binary Mergers Establishes Neutron Stars as the Central Engines of Short GRBs

We expand the theoretical framework by O. Gottlieb et al., which connects binary merger populations with long and short binary gamma-ray bursts (lbGRBs and sbGRBs, respectively), incorporating kilonovae (KNe) as a key diagnostic tool. We show that lbGRBs, powered by massive accretion disks around black holes (BHs), should be accompanied by bright, red KNe. In contrast, sbGRBs—if also powered by BHs—would produce fainter, red KNe, potentially biasing against their detection. However, magnetized hypermassive neutron star (HMNS) remnants that precede BH formation can produce jets with power (P NS ≈ 10 51 erg s –1 ) and Lorentz factor (Γ > 10) likely compatible with sbGRB observations, and would result in distinctly bluer KNe, offering a pathway to identifying the sbGRB central engine. Recent modeling by J. C. Rastinejad et al. found luminous red KNe consistently accompany lbGRBs, supporting their origin in BH-massive disk systems, likely following a short-lived HMNS phase. The preferential association of sbGRBs with comparably luminous KNe argues against the BH engine hypothesis for sbGRBs, while the bluer hue of these KNe provides additional support for an HMNS-driven mechanism. Within this framework, BH–NS mergers likely contribute exclusively to the lbGRB population with red KNe. Our findings suggest that GW170817 may, in fact, have been an lbGRB to on-axis observers. Finally, we discuss major challenges faced by alternative lbGRB progenitor models, such as white dwarf–NS or white dwarf–BH mergers and accretion-induced collapse forming magnetars, which fail to align with observed GRB timescales, energies, and KN properties.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Kilonova Emissions from Neutron Star Merger Remnants: Implications for the Nuclear Equation of State

Multimessenger observations of binary neutron star mergers can provide valuable information on the nuclear equation of state (EOS). Here, we investigate the extent to which electromagnetic observations of the associated kilonovae allow us to place constraints on the EOS. For this, we use state-of-the-art three-dimensional general-relativistic magnetohydrodynamics simulations and detailed nucleosynthesis modeling to connect properties of observed light curves to properties of the accretion disk, and hence, the EOS. Using our general approach, we use multimessenger observations of GW170817/AT2017gfo to study the impact of various sources of uncertainty on inferences of the EOS. We constrain the radius of a 1.4M ⊙ neutron star to lie within 10.30 ≤ R 1.4 ≤ 13.0 km and the maximum mass to be M TOV ≤ 3.06M ⊙ .

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

AT2025ulz and S250818k: Zooming in with the Hubble Space Telescope

AT2025ulz is an optical/near-infrared transient discovered during follow-up of the candidate gravitational wave (GW) event S250818k. Its young age (≲1 day), rapid decline, and strong color evolution over the first 48 hr classify it as a potential kilonova candidate. In this work, we present the results of our observing campaign, carried out with the Gran Telescopio Canarias and the Hubble Space Telescope (HST). Although the early time evolution of AT2025ulz resembles some aspects of a kilonova, its rapid onset (∼3 hr after the GW trigger) and luminosity (a factor of ∼5 brighter than AT2017gfo in g band) are difficult to reproduce. Only a small subset of our kilonova models matches its multicolor light curve, and the inferred ejecta mass is uncomfortably large given the low chirp mass (≲0.87 M ⊙ ) of the GW candidate. HST observations place the transient within a nearby (z = 0.08489) spiral galaxy with ongoing star formation and measure a color (F336W − F160W ≈ 1.4 mag) that is too blue to match with a kilonova. Our data support the classification of AT2025ulz as a supernova, initially undergoing a shock-cooling phase and later entering its photospheric phase and spectroscopically identified via its broad absorption features.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Do we owe our existence to gravitational waves?

Two heavy elements essential to human biology are thought to have been produced by the astrophysical r-process, which occurs in neutron-rich environments: iodine is a constituent of thyroid hormones that affect many physiological processes including growth and development, body temperature and heart rate, and bromine is essential for tissue development and architecture. Collisions of neutron stars (kilonovae) have been identified as sources of r-process elements including tellurium, which is adjacent to iodine in the periodic table, and lanthanides. Neutron-star collisions arise from energy loss due to gravitational-wave emission from binary systems, leading us to suggest that gravitational waves have played a key role in enabling human life by producing iodine and bromine. We propose probing this proposal by searching in lunar material for live 129 I deposited by a recent nearby kilonova explosion.

Human biochemistry↗
DOE OSTI2025

Black hole-neutron star binaries near neutron star disruption limit in the mass regime of event GW230529

In May 2023, the LIGO-Virgo-KAGRA (LVK) Collaboration reported the likely black hole-neutron star (BHNS) merger GW230529_181500. The signal was observed with high significance in only one detector, limiting constraints on the black hole spin and motivating our study of disruption in this mass regime. That event is expected to be the merger of a 2.5–4.5 M ⊙ primary with a secondary compact object of mass between 1.2–2.0 M ⊙ . This makes it the first BHNS merger with a significant potential for the production of electromagnetic (EM) counterparts, and provides further evidence for compact objects existing within the suspected lower mass gap. To produce post-merger EM transients, the component of the black hole spin aligned with the orbital angular momentum must be sufficiently high, allowing the neutron star to be tidally disrupted. The disrupting BHNS binary may then eject a few percent of a solar mass of matter, leading to an observable kilonova driven by radioactive decays in ejecta, and/or a compact-binary gamma-ray burst (cbGRB) resulting from the formation of an accretion disk and relativistic jet. Determining which mergers lead to disruption of the neutron star is necessary to predict the prevalence of EM signals from BHNS mergers, yet most BHNS simulations so far have been performed far from the minimum spin required for tidal disruption. Here, we use the Spectral Einstein Code to explore the behavior of BHNS mergers in a mass range consistent with GW230529_181500 close to that critical spin, and compare our results against the mass remnant model currently used by the LVK Collaboration to predict the probability of tidal disruption. Furthermore, our numerical results reveal the emergence of non-zero accretion disks even below the predicted NS disruption limit, of low mass but capable of powering cbGRBs. Our results also demonstrate that the remnant mass model underpredicts the disk mass for the DD2 equation of state, while they are within expected modeling errors for SFHo. The disruption limit itself, however, is not found to significantly differ from the predictions of the analytical model, unless remnant masses M rem ≲ 0.001 M ⊙ prove interesting observationally. In all of our simulations, any kilonova signal would be dim and most likely dominated by post-merger disk outflows.

GW230529_181500↗
DOE OSTI2026

Exploring the connection between compact object mergers and fast X-ray transients

Context. The connection between compact object mergers and some extragalactic fast X-ray transients (FXRTs) has long been hypothesised but never ultimately established. Aims. In this work, we investigate two FXRTs, the LEIA X-ray Transient LXT 240402A and the Einstein Probe EP 250207b, whose precise positions lie close to nearby (z ≲ 0.1) quiescent galaxies with a negligible probability of chance coincidence, identifying them as particularly promising cases of merger-driven explosions in the local Universe. Methods. We used Chandra to derive accurate localisations for both events and secure otherwise ambiguous associations with their optical counterparts. Deep optical and near-infrared observations with VLT, GTC, and LBT were performed to characterise the surrounding environment and search for kilonova emission, the hallmark of neutron star mergers. Complementary early-time X-ray monitoring with Swift and Einstein Probe was used to constrain the non-thermal afterglow. Results. We find that both FXRTs remain compatible with a compact binary merger progenitor, which produced low-mass ejecta and kilonova emission subdominant to the afterglow. However, alternative explanations such as a distant (z ≳ 1) core-collapse supernova cannot be conclusively ruled out.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

New Ideas on the Formation and Astrophysical Detection of Primordial Black Holes

Recently, a number of novel scenarios for primordial black hole (PBH) formation have been discovered. Some of them require very minimal new physics, some others require no new ingredients besides those already present in commonly considered models, such as supersymmetry. At the same time, new strategies have emerged for detection of PBHs. For example, an observation of an orphan kilonova unaccompanied by the gravitational waves signal of merging neutron stars, but associated with a fast radio burst, could be a smoking gun of PBH dark matter. We review some new ideas for PBH formation and detection.

Flores, Marcos M.↗
DOE OSTI2025

Assessing the difficulty of capturing the distribution function of neutrinos in neutron star merger simulations

The collision of two neutron stars is a rich source of information about nuclear physics. In particular, the kilonova signal following a merger can help us elucidate the role of neutron stars in nucleosynthesis, and informs us about the properties of matter above nuclear saturation. Approximate modeling of neutrinos remains an important limitation to our ability to make predictions for these observables. Part of the problem is the fermionic nature of neutrinos. By the exclusion principle, the expected value 𝑓 𝜈 for the number of neutrinos in a quantum state is at most 1. Any process producing neutrinos is suppressed by a blocking factor (1 −𝑓 𝜈 ). Recent simulations focused on neutrino physics mostly use a gray two-moment scheme to evolve neutrinos. This evolves integrals of 𝑓 𝜈 over momentum space, preventing direct calculations of blocking factors. Monte Carlo methods may be an attractive alternative, providing access to the full distribution of neutrinos. Their current implementation is, however, inadequate to estimate 𝑓 𝜈 : in our most recent simulations, a single Monte Carlo packet causes, in the worst cases, estimates of 𝑓 𝜈 to jump from 𝑓 𝜈 =0 to 𝑓 𝜈 ∼10 5 . While this is concerning, this brazen violation of the fermionic nature of neutrinos has been largely inconsequential, as the interactions used in simulations avoid direct calculations of 𝑓 𝜈 . We are, however, reaching a level of modeling at which this problem can no longer be ignored. Here, we discuss the relatively simple origin of this issue. We then show that very rough estimates of 𝑓 𝜈 can in theory be obtained in merger simulations, but that they will require a combination of unintuitive weighting schemes for Monte Carlo packets and smoothing of the neutrino distribution at coarser resolution than what the merger simulation uses.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Impact of Newly Measured 𝛽-Delayed Neutron Emitters around 78 Ni on Light Element Nucleosynthesis in the Neutrino Wind Following a Neutron Star Merger

Neutron emission probabilities and half-lives of 37 𝛽-delayed neutron emitters from 75 Ni to 92 Br were measured at the RIKEN Nishina Center in Japan, including 11 one-neutron and 13 two-neutron emission probabilities and six half-lives for the first time that supersede theoretical estimates. These nuclei lie in the path of the weak 𝑟 process occurring in neutrino-driven winds from the accretion disk formed after the merger of two neutron stars synthesizing elements in the 𝐴∼80 abundance peak. The presence of such elements dominates the accompanying kilonova emission over the first few days and have been identified in the AT2017gfo event, associated to the gravitational wave detection GW170817. Abundance calculations based on over 17,000 simulated trajectories describing the evolution of matter properties in the merger outflows show that the new data lead to an increase of 50%–70% in the abundance of Y, Zr, Nb, and Mo. This enhancement is large compared to the scatter of relative abundances observed in old very metal poor stars and thus is significant in the comparison with other possible astrophysical processes contributing to the light-element production. These results underline the importance of including experimental decay data for very neutron-rich 𝛽 -delayed neutron emitters into 𝑟 -process models.

59 ≤ A ≤ 89↗
DOE OSTI2025

From Existing and New Nuclear and Astrophysical Constraints to Stringent Limits on the Equation of State of Neutron-Rich Dense Matter

Through continuous progress in nuclear theory and experiment and an increasing number of neutron-star (NS) observations, a multitude of information about the equation of state (EOS) for matter at extreme densities is available. To constrain the EOS across its entire density range, this information needs to be combined consistently. However, the impact and model dependency of individual observations vary. Given their growing number, assessing the various methods is crucial to compare the respective effects on the EOS and discover potential biases. For this purpose, we present a broad compendium of different constraints and apply them individually to a large set of EOS candidates within a Bayesian framework. Specifically, we explore different ways of how chiral effective field theory and perturbative quantum chromodynamics can be used to place a likelihood on EOS candidates. We also investigate the impact of nuclear experimental constraints, as well as different radio and x-ray observations of NS masses and radii. This is augmented by reanalyses of the existing data from binary neutron star coalescences, in particular of GW170817, with improved models for the tidal waveform and kilonova light curves, which we also utilize to construct a tight upper limit of 2.39 M ⊙ on the TOV mass based on GW170817’s remnant. Our diverse set of constraints is eventually combined to obtain stringent limits on NS properties. We organize the combination in a way to distinguish between constraints where the systematic uncertainties are deemed small and those that rely on less conservative assumptions. For the former, we find the radius of the canonical 1.4 M ⊙ neutron star to be R 1.4 = 12.2 6 − 0.91 + 0.80 km and the TOV mass at M TOV = 2.2 5 − 0.22 + 0.42 M ⊙ (95% credibility). Including all the presented constraints yields R 1.4 = 12.2 0 − 0.48 + 0.50 km and M TOV = 2.3 0 − 0.20 + 0.07 M ⊙ . When comparing these limits to individual data points, we find that the quoted radius of HESS J1731-347 displays noticeable tension with other constraints. Constraining microphysical properties of the EOS proves more challenging. For instance, the symmetry energy slope is restricted to L sym = 48 − 25 + 21 MeV , where this constraint is mainly dominated by our reanalysis of the PREX-II and CREX experiment. Published by the American Physical Society 2025

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Towards Exascale Astrophysics of Mergers and Supernovae (TEAMS)

The TEAMS project brought together cutting-edge simulations, theoretical insights, and collaborative efforts to deepen our understanding of some of the universe’s most extreme phenomena—supernovae, neutron star mergers, and the powerful signals they emit. Using one of the largest suites of 3D supernova simulations ever conducted, researchers uncovered new insights into how massive stars explode, how those explosions vary by stellar mass, and what conditions lead to the birth of neutron stars or black holes. They also studied the radiation and gravitational wave signals emitted during these events, revealing how future observations can be used to uncover what happens deep inside collapsing stars. The team developed improved tools for modeling how light and neutrinos behave in such explosive environments, enabling more accurate predictions of what astronomers might observe. Work also explored how the chemical composition and geometry of kilonovae—the visible explosions that follow neutron star mergers—influence their signals and can reveal the origins of heavy elements like gold. These efforts not only advanced scientific knowledge, but also trained a new generation of researchers at the intersection of astrophysics, computational science, and nuclear theory.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

The Gravity Collective: A Comprehensive Analysis of the Electromagnetic Search for the Binary Neutron Star Merger GW190425

We present an ultraviolet to infrared search for the electromagnetic (EM) counterpart to GW190425, the second ever binary neutron star merger discovered by the LIGO-Virgo-KAGRA Collaboration. GW190425 was more distant and had a larger localization area than GW170817, so we use a new tool, Teglon, to redistribute the GW190425 localization probability in the context of galaxy catalogs within the final localization volume. We derive a 90th percentile area of 6688 deg 2 , a ∼1.5× improvement relative to the LIGO/Virgo map, and show how Teglon provides an order-of-magnitude boost to the search efficiency of small (≤1 deg 2 ) field-of-view instruments. We combine our data with a large, publicly reported imaging data set, covering 9078.59 deg 2 of unique area and 48.13% of the LIGO/Virgo-assigned localization probability, to calculate the most comprehensive kilonova (KN), short gamma-ray burst (sGRB) afterglow, and model-independent constraints on the EM emission from a hypothetical counterpart to GW190425 to date under the assumption that no counterpart was found in these data. If the counterpart were similar to AT 2017gfo, there would be a 28.4% chance of it being detected in the combined data set. We are relatively insensitive to an on-axis sGRB, and rule out a generic transient with a similar peak luminosity and decline rate as AT 2017gfo to 30% confidence. Finally, across our new imaging and publicly reported data, we find 28 candidate optical counterparts that we cannot rule out as being associated with GW190425, finding that four such counterparts discovered within the localization volume and within 5 days of merger exhibit luminosities consistent with a KN.

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 OSTI2024

Fragmentation in Gravitationally Unstable Collapsar Disks and Subsolar Neutron Star Mergers

Abstract Although stable neutron stars (NSs) can in principle exist down to massesM ns ≈ 0.1M ⊙ , standard models of stellar core-collapse predict a robust lower limitM ns ≳ 1.2M ⊙ , roughly commensurate with the Chandrasekhar massM Ch of the progenitor’s iron core (electron fractionY e ≈ 0.5). However, this limit may be circumvented in sufficiently dense neutron-rich environments (Y e < 0.5) for which M Ch ∝ Y e 2 is reduced to ≲1M ⊙ . Such physical conditions could arise in the black hole accretion disks formed from the collapse of rapidly rotating stars (“collapsars”), as a result of gravitational instabilities and cooling-induced fragmentation, similar to models for planet formation in protostellar disks. We confirm that the conditions to form subsolar-mass NS (ssNS) may be marginally satisfied in the outer regions of massive neutrino-cooled collapsar disks. If the disk fragments into multiple ssNSs, their subsequent coalescence offers a channel for precipitating subsolar mass LIGO/Virgo gravitational-wave mergers that does not implicate primordial black holes. The model makes several additional predictions: (1) ∼Hz frequency Doppler modulation of the ssNS-merger gravitational-wave signals due to the binary’s orbital motion in the disk; (2) at least one additional gravitational-wave event (coincident within ≲hours), from the coalescence of the ssNS-merger remnant(s) with the central black hole; (3) an associated gamma-ray burst and supernova counterpart, the latter boosted in energy and enriched withr-process elements from the NS merger(s) embedded within the exploding stellar envelope (“kilonovae inside a supernova”).

Astronomy & Astrophysics↗
DOE OSTI2025

Multimessengers from the Radioactive Decay of r -process Nuclei

The radioactive β-decay of nuclei synthesized in the rapid neutron capture process (r-process) releases a variety of particles, including electrons, γ-rays, neutrinos, and neutrons. These particles provide a rich set of multimessenger signals that carry information about the astrophysical environments where neutron-rich nucleosynthesis occurs. In this work, we calculate from first principles the emission spectra resulting from the β-decay of r-process nuclei. Our approach incorporates detailed nuclear structure and decay data to model the energy distributions of each particle species. We couple the spectra with a nuclear reaction network simulation to obtain the temporal evolution of these distributions. We find that the emission distributions vary significantly in time and are nonthermal, with substantial average energies. We investigate these nuclear signals as a direct probe of heavy element formation and show that they are complementary observables to kilonova.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗