O'Shaughnessy, R. | Search NASA

DOE OSTI2021

A Broad Grid of 2D Kilonova Emission Models

Depending upon the properties of their compact remnants and the physics included in the models, simulations of neutron star mergers can produce a broad range of ejecta properties. The characteristics of this ejecta, in turn, define the kilonova emission. To explore the effect of ejecta properties, we present a grid of 2-component 2D axisymmetric kilonova simulations that vary mass, velocity, morphology, and composition. The masses and velocities of each component vary, respectively, from 0.001 to 0.1 M ⊙ and 0.05 to 0.3c, covering much of the range of results from the neutron star merger literature. The set of 900 models is constrained to have a toroidal low electron fraction (Y e ) ejecta with a robust r-process composition and either a spherical or lobed high-Y e ejecta with two possible compositions. We simulate these models with the Monte Carlo radiative transfer code SuperNu using a full suite of lanthanide and 4th row element opacities. We examine the trends of these models with parameter variation, show how it can be used with statistical tools, and compare the model light curves and spectra to those of AT2017gfo, the electromagnetic counterpart of GW170817.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2014

Reconstructing the Sky Location of Gravitational-Wave Detected Compact Binary Systems: Methodology for Testing and Comparison

The problem of reconstructing the sky position of compact binary coalescences detected via gravitational waves is a central one for future observations with the ground-based network of gravitational-wave laser interferometers, such as Advanced LIGO and Advanced Virgo. Different techniques for sky localization have been independently developed. They can be divided in two broad categories: fully coherent Bayesian techniques, which are high latency and aimed at in-depth studies of all the parameters of a source, including sky position, and "triangulation-based" techniques, which exploit the data products from the search stage of the analysis to provide an almost real-time approximation of the posterior probability density function of the sky location of a detection candidate. These techniques have previously been applied to data collected during the last science runs of gravitational-wave detectors operating in the so-called initial configuration. Here, we develop and analyze methods for assessing the self consistency of parameter estimation methods and carrying out fair comparisons between different algorithms, addressing issues of efficiency and optimality. These methods are general, and can be applied to parameter estimation problems other than sky localization. We apply these methods to two existing sky localization techniques representing the two above-mentioned categories, using a set of simulated inspiralonly signals from compact binary systems with a total mass of equal to or less than 20M solar mass and nonspinning components. We compare the relative advantages and costs of the two techniques and show that sky location uncertainties are on average a factor approx. equals 20 smaller for fully coherent techniques than for the specific variant of the triangulation-based technique used during the last science runs, at the expense of a factor approx. equals 1000 longer processing time.

technique used↗