Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Analytical Balance”

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.

900 records · Page 50
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

N -Arylimide Molecular Balances: A Comprehensive Platform for Studying Aromatic Interactions in Solution

Noncovalent interactions of aromatic surfaces play a key role in many biological processes and in determining the properties and utility of synthetic materials, sensors, and catalysts. However, the study of aromatic interactions has been challenging because these interactions are usually very weak and their trends are modulated by many factors such as structural, electronic, steric, and solvent effects. Recently, N-arylimide molecular balances have emerged as highly versatile and effective platforms for studying aromatic interactions in solution. These molecular balances can accurately measure weak noncovalent interactions in solution via their influence on the folded–unfolded conformational equilibrium. The structure (i.e., size, shape, π-conjugation, and substitution) and nature (i.e., element, charge, and polarity) of the π-surfaces and interacting groups can be readily varied, enabling the study of a wide range of aromatic interactions. These include aromatic stacking, heterocyclic aromatic stacking, and alkyl–π, chalcogen–π, silver–π, halogen–π, substituent–π, and solvent–π interactions. The ability to measure a diverse array of aromatic interactions within a single model system provides a unique perspective and insights as the interaction energies, stability trends, and solvent effects for different types of interactions can be directly compared. Some broad conclusions that have emerged from this comprehensive analysis include: (1) The strongest aromatic interactions involve groups with positive charges such as pyridinium and metal ions which interact with the electrostatically negative π-face of the aromatic surface via cation–π or metal–π interactions. Attractive electrostatic interactions can also form between aromatic surfaces and groups with partial positive charges. (2) Electrostatic interactions involving aromatic surfaces can be switched from repulsive to attractive using electron-withdrawing substituents or heterocycles. These electrostatic trends appear to span many types of aromatic interactions involving a polar group interacting with a π-surface such as halogen–π, chalcogen–π, and carbonyl–π. (3) Nonpolar groups form weak but measurable stabilizing interactions with aromatic surfaces in organic solvents due to favorable dispersion and/or solvophobic effects. Furthermore, a good predictor of the interaction strength is provided by the change in solvent-accessible surface area. (4) Solvent effects modulate the aromatic interactions in the forms of solvophobic effects and competitive solvation, which can be modeled using solvent cohesion density and specific solvent–solute interactions.

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

The energetics and dynamics of confinement in flexible frameworks and molecular confinement

Porous frameworks form the chemical and structural basis for critical technologies in separations, catalysis, nuclear waste containment and biomedical applications. Hundreds of zeolites and metal organic frameworks (MOFs) have been synthesized, and their ability to separate and store hydrogen, methane and carbon dioxide has been investigated both experimentally and theoretically. Nevertheless, a fundamental and systematic molecular-level understanding of the thermodynamic and structural factors governing the stability and guest-host interactions in these materials lags behind focused studies of specific systems. Because the guest molecules interact with each other and with the host framework, molecular confinement is a finely balanced and complex phenomenon. The ability of the guest molecules to bind and diffuse through the pores is determined by the nature of the host framework which, in turn, responds to the nature and concentration of guest molecules and to pressure and temperature. The work explores how framework flexibility, tailored by structure, composition, temperature and pressure, is a general phenomenon, similar in nature but variable in extent, in both zeolites and MOFs and is part of a free energy landscape in which framework-guest interactions, pressure, and temperature result in changes in framework geometry and, in some cases, phase transitions. These subtle and/or pronounced changes in lattice geometry, energetics, and dynamics can play a decisive role in confinement and in differentiating the binding of molecules of similar size. The free energy landscape created by these structural changes links polymorphism, amorphization, “breathing,” “gate opening” and confinement. Specifically, the generality of such behavior arises from commonalities in lattice dynamics and energetics of frameworks containing a combination of strong rigid bonds and weaker more flexible deformation modes. Identifying and describing these common and collective phenomena is the focus of the research on a selected group of zeolites and MOFs. Structural studies using X-ray and neutron diffraction explore the mechanical functionality of these important materials, specifically how framework materials respond to changes in temperature, pressure and guest loading. These structural studies are combined with calorimetric measurements, using techniques uniquely developed in the participating laboratories, of heats of formation, heat capacities and entropies, and guest-host interactions. The experimental thermodynamic studies are complemented by inelastic neutron scattering studies of the lattice dynamics related both to framework vibrations and to guest-host interactions. The report below summarizes the work done at UC Davis, which is complemented by work done at BYU (Woodfield) and Virginia Tech (Ross).

36 MATERIALS SCIENCE↗
DOE OSTI2025

Identifying High Ionic Conductivity Compositions of Ionic Liquid Electrolytes Using Features of the Solvation Environment

Binary mixtures of ionic liquids with molecular solvents are gaining interest in electrochemical applications due to the improvement in their performance over neat ionic liquids. Dilution with suitable molecular solvents can reduce the viscosity and facilitate faster diffusion of ions, thereby yielding substantially higher ionic conductivity than that for a pure ionic liquid. Although viscosity and diffusion coefficients typically behave as monotonic functions of concentration, ionic conductivity often passes through a peak value at an optimum molar ratio of the molecular solvent to the ionic liquid. The ionic conductivity maximum is generally explained in terms of a balance between the ease of charge transport and the concentration of the charge carriers. In this work, fluctuation in the local environment surrounding an ion is invoked as a plausible explanation for the ionic conductivity mechanism with a binary mixture of 1-ethyl-3-methylimidazolium tetrafluoroborate and ethylene glycol as an example. The magnitude of the dynamism in the local environment is captured by measuring the spatial and temporal features of the solvation environment. Standard deviation in the number of ions in the solvation environment serves as a spatial feature, while the cage correlation lifetimes for oppositely charged ions within the first solvation shell serve as a temporal feature. Large standard deviations in the cluster ion population and short cage correlation lifetimes are indicators of highly dynamic ionic environment at the molecular level and consequently yield high ionic conductivity. Such compositions were found to be in good agreement with the optimum ionic liquid mole fractions obtained through experimental measurement. Short cage correlation lifetimes enable the identification of optimum mixture compositions using simulation trajectories significantly shorter than those required to implement the Nernst–Einstein or Einstein formalisms for calculating ionic conductivity. We validated the applicability of this approach across force fields and in six ionic liquid-molecular solvent electrolytes formed with combination of cations, anions, and solvents. We offer a computationally efficient approach of screening ionic liquid-molecular solvent binary mixture electrolytes to identify molar ratios that yield high ionic conductivity.

25 ENERGY STORAGE↗
DOE OSTI2024

HydraGNN_Predictive_GFM_2024 - Ensemble of predictive graph foundation models for ground state atomistic materials modeling

