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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

All sources NASA NTRS DOE OSTI

DOE OSTI2017

Infrared Emission from Kilonovae: The Case of the Nearby Short Hard Burst GRB 160821B

In this paper, we present constraints on Ks-band emission from one of the nearest short hard gamma-ray bursts, GRB 160821B, at z = 0.16, at three epochs. We detect a red relativistic afterglow from the jetted emission in the first epoch but do not detect any excess kilonova emission in the second two epochs. We compare upper limits obtained with Keck I/MOSFIRE to multi-dimensional radiative transfer models of kilonovae, that employ composition-dependent nuclear heating and LTE opacities of heavy elements. We discuss eight models that combine toroidal dynamical ejecta and two types of wind and one model with dynamical ejecta only. We also discuss simple, empirical scaling laws of predicted emission as a function of ejecta mass and ejecta velocity. Our limits for GRB 160821B constrain the ejecta mass to be lower than 0.03 M ⊙ for velocities greater than 0.1 c. At the distance sensitivity range of advanced LIGO, similar ground-based observations would be sufficiently sensitive to the full range of predicted model emission including models with only dynamical ejecta. Finally, the color evolution of these models shows that I–K color spans 7–16 mag, which suggests that even relatively shallow infrared searches for kilonovae could be as constraining as optical searches.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2017

Toward Rapid Transient Identification and Characterization of Kilonovae

With the increasing sensitivity of advanced gravitational-wave (GW) detectors, the first joint detection of an electromagnetic and GW signal from a compact binary merger will hopefully happen within this decade. However, current GW likelihood sky areas span ∼100--1000 deg{sup 2}, and thus it is a challenging task to identify which, if any, transient corresponds to the GW event. In this study, we make a comparison between recent kilonova/macronova light-curve models for the purpose of assessing potential light-curve templates for counterpart identification. We show that recent analytical and parameterized models for these counterparts result in qualitative agreement with more complicated radiative transfer simulations. Our analysis suggests that with improved light-curve models with smaller uncertainties it will become possible to extract information about ejecta properties and binary parameters directly from the light-curve measurement. Even tighter constraints are obtained in cases for which GW and kilonova parameter estimation results are combined. It will therefore be important to make comparisons and potentially combine parameter estimation with the kilonova and GW results. However, to be prepared for upcoming detections, more realistic kilonova models are needed. These will require numerical relativity with more detailed microphysics, better radiative transfer simulations, and a better understanding of the underlying nuclear physics.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2017

A Search for Kilonovae in the Dark Energy Survey

The coalescence of a binary neutron star pair is expected to produce gravitational waves (GW) and electromagnetic radiation, both of which may be detectable with currently available instruments. We describe a search for a predicted r-process optical transient from these mergers, dubbed the “kilonova” (KN), using griz broadband data from the Dark Energy Survey Supernova Program (DES-SN). Some models predict KNe to be redder, shorter-lived, and dimmer than supernovae (SNe), but the event rate of KNe is poorly constrained. We simulate KN and SN light curves with the Monte-Carlo simulation code SNANA to optimize selection requirements, determine search efficiency, and predict SN backgrounds. Our analysis of the first two seasons of DES-SN data results in 0 events, and is consistent with our prediction of 1.1 ± 0.2 background events based on simulations of SNe. From our prediction, there is a 33% chance of finding 0 events in the data. Assuming no underlying galaxy flux, our search sets 90% upper limits on the KN volumetric rate of 1.0 ×10{sup 7} Gpc{sup −3} yr{sup −1} for the dimmest KN model we consider (peak i-band absolute magnitude M{sub i}=−11.4 mag) and 2.4 ×10{sup 4} Gpc{sup −3} yr{sup −1} for the brightest (M{sub i}=−16.2 mag). Accounting for anomalous subtraction artifacts on bright galaxies, these limits are ∼3 times higher. This analysis is the first untriggered optical KN search and informs selection requirements and strategies for future KN searches. Our upper limits on the KN rate are consistent with those measured by GW and gamma-ray burst searches.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2017

