Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nuclear Equation of State”

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 OSTI2021

Analysis of the neutron matter equation of state and the symmetry energy up to fourth order of chiral effective field theory

We present predictions for the neutron matter equation of state, from leading to fourth order of chiral effective field theory, using recently developed, accurate chiral nucleon-nucleon potentials. For the many-body method, we employ the nonperturbative particle-particle ladder approximation, that is, we solve the G-matrix equation. Furthermore, we find the impact of subleading three-neutron forces to be mild and attractive. We also show order-by-order predictions for the symmetry energy, and discuss its density dependence in relation to empirical constraints. For the nuclear matter equation of state, in this work we adopt an empirical parametrization with good saturation properties. This is to highlight, specifically, the energy and pressure in neutron matter, particularly when comparing with empirical constraints.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Bayesian mixture model approach to quantifying the empirical nuclear saturation point

The equation of state (EOS) in the limit of infinite symmetric nuclear matter exhibits an equilibrium density, $n_0 \approx 0.16 \, \mathrm{fm}^{-3}$, at which the pressure vanishes and the energy per particle attains its minimum, $E_0 \approx -16 \, \mathrm{MeV}$. Although not directly measurable, the nuclear saturation point $(n_0,E_0)$ can be extrapolated by density functional theory (DFT), providing tight constraints for microscopic interactions derived from chiral effective field theory (EFT). However, when considering several DFT predictions for $(n_0,E_0)$ from Skyrme and Relativistic Mean Field (RMF) models together, a discrepancy between these model classes emerges at high confidence levels that each model prediction's uncertainty cannot explain. How can we leverage these DFT constraints to rigorously benchmark nuclear saturation properties of chiral interactions? To address this question, we present a Bayesian mixture model that combines multiple DFT predictions for $(n_0,E_0)$ using an efficient conjugate prior approach. The inferred posterior distribution for the saturation point's mean and covariance matrix follows a Normal-inverse-Wishart class, resulting in posterior predictives in the form of correlated, bivariate $t$-distributions. The DFT uncertainty reports are then used to mix these posteriors using an ordinary Monte Carlo approach. At the 95\% credibility level, we estimate $n_0 \approx 0.157 \pm 0.010 \, \mathrm{fm}^{-3}$ and $E_0 \approx -15.97 \pm 0.40 \, \mathrm{MeV}$ for the marginal (univariate) $t$-distributions. Combined with chiral EFT calculations of the pure neutron matter EOS, we obtain bivariate normal distributions for the nuclear symmetry energy and its slope parameter evaluated at $n_0$: $S_v \approx 32.0 \pm 1.1 \, \mathrm{MeV}$ and $L\approx 52.6\pm 8.1 \, \mathrm{MeV}$ (95\%), respectively. Furthermore, our Bayesian framework is publicly available, so practitioners can readily use and extend our results.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Influence of helium bubbles on the bulk equation of state of gold under static compression

Nuclear materials often evolve into two-phase systems comprising a bulk matrix with dispersed inert-gas bubbles. The presence of these bubbles can have consequences to the thermomechanical response of materials and is a key life-limiting factor in some nuclear fuel forms. Understanding the behavior of these two-phase, bubble-matrix systems is, thus, important to improved predictive models and frameworks for many nuclear materials applications. While temperature excursions of these two-phase systems have been characterized, fewer studies have focused on the evolution of inert-gas bubbles under pressure. Here, in this paper, we use x-ray tools to interrogate a He-implanted gold foil to determine the pressure-dependent evolution of the individual components (Au matrix + bubbles), and we compare that total pressure dependence to theoretical equation-of-state descriptions based on mixing rules.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
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 OSTI2021

Neutron Skin Measurement of208Pb Using Parity-Violating Electron Scattering

