Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reacting”

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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Influence of Nd(II) and Nd(III) Ions on the Speciation and Kinetic Dynamics of Radiolytic Transients in Molten LiCl-KCl Eutectic Salt Mixtures

In the pyrochemical reprocessing of used nuclear fuel, understanding the effects of ionizing radiation on the chemical behavior of fission products in molten chloride eutectic salt is crucial. Here, this study investigates the impact of radiation on the chemical and kinetic speciation of Nd ions in molten LiCl-KCl eutectic, which are challenging to separate due to their multivalent oxidation states and comproportionation reactions. Electron pulse radiolysis techniques were used to determine transient absorption spectra, chemical kinetics, and activation parameters for the reaction of Nd ions with the primary products of molten LiCl-KCl eutectic radiolysis: the solvated electron (e S - ) and dichlorine radical anion (Cl 2 •- ). Initially, Nd(III) reacted rapidly with e S - , forming Nd(II) with a second-order rate coefficient of k = (4.54 ± 0.07) × 10 10 M -1 s -1 at 673 K. The resulting Nd(II) ions then reacted slowly with Cl 2 •- , regenerating Nd(III) with an estimated rate of k = (1.72 ± 0.04) × 10 10 M -1 s -1 at 673 K. Additionally, the comproportionation reaction of Nd(III) and Nd(s) to form Nd(II) was monitored chemically and electrochemically. This work suggests that the comproportionation reaction of Nd(III) with Nd(s) to form Nd(II) is unstable in molten LiCl-KCl eutectic.

36 - MATERIALS SCIENCE↗
DOE OSTI2021

CyProduct: A software tool for accurately predicting the byproducts of human cytochrome P450 metabolism

In silico metabolism prediction is a cheminformatic task of autonomously predicting the set of metabolic byproducts produced from a specified molecule and a set of enzymes or reactions. Here we describe a novel machine-learned in silico cytochrome P450 (CYP450) metabolism prediction suite, called CyProduct, that accurately predicts metabolic byproducts for a specified molecule and a human CYP450 isoform. It includes three modules: (1) CypReact, a tool that predicts if the query compound reacts with a given CYP450 enzyme; (2) CypBoM, a tool that accurately predicts the “bond site” of the reaction (i.e., which specific bonds within the query molecule react with the CYP isoform); and (3) MetaboGen, a tool that generates the metabolic byproducts based on CypBoM’s bond-site prediction. CyProduct predicts metabolic biotransformation products for each of the nine most important human CYP450 enzymes. CypBoM uses an important new concept called “Bond of Metabolism” (BoM), which extends the traditional “Site of Metabolism" (SoM) by specifying the information about the set of chemical bonds that is modified or formed in a metabolic reaction (rather than the specific atom). We created a BoM database for 3487 CYP450-mediated Phase I reactions, then used this to train the CypBoM Predictor to predict the reactive bond locations on substrate molecules. CypBoM Predictor’s cross-validated Jaccard score for reactive bond prediction ranged from 0.380 to 0.452 over the nine CYP450 enzymes. Over variants of a test set of 72 known CYP450 substrates and 30 non-reactants, CyProduct outperformed the other packages -- including ADMET Predictor, BioTransformer and GLORY -- by an average of 200% (wrt Jaccard score) in terms of predicting metabolites. The CyProduct suite and the datasets are freely available at https://bitbucket.org/wishartlab/cyproduct/src/master/.

Machine learning, Cytochrome P450, Metabolism pred↗
DOE OSTI2024

Evaluating the Chemical Reactivity of DFT-Simulated Liquid Water with Hydrated Electrons via the Dual Descriptor

