The detailed observations of GW170817 proved for the first time directly that neutron star mergers are a major production site of heavy elements. The observations could be fit by a number of simulations that qualitatively agree, but can quantitatively differ (e.g., in total r-process mass) by an order of magnitude. We categorize kilonova ejecta into several typical morphologies motivated by numerical simulations, and apply a radiative transfer Monte Carlo code to study how the geometric distribution of the ejecta shapes the emitted radiation. We find major impacts on both spectra and light curves. The peak bolometric luminosity can vary by two orders of magnitude and the timing of its peak by a factor of five. These findings provide the crucial implication that the ejecta masses inferred from observations around the peak brightness are uncertain by at least an order of magnitude. Mixed two-component models with lanthanide-rich ejecta are particularly sensitive to geometric distribution. A subset of mixed models shows very strong viewing angle dependence due to lanthanide "curtaining," which persists even if the relative mass of lanthanide-rich component is small. The angular dependence is weak in the rest of our models, but different geometric combinations of the two components lead to a highly diverse set of light curves. We identify geometry-dependent P Cygni features in late spectra that directly map out strong lines in the simulated opacity of neodymium, which can help to constrain the ejecta geometry and to directly probe the r-process abundances.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

EREBUS Coupled Hydrodynamic Radiographic Dynamic Reconstruction & Limited View Reconstructions [Slides]

Abstract not provided.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

A tale of two mergers: constraints on kilonova detection in two short GRBs at z~0.5

Explore the source record for details and available documents.

O'Conner, B.↗

DOE OSTI2020

The Ristra Project: FY20/21 Milestone Report

The ASC Advanced Technology Development and Mitigation (ATDM) sub-program was established in 2014 to develop new simulation tools operating on exascale-class computers to serve NNSA (see Appendix B). Over the course of ATDM, LANL management have set a strategy for exascale-class application codes that follows two supportive and mutually risk-mitigating paths: evolution for established production integrated design codes (IDCs) – with a strong pedigree within the user community – based upon existing programming paradigms(MPI+X); and a new start ATDM project, Ristra, a high-risk/high-reward push for a next-generation multi-physics, multi-scale simulation toolkit based on emerging advanced programming systems(with an initial focus on data-flow task-based models exemplified by Legion). The role of Ristra as the high-risk/high-reward path for LANL’s codes was fully consistent with the goals of ATDM as described in Appendix B, in particular its emphasis on evolving ASC capabilities through novel computing programming models and computing technologies.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

A search for optical and near-infrared counterparts of the compact binary merger GW190814

Here, we report on our observing campaign of the compact binary merger GW190814, detected by the Advanced LIGO and Advanced Virgo detectors on 2019 August 14. This signal has the best localization of any observed gravitational wave (GW) source, with a 90 per cent probability area of 18.5 deg 2 , and an estimated distance of ≈240 Mpc. We obtained wide-field observations with the Deca-Degree Optical Transient Imager (DDOTI) covering 88 per cent of the probability area down to a limiting magnitude of w = 19.9 AB. Nearby galaxies within the high probability region were targeted with the Lowell Discovery Telescope (LDT), whereas promising candidate counterparts were characterized through multicolour photometry with the Reionization and Transients InfraRed (RATIR) and spectroscopy with the Gran Telescopio de Canarias (GTC). We use our optical and near-infrared limits in conjunction with the upper limits obtained by the community to constrain the possible electromagnetic counterparts associated with the merger. A gamma-ray burst seen along its jet’s axis is disfavoured by the multiwavelength data set, whereas the presence of a burst seen at larger viewing angles is not well constrained. Although our observations are not sensitive to a kilonova similar to AT2017gfo, we can rule out high-mass (>0.1 M ⊙ ) fast-moving (mean velocity ≥0.3c) wind ejecta for a possible kilonova associated with this merger.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Composition Effects on Kilonova Spectra and Light Curves. I

