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

Simulating Neutron Stars with Solid Quark Cores: Rotations, Oscillations, and Binary Mergers [Slides]

Neutron stars feature complexity from GR, QM, particle physics and nuclear physics effects which cannot be investigated in the lab up to the densities in their core. GWs from binary neutron star mergers are on the rise to provide more insight into the physics inside NSs. There is a possibility that involves crystalline, color-superconducting quark matter, with implications towards altering the GW signal if NSs harbor this exotic matter. There are no 3D dynamical simulations of solid quark matter inside Neutron stars. We perform the first 3D dynamical simulations of solid quark matter inside NSs. These are performed in the Newtonian Approximation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Signature of hadron-quark crossover in binary-neutron-star mergers

Here, we study observational signatures of the hadron-quark crossover in binary-neutron-star mergers by numerical-relativity simulations with various mass configurations. We employ two equations of state (EoSs) for matter consistent with inference from the observational data. In the crossover scenario the EoS is softened in a density realized in binary-neutron-star mergers and is smoothly continued to quark matter. In the phase transition scenario without crossover, the EoS remains stiff and a first-order phase transition takes place in a density out of reach of mergers. A GW170817-like system forms a remnant massive neutron star in both scenarios, and it collapses into a black hole only in the crossover scenario due to the softening while gravitational-wave emission is strong. This difference is clearly reflected in the sudden shutdown of gravitational waves. For a given EoS, the lifetime of the merger remnant is determined primarily by the total mass of the system. Identifying these features in a variety of future events with the next generation of ground-based gravitational-wave detectors will enable us to clarify details of hadron-quark transition. The mass of the accretion disk surrounding the remnant black hole is affected not only by the lifetime of the remnant but also by the mass ratio of the system. Electromagnetic emission associated with the disk outflow will also be useful for detailed investigation of the hadron-quark transition.

QCD phase transitions↗
DOE OSTI2024

Measuring neutron star radius with second and third generation gravitational wave detector networks

The next generation of ground-based interferometric gravitational wave detectors will observe mergers of black holes and neutron stars throughout cosmic time. A large number of the binary neutron star merger events will be observed with extreme high fidelity, and will provide stringent constraints on the equation of state of nuclear matter. In this paper, we investigate the systematic improvement in the measurability of the equation of state with increase in detector sensitivity by combining constraints obtained on the radius of a 1.4 M ⊙ neutron star from a simulated source population. Since the measurability of the equation of state depends on its stiffness, we consider a range of realistic equations of state that span the current observational constraints. We show that a single 40 km Cosmic Explorer detector can pin down the neutron star radius for a soft, medium and stiff equation of state with a precision of 10 m within a decade, whereas the current generation of ground-based detectors like the Advanced LIGO-Virgo network would take $\mathcal{O}(10_5)$ years to do so for a soft equation of state.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Limiting Rotation Rate of Neutron Stars from Crust Breaking and Gravitational Waves

Abstract Neutron stars are not observed to spin faster than about half their breakup rate. This limiting rotational frequency may be related to the strength of their crusts. As a star spins up from accretion, centrifugal forces stress the crust. We perform finite-element simulations of rotating neutron stars and find that the crust fails at rotation rates about half the breakup rate. Given uncertainties in microphysics, we have not determined the crust configuration after this failure. Instead, we argue that the crust may fail in an asymmetric way and could produce a configuration with a significant ellipticity (fractional difference in moments of inertia). If the ellipticity is large, a rotating star will radiate gravitational waves that may limit further spin-up. These stars may be promising sources for LIGO/VIRGO and next-generation gravitational-wave detectors.

Morales, J. A.↗
DOE OSTI2024

Modeling Competing Line-broadening Mechanisms in Neutron Star Atmospheres: Interference between Motional Stark and Ion Broadening

