Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Bounds”

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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Ameloblastin binding to biomimetic models of cell membranes – A continuum of intrinsic disorder

A 37-residue amino acid sequence corresponding to the segment encoded by exon-5 of murine ameloblastin (Ambn), AB2 (Y67-Q103), has been implicated with membrane association, ameloblastin self-assembly, and amelogenin-binding. Here, our aim was to characterize, at the residue level, the structural behavior of AB2 bound to chemical mimics of biological membranes using NMR spectroscopy. To better define the structure of AB2 using NMR-based methods, recombinant 13 C- and 15 N-labelled AB2 (*AB2) was prepared and data collected free in solution and with deuterated dodecylphosphocholine (dPC) micelles, deuterated bicelles, and both small and large unilamellar vesicles. Amide chemical shift and intensity perturbations observed in 1 H- 15 N HSQC spectra of *AB2 in the presence of bicelles and dPC micelles suggest that a region of *AB2, S6-E36 (murine Ambn S68 – E98), associates with the membrane biomimetics. A CSI-3 analysis of the NMR chemical shift assignments for *AB2 free in solution and bound to dPC micelles indicated the peptide remains disordered except for the adoption of a short, 12-residue α-helix, F10-G21 (murine Ambn F72-G83). In dPC micelles, the NOE NMR data was void of patterns characteristic of long-lived helical structure indicating this helix was transient in nature. A continuum of intrinsic disorder in the membrane-bound state may be responsible for ameloblastin’s ability to dynamically interact with multiple partners at the same site during amelogenesis.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Scaled boundary cubature scheme for numerical integration over planar regions with affine and curved boundaries

This paper introduces the scaled boundary cubature (SBC) scheme for accurate and efficient integration of functions over polygons and two-dimensional regions bounded by parametric curves. Over two-dimensional domains, the SBC method reduces integration over a region bounded by m curves to integration over m regions (referred to as curved triangular regions), where each region is bounded by two line segments and a curve. With proper (counterclockwise) orientation of the boundary curves, the scheme is applicable to convex and nonconvex domains. Additionally, for star-convex domains, a tensor-product cubature rule with positive weights and integration points in the interior of the domain is obtained. If the integrand is homogeneous, we show that this new method reduces to the homogeneous numerical integration scheme; however, the SBC scheme is more versatile since it is equally applicable to both homogeneous and non-homogeneous functions. This paper also introduces several methods for smoothing integrands with point singularities and near-singularities. When these methods are used, highly efficient integration of weakly singular functions is realized. Finally, the SBC method is applied to a number of benchmark problems, which reveal its broad applicability and superior performance (in terms of time to generate a rule and accuracy per cubature point) when compared to existing methods for integration.

42 ENGINEERING↗
DOE OSTI2020

MADHAT: Model-Agnostic Dark Halo Analysis Tool

We present the Model-Agnostic Dark Halo Analysis Tool (MADHAT), a numerical tool which implements a Fermi-LAT data-driven, model-independent analysis of gamma-ray emission from dwarf satellite galaxies and dwarf galaxy candidates due to dark matter annihilation, dark matter decay, or other nonstandard or unknown astrophysics. This tool efficiently provides statistical upper bounds on the number of observed photons in excess of the number expected, based on empirical determinations of foregrounds and backgrounds, using a stacked analysis of any selected set of dwarf targets. It also calculates the resulting bounds on the properties of dark matter under any assumptions the user makes regarding dark sector particle physics or astrophysics. As an application, we determine new bounds on Sommerfeld-enhanced dark matter annihilation in a set of eight dwarfs. MADHAT v1.0 includes 58 dwarfs and dwarf candidate targets, and we discuss future planned developments. MADHAT is available and will be maintained at https://github.com/MADHATdm.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Three-dimensional Skyrme Hartree-Fock-Bogoliubov solver in coordinate-space representation

The coordinate-space representation of the Hartree-Fock-Bogoliubov theory is the method of choice to study weakly bound nuclei whose properties are affected by the quasiparticle continuum space. To describe such systems, we developed a three-dimensional Skyrme-Hartree-Fock-Bogoliubov solver HFBFFT based on the existing, highly optimized and parallelized Skyrme-Hartree-Fock code Sky3D. The code does not impose any self-consistent spatial symmetries such as mirror inversions or parity. The underlying equations are solved in HFBFFT directly in the canonical basis using the fast Fourier transform. To remedy the problems with pairing collapse, we implemented the soft energy cutoff and pairing annealing. The convergence of HFB solutions was improved by a sub-iteration method. The Hermiticity violation of differential operators brought by Fourier-transform-based differentiation has also been solved. Furthermore, the accuracy and performance of HFBFFT were tested by benchmarking it against other HFB codes, both spherical and deformed, for a set of nuclei, both well-bound and weakly-bound.