We provide the ensemble of fifteen pre-trained graph foundation models (GFMs) for atomistic materials modeling applications. Each one of the fifteen GFMs has been trained on five open-source datasets that (once aggregated) amount to over 154 million atomistic structures, which cover over two-thirds of the natural elements of the periodic table and that comprises a broad set of organic and inorganic compounds. This vast set of atomistic structures comprises ground state configurations that are dynamically stable (i.e., equilibrated structures with atomic forces approximately close to zero values) as well as dynamically unstable structures (i.e., non-equilibrium structures with non-negligible non-zero values of atomic forces). The ensemble of datasets aggregated does NOT include excited states. The datasets have been curated to remove atomistic structures with spectral norm of the force tensor above 100 eV/angstrom. Moreover, a linear term of the energy was computed for each dataset using a linear regression model that uses the chemical concentration of each natural element as regressor. The linear term predicted by the linear regression model has been subtracted from each original energy value to perform a re-alignment of the energy values across different electronic structures approximation theories performed to generate the diverse multi-source, multi-fidelity datasets. The folder "ADIOS_files" contains the set of pre-processed datasets in Adaptable I/O System (ADIOS) format (https://www.exascaleproject.org/research-project/adios/) that have been used for the development and training of GFMs in this work. The "ADIOS_files" directory contains 6 sub-directories named as follows: - ANI1x-v3.bp - MPTrj-v3.bp - OC2020-20M-v3.bp - OC2020-v3.bp - OC2022-v3.bp - qm7x-v3.bp Each sub-directory contains the pre-processed datasets converted in Adaptable I/O System (ADIOS) format (https://www.exascaleproject.org/research-project/adios/) that have been used to the development, training, and performance testing of the ensemble go predictive graph foundation models. Each GFM was developed using HydraGNN (https://github.com/ORNL/HydraGNN) as underlying graph neural network (GNN) architecture. The multi-task learning (MTL) capability of HydraGNN was used to simultaneously train the GFMs on labeled values for direct predictions of energy (a total system property of an atomistic structure that measures the chemical stability) and atomic forces (an atomic level property of an atomistic structure that measures the dynamical stability). The hyper parameters of the GFM have been tuned using scalable hyperparameter optimization (HPO) algorithms implemented in the software DeepHyper (https://github.com/deephyper/deephyper). The pre-training of each HPO trial was performed using distributed data parallelism (DDP) to scale the training across 128 compute nodes of the exascale OLCF supercomputer Frontier. Each HPO trial was trained only for 10 epochs and an early stopping was performed to avoid wasting significant computational resources on GNN architectures that were clearly underperforming. For each HPO trial, the 'omnistat' tool developed by (AMD Research - Advanced Micro Device) was used to measure the total energy consumption in kWh. The ensemble of GFMs was obtained by selecting the fifteen best performing HPO trials. Four models have been selected for their clear advantage in accuracy, and these are the GFMs with IDs 229, 156, 147, 260. Additional eleven models have been selected based on judicious balance between accuracy and energy consumption needed for training, and these are the GFMs with IDs 165, 78, 137, 1, 175, 171, 181, 67, 179, 167, 351. Each selected GFM of the ensemble was continued to cumulate a total of at most 30 epochs. In some cases, the total number of epochs actually performed was les than 30 due to two combined factors: (1) the size of the GFM (i.e., the number of model parameters to train) and (2) the total wall-clock time for which the computational resources could be allocated on OLCF-Frontier. The "Ensemble_of_models" directory contains 15 sub-directories named as follows: - gfm_0.229 - gfm_0.156 - gfm_0.147 - gfm_0.260 - gfm_0.165 - gfm_0.78 - gfm_0.137 - gfm_0.1 - gfm_0.175 - gfm_0.171 - gfm_0.181 - gfm_0.67 - gfm_0.179 - gfm_0.167 - gfm_0.351 Each one of these sub-directories refers to one of the fifteen HPO trials that have been selected to continue the pre-training with at most 30 epochs. With each sub-directory associated with a specific HPO trial, the following files can be found: - config.json: file for argument parsing to develop and train an HydraGNN architecture - gfm_0.ID_epoch_N.pk: file with model parameters for HPO ID trial after N epochs of training The ensemble of fifteen GFM architectures was used for (1) ensemble averaging to stabilize the predictions of energy and atomic forces after pre-training for post-processing analysis and (2) ensemble uncertainty quantification (UQ). The code used to develop, pre-train, and load the pre-trained models for post-processing analysis is available on the ORNL-GitHub at the following link: https://github.com/ORNL/HydraGNN/tree/Predictive_GFM_2024

36 MATERIALS SCIENCE↗
DOE OSTI2023

Chemistry and Dynamics of Supercritical Carbon Dioxide and Methane in the Slit Pores of Layered Silicates

In the mid 2010s, high-pressure diffraction and spectroscopic tools opened a window into the molecular-scale behavior of fluids at the conditions of many CO 2 sequestration and shale/tight gas reservoirs - conditions where CO 2 and CH 4 are present as variably wet supercritical fluids. Integrating high-pressure spectroscopy and diffraction with molecular modeling has revealed much about the ways that supercritical CO 2 and CH 4 behave in reservoir components, particularly in the slit-shaped micro- and meso-pores of layered silicates (phyllosilicates) abundant in caprocks and shales. This Account summarizes how supercritical CO 2 and CH 4 behave in the slit pores of swelling phyllosilicates as functions of the H 2 O activity, framework structural features, and charge-balancing cation properties at 90 bar and 323 K – conditions similar to a reservoir at ~1 km depth. Slit pores containing cations with large radii, low hydration energy, and large polarizability readily interact with CO 2 , allowing CO 2 and H 2 O to adsorb and co-exist in these interlayer pores over a wide range of fluid humidities. In contrast, cations with small radii, high hydration energy, and low polarizability weakly interact with CO 2 leading to reduced CO 2 uptake and a tendency to exclude CO 2 from interlayers when H 2 O is abundant. The reorientation dynamics of confined CO 2 depends on the interlayer pore height, which is strongly influenced by the cation properties, framework properties, and fluid humidity. The silicate structural framework also influences CO 2 uptake and behavior - for example, smectites with increasing F-for-OH substitution in the framework take up greater quantities of CO 2 . Reactions that trap CO 2 in carbonate phases have been observed in thin H 2 O films near smectite surfaces, including a dissolution-reprecipitation mechanism when the edge surface area is high and an ion exchange-precipitation mechanism when the interlayer cation can form a highly insoluble carbonate. In contrast, supercritical CH 4 does not readily associate with cations, does not react with smectites, and is only incorporated into interlayer slit mesopores when (i) the pore has a $z$-dimension large enough to accommodate CH 4 , (ii) the smectite has low charge and (iii) the H 2 O activity is low. In conclusion, the adsorption and displacement of CH 4 by CO 2 and vice-versa have been studied on the molecular scale in one shale, but opportunities remain to examine behavioral details in this more complicated, slit-pore inclusive system.

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

Steric and Lewis Basicity Influence of the Second Coordination Sphere on Electrocatalytic CO 2 Reduction by Manganese Bipyridyl Complexes

This study aims to provide a greater insight into the balance between steric (bpy vs (Ph) 2 bpy vs mes 2 bpy ligands) and Lewis basic ((Ph) 2 bpy vs (MeOPh) 2 bpy vs (MeSPh) 2 bpy ligands) influence on the efficiencies of the protonation-first vs reduction-first CO 2 reduction mechanisms with [Mn I (R 2 bpy)(CO) 3 (CH 3 CN)] + precatalysts, and on their respective transition-state geometries/energies for rate-determining C–OH bond cleavage toward CO evolution. The presence of only modest steric bulk at the 6,6'-diphenyl-2,2'-bipyridyl ((Ph) 2 bpy) ligand has here allowed unique insight into the mechanism of catalyst activation and CO 2 binding by navigating a perfect medium between the nonsterically encumbered bpy-based and the highly sterically encumbered mes2bpy-based precatalysts. Cyclic voltammetry conducted in CO 2 -saturated electrolyte for the (Ph) 2 bpy-based precatalyst [2-CH 3 CN] + confirms that CO 2 binding occurs at the two-electron-reduced activated catalyst [2] – in the absence of an excess proton source, in contrast to prior assumptions that all manganese catalysts require a strong acid for CO 2 binding. This observation is supported by computed free energies of the parent–child reaction for [Mn–Mn] o dimer formation, where increased steric hindrance relative to the bpy-based precatalyst correlates with favorable CO 2 binding. A critical balance must be adhered to, however, as the absence of steric bulk in the bpy-based precatalyst [1-CH 3 CN] + maintains a lower overpotential than [2-CH 3 CN] + at the protonation-first pathway with comparable kinetic performance, whereas an ~2-fold greater TOF max is observed at its reduction-first pathway with an almost identical overpotential as [2-CH 3 CN] + . Notably, excessive steric bulk in the mes2bpy-based precatalyst [3-CH 3 CN] + results in increased activation free energies of the C–OH bond cleavage transition states for both the protonation-first and the reduction-first pathways relative to both [1-CH 3 CN] + and [2-CH 3 CN] + . In fact, [3-CH 3 CN] + requires a 1 V window beyond its onset potential to reach its peak catalytic current, which is in contrast to the narrower (<0.30 V) potential response window of the remaining catalysts here studied. Furthermore, voltammetry recorded under 1 atm of CO 2 with 2.8 M (5%) H 2 O establishes [2-CH 3 CN] + to have the lowest overpotential (η = 0.75 V) in the series here studied, attributed to its ability to lie “on the fence” when providing sufficient steric bulk to hinder (but not prevent) [Mn–Mn] o dimerization, while simultaneously having a limited steric impact on the free energy of activation for the rate-determining C–OH bond cleavage transition state. While the methoxyphenyl bpy-based precatalyst [4-CH 3 CN] + possesses an increased steric presence relative to [2-CH 3 CN] + , this is offset by its capacity to stabilize the C–OH bond cleavage transition states of both the protonation-first and the reduction-first pathways by facilitating second coordination sphere H-bonding stabilization.

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

Ca X ML: Chemistry‐informed machine learning explains mutual changes between protein conformations and calcium ions in calcium‐binding proteins using structural and topological features

Proteins' flexibility is a feature in communicating changes in cell signaling instigated by binding with secondary messengers, such as calcium ions, associated with the coordination of muscle contraction, neurotransmitter release, and gene expression. When binding with the disordered parts of a protein, calcium ions must balance their charge states with the shape of calcium-binding proteins and their versatile pool of partners depending on the circumstances they transmit. Accurately determining the ionic charges of those ions is essential for understanding their role in such processes. However, it is unclear whether the limited experimental data available can be effectively used to train models to accurately predict the charges of calcium-binding protein variants. Here, we developed a chemistry-informed, machine-learning algorithm that implements a game theoretic approach to explain the output of a machine-learning model without the prerequisite of an excessively large database for high-performance prediction of atomic charges. We used the ab initio electronic structure data representing calcium ions and the structures of the disordered segments of calcium-binding peptides with surrounding water molecules to train several explainable models. Network theory was used to extract the topological features of atomic interactions in the structurally complex data dictated by the coordination chemistry of a calcium ion, a potent indicator of its charge state in protein. Our design created a computational tool of Ca X ML, which provided a framework of explainable machine learning model to annotate ionic charges of calcium ions in calcium-binding proteins in response to the chemical changes in an environment. Our framework will provide new insights into protein design for engineering functionality based on the limited size of scientific data in a genome space.

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

The Influence of Nonlocal Metal–Dielectric Environments on Physical and Chemical Processes

Plasmonic nanolayers and laminar metallic/dielectric multilayers were originally developed for optical cloaking applications and lensing applications that could potentially image objects whose size was below the diffraction limit. These assemblies were initially formed from gold or silver nanorods grown within an alumina mesh. However, more recently, assemblies with similar properties have also been prepared by sequential thin-layer deposition of alternating layers of gold and magnesium fluoride (MgF 2 ). These metal/dielectric composite materials enable control of the dielectric constant in the directions perpendicular to the layers and balance the real and imaginary dielectric constants of the assembly such that the speed and the amplitude of the waves traveling through the assembly are not attenuated. In this Account, we will also focus on a few of the applications ranging from surface wetting to fluorescence quenching to enhancement of photochemical reactions. First, we will share an introduction to processes used to create these materials, which are combinations of low refractive index metals and transparent higher index materials arranged in a scalable repeating fashion. Two fabrication methods were employed: an electrochemical deposition of Ag nanorods into an anodized alumina matrix which produced materials with an anisotropic negative refractive index material within the plane of the film and lamellar metal/dielectric layers in which the negative index perpendicular to the growth direction. These alternating layers of plasmonic metals and dielectric materials were ultimately chosen to prepare films for further testing, because of their relative ease of fabrication. We will continue with a discussion of a few of the applications of both of these nonlocal dielectric composite materials including more specialized plasmonic, composite, and hyperbolic metamaterials including fluorescence quenching, photochemical reactions, and surface wetting. In each of these applications, the unique response caused by the enhancement of the electric field and the interface between hyperbolic materials and plasmonic materials as they interact photophysically with their near neighbors is presented. In each of the applications, the enhanced electric field extends from the composite substrate layer to interact with its near neighbors and beyond. Furthermore, the presence of this extended interaction can be observed in the form of decreased emission lifetime, enhancement of photochemical reaction rates, and changes in the surface energies measured by contact angle goniometry. In this Account, all of these situations will be addressed. Finally, we will conclude with a summary and vision for the future as well as a discussion of the unique challenges and opportunities available as research active faculty at an HBCU.

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

Investigation of design principles for metal-binding and conductive protein assemblies

Throughout the lifetime of this initiative, including renewals, we focused on understanding the fundamental principles of protein-protein interface design that enable predictable and modular spatial and kinetic control of multi-component protein self-assembly in 1D, 2D, and 3D, including the interface with inorganic materials, small molecules, and metal ions. We designed individual protein components that bind specific metal ions, including REEs and transport ions across lipid membranes. We created helical 1D filaments of repeating units with programmed periodicity, pitch, and multi-component environmentally responsive self-assembling protein fibers. We showed that these filaments reversibly assemble and disassemble under specific pH conditions and created end-specific caps that independently tune the balance of attachment and detachment rates at each terminus of the filament. Using similar filaments, we succeeded in binding arrays of heme and chlorophyll molecules and assembling patterned helical coatings around carbon nanotubes in efforts to create de novo conductive nanowires. By arraying REE binding sites in a large circular tandem array with a repeat protein-based cyclic oligomer, we created a molecular scaffold for superradiance and paramagnetic quantum sensing. We created a range of one-component and two-component self-assembling 2D arrays and showed that when designed to engage cell receptors, these arrays can control cell behavior from outside the cell signal to inside the cell. We designed helical repeat proteins with variable lengths displaying charged residues in a pattern matched to the cation lattice of mica. achieved a range of ordered states with an epitaxial match to the underlying crystal lattice. We further applied the learned principles of protein-induced biomineralization to design proteins with an interface lattice matching CaCO 3 and guide the formation of specific crystal forms of CaCO 3 from solution, a significant advance toward the global need to manage carbon. In all cases of mineral lattice matching and biomineralization, we followed assembly using molecularly resolved in situ AFM imaging and extracted information about assembly pathways and energetics, applying deep learning to quantify the dynamics of protein self-organization. We developed techniques for using dynamic metal-dependent interfaces on protein nanopores for discriminatively sensing dilute REEs in solution and demonstrated the use of strong metal-binding interfaces to drive nanocage disassembly for conditional nanocompartmentalization applications. This grant supported 11 people, including Asim Bera, Evans Brackenbrough, Andrew Borst, Nikita Hanikel, Timothy Huddy, Emily Joyce, Alex Young-Seug Kang, Ryan Kibler, Joshua Morris Lubner, Harley Pyles, and Shuai Zhang. The research effort culminated in the production of published papers and theses. Electronic Thesis/Dissertation are distributed by ProQuest/UMI Dissertation Publishing and made available on an open access basis through UW Libraries ResearchWorks Service.

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

Solar cell contacts: quantifying the impact of interfacial layers on selectivity, recombination, charge transfer, and V oc

Interfacial layers (IFLs) are ubiquitous in solar cells, but their precise impact on carrier transfer rates and the relation of these rates to performance metrics and the concepts of selectivity and recombination is lacking. We report the use of a well-defined interdigitated back-contact (IBC) silicon solar cell to determine the precise role of such IFLs. We characterize the action of the common IFL spiro-OMeTAD by making it a third contact to the IBC cell. This architecture creates three solar cells within a single structure that, with numerical simulation, provide the exchange current densities (i.e., charge transfer rates) for electrons ( J 0n ) and holes ( J 0p ) and the quasi Fermi-level splitting in the absorber, which measures the balance of generation and recombination. Further, we describe the relation of V oc to contact recombination, the asymmetry in electron/hole collection rates at a single contact (contact selectivity), and the asymmetry in collection rates of the same carrier at separate contacts (carrier selectivity). Relative to bare gold, neat spiro-OMeTAD reduces J 0n and J 0p (their geometric mean ( J 0n J 0p ) 0.5 decreases by 10 4 ), decreasing contact recombination. Addition of the common dopant Li-TFSI and air increase J 0p / J 0n by 10 6 with little effect on ( J 0n J 0p ) 0.5 , increasing contact hole selectivity. The significant increase in V oc observed by introducing spiro-OMeTAD/Li-TFSI IFLs into the cells studied, however, is due to an increase in the carrier selectivity rather than the contact selectivity or recombination of the spiro-OMeTAD-modified contact. Operando measurements further show voltage-dependent changes in the J 0 s, demonstrating that spiro-OMeTAD contributes to hysteresis. Broader context : Solar cells operate by photogenerating excess charge carriers in an absorber material and, in competition with recombination, asymmetrically extracting them at so-called carrier-selective contacts, one that ideally collects only electrons and the other, only holes. Particularly in emerging photovoltaics such as perovskites, thin layers of organic semiconductors or related materials are introduced between the absorber and contact to improve power conversion efficiency. In terms of interfacial charge transfer, a prevailing view is that such interfacial layers improve performance by helping block the collection of the undesired carrier, considered a form of recombination. We use a novel platform to study the simultaneous impact of spiro-OMeTAD, a common interfacial layer, on the collection of electrons, the collection of holes, and on the recombination of electrons and holes. We quantitatively demonstrate that spiro-OMeTAD layers indeed passivate the contact toward recombination, but that larger improvements in the open-circuit voltage, a key cell metric, can come not only from this but also from increasing the asymmetry of the collection of electrons in the system. Further, operando measurements show transient changes in the properties of spiro-OMeTAD which suggest that it contributes to hysteresis phenomena commonly observed in perovskite and other solar cells.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Composition, Activity, and Stability of IrO x Oxygen Evolution Reaction Electrocatalysts

The oxygen evolution reaction (OER) is integral to several electrochemical energy conversion and storage technologies, including carbon dioxide reduction to value added fuels, nitrogen reduction to ammonia, reversible fuel cells, rechargeable metal−air batteries, and water electrolysis to produce hydrogen. Iridium oxide (IrO x ) is widely recognized as the benchmark OER catalyst for acidic environments. Despite widespread use of IrO x catalysts, most notably in proton-exchange membrane water electrolyzers (PEMWEs), a comprehensive understanding of the physicochemical properties of commercial catalysts and the impact of these properties on both the activity and stability of these catalysts is lacking. Here, we study commercial IrO x catalysts with different physicochemical properties, three nominally considered amorphous and three rutile, to elucidate how structural and compositional variations affect OER activity and stability. Utilizing standardized aqueous electrochemical protocols, time-resolved dissolution quantification using inductively-coupled plasma mass spectrometry, and physicochemical characterization, including multiple synchrotron X-ray techniques, we systematically correlate catalyst properties with OER performance and degradation behavior aided by principal component analysis (PCA). Our results demonstrate the general trend of amorphous IrO x having higher intrinsic activity but limited stability and crystalline rutile IrO 2 having lower activity but enhanced stability against dissolution. The trends within the amorphous and rutile catalyst groups correlate with inherent material properties, including phase composition and structure, crystallinity, particle size, surface area, and surface structure/chemistry. Notably, we identify a rutile catalyst with the largest crystallite/ domain sizes, moderate surface area, a small fraction of hydrous phase, and a favorable pore structure (trimodal distributions of pore sizes ranging from 2−5 nm) that exhibits the best balance between activity and stability among the six catalysts studied here. These findings illustrate a fundamental structure-governed trade-off between activity and stability and highlight the critical role of surface chemistry modification and structure engineering in IrO x catalyst optimization.

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

Single-Step Conversion of Ethanol to n-Butene over Ag-ZrO2/SiO2 Catalysts

Ethanol is a promising platform molecule for production of a variety of fuels and chemicals. Of particular interest is producing middle distillate fuels (i.e., jet and diesel blendstock) from renewable ethanol feedstock. State-of-the-art alcohol-to-jet technology requires multiple process steps based on catalytic dehydration of ethanol to form ethylene, followed by sometimes a multi-step oligomerization, and then hydrotreatment and distillation. Here we report on a new catalytic route in which ethanol is directly converted to n-butene (1- and 2-butene mixtures) over Ag-ZrO2/SBA-16, thus offering the potential for a reduction in the number of required processing steps versus conventional alcohol-to-jet technology. This catalyst system provides the balanced metal and Lewis acid sites required to selectively facilitate a cascading sequence of reactions that includes dehydrogenation, aldol condensation, Meerwein–Ponndorf–Verley reduction, dehydration, and hydrogenation. High conversion and selectivity toward either n-butene or 1,3-butadiene is achieved by tuning the hydrogen feed partial pressure and other process/catalyst parameters. With sufficient hydrogen partial pressure 1,3-butadiene is completely and selectively hydrogenated to form n-butene. The reaction mechanism was elucidated through operando-nuclear magnetic resonance investigations coupled with reactivity measurements. Combined experimental-computational investigation reveals how changes in silver and zirconium composition and the silver oxidation state affects reactivity under controlled hydrogen partial pressures and after prolonged run times. Finally, catalyst effectiveness also was demonstrated when using wet ethanol feed, thus highlighting process flexibility in terms of feedstock purity requirements. This work was financially supported by the U.S. Department of Energy (DOE), Office of Energy Efficiency and Renewable Energy, Bioenergy Technologies Office, and was performed at the Pacific Northwest National Laboratory (PNNL) under Contract No. DE-AC05-76RL01830 and the National Renewable Energy Laboratory under Contract No. DE-AC36- 08GO28308. Part of the work conducted by S. A. Akhade was performed under the auspices of the U.S. DOE at Lawrence Livermore National Laboratory under Contract No. DE-AC52-07NA27344. This work was partly supported through the PNNL-WSU Distinguished Graduate Research Program for ADW. NMR and XPS experiments were performed using EMSL (grid.436923.9), a DOE Office of Science User Facility sponsored by the Office of Biological and Environmental Research

Dagle, Vanessa↗
DOE OSTI2021

Center for Gas Separations (CGS)

The total energy consumption in the U.S. has been rising steadily for decades, and it currently amounts to ~98,000 TBtu/yr, with approximately 30% of this total attributable to the industrial sector. Reasonable estimates indicate that 45–55% of total industry energy consumption derives from chemical separations, and for example, over 120 TBtu/yr alone is used in carrying out olefin/paraffin separations via energy-intensive cryogenic distillation. Therefore, the pursuit of new, even radically different approaches to some of the most energy-intensive industrial separations processes is an imperative scientific pursuit for reducing energy consumption toward a more sustainable future. Adsorbent and membrane-based separations can require a fraction of the energy needed for distillation methods, and as such are considered promising solutions for balancing increasing energy demand in the U.S. with the need for a massive reduction in energy consumption. Although considerable research effort has been devoted to the design of materials capable of carrying out various gas separations, usually operating through size-selective, chemisorptive, or physisorptive mechanisms, it remains a great challenge to design materials that function adequately for real-world applications. Indeed, the chemical and physical differences between molecules in gas mixtures of interest are often small, and therefore it is necessary, through the use of nanoscience and synthetic chemistry, to engineer unprecedented molecular-level control in adsorbate–adsorbent interactions. The overarching mission of the Center for Gas Separations (CGS) was to discover fundamental innovations that have the potential to dramatically reduce the energy associated with critical gas separations. In particular, the CGS developed novel synthetic routes, guided by molecular chemistry principles, as well as advanced characterization and computational methods, that have enabled the discovery of new materials and membranes tailor-made to exhibit exceptional performance for a range of gas separations processes, as required in the clean use of fossil fuels and in reducing CO 2 emissions from industry. A challenge of this magnitude required the collaboration and synergy of a large team of researchers with expertise in materials synthesis, characterization, and computations. During the 11-year project period, the CGS created a range of new materials within the family of highly-tunable, porous solids known as metal–organic frameworks (MOFs). These new frameworks demonstrate novel mechanisms for key industrial gas separations, including revolutionary new cooperative adsorption processes that enable low-energy CO 2 and CO capture, and are capable of efficiently separating olefins from paraffins, O 2 from air, and the shape-selective separation of alkane isomers. In addition, the CGS developed new strategies for incorporating these materials into composite membranes toward highly efficient and selective membrane-based separations. As a testament to the success of the CGS, two start-up companies, Mosaic Materials,4 Inc. and Flux Technology, Inc., grew out of these research efforts, and these companies are seeking to commercialize MOF and composite membranes materials for key separations in industry, including large-scale CO 2 capture and hydrocarbon separations, respectively. Another company, framergy, Inc., licensed IP resulting from CGS research toward the commercialization of adsorbents for various energy-relevant applications.

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

Defense Waste Processing Facility Nitric-Glycolic Flowsheet Chemical Process Cell Chemistry: Part 2

The conversions of nitrite to nitrate, the destruction of glycolate, and the conversion of glycolate to formate and oxalate were modeled for the Nitric-Glycolic flowsheet using data from Chemical Process Cell (CPC) simulant runs conducted by Savannah River National Laboratory (SRNL) from 2011 to 2016. The goal of this work was to develop empirical correlation models to predict these values from measurable variables from the chemical process so that these quantities could be predicted a-priori from the sludge or simulant composition and measurable processing variables. The need for these predictions arises from the need to predict the REDuction/OXidation (REDOX) state of the glass from the Defense Waste Processing Facility (DWPF) melter. This report summarizes the work on these correlations based on the aforementioned data. Previous work on these correlations was documented in a technical report covering data from 2011-2015. This current report supersedes this previous report. Further refinement of the models as additional data are collected is recommended. The glass REDOX depends on the concentrations of nitrate and manganese (oxidants), and of glycolate, formate, oxalate, carbon, and antifoam (reductants) in the melter feed. The waste sludge contains nitrite, nitrate, manganese (Mn), and oxalate. Virtually all of the nitrite is converted to nitrate or NO+NO 2 +N 2 O gases in the CPC. The portion of the nitrite converted to nitrate increases the amount of nitrate in the sludge. The amount of glycolate in the final melter feed depends on the amount of the glycolic acid feed that is destroyed. Similarly, the amounts of formate and oxalate formed during the decomposition of glycolic acid are required. The material balance on carbon was found to not close in most cases. Generally, there was less carbon at the end of testing compared to the inputs. The most uncertain product variable was glycolate, so material balances were performed where the glycolate concentration was adjusted, usually upward, to close the balance. Correlation versus the original, as-measured, data was generally poor, but correlation against the material balance adjusted values was greatly improved. It was also shown that the correlation of the measured REDOX versus the predicted REDOX was much better when the material balance adjusted glycolate values were used. Three data series were primarily used during the regressions of the data; these series were 1) Sludge Batch 9 NG flowsheet simulant runs NG51-62 (SB9-NG); 2) Scaled Runs + Bounding Hydrogen Runs (SR+BH); and 3) Runs GN43-50 and 57 (43-50,57). The glycolate destruction was found to correlate with acid stoichiometry (AS), percent reducing acid (PRA), and for some data series, headspace to simulant volume ratio (HSV), mercury (Hg), and nitrate. Although glycolate destruction for pairs of data series (e.g., [SB9-NG] and [SR+BH]) were found to depend on HSV, the combination of all three data series was not found to have significant dependence on this variable. The best model for glycolate destruction depended on AS, nitrate, and Hg. This model predicted the product glycolate compositions of the data to within 92-106%. The conversion of glycolate to formate was high when noble metals and Hg were not present, with values up to 100%. When noble metals and Hg were present, this conversion ranged from zero to 7%, and was dependent on AS. Lower AS gave higher conversions to formate. The conversion to oxalate was found to depend on the AS and the initial concentration of nitrite. An alternative fit versus AS and the form of ruthenium (Ru) used is a possible alternative. This fit was somewhat less statistically significant. This second model predicts that more oxalate is formed when Ru-nitrosyl nitrate is used rather than Ru chloride. The conversion of glycolate to oxalate ranged from zero to 6%. The conversion of nitrite to nitrate depended primarily on AS and PRA, with HSV and Hg being significant when these variables were varied. For multiple series of data, nitrite was also needed to SRNL-STI-2017-00172 5HYLVLRQ viL distinguish between data series, and the effect of HSV became insignificant. The best model for nitrite to nitrate conversion depended on AS, PRA, nitrite, and Hg. The 95% confidence intervals on the predicted values of glycolate destruction, glycolate to oxalate conversion, and nitrite to nitrate conversion were used to determine the uncertainty in the predicted REDOX when starting with only the composition of the sludge, AS, and PRA. Using the 95% confidences on an individual value (that is the confidence in getting a particular value for one single test as opposed to what the mean would be for multiple tests), the uncertainty in the predicted REDOX was calculated. The uncertainty in the actual product composition glycolate, oxalate, formate, and nitrate concentrations translated to an uncertainty in the REDOX value of ±0.1,which is approximately the uncertainty claimed in the REDOX model itself.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Multi-Sourced Collaboration for the Production and Refining of Rare Elements and Critical Metals (Final Technical Report)

