Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

General relativistic hydrodynamic simulations of binary strange star mergers

We perform fully general-relativistic simulations of binary strange star mergers considering two different approaches for thermal effects. The first uses a cold equation of state (EOS) derived from a modified version of the MIT bag model which is then supplemented by a Γ-law correction. The second approach employs a microphysical description of the finite-temperature effects. We describe results obtained with the two treatments, highlighting the influence of thermal effects. We find that the postmerger dynamics differs significantly in the two cases, leading to quantitative differences in the postmerger gravitational-wave spectrum and ejecta mass. The peak frequency of the postmerger gravitational-wave emission is consistent with the established quasi-universal relations for binary neutron star mergers and as a result, our simulations cannot distinguish between mergers of neutron stars and those of strange stars. Our models with realistic treatment of finite-temperature effects produce a significant amount of ejecta ≳0.02 M ⊙ ​. Here, the resulting flux of strangelets near the Earth, computed assuming that all neutron star mergers are in fact strange-stars mergers and that the binary considered here is representative, is in tension with experimental upper limits. As such, our results tentatively disfavor a scenario in which strange-quark matter is the lowest energy state of matter.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Exciting prospects for detecting late-time neutrinos from core-collapse supernovae

The importance of detecting neutrinos from a Milky Way core-collapse supernova is well known. An understudied phase is proto-neutron star cooling. For SN 1987A, this seemingly began at about 2 s and is thus probed by only 6 of the 19 events (and only the $\bar{ν}_e$ flavor) in the Kamiokande-II and IMB detectors. With the higher statistics expected for present and near-future detectors, it should be possible to measure detailed neutrino signals out to very late times. Here, we present the first comprehensive study of neutrino detection during the proto-neutron star cooling phase, considering a variety of outcomes, using all flavors, and employing detailed detector physics. For our nominal model, the event yields (at 10 kpc) after 10 s—the approximate duration of the SN 1987A signal—far exceed the entire SN 1987A yield, with ≃ 250 $\bar{ν}_e$ events (to 50 s) in Super-Kamiokande, ≃ 110 ν e events (40 s) in the Deep Underground Neutrino Experiment (DUNE), and ≃10 $ν_μ, ν_τ, \bar{ν}_μ, \bar{ν}_τ$ events (to 20 s) in the Jiangmen Underground Neutrino Observatory. These data would allow unprecedented probes of the proto-neutron star, including the onset of neutrino transparency and hence its transition to a neutron star. If a black hole forms, even at very late times, this can be clearly identified. But will the detectors fulfill their potential for this perhaps once-ever opportunity for an all-flavor, high-statistics detection of a core collapse? Maybe. Further work is urgently needed, especially for DUNE to thoroughly investigate and improve its MeV capabilities.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Impact of the nuclear equation of state on the formation of twin stars

Twin stars-two stable neutron stars (NSs) with the same mass but different radii have long been proposed to appear as a consequence of a possible first-order phase transition in NS matter. Within a meta-model for the EOS of hybrid stars, we revisit the viability of twin stars and its dependence on numerous parameters characterizing the EOS of nuclear matter, quark matter, and the phase transition between them. While essentially no experimental constraint exists for the last two, parameters characterizing the EOS of neutron-rich nucleonic matter have been constrained within various ranges by terrestrial experiments and astrophysical observations. Within these ranges, the impact of nuclear EOS and crust-core transition density on the formation of twin stars is studied. Here, it is found that the symmetry energy of neutron-rich nucleonic matter notably influences the formation of twin stars, particularly through its slope L and curvature K sym . Conversely, varying the EOS of symmetric nuclear matter within their currently known uncertainty ranges shows minimal influence on the formation of twin stars.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Framework for phase transitions between the Maxwell and Gibbs constructions

By taking the nucleon-to-quark phase transition within a neutron star as an example, we present a thermodynamically consistent method to calculate the equation of state of ambient matter so that transitions that are intermediate to those of the familiar Maxwell and Gibbs constructions can be described. This method does not address the poorly known surface tension between the two phases microscopically (as, for example, in the calculation of the core pasta phases via the Wigner-Seitz approximation) but instead combines the local and global charge neutrality conditions characteristic of the Maxwell and Gibbs constructions, respectively. Overall charge neutrality is achieved by dividing the leptons to those that obey local charge neutrality (Maxwell) and those that maintain global charge neutrality (Gibbs). The equation of state is obtained by using equilibrium constraints derived from minimizing the total energy density. The results of this minimization are then used to calculate neutron star mass-radius curves, tidal deformabilities, equilibrium and adiabatic sound speeds, and nonradial g-mode oscillation frequencies for several intermediate constructions. Various quantities of interest transform smoothly from their Gibbs structures to those of Maxwell as the local-to-total electron ratio η, introduced to mimic the hadron-to-quark interface tension from 0 (Gibbs) to ∞ (Maxwell), is raised from 0 to 1. As a result, a notable exception is the g-mode frequency for the specific case of η = 1 for which a gap appears between the quark and hadronic branches.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

