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 127 records · Page 7

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↗

Neutrino emission from binary neutron star mergers: characterising light curves and mean energies

Abstract Neutrinos are copiously emitted by neutron star mergers, due to the high temperatures reached by dense matter during the merger and its aftermath. Neutrinos influence the merger dynamics and shape the properties of the ejecta, including the resulting r -process nucleosynthesis and kilonova emission. In this work, we analyse neutrino emission from a large sample of binary neutron star merger simulations in Numerical Relativity, covering a broad range of initial masses, nuclear equation of state and viscosity treatments. We extract neutrino luminosities and mean energies, and compute quantities of interest such as the peak values, peak broadnesses, time averages and decrease time scales. We provide a systematic description of such quantities, including their dependence on the initial parameters of the system. We find that for equal-mass systems the total neutrino luminosity (several $$10^{53}{\hbox {erg}~{\hbox {s}}^{-1}}$$ 10 53 erg s - 1 ) decreases as the reduced tidal deformability increases, as a consequence of the less violent merger dynamics. Similarly, tidal disruption in asymmetric mergers leads to systematically smaller luminosities. Peak luminosities can be twice as large as the average ones. Electron antineutrino luminosities dominate (initially by a factor of 2-3) over electron neutrino ones, while electron neutrinos and heavy flavour neutrinos have similar luminosities. Mean energies are nearly constant in time and independent on the binary parameters. Their values reflect the different decoupling temperature inside the merger remnant. Despite present uncertainties in neutrino modelling, our results provide a broad and physically grounded characterisation of neutrino emission, and they can serve as a reference point to develop more sophisticated neutrino transport schemes.

79 ASTRONOMY AND ASTROPHYSICS↗

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↗

Markov Chain Monte Carlo Predictions of Neutron-rich Lanthanide Properties as a Probe of r -process Dynamics

Lanthanide element signatures are key to understanding many astrophysical observables, from merger kilonova light curves to stellar and solar abundances. To learn about the lanthanide element synthesis that enriched our solar system, we apply the statistical method of Markov Chain Monte Carlo to examine the nuclear masses capable of forming the r-process rare-earth abundance peak. We describe the physical constraints we implement with this statistical approach and demonstrate the use of the parallel chains method to explore the multidimensional parameter space. We apply our procedure to three moderately neutron-rich astrophysical outflows with distinct types of r-process dynamics. We show that the mass solutions found are dependent on outflow conditions and are related to the r-process path. We describe in detail the mechanism behind peak formation in each case. We then compare our mass predictions for neutron-rich neodymium and samarium isotopes to the latest experimental data from the CPT at CARIBU. Here, we find our mass predictions given outflows that undergo an extended (n,γ)⇌(γ,n) equilibrium to be those most compatible with both observational solar abundances and neutron-rich mass measurements.

79 ASTRONOMY AND ASTROPHYSICS↗

Igniting Weak Interactions in Neutron Star Postmerger Accretion Disks

The merger of two neutron stars or a neutron star and a black hole typically results in the formation of a postmerger accretion disk. Outflows from disks may dominate the overall ejecta from mergers and be a major source of r-process nuclei in our universe. We explore the parameter space of such disks and their outflows and r-process yields by performing 3D general-relativistic magnetohydrodynamic simulations with weak interactions and approximate neutrino transport. We discuss the mapping between the initial binary parameters and the parameter space of the resulting disks, chiefly characterized by their initial accretion rate. We demonstrate the existence of an ignition threshold for weak interactions at around ~10 -3 M⊙ s -1 for typical parameters by means of analytic calculations and numerical simulations. We find a degenerate, self-regulated, neutrino-cooled regime above the threshold and an advection-dominated regime below the threshold. Excess heating in the absence of neutrino cooling below the threshold leads to ≳60% of the initial disk mass being ejected in outflows, with typical velocities of ~(0.1–0.2)c, compared to ≲40% at ~(0.1–0.15)c above the threshold. While disks below the threshold show suppressed production of light r-process elements, disks above the threshold can produce the entire range of r-process elements, in good agreement with the observed solar system abundances. Disks below the ignition threshold may produce an overabundance of actinides seen in actinide-boost stars. As gravitational-wave detectors start to sample the neutron star merger parameter space, different disk realizations may be observable via their associated kilonova emission.

79 ASTRONOMY AND ASTROPHYSICS↗

Implementation of Monte Carlo Transport in the General Relativistic SpEC Code

