Source author record

Sergei N. Yurchenko

Sergei N. Yurchenko appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

51works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

51 published item(s)

preprint2026arXiv

ExoMol line lists -- LXIII: ExoMol line lists for 12 isotopologues of CO$_2$

Extensive rovibrational line lists are constructed for 12 isotopologues of carbon dioxide: $^{12}$C$^{16}$O$_2$, $^{13}$C$^{16}$O$_2$, $^{12}$C$^{17}$O$_2$, $^{13}$C$^{17}$O$_2$, $^{12}$C$^{18}$O$_2$, $^{13}$C$^{18}$O$_2$, $^{16}$O$^{12}$C$^{17}$O, $^{16}$O$^{12}$C$^{18}$O, $^{16}$O$^{13}$C$^{17}$O, $^{16}$O$^{13}$C$^{18}$O, $^{17}$O$^{12}$C$^{18}$O, and $^{17}$O$^{13}$C$^{18}$O. The variational program TROVE was employed together with an exact kinetic energy operator, accurate empirical potential energy surface (Ames-2) and the ab initio dipole moment surface Ames-2021-40K. Empirical energy levels from the most recent MARVEL analyses, as well as from the HITRAN and CDSD databases, are used to replace calculated values where available. The line lists are further supplemented by assigning AFGL quantum numbers using machine-learning based estimators. The resulting data were employed to generate opacities with four radiative transfer codes, TauREx, ARCiS, NEMESIS, and petitRADTRANS, both for individual isotopologues and for CO$_2$ at terrestrial isotopic natural abundance. All line lists and associated data are available at www.exomol.com.

preprint2025arXiv

ExoAtom: A Database of Atomic Spectra in ExoMol Format

We present the ExoAtom database, www.exomol.com/exoatom, an extension of the ExoMol database to provide atomic line lists in the ExoMol format. ExoAtom is designed for detailed astrophysical, planetary, and laboratory applications. ExoAtom currently includes atomic data for 80 neutral atoms and 74 singly charged ions. These data are extracted from both the NIST and Kurucz databases, with 79/71 atoms/ions sourced from NIST and 38/37 atoms/ions sourced from Kurucz. ExoAtom uses the file types .all, .def, .states, .trans and .pf as fundamental components for structuring atomic data in a consistent hierarchy. The .states file contains quantum numbers, uncertainties, lifetimes, etc. The .trans file specifies Einstein A coefficients and their associated wavenumbers. The .pf file provides partition functions over a wide grid of temperatures. Post-processing of the ExoAtom data is provided by the program PyExoCross. Future development of ExoAtom will include additional ionization stages.

preprint2024arXiv

The impact of spectral line wing cut-off: Recommended standard method with application to MAESTRO opacity database

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.

preprint2022arXiv

A method for the variational calculation of hyperfine-resolved rovibronic spectra of diatomic molecules

An algorithm for the calculation of hyperfine structure and spectra of diatomic molecules based on the variational nuclear motion is presented. Hyperfine coupling terms considered are Fermi-contact, nuclear spin-electron spin dipole-dipole, nuclear spin-orbit, nuclear spin-rotation and nuclear electric quadrupole interactions. Initial hyperfine-unresolved wavefunctions are obtained for given set of potential energy curves and associated couplings by variation solution of the nuclear-motion Schrödinger equation. Fully hyperfine-resolved parity-conserved rovibronic Hamiltonian matrices for a given final angular momentum, $\bm{F}$, are constructed and then diagonalized to give hyperfine-resolved energies and wavefunctions.Electric transition dipole moment curves can then be used to generate a hyperfine-resolved line list by applying rigorous selection rules.The algorithm is implemented in \textsc{Duo}, which is a general program for calculating spectra of diatomic molecules. This approach is tested for NO and MgH,and the results are compared to experiment and shown to be consistent with those given the well-used effective Hamiltonian code PGOPHER.

preprint2022arXiv

Cross-sections for heavy atmospheres: H$_2$O continuum

Most of the exoplanets detected up to now transit in front of their host stars, allowing for the generation of transmission spectra; the study of exoplanet atmospheres relies heavily upon accurate analysis of these spectra. Recent discoveries mean that the study of atmospheric signals from low-mass, temperate worlds are becoming increasingly common. The observed transit depth in these planets is small and more difficult to analyze. Analysis of simulated transmission spectra for two small, temperate planets (GJ 1214 b and K2-18 b) is presented, giving evidence for significant differences in simulated transit depth when the water vapor continuum is accounted for when compared to models omitting it. These models use cross-sections from the CAVIAR lab experiment for the water self-continuum up to 10,000 cm$^{-1}$; these cross-sections exhibit an inverse relationship with temperature, hence lower-temperature atmospheres are the most significantly impacted. Including the water continuum strongly affects transit depths, increasing values by up to 80 ppm, with the differences for both planets being detectable with the future space missions Ariel and JWST. It is imperative that models of exoplanet spectra move toward adaptive cross-sections, increasingly optimized for H$_2$O-rich atmospheres. This necessitates including absorption contribution from the water vapor continuum into atmospheric simulations.

preprint2022arXiv

Cross-sections for heavy atmospheres: H$_2$O self-broadening

The discovery of super-Earth and mini-Neptune exoplanets means that atmospheric signals from low-mass, temperate exoplanets are being increasingly studied. The signal acquired as the planet transits its host star, known as the transit depth, is smaller for these planets and, as such, more difficult to analyze. The launch of the space telescopes James Webb (JWST) & Ariel will give rise to an explosion in the quality and quantity of spectroscopic data available for an unprecedented number of exoplanets in our galaxy. Accurately extracting the information content, thereby permitting atmospheric science, of such data-sets will require robust models and techniques. We present here the analysis of simulated transmission spectra for water-rich atmospheres, giving evidence for non-negligible differences in simulated transit depths when self-broadening of H$_2$O is correctly accounted for, compared with the currently typically accepted standard of using H$_2$ and He-broadened cross-sections. Our case-study analysis is carried out on two super-Earths, focusing on water-based atmospheres, ranging from H$_2$-rich to H$_2$O-rich. The transit depth is considerably affected, increasing values by up to 60 ppm, which is shown to be detectable with JWST and Ariel. The differences are most pronounced for the lighter (i.e. $mu$ $\sim$ 4) atmospheres. Our work illustrates that it is imperative that the field of exoplanet spectroscopy moves toward adapted cross-sections, increasingly optimized for high-mu atmospheres for studies of super-Earths and mini-Neptunes.

preprint2022arXiv

ExoMol photodissociation cross sections I: HCl and HF

