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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Neutron Stars and the Nuclear Matter Equation of State

Neutron stars provide a window into the properties of dense nuclear matter. Several recent observational and theoretical developments provide powerful constraints on their structure and internal composition. Among these are the first observed binary neutron star merger, GW170817, whose gravitational radiation was accompanied by electromagnetic radiation from a short γ-ray burst and an optical afterglow believed to be due to the radioactive decay of newly minted heavy r-process nuclei. These observations give important constraints on the radii of typical neutron stars and on the upper limit to the neutron star maximum mass and complement recent pulsar observations that established a lower limit. Pulse-profile observations by the Neutron Star Interior Composition Explorer (NICER) X-ray telescope provide an independent, consistent measure of the neutron star radius. Theoretical many-body studies of neutron matter reinforce these estimates of neutron star radii. Studies using parameterized dense matter equations of state (EOSs) reveal several EOS-independent relations connecting global neutron star properties.

Physics↗
DOE OSTI2023

A Unified Picture of Short and Long Gamma-Ray Bursts from Compact Binary Mergers

Abstract The recent detections of the ∼10 s longγ-ray bursts (GRBs) 211211A and 230307A followed by softer temporally extended emission (EE) and kilonovae point to a new GRB class. Using state-of-the-art first-principles simulations, we introduce a unifying theoretical framework that connects binary neutron star (BNS) and black hole–NS (BH–NS) merger populations with the fundamental physics governing compact binary GRBs (cbGRBs). For binaries with large total masses,M tot ≳ 2.8M ⊙ , the compact remnant created by the merger promptly collapses into a BH surrounded by an accretion disk. The duration of the pre-magnetically arrested disk (MAD) phase sets the duration of the roughly constant power cbGRB and could be influenced by the disk mass,M d . We show that massive disks (M d ≳ 0.1M ⊙ ), which form for large binary mass ratiosq≳ 1.2 in BNS orq≲ 3 in BH–NS mergers, inevitably produce 211211A-like long cbGRBs. Once the disk becomes MAD, the jet power drops with the mass accretion rate as M ̇ ∼ t − 2 , establishing the EE decay. Two scenarios are plausible for short cbGRBs. They can be powered by BHs with less massive disks, which form for otherqvalues. Alternatively, for binaries withM tot ≲ 2.8M ⊙ , mergers should go through a hypermassive NS (HMNS) phase, as inferred for GW170817. Magnetized outflows from such HMNSs, which typically live for ≲1 s, offer an alternative progenitor for short cbGRBs. The first scenario is challenged by the bimodal GRB duration distribution and the fact that the Galactic BNS population peaks at sufficiently low masses that most mergers should go through an HMNS phase.

Astronomy & Astrophysics↗
DOE OSTI2024

Impact of dark matter on tidal signatures in neutron star mergers with the Einstein Telescope

If dark matter (DM) accumulates inside neutron stars (NS), it changes their internal structure and causes a shift of the tidal deformability from the value predicted by the dense-matter equation of state (EOS). In principle, this shift could be observable in the gravitational-wave (GW) signal of binary neutron star (BNS) mergers. We investigate the effect of fermionic, noninteracting DM when observing a large number of GW events from DM-admixed BNSs with the precision of the proposed Einstein telescope (ET). Specifically, we study the impact on the recovery of the baryonic EOS and whether DM properties can be constrained. For this purpose, we create event catalogs of BNS mock events with DM fraction up to 1%, from which we reconstruct the posterior uncertainties with the Fisher matrix approach. Using this data, we perform joint Bayesian inference on the baryonic EOS, DM particle mass, and DM particle fraction in each event. Here, our results reveal that when falsely ignoring DM effects, the EOS posterior is biased toward softer EOSs, though the offset is rather small. Further, we find that within our assumptions of our DM model and population, ET will likely not be able to test the presence of DM in BNSs, even when combining many events and adding Cosmic Explorer (CE) to the next-generation detector network. Likewise, the potential constraints on the DM particle mass will remain weak because of degeneracies with the fraction and EOS.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

What Can We Learn about the Unstable Equation-of-state Branch from Neutron Star Mergers?

Abstract The equation of state (EOS) of dense strongly interacting matter can be probed by astrophysical observations of neutron stars (NS), such as X-ray detections of pulsars or the measurement of the tidal deformability of NSs during the inspiral stage of NS mergers. These observations constrain the EOS at most up to the density of the maximum-mass configuration, n TOV , which is the highest density that can be explored by stable NSs for a given EOS. However, under the right circumstances, binary neutron star (BNS) mergers can create a postmerger remnant that explores densities above n TOV . In this work, we explore whether the EOS above n TOV can be measured from gravitational-wave or electromagnetic observations of the postmerger remnant. We perform a total of 25 numerical-relativity simulations of BNS mergers for a range of EOSs and find no case in which different descriptions of the matter above n TOV have a detectable impact on postmerger observables. Hence, we conclude that the EOS above n TOV can likely not be probed through BNS merger observations for the current and next generation of detectors.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Nucleosynthesis in outflows of compact objects and detection prospects of associated kilonovae

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.

