Compact Binary Mergers Research Articles

On non-detection of gamma-ray bursts in three compact binary merger events detected by LIGO

ABSTRACT The joint detection of the gravitational wave (GW) event GW170817 and the short-duration gamma-ray burst (SGRB) event GRB 170817A, marked the beginning of GW multimessenger astronomy and confirmed that binary neutron star mergers are progenitors of at least some SGRBs. An estimated joint detection rate of 0.3–1.7 per year between the LIGO-Hanford, LIGO-Livingston, and Virgo GW network at design sensitivity, and the Fermi Gamma-ray Burst Monitor was predicted. However, to date, the GW170817/GRB 170817A joint detection has been the only event of its kind so far. Taking into account that SGRBs are narrowly beamed and are emitted perpendicular to the orbital plane of the binary system, we propose that previous mergers involving neutron stars, were orientated such that observation of the emitted SGRB along this narrow jet was not possible. To support this hypothesis we have estimated the inclination of the binary systems for previously detected Binary Neutron Star (BNS) and Black Hole Neutron Star (BHNS) mergers through GW analysis. This analysis was performed using bilby, a python based Bayesian inference library, to estimate the inclination of the BNS events GW170817 and GW190425, and the BHNS events GW190917_114630 and GW200115_042309. The results obtained in this study indicate that these binaries may have had inclinations greater than 33° with respect to the line of sight from Earth, an upper limit on the viewing angle set from observations of GRB 170817A. This then suggests that the observation of the emitted SGRB from these past mergers might not have been possible.

Impact of Ejecta Temperature and Mass on the Strength of Heavy Element Signatures in Kilonovae

A kilonova, the electromagnetic emission produced by compact binary mergers, is formed through a delicate interplay of physical processes, involving r-process nucleosynthesis and interactions between heavy elements and photons through radiative transfer. This complexity makes it difficult to achieve a comprehensive understanding of kilonova spectra. In this study, we aim to enhance our understanding and establish connections between physical parameters and observables through radiative-transfer simulations. Specifically, we investigate how ejecta temperature and element mass influence the resulting kilonova spectrum. For each species, the strength of its line features depends on these parameters, leading to the formation of a distinct region in the parameter space, dubbed the resonance island, where the line signature of that species is notably evident in the kilonova spectrum. We explore its origin and applications. Among explored r-process elements (31 ≤ Z ≤ 92), we find that four species—SrII, YII, BaII, and CeII—exhibit large and strong resonance islands, suggesting their significant contributions to kilonova spectra at specific wavelengths. In addition, we discuss potential challenges and future perspectives in observable heavy elements and their masses in the context of the resonance island.

Amplification and Dissipation of Magnetic Fields in Accreting Compact Objects

Magnetic fields play a crucial role in shaping the dynamics of accreting compact objects. Whether we consider the formation of a proto-neutron star during the gravitational collapse of a massive star or the accretion disk around a black hole after a compact binary merger, a key process that remains challenging to include in large-scale simulations is the amplification and dissipation of magnetic fields driven by turbulent fluid motions. Despite the enormous increase in computational power currently available, the large separation between all the relevant spatial and temporal scales still poses severe limits to what can be achieved with ideal fluid simulations. One way to tackle such issue is to rely on sub-grid models, which however need to be appropriately tuned in light of models probing the small-scale dynamics. In this work we present the current state-of-the-art of dynamo models in proto-neutron stars, which aim at describing the amplification of magnetar-like magnetic fields during the gravitational collapse of a massive star. We also review some of the works from the past few years that included turbulent dynamos in accretion disks around a black holes, relying on a mean-field formalism. Finally, we will present a recent study on polar jets with explicit turbulent resistivity which showcases the importance of employing highly accurate numerical schemes.

Searches for neutrino counterparts of gravitational waves from the LIGO/Virgo third observing run with KM3NeT

