Ab initio calculations and Franck-Condon simulation of the absorption spectra of GeCl2 including anharmonicity.

Geometrical parameters, vibrational frequencies and relative electronic energies of the X1A1, ?3B1 and A1B1 states of GeCl2 have been calculated at the CCSD(T) and/or CASSCF/MRCI level with basis sets of up to aug-cc-pV5Z quality. Core electron correlation and relativistic contributions were also investigated. RCCSD(T)/ aug-cc-pVQZ potential energy functions (PEFs) of the X1A1 and ?3B, states, and a CASSCF/MRCl/aug-cc-pVQZ PEF of the A1B1 state of GeCl2 are reported. Anharmonic vibrational wavefunctions of these electronic states of GeCl2, obtained variationally using the computed PEFs, are employed to calculate the Franck-Condon factors (FCFs) of the ?-X and A-X transitions of GeCl2. Simulated absorption spectra of these transitions based on the computed FCFs are compared with the corresponding experimental laser-induced fluorescence (LIF) spectra of Karolczak et al. [J. Chem. Phys. 1993, 98, 60-70]. Excellent agreement is obtained between the simulated absorption spectrum and observed LIF spectrum of the ?-X transition of GeCl2, which confirms the molecular carrier, the electronic states involved and the vibrational assignments of the LIF spectrum. However, comparison between the simulated absorption spectrum and experimental LIF spectrum of the A-X transition of GeCl2 leads to a revision of vibrational assignments of the LIF spectrum and suggests that the X1A1 state of GeCl2 was prepared in the experimental work, with a non-Boltzmann vibrational population distribution. The X(0,0,1) level is populated over 4000 times more than expected from a Boltzmann distribution at 60 K, which is appropriate for the relative population of the other low-lying vibrational levels, such as the X(1,0,0) and X(0,1,0) levels.     

