Engineering topics

Bauschlicher, C. W., Jr.

Publications and source records attributed to Bauschlicher, C. W., Jr..

At least 37 records · Page 2

NASA NTRS1986

Ab initio calculations on the positive ions of the alkaline-earth oxides, fluorides, and hydroxides

Theoretical dissociation energies are presented for the alkaline-earth fluoride, hydroxide, and oxide positive ions that are considered to be accurate to 0.1-0.2 eV. The r(e) for the positive ions are found to be consistently shorter than the corresponding neutrals by 0.07 + or -0.02 A. The bonding in the ground states is demonstrated to be of predominantly M + 2 X - character. The a 3 Pi and A 1 Pi are found to lie considerably above the X 1 Sigma + ground states of the alkaline-earth fluoride and hydroxide positive ions. The overall agreement of the theoretical ionization potentials with the available experimental appearance potentials is satisfactory; these values should represent the most accurate and consistent set available.

Partridge, H.↗

NASA NTRS1986

Theoretical study of the X 2Pi and A 2Sigma(+) states of CuO and CuS

Theoretical spectroscopic constants and dipole moments are determined for the X 2Pi and A 2Sigma(+) states of CuO and CuS, using extended Gaussian basis sets and incorporating correlation using both configuration interaction (CI) and coupled pair (CPF) methods. Relativistic corrections (Darwin plus mass velocity), included using first-order perturbation theory, are found to be relatively small. At the CPF level, significant configuration mixing with the 3d(9)4s(1) occupation occurs for both states, resulting in some bond shortening and a significant increase in the dissociation energy. The best spectroscopic parameters are in excellent agreement with experiment for both CuO and CuS. In analogy with CuO, CuS is predicted to have a yet unobserved 2Sigma(+) state near 10,000 per cm.

Langhoff, S. R.↗

NASA NTRS1986

The structure of small metal clusters

One metal atom surrounded by its 12 nearest neighbors is considered for both D(3d) (face-centered cubic-like) and D(3h) (hexagonal close-packed-like) geometries. For Al and Be, the neutral cluster and the positive and negative ions are considered for idealized (all bonds equal) and distorted geometries. The D(3d) geometry is found to be the lowest for Be13, while the D(3h) geometry is lower for Al13. This is the reverse of what is expected based upon the bulk metal structures, Be(hcp) and Al(fcc). Al13 is found to have only small distortions, while Be13 shows large distortions for both the D(3d) and D(3h) geometries. The ions have geometries which are similar to those found for the neutral systems. Both all-electron and effective core potential calculations were carried out on the X13 clusters; the agreement is very good.

Bauschlicher, C. W., Jr.↗

NASA NTRS1986

Theoretical dissociation energies for the alkali and alkaline-earth monofluorides and monochlorides

Spectroscopic parameters are accurately determined for the alkali and alkaline-earth monofluorides and monochlorides by means of ab initio self-consistent field and correlated wave function calculations. Numerical Hartree-Fock calculations are performed on selected systems to ensure that the extended Slater basis sets employed are near the Hartree-Fock limit. Since the bonding is predominantly electrostatic in origin, a strong correlation exists between the dissociation energy (to ions) and the spectroscopic parameter r(e). By dissociating to the ionic limits, most of the differential correlation effects can be embedded in the accurate experimental electron affinities and ionization potentials.

Langhoff, S. R.↗

NASA NTRS1986

On the bonding in Fe2(CO)9

The present paper is concerned with an analysis of the bonding in Fe2(CO)9, taking into account the results of self-consistent-field (SCF) calculations. The basis sets used are the same as those employed in the study of Fe(CO)5 conducted by Bauschlicher and Bagus (1985). The Fe basis set starts with the 14s9p5d primitive set optimized by Wachters (1970), which is contracted to 8s4p5d. An SCF wave function, consisting of only the core electrons, was determined, and then electrons were slowly added to the virtual orbitals with the lowest eigenvalues. The conducted calculations show no direct Fe-Fe bonding, but rather bonding resulting from the bridging CO's. The bonding is found to involve both the 5 sigma and 2 pi(asterisk) orbitals of the bridging CO's.

Bauschlicher, C. W., Jr.↗

NASA NTRS1986

Ab initio study of the alkali and alkaline-earth monohydroxides

A systematic study of the structures and dissociation energies of all the alkali and alkaline-earth monohydroxides is conducted. A theoretical model for determining accurate dissociation energies of ionic molecules is discussed. The obtained theoretical structures and dissociation energies of the alkali and alkaline-earth monohydroxides, respectively, are compared with experimental data. It is found that the theoretical studies of the bending potentials of BeOH, MgOH, and CaOH reveal the different admixture of covalent character in these systems. The BeOH molecule with the largest degree of covalent character is found to be bent (theta equals 147 deg). The MgOH is also linear. The theoretical dissociation energies for the alkali and akaline-earth hydroxides are thought to be accurate to 0.1 eV.