The Electromagnetic Counterpart of the Binary Neutron Star Merger LIGO/Virgo GW170817. III. Optical and UV Spectra of a Blue Kilonova from Fast Polar Ejecta

We present optical and ultraviolet spectra of the first electromagnetic counterpart to a gravitational-wave (GW) source, the binary neutron star merger GW170817. Spectra were obtained nightly between 1.5 and 9.5 days post-merger, using the Southern Astrophysical Research and Magellan telescopes; the UV spectrum was obtained with the Hubble Space Telescope at 5.5 days. Our data reveal a rapidly fading blue component (T≈5500 K at 1.5 days) that quickly reddens; spectra later than ≳4.5 days peak beyond the optical regime. The spectra are mostly featureless, although we identify a possible weak emission line at ∼7900 Å at t≲4.5 days. The colors, rapid evolution, and featureless spectrum are consistent with a “blue” kilonova from polar ejecta comprised mainly of light r-process nuclei with atomic mass number A≲140. This indicates a sightline within θ{sub obs}≲45{sup ∘} of the orbital axis. Comparison to models suggests ∼0.03 M {sub ⊙} of blue ejecta, with a velocity of ∼0.3c. The required lanthanide fraction is ∼10{sup −4}, but this drops to <10{sup −5} in the outermost ejecta. The large velocities point to a dynamical origin, rather than a disk wind, for this blue component, suggesting that both binary constituents are neutron stars (as opposed to a binary consisting of a neutron star and a black hole). For dynamical ejecta, the high mass favors a small neutron star radius of ≲12 km. This mass also supports the idea that neutron star mergers are a major contributor to r-process nucleosynthesis.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2017

The Electromagnetic Counterpart of the Binary Neutron Star Merger LIGO/Virgo GW170817. VI. Radio Constraints on a Relativistic Jet and Predictions for Late-time Emission from the Kilonova Ejecta

We present Very Large Array (VLA) and Atacama Large Millimeter/submillimeter Array (ALMA) radio observations of GW170817, the first Laser Interferometer Gravitational-wave Observatory (LIGO)/Virgo gravitational wave (GW) event from a binary neutron star merger and the first GW event with an electromagnetic (EM) counterpart. Our data include the first observations following the discovery of the optical transient at both the centimeter (13.7 hr post-merger) and millimeter (2.41 days post-merger) bands. We detect faint emission at 6 GHz at 19.47 and 39.23 days after the merger, but not in an earlier observation at 2.46 days. We do not detect cm/mm emission at the position of the optical counterpart at frequencies of 10–97.5 GHz at times ranging from 0.6 to 30 days post-merger, ruling out an on-axis short gamma-ray burst (SGRB) for energies ≳10{sup 48} erg. For fiducial SGRB parameters, our limits require an observer viewer angle of ≳20°. The radio and X-ray data can be jointly explained as the afterglow emission from an SGRB with a jet energy of ∼10{sup 49}--10{sup 50} erg that exploded in a uniform density environment with n∼10{sup −4}--10{sup −2} cm{sup −3}, viewed at an angle of ∼20°–40° from the jet axis. Using the results of our light curve and spectral modeling, in conjunction with the inference of the circumbinary density, we predict the emergence of late-time radio emission from the deceleration of the kilonova (KN) ejecta on a timescale of ∼5–10 years that will remain detectable for decades with next-generation radio facilities, making GW170817 a compelling target for long-term radio monitoring.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2017

The Electromagnetic Counterpart of the Binary Neutron Star Merger LIGO/Virgo GW170817. II. UV, Optical, and Near-infrared Light Curves and Comparison to Kilonova Models

