Partial oxidation of propylene catalyzed by VO3 clusters: a density functional theory study

Density functional theory (DFT) calculations are carried out to investigate partial oxidation of propylene over neutral VO 3 clusters. C=C bond cleavage products CH 3CHO + VO 2CH 2 and HCHO + VO 2CHCH 3 can be formed overall barrierlessly from the reaction of propylene with VO 3 at room temperature. Formation of hydrogen transfer products H 2O + VO 2C 3H 4, CH 2=CHCHO + VO 2H 2, CH 3CH 2CHO + VO 2, and (CH 3) 2CO + VO 2 is subject to tiny (0.01 eV) or small (0.06 eV, 0.19 eV) overall free energy barriers, although their formation is thermodynamically more favorable than the formation of C=C bond cleavage products. These DFT results are in agreement with recent experimental observations. VO 3 regeneration processes at room temperature are also investigated through reaction of O 2 with the CC bond cleavage products VO 2CH 2 and VO 2CHCH 3. The following barrierless reaction channels are identified: VO 2CH 2 + O 2 --> VO 3 + CH 2O; VO 2CH 2 + O 2 --> VO 3C + H 2O, VO 3C + O 2 --> VO 3 + CO 2; VO 2CHCH 3 + O 2 --> VO 3 + CH 3CHO; and VO 2CHCH 3 + O 2 --> VO 3C + CH 3OH, VO 3C + O 2 --> VO 3 + CO 2. The kinetically most favorable reaction products are CH 3CHO, H 2O, and CO 2 in the gas phase model catalytic cycles. The results parallel similar behavior in the selective oxidation of propylene over condensed phase V 2O 5/SiO 2 catalysts.     

Guided-ion beam and theoretical study of the potential energy surface for activation of methane by W+

A guided-ion beam tandem mass spectrometer is used to study the reactions, W(+) + CH(4) (CD(4)) and [W,C,2H](+) + H(2) (D(2)), to probe the [W,C,4H](+) potential energy surface. The reaction W(+) + CH(4) produces [W,C,2H](+) in the only low-energy process. The analogous reaction in the CD(4) system exhibits a cross section with strong differences at the lowest energies caused by zero-point energy differences, demonstrating that this reaction is slightly exothermic for CH(4) and slightly endothermic for CD(4). The [W,C,2H](+) product ion reacts further at thermal energies with CH(4) to produce W(CH(2))(x)(+) (x = 2-4). At higher energies, the W(+) + CH(4) reaction forms WH(+) as the dominant ionic product with smaller amounts of WCH(3)(+), WCH(+), and WC(+) also formed. The energy dependent cross sections for endothermic formation of the various products are analyzed and allow the determination of D(0)(W(+)-CH(3)) approximately 2.31 +/- 0.10 eV, D(0)(W(+)-CH(2)) = 4.74 +/- 0.03 eV, D(0)(W(+)-CH) = 6.01 +/- 0.28 eV, and D(0)(W(+)-C) = 4.96 +/- 0.22 eV. We also examine the reverse reaction, [W,C,2H](+) + H(2) (D(2)) --> W(+) + CH(4) (CH(2)D(2)). Combining the cross sections for the forward and reverse processes yields an equilibrium constant from which D(0)(W(+)-CH(2)) = 4.72 +/- 0.04 eV is derived. Theoretical calculations performed at the B3LYP/HW+/6-311++G(3df,3p) level yield thermochemistry in reasonable agreement with experiment. These calculations help identify the structures and electronic states of the species involved and characterize the potential energy surface for the [W,C,4H](+) system.     

Reactions of hydrazoic acid and trimethylindium on TiO2 rutile (110) surface: a computational study on the formation of the first monolayer InN

