Search

Light-matter interactions lie at the heart of our exploration of exoplanetary atmospheres. Interpreting data obtained by remote sensing is enabled by meticulous, time- and resource-consuming work aiming at deepening our understanding of such interactions (i.e., opacity models). Recently, \citet{Niraula2022} pointed out that due primarily to limitations on our modeling of broadening and far-wing behaviors, opacity models needed a timely update for exoplanet exploration in the JWST era, and thus argued for a scalable approach. In this Letter, we introduce an end-to-end solution from ab initio calculations to pressure broadening, and use the perturbation framework to identify the need for precision to a level of $\sim$10\%. We focus on the CO$_2$-H$_2$ system as CO$_2$ presents a key absorption feature for exoplanet research (primarily driven by the observation of gas giants) at $\sim$4.3$\mu$m and yet severely lack opacity data. We compute elastic and inelastic cross-sections for the collision of {ortho-}H$_2$ ~with CO$_2$, in the ground vibrational state, and at the coupled-channel fully converged level. For scattering energies above $\sim$20~cm$^{-1}$, moderate precision inter-molecular potentials are indistinguishable from high precision ones in cross-sections. Our calculations agree with the currently available measurement within 7\%, i.e., well beyond the precision requirements. Our proof-of-concept introduces a computationally affordable way to compute full-dimensional interaction potentials and scattering quantum dynamics with a precision sufficient to reduce the model-limited biases originating from the pressure broadening and thus support instrument-limited science with JWST and future missions., Comment: Submitted to ApJL

The goal of this white paper is to provide a snapshot of the data availability and data needs primarily for the Ariel space mission, but also for related atmospheric studies of exoplanets and brown dwarfs. It covers the following data-related topics: molecular and atomic line lists, line profiles, computed cross-sections and opacities, collision-induced absorption and other continuum data, optical properties of aerosols and surfaces, atmospheric chemistry, UV photodissociation and photoabsorption cross-sections, and standards in the description and format of such data. These data aspects are discussed by addressing the following questions for each topic, based on the experience of the "data-provider" and "data-user" communities: (1) what are the types and sources of currently available data, (2) what work is currently in progress, and (3) what are the current and anticipated data needs. We present a GitHub platform for Ariel-related data, with the goal to provide a go-to place for both data-users and data-providers, for the users to make requests for their data needs and for the data-providers to link to their available data. Our aim throughout the paper is to provide practical information on existing sources of data whether in databases, theoretical, or literature sources., Comment: 58 pages, submitted to RAS Techniques and Instruments (RASTI). The authors welcome feedback: corresponding author emails can be found as footnotes on page 2

The sulfur dimer (S$_2$) is an important molecular constituent in cometary atmospheres and volcanic plumes on Jupiter's moon Io. It is also expected to play an important role in the photochemistry of exoplanets. The UV spectrum of S$_2$ contains transitions between vibronic levels above and below the dissociation limit, giving rise to a distinctive spectral signature. By using spectroscopic information from the literature, and the spectral simulation program PGOPHER, a UV line list of S$_2$ is provided. This line list includes the primary $B\,^{3}\Sigma^{-}_{u}-X\,^{3}\Sigma^{-}_{g}$ ($v'$=0-27, $v''$=0-10) electronic transition, where vibrational bands with $v'$$\geq$10 are predissociated. Intensities have been calculated from existing experimental and theoretical oscillator strengths, and semi-empirical strengths for the predissociated bands of S$_2$ have been derived from comparisons with experimental cross-sections. The S$_2$ line list also includes the $B''\,^{3}\Pi_{u}-X\,^{3}\Sigma^{-}_{g}$ ($v'$=0-19, $v''$=0-10) vibronic bands due to the strong interaction with the $B$ state. In summary, we present the new HITRAN-formatted S$_2$ line list and its validation against existing laboratory spectra. The extensive line list covers the spectral range 21700$-$41300~cm$^{-1}$ ($\sim$242$-$461~nm) and can be used for modeling both absorption and emission., Comment: Accepted to the Monthly Notices of the Royal Astronomical Society

