Search NASA⌕ Search

SEARCH · Search NASA

Results for “binary neutron stars”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Quantifying modeling uncertainties when combining multiple gravitational-wave detections from binary neutron star sources

Here, with the increasing sensitivity of gravitational-wave detectors, we expect to observe multiple binary neutron-star systems through gravitational waves in the near future. The combined analysis of these gravitational-wave signals offers the possibility to constrain the neutron-star radius and the equation of state of dense nuclear matter with unprecedented accuracy. However, it is crucial to ensure that uncertainties inherent in the gravitational-wave models will not lead to systematic biases when information from multiple detections is combined. To quantify waveform systematics, we perform an extensive simulation campaign of binary neutron-star sources and analyze them with a set of four different waveform models. For our analysis with 38 simulations, we find that statistical uncertainties in the neutron-star radius decrease to ± 250 m (2% at 90% credible interval) but that systematic differences between currently employed waveform models can be twice as large. Hence, it will be essential to ensure that systematic biases will not become dominant in inferences of the neutron-star equation of state when capitalizing on future developments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Second release of the CoRe database of binary neutron star merger waveforms

Abstract We present the second data release of gravitational waveforms from binary neutron star (BNS) merger simulations performed by the Computational Relativity ( CoRe ) collaboration. The current database consists of 254 different BNS configurations and a total of 590 individual numerical-relativity simulations using various grid resolutions. The released waveform data contain the strain and the Weyl curvature multipoles up to ℓ = m = 4 . They span a significant portion of the mass, mass-ratio, spin and eccentricity parameter space and include targeted configurations to the events GW170817 and GW190425. CoRe simulations are performed with 18 different equations of state, seven of which are finite temperature models, and three of which account for non-hadronic degrees of freedom. About half of the released data are computed with high-order hydrodynamics schemes for tens of orbits to merger; the other half is computed with advanced microphysics. We showcase a standard waveform error analysis and discuss the accuracy of the database in terms of faithfulness. We present ready-to-use fitting formulas for equation of state-insensitive relations at merger (e.g. merger frequency), luminosity peak, and post-merger spectrum.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Signatures of deconfined quark phases in binary neutron star mergers

We investigate the quark deconfinement phase transition in the context of binary neutron star (BNS) mergers. We employ a new finite-temperature composition-dependent equation of state (EOS) with a first order phase transition between hadrons and deconfined quarks to perform numerical relativity simulations of BNS mergers. Here, the softening of the EOS due to the phase transition causes the merger remnants to be more compact and to collapse to a black hole (BH) at earlier times. The phase transition is imprinted on the postmerger gravitational wave (GW) signal duration, amplitude, and peak frequency. However, this imprint is only detectable for binaries with sufficiently long-lived remnants. Moreover, the phase transition does not result in significant deviations from quasi-universal relations for the postmerger GW peak frequency. We also study the impact of the phase transition on dynamical ejecta, remnant accretion disk masses, r-process nucleosynthetic yields and associated electromagnetic (EM) counterparts. While there are differences in the EM counterparts and nucleosynthesis yields between the purely hadronic models and the models with phase transitions, these can be primarily ascribed to the difference in remnant collapse time between the two. An exception is the non-thermal afterglow caused by the interaction of the fastest component of the dynamical ejecta and the interstellar medium, which is systematically boosted in the binaries with phase transition as a consequence of the more violent merger they experience.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Thermal Effects in Binary Neutron Star Mergers

Abstract We study the impact of finite-temperature effects in numerical-relativity simulations of binary neutron star mergers with microphysical equations of state and neutrino transport in which we vary the effective nucleon masses in a controlled way. We find that, as the specific heat is increased, the merger remnants become colder and more compact due to the reduced thermal pressure support. Using a full Bayesian analysis, we demonstrate that this effect will be measurable in the postmerger gravitational wave signal with next-generation observatories at signal-to-noise ratios of 15.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Handing off the outcome of binary neutron star mergers for accurate and long-term postmerger simulations