GW190425: Observation of a Compact Binary Coalescence with Total Mass ~ 3.4 M ⊙

On 2019 April 25, the LIGO Livingston detector observed a compact binary coalescence with signal-to-noise ratio 12.9. The Virgo detector was also taking data that did not contribute to detection due to a low signal-to-noise ratio, but were used for subsequent parameter estimation. The 90% credible intervals for the component masses range from 1.12 to 2.52 M ⊙ (1.46-1.87 M ⊙ if we restrict the dimensionless component spin magnitudes to be smaller than 0.05). These mass parameters are consistent with the individual binary components being neutron stars. However, both the source-frame chirp mass 1.44 +0.02 -0.02 M ⊙ and the total mass 3.4 +0.3 -0.1 M ⊙ of this system are significantly larger than those of any other known binary neutron star (BNS) system. The possibility that one or both binary components of the system are black holes cannot be ruled out from gravitational-wave data. We discuss possible origins of the system based on its inconsistency with the known Galactic BNS population. Under the assumption that the signal was produced by a BNS coalescence, the local rate of neutron star mergers is updated to 250-2810 Gpc -3 yr -1 .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Towards Exascale Astrophysics of Mergers and Supernovae (TEAMS)

The TEAMS project brought together cutting-edge simulations, theoretical insights, and collaborative efforts to deepen our understanding of some of the universe’s most extreme phenomena—supernovae, neutron star mergers, and the powerful signals they emit. Using one of the largest suites of 3D supernova simulations ever conducted, researchers uncovered new insights into how massive stars explode, how those explosions vary by stellar mass, and what conditions lead to the birth of neutron stars or black holes. They also studied the radiation and gravitational wave signals emitted during these events, revealing how future observations can be used to uncover what happens deep inside collapsing stars. The team developed improved tools for modeling how light and neutrinos behave in such explosive environments, enabling more accurate predictions of what astronomers might observe. Work also explored how the chemical composition and geometry of kilonovae—the visible explosions that follow neutron star mergers—influence their signals and can reveal the origins of heavy elements like gold. These efforts not only advanced scientific knowledge, but also trained a new generation of researchers at the intersection of astrophysics, computational science, and nuclear theory.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Modeling Kilonova Light Curves: Dependence on Nuclear Inputs

The mergers of binary neutron stars, as well as black hole–neutron star systems, are expected to produce an electromagnetic counterpart that can be analyzed to infer the element synthesis that occurred in these events. We investigate one source of uncertainties pertinent to lanthanide-rich outflows: the nuclear inputs to rapid neutron capture nucleosynthesis calculations. We begin by examining 32 different combinations of nuclear inputs: eight mass models, two types of spontaneous fission rates, and two types of fission daughter product distributions. We find that such nuclear physics uncertainties typically generate at least one order of magnitude uncertainty in key quantities such as the nuclear heating (one and a half orders of magnitude at 1 day post-merger), the bolometric luminosity (one order of magnitude at 5 days post-merger), and the inferred mass of material from the bolometric luminosity (factor of 8 when considering the 8–10 day region). Since particular nuclear processes are critical for determining the electromagnetic signal, we provide tables of key nuclei undergoing β-decay, α-decay, and spontaneous fission important for heating at different times, identifying decays that are common among the many nuclear input combinations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Validity of a finite temperature expansion for dense nuclear matter

In this work we provide a new, well-controlled expansion of the equation of state of dense matter from zero to finite temperatures (𝑇) while covering a wide range of charge fractions (𝑌 𝑄 ), from pure neutron to isospin symmetric nuclear matter. Our expansion can be used to describe neutron star mergers using the equation of state inferred from neutron star observations. We discuss how knowledge from low-energy nuclear experiments and heavy-ion collisions can be directly incorporated into the expansion. We also suggest new thermodynamic quantities of interest that can be calculated from theoretical models or directly inferred by experimental data that can be used to infer the finite temperature equation of state. With our new method, we can quantify the uncertainty in our finite 𝑇 and 𝑌 𝑄 expansions without making assumptions about the underlying degrees of freedom. We can reproduce results from a microscopic equation of state up to 𝑇 = 100 MeV for baryon chemical potential 𝜇 𝐵 ≳ 1100 MeV [≈(1–2)⁢𝑛 sat ] within 5% error, with even better results for larger 𝜇 𝐵 and/or lower 𝑇. We investigate the sources of numerical and theoretical uncertainty and discuss future directions of study.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Gamma Rays from Kilonova: A Potential Probe of r-process Nucleosynthesis

