Search NASA⌕ Search

SEARCH · Search NASA

Results for “kilonovae”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Kilonovae and Long-duration Gamma-Ray Bursts

Recent detections of kilonova-like emission following long-duration gamma-ray bursts GRB 211211A and GRB 230307A have been interpreted as originating from the merger of two neutron stars. In this work, we demonstrate that these observations are also consistent with nucleosynthesis originating from a collapsar scenario. Our model is capable of reproducing the observed optical and infrared light curves using a single, weak r-process component. The absence of lanthanide-rich material in our model, consistent with the data, challenges the prevailing interpretation that a red evolution in such transients necessarily indicates the presence of heavy r-process elements.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Black hole–neutron star mergers: The first mass gap and kilonovae

Observations of X-ray binaries indicate a dearth of compact objects in the mass range from ~2 –5 M ⊙ . The existence of this (first mass) gap has been used to discriminate between proposed engines behind core-collapse supernovae. From LIGO/Virgo observations of binary compact remnant masses, several candidate first mass gap objects, either neutron stars (NSs) or black holes (BHs), were identified during the O3 science run. Motivated by these new observations, we study the formation of BH-NS mergers in the framework of isolated classical binary evolution, using population synthesis methods to evolve large populations of binary stars (Population I and II) across cosmic time. We present results on the NS to BH mass ratios (q = M NS /M BH ) in merging systems, showing that although systems with a mass ratio as low as q = 0.02 can exist, typically BH-NS systems form with moderate mass ratios q = 0.1 –0.2. If we adopt a delayed supernova engine, we conclude that ~30% of BH-NS mergers may host at least one compact object in the first mass gap (FMG). Even allowing for uncertainties in the processes behind compact object formation, we expect the fraction of BH-NS systems ejecting mass during the merger to be small (from ~0.6 –9%). In our reference model, we assume: (i) the formation of compact objects within the FMG, (ii) natal NS/BH kicks decreased by fallback, (iii) low BH spins due to Tayler-Spruit angular momentum transport in massive stars. We find that ≲1% of BH-NS mergers will have any mass ejection and about the same percentage will produce kilonova bright enough to have a chance of being detected with a large (Subaru-class) 8 m telescope. Interestingly, all these mergers will have both a BH and an NS in the FMG.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Bayesian inference of multi-messenger astrophysical data: Joint and coherent inference of gravitational waves and kilonovae

Multi-messenger observations of binary neutron star mergers can provide information on the neutron star’s equation of state (EOS) above the nuclear saturation density by directly constraining the mass-radius diagram. We present a Bayesian framework for joint and coherent analyses of multi-messenger binary neutron star signals. As a first application, we analyze the gravitational-wave GW170817 and the kilonova (kN) AT2017gfo data. These results are then combined with the most recent X-ray pulsar analyses of PSR J0030+0451 and PSR J0740+6620 to obtain new EOS constraints.We extend the bajes infrastructure with a joint likelihood for multiple datasets, support for various semi-analytical kN models, and numerical-relativity (NR)-informed relations for the mass ejecta, as well as a technique to include and marginalize over modeling uncertainties. The analysis of GW170817 used the TEOBResumS effective-one-body waveform template to model the gravitational-wave signal. The analysis of AT2017gfo used a baseline multicomponent spherically symmetric model for the kN light curves. Various constraints on the mass-radius diagram and neutron star properties were then obtained by resampling over a set of ten million parameterized EOSs, which was built under minimal assumptions (general relativity and causality).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Actinide opacities for modelling the spectra and light curves of kilonovae

ABSTRACT We extend previous ab initio calculations of lanthanide opacities to include a complete set of actinide opacities for use in the modelling of kilonova (KN) light curves and spectra. Detailed, fine-structure line features are generated using the configuration-interaction approach. These actinide opacities display similar trends to those observed for lanthanide opacities, such as the lighter actinides producing higher opacity than the heavier ones for relevant conditions in the dynamical ejecta. A line-binned treatment is employed to pre-compute opacity tables for 14 actinide elements (89 ≤ Z ≤ 102) over a grid of relevant temperatures and densities. These tabular opacities will be made publicly available for general usage in KN modelling. We demonstrate the usefulness of these opacities in KN simulations by exploring the sensitivity of light curves and spectra to different actinide abundance distributions that are predicted by different nuclear theories, as well as to different choices of ejecta mass and velocity. We find very little sensitivity to the two considered distributions, indicating that opacities for actinides with Z ≥ 99 do not contribute strongly. On the other hand, a single actinide element, protactinium, is found to produce faint spectral features in the far-infrared at late times (5–7 d post merger). More generally, we find that the choice of ejecta mass and velocity have the most significant effect on KN emission for this study.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Opacities of singly and doubly ionized neodymium and uranium for kilonova emission modeling

