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 127 records · Page 7

Computing the QRPA level density with the finite amplitude method

Here, we describe a new algorithm to calculate the vibrational nuclear level density of an atomic nucleus. Fictitious perturbation operators that probe the response of the system are generated by drawing their matrix elements from some probability distribution function. We use the Finite Amplitude Method to explicitly compute the response for each such sample. With the help of the Kernel Polynomial Method, we build an estimator of the vibrational level density and provide the upper bound of the relative error in the limit of infinitely many random samples. The new algorithm can give accurate estimates of the vibrational level density. Since it is based on drawing multiple samples of perturbation operators, its computational implementation is naturally parallel and scales like the number of available processing units.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

CSBCOO SPMM: Compressed Sparse Block Coordinate format sparse matrix dense vector multiply benchmark (CSBCOO SPMM) v1.0

The purpose of this microbenchmark is to provide a means to explore different programming methodologies using a simple, but not trivial, mathematical kernel. The kernel is based on the matrix vector product in the MFDn nuclear configuration interaction code's LOBPCG eigensolver. This particular mathematical operation is a key element in the solution of systems of equations and block eigensolvers.

Cook, Brandon↗

A Study of a Surrogate Nuclear Thermal Propulsion Fuel Element

We report a surrogate fuel, Zirconium Carbide/Titanium Nitride has been prepared using a simple method to obtain a more uniform distribution of titanium nitride in the zirconium carbide matrix. Powders were sintered using Spark Plasma Sintering at 1750oC and 1850oC using a modified sintering profile. This resulted in a fairly uniform distribution of titanium nitride in the zirconium carbide matrix with a density of 95% of theoretical at 1850oC. Scanning electron microscopy, indentation hardness and x-ray diffraction results are given.

36 MATERIALS SCIENCE↗

An IIPG-based finite element framework in MOOSE for modeling fiber reinforced composite failure governed by extrinsic cohesive laws

Fiber reinforced composites (FRC) are widely applied to many important engineering structures such as silicon carbide ceramic matrix composite (SiC CMC) of nuclear fuel claddings. Many material and structure failures are initiated at the interface between two distinctive materials, such as the debonding of interfaces of matrix/fiber of composites or bulk matrix/inclusion in metals at micro-scales. Accurate and robust modeling of these interface failures is critical to the development of novel material and structure with improved functionality. This paper presents a discontinuous Galerkin (DG) finite element method based on the Incomplete Interior Penalty scheme and its implementation in MOOSE for modeling material interface behaviors. The proposed method naturally accommodates important extrinsic cohesive zone material models and is capable of modeling the failures of material interfaces across which material properties are highly distinctive. We demonstrate a good performance of the proposed method using single fiber and multiple fiber debonding problems.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Multi-channel direct detection of light dark matter: theoretical framework

We present a unified theoretical framework for computing spin-independent direct detection rates via various channels relevant for sub-GeV dark matter — nuclear re- coils, electron transitions and single phonon excitations. Despite the very different physics involved, in each case the rate factorizes into the particle-level matrix element squared, and an integral over a target material- and channel-specific dynamic structure factor. We show how the dynamic structure factor can be derived in all three cases following the same procedure, and extend previous results in the literature in several aspects. For electron transitions, we incorporate directional dependence and point out anisotropic target materials with strong daily modulation in the scattering rate. For single phonon excitations, we present a new derivation of the rate formula from first principles for generic spin-independent couplings, and include the first calculation of phonon excitation through electron couplings. We also discuss the interplay between single phonon excitations and nuclear recoils, and clarify the role of Umklapp processes, which can dominate the single phonon production rate for dark matter heavier than an MeV. Our results highlight the complementarity between various search channels in probing different kinematic regimes of dark matter scattering, and provide a common reference to connect dark matter theories with ongoing and future direct detection experiments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Axions in neutron star mergers

Supernovae and cooling neutron stars have long been used to constrain the properties of axions, such as their mass and interactions with nucleons and other Standard Model particles. We investigate the prospects of using neutron star mergers as a similar location where axions can be probed in the future. We examine the impact axions would have on mergers, considering both the possibility that they free-stream through the dense nuclear matter and the case where they are trapped. We calculate the mean free path of axions in merger conditions, and find that they would free-stream through the merger in all thermodynamic conditions. In contrast to previous calculations, we integrate over the entire phase space while using a relativistic treatment of the nucleons, assuming the matrix element is momentum-independent. In particular, we use a relativistic mean field theory to describe the nucleons, taking into account the precipitous decrease in the effective mass of the nucleons as density increases above nuclear saturation density. We find that within current constraints on the axion-neutron coupling, axions could cool nuclear matter on timescales relevant to neutron star mergers. Our results may be regarded as first steps aimed at understanding how axions affect merger simulations and potentially interface with observations.

