A new hot line list for the main isotopologue of CO$_2$ is presented. The line list consists of almost 2.5 billion transitions between 3.5 million rotation-vibration states of CO$_2$ in its ground electronic state, covering the wavenumber range 0-20000 cm$^{-1}$ ($\lambda >0.5$ $\mu$m) with the upper and lower energy thresholds of 36 000 cm$^{-1}$ and 16 000 cm$^{-1}$, respectively. The ro-vibrational energies and wavefunctions are computed variationally using the Ames-2 accurate empirical potential energy surface. The ro-vibrational transition probabilities in the form of Einstein coefficients are computed using an accurate ab initio dipole moment surface using variational program TROVE. A new implementation of TROVE which uses an exact nuclear-motion kinetic energy operator is employed. Comparisons with the existing hot line lists are presented. The line list should be useful for atmospheric retrievals of exoplanets and cool stars. The UCL-4000 line list is available from the CDS and ExoMol databases.

Aims: We survey the transmission spectrum of WASP-121 b for line-absorption by metals and molecules at high spectral resolution, and elaborate on existing interpretations of the optical transmission spectrum observed with HST/STIS and WFC3. Methods: We use the cross-correlation technique and direct differential spectroscopy to search for sodium and other neutral and ionised atoms, TiO, VO and SH in high-resolution transit spectra obtained with the HARPS spectrograph. We inject models assuming chemical and hydrostatic equilibrium with varying temperature and composition to enable model comparison, and employ two bootstrap methods to test the robustness of our detections. Results: We detect neutral Mg, Na, Ca, Cr, Fe, Ni and V, which we predict exists in equilibrium with a significant quantity of VO, supporting earlier observations by HST/WFC3. Non-detections of Ti and TiO support the hypothesis that Ti is depleted via a cold-trap mechanism as has been proposed in the literature. Atomic line depths are under-predicted by hydrostatic models by a factor of 1.5 to 8, confirming recent findings that the atmosphere is extended. We predict the existence of significant concentrations of gas-phase TiO$_2$, VO$_2$ and TiS, which could be important absorbers at optical and NIR wavelengths in hot Jupiter atmospheres, but for which accurate line-list data is currently not available. We find no evidence for absorption by SH, and find that inflated atomic lines can plausibly explain the slope of the transmission spectrum observed in the NUV with HST/STIS. The Na D lines are significantly broadened and show a difference in their respective depths of 15 scale heights, which is not expected from isothermal hydrostatic theory.

The calcium monohydroxide radical (CaOH) is an important astrophysical molecule relevant to cool stars and rocky exoplanets, amongst other astronomical environments. Here, we present a consistent set of highly accurate rovibronic (rotation-vibration-electronic) energy levels for the five lowest electronic states ($\tilde{X}\,^2\Sigma^+$, $\tilde{A}\,^2\Pi$, $\tilde{B}\,^2\Sigma^+$, $\tilde{C}\,^2\Delta$, $\tilde{D}\,^2\Sigma^+$) of CaOH. A comprehensive analysis of the published spectroscopic literature on this system has allowed 1955 energy levels to be determined from 3204 rovibronic experimental transitions, all with unique quantum number labelling and measurement uncertainties. The dataset covers rotational excitation up to $J=62.5$ for molecular states below 29\,000~cm$^{-1}$. The analysis was performed using the MARVEL algorithm, which is a robust procedure based on the theory of spectroscopic networks. The dataset provided will significantly aid future interstellar, circumstellar and atmospheric detections of CaOH, as well as assisting in the design of efficient laser cooling schemes in ultracold molecule research and precision tests of fundamental physics.

Silicon dioxide (SiO$_2$) is expected to occur in the atmospheres of hot rocky super-Earth exoplanets but a lack of spectroscopic data is hampering its possible detection. Here, we present the first, comprehensive molecular line list for SiO$_2$. The line list, named OYT3, covers the wavenumber range 0\,--\,6000~cm$^{-1}$ (wavelengths $\lambda > 1.67$~$\mu$m) and is suitable for temperatures up to $T=3000$~K. Almost 33 billion transitions involving 5.69 million rotation-vibration states with rotational excitation up to $J=255$ have been computed using robust first-principles methodologies. The OYT3 line list is available from the ExoMol database at http://www.exomol.com.

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 (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~$\mu$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.

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.

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$ $\mu$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.

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}\Sigma^{+}$, $A\,{}^{2}\Pi$, $C\,{}^{2}\Pi$, $A'\,{}^{2}\Delta$, $B\,{}^{2}\Sigma^{+}$, $D\,{}^{2}\Sigma^{+}$ 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-$\zeta$ to quintuple-$\zeta$ 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\"{o}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.

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.

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.

The first overtone ($\Delta v$ = 2) bands of the monosubstituted isotopologues of CO at 2.3 $\mu$m in the spectrum of Arcturus (K2 III), and of the monosubstituted isotopologues of SiO at 4 $\mu$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 $\mu$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.

A new hot line list for $^{14}$NH$_3$ is presented. The line list CoYuTe was constructed using an accurate, empirically refined potential energy surface and a CCSD(T)/aug-cc-pVQZ ab initio dipole moment surface of ammonia, previously reported. The line list is an improvement of the ammonia line list BYTe [Yurchenko et al., Mon. Not. R. Astron. Soc., 413, 1828 (2011)]. The CoYuTe line list covers wavenumbers up to 20000 cm$^{-1}$, i.e. wavelengths beyond 0.5 $\mu$m for temperatures up to 1500 K. Comparisons with the high temperature experimental data from the literature show excellent agrement for wavenumbers below 6000 cm$^{-1}$. The CoYuTe line list contains 16.9 billion transitions and is available from the ExoMol website (www.exomol.com) and the CDS database.