ABSTRACT Even though the electromagnetic counterpart AT2017gfo to the binary neutron star merger GW170817 is powered by the radioactive decay of r-process nuclei, only few tentative identifications of light r-process elements have been made so far. One of the major limitations for the identification of heavy nuclei is incomplete or missing atomic data. While substantial progress has been made on lanthanide atomic data over the last few years, for actinides there has been less emphasis, with the first complete set of opacity data only recently published. We perform atomic structure calculations of neodymium (Z = 60) as well as the corresponding actinide uranium (Z = 92). Using two different codes [flexible atomic code (fac) and hartree–fock-relativistic (hfr)] for the calculation of the atomic data, we investigate the accuracy of the calculated data (energy levels and electric dipole transitions) and their effect on kilonova opacities. For the fac calculations, we optimize the local central potential and the number of included configurations and use a dedicated calibration technique to improve the agreement between theoretical and available experimental atomic energy levels (AELs). For ions with vast amounts of experimental data available, the presented opacities agree quite well with previous estimations. On the other hand, the optimization and calibration method cannot be used for ions with only few available AELs. For these cases, where no experimental nor benchmarked calculations are available, a large spread in the opacities estimated from the atomic data obtained with the various atomic structure codes is observed. We find that the opacity of uranium is almost double the neodymium opacity.

Flörs, A. (ORCID:0000000320242819)↗
DOE OSTI2023

Nucleosynthesis in outflows of compact objects and detection prospects of associated kilonovae

ABSTRACT We perform a comparative analysis of nucleosynthesis yields from binary neutron star (BNS) mergers, black hole-neutron star (BHNS) mergers, and core-collapse supernovae (CCSNe) with the goal of determining which are the most dominant sources of r-process enrichment observed in stars. We find that BNS and BHNS binaries may eject similar mass distributions of robust r-process nuclei post-merger (up to third peak and actinides, A ∼ 200−240), after accounting for the volumetric event rates. Magnetorotational (MR) CCSNe likely undergo a weak r-process (up to A ∼ 140) and contribute to the production of light element primary process (LEPP) nuclei, whereas typical thermal, neutrino-driven CCSNe only synthesize up to first r-process peak nuclei (A ∼ 80−90). We also find that the upper limit to the rate of MR CCSNe is $\lesssim 1~{{\ \rm per\ cent}}$ the rate of typical thermal CCSNe; if the rate was higher, then weak r-process nuclei would be overproduced. Although the largest uncertainty is from the volumetric event rate, the prospects are encouraging for confirming these rates in the next few years with upcoming surveys. Using a simple model to estimate the resulting kilonova light curve from mergers and our set of fiducial merger parameters, we predict that ∼7 BNS and ∼2 BHNS events will be detectable per year by the Vera C. Rubin Observatory (LSST), with prior gravitational wave (GW) triggers.

Astronomy & Astrophysics↗
DOE OSTI2024

Actinide signatures in low electron fraction kilonova ejecta

ABSTRACT Neutron star (NS) mergers are known to produce heavy elements through rapid neutron capture (r-process) nucleosynthesis. Actinides are expected to be created solely by the r-process in the most neutron-rich environments. Confirming if NS mergers provide the requisite conditions for actinide creation is therefore central to determining their origin in the Universe. Actinide signatures in kilonova (KN) spectra may yield an answer, provided adequate models are available in order to interpret observational data. In this study, we investigate actinide signatures in neutron-rich merger ejecta. We use three ejecta models with different compositions and radioactive power, generated by nucleosynthesis calculations using the same initial electron fraction ($Y_e = 0.15$) but with different nuclear physics inputs and thermodynamic expansion history. These are evolved from 10 to 100 d after merger using the sumo non-local thermodynamic equilibrium (NLTE) radiative transfer code. We highlight how uncertainties in nuclear properties, as well as choices in thermodynamic trajectory, may yield entirely different outputs for equal values of $Y_e$. We consider an actinide-free model and two actinide-rich models, and find that the emergent spectra and light-curve evolution are significantly different depending on the amount of actinides present, and the overall decay properties of the models. We also present potential key actinide spectral signatures, of which doubly ionized $_{89}$Ac and $_{90}$Th may be particularly interesting as spectral indicators of actinide presence in KN ejecta.

Pognan, Quentin (ORCID:0000000152550782)↗
DOE OSTI2026

Exploring composition mixing in kilonova ejecta with ray-by-ray simulations

