Engineering topics

Walch, Stephen P.

Publications and source records attributed to Walch, Stephen P..

At least 19 records

NASA NTRS2004

The Bonding of NO2, NH3, and CH2NY to Models of a (10,0) Carbon Nanotube

We have studied the bonding of NO2, NH3, and CH2NH to a (10,O) carbon nanotube using the MP2 and ONIOM methods with extended basis sets. We find bond strengths of 3.5, 3.6. and 6.3 kcal/mol for NO2, NH3, and CH2NH, respectively, using the ONIOM method with the high accuracy part treated at the MP2/aug-CC-pVTZ level and the remainder of the CNT approximated at the UFF level and including an estimate of basis set superposition error using the counterpoise method.

Walch, Stephen P.↗

NASA NTRS2002

Gas-Phase Reaction Pathways and Rate Coefficients for the Dichlorosilane-Hydrogen and Trichlorosilane-Hydrogen Systems

As part of NASA Ames Research Center's Integrated Process Team on Device/Process Modeling and Nanotechnology our goal is to create/contribute to a gas-phase chemical database for use in modeling microelectronics devices. In particular, we use ab initio methods to determine chemical reaction pathways and to evaluate reaction rate coefficients. Our initial studies concern reactions involved in the dichlorosilane-hydrogen (SiCl2H2--H2) and trichlorosilane-hydrogen (SiCl2H-H2) systems. Reactant, saddle point (transition state), and product geometries and their vibrational harmonic frequencies are determined using the complete-active-space self-consistent-field (CASSCF) electronic structure method with the correlation consistent polarized valence double-zeta basis set (cc-pVDZ). Reaction pathways are constructed by following the imaginary frequency mode of the saddle point to both the reactant and product. Accurate energetics are determined using the singles and doubles coupled-cluster method that includes a perturbational estimate of the effects of connected triple excitations (CCSD(T)) extrapolated to the complete basis set limit. Using the data from the electronic structure calculations, reaction rate coefficients are obtained using conventional and variational transition state and RRKM theories.

Dateo, Christopher E.↗

NASA NTRS2001

Reactions of SiCl2 and SiHCl with H and Cl Atoms

Calculations have been carried out for the reaction of SiCl2 and SiHCl with H and Cl atoms. In each case, the stationary point geometries and harmonic frequencies were characterized using CASSCF/derivative methods and the cc-pVDZ basis set. Accurate energetics were obtained by combining the CCSD(T) results using the a-cc-pVTZ basis set with an extrapolation to the basis set limit using the a-cc-pVDZ, a-cc-pVTZ, and a-cc-pVQZ basis sets at the MP2 level. The geometries, energetics, and harmonic frequencies were used to obtain rate constants using conventional transition state theory or a Gorin-like model. In each case we find direct abstraction pathways compete with an addition elimination pathway. In the case of SiClH + H the two direct pathways are H abstraction which is barrierless and Cl abstraction with a barrier of 13.5 kcal/mol, while the addition elimination process has a barrier of 26.9 kcal/mol. In the case of SiCl2 + H the direct pathway is Cl abstraction with a barrier of 16.4 kcal/mol, while the addition elimination pathway has a barrier of 29.6 kcal/mol. In the case of SiClH + Cl the direct pathway is H abstraction which is barrierless and the addition elimination pathway has a barrier of 2.0 kcal/mol.

Walch, Stephen P.↗

NASA NTRS2000

Atomistic Design and Simulations of Nanoscale Machines and Assembly

Over the three years of this project, we made significant progress on critical theoretical and computational issues in nanoscale science and technology, particularly in:(1) Fullerenes and nanotubes, (2) Characterization of surfaces of diamond and silicon for NEMS applications, (3) Nanoscale machine and assemblies, (4) Organic nanostructures and dendrimers, (5) Nanoscale confinement and nanotribology, (6) Dynamic response of nanoscale structures nanowires (metals, tubes, fullerenes), (7) Thermal transport in nanostructures.

Goddard, William A., III↗

NASA NTRS2000

On the reaction CH2O + NH3 Yields CH2NH + H2O

We study the energetics of CH2O + NH3-Yields CH2NH + H2O in the gas phase at the B3LYP and CCSD(T)levels. This reaction is shown to have a sizeable barrier. Ionization of NH3 reduces the barrier to about 5 kcal/mol. We also show that in water, a proton catalyzed mechanism yields no barriers in excess of the reaction endothermicity. Since this reaction has been proposed as one of the steps in interstellar synthesis of glycine, the simplest amino acid, this work suggests that the formation of amino acids is occurring in and/or on interstellar water ice grains, and not in the gas phase.