79 ASTRONOMY AND ASTROPHYSICS↗

Secondary γ-ray decays from the partial-wave $\textit{T}$ matrix with an $\textit{R}$-matrix application to 15 N($p, α_1γ$) 12 C

The secondary γ rays emitted following a nuclear reaction are frequently relatively straightforward to detect experimentally. Despite the large volume of such data, a practical formalism for describing these γ rays in terms of partial-wave $\textit{T}$-matrix elements has never been given. The partial-wave formalism is applicable when $\textit{R}$-matrix methods are used to describe the reaction in question. This paper supplies the needed framework, and it is demonstrated by the application to the 15 N($p, α_1γ$) 12 C reaction.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

The convergence of quadrupole rotational invariants from the nuclear shell model

Nuclei exhibit both single-particle and collective degrees of freedom, with the latter often subdivided into vibrational and rotational motions. Experimentally identifying the relative roles of these collective modes is extremely challenging, particularly in the face of possible shape coexistence. Model-independent, invariant quantities describing the deformation of a nucleus in the intrinsic frame have long been known but their determination potentially requires a large quantity of experimental data to achieve convergence. Through comparison with the nuclear shell model, the question of convergence is addressed. Shell-model calculations performed in the sd- and pf-shell model spaces are used to determine electric-quadrupole matrix elements for a multitude of low-lying states using the first 40 states of the relevant spins. Relative contributions to the rotationally invariant quantities from multiple states can therefore be determined. It is found that, on average, the inclusion of four intermediate states results in the leading-order invariant, $\langle\hat{Q^2}\rangle$, converging to within 10% of its true value and the triaxiality term, cos (3δ), converging to its true value, though some variance remains. Higher-order quantities relating to the softness of the nuclear shape are found to converge more slowly. The convergence of quadrupole rotationally invariant sum rules was quantified in the sd- and pf-shell model spaces and indicates the challenge inherent in a full determination of nuclear shape. The present study is limited to relatively small valence spaces. Finally, larger spaces, such as the rare-earth region, potentially offer faster convergence.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Weak Transitions in Light Nuclei

Nuclei are used for high-precision tests of the Standard Model and for studies of physics beyond the Standard Model. Without a thorough understanding of nuclei, we will not be able to meaningfully interpret the growing body of experimental data nor will we be able to disentangle new physics signals from underlying nuclear effects. This calls for accurate calculations of nuclear structure and reactions. In this work, we focus on electroweak decays in nuclei with mass number A ≤ 10 and report on ab initio Quantum Monte Carlo calculations of reduced matrix elements entering beta decays and electron captures in nuclei with mass number A ≤ 10. The many-body wave functions are calculated using selected Norfolk two- and three-nucleon potential models and associated one- and two-body axial currents at tree-level obtained from a chiral effective field theory with pions, nucleons, and Δ. The agreement with the experimental data is satisfactory except for transitions in A = 8 nuclei. In this specific case, the theory significantly underpredicts the experimental data, which indicates the need of further improvements in the corresponding nuclear wave functions. In this study, emphasis is placed on the contributions of two-body axial currents that are carefully analyzed using two-body transition densities. This allow us to study the spatial distribution and short-range behavior of two-body dynamics. In particular, the transition densities when scaled to peak at 1.0 exhibit universal short-range behavior across the considered nuclei, while they differ in the long-range tails.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Scattering Observables from One- and Two-body Densities: Formalism and Application to $\pmb \gamma $ ${}^3\hbox {He}$ Scattering

We introduce the transition-density formalism, an efficient and general method for calculating the interaction of external probes with light nuclei. One- and two-body transition densities that encode the nuclear structure of the target are evaluated once and stored. They are then convoluted with an interaction kernel to produce amplitudes, and hence observables. Here, by choosing different kernels, the same densities can be used for any reaction in which a probe interacts perturbatively with the target. The method therefore exploits the factorisation between nuclear structure and interaction kernel that occurs in such processes. We study in detail the convergence in the number of partial waves for matrix elements relevant in elastic Compton scattering on 3 He. The results are fully consistent with our previous calculations in Chiral Effective Field Theory. But the new approach is markedly more computationally efficient, which facilitates the inclusion of more partial-wave channels in the calculation. We also discuss the usefulness of the transition-density method for other nuclei and reactions. Calculations of elastic Compton scattering on heavier targets like 4 He are straightforward extensions of this study, since the same interaction kernels are used. And the generality of the formalism means that our 3 He densities can be used to evaluate any 3 He elastic-scattering observable with contributions from one- and two-body operators.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Ab initio estimation of E 2 strengths in Li 8 and its neighbors by normalization to the measured quadrupole moment

