Molecular-level descriptions of surface chemistry in kinetic models using density functional theory

On the mechanism of low-temperature water gas shift reaction on copper

Periodic, self-consistent density functional theory (DFT-GGA) calculations are used to investigate the water gas shift reaction (WGSR) mechanism on Cu(111). The thermochemistry and activation energy barriers for all the elementary steps of the commonly accepted redox mechanism, involving complete water activation to atomic oxygen, are presented. Through our calculations, we identify carboxyl, a new reactive intermediate, which plays a central role in WGSR on Cu(111). The thermochemistry and activation energy barriers of the elementary steps of a new reaction path, involving carboxyl, are studied. A detailed DFT-based microkinetic model of experimental reaction rates, accounting for both the previous and the new WGSR mechanism show that, under relevant experimental conditions, (1) the carboxyl-mediated route is the dominant path, and (2) the initial hydrogen abstraction from water is the rate-limiting step. Formate is a stable "spectator" species, formed predominantly through CO2 hydrogenation. In addition, the microkinetic model allows for predictions of (i) surface coverage of intermediates, (ii) WGSR apparent activation energy, and (iii) reaction orders with respect to CO, H2O, CO2, and H2.

Methane oxidation mechanism on Pt(111): a cluster model DFT study

The electronic energy barriers of surface reactions pertaining to the mechanism of the electrooxidation of methane on Pt (111) were estimated with density functional theory calculations on a 10-atom Pt cluster, using both the B3LYP and PW91 functionals. Optimizations of initial and transition states were performed for elementary steps that involve the conversion of CH(4) to adsorbed CO at the Pt/vacuum interface. As a first approximation we do not include electrolyte effects in our model. The reactions include the dissociative chemisorption of CH(4) on Pt, dehydrogenation reactions of adsorbed intermediates (*CH(x) --> *CH(x-1) + *H and *CH(x)O --> *CH(x-1)O + *H), and oxygenation reactions of adsorbed CH(x) species (*CH(x) + *OH --> *CH(x)OH). Many pathways were investigated and it was found that the main reaction pathway is CH(4) --> *CH(3) --> *CH(2) --> *CH --> *CHOH --> *CHO --> *CO. Frequency analysis and transition-state theory were employed to show that the methane chemisorption elementary step is rate-limiting in the above mechanism. This conclusion is in agreement with published experimental electrochemical studies of methane oxidation on platinum catalysts that have shown the absence of an organic adlayer at electrode potentials that allow the oxidation of adsorbed CO. The mechanism of the electrooxidation of methane on Pt is discussed.

Lattice strain effects on CO oxidation on Pt(111)

Surface strain plays a major role in determining the rate limiting step and catalytic activity of platinum for CO oxidation.

Ethene Epoxidation Selectivity Inhibited by Twisted Oxametallacycle: A DFT Study on Ag Surface-Oxide

Competitive ethene oxidation pathways are presented for a p(4 x 4) surface-oxide phase on Ag(111) obtained from density functional theory (DFT) calculations. Both parallel routes are found to proceed from a common oxametallacycle intermediate (OMME) in agreement with previous mechanistic studies on low coverage O adatom phase, although acetaldehyde (AcH) is favored by almost 2 kcal/mol. An even more striking difference with pure metal surface appears with the oxide regeneration pathways, which are found non-rate controlling. Furthermore, a kinetic model is developed on the basis of these DFT calculations and yields 96% selectivity in favor of AcH for a simulation in realistic catalytic conditions (600 K and respective partial pressures of 1 atm for ethene and oxygen reactants). As a key finding, this low ethene epoxide selectivity is proposed to be directly linked to the conformational barrier of the pivotal intermediate. In fact, the elasticity of the ultrathin oxide adlayer enables a twisted OMME structure as a true minimum, which agrees well with orbital prerequisite of the concerted H migration toward AcH. On the contrary, the desired selective ring closure forming ethene epoxide (EO) requires conformational inversion although the eclipsed form lies 2 kcal/mol above.

Hierarchical Multiscale Mechanism Development for Methane Partial Oxidation and Reforming and for Thermal Decomposition of Oxygenates on Rh

A thermodynamically consistent C1 microkinetic model is developed for methane partial oxidation and reforming and for oxygenate (methanol and formaldehyde) decomposition on Rh via a hierarchical multiscale methodology. Sensitivity analysis is employed to identify the important parameters of the semiempirical unity bond index quadratic exponential potential (UBI-QEP) method and these parameters are refined using quantum mechanical density functional theory. With adjustment of only two pre-exponentials in the CH4 oxidation subset, the C1 mechanism captures a multitude of catalytic partial oxidation (CPOX) and reforming experimental data as well as thermal decomposition of methanol and formaldehyde. We validate the microkinetic model against high-pressure, spatially resolved CPOX experimental data. Distinct oxidation and reforming zones are predicted to exist, in agreement with experiments, suggesting that hydrogen is produced from reforming of methane by H2O formed in the oxidation zone. CO is produced catalytically by partial oxidation up to moderately high pressures, with water-gas shift taking place in the gas-phase at sufficiently high pressures resulting in reduction of CO selectivity.

In Situ ATR-IR Spectroscopic and Reaction Kinetics Studies of Water−Gas Shift and Methanol Reforming on Pt/Al2O3 Catalysts in Vapor and Liquid Phases

Reaction kinetics measurements of the water-gas shift reaction were carried out at 373 K on Pt/Al2O3 in vapor phase to investigate the effects of CO, H2, and H2O partial pressures. Results of in situ ATR-IR studies conducted in vapor phase under similar conditions suggest that the Pt surface coverage by adsorbed CO is high (approximately 90% of the saturation coverage), leading to a negligible effect of the CO pressures on the rate of reaction. The negative reaction order with respect to the H2 pressure is caused by the increased coverage of adsorbed H atoms, and the fractional positive order with respect to the water pressure is consistent with non-equilibrated H2O dissociation on Pt. Results of in situ ATR-IR studies carried out at 373 K show that the presence of liquid water leads to a slight decrease in the Pt surface coverage by adsorbed CO (approximately 80% of the saturation coverage) when the CO partial pressure is the same as in the vapor-phase studies. The rate of the WGS reaction in the presence of liquid water is comparable to the rate under complete vaporization conditions when other factors (such as CO partial pressure) are held constant. Reaction kinetics measurements of methanol reforming were carried out at 423 K over a total pressure range of 1.36-5.84 bar. In situ ATR-IR studies were conducted at 423 K to determine the Pt surface coverage by adsorbed CO in completely vaporized methanol feeds and in aqueous methanol solutions. The decomposition of methanol is found to be slower during the reforming of methanol in liquid phase than in vapor phase, which leads to a lower rate of hydrogen production in liquid phase (0.08 min(-1) at 4.88 bar) than in vapor phase (0.23 min(-1) at 4.46 bar). The lower reaction order with respect to methanol concentration observed for vapor-phase versus liquid-phase methanol reforming (0.2 versus 0.8, respectively) is due to the higher extent of CO poisoning on Pt for reforming in vapor phase than in liquid phase, based on the higher coverage by adsorbed CO observed in completely vaporized methanol feeds (55-60% of the saturation coverage) than in aqueous methanol feed solutions (29-40% of the saturation coverage).