In this work, we perform binary neutron star (BNS) merger simulations in full dynamical general relativity with IllinoisGRMHD, on a Cartesian grid with adaptive-mesh refinement. After the remnant black hole has become nearly stationary, the evolution of the surrounding accretion disk on Cartesian grids over long timescales (~1s) is suboptimal, as Cartesian coordinates over-resolve the angular coordinates at large distances, and the accreting plasma flows obliquely across coordinate lines dissipating angular momentum artificially from the disk. To address this, we present the HandOff, a set of computational tools that enables the transfer of general relativistic magnetohydrodynamic (GRMHD) and spacetime data from IllinoisGRMHD to harm3d, a GRMHD code that specializes in modeling black hole accretion disks in static spacetimes over long timescales, making use of general coordinate systems with spherical topology. We demonstrate that the handoff allows for a smooth and reliable transition of GRMHD fields and spacetime data, enabling us to efficiently and reliably evolve BNS dynamics well beyond merger. We also discuss future plans, which involve incorporating advanced equations of state and neutrino physics into BNS simulations using the HandOff approach.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

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

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

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Testing claims of the GW170817 binary neutron star inspiral affecting β-decay rates

On August 17, 2017, the first gravitational wave signal from a binary neutron star inspiral (GW170817) was detected by Advanced LIGO and Advanced VIRGO. Here we present radioactive β-decay rates of three independent sources 44 Ti, 60 Co and 137 Cs, monitored during the same period by a precision experiment designed to investigate the decay of long-lived radioactive sources. We do not find any significant correlations between decay rates in a 5 h time interval following the GW170817 observation. This contradicts a previous claim published in this journal of an observed 2.5σ Pearson Correlation between fluctuations in the number of observed decays from two β-decaying isotopes ( 32 Si and 36 Cl) in the same time interval. By correcting for the choice of an arbitrary time interval, we find no evidence of a correlation above 1.5σ confidence. In addition, we argue that such analyses on correlations in arbitrary time intervals should always correct for the so-called Look-Elsewhere effect by quoting the global significance.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Magnetosphere Evolution and Precursor-driven Electromagnetic Signals in Merging Binary Neutron Stars

We detail new force-free simulations to investigate magnetosphere evolution and precursor electromagnetic (EM) signals from binary neutron stars. Our simulations fully follow a representative inspiral motion, capturing the intricate magnetospheric dynamics and their impact on EM outflows. We explore a range of stellar magnetic moment orientations and relative strengths, finding that the magnetospheres and Poynting flux evolution are strongly configuration dependent. The Poynting flux exhibits pulsations at twice the orbital frequency, 2Ω, and is highly anisotropic, following a power-law dependence on orbital frequency. The index ranges from 1 to 6, shaped by the intricate dynamics of the magnetospheres. Furthermore, we present the first computation of (1) the EM forces acting on the star surfaces, revealing the presence of torques that, for highly magnetized stars, could influence the orbital dynamics or break the crust; (2) the high-energy emission signals from these systems by adopting the established isolated pulsar theory. Assuming curvature radiation in the radiation-reaction limit, we find that photons could reach TeV–PeV energies in the last ∼ms for magnetic field strengths 10 10 –10 15 G. However, our analysis of single photon magnetic pair production suggests that these photons are unlikely to escape, with the MeV band emerging as a promising observational window for precursor high-energy emission. In this framework, we construct high-energy emission skymaps and light curves, exploring observational implications. Finally, we propose potential precursor radio emission and delayed afterglow echoes from magnetized outflows, which may contribute to late-time rebrightening in short gamma-ray bursts or to orphan afterglows.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Dynamical ejecta from binary neutron star mergers: Impact of a small residual eccentricity and of the equation of state implementation

Predicting the properties of the matter ejected during and after a neutron star merger is crucial to our ability to use electromagnetic observations of these mergers to constrain the masses of the neutron stars, the equation of state of dense matter, and the role of neutron star mergers in the enrichment of the Universe in heavy elements. Furthermore, our ability to reliably provide such predictions is however limited by a broad range of factors, including the finite resolution of numerical simulations, their treatment of magnetic fields, neutrinos, and neutrino-matter interactions, and the approximate modeling of the equation of state of dense matter. In this manuscript, we study specifically the role that a small residual eccentricity and different implementations of the same equation of state have on the matter ejected during the merger of a 1.3M ⊙ –1.4M ⊙ binary neutron star system. We find that a residual eccentricity e ~ 0.01, as measured ~ 4–6 orbits before merger, causes O(25%–30%) changes in the amount of ejected mass, mainly due to changes in the amount of matter ejected as a result of core bounces during merger. We note that O(1%) residual eccentricities have regularly been used in binary neutron star merger simulations as proxy for circular binaries, potentially creating an additional source of error in predictions for the mass of the dynamical ejecta.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Chasing Gamma-Ray Signals from Binary Neutron Star Coalescences with the Cherenkov Telescope Array: Prospects and Observing Strategies