For electric quadrupole (E2) observables, which depend on the large-distance tails of the nuclear wave function, ab initio no-core configuration interaction (NCCI) calculations converge slowly, making meaningful predictions challenging to obtain. Nonetheless, the calculated values for different E2 matrix elements, particularly those involving levels with closely-related structure (e.g., within the same rotational band) are found to be robustly proportional. This observation suggests that a known value for one observable may be used to determine the overall scale of E2 strengths, and thereby provide predictions for others. In particular, we demonstrate that meaningful predictions for E2 transitions may be obtained by calibration to the ground-state quadrupole moment. Here, we test this approach for well-measured low-lying E2 transitions in 7 Li and 9 Be, then provide predictions for transitions in 8 Li and 9 Li. In particular, we address the 2 + → 1 + transition in 8 Li, for which the reported measured strength exceeds ab initio Green's function Monte Carlo (GFMC) predictions by over an order of magnitude.

6 ≤ A ≤ 19↗

Nuclear structure and band mixing in Pt 194

We introduce a two-particle, two-hole mixed configuration scheme to fit E2 strengths for the 0 ↔ 2, 2 ↔ 4, and 4 ↔ 6 transitions in 194 Pt. The interaction includes two sets of pairing operators, {S ± (t), S 0 (t)} (t = s, d). Solutions within this framework are used to analyze energy spectra, E2 transitions, and band-mixing features of the model. The results confirm that mixing is small and similar for J = 2, 4, and 6, with the calculated energies and transition matrix elements in excellent agreement with experimental data.

190 ≤ A ≤ 219↗

Quasiclassical Gluon Fields and Low’s Soft Theorem at Small Momentum-Fraction x

In the high-energy limit, soft gluons can be approximately described by quasiclassical gluon fields. It is well known that the gluon field is a pure gauge field on the transverse plane at eikonal order. We derived the complete next-to-eikonal order solutions of the classical Yang-Mills equations for soft gluons in the dense nuclear regime. Utilizing these solutions, it is shown that Low’s soft theorem at small x can be obtained by considering off-diagonal matrix elements of quasiclassical chromoelectric field between single-gluon states in the dilute regime. We further propose on extending Low’s soft theorem at small x to incorporate the effects of gluon saturation in the dense regime.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Evaluation of Light-Element Reactions in the Resolved Resonance Region

Light-element reactions at low energies in the resolved resonance region are important for a range of applications in basic and applied sciences including nuclear reactors, nonproliferation, cultural heritage, forensics and environmental control, rare event investigations and nuclear astrophysics. In this paper, we report on an effort to evaluate charged-particle cross sections in the resolved resonance region and produce evaluated nuclear data files for further processing and inclusion in evaluated data libraries. We discuss the open issues in R-matrix calculations as we extend to higher energies, such as dealing with the rapidly growing number of open channels and merging with the regime of smooth cross sections described by the statistical model, and present attempts to address these issues in neutron-induced reactions relevant to nuclear reactor applications.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

3D TRISO particle-explicit compact meshing

The TRI-structural ISOtropic (TRISO) layered fuel particle is a robust nuclear fuel form offering enhanced safety and performance for advanced reactor concepts, including high-temperature gas-cooled reactors and other Generation IV designs. These poppy-seed-sized particles are embedded in a graphite matrix to form fuel elements that must withstand elevated temperatures and high burn-up levels. The heterogeneous nature of these fuel elements — comprising thousands of randomly distributed TRISO particles — produces complex stress fields and thermal gradients that one- and two-dimensional models cannot accurately capture. While three-dimensional modeling has improved predictions of dimensional changes, internal pressure buildup, and fission product transport under irradiation, current approaches rely on homogenized material properties that are known to have considerable divergence from experimental observations. This work presents a methodology for optimized random packing of TRISO fuel compacts and full three-dimensional mesh generation within the BISON fuel performance code, with each particle coating layer individually discretized. The resulting mesh was demonstrated through heat conduction simulations under representative in-reactor operating conditions, showing strong agreement with expected behavior. This capability enables detailed analysis of particle-to-particle interactions, matrix cracking mechanisms, and the statistical distribution of coating layer failures — all of which directly govern fuel performance and safety margins.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