Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Models & methods for nuclear reactions”

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 253 records · Page 14

High-Precision Spectroscopy of 20 O Benchmarking Ab Initio Calculations in Light Nuclei

The excited states of unstable 20 O were investigated via γ-ray spectroscopy following the 19 O(d, p) ⁢20 O reaction at 8 A⁢ MeV. By exploiting the Doppler shift attenuation method, the lifetimes of the $2^+_2$ and $3^+_1$ states were firmly established. From the γ-ray branching and E⁢2/M⁢1 mixing ratios for transitions deexciting the $2^+_2$ and $3^+_1$ states, the B⁡(E⁢2) and B⁡(M⁢1) were determined. Various chiral effective field theory Hamiltonians, describing the nuclear properties beyond ground states, along with a standard USDB interaction, were compared with the experimentally obtained data. Such a comparison for a large set of γ-ray transition probabilities with the valence space in medium similarity renormalization group ab initio calculations was performed for the first time in a nucleus far from stability. It was shown that the ab initio approaches using chiral effective field theory forces are challenged by detailed high-precision spectroscopic properties of nuclei. Finally, the reduced transition probabilities were found to be a very constraining test of the performance of the ab initio models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Artificial-intelligence-assisted analysis of 28 Si * → 7⁢𝛼 breakup data

Mid-weight 𝛼-conjugate nuclei are predicted to possess exotic toroid like resonances with high angular momenta. The search for these states in 28 Si* is the main point of two published experimental investigations of the peripheral 28 Si + 12 C reaction by Cao and collaborators and by Hannaman and collaborators. In this work, we develop a novel artificial intelligence (AI) based machine learning method utilizing the Gaussian Mixture Model (GMM) to analyze available experimental and theoretical data. Here, we additionally study the reaction with the Hybrid 𝛼-Cluster (H⁡𝛼⁢C) model. In all the examined data, our results suggest the presence of underlying structure which is close to that predicted for toroidal states.

Breakup reactions↗

Active species in chloroaluminate ionic liquids catalyzing low-temperature polyolefin deconstruction

Abstract Chloroaluminate ionic liquids selectively transform (waste) polyolefins into gasoline-range alkanes through tandem cracking-alkylation at temperatures below 100 °C. Further improvement of this process necessitates a deep understanding of the nature of the catalytically active species and the correlated performance in the catalyzing critical reactions for the tandem polyolefin deconstruction with isoalkanes at low temperatures. Here, we address this requirement by determining the nuclearity of the chloroaluminate ions and their interactions with reaction intermediates, combining in situ 27 Al magic-angle spinning nuclear magnetic resonance spectroscopy, in situ Raman spectroscopy, Al K-edge X-ray absorption near edge structure spectroscopy, and catalytic activity measurement. Cracking and alkylation are facilitated by carbenium ions initiated by AlCl 3 - tert -butyl chloride (TBC) adducts, which are formed by the dissociation of Al 2 Cl 7 − in the presence of TBC. The carbenium ions activate the alkane polymer strands and advance the alkylation cycle through multiple hydride transfer reactions. In situ 1 H NMR and operando infrared spectroscopy demonstrate that the cracking and alkylation processes occur synchronously; alkenes formed during cracking are rapidly incorporated into the carbenium ion-mediated alkylation cycle. The conclusions are further supported by ab initio molecular dynamics simulations coupled with an enhanced sampling method, and model experiments using n-hexadecane as a feed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Reduced-Order Modeling of Multigroup Neutron Cross Sections for High-Temperature Gas-cooled Reactors

