Three-dimensional potential energy surface of the Ar-OH(2Pi i) complex

Pure rotational transitions in the ground state for Ar-OH and Ar-OD [Y. Ohshima et al., J. Chem. Phys. 95, 7001 (1991) and Y. Endo et al., Faraday Discuss. 97, 341 (1994)], those in the excited states of the OH vibration, nu(s)=1 and 2, observed by Fourier-transform microwave spectroscopy in the present study, rotation-vibration transitions observed by infrared-ultraviolet double-resonance spectroscopy [K. M. Beck et al., Chem. Phys. Lett. 162, 203 (1989) and R. T. Bonn et al., J. Chem. Phys. 112, 4942 (2000)], and the P-level structure observed by stimulated emission pumping spectroscopy [M. T. Berry et al., Chem. Phys. Lett. 178, 301 (1991)] have been simultaneously analyzed to determine the potential energy surface of Ar-OH in the ground state. A Schrodinger equation, considering all the freedom of motions for an atom-diatom system in the Jacobi coordinate, R, theta, and r, was numerically solved to obtain energies of the rovibrational energy levels using the discrete variable representation method. A three-dimensional potential energy surface is determined by a least-squares fitting. In the analysis the potential parameters, obtained by ab initio calculations at the RCCSD(T) level of theory with a set of basis functions of aug-cc-pVTZ and midbond functions, are used as initial values. The determined intermolecular potential energy surface and its dependence on the OH monomer bond length are compared with those of an isovalent radical complex, Ar-SH.     

Intermolecular potential energy surface of Ar-NO

Rotational spectra of an open-shell complex, Ar-NO, in the electronic ground state have been analyzed by employing an analysis using a free-rotor model, where previously observed data by Mills et al. [J. Phys. Chem. 90, 3331 (1986); 90, 4961 (1986)] and additional transitions observed by Fourier-transform microwave spectroscopy in the present study are simultaneously analyzed with a standard deviation of the least-squares fit to be 27.5 kHz. A two-dimensional intermolecular potential energy surface for Ar-NO has been determined from the analysis. The determined potential energy surface is compared with those of Ar-OH and Ar-SH, which are also complexes containing an open-shell species with the 2Pi ground electronic state.     

Exploring the new three-dimensional ab initio interaction energy surface of the Ar-HF complex: rovibrational calculations for Ar-HF and Ar-DF with vibrationally excited diatoms

Several features and the performance of the recently published [P. Jankowski and M. Ziolkowski, Mol. Phys. 104, 2293 (2006)] three-dimensional intermolecular potential energy surface for the Ar-HF complex have been investigated. This full-dimensional surface has been obtained using the method of the local expansion of the exact interaction energy surface [P. Jankowski, J. Chem. Phys. 121, 1655 (2004)] in the Taylor series with respect to intramolecular coordinates. The interaction energies have been calculated with the coupled-cluster supermolecular method with single, double, and noniterative triple excitations. The convergence of the interaction energy with respect to the size of the basis set is discussed. The two-dimensional surfaces resulting from averaging of the full-dimensional surface over the intramolecular vibration of HF have been obtained and directly compared to the empirical H6(4,3,2) set of surfaces proposed by Hutson [J. Chem. Phys. 96, 6752 (1992)]. A very good agreement has been observed. The averaged potentials have been used to calculate the rovibrational energy levels of the Ar-HF and Ar-DF complexes and compared to the experimental data. The accuracy of rovibrational calculations achieved with the new surface is much better than with any of the ab initio surfaces available so far. Predictions of the rovibrational energy levels and spectroscopic constants have also been done for Ar-HF with HF in the v=4,5 vibrational states, and for Ar-DF and DF in the v=3,4 states. The full-dimensional surface studied in this paper is the first ab initio surface which is fully compatible with the empirical H6(4,3,2) surface proposed by Hutson.     

Rotational excitation cross sections of para-H2 + para-H2 collisions. A full-dimensional wave-packet propagation study using an exact form of the kinetic energy

A full-dimensional quantum dynamical study of the rotational excitation in para-para H2 + H2 collisions using the potential-energy surface of Boothroyd et al. [J. Chem. Phys. 116, 666 (2002)] is reported. The multiconfiguration time-dependent Hartree algorithm is utilized to propagate wave packets and the cross sections for collision energies up to 1.2 eV are determined by a flux analysis through the interaction of the wave packet with a complex absorbing potential. Calculations for a collection of total angular momenta up to J = 70 are performed; the missing channels are obtained with a J-interpolation algorithm.     

