Collisionally assisted spectroscopy of water from 27,000 to 34,000 cm(-1)

We report here an experimental approach that enables measurement of weak transitions to a wide range of rovibrational levels of water in the energy region 27,000-34,200 cm(-1). We have previously demonstrated the use of laser double-resonance overtone excitation to access highly excited vibrational levels from single rovibrational states. Although this approach simplifies the assignment of the spectra, it strongly reduces the number of observed transitions and hence our ability to test theoretical predictions. Here, we increase significantly the number of observed transitions by allowing rotational relaxation of H2O at intermediate levels of the double-resonance excitation scheme to the levels of the same nuclear spin (ortho or para). Our recently developed semiempirical potential energy surface PES12 enables assignment of the resulting complex spectra and reproduction of the measured transitions with accuracy better than 1 cm(-1).     

Rovibrational spectra of ammonia. I. Unprecedented accuracy of a potential energy surface used with nonadiabatic corrections

In this work, we build upon our previous work on the theoretical spectroscopy of ammonia, NH(3). Compared to our 2008 study, we include more physics in our rovibrational calculations and more experimental data in the refinement procedure, and these enable us to produce a potential energy surface (PES) of unprecedented accuracy. We call this the HSL-2 PES. The additional physics we include is a second-order correction for the breakdown of the Born-Oppenheimer approximation, and we find it to be critical for improved results. By including experimental data for higher rotational levels in the refinement procedure, we were able to greatly reduce our systematic errors for the rotational dependence of our predictions. These additions together lead to a significantly improved total angular momentum (J) dependence in our computed rovibrational energies. The root-mean-square error between our predictions using the HSL-2 PES and the reliable energy levels from the HITRAN database for J = 0-6 and J = 7∕8 for (14)NH(3) is only 0.015 cm(-1) and 0.020∕0.023 cm(-1), respectively. The root-mean-square errors for the characteristic inversion splittings are approximately 1∕3 smaller than those for energy levels. The root-mean-square error for the 6002 J = 0-8 transition energies is 0.020 cm(-1). Overall, for J = 0-8, the spectroscopic data computed with HSL-2 is roughly an order of magnitude more accurate relative to our previous best ammonia PES (denoted HSL-1). These impressive numbers are eclipsed only by the root-mean-square error between our predictions for purely rotational transition energies of (15)NH(3) and the highly accurate Cologne database (CDMS): 0.00034 cm(-1) (10 MHz), in other words, 2 orders of magnitude smaller. In addition, we identify a deficiency in the (15)NH(3) energy levels determined from a model of the experimental data.     

Approaching the full set of energy levels of water

We report here the measurements of rovibrational levels in the electronic ground state of water molecule at the previously inaccessible energies above 26,000 cm(-1). The use of laser double-resonance overtone excitation extends this limit to 34,200 cm(-1), which corresponds to 83% of the water dissociation energy. We use experimental data to generate a semiempirical potential energy surface that now allows prediction of water levels with sub-cm(-1) accuracy at any energy up to the new limit.     

CVRQD ab initio ground-state adiabatic potential energy surfaces for the water molecule

The high accuracy ab initio adiabatic potential energy surfaces (PESs) of the ground electronic state of the water molecule, determined originally by Polyansky et al. [Science 299, 539 (2003)] and called CVRQD, are extended and carefully characterized and analyzed. The CVRQD potential energy surfaces are obtained from extrapolation to the complete basis set of nearly full configuration interaction valence-only electronic structure computations, augmented by core, relativistic, quantum electrodynamics, and diagonal Born-Oppenheimer corrections. We also report ab initio calculations of several quantities characterizing the CVRQD PESs, including equilibrium and vibrationally averaged (0 K) structures, harmonic and anharmonic force fields, harmonic vibrational frequencies, vibrational fundamentals, and zero-point energies. They can be considered as the best ab initio estimates of these quantities available today. Results of first-principles computations on the rovibrational energy levels of several isotopologues of the water molecule are also presented, based on the CVRQD PESs and the use of variational nuclear motion calculations employing an exact kinetic energy operator given in orthogonal internal coordinates. The variational nuclear motion calculations also include a simplified treatment of nonadiabatic effects. This sophisticated procedure to compute rovibrational energy levels reproduces all the known rovibrational levels of the water isotopologues considered, H(2) (16)O, H(2) (17)O, H(2) (18)O, and D(2) (16)O, to better than 1 cm(-1) on average. Finally, prospects for further improvement of the ground-state adiabatic ab initio PESs of water are discussed.     

