Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nuclear matter in 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.

2,811 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Effect of magnetic fields on Urca rates in neutron star mergers

Isospin-equilibrating weak processes, called “Urca” processes, are of fundamental importance in astrophysical environments like (proto-)neutron stars, neutron star mergers, and supernovae. In these environments, matter can reach high temperatures of tens of MeVs and be subject to large magnetic fields. We thus investigate Urca rates at different temperatures and field strengths by performing the full temperature and magnetic-fielddependent rate integrals for different equations of state. We find that the magnetic fields play an important role at temperatures of a few MeV, especially close to or below the direct Urca threshold, which is softened by the magnetic field. At higher temperatures, the effect of the magnetic fields can be overshadowed by the thermal effects. Finally, we observe that the magnetic field influences the neutron decay rates more strongly than the electron-capture rates, leading to a shift in the flavor equilibrium.

electroweak interactions in nuclear 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 OSTI2025

Imprints of High-Density Nuclear Symmetry Energy on Crustal Fraction of Neutron Star Moment of Inertia

The density dependence of nuclear symmetry energy E sym (ρ) remains the most uncertain aspect of the equation of state (EOS) of supradense neutron-rich nucleonic matter. Utilizing an isospin-dependent parameterization of the nuclear EOS, we investigate the implications of the observational crustal fraction of the neutron star (NS) moment of inertia ΔI/I for the E sym (ρ). We find that symmetry energy parameters significantly influence the ΔI/I, while the EOS of symmetric nuclear matter has a negligible effect. In particular, an increase in the slope L and skewness J sym of symmetry energy results in a larger ΔI/I, whereas an increase in the curvature K sym leads to a reduction in ΔI/I. Moreover, the ΔI/I is shown to have the potential for setting a lower limit of symmetry energy at densities exceeding 3 ρ 0 , particularly when L is constrained to values less than 60 MeV, thereby enhancing our understanding of supradense NS matter.

equation of state↗
DOE OSTI2023

Dense nuclear matter equation of state from heavy-ion collisions

The nuclear equation of state (EOS) is at the center of numerous theoretical and experimental efforts in nuclear physics. With advances in microscopic theories for nuclear interactions, the availability of experiments probing nuclear matter under conditions not reached before, endeavors to develop sophisticated and reliable transport simulations to interpret these experiments, and the advent of multi-messenger astronomy, the next decade will bring new opportunities for determining the nuclear matter EOS, elucidating its dependence on density, temperature, and isospin asymmetry. Among controlled terrestrial experiments, collisions of heavy nuclei at intermediate beam energies (from a few tens of MeV/nucleon to about 25 GeV/nucleon in the fixed-target frame) probe the widest ranges of baryon density and temperature, enabling studies of nuclear matter from a few tenths to about 5 times the nuclear saturation density and for temperatures from a few to well above a hundred MeV, respectively. Collisions of neutron-rich isotopes further bring the opportunity to probe effects due to the isospin asymmetry. However, capitalizing on the enormous scientific effort aimed at uncovering the dense nuclear matter EOS, both at RHIC and at FRIB as well as at other international facilities, depends on the continued development of state-of-the-art hadronic transport simulations. Furthermore, this white paper highlights the essential role that heavy-ion collision experiments and hadronic transport simulations play in understanding strong interactions in dense nuclear matter, with an emphasis on how these efforts can be used together with microscopic approaches and neutron star studies to uncover the nuclear EOS.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Strong Kitaev Interaction in BaCo 2⁢ (AsO 4 ) 2

The inelastic neutron scattering results and their analysis unequivocally point to a dominant Kitaev interaction in the honeycomb-lattice cobaltate BaCo 2 ⁢(AsO 4 ) 2 . Our anisotropic-exchange model closely describes all available neutron scattering data in the material’s field-polarized phase. Furthermore, the density-matrix renormalization group results for our model are in close accord with the unusual double-zigzag magnetic order and the low in-plane saturation field of BaCo 2 ⁢(AsO 4 ) 2 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2026

X-Ray Polarization of Z-Type Neutron Star Low-Mass X-Ray Binaries II. Spectropolarimetric Analysis