This article reports the result of a computational study on the reaction of hydrazoic acid and trimethylindium (TMIn), coadsorbed on TiO2 rutile (110) surface. The adsorption geometries and energies of possible adsorbates including HN3-In(CH3)3(a) and its derivatives, HN3-In(CH3)2(a), N3-In(CH3)2(a), N3-In(CH3)(a), and N-In(a), have been predicted by first-principles calculations based on the density functional theory (DFT) and the pseudopotential method. The mechanisms of these surface reactions have also been explicitly elucidated with the computed potential energy surfaces. Starting from the interaction of three stable HN3 adsorbates, HN3-Ob(a), H(N2)N-Ob(a), and Ti-NN(H)N-Ob(a), where Ob is the bridged O site on the surface, with two stable intermediates from the adsorption and dissociative adsorption of TMIn, (H3C)3In-Ob(a) and (H3C)2In-Ob(a)+H3C-Ob(a), InN products can be formed exothermically via four reaction paths following the initial barrierless In-atom association with the N atom directly bonded to H, by CH4 elimination (with approximately 40 kcal/mol barriers), the InN-N bond breaking and the final CH3 elimination or migration (with <20 kcal/mol barriers). These Langmuir-Hinshelwood processes producing the two most stable InN(a) side-on adsorptions confirm that HN3 and TMIn are indeed very efficient precursors for the deposition of InN films on TiO2 nanoparticles. The result of similar calculations for the reactions occurring by the Rideal-Eley mechanism involving HN3(a)+TMIn(g) and HN3(g)+TMIn(a) indicates that they are energetically less favored and produce the less stable InN(a) with end-on configurations.     

Hafnium oxide and zirconium oxide atomic layer deposition: initial precursor and potential side-reaction product pathways with H/Si(100)-2 x 1

Hybrid density functional calculations have been carried out using cluster models of the H/Si(100)-2 x 1 surface to investigate the mechanistic details of the initial surface reactions occurring in the atomic layer deposition of hafnium and zirconium oxides (HfO2 and ZrO2). Reaction pathways involving the metal precursors ZrCl4, Zr(CH3)4, HfCl4, and Hf(CH3)4 have been examined. Pathways leading to the formation of a Zr-Si or Hf-Si linkage show a significant sensitivity to the identity of the leaving group, with chloride loss reactions being both kinetically and thermodynamically less favorable than reactions leading to the loss of a methyl group. The energetics of the Zr(CH3)4 and Hf(CH3)4 reactions are similar with an overall exothermicity of 0.3-0.4 eV and a classical barrier height of 1.1-1.2 eV. For the reaction between H2O and the H/Si(100)-2 x 1 surface, the activation energy and overall reaction enthalpy are 1.6 and -0.8 eV, respectively. Due to contamination, trace amounts of H2O may be encountered by metal precursors, leading to the formation of minor species that can lead to unanticipated side-reaction pathways. Such gas-phase reactions between the halogenated and alkylated metal precursors and H2O are exothermic with small or no reaction barriers, allowing for the possibility of metal precursor hydroxylation before the H/Si surface is encountered. Of the contaminant surface reaction pathways, the most kinetically favorable corresponds to the surface -OH deposition. Interestingly, for the hydroxylated metal precursors, a unique reaction pathway resulting in the direct formation of Si-O-Zr and Si-O-Hf linkages has been identified and found to be the most thermodynamically stable pathway available, being exothermic by approximately 1.0 eV.     

A theoretical investigation of the gas-phase oxidation reaction of the saturated tert-butyl radical.

The radical-radical reaction mechanisms and dynamics of ground-state atomic oxygen [O(3P)] with the saturated tert-butyl radical (t-C4H9) are investigated using the density functional method and the complete basis set model. Two distinctive reaction pathways are predicted to be in competition: addition and abstraction. The barrierless addition of O(3P) to t-C4H9 leads to the formation of an energy-rich intermediate (OC4H9) on the lowest doublet potential energy surface, which undergoes subsequent direct elimination or isomerization-elimination leading to various products: C3H6O + CH3, iso-C4H8O + H, C3H7O + CH2, and iso-C4H8 + OH. The respective microscopic reaction processes examined with the aid of statistical calculations, predict that the major addition pathway is the formation of acetone (C3H6O) + CH3 through a low-barrier, single-step cleavage. For the direct, barrierless H-atom abstraction mechanism producing iso-C4H8 (isobutene) + OH, which was recently reported in gas-phase crossed-beam investigations, the reaction is described in terms of both an abstraction process (major) and a short-lived addition dynamic complex (minor).     

Theoretical studies on cycloaddition reactions between the 2-aza-1,3-butadiene cation and olefins

Density functional (B3LYP) calculations, using the 6-31G basis set, have been employed to study the title reactions. For the model reaction (H(2)C=C-NH(+)=CH(2) + H(2)C=CH(2)), a complex has been formed with 6.2 kcal/mol of stabilization energy and the transition state is 4.0 kcal/mol above this complex, but 2.1 kcal/mol below the reactants. However, the substituent effects are quite remarkable. When ethene is substituted by electron-withdrawing group CN, the reaction could also yield six-membered-ring products, but the energy barriers are all more than 7 kcal/mol, which shows that CN group unfavors the reaction. The other substituents, such as CH(3)O and CH(3) groups, have also been considered in the present work, and the results show that they are favorable for the formation of six-membered-ring adducts. The calculated results have been rationalized with frontier orbital interaction and topological analysis.     

