Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nuclear matrix elements”

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 343 records · Page 19

Multistep Coulomb excitation of Ni 64 : Shape coexistence and nature of low-spin excitations

Here, the structure of 64 Ni, the heaviest stable Ni isotope, has been investigated via high-statistics, multistep safe Coulomb excitation to search for shape coexistence, a phenomenon recently observed in neutron-rich 66 Ni and 70 Ni as well as in doubly magic, N = 40, 68 Ni. The study was motivated by recent, state-of-the-art Monte Carlo shell-model calculations (MCSM), where a Hamiltonian with effective interactions incorporating the monopole tensor force predicts the existence of shape coexistence, also in the lower-mass 62,64 Ni isotopes. A set of transition and static E2 matrix elements for both yrast and near-yrast structures was extracted from the differential Coulomb excitation cross sections. From comparisons between the new results and MCSM as well as other shell-model calculations, a clearer picture of the structure of 64 Ni emerges. Specifically, the low-spin states are shown to be dominated by proton and neutron excitations mainly within the fp shell, with minimal contribution from the g 9/2 shape-driving neutron orbital. The agreement between experimental data and MCSM results indicates a small oblate deformation for the $0$$^{+}_{2}$ level and a spherical shape for the $0$$^{+}_{3}$ state. In addition, the small upper limit determined for the B(E2) probability of a transition associated with the decay of the recently observed 3463-keV, $0$$^{+}_{4}$ state agrees with its proposed assignment to a prolate shape, herewith providing first evidence for triple shape coexistence in a stable Ni isotope.

59 ≤ A ≤ 89↗

Evaluation of Sb-Nd and Te-Nd phases within the U-Zr fuel matrix and their interactions with HT9 alloy

Antimony (Sb) and tellurium (Te) were investigated as potential additives for U-10Zr (wt.%) metallic fuel to limit the fuel-cladding chemical interaction (FCCI) with HT-9 alloy. Neodymium (Nd) was utilized to simulate the formation of lanthanide-based solid fission products which are known to play a detrimental role in FCCI. Fuel alloys of U-Zr-Sb-Nd and U-Zr-Te-Nd were evaluated in their annealed condition and compared against their as-cast conditions. Isothermal diffusion couple experiments were performed between U-Zr-Nd, U-Zr-Sb-Nd, and U-Zr-Te-Nd against the cladding alloy HT9 to evaluate the effectiveness of the additives to stabilize Nd within the fuel alloys, as well as investigate the interaction regions that form between the different fuel alloys and HT9. Further, SbNd and Sb 3 Nd 4 , and TeNd are found to be the primary neodymium-based phases formed in the U-Zr-Sb-Nd and U-Zr-Te-Nd alloys, respectively. The zirconium-based phase, Zr 2 Sb, is also found to form within the former alloy. All phases were found to remain stable through the diffusion experiments and exhibited no interaction with HT9 constituent elements. Preferential interaction between Nd with additivities Te and Sb compared to constituting elements in HT9 was further verified based on density functional theory (DFT) calculated enthalpy of mixing.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

NEAMS Technical Area Support in MOOSE

The Multiphysics Object-Oriented Simulation Environment (MOOSE) framework is a foundational capability used by the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program to create over 15 different simulation tools for advanced nuclear reactors. Due to this ubiquity, improvements to the framework in support of modeling and simulation goals are critical to the program. These improvements can take many forms, including optimization, improved user experience, streamlined application programming interfaces (APIs), parallelism, and other new capabilities. The work described in this report was conducted in direct support of the simulation tools and has already been deployed. The capabilities outlined in this report include implementing hash table matrix assembly for efficient sparsity pattern construction for contact in BISON, developing re-step testing infrastructure for ensuring the viability of overlapping domain coupling between SAM and Pronghorn, allowing unique preconditioners for single-input multi-system solves, supporting multi-system in MOOSE’s workhorse executioners, and many more smaller feature enhancements and bug fixes.

