Self-Consistent Picture of the Mass Ejection from a One Second Long Binary Neutron Star Merger Leaving a Short-Lived Remnant in a General-Relativistic Neutrino-Radiation Magnetohydrodynamic Simulation

We perform a general-relativistic neutrino-radiation magnetohydrodynamic simulation of a one second-long binary neutron star merger on the Japanese supercomputer Fugaku using about 85 million CPU hours with 20 736 CPUs. We consider an asymmetric binary neutron star merger with masses of 1.2M_{⊙} and 1.5M_{⊙} and a "soft" equation of state SFHo. It results in a short-lived remnant with the lifetime of ≈0.017  s, and subsequent massive torus formation with the mass of ≈0.05M_{⊙} after the remnant collapses to a black hole. For the first time, we find that after the dynamical mass ejection, which drives the fast tail and mildly relativistic components, the postmerger mass ejection from the massive torus takes place due to the magnetorotational instability-driven turbulent viscosity in a single simulation and the two ejecta components are seen in the distributions of the electron fraction and velocity with distinct features.

Toward Calibration of the Global Network of Gravitational Wave Detectors with Sub-Percent Absolute and Relative Accuracy

The detection of gravitational-wave signals by the LIGO and Virgo observatories during the past few years has ushered us into the era of gravitational-wave astronomy, shifting our focus from detection to source parameter estimation. This has imposed stringent requirements on calibration in order to maximize the astrophysical information extracted from these detected signals. Current detectors rely on photon radiation pressure from auxiliary lasers to achieve required calibration accuracy. These photon calibrators have made significant improvements over the last few years, realizing fiducials displacements with sub-percent accuracy. This achieved accuracy is directly dependent on the laser power calibration. For the next observing campaign, scheduled to begin at the end of 2022, a new scheme is being implemented to achieve improved laser power calibration accuracy for all of the GW detectors in the global network. It is expected to significantly improve absolute and relative calibration accuracy for the entire network.

Inference-Optimized AI and High Performance Computing for Gravitational Wave Detection at Scale

We introduce an ensemble of artificial intelligence models for gravitational wave detection that we trained in the Summit supercomputer using 32 nodes, equivalent to 192 NVIDIA V100 GPUs, within 2 h. Once fully trained, we optimized these models for accelerated inference using NVIDIA TensorRT. We deployed our inference-optimized AI ensemble in the ThetaGPU supercomputer at Argonne Leadership Computer Facility to conduct distributed inference. Using the entire ThetaGPU supercomputer, consisting of 20 nodes each of which has 8 NVIDIA A100 Tensor Core GPUs and 2 AMD Rome CPUs, our NVIDIA TensorRT-optimized AI ensemble processed an entire month of advanced LIGO data (including Hanford and Livingston data streams) within 50 s. Our inference-optimized AI ensemble retains the same sensitivity of traditional AI models, namely, it identifies all known binary black hole mergers previously identified in this advanced LIGO dataset and reports no misclassifications, while also providing a 3X inference speedup compared to traditional artificial intelligence models. We used time slides to quantify the performance of our AI ensemble to process up to 5 years worth of advanced LIGO data. In this synthetically enhanced dataset, our AI ensemble reports an average of one misclassification for every month of searched advanced LIGO data. We also present the receiver operating characteristic curve of our AI ensemble using this 5 year long advanced LIGO dataset. This approach provides the required tools to conduct accelerated, AI-driven gravitational wave detection at scale.

The missing link in gravitational-wave astronomy: A summary of discoveries waiting in the decihertz range

