Engineering topics

Langhoff, S. R.

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

At least 19 records

NASA NTRS1999

Carbon chain abundance in the diffuse interstellar medium

Thanks to the mid-IR sensitivities of the ISO and IRTS orbiting spectrometers it is now possible to search the diffuse interstellar medium for heretofore inaccessible molecular emission. In view of the recent strong case for the presence of C(7-) (Kirkwood et al. 1998, Tulej et al. 1998),and the fact that carbon chains possess prominent infrared active modes in a very clean portion of the interstellar spectrum, we have analyzed the IRTS spectrum of the diffuse interstellar medium for the infrared signatures of these species. Theoretical and experimental infrared band frequencies and absolute intensities of many different carbon chain species are presented. These include cyanopolyynes, neutral and anionic linear carbon molecules, and neutral and ionized, even-numbered, hydrogenated carbon chains. We show that--as a family--these species have abundances in the diffuse ISM on the order of 10(-10) with respect to hydrogen, values consistent with their abundances in dense molecular clouds. Assuming an average length of 10 C atoms per C-chain implies that roughly a millionth of the cosmically available carbon is in the form of carbon chains and that carbon chains can account for a few percent of the visible to near-IR diffuse interstellar band (DIB) total equivalent width (not DIB number).

Carbon/analysis/chemistry↗

NASA NTRS1998

Infrared spectra of substituted polycylic aromatic hydrocarbons

Calculations are carried out using density functional theory (DFT) to determine the harmonic frequencies and intensities of 1-methylanthracene, 9-methylanthracene, 9-cyanoanthracene, 2-aminoanthracene, acridine, and their positive ions. The theoretical data are compared with matrix-isolation spectra for these species also reported in this work. The theoretical and experimental frequencies and relative intensities for the neutral species are in generally good agreement, whereas the positive ion spectra are only in qualitative agreement. Relative to anthracene, we find that substitution of a methyl or CN for a hydrogen does not significantly affect the spectrum other than to add the characteristic methyl C-H and C triple bond N stretches near 2900 and 2200 cm-1, respectively. However, addition of NH2 dramatically affects the spectrum of the neutral. Not only are the NH2 modes themselves strong, but this electron-withdrawing group induces sufficient partial charge on the ring to give the neutral molecule spectra characteristics of the anthracene cation. The sum of the absolute intensities is about four times larger for 2-aminoanthracene than those for 9-cyanoanthracene. Substituting nitrogen in the ring at the nine position (acridine) does not greatly alter the spectrum compared with anthracene.

NASA Discipline Exobiology↗

NASA NTRS1995

Minimum Energy Pathways for Chemical Reactions

Computed potential energy surfaces are often required for computation of such parameters as rate constants as a function of temperature, product branching ratios, and other detailed properties. We have found that computation of the stationary points/reaction pathways using CASSCF/derivative methods, followed by use of the internally contracted CI method to obtain accurate energetics, gives useful results for a number of chemically important systems. The talk will focus on a number of applications to reactions leading to NOx and soot formation in hydrocarbon combustion.

Walch, S. P.↗

NASA NTRS1994

Intermediates in the Formation of Aromatics in Hydrocarbon Combustion

