Engineering Papers⌕ Search

Engineering topics

Langhoff, Stephen R.

Publications and source records attributed to Langhoff, Stephen R..

At least 73 records · Page 4

Ab Initio Potential Energy Surfaces and the Calculation of Accurate Vibrational Frequencies

Due to advances in quantum mechanical methods over the last few years, it is now possible to determine ab initio potential energy surfaces in which fundamental vibrational frequencies are accurate to within plus or minus 8 cm(exp -1) on average, and molecular bond distances are accurate to within plus or minus 0.001-0.003 Angstroms, depending on the nature of the bond. That is, the potential energy surfaces have not been scaled or empirically adjusted in any way, showing that theoretical methods have progressed to the point of being useful in analyzing spectra that are not from a tightly controlled laboratory environment, such as vibrational spectra from the interstellar medium. Some recent examples demonstrating this accuracy will be presented and discussed. These include the HNO, CH4, C2H4, and ClCN molecules. The HNO molecule is interesting due to the very large H-N anharmonicity, while ClCN has a very large Fermi resonance. The ab initio studies for the CH4 and C2H4 molecules present the first accurate full quartic force fields of any kind (i.e., whether theoretical or empirical) for a five-atom and six-atom system, respectively.

Lee, Timothy J.↗

Ab Initio Studies of Chlorine Oxide and Nitrogen Oxide Species of Interest in Stratospheric Chemistry

The ability of modern state-of-the art ab initio quantum chemical techniques to characterize reliably the gas-phase molecular structure, vibrational spectrum, electronic spectrum, and thermal stability of chlorine oxide and nitrogen oxide species will be demonstrated by presentation of some example studies. In particular the geometrical structures, vibrational spectra, and heats of formation Of ClNO2, CisClONO, and trans-ClONO are shown to be in excellent agreement with the available experimental data, and where the experimental data are either not known or are inconclusive, the ab initio results are shown to fill in the gaps and to resolve the experimental controversy. In addition, ab initio studies in which the electronic spectra and the characterization of excited electronic states of ClONO2, HONO2, ClOOC17 ClOOH, and HOOH will also be presented. Again where available, the ab initio results are compared to experimental observations, and are used to aid in the interpretation of the experimental studies.

Lee, Timothy J.↗

Ab Initio Studies of Halogen and Nitrogen Oxide Species of Interest in Stratospheric Chemistry

The ability of modern state-of-the art ab initio quantum chemical techniques to characterize reliably the gas-phase molecular structure, vibrational spectrum, electronic spectrum, and thermal stability of fluorine, chlorine, bromine and nitrogen oxide species will be demonstrated by presentation of some example studies. The ab initio results are shown to be in excellent agreement with the available experimental data, and where the experimental data are either not known or are inconclusive, the theoretical results are shown to fill in the gaps and to resolve experimental controversies. In addition, ab initio studies in which the electronic spectra and the characterization of excited electronic states of halogen oxide species will also be presented. Again where available, the ab initio results are compared to experimental observations, and are used to aid in the interpretation of experimental studies.

Lee, Timothy J.↗

A Comparison of Correlation-Consistent and Pople-Type Basis Sets

Correlation-consistent basis sets are compared with Pople-type basis sets for the calculation of atomization energies using the G2(B3LYP/MP2/CC) approach. Without a higher-level correction, the smallest errors are found for the augmented- correlation-consistent polarized valence triple zeta plus tight d basis sets. However, after inclusion of a higher-level correction, the combination of the 6-311G** basis set for the high-level correlation treatment and 6-311+G(3df,2p) for the MP2 step yields the best results. At the spin-unrestricted CCSD(T) level and for the basis set considered in this work, the G2(B3LYP/MP2/CC) approach works better than the PCIX approach for the atomization energies of the 55 G2 molecules.

Bauschlicher, Charles W., Jr.↗

Vibrational and Electronic Energy Transfer and Dissociation of Diatomic Molecules by Electron Collisions

