The singlet electronic ground state isomers of dialuminum monoxide: AlOAl, AlAlO, and the transition state connecting them

The singlet electronic ground state isomers, X (1)Sigma(g) (+) (AlOAl D(infinityh)) and X (1)Sigma(+) (AlAlO C(infinitynu)), of dialuminum monoxide have been systematically investigated using ab initio electronic structure theory. The equilibrium structures and physical properties for the two molecules have been predicted employing self-consistent field (SCF) configuration interaction with single and double excitations (CISD), multireference CISD (MRCISD), coupled cluster with single and double excitations (CCSD), CCSD with perturbative triples [CCSD(T)], CCSD with iterative partial triple excitations (CCSDT-3 and CC3), and full triples (CCSDT) coupled cluster methods. Four correlation consistent polarized valence (cc-pVXZ) type basis sets were used. The AlAlO system is rather challenging theoretically. The two isomers are confirmed to have linear structures at all levels of theory. The symmetric isomer AlOAl is predicted to lie 81.9 kcal mol(-1) below the asymmetric isomer AlAlO at the cc-pV(Q+d)Z CCSD(T) level of theory. The predicted harmonic vibrational frequencies for the X (1)Sigma(g) (+) AlOAl molecule, omega(1)=517 cm(-1), omega(2)=95 cm(-1), and omega(3)=1014 cm(-1), are in good agreement with experimental values. The harmonic vibrational frequencies for the X (1)Sigma(+) AlAlO structure, omega(1)=1042 cm(-1), omega(2)=73 cm(-1), and omega(3)=253 cm(-1), presently have no experimental values with which to be compared. With the same methods the barrier heights for the isomerization AlOAl-->AlAlO and AlAlO-->AlOAl reactions were predicted to be 84.3 and 2.4 kcal mol(-1), respectively. The dissociation energies D(0) for AlOAl (X (1)Sigma(g) (+)) and AlAlO (X (1)Sigma(+))-->AlO (X (2)Sigma(+))+Al ((2)P) were determined to be 130.8 and 48.9 kcal mol(-1), respectively. Thus, both symmetric AlOAl (X (1)Sigma(g) (+)) and asymmetric AlAlO (X (1)Sigma(+)) isomers are expected to be thermodynamically stable with respect to the dissociation into AlO (X (2)Sigma(+)) + Al ((2)P) and kinetically stable for the isomerization reaction (AlAlO-->AlOAl) at sufficiently low temperatures.     

Laser induced fluorescence spectroscopy of the C3N radical

Electronic spectra of the C3N radical have been observed for the first time in the near ultraviolet wavelength region by laser induced fluorescence (LIF) spectroscopy. Seventeen vibronic bands of the B 2Pii-X 2Sigma+ electronic transition system of C3N were identified in LIF spectra of products in a discharge of HC3N. The origin of the B 2Pii state was determined to be 27,929.985(1) cm(-1) from rovibrational analyses. It was found that observations of two types of 2Sigma vibronic levels, which have 2Sigma+ and 2Sigma+/- symmetries originated from excitations of the nu4 trans-bending mode (omega4=369.1(20) cm(-1)) with a large Renner-Teller (RT) interaction (epsilon4=-0.1549(50)), and the nu5 cis-bending mode (omega5=163.24(84) cm(-1)) with a small Renner-Teller interaction (epsilon5=-0.0503(68)), respectively. Vibronic levels, with excitations of the C-C stretching (omega3=869.7 cm(-1)) mode, were also identified. The spin-orbit interaction constant was determined to be Aso=-36.7(50) cm(-1) from the RT analysis. In dispersed fluorescence spectra from B 2Pii, vibrational structures of the low-lying electronically excited A 2Pii state were clearly observed with a strong progression due to the nu3' mode, together with those of the X 2Sigma+ state with weak intensities. The origin of A 2Pii, T0=1844(3) cm(-1), and the vibrational frequencies, omega3'=883(3) cm(-1) and omega5'=121(3) cm(-1) for A 2Pii, and omega3"=1054(3) cm(-1), omega4"=405(3) cm(-1), and omega5"=131(3) cm(-1) for X 2Sigma+, were determined. Time profiles of fluorescence from B 2Pii have short (50-200 ns) and long (>1 micros) decay components with quantum beats, indicating that there is a competition between radiative decay and the nonradiative internal conversion to vibrationally highly excited A 2Pii and X 2Sigma+.     

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.     

