Monte Carlo simulation of catalytic CO oxidation

A simple Monte Carlo model of the CO oxidation on a single-crystal catalyst surface is presented. The simulation model considers the following elementary reaction steps:

1. (1) chemisorption of a CO molecule, its surface migration and possible desorption

2. (2) physisorption of an O₂ molecule to a precursor state and its subsequent dissociative chemisorption

3. (3) activated reaction of adsorbed O and CO (the Langmuir - Hinshelwood reaction mechanism), formation of CO₂ and its rapid desorption.

The changes in the activation energy of reaction and in the adsorption energy of CO resulting from the interactions between adsorbed species are also considered. The model makes possible to monitor temperature programmed reaction spectra or reaction spectra obtained during changes of the ratio of the partial pressures of CO and O₂. The results of simulations for a Pd(111) single-crystal plane are compared with experiment.     

Oxidation of CO by SO2: a theoretical study

The elementary reaction SO(2) + CO --> CO(2) + SO((3)Sigma) (1) and the subsequent reaction SO((3)Sigma) + CO --> CO(2) + S((3)P) (2) have been studied by the application of the Gaussian-3//B3LYP quantum chemical approach to characterize the potential energy surfaces and transition state kinetic analysis to derive rate coefficients. Reaction 1 is found to take place via two transition states (TS), a cis-OSOCO TS and a trans-OSOCO TS. Reaction via the cis-TS is concerted and takes place on a singlet surface. Intersystem crossing to the final products occurs after passage through the barrier on the singlet surface. The trans-TS leads to a very weakly bound singlet OSOCO intermediate that then passes through a second TS (on the triplet surface) to form the products. Reaction 2 takes place on triplet surfaces. There is a concerted reaction through a cis-SOCO TS and a weakly bound trans-SOCO has also been identified. Reaction 2 is analogous to the reaction CO + O(2)((3)Sigma) --> CO(2) + O((3)P) (3), and this reaction has been reinvestigated at a similar level of theory and the rate coefficient derived by quantum chemistry is compared with experiment. The sensitive effects of trace impurities such as H(2), H(2)O, and hydrocarbons on the accurate experimental determination of the rate coefficient of reaction 3 is discussed. Using rate coefficients for reactions 1 and 2 obtained via quantum chemical calculations, we have been unable to model the extent of decomposition of SO(2) measured in a shock tube study of reaction between SO(2) and CO [Bauer, S. H.; Jeffers, P.; Lifshitz, A.; Yadava, B. P. Proc. Combust. Inst. 1971, 13, 417]. In light of the known sensitivity of reaction 3 to trace impurities, we have incorporated trace amounts of H(2), CH(4), or H(2)O, together with our rate coefficients for (1) and (2), in a kinetic model of Alzueta et al. [Combust. Flame 2001, 127, 2234], which is then shown to be able to substantially model the SO(2) data of Bauer et al. In the course of this modeling study we also computed heats of formation for a number of sulfur-containing small molecules: HS, HSO, HSOH, HOSO, HS(2), HSO(2), HOSO(2), HOSOH, and HOSHO.     

Competing mechanisms in the optically activated functionalization of the hydrogen-terminated Si(111) surface

Recent experiments have shown that organic monolayers on silicon surfaces can be formed through the optically activated surface reaction of H-terminated Si surfaces with terminally unsaturated organic molecules (Eves et al. J. Am. Chem. Soc. 2004, 126, 14318; Sun et al. J. Am. Chem. Soc. 2005, 127, 2514). Possible mechanisms for the formation of this monolayer involve the abstraction of a H atom either at the same attachment site of the molecule (Path A) or from a neighboring site (Path B). Using periodic Density Functional Theory calculations together with an efficient method for finding reaction pathways, we examine both optically activated reaction mechanisms for an alkene and an aldehyde reacting with H-Si(111). Our results show that while Path A is energetically more favorable its significant barrier is likely to limit its viability. Path B on the other hand encounters a much lower H atom abstraction barrier and appears to be more viable.     

Experimental and theoretical study of reactivity trends for methanol on Co/Pt(111) and Ni/Pt(111) bimetallic surfaces