The formation of the first benzene ring is believed to be the rate limiting step in soot formation. Two different mechanisms have been proposed for formation of cyclic C6 species. The first involves the reaction of two acetylenes to give CH2CHCCH (vinyl acetylene), the loss of a H to give CHCHCCH (n-C41-13) or CH2CCCH (iso-C4H3), and addition of another acetylene to n-C4H3, followed by ring closure to give phenyl radical. Miller and Melius argue that only n-C4H3 leads to phenyl radical and since iso-C4H3 is more stable than n-C4H3 this mechanism is unlikely. An alternative mechanism proposed by them is formation of benzene from the dimerization of two CH2CCH (propargyl) radicals (formed by the reaction of singlet methylene with C2H2). We report reaction pathways and accurate energetics (from CASSCF/internally contracted CI calculations) for the reactions of CH(pi-2) and CH2-1 with acetylene, the reaction of vinylidene with acetylene, and the reaction of n-C4H3 and iso-C4H3 with acetylene. These calculations identify two new reactive intermediates CHCHCH ( a A"-2 ground state in Cs symmetry; spin coupling is a doublet from three singly occupied orbitals) and CHCCH (B-3 ground state in C2 symmetry) from the reaction of CH with acetylene. These species dimerize with no barrier to form benzene and para-benzyne, respectively. CHCCH is proposed as a reactive intermediate which can add to benzene to give higher polynuclear aromatic hydrocarbons or fullerenes. The addition of a C3H2 unit releases two C-C bond energies and thus the resulting addition product contains sufficient energy to break several CH bonds leading to a reduction in the H to C ratio as the cluster size increases. It is found that iso-C4H3 adds to acetylene to initially give a fulvene radical but that this species rearranges to phenyl radical. Thus, the reaction of acetylene with iso-C4H3 does lead to phenyl radical and the cyclization pathway may also contribute to formation of the initial benzene ring.

Walch, Stephen P.↗

NASA NTRS1994

Computed Potential Energy Surfaces and Minimum Energy Pathway for Chemical Reactions

Computed potential energy surfaces are often required for computation of such observables as rate constants as a function of temperature, product branching ratios, and other detailed properties. We have found that computation of the stationary points/reaction pathways using CASSCF/derivative methods, followed by use of the internally contracted CI method with the Dunning correlation consistent basis sets to obtain accurate energetics, gives useful results for a number of chemically important systems. Applications to complex reactions leading to NO and soot formation in hydrocarbon combustion are discussed.

Walch, Stephen P.↗

NASA NTRS1994

Characterization of the Minimum Energy Paths for the Reactions of CH(X(sup 2 Pi) and (1)CH2 with C2H2

The reactions of CH(sup 2 Pi) and singlet methylene (1)CH2 with acetylene lead to intermediates which may be important in soot formation. CH(sup 2 Pi) + acetylene leads to CHCHCH (C3H3), CHCCH (C3H2), and propargyl (CH2CCH). (1)CH2 + acetylene leads to cyclopropene and propargyl. All of these reaction products are formed with no barrier. Miller and Melius have previously discussed the dimerization of propargyl to give benzene. C3H3 and C3H2 can dimerize with no barrier to give benzene and para-benzyne, respectively. C3H3 and C3H2 can also add to smaller polynuclear aromatic hydrocarbons (PAH), and may be important species in forming larger PAH or fullerenes.

Walch, Stephen P.↗

NASA NTRS1994

Accurate Energetics for the Unimolecular Decomposition of HN2

We report internally contracted configuration interaction (ICCI) calculations along the constrained energy minimum (CEM) path for the addition of H to N2 with both CC-pVQZ and augmented CC-pVQZ basis sets. The results agree to within 0.5 kcal/mol with the larger basis set results of Walch, which were used to calibrate the HN2 potential energy surface used by Koizumi, Schatz, and Walch in coupled channel calculations to determine the lifetimes for unimolecular decay of the lowest vibrational levels of HN2. These small changes in the potential energy surface should have no significant effect on the computed lifetimes.

Walch, Stephen P.↗

NASA NTRS1994

Computed Potential Energy Surfaces and Minimum Energy Pathways for Chemical Reactions

Computed potential energy surfaces are often required for computation of such parameters as rate constants as a function of temperature, product branching ratios, and other detailed properties. For some dynamics methods, global potential energy surfaces are required. In this case, it is necessary to obtain the energy at a complete sampling of all the possible arrangements of the nuclei, which are energetically accessible, and then a fitting function must be obtained to interpolate between the computed points. In other cases, characterization of the stationary points and the reaction pathway connecting them is sufficient. These properties may be readily obtained using analytical derivative methods. We have found that computation of the stationary points/reaction pathways using CASSCF/derivative methods, followed by use of the internally contracted CI method to obtain accurate energetics, gives usefull results for a number of chemically important systems. The talk will focus on a number of applications including global potential energy surfaces, H + O2, H + N2, O(3p) + H2, and reaction pathways for complex reactions, including reactions leading to NO and soot formation in hydrocarbon combustion.

Walch, Stephen P.↗

NASA NTRS1994

Computed Minimum Energy Pathway for Isomerization in Ketene

A minimum energy pathway for interchange of the CH groups in ketene via a C2v structure has been obtained using complete active space self consistent field (CASSCF) derivative methods with a polarized valence double zeta basis set to define the reaction pathway followed by multi-reference internally contracted configuration interaction (ICCI) calculations with a [3s3p2d/3s2p] basis set to determine the energetics. Qualitatively, the C2v structure is found to be a shallow minimum on the potential energy surface separated from ketene by a small barrier (0.2 kcal/mol), a second minimum, and a larger barrier (3.0 kcal/mol). The minimum energy pathway leading from the C2v minimum to ketene starts by simultaneous rotation of the farther CH group out of the plane and away from the oxygen followed by increase of the CCO angle and subsequent 1,2-migration of the H of the nearer CH group toward the carbon of the farther CH group.

Walch, Stephen P.↗

NASA NTRS1986

On the electron affinity of the oxygen atom

The electron affinity of oxygen is computed to be 1.287 eV, at the full CI level using a 6s5p3d 2f Slater-type orbital basis and correlating only the 2p electrons. The best CASSCF-MRCI result including only 2p correlation is 1.263 eV. However, inclusion of 2s intrashell and 2s2p intershell coorelation increases the computed EA to 1.290 eV at the CASSCF-MRCI level. At the full CI basis set limit, the 2s contribution to the electron affinity is estimated to be as large as 0.1 eV. The higher excitation contribution to the electron affinity is found to increase substantially with basis set completeness, especially when the 2s electrons are correlated. Relativistic effects are shown to make a small (less than 0.01 eV) change in the EA.

Bauschlicher, C. W., Jr.↗

NASA NTRS1986

Theoretical study of the electric dipole moment function of the ClO molecule

The potential energy function and electric dipole moment function (EDMF) are computed for ClO X 2Pi 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. The band strength of the 1-0 fundamental transition is computed to be 12 + or - 2 per 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.

Pettersson, L. G. M.↗

NASA NTRS1986

Theoretical dipole moments for the first-row transition metal hydrides

Spectroscopic parameters (D sub e, r sub e, mu) are determined for the first-row transition metal hydrides using better than DZP basis sets at the modified coupled pair functional (MCPF) level. Extensive comparisons between MCPF and complete-active space self-consistent field (CASSCF)/MRCI calculations with natural orbital iterations, and studies with more extensive basis sets, show this level of treatment to supply an accurate and cost-effective treatment of these systems. For the transition metal hydrides, the bonding can arise from either the 3d(n)4s(2) or 3d(n + 1)4s(1) atomic asymptotes, or a mixture of both. Since the dipole moment arising from these two bonding mechanisms is very different, the dipole moment is found to be directly related to the 3d population. Thus, the magnitude of the dipole moments provides a sensitive test of the wave function, and gives insight into the nature of the bonding.

Chong, D. P.↗

NASA NTRS1986

Benchmark full configuration-interaction calculations on HF and NH2

Full configuration-interaction (FCI) calculations are performed at selected geometries for the 1-sigma(+) state of HF and the 2-B(1) and 2-A(1) states of NH2 using both DZ and DZP gaussian basis sets. Higher excitations become more important when the bonds are stretched and the self-consistent field (SCF) reference becomes a poorer zeroth-order description of the wave function. The complete active space SCF - multireference configuration-interaction (CASSCF-MRCI) procedure gives excellent agreement with the FCI potentials, especially when corrected with a multi-reference analog of the Davidson correction.

Bauschlicher, C. W., Jr.↗

NASA NTRS1986

Benchmark full configuration-interaction calculations on HF and NH2

Bauschlicher, C. W., Jr.↗

NASA NTRS1986

On the electron affinity of the oxygen atom

The electron affinity (EA) of oxygen is computed to be 1.287 eV, using 2p electron full configuration-interaction (CI) wave functions expanded in a 6s5p3d2f Slater-type orbital basis. The best complete active space self-consistent field - multireference CI (CASSCF-MRCI) result including only 2p correlation is 1.263 eV. However, inclusion of 2s intrashell and 2s2p intershell correlation increases the computed EA to 1.290 at the CASSCF-MRCI level. At the full CI basis set limit, the 2s contribution to the electron affinity is estimated to be as large as 0.1 eV. This study clearly establishes the synergistic effect between the higher excitations and basis set completeness on the electron affinity when the 2s electrons are correlated.

Bauschlicher, C. W., Jr.↗

NASA NTRS1986

Theoretical study of the 7Sigma(u)+ state of N2

Theoretical potentials for the 7Sigma(u)+ state are reported using both extended Slater and Gaussian basis sets. Electron correlation is included using the interacting correlated fragments (ICF), the singles plus doubles configuration-interaction (SDCI), and the coupled-pair functional (CPF) approaches. The best potential, corrected for basis set superposition errors, has a well depth of about 21/cm and an r(e) of 7.52 bohr. The inclusion of the nitrogen 2s correlation significantly increases the well depth. These results are in reasonably good agreement with the empirical potential of Ferrante and Stwalley (1983), and support the contention that spin-polarized atomic nitrogen should behave like a classical solid such as Ne. It is shown that the modified Buckingham potential used by Ferrante and Stwalley is consistent with the present ab initio potentials if smaller values are used for the dispersion coefficients of the nitrogen atom.

Partridge, H.↗

NASA NTRS1986

A full CI treatment of Ne atom - A benchmark calculation performed on the NAS CRAY 2

Full CI calculations are performed for Ne atom using Gaussian basis sets of up to triple-zeta plus double polarization quality. The total valence correlation energy through double, triple, quadruple and octuple excitations is compared for eight different basis sets. These results are expected to be an important benchmark for calibrating methods for estimating the importance of higher excitations.

Bauschlicher, C. W., Jr.↗

NASA NTRS1986

A modified coupled pair functional approach

A modified coupled pair functional (CPF) method is presented for the configuration interaction problem that dramatically improves properties for cases where the Hartree-Fock reference configuration is not a good zeroth-order wave function description. It is shown that the tendency for CPF to overestimate the effect of higher excitations arises from the choice of the geometric mean for the partial normalization denominator. The modified method is demonstrated for ground state dipole moment calculations of the NiH, CuH, and ZnH transition metal hydrides, and compared to singles-plus-doubles configuration interaction and the Ahlrichs et al. (1984) CPF method.

Chong, D. P.↗