The chemical environments of young planets are assumed to be largely influenced by impacts of bodies lingering on unstable trajectories after the dissolution of the protoplanetary disk. We explore the chemical consequences of impacts within the context of reducing planetary atmospheres dominated by carbon monoxide, methane and molecular nitrogen. A terawatt high-power laser was selected in order to simulate the airglow plasma and blast wave surrounding the impactor. The chemical results of these experiments are then applied to a theoretical atmospheric model. The impact simulation results in substantial volume mixing ratios within the reactor of 5% hydrogen cyanide (HCN), 8% acetylene (C$_2$H$_2$), 5% cyanoacetylene (HC$_3$N) and 1% ammonia (NH$_3$). These yields are combined with estimated impact rates for the Early Earth to predict surface boundary conditions for an atmospheric model. We show that impacts might have served as sources of energy that would have led to steady-state surface quantities of 0.4% C$_2$H$_2$, 400 ppm HCN and 40 ppm NH$_3$. We provide simulated transit spectra for an Earth-like exoplanet with this reducing atmosphere during and shortly after eras of intense impacts. We predict that acetylene is as observable as other molecular features on exoplanets with reducing atmospheres that have recently gone through their own `Heavy Bombardments', with prominent features at 3.05 $\mu$m and 10.5 $\mu$m.

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.

The GYT line list covering rotational, rovibrational and rovibronic transitions of the mercapto radical SH is presented. This work extends and replaces the SNaSH line list [Yurchenko et al., 2018, MNRAS, 478, 270] which covers the ground (electronic) $X$ $^{2}\Pi$ state only. This extension is prompted by the tentative identification of the ultra-violet features of SH as being of importance in the transmission spectrum of the ultra-hot Jupiter exoplanet WASP-121b [Evans et al., 2018, AJ., 156, 283]. This GYT line list model is generated by fitting empirical potential energy, spin-orbit and electronic angular momenta functions to experimentally measured wavelengths within the $X$ $^{2}\Pi$ and $A$ $^{2}\Sigma^{+}$ states and to the $A$ $^{2}\Sigma^{+}$ - $X$ $^{2}\Pi$ band system using ab initio curves as a starting reference point. The fits are compatible with the quoted uncertainty of the experimental data used of $\sim$ 0.03 - 0.3 cm$^{-1}$. The GYT line list covers wavelengths longer than 0.256 $\mu$m and includes 7686 rovibronic states and 572 145 transitions for $^{32}$SH. Line lists for the $^{33}$SH, $^{34}$SH, $^{36}$SH and $^{32}$SD isotopologues are generated including a consideration of non-Born-Oppenheimer effects for SD. The line lists are available from the CDS (http://cdsarc.u-strasbg.fr) and ExoMol (www.exomol.com) data bases.

Context: Reliable, directly measured optical properties of astrophysical ice analogs in the infrared (IR) and terahertz (THz) range are missing. These parameters are of great importance to model the dust continuum radiative transfer in dense and cold regions, here thick ice mantles are present, and are necessary for the interpretation of future observations planned in the far-IR region. Aims: Coherent THz radiation allows direct measurement of the complex dielectric function (refractive index) of astrophysically relevant ice species in the THz range. Methods: The time-domain waveforms and the frequency-domain spectra of reference samples of CO ice, deposited at a temperature of 28.5 K and annealed to 33 K at different thicknesses, have been recorded. A new algorithm is developed to reconstruct the real and imaginary parts of the refractive index from the time-domain THz data. Results: The complex refractive index in the wavelength range of 1 mm - 150 ${\mu}$m (0.3 - 2.0 THz) has been determined for the studied ice samples, and compared with available data found in the literature. Conclusions: The developed algorithm of reconstructing the real and imaginary parts of the refractive index from the time-domain THz data enables, for the first time, the determination of optical properties of astrophysical ice analogs without using the Kramers-Kronig relations. The obtained data provide a benchmark to interpret the observational data from current ground based facilities as well as future space telescope missions, and have been used to estimate the opacities of the dust grains in presence of CO ice mantles.

Accurate line lists are crucial for correctly modelling a variety of astrophysical phenomena, including stellar photospheres and the atmospheres of extra-solar planets. This paper presents a new line database Toto for the main isotopologues of titanium oxide (TiO): $^{46}$Ti$^{16}$O, $^{47}$Ti$^{16}$O, $^{48}$Ti$^{16}$O, $^{49}$Ti$^{16}$O and $^{50}$Ti$^{16}$O. The TiO line list contains transitions with wave-numbers up to 30,000 cm$^{-1}$ ie long-wards of 0.33 $\mu$m. The Toto line list includes all dipole-allowed transitions between 13 low-lying electronic states (X $^3\Delta$, a $^1\Delta$, d $^1\Sigma^+$, E $^3\Pi$, A $^3\Phi$ B $^3\Pi$, C $^3\Delta$, b $^1\Pi$, c $^1\Phi$, f $^1\Delta$, e $^1\Sigma^+$). Ab initio potential energy curves (PECs) are computed at the icMRCI level and combined with spin-orbit and other coupling curves. These PECs and couplings are iteratively refined to match known empirical energy levels. Accurate line intensities are generated using ab initio dipole moment curves. The Toto line lists are appropriate for temperatures below 5000 K and contain 30 million transitions for TiO; it is made available in electronic form via the CDS data centre and via www.exomol.com. Tests of the line lists show greatly improved agreement with observed spectra for objects such as M-dwarfs GJ876 and GL581.

A rovibronic line list for the ground ($X$ $^3\Sigma^-$) and first excited ($a$ $^1\Delta$) states of phosphinidene, $^{31}$PH, is computed. The line list is designed for studies of exoplanetary and cool stellar atmospheres with temperatures up to 4000 K. A combination of empirical and ab initio data are used to produce the line list: potential energy curves (PECs) are fitted using experimental transition frequencies; these transitions are reproduced with a root mean square error of 0.01 cm$^{-1}$. The nuclear Schr\"{o}dinger is solved using these PECs plus Born-Oppenheimer and spin splitting correction terms. Line intensities and Einstein $A$ coefficients are computed using ab initio Dipole Moment Curves (DMC) $X$-$X$ and $a$-$a$. The resulting LaTY line list, which contains 65055 transitions for 2528 rovibronic states up to 24500 cm$^{-1}$ and $J=80$ is used to simulate spectra in emission and absorption at a range of temperatures. The line list is made available in electronic form at the CDS and ExoMol databases.

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.

We present the first variational calculation of a hot temperature ab initio line list for the CH$_3$ radical. It is based on a high level ab initio potential energy surface and dipole moment surface of CH$_3$ in the ground electronic state. The ro-vibrational energy levels and Einstein $A$ coefficients were calculated using the general-molecule variational approach implemented in the computer program TROVE. Vibrational energies and vibrational intensities are found to be in very good agreement with the available experimental data. The line list comprises 9,127,123 ro-vibrational states ($J\le 40$) and 2,058,655,166 transitions covering the wavenumber range up to 10000 cm$^{-1}$ and should be suitable for temperatures up to $T= 1500$ K.

Growing astronomical interest in phosphorous (P) chemistry is stimulating the search for new interstellar P-bearing molecules; a task requiring detailed knowledge of the microwave and infrared molecular spectrum. In this work, we present comprehensive rotation-vibration line lists of the \cis- and \trans-isomers of diphosphene (P$_2$H$_2$). The line lists have been generated using robust, first-principles methodologies based on newly computed, high-level \ai\ potential energy and dipole moment surfaces. Transitions are considered between states with energies up to $8000$~cm$^{-1}$ and total angular momentum $J\leq25$. These are the first-ever line lists to be reported for P$_2$H$_2$ and they should significantly facilitate future spectroscopic characterization of this system. The deuterated species \trans-P$_2$HD and the effect of its dynamic dipole moment on the rovibrational spectrum is also discussed.

Line lists for magnesium oxide are computed and extensive comparisons are made with existing experimental spectra. The LiTY line lists cover all ro-vibration transitions within the five lowest-lying electronic states ($X\,{}^1\Sigma^+$, $a\,{}^3\Pi$, $A\,{}^1\Pi$, $B\,{}^1\Sigma^+$ and $b\,{}^3\Sigma^{+}$) and five isotopologues: $^{24}$Mg$^{16}$O, $^{25}$Mg$^{16}$O, $^{26}$Mg$^{16}$O, $^{24}$Mg$^{17}$O, $^{24}$Mg$^{18}$O, $^{24}$Mg$^{17}$O and $^{24}$Mg$^{18}$O. The calculation use potential energy cures, spin-orbit and electronic angular momentum couplings curves determined by fitting to empirical energy levels; these levels are reproduced to within 0.01 \cm\ in most cases. Computed nuclear-motion wavefunctions are combined with {\it ab initio} dipole moment curves to give transition intensities and excited state radiative lifetimes which are compared with laboraroty measurements. The $^{24}$Mg$^{16}$O line list comprises 186 842 ($J\le 320$) ro-vibronic states and 72 833 173 transitions with angular momenta, $J$, up to 300 and covering wavenumbers up to 33 000 cm$^{-1}$ ($\lambda > 0.3$ $\mu$m). The line lists are suitable for temperatures up to about 5000 K. They are relevant to astrophysical studies of exoplanet atmospheres, cool stars and brown dwarfs, and are made available in electronic form at the CDS and ExoMol databases.

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.

Calculations are reported on the rotation-vibration energy levels of the arsine molecule with associated transition intensities. A potential energy surface (PES) obtained from \textit{ab initio} electronic structure calculations is refined to experimental data, and the resulting energy levels display sub-wavenumber accuracy for all reliably known $J=0$ term values under 6500 cm$^{-1}$. After a small empirical adjustment of the band centres, our calculated ($J=1-6$) rovibrational states reproduce 578 experimentally derived energies with a root-mean-square error of 0.122 cm$^{-1}$. Absolute line intensities are computed using the refined PES and a new dipole moment surface (DMS) for transitions between states with energies up to 10~500 cm$^{-1}$ and rotational quantum number $J=30$. The computed DMS reproduces experimental line intensities to within 10\% uncertainty for the $\nu_1$ and $\nu_3$ bands. Furthermore, our calculated absorption cross-sections display good agreement with the main absorption features recorded in Pacific Northwest National Laboratory (PNNL) for the complete range of $600-6500$ cm$^{-1}$.

44325 measured and assigned transitions of H$_2^{32}$S, the parent isotopologue of the hydrogen sulfide molecule, are collated from 33 publications into a single database and reviewed critically. Based on this information, rotation-vibration energy levels are determined for the ground electronic state using the Measured Active Rotational-Vibrational Energy Levels (MARVEL) technique. The ortho and para principal components of the measured spectroscopic network of H$_2^{32}$S are considered separately. The verified set of 25293 ortho- and 18778 para- H$_2^{32}$S transitions determine 3969 ortho and 3467 para energy levels. The Marvel results are compared with alternative data compilations, including a theoretical variational linelist.

Accurate line lists for the carbon dimer, C$_2$, are presented. These line lists cover rovibronic transitions between the eight lowest electronic states: $X\,{}^{1}\Sigma_{g}^{+}$, $a\,{}^{3}\Pi_{u}$, $A\,{}^{1}\Pi_{u}$, $b\,^{3}\Sigma_{g}^{-}$, $b\,^{3}\Sigma_{g}^{-}$, $c\,^{3}\Sigma_{u}^{+}$, $d\,{}^{3}\Pi_{g}$, $B\,{}^{1}\Delta_{g}$, $B^\prime\,{}^{1}\Sigma_{g}^{+}$. Potential energy curves (PECs) and transition dipole moment curves are computed on a large grid of geometries using the aug-cc-pwCVQZ-DK/MRCI level of theory including core and core-valence correlations and scalar relativistic energy corrections. The same level of theory is used to compute spin-orbit and electronic angular momentum couplings. The PECs and couplings are refined by fitting to the empirical (MARVEL) energies of $^{12}$C$_2$ using the nuclear-motion program Duo. The transition dipole moment curves are represented as analytical functions to reduce the numerical noise when computing transition line strengths. Partition functions, full line lists, Land\'{e}-factors and lifetimes for three main isotopologues of C$_2$ ($^{12}$C$_2$,$^{13}$C$_2$ and $^{12}$C$^{13}$C) are made available in electronic form from the CDS (http://cdsarc.u-strasbg.fr) and ExoMol (www.exomol.com) databases.

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).

The ExoMol database (www.exomol.com) provides extensive line lists of molecular transitions which are valid over extended temperatures ranges. The status of the current release of the database is reviewed and a new data structure is specified. This structure augments the provision of energy levels (and hence transition frequencies) and Einstein $A$ coefficients with other key properties, including lifetimes of individual states, temperature-dependent cooling functions, Land\'e $g$-factors, partition functions, cross sections, $k$-coefficients and transition dipoles with phase relations. Particular attention is paid to the treatment of pressure broadening parameters. The new data structure includes a definition file which provides the necessary information for utilities accessing ExoMol through its application programming interface (API). Prospects for the inclusion of new species into the database are discussed.

The mass sensitivity of the vibration-rotation-inversion transitions of H$_3{}^{16}$O$^+$, H$_3{}^{18}$O$^+$, and D$_3{}^{16}$O$^+$ is investigated variationally using the nuclear motion program TROVE~\citep{TROVE:2007}. The calculations utilize new high-level \textit{ab initio} potential energy and dipole moment surfaces. Along with the mass dependence, frequency data and Einstein A coefficients are computed for all transitions probed. Particular attention is paid to the $\Delta|k|=3$ and $\Delta|k-l|=3$ transitions comprising the accidentally coinciding $|J,K\!=\!0,v_2\!=\!0^+\rangle$ and $|J,K\!=\!3,v_2\!=\!0^-\rangle$ rotation-inversion energy levels. The newly computed probes exhibit sensitivities comparable to their ammonia and methanol counterparts, thus demonstrating their potential for testing the cosmological stability of the proton-to-electron mass ratio. The theoretical TROVE results are in close agreement with sensitivities obtained using the nonrigid and rigid inverter approximate models, confirming that the \textit{ab initio} theory used in the present study is adequate.

A robust variational approach is used to investigate the sensitivity of the rotation-vibration spectrum of phosphine (PH$_3$) to a possible cosmological variation of the proton-to-electron mass ratio, $\mu$. Whilst the majority of computed sensitivity coefficients, $T$, involving the low-lying vibrational states acquire the expected values of $T\approx-1$ and $T\approx-1/2$ for rotational and ro-vibrational transitions, respectively, anomalous sensitivities are uncovered for the $A_1\!-\!A_2$ splittings in the $\nu_2/\nu_4$, $\nu_1/\nu_3$ and $2\nu_4^{\ell=0}/2\nu_4^{\ell=2}$ manifolds of PH$_3$. A pronounced Coriolis interaction between these states in conjunction with accidentally degenerate $A_1$ and $A_2$ energy levels produces a series of enhanced sensitivity coefficients. Phosphine is expected to occur in a number of different astrophysical environments and has potential for investigating a drifting constant. Furthermore, the displayed behaviour hints at a wider trend in molecules of ${\bf C}_{3\mathrm{v}}\mathrm{(M)}$ symmetry, thus demonstrating that the splittings induced by higher-order ro-vibrational interactions are well suited for probing $\mu$ in other symmetric top molecules in space, since these low-frequency transitions can be straightforwardly detected by radio telescopes.

A new line list for H$_2$$^{16}$O is presented. This line list, which is called POKAZATEL, includes transitions between rotation-vibrational energy levels up to 41000 cm$^{-1}$ in energy and is the most complete to date. The potential energy surface (PES) used for producing the line list was obtained by fitting a high-quality ab initio PES to experimental energy levels with energies of 41000 cm$^{-1}$ and for rotational excitations up to $J=5$. The final line list comprises all energy levels up to 41000 cm$^{-1}$ and rotational angular momentum $J$ up to 72. An accurate ab initio dipole moment surface (DMS) was used for the calculation of line intensities and reproduces high-precision experimental intensity data with an accuracy close to 1 %. The final line list uses empirical energy levels whenever they are available, to ensure that line positions are reproduced as accurately as possible. The POKAZATEL line list contains over 5 billion transitions and is available from the ExoMol website (www.exomol.com) and the CDS database.

SiS has long been observed in the circumstellar medium of the carbon-rich star IRC+10216 CW Leo. Comprehensive and accurate rotation-vibrational line lists and partition functions are computed for 12 isotopologues of silicon sulphide ($^{28}$Si$^{32}$S, $^{28}$Si$^{34}$S, $^{29}$Si$^{32}$S, $^{28}$Si$^{33}$S, $^{30}$Si$^{32}$S, $^{29}$Si$^{34}$S, $^{30}$Si$^{34}$S, $^{28}$Si$^{36}$S, $^{29}$Si$^{33}$S, $^{29}$Si$^{36}$S, $^{30}$Si$^{33}$S, $^{30}$Si$^{36}$S) in its ground (X $^1\Sigma^+$) electronic state. The calculations employ an existing spectroscopically accurate potential energy curve (PEC) derived from experimental measurements and a newly-computed ab initio dipole moment curve (DMC). The $^{28}$Si$^{32}$S line list includes 10104 states and 91715 transitions. These line lists are available from the ExoMol website www.exomol.com and the CDS database.

Comprehensive rotation-vibration line lists are presented for the two main isotopologues of methyl chloride, $^{12}$CH$_3{}^{35}$Cl and $^{12}$CH$_3{}^{37}$Cl. The line lists, OYT-35 and OYT-37, are suitable for temperatures up to $T=1200\,$K and consider transitions with rotational excitation up to $J=85$ in the wavenumber range $0$--$6400\,$cm$^{-1}$ (wavelengths $\lambda> 1.56\,\mu$m). Over 166 billion transitions between 10.2 million energy levels have been calculated variationally for each line list using a new empirically refined potential energy surface, determined by refining to 739 experimentally derived energy levels up to $J=5$, and an established {\it ab initio} dipole moment surface. The OYT line lists show excellent agreement with newly measured high-temperature infrared absorption cross-sections, reproducing both strong and weak intensity features across the spectrum. The line lists are available from the ExoMol database and the CDS database.

A new line list for ethylene, $^{12}$C$_2$$^1$H$_4$ is presented. The line list is based on high level ab initio potential energy and dipole moment surfaces. The potential energy surface is refined by fitting to experimental energies. The line list covers the range up to 7000 cm$^{-1}$ (1.43 $\mu$m) with all ro-vibrational transitions (50 billion) with the lower state below 5000 cm$^{-1}$ included and thus should be applicable for temperatures up to 700 K. A technique for computing molecular opacities from vibrational band intensities is proposed and used to provide temperature dependent cross sections of ethylene for shorter wavelength and higher temperatures. When combined with realistic band profiles (such as the proposed three-band model), the vibrational intensity technique offers a cheap but reasonably accurate alternative to the full ro-vibrational calculations at high temperatures and should be reliable for representing molecular opacities. The C$_2$H$_4$ line list, which is called MaYTY, is made available in electronic form from the CDS.

A new line list for AlH is produced. The WYLLoT line list spans two electronic states $X\,{}^1\Sigma^+$ and $A\,{}^1\Pi$. A diabatic model is used to model the shallow potential energy curve of the $A\,{}^1\Pi$ state, which has a strong pre-dissociative character with only two bound vibrational states. Both potential energy curves are empirical and were obtained by fitting to experimentally derived energies of the $X\,{}^1\Sigma^+$ and $A\,{}^1\Pi$ electronic states using the diatomic nuclear motion codes Level and Duo. High temperature line lists plus partition functions and lifetimes for three isotopologues $^{27}$AlH, $^{27}$AlD and $^{26}$AlH were generated using ab initio dipole moments. The line lists cover both the $X$--$X$ and $A$--$X$ systems and are made available in electronic form at the CDS and ExoMol databases.

The ExoMol project is dedicated to providing molecular line lists for exoplanet and other hot atmospheres. The ExoMol procedure uses a mixture of ab initio calculations and available laboratory data. The actual line lists are generated using variational nuclear motion calculations. These line lists form the input for opacity models for cool stars and brown dwarfs as well as for radiative transport models involving exoplanets. This paper is a collection of molecular opacities for 52 molecules (130 isotopologues) at two reference temperatures, 300 K and 2000 K, using line lists from the ExoMol database. So far, ExoMol line lists have been generated for about 30 key molecular species. Other line lists are taken from external sources or from our work predating the ExoMol project. An overview of the line lists generated by ExoMol thus far is presented and used to evaluate further molecular data needs. Other line lists are also considered. The requirement for completeness within a line list is emphasized and needs for further line lists discussed.

Line lists for the sulphur-containing molecules SH (the mercapto radical) and NS are computed as part of the ExoMol project. These line lists consider transitions within the $X$ ${}^2\Pi$ ground state for $^{32}$SH, $^{33}$SH, $^{34}$SH and $^{\text{32}}$SD, and $^{14}$N$^{32}$S, $^{14}$N$^{33}$S, $^{14}$N$^{34}$S, $^{14}$N$^{36}$S and $^{15}$N$^{32}$S. Ab initio potential energy (PEC) and spin-orbit coupling (SOC) curves are computed and then improved by fitting to experimentally observed transitions. Fully ab initio dipole moment curves (DMCs) computed at high level of theory are used to produce the final line lists. For SH, our fit gives a root-mean-square (rms) error of 0.03 cm$^{-1}$ between the observed ($v_{\rm max}=4$, $J_{\rm max} = 34.5$) and calculated transitions wavenumbers; this is extrapolated such that all $X$ $^2\Pi$ rotational-vibrational-electronic (rovibronic) bound states are considered. For $^{\text{32}}$SH the resulting line list contains about 81000 transitions and 2300 rovibronic states, considering levels up to $v_{\rm max} = 14$ and $J_{\rm max} = 60.5$. For NS the refinement used a combination of experimentally determined frequencies and energy levels and led to an rms fitting error of 0.002 cm$^{-1}$. Each NS calculated line list includes around 2.8 million transitions and 31000 rovibronic states with a vibrational range up to $v=53$ and rotational range to $J=235.5$, which covers up to 23000 cm$^{-1}$. Both line lists should be complete for temperatures up to 5000 K. Example spectra simulated using this line list are shown and comparisons made to the existing data in the CDMS database. The line lists are available from the CDS (http://cdsarc.u-strasbg.fr) and ExoMol (www.exomol.com) data bases.

ExoCross is a Fortran code for generating spectra (emission, absorption) and thermodynamic properties (partition function, specific heat etc.) from molecular line lists. Input is taken in several formats, including ExoMol and HITRAN formats. ExoCross is efficiently parallelized showing also a high degree of vectorization. It can work with several line profiles such as Doppler, Lorentzian and Voigt and support several broadening schemes. Voigt profiles are handled by several methods allowing fast and accurate simulations. Two of these methods are new. ExoCross is also capable of working with the recently proposed method of super-lines. It supports calculations of lifetimes, cooling functions, specific heats and other properties. ExoCross can be used to convert between different formats, such as HITRAN, ExoMol and Phoenix. It is capable of simulating non-LTE spectra using a simple two-temperature approach. Different electronic, vibronic or vibrational bands can be simulated separately using an efficient filtering scheme based on the quantum numbers.

We present here the analysis of 30 gaseous extrasolar planets, with temperatures between 600 and 2400 K and radii between 0.35 and 1.9 $R_\mathrm{Jup}$. The quality of the HST/WFC3 spatially scanned data combined with our specialized analysis tools allow us to study the largest and most self-consistent sample of exoplanetary transmission spectra to date and examine the collective behavior of warm and hot gaseous planets rather than isolated case-studies. We define a new metric, the Atmospheric Detectability Index (ADI) to evaluate the statistical significance of an atmospheric detection and find statistically significant atmospheres around 16 planets out of the 30 analysed. For most of the Jupiters in our sample, we find the detectability of their atmospheres to be dependent on the planetary radius but not on the planetary mass. This indicates that planetary gravity plays a secondary role in the state of gaseous planetary atmospheres. We detect the presence of water vapour in all of the statistically detectable atmospheres, and we cannot rule out its presence in the atmospheres of the others. In addition, TiO and/or VO signatures are detected with 4$\sigma$ confidence in WASP-76 b, and they are most likely present in WASP-121 b. We find no correlation between expected signal-to-noise and atmospheric detectability for most targets. This has important implications for future large-scale surveys.

A variationally computed $^{28}$SiH$_4$ rotation-vibration line list applicable for temperatures up to $T=1200\,$K is presented. The line list, called OY2T, considers transitions with rotational excitation up to $J=42$ in the wavenumber range $0$--$5000\,$cm$^{-1}$ (wavelengths $\lambda> 2\,\mu$m). Just under 62.7 billion transitions have been calculated between 6.1 million energy levels. Rovibrational calculations have utilized a new `spectroscopic' potential energy surface determined by empirical refinement to 1452 experimentally derived energy levels up to $J=6$, and a previously reported \textit{ab initio} dipole moment surface. The temperature-dependent partition function of silane, the OY2T line list format, and the temperature dependence of the OY2T line list are discussed. Comparisons with the PNNL spectral library and other experimental sources indicate that the OY2T line list is robust and able to accurately reproduce weaker intensity features. The full line list is available from the ExoMol database and the CDS database.

