Search NASA⌕ Search

SEARCH · Search NASA

Results for “kMC”

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

Microstructure-process relationships in monolithic U-10Mo fuel foil single-pass rolling: A parametric simulation study

In this work, a previously validated coupling of Kinetic Monte Carlo (KMC) Potts Model and finite element method (FEM) simulations was implemented to investigate the effects of microstructural features in as-cast and homogenized monolithic U-10Mo foils on the emergent microstructure after rolling and reheating. Parameters that could potentially affect recrystallization behavior of the rolled U-10Mo foils were considered: grain size distribution, uranium carbide (UC) size distribution, UC volume fraction, spatial distribution of UC, and rolling reduction magnitude. Grain structure and the magnitude of rolling reduction have the strongest influence on recrystallization kinetics and the fabricated grain size distribution. The UC distribution had only a weak effect on the recrystallization kinetics and final microstructures. While particle-stimulated nucleation (PSN) occurred in simulation more frequently as grain size increased, its incidence did not appear to considerably affect the recrystallization kinetics or grain size distribution.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Gas generation and intramolecular isotope study in laboratory pyrolysis of the Springfield coal from the Illinois Basin

Position-specific (PS) isotopes of propane have been proposed as a potential geochemical tool to decipher various geological processes (e.g., thermal cracking, biodegradation, H exchange) in natural reservoirs. The limited studies have been conducted on the PS isotopes of propane from the pyrolysis gases from marine shales, and natural gases sourced from lacustrine and marine kerogens, but little is known on gases produced from the humic kerogen. Here, this study investigated the PS δ 13 C of propane in the closed-system pyrolysis of the Springfield coal, Illinois Basin, Indiana, at 310 to 470 °C (Easy %Ro: 0.76 to 3.07). The C kinetic isotope effect (KIE) of CH 4 produced in both this study and previous low-temperature pyrolysis of the same coal indicates the cleavage of C—O bonds is the main generation pathway at the early kerogen cracking stage, followed by the breakdown of alkyl groups. At the wet-gas cracking stage, C 3 H 8 production from thermally stable compounds has a significant influence on the bulk and position-specific C KIE in the pyrolysis of marine Woodford kerogen and Springfield coal. According to the PS δ 13 C of propane, the central site is likely more enriched in 13 C and the δ 13 C of the terminal site is relatively heterogeneous within the propyl group attached to different functional groups of the gas-prone kerogen. Our findings based on the pyrolysis experiments and natural gas samples indicate thermal cracking and biodegradation appear to alter the δ 13 C cen values more significantly than the δ 13 C ter values of propane. The larger magnitude of ΔC c-t in the kMC simulations (Peterson et al., 2018) compared with those from the marine shale and coal possibly implies the non-random distributions of 13 C of propane precursors in the kerogens. As a new dimension of intramolecular isotopic information of propane, the PS δ 13 C values can contribute to fingerprinting the gas origins and identifying the various geological processes (e.g., kerogen cracking, wet-gas cracking, microbial activities) in sedimentary basins.

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

Bay capping via acetylene addition to polycyclic aromatic hydrocarbons: Mechanism and kinetics

The bay-capping mechanism on PAH armchair edges and the kinetics of acetylene addition to 6–6–5 and 5–6–5 bays have been explored by ab initio/RRKM-ME calculations. The bays on the edges were modeled by C 21 H 11 and C 20 H 9 radicals produced by H abstractions from 7H-benzo[c]cyclopenta[e]pyrene and dicyclopenta[cf]pyrene. The C 20 H 9 + C 2 H 2 reaction is shown to have a low entrance barrier and to rapidly form the capped product, indaceno[2,1,8,7-cdefg]pyrene, along with ethynyl substituted dicyclopenta[cf]pyrene at temperatures above 1400 K. The reactivity of C 21 H 11 is shown to be governed by the location of the unpaired electron; the π radical R1 formed by H abstraction from the CH 2 group in 7H-benzo[c]cyclopenta[e]pyrene reacts with C 2 H 2 very slowly owing to a high entrance barrier, with the bay-capping rate constant approaching 10 -16 cm 3 molecule -1 s -1 only at temperatures above 2000K. Further, this result reaffirms that the growth of π aryl radicals via acetylene addition is inefficient and reflects the generally low reactivity of such radicals where the spin density is highly delocalized over the entire polyaromatic system. Alternatively, the σ C 21 H 11 radical R2 produced by H abstraction from the five-membered ring at the bay rapidly reacts with C 2 H 2 forming the bay-capped product, with the rate constant on the order of 10 -12 cm 3 molecule -1 s -1 at T≥1500K. Rate constants for the capping reactions at the 6–6–5 and 5–6–5 bays are compared with those at the 6–0–6, 6–6–6, and 6–5–6 bays. The site-specific bay-capping rate constants have been utilized in kMC simulations of the PAH growth and the results showed measurable differences when the 6–6–5 and 5–6–5 bay-capping reactions are taken into account, including an increase of the growth rate and the formation of closed-shell PAH and a rise of the number of embedded five-membered rings accompanied with a slight decrease of their overall amount.

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