Neutrino transport and neutrino-matter interactions are known to play an important role in the evolution of neutron star mergers and of their post-merger remnants. Neutrinos cool remnants, drive post-merger winds, and deposit energy in the low-density polar regions where relativistic jets may eventually form. Neutrinos also modify the composition of the ejected material, impacting the outcome of nucleosynthesis in merger outflows and the properties of the optical/infrared transients that they power (kilonovae). So far, merger simulations have largely relied on approximate treatments of the neutrinos (leakage, moments) that simplify the equations of radiation transport in a way that makes simulations more affordable but also introduces unquantifiable errors in the results. To improve on these methods, we recently published a first simulation of neutron star mergers using a low-cost Monte Carlo algorithm for neutrino radiation transport. Our transport code limits costs in optically thick regions by placing a hard ceiling on the value of the absorption opacity of the fluid, yet all approximations made within the code are designed to vanish in the limit of infinite numerical resolution. We provide here an in-depth description of this algorithm, of its implementation in the SpEC merger code, and of the expected impact of our approximations in optically thick regions. We argue that the last is a subdominant source of error at the accuracy reached by current simulations and for the interactions currently included in our code. We also provide tests of the most important features of this code.

79 ASTRONOMY AND ASTROPHYSICS↗

Nuclear Physics Multimessenger Astrophysics Constraints on the Neutron Star Equation of State: Adding NICER’s PSR J0740+6620 Measurement

In the past few years, new observations of neutron stars (NSs) and NS mergers have provided a wealth of data that allow one to constrain the equation of state (EOS) of nuclear matter at densities above nuclear saturation density. However, most observations were based on NSs with masses of about 1.4 M ⊙ , probing densities up to ~three to four times the nuclear saturation density. Even higher densities are probed inside massive NSs such as PSR J0740+6620. Very recently, new radio observations provided an update to the mass estimate for PSR J0740+6620, and X-ray observations by the NICER and XMM telescopes constrained its radius. Based on these new measurements, we revisit our previous nuclear physics multimessenger astrophysics constraints and derive updated constraints on the EOS describing the NS interior. By combining astrophysical observations of two radio pulsars, two NICER measurements, the two gravitational-wave detections GW170817 and GW190425, detailed modeling of the kilonova AT 2017gfo, and the gamma-ray burst GRB 170817A, we are able to estimate the radius of a typical 1.4 M ⊙ NS to be ${11.94}_{-0.87}^{+0.76}\,\mathrm{km}$ at 90% confidence. Finally, our analysis allows us to revisit the upper bound on the maximum mass of NSs and disfavors the presence of a strong first-order phase transition from nuclear matter to exotic forms of matter, such as quark matter, inside NSs.

60 APPLIED LIFE SCIENCES↗

Hubble Space Telescope Observations of GW170817: Complete Light Curves and the Properties of the Galaxy Merger of NGC 4993

We present the complete set of Hubble Space Telescope imaging of the binary neutron star merger GW170817 and its optical counterpart AT 2017gfo. Including deep template imaging in F814W, F110W, F140W, and F160W at 3.4 yr post-merger, we reanalyze the full light curve of AT 2017gfo across 12 bands from 5 to 1273 rest-frame days after merger. We obtain four new detections of the short γ-ray burst 170817A afterglow from 109 to 170 rest-frame days post-merger. These detections are consistent with the previously observed β = -0.6 spectral index in the afterglow light curve with no evidence for spectral evolution. We also analyze our limits in the context of kilonova afterglow or IR dust echo emission but find that our limits are not constraining for these models. We use the new data to construct deep optical and IR stacks, reaching limits of M = -6.3 to -4.6 mag, to analyze the local environment around AT 2017gfo and low surface brightness features in its host galaxy NGC 4993. We rule out the presence of any globular cluster at the position of AT 2017gfo to 2.3 × 10$^{4}$L$_{⊙}$, including those with the reddest V - H colors. Finally, we analyze the substructure of NGC 4993 in deep residual imaging and find shell features that extend up to 71."8 (14.2 kpc) from NGC 4993. The shells have a cumulative stellar mass of 6.3 × 10$^{8}$M$_{⊙}$, roughly 2% of NGC 4993, and mass-weighted ages of >3 Gyr. We conclude that it was unlikely that the GW170817 progenitor system formed in the galaxy merger.

79 ASTRONOMY AND ASTROPHYSICS↗

Proposed Lunar Measurements of r-Process Radioisotopes to Distinguish the Origin of Deep-sea 244 Pu