Methanol was used as a probe molecule to examine the reforming activity of oxygenates on NiPt(111) and CoPt(111) bimetallic surfaces, utilizing density functional theory (DFT) modeling, temperature-programmed desorption, and high-resolution electron energy loss spectroscopy (HREELS). DFT results revealed a correlation between the methanol and methoxy binding energies and the surface d-band center of various NiPt(111) and CoPt(111) bimetallic surfaces. Consistent with DFT predictions, increased production of H2 and CO from methanol was observed on a Ni surface monolayer on Pt(111), designated as Ni-Pt-Pt(111), as compared to the subsurface monolayer Pt-Ni-Pt(111) surface. HREELS was used to verify the presence and subsequent decomposition of methoxy intermediates on NiPt(111) and CoPt(111) bimetallic surfaces. On Ni-Pt-Pt(111) the methoxy species decomposed to a formaldehyde intermediate below 300 K; this species reacted at approximately 300 K to form CO and H2. On Co-Pt-Pt(111), methoxy was stable up to approximately 350 K and decomposed to form CO and H2. Overall, trends in methanol reactivity on NiPt(111) bimetallic surfaces were similar to those previously determined for ethanol and ethylene glycol.     

Selective solar-driven reduction of CO2 to methanol using a catalyzed p-GaP based photoelectrochemical cell

With rising atmospheric CO2 levels, there has been increasing interest in artificial photosynthetic schemes for converting this greenhouse gas into valuable fuels and small organics. Photoelectrochemical schemes for activating the inert CO2 molecule, however, operate at excessive overpotentials and thus do not convert actual light energy to chemical energy. Here we describe the selective conversion of CO2 to methanol at a p-GaP semiconductor electrode with a homogeneous pyridinium ion catalyst, driving the reaction with light energy to yield faradaic efficiencies near 100% at potentials well below the standard potential.     

单原子Ge助剂修饰Cu(111)晶面上CO2加氢制甲醇的机理研究

针对CO₂热催化转化制甲醇过程中CO₂吸附、活化较困难及副产物较多的问题,提出采用单原子Ge助剂修饰Cu(111)晶面的解决思路,通过密度泛函理论(DFT)计算研究了CO₂在Ge-Cu(111)晶面上加氢合成甲醇的反应机理。结果表明,单原子Ge助剂的电子调控增加了与其相邻的 Cu 原子的电子云密度,使 CO₂分子在含 Ge 活性界面上的吸附能力显著增强:CO₂在 Ge-Cu(111)晶面上的吸附能约为Cu(111)晶面的1.5倍,约为Pd改性Cu(111)晶面的2.4倍,进而使逆水煤气变换(RWGS)反应路径速控步骤的活化能降低了近 20 kJ·mol^-1,同时衍生出 3条生成甲醇的 RWGS新路径;此外,Ge-Cu(111)晶面上甲酸盐路径由于速控步骤活化能大幅上升而被禁阻,进而CO及烃类等副产物选择性大幅降低,Ge-Cu(111)晶面上CO₂加氢制甲醇选择性升高。     

单原子Ge助剂修饰Cu(111)晶面上CO2加氢制甲醇的机理研究

针对CO₂热催化转化制甲醇过程中CO₂吸附、活化较困难及副产物较多的问题,提出采用单原子Ge助剂修饰Cu (111)晶面的解决思路,通过密度泛函理论(DFT)计算研究了CO₂在Ge-Cu(111)晶面上加氢合成甲醇的反应机理。结果表明,单原子Ge助剂的电子调控增加了与其相邻的Cu原子的电子云密度,使CO₂分子在含Ge活性界面上的吸附能力显著增强：CO₂在Ge-Cu(111)晶面上的吸附能约为Cu (111)晶面的1.5倍,约为Pd改性Cu(111)晶面的2.4倍,进而使逆水煤气变换(RWGS)反应路径速控步骤的活化能降低了近20 kJ·mol^-1,同时衍生出3条生成甲醇的RWGS新路径;此外,Ge-Cu(111)晶面上甲酸盐路径由于速控步骤活化能大幅上升而被禁阻,进而CO及烃类等副产物选择性大幅降低,Ge-Cu(111)晶面上CO₂加氢制甲醇选择性升高。     

Ab initio design of picosecond infrared laser pulses for controlling vibrational-rotational excitation of CO molecules

Optimal control of rovibrational excitations of the CO molecule using picosecond infrared laser pulses is described in the framework of the electric-nuclear Born-Oppenheimer approximation [G. G. Balint-Kurti et al., J. Chem. Phys. 122, 084110 (2005)]. The potential energy surface of the CO molecule in the presence of an electric field is calculated using coupled cluster theory with a large orbital basis set. The quantum dynamics of the process is treated using a full three dimensional treatment of the molecule in the laser field. The detailed mechanisms leading to efficient control of the selected excitation processes are discussed.     

