Localized and delocalized bound states of the main isotopologue 48O3 and of 18O-enriched 50O3 isotopomers of the ozone molecule near the dissociation threshold

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.

Ozone spectroscopy in the terahertz range from first high-resolution Synchrotron SOLEIL experiments combined with far-infrared measurements and ab initio intensity calculations

Ozone is one of the important molecules in terms of the impact on the atmospheric chemistry, climate changes, bio- and eco-systems and human health. It has a strong absorption in the microwave, terahertz and far-infrared spectral ranges where a large part of the Earth's outgoing longwave radiation to space is located. In this work, the observations, and analyses of the ozone high-resolution spectra in the THz range recorded using the Synchrotron light source of the SOLEIL CNRS equipment are reported for the first time. Thanks to the exceptional brightness of the Synchrotron radiation and to the signal/noise ratio, it was possible to observe many more ozone transitions of the cold rotational band and the hot ν2-ν2 band in the range 0.9-6 THz compared to the previous works. In addition, we have carried out new measurements and assignments for the ν2 band. The simultaneous fit of the rotational band GS-GS, the hot band ν2-ν2 and the FIR ν2 band yielded an overall weighted standard deviation of 0.68 for 13,466 line positions within the experimental accuracy. This includes all previously available MW (with the best uncertainty 0.1 - 10 kHz), FIR data and the original SOLEIL measurements that provided experimental accuracy of 0.00005 - 0.0001 cm-1 for the best lines. Significant deviations in new experimental spectra compared to available spectroscopic databases were evidenced, particularly for the line positions and energy levels at high J, Ka rotational quantum numbers that are the most pronounced in the 4.5 - 6 THz range. Accurate ab initio calculations of line intensities combined with empirically fitted line positions were used to create new linelists that permit theoretical modelling of the transmittance in a good agreement with the Synchrotron spectra in the entire range of observations for various pressures and optical paths. The region near 100 cm-1 and above appears to be more sensitive to the temperature conditions that should be considered in atmospheric observation for the currently operational and future ground based and space missions.

Ozone Formation in Ternary Collisions: Theory and Experiment Reconciled

The present Letter shows that the formation of ozone in ternary collisions O+O_{2}+M-the primary mechanism of ozone formation in the stratosphere-at temperatures below 200 K (for M=Ar) proceeds through a formation of a temporary complex MO_{2}, while at temperatures above ∼700  K, the reaction proceeds mainly through a formation of long-lived vibrational resonances of O_{3}^{*}. At intermediate temperatures 200-700 K, the process cannot be viewed as a two-step mechanism, often used to simplify and approximate collisions of three atoms or molecules. The developed theoretical approach is applied to the reaction O+O_{2}+Ar because of extensive experimental data available. The rate coefficients for the formation of O_{3} in ternary collisions O+O_{2}+Ar without using two-step approximations were computed for the first time as a function of collision energy. Thermally averaged coefficients were derived for temperatures 5-900 K. It is found that the majority of O_{3} molecules formed initially are weakly bound. Accounting for the process of vibrational quenching of the nascent population, a good agreement with available experimental data for temperatures 100-900 K is obtained.

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.

Efficient Method for an Approximate Treatment of the Coriolis Effect in Calculations of Quantum Dynamics and Spectroscopy, with Application to Scattering Resonances in Ozone

A numerical approach is developed to capture the effect of rotation-vibration coupling in a practically affordable way. In this approach only a limited number of adjacent rotational components are considered to be coupled, while the couplings to other rotational components are neglected. This partially coupled (PC) approach permits to reduce the size of Hamiltonian matrix significantly, which enables the calculations of ro-vibrational states above dissociation threshold (scattering resonances) for large values of total angular momentum. This method is employed here to reveal the role of the Coriolis effect in the ozone formation reaction at room temperature, dominated by large values of total angular momentum states, on the order of J = 24 and 28. We found that, overall, the effect of ro-vibrational coupling is not minor for large J. Compared to the results of symmetric top rotor approximation, where the ro-vibrational coupling is neglected, we found that the widths of scattering resonances, responsible for the lifetimes of metastable ozone states, remain nearly the same (on average), but the number of these states increases by as much as 20%. We also found that these changes are nearly the same in symmetric and asymmetric ozone isotopomers ¹⁶O¹⁸O¹⁶O and ¹⁶O¹⁶O¹⁸O. Therefore, based on the results of these calculations, the Coriolis coupling does not seem to favor the formation of asymmetric ozone molecules and thus cannot be responsible for symmetry-driven mass-independent fractionation of oxygen isotopes.

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.

Diagonal Born-Oppenheimer corrections to the ground electronic state potential energy surfaces of ozone: improvement of ab initio vibrational band centers for the 16O3, 17O3 and 18O3 isotopologues

Mass-dependent diagonal Born-Oppenheimer corrections (DBOCs) to the ab initio electronic ground state potential energy surface for the main ¹⁶O₃ isotopologue and for homogeneous isotopic substitutions ¹⁷O₃ and ¹⁸O₃ of the ozone molecule are reported for the first time. The system being of strongly multiconfigurational character, multireference configuration interaction wave function ansatz with different complete active spaces was used. The reliable DBOC calculations with the targeted accuracy were possible to carry out up to about half of the dissociation threshold D₀. The comparison with the experimental band centers shows a significant improvement of the accuracy with respect to the best Born-Oppenheimer (BO) ab initio calculations reducing the total root-mean-squares (calculated-observed) deviations by about a factor of two. For the set of ¹⁶O₃ vibrations up to five bending and four stretching quanta, the mean (calculated-observed) deviations drop down from 0.7 cm^-1 (BO) to about 0.1 cm^-1, with the most pronounced improvement seen for bending states and for mixed bending-stretching polyads. In the case of bending band centers directly observed under high spectral resolutions, the errors are reduced by more than an order of magnitude down to 0.02 cm^-1 from the observed levels, approaching nearly experimental accuracy. A similar improvement for heavy isotopologues shows that the reported DBOC corrections almost remove the systematic BO errors in vibrational levels below D₀/2, though the scatter increases towards higher energies. The possible reasons for this finding, as well as remaining issues are discussed in detail. The reported results provide an encouraging accuracy validation for the multireference methods of the ab initio theory. New sets of ab initio vibrational states can be used for improving effective spectroscopic models for analyses of the observed high-resolution spectra, particularly in the cases of accidental resonances with "dark" states requiring accurate theoretical predictions.

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.