A simple model for short-range ordering kinetics in multi-principal element alloys

Short-range ordering (SRO) in multi-principal element alloys influences material properties such as strength and corrosion. While some degree of SRO is expected at equilibrium, predicting the kinetics of its formation is challenging. We present a simplified isothermal concentration-wave (CW) model to estimate an effective relaxation time of SRO formation. Estimates from the CW model agree to within a factor of five with relaxation times obtained from kinetic Monte Carlo (kMC) simulations when above the highest ordering instability temperature. Further, the advantage of the CW model is that it only requires mobility and thermodynamic parameters, which are readily obtained from alloy mobility databases and Metropolis Monte Carlo simulations, respectively. The simple parameterization of the CW model and its analytical nature makes it an attractive tool for the design of processing conditions to promote or suppress SRO in multicomponent alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Maintaining a Flat Li Surface during the Li Stripping Process via Interface Design

Electroplating has been the main focus in mitigating the dendrite growth on the Li-metal electrode; however, the stripping process is equally critical, since the nonsmooth Li surface during stripping will lead to nonuniform local current density, planting the seeds for dendrite growth. In this paper, density functional theory (DFT) and kinetic Monte Carlo (KMC) techniques were combined to investigate the vacancy evolution in Li interfaced with different solid–electrolyte interphase (SEI) materials. It was found that the lithiophilic interface, such as Li/Li 2 O, repels vacancies into the bulk Li, so Li atoms can quickly fill the Li vacancies near the Li/Li 2 O interface and maintain a smooth Li surface. In contrast, the lithiophobic interface, such as Li/LiF, traps Li vacancies toward the interface, and the accumulated Li vacancies form voids and roughen the surface. The predicted critical stripping current density, below which a smooth Li surface will be maintained, is therefore much faster at the lithiophilic interface than that at the lithiophobic interface. It was further revealed that the lithiophilicity at different SEI or coating materials can be ranked as Li/Li 2 O > Li/LiPON > Li/Li 2 CO 3 > Li/LiF based on the calculated interfacial adhesion and accumulation of electron density at the interface. Furthermore, this suggests that interface and coating design at nanoscale can be effective for maintaining a smooth Li surface during the stripping process, solving another challenge to achieving a dendrite-free Li-metal electrode in both liquid and solid electrolytes.

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

Thermal Shape Stability of fcc Metal Nanocrystals Synthesized with Faceted Nonequilibrium Shapes

Highly refined capabilities of the shape-controlled solution-phase synthesis of metal nanocrystals (NCs) allow the generation of NCs with faceted nonequilibrium shapes, which optimize properties for target applications such as catalysis and plasmonics. Often, for such applications and also for TEM analysis, the NCs are removed from the solution-phase environment. We explore the postsynthesis evolution of these metastable NCs in a high-vacuum TEM environment. Specifically, here we analyze their reshaping toward the equilibrium Wulff shapes mediated by surface diffusion, where such reshaping degrades the above-mentioned optimized properties. Typical sizes for these NCs range from 5 to 30 nm or 10 3 –10 6 atoms, and reshaping often occurs on the time scale of minutes for temperatures around, say, 400 °C. We discuss the development of predictive stochastic atomistic-level models for NC evolution with a realistic description of surface diffusion. These models, in contrast to Molecular Dynamics, can naturally address the relevant time and length scales for these systems. KMC simulation results for the stochastic models are described, focusing on the reshaping of slightly elongated nanorods and of mildly truncated octahedra and nanocubes. In addition, we review appropriate theoretical formulations for reshaping, which involves the nucleation and growth on 2D islands or layers on outer facets of the NC. We note the limitations of classical nucleation theory in some scenarios and demonstrate the successes of a more fundamental and general master equation-based analysis.

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

