Multi-messenger Observations Research Articles

Disk-induced Binary Precession: Implications for Dynamics and Multimessenger Observations of Black Hole Binaries

Many studies have recently documented the orbital response of eccentric binaries accreting from thin circumbinary disks, characterizing the change in the binary semimajor axis and eccentricity. We extend these calculations to include the precession of the binary’s longitude of periapse induced by the circumbinary disk, and we characterize this precession continuously with binary eccentricity e b for equal mass components. This disk-induced apsidal precession is prograde with a weak dependence on the binary eccentricity when e b ≲ 0.4 and decreases approximately linearly for e b ≳ 0.4; yet at all e b binary precession is faster than the rates of change to the semimajor axis and eccentricity by an order of magnitude. We estimate that such precession effects are likely most important for subparsec separated binaries with masses ≲107 M ⊙, like LISA precursors. We find that accreting, equal-mass LISA binaries with M < 106 M ⊙ (and the most massive M ∼ 107 M ⊙ binaries out to z ∼ 3) may acquire a detectable phase offset due to the disk-induced precession. Moreover, disk-induced precession can compete with general relativistic precession in a vacuum, making it important for observer-dependent electromagnetic searches for accreting massive binaries—like Doppler boost and binary self-lensing models—after potentially only a few orbital periods.

Shining light on the hosts of the nano-Hertz gravitational wave sources: a theoretical perspective

ABSTRACT The formation of supermassive black holes (SMBHs) in the Universe and its role in the properties of the galaxies is one of the open questions in astrophysics and cosmology. Though, traditionally, electromagnetic waves have been instrumental in direct measurements of SMBHs, significantly influencing our comprehension of galaxy formation, gravitational waves (GW) bring an independent avenue to detect numerous binary SMBHs in the observable Universe in the nano-Hertz range using the pulsar timing array observation. This brings a new way to understand the connection between the formation of binary SMBHs and galaxy formation if we can connect theoretical models with multimessenger observations namely GW data and galaxy surveys. Along these lines, we present here the first paper on this series based on romulus25 cosmological simulation on the properties of the host galaxies of SMBHs and propose on how this can be used to connect with observations of nano-Hertz GW signal and galaxy surveys. We show that the most dominant contribution to the background will arise from sources with high chirp masses which are likely to reside in low-redshift early-type galaxies with high stellar mass, largely old stellar population, and low star formation rate, and that reside at centres of galaxy groups and manifest evidence of recent mergers. The masses of the sources show a correlation with the halo mass and stellar mass of the host galaxies. This theoretical study will help in understanding the host properties of the GW sources and can help in establishing a connection with observations.

Constraints on the proton fraction of cosmic rays at the highest energies and the consequences for cosmogenic neutrinos and photons

Over the last decade, observations have shown that the mean mass of ultra-high-energy cosmic rays (UHECRs) increases progressively toward the highest energies. However, the precise composition is still unknown and several theoretical studies hint at the existence of a subdominant proton component up to the highest energies. Motivated by the exciting prospect of performing charged-particle astronomy with ultra-high-energy (UHE) protons we quantify the level of UHE-proton flux that is compatible with present multimessenger observations and the associated fluxes of neutral messengers produced in the interactions of the protons. We study this scenario with numerical simulations of two independent populations of extragalactic sources and perform a fit to the combined UHECR energy spectrum and composition observables, constrained by diffuse gamma-ray and neutrino observations. We find that up to of order 10% of the cosmic rays at the highest energies can be UHE protons, although the result depends critically on the selected hadronic interaction model for the air showers. Depending on the maximum proton energy (E max p) and the redshift evolution of sources, the associated flux of cosmogenic neutrinos and UHE gamma rays can significantly exceed the multimessenger signal of the mixed-mass cosmic rays. Moreover, if E max p is above the GZK limit, we predict a large flux of UHE neutrinos above EeV energies that is absent in alternate scenarios for the origin of UHECRs. We present the implications and opportunities afforded by these UHE proton, neutrino and photon fluxes for future multimessenger observations.

Could a Kilonova Kill: A Threat Assessment