Ab initio calculations on the X (1)A(') and A (1)A(") states of HPO and Franck-Condon simulation of the single vibronic level emission spectra of HPO and DPO

Minimum-energy geometries and relative electronic energies of the X (1)A(') and A (1)A(") states of HPO have been computed employing the coupled-cluster single-double plus perturbative triple excitations {RCCSD(T)} and/or complete-active-space self-consistent-field (CASSCF) multireference internally contracted configuration interaction (MRCI) methods with basis sets of up to the augmented correlation-consistent polarized-valence quintuple-zeta (aug-cc-pV5Z) quality. In addition, RCCSD(T)/aug-cc-pVQZ and CASSCF/MRCI/aug-cc-pVQZ potential energy functions, anharmonic vibrational wave functions, and energies involving all three vibrational modes for both electronic states of HPO and DPO, and Franck-Condon factors between the two electronic states, which allow for Duschinsky rotation and anharmonicity, were computed. Computed Franck-Condon factors were then used to simulate single vibronic level (SVL) emission spectra recently reported by Tackett and Clouthier [J. Chem. Phys. 117, 10604 (2002)]. Excellent agreement between the simulated and observed spectra was obtained for the A (1)A(")(1,0,0)-->X (1)A(') SVL emission of HPO and DPO, when the best estimated ab initio geometries of the two states, which include contributions from core correlation and extrapolation to the complete basis set limit, were used in the simulation, suggesting that the best estimated ab initio geometry of the A (1)A(") state of HPO, particularly the bond angle of 94.5 degrees , is more reliable than the available experimentally derived geometry. A discussion on the geometrical parameters derived from rotational constants obtained from the rotational analysis of a high-resolution spectrum and from Franck-Condon simulation of the vibrational structure of an electronic spectrum is given.     

Ab initio calculations on SF2 and its low-lying cationic states: anharmonic Franck-Condon simulation of the UV photoelectron spectrum of SF2

Geometry optimization calculations were carried out on the (approximate)X(1)A(1) state of SF2 and the (approximate)X(2)B(1), (approximate)A(2)A(1), (approximate)B(2)B(2), (approximate)C(2)B(2), (approximate)D(2)A(1), and (approximate)E(2)A(2) states of SF2(+) employing the restricted-spin coupled-cluster single-double plus perturbative triple excitation [RCCSD(T)] method and basis sets of up to the augmented correlation-consistent polarized quintuple-zeta [aug-cc-pV(5+d)Z] quality. Effects of core electron (S 2s(2)2p(6) and F 1s(2) electrons) correlation and basis set extension to the complete basis set limit on the computed minimum-energy geometries and relative electronic energies (adiabatic and vertical ionization energies) were investigated. RCCSD(T) potential energy functions (PEFs) were calculated for the (approximate)X(1)A(1) state of SF2 and the low-lying states of SF2(+) listed above employing the aug-cc-pV(5+d)Z and aug-cc-pV5Z basis sets for S and F, respectively. Anharmonic vibrational wave functions of these neutral and cationic states of SF2, and Franck-Condon (FC) factors of the lowest four one-electron allowed neutral photoionizations were computed employing the RCCSD(T) PEFs. Calculated FC factors with allowance for Duschinsky rotation and anharmonicity were used to simulate the first four photoelectron bands of SF2. The agreement between the simulated and observed first bands in the He I photoelectron spectrum reported by de Leeuw et al. [Chem. Phys. 34, 287 (1978)] is excellent. Our calculations largely support assignments made by de Leeuw et al. on the higher ionization energy bands of SF2.     

Franck-Condon simulation of the single vibronic level emission spectra of HSiF and DSiF including anharmonicity

Potential energy functions (PEFs) of the X (1)A(') and A (1)A(") states of HSiF have been computed using the coupled-cluster single-double plus perturbative triple excitations and complete-active-space self-consistent-field multireference internally contracted configuration interaction methods, respectively, employing augmented correlation-consistent polarized-valence quadruple-zeta basis sets. For both electronic states of HSiF and DSiF, anharmonic vibrational wavefunctions and energies of all three modes have been calculated variationally with the ab initio PEFs and using Watson's Hamiltonian for nonlinear molecules. Franck-Condon factors between the two electronic states, allowing for Duschinsky rotation, were computed using the calculated anharmonic vibrational wavefunctions. These Franck-Condon factors were used to simulate the single vibronic level (SVL) emission spectra recently reported by Hostutler et al. in J. Chem. Phys. 114, 10728 (2001). Excellent agreement between the simulated and observed spectra was obtained for the A (1)A(")(1,0,0)-->X (1)A(') SVL emission of HSiF. Discrepancies between the simulated and observed spectra of the A (1)A(")(0,1,0) and (1,1,0) SVL emissions of HSiF have been found. These are most likely, partly due to experimental deficiencies and, partly to inadequacies in the ab initio levels of theory employed in the calculation of the PEFs. Based on the computed Franck-Condon factors, minor revisions of previous vibrational assignments are suggested. The calculated anharmonic wave functions of higher vibrational levels of the X (1)A(') state show strong mixings between the three vibrational modes of HSi stretching, bending, and SiF stretching.     

Ab initio calculations on SCl2 and low-lying cationic states of SCl2+: Franck-Condon simulation of the UV photoelectron spectrum of SCl2

Geometry optimization calculations were carried out on the (approximate)X (1)A(1) state of SCl(2) and the (approximate)X(2)B(1), (approximate)A(2)B(2), (approximate)B(2)A(1), (approximate)C(2)A(1), (approximate)D(2)A(2), and (approximate)E (2)B(2) states of SCl(2) (+) at the restricted-spin coupled-cluster single-double plus perturbative triple excitation [RCCSD(T)] level with basis sets of up to the augmented correlation-consistent polarized quintuple-zeta [aug-cc-pV(5+d)Z] quality. Effects of core electron correlation, basis set extension to the complete basis set limit, and relativistic contributions on computed minimum-energy geometrical parameters and/or relative electronic energies were also investigated. RCCSD(T) potential energy functions (PEFs) were calculated for the (approximate)X (1)A(1) state of SCl(2) and the low-lying states of SCl(2)(+) listed above employing the aug-cc-pV(5+d)Z basis set. Anharmonic vibrational wave functions of these neutral and cationic states of SCl(2), and Franck-Condon (FC) factors of the lowest four one-electron allowed neutral photoionizations were computed employing the RCCSD(T)aug-cc-pV(5+d)Z PEFs. Calculated FC factors with allowance for the Duschinsky rotation and anharmonicity were used to simulate the first four photoelectron (PE) bands of SCl(2). The agreement between simulated and observed He I PE spectra reported by Colton et al. [J. Electron Spectrosc. Relat. Phenom. 3, 345 (1974)] and Solouki et al. [Chem. Phys. Lett. 26, 20 (1974)] is excellent. However, our FC spectral simulations indicate that the first observed vibrational component in the first PE band of SCl(2) is a "hot" band arising from the SCl(2)(+)(approximate)X(2)B(1)(0,0,0)<--SCl(2)(approximate)X (1)A(1)(1,0,0) ionization. Consequently, the experimental adiabatic ionization energy of SCl(2) is revised to 9.55+/-0.01 eV, in excellent agreement with results obtained from state-of-the-art ab initio calculations in this work.     

A combined ab initio and Franck-Condon factor simulation study on the photodetachment spectrum of ScO2(-)

Restricted-spin coupled-cluster single-double plus perturbative triple excitation {RCCSD(T)} potential energy functions (PEFs) of the X(2)B2 state of ScO2 and the 1A1 state of ScO2(-) were computed, employing the augmented correlation-consistent polarized-weighted core-valence quadruple-zeta (aug-cc-pwCVQZ) basis set for Sc and augmented correlation-consistent polarized valence quadruple-zeta (aug-cc-pVQZ) basis set for O, and with the outer core Sc 3s(2)3p(6) electrons being explicitly correlated. Franck-Condon factors, which include allowance for Duschinsky rotation and anharmonicity, were calculated using the computed RCCSD(T) PEFs, and were used to simulate the first photodetachment band of ScO2(-). The simulated spectrum matches well with the corresponding experimental 355 nm photodetachment spectrum of Wu and Wang, J Phys Chem A 1998, 102, 9129, confirming the assignment of the photodetachment spectrum and the reliability of the RCCSD(T) PEFs used. Further calculations on low-lying electronic states of ScO2 gave adiabatic relative electronic energies (T(e)'s) of, and vertical excitation energies (T(v)'s) to, the 2A1, 2B1, and 2A2 states of ScO2 (from the X(2)B2 state of ScO2), as well as electron affinities (EAs) and vertical detachment energies (VDEs) to these neutral states from the 1A1 state of ScO2(-).     

An ab initio study of the low-lying electronic states of YO2 and Franck-Condon simulation of the first photodetachment band of YO2(-)

A variety of density functional theory and ab initio methods, including B3LYP, B98, BP86, CASSCF, CASSCF/RS2, CASSCF/MRCI, BD, BD(T), and CCSD(T), with ECP basis sets of up to the quintuple-zeta quality for Y, have been employed to study the X(2)B2 state of YO2 and the X(1)A1 state of YO2(-). Providing that the Y 4s(2)4p(6) outer-core electrons are included in the correlation treatment, the RCCSD(T) method gives the most consistent results and is concluded to be the most reliable and practical computational method for YO2 and YO2(-). In addition, RCCSD(T) potential energy functions (PEFs) of the X(2)B2 state of YO2 and the X(1)A1 state of YO2(-) were computed, employing the ECP28MDF_aug-cc-pwCVTZ and aug-cc-pVTZ basis sets for Y and O, respectively. Franck-Condon factors, which include allowance for Duschinsky rotation and anharmonicity, were calculated using the computed RCCSD(T) PEFs and were used to simulate the first photodetachment band of YO2(-). The simulated spectrum matches very well with the corresponding experimental 355 nm photodetachment spectrum of Wu, H.; Wang, L.-S. J. Phys. Chem. A 1998, 102, 9129, confirming the reliability of the RCCSD(T) PEFs used. Further calculations on low-lying electronic states of YO2 gave T(e)'s and T(vert)'s of the A(2)A1, B(2)B1, and C(2)A2 states of YO2, as well as EAs and VDEs to these states from the X(1)A1 state of YO2(-). On the basis of the ab initio VDEs obtained in the present study, previous assignments of the second and third photodetachment bands of YO2(-) have been revised.     

A new reaction path for the C + NO reaction: dynamics on the 4A' potential-energy surface

We present a new reaction path without significant barriers for the C + NO reaction, forming ground state N((4)S) and CO. Electronic structure (CASPT2) calculations have been performed for the two lowest (4)A' states of the CNO system. The lowest of these states shows no significant barriers against reaction in the C + NO or O + CN channels. This surface has been fitted to an analytical function using a many-body expansion. Using this surface, and the previously published (2)A' and (2)A' surfaces [Andersson et al., Phys. Chem. Chem. Phys., 2000, 2, 613; Andersson et al., Chem. Phys., 2000, 259, 99], we have performed quasiclassical trajectory (QCT) calculations, obtaining rate coefficients for the C((3)P) + NO(X(2)Pi) --> CO(X(1)Sigma(+)) + N((4)S,(2)D) and C((3)P) + NO(X(2)Pi) --> O((3)P) + CN(X(2)Sigma(+)) reactions. We have also simulated the crossed molecular beam experiments of Naulin et al. [Chem. Phys., 1991, 153, 519] The inclusion of the (4)A' surface in the QCT calculations gives excellent agreement with experiments. This is the first time an adiabatic pathway from C((3)P) + NO(X(2)Pi) to CO(X(1)Sigma(+))+N((4)S) has been reported.     

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.     

Stimulated emission pumping spectroscopy of the [X](1)A' state of CHF

We have recorded stimulated emission pumping (SEP) spectra of the A1A' ' 1A' system of CHF, which reveal rich detail concerning the rovibronic structure of the 1A' state up to approximately 7000 cm-1 above the vibrationless level. Using several intermediate A1A' ' state levels, we obtained rotationally resolved spectra for 16 of the 33 levels observed in our previous single vibronic level (SVL) emission study (Fan et al., J. Chem. Phys. 2005, 123, 014314), in addition to one new level. An anharmonic effective Hamiltonian model poorly reproduces the term energies even with the improved set of data because of the extensive interactions among levels in a given polyad (p) having combinations of nu1, nu2, nu3, which satisfy the relationship p = 2nu1 + nu2 + nu3. However, the precise A rotational constants determined from the SEP data were invaluable in clarifying the assignments for these strongly perturbed levels, and the data are well reproduced using a multiresonance effective Hamiltonian model. The derived vibrational parameters are in good agreement with high level ab initio calculations. The experimental frequencies were combined with those of CDF to derive a harmonic force field and average (rz,r(z)e) structures for the ground state.     

Nonadiabatic quantum dynamics in O(3P)+H2→OH+H: a revisited study

To investigate the extent of nonadiabatic effects in the title reaction, quasi-classical trajectory and nonadiabatic quantum scattering as well as the nonadiabatic quantum-classical trajectory calculations were performed on the accurate ab initio benchmark potential energy surfaces of the lowest (3)A' and (3)A" electronic states [Rogers et al., J Phys Chem A 2000, 104, 2308], together with the spin-orbit coupling matrix [Maiti and Schatz, J Chem Phys 2003, 119, 12360] and the lowest singlet (1) A' potential energy surface [Dobby and Knowles, Faraday Discuss 1998, 110, 247]. Comparison of the calculated total cross sections from both adiabatic and nonadiabatic calculations has demonstrated that for adiabatic channels including (3)A'→(3)A' and (3)A"→(3)A", difference does exist between the two kinds of adiabatic and nonadiabatic calculations, showing nonadiabatic effects to some extent. Such nonadiabatic effects tend to become more conspicuous at high collision energies and are found to be more pronounced with trajectories/quantum wave packet initiated on (3)A' than on (3)A". Furthermore, the present study also showed that nonadiabatic effects can bring the component of forward-scattering in the product angular distributions.     

A combined ab initio/Franck-Condon study of the A-X single-vibronic-level emission spectrum of CCl2 and the photodetachment spectrum of CCl2 -.

State-of-the-art ab initio calculations have been carried out on the X1A1, ?3B1 and A1B1 states of CCl2 and the X2B1 state of CCl2-. Franck-Condon factors including anharmonicity have been calculated, between the CCl2 states, and between the CCl2- X2B1 state and the CCl2 states. They are used to simulate the A-X single-vibronic-level (SVL) emission spectra of CCl2 determined by M.-L. Lui et al. [PCCP 2003, 5, 352] and the 364 nm laser photodetachment spectrum of CCl2- obtained by R. L. Schwartz et al. [J. Phys. Chem. A 1999, 103, 8213]. Comparison between simulated and observed spectra confirms the vibrational assignments of the X2B1 SVL emission spectra and the T0 position of the A1B1 state of CCl2. For the photodetachment spectrum of CCl2-, spectral simulation shows that the higher binding energy ?3B1(CCl2) <-- X2B1(CCl2-) band is well separated from the X1A1(CCl2) <-- X2B1(CCl2-) band. It is concluded that the observed second band, which overlaps heavily with the X1A1(CCl2) <-- X2B1(CCl2-) band in the photodetachment spectrum of CCl2- cannot be assigned to the CCl2(?3B1) + e --> CCl2-(X2B1) detachment process. Further ab initio calculations carried out in the present investigation support the suggestion that the second band in the 364 nm photodetachment spectrum of CCl2- is due to detachment from an excited state of CCl2-, a linear quartet state, to a triplet state of CCl2. These calculations identify the anionic state to be the lowest 4Sigmag- (4Sigma-) state, which photodetaches vertically to the 3Sigmag- (3Sigma-; adiabatically ?3B1) and/or 3Pi(u) (3Pi) states of CCl2 to give the second band observed in the 364 nm photodetachment spectrum of CCl2-.     

The double Renner effect in the X(2)A" and A(2)A' electronic states of HO(2)

A theoretical investigation of the X(2)A" and A(2)A' electronic states of the HO(2) radical is reported. Both electronic states have nonlinear equilibrium geometries and they correlate with a (2)Pi state at linear geometries so that they exhibit the Renner effect. In highly excited bending states, there is tunneling between two equivalent minima (with geometries where the H nucleus is bound to one, or the other, of the two O nuclei), and the two linear geometries H-O-O and O-O-H become accessible to the molecule. Thus, HO(2) affords an example of the so-called double Renner effect. Three-dimensional potential energy surfaces for the X(2)A" and A(2)A' electronic states of HO(2) have been calculated ab initio and the global potential energy surfaces for the states have been constructed. These surfaces have been used, in conjunction with the computer program DR [Odaka et al., J. Mol. Structure 795, 14 (2006); Odaka et al., J. Chem. Phys. 126, 094301 (2007)], for calculating HO(2) rovibronic energies in the "double-Renner"-degenerate electronic states X(2)A" and A(2)A'. The results of the ab initio calculations, the rovibronic energies obtained, and analyses of the wavefunctions for selected states are presented.     

Photodissociation of N2O: triplet states and triplet channel

The role of triplet states in the UV photodissociation of N(2)O is investigated by means of quantum mechanical wave packet calculations. Global potential energy surfaces are calculated for the lowest two (3)A' and the lowest two (3)A' states at the multi-reference configuration interaction level of electronic structure theory using the augmented valence quadruple zeta atomic basis set. Because of extremely small transition dipole moments with the ground electronic state, excitation of the triplet states has only a marginal effect on the far red tail of the absorption cross section. The calculations do not show any hint of an increased absorption around 280 nm as claimed by early experimental studies. The peak observed in several electron energy loss spectra at 5.4 eV is unambiguously attributed to the lowest triplet state 1(3)A'. Excitation of the 2(1)A' state and subsequent transition to the repulsive branch of the 2(3)A' state at intermediate NN-O separations, promoted by spin-orbit coupling, is identified as the main pathway to the N(2)((1)Σ(g)(+))+O((3)P) triplet channel. The yield, determined in two-state wave packet calculations employing calculated spin-orbit matrix elements, is 0.002 as compared to 0.005 ± 0.002 measured by Nishida et al. [J. Phys. Chem. A 108, 2451 (2004)].     

Hydrogen bridging in the compounds X2H (X=Al,Si,P,S)

X2H hydrides (X=Al, Si, P, and S) have been investigated using coupled cluster theory with single, double, and triple excitations, the latter incorporated as a perturbative correction [CCSD(T)]. These were performed utilizing a series of correlation-consistent basis sets augmented with diffuse functions (aug-cc-pVXZ, X=D, T, and Q). Al2H and Si2H are determined to have H-bridged C2v structures in their ground states: the Al2H ground state is of 2B1 symmetry with an Al-H-Al angle of 87.6 degrees, and the Si2H ground state is of 2A1 symmetry with a Si-H-Si angle of 79.8 degrees. However, P2H and S2H have nonbridged, bent Cs structures: the P2H ground state is of 2A' symmetry with a P-P-H angle of 97.0 degrees, and the S2H ground state is of 2A' symmetry with a S-S-H angle of 93.2 degrees. Ground state geometries, vibrational frequencies, and electron affinities have been computed at all levels of theory. Our CCSD(T)/aug-cc-pVQZ adiabatic electron affinity of 2.34 eV for the Si2H radical is in excellent agreement with the photoelectron spectroscopy experiments of Xu et al. [J. Chem. Phys. 108, 7645 (1998)], where the electron affinity was determined to be 2.31+/-0.01 eV.     

An ab initio investigation of the O(3P)-H2(1sigma(g)+) van der Waals well

We report an ab initio study of the van der Waals region of the O(3P)-H2 potential energy surface based on RCCSD(T) calculations with an aug-cc-pVQZ basis supplemented by bond functions. In addition, an open-shell implementation of symmetry-adapted perturbation theory (SAPT) is used to corroborate the RCCSD(T) calculations and to investigate the relative magnitudes of the various contributions to the van der Waals interaction. We also investigate the effect of the spin-orbit coupling on the position and depth of the van der Waals well. We predict the van der Waals minimum to occur in perpendicular geometry, and located at a closer distance than a secondary well in colinear geometry. The potentials obtained in the present study confirm the previous calculations of Alexander [M. H. Alexander, J. Chem. Phys., 1998, 108, 4467], but disagree with the earlier work of Harding and co-workers [Z. Li, V. A. Apkarian and L. B. Harding, J. Chem. Phys., 1997, 106, 942] as well as with recently refitted surfaces of Brand?o and coworkers [J. Brand?o, C. Mogo and B. C. Silva, J. Chem. Phys., 2004, 121, 8861]. Inclusion of spin-orbit coupling reduces the depth of the van der Waals minimum without causing a change in its position.     

Electronic spectroscopy of the A1A" <--> X1A' system of CDBr

We report fluorescence excitation and single vibronic level emission spectra of jet-cooled CDBr in the 450-750 nm region. A total of 32 cold bands involving the pure bending levels 2(0)n with n=3-10 and combination bands 2(0)n3(0)1 (n=2-10), 2(0)n3(0)2 (n=2-9), 1(0)(1)2(0)n (n=7-10), and 1(0)(1)2(0)n3(0)(1) (n=6,8-9) in the A1A" <-- X1A' system of this carbene were observed; most of these are reported and/or rotationally analyzed here for the first time. Rotational analysis yielded band origins and effective (B) rotational constants for both bromine isotopomers (CD79Br and CD81Br). The derived A1A" vibrational intervals are combined with results of Yu et al. [J. Chem. Phys. 115, 5433 (2001)] to derive barriers to linearity for the 2n, 2n3(1), and 2n3(2) progressions. The A1A" state C-D stretching frequency (2350 cm(-1)) is determined for the first time, in excellent agreement with theory, as are the 79Br-81Br isotope splittings in the excited state. Our emission spectra probe the vibrational structure of the X1A' and a3A" states up to approximately 9000 cm(-1) above the vibrationless level of the X1A' state; the total number of levels observed is around twice that previously reported. Unlike CHBr, where even the lowest bending levels are perturbed by spin-orbit interaction with the triplet origin, the term energy of every level save one below 3000 cm(-1) in CDBr is reproduced by a Dunham expansion to within a standard deviation of 1 cm(-1), and a spin-orbit coupling matrix element of approximately 330 cm(-1) is derived from a deperturbation analysis of the triplet origin. The multireference configuration interaction (MRCI) calculations of Yu et al. [J. Chem. Phys. 115, 5433 (2001)] well reproduce triplet perturbations in the pure bending manifold, and globally, the vibrational frequencies of X1A', a3A", and A1A" are in excellent agreement with theoretical predictions.     

Franck-Condon simulation, including anharmonicity, of the photodetachment spectrum of P2H(-): restricted-spin coupled-cluster single-double plus perturbative triple and unrestricted-spin coupled-cluster single-double plus perturbative triple -F12x potential energy functions of P2H and P2H(-)

Geometry optimization and harmonic vibrational frequency calculations have been carried out on the X?(2)A(') state of P(2)H and the X?(1)A(') state of P(2)H(-) using the restricted-spin coupled-cluster single-double plus perturbative triple excitation [RCCSD(T)] and explicitly correlated unrestricted-spin coupled-cluster single-double plus perturbative triple excitation [UCCSD(T)-F12x] methods. For RCCSD(T) calculations, basis sets of up to the augmented correlation-consistent polarized valence quintuple-zeta (aug-cc-pV5Z) quality were employed, and contributions from extrapolation to the complete basis set limit and from core correlation of the P 2s(2)2p(6) electrons were also included. For UCCSD(T)-F12x calculations, different atomic orbital basis sets of triple-zeta quality with different associated complementary auxiliary basis sets and different geminal Slater exponents were used. When the P 2s(2)2p(6) core electrons were correlated in these F12x calculations, appropriate core-valence basis sets were employed. In addition, potential energy functions (PEFs) of the X?(2)A(') state of P(2)H and the X?(1)A(') state of P(2)H(-) were computed at different RCCSD(T) and UCCSD(T)-F12x levels, and were used in variational calculations of anharmonic vibrational wavefunctions, which were then utilized to calculate Franck-Condon factors (FCFs) between these two states, employing a method which includes allowance for anharmonicity and Duschinsky rotation. The photodetachment spectrum of P(2)H(-) was then simulated using the computed FCFs. Simulated spectra obtained using the RCCSD(T)/aug-cc-pV5Z and UCCSD(T)-F12x(x = a or b)/aug-cc-pCVTZ PEFs are compared and found to be essentially identical. Based on the computed FCFs, a more detailed assignment of the observed vibrational structure than previously reported, which includes "hot bands," has been proposed. Comparison between simulated and available experimental spectra has been made, and the currently most reliable sets of equilibrium geometrical parameters for P(2)H and its anion have been derived. The photodetachment spectrum of P(2)D, yet to be recorded, has also been simulated.     

Ab initio potential energy surfaces and nonadiabatic interactions in the H+ +NO collision system

Ab initio calculations on the H(+)+NO system have been carried out in Jacobi coordinates at the multireference configuration interaction level employing Dunning's correlation-consistent polarized valence triple zeta basis set to analyze the role of low-lying electronic excited states in influencing the collision dynamics relevant to the experimental collision energy range of 9.5-30 eV. The lowest two adiabatic potential energy surfaces, asymptotically correlating to H(+)+NO(X (2)Pi) and H((2)S)+NO(+)(X (1)Sigma(+)), have been obtained. Using ab initio procedures, the (radial) nonadiabatic couplings and the mixing angle between the lowest two electronic states (1 (2)A' and 2 (2)A') have been obtained to yield the corresponding quasidiabatic potential energy matrix. The strengths of the computed vibrational coupling matrix elements reflect a similar trend, as has been observed experimentally in the magnitudes of the state-to-state transition probability for the inelastic vibrational excitations [J. Krutein and F. Linder, J. Chem. Phys. 71, 559 (1979); F. A. Gianturco et al., J. Phys. B 14, 667 (1981)].     

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.