Ab initio predictions and laboratory validation for consistent ozone intensities in the MW, 10 and 5 μm ranges

Ozone spectroscopy in the terahertz range from first high-resolution Synchrotron SOLEIL experiments combined with far-infrared measurements and ab initio intensity calculations

Ozone is one of the important molecules in terms of the impact on the atmospheric chemistry, climate changes, bio- and eco-systems and human health. It has a strong absorption in the microwave, terahertz and far-infrared spectral ranges where a large part of the Earth's outgoing longwave radiation to space is located. In this work, the observations, and analyses of the ozone high-resolution spectra in the THz range recorded using the Synchrotron light source of the SOLEIL CNRS equipment are reported for the first time. Thanks to the exceptional brightness of the Synchrotron radiation and to the signal/noise ratio, it was possible to observe many more ozone transitions of the cold rotational band and the hot ν2-ν2 band in the range 0.9-6 THz compared to the previous works. In addition, we have carried out new measurements and assignments for the ν2 band. The simultaneous fit of the rotational band GS-GS, the hot band ν2-ν2 and the FIR ν2 band yielded an overall weighted standard deviation of 0.68 for 13,466 line positions within the experimental accuracy. This includes all previously available MW (with the best uncertainty 0.1 - 10 kHz), FIR data and the original SOLEIL measurements that provided experimental accuracy of 0.00005 - 0.0001 cm-1 for the best lines. Significant deviations in new experimental spectra compared to available spectroscopic databases were evidenced, particularly for the line positions and energy levels at high J, Ka rotational quantum numbers that are the most pronounced in the 4.5 - 6 THz range. Accurate ab initio calculations of line intensities combined with empirically fitted line positions were used to create new linelists that permit theoretical modelling of the transmittance in a good agreement with the Synchrotron spectra in the entire range of observations for various pressures and optical paths. The region near 100 cm-1 and above appears to be more sensitive to the temperature conditions that should be considered in atmospheric observation for the currently operational and future ground based and space missions.

Convergence of series expansions in rovibrational configuration interaction (RVCI) calculations

Rotational and rovibrational spectra are a key in astrophysical studies, atmospheric science, pollution monitoring, and other fields of active research. The ab initio calculation of such spectra is fairly sensitive with respect to a multitude of parameters and all of them must be carefully monitored in order to yield reliable results. Besides the most obvious ones, i.e., the quality of the multidimensional potential energy surface and the vibrational wavefunctions, it is the representation of the μ-tensor within the Watson Hamiltonian, which has a significant impact on the desired line lists or simulated spectra. Within this work, we studied the dependence of high-resolution rovibrational spectra with respect to the truncation order of the μ-tensor within the rotational contribution and the Coriolis coupling operator of the Watson operator. Moreover, the dependence of the infrared intensities of the rovibrational transitions on an n-mode expansion of the dipole moment surface has been investigated as well. Benchmark calculations are provided for thioformaldehyde, which has already served as a test molecule in other studies and whose rovibrational spectrum was found to be fairly sensitive. All calculations rely on rovibrational configuration interaction theory and the discussed high-order terms of the μ-tensor are a newly implemented feature, whose theoretical basics are briefly discussed.

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.

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.

A Decade with VAMDC: Results and Ambitions

This paper presents an overview of the current status of the Virtual Atomic and Molecular Data Centre (VAMDC) e-infrastructure, including the current status of the VAMDC-connected (or to be connected) databases, updates on the latest technological development within the infrastructure and a presentation of some application tools that make use of the VAMDC e-infrastructure. We analyse the past 10 years of VAMDC development and operation, and assess their impact both on the field of atomic and molecular (A&M) physics itself and on heterogeneous data management in international cooperation. The highly sophisticated VAMDC infrastructure and the related databases developed over this long term make them a perfect resource of sustainable data for future applications in many fields of research. However, we also discuss the current limitations that prevent VAMDC from becoming the main publishing platform and the main source of A&M data for user communities, and present possible solutions under investigation by the consortium. Several user application examples are presented, illustrating the benefits of VAMDC in current research applications, which often need the A&M data from more than one database. Finally, we present our vision for the future of VAMDC.

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.

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.