IXPE has provided for the first time detailed energy- and time-resolved X-ray polarimetry of Z-type neutron star low-mass X-ray binaries (NS-LMXBs) as they move along their color-color diagrams (CCDs). These sources can reach the highest polarization observed for NS-LXMBs in the 2–8 keV range when they move along the horizontal branch. In a previous paper, we characterized the spectral state of a sample of Z-sources using the CCD and estimated the polarization with model-independent analysis. Here, we present detailed spectropolarimetric analysis for each source on each branch using data from IXPE, NICER, and NuSTAR. The continuum X-ray emission of all the sources is well described with a combination of thermal accretion disk emission plus a harder Comptonized component. In addition, reflection features, in particular the relativistically broadened Fe Kα line, are observed for our sources, except GX 5–1. For most of the sources and branches, the main contribution to the X-ray emission and polarization is due to Comptonization: moving from the horizontal branch (HB) to the normal branch (NB), the polarization degree (PD) in the 2–8 keV band varies from about 6% to 3–4%, while the PD is loosely constrained in the flaring branch (FB), due to the shorter exposures. These PD values are significantly higher than theoretical expectations for typical spreading or boundary layer configurations. The polarization of the disk is generally lower (below 3%) but still higher than predictions for an electron scattering-dominated, plane-parallel atmosphere above the disk observed at the corresponding inclination. Moreover, the polarization angle (PA) of the disk seems to be significantly misaligned and not perpendicular to that of Comptonization. We find no correlation between the polarization signal and the inclination, nor with the contribution of reflected photons throughout the Z-track.

accretion↗
NASA NTRS2026

X-Ray Polarization From the Atoll 4U 1735−44 Suggests a Low Inclination

X-ray polarimetry is a new tool capable of probing the geometry of accretion onto weakly magnetized neutron stars. Here we present the first X-ray spectro-polarimetric results from coordinated observations of the atoll source 4U 1735−44, conducted with the Imaging X-ray Polarimetry Explorer (IXPE), NICER, and NuSTAR . Over the 2–8 keV energy range, we obtained a marginal detection of polarization with the polarization degree of 1.4% ± 0.7% and polarization angle of −29◦ ± 14◦, corresponding to a 3σ upper limit on the polarization degree of 3.5%. The best-fit model to describe the spectrum comprises a thermal component associated with the accretion disk, a Comptonized black body component, and a relativistic reflection component. From the reflection model, we infer a disk inclination of ∼40◦. The spectroscopic and polarimetric properties of 4U 1735−44 are consistent with those observed in other atoll sources studied by IXPE, with its low polarization likely due to its low inclination

X-rays: binaries↗
NASA NTRS1962

Some Experiments on the Passage of High-Energy Protons in Dense Matter

A series of experiments designed to study the nuclear cascade resulting from the passage of 1 to 3 GeV protons in matter has been in progress at the Brookhaven Cosmotron. Preliminary results on the fluxes of fast neutrons (upper limits) and of strongly interacting particles above 50 MeV are summarized here (Fig. 1 and Table 1). The four cases studied are: 1-GeV protons on Fe, on chondritic material, and on C 5 H 8 O 2 , and 3-GeV protons on Fe.

Proton Irradiation↗
DOE OSTI2026

Simulation to a Newborn Supernova Remnant from a Low-mass Iron Core Star

Supernova remnant observations show a high degree of asymmetry, mixing, and inhomogeneity. These asymmetries are seeded during the early seconds of the explosion and are further enhanced and modified as the shock and ejecta move through the stellar progenitor and into the circumstellar medium. We present simulations of a 9.6 M⊙ zero-metallicity progenitor initialized after shock revival and evolved for several years when the ejecta is in the circumstellar medium. A suite of 1D and 2D simulations examines the effects of neutron-star wind and radioactive decay heating. In 1D, decay heating forms a low-density bubble that suppresses the reverse shock. While in 2D, the heating is localized to metal-rich pockets, inflating them and compressing the surrounding material into dense shells. In 3D, the neutron-star wind and decay heating modify the plume morphology, producing more large-scale structures. The extended plume morphology leads to an asymmetrical shock breakout. After breakout, the leading plumes cannot keep up with the shock front, resulting in deceleration and fragmentation by the reverse shock while retaining the large-scale asymmetry. The projected ejecta morphology and velocities are strongly viewing angle dependent. The relatively uniform metal-rich distribution does not resemble the strongly inhomogeneous ejecta structure of Cas A. The 160-isotope decay network shows that 24.4% of the radioactive heating comes from decay chains other than the canonical 56Ni chain. The low explosion energy, low 56Ni yield, and Ni/Fe ratio greater than unity suggest an observational signature similar to an electron capture supernova.

