Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nuclear structure functions”

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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Disordered interfaces of alkaline aluminate salt hydrates provide glimpses of Al 3+ coordination changes

Hypothesis: The precipitation and dissolution of aluminum-bearing mineral phases in aqueous systems often proceed via changes in both aluminum coordination number and connectivity, complicating molecular-scale interpretation of the transformation mechanism. Here, the thermally induced transformation of crystalline sodium aluminum salt hydrate, a phase comprised of monomeric octahedrally coordinated aluminate which is of relevance to industrial aluminum processing, has been studied. Because intermediate aluminum coordination states during melting have not previously been detected, it is hypothesized that the transition to lower coordinated aluminum ions occurs within a highly disordered quasi-two-dimensional phase at the solid-solution interface. Experiments and simulations: In this work, in situ X-ray diffraction (XRD), Raman and 27 Al nuclear magnetic resonance (NMR) spectroscopy were used to monitor the melting transition of nonasodium aluminate hydrate (NSA, Na 9 [Al(OH) 6 ] 2 ·3(OH)·6H 2 O). A mechanistic interpretation was developed based on complementary classical molecular dynamics (CMD) simulations including enhanced sampling. A reactive forcefield was developed to bridge speciation in the solution and in the solid phase. Findings: In contrast to classical dissolution, aluminum coordination change proceeds through a dynamically stabilized ensemble of intermediate states in a disordered layer at the solid-solution interface. In both melting and dissolution of NSA, octahedral, monomeric aluminum transition through an intermediate of pentahedral coordination. The intermediate dehydroxylates to form tetrahedral aluminate (Al(OH)$_4^–$) in the liquid phase. This coordination change is concomitant with a breaking of the ionic aluminate-sodium ion linkages. The solution phase Al(OH)$_4^–$ ions subsequently polymerize into polynuclear aluminate ions. However, there are some differences between bulk melting and interfacial dissolution, with the onset of the surface-controlled process occurring at a lower temperature (~30 °C) and the coordination change taking place more gradually as a function of temperature. This work to determine the local structure and dynamics of aluminum in the disordered layer provides a new basis to understand mechanisms controlling aluminum phase transformations in highly alkaline solutions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Detailed study of quark-hadron duality in spin structure functions of the proton and neutron

Background: The response of hadrons, the bound states of the strong force (QCD), to external probes can be described in two different, complementary frameworks: as direct interactions with their fundamental constituents, quarks and gluons, or alternatively as elastic or inelastic coherent scattering that leaves the hadrons in their ground state or in one of their excited (resonance) states. The former picture emerges most clearly in hard processes with high momentum transfer, where the hadron response can be described by the perturbative expansion of QCD, while at lower energy and momentum transfers, the resonant excitations of the hadrons dominate the cross section. The overlap region between these two pictures, where both yield similar predictions, is referred to as quark-hadron duality and has been extensively studied in reactions involving unpolarized hadrons. Some limited information on this phenomenon also exists for polarized protons, deuterons, and 3 He nuclei, but not yet for neutrons. Purpose: In this paper, we present comprehensive and detailed results on the correspondence between the extrapolated deep inelastic structure function g 1 of both the proton and the neutron with the same quantity measured in the nucleon resonance region. Thanks to the fine binning and high precision of our data, and using a well-controlled perturbative QCD (pQCD) fit for the partonic prediction, we can make quantitative statements about the kinematic range of applicability of both local duality and global duality. Method: We use the most updated QCD global analysis results at high x from the Jefferson Lab Angular Momentum Collaboration to extrapolate the spin structure function g 1 into the nucleon resonance region and then integrate over various intervals in the scaling variable x. We compare the results with the large data set collected in the quark-hadron transition region by the CLAS Collaboration, including, for the first time, deconvoluted neutron data, integrated over the same intervals. We present this comparison as a function of the momentum transfer Q 2 . Results: We find that, depending on the integration interval and the minimum momentum transfer chosen, a clear transition to quark-hadron duality can be observed in both nucleon species. Furthermore, we show, for the first time, the approach to scaling behavior for g 1 measured in the resonance region at sufficiently high momentum transfer. Conclusions: Here, our results can be used to quantify the deviations from the applicability of pQCD for data taken at moderate energies and can help with extraction of quark distribution functions from such data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Dramatic Conformer-Dependent Reactivity of the Acetaldehyde Oxide Criegee Intermediate with Dimethylamine Via a 1,2-Insertion Mechanism