Abstract – Deterministic neutronics calculations rely on multigroup neutron cross section libraries, which consist of databases of tabulated values, used to calculate the neutron cross sections through multivariate linear interpolation. However, interpolation of the multidimensional cross section data becomes memory inefficient and time consuming as the number of tabulations increases, significantly slowing down the neutronics calculation, especially in the case of microscopic cross section libraries where every isotope (on the order of hundreds) has its own set of specific reactions and cross sections. In order to address this challenge, this work constructs efficient and robust reduced-order models (ROMs) of the multi-group cross sections to support the Griffin simulation of high-temperature gas-cooled reactors (HTGRs). The first part of the study investigates the linearity of the multigroup cross section data across isotopes, reaction types, and energy groups on pre-generated datasets for the purpose of dimensionality reduction. Secondly, a down-selection of ROM techniques is presented on representative classical machine learning (ML) techniques, including variants of linear regression, kernel-based methods, tree-based algorithms, and artificial neural networks. The selection criteria jointly consider the memory efficiency, predictive accuracy, prediction speed, and scalability in comparison to the multidimensional interpolation. Among all the ML techniques, deep neural networks (DNNs) have proven to be the best selection with sufficient accuracy, high robustness, good memory efficiency, great scalability, and superior flexibility. DNNs have been trained for all isotopes in this work and systematic Griffin testing is ongoing to ensure the feasibility of this ROM technique for predicting cross section and reducing memory requirements without a significant sacrifice in computational performance.

42 - ENGINEERING↗

Nucleon Structure Studies at Jefferson Lab and the Electron Ion Collider

The research programs of Thomas Jefferson Laboratory (JLab) and the future Electron- Ion Collider (EIC) focus on one of the main goals of strong interaction studies : understanding the structure of nucleons in terms of the quarks and gluons composing them. Their structure is encoded in functions such as Generalized Parton Distributions (GPDs), which describe how quarks and gluons’ transverse position and longitudinal momentum are distributed inside nucleons. GPDs allow to obtain three-dimensional pictures of nucleons and to understand some of their fundamental properties, such as their internal pressure or the emergence of their spin from the dynamics of the partons composing them. At JLab and the EIC, electron beams are used to probe nucleons. Measurement of reactions such as Deeply Virtual Compton Scattering (DVCS) allows access to GPDs. The first longitudinally polarized-target experiment of the CLAS12 program at Jlab took place in 2022-2023. Combining polarized electron beams and nucleon targets, this experiment offers unique access to observables that allow the measurement of different types of GPDs. In particular, the DVCS beam- and target-spin asymmetries for protons and neutrons in deuterium will be measured for the first time. They give access to kinds of GPDs that are still poorly known, and the comparison between proton and neutron data will allow the extraction of the flavor dependence of the structure of nucleons. Specific analysis methods have been implemented to work with a polarized nuclear target and are presented in this thesis. These methods allow to obtain preliminary results for the asymmetries, waiting for the complete statistics to be available. In the long term, the experimental program for the EIC has been established with a strong emphasis on the measurement of the structure of nucleons at high energy. Measurements of reactions such as DVCS impose strict requirements on the electromagnetic calorimeter that will allow to measure the energy of the scattered electrons and photons. This calorimeter, which is under development, will be based on scintillating crystals read by Silicon Photomultipliers (SiPMs). A new type of glass-based scintillating material was tested, evaluating the possibilities to meet the technical requirements concerning their light yield and resistance to radiation damage in particular. Several models of SiPMs have been characterized, demonstrating they can operate over the vast energy range necessary to address the physics case at the EIC and providing guidelines for developing their readout electronics.

Pilleux, Noemie↗

Modeling quasielastic interactions of monoenergetic kaon decay-at-rest neutrinos

Monoenergetic muon neutrinos at 236 MeV are readily produced in intense medium-energy proton facilities (≳2–3GeV) when a positive kaon decays at rest (KDAR; K + →μ + ν μ ). These neutrinos provide a unique opportunity to both study the neutrino interaction and probe the nucleus with a monoenergetic weak-interaction-only tool. We present cross-section calculations for quasielastic scattering of these 236-MeV neutrinos off 12 C and 40 Ar, paying special attention to low-energy aspects of the scattering process. Our model takes the description of the nucleus in a mean-field approach as the starting point, where we solve Hartree-Fock equations using a Skyrme type nucleon-nucleon interaction. Thereby, we introduce long-range nuclear correlations by means of a continuum random-phase approximation (CRPA) framework where we solve the CRPA equations using a Green's function method. The model successfully describes (e,e') data on 12 C and 40 Ca in the kinematic region that overlaps with the KDAR ν μ phase space. In addition to these results, we present future prospects for precision KDAR cross-section measurements and applications of our calculations in current and future experiments that will utilize these neutrinos.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Vibrational properties of heme-nitrosoalkane complexes in comparison with those of their HNO analogs, and reactivity studies towards nitric oxide and Lewis acids