Nuclear stability requires large nuclei to have more neutrons than protons. While the central core of the nucleus is composed of symmetric nuclear matter, the extra neutrons in heavy nuclei are pushed out to the surface of the nucleus, forming a pure outer neutron skin layer. While the proton distribution inside a nucleus has been accurately measured using electromagnetic probes, previous measurements of neutron distributions in complex nuclei have suffered a lack of systematic precision. PREX-2 and CREX use Parity-Violating electron Scattering (PVeS), as an electroweak interaction probe, to measure the neutron distributions in208Pb and48Ca, respectively. Making use of elastic scattering of longitudinally polarized electrons from unpolarized (isotopically pure) targets, PREX-2 and CREX measure the asymmetry in scattering rates between opposite beam polarization states. Since both experiments share the same exact measurement concept and apparatus, much of this dissertation applies to both PREX-2 and CREX. However, the presented asymmetry analysis and neutron skin extraction are only for PREX-2. PREX-2 measures the parity-violating asymmetry, APV, for 953 MeV electrons scattered elastically from208Pb at?5°in the lab. The final result is APV= 550.00±16.09(stat.)±8.16(syst.)ppb. From this measurement, we report the neutron skin thickness in208Pb nuclei to be, Rskin= 0.278±0.078(exp.)±0.012(model) fm. This measurement has broad implications throughout the physics community, increasing our knowledge in neutron star structure, the equation of state of nuclear matter, nuclear baryon density, nuclear electric dipole polarizability, and more. Key Words: Neutron skin, electroweak interaction, PREX-2, CREX, parity-violating asymmetry, neutron star, quartz detector, weak charge radius, symmetry energy, Hall A, Jefferson Lab

Adhikari, Devi↗
NASA NTRS2016

The NuSTAR View of the Non-Thermal Emission from PSR J0437-4715

We present a hard X-ray Nuclear Spectroscopic Telescope Array (NuSTAR) observation of PSR J0437-4715, the nearest millisecond pulsar. The known pulsations at the apparent pulse period approximately 5.76 ms are observed with a significance of 3.7sigma, at energies up to 20 keV above which the NuSTAR background dominates. We measure a photon index gamma = 1.50 +/- 0.25(90 per cent confidence) for the power-law fit to the non-thermal emission. It had been shown that spectral models with two or three thermal components fit the XMM-Newton spectrum of PSR J0437-4715, depending on the slope of the power-law component, and the amount of absorption of soft X-rays. The new constraint on the high-energy emission provided by NuSTAR removes ambiguities regarding the thermal components of the emission below 3 keV. We performed a simultaneous spectral analysis of the XMM-Newton and NuSTAR data to confirm that three thermal components and a power law are required to fit the 0.3-20 keV emission of PSR J0437-4715. Adding a ROSAT-PSPC spectrum further confirmed this result and allowed us to better constrain the temperatures of the three thermal components. A phase resolved analysis of the NuSTAR data revealed no significant change in the photon index of the high-energy emission. This NuSTAR observation provides further impetus for future observations with the NICER mission (Neutron Star Interior Composition Explorer) whose sensitivity will provide much stricter constraints on the equation of state of nuclear matter by combining model fits to the pulsars phase-folded light curve with the pulsars well-defined mass and distance from radio timing observations.

stars; neutro↗
DOE OSTI2023

Secular outflows from 3D MHD hypermassive neutron star accretion disc systems

Magnetized hypermassive neutron stars (HMNSs) have been proposed as a way for neutron star mergers to produce high electron fraction, high-velocity ejecta, as required by kilonova models to explain the observed light curve of GW170817. The HMNS drives outflows through neutrino energy deposition and mechanical oscillations, and raises the electron fraction of outflows through neutrino interactions before collapsing to a black hole (BH). Here, in this study, we perform 3D numerical simulations of HMNS–torus systems in ideal magnetohydrodynamics, using a leakage/absorption scheme for neutrino transport, the nuclear APR equation of state, and Newtonian self-gravity, with a pseudo-Newtonian potential added after BH formation. Due to the uncertainty in the HMNS collapse time, we choose two different parametrized times to induce collapse. We also explore two initial magnetic field geometries in the torus, and evolve the systems until the outflows diminish significantly ($\sim\!\! 1\!\! - \!\!2\ \mathrm{s}$). We find bluer, faster outflows as compared to equivalent BH–torus systems, producing M ∼ 10 −3 M ⊙ of ejecta with Y e ≥ 0.25 and v ≥ 0.25c by the simulation end. Approximately half the outflows are launched in disc winds at times $t\lesssim 500 \ \mathrm{ms}$, with a broad distribution of electron fractions and velocities, depending on the initial condition. The remaining outflows are thermally driven, characterized by lower velocities and electron fractions. Nucleosynthesis with tracer particles shows patterns resembling solar abundances in all models. Although outflows from our simulations do not match those inferred from two-component modelling of the GW170817 kilonova, self-consistent multidimensional detailed kilonova models are required to determine whether our outflows can power the blue kilonova.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Horizon radiation reaction forces

