Theoretical study of the adsorption and dissociation of oxygen on Pt(111) in the presence of homogeneous electric fields

The effect of homogeneous electric fields on the adsorption energies of atomic and molecular oxygen and the dissociation activation energy of molecular oxygen on Pt(111) were studied by density functional theory (DFT). Positive electric fields, corresponding to positively charged surfaces, reduce the adsorption energies of the oxygen species on Pt(111), whereas negative fields increase the adsorption energies. The magnitude of the energy change for a given field is primarily determined by the static surface dipole moment induced by adsorption. On 10-atom Pt(111) clusters, the adsorption energy of atomic oxygen decreased by ca. 0.25 eV in the presence of a 0.51 V/A (0.01 au) electric field. This energy change, however, is heavily dependent on the number of atoms in the Pt(111) cluster, as the static dipole moment decreases with cluster size. Similar calculations with periodic slab models revealed a change in energy smaller by roughly an order of magnitude relative to the 10-atom cluster results. Calculations with adsorbed molecular oxygen and its transition state for dissociation showed similar behavior. Additionally, substrate relaxation in periodic slab models lowers the static dipole moment and, therefore, the effect of electric field on binding energy. The results presented in this paper indicate that the electrostatic effect of electric fields at fuel cell cathodes may be sufficiently large to influence the oxygen reduction reaction kinetics by increasing the activation energy for dissociation.     

Silver as an electron source for photodissociation of hydronium

The photochemistry of a solvated hydronium ion near a silver surface is investigated using ab initio self-consistent field and configuration interaction theory. Photoinduced electron attachment can occur at energies in the range of 1.1-1.2 eV depending upon the initial orientation of the hydronium relative to the silver surface. Rearrangement of solvating waters considerably reduces transition state barriers to dissociation on the excited-state potential energy surface, such that fast dissociation of the neutralized hydronium would occur with no barrier. Both the H and H(2) product channels are exothermic pathways on the excited state surface and in several instances exothermic compared to the energy of the initial structure.     

CO2在Cu表面还原成碳氢化合物的DFT计算研究

应用密度泛函理论（DFT）反应能计算及最小能量路径分析研究了CO2在气相和电化学环境中于Cu(111)单晶表面的还原过程。气相中,CO2还原为碳氢化合物的反应路径可能为：CO2(g) + H* → COOH* → (CO +OH)* → CHO*;CHO + H* → CH2O* → (CH2 + O)*;CH2* + 2H* → CH4或2CH2* → C2H4。整个反应由CO2(g) + H* → COOH* → (CO +OH)*,(CO + H)* → CHO*和CH2O* → (CH2 + O)*等几个步骤联合控制。在-0.50V (vs RHE) 以正的电势下,CO2在Cu(111)表面电化学还原主要形成HCOO-和CO吸附物;随着电势逐渐负移,CO2加氢解离形成CO的反应越来越容易,CO成为主要产物;随电势进一步变负,形成碳氢化合物的趋势逐渐变强。与CO2的气相化学还原不同的是,电化学环境下CO质子化形成的CHO中间体倾向于解离形成CH,而在气相中CHO中间体则倾向于进一步质子化形成CH2O中间体。     

Coverage dependence and hydroperoxyl-mediated pathway of catalytic water formation on Pt (111) surface

Hydrogen oxidation on Pt (111) surface is modeled by density functional theory (DFT). Previous DFT calculations showed too large O2 dissociation barriers, but we find them highly coverage dependent: when the coverage is low, dissociation barriers close to experimental values (approximately 0.3 eV) are obtained. For the whole reaction, a new pathway involving hydroperoxyl (OOH) intermediate is found, with the highest reaction barrier of only approximately 0.4 eV. This may explain the experimental observation of catalytic water formation on Pt (111) surface above the H2O desorption temperature of 170 K, despite that the direct reaction between chemisorbed O and H atoms is a highly activated process with barrier approximately 1 eV as previous calculations showed.     

Mechanism and Tafel lines of electro-oxidation of water to oxygen on RuO2(110)

