Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “explainable”

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

A Detailed Reaction Mechanism for Thiosulfate Oxidation by Ozone in Aqueous Environments

The ozone oxidation, or ozonation, of thiosulfate is an important reaction for wastewater processing, where it is used for remediation of mining effluents, and for studying aerosol chemistry, where its fast reaction rate makes it an excellent model reaction. Although thiosulfate ozonation has been studied since the 1950s, challenges remain in developing a realistic reaction mechanism that can satisfactorily account for all observed products with a sequence of elementary reaction steps. Here, we present novel measurements using trapped microdroplets to study the pH-dependent thiosulfate ozonation kinetics. We detect known products and intermediates, including SO32-, SO42-, S3O62-, and S4O62-, establishing agreement with the literature. However, we identify S2O42- as a new reaction intermediate and find that the currently accepted mechanism does not directly explain observed pH effects. Thus, we develop a new mechanism, which incorporates S2O42- as an intermediate and uses elementary steps to explain the pH dependence of thiosulfate ozonation. The proposed mechanism is tested using a kinetic model benchmarked to the experiments presented here, then compared to literature data. We demonstrate good agreement between the proposed thiosulfate ozonation mechanism and experiments, suggesting that the insights in this paper can be leveraged in wastewater treatment and in understanding potential climate impacts.

Deal, Alexandra M↗
DOE OSTI2025

Integrating Contaminant Source Indicators, Water Quality Measures, and Ecotoxicity to Characterize Contaminant Mixtures and Per- and Polyfluoroalkyl Substance (PFAS) Variability in an Urban Watershed

Thousands of chemical contaminants threaten watersheds but are time and cost prohibitive to monitor. Identifying their sources, transport, and ecological risk is limited in heterogeneous urban watersheds. We present an integrative watershed approach using source-specific indicator compounds, common water quality measures, and ecotoxicity assays to examine the distribution of contaminant mixtures in an urbanized watershed. Indicator compound concentrations were temporally and spatially distributed for treated/untreated sewage (sucralose, artificial sweetener), road runoff (diphenyl-guanidine [DPG] and 6PPD-quinone [6PPD-Q], automobile tire additives), and lawncare runoff (aminomethanephosphonic acid (AMPA), major degradant of the herbicide glyphosate). Sucralose was predominately sourced from treated wastewater; measurable concentrations in tributaries indicated raw sewage inputs. DPG and 6PPD-Q concentrations correlated to road density during base flow and were elevated during stormflow. AMPA was measurable spring through fall, especially where lawns were dense. When specific sources dominated flow, water quality measures correlated with wastewater (sulfate, potassium, chloride, and sodium) and road runoff (chromium and lead) indicators. The limited behavioral toxicity observed in exposed zebrafish (Danio rerio) (18%) was not well explained by source-indicators. PFAS concentrations were highly variable spatially but not well explained by our source-specific indicator compounds. Here, more costly compound-specific monitoring may be necessary when multiple sources exist or when unexpected toxicity trends occur.

computer simulations↗
DOE OSTI2020

Practical considerations for atomistic structure modeling with Cryo-EM maps

We present our experiences and opinions with selected software with respect to atomistic structure modeling with single particle analysis derived maps. These are atomistic model building and fitting method selection criteria and implementation procedures. Since we explain common concepts and caveats that we have explained to our collaborators, we hope that this short perspective will help facilitate active discussion between scientists at different levels and backgrounds.

Kim, Doo Nam↗
DOE OSTI2025

Energetic and Electronic Properties of AcX and LaX (X = O and F)

