Search

The 1Sigma(+), 3Sigma(+), 1Pi, and 3Pi states of the negative ions and the 2Pi and 2Sigma(+) states of the neutral alkali oxides are studied at high levels of theory. The calculations show that ground state of the negative ions changes from 3Pi for LiO(-) to 1Sigma(+) for KO(-). Although the calculations give a 3Pi ground state for NaO(-), we cannot rule out the possibility that the very low-lying 1Sigma(+) state is the true ground state. The Franck-Condon factors for photodetachment of an electron from the 1Sigma(+) and/or 3Pi states of the negative ion are presented to help interpret photodetachment experiments. Our best results for the A 2Sigma(+) - X 2Pi separations in LiO and NaO are 2496 and 2061/cm, which are in excellent agreement with that deduced (2516 and 2018/cm) from experiment.

The bonding in the X 1Sigma(+)g state of Mg2 is investigated using near-complete valence one-particle Slater and Gaussian basis sets containing up to h functions. It is shown that the four-electron complete CI limit can be approached using a sequence of either second-order CI (SOCI) or interacting correlated fragment (ICF) calculations. At the valence level, the best estimate of the dissociation energy D(e) was 464/cm. This is a lower limit and is probably within 5/cm of the complete basis value.

Spectroscopic results were determined for the ground and low-lying states of the oxides and sulfides of Mo and Tc, using the single-reference-based modified coupled pair functional method of Ahlrichs et al. (1985) and Chong et al. (1986) and the multireference-based state-averaged CASSCF/MRCI method. Spectroscopic constants, dipole moments, Mulliken populations, and radiative lifetimes were calculated for selected low-lying states of these molecular systems. The spectroscopy of the MoS and TcS molecules was found to be quite analogous to the corresponding oxides.

The low-lying (up to about 40,000/cm) doublet and quartet states of AlC have been studied at the CASSCF and MRCI levels of theory. The bonding in the X4Sigma(-) state is shown to involve three one-electron bonds. The first excited state, a 2Pi, lies about 8700/cm higher and like the X4Sigma(-) state is derived from ground state Al and C, but has a two-electron pi bond leading to a much shorter bond length.

The equilibrium bond distance in methane has been optimized using coupled-pair functional and contracted CI wave functions, and a Gaussian basis that includes g-type functions on carbon and d-type functions on hydrogen. The resulting bond distance, when corrected for core-valance correlation effects, agrees with the experimental value of 2.052 a(0) to within the experimental uncertainty of 0.002 a(0). The main source of error in the best previous studies, which showed discrepancies with experiment of 0.007 a(0) is shown to be basis set incompleteness. In particular, it is important that the basis set be close to saturation, at least for the lower angular quantum numbers.

Theoretical dissociation energies for the first- and second-row transition metal hydride positive ions are critically compared against recent experimental values obtained from ion beam reactive scattering methods. Theoretical spectroscopic parameters and dipole moments are presented for the ground and several low-lying excited states. The calculations employ large Gaussian basis sets and account for electron correlation using the single-reference single- and double-excitation configuration interaction and coupled-pair-functional methods. The Darwin and mass-velocity contributions to the relativistic energy are included in the all-electron calculations on the first-row systems using first-order perturbation theory, and in the second-row systems using the Hay and Wadt relativistic effective core potentials. The theoretical D(0) values for the second-row transition metal hydride positive ions should provide a critical measure of the experimental values, which are not as refined as many of those in the first transition row.

A detailed study is presented of the basis set requirements and effects of correlation on the geometry and structure of small Al(n) clusters n = 2, 4, and 13. An effective core potential (ECP) is developed from the Huzinaga basis which accurately reproduces the all-electron results. It is found that depolarization functions are very important in computing the bond length, and that the bond contraction obtained is about as large for a 13-atom cluster as for the Al(4) rhombus structure. With d functions on each center in Al(13), a bond distance shorter than the bulk is obtained, as expected. The inclusion of correlation in Al(4) is found to be less important than the addition of polarization functions for determining the bond length, but increases substantially the calculated atomization energy. These calibration calculations allow an accurate treatment for Al(n) clusters to be designed.

Qualitative calculations are performed to elucidate the nature of the side-on interaction of both N2 and CO with a single Fe atom. The systems are found to be quite similar, with bonding leading to an increase in the CO or N2 bond length and a decrease in the vibrational frequency. The CO or N2 stretching modes lead to a large dipole derivative along the metal-ligand bond axis. The populations show an almost identical, large donation from the Fe 3d orbitals into the CO or N2 Pi-asterisk. The larger system Fe2N2 is then considered, with the N2 bridging the Fe2, both parallel and perpendicular to the Fe2 bond axis for two different Fe-Fe distances. For FeN2, the shift in the observed N2 frequency is smaller than observed for the alpha state of N2/Fe(111). The shift in the N2 vibrational frequency increases when the N2 interacts with two Fe atoms, either at the Fe-Fe nearest neighbor distance or at the first layer Fe-Fe distance, when the side-on N2 axis is oriented perpendicular to an Fe-Fe bond.

The structure and stability of aluminum clusters containing up to six atoms have been studied using correlated wave functions and extended basis sets. The lowest energy structure is planar for Al4 and Al5, but three dimensional for Al6. The icosahedral, hcp, fcc, and two planar structures of Al13 were considered at the SCF level. The lowest energy structure is the icosahedron, but the planar structures are fairly low lying even in this case. A simplified description using two- and three-body interaction potentials is found to agree well with the ab initio structures and binding energies.

