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 145 records · Page 8

An Alternative Conceptual Model for the Spent Nuclear Fuel–Water Interaction in Deep Geologic Disposal Conditions

For the long-term safety assessment of direct disposal of spent nuclear fuel in deep geologic repositories, knowledge on the radionuclide release rate from the UO 2 matrix is essential. This work provides a conceptual model to explain the results of leaching experiments involving used nuclear fuel or simulant materials in confirmed reducing conditions. Key elements of this model are: direct effect of radiation from radiolytic species (including defects and excited states) in the solid and in the first water layers in contact with its surface; and excess H 2 may be produced due to processes occurring at the surface of the spent fuel and in confined water volumes, which may also play a role in keeping the spent fuel surface in a reduced state. The implication is that the fractional radionuclide release rate used in most long-term safety assessments (10 -7 year -1 ) is over estimated because it assumes that there is net UO 2 oxidation caused by radiolysis, in contrast with the alternative conceptual model presented here. Furthermore, conventional water radiolysis models and radiation chemical yields published in the literature are not directly applicable to a heterogeneous system such as the spent fuel–water interface. Suggestions are provided for future work to develop more reliable models for the long-term safety assessment of spent nuclear fuel disposal.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Comparative study of He 6 β -decay based on different similarity-renormalization-group evolved chiral interactions

We report on a study of the Gamow-Teller matrix element contributing to ${}^6$He $\beta$-decay with similarity renormalization group (SRG) versions of momentum- and configuration-space two-nucleon interactions. These interactions are derived from two different formulations of chiral effective field theory ($\chi$EFT) -- without and with the explicit inclusion of $\Delta$-isobars. We consider evolution parameters $\Lambda_{\rm SRG}$ in the range between 1.2 and 2.0 fm$^{-1}$ and, for the $\Delta$-less case, also the unevolved (bare) interaction. The axial current contains one- and two-body terms, consistently derived at tree level (no loops) in the two distinct $\chi$EFT formulations we have adopted here. The ${}^6$He and ${}^6$Li ground-state wave functions are obtained from hyperspherical-harmonics (HH) solutions of the nuclear many-body problem. In $A\,$=$\,6$ systems, the HH method is limited at present to treat only two-body interactions and non-SRG evolved currents. Our results exhibit a significant dependence on $\Lambda_{\text{SRG}}$ of the contributions associated with two-body currents, suggesting that a consistent SRG-evolution of these is needed in order to obtain reliable estimates. We also show that the contributions from one-pion-exchange currents depend strongly on the model (chiral) interactions and on the momentum- or configuration-space cutoffs used to regularize them. These results might prove helpful in clarifying the origin of the sign difference recently found in No-Core-Shell-Model and Quantum Monte Carlo calculations of the ${}^6$He Gamow-Teller matrix element.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Extraction of low-energy constants of single- and double-β decays from lattice QCD: A sensitivity analysis

Lattice quantum chromodynamics (LQCD) has the promise of constraining low-energy constants (LECs) of nuclear effective field theories (EFTs) from first-principles calculations that incorporate the dynamics of quarks and gluons. Given the Euclidean and finite-volume nature of LQCD outputs, complex mappings are developed in recent years to obtain the Minkowski and infinite-volume counterparts of LQCD observables. In particular, as LQCD is moving toward computing a set of important few-nucleon matrix elements at the physical values of the quark masses, it is important to investigate whether the anticipated precision of LQCD spectra and matrix elements will be sufficient to guarantee tighter constraints on the relevant LECs than those already obtained from phenomenology, considering the nontrivial mappings involved. With a focus on the leading-order LECs of the pionless EFT, L 1,A and $g$$^{NN}_{v}$, which parametrize, respectively, the strength of the isovector axial two-body current in a single-β decay (and other related processes such pp fusion), and of the isotensor contact two-body operator in the neutrinoless double-β decay within the light neutrino exchange scenario, the expected uncertainty on future extractions of L 1,A and $g$$^{NN}_{v}$ are examined using synthetic data at the physical values of the quark masses. It is observed that achieving small uncertainties in L 1,A will be challenging, and (sub)percent-level precision in the two-nucleon spectra and matrix elements is essential in reducing the uncertainty on this LEC compared to the existing constraints. On the other hand, the short-distance coupling of the neutrinoless double-β decay, $g$$^{NN}_{v}$, is shown to be less sensitive to uncertainties on both LQCD energies and the matrix element, and can likely be constrained with percent-level precision in the upcoming LQCD calculations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Summary of light-element evaluation work: nuclear data uncertainty quantification

