Communication: An accurate global potential energy surface for the ground electronic state of ozone

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.

Classical and quantum trial wave functions in auxiliary-field quantum Monte Carlo applied to oxygen allotropes and a CuBr2 model system

In this work, we test a recently developed method to enhance classical auxiliary-field quantum Monte Carlo (AFQMC) calculations with quantum computers against examples from chemistry and material science, representative of classes of industry-relevant systems. As molecular test cases, we calculate the energy curve of H4 and the relative energies of ozone and singlet molecular oxygen with respect to triplet molecular oxygen, which is industrially relevant in organic oxidation reactions. We find that trial wave functions beyond single Slater determinants improve the performance of AFQMC and allow it to generate energies close to chemical accuracy compared to full configuration interaction or experimental results. In the field of material science, we study the electronic structure properties of cuprates through the quasi-1D Fermi-Hubbard model derived from CuBr2, where we find that trial wave functions with both significantly larger fidelities and lower energies over a mean-field solution do not necessarily lead to AFQMC results closer to the exact ground state energy.

Semiclassical Dynamics on Machine-Learned Coupled Multireference Potential Energy Surfaces: Application to the Photodissociation of the Simplest Criegee Intermediate

Determination of high-dimensional potential energy surfaces (PESs) and nonadiabatic couplings have always been quite challenging. To this end, machine learning (ML) models, trained with a finite set of ab initio data, allow accurate prediction of such properties. To express the PESs in terms of atomic contributions is the cornerstone of any ML based technique because it can be easily scaled to large systems. In this work, we have constructed high fidelity PESs and nonadiabatic coupling terms at the CASSCF level of ab initio data using a machine learning technique, namely, kernel-ridge regression. Additional MRCI-level calculations were carried out to assess the quality of the PESs. We use these machine-learned PESs and nonadiabatic couplings to simulate excited-state molecular dynamics based on Tully's fewest-switches surface hopping method (FSSH). FSSH is a semiclassical method in which nuclei move on the PESs due to the electrons according to the laws of classical mechanics. Nonadiabatic effects are taken into account in terms of transitions between PESs. We apply this scheme to study the O-O photodissociation of the simplest Criegee intermediate (CH₂OO). The FSSH trajectories were initiated on the lowest optically bright singlet excited state (S₂) and propagated along the three most important internal coordinates, namely, O-O and C-O bond distances and the COO bond angle. Some of the trajectories end up on energetically lower PESs as a result of radiationless transfer through conical intersections. All of the trajectories lead to the dissociation of the O-O bond due to the dissociative nature of the excited PESs through one of the two dissociative channels. The simulation reveals that there is about 88.4% probability of dissociation through the lower channel leading to the H₂CO (X¹A₁) and O (¹D) products, whereas there is only 11.6% probability of dissociation through the upper channel leading to H₂CO (a³A″) and O (³P) products.

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.

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.

Quantum Proton Effects from Density Matrix Renormalization Group Calculations

We recently introduced [J. Chem. Phys. 2020, 152, 204103] the nuclear-electronic all-particle density matrix renormalization group (NEAP-DMRG) method to solve the molecular Schrödinger equation, based on a stochastically optimized orbital basis, without invoking the Born-Oppenheimer approximation. In this work, we combine the DMRG method with the nuclear-electronic Hartree-Fock (NEHF-DMRG) approach, treating nuclei and electrons on the same footing. Inter- and intraspecies correlations are described within the DMRG method without truncating the excitation degree of the full configuration interaction wave function. We extend the concept of orbital entanglement and mutual information to nuclear-electronic wave functions and demonstrate that they are reliable metrics to detect strong correlation effects. We apply the NEHF-DMRG method to the HeHHe⁺ molecular ion, to obtain accurate proton densities, ground-state total energies, and vibrational transition frequencies by comparison with state-of-the-art data obtained with grid-based approaches and modern configuration interaction methods. For HCN, we improve on the accuracy of the latter approaches with respect to both the ground-state absolute energy and proton density, which is a major challenge for multireference nuclear-electronic state-of-the-art methods.

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.

The Low-Lying Electronic States of NO2: Potential Energy and Dipole Surfaces, Bound States, and Electronic Absorption Spectrum