Binary neutron star mergers produce high-energy emissions from several physically different sources, including a gamma-ray burst (GRB) and its afterglow, a kilonova (KN), and, at late times, a remnant many parsecs in size. Ionizing radiation from these sources can be dangerous for life on Earth-like planets when located too close. Work to date has explored the substantial danger posed by the GRB to on-axis observers; here we focus instead on the potential threats posed to nearby off-axis observers. Our analysis is based largely on observations of the GW170817/GRB 170817A multi-messenger event, as well as theoretical predictions. For baseline KN parameters, we find that the X-ray emission from the afterglow may be lethal out to ∼1 pc and the off-axis gamma-ray emission may threaten a range out to ∼4 pc, whereas the greatest threat comes years after the explosion, from the cosmic rays accelerated by the KN blast, which can be lethal out to distances up to ∼11 pc. The distances quoted here are typical, but the values have significant uncertainties and depend on the viewing angle, ejected mass, and explosion energy in ways we quantify. Assessing the overall threat to Earth-like planets, KNe have a similar kill distance to supernovae, but are far less common. However, our results rely on the scant available KN data, and multi-messenger observations will clarify the danger posed by such events.

On the Source Contribution to the Galactic Diffuse Gamma Rays above 398 TeV Detected by the Tibet ASγ Experiment

Potential contribution from gamma-ray sources to the Galactic diffuse gamma rays observed above 100 TeV (sub-PeV energy range) by the Tibet ASγ experiment is an important key to interpreting recent multimessenger observations. This paper reveals a surprising fact: none of the 23 Tibet ASγ diffuse gamma-ray events above 398 TeV within the Galactic latitudinal range of ∣b∣ < 10° come from the 43 sub-PeV gamma-ray sources reported in the 1LHAASO catalog, which proves that these sources are not the origins of the Tibet ASγ diffuse gamma-ray events. No positional overlap between the Tibet ASγ diffuse gamma-ray events and the sub-PeV LHAASO sources currently supports the diffusive nature of the Tibet ASγ diffuse gamma-ray events, although their potential origin in the gamma-ray sources yet unresolved in the sub-PeV energy range cannot be ruled out.

Real-time monitoring for the next core-collapse supernova in JUNO

The core-collapse supernova (CCSN) is considered one of the most energetic astrophysical events in the universe. The early and prompt detection of neutrinos before (pre-SN) and during the supernova (SN) burst presents a unique opportunity for multi-messenger observations of CCSN events. In this study, we describe the monitoring concept and present the sensitivity of the system to pre-SN and SN neutrinos at the Jiangmen Underground Neutrino Observatory (JUNO), a 20 kton liquid scintillator detector currently under construction in South China.The real-time monitoring system is designed to ensure both prompt alert speed and comprehensive coverage of progenitor stars. It incorporates prompt monitors on the electronic board as well as online monitors at the data acquisition stage.Assuming a false alert rate of 1 per year, this monitoring system exhibits sensitivity to pre-SN neutrinos up to a distance of approximately 1.6 (0.9) kiloparsecs and SN neutrinos up to about 370 (360) kiloparsecs for a progenitor mass of 30 solar masses, considering both normal and inverted mass ordering scenarios.The pointing ability of the CCSN is evaluated by analyzing the accumulated event anisotropy of inverse beta decay interactions from pre-SN or SN neutrinos. This, along with the early alert, can play a crucial role in facilitating follow-up multi-messenger observations of the next galactic or nearby extragalactic CCSN.

Premerger Sky Localization of Gravitational Waves from Binary Neutron Star Mergers Using Deep Learning