How to efficiently oxidize H(2)O to O(2) (H(2)O → 1/2O(2) + 2H(+) + 2e(-)) is a great challenge for electrochemical/photo water splitting owing to the high overpotential and catalyst corrosion. Here extensive periodic first-principles calculations integrated with modified-Poisson-Boltzmann electrostatics are utilized to reveal the physical origin of the high overpotential of the electrocatalytic oxygen evolution reaction (OER) on RuO(2)(110). By determining the surface phase diagram, exploring the possible reaction channels, and computing the Tafel lines, we are able to elucidate some long-standing puzzles on the OER kinetics from the atomic level. We show that OER occurs directly on an O-terminated surface phase above 1.58 V vs NHE, but indirectly on a OH/O mixed phase below 1.58 V by converting first the OH/O mixed phase to the O-terminated phase locally. The rate-determining step of OER involves an unusual water oxidation reaction following a Eley-Rideal-like mechanism, where a water molecule from solution breaks its OH bond over surface Os with concurrent new O-OH bond formation. The free energy barrier is 0.74 eV at 1.58 V, and it decreases linearly with the increase of potential above 1.58 V (a slope of 0.56). In contrast, the traditionally regarded surface oxygen coupling reaction with a Langmuir-Hinshelwood mechanism is energetically less favored and its barrier is weakly affected by the potential. Fundamentally, we show that the empirical linear barrier~potential relation is caused by the linear structural response of the solvated transition state to the change of potential. Finally, the general strategy for finding better OER anode is also presented.     

Density functional theory model study of size and structure effects on water dissociation by platinum nanoparticles

Size and structure effects on the homolytic water dissociation reaction mediated by Pt nanoparticles have been investigated through density functional theory calculations carried out on a series of cubooctahedral Pt(n) nanoparticles of increasing sizes (n = 13, 19, 38, 55, 79, and 140). Water adsorption energy is not significantly influenced by the nanoparticle size. However, activation energy barrier strongly depends on the particle size. In general, the activation energy barrier increases with nanoparticles size, varying from 0.30 eV for Pt(19) to 0.70 eV for Pt(140). For the largest particle the calculated barrier is very close to that predicted for water dissociation on Pt(111) (0.78 eV) even though the reaction mediated by the Pt nanoparticles involves adsorption sites not present on the extended surface.     

A first principles study of oxygen reduction reaction on a Pt(111) surface modified by a subsurface transition metal M (M = Ni, Co, or Fe)

We have performed first-principle density functional theory calculations to investigate how a subsurface transition metal M (M = Ni, Co, or Fe) affects the energetics and mechanisms of oxygen reduction reaction (ORR) on the outermost Pt mono-surface layer of Pt/M(111) surfaces. In this work, we found that the subsurface Ni, Co, and Fe could down-shift the d-band center of the Pt surface layer and thus weaken the binding of chemical species to the Pt/M(111) surface. Moreover, the subsurface Ni, Co, and Fe could modify the heat of reaction and activation energy of various elementary reactions of ORR on these Pt/M(111) surfaces. Our DFT results revealed that, due to the influence of the subsurface Ni, Co, and Fe, ORR would adopt a hydrogen peroxide dissociation mechanism with an activation energy of 0.15 eV on Pt/Ni(111), 0.17 eV on Pt/Co(111), and 0.16 eV on Pt/Fe(111) surface, respectively, for their rate-determining O2 protonation reaction. In contrast, ORR would follow a peroxyl dissociation mechanism on a pure Pt(111) surface with an activation energy of 0.79 eV for its rate-determining O protonation reaction. Thus, our theoretical study explained why the subsurface Ni, Co, and Fe could lead to multi-fold enhancement in catalytic activity for ORR on the Pt mono-surface layer of Pt/M(111) surfaces.     

Mechanistic study of the electrochemical oxygen reduction reaction on Pt(111) using density functional theory

Density functional theory (DFT) was used to study the electrolyte solution effects on the oxygen reduction reaction (ORR) on Pt(111). To model the acid electrolyte, an H(5)O(2)(+) cluster was used. The vibrational proton oscillation modes for adsorbed H(5)O(2)(+) computed at 1711 and 1010 cm(-1), in addition to OH stretching and H(2)O scissoring modes, agree with experimental vibrational spectra for proton formation on Pt surfaces in ultrahigh vacuum. Using the H(5)O(2)(+) model, protonation of adsorbed species was found to be facile and consistent with the activation barrier of proton transfer in solution. After protonation, OOH dissociates with an activation barrier of 0.22 eV, similar to the barrier for O(2) dissociation. Comparison of the two pathways suggests that O(2) protonation precedes dissociation in the oxygen reduction reaction. Additionally, an OH diffusion step following O protonation inhibits the reaction, which may lead to accumulation of oxygen on the electrode surface.     

Cu_2O(111)催化水煤气变换反应机理的理论研究