The KM3NeT neutrino telescope is currently being deployed at two different sites in the Mediterranean Sea. First searches for astrophysical neutrinos have been performed using data taken with the partial detector configuration already in operation. The paper presents the results of two independent searches for neutrinos from compact binary mergers detected during the third observing run of the LIGO and Virgo gravitational wave interferometers. The first search looks for a global increase in the detector counting rates that could be associated with inverse beta decay events generated by MeV-scale electron anti-neutrinos. The second one focuses on upgoing track-like events mainly induced by muon (anti-)neutrinos in the GeV–TeV energy range. Both searches yield no significant excess for the sources in the gravitational wave catalogs. For each source, upper limits on the neutrino flux and on the total energy emitted in neutrinos in the respective energy ranges have been set. Stacking analyses of binary black hole mergers and neutron star-black hole mergers have also been performed to constrain the characteristic neutrino emission from these categories.

Cosmology and Astrophysics with Standard Sirens and Galaxy Catalogs in View of Future Gravitational Wave Observations

With the growing number of gravitational wave (GW) detections and the advent of large galaxy redshift surveys, a new era in cosmology is emerging. This study explores the synergies between GWs and galaxy surveys to jointly constrain cosmological and GW population parameters. We introduce CHIMERA, a novel code for GW cosmology combining information from the population properties of compact binary mergers and galaxy catalogs. We study constraints for scenarios representative of the LIGO-Virgo-KAGRA O4 and O5 observing runs, assuming to have a complete catalog of potential host galaxies with either spectroscopic or photometric redshift measurements. We find that a percent-level measurement of H 0 could be achieved with the best 100 binary black holes (BBHs) in O5 using a spectroscopic galaxy catalog. In this case, the intrinsic correlation that exists between H 0 and the BBH population mass scales is broken. Instead, by using a photometric catalog the accuracy is degraded up to a factor of ∼9, leaving a significant correlation between H 0 and the mass scales that must be carefully modeled to avoid bias. Interestingly, we find that using spectroscopic redshift measurements in the O4 configuration yields a better constraint on H 0 compared to the O5 configuration with photometric measurements. In view of the wealth of GW data that will be available in the future, we argue the importance of obtaining spectroscopic galaxy catalogs to maximize the scientific return of GW cosmology.

3D NLTE modelling of Y and Eu

Context. Abundances of s- and r-process elements in Sun-like stars constrain nucleosynthesis in extreme astrophysical events, such as compact binary mergers and explosions of highly magnetised rapidly rotating massive stars. Aims. We measure solar abundances of yttrium (Y) and europium (Eu) using 3D non-local thermal equilibrium (NLTE) models. We use the model to determine the abundance of Y, and also explore the model’s ability to reproduce the solar centre-to-limb variation of its lines. In addition, we determine the Eu abundance using solar disc-centre and integrated flux spectra. Methods. We developed an NLTE model of Eu and updated our model of Y with collisional data from detailed quantum-mechanical calculations. We used the IAG spatially resolved high-resolution solar spectra to derive the solar abundances of Y across the solar disc and of Eu for integrated flux and at disc centre using a set of carefully selected lines and a 3D radiation-hydrodynamics model of the solar atmosphere. Results. We find 3D NLTE solar abundances of A(Y)3D NLTE = 2.30 ± 0.03stat ± 0.07syst dex based on observations at all angles and A(Eu) = 0.57 ± 0.01stat ± 0.06syst dex based on the integrated flux and disc-centre intensity. 3D NLTE modelling offers the most consistent abundances across the solar disc, and resolves the problem of severe systematic bias in Y and Eu abundances inherent to 1D LTE, 1D NLTE, and 3D LTE modelling.

Template bank for compact binary mergers in the fourth observing run of Advanced LIGO, Advanced Virgo, and KAGRA

Template bank for compact binary mergers in the fourth observing run of Advanced LIGO, Advanced Virgo, and KAGRA

Xkn: a semi-analytic framework for the modelling of kilonovae