Photon initiated chemistry, \textit{i.e.} the interaction of light with chemical species, is a key factor in the evolution of the atmosphere of exoplanets. For planets orbiting stars in UV-rich environments, photodissociation induced by high energy photons dominates the atmosphere composition and dynamics. The rate of photodissociation can be highly dependent on atmospheric temperature, as increased temperature leads to increased population of vibrational excited states and the consequent lowering of the photodissociation threshold. This paper inaugurates a new series of papers presenting computed temperature-dependent photodissociation cross sections with rates generated for different stellar fields. Cross sections calculations are performed by solving the time-independent Schrödinger equation for each electronic state involved in the process. Here photodissociation cross sections for hydrogen chloride and hydrogen fluoride are computed for a grid of 34 temperatures between 0 and 10~000 K. Use of different radiation fields shows that for the Sun and cooler stars the photodissociation rate can increase exponentially for molecular temperatures above 1000 K; conversely the photodissociation rates in UV rich fields instead are almost insensitive to the temperature of the molecule. Furthermore, these rates show extreme sensitivity to the radiation model used for cool stars, suggesting that further work on these may be required. The provision of an ExoMol database of cross sections is discussed.

preprint2022arXiv

Non-Local thermal equilibrium spectra of atmospheric molecules for exoplanets

Here we present a study of non-LTE effects on the exoplanetary spectra of a collection of molecules which are key in the investigation of exoplanet atmospheres: water, methane, carbon monoxide and titanium oxide. These molecules are chosen as examples of different spectral ranges (IR and UV), molecular types (diatomics and polyatomics) and spectral types (electronic and ro-vibrational); the importance of different vibrational bands in forming distinct non-LTE spectral features are investigated. Most notably, such key spectral signatures for distinguishing between the LTE and non-LTE cases include: for CH4 the 3.15 $μ$m band region; for H2O the 2.0 $μ$m and 2.7 $μ$m band regions; for TiO, a strong variation in intensity in the bands between 0.5 and 0.75 $μ$m; and a sole CO signature between 5 and 6 $μ$m. The analysis is based on the ExoMol cross sections and takes advantage of the extensive vibrational assignment of these molecular line lists in the ExoMol database. We examine LTE and non-LTE cross sections under conditions consistent with those on WASP-12b and WASP-76b using the empirically motivated bi-temperature Treanor model. In addition, we make a simplistic forward model simulation of transmission spectra for H2O in the atmosphere of WASP-12b using the TauREx 3 atmospheric modelling code.

preprint2020arXiv

An update to the MARVEL dataset and ExoMol line list for 12C2

The spectrum of dicarbon (C2) is important in astrophysics and for spectroscopic studies of plasmas and flames. The C2 spectrum is characterized by many band systems with new ones still being actively identified; astronomical observations involve eight of these bands. Recently, Furtenbacher et al. (2016, Astrophys. J. Suppl., 224, 44) presented a set of 5699 empirical energy levels for 12C2, distributed among 11 electronic states and 98 vibronic bands, derived from 42 experimental studies and obtained using the MARVEL (Measured Active Rotational-Vibrational Energy Levels) procedure. Here, we add data from 13 new sources and update data from 5 sources. Many of these data sources characterize high-lying electronic states, including the newly detected 3Pig state. Older studies have been included following improvements in the MARVEL procedure which allow their uncertainties to be estimated. These older works in particular determine levels in the C1Pig state, the upper state of the insufficiently characterized Deslandres-d'Azambuja (C1Pig-A1Piu) band. The new compilation considers a total of 31323 transitions and derives 7047 empirical(MARVEL) energy levels spanning 20 electronic and 142 vibronic states. These new empirical energy levels are used here to update the 8states C2 ExoMol line list. This updated line list is highly suitable for high-resolution cross-correlation studies in astronomical spectroscopy of, for example, exoplanets, as 99.4% of the transitions with intensities over 10^(-18) cm/molecule at 1000 K have frequencies determined by empirical energy levels.

preprint2020arXiv

Empirical Line Lists in the ExoMol Database

The ExoMol database aims to provide comprehensive molecular line lists for exoplanetary and other hot atmospheres. The data are expanded by inclusion of empirically derived line lists taken from the literature for a series of diatomic molecules, namely CH, NH, OH, AlCl, AlF, OH$^+$, CaF, MgF, KF, NaF, LiCl, LiF, MgH, TiH, CrH, FeH, C$_2$, CP, CN, CaH, and triplet N$_2$. Generally, these line lists are constructed from measured spectra using a combination of effective rotational Hamiltonian models for the line positions and ab initio (transition) dipole moments to provide intensities. This work results in the inclusion of 22 new molecules (36 new isotopologues) in the ExoMol database.

preprint2020arXiv

First observation of electric-quadrupole infrared transitions in water vapour

Molecular absorption of infrared radiation is generally due to ro-vibrational electric-dipole transitions. Electric-quadrupole transitions may still occur, but they are typically a million times weaker than electric-dipole transitions, rendering their observation extremely challenging. In polyatomic or polar diatomic molecules, ro-vibrational quadrupole transitions have never been observed. Here, we report the first direct detection of quadrupole transitions in water vapor. The detected quadrupole lines have intensity largely above the standard dipole intensity cut-off of spectroscopic databases and thus are important for accurate atmospheric and astronomical remote sensing.

preprint2020arXiv

Molecular Cross Sections for High Resolution Spectroscopy of Super Earths, Warm Neptunes and Hot Jupiters

High resolution spectroscopy (HRS) has been used to detect a number of species in the atmospheres of hot Jupiters. Key to such detections is accurately and precisely modelled spectra for cross-correlation against the R$\gtrsim$20,000 observations. There is a need for the latest generation of opacities which form the basis for high signal-to-noise detections using such spectra. In this study we present and make publicly available cross sections for six molecular species, H$_2$O, CO, HCN, CH$_4$, NH$_3$ and CO$_2$ using the latest line lists most suitable for low- and high-resolution spectroscopy. We focus on the infrared (0.95-5~$μ$m) and between 500-1500~K where these species have strong spectral signatures. We generate these cross sections on a grid of pressures and temperatures typical for the photospheres of super Earth, warm Neptunes and hot Jupiters using the latest H$_2$ and He pressure broadening. We highlight the most prominent infrared spectral features by modelling three representative exoplanets, GJ~1214~b, GJ~3470~b and HD~189733~b, which encompass a wide range in temperature, mass and radii. In addition, we verify the line lists for H$_2$O, CO and HCN with previous high resolution observations of hot Jupiters. However, we are unable to detect CH$_4$ with our new cross sections from HRS observations of HD~102195~b. These high accuracy opacities are critical for atmospheric detections with HRS and will be continually updated as new data becomes available.

preprint2020arXiv

Spectroscopy of YO from first principles