Binary neutron star merger (BNSM) ejecta are considered a primary repository of r-process nucleosynthesis and a source of the observed heavy-element abundances. We implement composition mixing into ray-by-ray radiation-hydrodynamic simulations of BNSM ejecta, coupled with an online nuclear network (NN). We model mixing via a gradient-based mixing approximation that evolves simultaneously with the hydrodynamics. Here, we find that mixing occurs in regions where the electron fraction changes rapidly. While mixing smooths composition gradients in transition regions, it has a negligible impact on the heavy-element yields. This is because the primary r-process site (the equatorial ejecta) is initially homogeneous in free neutrons, leaving no strong gradients for mixing to act upon. In each angular ray, the abundances of the most produced elements are robust under mixing, while the less abundant ones are more affected. The total global abundances change only slightly from mixing, since each angular ray contributes its most abundant elements. Furthermore, the predicted kilonova light curves show only minor reddening, with differences below the detectability of state-of-the-art telescopes. In general, we do not observe significant effects from mixing in the time span of the r-process. Consequently, mixing only leads to minor variations in abundances and light curves in ray-by-ray simulations.

Explosive burning↗
DOE OSTI2022

Kilonova Detectability with Wide-Field Instruments

A diverse set of instruments increase the chances of detecting and identifying a kilonova. Early observations increase the probability of detection. More sensitive wide-field ultraviolet instruments are needed as GW detectors reach design sensitivity. There is a dearth of wide-field infrared instruments.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Light Curves and Spectra from Kilonova Models [Poster]

White dwarf-neutron star merger simulations show a Ca IR line feature in spectra, pointing to Ca-rich transients as a possible observation (Kaltenborn et al 2024, Gillanders et al 2024). Supernova and new/old kilonova (KN) models have been compared to recent peculiar long gamma-ray burst (GRB) with KN-like signal (Troja et al 2022).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

The X-Ray Fundamental Plane of the Platinum Sample, the Kilonovae, and the SNe Ib/c Associated with GRBs

A large fraction of gamma-ray burst (GRB) lightcurves (LCs) show X-ray plateaus. We analyze all GRBs with known redshifts presenting plateaus observed by The Neil Gehrels Swift Observatory from its launch until 2019 August. The fundamental plane relation between the rest-frame time and X-ray luminosity at the end of the plateau emission and the peak prompt luminosity holds for all the GRB classes when selection biases and cosmological evolutions are applied. We have discovered two important findings: (1) a new class of long GRBs with good data coverage: the platinum sample; and (2) the platinum, the SNe-LGRB, and the KN-SGRB samples, yield the smallest intrinsic scatter with σ platinum , GRB-SNe = 0.22 ± 0.10 and σ KN-SGRB = 0.24 ± 0.12. The SNe-LGRBs are composed of GRBs associated spectroscopically with the SNe Ib,c, the KN-SGRBs are composed by eight GRBs associated with kilonovae or where there could have been such an association. The highest correlation coefficients are yielded for the SN-LGRB-ABC sample, which includes GRBs spectroscopically associated with SNe Ib/c or with a clear optical bump in the LC resembling the SNe Ib/c, ( ), the SN-LGRBs ( ), and the KN-SGRBs ( ) when the redshift evolution is considered. We find that these category planes are reliable candidates to use as cosmological tools. Furthermore, the distance from the gold fundamental plane is a crucial discriminant among classes. In fact, we find that the distributions of the distances of the SNe-LGRB, SNe-LGRB-ABC, KN-SGRB, and SGRB samples from the gold fundamental plane are statistically different from the distribution of the gold GRBs' distances from the gold fundamental plane with and without considering evolution cases. <

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

The Challenges Ahead for Multimessenger Analyses of Gravitational Waves and Kilonova: A Case Study on GW190425

In recent years, there have been significant advances in multimessenger astronomy due to the discovery of the first, and so far only confirmed, gravitational wave event with a simultaneous electromagnetic (EM) counterpart, as well as improvements in numerical simulations, gravitational wave (GW) detectors, and transient astronomy. This has led to the exciting possibility of performing joint analyses of the GW and EM data, providing additional constraints on fundamental properties of the binary progenitor and merger remnant. Here, in this work, we present a new Bayesian framework that allows inference of these properties, while taking into account the systematic modeling uncertainties that arise when mapping from GW binary progenitor properties to photometric light curves. We extend the relative binning method presented in Zackay et al. to include extrinsic GW parameters for fast analysis of the GW signal. The focus of our EM framework is on light curves arising from r-process nucleosynthesis in the ejected material during and after merger, the so-called kilonova, and particularly on black hole–neutron star systems. As a case study, we examine the recent detection of GW190425, where the primary object is consistent with being either a black hole or a neutron star. We show quantitatively how improved mapping between binary progenitor and outflow properties, and/or an increase in EM data quantity and quality are required in order to break degeneracies in the fundamental source parameters.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

