Neutron Star Mergers Research Articles

Nuclear physics constraints from binary neutron star mergers in the Einstein Telescope era

Nuclear physics constraints from binary neutron star mergers in the Einstein Telescope era

Effects of onset of phase transition on binary neutron star mergers

ABSTRACT Quantum Chromodynamics predicts phase transition from hadronic matter to quark matter at high density, which is highly probable in astrophysical systems like binary neutron star mergers. To explore the critical density where such phase transition can occur, we performed numerical relativity simulations of binary neutron star mergers with various masses (equal and unequal binaries). We aim to understand the effect of the onset of phase transition on the merger dynamics and gravitational wave spectra. We generated a set of equations of states by agnostically changing the onset of phase transition, having the hadronic matter part and quark matter part fixed. This particular arrangement of the equation of states explores the scenario of mergers where mixed phases of matter are achieved before or during the merger. Under these circumstances, if the matter properties with hadronic and quark degrees differ significantly, it is reflected in the stability of the final merger product for the intermediate mass binary. We performed a case study on mixed species merger, where one of the binary companions is hybrid star. If quark matter appears at low densities, we observe significant change in post-merger gravitational wave analysis in terms of higher peak frequencies and post-merger frequencies in power spectral density. We report indications expressed as spikes in phase difference plots at merger time for mixed mergers. We found that the expression of phase transition in post-merger gravitational wave signals is more significant for unequal mass binary than for equal mass binary having the same total baryonic mass.

Element abundance patterns in stars indicate fission of nuclei heavier than uranium

The heaviest chemical elements are naturally produced by the rapid neutron-capture process ( r -process) during neutron star mergers or supernovae. The r -process production of elements heavier than uranium (transuranic nuclei) is poorly understood and inaccessible to experiments so must be extrapolated by using nucleosynthesis models. We examined element abundances in a sample of stars that are enhanced in r -process elements. The abundances of elements ruthenium, rhodium, palladium, and silver (atomic numbers Z = 44 to 47; mass numbers A = 99 to 110) correlate with those of heavier elements (63 ≤ Z ≤ 78, A > 150). There is no correlation for neighboring elements (34 ≤ Z ≤ 42 and 48 ≤ Z ≤ 62). We interpret this as evidence that fission fragments of transuranic nuclei contribute to the abundances. Our results indicate that neutron-rich nuclei with mass numbers >260 are produced in r -process events.

Jets from Neutron-Star Merger Remnants and Massive Blue Kilonovae.

We perform three-dimensional general-relativistic magnetohydrodynamic simulations with weak interactions of binary neutron-star (BNS) mergers resulting in a long-lived remnant neutron star, with properties typical of galactic BNS and consistent with those inferred for the first observed BNS merger GW170817. We demonstrate self-consistently that within ≲30 ms postmerger magnetized (σ∼5-10) incipient jets emerge with asymptotic Lorentz factor Γ∼5-10, which successfully break out from the merger debris within ≲20 ms. A fast (v≲0.6c), magnetized (σ∼0.1) wind surrounds the jet core and generates a UV/blue kilonova precursor on timescales of hours, similar to the precursor signal due to free neutron decay in fast dynamical ejecta. Postmerger ejecta are quickly dominated by magnetohydrodynamically driven outflows from an accretion disk. We demonstrate that, within only 50ms postmerger, ≳2×10^{-2}M_{⊙} of lanthanide-free, quasispherical ejecta with velocities ∼0.1-0.2c is launched, yielding a kilonova signal consistent with GW170817 on timescales of ≲5 d.

Impact of a mean field dynamo on neutron star mergers leading to magnetar remnants

Impact of a mean field dynamo on neutron star mergers leading to magnetar remnants

Big-bang Nucleosynthesis on a bubble universe nucleated in Kerr-AdS5.

Modelling the Galactic chemical evolution (GCE) of short-lived radioisotopes (SLRs, with half-lives of the order of million years) can provide timing information on recent nucleosynthesis. The knowledge of their spatial distribution throughout the interstellar medium (ISM) is crucial. We are using a three-dimensional GCE model to investigate the evolution of four SLRs: \(^{53}\)Mn from supernovae of type Ia (SNeIa), \(^{60}\)Fe from core-collapse supernovae (CCSNe), \(^{182}\)Hf from stellar winds from intermediate mass stars (IMSs), and \(^{244}\)Pu from neutron star mergers (NSMs) to explain the recent (within the last $\approx 1 -- 20$ Myr) deposition of live SLRs into deep-sea reservoirs. We find that although these SLRs are produced at very different nucleosynthetic sites, they are deposited on Earth conjointly.

