1
|
Kokoouline V, Alijah A, Tyuterev V. Lifetimes and decay mechanisms of isotopically substituted ozone above the dissociation threshold: matching quantum and classical dynamics. Phys Chem Chem Phys 2024; 26:4614-4628. [PMID: 38251711 DOI: 10.1039/d3cp04286c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2024]
Abstract
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.
Collapse
Affiliation(s)
| | - Alexander Alijah
- Groupe de Spectrometrie Moléculaire et Atmospherique, UMR CNRS 7331, University of Reims Champagne-Ardenne, Reims Cedex 2, F-51687, France
| | - Vladimir Tyuterev
- Laboratory of Molecular Quantum Mechanics and Radiative transfer, Tomsk State University, Tomsk, Russia
- Laboratory of Theoretical Spectroscopy, V.E. Zuev Institute of Atmospheric Optics, Russian Academy of Sciences, Tomsk, 634055, Russia
| |
Collapse
|
2
|
Gayday I, Babikov D. Efficient Method for an Approximate Treatment of the Coriolis Effect in Calculations of Quantum Dynamics and Spectroscopy, with Application to Scattering Resonances in Ozone. J Phys Chem A 2021; 125:5661-5669. [PMID: 34156247 DOI: 10.1021/acs.jpca.1c03350] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
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.
Collapse
Affiliation(s)
- Igor Gayday
- Department of Chemistry, Marquette University, Wehr Chemistry Building, Milwaukee, Wisconsin 53201-1881, United States
| | - Dmitri Babikov
- Department of Chemistry, Marquette University, Wehr Chemistry Building, Milwaukee, Wisconsin 53201-1881, United States
| |
Collapse
|
3
|
Sarka J, Poirier B. Comment on "Calculated vibrational states of ozone up to dissociation" [J. Chem. Phys. 144, 074302 (2016)]. J Chem Phys 2020; 152:177101. [PMID: 32384842 DOI: 10.1063/5.0002762] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- János Sarka
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409, USA
| | - Bill Poirier
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409, USA
| |
Collapse
|
4
|
Gayday I, Teplukhin A, Kendrick BK, Babikov D. The role of rotation-vibration coupling in symmetric and asymmetric isotopomers of ozone. J Chem Phys 2020; 152:144104. [PMID: 32295370 DOI: 10.1063/1.5141365] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
A theoretical framework and a computer code (SpectrumSDT) are developed for accurate calculations of coupled rotational-vibrational states in triatomic molecules using hyper-spherical coordinates and taking into account the Coriolis coupling effect. Concise final formulas are derived for the construction of the Hamiltonian matrix using an efficient combination of the variational basis representation and discrete variable representation methods with locally optimized basis sets and grids. First, the new code is tested by comparing its results with those of the APH3D program of Kendrick et al. [Kendrick, Pack, Walker, and Hayes, J. Chem. Phys. 110, 6673 (1999)]. Then, accurate calculations of the rovibrational spectra are carried out for doubly substituted symmetric (18O16O18O) and asymmetric (18O18O16O) ozone isotopomers for the total angular momentum up to J = 5. Together with similar data recently reported for the singly substituted symmetric (16O18O16O) and asymmetric (16O16O18O) ozone isotopomers, these calculations quantify the role of the Coriolis coupling effect in the large mass-independent isotopic enrichment of ozone, observed in both laboratory experiments and the atmosphere of the Earth. It is found that the Coriolis effect in ozone is relatively small, as evidenced by deviations of its rotational constants from the symmetric-top-rotor behavior, magnitudes of parity splittings (Λ-doubling), and ratios of rovibrational partition functions for asymmetric vs symmetric ozone molecules. It is concluded that all of these characteristics are influenced by the isotopic masses as much as they are influenced by the overall symmetry of the molecule. It is therefore unlikely that the Coriolis coupling effect could be responsible for symmetry-driven mass-independent fractionation of oxygen isotopes in ozone.
Collapse
Affiliation(s)
- Igor Gayday
- Department of Chemistry, Wehr Chemistry Building, Marquette University, Milwaukee, Wisconsin 53201-1881, USA
| | - Alexander Teplukhin
- Theoretical Division (T-1, MS B221), Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| | - Brian K Kendrick
- Theoretical Division (T-1, MS B221), Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| | - Dmitri Babikov
- Department of Chemistry, Wehr Chemistry Building, Marquette University, Milwaukee, Wisconsin 53201-1881, USA
| |
Collapse
|
5
|
Gayday I, Teplukhin A, Kendrick BK, Babikov D. Theoretical Treatment of the Coriolis Effect Using Hyperspherical Coordinates, with Application to the Ro-Vibrational Spectrum of Ozone. J Phys Chem A 2020; 124:2808-2819. [PMID: 32227893 DOI: 10.1021/acs.jpca.0c00893] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
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 16O18O16O and 16O16O18O 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 16O18O16O and 16O16O18O 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.
Collapse
Affiliation(s)
- Igor Gayday
- Department of Chemistry, Wehr Chemistry Building, Marquette University, Milwaukee, Wisconsin 53201-1881, United States
| | - Alexander Teplukhin
- Theoretical Division (T-1, MS B221), Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Brian K Kendrick
- Theoretical Division (T-1, MS B221), Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Dmitri Babikov
- Department of Chemistry, Wehr Chemistry Building, Marquette University, Milwaukee, Wisconsin 53201-1881, United States
| |
Collapse
|
6
|
Gayday I, Grushnikova E, Babikov D. Influence of the Coriolis effect on the properties of scattering resonances in symmetric and asymmetric isotopomers of ozone. Phys Chem Chem Phys 2020; 22:27560-27571. [PMID: 33236748 DOI: 10.1039/d0cp05060a] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
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.
Collapse
Affiliation(s)
- Igor Gayday
- Department of Chemistry, Marquette University, Wehr Chemistry Building, Milwaukee, Wisconsin 53201-1881, USA.
| | | | | |
Collapse
|
7
|
Kokoouline V, Lapierre D, Alijah A, Tyuterev V. Localized and delocalized bound states of the main isotopologue 48O 3 and of 18O-enriched 50O 3 isotopomers of the ozone molecule near the dissociation threshold. Phys Chem Chem Phys 2020; 22:15885-15899. [PMID: 32642747 DOI: 10.1039/d0cp02177f] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
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.
Collapse
Affiliation(s)
| | - David Lapierre
- Groupe de Spectrometrie Moléculaire et Atmospherique, UMR CNRS 7331, University of Reims Champagne-Ardenne, F-51687, Reims Cedex 2, France.
| | - Alexander Alijah
- Groupe de Spectrometrie Moléculaire et Atmospherique, UMR CNRS 7331, University of Reims Champagne-Ardenne, F-51687, Reims Cedex 2, France.
| | - Vladimir Tyuterev
- Groupe de Spectrometrie Moléculaire et Atmospherique, UMR CNRS 7331, University of Reims Champagne-Ardenne, F-51687, Reims Cedex 2, France. and Quamer Laboratory, Tomsk State University, Tomsk, Russia
| |
Collapse
|
8
|
Sarka J, Petty C, Poirier B. Exact bound rovibrational spectra of the neon tetramer. J Chem Phys 2019; 151:174304. [PMID: 31703493 DOI: 10.1063/1.5125145] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Exact quantum dynamics calculations are performed for the bound rovibrational states of the neon tetramer (Ne4) in its ground electronic state, using pair-wise Lennard-Jones potentials and the ScalIT suite of parallel codes. The vibrational states separate into a low-lying group mostly localized to a single potential well and a higher-energy delocalized group lying above the isomerization threshold-with such a structure serving as a testament to the "intermediate" quantum nature of the Ne4 system. To accurately and efficiently represent both groups of states, the phase-space optimized discrete variable representation (PSO-DVR) approach was used, as implemented in the ScalIT code. The resultant 1D PSO effective potentials also shed significant light on the quantum dynamics of the system. All vibrational states were computed well up into the isomerization band and labeled up to the classical isomerization threshold-defined as the addition of the classical energy of a single bond, ϵ = 24.7 cm-1, to the quantum ground state energy. Rovibrational energy levels for all total angular momentum values in the range J = 1-5 were also computed, treating all Coriolis coupling exactly.
Collapse
Affiliation(s)
- János Sarka
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409-1061, USA
| | - Corey Petty
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409-1061, USA
| | - Bill Poirier
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409-1061, USA
| |
Collapse
|
9
|
Teplukhin A, Kendrick BK, Babikov D. Calculation of Molecular Vibrational Spectra on a Quantum Annealer. J Chem Theory Comput 2019; 15:4555-4563. [PMID: 31314517 DOI: 10.1021/acs.jctc.9b00402] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
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: O2 (oxygen) and O3 (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.
Collapse
Affiliation(s)
- Alexander Teplukhin
- Theoretical Division (T-1, MS B221) , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States
| | - Brian K Kendrick
- Theoretical Division (T-1, MS B221) , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States
| | - Dmitri Babikov
- Department of Chemistry , Marquette University , Milwaukee , Wisconsin 53021 , United States
| |
Collapse
|