The project objective was to develop a feasible and cost-effective method for recovering rare earth elements (REEs) and critical materials (CMs) from coal and coal byproducts, resulting in high-purity individually separated REEs and CMs. The targeted REEs included Y, Pr, Nd, Gd, Dy, and Sm, with a purity of over 99.5%, while the CMs included Co, Mn, Ga, Sr, Li, Ni, Zn, and Ge, with a purity of over 90%. The project aimed to design a prototype facility capable of producing 1-3 tonnes/day of high-purity REO mixes. The work was divided into four designated circuits: 1) REE extraction and concentration, 2) REE separation and purification, 3) RE metal production, and 4) CM production. To achieve these goals, the project involved 11 tasks, including technology reviews, research, process flow diagram development, mass balance estimation, and preliminary technical-economic analysis. The project team included researchers from the University of Kentucky, University of Alabama and Virginia Tech as well as process specialists from Argonne National Laboratory. MP Materials provided technical support regarding rare earth markets and processing while Alliance Coal performed resource assessment. The project included a market analysis for Nd/Pr, Tb, Dy, Gd, Y, Co, Mn, Li, Sr, Ga, Ni, Zn, and Ge. These analyses provided insights into the supply and demand trends as well as historic and future projections of market price relative to purity requirements for these elements. Two coal resources were selected for the project: the West Kentucky No. 13 (Baker) Seam and an undisclosed lignite resource in the Illinois coal basin. The estimated quantities of REEs in these resources were calculated based on production samples and drilling data. It was estimated that there is adequate supply for an operation producing one metric ton daily of higher purity mixed rare earth oxides (MREO) for approximately 20 years at a site located in western Kentucky. In Circuit 1, project data was obtained from a pilot heap leach and REE concentration facility. It was concluded that the existing circuit, which generated a MREO concentrate, two types of CM mixed products, and Li- and Sr-containing waters, would be suitable feed for circuits 2-4. Data from the first-of-its-kind coal coarse refuse heap leach pilot pad played a crucial role in estimating reliable elemental concentrations of the pregnant leaching solution (PLS). The average total REE concentration in the PLS was found to be 28.6 ppm. In Circuit 2, several concepts were explored including a novel process referred to as solvent-assisted chromatography (SAC). This concept involved a novel columnar reactor that incorporated multiple mixer/settlers, thereby enabling the operation of counter-flowing aqueous and organic phases. Unfortunately, due to project time constraints, a complete fundamental modeling analysis could not be completed to fully evaluate the technology. Molten salt electrowinning was considered as an alternative for circuit 3 following circuit 2 purification circuit utilizing the novel SAC process. A mass and energy balance of Nd reduction to metal in a fluoride containing molten salt electrolyte was conducted. Comparisons were made with the current state of Asian molten salt electrorefining, and potential improvements in siphoning rare earth metals (REM) from the reactor were presented. A cost estimate was performed for the production of 1 tonne per day, which yielded a total of $2.29 million for the nine electrowinning (EW) cells required. The selected option for circuits 2 and 3 was a plasma distillation process, which initially separates rare earth elements (REEs) from other elements. This is followed by selective electrowinning in various ionic liquids. The selection was made on the basis of thermodynamic modeling and experimental data previously published by a project partner. The combination offers an innovative approach to integrated refining and RE metal production. For Circuit 4, an extensive literature review was conducted for the processing of the CMs. The ultimate decision was to utilize a combined plasma and ionic liquid process as well to produce individual high-purity concentrates of Zn, Ni, Co, Mn, and Mg. A separate flowsheet for Li and Sr was recommended, which would yield carbonates of these elements. Due to the lack of suitable experimental data at this time, a process recommendation could not be provided but several methods have been proposed for consideration. Lastly, a techno-economic analysis (TEA) was conducted to assess the effectiveness of the proposed process for further investigation. The TEA results revealed a capital expense (CapEx) of $737 million and an annual operational expense (OpEx) of $220 million. Due to the selected elements, the hypothetical heap leach pad can produce 1 metric tonne per day of REO equivalent, but a conscious decision was made to only treat targeted REEs, resulting in the production of 0.4 metric tonne of REM. An estimated annual revenue of $90.87 million was projected based on standard market pricing information provided by the funding agency. During the TEA, ten different modules were evaluated for costing purposes. The precipitation circuit was identified as the largest single operational expense, followed by the Mg/Mn process due to the amount of treated metal. In terms of capital expenditures, the heap leach process incurred the highest cost, followed by the Mg/Mn process. The scalability of the plasma process is a crucial consideration since the reactors cannot be scaled beyond the largest demonstrated size due to their reliance on surface area of the slag and vapor phase. The purity estimate for the REEs are generally 98%±2% to produce a metal. The purity level being lower than the project objective was due to the lack of specific experimental data needed to tighten the tolerance of the estimates. Based on literature and previous experience, the CMs are estimated as follows; Ga (95%+, metal), Sr (95%+, carbonate), Li (95%+, carbonate), Ni (98%±2%, metal), Zn (95%+, metal sponge), Ge (95%+, metal), Co (98%±2%, metal), and Mn (98%±2%, metal).

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2021

