Abstract
General relativistic radiation hydrodynamic simulations are necessary to accurately model a number of astrophysical systems involving black holes and neutron stars. Photon transport plays a ...crucial role in radiatively dominated accretion discs, while neutrino transport is critical to core-collapse supernovae and to the modelling of electromagnetic transients and nucleosynthesis in neutron star mergers. However, evolving the full Boltzmann equations of radiative transport is extremely expensive. Here, we describe the implementation in the general relativistic spec code of a cheaper radiation hydrodynamic method that theoretically converges to a solution of Boltzmann's equation in the limit of infinite numerical resources. The algorithm is based on a grey two-moment scheme, in which we evolve the energy density and momentum density of the radiation. Two-moment schemes require a closure that fills in missing information about the energy spectrum and higher order moments of the radiation. Instead of the approximate analytical closure currently used in core-collapse and merger simulations, we complement the two-moment scheme with a low-accuracy Monte Carlo evolution. The Monte Carlo results can provide any or all of the missing information in the evolution of the moments, as desired by the user. As a first test of our methods, we study a set of idealized problems demonstrating that our algorithm performs significantly better than existing analytical closures. We also discuss the current limitations of our method, in particular open questions regarding the stability of the fully coupled scheme.
Of the three main types of binaries detectable through ground-based gravitational wave observations, black hole-neutron star (BHNS) mergers remain the most elusive. While candidates BHNS exist in the ...triggers released during the third observing run of the Advanced LIGO/Virgo collaboration, no detection has been confirmed so far. As for binary neutron star systems, BHNS binaries allow us to explore a wide range of physical processes, including the neutron star equation of state, nucleosynthesis, stellar evolution, high-energy astrophysics, and the expansion of the Universe. Here, we review some of the main features of BHNS systems: the distinction between disrupting and non-disrupting binaries, the types of outflows that BHNS mergers can produce, and the information that can be extracted from the observation of their gravitational wave and electromagnetic signals. We also emphasize that for the most likely binary parameters, BHNS mergers seem less likely to power electromagnetic signals than binary neutron star systems. Finally, we discuss some of the issues that still limit our ability to model and interpret electromagnetic signals from BHNS binaries.
Gravitational-wave (GW) and electromagnetic (EM) signals from the merger of a neutron star (NS) and a black hole (BH) are a highly anticipated discovery. We present a simple formula, validated with ...75 simulations, that distinguishes between potential merger outcomes and predicts the baryon mass left outside of the BH after merger. Our formula describes critical unexplored regimes: comparable masses with nonspinning BHs, and higher BH spins, and is essential in assessing whether events such as GW170817 could be NS-BH systems instead of NS-NS mergers.
There is irresistible observational evidence that binary systems of compact objects with at least one neutron star are progenitors of short gamma-ray bursts, as well as a production site for r ...-process elements, at least when some matter is ejected by the merger and an accretion disk is formed. The recent observations of gravitational waves in conjunction with electromagnetic counterparts fuel the need for models predicting the outcome of a given merger and the properties of the associated matter outflows as a function of the initial parameters of the binary. In this manuscript, we provide updated fitting formulas that estimate the disk mass for double neutron star binaries and ejecta masses for black hole–neutron star and double neutron star binaries, fitted to the results of numerical simulations. Our proposed fitting formulas improve on existing models by aiming for analytical simplicity, by covering a larger region of parameter space, and by accounting for regions of parameter space not covered by numerical simulations but with physically manifest merger outcomes.
ABSTRACT
GW170817 showed that neutron star mergers not only emit gravitational waves but also can release electromagnetic signatures in multiple wavelengths. Within the first half of the third ...observing run of the Advanced LIGO and Virgo detectors, there have been a number of gravitational wave candidates of compact binary systems for which at least one component is potentially a neutron star. In this article, we look at the candidates S190425z, S190426c, S190510g, S190901ap, and S190910h, predicted to have potentially a non-zero remnant mass, in more detail. All these triggers have been followed up with extensive campaigns by the astronomical community doing electromagnetic searches for their optical counterparts; however, according to the released classification, there is a high probability that some of these events might not be of extraterrestrial origin. Assuming that the triggers are caused by a compact binary coalescence and that the individual source locations have been covered during the EM follow-up campaigns, we employ three different kilonova models and apply them to derive possible constraints on the matter ejection consistent with the publicly available gravitational-wave trigger information and the lack of a kilonova detection. These upper bounds on the ejecta mass can be related to limits on the maximum mass of the binary neutron star candidate S190425z and to constraints on the mass-ratio, spin, and NS compactness for the potential black hole–neutron star candidate S190426c. Our results show that deeper electromagnetic observations for future gravitational wave events near the horizon limit of the advanced detectors are essential.
Searches for electromagnetic counterparts of gravitational-wave signals have redoubled since the first detection in2017 of a binary neutron star merger with a gamma-ray burst, optical/infrared ...kilonova, and panchromatic after glow. Yet, one LIGO/Virgo observing run later, there has not yet been a second, secure identification of an electromagnetic counterpart. This is not surprising given that the localization uncertainties of events in LIGO and Virgo’s third observing run, O3, were much larger than predicted. We explain this by showing that improvements in data analysis that now allow LIGO/Virgo to detect weaker and hence more poorly localized events have increased the overall number of detections, of which well-localized, gold-plated events make up a smaller proportion overall. We present simulations of the next two LIGO/Virgo/KAGRA observing runs, O4 and O5, that are grounded in the statistics ofO3 public alerts. To illustrate the significant impact that the updated predictions can have, we study the follow-up strategy for the Zwicky Transient Facility. Realistic and timely forecasting of gravitational-wave localization accuracy is paramount given the large commitments of telescope time and the need to prioritize which events are followed up. We include a data release of our simulated localizations as a public proposal planning resource for astronomers
Gravitational waves from neutron star binary inspirals contain information on strongly interacting matter in unexplored, extreme regimes. Extracting this requires robust theoretical models of the ...signatures of matter in the gravitational-wave signals due to spin and tidal effects. In fact, spins can have a significant impact on the tidal excitation of the quasinormal modes of a neutron star, which is not included in current state-of-the-art waveform models. We develop a simple approximate description that accounts for the Coriolis effect of spin on the tidal excitation of the neutron star's quadrupolar and octupolar fundamental quasinormal modes and incorporate it in the SEOBNRv4T waveform model. We show that the Coriolis effect introduces only one new interaction term in an effective action in the corotating frame of the star, and fix the coefficient by considering the spin-induced shift in the resonance frequencies that has been computed numerically for the mode frequencies of rotating neutron stars in the literature. We investigate the impact of relativistic corrections due to the gravitational redshift and frame-dragging effects, and identify important directions where more detailed theoretical developments are needed in the future. Comparisons of our model to numerical-relativity simulations of double neutron star and neutron star black-hole binaries show improved consistency in the agreement compared to current models used in data analysis.
Abstract
Black hole—neutron star (BH–NS) mergers are a major target for ground-based gravitational wave observatories. A merger can also produce an electromagnetic counterpart (a kilonova) if it ...ejects neutron-rich matter that assembles into heavy elements through
r
-process nucleosynthesis. We study the kilonova signatures of the unbound dynamical ejecta of a BH–NS merger. We take as our initial state the results from a numerical relativity simulation and then use a general relativistic hydrodynamics code to study the evolution of the ejecta with parameterized
r
-process heating models. The unbound dynamical ejecta is initially a flattened, directed tidal tail largely confined to a plane. Heating from the
r
-process inflates the ejecta into a more spherical shape and smooths its small-scale structure, though the ejecta retains its bulk directed motion. We calculate the electromagnetic signatures using a 3D radiative transfer code and a parameterized opacity model for lanthanide-rich matter. The light curve varies with viewing angle because of two effects: asphericity results in brighter emission for orientations with larger projected areas, while Doppler boosting results in brighter emission for viewing angles more aligned with the direction of bulk motion. For typical
r
-process heating rates, the peak bolometric luminosity varies by a factor of ∼3 with orientation while the peak in the optical bands varies by ∼3 magnitudes. The spectrum is blueshifted at viewing angles along the bulk motion, which increases the
V
-band peak magnitude to ∼−14 despite the lanthanide-rich composition.
Binary neutron star mergers are promising sources of gravitational waves for ground-based detectors such as Advanced LIGO. Neutron-rich material ejected by these mergers may also be the main source ...of r-process elements in the Universe, while radioactive decays in the ejecta can power bright electromagnetic postmerger signals. Neutrino-matter interactions play a critical role in the evolution of the composition of the ejected material, which significantly impacts the outcome of nucleosynthesis and the properties of the associated electromagnetic signal. In this work, we present a simulation of a binary neutron star merger using an improved method for estimating the average neutrino energies in our energy-integrated neutrino transport scheme. These energy estimates are obtained by evolving the neutrino number density in addition to the neutrino energy and flux densities. We show that significant changes are observed in the composition of the polar ejecta when comparing our new results with earlier simulations in which the neutrino spectrum was assumed to be the same everywhere in optically thin regions. In particular, we find that material ejected in the polar regions is less neutron rich than previously estimated. Our new estimates of the composition of the polar ejecta make it more likely that the color and time scale of the electromagnetic signal depend on the orientation of the binary with respect to an observer’s line of sight. These results also indicate that important observable properties of neutron star mergers are sensitive to the neutrino energy spectrum, and may need to be studied through simulations including a more accurate, energy-dependent neutrino transport scheme.
We investigate the ejecta from black hole-neutron star mergers by modeling the formation and interaction of mass ejected in a tidal tail and a disk wind. The outflows are neutron-rich, giving rise to ...optical/infrared emission powered by the radioactive decay of r-process elements (a kilonova). Here we perform an end-to-end study of this phenomenon, where we start from the output of a fully-relativistic merger simulation, calculate the post-merger hydrodynamical evolution of the ejecta and disk winds including neutrino physics, determine the final nucleosynthetic yields using post-processing nuclear reaction network calculations, and compute the kilonova emission with a radiative transfer code. We study the effects of the tail-to-disk mass ratio by scaling the tail density. A larger initial tail mass results in fallback matter becoming mixed into the disk and ejected in the subsequent disk wind. Relative to the case of a disk without dynamical ejecta, the combined outflow has lower mean electron fraction, faster speed, larger total mass, and larger absolute mass free of high-opacity Lanthanides or Actinides. In most cases, the nucleosynthetic yield is dominated by the heavy r-process contribution from the unbound part of the dynamical ejecta. A Solar-like abundance distribution can however be obtained when the total mass of the dynamical ejecta is comparable to the mass of the disk outflows. The kilonova has a characteristic duration of 1 week and a luminosity of 1041 erg s−1, with orientation effects leading to variations of a factor 2 in brightness. At early times (<1 d) the emission includes an optical component from the (hot) Lanthanide-rich material, but the spectrum evolves quickly to the infrared thereafter.