When computing cross-sections from a line list, the result depends not only on the line strength, but also the line shape, pressure-broadening parameters, and line wing cut-off (i.e., the maximum distance calculated from each line centre). Pressure-broadening can be described using the Lorentz lineshape, but it is known to not represent the true absorption in the far wings. Both theory and experiment have shown that far from the line centre, non-Lorentzian behaviour controls the shape of the wings and the Lorentz lineshape fails to accurately characterize the absorption, leading to an underestimation or overestimation of the opacity continuum depending on the molecular species involved. The line wing cut-off is an often overlooked parameter when calculating absorption cross sections, but can have a significant effect on the appearance of the spectrum since it dictates the extent of the line wing that contributes to the calculation either side of every line centre. Therefore, when used to analyse exoplanet and brown dwarf spectra, an inaccurate choice for the line wing cut-off can result in errors in the opacity continuum, which propagate into the modeled transit spectra, and ultimately impact/bias the interpretation of observational spectra, and the derived composition and thermal structure. Here, we examine the different methods commonly utilized to calculate the wing cut-off and propose a standard practice procedure (i.e., absolute value of 25~cm$^{-1}$ for $P\leqslant$~200~bar and 100~cm$^{-1}$ for $P >$ ~200~bar) to generate molecular opacities which will be used by the open-access {\tt MAESTRO} (Molecules and Atoms in Exoplanet Science: Tools and Resources for Opacities) database. The pressing need for new measurements and theoretical studies of the far-wings is highlighted., Comment: 13 pages, 6 figures

Advancements in space telescopes have opened new avenues for gathering vast amounts of data on exoplanet atmosphere spectra. However, accurately extracting chemical and physical properties from these spectra poses significant challenges due to the non-linear nature of the underlying physics. This paper presents novel machine learning models developed by the AstroAI team for the Ariel Data Challenge 2023, where one of the models secured the top position among 293 competitors. Leveraging Normalizing Flows, our models predict the posterior probability distribution of atmospheric parameters under different atmospheric assumptions. Moreover, we introduce an alternative model that exhibits higher performance potential than the winning model, despite scoring lower in the challenge. These findings highlight the need to reevaluate the evaluation metric and prompt further exploration of more efficient and accurate approaches for exoplanet atmosphere spectra analysis. Finally, we present recommendations to enhance the challenge and models, providing valuable insights for future applications on real observational data. These advancements pave the way for more effective and timely analysis of exoplanet atmospheric properties, advancing our understanding of these distant worlds., Comment: Conference proceeding for the ECML PKDD 2023

We use optical frequency comb Fourier transform spectroscopy to record high-resolution, low-pressure, room-temperature spectra of formaldehyde (H$_2$$^{12}$C$^{16}$O) in the range of 1250 to 1390 cm$^{-1}$. Through line-by-line fitting, we retrieve line positions and intensities of 747 rovibrational transitions: 558 from the ${\nu}_6$ band, 129 from the ${\nu}_4$ band, and 14 from the ${\nu}_3$ band, as well as 46 from four different hot bands. We incorporate the accurate and precise line positions (0.4 MHz median uncertainty) into the MARVEL (measured active vibration-rotation energy levels) analysis of the H$_2$CO spectrum. This increases the number of MARVEL-predicted energy levels by 82 and of rovibrational transitions by 5382, and substantially reduces uncertainties of MARVEL-derived H$_2$CO energy levels over a large range: from pure rotational levels below 200 cm$^{-1}$ up to multiply excited vibrational levels at 6000 cm$^{-1}$. This work is an important step toward filling the gaps in formaldehyde data in the HITRAN database.

As the James Webb Space Telescope (JWST) came online last summer, we entered a new era of astronomy. This new era is supported by data products of unprecedented information content that require novel reduction and analysis techniques. Recently, Niraula et al. 2022 (N22) highlighted the need for upgraded opacity models to prevent facing a model-driven accuracy wall when interpreting exoplanet transmission spectra. Here, we follow the same approach as N22 to explore the sensitivity of inferences on the atmospheric properties of WASP-39 b to the opacity models used. We find that the retrieval of the main atmospheric properties from this first JWST exoplanet spectrum is mostly unaffected by the current state of the community's opacity models. Abundances of strong opacity sources like water and carbon dioxide are reliably constrained within $\sim$0.30 dex, beyond the 0.50 dex accuracy wall reported in N22. Assuming the completeness and accuracy of line lists, N22's accuracy wall is primarily driven by model uncertainties on broadening coefficients and far-wing behaviors, which we find to have marginal consequences for interpreting the transmission spectra of large, hot, high-metallicity atmospheres such as WASP-39 b, in opposition to emission spectra and climate modeling which depend on deeper parts of a planetary atmosphere. The origin of the opacity challenge in the retrieval of metal-rich hot Jupiters via transmission spectroscopy will thus mostly stem from the incompleteness and inaccuracy of line lists., Comment: 10 Pages, 6 Figures

