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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

All-sky search for continuous gravitational waves from isolated neutron stars using Advanced LIGO and Advanced Virgo O3 data

We present results of an all-sky search for continuous gravitational waves which can be produced by spinning neutron stars with an asymmetry around their rotation axis, using data from the third observing run of the Advanced LIGO and Advanced Virgo detectors. Four different analysis methods are used to search in a gravitational-wave frequency band from 10 to 2048 Hz and a first frequency derivative from -10 -8 to 10 -9 Hz/s. No statistically significant periodic gravitational-wave signal is observed by any of the four searches. As a result, upper limits on the gravitational-wave strain amplitude ℎ 0 are calculated. The best upper limits are obtained in the frequency range of 100 to 200 Hz and they are ~1.1 ×10 -25 at 95% confidence level. The minimum upper limit of 1.10 ×10 -25 is achieved at a frequency 111.5 Hz. In conclusion, we also place constraints on the rates and abundances of nearby planetary- and asteroid-mass primordial black holes that could give rise to continuous gravitational-wave signals.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Central speed of sound, the trace anomaly, and observables of neutron stars from a perturbative analysis of scaled Tolman-Oppenheimer-Volkoff equations

The central speed of sound (SS) measures the stiffness of the equation of state (EOS) of superdense neutron star (NS) matter. Its variations with density and radial coordinate in NSs in conventional analyses often suffer from uncertainties of the specific nuclear EOS used. Using the central SS and NS mass/radius scaling obtained from solving perturbatively the scaled Tolman-Oppenheimer-Volkoff (TOV) equations, we study the variations of SS, trace anomaly and several closely related properties of NSs in an EOSmodel- independent manner. We find that the SS increases with the reduced central pressure Pˆ c ≡ P c =ε c (scaled by the central energy density εc), and the conformal bound for SS tends to break down for NSs with masses higher than about 1.9M ⊙ . The ratio P=ε is upper bounded as P=ε ≲ 0.374 around the centers of stable NSs. We demonstrate that it is an intrinsic property of strong field gravity and is more relevant than the perturbative QCD bound on it. While a sharp phase transition at high densities characterized by a sudden vanishing of SS in cores of massive NSs are basically excluded, the probability for a continuous crossover signaled by a peaked radial profile of SS is found to be enhanced as Pˆ c decreases, implying it likely happens near the centers of massive NSs. Moreover, a new and more stringent causality boundary as R max =km ≳ 4.73$M$$^{max}_{NS}$ =M ⊙ þ 1.14 for the NS mass-radius curve is found to be excellently consistent with observational data on NS masses and radii. Here, new constraints on the ultimate energy density and pressure allowed in NSs before collapsing into black holes are obtained and compared with earlier predictions in the literature

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

AT2017gfo: Bayesian inference and model selection of multicomponent kilonovae and constraints on the neutron star equation of state

The joint detection of the gravitational wave GW170817, of the short γ-ray burst GRB170817A and of the kilonova AT2017gfo, generated by the the binary neutron star (NS) merger observed on 2017 August 17, is a milestone in multimessenger astronomy and provides new constraints on the NS equation of state. We perform Bayesian inference and model selection on AT2017gfo using semi-analytical, multicomponents models that also account for non-spherical ejecta. Observational data favour anisotropic geometries to spherically symmetric profiles, with a log-Bayes’ factor of ~104, and favour multicomponent models against single-component ones. The best-fitting model is an anisotropic three-component composed of dynamical ejecta plus neutrino and viscous winds. Using the dynamical ejecta parameters inferred from the best-fitting model and numerical–relativity relations connecting the ejecta properties to the binary properties, we constrain the binary mass ratio to q < 1.54 and the reduced tidal parameter to $120\lt \tilde{\Lambda }\lt 1110$. Finally, we combine the predictions from AT2017gfo with those from GW170817, constraining the radius of a NS of 1.4 M ⊙ to 12.2 ± 0.5 km (1σ level). This prediction could be further strengthened by improving kilonova models with numerical-relativity information.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Neutron star mass estimates from gamma-ray eclipses in spider millisecond pulsar binaries

