Abstract:A new time-dependent approach to the cumulative reaction probability, N(E), has been developed based on the famous formulation given by Miller and co-workers [J. Chem. Phys. 79, 4889 (1983)], N(E)=[(2π)2/2] tr[δ(E−H)Fδ(E−H)F]. Taking advantage of the fact that the flux operator has only two nonzero eigenvalues, we evaluate the trace efficiently in a direct product basis of the first flux operator eigenstates and the Hamiltonian eigenstates on the dividing surface (internal states). Because the microcanonical density operator, δ(E−H), will eliminate contributions to N(E) from an internal state with the energy much higher than the total energy E, we can minimize the number of internal states required by choosing a dividing surface with the lowest density of internal states. If the dividing surface is located in an asymptotic region, one just needs to include all the open channels, i.e., with internal energy lower than the total energy. Utilizing the Fourier transform for δ(E−H), we can obtain the information for all the energies desired by propagating these wave packets once. Thus the present approach will be much more efficient than the initial state selected wave packet (ISSWP) approach to N(E) for systems with many rotation degrees of freedom because the density of states in asymptotic region for such systems is much higher than that in the transition state region. With the present method one can also calculate the cumulative reaction probability from an initial state (or to a final state) by locating the second flux operator in the corresponding asymptotic region. This provides an alternative to the ISSWP approach which may be more efficient if the reaction probabilities from a large number of initial states are desired. The method is applied to the 3D H + H2 (even rotation) reaction for J=0 by locating the first dividing surface in the transition state region. The demonstration also shows an aspect less than ideal; the contribution to N(E) from a wave packet may be slightly larger than 1 or slightly smaller than 0, making it improper to interpret the contribution as a probability.

Direct Calculation of Cumulative Reaction Probabilities from Chebyshev Correlation Functions

Calculation of Multiple Initial State Selected Reaction Probabilities from Chebyshev Flux-Flux Correlation Functions: Influence of Reactant Internal Excitations on H + H2O → OH + H2.

Cumulative Reaction Probability Via Transition State Wave Packets

Quantum Wavepacket Method for State-to-state Reactive Cross Sections in Hyperspherical Coordinates

State-to-state Dynamics of the C+(2P) + SH(X2Π) → H(2S) + CS+(X2Σ+) Reaction Using a Time-Dependent Wave Packet and Quasi-Classical Trajectory Methods

The Cumulative Reaction Probability for the H-2+Oh Reaction

Accurate Potential Energy Surface and Quantum Reaction Rate Calculations for the H+Ch4 -> H-2+Ch3 Reaction

Accurate Quantum Calculation for the Benchmark Reaction H2+OH→H2O +H in Five-Dimensional Space: Reaction Probabilities for J=0

Accuracy of the Centrifugal Sudden Approximation in the H + CHD3 → H2 + CD3 Reaction

Comparison of Second-Order Split Operator and Chebyshev Propagator in Wave Packet Based State-to-state Reactive Scattering Calculations

The Three-Dimensional Nonadiabatic Dynamics Calculation of Dh2+ and Hd2+ Systems by Using the Trajectory Surface Hopping Method Based on the Zhu-Nakamura Theory

Quantum Calculations of Reaction Probabilities for HO + CO→ H + CO2 and Bound States of HOCO

Accurate Quantum Calculation for the Benchmark Reaction H-2+Oh-]H2o+H in 5-Dimensional Space - Reaction Probabilities for J=O

Calculation Of State-To-State Cross Sections For Triatomic Reaction By The Multi-Configuration Time-Dependent Hartree Method

The quantum transition state wavepacket method

A Reactant-Coordinate-based Time-Dependent Wave Packet Method for Triatomic State-to-state Reaction Dynamics: Application to the H + O2 Reaction.

Non-Born-Oppenheimer Dynamics Calculations Using the Coherent Switching with Decay of Mixing Method

Quantum Dynamics Calculation of Reaction Probability for H + Cl2 → HCl + Cl

An Interaction-Asymptotic Region Decomposition Method for General State-to-state Reactive Scatterings.

A Full-Dimensional Time-Dependent Wave Packet Study of the OH + CO → H + CO2 Reaction

Interaction-Asymptotic Region Decomposition Method for an Insertion Reaction: Application to the S(1D) + H2 Reaction.