Hale & Paris (T-2) have evaluated the 15 N system to provide nuclear cross section and covariance information for nuclear data uncertainty quantification. This memo gives a brief discussion of the R-matrix method used to generate this nuclear scattering and reaction data in the ENDF-6 format.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Advanced Ab Initio Methods for Nuclear Structure (Final Report)

Over the past decade, there has been enormous progress in the description of nuclear structure from first principles, using interactions from Chiral Effective Field Theory that are rooted in Quantum Chromodynamics, the fundamental theory of strong interactions, and many-body methods that solve the Schrödinger equation with systematically improvable approximations. Amongst those are the family of In-Medium Similarity Renormalization Group (IMSRG) framework developed by the PI and his co-workers. While these methods scale polynomially in the size N of the single-particle basis, computational efforts still grows dramatically as we increase the degrees of freedom for the nucleons by relaxing symmetries or introducing continuum couplings (see below), or as we push to improved truncations to provide precise inputs for experimental efforts, in particular in fundamental symmetry searches. In order to address this growing computational cost, one focus area of this award was the exploration of compression and factorization methods. The key to success or failure is the presence of low-rank structures within the matrix elements of NN and 3N interactions or the IMSRG evolution operator, and a means to reformulate the method that will let us exploit them.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Diabatic Hamiltonian matrix elements made simple

With a view to applying the generator coordinate method to large configuration spaces, we propose a simple approximate formula to compute diabatic many-body matrix elements without having to evaluate two-body interaction matrix elements. Here, the method is illustrated with two analytically solvable Hamiltonians based on the harmonic oscillator.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Statistical uncertainty of fission matrix eigenvalues using perturbation theory

Eigenvalue search of high dominance ratio systems may be slow to converge. The fission matrix element is defined by its element (FM){sub ij}, which are the probability for a neutron born in cell i to create a fission in cell j for a spatial mesh of n{sub i}*n{sub j} elements. Fission matrices are used in Monte Carlo criticality simulations to enhance computing speed, but also to find higher order eigenvalues. However, few studies have been made on the link between statistical uncertainties of fission matrix elements and eigenvalues uncertainties. Thus, dominance ratio statistical uncertainties remain unknown. This paper uses a new generalized perturbation theory (GPT) method to estimate sensitivities of eigenvalues to fission matrix elements and then to calculate dominance ratio uncertainties.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Uncertainty quantification in elastic constants of SiC f /SiC m tubular composites using global sensitivity analysis

Silicon carbide fiber and silicon carbide matrix (SiC f /SiC m ) tubes produced through the chemical vapor infiltration process have become a candidate cladding material in nuclear applications. The performance of this composite is influenced by many variables such as braiding angle, porosity, material properties, etc., which vary over a range of values due to the inherent fluctuations in the manufacturing process. In this study, the variability in elastic constants of SiC f /SiC m composite has been quantified through multiscale finite element (FE) simulations, variable screening, and high-fidelity surrogate modeling. The key variables dominantly affecting the elastic constants of SiC f /SiC m tubes were identified using global sensitivity analysis. A surrogate to the high-fidelity FE-based model was used in Monte Carlo simulations to generate a hundred thousand samples from which the uncertainty in elastic constants was assessed. It turned out that the coefficient of variation was less than 10%.

Materials Science↗

Suppressed electric quadrupole collectivity in Si 32

Lying between O 16 and Ca 40 , the s d shell is well described by robust phenomenological and nuclear theories. In this work, however, we highlight an unexplained reduction in electric-quadrupole strength in the rare isotope Si 32 , studied through sub-barrier Coulomb excitation. It is found that the oblate nature of the deformation is well reproduced, while the absolute scale of quadrupole deformation, however, is inhibited by approximately a factor of 2 compared to theoretical predictions. Through comparison with shell-model and calculations, we present a number of possible explanations for this inhibited E 2 strength. By comparing the results of these calculations to multiple observables, we conclude that there is a reduced role for out-of-space excitations in Si 32 , resulting in a reduction in the corrections normally applied to both models. Published by the American Physical Society 2024

20 ≤ A ≤ 38↗

Spectroscopy of states in 136 Ba using the 138 Ba(p,t) reaction