Differential rotation in compact objects with hyperons and delta isobars

AbstractNeutron stars may experience differential rotation on short, dynamical timescales following extreme astrophysical events like binary neutron star mergers. In this work, the masses and radii of differentially rotating neutron star models are computed. We employ a set of equations of states for dense hypernuclear and ‐admixed‐hypernuclear matter obtained within the framework of CDF theory in the relativistic Hartree‐Fock (RHF) approximation. Results are shown for varying meson‐ couplings, or equivalently the ‐potential in nuclear matter. A comparison of our results with those obtained for nonrotating stars shows that the maximum mass difference between differentially rotating and static stars is independent of the underlying particle composition of the star. We further find that the decrease in the radii and increase in the maximum masses of stellar models when ‐isobars are added to hyperonuclear matter (as initially observed for static and uniformly rotating stars) persist also in the case of differentially rotating neutron stars.

SPLUS J142445.34–254247.1: An r-process–enhanced, Actinide-boost, Extremely Metal-poor Star Observed with GHOST

We report on a chemo-dynamical analysis of SPLUS J142445.34−254247.1 (SPLUS J1424−2542), an extremely metal-poor halo star enhanced in elements formed by the rapid neutron-capture process (r-process). This star was first selected as a metal-poor candidate from its narrowband S-PLUS photometry and followed up spectroscopically in medium resolution with Gemini-South/GMOS, which confirmed its low-metallicity status. High-resolution spectroscopy was gathered with GHOST at Gemini-South, allowing for the determination of the chemical abundances for 36 elements, from carbon to thorium. At [Fe/H] = −3.39, SPLUS J1424−2542 is one of the lowest-metallicity stars with measured Th and has the highest observed to date, making it part of the “actinide-boost” category of r-process–enhanced stars. The analysis presented here suggests that the gas cloud from which SPLUS J1424−2542 formed must have been enriched by at least two progenitor populations. The light-element (Z ≤ 30) abundance pattern is consistent with the yields from a supernova explosion of metal-free stars with 11.3–13.4 M ⊙, and the heavy-element (Z ≥ 38) abundance pattern can be reproduced by the yields from a neutron star merger (1.66 M ⊙ and 1.27 M ⊙) event. A kinematical analysis also reveals that SPLUS J1424−2542 is a low-mass, old halo star with a likely in situ origin, not associated with any known early merger events in the Milky Way.

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.

Predictions for electromagnetic counterparts to Neutron Star mergers discovered during LIGO-Virgo-KAGRA observing runs 4 and 5

ABSTRACT We present a comprehensive, configurable open-source software framework for estimating the rate of electromagnetic detection of kilonovae (KNe) associated with gravitational wave detections of binary neutron star (BNS) mergers. We simulate the current LIGO-Virgo-KAGRA (LVK) observing run (O4) using current sensitivity and uptime values as well as using predicted sensitivites for the next observing run (O5). We find the number of discoverable kilonovae during LVK O4 to be ${ 1}_{- 1}^{+ 4}$ or ${ 2 }_{- 2 }^{+ 3 }$, (at 90 per cent confidence) depending on the distribution of NS masses in coalescing binaries, with the number increasing by an order of magnitude during O5 to ${ 19 }_{- 11 }^{+ 24 }$. Regardless of mass model, we predict at most five detectable KNe (at 95 per cent confidence) in O4. We also produce optical and near-infrared light curves that correspond to the physical properties of each merging system. We have collated important information for allocating observing resources for search and follow-up observations, including distributions of peak magnitudes in several broad-bands and time-scales for which specific facilities can detect each KN. The framework is easily adaptable, and new simulations can quickly be produced in response to updated information such as refined merger rates and NS mass distributions. Finally, we compare our suite of simulations to the thus-far completed portion of O4 (as of 2023, October 14), finding a median number of discoverable KNe of 0 and a 95 percentile upper limit of 2, consistent with no detections so far in O4.

