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 325 records · Page 18

From suppressed void growth to significant void swelling in NiCoFeCr complex concentrated solid-solution alloy

Void swelling can result in dimensional instability and undermine the safe operation of nuclear reactors. Current strategies to inhibit void swelling mainly focus on enhancing defect absorption and recombination by introducing high-density defect sinks. Complex concentrated solid-solution alloys (CSAs), including high-entropy alloys, can withstand severe radiation damage due to their inherent chemical complexity without interfaces. However, the underlying mechanisms for void suppression in CSAs are far from clear. In this research, we studied the void evolution with respect to irradiation depths, doses, and temperatures in equiatomic NiCoFeCr under 3 MeV Ni ion irradiations. At relatively low doses (16 and 54 displacements per atom, dpa), voids form mainly outside of the ion-damaged region, and void formation in the peak damage region is suppressed, leading to negligible swelling. However, with further increase of dose (86 up to 250 dpa), significant void growth occurs in the peak damage region and extended dislocation lines dominate instead of short dislocation lines and loops formed at lower doses. From 500 to 700 °C, the dislocation density decreases while dislocations grow. Although the overall void swelling increases dramatically from 500 to 580 °C at 54 dpa, void growth in the peak damage region is still suppressed. The transition from suppressed void growth to significant void swelling is attributed to dislocation evolution and local chemical inhomogeneity (enrichment of Fe/Cr in the matrix) at higher doses. Our study shows that controlling element diffusion and defect evolution through tuning chemical complexity can further enhance the swelling resistance of CSAs.

36 MATERIALS SCIENCE↗

Detailed analysis of excited-state systematics in a lattice QCD calculation of 𝑔 𝐴

Excited state contamination remains one of the most challenging sources of systematic uncertainty to control in lattice QCD calculations of nucleon matrix elements and form factors: early time separations are contaminated by excited states and late times suffer from an exponentially bad signal-to-noise problem. High-statistics calculations at large time separations ≳ 1 fm are commonly used to combat these issues. In this work, focusing on g A , we explore the alternative strategy of utilizing a large number of relatively low-statistics calculations at short to medium time separations (0.2–1 fm), combined with a multistate analysis. On an ensemble with a pion mass of approximately 310 MeV and a lattice spacing of approximately 0.09 fm, we find this provides a more robust and economical method of quantifying and controlling the excited state systematic uncertainty. A quantitative separation of various types of excited states enables the identification of the transition matrix elements as the dominant contamination. The excited state contamination of the Feynman-Hellmann correlation function is found to reduce to the 1% level at approximately 1 fm while, for the more standard three-point functions, this does not occur until after 2 fm. Critical to our findings is the use of a global minimization, rather than fixing the spectrum from the two-point functions and using them as input to the three-point analysis. We find that the ground state parameters determined in such a global analysis are stable against variations in the excited state model, the number of excited states, and the truncation of early-time or late-time numerical data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Investigation of post-breakup Coulomb acceleration using a trajectory model

Intermediate mass fragments ejected during the deexcitation of excited projectilelike fragments may promptly decay following ejection; the daughter particles that are subsequently produced are subject to interactions with the residual nucleus that affect final-state observables, a process herein referred to as post-breakup Coulomb acceleration. A simple classical Coulomb interaction model was used to study modification of 8 Be (2 + ), 5 Li (3/2 – ), 7 Li (7/2 – ), 7 Be (7/2 – ), and states in 12 B emitted following heavy-ion collisions of 28 Si + 12 C at 35 MeV/nucleon. Here, in contrast to previous work studying 8 Be (2 + ), excellent agreement between simulation and experiment was obtained using only Coulomb forces when either a Lorentzian or R-matrix line shape was used to describe the initial relative energy rather than a Gaussian. In consideration of the obtained results, improvements to the model and evaluation of experimental data are discussed as future directions, but it was concluded that the effects observed in the present data can be accurately described using only elements of classical mechanics and that the process is largely understood for a wide range of state lifetimes and mass and charge (a)symmetries.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Direct and Indirect Measurements for a Better Understanding of the Primordial Nucleosynthesis

