Ab initio potential energy surface, variational transition state theory, and quasiclassical trajectory studies of the F+CH4-->HF+CH3 reaction

In this work we present a study of the F+CH(4)-->HF+CH(3) reaction (DeltaHdegrees(298 K)=-32.0 kcal mol(-1)) using different methods of the chemical reaction theory. The ground potential energy surface (PES) is characterized using several ab initio methods. Full-dimensional rate constants have been calculated employing the variational transition state theory and using directly ab initio data. A triatomic analytical representation of the ground PES was derived from ab initio points calculated at the second- and fourth-order M?ller-Plesset levels with the 6-311+G(2df,2pd) basis set, assuming the CH(3) fragment to be a 15 a.m.u. pseudoatom in the fitting process. This is suggested from experiments that indicate that the methyl group is uncoupled to the reaction coordinate. A dynamics study by means of the quasiclassical trajectory (QCT) method and employing this analytical surface was also carried out. The experimental data available on the HF internal states distributions are reproduced by the QCT results. Very recent experimental information about the reaction stereodynamics is also borne out by our QCT calculations. Comparisons with the benchmark F+H(2) and analogous Cl+CH(4) reactions are established throughout.     

Quasiclassical trajectory study of the F + CH4 reaction dynamics on a dual-level interpolated potential energy surface

An ab initio interpolated potential energy surface (PES) for the F + CH4 reactive system has been constructed using the interpolation method of Collins and co-workers. The ab initio calculations have been performed using second-order M?ller-Plesset (MP2) perturbation theory to build the initial PES. Scaling all correlation (SAC) methodology has been employed to improve the ab initio calculations and to construct a dual-level PES. Using this PES, a detailed quasiclassical trajectory study of integral and differential cross sections, product rovibrational populations and internal energy distributions has been carried out for the F + CH4 and F + CD4 reactions and the theoretical results have been compared with the available experimental data.     

Potential energy surface, kinetics, and dynamics study of the Cl+CH4-->HCl+CH3 reaction

A modified and recalibrated potential energy surface for the gas-phase Cl+CH4-->HCl+CH3 reaction is reported and tested. It is completely symmetric with respect to the permutation of the four methane hydrogen atoms and is calibrated with respect to updated experimental and theoretical stationary point properties and experimental forward thermal rate constants. From the kinetics point of view, the forward and reverse thermal rate constants and the activation energies were calculated using the variational transition-state theory with semiclassical transmission coefficients over a wide temperature range of 150-2500 K. The theoretical results reproduce the available experimental data, with a small curvature of the Arrhenius plot which indicates the role of tunneling in this hydrogen abstraction reaction. A dynamics study was also performed on this PES using quasiclassical trajectory (QCT) calculations, including corrections to avoid zero-point energy leakage along the trajectories. First, we found a noticeable internal energy in the coproduct methyl radical, both in the ground-state [CH4 (v=0)] and vibrationally excited [CH4 (v=1)] reactions. This CH3 internal energy was directly precluded in some experiments or oversimplified in previous theoretical studies using pseudotriatomic models. Second, our QCT calculations give HCl rotational distributions slightly hotter than those in experiment, but correctly describing the experimental trend of decreasing the HCl product rotation excitation in going from HCl (v'=0) to HCl (v'=1) for the CH4 (v=1) reaction. Third, the state specific scattering distributions present a reasonable agreement with experiment, although they tend to make the reaction more forward and backward scattered than found experimentally probably because of the hotter rotational distribution and the deficiencies of the QCT methods.     

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.     

辛准经典轨迹法研究Cl+H2反应

用辛准经典轨迹法模拟了Cl+H2反应在mBW2势能面上的动力学行为。研究了各种初始条件下的反应碰撞截面,产物的能量分配,角度分布和态分布。另外,我们还比较了反应物的三种能量形式(平动能,转动能和振动能)对反庆的有效性。     

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).     

Theoretical study of the dynamics of Cl + O3 reaction I. Ab initio potential energy surface and quasiclassical trajectory results