97 - MATHEMATICS AND COMPUTING↗

Direct quarkonium production in DIS from a joint CGC and NRQCD framework

We compute the differential cross section for direct quarkonium production in high-energy electron-nucleus collisions at small 𝑥. Our computation is performed within the nonrelativistic QCD factorization formalism that separates the calculation into short distance coefficients and long distance matrix elements that depend on the color and spin of the state. We obtain the short distance coefficients of the production of the heavy quark pair within the framework of the color glass condensate effective field theory, which resums coherent multiple interactions of the heavy quark pair with the nucleus to all orders. Our results are expressed as the convolution of perturbatively calculable functions with multipoint lightlike Wilson line correlators. In the correlation limit, we establish the correspondence between our color glass condensate formulation with calculations employing the transverse momentum dependent (TMD) framework. We extend this correspondence by resumming kinematic power corrections within the improved TMD framework, which interpolates between the TMD formalism and 𝑘 ⊥ -factorization formalism. We present a detailed numerical analysis, focusing on 𝐽/𝜓 production in the kinematics accessible at the future Electron-Ion Collider, highlighting the importance of genuine higher-order saturation contributions when the electron collides with a large nucleus. Our results are also valid in the photoproduction limit where we expect the largest contribution from genuine higher-order saturation contributions which could be accessed in ultraperipheral collisions of relativistic heavy ions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Lattice QCD calculation of 𝑥-dependent meson distribution amplitudes at physical pion mass with threshold logarithm resummation

We present a lattice quantum chromodynamics (QCD) calculation of the 𝑥-dependent pion and kaon distribution amplitudes (DA) in the framework of large momentum effective theory. This calculation is performed on a fine lattice of 𝑎 = 0.076 fm at physical pion mass, with the pion boosted to 1.8 GeV and kaon boosted to 2.3 GeV. We renormalize the matrix elements in the hybrid scheme and match to Math output error with a subtraction of the leading renormalon in the Wilson-line mass. The perturbative matching is improved by resumming the large logarithms related to the small quark and gluon momenta in the soft-gluon limit. After resummation, we demonstrate that we are able to calculate a range of 𝑥 ∈[𝑥 0 ,1 − 𝑥 0 ] with 𝑥 0 = 0.25 for pion and 𝑥 0 = 0.2 for kaon with theoretical systematic errors under control. The kaon DA is shown to be slighted skewed, and narrower than pion DA. Although the 𝑥-dependence cannot be direct calculated beyond these ranges, we estimate higher moments of the pion and kaon DAs by complementing our calculation with short-distance factorization.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Tensor network simulations of quasi-GPDs in the massive Schwinger model

Generalized parton distribution functions (GPDs) are off-diagonal light-cone matrix elements that encode the internal structure of hadrons in terms of quark and gluon degrees of freedom. In this work, we present the first nonperturbative study of quasi-GPDs in the massive Schwinger model, quantum electrodynamics in 1+1 dimensions (QED 2 ), within the Hamiltonian formulation of lattice field theory. Quasidistributions are spatial correlation functions of boosted states, which approach the relevant light-cone distributions in the luminal limit. Using tensor networks, we prepare the first excited state in the strongly coupled regime and boost it to close to the light-cone on lattices of up to 400 lattice sites. We compute both quasiparton distribution functions and, for the first time, quasi-GPDs, and study their convergence for increasingly boosted states. In addition, we perform analytic calculations of GPDs in the two-particle Fock-space approximation and in the Reggeized limit, providing qualitative benchmarks for the tensor network results. Our analysis establishes computational benchmarks for accessing partonic observables in low-dimensional gauge theories, offering a starting point for future extensions to higher dimensions, non-Abelian theories, and quantum simulations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Proton and neutron contributions to the quadrupole transition strengths in 39 Ca and 39 K studied by lifetime measurements of mirror transitions

