An isotopic-independent highly accurate potential energy surface for CO2 isotopologues and an initial 12C16O2 infrared line list

A two-step quadrature-based variational calculation of ro-vibrational levels and wavefunctions of CO2 using a bisector-x molecule-fixed frame

In this paper, we propose a new two-step strategy for computing ro-vibrational energy levels and wavefunctions of a triatomic molecule and apply it to CO2. A two-step method [J. Tennyson and B. T. Sutcliffe, Mol. Phys., 1986, 58, 1067] uses a basis whose functions are products of K-dependent "vibrational" functions and symmetric top functions. K is the quantum number for the molecule-fixed z component of the angular momentum. For a linear molecule, a two-step method is efficient because the Hamiltonian used to compute the basis functions includes the largest coupling term. The most important distinguishing feature of the two-step method we propose is that it uses an associated Legendre basis and quadrature rather than a K-dependent discrete variable representation. This reduces the cost of the calculation and simplifies the method. We have computed ro-vibrational energy levels with J up to 100 for CO2, on an accurate available potential energy surface which is known as the AMES-2 PES and present a subset of those levels. We have converged most levels up to 20 000 cm-1 to 0.0001 cm-1.

Methane throughout the atmosphere of the warm exoplanet WASP-80b

The abundances of main carbon- and oxygen-bearing gases in the atmospheres of giant exoplanets provide insights into atmospheric chemistry and planet formation processes1,2. Thermochemistry suggests that methane (CH4) should be the dominant carbon-bearing species below about 1,000 K over a range of plausible atmospheric compositions3; this is the case for the solar system planets4 and has been confirmed in the atmospheres of brown dwarfs and self-luminous, directly imaged exoplanets5. However, CH4 has not yet been definitively detected with space-based spectroscopy in the atmosphere of a transiting exoplanet6-11, but a few detections have been made with ground-based, high-resolution transit spectroscopy12,13 including a tentative detection for WASP-80b (ref. 14). Here we report transmission and emission spectra spanning 2.4-4.0 μm of the 825 K warm Jupiter WASP-80b taken with the NIRCam instrument of the JWST, both of which show strong evidence of CH4 at greater than 6σ significance. The derived CH4 abundances from both viewing geometries are consistent with each other and with solar to sub-solar C/O and around five times solar metallicity, which is consistent with theoretical predictions15-17.

A Spectroscopic Description of Asymmetric Isotopologues of CO2

A polyad-conserving algebraic model applied to vibrational excitations of asymmetric isotopologues of CO2 is presented. First, the problem of vibrational excitations is studied by taking into account only the minimum subspace of states to characterize the Fermi interaction. This analysis allows an estimation of the force constants as well as the feasibility of describing the system in a local mode scheme, in terms of SU(2) operators associated with Morse ladder operators for the stretches. This description together with the algebraic U(3) for the bends establishes the dynamical group SU1(2) × U(3) × SU2(2) for a series of isotopologues. Six isotopologues are considered, namely, 16O12C18O, 16O12C17O, 16O13C18O, 16O13C17O, 17O13C18O, and 17O12C18O in their electronic ground states. For isotopologues 16O12C18O, 16O12C17O, 17O12C18O, and 16O13C18O, the vibrational description was carried out using a Hamiltonian involving 14 parameters. For this series of isotopologues with a number of energy terms 90, 57, 42, and 40, the deviations obtained were rms = 0.15, 0.10, 0.06, and 0.07 cm-1, respectively. For 16O13C17O, with 28 experimental energies and involving 13 parameters, the deviation was rms = 0.05 cm-1, while for 17O13C18O, a different strategy was proposed since only 12 experimental energy levels. In all cases, the polyad scheme P212 = 2(ν1 + ν3) + ν2 was considered. In addition, a new criterion of locality/normality degree is proposed, embracing the case of molecules with normal mode behavior, in particular, the isotopologues of CO2.

Ames-2021 CO2 Dipole Moment Surface and IR Line Lists: Toward 0.1% Uncertainty for CO2 IR Intensities