The simultaneous observation of gravitational waves (GW) and prompt electromagnetic counterparts from the merger of two neutron stars can help reveal the properties of extreme matter and gravity during and immediately after the final plunge. Rapid sky localization of these sources is crucial to facilitate such multimessenger observations. As GWs from binary neutron star (BNS) mergers can spend up to 10–15 minutes in the frequency bands of the detectors at design sensitivity, early-warning alerts and premerger sky localization can be achieved for sufficiently bright sources, as demonstrated in recent studies. In this work, we present premerger BNS sky localization results using GW-SkyLocator, a deep-learning model capable of inferring sky location posterior distributions of GW sources at orders of magnitude faster speeds than standard Markov Chain Monte Carlo methods. We test our model’s performance on a catalog of simulated injections from Sachdev, recovered at 0–60 s before the merger, and obtain comparable sky localization areas to the rapid localization tool BAYESTAR. These results show the feasibility of our model for premerger sky localization and the possibility of follow-up observations for precursor emissions from BNS mergers.

Neutrino production in blazar radio cores

Models of the origin of astrophysical neutrinos with energies from TeVs to PeVs are strongly constrained by multimessenger observations and population studies. Recent results point to statistically significant associations between these neutrinos and active galactic nuclei (AGN) selected by their radio flux observed with very-long-baseline interferometry (VLBI). This suggests that the neutrinos are produced in central parsecs of blazars, AGN with relativistic jets pointing to the observer. However, conventional AGN models tend to explain only the highest-energy part of the neutrino flux observationally associated with blazars. Here we discuss in detail how the neutrinos can be produced in the part of an AGN giving the dominant contribution to the VLBI radio flux, the radio core located close to the jet base. Physical conditions there differ both from the immediate environment of the central black hole and from the plasma blobs moving along the jet. Required neutrino fluxes, considerably smaller than those of photons, can be produced in interactions of relativistic protons, accelerated closer to the black hole, with radiation in the core.

High-energy neutrinos from merging stellar-mass black holes in active galactic nuclei accretion disc

ABSTRACT A population of binary stellar-mass black hole (BBH) mergers are believed to occur embedded in the accretion disc of active galactic nuclei (AGNs). In this Letter, we demonstrate that the jets from these BBH mergers can propagate collimatedly within the disc atmosphere along with a forward shock and a reverse shock forming at the jet head. Efficient proton acceleration by these shocks is usually expected before the breakout, leading to the production of TeV−PeV neutrinos through interactions between these protons and electron-radiating photons via photon–meson production. AGN BBH mergers occurring in the outer regions of the disc are more likely to produce more powerful neutrino bursts. Taking the host AGN properties of the potential GW190521 electromagnetic (EM) counterpart as an example, one expects ≳1 neutrino events detectable by IceCube if the jet is on-axis and the radial location of the merger is R ≳ 105Rg, where Rg is the gravitational radius of the supermassive BH. Neutrino bursts from AGN BBH mergers could be detected by IceCube following the observation of gravitational waves (GWs), serving as precursor signals before the detection of EM breakout signals. AGN BBH mergers are potential target sources for future joint GW, neutrino, and EM multi-messenger observations.

What Is the Nature of the HESS J1731-347 Compact Object?

Once further confirmed in future analyses, the radius and mass measurement of HESS J1731-347 with and will be among the lightest and smallest compact objects ever detected. This raises many questions about its nature and opens up the window for different theories to explain such a measurement. In this article, we use the information from Doroshenko et al. on the mass, radius, and surface temperature together with the multimessenger observations of neutron stars to investigate the possibility that HESS J1731-347 is one of the lightest observed neutron star, a strange quark star, a hybrid star with an early deconfinement phase transition, or a dark matter–admixed neutron star. The nucleonic and quark matter are modeled within realistic equation of states (EOSs) with a self-consistent calculation of the pairing gaps in quark matter. By performing the joint analysis of the thermal evolution and mass–radius constraint, we find evidence that within a 1σ confidence level, HESS J1731-347 is consistent with the neutron star scenario with the soft EOS as well as with a strange and hybrid star with the early deconfinement phase transition with a strong quark pairing and neutron star admixed with dark matter.

Observing white dwarf tidal stripping with TianQin gravitational wave observatory