The bonding and spectroscopic properties of LaX and AcX (X = O and F) diatomic molecules were studied by highlevel ab initio CCSD(T) and SO-CASPT2 electronic structure calculations. Bond dissociation energies (BDEs) were calculated at the Feller−Peterson−Dixon (FPD) level. Potential energy curves and spectroscopic constants for the lowest-lying spin−orbit Ω states were obtained at the SO-CASPT2/aQ-DK level. A dense manifold of excited states was described for the monofluorides with the ground states well separated from the excited states. The spectroscopic parameters were in good agreement with those reported experimentally for LaO and LaF. For the diatomic molecules containing actinides, no experimental data of these parameters was found, but the results were consistent with other high-level calculations. The BDEs calculated at the FPD level were 791.3 (LaO), 705.2 (AcO), 650.0 (LaF), and 678.6 (AcF) kJ/ mol. The NBO analysis showed that the monofluorides are essentially ionic, which explains why the BDE(AcF) is higher than BDE(LaF); for the monoxides, covalent contributions involving the d orbitals of the metal and the p orbitals of the oxygen are stronger for LaO than AcO, which explains the higher BDE for LaO. Finally, the bond orders are predicted to be 2 for LaF and AcF, 3 for AcO, and higher than 3 for LaO.

Basis sets↗
DOE OSTI2025

Modeling Pb(II) Adsorption on Mineral Surfaces: Bridging Density Functional Theory and Experiment with Thermodynamic Insights

Despite decades of work on aqueous lead (Pb) adsorption on a-Fe2O3 (hematite) and a-Al2O3 (alumina), gaps between measurements and modeling obscure molecular-level understanding. Achieving well-matched geometries between theory and experimental for mineral-water interfaces is a hurdle, as surface functional group type and distribution must be accounted for in determining mechanisms. Additionally, computational methods that can describe the substrate are often not appropriate to capture aqueous effects. Progress requires focusing on well-studied and relevant systems, such as key facets (001),(012), and (110) of hematite and alumina, and ubiquitous contaminants such as aqueous Pb. In the past, bulk-parametrized bond-valence principles were used to rationalize Pb(II) adsorption trends. These approaches can break down at surfaces, where flexible bonding environments and adsorption-induced surface relaxations play a critical role. Here, we adapt and apply a density functional theory (DFT) and thermodynamics framework, integrating DFT-calculated energies with experimental data and electrochemical principles, to predict Pb(II) adsorption. Our model results capture trends across the full set of surfaces and predict that inner-sphere Pb(II) sorption on (001) alumina varies from unfavorable to weakly favorable across a range of pH conditions. This aligns with experiment insights that Pb(II) interacts at that surface through outer-sphere interactions. Extending to Fe(II) adsorption, we demonstrate a coverage-dependent site preference, potentially explaining disorder in overlayers grown by the oxidative adsorption of Fe(II) on hematite (001).

lead contamination↗
DOE OSTI2020

Surface Enhanced Raman Scattering Selectivity in Proteins Arises from Electron Capture and Resonant Enhancement of Radical Species

Plasmon-enhanced Raman scattering is a powerful approach to detecting and characterizing proteins in live and dynamic biological systems. However, the selective detection/enhancement of specific residues as well as spectral diffusion and fluctuations have complicated the interpretation of enhanced Raman spectra and images of biological matter. In this paper, we therefore demonstrate that the amino acid tryptophan (Trp) can capture an electron from an excited plasmon, which generates a radical anion that is resonantly enhanced: a visible excited electronic state slides into resonance upon charging. This surface enhanced resonance Raman scattering (SERRS) mechanism explains the persistence of Trp signatures in the SERS and TERS spectra of proteins. Evidence for this picture includes the observation of visible resonances in the UV-Vis extinction spectrum, changes in the ground state vibrational spectrum, and plasmon-resonance dependent behavior. DFT calculations support the experimental observations. The behavior observed from the free Trp molecule is shown to explain the SERS spectrum of the Trp-cage protein. In effect, resonant Raman scattering from radicals formed through plasmonic excitation represents an under-investigated mechanism that may be exploited for chemical sensing applications. This work was supported by the National Institutes of Health Award R01-GM109988. The authors thank Hannah Shafaat for helpful discussions and providing an azurin sample for comparison. The authors also thank Joel Harris for providing constructive comments on our manuscript.

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

Understanding Fragmentation of Organic Small Molecules in Atom Probe Tomography