Neutron star surfaces have extremely high magnetic fields. In the atmosphere, the broadening of spectral lines will be substantial from the dense plasma as well as from the magnetic field. One broadening mechanism of note is due to the motional Stark effect (MSE)—an additional electric field that arises from the motion of the atom in the magnetic field. However, approximate formulae are often used to construct atmosphere models, and the MSE is assumed to be the dominant line-broadening mechanism even in ions. Detailed pressure-broadening models in these extreme magnetic fields are now currently being developed. In these more detailed models, it was suggested that the MSE may not be as large as previously predicted. If correct, this hypothesis implies that neutron star line widths might be dominated by pressure broadening rather than by motional Stark broadening. We find that, in the absence of plasma perturbations, for typical magnetic fields (B = 10 12 G), mid-Z elements, such as oxygen, have motional Stark widths of order 1 eV for transitions between dipole-allowed transitions from the ground state, though higher temperatures and transitions to higher-energy states are expected to have more broadening. The MSE also breaks down selection rules, giving rise to forbidden transitions, which have much larger widths. When plasma perturbations are included, we find that the plasma perturbation and motional Stark processes are not independent and, as a result, the spectral lines become narrow in a nontrivial way and display harmonics of the ion cyclotron frequency.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Combining Electromagnetic and Gravitational-Wave Constraints on Neutron-Star Masses and Radii

In this work, we perform a joint Bayesian inference of neutron-star mass and radius constraints based on GW170817, observations of quiescent low-mass x-ray binaries (QLMXBs), photospheric radius expansion x-ray bursting sources, and x-ray timing observations of J0030+0451. With this dataset, the form of the prior distribution still has an impact on the posterior mass-radius curves and equation of state (EOS), but this impact is smaller than recently obtained when considering QLMXBs alone. We analyze the consistency of the electromagnetic data by including an “intrinsic scattering” contribution to the uncertainties, and find only a slight broadening of the posteriors. This suggests that the gravitational-wave and electromagnetic observations of neutron-star structure are providing a consistent picture of the neutron-star mass-radius curve and the EOS.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Impact of Pycnonuclear Fusion Uncertainties on the Cooling of Accreting Neutron Star Crusts

Abstract The observation of X-rays during quiescence from transiently accreting neutron stars provides unique clues about the nature of dense matter. This, however, requires extensive modeling of the crusts and matching the results to observations. The pycnonuclear fusion reaction rates implemented in these models are theoretically calculated by extending phenomenological expressions and have large uncertainties spanning many orders of magnitude. We present the first sensitivity studies of these pycnonuclear fusion reactions in realistic network calculations. We also couple the reaction network with the thermal evolution code dStar to further study their impact on the neutron star cooling curves in quiescence. Varying the pycnonuclear fusion reaction rates alters the depth at which nuclear heat is deposited although the total heating remains constant. The enhancement of the pycnonuclear fusion reaction rates leads to an overall shallower deposition of nuclear heat. The impurity factors are also altered depending on the type of ashes deposited on the crust. These total changes correspond to a variation of up to 9 eV in the modeled cooling curves. While this is not sufficient to explain the shallow heat source, it is comparable to the observational uncertainties and can still be important for modeling the neutron star crust.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

A Theory for Neutron Star and Black Hole Kicks and Induced Spins

Abstract Using 20 long-term 3D core-collapse supernova simulations, we find that lower compactness progenitors that explode quasi-spherically due to the short delay to explosion experience smaller neutron star recoil kicks in the ∼100−200 km s −1 range, while higher compactness progenitors that explode later and more aspherically leave neutron stars with kicks in the ∼300−1000 km s −1 range. In addition, we find that these two classes are correlated with the gravitational mass of the neutron star. This correlation suggests that the survival of binary neutron star systems may in part be due to their lower kick speeds. We also find a correlation between the kick and both the mass dipole of the ejecta and the explosion energy. Furthermore, one channel of black hole birth leaves masses of ∼10 M ⊙ , is not accompanied by a neutrino-driven explosion, and experiences small kicks. A second channel is through a vigorous explosion that leaves behind a black hole with a mass of ∼3.0 M ⊙ kicked to high speeds. We find that the induced spins of nascent neutron stars range from seconds to ∼10 ms, but do not yet see a significant spin/kick correlation for pulsars. We suggest that if an initial spin biases the explosion direction, a spin/kick correlation would be a common byproduct of the neutrino mechanism of core-collapse supernovae. Finally, the induced spin in explosive black hole formation is likely large and in the collapsar range. This new 3D model suite provides a greatly expanded perspective and appears to explain some observed pulsar properties by default.

Astronomy & Astrophysics↗
DOE OSTI2024

Gravitational waves from binary neutron star mergers with a spectral equation of state