An improved effective core potential (ECP) technique is described and used to give ECP parameters for the atoms of the first two rows of the periodic table. A given basis set is parametrized which allows for a direct comparison with all-electron calculations. Extensive test calculations on first- and second-row molecules using the ECP have been performed, giving excellent agreement with the all-electron results at the SCF, CASSCF, and CI levels of treatment. Correlating and diffuse functions may be added without modifying the ECP parameters. The present ECP descriptions result in CPU time reductions of the order of 50 percent in addition to reduced disk storage.

The equilibrium bond distance in methane was optimized using coupled-pair functional and contracted CI wave functions, and a Gaussian basis that includes g-type functions on carbon and d-type functions on hydrogen. The resulting bond distance, when corrected for core-valence correlation effects, agrees with the experimental value of 2.052 a(0) to within the experimental uncertainty of 0.002 a(0). The main source of error in the best previous studies, which showed discrepancies with experiment of 0.007 a(0) is shown to be basis set incompleteness. In particular, it is important that the basis set be close to saturation, at least for the lower angular quantum numbers.

A systematic analysis of the low-lying states of all of the second-row transition metal (TM) hydrides except CdH is reported. The calculations included the dominant relativistic contributions through the use of the relativistic effective core potentials of Hay and Wadt (1985). Electron correlation was incorporated, using single-plus-double configuration interaction, the coupled pair functional (CPF) formalism of Ahlrichs et al. (1985), and the Chong and Langhoff (1986) modified version of the CPF method. The spectroscopic parameters D(e), r(e), and mu(e) determined for the low-lying states are compared with the available experimental data and previous theoretical results. In contrast to the first-row TM hydrides studied earlier (Chong et al., 1986), the spectroscopic constants for the second-row TM hydrides were found to be much less sensitive to the level of correlation treatment.

Reliable experimental diople moments are available for the ground states of SeH and BrH whereas no values have been reported for GaH and AsH a recently reported experimental dipole moment for GeH of 1.24 + or -0.01 D has been seriously questioned, and a much lower value of, 0.1 + or - 0.05 D, suggested. In this work, we report accurate theoretical dipole moments, dipole derivatives, dissociation energies, and spectroscopic constants (tau(sub e), omega(sub e)) for the ground states of GaH through BrH.

The potential energy function and electric dipole moment function (EDMF) are computed for CIO Chi(sup 2)Pi using several different techniques to include electron correlation. The EDMF is used to compute Einstein coefficients, vibrational lifetimes, and dipole moments in higher vibrational levels. Remaining questions concerning the position of the maximum of the EDMF may be resolved through experimental measurement of dipole moments of higher vibrational levels. The band strength of the 1-0 fundamental transition is computed to be 12 +/- 2 /sq cm atm in good agreement with three experimental values, but larger than a recent value of 5 /sq cm atm determined from infrared heterodyne spectroscopy. The theoretical methods used include SCF, CASSCF, multireference singles plus doubles configuration interaction (MRCI) and contracted CI, coupled pair functional (CPF), and a modified version of the CPF method. The results obtained using the different methods are critically compared.

Face-centered cubic (fcc) and hexagonally close-packed (hcp) structures are compared for Be(13) and Be(55) clusters. Both Be(13) and Be(55) prefer the fcc structure over the bulk hcp structure, but the energy difference per atom decreases for Be(55) relative to Be(13). The binding energy per atom, 1.3 eV for Be(55) and 0.8-0.9 eV for Be(13), reflects the greater total number of bonds in the larger cluster rather than a difference in bonding. The energies per bond are much more similar, in the range of 0.30-0.34 eV for both clusters. The size of the p-basis set used influences both stability and ionization potentials strongly.

Accurate theoretical dipole moments (mu(sub c) have been computed for the X(exp 2)Pi ground states of Si(-)H(+)(0.118 D), Ge(+)H(-)(0.085 D) and Sn(+)H(-)(0.357 D). The trend down the periodic table is regular and follows that expected from the electronegativities of the group IV atoms. The dipole moment of 1.24 +/- 0.1 D for GeH recently derived by Brown, Evenson and Sears from the relative intensities of electric and magnetic dipole transitions in the 10 microns spectrum of the X(exp 2)Pi state is seriously questioned.

Recently, the TR frequencies for NH3 molecularly adsorbed on Cu(100) have been measured by infrared laser-induced photodesorption using a tunable laser operating in the group frequency range. Additionally, photodesorption, photo-fragmentation, and photoionization has been observed using a CO2 laser in the 10 microns range. In this contribution we report on work in progress where we use quantum chemical methods to analyze the shifts in the NH, frequencies relative to the gas phase upon chemisorption. We also compute the Cu- N stretch vibration, which has been observed experimentally,and the change in the IR intensities associated with chemisorption; this change is insufficient by itself to explain the increased photoionization efficiency observed.

An accurate effective core potential (ECP), including frozen 3s, 3p orbitals and a single-zeta contracted 3d orbital, has been developed for germanium. The ECP with associated valence basis set reproduces the corresponding all-electron results for the atomic excitations and the geometry and excitation energies of GeH2. At the SCF, CAS SCF and CI levels the maximum difference from the all-electron results is 0.5 kcal/mol in the (sup 1)A(sub 1) - (sup 3)B(sub 1) excitation energy. Finally, the ECP description is used with an extended valence basis set and large-scale CAS SCF and multi-reference CI wavefunctions to compute the singlet-triplet separation; the final CI result including the Davidson correction is 22.5 kcal/mol.