Binary Mergers Research Articles

Host Galaxy Properties of Gamma-Ray Bursts Involving Neutron Star Binary Mergers and Their Impact on Kilonovae Rates

In the upcoming gravitational-wave (GW) observing runs, identifying host galaxies is crucial as it provides essential redshift information and enables the use of GW events as standard sirens. However, pinpointing host galaxies remains challenging due to the large localization uncertainties and the rapidly fading nature of their optical counterparts. Analyzing the host galaxies of short gamma-ray bursts (sGRBs) offers an alternative approach to deepen our understanding of the environments where binary neutron stars (BNS) primarily merge. This study compiles archival photometric data for the host galaxies of 76 sGRBs and four hybrid GRBs that are long GRBs accompanied by kilonova-like signals. We use this data to evaluate their physical properties through spectral energy distribution fitting. To assess the characteristics of the host galaxies, we utilized a volume-limited sample (z < 0.5) from the COSMOS field as a control group. Contrary to expectations that the BNS merger rate is proportional to host stellar mass, the short and hybrid GRB population appears less massive than the mass-weighted distribution of the control sample. Instead, we propose a formulation for the expected BNS merger rate from a galaxy as log(nBNSGyr-1)=0.86×log(M*/M⊙) + 0.44×log(sSFRyr-1)+0.857 , which optimally explains the deviation between the stellar mass distributions of the GRB host galaxies and the control sample. These insights provide a strategic framework for targeted GW follow-ups and enhance our ability to identify potential host galaxies for future GW events.

Inferring Binary Properties from Gravitational-Wave Signals

This review provides a conceptual and technical survey of methods for parameter estimation of gravitational-wave signals in ground-based interferometers such as Laser Interferometer Gravitational-Wave Observatory (LIGO) and Virgo. We introduce the framework of Bayesian inference and provide an overview of models for the generation and detection of gravitational waves from compact binary mergers, focusing on the essential features that are observable in the signals. Within the traditional likelihood-based paradigm, we describe various approaches for enhancing the efficiency and robustness of parameter inference. This includes techniques for accelerating likelihood evaluations, such as heterodyne/relative binning, reduced-order quadrature, multibanding, and interpolation. We also cover methods to simplify the analysis to improve convergence, via reparameterization, importance sampling, and marginalization. We end with a discussion of recent developments in the application of likelihood-free (simulation-based) inference methods to gravitational-wave data analysis.

Navigating Unknowns: Deep Learning Robustness for Gravitational-wave Signal Reconstruction

We present a rapid and reliable deep-learning-based method for gravitational-wave (GW) signal reconstruction from elusive, generic binary black hole mergers in LIGO data. We demonstrate that our model, AWaRe, effectively recovers GWs with parameters it has not encountered during training. This includes features like higher black hole masses, additional harmonics, eccentricity, and varied waveform systematics, which introduce complex modulations in the waveform’s amplitudes and phases. The accurate reconstructions of these unseen signal characteristics demonstrate AWaRe's ability to handle complex features in the waveforms. By directly incorporating waveform reconstruction uncertainty estimation into the AWaRe framework, we show that for real GW events, the uncertainties in AWaRe's reconstructions align closely with those achieved by benchmark algorithms like BayesWave and coherent WaveBurst. The robustness of our model to real GW events and its ability to extrapolate to unseen data open new avenues for investigations in various aspects of GW astrophysics and data analysis, including tests of general relativity and the enhancement of current GW search methodologies.

Detection of the Fe K lines from the binary AGN in 4C+37.11