Energetics and mechanisms of C-H bond activation by a doubly charged metal ion: guided ion beam and theoretical studies of Ta2+ + CH4

A guided ion beam tandem mass spectrometer is used to study the kinetic-energy dependence of doubly charged atomic tantalum cations (Ta(2+)) reacting with CH4 and CD4. As for the analogous singly charged system, the dehydrogenation reaction to form TaCH2(2+) + H2 is exothermic. The charge-transfer reaction to form Ta(+) + CH4(+) and the charge-separation reaction to form TaH(+) + CH3(+) are also observed at low energies in exothermic processes, as is a secondary reaction of TaCH2(2+) to form TaCH3(+) + CH3(+). At higher energies, other doubly charged products, TaC(2+) and TaCH3(2+), are observed, although no formation of TaH(2+) was observed. Modeling of the endothermic cross sections provides 0 K bond dissociation energies (in electronvolts) of D0(Ta(2+)-C) = 5.42 +/- 0.19 and D0(Ta(2+)-CH3) = 3.40 +/- 0.16. These experimental bond energies are in poor agreement with density functional calculations at the B3LYP/HW+/6-311++G(3df,3p) level of theory. However, the Ta(2+)-C bond energy is in good agreement with calculations at the QCISD(T) level of theory, and the Ta(2+)-CH3 bond energy is in good agreement with density functional calculations at the BHLYP level of theory. Theoretical calculations reveal the geometric and electronic structures of all product ions and are used to map the potential energy surface, which describes the mechanism of the reaction and key intermediates. Both experimental and theoretical results suggest that TaH(+), TaCH2(2+), and TaCH3(2+) are formed through a H-Ta(2+)-CH3 intermediate.     

Trajectory surface hopping study of the O(3P) + ethylene reaction dynamics

Spin-orbit coupling (SOC) induced intersystem crossing (ISC) has long been believed to play a crucial role in determining the product distributions in the O(3P) + C2H4 reaction. In this paper, we present the first nonadiabatic dynamics study of the title reaction at two center-of-mass collision energies: 0.56 eV, which is barely above the H-atom abstraction barrier on the triplet surface, and 3.0 eV, which is in the hyperthermal regime. The calculations were performed using a quasiclassical trajectory surface hopping (TSH) method with the potential energy surface generated on the fly at the unrestricted B3LYP/6-31G(d,p) level of theory. To simplify our calculations, nonadiabatic transitions were only considered when the singlet surface intersects the triplet surface. At the crossing points, Landau-Zener transition probabilities were computed assuming a fixed spin-orbit coupling parameter, which was taken to be 70 cm-1 in most calculations. Comparison with a recent crossed molecular beam experiment at 0.56 eV collision energy shows qualitative agreement as to the primary product branching ratios, with the CH3 + CHO and H + CH2CHO channels accounting for over 70% of total product formation. However, our direct dynamics TSH calculations overestimate ISC so that the total triplet/singlet ratio is 25:75, compared to the observed 43:57. Smaller values of SOC reduce ISC, resulting in better agreement with the experimental product relative yields; we demonstrate that these smaller SOC values are close to being consistent with estimates based on CASSCF calculations. As the collision energy increases, ISC becomes much less important and at 3.0 eV, the triplet to singlet branching ratio is 71:29. As a result, the triplet products CH2 + CH2O, H + CH2CHO and OH + C2H3 dominate over the singlet products CH3 + CHO, H2 + CH2CO, etc.     

Ni+ reactions with aminoacrylonitrile, a species of potential astrochemical relevance

