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

A Parameterized Neutrino Emission Model to Study Mass Ejection in Failed Core-collapse Supernovae

Some massive stars end their lives as failed core-collapse supernovae (CCSNe) and become black holes (BHs). Although in this class of phenomena the stalled supernova (SN) shock is not revived, the outer stellar envelope can still be partially ejected. This occurs because the hydrodynamic equilibrium of the star is disrupted by the gravitational mass loss of the protoneutron star (PNS) due to neutrino emission. We develop a simple parameterized model that emulates PNS evolution and its neutrino emission and use it to simulate failed CCSNe in spherical symmetry for a wide range of progenitor stars. Our model allows us to study mass ejection of failed CCSNe where the PNS collapses into a BH within ~100 ms and up to ~10 6 s. We perform failed CCSNe simulations for 262 different pre-SN progenitors and determine how the energy and mass of the ejecta depend on progenitor properties and the equation of state (EOS) of dense matter. In the case of a future failed CCSN observation, the trends obtained in our simulations can be used to place constraints on the pre-SN progenitor characteristics, the EOS, and on PNS properties at BH formation time.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Hydrodynamic simulations of electron-capture supernovae: progenitor and dimension dependence

ABSTRACT We present neutrino-transport hydrodynamic simulations of electron-capture supernovae (ECSNe) in flash with new two-dimensional (2D) collapsing progenitor models. These progenitor models feature the 2D modelling of oxygen-flame propagation until the onset of core collapse. We perform axisymmetric simulations with six progenitor models that, at the time of collapse, span a range of propagating flame front radii. For comparison, we also perform a simulation with the same set-up using the canonical, spherically symmetrical progenitor model n8.8. We found that the variations in the progenitor models inherited from simulations of stellar evolution and flame propagation do not significantly alter the global properties of the neutrino-driven ECSN explosion, such as the explosion energy (∼1.36–1.48 × 1050 erg) and the mass (∼0.017–0.018 M⊙) and composition of the ejecta. Due to aspherical perturbations induced by the 2D flame, the ejecta contains a small amount (≲1.8 × 10−3 M⊙) of low-Ye (0.35 < Ye < 0.4) component. The baryonic mass of the protoneutron star is ∼1.34 M⊙ (∼1.357 M⊙) with the new (n8.8) progenitor models when simulations end at ∼400 ms and the discrepancy is due to updated weak-interaction rates in the progenitor evolutionary simulations. Our results reflect the nature of ECSN progenitors containing a strongly degenerate oxygen–neon–magnesium (ONeMg) core and suggest a standardized ECSN explosion initialized by ONeMg core collapse. Moreover, we carry out a rudimentary three-dimensional simulation and find that the explosion properties are fairly compatible with the 2D counterpart. Our paper facilitates a more thorough understanding of ECSN explosions following the ONeMg core collapse, though more three-dimensional simulations are still needed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

New Covariant Density Functionals of Nuclear Matter for Compact Star Simulations

We generate three families of extended covariant density functionals of nuclear matter that have varying slope of symmetry energy and skewness at nuclear saturation density, but otherwise share the same basic parameters (symmetry energy, compressibility, saturation parameters, etc.) with the standard DDME2, DD2, and MPE functionals. Tables of the parameters of these new density functionals are given, which can be straightforwardly used in DDME2, DD2, and MPE parameterization-based codes. Furthermore, we provide tables of a large number of equations of state (81 for each family) that can be used in astrophysical simulations to assess the impact of variations of the not-well-known slope of symmetry energy and skewness of nuclear systems on the astrophysics of compact objects. We also provide tables of computed integral parameters (mass, radius, and tidal deformability) that can be used, e.g., for modeling gravitational waveforms. Finally, for the extended DDME2-based parameterization, we implement a first-order phase transition to quark matter to obtain a family of equations of state that accommodates a phase transition to quark matter. Analogous tables of the equations of state and integral parameters are provided for this case as well.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

The Breakthrough Listen Search For Intelligent Life Near the Galactic Center. I