3D coordinate-space representation↗
DOE OSTI2020

Bio-reduction of ferrihydrite-montmorillonite-organic matter complexes: Effect of montmorillonite and fate of organic matter

Organic matter (OM) is often associated with Fe (hydr)oxides such as ferrihydrite (Fh) in soils and sediments, forming binary Fh-OM complexes. Microbial reduction of Fh results in destabilization of the complexes and mineral/OM transformation. However, little is known about the role of clay minerals in such processes, despite their common co-existence with Fh and OM in natural environments. Here Fh-OM complexes were synthesized in the presence of montmorillonite (SWy-2), forming ternary Fh-(SWy-2)-OM complexes. A metal-reducing bacterium Geobacter sulfurreducens was used to reduce Fh in the complexes under circumneutral pH and anoxic conditions with or without H 2 as extra electron donor. Various spectroscopy and mass spectrometry methods were used to monitor the progress of Fh bio-reduction and mineral/OM transformation. Results showed that G. sulfurreducens utilized mineral-bound OM as electron donor and/or carbon source to couple with Fh reduction. Relative to Fh-OM complex, Addition of SWy-2 to Fh-OM complex enhanced the bio-reduction extent of Fh by increasing the proportion of bioavailable OM that was weakly bound to SWy-2. However, its effect on the bio-reduction rate was variable. SWy-2 initially decreased the rate, because it spatially separated OM (electron donor) from Fh (electron acceptor). During later incubation, SWy-2 increased the reduction rate by sorbing biogenic Fe 2+ that would otherwise passivate the Fh and cell surfaces. Bio-reduction transformed mineral-bound OM to microbial products (e.g. necromass, extracellular polymeric substances), but organic compounds with aromatic structures, carboxyl groups and large molecular weight were more resistant to desorption and oxidation. The persistence of these compounds against bio-reduction induced transformation is likely due to their stronger binding with minerals and/or lower nominal oxidation states of carbon relative to other compounds. Our results provide new insights into the role of clay minerals in regulating biogeochemical cycling of solid-phase Fe and transformation of mineral-associated OM in anoxic soil environments.

58 GEOSCIENCES↗
DOE OSTI2024

The STAG code: A fully relativistic super transition array calculation using Green’s functions

Calculating opacities for a wide range of plasma conditions (i.e. temperature, density, element) requires detailed knowledge of the plasma configuration space and electronic structure. For plasmas composed of heavier elements, relativistic effects are important in both the electronic structure and the details of opacity spectra. We extend our previously described superconfiguration and super transition array capabilities (Gill et al., 2023) to include a fully relativistic formalism. Further, the use of hybrid bound-continuum supershells in our superconfigurations demonstrates the importance of a consistent treatment of bound and continuum electrons in dense plasma opacities, and we expand the discussion of these consequences to include issues associated with equation of state and electron correlations between bound and continuum electrons.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Losses of CO and CO 2 upon collision-activated dissociation of substituted 2-methoxyphenoxides after methyl radical loss

Some deprotonated, substituted 2-methoxyphenols fragment by methyl radical loss followed by competing losses of CO and CO 2 upon collisional activation in linear quadrupole ion trap mass spectrometers. These reactions were examined experimentally and computationally in order to determine their mechanisms. For deprotonated vanillin, the CO loss was found to involve ring contraction, with the highest free energy barrier of 58.0 kcal mol -1 (for the syn rotamer; here, syn refers to the relationship between the aldehyde oxygen atom and the oxygen atom of the methoxy group) for the entire process. The atoms lost in this fragmentation are the oxygen atom that was bound to the first eliminated methyl group and the aromatic carbon atom bound to this oxygen. Examination of carbon-13 labeled vanillin supports these assignments. Examination of several model compounds revealed that this reaction requires the presence of an electron-withdrawing substituent in the para-position relative to the phenol moiety. In contrast, the CO 2 loss from deprotonated vanillin occurs via ring opening followed by re-cyclization and then ring contraction, leading to the loss of CO 2 in a process wherein the highest free energy barrier is 84.5 kcal mol -1 (for the syn isomer). The atoms lost in this fragmentation are the carbon and oxygen atoms from the phenoxide group and the oxygen atom that was bound to the first eliminated methyl group, which is supported by examination of carbon-13 labeled vanillin. Despite the higher total free energy requirement, the CO 2 loss is competitive with CO loss, possibly due to favorable entropy and low activation energy (free energy barrier of 48.9 kcal mol -1 ) for the first reaction step (the analogous value for CO loss is 58.0 kcal mol -1 ). The extent of CO 2 loss is strongly affected by substituents – it either is not observed or is very slow for the other compounds studied here. For example, it is substantially less favorable than CO loss (highest free energy barrier 60.9 kcal mol -1 ) for deprotonated acetovanillone for which the first reaction step for CO 2 loss has a substantially greater free energy barrier (60.6 kcal mol-1) than for vanillin.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Spectroscopy and dynamics of the v = 1 levels derived from the OH(D) stretching modes in the I‾∙HDO complex using CW and time-resolved, resonant two-photon infrared excitation of the cryogenically cooled ions

