Several levels of theory for description of isotope effects in ozone: Effect of resonance lifetimes and channel couplings

In this paper, two levels of theory are developed to determine the role of scattering resonances in the process of ozone formation. At the lower theory level, we compute resonance lifetimes in the simplest possible way, by neglecting all couplings between the diabatic vibrational channels in the problem. This permits to determine the effect of "shape" resonances, trapped behind the centrifugal barrier and populated by quantum tunneling. At the next level of theory, we include couplings between the vibrational channels, which permits to determine the role of Feshbach resonances and interaction of different reaction pathways on the global PES of ozone. Pure shape resonances are found to contribute little to the overall recombination process since they occur rather infrequently in the spectrum, in the vicinity of the top of the centrifugal barrier only. Moreover, the associated isotope effects are found to disagree with experimental data. By contrast, Feshbach-type resonances are found to make dominant contribution to the process. They occur in a broader range of spectrum, and their density of states is much higher. The properties of Feshbach resonances are studied in detail. They explain the isotopic ζ -effect, giving theoretical prediction in good agreement with experiments for both singly and doubly substituted ozone molecules. Importantly, Feshbach resonances also contribute to the isotopic η -effect, moving theoretical predictions in the right direction. Some differences with experimental data remain, which indicates that there may be another additional source of the η -effect.

Quantum, Statistical, and Quasiclassical Trajectory Studies For the Ne + HeH(+) → NeH(+) + He Reaction on the Ground Electronic State

Real wave packet, statistical quantum, and quasiclassical trajectory methods were employed to study the dynamics of Ne + HeH(+)(v0,j0) → He + NeH(+) reaction on an ab initio potential energy surface [J. Phys. Chem. A 2013, 117, 13070-13078]. Quantum and statistical quantum calculations were performed within the centrifugal sudden (CS) approximation as well as including the Coriolis coupling (CC). Dense oscillatory structures of the quantum reaction probabilities and fair agreement between quantum and statistical cross sections suggest a complex forming mechanism for the reaction. No significant differences between cross sections obtained within the CS and CC approaches are observed. Quasiclassical trajectory results give an excellent average description of the quantum CC results. At low collision energies, there is a substantial decrease in reactivity for the reaction upon rovibrational excitation. Initial state selected rate constants for the title reaction are calculated between 20 and 1000 K, and the calculated value at 300 K agrees quite well with the available experimental result. Reaction cross sections and rate constants are also compared with those calculated via the Langevin capture model for exothermic reactions.

Wave packet and statistical quantum calculations for the He + NeH⁺ → HeH⁺ + Ne reaction on the ground electronic state

A real wave packet based time-dependent method and a statistical quantum method have been used to study the He + NeH(+) (v, j) reaction with the reactant in various ro-vibrational states, on a recently calculated ab initio ground state potential energy surface. Both the wave packet and statistical quantum calculations were carried out within the centrifugal sudden approximation as well as using the exact Hamiltonian. Quantum reaction probabilities exhibit dense oscillatory pattern for smaller total angular momentum values, which is a signature of resonances in a complex forming mechanism for the title reaction. Significant differences, found between exact and approximate quantum reaction cross sections, highlight the importance of inclusion of Coriolis coupling in the calculations. Statistical results are in fairly good agreement with the exact quantum results, for ground ro-vibrational states of the reactant. Vibrational excitation greatly enhances the reaction cross sections, whereas rotational excitation has relatively small effect on the reaction. The nature of the reaction cross section curves is dependent on the initial vibrational state of the reactant and is typical of a late barrier type potential energy profile.

Laser-driven isomerization of HCN → HNC: the importance of rotational excitation

We report a time-dependent quantum wave packet theory, which is employed to interpret the isomerization dynamics of HCN molecules induced by an intense picosecond infrared laser field. Considering the molecular rotational degrees of the freedom, the wave functions are expanded in terms of molecular rotational bases. Our full-dimensional quantum model includes the full Coriolis coupling in the molecular kinetic energy Hamiltonian and dipole approximation in interaction terms. The numerical results show that the field-induced molecule rotational excitation plays an important role in the isomerization dynamical process. Some phenomena appear such as two-step two-photon absorption and highly oscillatory structure in rotational state distributions. The centrifugal sudden (CS) approximation calculation is also carried out and compared in this work; it is shown that the Coriolis couplings may lead to a significant decrease in the isomerization rate but highly enhanced molecular rotational excitation.

Quantum dynamical study of the He + NeH+ reaction on a new analytical potential energy surface

An analytical potential energy surface (PES) for the ground state of the [HeHNe](+) system has been constructed from a set of 19,605 ab initio data points, obtained from coupled cluster singles and doubles with perturbative triples correction calculations and the aug-cc-pVQZ basis set. The PES is based on the many-body expansion form proposed by Aguado and Paniagua (J. Chem. Phys. 1992, 96, 1265), and it has a root-mean-square error of 0.03 kcal/mol. The minimum energy pathways (MEPs) for different Ne-H-He angles are calculated, and it is found that the MEP for 180° (linear) goes through the deepest potential energy well. Preliminary quantum dynamical studies are performed for the He + NeH(+) (v = 0-2, j = 0-3) → HeH(+) + Ne reaction in the 0.0-0.5 eV collision energy range. Quantum calculations are carried out using a time-dependent wave packet method within the centrifugal sudden approximation. Reaction probabilities exhibit strong oscillatory behavior arising because of the metastable [HeHNe](+). Vibrational excitation has been found to enhance the reaction cross sections.

Time-dependent quantum wave packet study of the Ar+H2+→ArH(+)+H reaction on a new ab initio potential energy surface for the ground electronic state (1(2)A')

A new global potential energy surface for the ground electronic state (1(2)A') of the Ar+H2(+)→ArH(+)+H reaction has been constructed by multi-reference configuration interaction method with Davidson correction and a basis set of aug-cc-pVQZ. Using 6080 ab initio single-point energies of all the regions for the dynamics, a many-body expansion function form has been used to fit these points. The quantum reactive scattering dynamics calculations taking into account the Coriolis coupling (CC) were carried out on the new potential energy surface over a range of collision energies (0.03-1.0 eV). The reaction probabilities and integral cross sections for the title reaction were calculated. The significance of including the CC quantum scattering calculation has been revealed by the comparison between the CC and the centrifugal sudden approximation calculation. The calculated cross section is in agreement with the experimental result at collision energy 1.0 eV.

Accurate study on the quantum dynamics of the He + HeH(+) (X1Σ+) reaction on a new ab initio potential energy surface for the lowest 1(1)A' electronic singlet state

A time-dependent quantum wave packet method is used to investigate the dynamics of the He + HeH(+)(X(1)Σ(+)) reaction based on a new potential energy surface [Liang et al., J. Chem. Phys.2012, 136, 094307]. The coupled channel (CC) and centrifugal-sudden (CS) reaction probabilities as well as the total integral cross sections are calculated. A comparison of the results with and without Coriolis coupling revealed that the number of K states N(K) (K is the projection of the total angular momentum J on the body-fixed z axis) significantly influences the reaction threshold. The effective potential energy profiles of each N(K) for the He + HeH(+) reaction in a collinear geometry indicate that the barrier height gradually decreased with increased N(K). The calculated time evolution of CC and CS probability density distribution over the collision energy of 0.27-0.36 eV at total angular momentum J = 50 clearly suggests a lower reaction threshold of CC probabilities. The CC cross sections are larger than the CS results within the entire energy range, demonstrating that the Coriolis coupling effect can effectively promote the He + HeH(+) reaction.