The reactivity of carbonyl oxides has previously been shown to exhibit strong conformer and substituent dependencies. Through a combination of synchrotron multiplexed photoionization mass spectrometry experiments (298 K, 4 Torr) and high-level theory (CCSD(T)-F12/cc-pVTZ-F12//B2PLYP-D3/cc-pVTZ with an added CCSDT(Q) correction), we explore the conformer dependence of the reaction of acetaldehyde oxide (CH 3 CHOO) with dimethyl amine (DMA). The experimental data supports the theoretically predicted 1,2-insertion mechanism and the formation of an amine-functionalized hydroperoxide reaction product. Tunable-VUV photoionization probing of anti- or anti- + syn-CH 3 CHOO reveals a strong conformer dependence of the title reaction. Here, the rate coefficient of DMA with anti-CH 3 CHOO is predicted to exceed that for the reaction with syn-CH 3 CHOO by a factor of ~34,000, which is attributed to submerged barrier (syn) vs. barrierless (anti) mechanisms for energetically downhill reactions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Classical and quantum computing of shear viscosity for ( 2 + 1 ) D SU(2) gauge theory

We perform a nonperturbative calculation of the shear viscosity for ( 2 + 1 )-dimensional SU(2) gauge theory by using the lattice Hamiltonian formulation. The retarded Green’s function of the stress-energy tensor is calculated from real time evolution via exact diagonalization of the lattice Hamiltonian with a local Hilbert space truncation, and the shear viscosity is obtained via the Kubo formula. When taking the continuum limit, we account for the renormalization group flow of the coupling but no additional operator renormalization. We find the ratio of the shear viscosity and the entropy density η s is consistent with a well-known holographic result 1 4 π at several temperatures on a 4 × 4 honeycomb lattice with the local electric representation truncated at j max = 1 2 . We also find the ratio of the spectral function and frequency ρ x y ( ω ) ω exhibits a peak structure when the frequency is small. Both the exact diagonalization method and simple matrix product state classical simulation method beyond j max = 1 2 on bigger lattices require exponentially growing resources. So we develop a quantum computing method to calculate the retarded Green’s function and analyze various systematics of the calculation including j max truncation and finite size effects, Trotter errors and the thermal state preparation efficiency. Our thermal state preparation method still requires resources that grow exponentially with the lattice size, but with a very small prefactor at high temperature. We test our quantum circuit on both the Quantinuum emulator and the IBM simulator for a small lattice and obtain results consistent with the classical computing ones. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Superradiance in alpha clustered mirror nuclei

Resonances in unstable quantum systems are radiating states that despite decaying overall normalization have a well-defined structure which is being balanced by outgoing radiation. Such an interplay between outgoing wave and internal quantum many-body dynamics leads to several unique effects. One of those is known as superradiance, or alignment, where due to decay or virtual coupling to the continuum the states undergo restructuring so that their wave functions align towards the decay channels thus facilitating the decay. This effect is well understood theoretically and is closely related to the fundamental properties of reaction physics. Direct observation of superradiance in open quantum many-body systems is difficult because it is hard to find identical complex quantum systems that are different only in their coupling to the continuum of reaction states describing the decay. Here we report this phenomenon in alpha cluster decays of mirror nuclei 18 O and 18 Ne.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Dynamics of McMillan mappings II. axially symmetric map

Here, in this article, we investigate the transverse dynamics of a single particle in a model integrable accelerator lattice, based on a McMillan axially-symmetric electron lens. Although the McMillan e-lens has been considered as a device potentially capable of mitigating collective space charge forces, some of its fundamental properties have not been described yet. The main goal of our work is to close this gap and understand the limitations and potentials of this device. It is worth mentioning that the McMillan axially symmetric map provides the first-order approximations of dynamics for a general linear lattice plus an arbitrary thin lens with motion separable in polar coordinates. Therefore, advancements in its understanding should give us a better picture of more generic and not necessarily integrable round beams. In the first part of the article, we classify all possible regimes with stable trajectories and find the canonical action-angle variables. This provides an evaluation of the dynamical aperture, Poincaré rotation numbers as functions of amplitudes, and thus determines the spread in nonlinear tunes. Also, we provide a parameterization of invariant curves, allowing for the immediate determination of the map image forward and backward in time. The second part investigates the particle dynamics as a function of system parameters. We show that there are three fundamentally different configurations of the accelerator optics causing different regimes of nonlinear oscillations. Each regime is considered in great detail, including the limiting cases of large and small amplitudes. In addition, we analyze the dynamics in Cartesian coordinates and provide a description of observable variables and corresponding spectra.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Isospin symmetry holds beyond the proton drip line: First observation of excited states in 55Cu

Isospin symmetry, combined with shell-model description, has long been successful in reproducing nuclear structure and predicting mirror energy differences. However, these descriptions have been primarily limited to mirror nuclei with well-bound ground states. As one approaches the proton drip line, coupling to the continuum becomes important. Despite this, little is known about medium-mass nuclei at the drip line. Here we report the measurement of the excitation energies of several unbound states in a nucleus beyond the proton drip line, 55Cu, performed via γ-ray spectroscopy at the Radioactive Isotope Beam Factory (RIKEN). Several low-lying yrast states are observed to decay predominantly by γ-ray emission. In contrast, the state–predicted to be strongly populated in the one-proton knockout reaction–is notably absent, suggesting that it decays primarily via prompt proton emission. We show that even for such an unbound system, shell-model calculations based on a harmonic oscillator potential reproduce the observed mirror energy differences remarkably well in all the observed states. Based on this result, we propose that a general pattern might emerge where, in medium-mass systems at the proton drip line, the stronger Coulomb barrier together with the centrifugal barrier effectively confines the wave functions, so that the shell model is able to describe low-lying states, enabling accurate predictions of mirror energy differences between isobaric analogue states.

Pigliapoco, S [Dipartimento di Fisica e Astronomia↗
DOE OSTI2020

Hadronic structure in high-energy collisions

Here, Parton distribution functions (PDFs) describe the structure of hadrons as composed of quarks and gluons. They are needed to make predictions for short-distance processes in high-energy collisions and are determined by fitting to cross-section data. Definitions of the PDFs and their relations to high-energy cross sections are reviewed. The focus is on the PDFs in protons, but PDFs in nuclei are also discussed. The standard statistical treatment needed to fit the PDFs to data using the Hessian method is reviewed in some detail. Tests are discussed that critically examine whether the needed assumptions are indeed valid. Also presented are some ideas of what one can do in case tests indicate that the assumptions fail.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Unsupervised Resource Allocation with Graph Neural Networks

We present an approach for maximizing a global utility function by learning how to allocate resources in an unsupervised way. We expect interactions between allocation targets to be important and therefore propose to learn the reward structure for near-optimal allocation policies with a GNN. By relaxing the resource constraint, we can employ gradient-based optimization in contrast to more standard evolutionary algorithms. Our algorithm is motivated by a problem in modern astronomy, where one needs to select-based on limited initial information-among $10^9$ galaxies those whose detailed measurement will lead to optimal inference of the composition of the universe. Our technique presents a way of flexibly learning an allocation strategy by only requiring forward simulators for the physics of interest and the measurement process. We anticipate that our technique will also find applications in a range of resource allocation problems.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Global polarization of hyperons and spin alignment of vector mesons in quark matters

Relativistic heavy ion collider (RHIC) as a dedicated nuclear facility has made a few major discoveries in physics. This year marks the 30th year STAR Collaboration formation and the 23th year of STAR detector operation and data collection at RHIC. In the last two decades, STAR has collected many datasets, exhibiting scientific versatility and flexibility of the RHIC facility. The total dataset in the first year is less than 1 million good events, and currently there are about 1 billion events per dataset. The Global Hyperon Polarization was proposed in 2004. This immediately prompted the STAR Collaboration to search for this phenomenon from the early datasets. The null results were presented at Quark Matter Conference in Shanghai in 2006 and subsequently published. Although there were peripheral and continuous efforts in the following decade, no positive result has been observed experimentally. This situation changed in the following decade with the upgrade of high data rate and time-of-flight (TOF) detector and the progress of the Beam Energy Scan Phase I (BES-I). The experimental discoveries of the global polarization of hyperons in 2017 and the spin alignment of vector mesons in 2023 at RHIC-STAR confirm the theory which was established nearly twenty years ago. The theory and these measurements open the way to studying the properties of the hot and dense nuclear matter created in high-energy heavy ion collisions from a new degree of freedom, spin. We briefly review these discoveries from the proposals of theory to the experimental measurements, and summarize the related measurements at the existing facilities and the theoretical explanations to the original proposal. The basic understanding and the original proposal are still valid and fundamental, that is, the angular momentum of system can transform into a spin effect observable in experiment. However, it appears that in each case a new model is needed to explain the new experimental observation. We need a more basic theory to help us unify all these spin related phenomena. Over the past five years, STAR has successfully installed 3 new detectors and we have begun to see the physical analysis results from datasets with those new functions. What makes the STAR detector viable after 20 years of operation is its continuous evolution through successful upgrades, with new scientific programs added year by year. The next big thing is to forward upgrade a tracking system (3 layers of silicon strips and 4 layers of sTGC chambers) and a calorimetry system (electromagnetic and hadronic calorimeters). In addition to studying the spin structure of protons by using the polarized proton beams at RHIC, the upgrades also provide a unique ability to investigate the origin of Λ Global Polarization as a function of rapidity and rapidity (de-)correlations in Au+Au collisions.

Physics↗
DOE OSTI2021

Functionalized Hydroperoxide Formation from the Reaction of Methacrolein-Oxide, an Isoprene-Derived Criegee Intermediate, with Formic Acid: Experiment and Theory

Methacrolein oxide (MACR-oxide) is a four-carbon, resonance-stabilized Criegee intermediate produced from isoprene ozonolysis, yet its reactivity is not well understood. This study identifies the functionalized hydroperoxide species, 1-hydroperoxy-2-methylallyl formate (HPMAF), generated from the reaction of MACR-oxide with formic acid using multiplexed photoionization mass spectrometry (MPIMS, 298 K = 25 °C, 10 torr = 13.3 hPa). Electronic structure calculations indicate the reaction proceeds via an energetically favorable 1,4-addition mechanism. The formation of HPMAF is observed by the rapid appearance of a fragment ion at m/z 99, consistent with the proposed mechanism and characteristic loss of HO2 upon photoionization of functional hydroperoxides. The identification of HPMAF is confirmed by comparison of the appearance energy of the fragment ion with theoretical predictions of its photoionization threshold. The results are compared to analogous studies on the reaction of formic acid with methyl vinyl ketone oxide (MVK-oxide), the other four-carbon Criegee intermediate in isoprene ozonolysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Theoretical Prediction and Interpretation of 237 Np Mössbauer Isomer Shifts

Different theoretical approaches for the calculation of 237 Np Mössbauer isomer shifts are investigated. The traditional contact density route is compared to a previously proposed alternative approach that uses energy derivatives with respect to the nuclear radius. Both approaches yield similar results as long as suitable basis sets augmented with large exponents and relativistic methods are used. Density functional theory (DFT) calculations do not show a strong dependency of the 237 Np isomer shift on the chosen functional. Wavefunction calculations show that dynamic electron correlation can be important when covalent bonding influences the isomer shift. Effects from spin–orbit coupling are small. The isomer shifts of ionic solids and Np(III) organometallic complexes are largely governed by the oxidation state of Np. Isomer shifts of organometallic Np(IV) complexes are strongly affected by donation bonding. Furthermore, detailed analysis of the wavefunction results with different active spaces demonstrates that correlation among the outer core Np and occupied ligand frontier orbitals contributes significantly to isomer shifts of Np(IV) compounds.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

PDB ‐101: Educational resources supporting molecular explorations through biology and medicine

Abstract The Protein Data Bank (PDB) archive is a rich source of information in the form of atomic‐level three‐dimensional (3D) structures of biomolecules experimentally determined using macromolecular crystallography, nuclear magnetic resonance (NMR) spectroscopy, and electron microscopy (3DEM). Originally established in 1971 as a resource for protein crystallographers to freely exchange data, today PDB data drive research and education across scientific disciplines. In 2011, the online portal PDB‐101 was launched to support teachers, students, and the general public in PDB archive exploration ( pdb101.rcsb.org ). Maintained by the Research Collaboratory for Structural Bioinformatics PDB, PDB‐101 aims to help train the next generation of PDB users and to promote the overall importance of structural biology and protein science to nonexperts. Regularly published features include the highly popular Molecule of the Month series, 3D model activities, molecular animation videos, and educational curricula. Materials are organized into various categories (Health and Disease, Molecules of Life, Biotech and Nanotech, and Structures and Structure Determination) and searchable by keyword. A biennial health focus frames new resource creation and provides topics for annual video challenges for high school students. Web analytics document that PDB‐101 materials relating to fundamental topics (e.g., hemoglobin, catalase) are highly accessed year‐on‐year. In addition, PDB‐101 materials created in response to topical health matters (e.g., Zika, measles, coronavirus) are well received. PDB‐101 shows how learning about the diverse shapes and functions of PDB structures promotes understanding of all aspects of biology, from the central dogma of biology to health and disease to biological energy.

Zardecki, Christine↗
DOE OSTI2024

An extreme thermal cycling reliability test of ATLAS ITk Strips barrel modules

At the end of Run 3 of the Large Hadron Collider (LHC), the accelerator complex will be upgraded to the High-Luminosity LHC (HL-LHC) in order to increase the total amount of data provided to its experiments. To cope with the increased rates of data, radiation, and pileup, the ATLAS detector will undergo a substantial upgrade, including a replacement of the Inner Detector with a future Inner Tracker, called the ITk. The ITk will be composed of pixel and strip sub-detectors, where the strips portion will be composed of 17,888 silicon strip detector modules. During the HL-LHC running period, the ITk will be cooled and warmed a number of times from about -35°C to room temperature as part of the operational cycle, including warm-ups during yearly shutdowns. To ensure ITk Strips modules are functional after these expected temperature changes, and to ensure modules are mechanically robust, each module must undergo ten thermal cycles and pass a set of electrical and mechanical criteria before it is placed on a local support structure. This paper describes the thermal cycling Quality Control (QC) procedure, and results from the barrel pre-production phase (about 5% of the production volume). Additionally, in order to assess the headroom of the nominal QC procedure of 10 cycles and to ensure modules don't begin failing soon after, four representative ITk Strips barrel modules were thermally cycled 100 times — this study is also described.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Correlated Anion Disorder in Heteroanionic Cubic TiOF 2

Resolving anion configurations in heteroanionic materials is crucial for understanding and controlling their properties. For anion-disordered oxyfluorides, conventional Bragg diffraction cannot fully resolve the anionic structure, necessitating alternative structure determination methods. We have investigated the anionic structure of anion-disordered cubic (ReO 3 -type) TiOF 2 using X-ray pair distribution function (PDF), 19 F MAS NMR analysis, density functional theory (DFT), cluster expansion modeling, and genetic-algorithm structure prediction. Our computational data predict short-range anion ordering in TiOF 2 , characterized by predominant cis-[O 2 F 4 ] titanium coordination, resulting in correlated anion disorder at longer ranges. To validate our predictions, we generated partially disordered supercells using genetic-algorithm structure prediction and computed simulated X-ray PDF data and 19 F MAS NMR spectra, which we compared directly to experimental data. To construct our simulated 19 F NMR spectra, we derived new transformation functions for mapping calculated magnetic shieldings to predicted magnetic chemical shifts in titanium (oxy)fluorides, obtained by fitting DFT-calculated magnetic shieldings to previously published experimental chemical shift data for TiF 4 . We find good agreement between our simulated and experimental data, which supports our computationally predicted structural model and demonstrates the effectiveness of complementary experimental and computational techniques in resolving anionic structure in anion-disordered oxyfluorides. From additional DFT calculations, we predict that increasing anion disorder makes lithium intercalation more favorable by, on average, up to 2 eV, highlighting the significant effect of variations in short-range order on the intercalation properties of anion-disordered materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

TRISO Particle Isolation and Graphite Removal Using SRNL Vapor Digestion Technology

Tri-structural isotropic fueled reactors are planned to come online in the next decade, but there is not currently a widely agreed upon disposition pathway for this new fuel stream. For long term used nuclear fuel storage, there would be a significant benefit if the waste volume could be reduced or mitigated. Most of the volume from used TRISO fuel contains graphite. SRNL currently has an active patent which describes one pathway to digest the graphite. One benefit in the utilization of this pathway is that multiple cylindrical pieces can be stacked end-to-end in a reaction tube, then nitric acid and water vapor can be flowed through the tube, and the weight change of the individual cores can provide both a reaction profile and overall oxidant use efficiency. The reaction of nitric acid and graphite passes through a series of intermediate products – NO 2 , NO, and N 2 O – from reaction and decomposition to eventually form N2. The data shows that NO2 grows in with increasing temperature between 500-600°C and then decreases by way of either reaction with graphite or decomposition. Nitric oxide, a reaction and thermal decomposition product of NO 2 , exhibits a consistent decline as a function of temperature, which is consistent with the literature. Similarly, the concentrations of N 2 O and N 2 increase as a function of temperature as their reactions with graphite become more favorable.

Advanced Reactors↗
DOE OSTI2022

Determination of Molecular Structure and Dynamics of Molten Salts by Advanced Neutron and X-ray Scattering Measurements and Computer Modeling

The design and development of fully functional Molten-Salt Reactors (MSR) require detailed knowledge of the molten salt properties in order to understand and predict the salt’s behavior. Fundamental properties of interest include molecular structure, speciation, and dynamics (such as diffusion coefficients) of salt components and dissolved corrosion and fission products. Computer modeling is necessary to predict changes in physical and chemical properties due to irradiation, burning of dissolved fuel, and corrosion. The modeling requires experimental data, and advanced neutron and x-ray scattering and spectroscopy provide the most reliable and direct determination of the structure (Pair-Distribution Functions, PDF), and dynamics of ions in the melt. This project dealt with both fluoride and chloride salts. The PDFs have been measured by a combination of neutron and x-ray diffraction. We utilized the techniques of isotope substitutions, a very powerful tool available for neutron-scattering, to extract the details of the liquid structure. Although similar measurements have been done before, modern advanced neutron and x-ray-scattering techniques allow collecting the data at much higher resolution and in a wider range of temperatures. Importantly, we were among the first to study fluoride salts by neutron scattering. The importance of impurities and their effects on salt properties have become apparent recently and so new methods of salt purification were developed. We took advantage of these developments to produce reliable data, which have been used for computer simulations of both clean salts and those with added fission and corrosion products most relevant for MSRs. Ab initio molecular dynamics simulations have been performed to understand the multi-component liquid solution, in particular solubility of impurities and thermodynamic interactions in relation to the ionic-cluster structure of the fluid. We applied machine learning to regress from the simulation and experimental data in order to develop a fast-acting model that can handle molten salt with an arbitrary (≥ 10) number of chemical elements and be able to predict chemical potential as a function of composition and temperature. This project resulted in a number of experimental and computer-simulation publications, a patent application, and numerous conference presentations (American Physical Society, American Chemical Society, and The Electrochemical Society among others). Multiple students and postdocs participated and collaborated on aspects of this project. This project seeded new collaborations between MIT and other institutions, such as the University of Massachusetts Lowell, the University of Illinois Urbana-Champaign, the University of California Berkeley, and Oak Ridge and Los Alamos National Labs. As such, this project has had a broad and lasting impact beyond its original scientific scope.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Optimizing group structures using hierarchical division

Creating group structures with few groups that give low errors is a difficult problem in reactor analysis. In recent years, automated optimization techniques have been applied to this task. We continue this trend by applying the hierarchical division algorithm to generate optimized group structures that minimize a cost function. At each stage, the algorithm adds a single group boundary to an existing group structure, dividing one group into two to increase the resolution of the group structure. The location of the added boundary is the one that gives the lowest error over all possible new boundary locations. Our implementation requires a beginning group structure, a set of candidate new boundary locations, and a set of reference reaction rates. As a proof of concept, we used WIMS-69 as the initial group structure, XMAS-172 as the ending group structures, and a 344-group reference group structure. Testing on two simple, homogenized reactor problems, we found that hierarchical division was able to reduce the error by a factor of around 5 with an increase of only 15% in the number of groups. Because hierarchical division can get stuck in local minima, it often reaches a plateau in its error reduction capability as many groups are added. Nevertheless, we find hierarchical division has strong potential to make good group structures into great group structures at a modest increase in computational cost. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