SiH has long been observed in the spectrum of our Sun and other cool stars. Computed line lists for the main isotopologues of silicon monohydride, $^{28}$SiH, $^{29}$SiH, $^{30}$SiH and $^{28}$SiD are presented. These line lists consider rotation-vibration transitions within the ground $X$ $^{2}\Pi$ electronic state as well as transitions to the low-lying $A$ $^{2}\Delta$ and $a$ $^{4}\Sigma^-$ states. Ab initio potential energy (PECs) and dipole moment curves (DMCs) along with spin-orbit and electronic-angular-momentum couplings between them are calculated using the MRCI level of theory with the MOLPRO package. The PEC for the ground X $^2\Pi$ state is refined to available experimental data with a typical accuracy of around 0.01 cm$^{-1}$ or better. The $^{28}$SiH line list includes 11,785 rovibronic states and 1,724,841 transitions with associated Einstein-A coefficients for angular momentum $J$ up to $82.5$ and covering wavenumbers up to 31340 cm$^{-1}$ ($\lambda$ $<$ 0.319 $\mu$m). Spectra are simulated using the new line list and comparisons made with various experimental spectra. These line lists are applicable up to temperatures of ~5000 K, making them relevant to astrophysical objects such as exoplanetary atmospheres and cool stars and opening up the possibility of detection in the interstellar medium. These line lists are available at the ExoMol www.exomol.com and CDS database websites.