Modeling the various properties of liquid water, particularly its reactivity, has been a longstanding challenge for simulation methods. Recently, ab initio simulations based on density functional theory (DFT) have come to the fore as tenable methods for calculating the properties and reactivity of water, with varying degrees of success for different exchange-correlation functionals. In particular, hybrid-GGA and meta-GGA functionals have been shown to reproduce many of the structural, dynamical, and energetic properties of water to a high degree of accuracy relative to their computational cost. Here, we show that the dual descriptor (DD) measure of nucleophilicity and electrophilicity, which is sometimes used to elucidate organic chemistry reaction mechanisms, can also be used to characterize the reactivity of DFTsimulated liquid water. The DD is especially apt for understanding the reactivity of excess electrons with water as its calculation explicitly involves adding and removing an excess electron from a reference system. We use the DD to explore the reactivity of water simulated using three different DFT functionals: the LDA functional (LDA), a hybrid-GGA functional (PBE0), and a hybrid meta- GGA functional (SCAN0). Using the DD, we show that the SCAN0 functional with the standard 25% Hartree–Fock exchange produces simulated liquid water with many regions that are far more reactive than either PBE0 or LDA. To understand the implications of these highly reactive regions, we then add a strong nucleophile in the form of an excess electron and find that although PBE0 and LDA predict stable hydrated electrons, the excess electron reacts nearly instantaneously with SCAN0 water via proton abstraction to form a hydrogen atom and hydroxide ion. We show that the DD provides the ability to not only predict whether or not liquid water will react with a hydrated electron but also which particular waters will be involved solely from analyzing pure water configurations generated with each functional. We rationalize this result in terms of the known trap-seeking behavior of injected hydrated electrons, which are able to find the most electronegative region in bulk water. These results highlight the utility of the dual descriptor as a fast and interpretable method for investigating condensed-phase reactivity with excess electrons.

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

Coupled Interfacial and Bulk Kinetics Govern the Timescales of Multiphase Ozonolysis Reactions

Chemical transformations in aerosols impact the lifetime of particle phase species, the fate of atmospheric pollutants, and both climate- and health-relevant aerosol properties. Timescales for multiphase reactions of ozone in atmospheric aqueous phases are governed by coupled kinetic processes between the gas phase, the particle interface, and its bulk, which respond dynamically to reactive consumption of O 3 . However, models of atmospheric aerosol reactivity often do not account for the coupled nature of multiphase processes. Herein, to examine these dynamics, we use new and prior experimental observations of aqueous droplet reaction kinetics, including three systems with a range of surface affinities and ozonolysis rate coefficients ( trans -aconitic acid (C 6 H 6 O 6 ), maleic acid (C 4 H 4 O 4 ), and sodium nitrite (NaNO 2 )). Using literature rate coefficients and thermodynamic properties, we constrain a simple two-compartment stochastic kinetic model which resolves the interface from the particle bulk and represents O 3 partitioning, diffusion, and reaction as a coupled kinetic system. Our kinetic model accurately predicts decay kinetics across all three systems, demonstrating that both the thermodynamic properties of O 3 and the coupled kinetic and diffusion processes are key to making accurate predictions. An enhanced concentration of adsorbed O 3 , compared to gas and bulk phases is rapidly maintained and remains constant even as O 3 is consumed by reaction. Multiphase systems dynamically seek to achieve equilibrium in response to reactive O 3 loss, but this is hampered at solute concentrations relevant to aqueous aerosol by the rate of O 3 arrival in the bulk by diffusion. As a result, bulk-phase O 3 becomes depleted from its Henry's law solubility. This bulk-phase O 3 depletion limits reaction timescales for relatively slow-reacting organic solutes with low interfacial affinity (i.e., trans-aconitic and maleic acids, with k rxn ≈ 10 3 -10 4 M -1 s -1 ), which is in contrast to fast-reacting solutes with higher surface affinity (i.e., nitrite, with k rxn ≈ 10 5 M -1 s -1 ) where surface reactions strongly impact the observed decay kinetics.

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

Diffusion-Limited Kinetics in Reactive Systems