Theoretical Perspective of Promoting Direct Methane-to-Methanol Conversion at Complex Metal Oxide–Metal Interfaces

Direct methane conversion to methanol has been considered as an effective and economic way to release greenhouse effects and current high demand on methanol in industry. However, the process has long been challenging due to lack of viable catalysts to compromise the activation of methane that typically occurs high temperatures and retaining of produced methanol that requires mild conditions. Here, this perspective demonstrates an effective strategy to promote direct methane to methanol conversion by engineering the active sites and chemical environments at complex metal oxide - copper oxide - copper interfaces. Such effort strongly depends on the extensive theoretical studies by combining density functional theory (DFT) calculations and kinetic Monte Carlo (KMC) simulations to provide in-depth understanding of reaction mechanism and active sites, which built the strong basis to enable the identification of design principles and advance the catalyst optimization for selective CH 4 -to-CH 3 OH conversion.

03 NATURAL GAS↗
DOE OSTI2023

Correlating Experimentally Determined Hydrogen Binding Energy with Hydrogen Evolution Activity over Metal Monolayers on Molybdenum Nitride

It is well established that hydrogen binding energy (HBE) is a key descriptor for hydrogen evolution reaction (HER) activity, and such a relationship is a useful tool for searching efficient and cost-effective HER catalysts. However, in almost all cases, the HBE values are obtained from density functional theory (DFT) calculations. Here in this study, temperature programmed desorption (TPD) was used to experimentally determine the HBE values of metal monolayers supported on molybdenum nitride (Mo 2 N), and electrochemical measurements were performed on the same surfaces. Combined DFT and kinetic Monte Carlo (kMC) simulations were used to validate the trend observed with TPD and the electrochemical HER activity. Depositing one monolayer of Pt on Mo 2 N led to similar HBE values seen for bulk Pt, and electrochemical measurements showed that monolayer Pt on Mo 2 N had similar HER activity to that of bulk Pt. Similar studies were also performed for monolayer Pd on Mo 2 N to confirm the correlation. This work also demonstrates that Mo 2 N is a promising support to reduce precious metal loading in HER catalysts.

08 HYDROGEN↗
DOE OSTI2020

Methanol Synthesis from CO 2 Hydrogenation over a Potassium-Promoted Cu x O/Cu(111) ( x ≤ 2) Model Surface: Rationalizing the Potential of Potassium in Catalysis

Alkalis have been reported as a promotor in the heterogeneous catalysis, being able to enhance the activity and selectivity of catalysts. The effective utilization of alkalis in catalyst optimization requires the fundamental understanding of the underlying mechanism. In this work, we take a potassium (K)-modified Cu x O/Cu(111) (x ≤ 2) model surface as a case study to rationalize the nature of K during the carbon dioxide hydrogenation using combined density functional theory (DFT) calculation and the kinetic Monte Carlo (KMC) simulation. Our result demonstrates the significant tuning of selectivity from carbon monoxide to methanol on going from Cu(111) to K-modified Cu x O/Cu(111). The deposited K + stabilizes the Cu x O thin film under the reducing condition of carbon dioxide hydrogenation. More importantly, our study reveals that K + acts as an active center for selective tuning in the binding, an accelerator for charge transfer, and a mediator for the electron tunneling. As a result, the K-modified Cu x O/Cu(111) opens a methanediol [H 2 C(OH) 2 ]-mediated formate pathway to facilitate the selective conversion of carbon dioxide to methanol. Our study develops the intrinsic rules of design to tune the catalytic performance using alkali metals.

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

SEI Formation and Lithium-Ion Electrodeposition Dynamics in Lithium Metal Batteries via First-Principles Kinetic Monte Carlo Modeling