C-Nitroso compounds (RNO, R = alkyl and aryl) are byproducts of drug metabolism and bind to heme proteins, and their heme-RNO adducts are isoelectronic to ferrous nitroxyl (NO-/HNO) complexes. Importantly, heme-HNO compounds are key intermediates in the reduction of NO to N 2 O and nitrite to ammonium in the nitrogen cycle. Ferrous heme-RNO complexes act as stable analogs of these species, potentially allowing for the investigation of the vibrational and electronic properties of unstable heme-HNO intermediates. In this paper, a series of six-coordinate ferrous heme-RNO complexes (where R = iPr and Ph) were prepared using the TPP 2- and 3,5-Me-BAFP 2- co-ligands, and tetrahydrofuran, pyridine, and 1-methylimidazole as the axial ligands (bound trans to RNO). These complexes were characterized using different spectroscopic methods and X-ray crystallography. The complex [Fe(TPP)(THF)(iPrNO)] was further utilized for nuclear resonance vibrational spectroscopy (NRVS), allowing for the detailed assignment of the Fe–N(R)O vibrations of a heme-RNO complex for the first time. The vibrational properties of these species were then correlated with those of their HNO analogs, using DFT calculations. Our studies support previous findings that RNO ligands in ferrous heme complexes do not elicit a significant trans effect. In addition, the complexes are air-stable, and do not show any reactivity of their RNO ligands towards NO. So although ferrous heme-RNO complexes are suitable structural and electronic models for their HNO analogs, they are unsuitable to model the reactivity of heme-HNO complexes. We further investigated the reaction of our heme-RNO complexes with different Lewis acids. Here, [Fe(TPP)(THF)(iPrNO)] was found to be unreactive towards Lewis acids. In contrast, [Fe(3,5-Me-BAFP)(iPrNO) 2 ] is reactive towards all of the Lewis acids investigated here, but in most cases the iron center is simply oxidized, resulting in the loss of the iPrNO ligand. In the case of the Lewis acid B 2 (pin) 2 , the reduced product [Fe(3,5-Me-BAFP)(iPrNH 2 )(iPrNO)] was identified by X-ray crystallography.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Statistical analysis for the neutrinoless double-β-decay matrix element of 48 Ca

Neutrinoless double-β-decay (0⁢vββ) nuclear matrix elements (NME) are the object of many theoretical calculation methods, and are very important for analysis and guidance of a large number of experimental efforts. However, there are large discrepancies between the NME values provided by different methods. Here, in this paper, we propose a statistical analysis of the 48 Ca 0v⁢ββ NME using the interacting shell model, emphasizing the range of the NME probable values and their correlations with observables that can be obtained from the existing nuclear data. Based on this statistical analysis with three independent effective Hamiltonians, we propose a common probability distribution function for the 0⁢vββ NME, which has a range of (0.45–0.95) at 90% confidence level, and a mean value of 0.68.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Comparison of model predictions with measured proton-induced production of nickel and iridium isotopes

Calculations were performed to support method development for simultaneous production of Ni and Ir isotopes. Here, the work scope included development of the physical target, irradiation configuration, and post-irradiation radiochemical separations. The assumptions of predictive models previously developed were refined based upon the precise experimental configuration selected, including Os target material and a set of stacked targets for concurrent production of Ni and Ir isotopes. Model predictions of the reaction cross sections with EMPIRE 3.21 and implied isotopic yields are compared with proton beam irradiation measurements to refine model parameters and guide future experiments.

Isotope Production↗

Fatty acid capped, metal oxo clusters as the smallest conceivable nanocrystal prototypes