Reliable neutron star mass measurements are key to determining the equation of state of cold nuclear matter, but such measurements are rare. Black widows and redbacks are compact binaries consisting of millisecond pulsars and semi-degenerate companion stars. Spectroscopy of the optically bright companions can determine their radial velocities, providing inclination-dependent pulsar mass estimates. Although inclinations can be inferred from subtle features in optical light curves, such estimates may be systematically biased due to incomplete heating models and poorly understood variability. Using data from the Fermi Large Area Telescope, we have searched for gamma-ray eclipses from 49 spider systems, discovering significant eclipses in 7 systems, including the prototypical black widow PSR B1957+20. Gamma-ray eclipses require direct occultation of the pulsar by the companion, and so the detection, or significant exclusion, of a gamma-ray eclipse strictly limits the binary inclination angle, providing new robust, model-independent pulsar mass constraints. For PSR B1957+20, the eclipse implies a much lighter pulsar (1.81 ± 0.07 solar masses) than inferred from optical light curve modelling.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Thermalization and annihilation of dark matter in neutron stars

The capture of dark matter, and its subsequent annihilation, can heat old, isolated neutron stars. In order for kinetic heating to be achieved, the captured dark matter must undergo sufficient scattering to deposit its kinetic energy in the star. We find that this energy deposit typically occurs quickly, for most of the relevant parameter space. In order for appreciable annihilation heating to also be achieved, the dark matter must reach a state of capture-annihilation equilibrium in the star. We show that this can be fulfilled for all types of dark matter - baryon interactions. This includes cases where the scattering or annihilation cross sections are momentum or velocity suppressed in the non-relativistic limit. Importantly, we find that capture-annihilation equilibrium, and hence maximal annihilation heating, can be achieved without complete thermalization of the captured dark matter. For scattering cross sections that saturate the capture rate, we find that capture-annihilation equilibrium is typically reached on a timescale of less than 1 year for vector interactions and 10 4 years for scalar interactions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Mirror neutron stars

The fundamental nature of dark matter is entirely unknown. A compelling candidate is twin Higgs mirror matter, invisible hidden-sector cousins of the Standard Model particles and forces. This predicts mirror neutron stars made entirely of mirror nuclear matter. Here we find their structure using realistic equations of state, robustly modified based on first-principle quantum chromodynamic calculations, for the first time. This allows us to predict their gravitational wave signals, demonstrating an impressive discovery potential and ability to probe dark sectors connected to the hierarchy problem.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Empirical radius formulas for canonical neutron stars from bidirectionally selecting features of equations of state in extended Bayesian analyses of observational data

Significant advancement in Bayesian inference of nuclear equation of state (EOS) from gravitational wave and x-ray observations of neutron stars (NSs) has been made by the nuclear astrophysics community especially since GW170817. By extending the traditional Bayesian analysis which normally ends at presenting the marginalized posterior probability distribution functions (PDFs) of individual EOS parameters and their correlations (or sometimes only the Pearson correlation coefficients which are only reliably useful when the variables are linearly correlated while they are actually often not), we search for a data-driven and robust empirical formula for the radius 𝑅 1.4 of canonical NSs in terms of the characteristic EOS parameters (features). We also identify the single most important but currently poorly known EOS parameter for determining the 𝑅 1.4 . Using three regression-model-building methodologies: bidirectional stepwise feature selection, least absolute shrinkage selection operator (LASSO) regression, and neural network regression on a large set of posterior EOSs and the corresponding 𝑅 1.4 values inferred from earlier comprehensive Bayesian analyses of NS observational data, we systematically and rigorously develop the most probable 𝑅 1.4 formulas with varying statistical accuracy and technical complexity. Here, the most important EOS parameters for determining 𝑅 1.4 are found consistently in each of the feature selection processes to be (in order of decreasing importance): curvature 𝐾 sym , slope 𝐿, skewness 𝐽 sym of nuclear symmetry energy, skewness 𝐽 0 , incompressibility 𝐾 0 of symmetric nuclear matter, and the magnitude 𝐸 sym ⁡(𝜌 0 ) of symmetry energy at the saturation density 𝜌 0 of nuclear matter.

Bayesian methods↗
DOE OSTI2023

Deducing neutron star equation of state parameters directly from telescope spectra with uncertainty-aware machine learning