Background: The 136 Ba isotope is the daughter nucleus in 136 Xe $ββ$ decay. It also lies in a shape transitional region of the nuclear chart, making it a suitable candidate to test a variety of nuclear models. Purpose: To obtain spectroscopic information on states in 136 Ba, which will allow a better understanding of its low-lying structure. These data may prove useful to constrain future 136 Xe → 136 Ba neutrinoless $ββ$ decay matrix element calculations. Methods: Here, a 138 Ba(p, t) reaction was used to populate states in 136 Ba up to approximately 4.6 MeV in excitation energy. The tritons were detected using a high-resolution Q3D magnetic spectrograph. A distorted wave Born approximation analysis was performed for the measured triton angular distributions. Results: 102 excited states in 136 Ba were observed, out of which 52 are reported for the first time. Definite spin-parity assignments are made for 26 newly observed states, while previously ambiguous assignments for ten other states are resolved. Together with other available data, the results are used to determine level densities in 136 Ba. These were compared with theory predictions, obtained using shell model calculations with Hamiltonians previously used for 136 Xe neutrinoless $ββ$ decay matrix element evaluations. Conclusions: The shell model predicted level densities agree reasonably well for the two Hamiltonians. However the results for theory and experiment are found to agree only at lower energies, diverging from one another for the higher lying states, with the discrepancy increasing with energy. This is presumably because of lower production cross sections for a majority of the higher-lying predicted states and the experimental limitations in resolving a large number of nearly degenerate states predicted by the theory.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Absorber Clamp for Microcalorimeter Decay Energy Spectrometry

Microcalorimeter Decay Energy Spectrometry (DES) is of interest to nuclear safeguards due to its ability to provide high precision isotopic compositions of nanogram-to microgram-scale samples of Pu and U and related daughter products. The DES method is able to record decay energy of each alpha-decay event in a sample that is embedded in a metal matrix (absorber) and thermally linked to a microcalorimeter detector. This work optimizes the DES technique used to thermally link the absorber and microcalorimeter detector element to allow for more rapid assembly and to increase detector performance and operating life. Optimized attachment methods are crucial for enhancing the viability of DES in high-sample-throughput facilities, such as those that support international nuclear safeguards measurements. Here, in this study, we designed and implemented a pressure-based absorber clamp and evaluated the performance of this new attachment method relative to pressed indium bond attachment. Results using the absorber clamp demonstrate a streamlined detector assembly procedure that minimizes accidental damage to detectors, as well as increasing detector pulse speeds by 57%. Spectral comparison shows the clamp preserves detector performance relative to indium attachment.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Impact of heavy quark and quarkonium data on nuclear gluon PDFs

A clear understanding of nuclear parton distribution functions (nPDFs) plays a crucial role in the interpretation of collider data taken at the Relativistic Heavy Ion Collider, the Large Hadron Collider (LHC), and in the near future at the Electron-Ion Collider. Even with the recent inclusions of vector boson and light meson production data, the uncertainty of the gluon PDF remains substantial and limits the interpretation of heavy ion collision data. To obtain new constraints on the nuclear gluon PDF, we extend our recent nCTEQ15WZ+SIH analysis to inclusive quarkonium and open heavy-flavor meson production data from the LHC. This vast new data set covers a wide kinematic range and puts strong constraints on the nuclear gluon PDF down to x ≲ 10 -5 . The theoretical predictions for these data sets are obtained from a data-driven approach, where proton-proton data are used to determine effective scattering matrix elements. This approach is validated with detailed comparisons to existing next-to-leading order calculations in nonrelativistic QCD for quarkonia and in the general-mass variable-flavor-number scheme for the open heavy-flavored mesons. In addition, the uncertainties from the data-driven approach are determined using the Hessian method and accounted for in the PDF fits. This extension of our previous analyses represents an important step toward the next generation of PDFs not only by including new data sets, but also by exploring new methods for future analyses.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

One-loop matching for quark dipole operators in a gradient-flow scheme

The quark chromoelectric dipole (qCEDM) operator is a CP-violating operator describing, at hadronic energies, beyond-the-standard-model contributions to the electric dipole moment of particles with nonzero spin. In this paper we define renormalized dipole operators in a regularization-independent scheme using the gradient flow, and we perform the matching at one loop in perturbation theory to renormalized operators of the same and lower dimension in the more familiar MS scheme. We also determine the matching coefficients for the quark chromo-magnetic dipole operator (qCMDM), which contributes for example to matrix elements relevant to CP-violating and CP-conserving kaon decays. The calculation provides a basis for future lattice QCD computations of hadronic matrix elements of the qCEDM and qCMDM operators.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Understanding the structural and morphological effects of synthesis route on NpO 2