Nitrogen dioxide, NO₂, is a free radical composed of the two most abundant elements in Earth's atmosphere, nitrogen and oxygen, and is relevant to atmospheric and combustion chemistry. The electronic structure of even its lowest-lying states is remarkably complex, with various conical intersections and Renner-Teller pairings, giving rise to complex and perturbed vibronic states. Here we report some analysis of the 18 molecular states of doublet spin-multiplicity formed by combining ground-state N(⁴S_u) and O(³P_g) atoms. Three-dimensional potential energy surfaces were fit at the MRCI(Q)-F12/VTZ-F12 level, describing the lowest four (X̃, Ã, B̃, and C̃) electronic states. A properties-based diabatization procedure was applied to accommodate the intersections, producing energies in a quasidiabatic representation and yielding couplings that were also fit into surfaces. The low-lying vibrational levels on the ground X̃ state were computed and compared with experimental measurements. Compared to experiment, the lowest 125 calculated vibrational levels (up to 8500 cm^-1 above the zero-point energy) have a root-mean-squared error of 16.5 cm^-1. In addition, dipole moments for each of the lowest four electronic states-and the transition dipoles between them-were also computed and fit. With the coupled energy and dipole surfaces, the electronic spectrum was calculated in absolute intensity and compared with experimental measurements. Detailed structure in the experimental spectrum was successfully reproduced, and the total integrated intensity matches experiment to an accuracy of ∼1.5% with no empirical adjustments.

On neglecting Coriolis and related couplings in first-principles rovibrational spectroscopy: Considerations of symmetry, accuracy, and simplicity. II. Case studies for H2O isotopologues, H3+, O3, and NH3

For centuries, it has been known that vibrational and rotational degrees of freedom are in general not separable. Nevertheless, surprisingly little is known about the best strategies for approximately separating these degrees of freedom in practice-even in the case of semirigid molecules, where the separation is most meaningful. There is also some confusion in the literature about the proper way to quantify the magnitude of the Coriolis (i.e., rotation-vibration) coupling in rovibrational Hamiltonians or its effect on the rovibrational eigenenergies. In this study, a vibrational-coordinate-independent metric is proposed to quantify the magnitude of the Coriolis contribution to the rovibrational Hamiltonian. The impact of Coriolis coupling on the rovibrational eigenenergies is computed numerically exactly, using both full and various truncated Hamiltonians. The role played by the choice of the vibrational coordinate system-and especially by the choice of "embedding" or body-fixed frame-is examined extensively, both numerically and analytically. This investigation targets several molecular prototypes, all of which serve as important benchmarks for the high-resolution spectroscopic community. Most of these are triatomic molecules, including water (H₂¹⁶O), its deuterated isotopologues (D₂¹⁶O and HD¹⁶O), H₃⁺, and ozone (¹⁶O₃), but the tetratomic ammonia molecule (¹⁴NH₃) is also investigated. These studies provide important insight into the nature of Coriolis coupling under various circumstances. The findings of this study also have significant practical ramifications, vis-à-vis the use of simplifying numerical approximation techniques in nuclear-motion computations.

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.

Nuclear-electronic all-particle density matrix renormalization group

We introduce the Nuclear-Electronic All-Particle Density Matrix Renormalization Group (NEAP-DMRG) method for solving the time-independent Schrödinger equation simultaneously for electrons and other quantum species. In contrast to the already existing multicomponent approaches, in this work, we construct from the outset a multi-reference trial wave function with stochastically optimized non-orthogonal Gaussian orbitals. By iterative refining of the Gaussians' positions and widths, we obtain a compact multi-reference expansion for the multicomponent wave function. We extend the DMRG algorithm to multicomponent wave functions to take into account inter- and intra-species correlation effects. The efficient parameterization of the total wave function as a matrix product state allows NEAP-DMRG to accurately approximate the full configuration interaction energies of molecular systems with more than three nuclei and 12 particles in total, which is currently a major challenge for other multicomponent approaches. We present the NEAP-DMRG results for two few-body systems, i.e., H₂ and H₃ ⁺, and one larger system, namely, BH₃.

The generality of the GUGA MRCI approach in COLUMBUS for treating complex quantum chemistry

The core part of the program system COLUMBUS allows highly efficient calculations using variational multireference (MR) methods in the framework of configuration interaction with single and double excitations (MR-CISD) and averaged quadratic coupled-cluster calculations (MR-AQCC), based on uncontracted sets of configurations and the graphical unitary group approach (GUGA). The availability of analytic MR-CISD and MR-AQCC energy gradients and analytic nonadiabatic couplings for MR-CISD enables exciting applications including, e.g., investigations of π-conjugated biradicaloid compounds, calculations of multitudes of excited states, development of diabatization procedures, and furnishing the electronic structure information for on-the-fly surface nonadiabatic dynamics. With fully variational uncontracted spin-orbit MRCI, COLUMBUS provides a unique possibility of performing high-level calculations on compounds containing heavy atoms up to lanthanides and actinides. Crucial for carrying out all of these calculations effectively is the availability of an efficient parallel code for the CI step. Configuration spaces of several billion in size now can be treated quite routinely on standard parallel computer clusters. Emerging developments in COLUMBUS, including the all configuration mean energy multiconfiguration self-consistent field method and the graphically contracted function method, promise to allow practically unlimited configuration space dimensions. Spin density based on the GUGA approach, analytic spin-orbit energy gradients, possibilities for local electron correlation MR calculations, development of general interfaces for nonadiabatic dynamics, and MRCI linear vibronic coupling models conclude this overview.

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.