The Big Bang Nucleosynthesis (BBN) model is a great success of nuclear astrophysics due to the outstanding agreement between observational and predicted light elements abundances. One exception, however, is the so-called “lithium problem.” In this context, experimental efforts to measure the relevant reactions have been brought to an increased level of accuracy in measuring primordial abundances, and the introduction of indirect methods has allowed researchers to overcome the natural limitations of direct measurements in the energy range of interest for BBN. Here we review the results obtained from the application of the Trojan Horse Method to some of the most influential reactions of the standard network, such as 2 H(d,p) 3 H, 2 H(d,n) 3 He, 3 He(d,p) 4 He, 7 Li(p,α) 4 He, and 7 Be(n,α) 4 He. The relevant cross sections have been then used as new inputs to a classical BBN code, resulting in important constraints that make suggestions for a possible solution for the lithium problem outside of nuclear physics.

79 ASTRONOMY AND ASTROPHYSICS↗

Microstructural analysis of tristructural isotropic particles in high-temperature steam mixed gas atmospheres

High-temperature gas-cooled reactors (HTGRs) use tristructural isotropic (TRISO) particles embedded in a graphitic matrix material to form the integral fuel element. Potential off-normal reactor conditions for HTGRs include steam ingress with temperatures above 1,000 °C. Fuel element exposure to steam can cause the graphitic matrix material to evolve, forming an atmosphere composed of oxidants and oxidation products and potentially exposing the TRISO particles to these conditions. Investigating the oxidation response of TRISO particles exposed to a mixed gas atmosphere will provide insight into the stability under off-normal conditions. In this study, surrogate TRISO particles were exposed to high temperatures (T = 1,200 °C) in flowing steam (3% < pH 2 O < 21%) and CO (pCO < 1%) to determine the oxidation behavior of the SiC layer when exposed to various mixed gas atmospheres. Scanning electron microscopy, x-ray diffraction, and focused ion beam milling was used to determine the impact of CO and steam on the oxidation behavior of the SiC layer. Therefore, the data presented demonstrates how the SiC layer showed strong oxidation resistance due to limited SiO 2 growth and maintained its structural integrity under these off-normal conditions.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Development of Direct Injection/Ionization Mass Spectrometry Methods for Whole Molecule Characterization

The objective of this work is to adapt ambient ionization mass spectrometry (AMS) techniques for the rapid analysis of intact uranium complexes, stable strontium, cerium, and explosive compounds. The methods used were “soft ionization” techniques, which facilitate the detection of whole molecule complexes. The soft ionization mass spectrometry (MS) techniques that were investigated include paper spray ionization (PSI), matrix-assisted ionization (MAI), electrospray ionization (ESI) and direct analysis in real time (DART). For the first time, PSI-MS methods were successfully developed for whole molecule uranium-containing analytes (uranyl acetate, uranyl nitrate, and uranyl-tributylphosphate complexes). This was also the first demonstration of uranium complex detection and characterization and one of the few examples of inorganic analysis using MAI techniques. Proof of concept experiments also putatively identified matrix-derived ions and ion complexes that have not previously been described in the literature. Additionally, PSI-MS on cotton swipe samples doped with a multi-element standard containing μg levels of U, Bi, Pb, Cd, Fe, and Zn were directly analyzed without purification, representing a major improvement over existing methods. Both PSI and MAI methods demonstrated limits of detection (LODs) in the 10 - 100’s ng for various uranyl species within a range of 10’s ppm - 100’s ppb, dependent on analytical method and analyte species. AMS methods were also developed for other inorganics, including Ce and Sr, and organic explosive residues to address specific challenges in environmental monitoring and forensics. Further refinement and qualification of the AMS techniques developed within this effort would lead to significant cost reduction and timeliness by facilitating the triage and queueing of samples for subsequent more sensitive and time-consuming analyses.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Multiscale Modeling of Silicon Carbide Cladding for Nuclear Applications: Thermal Performance Modeling