Using Effective Field Theory (EFT) methods, we compute the effects of horizon dissipation on the gravitational interactions of relativistic binary black hole systems. We assume that the dynamics is perturbative, i.e it admits an expansion in powers of Newton’s constant (post-Minkowskian, or PM, approximation). As applications, we compute corrections to the scattering angle in a black hole collision due to dissipative effects to leading PM order, as well as the post-Newtonian (PN) corrections to the equations of motion of binary black holes in non-relativistic orbits, which represents the leading order finite size effect in the equations of motion. The methods developed here are also applicable to the case of more general compact objects, eg. neutron stars, where the magnitude of the dissipative effects depends on non-gravitational physics (e.g, the equation of state for nuclear matter).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Efficient learning of accurate surrogates for simulations of complex systems

Machine learning methods are increasingly deployed to construct surrogate models for complex physical systems at a reduced computational cost. However, the predictive capability of these surrogates degrades in the presence of noisy, sparse or dynamic data. Here, we introduce an online learning method empowered by optimizer-driven sampling that has two advantages over current approaches: it ensures that all local extrema (including endpoints) of the model response surface are included in the training data, and it employs a continuous validation and update process in which surrogates undergo retraining when their performance falls below a validity threshold. We find, using benchmark functions, that optimizer-directed sampling generally outperforms traditional sampling methods in terms of accuracy around local extrema even when the scoring metric is biased towards assessing overall accuracy. Finally, the application to dense nuclear matter demonstrates that highly accurate surrogates for a nuclear equation-of-state model can be reliably autogenerated from expensive calculations using few model evaluations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Magnetars and axion-like particles: probes with the hard X-ray spectrum

Quiescent hard X-ray and soft gamma-ray emission from neutron stars constitute a promising frontier to explore axion-like-particles (ALPs). ALP production in the core peaks at energies of a few keV to a few hundreds of keV; subsequently, the ALPs escape and convert to photons in the magnetosphere. The emissivity goes as ~T 6 while the conversion probability is enhanced for large magnetic fields, making magnetars, with their high core temperatures and strong magnetic fields, ideal targets for probing ALPs. We compute the energy spectrum of photons resulting from conversion of ALPs in the magnetosphere and then compare it against hard X-ray data from NuSTAR, INTEGRAL, and XMM-Newton, for a set of eight magnetars for which such data exists. Upper limits are placed on the product of the ALP-nucleon and ALP-photon couplings. For the production in the core, we perform a calculation of the ALP emissivity in degenerate nuclear matter modeled by a relativistic mean field theory. The reduction of the emissivity due to improvements to the one-pion exchange approximation is incorporated, as is the suppression of the emissivity due to proton superfluidity in the neutron star core. A range of core temperatures is considered, corresponding to different models of the steady heat transfer from the core to the stellar surface. For the subsequent conversion, we solve the coupled differential equations mixing ALPs and photons in the magnetosphere. The conversion occurs due to a competition between the dipolar magnetic field and the photon refractive index induced by the external magnetic field. Semi-analytic expressions are provided alongside the full numerical results. Furthermore, we also present an analysis of the uncertainty on the axion limits we derive due to the uncertainties in the magnetar masses, nuclear matter equation of state, and the proton superfluid critical temperature.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Intrinsic correlations among characteristics of neutron-rich matter imposed by the unbound nature of pure neutron matter