The short gamma-ray burst population in a quasi-universal jet scenario

We present a model of the short gamma-ray burst (SGRB) population under a ‘quasi-universal jet’ scenario in which jets can differ somewhat in their on-axis peak prompt emission luminosity, Lc, but share a universal angular luminosity profile, ℓ(θv) = L(θv)/Lc, as a function of the viewing angle, θv. The model was fitted, through a Bayesian hierarchical approach inspired by gravitational wave (GW) population analyses, to three observed SGRB samples simultaneously: the Fermi/GBM sample of SGRBs with spectral information available in the catalogue (367 events); a flux-complete sample of 16 Swift/BAT SGRBs that are also detected by the GBM and have a measured redshift; and a sample of SGRBs with a binary neutron star (BNS) merger counterpart, which only includes GRB 170817A at present. Particular care was put into modelling selection effects. The resulting model, which reproduces the observations, favours a narrow jet ‘core’ with half-opening angle θc = 2.1−1.4+2.4 deg (uncertainties hereon refer to 90% credible intervals from our fiducial ‘full sample’ analysis) whose peak luminosity, as seen on-axis, is distributed as a power law, p(Lc) ∝ Lc−A with A = 3.2−0.4+0.7, above a minimum isotropic-equivalent luminosity, Lc⋆ = 5−2+11 × 1051 erg s−1. For viewing angles larger than θc, the luminosity profile scales as a single power law, l ∝ θv−αL with αL = 4.7−1.4+1.2, with no evidence of a break, despite the model allowing for it. While the model implies an intrinsic ‘Yonetoku’ correlation between L and the peak photon energy, Ep, of the spectral energy distribution, its slope is somewhat shallower, Ep ∝ L0.4 ± 0.2, than the apparent one, and the normalisation is offset towards larger Ep due to selection effects. The implied local rate density of SGRBs (regardless of the viewing angle) is between about one hundred up to several thousand events per cubic gigaparsec per year, in line with the BNS merger rate density inferred from GW observations. Based on the model, we predict 0.2 to 1.3 joint GW+SGRB detections per year by the advanced GW detector network and Fermi/GBM during the O4 observing run.

Ab-initio General-relativistic Neutrino-radiation Hydrodynamics Simulations of Long-lived Neutron Star Merger Remnants to Neutrino Cooling Timescales

We perform the first 3D ab-initio general-relativistic neutrino-radiation hydrodynamics of a long-lived neutron star merger remnant spanning a fraction of its cooling timescale. We find that neutrino cooling becomes the dominant energy loss mechanism after the gravitational-wave dominated phase (∼20 ms postmerger). Electron flavor antineutrino luminosity dominates over electron flavor neutrino luminosity at early times, resulting in a secular increase of the electron fraction in the outer layers of the remnant. However, the two luminosities become comparable ∼20–40 ms postmerger. A dense gas of electron antineutrinos is formed in the outer core of the remnant at densities ∼1014.5 g cm−3, corresponding to temperature hot spots. The neutrinos account for ∼10% of the lepton number in this region. Despite the negative radial temperature gradient, the radial entropy gradient remains positive, and the remnant is stably stratified according to the Ledoux criterion for convection. A massive accretion disk is formed from the material squeezed out of the collisional interface between the stars. The disk carries a large fraction of the angular momentum of the system, allowing the remnant massive neutron star to settle to a quasi-steady equilibrium within the region of possible, stable, rigidly rotating configurations. The remnant is differentially rotating, but it is stable against the magnetorotational instability. Other MHD mechanisms operating on longer timescales are likely responsible for the removal of the differential rotation. Our results indicate the remnant massive neutron star is thus qualitatively different from a protoneutron stars formed in core-collapse supernovae.

Stellar Properties of Observed Stars Stripped in Binaries in the Magellanic Clouds