Rotation-vibration energy levels are determined for the electronic ground state of the acetylene molecule, $^{12}$C$_2$H$_2$, using the Measured Active Rotational-Vibrational Energy Levels (MARVEL) technique. 37,813 measured transitions from 61 publications are considered. The distinct components of the spectroscopic network linking ortho and para states are considered separately. The 20,717 ortho and 17,096 para transitions measured experimentally are used to determine 6013 ortho and 5200 para energy levels. The MARVEL results are compared with alternative compilations based on the use of effective Hamiltonians.

Comprehensive line lists for phosphorus monoxide ($^{31}$P$^{16}$O) and phosphorus monosulphide ($^{31}$P$^{32}$S) in their $X$ $^2\Pi$ electronic ground state are presented. The line lists are based on new ab initio potential energy (PEC), spin-orbit (SOC) and dipole moment (DMC) curves computed using the MRCI+Q-r method with aug-cc-pwCV5Z and aug-cc-pV5Z basis sets. The nuclear motion equations (i.e. the rovibronic Schr\"odinger equations for each molecule) are solved using the program Duo. The PECs and SOCs are refined in least-squares fits to available experimental data. Partition functions, $Q(T)$, are computed up to $T=$ 5000 K, the range of validity of the line lists. These line lists are the most comprehensive available for either molecule. The characteristically sharp peak of the $Q$-branches from the spin-orbit split components give useful diagnostics for both PO and PS in spectra at infrared wavelengths. These line lists should prove useful for analysing observations and setting up models of environments such as brown dwarfs, low-mass stars, O-rich circumstellar regions and potentially for exoplanetary retrievals. Since PS is yet to be detected in space, the role of the two lowest excited electronic states ($a$ $^4\Pi$ and $B$ $^2\Pi$) are also considered. An approximate line list for the PS $X - B$ electronic transition, which predicts a number of sharp vibrational bands in the near ultraviolet, is also presented. he line lists are available from the CDS (cdsarc.u-strasbg.fr) and ExoMol (www.exomol.com) databases.

Molecular line lists are important for modelling absorption and emission processes in atmospheres of different astronomical objects, such as cool stars and exoplanets. In order to be applicable for high temperatures, line lists for molecules like methane must contain billions of transitions, which makes their direct (line-by-line) application in radiative transfer calculations impracticable. Here we suggest a new, hybrid line list format to mitigate this problem, based on the idea of temperature-dependent absorption continuum. The line list is partitioned into a large set of relatively weak lines and a small set of important, stronger lines. The weaker lines are then used either to construct a temperature-dependent (but pressure-independent) set of intensity cross sections or are blended into a greatly reduced set of super-lines. The strong lines are kept in the form of temperature independent Einstein A coefficients. A line list for methane is constructed as a combination of 17 million strong absorption lines relative to the reference absorption spectra and a background methane continuum in two temperature-dependent forms, of cross sections and super-lines. This approach eases the use of large high temperature line lists significantly as the computationally expensive calculation of pressure dependent profiles only need to be performed for a relatively small number of lines. Both the line list and cross sections were generated using a new 34 billion methane line list (34to10), which extends the 10to10 line list to higher temperatures (up to 2000 K). The new hybrid scheme can be applied to any large line lists containing billions of transitions. We recommend to use super-lines generated on a high resolution grid based on resolving power (R = 1,000,000) to model the molecular continuum as a more flexible alternative to the temperature dependent cross sections.

The majority of stars are now thought to support exoplanets. Many of those exoplanets discovered thus far are categorized as rocky objects with an atmosphere. Most of these objects are however hot due to their short orbital period. Models suggest that water is the dominant species in their atmospheres. The hot temperatures are expected to turn these atmospheres into a (high pressure) steam bath containing remains of melted rock. The spectroscopy of these hot rocky objects will be very different from that of cooler objects or hot gas giants. Molecules suggested to be important for the spectroscopy of these objects are reviewed together with the current status of the corresponding spectroscopic data. Perspectives of building a comprehensive database of linelist/cross sections applicable for atmospheric models of rocky super-Earths as part of the ExoMol project are discussed. The quantum-mechanical approaches used in linelist productions and their challenges are summarized.

Line lists for the ground electronic ground state for six major isotopologues of nitric oxide are presented. The line lists are constructed using empirical energy levels (and line positions) and high-level {\it ab inito} intensities. The energy levels were obtained using a combination of two approaches, from an effective Hamiltonian and from solving the rovibronic Schr\"{o}dinger equation variationally. The effective hamiltonian model was obtained through a fit to the experimental line positions of NO available in the literature for all six isotopologues using the programs SPFIT and SPCAT. The variational model was built through a least squares fit of the \textit{ab inito} potential and spin-orbit curves to the experimentally derived energies and experimental line positions of the main isotopologue only using the Duo program. The \textit{ab inito} potential energy, spin-orbit and dipole moment curves (PEC, SOC and DMC) are computed using high-level {\it ab inito} methods and the MARVEL method is used to obtain energies of NO from experimental transition frequencies. Each line list covers a wavenumber range from 0 - 40,000 \cm with approximately 22,000 rovibronic states and 2.3-2.6 million transitions extending to $J_{max} = 184.5$ and $v_{max} = 51$. Partition functions are also calculated up to a temperature of 5000 K. The calculated absorption line intensities at 296 K using these line lists show excellent agreement with those included in the HITRAN and HITEMP databases. The computed NO line lists are the most comprehensive to date, covering a wider wavenumber and temperature range compared to both the HITRAN and HITEMP databases. These line lists are also more accurate than those used in HITEMP. The full line lists are available from the CDS and ExoMol databases; data will also be available from CDMS.

We present a model for lightning shock induced chemistry that can be applied to atmospheres of arbitrary H/C/N/O chemistry, hence for extrasolar planets and brown dwarfs. The model couples hydrodynamics and the STAND2015 kinetic gas-phase chemistry. For an exoplanet analogue to the contemporary Earth, our model predicts NO and NO2 yields in agreement with observation. We predict height-dependent mixing ratios during a storm soon after a lightning shock of NO ~ 1e-3 at 40 km and NO2 ~ 1e-4 below 40 km, with O3 reduced to trace quantities (<< 1e-10). For an Earth-like exoplanet with a CO2/N2 dominated atmosphere and with an extremely intense lightning storm over its entire surface, we predict significant changes in the amount of NO, NO2, O3, H2O, H2, and predict significant abundance of C2N. We find that, for the Early Earth, O2 is formed in large quantities by lightning but is rapidly processed by the photochemistry, consistent with previous work on lightning. The effect of persistent global lightning storms are predicted to be significant, primarily due to NO2, with the largest spectral features present at ~3.4 {\mu}m and ~6.2 {\mu}m. The features within the transmission spectrum are on the order of 1 ppm and therefore are not likely detectable with JWST. Depending on its spectral properties, C2N could be a key tracer for lightning on Earth-like exoplanets with a N2/CO2 bulk atmosphere, unless destroyed by yet unknown chemical reactions.

H$_3^+$ is a ubiquitous and important astronomical species whose spectrum has been observed in the interstellar medium, planets and tentatively in the remnants of supernova SN1897a. Its role as a cooler is important for gas giant planets and exoplanets, and possibly the early Universe. All this makes the spectral properties, cooling function and partition function of H$_3^+$ key parameters for astronomical models and analysis. A new high-accuracy, very extensive line list for H$_3^+$ called MiZATeP was computed as part of the ExoMol project alongside a temperature-dependent cooling function and partition function as well as lifetimes for %individual excited states. These data are made available in electronic form as supplementary data to this article and at http://www.exomol.com

We describe the recent updates to the ExoMol database regarding the molecular spectral line shapes. ExoMol provides comprehensive molecular line lists with a special emphasis on the applications involving characterization of atmospheres of exoplanets and cool stars. Among important requirements of such applications are (i) the broadening parameters for hydrogen and helium dominating atmospheres and (ii) very broad ranges of temperature and pressures. The current status of the available line shape data in the literature, demands from the exoplanetary community and their specific needs are discussed.

Hot line lists for two isotopologues of water, \octo\ and \heto, are presented. The calculations employ newly constructed potential energy surfaces (PES) which take advantage of a novel method for using the large set of experimental energy levels for \hato\ to give high quality predictions for \octo\ and \heto. This procedure greatly extends the energy range for which a PES can be accurately determined, allowing accurate prediction of higher-lying energy levels than are currently known from direct laboratory measurements. This PES is combined with a high-accuracy, {\it ab initio} dipole moment surface of water in the computation of all energy levels, transition frequencies and associated Einstein A coefficients for states with rotational excitation up to $J=50$ and energies up to 30~000 \cm. The resulting HotWat78 line lists complement the well-used BT2 \hato\ line list (Barber et.al, 2006, MNRAS, {\bf 368}, 1087). Full line lists are made available in the electronic form as supplementary data to this article and at \url{www.exomol.com}.

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 ($\Delta \tilde{\nu}$ = 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.

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 $\mu$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.

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.

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 $\mu$m ($\nu < 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}.

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}$ ($\lambda < 1.67$~$\mu$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

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 $\mu$m) for temperatures up to 2000 K. Room temperature spectra can be simulated up to 20000 \cm\ (0.5 $\mu$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}.

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.

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 $\mu$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.

We report the analysis of two new spectroscopic observations of the super-Earth 55 Cancri e, in the near infrared, obtained with the WFC3 camera onboard the HST. 55 Cancri e orbits so close to its parent star, that temperatures much higher than 2000 K are expected on its surface. Given the brightness of 55 Cancri, the observations were obtained in scanning mode, adopting a very long scanning length and a very high scanning speed. We use our specialized pipeline to take into account systematics introduced by these observational parameters when coupled with the geometrical distortions of the instrument. We measure the transit depth per wavelength channel with an average relative uncertainty of 22 ppm per visit and find modulations that depart from a straight line model with a 6$\sigma$ confidence level. These results suggest that 55 Cancri e is surrounded by an atmosphere, which is probably hydrogen-rich. Our fully Bayesian spectral retrieval code, T-REx, has identified HCN to be the most likely molecular candidate able to explain the features at 1.42 and 1.54 $\mu$m. While additional spectroscopic observations in a broader wavelength range in the infrared will be needed to confirm the HCN detection, we discuss here the implications of such result. Our chemical model, developed with combustion specialists, indicates that relatively high mixing ratios of HCN may be caused by a high C/O ratio. This result suggests this super-Earth is a carbon-rich environment even more exotic than previously thought.

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.

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\"{o}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}.

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.

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 $\mu$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 $\mu$m that are capable of providing a clear signature for HNO$_3$; the feature at 11.25 $\mu$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.

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}.

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}.

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\Sigma^+$) and first excited ($A$~$^1\Sigma^+$) 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}.

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 $\Delta k=\pm 3$ transitions between the accidentally coinciding rotation-inversion energy levels of the $\nu_2=0^+,0^-,1^+$ and $1^-$ states, and the inversion transitions in the $\nu_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 $\mu$. Moreover, a simultaneous comparison of the calculated sensitivities reveals a sizeable isotopic dependence which could aid an exclusive ammonia detection.

