Search NASA⌕ Search

SEARCH · Search NASA

Results for “Neutron star mergers”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Impact of Newly Measured 𝛽-Delayed Neutron Emitters around 78 Ni on Light Element Nucleosynthesis in the Neutrino Wind Following a Neutron Star Merger

Neutron emission probabilities and half-lives of 37 𝛽-delayed neutron emitters from 75 Ni to 92 Br were measured at the RIKEN Nishina Center in Japan, including 11 one-neutron and 13 two-neutron emission probabilities and six half-lives for the first time that supersede theoretical estimates. These nuclei lie in the path of the weak 𝑟 process occurring in neutrino-driven winds from the accretion disk formed after the merger of two neutron stars synthesizing elements in the 𝐴∼80 abundance peak. The presence of such elements dominates the accompanying kilonova emission over the first few days and have been identified in the AT2017gfo event, associated to the gravitational wave detection GW170817. Abundance calculations based on over 17,000 simulated trajectories describing the evolution of matter properties in the merger outflows show that the new data lead to an increase of 50%–70% in the abundance of Y, Zr, Nb, and Mo. This enhancement is large compared to the scatter of relative abundances observed in old very metal poor stars and thus is significant in the comparison with other possible astrophysical processes contributing to the light-element production. These results underline the importance of including experimental decay data for very neutron-rich 𝛽 -delayed neutron emitters into 𝑟 -process models.

59 ≤ A ≤ 89↗
DOE OSTI2025

Impact of Magnetic-field-driven Anisotropies on the Equation of State Probed in Neutron Star Mergers

Binary neutron star mergers can produce extreme magnetic fields, some of which can lead to strong magnetar-like remnants. While strong magnetic fields have been shown to affect the dynamics of outflows and angular momentum transport in the remnant, they can also crucially alter the properties of nuclear matter probed in the merger. In this work, we provide a first assessment of the latter, determining the strength of the pressure anisotropy caused by Landau-level quantization and the anomalous magnetic moment. To this end, we perform the first numerical relativity simulation with a magnetic polarization tensor and a magnetic-field-dependent equation of state using a new algorithm we present here, which also incorporates a mean-field dynamo model to control the magnetic field strength present in the merger remnant. Our results show that—in the most optimistic case—corrections to the anisotropy can be in excess of 10% and are potentially largest in the outer layers of the remnant. This work paves the way for a systematic investigation of these effects.

General relativity↗
NASA NTRS1993

Gamma-ray bursts from neutron star mergers

Binary neutron stars merger (NS(sup 2)M) at cosmological distances is probably the only gamma-ray bursts model based on an independently observed phenomenon which is known to be taking place at a comparable rate. We describe this model, its predictions and some open questions.

Piran, Tsvi↗
NASA NTRS2021

HARM3D+NUC: A New Method for Simulating the Post-merger Phase of Binary Neutron Star Mergers with GRMHD, Tabulated EOS, and Neutrino Leakage

The first binary neutron star merger has already been detected in gravitational waves. The signal was accompanied by an electromagnetic counterpart including a kilonova component powered by the decay of radioactive nuclei, as well as a short γ-ray burst. In order to understand the radioactively powered signal, it is necessary to simulate the outflows and their nucleosynthesis from the post-merger disk. Simulating the disk and predicting the composition of the outflows requires general relativistic magnetohydrodynamical (GRMHD) simulations that include a realistic, finite-temperature equation of state (EOS) and self-consistently calculating the impact of neutrinos. In this work, we detail the implementation of a finite-temperature EOS and the treatment of neutrinos in the GRMHD code HARM3D+NUC, based on HARM3D. We include formal tests of both the finite-temperature EOS and the neutrino-leakage scheme. We further test the code by showing that, given conditions similar to those of published remnant disks following neutron star mergers, it reproduces both recombination of free nucleons to a neutron-rich composition and excitation of a thermal wind.

Ariadna Murguia-Berthier↗
DOE OSTI2024

Turbulence modelling in neutron star merger simulations

Observations of neutron star mergers have the potential to unveil detailed physics of matter and gravity in regimes inaccessible by other experiments. Quantitative comparisons to theory and parameter estimation require nonlinear numerical simulations. However, the detailed physics of energy and momentum transfer between different scales, and the formation and interaction of small scale structures, which can be probed by detectors, are not captured by current simulations. This is where turbulence enters neutron star modelling. This review will outline the theory and current status of turbulence modelling for relativistic neutron star merger simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Connecting GRBs from binary neutron star mergers to nuclear properties of neutron stars