The stabilization and enhanced performance of lithium metal batteries (LMBs) depend on the formation and evolution of the Solid Electrolyte Interphase (SEI) layer as a critical component for regulating the Li metal electrodeposition processes. This study employs a first-principles kinetic Monte Carlo (kMC) model to simulate the SEI formation and Li + electrodeposition processes on a lithium metal anode, integrating both the electrochemical electrolyte reduction reactions and the diffusion events giving place to the SEI aggregation processes during battery charge and discharge processes. The model replicates the competitive interactions between organic and inorganic SEI components, emphasizing the influence of the cycling regime. Results indicate that grain boundaries within the SEI facilitate faster lithium-ion transport compared to crystalline regions, crucial for improving the performance and stability of LMBs. The findings underscore the importance of dynamic SEI modeling for further development of next-generation high-energy-density batteries.

25 ENERGY STORAGE↗
DOE OSTI2025

Unveiling the Origin of Morphological Instability in Topologically Complex Electrocatalytic Nanostructures

Coarsening and degradation phenomena in metals have largely focused on thermally driven processes, such as bulk and surface diffusion. However, dramatic coarsening has been reported in high-surface-area, nanometer-sized Pt-based catalysts during potential cycling in an electrolyte at room temperature─a temperature too low for the process to be explained purely by surface mobility values measured in both vacuum and electrolytes (∼10 –22 and ∼10 –18 cm 2 /s, respectively). This morphological evolution must be due to a different mechanism for mass transport that is sensitive to electrochemical conditions (e.g., electrolyte composition, potential limits, and scan rate). However, there have been no notable studies of electrochemically induced coarsening in nanometer-sized electrocatalysts. Here, we unveil the origins of coarsening in an electrolyte through coupled in situ experiments and atomistic kinetic Monte Carlo (kMC) simulations. Our work demonstrates electrochemical coarsening is driven by two concurrent mechanisms that can be explained at the atomistic level: (i) dissolution/redeposition during the reduction of an oxidized species and (ii) rapid surface diffusion of undercoordinated atoms.

gold↗
DOE OSTI2020

Final Technical Report for project entitled Highly Active, Durable, and Ultra-Low PGM NSTF Thin Film ORR Catalysts and Support