Astronomy & Astrophysics↗
DOE OSTI2025

Effects of Landau quantization on neutrino emission and absorption

Some neutron stars known as magnetars possess very strong magnetic fields, with surface fields as large as 10 15 G and internal fields that are possibly stronger. Recent observations of the radio pulsar GLEAM-X J1627 suggest it may have a surface field as strong as 10 16 G. In the presence of a strong magnetic field, Leinson and Perez showed and Baiko and Yakovlev confirmed that the energy levels of electrons and protons are quantized and the direct Urca process allows neutron stars to cool rapidly, even at low density. For the case of magnetic fields 𝐵 ≳ 10 16 G, we find features in the emissivity due to energy quantization that are not captured by the frequently employed quasiclassical approximation where energy levels are treated as nearly continuous. Resonances can result in amplification of the neutrino emissivity at specific densities compared to a calculation that neglects quantization, particularly at low temperature. These effects are not important for the thermal evolution of an entire neutron star, but may be relevant for phenomena that depend on behavior at specific densities. Here, we present a fully relativistic calculation of the direct Urca rate in a strong magnetic field using the standard V-A weak Lagrangian incorporating mean field nuclear effects and discuss approaches to the numerical challenge the modified wave functions present and a new semianalytic approximation. These tools are also applicable to calculating neutrino opacities in strong magnetic fields in the ejecta of binary neutron star mergers. We calculate the opacities for neutrinos capturing on free nucleons at subsaturation densities and temperatures exceeding an MeV. We find an enhancement to capture processes of the lowest energy neutrinos by an order of magnitude or more due to suppression of electron Pauli blocking in the case of capture on neutrons, and from the effect of the nucleon magnetic moments in the case of capture on protons.

astrophysical electromagnetic fields↗
DOE OSTI2024

Do black holes remember what they are made of?

We study the ringdown signal of black holes formed in prompt-collapse binary neutron star mergers. We analyze data from 47 numerical relativity simulations. We show that the ($l$ = 2, $m$ = 2) and ($l$ = 2, $m$ = 1) multipoles of the gravitational wave signal are well fitted by decaying damped exponentials, as predicted by black-hole perturbation theory. We show that the ratio of the amplitude in the two modes depends on the progenitor binary mass ratio q and reduced tidal parameter $\overline{Λ}$. Unfortunately, the numerical uncertainty in our data is too large to fully quantify this dependency. If confirmed, these results will enable novel tests of general relativity in the presence of matter with next-generation gravitational-wave observatories.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Quantum maximum entropy closure for small flavor coherence

Quantum angular moment transport schemes are an important avenue toward describing neutrino flavor mixing phenomena in dense astrophysical environments such as supernovae and merging neutron stars. Successful implementation will require new closure relations that go beyond those used in classical transport. In this paper, we derive the first analytic expression for a quantum M1 closure, valid in the limit of small flavor coherence, based on the maximum entropy principle. We verify that the resulting closure relation has the appropriate limits and characteristic speeds in the diffusive and free-streaming regimes. We then use this new closure in a moment linear stability analysis to search for fast flavor instabilities in a binary neutron star merger simulation and find better results as compared with previously designed, ad hoc , semiclassical closures.

astrophysical & cosmological simulations↗
DOE OSTI2021

Projecting the likely importance of weak-interaction-driven bulk viscosity in neutron star mergers

ABSTRACT In this work, we estimate how much bulk viscosity driven by Urca processes is likely to affect the gravitational wave signal of a neutron star coalescence. In the late inspiral, we show that bulk viscosity affects the binding energy at fourth post-Newtonian order. Even though this effect is enhanced by the square of the gravitational compactness, the coefficient of bulk viscosity is likely too small to lead to observable effects in the waveform during the late inspiral, when only considering the orbital motion itself. In the post-merger, however, the characteristic time-scales and spatial scales are different, potentially leading to the opposite conclusion. We post-process data from a state-of-the-art equal-mass binary neutron star merger simulation to estimate the effects of bulk viscosity (which was not included in the simulation itself). In that scenario, we find that bulk viscosity can reach high values in regions of the merger. We compute several estimates of how much it might directly affect the global dynamics of the considered merger scenario, and find that it could become significant. Even larger effects could arise in different merger scenarios or in simulations that include non-linear effects. This assessment is reinforced by a quantitative comparison with relativistic heavy-ion collisions where such effects have been explored extensively.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

