Born–Oppenheimer and Renner–Teller Quantum Dynamics of CH(X2Π) + D(2S) Reactions on Three CHD Potential Surfaces

Quantum Dynamics of Nonadiabatic Renner-Teller Effects in Atom + Diatom Collisions

We review the quantum nonadiabatic dynamics of atom + diatom collisions due to the Renner-Teller (RT) effect, i.e., to the Hamiltonian operators that contain the total spinless electronic angular momentum L̂. As is well-known, this rovibronic effect is large near collinear geometries when at least one of the interacting states is doubly degenerate. In general, this occurs in insertion reactions and at short-range, where the potential wells exhibit deep minima and support metastable complexes. Initial-state-resolved reaction probabilities, integral cross sections, and thermal rate constants are calculated via the real wavepacket method, solving the equation of motion with an approximated or with an exact spinless RT Hamiltonian. We present the dynamics of 10 single-channel or multichannel reactions showing how RT effects depend on the product channels and comparing with the Born-Oppenheimer (BO) approximation or coexisting conical-intersection (CI) interactions. RT effects not only can significantly modify the adiabatic dynamics or correct purely CI results, but also they can be very important in opening collision channels which are closed at the BO or CI level, as in electronic-quenching reactions. In the OH(A²Σ⁺) + Kr electronic quenching, where both nonadiabatic effects (CI and RT) coexist, they are in competition because CI dominates the reactivity but RT couplings reduce the large CI cross section and open a CI-forbidden evolution toward products, so that CI + RT quantum results are in good agreement with experimental or semiclassical findings. The different roles of these couplings are due to the unlike nuclear geometries where they are large: rather far from or near to linearity for CI or RT, respectively. The OH(A²Σ⁺) + Kr electronic quenching was investigated with the exact RT Hamiltonian, validating the approximated one, which was employed for all other collisions.

Non-adiabatic Quantum Dynamics of the Dissociative Charge Transfer He++H2 → He+H+H. Front Chem 2019;7:249. [PMID: 31041310 PMCID: PMC6477054 DOI: 10.3389/fchem.2019.00249]

We present the non-adiabatic, conical-intersection quantum dynamics of the title collision where reactants and products are in the ground electronic states. Initial-state-resolved reaction probabilities, total integral cross sections, and rate constants of two H₂ vibrational states, v₀ = 0 and 1, in the ground rotational state (j₀ = 0) are obtained at collision energies E_coll ≤ 3 eV. We employ the lowest two excited diabatic electronic states of HeH2+ and their electronic coupling, a coupled-channel time-dependent real wavepacket method, and a flux analysis. Both probabilities and cross sections present a few groups of resonances at low E_coll, whose amplitudes decrease with the energy, due to an ion-induced dipole interaction in the entrance channel. At higher E_coll, reaction probabilities and cross sections increase monotonically up to 3 eV, remaining however quite small. When H₂ is in the v₀ = 1 state, the reactivity increases by ~2 orders of magnitude at the lowest energies and by ~1 order at the highest ones. Initial-state resolved rate constants at room temperature are equal to 1.74 × 10⁻¹⁴ and to 1.98 × 10⁻¹² cm³s⁻¹ at v₀ = 0 and 1, respectively. Test calculations for H₂ at j₀ = 1 show that the probabilities can be enhanced by a factor of ~1/3, that is ortho-H₂ seems ~4 times more reactive than para-H₂.

Dynamics and resonances of the H(2S) + CH+(X1Σ+) reaction in the electronic ground state: a detailed quantum wavepacket study

Initial state-selected and energy resolved channel-specific reaction probabilities, integral cross sections and thermal rate constants of the H(²S) + CH⁺(X¹Σ⁺) reaction are calculated within the coupled states approximation by a time-dependent wave packet propagation method. The new ab initio global potential energy surface (PES) of the electronic ground state (1 ²A') of the system, recently reported by Li et al. [J. Chem. Phys., 2015, 142, 124302], is employed for this purpose. All partial wave contributions up to the total angular momentum J = 60 are considered to obtain the converged integral reaction cross section up to a collision energy of 1.0 eV. Thermal rate constants are calculated by averaging the reaction cross sections over the Boltzmann distribution of energies and compared with the available theoretical and experimental results for the temperature range 10-1000 K. Investigation of the channel-specific reaction attributes shows that the H abstraction (CH⁺ destruction) channel is highly favored over the H exchange channel. The effect of rotational and vibrational excitations of the CH⁺ reagent on the dynamics is also studied. The resonances formed during the course of the reaction are also identified by calculating the transition state spectrum and characterized in terms of the eigenfunctions and lifetimes. More than 260 vibrational levels are obtained and their eigenfunctions are calculated, which are represented in terms of the nodal assignments and the eigenenergies. They reveal both the local and hyperspherical behavior for the bound and quasibound states of the CH₂⁺ complex in the ground 1 ²A' surface. The lifetime analysis of the quasibound states indicates that the CH₂⁺ resonances survive for as long as ∼400 fs at high energies (E ∼ 2.0 eV) and are expected to decay faster with further increasing energy. Finally, the type of mechanism for the formation of the product (C⁺ + H₂) is elucidated.

Theoretical investigation of rotationally inelastic collisions of CH(X2Π) with hydrogen atoms

We report calculations of state-to-state cross sections for collision-induced rotational transitions of CH(X²Π) with atomic hydrogen. These calculations employed the four adiabatic potential energy surfaces correlating CH(X²Π) + H(²S), computed in this work through the multi-reference configuration interaction method [MRCISD + Q(Davidson)]. Because of the presence of deep wells on three of the potential energy surfaces, the scattering calculations were carried out using the quantum statistical method of Manolopoulos and co-workers [Chem. Phys. Lett. 343, 356 (2001)]. The computed cross sections included contributions from only direct scattering since the CH₂ collision complex is expected to decay predominantly to C + H₂. Rotationally energy transfer rate constants were computed for this system since these are required for astrophysical modeling.