At high altitudes and velocities equal to or greater than the geosynchronous return velocity (10 kilometers per second), the shock layer of a hypersonic flight will be in thermochemical nonequilibrium and partially ionized. The amount of ionization is determined by the velocity. For a trans atmospheric flight of 10 kilometers per second and at an altitude of 80 kilometers, a maximum of 1% ionization is expected. At a velocity of 12 - 17 kilometer per second, such as a Mars return mission, up to 30% of the atoms and molecules in the flow field will be ionized. Under those circumstances, electrons play an important role in determining the internal states of atoms and molecules in the flow field and hence the amount of radiative heat load and the distance it takes for the flow field to re-establish equilibrium. Electron collisions provide an effective means of transferring energy even when the electron number density is as low as 1%. Because the mass of an electron is 12,760 times smaller than the reduced mass of N2, its average speed, and hence its average collision frequency, is more than 100 times larger. Even in the slightly ionized regime with only 1% electrons, the frequency of electron-molecule collisions is equal to or larger than that of molecule-molecule collisions, an important consideration in the low density part of the atmosphere. Three electron-molecule collision processes relevant to hypersonic flows will be considered: (1) vibrational excitation/de-excitation of a diatomic molecule by electron impact, (2) electronic excitation/de-excitation, and (3) dissociative recombination in electron-diatomic ion collisions. A review of available data, both theory and experiment, will be given. Particular attention will be paid to tailoring the molecular physics to the condition of hypersonic flows. For example, the high rotational temperatures in a hypersonic flow field means that most experimental data carried out under room temperatures are not applicable. Also, the average electron temperature is expected to be between 10,000 and 20,000 K. Thus only data for low energy electrons are relevant to the model.

Huo, Winifred M.↗

Rotational Energy Transfer and Collisional Induced Raman Linewidths in N2 Gas: Energy Transfer Rates - 1

Rotationally inelastic transitions of N2 have been studied in the coupled state (CS) and infinite-order-sudden (IOS) approximations, using the N2-N2 rigidrotor potential of van der Avoird et al. For benchmarking purposes, close coupling (CC) calculations have also been carried out over a limited energy range. The CC and CS cross sections have been obtained both with and without identical molecule exchange symmetry, whereas exchange was neglected in the IOS calculations. The CS results track the CC cross sections rather well; between 113 - 219 cm(exp -1) the average deviation is 14%. Comparison between the CS and IOS cross sections at the high energy end of the CS calculation, 500 - 680 cm(exp -1), shows that IOS is sensitive to the amount of inelasticity and the results for large DELTA J transitions are subject to larger errors. It is found that the state-to-state cross sections with even and odd exchange symmetry agree to better than 2% and are well represented as a sum of direct and exchange cross sections for distinguishable molecules, an indication of the applicability of a classical treatment for this system. This result, however, does not apply to partial cross sections for given total J, but arises from a near cancellation in summing over partial waves. In order to use rigid-rotor results for the calculation of effective rotational excitation rates of N2 in the v=1 vibrational level colliding with bath N2 molecules in the v=0 level, it is assumed that exchange scattering between molecules in different vibrational levels is negligible and direct scattering is independent of Y. Good agreement with room temperature experimental data is obtained. The effective rates determined using the IOS and energy corrected sudden (ECS) approximations are also in reasonable agreement with experiment, with the ECS results being somewhat better. The problem with a degeneracy factor in earlier cross section expressions for collisions between identical molecules is pointed out and corrected.

Huo, Winifred M.↗

Rotational Energy Transfer Cross Sections in N2-N2 Collisions

Rotational inelastic transitions of N2 have been studied in the coupled state (CS, also called centrifugal sudden) and infinite-order-sudden (IOS) approximations, using the N2-N2 rigid-rotor potential of van der Avoird et al. For benchmarking purposes, close coupling (CC) calculations have also been carried out over a limited energy range and for even j - even j collisions only. Both the CC and CS cross sections have been obtained with and without exchange symmetry, whereas exchange is neglected in the IOS calculations. The CS results track the CC cross sections rather well. At total energies between 113 to 219 cm(exp -1) the average deviation is 14%. The deviation decrease with increasing energy, indicating that the CS approximation can be used as a substitute at higher energies when the CC calculations become impractical. Comparison between the CS and IOS cross sections at the high energy end of the CS calculation, 500 - 680 cm(exp-1), shows significant differences between the two. In addition, the IOS results exhibits sensitivity to the amount of inelasticity and the results for large DELTA J transitions are subjected to bigger errors. At total energy 113 cm(exp -1) and above, the average deviation between state-to-state cross sections calculated with even and odd exchange symmetries is 1.5%.

