Neutron Star Mergers Research Articles

Can we identify primordial black holes? Tidal tests for subsolar-mass gravitational-wave observations

The detection of a subsolar object in a compact binary merger is regarded as one of the smoking gun signatures of a population of primordial black holes (PBHs). We critically assess whether these systems could be distinguished from stellar binaries, for example composed of white dwarfs or neutron stars, which could also populate the subsolar mass range. At variance with PBHs, the gravitational-wave signal from stellar binaries is affected by tidal effects, which dramatically grow for moderately compact stars as those expected in the subsolar range. We forecast the capability of constraining tidal effects of putative subsolar neutron star binaries with current and future LIGO-Virgo-KAGRA (LVK) sensitivities as well as next-generation experiments. We show that, should LVK O4 run observe subsolar neutron-star mergers, it could measure the (large) tidal effects with high significance. In particular, for subsolar neutron-star binaries, O4 and O5 projected sensitivities would allow measuring the effect of tidal disruption on the waveform in a large portion of the parameter space, also constraining the tidal deformability at O(10%) level, thus excluding a primordial origin of the binary. Vice versa, for subsolar PBH binaries, model-agnostic upper bounds on the tidal deformability can rule out neutron stars or more exotic competitors. Assuming events similar to the subthreshold candidate SSM200308 reported in LVK O3b data are PBH binaries, O4 projected sensitivity would allow ruling out the presence of neutron-star tidal effects at ≈3σ CL, thus strengthening the PBH hypothesis. Future experiments would lead to even stronger (>5σ) conclusions on potential discoveries of this kind. Published by the American Physical Society 2024

Comparative study of kilonova opacities for three elements of the sixth period (hafnium, osmium, and gold) from new atomic structure calculations in Hf I–IV, Os I–IV, and Au I–IV

Aims. It is now well established that a large amount of heavy (trans-iron) elements are produced during neutron star (NS) mergers. These elements can be detected in the spectra of the kilonova emitted from the post-merger ejected materials. Due to the high level densities that characterize the complex configurations belonging to heavy elements, thus giving rise to millions of absorption lines, the kilonova ejecta opacity is of significant importance. The elements that contribute the most to the latter are those with an unfilled nd subshell belonging to the fifth and the sixth rows of the periodic table, and those with an unfilled nf subshell belonging to the lanthanide and actinide groups. The aim of the present work is to make a new contribution to this field by performing large-scale atomic structure calculations in three specific sixth-row 5d elements, namely hafnium, osmium, and gold, in the first four charge stages (I–IV), and by computing the corresponding opacities, while focusing on the importance of the atomic models used. Methods. The pseudo-relativistic Hartree–Fock (HFR) method, including extended sets of interacting configurations, was used for the atomic structure and radiative parameter calculations, while the expansion formalism was used to estimate the opacities. Results. Theoretical energy levels, wavelengths, and oscillator strengths were computed for millions of spectral lines in Hf I–IV, Os I–IV, and Au I–IV ions, the reliability of these parameters being assessed through detailed comparisons with previously published experimental and theoretical results. The newly obtained atomic data were then used to calculate expansion opacities for typical kilonova conditions expected one day after the NS merger; these are a density of ρ = 10−13 g cm−3 and temperatures ranging from T = 5000 K to T = 15 000 K. Some agreements and differences were found when comparing our results with available data, highlighting the importance of using sufficiently complete atomic models for the determination of opacities.

Inference of neutron-star properties with unified crust-core equations of state for parameter estimation

Context. Relating different global neutron-star (NS) properties, such as tidal deformability and radius, or mass and radius, requires an equation of state (EoS). Determining the NS EoS is therefore not only the science goal of a variety of observational projects, but it also enters in the analysis process; for example, to predict a NS radius from a measured tidal deformability via gravitational waves (GW) during the inspiral of a binary NS merger. To this aim, it is important to estimate the theoretical uncertainties on the EoS, one of which is the possible bias coming from an inconsistent treatment of the low-density region; that is, the use of a so called non-unified NS crust. Aims. We propose a numerical tool allowing the user to consistently match a nuclear-physics informed crust to an arbitrary high-density EoS describing the core of the star. Methods. We introduce an inversion procedure of the EoS close to saturation density that allows users to extract nuclear-matter parameters and extend the EoS to lower densities in a consistent way. For the treatment of inhomogeneous matter in the crust, a standard approach based on the compressible liquid-drop (CLD) model approach was used in our work. A Bayesian analysis using a parametric agnostic EoS representation in the high-density region is also presented in order to quantify the uncertainties induced by an inconsistent treatment of the crust. Results. We show that the use of a fixed, realistic-but-inconsistent model for the crust causes small but avoidable errors in the estimation of global NS properties and leads to an underestimation of the uncertainties in the inference of NS properties. Conclusions. Our results highlight the importance of employing a consistent EoS in inference schemes. The numerical tool that we developed to reconstruct such a thermodynamically consistent EoS, CUTER, has been tested and validated for use by the astrophysical community.