Relativistic DFT study on the reaction mechanism of second-row transition metal Ru with CO2

To estimate the importance of relativistic effects on the reaction mechanisms between Ru and CO2, the potential energy surfaces have been performed in the triplet and quintet electronic states using quasi-relativistic (Pauli), zero-order regularly approximated (ZORA), and nonrelativistic (NR) density functional theory (DFT) at the PW91/TZP level. The results demonstrate that there are two rival reaction mechanisms: one is an addition mechanism and the other is an insertion mechanism in the triplet state. The only mechanism in the quintet state is the insertion mechanism. The most favored reaction mechanism in Ru + CO2 is that the Ru atom in its ground state first attacks the CO bond of CO2, forming q-Ru(CO)O (5A') with the insertion mechanism, and then undergoes an intersystem crossing to t-Ru(CO)O (3A'). Then it crosses t-TS3 to produce t-ORuCO molecule. The relativistic effects are important for reactivity of the second-row transition metal to CO2. In the key step of t-Ru(CO)O via t-TS3 to t-ORuCO, relativistic effects reduce the barrier energy by 10.3 kcal/mol, which is nearly half the nonrelativistic barrier energy.     

From CO2 to dimethyl carbonate with dialkyldimethoxystannanes: the key role of monomeric species

The formation of dimethyl carbonate (DMC) from CO(2) and methanol with the dimer [n-Bu(2)Sn(OCH(3))(2)](2) was investigated by experimental kinetics in support of DFT calculations. Under the reaction conditions (357-423 K, 10-20 MPa), identical initial rates are observed with three different reacting mixtures, CO(2)/toluene, supercritical CO(2), and CO(2)/methanol, and are consistent with the formation of monomeric di-n-butyltin(iv) species. An intramolecular mechanism is, therefore, proposed with an Arrhenius activation energy amounting to 104 ± 10 kJ mol(-1) for DMC synthesis. DFT calculations on the [(CH(3))(2)Sn(OCH(3))(2)](2)/CO(2) system show that the exothermic insertion of CO(2) into the Sn-OCH(3) bond occurs by a concerted Lewis acid-base interaction involving the tin center and the oxygen atom of the methoxy ligand. The Gibbs energy diagrams highlight that, under the reaction conditions, the dimer-monomer equilibrium is significantly shifted towards monomeric species, in agreement with the experimental kinetics. Importantly, the two Sn-OCH(3) bonds are prompt to insert CO(2). These results provide new insight into the reaction mechanism and catalyst design to enhance the turnover numbers.     

Examination of the photophysical processes of chlorophyll d leading to a clarification of proposed uphill energy transfer processes in cells of Acaryochloris marinas

A comprehensive study of the photophysical properties of chlorophyll (Chl) d in 1:40 acetonitrile-methanol solution is performed over the temperature range 170-295 K. From comparison of absorption and emission spectra, time-dependent density-functional calculations and homologies with those of Chl a, we assign the key features of the absorption and fluorescence spectra. Possible photophysical energy relaxation mechanisms are summarized, and thermal equilibration processes are studied in detail by monitoring the observed emission profiles and quantum yields as a function of excitation energy. In particular, we concentrate on emission subsequent to excitation in the extreme far-red tail of the Qy absorption spectrum, with this emission partitioned into contributions from hot-band absorptions as well as uphill energy transfer processes that occur subsequent to absorption. No unusual photophysical processes are detected for Chl d; it appears that all intramolecular relaxation processes reach thermal equilibration on shorter timescales than the fluorescence lifetime even at 170 K. The results from these studies are used to reinterpret a previous study of photochemical processes observed in intact cells and their acetone extracts of the photosynthetic system of Acaryochloris marina. In the study of Mimuro et al., light absorbed by Chl d at 736 nm is found to give rise to emission by another species, believed to also be Chl d, at 703 nm; this uphill energy transfer process is easily rationalized in terms of the thermal equilibration processes that we deduced for Chl d. However, no evidence is found in the experimental results of Mimuro et al. to support claims that (nonequilibrium) uphill energy transfer is additionally observed to Chl a species that emit at 670-680 nm. This finding is relevant to broader issues concerning the nature of the special pair in photosystem II of A. marina because suggestions that it is comprised of Chl a can only be correct if nonthermal uphill energy transfer processes from Chl d are operative.     