Intermolecular potential energy surface and spectra of He-HCl with generalization to other rare gas-hydrogen halide complexes

A two-dimensional (rigid monomer) intermolecular potential energy surface (PES) of the He-HCl complex has been obtained from ab initio calculations utilizing the symmetry-adapted perturbation theory (SAPT) and an spdfg basis set including midbond functions. The bond length in HCl was chosen to be equal to the expectation value in the ground vibrational state of isolated HCl. The rigid-monomer potential should be a very good approximation to the complete (three-dimensional) potential for H-Cl distances corresponding to the lowest vibrational levels of the monomer since the He-HCl interaction energy was found to be only weakly dependent on the HCl bond length in this region, at least as compared to systems such as Ar-HF. The calculated points were fitted using an analytic function with ab initio computed asymptotic coefficients. As expected, the complex is loosely bound, with the dispersion energy providing the majority of the attraction. Our SAPT PES agrees with the semiempirical PES of Willey et al. [J. Chem. Phys. 96, 898 (1992)], in finding that, atypically for rare gas-hydrogen halide complexes including the lighter halide atoms, the global minimum is on the Cl side (with intermonomer separation 3.35 A and depth of 32.8 cm(-1)), rather than on the H side, where there is only a local minimum (3.85 A, 30.8 cm(-1)). The ordering of the minima was confirmed by single-point calculations in larger basis sets and complete basis set extrapolations, and also using higher levels of theory. We show that the opposite findings in the recent calculations of Zhang and Shi [J. Mol. Struct: THEOCHEM 589, 89 (2002)] are due to the lack of midbond functions in their basis set. Despite the closeness in depth of the two linear minima, the existence of a relatively high barrier between them invalidates the assumption of isotropy, a feature of some literature potentials. The trends concerning the locations of minima within the family of rare gas-hydrogen halide complexes are rationalized in terms of the physical components of the intermolecular forces and related to monomer properties. The accuracy of the SAPT PES was tested by performing calculations of rovibrational levels. The transition frequencies obtained were found to be in excellent agreement (to within 0.02 cm(-1)) with the measurements of Lovejoy and Nesbitt [J. Chem. Phys. 93, 5387 (1990)]. The SAPT PES predicts a dissociation energy for the complex of 7.74 cm(-1) which is probably more accurate than the experimental value of 10.1+/-1.2 cm(-1). Our analysis of the ground-state rovibrational wave function shows that the He-HCl configuration is favored over the He-ClH configuration despite the ordering of minima. This is due to the greater volume of the well in the former case. We have also determined positions and widths of three low-lying resonance states through scattering calculations. These predictions are expected to be more accurate than values derived from experiment.     

Calibration-quality adiabatic potential energy surfaces for H3(+) and its isotopologues