Realistic models of general-relativistic differentially rotating stars

ABSTRACT General-relativistic equilibria of differentially rotating stars are expected in a number of astrophysical scenarios, from core-collapse supernovae to the remnant of binary neutron-star mergers. The latter, in particular, have been the subject of extensive studies where they were modelled with a variety of laws of differential rotation with varying degree of realism. Starting from accurate and fully general-relativistic simulations of binary neutron-star mergers with various equations of state and mass ratios, we establish the time when the merger remnant has reached a quasi-stationary equilibrium and extract in this way realistic profiles of differential rotation. This allows us to explore how well traditional laws reproduce such differential-rotation properties and to derive new laws of differential rotation that better match the numerical data in the low-density Keplerian regions of the remnant. In this way, we have obtained a novel and somewhat surprising result: the dynamical stability line to quasi-radial oscillations computed from the turning-point criterion can have a slope that is not necessarily negative with respect to the central rest-mass density, as previously found with traditional differential-rotation laws. Indeed, for stellar models reproducing well the properties of the merger remnants, the slope is actually positive, thus reflecting remnants with angular momentum at large distances from the rotation axis, and hence with cores having higher central rest-mass densities and slower rotation rates.

Location of the neutron drip line for Sn and its impact on r-process abundances

Nuclear physical inputs are critical for predicting the abundance of r-process elements. Extensive sensitivity studies have recently been performed to gauge the impact of the individual properties of nuclei on the r-process. In this work, we investigate the impact of the large uncertainties in the theoretical predictions of the masses of neutron-rich Sn isotopes and the location of the neutron drip line on the abundance of r-process abundances. The uncertainties in the predicted r-process abundances are obtained through large-scale network calculations by simultaneously varying the masses and reaction rates of Sn within the predicted mass uncertainties. The calculations use a generally indicative astrophysical trajectory of neutron-star mergers and are based on three different mass models. The results indicate that the large uncertainty in the location of the neutron drip line for Sn can significantly and asymmetrically affect the predicted abundances of nuclei after the second peak, in particular around A=140−155, sometimes accompanied by a delayed freeze-out.

Prospects for weighing neutrinos in interacting dark energy models using joint observations of gravitational waves and γ-ray bursts* *Supported by the National Natural Science Foundation of China (12305069, 11947022, 11975072, 11875102, 11835009), the National SKA Program of China (2022SKA0110200, 2022SKA0110203), the National 111 Project (B16009), and the Program of the Education Department of Liaoning Province, China (JYTMS20231695)

Cosmological observations can be used to weigh neutrinos, but this method is model-dependent, with results relying on the cosmological model considered. If we consider interactions between dark energy and dark matter, the neutrino mass constraints differ from those derived under the standard model. On the contrary, gravitational wave (GW) standard siren observations can measure absolute cosmological distances, helping to break parameter degeneracies inherent in traditional cosmological observations, thereby improving constraints on neutrino mass. This paper examines the constraints on neutrino mass within interacting dark energy (IDE) models and explores how future GW standard siren observations could enhance these results. For multi-messenger GW observations, we consider the joint observations of binary neutron star mergers by third-generation ground-based GW detectors and short γ-ray burst observations by missions similar to the THESEUS satellite project. Using current cosmological observations (CMB+BAO+SN), we obtain an upper limit on the neutrino mass in the IDE models of 0.15 (or 0.16) eV. With the inclusion of GW data, the upper limit on the neutrino mass improves to 0.14 eV. This indicates that in the context of IDE models, the improvement in neutrino mass constraints from GW observations is relatively limited. However, GW observations significantly enhance the constraints on other cosmological parameters, such as matter density parameter, the Hubble constant, and coupling strength between dark energy and dark matter.

