skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.
Attention:The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 7:00 AM ET to 7:30 AM ET on Friday, April 24 due to maintenance. We apologize for the inconvenience.


Title: Insights into binary neutron star merger simulations: a multi-code comparison
Gravitational Wave (GW) signals from binary neutron star (BNS) mergers provide critical insights into the properties of matter under extreme conditions. Due to the scarcity of observational data, numerical relativity (NR) simulations are indispensable for exploring these phenomena, without replacing the need for observational confirmation. However, simulating BNS mergers is a formidable challenge, and ensuring the consistency, reliability or convergence, especially in the post-merger, remains a work in progress. In this paper we assess the performance of current BNS merger simulations by analyzing open-source GW waveforms from five leading NR codes –SACRA,BAM,THC,Whisky, andSpEC. We focus on the accuracy of these simulations and on the effect of the equation of state on waveform predictions. We first check if different codes give similar results for similar initial data, then apply two methods to calculate convergence and quantify discretization errors. Lastly, we perform a thorough investigation into the effect of tidal interactions on key frequencies in the GW spectrum. We introduce a novel quasi-universal relation for the transient post-merger time, enhancing our understanding of remnant dynamics in this region. This detailed analysis clarifies agreements and discrepancies between these leading NR codes, and highlights necessary improvements for the advanced accuracy requirements of future GW detectors.  more » « less
Award ID(s):
2349289
PAR ID:
10667070
Author(s) / Creator(s):
;
Publisher / Repository:
IOP Publishing
Date Published:
Journal Name:
Classical and Quantum Gravity
Volume:
42
Issue:
9
ISSN:
0264-9381
Page Range / eLocation ID:
095012
Subject(s) / Keyword(s):
binary neutron star merger simulations numerical relativity gravitational waves multi-code comparison code convergence quasi-universal relations
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT The association of GRB170817A with a binary neutron star (BNS) merger has revealed that BNSs produce at least a fraction of short gamma-ray bursts (SGRBs). As gravitational wave (GW) detectors push their horizons, it is important to assess coupled electromagnetic (EM)/GW probabilities and maximize observational prospects. Here, we perform BNS population synthesis calculations with the code mobse, seeding the binaries in galaxies at three representative redshifts, $$z$$ = 0.01, 0.1, and 1 of the Illustris TNG50 simulation. The binaries are evolved and their locations numerically tracked in the host galactic potentials until merger. Adopting the microphysics parameters of GRB170817A, we numerically compute the broad-band light curves of jets from BNS mergers, with the afterglow brightness dependent on the local medium density at the merger site. We perform Monte Carlo simulations of the resulting EM population assuming either a random viewing angle with respect to the jet, or a jet aligned with the orbital angular momentum of the binary, which biases the viewing angle probability for GW-triggered events. We find a gamma-ray detection probability of $$\sim\!2{{\rm per\ cent}},10{{\rm per\ cent}},\mathrm{and}\ 40{{\rm per\ cent}}$$ for BNSs at $$z$$ = 1, 0.1, and 0.01, respectively, for the random case, rising to $$\sim\!75{{\rm per\ cent}}$$ for the $$z$$ = 0.01, GW-triggered aligned case. Afterglow detection probabilities of GW-triggered BNS mergers vary in the range of $$\sim \! 0.3 \!-\! 0.5{{\rm per\ cent}}$$, with higher values for aligned jets, and are comparable across the high- and low-energy bands, unlike gamma-ray-triggered events (cosmological SGRBs) which are significantly brighter at higher energies. We further quantify observational biases with respect to host galaxy masses. 
    more » « less
  2. Black holes (BHs) with masses between 3 5 M , produced by a binary neutron star (BNS) merger, can further pair up with a neutron star or BH and merge again within a Hubble time. However, the astrophysical environments in which this can happen and the rate of such mergers are open questions in astrophysics. Gravitational waves may play an important role in answering these questions. In this context, we discuss the possibility that the primary of the recent LIGO-Virgo-KAGRA binary GW230529_181500 (GW230529, in short) is the product of a previous BNS merger. Invoking numerical relativity (NR)-based fitting formulas that map the binary constituents’ masses and tidal deformabilities to the mass, spin, and kick velocity of the remnant BH, we investigate the potential parents of GW230529’s primary. Our calculations using NR fits based on BNS simulations reveal that the remnant of a high-mass BNS merger similar to GW190425 is consistent with the primary of GW230529. This argument is further strengthened by the gravitational wave-based merger rate estimation of GW190425-like and GW230529-like populations. We show that around 18% (median) of the GW190425-like remnants could become the primary component in GW230529-like mergers. The dimensionless tidal deformability parameter of the heavier neutron star in the parent binary is constrained to 67 61 + 163 at 90% credibility. Using estimates of the gravitational-wave kick imparted to the remnant, we also discuss the astrophysical environments in which these types of mergers can take place and the implications for their future observations. 
    more » « less
  3. Abstract For the first ∼3 yrs after the binary neutron star merger event GW 170817, the radio and X-ray radiation has been dominated by emission from a structured relativistic off-axis jet propagating into a low-density medium with n < 0.01 cm −3 . We report on observational evidence for an excess of X-ray emission at δt > 900 days after the merger. With L x ≈ 5 × 10 38 erg s −1 at 1234 days, the recently detected X-ray emission represents a ≥3.2 σ (Gaussian equivalent) deviation from the universal post-jet-break model that best fits the multiwavelength afterglow at earlier times. In the context of JetFit afterglow models, current data represent a departure with statistical significance ≥3.1 σ , depending on the fireball collimation, with the most realistic models showing excesses at the level of ≥3.7 σ . A lack of detectable 3 GHz radio emission suggests a harder broadband spectrum than the jet afterglow. These properties are consistent with the emergence of a new emission component such as synchrotron radiation from a mildly relativistic shock generated by the expanding merger ejecta, i.e., a kilonova afterglow. In this context, we present a set of ab initio numerical relativity binary neutron star (BNS) merger simulations that show that an X-ray excess supports the presence of a high-velocity tail in the merger ejecta, and argues against the prompt collapse of the merger remnant into a black hole. Radiation from accretion processes on the compact-object remnant represents a viable alternative. Neither a kilonova afterglow nor accretion-powered emission have been observed before, as detections of BNS mergers at this phase of evolution are unprecedented. 
    more » « less
  4. ABSTRACT We perform a comparative analysis of nucleosynthesis yields from binary neutron star (BNS) mergers, black hole-neutron star (BHNS) mergers, and core-collapse supernovae (CCSNe) with the goal of determining which are the most dominant sources of r-process enrichment observed in stars. We find that BNS and BHNS binaries may eject similar mass distributions of robust r-process nuclei post-merger (up to third peak and actinides, A ∼ 200−240), after accounting for the volumetric event rates. Magnetorotational (MR) CCSNe likely undergo a weak r-process (up to A ∼ 140) and contribute to the production of light element primary process (LEPP) nuclei, whereas typical thermal, neutrino-driven CCSNe only synthesize up to first r-process peak nuclei (A ∼ 80−90). We also find that the upper limit to the rate of MR CCSNe is $$\lesssim 1~{{\ \rm per\ cent}}$$ the rate of typical thermal CCSNe; if the rate was higher, then weak r-process nuclei would be overproduced. Although the largest uncertainty is from the volumetric event rate, the prospects are encouraging for confirming these rates in the next few years with upcoming surveys. Using a simple model to estimate the resulting kilonova light curve from mergers and our set of fiducial merger parameters, we predict that ∼7 BNS and ∼2 BHNS events will be detectable per year by the Vera C. Rubin Observatory (LSST), with prior gravitational wave (GW) triggers. 
    more » « less
  5. Abstract Gravitational-wave (GW) detectors are observing compact object mergers from increasingly far distances, revealing the redshift evolution of the binary black hole (BBH)—and soon the black hole–neutron star (BHNS) and binary neutron star (BNS)—merger rate. To help interpret these observations, we investigate the expected redshift evolution of the compact object merger rate from the isolated binary evolution channel. We present a publicly available catalog of compact object mergers and their accompanying cosmological merger rates from population synthesis simulations conducted with the COMPAS software. To explore the impact of uncertainties in stellar and binary evolution, our simulations use two-parameter grids of binary evolution models that vary the common-envelope efficiency with mass transfer accretion efficiency and supernova (SN) remnant mass prescription with SN natal kick velocity, respectively. We quantify the redshift evolution of our simulated merger rates using the local (z∼ 0) rate, the redshift at which the merger rate peaks, and the normalized differential rates (as a proxy for slope). We find that although the local rates span a range of ∼103across our model variations, their redshift evolutions are remarkably similar for BBHs, BHNSs, and BNSs, with differentials typically within a factor 3 and peaks ofz≈ 1.2–2.4 across models. Furthermore, several trends in our simulated rates are correlated with the model parameters we explore. We conclude that future observations of the redshift evolution of the compact object merger rate can help constrain binary models for stellar evolution and GW formation channels. 
    more » « less