Walch, Stephen P.↗

NASA NTRS1998

Computational Studies of the Interaction of H/H2 with Diamond and Silicon Surfaces

The interaction of hydrogen atoms and molecules with diamond and silicon surfaces is important in several important applications. Two areas that we are interested are: 1) tribology (molecular level friction) and 2) the role of H atoms in silicon chemical vapor deposition (CVD). In the tribology area, H atoms can be used to tie off dangling bonds, which otherwise form bonds between adjacent surfaces, and lead to resistance to sliding the surfaces by each other. Processes which are important in understanding molecular level friction include barriers to addition of H/H2 to the surface and barriers to migration of H atoms on the surface. In the silicon CVD area, we have studied the process of H2 elimination from the 100 surface of silicon. Cluster models for the dime surfaces of diamond are presented. The unrelaxed 100 surface has carbene like surface carbon atoms; however, for the relaxed surface these dimerize to give rows of surface dimers and there is a significant amount of p bonding between the radical orbitals of the dimer. The 110 surface has zig-zag rows of carbon atoms with a dangling bond on each carbon atom. These dangling bonds are hybridized away from each other and thus interact less strongly than for the 100 surface. Finally, the 111 surface has surface C atoms arranged in a triangular pattern and the surface dangling bonds are well separated from each other (second nearest neighbor distance) leading to almost no interaction between adjacent dangling bonds. These qualitative features may be quantified by computing the overlap of adjacent dangling bonds in a GVB(pp) calculation. The overlaps are 0.462, 0.292, and 0.016 for the diamond 100, 110, and 111 surfaces, respectively.

Walch, Stephen P.↗

NASA NTRS1997

The Role of Conical Intersections in Electronic Quenching of a State OH by H2 and N2

Lester and coworkers have experimentally characterized complexes of OH in the X and A states with H2 and N2. Recently, we have carried out ab initio calculations of relevant portions of the ground state and excited state potential energy surfaces for these systems, including the conical intersection regions, which are responsible for electronic quenching of the A state of OH by H2 and N2. Both of these systems have weakly bound complexes in the X state and strongly bound complexes in the A state. The OH-H2 complex is T-shaped, while the OH-N2 complex is collinear. In both cases the H end of OH is oriented toward the H2 or N2 molecule, respectively. Rotation of the OH about its center of mass involves only a small barrier and rotation by 1800 so that the 0 end of OH is oriented toward the H2 or N2 molecule leads to conical intersections with the ground state surface. Since there is about 95 kcal/mol of available energy after crossing to the ground state surface, chemical reactions on the ground state surface are possible, in addition to electronic quenching. In the case of OH-H2, the conical intersection is much lower in energy than the OH (A state) + H2 asymptotic energy and the A state complex can be characterized as a hanging well on the upper cone of the conical intersection. Passage through the conical intersection places the system on the ground state potential energy surface with the possibility of going to OH (X state) + H2 (i.e. electronic quenching) or to H2O + H (reaction). For OH-N2, the A state complex is also a hanging well on the upper cone of the conical intersection, but the conical intersection is only slightly below the OH (A state) + H2 asymptotic energy. Passage through the conical intersection can lead to OH (X state) + N2 (i.e. electronic quenching), but so far reactive pathways have not been found.

Walch, Stephen P.↗

NASA NTRS1996

Theoretical Studies of the HO/HO2 Catalytic Cycle for Ozone Destruction

Recently it has been determined that the HO/HO2 catalytic cycle accounts for nearly one-half of the total ozone depletion in the lower stratosphere. The catalytic cycle is: (1) HO + O3 yields HO2 + O2; (2) HO2 + O3 yields HO + O2 + O2. The net reaction is 2O3 yields 3O2. The rate limiting step in this process is the reaction of HO2 with ozone. There is a problem extending the experimental measurement of the rate of this reaction over the range 233-400 K down to stratospheric temperatures of 210-220 K. Therefore we have undertaken a project to determine the temperature dependence of the rate constant for this reaction in the low temperature region. The first step in this project, which is described in this poster, is the determination of the relevant potential energy surfaces. The calculations use CASSCF/derivative methods to define the pathways followed by CASSCF/ACPF to determine the energetics. The HO + O3 reaction is found to proceed through an HO4 complex, which is unstable with respect to HO2 + O2. The HO2 +O3 reaction is more complex. One pathway, which has been characterized, is the formation of an HO5 complex which decomposes to HO3 + O2 and subsequently to HO + O2 + O2. Another pathway, which is believed to also play a role, is hydrogen abstraction to give O2 + HO3 and subsequent decomposition of HO3 to HO + O2. Isotopic labeling experiments indicate that the later pathway is dominant. However, so far attempts to locate the saddle point for this pathway have not been successful. We have also characterized the potential energy surfaces for a number of species involved in these reactions, including HO3 and triplet O4. The triplet O4 species is probably involved in the reaction of vibrationally excited O2 with ground state O2 leading to O3 + O. The latter reaction is believed to be important as an additional source of stratospheric ozone.

