Propene Oxidation with the Anionic Cluster V4o11-: Selective Epoxidation
Hai-Bei Li,Shan Xi Tian,Jinlong Yang
DOI: https://doi.org/10.1002/chem.200901165
2009-01-01
Chemistry - A European Journal
Abstract:Vanadium cage catalysis: High selectivity of the epoxidation of propene with anionic cluster V4O11− is demonstrated on the basis of the first-principle calculations. Oxygen, proton, and negative charge transfers are the general reaction mechanisms (see figure) of the hydrocarbon oxidations by vanadium oxides. Vanadium oxides have been extensively studied, owing to the diverse applications as bolometric detectors,1 cathode materials in lithium batteries,2 spintronic devices,3 and the supported catalysts.4 Toward an understanding of the operation mode of the catalysts, the reaction mechanism of the nano-sized and subnano-sized vanadium oxide (VxOy) is needed to be established at the molecular level. Recently, a lot of efforts both of experiments and theoretical calculations have been made to reveal the chemical reactivity of the neutral and ionic VxOy+/0/− clusters.5–7 Oxidation of the hydrocarbons with the VxOy solid surface4 or clusters5–7 is such reaction in which the oxygen atom absorbed on VxOy is transferred to hydrocarbons then the oxide products are produced. Propylene oxide (PO), as an important hydrocarbon oxide, is widely used in manufacturing polyurethane and glycols. The oxidation of propene coabsorbed on a active site of VxOy in heterogeneous catalysis is promising route, but it is currently a challenging goal of the production of PO with both the high selectivity and conversion.5–7 Molecular association and oxygen transfer are frequently observed in the reactions between the cationic VxOy+ cluster and alkyl or allyl hydrocarbons, which is extremely sensitive to the stoichiometry of the cationic clusters;5, 6 Propene can be oxidated by the neutral oxygen-rich cluster VO3(V2O5)n through the CC bond cleavage.6b Recently, Li and co-workers found the unique chemical reactivity between V4O11− cluster and propene in the photocatalysis experiments.7 The predominated products are V4O10− and C3H6O, while V4O10H2− and C3H4O are affiliated. The former oxide products C3H6O, possibly, PO, propanal (PR), and acetone (AC), are not further identified.7 With respect to PO, PR and AC can be regarded as the by-products suffering the PO outcome in the industry.4 Herein the association complex [V4O11⋅C3H6]−, as an excellent cluster model of the coadsorption of propene on a VxOy active site in heterogeneous catalysis, is studied by the first-principle calculations. A high selectivity to produce PO via the efficient oxygen transfer is demonstrated. Two low-lying Cs symmetrical V4O11− clusters, V4O11−(1) and V4O11−(2), as shown in Figure 1, are selected because they have been proved to be the most stable ones (nearly degenerate in energy) in gas phase8 and their low-symmetrical structures having both the bridged peroxo and the open dioxo units mimic well the reactants adsorbed and activated on the solid surface. After the oxygen transfer to propene, two possible clusters V4O10−(1) (D2d symmetry) and V4O10−(2) (Cs symmetry) are formed, and the corresponding oxide products C3H6O may be PO, PR, and AC. The highly symmetrical V4O10−(1) can be served to approximate vanadium oxide surface structure, while V4O10− (2) can represent a model of some suggested defect sites of the solid surface.9 The structures of the reactants, products, and the association complexes [V4O11⋅C3H6]− are fully optimized by the first-principle hybrid density functional B3LYP method10 together with the TZVP basis set.11 According to the frontier molecular orbitals of V4O11−(1) and V4O11−(2),8 these two structures have the different active sites, implying that there may be the different reactions pathways after the adsorption of the propene. Two most stable structures of V4O11− and the transformation from V4O11− to V4O10− after the atomic oxygen transfer. The relative energies (in kcal mol−1) are given in the parenthesis. a–d, a′–d′ denote the potential active sites. Our trial calculations indicate that the dioxo units (Oa−d, Oa′−d′) are much more active than the bridged peroxo units for the associations with propene. For the reactions between V4O11−(1) and propene, as depicted in Figure 2, the initial step of five parallel reaction pathways is the formation of a very stable association complex. In the paths 1B–1E, one covalent bond is formed between oxygen and carbon (any one in CC) atoms; while two bonds formed in the path 1 A which leads to the most stable association complex, 1 A-R. A comparison among the five parallel reaction energy barriers reveals that the path 1 E (the lowest barrier 3.48 kcal mol−1) is the most favourable for the oxygen transfer process, and the products of this path spanning over the transition state 1 E-TS are V4O10−(1) and PO. The energy barriers for the paths 1 A–1D are 42.87, 20.97, 11.54, and 8.48 kcal mol−1, respectively. PR is produced in the paths 1 A, while AC is the outcome in the paths 1 B and 1 D. Only two oxides, PO and AC, are produced in the reactions between V4O11−(2) and propene, as shown in Figure 3. Two reaction paths, 2 D and 2 C, exhibit the comparable energy barriers 12.64 and 12.85 kcal mol−1 that are much lower than those of paths 2 A-1, 2 A-2, and 2 B (the barriers are 25.19 kcal mol−1, 17.77 kcal mol−1, and 18.05 kcal mol−1, respectively). In the reactions between V4O11−(2) and propene, AC is produced in the path 2 A-1, while the others yield PO. Structures and relative energetic profiles of the five parallel reactive pathways for the reaction between V4O11−(1) and propene. All energies (in kcal mol−1) are corrected with zero-point vibrational energy. Structures and relative energetic profiles of the five parallel reactive pathways for the reaction between V4O11−(2) and propene. All energies (in kcal mol−1) are corrected with zero-point vibrational energy. In order to assess the relative importance of the different reaction channels, we performed the Rice-Ramsperger–Kassel–Marcus (RRKM) calculations to predict the rate constants (k), over a grid of energies up to 200 kcal mol−1 of internal energy. No tunnelling correction was included because: Firstly, in the most cases the rates are extremely high when the internal energy is above the reaction threshold; Secondly, the most barrier heights are less than 20 kcal mol−1 and do not warrant the attempt to achieve the accuracy addressed by a tunnelling correction. The rate constants k as a function of the internal energies are plotted for the reactions of V4O11−(1) in Figure 4a and for V4O11−(2) in Figure 4b. The above direct oxidation (or oxygen transfer) processes are treated as the first-order reactions. Evidently, the paths 1 E and 2 D with the lowest-energy barriers show the highest-rate constants over the whole energy range. Both of these two paths yield the common product PO, indicating that the epoxidation of the propene is the overwhelming process among all reaction paths. In the experiment,6 ultraviolet femtosecond laser (λcenter=272 nm) was used to assist the oxygen transfer from V4O11− to propene. This energy is high enough to produce the hydrocarbon oxides, PO, PR, and AC. The energy barrier in the path 1E 3.48 kcal mol−1 is as low as the energy of the fundamental dioxo VO stretching vibration (≈1000 cm−1) of the isolated V4O11−,8 implying that an infrared absorption may touch off the reaction. Rate constants for the five parallel reaction pathways of the reaction between V4O11− [a) 1; b) 2] and propene. The physicochemical nature of the parallel reactions between V4O11−(1, 2) and propene are further explored by inspection of the charge-density difference maps (as shown in the supporting Figure S1 and S2, in which the site-changes of the negative charge are plotted). These charge transfers are succinctly illustrated in Scheme 1 to grasp the distinct differences in each reaction paths. In the isolated cluster V4O11−(1), the negative charge is localized at one vanadyl dioxo group (the most on Ob atom). Along the interaction pathway 1 E, the ethylene-end carbon atom C1 (in the propene C1C2C3) initially attacks Oa atom, then propene aggregates to form the association complex (1 E-R). This is similar to that in the reaction between VxOy+ and the ethylene.6 In the transition state structure 1 E-TS, the delocalized negative charge is partially at both the second carbon C2 and V atoms. Finally, the negative charge is completely localized at the V atom in the V4O10− moiety of 1 E-P (as shown in the left side of Scheme 1a, then it is delocalized at the isolated V4O10−(1). As shown in Figure S2, the charge transfer of the path 2 B is in a similar scenario. When the middle carbon atom C2 in propene attacks the O site, e.g., along the paths 1 D, 2 A-2, and 2 C, the negative charge is transferred to the methylene end (as shown in the middle and right panels of Scheme 1a. AC is produced in the path 1D where a proton connecting the C2 atom shifts to the methylene. This concerted oxygen/proton transfer process is illustrated in the right panel of Scheme 1a. Only in the path 1 A the negative charge is localized at the V atom all long, see Scheme 1b and Figure S1, where the proton of methylene group is transferred to the middle site and PR is finally produced. Mechanism of the reaction between V4O11− and propene. PR and AC were predicted to be produced in the reactions between VO2+ and propene in which propene also combined with the active oxygen site.6c The oxygen/proton transfer mechanism was proposed for the above reactions.5, 6 Dibenedetto et al. found AC was the absolutely predominated product in the V2O5-propene reaction.5a Up to date, there are no reports that PO can be the main product in the reactions with VxOy0/+. A previous study indicated that neither the association clusters nor oxygen transfer products were found for the anionic VxOy− and C2-hydrocarbons (e.g., ethane, ethylene), hypothesizing the reasons in the negative charge and oxygen site of the anionic VxOy−, or the decreased polarity of C2-hydrocarbons.6h On the basis of the experimental findings by Li et al.,7 our comprehensive theoretical study suggests that the methyl group in propene should be in charge of the high reactivity with V4O11−. The methyl group has strong electron withdrawing effect, which weakens the double bond CC and makes propene more active than ethylene. The more polarized propene can depress the activation barriers for reaction. Both of the present study and the previous work on the VxOy+ reactions6a indicate that the electron transfer induces the oxygen transfer from vanadium oxide to hydrocarbon. This mechanism is obviously practicable and efficient for the negatively charged V4O11−. Owing to the high electron affinity (≈4.16 eV) of V4O10−(1),12 the negative charge transfers back to the vanadium oxide moiety after spanning the transition states [V4O11⋅C3H6]−≠. Thus, in the thermodynamical scheme, a direct gas-phase oxidation of propene by oxygen atom is readily achieved in the reactions with V4O11−. The interaction between propene and oxygen atom yields PO without any energetic barrier. As shown in the supporting Figure S3, we also explore the reaction paths for O−-propene system. PR and AC, not PO, are the more favourable products. This implies that the partial even unit negative charge on the oxygen of the dioxo group can not accelerate the reactions to PO. On the other hand, the transferred partial negative charge on the C1 or C2 atom at the transition states of the association complexes (see Scheme 1) may boost the epoxidation or proton transfer. In summary, the reactions of the vanadium oxide anion clusters V4O11− with propene have been clarified by the first-principle calculations. The main C3H6O products observed in the photocatalytic experiment7 are identified as PO, PR, and AC, and further rationalized with the mechanism involving oxygen, proton, and the negative charge transfers. An excellent cluster model of the heterogeneous catalysis is established for the activation of the molecular oxygen and the coadsorption of propene on vanadium oxide anion clusters, in which the higher selectivity to produce PO is demonstrated, moreover, the much higher conversion to PO is hopefully archived. In which ρA and ρN represented the electron density distributions of the anionic and neutral species, respectively. The footnote anion denoted that the anionic structures were used through the above calculations. This work is supported by NNSFC (Grant. 20673105), MOST (Grants. 2006CB92204, 2007CB815204), and NCET program. Detailed facts of importance to specialist readers are published as ”Supporting Information”. Such documents are peer-reviewed, but not copy-edited or typeset. They are made available as submitted by the authors. Please note: The publisher is not responsible for the content or functionality of any supporting information supplied by the authors. Any queries (other than missing content) should be directed to the corresponding author for the article.