Neopane, Sudarshan [University of Tennessee (UT)]↗
NASA NTRS2020

The Slow Demise of the Long-Lived Sn 2005ip

The Type IIn supernova (SN IIn) 2005ip is one of the most well-studied and long-lasting examples of an SN interacting with its circumstellar environment. The optical light curve plateaued at a nearly constant level for more than five years, suggesting ongoing shock interaction with an extended and clumpy circumstellar medium (CSM). Here, we present continued observations of the SN from ∼1000 to 5000 d post-explosion at all wavelengths, including X-ray, ultraviolet, near-infrared (NIR), and mid-infrared. The UV spectra probe the pre-explosion mass loss and show evidence for CNO processing. From the bolometric light curve, we find that the total radiated energy is in excess of 10 50 erg, the progenitor star’s pre-explosion mass-loss rate was >~ 1 × 10 −2 Mꙩ yr −1 , and the total mass lost shortly before explosion was >~1Mꙩ, though the mass lost could have been considerably larger depending on the efficiency for the conversion of kinetic energy to radiation. The ultraviolet through NIR spectrum is characterized by two high-density components, one with narrow high-ionization lines, and one with broader low-ionization H I, He I, [OI], Mg II, and Fe II lines. The rich Fe II spectrum is strongly affected by Lyα fluorescence, consistent with spectral modelling. Both the Balmer and He I lines indicate a decreasing CSM density during the late interaction period. We find similarities to SN 1988Z, which shows a comparable change in spectrum at around the same time during its very slow decline. These results suggest that, at long last, the shock interaction in SN 2005ip may finally be on the decline.

stars↗
DOE OSTI2026

Collapse of magnetized white dwarfs as site of heavy-element formation and kilonova signal

We present the first end-to-end calculation connecting the accretion-induced collapse (AIC) of a magnetized, rapidly rotating white dwarf to observable kilonova signatures, combining two-dimensional (2D) general-relativistic neutrino-magnetohydrodynamic simulations, followed by radiation hydrodynamics with in-situ nuclear network and 2D Monte Carlo radiative transfer with spatially resolved heating rates. Unlike all previous unmagnetized AIC models – which predicted proton-rich, $^{56}$Ni-dominated ejecta – strong magnetic fields eject ${\approx }\, 0.2\, \mathrm{ M}_\odot$ of neutron-rich material ($\langle Y_e \rangle \sim 0.24$) on dynamical time-scales, before neutrino irradiation can raise the electron fraction, enabling strong r-process nucleosynthesis up to and beyond the third peak. The resulting kilonova is lanthanide-rich ($X_{\rm lan} \approx 8~{{\ \rm per\ cent}}$) and dominated by near-infrared emission. We compute synthetic light curves in the Large Synoptic Survey Telescope and J ames Webb Space Telescope bands and find striking agreement, without parameter tuning, between the observations of AT 2023vfi/GRB 230307A and our broadband light curves for polar viewing angles. These results establish magnetized AIC as a viable channel for heavy r-process element production and a compelling progenitor candidate for long-duration gamma-ray bursts with kilonova signatures.

MHD↗
DOE OSTI2023

TRTR NRAD NRS Beamline Abstract