The E2 transition matrix elements of isobaric multiplets are expected to follow a linear trend as a function of isospin projection. However, measurements of the 2 + → 0 + transitions in the A = 38 triplet of Ca, K, and Ar show a deviation from this trend with an enhanced transition strength in 38 Ca with respect to its mirror 38 Ar. We have studied analogue 11/2 – → 7/2 – E2 transitions in 39 Ca and its mirror partner 39 K to determine if this enhancement persists in neighboring Ca isotopes. Recoil-distance lifetime measurements of 39 Ca and 39 K were performed utilizing a 42 Sc secondary beam, the TRIPLEX plunger, the GRETINA array, and the S800 spectrograph. Our data provide a lifetime measurement of the (11/2 – ) state in 39 Ca as well as an improved lifetime result for the (9/2 – ) state, while the 39 K data are used to validate the present analysis. Furthermore, a comparison of the present data to shell-model calculations suggests an enhanced transition strength in 39 Ca, pointing to both proton and neutron contributions to core excitations across the Z = N = 20 shell gaps in close proximity to 40 Ca.

39 ≤ A ≤ 58↗

Low-Activity Waste Glass Standards Preparation and Characterization for Calibration of Analytical Instruments

A matrix of seven low-activity waste (LAW) glasses was fabricated and characterized to be used by the Hanford Waste Treatment and Immobilization Plant (WTP) as analytical instrument calibration and matrix interference standards when measuring LAW glass compositions. The matrix of glasses was designed to represent the range of LAW glass compositions that will be generated by the WTP. Samples from each of the seven glasses were measured with electron probe microanalysis (EPMA), ion chromatography (IC), and bulk inductively coupled plasma mass spectroscopy (ICP-MS) to determine the variability and accuracy across batches and within individually poured glass bar samples. High relative percent differences and standard deviations were calculated for components with low concentrations due to instrument detection limits. IC and ICP-MS percent differences and standard deviations were higher at low component concentrations. Otherwise, variability in measured compositions from IC and ICP-MS was generally lower than measurements from EPMA, with the exclusion of Cl, F, Si, and B. Overall, the two methods proved comparable and complementary. From EPMA, IC, and ICP-MS data, volatile elements were underrepresented compared to batched compositions, which suggested loss of these elements during melting. EPMA data was statistically analyzed to evaluate homogeneity within single bars and whole compositions for each of the seven glasses. Most of the variability within glass compositions occurred at the bar-to-bar level; however, no clear systematic trends were observed. Based on these analyses, when using these glasses as reference material, it is recommended that the overall mean and variance of each composition be used for performing instrument calibrations and analytical corrections.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Moments of axial-vector GPD from lattice QCD: quark helicity, orbital angular momentum, and spin-orbit correlation

In this work, we present a lattice QCD calculation of the Mellin moments of the twist-2 axial-vector generalized parton distribution (GPD), $\overset{\sim }{H}\left(x,\xi, t\right)$ , at zero skewness, ξ, with multiple values of the momentum transfer, t. Our analysis employs the short-distance factorization framework on ratio-scheme renormalized quasi-GPD matrix elements. The calculations are based on an N f = 2 + 1 + 1 twisted mass fermions ensemble with clover improvement, a lattice spacing of a = 0.093 fm, and a pion mass of m π = 260 MeV. We consider both the iso-vector and iso-scalar cases, utilizing next-to-leading-order perturbative matching while omitting the disconnected contributions and gluon mixing in the iso-scalar case. For the first time, we determine the Mellin moments of $\overset{\sim }{H}$ up to the fifth order. From these moments, we discuss the quark helicity and orbital angular momentum contributions to the nucleon spin, as well as the spin-orbit correlations of the quarks. Additionally, we perform a Fourier transform over the momentum transfer, which allows us to explore the spin structure in the impact-parameter space.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Measurement of the branching fraction of the ${{B} ^0} {\rightarrow }{{D} ^+_{s}} {{\pi } ^-} $ decay