Walch, Stephen P.↗

NASA NTRS1996

Theoretical Studies of Some HEDM Species: Cyclic O4, Cyclic O3 and Cubane

Calculations have been carried out for the HEDM species (cyclic O4, cyclic O3, and cubane) using CASSCF/derivative and CASSCF/ICCI methods. Cyclic O4 is of interest both as a potential HEDM species and because of its possible role in the ozone deficit problem in atmospheric chemistry. We have studied the pathway for decomposition from the D(2d) minimum and also have found the approximate location of the singlet triplet crossing. The barrier to decomposition is found to be about 9 kcal/mol and is not limited by the singlet triplet crossing. For cyclic O3 we have focused on the crossings between the lowest five surfaces (X(1)A(1), s(1)A(1), (1)A(2), (1)B(1), and (1)B(2)) to provide some insight into ways to form cyclic O3 photochemically. The crossing region between the X(1)A(1) and 2(1)A(1) surfaces is in agreement with the work of Xantheas et al. The calculations show that vertical excitation from the ground state to the (1)A(2) state leads to a crossing with the (1)A(1) manifold near the crossing region of the X(1)A(1) and 2(1)A(1) surfaces. We have studied the decomposition pathways for cubane to benzene plus acetylene and to cyclooctatetraene. We have also studied the ground and excited states for the photochemical ring closure step. The state which closes to cubane can be described as a double triplet pi to pi* excitation with respect to the ground state. Thus, this state has only a small oscillator strength with respect to the ground state. However, there is a singlet pi to pi* state at nearly the same energy and excitation to this state followed by intersystem crossing could lead to the triplet pi to pi* state.

Walch, Stephen P.↗

NASA NTRS1995

Characterization of the Minimum Energy Paths for the Ring Closure Reactions of C4H3 with Acetylene

The ring closure reaction of C4H3 with acetylene to give phenyl radical is one proposed mechanism for the formation of the first aromatic ring in hydrocarbon combustion. There are two low-lying isomers of C4H3; 1-dehydro-buta-l-ene-3-yne (n-C4H3) and 2-dehydro-buta-l-ene-3-yne (iso-C4H3). It has been proposed that only n-C4H3 reacts with acetylene to give phenyl radical, and since iso-C4H3 is more stable than n-C4H3, formation of phenyl radical by this mechanism is unlikely. We report restricted Hartree-Fock (RHF) plus singles and doubles configuration interaction calculations with a Davidson's correction (RHF+1+2+Q) using the Dunning correlation consistent polarized valence double zeta basis set (cc-pVDZ) for stationary point structures along the reaction pathway for the reactions of n-C4H3 and iso-C4H3 with acetylene. n-C4H3 plus acetylene (9.4) has a small entrance channel barrier (17.7) (all energetics in parentheses are in kcal/mol with respect to iso-C4H3 plus acetylene) and the subsequent closure steps leading to phenyl radical (-91.9) are downhill with respect to the entrance channel barrier. Iso-C4H3 Plus acetylene also has an entrance channel barrier (14.9) and there is a downhill pathway to 1-dehydro-fulvene (-55.0). 1-dehydro-fulvene can rearrange to 6-dehydro-fulvene (-60.3) by a 1,3-hydrogen shift over a barrier (4.0), which is still below the entrance channel barrier, from which rearrangement to phenyl radical can occur by a downhill pathway. Thus, both n-C4H3 and iso-C4H3 can react with acetylene to give phenyl radical with small barriers.

Walch, Stephen P.↗

NASA NTRS1995

Characterization of the Minimum Energy Paths and Energetics for the reaction of Vinylidene with Acetylene