On the Minimum Radius of Very Massive Neutron Stars

Prospects of establishing the radii of massive neutron stars in PSR J1614 – 2230 and PSR J0740 + 6620 from Neutron Star Interior Composition ExploreR and Chandra observatories hold the potential to constrain the equation of state (EoS) of matter to densities well beyond those encountered in canonical stars of mass $\sim 1.4\,{M}_{\odot }$. In this work, we investigate the relation between the radii of very massive neutron stars up to the maximum mass, ${M}_{\max }$, supported by dense matter EoSs. Results from models with hadronic matter are contrasted with those that include a first-order hadron-to-quark phase transition. We find that a lower bound on ${M}_{\max }$ with an upper bound on the radius of massive pulsars serves to rule out quark matter that is too soft, and an upper bound on ${M}_{\max }$ with a lower bound on the radius of massive pulsars strongly disfavors a transition into quark matter that is too stiff appearing at low densities. In conclusion, the complementary role played by radius inferences from future gravitational-wave events of inspiraling binary neutron stars is also briefly discussed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Exploring dark forces with multimessenger studies of extreme mass ratio inspirals

Abstract The exploration of dark sector interactions via gravitational waves (GWs) from binary inspirals has been a subject of recent interest. We study dark forces using extreme mass ratio inspirals (EMRIs), pointing out two issues of interest. Firstly, the innermost stable circular orbit (ISCO) of the EMRI, which sets the characteristic length scale of the system and hence the dark force range to which it exhibits enhanced sensitivity, probes force mediator masses that complement those studied with supermassive black hole (SMBH) or neutron star binaries. The LISA mission (the proposedμAres detector) will probe mediators with massesm V ∼ 10 -16 eV (m V ∼ 10 -18 eV), corresponding to ISCOs of 10 6 M ⊙ (10 8 M ⊙ ) central SMBHs. Secondly, while the sensitivity to dark couplings is typically limited by the uncertainty in the binary component masses, independent mass measurements of the central SMBH through reverberation mapping campaigns or the motion of dynamical tracers enable one to break this degeneracy. Our results therefore highlight the necessity for coordinated studies, loosely referred to as “multimessenger”, between futureμHz- mHz GW observatories and ongoing and forthcoming SMBH mass measurement campaigns, including OzDES-RM, SDSS-RM, and SDSS-V Black Hole Mapper.

Astronomy & Astrophysics↗
DOE OSTI2020

The Binary–Host Connection: Astrophysics of Gravitational-Wave Binaries from Host Galaxy Properties

Gravitational waves from the merger of binary neutron stars (BNSs) are accompanied by electromagnetic counterparts, making it possible to identify the associated host galaxy. In this work, we explore how properties of the hosts relate to the astrophysical processes leading to the mergers. It is thought that the BNS merger rate within a galaxy at a given epoch depends primarily on the galaxy's star formation history, as well as the underlying merger time-delay distribution of the binary systems. The stellar history of a galaxy, meanwhile, depends on the cosmological evolution of the galaxy through time, and is tied to the growth of structure in the universe. We study the hosts of BNS mergers in the context of structure formation by populating the UniverseMachine simulations with gravitational wave (GW) events, based on a simple time-delay model. We find that different time-delay distributions predict different properties of the associated host galaxies, including the distributions of stellar mass, star formation rate, halo mass, and local and large-scale clustering of hosts. Moreover, BNSs merging today with short delay times occur preferentially in hosts with high star formation rates, while those with long delay times live in dense regions within massive halos that have low star formation. Here, we show that with ${ \mathcal O }(10)$ events from current GW detector networks, it is possible to make preliminary distinctions between formation channels which trace stellar mass, halo mass, or star formation rate. We also find that strategies to follow-up GW events with electromagnetic telescopes can be significantly optimized using the clustering properties of their hosts.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Mass measurements of neutron-rich gallium isotopes refine production of nuclei of the first r -process abundance peak in neutron-star merger calculations

We report mass measurements of neutron-rich Ga isotopes 80–85 Ga with TRIUMF's Ion Trap for Atomic and Nuclear science. The measurements determine the masses of 80–83 Ga in good agreement with previous measurements. Here, the masses of 84 Ga and 85 Ga were measured for the first time. Uncertainties between 25 and 48 keV were reached. The new mass values reduce the nuclear uncertainties associated with the production of A ≈ 84 isotopes by the r-process for astrophysical conditions that might be consistent with a binary neutron star (BNS) merger producing a blue kilonova. Our nucleosynthesis simulations confirm that BNS merger may contribute to the first abundance peak under moderate neutron-rich conditions with electron fractions Y e = 0.35–0.38.