采用密度泛函理论结合周期平板模型方法系统地研究了水煤气变换反应在Cu_2O(111)表面上的反应机理,包括氧化还原机理、羧基机理和甲酸根机理.结果表明,在Cu_2O(111)表面,羧基机理和甲酸根机理均可行,且甲酸根机理更为有利,其最佳反应途径为H_2O~*→H~*+OH~*;CO(g)+H~*+OH~*→trans-HCOOH~*(1)→cis-HCOOH~*→CO_2~*+H_2(g).其中trans-HCOOH~*(1)→cis-HCOOH~*为其决速步,该基元反应的能垒仅为59 kJ·mol~(-1).羧基机理的最优反应路径同样是以H_2O的解离反应开始,随后CO(g)+OH~*→cis-COOH~*→trans-COOH~*→CO_2(g)+H~*,最后产生的两个吸附的H原子先迁移再结合生成H_2,整个反应的控速步骤为H原子的迁移,迁移能垒为96 kJ·mol~(-1).氧化还原机理则由于OH解离需要越过一个很高的能垒(254 vs.187 kJ·mol~(-1))而不可行.     

The reaction pathways for HSCH3 adsorption on Au(111): a density functional theory study

Density functional theory was used to investigate the reaction pathways for HSCH(3) adsorption on Au(111) at low coverage. A molecular adsorbed state was found with the S atom bond on Top sites (E approximately -0.38 eV) and molecular adsorption is nonactivated. The H-SCH(3) dissociation process is energetically less favorable and becomes slightly exothermic only when surface relaxation is considered (DeltaE approximately -0.2 eV). All the reaction pathways present a sizable activation energy barrier, with the lowest being approximately 0.52 eV (0.41 eV taking into account slab relaxation). In the corresponding saddle point of the potential energy surface, the S atom of the methylthiolate molecule is placed on Top sites and the H near a Bridge site. The high barrier obtained explains the complete absence of reactive methanethiol dissociation found in recent experiments.     

The mechanism of N2O formation via the (NO)2 dimer: a density functional theory study

Catalytic formation of N(2)O via a (NO)(2) intermediate was studied employing density functional theory with generalized gradient approximations. Dimer formation was not favored on Pt(111), in agreement with previous reports. On Pt(211) a variety of dimer structures were studied, including trans-(NO)(2) and cis-(NO)(2) configurations. A possible pathway involving (NO)(2) formation at the terrace near to a Pt step is identified as the possible mechanism for low-temperature N(2)O formation. The dimer is stabilized by bond formation between one O atom of the dimer and two Pt step atoms. The overall mechanism has a low barrier of approximately 0.32 eV. The mechanism is also put into the context of the overall NO + H(2) reaction. A consideration of the step-wise hydrogenation of O(ads) from the step is also presented. Removal of O(ads) from the step is significantly different from O(ads) hydrogenation on Pt(111). The energetically favored structure at the transition state for OH(ads) formation has an activation energy of 0.63 eV. Further hydrogenation of OH(ads) has an activation energy of 0.80 eV.     

Adsorption and dissociation of H2O on Cu2O(100):A computational study

The adsorption and dissociation of water on Cu2O(100) have been investigated by the density functional theory-generalized gradient approximation (DFT-GGA) method. The corresponding reaction energies, the structures of the transition states and the activation energies were determined. Calculations with and without dipole correction were both studied to get an understanding of the effect of the dipole moment on the adsorption and reaction of water on dipole surface Cu2O(100). When dipole correction was added, the adsorption energies of H2O on different sites generally decreased. The calculated activation barriers for HxO (x = 1, 2) dehydrogenation are 0.42 eV (1.01 eV without the dipole correction) and 1.86 eV, respectively, including the zero point energy correction. The first dehydrogenation outcome is energetically the most stable product.     

Catalytic oxidation activity of Pt3O4 surfaces and thin films

The catalytic oxidation activity of platinum particles in automobile catalysts is thought to originate from the presence of highly reactive superficial oxide phases which form under oxygen-rich reaction conditions. Here we study the thermodynamic stability of platinum oxide surfaces and thin films and their reactivities toward oxidation of carbon compounds by means of first-principles atomistic thermodynamics calculations and molecular dynamics simulations based on density functional theory. On the Pt(111) surface the most stable superficial oxide phase is found to be a thin layer of alpha-PtO2, which appears not to be reactive toward either methane dissociation or carbon monoxide oxidation. A PtO-like structure is most stable on the Pt(100) surface at oxygen coverages of one monolayer, while the formation of a coherent and stress-free Pt3O4 film is favored at higher coverages. Bulk Pt3O4 is found to be thermodynamically stable in a region around 900 K at atmospheric pressure. The computed net driving force for the dissociation of methane on the Pt3O4(100) surface is much larger than that on all other metallic and oxide surfaces investigated. Moreover, the enthalpy barrier for the adsorption of CO molecules on oxygen atoms of this surface is as low as 0.34 eV, and desorption of CO2 is observed to occur without any appreciable energy barrier in molecular dynamics simulations. These results, combined, indicate a high catalytic oxidation activity of Pt3O4 phases that can be relevant in the contexts of Pt-based automobile catalysts and gas sensors.     