We report an ab initio study on the spectroscopy of the open-shell diatomic molecule yttrium oxide, YO. The study considers the six lowest doublet states, $X\,{}^{2}Σ^{+}$, $A\,{}^{2}Π$, $C\,{}^{2}Π$, $A'\,{}^{2}Δ$, $B\,{}^{2}Σ^{+}$, $D\,{}^{2}Σ^{+}$ and a few higher-lying quartet states using high levels of electronic structure theory and accurate nuclear motion calculations. The coupled cluster singles, doubles, and perturbative triples, CCSD(T), and multireference configuration interaction (MRCI) methods are employed in conjunction with a relativistic pseudopotential on the yttrium atom and a series of correlation-consistent basis sets ranging in size from triple-$ζ$ to quintuple-$ζ$ quality. Core-valence correlation effects are taken into account and complete basis set limit extrapolation is performed for CCSD(T). Spin-orbit coupling is included through the use of both MRCI state-interaction with spin-orbit (SI-SO) approach and four-component relativistic equation-of-motion CCSD calculations. Using the ab initio data for bond lengths ranging from 1.0 to 2.5 A, we compute 6 potential energy, 12 spin-orbit, 8 electronic angular momentum, 6 electric dipole moment and 12 transition dipole moment (4 parallel and 8 perpendicular) curves which provide a complete description of the spectroscopy of the system of six lowest doublet states. The Duo nuclear motion program is used to solve the coupled nuclear motion Schrödinger equation for these six electronic states. The spectra of $^{89}$Y$^{16}$O simulated for different temperatures are compared with several available high resolution experimental studies; good agreement is found once minor adjustments are made to the electronic excitation energies.

preprint2020arXiv

The 2020 release of the ExoMol database: molecular line lists for exoplanet and other hot atmospheres

The ExoMol database (www.exomol.com) provides molecular data for spectroscopic studies of hot atmospheres. While the data is intended for studies of exoplanets and other astronomical bodies, the dataset is widely applicable. The basic form of the database is extensive line lists; these are supplemented with partition functions, state lifetimes, cooling functions, Landé g-factors, temperature-dependent cross sections, opacities, pressure broadening parameters, $k$-coefficients and dipoles. This paper presents the latest release of the database which has been expanded to consider 80 molecules and 190 isotopologues totaling over 700 billion transitions. While the spectroscopic data is concentrated at infrared and visible wavelengths, ultraviolet transitions are being increasingly considered in response to requests from observers. The core of the database comes from the ExoMol project which primarily uses theoretical methods, albeit usually fine-tuned to reproduce laboratory spectra, to generate very extensive line lists for studies of hot bodies. The data has recently been supplemented by line lists deriving from direct laboratory observations, albeit usually with the use of ab initio transition intensities. A major push in the new release is towards accurate characterisation of transition frequencies for use in high resolution studies of exoplanets and other bodies.

preprint2020arXiv

The ExoMolOP Database: Cross-sections and k-tables for Molecules of Interest in High-Temperature Exoplanet Atmospheres

A publicly available database of opacities for molecules of astrophysical interest, ExoMolOP, has been compiled for over 80 species, based on the latest line list data from the ExoMol, HITEMP and MoLLIST databases. These data are generally suitable for characterising high temperature exoplanet or cool stellar/substellar atmospheres, and have been computed at a variety of pressures and temperatures, with a few molecules included at room-temperature only from the HITRAN database. The data are formatted in different ways for four different exoplanet atmosphere retrieval codes; ARCiS, TauREx, NEMESIS and petitRADTRANS, and include both cross-sections (at R~=~$\fracλ{Δλ}$~=~15,000) and k-tables (at R~=~$\fracλ{Δλ}$~=~1000) for the 0.3~-~50$μ$m wavelength region. Opacity files can be downloaded and used directly for these codes. Atomic data for alkali metals Na and K are also included, using data from the NIST database and the latest line shapes for the resonance lines. Broadening parameters have been taken from the literature where available, or from those for a known molecule with similar molecular properties where no broadening data are available. The data are available from www.exomol.com.

preprint2020arXiv

The high-temperature rotation-vibration spectrum and rotational clustering of silylene (SiH$_2$)

A rotation-vibration line list for the electronic ground state ($\tilde{X}^{1}A_{1}$) of SiH$_2$ is presented. The line list, named CATS, is suitable for temperatures up to 2000 K and covers the wavenumber range 0 - 10,000 cm$^{-1}$(wavelengths $>1.0$ $μ$m) for states with rotational excitation up to $J=52$. Over 310 million transitions between 593 804 energy levels have been computed variationally with a new empirically refined potential energy surface, determined by refining to 75 empirical term values with $J\leq 5$ and a newly computed high-level ab initio dipole moment surface. This is the first, comprehensive high-temperature line list to be reported for SiH$_2$ and it is expected to aid the study of silylene in plasma physics, industrial processes and possible astronomical detection. Furthermore, we investigate the phenomenon of rotational energy level clustering in the spectrum of SiH$_2$. The CATS line list is available from the ExoMol database (www.exomol.com) and the CDS database.

preprint2019arXiv

Analysis of first overtone bands of isotopologues of CO and SiO in stellar spectra