Huo, Winifred M.↗

The FEM-R-Matrix Approach: Use of Mixed Finite Element and Gaussian Basis Sets for Electron Molecule Collisions

For the calculation of electron molecule collision cross sections R-matrix methods automatically take advantage of the division of configuration space into an inner region (I) bounded by radius tau b, where the scattered electron is within the molecular charge cloud and the system is described by an correlated Configuration Interaction (CI) treatment in close analogy to bound state calculations, and an outer region (II) where the scattered electron moves in the long-range multipole potential of the target and efficient analytic methods can be used for solving the asymptotic Schroedinger equation plus boundary conditions.

Thuemmel, Helmar T.↗

The MCH2+ Systems: Do ScCH2+ and TiCH2+ Have C(sub s) or C(sub 2v) Symmetry and A Comparison of the B3LYP Method to Other Approaches

ScCH2(+) and TiCH2(+) are found to have C(sub s) symmetry, while the remaining systems have the expected C(sub 2v) symmetry. The C(sub s) symmetry structure is favored for ScCH2(+) and TiCH2(+) because it allows donation from one of the CH bonds into an empty 3d orbital. This distortion stabilizes the systems by less than 2 kcal/mol at the B3LYP, CASSCF, and CCSD(T) levels of theory. The SCF and B3LYP approaches find the C2(sub v) structure to be a saddle point, while the CASSCF approach finds the C(sub 2v) structure to be a stationary point. For V-Cu there are no empty 3d orbitals of the correct symmetry, and therefore the distortion to C(sub s), is unfavorable. The B3LYP binding energies are in good agreement with experiment and our previous best estimates.

Ricca, Alessandra↗

Rotational And Rovibrational Energy Transfer In Electron Collisions With Molecules

Air flows around a hypervelocity reentry vehicle undergo dissociation, rovibrational excitation and ionization. More specifically the air, initially 80% N2 and 20% O2, in the shock layer consists of species such as N, O, N2, O2, NO, N+, O+, N+, O+, NO+ and 2 free electrons. It was pointed out in multi temperature models'' that the temperature of the rotational energy modes and the gas-kinetic translational temperature are quickly equilibrated by a few collisions and rise rapidly to high temperatures as 50000K before falling off to equilibrium value of 10000K. Contrary, the electronic and vibrational temperatures state energy distributions remain low (less than 15000K) because of the slow equilibration. Electron vibrational energy transfer is thought to play a crucial role in such a ionizing flow regime since chemical reaction rates and dissociation depend strongly on the vibrational temperatures. Modeling of these flowfields in principle require the rovibrational excitation and de-excitation cross section data for average electron energies from threshold up to several eV (leV=11605.4 K). In this lecture we focus on theoretical description of rotational effects i.e. energy transfer of electrons to molecules such that the molecular rotational (vojo goes to voj) or vibrational and rotational (v(sub 0)j(sub 0) goes to vj) states are changed. Excitation and de-excitation of electronic states was discussed in a previous talk at this conference.

Thuemmel, Helmar T.↗

Potential Energy Curves and Collisions Integrals of Air Components: Interactions Involving Ionized Atoms - 2

Collision integrals are fundamental quantities required to determine the transport properties of the environment surrounding aerospace vehicles in the upper atmosphere. These collision integrals can be determined as a function of temperature from the potential energy curves describing the atomic and molecular collisions. Ab initio calculations provide a practical method of computing the required interaction potentials. In this work we will discuss recent advances in scattering calculations with an emphasis on the accuracy that is obtainable. Results for interactions of the atoms and ionized atoms of nitrogen and oxygen will be reviewed and their application to the determination of transport properties, such as diffusion and viscosity coefficients, will be examined.

Stallcop, James R.↗

Variational Calculations of Ro-Vibrational Energy Levels and Transition Intensities for Tetratomic Molecules

A description is given of an algorithm for computing ro-vibrational energy levels for tetratomic molecules. The expressions required for evaluating transition intensities are also given. The variational principle is used to determine the energy levels and the kinetic energy operator is simple and evaluated exactly. The computational procedure is split up into the determination of one dimensional radial basis functions, the computation of a contracted rotational-bending basis, followed by a final variational step coupling all degrees of freedom. An angular basis is proposed whereby the rotational-bending contraction takes place in three steps. Angular matrix elements of the potential are evaluated by expansion in terms of a suitable basis and the angular integrals are given in a factorized form which simplifies their evaluation. The basis functions in the final variational step have the full permutation symmetries of the identical particles. Sample results are given for HCCH and BH3.