In atom probe tomography of molecular organic materials, field ionization of either entire molecules or molecular fragments can occur, but the mechanism governing this behavior was not previously understood. This work explains when a doubly ionized small molecule organic material is expected to undergo fragmentation. Here, we find that multiple detection events arising from post-ionization fragmentation of a parent molecular dication into two daughter ions is well explained by the free energy and geometries of the molecules computed using density functional theory. Of the systems studied, exergonic free energies for formation of the daughter ions, smaller activation energies for dissociation, and increases in bond length are all found to be quantitative predictors for ion fragmentation. This work expands the applicability of atom probe tomography to organic materials by increasing the fundamental understanding of processes occurring during this analysis technique.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Role of Local Conductivities in the Plasmon Reflections at the Edges and Stacking Domain Boundaries of Trilayer Graphene

Here, we employed infrared scattering-type scanning near-field optical microscopy (IR-sSNOM) to study surface plasmon polaritons (SPPs) in trilayer graphene (TLG). Our study reveals systematic differences in near-field IR spectra and SPP wavelengths between Bernal (ABA) and rhombohedral (ABC) TLG domains on SiO 2 , which can be explained by stacking-dependent intraband conductivities. We also observed that the SPP reflection profiles at ABA-ABC boundaries could be mostly accounted for by an idealized domain boundary defined by the conductivity discontinuity. However, we identified distinct shapes in the SPP profiles at the edges of the ABA and ABC TLG, which cannot be solely attributed to idealized edges with stacking-dependent conductivities. Instead, this can be explained by the presence of various edge structures with local conductivities differing from those of bulk TLGs. Our findings unveil a new structural element that can control SPP and provide insights into the structures and electronic states of the edges of few-layer graphene.

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

Wetting Behavior of A -block- (B- random -C) Copolymers with Equal Block Surface Energies on Surfaces Functionalized with B- random -C Copolymers

To form nanopatterns with self-assembled block copolymers (BCPs), it is desirable to have through-film domains that are oriented perpendicular to the substrate. The domain orientation is determined by the interfacial interactions of the BCP domains with the substrate and with the free surface. Here, we use thin films of two different sets of BCPs with A-block-(B-random-C) architecture matched with a corresponding B-random-C copolymer nanocoating on the substrate to demonstrate two distinct wetting behaviors. The two sets of A-b-(B-r-C) BCPs are made by using thiol–epoxy click chemistry to functionalize polystyrene-block-poly(glycidyl methacrylate) with trifluoroethanethiol (TFET) and either 2-mercaptopyridine (2MP) or methyl thioglycolate (MTG). For each set of BCPs, the composition ratio of the two thiols in the BCP (φ 1 ) is found that results in the two blocks of the modified BCP having equal surface energies (Δγ air = 0). The corresponding B-r-C random copolymers were synthesized and used to modify the substrate, and the composition ratio (φ 2 ) values that resulted in the two blocks of the BCP having equal interfacial energy with the substrate (Δγ sub = 0) were determined with scanning electron microscopy. The correlation between each block’s γ sub value and the interaction parameter, χ, is employed to explain the different wetting behaviors of the two sets of BCPs. For the thiol pair 2MP and TFET, the values of φ 1 and φ 2 that lead to Δγ air = 0 and Δγ sub = 0, respectively, are significantly different. A similar difference was observed between the φ 1 and φ 2 values that lead to Δγ air = 0 and Δγ sub = 0 for the BCPs made with the thiol pair MTG and TFET. In the latter case, for Δγ sub = 0 two windows of φ 2 are identified, which can be explained by the thermodynamic interactions of the specific thiol pair and the A-b-(B-r-C) architecture.

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

Nonlinear Shear Rheology of Entangled Polymer Rings