A comprehensive kinetic mechanism for CO,CH2O,and CH3OH combustion

New experimental profiles of stable species concentrations are reported for formaldehyde oxidation in a variable pressure flow reactor at initial temperatures of 850–950 K and at constant pressures ranging from 1.5 to 6.0 atm. These data, along with other data published in the literature and a previous comprehensive chemical kinetic model for methanol oxidation, are used to hierarchically develop an updated mechanism for CO/H₂O/H₂/O₂, CH₂O, and CH₃OH oxidation. Important modifications include recent revisions for the hydrogen–oxygen submechanism (Li et al., Int J Chem Kinet 2004, 36, 565), an updated submechanism for methanol reactions, and kinetic and thermochemical parameter modifications based upon recently published information. New rate constant correlations are recommended for CO + OH = CO₂ + H ( R23 ) and HCO + M = H + CO + M ( R24 ), motivated by a new identification of the temperatures over which these rate constants most affect laminar flame speed predictions (Zhao et al., Int J Chem Kinet 2005, 37, 282). The new weighted least‐squares fit of literature experimental data for ( R23 ) yields k₂₃ = 2.23 × 10⁵T^1.89exp(583/T) cm³/mol/s and reflects significantly lower rate constant values at low and intermediate temperatures in comparison to another recently recommended correlation and theoretical predictions. The weighted least‐squares fit of literature results for ( R24 ) yields k₂₄ = 4.75 × 10¹¹T^0.66exp(?7485/T) cm³/mol/s, which predicts values within uncertainties of both prior and new (Friedrichs et al., Phys Chem Chem Phys 2002, 4, 5778; DeSain et al., Chem Phys Lett 2001, 347, 79) measurements. Use of either of the data correlations reported in Friedrichs et al. (2002) and DeSain et al. (2001) for this reaction significantly degrades laminar flame speed predictions for oxygenated fuels as well as for other hydrocarbons. The present C₁/O₂ mechanism compares favorably against a wide range of experimental conditions for laminar premixed flame speed, shock tube ignition delay, and flow reactor species time history data at each level of hierarchical development. Very good agreement of the model predictions with all of the experimental measurements is demonstrated. © 2007 Wiley Periodicals, Inc. 39: 109–136, 2007     

Active surfaces for CO oxidation on palladium in the hyperactive state

Hyperactivity was previously observed for CO oxidation over palladium, rhodium, and platinum surfaces under oxygen-rich conditions, characterized by reaction rates 2-3 orders higher than those observed under stoichiometric reaction conditions [Chen et al. Surf. Sci. 2007, 601, 5326]. In the present study, the formation of large amounts of CO(2) and the depletion of CO at the hyperactive state on both Pd(100) and polycrystalline Pd foil were evidenced by the infrared intensities of the gas phase CO(2) and CO, respectively. The active surfaces at the hyperactive state for palladium were characterized using infrared reflection absorption spectroscopy (IRAS, 450-4000 cm(-1)) under the realistic catalytic reaction condition. Palladium oxide on a Pd(100) surface was reduced eventually by CO at 450 K, and also under CO oxidation conditions at 450 K. In situ IRAS combined with isotopic (18)O(2) revealed that the active surfaces for CO oxidation on Pd(100) and Pd foil are not a palladium oxide at the hyperactive state and under oxygen-rich reaction conditions. The results demonstrate that a chemisorbed oxygen-rich surface of Pd is the active surface corresponding to the hyperactivity for CO oxidation on Pd. In the hyperactive region, the CO(2) formation rate is limited by the mass transfer of CO to the surface.     

A comprehensive and compact n-heptane oxidation model derived using chemical lumping

