Semiempirical calculation of the CH4 + CH3 · → CH3 · + CH4 radical-reaction potential surfaces in a basis of frontier orbitals

A semiempirical approach is suggested to describe potential-energy surfaces (PESs) of some radical reactions in the zero-differential overlap (ZDO) approximation. An incomplete basis set is used including only frontier (single-filled) radical molecular orbitals (MOs) depending on geometrical parameters. All possible configurations are taken into account. The parameter selection techniques are analyzed. The approach is applied to PES calculation of the CH₄ + CH₃ · CH₃ · + CH₄ reaction. Some points of the PES are verified by a nonempirical method using the perturbation theory and taking into account the correlation energy. The relaxation energies are calculated. The one-center parameters are determined nonempirically from the CH₃ + CH₃ ·, and CH₃ ^– energetics. The two-center parameters are found by modeling the CH₄ CH₃ · + H and C₂H₆ 2CH₃ reactions in the same single-orbital approximation. The energy parameters of the reactions considered are overestimated by 10%, whereas the geometrical parameters are under-estimated by 15%. Further, a comparative analysis of the Hartree-Fock solutions and those including correlation interactions (CIs) is given. The variations in the spin and charge densities on the reaction centers are considered.Institute of Chemical Kinetics and Combustion, Academy of Sciences of the USSR, Siberian Branch, Novosibirsk. Translated from Teoreticheskaya i Éksperimental'naya Khimiya, Vol. 27, No. 4, pp. 499–506, July–August, 1991. Original article submitted February 10, 1989.     

Potential energy surface and reaction mechanism for the ion-molecule reaction: CH4 + CH4+ → CH3 + CH5+

The potential energy surface for the CH₄+CH₄⁺ reaction system has been calculated with the ab initio method. A stable complex, responsible for the complex mechanism, has been found but is hard to reach. Each of the two direct mechanisms, hydrogen transfer and proton transfer, has been shown to consist of a combination of electron transfer and hydrogen atom transfer processes.     

A photoionization study of the ion-neutral complexes CH3CH +CH 3 · CH2CH3] and CH3CH2CH+CH 3 · CH3 in the gas phase: Formation,H-transfer and C—C bond formation between partners,and channeling of energy into dissociation

Photoionization mass spectrometry was used to investigate the dynamics of ion-neutral complex-mediated dissociations of the n-pentane ion (1). Reinterpretation of previous data demonstrates that a fraction of ions 1 isomerizes to the 2-methylbutane ion (2) through the complex CH₃CH⁺CH ₃ ^· CH₂CH₃ (3), but not through CH₃CH⁺CH₂CH ₃ ^· CH₃ (4). The appearance energy for C₃Hin ₇ ⁺ formation from 1 is 66 kJ mol^?1 below that expected for the formation of n-C₃H ₇ ⁺ and just above that expected for formation of i-C₃H ₇ ⁺ . This demonstrates that the H shift that isomerizes C₃H ₇ ⁺ is synchronized with bond cleavage at the threshold for dissociation to that product. It is suggested that ions that contain n-alkyl chains generally dissociate directly to more stable rearranged carbenium ions. Ethane elimination from 3 is estimated to be about seven times more frequent than is C-C bond formation between the partners in that complex to form 2, which demonstrates a substantial preference in 3 for H abstraction over C-C bond formation. In 1 → CH₃CH⁺CH₂CH₃ + CH₃ by direct cleavage of the C1–C2 bond, the fragments part rapidly enough to prevent any reaction between them. However, 1 → 2 → 4 → C4H ₈ ⁺ + CH₄ occurs in this same energy range. Thus some of the potential energy made available by the isomerization of n-C₄H₉ in 1 is specifically channeled into the coordinate for dissociation. In contrast, analogous formation of 3 by 1 → 3 is predominantly followed by reaction between the electrostatically bound partners.     

Capture and dissociation in the complex-forming CH + H2 → CH2 + H, CH + H2 reactions

The rate coefficients for the capture process CH + H(2)→ CH(3) and the reactions CH + H(2)→ CH(2) + H (abstraction), CH + H(2) (exchange) have been calculated in the 200-800 K temperature range, using the quasiclassical trajectory (QCT) method and the most recent global potential energy surface. The reactions, which are of interest in combustion and in astrochemistry, proceed via the formation of long-lived CH(3) collision complexes, and the three H atoms become equivalent. QCT rate coefficients for capture are in quite good agreement with experiments. However, an important zero point energy (ZPE) leakage problem occurs in the QCT calculations for the abstraction, exchange and inelastic exit channels. To account for this issue, a pragmatic but accurate approach has been applied, leading to a good agreement with experimental abstraction rate coefficients. Exchange rate coefficients have also been calculated using this approach. Finally, calculations employing QCT capture/phase space theory (PST) models have been carried out, leading to similar values for the abstraction rate coefficients as the QCT and previous quantum mechanical capture/PST methods. This suggests that QCT capture/PST models are a good alternative to the QCT method for this and similar systems.     