A proper representation of chemical kinetics is vital to understanding, modeling, and optimizing many important chemical processes. In liquid and surface phases, where diffusion is slow, the rate at which the reactants diffuse together limits the overall rate of many elementary reactions. Commonly, the textbook Smoluchowski theory is utilized to estimate effective rate coefficients in the liquid phase. On surfaces, modelers commonly resort to much more complex and expensive Kinetic Monte Carlo (KMC) simulations. Here, in this study, we extend the Smoluchowski model to allow the diffusing species to undergo chemical reactions and derive analytical formulas for the diffusion-limited rate coefficients for 3D, 2D, and 2D/3D interface cases. With these equations, we are able to demonstrate that when species react faster than they diffuse they can react orders of magnitude faster than predicted by Smoluchowski theory, through what we term “the reactive transport effect”. We validate the derived steady-state equations against particle Monte Carlo (PMC) simulations, KMC simulations, and non-steady-state solutions. Furthermore, using PMC and KMC simulations, we propose corrections that agree with all limits and the computed data for the 2D and 2D/3D interface steady-state equations, accounting for unique limitations in the associated derived equations. Additionally, we derive equations to handle couplings between diffusion-limited rate coefficients in reaction networks. We believe these equations should make it possible to run much more accurate mean-field simulations of liquids, surfaces, and liquid–surface interfaces accounting for diffusion limitations and the reactive transport effect.

Johnson, Matthew Sean↗
DOE OSTI2020

Investigating the Chemical Reactivity of Lithium Silicate Model SEI Layers

Silicon anodes suffer from an unstable solid electrolyte interphase (SEI) layer that contributes to undesirable capacity fade with cycling. A key part to addressing this unstable SEI formation is to examine how certain components of the SEI react with the electrolyte over time. One SEI component that has not been thoroughly studied in the context of the chemical reactivity against the electrolyte is lithiated silicate. Four model silicate thin films with increasing lithium content were deposited by radio frequency (RF) magnetron sputtering to study how the lithiation of the native oxide on a silicon anode affects the chemical stability of the anode surface. SiO2, Li2Si2O5, Li2SiO3, and Li3SiOx films were exposed to 1.2 M LiPF6 in the 3:7 wt% ethylene carbonate/ethyl methyl carbonate (EC/EMC) electrolyte for periods of time that are representative of the amount of time it takes to undergo cell formations. Soaked samples were rinsed, dried, and characterized by a combination of attenuated total reflectance-infrared spectroscopy (ATR-IR), focused ion beam-secondary electron microscopy (FIB-SEM), and X-ray photoelectron spectroscopy (XPS) depth profiling. It was found that the rate of the decrease in film thickness of the silicates exposed to the electrolyte over time increases as a function of the lithium content in the thin film. This reaction involves HF etching and LiPF6 salt degradation leading to silicate loss and fluorination throughout the bulk. Understanding this chemical instability is critical to determining the overall mechanism of SEI degradation over time.

ADVANCED PROPULSION SYSTEMS,ENERGY STORAGE↗
DOE OSTI2020

An Unexpected Role of H During SiC Corrosion in Water

During aqueous corrosion, atoms in the solid react chemically with oxygen, leading either to the formation of an oxide film or to the dissolution of the host material. Commonly, the first step in corrosion involves an oxygen atom from the dissociated water that reacts with the surface atoms and breaks near-surface bonds. In contrast, hydrogen on the surface often functions as a passivating species. Here, we discovered that the roles of O and H are reversed in the early corrosion stages on a Si-terminated SiC surface. O forms stable species on the surface, and chemical attack occurs by H that breaks the Si–C bonds. Here, this so-called hydrogen scission reaction is enabled by a newly discovered metastable bridging hydroxyl group that can form during water dissociation. The Si atom that is displaced from the surface during water attack subsequently forms H 2 SiO 3 , which is a known precursor to the formation of silica and silicic acid. This study suggests that the roles of H and O in oxidation need to be reconsidered.

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

Closing the Gap Between Experiment and Theory: Reactive Scattering of HCl from Au(111)

