Kinetic isotope effect of the 16O + 36O2 and 18O + 32O2 isotope exchange reactions: Dominant role of reactive resonances revealed by an accurate time-dependent quantum wavepacket study

Diabatic Potential Energy Surfaces of SrH2+ and Dynamics Studies of the Sr+(5s2S) + H2 Reaction

Based on the ab initio energy points of both ground and excited states, a neural network fitting method combined with a specific function was successfully used to construct the diabatic potential energy surfaces (PESs) of the SrH2+ system. The topographical features of the diabatic PESs were examined in detail. The results indicate that the nonadiabatic transition characteristics between ground and excited states are accurately described by the newly constructed diabatic PESs. To verify the validity and applicability of the diabatic PESs, as well as the nonadiabatic effects during the reaction process, the quantum dynamics studies of the Sr+(5s2S) + H2 reaction were performed based on both adiabatic and diabatic PESs. The dynamics results indicate that adiabatic dynamics results are dozens of times larger than those of nonadiabatic. This illustrates the significant effect of nonadiabaticity, indicating that adiabatic dynamics results often overestimate the actual values. The integral cross sections (ICSs) were calculated and compared with the experimental data. The diabatic ICSs are in good agreement with the experimental results. The reasonable dynamics results indicate that the newly constructed diabatic PESs are suitable for the relevant dynamics studies.

Construction of diabatic potential energy surfaces for the SiH2+ system and dynamics studies of the Si+(2P1/2, 3/2) + H2 reaction

The construction of diabatic potential energy surfaces (PESs) for the SiH2+ system, related to the ground (12A') and excited states (22A'), has been successfully achieved. This was accomplished by utilizing high-level ab initio energy points, employing a neural network fitting method in conjunction with a specifically designed function. The newly constructed diabatic PESs are carefully examined for dynamics calculations of the Si+(2P1/2, 3/2) + H2 reaction. Through time-dependent quantum wave packet calculations, the reaction probabilities, integral cross sections (ICSs), and differential cross sections (DCSs) of the Si+(2P1/2, 3/2) + H2 reaction were reported. The dynamics results indicate that the total ICS is in excellent agreement with experimental data within the collision energy range studied. The results also indicate that the SiH+ ion is hardly formed via the Si+(2P3/2) + H2 reaction. The results from the DCSs suggest that the "complex-forming" reaction mechanism predominates in the low collision energy region. Conversely, the forward abstraction reaction mechanism is dominant in the high collision energy region.

New Diabatic Potential Energy Surfaces for the Li + H2 Reaction and Time-Dependent Quantum Wave Packet Studies

New global diabatic potential energy surfaces (DPESs) for the ground (12A') and first excited (22A') states for the Li + H2 system were developed, with more than 30,000 energy points at the IC-MRCI+Q level of theory, utilizing the aug-cc-pV5Z basis set for the H atoms and the cc-pCV5Z basis set for the Li atom, fitted by a single neural network (NN) with symmetry. Product state-resolved quantum dynamics calculations of the nonadiabatic reaction Li (2P) + H2 (X 1 ∑g+, v0 = 0, j0 = 0) → LiH (X 1∑+) + H(2S) were carried out using these new DPESs and also the previous HYLC-DPESs. The numerical results suggested that our newly constructed DPESs provided an accurate description of the LiH2 system.

New Potential Energy Surface for the H + Cl2 Reaction and Quantum Dynamics Studies

The reaction of H + Cl2 → HCl + Cl plays a crucial role in various fields. However, no previous study has investigated this reaction using accurate quantum mechanical methods. In this paper, we construct a global potential energy surface (PES) using the neural network method with more than 20,000 ab initio energies obtained by the MRCI-F12+Q method with the aug-cc-pV5Z basis and extrapolated to the complete basis set limit. The spin-orbit coupling of the Cl atom has been considered in the PES. With this new PES, product state-resolved quantum dynamics calculations for the H + Cl2 (v0 = 0, j0 = 0-2) → HCl + Cl reaction was carried out. Numerical results show that the initial rotational excitation of the Cl2 has negligible effects on the reactivity. Product state-resolved integral cross sections (ICS) and rate constants reveal that the HCl is most favorably produced in its v' = 2 vibrational state. The calculated product vibrational state-resolved and total reaction rate constants suggest that the new global PES is accurate enough, as compared with the available experimental measurements.

Product State-Resolved Reactive Scattering Studies of the H + Cl2 (v0 = 1-3, j0 = 0) → HCl + Cl Reaction by the Time-Dependent Wave Packet Method

