Search

Multi-reference configuration interaction (MRCI) potential energy curves (PECs) and spin-orbit couplings for the X $^2\Pi$, A $^2 \Sigma^+$, 1 $^2 \Sigma^-$, 1 $^4 \Sigma^-$, and 1 $^4 \Pi$ states of OH are computed and refined against empirical energy levels and transitions to produce a spectroscopic model. Predissociation lifetimes are determined by discretising continuum states in the variational method nuclear motion calculation by restricting the calculation to finite range of internuclear separations. Varying this range give a series of avoided crossings between quasi-bound states associated with the A $^2 \Sigma^+$ and continuum states, from which predissociation lifetimes are extracted. 424 quasi-bound A $^2 \Sigma^+$ state rovibronic energy levels are analysed and 374 predissociation lifetimes are produced, offering good coverage of the predissociation region. Agreement with measured lifetimes is satisfactory and a majority of computed results were within experimental uncertainty. A previously unreported A $^2 \Sigma^+$ state predissociation channel which goes via the X $^2\Pi$ is identified in the calculations. A python package, binslt, is produced to calculate predissociation lifetimes, associated line broadening parameters, and uncertainties from Duo *.states files is made available. The PECs and other curves from this work will be used to produce a rovibronic ExoMol linelist and temperature-dependent photodissociation cross sections for the hydroxyl radical., Comment: 14 Pages, 11 figures, 11 tables, 3 supplementary files. To be published in Journal of Chemical Physics

We present an ab initio study of the rovibronic spectra of sulfur monoxide ($^{32}$S$^{16}$O) using internally contracted multireference confoguration interaction (ic-MRCI) method and aug-cc-pV5Z basis sets. It covers 13 electronic states $X^{3}\Sigma^{-}$, $a^{1}\Delta$, $b^{1}\Sigma^{+}$, $c^{1}\Sigma^{-}$, $A^{\prime\prime 3}\Sigma^{+}$, $A^{\prime 3}\Delta$, $A^{3}\Pi$, $B^{3}\Sigma^{-}$, $C^{3}\Pi$, $d^{1}\Pi$, $e^{1}\Pi$, $C^{\prime 3}\Pi$, and $(3)^{1}\Pi$ ranging up to 66800 cm$^{-1}$. The ab initio spectroscopic model includes 13 potential energy curves, 23 dipole and transition dipole moment curves, 23 spin-orbit curves, and 14 electronic angular momentum curves. A diabatic representation is built by removing the avoided crossings between the spatially degenerate pairs $C^{3}\Pi - C^{\prime 3}\Pi$ and $e^{1}\Pi - (3)^{1}\Pi$ through a property-based diabatisation method. We also present non-adiabatic couplings and diabatic couplings for these avoided crossing systems. All phases for our coupling curves are defined, and consistent, providing the first fully reproducible spectroscopic model of SO covering the wavelength range longer than 147 nm. Finally, an ab initio rovibronic spectrum of SO is computed., Comment: 15 pages, 14 figures

We present a unified variational treatment of the electric quadrupole (E2) matrix elements, Einstein coefficients, and line strengths for general open-shell diatomic molecules in the general purpose diatomic code \Duo. Transformation relations between the Cartesian representation (typically used in electronic structure calculations) to the tensorial representation (required for spectroscopic applications) of the electric quadrupole moment components are derived. The implementation has been validated against accurate theoretical calculations and experimental measurements of quadrupole intensities of $^1\text{H}_2$ available in the literature. We also present accurate electronic structure calculations of the electric quadrupole moment functions for the $X^1\Sigma^+$ electronic states of $\text{CO}$ and $\text{HF}$ at the CCSD(T) and MRCI levels of theory, respectively, as well for the $a^1\Delta_g$ -- $b^1\Sigma_g^+$ quadrupole transition moment of $\text{O}_2$ with MRCI level of theory. Accurate infrared E2 line lists for $^{12}\text{C}^{16}\text{O}$ and $^1\text{H}^{19}\text{F}$ are provided. A demonstration of spectroscopic applications is presented by simulating E2 spectra for $^{12}\text{C}^{16}\text{O}$, $\text{H}^{19}\text{F}$ and $^{16}\text{O}_2$ (Noxon $a^1\Delta_g$ -- $b^1\Sigma_g^+$ band)., Comment: 38 pages, 12 figures