A branching fraction measurement of the \({{B} ^0} {\rightarrow }{{D} ^+_{s}} {{\pi } ^-} \) decay is presented using proton–proton collision data collected with the LHCb experiment, corresponding to an integrated luminosity of \(5.0\,\text {fb} ^{-1} \) . The branching fraction is found to be \({\mathcal {B}} ({{B} ^0} {\rightarrow }{{D} ^+_{s}} {{\pi } ^-} ) =(19.4 \pm \) \(1.8\pm 1.3 \pm 1.2)\times 10^{-6}\) , where the first uncertainty is statistical, the second systematic and the third is due to the uncertainty on the \({{B} ^0} {\rightarrow }{{D} ^-} {{\pi } ^+} \) , \({{D} ^+_{s}} {\rightarrow }{{K} ^+} {{K} ^-} {{\pi } ^+} \) and \({{D} ^-} {\rightarrow }{{K} ^+} {{\pi } ^-} {{\pi } ^-} \) branching fractions. This is the most precise single measurement of this quantity to date. As this decay proceeds through a single amplitude involving a \(b{\rightarrow }u\) charged-current transition, the result provides information on non-factorisable strong interaction effects and the magnitude of the Cabibbo–Kobayashi–Maskawa matrix element \(V_{ub}\) . Additionally, the collision energy dependence of the hadronisation-fraction ratio \(f_s/f_d\) is measured through \({{\overline{B}} {}^0_{s}} {\rightarrow }{{D} ^+_{s}} {{\pi } ^-} \) and \({{B} ^0} {\rightarrow }{{D} ^-} {{\pi } ^+} \) decays.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Enabling the Broader Use of MOOSE for Nuclear Energy and Other Simulation

This Final Scientific and Technical Report summarizes work performed under the Phase IIA SBIR project “Enabling the Broader Use of MOOSE for Nuclear Energy and Other Simulation” (DE-SC0020906) from August 2023 through August 2025. The objective of the Phase IIA effort was to mature and harden capabilities developed during Phase II, with the goal of enabling practical interoperability between Coreform’s isogeometric analysis (IGA) technologies and the Multiphysics Object-Oriented Simulation Environment (MOOSE), while improving robustness, performance, and scalability for complex, nuclear-relevant geometries. Over the course of Phase IIA, the project established and validated an extraction-based interoperability pathway between Coreform tools and MOOSE. A combined mesh and matrix format was defined collaboratively with MOOSE developers and integrated into the solver, enabling standard MOOSE workflows to operate on data exported from Coreform’s IGA and Flex Representation Method (FRM) pipelines. Early demonstrations validated architectural compatibility using linear solid mechanics problems, while later efforts focused on benchmark testing and external use. By the end of the project period, engineers at BWXT were able to independently set up and execute a simulation using the Coreform–MOOSE workflow and provide direct feedback that informed further refinement. In parallel, substantial effort was devoted to improving the robustness of trimmed U-spline construction for complex CAD geometries. A growing test suite of nuclear-relevant models was compiled through collaboration with multiple stakeholders and used to drive extensive bug fixing and reliability improvements. These efforts resulted in improved robustness and performance, including the addition of fallback capabilities that enhance reliability when the underlying commercial CAD kernel fails. Performance-oriented work progressed later in the project, with the development and demonstration of methods to decompose complex geometries into structured subregions and updated data representations to support more efficient solver processing. Additionally, extensive enhancements to threadsafe parallel data structures and trimming operations established a foundation for scalable processing of large assemblies. Collaboration with Sandia National Laboratories on the SGM geometric modeling kernel advanced to a functioning interface test case, positioning the workflow for future kernel integration. Overall, the Phase IIA effort successfully transitioned the project from architectural proof-of-concept to externally exercised, solver-integrated capability, while clarifying remaining technical challenges related to standardization, performance optimization, and kernel integration.

42 ENGINEERING↗

Proton spin structure and generalized polarizabilities in the strong quantum chromodynamics regime