Accurate simulation of molecules reacting on metal surfaces, which can help in improving heterogeneous catalysts, remains out of reach for several reactions. For example, a large disagreement between theory and experiment for HCl reacting on Au(111) still remains, despite many efforts. In this work, the dissociative chemisorption of HCl on Au(111) is investigated with a recently developed MGGA density functional (MS-RPBEl) and a high-dimensional neural network potential. Additionally, previous experimental sticking probabilities are re-examined. A considerably improved agreement between experiment and theory is obtained, although theory still overestimates experimental sticking probabilities by a factor of 2-7 at the highest incidence energy. Computed and measured vibrational transition probabilities are also in improved agreement. Several dynamical effects such as angular steering and energy transfer from the molecule to the surface are found to play an important role.

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

Detecting Reactive Products in Carbon Capture Polymers with Chemical Shift Anisotropy and Machine Learning

Aminopolymers are attractive sorbents for CO 2 direct air capture applications due to their high density of amine groups, which can readily react with atmospheric levels of CO 2 to form chemisorbed species. The identity of these chemisorbed species and the functional groups that form upon oxidative degradation depends on both material properties and processing conditions, forming a variety of carbonyl-type sites such as ammonium carbamates, bicarbonates, carbonates, carbamic acids, ureas, and amides. 13 C solid-state nuclear magnetic resonance (NMR) is often used to help elucidate the identity of these reacted species, but it is challenging due to the narrow chemical shift range of carbonyl sites. Herein, we demonstrate the application of a two-dimensional (2D) chemical shift anisotropy (CSA) recoupling pulse sequence (ROCSA) to obtain CSA tensor values at each isotropic chemical shift, overcoming limitations of isotropic peak resolution. CSA tensor values describe the local chemical environment and can readily differentiate between the chemisorbed and degradation products. To aid identification, we also developed a k-nearest neighbor (kNN) classification model to distinguish the functional groups via their CSA tensor parameters. This methodology was demonstrated on poly(ethylenimine) in γ-Al 2 O 3 exposed to CO 2 and showed that the chemisorbed products are ammonium carbamate and a mixed carbamate–carbamic acid species. The sample was analyzed again after desorption at 100 °C inducing mild degradation, and the remaining products were strongly bound carbamate and urea species. In conclusion, the combination of 2D CSA measurements coupled with a kNN classification model enhances the ability to accurately identify chemisorbed or degradation products in complex carbon capture materials.

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

Reactivity of [Tism Pr i Benz ]MgH and [Tism Pr i Benz ]MgMe towards Carbonyl Compounds: Access to Terminal Alkoxide and Enolate Complexes

The hydride and methyl compounds [Tism Pr i Benz ]MgH and [Tism Pr i Benz ]MgMe undergo insertion of the carbonyl moieties of non-enolizable aldehydes and ketones such as PhCHO and Ph 2 CO into the Mg–H and Mg–Me bonds to form alkoxide compounds, namely [Tism Pr i Benz ]MgOCH 2 Ph, [Tism Pr i Benz ]MgOCHPh 2 , [Tism Pr i Benz ]MgMOCH(Me)Ph and [Tism Pr i Benz ]MgOCMePh 2 . In contrast to the insertion of the carbonyl moiety, the reactions of the enolizable ketones Me 2 CO and PhC(O)Me with [Tism Pr i Benz ]MgMe afford the enolate complexes, [Tism Pr i Benz ]MgOC(Me)=CH 2 and [Tism Pr i Benz ]MgOC(Ph)=CH 2 . The formation of [Tism Pr i Benz ]MgOC(Me)=CH 2 is of note because methyl Grignard reagents preferentially react with acetone to form t-butoxide derivatives. The hydride compound, [Tism Pr i Benz ]MgH, also reacts with acetone to yield the enolate compound, [Tism Pr i Benz ]MgOC(Me)=CH 2 , but while the overall transformation is similar to that of the methyl derivative, [Tism Pr i Benz ]MgMe, the enolate compound is not the initially formed product. Specifically, acetone undergoes preferential insertion into the Mg–H bond to generate the corresponding alkoxide, [Tism Pr i Benz ]MgOPr i , which subsequently converts to the respective enolate in the presence of excess acetone. Furthermore, the relative ability of the hydride and methyl compounds to undergo insertion of carbonyl compounds into the Mg–H and Mg–Me bonds has been addressed computationally, which indicates that the barrier for insertion of the carbonyl group into the Mg–H bond is lower than that for insertion into the Mg–Me bond. The molecular structures of [Tism Pr i Benz ]MgOCH 2 Ph, [Tism Pr i Benz ]MgOCHPh 2 , [Tism Pr i Benz ]MgOCMePh 2 , [Tism Pr i Benz ]MgOC(Me)=CH 2 and [Tism Pr i Benz ]MgOC(Ph)=CH 2 have been determined by X-ray diffraction.

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