Minimizing Cr-Evaporation from Balance Of Plant Components By Utilizing Cost-Effective Alumina-Forming Austenitic Steels

A solid oxide fuel cell (SOFC) is a clean and efficient energy conversion device. The development of intermediate-temperature SOFCs has made it preferable to use metallic interconnects (MICs) to greatly reduce the cost and significantly increase the efficiency compared to ceramic interconnect materials. However, gaseous chromium species will evaporate from the chromium-containing layer formed on the surface of commonly used MICs and balance of plant (BoP) components. Volatile chromium species have been shown to form solid deposits which poison the cathodes of SOFCs, causing drastic cell performance degradation and thereby limiting commercialization. In order to alleviate the Cr poisoning and achieve long-term high performance of SOFC stacks, various Al 2 O 3 -forming austenitic (AFA) stainless steels applied at different temperatures are evaluated in this project. Based on our Phase I results (500 h operation), it is shown that on the AFAs, an alumina-based protective layer forms under high temperature that is invulnerable to water vapor effects and suppresses the diffusion of chromium and manganese which can prevent the generation of spinels on the alloy surface. Therefore, 310S, OC4 and OC5 at 800 °C and 625, OC11 and OC11LZ at 900 °C are selected to be tested in the Phase II long-term operation. From the results, we find that besides the lower Cr evaporation rates and better oxidation resistance of AFAs than benchmark alloys after short-term (500 h) operation, AFAs also possess the sturdy and compact alumina layer after a long-term operation (5000 h). A protective oxide layer is of great importance for the long-term high-temperature operation of structural materials. The high-temperature oxidation behavior of AFAs and commercial 310S and 625 alloys in Air + 10% H 2 O at 800 °C and 900 °C is systematically investigated. Severe breakaway oxidation and minor spallation are observed for the 310S and 625 after the long-term operation, while the AFAs show high oxidation resistance. It is found that the formation of a continuous alumina layer could greatly prevent the volatilization of chromium vapors. In addition, the corresponding evolving models are discussed. Chromium evaporation from BoP components in high-temperature environment could severely deteriorate the electrochemical performance of SOFC. Several methods were applied to evaluate the Cr evaporation rates of BoP after 500 h exposure at 800 °C to 900 °C in air with 10% H 2 O. An optimal method was designed to exclude the effect of silicon (Si) deposits from quartz tube and sodium (Na) deposits from the sodium carbonate on the oxidation process and the chemical interaction between Cr gaseous species and alumina tube which could provide further quantitative correlation of the evaporated Cr species quantities and degradation rates of SOFC. Based on the great performance of AFAs after long-term operation, AFAs are assembled with Anode-supported cells (ASC) to investigate the the Cr deposition of anode-supported cell under a constant current density of 0.5 A cm -2 at 800 °C with AFA alloys compared with commercial alloys. In addition, the anodic and cathodic processes are deconvoluted by distribution of relaxation times (DRT) method which are comprehensively discussed. It is found out that the voltage exhibited a slight decrease of 5.09 % and 1.54 % in the presence of OC11 and OC11LZA alloy, respectively. However, a considerable decrease of 22.14 % and 12.06 % was determined in the presence of 310S and 625 alloy, respectively. AFAs possessing the low-cost, low Cr evaporation rates and the high oxidation resistance will make it of great potential to replace the existing BoP components.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2020