A line of sight toward the Galactic Center (GC) offers the largest number of potentially habitable systems of any direction in the sky. The Breakthrough Listen program is undertaking the most sensitive and deepest targeted SETI surveys toward the GC. Here, we outline our observing strategies with Robert C. Byrd Green Bank Telescope (GBT) and Parkes telescope to conduct 600 hr of deep observations across 0.7–93 GHz. We report preliminary results from our survey for extraterrestrial intelligence (ETI) beacons across 1–8 GHz with 7.0 and 11.2 hr of observations with Parkes and GBT, respectively. With our narrowband drifting signal search, we were able to place meaningful constraints on ETI transmitters across 1–4 GHz and 3.9–8 GHz with EIRP limits of ≥4 × 10{sup 18} W among 60 million stars and ≥5 × 10{sup 17} W among half a million stars, respectively. For the first time, we were able to constrain the existence of artificially dispersed transient signals across 3.9–8 GHz with EIRP ≥1 × 10{sup 14} W/Hz with a repetition period ≤4.3 hr. We also searched our 11.2 hr of deep observations of the GC and its surrounding region for Fast Radio Burst–like magnetars with the DM up to 5000 pc cm{sup -3} with maximum pulse widths up to 90 ms at 6 GHz. We detected several hundred transient bursts from SGR J1745–2900, but did not detect any new transient bursts with the peak luminosity limit across our observed band of ≥10{sup 31} erg s{sup -1} and burst rate of ≥0.23 burst hr{sup -1}. These limits are comparable to bright transient emission seen from other Galactic radio-loud magnetars, constraining their presence at the GC.

47 OTHER INSTRUMENTATION↗
DOE OSTI2026

Uncertainties in the production of iron-group nuclides in core-collapse supernovae from Monte Carlo variations of reaction rates

Core-collapse supernovae, occurring at the end of massive star evolution, produce heavy elements, including those in the iron peak. Although the explosion mechanism is not yet fully understood, theoretical models can reproduce optical observations and observed elemental abundances. However, many nuclear reaction rates involved in explosive nucleosynthesis have large uncertainties, impacting the reliability of abundance predictions. To address this, we have previously developed a Monte Carlo-based nucleosynthesis code that accounts for reaction rate uncertainties and has been applied to nucleosynthesis processes beyond iron. Our framework is also well suited for studying explosive nucleosynthesis in supernovae. In this paper, we investigate 1D explosion models using the ‘PUSH method’ , focusing on progenitors with varying metallicities and initial masses around $M_{\rm ZAMS} = 16\, {\rm M}_{\odot }$. Detailed post-process nucleosynthesis calculations and Monte Carlo analyses are used to explore the effects of reaction rate uncertainties and to identify key reaction rates in explosive nucleosynthesis. We find that many reactions have little impact on the production of iron-group nuclei, as these elements are primarily synthesized in the nuclear statistical equilibrium. However, we identify a few ‘key reactions’ that significantly influence the production of radioactive nuclei, which may affect astrophysical observables. In particular, for the production of ${}^{44}{\rm Ti}$, we confirm that several traditionally studied nuclear reactions have a strong impact. However, determining a single reaction rate is insufficient to draw a definitive conclusion.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Exploring the role of 𝑑* hexaquarks on quark deconfinement and hybrid stars

We investigate the impact of the d*( 2380) hexaquark on the equation of state (EOS) of dense matter within hybrid stars (HSs) using the chiral mean-field model (CMF). The hexaquark is included as a new degree of freedom in the hadronic phase, and its influence on the deconfinement transition to quark matter is explored. We reparametrize the CMF model to ensure compatibility with recent astrophysical constraints, including the observation of massive pulsars and gravitational wave events. Our results show that the presence of d* significantly modifies the EOS, leading to a softening at high densities and a consequent reduction in the predicted maximum stellar masses. Furthermore, we examine the possibility of a first-order deconfinement phase transition within the context of the extended stability branch of slow stable HSs (SSHSs). We find that the presence of hexaquarks can delay the deconfinement phase transition and reduce the associated energy density gap, affecting the structure and stability of HSs. Our results suggest that, as the hexaquark appearance tends to destabilize stellar configurations, fine-tuning of model parameters is required to obtain both the presence of hexaquarks and quark deconfinement in these systems. In this scenario, the SSHS branch plays a crucial role in obtaining HSs with hexaquarks that satisfy current astrophysical constraints. Our work provides new insights into the role of exotic particles like d* in dense matter and the complex interplay between hadronic and quark degrees of freedom inside compact stellar objects.

Nuclear astrophysics↗
DOE OSTI2022

Equation of state of neutron-rich matter in d -dimensions