Rotationally inelastic scattering of O3–Ar: state-to-state rates with the multiconfigurational time dependent Hartree method

The rates of state-changing collisions are compared for different isotopologues of ozone from quantum scattering calculations with the MCTDH method.

Mechanism Change in the Dynamics of the O' + O2 → O'O + O Atom Exchange Reaction at High Collision Energies

The extreme velocity and the large available energy of atoms with hyperthermal kinetic energies can give rise to novel mechanisms and behavior of chemical reactions unseen at thermal conditions. Crossed-molecular-beams experiments combined with isotope labeling on the reaction of hyperthermal O atoms with O₂ molecules have provided an example of the arising complexity of such systems. Quasiclassical trajectory (QCT) calculations proved to be instructive in the exploration of the microscopic mechanism of the reactive and inelastic scattering observed, and a new mechanism has been identified: there are reactive collisions in which the potential energy remains repulsive during the entire encounter ("direct" reactions in which, in a sense, no complex is formed). In this work, the effect of the magnitude of the collision energy on this mechanism is explored. At hyperthermal collision energies, the reaction is characterized by a unique impact parameter window favorable for reaction through complex formation, while the direct collisions take place exclusively at small impact parameters. In direct reactive collisions, contributing as much as 12% to the reaction cross section, first the existing bond is broken, and the new bond is formed afterward. This kind of collision is unique to extremely high collision energies. Analysis of various correlations was used to find out the details of the reaction dynamics. The observed phenomena indicate that when the collision energy is extremely high, one can expect deviation from what an extrapolation from the more familiar energy ranges would predict.

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.

Calculation of Molecular Vibrational Spectra on a Quantum Annealer

Until recently molecular energy calculations using quantum computing hardware have been limited to gate-based quantum computers. In this paper, a new methodology is presented to calculate the vibrational spectrum of a molecule on a quantum annealer. The key idea of the method is a mapping of the ground state variational problem onto an Ising or quadratic unconstrained binary optimization (QUBO) problem by expressing the expansion coefficients using spins or qubits. The algorithm is general and represents a new revolutionary approach for solving the real symmetric eigenvalue problem on a quantum annealer. The method is applied to two chemically important molecules: O₂ (oxygen) and O₃ (ozone). The lowest two vibrational states of these molecules are computed using both a hardware quantum annealer and a software based classical annealer. Extension of the algorithm to higher dimensions is explicitly demonstrated for an N-dimensional harmonic oscillator (N ≤ 5). The algorithm scales exponentially with dimensionality if a direct product basis is used but will exhibit polynomial scaling for a nondirect product basis.

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.

Anab initiobased full-dimensional potential energy surface for OH + O2⇄ HO3and low-lying vibrational levels of HO3

A full-dimensional potential energy surface for HO₃, including the HO + O₂dissociation asymptote, is developed and rigorous quantum dynamics calculations based on this PES have been carried out to compute the vibrational energy levels of HO₃.

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.

Spectroscopic studies of ozone in cryosolutions: FT-IR spectra of 16O3 in liquid nitrogen, oxygen, argon and krypton

We have measured and interpreted the IR spectra of ozone dissolved in liquid nitrogen, oxygen, argon, and krypton in the 650-4700cm^-1 spectral region at 79-117K. Frequency shifts, band intensities and bandshapes of 22 spectral features of soluted ozone were analyzed. The bands of the А₁ symmetry have a complex contour and possess an excess intensity with respect to the value of the purely vibrational transition moment. It was found that this effect is related to the manifestation of the Coriolis interaction. The bandshape distortion manifests itself as an additional intensity from the side of the В₁ symmetry band being an intensity source in the case of the Coriolis interaction.

Global potential energy surface of ground state singlet spin O4

A new global potential energy for the singlet spin state O₄ system is reported using CASPT2/aug-cc-pVTZ ab initio calculations. The geometries for the six-dimensional surface are constructed using a novel point generation scheme that employs randomly generated configurations based on the beta distribution. The advantage of this scheme is apparent in the reduction of the number of required geometries for a reasonably accurate potential energy surface (PES) and the consequent decrease in the overall computational effort. The reported surface matches well with the recently published singlet surface by Paukku et al. [J. Chem. Phys. 147, 034301 (2017)]. In addition to the O₄ PES, the ground state N₄ PES is also constructed using the point generation scheme and compared with the existing PES [Y. Paukku et al., J. Chem. Phys. 139, 044309 (2013)]. The singlet surface is constructed with the aim of studying high energy O₂-O₂ collisions and predicting collision induced dissociation cross section to be used in simulating non-equilibrium aerothermodynamic flows.

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.