Steady-state shear viscosity η($\bar{γ}$) of unconcatenated ring polymer melts as a function of the shear rate $\bar{γ}$ is studied by a combination of experiments, simulations, and theory. Experiments using polystyrenes with Z ≈ 5 and Z ≈ 11 entanglements indicate weaker shear thinning for rings compared to linear polymers exhibiting power law scaling of shear viscosity η ~ $\bar{γ}$ –0.56 ± 0.02 , independent of chain length, for Weissenberg numbers up to about 10 2 . Nonequilibrium molecular dynamics simulations using the bead-spring model reveal a similar behavior with η ~ $\bar{γ}$ -0.57 ± 0.08 for 4 ≤ Z ≤ 57. Viscosity decreases with chain length for high $\bar{γ}$. In our experiments, we see the onset of this regime, and in simulations, which we extended to Wi ~ 10 4 , the nonuniversality is fully developed. In addition to a naive scaling theory yielding for the universal regime η ~ $\bar{γ}$ –0.57 , we developed a novel shear slit model explaining many details of observed conformations and dynamics as well as the chain length-dependent behavior of viscosity at large $\bar{γ}$. The signature feature of the model is the presence of two distinct length scales: the size of tension blobs and much larger thickness of a shear slit in which rings are self-consistently confined in the velocity gradient direction and which is dictated by the size of a chain section with relaxation time 1/$\bar{γ}$. These two length scales control the two normal stress differences. In this model, the chain length-dependent onset of nonuniversal behavior is set by tension blobs becoming as small as about one Kuhn segment. This model explains the approximate applicability of the Cox–Merz rule for ring polymers.

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

Entropic Penalty Switches Li + Solvation Site Formation and Transport Mechanisms in Mixed Polarity Copolymer Electrolytes

Emerging solid polymer electrolyte (SPE) designs for efficient Li-ion (Li + ) conduction have relied on polarity and mobility contrast to improve conductivity. To further develop this concept, we employ simulations to examine Li + solvation and transport in poly(oligo ethylene methacrylate) (POEM) and its copolymers with poly(glycerol carbonate methacrylate) (PGCMA). We find that Li + is solvated by ether oxygens instead of the highly polar PGCMA, due to lower entropic penalties. The presence of PGCMA promotes single-chain solvation, thereby suppressing interchain Li + hopping. The conductivity difference between random copolymer PGCMA-r-POEM and block copolymer PGCMA-b-POEM is explained in terms of a hybrid solvation site mechanism. With diffuse microscopic interfaces between domains, PGCMA near the POEM contributes to Li + transport by forming hybrid solvation sites. The formation of such sites is hindered when PGCMA is locally concentrated. These findings help explain how thermodynamic driving forces govern Li + solvation and transport in mixed SPEs.

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

Mechanistic Advantages of Organotin Molecular EUV Photoresists

Here, extreme ultraviolet (EUV)-induced radiation exposure chemistry in organotin-oxo systems, represented by the archetypal [(R-Sn) 12 O 14 (OH) 6 ](A) 2 cage, has been investigated with density functional theory. Upholding existing experimental evidence of Sn-C cleavage-dominant chemistry, computations have revealed that either electron attachment or ionization can single-handedly trigger tin-carbon bond cleavage, partially explaining the current EUV sensitivity advantage of metal oxide systems. We have revealed that tin atoms at different parts of the molecule react differently to ionization and electron attachment and have identified such selectivity as a result of local coordination chemistry instead of the macro geometry of the molecule. An ionization-deprotonation pathway has also been identified to explain the observed evolution of an anion conjugate acid upon exposure and anion mass dependence in resist sensitivity.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2020

The Role of Protons and Hydrides in the Catalytic Hydrogenolysis of Guaiacol at the Ruthenium Nanoparticle-water Interface