Nuclear systems under constraints, with high degrees of symmetries and/or collectivities may be considered as moving effectively in spaces with reduced spatial dimensions. We first derive analytical expressions for the nucleon specific energy E 0 (ρ), pressure P 0 (ρ), incompressibility coefficient K 0 (ρ) and skewness coefficient J 0 (ρ) of symmetric nucleonic matter (SNM), the quadratic symmetry energy E sym (ρ), its slope parameter L(ρ) and curvature coefficient K sym (ρ) as well as the fourth-order symmetry energy E sym,4 (ρ) of neutron-rich matter in general d spatial dimensions (abbreviated as "dD") in terms of the isoscalar and isovector parts of the isospin-dependent single-nucleon potential according to the generalized Hugenholtz-Van Hove (HVH) theorem. The equation of state (EOS) of nuclear matter in dD can be linked to that in the conventional 3-dimensional (3D) space by the ε-expansion which is a perturbative approach successfully used previously in treating second-order phase transitions and related critical phenomena and more recently in studying the EOS of cold atoms. The ε-expansion of nuclear EOS in dD based on a reference dimension d f = d – ε is shown to be effective with –1 ≲ ε ≲ 1 starting from 1 ≲ d f ≲ 3 in comparison with the exact expressions derived using the HVH theorem. Moreover, the EOS of SNM (with/without considering its potential part) is found to be reduced (enhanced) in lower (higher) dimensions, indicating in particular that the many-nucleon system tends to be deeper bounded but saturate at higher densities in spaces with lower dimensions. Here, the links between the EOSs in 3D and dD spaces from the ε-expansion provide new perspectives to the EOS of neutron-rich matter.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Atomic masses with machine learning for the astrophysical r process

The astrophysical r process plays a vital role in the production of heavy elements. Modeling of the r process is sensitive to masses and further requires knowledge of masses beyond current experimental reach. Therefore, simulations of the r process offer a unique test bed for predicting mass extrapolations. We take a Machine-Learning (ML) approach to model the masses across the entire chart of nuclides. For the first time, we simulate r-process nucleosynthesis with a physics-based ML mass model. We compare simulated abundances to solar data in order to evaluate the model's performance far from stability. The resulting r-process abundances up to thorium and uranium qualitatively match those of the observed solar system abundance pattern, with the characteristic peaks well positioned. We propagate the mass uncertainties obtained from the ML model to r-process abundance yields to estimate an uncertainty band associated with our approach. The size of the uncertainty band is approximately one order of magnitude which aligns with the uncertainty reported using alternative techniques.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Implications of a weakening N = 126 shell closure away from stability for r -process astrophysical conditions

The formation of the third r-process abundance peak near A ∼ 195 is highly sensitive to both nuclear structure far from stability and the astrophysical conditions that produce the heaviest elements. In particular, the N = 126 shell closure plays a crucial role in shaping this peak. Experimental data hints that the shell weakens as proton number departs from Z = 82, a trend largely missed by global mass models. To investigate its impact on r-process nucleosynthesis, we employ both standard global models with strong closures and modified Duflo-Zuker (DZ) models that reproduce the weakening, combined with three sets of β − -decay rates. Strong shell closures generate sharply peaked abundances, whereas weakened closures consistent with the experimental trend produce broader, flatter patterns. Accurately reproducing the solar third peak under weakened shell strength requires sufficiently neutron-rich conditions that significant fission occurs, and slower decay rates. These results demonstrate that a weakening N = 126 shell closure away from stability imposes significant constraints on the astrophysical environments of the r-process and underscores the need for precise mass measurements and improved characterization of β − -decay properties in this region.

N = 126 closed shells↗
DOE OSTI2024

Pulsar timing anomalies: a window into baryon number violation

Here, we investigate the influence of a specific class of slow Baryon Number Violation (BNV)—one that induces quasi-equilibrium evolution — on pulsar spin characteristics. This work reveals how BNV can potentially alter observable parameters, including spin-down rates, the second derivative of spin frequency, and braking indices of pulsars. Moreover, we demonstrate that BNV could lead to anomalies in pulsar timing, along with a wide array of braking indices, both positive and negative. In addition, we examine the possibility of pulsar spin-up due to BNV, which may result in a novel mechanism for the revival of “dead” pulsars. We conclude by assessing the sensitivity required for future pulsar timing efforts to detect such BNV effects, thus highlighting the potential for pulsars to serve as laboratories for testing fundamental physics.

cosmology↗
DOE OSTI2025

On the motion of compact objects in relativistic viscous fluids

We present a world-line effective field theory of compact objects moving relativistically through a viscous fluid. The theory is valid when velocity gradients are small compared to the inverse size of the object. Working within the EFT eliminates the need to solve a boundary value problem by turning all interactions between the fluid and the object into a source term in the action. We use the EFT to derive the relativistic equations of motion for a compact object immersed in a viscous fluid in a curved background, when the relative velocity of the object and the fluid is small compared to the speed of light.

astrophysical black holes↗
DOE OSTI2020

