Collisional stabilization of van der Waals states of ozone

Mixed quantum/classical theory (MQCT) approach to the dynamics of molecule-molecule collisions in complex systems

We developed a general theoretical approach and a user-ready computer code that permit study of the dynamics of collisional energy transfer and ro-vibrational energy exchange in complex molecule-molecule collisions. The method is a mixture of classical and quantum mechanics. The internal ro-vibrational motion of collision partners is treated quantum mechanically using a time-dependent Schrödinger equation that captures many quantum phenomena including state quantization and zero-point energy, propensity and selection rules for state-to-state transitions, quantum symmetry and interference phenomena. A significant numerical speed up is obtained by describing the translational motion of collision partners classically, using the Ehrenfest mean-field trajectory approach. Within this framework a family of approximate methods for collision dynamics is developed. Several benchmark studies for diatomic and triatomic molecules, such as H2O and ND3 collided with He, H2 and D2, show that the results of MQCT are in good agreement with full-quantum calculations in a broad range of energies, especially at high collision energies where they become nearly identical to the full quantum results. Numerical efficiency of the method and massive parallelism of the MQCT code permit us to embrace some of the most complicated collisional systems ever studied, such as C6H6 + He, CH3COOH + He and H2O + H2O. Application of MQCT to the collisions of chiral molecules such as CH3CHCH2O + He, and to molecule-surface collisions is also possible and will be pursued in the future.

Four Isotope-Labeled Recombination Pathways of Ozone Formation

A theoretical approach is developed for the description of all possible recombination pathways in the ozone forming reaction, without neglecting any process a priori, and without decoupling the individual pathways one from another. These pathways become physically distinct when a rare isotope of oxygen is introduced, such as ¹⁸O, which represents a sensitive probe of the ozone forming reaction. Each isotopologue of O₃ contains two types of physically distinct entrance channels and two types of physically distinct product wells, creating four recombination pathways. Calculations are done for singly and doubly substituted isotopologues of ozone, eight rate coefficients total. Two pathways for the formation of asymmetric ozone isotopomer exhibit rather different rate coefficients, indicating large isotope effect driven by ΔZPE-difference. Rate coefficient for the formation of symmetric isotopomer of ozone (third pathway) is found to be in between of those two, while the rate of insertion pathway is smaller by two orders of magnitude. These trends are in good agreement with experiments, for both singly and doubly substituted ozone. The total formation rates for asymmetric isotopomers are found to be somewhat larger than those for symmetric isotopomers, but not as much as in the experiment. Overall, the distribution of lifetimes is found to be very similar for the metastable states in symmetric and asymmetric ozone isotopomers.

Influence of the Coriolis effect on the properties of scattering resonances in symmetric and asymmetric isotopomers of ozone

Scattering resonances above dissociation threshold are computed for four isotopically substituted ozone species: 16O18O16O, 16O16O18O, 18O16O18O and 16O18O18O, using a variational method with accurate treatment of the rotation-vibration coupling terms (Coriolis effect) for all values of the total angular momentum J from 0 to 4. To make these calculations numerically affordable, a new approach was developed which employs one vibrational basis set optimized for a typical rotational excitation (J,Λ), to run coupled rotation-vibration calculations at several desired values of J. In order to quantify the effect of Coriolis coupling, new data are contrasted with those computed using the symmetric-top rotor approximation, where the rotation-vibration coupling terms are neglected. It is found that, overall, the major properties of scattering resonances (such as their lifetimes, the number of these states, and their cumulative partition function Q) are all influenced by the Coriolis effect and this influence grows as the angular momentum J is raised. However, it is found that the four isotopically substituted ozone molecules are affected roughly equally by the Coriolis coupling. When the ratio η of partition functions for asymmetric over symmetric ozone molecules is computed, the Coriolis effect largely cancels, and this cancelation seems to occur for all values of J. Therefore, it does not seem grounded to attribute any appreciable mass-independent symmetry-driven isotopic fractionation to the Coriolis coupling effect.

Development of a potential energy surface for the O3–Ar system: rovibrational states of the complex

Collisional stabilization is an important step in the process of atmospheric formation of ozone.

Several Levels of Theory for Description of Isotope Effects in Ozone: Symmetry Effect and Mass Effect