The electronic structure of six low-lying electronic states of scandium hydride, $X\,{}^{1}\Sigma^+$, $a\,{}^{3}\Delta$, $b\,{}^{3}\Pi$, $A\,{}^{1}\Delta$ $c\,{}^{3}\Sigma^+$, and $B\,{}^{1}\Pi$, 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.

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}.

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.

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.

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.

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.

The discovery of almost 2000 exoplanets has revealed an unexpectedly diverse planet population. Observations to date have shown that our Solar System is certainly not representative of the general population of planets in our Milky Way. The key science questions that urgently need addressing are therefore: What are exoplanets made of? Why are planets as they are? What causes the exceptional diversity observed as compared to the Solar System? EChO (Exoplanet Characterisation Observatory) has been designed as a dedicated survey mission for transit and eclipse spectroscopy capable of observing a large and diverse planet sample within its four-year mission lifetime. EChO can target the atmospheres of super-Earths, Neptune-like, and Jupiter-like planets, in the very hot to temperate zones (planet temperatures of 300K-3000K) of F to M-type host stars. Over the next ten years, several new ground- and space-based transit surveys will come on-line (e.g. NGTS, CHEOPS, TESS, PLATO), which will specifically focus on finding bright, nearby systems. The current rapid rate of discovery would allow the target list to be further optimised in the years prior to EChO's launch and enable the atmospheric characterisation of hundreds of planets. Placing the satellite at L2 provides a cold and stable thermal environment, as well as a large field of regard to allow efficient time-critical observation of targets randomly distributed over the sky. A 1m class telescope is sufficiently large to achieve the necessary spectro-photometric precision. The spectral coverage (0.5-11 micron, goal 16 micron) and SNR to be achieved by EChO, thanks to its high stability and dedicated design, would enable a very accurate measurement of the atmospheric composition and structure of hundreds of exoplanets.

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~$\mu$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\"odinger 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}.

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.

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.

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.

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 $\mu$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\"{o}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.