We present a global full dimensional potential energy surface (PES) for the Cl + O(3)→ ClO + O(2) reaction, which is an elementary step in a catalytic cycle that leads to the destruction of ozone in the stratosphere. The PES is constructed by interpolation of quantum chemistry data using the method developed by Collins and co-workers. Ab initio data points (energy, gradients and Hessian matrix elements) have been calculated at the UQCISD/aug-cc-pVDZ (unrestricted quadratic configuration interaction with single and double excitations) level of theory. The ab initio calculations predict a markedly non-coplanar (dihedral angle of 80°) transition state for the reaction, located very early in the reactant valley and slightly below the energy of the reactants as long as the spin-orbit splitting is neglected. Quasiclassical trajectory (QCT) calculations have been carried out at several collision energies to investigate the reaction dynamics. The QCT excitation function shows no threshold, displays a minimum at a collision energy of 2.5 kcal mol(-1), and then increases monotonically at larger collision energies. This behaviour is consistent with a barrierless reaction dominated by an oxygen-abstraction mechanism. The calculated product vibrational distributions (strongly inverted for ClO) and rate constants are compared with experimental determinations. Differential cross sections (DCS) summed over all final states are found to be in fairly good agreement with those derived from crossed molecular beam experiments.     

Exit interaction effect on nascent product state distribution of O(1D)+N2O-->NO+NO