Since 2015 the gravitational-wave observations of LIGO and Virgo have transformed our understanding of compact-object binaries. In the years to come, ground-based gravitational-wave observatories such as LIGO, Virgo, and their successors will increase in sensitivity, discovering thousands of stellar-mass binaries. In the 2030s, the space-based LISA will provide gravitational-wave observations of massive black holes binaries. Between the ∼ 10 -103 Hz band of ground-based observatories and the ∼ 1 0 - 4 -10- 1 Hz band of LISA lies the uncharted decihertz gravitational-wave band. We propose a Decihertz Observatory to study this frequency range, and to complement observations made by other detectors. Decihertz observatories are well suited to observation of intermediate-mass ( ∼ 1 0 2 -104 M ⊙) black holes; they will be able to detect stellar-mass binaries days to years before they merge, providing early warning of nearby binary neutron star mergers and measurements of the eccentricity of binary black holes, and they will enable new tests of general relativity and the Standard Model of particle physics. Here we summarise how a Decihertz Observatory could provide unique insights into how black holes form and evolve across cosmic time, improve prospects for both multimessenger astronomy and multiband gravitational-wave astronomy, and enable new probes of gravity, particle physics and cosmology.

Multimessenger Binary Mergers Containing Neutron Stars: Gravitational Waves, Jets, and γ-Ray Bursts

Neutron stars (NSs) are extraordinary not only because they are the densest form of matter in the visible Universe but also because they can generate magnetic fields ten orders of magnitude larger than those currently constructed on earth. The combination of extreme gravity with the enormous electromagnetic (EM) fields gives rise to spectacular phenomena like those observed on August 2017 with the merger of a binary neutron star system, an event that generated a gravitational wave (GW) signal, a short γ -ray burst (sGRB), and a kilonova. This event serves as the highlight so far of the era of multimessenger astronomy. In this review, we present the current state of our theoretical understanding of compact binary mergers containing NSs as gleaned from the latest general relativistic magnetohydrodynamic simulations. Such mergers can lead to events like the one on August 2017, GW170817, and its EM counterparts, GRB 170817 and AT 2017gfo. In addition to exploring the GW emission from binary black hole-neutron star and neutron star-neutron star mergers, we also focus on their counterpart EM signals. In particular, we are interested in identifying the conditions under which a relativistic jet can be launched following these mergers. Such a jet is an essential feature of most sGRB models and provides the main conduit of energy from the central object to the outer radiation regions. Jet properties, including their lifetimes and Poynting luminosities, the effects of the initial magnetic field geometries and spins of the coalescing NSs, as well as their governing equation of state, are discussed. Lastly, we present our current understanding of how the Blandford-Znajek mechanism arises from merger remnants as the trigger for launching jets, if, when and how a horizon is necessary for this mechanism, and the possibility that it can turn on in magnetized neutron ergostars, which contain ergoregions, but no horizons.

Black hole-neutron star coalescence: Effects of the neutron star spin on jet launching and dynamical ejecta mass

Black hole-neutron star (BHNS) mergers are thought to be sources of gravitational waves (GWs) with coincident electromagnetic (EM) counterparts. To further probe whether these systems are viable progenitors of short gamma-ray bursts (SGRBs) and kilonovas, and how one may use (the lack of) EM counterparts associated with LIGO/Virgo candidate BHNS GW events to sharpen parameter estimation, we study the impact of neutron star spin in BHNS mergers. Using dynamical spacetime magnetohydrodynamic simulations of BHNSs initially on a quasicircular orbit, we survey configurations that differ in the BH spin (a BH/M BH = 0 and 0.75), the NS spin (a NS/M NS = -0.17, 0, 0.23, and 0.33), and the binary mass ratio (q = M BH:M NS = 3:1 and 5:1). The general trend we find is that increasing the NS prograde spin increases both the rest mass of the accretion disk onto the remnant black hole, and the rest mass of dynamically ejected matter. By a time Δt ~ 3500-5500M ~ 88-138(M NS/1.4 M ⊙) ms after the peak gravitational-wave amplitude, a magnetically driven jet is launched only for q = 3:1 regardless of the initial NS spin. The lifetime of the jets [Δt ~ 0.5-0.8(M NS/1.4 M ⊙) s] and their outgoing Poynting luminosity [L Poyn ~ 1051.5±0.5 erg/s] are consistent with typical SGRBs' luminosities and expectations from the Blandford-Znajek mechanism. By the time we terminate our simulations, we do not observe either an outflow or a large-scale magnetic-field collimation for the other systems we consider. The mass range of dynamically ejected matter is 10-4.5-10-2(M NS/1.4 M ⊙) M ⊙, which can power kilonovas with peak bolometric luminosities L knova ~ 1040-1041.4 erg/s with rise times ≲6.5 h and potentially detectable by the LSST.