59 ≤ A ≤ 89↗
DOE OSTI2025

From Existing and New Nuclear and Astrophysical Constraints to Stringent Limits on the Equation of State of Neutron-Rich Dense Matter

Through continuous progress in nuclear theory and experiment and an increasing number of neutron-star (NS) observations, a multitude of information about the equation of state (EOS) for matter at extreme densities is available. To constrain the EOS across its entire density range, this information needs to be combined consistently. However, the impact and model dependency of individual observations vary. Given their growing number, assessing the various methods is crucial to compare the respective effects on the EOS and discover potential biases. For this purpose, we present a broad compendium of different constraints and apply them individually to a large set of EOS candidates within a Bayesian framework. Specifically, we explore different ways of how chiral effective field theory and perturbative quantum chromodynamics can be used to place a likelihood on EOS candidates. We also investigate the impact of nuclear experimental constraints, as well as different radio and x-ray observations of NS masses and radii. This is augmented by reanalyses of the existing data from binary neutron star coalescences, in particular of GW170817, with improved models for the tidal waveform and kilonova light curves, which we also utilize to construct a tight upper limit of 2.39 M ⊙ on the TOV mass based on GW170817’s remnant. Our diverse set of constraints is eventually combined to obtain stringent limits on NS properties. We organize the combination in a way to distinguish between constraints where the systematic uncertainties are deemed small and those that rely on less conservative assumptions. For the former, we find the radius of the canonical 1.4 M ⊙ neutron star to be R 1.4 = 12.2 6 − 0.91 + 0.80 km and the TOV mass at M TOV = 2.2 5 − 0.22 + 0.42 M ⊙ (95% credibility). Including all the presented constraints yields R 1.4 = 12.2 0 − 0.48 + 0.50 km and M TOV = 2.3 0 − 0.20 + 0.07 M ⊙ . When comparing these limits to individual data points, we find that the quoted radius of HESS J1731-347 displays noticeable tension with other constraints. Constraining microphysical properties of the EOS proves more challenging. For instance, the symmetry energy slope is restricted to L sym = 48 − 25 + 21 MeV , where this constraint is mainly dominated by our reanalysis of the PREX-II and CREX experiment. Published by the American Physical Society 2025

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Photon temporal-mode readout for inference of neutron star merger remnant gravitational waves

Gravitational waves emitted after neutron star binary coalescences and the information they carry about dense matter are a high-priority target for next-generation detectors. Even though such detectors are expected to observe millions of signals, detectable postmerger emission will remain rare. Here, in this work, we explore postmerger detectability and inference through an alternative detector readout scheme for data dominated by quantum-noise, which is the case above 1 kHz; photon-counting. In such a readout, signals and noise become quantized into discrete distributions corresponding to the detection of single photons measured in a chosen basis of modes. Through simulated data, we demonstrate that photon counting can be efficient even for weak signals. We find ∼1 in 100 signals with a postmerger signal-to-noise ratio of 0.2 can result in a single photon and thus be detected. Furthermore, after 2 ×10 4 signals—equivalent to 10 −2 to 1.5 years of observation—photon counting results in a twofold improvement in the measurement of the radius of a 1.6⁢𝑀 ⊙ neutron star. Constraints can be further tightened if the detector classical noise is reduced. Photon counting offers a promising alternative to traditional homodyne readout techniques for extracting information from low signal-to-noise ratio postmerger signals.

gravitational wave detection↗
DOE OSTI2023

Modelling populations of kilonovae

Abstract The 2017 detection of a kilonova coincident with gravitational-wave emission has identified neutron star mergers as the major source of the heaviest elements and dramatically constrained alternative theories of gravity. Observing a population of such sources has the potential to transform cosmology, nuclear physics, and astrophysics. However, with only one confident multi-messenger detection currently available, modelling the diversity of signals expected from such a population requires improved theoretical understanding. In particular, models that are quick to evaluate and are calibrated with more detailed multi-physics simulations are needed to design observational strategies for kilonovae detection and to obtain rapid-response interpretations of new observations. We use grey-opacity models to construct populations of kilonovae, spanning ejecta parameters predicted by numerical simulations. Our modelling focuses on wavelengths relevant for upcoming optical surveys, such as the Rubin Observatory Legacy Survey of Space and Time (LSST). In these simulations, we implement heating rates that are based on nuclear reaction network calculations. We create a Gaussian-process emulator for kilonova grey opacities, calibrated with detailed radiative transfer simulations. Using recent fits to numerical relativity simulations, we predict how the ejecta parameters from binary neutron star (BNS) mergers shape the population of kilonovae, accounting for the viewing-angle dependence. Our simulated population of BNS mergers produce peak i-band absolute magnitudes of −20 ≤ Mi ≤ −11. A comparison with detailed radiative transfer calculations indicates that further improvements are needed to accurately reproduce spectral shapes over the full light curve evolution.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