Schwenke, David W.↗

Theoretical Studies of the Quasi-Bound Vibrational Levels in the B State of CH3

Consider the calculation of the oscillator strength of th B yielded by X electronic transition in CH3. The X state is strongly bound, but in the B state, CH3 is pre-dissociated by H atom tunneling through a potential energy barrier to form CH2+H. If we turn off pre-dissociation by making the potential outside of the barrier repulsive, the calculated 0-0 oscillator strength for CD3 from high quality ab initio potential and transition moment data along with accurate ro-vibrational wavefunctions is 0.0212, which is a factor of two greater than experiment. This discrepancy is largely due to the neglect of pre-dissociation. In this poster we explore techniques for including pre-dissociation in the calculations.

Schwenke, David W.↗

A Comparison of the Accuracy of Different Functionals

The errors in the computed geometries, zero-point energies, and atomization energies of molecules containing only first and second row atoms are compared for several levels of theory, including Hartree-Fock, second-order Moller-Plesset perturbation theory (MP2), and density functional theory (DFT) using five different functionals, including two hybrid functionals. The 6-31G* and 6-311+G(3df,2p) basis set are used. Overall, the B3LYP hybrid functional yields the best results.

Bauschlicher, Charles W., Jr.↗

The Structure of Li7 and K7

The self-consistent-field (SCF) approach and density functional theory, using the B3LYP hybrid functional, yield three low-lying structures for Li7(-). The relative separations differ for the SCF and B3LYP approaches, however the B3LYP results are in good agreement with the coupled cluster results. For K7(-), only an octahedron with one face capped is found to be a minimum; this the second most stable structure for Li7(-). A comparison of the computed separations between the low-lying states of K7 and the photoelectron detachment spectra does not allow an unambiguous assignment of the structure of K7(-).

Bauschlicher. Charles W., Jr.↗

The low-lying electronic states of LiC

The spectroscopic constants for the doublet and quartet states of LiC below about 30,000/cm are determined using an internally contracted multireference configuration-interaction approach in conjunction with a [6s 5p 3d 2f] atomic natural orbital basis sets. All of the strongly bound states, X(sup 4)(SIGMA)(sup -),(1)(sup 2)(DELTA), (1)(sup 2)(SIGMA)(sup +), and (2)(sup 2) II, very ionic in character. The only bound-bound quartet transition in this energy range is (2)(sup 4)SIGMA(sup -) and Franck-Condon factors, Einstein A values, and lifetimes are reported for this transition.

Ricca, Alessandra↗

Theoretical Study of Low-Energy Electron-Molecule Collisions

Polarized studies of electron collisions with carbon tetrafluoride and beryllium carbonyl are reported. In the e-CF4 study, the resonance structure between 8-9 eV is shown to be a superposition of the T2- and A1-symmetry shape resonances. The unequivocal assignment of this feature as a double-resonance structure clarifies certain discrepancies in previous attempts to ascribe the 8-9 eV feature to a shape resonance in one or another partial channel and explains the sensitivities found in the fragment ion production in the resonance region. The Be-CO system is chosen as a prototypical example of an absorbate-substrate interaction. To understand how bonding with the substrate affects the CO resonance, electron collision with two states of BeCO are studied, one with a very weak van der Waals bond and the other, with a normal chemical bond. It is shown that the nature of the Be-CO bond has strong effects on the cross section feature.

Huo, Winfred M.↗

The Dissociation Energies of CH4 and C2H2 Revisited

The bond dissociation energies of CH4 and C2H2 and their fragments are investigated using basis set extrapolations and high levels of correlation. The computed bond dissociation energies (D(sub e)) are accurate to within 0.2 kcal/mol. The agreement with the experimental (D(sub 0)) values is excellent if we assume that the zero-point energy of C2H is 9.18 kcal/mol. The effect of core (1s) correlation on the bond dissociation energies of C-H bonds is shown to vary from 0.2 to 0.7 kcal/mol and that for C-C bonds varies from 0.4 to 2.2 kcal/mol.

Partridge, Harry↗