The fate of the binary neutron star (NS) merger remnants hinges sensitively upon the NS equation of state and the threshold mass, M ls , that separates a long-lived from a short-lived NS remnant. The nature of the electromagnetic counterparts is also influenced by the remnant type, particularly in determining whether a gamma-ray burst from a compact binary merger (cbGRB) is of short or long duration. We propose a novel approach to probe Mls by linking it to the estimated observed ratio of long to short cbGRBs. We find that current observations broadly favor a relatively high value for this transition, M ls ≃ 1.3M TOV , for which M TOV ≲ 2.6M ⊙ , consistent with numerical simulations, as also shown here. Furthermore, our results disfavor nuclear physics scenarios that would lead to catastrophic pressure loss at a few times nuclear density and temperatures of tens of MeV, leading to a rapid gravitational collapse of binaries with total mass M ≲ 1.3M TOV . Future individual gravitational wave events with on-axis cbGRBs can further bound Mls.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Fragmentation in Gravitationally Unstable Collapsar Disks and Subsolar Neutron Star Mergers

Abstract Although stable neutron stars (NSs) can in principle exist down to massesM ns ≈ 0.1M ⊙ , standard models of stellar core-collapse predict a robust lower limitM ns ≳ 1.2M ⊙ , roughly commensurate with the Chandrasekhar massM Ch of the progenitor’s iron core (electron fractionY e ≈ 0.5). However, this limit may be circumvented in sufficiently dense neutron-rich environments (Y e < 0.5) for which M Ch ∝ Y e 2 is reduced to ≲1M ⊙ . Such physical conditions could arise in the black hole accretion disks formed from the collapse of rapidly rotating stars (“collapsars”), as a result of gravitational instabilities and cooling-induced fragmentation, similar to models for planet formation in protostellar disks. We confirm that the conditions to form subsolar-mass NS (ssNS) may be marginally satisfied in the outer regions of massive neutrino-cooled collapsar disks. If the disk fragments into multiple ssNSs, their subsequent coalescence offers a channel for precipitating subsolar mass LIGO/Virgo gravitational-wave mergers that does not implicate primordial black holes. The model makes several additional predictions: (1) ∼Hz frequency Doppler modulation of the ssNS-merger gravitational-wave signals due to the binary’s orbital motion in the disk; (2) at least one additional gravitational-wave event (coincident within ≲hours), from the coalescence of the ssNS-merger remnant(s) with the central black hole; (3) an associated gamma-ray burst and supernova counterpart, the latter boosted in energy and enriched withr-process elements from the NS merger(s) embedded within the exploding stellar envelope (“kilonovae inside a supernova”).

Astronomy & Astrophysics↗
DOE OSTI2025

Nuclear Astrophysics through simulations of neutron star mergers using Monte‐Carlo neutrino radiation transport (DE-SC0020435 Final Technical Report)

Neutron star mergers are an important source of information for nuclear physics. Except for black holes, neutron stars are the densest macroscopic objects known to exist in the Universe. They provide us with a remarkable laboratory to study the poorly understood high-density states of matter, complementing the study of heavy nuclei performed by nuclear physics experiments on Earth. Collisions of neutron stars are of particular interest. They are among the most energetic events observable in the Universe, powering a broad range of signals across the electromagnetic spectrum as well as gravitational wave signals. These signals contain important information about the properties of extremely dense matter. Colliding neutron stars additionally eject large amounts of neutron-rich material into the surrounding interstellar medium -- material that then undergoes rapid neutron-capture (r-process) nucleosynthesis, the mechanism thought to be responsible for the production of about half of the heavy nuclei. In order to understand current and future observations of neutron star mergers, we need reliable models for the signals that they power. An important component in the construction of such models are numerical simulations of colliding neutron stars. These simulations are costly, running for multiple months on supercomputers, and require the inclusion of complex physics (general relativity, magnetohydrodynamics, neutrino physics, nuclear reactions). In this document, we describe the development of new methods for the treatment of neutrinos in merger simulations, as well as the first simulations capable of evolving the equations of neutrino radiation transport directly and their impact on our understanding of neutron star mergers and more broadly in nuclear astrophysics. We also discuss additional work performed in the study of nucleosynthesis and neutrino physics in neutron star mergers as part of Early Career Award DE-SC0020435.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Bound Debris Expulsion from Neutron Star Merger Remnants

Abstract 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 ofr-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.