Pd掺杂对ZnO(1120)面上水解离的影响

采用密度泛函理论计算研究了清洁的以及Pd掺杂的ZnO(1120)面上水分子的吸附和解离.结果表明,在清洁ZnO(1120)上,水分子倾向于分子吸附,解离吸附较为困难.在Pd掺杂的ZnO上,水分子仍倾向吸附在Zn原子上,且吸附能与其在清洁ZnO表面的相当.然而,Pd的掺杂可增强水解离产物OH和H的吸附,从而显著提高了水的解离活性,相应的水解离能垒为0.36eV,放热0.21eV.     

Potential energy surface for activation of methane by Pt(+): a combined guided ion beam and DFT study

A guided-ion beam tandem mass spectrometer is used to study the reactions of Pt(+) with methane, PtCH(2)(+) with H(2) and D(2), and collision-induced dissociation of PtCH(4)(+) and PtCH(2)(+) with Xe. These studies experimentally probe the potential energy surface for the activation of methane by Pt(+). For the reaction of Pt(+) with methane, dehydrogenation to form PtCH(2)(+) + H(2) is exothermic, efficient, and the only process observed at low energies. PtH(+), formed in a simple C-H bond cleavage, dominates the product spectrum at high energies. The observation of a PtH(2)(+) product provides evidence that methane activation proceeds via a (H)(2)PtCH(2)(+) intermediate. Modeling of the endothermic reaction cross sections yields the 0 K bond dissociation energies in eV (kJ/mol) of D(0)(Pt(+)-H) = 2.81 +/- 0.05 (271 +/- 5), D(0)(Pt(+)-2H) = 6.00 +/- 0.12 (579 +/- 12), D(0)(Pt(+)-C) = 5.43 +/- 0.05 (524 +/- 5), D(0)(Pt(+)-CH) = 5.56 +/- 0.10 (536 +/- 10), and D(0)(Pt(+)-CH(3)) = 2.67 +/- 0.08 (258 +/- 8). D(0)(Pt(+)-CH(2)) = 4.80 +/- 0.03 eV (463 +/- 3 kJ/mol) is determined by measuring the forward and reverse reaction rates for Pt(+) + CH(4) right harpoon over left harpoon PtCH(2)(+) + H(2) at thermal energy. We find extensive hydrogen scrambling in the reaction of PtCH(2)(+) with D(2). Collision-induced dissociation (CID) of PtCH(4)(+), identified as the H-Pt(+)-CH(3) intermediate, with Xe reveals a bond energy of 1.77 +/- 0.08 eV (171 +/- 8 kJ/mol) relative to Pt(+) + CH(4). The experimental thermochemistry is favorably compared with density functional theory calculations (B3LYP using several basis sets), which also establish the electronic structures of these species and provide insight into the reaction mechanism. Results for the reaction of Pt(+) with methane are compared with those for the analogous palladium system and the differences in reactivity and mechanism are discussed.     

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.     

Ion-molecule rate constants and branching ratios for the reaction of N(3) (+)+O(2) from 120 to 1400 K

The kinetics of the reaction of N(3) (+) with O(2) has been studied from 120 to 1400 K using both a selected ion flow tube and high-temperature flowing afterglow. The rate constant decreases from 120 K to approximately 1200 K and then increases slightly up to the maximum temperature studied, 1400 K. The rate constant compares well to most of the previous measurements in the overlapping temperature range. Comparing the results to drift tube data shows that there is not a large difference between increasing the translational energy available for reaction and increasing the internal energy of the reactants over much of the range, i.e., all types of energies drive the reactivity equally. The reaction produces both NO(+) and NO(2) (+), the latter of which is shown to be the higher energy NOO(+) linear isomer. The ratio of NOO(+) to NO(+) decreases from a value of over 2 at 120 K to less than 0.01 at 1400 K because of dissociation of NOO(+) at the higher temperatures. This ratio decreases exponentially with increasing temperature. High-level theoretical calculations have also been performed to compliment the data. Calculations using multi-reference configuration interaction theory at the MRCISD(Q)/cc-pVTZ level of theory show that singlet NOO(+) is linear and is 4.5 eV higher in energy than ONO(+). A barrier of 0.9 eV prevents dissociation into NO(+) and O((1)D); however, a crossing to a triplet surface connects to NO(+) and O((3)P) products. A singlet and a triplet potential energy surface leading to products have been determined using coupled cluster theory at the CCSD(T)/aug-cc-pVQZ level on structures optimized at the Becke3-Lee, Yang, and Parr (B3LYP)/aug-cc-pVTZ level of theory. The experimental results and reaction mechanism are evaluated using these surfaces.     