Neutron stars provide a unique laboratory for studying matter at extreme pressures and densities. While there is no direct way to explore their interior structure, X-rays emitted from these stars can indirectly provide clues to the equation of state (EOS) of the superdense nuclear matter through the inference of the star's mass and radius. However, inference of EOS directly from a star's X-ray spectra is extremely challenging and is complicated by systematic uncertainties. The current state of the art is to use simulation-based likelihoods in a piece-wise method which relies on certain theoretical assumptions and simplifications about the uncertainties. It first infers the star's mass and radius to reduce the dimensionality of the problem, and from those quantities infer the EOS. Here we demonstrate a series of enhancements to the state of the art, in terms of realistic uncertainty quantification and a path towards circumventing the need for theoretical assumptions to infer physical properties with machine learning. We also demonstrate novel inference of the EOS directly from the high-dimensional spectra of observed stars, avoiding the intermediate mass-radius step. Our network is conditioned on the sources of uncertainty of each star, allowing for natural and complete propagation of uncertainties to the EOS.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

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

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

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

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

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

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Thermodynamics conditions of matter in the neutrino decoupling region during neutron star mergers

In this work we investigate the thermodynamics conditions at which neutrinos decouple from matter in neutron star merger remnants by post-processing results of merger simulations. We find that the matter density and the neutrino energies are the most relevant quantities in determining the decoupling surface location. For mean energy neutrinos (~9, 15 and 24 MeV for ν e , $\bar{ν}$ e and ν μ,τ , respectively) the transition between diffusion and free-streaming conditions occurs around 10 11 g cm -3 for all neutrino species. Weak and thermal equilibrium freeze-out occurs deeper (several 10 12 g cm -3 ) for heavy-flavor neutrinos than for $\bar{ν}$ e and ν e (≳10 11 g cm -3 ). Decoupling temperatures are broadly in agreement with the average neutrino energies, with softer equations of state characterized by ~1 MeV larger decoupling temperatures. Neutrinos streaming at infinity with different energies come from different remnant parts. While low-energy neutrinos (~3 MeV) decouple at ρ~10 13 g cm -3 , T~10 MeV and Y e ≲0.1 close to weak equilibrium, high-energy ones (~50 MeV ) decouple from the disk at ρ~10 9 g cm -3 , T~2 MeV and Y e ≳0.25. The presence of a massive NS or a BH influences the neutrino thermalization. While in the former case decoupling surfaces are present for all relevant energies, the lower maximum density ( ≲10 12 g cm -3 ) in BH-torus systems does not allow softer neutrinos to thermalize and diffuse.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Evidence from Disrupted Halo Dwarfs that r-process Enrichment via Neutron Star Mergers is Delayed by ≳500 Myr

The astrophysical origins of r-process elements remain elusive. Neutron star mergers (NSMs) and special classes of core-collapse supernovae (rCCSNe) are leading candidates. Due to these channels' distinct characteristic timescales (rCCSNe: prompt, NSMs: delayed), measuring r-process enrichment in galaxies of similar mass but differing star formation durations might prove informative. Two recently discovered disrupted dwarfs in the Milky Way's stellar halo, Kraken and Gaia-Sausage Enceladus (GSE), afford precisely this opportunity: Both have M * ≈ 10 8 M ⊙ but differing star formation durations of ≈2 Gyr and ≈3.6 Gyr. Here we present R ≈ 50,000 Magellan/MIKE spectroscopy for 31 stars from these systems, detecting the r-process element Eu in all stars. Stars from both systems have similar [Mg/H] ≈ -1, but Kraken has a median [Eu/Mg] ≈ -0.1 while GSE has an elevated [Eu/Mg] ≈ 0.2. With simple models, we argue NSM enrichment must be delayed by 500–1000 Myr to produce this difference. rCCSNe must also contribute, especially at early epochs, otherwise stars formed during the delay period would be Eu free. In this picture, rCCSNe account for ≈50% of the Eu in Kraken, ≈25% in GSE, and ≈15% in dwarfs with extended star formation durations like Sagittarius. The inferred delay time for NSM enrichment is 10×–100× longer than merger delay times from stellar population synthesis—this is not necessarily surprising because the enrichment delay includes time taken for NSM ejecta to be incorporated into subsequent generations of stars. For example, this may be due to natal kicks that result in r-enriched material deposited far from star-forming gas, which then takes ≈10 8 –10 9 yr to cool in these galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Superfluidity in Disordered Neutron Stars Crusts