We present UV, optical, and near-infrared (NIR) photometry of the first electromagnetic counterpart to a gravitational wave source from Advanced Laser Interferometer Gravitational-wave Observatory (LIGO)/Virgo, the binary neutron star merger GW170817. Our data set extends from the discovery of the optical counterpart at 0.47–18.5 days post-merger, and includes observations with the Dark Energy Camera (DECam), Gemini-South/FLAMINGOS-2 (GS/F2), and the Hubble Space Telescope (HST). The spectral energy distribution (SED) inferred from this photometry at 0.6 days is well described by a blackbody model with T≈8300 K, a radius of R≈4.5×10{sup 14} cm (corresponding to an expansion velocity of v≈0.3c), and a bolometric luminosity of L{sub bol}≈5×10{sup 41} erg s{sup −1}. At 1.5 days we find a multi-component SED across the optical and NIR, and subsequently we observe rapid fading in the UV and blue optical bands and significant reddening of the optical/NIR colors. Modeling the entire data set, we find that models with heating from radioactive decay of {sup 56}Ni, or those with only a single component of opacity from r-process elements, fail to capture the rapid optical decline and red optical/NIR colors. Instead, models with two components consistent with lanthanide-poor and lanthanide-rich ejecta provide a good fit to the data; the resulting “blue” component has M{sub ej}{sup blue}≈0.01 M{sub ⊙} and v{sub ej}{sup blue}≈0.3 c, and the “red” component has M{sub ej}{sup red}≈0.04 M{sub ⊙} and v{sub ej}{sup red}≈0.1 c. These ejecta masses are broadly consistent with the estimated r-process production rate required to explain the Milky Way r-process abundances, providing the first evidence that binary neutron star (BNS) mergers can be a dominant site of r-process enrichment.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Kilonova Luminosity Function Constraints Based on Zwicky Transient Facility Searches for 13 Neutron Star Merger Triggers during O3

We present a systematic search for optical counterparts to 13 gravitational wave (GW) triggers involving at least one neutron star during LIGO/Virgo's third observing run (O3). We searched binary neutron star (BNS) and neutron star black hole (NSBH) merger localizations with the Zwicky Transient Facility (ZTF) and undertook follow-up with the Global Relay of Observatories Watching Transients Happen (GROWTH) collaboration. The GW triggers had a median localization area of 4480 deg 2 , a median distance of 267 Mpc, and false-alarm rates ranging from 1.5 to 10 -25 yr -1 . The ZTF coverage in the g and r bands had a median enclosed probability of 39%, median depth of 20.8 mag, and median time lag between merger and the start of observations of 1.5 hr. The O3 follow-up by the GROWTH team comprised 340 UltraViolet/Optical/InfraRed (UVOIR) photometric points, 64 OIR spectra, and three radio images using 17 different telescopes. We find no promising kilonovae (radioactivity-powered counterparts), and we show how to convert the upper limits to constrain the underlying kilonova luminosity function. Initially, we assume that all GW triggers are bona fide astrophysical events regardless of false-alarm rate and that kilonovae accompanying BNS and NSBH mergers are drawn from a common population; later, we relax these assumptions. Assuming that all kilonovae are at least as luminous as the discovery magnitude of GW170817 (-16.1 mag), we calculate that our joint probability of detecting zero kilonovae is only 4.2%. If we assume that all kilonovae are brighter than-16.6 mag (the extrapolated peak magnitude of GW170817) and fade at a rate of 1 mag day -1 (similar to GW170817), the joint probability of zero detections is 7%. If we separate the NSBH and BNS populations based on the online classifications, the joint probability of zero detections, assuming all kilonovae are brighter than -16.6 mag, is 9.7% for NSBH and 7.9% for BNS mergers. Moreover, no more than <57% (<89%) of putative kilonovae could be brighter than -16.6 mag assuming flat evolution (fading by 1 mag day -1 ) at the 90% confidence level. If we further take into account the online terrestrial probability for each GW trigger, we find that no more than <68% of putative kilonovae could be brighter than -16.6 mag. Comparing to model grids, we find that some kilonovae must have M ej < 0.03 M ⊙ , X lan > 10 -4 , or φ > 30° to be consistent with our limits. We look forward to searches in the fourth GW observing run; even 17 neutron star mergers with only 50% coverage to a depth of-16 mag would constrain the maximum fraction of bright kilonovae to <25%.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2019

A luminosity distribution for kilonovae based on short gamma-ray burst afterglows