Astronomy & Astrophysics↗
DOE OSTI2023

Neutrino fast flavor instability in three dimensions for a neutron star merger

The flavor evolution of neutrinos in core collapse supernovae and neutron star mergers is a critically important unsolved problem in astrophysics. Following the electron flavor evolution of the neutrino system is essential for calculating the thermodynamics of compact objects as well as the chemical elements they produce. Accurately accounting for flavor transformation in these environments is challenging for a number of reasons, including the large number of neutrinos involved, the small spatial scale of the oscillation, and the nonlinearity of the system. We take a step in addressing these issues by presenting a method which describes the neutrino fields in terms of angular moments. We apply our moment method to neutron star merger conditions and show it simulates fast flavor neutrino transformation in a region where this phenomenon is expected to occur. By comparing with particle-in-cell calculations we show that the moment method is able to capture the three phases of growth, saturation, and decoherence, and correctly predicts the lengthscale of the fastest growing fluctuations in the neutrino field.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Impact of in situ nuclear networks and atomic opacities on neutron star merger ejecta dynamics, nucleosynthesis, and kilonovae

Context. Binary neutron star merger (BNSM) ejecta are key sites of rapid neutron capture (r-process) nucleosynthesis and they produce kilonovae powered by the radioactive decay of freshly synthesized nuclei. Modeling their evolution requires multi-physics simulations involving hydrodynamics, nuclear reactions, and radiative processes. The impact of nuclear burning and atomic opacity is poorly understood and often treated with simplified prescriptions. Aims. We systematically investigate different treatments of nuclear heating, particle thermalization, and atomic opacities in radiation-hydrodynamics simulations of BNSM ejecta and kilonova light curves. Methods. Ejecta profiles from long-term numerical-relativity simulations of asymmetric neutron star binaries with a massive neutron star remnant were evolved to ∼30 days using a 2D ray-by-ray approach. We compared simplified heating-rate and thermalization prescriptions with in situ Nuclear reaction Network (NN) calculations that track nuclear energy deposition and include a composition-dependent thermalization scheme. We also contrasted various gray opacity models with a frequency-dependent treatment based on atomic calculations. Results. Coupling NN and hydrodynamics significantly affects nucleosynthesis and kilonova emission. Assuming homologous expansion alters abundance evolution and produces a narrower, less populated second r-process peak and a third peak shifted to higher mass numbers. The back-reaction of nuclear heating affects the temperature evolution enough to delay and redden the early (t∼ hours) kilonova peaks. A constant thermalization efficiency underestimates and reddens the early emission while overestimating the late-time luminosity compared to the composition-dependent treatment. Analytical opacity prescriptions yield a more extended, colder photosphere, resulting in dimmer, redder kilonovae at early times (t≲ hour), while the delayed recession of the photosphere prolongs the red emission at t ≳ 5 days. Conclusions. Coupling hydrodynamics to an in situ NN is crucial for reliable nucleosynthesis and kilonova predictions. Resolving the first several hundred milliseconds of the hydrodynamics is essential for robust nucleosynthesis calculations. Composition-dependent thermalization and frequency-dependent, atomic-physics-based opacities are needed to accurately capture the temperature evolution of the ejecta and the brightness and color evolution of the kilonova. Calibrated analytic nuclear-power fits with simplified thermalization and opacity prescriptions can still reproduce the density and temperature evolution of the ejecta.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Neutrino fast flavor oscillations with moments: Linear stability analysis and application to neutron star mergers

Providing an accurate modeling of neutrino physics in dense astrophysical environments such as binary neutron star mergers presents a challenge for hydrodynamic simulations. Nevertheless, understanding how flavor transformation can occur and affect the dynamics, the mass ejection, and the nucleosynthesis will need to be achieved in the future. Computationally expensive, large-scale simulations frequently evolve the first classical angular moments of the neutrino distributions. By promoting these quantities to matrices in flavor space, we develop a linear stability analysis of fast flavor oscillations using only the first two “quantum” moments, which notably requires generalizing the classical closure relations that appropriately truncate the hierarchy of moment equations in order to treat quantum flavor coherence. After showing the efficiency of this method on a well-understood test situation, we perform a systematic search of the occurrence of fast flavor instabilities in a neutron star merger simulation. Here, we discuss the successes and shortcomings of moment linear stability analysis, as this framework provides a time-efficient way to design and study better closure prescriptions in the future.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Secular Outflows from Long-Lived Neutron Star Merger Remnants