A highly accurate CO₂ ab initio dipole moment surface (DMS), Ames-2021, is reported along with ¹²C¹⁶O₂ infrared (IR) intensity comparisons approaching a 1-4‰ level of agreement and uncertainty. The Ames-2021 DMS was accurately fitted from CCSD(T) finite-field dipoles computed with the aug-cc-pVXZ (X = T, Q, 5) basis for C atom and the d-aug-cc-pVXZ (X = T, Q, 5) basis for O atoms, and extrapolated to the one particle basis set limit. Fitting σ_rms is 3.8 × 10^-7 au for 4443 geometries below 15 000 cm^-1. The corresponding IR intensity, S_Ames-2021, are computed using the Ames-2 potential energy surface (PES), which is the best PES available for CO₂. Compared to high accuracy IR studies for 2001i-00001 and 3001i-00001 bands, S_Ames-2021 matches NIST experiment-based intensities [S_NIST-HIT16 or S_HIT20] to -1.0 ± 1.3‰, or matches DLR experiment-based intensities [S_{DLR-HIT16/UCL/Ames}] to 1.9 ± 3.7‰. This indicates the systematic deviations and uncertainties have been significantly reduced in S_Ames-2021. The S_UCL2015 (or S_HITRAN2016) have larger deviations (vs S_DLR) and uncertainties (vs S_DLR, S_NIST) which are attributed to the less accurate Ames-1 PES adopted in UCL-296 line list calculation. The S_Ames-2021 intensity of ¹²C¹⁶O₂ and ¹³C¹⁶O₂ is utilized to derive new absolute ¹³C/¹²C ratios for Vienna PeeDee Belemnite (VPDB) with uncertainty reduced by ¹/₃ or ²/₃. Further evaluation of S_Ames-2021 intensities are carried out on those CO₂ bands discussed in the HITRAN2020 update paper. Consistent improvements and better accuracies are found in band-by-band analysis, except for those bands strongly affected by Coriolis couplings, or very weak bands measured with relatively larger experimental uncertainties. The Ames-2021 296 K IR line lists are generated for 13 CO₂ isotopologues, with 18 000 cm^-1 and S_296 K > 1 × 10^-31 cm/molecule cutoff and then combined with CDSD line positions (except ¹⁴C¹⁶O₂). The Ames-2021 DMS and 296 K IR line lists represent a major improvement over previous CO₂ theoretical IR intensity studies, including Ames-2016, UCL-296, and recent UCL DMS 2021 update. A real 1 permille level of agreement and uncertainty will definitely require both more accurate PES and more accurate DMS.

What It Takes to Compute Highly Accurate Rovibrational Line Lists for Use in Astrochemistry

ConspectusWe review the Best Theory + Reliable High-Resolution Experiment (BTRHE) strategy for obtaining highly accurate molecular rovibrational line lists with InfraRed (IR) intensities. The need for highly accurate molecular rovibrational line lists is twofold: (a) assignment of the many rovibrational lines for common stable molecules especially those that exhibit a large amplitude motion, such as NH₃, or have a high density of states such as SO₂; and (b) characterization of the atmospheres of exoplanets, which will be one of the main areas of research in astronomy in the coming decades. The first motivation arises due to the need to eliminate lines due to common molecules in an astronomical observation in order to identify lines from new molecules, while the second motivation arises due to the need to obtain accurate molecular opacities in order to characterize the atmosphere of an exoplanet. The BTRHE strategy first consists of using high-quality ab initio quantum-chemical methods to obtain a global potential energy surface (PES) and dipole moment surface (DMS) that contains the proper physics. The global PES is then refined using a subset of the reliable high-resolution experimental data. The refined PES then gives energy-level predictions to an accuracy similar to the reproduction accuracy of the experimental data used in the refinement step in the interpolation region (i.e., within the range of the experimental data used in the refinement step). The accuracy of the energy levels will slowly degrade as they are extrapolated to spectral regions beyond the high-resolution experimental data used in the refinement step. However, because the degradation is slow, the predicted energy levels can be used to assign new high-resolution experiments, and the data from these can then be used in a subsequent refinement step. In this way, the global PES eventually can yield highly accurate energy levels for all desired spectral regions including to very high energies and high J values. We show that IR intensities computed with the BTRHE rovibrational wave functions and the DMS can be very accurate provided one has minimized the fitting error of the DMS and tested the completeness of the DMS. Some examples of our work on NH₃, CO₂, and SO₂ are given to highlight the usefulness of the BTRHE strategy and to provide ideas on how to further improve its predictive power in the future. In particular, it is shown how successive refinement steps, once new high-resolution data are available, can lead to PESs that yield highly accurate transition energies to larger spectral regions. The importance of including nonadiabatic corrections to reduce the J-dependence of errors for H-containing molecules is shown with work on NH₃. Another very important aspect of the BTRHE approach is the consistency across isotopologues, which allows for highly accurate line lists for any isotopologue once one is obtained for the main isotopologue (which has more high-resolution data available for refinement).

Heavy Particle Impact Vibrational Excitation and Dissociation Processes in CO2