The merger of neutron star binaries is believed to eject a wide range of heavy elements into the universe. By observing the emission from this ejecta, scientists can probe the ejecta properties (mass, velocity, and composition distributions). The emission (a.k.a. kilonova) is powered by the radioactive decay of the heavy isotopes produced in the merger and this emission is reprocessed by atomic opacities to optical and infrared wavelengths. Understanding the ejecta properties requires calculating the dependence of this emission on these opacities. The strong lines in the optical and infrared in lanthanide opacities have been shown to significantly alter the light curves and spectra in these wavelength bands, arguing that the emission in these wavelengths can probe the composition of this ejecta. In this work, we study variations in the kilonova emission by varying individual lanthanide (and the actinide uranium) concentrations in the ejecta. The broad forest of lanthanide lines makes it difficult to determine the exact fraction of individual lanthanides. Nd is an exception. Its opacities above 1 μm are higher than other lanthanides and observations of kilonovae can potentially probe increased abundances of Nd. Similarly, at early times when the ejecta is still hot (first day), the U opacity is strong in the 0.2–1 μm wavelength range and kilonova observations may also be able to constrain these abundances.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Time-dependent density-functional theory of high-order harmonic generation from noble-gas atoms driven by orthogonally polarized two-color laser fields

An orthogonally polarized two-color (OTC) field of frequency ratio 1 : 2 is known to generate both odd and even harmonic spectra in the mutually orthogonal directions of the fundamental and the second-harmonic fields, respectively. Recollision angles of an ionized electron driven back to its parent ion by an OTC field are inferred from the intensity ratio of odd and even harmonics, where the relative phase delays between the two-color fields can be used in time-resolved measurements on a subcycle time scale. Here, we present the ab initio, many-electron calculation of time-resolved recollision angles of noble-gas atoms driven by an OTC field, based on the time-dependent density-functional theory. Our results agree well with experimental data by Shafir et al. in Nat. Phys. 5, 412 (2009) of helium atom and New J. Phys. 12, 073032 (2010) of neon atom.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2020

Gamma-Rays from Kilonovae and the Cosmic Gamma-Ray Background

The recent detection of the gravitational-wave event GW170817, produced by the coalescence of two neutron stars, and of its optical–infrared counterpart, powered by the radioactive decay of r-process elements, has opened a new window onto gamma-ray astronomy: the direct detection of photons coming from such decays. In this work, we calculate the contribution of kilonovae to the diffuse gamma-ray background in the MeV range, using recent results on the spectra of the gamma-rays emitted in individual events, and we compare it with that from other sources. We find that the contribution from kilonovae is not dominant in such an energy range, but within current uncertainties, and its addition to other sources might help to fit the observational data.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2020

Gamma Rays from Kilonova: A Potential Probe of r-process Nucleosynthesis

The mergers of compact binaries with at least one neutron star component are the potential leading sites of the production and ejection of r-process elements. Discoveries of galactic binary pulsars, short gamma-ray bursts, and gravitational-wave detections have all been constraining the rate of these events, while the gravitational wave plus broadband electromagnetic coverage of binary neutron star merger (GW170817) has also placed constraints on the properties (mass and composition) of the merger ejecta. But uncertainties and ambiguities in modeling the optical and infrared emission make it difficult to definitively measure the distribution of heavy isotopes in these mergers. In contrast, gamma rays emitted in the decay of these neutron-rich ejecta may provide a more direct measurement of the yields. We calculate the gamma production in remnants of neutron star mergers, considering two epochs: a kilonova epoch, lasting about two weeks, and a much later epoch of tens and hundreds of thousands of years after the merger. For the kilonova epoch, when the expanding ejecta is still only partially transparent to gamma radiation, we use 3D radiative transport simulations to produce the spectra. Futhermore, we show that the gamma-ray spectra associated with beta- and alpha-decay provide a fingerprint of the ejecta properties and, for a sufficiently nearby remnant, may be detectable, even for old remnants. We compare our gamma spectra with the potential detection limits of next generation detectors, including the Lunar Occultation Explorer (LOX), the All-sky Medium Energy Gamma-ray Observatory (AMEGO), and the Compton Spectrometer and Imager (COSI). We show that fission models can be discriminated via the presence of short-lived fission fragments in the remnant spectra.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2019

Nucleosynthesis Probes of Cosmic Explosions [Slides]

GW170817 demonstrated the potential of GW+EM observations of NS mergers: most direct observation of r-process production to date. But what can we really learn from the observations. Getting from the observations to constraining properties of the merger ejecta requires multiple aspects of simulation (merger, disk calculations, and radiation flow) and physics (atomic, nuclear cross-section, dense nuclear matter, neutrino, …). Can we maximize what we can do with observations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