The vibrational energy levels of the two isotopomers adopted by the I‾∙HDO ion-molecule complex occur such that the OH(D) stretching fundamentals span its dissociation energy, thus enabling a spectroscopic investigation of the dynamics displayed by a system prepared in the vicinity of the dissociation threshold. This regime is explored using infrared photoexcitation of the mass-selected complexes cooled in a cryogenic radiofrequency (Paul) ion trap. Survey spectra are obtained at modest resolution using two-color, IR-IR photodissociation with nanosecond lasers to establish the level structure and unimolecular decay dynamics of the v = 1 and 2 levels of the bound OH(D) oscillator. The v = 1 levels are prepared by fixed frequency excitation in the trap and the absorption spectra arising from this excited state are probed by photofragmentation of the complex with a second pulsed IR laser after a variable delay time (0 to 20 ms). The transitions to levels above the dissociation threshold for I‾ + HDO formation are observed to be sharp (~5 cm -1 FWHM). At low pressure, the bound OD (v = 1) population relaxes very slowly (~3 ms), consistent with resonant fluorescence in the IR. The collisional quenching rate constants of this level by the He buffer gas were estimated to be on the order 3 x 10 -12 cm 3 /s based on a crude Stern-Volmer analysis. Here, the rotational fine structure and linewidths (≲0.01 cm -1 FWHM) arising from transitions of the non-bonded OH stretch fundamental of the OD-bound isotopomer that lies just above the dissociation energy are determined using single photon photodissociation by excitation of the 10 K ions with a single-frequency, CW IR laser in the ion trap.

Cryogenic ion infrared spectroscopy↗
DOE OSTI2023

A human antibody epitope map of Pfs230D1 derived from analysis of individuals vaccinated with a malaria transmission-blocking vaccine

Pfs230 domain 1 (Pfs230D1) is an advanced malaria transmission-blocking vaccine antigen demonstrating high functional activity in clinical trials. However, the structural and functional correlates of transmission-blocking activity are not defined. Here, we characterized a panel of human monoclonal antibodies (hmAbs) elicited in vaccinees immunized with Pfs230D1. These hmAbs exhibited diverse transmission-reducing activity, yet all bound to Pfs230D1 with nanomolar affinity. We compiled epitope-binning data for seventeen hmAbs and structures of nine hmAbs complexes to construct a high-resolution epitope map and revealed that potent transmission-reducing hmAbs bound to one face of Pfs230D1, while non-potent hmAbs bound to the opposing side. The structure of Pfs230D1D2 revealed that non-potent transmission-reducing epitopes were occluded by the second domain. The hmAb epitope map delineated binary hmAb combinations that synergized for extremely high-potency, transmission-reducing activity. Finally, this work provides a high-resolution guide for structure-based design of enhanced immunogens and informs diagnostics that measure the transmission-reducing response.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Structure of 3-mercaptopropionic acid dioxygenase with a substrate analog reveals bidentate substrate binding at the iron center