The unbound nature of pure neutron matter (PNM) requires intrinsic correlations between the symmetric nuclear matter (SNM) EOS parameters (incompressibility K 0 , skewness J 0 and kurtosis I 0 ) and those (slope L, curvature K sym and skewness J sym ) characterizing the symmetry energy independent of any nuclear many-body theory. We investigate these intrinsic correlations and their applications in better constraining the poorly known high-density behavior of nuclear symmetry energy. Several novel correlations connecting the characteristics of SNM EOS with those of nuclear symmetry energy are found. In particular, at the lowest-order of approximations, the bulk parts of the slope L, curvature K sym and skewness J sym of the symmetry energy are found to be L ≈ K 0 /3, K sym ≈ LJ 0 /2K 0 and J sym ≈ I 0 L/3K 0 , respectively. High-order corrections to these simple relations can be written in terms of the small ratios of high-order EOS parameters. The resulting intrinsic correlations among some of the EOS parameters reproduce very nicely their relations predicted by various microscopic nuclear many-body theories and phenomenological models constrained by available data of terrestrial experiments and astrophysical observations in the literature. Here, the unbound nature of PNM is fundamental and the required intrinsic correlations among the EOS parameters characterizing both the SNM EOS and symmetry energy are universal. These intrinsic correlations provide a novel and model-independent tool not only for consistency checks but also for investigating the poorly known high-density properties of neutron-rich matter by using those with smaller uncertainties.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Nucleon-nucleon potentials from Δ-full chiral effective-field-theory and implications

We closely investigate NN potentials based upon the Δ-full version of chiral effective-field-theory. We find that recently constructed NN potentials of this kind, which (when applied together with three-nucleon forces) were presented as predicting accurate binding energies and radii for a range of nuclei from A = 16 to A = 132 and providing accurate equations of state for nuclear matter, yield a χ 2 /datum of 60 for the reproduction of the pp data below 100 MeV laboratory energy. This χ 2 is more than three times what the Hamada-Johnston potential of the year of 1962 achieved already some 60 years ago. We perceive this historical fact as concerning in view of the current emphasis on precision. We are able to trace the very large χ 2 as well as the apparent success of the potentials in nuclear structure to unrealistic predictions for P-wave states, in which the Δ-full next-to-next-to-leading order (NNLO) potentials are off by up to 40 times the NNLO truncation errors. In fact, we show that the worse the description of the P-wave states, the better the predictions in nuclear structure. Thus, these potentials cannot be seen as the solution to the outstanding problems in current microscopic nuclear structure physics.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Accurate Determination of the Neutron Skin Thickness of Pb 208 through Parity-Violation in Electron Scattering

We report a precision measurement of the parity-violating asymmetry A PV in the elastic scattering of longitudinally polarized electrons from 208 Pb. We measure A PV = 550 ± 16 (stat) ±8 (syst) parts per billion, leading to an extraction of the neutral weak form factor F W (Q 2 = 0.00616 GeV 2 ) = 0.368 ± 0.013. Combined with our previous measurement, the extracted neutron skin thickness is R n -R p = 0.283 ± 0.071 fm. The result also yields the first significant direct measurement of the interior weak density of 208 Pb: ρ$^0_W$ = -0.0796 ± 0.0036(exp) ± 0.0013(theo) fm -3 leading to the interior baryon density ρ$^0_b$ = 0.1480 ± 0.0036(exp) ± 0.0013(theo) fm -3 . Finally, the measurement accurately constrains the density dependence of the symmetry energy of nuclear matter near saturation density, with implications for the size and composition of neutron stars.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Charge Radius of Neutron-Deficient Ni 54 and Symmetry Energy Constraints Using the Difference in Mirror Pair Charge Radii

he nuclear root-mean-square charge radius of Ni 54 was determined with collinear laser spectroscopy to be R ( Ni 54 ) = 3.737 ( 3 ) fm . In conjunction with the known radius of the mirror nucleus Fe 54 , the difference of the charge radii was extracted as Δ R ch = 0.049 ( 4 ) fm . Based on the correlation between Δ R ch and the slope of the symmetry energy at nuclear saturation density ( L ), we deduced 21 ≤ L ≤ 88 MeV . The present result is consistent with the L from the binary neutron star merger GW170817, favoring a soft neutron matter EOS, and barely consistent with the PREX-2 result within 1 σ error bands. Our result indicates the neutron-skin thickness of Ca 48 as 0.15–0.21 fm.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Probing the Incompressibility of Nuclear Matter at Ultrahigh Density through the Prompt Collapse of Asymmetric Neutron Star Binaries