The first overtone ($Δv$ = 2) bands of the monosubstituted isotopologues of CO at 2.3 $μ$m in the spectrum of Arcturus (K2 III), and of the monosubstituted isotopologues of SiO at 4 $μ$m in the spectrum of the red giant HD196610 (M6 III) are modelled. To investigate problems involving the computation of the first overtone bands of isotopologues of CO and SiO in spectra of late-type stars and to determine isotopic abundances. We use fits of theoretical synthetic spectra to the observed stellar molecular bands of CO and SiO to determine abundances of isotopes of C, O and Si. Fits of synthetic spectra of the \CO first overtone bands at 2.3 $μ$m computed with three available line lists (Goorvitch, HITEMP2010 and Li et al.) to the observed spectrum of Arcturus provide the same carbon abundance [C]=$-0.6$ and isotopic ratio of carbon $^{12}$C/$^{13}$C=10 $\pm$ 2. However, the quality of fits to the observed spectrum differ for three line lists used. Furthermore, the derived oxygen isotopic ratio $^{16}$O/$^{18}$O = 2000 $\pm$ 500 is larger than that known in the solar system where $^{16}$O/$^{18}$O = 500. The silicon isotopic ratio in the atmosphere of the red giant HD196610 is revised. Using the ExoMol SiO line list with appropriate statistical weights for the SiO isotopologues the `non-solar' ratio $^{28}$Si:$^{29}$Si:$^{30}$Si = 0.86$\pm$0.03:0.12$\pm$0.02:0.02$\pm$0.01 is obtained. We found that a) the computed isotopic carbon and silicon ratios determined by the fits to the observed spectrum depend on the adopted abundance of C and Si, respectively; b) Correct treatment of the nuclear spin degeneracies parameter is of crucial importance for the use of nowadays HITRAN/ExoMol line lists in the astrophysical computations.

preprint2016arXiv

Ab initio calculations to support accurate modelling of the rovibronic spectroscopy calculations of vanadium monoxide (VO)

Accurate knowledge of the rovibronic near-infrared and visible spectra of vanadium monoxide (VO) is very important for studies of cool stellar and hot planetary atmospheres. Here, the required ab initio dipole moment and spin-orbit coupling curves for VO are produced. This data forms the basis of a new VO line list considering 13 different electronic states and containing over 277 million transitions. Open shell transition, metal diatomics are challenging species to model through ab initio quantum mechanics due to the large number of low-lying electronic states, significant spin-orbit coupling and strong static and dynamic electron correlation. Multi-reference configuration interaction methodologies using orbitals from a complete active space self-consistent-field (CASSCF) calculation are the standard technique for these systems. We use different state-specific or minimal-state CASSCF orbitals for each electronic state to maximise the calculation accuracy. The off-diagonal dipole moment controls the intensity of electronic transitions. We test finite-field off-diagonal dipole moments, but found that (1) the accuracy of the excitation energies were not sufficient to allow accurate dipole moments to be evaluated and (2) computer time requirements for perpendicular transitions were prohibitive. The best off-diagonal dipole moments are calculated using wavefunctions with different CASSCF orbitals.

preprint2016arXiv

Duo: a general program for calculating spectra of diatomic molecules

Duo is a general, user-friendly program for computing rotational, rovibrational and rovibronic spectra of diatomic molecules. Duo solves the Schrödinger equation for the motion of the nuclei not only for the simple case of uncoupled, isolated electronic states (typical for the ground state of closed-shell diatomics) but also for the general case of an arbitrary number and type of couplings between electronic states (typical for open-shell diatomics and excited states). Possible couplings include spin-orbit, angular momenta, spin-rotational and spin-spin. Corrections due to non-adiabatic effects can be accounted for by introducing the relevant couplings using so-called Born-Oppenheimer breakdown curves. Duo requires user-specified potential energy curves and, if relevant, dipole moment, coupling and correction curves. From these it computes energy levels, line positions and line intensities. Several analytic forms plus interpolation and extrapolation options are available for representation of the curves. Duo can refine potential energy and coupling curves to best reproduce reference data such as experimental energy levels or line positions. Duo is provided as a Fortran 2003 program and has been tested under a variety of operating systems.

preprint2016arXiv

ExoMol line lists XV: A new hot line list for hydrogen peroxide

A computed line list for hydrogen peroxide, H$_2{}^{16}$O$_2$, applicable to temperatures up to $T=1250$~K is presented. A semi-empirical high accuracy potential energy surface is constructed and used with an {\it ab initio} dipole moment surface as input TROVE to compute 7.5 million rotational-vibrational states and around 20 billion transitions with associated Einstein-$A$ coefficients for rotational excitations up to $J=85$. The resulting APTY line list is complete for wavenumbers below 6~000 cm$^{-1}$ ($λ< 1.67$~$μ$m) and temperatures up to 1250~K. Room-temperature spectra are compared with laboratory measurements and data currently available in the HITRAN database and literature. Our rms with line positions from the literature is 0.152 \cm\ and our absolute intensities agree better than 10\%. The full line list is available from the CDS databas

preprint2016arXiv

ExoMol line lists XVIII. The high temperature spectrum of VO

An accurate line list, VOMYT, of spectroscopic transitions is presented for hot VO. The 13 lowest electronic states are considered. Curves and couplings are based on initial {\it ab initio} electronic structure calculations and then tuned using available experimental data. Dipole moment curves, used to obtain transition intensities, are computed using high levels of theory (e.g. MRCI/aug-cc-pVQZ using state-specific or minimal-state CAS for dipole moments). This line list contains over 277 million transitions between almost 640,000 energy levels. It covers the wavelengths longer than 0.29 $μ$m and includes all transitions from energy levels within the lowest nine electronic states which have energies less than 20,000 \cm{} to upper states within the lowest 13 electronic states which have energies below 50,000 \cm{}. The line lists give significantly increased absorption at infrared wavelengths compared to currently available VO line lists. The full line lists is made available in electronic form via the CDS database and at www.exomol.com.

preprint2016arXiv

ExoMol molecular line lists - XIV: The rotation-vibration spectrum of hot SO$_2$

Sulphur dioxide is well-known in the atmospheres of planets and satellites, where its presence is often associated with volcanism, and in circumstellar envelopes of young and evolved stars as well as the interstellar medium. This work presents a line list of 1.3 billion $^{32}$S$^{16}$O$_2$ vibration-rotation transitions computed using an empirically-adjusted potential energy surface and an ab initio dipole moment surface. The list gives complete coverage up to 8000 cm$^{-1}$ (wavelengths longer than 1.25 $μ$m) for temperatures below 2000 K. Infrared absorption cross sections are recorded at 300 and 500 C are used to validated the resulting ExoAmes line list. The line list is made available in electronic form as supplementary data to this article and at www.exomol.com.

preprint2016arXiv

ExoMol molecular line lists - XVI: The rotation-vibration spectrum of hot H$_2$S

This work presents the AYT2 line list: a comprehensive list of 114 million $^{1}$H$_2$$^{32}$S vibration-rotation transitions computed using an empirically-adjusted potential energy surface and an {\it ab initio} dipole moment surface. The line list gives complete coverage up to 11000 \cm\ (wavelengths longer than 0.91 $μ$m) for temperatures up to 2000 K. Room temperature spectra can be simulated up to 20000 \cm\ (0.5 $μ$m) but the predictions at visible wavelengths are less reliable. AYT2 is made available in electronic form as supplementary data to this article and at \url{www.exomol.com}.

preprint2016arXiv

ExoMol molecular line lists - XVII The rotation-vibration spectrum of hot SO$_3$

Sulphur trioxide (SO$_3$) is a trace species in the atmospheres of the Earth and Venus, as well as well as being an industrial product and an environmental pollutant. A variational line list for $^{32}$S$^{16}$O$_{3}$, named UYT2, is presented containing 21 billion vibration-rotation transitions. UYT2 can be used to model infrared spectra of SO$_3$ at wavelengths longwards of 2 $μ$m ($ν< 5000$ cm$^{-1}$) for temperatures up to 800 K. Infrared absorption cross sections are also recorded at 300 and 500 C are used to validate the UYT2 line list. The intensities in UYT2 are scaled to match the measured cross sections. The line list is made available in electronic form as supplementary data to this article and at \url{www.exomol.com}.

preprint2016arXiv

Experimental energy levels and partition function of the $^{12}$C$_2$ molecule

The carbon dimer, the $^{12}$C$_2$ molecule, is ubiquitous in astronomical environments. Experimental-quality rovibronic energy levels are reported for $^{12}$C$_2$, based on rovibronic transitions measured for and among its singlet, triplet, and quintet electronic states, reported in 42 publications. The determination utilizes the Measured Active Rotational-Vibrational Energy Levels (MARVEL) technique. The 23,343 transitions measured experimentally and validated within this study determine 5,699 rovibronic energy levels, 1,325, 4,309, and 65 levels for the singlet, triplet, and quintet states investigated, respectively. The MARVEL analysis provides rovibronic energies for six singlet, six triplet, and two quintet electronic states. For example, the lowest measurable energy level of the \astate\ state, corresponding to the $J=2$ total angular momentum quantum number and the $F_1$ spin-multiplet component, is 603.817(5) \cm. This well-determined energy difference should facilitate observations of singlet--triplet intercombination lines which are thought to occur in the interstellar medium and comets. The large number of highly accurate and clearly labeled transitions that can be derived by combining MARVEL energy levels with computed temperature-dependent intensities should help a number of astrophysical observations as well as corresponding laboratory measurements. The experimental rovibronic energy levels, augmented, where needed, with {\it ab initio} variational ones based on empirically adjusted and spin-orbit coupled potential energy curves obtained using the \Duo\ code, are used to obtain a highly accurate partition function, and related thermodynamic data, for $^{12}$C$_2$ up to 4,000 K.

preprint2016arXiv

Pressure-dependent water absorption cross sections for exoplanets and other atmospheres

Many atmospheres (cool stars, brown dwarfs, giant planets, extrasolar planets) are predominately composed of molecular hydrogen and helium. H$_2{}^{16}$O is one of the best measured molecules in extrasolar planetary atmospheres to date and a major compound in the atmospheres of brown-dwarfs and oxygen-rich cool stars, yet the scope of experimental and theoretical studies on the pressure broadening of water vapour lines by collision with hydrogen and helium remains limited. Theoretical H$_2$- and He-broadening parameters of water vapour lines (rotational quantum number $J$ up to 50) are obtained for temperatures in the range 300 - 2000 K. Two approaches for calculation of line widths were used: (i) the averaged energy difference method and (ii) the empirical expression for $J$\p $J$\pp-dependence. Voigt profiles based on these widths and the BT2 line list are used to generate high resolution ($Δ\tildeν$ = 0.01 \cm) pressure broadened cross sections for a fixed range of temperatures and pressures between 300 - 2000 K and 0.001 - 10 bar. An interpolation procedure which can be used to determine cross sections at intermediate temperature and pressure is described. Pressure broadening parameters and cross sections are presented in new ExoMol format.

preprint2016arXiv

Radiative cooling of H3O+ and its deuterated isotopologues

In conjunction with ab initio potential energy and dipole moment surfaces for the electronic ground state, we have made a theoretical study of the radiative lifetimes for the hydronium ion H$_3$O$^{+}$ and its deuterated isotopologues. We compute the ro-vibrational energy levels and their associated wavefunctions together with Einstein coefficients for the electric dipole transitions. A detailed analysis of the stability of the ro-vibrational states have been carried out and the longest-living states of the hydronium ions have been identified. We report estimated radiative lifetimes and cooling functions for temperatures $<$ 200 K. A number of long-living meta-stable states are identified, capable of population trapping.

preprint2016arXiv

Radiative lifetimes and cooling functions for astrophysically important molecules

Extensive line lists generated as part of the ExoMol project are used to compute lifetimes for individual rotational, rovibrational and rovibronic excited states, and temperature-dependent cooling functions by summing over all dipole-allowed transitions for the states concerned. Results are presented for SiO, CaH, AlO, ScH, H$_2$O and methane. The results for CH$_4$ are particularly unusual with 4 excited states with no dipole-allowed decay route and several others where these decays lead to exceptionally long lifetimes. These lifetime data should be useful in models of masers and estimates of critical densities, and can provide a link with laboratory measurements. Cooling functions are important in stellar and planet formation.

preprint2015arXiv

A theoretical room-temperature line list for 15NH3

A new room temperature line list for $^{15}$NH$_3$ is presented. This line list comprised of transition frequencies and Einstein coefficients has been generated using the `spectroscopic' potential energy surface NH3-Y2010 and an ab initio dipole moment surface. The $^{15}$NH$_3$ line list is based on the same computational procedure used for the line list for $^{14}$NH$_3$ BYTe reported recently and should be as accurate. Comparisons with experimental frequencies and intensities are presented. The synthetic spectra show excellent agreement with experimental spectra.