The combined detection of a gravitational-wave signal, kilonova, and short gamma-ray burst (sGRB) from GW170817 marked a scientific breakthrough in the field of multimessenger astronomy. But even before GW170817, there have been a number of sGRBs with possible associated kilonova detections. In this work, we re-examine these ‘historical’ sGRB afterglows with a combination of state-of-the-art afterglow and kilonova models. This allows us to include optical/near-infrared synchrotron emission produced by the sGRB as well as ultraviolet/optical/near-infrared emission powered by the radioactive decay of r-process elements (i.e. the kilonova). Fitting the light curves, we derive the velocity and the mass distribution as well as the composition of the ejected material. The posteriors on kilonova parameters obtained from the fit were turned into distributions for the peak magnitude of the kilonova emission in different bands and the time at which this peak occurs. From the sGRB with an associated kilonova, we found that the peak magnitude in H bands falls in the range [-16.2, -13.1] (95 percent of confidence) and occurs within 0.8-3.6d after the sGRB prompt emission. In g band instead we obtain a peak magnitude in range [-16.8, -12.3] occurring within the first 18 h after the sGRB prompt. From the luminosity distributions of GW170817/AT2017gfo, kilonova candidates GRB130603B, GRB050709, and GRB060614 (with the possible inclusion of GRB150101B, GRB050724A, GRB061201, GRB080905A, GRB150424A, and GRB160821B) and the upper limits from all the other sGRBs not associated with any kilonova detection we obtain for the first time a kilonova luminosity distribution in different bands.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Kilonova Detectability with Wide-field Instruments

Abstract Kilonovae are ultraviolet, optical, and infrared transients powered by the radioactive decay of heavy elements following a neutron star merger. Joint observations of kilonovae and gravitational waves can offer key constraints on the source of Galactic r -process enrichment, among other astrophysical topics. However, robust constraints on heavy element production require rapid kilonova detection (within ∼1 day of merger) as well as multiwavelength observations across multiple epochs. In this study, we quantify the ability of 13 wide-field-of-view instruments to detect kilonovae, leveraging a large grid of over 900 radiative transfer simulations with 54 viewing angles per simulation. We consider both current and upcoming instruments, collectively spanning the full kilonova spectrum. The Roman Space Telescope has the highest redshift reach of any instrument in the study, observing kilonovae out to z ∼ 1 within the first day post-merger. We demonstrate that BlackGEM, DECam, GOTO, the Vera C. Rubin Observatory’s LSST, ULTRASAT, VISTA, and WINTER can observe some kilonovae out to z ∼ 0.1 (∼475 Mpc), while DDOTI, MeerLICHT, PRIME, Swift/UVOT, and ZTF are confined to more nearby observations. Furthermore, we provide a framework to infer kilonova ejecta properties following nondetections and explore variation in detectability with these ejecta parameters.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2017

Lanthanides or Dust in Kilonovae: Lessons Learned from GW170817

The unprecedented optical and near-infrared lightcurves of the first electromagnetic counterpart to a gravitational-wave source, GW170817, a binary neutron star merger, exhibited a strong evolution from blue to near-infrared (a so-called “kilonova” or “macronova”). The emerging near-infrared component is widely attributed to the formation of r-process elements that provide the opacity to shift the blue light into the near-infrared. An alternative scenario is that the light from the blue component gets extinguished by dust formed by the kilonova and subsequently is re-emitted at near-infrared wavelengths. We test here this hypothesis using the lightcurves of AT 2017gfo, the kilonova accompanying GW170817. We find that of the order of 10{sup −5} M{sub ⊙} of carbon is required to reproduce the optical/near-infrared lightcurves as the kilonova fades. This putative dust cools from ∼2000 K at ∼4 days after the event to ∼1500 K over the course of the following week, thus requiring dust with a high condensation temperature, such as carbon. We contrast this with the nucleosynthetic yields predicted by a range of kilonova wind models. These suggest that at most 10{sup −9} M{sub ⊙} of carbon is formed. Moreover, the decay in the inferred dust temperature is slower than that expected in kilonova models. We therefore conclude that in current models of the blue component of the kilonova, the near-infrared component in the kilonova accompanying GW170817 is unlikely to be due to dust.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