Calibration-quality ab initio adiabatic potential energy surfaces (PES) have been determined for all isotopologues of the molecular ion H(3)(+). The underlying Born-Oppenheimer electronic structure computations used optimized explicitly correlated shifted Gaussian functions. The surfaces include diagonal Born-Oppenheimer corrections computed from the accurate electronic wave functions. A fit to the 41,655 ab initio points is presented which gives a standard deviation better than 0.1 cm(-1) when restricted to the points up to 6000 cm(-1) above the first dissociation asymptote. Nuclear motion calculations utilizing this PES, called GLH3P, and an exact kinetic energy operator given in orthogonal internal coordinates are presented. The ro-vibrational transition frequencies for H(3)(+), H(2)D(+), and HD(2)(+) are compared with high resolution measurements. The most sophisticated and complete procedure employed to compute ro-vibrational energy levels, which makes explicit allowance for the inclusion of non-adiabatic effects, reproduces all the known ro-vibrational levels of the H(3)(+) isotopologues considered to better than 0.2 cm(-1). This represents a significant (order-of-magnitude) improvement compared to previous studies of transitions in the visible. Careful treatment of linear geometries is important for high frequency transitions and leads to new assignments for some of the previously observed lines. Prospects for further investigations of non-adiabatic effects in the H(3)(+) isotopologues are discussed. In short, the paper presents (a) an extremely accurate global potential energy surface of H(3)(+) resulting from high accuracy ab initio computations and global fit, (b) very accurate nuclear motion calculations of all available experimental line data up to 16,000 cm(-1), and (c) results suggest that we can predict accurately the lines of H(3)(+) towards dissociation and thus facilitate their experimental observation.     

The Al+ -H(2) cation complex: rotationally resolved infrared spectrum, potential energy surface, and rovibrational calculations

The infrared spectrum of the Al(+)-H(2) complex is recorded in the H-H stretch region (4075-4110 cm(-1)) by monitoring Al(+) photofragments. The H-H stretch band is centered at 4095.2 cm(-1), a shift of -66.0 cm(-1) from the Q(1)(0) transition of the free H(2) molecule. Altogether, 47 rovibrational transitions belonging to the parallel K(a)=0-0 and 1-1 subbands were identified and fitted using a Watson A-reduced Hamiltonian, yielding effective spectroscopic constants. The results suggest that Al(+)-H(2) has a T-shaped equilibrium configuration with the Al(+) ion attached to a slightly perturbed H(2) molecule, but that large-amplitude intermolecular vibrational motions significantly influence the rotational constants derived from an asymmetric rotor analysis. The vibrationally averaged intermolecular separation in the ground vibrational state is estimated as 3.03 A, decreasing by 0.03 A when the H(2) subunit is vibrationally excited. A three-dimensional potential energy surface for Al(+)-H(2) is calculated ab initio using the coupled cluster CCSD(T) method and employed for variational calculations of the rovibrational energy levels and wave functions. Effective dissociation energies for Al(+)-H(2)(para) and Al(+)-H(2)(ortho) are predicted, respectively, to be 469.4 and 506.4 cm(-1), in good agreement with previous measurements. The calculations reproduce the experimental H-H stretch frequency to within 3.75 cm(-1), and the calculated B and C rotational constants to within approximately 2%. Agreement between experiment and theory supports both the accuracy of the ab initio potential energy surface and the interpretation of the measured spectrum.     

He-ThO(1Σ+) interactions at low temperatures: elastic and inelastic collisions, transport properties, and complex formation in cold 4He gas

We present an ab initio study of cold (4)He + ThO((1)Σ(+)) collisions based on an accurate potential energy surface (PES) evaluated by the coupled cluster method with single, double, and noniterative triple excitations using an extended basis set augmented by bond functions. Variational calculations of rovibrational energy levels show that the (4)He-ThO van der Waals complex has a binding energy of 10.9 cm(-1) in its ground J = 0 rotational state. The calculated energy levels are used to obtain the temperature dependence of the chemical equilibrium constant for the formation of the He-ThO complex. We find that complex formation is thermodynamically favored at temperatures below 1 K and predict the maximum abundance of free ground-state ThO(v = 0, j = 0) molecules between 2 and 3 K. The calculated cross sections for momentum transfer in elastic He + ThO collisions display a rich resonance structure below 5 cm(-1) and decline monotonically above this collision energy. The cross sections for rotational relaxation accompanied by momentum transfer decline abruptly to zero at low collision energies (<0.1 cm(-1)). We find that Stark relaxation in He + ThO collisions can be enhanced by applying an external dc electric field of less than 100 kV∕cm. Finally, we present calculations of thermally averaged diffusion cross sections for ThO in He gas, and find these to be insensitive to small variations of the PES at temperatures above 1 K.     

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.     

