A hybrid line list for CH4 and hot methane continuum

Sulfur dioxide in the mid-infrared transmission spectrum of WASP-39b

The recent inference of sulfur dioxide (SO2) in the atmosphere of the hot (approximately 1,100 K), Saturn-mass exoplanet WASP-39b from near-infrared JWST observations1-3 suggests that photochemistry is a key process in high-temperature exoplanet atmospheres4. This is because of the low (<1 ppb) abundance of SO2 under thermochemical equilibrium compared with that produced from the photochemistry of H2O and H2S (1-10 ppm)4-9. However, the SO2 inference was made from a single, small molecular feature in the transmission spectrum of WASP-39b at 4.05 μm and, therefore, the detection of other SO2 absorption bands at different wavelengths is needed to better constrain the SO2 abundance. Here we report the detection of SO2 spectral features at 7.7 and 8.5 μm in the 5-12-μm transmission spectrum of WASP-39b measured by the JWST Mid-Infrared Instrument (MIRI) Low Resolution Spectrometer (LRS)10. Our observations suggest an abundance of SO2 of 0.5-25 ppm (1σ range), consistent with previous findings4. As well as SO2, we find broad water-vapour absorption features, as well as an unexplained decrease in the transit depth at wavelengths longer than 10 μm. Fitting the spectrum with a grid of atmospheric forward models, we derive an atmospheric heavy-element content (metallicity) for WASP-39b of approximately 7.1-8.0 times solar and demonstrate that photochemistry shapes the spectra of WASP-39b across a broad wavelength range.

Sub-Doppler optical-optical double-resonance spectroscopy using a cavity-enhanced frequency comb probe

Accurate parameters of molecular hot-band transitions, i.e., those starting from vibrationally excited levels, are needed to accurately model high-temperature spectra in astrophysics and combustion, yet laboratory spectra measured at high temperatures are often unresolved and difficult to assign. Optical-optical double-resonance (OODR) spectroscopy allows the measurement and assignment of individual hot-band transitions from selectively pumped energy levels without the need to heat the sample. However, previous demonstrations lacked either sufficient resolution, spectral coverage, absorption sensitivity, or frequency accuracy. Here we demonstrate OODR spectroscopy using a cavity-enhanced frequency comb probe that combines all these advantages. We detect and assign sub-Doppler transitions in the spectral range of the 3ν3 ← ν3 resonance of methane with frequency precision and sensitivity more than an order of magnitude better than before. This technique will provide high-accuracy data about excited states of a wide range of molecules that is urgently needed for theoretical modeling of high-temperature data and cannot be obtained using other methods.

Early Release Science of the exoplanet WASP-39b with JWST NIRISS

The Saturn-mass exoplanet WASP-39b has been the subject of extensive efforts to determine its atmospheric properties using transmission spectroscopy1-4. However, these efforts have been hampered by modelling degeneracies between composition and cloud properties that are caused by limited data quality5-9. Here we present the transmission spectrum of WASP-39b obtained using the Single-Object Slitless Spectroscopy (SOSS) mode of the Near Infrared Imager and Slitless Spectrograph (NIRISS) instrument on the JWST. This spectrum spans 0.6-2.8 μm in wavelength and shows several water-absorption bands, the potassium resonance doublet and signatures of clouds. The precision and broad wavelength coverage of NIRISS/SOSS allows us to break model degeneracies between cloud properties and the atmospheric composition of WASP-39b, favouring a heavy-element enhancement ('metallicity') of about 10-30 times the solar value, a sub-solar carbon-to-oxygen (C/O) ratio and a solar-to-super-solar potassium-to-oxygen (K/O) ratio. The observations are also best explained by wavelength-dependent, non-grey clouds with inhomogeneous coverageof the planet's terminator.

Early Release Science of the exoplanet WASP-39b with JWST NIRSpec G395H

Measuring the abundances of carbon and oxygen in exoplanet atmospheres is considered a crucial avenue for unlocking the formation and evolution of exoplanetary systems1,2. Access to the chemical inventory of an exoplanet requires high-precision observations, often inferred from individual molecular detections with low-resolution space-based3-5 and high-resolution ground-based6-8 facilities. Here we report the medium-resolution (R ≈ 600) transmission spectrum of an exoplanet atmosphere between 3 and 5 μm covering several absorption features for the Saturn-mass exoplanet WASP-39b (ref. 9), obtained with the Near Infrared Spectrograph (NIRSpec) G395H grating of JWST. Our observations achieve 1.46 times photon precision, providing an average transit depth uncertainty of 221 ppm per spectroscopic bin, and present minimal impacts from systematic effects. We detect significant absorption from CO2 (28.5σ) and H2O (21.5σ), and identify SO2 as the source of absorption at 4.1 μm (4.8σ). Best-fit atmospheric models range between 3 and 10 times solar metallicity, with sub-solar to solar C/O ratios. These results, including the detection of SO2, underscore the importance of characterizing the chemistry in exoplanet atmospheres and showcase NIRSpec G395H as an excellent mode for time-series observations over this critical wavelength range10.

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.

Efficient and Automated Quantum Chemical Calculation of Rovibrational Nonresonant Raman Spectra

An outline of a newly developed program for the simulation of rovibrational nonresonant Raman spectra is presented. This program is an extension of our recently developed code for rovibrational infrared spectra [J. Chem Phys. 152 (2020) 244104] and relies on vibrational wavefunctions from variational configuration interaction theory to allow for an almost fully automated calculation of such spectra in pure ab initio fashion. Due to efficient contraction schemes this program requires modest computational resources and it can be controlled by only a few lines of input. As the required polarizability surfaces are also computed in an automated fashion, this implementation enables the routine application to small molecules. For demonstrating its capabilities, benchmark calculations for water H216O are compared to reference data and spectra for the beryllium dihydride dimer, Be2H4 (D2h), are predicted. The inversion symmetry of the D2h systems lead to complementary infrared and Raman spectra, which are needed both for a comprehensive investigation of this system.

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.

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.

High enthalpy source dedicated to quantitative infrared emission spectroscopy of gas flows at elevated temperatures

The High Enthalpy Source (HES) is a novel high temperature source developed to measure infrared line-by-line integrated absorption cross sections of flowing gases up to 2000 K. The HES relies on a porous graphite furnace designed to uniformly heat a constant flow of gas. The flow compensates thermal dissociation by renewing continuously the gas sample and eliminating dissociation products. The flowing characteristics have been investigated using computational fluid dynamics simulation confirming good temperature uniformity. The HES has been coupled to a high-resolution Fourier transform spectrometer to record emission spectra of methane at temperatures ranging between 700 and 1400 K. A radiative model has been developed to extract absolute line intensities from the recorded spectra.

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.