We report the discovery of the Fe K line emission at $ $ keV with a width of $ $ keV using two epochs of Chandra archival data for the nucleus of the galaxy 4C+37.11, which is known to host a binary supermassive black hole (BSMBH) system where the SMBHs are separated by $ mas or sim 7pc. Our study reports the first detection of the Fe K line from a known binary AGN, which has an F-statistic value of 20.98 and a probability of $2.47 $. Stacking two spectra reveals another Fe K line component at $ $ keV. Different model scenarios indicate that the lines originate from the combined effects of accretion disk emission and circumnuclear collisionally ionized medium. The observed low column density favors a gas-poor merger scenario, where the high temperature of the hot ionized medium may be associated with the shocked gas in the binary merger and not with star formation activity. The estimated total BSMBH mass and disk inclination are $ $ M$_ and $ indicating that the BSMBH is probably a high-inclination system. We were not able to tightly constrain the spin parameter using the present data sets. Our results draw attention to the fact that detecting the Fe K line emissions from BSMBHs is important for estimating the individual SMBH masses and the spins of the binary SMBHs, as well as for exploring their emission regions.

The Detection Prospect of the Counter Jet Radiation in the Late Afterglow of GRB 170817A

The central engine of a gamma-ray burst (GRB) is widely believed to launch a pair of oppositely moving jets, i.e., the forward jet moving toward us and the counter jet regressing away. The forward jet generates the radiation typically observed in GRBs, while the counter jet has not been detected yet due to its dimness. GRB 170817A, a short burst associated with a binary neutron star merger event, is a nearby event (z = 0.0097) with an off-axis structured energetic forward jet and hence probably the most suitable target for searching the counter jet radiation. Assuming the same properties for the forward and counter jet components as well as the shock parameters, the fit to the multiwavelength afterglow emission of GRB 170817A suggests a peak time ∼quite a few ×103 day of the counter jet radiation, but the detection prospect of this new component is not promising. Anyhow, if the shock parameters (ϵ e and ϵ B) of the counter jet component are (a few times) higher than that of the forward shock, e.g., the posterior results of the magnetic energy fraction of the forward shock and the counter jet are log10ϵB =−4.23−0.69+1.42 and log10ϵB,cj =−3.65−2.11+3.06 , respectively, as still allowed by the current data, the counter jet afterglow emission will be enhanced and hence may be detected. A few hours of exposure by JWST in the F356W band will stringently test such a scenario.

Waveform systematics in gravitational-wave inference of signals from binary neutron star merger models incorporating higher-order modes information

Waveform systematics in gravitational-wave inference of signals from binary neutron star merger models incorporating higher-order modes information

Probing general relativistic spin–orbit coupling with gravitational waves from hierarchical triple systems

ABSTRACT Wave packets propagating in inhomogeneous media experience a coupling between internal and external degrees of freedom and, as a consequence, follow spin-dependent trajectories. These phenomena, well known in optics and condensed matter physics, are referred to as spin Hall effects. Similarly, the gravitational spin Hall effect is expected to affect the propagation of gravitational waves on curved spacetimes. In this general-relativistic setup, the curvature of spacetime acts as impurities in a semiconductor or inhomogeneities in an optical medium, leading to a frequency- and polarization-dependent propagation of wave packets. In this letter, we study this effect for strong-field lensed gravitational waves generated in hierarchical triple black hole systems in which a stellar-mass binary merges near a more massive black hole. We calculate how the gravitational spin Hall effect modifies the gravitational waveforms and show its potential for experimental observation. If detected, these effects will have profound implications for astrophysics and tests of general relativity.

Using equation of state constraints to classify low-mass compact binary mergers

Using equation of state constraints to classify low-mass compact binary mergers

MURCA driven bulk viscosity in neutrino trapped baryonic matter

