SEARCH · Engineering Papers
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.
The nuclear Thomas-Fermi model which is based on nuclear matter calculations has been successfully applied to the study of the bulk properties of nuclei. It is ideally suited for extrapolation into the region of very neutron-rich and of superheavy nuclei. It is therefore a valuable approach for r-process calculations as well as for the study of neutron star matter at subnuclear densities.
Here, this paper compiles the model parameters and zero-temperature properties of an extensive collection of published theoretical nuclear interactions, including 255 nonrelativistic (Skyrme-like) forces, 270 relativistic mean field (RMF) and point-coupling (RMF-PC) forces, and 13 Gogny-like forces. This forms the most exhaustive tabulation of model parameters to date. The properties of uniform symmetric matter and pure neutron matter at the saturation density are determined. Symmetry properties found from the second-order term of a Taylor expansion in neutron excess are compared with the energy difference of pure neutron and symmetric nuclear matter at the saturation density. Selected liquid-droplet model parameters, including the surface tension and surface symmetry energy, are determined for semi-infinite surfaces. Theoretical liquid droplet model neutron skin thicknesses of the neutron-rich closed-shell nuclei 48 Ca and 208 Pb are compared with published theoretical Hartree-Fock and experimental results. A similar comparison is made between theoretical liquid-droplet and experimental values of dipole polarizabilities. In addition, radii, binding energies, moments of inertia and tidal deformabilities of 1.2M ⊙ , 1.4M ⊙ and 1.6M ⊙ neutron stars are computed. An extensive correlation analysis of bulk matter, nuclear structure, and low-mass neutron star properties is performed and compared with nuclear experiments and astrophysical observations.
We show that, for physically reasonable bulk and surface properties, the lowest energy state of dense matter consists of quark matter coexisting with nuclear matter in the presence of an essentially uniform background of electrons. We estimate the size and nature of spatial structure in this phase, and show that at the lowest densities the quark matter forms droplets embedded in nuclear matter, whereas at higher densities it can exhibit a variety of different topologies. A finite fraction of the interior of neutron stars could consist of matter in this new phase, which would provide new mechanisms for glitches and cooling.
We construct efficient emulators for the computation of the infinite nuclear matter equation of state. These emulators are based on the subspace-projected coupled-cluster method for which we here develop a new algorithm called small-batch voting to eliminate spurious states that might appear when emulating quantum many-body methods based on a non-Hermitian Hamiltonian. The efficiency and accuracy of these emulators facilitate a rigorous statistical analysis within which we explore nuclear matter predictions for > 10 6 different parametrizations of a chiral interaction model with explicit Δ -isobars at next-to-next-to leading order. Constrained by nucleon-nucleon scattering phase shifts and bound-state observables of light nuclei up to He 4 , we use history matching to identify nonimplausible domains for the low-energy coupling constants of the chiral interaction. Within these domains we perform a Bayesian analysis using sampling and importance resampling with different likelihood calibrations and study correlations between interaction parameters, calibration observables in light nuclei, and nuclear matter saturation properties. Published by the American Physical Society 2024
Strongly-interacting dark matter can be accumulated in large quantities inside the Earth, and for dark matter particles in a few GeV mass range, it can exist in large quantities near the Earth’s surface. We investigate the constraints imposed on such dark matter properties by its upscattering by fast neutrons in nuclear reactors with subsequent scattering in nearby well-shielded dark matter detectors, schemes which are already used for searches of the coherent reactor neutrino scattering. We find that the existing experiments cover new parameter space on the spin-dependent interaction between dark matter and the nucleon. Similar experiments performed with research reactors, and lesser amount of shielding, may provide additional sensitivity to strongly-interacting dark matter.
Properties of dense matter in nuclear statistical equilibrium analyzed at high densities and temperatures in connection with stellar evolution
Absorption of fermionic dark matter leads to a range of distinct and novel signatures at dark matter direct detection and neutrino experiments. We study the possible signals from fermionic absorption by nuclear targets, which we divide into two classes of four Fermi operators: neutral and charged current. In the neutral current signal, dark matter is absorbed by a target nucleus and a neutrino is emitted. This results in a characteristically different nuclear recoil energy spectrum from that of elastic scattering. The charged current channel leads to induced $\beta$ decays in isotopes which are stable in vacuum as well as shifts of the kinematic endpoint of $ \beta$ spectra in unstable isotopes. To confirm the possibility of observing these signals in light of other constraints, we introduce UV completions of example higher dimensional operators that lead to fermionic absorption signals and study their phenomenology. Most prominently, dark matter which exhibits fermionic absorption signals is necessarily unstable leading to stringent bounds from indirect detection searches. Nevertheless, we find a large viable parameter space in which dark matter is sufficiently long lived and detectable in current and future experiments.
In the past few years, new observations of neutron stars (NSs) and NS mergers have provided a wealth of data that allow one to constrain the equation of state (EOS) of nuclear matter at densities above nuclear saturation density. However, most observations were based on NSs with masses of about 1.4 M ⊙ , probing densities up to ~three to four times the nuclear saturation density. Even higher densities are probed inside massive NSs such as PSR J0740+6620. Very recently, new radio observations provided an update to the mass estimate for PSR J0740+6620, and X-ray observations by the NICER and XMM telescopes constrained its radius. Based on these new measurements, we revisit our previous nuclear physics multimessenger astrophysics constraints and derive updated constraints on the EOS describing the NS interior. By combining astrophysical observations of two radio pulsars, two NICER measurements, the two gravitational-wave detections GW170817 and GW190425, detailed modeling of the kilonova AT 2017gfo, and the gamma-ray burst GRB 170817A, we are able to estimate the radius of a typical 1.4 M ⊙ NS to be ${11.94}_{-0.87}^{+0.76}\,\mathrm{km}$ at 90% confidence. Finally, our analysis allows us to revisit the upper bound on the maximum mass of NSs and disfavors the presence of a strong first-order phase transition from nuclear matter to exotic forms of matter, such as quark matter, inside NSs.
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
The core collapse of a massive star results in the formation of a proto-neutron star (PNS). If enough material is accreted onto a PNS, it will become gravitationally unstable and further collapse into a black hole (BH). We perform a systematic study of failing core-collapse supernovae in spherical symmetry for a wide range of pre-supernova progenitor stars and equations of state (EOSs) of nuclear matter. We analyze how variations in progenitor structure and the EOS of dense matter above nuclear saturation density affect the PNS evolution and subsequent BH formation. Comparisons of core collapse for a given progenitor star and different EOSs show that the path traced by the PNS in mass-specific entropy phase space ${M}_{\mathrm{grav}}^{\mathrm{PNS}}-\tilde{s}$ is well correlated with the progenitor compactness and is almost EOS independent, apart from the final end point. Furthermore, BH formation occurs, to a very good approximation, soon after the PNS overcomes the maximum gravitational mass supported by a hot NS with constant specific entropy equal to $\tilde{s}$. These results show a path to constraining the temperature dependence of the EOS through the detection of neutrinos from a failed galactic supernova.
Neutrino scattering and absorption rates of relevance to supernovae and neutron star mergers are obtained from nuclear matter dynamical structure functions that encode many-body effects from nuclear mean fields and correlations. In this study, we employ nuclear interactions from chiral effective field theory to calculate the density, spin, isospin, and spin-isospin response functions of warm β-equilibrium nuclear matter. We include corrections to the single-particle energies in the mean field approximation as well as vertex corrections resummed in the random phase approximation (RPA), including both direct and exchange diagrams. We find that correlations included through the RPA redistribute the strength of the response to higher energy for neutrino absorption and lower energy for antineutrino absorption. This tends to suppress the absorption rate of electron neutrinos across all relevant energy scales. In contrast, the inclusion of RPA correlations enhances the electron antineutrino absorption rate at low energy and suppresses the rate at high energy. These effects are especially important at high-density and in the vicinity of the neutrino decoupling region. Implications for heavy element nucleosynthesis, electromagnetic signatures of compact object mergers, supernova dynamics, and neutrino detection from galactic supernovae are discussed briefly.
Abstract We construct four equation of state (EoS) tables, tabulated over a range of temperatures, densities, and charge fractions, relevant for neutron star applications such as simulations of neutron star mergers. The EoS are computed from a relativistic mean-field theory constrained by the pure neutron matter EoS from chiral effective field theory, inferred properties of isospin-symmetric nuclear matter, and astrophysical observations of neutron star structure. To model nuclear matter at low densities, we attach an EoS that models inhomogeneous nuclear matter at arbitrary temperatures and charge fractions. The four EoS tables we develop are available from the CompOSE EoS repository https://compose.obspm.fr/eos/297 and gitlab.com/ahaber/qmc-rmf-tables .
A clear connection can be established between properties of nuclear matter and finite-nuclei observables, such as the correlation between the slope of the symmetry energy and the dipole polarizability, or between compressibility and the isoscalar monopole giant resonance excitation energy. Establishing a connection between realistic atomic nuclei and an idealized infinite nuclear matter leads to a better understanding of underlying physical mechanisms that govern nuclear dynamics. In this work, we aim to study the dependence of the binding energies and related quantities (e.g., location of drip lines, the total number of bound even-even nuclei) on the symmetry energy S 2 (ρ). The properties of finite nuclei are calculated by employing the relativistic Hartree-Bogoliubov model, assuming even-even axial and reflection symmetric nuclei. Calculations are performed by employing two families of relativistic energy density functionals, based on different effective Lagrangians, constrained to a specific symmetry energy at the saturation density J within the interval of 30–36 MeV. Nuclear binding energies and related quantities of bound nuclei are calculated between 8 ≤ Z ≤ 104 from the two-proton to the two-neutron drip line. As the neutron drip line is approached, the interactions with stiffer J tend to predict more bound nuclei, resulting in a systematic shift of the two-neutron drip line towards more neutron-rich nuclei. Consequentially, a correlation between the number of bound nuclei N nucl and S 2 (ρ) is established for a set of functionals constrained using the similar optimization procedures. Furthermore the direction of the relationship between the number of bound nuclei and the symmetry energy highly depends on the density under consideration.
In this work, using as references the posterior probability distribution functions of the equation of state (EOS) parameters inferred from the radii of canonical neutron stars (NSs) reported by the LIGO/VIRGO and NICER Collaborations based on their observations of GW170817 and PSR J0030+0451, we investigate how future radius measurements of more massive NSs will improve our current knowledge about the EOS of superdense neutron-rich nuclear matter, especially its symmetry energy term. Within the Bayesian statistical approach using an explicitly isospin-dependent parametric EOS for the core of NSs, we infer the EOS parameters of superdense neutron-rich nuclear matter from three sets of imagined mass–radius correlation data representing typical predictions by various nuclear many-body theories, that is, the radius stays the same, decreases, or increases with increasing NS mass within ±15% between 1.4 and 2.0 M$_\odot$. The corresponding NS average density increases quickly or slowly or slightly decreases as the NS mass increases from 1.4 to 2.0 M$_\odot$. While the EOSs of symmetric nuclear matter (SNM) inferred from the three data sets are approximately the same, the corresponding symmetry energies above about twice the saturation density of nuclear matter are very different, indicating that the radii of massive NSs carry important information about the high-density behavior of nuclear symmetry energy with little influence from the remaining uncertainties of the SNM EOS at suprasaturation densities.
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 model mixing (BMM) is a statistical technique that can combine constraints from different regions of an input space in a principled way. Here we extend our BMM framework for the equation of state (EOS) of strongly interacting matter from symmetric nuclear matter to asymmetric matter, specifically focusing on zero-temperature, charge-neutral, 𝛽-equilibrated matter. We use Gaussian processes (GPs) to infer constraints on the neutron-star matter EOS at intermediate densities from two different microscopic theories: chiral effective-field theory (𝜒EFT) at baryon densities around nuclear saturation, 𝑛 𝐵 ∼ 𝑛 0 , and perturbative QCD at asymptotically high baryon densities, 𝑛 𝐵 ⩾ 20𝑛 0 . The uncertainties of the 𝜒EFT and pQCD EOSs are obtained using the BUQEYE truncation error model. We demonstrate the flexibility of our framework through the use of two categories of GP kernels: conventional stationary kernels and a nonstationary changepoint kernel. We use the latter to explore potential constraints on the dense matter EOS by including exogenous data representing theory predictions and heavy-ion collision measurements at densities ⩾ 2𝑛 0 . We also use our EOSs to obtain neutron-star mass-radius relations and their uncertainties. Finally, our framework, whose implementation will be available through a GitHub repository, provides a prior distribution for the EOS that can be used in large-scale neutron-star inference frameworks.
The development of systematic effective field theories (EFTs) for nuclear forces and advances in solving the nuclear many-body problem have greatly improved our understanding of dense nuclear matter and the structure of finite nuclei. For global nuclear calculations, density functional theories (DFTs) have been developed to reduce the complexity and computational cost required in describing nuclear systems. However, DFT often makes approximations and assumptions about terms included in the functional, which may introduce systematic uncertainties compared to microscopic calculations using EFTs. In this work, we investigate possible avenues of improving nuclear DFT using nonlinear relativistic mean-field (RMF) theory. We explore the impact of RMF model extensions by fitting the nonlinear RMF model to predictions of nuclear matter and selected closed-shell nuclei using four successful chiral EFT Hamiltonians. We find that these model extensions are impactful and important in capturing the physics present within chiral Hamiltonians, particularly for charge radii and neutron skins of closed-shell nuclei. However, there are additional effects that are not captured within the RMF model, particularly within the isoscalar sector of RMF theory. Additional model extensions and the reliability of the nonlinear RMF model are discussed.
In this work, we construct the equation of state for high-density neutron star matter at zero temperature using the two-flavor Nambu–Jona-Lasinio (NJL) model as an effective theory of QCD. We build nuclear matter, quark matter, and the mixed phases from the same NJL Lagrangian, which has been used to model free and in-medium hadrons as well as nuclear systems. A focus here is to determine if the same coupling constants in the scalar diquark and vector meson channels, which give a good description of nucleon structure and nuclear matter, can also be used for the color superconducting high-density quark matter phase. We find that this is possible for the scalar diquark (pairing) interaction, but the vector meson interaction has to be reduced so that superconducting quark matter becomes the stable phase at high densities. We compare our equation of state with recent phenomenological parametrizations based on generic stability conditions for neutron stars. We find that the maximum mass of a hybrid star, with a color superconducting quark matter core, exceeds 2.01 ± 0.04 M⊙, which is the value of the recently observed massive neutron star PSR J0348+0432. The mass-radius relation is also consistent with gravitational wave observations (GW170817).