The complex multiscale and anisotropic nature of silicon carbide (SiC) ceramic matrix composite (CMC) makes it difficult to accurately model its performance in nuclear applications. The existing models for nuclear grade composite SiC do not account for the microstructural features and how these features can affect the thermal and structural behavior of the cladding and its anisotropic properties. In addition to the microstructural features, the properties of individual constituents of the composites and fiber tow architecture determine the bulk properties. Models for determining the relationship between the individual constituents’ properties and the bulk properties of SiC composites for nuclear applications are absent, although empirical relationships exist in the literature. Here, a hierarchical multiscale modeling approach was presented to address this challenge. This modular approach addressed this difficulty by dividing the various aspects of the composite material into separate models at different length scales, with the evaluated property from the lower-length-scale model serving as an input to the higher-length-scale model. The multiscale model considered the properties of various individual constituents of the composite material (fiber, matrix, and interphase), the porosity in the matrix, the fiber volume fraction, the composite architecture, the tow thickness, etc. By considering inhomogeneous and anisotropic contributions intrinsically, our bottom-up multiscale modeling strategy is naturally physics-informed, bridging constitutive law from micromechanics to meso-mechanics and structural mechanics. The effects that these various physical attributes and thermo-physical properties have on the composite’s bulk thermal properties were easily evaluated and demonstrated through the various analyses presented herein. Since silicon carbide fiber-reinforced SiC CMCs are also promising thermal–structural materials with a broad range of high-end technology applications beyond nuclear applications, we envision that the multiscale modeling method we present here may prove helpful in future efforts to develop and construct reinforced CMCs and other advanced composite nuclear materials, such as MAX phase materials, that can service under harsh environments of ultrahigh temperatures, oxidation, corrosion, and/or irradiation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Bayesian Optimization Framework for Imperfect Data or Models

Conventional Bayesian optimization methods implicitly assume that the data and model being optimized are “perfect.” This assumption leads to inaccurate posterior probability distribution functions (PDFs) when applied to “imperfect” data or models. The new Bayesian optimization framework presented in this report provides a way to parameterize the effect of imperfections usually encountered in a prior PDF of generalized data or a model on the posterior PDF. The effects of imperfections are parameterized by a set of constraints imposed on the posterior expectation values of deviations between the data and the model and on their covariance matrix elements. A particular set of values for these constraints conveys an evaluator’s best estimate of the effect of imperfections on the corresponding posterior expectation values. When a prior PDF of generalized data is assumed to be normal, an expression for a posterior PDF satisfying an arbitrary set of constraints is derived analytically for linear models. An analogous iterative algorithm is given for nonlinear models. The corresponding posterior PDF should be used to estimate any posterior expectation values in the presence of imperfections parameterized by that set of constraints. A posterior PDF of a conventional Bayesian optimization method is recovered analytically when all evaluator-specified constraints are set to zero (i.e., in the absence of any imperfections). The analytical expressions derived in this report for normal PDFs and linear models were verified numerically by a Metropolis–Hastings Monte Carlo method. The methods presented herein could be applied to any kind of data or models, including differential cross-section data or integral benchmark experiments.

97 MATHEMATICS AND COMPUTING↗

Structure of Low-spin States in \(^{45}\)Sc Studied via Coulomb Excitation

The electromagnetic structure of 45 Sc at low excitation energy was investigated via low-energy Coulomb excitation at the Heavy Ion Laboratory (HIL) of the University of Warsaw and at the Inter-University Accelerator Centre (IUAC) in New Delhi. A set of reduced E2, E3, and M1 matrix elements was extracted from the collected data using the GOSIA code. The reduced transition probability B(E2; 11/2 – → 7/2 – ) has been determined, allowing us to deduce the lifetime of the 11/2 – state at 1237 keV. In addition, the upper limit on the reduced transition probability B(E3; 7/2 – → 5/2 + ) has been determined for the first time. Finally, new large-scale shellmodel and beyond-mean-field calculations were performed to interpret the structure of this nucleus.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Harvesting 88 Zr from heavy-ion beam irradiated tungsten at the National Superconducting Cyclotron Laboratory