Nonequilibrium conditions imposed by neutrino cooling through the liquid-solid transition lead to disorder in the solid crust of neutron stars. Disorder reduces the superfluid fraction, $\rho_s/\rho$, at densities above that of neutron drip, $\rho_d \approx 4\times 10^{11}\,g/cm^3$. For an amorphous solid crust the suppression of $\rho_s$ is small, except in the highest density regions of the crust. In contrast to the strong reduction in neutron conduction predicted for coherent Bragg scattering in a crystalline crust, the disordered solid crust supports sufficient neutron superfluid density to account for pulsar glitches.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Unraveling trace anomaly of supradense matter via neutron star compactness scaling

The trace anomaly Δ ≡ 1/3 −𝑃/𝜖 =1/3 −𝜙 quantifies the possibly broken conformal symmetry in supradense matter under pressure 𝑃 at energy density 𝜖. Perturbative QCD (pQCD) predicts a vanishing Δ at extremely high energy or baryon densities when the conformal symmetry is realized but its behavior at intermediate densities reachable in neutron stars (NSs) is still very uncertain. The extraction of Δ from NS observations strongly depends on the employed model for nuclear equation of state (EOS). Using the IPAD-TOV method based on an intrinsic and perturbative analysis of the dimensionless (IPAD) Tolman-Oppenheimer-Volkoff (TOV) equations that are further verified numerically by using 10 5 EOSs generated randomly with a metamodel in a very broad EOS parameter space constrained by terrestrial nuclear experiments and astrophysical observations, here we first show that the compactness 𝜉 ≡ 𝐺⁡𝑀 NS /𝑅⁢𝑐 2 ≡ 𝑀 NS /𝑅 of a NS with mass 𝑀 NS and radius 𝑅 scales very accurately with $\bar{Π}$ c ≡ $Π$ c · (1 +18⁢X/25) ≡ X/(1 +3⁢X 2 +4⁢X) · (1 +18⁢X/25) where X ≡ 𝜙 c = 𝑃 c /𝜖 c is the ratio of pressure over energy density at NS centers. The scaling of NS compactness thus enables one to readily read off the central trace anomaly Δ c = 1/3 −X directly from the observational data of either the mass-radius or red-shift measurements. Finally, we then demonstrate indeed that the available NS data themselves from recent X-ray and gravitational wave observations can determine model insensitively the trace anomaly as a function of energy density in NS cores, providing a stringent test of existing NS models and a clear guidance in a new direction for further understanding the nature and EOS of supradense matter.

nuclear astrophysics↗
DOE OSTI2022

A Systematic Exploration of Kilonova Candidates from Neutron Star Mergers during the Third Gravitational-wave Observing Run

We present a comprehensive analysis of 653 optical candidate counterparts reported during the third gravitational-wave (GW) observing run. Our sample concentrates on candidates from the 15 events (published in GWTC-2, GWTC-3, or not retracted on GraceDB) that had a >1% chance of including a neutron star in order to assess their viability as true kilonovae. In particular, we leverage tools available in real time, including pre-merger detections and cross-matching with catalogs (i.e., point-source, variable-star, quasar and host-galaxy redshift data sets), to eliminate 65% of candidates in our sample. We further employ spectroscopic classifications, late-time detections, and light-curve behavior analyses and conclude that 66 candidates remain viable kilonovae. These candidates lack sufficient information to determine their classifications, and the majority would require luminosities greater than that of AT 2017gfo. Pre-merger detections in public photometric survey data and comparison of cataloged host-galaxy redshifts with the GW event distances are critical to incorporate into vetting procedures, as these tools eliminated >20% and >30% of candidates, respectively. We expect that such tools that leverage archival information will significantly reduce the strain on spectroscopic and photometric follow-up resources in future observing runs. Finally, we discuss the critical role prompt updates from GW astronomers to the EM community play in reducing the number of candidates requiring vetting.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

