Search NASA⌕ Search

Engineering topics

Neuhauser, Daniel

Publications and source records attributed to Neuhauser, Daniel.

Sparse-Stochastic Fragmented Exchange for Large-Scale Hybrid Time-Dependent Density Functional Theory Calculations

Here we extend our recently developed sparse-stochastic fragmented exchange formalism for ground-state near-gap hybrid DFT to calculate absorption spectra within linear-response time-dependent generalized Kohn-Sham DFT (LR-GKS-TDDFT) for systems consisting of thousands of valence electrons within a grid-based/plane-wave representation. A mixed deterministic/fragmented-stochastic compression of the exchange kernel, here using long-range explicit exchange functionals, provides an efficient method for accurate optical spectra. Both real-time propagation as well as frequency-resolved Casida-equation-type approaches for spectra are presented, and the method is applied to large molecular dyes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Optimized attenuated interaction: Enabling stochastic Bethe–Salpeter spectra for large systems

We develop an improved stochastic formalism for the Bethe–Salpeter equation (BSE), based on an exact separation of the effective-interaction W into two parts, W = (W – vW) + vW, where the latter is formally any translationally invariant interaction, vW(r – r'). When optimizing the fit of the exchange kernel vW to W, using a stochastic sampling W, the difference W – vW becomes quite small. Then, in the main BSE routine, this small difference is stochastically sampled. Furthermore, the number of stochastic samples needed for an accurate spectrum is then largely independent of system size. While the method is formally cubic in scaling, the scaling prefactor is small due to the constant number of stochastic orbitals needed for sampling W.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

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↗

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↗

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↗

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↗

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↗

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

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↗