Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Transition state”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

HȮ 2 + HȮ 2 : High level theory and the role of singlet channels

At low temperatures and high pressures, the HȮ 2 + HȮ 2 reaction is an important reaction in combustion. There are significant unresolved discrepancies between prior theoretical and experimental studies of this reaction. It has generally been presumed to occur as an abstraction on the triplet surface to produce H 2 O 2 + 3 O 2 . Here, we employ a combination of high-level electronic structure theory (the ANL0 composite method or multi-reference methods as appropriate), sophisticated transition state theory (vibrationally adiabatic torsions, variational, and variable reaction coordinate), and master equation analyses to predict the thermal kinetics on the H 2 O 4 surface. Notably, this analysis suggests a significant branching to HȮ 3 + ȮH near 1000 K via reaction on the singlet surface. This channel does not appear to have been considered in prior combustion models. HȮ 3 itself is a metastable complex (bound by only 3 kcal mol -1 ) that rapidly dissociates to ȮH + 3 O 2 . Thus, the net reaction for this channel, HȮ 2 + HȮ 2 → ȮH + ȮH + O 2 , converts two low reactivity HȮ 2 radicals into two highly reactive ȮH radicals. The ramifications of these newly derived rate expressions are highlighted through kinetic modeling studies of H 2 , CH 3 OH, $\textit{n}$-heptane, and isooctane oxidation; all at the high pressures of relevance to combustion devices.

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

An experimental and kinetic modeling study of the pyrolysis of isoprene, a significant biogenic hydrocarbon in naturally occurring vegetation fires

Isoprene dominates the carbon flux emitted by vegetation and constitutes 40% of non-methane biogenic emissions worldwide. Despite pyrolysis experiments at temperatures above 1000 K showing a link between isoprene combustion and aromatic species formation, comprehensive mechanistic research on isoprene is scarce in the literature. Here, in this work, we carry out an experimental and theoretical study to build, for the first time, a chemical kinetic model describing isoprene pyrolysis. The formation of polycyclic aromatic hydrocarbon (PAH) precursor species, often observed in vegetation fire plumes, is partially explained by isoprene pyrolysis experiments and theoretical modeling. Molecular dynamics (MD) simulations unveil reaction pathways from allylic isoprenyl radicals to allene and cyclopentadiene (CPD) intermediates, two relevant species detected in the experiments. Rate constants for these identified pathways are calculated using variational transition state theory to update the kinetic model, which is validated against single-pulse shock tube (SPST), and jet-stirred reactor (JSR) experimental data in the temperature range of 850–1690 K. The kinetic model presents satisfactory agreement with the SPST experimental data, and a reaction pathway analysis shows that association of propargyl radicals results in benzene formation. The JSR pathway analysis also identifies the prominent reactions for CPD, benzene, styrene, and toluene formation. Our model does not reproduce the CPD experimental profiles, indicating that additional studies are necessary. Overall, our findings advance the understanding of isoprene pyrolysis and its related atmospheric pollutants in naturally occurring vegetation fires where smoldering and oxygen-deficient combustion processes are present.

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

From electronic structure to model application of key reactions for gasoline/alcohol combustion: Hydrogen-atom abstractions by $CH_3\dot{O}$ radical

H-atom abstraction by methoxy radical ($CH_3\dot{O}$) plays an important role in capturing the kinetics of reactions between gasoline components and alcohols. This study focuses on determining the reaction rates and thermodynamic properties of methoxy radical reactions with five gasoline fuel components: n-heptane, iso-octane, 1-hexene, cyclopentane and toluene. Electronic structure calculations were per-formed for all the stationary points with M06-2X/6 -311 ++ g(d,p) method. G3 composite method with atomization method is used for determining Δ f H 0 of all the closed shell and radical species, using which the necessary thermodynamic data of all the species was determined. Coupled cluster theory QCISD(T)/cc-pVXZ (where X = D and T) and Møller-Plesset perturbation theory MP2/cc-pVXZ (where X = D, T and Q) were used to calculate single point energies. Subsequently, rate constants for all hydrogen atom ab-straction channels have been performed using conventional transition state theory with unsymmetric tunneling corrections. A systematic comparison of rates for abstraction from different sites within the same species and same site from different species is done in order to get insights into this reaction class. Here, the computed thermodynamic properties and rate constants were incorporated into a recent gasoline mechanism to investigate the impact of the calculations performed in this work. A shift in predicted NTC (negative temperature coefficient) behavior and significant reduction in model reactivity is observed upon incorporating the rates calculated herein.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2025

