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

The Ratio of The Number of States in Asymmetric and Symmetric Ozone Molecules Deviates from The Statistical Value Of 2

Accurate calculations of vibrational states in singly and doubly substituted ozone molecules are carried out, up to dissociation threshold. Analysis of these spectra reveals noticeable deviations from the statistical factor of 2 for the ratio between the number of states in asymmetric and symmetric ozone molecules. It is found that, for the lower energy parts of spectra, the ratio is less than 2 in the singly substituted ozone molecules, but it is more than 2 in the doubly substituted ozone molecules. However, the upper parts of spectra, just below dissociation thresholds, exhibit a different behavior. In this energy range, the singly and doubly substituted ozone molecules behave similar, with the ratio of states in asymmetric and symmetric ozone molecules being more than 2 in both cases. This property may contribute to an explanation of the mysterious η-effect in the ozone forming reaction that favors the formation of the asymmetric ozone molecules

Four Isotope-Labeled Recombination Pathways of Ozone Formation

A theoretical approach is developed for the description of all possible recombination pathways in the ozone forming reaction, without neglecting any process a priori, and without decoupling the individual pathways one from another. These pathways become physically distinct when a rare isotope of oxygen is introduced, such as 18O, which represents a sensitive probe of the ozone forming reaction. Each isotopologue of O3 contains two types of physically distinct entrance channels and two types of physically distinct product wells, creating four recombination pathways. Calculations are done for singly and doubly substituted isotopologues of ozone, eight rate coefficients total. Two pathways for the formation of asymmetric ozone isotopomer exhibit rather different rate coefficients, indicating large isotope effect driven by -difference. Rate coefficient for the formation of symmetric isotopomer of ozone (third pathway) is found to be in between of those two, while the rate of insertion pathway is smaller by two orders of magnitude. These trends are in good agreement with experiments, for both singly and doubly substituted ozone. The total formation rates for asymmetric isotopomers are found to be somewhat larger than those for symmetric isotopomers, but not as much as in the experiment. Overall, the distribution of lifetimes is found to be very similar for the metastable states in symmetric and asymmetric ozone isotopomers

Molecular Dynamics on Quantum Annealers

In this work we demonstrate a practical prospect of using quantum annealers for simulation of molecular dynamics. A methodology developed for this goal, dubbed Quantum Differential Equations (QDE), is applied to propagate classical trajectories for the vibration of the hydrogen molecule in several regimes: nearly harmonic, highly anharmonic, and dissociative motion. The results obtained using the D-Wave 2000Q quantum annealer are all consistent and quickly converge to the analytical reference solution. Several alternative strategies for such calculations are explored and it was found that the most accurate results and the best efficiency are obtained by combining the quantum annealer with classical post-processing (greedy algorithm). Importantly, the QDE framework developed here is entirely general and can be applied to solve any system of first-order ordinary nonlinear differential equations using a quantum annealer

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 16O18O16O and 16O16O18O. 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

Computational Analysis of Vibrational Modes in Tetra-Sulfur Using Dimensionally Reduced Potential Energy Surface

Electronic structure calculations are carried out for the S4 molecule at the CCSD(T)-F12a/VTZ-F12 level of theory to map out its potential energy surface, which possesses a double-well shape with a low-energy barrier. Two degrees of freedom are considered: the distance R and the gearing motion angle α between the two weakly-perturbed S2 dimers, which form S4. Vibrational states are computed on this 2D-surface and assigned quantum numbers based on their energies and the shapes of their wave functions. Two progressions of vibrational states are identified: a long progression of easily assignable states that develop nodes along the ‘channels’ on the surface, and a shorter progression of states that develop nodes across the ‘channels’ and are much harder to assign, due to the double-well effect. Normal mode analysis indicates that these two modes in S4 represent a significant mixture of conventional bending and stretching motions. When the angle α is increased, the lower frequency mode corresponds to stretching of the distance R, while the higher frequency mode corresponds to compression of R. Frequencies of the modes, ∼180 and ∼420 cm−1, are in a qualitative agreement with earlier ab initio studies of tetra-sulfur, and with sparse experimental data

Influence of the Coriolis Effect on the Properties of Scattering Resonances in Symmetric and Asymmetric Isotopomers of Ozone

Scattering resonances above dissociation threshold are computed for four isotopically substituted ozone species: 16O18O16O, 16O16O18O, 18O16O18O and 16O18O18O, using a variational method with accurate treatment of the rotation–vibration coupling terms (Coriolis effect) for all values of the total angular momentum J from 0 to 4. To make these calculations numerically affordable, a new approach was developed which employs one vibrational basis set optimized for a typical rotational excitation (J,Λ), to run coupled rotation–vibration calculations at several desired values of J. In order to quantify the effect of Coriolis coupling, new data are contrasted with those computed using the symmetric-top rotor approximation, where the rotation–vibration coupling terms are neglected. It is found that, overall, the major properties of scattering resonances (such as their lifetimes, the number of these states, and their cumulative partition function Q) are all influenced by the Coriolis effect and this influence grows as the angular momentum J is raised. However, it is found that the four isotopically substituted ozone molecules are affected roughly equally by the Coriolis coupling. When the ratio η of partition functions for asymmetric over symmetric ozone molecules is computed, the Coriolis effect largely cancels, and this cancelation seems to occur for all values of J. Therefore, it does not seem grounded to attribute any appreciable mass-independent symmetry-driven isotopic fractionation to the Coriolis coupling effect

Methods of instrument testing of smoke detectors performance

This article is devoted to the currently relevant task - determining the real operability of the smoke detector and developing a method for the rapid diagnosis of fire detectors included in the fire alarm loop. In the process, the sensitivity of the smoke optoelectronic fire detector was checked

SpectrumSDT: A Program for Parallel Calculation of Coupled Rotational-Vibrational Energies and Lifetimes of Bound States and Scattering Resonances in Triatomic Systems

We present SpectrumSDT – a program for calculations of energies and lifetimes of bound rotational-vibrational states below and scattering resonances above the dissociation threshold on a global potential energy surface of a triatomic system, which may include stable molecules, weekly-bound van-der-Waals complexes, and unbound atom + diatom scattering systems. Large-amplitude vibrational motion is treated explicitly using hyper-spherical coordinates. Three options for the rotational-vibrational interaction are supported: uncoupled (symmetric top rotor), partially coupled (to include interaction between several nearest states only) and full-coupled (vibrating asymmetric-top rotor). In addition to energies and lifetimes, SpectrumSDT is able to integrate ro-vibrational wave functions over the user-defined regions of potential energy surface, which helps to classify these states. In this release of the code, SpectrumSDT is limited to ABA-type molecules with wave functions that do not extend into the regions near Eckart singularities