The reaction of aminoacrylonitrile, a species of astrochemical interest, with Ni(+)((2)D(5/2)) was investigated by means of mass spectrometry techniques and density functional theory calculations. The dominant fragmentations in the MIKE spectrum correspond to the loss of [C2,N,H3], HCN, and NH3, the loss of H2 being very minor. The structure and bonding of the different aminoacrylonitrile-Ni(+) complexes were investigated at the B3LYP/6-311G(d,p) level of theory. The same approach was employed in our survey of the corresponding potential energy surface. This survey indicates that the [C2,N,H3] neutral product can be formed either as ketenimine (CH2CNH) or acetonitrile. The formation of the latter is significantly more exothermic but involves slightly higher activation barriers; so very likely, both isomers are produced along the reaction process. The lost of HNC is not competitive with the loss of HCN, because when the former is formed the products lie higher in energy and the corresponding mechanisms involve energy barriers above the entrance channel. The loss of NH3 is associated with the formation of a complex between cyanoacetylene, HCCCN, which is very abundant in the interstellar media, and Ni(+).     

Internal proton transfer and H2 rotations in the H5(+) cluster: a marked influence on its thermal equilibrium state

Classical and path integral Monte Carlo (CMC, PIMC) "on the fly" calculations are carried out to investigate anharmonic quantum effects on the thermal equilibrium structure of the H5(+) cluster. The idea to follow in our computations is based on using a combination of the above-mentioned nuclear classical and quantum statistical methods, and first-principles density functional (DFT) electronic structure calculations. The interaction energies are computed within the DFT framework using the B3(H) hybrid functional, specially designed for hydrogen-only systems. The global minimum of the potential is predicted to be a nonplanar configuration of C(2v) symmetry, while the next three low-lying stationary points on the surface correspond to extremely low-energy barriers for the internal proton transfer and to the rotation of the H2 molecules, around the C2 axis of H5(+), connecting the symmetric C(2v) minima in the planar and nonplanar orientations. On the basis of full-dimensional converged PIMC calculations, results on the quantum vibrational zero-point energy (ZPE) and state of H5(+) are reported at a low temperature of 10 K, and the influence of the above-mentioned topological features of the surface on its probability distributions is clearly demonstrated.     

Reactions of C2H with 1- and 2-butynes: an ab initio/RRKM study of the reaction mechanism and product branching ratios

Ab initio CCSD(T)/cc-pVTZ(CBS)//B3LYP/6-311G** calculations of the C(6)H(7) potential energy surface are combined with RRKM calculations of reaction rate constants and product branching ratios to investigate the mechanism and product distribution in the C(2)H + 1-butyne/2-butyne reactions. 2-Ethynyl-1,3-butadiene (C(6)H(6)) + H and ethynylallene (C(5)H(4)) + CH(3) are predicted to be the major products of the C(2)H + 1-butyne reaction. The reaction is initiated by barrierless ethynyl additions to the acetylenic C atoms in 1-butyne and the product branching ratios depend on collision energy and the direction of the initial C(2)H attack. The 2-ethynyl-1,3-butadiene + H products are favored by the central C(2)H addition to 1-butyne, whereas ethynylallene + CH(3) are preferred for the terminal C(2)H addition. A relatively minor product favored at higher collision energies is diacetylene + C(2)H(5). Three other acyclic C(6)H(6) isomers, including 1,3-hexadiene-5-yne, 3,4-hexadiene-1-yne, and 1,3-hexadiyne, can be formed as less important products, but the production of the cyclic C(6)H(6) species, fulvene, and dimethylenecyclobut-1-ene (DMCB), is predicted to be negligible. The qualitative disagreement with the recently measured experimental product distribution of C(6)H(6) isomers is attributed to a possible role of the secondary 2-ethynyl-1,3-butadiene + H reaction, which may generate fulvene as a significant product. Also, the photoionization energy curve assigned to DMCB in experiment may originate from vibrationally excited 2-ethynyl-1,3-butadiene molecules. For the C(2)H + 2-butyne reaction, the calculations predict the C(5)H(4) isomer methyldiacetylene + CH(3) to be the dominant product, whereas very minor products include the C(6)H(6) isomers 1,1-ethynylmethylallene and 2-ethynyl-1,3-butadiene.     

Hydrogen atom transfer reactions of C2-, C4-, and C6-: bond dissociation energies of linear H-C2n- and H-C2n (n = 1, 2, 3)