Experimental and modeling study of the autoignition behavior of a saturated heterocycle: Pyrrolidine

Experiments are conducted in both rapid compression machine (RCM) and shock tube (ST) to better quantify autoignition behavior (e.g., ignition delay, heat release) and understand heteroatomic effects in heterocyclic compounds, which are important reference components for the combustion of biomass-derived liquid fuels. These tests focus on the nitrogen-containing, five-membered saturated ring, pyrrolidine, at diluted conditions covering pressures of 20 and 50 bar, temperatures of 720–1450 K and a range of stoichiometries (ϕ = 0.5–2). A chemical kinetic model is developed and coupled to an existing combustion kinetics framework describing key nitrogen containing intermediates (e.g. pyrrole, ammonia and NOx). H-abstraction reactions by OH, H, CH 3 and HO 2 , are determined using ab-initio transition state theory methods, while analogies to cyclopentane are adopted for many other reactions, such as ring-opening. Here, the autoignition measurements reveal the lack of negative temperature coefficient (NTC) behavior and low-temperature chemistry for pyrrolidine, as opposed to its saturated hydrocarbon analogue, cyclopentane. Interestingly, at the lowest temperatures (T < 750 K), the reactivity of cyclopentane is greater than pyrrolidine, while at higher temperatures, pyrrolidine becomes more reactive. Agreement between the experimental measurements and the model is good, and it is found that H-abstraction reactions by HO 2 and ensuing chemistry play key roles in controlling the reactivity of this cyclic amine. Most of the fluxes, i.e., >70 %, are predicted to move through 1- or 2-pyrroline (C 4 H 7 N) and then the cyclic C 4 H 6 N radical, at both lower and higher temperatures, to form either CH 2 CHCHCHNH via ring-opening or pyrrole via β-scission. It appears that the ring opens more easily at lower temperature whereas the C–H β-scission dominates at higher temperature and lower pressure, such that the reaction of the fuel radical intermediate carrying an unpaired electron on the nitrogen atom with HO 2 is the next most notable in promoting oxidation. When comparing pyrrolidine and cyclopentane, which exhibits distinct pathways in different temperature regimes, the pyrrolidine pathways and sensitivity analysis align more closely to the high temperature case of cyclopentane where the important role of HO 2 radicals is seen to provide chain branching through HO 2 reaction with the fuel, accompanied by H 2 O 2 formation and decomposition to OH. The formation of 5-membered diene rings and ring opening reactions are also found to be highly relevant. Of particular note, it is found that there is little influence of small molecule nitrogen-chemistry, e.g., NH 2 , HCN, NO/NO 2 on the reactivity of the pyrrolidine mixtures investigated here where no recirculated combustion gases are included.

Autoignition↗
DOE OSTI2020

An automated cluster surface scanning method for exploring reaction paths on metal-cluster surfaces