Close-coupling study of rotational energy transfer of CO (upsilon=2) by collisions with He atoms

Quantum close-coupling scattering calculations of rotational energy transfer in the vibrationally excited CO due to collisions with He atom are presented for collision energies between 10(-5) and approximately 1000 cm-1 with CO being initially in the vibrational level upsilon=2 and rotational levels j=0,1,4, and 6. The He-CO interaction potential of Heijmen et al. [J. Chem. Phys. 107, 9921 (1997)] was adopted for the calculations. Cross sections for rovibrational transitions and state-to-state rotational energy transfer from selected initial rotational levels were computed and compared with recent measurements of Carty et al. [J. Chem. Phys. 121, 4671 (2004)] and available theoretical results. Comparison in all cases is found to be excellent, providing a stringent test for the scattering calculations as well as the reliability of the He-CO interaction potential by Heijmen et al.     

Characterization of silicon-carbon clusters by infrared laser spectroscopy: the nu 3(sigma u) band of linear Si2C3

The nu 3(sigma u) fundamental vibration of 1 sigma g+ Si2C3 has been observed using a laser vaporization-supersonic cluster beam-diode laser spectrometer. Forty rovibrational transitions were measured in the range of 1965.8 to 1970.9 cm-1 with a rotational temperature of 10-15 K. A least-squares fit of these transitions yielded the following molecular constants: nu 3(sigma u)=1968.188 31(18) cm-1, B"=0.031 575 1(60) cm-1, and B'=0.031 437 4(57) cm-1. These results are in excellent agreement with recent Fourier transform infrared (FTIR) measurements of Si2C3 trapped in a solid Ar matrix [J. Chem. Phys. 100, 181(1994)] and with ab initio calculations [J. Chem. Phys. 100, 175 (1994)] which suggest cumulenic-like bonding for Si2C3, analogous to the isovalent C5 carbon cluster.     

Reactivity enhancement of ultracold O(3P)+H2 collisions by van der Waals interactions