r-Process Radioisotopes from Near-Earth Supernovae and Kilonovae

The astrophysical sites where r-process elements are synthesized remain mysterious: it is clear that neutron star mergers (kilonovae (KNe)) contribute, and some classes of core-collapse supernovae (SNe) are also possible sources of at least the lighter r-process species. The discovery of 60Fe on the Earth and Moon implies that one or more astrophysical explosions have occurred near the Earth within the last few million years, probably SNe. Intriguingly, 244 Pu has now been detected, mostly overlapping with 60 Fe pulses. However, the 244 Pu flux may extend to before 12 Myr ago, pointing to a different origin. Here, motivated by these observations and difficulties for r-process nucleosynthesis in SN models, we propose that ejecta from a KN enriched the giant molecular cloud that gave rise to the Local Bubble, where the Sun resides. Accelerator mass spectrometry (AMS) measurements of 244 Pu and searches for other live isotopes could probe the origins of the r-process and the history of the solar neighborhood, including triggers for mass extinctions, e.g., that at the end of the Devonian epoch, motivating the calculations of the abundances of live r-process radioisotopes produced in SNe and KNe that we present here. Given the presence of 244 Pu, other r-process species such as 93 Zr, 107 Pd, 129 I, 135 Cs, 182 Hf, 236 U, 237 Np, and 247 Cm should be present. Their abundances and well-resolved time histories could distinguish between the SN and KN scenarios, and we discuss prospects for their detection in deep-ocean deposits and the lunar regolith. We show that AMS 129 I measurements in Fe–Mn crusts already constrain a possible nearby KN scenario.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Could a Kilonova Kill: A Threat Assessment

Abstract Binary neutron star mergers produce high-energy emissions from several physically different sources, including a gamma-ray burst (GRB) and its afterglow, a kilonova (KN), and, at late times, a remnant many parsecs in size. Ionizing radiation from these sources can be dangerous for life on Earth-like planets when located too close. Work to date has explored the substantial danger posed by the GRB to on-axis observers; here we focus instead on the potential threats posed to nearby off-axis observers. Our analysis is based largely on observations of the GW170817/GRB 170817A multi-messenger event, as well as theoretical predictions. For baseline KN parameters, we find that the X-ray emission from the afterglow may be lethal out to ∼1 pc and the off-axis gamma-ray emission may threaten a range out to ∼4 pc, whereas the greatest threat comes years after the explosion, from the cosmic rays accelerated by the KN blast, which can be lethal out to distances up to ∼11 pc. The distances quoted here are typical, but the values have significant uncertainties and depend on the viewing angle, ejected mass, and explosion energy in ways we quantify. Assessing the overall threat to Earth-like planets, KNe have a similar kill distance to supernovae, but are far less common. However, our results rely on the scant available KN data, and multi-messenger observations will clarify the danger posed by such events.

Astronomy & Astrophysics↗
DOE OSTI2024

On a Spectral Method for β -particle Bound Excitation Collisions in Kilonovae

The interaction of β-particles with the weakly ionized plasma background is an important mechanism for powering the kilonova (KN) transient signal from neutron star mergers. For this purpose, we present an implementation of the approximate fast-particle collision kernel, described by Inokuti following the seminal formulation of Bethe, in a spectral solver of the Vlasov–Maxwell–Boltzmann equation. In particular, we expand the fast-particle plane-wave atomic excitation kernel into coefficients of the Hermite basis, and derive the relevant discrete spectral system. In this fast-particle limit, the approach permits the direct use of atomic data, including optical oscillator strengths, normally applied to photon–matter interaction. The resulting spectral matrix is implemented in the MASS-APP spectral solver framework, in a way that avoids full matrix storage per spatial zone. We numerically verify aspects of the matrix construction, and present a proof-of-principle 3D simulation of a 2D axisymmetric KN ejecta snapshot. Our preliminary numerical results indicate that a reasonable choice of Hermite basis parameters for β-particles in the KN is a bulk velocity parameter u = 0, a thermal velocity parameter α = 0.5c, and a 9 × 9 × 9 mode velocity basis set (Hermite orders of 0–8 in each dimension). For interior-ejecta sample zones, we estimate that the ratio of thermalization from large-angle (≳2fdg5) bound excitation scattering to total thermalization is ~0.002–0.003.

79 ASTRONOMY AND ASTROPHYSICS↗
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↗