The Neutron Radiography (NRAD) Reactor is a 250kW TRIGA housed under the largest hot cell in the United States, making it the foremost location to perform neutron imaging of irradiated nuclear fuels and materials, including those intended for use with advanced reactors. Currently, NRAD has two radial beamlines that are used for neutron radiography and tomography. Upgrades to the North neutron beamline include the replacement of the in-tank beam chamber, through-the-wall collimator, and neutron shutter. These beamline upgrades will not only improve image quality for current capabilities but will also condition the beam to be more suitable for advanced methods such as neutron powder diffraction. These upgrades will increase the excess reactivity of the core, reduce unnecessary activation and exposure to workers, and significantly reduce the amount of shielding required around the beam. This presentation will describe how these upgrades improve beam quality, increase utilization of the reactor, and reduce radiation exposure to personnel.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗
NASA NTRSDate not supplied

Lunar Nuclear Reactor Neutron Fluence and Gamma Dose Estimates

NASA aims to deploy a fission reactor on the lunar surface by 2030 to generate 100 kWe of power for at least 10 years without maintenance or refueling. To assess system and subsystem survivability according to those requirements, a model of the reactor pallet was created in the radiation transport code MCNP to determine the neutron fluences and gamma doses that the materials and systems would be exposed to. The model was based on the 2025 government design. The model indicated that the control drum motors, Brayton cycle engines, and pipes and valves would be exposed to neutron fluences of 1×10 15 – 1×10 16 n/cm 2 per year and gamma doses ranging from 10 to 500 Mrad per year. The electronics box housing the instrumentation and control elements is expected to receive 1×10 14 – 1×10 15 n/cm 2 per year and gamma doses in the range of 1–5 Mrad per year. These doses necessitate that the control drum motors, Brayton cycle engines, valve components, control electronics, and other components be radiation-hardened (>10 16 n/cm 2 and >10 Mrad gamma dose) or additionally shielded to survive the 10-year lunar mission.

Radiation-Hardened Materials↗
DOE OSTI1969

Application of Nuclear Technology to Art Identification Problems: First Annual Report

Throughout the centuries, as man has become increasingly affluent, his interest in antiquities and objects of art, and the price he has been willing to pay for them, has increased. Inevitably, one result has been to attract the unscrupulous to the production and sale of counterfeit paintings, sculptures, and other works of art. The skill of forgers varies but often is highly sophisticated. However, dealers, museum curators, and other experts have also become more skillful in the detection of forgeries. Scientific tools of increasing sensitivity and sophistication have gradually been applied to the examination of the materials of art and archaeology. Such tools frequently support and render more reliable the judgment of the experts. The quality of forgeries may improve further as their makers become acquainted with the new methods of detection and in turn learn to circumvent or confound the methods. Therefore, the development of still more advanced methods of examination has great utility in the art world. This is particularly true if the ultimate effect is to make circumvention of the methods of examination so costly that the economic incentives for producing forgeries may be substantially reduced, if not eliminated.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Integrated Neutronics Modeling for Inertial Fusion Energy Systems: Development and Application to LD-FIRST

Lawrence Livermore National Laboratory (LLNL) is proposing a new Laser Driven Fusion Integration Research and Science Test Facility (LD-FIRST) with the goal of providing an experimental testbed for future Inertial Fusion Energy (IFE) systems. However, IFE systems require detailed and accurate multiphysics modeling to quantify material damage, thermal loading, and tritium breeding within complex chamber environments. This article presents the first step in an integrated multiphysics framework that couples meshed CAD-based geometry within Monte Carlo neutronic simulations to enable high-fidelity analysis of IFE chamber concepts, with future coupling to external codes. The neutronics workflow utilizes OpenMC and its third-party capability to use CAD-based geometries through DAGMC and tally on unstructured meshes with Libmesh to evaluate neutron transport behavior, geometric fidelity, and material performance under reactor-relevant conditions. The use of tailored tallies on unstructured meshes in this framework allows direct transfer without interpolating to CFD simulation tools. Two IFE chambers were evaluated, both conceived by LLNL: HYLIFE-II and Laser IFE (LIFE). This work produced high-fidelity conformal surface and volumetric meshes of the HYLIFE-II and LIFE chambers with mapped spatial insight into material damage, thermal loading, and tritium breeding. The HYLIFE-II model was built utilizing available resources and used as a test case to verify that the neutronics framework can handle complex geometries. The LIFE chamber CAD was provided by LLNL and was the main focus of this work. This work analyzes multiple ternary alloy breeding materials for the LIFE chamber, across different 6 Li enrichments to produce data relevant to the LD-FIRST project. This work also investigates the level of model fidelity for the LIFE chamber, and results show that inclusion of detailed first wall and coolant structures increased the predicted tritium breeding ratio (TBR) by ~30%, highlighting the sensitivity of tritium breeding and the need for a high-fidelity simulation framework for IFE chambers. These developments provide a scalable toolset for the design and optimization of next-generation IFE chambers, forming a solid foundation for future coupled multiphysics analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRSDate not supplied