Massive propagating modes of torsion

The dynamics of the torsion field is analyzed in the framework of the Covariant Canonical Gauge Theory of Gravity (CCGG), a De Donder–Weyl Hamiltonian formulation of gauge gravity. The action is quadratic in both, the torsion and the Riemann–Cartan tensor. Since the latter adds the derivative of torsion to the equations of motion, torsion is no longer identical to spin density, as in the Einstein–Cartan theory, but an additional propagating degree of freedom. As torsion turns out to be totally anti-symmetric, it can be parametrised via a single axial vector. It is shown in this paper that, in the weak torsion limit, the axial vector obeys a wave equation with an effective mass term which is partially dependent on the scalar curvature. The source of torsion is thereby given by the fermion axial current which is the net fermionic spin density of the system. Possible measurable effects and approaches to experimental analysis are addressed. For example, neutron star mergers could act as a dipoles or quadrupoles for torsional radiation, and an analysis of radiation of pulsars could lead to a detection of torsion wave background radiation.

Probing ultralight isospin-violating mediators at GW170817

Gravitational wave (GW) signals arising from binary neutron star mergers offer new, sensitive probes to ultralight mediators. Here we analyze the GW signals in the GW170817 event detected by the LIGO/Virgo collaboration to impose constraints on the ultralight isospin-violating mediator that has different couplings to protons and neutrons. Neutron stars, which primarily consist of neutrons, are the ideal places to probe the isospin-violating mediator. Such a mediator can significantly alter the dynamics of the binary neutron star mergers, through both the long-range Yukawa force and the new dipole radiation. We compute the gravitational waveform by taking into account the new physics effects due to the isospin-violating mediator and use the Bayesian inference to analyze the gravitational wave data in the GW170817 event. We find that although the current fifth force experiments (including MICROSCOPE and EW) often provide more stringent constraints than the GW170817 data, in the parameter space where the isospin-violating force is completely screened by the Earth (namely, the Earth is charge neutral under this force), the GW170817 data offer the leading constraints: the upper bound on the neutron coupling is fn ≲ 10−19 in the mediator mass range of ≃ (3 × 10−16, 5 × 10−14) eV.

Extension of semi-analytic kilonova model and the application in multicolour light curves fitting

ABSTRACT Observations of AT2017gfo, the counterpart of the double neutron star merger, confirmed the existence of kilonovae that are suggested to originate from the decay of the very heavy elements. However, the different simplified assumptions used in analytical and semi-analytical models can lead to a variety of fitting results. Previous semi-analytical models have assumed a grey body because of the challenge in calculating the detailed opacity of ejecta and to save the calculation time. However, the variation of opacity with wavelength cannot be ignored when calculating the multicolour flux. Though numerical simulations can calculate kilonova details, they require an excessive amount of time for data modelling. In this work, we propose an extension of previous semi-analytical models by not assuming the grey-body opacity and constant thermalization efficiency assumptions. We demonstrate that the impact of opacity changing with wavelength is non-negligible and then introduce an empirical analytical opacity in optical and near-infrared bands to calculate multicolour light curves of kilonovae with a three-dimensional, three-component model. Additionally, we consider the evolution of thermalization efficiency at later times, enabling application of our model to both the early and later phases of kilonovae. To enable efficient fitting of models with many parameters (23 parameters in the case of double neutron star merger), a quick method based on the greedy algorithm is developed. We can well reproduce the infrared/optical/ultraviolet observations of AT2017gfo. This model can also be extended to black-hole and neutron star mergers as well.

Finite-temperature equations of state of compact stars with hyperons: three-dimensional tables

We construct tables of finite temperature equation of state (EoS) of hypernuclear matter in the range of densities, temperatures, and electron fractions that are needed for numerical simulations of supernovas, proto-neutron stars, and binary neutron star mergers and cast them in the format of CompOSE database. The tables are extracted from a model that is based on covariant density functional (CDF) theory that includes the full JP=1/2+\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$J^P=1/2^+$$\\end{document} baryon octet in a manner that is consistent with the current astrophysical and nuclear constraints. We employ a parameterization with three different values of the slope of the symmetry energy Lsym=30\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$L_\ extrm{sym}=30$$\\end{document}, 50 and 70 MeV and fixed skewness Qsat=400\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$Q_\ extrm{sat}= 400$$\\end{document} MeV for above saturation matter. A model for the EoS of inhomogeneous matter is matched at sub-saturation density to the high-density hypernuclear EoS. We discuss the generic features of the resulting EoS and the composition of matter as a function of density, temperature, and electron fraction. The nuclear characteristics and strangeness fraction of these models are compared to the alternatives from the literature. The integral properties of static and rapidly rotating compact stars in the limit of zero temperature are discussed and confronted with the multimessenger astrophysical constraints.

