Engineering topics
Publications and source records attributed to Langhoff, Stephen R..
The D(sub 0) values for AlH2 and AlAr are computed using the coupled cluster approach in conjunction with large basis sets. Basis set superposition and spin-orbit effects are accounted for as they are sizeable due to the small binding energy. The computed dissociation energy for AlAr is 101 /cm , which is 83% of the experimental value (122.4/ cm). Our best estimate for the H2 binding energy in AlH2 is 40 +/- 28 /cm.
Multireference configuration-interaction (MRCI) wavefunctions and potential energy curves have been calculated for the X(sup 3)Sigma(sub g,sup -), B(sup 3)Sigma(sub u, Sup -) and B"(sup 3)Pi((sub u) states of S(sub 2) using correlation consistent Gaussian basis sets. These wavefunctions are utilized to compute the the transition dipole moments of the B(sup 3)Sigma(sub g, sup -) - X(sup 3) Sigma(sub g, sup -) and B"(sup 3)Pi(sub u) - X(sup 3)Sigma(sub g, sup -) systems. Oscillator strengths, transition probabilities, and radiative lifetimes are computed for the X-B system and comparison is made with experimental data.
In 1986, we reported a dissociation energy (D(sub 0) of 2.31 eV for the X(sup 2)Pi ground state of MgO(+). This value was determined by computing the dissociation energy to the Mg(2+) + O(-) limit and adjusting the value to the Mg(+) + O limit using the experimental Ionization Potential (IP) of Mg(+) and the Electron Affinity (EA) of O. The success of this method relies on the assumption that there is little covalent contribution to the bonding. The very small (0.04 eV) correlation contribution to the binding energy was taken as corroboration for the validity of this approach. Our earlier theoretical value was estimated to be accurate to at least 0.2 eV. It is in excellent agreement with the subsequent value of 2.30 +/- 0.13 eV determined by Freiser and co-workers from photodissociation experiments. It is also consistent with the upper (less than 3.1 eV) and lower (greater than 1.1 eV) bounds determined by Rowe obtained by studying the reactions of Mg(+) with 03 and NO2. However, it is inconsistent with an upper bound of 1.7 eV reported by Kappes and Staley based on their failure to observe MgO(+) in the reaction of Mg(+) with N2O. The picture became somewhat clouded, however, by the recent guided-ion beam mass spectrometric studies of Dalleska and Armentrout. Their initial analysis of the reaction data for Mg(+) + O2 lead to a bond dissociation energy of 2.92 +/- 0.25 eV, which is considerably larger than the value of 2.47 +/- 0.06 eV deduced from their studies of the Mg(+)+NO2 reaction.
The singlet and triplet states of AlCu and AlAg below about 32 000/cm are studied using the internally contracted multireference configuration-interaction method. A more elaborate study of the X(sup 1)Sum(sup +) ground state of AlCu is undertaken using extended Gaussian basis sets, including the effect of inner-shell correlation and including a perturbational estimate of relativistic effects. Our best estimate of the spectroscopic constants (r(sub 0), DeltaG(sub 1/2), and D(sub 0)) for the X(sup 1)Sum(sup+) state with the experimental values in parentheses are: 4.416(4.420) a(sub 0), 295 (294) /cm, and 2.318 (2.315) eV. The calculations definitively assign the upper state in the observed transition at 14 892/cm to the lowest (sup 1)Prod state. The calculated spectroscopic constants and radiative lifetime for the (sup 1)Prod state are in good agreement with experiment. The calculations support the tentative assignments of Behm et al. for three band systems observed in the visible region between 25 000 and 28 000 / cm. However, the computed spectroscopic constants are in very poor agreement with those deduced from an analysis of the spectra. Analogous theoretical results for AlAg suggest that the (2)(sup 3)Prod, (3)(sup 3)Prod, and (3)(sup 1)Sum(sup +) states account for the bands observed, but not assigned, by Duncan and co-workers.
The dipole moment function for the X(sup 2)II state of NO is studied as a function of the completeness in both the one- and n-particle spaces. Einstein coefficients are presented that are significantly more accurate than previous tabulations for the higher vibrational levels. The theoretical values give considerable insight into the limitations of recently published ratios of Einstein coefficients measured by spectrally resolved infrared chemiluminescence.
Due to recent advances in methods and computers, the accuracy of ab calculations has reached a point where these methods can be used to provide accurate spectroscopic constants for small molecules; this will be illustrated with several examples. We will show how ab initio calculations where used to identify the Hermann infrared system in N2 and two band systems in CO. The identification of all three of these band systems relied on very accurate calculations of quintet states. The analysis of the infrared spectra of cool stars requires knowledge of the intensity of vibrational transitions in SiO for high nu and J levels. While experiment can supply very accurate dipole moments for nu = 0 to 3, this is insufficient to construct a global dipole moment function. We show how theory, combined by the experiment, can be used to generate the line intensities up to nu = 40 and J = 250. The spectroscopy of transition metal containing systems is very difficult for both theory and experiment. We will discuss the identification of the ground state of Ti2 and the spectroscopy of AlCu as examples of how theory can contribute to the understanding of these complex systems.
Density functional theory (DFT) is found to give a better description of the geometries and vibrational frequencies of FeL and FeL(sup +) systems than second order Moller Plesset perturbation theory (MP2). Namely, the DFT correctly predicts the shift in the CO vibrational frequency between free CO and the Sigma(sup -) state of FeCO and yields a good result for the Fe-C distance in the quartet states of FeCH4(+) 4 These are properties where the MP2 results are unsatisfactory. Thus DFT appears to be an excellent approach for optimizing the geometries and computing the zero-point energies of systems containing first transition row atoms. Because the DFT approach is biased in favor of the 3d(exp 7) occupation, whereas the more traditional approaches are biased in favor of the 3d(exp 6) occupation, differences are found in the relative ordering of states. It is shown that if the dissociation is computed to the most appropriate atomic asymptote and corrected to the ground state asymptote using the experimental separations, the DFT results are in good agreement with high levels of theory. The energetics at the DFT level are much superior to the MP2 and in most cases in good agreement with high levels of theory.
BH3CO, BF3CO, BH3NH3, and BF3NH3 have been studied using density functional theory (DFT). The computed geometries are in reasonable agreement with experiment. BF3CO is found to be weakly bound. The binding energies of BH3CO, BH3NH3 , and BF3NH3 have also been computed using the coupled cluster singles and doubles level of theory, including a perturbational estimate of the connected triples [CCSD(T)]. These CCSD(T) calculations are performed using correlation consistent polarized-valence triple-zeta basis sets. These are expected to be the most accurate binding energies to date. The DFT binding energies are in good agreement with those obtained at the CCSD(T) level. The bonding is analyzed using the constrained space orbital variation (CSOV) technique. In addition to analyzing the bonding in BH3CO, BH3NH3, and BF3NH3, the CSOV approach is used to study why the bonding in BF3CO is so weak.
Potential-energy curves, transition-moment functions, Einstein Alpha coefficients and radiative lifetimes associated with the Rydberg state of nitric oxide are calculated in large Gaussian basis sets employing internally contracted MRCI from state-averaged CASSCF wave functions. Computationally, a strong sensitivity of the results to the choice of the CASSCF active space is seen, and careful examination of this behavior gives insight into matters related to the treatment of weakly occupied orbitals, which is a consideration in any study Rydberg excited states. The epsilon system transition-moment function and radiative lifetimes are in disagreement with another calculation in the literature, but reasons for the discrepancies are examined, and the present results are determined to be more accurate.
Abstract Cr2 is studied at the coupled cluster singles and doubles level of theory, including a perturbational estimate of the connected triples [CCSD(T)]. Both spin restricted, RCCSD(T), and spin unrestricted, UCCSD(T), approaches are used. The UCCSD(T) bond length is much longer than the experimental value, while the RCCSD(T) bond length is shorter. This difference is due to the inability of the spin restricted Hartree-Fock reference to describe the antiferromagnetic coupling at larger r values. In contrast with the results obtained using density functional theory in conjunction with older functionals, the results obtained with recently developed functionals are in agreement with the UCCSD(T) results.
Electronic structure calculations can provide accurate spectroscopic data (such as molecular structures) vibrational frequencies, binding energies, etc.) that have been very useful in explaining trends in experimental data and in identifying incorrect experimental measurements. In addition, ab initio calculations. have given considerable insight into the many interactions that make the chemistry of transition metal systems so diverse. In this review we focus on cases where calculations and experiment have been used to solve interesting chemical problems involving metal ions. The examples include cases where theory was used to differentiate between disparate experimental values and cases where theory was used to explain unexpected experimental results.
This article provides an overview of the principal computational approaches and their accuracy for the study of electronic spectroscopy of diatomic molecules. We include a number of examples from our work that illustrate the range of application. We show how full configuration interaction benchmark calculations were instrumental in improving the understanding of the computational requirements for obtaining accurate results for diatomic spectroscopy. With this understanding it is now possible to compute radiative lifetimes accurate to within 10% for systems involving first- and second-row atoms. We consider the determination of the infrared vibrational transition probabilities for the ground states of SiO and NO, based on a globally accurate dipole moment function. We show how we were able to assign the a(sup "5)II state of CO as the upper state in the recently observed emission bands of CO in an Ar matrix. We next discuss the assignment of the photoelectron detachment spectra of NO and the alkali oxide negative ions. We then present several examples illustrating the state-of-the-art in determining radiative lifetimes for valence-valence and valence-Rydberg transitions. We next compare the molecular spectroscopy of the valence isoelectronic B2, Al2, and AlB molecules. The final examples consider systems involving transition metal atoms, which illustrate the difficulty in describing states with different numbers of d electrons.
The bond dissociation energies of 1-butene and several calibration systems are computed using the G2(MP2) approach. The agreement between the calibration systems and experiment is very good. The computed values for 1-butene are compared with calibration systems and the agreement between the computed results for 1-butene and the "rule of thumb" values from the smaller systems is remarkably good.
Large atomic natural orbital (ANO) basis sets are tabulated for the Sc to Cu. The primitive sets are taken from the large sets optimized by Partridge, namely (21s 13p 8d) for Sc and Ti and (20s 12p 9d) for V to Cu. These primitive sets are supplemented with three p, one d, six f, and four g functions. The ANO sets are derived from configuration interaction density matrices constructed as the average of the lowest states derived from the 3d(sup n)4s(sup 2) and 3d(sup n+1)4s(sup 1) occupations. For Ni, the 1S(3d(sup 10)) state is included in the averaging. The choice of basis sets for molecular calculations is discussed.
The molecular structure of cis-FONO has been determined with the CCSD(T) correlation method using an spdf quality basis set. In agreement with previous coupled-cluster calculations but in disagreement with density functional theory, cis-FONO is found to exhibit normal bond distances. The quadratic and cubic force fields of cis-FONO have also been determined in order to evaluate the effect of vibrational averaging on the molecular geometry. Vibrational averaging is found to increase bond distances, as expected, but it does not affect the qualitative nature of the bonding. The CCSD(T)/spdf harmonic frequencies of cis-FONO support our previous assertion that a band observed at 1200 /cm is a combination band (upsilon(sub 3) + upsilon(sub 4)), and not a fundamental.
The XCN and XNC (X=F, Cl) isomers have been investigated using the CCSD and CCSD(T) methods in conjunction with a TZ2P basis set. Equilibrium geometries, dipole moments, harmonic frequencies, IR intensities and relative energies have been evaluated. The CCSD(T) geometries and vibrational frequencies for the XCN isomers are in good agreement with the available experimental data. The CCSD(T) results for FCN and FNC are in good agreement with the CEPA calculations of Botshwina et al., with the exception of the energy difference, which the CEPA method underestimates by about 1.2 kcal/mol. FCN and CICN are shown to be lower in energy than the FNC and ClNC isomers by 69.511.0 and 42.711.0 kcal/mol (0 K), respectively.
Calculations performed using both density functional theory (DFT) and the modified coupled-pair functional (MCPF) approach show that FeC5H6(+) is more stable than HFeC5H5(+) by about 10 kcal/mol. The ground state of FeC5H6(+) is a quartet state derived from the 3d(sup 7) occupation of Fe+. For HFeC5H6(+), the MCPF 5 approach yields a sextet ground state while DFT yields a quartet; however, these two states are close in energy at both levels of theory.
The C-H bond energy in C2H6 is computed to be 99.76 +/- 0.35 kcal/mol, which is in excellent agreement with the most recent experimental values. The calculation of the C-H bond energy by direct dissociation and by an isodesmic reaction is discussed.