Abstract After GW170817, kilonovae have become of great interest for the astronomical, astrophysics and nuclear physics communities, due to their potential in revealing key information on the compact binary merger from which they emerge, such as the fate of the central remnant or the composition of the expelled material. Therefore, the landscape of models employed for their analysis is rapidly evolving, with multiple approaches being used for different purposes. In this paper, we present xkn, a semi-analytic framework which predicts and interprets the bolometric luminosity and the broadband light curves of such transients. xkn models the merger ejecta structure accounting for different ejecta components and non-spherical geometries. In addition to light curve models from the literature based on time scale and random-walk arguments, it implements a new model, xkn-diff, which is grounded on a solution of the radiative transfer equation for homologously expanding material. In order to characterize the variety of the ejecta conditions, it employs time and composition dependent heating rates, thermalization efficiencies and opacities. We compare xkn light curves with reference radiative transfer calculations, and we find that xkn-diff significantly improves over previous semi-analytic prescriptions. We view xkn as an ideal tool for extensive parameter estimation data analysis applications.

Perturbing fast neutrino flavor conversion

The flavor evolution of neutrinos in dense astrophysical sources, such as core-collapse supernovae or compact binary mergers, is non-linear due to the coherent forward scattering of neutrinos among themselves. Recent work in this context has been addressed to figure out whether flavor equipartition could be a generic flavor outcome of fast flavor conversion. We investigate the flavor conversion physics injecting random perturbations in the neutrino field in two simulation setups: 1. a spherically symmetric simulation shell without periodic boundaries, with angular distributions evolving dynamically thanks to non-forward scatterings of neutrinos with the background medium, and neutrino advection; 2. a periodic simulation shell, with angular distributions of neutrinos defined a priori and neutrino advection. We find that, independent of the exact initial flavor configuration and type of perturbations, flavor equipartition is generally achieved in the system with periodic boundaries; in this case, perturbations aid the diffusion of flavor structures to smaller and smaller scales. However, flavor equipartition is not a general outcome in the simulation shell without periodic boundaries, where the inhomogeneities induced perturbing the neutrino field affect the flavor evolution, but do not facilitate the diffusion of flavor waves. This work highlights the importance of the choice of the simulation boundary conditions in the exploration of fast flavor conversion physics.

New and robust gravitational-waveform model for high-mass-ratio binary neutron star systems with dynamical tidal effects

For the analysis of gravitational-wave signals, fast and accurate gravitational-waveform models are required. These enable us to obtain information on the system properties from compact binary mergers. In this article, we introduce the model, which contains a closed-form expression that describes tidal effects, focusing on the description of binary neutron star systems. The model improves upon previous versions by employing a larger set of numerical-relativity data for its calibration, by including high-mass ratio systems covering also a wider range of equations of state. It also takes into account dynamical tidal effects and the known post-Newtonian mass-ratio dependence of individual calibration parameters. We implemented the model in the publicly available uite software library by augmenting different binary black hole waveform models (, , and ). We test the validity of by comparing it with numerical-relativity waveforms, as well as other tidal models. Finally, we perform parameter estimation for GW170817 and GW190425 with the new tidal approximant and find overall consistent results with respect to previous studies. Published by the American Physical Society 2024

Targeted search for gravitational waves from highly spinning light compact binaries

Abstract Searches for gravitational waves from compact-binary mergers, which to date have reported ∼100 observations, have previously ignored binaries whose components are both consistent with the mass of neutron stars (1M⊙ to 2M⊙) and have high dimensionless spin >0.05. While previous searches targeted sources that are representative of observed neutron star binaries in the galaxy, it is already known that neutron stars can regularly be spun up to a dimensionless spin of ∼0.4, and in principle reach up to ∼0.7 before breakup would occur. Furthermore, there may be primordial black hole binaries or exotic formation mechanisms to produce light black holes. In these cases, it is possible for the binary constituent to be spun up beyond that achievable by a neutron star. A single detection of this type of source would reveal a novel formation channel for compact-binaries. To determine if there is evidence for any such sources, we use PyCBC to conduct a targeted search of LIGO and Virgo data for light compact objects with high spin. Our analysis detects previously known observations GW170817 and GW200115; however, we report no additional mergers. The most significant candidate, not previously known, is consistent with the noise distribution, and so we constrain the merger rate of spinning light binaries.

