Effect of internal excitations of reagent diatom on initial state-selected dynamics of C + OH reaction on its second excited (14A″) electronic state

Quantum and quasi-classical dynamics of the C(3P) + O2(3Σ-g) → CO(1Σ+) + O(1D) reaction on its electronic ground state

The dynamics of the C(³P) + O₂(³Σ⁻_g) → CO(¹Σ⁺) + O(¹D) reaction on its electronic ground state is investigated by using time-dependent wave packet propagation (TDWP) and quasi-classical trajectory (QCT) simulations. For the moderate collision energies considered (E_c = 0.001 to 0.4 eV, corresponding to a range from 10 K to 4600 K) the total reaction probabilities from the two different treatments of the nuclear dynamics agree very favourably. The undulations present in P(E) from the quantum mechanical treatment can be related to stabilization of the intermediate CO₂ complex with lifetimes on the 0.05 ps time scale. This is also confirmed from direct analysis of the TDWP simulations and QCT trajectories. Product diatom vibrational and rotational level resolved state-to-state reaction probabilities from TDWP and QCT simulations agree well except for the highest product vibrational states (v′ ≥ 15) and for the lowest product rotational states (j′ ≤ 10). Opening of the product vibrational level CO(v′ = 17) requires ∼0.2 eV from QCT and TDWP simulations with O₂(j = 0) and decreases to 0.04 eV if all initial rotational states are included in the QCT analysis, compared with E_c > 0.04 eV obtained from experiments. It is thus concluded that QCT simulations are suitable for investigating and realistically describe the C(³P) + O₂(³Σ⁻_g) → CO(¹Σ⁺) + O(¹D) reaction down to low collision energies when compared with results from a quantum mechanical treatment using TDWPs.

The dynamics of the C(³P) + O₂(³Σ⁻_g) → CO(¹Σ⁺) + O(¹D) reaction on its electronic ground state is investigated by using time-dependent wave packet propagation (TDWP) and quasi-classical trajectory (QCT) simulations.

Time-dependent quantum mechanical wave packet dynamics

Starting from a model study of the collinear (H, H₂) exchange reaction in 1959, the time-dependent quantum mechanical wave packet (TDQMWP) method has come a long way in dealing with systems as large as Cl + CH₄. The fast Fourier transform method for evaluating the second order spatial derivative of the wave function and split-operator method or Chebyshev polynomial expansion for determining the time evolution of the wave function for the system have made the approach highly accurate from a practical point of view. The TDQMWP methodology has been able to predict state-to-state differential and integral reaction cross sections accurately, in agreement with available experimental results for three dimensional (H, H₂) collisions, and identify reactive scattering resonances too. It has become a practical computational tool in predicting the observables for many A + BC exchange reactions in three dimensions and a number of larger systems. It is equally amenable to determining the bound and quasi-bound states for a variety of molecular systems. Just as it is able to deal with dissociative processes (without involving basis set expansion), it is able to deal with multi-mode nonadiabatic dynamics in multiple electronic states with equal ease. We present an overview of the method and its strength and limitations, citing examples largely from our own research groups.

Effect of Reagent Vibration and Rotation on the State-to-State Dynamics of the Hydrogen Exchange Reaction, H + H2 → H2 + H

State-to-state dynamics of the benchmark hydrogen exchange reaction H + H₂ (v = 0-4, j = 0-3) → H₂ (v', j') + H is investigated with the aid of the real wave packet approach of Gray and Balint-Kurti (J. Chem. Phys. 1998, 108, 950-962) and electronic ground BKMP2 potential energy surface of Boothroyd et al. (J. Chem. Phys. 1996, 104, 7139-7152). Initial state-selected and product state-resolved reaction probabilities, integral cross section, and product diatom vibrational and rotational level populations at a few collision energies are reported to elucidate the energy disposal mechanism. State-specific thermal rate constants are also calculated and compared with the available literature results. Coriolis coupling terms of the nuclear Hamiltonian are included, and calculations are parallelized over the helicity quantum number, Ω'. Attempts are made, in particular, to study the effect of reagent vibrational and rotational excitations on the dynamical attributes. It is found that the calculations become computationally expensive with reagent vibrational and rotational excitation. Reagent vibrational excitation is found to enhance the reactivity and has significant impact on the energy disposal to the vibrational and rotational degrees of freedom of the product. The interplay of reagent translational and vibrational energy on the product vibrational distribution unfolds an important aspect of the energy disposal mechanism. The effect of reagent rotation on the state-to-state dynamics is found not to be very significant, and the weak effect turns out to be specific to v'.

Quasi-Classical Trajectory Study of NH(3∑-) + NH(3∑-) Reactive Collisions

A full-dimension quasi-classical trajectory study of collisions between two NH radicals is presented. Interatomic interactions are represented by a previously reported global six-dimensional potential energy surface for singlet electronic state of the N₂H₂ system. This study suggests that the formation of N₂ from the collision of two NH radicals may occur via a one-step (NH + NH → N₂ + H + H) or two-step (NH + NH → N₂H + H → N₂ + H + H) microscopic reaction mechanism. A fast vibrational energy redistribution is observed in the four-body complex in the latter mechanism. Excitation functions are presented and discussed. A variant of the vibrational energy quantum mechanical threshold method was used to correct the zero-point energy leakage in the classical calculations. The influence of reactant's rotational and vibrational energy on reactivity was investigated by state-specific calculations with one or both reactants excited. Reaction rate constants for the ground and some rotationally excited states are presented using an Arrhenius-Kooij-like functional form.