We examine bulk viscosity, taking into account trapped neutrinos in baryonic matter, in the context of binary neutron star mergers. Following the merging event, the binary star can yield a remnant compact object with densities up to 5 nuclear saturation density and temperature up to 50 MeV resulting in the retention of neutrinos. We employ two relativistic mean field models, NL3 and DDME2, to describe the neutrino-trapped baryonic matter. The dissipation coefficient is determined by evaluating the Modified URCA interaction rate in the dense baryonic medium, and accounting for perturbations caused by density oscillations. We observe the resonant behavior of bulk viscosity as it varies with the temperature of the medium. The bulk viscosity peak remains within the temperature range of ∼13\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\sim 13$$\\end{document}–50 MeV, depending upon the underlying equation of states and lepton fractions. This temperature range corresponds to the relevant domain of binary neutron star mergers. We also note that in presence of neutrinos in the medium the bulk viscosity peak shifts towards higher temperature and the peak value of bulk viscosity also changes. The time scale of viscous dissipation is dictated by the beta-off-equilibrium susceptibilities derived from the nuclear equation of state. The resulting viscous decay time scale ranges from 32–100 ms, which aligns with the order of magnitude of the post-merger object’s survival time in some specific scenarios.

New results on the gamma-ray burst variability--luminosity relation

At the dawn of the gamma--ray burst (GRB) afterglow era, a Cepheid-like correlation was discovered between the time variability $V$ and the isotropic-equivalent peak luminosity $L_ iso $ of the prompt emission of about a dozen long GRBs with measured redshift available at that time. Soon afterwards, the correlation was confirmed in a sample of about 30 GRBs, even though it was affected by significant scatter. Unlike the minimum variability timescale (MVT), $V$ measures the relative power of short-to-intermediate timescales. We aim to test the correlation using about 200 long GRBs with spectroscopically measured redshift, detected by Swift Fermi and Konus/ WIND for which both observables can be accurately estimated. The variability for all selected GRBs was calculated according to the original definition using the 64 ms background-subtracted light curves of Swift /BAT ( Fermi /GBM) in the 15--150 (8--900) keV energy passband. Peak luminosities were either taken from the literature or derived from modelling broad-band spectra acquired with either Konus/ WIND or Fermi /GBM. The statistical significance of the correlation has weakened to $ 2$&lt;!PCT!&gt;, mostly due to the appearance of a number of smooth and luminous GRBs that are characterised by a relatively small $V$. At odds with most long GRBs, three out of four long-duration merger candidates have high $V$ and low iso The luminosity is more tightly connected with shortest timescales measured by MVT than the short to intermediate timescales measured by $V$. We discuss the implications for internal dissipation models and the role of the $e^ pm $ photosphere. We identified a few smooth GRBs with a single broad pulse and low $V$ that might have an external shock origin, in contrast with most GRBs. The combination of high variability $(V 0.1)$, low luminosity iso $ erg s$^ $, and short MVT ($ 0.1$ s) could be a good indicator for a compact binary merger origin.

Binary Mergers in the Centers of Galaxies: Synergy between Stellar Flybys and Tidal Fields

Galactic centers (GCs) are very dynamically active environments, often harboring a nuclear star cluster and supermassive black hole at their cores. Binaries in these environments are subject to strong tidal fields that can efficiently torque its orbit, exciting near-unity eccentricities that ultimately lead to their merger. In turn, frequent close interactions with passing stars impulsively perturb the orbit of the binary, generally softening their orbit until their evaporation, potentially hindering the role of tides to drive these mergers. In this work, we study the evolution of compact object binaries in the GC and their merger rates, focusing for the first time on the combined effect of the cluster’s tidal field and flyby interactions. We find a significant synergy between both processes, where merger rates increase by a factor of ∼10−30 compared to models in which only flybys or tides are taken into account. This synergy is a consequence of the persistent tides-driven eccentricity excitation that is enhanced by the gradual diffusion of j z driven by flybys. The merger efficiency peaks when the diffusion rate is ∼10−100 slower than the tides-driven torquing. Added to this synergy, we also find that the gradual softening of the binary can lift the relativistic quenching of initially tight binaries, otherwise unable to reach extreme eccentricities, and thus expanding the available phase space for mergers. Cumulatively, we conclude that despite the gradual softening of binaries due to flybys, these greatly enhance their merger rates in GCs by promoting the tidal-field-driven eccentricity excitation.