Distribution Development for the RDX Regional Model at Los Alamos National Laboratory - 20385

Representing uncertainty in model inputs often means finding a balance between uncertainty and physical reality. Developing wide distributions may seem conservative in principle, but this approach may lead to unrealistic model results. Characterizing the current state of knowledge of stochastic inputs presents many challenges, especially if the data available are limited or have limited relevance to the site. Relationships among these inputs may also be important to represent but are typically complex or difficult to define. Often special adjustments must be made to account for reduced credibility in particular data. If parameters are strongly related to one another, a correlation structure may be developed for input to the model. Other techniques such as regression models may be used to incorporate relationships between the information available and the desired parameters. The process of developing distributions must consider details of the model in terms of what the distribution is meant to represent. This paper uses the example of a probabilistic fate and transport model for hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX) in the regional aquifer at Los Alamos National Laboratory (LANL). For many parameters, a single draw is applied to all space and time over which the model is run, for a single iteration. This simplification is often made for many reasons, and can often be beneficial, but also adds additional complexity in the distribution development process. Defining the distributional goals as they relate to the modeling process is an important step, which should take place prior to evaluation of the data. Usually, the distributions developed are meant to characterize the average value of the parameter over the spatial and temporal domain of the model. Distribution development requires consideration of many sources of information on the parameter where available, ideally from multiple references. Examples of different sources include data from different references but also from different conditions, measurement methods, or experimental types. Depending on these conditions and the reliability or relevance of particular references, different sources of data may each contribute valuable information but have varying relevance to the site. In these cases, weighting data unequally is a useful way to incorporate this information. As an example, aqueous dispersivity data are available for a variety of rock types. Only a few values are available for the desired rock type, and this is not enough to develop a distribution. Therefore, dispersivity values from other rock types are included in distribution development but are down-weighted such that the best data have the most influence on the distribution developed. In another case, K{sub d} distributions in the model are meant to represent a known composition of multiple soil types. Data from these materials are weighted accordingly to develop a distribution for the weighted average K{sub d} across soil types. Other cases include varying reliability of different sources, and weighting data according to the confidence in these sources. In some cases, input parameters are correlated with one another. An example is advective porosity, which is positively correlated with total porosity and must be less than total porosity. Paired data with both parameters must exist to discern the relationship between the two parameters and if it is necessary to build a correlation structure into the model. In general, correlated parameters may be represented in the model by a multivariate distribution, or perhaps more desirably, capturing the correlation within the developed distributions which may be treated as independent from one another. In the example of porosity, this can be done by transforming advective porosity into a proportion of total porosity which may be drawn independently from the distribution of total porosity. This paper explores how complex data and correlations can be incorporated to meet the distributional goals of the model, using the RDX regional model as a detailed example. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Degradation of Poly- and Perfluoroalkyl Substances (PFAS) in Water via High Power, Energy-Efficient Electron Beam Accelerator