Tungsten is a commonly used material at many heavy-ion beam facilities, and it often becomes activated due to interactions with a beam. Many of the activation products are useful in basic and applied sciences if they can be recovered efficiently. In order to develop the radiochemistry for harvesting group (IV) elements from irradiated tungsten, a heavy-ion beam containing 88 Zr was embedded into a stack of tungsten foils at the National Superconducting Cyclotron Laboratory and a separation methodology was devised to recover the 88 Zr. The foils were dissolved in 30% hydrogen peroxide, and the 88 Zr was chemically purified from the tungsten matrix and from other co-implanted radionuclides (such as 85 Sr and 88 Y) using strong cation-exchange (AG MP-50) chromatographic resin in sulfuric acid media. The procedure provided 88 Zr in approximately 60 mL 0.5 M sulfuric acid with no detectable radio-impurities. The overall recovery yield for 88 Zr was (92.3 ± 1.2)%. In conclusion, this proof-of-concept experiment has facilitated the development of methodologies to harvest from tungsten and tungsten-alloy parts that are regularly irradiated at heavy-ion beam facilities.

07 ISOTOPE AND RADIATION SOURCES↗

JIMWLK on a quantum computer

We propose a method for solving the Jalilian-Marian-Iancu-McLerran-Weigert-Leonidov-Kovner (JIMWLK) evolution equation on quantum computers. Our approach exploits the reformulation of the JIMWLK equation as a Lindblad master equation governing the rapidity evolution of the hadronic density matrix, as established in prior work. To render the problem tractable for quantum simulation, we introduce several approximations: the two-dimensional transverse plane is reduced to a one-dimensional radial lattice by assuming azimuthal symmetry of the jump operators; the gauge group is restricted to SU(2); and the infinite Wilson lines of the JIMWLK equation are replaced by finite Wilson links along the light-cone direction. The resulting bosonic Hilbert space is truncated using the electric field basis familiar from Hamiltonian lattice gauge theory, with states restricted to angular momenta 𝑗 ≤ 𝑗 max . We derive the matrix elements of the JIMWLK Lindblad jump operators in this basis. As a benchmark, we demonstrate rapid convergence of the fundamental dipole expectation value with 𝑗 max for both pure and mixed Gaussian initial density matrices. For the simplest truncation, 𝑗 max =1/2, we implement the Lindblad evolution using a quantum simulation algorithm verified with the Qiskit statevector simulator by decomposing the non-unitary evolution operator into a linear combination of unitaries. This work establishes a concrete pathway toward quantum simulation of high-energy quantum chromodynamics evolution equations, with direct relevance to the physics program of the Electron-Ion Collider.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Complete set of bound negative-parity states in the neutron-rich nucleus 18 N

High-resolution γ -ray spectroscopy of N 18 is performed with the Advanced GAmma Tracking Array, following deep-inelastic processes induced by an O 18 beam on a Ta 181 target. Six states are newly identified, which together with the three known excitations exhaust all negative-parity excited states expected in N 18 below the neutron threshold. Spin and parities are proposed for all located states on the basis of decay branchings and comparison with large-scale shell-model calculations performed in the p-sd space, with the YSOX interaction. Of particular interest is the location of the 0 1 - and 1 2 - excitations, which provide strong constrains for cross-shell p-sd matrix elements based on realistic interactions and help to simultaneously reproduce the ground and first-excited states in N 16 and N 18 , for the first time. Overall, understanding the N 18 structure may also have significant impact on neutron-capture cross-section calculations in r-process modeling including light neutron-rich nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Coulomb excitation of 80,82 Kr and a change in structure approaching N = Z = 40