The mergers of compact binaries with at least one neutron star component are the potential leading sites of the production and ejection of r-process elements. Discoveries of galactic binary pulsars, short gamma-ray bursts, and gravitational-wave detections have all been constraining the rate of these events, while the gravitational wave plus broadband electromagnetic coverage of binary neutron star merger (GW170817) has also placed constraints on the properties (mass and composition) of the merger ejecta. But uncertainties and ambiguities in modeling the optical and infrared emission make it difficult to definitively measure the distribution of heavy isotopes in these mergers. In contrast, gamma rays emitted in the decay of these neutron-rich ejecta may provide a more direct measurement of the yields. We calculate the gamma production in remnants of neutron star mergers, considering two epochs: a kilonova epoch, lasting about two weeks, and a much later epoch of tens and hundreds of thousands of years after the merger. For the kilonova epoch, when the expanding ejecta is still only partially transparent to gamma radiation, we use 3D radiative transport simulations to produce the spectra. Futhermore, we show that the gamma-ray spectra associated with beta- and alpha-decay provide a fingerprint of the ejecta properties and, for a sufficiently nearby remnant, may be detectable, even for old remnants. We compare our gamma spectra with the potential detection limits of next generation detectors, including the Lunar Occultation Explorer (LOX), the All-sky Medium Energy Gamma-ray Observatory (AMEGO), and the Compton Spectrometer and Imager (COSI). We show that fission models can be discriminated via the presence of short-lived fission fragments in the remnant spectra.

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 OSTI2021

Bulk viscosity from Urca processes: npeμ matter in the neutrino-trapped regime

In this work, we extend our previous study of the bulk viscosity of hot and dense npe matter induced by the Urca process in the neutrino trapped regime to npeμ matter by adding the muonic Urca processes as well as the purely leptonic electroweak processes involving electron-muon transition. Here, the nuclear matter is modeled in a relativistic density functional approach with two different parametrizations which predict neutrino dominated matter (DDME2 model) and anti-neutrino dominated matter (NL3 model) at temperatures for which neutrinos/anti-neutrinos are trapped. In the case of neutrino-dominated matter, the main equilibration mechanism is lepton capture, whereas in the case of antineutrino-dominated matter this is due to neutron decay. We find that the equilibration rates of Urca processes are higher than that of the pure leptonic processes, which implies that the Urca bulk viscosity can be computed with the leptonic reactions assumed to be frozen. We find that the bulk viscosity decreases with temperature as ζ~T –2 at moderate temperatures. At high temperatures this scaling breaks down by sharp drops of the bulk viscosity close to the temperature where the proton fraction is density-independent and the matter becomes scale-invariant. This occurs also when the matter undergoes a transition from the antineutrino-dominated regime to the neutrino-dominated regime where the bulk viscosity attains a local maximum. We also estimate the bulk viscous dissipation timescales and find that these are in the range ≳ 1 s for temperatures above the neutrino trapping temperature. These timescales would be relevant only for long-lived objects formed in binary neutron star mergers and hot proto-neutron stars formed in core-collapse supernovas.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Excluded-volume model for quarkyonic matter. II. Three-flavor shell-like distribution of baryons in phase space

We extend the excluded-volume model of isospin symmetric two-flavor dense quarkyonic matter [Phys. Rev. C 101, 035201 (2020)] including strange particles and address its implications for neutron stars. The effective sizes of baryons are defined from the diverging hard-core potentials in the short interdistance regime. Around the hard-core density, the repulsive core between baryons at short distances leads to a saturation in the number density of baryons and generates perturbative quarks from the lower phase space, which leads to the shell-like distribution of baryons by the Pauli exclusion principle. The strange-quark Fermi sea always appears at high densities but the Λ hyperon shell only appears when the effective size of the Λ hyperon is smaller than the effective size of nucleons. We find that the pressure of strange quarkyonic matter can be large enough to support neutron stars with two times solar mass and can have a large sound speed, $c^2_s ≃ 0.7$. Finally, the fraction of the baryon number carried by perturbative quarks is about 30% at the inner core of most massive neutron stars.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Constraints on Nuclear Symmetry Energy Parameters