Ground state potential energy curve and dissociation energy of MgH

New high-resolution visible emission spectra of the MgH molecule have been recorded with high signal-to-noise ratios using a Fourier transform spectrometer. Many bands of the A 2Pi-->X 2Sigma+ and B' 2Sigma+-->X 2Sigma+ electronic transitions of 24MgH were analyzed; the new data span the v' = 0-3 levels of the A 2Pi and B'2Sigma+ excited states and the v'=0-11 levels of the X 2Sigma+ ground electronic state. The vibration-rotation energy levels of the perturbed A 2Pi and B' 2Sigma+ states were fitted as individual term values, while those of the X 2Sigma+ ground state were fitted using the direct-potential-fit approach. A new analytic potential energy function that imposes the theoretically correct attractive potential at long-range, and a radial Hamiltonian that includes the spin-rotation interaction were employed, and a significantly improved value for the ground state dissociation energy of MgH was obtained. The v'=11 level of the X 2Sigma+ ground electronic state was found to be the highest bound vibrational level of 24MgH, lying only about 13 cm(-1) below the dissociation asymptote. The equilibrium dissociation energy for the X 2Sigma+ ground state of 24MgH has been determined to be De=11104.7+/-0.5 cm(-1) (1.37681+/-0.00006 eV), whereas the zero-point energy (v'=0) is 739.11+/-0.01 cm(-1). The zero-point dissociation energy is therefore D0=10365.6+/-0.5 cm(-1) (1.28517+/-0.00006 eV). The uncertainty in the new experimental dissociation energy of MgH is more than 2 orders of magnitude smaller than that for the best value available in the literature. MgH is now the only hydride molecule other than H2 itself for which all bound vibrational levels of the ground electronic state are observed experimentally and for which the dissociation energy is determined with subwavenumber accuracy.     

Experimental long range potential of the B 1Pi state in NaRb

The long range potential of the B 1Pi state in NaRb has been investigated by observation of rovibrational levels that it supports, including the high lying ones, with the technique of polarization labeling spectroscopy. This has allowed us to characterize the potential energy curve up to 1.9 cm(-1) from the dissociation limit. The highest observed rovibrational level v=49, J=10 has the outer turning point at R=16.48 A.     

Rovibrational spectra of ammonia. II. Detailed analysis, comparison, and prediction of spectroscopic assignments for 14NH3, 15NH3, and 14ND3

Several aspects of ammonia rovibrational spectra have been investigated using the new HSL-2 potential energy surface that includes an approximate correction for nonadiabatic effects. The unprecedented accuracy of rovibrational energy levels and transition energies computed using HSL-2 was demonstrated in Part I of this study. For (14)NH(3), new assignments for a few ν(3) + ν(4) band transitions and energy levels are suggested, and discrepancies between computed and HITRAN energy levels in the 2ν(4) band are analyzed (2ν(4) is the most difficult band below 5000 cm(-1)). New assignments are suggested for existing or missing 2ν(4) levels. Several new vibrational bands are identified from existing, unassigned HITRAN data, including 2ν(2) + ν(4), (ν(3) + ν(4)) -A(')∕A("), ν(1) + 2ν(2), and 2ν(2) + 2ν(4). The strong mixing between the 2ν(4) and 2ν(2) + ν(4) bands is carefully examined and found to be the source of the difficulties in the experimental modeling of 2ν(4). Discussion is presented for preliminary J = 10 results, where the overall root-mean-square error is estimated to be less than 0.039 cm(-1). The analysis of the 4ν(2) band demonstrates both the reliability and the accuracy of predictions from HSL-2. The full list of computed J = 0 band origins (with assignments) and the inversion splittings up to 7000-8000 cm(-1) above the zero-point energy are presented. J = 0-2 levels are reported for those bands below 5100 cm(-1) that are missing from the HITRAN database. For (15)NH(3), excellent agreement is found for the available ν(2) and ν(3) + ν(4)(E) transition energies, but significant deficiencies are shown for HITRAN levels and several corrections are suggested. The (15)N isotopic effects are presented for the J = 0-6 levels of 13 HITRAN bands. For (14)ND(3), we reproduce the pure rotational inversion spectra line frequencies with an accuracy similar to that for (14)NH(3). However, it is not possible to reproduce simultaneously all four pairs of inversion-split vibrational fundamentals to better than 0.05 cm(-1) uncertainty. It is suggested that a reanalysis of some suspicious (14)ND(3) fundamental bands is required. The analyses presented here and in Part I show that rovibrational energy levels and transition frequencies computed with HSL-2 (with nonadiabatic corrections) remain highly accurate well beyond the experimental data used in the refinement procedure. Calculations using HSL-2 are capable of revealing many deficiencies in experimental analyses of ammonia spectra and provide reliable predictions with similar accuracy. It is expected that the results of this study will be useful in the future interpretation of high-resolution spectra from laboratory experiments or from astronomical observations. The present work represents a very significant advance in the state of our knowledge of the spectroscopy of ammonia and its isotopologues.     

An isotopic-independent highly accurate potential energy surface for CO2 isotopologues and an initial (12)C(16)O2 infrared line list

An isotopic-independent, highly accurate potential energy surface (PES) has been determined for CO(2) by refining a purely ab initio PES with selected, purely experimentally determined rovibrational energy levels. The purely ab initio PES is denoted Ames-0, while the refined PES is denoted Ames-1. Detailed tests are performed to demonstrate the spectroscopic accuracy of the Ames-1 PES. It is shown that Ames-1 yields σ(rms) (root-mean-squares error) = 0.0156 cm(-1) for 6873 J = 0-117 (12)C(16)O(2) experimental energy levels, even though less than 500 (12)C(16)O(2) energy levels were included in the refinement procedure. It is also demonstrated that, without any additional refinement, Ames-1 yields very good agreement for isotopologues. Specifically, for the (12)C(16)O(2) and (13)C(16)O(2) isotopologues, spectroscopic constants G(v) computed from Ames-1 are within ±0.01 and 0.02 cm(-1) of reliable experimentally derived values, while for the (16)O(12)C(18)O, (16)O(12)C(17)O, (16)O(13)C(18)O, (16)O(13)C(17)O, (12)C(18)O(2), (17)O(12)C(18)O, (12)C(17)O(2), (13)C(18)O(2), (13)C(17)O(2), (17)O(13)C(18)O, and (14)C(16)O(2) isotopologues, the differences are between ±0.10 and 0.15 cm(-1). To our knowledge, this is the first time a polyatomic PES has been refined using such high J values, and this has led to new challenges in the refinement procedure. An initial high quality, purely ab initio dipole moment surface (DMS) is constructed and used to generate a 296 K line list. For most bands, experimental IR intensities are well reproduced for (12)C(16)O(2) using Ames-1 and the DMS. For more than 80% of the bands, the experimental intensities are reproduced with σ(rms)(ΔI) < 20% or σ(rms)(ΔI∕δ(obs)) < 5. A few exceptions are analyzed and discussed. Directions for future improvements are discussed, though it is concluded that the current Ames-1 and the DMS should be useful in analyzing and assigning high-resolution laboratory or astronomical spectra.     

The Role of High Excitations in Constructing Sub-spectroscopic Accuracy Intermolecular Potential of He-HCN: Critically Examined by the High-Resolution Spectra with Resonance States

Interpreting high-resolution rovibrational spectra of weakly bound complexes commonly requires spectroscopic accuracy (<1 cm^-1) potential energy surfaces (PES). Constructing high-accuracy ab initio PES relies on the high-level electronic structure approaches and the accurate physical models to represent the potentials. The coupled cluster approaches including single and double excitations with a perturbational estimate of triple excitations (CCSD(T)) have been termed the "gold standard" of electronic structure theory, and widely used in generating intermolecular interaction energies for most van der Waals complexes. However, for HCN-He complex, the observed millimeter-wave spectroscopy with high-excited resonance states has not been assigned and interpreted even on the ab initio PES computed at CCSD(T) level of theory with the complete basis set (CBS) limit. In this work, an effective three-dimensional ab initio PES for HCN-He, which explicitly incorporates dependence on the Q₁ (C-H) normal-mode coordinate of the HCN monomer has been calculated at the CCSD(T)/CBS level. The post-CCSD(T) interaction energy has been examined and included in our PES. Analytic two-dimensional PESs are obtained by least-squares fitting vibrationally averaged interaction energies for v₁(C-H)=0, and 1 to the Morse/Long-Range potential function form with root-mean-square deviations (RMSD) smaller than 0.011 cm^-1. The role and significance of the post-CCSD(T) interaction energy contribution are clearly illustrated by comparison with the predicted rovibrational energy levels. With or without post-CCSD(T) corrections, the value of dissociation limit (D₀) is 8.919 or 9.403 cm^-1, respectively. The predicted millimeter-wave transitions and intensities from the PES with post-CCSD(T) excitation corrections are in good agreement with the available experimental data with RMS discrepancy of 0.072 cm^-1. Moreover, the infrared spectrum for HCN-He complex is predicted for the first time. These results will serve as a good starting point and provide reliable guidance for future infrared studies of HCN doped in (He)_n clusters.     

Close-coupling calculations of low-energy inelastic and elastic processes in 4He collisions with H2: A comparative study of two potential energy surfaces

The two most recently published potential energy surfaces (PESs) for the HeH2 complex, the so-called MR (Muchnick and Russek) and BMP (Boothroyd, Martin, and Peterson) surfaces, are quantitatively evaluated and compared through the investigation of atom-diatom collision processes. The BMP surface is expected to be an improvement, approaching chemical accuracy, over all conformations of the PES compared to that of the MR surface. We found significant differences in inelastic rovibrational cross sections computed on the two surfaces for processes dominated by large changes in target rotational angular momentum. In particular, the H2(nu=1,j=0) total quenching cross section computed using the BMP potential was found to be a factor of 1000 larger than that obtained with the MR surface. A lesser discrepancy persists over a large range of energies from the ultracold to thermal and occurs for other low-lying initial rovibrational levels. The MR surface was used in previous calculations of the H2(nu=1,j=0) quenching rate coefficient and gave results in close agreement with the experimental data of Audibert et al. which were obtained for temperatures between 50 and 300 K. Examination of the rovibronic coupling matrix elements, which are obtained following a Legendre expansion of the PES, suggests that the magnitude of the anisotropy of the BMP potential is too large in the interaction region. However, cross sections for elastic and pure rotational processes obtained from the two PESs differ typically by less than a factor of 2. The small differences may be ascribed to the long-range and anharmonic components of the PESs. Exceptions occur for (nu=10,j=0) and (nu=11,j=1) where significant enhancements have been found for the low-energy quenching and elastic cross sections due to zero-energy resonances in the BMP PES which are not present in the MR potential.     

Kinetic study of the gas-phase reaction of OH with Br2

An experimental, temperature-dependent kinetic study of the gas-phase reaction of the hydroxyl radical with molecular bromine (reaction 1) has been performed by using a pulsed laser photolysis/pulsed-laser-induced fluorescence technique over a wide temperature range of 297-766 K, and at pressures between 6.68 and 40.29 kPa of helium. The experimental rate coefficients for reaction 1 demonstrate no correlation with pressure and exhibit a negative temperature dependence with a slight negative curvature in the Arrhenius plot. A nonlinear least-squares fit with two floating parameters of the temperature-dependent k(1)(T) data set using an equation of the form k(1)(T) = AT(n) yields the recommended expression k(1)(T) = (1.85 x 10(-9))T(-0.66) cm(3) molecule(-1) s(-1) for the temperature dependence of the reaction 1 rate coefficient. The potential energy surface (PES) of reaction 1 was investigated with use of quantum chemistry methods. The reaction proceeds through formation of a weakly bound OH...Br(2) complex and a PES saddle point with an energy below that of the reactants. Temperature dependence of the reaction rate coefficient was modeled by using the RRKM method on the basis of the calculated PES.     

Close-coupling study of rotational energy transfer and differential scattering in H2O collisions with He atoms

Quantum close-coupling scattering calculations of rotational energy transfer (RET) of rotationally excited H(2)O due to collisions with He are presented for collision energies between 10(-6) and 1000 cm(-1) with para-H(2)O initially in levels 1(1,1), 2(0,2), 2(1,1), and 2(2,0) and ortho-H(2)O in levels 1(1,0), 2(1,2), and 2(2,1). Quenching cross sections and rate coefficients for state-to-state RET were computed. Both elastic and inelastic differential cross sections are also calculated and compared with relative experimental results giving generally good agreement in all cases, but less so for inelastic results. Significant differences in the computed collisional parameters, obtained on three different potential energy surfaces (PESs), were found particularly in the ultracold regime. In the thermal regime, the rate coefficients calculated on each of the surfaces are generally in better agreement and comparable, but typically larger, than those obtained in a previous calculation. Unfortunately, a lack of absolute differential or integral inelastic experimental data prevents firm determination of a preferred PES.     

Highly accurate potential-energy and dipole moment surfaces for vibrational state calculations of methane

Full-dimensional ab initio potential-energy surface (PES) and dipole moment surface are constructed for a methane molecule at the CCSD(T)/cc-pVTZ and MP2/cc-pVTZ levels of theory, respectively, by the modified Shepard interpolation method based on the fourth-order Taylor expansion [MSI(4th)]. The reference points for the interpolation have been set in the coupling region of CH symmetric and antisymmetric stretching modes so as to reproduce the vibrational energy levels related to CH stretching vibrations. The vibrational configuration-interaction calculations have been performed to obtain the energy levels and the absorption intensities up to 9000 cm(-1) with the use of MSI(4th)-PES. The calculated fundamental frequencies and low-lying vibrational energy levels show that MSI(4th) is superior to the widely employed quartic force field, giving a better agreement with the experimental values. The absorption bands of overtones as well as combination bands, which are caused by purely anharmonic effects, have been obtained up to 9000 cm(-1). Strongly coupled states with visible intensity have been found in the 6500-9000 cm(-1) region where the experimental data are still lacking.     

Effect of the overall rotation on the cis-trans isomerization of HONO induced by an external field

Rovibrational eigenenergies of HONO are computed and compared to experimental energies available in the literature. For their computation, we use a previously developed potential energy surface (PES) and a newly derived exact kinetic energy operator (KEO) including the overall rotation for a tetra-atomic molecule in non-orthogonal coordinates. In addition, we use the Heidelberg Multi-Configuration Time-Dependent Hartree (MCTDH) package. We compare the experimental rovibrational eigenvalues of HONO available in the literature with those obtained with MCTDH and a previously developed potential energy surface (PES) [F. Richter et al., J. Chem. Phys., 2004, 120, 1306.] for the cis geometry. The effect of the overall rotation on the process studied in our previous work on HONO [F. Richter et al., J. Chem. Phys., 2007, 127, 164315.] leading to the cis→trans isomerization of HONO is investigated. This effect on this process is found to be weak.     

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.