Chloride, Alkoxide, or Silicon: The Bridging Ligand Dictates the Spin State in Dicobalt Expanded Pincer Complexes

We report the synthesis and characterization of a series of high- and low-spin dicobalt complexes of the tBu PNNP expanded pincer ligand. Reacting this dinucleating ligand in its neutral form with two equiv of CoCl 2 (tetrahydrofuran) 1.5 yields a high-spin dicobalt complex featuring one Co inside and one Co outside of the dinucleating pocket. Performing the same reaction in the presence of two equivalents of KOtBu provides access to a high-spin dicobalt complex wherein both Co centers are bound within the PNNP pocket, and this complex also features a bridging OtBu ligand. Reacting either of the high-spin complexes with excess diethyl silane affords a low-spin dicobalt complex containing two unusual bridging Si-based ligands. These complexes were investigated using NMR spectroscopy, XAS, single crystal X-ray structure determination, and computational methods, showing that the Si-based ligands are best described as base-stabilized silylenes.

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

Theory of Cation Solvation in the Helmholtz Layer of Li-Ion Battery Electrolytes

The solvation environments of Li + in conventional nonaqueous battery electrolytes, such as LiPF 6 in mixtures of ethylene carbaronate (EC) and ethyl methyl carbonate (EMC), are often used to rationalize transport properties and solid electrolyte interphase (SEI) formation. Solvation environments in the compact electrical double layer (EDL) next to the electrode, also known as the Helmholtz layer, determine (partially) what species can react to form the SEI, with bulk solvation environments often being used as a proxy. Here, we develop and test a theory of cation solvation in the Helmholtz layer of nonaqueous Li-ion battery electrolytes. First, we validate the theory against bulk and diffuse EDL atomistic molecular dynamics (MD) simulations of LiPF 6 EC/EMC mixtures as a function of surface charge, where we find the theory can qualitatively capture the solvation environments. Next, we turn to the Helmholtz layer, where we find the main effect of the solvation structures next to the electrode is an apparent reduction in the number of binding sites between Li + and the solvents, again where we find reasonable agreement with our developed theory. Finally, by solving a simplified version of the theory, we find that the probability of Li + binding to each solvent remains equal to the bulk probability, suggesting that the bulk solvation environments are a reasonable place to start when understanding battery electrolytes. Our developed formalism can be parametrized from bulk MD simulations and used to predict the solvation environments in the Helmholtz layer through reducing the number of available coordination sites, which can be used to determine what could react and form the SEI.

Helmholtz↗
DOE OSTI2023

Experimental and Computational Mechanisms that Govern Long-Term Stability of CO 2 -Adsorbed ZIF-8-Based Porous Liquids