The detection of gravitational waves (GWs) from a binary neutron star (BNS) merger by Advanced LIGO and Advanced Virgo (GW170817), together with its electromagnetic counterpart, the short gamma-ray burst GRB 170817A, heralded the birth of multimessenger astronomy. The detection of TeV emission from GRBs motivates follow-up observations with the Cherenkov Telescope Array Observatory (CTAO), which is ideal for detecting such signals due to its unprecedented sensitivity, rapid response, and wide-field survey capabilities. The aim of this work is to evaluate GeV–TeV GW follow-up strategies for CTAO using a multistep simulation pipeline and to estimate the expected rate of joint GW–GRB detections during observing run O5. Using a simulated sample of BNS systems with corresponding GW detections, gamma-ray emission is simulated through phenomenological prescriptions based on the observed population of short GRBs, including off-axis jet scenarios. CTAO observations are simulated to account for instrument response, sky tiling strategies, integration times, and varying observing conditions. Strategies with variable and constant integration times are investigated. We find that, via an optimized follow-up strategy, about 5% of simulated GW-associated short GRBs produce GeV–TeV radiation detectable by CTAO. Detectability is strongly influenced by the jet opening angle and viewing angle, suggesting that even rough estimates of the viewing angle in GW alerts could enhance targeting. This framework motivates future follow-ups of GW-detectable events, including neutron star–black hole mergers, and further supports the development of advanced strategies incorporating galaxy distributions and synergies with future detectors such as the Einstein Telescope.

