Engineering Papers⌕ Search

Engineering topics

Walch, Stephen P.

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

At least 37 records · Page 2

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.↗

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.↗

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.↗

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.↗

Computed potential energy surfaces for chemical reactions

A new global potential energy surface (PES) is being generated for O(P-3) + H2 yields OH + H. This surface is being fit using the rotated Morse oscillator method, which was used to fit the previous POL-CI surface. The new surface is expected to be more accurate and also includes a much more complete sampling of bent geometries. A new study has been undertaken of the reaction N + O2 yields NO + O. The new studies have focused on the region of the surface near a possible minimum corresponding to the peroxy form of NOO. A large portion of the PES for this second reaction has been mapped out. Since state to state cross sections for the reaction are important in the chemistry of high temperature air, these studies will probably be extended to permit generation of a new global potential for reaction.

Walch, Stephen P.↗

Theoretical characterization of the reaction NH2+NO - products

The potential energy surface (PES) for NH2+NO is characterized using complete active space self-consistent-field (CASSCF)/derivative methods to locate the stationary points, followed by internally contracted configuration interaction (ICCI) calculations to determine the energetics. Production of N2+H2O is found to involve a complex mechanism, which, however, has no barrier with respect to NH2+NO. This pathway is exothermic by 124.5 kcal/mol. Production of NH2+OH can occur with no barrier other than the exothermicity from any of three isomers of HNNOH.

Walch, Stephen P.↗

Theoretical characterization of the reaction NH2 + O yields products

The potential energy surface for NH2+O has been characterized using complete active space self-consistent field (CASSCF)/derivative calculations to determine stationary point geometries and frequencies followed by internally contracted configuration interaction (ICCI) calculations to determine the energetics. The calculations predict a NO bond strength of 85.8 kcal/mol for NH2O. The barrier for isomerization of NH2O to trans-HNOH is predicted to be 48.0 kcal/mol and the barriers for H+HNO forming NH2O and NHOH are predicted to be 2.1 and 8.3 kcal/mol, respectively (all corrected for zero-point energy). The computed heats of formation for NH2O and cis- and trans-HNOH are in good agreement with the present results. The barrier for H + HNO yields H2 + NO is computed to be about 0.3 kcal/mol.

Walch, Stephen P.↗

A global potential energy surface for ArH2

We describe a simple analytic representation of the ArH2 potential energy surface which well reproduces the results of extensive ab initio electronic structure calculations. The analytic representation smoothly interpolates between the dissociated H2 and strong bonding limits. In the fitting process, emphasis is made on accurately reproducing regions of the potential expected to be important for high temperature (ca. 3000 K) collision processes. Overall, the anisotropy and H2 bond length dependence of the analytic representation well reproduce the input data.

Schwenke, David W.↗

Computed barrier heights for H + CH2O yields CH3O yields CH2OH

The barrier heights (including zero-point effects) for H + CH2O yields CH3O and CH3O yields CH2OH have been computed using complete active space self consistent field (CASSCF)/gradient calculations to define the stationary point geometries and harmonic frequencies and internally contracted configuration-interaction (CCI) to refine the energetics. The computed barrier heights are 5.6 kcal/mol and 30.1 kcal/mol, respectively. The former barrier height compares favorably to an experimental activation energy of 5.2 kcal/mol.

Walch, Stephen P.↗

Theoretical characterization of the reaction CH3 + OH yields CH3OH yields products - The (1)CH2 + H2O, H2 + HCOH, and H2 + H2CO channels

The potential energy surface (PES) for the CH3OH system has been characterized for the (1)CH2 + H2O, H2 + HCOH, and H2 + H2CO product channels using complete-active-space self-consistent-field (CASSCF) gradient calculations to determine the stationary point geometries and frequencies followed by CASSCF/internally contracted configuration-interaction (CCI) calculations to refine the energetics. The (1)CH2 + H2O channel is found to have no barrier. The long range interaction is dominated by the dipole-dipole term, which orients the respective dipole moments parallel to each other but pointing in opposite directions. At shorter separations there is a dative bond structure in which a water lone pair donates into the empty 'a' orbital of CH2. Subsequent insertion of CH2 into an OH bond of water have barriers located at -5.2 kcal/mol and 1.7 kcal/mol, respectively, with respect to CH3 + OH. From comparison of the computed energetics of the reactants and products to known thermochemical data it is estimated that the computed PES is accurate to plus or minus 2 kcal/mol.

Walch, Stephen P.↗

Computed Potential Energy Surfaces for Chemical Reactions

A manuscript describing the calculations on the (1)CH2 + H2O, H2 + HCOH, and H2 + H2CO product channels in the CH3 + OH reaction, which were described in the last progress report, has been accepted for publication in J. Chem. Phys., and a copy of the manuscript is included in the appendix. The production of (1)CH2 in this reaction is important in hydrocarbon combustion since (1)CH2 is highly reactive and would be expected to insert into N2, possibly leading to a new source for prompt NO(x) (vide infra). During the last six months new calculations have been carried out for the NH2 + NO system, which is important in the thermal de-NO(x) process.

Heinemann, K.↗