We have determined the rotational state distributions of NO(v'=0,1,2) products produced from the reaction O(1D)+N2O. This is the first full characterization of the product rotational distribution of this reaction. The main part of each rotational distribution (up to j' approximately 80) has rotational temperature approximately 20,000 K and all these distributions are quite near to those predicted by the phase space theory (PST). This observation and previously reported vibrational distribution indicate that the most part of the energy partitioning of the reaction products is at least apparently statistical although the intermediate of this reaction is not so stable as to ensure the long lifetime. On the other hand, the distributions in the high rotational levels (j'=80-100) are found to decrease more sharply as j' increases than the PST predictions. The origin of the observed decrease of the distribution is discussed with quasiclassical trajectory (QCT) calculations on a five-dimensional ab initio potential energy surface (PES). The observed near-statistical distribution and the sharp decrease in the high-j' levels are well reproduced by a "half-collision" QCT calculation, where statistical distribution at the reaction intermediate is assumed. This agreement shows the rotation-translation interaction in the exit region has an effect of yielding small high-j' populations. However, a little bias of the calculated distribution toward lower rotational excitation than the observed one indicates that the combination of the statistical intermediate and the exit interaction on the current PES does not completely describe the real system. It is suggested that the reaction intermediate is generated with the distribution which is close to statistical but a little biased toward yielding high-j' products, and that the interaction in the exit region of the PES results in the sharp decrease in the high-j' levels.     

Ab initio and direct quasiclassical-trajectory study of the Cl + CH4-->HCl + CH3 reaction

We present an electronic structure and dynamics study of the Cl + CH(4)--> HCl + CH(3) reaction. We have characterized the stationary points of the ground-state potential-energy surface using various electronic structure methods and basis sets. Our best calculations, CCSD(T) extrapolated to the complete basis-set limit based on geometries and harmonic frequencies obtained at the CCSD(T)/aug-cc-pvtz level, are in agreement with the experimental reaction energy and indirect measurements of the barrier height. Using ab initio information, we have reparametrized a semiempirical Hamiltonian so that the predictions of the improved Hamiltonian agree with the higher-level calculations in various regions of the potential-energy surface. This improved semiempirical Hamiltonian is then used to propagate quasiclassical trajectories and characterize the reaction dynamics. The good agreement of the calculated HCl rotational and angular distributions with the experiment indicates that reparametrizing semiempirical Hamiltonians is a promising approach to derive accurate potential-energy surfaces for polyatomic reactions. However, excessive energy leakage from the initial vibrational energy of the CH(4) molecule to the reaction coordinate in the trajectory calculations calls into question the suitability of the standard quasiclassical-trajectory method to describe energy partitioning in polyatomic reactions.     

Quantum mechanical and quasiclassical trajectory scattering calculations for the C(1D) + H2 reaction on the second excited 1 1A" potential energy surface

Time-independent quantum mechanical (QM) and quasiclassical trajectory (QCT) scattering calculations have been carried out for the C(1D) + H2 --> CH + H reaction at a collision energy of 80 meV on a newly developed ab initio potential energy surface [B. Bussery-Honvault et al., Phys. Chem. Chem. Phys. 7, 1476 (2005)] of 1 1A" symmetry, corresponding to the second singlet state 1 1B1 of CH2. A general good agreement has been found between the QM and QCT rotational distributions and differential cross sections (DCSs). In both cases, DCSs are strongly peaked in the forward direction with a small contribution in the backward direction in contrast with those obtained on the 1 1A' surface, which are nearly symmetric. Rotational distributions obtained on the 1 1A" surface are somewhat colder than those calculated on the 1 1A' surface. The specific dynamics and the contribution of the 1 1A" surface to the overall reactivity of this system are discussed.     

Probing the structure of CH5+ by dissociative charge exchange

Dissociative charge exchange of CH5+ with Cs, coupled with quasiclassical trajectory calculations on an ab initio PES for CH5, has been used to probe the structure of the CH5+ cation. Product kinetic energy release distributions and branching ratios for CH5 --> CH4 + H and CH5 --> CH3 + H2 have been compared. The agreement of the product branching ratios provides evidence for the fluxional nature of CH5+.     

Quasi-classical Trajectory Study of F+H2O→HF+OH Reaction: Influence of Barrier Height,Reactant Rotational Excitation,and Isotopic Substitution (cited: 4)

The reaction dynamics of the F+H₂O/D₂O→HF/DF+OH/OD are investigated on an ac-curate potential energy surface (PES) using a quasi-classical trajectory method. For bothisotopomers, the hydrogen/deuterium abstraction reaction is dominated by a direct rebound mechanism over a very low “reactant-like” barrier, which leads to a vibrationally hot HF/DF product with an internally cold OH/OD companion. It is shown that the lowered reaction barrier on this PES, as suggested by high-level ab initio calculations, leads to a much better agreement with the experimental reaction cross section, but has little impact on the product state distributions and mode selectivity. Our results further indicate that rotational exci-tation of the H₂O reactant leads to significant enhancement of the reactivity, suggesting a strong coupling with the reaction coordinate.     

Searching for resonances in the reaction Cl+CH4-->HCl+CH3: quantum versus quasiclassical dynamics and comparison with experiments

A quantum-mechanical (QM) and quasiclassical trajectory (QCT) study was performed on the title reaction, using a pseudotriatomic ab initio based surface. Probabilities and integral cross sections present some clear peaks versus the collision energy E(col), which we assign to Feshbach resonances of the transition state, where the light H atom oscillates between the heavy Cl and CH(3) groups. For ground-state reactants, reactivity is essentially of quantum origin (QCT observables and oscillations are smaller, or much smaller, than QM ones), and the calculated integral cross section and product distributions are in reasonable agreement with the experiment. The reaction occurs through an abstraction mechanism, following both a direct and an indirect mechanism. The quasiclassical trajectory calculations show the participation of a short-lived collision complex in the microscopic reaction mechanism. Finally, QCT differential cross sections of Cl+CH(4)-->HCl (nu(')=0 and 1)+CH(3) oscillate versus E(col), whereas experimentally this only occurs for HCl (nu(')=1). This theoretical result and other oscillating properties found here could, however, be related to the existence of a Feshbach resonance for the production of HCl (nu(')=1), as suggested by experimentalists.     

Comparison of experimental time-of-flight spectra of the HF products from the F+H2 reaction with exact quantum mechanical calculations

High resolution HF product time-of-flight spectra measured for the reactive scattering of F atoms from n-H2(p-H2) molecules at collision energies between 69 and 81 meV are compared with exact coupled-channel quantum mechanical calculations based on the Stark-Werner ab initio ground state potential energy surface. Excellent agreement between the experimental and computed rotational distributions is found for the HF product vibrational states v'=1 and v'=2. For the v'=3 vibrational state the agreement, however, is less satisfactory, especially for the reaction with p-H2. The results for v'=1 and v'=2 confirm that the reaction dynamics for these product states is accurately described by the ground electronic state 1 (2)A' potential energy surface. The deviations for HF(v'=3, j' > or =2) are attributed to an enhancement of the reaction resulting from the 25% fraction of excited ((2)P(12)) fluorine atoms in the reactant beam.     

Evidence for excited spin-orbit state reaction dynamics in F+H2: theory and experiment

We describe fully quantum, time-independent scattering calculations of the F+H2-->HF+H reaction, concentrating on the HF product rotational distributions in v'=3. The calculations involved two new sets of ab initio potential energy surfaces, based on large basis set, multireference configuration-interaction calculations, which are further scaled to reproduce the experimental exoergicity of the reaction. In addition, the spin-orbit, Coriolis, and electrostatic couplings between the three quasidiabatic F+H2 electronic states are included. The calculated integral cross sections are compared with the results of molecular beam experiments. At low collision energies, a significant fraction of the reaction is due to Born-Oppenheimer forbidden, but energetically allowed reaction of F in its excited (2P 1/2) spin-orbit state. As the collision energy increases, the Born-Oppenheimer allowed reaction of F in its ground (2P 3/2) spin-orbit state rapidly dominates. Overall, the calculations agree reasonably well with the experiment, although there remains some disagreement with respect to the degree of rotational excitation of the HF(v'=3) products as well as with the energy dependence of the reactive cross sections at the lowest collision energies.     

Experimental and theoretical differential cross sections for the N(2D) + H2 reaction

In this paper, we report a combined experimental and theoretical study on the dynamics of the N(2D) + H2 insertion reaction at a collision energy of 15.9 kJ mol(-1). Product angular and velocity distributions have been obtained in crossed beam experiments and simulated by using the results of quantum mechanical (QM) scattering calculations on the accurate ab initio potential energy surface (PES) of Pederson et al. (J. Chem. Phys. 1999, 110, 9091). Since the QM calculations indicate that there is a significant coupling between the product angular and translational energy distributions, such a coupling has been explicitly included in the simulation of the experimental results. The very good agreement between experiment and QM calculations sustains the accuracy of the NH2 ab initio ground state PES. We also take the opportunity to compare the accurate QM differential cross sections with those obtained by two approximate methods, namely, the widely used quasiclassical trajectory calculations and a rigorous statistical method based on the coupled-channel theory.     

Quasiclassical trajectory study of the O(3P) + CH4 --> OH + CH3 reaction with a specific reaction parameters semiempirical Hamiltonian

We present a theoretical study of the O(3P) + CH4 --> OH + CH3 reaction using electronic structure, kinetics, and dynamics calculations. We calculate a grid of ab initio points at the PMP2/AUG-cc-pVDZ level to characterize the potential energy surface in regions of up to 1.3 eV above reagents. This grid of ab initio points is used to derive a set of specific reaction parameters (SRP) for the MSINDO semiempirical Hamiltonian. The resulting SRP-MSINDO Hamiltonian improves the quality of the standard Hamiltonian, particularly in regions of the potential energy surface beyond the minimum-energy reaction path. Quasiclassical-trajectory calculations are used to study the reaction dynamics with the original and the improved MSINDO semiempirical Hamiltonians, and a prior surface. The SRP-MSINDO semiempirical Hamiltonian yields OH rotational distributions in agreement with experimental results, improving over the results of the other surfaces. Thermal rate constants estimated with Variational Transition State Theory using the SRP-MSINDO Hamiltonian are also in agreement with experiments. Our results indicate that reparametrized semiempirical Hamiltonians are a good alternative to generating potential energy surfaces for accurate dynamics studies of polyatomic reactions.     

用神经网络方法拟合的反应体系H＋CH4←→H2＋CH3的一个全域势能面

利用神经网络力法,基于47783个高精度从头算能量点构建了反应体系H＋CH4←→H2＋CH3的一个全域势能面．通过大最的准经典轨线以及量子散射计算测试了势能面的收敛性质．这个势能面对于拟合过程以及从头算点的数目都已经完全收敛,拟合误差很小县比Shepard插值的势能面计算速度更快,代表了此标志性多原子反应体系最好的势能面．     

Understanding the complex dissociation dynamics of energy selected dichloroethylene ions: neutral isomerization energies and heats of formation by imaging photoelectron-photoion coincidence

The dissociative photoionization of 1,1-C(2)H(2)Cl(2), (E)-1,2-C(2)H(2)Cl(2), and (Z)-1,2-C(2)H(2)Cl(2) has been investigated at high energy and mass resolution using the imaging photoelectron photoion coincidence instrument at the Swiss Light Source. The asymmetric Cl-atom loss ion time-of-flight distributions were fitted to obtain the dissociation rates in the 10(3) s(-1) < k < 10(7) s(-1) range as a function of the ion internal energy. The results, supported by ab initio calculations, show that all three ions dissociate to the same C(2v) symmetry ClC═CH(2)(+) product ion. The 0 K onset energies thus establish the relative heats of formation of the neutral isomers, that is, the isomerization energies. The experimental rate constants, k(E), as well as ab initio calculations indicate an early isomerization transition state and no overall reverse barrier to dissociation. The major high energy channels are the parallel HCl loss and the sequential ClC═CH(2)(+) → HCCH(+) + Cl process, the latter in competition with a ClC═CH(2)(+) → ClCCH(+) + H reaction. A parallel C(2)H(2)Cl(2)(+) → C(2)HCl(2)(+) + H channel also weakly asserts itself. The 0 K onset energy for the sequential Cl loss reaction suggests no barrier to the production of the most stable acetylene ion product; thus the sequential Cl-atom loss is preceded by a ClC═CH(2)(+) → HC(Cl)CH(+) reorganization step with a barrier lower than that of the second Cl-atom loss. The breakdown diagram corresponding to this sequential dissociation reveals the internal energy distribution of the first C(2)H(2)Cl(+) daughter ion, which is determined by the kinetic energy release in the first, Cl loss reaction at high excess energies. At low kinetic energy release, this distribution corresponds to the predicted two translational degrees of freedom, whereas at higher energies, the excess energy partitioning is characteristic of only one translational degree of freedom. New Δ(f)H(o)(298K) of 3.7, 2.5, and 0.2 ± 1.75 kJ mol(-1) are proposed for 1,1-C(2)H(2)Cl(2), (E)-1,2-C(2)H(2)Cl(2), and (Z)-1,2-C(2)H(2)Cl(2), respectively, and the proton affinity of ClCCH is found to be 708.6 ± 2.5 kJ mol(-1).     

An ab initio spin-orbit-corrected potential energy surface and dynamics for the F + CH4 and F + CHD3 reactions

We report an analytical ab initio three degrees of freedom (3D) spin-orbit-correction surface for the entrance channel of the F + methane reaction obtained by fitting the differences between the spin-orbit (SO) and non-relativistic electronic ground state energies computed at the MRCI+Q/aug-cc-pVTZ level of theory. The 3D model surface is given in terms of the distance, R(C-F), and relative orientation, Euler angles ? and θ, of the reactants treating CH(4) as a rigid rotor. The full-dimensional (12D) "hybrid" SO-corrected potential energy surface (PES) is obtained from the 3D SO-correction surface and a 12D non-SO PES. The SO interaction has a significant effect in the entrance-channel van der Waals region, whereas the effect on the energy at the early saddle point is only ～5% of that at the reactant asymptote; thus, the SO correction increases the barrier height by ～122 cm(-1). The 12D quasiclassical trajectory calculations for the F + CH(4) and F + CHD(3) reactions show that the SO effects decrease the cross sections by a factor of 2-4 at low collision energies and the effects are less significant as the collision energy increases. The inclusion of the SO correction in the PES does not change the product state distributions.