Porous liquids (PLs) based on the zeolitic imidazole framework ZIF-8 are attractive systems for carbon capture since the hydrophobic ZIF framework can be solvated in aqueous solvent systems without porous host degradation. However, solid ZIF-8 is known to degrade when exposed to CO 2 in wet environments, and therefore the long-term stability of ZIF-8-based PLs is unknown. Here, through aging experiments, the long-term stability of a ZIF-8 PL formed using the water, ethylene glycol, and 2-methylimidazole solvent system was systematically examined, and the mechanisms of degradation were elucidated. The PL was found to be stable for several weeks, with no ZIF framework degradation observed after aging in N 2 or air. However, for PLs aged in a CO 2 atmosphere, formation of a secondary phase occurred within 1 day from the degradation of the ZIF-8 framework. From the computational and structural evaluation of the effects of CO 2 on the PL solvent mixture, it was identified that the basic environment of the PL caused ethylene glycol to react with CO 2 forming carbonate species. These carbonate species further react within the PL to degrade ZIF-8. The mechanisms governing this process involves a multistep pathway for PL degradation and lays out a long-term evaluation strategy of PLs for carbon capture. Additionally, it clearly demonstrates the need to examine the reactivity and aging properties of all components in these complex PL systems in order to fully assess their stabilities and lifetimes.

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

Reaction Layer Formation on MgO in the Presence of Humidity

Mineralization by MgO is an attractive potential strategy for direct air capture (DAC) of CO 2 due to its tendency to form carbonate phases upon exposure to water and CO 2 . Hydration of MgO during this process is typically assumed to not be rate limiting, even at ambient temperatures. However, surface passivation by hydrated phases likely reduces the CO 2 capture capacity. Here, we examine the initial hydration reactions that occur on MgO(100) surfaces to determine whether they could potentially impact CO 2 uptake. We first used atomic force microscopy (AFM) to explore changes in reaction layers in water (pH = 6 and 12) and MgO-saturated solution (pH = 11) and found the reaction layers on MgO are heterogeneous and nonuniform. To determine how relative humidity (R.H.) affects reactivity, we reacted samples at room temperature in nominally dry N 2 (~11–12% R.H.) for up to 12 h, in humid (>95% R.H.) N 2 for 5, 10, and 15 min, and in air at 33 and 75% R.H. for 8 days. X-ray reflectivity and electron microscopy analysis of the samples reveal that hydrated phases form rapidly upon exposure to humid air, but the growth of the hydrated reaction layer slows after its initial formation. Reaction layer thickness is strongly correlated with R.H., with denser reaction layers forming in 75% R.H. compared with 33% R.H. or nominally dry N 2 . The reaction layers are likely amorphous or poorly crystalline based on grazing incidence X-ray diffraction measurements. After exposure to 75% R.H. in air for 8 days, the reaction layer increases in density as compared to the sample reacted in humid N 2 for 5–15 min. Further, this may represent an initial step toward the crystallization of the reaction layer. Overall, high R.H. favors the formation of a hydrated, disordered layer on MgO. Based on our results, DAC in a location with a higher R.H. will be favorable, but growth may slow significantly from initial rates even on short timescales, presumably due to surface passivation.

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

Reactive Species and Reaction Pathways for the Oxidative Cleavage of 4-Octene and Oleic Acid with H 2 O 2 over Tungsten Oxide Catalysts

The oxidative cleavage of carbon-carbon double bonds (C=C) in alkenes and fatty acids produces aldehydes and acids valued as chemical intermediates. Solid tungsten oxide catalysts are low cost, non-toxic, and selective for the oxidative cleavage of C=C bonds with hydrogen peroxide (H 2 O 2 ) and are, therefore, a promising option for continuous processes. Despite the relevance of these materials, the elementary steps involved and their sensitivity to the form of W sites present on surfaces have not been described. In this work, we combine in situ spectroscopy and rate measurements to identify significant steps in the reaction, the reactive species present on catalysts, and examine differences between the kinetics of this reaction on isolated W atoms grafted to alumina and on those exposed on crystalline WO 3 nanoparticles. Raman spectroscopy shows that W-peroxo complexes (W-(η 2 -O 2 ) formed from H 2 O 2 ) react in a kinetically relevant step to produce epoxides, which undergo hydrolysis at protic surface sites. Subsequently, the CH 3 CN solvent deprotonates diols to form alpha hydroxy ketones that react to form aldehydes and water following nucleophilic attack of H 2 O 2 . Turnover rates for oxidative cleavage, determined by in situ site titrations, on WO x -Al 2 O 3 are 75% greater than those on WO 3 at standard conditions. These differences reflect activation enthalpies (ΔH ‡ ) for oxidative cleavage of 4-octene that are much lower for the isolated WO x sites (36±3 and 60±6 kJ∙mol -1 for WO x -Al 2 O 3 and WO 3 , respectively) and correlate strongly with the difference between enthalpies of adsorption for epoxyoctane (ΔH ads,epox ), which resembles the transition state for epoxidation. The WO x -Al 2 O 3 catalysts mediate oxidative cleavage of oleic acid with H 2 O 2 following a mechanism comparable to that for 4-octene. The WO 3 materials, however, form only the epoxide and do not cleave the C-C bond or produce aldehydes and acids. These differences reflect the distinct site requirements for these reaction pathways and indicate that acid sites required for diol formation are strongly inhibited by oleic acids and epoxides on WO 3 whereas the Al 2 O 3 support provides sites competent for this reaction and increase the yield of the oxidative cleavage products.

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

