Engineering PapersSearch

SEARCH · Engineering Papers

Results for “potential energy surfaces”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

Theoretical characterization of the potential energy surface for the reversible reaction H + O2 yields HO2(asterisk) yields OH + O. III - Computed points to define a global potential energy surface

Computed energies and geometries are reported which, combined with previously published calculations, permit a global representation of the potential energy surface for the reaction H + O2 yields HO2(asterisk) yields OH + O. These new calculations characterize the potential energy surface (PES) for all H atom angles of approach to O2 and for the region of the inner repulsive wall. The region of the T-shaped H-O2 exchange saddle point is connected with the constrained energy minimum (CEM) path, and a new collinear H-O2 exchange saddle point is characterized which lies only 9 kcal/mol above the H + O2 asymptote. A vibrational analysis which utilizes local cubic and quartic polynomial representations of the PES along the CEM path has been carried out. Optimal geometries, energies, and harmonic frequencies are reported along with anharmonic analyses for the O2 and OH asymptotes and for the HO2 minimum region of the PES.

Walch, Stephen P.

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.

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.

Classical trajectory studies of gas phase reaction dynamics and kinetics using ab initio potential energy surfaces

Strategies for constructing global potential energy surfaces from a limited number of accurate ab initio electronic energy calculations are discussed. Generally, these data are concentrated in small regions of configuration space (e.g., in the vicinity of saddle points and energy minima) and difficulties arise in generating a potential function that is globally well-behaved. Efficient computer codes for carrying out classical trajectory calculations on vector and parallel processors are also described. Illustrations are given from recent work on the following chemical systems: Ca + HF yields CaF + H, H + H + H2 yields H2 + H2, N + O2 yields NO + O and O + N2 yields NO + N. The dynamics and kinetics of metathesis, dissociation, recombination, energy transfer and complex formation processes will be discussed.

Jaffe, Richard L.

An Ab Initio Based Potential Energy Surface for Water

We report a new determination of the water potential energy surface. A high quality ab initio potential energy surface (PES) and dipole moment function of water have been computed. This PES is empirically adjusted to improve the agreement between the computed line positions and those from the HITRAN 92 data base. The adjustment is small, nonetheless including an estimate of core (oxygen 1s) electron correlation greatly improves the agreement with experiment. Of the 27,245 assigned transitions in the HITRAN 92 data base for H2(O-16), the overall root mean square (rms) deviation between the computed and observed line positions is 0.125/cm. However the deviations do not correspond to a normal distribution: 69% of the lines have errors less than 0.05/cm. Overall, the agreement between the line intensities computed in the present work and those contained in the data base is quite good, however there are a significant number of line strengths which differ greatly.

Partridge, Harry

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

Ab Initio Potential Energy Surface for H-H2

Ab initio calculations employing large basis sets are performed to determine an accurate potential energy surface for H-H2 interactions for a broad range of separation distances. At large distances, the spherically averaged potential determined from the calculated energies agrees well with the corresponding results determined from dispersion coefficients; the van der Waals well depth is predicted to be 75 +/- 3 micro E(h). Large basis sets have also been applied to reexamine the accuracy of theoretical repulsive potential energy surfaces (25-70 kcal/mol above the H-H2 asymptote) at small interatomic separations; the Boothroyd, Keogh, Martin, and Peterson (BKMP) potential energy surface is found to agree with results of the present calculations within the expected uncertainty (+/- 1 kcal/mol) of the fit. Multipolar expansions of the computed H-H2 potential energy surface are reported for four internuclear separation distances (1.2, 1.401, 1.449, and 1.7a(0)) of the hydrogen molecule. The differential elastic scattering cross section calculated from the present results is compared with the measurements from a crossed beam experiment.

Patridge, Harry

Computed potential energy surfaces for chemical reactions