Measuring the Hubble constant with a sample of kilonovae

Kilonovae produced by the coalescence of compact binaries with at least one neutron star are promising standard sirens for an independent measurement of the Hubble constant (H₀). Through their detection via follow-up of gravitational-wave (GW), short gamma-ray bursts (sGRBs) or optical surveys, a large sample of kilonovae (even without GW data) can be used for H₀ contraints. Here, we show measurement of H₀ using light curves associated with four sGRBs, assuming these are attributable to kilonovae, combined with GW170817. Including a systematic uncertainty on the models that is as large as the statistical ones, we find \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${H}_{0}=73.{8}_{-5.8}^{+6.3}\ {\rm{km}}\ {{\rm{s}}}^{-1}\ {{\rm{Mpc}}}^{-1}$$\end{document}H0=73.8−5.8+6.3kms−1Mpc−1 and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${H}_{0}=71.{2}_{-3.1}^{+3.2}\ {\rm{km}}\ {{\rm{s}}}^{-1}\ {{\rm{Mpc}}}^{-1}$$\end{document}H0=71.2−3.1+3.2kms−1Mpc−1 for two different kilonova models that are consistent with the local and inverse-distance ladder measurements. For a given model, this measurement is about a factor of 2-3 more precise than the standard-siren measurement for GW170817 using only GWs.

Kilonovae observations can be used to out constraints on the Hubble constant (H0). Here, the authors show H0 measurements by combining light curves of four short gamma-ray burts with GW170817 are about a factor of 2-3 more precise than the standard-siren measurements using only gravitational-waves.

Asymmetric mass ratios for bright double neutron-star mergers

The discovery of a radioactively powered kilonova associated with the binary neutron-star merger GW170817 remains the only confirmed electromagnetic counterpart to a gravitational-wave event^1,2. Observations of the late-time electromagnetic emission, however, do not agree with the expectations from standard neutron-star merger models. Although the large measured ejecta mass^3,4 could be explained by a progenitor system that is asymmetric in terms of the stellar component masses (that is, with a mass ratio q of 0.7 to 0.8)⁵, the known Galactic population of merging double neutron-star systems (that is, those that will coalesce within billions of years or less) has until now consisted only of nearly equal-mass (q > 0.9) binaries⁶. The pulsar PSR J1913+1102 is a double system in a five-hour, low-eccentricity (0.09) orbit, with an orbital separation of 1.8 solar radii⁷, and the two neutron stars are predicted to coalesce in [Formula: see text] million years owing to gravitational-wave emission. Here we report that the masses of the pulsar and the companion neutron star, as measured by a dedicated pulsar timing campaign, are 1.62 ± 0.03 and 1.27 ± 0.03 solar masses, respectively. With a measured mass ratio of q = 0.78 ± 0.03, this is the most asymmetric merging system reported so far. On the basis of this detection, our population synthesis analysis implies that such asymmetric binaries represent between 2 and 30 per cent (90 per cent confidence) of the total population of merging binaries. The coalescence of a member of this population offers a possible explanation for the anomalous properties of GW170817, including the observed kilonova emission from that event.

Continuous Gravitational Waves from Neutron Stars: Current Status and Prospects

Gravitational waves astronomy allows us to study objects and events invisible in electromagnetic waves. It is crucial to validate the theories and models of the most mysterious and extreme matter in the Universe: the neutron stars. In addition to inspirals and mergers of neutrons stars, there are currently a few proposed mechanisms that can trigger radiation of long-lasting gravitational radiation from neutron stars, such as e.g., elastically and/or magnetically driven deformations: mountains on the stellar surface supported by the elastic strain or magnetic field, free precession, or unstable oscillation modes (e.g., the r-modes). The astrophysical motivation for continuous gravitational waves searches, current LIGO and Virgo strategies of data analysis and prospects are reviewed in this work.