Dynamics of the O(3P) + CH4 → OH + CH3 reaction is similar to that of a triatomic reaction

The O((3)P) + CH(4) reaction has been investigated using the quasi-classical trajectory (QCT) method and an ab initio pseudotriatomic potential energy surface (PES). This has been mainly motivated by very recent experiments which support the reliability of the triatomic modeling even at high collision energy ( = 64 kcal mol(-1)). The QCT results agree rather well with the experiments (translational and angular distributions of products); i.e., the ab initio pseudotriatomic modeling "captures" the essence of the reaction dynamics, although the PES was not optimized for high E(col). Furthermore, similar experiments on the O((3)P) + CD(4) reaction at moderate E(col) (12.49 kcal mol(-1)) have also been of a large interest here and, under these softer reaction conditions, the QCT method leads to results which are almost in quantitative agreement with experiments. The utility of the ab initio pseudotriatomic modeling has also been recognized for other analogous systems (X + CH(4)) but with very different PESs.     

Crossed-beams studies of the dynamics of the H-atom abstraction reaction, O(3P) + CH4 → OH + CH3, at hyperthermal collision energies

The H-atom abstraction reaction, O((3)P) + CH(4) → OH + CH(3), has been studied at a hyperthermal collision energy of 64 kcal mol(-1) by two crossed-molecular-beams techniques. The OH products were detected with a rotatable mass spectrometer employing electron-impact ionization, and the CH(3) products were detected with the combination of resonance-enhanced multiphoton ionization (REMPI) and time-sliced ion velocity-map imaging. The OH products are mainly formed through a stripping mechanism, in which the reagent O atom approaches the CH(4) molecule at large impact parameters and the OH product is scattered in the forward direction: roughly the same direction as the reagent O atoms. Most of the available energy is partitioned into product translation. The dominance of the stripping mechanism is a unique feature of such H-atom abstraction reactions at hyperthermal collision energies. In the hyperthermal reaction of O((3)P) with CH(4), the H-atom abstraction reaction pathway accounts for 70% of the reactive collisions, while the H-atom elimination pathway to produce OCH(3) + H accounts for the other 30%.     

RISM-SCF study of the free-energy profile of the Menshutkin-type reaction NH3+CH3Cl→NH3CH3++Cl− in aqueous solution

The free-energy profile for the Menshutkin-type reaction NH₃ + CH₃Cl → NH₃CH₃ ⁺ + Cl⁻ in aqueous solution is studied using the RISM-SCF method. The effect of electron correlation on the free-energy profile is estimated by the RISM-MP2 method at the HF optimized geometries along the reaction coordinate. Solvation was found to have a large influence on the vibrational frequencies at the reactant, transition state and product; these vibrational frequencies are utilized to calculate the zero-point energy correction of the free-energy profile. The computed barrier height and reaction exothermicity are in reasonable agreement with those of experiment and previous calculations. The change of solvation structure along the reaction path is represented by radial distribution functions between solute-solvent atomic sites. The mechanisms of the reaction are discussed from the view points of solute electronic and solvation structures. Received: 26 June 1998/Accepted: 28 August 1998 / Published online: 2 November 1998     

Reduced dimensionality spin-orbit dynamics of CH3 + HCl ⇌ CH4 + Cl on ab initio surfaces

A reduced dimensionality quantum scattering method is extended to the study of spin-orbit nonadiabatic transitions in the CH(3) + HCl ? CH(4) + Cl((2)P(J)) reaction. Three two-dimensional potential energy surfaces are developed by fitting a 29 parameter double-Morse function to CCSD(T)/IB//MP2/cc-pV(T+d)Z-dk ab initio data; interaction between surfaces is described by geometry-dependent spin-orbit coupling functions fit to MCSCF/cc-pV(T+d)Z-dk ab initio data. Spectator modes are treated adiabatically via inclusion of curvilinear projected frequencies. The total scattering wave function is expanded in a vibronic basis set and close-coupled equations are solved via R-matrix propagation. Ground state thermal rate constants for forward and reverse reactions agree well with experiment. Multi-surface reaction probabilities, integral cross sections, and initial-state selected branching ratios all highlight the importance of vibrational energy in mediating nonadiabatic transition. Electronically excited state dynamics are seen to play a small but significant role as consistent with experimental conclusions.     