The strong interaction is not well understood at low energies or for interactions with low momentum transfer. Chiral perturbation theory gives testable predictions for the nucleonic generalized polarizabilities, which are fundamental quantities describing the nucleon’s response to an external field. We report a measurement of the proton’s generalized spin polarizabilities extracted with a polarized electron beam and a polarized solid ammonia target in the region where chiral perturbation theory is expected to be valid. The investigated structure function g2 characterizes the internal spin structure of the proton. From its moments, we extract the longitudinal–transverse spin polarizability δLT and twist-3 matrix element and polarizability d2¯. Our results provide discriminating power between existing chiral perturbation theory calculations and will help provide a better understanding of this strong quantum chromodynamics regime.

Ruth, D↗

Kaon mixing beyond the standard model with physical masses

We present nonperturbative results for beyond the standard model kaon mixing matrix elements in the isospin symmetric limit ( m u = m d ) of QCD, including a complete estimate of all dominant sources of systematic error. Our results are obtained from numerical simulations of lattice QCD with N f = 2 + 1 flavors of dynamical domain wall fermions. For the first time, these quantities are simulated directly at the physical pion mass m π ∼ 139 MeV for two different lattice spacings. We include data at three lattice spacings in the range a = 0.11 – 0.07 fm and with pion masses ranging from the physical value up to 450 MeV. Compared to our earlier work, we have added both direct calculations at physical quark masses and a third lattice spacing making the removal of discretization effects significantly more precise and eliminating the need for any significant mass extrapolation beyond the range of simulated data. We renormalize the lattice operators nonperturbatively using RI-SMOM off-shell schemes. These schemes eliminate the need to model and subtract nonperturbative pion poles that arises in the RI-MOM scheme and, since the calculations are performed with domain wall fermions, the unphysical mixing between chirality sectors is suppressed. Our results for the bag parameters in the MS ¯ scheme at 3 GeV are B K ≡ B 1 = 0.5240 ( 17 ) ( 54 ) , B 2 = 0.4794 ( 25 ) ( 35 ) , B 3 = 0.746 ( 13 ) ( 17 ) , B 4 = 0.897 ( 02 ) ( 10 ) and B 5 = 0.6882 ( 78 ) ( 94 ) , where the first error is from lattice uncertainties and the second is the uncertainty due to the perturbative matching to MS ¯ . Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Multiscale Modeling of the Mechanical Response of Silicon Carbide Composite Within the Accelerated Fuel Qualification Framework

The accelerated fuel qualification (AFQ) framework has been used for the initial development of multiscale modeling of silicon carbide (SiC) fiber reinforced composite (SiC-SiC). The AFQ framework provides a methodology to leverage physics-informed multiscale modeling along with a reduced set of empirical test data to reduce the time and cost of licensing and qualification of new nuclear fuel systems while maintaining the overall nuclear power plant safety case. SiC-SiC is being proposed for in-core applications, most notably fuel cladding, for current and next-generation nuclear reactors because of its high temperature stability, irradiation tolerance, and ability to withstand many accident conditions. As these composites exhibit multiscale architectures and complex microstructure-based fracture mechanics, it is an appealing use case for the AFQ methodology. While the end goal of this work is a single multiscale model that can be used for predictive in-core performance, current focus is on the individual various length scale models. Four individual models have been initially developed from microscale to engineering system level to capture key physics-based effects across different length scales. These models include a microscale homogenized tow model, a mesoscale fast Fourier transform–based weave model that integrates the homogenized tow model, a mesoscale finite element–based weave model, and a system-level BISON fuel performance model. Results of these models have undergone an initial comparison with separate-effects test data showing a good match to experimental results. By using the AFQ framework during model development, several near-term benefits have been secured including a reduction in development time for the SiC-SiC cladding, more targeted irradiation testing, and a better understanding of uncertainty.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Search for heavy neutral leptons decaying into muon-pion pairs in the MicroBooNE detector