Influence of the Molecular Structure on the Electrocatalytic Hydrogenation of Carbonyl Groups and H 2 Evolution on Pd

We investigated the electrocatalytic hydrogenation (ECH) of model aldehydes and ketones over carbon supported Pd in aqueous phase. We propose reaction mechanisms based on kinetic measurements and on spectroscopic and electrochemical characterization of the working catalyst. The reaction rates of ECH and of the H 2 evolution reaction (HER) vary with the applied electric potential following trends that strongly depend on the organic substrate. The intrinsic rates of hydrogenation and H 2 evolution are influenced, in opposing ways, by the sorption of the reacting organic substrate. Strong interactions, i.e., higher standard free energies of adsorption of the organic compound, induce high hydrogenation rates but decrease the overall currents. The fast hydrogenation kinetics produces a hydrogen-depleted environment that kinetically hinders the HER and the bulk phase transition of Pd to a H-rich bulk Pd hydride, which is triggered by the applied potential in the absence of reacting organic compounds. As consequence of strong organic-metal interactions, hydrogenation dominates at low overpotential. However, the coverages of organic substrates on the metal surface decrease and the rates of H 2 evolution surpass those of hydrogenation with increasingly negative electric potential. We determined the range of electric potential favoring hydrogenation on Pd and quantitatively deconvoluted the effects of the sorption of the organic compound, and of the rates of proton coupled electron transfers, on the kinetics of both ECH and HER. The results indicate that electrocatalysis offers hydrogenation pathways for polar molecules that are different and, in some cases, faster than those dominating in the absence of an external electric potential.

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

Polypropylene Degradation Catalyzed by Tantalum Hydrides Supported on Sulfated Alumina

Ta–H + sites supported on sulfated aluminum oxide (SAO), [Ta(H) 2 (O–) 2 ][SAO] (1), catalyze the hydrogenolysis of isotactic polypropylene (iPP, M n = 13.3 kDa; D = 2.4; mmmm = 94%) to form low-molecular-weight branched alkanes (C 11 –C 30 ) in good yields (70%). Here, the alkanes formed lose stereochemical information originating from iPP, but residual iPP remains highly isotactic. In the presence of D 2 , similar mixtures of alkanes are formed containing –CH 3–x D x , –CHD, and –CD–. Residual iPP maintains high tacticity and incorporates deuterium primarily into –CH 3 groups of the polymer (–CH 3–x D x /–CHD–/–CD– ~10:1:0). 1 reacts with pinacolborane to form [TaH(κ 2 -H 2 BPin)(O–) 2 ][SAO] (2). At 200 °C in an iPP melt, 2 reacts to form products indicative of C–H bond activation at –CH 3 groups and internal –CH 2 – groups in iPP in a ~1:3 ratio, indicating a slight kinetic preference for C–H bond activation at internal –CH 2 – groups. 1 is more reactive toward iPP than high-density polyethylene in hydrogenolysis reactions

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