Search NASA⌕ Search

Engineering topics

Kouri, Donald J.

Publications and source records attributed to Kouri, Donald J..

15 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

State-to-state rates for the D + H2(v = 1, j = 1) yield HD(v-prime, j-prime) + H reaction - Predictions and measurements

A fully quantal wavepacket approach to reactive scattering in which the best available H3 potential energy surface was used enabled a comparison with experimentally determined rates for the D + H2(v = 1, j = 1) yield HD(v-prime = 0, 1, 2; j-prime) + H reaction at significantly higher total energies (1.4 to 2.25 electron volts) than previously possible. The theoretical results are obtained over a sufficient range of conditions that a detailed simulation of the experiment was possible, thus making this a definitive comparison of experiment and theory. Good to excellent agreement is found for the vibrational branching ratios and for the rotational distributions within each product vibrational level. However, the calculated rotational distributions are slightly hotter than the experimentally measured ones.

Neuhauser, Daniel↗
NASA NTRS1991

Converged three-dimensional quantum mechanical reaction probabilities for the F + H2 reaction on a potential energy surface with realistic entrance and exit channels and comparisons to results for three other surfaces

Accurate three-dimensional quantum mechanical reaction probabilities are presented for the reaction F + H2 yields HF + H on the new global potential energy surface 5SEC for total angular momentum J = 0 over a range of translational energies from 0.15 to 4.6 kcal/mol. It is found that the v-prime = 3 HF vibrational product state has a threshold as low as for v-prime = 2.

Lynch, Gillian C.↗
NASA NTRS1991

A comparison of three time-dependent wave packet methods for calculating electron-atom elastic scattering cross sections

Three time-dependent wave packet methods for performing elastic scattering calculations from screened Coulomb potentials are compared. The three methods are the time-dependent amplitude density method (TDADM), a Cayley-transform method (CTM), and the Chebyshev propagation method of Tal-Ezer and Kosloff. Both the TDADM and the CTM are based on a time-dependent integral equation for the wave function. In the first, the time-dependent amplitude density is propagated, while in the other two, the wave function is propagated. As a numerical example, phase shifts and cross sections are calculated using a screened Coulomb, Yukawa type potential over the range 200-1000 eV. It is found that, in most cases, all three methods yield comparable accuracy and are about equally efficient computationally. However for l = 0, where the Coulomb well is not screened by the centrifugal potential, the TDADM requires smaller grid spacings to maintain accuracy.

Judson, Richard S.↗
NASA NTRS1991

Total integral reactive cross sections for F + H2 yielding HF + H - Comparison of converged quantum, quasiclassical trajectory and experimental results

The paper reports converged quantum total integral reactive cross sections for the reaction F + H2 yielding HF + H, for initial rotational states j sub i = 0 and 1, using a time-dependent method. The results are compared to classical results and to the experimental results of Neumark et al. (1985). Strong quantum effects are found in the threshold region for both initial states (i.e., in the dependence of the reaction on initial state for low energies). The classical results agree better with experiment than do the quantum results; this appears to be due to errors in the potential used.

Neuhauser, Daniel↗
NASA NTRS1991

Theoretical studies of molecular collisions

The following subject areas are covered: (1) total integral reactive cross sections and vibrationally resolved reaction probabilities for F + H2 = HF + H; (2) a theoretical study of inelastic O + N2 collisions; (3) body frame close coupling wave packet approach to gas phase atom-rigit rotor inelastic collisions; (4) wave packet study of gas phase atom-rigit motor scattering; (5) the application of optical potentials for reactive scattering; (6) time dependent, three dimensional body frame quantal wave packet treatment of the H + H2 exchange reaction; (7) a time dependent wave packet approach to atom-diatom reactive collision probabilities; (8) time dependent wave packet for the complete determination of s-matrix elements for reactive molecular collisions in three dimensions; (9) a comparison of three time dependent wave packet methods for calculating electron-atom elastic scattering cross sections; and (10) a numerically exact full wave packet approach to molecule-surface scattering.

Kouri, Donald J.↗
NASA NTRS1991

The application of time-dependent wavepacket methods to reactive scattering

Several methods for performing numerically exact reactive scattering calculations using time-dependent wavepackets are reviewed. The basic idea is to take the multiarrangement reactive problem and reformulate it as one or more inelastic ones. In the simplest method, total reaction probabilities are extracted by calculating the flux of the wavepacket as it leaves the interaction region in the direction of the reactive arrangement. To make this practical, complex potentials that absorb the wavepacket before it reaches the numerical grid boundary are used. Methods that generate observables ranging from total, energy-averaged reaction probabilities up to energy- and state-resolved S-matrix elements are used. Techniques for efficiently performing the necessary inelastic wavepacket propagation are also reviewed.

Neuhauser, Daniel↗
NASA NTRS1990

Time-dependent treatment of scattering. II - Novel integral equation approach to quantum wave packets