Thiol dioxygenases are a subset of nonheme iron oxygenases that catalyze the formation of sulfinic acids from sulfhydryl-containing substrates and dioxygen. Among this class, cysteine dioxygenases (CDOs) and 3-mercaptopropionic acid dioxygenases (3MDOs) are the best characterized, and the mode of substrate binding for CDOs is well understood. However, the manner in which 3-mercaptopropionic acid (3MPA) coordinates to the nonheme iron site in 3MDO remains a matter of debate. A model for bidentate 3MPA coordination at the 3MDO Fe-site has been proposed on the basis of computational docking, whereas steady-state kinetics and EPR spectroscopic measurements suggest a thiolate-only coordination of the substrate. To address this gap in knowledge, we determined the structure of Azobacter vinelandii 3MDO (Av3MDO) in complex with the substrate analog and competitive inhibitor, 3-hydroxypropionic acid (3HPA). The structure together with DFT computational modeling demonstrates that 3HPA and 3MPA associate with iron as chelate complexes with the substrate-carboxylate group forming an additional interaction with Arg168 and the thiol bound at the same position as in CDO. A chloride ligand was bound to iron in the coordination site assigned as the O 2 -binding site. Supporting HYSCORE spectroscopic experiments were performed on the (3MPA/NO)-bound Av3MDO iron nitrosyl (S = 3/2) site. In combination with spectroscopic simulations and optimized DFT models, this work provides an experimentally verified model of the Av3MDO enzyme–substrate complex, effectively resolving a debate in the literature regarding the preferred substrate-binding denticity. These results elegantly explain the observed 3MDO substrate specificity, but leave unanswered questions regarding the mechanism of substrate-gated reactivity with dioxygen.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A broad specificity β-propeller enzyme from Rhodopseudomonas palustris that hydrolyzes many lactones including γ-valerolactone

Lactones are prevalent in biological and industrial settings, yet there is a lack of information regarding enzymes used to metabolize these compounds. One compound, γ-valerolactone (GVL), is used as a solvent to dissolve plant cell walls into sugars and aromatic molecules for subsequent microbial conversion to fuels and chemicals. Despite the promise of GVL as a renewable solvent for biomass deconstruction, residual GVL can be toxic to microbial fermentation. Here, we identified a Ca 2+ -dependent enzyme from Rhodopseudomonas palustris (Rpa3624) and showed that it can hydrolyze aliphatic and aromatic lactones and esters, including GVL. Maximum-likelihood phylogenetic analysis of other related lactonases with experimentally determined substrate preferences shows that Rpa3624 separates by sequence motifs into a subclade with preference for hydrophobic substrates. Additionally, we solved crystal structures of this β-propeller enzyme separately with either phosphate, an inhibitor, or a mixture of GVL and products to define an active site where calcium-bound water and calcium-bound aspartic and glutamic acid residues make close contact with substrate and product. Our kinetic characterization of WT and mutant enzymes combined with structural insights inform a reaction mechanism that centers around activation of a calcium-bound water molecule promoted by general base catalysis and close contacts with substrate and a potential intermediate. Similarity of Rpa3624 with other β-propeller lactonases suggests this mechanism may be relevant for other members of this emerging class of versatile catalysts.

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

Characterization of two distinct neutrophil serine protease–binding modes within a Staphylococcus aureus innate immune evasion protein family

Extracellular adherence protein domain (EAPs) proteins are a class of innate immune evasion proteins secreted by the human pathogen Staphylococcus aureus. EAPs are potent and selective inhibitors of cathepsin-G (CG) and neutrophil elastase (NE), which are the two most abundant neutrophil serine proteases (NSPs). Previous work from our group has shown that the prototypical EAP, EapH1, relies on plasticity within a single inhibitory site to block the activities of CG and NE. However, whether other EAPs follow similar structure–function relationships is unclear. To address this question, we studied the inhibitory properties of the first (Eap1) and second (Eap2) domains of the modular extracellular adherence protein of S. aureus and determined their structures when bound to CG and NE, respectively. We observed that both Eap1 and Eap2 displayed time-dependent inhibition of CG (on the order of 10 -9 M) and of NE (on the order of 10 -10 M). We also found that whereas the structures of Eap1 and Eap2 bound to CG showed an overall inhibitory mode like that seen previously for EapH1, the structures of Eap1 and Eap2 bound to NE revealed a new inhibitory mode involving a distal region of the EAP domain. Using site-directed mutagenesis of Eap1 and Eap2, along with enzyme assays, we confirmed the roles of interfacial residues in NSP inhibition. Taken together, our work demonstrates that EAPs can form structurally divergent complexes with two closely related serine proteases and further suggests that certain EAPs may be capable of inhibiting two NSPs simultaneously.

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

Identification and quantification of distinct active sites in Hf-Beta zeolites for transfer hydrogenation catalysis