preprint2015arXiv

Accurate Prediction of the Ammonia Probes of a Variable Proton-to-Electron Mass Ratio

A comprehensive study of the mass sensitivity of the vibration-rotation-inversion transitions of $^{14}$NH$_3$, $^{15}$NH$_3$, $^{14}$ND$_3$, and $^{15}$ND$_3$ is carried out variationally using the TROVE approach. Variational calculations are robust and accurate, offering a new way to compute sensitivity coefficients. Particular attention is paid to the $Δk=\pm 3$ transitions between the accidentally coinciding rotation-inversion energy levels of the $ν_2=0^+,0^-,1^+$ and $1^-$ states, and the inversion transitions in the $ν_4=1$ state affected by the "giant" $l$-type doubling effect. These transitions exhibit highly anomalous sensitivities, thus appearing as promising probes of a possible cosmological variation of the proton-to-electron mass ratio $μ$. Moreover, a simultaneous comparison of the calculated sensitivities reveals a sizeable isotopic dependence which could aid an exclusive ammonia detection.

preprint2015arXiv

ExoMol line lists VIII: A variationally computed line list for hot formaldehyde

A computed line list for formaldehyde, H$_2{}^{12}$C$^{16}$O, applicable to temperatures up to $T=1500$~K is presented. An empirical potential energy and {\it ab initio} dipole moment surfaces are used as the input to nuclear motion program TROVE. The resulting line list, referred to as \textit{AYTY}, contains 10.3 million rotational-vibrational states and around 10 billion transition frequencies. Each transition includes associated Einstein-$A$ coefficients and absolute transition intensities, for wavenumbers below 10~000 cm\(^{-1}\) and rotational excitations up to \(J=70\). Room-temperature spectra are compared with laboratory measurements and data currently available in the HITRAN database. These spectra show excellent agreement with experimental spectra and highlight the gaps and limitations of the HITRAN data. The full line list is available from the CDS database as well as at \url{www.exomol.com}.

preprint2015arXiv

ExoMol molecular line lists - XIII: The spectrum of CaO