A detailed reaction mechanism for n-heptane oxidation has been compiled and subsequently simplified. The model is based on a kinetic model for C1-C4 fuel oxidation of Hoyermann et al. [Phys. Chem. Chem. Phys., 2004, 6, 3824] and a detailed mechanism for n-heptane oxidation developed by Curran et al. [Combust. Flame, 1998, 114, 149]. The generated mechanism is kept compact by limiting the application of the low temperature oxidation pathways to the fuel molecule. The first reaction steps and the complex low temperature paths in the oxidation process have been simplified and reorganized by linear chemical lumping. The reported procedure allows a decrease in number of species and reactions with only a minor loss of model accuracy. The simplified model is of very compact size and gives an advantageous starting point for further model reduction. By this chemically lumped general mechanism without further adjustments the large set of experimental data for the high and low temperature oxidation (ignition delay times, species concentration profiles, heat release and engine pressure profiles, flame speeds and flame structure data) for conditions ranging from very low to high temperatures (550-2300 K), very lean to extremely fuel rich (0.22 < phi < 3) mixtures and pressures between 1 and 42 bar is consistently described providing a basis for reliable predictions for future applications, (i) building reaction mechanisms for similar but chemically more complex fuels (e.g. iso-octane, n-decane,...) and (ii) calculating complex flow fields ("fluid dynamics") after further simplification with advanced reduction tools.     

Influence of the cluster dimensionality on the binding behavior of CO and O2 on Au13

We present an ab initio density functional theory study of the binding behavior of CO and O(2) molecules to two- and three-dimensional isomers of Au(13) in order to investigate the potential catalytic activity of this cluster towards low-temperature CO oxidation. First, we scanned the potential energy surface of Au(13) and studied the effect of spin-orbit coupling on the relative stabilities of the 21 isomers we identified. While spin-orbit coupling increases the stability of the three-dimensional more than the two-dimensional isomers, the ground state structure at 0 K remains planar. Second, we systematically studied the binding of CO and O(2) molecules onto the planar and three-dimensional structures lowest in energy. We find that the isomer dimensionality has little effect on the binding of CO to Au(13). O(2), on the other hand, binds significantly to the three-dimensional isomer only. The simultaneous binding of multiple CO molecules decreases the binding energy per molecule. Still, the CO binding remains stronger than the O(2) binding. We did not find a synergetic effect due to the co-adsorption of both molecular species. On the three-dimensional isomer, we find O(2) dissociation to be exothermic with an dissociation barrier of 1.44 eV.     

A DFT study on the mechanism of hydrosilylation of unsaturated compounds with neutral hydrido(hydrosilylene)tungsten complex

Recently, Tobita et al. reported stoichiometric hydrosilylation reactions of acetone and acetonitrile with neutral hydrido(hydrosilylene)tungsten complexes Cp'(CO)2(H)W=Si(H)[C(SiMe(3))(3)] (Cp' = Cp*, C(5)Me(4)Et). The mechanisms of the hydrosilylation reactions of unsaturated compounds (ketone and nitrile) with the tungsten complexes have been investigated with the B(3)LYP density functional theory method. Four possible reaction mechanisms were studied. The results of the calculations indicate that the hydrosilylation of acetone proceeds via a metal hydride migration mechanism proposed by Tobita et al., while the hydrosilylation of nitrile occurs through a silyl migration mechanism, analogous to the modified Chalk-Harrod mechanism. The [2(sigma)+2(pi)] additions of various CX (CX = C=O or CN) multiple bonds with the Si-H bonds in the neutral complexes have very high barriers although similar additions were found feasible in other related cationic complexes. All the hydrosilylation reactions studied here give stable tungsten-silylene or tungsten-silyl products, which are not easily converted into the starting hydrido(hydrosilylene)tungsten complexes when reacting with a hydrosilane substrate molecule. Therefore, we predict that hydrosilylation of acetonitrile and acetone catalyzed by these tungsten complexes is difficult to achieve.     

Oxidation of methanol by FeO2+ in water: DFT calculations in the gas phase and ab initio MD simulations in water solution

We investigate the mechanism of methanol oxidation to formaldehyde by ironoxido ([Fe(IV)O]2+), the alleged active intermediate in the Fenton reaction. The most likely reaction mechanisms are explored with density functional theory (DFT) calculations on microsolvated clusters in the gas phase and, for a selected set of mechanisms, with constrained Car-Parrinello molecular dynamics (CPMD) simulations in water solution. Helmholtz free energy differences are calculated using thermodynamic integration in a simulation box with 31 water molecules at 300 K. The mechanism of the reaction is investigated with an emphasis on whether FeO2+ attacks methanol at a C-H bond or at the O-H bond. We conclude that the most likely mechanism is attack by the oxido oxygen at the C-H bond ("direct CH mechanism"). We calculate an upper bound for the reaction Helmholtz free energy barrier in solution of 50 kJ/mol for the C-H hydrogen transfer, after which transfer of the O-H hydrogen proceeds spontaneously. An alternative mechanism, starting with coordination of methanol directly to Fe ("coordination OH mechanism"), cannot be ruled out, as it involves a reaction Helmholtz free energy barrier in solution of 44 +/- 10 kJ/mol. However, this coordination mechanism has the disadvantage of requiring a prior ligand substitution reaction, to replace a water ligand by methanol. Because of the strong acidity of [FeO(H2O)5]2+, we also investigate the effect of deprotonation of a first-shell water molecule. However, this is found to increase the barriers for all mechanisms.     

