Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nuclear interactions”

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 73 records · Page 4

Nuclear-matter saturation and symmetry energy within Δ -full chiral effective field theory

Nuclear saturation and the symmetry energy are key properties of low-energy nuclear physics that depend on fine details of the nuclear interaction. The equation of state around saturation is also an important anchor for extrapolations to higher densities and studies of neutron stars. Here we develop a unified statistical framework that uses realistic nuclear forces to link the theoretical modeling of finite nuclei and infinite nuclear matter. We construct fast and accurate emulators for nuclear-matter observables and employ an iterative history-matching approach to explore and reduce the enormous parameter domain of Δ -full chiral interactions. We perform rigorous uncertainty quantification and find that model calibration including O 16 observables gives saturation predictions that are more precise than those that only use few-body data. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Chiral effective field theory calculation of neutrino reactions in warm neutron-rich matter

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.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Quantifying uncertainties in α -nucleus reaction dynamics informed from first principles

The ab initio symmetry-adapted no-core shell model is a microscopic many-body method which naturally describes challenging collective and clustering features of atomic nuclei. Wave functions and observables computed with realistic nucleon-nucleon forces in this framework are tied to first principles, and are hence well-suited for rigorous uncertainty quantification. We discuss α-deuteron and α- 12 C cluster potentials informed by symmetry-adapted calculations, and propagate uncertainties in the effective binary cluster method as well as those in the nuclear interaction to reaction observables, namely scattering phase shifts, cross sections, partial widths and resonance energies. Here, we find that the overall uncertainties are dominated by those originating in the underlying nuclear force, speaking to the need for tighter constraints on realistic nucleon-nucleon interactions.

Ab initio↗

Windowed multipole representation of R -matrix cross sections

Nuclear cross sections are basic inputs to any nuclear computation. Campaigns of experiments are fitted with the parametric R-matrix model of quantum nuclear interactions, and the resulting cross sections are documented—both pointwise and as resonance parameters (with uncertainties)—in standard evaluated nuclear data libraries (ENDF, JEFF, BROND, JENDL, CENDL, TENDL): these constitute our common knowledge of fundamental low-energy nuclear cross sections. In the past decade, a collaborative effort has been deployed to establish a new nuclear cross-section library format—the Windowed Multipole Library—with the goal of considerably reducing the computational cost of cross-section calculations in nuclear transport simulations. This work lays the theoretical foundations underpinning these efforts. From general R-matrix scattering theory, we derive the windowed multipole representation of nuclear cross sections. Though physically and mathematically equivalent to R-matrix cross sections, the windowed multipole representation is particularly well suited for subsequent temperature treatment of angle-integrated cross sections, in particular Doppler broadening, which is the averaging of cross sections over the thermal motion of the target atoms. Doppler broadening is of critical importance in neutron transport applications, as it ensures the stability of many nuclear reactors (negative thermal reactivity). Yet, Doppler broadening of nuclear cross sections has been a considerable bottleneck for nuclear transport computations, often requiring memory-costly pretabulations. We show that the windowed multipole representation can perform accurate Doppler broadening analytically (up to the first reaction threshold), from which we derive cross-section temperature derivatives to any order—all computable on the fly (without precalculations stored in memory). Furthermore, we here establish a way of converting the R-matrix resonance parameters uncertainty (covariance matrices) into windowed multipole parameters uncertainty. We show that generating stochastic nuclear cross sections by sampling from the resulting windowed multipole covariance matrix can reproduce the cross-section uncertainty in the original nuclear data file. The windowed multipole representation is therefore a novel nuclear physics formalism able to generate Doppler broadened stochastic nuclear cross sections on the fly, unlocking breakthrough computational gains for nuclear computations. Through this foundational paper, we hope to make the windowed multipole representation accessible, reproducible, and usable for the nuclear physics community, as well as provide the theoretical basis for future research on expanding its capabilities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NUCLEI (SciDAC-4) effort at The University of Tennessee, Knoxville

This Report summarizes the accomplishments of NUCLEI researchers at The University of Tennessee, Knoxville (UTK) and those supported via subcontracts at the University of North Carolina at Chapel Hill (UNC) and Iowa State University (ISU). The SciDAC-4 NUCLEI collaboration had the goal to advance the computation of atomic nuclei and to make reliable predictions including uncertainty estimates for observables of interest. These goals were achieved by developing improved nuclear interactions and currents that are input to the computations, and by advancing the algorithms and techniques used for computing nuclear observables of interest. The computed results advanced our understanding of how nuclei work; they helped to guide and interpret experiments at existing low-energy nuclear physics facilities, the Facility for Rare Isotope Beams, Jefferson Laboratory, and neutrino experiments. In addition, predictions were made that await confrontation with experimental results.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Progress in Constraining Nuclear Symmetry Energy Using Neutron Star Observables Since GW170817