The reactions of C2-, C4-, and C6- with D2O and ND3 and of C4- with CH3OH, CH4, and C2H6 have been investigated using guided ion beam tandem mass spectrometry. Hydrogen (or deuterium) atom transfer is the major product channel for each of the reactions. The reaction threshold energies for collisional activation are reported. Several of the reactions exhibit threshold energies in excess of the reaction endothermicity. Potential energy calculations using density functional theory show energy barriers for some of the reactions. Dynamic restrictions related to multiple wells along the reaction path may also contribute to elevated threshold energies. The results indicate that the reactions with D2O have the smallest excess threshold energies, which may therefore be used to derive lower limits on the C-H bond dissociation energies of the C2nH- and C2nH (n = 1-3) linear species. The experimental lower limits for the bond dissociation energies of the neutral radicals to linear products are D0(C2-H) >or= 460 +/- 15 kJ/mol, D0(C4-H) >or= 427 +/- 12 kJ/mol, and D0(C6-H) >or= 405 +/- 11 kJ/mol.     

Overtone-induced dissociation and isomerization dynamics of the hydroxymethyl radical (CH2OH and CD2OH). I. A theoretical study

The dissociation of the hydroxymethyl radical, CH(2)OH, and its isotopolog, CD(2)OH, following the excitation of high OH stretch overtones is studied by quasi-classical molecular dynamics calculations using a global potential energy surface (PES) fitted to ab initio calculations. The PES includes CH(2)OH and CH(3)O minima, dissociation products, and all relevant barriers. Its analysis shows that the transition states for OH bond fission and isomerization are both very close in energy to the excited vibrational levels reached in recent experiments and involve significant geometry changes relative to the CH(2)OH equilibrium structure. The energies of key stationary points are refined using high-level electronic structure calculations. Vibrational energies and wavefunctions are computed by coupled anharmonic vibrational calculations. They show that high OH-stretch overtones are mixed with other modes. Consequently, trajectory calculations carried out at energies about ~3000 cm(-1) above the barriers reveal that despite initial excitation of the OH stretch, the direct OH bond fission is relatively slow (10 ps) and a considerable fraction of the radicals undergoes isomerization to the methoxy radical. The computed dissociation energies are: D(0)(CH(2)OH → CH(2)O + H) = 10,188 cm(-1), D(0)(CD(2)OH → CD(2)O + H) = 10,167 cm(-1), D(0)(CD(2)OH → CHDO + D) = 10,787 cm(-1). All are in excellent agreement with the experimental results. For CH(2)OH, the barriers for the direct OH bond fission and isomerization are: 14,205 and 13,839 cm(-1), respectively.     

Ab initio study of C + H3+ reactions

The reaction C + H3+ --> CH(+) + H2 is frequently used in models of dense interstellar cloud chemistry with the assumption that it is fast, i.e. there are no potential energy barriers inhibiting it. Ab initio molecular orbital study of the triplet CH3+ potential energy surface (triplet because the reactant carbon atom is a ground state triplet) supports this hypothesis. The reaction product is 3 pi CH+; the reaction is to exothermic even though the product is not in its electronic ground state. No path has been found on the potential energy surface for C + H3+ --> CH2(+) + H reaction.     

Infrared spectrum and bonding in uranium methylidene dihydride, CH2=UH2

Uranium atoms activate methane upon ultraviolet excitation to form the methyl uranium hydride CH3-UH, which undergoes alpha-H transfer to produce uranium methylidene dihydride, CH2=UH2. This rearrangement most likely occurs on an excited-quintet potential-energy surface and is followed by relaxation in the argon matrix. These simple U+CH4 reaction products are identified through isotopic substitution (13CH4, CD4, CH2D2) and density functional theory frequency and structure calculations for the strong U-H stretching modes. Relativistic multiconfiguration (CASSCF/CASPT2) calculations substantiate the agostic distorted C1 ground-state structure for the triplet CH2=UH2 molecule. We find that uranium atoms are less reactive in methane activation than thorium atoms. Our calculations show that the CH2=UH2 complex is distorted more than CH2=ThH2. A favorable interaction between the low energy open-shell U(5f) sigma orbital and the agostic hydrogen contributes to the distortion in the uranium methylidene complexes.     

Ab initio investigations of the radical-radical reaction of O(3P) + C3H3

