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.

Semiclassical treatment of Feshbach resonances by transfer matrices

A semiclassical method is presented for the calculation of Feshbach resonance positions and widths. This approach, based on semiclassical transfer matrices, relies only on relatively short trajectory fragments, thus avoiding problems associated with the long trajectories needed in more straightforward semiclassical techniques. Complex resonance energies are obtained from an implicit equation that is developed to compensate for the inaccuracy of the stationary phase approximation underlying the semiclassical transfer matrix applications. Although this treatment requires calculation of transfer matrices for complex energies, an initial value representation method makes it possible to extract such quantities from ordinary real-valued classical trajectories. This treatment is applied to obtain positions and widths for resonances in a model two-dimensional system, and the results are compared to those obtained from accurate quantum mechanical calculations. The semiclassical method successfully captures the irregular energy dependence of resonance widths that vary over a range of more than two orders of magnitude. An explicit semiclassical expression for the width of narrow resonances is also presented and serves as a simpler, useful approximation for many cases.

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.

Oxygen and magnesium mass-independent isotopic fractionation induced by chemical reactions in plasma

Both the physical effect and the chemical conditions at the origin of the oxygen isotope variations in the solar system have been puzzling questions for 50 y. The data reported here bring the MIF effect (the mass-independent fractionation originally identified on ozone) back to the center of the debate. Similar to Ti isotopes, we observe that the MIF effect for O and Mg is triggered by redox reactions in plasmas. These observations reinforce the idea of a universal mechanism observable in photochemical reactions when molecular collisions involving indistinguishable isotopes yield a symmetrical complex stabilized as a chemical product.

Enrichment or depletion ranging from −40 to +100% in the major isotopes ¹⁶O and ²⁴Mg were observed experimentally in solids condensed from carbonaceous plasma composed of CO₂/MgCl₂/Pentanol or N₂O/Pentanol for O and MgCl₂/Pentanol for Mg. In NanoSims imaging, isotope effects appear as micrometer-size hotspots embedded in a carbonaceous matrix showing no isotope fractionation. For Mg, these hotspots are localized in carbonaceous grains, which show positive and negative isotopic effects so that the whole grain has a standard isotope composition. For O, no specific structure was observed at hotspot locations. These results suggest that MIF (mass-independent fractionation) effects can be induced by chemical reactions taking place in plasma. The close agreement between the slopes of the linear correlations observed between δ²⁵Mg versus δ²⁶Mg and between δ¹⁷O versus δ¹⁸O and the slopes calculated using the empirical MIF factor η discovered in ozone [M. H. Thiemens, J. E. Heidenreich, III. Science 219, 1073–1075; C. Janssen, J. Guenther, K. Mauersberger, D. Krankowsky. Phys. Chem. Chem. Phys. 3, 4718–4721] attests to the ubiquity of this process. Although the chemical reactants used in the present experiments cannot be directly transposed to the protosolar nebula, a similar MIF mechanism is proposed for oxygen isotopes: at high temperature, at the surface of grains, a mass-independent isotope exchange could have taken place between condensing oxides and oxygen atoms originated form the dissociation of CO or H₂O gas.

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.

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.

The role of rotation-vibration coupling in symmetric and asymmetric isotopomers of ozone

A theoretical framework and a computer code (SpectrumSDT) are developed for accurate calculations of coupled rotational-vibrational states in triatomic molecules using hyper-spherical coordinates and taking into account the Coriolis coupling effect. Concise final formulas are derived for the construction of the Hamiltonian matrix using an efficient combination of the variational basis representation and discrete variable representation methods with locally optimized basis sets and grids. First, the new code is tested by comparing its results with those of the APH3D program of Kendrick et al. [Kendrick, Pack, Walker, and Hayes, J. Chem. Phys. 110, 6673 (1999)]. Then, accurate calculations of the rovibrational spectra are carried out for doubly substituted symmetric (¹⁸O¹⁶O¹⁸O) and asymmetric (¹⁸O¹⁸O¹⁶O) ozone isotopomers for the total angular momentum up to J = 5. Together with similar data recently reported for the singly substituted symmetric (¹⁶O¹⁸O¹⁶O) and asymmetric (¹⁶O¹⁶O¹⁸O) ozone isotopomers, these calculations quantify the role of the Coriolis coupling effect in the large mass-independent isotopic enrichment of ozone, observed in both laboratory experiments and the atmosphere of the Earth. It is found that the Coriolis effect in ozone is relatively small, as evidenced by deviations of its rotational constants from the symmetric-top-rotor behavior, magnitudes of parity splittings (Λ-doubling), and ratios of rovibrational partition functions for asymmetric vs symmetric ozone molecules. It is concluded that all of these characteristics are influenced by the isotopic masses as much as they are influenced by the overall symmetry of the molecule. It is therefore unlikely that the Coriolis coupling effect could be responsible for symmetry-driven mass-independent fractionation of oxygen isotopes in ozone.