Bayesian inference of multi-messenger astrophysical data: Joint and coherent inference of gravitational waves and kilonovae

Context. Multi-messenger observations of binary neutron star mergers can provide information on the neutron star’s equation of state (EOS) above the nuclear saturation density by directly constraining the mass-radius diagram. Aims. We present a Bayesian framework for joint and coherent analyses of multi-messenger binary neutron star signals. As a first application, we analyze the gravitational-wave GW170817 and the kilonova (kN) AT2017gfo data. These results are then combined with the most recent X-ray pulsar analyses of PSR J0030+0451 and PSR J0740+6620 to obtain new EOS constraints. Methods. We extend the bajes infrastructure with a joint likelihood for multiple datasets, support for various semi-analytical kN models, and numerical-relativity (NR)-informed relations for the mass ejecta, as well as a technique to include and marginalize over modeling uncertainties. The analysis of GW170817 used the TEOBResumS effective-one-body waveform template to model the gravitational-wave signal. The analysis of AT2017gfo used a baseline multicomponent spherically symmetric model for the kN light curves. Various constraints on the mass-radius diagram and neutron star properties were then obtained by resampling over a set of ten million parameterized EOSs, which was built under minimal assumptions (general relativity and causality). Results. We find that a joint and coherent approach improves the inference of the extrinsic parameters (distance) and, among the intrinsic parameters, the mass ratio. The inclusion of NR-informed relations marks a strong improvement over the case in which an agnostic prior is used on the intrinsic parameters. Comparing Bayes factors, we find that the two observations are better explained by the common source hypothesis only by assuming NR-informed relations. These relations break some of the degeneracies in the employed kN models. The EOS inference folding-in PSR J0952-0607 minimum-maximum mass, PSR J0030+0451 and PSR J0740+6620 data constrains, among other quantities, the neutron star radius to R1.4TOV = 12.30− 0.56+ 0.81 km(R1.4TOV = 13.20− 0.90+ 0.91 km) and the maximum mass to MmaxTOV = 2.28− 0.17+ 0.25M⊙(MmaxTOV = 2.32− 0.19+ 0.30M⊙), where the ST+PDT (PDT-U) analysis of Vinciguerra et al. (2024, ApJ, 961, 62) for PSR J0030+0451 was employed. Hence, the systematics on the PSR J0030+0451 data reduction currently dominate the mass-radius diagram constraints. Conclusions. We conclude that bajes delivers robust analyses in line with other state-of-the-art results in the literature. Strong EOS constraints are provided by pulsars observations, albeit with large systematics in some cases. Current gravitational-wave constraints are compatible with pulsar constraints and can further improve the latter.

Close Encounters of Wide Binaries Induced by the Galactic Tide: Implications for Stellar Mergers and Gravitational-wave Sources

A substantial fraction of stars can be found in wide binaries with projected separations between ∼102 and 105 au. In the standard lore of binary physics, these would evolve as effectively single stars that remotely orbit one another on stationary Keplerian ellipses. However, embedded in their Galactic environment, the low binding energy of wide binaries makes them exceptionally prone to perturbations from the gravitational potential of the Milky Way and encounters with passing stars. Employing a fully relativistic N-body integration scheme, we study the impact of these perturbations on the orbital evolution of wide binaries along their trajectory through the Milky Way. Our analysis reveals that the torques exerted by the Galaxy can cause large-amplitude oscillations of the binary eccentricity to 1 − e ≲ 10−8. As a consequence, the wide binary members pass close to each other at periapsis, which, depending on the type of binary, potentially leads to a mass transfer or collision of stars or to an inspiral and subsequent merger of compact remnants due to gravitational-wave radiation. Based on a simulation of 105 wide binaries across the Galactic field, we find that this mechanism could significantly contribute to the rate of stellar collisions and binary black hole mergers as inferred from observations of luminous red novae and gravitational-wave events by LIGO/Virgo/Kagra. We conclude that the dynamics of wide binaries, despite their large mean separation, can give rise to extreme interactions between stars and compact remnants.