Ten-Dimensional Quantum Dynamics Study of H+CH3D→H2+CH2D Reaction

The mode selectivity of the H+CH3D→H2+CH2D reaction was studied using a recently developed ten-dimensional time-dependent wave packet method.The reac-tion dynam...     

Mechanistic and kinetic study of CF3CH═CH2 + OH reaction

The potential energy surfaces of the CF(3)CH═CH(2) + OH reaction have been investigated at the BMC-CCSD level based on the geometric parameters optimized at the MP2/6-311++G(d,p) level. Various possible H (or F)-abstraction and addition/elimination pathways are considered. Temperature- and pressure-dependent rate constants have been determined using Rice-Ramsperger-Kassel-Marcus theory with tunneling correction. It is shown that IM1 (CF(3)CHCH(2)OH) and IM2 (CF(3)CHOHCH(2)) formed by collisional stabilization are major products at 100 Torr pressure of Ar and in the temperature range of T < 700 K (at P = 700 Torr with N(2) as bath gas, T ≤ 900 K), whereas CH(2)═CHOH and CF(3) produced by the addition/elimination pathway are the dominant end products at 700-2000 K. The production of CF(3)CHCH and CF(3)CCH(2) produced by hydrogen abstractions become important at T ≥ 2000 K. The calculated results are in good agreement with available experimental data. The present theoretical study is helpful for the understanding the characteristics of the reaction of CF(3)CH═CH(2) + OH.     

Ab initio potential energy surface and quantum dynamics for the H + CH4 → H2 + CH3 reaction

A new full-dimensional potential energy surface for the title reaction has been constructed using the modified Shepard interpolation scheme. Energies and derivatives were calculated using the UCCSD(T) method with aug-cc-pVTZ and 6-311++G(3df,2pd) basis sets, respectively. A total number of 30,000 data points were selected from a huge number of molecular configurations sampled by trajectory method. Quantum dynamical calculations showed that the potential energy surface is well converged for the number of data points for collision energy up to 2.5 eV. Total reaction probabilities and integral cross sections were calculated on the present surface, as well as on the ZBB3 and EG-2008 surfaces for the title reaction. Satisfactory agreements were achieved between the present and the ZBB3 potential energy surfaces, indicating we are approaching the final stage to obtain a global potential energy surface of quantitative accuracy for this benchmark polyatomic system. Our calculations also showed that the EG-2008 surface is less accurate than the present and ZBB3 surfaces, particularly in high energy region.     

Communication: probing the entrance channels of the X+CH4→HX+CH3 (X = F, Cl, Br, I) reactions via photodetachment of X(-)-CH4

The entrance channel potentials of the prototypical polyatomic reaction family X + CH(4) → HX + CH(3) (X = F, Cl, Br, I) are investigated using anion photoelectron spectroscopy and high-level ab initio electronic structure computations. The pre-reactive van der Waals (vdW) wells of these reactions are probed for X = Cl, Br, I by photodetachment spectra of the corresponding X(-)-CH(4) anion complex. For F-CH(4), a spin-orbit splitting (～1310 cm(-1)) much larger than that of the F atom (404 cm(-1)) was observed, in good agreement with theory. This showed that in the case of the F-CH(4) system the vertical transition from the anion ground state to the neutral potentials accesses a region between the vdW valley and transition state of the early-barrier F + CH(4) reaction. The doublet splittings observed in the other halogen complexes are close to the isolated atomic spin-orbit splittings, also in agreement with theory.     

The analysis of high-resolution spectra of the ν2 and ν5 bands of CH335Cl and CH337Cl

The ν₂ and ν₅ bands of CH³₃⁵Cl and CH³₃⁷Cl between 1300 and 1600 cm⁻¹ have been analysed using a Fourier transform spectrum with 0.006 cm⁻¹ resolution. For CH³₃⁵Cl, the microwave data and 1200 lines from the IR spectrum with J⩽ 50 were fitted with an overall r.m.s. error of 0.00079 cm⁻¹ using the method of predicative observations. A similar fit for 900 lines of CH³₃⁷Cl gave an overall r.m.s. error of 0.00055 cm⁻¹, providing erroneous microwave data on the ν₅ level are omitted. Improved molecular constants are reported for both isotopic species. As expected, the values for ν₂ and ν₅ are little affected by chlorine isotopic substitution.     