Despite the significant progress made in characterizing different framework heteroatom sites that exist in Lewis acidic zeolites with probe molecule adsorption and spectroscopy, methods to reliably quantify site counts remain indefinite and have been primarily limited to Sn and Ti Lewis acid sites. Here, methods to quantify framework Lewis acidic Hf 4+ sites in zeolite Beta (Hf-Beta) with two Lewis base titrants (pyridine, deuterated acetonitrile) were developed using infrared (IR) spectroscopy. Lewis acid site counts for Hf-Beta zeolites were validated by measuring integrated molar extinction coefficients (IMECs; ε, cm μmol –1 ) on Sn-Beta zeolites using identical Lewis base titrants to benchmark site counts with established literature procedures to quantify Lewis acid sites in Beta zeolites from IR spectra. Highlighting the importance of benchmarking active site counts against well-established experimental protocols, IMECs of CD 3 CN bound to open (ε(Sn; 2316 cm –1 ): 1.80 ± 0.25) and closed (ε(Sn; 2308 cm –1 ): 3.76 ± 0.33) Sn sites were ~1.8x larger than those previously reported while total Lewis acid site counts agreed with those measured by pyridine (ε(Sn; 1451 cm –1 ): 1.58 ± 0.16) on six different Sn-Beta zeolites. IMECs measured for IR peaks reflecting pyridine bound to Lewis acidic Hf sites (ε(Hf; 1448 cm –1 ): 1.54 ± 0.21) and CD 3 CN bound to open (ε(Hf; 2313 cm –1 ): 2.40 ± 0.22) and closed (ε(Hf; 2307 cm –1 ): 3.55 ± 0.41) Hf sites, gave similar counts for the total number of Lewis acidic sites across six Hf-Beta zeolites (Si/Hf = 100–413). Consistent with previous reports with Sn-Beta catalysts where open Sn sites are responsible for catalytic turnover, apparent first and zero-order MPVO rate constants (0.01–1 M cyclohexanone in 2-butanol; per total Hf, 373 K) correlated with the total number of open Hf sites, per total Hf, but not with the total number of closed Hf sites or total Lewis acid site counts. Measured initial MPVO rates (0.1 M cyclohexanone in 2-butanol, per open Hf, 373 K) were ~25x higher on hydrophobic Hf-Beta-F than on hydrophilic Hf-Beta-OH zeolites. Overall, the apparent first-order MPVO rate constants (2-butanol solvent, per open Hf, 373 K) were ~6x higher on Hf-Beta-F than on Hf-Beta-OH zeolites. The characterization methods reported here enable normalization of MPVO turnover rates on Sn- and Hf-Beta zeolites by their number of open sites. Finally, this enables performing quantitative rate comparisons across Lewis acid zeolites of varying active site identity, solvation, and pore topology used in liquid-phase catalysis.

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

Accelerating high-order continuum kinetic plasma simulations using multiple GPUs

Kinetic plasma simulations solve the Vlasov-Poisson or Vlasov-Maxwell equations to evolve scalar-variable distribution functions in position-velocity phase space and vector-variable electromagnetic fields in configuration space. The immense computational cost of evolving high-dimensional variables, and their large number of degrees of freedom, often limits the utility of continuum kinetic simulations and presents a challenge when it comes to accurately simulating real-world physical phenomena. To address this challenge, we present techniques that accelerate and minimize the computational work required for a scalable Vlasov-Poisson solver. We show theoretical hardware compute and communication bounds for solving a fourth-order finite-volume Vlasov-Poisson system. These bounds are then used to inform and evaluate the design of performance portable algorithms for a multiple graphics processing unit (GPU) accelerated version of the Vlasov-Poisson solver VCK-CPU [1]. We demonstrate that the multi-GPU Vlasov solver implementation, VCK-GPU, simultaneously minimizes required inter-process data transfer while also being bounded by the machine network performance limits. This results in an overall strong scaling speedup per timestep of up to 40x in three-dimensional phase space (one position, two velocity coordinates) and 54x in four dimensional phase space (two position, two velocity coordinates) and a 341x increase in simulation throughput of the GPU accelerated code over the existing CPU code. The GPU code is also able to weak scale up to 256 compute nodes and 1024 GPUs. In conclusion, we demonstrate that the improved compute performance enables exploring configurations which were previously computationally infeasible, including resolving fine-scale distribution function filamentation and multi-species dynamics with realistic electron-proton mass ratios.

Continuum kinetics↗
DOE OSTI2023

Photoluminescence of CdTe and CdZnTe compounds doped with 2% selenium