Metal-cluster surfaces present a wide variety of unique coordination environments. This complexity makes it difficult to manually probe the surface reactivity of such clusters. Here, we present a simple and automated method to systematically discover reaction pathways on cluster surfaces, based on the automated cluster surface scanning (ACSS) technique for mapping out potential energy surfaces. We showcase our method on 55-atom icosahedral Cu and Ag clusters, where we determine the activation energies of four elementary steps common in heterogeneous catalysis – hydrogen recombination (H* + H* → H 2 * + *), oxygen recombination (O* + O* → O 2 * + *), water formation (OH* + H* → H 2 O* + *), and CO oxidation (CO* + O* → CO 2 * + *) – with density functional theory calculations (DFT-PBE + D3). We show that the ACSS method requires significantly less human effort than the established manually performed climbing-image nudged elastic band (MP + CI-NEB) technique and locates transition states with comparable accuracy (root-mean-squared error of 0.10 eV) and similar computational cost. Rigorous sampling of the potential energy surface with the ACSS method allows one to locate all lowest-energy reaction pathways obtained via the MP+CI-NEB approach, as well as alternative pathways that one may have missed with the MP+CI-NEB approach due to the many possible pathways available on these clusters. The accuracy and efficiency afforded by the ACSS method could enable high-throughput exploration of the diverse reactivity of metal clusters.

36 MATERIALS SCIENCE↗
DOE OSTI2023

QMMM 2023: A program for combined quantum mechanical and molecular mechanical modeling and simulations

Combined quantum mechanical and molecular mechanical (QM/MM) methods play an important role in multiscale modeling and simulations. QMMM 2023 is a general-purpose program for single-point calculations, geometry optimizations, transition-state optimizations, and molecular dynamics (MD) at the QM/MM level. It calls a QM package and an MM package to perform the required single-level calculations and combines them into a QM/MM energy by a variety of schemes. QMMM 2023 supports GAMESS-US, Gaussian, and ORCA as QM packages and Tinker as the MM package. Four types of treatments are available for embedding the QM subsystem in the MM environment: mechanical embedding with gas-phase calculations of the QM region, electronic embedding that allows polarization of the QM region by the MM environment, polarizable embedding for mutual polarization of the QM and MM regions, and flexible embedding for both mutual polarization and partial charge transfer between the QM and MM regions. Boundaries between QM and MM regions that pass through covalent bonds can be treated by several methods, including the redistributed charge (RC) scheme, redistributed charge and dipole (RCD) scheme, balanced-RC scheme, balanced-RCD scheme, screened charge scheme that takes account of charge penetration effects, and smeared charge scheme that delocalizes the MM charges near the QM–MM boundary. Geometry optimization can be done using the optimizer implemented in QMMM 2023 or the Berny optimizer in Gaussian through external calls to Gaussian. Molecular dynamics simulations can be performed at the pure-MM level, pure-QM level, fixed-partitioning QM/MM level, and adaptive-partitioning QM/MM level. As a result, the adaptive-partitioning treatments permit on-the-fly relocation of the QM–MM boundary by dynamically reclassifying atoms or groups into the QM or MM subsystems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Preparation of a uranium monocarbide anode and electrochemical characterization in molten LiCl-KCl-UCl 3

Porous uranium carbide (UC) pellets possessing moderate electrical conductivity were synthesized by reaction of UO 2 with graphite at temperatures up to 1550°C under rough vacuum. Conversions as high as 98% were achieved at soak times of 2-4 hours. The electrochemistry of the UC pellets in molten LiCl-KCl-6.5 wt% UCl 3 was explored using a variety of techniques including DC polarization methods, cyclic voltammetry, chronopotentiometry and bulk electrolysis. Here, the electrode reaction for anodic dissolution was found to be kinetically controlled by dissociation of UC to a transition state complex that was hypothesized to consist of a uranium atom partially complexed by chloride ions. Precise measurements of current efficiencies using chronopotentiometry indicated upper limits of 90.9 ± 3.4% and 98.3 +1.7/-3.7% for anode and cathode, respectively, when operating at anodic overpotentials near +300 mV. Bulk electrolysis of a UC pellet performed by passing 98% of the theoretical charge resulted in nearly complete recovery of its uranium content as highly pure metal at the cathode.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Fission yeast Duf89 and Duf8901 are cobalt/nickel-dependent phosphatase–pyrophosphatases that act via a covalent aspartyl–phosphate intermediate