A revised rotation-vibration line list for the combined hydrogen cyanide (HCN) / hydrogen isocyanide (HNC) system is presented. The line list uses {\it ab initio} transition intensities calculated previously (Harris et al., ApJ, 2002, 578, 657) and extensive datasets of recently measured experimental energy levels (Mellau, J. Chem. Phys. and J. Mol. Spectrosc. 2010-2011). The resulting line list has significantly more accurate wavelengths than previous ones for these systems. An improved value for the separation between HCN and HNC is adopted leading to an approximately 25%\ lower predicted thermal population of HNC as a function of temperature in the key 2000 to 3000 K region. Temperature-dependent partition functions and equilibrium constants are presented. The line lists are validated by comparison with laboratory spectra and are presented in full as supplementary data to the article and at \url{www.exomol.com}.

Accurate rotation-vibration line lists are calculated for silicon monoxide. Line lists are presented for the main isotopologue, $^{28}$Si$^{16}$O, and for four monosubsituted isotopologues ($^{29}$Si$^{16}$O, $^{30}$Si$^{16}$O, $^{28}$Si$^{18}$O and $^{28}$Si$^{17}$O), in their ground electronic states. These line lists are suitable for high temperatures (up to 9000 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 are determined to high accuracy by fitting to extensive data from the analysis of both laboratory and sunspot spectra; a high quality {\it ab initio} dipole moment curve is calculated at the large basis set, multi-reference configuration interaction level. A partition function plus full line lists of rotation-vibration transitions are made available in an electronic form as supplementary data to this article and at \url{www.exomol.com}.

We report near-infrared measurements of the terminator region transmission spectrum and dayside emission spectrum of the exoplanet WASP-12b obtained using the HST WFC3 instrument. The disk-average dayside brightness temperature averages about 2900 K, peaking to 3200 K around 1.46 microns. We modeled a range of atmospheric cases for both the emission and transmission spectrum and confirm the recent finding by Crossfield et al. (2012b) that there is no evidence for C/O >1 in the atmosphere of WASP-12b. Assuming a physically plausible atmosphere, we find evidence that the presence of a number of molecules is consistent with the data, but the justification for inclusion of these opacity sources based on the Bayesian Information Criterion (BIC) is marginal. We also find the near-infrared primary eclipse light curve is consistent with small amounts of prolate distortion. As part of the calibration effort for these data, we conducted a detailed study of instrument systematics using 65 orbits of WFC3-IR grims observations. The instrument systematics are dominated by detector-related affects, which vary significantly depending on the detector readout mode. The 256x256 subarray observations of WASP 12 produced spectral measurements within 15% of the photon-noise limit using a simple calibration approach. Residual systematics are estimated to be less than 70 parts per million.

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.

{\it Ab initio} potential energy and dipole moment surfaces are computed for sulfur trioxide (SO$_3$) at the CCSD(T)-F12b level of theory and appropriate triple-zeta basis sets. These surfaces are fitted and used, with a slight correction for the equilibrium S--O distance, to compute pure rotational and rotation-vibraton spectra of \sothree\ using the variational nuclear motion program TROVE. The calculations considered transitions in the region 0--4000 \cm\ with rotational states up to J=85. The resulting line list of 349~348~513 transitions is appropriate for modelling room temperature SO3 spectra. Good agreement is found with the observed infrared absorption spectra and the calculations are used to place the measured relative intensities on an absolute scale. A list of 10~878 experimental transitions is provided in a form suitable for inclusion in standard atmospheric and planetary spectroscopic databases.

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$.

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$.

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.

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.

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.

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}.