The dragon-II simulations – III. Compact binary mergers in clusters with up to 1 million stars: mass, spin, eccentricity, merger rate, and pair instability supernovae rate

ABSTRACT Compact binary mergers forming in star clusters may exhibit distinctive features that can be used to identify them among observed gravitational-wave sources. Such features likely depend on the host cluster structure and the physics of massive star evolution. Here, we dissect the population of compact binary mergers in the dragon-II simulation data base, a suite of 19 direct N-body models representing dense star clusters with up to 106 stars and $\lt 33~{{\ \rm per\ cent}}$ of stars in primordial binaries. We find a substantial population of black hole binary (BBH) mergers, some of them involving an intermediate-mass BH (IMBH), and a handful mergers involving a stellar BH and either a neutron star (NS) or a white dwarf (WD). Primordial binary mergers, $\sim 30~{{\ \rm per\ cent}}$ of the whole population, dominate ejected mergers. Dynamical mergers, instead, dominate the population of in-cluster mergers and are systematically heavier than primordial ones. Around 20 per cent of dragon-II mergers are eccentric in the Laser Interferometer Space Antenna (LISA) band and 5 per cent in the LIGO band. We infer a mean cosmic merger rate of $\mathcal {R}\sim 30(4.4)(1.2)$ yr−1 Gpc−3 for BBHs, NS–BH, and WD–BH binary mergers, respectively, and discuss the prospects for multimessenger detection of WD–BH binaries with LISA. We model the rate of pair-instability supernovae (PISNe) in star clusters and find that surveys with a limiting magnitude mbol = 25 can detect ∼1–15 yr−1 PISNe. Comparing these estimates with future observations could help to pin down the impact of massive star evolution on the mass spectrum of compact stellar objects in star clusters.

Search technique to observe precessing compact binary mergers in the advanced detector era

Search technique to observe precessing compact binary mergers in the advanced detector era

Identifying the Physical Origin of Gamma-Ray Bursts with Supervised Machine Learning

The empirical classification of gamma-ray bursts (GRBs) into long and short GRBs based on their durations is already firmly established. This empirical classification is generally linked to the physical classification of GRBs originating from compact binary mergers and GRBs originating from massive star collapses, or Type I and II GRBs, with the majority of short GRBs belonging to Type I and the majority of long GRBs belonging to Type II. However, there is a significant overlap in the duration distributions of long and short GRBs. Furthermore, some intermingled GRBs, i.e., short-duration Type II and long-duration Type I GRBs, have been reported. A multiparameter classification scheme of GRBs is evidently needed. In this paper, we seek to build such a classification scheme with supervised machine-learning methods, chiefly XGBoost. We utilize the GRB Big Table and Greiner’s GRB catalog and divide the input features into three subgroups: prompt emission, afterglow, and host galaxy. We find that the prompt emission subgroup performs the best in distinguishing between Type I and II GRBs. We also find the most important distinguishing features in prompt emission to be T 90, the hardness ratio, and fluence. After building the machine-learning model, we apply it to the currently unclassified GRBs to predict their probabilities of being either GRB class, and we assign the most probable class of each GRB to be its possible physical class.

Probing Cosmology with Baryon Acoustic Oscillations Using Gravitational Waves