Application of equation-of-motion coupled-cluster methods to low-lying singlet and triplet electronic states of HBO and BOH

The equilibrium structures and physical properties of the X (1)sigma(+) linear electronic states, linear excited singlet and triplet electronic states of hydroboron monoxide (HBO) (A (1)sigma(-), B (1)delta, a (3)sigma(+), and b (3)delta) and boron hydroxide (BOH) (A (1)sigma(+), B (1)Pi, and b (3)Pi), and their bent counterparts (HBO a (3)A('), b (3)A("), A (1)A("), B (1)A(') and BOH X (1)A('), b (3)A('), c (3)A("), A (1)A('), B (1)A('), C (1)A(")) are investigated using excited electronic state ab initio equation-of-motion coupled-cluster (EOM-CC) methods. A new implementation of open-shell EOM-CC including iterative partial triple excitations (EOM-CC3) was tested. Coupled-cluster wave functions with single and double excitations (CCSD), single, double, and iterative partial triple excitations (CC3), and single, double, and full triple excitations (CCSDT) are employed with the correlation-consistent quadruple and quintuple zeta basis sets. The linear HBO X (1)sigma(+) state is predicted to lie 48.3 kcal mol(-1) (2.09 eV) lower in energy than the BOH X (1)sigma(+) linear stationary point at the CCSDT level of theory. The CCSDT BOH barrier to linearity is predicted to lie 3.7 kcal mol(-1) (0.16 eV). With a harmonic zero-point vibrational energy correction, the HBO X (1)sigma(+)-BOH X (1)A(') energy difference is 45.2 kcal mol(-1) (1.96 eV). The lowest triplet excited electronic state of HBO, a (3)A('), has a predicted excitation energy (T(e)) of 115 kcal mol(-1) (4.97 eV) from the HBO ground state minimum, while the lowest-bound BOH excited electronic state, b (3)A('), has a T(e) of 70.2 kcal mol(-1) (3.04 eV) with respect to BOH X (1)A('). The T(e) values predicted for the lowest singlet excited states are A (1)A(")<--X (1)sigma(+)=139 kcal mol(-1) (6.01 eV) for HBO and A (1)A(')<--X (1)A(')=102 kcal mol(-1) (4.42 eV) for BOH. Also for BOH, the triplet vertical transition energies are b (3)A(')<--X (1)A(')=71.4 kcal mol(-1) (3.10 eV) and c (3)A(")<--X (1)A(')=87.2 kcal mol(-1) (3.78 eV).     

The d (3)Pi(g)-c (3)Sigma(u) (+) band system of C2

A two-dimensional fluorescence (excitation/emission) spectrum of C2 produced in an acetylene discharge was used to identify and separate emission bands from the d (3)Pi(g)<--c (3)Sigma(u) (+) and d (3)Pi(g)<--a (3)Pi(u) excitations. Rotationally resolved excitation spectra of the (4<--1), (5<--1), (5<--2), and (7<--3) bands in the d (3)Pi(g)<--c (3)Sigma(u) (+) system of C2 were observed by laser-induced fluorescence spectroscopy. The molecular constants of each vibrational level, determined from rotational analysis, were used to calculate the spectroscopic constants of the c (3)Sigma(u) (+) state. The principal molecular constants for the c (3)Sigma(u) (+) state are B(e)=1.9319(19) cm(-1), alpha(e)=0.018 55(69) cm(-1), omega(e)=2061.9 cm(-1), omega(e)x(e)=14.84 cm(-1), and T(0)(c-a)=8662.925(3) cm(-1). We report also the first experimental observations of dispersed fluorescence from the d (3)Pi(g) state to the c (3)Sigma(u) (+) state, namely, d (3)Pi(g)(v=3)-->c (3)Sigma(u) (+)(v=0,1).     

The lowest singlet-triplet excitation energy of BN: A converged coupled cluster perspective

The notoriously small X 3Pi-a 1Sigma+ excitation energy of the BN diatomic has been calculated using high-order coupled cluster methods. Convergence has been established in both the one-particle basis set and the coupled cluster expansion. Explicit inclusion of connected quadruple excitations T4 is required for even semiquantitative agreement with the limit value, while connected quintuple excitations T5 still have an effect of about 60 cm(-1). Still higher excitations only account for about 10 cm(-1). Inclusion of inner-shell correlation further reduces Te by about 60 cm(-1) at the CCSDT, and 85 cm(-1) at the CCSDTQ level. Our best estimate, Te = 183+/-40 cm(-1), is in excellent agreement with earlier calculations and experiment, albeit with a smaller (and conservative) uncertainty. The dissociation energy of BN(X 3Pi) is De = 105.74+/-0.16 kcal/mol and D0 = 103.57+/-0.16 kcal/mol.     

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.     

Renner-Teller effect in C2H2+(X2Pi u) studied by rotationally resolved zero kinetic energy photoelectron spectroscopy

The Renner-Teller effect in C(2)H(2)(+)(X(2)Pi(u)) has been studied by using zero kinetic energy (ZEKE) photoelectron spectroscopy and coherent extreme ultraviolet (XUV) radiation. The rotationally resolved vibronic spectra have been recorded for energies up to 2000 cm(-1) above the ground vibrational state. The C triple bond C symmetric stretching (upsilon(2)), the CCH trans bending (upsilon(4)), and the CCH cis bending (upsilon(5)) vibrational excitations have been observed. The assigned vibronic bands are 4(1)(1)(kappa(2)Sigma(u)(+))(hot band), 4(1)(0)(mu/kappa(2)Sigma (u)(-/+)), 5(1)(0)(mu/kappa(2)Sigma (g)(+/-)), and 4(2)(0)(mu(2)Pi(u)), 4(2)(0)(kappa(2)Pi(u)), 4(1)(0)5(1)(0) (mu(2)Pi(g)), 0(0)(0)(X(2)Pi(u)), and 2(1)(0)(X(2)Pi(u)). The Renner-Teller parameters, the harmonic frequencies, the spin-orbit coupling constants, and the rotational constants for the corresponding vibronic bands have been determined by fitting the spectra with energy eigenvalues from the Hamiltonian that considers simultaneously Renner-Teller coupling, vibrational energies, rotational energies, and spin-orbit coupling interaction.     

The Beryllium tetramer: profiling an elusive molecule

The structure and energetics of Be(4) are investigated using state-of-the-art coupled-cluster methods. We compute the optimized bond length, dissociation energy, and anharmonic vibrational frequencies. A composite approach is employed, starting from coupled-cluster theory with single, double, and perturbative triple excitations extrapolated to the complete basis set (CBS) limit using Dunning's correlation consistent cc-pCVQZ and cc-pCV5Z basis sets. A correction for full triple and connected quadruple excitations in the smaller cc-pCVDZ basis set is then added, yielding an approximation to CCSDT(Q)/CBS denoted c～CCSDT(Q). Corrections are included for relativistic and non-Born-Oppenheimer effects. We obtain D(e) = 89.7 kcal mol(-1), D(0) = 84.9 kcal mol(-1), and r(e) = 2.043 A?. Second-order vibrational perturbation theory (VPT2) is applied to a full quartic force field computed at the c～CCSDT(Q) level of theory, yielding B(e) = 0.448 cm(-1) and fundamental frequencies of 666 (a(1)), 468 (e), and 571 (t(2)) cm(-1). Computations on the spectroscopically characterized Be(2) molecule are reported for the purpose of benchmarking our methods. Perturbative estimates of the effect of quadruple excitations are found to be essential to computing accurate parameters for Be(2); however, they seem to exert a much smaller influence on the structure and energetics of Be(4). Our extensive characterization of the Be(4) bonding potential energy surface should aid in the experimental identification of this thermodynamically viable but elusive molecule.     

Ab initio investigation of potential energy curves of the 23 electronic states of IBr correlating to neutral (2)P atoms

Potential energy surfaces for all Born-Oppenheimer electronic states of IBr molecule correlating to the neutral (2)P ((2)P(3/2) and (2)P(1/2)) iodine and bromine are calculated for the first time. Electric dipole and polarizability curves (static and transition) are also determined. Calculations include scalar and spin-orbit relativistic effects within all-electron Douglas-Kroll two-component Hamiltonian. Electron correlation is treated with quasi-degenerate multi-reference second-order perturbation theory. Seven adiabatic electronic states (X (1)Sigma(+), A'(3)Pi(2), A (3)Pi(1), 1 (3)Pi(0-), B (3)Pi(0+), B'(3)Sigma, and 2 (3)Pi(0+)) exhibit significant covalent bonding, and can support vibrational states. Calculated spectroscopic parameters agree with experiment to better than 1000 cm(-1) (T(e)), 10 cm(-1) (omega(e)), and 0.05 Angstrom (r(e)). A new 1 (3)Pi(0-) state correlating to ground-state atoms is predicted at T(e) approximately 14 000 cm(-1), omega(e) approximately 80 cm(-1), and r(e) approximately 3.0 Angstrom. The second new state (2 (3)Pi(0+)) correlates to excited iodine atom, with T(e) approximately 20 000 cm(-1), omega(e) approximately 115 cm(-1), and r(e) approximately 3.3 Angstrom. Non-adiabatic coupling parameters are calculated for the four avoided crossings, which arise due to electronic spin-orbit interaction. Estimated parameters of the B (3)Pi(0+)/B'(3)Sigma crossing (R(c) approximately 3.32 Angstrom; V approximately 120 cm(-1)) agree with experimental values. The previously unsuspected 2 (3)Pi(0-)/1 (1)Sigma(-) crossing of two repulsive surfaces provides a new collisional deactivation channel for Br* atoms at relative velocities above 1000 m s(-1). Several repulsive states (including 1 (1)Pi(1) and 2 (3)Pi(1)) intersect the B/B' system near the avoided crossing point, and may affect dynamics of IBr in strong laser fields.     

In search of definitive signatures of the elusive NCCO radical

Previous experimental assignments of the fundamental vibrational frequencies of NCCO have been brought into question by subsequent unsuccessful attempts to observe IR signatures of this radical at these frequencies. Here we compute the fundamental vibrational frequencies by applying second-order vibrational perturbation theory to the complete quartic force field computed at the all-electron (AE) coupled cluster singles, doubles, and perturbative triples level [CCSD(T)] with the correlation-consistent, polarized core-valence quadruple-zeta (cc-pCVQZ) basis set, which has tight functions to correctly describe core correlation. The AE-CCSD(T)/cc-pCVQZ geometric parameters are r(e)(N-C)=1.1623 A, r(e)(C-C)=1.4370 A, r(e)(C-O)=1.1758 A, theta(e)(N-C-C)=168.55 degrees , and theta(e)(C-C-O)=132.22 degrees . Our CCSD(T)/cc-pCVQZ values of the characteristic stretching frequencies nu(1) and nu(2) are 2171 and 1898 cm(-1), respectively, in stark contrast to the experimentally derived values of 2093 and 1774 cm(-1). Finally, focal-point extrapolations using correlation-consistent basis sets cc-pVXZ (X=D,T,Q,5,6) and electron correlation treatments as extensive as full coupled cluster singles, doubles, and triples (CCSDT) with perturbative accounting of quadruple excitations [CCSDT(Q)] determine the vibrationless barrier to linearity of NCCO and the dissociation energy (D(0)) of NCCO-->NC+CO to be 8.4 and 26.5 kcal mol(-1), respectively. Using our precisely determined dissociation energy, we recommend a new 0 K enthalpy of formation for NCCO of 50.9+/-0.3 kcal mol(-1).     

Born-Oppenheimer symmetry breaking in the C state of NO2: importance of static and dynamic correlation effects

A systematic theoretical treatment is performed with highly correlated ab initio theoretical methods to establish the structural nature of the C state of NO2. We predict the C state to have an asymmetric structure (point group C(s)). Spin-restricted and spin-unrestricted configuration interaction (CISD), coupled cluster [CCSD and CCSD(T)], multireference complete active space self-consistent field (CASSCF), and internally contracted multireference configuration interaction (ICMRCI) methods were used in conjunction with very large correlation-consistent polarized valence zeta cc-pVXZ and aug-cc-pVXZ [X = T, Q, 5] basis sets. The asymmetric C 2A' state is predicted to lie T(e) = 47.5 kcal/mol (2.06 eV, 16,600 cm(-1)) above the X 2A1 state at the aug-cc-pV5Z/UCCSD(T) level of theory, with T0 = 46.0 kcal/mol (2.00 eV, 16,100 cm(-1)), in good agreement with the experimental values of 46.77 kcal/mol (2.028 eV, 16,360 cm(-1)) by Weaver and 46.42 kcal/mol (2.013 eV, 16,234 cm(-1)) by Aoki. The symmetric structure (in C(2v) symmetry) with re(NO) = 1.274 A and theta(e) (ONO) = 109.9 degrees is a transition state between the two equivalent asymmetric (in C(s) symmetry) structures and is located only 1.53 kcal/mol (0.066 eV, 540 cm(-1)) above the asymmetric structure. The asymmetric structure is predicted to have structural parameters r(e)(NOl) = 1.489 A, r(e)(NOs) = 1.169 A, and theta(e)(ONO) = 109.7 degrees with the same method, aug-cc-pV5Z/UCCSD(T). The averaged NO bond distance is 1.329 A, and the difference between the two NO bond distances is 0.320 A. The three harmonic vibrational frequencies for the C 2A' state are 1656 (in-phase stretch), 759 (bend), and 378 (out- of-phase stretch) cm(-1). While these theoretical results further corroborate the previous predictions concerning the asymmetric nature of the C state, there remains discrepancy between the theoretical and experimental symmetric stretching mode omega1 (1656 and 923 cm(-1), respectively). It is possible, however, that this disagreement could be resolved by a reassignment of the corresponding lines in the experimental spectrum, though additional vibronic simulations of the spectrum are required to confirm this proposition.     

1 Pi<--X1 Sigma+ band systems of jet-cooled ScCo and YCo

Rotationally resolved resonant two-photon ionization (R2PI) spectra of ScCo and YCo are reported. The measured spectra reveal that these molecules possess ground electronic states of (1)Sigma(+) symmetry, as previously found in the isoelectronic Cr(2) and CrMo molecules. The ground state rotational constants for ScCo and YCo are B(0)(")=0.201 31(22) cm(-1) and B(0) (")=0.120 96(10) cm(-1), corresponding to ground state bond lengths of r(0) (")=1.812 1(10) A and r(0) (")=1.983 0(8) A, respectively. A single electronic band system, assigned as a (1)Pi<--X (1)Sigma(+) transition, has been identified in both molecules. In ScCo, the (1)Pi state is characterized by T(0)=15,428.8, omega(e)(')=246.7, and omega(e)(')x(e)(')=0.73 cm(-1). In YCo, the (1)Pi state has T(0)=13 951.3, omega(e)(')=231.3, and omega(e)(')x(e) (')=2.27 cm(-1). For YCo, hot bands originating from levels up to v(")=3 are observed, allowing the ground state vibrational constants omega(e)(")=369.8, omega(e)(")x(e)(")=1.47, and Delta G(12)(")=365.7 cm(-1) to be deduced. The bond energy of ScCo has been measured as 2.45 eV from the onset of predissociation in a congested vibronic spectrum. A comparison of the chemical bonding in these molecules to related molecules is presented.     

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.     

Vibronic spectroscopy of unsaturated transition metal complexes: CrC2H, CrCH3, and NiCH3

Vibronically resolved resonant two-photon ionization and dispersed fluorescence spectra of the organometallic radicals CrC(2)H, CrCH(3), and NiCH(3) are reported in the visible and near-infrared wavelength regions. For CrC(2)H, a complicated vibronic spectrum is found in the 11 100-13 300 cm(-1) region, with a prominent vibrational progression having omega(e) (')=426.52+/-0.84 cm(-1), omega(e) (')x(e) (')=0.74+/-0.13 cm(-1). Dispersed fluorescence reveals a v(")=1 level of the ground state with DeltaG(1/2) (")=470+/-20 cm(-1). These vibrational frequencies undoubtedly pertain to the Cr-C(2)H stretching mode. It is suggested that the spectrum corresponds to the A (6)Sigma(+)<--X (6)Sigma(+) band system, with the CrC(2)H molecule being linear in both the ground and the excited state. The related CrCH(3) molecule displays a vibronic spectrum in the 11 500-14 000 cm(-1) region. The upper state of this system displays six sub-bands that are too closely spaced to be vibrational structure, but too widely separated to be K structure. It is suggested that the observed spectrum is a (6)E<--X (6)A(1) band system, analogous to the well-known B (6)Pi<--X (6)Sigma(+) band systems of CrF and CrCl. The ground state Cr-CH(3) vibration is characterized by omega(e) (")=525+/-17 cm(-1) and omega(e) (")x(e) (")=7.9+/-6 cm(-1). The spectrum of NiCH(3) lies in the 16 100-17 400 cm(-1) range and has omega(e) (')=455.3+/-0.1 cm(-1) and omega(e) (')x(e) (')=6.60+/-0.03 cm(-1). Dispersed fluorescence studies provide ground state vibrational constants of omega(e) (")=565.8+/-1.6 cm(-1) and omega(e) (")x(e) (")=1.7+/-3.0 cm(-1). Again, these values correspond to the Ni-CH(3) stretching motion. (c) 2004 American Institute of Physics.     

The ground and two lowest-lying singlet excited electronic states of copper hydroxide (CuOH)

Various ab initio methods, including self-consistent field (SCF), configuration interaction, coupled cluster (CC), and complete-active-space SCF (CASSCF), have been employed to study the electronic structure of copper hydroxide (CuOH). Geometries, total energies, dipole moments, harmonic vibrational frequencies, and zero-point vibrational energies are reported for the linear 1Sigma+ and 1Pi stationary points, and for the bent ground-state X 1A', and excited-states 2 1A' and 1 1A". Six different basis sets have been used in the study, Wachters/DZP being the smallest and QZVPP being the largest. The ground- and excited-state bending modes present imaginary frequencies for the linear stationary points, indicating that bent structures are more favorable. The effects of relativity for CuOH are important and have been considered using the Douglas-Kroll approach with cc-pVTZ/cc-pVTZ_DK and cc-pVQZ/cc-pVQZ_DK basis sets. The bent ground and two lowest-lying singlet excited states of the CuOH molecule are indeed energetically more stable than the corresponding linear structures. The optimized geometrical parameters for the X 1A' and 1 1A" states agree fairly well with available experimental values. However, the 2 1A' structure and rotational constants are in poor agreement with experiment, and we suggest that the latter are in error. The predicted adiabatic excitation energies are also inconsistent with the experimental values of 45.5 kcal mol(-1) for the 2 1A' state and 52.6 kcal mol(-1) for the 1 1A" state. The theoretical CC and CASSCF methods show lower adiabatic excitation energies for the 1 1A" state (53.1 kcal mol(-1)) than those for the corresponding 2 1A' state (57.6 kcal mol(-1)), suggesting that the 1 1A" state might be the first singlet excited state while the 2 1A' state might be the second singlet excited state.     

1-Germavinylidene (Ge═CH2), germyne (HGeCH), and 2-germavinylidene (H2Ge═C) molecules and isomerization reactions among them: anharmonic rovibrational analyses

Theoretical investigations of three equilibrium structures and two associated isomerization reactions of the GeCH(2) - HGeCH - H(2)GeC system have been systematically carried out. This research employed ab initio self-consistent-field (SCF), coupled cluster (CC) with single and double excitations (CCSD), and CCSD with perturbative triple excitations [CCSD(T)] wave functions and a wide variety of correlation-consistent polarized valence cc-pVXZ and cc-pVXZ-DK (where X = D, T, Q) basis sets. For each structure, the total energy, geometry, dipole moment, harmonic vibrational frequencies, and infrared intensities are predicted. Complete active space SCF (CASSCF) wave functions are used to analyze the effects of correlation on physical properties and energetics. For each of the equilibrium structures, vibrational second-order perturbation theory (VPT2) has been utilized to obtain the zero-point vibration corrected rotational constants, centrifugal distortion constants, and fundamental vibrational frequencies. The predicted rotational constants and anharmonic vibrational frequencies for 1-germavinylidene are in good agreement with available experimental observations. Extensive focal point analyses, including CCSDT and CCSDT(Q) energies and basis sets up to quintuple zeta, are used to obtain complete basis set (CBS) limit energies. At all levels of theory employed in this study, the global minimum of the GeCH(2) potential energy surface (PES) is confirmed to be 1-germavinylidene (GeCH(2), 1). The second isomer, germyne (HGeCH, 2) is predicted to lie 40.4(41.1) ± 0.3 kcal mol(-1) above the global minimum, while the third isomer, 2-germavinylidene (H(2)GeC, 3) is located 92.3(92.7) ± 0.3 kcal mol(-1) above the global minimum; the values in parentheses indicate core-valence and zero-point vibration energy (ZPVE) corrected energy differences. The barriers for the forward (1→2) and reverse (2→1) isomerization reactions between isomers 1 and 2 are 48.3(47.7) ± 0.3 kcal mol(-1) and 7.9(6.6) ± 0.3 kcal mol(-1), respectively. On the other hand, the barriers of the forward (2→3) and reverse (3→2) isomerization reactions between isomers 2 and 3 are predicted to be 55.2(53.2) ± 0.3 kcal mol(-1) and 3.3(1.6) ± 0.3 kcal mol(-1), respectively.     

The low-lying electronic excited states of NiCO

Highly correlated coupled cluster methods with single and double excitations (CSSD) and CCSD with perturbative triple excitations were used to predict molecular structures and harmonic vibrational frequencies for the electronic ground state X 1Sigma+, and for the 3Delta, 3Sigma+, 3Phi, 1 3Pi, 2 3Pi, 1Sigma+, 1Delta, and 1Pi excited states of NiCO. The X 1Sigma+ ground state's geometry is for the first time compared with the recently determined experimental structure. The adiabatic excitation energies, vertical excitation energies, and dissociation energies of these excited states are predicted. The importance of pi and sigma bonding for the Ni-C bond is discussed based on the structures of excited states.     

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.     

Low-lying electronic states of FeNC and FeCN: a theoretical journey into isomerization and quartet/sextet competition

With several levels of multireference and restricted open-shell single-reference electronic structure theory, optimum structures, relative energetics, and spectroscopic properties of the low-lying (6)Delta, (6)Pi, (4)Delta, (4)Pi, and (4)Sigma(-) states of linear FeNC and FeCN have been investigated using five contracted Gaussian basis sets ranging from Fe[10s8p3d], C/N[4s2p1d] to Fe[6s8p6d3f2g1h], C/N[6s5p4d3f2g]. Based on multireference configuration interaction (MRCISD+Q) results with a correlation-consistent polarized valence quadruple-zeta (cc-pVQZ) basis set, appended with core correlation and relativistic corrections, we propose the relative energies: T(e)(FeNC), (6)Delta(0)<(6)Pi (2300 cm(-1))<(4)Delta (2700 cm(-1))<(4)Pi (4200 cm(-1))<(4)Sigma(-); and T(e)(FeCN), (6)Delta(0)<(6)Pi (1800 cm(-1))<(4)Delta (2500 cm(-1))<(4)Pi (2900 cm(-1))<(4)Sigma(-). The (4)Delta and (4)Pi states have massive multireference character, arising mostly from 11sigma-->12sigma promotions, whereas the sextet states are dominated by single electronic configurations. The single-reference CCSDT-3 (coupled cluster singles and doubles with iterative partial triples) method appears to significantly overshoot the stabilization of the quartet states provided by both static and dynamical correlation. The (4,6)Delta and (4,6)Pi states of both isomers are rather ionic, and all have dipole moments near 5 D. On the ground (6)Delta surface, FeNC is predicted to lie 0.6 kcal mol(-1) below FeCN, and the classical barrier for isocyanide/cyanide isomerization is about 6.5 kcal mol(-1). Our data support the recent spectroscopic characterization by Lei and Dagdigian [J. Chem. Phys. 114, 2137 (2000)] of linear (6)Delta FeNC as the first experimentally observed transition-metal monoisocyanide. Their assignments for the ground term symbol, isotopomeric rotational constants, and the Fe-N omega(3) stretching frequency are confirmed; however, we find rather different structural parameters for (6)Delta FeNC:r(e)(Fe-N)=1.940 A and r(N-C)=1.182 A at the cc-pVQZ MRCISD+Q level. Our results also reveal that the observed band of FeNC originating at 27 236 cm(-1) should have an analog in FeCN near 23 800 cm(-1) of almost equal intensity. Therefore, both thermodynamic stability and absorption intensity factors favor the eventual observation of FeCN via a (6)Pi<--(6)Delta transition in the near-UV.