The reaction of vinylidene (CH2C) with acetylene may be an initiating reaction in soot formation. We report minimum energy paths and accurate energetics for a pathway leading to vinylacetylene and for a number of isomers Of C4H4. The calculations use complete active space self-consistent field (CASSCF) derivative methods to characterize the stationary points and internally contacted configuration interaction (ICCI) and/or coupled cluster singles and doubles with a perturbational estimate of triple excitations (CCSD(T)) to determine the energetics. We find an entrance channel barrier of about 5 kcal/mol for the addition of vinylidene to acetylene, but no barriers above reactants for the reaction pathway leading to vinylacetylene.

Walch, Stephen P.↗

NASA NTRS1995

Characterization of the Minimum Energy Path for the Reaction of Singlet Methylene with N2: The Role of Singlet Methylene in Prompt NO

We report calculations of the minimum energy pathways connecting CH2 + N2 to diazomethane and diazirine, for the rearrangement of diazirine to diazomethane, for the dissociation of diazirine to HCN2+H, and of diazomethane to CH2N+N. The calculations use Complete Active Space Self-Consistent Field (CASSCF) derivative methods to characterize the stationary points and Internally Contracted Configuration Interaction (ICCI) to determine the energetics. The calculations suggest a potential new source of prompt NO from the reaction CH2 with N2 to give diazirine, and subsequent reaction of diazirine with hydrogen abstracters to form doublet HCN2, which leads to HCN+N(S-4) on the previously studied CH+N2 surface. The calculations also predict accurate 0 K heats of formation of 77.7 kcal/mol and 68.0 kcal/mol for diazirine and diazomethane, respectively.

Walch, Stephen P.↗

NASA NTRS1995

Characterization of the Minimum Energy Path for the Reaction of Singlet Methylene with N2: The Role of Singlet Methylene in Prompt NO

We report calculations of the minimum energy pathways connecting (1)CH2+N2 to diazomethane and diazirine, for the rearrangement of diazirine to diazomethane, for the dissociation of diazirine to HCN2+H, and of diazomethane to CH2N+N. The calculations use complete active space self-consistent field (CASSCF) derivative methods to characterize the stationary points and internally contracted configuration interaction (ICCI) to determine the energetics. The calculations suggest a potential new source of prompt NO from the reaction of (1)CH2 with N2 to give diazirine, and subsequent reaction of diazirine with hydrogen abstracters to form doublet HCN2, which leads to HCN+N(S-4) on the previously studied CH+N2 Surface. The calculations also predict accurate 0 K heats of formation of 77.7 kcal/mol and 68.0 kcal/mol for diazirine and diazomethane, respectively.

Walch, Stephen P.↗

NASA NTRS1995

A New Potential Energy Surface for N+O2: Is There an NOO Minimum?

We report a new calculation of the N+02 potential energy surface using complete active space self-consistent field internally contracted configuration interaction with the Dunning correlation consistent basis sets. The peroxy isomer of N02 is found to be a very shallow minimum separated from NO+O by a barrier of only 0.3 kcal/mol (excluding zero-point effects). The entrance channel barrier height is estimated to be 8.6 kcal/mol for ICCI+Q calculations correlating all but the Ols and N1s electrons with a cc-p VQZ basis set.

Walch, Stephen P.↗

NASA NTRS1995

Characterization of the Minimum Energy Paths and Energetics for the Reaction of Vinylidene with Acetylene

The reaction of vinylidene (CH2C) with acetylene may be an initiating reaction in soot formation. We report minimum energy paths and accurate energetics for a pathway leading to vinyl-acetylene and for a number of isomers of C4H4. The calculations use complete active space self-consistent field (CASSCF) derivative methods to characterize the stationary points and internally contacted configuration interaction (ICCI) and/or coupled cluster singles and doubles with a perturbational estimate of triple excitations (CCSD(T)) to determine the energetics. We find an entrance channel barrier of about 5 kcal/mol for the addition of vinylidene to acetylene, but no barriers above reactants for the reaction pathway leading to vinyl-acetylene.

Walch, Stephen P.↗

NASA NTRS1994

Computed potential energy surfaces for chemical reactions

Quantum mechanical methods have been used to compute potential energy surfaces for chemical reactions. The reactions studied were among those believed to be important to the NASP and HSR programs and included the recombination of two H atoms with several different third bodies; the reactions in the thermal Zeldovich mechanism; the reactions of H atom with O2, N2, and NO; reactions involved in the thermal De-NO(x) process; and the reaction of CH(squared Pi) with N2 (leading to 'prompt NO'). These potential energy surfaces have been used to compute reaction rate constants and rates of unimolecular decomposition. An additional application was the calculation of transport properties of gases using a semiclassical approximation (and in the case of interactions involving hydrogen inclusion of quantum mechanical effects).

Walch, Stephen 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.↗