Electronic states and potential energy curves of molybdenum carbide and its ions

The potential energy curves and spectroscopic constants of the ground and 29 low-lying excited states of MoC with different spin and spatial symmetries within 48 000 cm(-1) have been investigated. We have used the complete active space multiconfiguration self-consistent field methodology, followed by multireference configuration interaction (MRCI) methods. Relativistic effects were considered with the aid of relativistic effective core potentials in conjunction with these methods. The results are in agreement with previous studies that determined the ground state as X (3)Sigma(-). At the MRCISD+Q level, the transition energies to the 1 (3)Delta and 4 (1)Delta states are 3430 and 8048 cm(-1), respectively, in fair agreement with the results obtained by DaBell et al. [J. Chem. Phy. 114, 2938 (2001)], namely, 4003 and 7834 cm(-1), respectively. The three band systems located at 18 611, 20 700, and 22 520 cm(-1) observed by Brugh et al. [J. Chem. Phy. 109, 7851 (1998)] were attributed to the excited 11 (3)Sigma(-), 14 (3)Pi, and 15 (1)Pi states respectively. At the MRCISD level, these states are 17 560, 20 836, and 20 952 cm(-1) above the ground state respectively. We have also identified a (3)Pi state lying 14 309 cm(-1) above the ground state. The ground states of the molecular ions are predicted to be (4)Sigma(-) and (2)Delta for MoC(-) and MoC(+), respectively.     

Theoretical study of electronic states of N22+in an intense radiation field

The Floquet states of N(2) (2+) created by the interactions of the six lowest singlet (1 (1)Sigma(g) (+), 1 (1)Delta(g), 2 (1)Sigma(g) (+), 1 (1)Pi(u), 1 (1)Pi(g), and 1 (1)Sigma(u) (-)) states of the dication with intense (0.4 x 10(13) Wcm(2)) radiation have been studied using the recently developed multireference configuration interaction method with single and double excitations (MRCISD)-based approach. The adiabatic Floquet state coinciding near its minimum with the initial X (1)Sigma(g) (+) ground state and asymptotically correlating with A (1)Pi(u) (m = -1), i.e., with one less photon in the dressed state, is expected to be metastable, as is the ground state in the absence of a field, but to support up to the v(max) = 12 quasibound vibrational level in comparison with v(max) = 11 in the parent field-free X (1)Sigma(g) (+) ground state. The tunneling lifetimes of the highest vibrational levels in this adiabatic Floquet state are predicted to be several orders longer than those in the parent field-free state. Analysis of the complete basis set limit extrapolated MRCISD potential energy curve of the field-free X (1)Sigma(g) (+) state of N(2) (2+) calculated in the present work (R(e) = 1.130 A, omega(e) = 2011 cm(-1), omega(e)x(e) = 26.1 cm(-1)) is in good agreement with spectroscopic experimental data. Calculations on the field-free A (1)Pi(u) state (T(e) = 12 106 cm(-1), R(e) = 1.252 A, omega(e) = 1438 cm(-1), omega(e)x(e) = 23.5 cm(-1)) generally support earlier theoretical work and do not support reported experimental values.     

Multireference configuration interaction study of the electronic states of ZrC

The potential energy curves and spectroscopic constants of the ground and 32 low-lying electronic states of ZrC have been studied by employing multireference configuration interaction methods, in conjunction with relativistic effective core potentials and 5s3p3d1f, 3s3p1d basis sets con Zr and C, respectively. We have determined that the ground state is (3)Sigma(+). However there are two low-lying (1)Sigma(+) states (below 5000 cm(-1)) which strongly interact resulting in avoided crossings. The lowest (1)Sigma(+) state corresponds to a combination of 1sigma(2) Xsigma(2) 1pi(4) configurations whereas the second is an open shell singlet 1sigma(2) 2sigma(1) 3sigma(1) 1pi(4). Several avoided crossings were observed, for (1)Pi, (3)Pi, (1)Delta, (3)Sigma(+), and (3)Delta states. We have identified (3)Pi and (1)Pi lying at 4367 and 5797 cm(-1), respectively. The results are in good agreement with the recent experimental findings of Rixon et al. [J. Mol. Spectrosc. 228, 554 (2004)], and indicate that the (3)Pi-(3)Sigma(+), and (1)Pi-(1)Sigma(+), bands located between 16 000-19 000 cm(-1) are extremely complex due to near degeneracy of several (1)Pi and (3)Pi states. We also have identified a (1)Sigma(+) state in the same region that may interfere with the (1)Pi emission bands. The present results not only shed further light into the spectra of ZrC but also predict yet to be observed systems.     