ABSTRACT Recently discovered regular X-ray bursts known as quasi-periodic eruptions have a proposed model that suggests a tidal stripping white dwarf inspiralling into the galaxy’s central black hole on an eccentric orbit. According to this model, the interaction of the stripping white dwarf with the central black hole would also emit gravitational wave signals, their detection can help explore the formation mechanism of quasi-periodic eruptions and facilitate multimessenger observations. In this paper, we investigated the horizon distance of TianQin on this type of gravitation wave signal and found it can be set to 200 Mpc. We also find that those stripping white dwarf model sources with central black hole mass within $10^4 \!-\! 10^{5.5}\, \mathrm{M}_\odot$ are more likely to be detected by TianQin. We assessed the parameter estimation precision of TianQin on those stripping white dwarf model sources. Our result shows that, even in the worst case, TianQin can determine the central black hole mass, the white dwarf mass, the central black hole spin, and the orbital initial eccentricity with a precision of 10−2. In the optimistic case, TianQin can determine the central black hole mass and the white dwarf mass with a precision of 10−7, determine the central black hole spin with a precision of 10−5, and determine the orbital initial eccentricity with a precision of 10−8. Moreover, TianQin can determine the luminosity distance with a precision of 10−1 and determine the sky localization with a precision of 10−2–10 $\rm deg^2$.

Decomposing the Origin of TeV–PeV Emission from the Galactic Plane: Implications of Multimessenger Observations

High-energy neutrino and γ-ray emission has been observed from the Galactic plane, which may come from individual sources and/or diffuse cosmic rays. We evaluate the contribution of these two components through the multimessenger connection between neutrinos and γ-rays in hadronic interactions. We derive maximum fluxes of neutrino emission from the Galactic plane using γ-ray catalogs, including 4FGL, HGPS, 3HWC, and 1LHAASO, and measurements of the Galactic diffuse emission by Tibet ASγ and LHAASO. We find that the IceCube Galactic neutrino flux is larger than the contribution from all resolved sources when excluding promising leptonic sources such as pulsars, pulsar wind nebulae, and TeV halos. Our result indicates that the Galactic neutrino emission is likely dominated by the diffuse emission by the cosmic-ray sea and unresolved hadronic γ-ray sources. In addition, the IceCube flux is comparable to the sum of the flux of nonpulsar sources and the LHAASO diffuse emission especially above ∼30 TeV. This implies that the LHAASO diffuse emission may dominantly originate from hadronic interactions, either as the truly diffuse emission or unresolved hadronic emitters. Future observations of neutrino telescopes and air-shower γ-ray experiments in the Southern hemisphere are needed to accurately disentangle the source and diffuse emission of the Milky Way.

What constraints can one pose on the maximum mass of neutron stars from multimessenger observations?

ABSTRACT The maximum mass of neutron stars (MTOV) plays a crucial role in understanding their equation of state (EoS). Previous studies have used the measurements for the compactness of massive pulsars and the tidal deformability of neutron stars in binary neutron star (BNS) mergers to constrain the EoS and thus the MTOV. The discovery of the most massive pulsar, PSR J0952−0607, with a mass $\sim 2.35\, {\rm M}_{\odot }$, has provided a valuable lower limit for MTOV. Another efficient method to constrain MTOV is by examining the type of central remnant formed after a BNS merger. Gravitational wave (GW) data can provide the total mass of the system, while accompanying electromagnetic signals can help infer the remnant type. In this study, we combine all the previous constraints and utilize the observational facts that about 24 per cent of the short gamma-ray bursts are followed by an X-ray internal plateau, which indicate that roughly this fraction of BNS mergers yield supermassive neutron stars, to perform (Markov Chain) Monte Carlo simulations. These simulations allow us to explore the probability density distribution of MTOV and other parameters related to BNS mergers. Our findings suggest that MTOV is likely around $2.49\!-\!2.52\, {\rm M}_{\odot }$, with an uncertainty range of approximately [$-0.16$, $0.15\, {\rm M}_{\odot }$] ([$-0.28$, $0.26\, {\rm M}_{\odot }$]) at 1σ (2σ) confidence level. Furthermore, we examine the type of merger remnants in specific events like GW170817 and GW190425 to further constrain MTOV and other relevant parameters, which can help to understand the physical processes involved in BNS mergers.

Impact of modelling galaxy redshift uncertainties on the gravitational-wave dark standard siren measurement of the Hubble constant