Impact of in situ nuclear networks and atomic opacities on neutron star merger ejecta dynamics, nucleosynthesis, and kilonovae

Context. Binary neutron star merger (BNSM) ejecta are key sites of rapid neutron capture (r-process) nucleosynthesis and they produce kilonovae powered by the radioactive decay of freshly synthesized nuclei. Modeling their evolution requires multi-physics simulations involving hydrodynamics, nuclear reactions, and radiative processes. The impact of nuclear burning and atomic opacity is poorly understood and often treated with simplified prescriptions. Aims. We systematically investigate different treatments of nuclear heating, particle thermalization, and atomic opacities in radiation-hydrodynamics simulations of BNSM ejecta and kilonova light curves. Methods. Ejecta profiles from long-term numerical-relativity simulations of asymmetric neutron star binaries with a massive neutron star remnant were evolved to ∼30 days using a 2D ray-by-ray approach. We compared simplified heating-rate and thermalization prescriptions with in situ Nuclear reaction Network (NN) calculations that track nuclear energy deposition and include a composition-dependent thermalization scheme. We also contrasted various gray opacity models with a frequency-dependent treatment based on atomic calculations. Results. Coupling NN and hydrodynamics significantly affects nucleosynthesis and kilonova emission. Assuming homologous expansion alters abundance evolution and produces a narrower, less populated second r-process peak and a third peak shifted to higher mass numbers. The back-reaction of nuclear heating affects the temperature evolution enough to delay and redden the early (t∼ hours) kilonova peaks. A constant thermalization efficiency underestimates and reddens the early emission while overestimating the late-time luminosity compared to the composition-dependent treatment. Analytical opacity prescriptions yield a more extended, colder photosphere, resulting in dimmer, redder kilonovae at early times (t≲ hour), while the delayed recession of the photosphere prolongs the red emission at t ≳ 5 days. Conclusions. Coupling hydrodynamics to an in situ NN is crucial for reliable nucleosynthesis and kilonova predictions. Resolving the first several hundred milliseconds of the hydrodynamics is essential for robust nucleosynthesis calculations. Composition-dependent thermalization and frequency-dependent, atomic-physics-based opacities are needed to accurately capture the temperature evolution of the ejecta and the brightness and color evolution of the kilonova. Calibrated analytic nuclear-power fits with simplified thermalization and opacity prescriptions can still reproduce the density and temperature evolution of the ejecta.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2016

JAMES WEBB SPACE TELESCOPE CAN DETECT KILONOVAE IN GRAVITATIONAL WAVE FOLLOW-UP SEARCH

Kilonovae represent an important electromagnetic counterpart for compact binary mergers, which could become the most commonly detected gravitational-wave (GW) source. Follow-up observations of kilonovae, triggered by GW events, are nevertheless difficult due to poor localization by GW detectors and due to their faint near-infrared peak emission, which has limited observational capability. We show that the Near-Infrared Camera (NIRCam) on the James Webb Space Telescope will be able to detect kilonovae within the relevant GW-detection range of ∼200 Mpc in short (≲12-s) exposure times for a week following the merger. Despite this sensitivity, a kilonova search fully covering a fiducial localized area of 10 deg{sup 2} will not be viable with NIRCam due to its limited field of view. However, targeted surveys may be developed to optimize the likelihood of discovering kilonovae efficiently within limited observing time. We estimate that a survey of 10 deg{sup 2} focused on galaxies within 200 Mpc would require about 13 hr, dominated by overhead times; a survey further focused on galaxies exhibiting high star formation rates would require ∼5 hr. The characteristic time may be reduced to as little as ∼4 hr, without compromising the likelihood of detecting kilonovae, by surveying sky areas associated with 50%, rather than 90%, confidence regions of 3 GW events, rather than a single event. Upon the detection and identification of a kilonova, a limited number of NIRCam follow-up observations could constrain the properties of matter ejected by the binary and the equation of state of dense nuclear matter.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

Radioactivity and thermalization in the ejecta of compact object mergers and their impact on kilonova light curves