Nuclei approaching N = Z = 40 are known to exhibit strongly deformed structures and are thought to be candidates for shape coexistence. In the krypton isotopes, 80,82 Kr are poorly characterized, preventing an understanding of evolving deformation approaching N = 40. The present work aims to determine electric quadrupole transition strengths and quadrupole moments of 80,82 Kr in order to better characterize their deformation. Sub-barrier Coulomb excitation was employed, impinging the isotopes of krypton on 196 Pt and 208 Pb targets. Utilizing a semiclassical description of the safe Coulomb-excitation process E 2 matrix elements could then be determined. Eleven new or improved matrix elements are determined in 80 Kr and seven in 82 Kr . The new B(E2;0$^+_1$ → 2$^+_1$ ) value in 82 Kr disagrees with the evaluated value by 3σ , which can be explained in terms of deficiencies in a previous Coulomb-excitation analysis. Conclusions: Comparison of measured Q s ( 2$^+_1$ ) and B(E2; 0$^+_1$ → 2$^+_1$ ) values indicates that neutron-deficient ( N ≤ 42 ) isotopes of krypton are closer to axial deformation than other isotopic chains in the mass region. A continuation of this trend to higher Z may result in Sr and Zr isotopes exhibiting near-axial prolate deformation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Factorization at subleading power and endpoint divergences in h → γγ decay. Part II. Renormalization and scale evolution

Building on the recent derivation of a bare factorization theorem for the b-quark induced contribution to the h → γγ decay amplitude based on soft-collinear effective theory, we derive the first renormalized factorization theorem for a process described at subleading power in scale ratios, where λ = m b /M h « 1 in our case. We prove two refactorization conditions for a matching coefficient and an operator matrix element in the endpoint region, where they exhibit singularities giving rise to divergent convolution integrals. The refactorization conditions ensure that the dependence of the decay amplitude on the rapidity regulator, which regularizes the endpoint singularities, cancels out to all orders of perturbation theory. We establish the renormalized form of the factorization formula, proving that extra contributions arising from the fact that “endpoint regularization” does not commute with renormalization can be absorbed, to all orders, by a redefinition of one of the matching coefficients. We derive the renormalization-group evolution equation satisfied by all quantities in the factorization formula and use them to predict the large logarithms of order $αα^2_sL^k$ in the three-loop decay amplitude, where $L=ln(-M^2_h/m^2_b)$ and k = 6, 5, 4, 3. We find perfect agreement with existing numerical results for the amplitude and analytical results for the three-loop contributions involving a massless quark loop. On the other hand, we disagree with the results of previous attempts to predict the series of subleading logarithms $~αα^n_sL^{2n+1}$

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Mach Conference, 2022.

There has been considerable interest in chemical vapor infiltration (CVI) manufactured silicon carbide fiber and silicon carbide matrix (SiC/SiC) composite tubes due to their superior mechanical properties. The SiC/SiC tubes are known to remain stable even after prolonged exposure to radiation. However, there are many variables in the manufacturing of SiC/SiC composite tubes such as elastic constants of constituent SiC fiber and SiC matrix, braiding angle, porosity, and others whose values vary over a range. The intention is to identify the distribution of elastic constants of SiC/SiC tubes which would help model the performance of the novel cladding material. To do so, a sensitivity analysis will be carried out to determine the dominant variables that influence the elastic constants of SiC/SiC tubes. A chosen number of combinations of dominant variables will be sampled through the Latin hypercube sampling (LHS) method. The elastic constants of composite tube namely Young's moduli in the circumferential and longitudinal direction, Poisson's ratio, and shear modulus will be calculated at sampled points through finite element (FE) analysis. A polynomial response surface (PRS) will be built for each of the elastic constants to be used as a surrogate to FE analysis. Once the predictive accuracy of PRS is verified at the validation set of data points, the PRS will be invoked 105 times in Monte Carlo simulations (MCS). The uncertainty can be quantified by calculating the coefficient of variation (CV) based on MCS-generated data which will be presented at the conference.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Measurement of 𝐽∕𝜓and 𝜓(2𝑆)production in 𝑝 +𝑝and 𝑝 +𝑑interactions at 120 GeV