ABSTRACT Gravitational wave science is a new and rapidly expanding field of observational astronomy. Multimessenger observations of the binary neutron star merger GW170817 have provided some iconic results including the first gravitational-wave standard siren measurement of the Hubble constant, opening up a new way to probe cosmology. The majority of the compact binary sources observed in gravitational waves are, however, without bright electromagnetic counterparts. In these cases, one can fall back on the ‘dark standard siren’ approach to include information statistically from potential host galaxies. For such a measurement, we need to be cautious about all possible sources of systematic errors. In this paper, we begin to study the possible errors coming from the galaxy catalogue sector, and in particular, look into the effect of galaxy redshift uncertainties for the cases where these are photometry-based. We recalculate the dark standard siren Hubble constant using the latest third Gravitational-Wave Transient Catalog (GWTC-3) events and associated galaxy catalogues, with different galaxy redshift uncertainty models, namely, the standard Gaussian, a modified Lorentzian, and no uncertainty at all. We find that not using redshift uncertainties at all can lead to a potential bias comparable with other potential systematic effects previously considered for the GWTC-3 H0 measurement (however still small compared to the overall statistical error in this measurement). The difference between different uncertainty models leads to small differences in the results for the current data; their impact is much smaller than the current statistical errors and other potential sources of systematic errors which have been considered in previous robustness studies.

Pre-merger alert to detect prompt emission in very-high-energy gamma-rays from binary neutron star mergers: Einstein Telescope and Cherenkov Telescope Array synergy

The current generation of very-high-energy gamma-ray (VHE; E &gt; 30 GeV) detectors (MAGIC and H.E.S.S.) have recently demonstrated the ability to detect the afterglow emission of gamma-ray bursts (GRBs). However, the GRB prompt emission, typically observed in the 10 keV–10 MeV band, is still undetected at higher energies. Here, we investigate the perspectives of multi-messenger observations to detect the earliest VHE emission from short GRBs. Considering binary neutron star mergers as progenitors of short GRBs, we evaluate the joint detection efficiency of the Cherenkov Telescope Array (CTA) observing in synergy with the third generation of gravitational-wave detectors, such as the Einstein Telescope (ET) and Cosmic Explorer (CE). In particular, we evaluate the expected capabilities to detect and localize gravitational-wave events in the inspiral phase and to provide an early warning alert able to drive the VHE search. We compute the amount of possible joint detections by considering several observational strategies, and demonstrate that the sensitivity of CTA make the detection of the VHE emission possible even if it is several orders fainter than that observed at 10 keV–10 MeV. We discuss the results in terms of possible scenarios of the production of VHE photons from binary neutron star mergers.

Measuring the Hubble Constant with Dark Neutron Star–Black Hole Mergers

Detection of gravitational waves (GWs) from neutron star-black hole (NSBH) standard sirens provides local measurements of the Hubble constant (H 0), regardless of the detection of an electromagnetic (EM) counterpart, given that matter effects can be exploited to break the redshift degeneracy of the GW waveforms. The distinctive merger morphology and the high-redshift detectability of tidally disrupted NSBH make them promising candidates for this method. Also, the detection prospects of an EM counterpart for these systems will be limited to z < 0.8 in the optical, in the era of future GW detectors. Using recent constraints on the equation of state of NSs from multi-messenger observations of NICER and LIGO/Virgo/KAGRA, we show the prospects of measuring H 0 solely from GW observation of NSBH systems, achievable by the Einstein telescope (ET) and Cosmic Explorer (CE) detectors. We first analyze individual events to quantify the effect of high-frequency (≥500 Hz) tidal distortions on the inference of NS tidal deformability parameter (Λ) and hence on H 0. We find that disruptive mergers can constrain Λ up to more precisely than nondisruptive ones. However, this precision is not sufficient to place stringent constraints on the H 0 from individual events. By performing Bayesian analysis on simulated NSBH data (up to N = 100 events, corresponding to a day of observation) in the ET+CE detectors, we find that NSBH systems enable unbiased 4%–13% precision on the estimate of H 0 (68% credible interval). This is a similar measurement precision found in studies analyzing NSBH mergers with EM counterparts in the LVKC O5 era.