In numerical simulations of binary neutron star systems, the equation of state of the dense neutron star matter is an important factor in determining both the physical realism and the numerical accuracy of the simulations. Some equations of state used in simulations are C 2 or smoother in the pressure/density relationship function, such as a polytropic equation of state, but may not have the flexibility to model stars or remnants of different masses while keeping their radii within known astrophysical constraints. Other equations of state, such as tabular or piece-wise polytropic, may be flexible enough to model additional physics and multiple stars' masses and radii within known constraints, but are not as smooth, resulting in additional numerical error. We will study in this paper a recently developed family of equation of state, using a spectral expansion with sufficient free parameters to allow for a larger flexibility than current polytropic equations of state, and with sufficient smoothness to reduce numerical errors compared to tabulated or piece-wise polytropic equations of state. We perform simulations at three mass ratios with a common chirp mass, using two distinct spectral equations of state, and at multiple numerical resolutions. We evaluate the gravitational waves produced from these simulations, comparing the phase error between resolutions and equations of state, as well as with respect to analytical models. From our simulations we estimate that the phase difference at merger for binaries with a dimensionless weighted tidal deformability difference greater than Δ$\tilde{Λ}$=55 can be captured by the SpEC code for these equations of state.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Probing Dark Sectors with Neutron Stars

Tensions in the measurements of neutron and kaon weak decays, such as of the neutron lifetime, may speak to the existence of new particles and dynamics not present in the Standard Model (SM). In scenarios with dark sectors, particles that couple feebly to those of the SM appear. We offer a focused overview of such possibilities and describe how the observations of neutron stars, which probe either their structure or dynamics, limit them. In realizing these constraints, we highlight how the assessment of particle processes within dense baryonic matter impacts the emerging picture—and we emphasize both the flavor structure of the constraints and their broader connections to cogenesis models of dark matter and baryogenesis.

baryon number↗
DOE OSTI2020

Curvature-slope correlation of nuclear symmetry energy and its imprints on the crust-core transition, radius, and tidal deformability of canonical neutron stars

Background: The nuclear symmetry energy E sym (ρ) encodes information about the energy necessary to make nuclear systems more neutron-rich. While its slope parameter L at the saturation density ρ 0 of nuclear matter has been relatively well constrained by recent astrophysical observations and terrestrial nuclear experiments, its curvature K sym characterizing the E sym (ρ) around 2ρ 0 remains largely unconstrained. Over 520 calculations for E sym (ρ) using various nuclear theories and interactions in the literature have predicted several significantly different K sym –L correlations. Purpose: If a unique K sym –L correlation of E sym (ρ) can be firmly established, it will enable us to progressively better constrain the high-density behavior of E sym (ρ) using the available constraints on its slope parameter L. Here, we investigate if and by how much the different K sym –L correlations may affect neutron star observables. Method: A meta-model of nuclear Equation of States (EOSs) with three representative K sym –L correlation functions is used to generate multiple EOSs for neutron stars. We then examine effects of the K sym –L correlation on the crust-core transition density and pressure as well as the radius and tidal deformation of canonical neutron stars. Results: The K sym –L correlation affects significantly both the crust-core transition density and pressure. It also has strong imprints on the radius and tidal deformability of canonical neutron stars especially at small L values. The available data from LIGO/VIRGO and NICER set some useful limits for the slope L but can not distinguish the three representative K sym –L correlations considered.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
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 OSTI2026

Magnetic flux emergence in binary neutron star remnants

Using high-resolution AthenaK simulations of a twisted toroidal flux tube, we study the flux emergence of magnetic structures in the shear layer of a hot massive neutron star typical of a binary neutron star remnant. High-resolution simulations demonstrate that magnetic buoyant instabilities allow for emergence only for extremely large magnetic fields significantly exceeding 10 17 G, and more typical fields around 10 16 G are instead dominated by hydrodynamic effects. Because merger remnants tend to be stable against hydrodynamic convection, our work places strong limitations on the mechanisms by which massive binary neutron star remnants can produce the magnetically-driven outflows needed to power jets.

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 OSTI2023

Deducing the EOS of dense neutron star matter with machine learning