The addition of 2 % of selenium in CdTe and CdZnTe (CZT) exhibited an additional acceptor-bound exciton in photoluminescence (PL). From the PL and inductively coupled plasma mass spectrometer (ICP/MS) analyses, the additional acceptor-bound exciton emission is associated with a Cu impurity in the relatively low-purity (5N) precursor CdSe starting material. Additionally, thermal quenching of the acceptor/donor bound exciton emissions do not require a shift to any higher excited states or to the free exciton state. ICP/MS analysis also showed a higher concentration of Sn in CdZnTeSe (CZTS) compared to CdTe and CZT which is known to form a deep level trap at around E C – 0.83 eV. The trapping and de-trapping times indicate that the Sn-related deep trap acts as an electron trapping center. The relatively low electron mobility-lifetime product of CZTS might be correlated with the concentration of Sn in the impure CdSe material. Pulse height spectra using an Eu-152 gamma-ray source was taken for a 2-mm-thick planar CZTS detector. The spectra clearly showed all the expected gamma peaks with energy lower than 344.3 keV.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Small uranium and oxygen interstitial clusters in UO 2 : An empirical potential study

The aggregation of defect clusters in UO 2 under irradiations creates dislocation loops and cavities, which in turn contribute to swelling of the fuels, leading to detrimental pellet-clad mechanical interactions. The kinetics of the defect clusters, however, have received less attention and are not well understood. As such, using two interatomic potentials, the Morelon potential and a many-body model by Cooper, Rushton, and Grimes (CRG), the energetics and kinetics of small uranium interstitial clusters, with and without bound oxygen interstitials, are investigated with molecular dynamics and accelerated molecular dynamics simulations. For the single bound anti-Schottky interstitial cluster (U i 2O i ), both potentials find the interstitialcy migration mechanism to be a low barrier migration mechanism. However, for the CRG potential, an oxygen Frenkel pair assisted interstitial cluster migration mechanism with a lower migration barrier is identified during accelerated molecular dynamics simulations. For the Morelon potential, the interstitial clusters have geometrically compact structures, and small interstitial clusters (two and three anti-Schottky clusters) are more mobile than the single bound anti-Schottky interstitial cluster; for the CRG potential, the interstitial clusters have geometrically non-compact structures, yielding a more complex energy landscape. The rapid migration of interstitials and interstitial clusters are important for the irradiation response of UO 2 and correlates well with the very low barriers reported in a few historic experimental studies.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2021

A review of rare earth elements and yttrium in coal ash: Content, modes of occurrences, combustion behavior, and extraction methods

Rare earth elements and yttrium (REY) have attracted considerable attention over the last decade because of their vital roles in clean energy, consumer product, national defense and security applications, among other uses. Due to the retention of REY during coal burning, coal combustion ash is considered as potential alternative sources for REY. Understanding the content, speciation, retention and/or transformation behavior of REY during coal combustion not only expands our knowledge of the combustion behavior of the trace elements in coal, but also provides basis for modeling REY partitioning during coal combustion and for developing economically viable REY recovery technologies. This review makes a critical summary of recent progress in the study of REY in coal ash. The contents and the extraction potentials of REY in coal ash derived from 15 major coal-producing countries worldwide were summarized and evaluated. Various analytical methods for determining REY bulk contents and speciation, together with the solid sample pretreatment, analytical accuracy and precision, advantages and disadvantages were summarized and compared. Modern analytical approaches combined indirect methods (e.g., sequential extraction) shed light on the physical distribution, mineralogy, and the chemical state of REY in coal ash. Three types of REY occurrences in coal ash, including Si-Al glassy association, discrete minerals or compounds, and organic association (bound with unburned carbon) were defined in the review. The glassy association can be further divided into REY minerals closely bound to glass phases and dispersed throughout the glassy structure. REY partitioning in various emission streams, the size distribution, and their enrichment behavior in coal ash were discussed. Additionally, thermal behavior and transformation of various REY forms in coal during combustion process, including organic-associated REY, REY phosphates, REY carbonates, clay-bound REY and among others were summarized. Two possible retention mechanisms of REY by aluminosilicate glass at boiler temperature were proposed: the incorporation of the individual REY phases into the glass as inclusions and the diffusion of REY phases throughout Si-Al glass structures in the melting process. Feed coal mineral types, mineral-mineral associations, boiler conditions, and other factors control the retention process. After coal combustion, the speciation of REY in fly ash may be modified by the reactions of REY phases with flue gas components. Further, an overview of REY transformation mechanisms during coal combustion was deeply discussed. Finally, current extraction techniques for REY recovery from coal combustion ash were introduced. Future outlooks and research problems were also identified.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