We present ab initio calculations of the reaction of ground-state atomic oxygen [O((3)P)] with a propargyl (C(3)H(3)) radical based on the application of the density-functional method and the complete basis-set model. It has been predicted that the barrierless addition of O((3)P) to C(3)H(3) on the lowest doublet potential-energy surface produces several energy-rich intermediates, which undergo subsequent isomerization and decomposition steps to generate various exothermic reaction products: C(2)H(3)+CO, C(3)H(2)O+H, C(3)H(2)+OH, C(2)H(2)+CHO, C(2)H(2)O+CH, C(2)HO+CH(2), and CH(2)O+C(2)H. The respective reaction pathways are examined extensively with the aid of statistical Rice-Ramsperger-Kassel-Marcus calculations, suggesting that the primary reaction channel is the formation of propynal (CHCCHO)+H. For the minor C(3)H(2)+OH channel, which has been reported in recent gas-phase crossed-beam experiments [H. Lee et al., J. Chem. Phys. 119, 9337 (2003); 120, 2215 (2004)], a comparison on the basis of prior statistical calculations is made with the nascent rotational state distributions of the OH products to elucidate the mechanistic and dynamic characteristics at the molecular level.     

An ab initio based global potential energy surface describing CH5+ --> CH3+ + H2

A full-dimensional, ab initio based potential energy surface (PES) for CH(5)(+), which can describe dissociation is reported. The PES is a precise fit to 36173 coupled-cluster [CCSD(T)] calculations of electronic energies done using an aug-cc-pVTZ basis. The fit uses a polynomial basis that is invariant with respect to permutation of the five H atoms, and thus describes all 120 equivalent minima. The rms fitting error is 78.1 cm(-1) for the entire data set of energies up to 30,000 cm(-1) and a normal-mode analysis of CH(5)(+) also verifies the accuracy of the fit. Two saddle points have been located on the surface as well and compared with previous theoretical work. The PES dissociates correctly to the fragments CH(3)(+) + H(2) and the equilibrium geometry and normal-mode analyses of these fragments are also presented. Diffusion Monte Carlo calculations are done for the zero-point energies of CH(5)(+) (and some isotopologs) as well as for the separated fragments of CH(5)(+), CH(3)(+) + H(2) and those of CH(4)D(+), CH(3)(+) + HD and CH(2)D(+) + H(2). Values of D(0) are reported for these dissociations. A molecular dynamics calculation of CH(4)D(+) dissociation at one total energy is also performed to both validate the applicability of the PES for dynamics studies as well as to test a simple classical statistical prediction of the branching ratio of the dissociation products.     

Experimental and theoretical study of the reactions between neutral vanadium oxide clusters and ethane, ethylene, and acetylene

Reactions of neutral vanadium oxide clusters with small hydrocarbons, namely C2H6, C2H4, and C2H2, are investigated by experiment and density functional theory (DFT) calculations. Single photon ionization through extreme ultraviolet (EUV, 46.9 nm, 26.5 eV) and vacuum ultraviolet (VUV, 118 nm, 10.5 eV) lasers is used to detect neutral cluster distributions and reaction products. The most stable vanadium oxide clusters VO2, V2O5, V3O7, V4O10, etc. tend to associate with C2H4 generating products V(m)O(n)C2H4. Oxygen-rich clusters VO3(V2O5)(n=0,1,2...), (e.g., VO3, V3O8, and V5O13) react with C2H4 molecules to cause a cleavage of the C=C bond of C2H4 to produce (V2O5)(n)VO2CH2 clusters. For the reactions of vanadium oxide clusters (V(m)O(n)) with C2H2 molecules, V(m)O(n)C2H2 are assigned as the major products of the association reactions. Additionally, a dehydration reaction for VO3 + C2H2 to produce VO2C2 is also identified. C2H6 molecules are quite stable toward reaction with neutral vanadium oxide clusters. Density functional theory calculations are employed to investigate association reactions for V2O5 + C2H(x). The observed relative reactivity of C2 hydrocarbons toward neutral vanadium oxide clusters is well interpreted by using the DFT calculated binding energies. DFT calculations of the pathways for VO3+C2H4 and VO3+C2H2 reaction systems indicate that the reactions VO3+C2H4 --> VO2CH2 + H2CO and VO3+C2H2 --> VO2C2 + H2O are thermodynamically favorable and overall barrierless at room temperature, in good agreement with the experimental observations.     

Dissociation of benzene dication [C6H6]2+: exploring the potential energy surface

The singlet potential energy surface for the dissociation of benzene dication has been explored, and its three major dissociation channels have been studied: C6H6(2+) --> C3H3(+) + C3H3(+), C4H3(+) + C2H3(+), and C5H3(+) + CH3(+). The calculated energetics suggest that the products will be formed with considerable translational energy because of the Coulomb repulsion between the charged fragments. The calculated energy release in the three channels shows a qualitative agreement with the experimentally observed kinetic energy release. The formation of certain intermediates is found to be common to the three dissociation channels.