We report the p + p and p + d differential cross sections measured in the SeaQuest experiment for J/psi and psi(2S) production at 120GeV beam energy covering the forward x-Feynman (x(F)) range of 0.5 < x(F) < 0.9. The measured cross sections are in good agreement with theoretical calculations based on the nonrelativistic QCD (NRQCD) using the long-distance matrix elements deduced from a recent global analysis of proton- and pion-induced charmonium production data. The sigma(psi(2S))/sigma(J/psi) cross section ratios are found to increase as x(F) increases, indicating that the q (q) over bar. annihilation process has larger contributions in the psi(2S) production than the J/psi production. The sigma(pd)/2 sigma(pp) cross section ratios are observed to be significantly different for the DrellYan process and J/psi production, reflecting their different production mechanisms. We find that the sigma(pd)/2 sigma(pp) ratios for J/psi production at the forward x(F) region are sensitive to the (d) over bar/(u) over bar flavor asymmetry of the proton sea, analogous to the Drell-Yan process. The transverse momentum (p(T)) distributions for J/psi and psi(2S) production are also presented and compared with data collected at higher center-of-mass energies.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Engineering defect energy landscape of CoCrFeNi high-entropy alloys by the introduction of additional dopants

The concept of high-entropy alloys (HEAs) focusing on tuning the overall chemical complexity represents a novel alloy design strategy. In contrast, alloying of a metallic matrix with minor doping elements with limited and localized tunability has been a common practice to improve material performance. Combining the idea of globally engineering defect energy landscape in HEAs and the localized doping strategy in dilute alloys, in this work, we explore doping effects of minor elements in a HEA matrix to further enhance the overall and localized chemical tunability, aiming to improve its irradiation resistance. Specifically, we study the influence of minor Al, Cu, Ti, and Pd substitutional doping elements on defect energetics in a CoCrFeNi model HEA based on density-functional theory (DFT) calculations. The DFT results indicate that the formation and migration energies of vacancies can be strongly influenced when a dopant is introduced at the first nearest neighbor shells around a vacancy. On the other hand, interstitial energetics are only slightly affected. Among the four elements, Ti and Pd generally decrease vacancy formation energies and increase vacancy migration energies more significantly than Al and Cu. The doping effects become more pronounced when the concentration of the substitutional dopants increases. Based on the energy distributions obtained from DFT, we build a kinetic Monte Carlo (kMC) model to assess the impact of dopants on vacancy-mediated diffusivity in the doped HEAs. Our results suggest that Ti and Pd can lower the tracer diffusivity in the considered HEAs and act as trapping sites, whereas Cu may enhance the atomic transport. This work indicates that substitutional doping in HEAs is an effective strategy in metallurgy to further tune the defect and transport properties of complex alloys.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Triaxiality and the nature of low-energy excitations in Ge 76

The deformation properties of the low-lying states in 76 Ge have been investigated following a safe energy Coulomb excitation measurement with the GRETINA tracking array and CHICO2 heavy-ion counter at the ATLAS accelerator facility at Argonne National Laboratory. A comprehensive set of transition and static E2 matrix elements were extracted from the measured differential Coulomb cross-sections, and compared with results of configuration interaction shell-model calculations and computations carried out within the framework of the generalized triaxial rotor model. The remarkable agreement between the calculated and experimental data supports a near-maximum triaxial deformation for the ground state of 76 Ge. Additionally, the degree of softness of the asymmetry in 76 Ge and 76 Se was investigated using rotational invariants generated from configuration interaction shell-model wave functions computed with the jj44b and JUN45 effective interactions. The resulting invariants are shown to be consistent with a stiff triaxial deformation in 76 Ge and a predominantly soft triaxial potential for 76 Se, in agreement with the conclusions of recent works by this collaboration.

59 ≤ A ≤ 89↗