Bauschlicher, C. W., Jr.↗

NASA NTRS1986

Algorithms versus architectures for computational chemistry

The algorithms employed are computationally intensive and, as a result, increased performance (both algorithmic and architectural) is required to improve accuracy and to treat larger molecular systems. Several benchmark quantum chemistry codes are examined on a variety of architectures. While these codes are only a small portion of a typical quantum chemistry library, they illustrate many of the computationally intensive kernels and data manipulation requirements of some applications. Furthermore, understanding the performance of the existing algorithm on present and proposed supercomputers serves as a guide for future programs and algorithm development. The algorithms investigated are: (1) a sparse symmetric matrix vector product; (2) a four index integral transformation; and (3) the calculation of diatomic two electron Slater integrals. The vectorization strategies are examined for these algorithms for both the Cyber 205 and Cray XMP. In addition, multiprocessor implementations of the algorithms are looked at on the Cray XMP and on the MIT static data flow machine proposed by DENNIS.

Partridge, H.↗

NASA NTRS1986

Theoretical dissociation energies for ionic molecules

Ab initio calculations at the self-consistent-field and singles plus doubles configuration-interaction level are used to determine accurate spectroscopic parameters for most of the alkali and alkaline-earth fluorides, chlorides, oxides, sulfides, hydroxides, and isocyanides. Numerical Hartree-Fock (NHF) calculations are performed on selected systems to ensure that the extended Slater basis sets employed for the diatomic systems are near the Hartree-Fock limit. Extended Gaussian basis sets of at least triple-zeta plus double polarization equality are employed for the triatomic system. With this model, correlation effects are relatively small, but invariably increase the theoretical dissociation energies. The importance of correlating the electrons on both the anion and the metal is discussed. The theoretical dissociation energies are critically compared with the literature to rule out disparate experimental values. Theoretical (sup 2)Pi - (sup 2)Sigma (sup +) energy separations are presented for the alkali oxides and sulfides.

Langhoff, S. R.↗

NASA NTRS1986

Coverage dependent effects on metal surfaces - O, S, F, and Cl on Ni

The c(2 x 2) and p(2 x 2) vibrational frequencies for a rigid Ni lattice are computed for O, S, F, and Cl on Ni(100). It is found that only O/Ni(100) shows a large difference in rigid lattice vibrational frequency. A large shift in rigid lattice vibrational frequency is also found for O/Ni(111) between 0.25 and 0.5 monolayer coverages. Also associated with the change in coverage, is a change in the height above the surface R(perpendicular). This varies from a 0.11 A reduction for O/Ni(100), to an expansion of 0.06 A for F/Ni(100).

Bauschlicher, C. W., Jr.↗

NASA NTRS1986

Transition metal-ligand bonding. II

The nature of the bonding of CO, H2O, and NH3 to transition metal atoms is analyzed using the constrained-space-orbital-variation (CSOV) technique. The cooperative effects for Ni(CO)2 are found to be different than those for Ni(H2O)2. The bonding between neutral systems and the positive ions is found to be quite different; NiCO(+) has little pi bonding, while NiCO has strong pi bonding. The positive ion of NiH2O is far more strongly bound than the neutral, while for NiCO the positive ion and neutral are bound by about the same energy.

Bauschlicher, C. W., Jr.↗

NASA NTRS1986

On the uniqueness of the constrained space orbital variation (CSOV) technique

Several CSOV analyses are performed for the 1Sigma(+) state of NiCO, and it is shown that the importance of the CO sigma donation, Ni pi back donation, and interunit polarizations are virtually independent of the order of the CSOV steps, provided that the open-shell 3d sigma and 4s Ni orbitals are orthogonalized to the CO. This order of orthogonalization is consistent with the polarization of the Ni observed in the unconstrained SCF wavefunction. If instead the CO is orthogonalized to the open-shell Ni orbitals, the frozen orbital repulsion and entire CSOV analysis becomes unphysical. A comparison of the SCF and CAS SCF descriptions for the NiCO 1Sigma(+) state shows the importance of the s to d promotion and sd hybridization in reducing the repulsion and increasing the Ni to CO pi bonding. For LiF, CSOV analyses starting from both the neutral and ionic asymptotes show the bonding to be predominantly Li(+) - F(-). These examples show the uniqueness of the CSOV decomposition.

Bauschlicher, C. W., Jr.↗

NASA NTRS1985

The effect of bond functions on dissociation energies