One promising electromagnetic signature of compact object mergers are kilonovae: approximately isotropic radioactively powered transients that peak days to weeks post-merger. Key uncertainties in kilonova modeling include the emission profiles of the radioactive decay products—non-thermal β-particles, α-particles, fission fragments, and γ-rays—and the efficiency with which their kinetic energy is absorbed by the ejecta. The radioactive energy emitted, along with its thermalization efficiency, sets the luminosity budget and is therefore crucial for predicting kilonova light curves. We outline uncertainties in the radioactivity, describe the processes by which the decay products transfer energy to the ejecta, and calculate time-dependent thermalization efficiencies for each particle type. We determine the net thermalization efficiency and explore its dependence on r-process yields—in particular, the production of α-decaying translead nuclei—and on ejecta mass, velocity, and magnetic fields. We incorporate our results into detailed radiation transport simulations, and calculate updated kilonova light curve predictions. Thermalization effects reduce kilonova luminosities by a factor of roughly 2 at peak, and by an order of magnitude at later times (15 days or more after explosion). We present analytic fits to time-dependent thermalization efficiencies, which can be used to improve light curve models. We revisit the putative kilonova that accompanied gamma-ray burst 130603B, and estimate the mass ejected in that event. We find later time kilonova light curves can be significantly impacted by α-decay from translead isotopes; data at these times may therefore be diagnostic of ejecta abundances.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Designing an Optimal Kilonova Search Using DECam for Gravitational-wave Events

We address the problem of optimally identifying all kilonovae detected via gravitational-wave emission in the upcoming LIGO/Virgo/KAGRA observing run, O4, which is expected to be sensitive to a factor of ~7 more binary neutron star (BNS) alerts than previously. Electromagnetic follow-up of all but the brightest of these new events will require >1 m telescopes, for which limited time is available. We present an optimized observing strategy for the DECam during O4. We base our study on simulations of gravitational-wave events expected for O4 and wide-prior kilonova simulations. We derive the detectabilities of events for realistic observing conditions. We optimize our strategy for confirming a kilonova while minimizing telescope time. For a wide range of kilonova parameters, corresponding to a fainter kilonova compared to GW170817/AT 2017gfo, we find that, with this optimal strategy, the discovery probability for electromagnetic counterparts with the DECam is ~80% at the nominal BNS gravitational-wave detection limit for O4 (190 Mpc), which corresponds to an ~30% improvement compared to the strategy adopted during the previous observing run. For more distant events (~330 Mpc), we reach an ~60% probability of detection, a factor of ~2 increase. For a brighter kilonova model dominated by the blue component that reproduces the observations of GW170817/AT 2017gfo, we find that we can reach ~90% probability of detection out to 330 Mpc, representing an increase of ~20%, while also reducing the total telescope time required to follow up events by ~20%.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2023

Modelling populations of kilonovae

Abstract The 2017 detection of a kilonova coincident with gravitational-wave emission has identified neutron star mergers as the major source of the heaviest elements and dramatically constrained alternative theories of gravity. Observing a population of such sources has the potential to transform cosmology, nuclear physics, and astrophysics. However, with only one confident multi-messenger detection currently available, modelling the diversity of signals expected from such a population requires improved theoretical understanding. In particular, models that are quick to evaluate and are calibrated with more detailed multi-physics simulations are needed to design observational strategies for kilonovae detection and to obtain rapid-response interpretations of new observations. We use grey-opacity models to construct populations of kilonovae, spanning ejecta parameters predicted by numerical simulations. Our modelling focuses on wavelengths relevant for upcoming optical surveys, such as the Rubin Observatory Legacy Survey of Space and Time (LSST). In these simulations, we implement heating rates that are based on nuclear reaction network calculations. We create a Gaussian-process emulator for kilonova grey opacities, calibrated with detailed radiative transfer simulations. Using recent fits to numerical relativity simulations, we predict how the ejecta parameters from binary neutron star (BNS) mergers shape the population of kilonovae, accounting for the viewing-angle dependence. Our simulated population of BNS mergers produce peak i-band absolute magnitudes of −20 ≤ Mi ≤ −11. A comparison with detailed radiative transfer calculations indicates that further improvements are needed to accurately reproduce spectral shapes over the full light curve evolution.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Kilonovae Across the Nuclear Physics Landscape: The Impact of Nuclear Physics Uncertainties on r-process-powered Emission

