Search NASA⌕ Search

SEARCH · Search NASA

Results for “quasiclassical trajectory calculations”

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.

On the Development of a New Nonequilibrium Chemistry Model for Mars Entry

This paper represents a summary of results to date of an on-going effort at NASA Ames Research Center to develop a physics-based non-equilibrium model for hypersonic entry into the Martian atmosphere. Our approach is to first compute potential energy surfaces based on accurate solutions of the electronic Schroedinger equation and then use quasiclassical trajectory calculations to obtain reaction cross sections and rate coefficients based on these potentials. We have presented new rate coefficients for N2 dissociation and CO dissociation and exchange reactions. These results illustrate shortcomings with some of the rate coefficients in Parks original T-Tv model for Mars entries and with some of the 30-45 year old shock tube data. We observe that the shock tube experiments of CO + O dissociation did not adequately account for the exchange reaction that leads to formation of C + O2. This reaction is actually the primary channel for CO removal in the shock layer at temperatures below 10,000 K, because the reaction enthalpy for exchange is considerably lower than the comparable value for dissociation.

quasiclassical trajectory calculation↗

Calculations of rate constants for the three-body recombination of H2 in the presence of H2

A new global potential energy hypersurface for H2 + H2 is constructed and quasiclassical trajectory calculations performed using the resonance complex theory and energy transfer mechanism to estimate the rate of three body recombination over the temperature range 100 to 5000 K. The new potential is a faithful representation of ab initio electron structure calculations, is unchanged under the operation of exchanging H atoms, and reproduces the accurate H3 potential as one H atom is pulled away. Included in the fitting procedure are geometries expected to be important when one H2 is near or above the dissociation limit. The dynamics calculations explicitly include the motion of all four atoms and are performed efficiently using a vectorized variable-stepsize integrator. The predicted rate constants are approximately a factor of two smaller than experimental estimates over a broad temperature range.

Schwenke, David W.↗

Calculations of rate constants for the three-body recombination of H2 in the presence of H2

A new global potential energy hypersurface for H2 + H2 is constructed and quasiclassical trajectory calculations performed using the resonance complex theory and energy transfer mechanism to estimate the rate of three body recombination over the temperature range 100 to 5000 K. The new potential is a faithful representation of ab initio electron structure calculations, is unchanged under the operation of exchanging H atoms, and reproduces the accurate H3 potential as one H atom is pulled away. Included in the fitting procedure are geometries expected to be important when one H2 is near or above the dissociation limit. The dynamics calculations explicitly include the motion of all four atoms and are performed efficiently using a vectorized variable-stepsize integrator. The predicted rate constants are approximately a factor of two smaller than experimental estimates over a broad temperature range.

Schwenke, David W.↗

A New Global Potential Energy Surface for the Hydroperoxyl Radical, HO2: Reaction Coefficients for H + O2 and Vibrational Splittings for H Atom Transfer

A new analytic global potential energy surface describing the hydroperoxyl radical system H((sup 2)S) + O2(X (sup 3)Sigma((sup -)(sub g))) (reversible reaction) HO2 ((X-tilde) (sup 2)A'') (reversible reaction) O((sup 3)P) + O H (X (sup 2)Pi) has been fitted using the ab initio complete active space SCF (self-consistent-field)/externally contracted configuration interaction (CASSCF/CCI) energy calculations of Walch and Duchovic. Results of quasiclassical trajectory studies to determine the rate coefficients of the forward and reverse reactions at combustion temperatures will be presented. In addition, vibrational energy levels were calculated using the quantum DVR-DGB (discrete variable representation-distributed Gaussian basis) method and the splitting due to H atom migration is investigated. The material of the proposed presentation was reviewed and the technical content will not reveal any information not already in the public domain and will not give any foreign industry or government a competitive advantage.

Dateo, Christopher E.↗

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↗

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↗

A new semiclassical decoupling scheme for electronic transitions in molecular collisions - Application to vibrational-to-electronic energy transfer

A new semiclassical decoupling scheme (the trajectory-based decoupling scheme) is introduced in a computational study of vibrational-to-electronic energy transfer for a simple model system that simulates collinear atom-diatom collisions. The probability of energy transfer (P) is calculated quasiclassically using the new scheme as well as quantum mechanically as a function of the atomic electronic-energy separation (lambda), with overall good agreement between the two sets of results. Classical mechanics with the new decoupling scheme is found to be capable of predicting resonance behavior whereas an earlier decoupling scheme (the coordinate-based decoupling scheme) failed. Interference effects are not exhibited in P vs lambda results.

Lee, H.-W.↗