Elementary constituents of microdevices: the Ge2H fragment

Highly correlated ab initio electronic structure theory has been used to systematically investigate the linear (X 2Pi) GeGeH and H-bridged (X 2B1 and A 2A1) GeHGe structures and the isomerization transition state (A 2A') connecting X 2Pi with A 2A1. The equilibrium structures and physical properties have been predicted employing self-consistent field, configuration interaction with single and double excitations, coupled cluster with single and double excitations (CCSD), and CCSD with perturbative triple excitations [CCSD(T)]. Four correlation-consistent polarized valence-[cc-pVXZ and aug-cc-pVXZ (X=T and Q)] type basis sets were used in the study. For the H-bridged GeHGe isomer, the X 2B1 ground state is predicted to lie only 0.74 [0.95 with zero-point vibrational energy (ZPVE) corrections] kcal mol-1 below the A 2A1 excited state at the CCSD(T) level of theory with the augmented correlation-consistent polarized valence quadruple zeta (aug-cc-pVQZ) basis set. The X 2B1 state of the H-bridged GeHGe isomer lies 8.6 kcal mol-1 below the X 2Pi ground state of the linear GeGeH isomer. The forward isomerization barrier from the linear ground state to the A 2A' transition state is predicted to be 3.7 kcal mol-1. The reverse isomerization barrier between the A 2A1 GeHGe structure and the X 2Pi GeGeH structure is predicted to be 11.6 (10.8 with the ZPVE corrections) kcal mol-1 at the aug-cc-pVQZ CCSD(T) level of theory.     

Characterization of singlet ground and low-lying electronic excited states of phosphaethyne and isophosphaethyne

The singlet ground ((approximate)X(1)Sigma1+) and excited (1Sigma-,1Delta) states of HCP and HPC have been systematically investigated using ab initio molecular electronic structure theory. For the ground state, geometries of the two linear stationary points have been optimized and physical properties have been predicted utilizing restricted self-consistent field theory, coupled cluster theory with single and double excitations (CCSD), CCSD with perturbative triple corrections [CCSD(T)], and CCSD with partial iterative triple excitations (CCSDT-3 and CC3). Physical properties computed for the global minimum ((approximate)X(1)Sigma+HCP) include harmonic vibrational frequencies with the cc-pV5Z CCSD(T) method of omega1=3344 cm(-1), omega2=689 cm(-1), and omega3=1298 cm(-1). Linear HPC, a stationary point of Hessian index 2, is predicted to lie 75.2 kcal mol(-1) above the global minimum HCP. The dissociation energy D0[HCP((approximate)X(1)Sigma+)-->H(2S)+CP(X2Sigma+)] of HCP is predicted to be 119.0 kcal mol(-1), which is very close to the experimental lower limit of 119.1 kcal mol(-1). Eight singlet excited states were examined and their physical properties were determined employing three equation-of-motion coupled cluster methods (EOM-CCSD, EOM-CCSDT-3, and EOM-CC3). Four stationary points were located on the lowest-lying excited state potential energy surface, 1Sigma- -->1A", with excitation energies Te of 101.4 kcal mol(-1) (1A"HCP), 104.6 kcal mol(-1)(1Sigma-HCP), 122.3 kcal mol(-1)(1A" HPC), and 171.6 kcal mol(-1)(1Sigma-HPC) at the cc-pVQZ EOM-CCSDT-3 level of theory. The physical properties of the 1A" state with a predicted bond angle of 129.5 degrees compare well with the experimentally reported first singlet state ((approximate)A1A"). The excitation energy predicted for this excitation is T0=99.4 kcal mol(-1) (34 800 cm(-1),4.31 eV), in essentially perfect agreement with the experimental value of T0=99.3 kcal mol(-1)(34 746 cm(-1),4.308 eV). For the second lowest-lying excited singlet surface, 1Delta-->1A', four stationary points were found with Te values of 111.2 kcal mol(-1) (2(1)A' HCP), 112.4 kcal mol(-1) (1Delta HPC), 125.6 kcal mol(-1)(2(1)A' HCP), and 177.8 kcal mol(-1)(1Delta HPC). The predicted CP bond length and frequencies of the 2(1)A' state with a bond angle of 89.8 degrees (1.707 A, 666 and 979 cm(-1)) compare reasonably well with those for the experimentally reported (approximate)C(1)A' state (1.69 A, 615 and 969 cm(-1)). However, the excitation energy and bond angle do not agree well: theoretical values of 108.7 kcal mol(-1) and 89.8 degrees versus experimental values of 115.1 kcal mol(-1) and 113 degrees. of 115.1 kcal mol(-1) and 113 degrees.     

Characterization of the X2A1, A2B1, and X2Pi electronic states of the Ga2H molecule and the X2A' and A 2A" isomerization transition states connecting the three minima

A wide range of highly correlated ab initio methods has been used to predict the geometrical parameters of the linear (X (2)Pi) and H-bridged (X (2)A(1) and A (2)B(1)) Ga(2)H isomers and two isomerization transition states (X (2)A(') and A (2)A(")) connecting the three minima. Dipole moments and vibrational frequencies are also obtained. The global minimum X (2)A(1) ground state of the H-bridged GaHGa isomer is predicted to lie only 1.6 [1.9 with the zero-point vibrational energy (ZPVE) corrections] kcal mol(-1) below the A (2)B(1) state. The X (2)A(1) state lies 5.4 kcal mol(-1) below the X (2)Pi ground state of the linear GaGaH isomer at the coupled-cluster with single, double, and perturbative triple excitations [CCSD(T)] level of theory with the augmented correlation-consistent polarized valence quadruple-zeta (aug-cc-pVQZ) basis set. The full triples coupled-cluster method is found to alter these CCSD(T) predictions by as much as 0.3 kcal mol(-1). The forward isomerization barriers from the linear ground state to the X (2)A(') and A (2)A(") transition states are determined to be 3.3 and 5.3 kcal mol(-1), respectively. The reverse isomerization barrier between the X (2)A(1) GaHGa structure and the X (2)Pi GaGaH structure is predicted to be 8.6 (8.2 with the ZPVE corrections) kcal mol(-1) at the aug-cc-pVQZ CCSD(T) level of theory.     

Assessment of recently developed exchange-correlation functionals for the description of torsion potentials in pi-conjugated molecules

Newly developed exchange-correlation functionals in density functional theory (DFT) have been applied to describe conjugation effects in organic molecules. The performance of the various approaches is assessed through the calculation of torsion energy profiles and their critical comparison with available experimental data. Our results indicate that the OPTX-B95 exchange-correlation functional as well as its corresponding hybrid versions perform better than the well-established BLYP or B3LYP schemes when dealing with pi-conjugated molecules. In contrast, the recently introduced VSXC functional is not as reliable as other DFT methods for the systems examined here.     

Study of the C(3P) + OH(X2Pi) --> CO(X1Sigma(g)+) + H(2S) reaction: a fully global ab initio potential energy surface of the X2A' state

The C((3)P) + OH(X (2)Pi) --> CO(X (1)Sigma(g)(+)) + H((2)S) reaction has been investigated by ab initio electronic structure calculations of the X(2)A' state based on the multireference (MR) internally contracted single and double configuration interaction (SDCI) method plus Davidson correction (+Q) using Dunning aug-cc-pVQZ basis sets. In particular, the multireference space is taken to be a complete active space (CAS). Improvement over previously proposed potential energy surfaces for HCO/COH is obtained in the sense that present surface describes also the potential part where the CO interatomic distance is large. A large number of geometries (around 2000) have been calculated and analytically fitted using the reproducing kernel Hilbert space (RKHS) method of Ho and Rabitz both for the two-body and three-body terms following the many-body decomposition of the total electronic energies. Results show that the global reaction is highly exothermic ( approximately 6.4 eV) and barrierless (relative to the reactant channel), while five potential barriers are located on this surface. The three minima and five saddle points observed are characterized and found to be in good agreement with previous work. The three minima correspond to the formation of HCO and COH complexes and to the CO + H products, with the COH complex being a metastable minimum relative to the product channel. The five saddle points correspond to potential barriers for both the dissociation/formation of HCO and COH into/from CO + H, to barriers for the isomerization of HCO into COH and to barriers for the inversion of HCO and COH through their respective linear configuration.     

PdO0,±1, PdH0,±1及PdOH分子结构和电子性质的从头计算

利用多种从头计算方法对PdO0,±1, PdH0,±1以及PdOH的结构和电子性质进行了理论研究. 优化了PdO分子两个能量相近的态(3∏和3∑-), 结果表明PdO 的3∏态较3∑-态更为稳定, 在单双激发耦合簇理论(CCSD)水平下3∑-态的总能量比3∏态的高出0.286 eV. 在结构优化的基础上计算了PdO以及PdH的绝热电离能(AIE)和电子亲和能(AEA), 计算结果与实验值符合得很好. PdOH的基态为Cs对称性的角型结构Pd—OH(2A'态), 另外还优化得到两个2A'态的亚稳结构, 分别对应于Pd—OH和O—Pd—H, CCSD水平下两个亚稳态的总能量较基态分别高0.405和2.284 eV, 优化得到了连接这两个2A'态的过渡态, 并计算了相应的反应能垒.     

Ab initio and DFT benchmark study for nucleophilic substitution at carbon (SN2@C) and silicon (SN2@Si)

To obtain a set of consistent benchmark potential energy surfaces (PES) for the two archetypal nucleophilic substitution reactions of the chloride anion at carbon in chloromethane (S(N)2@C) and at silicon in chlorosilane (S(N)2@Si), we have explored these PESes using a hierarchical series of ab initio methods [HF, MP2, MP4SDQ, CCSD, CCSD(T)] in combination with a hierarchical series of six Gaussian-type basis sets, up to g polarization. Relative energies of stationary points are converged to within 0.01 to 0.56 kcal/mol as a function of the basis-set size. Our best estimate, at CCSD(T)/aug-cc-pVQZ, for the relative energies of the [Cl(-), CH(3)Cl] reactant complex, the [Cl-CH(3)-Cl](-) transition state and the stable [Cl-SiH(3)-Cl](-) transition complex is -10.42, +2.52, and -27.10 kcal/mol, respectively. Furthermore, we have investigated the performance for these reactions of four popular density functionals, namely, BP86, BLYP, B3LYP, and OLYP, in combination with a large doubly polarized Slater-type basis set of triple-zeta quality (TZ2P). Best overall agreement with our CCSD(T)/aug-cc-pVQZ benchmark is obtained with OLYP and B3LYP. However, OLYP performs better for the S(N)2@C overall and central barriers, which it underestimates by 2.65 and 4.05 kcal/mol, respectively. The other DFT approaches underestimate these barriers by some 4.8 (B3LYP) to 9.0 kcal/mol (BLYP).     

Describing weak interactions of biomolecules with dispersion-corrected density functional theory

Interaction energies of the biomolecules in the JSCH-2005 database are calculated with density functional theory using the exchange-correlation functional BLYP augmented with dispersion-corrected atom-centered potentials (DCACPs). The results are in excellent agreement with extrapolated CCSD(T) complete basis set limit references with unsigned mean errors of less than 1.6 kcal mol(-1). Geometry optimisations all reach stable configurations that are close to the MP2-optimised reference geometries.     

Microsolvation of protonated methane: structures and energetics of CH5(+)(H2)n

Effects of microsolvating CH(5)(+) with up to four H(2) molecules have been investigated in terms of structures and energies. For the smaller complexes, benchmark calculations have been carried out using MP2 and CCSD(T) with basis sets up to aug-cc-pV5Z quality and energies have been extrapolated to the infinite basis set limit. It is found that MP2 calculations using the aug-cc-pVQZ basis set or better yield robust reference data for both structures and energies. More than 30 stationary points including minima and first-order as well as second-order stationary points have been characterized by this method and are discussed in terms of solvation motifs. Finally, the performance of several density functionals has been assessed for this very demanding case. Popular GGA functionals such as BLYP and PBE fail, whereas the TPSS meta-GGA functional captures many structural and energetic aspects of microsolvation satisfactorily.     

ClO/ClO~-耦合体系结构性质及反应通道研究

利用从头算法和密度泛函理论对ClO/ClO~-体系进行了研究。首先应用密度泛 函理论的六种方法(B3LYP,BLYP,B3P86,BP86,BHLYP,LSDA)和从头算理论的CCSD方法 在6-3+G~*,6-311+G~*及6-311+G(3df)基组水平上对单体进行了优化,选出最适合 该体系的方法和基组B3LYP/6-311+G(3df);然后在B3LYP/6-311+G(3df)水平上计算 了沿各种反应通道的离解能,并且进行了校正。找出了存在的中间体及中间体异构 化的过渡态,并进行了IRC路径解析。对各稳定体进行了频率分析和成键分析。结 果表明,单体ClO和单体ClO~-结合为稳定的中间体后,其离解方式主要是向着 ClOO+Cl~-和ClOO~-+Cl两种方式进行,两种离解方式的离解能分别为-33.39和82. 88kJ/mol,并且前者是一个离解性电子转移过程,经历一个电子转移过渡态。     

Static dipole polarizability and binding energy of sodium clusters Nan (n=1-10): A critical assessment of all-electron based post Hartree-Fock and density functional methods

A systematic all electron post Hartree-Fock as well as density functional theory (DFT) based calculations for the polarizability and binding energy of sodium metal clusters have been performed and an in-depth analysis of the discrepancy between the experimental and theoretical results is presented. A systematic investigation for the assessment of different DFT exchange-correlation functionals in predicting the polarizability values has also been reported. All the pure DFT functionals have been found to considerably underestimate the calculated polarizability values as compared to the MP2 results. DFT calculations using the full Hartree-Fock exchange along with one-parameter progressive correlation functional have, however, been shown to yield results in good agreement with the MP2 and experimental results. The possible sources of error present in the experimental measurements as well as in the different theoretical methods have also been analyzed. One of the most important conclusions of the present study is that the effect of electron correlation plays a significant role in determining the polarizability of the clusters and the MP2 method can be considered to be one of the most reliable methods for their prediction. It has also been noted that the polarizability value of the lower member clusters (Na2 and Na4) calculated by highly sophisticated methods such as, CCSD and CCSD(T) are found to be very close to the corresponding MP2 values. The polarizability and the binding energy of the clusters are found to be inversely related to each other and their correlation is rationalized by invoking the minimum polarizability principle. A good linear correlation between the polarizability and volume of the cluster has also been found to exist.     

On the accuracy of density-functional theory exchange-correlation functionals for H bonds in small water clusters: benchmarks approaching the complete basis set limit

The ability of several density-functional theory (DFT) exchange-correlation functionals to describe hydrogen bonds in small water clusters (dimer to pentamer) in their global minimum energy structures is evaluated with reference to second order Moller-Plesset perturbation theory (MP2). Errors from basis set incompleteness have been minimized in both the MP2 reference data and the DFT calculations, thus enabling a consistent systematic evaluation of the true performance of the tested functionals. Among all the functionals considered, the hybrid X3LYP and PBE0 functionals offer the best performance and among the nonhybrid generalized gradient approximation functionals, mPWLYP and PBE1W perform best. The popular BLYP and B3LYP functionals consistently underbind and PBE and PW91 display rather variable performance with cluster size.     

Low-lying quartet electronic states of nitrogen dioxide

The environmentally active molecule nitrogen dioxide (NO2) has been systematically studied using high level theoretical methods. The electronic ground state and the low-lying quartet states of NO2 have been investigated. Single reference restricted open-shell self-consistent field (SCF), complete active space SCF (CASSCF), spin-restricted (R) and spin-unrestricted (U) configuration interaction with single and double excitations (CISD), coupled cluster with single and double excitations (CCSD), CCSD with perturbative triple excitations [CCSD(T)], and internally contracted multireference configuration interaction (ICMRCI) methods along with Dunning's correlation consistent polarized valence cc-pVXZ and augmented cc-pVXZ (where X=T,Q,5) basis sets were used in this research. At the aug-cc-pV5Z/UCCSD(T) level the classical adiabatic excitation energies (Te values) of the three lowest-lying quartet excited states were predicted to be 83.3 kcalmol (3.61 eV, 29 200 cm(-1)) for the ? 4A2 state, 93.3 kcalmol (4.05 eV, 32 600 cm(-1)) for the b 4B2 state, and 100.8 kcalmol (4.37 eV, 35 300 cm(-1)) for the c 4A1 state. The quantum mechanical excitation energies (T 0 values) were determined to be 81.6 kcalmol (3.54 eV, 28 500 cm(-1)) for the a 4A2 state and 90.7 kcalmol (3.93 eV, 31 700 cm(-1)) for the b 4B2 state. The lowest quartet linear Renner-Teller 4Pi state gives rise to the a 4A2 state with 112.8 degrees and the b 4B2 state with 124.4 degrees <(ONO) bond angles upon bending. The b state shows some peculiar behavior. Although CASSCF, RCISD, UCISD, RCCSD, UCCSD, and RCCSD(T) methods predicted the presence of a Cs equilibrium geometry (a double minimum 4A' state), SCF, UCCSD(T), and ICMRCI wave functions predicted the C2v structure for the b 4B2 state. The importance of both dynamical and nondynamical correlation treatments for the energy difference between C2v and Cs structures of b state is highlighted in this context. The c 4A1 state is predicted to have a very small bond angle of 85.8 degrees . Potential energy diagrams with respect to the bond angles of the ground state and four quartet states are presented.     

A CASSCF/MR-CI study toward the understanding of wavelength-dependent and geometrically memorized photodissociation of formic acid

The S(0), T(1), and S(1) potential energy surfaces for the HCOOH dissociation and isomerization processes have been mapped with different ab initio methods. The wavelength-dependent mechanism for the HCOOH dissociation was elucidated through the computed potential energy surfaces and the surface crossing points. The HCOOH molecules in S(1) by excitation at 248 nm mainly decay to the ground state via the S(0) and S(1) vibronic interaction, followed by molecular eliminations in the ground state. The S(1) direct dissociation to HCO((2)A') + OH((2)Pi) is the dominant pathway upon photoexcitation at 240-210 nm. Meanwhile, there is a slight probability that the system relaxes to the ground state via the S(0) and S(1) vibronic interaction at these wavelengths. After irradiation of HCOOH at 193 nm, the S(1) direct dissociation into HCO((2)A') + OH((2)Pi) is energetically the most favorable pathway. In view of high IC efficiency at the S(0)/S(1) conical crossing, the S(1) --> S(0) internal conversion via the S(0)/S(1) point can occur with considerable efficiency. In addition, the S(1) isomerization probably plays a dominant role in the partially conformational memory of the HCOOH photodissociation, which has been discussed in detail.     

Performance of multi-configurational calculations for a 1,4-bis(phenylethynyl)benzene derivative conjugated molecule.

The theoretical challenge of finding a single method that quantitatively reproduces both the experimental low-lying excitation energies and the torsional barrier of a prototypical conjugated molecule, which could act as a molecular wire, has been addressed here. The results indicate that this goal can be reasonably achieved when multi-reference perturbation theory up to second order (MRMP2) based on a complete active space self-consistent field (CASSCF) wave function using large active spaces is used. The results obtained were also used to compare with less expensive Kohn-Sham (KS) density functional theory (DFT) calculations when applied to these properties. The results obtained with BLYP and B3LYP exchange-correlation functionals indicate that quantitative agreement with all the experimental data cannot be obtained with this methodology, with a clear dependence on the exchange-correlation form selected. We thus encourage a careful testing of pure and hybrid density functionals whenever KS DFT is used for the rational design of conjugated materials for charge conduits.     

Charge-transfer excitations in uranyl tetrachloride ([UO2Cl4]2-): how reliable are electronic spectra from relativistic time-dependent density functional theory?

Four-component relativistic time-dependent density functional theory (TD-DFT) is used to study charge-transfer (CT) excitation energies of the uranyl molecule as well as the uranyl tetrachloride complex. Adiabatic excitation energies and vibrational frequencies of the excited states are calculated for the lower energy range of the spectrum. The results for TD-DFT with the CAM-B3LYP exchange-correlation functional for the [UO(2)Cl(4)](2-) system are in good agreement with the experimentally observed spectrum of this species and agree also rather well with other theoretical data. Use of the global hybrid B3LYP gives qualitatively correct results, while use of the BLYP functional yields results that are qualitatively wrong due to the too low CT states calculated with this functional. The applicability of the overlap diagnostic of Peach et al. (J. Chem. Phys.2008, 128, 044118) to identify such CT excitations is investigated for a wide range of vertical transitions using results obtained with three different approximate exchange-correlation functionals: BLYP, B3LYP, and CAM-B3LYP.     

Structures and stability of medium-sized silicon clusters. III. Reexamination of motif transition in growth pattern from Si15 to Si20

It has been established from experiments that stable medium-sized ionic clusters Si15-Si20 are prolate in shape. Density-functional theories (DFTs) also predict that nearly all low-lying neutral clusters in this size range are prolate in shape. Moreover, most of them are built onto two generic structural motifs, either the tricapped-trigonal-prism (TTP) Si9 motif or the six/six Si6Si6 (sixfold-puckered hexagonal ring Si6 plus six-atom tetragonal bipyramid Si6) motif. However, it appears that the exact location of the TTP-to-six/six motif transition is dependent on the functional (e.g., PBE or BLYP) used in the DFT calculations. Here, we present total-energy calculations for two series of clusters (one series containing six/six motif and the other containing the TTP motif) in the size range of Si16-Si20. The calculations were based on all-electron DFT methods with a medium [6-311G (2d)] and a large (cc-pVTZ) basis sets, as well as coupled-cluster single and double substitutions (including triple excitations) [CCSD(T)] method with a modest (cc-pVDZ) basis set. In the DFT calculations, two popular hybrid density functionals, the B3LYP and PBE1PBE, were selected. It is found that the B3LYP total-energy calculations slightly favor the six/six motif, whereas the PBE1PBE calculations slightly favor the TTP motif. The CCSD(T) total-energy calculations, however, show that isomers based on the six/six motif are energetically slightly favorable in the size range of Si16-Si20. Hence, the TTP-to-six/six motif transition is more likely to occur at Si16.     

Pi and sigma-phenylethynyl radicals and their isomers o-, m-, and p-ethynylphenyl: structures, energetics, and electron affinities

Molecular structures, energetics, vibrational frequencies, and electron affinities are predicted for the phenylethynyl radical and its isomers. Electron affinities are computed using density functional theory, -namely, the BHLYP, BLYP, B3LYP, BP86, BPW91, and B3PW91 functionals-, employing the double-zeta plus polarization DZP++ basis set; this level of theory is known to perform well for the computation of electron affinities. Furthermore, ab initio computations employing perturbation theory, coupled cluster with single and double excitations [CCSD], and the inclusion of perturbative triples [CCSD(T)] are performed to determine the relative energies of the isomers. These higher level computations are performed with the correlation consistent family of basis sets cc-pVXZ (X = D, T, Q, 5). Three electronic states are probed for the phenylethynyl radical. In C2v symmetry, the out-of-plane (2B1) radical is predicted to lie about 10 kcal/mol below the in-plane (2B2) radical by DFT methods, which becomes 9.4 kcal/mol with the consideration of the CCSD(T) method. The energy difference between the lowest pi and sigma electronic states of the phenylethynyl radical is also about 10 kcal/mol according to DFT; however, CCSD(T) with the cc-pVQZ basis set shows this energy separation to be just 1.8 kcal/mol. The theoretical electron affinities of the phenylethynyl radical are predicted to be 3.00 eV (B3LYP/DZP++) and 3.03 eV (CCSD(T)/DZP++//MP2/DZP++). The adiabatic electron affinities (EAad) of the three isomers of phenylethynyl, that is, the ortho-, meta-, and para-ethynylphenyl, are predicted to be 1.45, 1.40, and 1.43 eV, respectively. Hence, the phenylethynyl radical binds an electron far more effectively than the three other radicals studied. Thermochemical predictions, such as the bond dissociation energies of the aromatic and ethynyl C-H bonds and the proton affinities of the phenylethynyl and ethynylphenyl anions, are also reported.     

Protonated carbonyl sulfide: prospects for the spectroscopic observation of the elusive HSCO+ isomer

Spurred by the apparent conflict between ab initio predictions and infrared spectroscopic evidence regarding the relative stability of isomers of protonated carbonyl sulfide, key stationary points on the isomerization surface of HOCS(+) have been examined via systematic extrapolations of ab initio energies. Electron correlation has been accounted for using second-order M?ller-Plesset perturbation theory and coupled cluster theory through triple excitations [CCSD, CCSD(T), and CCSDT] in conjunction with the correlation consistent hierarchy of basis sets, cc-pVXZ (X=D,T,Q,5,6). HSCO(+) is predicted to lie lower in energy than HOCS(+) by 4.86 kcal mol(-1), computed using the focal point extrapolation scheme of Allen and co-workers [J. Chem. Phys. 99, 4638 (1993)] with corrections for anharmonic zero-point vibrational energy, core correlation, non-Born-Oppenheimer, and scalar relativistic effects. A transition state has been located, constituting the barrier to isomerization of HSCO(+) to HOCS(+), lying 68.9 kcal mol(-1) higher in energy than HSCO(+). This is well above predicted exothermicity [DeltaH(r) (o)(0 K)=48.1 kcal mol(-1), cc-pVQZ CCSD(T)] for the reaction considered in the experiments (HSCO(+)+H(2)-->OCS+H(3) (+)). Though proton tunneling will lead to a lower effective barrier, this prediction is consistent with the lack of HSCO(+) in electrical discharges in H(2)OCS, since the relative populations of HOCS(+) and HSCO(+) will depend on the experimental details of the protonation route rather than the relative thermodynamic stability of the isomers. Anharmonic vibrational frequencies and vibrationally corrected rotational constants from cc-pVTZ CCSD(T) cubic and quartic force constants are provided, to aid in the spectroscopic observation of the energetically favorable but apparently elusive HSCO(+) isomer.