Domain of Unknown Function 89 (DUF89) proteins are metal-dependent phosphohydrolases. Exemplary DUF89 enzymes differ in their metal and phosphosubstrate preferences. Here, we interrogated the activities and structures of two DUF89 paralogs from fission yeast—Duf89 and Duf8901. We find that Duf89 and Duf8901 are cobalt/nickel-dependent phosphohydrolases adept at hydrolyzing p-nitrophenylphosphate and PP i . Crystal structures of metal-free Duf89 and Co 2+ -bound Duf8901 disclosed two enzyme conformations that differed with respect to the position of a three-helix module, which is either oriented away from the active site in Duf89 or forms a lid over the active site in Duf8901. Lid closure results in a 16 Å movement of Duf8901 Asp195, vis-à-vis Asp199 in Duf89, that brings Asp195 into contact with an octahedrally coordinated cobalt. Reaction of Duf8901 with BeCl 2 and NaF in the presence of divalent cations Co 2+ , Ni 2+ , or Zn 2+ generated covalent Duf8901-(Asp248)–beryllium trifluoride (BeF 3 )•Co 2+ , Duf8901-(Asp248)–BeF 3 •Ni 2+ , or Duf8901-(Asp248)–BeF 3 •Zn 2+ adducts, the structures of which suggest a two-step catalytic mechanism via formation and hydrolysis of an enzyme-(aspartyl)–phosphate intermediate. Alanine mutations of Duf8901 Asp248, Asn249, Lys401, Asp286, and Asp195 that interact with BeF 3 •Co 2+ squelched p-nitrophenylphosphatase activity. A 1.8 Å structure of a Duf8901-(Asp248)–AlF 4 –OH 2 •Co 2+ transition-state mimetic suggests an associative mechanism in which Asp195 and Asp363 orient and activate the water nucleophile. Whereas deletion of the duf89 gene elicited a phenotype in which expression of phosphate homeostasis gene pho1 was derepressed, deleting duf8901 did not, thereby hinting that the DUF89 paralogs have distinct functional repertoires in vivo.

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

Acetone to isobutene conversion on Zn x Ti y O z : Effects of TiO 2 facet

In this study, liquid-phase chemical grafting method was used to graft Zn onto TiO 2 with preferentially exposed (1 0 1) or (0 0 1) facet. The obtained Zn x Ti y O z materials were characterized using various techniques (e.g., XRD, Raman, DRIFTS etc.) and evaluated for the acetone-to-isobutene reaction. It was found that over TiO 2 (0 0 1), both terminal and bridging hydroxyls were readily titrated by Zn deposition, whereas a substantial amount of bridging hydroxyls on TiO 2 (1 0 1) remained. Although dominant Zn-O-terminated surface was obtained on two Zn x Ti y O z samples, bridging hydroxyls with high H-D exchange reactivity were observed on Zn x Ti y O z (1 0 1) compared with Zn x Ti y O z (0 0 1). The bridging hydroxyls showing rapid proton transfer efficiently stabilizes a transition state of diacetone alcohol intramolecular rearrangement for isobutene production as opposed to the diacetone alcohol dehydration.

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

Understanding the influence of solvents on the Pt-catalyzed hydrodeoxygenation of guaiacol

Bio-oil derived from fast pyrolysis of lignocellulosic biomass needs to be deoxygenated to become a substitute for petroleum fuels. Here, we study the hydrodeoxygenation mechanism of guaiacol, a bio-oil model compound derived from the lignin fraction of biomass, on Pt(1 1 1) terrace sites in the presence of water, diethyl ether, 1-butanol, and n-hexane as solvent. Using first-principles periodic density functional theory (DFT) calculations and mean-field microkinetic reactor modeling, a detailed reaction mechanism is investigated targeting various products such as catechol, phenol, anisole, benzene, cyclohexanone, and cyclohexanol. Solvent phase DFT outcomes are mostly similar to that of the vapor phase; however, microkinetic modeling results suggest that rate controlling species and transition states differ somewhat in the various reaction environments. Catechol was found to be the major aromatic product across all reaction environments. Over Pt(1 1 1), unsaturated monooxygenate production from catechol is unlikely.

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