Primordial black hole interpretation in subsolar mass gravitational wave candidate SSM200308

In the recent second part of the third observation run by the LIGO-Virgo-KAGRA collaboration, a candidate with sub-solar mass components was reported, which we labelled as SSM200308. This study investigates the premise that primordial black holes (PBHs), arising from Gaussian perturbation collapses, could explain SSM200308. Through Bayesian analysis, we obtain the primordial curvature power spectrum that leads to the merger rate of PBHs aligning with observational data as long as they constitute f PBH = 5.66+58.68 -5.44 × 10-2 of the dark matter. However, while the gravitational wave (GW) background from binary PBH mergers is within current observational limits, the scalar-induced GWs associated with PBH formation exceed the constraints imposed by pulsar timing arrays, challenging the Gaussian perturbation collapse PBH model as the source of SSM200308.

Unequal-mass highly spinning binary black hole mergers in the stable mass transfer formation channel

Unequal-mass highly spinning binary black hole mergers in the stable mass transfer formation channel

Gauss-Bonnet Cosmology: large-temperature behaviour and bounds from Gravitational Waves

We provide a transparent discussion of the high temperature asymptotic behaviour of Cosmology in a dilaton-Einstein-Gauss-Bonnet (dEGB) scenario of modified gravity with vanishing scalar potential. In particular, we show that it has a clear interpretation in terms of only three attractors (stable critical points) of a set of autonomous differential equations: w = -1/3, w = 1 and 1 < w < 7/3, where w ≡ p/ρ is the equation of state, defined as the ratio of the total pressure and the total energy density. All the possible different high-temperature evolution histories of the model are exhausted by only eight paths in the flow of the set of the autonomous differential equations. Our discussion clearly explains why five out of them are characterized by a swift transition of the system toward the attractor, while the remaining three show a more convoluted evolution, where the system follows a meta-stable equation of state at intermediate temperatures before eventually jumping to the real attractor at higher temperatures. Compared to standard Cosmology, the regions of the dEGB parameter space with w = -1/3 show a strong enhancement of the expected Gravitational Wave stochastic background produced by the primordial plasma of relativistic particles of the Standard Model. This is due to the very peculiar fact that dEGB allows to have an epoch when the energy density ρ rad of the relativistic plasma dominates the energy of the Universe while at the same time the rate of dilution with T of the total energy density is slower than what usually expected during radiation dominance. This allows to use the bound from Big Bang Nucleosynthesis (BBN) to put in dEGB a constraint T RH ≲ (108 – 109) GeV on the reheating temperature of the Universe T RH. Such BBN bound is complementary to late-time constraints from compact binary mergers.

The Pristine Inner Galaxy Survey (PIGS)

The most metal-poor stars provide valuable insights into the early chemical enrichment history of a system, carrying the chemical imprints of the first generations of supernovae. The most metal-poor region of the Sagittarius dwarf galaxy remains inadequately observed and characterised. To date, only ∼4 stars with [Fe/H] &lt; −2.0 have been chemically analysed with high-resolution spectroscopy. In this study, we present the most extensive chemical abundance analysis of 12 low-metallicity stars with metallicities down to [Fe/H] = −3.26 and located in the main body of Sagittarius. These targets, selected from the Pristine Inner Galaxy Survey, were observed using the MIKE high-resolution spectrograph at the Magellan-Clay telescope, which allowed us to measure up to 17 chemical species. The chemical composition of these stars reflects the imprint of a variety of type II supernovae (SNe II). A combination of low- to intermediate-mass high-energy SNe and hypernovae (∼10 − 70 M⊙) is required to account for the abundance patterns of the lighter elements up to the Fe-peak. The trend of the heavy elements suggests the involvement of compact binary merger events and fast-rotating (up to ∼300 km s−1) intermediate-mass to massive metal-poor stars (∼25 − 120 M⊙) that are the sources of rapid and slow processes, respectively. Additionally, asymptotic giant branch stars contribute to a wide dispersion of [Ba/Mg] and [Ba/Eu]. The absence of an α−knee in our data indicates that type Ia supernovae did not contribute in the very metal-poor region ([Fe/H] ≤ −2.0). However, they might have started to pollute the interstellar medium at [Fe/H] &gt; −2.0, given the relatively low [Co/Fe] in this metallicity region.