Research on Binary Black Hole Systems by Analysis of Gravitational Waves

The mass ratio, effective spin, and chirp mass are three important properties to describe the binary systems of black holes. The gravitational wave sources can be a breakthrough in the investigation of black holes. The formation of black holes is a significant problem to solve as the electromagnetic waves (light) cannot escape from the black hole, and the gravitational waves can be emitted when the binary systems of black holes or neutron stars merge. Based on the detection of gravitational waves, the process of the formation of black holes can be traced and studied in the binary black holes and primordial black holes formed in the early universe. The mass ratio distribution is plotted to investigate the difference between the binary black hole and neutron star black hole and the mass ratio which makes it easier for the merger of the black holes to happen. The investigation of the effective spin helps classify the model of the formation of black holes as there are five models of the black hole. The statistical model can be tested with the chosen data, and the cumulative standard Gaussian distribution is used in this paper. The result states that the model is not suitable for the chosen data and another statistical model should be introduced to analyse the effective spin. These parameters are all important for determining the stage and formation process of the binary systems.

HERMES: Gamma-ray burst and gravitational wave counterpart hunter

Gamma-ray bursts (GRBs) bridge relativistic astrophysics and multimessenger astronomy. Space--based gamma - and X-ray wide-field detectors have proven essential for detecting and localizing the highly variable GRB prompt emission, which is also a counterpart of gravitational wave events. We studied the capability of detecting long and short GRBs with the High Energy Rapid Modular Ensemble of Satellites (HERMES) Pathfinder (HP) and SpIRIT, namely a swarm of six 3U CubeSats to be launched in early 2025, and a 6U CubeSat launched on December 1 2023. We also studied the capabilities of two advanced configurations of swarms of more than eight satellites with improved detector performances (HERMES Constellations). The HERMES detectors, sensitive down to sim \,2-3\,keV, will be able to detect faint and soft GRBs, which comprise X-ray flashes and high-redshift bursts. By combining state-of-the-art long- and short-GRB population models with a description of the single module performance, we estimate that HP will detect sim $ long GRBs ($3.4^ $ short GRBs per year. The larger HERMES Constellations under study can detect between sim \,1300 and sim \,3000 long GRBs per year and between sim \,160 and sim \,400 short GRBs per year, depending on the chosen configuration, with a rate of long GRBs above $z>6$ of between 30 and 75 per year. Finally, we explored the capability of HERMES to detect short GRBs as electromagnetic counterparts of binary neutron star (BNS) mergers detected as gravitational signals by current and future ground--based interferometers. Under the assumption that the GRB jets are structured, we estimate that HP can provide up to sim 1 (14) $ joint detections during the fifth LIGO-Virgo-KAGRA observing run ( Einstein Telescope single triangle 10 km arm configuration). These numbers become sim \,4 (100) $, respectively, for the HERMES Constellation configuration.

Revealing phase transition in dense matter with gravitational wave spectroscopy of binary neutron star mergers

Revealing phase transition in dense matter with gravitational wave spectroscopy of binary neutron star mergers

Detecting short-term gravitational waves from post-merger hyper-massive neutron stars with a kilohertz detector

Abstract Gravitational waves emanating from binary neutron star inspirals, alongside electromagnetic transients resulting from the aftermath of the GW170817 merger, have been successfully detected. However, the intricate post-merger dynamics that bridge these two sets of observables remain enigmatic. This includes if, and when, the post-merger remnant star collapses to a black hole, and what are the necessary conditions to power a short gamma-ray burst, and other observed electromagnetic counterparts. Our focus is on the detection of gravitational wave (GW) emissions from hyper-massive neutron stars (NSs) formed through binary neutron star (BNS) mergers. Utilizing several kilohertz GW detectors, we simulate BNS mergers within the detection limits of LIGO-Virgo-KARGA O4. Our objective is to ascertain the fraction of simulated sources that may emit detectable post-merger GW signals. For kilohertz detectors equipped with a new cavity design, we estimate that approximately 1.1%–32% of sources would emit a detectable post-merger GW signal. This fraction is contingent on the mass converted into gravitational wave energy, ranging from 0.01M sun to 0.1M sun. Furthermore, by evaluating other well-regarded proposed kilohertz GW detectors, we anticipate that the fraction can increase to as much as 2.1%–61% under optimal performance conditions.

Multi-messenger prospects for black hole – neutron star mergers in the O4 and O5 runs

The existence of merging black hole-neutron star (BHNS) binaries has been ascertained through the observation of their gravitational wave (GW) signals. However, to date, no definitive electromagnetic (EM) emission has been confidently associated with these mergers. Such an association could help unravel crucial information on these systems, for example, their BH spin distribution, the equation of state (EoS) of the neutron star and the rate of heavy element production. We modeled the multi-messenger (MM) emission from BHNS mergers detectable during the fourth (O4) and fifth (O5) observing runs of the LIGO-Virgo-KAGRA (LVK) GW detector network in order to provide detailed predictions that can help enhance the effectiveness of observational efforts and extract the highest possible scientific information from such remarkable events. Our methodology is based on a population synthesis approach, which includes the modeling of the signal-to-noise ratio of the GW signal in the detectors, the GW-inferred sky localization of the source, the kilonova (KN) optical and near-infrared light curves, the relativistic jet gamma-ray burst (GRB) prompt emission peak photon flux, and the GRB afterglow light curves in the radio, optical, and X-ray bands. The resulting prospects for BHNS MM detections during O4 are not promising, with an LVK GW detection rate of 15.0−8.8+15.4 yr−1, but joint MM rates of ∼10−1 yr−1 for the KN and ∼10−2 yr−1 for the jet-related emission. In O5, we found an overall increase in expected detection rates by around an order of magnitude, owing to both the enhanced sensitivity of the GW detector network and the coming online of future EM facilities. Considering variations in the NS EoS and BH spin distribution, we find that the detection rates can increase further by up to a factor of several tens. Finally, we discuss direct searches for the GRB radio afterglow with large field-of-view instruments during O5 and beyond as a new possible follow-up strategy in the context of ever-dimming prospects for KN detection due to the recession of the GW horizon.

Diversity of Early Kilonova with the Realistic Opacities of Highly Ionized Heavy Elements

We investigate early (t < 1 day) kilonova from a neutron star merger by deriving atomic opacities for all the elements from La to Ra (Z = 57–88) ionized to the states V–XI. The opacities at high temperatures for the elements with open f-shells (e.g., lanthanides) are exceptionally high, reaching κexp∼3×103cm2g−1 at λ ≤ 1000 Å at T ∼ 70,000 K, whereas the opacities at the same temperature and wavelengths for the elements with open d-, p-, and s-shells reach κexp∼1 , 0.1, and 0.01 cm2 g−1, respectively. Using the new opacity data set, we derive early kilonovae for various compositions and density structures expected for neutron star merger ejecta. The bolometric luminosity of the lanthanide-rich ejecta shows distinct signatures and is fainter than that of the lanthanide-free ejecta. Early luminosity is suppressed by the presence of a thin outer layer, agreeing with the results of Kasen et al. and Banerjee et al. The early brightness in the Swift UVOT filters and in the optical g, r, i, and z filters for a source at 100 Mpc are about ∼22–19.5 and ∼21–20 mag, respectively, at t ∼ 0.1 day. Such kilonovae are ideal targets for the upcoming UV satellites, such as ULTRASAT, UVEX, and DORADO, and the upcoming surveys, e.g., the Vera Rubin Observatory. We suggest that the gray opacities that reproduce the bolometric light curves with and without lanthanides are ∼1–10 and ∼0.8 cm2 g−1.

Revisiting the Constraint on the Equation of State of Neutron Stars Based on Binary Neutron Star Mergers

The merger of a neutron star (NS)–NS binary can form different productions of compact remnants, among which a supramassive NS (SMNS) can create an internal plateau, and the following steep decay marks the collapse of the SMNS. The proportion of the SMNS and the corresponding collapse time are often used to constrain the NS equation of state (EOS). This paper revisits this topic by considering the effect of an accretion disk on a compact remnant, which is not considered in previous works. Compared with previous works, the collapse-time distribution (peaks ∼100 s) of SMNSs formed from an NS–NS merger is almost unaffected by the initial surface magnetic field (B s,i ) of the NS, but the total energy output of the magnetic dipole radiation from the SMNSs depends on B s,i significantly. Coupling the constraints from the SMNS fraction, we exclude some EOSs and obtain three candidate EOSs, i.e., DD2, ENG, and MPA1. By comparing the distributions of the collapse time and the luminosity of the internal plateau (in the short gamma-ray bursts) for observations obtained based on the three candidate EOSs, it is shown that only the EOS of ENG is favored. Our sample, based on the ENG EOS and a mass distribution motivated by Galactic systems, suggests that approximately 99% of NS–NS mergers collapse to form a black hole (BH) within 107s. This includes scenarios promptly forming a BH (36.5%), an SMNS (60.7%), or a stable NS that transitions into a BH or an SMNS following accretion (2.1%). It also indicates that the remnants for GW170817 and GW190425, and the second object of GW190814, are more likely to be BHs.

Current Sheet Alignment in Oblique Black Hole Magnetospheres: A Black Hole Pulsar?

We study the magnetospheric evolution of a nonaccreting spinning black hole (BH) with an initially inclined split monopole magnetic field by means of 3D general relativistic magnetohydrodynamic simulations. This serves as a model for a neutron star (NS) collapse or a BH–NS merger remnant after the inherited magnetosphere has settled into a split monopole field creating a striped wind. We show that the initially inclined split monopolar current sheet aligns over time with the BH equatorial plane. The inclination angle evolves exponentially toward alignment, with an alignment timescale that is inversely proportional to the square of the BH angular velocity, where higher spin results in faster alignment. Furthermore, magnetic reconnection in the current sheet leads to exponential decay of event-horizon-penetrating magnetic flux with nearly the same timescale for all considered BH spins. In addition, we present relations for the BH mass and spin in terms of the period and alignment timescale of the striped wind. The explored scenario of a rotating, aligning, and reconnecting current sheet can potentially lead to multimessenger electromagnetic counterparts to a gravitational-wave event due to the acceleration of particles powering high-energy radiation, plasmoid mergers resulting in coherent radio signals, and pulsating emission due to the initial misalignment of the BH magnetosphere.

Determining the Hubble constant with AGN-assisted black hole mergers

ABSTRACT Gravitational waves from neutron star mergers have long been considered a promising way to measure the Hubble constant, H0, which describes the local expansion rate of the Universe. While black hole mergers are more abundantly observed, their expected lack of electromagnetic emission and poor gravitational-wave localization make them less well suited for measuring H0. Black hole mergers within the discs of Active Galactic Nuclei (AGN) could be an exception. Accretion from the AGN disc may produce an electromagnetic signal, pointing observers to the host galaxy. Alternatively, the low number density of AGNs could help identify the host galaxy of $1{\!-\!}5~{{\ \rm per\ cent}}$ of mergers. Here we show that black hole mergers in AGN discs may be a sensitive way to determine H0 with gravitational waves. If 1 per cent (10 per cent) of LIGO’s observations occur in AGN discs with identified host galaxies, we could measure H0 with 12 per cent (4 per cent) uncertainty in five years, possibly comparable to the sensitivity of neutron star mergers and set to considerably improve current gravitational wave measurements.

Chaos in inhomogeneous neutrino fast flavor instability

In dense neutrino gases, the neutrino-neutrino coherent forward scattering gives rise to a complex flavor oscillation phenomenon not fully incorporated in simulations of neutron star mergers (NSM) and core collapse supernovae (CCSNe). Moreover, it has been proposed to be chaotic, potentially limiting our ability to predict neutrino flavor transformations in simulations. To address this issue, we explore how small flavor perturbations evolve in the nonlinear regime of the neutrino quantum kinetic equation within a narrow centimeter-scale region inside a NSM and a toy neutrino distribution. Our findings reveal that paths in the flavor state space of solutions with similar initial conditions diverge exponentially, exhibiting chaos. This inherent chaos makes the microscopic scales of neutrino flavor transformations unpredictable. However, the domain-averaged neutrino density matrix remains relatively stable, with chaos minimally affecting it. This particular property suggests that domain-averaged quantities remain reliable despite the exponential amplification of errors. Published by the American Physical Society 2024