The mechanistic roles of free hydronium ions, surface hydrides, and interfacial protons during guaiacol hydrodeoxygenation (HDO) on ruthenium nanoparticles are established. As guaiacol adsorbs on Ru, it loses its strong aromaticity and undergoes a rapid H-shift from its hydroxyl to meta carbons (in relation to its hydroxyl group), leading enol and keto surface isomers to exist in chemical equilibrium. HDO occurs via a hydridic H-adatom (H*) attack to the enol, followed by a kinetically relevant C-O bond rupture step, during which water shuttles the hydroxyl proton, enabling its intramolecular attack to the methoxy, evolving a high charged [Ru(s)-(C6H5O-)…(H+)…OCH3]† transition state. The competing HYD begins with a rapid H* attack to the keto, before a second, kinetically relevant H* attack, without proton involvement. Water, despite shifting the thermodynamics towards the more polar surface keto, promotes HDO to a much greater extent than HYD, because of its dual catalytic roles—it mobilizes hydroxyl proton (Brønsted acid) to cleave the strong C-O bond, synchronizing with the Ru metal surface (base) function that stabilizes the resulting [Ru(s)-(C6H5O-)…(H+)…OCH3]† transition state, and the water layers solvate this charged transition state, further reducing its free energy. Free hydronium ions do catalyze a separate homogeneous enol-keto isomerization, but this reaction is kinetically unrelated to HDO catalysis. This mechanistic picture explains the strong effects of polar protic solvent in hydrodeoxygenation, highlighting (i) the requirements of surface hydrides and interfacial protons acting in tandem to complete a HDO turnover and (ii) the cooperative role of protic solvent and metal surface in breaking the aromaticity and stabilizing charged reactive precursors and transition states.

Shangguan, Junnan↗
DOE OSTI2021

Unraveling the Intermediate Reaction Complexes and Critical Role of Support-Derived Oxygen Atoms in CO Oxidation on Single-Atom Pt/CeO 2

CeO 2 -supported Pt single-atom catalysts have been extensively studied due to their relevance in automobile emission control and for the fundamental understanding of CeO 2 -based catalysts. Though CeO 2 -supported Pt nanoparticles are often more active than their single-atom counterparts, the former could easily redisperse to Pt single atom under oxidizing diesel conditions. Therefore, to maximize the reactivity of every Pt atom, it is important to fully understand the reaction mechanism of CeO 2 -supported Pt single atoms. Here, we report a CO oxidation study on a Pt/CeO 2 single-atom catalyst, where we can account for all of the neighbors using in situ and operando spectroscopy techniques and microcalorimetric measurements. Coupled with density functional theory calculations, we present a comprehensive picture of the dynamics of the surface species, the role of surface intermediates, and explain the observed reaction kinetics. We started with a catalyst containing exclusively single atoms and used in situ/operando spectroscopy to provide evidence for their stability during the reaction and to identify the Pt 1 complexes before and during the reaction and their binding to CO. The results reveal that in the precatalyst, Pt is present as Pt(O) 4 on the CeO 2 (111) step edge sites, but during CO oxidation, we find that two Pt 1 complexes coexist, representing two states of the same active site in the reaction cycle. The dominant state/complex remains Pt(O) 4 , which adsorbs CO very weakly as shown by CO microcalorimetry. The second, minority state/complex, Pt(CO)(O) 3 is generated through the reaction of Pt(O) 4 with CO, and CO is bound strongly to Pt 1 . Labile oxygen adatoms from the CeO 2 surface play a major role in the regeneration of Pt(O) 4 either directly from Pt(O) 3 or by reaction with the strongly adsorbed CO in Pt(CO)(O) 3 . We show that the formation of an oxygen vacancy and generation of a labile O* are not barrierless, which explains the long lifetime of Pt(CO)(O) 3 and its detectability despite being a minority complex. The results help to develop a comprehensive view of the dynamic evolution of Pt 1 complexes along the reaction cycle and provide mechanistic insights to guide the design of Pt-based single-atom catalysts.

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

Electrochemical CO2 Reduction over Metal-/Nitrogen-Doped Graphene Single-Atom Catalysts Modeled Using the Grand-Canonical Density Functional Theory