Abstract 244 Pu has recently been discovered in deep-sea deposits spanning the past 10 Myr, a period that includes two 60 Fe pulses from nearby supernovae. 244 Pu is among the heaviest r -process products, and we consider whether it was created in supernovae, which is disfavored by nucleosynthesis simulations, or in an earlier kilonova event that seeded the nearby interstellar medium with 244 Pu that was subsequently swept up by the supernova debris. We discuss how these possibilities can be probed by measuring 244 Pu and other r -process radioisotopes such as 129 I and 182 Hf, both in lunar regolith samples returned to Earth by missions such as Chang’e and Artemis, and in deep-sea deposits.

79 ASTRONOMY AND ASTROPHYSICS↗

Bound Debris Expulsion from Neutron Star Merger Remnants

Abstract Many studies have found that neutron star mergers leave a fraction of the stars’ mass in bound orbits surrounding the resulting massive neutron star or black hole. This mass is a site ofr-process nucleosynthesis and can generate a wind that contributes to a kilonova. However, comparatively little is known about the dynamics determining its mass or initial structure. Here we begin to investigate these questions, starting with the origin of the disk mass. Using tracer particle as well as discretized fluid data from numerical simulations, we identify where in the neutron stars the debris came from, the paths it takes in order to escape from the neutron stars’ interiors, and the times and locations at which its orbital properties diverge from those of neighboring fluid elements that end up remaining in the merged neutron star.

Astronomy & Astrophysics↗

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↗

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↗

The Impact of Observing Strategy on Cosmological Constraints with LSST

The generation-defining Vera C. Rubin Observatory will make state-of-the-art measurements of both the static and transient universe through its Legacy Survey for Space and Time (LSST). With such capabilities, it is immensely challenging to optimize the LSST observing strategy across the survey's wide range of science drivers. Many aspects of the LSST observing strategy relevant to the LSST Dark Energy Science Collaboration, such as survey footprint definition, single-visit exposure time, and the cadence of repeat visits in different filters, are yet to be finalized. Here, we present metrics used to assess the impact of observing strategy on the cosmological probes considered most sensitive to survey design; these are large-scale structure, weak lensing, type Ia supernovae, kilonovae, and strong lens systems (as well as photometric redshifts, which enable many of these probes). We evaluate these metrics for over 100 different simulated potential survey designs. Our results show that multiple observing strategy decisions can profoundly impact cosmological constraints with LSST; these include adjusting the survey footprint, ensuring repeat nightly visits are taken in different filters, and enforcing regular cadence. We provide public code for our metrics, which makes them readily available for evaluating further modifications to the survey design. We conclude with a set of recommendations and highlight observing strategy factors that require further research.

79 ASTRONOMY AND ASTROPHYSICS↗

The R -process Alliance: A Nearly Complete R -process Abundance Template Derived from Ultraviolet Spectroscopy of the R -process-enhanced Metal-poor Star HD 222925 *

We present a nearly complete rapid neutron-capture process (r-process) chemical inventory of the metal-poor ([Fe/H] = -1.46 ± 0.10) r-process-enhanced ([Eu/Fe] = +1.32 ± 0.08) halo star HD 222925. This abundance set is the most complete for any object beyond the solar system, with a total of 63 metals detected and seven with upper limits. It comprises 42 elements from 31 ≤ Z ≤ 90, including elements rarely detected in r-process-enhanced stars, such as Ga, Ge, As, Se, Cd, In, Sn, Sb, Te, W, Re, Os, Ir, Pt, and Au. We derive these abundances from an analysis of 404 absorption lines in ultraviolet spectra collected using the Space Telescope Imaging Spectrograph on the Hubble Space Telescope and previously analyzed optical spectra. A series of appendices discusses the atomic data and quality of fits for these lines. The r-process elements from Ba to Pb, including all elements at the third r-process peak, exhibit remarkable agreement with the solar r-process residuals, with a standard deviation of the differences of only 0.08 dex (17%). In contrast, deviations among the lighter elements from Ga to Te span nearly 1.4 dex, and they show distinct trends from Ga to Se, Nb through Cd, and In through Te. The r-process contribution to Ga, Ge, and As is small, and Se is the lightest element whose production is dominated by the r-process. The lanthanide fraction, log X La = -1.39 ± 0.09, is typical for r-process-enhanced stars and higher than that of the kilonova from the GW170817 neutron-star merger event. We advocate adopting this pattern as an alternative to the solar r-process-element residuals when confronting future theoretical models of heavy-element nucleosynthesis with observations.

79 ASTRONOMY AND ASTROPHYSICS↗

Results of the Photometric LSST Astronomical Time-series Classification Challenge (PLAsTiCC)