Using 250 neutron star merger simulations with microphysics, we explore for the first time the role of nuclear incompressibility in the prompt collapse threshold for binaries with different mass ratios. We demonstrate that observations of prompt collapse thresholds, either from binaries with two different mass ratios or with one mass ratio but combined with the knowledge of the maximum neutron star mass or compactness, will constrain the incompressibility at the maximum neutron star density K max to within tens of percent. This otherwise inaccessible measure of K max can potentially reveal the presence of hyperons or quarks inside neutron stars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Advancements in High-Q development for novel medium-velocity 644 MHz 5-cell elliptical superconducting RF cavities for continuous-wave operation in heavy-ion linacs

The Facility for Rare Isotope Beams is a world-leading center for experimental nuclear physics research, and relies on a first-in-kind superconducting RF driver linac to supply 400 kW beams of a uniquely large range of particles, from protons through the heaviest uranium ions. The FRIB400 project proposes to double the end energy of the current superconducting FRIB linac for the heaviest uranium ions from 200 MeV per nucleon (MeV/u) to 400 MeV/u, which equates to approximately 1 GeV for protons. The increased rare isotope production from higher-energy drive beams would deliver new and exciting capabilities to FRIB users, which include extending the facility's reach along the neutron drip-line, greatly increasing the facility's rare isotope yield, and significantly improving the precision of parameter measurements for the nuclear matter equation of state, which is a particularly timely complement to the recent advent of multi-messenger astronomy, and the study of neutron star mergers. The unique set of design parameters defined by the 400 MeV/u (uranium) energy goal, the physical size, continuous-wave operation, and cryogenic capacity of the extant facilities at FRIB, requires development of a first-in-kind beta = v/c = 0.65, 644 MHz 5-cell elliptical superconducting RF (SRF) cavity, with field-leading efficiency, measured by the cavity's intrinsic quality factor, Q_0. The field of SRF has recently experienced a very active period of development, which demonstrated the potential to raise the Q_0 of 1.3 GHz cavities by factors of two to three or more. In addition to uncovering fascinating new properties of the fundamental physics of superconducting RF phenomena, these improvements suggest substantial improvements could be made over the minimum Q_0 of the proposed FRIB400 upgrade SRF cavities, if these new techniques can be adapted to the beta = 0.65, 644 MHz regime. This work presents the first validation of the the novel beta = 0.65, 644 MHz 5-cell elliptical superconducting RF cavity design, and the first results of both conventional and advanced RF surface preparation techniques. In particular, this work focuses on two recently-developed high-Q_0 recipes: N-doping, and furnace baking, setting the current world-record Q_0 for this type of resonant cavity at 3.8 x 10^10 at the FRIB400 operating gradient of 17.5 MV/m. We also find evidence that these results may be further improved upon, and make the case for further RF surface preparation recipe optimization. We then pair these high-power RF investigations with techniques borrowed from material science to try to better understand the ways in which specific properties of the niobium material, such as grain structure and superconducting flux pinning force, can be related to its superconducting RF performance. These experiments better elucidate the effects of high-temperature annealing on niobium samples, which suggest methods of further optimizing cavity flux expulsion and thus SRF performance.Read

McGee, Kellen E.↗
DOE OSTI2025

Illuminating Dark Matter Admixed in Neutron Stars with Simultaneous Mass–Radius Constraints

We investigate how simultaneous mass and radius measurements of massive neutron stars can help constrain the properties of dark matter possibly admixed in them. Within a fermionic dark matter model that interacts only through gravitation, along with a well-constrained nuclear matter equation of state, we show that the simultaneous mass and radius measurement of PSRJ0740+6620 reduces the uncertainty of dark matter central energy density by more than 50% compared to the results obtained from using the two observables independently, while other dark matter parameters remain unconstrained. Additionally, we find that the dark matter fraction fD should be smaller than 2% when constrained by the observed neutron star maximum mass alone, and it could be even smaller than 0.3% with the simultaneous measurement of mass and radius, supporting the conclusion that only a small amount of dark matter exists in dark matter admixed neutron stars (DANSs).

NICER↗