Computing and using opacities is a key part of modeling and interpreting data of exoplanetary atmospheres. Since the underlying spectroscopic line lists are constantly expanding and currently include up to ~ 10^10 - 10^11 transition lines, the opacity calculator codes need to become more powerful. Here we present major upgrades to the HELIOS-K GPU-accelerated opacity calculator and describe the necessary steps to process large line lists within a reasonable amount of time. Besides performance improvements, we include more capabilities and present a toolbox for handling different atomic and molecular data sets: from downloading and pre-processing the data to performing the opacity calculations in a user-friendly way. HELIOS-K supports line lists from ExoMol, HITRAN, HITEMP, NIST, Kurucz and VALD3. By matching the resolution of 0.1 cm^-1 and cutting length of 25 cm^-1 used by the ExoCross code for timing performance (251 seconds excluding data read-in time), HELIOS-K can process the ExoMol BT2 water line list in 12.5 seconds. Using a resolution of 0.01 cm^-1, it takes 45 seconds - equivalent to about 10^7 lines per second. As a wavenumber resolution of 0.01 cm^-1 suffices for most exoplanetary atmosphere spectroscopic calculations, we adopt this resolution in calculating opacity functions for several hundred atomic and molecular species, and make them freely available on the open-access DACE database. For the opacity calculations of the database, we use a cutting length of 100 cm^-1 for molecules and no cutting length for atoms. Our opacities are available for downloading from https://dace.unige.ch/opacityDatabase and may be visualized using https://dace.unige.ch/opacity., Comment: Published in The Astrophysical Journal Supplement Series

With the upcoming launch of space telescopes dedicated to the study of exoplanets, the \textit{Atmospheric Remote-Sensing Infrared Exoplanet Large-survey} (ARIEL) and the \textit{James Webb Space Telescope} (JWST), a new era is opening in the exoplanetary atmospheric explorations. However, especially in relatively cold planets around later-type stars, photochemical hazes and clouds may mask the composition of the lower part of the atmosphere, making it difficult to detect any chemical species in the troposphere or to understand whether there is a surface or not. This issue is particularly exacerbated if the goal is to study the habitability of said exoplanets and to search for biosignatures.\par This work combines innovative laboratory experiments, chemical modeling and simulated observations at ARIEL and JWST resolutions. We focus on the signatures of molecular ions that can be found in upper atmospheres above cloud decks. Our results suggest that H$_3^+$ along with H$_3$O$^+$ could be detected in the observational spectra of sub-Neptunes based on realistic mixing ratio assumption. This new parametric set may help to distinguish super-Earths with a thin atmosphere from H$_2$-dominated sub-Neptunes, to address the critical question whether a low-gravity planet around a low-mass active star is able to retain its volatile components. These ions may also constitute potential tracers to certain molecules of interest such as H$_2$O or O$_2$ to probe the habitability of exoplanets. Their detection will be an enthralling challenge for the future JWST and ARIEL telescopes.

High resolution spectroscopy (R > 20,000) is currently the only known method to constrain the orbital solution and atmospheric properties of non-transiting hot Jupiters. It does so by resolving the spectral features of the planet into a forest of spectral lines and directly observing its Doppler shift while orbiting the host star. In this study, we analyse VLT/CRIRES (R = 100,000) L-band observations of the non-transiting giant planet HD 179949 b centred around 3.5 microns. We observe a weak (3.0 sigma, or S/N = 4.8) spectral signature of H2O in absorption contained within the radial velocity of the planet at superior-conjunction, with a mild dependence on the choice of line list used for the modelling. Combining this data with previous observations in the K-band, we measure a detection significance of 8.4 sigma for an atmosphere that is most consistent with a shallow lapse-rate, solar C/O ratio, and with CO and H2O being the only major sources of opacity in this wavelength range. As the two sets of data were taken three years apart, this points to the absence of strong radial-velocity anomalies due, e.g., to variability in atmospheric circulation. We measure a projected orbital velocity for the planet of KP = (145.2 +- 2.0)kms^{-1} (1 sigma) and improve the error bars on this parameter by ~70%. However, we only marginally tighten constraints on orbital inclination (66.2 +3.7 -3.1 degrees) and planet mass (0.963 +0.036 -0.031 Jupiter masses), due to the dominant uncertainties of stellar mass and semi-major axis. Follow ups of radial-velocity planets are thus crucial to fully enable their accurate characterisation via high resolution spectroscopy., Comment: Accepted for publication in MNRAS; 13 pages, 10 figures

A new ro-vibrational line list for the ground electronic state of the main isotopologue of acetylene, $^{12}$C$_2$H$_2$, is computed as part of the ExoMol project. The aCeTY line list covers the transition wavenumbers up to 10,000 cm$^{-1}$ ($ \lambda >1$ $\mu$m), with lower and upper energy levels up to 12,000 cm$^{-1}$ and 22,000 cm$^{-1}$ considered, respectively. The calculations are performed up to a maximum value for the vibrational angular momentum, $K_{\rm max}=L_{\rm max}$ = 16, and maximum rotational angular momentum, $J$ = 99. Higher values of $J$ were not within the specified wavenumber window. The aCeTY line list is considered to be complete up to 2200 K, making it suitable for use in characterising high-temperature exoplanet or cool stellar atmospheres. Einstein-A coefficients, which can directly be used to calculate intensities at a particular temperature, are computed for 4.3 billion (4,347,381,911) transitions between 5 million (5,160,803) energy levels. We make comparisons against other available data for $^{12}$C$_2$H$_2$, and demonstrate this to be the most complete line list available. The line list is available in electronic form from the online CDS and ExoMol databases.