Energy of charged states in the RDX crystal: trapping of charge-transfer pairs as a possible mechanism for initiating detonation

Calculations for the crystalline energetic material RDX (1,3,5-trinitro-1,3,5-triazacyclohexane) yield the effective polarizability (17.2 angstroms3), local electric field tensor, effective dipole moment (9.40 D), and dipole-dipole energy (-27.2 kJ/mol). Fourier-transform techniques give the polarization energy P for a single charge in the perfect crystal as -1.14 eV; the charge-dipole energy W(D) is zero if the crystal carries no bulk dipole moment. Polarization energies for charge-transfer (CT) pairs combine with the Coulomb energy E(C) to give the screened Coulomb energy E(scr); screening is nearly isotropic with E(scr) approximately = E(C)2.6. For CT pairs W(D) reduces to a term deltaW(D) arising from the interaction of the charge on each ion with the change in dipole moment on the other ion relative to the neutral molecule. The dipole moments are calculated as 7.40 D for the neutral molecule and 6.84 D and 7.44 D for the anion and cation, giving the lowest two CT pairs at -1.34 eV and -0.94 eV. The changes in P and W(D) near a molecular vacancy yield traps with depths that reach 400 meV for single charges and 185 meV for the nearest-neighbor CT pair. Divacancies yield traps with depths nearly equal to the sum of those produced by the separate vacancies. These results are consistent with a mechanism in which detonation of RDX is initiated by mechanical generation of CT pairs that localize at vacancies, recombine, and release energy sufficient to break bonds; crystals of molecules with lower dipole moments should be less sensitive.     

Ionization induced relaxation in solvation structure: a comparison between Na(H2O)n and Na(NH3)n

The constant ionization potential for hydrated sodium clusters Na(H2O)n just beyond n=4, as observed in photoionization experiments, has long been a puzzle in violation of the well-known (n+1)(-1/3) rule that governs the gradual transition in properties from clusters to the bulk. Based on first principles calculations, a link is identified between this puzzle and an important process in solution: the reorganization of the solvation structure after the removal of a charged particle. Na(H2O)n is a prototypical system with a solvated electron coexisting with a solvated sodium ion, and the cluster structure is determined by a balance among three factors: solute-solvent (Na+-H2O), solvent-solvent (H2O-H2O), and electron-solvent (OH{e}HO) interactions. Upon the removal of an electron by photoionization, extensive structural reorganization is induced to reorient OH{e}HO features in the neutral Na(H2O)n for better Na+-H2O and H2O-H2O interactions in the cationic Na+(H2O)n. The large amount of energy released, often reaching 1 eV or more, indicates that experimentally measured ion signals actually come from autoionization via vertical excitation to high Rydberg states below the vertical ionization potential, which induces extensive structural reorganization and the loss of a few solvent molecules. It provides a coherent explanation for all the peculiar features in the ionization experiments, not only for Na(H2O)n but also for Li(H2O)n and Cs(H2O)n. In addition, the contrast between Na(H2O)n and Na(NH3)n experiments is accounted for by the much smaller relaxation energy for Na(NH3)n, for which the structures and energetics are also elucidated.     

Density-functional study for the NOx (x = 1, 2) dissociation mechanism on the Cu(1 1 1) surface

Spin-polarized density functional theory calculation is employed to study the adsorption and dissociation of NO₂ molecule on Cu(1 1 1) surface. It is shown that the most favorable adsorption structure is the NO₂ (T,T-O-,O′-nitrito) configuration which has an adsorption energy of −1.49 eV. The barriers for step-wise NO₂ dissociation reaction, NO_2(g) → N_(a) + 2O_(a), are 1.05 (for O–N–O bond activation), and 2.08 eV (for N–O bond activation), respectively, and the entire process is 0.6 eV exothermic. The energetics of single N–O dissociation with and without the presence of N atom or O atom on the surface are also calculated. The results indicate that in the presence of O atom on Cu(1 1 1) surface would raise the N–O dissociation barrier, whereas in the presence of N atom decrease it. The interaction nature between adsorbates and substrate is analyzed by the local density of states (LDOS) calculation.