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.

Ozone Formation in Ternary Collisions: Theory and Experiment Reconciled.

The present Letter shows that the formation of ozone in ternary collisions O+O_{2}+M-the primary mechanism of ozone formation in the stratosphere-at temperatures below 200 K (for M=Ar) proceeds through a formation of a temporary complex MO_{2}, while at temperatures above ∼700 K, the reaction proceeds mainly through a formation of long-lived vibrational resonances of O_{3}^{*}. At intermediate temperatures 200-700 K, the process cannot be viewed as a two-step mechanism, often used to simplify and approximate collisions of three atoms or molecules. The developed theoretical approach is applied to the reaction O+O_{2}+Ar because of extensive experimental data available. The rate coefficients for the formation of O_{3} in ternary collisions O+O_{2}+Ar without using two-step approximations were computed for the first time as a function of collision energy. Thermally averaged coefficients were derived for temperatures 5-900 K. It is found that the majority of O_{3} molecules formed initially are weakly bound. Accounting for the process of vibrational quenching of the nascent population, a good agreement with available experimental data for temperatures 100-900 K is obtained.

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.

Towards a complete elucidation of the ro-vibrational band structure in the SF6 infrared spectrum from full quantum-mechanical calculations.

The first accurate and complete theoretical room-temperature rotationally resolved spectra in the range 300-3000 cm-1 are reported for the three most abundant isotopologues (32SF6, 33SF6 and 34SF6) of the sulfur hexafluoride molecule. The literature reports that SF6 is widely used as a prototype molecule for studying the multi-photon excitation processes with powerful lasers in the infrared range. On the other hand, SF6 is an important greenhouse molecule with a very long lifetime in the atmosphere. Because of relatively low vibrational frequencies, the hot bands of this molecule contribute significantly to the absorption infrared spectra even at room temperature. This makes the calculation of complete rovibrational line lists required for fully converged opacity modeling extremely demanding. In order to reduce the computational costs, symmetry was exploited at all stages of the first global variational nuclear motion calculations by means of irreducible tensor operators. More than 2600 new vibrational band centers were predicted using our empirically refined ab initio potential energy surface. Highly excited rotational states were calculated up to J = 121, resulting in 6 billion transitions computed from an ab initio dipole moment surface and distributed over more than 500 cold and hot bands. The final line lists are made available through the TheoReTS information system (http://theorets.univ-reims.fr, http://theorets.tsu.ru). For the first time, the major (ro)vibrational band structures in the wavenumber range corresponding to the strongest absorption in the infra-red are completely elucidated for a seven-atom molecule, providing excellent agreement with the observed spectral patterns. It is shown that the obtained results are more complete than all available line lists, permitting reliable modelling of the temperature dependence of the molecular opacity.

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 16O3 isotopologue and for homogeneous isotopic substitutions 17O3 and 18O3 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 D0. 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 16O3 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 D0/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.

Derivation of ρ-dependent coordinate transformations for nonrigid molecules in the Hougen-Bunker-Johns formalism.

In this paper, we report a series of transformations for the construction of a Hamiltonian model for nonrigid polyatomic molecules in the framework of the Hougen-Bunker-Johns formalism (HBJ). This model is expressed in normal mode coordinates for small vibrations and in a specific coordinate ρ to describe the large amplitude motion. For the first time, a general procedure linking the "true" curvilinear coordinates to ρ is proposed, allowing the expression of the potential energy part in the same coordinate representation as the kinetic energy operator, whatever the number of atoms. A Lie group-based method is also proposed for the derivation of the reference configuration in the internal axis system. This work opens new perspectives for future high-resolution spectroscopy studies of nonrigid, medium-sized molecules using HBJ-type Hamiltonians. Illustrative examples and computation of vibrational energy levels on semirigid and nonrigid molecules are given to validate this method.

First Full-Dimensional Potential Energy and Dipole Moment Surfaces of SF6.

A 15-dimensional analytical form for the potential energy and dipole moment surfaces of the SF6 molecule in the ground electronic state is obtained using ab initio methods. In order to calculate the equilibrium S-F distance, we applied the coupled cluster CCSD(T) method and several versions of the correlation-consistent basis sets from valence triple-zeta (VTZ) and augmented valence triple-zeta (AVTZ) to core-valence quadruple-zeta (CVQZ) with Douglas-Kroll (DK) relativistic corrections that provided good agreement with an empirical equilibrium value. Ab initio electronic energies on 15D grids of nuclear geometries are computed using the CCSD(T) method with VTZ and CVQZ-DK basis sets. The analytical representation of the potential energy surface is determined through an expansion in symmetry-adapted products of nonlinear coordinates up to the 5th order. The influence of additional redundant coordinates on the quality of the fit was investigated. Parameters of full-dimensional dipole moment surfaces are determined up to the 4th order expansion in normal mode coordinates. For validation of ab initio results, the fundamental vibration frequencies and absorption cross sections of the principal sulfur hexafluoride isotopologue are calculated, giving good agreement with cold (180 K) and room temperature (296 K) experimental spectra. Absorption cross sections calculated from our preliminary line list agree much better with these observations than the simulations using SF6 line-by-line lists constructed from effective models included in currently available spectroscopic databases.

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: 16O16O + 18O → 16O18O + 16O. 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.

Symmetry effects in rotationally resolved spectra of bi-deuterated ethylene: Theoretical line intensities of cis, trans, and as-C2H2D2 isotopomers.

In this paper, we report accurate first-principles variational rovibrational spectra predictions for the three double deuterated ethylene isotopologs denoted as cis, trans, and as-12C2H2D2. Calculations were performed in the framework of the normal-mode approach using our ab initio12C2H4 (D2h) Born-Oppenheimer potential energy and dipole moment surfaces. Symmetry breaking effects under bideuterated H → D substitutions (D2h → C2v/C2h) and their impact on infrared spectra are studied from normal mode transformations. All theoretical spectra simulated at 296 K up to J = 38 are in good agreement, both for line positions and in absorption cross sections, with the experiment. Accurate theoretical line lists providing for the first time intensities of rovibrational transitions are computed for the three species in the range (0-4500) cm-1 and will be available on our TheoReTS information system (http://theorets.univ-reims.fr and http://theorets.tsu.ru). These results can be used for remote sensing retrieval of isotopic species using predicted line intensities and experimentally optimized line positions.

Isotopic and symmetry breaking effects on phosphine spectra under H → D substitutions from ab initio variational calculations.

Variationally computed infrared spectra in the range [0-5000] cm-1 are reported for the deuterated PH2D and PHD2 molecules from accurate potential energy and dipole moment surfaces initially derived for the major isotopologue PH3( C 3 v ). Energy level and line intensity calculations were performed by using a normal-mode model combined with isotopic and symmetry transformations for the H → D substitutions. Theoretical spectra were computed at 296 K up to J max = 30 and will be made available through the TheoReTS information system (http://theorets.univ-reims.fr, http://theorets.tsu.ru). For the very first time, ab initio intensity predictions of PH2D/PHD2 are in good qualitative agreement with the literature. This work will be useful for spectral intensity analysis for which accurate spectral intensity data are still missing.

Highly excited vibrational levels of methane up to 10 300 cm-1: Comparative study of variational methods.

In this work, we report calculated vibrational energy levels of the methane molecule up to 10 300 cm-1. Two potential energy surfaces constructed in quite different coordinate systems with different analytical representations are employed in order to evaluate the uncertainty of vibrational predictions. To calculate methane energy levels, we used two independent techniques of the variational method. One method uses an exact kinetic energy operator in internal curvilinear coordinates. Another one uses an expansion of Eckart-Watson nuclear motion Hamiltonian in rectilinear normal coordinates. In the Icosad range (up to five vibrational quanta bands-below 7800 cm-1), the RMS standard deviations between calculated and observed energy levels were 0.22 cm-1 and 0.41 cm-1 for these two quite different approaches. For experimentally well-known 3v3 sub-levels, the calculation accuracy is estimated to be ∼1 cm-1. In the Triacontad range (7660-9188 cm-1), the average error of the calculation is about 0.5 cm-1. The accuracy and convergence issues for higher energy ranges are discussed.

Understanding global infrared opacity and hot bands of greenhouse molecules with low vibrational modes from first-principles calculations: the case of CF4.

Fluorine containing molecules have a particularly long atmospheric lifetime and their very big estimated global warming potentials are expected to rapidly increase in the future. This work is focused on the global theoretical prediction of infrared spectra of the tetrafluoromethane molecule that is considered as a potentially powerful greenhouse gas having the largest estimated lifetime of over 50 000 years in the atmosphere. The presence of relatively low vibrational frequencies makes the Boltzmann population of the excited levels important. Consequently, the "hot bands" corresponding to transitions among excited rovibrational states contribute significantly to the CF4 opacity in the infrared even at room temperature conditions but the existing laboratory data analyses are not sufficiently complete. In this work, we construct the first accurate and complete ab initio based line lists for CF4 in the range 0-4000 cm-1, containing rovibrational bands that are the most active in absorption. An efficient basis set compression method was applied to predict more than 700 new bands and subbands via variational nuclear motion calculations. We show that already at room temperature a quasi-continuum of overlapping weak lines appears in the CF4 infrared spectra due to the increasing density of bands and transitions. In order to converge the infrared opacity at room temperature, it was necessary to include a high rotational quantum number up to J = 80 resulting in 2 billion rovibrational transitions. In order to make the cross-section simulation faster, we have partitioned our data into two parts: (a) strong & medium line lists with lower energy levels for calculation of selective absorption features that can be used at various temperatures and (b) compressed "super-line" libraries of very weak transitions contributing to the quasi-continuum modelling. Comparisons with raw previously unassigned experimental spectra showed a very good accuracy for integrated absorbance in the entire range of the reported spectra predictions. The data obtained in this work will be made available through the TheoReTS information system (http://theorets.univ-reims.fr, http://theorets.tsu.ru) that contains ab initio born line lists and provides a user-friendly graphical interface for a fast simulation of the CF4 absorption cross-sections and radiance under various temperature conditions from 80 K to 400 K.

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, 18O + 32O2, 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.

Accurate ab initio dipole moment surfaces of ozone: First principle intensity predictions for rotationally resolved spectra in a large range of overtone and combination bands.

Ab initio dipole moment surfaces (DMSs) of the ozone molecule are computed using the MRCI-SD method with AVQZ, AV5Z, and VQZ-F12 basis sets on a dense grid of about 1950 geometrical configurations. The analytical DMS representation used for the fit of ab initio points provides better behavior for large nuclear displacements than that of previous studies. Various DMS models were derived and tested. Vibration-rotation line intensities of 16O3 were calculated from these ab initio surfaces by the variational method using two different potential functions determined in our previous works. For the first time, a very good agreement of first principle calculations with the experiment was obtained for the line-by-line intensities in rotationally resolved ozone spectra in a large far- and mid-infrared range. This includes high overtone and combination bands up to ΔV = 6. A particular challenge was a correct description of the B-type bands (even ΔV3 values) that represented major difficulties for the previous ab initio investigations and for the empirical spectroscopic models. The major patterns of various B-type bands were correctly described without empirically adjusted dipole moment parameters. For the 10 µm range, which is of key importance for the atmospheric ozone retrievals, our ab initio intensity results are within the experimental error margins. The theoretical values for the strongest lines of the ν3 band lie in general between two successive versions of HITRAN (HIgh-resolution molecular TRANsmission) empirical database that corresponded to most extended available sets of observations. The overall qualitative agreement in a large wavenumber range for rotationally resolved cold and hot ozone bands up to about 6000 cm-1 is achieved here for the first time. These calculations reveal that several weak bands are yet missing from available spectroscopic databases.

Ab initio variational predictions for understanding highly congested spectra: rovibrational assignment of 108 new methane sub-bands in the icosad range (6280-7800 cm(-1)).

A detailed study of methane spectra in the highly congested icosad range of 6280-7800 cm(-1) has been performed using global variational calculations derived from accurate ab initio potential energy and dipole moment surfaces. About 13,000 (12)CH4 lines of the WKLMC line lists recorded at 80 and 296 K using very sensitive laser techniques (DAS, CRDS) have been rovibrationally assigned from first principles predictions. Overall, a total of 7436 upper energy levels were determined. Among the 20 bands and the 134 sub-levels contained in the icosad system, 20 and 108 have been identified for the first time, respectively. The assigned transitions represent 98% of the sum of the experimental intensity at 80 K. This work demonstrates for the first time how accurate first principles global calculations allow assigning complicated spectra of a molecule with more than 4 atoms.

First predictions of rotationally resolved infrared spectra of dideuteromethane ((12)CH2D2) from potential energy and dipole moment surfaces.

We report the variationally computed infrared spectrum of (12)CH2D2 using our recent potential energy and dipole moment methane surfaces, which have been initially derived in the irreducible tensor representation adapted to the tetrahedral symmetry of the major isotopologue (12)CH4. The nuclear motion calculations are accomplished by combining the normal-mode Eckart-Watson Hamiltonian with isotopic and symmetry transformations. Our direct vibrational calculations are compared to the 93 observed band centers up to 6300 cm(-1). Except for two outliers the root-mean-square deviation is 0.22 cm(-1) and the maximum error is 0.7 cm(-1) without empirical adjustment of parameters. The work aims at filling the gap concerning missing line strength information for this molecule. Theoretical spectra predictions are given up to J = 25 and, for the very first time, ab initio intensity predictions for rovibrational line transitions are in good qualitative agreement with available experimental spectra.

A new accurate ground-state potential energy surface of ethylene and predictions for rotational and vibrational energy levels.

In this paper we report a new ground state potential energy surface for ethylene (ethene) C2H4 obtained from extended ab initio calculations. The coupled-cluster approach with the perturbative inclusion of the connected triple excitations CCSD(T) and correlation consistent polarized valence basis set cc-pVQZ was employed for computations of electronic ground state energies. The fit of the surface included 82,542 nuclear configurations using sixth order expansion in curvilinear symmetry-adapted coordinates involving 2236 parameters. A good convergence for variationally computed vibrational levels of the C2H4 molecule was obtained with a RMS(Obs.-Calc.) deviation of 2.7 cm(-1) for fundamental bands centers and 5.9 cm(-1) for vibrational bands up to 7800 cm(-1). Large scale vibrational and rotational calculations for (12)C2H4, (13)C2H4, and (12)C2D4 isotopologues were performed using this new surface. Energy levels for J = 20 up to 6000 cm(-1) are in a good agreement with observations. This represents a considerable improvement with respect to available global predictions of vibrational levels of (13)C2H4 and (12)C2D4 and rovibrational levels of (12)C2H4.