A heavy particle impact vibrational excitation and dissociation model for CO₂ is presented. This state-to-state model is based on the forced harmonic oscillator (FHO) theory, which is more accurate than current state-of-the-art kinetic models of CO₂ based on first-order perturbation theory. The first excited triplet state ³B₂ of CO₂, including its vibrational structure, is considered in our model, and a more consistent approach to CO₂ dissociation is also proposed. The model is benchmarked against a few academic zero-dimensional (0D) cases and compared to decomposition time measurements in a shock tube. Our model is shown to have reasonable predictive capabilities, and the CO₂ + O ↔ CO + O₂ reaction is found to have a key influence on the dissociation dynamics of CO₂ shocked flows, warranting further theoretical studies. We conclude this study with a discussion on the theoretical improvements that are still required for a more consistent analysis of the vibrational/dissociation dynamics of CO₂.

Absolute 13C/12C Isotope Amount Ratio for Vienna Pee Dee Belemnite from Infrared Absorption Spectroscopy

Measurements of isotope ratios are predominantly made with reference to standard specimens that have been characterized in the past. In the 1950s, the carbon isotope ratio was referenced to a belemnite sample collected by Heinz Lowenstam and Harold Urey¹ in South Carolina's Pee Dee region. Due to the exhaustion of the sample since then, reference materials that are traceable to the original artefact are used to define the Vienna Pee Dee Belemnite (VPDB) scale for stable carbon isotope analysis². However, these reference materials have also become exhausted or proven to exhibit unstable composition over time³, mirroring issues with the international prototype of the kilogram that led to a revised International System of Units⁴. A campaign to elucidate the stable carbon isotope ratio of VPDB is underway⁵, but independent measurement techniques are required to support it. Here we report an accurate value for the stable carbon isotope ratio inferred from infrared absorption spectroscopy, fulfilling the promise of this fundamentally accurate approach⁶. Our results agree with a value recently derived from mass spectrometry⁵, and therefore advance the prospects of SI-traceable isotope analysis. Further, our calibration-free method could improve mass balance calculations and enhance isotopic tracer studies in CO₂ source apportionment.

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.

Comb-assisted, Pound-Drever-Hall locked cavity ring-down spectrometer for high-performance retrieval of transition parameters

Fast and high-performance cavity ring-down spectrometer (CRDS) is highly desired to precisely extract spectral parameters. In this paper, we present our comb-assisted Pound-Drever-Hall (PDH) locked CRDS setup, aiming to retrieve molecular parameters. In the setup, a dynamic feedback is used to keep the tight PDH locking even under strong absorption in the spectral measurement. PDH light and probing light enter the ring-down cavity simultaneously under orthogonal polarization, which enables a fast acquisition of ring-down events without interrupting PDH locking. Ultra-stable cavity temperature is realized, which has an accuracy below 0.5 mK in 27 minutes. The optical frequency comb (OFC) system is developed to rapidly and automatically measure the frequency axis with a relatively wide beat-note range. The minimum detectable absorption coefficient and noise-equivalent absorption coefficient (NEA) are 7.6×10^-12cm^-1 and 5.3×10^-12cm^-1Hz^-1/2, respectively. The spectrometer is implemented to measure CO₂ transition and extract line parameters. The uncertainty for line position is evaluated to be 120 kHz. An accuracy of 0.31% for line intensity is beneficial to the precise determination of CO₂ content for the purpose of environment protection and other applications.

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 3v₃ 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.

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)

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.

High-Accuracy CO(2) Line Intensities Determined from Theory and Experiment

Atmospheric CO(2) concentrations are being closely monitored by remote sensing experiments which rely on knowing line intensities with an uncertainty of 0.5% or better. Most available laboratory measurements have uncertainties much larger than this. We report a joint experimental and theoretical study providing rotation-vibration line intensities with the required accuracy. The ab initio calculations are extendible to all atmospherically important bands of CO(2) and to its isotologues. As such, they will form the basis for detailed CO(2) spectroscopic line lists for future studies.

Limited rotational and rovibrational line lists computed with highly accurate quartic force fields and ab initio dipole surfaces

In this work, computational procedures are employed to compute the rotational and rovibrational spectra and line lists for H2O, CO2, and SO2. Building on the established use of quartic force fields, MP2 and CCSD(T) Dipole Moment Surfaces (DMSs) are computed for each system of study in order to produce line intensities as well as the transition energies. The computed results exhibit a clear correlation to reference data available in the HITRAN database. Additionally, even though CCSD(T) DMSs produce more accurate intensities as compared to experiment, the use of MP2 DMSs results in reliable line lists that are still comparable to experiment. The use of the less computationally costly MP2 method is beneficial in the study of larger systems where use of CCSD(T) would be more costly.