Renewably driven, electrochemical conversion of carbon dioxide into value-added products is expected to be a critical tool in global decarbonization. However, theoretical studies based on the computational hydrogen electrode largely ignore the nonlinear effects of the applied potential on the calculated results, leading to inaccurate predictions of catalytic behavior or mechanistic pathways. Here, we use grand canonical density functional theory (GC-DFT) to model electrochemical CO2 reduction (CO2R) over metal- and nitrogen-doped graphene catalysts (MNCs) and explicitly include the effects of the applied potential. We used GC-DFT to compute the CO2 to CO reaction intermediate energies at -0.3, -0.7, and -1.2 VSHE catalyzed by MNCs each doped with 1 of the 10 3d block metals coordinated by four pyridinic nitrogen atoms. Our results predict that Sc-, Ti-, Co-, Cu-, and Zn-N4Cs effectively catalyze CO2R at moderate to large reducing potentials (-0.7 to -1.2 VSHE). ZnN4C is a particularly promising electrocatalyst for CO2R to CO both at low and moderate applied potentials based on our thermodynamic analysis. Our findings also explain the observed pH independence of CO production over FeN4C and predict that the rate-determining step of CO2R over FeN4C is not *CO2- formation but rather *CO desorption. Additionally, the GC-DFT-computed density of states analysis illustrates how the electronic states of MNCs and adsorbates change non-uniformly with applied potential, resulting in a significantly increased *CO2- stability relative to other intermediates and demonstrating that the formation of the adsorbed *CO2- anion is critical to CO2R activation. This work demonstrates how GC-DFT paves the way for physically realistic and accurate theoretical simulations of reacting electrochemical systems.

CO2 reduction↗
DOE OSTI2022

Hydrogen Dissociation Controls 1-Hexyne Selective Hydrogenation on Dilute Pd-in-Au Catalysts

Increasing the selectivity of the catalytic hydrogenation of alkynes to alkenes is of major importance for the processing of petrochemicals and the production of fine chemicals. Achieving high selectivity for alkene formation at high conversions, however, remains a long-standing challenge in heterogeneous catalysis. Here, the mechanism and origin of the high selectivity of dilute Pd-in-Au catalysts has been studied by a combination of first-principles calculations, microkinetic simulations, and isotopic exchange hydrogenation experiments using Pd 0.04 Au 0.96 nanoparticles embedded in raspberry colloid-templated silica. The Pd is predominantly in the form of isolated atoms, only surrounded by Au atoms, based on prior studies. The simulations indicate that the rate-limiting process for 1-hexyne hydrogenation on Pd monomers in Au(111) is H 2 dissociation, which has a large free energy barrier of 0.86 eV at 363 K and 0.2 bar of H 2 . The C–H bond formation steps, on the other hand, proceed with lower barriers, which contrasts with previous studies of extended Pd catalysts. The microkinetic simulations identify the sizable H 2 dissociation barrier and the small barrier for the hydrogenation of 1-hexyne as key factors that lead to a high selectivity for the production of 1-hexene from 1-hexyne, even at high conversion. The unconventional H 2 dissociation limiting process in combination with the low coverage of weakly bound hydrocarbon intermediates explains the near-zero order of 1-hexyne found experimentally. Furthermore, the partial hydrogenation of 1-hexyne to form 1-hexene is shown to be an irreversible process from our isotopic exchange hydrogenation experiments and is explained by the strongly exothermic nature of the reaction. Diluting active species, like Pd, in a less active host metal, like Au, hence appears promising as a means of tuning the binding energy of reactants and altering reaction profiles, leading to distinct kinetic behavior for an optimal catalytic activity and selectivity. Furthermore, the combination of microkinetic modeling, density functional theory calculations, and isotopic exchange experiments is thus demonstrated to be an effective approach to modeling important catalytic phenomena.

1-hexyne↗
DOE OSTI2026

Ripening of Rh Nanoparticle Catalysts in Reverse Water–Gas Shift via a Data-Driven Model Combining Physics, Theory, and Experiment