The typical hydrogen atom plus halogen molecule reaction H + Cl2 → HCl + Cl has implications across many fields. In this paper, product state-resolved quantum dynamics calculations for the vibrationally excited reaction H + Cl2 (v0 = 1-3, j0 = 0) → HCl + Cl were conducted using the time-dependent wave packet method on a newly developed accurate potential energy surface. Numerical results indicate that the initial vibrational excitation of Cl2 does enhance the reactivity for this early barrier reaction, although less than the enhancement of the translational energy. The calculated product vibrational state-resolved integral cross sections and rate constants reveal that the product vibrational state distribution and the initial vibrational state of Cl2 are highly correlated. The thermal rate constant in the temperature range from 100 to 1000 K was given and is found to be in reasonable agreement with the experimental measurements.

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.

Nonadiabatic dynamics studies of the H(2S) + RbH(X1Σ+) reaction: based on new diabatic potential energy surfaces

The global diabatic potential energy surfaces (PESs) that correspond to the ground (1²A') and first excited states (2²A') of the RbH₂ system PES are constructed based on 17 786 ab initio points. The neural network method is used to fit the PESs and the topographic features of the new diabatic PESs are discussed in detail. Based on the newly constructed diabatic PESs, the dynamics calculations of the H(²S) + RbH(X¹Σ⁺) → Rb(5²S) + H₂(X¹Σ_g ⁺)/Rb(5²P) + H₂(X¹Σ_g ⁺) reactions are performed using the time-dependent wave packet method. The dynamics properties of these two channels such as the reaction probabilities, integral cross sections, and differential cross sections (DCSs) are calculated at state-to-state level of theory. The nonadiabatic effects are discussed in detail, and the results indicate that the adiabatic results are overestimated from the dynamics values. The DCSs of these two channels are forward biased, which indicates that the abstraction mechanism plays a dominant role in the reaction.

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.

Inelastic scattering in isotopologues of O2-Ar: the effects of mass, symmetry, and density of states

The two species considered here, O2 (oxygen molecule) and Ar (argon-atom), are both abundant components of Earth's atmosphere and hence familiar collision partners in this medium. O2 is quite reactive and extensively involved in atmospheric chemistry, including Chapman's cycle of the formation and destruction of ozone; while Ar, like N2, typically plays the nevertheless crucial role of inert collider. Inert species can provide stabilization to metastable encounter-complexes through the energy transfer associated with inelastic collisions. The interplay of collision frequency and energy transfer efficiency, with state lifetimes and species concentrations, contributes to the rich and varied chemistry and dynamics found in diverse environments ranging from planetary atmospheres to the interstellar and circumstellar media. The nature and density of bound and resonance states, coupled electronic states, symmetry, and nuclear spin-statistics can all play a role. Here, we systematically investigate some of those factors by looking at the O2-Ar system, comparing rigorous quantum-scattering calculations for the 16O16O-40Ar, 18O16O-40Ar, and 18O18O-40Ar isotope combinations. A new accurate potential energy surface was constructed for this purpose holding the O2 bond distance at its vibrationally averaged distance.

A new global potential energy surface of the SH2+(X4A'') system and quantum calculations for the S+ + H2(v = 0-3, j = 0) reaction

A new global potential energy surface (PES) for the ground state of the SH2+(X4A'') system is constructed using a permutation invariant polynomial neural network method. In ab initio calculations, the MRCI-F12 method with the AVTZ basis set is used. Furthermore, the dynamics calculations of the S+ + H2(v = 0-3, j = 0) → SH+ + H reaction are carried out based on the new PES. The reaction probabilities and integral cross sections are compared with available theoretical calculations. Present values are in general good agreement with the previous theoretical studies. However, some discrepancies can still be found due to different PESs used in the calculation. Furthermore, the vibrational energy of the reactant molecule can significantly enhance the reactivity compared to the translational energy. The differential cross sections indicated that the reaction mechanism is changed from the "head-on" rebound mechanism to the tripping mechanism with the increasing number of initial vibrational excitation state.

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.

Time-dependent quantum mechanical wave packet dynamics

Starting from a model study of the collinear (H, H₂) exchange reaction in 1959, the time-dependent quantum mechanical wave packet (TDQMWP) method has come a long way in dealing with systems as large as Cl + CH₄. The fast Fourier transform method for evaluating the second order spatial derivative of the wave function and split-operator method or Chebyshev polynomial expansion for determining the time evolution of the wave function for the system have made the approach highly accurate from a practical point of view. The TDQMWP methodology has been able to predict state-to-state differential and integral reaction cross sections accurately, in agreement with available experimental results for three dimensional (H, H₂) collisions, and identify reactive scattering resonances too. It has become a practical computational tool in predicting the observables for many A + BC exchange reactions in three dimensions and a number of larger systems. It is equally amenable to determining the bound and quasi-bound states for a variety of molecular systems. Just as it is able to deal with dissociative processes (without involving basis set expansion), it is able to deal with multi-mode nonadiabatic dynamics in multiple electronic states with equal ease. We present an overview of the method and its strength and limitations, citing examples largely from our own research groups.

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.

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.