A search for optical and near-infrared counterparts of the compact binary merger GW190814

Here, we report on our observing campaign of the compact binary merger GW190814, detected by the Advanced LIGO and Advanced Virgo detectors on 2019 August 14. This signal has the best localization of any observed gravitational wave (GW) source, with a 90 per cent probability area of 18.5 deg 2 , and an estimated distance of ≈240 Mpc. We obtained wide-field observations with the Deca-Degree Optical Transient Imager (DDOTI) covering 88 per cent of the probability area down to a limiting magnitude of w = 19.9 AB. Nearby galaxies within the high probability region were targeted with the Lowell Discovery Telescope (LDT), whereas promising candidate counterparts were characterized through multicolour photometry with the Reionization and Transients InfraRed (RATIR) and spectroscopy with the Gran Telescopio de Canarias (GTC). We use our optical and near-infrared limits in conjunction with the upper limits obtained by the community to constrain the possible electromagnetic counterparts associated with the merger. A gamma-ray burst seen along its jet’s axis is disfavoured by the multiwavelength data set, whereas the presence of a burst seen at larger viewing angles is not well constrained. Although our observations are not sensitive to a kilonova similar to AT2017gfo, we can rule out high-mass (>0.1 M ⊙ ) fast-moving (mean velocity ≥0.3c) wind ejecta for a possible kilonova associated with this merger.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Actinide opacities for modelling the spectra and light curves of kilonovae

ABSTRACT We extend previous ab initio calculations of lanthanide opacities to include a complete set of actinide opacities for use in the modelling of kilonova (KN) light curves and spectra. Detailed, fine-structure line features are generated using the configuration-interaction approach. These actinide opacities display similar trends to those observed for lanthanide opacities, such as the lighter actinides producing higher opacity than the heavier ones for relevant conditions in the dynamical ejecta. A line-binned treatment is employed to pre-compute opacity tables for 14 actinide elements (89 ≤ Z ≤ 102) over a grid of relevant temperatures and densities. These tabular opacities will be made publicly available for general usage in KN modelling. We demonstrate the usefulness of these opacities in KN simulations by exploring the sensitivity of light curves and spectra to different actinide abundance distributions that are predicted by different nuclear theories, as well as to different choices of ejecta mass and velocity. We find very little sensitivity to the two considered distributions, indicating that opacities for actinides with Z ≥ 99 do not contribute strongly. On the other hand, a single actinide element, protactinium, is found to produce faint spectral features in the far-infrared at late times (5–7 d post merger). More generally, we find that the choice of ejecta mass and velocity have the most significant effect on KN emission for this study.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

The invisible black widow PSR J1720−0534: implications for the electron density towards the North Polar Spur

Radio emission from pulsars can be used to map out their distances through dispersion measure (DM), which quantifies the amount of radio pulse dispersion. However, this method relies on accurately modelling the free electron density in the line of sight. Here, we present a detailed study of the multiwavelength emission from PSR J1720−0534, a black widow compact binary millisecond pulsar discovered in 2021, which the latest electron density model of the Galaxy places at only 191 pc. We obtained and analysed deep multiwavelength observations in the γ-ray ( Fermi -Large Area Telescope, 2008–2022), optical (Las Cumbres Observatory, 2.7 h), near-infrared (Nordic Optical Telescope, 3.5 h), and X-ray ( Swift -X-Ray Telescope, 10 ks) bands. We found no significant detection of γ-ray, optical, near-infrared, or X-ray counterparts around the radio-timing position of PSR J1720−0534, which we thus nickname ‘the invisible black widow’. Employing the most constraining near-infrared limit ( J > 23.4 mag), we established a lower limit on the source distance, d > 1.1 kpc, assuming conservative properties for the black widow companion star. This distance lower limit differs drastically (by a factor of more than 5) from the Yao et al. DM distance estimate. We attribute this difference to the inclusion in the Yao et al. model of a large and dense component towards the North Polar Spur. Considering our results and recent parallax distances to other pulsars in this direction, we argue that such a local and large component in the electron density model of the Galaxy is unnecessary.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Plasma instabilities dominate radioactive transients magnetic fields: the self-confinement of leptons in Type Ia and core-collapse supernovae, and kilonovae