Metal oxo clusters of the type M 6 O 4 (OH) 4 (OOCR) 12 (M = Zr or Hf) are valuable building blocks for materials science. Here, we synthesize a series of zirconium and hafnium oxo clusters with ligands that are typically used to stabilize oxide nanocrystals (fatty acids with long and/or branched chains). The fatty acid capped oxo clusters have a high solubility but do not crystallize, precluding traditional purification and single-crystal XRD analysis. We thus develop alternative purification strategies and we use X-ray total scattering and Pair Distribution Function (PDF) analysis as our main method to elucidate the structure of the cluster core. We identify the correct structure from a series of possible clusters (Zr3, Zr4, Zr6, Zr12, Zr10, and Zr26). Excellent refinements are only obtained when the ligands are part of the structure model. Further evidence for the cluster composition is provided by nuclear magnetic resonance (NMR), infrared spectroscopy (FTIR), thermogravimetry analysis (TGA), and mass spectrometry (MS). We find that hydrogen bonded carboxylic acid is an intrinsic part of the oxo cluster. Using our analytical tools, we elucidate the conversion from a Zr6 monomer to a Zr12 dimer (and vice versa), induced by carboxylate ligand exchange. Finally, we compare the catalytic performance of Zr12-oleate clusters with oleate capped, 5.5 nm zirconium oxide nanocrystals in the esterification of oleic acid with ethanol. The oxo clusters present a five times higher reaction rate, due to their higher surface area. Since the oxo clusters are the lower limit of downscaling oxide nanocrystals, we present them as appealing catalytic materials, and as atomically precise model systems. In addition, the lessons learned regarding PDF analysis are applicable to other areas of cluster science as well, from semiconductor and metal clusters, to polyoxometalates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Microcanonical rates from ring-polymer molecular dynamics: Direct-shooting, stationary-phase, and maximum-entropy approaches

We address the calculation of microcanonical reaction rates for processes involving significant nuclear quantum effects using ring-polymer molecular dynamics (RPMD), both with and without electronically non-adiabatic transitions. After illustrating the shortcomings of the naive free-particle direct-shooting method, in which the temperature of the internal ring-polymer modes is set to the translational energy scale, we investigate alternative strategies based on the expression for the microcanonical rate in terms of the inverse Laplace transform of the thermal reaction rate. It is shown that simple application of the stationary-phase approximation (SPA) dramatically improves the performance of the microcanonical rates using RPMD, particularly in the low-energy region where tunneling dominates. Using the SPA as a Bayesian prior, numerically exact RPMD microcanonical rates are then obtained using maximum entropy inversion of the thermal reaction rates for both electronically adiabatic and non-adiabatic model systems. Finally, the direct-shooting method is revisited using the SPA-determined temperature for the internal ring-polymer modes, leading to a simple, direct-simulation method with improved accuracy in the tunneling regime. This work suggests a general strategy for the extraction of microcanonical dynamical quantities from RPMD (or other approximate thermal) simulations.

Tao, Xuecheng↗

The LLNL nuclear data infrastructure for the GNDS data format

The next generation of nuclear data infrastructure tools at the Livermore National Laboratory (LLNL) consists of pipeline of codes that read and process nuclear data from evaluated files saved in the new GNDS (Generalised Nuclear Data Structure) nuclear data format. The processing code FUDGE (For Updating Data and Generating Evaluations) is at the front-end of this pipeline as it reads and process the evaluated data for use in downstream transport codes. FUDGE is Python based with C and C++ extensions for computationally intensive tasks. As is the case for the evaluated data, the processed output is also saved in the GNDS format and the GIDI+ API is provided as the interface between the processed data and the transport codes. GIDI+ is a C++ based suite of codes and it includes GIDI (General Interaction Data Interface), a library for reading and writing GNDS data, and MCGIDI which is the cross section lookup, and reaction and product distribution sampling interface between Monte Carlo transport codes and the GNDS data. GIDI provides methods for easy access to the multi-group processed GNDS data and this is demonstrated through its implementation in ARDRA, the LLNL deterministic transport code. The evaluation and sampling methods in MCGIDI are available as both CPU and GPU methods which facilitates the use of MCGIDI in both traditional CPU-based as well as the next generation mixed model computational architectures. This is demonstrated through the GIDI+ implementation in MERCURY, the LLNL Monte Carlo transport code. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Evolution of collectivity in 118 Xe