The third-generation (3G) gravitational wave (GW) detectors such as the Einstein Telescope or Cosmic Explorer are expected to play an important role in cosmology. With the help of 3G detectors, we will be able to probe large-scale structure features such as baryon acoustic oscillations (BAO), galaxy bias, etc. We explore the possibility to do precision cosmology, with the 3G GW detectors by measuring the angular BAO scale using localization volumes of compact binary merger events. Through simulations, we show that with a 3G detector network, by probing the angular BAO scale using purely GW observations, we can constrain the Hubble constant for the standard model of cosmology (ΛCDM) with 90% credible regions as . When combined with BAO measurements from galaxy surveys, we show that it can be used to constrain various models of cosmology such as parameterized models for dark energy equations of state. We also show how cosmological constraints using BAO measurements from GW observations in the 3G era will complement the same from spectroscopic surveys.

A Unified Picture of Short and Long Gamma-Ray Bursts from Compact Binary Mergers

The recent detections of the ∼10 s long γ-ray bursts (GRBs) 211211A and 230307A followed by softer temporally extended emission (EE) and kilonovae point to a new GRB class. Using state-of-the-art first-principles simulations, we introduce a unifying theoretical framework that connects binary neutron star (BNS) and black hole–NS (BH–NS) merger populations with the fundamental physics governing compact binary GRBs (cbGRBs). For binaries with large total masses, M tot ≳ 2.8 M ⊙, the compact remnant created by the merger promptly collapses into a BH surrounded by an accretion disk. The duration of the pre-magnetically arrested disk (MAD) phase sets the duration of the roughly constant power cbGRB and could be influenced by the disk mass, M d . We show that massive disks (M d ≳ 0.1 M ⊙), which form for large binary mass ratios q ≳ 1.2 in BNS or q ≲ 3 in BH–NS mergers, inevitably produce 211211A-like long cbGRBs. Once the disk becomes MAD, the jet power drops with the mass accretion rate as , establishing the EE decay. Two scenarios are plausible for short cbGRBs. They can be powered by BHs with less massive disks, which form for other q values. Alternatively, for binaries with M tot ≲ 2.8 M ⊙, mergers should go through a hypermassive NS (HMNS) phase, as inferred for GW170817. Magnetized outflows from such HMNSs, which typically live for ≲1 s, offer an alternative progenitor for short cbGRBs. The first scenario is challenged by the bimodal GRB duration distribution and the fact that the Galactic BNS population peaks at sufficiently low masses that most mergers should go through an HMNS phase.

General relativistic simulations of high-mass binary neutron star mergers: rapid formation of low-mass stellar black holes

ABSTRACT Almost a hundred compact binary mergers have been detected via gravitational waves by the LIGO-Virgo-KAGRA collaboration in the past few years providing us with a significant amount of new information on black holes and neutron stars. In addition to observations, numerical simulations using newly developed modern codes in the field of gravitational wave physics will guide us to understand the nature of single and binary degenerate systems and highly energetic astrophysical processes. We here presented a set of new fully general relativistic hydrodynamic simulations of high-mass binary neutron star systems using the open-source EINSTEIN TOOLKIT and LORENE codes. We considered systems with total baryonic masses ranging from 2.8$\, \mathrm{M}_\odot$ to 4.0$\, \mathrm{M}_\odot$ and used the SLy equation of state. We analysed the gravitational wave signal for all models and report potential indicators of systems undergoing rapid collapse into a black hole that could be observed by future detectors like the Einstein Telescope and the Cosmic Explorer. The properties of the post-merger black hole, the disc and ejecta masses, and their dependence on the binary parameters were also extracted. We also compared our numerical results with recent analytical fits presented in the literature and provided parameter-dependent semi-analytical relations between the total mass and mass ratio of the systems and the resulting black hole masses and spins, merger frequency, BH formation time, ejected mass, disc mass, and radiated gravitational wave energy.

Gamma-ray burst precursors from tidally resonant neutron star oceans: potential implications for GRB 211211A