Abe, S. [University of Tokyo] (ORCID:0000000172503↗
DOE OSTI2025

The Gravity Collective: A Comprehensive Analysis of the Electromagnetic Search for the Binary Neutron Star Merger GW190425

We present an ultraviolet to infrared search for the electromagnetic (EM) counterpart to GW190425, the second ever binary neutron star merger discovered by the LIGO-Virgo-KAGRA Collaboration. GW190425 was more distant and had a larger localization area than GW170817, so we use a new tool, Teglon, to redistribute the GW190425 localization probability in the context of galaxy catalogs within the final localization volume. We derive a 90th percentile area of 6688 deg 2 , a ∼1.5× improvement relative to the LIGO/Virgo map, and show how Teglon provides an order-of-magnitude boost to the search efficiency of small (≤1 deg 2 ) field-of-view instruments. We combine our data with a large, publicly reported imaging data set, covering 9078.59 deg 2 of unique area and 48.13% of the LIGO/Virgo-assigned localization probability, to calculate the most comprehensive kilonova (KN), short gamma-ray burst (sGRB) afterglow, and model-independent constraints on the EM emission from a hypothetical counterpart to GW190425 to date under the assumption that no counterpart was found in these data. If the counterpart were similar to AT 2017gfo, there would be a 28.4% chance of it being detected in the combined data set. We are relatively insensitive to an on-axis sGRB, and rule out a generic transient with a similar peak luminosity and decline rate as AT 2017gfo to 30% confidence. Finally, across our new imaging and publicly reported data, we find 28 candidate optical counterparts that we cannot rule out as being associated with GW190425, finding that four such counterparts discovered within the localization volume and within 5 days of merger exhibit luminosities consistent with a KN.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Impact of moment-based, energy integrated neutrino transport on microphysics and ejecta in binary neutron star mergers

We present an extensive study of the effects of neutrino transport in three-dimensional general relativistic radiation hydrodynamics (GRHD) simulations of binary neutron star (BNS) mergers using our moment-based, energy-integrated neutrino radiation transport (M1) scheme. Here, we consider a total of eight BNS configurations, while varying equation of state models, mass ratios, and grid resolutions, for a total of 16 simulations. We find that M1 neutrino transport is crucial in modeling the local absorption of neutrinos and the deposition of lepton number throughout the medium. We provide an in-depth look at the effects of neutrinos on the fluid dynamics and luminosity during the late inspiral and postmerger phases, the properties of ejecta and outflow, and the postmerger nucleosynthesis. The simulations presented in this work comprise an extensive study of the combined effect of the equation of state and M1 neutrino transport in GRHD simulations of BNS mergers, and establish that the solution provided by our M1 scheme is robust across system properties.

150 ≤ A ≤ 18959 ≤ A ≤ 8990 ≤ A ≤ 149↗
DOE OSTI2022

Binary neutron star mergers as a probe of quark-hadron crossover equations of state

It is anticipated that the gravitational radiation detected in future gravitational wave (GW) detectors from binary neutron star (NS) mergers can probe the high-density equation of state (EOS). In this work, we perform the first simulations of binary NS mergers which adopt various parametrizations of the quark-hadron crossover (QHC) EOS. These are constructed from combinations of a hadronic EOS ($n_b$<2 $n_0$) and a quark-matter EOS ($n_b$>5 $n_0$), where $n_b$ and $n_0$ are the baryon number density and the nuclear saturation density, respectively. At the crossover densities (2 $n_0$ < $n_b$ < 5 $n_0$) the QHC EOSs continuously soften, while remaining stiffer than hadronic and first-order phase transition EOSs, achieving the stiffness of strongly correlated quark matter. This enhanced stiffness leads to significantly longer lifetimes of the postmerger NS than that for a pure hadronic EOS. We find a dual nature of these EOSs such that their maximum chirp GW frequencies $f_{\text{max}}$ fall into the category of a soft EOS while the dominant peak frequencies ($f_{\text{peak}}$) of the postmerger stage fall in between that of a soft and stiff hadronic EOS. An observation of this kind of dual nature in the characteristic GW frequencies will provide crucial evidence for the existence of strongly interacting quark matter at the crossover densities for QCD.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

A DESGW Search for the Electromagnetic Counterpart to the LIGO/Virgo Gravitational-wave Binary Neutron Star Merger Candidate S190510g

In this work, we present the results from a search for the electromagnetic counterpart of the LIGO/Virgo event S190510g using the Dark Energy Camera (DECam). S190510g is a binary neutron star (BNS) merger candidate of moderate significance detected at a distance of 227 ± 92 Mpc and localized within an area of 31 (1166) square degrees at 50% (90%) confidence. While this event was later classified as likely nonastrophysical in nature within 30 hours of the event, our short latency search and discovery pipeline identified 11 counterpart candidates, all of which appear consistent with supernovae following offline analysis and spectroscopy by other instruments. Later reprocessing of the images enabled the recovery of six more candidates. Additionally, we implement our candidate selection procedure on simulated kilonovae and supernovae under DECam observing conditions (e.g., seeing and exposure time) with the intent of quantifying our search efficiency and making informed decisions on observing strategy for future similar events. This is the first BNS counterpart search to employ a comprehensive simulation-based efficiency study. We find that using the current follow-up strategy, there would need to be 19 events similar to S190510g for us to have a 99% chance of detecting an optical counterpart, assuming a GW170817-like kilonova. We further conclude that optimization of observing plans, which should include preference for deeper images over multiple color information, could result in up to a factor of 1.5 reduction in the total number of follow-ups needed for discovery.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Performance-portable Binary Neutron Star Mergers with AthenaK

We introduce an extension to the AthenaK code for general-relativistic magnetohydrodynamics (GRMHD) in dynamical spacetimes using a 3+1 conservative Eulerian formulation. Like the fixed-spacetime GRMHD solver, we use standard finite-volume methods to evolve the fluid and a constrained-transport scheme to preserve the divergence-free constraint for the magnetic field. We also utilize a first-order flux correction (FOFC) scheme to reduce the need for an artificial atmosphere and optionally enforce a maximum principle to improve robustness. We demonstrate the accuracy of AthenaK using a set of standard tests in flat and curved spacetimes. Using a SANE accretion disk around a Kerr black hole, we compare the new solver to the existing solver for stationary spacetimes using the so-called "HARM-like" formulation. We find that both formulations converge to similar results. We also include the first published binary neutron star (BNS) mergers performed on graphical processing units (GPUs). Thanks to the FOFC scheme, our BNS mergers maintain a relative error of $\mathcal{O}$(10 –11 ) or better in baryon mass conservation up to collapse. Finally, we perform scaling tests of AthenaK on OLCF Frontier, where we show excellent weak scaling of ≥80% efficiency up to 32,768 GPUs and 74% up to 65,536 GPUs for a GRMHD problem in dynamical spacetimes with six levels of mesh refinement. AthenaK achieves an order-of-magnitude speedup using GPUs compared to CPUs, demonstrating that it is suitable for performing numerical relativity problems on modern exascale resources.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Binary neutron star mergers as potential sources for ultrahigh-energy cosmic rays and high-energy neutrinos

Recent studies suggest that the most energetic cosmic rays, exceeding 100 EeV, may primarily consist of 𝑟-process nuclei. This highlights binary neutron star mergers and collapsars as promising sources of ultrahigh-energy cosmic rays (UHECRs). Building on these insights, we examine the conditions that facilitate the efficient production of UHE 𝑟-process nuclei during the prompt radiation (PR), extended emission (EE), and plateau emission phases of short gamma-ray bursts (sGRBs) following neutron star mergers. Our study reveals that jets associated with the PR phase, characterized by typical bulk Lorentz factors (≳ 400–500), dissipation radii, and magnetic field strengths, can accelerate 𝑟-process nuclei to energies ≳ 100 EeV, while preserving them during propagation within the source. Additionally, we investigate the production of HE neutrinos from photomeson and hadronic interactions, as well as from the 𝛽 decay of accelerated 𝑟-process nuclei. We find that the HE neutrino fluxes from sGRBs, mainly produced via photomeson interactions, are significantly limited to preserve the accelerated heavy nuclei, leading to lower fluxes than the predictions without allowing for contributions to UHECRs. Our results suggest that sGRBs may potentially contribute to UHECRs during the PR phase and to HE neutrinos during the EE phase—a scenario that can be tested by future neutrino observatories.

Cosmic ray sources↗
DOE OSTI2024

Toward Accelerated Nuclear-physics Parameter Estimation from Binary Neutron Star Mergers: Emulators for the Tolman–Oppenheimer–Volkoff Equations

Abstract Gravitational-wave observations of binary neutron-star (BNS) mergers have the potential to revolutionize our understanding of the nuclear equation of state (EOS) and the fundamental interactions that determine its properties. However, Bayesian parameter estimation frameworks do not typically sample over microscopic nuclear-physics parameters that determine the EOS. One of the major hurdles in doing so is the computational cost involved in solving the neutron-star structure equations, known as the Tolman–Oppenheimer–Volkoff (TOV) equations. In this paper, we explore approaches to emulating solutions for the TOV equations: multilayer perceptrons (MLPs), Gaussian processes, and a data-driven variant of the reduced basis method (RBM). We implement these emulators for three different parameterizations of the nuclear EOS, each with a different degree of complexity represented by the number of model parameters. We find that our MLP-based emulators are generally more accurate than the other two algorithms, whereas the RBM results in the largest speedup with respect to the full high-fidelity TOV solver. We employ these emulators for a simple parameter inference using a potentially loud BNS observation and show that the posteriors predicted by our emulators are in excellent agreement with those obtained from the full TOV solver.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Study of the agreement between binary neutron star ejecta models derived from numerical relativity simulations

Neutron star mergers have recently become a tool to study extreme gravity, nucleosynthesis, and the chemical composition of the Universe. To date, there has been one joint gravitational and electromagnetic observation of a binary neutron star merger, GW170817, as well as a solely gravitational observation, GW190425. In order to accurately identify and interpret electromagnetic signals of neutron star mergers, better models of the matter outflows generated by these mergers are required. We compare a series of ejecta models to see where they provide strong constraints on the amount of ejected mass expected from a system, and where systematic uncertainties in current models prevent us from reliably extracting information from observed events. We also examine 2396 neutron star equations of state compatible with GW170817 to see whether a given ejecta mass could be reasonably produced with a neutron star of said equation of state, and whether different ejecta models provide consistent predictions. We find that the difference between models is often comparable to or larger than the error generally assumed for these models, implying better constraints on the models are needed. In conclusion, we also note that the extrapolation of outflow models outside of their calibration window, while commonly needed to analyze gravitational wave events, is extremely unreliable and occasionally leads to completely unphysical results.

79 ASTRONOMY AND ASTROPHYSICS↗