Massive stars (∼8–25 M ⊙) stripped of their hydrogen-rich envelopes via binary interaction are thought to be the main progenitors for merging neutron stars and stripped-envelope supernovae. We recently presented the discovery of the first set of such stripped stars in a companion paper. Here, we fit the spectra of 10 stars with new atmosphere models in order to constrain their stellar properties precisely. We find that the stellar properties align well with the theoretical expectations from binary evolution models for helium-core burning envelope-stripped stars. The fits confirm that the stars have high effective temperatures (T eff ∼ 50–100 kK), high surface gravities ( 5), and hydrogen-poor/helium-rich surfaces (X H,surf ∼ 0–0.4) while showing for the first time a range of bolometric luminosities (103–105 L ⊙), small radii (∼0.5–1 R ⊙), and low Eddington factors (Γ e ∼ 0.006–0.4). Using these properties, we derive intermediate current masses (∼1–8 M ⊙), which suggest that their progenitors were massive stars (∼5–25 M ⊙) and that a subset will reach core-collapse, leaving behind neutron stars or black holes. Using the model fits, we also estimate the emission rates of ionizing photons for these stars, which agree well with previous model expectations. Further, by computing models for a range of mass-loss rates, we find that the stellar winds are weaker than predicted by any existing scheme ( M ⊙ yr−1). The properties of this first sample of intermediate-mass helium stars suggest they both contain progenitors of type Ib and IIb supernovae, and provide important benchmarks for binary evolution and population synthesis models.

The Jet Opening Angle and Event Rate Distributions of Short Gamma-Ray Bursts from Late-time X-Ray Afterglows

We present a comprehensive study of 29 short gamma-ray bursts (SGRBs) observed ≈0.8−60 days postburst using Chandra and XMM-Newton. We provide the inferred distributions of the SGRB jet opening angles and true event rates to compare against neutron star merger rates. We perform a uniform analysis and modeling of their afterglows, obtaining 10 opening angle measurements and 19 lower limits. We report on two new opening angle measurements (SGRBs 050724A and 200411A) and eight updated values, obtaining a median value of 〈θ j〉 ≈ 6.°1 [−3.°2, +9.°3] (68% confidence on the full distribution) from jet measurements alone. For the remaining events, we infer θ j ≳ 0.°5–26°. We uncover a population of SGRBs with wider jets of θ j ≳ 10° (including two measurements of θ j ≳ 15°), representing ∼28% of our sample. Coupled with multiwavelength afterglow information, we derive a total true energy of 〈E true,tot〉 ≈ 1049–1050 erg, which is consistent with magnetohydrodynamic jet launching mechanisms. Furthermore, we determine a range for the beaming-corrected event rate of Gpc−3 yr−1, set by the inclusion of a population of wide jets on the low end, and the jet measurements alone on the high end. From a comparison with the latest merger rates, our results are consistent with the majority of SGRBs originating from binary neutron star mergers. However, our inferred rates are well above the latest neutron star–black hole merger rates, consistent with at most a small fraction of SGRBs originating from such mergers.

Modelling the spectra of the kilonova AT2017gfo – II. Beyond the photospheric epochs

ABSTRACT Binary neutron star mergers are the first confirmed site of element nucleosynthesis by the rapid neutron-capture process (r-process). The kilonova AT2017gfo is the only electromagnetic counterpart of a neutron star merger spectroscopically observed. We analyse the entire spectral sequence of AT2017gfo (from merger to +10.4 d) and identify seven emission-like features. We confirm that the prominent 1.08 $\mu{\text{m}}$ feature can be explained by the Sr ii near-infrared triplet evolving from a P-Cygni profile through to pure emission. We calculate the expected strength of the [Sr ii] doublet and show that its absence requires highly clumped ejecta. Near-infrared features at 1.58 and 2.07 $\mu {\text{m}}$ emerge after three days and become more prominent as the spectra evolve. We model these as optically thick P-Cygni profiles and alternatively as pure emission features (with FWHM ≃ 35 600 ± 6600 km s−1) and favour the latter interpretation. The profile of the strong 2.07 $\mu {\text{m}}$ emission feature is best reproduced with two lines, centred at 2.059 and 2.135 $\mu {\text{m}}$. We search for candidate ions for all prominent features in the spectra. Strong, permitted transitions of La iii, Ce iii, Gd iii, Ra ii, and Ac i are plausible candidates for the emission features. If any of these features are produced by intrinsically weak, forbidden transitions, we highlight candidate ions spanning the three r-process peaks. The second r-process peak elements Te and I have plausible matches to multiple features. We highlight the need for more detailed and quantitative atomic line transition data.

Effects of nuclear matter properties in neutron star mergers