The goal of the 2-year workplan was to see if electron beam (EB) could be used to break down a sub-set of the larger chemical family of per and polyfluoroalkylated substances (PFAS) in an energy efficient and economical manner when compared to conventional water treatment technologies. Year one (Y1) work focused on sample EB treatment work in the Fermi National Accelerator Laboratory’s (FNALs) Accelerator Applications Demonstration and Development (A2D2) EB accelerator. While there are reportedly thousands of types of PFAS, for the point of most of the work herein, a small subset was examined, typically perfluorooctane sulfonate (PFOS) and perfluorooctanoate (PFOA). PFOA and PFOS are two of the most well studied PFAS and are studied for baseline evaluations and are considered most useful. The work from Y1 provided information about the optimal operating parameters and additives to use when treating PFOS and PFOA via EB. The data were then used to see where in a water treatment system an EB accelerator would be best suited to treat PFAS. A conventional water treatment technology, GAC, was then compared to e-beam treatment technology with respect to energy and costs for treatment. In year two (Y2), several conventional e-beam accelerator designs, and FNAL’s developmental compact SRF accelerator design, were evaluated for their suitability in PFAS treatment, from an energy efficiency and cost standpoint. Several EB parameters were evaluated and optimized for the removal of PFOA and PFOS from water at normal pressure and temperature, measured as total PFAS removal. Under the optimized test conditions both PFOA showed complete destruction to inorganic fluoride, and PFOS to inorganic fluoride and sulfate, with mass balance. The effect on PFAS removal relative to solution pH, total EB dose, EB dose rate, dissolved oxygen concentration (DO), temperature, and initial PFAS concentration were evaluated. In general, PFOA was easier to destroy than PFOS. Degradation products, typically observed under less-than-optimal EB conditions, provided insight to degradation mechanisms. Products were identified to rule out possible deleterious biproduct formation. The water radiolysis radical reaction kinetics with PFOS and PFOA were not dependent on the initial concentration over 5-orders of magnitude from 2 μg/L to 20 mg/L. This is thought to be because there was an overabundance of the reactive water radiolysis radicals relative to PFAS molecules and largely attributed to aqueous electrons. The reaction rates appeared to be diffusion limited. Testing at higher concentrations (100-200 mg/L) showed a decrease in removal efficiency, suggesting alternative kinetics, possibly second order rates, at higher concentrations. In all, we successfully defined a set of optimal EB parameters to treat PFOA and PFOS at concentrations of 20 mg/L in water with destruction efficiencies near 100%. We further tested the optimized EB parameters with other types of PFAS, including shorter and longer fluorocarbon chain homologs of PFOA and PFOS, and PFAS with alternative functional groups such as sulfonamides. Based on our results EB can be optimized as an effective destructive technology for removing PFAS from water. The conditions optimized for PFOA and PFOS were less effective with ultra-short fluorocarbon compounds like TFMS, PFES, PFPS and PFBS, and likely require re-optimization of parameters to them. In all, it was determined that from a cost and energy efficiency standpoint, EB would be best applied to waste streams with relatively high concentrations of PFOS and PFOA and is not as cost effective as GAC treatment for removing low concentrations of PFAS from water. Higher concentrations of PFAS can be found in the wastewater of conventional treatment processes such as RO and IE and therefore EB may be used to supplement such treatment technologies. Some real-world IE regeneration wash water and RO reject water containing higher concentrations of PFAS and obtained from pilot scale industrial wastewater treatment system at a fluorochemical manufacturing facility, showed that EB could remove PFAS from such types of wastewaters. The IE regenerant wash water appeared to be the most efficient of the two types of wastewaters tested. However, some further optimization of the EB parameters for the specific PFAS types present in those wastewaters may be required. Also, the effects of co-present TOC and mineral salts should be considered during such optimization efforts. From the experimental Y1 results it was seen that the aqueous electron drives degradation of the PFAS. In a hypothetical water treatment skid using EB for PFAS destruction the parameters of the system should be optimized to promote aqueous electron production. Before EB treatment, the PFAS should be preconcentrated when possible, the pH should be raised to pH 10 or higher to enhance aqueous electron production, and the water should be nitrogen purged to remove dissolved oxygen to minimize aqueous electron scavenging. An excel spreadsheet was created that calculates optimal conditions based on inlet PFAS concentration and desired outlet concentration, by optimizing the accelerator power, dose rate, water treatment rate, pH and dissolved oxygen levels to reach the desired endpoint. Given this information on accelerator operating conditions five different EB accelerator systems were compared. One EB system was a continuous-wave, linear superconducting accelerator being designed at Fermilab. Three other EB systems (IMPELA at 5% and 25% duty factor and the ILU-14) were normal conducting pulsed linear accelerators. The fifth system was an IBA Rhodotron which is a normal conducting, circular, continuous-wave accelerator. The accelerator efficiency (% of the incoming power that is used in water treatment) was the dominating factor in accelerator choice. The radio frequency (RF) power supply and the accelerator design (superconducting versus warm technology) drive the accelerator efficiency. The IBA Rhodotron was seen to be the most energy efficient commercially available technology with a wall-plug (total) power efficiency of 43% at 400 kW. The Fermilab design, with a prototype for a different application currently being fabricated, was the most energy efficient at 55% when driven by a Klystron RF power supply and as high as 77% when powered by a magnetron. As the Fermilab design was the most energy efficient by approximately 10-30%, further design work was done on the accelerator and beam delivery system specific to the destruction of PFAS in water. The Fermilab design is unique from industrial accelerators in that is superconducting. Superconducting technology allows for the acceleration of electrons without losses. The accelerator must be cooled to below the point where it is superconducting and is operated around 4 degrees Kelvin. The bulk of the design work for the accelerator is on making the accelerator as energy efficient as possible so that it does not require liquid helium and can be cooled with conduction cooling via cryocoolers. Final design work resulted in an EB accelerator that would operate at minimally 200 kW and 10 MeV. Prototype construction would cost $\$ $7.8 million dollars when driven by a Klystron power supply. A second version of the same accelerator would cost $\$ $5.5 million dollars when driven by a magnetron that is still under development. The commercially available 300 kW IBA Rhodotron cost was estimated at approximately $\$ $9 million. While it is hard to directly compare, an operational GAC system used by 3M for groundwater treatment capital cost (2022 dollars) was estimated to cost $\$ $3.3 million. While the capital expense of the EB accelerator systems was higher than GAC, the accelerator EB treatment would result in destruction of the PFAS and not just sequestration of PFAS to form a new waste stream that requires further treatment or disposal. The operating cost to destroy the PFAS via 400 kw EB system was less than $\$ $1000/kg of PFAS destroyed when treating at a 20 mg/L PFAS concentration, compared to GAC with operating costs that calculated at $\$ $27,530 per kg of PFAS sequestered when treating 100 μg/L PFOA and PFOS combined concentration.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