Theoretical Treatment of the Coriolis Effect Using Hyperspherical Coordinates, with Application to the Ro-Vibrational Spectrum of Ozone

Several alternative methods for the description of the interaction between rotation and vibration are compared and contrasted using hyperspherical coordinates for a triatomic molecule. These methods differ by the choice of the z-axis and by the assumption of a prolate or oblate rotor shape of the molecule. For each case, a block-structure of the rotational-vibrational Hamiltonian matrix is derived and analyzed, and the advantages and disadvantages of each method are made explicit. This theory is then employed to compute ro-vibrational spectra of singly substituted ozone; roughly, 600 vibrational states of ¹⁶O¹⁸O¹⁶O and ¹⁶O¹⁶O¹⁸O isomers combined, with rotational excitations up to J = 5 and both inversion parities (21600 coupled ro-vibrational states in total). Splittings between the states of different parities, so-called K-doublings, are calculated and analyzed. The roles of the asymmetric-top rotor term and the Coriolis coupling term are determined individually, and it is found that they both affect these splittings, but in the opposite directions. Thus, the two effects partially cancel out, and the residual splittings are relatively small. Energies of the ro-vibrational states reported in this work for ¹⁶O¹⁸O¹⁶O and ¹⁶O¹⁶O¹⁸O are in excellent agreement with literature (available for low-vibrational excitation). New data obtained here for the highly excited vibrational states enable the first systematic study of the Coriolis effect in symmetric and asymmetric isotopomers of ozone.

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

Knowledge of highly excited rovibrational states of ozone isotopologues is of key importance for modelling the dynamics of exchange reactions, for understanding longstanding problems related to isotopic anomalies of the ozone formation, and for analyses of extra-sensitive laser spectral experiments currently in progress. This work is devoted to new theoretical study of high-energy states for the main isotopologue 48O3 = 16O16O16O and for the family of 18O-enriched isotopomers 50O3 = {16O16O18O, 16O18O16O, 18O16O16O} of the ozone molecule considered using a full-symmetry approach. Energies and wave functions of bound states near the dissociation threshold are computed in hyperspherical coordinates accounting for the permutation symmetry of three identical nuclei in 48O3 and of two identical nuclei in 50O3, using the most accurate potential energy surface available now. The obtained vibrational band centers agree with observed ones with the root-mean-squares deviation of about 1 cm-1, making the results appropriate for assignments and analyses of future experimental spectra. The levels delocalized between the three potential wells of ozone isomers are computed and analyzed. The states situated deep in the three (for 48O3) or two (for 50O3) equivalent potential wells have similar energies with negligible splitting. However, the states situated just below the potential barriers separating the wells, are split due to the tunneling between the wells resulting in the splitting of rovibrational sub-bands. We evaluate the amplitudes of the corresponding effects and consider possible perturbations in vibration-rotation bands due to interactions between three potential wells. Theoretical predictions for the splitting of observable band centers are provided for the first time.

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.

Use of Isotope Effects To Understand the Present and Past of the Atmosphere and Climate and Track the Origin of Life

Stable isotope ratio measurements have been used as a measure of a wide variety of processes, including solar system evolution, geological formational temperatures, tracking of atmospheric gas and aerosol chemical transformation, and is the only means by which past global temperatures may be determined over long time scales. Conventionally, isotope effects derive from differences of isotopically substituted molecules in isotope vibrational energy, bond strength, velocity, gravity, and evaporation/condensation. The variations in isotope ratio, such as ¹⁸ O/¹⁶ O (δ¹⁸ O) and ¹⁷ O/¹⁶ O (δ¹⁷ O) are dependent upon mass differences with δ¹⁷ O/δ¹⁸ O=0.5, due to the relative mass differences (1 amu vs. 2 amu). Relations that do not follow this are termed mass independent and are the focus of this Minireview. In chemical reactions such as ozone formation, a δ¹⁷ O/δ¹⁸ O=1 is observed. Physical chemical models capture most parameters but differ in basic approach and are reviewed. The mass independent effect is observed in atmospheric species and used to track their chemistry at the modern and ancient Earth, Mars, and the early solar system (meteorites).

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.