Communication: Rigorous quantum dynamics of O + O2 exchange reactions on an ab initio potential energy surface substantiate the negative temperature dependence of rate coefficients

Lifetimes and decay mechanisms of isotopically substituted ozone above the dissociation threshold: matching quantum and classical dynamics

Energies and lifetimes of vibrational resonances were computed for 18O-enriched isotopologue 50O3 = {16O16O18O and 16O18O16O} of the ozone molecule using hyperspherical coordinates and the method of complex absorbing potential. Various types of scattering resonances were identified, including roaming OO-O rotational states, the series corresponding to continuation of bound vibrational resonances of highly excited bending or symmetric stretching vibrational modes. Such a series become metastable above the dissociation limit. The coupling between the vibrationally excited O2 fragment and rotational roaming gives rise to Feshbach type resonances in ozone. Different paths for the formation and decay of symmetric 16O18O16O and asymmetric species 16O16O18O were also identified. The symmetry properties of the total rovibronic wave functions of the 18O-enriched isotopologues are discussed in the context of allowed dissociation channels.

Ab initio study of the O3-N2 complex: Potential energy surface and rovibrational states

The formation and destruction of O₃ within the Chapman cycle occurs as a result of inelastic collisions with a third body. Since N₂ is the most abundant atmospheric molecule, it can be considered as the most typical candidate when modeling energy-transfer dynamics. We report a new ab initio potential energy surface (PES) of the O₃-N₂ van der Waals complex. The interaction energies were calculated using the explicitly correlated single- and double-excitation coupled cluster method with a perturbative treatment of triple excitations [CCSD(T)-F12a] with the augmented correlation-consistent triple-zeta aug-cc-pVTZ basis set. The five-dimensional PES was analytically represented by an expansion in spherical harmonics up to eighth order inclusive. Along with the global minimum of the complex (D_e = 348.88 cm^-1), with N₂ being perpendicular to the O₃ plane, six stable configurations were found with a smaller binding energy. This PES was employed to calculate the bound states of the O₃-N₂ complex with both ortho- and para-N₂ for total angular momentum J = 0 and 1, as well as dipole transition probabilities. The nature of the bound states of the O₃-oN₂ and O₃-pN₂ species is discussed based on their rovibrational wave functions.

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.

Rotationally inelastic scattering of O3–Ar: state-to-state rates with the multiconfigurational time dependent Hartree method

The rates of state-changing collisions are compared for different isotopologues of ozone from quantum scattering calculations with the MCTDH method.

The Role of Ozone Vibrational Resonances in the Isotope Exchange Reaction 16O16O + 18O → 18O16O + 16O: The Time-Dependent Picture

We consider the time-dependent dynamics of the isotope exchange reaction in collisions between an oxygen molecule and an oxygen atom: ¹⁶O¹⁶O + ¹⁸O → ¹⁶O¹⁸O + ¹⁶O. A theoretical approach using the multiconfiguration time-dependent Hartree method was employed to model the time evolution of the reaction. Two potential surfaces available in the literature were used in the calculations, and the results obtained with the two surfaces are compared with each other as well as with results of a previous theoretical time-independent approach. A good agreement for the reaction probabilities with the previous theoretical results is found. Comparing the results obtained using two potential energy surfaces allows us to understand the role of the reef/shoulder-like feature in the minimum energy path of the reaction in the isotope exchange process. Also, it was found that the distribution of final products of the reaction is highly anisotropic, which agrees with experimental observations and, at the same time, suggests that the family of approximated statistical approaches, assuming a randomized distribution over final exit channels, is not applicable to this case.

First-Principles Computed Rate Constant for the O + O2 Isotopic Exchange Reaction Now Matches Experiment

We show, by performing exact time-independent quantum molecular scattering calculations, that the quality of the ground electronic state global potential energy surface appears to be of utmost importance in accurately obtaining even as strongly averaged quantities as kinetic rate constants. The oxygen isotope exchange reaction, ¹⁸O + ³²O₂, motivated by the understanding of a complex long-standing problem of isotopic ozone anomalies in the stratosphere and laboratory experiments, is explored in this context. The thermal rate constant for this key reaction is now in quantitative agreement with all experimental data available to date. A significant recent progress at the frontier of three research domains, advanced electronic structure calculations, ultrasensitive spectroscopy, and quantum scattering calculations, has therefore permitted a breakthrough in the theoretical modeling of this crucial collision process from first principles.

Properties of Feshbach and “shape”-resonances in ozone and their role in recombination reactions and anomalous isotope effects

Three reaction pathways for formation of symmetric and asymmetric isotopologues of ozone.

On stabilization of scattering resonances in recombination reaction that forms ozone