The density dependence of nuclear symmetry energy is among the most uncertain parts of the Equation of State (EOS) of dense neutron-rich nuclear matter. It is currently poorly known especially at suprasaturation densities partially because of our poor knowledge about isovector nuclear interactions at short distances. Because of its broad impacts on many interesting issues, pinning down the density dependence of nuclear symmetry energy has been a longstanding and shared goal of both astrophysics and nuclear physics. New observational data of neutron stars including their masses, radii, and tidal deformations since GW170817 have helped improve our knowledge about nuclear symmetry energy, especially at high densities. Based on various model analyses of these new data by many people in the nuclear astrophysics community, while our brief review might be incomplete and biased unintentionally, we learned in particular the following: (1) The slope parameter L of nuclear symmetry energy at saturation density ρ0 of nuclear matter from 24 new analyses of neutron star observables was about L≈57.7±19 MeV at a 68% confidence level, consistent with its fiducial value from surveys of over 50 earlier analyses of both terrestrial and astrophysical data within error bars. (2) The curvature Ksym of nuclear symmetry energy at ρ0 from 16 new analyses of neutron star observables was about Ksym≈−107±88 MeV at a 68% confidence level, in very good agreement with the systematics of earlier analyses. (3) The magnitude of nuclear symmetry energy at 2ρ0, i.e., Esym(2ρ0)≈51±13 MeV at a 68% confidence level, was extracted from nine new analyses of neutron star observables, consistent with the results from earlier analyses of heavy-ion reactions and the latest predictions of the state-of-the-art nuclear many-body theories. (4) While the available data from canonical neutron stars did not provide tight constraints on nuclear symmetry energy at densities above about 2ρ0, the lower radius boundary R2.01=12.2 km from NICER’s very recent observation of PSR J0740+6620 of mass 2.08±0.07M⊙ and radius R=12.2–16.3 km at a 68% confidence level set a tight lower limit for nuclear symmetry energy at densities above 2ρ0. (5) Bayesian inferences of nuclear symmetry energy using models encapsulating a first-order hadron–quark phase transition from observables of canonical neutron stars indicated that the phase transition shifted appreciably both L and Ksym to higher values, but with larger uncertainties compared to analyses assuming no such phase transition. (6) The high-density behavior of nuclear symmetry energy significantly affected the minimum frequency necessary to rotationally support GW190814’s secondary component of mass (2.50–2.67) M⊙ as the fastest and most massive pulsar discovered so far. Overall, thanks to the hard work of many people in the astrophysics and nuclear physics community, new data of neutron star observations since the discovery of GW170817 have significantly enriched our knowledge about the symmetry energy of dense neutron-rich nuclear matter.

79 ASTRONOMY AND ASTROPHYSICS↗

The Unitarity-Limit Expansion for Two Nucleons with Perturbative Pions: Digest and Ideas

Theorists love nontrivial fixed points. In the Unitarity Limit, the NN 𝑆-wave binding energies are zero, the scattering lengths infinite, Physics is universal, i.e. insensitive to details of the interactions, and observables display richer symmetries, namely invariance under both scaling and Wigner’s combined SU(4) transformation of spin and isospin. In “Pionless” EFT, both are explicitly but weakly broken and hence perturbative in the Unitarity Window (phase shifts 45° ≲ δ(k) ≲ 135°, i.e. momenta k≈mπ). This Unitarity Expansion provides strong hints that Nuclear Physics resides indeed in a sweet spot: bound weakly enough to be insensitive to the details of the nuclear interaction; and therefore interacting strongly enough that the NN scattering lengths are perturbatively close to the Unitarity Limit. In this paradigm change, NN details are less important than NNN interactions to explain the complexity and patterns of the nuclear chart. This presentation is a digest of the first quantitative exploration of corrections to this picture when pions are included [1] (see there for a more comprehensive list of references). Since the pion mass and decay constant introduce dimensionful scales in the NN system, they explicitly break the symmetries of the Unitarity fixed point. In χEFT, these symmetries must therefore be hidden and instead be classified as emergent.