The procedure employing bond functions recently suggested by Wright and Buenker has been applied to the N2 X 1 Sigma g + potential curve within the CAS SCF + MRSD CI treatment of electron correlation. The basis set used herein is identical to that employed by these authors in their SCF + CI calculations. The De and and the shape of the resulting potential curve, as judged by the computed vibrational levels, is not so accurate as would be expected from the results reported by Wright and Buenker (1984). The results indicate that using the CI superposition errors associated with bond functions to cancel basis set incompleteness depends on the treatment of the electron correlation.

Bauschlicher, C. W., Jr.↗

NASA NTRS1985

On 3d bonding in the transition metal trimers - The electronic structure of equilateral triangle Ca3, Sc3, Sc3(+), and Ti3(+)

It is pointed out that transition metals and transition metal (TM) compounds are currently of considerable interest because of their relevance to catalysis and to materials science problems such as hydrogen embrittlement and crack propagation in metals. The present paper is concerned with complete active space Self-Consistent Field (SCF) externally contracted configuration interaction (CASSCF/CCI) calculations for the low-lying states of Sc3 and Sc3(+). A comparison is conducted regarding the bonding in the Ca3, Sc3, and Cu3 molecules. This comparison makes it possible to predict general trends for the TM trimers. Attention is given to the qualitative features of the bonding in the TM trimers, the basis sets and other technical details of the calculations, the calculated results for Sc3 and Sc3(+), and conclusions from this work.

Walch, S. P.↗

NASA NTRS1985

Calculated electric dipole moment of NiH X2Delta

A calculated dipole moment of 2.39 D at R sub e = 2.79 a sub 0 is reported, obtained from complete active space SCF/configuration interaction calculations plus one natural orbital iteration. The calculation is in good agreement with the experimental value of 2.4 + or - 0.1 D measured for the lowest vibrational level. In agreement with Gray et al. (1985), it is found that the dipole moment is strongly correlated with the 3d electron population; the good agreement with experiment thus provides verification of the mixed state model of NiH. It is concluded that the electric dipole moment of NiH is a sensitive test of the quality of the NiH wave function.

Walch, S.↗

NASA NTRS1985

H2O/Ni(100) and NH3/Ni(100) - A computational approach

The adsorption of NH3 and H2O on the Ni(100) is treated using a cluster model. The adsorption is found to have only a small effect on the HXH angle. Rotation about the principal ligand axis requires virtually no energy. Ligand tilts of 15 deg are found to require small amounts of energy, leading to the suggestion that the ESDIAD signal is a result of excited tilting modes, not a change in the HXH angle. For H2O additional bending modes are considered and all are found to be quite flat in energy.

Bauschlicher, C. W., Jr.↗

NASA NTRS1985

Bonding of NH3 to Cu

Bagus et al. (1984) have studied the chemisorption of NH3 at the on-top site of Cu(100), using a cluster model, Cu5(1,4)NH3. On the basis of the value of the NH3 chemisorption energy obtained, it was concluded that NH3 chemisorbs directly above a Cu atom. In order to consider also other sites, Al(111) was substituted for Cu(100). Chemisorption on Al can be quite different from a transition metal, however, and questions arise regarding the suitability of the choice of the Al model. Thus, in this note, NH3 into hollows is considered for Cu(111). Attention is given to the conduction of self-consistent field calculations, problems of site determination, and the consideration of electrostatic interaction.

Bauschlicher, C. W., Jr.↗

NASA NTRS1985

The coverage-dependent effect on the barrier for oxygen penetration into a metal lattice

Model calculations for O/Cu(111) show that the rigid lattice penetration barrier is reduced from 8.2 to 5.3 eV as the O coverage is increased from 1/4 to 1/2 monolayer; it is expected that this large reduction will also apply to the relaxed lattice case. An increase of the distance between the first and second Cu layers, without changing the intralayer spacing, does not affect the barrier.

Bauschlicher, C. W., Jr.↗

NASA NTRS1985

The lower electronic states of ClOO - A computational investigation

Eight doublet and eight quartet states of ClOO were investigated by ab initio CI techniques. The potential energy surfaces of the four lowest energy doublet states of both A-double-prime and A-prime symmetry indicate that only the 1 2A-double-prime state is bound. In contrast to the model provided by the HO2 radical, all of the excited doublet states investigated were repulsive with respect to dissociation to Cl + O2 and metastable or bound with respect to dissociation to ClO + O. The transitions to the excited states investigated span the visible and near UV spectral regions, but the transition moments indicate that they are very weak. Since the photolysis products are the same as those of the rapid thermal dissociation, photolysis is not expected to be an important atmospheric process. The soft bending potential for the 1 2A-prime state and the shape of the 1 4A-double-prime state in the entrance channel of the ClO + O yields Cl + O2 reaction provide a qualitative explanation for the underprediction of the low temperature reaction rate by previous trajectory calculations.

Jafri, J. A.↗