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

Towards Understanding the Roaming Mechanism in H + MgH → Mg + HH Reaction

The roaming mechanism in the reaction H + MgH →Mg + HH is investigated by classical and quantum dynamics employing an accurate ab initio three-dimensional ground electronic state potential energy surface. The reaction dynamics are explored by running trajectories initialized on a four-dimensional dividing surface anchored on three-dimensional normally hyperbolic invariant manifold associated with a family of unstable orbiting periodic orbits in the entrance channel of the reaction (H + MgH). By locating periodic orbits localized in the HMgH well or involving H orbiting around the MgH diatom, and following their continuation with the total energy, regions in phase space where reactive or nonreactive trajectories may be trapped are found. In this way roaming reaction pathways are deduced in phase space. Patterns similar to periodic orbits projected into configuration space are found for the quantum bound and resonance eigenstates. Roaming is attributed to the capture of the trajectories in the neighborhood of certain periodic orbits. The complex forming trajectories in the HMgH well can either return to the radical channel or “roam” to the MgHH minimum from where the molecule may react

Dynamical Outcomes of Quenching: Reflections on a Conical Intersection

This review focuses on experimental studies of the dynamical outcomes following collisional quenching of electronically excited OH A² ∑ ⁺ radicals by molecular partners. The experimental observables include the branching between reactive and nonreactive decay channels, kinetic energy release and quantum state distributions of the products. Complementary theoretical investigations reveal regions of strong nonadiabatic coupling, known as conical intersections, which facilitate the quenching process. The dynamical outcomes observed experimentally are connected to the local forces and geometric properties of the nuclei in the conical intersection region. Dynamical calculations for the benchmark OH-H₂ system are in good accord with experimental observations, demonstrating that the outcomes reflect the strong coupling in the conical intersection region as the system evolves from the excited electronic state to quenched products

OH+ in astrophysical media: state-to-state formation rates, Einstein coefficients and inelastic collision rates with He

The rate constants required to model the OH$^+$ observations in different regions of the interstellar medium have been determined using state of the art quantum methods. First, state-to-state rate constants for the H$_2(v=0,J=0,1)$+ O$^+$($^4S$) $\rightarrow$ H + OH$^+(X ^3\Sigma^-, v', N)$ reaction have been obtained using a quantum wave packet method. The calculations have been compared with time-independent results to asses the accuracy of reaction probabilities at collision energies of about 1 meV. The good agreement between the simulations and the existing experimental cross sections in the $0.01-$1 eV energy range shows the quality of the results. The calculated state-to-state rate constants have been fitted to an analytical form. Second, the Einstein coefficients of OH$^+$ have been obtained for all astronomically significant ro-vibrational bands involving the $X^3\Sigma^-$ and/or $A^3\Pi$ electronic states. For this purpose the potential energy curves and electric dipole transition moments for seven electronic states of OH$^+$ are calculated with {\it ab initio} methods at the highest level and including spin-orbit terms, and the rovibrational levels have been calculated including the empirical spin-rotation and spin-spin terms. Third, the state-to-state rate constants for inelastic collisions between He and OH$^+(X ^3\Sigma^-)$ have been calculated using a time-independent close coupling method on a new potential energy surface. All these rates have been implemented in detailed chemical and radiative transfer models. Applications of these models to various astronomical sources show that inelastic collisions dominate the excitation of the rotational levels of OH$^+$. In the models considered the excitation resulting from the chemical formation of OH$^+$ increases the line fluxes by about 10 % or less depending on the density of the gas

Semiclassical Nonadiabatic Dynamics Based on Quantum Trajectories for the O(\u3csup\u3e3\u3c/sup\u3eP,\u3csup\u3e1\u3c/sup\u3eD)+H\u3csub\u3e2\u3c/sub\u3e System

The O(3P,1D)+H2→OH+H reaction is studied using trajectory dynamics within the approximate quantum potential approach. Calculations of the wave-packet reaction probabilities are performed for four coupled electronic states for total angular momentum J = 0 using a mixed coordinate/polar representation of the wave function. Semiclassical dynamics is based on a single set of trajectories evolving on an effective potential-energy surface and in the presence of the approximate quantum potential. Population functions associated with each trajectory are computed for each electronic state. The effective surface is a linear combination of the electronic states with the contributions of individual components defined by their time-dependent average populations. The wave-packet reaction probabilities are in good agreement with the quantum-mechanical results. Intersystem crossing is found to have negligible effect on reaction probabilities summed over final electronic states

Dynamics on Multiple Potential Energy Surfaces: Quantitative Studies of Elementary Processes Relevant to Hypersonics

The determination of thermal and vibrational relaxation rates of triatomic systems suitable for application in hypersonic model calculations is discussed. For this, potential energy surfaces for ground and electronically excited state species need to be computed and represented with high accuracy and quasiclassical or quantum nuclear dynamics simulations provide the basis for determining the relevant rates. These include thermal reaction rates, state-to-state cross-sections, or vibrational relaxation rates. For exemplary systems - [NNO], [NOO], and [CNO] - all individual steps are described and a literature overview for them is provided. Finally, as some of these quantities involve considerable computational expense, for the example of state-to-state cross sections the construction of an efficient model based on neural networks is discussed. All such data is required and being used in more coarse-grained computational fluid dynamics simulations.Comment: Review article, 46 pages, 8 figure

The hyperbolic matrix method for wave packet treatment of atomic and molecular dynamics

The hyperbolic matrix method for the treatment of atomic and molecular nuclear dynamics is derived by means of the wave packet technique within the Born–Oppenheimer approach formalism. This method allows one to calculate the evolution of a wave packet by means of the product of usual matrices instead of the time propagation of matrix exponentials. We provide a detailed description of this method, considering Tully’s model as an example. We also show good agreement with the Landau–Zener model application