Abstract The interior of a neutron star is a unique astrophysical laboratory for studying matter at extreme densities and pressures beyond what is replicable in terrestrial experiments. While there is no direct way to simulate the interior of these stars, one promising avenue to learning more about the equation of state (EOS) of such matter is through X‐rays emitted from the star's surface. The current state‐of‐the‐art method for inference of EOS from a star's X‐ray spectra uses piece‐wise, simulation‐based likelihoods that rely on theoretical assumptions complicated by systematic uncertainties. To reduce the dimensionality of the problem, this method infers macroscopic properties of the star (mass and radius) from emitted X‐ray spectra, and from those quantities infers the EOS. This work approaches the same problem using machine learning techniques, demonstrating a series of enhancements to the current state‐of‐the‐art by realistic uncertainty quantification and reducing the need for theoretical assumptions. We also demonstrate novel inference of the EOS directly from high‐dimensional simulated X‐ray spectra from neutron stars that negate the need for a piece‐wise approach. This inference allows for a natural propagation of uncertainties from the X‐ray spectra by conditioning the discussed networks on realistic sources of uncertainty for each star.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Neural network reconstruction of the dense matter equation of state from neutron star observables

The Equation of State (EoS) of strongly interacting cold and hot ultra-dense QCD matter remains a major challenge in the field of nuclear astrophysics. With the advancements in measurements of neutron star masses, radii, and tidal deformabilities, from electromagnetic and gravitational wave observations, neutron stars play an important role in constraining the ultra-dense QCD matter EoS. Here, in this work, we present a novel method that exploits deep learning techniques to reconstruct the neutron star EoS from mass-radius (M-R) observations. We employ neural networks (NNs) to represent the EoS in a model-independent way, within the range of ~1-7 times the nuclear saturation density. The unsupervised Automatic Differentiation (AD) framework is implemented to optimize the EoS, so as to yield through TOV equations, an M-R curve that best fits the observations. We demonstrate that this method works by rebuilding the EoS on mock data, i.e., mass-radius pairs derived from a randomly generated polytropic EoS. The reconstructed EoS fits the mock data with reasonable accuracy, using just 11 mock M-R pairs observations, close to the current number of actual observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Reconstructing Masses of Merging Neutron Stars from Stellar r -process Abundance Signatures

We report that neutron star mergers (NSMs) are promising astrophysical sites for the rapid neutron-capture ("r") process, but can their integrated yields explain the majority of heavy-element material in the Galaxy? One method to address this question implements a forward approach that propagates NSM rates and yields along with stellar formation rates and compares those results with observed chemical abundances of r-process-rich, metal-poor stars. In this work, we take the inverse approach by utilizing r-process-element abundance ratios of metal-poor stars as input to reconstruct the properties -especially the masses - of their neutron star (NS) binary progenitors. This novel analysis provides an independent avenue for studying the population of the original NS binary systems that merged and produced the r-process material now incorporated in Galactic metal-poor halo stars. We use ratios of elements typically associated with the limited-r-process and the actinide region to those in the lanthanide region (i.e., Zr/Dy and Th/Dy) to probe the NS masses of the progenitor merger. We find that NSMs can account for all r-process material in metal-poor stars that display r-process signatures, while simultaneously reproducing the present-day distribution of double-NS systems. Notably, with our model assumptions and the studied stellar sample, we postulate that the most r-process enhanced stars (the r–II stars) on their own would require progenitor NSMs of asymmetric systems that are distinctly different from present ones in the Galaxy. We also explore variations to the model and find that the predicted degree of asymmetry is most sensitive to the electron fraction of the remnant disk wind.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Probing axionlike particles with multimessenger observations of neutron star mergers

Axion-like particles (ALPs) can be copiously produced in binary neutron star (BNS) mergers through nucleon-nucleon bremsstrahlung if the ALP-nucleon couplings 𝑔 𝑎⁢𝑁 are sizable. The ALP-photon coupling 𝑔 𝑎⁢𝛾 may trigger conversions of ultralight ALPs into photons in the magnetic fields of the merger remnant and of the Milky Way. This effect would lead to a potentially observable short gamma-ray signal, in coincidence with the gravitational-wave signal produced during the merging process. This event could be detected through multimessenger observation of BNS mergers employing the synergy between gravitational-wave detectors and gamma-ray telescopes. Here, we study the sensitivity of current and proposed MeV gamma-ray experiments to detect such a signal. As an explicit example, we consider ALP couplings related as in the Kim-Shifman-Vainshtein-Zakharov axion model, and show that in this case the proposed instruments can reach a sensitivity down to 𝑔 𝑎⁢𝛾 ≳ few ×10 −13 GeV −1 for 𝑚 𝑎 ≲ 10 −9 eV, comparable with the SN 1987A limit.

Axion-like particles↗