The light curves of radioactive transients, such as supernovae and kilonovae, are powered by the decay of radioisotopes, which release high-energy leptons through $\beta ^+$ and $\beta ^-$ decays. These leptons deposit energy into the expanding ejecta. As the ejecta density decreases during expansion, the plasma becomes collisionless, with particle motion governed by electromagnetic forces. In such environments, strong or turbulent magnetic fields are thought to confine particles, though the origin of these fields and the confinement mechanism have remained unclear. Using fully kinetic particle-in-cell (PIC) simulations, we demonstrate that plasma instabilities can naturally confine high-energy leptons. These leptons generate magnetic fields through plasma streaming instabilities, even in the absence of pre-existing fields. The self-generated magnetic fields slow lepton diffusion, enabling confinement, and transferring energy to thermal electrons and ions. Our results naturally explain the positron trapping inferred from late-time observations of thermonuclear and core-collapse supernovae. Furthermore, they suggest potential implications for electron dynamics in the ejecta of kilonovae. We also estimate synchrotron radio luminosities from positrons for Type Ia supernovae and find that such emission could only be detectable with next-generation radio observatories from a Galactic or local-group supernova in an environment without any circumstellar material.

instabilities↗
DOE OSTI2025

Outflowing photoionized plasma in Circinus X-1 using the high-resolution X-ray spectrometer Resolve onboard XRISM and the radiative transfer code cloudy

High-resolution X-ray spectroscopy is a key to understanding the mass inflow and outflow of compact objects. Spectral lines carry information about the ionization, density, and velocity structures through their intensity ratios and profiles. They are formed in non-local thermodynamic equilibrium conditions under the intense radiation field from the compact objects; thus radiative transfer (RT) calculation is a requisite for proper interpretations. We present such a study for a low-mass X-ray binary, Circinus X-1, from which the P Cygni profile was discovered using the X-ray grating spectrometer onboard Chandra. We observed the source using the X-ray microcalorimeter onboard XRISM at an orbital phase of 0.93–0.97 and revealed many spectral features unidentified before: the higher series transitions ($n{\rightarrow }1$; $n>2$) of highly ionized (H- and He-like) S, Ca, Ar, and Fe in emission and absorption, the Fe K$\alpha$ and K$\beta$ inner-shell excitation absorption of mildly ionized (O- to Li-like) Fe, and resolved fine-structure level transitions in the Fe Ly$\alpha$ and He$\alpha$ complexes. They blend with each other at different velocity shifts on top of apparently variable continuum emission that changed its flux by an order of magnitude within a 70 ks telescope time. Despite such complexity in the observed spectra, most of them can be explained by a simple model consisting of the photoionized plasma outflowing at $\sim$300 km s$^{-1}$ and the variable blocking material in the line of sight of the incident continuum emission from the accretion disk. We demonstrate this with the aid of the RT code cloudy for the line ratio diagnostics and spectral fitting. We further constrain the physical parameters of the outflow and argue that the outflow is launched close to the outer edge of the accretion disk and can be driven radiatively by being assisted by the line force calculated using the RT simulation.

X-rays: binaries↗
DOE OSTI2026

Energy shift of Fe-K fluorescence lines due to low ionization demonstrated with XRISM in Centaurus X-3

The Fe K$\alpha$ fluorescence line at 6.4 keV is a powerful probe of cold matter surrounding X-ray sources and has been widely used in various astrophysical contexts. The X-ray microcalorimeter spectrometer onboard XRISM can measure line shifts with unprecedented precision of ${\sim }$0.2 eV, equivalent to a line-of-sight velocity of ${\sim }$10 km s$^{-1}$. At this level of accuracy, however, several factors that influence the line energy must be carefully considered prior to astrophysical interpretation. One such important factor is the ionization degree, Fe$^{q+}$. The K$\alpha$ line shifts redward by ${\sim }$4 eV as q increases from 0 (neutral) to 8 (Ar-like). Additionally, the accompanying Fe K$\beta$ line at 7.06 keV shifts blueward by ${\sim }$30 eV from $q=0$ to 8. We demonstrate that this effect is actually observable in the XRISM data of the high-mass X-ray binary Centaurus X-3 (Cen X-3). We advocate that the differential energy shift between the K$\alpha$ and K$\beta$ line provides a robust estimate of q by decoupling from other effects that shift the two lines in the same direction. We derived $q \sim 5$ (Sc-like) for the fluorescing matter by comparing the observation with atomic structure calculations of our own and in the literature. By accounting for the derived charge state and the corresponding shift in the rest-frame line energy, we made corrections for this effect and reached a consistent residual shift among the K$\alpha$, K$\beta$, and the optical measurement attributable to the systemic velocity of the system. Consequently, we obtained a new constraint on the location of the cold matter. This ionization effect needs to be assessed in all use cases of the Fe K$\alpha$ line shift beyond Cen X-3, and the proposed metric is generally applicable to all of them.

Astronomy and AstroPhysics↗