ABSTRACT The dynamics in mergers of binary neutron star (BNS) systems depend sensitively on the equation of state (EOS) of dense matter. This has profound implications on the emission of gravitational waves (GWs) and the ejection of matter in the merger and post-merger phases and is thus of high interest for multimessenger astronomy. Today, a variety of nuclear EOSs are available with various underlying microphysical models. This calls for a study to focus on EOS effects from different physical nuclear matter properties and their influence on BNS mergers. We perform simulations of equal-mass BNS mergers with a set of nine different EOSs based on Skyrme density functionals. In the models, we systematically vary the effective nucleon mass, incompressibility, and symmetry energy at saturation density. This allows us to investigate the influence of specific nuclear matter properties on the dynamics of BNS mergers. We analyse the impact of these properties on the merger dynamics, the fate of the remnant, disc formation, ejection of matter, and GW emission. Our results indicate that some aspects of the merger, such as the frequencies of the post-merger GW spectrum and the shock-heated ejecta mass, are sensitive to the EOS around saturation density while others, such as the contraction of the remnant and the tidal ejecta mass, are sensitive to the behaviour towards higher densities, e.g. characterized by the slope of the pressure versus density. The detailed density dependence of the EOS thus needs to be taken into account to describe its influence on BNS mergers.

Impact of anisotropic ejecta on jet dynamics and afterglow emission in binary neutron-star mergers

ABSTRACT Binary neutron-stars mergers widely accepted as potential progenitors of short gamma-ray bursts. After the remnant of the merger has collapsed to a black hole, a jet is powered and may breakout from the the matter expelled during the collision and the subsequent wind emission. The interaction of the jet with the ejecta may affect its dynamics and the resulting electromagnetic counterparts. We here examine how an inhomogeneous and anisotropic distribution of ejecta affects such dynamics, dictating the properties of the jet-ejecta cocoon and of the afterglow radiated by the jet upon deceleration. More specifically, we carry out general-relativistic hydrodynamical simulations of relativistic jets launched within a variety of geometrically inhomogeneous and anisotropic distributions of ejected matter. We find that different anisotropies impact the variance of the afterglow light curves as a function of the jet luminosity and ejected mass. A considerable amount of the jet energy is deposited in the cocoon through the jet-ejecta interaction with a small but important dependence on the properties of the ejecta. Furthermore, all configurations show a two-component behaviour for the polar structure of the jet, with a narrow core at large energies and Lorentz factors and a shallow segment at high latitudes from the jet axis. Hence, afterglows measured on off-axis lines of sight could be used to deduce the properties of the ejected matter, but also that the latter need to be properly accounted for when modelling the afterglow signal and the jet-launching mechanisms.

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.

Bound Debris Expulsion from Neutron Star Merger Remnants

Many studies have found that neutron star mergers leave a fraction of the stars’ mass in bound orbits surrounding the resulting massive neutron star or black hole. This mass is a site of r-process nucleosynthesis and can generate a wind that contributes to a kilonova. However, comparatively little is known about the dynamics determining its mass or initial structure. Here we begin to investigate these questions, starting with the origin of the disk mass. Using tracer particle as well as discretized fluid data from numerical simulations, we identify where in the neutron stars the debris came from, the paths it takes in order to escape from the neutron stars’ interiors, and the times and locations at which its orbital properties diverge from those of neighboring fluid elements that end up remaining in the merged neutron star.

Updated Observing Scenarios and Multimessenger Implications for the International Gravitational-wave Networks O4 and O5

An advanced LIGO and Virgo’s third observing run brought another binary neutron star merger (BNS) and the first neutron-star black hole mergers. While no confirmed kilonovae were identified in conjunction with any of these events, continued improvements of analyses surrounding GW170817 allow us to project constraints on the Hubble Constant (H 0), the Galactic enrichment from r-process nucleosynthesis, and ultra-dense matter possible from forthcoming events. Here, we describe the expected constraints based on the latest expected event rates from the international gravitational-wave network and analyses of GW170817. We show the expected detection rate of gravitational waves and their counterparts, as well as how sensitive potential constraints are to the observed numbers of counterparts. We intend this analysis as support for the community when creating scientifically driven electromagnetic follow-up proposals. During the next observing run O4, we predict an annual detection rate of electromagnetic counterparts from BNS of () for the Zwicky Transient Facility (Rubin Observatory).