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.     

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.     

Intermolecular interaction potentials for the Ar-C2H2, Kr-C2H2, and Xe-C2H2 weakly bound complexes: information from molecular beam scattering, pressure broadening coefficients, and rovibrational spectroscopy

Integral cross sections and pressure broadening coefficients have been measured for the acetylene-krypton complex, by molecular beam scattering and by high resolution IR spectroscopy, respectively. A new potential energy surface (PES) is proposed to describe structure and dynamical properties of this prototypical weakly bound complex. The PES has been parametrized exploiting a novel atom-bond pairwise additive scheme and has been fitted to the experimental data. A similar PES has been obtained for the acetylene-xenon system by a proper scaling of the interaction parameters of the krypton case, based on empirical considerations. These PESs together with that recently proposed by the same authors [J. Phys. Chem. 109, 8471 (2005)] for the acetylene-argon case have been employed for close coupling calculations of the pressure broadening cross sections and for a characterization of the rovibrational structure of the complexes.     

Exploring the accuracy level of new potential energy surfaces for the F + HD reactions: from exact quantum rate constants to the state-to-state reaction dynamics

Exact quantum reactive scattering calculations in the collision energy range 1-250 meV have been carried out for both the isotopic product channels of the title system. The dynamical studies compares an ab initio potential energy surface (PES) recently appeared in the literature (J. Chem. Phys., 2008, 129, 011103) with other phenomenological PESs. Vibrational branching ratios, cross sections and rate constants are presented and compared with molecular beam scattering experiments as well as with chemical kinetics data. In particular, the agreement of the vibrational branching ratios with experimental measurements is improved with respect to previous studies on other PESs, mainly because of the presence of a broad peak in the HF(v' = 3) integral cross section completely absent in the previous simulations. This feature, observed by molecular beam experiments, is the fingerprint of a new reaction mechanism operative in the dynamics described by the new PES. A conjecture for its origin, able to explain many of its characteristic aspects, is analyzed and discussed.     

Quasiclassical trajectory study of the Cl+CH4 reaction dynamics on a quadratic configuration interaction with single and double excitation interpolated potential energy surface

An ab initio interpolated potential energy surface (PES) for the Cl+CH(4) reactive system has been constructed using the interpolation method of Collins and co-workers [J. Chem. Phys. 102, 5647 (1995); 108, 8302 (1998); 111, 816 (1999); Theor. Chem. Acc. 108, 313 (2002)]. The ab initio calculations have been performed using quadratic configuration interaction with single and double excitation theory to build the PES. A simple scaling all correlation technique has been used to obtain a PES which yields a barrier height and reaction energy in good agreement with high level ab initio calculations and experimental measurements. Using these interpolated PESs, a detailed quasiclassical trajectory study of integral and differential cross sections, product rovibrational populations, and internal energy distributions has been carried out for the Cl+CH(4) and Cl+CD(4) reactions, and the theoretical results have been compared with the available experimental data. It has been shown that the calculated total reaction cross sections versus collision energy for the Cl+CH(4) and Cl+CD(4) reactions is very sensitive to the barrier height. Besides, due to the zero-point energy (ZPE) leakage of the CH(4) molecule to the reaction coordinate in the quasiclassical trajectory (QCT) calculations, the reaction threshold falls below the barrier height of the PES. The ZPE leakage leads to CH(3) and HCl coproducts with internal energy below its corresponding ZPEs. We have shown that a Gaussian binning (GB) analysis of the trajectories yields excitation functions in somehow better agreement with the experimental determinations. The HCl(v'=0) and DCl(v'=0) rotational distributions are as well very sensitive to the ZPE problem. The GB correction narrows and shifts the rotational distributions to lower values of the rotational quantum numbers. However, the present QCT rotational distributions are still hotter than the experimental distributions. In both reactions the angular distributions shift from backward peaked to sideways peaked as collision energy increases, as seen in the experiments and other theoretical calculations.     

New exchange-Coulomb N2-Ar potential-energy surface and its comparison with other recent N2-Ar potential-energy surfaces

The reliability of five N2-Ar potential-energy surfaces in representing the N2-Ar interaction has been investigated by comparing their abilities to reproduce a variety of experimental results, including interaction second viral coefficients, bulk transport properties, relaxation phenomena, differential scattering cross sections, and the microwave and infrared spectra of the van der Waals complexes. Four of the surfaces are the result of high-level ab initio quantal calculations; one of them utilized fine tuning by fitting to microwave data. To date, these four potential-energy surfaces have only been tested against experimental microwave data. The fifth potential-energy surface, based upon the exchange-Coulomb potential-energy model for the interaction of closed-shell species, is developed herein: it is a combination of a damped dispersion energy series and ab initio calculations of the Heitler-London interaction energy, and has adjustable parameters determined by requiring essentially simultaneous agreement with selected quality interaction second viral coefficient and microwave data. Comparisons are also made with the predictions of three other very good literature potential-energy surfaces, including the precursor of the new exchange-Coulomb potential-energy surface developed here. Based upon an analysis of a large body of information, the new exchange-Coulomb and microwave-tuned ab initio potential-energy surfaces provide the best representations of the N2-Ar interaction; nevertheless, the other potential-energy surfaces examined still have considerable merit with respect to the prediction of specific properties of the N2-Ar van der Waals complex. Of the two recommended surfaces, the new exchange-Coulomb surface is preferred on balance due to its superior predictions of the effective cross sections related to various relaxation phenomena, and to its reliable, and relatively simple, representation of the long-range part of the potential-energy surface. Moreover, the flexibility still inherent in the exchange-Coulomb potential form can be further exploited, if required, in future studies of the N2-Ar interaction.     

Quantum scattering calculation for reaction Br + H2 on two potential energy surfaces

Three‐dimensional time‐dependent quantum wave packet calculations have been carried out for Br + H₂ on a new global ab initio and a semi‐empirical extended London–Eyring–Polanyi–Sato potential energy surface. It is shown that on the ab initio surface, the threshold energy is much lower, and the reaction probabilities, cross sections, and rate constants are much larger. The effects of the initial rovibrational excitation have also been studied. Comparison of rate constants with experimental measurement implies that the ab initio surface is more suitable for quantum dynamic calculation. The possible reasons and mechanism for the dynamical difference on the two PES are analyzed and discussed. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Quantum dynamics of inelastic excitations and charge transfer processes in the H+ +NO collisions

State-resolved differential cross section, integral cross section, average vibrational energy transfer, and the relative transition probability are computed for the H(+)+NO system using our newly obtained ab initio potential energy surfaces (PES) at the multireference configuration interaction level of accuracy employing the correlation consistent polarized valence triple zeta basis set. The quantum dynamics is treated within the vibrational close-coupling rotational infinite-order sudden approximation using the coupled ground state and first excited state ab initio quasidiabatic PES. The computed collision attributes for the inelastic vibrational excitation are compared with the state-to-state scattering data available at E(c.m.)=9.5 eV and E(c.m.)=29.03 eV and are found to be in overall good agreement with those of the experiments. The results for the vibrational charge transfer processes at these collision energies are also presented.     

Stark coefficients for highly excited rovibrational states of H2O

Quantum beat spectroscopy is combined with triple-resonance vibrational overtone excitation to measure the Stark coefficients (SCs) of the water molecule for 28 rovibrational levels lying from 27,600 to 41,000 cm(-1). These data provide a stringent test for assessing the accuracy of the available potential energy surfaces (PESs) and dipole moment surfaces (DMSs) of this benchmark molecule in this energy region, which is inaccessible by direct absorption. SCs, calculated using the combination of a high accuracy, spectroscopically determined PES and a recent ab initio DMS, are within the 1% accuracy of available experimental data for levels below 25,000 cm(-1), and within 4.5% for coefficients associated with levels up to 35,000 cm(-1). However, the error in the computed coefficients is over 60% for the very high rovibrational states lying just below the lowest dissociation threshold, due, it seems, to lack of a high accuracy PES in this region. The comparative analysis suggests further steps, which may bring the theoretical predictions closer to the experimental accuracy.     

A new "spectroscopic" potential energy surface for formaldehyde in its ground electronic state

We report a new "spectroscopic" potential energy surface (PES) of formaldehyde (H(2)(12)C(16)O) in its ground electronic state, obtained by refining an ab initio PES in a least-squares fitting to the experimental spectroscopic data for formaldehyde currently available in the literature. The ab initio PES was computed using the CCSD(T)/aug-cc-pVQZ method at 30 840 geometries that cover the energy range up to 44 000 cm(-1) above equilibrium. Ro-vibrational energies of formaldehyde were determined variationally for this ab initio PES by means of the program TROVE [Theoretical ROtation-Vibration Energies; S. N. Yurchenko, W. Thiel, and P. Jensen, J. Mol. Spectrosc. 245, 126 (2007)]. The parameter values in the analytical representation of the PES were optimized in fittings to 319 ro-vibrational energies with J = 0, 1, 2, and 5. The initial parameter values in the fittings were those of the ab initio PES, the ro-vibrational eigenfunctions obtained from this PES served as a basis set during the fitting process, and constraints were imposed to ensure that the refined PES does not deviate unphysically from the ab initio one in regions of configuration space not sampled by the experimental data. The resulting refined PES, referred to as H(2)CO-2011, reproduces the available experimental J ≤ 5 data with a root-mean-square error of 0.04 cm(-1).     

Molecular-beam study of the water-helium system: features of the isotropic component of the intermolecular interaction and a critical test for the available potential-energy surfaces

We report molecular-beam measurements of the total integral cross sections for the scattering of water molecules by helium atoms. A combined analysis of the new experimental data together with available differential cross section results has allowed an accurate determination of the isotropic component of the interaction potential for this prototypical system. The potential well shows a depth of 0.265 +/- 0.010 kJ/mol at a distance between He and the center of mass of the water molecule of 0.345 +/- 0.02 nm. An effective isotropic long-range attraction constant C(LR) = (6.3+/-0.3) x 10(-4) kJ mol(-1) nm(-6), including both dispersion and induction contributions, has also been determined. The most recent and accurate ab initio potential-energy surfaces have been tested against these new experimental results.     

Potential energy surface of HDO up to 25,000 cm-1

A new spectroscopically determined potential energy surface (PES) for HD(16)O is presented. This surface is constructed by adjusting the high accuracy ab initio PES of Polyansky et al. [Science 299, 539 (2003)] by fitting to both published experimental data and our still unpublished data. This refinement used experimentally derived term values up to 25,000 cm(-1) and with J< or =8: a data set of 3478 energy levels once some levels with ambiguous assignment is excluded. To improve the extrapolation properties of the empirical PES, the restraint that the resulting PESs remain close to the ab initio surface was imposed. The new HDO_07 PES reproduces the experimental data, including high J levels not included in the fit, with a root mean square error of 0.035 cm(-1). Predictions for rotation-vibration term values up to J=12 are made.     

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

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

Accuracy of recent potential energy surfaces for the He-N2 interaction. I. Virial and bulk transport coefficients

A new exchange-Coulomb semiempirical model potential energy surface for the He-N2 interaction has been developed. Together with two recent high-level ab initio potential energy surfaces, it has been tested for the reliability of its predictions of second-virial coefficients and bulk transport phenomena in binary mixtures of He and N2. The agreement with the relevant available measurements is generally within experimental uncertainty for the exchange-Coulomb surface and the ab initio surface of Patel et al. [J. Chem. Phys. 119, 909 (2003)], but with slightly poorer agreement for the earlier ab initio surface of Hu and Thakkar [J. Chem. Phys. 104, 2541 (1996)].     

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

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

Ab initio quantum dynamical study of the vinylidene-acetylene isomerization

The dynamics of vinylidene-acetylene isomerization has been studied quantum mechanically, based on a new ab initio calculated potential energy surface (PES). The grid underlying the PES and the dynamic treatment rest on planar 4-atom Jacobi-like coordinates. Three degrees of freedom, the C-C bond length and the two Jacobi-like angles are treated explicitly, while the two other coordinates relating to the C-H bond lengths have been optimized in the ab initio treatment. The energies of the resulting 3D grid have been computed with the CCSD(T) method, using a cc-pVTZ basis set, while selected stationary points have been characterized at a higher level of theory. The subsequent dynamic treatment reproduces very satisfactorily the experimental photodetachment spectrum of Lineberger and coworkers. Preliminary results are also reported for the lifetime of vinylidene. Received: 5 June 1998/Accepted: 23 July 1998 / Published online: 9 October 1998     

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.     

Ultracold O2 + O2 collisions in a magnetic field: on the role of the potential energy surface

The collision dynamics of (17)O(2)((3)Σ(g)(-)) + (17)O(2)((3)Σ(g)(-)) in the presence of a magnetic field is studied within the close-coupling formalism in the range between 10 nK and 50 mK. A recent global ab initio potential energy surface (PES) is employed and its effect on the dynamics is analyzed and compared with previous calculations where an experimentally derived PES was used [T. V. Tscherbul et al., New J. Phys 134, 055021 (2009)]. Compared to the results using the older PES, magnetic-field dependence of the low-field-seeking state in the ultracold regime is characterized by a very large background scattering length, a(bg), and cross sections exhibit broader and more pronounced Feshbach resonances. The marked resonance structure is somewhat surprising considering the influence of inelastic scattering but it can be explained by resorting to the analytical van der Waals theory, where the short-range amplitude of the entrance channel wavefunction is enhanced by the large a(bg). This strong sensitivity to the short range of the ab initio PES persists up to relatively high energies (10 mK). After this study and despite quantitative predictions are very difficult, it can be concluded that the ratio between elastic and spin relaxation scattering is generally small, except for magnetic fields which are either low or close to an asymmetric Fano-type resonance. Some general trends found here, such as a large density of quasibound states and a propensity toward large scattering lengths, could be also characteristic of other anisotropic molecule-molecule systems.     

An Ab initio Theoretical Study on the Nonadiabatic Coupling for Na＋I2 Collision System

The ionic and neutral state potential energy surfaces (PESs) of Na I2 collision system have been calculated on QCISD(T) level by using ab initio method.The location and depth of the potential well,the collision radius and their fine structures have been analyzed at the equilibrium geometry of I2 molecule.The electronic transfer probabilities are also calculated in terms of Landau-Zener model.The lifetime of scattering resonance state is evaluated by the uncertainty principle.All the results have been compared with those obtained according to the Aten-Lanting-Los PES and Feng‘s PES.     

Vibrational and rotational cooling of NO+ in collisions with He

A quantum mechanical investigation of the vibrational and rotational deactivation of NO(+) in collisions with He atoms in the cold and ultracold regime is presented. Ab initio potential energy calculations are carried out at BCCD(T) level and a new global 3D potential energy surface (PES) is obtained by fitting ab initio points within the reproducing kernel Hilbert space method. As a first test of this PES the bound state energies of the (3)He-NO(+) and (4)He-NO(+) complexes are calculated and compared to previous rigid rotor calculations. The efficiency of the vibrational and the rotational cooling of this molecular ion using a buffer gas of helium is then investigated by performing close coupling scattering calculations for collision energy ranging from 10(-6) to 2000 cm(-1). The calculations are performed for the two isotopes (3)He and (4)He and the results are compared to the available experimental data.