How Well Do We Know the Neutron-Matter Equation of State at the Densities Inside Neutron Stars? A Bayesian Approach with Correlated Uncertainties

Here, we introduce a new framework for quantifying correlated uncertainties of the infinite-matter equation of state derived from chiral effective field theory (𝜒⁢EFT ). Bayesian machine learning via Gaussian processes with physics-based hyperparameters allows us to efficiently quantify and propagate theoretical uncertainties of the equation of state, such as 𝜒⁢EFT truncation errors, to derived quantities. We apply this framework to state-of-the-art many-body perturbation theory calculations with nucleon-nucleon and three-nucleon interactions up to fourth order in the 𝜒⁢EFT expansion. This produces the first statistically robust uncertainty estimates for key quantities of neutron stars. We give results up to twice nuclear saturation density for the energy per particle, pressure, and speed of sound of neutron matter, as well as for the nuclear symmetry energy and its derivative. At nuclear saturation density, the predicted symmetry energy and its slope are consistent with experimental constraints.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

The landscape of disc outflows from black hole–neutron star mergers

ABSTRACT We investigate mass ejection from accretion discs formed in mergers of black holes (BHs) and neutron stars (NSs). The third observing run of the LIGO/Virgo interferometers provided BH–NS candidate events that yielded no electromagnetic (EM) counterparts. The broad range of disc configurations expected from BH–NS mergers motivates a thorough exploration of parameter space to improve EM signal predictions. Here we conduct 27 high-resolution, axisymmetric, long-term hydrodynamic simulations of the viscous evolution of BH accretion discs that include neutrino emission/absorption effects and post-processing with a nuclear reaction network. In the absence of magnetic fields, these simulations provide a lower limit to the fraction of the initial disc mass ejected. We find a nearly linear inverse dependence of this fraction on disc compactness (BH mass over initial disc radius). The dependence is related to the fraction of the disc mass accreted before the ouflow is launched, which depends on the disc position relative to the innermost stable circular orbit. We also characterize a trend of decreasing ejected fraction and decreasing lanthanide/actinide content with increasing disc mass at fixed BH mass. This trend results from a longer time to reach weak freezout and an increasingly dominant role of neutrino absorption at higher disc masses. We estimate the radioactive luminosity from the disc outflow alone available to power kilonovae over the range of configurations studied, finding a spread of two orders of magnitude. For most of the BH–NS parameter space, the disc outflow contribution is well below the kilonova mass upper limits for GW190814.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

A Unified Model of Kilonovae and Gamma-Ray Bursts in Binary Mergers Establishes Neutron Stars as the Central Engines of Short GRBs

We expand the theoretical framework by O. Gottlieb et al., which connects binary merger populations with long and short binary gamma-ray bursts (lbGRBs and sbGRBs, respectively), incorporating kilonovae (KNe) as a key diagnostic tool. We show that lbGRBs, powered by massive accretion disks around black holes (BHs), should be accompanied by bright, red KNe. In contrast, sbGRBs—if also powered by BHs—would produce fainter, red KNe, potentially biasing against their detection. However, magnetized hypermassive neutron star (HMNS) remnants that precede BH formation can produce jets with power (P NS ≈ 10 51 erg s –1 ) and Lorentz factor (Γ > 10) likely compatible with sbGRB observations, and would result in distinctly bluer KNe, offering a pathway to identifying the sbGRB central engine. Recent modeling by J. C. Rastinejad et al. found luminous red KNe consistently accompany lbGRBs, supporting their origin in BH-massive disk systems, likely following a short-lived HMNS phase. The preferential association of sbGRBs with comparably luminous KNe argues against the BH engine hypothesis for sbGRBs, while the bluer hue of these KNe provides additional support for an HMNS-driven mechanism. Within this framework, BH–NS mergers likely contribute exclusively to the lbGRB population with red KNe. Our findings suggest that GW170817 may, in fact, have been an lbGRB to on-axis observers. Finally, we discuss major challenges faced by alternative lbGRB progenitor models, such as white dwarf–NS or white dwarf–BH mergers and accretion-induced collapse forming magnetars, which fail to align with observed GRB timescales, energies, and KN properties.

79 ASTRONOMY AND ASTROPHYSICS↗