Multimessenger Bayesian parameter inference of a binary neutron star merger

The combined detection of a binary neutron star merger in both gravitational waves (GWs) and electromagnetic (EM) radiation spanning the entire spectrum – GW170817/AT2017gfo/GRB170817A – marks a breakthrough in the field of multimessenger astronomy. Between the plethora of modelling and observations, the rich synergy that exists among the available data sets creates a unique opportunity to constrain the binary parameters, the equation of state of supranuclear density matter, and the physical processes at work during the kilonova and gamma-ray burst. We report, for the first time, Bayesian parameter estimation combining information from GW170817, AT2017gfo, and GRB170817 to obtain truly multimessenger constraints on the tidal deformability $\tilde{\Lambda } \in [302,860]$, total binary mass M ∈ [2.722, 2.751] M⊙, the radius of a 1.4 solar mass neutron star $R \in [11.3,13.5] \,\,\rm km$ (with additional $0.2\ \rm km$ systematic uncertainty), and an upper bound on the mass ratio of q ≤ 1.27, all at 90 per cent confidence. Our joint novel analysis uses new phenomenological descriptions of the dynamical ejecta, debris disc mass, and remnant black hole properties, all derived from a large suite of numerical relativity simulations.

Constraining the p-Mode-g-Mode Tidal Instability with GW170817

We analyze the impact of a proposed tidal instability coupling p modes and g modes within neutron stars on GW170817. This nonresonant instability transfers energy from the orbit of the binary to internal modes of the stars, accelerating the gravitational-wave driven inspiral. We model the impact of this instability on the phasing of the gravitational wave signal using three parameters per star: an overall amplitude, a saturation frequency, and a spectral index. Incorporating these additional parameters, we compute the Bayes factor (lnB_{!pg}^{pg}) comparing our p-g model to a standard one. We find that the observed signal is consistent with waveform models that neglect p-g effects, with lnB_{!pg}^{pg}=0.03_{-0.58}^{+0.70} (maximum a posteriori and 90% credible region). By injecting simulated signals that do not include p-g effects and recovering them with the p-g model, we show that there is a ≃50% probability of obtaining similar lnB_{!pg}^{pg} even when p-g effects are absent. We find that the p-g amplitude for 1.4  M_{⊙} neutron stars is constrained to less than a few tenths of the theoretical maximum, with maxima a posteriori near one-tenth this maximum and p-g saturation frequency ∼70  Hz. This suggests that there are less than a few hundred excited modes, assuming they all saturate by wave breaking. For comparison, theoretical upper bounds suggest ≲10^{3} modes saturate by wave breaking. Thus, the measured constraints only rule out extreme values of the p-g parameters. They also imply that the instability dissipates ≲10^{51}  erg over the entire inspiral, i.e., less than a few percent of the energy radiated as gravitational waves.

Search for Neutrinos in Super-Kamiokande Associated with Gravitational Wave Events

We report the results from a search in Super-Kamiokande for neutrino signals coincident with gravitational-wave events using a neutrino energy range from 3.5 MeV–100 PeV. We searched for coincident neutrino events within a time window of ± 500 s around the gravitational-wave detection time. In this work, we report the number of events within the search-window and the 90 % confidence level upper limits on the neutrino fluence for each gravitational-wave event.

Convective Excitation of Inertial Modes in Binary Neutron Star Mergers

We present the first very long-term simulations (extending up to ∼140  ms after merger) of binary neutron star mergers with piecewise polytropic equations of state and in full general relativity. Our simulations reveal that, at a time of 30-50 ms after merger, parts of the star become convectively unstable, which triggers the excitation of inertial modes. The excited inertial modes are sustained up to several tens of milliseconds and are potentially observable by the planned third-generation gravitational-wave detectors at frequencies of a few kilohertz. Since inertial modes depend on the rotation rate of the star and they are triggered by a convective instability in the postmerger remnant, their detection in gravitational waves will provide a unique opportunity to probe the rotational and thermal state of the merger remnant. In addition, our findings have implications for the long-term evolution and stability of binary neutron star remnants.