Search NASASearch

Engineering topics

Truhlar, Donald G.

Publications and source records attributed to Truhlar, Donald G..

Kinetics and Mechanism of the Singlet Oxygen Atom Reaction with Dimethyl Ether

Here, we combine in situ laser spectroscopy, quantum chemistry, and kinetic calculations to study the reaction of a singlet oxygen atom with dimethyl ether. Infrared laser absorption spectroscopy and Faraday rotation spectroscopy are used for the detection and quantification of the reaction products OH, H 2 O, HO 2 , and CH 2 O on submillisecond time scales. Fitting temporal profiles of products with simulations using an in-house reaction mechanism allows product branching to be quantified at 30, 60, and 150 Torr. The experimentally determined product branching agrees well with master equation calculations based on electronic structure data and transition state theory. The calculations demonstrate that the dimethyl peroxide (CH 3 OOCH 3 ) generated via O-insertion into the C–O bond undergoes subsequent dissociation to CH 3 O + CH 3 O through energetically favored reactions without an intrinsic barrier. This O-insertion mechanism can be important for understanding the fate of biofuels leaking into the atmosphere and for plasma-based biofuel processing technologies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

The quantum dynamics of electronically nonadiabatic chemical reactions

Considerable progress was achieved on the quantum mechanical treatment of electronically nonadiabatic collisions involving energy transfer and chemical reaction in the collision of an electronically excited atom with a molecule. In the first step, a new diabatic representation for the coupled potential energy surfaces was created. A two-state diabatic representation was developed which was designed to realistically reproduce the two lowest adiabatic states of the valence bond model and also to have the following three desirable features: (1) it is more economical to evaluate; (2) it is more portable; and (3) all spline fits are replaced by analytic functions. The new representation consists of a set of two coupled diabatic potential energy surfaces plus a coupling surface. It is suitable for dynamics calculations on both the electronic quenching and reaction processes in collisions of Na(3p2p) with H2. The new two-state representation was obtained by a three-step process from a modified eight-state diatomics-in-molecules (DIM) representation of Blais. The second step required the development of new dynamical methods. A formalism was developed for treating reactions with very general basis functions including electronically excited states. Our formalism is based on the generalized Newton, scattered wave, and outgoing wave variational principles that were used previously for reactive collisions on a single potential energy surface, and it incorporates three new features: (1) the basis functions include electronic degrees of freedom, as required to treat reactions involving electronic excitation and two or more coupled potential energy surfaces; (2) the primitive electronic basis is assumed to be diabatic, and it is not assumed that it diagonalizes the electronic Hamiltonian even asymptotically; and (3) contracted basis functions for vibrational-rotational-orbital degrees of freedom are included in a very general way, similar to previous prescriptions for locally adiabatic functions in various quantum scattering algorithms.

Truhlar, Donald G.

The nature and role of quantized transition states in the accurate quantum dynamics of the reaction O + H2 yields OH + H

Accurate quantum mechanical dynamics calculations are reported for the reaction probabilities of O(3P) + H2 yields OH + H with zero total angular momentum on a single potential energy surface. The results show that the reactive flux is gated by quantized transition states up to the highest energy studied, which corresponds to a total energy of 1.90 eV. The quantized transition states are assigned and compared to vibrationally adiabatic barrier maxima; their widths and transmission coefficients are determined; and they are classified as variational, supernumerary of the first kind, and supernumerary of the second kind. Their effects on state-selected and state-to-state reactivity are discussed in detail.

Chatfield, David C.

Complex generalized minimal residual algorithm for iterative solution of quantum-mechanical reactive scattering equations

Complex dense matrices corresponding to the D + H2 and O + HD reactions were solved using a complex generalized minimal residual (GMRes) algorithm described by Saad and Schultz (1986) and Saad (1990). To provide a test case with a different structure, the H + H2 system was also considered. It is shown that the computational effort for solutions with the GMRes algorithm depends on the dimension of the linear system, the total energy of the scattering problem, and the accuracy criterion. In several cases with dimensions in the range 1110-5632, the GMRes algorithm outperformed the LAPACK direct solver, with speedups for the linear equation solution as large as a factor of 23.

Chatfield, David C.

State-selected chemical reaction dynamics at the S matrix level - Final-state specificities of near-threshold processes at low and high energies

State-to-state reaction probabilities are found to be highly final-state specific at state-selected threshold energies for the reactions O + H2 yield OH + H and H + H2 yield H2 + H. The study includes initial rotational states with quantum numbers 0-15, and the specificity is especially dramatic for the more highly rotationally excited reactants. The analysis is based on accurate quantum mechanical reactive scattering calculations. Final-state specificity is shown in general to increase with the rotational quantum number of the reactant diatom, and the trends are confirmed for both zero and nonzero values of the total angular momentum.

Chatfield, David C.

A new diabatic representation of the coupled potential energy surfaces for Na(3p P-2) + H2 yields Na(3s S-2) + H2 or NaH + H

A diabatic representation is presented of the coupled potential-energy surfaces for Na(3p P-2) + H2 yields Na (3s S-2) + H2 or NaH + H. The representation is designed to yield, upon diagonalization, realistic values for the two lowest energy adiabatic states at both asymptotes of the chemical reaction as well as near the conical intersection in the three-body interaction region. It is economical to evaluate and portable. It is suitable for dynamics calculations on both the quenching process and the electronically nonadiabatic chemical reaction.

Halvick, Philippe

Improved techniques for outgoing wave variational principle calculations of converged state-to-state transition probabilities for chemical reactions