An accurate line list for calcium oxide is presented covering transitions between all bound ro-vibronic levels from the five lowest electronic states \X, \Ap, \A, \astate, and \bstate. The ro-vibronic energies and corresponding wavefunctionts were obtained by solving the fully coupled Schrödinger equation. \textit{Ab initio} potential energy, spin-orbit, and electronic angular momentum curves were refined by fitting to the experimental frequencies and experimentally derived energies available in the literature. Using our refined model we could (i) reassign the vibronic states for a large portion of the experimentally derived energies [van Groenendael A., Tudorie M., Focsa C., Pinchemel B., Bernath P. F., 2005, J. Mol. Spectrosc., 234, 255], (ii) extended this list of energies to $J=79-118$ and (iii) suggest a new description of the resonances from the \A--\X\ system. We used high level \textit{ab initio} electric dipole moments reported previously [Khalil H., Brites V., Le Quere F., Leonard C., 2011, Chem. Phys., 386, 50] to compute the Einstein A coefficients. Our work is the first fully coupled description of this system. Our line list is the most complete catalogue of spectroscopic transitions available for $^{40}$Ca$^{16}$O and is applicable for temperatures up to at least 5000~K. CaO has yet to be observed astronomically but its transitions are characterised by being particularly strong which should facilitate its detection. The CaO line list is made available in an electronic form as supplementary data to this article and at \url{www.exomol.com}.

preprint2015arXiv

Exomol molecular line lists VI: A high temperature line list for Phosphorus Nitride

Accurate rotational-vibrational line lists for $^{31}$P$^{14}$N and $^{31}$P$^{15}$N in their ground electronic states are computed. The line lists are produced using an empirical potential energy curve obtained by fitting to the experimental transition frequencies available in the literature in conjunction with an accurate, high level \textit{ab initio} dipole moment curve. In these calculations the programs DPotFit and LEVEL~8.0 were used. The new line lists reproduce the experimental wavenumbers with a root-mean-square error of 0.004~cm$^{-1}$. The line lists cover the frequency range 0--51000 cm$^{-1}$, contain almost 700~000 lines each and extend up to a maximum vibrational level of $v$=66 and a maximum rotational level of $J$=357. They should be applicable for a large range of temperature up to, at least, 5000~K. These new line lists are used to simulate spectra for PN at a range of temperatures and are deposited in the Strasbourg data centre. This work is performed as part of the ExoMol project.

preprint2015arXiv

ExoMol molecular line lists X: The spectrum of sodium hydride

Accurate and complete rotational, rotational-vibrational and rotational-vibrational-electronic line lists are calculated for sodium hydride: both the NaH and NaD isotopologues are considered. These line lists cover all ro-vibrational states of the ground ($X$~$^1Σ^+$) and first excited ($A$~$^1Σ^+$) electronic states. The calculations use available spectroscopically-determined potential energy curves and new high-quality, \textit{ab initio} dipole moment curves. Partition functions for both isotopologues are calculated and the effect of quasibound states is considered. The resulting line lists are suitable for temperatures up to about 7000~K and are designed for studies of exoplanet atmospheres, brown dwarfs and cool stars. In particular, the NaH $A-X$ band is found to show a broad absorption feature at about 385 nm which should provide a signature for the molecule. All partition functions, lines and transitions are available as Supplementary Information to this article and at \url{www.exomol.com}.

preprint2015arXiv

ExoMol molecular line lists XII: Line Lists for 8 isotopologues of CS

Comprehensive vibration-rotation line lists for eight isotopologues of carbon monosulphide (CS) ($^{12}$C$^{32}$S, $^{12}$C$^{33}$S, $^{12}$C$^{34}$S, $^{12}$C$^{36}$S, $^{13}$C$^{32}$S, $^{13}$C$^{33}$S, $^{13}$C$^{34}$S, $^{13}$C$^{36}$S) in their ground electronic states are calculated. These line lists are suitable for temperatures up to 3000 K. A spectroscopically-determined potential energy curve (PEC) and dipole moment curve (DMC) are taken from literature. This PEC is adapted to suit our method prior to the computation of ro-vibrational energies. The calculated energies are then substituted by experimental energies, where available, to improve the accuracy of the line lists. The {\it ab initio} DMC is used without refinement to generate Einstein A coefficients. Full line lists of vibration-rotation transitions and partition functions are made available in an electronic form as supporting information to this paper and at \url{www.exomol.com}.

preprint2015arXiv

ExoMol molecular line lists: IX The spectrum of AlO

Accurate line lists are calculated for aluminium monoxide covering the pure rotation, rotation-vibration and electronic (B -- X blue-green and A -- X infrared bands) spectrum. Line lists are presented for the main isotopologue, $^{27}$Al$^{16}$O, as well as for $^{27}$Al$^{17}$O, $^{27}$Al$^{18}$O and $^{26}$Al$^{16}$O. These line lists are suitable for high temperatures (up to 8000 K) including those relevant to exoplanetary atmospheres and cool stars. A combination of empirical and \textit{ab initio} methods is used: the potential energy curves were previously determined to high accuracy by fitting to extensive data from analysis of laboratory spectra; a high quality {\it ab initio} dipole moment curve is calculated using quadruple zeta basis set and the multi-reference configuration interaction (MRCI) method. Partition functions plus full line lists of transitions are made available in an electronic form as supplementary data to this article and at \url{www.exomol.com}.

preprint2015arXiv

ExoMol molecular line lists: XI The spectrum of nitric acid

Nitric acid is a possible biomarker in the atmospheres of exoplanets. An accurate line list of rotational and rotational-vibrational transitions is computed for nitric acid (HNO$_3$). This line list covers wavelengths longer than 1.42 $μ$m (0 - 7000 cm$^{-1}$) and temperatures up to 500 K. The line list is computed using a hybrid variational -- perturbation theory and empirically tuned potential energy and dipole surfaces. It comprises almost 7 billion transitions involving rotations up to $J=100$. Comparisons with spectra from the HITRAN and PNNL databases demonstrate the accuracy of our calculations. Synthetic spectra of water - nitric acid mixturessuggest that nitric acid has features at 7.5 and 11.25 $μ$m that are capable of providing a clear signature for HNO$_3$; the feature at 11.25 $μ$m is particularly promising. Partition functions plus full line lists of transitions are made available in an electronic form as supplementary data to the article and at www.exomol.com.

preprint2015arXiv

Methane and Ammonia in the near-infrared spectra of late T dwarfs

Analysis of T dwarfs using model atmospheres has been hampered by the absence of reliable line lists for methane and ammonia. Newly computed high temperature line lists for both of these important molecules are now available, so it is timely to investigate the appearance of the various absorption features in T dwarfs in order to better understand their atmospheres and validate the new line lists. We present high quality R~5000 Gemini/NIFS 1.0-2.4 microns spectra of the T8 standard 2MASS 0415-0935 and the T9 standard UGPS 0722-0540. We use these spectra to identify numerous methane and ammonia features not previously seen and we discuss the implications for our understanding of T dwarf atmospheres. Among our results, we find that ammonia is the dominant opacity source between ~1.233-1.266 microns in UGPS 0722-0540, and we tentatively identify several absorption features in this wavelength range in the T9's spectrum which may be due entirely to ammonia opacity. Our results also suggest that water rather than methane is the dominant opacity source in the red half of the J-band of the T8 dwarf. Water appears to be the main absorber in this wavelength region in the T9 dwarf until ~1.31 microns, when methane starts to dominate.