With a new generation of observatories coming online this decade, the process of characterizing exoplanet atmospheres will need to be reinvented. Currently mostly on the instrumental side, characterization bottlenecks will soon stand by the models used to translate spectra into atmospheric properties. Limitations stemming from our stellar and atmospheric models have already been highlighted. Here, we show that the current limitations of the opacity models used to decode exoplanet spectra propagate into an accuracy wall at ~0.5-1.0 dex (i.e., 3 to 10x) on the atmospheric properties, which is an order of magnitude above the precision targeted by JWST Cycle 1 programs and needed for, e.g., meaningful C/O-ratio constraints and biosignatures identification. We perform a sensitivity analysis using nine different opacity models and find that most of the retrievals produce harmonious fits owing to compensations in the form of >5$\sigma$ biases on the derived atmospheric parameters translating in the aforementioned accuracy wall. We suggest a two-tier approach to alleviate this problem involving a new retrieval procedure and guided improvements in opacity data, their standardization and optimal dissemination., Comment: 16 Pages, 12 Figures

The latest version of the HITRAN molecular spectroscopic database, HITRAN2020, has recently been released featuring many updates, including line-by-line broadening parameters (and their temperature dependence) appropriate for the dominant constituents of planetary atmospheres. In this work, line shape codes suitable for calculating microwave spectra have been implemented within the HITRAN Application Programming Interface (HAPI). These new additions allow for spectroscopic calculations of microwave absorbing species pertinent to current and future studies of the atmospheres of Jupiter and Venus, and more generally for the atmospheres of gas giants and rocky planets. The inversion spectrum of the NH$_3$ molecule broadened by H$_2$, He and H$_2$O dominates the microwave region of Jupiter. Whereas for Venus, accurate spectroscopic data of SO$_2$ broadened by CO$_2$ is necessary in order to determine its significance, if any, on the reported detection of PH$_3$ in the Venusian upper cloud deck. Comparisons have been made to available microwave laboratory opacities and the following results illustrate that HITRAN data can be used in conjunction with HAPI to reproduce the existing experimental measurements and provide reliable calculation of planetary opacities. Users should be mindful regarding selection of appropriate parameters in HITRAN and selecting suitable line shape functions in HAPI, depending on the spectral region, target molecular species, as well as ambient chemical and thermodynamic conditions., Comment: Accepted to the Monthly Notices of the Royal Astronomical Society (MNRAS) June 2022

We report the results of the trajectory-based simulation of far-infrared collision-induced absorption (CIA) due to CH$_4-$N$_2$ pairs at temperatures between 70 and 400 K. Our analysis utilizes recently calculated high-level potential energy (PES) and induced dipole surfaces (IDS) [Finenko, A. A., Chistikov, D. N., Kalugina, Y. N., Conway E. K., Gordon, I. E., Phys. Chem. Chem. Phys., 2021, doi: 10.1039/d1cp02161c]. Treating collision partners as rigid rotors, the time evolution of interaction-induced dipole is accumulated over a vast ensemble of classical trajectories and subsequently transformed into CIA spectrum via Fourier transform. In our calculations, both bound and unbound states are properly accounted for, and the rigorous theory of lower-order spectral moments is addressed to check the accuracy of simulated profiles. Classically derived trajectory-based profiles are subject to two approximate desymmetrization procedures so that resulting profiles conform to the quantum principle of detailed balance. The simulated profiles are compared to laboratory measurements and employed for modeling Titan's spectra in the 50-500 cm$^{-1}$ range. Based on the desymmetrized simulated profiles, a new semi-empirical model for CH$_4-$N$_2$ CIA is proposed for modeling Titan's infrared spectra. Synthetic spectra derived using this model yield an excellent agreement with the data recorded by the Composite Infrared Spectrometer (CIRS) aboard the Cassini spacecraft at low and high emission angles., Comment: 18 pages, 11 figures, accepted for publication in ApJS . The semiempirical CH4-N2 profiles tabulated in the HITRAN CIA format are available on the ApJS website or upon reasonable request