The objective was to obtain accurate potential energy surfaces (PES's) for a number of reactions which are important in the H/N/O combustion process. The interest in this is centered around the design of the SCRAM jet engine for the National Aerospace Plane (NASP), which was envisioned as an air-breathing hydrogen-burning vehicle capable of reaching velocities as large as Mach 25. Preliminary studies indicated that the supersonic flow in the combustor region of the scram jet engine required accurate reaction rate data for reactions in the H/N/O system, some of which was not readily available from experiment. The most important class of combustion reactions from the standpoint of the NASP project are radical recombinaton reactions, since these reactions result in most of the heat release in the combustion process. Theoretical characterizations of the potential energy surfaces for these reactions are presented and discussed.

Walch, Stephen P.

Theoretical characterization of the potential energy surface for H + N2 yields HN2. II - Computed points to define a global potential

A previous calculation for H + N2 (Walch et al., 1989) focused on the minimum energy path (MEP) region of the potential energy surface and on estimates of the lifetime of the HN2 species. In this paper, energies computed at geometries selected to permit a global representation of the potential energy surface (PES) are reported. As in the previous work, the calculations were performed using the complete active space self-consistent field/externally contracted configuration interaction method. The surface was characterized using the same basis set as in the previous paper except that an improved contraction of the H s-basis is used. Calculations with a larger basis set were carried out along an approximate MEP obtained with the smaller basis set. The new PES exhibits a sharp curvature, which was not present in the previous calculations, and has a slightly narrower and smaller barrier to dissociation. Saddle points for H atom exchange via collinear and T shaped HN2 complexes are also reported.

Walch, Stephen P.

Ab initio Potential Energy Surface for H-H2

Ab initio calculations employing large basis sets are performed to determine an accurate potential energy surface for H-H2 interactions for a broad range of separation distances. At large distances, the spherically averaged potential determined from the calculated energies agrees well with the corresponding results determined from dispersion coefficients; the van der Waals well depth is predicted to be 75 +/- (mu)E(sub h). Large basis sets have also been applied to reexamine the accuracy of theoretical repulsive potential energy surfaces. Multipolar expansions of the computed H-H2 potential energy surface are reported for four internuclear separation distances (1.2, 1.401, 1.449, and 1.7a(sub 0) of the hydrogen molecule. The differential elastic scattering cross section calculated from the present results is compared with the measurements from a crossed beam experiment.

Partridge, Harry

MCSCF potential energy surface for photodissociation of formaldehyde

The ground state potential energy surface for the dissociation of formaldehyde (H2CO to H2 and CO) is calculated with the ab initio MCSCF method with an extended (4-31G) basis set. The location, barrier height, and force constants of the transition state are determined, and the normal coordinate analysis is carried out. The calculated barrier height is 4.5 eV. Based on the calculated quantities, the detailed mechanism of the photochemical dissociation is discussed.

Jaffe, R. L.

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.

Collisional Dissociation of CO: ab initio Potential Energy Surfaces and Quasiclassical Trajectory Rate Coefficients

We have generated accurate global potential energy surfaces for CO+Ar and CO+O that correlate with atom-diatom pairs in their ground electronic states based on extensive ab initio electronic structure calculations and used these potentials in quasi-classical trajectory nuclear dynamics calculations to predict the thermal dissociation rate coefficients over 5000- 35000 K. Our results are not compatible with the 20-45 year old experimental results. For CO + Ar we obtain fairly good agreement with the experimental rate coefficients of Appleton et al. (1970) and Mick and Roth (1993), but our computed rate coefficients exhibit a stronger temperature dependence. For CO + O our dissociation rate coefficient is in close agreement with the value from the Park model, which is an empirical adjustment of older experimental results. However, we find the rate coefficient for CO + O is only 1.5 to 3.3 times larger than CO + Ar over the temperature range of the shock tube experiments (8000-15,000 K). The previously accepted value for this rate coefficient ratio is 15, independent of temperature. We also computed the rate coefficient for the CO + O ex- change reaction which forms C + O2. We find this reaction is much faster than previously believed and is the dominant process in the removal of CO at temperatures up to 16,000 K. As a result, the dissociation of CO is accomplished in two steps (react to form C+O2 and then O2 dissociates) that are endothermic by 6.1 and 5.1 eV, instead of one step that requires 11.2 eV to break the CO bond.

exchange rate coefficient

Comparison of Quantum Mechanical and Empirical Potential Energy Surfaces and Computed Rate Coefficients for N2 Dissociation

Comparisons are made between potential energy surfaces (PES) for N2 + N and N2 + N2 collisions and between rate coefficients for N2 dissociation that were computed using the quasiclassical trajectory method (QCT) on these PESs. For N2 + N we compare the Laganas empirical LEPS surface with one from NASA Ames Research Center based on ab initio quantum chemistry calculations. For N2 + N2 we compare two ab initio PESs (from NASA Ames and from the University of Minnesota). These use different methods for computing the ground state electronic energy for N4, but give similar results. Thermal N2 dissociation rate coefficients, for the 10,000K-30,000K temperature range, have been computed using each PES and the results are in excellent agreement. Quasi-stationary state (QSS) rate coefficients using both PESs have been computed at these temperatures using the Direct Molecular Simulation of Schwartentruber and coworkers. The QSS rate coefficients are up to a factor of 5 lower than the thermal ones and the thermal and QSS values bracket the results of shock-tube experiments. We conclude that the combination of ab initio quantum chemistry PESs and QCT calculations provides an attractive approach for the determination of accurate high-temperature rate coefficients for use in aerothermodynamics modeling.

quasiclassical trajectory method

A potential energy surface for the process H2 + H2O yielding H + H + H2O - Ab initio calculations and analytical representation

Extensive ab initio calculations on the ground state potential energy surface of H2 + H2O were performed using a large contracted Gaussian basis set and a high level of correlation treatment. An analytical representation of the potential energy surface was then obtained which reproduces the calculated energies with an overall root-mean-square error of only 0.64 mEh. The analytic representation explicitly includes all nine internal degrees of freedom and is also well behaved as the H2 dissociates; it thus can be used to study collision-induced dissociation or recombination of H2. The strategy used to minimize the number of energy calculations is discussed, as well as other advantages of the present method for determining the analytical representation.

Schwenke, David W.

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.