The novel wave-packet propagation scheme presented is based on the time-dependent form of the Lippman-Schwinger integral equation and does not require extensive matrix inversions, thereby facilitating application to systems in which some degrees of freedom express the potential in a basis expansion. The matrix to be inverted is a function of the kinetic energy operator, and is accordingly diagonal in a Bessel function basis set. Transition amplitudes for various orbital angular momentum quantum numbers are obtainable via either Fourier transform of the amplitude density from the time to the energy domain, or the direct analysis of the scattered wave packet.

Sharafeddin, Omar A.↗
NASA NTRS1990

The application of optical potentials for reactive scattering - A case study

A precise method is presented for measuring decay fractions of metastable cluster ions with corrections concerning instrumental artifacts and ion trajectory of parents and daughters. Experimental data are used to derive the Gspann parameter and heat capacity of clusters as described in Klots' evaporative ensemble model. The model is applied to obtain binding energies of ammonia cluster ions. The deduced binding energy values are in very good agreement with both thermochemical data and Engelking's (1986) modified statistical theory.

Neuhauser, Daniel↗
NASA NTRS1990

A time-dependent wave packet approach to atom-diatom reactive collision probabilities - Theory and application to the H + H2(J = 0) system

This paper describes a new approach to the study of atom-diatom reactive collisions in three dimensions employing wave packets and the time-dependent Schroedinger equation. The method uses a projection operator approach to couple the inelastic and reactive portions of the total wave function and optical potentials to circumvent the necessity of using product arrangement coordinates. Reactive transition probabilities are calculated from the state resolved flux of the wave packet as it leaves the interaction region in the direction of the reactive arrangement channel. The present approach is used to obtain such vibrationally resolved probabilities for the three-dimensional H + H2 (J = 0) hydrogen exchange reaction, using a body-fixed system of coordinates.

Neuhauser, Daniel↗
NASA NTRS1990

Time-dependent wave-packet method for the complete determination of S-matrix elements for reactive molecular collisions in three dimensions

An alternative time-dependent wave-packet method for treating three-dimensional gas phase reactive atom-diatom collisions is presented. The method employs a nonreactive body-frame wave packet propagation procedure, made possible by judicious use of absorbing optical potentials, a novel scheme for interpolating the wave function from coordinates in one arrangement to those in another and the fact that the time-dependent Schroedinger equation is an initial-value problem. The last feature makes possible a computationally viable and accurate procedure for changing from one arrangement's coordinates to another. In addition, the method allows the determination of S-matrix elements over a wide range of energies from a single wave-packet propagation. The method is illustrated by carrying out detailed calculations of inelastic and reactive scattering in the H + H2 system using the Liu-Siegbahn-Truhlar-Horowitz potential surface.

Judson, Richard S.↗
NASA NTRS1990

Dynamical basis sets for algebraic variational calculations in quantum-mechanical scattering theory

New basis sets are proposed for linear algebraic variational calculations of transition amplitudes in quantum-mechanical scattering problems. These basis sets are hybrids of those that yield the Kohn variational principle (KVP) and those that yield the generalized Newton variational principle (GNVP) when substituted in Schlessinger's stationary expression for the T operator. Trial calculations show that efficiencies almost as great as that of the GNVP and much greater than the KVP can be obtained, even for basis sets with the majority of the members independent of energy.

Sun, Yan↗
NASA NTRS1990

Time-dependent treatment of scattering - Integral equation approaches using the time-dependent amplitude density

The time-dependent form of the Lippmann-Schwinger integral equation is used as the basis of several new wave packet propagation schemes. These can be formulated in terms of either the time-dependent wave function or a time-dependent amplitude density. The latter is nonzero only in the region of configuratiaon space for which the potential is nonzero, thereby in principle obviating the necessity of large grids or the use of complex absorbing potentials when resonances cause long collision times (leading, consequently, to long propagation times). Transition amplitudes are obtained in terms of Fourier transforms of the amplitude density from the time to the energy domain. The approach is illustrated by an application to a standard potential scattering model problem where, as in previous studies, the action of the kinetic energy operator is evaluated by fast Fourier transform (FFT) techniques.

Hoffman, David K.↗
NASA NTRS1989

Direct calculation of the reactive transition matrix by L-squared quantum mechanical variational methods with complex boundary conditions

A new formalism of the generalized Newton variational principle for the calculation of quantum mechanical state-to-state reaction probabilities is presented. The reformulation involves solving directly for the transition matrix rather than the reactance mtrix so that calculations may be carried out for individual columns of the transition matrix without obtaining solutions for all possible initial channels. The convergence of calculations with real and complex boundary conditions are compared for H + H2 - H2 + H, O + H2 - OH + H, and O + HD - OH + D and OD + H.

Sun, Yan↗
NASA NTRS1989

Time dependent three-dimensional body frame quantal wave packet treatment of the H + H2 exchange reaction on the Liu-Siegbahn-Truhlar-Horowitz (LSTH) surface

The first successful application of the three-dimensional quantum body frame wave packet approach to reactive scattering is reported for the H + H2 exchange reaction on the LSTH potential surface. The method used is based on a procedure for calculating total reaction probabilities from wave packets. It is found that converged, vibrationally resolved reactive probabilities can be calculated with a grid that is not much larger than required for the pure inelastic calculation. Tabular results are presented for several energies.

Neuhauser, Daniel↗