We study mass ejection from a binary neutron star merger producing a long-lived massive neutron star remnant with general-relativistic neutrino-radiation hydrodynamics simulations. In addition to outflows generated by shocks and tidal torques during and shortly after the merger, we observe the appearance of a wind driven by spiral density waves in the disk. This spiral-wave-driven outflow is predominantly located close to the disk orbital plane and have a broad distribution of electron fractions. At higher latitudes, a high electron-fraction wind is driven by neutrino radiation. The combined nucleosynthesis yields from all the ejecta components is in good agreement with Solar abundance measurements.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Numerical relativity estimates of the remnant recoil velocity in binary neutron star mergers

We present, for the first time, recoil velocity estimates for binary neutron star mergers using data from numerical relativity simulations. We find that binary neutron star merger remnants can have recoil velocity of the order of a few tens of kilometers per second and as high as 150 kilometers per second in our dataset. These recoils are attained due to equivalent contributions from the anisotropic gravitational wave emission as well as the asymmetric ejection of dynamical matter during the merger. Here, we provide fits for net recoil velocity as well as its ejecta component as a function of the amount of ejected matter, which may be useful when constraints on the ejected matter are obtained through electromagnetic observations. We also estimate the mass and spin of the remnants and find them to be in the range [2.34,3.38]M ⊙ and [0.63, 0.82], respectively, for our dataset.

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↗
NASA NTRS2019

GROWTH on S190510g: DECam Observation Planning and Follow-up of a Distant Binary Neutron Star Merger Candidate

The first two months of the third Advanced LIGO and Virgo observing run (2019 April–May) showed that distant gravitational-wave (GW) events can now be readily detected. Three candidate mergers containing neutron stars (NS) were reported in a span of 15 days, all likely located more than 100 Mpc away. However, distant events such as the three new NS mergers are likely to be coarsely localized, which highlights the importance of facilities and scheduling systems that enable deep observations over hundreds to thousands of square degrees to detect the electromagnetic counterparts. On 2019 May 10 02:59:39.292 UT the GW candidate S190510g was discovered and initially classified as a binary neutron star (BNS) merger with 98% probability. The GW event was localized within an area of 3462 deg2, later refined to 1166 deg2 (90%) at a distance of 227 ± 92 Mpc. We triggered Target-of-Opportunity observations with the Dark Energy Camera (DECam), a wide-field optical imager mounted at the prime focus of the 4 m Blanco Telescope at Cerro Tololo Inter-American Observatory in Chile. This Letter describes our DECam observations and our real-time analysis results, focusing in particular on the design and implementation of the observing strategy. Within 24 hr of the merger time, we observed 65% of the total enclosed probability of the final skymap with an observing efficiency of 94%. We identified and publicly announced 13 candidate counterparts. S190510g was reclassified 1.7 days after the merger, after our observations were completed, with a "BNS merger" probability reduced from 98% to 42% in favor of a "terrestrial classification.

Igor Andreoni↗
DOE OSTI2024

First Constraints on the Photon Coupling of Axionlike Particles from Multimessenger Studies of the Neutron Star Merger GW170817

We use multimessenger observations of the neutron star merger event GW170817 to derive new constraints on axionlike particles (ALPs) coupling to photons. ALPs are produced via Primakoff and photon coalescence processes in the merger, escape the remnant, and decay back into two photons, giving rise to a photon signal approximately along the line of sight to the merger. We analyze the spectral and temporal information of the ALP-induced photon signal and use the Fermi Large Area Telescope (Fermi-LAT) observations of GW170817 to derive our new ALP constraints. We also show the improved prospects with future MeV 𝛾-ray missions, taking the spectral and temporal coverage of Fermi-LAT as an example.

Physics↗
DOE OSTI2024

Neutrino Trapping and Out-of-Equilibrium Effects in Binary Neutron-Star Merger Remnants

We study out-of-thermodynamic-equilibrium effects in neutron-star mergers with 3D general-relativistic neutrino-radiation large-eddy simulations. During mergers, the cores of the neutron stars remain cold (T ~ a few MeV) and out of thermodynamic equilibrium with trapped neutrinos originating from the hot collisional interface between the stars. However, within ~2 to 3 ms matter and neutrinos reach equilibrium everywhere in the remnant massive neutron star. Furthermore, our results show that dissipative effects, such as bulk viscosity, if present, are only active for a short window of time after the merger.

79 ASTRONOMY AND ASTROPHYSICS↗