Molecular line lists, particularly those computed for high temperature applications, often have very few states assigned local quantum numbers. These are often important components for accurately determining line shape parameters required for radiative transfer simulations. The projection of the total angular momentum onto the molecule fixed axis ($k$) is investigated in the Radau internal coordinate system to determine when it can be considered a good quantum number. In such a coordinate system, when the square of the $k^{th}$ component of the wavefunction is greater than one half, then we can classify $k$ as a good quantum number in accordance with the theorem of Hose and Taylor. Furthermore, it is demonstrated that when this holds true, oblate and prolate quantum labels $K_{a}$ and $K_{c}$ can reliably be predicted. This is demonstrated for the water and ozone molecules.

The application described has been designed to create bibliographic entries in large databases with diverse sources automatically, which reduces both the frequency of mistakes and the workload for the administrators. This new system uniquely identifies each reference from its digital object identifier (DOI) and retrieves the corresponding bibliographic information from any of several online services, including the SAO/NASA Astrophysics Data Systems (ADS) and CrossRef APIs. Once parsed into a relational database, the software is able to produce bibliographies in any of several formats, including HTML and BibTeX, for use on websites or printed articles. The application is provided free-of-charge for general use by any scientific database. The power of this application is demonstrated when used to populate reference data for the HITRAN and AMBDAS databases as test cases. HITRAN contains data that is provided by researchers and collaborators throughout the spectroscopic community. These contributors are accredited for their contributions through the bibliography produced alongside the data returned by an online search in HITRAN. Prior to the work presented here, HITRAN and AMBDAS created these bibliographies manually, which is a tedious, time-consuming and error-prone process. The complete code for the new referencing system can be found at \url{https://github.com/hitranonline/refs}., Comment: 11 pages, 5 figures, already published online at https://doi.org/10.3390/atoms8020016

Theoretical calculations are contributing a significantly higher proportion of data to contemporary spectroscopic databases, which have traditionally relied on experimental observations and semi-empirical models. It is now a common procedure to extend calculated line lists to include ro-vibrational transitions between all bound states of the ground electronic state up to the dissociation limit. Advanced ab initio methods are utilized to calculate the potential energy and dipole moment surfaces (PESs and DMSs), and semi-empirical PESs are then obtained by combining ab initio and experimental data. The objective is to reach high accuracy in the calculated transition intensities for all parts of spectrum, i.e. to increase the predictive power of the model. We show that in order to perform this task, one needs, in addition to the standard improvements of the PES and DMS in the spectroscopically accessible regions, to extend the ab initio calculations of the PES towards the united-atom limit along the stretching coordinates. The argument is based on the correlation between the intensities of high-overtone transitions and the repulsive potential wall that has previously been theoretically established for diatomic molecules and is empirically extended here to linear and nonlinear triatomic molecules. We generate partial line lists for water and ozone, and together with an already available line list for carbon dioxide, we derive the normal intensity distribution, which is a direct consequence of this correlation. The normal distribution is not an instrument to compute highly accurate intensities, rather it is a means to analyze the intensities computed by the traditional methods.

We present an evaluation of airborne Intensity-Modulated Continuous-Wave (IM-CW) lidar measurements of atmospheric column CO2 mole fractions during the ACT-America project. This lidar system transmits online and offline wavelengths simultaneously on the 1.57111-um CO2 absorption line, with each modulated wavelength using orthogonal swept frequency waveforms. After the spectral characteristics of this system were calibrated through short-path measurements, we used the HITRAN spectroscopic database to derive the average-column CO2 mixing ratio (XCO2) from the lidar measured optical depths. Based on in situ measurements of meteorological parameters and CO2 concentrations for calibration data, we demonstrate that our lidar CO2 measurements were consistent from season to season and had an absolute calibration error (standard deviation) of 0.80 ppm when compared to XCO2 values derived from in situ measurements. By using a 10-second or longer moving average, a long-term stability of 1 ppm or better was obtained. The estimated CO2 measurement precision for 0.1-s, 1-s, 10-s, and 60-s averages were determined to be 3.4 ppm (0.84%), 1.2 ppm (0.30%), 0.43 ppm (0.10%), and 0.26 ppm (0.063%), respectively. These correspond to measurement signal-to-noise ratios of 120, 330, 950, and 1600, respectively. The drift in XCO2 over one-hour of flight time was found to be below our detection limit of about 0.1 ppm. These analyses demonstrate that the measurement stability, precision and accuracy are all well below the thresholds needed to study synoptic-scale variations in atmospheric XCO2., Comment: 20 pages, 5 figures. Submitted to Earth, Space, and Science (AGU)

New line lists are presented for the two most abundant water isotopologues; H$_{2}$$^{16}$O and H$_{2}$$^{18}$O. The H$_{2}$$^{16}$O line list extends to 25710 cm$^{-1}$ with intensity stabilities provided via ratios of calculated intensities obtained from two different semi-empirical potential energy surfaces. The line list for H$_{2}$$^{18}$O extends to 20000 cm$^{-1}$. The minimum intensity considered for all is $10^{-30}$ cm molecule$^{-1}$ at 296~K, assuming 100\% abundance for each isotopologue. Fluctuation of calculated intensities caused by changes in the underlying potential energy are found to be significant, particularly for weak transitions. Direct comparisons are made against eighteen different sources of line intensities, both experimental and theoretical, many of which are used within the HITRAN2016 database. With some exceptions, there is excellent agreement between our line lists and the experimental intensities in HITRAN2016. In the infrared region, many H$_{2}$$^{16}$O bands which exhibit intensity differences of 5-10\% between to the most recent 'POKAZATEL' line list (Polyansky \textit{et al.}, [Mon. Not. Roy. Astron. Soc. \textbf{480}, 2597 (2018)] and observation, are now generally predicted to within 1\%. For H$_{2}$$^{18}$O, there are systematic differences in the strongest intensities calculated in this work versus those obtained from semi-empirical calculations. In the visible, computed cross sections show smaller residuals between our work and both HITRAN2016 and HITEMP2010 than POKAZATEL. While our line list accurately reproduces HITEMP2010 cross sections in the observed region, residuals produced from this comparison do however highlight the need to update line positions in the visible spectrum of HITEMP2010. These line lists will be used to update many transition intensities and line positions in the HITRAN2016 database.

A methane line list for the HITEMP spectroscopic database, covering 0-13,400 cm$^{-1}$ ($>$746 nm), is presented. To create this compilation, ab initio line lists of $^{12}$CH$_{4}$ from Rey et al. (2017) ApJ, 847, 105 (provided at separate temperatures in the TheoReTS information system), are now combined with HITRAN2016 methane data to produce a single line list suitable for high-temperature line-by-line calculations up to 2000 K. An effective-temperature interpolation model was created in order to represent continuum-like features at any temperature of interest. This model is advantageous to previously-used approaches that employ so-called ``super-lines'', which are suitable only at a given temperature and require separate line lists for different temperatures. The resultant HITEMP line list contains $\sim$32 million lines and is significantly more flexible than alternative line lists of methane, while accuracy required for astrophysical or combustion applications is retained. Comparisons against experimental observations of methane absorption at high temperatures have been used to demonstrate the accuracy of the new work. The line list includes both strong lines and quasi-continuum features and is provided in the common user-friendly HITRAN/HITEMP format, making it the most practical methane line list for radiative transfer modeling at high-temperature conditions., Comment: 25 pages, 12 figures, submitted to the Astrophysical Journal Supplement Series

Calculating dipole moments with high-order basis sets is generally only possible for the light molecules, such as water. A simple, yet highly effective strategy of obtaining high-order dipoles with small, computationally less expensive basis sets is described. Using the finite field method for computing dipoles, energies calculated with small basis sets can be extrapolated to produce dipoles that are comparable to those obtained in high order calculations. The method reduces computational resources by approximately 50% (allowing the calculation of reliable dipole moments for larger molecules) and simultaneously improves the agreement with experimentally measured infrared transition intensities. For atmospherically important molecules which are typically too large to consider the use of large basis sets, this procedure will provide the necessary means of improving calculated spectral intensities by several percent.

The amount of water vapor in the terrestrial atmosphere is highly variable both spatially and temporally. In the tropics it sometimes constitutes 4-5% of the atmosphere. At the same time collisional broadening of spectral lines by water vapor is much larger than that by nitrogen and oxygen. Therefore, in order to accurately characterize and model spectra of the atmospheres with significant amounts of water vapor, the line-shape parameters for spectral lines broadened by water vapor are required. In this work, the line-broadening coefficients (and their temperature dependence exponents) due to the pressure of water vapor for lines of CO2, N2O, CO, CH4, O2, NH3, and H2S from both experimental and theoretical studies were collected and carefully reviewed. A set of semi-empirical models based on these collected data was created and then used to estimate water broadening and its temperature dependence for all transitions of selected molecules in the HITRAN2016 database., Comment: Submitted to the Journal of Geophysical Research: Atmospheres on May 23, 2019. 3 Tables, 11 Figures

We are now on a clear trajectory for improvements in exoplanet observations that will revolutionize our ability to characterize their atmospheric structure, composition, and circulation, from gas giants to rocky planets. However, exoplanet atmospheric models capable of interpreting the upcoming observations are often limited by insufficiencies in the laboratory and theoretical data that serve as critical inputs to atmospheric physical and chemical tools. Here we provide an up-to-date and condensed description of areas where laboratory and/or ab initio investigations could fill critical gaps in our ability to model exoplanet atmospheric opacities, clouds, and chemistry, building off a larger 2016 white paper, and endorsed by the NAS Exoplanet Science Strategy report. Now is the ideal time for progress in these areas, but this progress requires better access to, understanding of, and training in the production of spectroscopic data as well as a better insight into chemical reaction kinetics both thermal and radiation-induced at a broad range of temperatures. Given that most published efforts have emphasized relatively Earth-like conditions, we can expect significant and enlightening discoveries as emphasis moves to the exotic atmospheres of exoplanets., Comment: Submitted as an Astro2020 Science White Paper

This work describes the update of NO along with the incorporation of NO$_{2}$ and N$_{2}$O to the HITEMP database. Where appropriate, the HITRAN line lists for the same molecules have also been updated. This work brings the current number of molecules provided by HITEMP to seven. The initial line lists originating from \textit{ab initio} and semi-empirical methods for each molecule have been carefully validated against available observations and, where necessary, adjustments have been made to match observations. We anticipate this work will be applied to a variety of high-temperature environments including astronomical applications, combustion monitoring, and non-local thermodynamic equilibrium conditions.

A new line list including positions and absolute intensities (in the form of Einstein $A$ values and oscillator strengths) has been produced for the OH ground X\DP\ state rovibrational (Meinel system) and pure rotational transitions. All possible transitions are included with v$\primed$ and v$\Dprimed$ up to 13, and $J$ up to between 9.5 and 59.5, depending on the band. An updated fit to determine molecular constants has been performed, which includes some new rotational data and a simultaneous fitting of all molecular constants. The absolute line intensities are based on a new dipole moment function, which is a combination of two high level ab initio calculations. The calculations show good agreement with an experimental v=1 lifetime, experimental $\mu_\mathrm{v}$ values, and $\Delta$v=2 line intensity ratios from an observed spectrum. To achieve this good agreement, an alteration in the method of converting matrix elements from Hund's case (b) to (a) was made. Partitions sums have been calculated using the new energy levels, for the temperature range 5-6000 K, which extends the previously available (in HITRAN) 70-3000 K range. The resulting absolute intensities have been used to calculate O abundances in the Sun, Arcturus, and two red giants in the Galactic open and globular clusters M67 and M71. Literature data based mainly on [O I] lines are available for the Sun and Arcturus, and excellent agreement is found., Comment: 17 pages, 8 figues. 7 supplementary files: dipole moment functions (OH-X-DMFs.txt), equilibrium constants (OH-X-Equilibrium_Constants.txt), partition function (OH-X-Q_5-6000K.dat), PGOPHER file with molecular constants and transition matric elements (OH-XX.pgo), vibrational Einstein A and f values (OH-XX-Avv_fvv.txt), line list (OH-XX-Line_list.txt), and OH-Transformation_Equation_Extra.docx