Improved techniques and well-optimized basis sets are presented for application of the outgoing wave variational principle to calculate converged quantum mechanical reaction probabilities. They are illustrated with calculations for the reactions D + H2 yields HD + H with total angular momentum J = 3 and F + H2 yields HF + H with J = 0 and 3. The optimization involves the choice of distortion potential, the grid for calculating half-integrated Green's functions, the placement, width, and number of primitive distributed Gaussians, and the computationally most efficient partition between dynamically adapted and primitive basis functions. Benchmark calculations with 224-1064 channels are presented.

Mielke, Steven L.

Use of scaled external correlation, a double many-body expansion, and variational transition state theory to calibrate a potential energy surface for FH2

New ab initio results and a double many-body expansion formalism have been used to parameterize a new FH2 potential energy surface with improved properties near the saddle point and in the region of long-range attraction. The functional form of the new surface includes dispersion forces by a double many-body expansion. Stationary point properties for the new surface are calculated along with the product-valley barrier maxima of vibrationally adiabatic potential curves for F + H2 - HF(nu-prime = 3) + H, F + HD - HF(nu-prime = 3) + D, and F + D2 - DF(nu-prime = 4) + D. The new surface should prove useful for studying the effect on dynamics of a low, early barrier with a wide, flat bend potential.

Lynch, Gillian C.

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.

Benchmark calculations of thermal reaction rates. I - Quantal scattering theory

The thermal rate coefficient for the prototype reaction H + H2 yields H2 + H with zero total angular momentum is calculated by summing, averaging, and numerically integrating state-to-state reaction probabilities calculated by time-independent quantum-mechanical scattering theory. The results are very carefully converged with respect to all numerical parameters in order to provide high-precision benchmark results for confirming the accuracy of new methods and testing their efficiency.

Chatfield, David C.

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

Localized basis functions and other computational improvements in variational nonorthogonal basis function methods for quantum mechanical scattering problems involving chemical reactions

The Generalized Newton Variational Principle for 3D quantum mechanical reactive scattering is briefly reviewed. Then three techniques are described which improve the efficiency of the computations. First, the fact that the Hamiltonian is Hermitian is used to reduce the number of integrals computed, and then the properties of localized basis functions are exploited in order to eliminate redundant work in the integral evaluation. A new type of localized basis function with desirable properties is suggested. It is shown how partitioned matrices can be used with localized basis functions to reduce the amount of work required to handle the complex boundary conditions. The new techniques do not introduce any approximations into the calculations, so they may be used to obtain converged solutions of the Schroedinger equation.

Schwenke, David W.

Spectroscopic analysis of transition state energy levels - Bending-rotational spectrum and lifetime analysis of H3 quasibound states

Converged quantum mechanical calculations of scattering matrices and transition probabilities are reported for the reaction of H with H2 with total angular momentum 0, 1, and 4 as functions of total energy in the range 0.85-1.15 eV on an accurate potential energy surface. The resonance structure is illustrated with Argand diagrams. State-to-state reactive collision delay times and lifetimes are presented. For J = 0, 1, and 4, the lowest-energy H3 resonance is at total energies of 0.983, 0.985, and 1.01 eV, respectively, with lifetimes of about 16-17 fs. For J = 1 and 4 there is a higher-energy resonance at 1.10-1.11 eV. For J = 1 the lifetime is about 4 fs and for J = 4 it is about 1 fs.

Zhao, Meishan

Computational strategies and improvements in the linear algebraic variational approach to rearrangement scattering

The computational steps in calculating quantum mechanical reactive scattering amplitudes by the L2 generalized Newton variational principle are discussed with emphasis on computational strategies and recent improvements that make the calculations more efficient. Special emphasis is placed on quadrature techniques, storage management strategies, use of symmetry, and boundary conditions. It is concluded that an efficient implementation of these procedures provides a powerful algorithm for the accurate solution of the Schroedinger equation for rearrangements.

Schwenke, David W.

Effect of orbital and rotational angular momentum averaging on branching ratios of dynamical resonances in the reaction H + p-H2 yields o-H2 + H

The paper reports extensive quantum mechanical calculations of the product vibrational branching ratios in the reaction H + p-H2 yields o-H2 + H. The calculations involve total angular momentum up to 2 and excited as well as ground initial rotational states, and they are completely converged with up to 513 channels in individual total angular momentum/parity blocks. Comparisons are made with recent experiments by Nieh and Valentini.

Mladenovic, Mirjana

Supercomputer chemistry at the University of Minnesota

The focus of cluster chemistry studies at the University of Minnesota is on the adsorption of small hydrocarbon molecules on metal and graphite surfaces, and on the substitution defects in diamond and silicon. Of particular interest is the study of very large systems such as the C60 molecule and the C150H30 molecule. Many of the properties of such clusters converge extremely slowly toward bulk values. Ultrathin films and microporous media are discussed as well as the Brownian dynamics of polymers.

Almlof, Jan

Quantum mechanical algebraic variational methods for inelastic and reactive molecular collisions

The quantum mechanical problem of reactive or nonreactive scattering of atoms and molecules is formulated in terms of square-integrable basis sets with variational expressions for the reactance matrix. Several formulations involving expansions of the wave function (the Schwinger variational principle) or amplitude density (a generalization of the Newton variational principle), single-channel or multichannel distortion potentials, and primitive or contracted basis functions are presented and tested. The test results, for inelastic and reactive atom-diatom collisions, suggest that the methods may be useful for a variety of collision calculations and may allow the accurate quantal treatment of systems for which other available methods would be prohibitively expensive.

Schwenke, David W.