Dynamics studies of the H + HBr reaction: Based on a new potential energy surface

The initial state specific quantum wave packet dynamics studies of the H + HBr (v₀ = 0, j₀ = 0-2) reaction were performed using a new global potential energy surface (PES) of the ground state of the BrH₂ system for the collision energy ranging from 0.01 to 2.0 eV. The PES was constructed using the permutation invariant polynomial neural network method based on approximately 63 000 ab initio points, which were calculated by the multireference configuration interaction method with AVTZ and AVQZ basis sets. To improve the accuracy of the PES, Davidson's correction and spin-orbit coupling effects were considered in the ab initio calculation and the basis set was extrapolated to complete basis set limit. The new PES was compared with the previous ones and also the available experimental data, which suggests that the new PES is more accurate. The state-to-state quantum wave packet dynamics was carried out using the reactant-coordinate based approach. The reaction probabilities, integral and differential cross sections, rovibrational state distributions of product and rate constants, etc., were compared with the available theoretical and experimental studies. In general, the present work is in better agreement with the available experimental data. The quantum dynamics studies suggest that the rotational excitation of HBr has little effect on the reaction.

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.

Fully coupled (J > 0) time-dependent wave-packet calculations using hyperspherical coordinates for the H + O2 reaction on the CHIPR potential energy surface

ICS calculation by time dependent wavepacket approach for H + O₂ reaction using non-zero J values.

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.

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.

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.

A new potential energy surface of the OH2+ system and state-to-state quantum dynamics studies of the O+ + H2 reaction

A new global potential energy surface of the O⁺ + H₂ system was constructed with neural network method, using about 63000 ab initio points, which were calculated by employing the multi-reference configuration interaction method with aug-cc-pVTZ and aug-cc-pVQZ basis sets.

Recombination reactions as a possible mechanism of mass-independent fractionation of sulfur isotopes in the Archean atmosphere of Earth

A hierarchy of isotopically substituted recombination reactions is formulated for production of sulfur allotropes in the anoxic atmosphere of Archean Earth. The corresponding system of kinetics equations is solved analytically to obtain concise expressions for isotopic enrichments, with focus on mass-independent isotope effects due to symmetry, ignoring smaller mass-dependent effects. Proper inclusion of atom-exchange processes is shown to be important. This model predicts significant and equal depletions driven by reaction stoichiometry for all rare isotopes: ³³S, ³⁴S, and ³⁶S. Interestingly, the ratio of capital [Formula: see text] values obtained within this model for ³³S and ³⁶S is -1.16, very close to the mass-independent fractionation line of the Archean rock record. This model may finally offer a mechanistic explanation for the striking mass-independent fractionation of sulfur isotopes that took place in the Archean atmosphere of Earth.

On stabilization of scattering resonances in recombination reaction that forms ozone

Calculations of energy transfer in the recombination reaction that forms ozone are carried out within the framework of the mixed quantum/classical theory and using the dimensionally reduced 2D-model of ozone molecule, with bending motion neglected. Recombination rate coefficients are obtained at room temperature for symmetric and asymmetric isotopomers of singly and doubly substituted isotopologues. The processes of resonance formation, spontaneous decay, collisional dissociation, and stabilization by bath gas (Ar) are all characterized and taken into account within the steady-state approximation for kinetics. The focus is on stabilization step, where the mysterious isotopic η-effect was thought to originate from. Our results indicate no difference in cross sections for stabilization of scatteringresonances in symmetric and asymmetric isotopomers. As practical results, the general and simple analytic models for stabilization and dissociation cross sections are presented, which can be applied to resonances in any ozone molecule, symmetric or asymmetric, singly or doubly substituted. Present calculations show some isotope effect that looks similar to the experimentally observed η-effect, and the origin of this phenomenon is in the rates of formation/decay of scatteringresonances, determined by their widths, that are somewhat larger in asymmetric isotopomers than in their symmetric analogues. However, the approximate two-dimensional model used here is insufficient for consistent and reliable description of all features of the isotopic effect in ozone. Calculations using an accurate 3D model are still needed.

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.

Quantum Dynamics of the (18)O + (36)O2 Collision Process

We report full quantum cross sections and rate constants for the (18)O + (36)O2 → (36)O2 + (18)O collision process. This constitutes to the best of our knowledge the first dynamical study of the (18)O(18)O(18)O system, with three identical (18)O oxygen atoms. We emphasize the comparison with the (16)O + (32)O2 collision as this latter presents the exact same features as the one treated here, except the consistent change of mass for all three atoms. We find very similar behaviors in the cross sections, and we confirm that the rates are faster when three identical nuclei are involved. In particular, we cannot dynamically study this system with classical trajectory methods, and we have to include properly the indistinguishability of the three (18)O nuclei; however, we note some slight differences with the (16)O(16)O(16)O benchmark system, and we focus our analysis on their origin.