We present upper limits on the production of heavy neutral leptons (HNLs) decaying to $\mu \pi$ pairs using data collected with the MicroBooNE liquid-argon time projection chamber (TPC) operating at Fermilab. This search is the first of its kind performed in a liquid-argon TPC. We use data collected in 2017 and 2018 corresponding to an exposure of $2.0 \times 10^{20}$ protons on target from the Fermilab Booster Neutrino Beam, which produces mainly muon neutrinos with an average energy of $\approx 800$ MeV. HNLs with higher mass are expected to have a longer time-of-flight to the liquid-argon TPC than Standard Model neutrinos. The data are therefore recorded with a dedicated trigger configured to detect HNL decays that occur after the neutrino spill reaches the detector. We set upper limits at the $90\%$ confidence level on the element $\lvert U_{\mu4}\rvert^2$ of the extended PMNS mixing matrix in the range $\lvert U_{\mu4}\rvert^2<(6.6$-$0.9)\times 10^{-7}$ for Dirac HNLs and $\lvert U_{\mu4}\rvert^2<(4.7$-$0.7)\times 10^{-7}$ for Majorana HNLs, assuming HNL masses between $260$ and $385$ MeV and $\lvert U_{e 4}\rvert^2 = \lvert U_{\tau 4}\rvert^2 = 0$.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Rethinking 𝛼−RuCl 3 : Parameters, models, and phase diagram

RuCl 3 was likely the first ever deliberately synthesized ruthenium compound, following the discovery of the 44 Ru element in 1844. For a long time it was known as an oxidation catalyst, with its physical properties being discrepant and confusing, until a decade ago when its allotropic form 𝛼−RuCl 3 rose to exceptional prominence. This “rediscovery” of 𝛼−RuCl 3 has not only reshaped the hunt for a material manifestation of the Kitaev spin liquid, but it has opened the floodgates of theoretical and experimental research in the many unusual phases and excitations that the anisotropic-exchange magnets as a class of compounds have to offer. Given its importance for the field of Kitaev materials, it is astonishing that the low-energy spin model that describes this compound and its possible proximity to the much-desired spin-liquid state is still a subject of significant debate ten years later. In the present study, we argue that the existing key phenomenological observations put strong natural constraints on the effective microscopic spin model of 𝛼−RuCl 3 , and specifically on its spin-orbit-induced anisotropic-exchange parameters that are responsible for the nontrivial physical properties of this material. These constraints allow one to focus on the relevant region of the multidimensional phase diagram of the 𝛼−RuCl 3 model, suggest an intuitive description of it via a different parametrization of the exchange matrix, offer a unifying view on the earlier assessments of its parameters, and bring closer together several approaches to the derivation of anisotropic-exchange models. We explore extended phase diagrams relevant to the 𝛼−RuCl 3 parameter space using quasiclassical, Luttinger-Tisza, exact diagonalization, and density-matrix renormalization-group methods, demonstrating a remarkably close quantitative accord between them on the general structure and hierarchy of the phases, with the zigzag, ferromagnetic, and incommensurate phases that are proximate to each other. As a result, one of the highlights is the detailed agreement on the nature of the incommensurate phases that realize two distinct counterrotating helical states.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Measurement of spin density matrix elements in Λ ( 1520 ) photoproduction at 8.2–8.8 GeV

We report on the measurement of spin density matrix elements of the $\Lambda(1520)$ in the photoproduction reaction $\gamma p\rightarrow \Lambda(1520)K^+$, via its subsequent decay to $K^{-}p$. The measurement was performed as part of the GlueX experimental program in Hall D at Jefferson Lab using a linearly polarized photon beam with $E_\gamma =$ 8.2-8.8 GeV. These are the first such measurements in this photon energy range. Results are presented in bins of momentum transfer squared, $-(t-t_\text{0})$. We compare the results with a Reggeon exchange model and determine that natural exchange amplitudes are dominant in $\Lambda(1520)$ photoproduction.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