preprint2015arXiv

The calculated rovibronic spectrum of scandium hydride, ScH

The electronic structure of six low-lying electronic states of scandium hydride, $X\,{}^{1}Σ^+$, $a\,{}^{3}Δ$, $b\,{}^{3}Π$, $A\,{}^{1}Δ$ $c\,{}^{3}Σ^+$, and $B\,{}^{1}Π$, is studied using multi-reference configuration interaction as a function of bond length. Diagonal and off-diagonal dipole moment, spin-orbit coupling and electronic angular momentum curves are also computed. The results are benchmarked against experimental measurements and calculations on atomic scandium. The resulting curves are used to compute a line list of molecular ro-vibronic transitions for $^{45}$ScH.

preprint2014arXiv

ExoMol line lists IV: The rotation-vibration spectrum of methane up to 1500 K

A new hot line list is calculated for $^{12}$CH$_4$ in its ground electronic state. This line list, called 10to10, contains 9.8 billion transitions and should be complete for temperatures up to 1500 K. It covers the wavelengths longer than 1 $μ$m and includes all transitions to upper states with energies below $hc \cdot 18\,000$ cm$^{-1}$ and rotational excitation up to $J=39$. The line list is computed using the eigenvalues and eigenfunctions of CH$_4$ obtained by variational solution of the Schrödinger equation for the rotation-vibration motion of nuclei employing program TROVE. An ab initio dipole moment surface and a new 'spectroscopic' potential energy surface are used. Detailed comparisons with other available sources of methane transitions including HITRAN, experimental compilations and other theoretical line lists show that these sources lack transitions both higher temperatures and near infrared wavelengths. This line list is suitable for modelling atmospheres of cool stars and exoplanets. It is available from the CDS database as well as at www.exomol.com.

preprint2014arXiv

ExoMol line lists VII: The rotation-vibration spectrum of phosphine up to 1500 K

A comprehensive hot line list is calculated for $^{31}$PH$_3$ in its ground electronic state. This line list, called SAlTY, contains almost 16.8 billion transitions between 7.5 million energy levels and it is suitable for simulating spectra up to temperatures of 1500~K. It covers wavelengths longer than 1~$μ$m and includes all transitions to upper states with energies below $hc \cdot 18\,000$~cm$^{-1}$ and rotational excitation up to $J=46$. The line list is computed by variational solution of the Schrödinger equation for the rotation-vibration motion employing the nuclear-motion program TROVE. A previously reported {\it ab initio} dipole moment surface is used as well as an updated `spectroscopic' potential energy surface (PES), obtained by refining an existing \textit{ab initio} surface through least-squares fitting to the experimentally derived energies. Detailed comparisons with other available sources of phosphine transitions confirms SAlTY's accuracy and illustrates the incompleteness of previous experimental and theoretical compilations for temperatures above 300 K. Atmospheric models are expected to severely underestimate the abundance of phosphine in disequilibrium environments, and it is predicted that phosphine will be detectable in the upper troposphere of many substellar objects. This list is suitable for modelling atmospheres of many astrophysical environments, namely carbon stars, Y dwarfs, T dwarfs, hot Jupiters and solar system gas giant planets. It is available in full as supplementary data to the article and at \url{www.exomol.com}.

preprint2014arXiv

ExoMol molecular line lists V: The ro-vibrational spectra of NaCl and KCl

Accurate rotation-vibration line lists for two molecules, NaCl and KCl, in their ground electronic states are presented. These line lists are suitable for temperatures relevant to exoplanetary atmospheres and cool stars (up to 3000 K). Isotopologues $^{23}$Na$^{35}$Cl, $^{23}$Na$^{37}$Cl, $^{39}$K$^{35}$Cl, $^{39}$K$^{37}$Cl, $^{41}$K$^{35}$Cl and $^{41}$K$^{37}$Cl are considered. Laboratory data was used to refine ab initio potential energy curves in order to compute accurate ro-vibrational energy levels. Einstein A coefficients are generated using newly determined ab initio dipole moment curves calculated using the CCSD(T) method. New Dunham Y$_{ij}$ constants for KCl are generated by a reanalysis of a published Fourier transform infrared emission spectra. Partition functions plus full line lists of ro-vibration transitions are made available in an electronic form as supplementary data to this article and at www.exomol.com.

preprint2014arXiv

The spectrum of hot methane in astronomical objects using a comprehensive computed line list

Hot methane spectra are important in environments ranging from flames to the atmospheres of cool stars and exoplanets. A new spectroscopic line list, 10to10, for $^{12}$CH$_4$ containing almost 10 billion transitions is presented. This comprehensive line list covers a broad spectroscopic range and is applicable for temperatures up to 1 500 K. Previous methane data are incomplete leading to underestimated opacities at short wavelengths and elevated temperatures. Use of 10to10 in models of the bright T4.5 brown dwarf 2MASS 0559-14 leads to significantly better agreement with observations and in studies of the hot Jupiter exoplanet HD 189733b leads to up to a twentifold increase in methane abundance. It is demonstrated that proper inclusion of the huge increase in hot transitions which are important at elevated temperatures is crucial for accurate characterizations of atmospheres of brown dwarfs and exoplanets, especially when observed in the near-infrared.

preprint2014arXiv

The Status of Spectroscopic Data for the Exoplanet Characterisation Missions

The status of laboratory spectroscopic data for exoplanet characterisation missions such as EChO is reviewed. For many molecules (eg H2O, CO, CO2, H3+, O2, O3) the data are already available. For the other species work is actively in progress constructing this data. Much of the is work is being undertaken by ExoMol project (www.exomol.com). This information will be used to construct and EChO-specific spectroscopic database.

preprint2013arXiv

A Computed Room Temperature Line List for Phosphine

An accurate and comprehensive room temperature rotation-vibration transition line list for phosphine (31PH3) is computed using a newly refined potential energy surface and a previously constructed ab initio electric dipole moment surface. Energy levels, Einstein A coefficients and transition intensities are computed using these surfaces and a variational approach to the nuclear motion problem as implemented in the program TROVE. A ro-vibrational spectrum is computed, covering the wavenumber range 0 to 8000 cm-1. The resulting line list, which is appropriate for temperatures up to 300 K, consists of a total of 137 million transitions between 5.6 million energy levels. Several of the band centres are shifted to better match experimental transition frequencies. The line list is compared to the most recent HITRAN database and other laboratorial sources. Transition wavelengths and intensities are generally found to be in good agreement with the existing experimental data, with particularly close agreement for the rotational spectrum. An analysis of the comparison between the theoretical data created and the existing experimental data is performed, and suggestions for future improvements and assignments to the HITRAN database are made.

preprint2013arXiv

Variational calculation of highly excited rovibrational energy levels of H2O2

Results are presented for highly accurate ab initio variational calculation of the rotation - vibration energy levels of H2O2 in its electronic ground state. These results use a recently computed potential energy surface and the variational nuclear-motion programmes WARV4, which uses an exact kinetic energy (EKE) operator, and TROVE, which uses a numerical expansion for the kinetic energy. The TROVE calculations are performed for levels with high values of rotational excitation, $J$ up to 35. The purely \ai\ calculations of the rovibrational energy levels reproduce the observed levels with a standard deviation of about 1 \cm, similar to that of the $J = 0$ calculation as the discrepancy between theory and experiment for rotational energies within a given vibrational state is substantially determined by the error in the vibrational band origin. Minor adjustments are made to the ab initio equilibrium geometry and to the height of the torsional barrier. Using these and correcting the band origins using the error in $J = 0$ states lowers the standard deviation of the observed $-$ calculated energies to only 0.002 \cm\ for levels up to $J = 10$ and 0.02 \cm\ for all experimentally know energy levels, which extend up to $J = 35$.

preprint2013arXiv

Vibrational transition moments of CH$_4$ from first principles

A new nine-dimensional (9D), \textit{ab initio} electric dipole moment surface (DMS) of methane in its ground electronic state is presented. The DMS is computed using an explicitly correlated coupled cluster CCSD(T)-F12 method in conjunction with an F12-optimized correlation consistent basis set of the TZ-family. A symmetrized molecular bond representation is used to parameterise the 9D DMS in terms of sixth-order polynomials. Vibrational transition moments as well as band intensities for a large number of IR-active vibrational bands of $^{12}$CH$_4$ are computed by vibrationally averaging the \textit{ab initio} dipole moment components. The vibrational wavefunctions required for these averages are computed variationally using the program TROVE and a new `spectroscopic' $^{12}$CH$_4$ potential energy surface. The new DMS will be used to produce a hot line list for $^{12}$CH$_4$.