Degradation via sintering is an ongoing challenge that impedes the broad commercial success of supported metallic nanoparticle catalysts. To mitigate degradation via informed catalyst design and process operations, here we aim to disambiguate the underlying mechanisms of sintering by combining theory and experiment in a quantitative framework. While mechanistic sintering models exist, they only model a single sintering pathway, even though multiple sintering mechanisms can occur simultaneously or dominate at different stages of the process. Data-driven machine learning models have emerged as a means to represent complex processes through data regression. However, machine learning models have very large data needs and lack mechanistic insights due to their black-box encoding. To develop an interpretive model of catalyst degradation via sintering, we constructed a hybrid model combining mechanistic “physics-based” models and data-driven methods to obtain both reliable predictions and mechanistic insights regarding experimentally observed sintering phenomena. Focusing on nanoparticle sintering in the Rh–TiO 2 catalyst for the reverse water–gas shift (RWGS) reaction, the hybrid model couples a mechanistic term for Ostwald ripening with energy values calculated via density functional theory (DFT) with a parametric, data-driven discrepancy function term for unmodeled mechanisms. The hybrid model is trained using Bayesian inference with data collected from small-angle X-ray scattering (SAXS) in situ experiments wherein average nanoparticle diameter versus time was measured at three relevant operating temperatures. The calibrated hybrid model results show that an Ostwald ripening-only model parameterized with fixed DFT energies does not fully capture the time and temperature dependence of the SAXS-observed sintering kinetics, and that an additional functional contribution, or DFT energy calibration, is required to reconcile simulation and experiment. Analysis of the hybrid-model error confirms that the hybrid model outperforms both the purely mechanistic and purely data-driven alternatives in terms of expected predictive accuracy for time-evolving average particle sizes. Furthermore, the results support the hypothesis that the Ostwald ripening mechanism is less important for explaining the sintering phenomena as operating temperature increases under an assumed fixed DFT parameterization. This could be explained in one of two ways: either latent, unmodeled sintering mechanisms dominate at higher temperatures, or the DFT uncertainty increases with temperature. The proposed modeling approach directly links theory to experiments and simulations via a statistical hybrid modeling framework and can be extended to other catalytic systems to improve predictive models and mechanistic understanding.

Bayesian hybrid modeling↗
DOE OSTI2021

Loss and Isotopic Fractionation of Alkali Elements during Diffusion-Limited Evaporation from Molten Silicate: Theory and Experiments

Moderately volatile elements (MVEs) are variably depleted in planetary bodies, reflecting the imprints of nebular and planetary processes. Among MVEs, Na, K, and Rb are excellent tracers for unraveling the history of MVE depletion in planetary bodies because they have similar geochemical behaviors but can be chemically fractionated by evaporation and condensation processes. Furthermore, K and Rb are amenable to high-precision isotopic analyses, which can help constrain the conditions of evaporation and condensation. To quantitatively understand why Na, K, and Rb are depleted in planetary bodies, we have carried out vacuum evaporation experiments from basaltic melt at 1200 and 1400 °C to study their evaporation kinetics and isotopic fractionations. We chose this composition because it is relevant to evaporation from small differentiated planetesimals. The Rb isotopic compositions of the evaporation residues were measured by multicollector inductively coupled plasma mass spectrometry (MC-ICPMS), and the K isotopic compositions were measured along profiles across the residues by secondary ion mass spectrometry (SIMS). In the 1400 °C run products, we found that the concentrations of both K and Rb in the run products decreased from core to rim, which was accompanied by a heavy K isotope enrichment near the surface. This indicates that, in this run, evaporation was limited by diffusion. To use those data quantitatively, we derive analytical equations that describe the evaporation rate and isotopic fractionation associated with diffusion-limited evaporation from a sphere, slab, and cylinder in transient and quasi-steady state regimes. This model is used to tease out the roles that diffusive transport in the melt and evaporation at the melt/gas interface play in setting the elemental depletion and isotopic composition of the residue. Under our experimental conditions, volatility decreases in the order of Na, Rb, and K. Using our experimental results in a thermodynamic model, we have estimated the product γΓ of activity coefficients × evaporation coefficients of Na, Rb, and K. The measured isotopic compositions of the residues are well explained using Rayleigh distillations, whereby the relative volatilities of K and Rb isotopes are given by the square root of their masses. We use our results and previously published data to predict how K and Rb could have been lost as a function of temperature, melt composition, oxygen fugacity, and saturation degree relevant to Vesta’s building blocks. We find that the K and Rb depletions, K/Rb elemental fractionation, and δ 41 K and δ 87 Rb isotopic fractionations of Vesta (as sampled by howardite-eucrite-diogenite (HED) meteorites) are best explained by evaporation of submillimeter size objects for 0.1-10 years at moderate temperatures (~1050 °C) in a medium ~98.8% saturated.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