The availability of actinide standard materials for use in nuclear safeguard applications is critical, as is thorough characterization thereof. Although accurate trace element compositions and isotopic considerations are paramount for deployment of reference standards, structural characterization is also essential towards accurately describing the chemical form and potential matrix effects in candidate materials. Here, to this end, samples of NpO 2 were synthesized via a direct denitration (DD) method and probed with powder X-ray diffraction (PXRD), Raman spectroscopy, and scanning electron microscopy (SEM) for structural and morphological characterization and comparison with NpO 2 materials produced via modified direct denitration (MDD). PXRD confirmed the bulk identity of NpO 2 , and no additional phases were identified using this method. Analysis of Raman data collected using a 532 nm excitation wavelength indicates that samples are mostly phase pure; however, some variability in spectral features is observed. Analysis of additional spectroscopic data collected with a 785 nm excitation wavelength revealed variability in the relative intensity of spectral features. Raman spectroscopy indicates that the sample is primarily NpO 2 ; however, additional signals indicate possible structural disorder, oxidized species, or potential contributions from other Np phases. To further investigate the possibility of additional phase contributions within the sample of NpO 2 , Raman spectroscopic mapping was employed to examine the homogeneity of the sample produced via DD. From this analysis, we determined that despite variability in the intensity of Raman-active vibrational modes, consistent spectra are obtained throughout the area of the sample investigated. SEM images show aggregates with variable sizes and shapes, with rounded, primary particles possessing an average diameter of approximately 100 nm. Comparison of the results of these multimodal analyses to the literature indicates that the crystal chemical, spectroscopic, and microstructural properties of NpO 2 vary based on synthesis method, even if X-ray diffraction data indicate that the bulk phase is NpO 2 .

Direct denitration↗

Construction of St. Benedict

The Superallowed Transition Beta-Neutrino Decay Ion Coincidence Trap (St. Benedict) is currently under construction at the University of Notre Dame Nuclear Science Laboratory. It is designed to measure the beta-neutrino angular correlation parameter in superallowed mixed beta-decay transitions between mirror nuclei in order to extract the Fermi-to-Gamow Teller mixing ratio and test theoretical corrections entering in the determination of the V ud element of the Cabibbo–Kobayashi–Maskawa quark mixing matrix. St. Benedict includes a large volume gas catcher to thermalize the radioactive ion beam from TwinSol, a differentially-pumped extraction system, a radiofrequency quadrupole ion cooler and buncher, and a Paul trap for the observation and detection of the β decay products. Finally, the status of each of these components will be presented.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Detection of Br and I as atomic anions using liquid sampling – atmospheric pressure glow discharge/Orbitrap mass spectrometry

The quantification and determination of halogen isotope composition is essential in many fields including geochemistry and nuclear forensics. Detection of the halogen elements is commonly attempted with inductively coupled plasma mass spectrometry (ICP-MS) however, there are multiple challenges faced. Most significantly this includes very poor ionization efficiency and the potential for isobaric interferences from either matrix or plasma species. Thus, sensitivity and accuracy are often limiting issues. Here a liquid sampling-atmospheric pressure glow discharge (LS-APGD) ionization source is coupled with an ultrahigh resolution Orbitrap MS to perform halogen detection of bromine and iodine as initial analytes. This facilitates simple and sensitive detection of these elements as atomic anions from simple aqueous salt solutions. Collision-induced dissociation (CID) and higher-energy collisional dissociation (HCD) energies were optimized to produce the maximum response of Br − and I − . The LS-APGD conditions were optimized using a design of experiments (DOE) approach for Br − and I − concurrently. Response curves for Br − and I − solutions determined limit of detection (LOD) values of 50 pg and 5 pg, respectively, in 20 μL aliquots. The curves indicated response factors of 0.67 for 79 Br − and 0.90 for 127 I − . Br isotope ratios were determined with precision between 0.7 and 7.3% RSD, with the isotope ratios determined with precision 0.8% RSD at concentrations above the limit of quantification. Preliminary tests were conducted to evaluate the effect of different cations (Na + , K + , and Mg 2+ ) on Br − responses, with no discernible impacts observed on the halogen signal responses. This study demonstrates the potential use of the LS-APGD-Orbitrap-MS for the detection of multiple halogens while avoiding interferences, minimizing sample preparation, and overcoming ionization barriers.

Lang, Dehlia A. [Clemson University, SC (United St↗

Gamma-ray spectrometry analysis methods for radioisotope dissolution and mixing for nuclear forensics applications

We report the production of surrogate nuclear explosive debris involves the dissolution of irradiated fissile material, elemental isolation, and recombination based on required fission product ratios. Gamma-ray spectrometry is used following each step to assess the radioactive fission product inventory. A high-purity germanium detector is used to measure each sample, and the detector efficiencies are determined prior to sample measurements. Through radiochemical techniques, elements are isolated and recombined into samples with specific fission product ratios. As a result, a matrix-based approach can account for imperfect chemical separations to calculate the proper mixture volumes to obtain intended volatile-to-refractory fission product ratios, or R-values.

07 ISOTOPE AND RADIATION SOURCES↗