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

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.

High-Energy Quantum Dynamics of the 15N + o-14N14N Rovibrational Activation and Isotope Exchange Processes

We report full quantum reaction probabilities, computed within the framework of time-independent quantum mechanics using hyperspherical coordinates, for the 15N + 14N14N inelastic and reactive collision processes, restricted to total angular momentum J = 0, for kinetic energies up to 4.5 eV. We take advantage of the nonzero (i = 1) nuclear spin of 14N, leading to the existence of two nuclear spin isomers of 14N14N, namely, ortho- and para-14N14N, to restrict the study to the ortho molecular nitrogen species, with even rotational quantum number j = 0, 2, ... states. Specifically, we start with diatomic reagents ortho-14N14N in the initial rotational state j = 0. A comparison with similar works previously published by other groups using time-dependent wave packet and quasi-classical trajectory methods for the 14N + 14N14N fully symmetric collision is given. We find that reactive processes 15N + 14N14N involving atom exchange do not happen for collision energies less than 2.2 eV. Collisions at energies of around 2.0 eV are most effective for populating reactants' rovibrational states, that is, for inelastic scattering, whereas those at energies close to 5.0 eV yield a newly formed 14N15N isotopologue in a wide variety of excited vibrational levels.

High Resolution Infrared Spectroscopy in Support of Ozone Atmospheric Monitoring and Validation of the Potential Energy Function

The first part of this review is a brief reminder of general information concerning atmospheric ozone, particularly related to its formation, destruction, observations of its decrease in the stratosphere, and its increase in the troposphere as a result of anthropogenic actions and solutions. A few words are said about the abandonment of the Airbus project Alliance, which was expected to be the substitute of the supersonic Concorde. This project is over due to the theoretical evaluation of the impact of a fleet in the stratosphere and has been replaced by the A380, which is now operating. The largest part is devoted to calculations and observations of the transitions in the infrared range and their applications for the atmosphere based both on effective models (Hamiltonian, symmetry rules, and dipole moments) and ab initio calculations. The complementarities of the two approaches are clearly demonstrated, particularly for the creation of an exhaustive line list consisting of more than 300,000 lines reaching experimental accuracies (from 0.00004 to 0.001 cm-1) for positions and a sub percent for the intensities in the 10 microns region. This contributes to definitively resolving the issue of the observed discrepancies between line intensity data in different spectral regions: between the infrared and ultraviolet ranges, on the one hand, and between 10 and 5 microns on the other hand. The following section is devoted to the application of recent work to improve the knowledge about the behavior of potential function at high energies. A controversial issue related to the shape of the potential function in the transition state range near the dissociation is discussed.

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.

Direct time delay computation applied to the O + O2 exchange reaction at low energy: Lifetime spectrum of O3 * species

We report full quantum dynamical calculations for lifetimes of scattering resonances, among which are true metastable states, of the intermediate heavy ozone complex ⁵⁰O₃ ^* of the ¹⁸O + ¹⁶O¹⁶O reaction, for any value of the total angular momentum quantum number J. We show that computations for nonzero values of J are mandatory in order to properly analyze resonances and time delays, with a view to establish a somewhat comprehensive eigenlife spectrum of the complex O₃ ^*. Calculations have been performed in a given low to moderate energy range, including the interval between zero-point energies (ZPEs) of reagents and product species. Quasi-bound states tend to be more numerous, and eigenlifetimes themselves are seen to increase with J, reaching unusually large values for J = 30. A very dense forest of O₃ ^* species is pictured already for J greater than 20, especially at the highest energies considered, leading to a quasi-continuum of metastable states. On the contrary, they appear as rather sparse and isolated at J = 0 and lower energies, including the domain between ¹⁸O¹⁶O and ¹⁶O¹⁶O ZPEs, embedded among many overlapping resonances that turn out to be not long-lived enough to be associated with genuine metastable states.

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.

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.