A recoil-distance Doppler shift experiment has been performed using the Pd 102 ( F 19 , p 2 n ) reaction at a beam energy of 73 MeV to measure the lifetime of excited states in Xe 118 . Additionally, the differential decay-curve method using γ γ coincidences and a gating procedure that allows to extract the lifetime without feeding assumptions has been employed. The lifetimes obtained for the yrast states up to spin-parity 8 + are compared with interacting boson model calculations and Xe 118 can be classified as a transitional nucleus between the spherical and a deformed shape. Systematics of the B ( E 2 ) values for the 2 + → 0 + and 4 + → 2 + transitions in the isotopic chains of tin, tellurium and xenon are presented. It is proposed that a “critical point” exists at which the B ( E 2 ; 4 + → 2 + ) / B ( E 2 ; 2 + → 0 + ) ratio drops to unity for lower neutron numbers within the isotopic chain. The position of the “critical point” varies with proton number, i.e., it is presumed to be located at the same mass number A = 114 in the Sn, Te, and Xe isotopes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Exploring Resonance Structures in the Partial-Wave Analysis of ¿p0 Photoproduction at GlueX

This thesis studies what happens when a photon (¿) collides with a proton (p) and produces a neutral omega (¿) and pion (p0) pair, with a recoiling proton (p'), expressed as ¿p ¿ ¿p0p'. We study this and other reactions to better understand the strong nuclear force; one of the four fundamental forces that govern all the physics of the universe. This force is specifically responsible for the binding and decay of subatomic particles, such as the ones here. While we understand the ¿ and p0, what we are actually interested in is a short-lived unknown particle X that decays via X ¿ ¿p0. There are a multitude of possible particles X can be, and so our focus in this work is to find out what X is by determining its properties from the particles we measure. We do this via an intricate analysis procedure known as “partial-wave analysis”. By analogy, one can think of our particle detector as a buoy, and the particles we want to analyze (X) as pebbles hitting a pond. The waves created by the pebble will move our buoy, giving us information about the pebble that produced the wave. However, when multiple pebbles hit our pond, the waves overlap and interfere with each other. Our buoy only can measure the complicated interfering result of all the waves. To disentangle this, our partial-wave analysis works by modeling this interference pattern so that we may infer the properties of the pebbles (particles) we produced. In this thesis, we review the relevant experimental history in photoproduction and related production mechanisms, as well as the theoretical foundations that motivate our measurement. We describe the methods used to collect our data at the GlueX experiment stationed at the Jefferson Lab accelerator facility in Newport News, Virginia. We then detail the selections we apply to ensure our events are almost exclusively ¿p ¿ ¿p0p'. We cover the intensity model, how we verify its capabilities, and finally present our results together with systematic studies. Our primary result is the detection of a b1(1235) meson interfering with a wide JPC = 1-- vector state, measured via a mass-independent partial-wave analysis. It provides precise experimental results that can be used as input for theoretical models of the reaction, yielding conclusions about the procedures responsible for how our universe behaves at its most basic level.

Scheuer, Kevin [College of William and Mary, Willi↗

Propagating neutronic uncertainties for FFTF LOFWOS Test #13

The safety evaluation conducted for licensing a Sodium-cooled Fast Reactor (SFR) may require modeling transients with best-estimate calculation tools that must first be validated against real-world measurements. To provide the community with a valuable benchmarking opportunity for validating SFR analysis tools and methods, the International Atomic Energy Agency (IAEA) initiated a coordinated research project (CRP) in 2018 for the analysis of the Fast Flux Test Facility (FFTF) Loss of Flow Without Scram (LOFWOS) Test #13. The impact of nuclear data uncertainties on neutronics parameters was previously investigated based on the COMMARA-2.0 covariance matrix. Since the transient simulation results are very sensitive to certain reactivity coefficients, it was decided to employ rigorous uncertainty propagation methods to quantify the impact of nuclear data uncertainties on the best-estimate predication of FFTF LOFWOS Test #13. The DAKOTA code is used to propagate neutronic uncertainties through SAS4A/SASSYS-1 transient simulations, while taking into account spatial and reaction-wise correlations within these uncertainties. This study shows that the remaining discrepancies observed between the Argonne National Laboratory (ANL) best-estimate results and the experimental measurements can be partly explained by the uncertainty associated with Gas Expansion Module (GEM) worth, which contributes the majority of the overall nuclear data uncertainty on the output from the FFTF LOFWOS Test #13 transient simulation. This study also confirmed the importance of including spatial and reaction-wise correlations of nuclear data uncertainties on feedback coefficients in the uncertainty propagation to avoid under-estimating their impact during the transient simulations.

42 ENGINEERING↗

Accurate determination of production data of the non-standard positron emitter 86 Y via the 86 Sr(p,n)-reaction

In view of several significant discrepancies in the excitation function of the 86 Sr(p,n) 86g+xm Y reaction which is the method of choice for the production of the non-standard positron emitter 86 Y for theranostic application, we carried out a careful measurement of the cross sections of this reaction from its threshold up to 16.2 MeV at Forschungszentrum Jülich (FZJ) and from 14.3 to 24.5 MeV at LBNL. Thin samples of 96.4% enriched 86 SrCO 3 were prepared by sedimentation and, after irradiation with protons in a stacked-form, the induced radioactivity was measured by high-resolution γ -ray spectrometry. The projectile flux was determined by using the monitor reactions nat Cu(p,xn) 62,63,65 Zn and nat Ti(p,x) 48 V, and the calculated proton energy for each sample was verified by considering the ratios of two reaction products of different thresholds. Additionally, the experimental cross section data obtained agreed well with the results of a nuclear model calculation based on the code TALYS. From the cross section data, the integral yield of 86 Y was calculated. Over the optimum production energy range E p = 14 → 7 MeV the yield of 86 Y amounts to 291 MBq/μA for 1 h irradiation time. This value is appreciably lower than the previous literature values calculated from measured and evaluated excitation functions. It is, however, more compatible with the experimental yields of 86 Y obtained in clinical scale production runs. The levels of the isotopic impurities 87m Y, 87g Y, and 88 Y were also estimated and found to be <2% in sum.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Preparation of excited states for nuclear dynamics on a quantum computer

We study two different methods to prepare excited states on a quantum computer, a key initial step to study dynamics within linear response theory. The first method uses unitary evolution for a short time T = O(√1 - F) to approximate the action of an excitation operator Ô with fidelity F and success probability P ≈ 1 – F. The second method probabilistically applies the excitation operator using the Linear Combination of Unitaries (LCU) algorithm. We benchmark these techniques on emulated and real quantum devices, using a toy model for thermal neutron-proton capture. Despite its larger memory footprint, the LCU-based method is efficient even on current generation noisy devices and can be implemented at a lower gate cost than a naive analysis would suggest. Here, these findings show that quantum techniques designed to achieve good asymptotic scaling on fault tolerant quantum devices might also provide practical benefits on devices with limited connectivity and gate fidelity.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Update on NNL TSL Evaluations and Validation [Slides]

12 new/revised TSL evaluations were contributed to ENDF/B-VIII.1. The thermal diffusion length $L$ is an integral property of a single material’s absorption and scattering cross sections (both differential and integral). No other neutron reactions or materials are involved. The MC21-calculated $L$ for water is consistent with the spread of experimental data and is sufficiently sensitive to different H-H 2 O TSL physics models to use the method as a TSL integral performance benchmark. Modern high-quality diffusion experiments at elevated $T$ would allow direct low-cost physics benchmarking of water TSLs when public elevated-$T$ critical benchmarks are limited. NNL is working with LLNL and RPI to develop experimental PNDA capability. An ICSBEP Fundamental Physics benchmark of the Nassar and Murphy PNDA experiment is being developed to provide an example of how to evaluate these type of benchmarks.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