In this project, the objective was to develop new oxygen reduction reaction (ORR) electrocatalysts for proton exchange membrane fuel cells (PEMFCs) which could exceed all of the Department of Energy (DOE) 2020 targets listed in DE-FOA-0001224, Subtopic 1b, Table I. The expected outcome was development of one or more electrocatalysts which are substantially improved over the current state-of-the-art in terms of overall activity, durability, and cost, and are suitable for automotive traction and stationary fuel cell applications. This report summarizes this project’s progress towards meeting the stated objectives and includes a summary of key findings and conclusions. The development activities towards new highly active and durable thin film electrocatalysts has led to the discovery of several electrocatalysts which approached or exceeded several DOE 2020 targets for activity, durability, and performance in proton exchange membrane fuel cell (membrane electrode assemblies (MEAs). In this project, electrocatalyst development focused on systematic physical and electrochemical characterization of electrocatalyst activity, durability, and performance in MEAs as a function of electrocatalyst fabrication, compositional and structural variables. The project approach was to develop relationships between the catalyst functional responses (activity, durability and performance) with catalyst physical properties and catalyst fabrication methods. Additionally, the thin film catalysts were integrated onto the unique and durable 3M Nanostructured Thin Film (NSTF) support, consisting of arrays of self-assembled organic crystalline whiskers. The catalysts were evaluated against the project targets via extensive electrochemical characterization in MEA format at 3M, advanced structural and compositional microscopy at Oak Ridge National Laboratory (ORNL), and atomic structure analysis via XAFS at Argonne National Laboratory (ANL). The extensive electrochemical and physical characterization resulted in development of several trends which correlate the catalysts’ electrochemical properties to their physical properties, including composition, structure, and method of fabrication. The development was guided by density functional theory (DFT) modeling at Purdue University and kinetic Monte Carlo (kMC) modeling at Johns Hopkins University. The catalyst simulations provided key insights into the observed experimental activity and durability trends, and additionally were utilized to assess new electrocatalyst concepts prior to or coincident with physical catalyst development. This development has led to several electrocatalyst candidates with activity and durability which approach or exceed DOE targets. One class of catalysts, based on nm-scale thin layers of Pt on Ir, met or exceeded 6 of the 6 DOE targets the project addressed.

08 HYDROGEN↗
DOE OSTI2023

Phase-Field Modeling of Materials Interfaces and Nanostructures

Nanostructured materials offer unique properties for a wide range of energy applications. Among various known processing methods, liquid metal dealloying (LMD)—the selective dissolution of a base alloy element into a metallic melt—has emerged as a powerful technique to produce a new class of nano-/meso-scale open porous and bicontinuous composite structures with ultra-high interfacial area. LMD has recently been complemented by the advent of vapor phase dealloying (VPD), a novel technique that exploits the selective evaporation of one element from a parent alloy containing elements with very different vapor pressures, thereby enabling to fabricate open nanoporous structures of various elements from less-noble metals to inorganic elements regardless of their chemical activity without requiring high LMD temperatures and chemical etching. Together LMD and VPD have greatly expanded the scope of dealloying techniques, limited by traditional electrochemical means to noble metals, and boosted the design of new functional and structural materials that combine a wide variety of elements. Topologically-connected open porous structures with ultra-high surface area allow mass transport within the structure while preserving structural integrity, enabling them to serve as catalysts, fuel cells, supercapacitors, or high-capacity battery materials. Bicontinuous composite structures in turn can display high strength and high ductility or superior radiation-damage resistance due to the ultra-high interface area between interpenetrating solid phases. This research program makes use of state-of-the-art computational methods to understand at a basic level the self-organizing dealloying process with main focus on dealloying kinetics and interfacial pattern formation at the dealloying front controlling initial structure size, topology, and phase compositions. Phase-field simulation studies of LMD focus on solid solutions, line compounds, and intermetallic systems that can form ternary composites by nucleation and growth of a new phase. Studies of VPD employ phase-field modeling and a hybrid method combining a kinetic Monte Carlo (KMC) model of evaporation and surface diffusion with molecular dynamics for vapor-phase transport inside nanopores. Simulations explore mechanisms of interface- and diffusion-controlled dealloying kinetics, both observed in VPD but not fundamentally understood. In addition, our newly developed multi-physics phase-field approach of large-volume-change phase transformations is being used to model novel 3D anode geometries based on dealloyed nanoporous structures including novel sandwiched graphene/Si/silica for highrate Li ion battery. Those studies are aimed at elucidating geometric design principles that improve mechanical stability. We expect this research to enhance the capability to tailor nano-/mesoscale structures for a wide range of energy-related materials applications and to yield further advances in computational methodologies that benefit a broad materials research community.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2024

Dynamic Binary Complexes (DBC) as Super-Adjustable Viscosity Modifiers for Hydraulic Fracturing Fluids

In the preceding project year two, we refined three DBC formulations from a selection of over 50 different chemistries. The optimization study primarily encompassed testing for (i) reversibility extent, (ii) performance in the presence of chemical additives, (iii) adhesion and friction behavior during displacement in wellbores and pipelines, (iv) corrosion protection performance, and (v) injection performance with model fracture systems at the laboratory scale. Highly promising results obtained from all these tests signify the significant potential of DBCs in enhancing hydrocarbon recovery from unconventional reservoirs. The primary activities in the third project year included publishing experimental findings across multiple articles and conducting outreach initiatives. Throughout the year, we undertook tasks such as replicating experimental results, further optimizing various formulations and their associated experimental sets, and conducting additional tests to address missing components based on reviewer feedback and suggestions. We also explored the surfactant and friction-reduction aspects of selected formulations through drag reduction tests. In addition, we constructed an improved fracturing performance setup and performed flow injection tests. The specific DBC formulations focused on during this project period were A8/B1, A12/B5, and A10/B12. We also obtained results for additional DBC formulations and a select few commercial fracturing fluids for the purpose of comparison. Within the project's scope, we aim to enhance the experimental findings with the development of various models. The first two years focused on two key aspects: (i) the creation of a high-fidelity hydraulic fracturing model for non-Newtonian fluids to gain insights into the implementation of DBC fluids in fracking environments, and (ii) the development of a multiphase flow simulator for estimating total production, fluid saturation in the reservoir, and the creation of a fracture propagation model and kinetic Monte Carlo (kMC) models for diverse applications. In the third year, we delved into the fundamental nanostructural properties of DBCs, exploring aspects such as material chemistry, pH tunability, and control of DBC formation and stability. Subsequently, in the extension year, we conducted a systematic investigation of various building blocks containing primary, secondary, and tertiary amine functional groups to understand their impact on rheological and viscoelastic properties. Furthermore, we explored a Dissipative Particle Dynamics (DPD) model to simulate self-assembly processes with precision, creating a high-fidelity representation of relevant nanostructures. The tasks performed this year with the significant results obtained have been discussed in Section 2.

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

Self-Learning Kinetic Monte Carlo Simulations of Radiation Damage in Nuclear Fuels

Understanding how irradiation affects the thermo-physical and mechanical properties of nuclear materials, such as thermal conductivity degradation in fuels and embrittlement of structural components, is critical to the safety and efficiency of nuclear reactors. These effects are largely governed by the formation and evolution of atomic-scale point defects and defect clusters. Due to their small sizes, however, these defects are invisible under high-resolution scanning transmission electron microscopy. This project aims to fill this experimental knowledge gap by integrating density functional theory (DFT), machine learning interatomic potential (MLIP), and kinetic Monte Carlo (KMC) techniques to predict longtime evolution of irradiation-induced defects in nuclear fuels.

36 - MATERIALS SCIENCE↗
DOE OSTI2024

New insights on the effects of chemistry and temperature on α’ precipitation during aging of FeCrAl alloys

Here, the effect of Al on α' precipitation in FeCrAl alloys was studied through thermal aging of several model binary FeCr and FeCrAl alloys with Cr content of 13 wt.%, 17 wt.%, and 25 wt.% with and without 5.5 wt.% Al. Aging was performed up to 1,000h at varying temperatures from 400-500°C. At 400°C no age hardening was observed due to slow kinetics at this temperature. For both the 17Cr and 25Cr alloys, the addition of Al shows a lowering of the miscibility gap, consistent with other reports in literature. Interestingly, however, for the 25Cr alloys the addition of Al in the FeCrAl ternary alloy accelerated α' precipitation below the miscibility gap. Such enhanced precipitation has also been predicted by our atomistic kinetic Monte Carlo (AKMC) simulations. While previous literature has often focused on Al suppressing α', here we show that while Al can lower the miscibility gap in the Fe-Cr ferritic system thermodynamically, it may also enhance the kinetics of precipitation.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Ethanol Conversion over La 0.7 Sr 0.3 MnO 3–x (100): Autocatalysis, Adjacent O-Vacancies, Disproportionation, and Dehydrogenation

The mechanism for catalytic conversion of ethanol over La 0.7 Sr 0.3 MnO 3– x (100) surface to acetaldehyde and ethene was investigated. Pre-exposure temperature-programmed reaction (PE-TPR) experiments were performed in which ethanol was introduced to oxidized or reduced surfaces followed by heating. In particular, sequential PE-TPR experiments were conducted to incrementally and gradually reduce the surface. The products and their ratios were investigated as a function of surface reduction. The data show that acetaldehyde and ethene production is catalyzed with hydrogen abstraction and oxygen abstraction reactions occurring by intermediates in vacancies at various temperatures >400 K. Adsorption of acetaldehyde followed by a temperature-programmed reaction does not produce ethene, indicating that acetaldehyde is not an intermediate to ethene and that the hydrogen and oxygen abstraction from ethanol to ethene are decoupled steps. Further evidence for this mechanistic nuance was obtained using isotopically labeled ethanol (CD 3 CH 2 OH), which produces CD 3 CHO and CD 2 CH 2 . Additionally, the ratio of aldehyde production to alkene production increases with reduction, suggesting that aldehyde is produced from a disproportionation reaction between ethoxy species in adjacent O-vacancies, while ethene is produced from a dehydrogenation reaction with ethoxy species in vacancies without requiring adjacent O-vacancies. Counterintuitively, this finding indicates that the more oxygenated product (aldehyde vs ethene) is favored with more vacancies and that the net alcohol conversion is autocatalytic. Density functional theory calculations were able to find the previously unknown disproportionation pathway between ethoxies in adjacent O-vacancies, and kinetic Monte Carlo simulations support this interpretation by reproducing experimental selectivities. The activation energies for these pathways are estimated as 132 ± 10 kJ/mol for the disproportionation reaction (when occurring between ethoxies in adjacent vacancies) and as 148 ± 11 kJ/mol for the direct dehydrogenation reaction of an ethoxy in a vacancy. Based on these results, a mechanism with operative pathways based on elementary steps in O-vacancies is reported.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