A dedicated mission to investigate exoplanetary atmospheres represents a major milestone in our quest to understand our place in the universe by placing our Solar System in context and by addressing the suitability of planets for the presence of life. EChO -the Exoplanet Characterisation Observatory- is a mission concept specifically geared for this purpose. EChO will provide simultaneous, multi-wavelength spectroscopic observations on a stable platform that will allow very long exposures. EChO will build on observations by Hubble, Spitzer and groundbased telescopes, which discovered the first molecules and atoms in exoplanetary atmospheres. EChO will simultaneously observe a broad enough spectral region -from the visible to the mid-IR- to constrain from one single spectrum the temperature structure of the atmosphere and the abundances of the major molecular species. The spectral range and resolution are tailored to separate bands belonging to up to 30 molecules to retrieve the composition and temperature structure of planetary atmospheres. The target list for EChO includes planets ranging from Jupiter-sized with equilibrium temperatures Teq up to 2000 K, to those of a few Earth masses, with Teq ~300 K. We have baselined a dispersive spectrograph design covering continuously the 0.4-16 micron spectral range in 6 channels (1 in the VIS, 5 in the IR), which allows the spectral resolution to be adapted from several tens to several hundreds, depending on the target brightness. The instrument will be mounted behind a 1.5 m class telescope, passively cooled to 50 K, with the instrument structure and optics passively cooled to ~45 K. EChO will be placed in a grand halo orbit around L2. We have also undertaken a first-order cost and development plan analysis and find that EChO is easily compatible with the ESA M-class mission framework.