BSE versus StarTrack: Implementations of new wind, remnant-formation, and natal-kick schemes in NBODY7 and their astrophysical consequences

As a result of their formation via massive single and binary stellar evolution, the masses of stellar-remnant black holes (BH) are subjects of great interest in this era of gravitational-wave detection from binary black hole (BBH) and binary neutron star merger events. In this work, we present new developments in the stellar-remnant formation and related schemes of the current N-body evolution program NBODY7. We demonstrate that the newly implemented stellar-wind and remnant-formation schemes in the stellar-evolutionary sector or BSE of the NBODY7 code, such as the “rapid” and the “delayed” supernova (SN) schemes along with an implementation of pulsational-pair-instability and pair-instability supernova (PPSN/PSN), now produce neutron star (NS) and BH masses that agree nearly perfectly, over large ranges of zero-age-main-sequence (ZAMS) mass and metallicity, with those from the widely recognised StarTrack population-synthesis program. We also demonstrate the new, recipe-based implementations of various widely debated mechanisms of natal kicks on NSs and BHs, such as “convection-asymmetry-driven”, “collapse-asymmetry-driven”, and “neutrino-emission-driven” kicks, in addition to a fully consistent implementation of the standard, fallback-dependent, momentum-conserving natal kick. All the above newly implemented schemes are also shared with the standalone versions of SSE and BSE. Furthermore, all these demonstrations are performed with both the updated standalone BSE and the updated NBODY7/BSE. When convolved with stellar and primordial-binary populations as observed in young massive clusters, such remnant-formation and natal-kick mechanisms crucially determine the accumulated number, mass, and mass distribution of the BHs retained in young massive, open, and globular clusters (GCs); these BHs would eventually become available for long-term dynamical processing. Among other conclusions, we find that although the newer, delayed SN remnant formation model gives birth to the largest number (mass) of BHs, the older remnant-formation schemes cause the largest number (mass) of BHs to survive in clusters, when incorporating SN material fallback onto the BHs. The SN material fallback also causes the convection-asymmetry-driven SN kick to effectively retain similar numbers and masses of BHs in clusters as for the standard, momentum-conserving kick. The collapse-asymmetry-driven SN kick would cause nearly all BHs to be retained in clusters irrespective of their mass, remnant-formation model, and metallicity, whereas the inference of a large population of BHs in GCs would potentially rule out the neutrino-driven SN kick mechanism. Pre-SN mergers of massive primordial binaries would potentially cause BH masses to deviate from the theoretical, single-star ZAMS to mass-remnant mass relation unless a substantial of the total merging stellar mass of up to ≈40% is lost during a merger process. In particular, such mergers, at low metallicities, have the potential to produce low-spinning BHs within the PSN mass gap that can be retained in a stellar cluster and be available for subsequent dynamical interactions. As recent studies indicate, the new remnant-formation modelling reassures us that young massive and open clusters would potentially contribute to the dynamical BBH merger detection rate to a similar extent as their more massive GC counterparts.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Quantum magic and computational complexity in the neutrino sector

We consider the quantum magic in systems of dense neutrinos undergoing coherent flavor transformations, relevant for supernova and neutron-star binary mergers. Mapping the three-flavor-neutrino system to qutrits, the evolution of quantum magic is explored in the single scattering angle limit for a selection of initial tensor-product pure states for 𝑁 𝜈 ≤ 8 neutrinos. For |𝜈𝑒⟩ ⊗𝑁𝜈 initial states, the magic, as measured by the 𝛼 = 2 stabilizer Renyi entropy ℳ 2 , is found to decrease with radial distance from the neutrino sphere, reaching a value that lies below the maximum for tensor-product qutrit states. Further, the asymptotic magic per neutrino, ℳ 2 /𝑁 𝜈 , decreases with increasing 𝑁 𝜈 . In contrast, the magic evolving from states containing all three flavors reaches values only possible with entanglement, with the asymptotic ℳ 2 /𝑁 𝜈 increasing with 𝑁 𝜈 . These results highlight the connection between the complexity in simulating quantum physical systems and the parameters of the Standard Model.

computational complexity↗