The role of van der Waals forces in O((3)P)+H(2)(upsilon=1,j=0) collisions is investigated theoretically at low and ultralow temperatures. Quantum scattering calculations have been performed for zero total angular momentum using the lowest London-Eyring-Polanyi-Sato double-polynomial (3)A(") potential-energy surface reported by [Rogers et al., J. Phys. Chem. A 104, 2308 (2000)] and its recent BMS1 and BMS2 extensions developed by [Brandao et al., J. Chem. Phys. 121, 8861 (2004)] which provide a more accurate treatment of the van der Waals interaction. Our calculations show that van der Waals forces strongly influence chemical reactivity at ultracold translational energies. The presence of a zero-energy resonance for the BMS1 surface is found to enhance reactivity in the ultracold regime and shift the Wigner threshold to lower temperatures.     

Rovibrational spectra of LiH2+, LiHD+ and LiD2+ determined from FCI property surfaces

Full configuration interaction (FCI) has been used in conjunction with the lithium [6s5p3d1f] (Iron, M. A.; et al. Mol. Phys. 2004, 101, 1345) and hydrogen aug-cc-pVTZ basis sets to construct an 83-point potential energy surface of the 1A1 ground state of 7LiH2+. Vibrational and rovibrational wave functions of the (6,7)LiH2+, (6,7)LiHD+, and (6,7)LiD2+ ground states were calculated variationally using an Eckart-Watson Hamiltonian. For (7)LiD2+, rovibrational transition frequencies for K = 0, 1, 2 and J < or = 10 are within ca. 0.1% of recent experimental values (Thompson, C. D.; et al. J. Chem. Phys. 2006, 125, 044310). A 47-point FCI dipole moment surface was embedded in the rovibrational Hamiltonian to calculate vibrational and rovibrational radiative properties. At 296 K, with v < or = 4 and J < or = 4, the 2(02) <-- 3(03) rotational transition in the |001> band was found to have the greatest spectral intensity with respect to the ground electronic states of (6,7)LiH2+, (6,7)LiHD+, and (6,7)LiD2+. In each case, the most intense rovibrational transitions have been assigned unequivocally using the J, Ka, Kc assignment scheme.     

A calculation of the rovibronic energies and spectrum of the B1A1 electronic state of SiH2

The B1A1 electronic state of silylene (SiH2) is the second excited singlet state of the molecule and, like the analogous c state of methylene (CH2), it is quasilinear with symmetry 1sigmag+ at linearity. This state dissociates to Si(1D) + H2(1sigmag+). At equilibrium, the B state of SiH2 has an energy that we calculate to be 0.71 eV above that of the dissociation products. However, there is a barrier to dissociation that allows quasibound rovibrational levels to occur, and some have been observed recently [Y. Muramoto et al., J. Chem. Phys. 122, 154302 (2005)]. Starting with our analytical ab initio potential-energy surface, we adjusted it in a fitting to the experimental term values in order to determine the optimum potential-energy function in the bound region. This potential has a C2v equilibrium structure with a SiH bond length of 1.459 angstroms and a bond angle of 165.4 degrees; the barrier to linearity is only 129 cm(-1). Using the optimized potential-energy surface we calculate B-state term values, and using our calculated y and z dipole moment surfaces, we simulate the rotation-vibration spectrum of the state in order to assist in the detection of the matrix isolation spectrum.     

Supersonically cooled hydronium ions in a slit-jet discharge: high-resolution infrared spectroscopy and tunneling dynamics of HD2O+

Jet-cooled high-resolution infrared spectra of partially deuterated hydronium ion (HD2O+) in the O-H stretch region (nu3 band) are obtained for the first time, exploiting the high ion densities, long absorption path lengths, and concentration modulation capabilities of the slit-jet discharge spectrometer. Least-squares analysis with a Watson asymmetric top Hamiltonian yields rovibrational constants and provides high level tests of ab initio molecular structure predictions. Transitions out of both the lower (nu3(+)<--0(+)) and the upper (nu3(-)<--0(-)) tunneling levels, as well as transitions across the tunneling gap (nu3(-)<--0(+)) are observed. The nu3(-)<--0(+) transitions in HD2O+ acquire oscillator strength by loss of D(3h) symmetry, and permit both ground-state-[27.0318(72) cm(-1)] and excited-state-[17.7612(54) cm(-1)]-tunneling splittings to be determined to spectroscopic precision from a single rovibrational band. The splittings and band origins calculated with recent high level ab initio six-dimensional potential surface predictions for H3O+ and isotopomers [X. C. Huang, S. Carter, and J. M. Bowman, J. Chem. Phys. 118, 5431 (2003); T. Rajamaki, A. Miani, and L. Halonen, J. Chem. Phys. 118, 10929 (2003)] are in very good agreement with the current experimental results.     

Rotational spectrum, potential energy surface, and bound states of the weakly bound complex He-N2O

Pure rotational transitions of the weakly bound complex He-N(2)O and three minor isotopomers (He-(14)N(15)NO, He-(15)N(14)NO, and He-(15)N(15)NO) were measured in the frequency region from 6 to 20 GHz. Predictions for the microwave transition frequencies were based on the infrared work by Tang and McKellar [J. Chem. Phys. 117, 2586 (2002)]. In the case of (14)N containing isotopomers, nuclear quadrupole hyperfine structure of the rotational transitions was observed and analyzed. The resulting spectroscopic parameters were used to determine geometrical and dynamical information about the complex. An ab initio potential energy surface was calculated at the coupled cluster level of theory with single and double excitations and perturbative inclusion of triple excitations. This surface was constructed using the augmented correlation consistent polarized valence triple zeta basis set for all atoms with the inclusion of bond functions for the van der Waals bond. Bound state calculations were done to determine the energies of low-lying rovibrational levels that are supported by the potential energy surface. The resulting transition energies agree with the experimental values to 1% or better.     

Ab initio potential energy surfaces, bound states, and electronic spectrum of the Ar-SH complex

New ab initio potential energy surfaces for the (2)Pi ground electronic state of the Ar-SH complex are presented, calculated at the RCCSD(T)/aug-cc-pV5Z level. Weakly bound rotation-vibration levels are calculated using coupled-channel methods that properly account for the coupling between the two electronic states. The resulting wave functions are analyzed and a new adiabatic approximation including spin-orbit coupling is proposed. The ground-state wave functions are combined with those obtained for the excited (2)Sigma(+) state [D. M. Hirst, R. J. Doyle, and S. R. Mackenzie, Phys. Chem. Chem. Phys. 6, 5463 (2004)] to produce transition dipole moments. Modeling the transition intensities as a combination of these dipole moments and calculated lifetime values [A. B. McCoy, J. Chem. Phys. 109, 170 (1998)] leads to a good representation of the experimental fluorescence excitation spectrum [M.-C. Yang, A. P. Salzberg, B.-C. Chang, C. C. Carter, and T. A. Miller, J. Chem. Phys. 98, 4301 (1993)].     

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.     

A new potential energy surface and predicted infrared spectra of He-CO2: dependence on the antisymmetric stretch of CO2

A new potential energy surface involving the antisymmetric Q(3) normal mode of CO(2) for the He-CO(2) van der Waals complex is constructed at the coupled-cluster singles and doubles with noniterative inclusion of connected triple [CCSD(T)] level with augmented correlation-consistent quadruple-zeta (aug-cc-pVQZ) basis set plus bond functions. Two vibrationally adiabatic potentials with CO(2) at both the ground and the first excited vibrational states are generated from the integration of the three-dimensional potential over the Q(3) coordinate. The potential has a T-shaped global minimum and two equivalent linear local minima. The bound rovibrational energy levels are obtained using the radial discrete variable representation/angular finite basis representation method and the Lanczos algorithm. The observed band origin shift of the complex (0.0946 cm(-1)) is successfully reproduced by our calculation (0.1034 cm(-1)). The infrared spectra of the complex are also predicted. The fundamental band is in excellent agreement with the experiment. Most of the transitions corresponding to the observed hot band [M. J. Weida et al., J. Chem. Phys. 101, 8351 (1994)] are assigned reasonably.     

State-to-state rotational transitions in H2+H2 collisions at low temperatures

We present quantum mechanical close-coupling calculations of collisions between two hydrogen molecules over a wide range of energies, extending from the ultracold limit to the superthermal region. The two most recently published potential energy surfaces for the H(2)-H(2) complex, the so-called Diep-Johnson (DJ) [J. Chem. Phys. 112, 4465 (2000); 113, 3480 (2000)] and Boothroyd-Martin-Keogh-Peterson (BMKP) [J. Chem. Phys. 116, 666 (2002)] surfaces, are quantitatively evaluated and compared through the investigation of rotational transitions in H(2)+H(2) collisions within rigid rotor approximation. The BMKP surface is expected to be an improvement, approaching chemical accuracy, over all conformations of the potential energy surface compared to previous calculations of H(2)-H(2) interaction. We found significant differences in rotational excitation/deexcitation cross sections computed on the two surfaces in collisions between two para-H(2) molecules. The discrepancy persists over a large range of energies from the ultracold regime to thermal energies and occurs for several low-lying initial rotational levels. Good agreement is found with experiment B. Mate et al., [J. Chem. Phys. 122, 064313 (2005)] for the lowest rotational excitation process, but only with the use of the DJ potential. Rate coefficients computed with the BMKP potential are an order of magnitude smaller.     

Stabilization and rovibronic spectra of the T-shaped and linear ground-state conformers of a weakly bound rare-gas-homonuclear dihalogen complex: He...Br2

Laser-induced fluorescence spectra of Br(2) entrained in a He supersonic expansion have been recorded in the Br(2) B-X, 8-0, 12-0, and 21-0 spectral regions at varying downstream distances, and thus different temperature regimes. Features associated with transitions of the T-shaped and linear He...Br(2)(X,nu(") = 0) complexes are identified. The changes in the relative intensities of the T-shaped and linear features with cooling in the expansion indicate that the linear conformer is energetically more stable than the T-shaped conformer. A He + Br(2)(X,nu(") = 0) ab initio potential-energy surface, computed at the coupled cluster level of theory with a large, flexible basis set, is used to calculate the binding energies of the two conformers, 15.8 and 16.5 cm(-1) for the T-shaped and linear complexes, respectively. This potential and an excited-state potential [M. P. de Lara-Castells, A. A. Buchachenko, G. Delgado-Barrio, and P. Villareal, J. Chem. Phys. 120, 2182 (2004)] are used to calculate the excitation spectra of He...(79)Br(2)(X,nu(") = 0) in the Br(2) B-X, 12-0 region. The calculated spectra are used to make spectral assignments and to determine the energies of the excited-state intermolecular vibrational levels accessed in the observed transitions. Temperature-dependent laser-induced fluorescence spectra and a simple thermodynamic model [D. S. Boucher, J. P. Darr, M. D. Bradke, R. A. Loomis, and A. B. McCoy, Phys. Chem. Chem. Phys. 6, 5275 (2004)] are used to estimate that the linear conformer is 0.4(2) cm(-1) more strongly bound than the T-shaped conformer. Two-laser action spectroscopy experiments reveal that the binding energy of the linear He...(79)Br(2)(X,nu(") = 0) conformer is 17.0(8) cm(-1), and that of the T-shaped He...(79)Br(2)(X,nu(") = 0) conformer is then 16.6(8) cm(-1), in good agreement with the calculated values.     

Rovibrational energy levels of hydrogen peroxide, studied by MULTIMODE with a reaction path Hamiltonian

Recently, Carter and Handy [J. Chem. Phys. 113 (2000) 987] have introduced the theory of the reaction path Hamiltonian (RPH) [J. Chem. Phys. 72 (1980) 99] into the variational scheme MULTIMODE, for the calculation of the J=0 vibrational levels of polyatomic molecules, which have a single large-amplitude motion. In this theory the reaction path coordinate s becomes the fourth dimension of the moment-of-inertia tensor, and must be treated separately from the remaining 3N-7 normal coordinates in the MULTIMODE program. In the modified program, complete integration is performed over s, and the M-mode MULTIMODE coupling approximation for the evaluation of the matrix elements applies only to the 3N-7 normal coordinates. In this paper the new algorithm is extended to the calculation of rotational-vibration energy levels (i.e. J>0) with the RPH, following from our analogous implementation for rigid molecules [Theoret. Chem. Acc. 100 (1998) 191]. The full theory is given, and all extra terms have been included to give the exact kinetic energy operator. In order to validate the new code, we report studies on hydrogen peroxide (H2O2), where the reaction path is equivalent to torsional motion. H2O2 has previously been studied variationally using a valence coordinate Hamiltonian; complete agreement for calculated rovibrational levels is obtained between the previous results and those from the new code, using the identical potential surface. MULTIMODE is now able to calculate rovibrational levels for polyatomic molecules which have one large-amplitude motion.     

Rovibrational energy transfer in ortho-H2+para-H2 collisions

We present the results of a full-dimensional quantum mechanical study of the rovibrational energy transfer in the collision between ortho-H2 and para-H2 in the energy range of 0.1-1.0 eV. The multiconfiguration time-dependent Hartree algorithm has been used to propagate the wave packets on the global potential energy surface by Boothroyd et al. [J. Chem. Phys. 116, 666 (2002)] and on a modification of this surface where the short range anisotropy is reduced. State-to-state attributes such as probabilities or integral cross sections are obtained using the formalism of Tannor and Weeks [J. Chem. Phys. 98, 3884 (1993)] by Fourier transforming the correlation functions. The effect of initial rotation of the diatoms on the inelastic and de-excitation processes is investigated.     

Accurate quantum calculation of the bound and resonant rovibrational states of Li-(H2)

In a recent paper [B. Poirier, Chem. Phys. 308, 305 (2005)] a full-dimensional quantum method for computing the rovibrational dynamics of triatomic systems was presented, incorporating three key features: (1) exact analytical treatment of Coriolis coupling, (2) three-body "effective potential," and (3) a single bend angle basis for all rotational states. In this paper, these ideas are applied to the Li-(H2) electrostatic complex, to compute all of the rovibrational bound state energies, and a number of resonance energies and widths, to very high accuracy (thousandths of a wave number). This application is very challenging, owing to the long-range nature of the interaction and to narrow level spacings near dissociation. Nevertheless, by combining the present method with a G4 symmetry-adapted phase-space-optimized representation, only modest basis sizes are required for which the matrices are amenable to direct diagonalization. Several new bound levels are reported, as compared with a previous calculation [D. T. Chang, G. Surratt, G. Ristroff, and G. I. Gellene, J. Chem. Phys. 116, 9188 (2002)]. The resonances exhibit a clear-cut separation into shape and Feshbach varieties, with the latter characterized by extremely long lifetimes (microseconds or longer).