Revise Thermal Winds of Remnant Neutron Stars in Gamma-Ray Bursts

It seems that the wealth of information revealed by the multi-messenger observations of the binary neutron star (NS) merger event, GW170817/GRB 170817A/kilonova AT2017gfo, places irreconcilable constraints to models of the prompt emission of this gamma-ray burst (GRB). The observed time delay between the merger of the two NSs and the trigger of the GRB and the thermal tail of the prompt emission can hardly be reproduced by these models simultaneously. We argue that the merger remnant should be an NS (last for, at least, a large fraction of 1 s), and that the difficulty can be alleviated by the delayed formation of the accretion disk due to the absorption of high-energy neutrinos emitted by the NS and the delayed emergence of effective viscosity in the disk. Further, we extend the consideration of the effect of the energy deposition of neutrinos emitted from the NS. If the NS is the central object of a GRB with a distance and duration similar to that of GRB 170817A, thermal emission of the thermal bubble inflated by the NS after the termination of accretion may be detectable. If our scenario is verified, it would be of interest to investigate the cooling of nascent NSs.

Baryonic models of ultra-low-mass compact stars for the central compact object in HESS J1731-347

The recent attempt on mass and radius inference of the central compact object within the supernova remnant HESS J1731-347 suggests for this object an unusually low mass of M=0.77−0.17+0.20M⊙ and a small radius of R=10.4−0.78+0.86km. We explore the ways such a result can be accommodated within models of dense matter with heavy baryonic degrees of freedom which are constrained by the multi-messenger observations. We find that to do so using only purely nucleonic models, one needs to assume a rather small value of the slope of symmetry energy Lsym. Once heavy baryons are included higher values of the slope Lsym become acceptable at the cost of a slightly reduced maximum mass of static configuration. These two scenarios are distinguished by the particle composition and will undergo different cooling scenarios. In addition, we show that the universalities of the I-Love-Q relations for static configurations can be extended to very low masses without loss in their accuracy.

Higgs Field-Induced Triboluminescence in Binary Black Hole Mergers

We conjecture that the Higgs potential can be significantly modified when it is in close proximity to the horizon of an astrophysical black hole, leading to the destabilization of the electroweak vacuum. In this situation, the black hole should be encompassed by a shell consisting of a “bowling substance” of the nucleating new-phase bubbles. In a binary black-hole merger, just before the coalescence, the nucleated bubbles can be prevented from falling under their seeding horizons, as they are simultaneously attracted by the gravitational potential of the companion. For a short time, the unstable vacuum will be “sandwiched” between two horizons of the binary black hole, and therefore the bubbles may collide and form micro-black holes, which are rapidly evaporated by thermal emission of Hawking radiation of all Standard Model species. This evaporation, being triggered by a gravitational wave signal from the binary black-hole merger, can manifest itself in observations of gamma rays and very-high-energy neutrinos, which makes it a perfect physics case for multi-messenger astronomical observations.

Simulating neutron stars with a flexible enthalpy-based equation of state parametrization in spectre

Numerical simulations of neutron star mergers represent an essential step toward interpreting the full complexity of multimessenger observations and constraining the properties of supranuclear matter. Currently, simulations are limited by an array of factors, including computational performance and input physics uncertainties, such as the neutron star equation of state. In this work, we expand the range of nuclear phenomenology efficiently available to simulations by introducing a new analytic parametrization of cold, beta-equilibrated matter that is based on the relativistic enthalpy. We show that the new enthalpy parametrization can capture a range of nuclear behavior, including strong phase transitions. We implement the enthalpy parametrization in the spectre code, simulate isolated neutron stars, and compare performance to the commonly used spectral and polytropic parametrizations. We find comparable computational performance for nuclear models that are well represented by either parametrization, such as simple hadronic equations of state. We show that the enthalpy parametrization further allows us to simulate more complicated hadronic models or models with phase transitions that are inaccessible to current parametrizations.