ABSTRACT Precursors have been observed seconds to minutes before some short gamma-ray bursts. While the precursor origins remain unknown, one explanation relies on the resonance of neutron star pulsational modes with the tidal forces during the inspiral phase of a compact binary merger. In this paper, we present a model for short gamma-ray burst precursors that relies on tidally resonant neutron star oceans. In this scenario, the onset of tidal resonance in the crust–ocean interface mode ignites the precursor flare, possibly through the interaction between the excited neutron star ocean and the surface magnetic fields. From just the precursor total energy, the time before the main event, and a detected quasi-periodic oscillation frequency, we may constrain the binary parameters and neutron star ocean properties. Our model can immediately distinguish neutron star–black hole mergers from binary neutron star mergers without gravitational wave detection. We apply our model to GRB 211211A, the recently detected long duration short gamma-ray burst with a quasi-periodic precursor, and explore the parameters of this system. The precursor of GRB 211211A is consistent with a tidally resonant neutron star ocean explanation that requires an extreme mass ratio neutron star–black hole merger and a high-mass neutron star. While difficult to reconcile with the main gamma-ray burst and associated kilonova, our results constrain the possible precursor mechanisms in this system. A systematic study of short gamma-ray burst precursors with the model presented here can test precursor origin and probe the possible connection between gamma-ray bursts and neutron star–black hole mergers.

Towards a robust and reliable deep learning approach for detection of compact binary mergers in gravitational wave data

The ability of deep learning (DL) approaches to learn generalised signal and noise models, coupled with their fast inference on GPUs, holds great promise for enhancing gravitational-wave (GW) searches in terms of speed, parameter space coverage, and search sensitivity. However, the opaque nature of DL models severely harms their reliability. In this work, we meticulously develop a DL model stage-wise and work towards improving its robustness and reliability. First, we address the problems in maintaining the purity of training data by deriving a new metric that better reflects the visual strength of the ‘chirp’ signal features in the data. Using a reduced, smooth representation obtained through a variational auto-encoder (VAE), we build a classifier to search for compact binary coalescence (CBC) signals. Our tests on real LIGO data show an impressive performance of the model. However, upon probing the robustness of the model through adversarial attacks, its simple failure modes were identified, underlining how such models can still be highly fragile. As a first step towards bringing robustness, we retrain the model in a novel framework involving a generative adversarial network (GAN). Over the course of training, the model learns to eliminate the primary modes of failure identified by the adversaries. Although absolute robustness is practically impossible to achieve, we demonstrate some fundamental improvements earned through such training, like sparseness and reduced degeneracy in the extracted features at different layers inside the model. We show that these gains are achieved at practically zero loss in terms of model performance on real LIGO data before and after GAN training. Through a direct search on days of LIGO data, we recover two significant CBC events from GWTC-2.1 (Abbott et al 2022 2108.01045 [gr-qc]), GW190519_153544 and GW190521_074359. We also report the search sensitivity obtained from an injection study.

Associating fast radio bursts with compact binary mergers via gravitational lensing

ABSTRACT The origin of fast radio bursts (FRBs) is currently an open question with several proposed sources and corresponding mechanisms for their production. Among them are compact binary coalescences (CBCs) that also generate gravitational waves (GWs). Spatial and temporal coincidences between GWs and FRBs have so far been used to search for potential FRB counterparts to GWs from CBCs. However, such methods suffer from relatively poor sky-localization of the GW sources, and similarly poor luminosity distance estimates of both GW and FRB sources. The expected time-delay between the GW and radio emission is also poorly understood. In this work, we propose an astrophysical scenario that could potentially provide an unambiguous association between CBCs and FRBs, if one exists, or unambiguously rule out FRB counterparts to a given CBC GW event. We demonstrate that if a CBC that emitted both GWs and FRBs is gravitationally lensed, we can make a >5σ association using time-delay estimates of the lensed GW and FRB images (in strong lensing), which are expected to be measured with millisecond (for GW) and nano-second (FRB) precisions. We also demonstrate that the CBC–FRB association can be made in the microlensing regime as well where wave-optics effects modulate the GW waveform. We further investigate the rate of such detected associations in future observing scenarios of both GW and radio detectors.