Probing the black holes in a dark matter halo of M87 using gravitational wave echoes

Variations at the event horizon structure of a black hole will emit the signals of the gravitational wave echoes associated with the ringdown of the binary black hole merger. In this work, combining mass model of M87 and Einasto profile for dark matter halo, we construct one formal solution of black holes in dark matter halo, and this solution includes the regular black hole for geometric parameter a<2M\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$a < 2M$$\\end{document} and the wormhole for a≥2M\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$a \\ge 2M$$\\end{document}. Then, we test this solution under axial gravitational perturbation and calculate their quasinormal modes (QNM). Our results show that when geometric parameter a>2M\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$a>2M$$\\end{document}, a series of gravitational wave echoes appear after the QNM, and one distinctive feature of gravitational wave echoes is the double barrier. Besides, we also study the impacts of shape parameters α\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\alpha $$\\end{document} of Einasto profile both on the QNM and gravitational wave echoes, and extract their frequencies using WKB method and Prony method. In principle, with the increasing of shape parameter α\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\alpha $$\\end{document}, the frequencies of the QNM and gravitational wave echoes between the different shape parameters must be different. But now there is zero difference of the frequencies between shape parameter α=0.18\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\alpha =0.18$$\\end{document} and α=0.20\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\alpha =0.20$$\\end{document}. Zero difference is not allowed, so we give an upper limit for shape parameter α\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\alpha $$\\end{document}, which is approximately α<0.18\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\alpha <0.18$$\\end{document}. These signals of gravitational wave echoes in a dark matter halo may be detected in the near future, and these characteristics of gravitational wave echoes may serve as local measurements of dark matter.

Afterglows from Binary Neutron Star Postmerger Systems Embedded in Active Galactic Nuclei Disks

The observability of afterglows from binary neutron star mergers occurring within active galactic nuclei (AGN) disks is investigated. We perform 3D GRMHD simulations of a postmerger system and follow the jet launched from the compact object. We use semianalytic techniques to study the propagation of the blast wave powered by the jet through an AGN disk-like external environment, extending to distances beyond the disk scale height. The synchrotron emission produced by the jet-driven forward shock is calculated to obtain the afterglow emission. The observability of this emission at different frequencies is assessed by comparing it to the quiescent AGN emission. In the scenarios where the afterglow could temporarily outshine the AGN, we find that detection will be more feasible at higher frequencies (≳1014 Hz) and the electromagnetic counterpart could manifest as a fast variability in the AGN emission, on timescales less than a day.

Impact of lensing of gravitational waves on the observed distribution of neutron star masses

The distribution of masses of neutron stars, particularly the maximum mass value, is considered a probe of their formation, evolution and internal physics (i.e., equation of state). This mass distribution could in principle be inferred from the detection of gravitational waves from binary neutron star mergers. Using mock catalogues of 105 dark sirens events, expected to be detected by Einstein Telescope over an operational period of , we show how the biased luminosity distance measurement induced by gravitational lensing affects the inferred redshift and mass of the merger. This results in higher observed masses than expected. Up to 2% of the events are predicted to fall above the maximum allowed neutron star mass depending on the intrinsic mass distribution and signal-to-noise ratio threshold adopted. The underlying true mass distribution and maximum mass could still be approximately recovered in the case of bright standard sirens.