preprint2012arXiv

ExoMol: molecular line lists for exoplanet and other atmospheres

The discovery of extrasolar planets is one of the major scientific advances of the last two decades. Hundreds of planets have now been detected and astronomers are beginning to characterise their composition and physical characteristics. To do this requires a huge quantity of spectroscopic data most of which is not available from laboratory studies. The ExoMol project will offer a comprehensive solution to this problem by providing spectroscopic data on all the molecular transitions of importance in the atmospheres of exoplanets. These data will be widely applicable to other problems and will be used for studies on cool stars, brown dwarfs and circumstellar environments. This paper lays out the scientific foundations of this project and reviews previous work in this area. A mixture of first principles and empirically-tuned quantum mechanical methods will be used to compute comprehensive and very large rotation-vibration and rotation-vibration-electronic (rovibronic) line lists. Methodologies will be developed for treating larger molecules such as methane and nitric acid. ExoMol will rely on these developments and the use of state-of-the-art computing.

preprint2012arXiv

Temperature-dependent molecular absorption cross sections for exoplanets and other atmospheres

Exoplanets, and in particular hot ones such as hot Jupiters, require a very significant quantities of molecular spectroscopic data to model radiative transport in their atmospheres or to interpret their spectra. This data is commonly provided in the form of very extensive transition line lists. The size of these line lists is such that constructing a single model may require the consideration of several billion lines. We present a procedure to simplify this process based on the use of cross sections. Line lists for water, H$_3^+$, HCN /HNC and ammonia have been turned into cross sections on a fine enough grid to preserve their spectroscopic features. Cross sections are provided at a fixed range of temperatures and an interpolation procedure which can be used to generate cross sections at arbitrary temperatures is described. A web-based interface (www.exomol.com/xsec) has been developed to allow astronomers to download cross sections at specified temperatures and spectral resolution. Specific examples are presented for the key water molecule.

preprint2012arXiv

The rovibrational spectrum of BeH, MgH and CaH at high temperatures in the $X\,{}^2Σ^+$ state: a theoretical study

Accurate line lists for three molecules, BeH, MgH and CaH, in their ground electronic states are presented. These line lists are suitable for temperatures relevant to exoplanetary atmospheres and cool stars (up to 2000K). A combination of empirical and \textit{ab initio} methods is used. The rovibrational energy levels of BeH, MgH and CaH are computed using the programs Level and DPotFit in conjunction with `spectroscopic' potential energy curves (PECs). The PEC of BeH is taken from the literature, while the PECs of CaH and MgH are generated by fitting to the experimental transition energy levels. Both spin-rotation interactions (except for BeH, for which it is negligible) and non-adiabatic corrections are explicitly taken into account. Accurate line intensities are generated using newly computed \textit{ab initio} dipole moment curves for each molecule using high levels of theory. Full line lists of rotation-vibration transitions for $^9$BeH, $^{24}$MgH, $^{25}$MgH, $^{26}$MgH and $^{40}$CaH are made available in an electronic form as supplementary data to this article and at \url{www.exomol.com}.

preprint2010arXiv

The discovery of a very cool, very nearby brown dwarf in the Galactic plane

We report the discovery of a very cool, isolated brown dwarf, UGPS 0722-05, with the UKIDSS Galactic Plane Survey. The near-infrared spectrum displays deeper H2O and CH4 troughs than the coolest known T dwarfs and an unidentified absorption feature at 1.275 um. We provisionally classify the object as a T10 dwarf but note that it may in future come to be regarded as the first example of a new spectral type. The distance is measured by trigonometric parallax as d=4.1{-0.5}{+0.6} pc, making it the closest known isolated brown dwarf. With the aid of Spitzer/IRAC we measure H-[4.5] = 4.71. It is the coolest brown dwarf presently known -- the only known T dwarf that is redder in H-[4.5] is the peculiar T7.5 dwarf SDSS J1416+13B, which is thought to be warmer and more luminous than UGPS 0722-05. Our measurement of the luminosity, aided by Gemini/T-ReCS N band photometry, is L = 9.2 +/- 3.1x10^{-7} Lsun. Using a comparison with well studied T8.5 and T9 dwarfs we deduce Teff=520 +/- 40 K. This is supported by predictions of the Saumon & Marley models. With apparent magnitude J=16.52, UGPS 0722-05 is the brightest T dwarf discovered by UKIDSS so far. It offers opportunities for future study via high resolution near-infrared spectroscopy and spectroscopy in the thermal infrared.