Griesshammer, Harald W. [The George Washington Uni↗

R -mode Stability of GW190814's Secondary Component as a Supermassive and Superfast Pulsar

The nature of GW190814's secondary component m 2 of mass (2.50–2.67)M ⊙ in the mass gap between the currently known maximum mass of neutron stars and the minimum mass of black holes is currently under hot debate. Among the many possibilities proposed in the literature, the m 2 was suggested as a superfast pulsar while its r-mode stability against the run-away gravitational radiation through the Chandrasekhar-Friedman-Schutz mechanism is still unknown. Previously, Fortin et al. constructed a sample of 33 unified equations of state using the same nuclear interactions from the crust to the core consistently; from that sample we use those equations that fulfill all currently known astrophysical and nuclear physics constraints to compare the minimum frequency required for m 2 to rotationally sustain a mass greater than 2.50 M ⊙ with the critical frequency above which the r-mode instability occurs. Here, we use two extreme damping models assuming that the crust is either perfectly rigid or elastic. Using the stability of 19 observed low-mass X-ray binaries as an indication that the rigid crust damping of the r-mode dominates within the models studied, we find that m 2 is r-mode-stable while rotating with a frequency higher than 870.2 Hz (0.744 times its Kepler frequency of 1169.6 Hz) as long as its temperature is lower than about 3.9 × 10 7 K, further supporting the proposal that GW190814's secondary component is a supermassive and superfast pulsar.

79 ASTRONOMY AND ASTROPHYSICS↗

Converged ab initio calculations of heavy nuclei

We propose a novel storage scheme for three-nucleon (3N) interaction matrix elements relevant for the normal-ordered two-body approximation used extensively in ab initio calculations of atomic nuclei. This scheme reduces the required memory by approximately two orders of magnitude, which allows the generation of 3N interaction matrix elements with the standard truncation of E 3max =28, well beyond the previous limit of 18. We demonstrate that this is sufficient to obtain the ground-state energy of 132 Sn converged to within a few MeV with respect to the E 3max truncation. In addition, we study the asymptotic convergence behavior and perform extrapolations to the un-truncated limit. Finally, we investigate the impact of truncations made when evolving free-space 3N interactions with the similarity renormalization group. We find that the contribution of blocks with angular momentum J rel > 9/2 to the ground-state energy is dominated by a basis-truncation artifact, which vanishes in the large-space limit, so these computationally expensive components can be neglected. For the two sets of nuclear interactions employed in this work, the resulting binding energy of 132 Sn agrees with the experimental value within theoretical uncertainties. This work enables converged ab initio calculations of heavy nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Electrical conductivity of a warm neutron star crust in magnetic fields: Neutron-drip regime

We compute the anisotropic electrical conductivity tensor of the inner crust of a compact star at nonzero temperature by extending a previous work on the conductivity of the outer crust. The physical scenarios, where such crust is formed, involve protoneutron stars born in supernova explosions, binary neutron star mergers, and accreting neutron stars. The temperature-density range studied covers the transition from a semidegenerate to a highly degenerate electron gas and assumes that the nuclei form a liquid, i.e., the temperature is above the melting temperature of the lattice of nuclei. The electronic transition probabilities include (i) the screening of electron-ion interaction in the hard-thermal-loop approximation for the QED plasma, (ii) the correlations of the ionic component in a one-component plasma, and (iii) finite nuclear size effects. The conductivity tensor is obtained from the Boltzmann kinetic equation in relaxation time approximation accounting for the anisotropy introduced by a magnetic field. The sensitivity of the results towards the matter composition of the inner crust is explored by using several compositions of the inner crust, which were obtained using different nuclear interactions and methods of solving the many-body problem. The standard deviations of relaxation time and components of the conductivity tensor from the average are below ≤25% except close to crust-core transition, where nonspherical nuclear structures are expected. Finally, our results can be used in dissipative magnetohydrodynamics simulations of warm compact stars.

Physics↗

Partial Wave Analysis of Strange Mesons Decaying to K+p+p- in the Reaction ¿p ¿ K+p+p-¿(1520) and the Commissioning of the GlueX DIRC Detector

Hadron spectroscopy is a cornerstone of our understanding of the strong nuclear interactions. Studying the hadron spectrum led to the postulation of quarks and gluons, and the development of Quantum Chromodynamics (QCD), the theory of the strong nuclear force. Today hadron spectroscopy provides an important test of QCD, particularly in the non-perturbative energy regime. One such test is the existence of hybrid hadrons that have gluonic degrees of freedom, e.g. q ?qg states, that are allowed by QCD but have remained elusive in experimental searches. The GlueX experiment located at Thomas Jefferson National Accelerator Facility, is designed to map the light meson spectrum, including hybrid mesons, produced by a linearly polarized photon beam that is scattered off of a liquid hydrogen target. This dissertation concerns the installation and commissioning of the GlueX DIRC upgrade and the analysis of three reactions measured by the GlueX experiment. The GlueX DIRC is a particle identification detector installed as an upgrade to the GlueX experiment in 2019 to improve the experiment?s ability to distinguish pions and kaons. The DIRC uses Cherenkov photons radiated by charged particles that travel faster than the speed of light in the detector medium to determine the velocity of said particles. The velocity measurement is used in conjunction with momentum measurements from other detectors to determine the mass of the particle. The DIRC will be essential to the future of strange meson spectroscopy in GlueX. The reactions ?p ? K+K??+??p, ?p ? K+K??0?0p, ?p ? K+K??0?p were initially chosen for this analysis to search for the hybrid meson candidate ?(2170) and related resonances. The ?(2170) has been observed in e+e? collider experiments decaying to K+K??+?? and K+K??0?0, among other final states. The ?(2170) was not observed in this analysis. The reaction ?p ? K+?+???(1520), a subset of the ?p ? K+K??+??p reaction, was analyzed with a Partial Wave Analysis (PWA) in order to study strange mesons that decay to K+?+??. In the PWA, a useful model for the angular distributions was found but no significant structures could be observed due to statistical limitations.

Hurley, Andrew↗

Lattice QCD Inputs for nuclear double beta decay

Second order β -decay processes with and without neutrinos in the final state are key probes of nuclear physics and of the nature of neutrinos. Neutrinoful double- β decay is the rarest Standard Model process that has been observed and provides a unique test of the understanding of weak nuclear interactions. Observation of neutrinoless double- β decay would reveal that neutrinos are Majorana fermions and that lepton number conservation is violated in nature. While significant progress has been made in phenomenological approaches to understanding these processes, establishing a connection between these processes and the physics of the Standard Model and beyond is a critical task as it will provide input into the design and interpretation of future experiments. The strong-interaction contributions to double- β decay processes are non-perturbative and can only be addressed systematically through a combination of lattice Quantum Chromoodynamics (LQCD) and nuclear many-body calculations. In this review, current efforts to establish the LQCD connection are discussed for both neutrinoful and neutrinoless double- β decay. LQCD calculations of the hadronic contributions to the neutrinoful process $nn → ppe^-e^-\bar{v}_e\bar{v}_e$ and to various neutrinoless pionic transitions are reviewed, and the connections of these calculations to the phenomenology of double- β decay through the use of effective field theory (EFTs) is highlighted. At present, LQCD calculations are limited to small nuclear systems, and to pionic subsystems, and require matching to appropriate EFTs to have direct phenomenological impact. However, these calculations have already revealed qualitatively that there are terms in the EFTs that can only be constrained from double- β decay processes themselves or using inputs from LQCD. Finally, future prospects for direct calculations in larger nuclei are also discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Expected performance of the FOOT experiment

The main objective of the FOOT (FragmentatiOn of Target) experimentis to measure the double differential cross-sections with respect to the kineticenergy and the angle of emission of fragments produced in nuclear interactions withan accuracy of 5%. Measurements will use the inverse kinematics approach to studythe interaction of heavy ions with proton-rich targets and will be performed in twoenergy regions, one up to 400 MeV/u, to assess the effects of nuclear fragmentationin hadrontherapy treatment with proton beams, and up to 700 MeV/u, to improvethe knowledge of the fragmentation of ions on nuclei to help with the design of spacevehicles shielding.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Exploring the Strong Interaction of Three-Body Systems at the LHC

Deuterons are atomic nuclei composed of a neutron and a proton held together by the strong interaction. Unbound ensembles composed of a deuteron and a third nucleon have been investigated in the past using scattering experiments, and they constitute a fundamental reference in nuclear physics to constrain nuclear interactions and the properties of nuclei. In this work, K + -d and p-d femtoscopic correlations measured by the ALICE Collaboration in proton-proton (pp) collisions at $\sqrt{s}$ = 13 TeV at the Large Hadron Collider (LHC) are presented. It is demonstrated that correlations in momentum space between deuterons and kaons or protons allow us to study three-hadron systems at distances comparable with the proton radius. The analysis of the K + -d correlation shows that the relative distances at which deuterons and protons or kaons are produced are around 2 fm. The analysis of the p-d correlation shows that only a full three-body calculation that accounts for the internal structure of the deuteron can explain the data. In particular, the sensitivity of the observable to the short-range part of the interaction is demonstrated. These results indicate that correlations involving light nuclei in pp collisions at the LHC will also provide access to any three-body system in the strange and charm sectors.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Multiphysics FISPACT-II and TENDL-2019 simulation: neutron-induced damage metrics

Nuclear interactions can be the source of atomic displacement, irradiation-induced defects and transmutation in structural materials. Such quantities are derived from, or can be correlated to, nuclear kinematic simulations of primary atomic energy distributions spectra and the quantification of the numbers of secondary defects produced per primary as a function of the available recoils, residual and emitted, energies. Recoil kinematics of neutral, residual, gas, proton, alpha particle emissions are now more rigorously treated based on recent, complete and enhanced nuclear data parsed in state-of-the-art processing tools. Defect production metrics are the starting point in the complex problem of correlating and simulating the behaviour of materials under irradiation, as experimental information is rare or scattered. Detailed, segregated primary knock-on-atom metrics are now becoming available as the starting point of further simulation processes of isolated and clustered defects in material lattices. This allows more materials, lattices, neutron incident energy ranges, and irradiation conditions to be explored with sufficient data to adequately cover both standard and novel applications and materials: the broader reactor applications landscape. The damage metrics of of materials are systematically explored under typical but different reactor's type environment. The inventory code FISPACT-II combined with the enhanced nuclear data forms of the TENDL-2019 libraries allow one to not only calculate dpa from mostly scattering events but to also properly predict gas production, nuclear heating and transmutation under the same conditions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

From Spin to Pseudospin Symmetry: The Origin of Magic Numbers in Nuclear Structure

Magic numbers lie at the heart of nuclear structure, reflecting enhanced stability in nuclei with closed shells. While the emergence of magic numbers beyond 20 is commonly attributed to strong spin-orbit coupling, the microscopic origin of the spin-orbit potential remains elusive, owing to its dependence on the resolution scale and renormalization scheme of nuclear forces. Here, we investigate the evolution of nuclear shell structure with varying momentum resolution in nuclear interactions derived from chiral effective field theory, using the similarity renormalization group to link different scales. We uncover a novel transition from spin symmetry to pseudospin symmetry as the resolution scale decreases, during which magic numbers emerge naturally. A similar pattern is found in calculations using relativistic one-boson-exchange potentials, underscoring the robustness of the phenomenon. This establishes a direct connection between realistic nuclear forces with a high resolution scale and effective nuclear forces at coarse-grained scales, offering a first-principles explanation for the origin of magic numbers and pseudospin symmetry in nuclear shell structure and new insights into the structure of exotic nuclei far from stability

Energy levels↗

Computational Modeling and Simulation for Nonproliferation: The History (and Present) of Monte Carlo and the MCNP(R) Code at Los Alamos [Slides]

The emergence of the Monte Carlo method as a research tool springs from work done at Los Alamos in the 1940s.Monte Carlo and the MCNP Code have been and continue to be developed at Los Alamos for many decades. From basic science in support of understanding nuclear interaction physics to global security applications, application uses of the code are extensive. Recent R&D projects and code modernization efforts make the MCNP code a great tool for nuclear nonproliferation applications. In collaboration with nuclear safeguards experts, new training has just recently been developed to help new practitioners learn how to use the code for nuclear safeguards applications.

97 MATHEMATICS AND COMPUTING↗

Optical potential uncertainties on exclusive tens-of-MeV Ar cross sections

Large liquid argon time projection chamber (LArTPC) neutrino detectors, such as those planned for the Deep Underground Neutrino Experiment (DUNE), show considerable promise as a platform for next-generation measurements of supernova neutrinos. Thanks to the neutron excess in $^{40}$Ar as well as the detailed tracking possible with LArTPCs, these detectors are expected to be uniquely capable of measuring supernova electron neutrinos with high statistics and minimal backgrounds. However, these technological advantages come at the price of complexities in data interpretation; reconstruction of the incident energies of supernova neutrinos in a future LArTPC-based analysis will be subject to a variety of systematic uncertainties related to nuclear interaction modeling. In this talk, we present a study of a subset of these uncertainties as implemented in the MARLEY event generator used by DUNE and other LArTPC neutrino experiments. Variations to the optical potential used in the MARLEY nuclear de-excitation model are applied to the calculation of exclusive tens-of-MeV neutrino-nucleus cross sections for the first time, and the impact on observables of interest for supernova neutrino detection is examined. We also discuss prospects for quantifying other sources of uncertainty on the MARLEY neutrino interaction model.

El-Haj, Luca Abu↗