Abstract Next-generation surveys like the Legacy Survey of Space and Time (LSST) on the Vera C. Rubin Observatory (Rubin) will generate orders of magnitude more discoveries of transients and variable stars than previous surveys. To prepare for this data deluge, we developed the Photometric LSST Astronomical Time-series Classification Challenge (PLAsTiCC), a competition that aimed to catalyze the development of robust classifiers under LSST-like conditions of a nonrepresentative training set for a large photometric test set of imbalanced classes. Over 1000 teams participated in PLAsTiCC, which was hosted in the Kaggle data science competition platform between 2018 September 28 and 2018 December 17, ultimately identifying three winners in 2019 February. Participants produced classifiers employing a diverse set of machine-learning techniques including hybrid combinations and ensemble averages of a range of approaches, among them boosted decision trees, neural networks, and multilayer perceptrons. The strong performance of the top three classifiers on Type Ia supernovae and kilonovae represent a major improvement over the current state of the art within astronomy. This paper summarizes the most promising methods and evaluates their results in detail, highlighting future directions both for classifier development and simulation needs for a next-generation PLAsTiCC data set.

79 ASTRONOMY AND ASTROPHYSICS↗

Constraints on the Dense Matter Equation of State and Neutron Star Properties from NICER’s Mass–Radius Estimate of PSR J0740+6620 and Multimessenger Observations

In recent years our understanding of the dense matter equation of state (EOS) of neutron stars has significantly improved by analyzing multimessenger data from radio/X-ray pulsars, gravitational wave events, and from nuclear physics constraints. Here we study the additional impact on the EOS from the jointly estimated mass and radius of PSR J0740+6620, presented in Riley et al. by analyzing a combined data set from X-ray telescopes NICER and XMM-Newton. We employ two different high-density EOS parameterizations: a piecewise-polytropic (PP) model and a model based on the speed of sound in a neutron star (CS). Here, at nuclear densities these are connected to microscopic calculations of neutron matter based on chiral effective field theory (EFT) interactions. In addition to the new NICER data for this heavy neutron star, we separately study constraints from the radio timing mass measurement of PSR J0740+6620, the gravitational wave events of binary neutron stars GW190425 and GW170817, and for the latter the associated kilonova AT2017gfo. By combining all these, and the NICER mass-radius estimate of PSR J0030+0451, we find the radius of a 1.4 M ⊙ neutron star to be constrained to the 95% credible ranges 12.33$^{+0.76}_{–0.81}$km (PP model) and 12.18$^{+0.56}_{–0.79}$km (CS model). In addition, we explore different chiral EFT calculations and show that the new NICER results provide tight constraints for the pressure of neutron star matter at around twice saturation density, which shows the power of these observations to constrain dense matter interactions at intermediate densities.

79 ASTRONOMY AND ASTROPHYSICS↗

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↗

Hours-long Near-UV/Optical Emission from Mildly Relativistic Outflows in Black Hole–Neutron Star Mergers

The ongoing LIGO–Virgo–KAGRA observing run O4 provides an opportunity to discover new multimessenger events, including binary neutron star (BNS) mergers such as GW170817 and the highly anticipated first detection of a multimessenger black hole–neutron star (BH–NS) merger. While BNS mergers were predicted to exhibit early optical emission from mildly relativistic outflows, it has remained uncertain whether the BH–NS merger ejecta provides the conditions for similar signals to emerge. We present the first modeling of early near-ultraviolet/optical emission from mildly relativistic outflows in BH–NS mergers. Adopting optimal binary properties, a mass ratio of q = 2, and a rapidly rotating BH, we utilize numerical relativity and general relativistic magnetohydrodynamic (GRMHD) simulations to follow the binary's evolution from premerger to homologous expansion. We use an M1 neutrino transport GRMHD simulation to self-consistently estimate the opacity distribution in the outflows and find a bright near-ultraviolet/optical signal that emerges due to jet-powered cocoon cooling emission, outshining the kilonova emission at early time. The signal peaks at an absolute magnitude of ~–15 a few hours after the merger, longer than previous estimates, which did not consider the first principles–based jet launching. By late 2024, the Rubin Observatory will have the capability to track the entire signal evolution or detect its peak up to distances of ≳1 Gpc. In 2026, ULTRASAT will conduct all-sky surveys within minutes, detecting some of these events within ~200 Mpc. The BH–NS mergers with higher mass ratios or lower BH spins would produce shorter and fainter signals.

79 ASTRONOMY AND ASTROPHYSICS↗