Merging neutron stars produce “kilonovae”—electromagnetic transients powered by the decay of unstable nuclei synthesized via rapid neutron capture (the r-process) in material that is gravitationally unbound during inspiral and coalescence. Kilonova emission, if accurately interpreted, can be used to characterize the masses and compositions of merger-driven outflows, helping to resolve a long-standing debate about the origins of r-process material in the Universe. We explore how the uncertain properties of nuclei involved in the r-process complicate the inference of outflow properties from kilonova observations. Using r-process simulations, we show how nuclear physics uncertainties impact predictions of radioactive heating and element synthesis. For a set of models that span a large range in both predicted heating and final abundances, we carry out detailed numerical calculations of decay product thermalization and radiation transport in a kilonova ejecta with a fixed mass and density profile. The light curves associated with our models exhibit great diversity in their luminosities, with peak brightness varying by more than an order of magnitude. We also find variability in the shape of the kilonova light curves and their color, which in some cases runs counter to the expectation that increasing levels of lanthanide and/or actinide enrichment will be correlated with longer, dimmer, redder emission.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

2D end-to-end modelling of kilonovae from binary neutron star merger remnants

We investigate the kilonova emission resulting from outflows produced in a 3D general-relativistic magnetohydrodynamic (GRMHD) simulation of a hypermassive neutron star (HMNS) remnant. We map the outflows into the flash hydrodynamics code to model their expansion in axisymmetry, and study the effects of employing different r-process heating rates. Except for the highest heating rate prescription, we find no significant differences with respect to overall ejecta dynamics and morphology compared to the simulation without heating. Once homologous expansion is attained, typically after ~2s for these ejecta, we map the outflows to the sedona radiative transfer code and compute the spectral evolution of the kilonova and broad-band light curves in various Legacy Survey of Space and Time (LSST) bands. The kilonova properties depend on the remnant lifetime, with peak luminosities and peak time-scales increasing for longer lived remnants that produce more massive ejecta. For all models, there is a strong dependence of both the bolometric and broad-band light curves on the viewing angle. While the short-lived (12 ms) remnant produces higher luminosities when viewed from angles closer to the pole, longer lived remnants (240 ms and 2.5 s) are more luminous when viewed from angles closer to the equator. Our results highlight the importance of self-consistent, long-term modelling of merger ejecta, and taking viewing-angle dependence into account when interpreting observed kilonova light curves. We find that magnetized outflows from an HMNS – if it survives long enough – could explain blue kilonovae, such as the blue emission seen in AT2017gfo.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Superheavy Elements in Kilonovae

Abstract As LIGO-Virgo-KAGRA enters its fourth observing run, a new opportunity to search for electromagnetic counterparts of compact object mergers will also begin. The light curves and spectra from the first “kilonova” associated with a binary neutron star merger (NSM) suggests that these sites are hosts of the rapid neutron capture (“ r ”) process. However, it is unknown just how robust elemental production can be in mergers. Identifying signposts of the production of particular nuclei is critical for fully understanding merger-driven heavy-element synthesis. In this study, we investigate the properties of very neutron-rich nuclei for which superheavy elements ( Z ≥ 104) can be produced in NSMs and whether they can similarly imprint a unique signature on kilonova light-curve evolution. A superheavy-element signature in kilonovae represents a route to establishing a lower limit on heavy-element production in NSMs as well as possibly being the first evidence of superheavy-element synthesis in nature. Favorable NSM conditions yield a mass fraction of superheavy elements X Z ≥104 ≈ 3 × 10 −2 at 7.5 hr post-merger. With this mass fraction of superheavy elements, we find that the component of kilonova light curves possibly containing superheavy elements may appear similar to those arising from lanthanide-poor ejecta. Therefore, photometric characterizations of superheavy-element rich kilonova may possibly misidentify them as lanthanide-poor events.

79 ASTRONOMY AND ASTROPHYSICS↗