Calculations of energy transfer in the recombination reaction that forms ozone are carried out within the framework of the mixed quantum/classical theory and using the dimensionally reduced 2D-model of ozone molecule, with bending motion neglected. Recombination rate coefficients are obtained at room temperature for symmetric and asymmetric isotopomers of singly and doubly substituted isotopologues. The processes of resonance formation, spontaneous decay, collisional dissociation, and stabilization by bath gas (Ar) are all characterized and taken into account within the steady-state approximation for kinetics. The focus is on stabilization step, where the mysterious isotopic η-effect was thought to originate from. Our results indicate no difference in cross sections for stabilization of scatteringresonances in symmetric and asymmetric isotopomers. As practical results, the general and simple analytic models for stabilization and dissociation cross sections are presented, which can be applied to resonances in any ozone molecule, symmetric or asymmetric, singly or doubly substituted. Present calculations show some isotope effect that looks similar to the experimentally observed η-effect, and the origin of this phenomenon is in the rates of formation/decay of scatteringresonances, determined by their widths, that are somewhat larger in asymmetric isotopomers than in their symmetric analogues. However, the approximate two-dimensional model used here is insufficient for consistent and reliable description of all features of the isotopic effect in ozone. Calculations using an accurate 3D model are still needed.

Dissociative chemisorption of methane on Ni(111) using a chemically accurate fifteen dimensional potential energy surface

A new chemically accurate potential energy surface for the dissociative chemisorption of methane on the rigid Ni(111) surface.

Testing the Palma-Clary Reduced Dimensionality Model Using Classical Mechanics on the CH4 + H → CH3 + H2 Reaction

Application of exact quantum scattering methods in theoretical reaction dynamics of bimolecular reactions is limited by the complexity of the equations of nuclear motion to be solved. Simplification is often achieved by reducing the number of degrees of freedom to be explicitly handled by freezing the less important spectator modes. The reaction cross sections obtained in reduced-dimensionality (RD) quantum scattering methods can be used in the calculation of rate coefficients, but their physical meaning is limited. The accurate test of the performance of a reduced-dimensionality method would be a comparison of the RD cross sections with those obtained in accurate full-dimensional (FD) calculations, which is not feasible because of the lack of complete full-dimensional results. However, classical mechanics allows one to perform reaction dynamics calculations using both the RD and the FD model. In this paper, an RD versus FD comparison is made for the 8-dimensional Palma-Clary model on the example of four isotopologs of the CH4 + H → CH3 + H2 reaction, which has 12 internal dimensions. In the Palma-Clary model, the only restriction is that the methyl group is confined to maintain C3v symmetry. Both RD and FD opacity and excitation functions as well as differential cross sections were calculated using the quasiclassical trajectory method. The initial reactant separation has been handled according to our one-period averaging method [ Nagy et al. J. Chem. Phys. 2016, 144, 014104 ]. The RD and FD excitation functions were found to be close to each other for some isotopologs, but in general, the RD reactivity parameters are lower than the FD reactivity parameters beyond statistical error, and for one of the isotopologs, the deviation is significant. This indicates that the goodness of RD cross sections cannot be taken for granted.

Rovibrational energy transfer and dissociation in O2-O collisions

A set of state-specific transition rates for each rovibrational level is generated for the O2(X(3)Σ(g)(-))-O(3)P system using the quasi-classical trajectory method at temperatures observed in hypersonic flows. A system of master equations describes the relaxation of the rovibrational ensemble to thermal equilibrium under ideal heat bath conditions at a constant translational temperature. Vibrational and rotational relaxation times, obtained from the average internal energies, exhibit a pattern inherent in a chemically reactive collisional pair. An intrinsic feature of the O3 molecular system with a large attractive potential is a weak temperature dependence of the rovibrational transition rates. For this reason, the quasi-steady vibrational and rotational temperatures experience a maximum at increasing translational temperature. The energy rate coefficients, that characterize the average loss of internal energy due to dissociation, quickly diminish at high temperatures, compared to other molecular systems.

A full-dimensional model of ozone forming reaction: the absolute value of the recombination rate coefficient, its pressure and temperature dependencies

Rigorous calculations of scattering resonances in ozone are carried out for a broad range of rotational excitations with a detailed analysis of their properties and contribution into recombination process.

Quantum Dynamics of the (18)O + (36)O2 Collision Process

We report full quantum cross sections and rate constants for the (18)O + (36)O2 → (36)O2 + (18)O collision process. This constitutes to the best of our knowledge the first dynamical study of the (18)O(18)O(18)O system, with three identical (18)O oxygen atoms. We emphasize the comparison with the (16)O + (32)O2 collision as this latter presents the exact same features as the one treated here, except the consistent change of mass for all three atoms. We find very similar behaviors in the cross sections, and we confirm that the rates are faster when three identical nuclei are involved. In particular, we cannot dynamically study this system with classical trajectory methods, and we have to include properly the indistinguishability of the three (18)O nuclei; however, we note some slight differences with the (16)O(16)O(16)O benchmark system, and we focus our analysis on their origin.