Accurate ab initio potential energy surface, thermochemistry, and dynamics of the Cl(2P, 2P(3/2) + CH4 → HCl + CH3 and H + CH3Cl reactions

We report a high-quality, ab initio, full-dimensional global potential energy surface (PES) for the Cl((2)P, (2)P(3/2)) + CH(4) reaction, which describes both the abstraction (HCl + CH(3)) and substitution (H + CH(3)Cl) channels. The analytical PES is a least-squares fit, using a basis of permutationally invariant polynomials, to roughly 16,000 ab initio energy points, obtained by an efficient composite method, including counterpoise and spin-orbit corrections for the entrance channel. This composite method is shown to provide accuracy almost equal to all-electron CCSD(T)/aug-cc-pCVQZ results, but at much lower computational cost. Details of the PES, as well as additional high-level benchmark characterization of structures and energetics are reported. The PES has classical barrier heights of 2650 and 15,060 cm(-1) (relative to Cl((2)P(3/2)) + CH(4)(eq)), respectively, for the abstraction and substitution reactions, in good agreement with the corresponding new computed benchmark values, 2670 and 14,720 cm(-1). The PES also accurately describes the potential wells in the entrance and exit channels for the abstraction reaction. Quasiclassical trajectory calculations using the PES show that (a) the inclusion of the spin-orbit corrections in the PES decreases the cross sections by a factor of 1.5-2.5 at low collision energies (E(coll)); (b) at E(coll) ≈ 13,000 cm(-1) the substitution channel opens and the H/HCl ratio increases rapidly with E(coll); (c) the maximum impact parameter (b(max)) for the abstraction reaction is ~6 bohr; whereas b(max) is only ~2 bohr for the substitution; (d) the HCl and CH(3) products are mainly in the vibrational ground state even at very high E(coll); and (e) the HCl rotational distributions are cold, in excellent agreement with experiment at E(coll) = 1280 cm(-1).     

Quasi-classical trajectory studies of the stereodynamics of the reaction O + HCl → ClO + H

The stereodynamics of the O + HCl → ClO + H reaction are investigated by quasi-classical trajectory (QCT) method. The calculations are carried out on the ground 1 ¹ A′ potential energy surface (PES). The orientation and alignments of the product rotational angular momentum for the title reaction are reported. The influence of collision energy on the product vector properties is also studied in the present work. Four (2π/σ)(dσ₀₀/dω _t ), (2π/σ)(dσ₂₀/dω _t ), (2π/σ)(dσ₂₂₊/dω _t ), and (2π / σ)(dσ₂₁₋/dω _t ), and have been calculated in the center of mass frame.     

Impact of Water Molecules on the Isomerization of CH3S（OH）CH2 to CH3S（O）CH3： A Computational Investigation

The isomerization of CH3S（OH）CH2 to CH3S（O）CH3 in the absence and presence of water has been investigated at the G3XMP2//B3LYP/6-311 ＋ G（2df, p） level. The naked isomerization, the reaction without water, gives the high barrier height （21.56 kcal.mol^-1）. Three models are constructed to describe the water influence on the isomerization, that is, water molecules are the catalyst and the microsolvation, and water molecules act as the catalyst and microsolvation simultaneously. Our results show that the isomerization barrier heights of CH3S（OH）CH2 to CH3S（O）CH3 are reduced by 12.32, 11.04, and 7.80 kcal.mol^-1, respectively, when one, two, and three water molecules are performed as catalyst, in contrast to the naked isomerization. Moreover, the rate constants of the isomerization are calculated using the transition state theory with the Wigner tunneling correction over the temperature range of 240-425 K. We find that the rate constant of a single water molecule as the catalyst is 1.58 times larger than the naked isomerization at 325 K, whereas it is slower by 6 orders of magnitude when water molecule serves as the microsolvation at 325 K, compared to naked reaction. So the water-catalyzed isomerization of CH3S（OH）CH2 to CH3S（O）CH3 is predicted to be the key role in lowering the activation energy. The isomerization involving water molecules acting as mierosolvation is unfavorable under atmospheric conditions.     

Theoretical Studies of the Reaction Mechanisms of CH3S＋NO2

The potential energy surface for the CH3S NO2 reaction has been studied using the ab initio G3(MP2) method. A variety of possible complexes and saddle points along the minimum energy reaction paths have been characterized at UMP2 (full)/6-31G(d) level. The calculations reveal dominating reaction mechanisms of the title reaction: CH3S NO2 firstly produce intermediate CH3SONO,then break up into CH3SO NO. The results are valuable to understand the atmospheric sulfur compounds oxidation mechanism.