We present an analysis of seven primary transit observations of the hot Neptune GJ436b at 3.6, 4.5 and $8~\mu$m obtained with the Infrared Array Camera (IRAC) on the Spitzer Space Telescope. After correcting for systematic effects, we fitted the light curves using the Markov Chain Monte Carlo technique. Combining these new data with the EPOXI, HST and ground-based $V, I, H$ and $K_s$ published observations, the range $0.5-10~\mu$m can be covered. Due to the low level of activity of GJ436, the effect of starspots on the combination of transits at different epochs is negligible at the accuracy of the dataset. Representative climate models were calculated by using a three-dimensional, pseudo-spectral general circulation model with idealised thermal forcing. Simulated transit spectra of GJ436b were generated using line-by-line radiative transfer models including the opacities of the molecular species expected to be present in such a planetary atmosphere. A new, ab-initio calculated, linelist for hot ammonia has been used for the first time. The photometric data observed at multiple wavelengths can be interpreted with methane being the dominant absorption after molecular hydrogen, possibly with minor contributions from ammonia, water and other molecules. No clear evidence of carbon monoxide and dioxide is found from transit photometry. We discuss this result in the light of a recent paper where photochemical disequilibrium is hypothesised to interpret secondary transit photometric data. We show that the emission photometric data are not incompatible with the presence of abundant methane, but further spectroscopic data are desirable to confirm this scenario.

We present 'BYTe', a comprehensive 'hot' line list for the ro-vibrational transitions of ammonia, 14NH3, in its ground electronic state. This line list has been computed variationally using the program suite TROVE, a new spectroscopically-determined potential energy surface and an ab initio dipole moment surface. BYTe, is designed to be used at all temperatures up to 1500K. It comprises 1137650964 transitions in the frequency range from 0 to 12000 cm-1, constructed from 1366519 energy levels below 18000 cm-1 having J values below 36. Comparisons with laboratory data confirm the accuracy of the line list which is suitable for modelling a variety of astrophysical problems including the atmospheres of extrasolar planets and brown dwarfs.

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.