On the ground and some low-lying excited states of ScB: a multiconfigurational study

The electronic structure of a series of low-lying excited triplet and quintet states of scandium boride (ScB) was examined using multireference configuration interaction (including Davidson's correction for quadruple excitations) and single-reference coupled cluster (CC) methods with averaged natural orbital (ANO) basis sets. The CC approach was used only for the lowest quintet state. The authors have analyzed eight low-lying triplets 3Sigma-(2), 3Sigma+, 3Pi(3), and 3Delta(2) dissociating to Sc(2D)/B(2P) atoms and eight low-lying quintet states 5Sigma-, 5Sigma+, 5Pi(2), 5Phi, and 5Delta(3) dissociating to Sc(4F)/B(2P) atoms. They report the potential energy curves and spectroscopic parameters of ScB obtained with the multireference configuration interaction (MRCI) technique including all singly and doubly excited configurations obtained with the ANO-S basis set. For the two lowest states they obtained also improved ANO-L spectroscopic constants, dipole and quadrupole moments as well as scalar relativistic effects based on the Douglas-Kroll-Hess Hamiltonian. They provide the analysis of the bonding based on Mulliken populations and occupation numbers. Since the two lowest states, 3Sigma- and 5Sigma-, lie energetically very close, their principal goal was to resolve the nature of the ground state of ScB. Their nonrelativistic MRCI(Q) (including Davidson correction) results indicate that the quintet is more stable than the triplet by about 800 cm(-1). Inclusion of scalar relativistic effects reduces this difference to about 240 cm(-1). The dissociation energies for 5Sigma- ScB range from 3.20 to 3.30 eV while those for the 3Sigma- range from 1.70 to 1.80 eV.     

Interactions of transition metal atoms in high-spin states: Cr2, Sc-Cr, and Sc-Kr

The high-spin van der Waals states are examined for the following dimers: Cr(2) ((13)Sigma(g)(+)), Sc-Cr ((8)Sigma(+), (8)Pi, (8)Delta), and Sc-Kr ((2)Sigma(+), (2)Pi, (2)Delta). These three systems offer a wide range of van der Waals interactions: anomalously strong, intermediate, and typically weak. The single-reference [coupled cluster with single, double, and noniterative triple excitations, RCCSD(T)] method is used in the calculations for all three systems. In addition, a range of configuration-interaction based methods is applied in Cr(2) and Sc-Cr. The three dimers are shown to be bound by the dispersion interaction of varying strength. In a related effort, the dispersion energy and its exchange counterpart are calculated using the newly developed open-shell variant of the symmetry-adapted perturbation theory (SAPT). The restricted open-shell time-dependent Hartree-Fock linear response function is used in the calculations of the dispersion energy in Sc-Cr and Sc-Kr calculations, while the restricted open-shell time-dependent density functional linear response function is used for Cr(2). A hybrid method combining the repulsive restricted open-shell Hartree-Fock (or complete active space self-consistent field) interaction energy with the dispersion and exchange-dispersion terms is tested against the RCCSD(T) results for the three complexes. The Cr(2) ((13)Sigma(g)(+)) complex has the well depth of 807.8 cm(-1) at the equilibrium distance of 6.18a(0) and the dissociation energy of 776.8 cm(-1). The octet-state Sc-Cr is about four times more strongly bound with the order of well depths of (8)Delta>(8)Pi>(8)Sigma(+) and a considerable anisotropy. The enhanced bonding is attributed to the unusually strong dispersion interaction. Sc-Kr ((2)Sigma(+), (2)Pi, (2)Delta) is a typical van der Waals dimer with well depths in the range of 81 cm(-1) ((2)Delta), 84 cm(-1) ((2)Sigma(+)), and 86 cm(-1) ((2)Pi). The hybrid model based on SAPT leads to results which are in excellent qualitative agreement with RCCSD(T) for all three interactions.     

A theoretical study of the excited states of CrH: potential energies, transition moments, and lifetimes

Ab initio calculations of low-lying electronic states of CrH are presented, including potential energies, dipole and transition dipole moment (TDM) functions, and radiative lifetimes for X (6)Sigma(+), A (6)Sigma(+), 3 (6)Sigma(+), 1 (6)Pi, 2 (6)Pi, 3 (6)Pi, and (6)Delta. Calculation of dynamic correlation effects was performed using the multistate complete active space second-order perturbation method, based on state-averaged complete active space self-consistent-field reference wave functions obtained with seven active electrons in an active space of 16 molecular orbitals. A relativistic atomic natural orbital-type basis set from the MOLCAS library was used for Cr. Good agreement is found between the current calculations and experiment for the lowest two (6)Sigma(+) states, the only states for which spectroscopic data are available. Potential curves for the 3 (6)Sigma(+) and 2 (6)Pi states are complicated by avoided crossings with higher states of the same symmetry, thus resulting in double-well structures for these two states. The measured bandhead T(0)=27 181 cm(-1), previously assigned to a (6)Pi<--X (6)Sigma(+) transition, is close to our value of T(0)=28 434 cm(-1) for the 2 (6)Pi state. We tentatively assign the ultraviolet band found experimentally at 30 386 cm(-1) to the 3 (6)Pi<--X (6)Sigma(+) transition for which the computed value is 29 660 cm(-1). The A (6)Sigma(+)<--X (6)Sigma(+) TDM and A (6)Sigma(+) lifetimes are found to be in reasonable agreement with previous calculations.     

Ab initio study of the KrH+ photodissociation

The multireference spin-orbit configuration interaction method is employed to calculate potential energy curves for the ground and low-lying excited states of the KrH(+) cation. For the first time, the spin-orbit interaction is taken into account and electric dipole moments are computed for transitions to the states responsible for the first absorption continuum (A band) of KrH(+). On this basis, the partial and total absorption spectra in this energy range are obtained. It is shown that the A-band absorption is dominated by the parallel A (1)Sigma(+)<--X (1)Sigma(+) transition. In the low-energy part of the band (<83x10(3) cm(-1)) the absorption is mainly caused by the spin-forbidden b (3)Pi(0(+) )<--X (1)Sigma(+) excitation, while perpendicular transitions to the B (1)Pi and b (3)Pi(1) states are significantly weaker. The branching ratio Gamma for the photodissociation products is calculated and it is shown to increase smoothly from 0 in the red tail of the band to 1 at E>or=90x10(3) cm(-1). The latter value corresponds to the exclusive formation of the spin-excited Kr(+)((2)P(12)) ions, which may be used to obtain laser generation on the Kr(+)((2)P(12)-(2)P(32)) transition.     

Quadrupole, octopole, and hexadecapole electric moments of Sigma, Pi, Delta, and Phi electronic states: cylindrically asymmetric charge density distributions in linear molecules with nonzero electronic angular momentum

The number of independent components, n, of traceless electric 2(l)-multipole moments is determined for C(infinity v) molecules in Sigma(+/-), Pi, Delta, and Phi electronic states (Lambda=0,1,2,3). Each 2(l) pole is defined by a rank-l irreducible tensor with (2l+1) components P(m)((l)) proportional to the solid spherical harmonic r(l)Y(m)(l)(theta,phi). Here we focus our attention on 2(l) poles with l=2,3,4 (quadrupole Theta, octopole Omega, and hexadecapole Phi). An important conclusion of this study is that n can be 1 or 2 depending on both the multipole rank l and state quantum number Lambda. For Sigma(+/-)(Lambda=0) states, all 2(l) poles have one independent parameter (n=1). For spatially degenerate states--Pi, Delta, and Phi (Lambda=1,2,3)--the general rule reads n=1 for l<2/Lambda/ (when the 2(l)-pole rank lies below 2/Lambda/ but n=2 for higher 2(l) poles with l>or=2/Lambda/. The second nonzero term is the off-diagonal matrix element [formula: see text]. Thus, a Pi(Lambda=1) state has one dipole (mu(z)) but two independent 2(l) poles for l>or=2--starting with the quadrupole [Theta(zz),(Theta(xx)-Theta(yy))]. A Delta(Lambda=2) state has n=1 for 2((1,2,3)) poles (mu(z),Theta(zz),Omega(zzz)) but n=2 for higher 2((l>or=4)) poles--from the hexadecapole Phi up. For Phi(Lambda=3) states, it holds that n=1 for 2(1) to 2(5) poles but n=2 for all 2((l>or=6)) poles. In short, what is usually stated in the literature--that n=1 for all possible 2(l) poles of linear molecules--only applies to Sigma(+/-) states. For degenerate states with n=2, all Cartesian 2(l)-pole components (l>or=2/Lambda/) can be expressed as linear combinations of two irreducible multipoles, P(m=0)((l)) and P/m/=2 Lambda)((l)) [parallel (z axis) and anisotropy (xy plane)]. Our predictions are exemplified by the Theta, Omega, and Phi moments calculated for Lambda=0-3 states of selected diatomics (in parentheses): X (2)Sigma(+)(CN), X (2)Pi(NO), a (3)Pi(u)(C(2)), X (2)Delta(NiH), X (3)Delta(TiO), X (3)Phi(CoF), and X (4)Phi(TiF). States of Pi symmetry are most affected by the deviation from axial symmetry.     

Spectroscopy of the ground and low-lying excited states of ThO+

The ThO(+) cation is of interest as it is a useful prototype for experimental and theoretical studies of bonding in a simple actinide compound. Formally the ground state of ThO(+) has the configuration Th(3+)(7s)O(2-), where there is a single unpaired electron associated with a closed-shell Th(4+)-ion core. The first tier of excited states above the X (2)Sigma(+) ground state is expected to be 1 (2)Delta, 1 (2)Pi, and 2 (2)Sigma(+) derived from the Th(3+)(6d)O(2-) configuration. Spectroscopic observations of ThO(+) using the pulsed field ionization-zero kinetic-energy photoelectron technique are reported here. Rotationally resolved spectra were recorded for the X (2)Sigma(+), 1 (2)Delta, and 1 (2)Pi states. Extensive vibrational progressions were observed. Surprisingly, it was found that ionization of ThO decreases the dissociation energy, while increasing the vibrational frequency and decreasing the bond length. Accurate values for the ionization energies of ThO [53 253.8(2) cm(-1)] and Th [50 868.71(8) cm(-1)] were determined as part of this investigation.     

Ab initio study of KN

The potential energy curves for the lowest (3)Sigma(-), (3)Pi, and (5)Sigma(-) states of the KN molecule have been calculated by the multireference singles and doubles configuration interaction method, including Davidson's corrections for quadruple excitations [MRCI(+Q)]. It is shown that the former two are bound, while the last one is repulsive. The electronic ground state of KN is predicted as (3)Sigma(-) state, although the term energy of the (3)Pi state is very small, 177.3 cm(-1). The binding energy for the (3)Sigma(-) state is evaluated as 0.838 eV, the rotational constant B(0) as 0.250 63 cm(-1), and harmonic frequency as 324.4 cm(-1). The spin-orbit coupling effects between the (3)Sigma(-) and (3)Pi states of KN are evaluated and discussed. The same MRCI(+Q) computational procedures are applied to the isovalent LiN, KC, KO, and KCl to confirm the accuracy of present calculations. Theoretical spectroscopic constants presented here will inspire experimental studies of KN.     

Ab initio investigation of the ground and low-lying states of the diatomic fluorides TiF, VF, CrF, and MnF

The electronic structure of the ground and low-lying states of the diatomic fluorides TiF, VF, CrF, and MnF was examined by multireference and coupled cluster methods in conjunction with extended basis sets. For a total of 34 states we report binding energies, spectroscopic constants, dipole moments, separation energies, and charge distributions. In addition, for all states we have constructed full potential curves. The suggested ground state binding energies of TiF(X (4)Phi), VF(X (5)Pi), CrF(X (6)Sigma(+)), and MnF(X (7)Sigma(+)) are 135, 130, 110, and 108 kcal/mol, respectively, with first excited states A (4)Sigma(-), A (5)Delta, A (6)Pi, and a (5)Sigma(+) about 2, 3, 23, and 19 kcal/mol higher. In essence all our numerical findings are in harmony with experimental results. For all molecules and states studied it is clear that the in situ metal atom (M) shows highly ionic character, therefore the binding is described realistically by M(+)F(-).     

Ab initio calculations and vibrational energy level fits for the lower singlet potential-energy surfaces of C3

Ab initio multireference configuration interaction potential energy surfaces are computed for the eight lowest singlet surfaces of C(3). These reveal several important features, including several conical intersections in linear, nonlinear, and equilateral triangle geometries. These intersections are important because, particularly for the excited A (1)Pi(u) state, reasonable ab initio results could only be obtained by including nearby, near degenerate, (1)Sigma(u) (-) and (1)Delta(u) states that cross the A (1)Pi(u) state around 4500 cm(-1) above the equilibrium geometry, and a (1)Pi(g) state whose potential in turn crosses the other states about 2000 cm(-1) further up. These states are probably responsible for the complexity of the shorter wavelength UV absorption spectrum of C(3). The computed potential energy surface for the ground, X (1)Sigma(g) (+), state and for the lowest two excited singlet surfaces (which both correlate with the A (1)Pi(u) state in a collinear geometry) are fitted to analytic functional forms. Vibrational energy levels are calculated for both states, taking account of the Renner-Teller coupling in the excited A (1)Pi(u) state. The potential parameters for both states are then least-squares fitted to experimental data. The ground-state fit covers a range of approximately 8500 cm(-1) above the lowest level, and reproduces 100 observed vibrational levels with an average error of 2.8 cm(-1). The A (1)Pi(u) state surfaces cover a range of 3250 cm(-1) above the zero-point level, and reproduce the 44 observed levels in this range with an average error of 2.8 cm(-1).     

van der Waals interactions and dipole polarizabilities of lanthanides: Tm(2F)-He and Yb(1S)-He potentials

Anisotropic dipole polarizabilities of Tm(2F), Tm+2(2F), and Yb(1S) are calculated using the finite-field multireference averaged quadratic coupled cluster (MR-AQCC) (Tm and Tm+2) and RCCSD(T) (Yb) methods with small-core relativistic pseudopotentials ECP28MWB combined with the augmented ANO basis sets. The lanthanide atoms are strongly polarizable with the scalar part originating from the 6s electrons and the tensorial part from the open 4f shells. The adiabatic interaction potentials 2Sigma+, 2Pi, 2Delta, and 2Phi of Tm(2F)-He and Tm+2(2F)-He were examined by the multireference approaches, multireference configuration interaction and MR-AQCC, using the basis sets designed in the polarizability calculations. A closed-shell lanthanide system Yb(1S)-He was included for comparison. The Tm-He 2Sigma+, 2Pi, 2Delta, and 2Phi interaction potentials are very shallow and nearly degenerate (within 0.01 cm(-1)), with the well depths in the range of 2.35-2.36 cm(-1) at R=6.17 A. The basis-set saturated well depths are expected to be larger by ca. 25%, as estimated using the bond-function augmented basis set. The interactions of lanthanide atoms with He are one order of magnitude less anisotropic than those involving first-row transition metal atoms. The suppression of anisotropy is chiefly attributed to the screening effected by the 6s shell. When these electrons are removed as in the di-cation complex Tm+2(2F)-He, the potentials deepen to a thousand wave number range and their anisotropy is enhanced 500-fold.     

Ab initio configuration-interaction study of the ground and low-lying electronic states of NiI

For the first time, we have studied the potential-energy curves, spectroscopic terms, vibrational levels, and the spectroscopic constants of the ground and low-lying excited states of NiI by employing the complete active space self-consistent-field method with relativistic effective core potentials followed by multireference configuration-interaction calculations. We have identified six low-lying electronic states of NiI with doublet spin multiplicities, including three states of Delta symmetry and three states of Pi symmetry of the molecule within 15 000 cm(-1). The lowest (2)Delta state is identified as the ground state of NiI, and the lowest (2)Pi state is found at 2174.56 cm(-1) above it. These results fully support the previous conclusion of the observed spectra although our computational energy separation of the two states is obviously larger than that of the experimental values. The present calculations show that the low-lying excited states [13.9] (2)Pi and [14.6] (2)Delta are 3 (2)Pi and 3 (2)Delta electronic states of NiI, respectively. Our computed spectroscopic terms, vibrational levels, and spectroscopic constants for them are in good agreement with the experimental data available at present. In the present work we have not only suggested assignments for the observed states but also computed more electronic states that are yet to be observed experimentally.     

Theoretical study of the ArH+ photodissociation

The multireference Spin-Orbit (SO) Configuration Interaction (CI) method in its Lambda-S Contracted SO-CI (LSC-SO-CI) version is employed to calculate potential energy curves for the ground and low-lying excited states of the ArH(+) cation. For the first time, electric dipole moments are also computed in the approach, including SO coupling for transitions to the states responsible for the first absorption continuum (A-band) of ArH(+). On this basis, the partial and total absorption spectra in this energy range are obtained. It is shown that absorption in the A-band is dominated by the parallel A(1)Sigma(+)<--X(1)Sigma(+) transition. In the low-energy part of the band (<95 x 10(3) cm(-1)) the absorption is caused by the perpendicular B(1)Pi<--X(1)Sigma(+) excitation, but transitions to the b(3)Pi(0(+),1) states are also not negligible. The branching ratio Gamma for the final photodissociation products is calculated and it is shown to increase smoothly from 0 in the red tail of the band to 1 at E>or= 10(5) cm(-1). The latter value corresponds to the exclusive formation of the spin-excited Ar(+)((2)P(1/2)) ions, and thus leads to the inverse population of the Ar(+)((2)P(1/2)-(2)P(3/2)) ion states.     

Ab initio investigation of the electronic structure and bonding of BH, BH(-), and HBBH molecules

By correlating all electrons and employing core-tuned correlation consistent basis sets of quintuple-zeta quality, we applied multireference and coupled-cluster methods to study 32 electronic states of the diatomic BH molecule, two bound states of BH(-), and three states of the linear HBBH molecule. We have constructed full potential energy curves and profiles, reporting binding energies, geometries, spectroscopic parameters, dipole moments, and energy separations, whereas our numerical results are in excellent agreement with available experimental numbers. We are trying as well to interpret the binding modes of a large number of the examined states. 18 states of BH are of Rydberg character, with the BH(-) anion revealing similar structural characteristics to the isoelectronic CH species. The first three states of HBBH X 3Sigma g (-), a 1Delta g, and b 1Sigma g + diabatically correlate to two a 3Pi BH fragments, they are similar to the states b 3Sigma g (-), B 1Delta g, and B' 1Sigma g + of the isoelectronic molecule C2, however, their ordering follows that of the first three states of the O2 molecule.     

Insight into the Rydberg states of CH

Ab initio electronic structure calculations of a relatively large number of Rydberg states of the CH radical were carried out employing the multireference single and double excitation configuration interaction (MRD-CI) method. A Gaussian basis set of cc-pV5Z quality augmented with 12 diffuse functions was used together with an extensive treatment of electron correlation. The main focus of this contribution is to investigate the 3d Rydberg complex assigned by Watson [Astrophys. J. 555, 472 (2001)] to three unidentified interstellar bands. The authors' calculations reproduce quite well the absolute excitation energies of the three components of the 3d complex, i.e., 2Sigma+(3dsigma), 2Pi(3dpi), and 2Delta(3ddelta), but not the energy ordering inferred from a rotational assignment of the 3d<--X 2Pi laboratory spectrum. The computation of the 4d complex is reported for the first time along with a number of other higher lying Rydberg species with an X 1Sigma+ core. The lowest Rydberg states belonging to series converging to the a 3Pi and A 1Pi excited states of CH+ are also calculated.     

Slow electron velocity-map imaging spectroscopy of the C4H- and C4D- anions

High resolution photodetachment spectra of C4H- and C4D- obtained via slow electron velocity-map imaging (SEVI) are presented. The spectra reveal closely spaced transitions to the neutral 2Sigma+ and 2Pi states which can be distinguished based on the corresponding photoelectron angular distributions. The C4H ground state is confirmed as the X2Sigma+ state, with the excited A2Pi state lying only 213 cm(-1) higher (201 cm(-1) for C4D). The electron affinities (EAs) are slightly revised to EA (C4H)=28,497+/-8 cm(-1) and EA (C4D)=28,478+/-10 cm(-1). Progressions in low frequency bending vibrations are observed in both states, yielding experimental frequencies of nu7=179(169) cm(-1) and nu6=408(392) cm(-1) for the X2Sigma+ state of C4H (C4D), and nu7=220(215)cm(-1) and nu6=446(437) cm(-1) for the A2Pi state.     

Dissociation of the OCS+ ion in low-lying electronic states studied using multiconfiguration second-order perturbation theory

Complete active space self-consistent-field (CASSCF) and multiconfiguration second-order perturbation theory (CASPT2) calculations with atomic natural orbital basis sets were performed to investigate the S-loss direct dissociation of the 1 2Pi(X 2Pi), 2 2Pi(A 2Pi), 1 2Sigma+(B 2Sigma+), 1 4Sigma-, 1 2Sigma-, and 1 2Delta states of the OCS+ ion and the predissociations of the 1 2Pi, 2 2Pi, and 1 2Sigma+ states. Our calculations indicate that the S-loss dissociation products of the OCS(+) ion in the six states are the ground-state CO molecule plus the S+ ion in different electronic states. The CASPT2//CASSCF potential energy curves were calculated for the S-loss dissociation from the six states. The calculations indicate that the dissociation of the 1 4Sigma- state leads to the CO + S+ (4Su) products representing the first dissociation limit; the dissociations of the 1 2Pi, 1 2Sigma-, and 1 2Delta states lead to the CO + S+(2Du) products representing the second dissociation limit; and the dissociations of the 2 2Pi and 1 2Sigma+ states lead to the CO + S+(2Pu) products representing the third dissociation limit. Seams of the 1 2Pi-1 4Sigma-, 2 2Pi-1 4Sigma-, 2 2Pi-1 2Sigma-, 2 2Pi-1 2Delta, and 1 2Sigma(+)-1 4Sigma- potential energy surface intersections were calculated at the CASPT2 level, and the minima along the seams were located. The calculations indicate that within the experimental energy range (15.07-16.0 eV) the 2 2Pi(A 2Pi) state can be predissociated by 1 4Sigma- forming the S+(4Su) ion and can undergo internal conversion to 1 2Pi followed by the direct dissociation of 1 2Pi forming S+(2Du) and that within the experimental energy range (16.04-16.54 eV) the 1 2Sigma+(B 2Sigma+) state can be predissociated by 1 4Sigma- forming the S+(4Su) ion and can undergo internal conversion to 2 2Pi followed by the predissociation of 2 2Pi by 1 2Sigma- and 1 2Delta forming the S+(2Du) ion. These indications are in line with the experimental fact that both the 4Su and 2Du states of the S+ ion can be formed from the 2 2Pi and 1 2Sigma+ states of the OCS+ ion.