The essential components of theory for the description of isotope effects in recombination reaction that forms ozone are presented, including the introduction of three reaction pathways for symmetric and asymmetric isotopomers, a brief review of relevant experimental data for singly- and doubly substituted isotopologues, the definitions of ζ-effect and η-effect, and the introduction of isotopic enrichment δ. Two levels of theory are developed to elucidate the role of molecular symmetry, atomic masses, vibrational zero-point energies, and rotational excitations in the recombination process. The issue of symmetry is not trivial, since the important factors, such as 1/2 and 2, appear in seven different places in the formalism. It is demonstrated that if all these effects are taken into account properly, then no anomalous isotope effects emerge. At the next level of theory, a model is considered in which one scattering resonance (sitting right at the top of centrifugal barrier) is introduced per ro-vibrational channel. It is found that this approach is equivalent to statistical treatment with partition functions at the transition state. Accurate calculations using hyper-spherical coordinates show that no isotope effects come from difference in the number of states. In contrast, differences in vibrational and rotational energies lead to significant isotope effects. However, those effects appear to be local, found for the rather extreme values of rotational quantum numbers. They largely cancel when rate coefficients are computed for the thermal distribution of rotational excitations. Although large isotope effects (observed in experiments) are not reproduced here, this level of theory can be used as a foundation for more detailed computational treatment, with accurate information about resonance energies and lifetimes computed and included.

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.

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.

Three sources of errors in the Ehrenfest treatment of inelastic scattering and possible ways of resolving them

In order to identify the origin of possible errors in the mixed quantum/classical approach to inelastic scattering [A. Semenov and D. Babikov, J. Chem. Phys. 140, 044306 (2014) and A. Semenov, M.-L. Dubernet, and D. Babikov, J. Chem. Phys. 141, 114304 (2014)], a simplified model is considered that consists of one intermolecular degree of freedom and two intramolecular states, coupled by a simple potential. For this system, analytic derivations are carried out to determine (i) the exact quantum mechanical solution of the inelastic scattering problem, (ii) a simplified version of it with all oscillatory terms neglected, and (iii) the Ehrenfest solution in which the translational motion is described by the mean-field trajectory while the internal molecular motion is treated by the time-dependent Schrodinger equation. It is shown that the appropriate choice of velocity for the mean-field trajectory permits to enforce microscopic reversibility and gives results in excellent agreement with full-quantum results. The average velocity method of Billing is rigorously derived as a limiting case (of this more general approach), when reversibility is enforced locally, at the initial moment of time only. It is demonstrated that errors of state-to-state transition probabilities in the Ehrenfest approach occur at lower values of total energy E if the magnitudes of excitation energy ΔE, potential energy difference between the two states ΔV, and coupling of two states V₁₂ are large. Possible ways of applying this concept to rotational transitions in real molecules are explored, using examples from CO + CO inelastic scattering.

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.

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.

Rotational quenching of H2O by He: mixed quantum/classical theory and comparison with quantum results

The mixed quantum/classical theory (MQCT) formulated in the space-fixed reference frame is used to compute quenching cross sections of several rotationally excited states of water molecule by impact of He atom in a broad range of collision energies, and is tested against the full-quantum calculations on the same potential energy surface. In current implementation of MQCT method, there are two major sources of errors: one affects results at energies below 10 cm(-1), while the other shows up at energies above 500 cm(-1). Namely, when the collision energy E is below the state-to-state transition energy ΔE the MQCT method becomes less accurate due to its intrinsic classical approximation, although employment of the average-velocity principle (scaling of collision energy in order to satisfy microscopic reversibility) helps dramatically. At higher energies, MQCT is expected to be accurate but in current implementation, in order to make calculations computationally affordable, we had to cut off the basis set size. This can be avoided by using a more efficient body-fixed formulation of MQCT. Overall, the errors of MQCT method are within 20% of the full-quantum results almost everywhere through four-orders-of-magnitude range of collision energies, except near resonances, where the errors are somewhat larger.

On molecular origin of mass-independent fractionation of oxygen isotopes in the ozone forming recombination reaction

Theoretical treatment of ozone forming reaction is developed within the framework of mixed quantum/classical dynamics. Formation and stabilization steps of the energy transfer mechanism are both studied, which allows simultaneous capture of the delta zero-point energy effect and η-effect and identification of the molecular level origin of mass-independent isotope fractionation. The central role belongs to scattering resonances; dependence of their lifetimes on rotational excitation, asymmetry; and connection of their vibrational wave functions to two different reaction channels. Calculations, performed within the dimensionally reduced model of ozone, are in semiquantitative agreement with experiment.