Global permutationally invariant potential energy surface for ozone forming reaction

Rovibrational internal energy transfer and dissociation of high-temperature oxygen mixture

This work constructs a rovibrational state-to-state model for the O₂ + O₂ system leveraging high-fidelity potential energy surfaces and quasi-classical trajectory calculations. The model is used to investigate internal energy transfer and nonequilibrium reactive processes in a dissociating environment using a master equation approach, whereby the kinetics of each internal rovibrational state is explicitly computed. To cope with the exponentially large number of elementary processes that characterize reactive bimolecular collisions, the internal states of the collision partner are assumed to follow a Boltzmann distribution at a prescribed internal temperature. This procedure makes the problem tractable, reducing the computational cost to a comparable scale with the O₂ + O system. The constructed rovibrational-specific kinetic database covers the temperature range of 7500-20 000 K. The reaction rate coefficients included in the database are parameterized in the function of kinetic and internal temperatures. Analysis of the energy transfer and dissociation process in isochoric and isothermal conditions reveals that significant departure from the equilibrium Boltzmann distribution occurs during the energy transfer and dissociation phase. Comparing the population distribution of the O₂ molecules against the O₂ + O case demonstrates a more significant extent of nonequilibrium characterized by a more diffuse distribution whereby the vibrational strands are more clearly identifiable. This is partly due to less efficient mixing of the rovibrational states, which results in more diffuse rovibrational distributions in the quasi-steady-state distribution of O₂ + O₂. A master equation analysis for the combined O₂ + O and O₂ + O₂ system reveals that the O₂ + O₂ system governs the early stage of energy transfer, whereas the O₂ + O system takes control of the dissociation dynamics. The findings of the present work will provide a strong physical foundation that can be exploited to construct an improved reduced-order model for oxygen chemistry.

Spectroscopy and Scattering Studies Using Interpolated Ab Initio Potentials

The Born-Oppenheimer potential energy surface (PES) has come a long way since its introduction in the 1920s, both conceptually and in predictive power for practical applications. Nevertheless, nearly 100 years later-despite astonishing advances in computational power-the state-of-the-art first-principles prediction of observables related to spectroscopy and scattering dynamics is surprisingly limited. For example, the water dimer, (H₂O)₂, with only six nuclei and 20 electrons, still presents a formidable challenge for full-dimensional variational calculations of bound states and is considered out of reach for rigorous scattering calculations. The extremely poor scaling of the most rigorous quantum methods is fundamental; however, recent progress in development of approximate methodologies has opened the door to fairly routine high-quality predictions, unthinkable 20 years ago. In this review, in relation to the workflow of spectroscopy and/or scattering studies, we summarize progress and challenges in the component areas of electronic structure calculations, PES fitting, and quantum dynamical calculations.

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.

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.

Accurate CHIPR Potential Energy Surface for the Lowest Triplet State of C3

We report the first global ab initio-based potential energy surface (PES) for ground-state triplet C₃(³A') based on accurate energies extrapolated to the complete basis set (CBS) limit, and using the combined-hyperbolic-inverse-power-representation method for the analytical modeling. By relying on a cost-effective CBS(D,T) protocol, we ensure that the final form reproduces all topographical features of the PES, including its cyclic-linear isomerization barrier, with CBS(5,6)-quality. To partially account for the incompleteness of the N-electron basis and other minor effects, the available accurate experimental data on the relevant diatomics were used to obtain direct-fit curves that replace the theoretical ones in the many-body expansion. Besides describing properly long-range interactions at all asymptotic channels and permutational symmetry by built-in construction, the PES reported here reproduces the proper exothermicities at dissociation regions as well as the spectroscopy of the diatomic fragments. Bound vibrational state calculations in both linear and cyclic isomers have also been carried out, unveiling a good match of the available data on C₃(ã ³Π_u), while assisting with IR band positions for C₃(³A₂^') that may serve as a guide for its laboratory and astronomical detection.

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).

Data sampling scheme for reproducing energies along reaction coordinates in high-dimensional neural network potentials

We propose a data sampling scheme for high-dimensional neural network potentials that can predict energies along the reaction pathway calculated using the hybrid density functional theory. We observed that a data sampling scheme that combined partial geometry optimization of intermediate structures with random displacement of atoms successfully predicted the energies along the reaction path with respect to five chemical reactions: Claisen rearrangement, Diels-Alder reaction, [1,5]-sigmatropic hydrogen shift, concerted hydrogen transfer in the water hexamer, and Cornforth rearrangement.

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.

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.

Dissociation cross section for high energy O2-O2 collisions

Collision-induced dissociation cross section database for high energy O₂-O₂ collisions (up to 30 eV) is generated and published using the quasiclassical trajectory method on the singlet, triplet, and quintet spin ground state O₄ potential energy surfaces. At equilibrium conditions, these cross sections predict reaction rate coefficients that match those obtained experimentally. The main advantage of the cross section database based on ab initio computations is in the study of complex flows with high degree of non-equilibrium. Direct simulation Monte Carlo simulations using the reactive cross section databases are carried out for high enthalpy hypersonic oxygen flow over a cylinder at rarefied ambient conditions. A comparative study with the phenomenological total collision energy chemical model is also undertaken to point out the difference and advantage of the reported ab initio reaction model.

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.

Dissociation cross sections for N2 + N → 3N and O2 + O → 3O using the QCT method

Cross sections for the homo-nuclear atom-diatom collision induced dissociations (CIDs): N₂ + N and O₂ + O are calculated using Quasi-Classical Trajectory (QCT) method on ab initio Potential Energy Surfaces (PESs). A number of studies for these reactions carried out in the past focused on the CID cross section values generated using London-Eyring-Polanyi-Sato PES and seldom listed the CID cross section data. A highly accurate CASSCF-CASPT2 N₃ and a new O₃ global PES are used for the present QCT analysis and the CID cross section data up to 30 eV relative energy are also published. In addition, an interpolating scheme based on spectroscopic data is introduced that fits the CID cross section for the entire ro-vibrational spectrum using QCT data generated at chosen ro-vibrational levels. The rate coefficients calculated using the generated CID cross section compare satisfactorily with the existing experimental and theoretical results. The CID cross section data generated will find an application in the development of a more precise chemical reaction model for Direct Simulation Monte Carlo code simulating hypersonic re-entry flows.

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.

Does the "reef structure" at the ozone transition state towards the dissociation exist? New insight from calculations and ultrasensitive spectroscopy experiments

Since the discovery of anomalies in ozone isotope enrichment, several fundamental issues in the dynamics linked to the shape of the potential energy surface in the transition state region have been raised. The role of the reeflike structure on the minimum energy path is an intricate question previously discussed in the context of chemical experiments. In this Letter, we bring strong arguments in favor of the absence of a submerged barrier from ultrasensitive laser spectroscopy experiments combined with accurate predictions of highly excited vibrations up to nearly 95% of the dissociation threshold.

High-Level, First-Principles, Full-Dimensional Quantum Calculation of the Ro-vibrational Spectrum of the Simplest Criegee Intermediate (CH2OO)

The ro-vibrational spectrum of the simplest Criegee intermediate (CH2OO) has been determined quantum mechanically based on nine-dimensional potential energy and dipole surfaces for its ground electronic state. The potential energy surface is fitted to more than 50 000 high-level ab initio points with a root-mean-square error of 25 cm(-1), using a recently proposed permutation invariant polynomial neural network method. The calculated rotational constants, vibrational frequencies, and spectral intensities of CH2OO are in excellent agreement with experiment. The potential energy surface provides a valuable platform for studying highly excited vibrational and unimolecular reaction dynamics of this important molecule.