Applied reaction dynamics: efficient synthesis gas production via single collision partial oxidation of methane to CO on Rh111

Supersonic molecular beams have been used to determine the yield of CO from the partial oxidation of CH4 on a Rh111 catalytic substrate, CH4+12O2-->CO+2H2, as a function of beam kinetic energy. These experiments were done under ultrahigh vacuum conditions with concurrent molecular beams of O2 and CH4, ensuring that there was only a single collision for the CH4 to react with the surface. The fraction of CH4 converted is strongly dependent on the normal component of the incident beam's translational energy, and approaches unity for energies greater than approximately 1.3 eV. Comparison with a simplified model of the methane-Rh111 reactive potential gives insight into the barrier for methane dissociation. These results demonstrate the efficient conversion of methane to synthesis gas, CO+2H2, are of interest in hydrogen generation, and have the optimal stoichiometry for subsequent utilization in synthetic fuel production (Fischer-Tropsch or methanol synthesis). Moreover, under the reaction conditions explored, no CO2 was detected, i.e., the reaction proceeded with the production of very little, if any, unwanted greenhouse gas by-products. These findings demonstrate the efficacy of overcoming the limitations of purely thermal reaction mechanisms by coupling nonthermal mechanistic steps, leading to efficient C-H bond activation with subsequent thermal heterogeneous reactions.     

Catalytic involvement of CO2 in the mutagenesis caused by reactions of ONOO(-) with guanine

The catalytic role of CO2 in reactions of ONOO- with guanine, leading to the formation of the mutagenic species 8-oxoguanine (8-oxoG) and 8-nitroguanine anion (8-nitroG-), was investigated by considering the reactions of nitrosoperoxycarbonate anion (ONOOCO2-), an adduct of ONOO- and CO2, with guanine at the B3LYP/6-31G** and B3LYP/AUG-cc-pVDZ levels of density functional theory in gas phase. In order to study bulk solvent effect, single-point energy calculations in aqueous media were carried out for all the species occurring in the reactions at the B3LYP/AUG-cc-pVDZ level of theory, by use of the polarizable continuum model (PCM). Vibrational frequency analysis was performed, and zero-point-energy (ZPE)-corrected total energies and Gibbs free energy changes at 298.15 K were obtained. The genuineness of the calculated transition states was confirmed by visually examining the vibrational modes and also by intrinsic reaction coordinate (IRC) calculations. The reaction between ONOOCO2- and guanine occurring through four different mechanisms leads to the formation of 8-oxoG or its anion, while the reaction between the same two species occurring through a different scheme leads to the formation of 8-nitroG-. It has been shown that the presence of a water molecule along with ONOOCO2- would not affect the reaction mechanisms significantly. Structures of the reactant complexes, product complexes and barrier energies involved in the reactions reveal that CO2 acts as a catalyst for the reaction between ONOO- and guanine. The cause of the catalytic action of CO2 is mainly due to intermediacy of the CO3 radical anion and NO2 radical into which ONOOCO2- is fragmented while reacting with guanine. The relative stabilities of the different product complexes suggest that the mutation caused by ONOO- in the presence of CO2 would mainly involve 8-oxoG.     

Reaction mechanism of methanol decomposition on Pt‐based model catalysts: A theoretical study

The decomposition mechanisms of methanol on five different Pt surfaces, the flat surface of Pt(111), Pt‐defect, Pt‐step, Pt(110)(1 × 1), and Pt(110)(2 × 1), have been studied with the DFT‐GGA method using the repeated slab model. The adsorption energies under the most stable configuration of the possible species and the activation energy barriers of the possible elementary reactions involved are obtained in this work. Through systematic calculations for the reaction mechanism of methanol decomposition on these surfaces, we found that such a reaction shows the same reaction mechanism on these Pt‐based model catalysts, that is, the final products are all H (H_ads) and CO (CO_ads) via O? H bond breaking in methanol and C? H bond scission in methoxy. These results are in general agreement with the previous experimental observations. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010.