Theoretical characterization of the reaction CH3 +OH yields CH3OH yeilds products: The (1)CH2 + H2O, H2 + HCOH, and H2 + H2CO channels

The potential energy surface (PES) for the CH3OH system has been characterized for the (1)CH2 + H2O, H2 + HCOH, and H2 + H2CO product channels using complete-active-space self-consistent-field (CASSCF) gradient calculations to determine the stationary point geometries and frequencies followed by CASSCF/internally contracted configuration-interaction (CCI) calculations to refine the energetics. The (1)CH2 + H2O channel is found to have no barrier. The long range interaction is dominated by the dipole-dipole term, which orients the respective dipole moments parallel to each other but pointing in opposite directions. At shorter separations there is a dative bond structure in which a water lone pair donates into the empty a" orbital of CH2. Subsequent insertion of CH2 into an OH bond of water have barriers located at -5.2 kcal/mol and 1.7 kcal/mol, respectively, with respect to CH3 + OH. From comparison of the computed energetics of the reactants and products to known thermochemical data it is estimated that the computed PES is accurate to plus or minus 2 kcal/mol.

Walch, Stephen P.↗

Computed barrier heights for H + CH2O yields CH3O yields CH2OH

The barrier heights (including zero-point effects) for H + CH2O yields CH3O and CH3O yields CH2OH have been computed using complete active space self consistent field (CASSCF)/gradient calculations to define the stationary point geometries and harmonic frequencies and internally contracted configuration-interaction (CCI) to refine the energetics. The computed barrier heights are 5.6 kcal/mol and 30.1 kcal/mol, respectively. The former barrier height compares favorably to an experimental activation energy of 5.2 kcal/mol.

Walch, Stephen P.↗

Theoretical characterization of the potential energy surface for NH + NO

The potential energy surface for NH + NO was characterized using complete active space self-consistent field (CASSCF) gradient calculation to determine the stationary point geometries and frequencies followed by CASSCF/internally contracted configuration interaction calculations to refine the energetics. The present results are in qualitative accord with the BAC-MP4 calculations, but there are differences as large as 8 kcal/mol in the detailed energetics.

Walch, Stephen P.↗

H-N2 interaction energies, transport cross sections, and collision integrals

The energies for the interaction of a hydrogen atom with a nitrogen molecule have been calculated for large separation distances using a complete-active-space self-consistent-field/externally contracted configuration interaction method. H-N2 transport cross sections and collision integrals have been calculated using sudden approximations and a semiclassical description of the scattering. The values of these quantities are found to be close to the corresponding values determined from the average (isotropic) potential energy. The collision integrals are applied to determine diffusion and viscosity coefficients; the theoretical diffusion agrees well with the measured data available from experiments at low temperatures.

Stallcop, James R.↗

Computed potential energy surfaces for chemical reactions

The work on the NH + NO system which was described in the last progress report was written up and a draft of the manuscript is included in the appendix. The appendix also contains a draft of a manuscript on an Ar + H + H surface. New work which was completed in the last six months includes the following: (1) calculations on the (1)CH2 + H2O, H2 + HCOH, and H2 + H2CO product channels in the CH3 + OH reaction; (2) calculations for the NH2 + O reaction; (3) calculations for the CH3 + O2 reaction; and (4) calculations for CH3O and the two decomposition channels--CH2OH and H + H2CO. Detailed descriptions of this work will be given in manuscripts; however, brief descriptions of the CH3 + OH and CH3 + O2 projects are given.

Heinemann, K.↗

Theoretical characterization of the potential energy surface for NH + NO

The potential energy surface (PES) for NH + NO was characterized using complete active space self-consistent field (CASSCF) gradient calculations to determine the stationary point geometries and frequencies followed by CASSCF/internally contracted configuration interaction (CCI) calculations to refine the energetics. The present results are in qualitative accord with the BAC-MP4 calculations, but there are differences as large as 8 kcal/mol in the detailed energetics. Addition of NH to NO on a (2)A' surface, which correlated with N2 + OH or H + N2O products, involves barriers of 3.2 kcal/mol (trans) and 6.3 kcal/mol (cis). Experimental evidence for these barriers is found in earlier works. The (2)A' surface has no barrier to addition, but does not correlate with products. Surface crossings between the barrierless (2)A' surface and the (2)A' surface may be important. Production of N2 + OH products is predicted to occur via a planar saddle point of (2)A' symmetry. This is in accord with the preferential formation of II(A') lambda doublet levels of OH in earlier experiments. Addition of NH (1)delta to NO is found to occur on an excited state surface and is predicted to lead to N2O product as observed in earlier works.

Walch, Stephen P.↗

A global potential energy surface for ArH2

We describe a simple analytic representation of the ArH2 potential energy surface which well reproduces the results of extensive ab initio electronic structure calculations. The analytic representation smoothly interpolates between the dissociated H2 and strong bonding limits. In the fitting process, emphasis is made on accurately reproducing regions of the potential expected to be important for high temperature (ca. 3000 K) collision processes. Overall, the anisotropy and H2 bond length dependence of the analytic representation well reproduce the input data.

Schwenke, David W.↗