A constitutive framework for rocks undergoing solid dissolution

Here, we formulate a time-dependent damage theory for rocks subjected to mechanical deformation and solid dissolution. The constitutive description is inspired by the transition state theory, which states that the rate of dissolution is a function of the reactive surface area measured through the crack density in the volume. We use a gradient-enhanced damage framework in which damage depends on the deformation of the material as well as on the amount of solid mass dissolved over time. The gradient-enhanced formulation is characterized by a three-field variational formulation with the solid displacement, nonlocal equivalent strain, and nonlocal rate of solid dissolution as the basic state variables. Traditionally, time-independent damage theories have only allowed damage to increase with increasing external load. In the proposed framework, the degree of damage may increase due to solid dissolution even when the external load is held fixed. In this way, solid dissolution is viewed as a process that is responsible for bringing about rate-dependent effects such as creep and stress-relaxation, which are two common features of geomaterial behavior.

42 ENGINEERING↗
DOE OSTI2025

Flexibility of oxygen sublattice and hydrogen bond length predict proton mobility in ternary metal oxides

Discovery of fast proton conductors is important for advancing clean energy technologies. This requires a better understanding of proton migration mechanisms. While structural and chemical traits of ternary metal oxides have been related to proton migration barriers, lattice dynamical effects have not been resolved quantitatively. Here, in this work, we introduce a phonon-based dynamic descriptor, termed ‘‘thermal O…O fluctuation,’’ quantifying the flexibility of donor-acceptor oxide-ion pairs. This enables direct comparison of O-sublattice flexibility across diverse metal oxides. Using regression models, we ranked physical descriptors as predictors of proton mobility, finding that H-bond length and thermal O…O fluctuation were the strongest descriptors. Further analysis revealed a critical O…O spacing of 2.4 A˚ at the transition state, which is easier to reach by more flexible donor-acceptor pairs, enabling facile proton transfer. Our results demonstrate oxygen sublattice flexibility as a dynamic descriptor and provide guiding principles for enhancing proton mobility in ternary metal oxides.

diffusivity↗
DOE OSTI2026

Synthesis and application of thermally responsive nanofiber coatings for overtemperature monitoring

This study presents a one-pot synthesis route to organometallic nanofibers based on copper thiolate, exhibiting distinctive chemical and physical characteristics. Electron microscopy analysis of morphology and composition revealed 2-10 μm-long, 50-90 nm-diameter hollow and non-hollow fibers composed of copper, sulfur, oxygen, hydrocarbon, and chlorine. Thermogravimetric analysis showed a pronounced mass loss within 120°C-135°C. To elucidate the thermal responsive pathways, the nanofibers were characterized before and after heating. X-ray photoelectron spectroscopy indicates that an initially mixed Cu(I)/Cu(II) oxidation states transition to predominantly Cu(I) upon heating. A layer of nanofiber was coated on battery pouch foil and evaluated as a candidate thermally sensitive coating. At elevated temperature (100-130°C), nanofiber coating released volatile organic compounds, sulfide and sulfur dioxide as detected using multiple gas sensors. This thermally responsive gas release/sensing approach provides a potential large-area temperature monitoring strategy, which is particularly relevant where direct temperature measurements of individual batteries is impractical. The results established proof of concept for nanofiber-coated battery pouch foil as overtemperature warning platform that can provide alerts when surface temperatures exceed a critical threshold. More broadly, the ability to form interconnected fiber networks positions copper thiolate nanofiber coatings as promising materials for advanced applications.