A review is made of constraints on the nuclear symmetry energy parameters arising from nuclear binding energy measurements, theoretical chiral effective field predictions of neutron matter properties, the unitary gas conjecture, and measurements of neutron skin thicknesses and dipole polarizabilities. While most studies have been confined to the parameters 𝑆 𝑉 and L, the important roles played by, and constraints on 𝐾 sym , or, equivalently, the neutron matter incompressibility 𝐾 𝑁 , are discussed. Strong correlations among 𝑆 𝑉 ,𝐿, and 𝐾 𝑁 are found from both nuclear binding energies and neutron matter theory. However, these correlations somewhat differ in the two cases, and those from neutron matter theory have smaller uncertainties. To 68% confidence, it is found from neutron matter theory that 𝑆 𝑉 = 32.0 ± 1.1 MeV, 𝐿 = 51.9 ± 7.9 MeV and 𝐾 𝑁 = 152.2 ± 38.1 MeV. Theoretical predictions for neutron skin thickness and dipole polarizability measurements of the neutron-rich nuclei 48 Ca, 120 Sn, and 208 Pb are compared to recent experimental measurements, most notably the CREX and PREX neutron skin experiments from Jefferson Laboratory. By themselves, PREX I+II measurements of 208 Pb and CREX measurement of 48 Ca suggest 𝐿 = 121 ± 47 MeV and 𝐿 = −5 ± 40 MeV, respectively, to 68% confidence. However, we show that nuclear interactions optimally satisfying both measurements imply 𝐿 = 53 ± 13 MeV, nearly the range suggested by either nuclear mass measurements or neutron matter theory, and is also consistent with nuclear dipole polarizability measurements. This small parameter range implies 𝑅 1.4 = 11.6 ± 1.0 km and Λ 1.4 = 228$^{+148}_{−90}$, which are consistent with NICER X-ray and LIGO/Virgo gravitational wave observations of neutron stars.

neutron skins↗
DOE OSTI2024

Astrophysical Equation-of-State Constraints on the Color-Superconducting Gap

We demonstrate that astrophysical constraints on the dense-matter equation of state place an upper bound on the color-superconducting gap in dense matter above the transition from nuclear matter to quark matter. Pairing effects in the color-flavor locked quark matter phase increase the pressure at high density, and if this effect is sufficiently large then the requirements of causality and mechanical stability make it impossible to reach such a pressure in a way that is consistent with what is known at lower densities. The intermediate-density equation of state is inferred by considering extensions of chiral effective field theory to neutron star densities, and conditioning these using current astrophysical observations of neutron star radius, maximum mass, and tidal deformability (PSR J⁢0348+0432, PSR J1624-2230, PSR J⁢0740+6620, GW170817). At baryon number chemical potential μ = 2.6 GeV we find a 95% upper limit on the color-flavor locked pairing gap Δ of 457 MeV using overly conservative assumptions and 216 MeV with more reasonable assumptions. Furthermore, this constraint may be strengthened by future astrophysical measurements as well as by future advances in high-density QCD calculations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Axion-photon conversion in strongly magnetised plasmas

Axion dark matter can resonantly convert to photons in the magnetosphere of neutron stars, possibly giving rise to radio signals observable on Earth. This method for the indirect detection of axion dark matter has recently received significant attention in the literature. The calculation of the radio signal is complicated by a number of effects; most importantly, the gravitational infall of the axions onto the neutron star accelerates them to semi-relativistic speed, and the neutron star magnetosphere is highly anisotropic. Both of these factors complicate the calculation of the conversion of axions to photons. In this work, we present the first fully three-dimensional calculation of the axion-photon conversion in highly magnetised anisotropic media. Depending on the axion trajectory, this calculation leads to orders-of-magnitude differences in the conversion compared to the simplified one-dimensional calculation used so far in the literature, altering the directionality of the produced photons. Our results will have important implications for the radio signal one would observe in a telescope.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Nonparametric extensions of nuclear equations of state: Probing the breakdown scale of relativistic mean-field theory

Phenomenological calculations of the properties of dense matter, such as relativistic mean-field theories, represent a pathway to predicting the microscopic and macroscopic properties of neutron stars. However, such theories do not generically have well-controlled uncertainties and may break down within neutron stars. To faithfully represent the uncertainty in this breakdown scale, we develop a hybrid representation of the dense-matter equation of state, which assumes the form of a relativistic mean-field theory at low densities, while remaining agnostic to any nuclear theory at high densities. To achieve this, we use a nonparametric equation of state model to incorporate the correlations of the underlying relativistic mean-field theory equation of state at low pressures and transition to more flexible correlations above some chosen pressure scale. We perform astrophysical inference under various choices of the transition pressure between the theory-informed and theory-agnostic models. Here, we further study whether the chosen relativistic mean-field theory breaks down above some particular pressure and find no such evidence. Using simulated data for future astrophysical observations at about two-to-three times the precision of current constraints, we show that our method can identify the breakdown pressure associated with a potential strong phase transition.

Equations of state of nuclear matter↗
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 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↗