In the past decade, observations from space and ground have found H$_2$O to be the most abundant molecular species, after hydrogen, in the atmospheres of hot, gaseous, extrasolar planets. Being the main molecular carrier of oxygen, H$_2$O is a tracer of the origin and the evolution mechanisms of planets. For temperate, terrestrial planets, the presence of H$_2$O is of great significance as an indicator of habitable conditions. Being small and relatively cold, these planets and their atmospheres are the most challenging to observe, and therefore no atmospheric spectral signatures have so far been detected. Super-Earths -- planets lighter than ten M$_\oplus$ -- around later-type stars may provide our first opportunity to study spectroscopically the characteristics of such planets, as they are best suited for transit observations. Here we report the detection of an H$_2$O spectroscopic signature in the atmosphere of \planet\ -- an eight M$_\oplus$ planet in the habitable-zone of an M-dwarf -- with high statistical confidence (ADI = 5.0, $\sim$3.6$\sigma$). In addition, the derived mean molecular weight suggests an atmosphere still containing some hydrogen. The observations were recorded with the Hubble Space Telescope/WFC3 camera, and analysed with our dedicated, publicly available, algorithms. While the suitability of M-dwarfs to host habitable worlds is still under discussion, \planet\ offers an unprecedented opportunity to get insight into the composition and climate of habitable-zone planets., Comment: Published in Nature Astronomy

In the present work, we report a detailed description of the symmetry properties of the eight-atomic molecule ethane, with the aim of facilitating the variational calculations of rotation-vibration spectra of ethane and related molecules. Ethane consists of two methyl groups $\text{CH}_3$ where the internal rotation (torsion) of one $\text{CH}_3$ group relative to the other is of large amplitude and involves tunneling between multiple minima of the potential energy function. The molecular symmetry group of ethane is the 36-element group $G_{36}$ but the construction of symmetrized basis functions is most conveniently done in terms of the 72-element extended molecular symmetry group $G_{36}\text{(EM)}$. This group can subsequently be used in the construction of block-diagonal matrix representations of the ro-vibrational Hamiltonian for ethane. The derived transformation matrices associated with $G_{36}\text{(EM)}$ have been implemented in the variational nuclear motion program TROVE (Theoretical ROVibrational Energies). TROVE variational calculations will be used as a practical example of a $G_{36}\text{(EM)}$ symmetry adaptation for large systems with a non-rigid, torsional degree of freedom. We present the derivation of irreducible transformation matrices for all 36 (72) operations of $G_{36}\text{(M)}$ ($G_{36}\text{(EM)}$) and also describe algorithms for a numerical construction of these matrices based on a set of four (five) generators. The methodology presented is illustrated on the construction of the symmetry-adapted representations both of the potential energy function of ethane and of the rotation, torsion and vibration basis set functions., Comment: Main body: 23 pages, 11 figures, 5 tables

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.

An analysis of currently available ammonia (NH$_3$) visible-to-near-infrared gas absorption data was recently undertaken by Irwin et al. (Icarus, 302 (2018) 426) to help interpret Very Large Telescope (VLT) MUSE observations of Jupiter from 0.48 - 0.93 $\mu$m, made in support of the NASA/Juno mission. Since this analysis a newly revised set of ammonia line data, covering the previously poorly constrained range 0.5 - 0.833 $\mu$m, has been released by the ExoMol project, "C2018" (Coles et al., JQSRT 219, 199 - 122, 2018), which demonstrates significant advantages over previously available data sets, and providing for the first time complete line data for the previously poorly constrained 5520- and 6475-\AA\ bands of NH$_3$. In this paper we compare spectra calculated using the ExoMol-C2018 data set (Coles et al., JQSRT 219, 199 - 122, 2018) with spectra calculated from previous sources to demonstrate its advantages. We conclude that at the present time the ExoMol-C2018 dataset provides the most reliable ammonia absorption source for analysing low- to medium-resolution spectra of Jupiter in the visible/near-IR spectral range, but note that the data are less able to model high-resolution spectra owing to small, but significant inaccuracies in the line wavenumber estimates. This work is of significance not only for solar system planetary physics, but for future proposed observations of Jupiter-like planets orbiting other stars, such as with NASA's planned Wide-Field Infrared Survey Telescope (WFIRST)., Comment: 12 Figures