Ihala Gamaralalage, Chanaka [ORNL] (ORCID:00000002↗
DOE OSTI2022

Experimental and modeling study of water time histories during H 2 S-N 2 O combustion in a shock tube

Here, with the interest in directly burning sour gas in gas turbines, and the fact that even small amounts of H 2 S or its combustion products can alter combustion characteristics, many research studies have been performed to better understand the combustion chemistry of H 2 S. In the present study, the water formation was followed by laser absorption with N 2 O as an oxidant, instead of O 2 . Nitrous Oxide being essentially consumed via N 2 O (+M) ⇌ N 2 + O (+M), the water formation via the H 2 S + O route can then be probed to further validate the models. Three H 2 S/N 2 O mixtures diluted in 98% Ar were studied to cover the following range of equivalence ratios: 0.5, 1.0, and 2.0, over a wide range of temperatures (1580–1940 K) around atmospheric pressure. A chemical kinetic model was then developed first by validating the base N 2 O kinetics mechanism, then by investigating the experimental results presented herein. The N 2 O kinetics mechanism was updated with recent work on the low- and high-pressure limits for N 2 O decomposition (N 2 O (+M) ⇌ N 2 + O (+M)) as well as a review of rate coefficients for N 2 O + H ⇌ N 2 + OH from the literature. Good agreement is shown for H 2 /N 2 O mixtures. Updates were then made to the H 2 S kinetics mechanism, specifically, an update from the literature on SO 2 + H (+M) ⇌ SO + OH (+M) and an adjustment to SO + SH ⇌ S 2 + OH. Additionally, reactions between SH and N 2 O were determined using W1BD data and transition state theory which required the addition of an NNS sub-mechanism. With these updates, the mechanism provides good agreement with the H 2 S/N 2 O experiments.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

From electronic structure to model application of key reactions for gasoline/alcohol combustion: Hydrogen-atom abstraction by $CH_3O\dot{O}$ radicals

Hydrogen atom abstraction by methyl peroxy ($CH_3O\dot{O}$) radicals can play an important role in gasoline/ethanol interacting chemistry for fuels that produce high concentrations of methyl radicals. Detailed kinetic reactions for hydrogen atom abstraction by $CH_3O\dot{O}$ radicals from the components of FGF-LLNL (a gasoline surrogate) including cyclopentane, toluene, 1-hexene, n -heptane, and isooctane have been systematically studied in this work. Here, the M06-2X/6-311 ++ G(d,p) level of theory was used to obtain the optimized structure and vibrational frequency for all stationary points and the low-frequency torsional modes. The 1-D hindered rotor treatment for low-frequency torsional modes was treated at M06-2X/6-31G level of theory. The UCCSD(T)-F12a/cc-pVDZ-F12 and QCISD(T)/CBS level of theory were used to calculate single point energies for all species. High pressure limiting rate constants for all hydrogen atom abstraction channels were performed using conventional transition state theory with unsymmetric tunneling corrections. Individual rate constants are reported in the temperature range from 298.15 to 2000 K. Our computed results show that the abstraction of allylic hydrogen atoms from 1-hexene is the fastest at low temperatures. When the temperature increases, the hydrogen atom abstraction reaction channel at the primary alkyl site gradually becomes dominant. Thermodynamics properties for all stable species and high-pressure limiting rate constants for each reaction pathway obtained in this work were incorporated into the latest gasoline surrogate/ethanol model to investigate the influence of the rate constants calculated here on model predicted ignition delay times.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2026

Pressure-dependent kinetics of cyclopentene consumption

Recent experiments have uncovered deficiencies in the ability of current kinetic models for cyclopentane combustion to predict species profiles as a function of O2 concentration. One likely source of these discrepancies flows from the limited availability of relevant kinetic data for key functionalized intermediates, requiring the use of distant analogies to estimate their rates. This work fills this gap for cyclopentene, the most prominent C5 intermediate in low-temperature cyclopentane combustion. Ab initio transition-state-theory based master equation methods are used to calculate pressure-dependent rate constants across four potential energy surfaces (C5H8 + OH, C5H8 + HO2, C5H7, and C5H7 + O2). Potential energy surfaces are characterized at the CCSD(T)-F12/cc-pVDZ-F12//revDSD-PBEP86-D3BJ/def2-TZVPP level of theory. Pressure-dependent rate constants and branching fractions are determined via 1D master equation calculations. The results are used to analyze the competition for radicals between cyclopentene and cyclopentane in the intermediate temperature region, and to identify the dominant pathways flowing from these initiation channels. On the OH initiation surface, we find good agreement with recently published experimental rate constants at high temperature, and we elucidate the complex temperature dependence below 800 K, where 𝜋𝜋 addition becomes prominent. On the HO2 initiation surface, we find an interesting alternative route to bicyclic ether formation via well-skipping over the HO2 adduct well. On the cyclopentenyl radical surface, the major product is cyclopentadiene plus H. On the cyclopentenyl + O2 surface, we find that resonance stabilization largely prevents the allylic radical from undergoing oxidation, whereas the alkylic radical has a somewhat elevated propensity for oxidation pathways relative to cyclopentyl.

ab initio calculations↗
DOE OSTI2021

On the structure sensitivity of and CO coverage effects on formic acid decomposition on Pd surfaces

Using density functional theory calculations, the Pd-catalyzed vapor-phase formic acid decomposition was studied, with a focus on the structure sensitivity and CO coverage effects. A comprehensive reaction network was developed on both the (111) and (100) facets of Pd, at CO coverages of 0 and 5/9 monolayer (ML). Pd(100) was determined to be more reactive than Pd(111) at both CO coverages. The introduction of 5/9 ML CO decreased the activity of both facets significantly, due to destabilization of the surface intermediates and transition states on the CO-decorated surfaces. Three reaction pathways were explored on the clean surfaces: the formate (HCOO) pathway, the carboxyl (COOH) pathway leading to the formation of CO 2 , and the COOH pathway leading to the formation of CO (COOH→CO). Based on the DFT-derived energetics alone, it appears that all three pathways contribute to the reaction on clean Pd, whereas the presence of 5/9 ML of CO inhibits the HCOO pathway on both facets and favors the COOH→CO pathway on the (111) facet, but the COOH→CO 2 one on the (100) facet. Moreover, at high CO coverages, alternative spectator CO-assisted adsorbate decomposition pathways were discovered, which could potentially play a role in formic acid decomposition on Pd catalysts under realistic reaction conditions.

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

Effects of co-adsorbed water on different bond cleavages involved in acetic acid decomposition on Pt (111)

Acetic acid decomposition on Pt (111) in the presence of co-adsorbed water is a good model system for oxygenate decomposition on Pt (111) in aqueous phase, with application in hydrogen production and biomass conversion. Here, in this study, we present a density functional theory (DFT) theory calculation of how co-adsorbed water affects different bond cleavages of acetic acid decomposition on Pt (111). The presence of co-adsorbed water generally enhances O$-$H bond cleavage while inhibiting OC$-$O and OC$-$OH bond cleavage. The influence of co-adsorbed water on C$-$H bond cleavage varies the most and depends on the nature of the transition state and how co-adsorbed water stabilizes the initial and final state. Although these trends are useful as general guidance, they are not sufficient to predict the effect on a complex reaction network such as acetic acid decomposition on Pt (111) which has several parallel reaction paths with similar energies. In the absence of co-adsorbed water, the two lowest energy pathways are decarboxylation (DCX) and decarbonylation (DCN) pathways through a common CH 2 COO intermediate, in which the DCX pathway (C$-$H bond cleavage of CH 2 COO) is more favorable than DCN pathway (OC$-$O bond cleavage of CH 2 COO). In the presence of co-adsorbed water, the energy difference between Csingle bondH bond cleavage and OC$-$O bond cleavage of CH 2 COO increases, suggesting an increase in the favorability of acetic acid decarboxylation (formation of carbon dioxide) over acetic acid decarbonylation (formation of carbon monoxide) on Pt (111).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