High-Fidelity Modeling Approach for Coupling Thermal, Fluid, and Neutronic Analysis in an NTP Engine

Nuclear thermal propulsion (NTP), with its increased efficiency compared to traditional rocket engines, is a potentially enabling technology for space missions of higher complexity, such as delivering heavy machinery or crew to the moon or other planets. Thus, exploring modeling approaches that can effectively couple the thermal, fluid, and neutronic behaviors of NTP reactors is essential to mitigate cost and schedule challenges with hardware development. This paper details one such approach which utilizes Cardinal to unify reactor physics analysis performed using OpenMC, heat transfer analysis performed using the MOOSE framework, and thermal hydraulic analysis performed using the MOOSE Thermal Hydraulic Module. This approach is demonstrated with a proof of concept that focuses on a legacy reactor geometry from the historic NERVA program. This preliminary model found an overprediction in, when compared to the reference data, axial power peaking and the temperatures of both solids and fluids by roughly 20%. It is clear, however, that the model conserved a number of broad trends, and a few model refinements should reduce the observed error. This approach has definite merit and should enable the production of informative results for more modern NTP systems. This work is a first step in creating a transient reactor model that can be updated across a sequence of time steps based on continually fluctuating, time-dependent boundary conditions.

Nuclear Thermal Propulsion↗
DOE OSTI2026

Shakeup and shakeoff spectra in the electron capture decay of atomic 7 Be

The most stringent laboratory-based experimental limits on the existence of submegaelectronvolt sterile neutrinos are currently set by decay spectroscopy of radioactive 7 Be embedded into superconducting sensors. The systematic uncertainties are dominated by the modeling of the electron shakeup and shakeoff spectra that are not based on state-of-the-art atomic theory and do not include electron correlations or relativistic effects. We have used the multiconfiguration Dirac-Fock formalism to obtain correlated wave functions ab initio and compute all single and double shake processes in the electron capture decay of atomic 7 Be . The simulations can explain some but not all of the observed spectral features, likely because the wave functions are modified by the Ta sensor material that the 7 Be is embedded into. The new models also show that the 𝐿/𝐾 electron capture ratio of 7 Be in Ta has previously been slightly underestimated revising the previous value of 0.070(7) to a new value of 0.0756(20).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Emergence of low-energy spin waves in superconducting electron-doped cuprates

In order to fully utilize the technological potential of unconventional superconductors, an enhanced understanding of the superconducting mechanism is necessary. In the best performing superconductors, the cuprates, superconductivity is intimately linked with magnetism, although the details of this coupling remain elusive. Here, we address this gap by studying the electron-doped cuprate Nd 1.85 Ce 0.15 CuO 4−δ that has an antiferromagnetic ground state when synthesized and only becomes superconducting after a reductive annealing process. Using neutron spectroscopy, we show that the as-grown crystal exhibits a large spin pseudogap in the magnetic fluctuation spectrum. Annealing removes defects introduced by the commonly employed synthesis method and significantly reduces the spin pseudogap. While the spin pseudogap in the annealed sample likely arises from superconductivity, in the as-grown sample it results from the absence of long-wavelength spin waves. These results reveal a direct connection between defects, magnetism, and superconductivity, offering new insight into the mechanisms underlying high-temperature superconductivity and guiding the design of improved superconducting materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