A new nine-dimensional potential energy surface (PES) for methane has been generated using state-of-the-art \textit{ab initio} theory. The PES is based on explicitly correlated coupled cluster calculations with extrapolation to the complete basis set limit and incorporates a range of higher-level additive energy corrections. These include: core-valence electron correlation, higher-order coupled cluster terms beyond perturbative triples, scalar relativistic effects and the diagonal Born-Oppenheimer correction. Sub-wavenumber accuracy is achieved for the majority of experimentally known vibrational energy levels with the four fundamentals of $^{12}$CH$_4$ reproduced with a root-mean-square error of $0.70{\,}$cm$^{-1}$. The computed \textit{ab initio} equilibrium C{--}H bond length is in excellent agreement with previous values despite pure rotational energies displaying minor systematic errors as $J$ (rotational excitation) increases. It is shown that these errors can be significantly reduced by adjusting the equilibrium geometry. The PES represents the most accurate \textit{ab initio} surface to date and will serve as a good starting point for empirical refinement.

Splitting due to tunnelling via the potential energy barrier has played a significant role in the study of molecular spectra since the early days of spectroscopy. The observation of the ammonia doublet led to attempts to find a phosphine analogous, but these have so far failed due to its considerably higher barrier. Full dimensional, variational nuclear motion calculations are used to predict splittings as a function of excitation energy. Simulated spectra suggest that such splittings should be observable in the near infrared via overtones of the $\nu_2$ bending mode starting with $4\nu_2$.

Tau-REx (Tau Retrieval of Exoplanets) is a novel, fully Bayesian atmospheric retrieval code custom built for extrasolar atmospheres. In Waldmann et al. (2015) the transmission spectroscopic case was introduced, here we present the emission spectroscopy spectral retrieval for the Tau-REx framework. Compared to transmission spectroscopy, the emission case is often significantly more degenerate due to the need to retrieve the full atmospheric temperature-pressure (TP) profile. This is particularly true in the case of current measurements of exoplanetary atmospheres, which are either of low signal-to-noise, low spectral resolution or both. Here we present a new way of combining two existing approaches to the modelling of the said TP profile: 1) the parametric profile, where the atmospheric TP structure is analytically approximated by a few model parameters, 2) the Layer-by-Layer approach, where individual atmospheric layers are modelled. Both these approaches have distinct advantages and disadvantages in terms of convergence properties and potential model biases. The Tau-REx hybrid model presented here is a new two-stage TP profile retrieval, which combines the robustness of the analytic solution with the accuracy of the Layer-by-Layer approach. The retrieval process is demonstrated using simulations of the hot-Jupiter WASP-76b and the hot SuperEarth 55 Cnc e, as well as on the secondary eclipse measurements of HD189733b., Comment: ApJ accepted

Spectroscopy of exoplanetary atmospheres has become a well established method for the characterisation of extrasolar planets. We here present a novel inverse retrieval code for exoplanetary atmospheres. TauRex (Tau Retrieval for Exoplanets) is a line-by-line radiative transfer fully Bayesian retrieval framework. TauRex includes the following features: 1) the optimised use of molecular line-lists from the Exomol project; 2) an unbiased atmospheric composition prior selection, through custom built pattern recognition software; 3) the use of two independent algorithms to fully sample the Bayesian likelihood space: nested sampling as well as a more classical Markov Chain Monte Carlo approach; 4) iterative Bayesian parameter and model selection using the full Bayesian Evidence as well as the Savage-Dickey Ratio for nested models, and 5) the ability to fully map very large parameter spaces through optimal code parallelisation and scalability to cluster computing. In this publication we outline the TauRex framework and demonstrate, using a theoretical hot-Jupiter transmission spectrum, the parameter retrieval and model selection. We investigate the impact of Signal-to-Noise and spectral resolution on the retrievability of individual model parameters, both in terms of error bars on the temperature and molecular mixing ratios as well as its effect on the model's global Bayesian evidence., Comment: ApJ published

Blind-source separation techniques are used to extract the transmission spectrum of the hot-Jupiter HD189733b recorded by the Hubble/NICMOS instrument. Such a 'blind' analysis of the data is based on the concept of independent component analysis. The de-trending of Hubble/NICMOS data using the sole assumption that nongaussian systematic noise is statistically independent from the desired light-curve signals is presented. By not assuming any prior, nor auxiliary information but the data themselves, it is shown that spectroscopic errors only about 10 - 30% larger than parametric methods can be obtained for 11 spectral bins with bin sizes of ~0.09 microns. This represents a reasonable trade-off between a higher degree of objectivity for the non-parametric methods and smaller standard errors for the parametric de-trending. Results are discussed in the light of previous analyses published in the literature. The fact that three very different analysis techniques yield comparable spectra is a strong indication of the stability of these results., Comment: ApJ accepted