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.

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.

We report follow-up observations of transiting exoplanets that have either large uncertainties (>10 minutes) in their transit times or have not been observed for over three years. A fully robotic ground-based telescope network, observations from citizen astronomers and data from TESS have been used to study eight planets, refining their ephemeris and orbital data. Such follow-up observations are key for ensuring accurate transit times for upcoming ground and space-based telescopes which may seek to characterise the atmospheres of these planets. We find deviations from the expected transit time for all planets, with transits occurring outside the 1 sigma uncertainties for seven planets. Using the newly acquired observations, we subsequently refine their periods and reduce the current predicted ephemeris uncertainties to 0.28 - 4.01 minutes. A significant portion of this work has been completed by students at two high schools in London as part of the Original Research By Young Twinkle Students (ORBYTS) programme.

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.

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.

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

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.

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.

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

Zirconium oxide(ZrO) is an important astrophysical molecule that defines the S-star classification class for cool giant stars. Accurate, empirical rovibronic energy levels, with associated labels and uncertainties, are reported for 9 low-lying electronic states of the diatomic 90Zr16O molecule. These 8088 empirical energy levels are determined using the Marvel (Measured Active Rotational-Vibrational Energy Levels) algorithm with 23 317 input assigned transition frequencies, 22 549 of which were validated. A temperature-dependent partition function is presented alongside updated spectroscopic constants for the 9 low-lying electronic states.

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.

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.

Cross sections are presented for dissociative recombination and electron-impact vibrational excitation of the ArH+ molecular ion at electron energies appropriate for the interstellar environment. The R-matrix method is employed to determine the molecular structure data, i.e. the position and width of the resonance states. The cross sections and the corresponding Maxwellian rate coefficients are computed using a method based on the Multichannel Quantum Defect Theory. The main result of the paper is the very low dissociative recombination rate found at temperatures below 1000K. This is in agreement with the previous upper limit measurement in merged beams and offers a realistic explanation to the presence of ArH+ in exotic interstellar conditions.

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.

Availability of reliable atomic and molecular opacity tables is essential in a wide variety of astronomical modeling: the solar and stellar interiors, stellar and planetary atmospheres, stellar evolution, pulsating stars, and protoplanetary disks, to name a few. With the advancement of powerful research techniques such as helio-seismology and asteroseismology, solar neutrino-flux measurements, exoplanet survey satellites, three-dimensional hydrodynamic atmospheric simulations (including non-LTE and granulation effects), high-performance computing of atomic and molecular data, and innovative plasma experiments the accuracy and completeness of opacity tables is being taken to an unprecedented level. The goal of the second Workshop on Astrophysical Opacities was to gather opacity data producers and consumers from both the atomic and molecular sectors to contribute to solving outstanding problems and to develop more effective and integrated interfaces. In this review we attempt to summa- rize the discussion at the workshop and propose future directions for opacity research.

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.

R-matrix calculations combined with the adiabatic-nuclei-rotation and Coulomb-Born approximations are used to compute electron-impact rotational rate coefficients for two open-shell diatomic cations of astrophysical interest: the hydoxyl and sulphanyl ions, OH$^+$ and SH$^+$. Hyperfine resolved rate coefficients are deduced using the infinite-order-sudden approximation. The propensity rule $\Delta F=\Delta j=\Delta N=\pm 1$ is observed, as is expected for cations with a large dipole moment. A model for OH$^+$ excitation in the Orion Bar photon-dominated region (PDR) is presented which nicely reproduces Herschel observations for an electron fraction $x_e=10^{-4}$ and an OH$^+$ column density of $3\times 10^{13}$~cm$^{-2}$. Electron impact electronic excitation cross sections and rate coefficients for the ions are also presented.

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.

We present a detailed theoretical study of the rotational excitation of CH$^+$ due to reactive and nonreactive collisions involving C$^+(^2P)$, H$_2$, CH$^+$, H and free electrons. Specifically, the formation of CH$^+$ proceeds through the reaction between C$^+(^2P)$ and H$_2(\nu_{\rm H_2}=1, 2)$, while the collisional (de)excitation and destruction of CH$^+$ is due to collisions with hydrogen atoms and free electrons. State-to-state and initial-state-specific rate coefficients are computed in the kinetic temperature range 10-3000~K for the inelastic, exchange, abstraction and dissociative recombination processes using accurate potential energy surfaces and the best scattering methods. Good agreement, within a factor of 2, is found between the experimental and theoretical thermal rate coefficients, except for the reaction of CH$^+$ with H atoms at kinetic temperatures below 50~K. The full set of collisional and chemical data are then implemented in a radiative transfer model. Our Non-LTE calculations confirm that the formation pumping due to vibrationally excited H$_2$ has a substantial effect on the excitation of CH$^+$ in photon-dominated regions. In addition, we are able to reproduce, within error bars, the far-infrared observations of CH$^+$ toward the Orion Bar and the planetary nebula NGC~7027. Our results further suggest that the population of $\nu_{\rm H_2}=2$ might be significant in the photon-dominated region of NGC~7027.

Accurate, experimental rovibronic energy levels, with associated labels and uncertainties, are reported for 11 low-lying electronic states of the diatomic \TiO\ molecule, determined using the \Marvel\ (Measured Active Rotational-Vibrational Energy Levels) algorithm. All levels are based on lines corresponding to critically reviewed and validated high-resolution experimental spectra taken from 24 literature sources. The transition data are in the 2 $-$ 22,160 \cm{} region. Out of the 49,679 measured transitions, 43,885 are triplet-triplet, 5710 are singlet-singlet and 84 are triplet-singlet transitions. A careful analysis of the resulting experimental spectroscopic network (SN) allows 48,590 transitions to be validated. The transitions determine 93 vibrational band origins of \TiO\, including 71 triplet and 22 singlet ones. There are 276 (73) triplet-triplet (singlet-singlet) band-heads derived from \Marvel\ experimental energies, 123 (38) of which have never been assigned in low or high resolution experiments. The highest $J$ value, where $J$ stands for the total angular momentum, for which an energy level is validated is 163. The number of experimentally-derived triplet and singlet \TiO\ rovibrational energy levels is 8682 and 1882, respectively. The lists of validated lines and levels for \TiO\ are deposited in the Supporting Information to this paper.

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

Remote sensing experiments require high-accuracy, preferably sub-percent, line intensities and in response to this need we present computed room temperature line lists for six symmetric isotopologues of carbon dioxide: $^{13}$C$^{16}$O$_2$, $^{14}$C$^{16}$O$_2$, $^{12}$C$^{17}$O$_2$, $^{12}$C$^{18}$O$_2$, $^{13}$C$^{17}$O$_2$ and $^{13}$C$^{18}$O$_2$, covering the range 0-8000 \cm. Our calculation scheme is based on variational nuclear motion calculations and on a reliability analysis of the generated line intensities. Rotation-vibration wavefunctions and energy levels are computed using the DVR3D software suite and a high quality semi-empirical potential energy surface (PES), followed by computation of intensities using an \abinitio\ dipole moment surface (DMS). Four line lists are computed for each isotopologue to quantify sensitivity to minor distortions of the PES/DMS. Reliable lines are benchmarked against recent state-of-the-art measurements and against the HITRAN2012 database, supporting the claim that the majority of line intensities for strong bands are predicted with sub-percent accuracy. Accurate line positions are generated using an effective Hamiltonian. We recommend the use of these line lists for future remote sensing studies and their inclusion in databases.

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

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

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

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

{\bf R}-matrix calculations combined with the adiabatic nuclei approximation are used to compute electron-impact rotiational excitation rates for three closed-shell diatomic cations, HeH$^+$, CH$^+$, ArH$^+$. Comparisons with previous studies show that an improved treatment of threshold effects leads to significant changes in the low temperature rates, furthermore the new calculations suggest that excitation of CH$^+$ is dominated by $\Delta J =1$ transitions as is expected for cations with a large dipole moment. A model for ArH$^+$ excitation in the Crab Nebula is presented which gives results consistent with the observations for electron densities in the range $2-3\times 10^3$~cm$^{-3}$.

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

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

Chemical models used to study the chemical composition of the gas and the ices in the interstellar medium are based on a network of chemical reactions and associated rate coefficients. These reactions and rate coefficients are partially compiled from data in the literature, when available. We present in this paper kida.uva.2014, a new updated version of the kida.uva public gas-phase network first released in 2012. In addition to a description of the many specific updates, we illustrate changes in the predicted abundances of molecules for cold dense cloud conditions as compared with the results of the previous version of our network, kida.uva.2011.

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 quantum electrodynamics (QED) correction surface for the simplest polyatomic and polyelectronic system H$_3^+$ is computed using an approximate procedure. This surface is used to calculate the shifts to vibration-rotation energy levels due to QED; such shifts have a magnitude of up to 0.25 cm$^{-1}$ for vibrational levels up to 15~000 cm$^{-1}$ and are expected to have an accuracy of about 0.02 cm$^{-1}$. Combining the new H$_3^+$ QED correction surface with existing highly accurate Born-Oppenheimer (BO), relativistic and adiabatic components suggests that deviations of the resulting {\it ab initio} energy levels from observed ones are largely due to non-adiabatic effects.

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

Utilising previous work by the authors on the spin-coupled rotational cross-sections for electron-CN collisions, data for the associated rate coefficients is presented. Data on rotational, fine-structure and hyperfine-structure transition involving rotational levels up to $N$=20 are computed for temperatures in the range 10 -- 1000~K. Rates are calculated by combining Born-corrected R-matrix calculations with the infinite-order-sudden (IOS) approximation. The dominant hyperfine transitions are those with $\Delta N=\Delta j= \Delta F=1$. For dipole-allowed transitions, electron-impact rates are shown to exceed those for excitation of CN by para-H$_2$($j=0$) by five orders of magnitude. The role of electron collisions in the excitation of CN in diffuse clouds, where local excitation competes with the cosmic microwave background (CMB) photons, is considered. Radiative transfer calculations are performed and the results compared to observations. These comparisons suggest that electron density lies in the range $n(e)\sim 0.01-0.06$~cm$^{-3}$ for typical physical conditions present in diffuse clouds.

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

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

We present a novel chemical database for gas-phase astrochemistry. Named the KInetic Database for Astrochemistry (KIDA), this database consists of gas-phase reactions with rate coefficients and uncertainties that will be vetted to the greatest extent possible. Submissions of measured and calculated rate coefficients are welcome, and will be studied by experts before inclusion into the database. Besides providing kinetic information for the interstellar medium, KIDA is planned to contain such data for planetary atmospheres and for circumstellar envelopes. Each year, a subset of the reactions in the database (kida.uva) will be provided as a network for the simulation of the chemistry of dense interstellar clouds with temperatures between 10 K and 300 K. We also provide a code, named Nahoon, to study the time-dependent gas-phase chemistry of 0D and 1D interstellar sources.

Energy exchange processes play a crucial role in the early Universe, affecting the thermal balance and the dynamical evolution of the primordial gas. In the present work we focus on the consequences of a non-thermal distribution of the level populations of H$_2$: first, we determine the excitation temperatures of vibrational transitions and the non-equilibrium heat transfer; second, we compare the modifications to chemical reaction rate coefficients with respect to the values obtained assuming local thermodynamic equilibrium; third, we compute the spectral distortions to the cosmic background radiation generated by the formation of H$_2$ in vibrationally excited levels. We conclude that non-equilibrium processes cannot be ignored in cosmological simulations of the evolution of baryons, although their observational signatures remain below current limits of detection. New fits to the equilibrium and non-equilibrium heat transfer functions are provided.

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.

Exoplanetary science is among the fastest evolving fields of today's astronomical research. Ground-based planet-hunting surveys alongside dedicated space missions (Kepler, CoRoT) are delivering an ever-increasing number of exoplanets, now numbering at ~690, with ESA's GAIA mission planned to bring this number into the thousands. The next logical step is the characterisation of these worlds: what is their nature? Why are they as they are? The use of the HST and Spitzer Space Telescope to probe the atmospheres of transiting hot, gaseous exoplanets has demonstrated that it is possible with current technology to address this ambitious goal. The measurements have also shown the difficulty of understanding the physics and chemistry of these environments when having to rely on a limited number of observations performed on a handful of objects. To progress substantially in this field, a dedicated facility for exoplanet characterization with an optimised instrument design (detector performance, photometric stability, etc.), able to observe through time and over a broad spectral range a statistically significant number of planets, will be essential. We analyse the performances of a 1.2/1.4m space telescope for exoplanet transit spectroscopy from the visible to the mid IR, and present the SNR ratio as function of integration time and stellar magnitude/spectral type for the acquisition of spectra of planetary atmospheres in a variety of scenarios: hot, warm, and temperate planets, orbiting stars ranging in spectral type from hot F to cool M dwarfs. We include key examples of known planets (e.g. HD 189733b, Cancri 55 e) and simulations of plausible terrestrial and gaseous planets, with a variety of thermodynamical conditions. We conclude that even most challenging targets, such as super-Earths in the habitable-zone of late-type stars, are within reach of a M-class, space-based spectroscopy mission.

Cooling of primordial gas plays a crucial role in the birth of the first structures in our Universe. Due to the low fractional abundance of molecular species at high redshifts, spontaneous emission rather than collisions represents the most efficient way to cool the pristine plasma. In the present work, radiative cooling functions are evaluated for the diatomic species HD, HD$^+$, HeH$^+$, LiH and LiH$^+$. Cooling functions for the triatomic ions H$_3^+$ and H$_2$D$^+$ are also considered. Analytic fits as functions of temperature are provided.

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.

Shock modelling predicts an electron density enhancement within the magnetic precursor of C-shocks. Previous observations of SiO, H13CO+, HN13C and H13CN toward the young L1448-mm outflow showed an over-excitation of the ion fluid that was attributed to an electron density enhancement in the precursor. We re-visit this interpretation and test if it still holds when we consider different source morphologies and kinetic temperatures for the observed molecules, and also give some insight on the spatial extent of the electron density enhancement around L1448-mm. We estimate the opacities of H13CO+ and HN13C by observing the J=3\to2 lines of rarer isotopologues to confirm that the emission is optically thin. To model the excitation of the molecules, we use the large velocity gradient (LVG) approximation with updated collisional coefficients to i) re- analyse the observations toward the positions where the over-excitation of H13CO+ has previously been observed [i.e. toward L1448- mm at offsets (0,0) and (0,-10)], and ii) to investigate if the electron density enhancement is still required for the cases of extended and compact emission, and for kinetic temperatures of up to 400 K. We also report several lines of SiO, HN13C and H13CO+ toward new positions around this outflow, to investigate the spatial extent of the over-excitation of the ions in L1448-mm. From the isotopologue observations, we find that the emission of H13CO+ and HN13C from the precursor is optically thin if this emission is extended. Using the new collisional coefficients, an electron density enhancement is still needed to explain the excitation of H13CO+ for extended emission and for gas temperatures of\le 400 K toward L1448-mm (0,-10), and possibly also toward L1448-mm (0,0). For compact emission the data cannot be fitted. We do not find any evidence for the over-excitation of the ion fluid toward the newly observed positions around L1448-mm. The observed line emission of SiO, H13CO+ and HN13C toward L1448-mm (0,0) and (0,-10) is consistent with an electron density enhancement in the precursor component, if this emission is spatially extended. This is also true for the case of high gas temperatures (\le400 K) toward the (0,-10) offset. The electron density enhancement seems to be restricted to the southern, redshifted lobe of the L1448-mm outflow. Interferometric images of the line emission of these molecules are needed to confirm the spatial extent of the over-excitation of the ions and thus, of the electron density enhancement in the magnetic precursor of L1448-mm.

The hot Jupiter HD 209458b was observed during primary transit at 3.6, 4.5, 5.8 and 8.0 microns using the Infrared Array Camera (IRAC) on the Spitzer Space Telescope. We detail here the procedures we adopted to correct for the systematic trends present in the IRAC data. The light curves were fitted including limb darkening effects and fitted using Markov Chain Monte Carlo and prayer-bead Monte Carlo techniques, finding almost identical results. The final depth measurements obtained by a combined Markov Chain Monte Carlo fit are at 3.6 microns, 1.469 +- 0.013 % and 1.448 +- 0.013 %; at 4.5 microns, 1.478 +- 0.017 % ; at 5.8 microns, 1.549 +- 0.015 % and at 8.0 microns 1.535 +- 0.011 %. Our results clearly indicate the presence of water in the planetary atmosphere. Our broad band photometric measurements with IRAC prevent us from determining the additional presence of other other molecules such as CO, CO2 and methane for which spectroscopy is needed. While water vapour with a mixing ratio of 10^-4-10^-3 combined with thermal profiles retrieved from the day-side may provide a very good fit to our observations, this data set alone is unable to resolve completely the degeneracy between water abundance and atmospheric thermal profile.

High resolution spectra of Comet 8P/Tuttle were obtained in the frequency range 3440.6-3462.6 cm-1 on 3 January 2008 UT using CGS4 with echelle grating on UKIRT. In addition to recording strong solar pumped fluorescent (SPF) lines of H2O, the long integration time (152 miutes on target) enabled eight weaker H2O features to be assigned, most of which had not previously been identified in cometary spectra. These transitions, which are from higher energy upper states, are similar in character to the so-called 'SH' lines recorded in the post Deep Impact spectrum of comet Tempel 1 (Barber et al., 2007). We have identified certain characteristics that these lines have in common, and which in addition to helping to define this new class of cometary line, give some clues to the physical processes involved in their production. Finally, we derive an H2O rotational temperature of 62+/- K and a water production rate of (1.4+/-0.3)E28 molecules/s.

The molecular gas composition in the inner 1 kpc disk of the starburst galaxy M82 resembles that of Galactic Photon Dominated Regions (PDRs). In particular, large abundances of the reactive ions HOC+ and CO+ have been measured in the nucleus of this galaxy. To investigate the origin of the large abundances of reactive ions in M82, we have completed our previous 30m HOC+ J=1-0 observations with the higher excitation HCO+ and HOC+ J=4-3 and 3-2 rotational lines. In addition, we have obtained with the IRAM Plateau de Bure Interferometer (PdBI) a 4" resolution map of the HOC+ 1-0 emission, the first ever obtained in a Galactic or extragalactic source. Our HOC+ interferometric image shows that the emission of the HOC+ 1-0 line is mainly restricted to the nuclear disk, with the maxima towards the E. and W. molecular peaks. In addition, line excitation calculations imply that the HOC+ emission arises in dense gas. Therefore, the HOC+ emission is arising in the dense PDRs embedded in the M82 nuclear disk, rather than in the intercloud phase and/or wind. We have improved our previous chemical model of M82 by (i) using the new version of the Meudon PDR code, (ii) updating the chemical network, and (iii) considering two different types of clouds (with different thickness). Most molecular observations (HCO+, HOC+, CO+, CN, HCN, H3O+) are well explained assuming that ~ 87 % of the mass of the molecular gas is forming small clouds (Av=5 mag) while only ~ 13 % of the mass is in large molecular clouds (Av=50 mag). Such small number of large molecular clouds suggests that M82 is an old starburst, where star formation has almost exhausted the molecular gas reservoir.

We present the 3.5m SPICA space telescope, a proposed Japanese-led JAXA-ESA mission scheduled for launch around 2017. The actively cooled (<5 K), single aperture telescope and monolithic mirror will operate from ~3.5 to ~210 um and will provide superb sensitivity in the mid- and far-IR spectral domain (better than JWST at lambda > 18 um). SPICA is one of the few space missions selected to go to the next stage of ESA's Cosmic Vision 2015-2025 selection process. In this White Paper we present the main specifications of the three instruments currently baselined for SPICA: a mid-infrared (MIR) coronagraph (~3.5 to ~27 um) with photometric and spectral capabilities (R~200), a MIR wide-field camera and high resolution spectrometer (R~30,000), and a far-infrared (FIR ~30 to ~210 um) imaging spectrometer - SAFARI - led by a European consortium. We discuss their capabilities in the context of MIR direct observations of exo-planets (EPs) and multiband photometry/high resolution spectroscopy observations of transiting exo-planets. We conclude that SPICA will be able to characterize the atmospheres of transiting exo-planets down to the super-Earth size previously detected by ground- or space-based observatories. It will also directly detect and characterize Jupiter/Neptune-size planets orbiting at larger separation from their parent star (>5-10 AU), by performing quantitative atmospheric spectroscopy and studying proto-planetary and debris disks. In addition, SPICA will be a scientific and technological precursor for future, more ambitious, IR space missions for exo-planet direct detection as it will, for example, quantify the prevalence exo-zodiacal clouds in planetary systems and test coronographic techniques, cryogenic systems and lightweight, high quality telescopes. (abridged)

A line list of vibration-rotation transitions for 13C substituted HCN is presented. The line list is constructed using known experimental levels where available, calculated levels and ab initio line intensities originally calculated for the major isotopologue. Synthetic spectra are generated and compared with observations for cool carbon star WZ Cas. It is suggested that high resolution HCN spectra recorded near 14 micron should be particularly sensitive to the 13C to 12C ratio.

We present high resolution (R = 90,000) mid-infrared spectra of M dwarfs. The mid infrared region of the spectra of cool low mass stars contain pure rotational water vapour transitions that may provide us with a new methodology in the determination of the effective temperatures for low mass stars. We identify and assign water transitions in these spectra and determine how sensitive each pure rotational water transition is to small (25 K) changes in effective temperature. We find that, of the 36 confirmed and assigned pure rotational water transitions, at least 10 should be sensitive enough to be used as temperature indicators.

We compute opacities for the electronic molecular band systems A 6Sigma+ -- X 6Sigma+ of CrH and CrD, and A 2Pi -- X 2Sigma+ of MgH and MgD. The opacities are computed by making use of existing spectroscopic constants for MgH and CrH. These constants are adjusted for the different reduced masses of MgD and CrD. Frank-Condon factors are used to provide intensities for the individual vibronic bands. These results are used in the computation of synthetic spectra between Tef = 1800 and 1200 K with an emphasis on the realisation of ``deuterium test'', first proposed by Bejar et al. (1999) to distinguish brown dwarfs from planetary mass objects. We discuss the possible use of CrD and MgD electronic bands for the "deuterium test". We find CrD to be the more promising of the two deuterides, potentially, the most useful bands of CrH/CrD are the Delta v = +1 and Delta v = -1 at 0.795 and 0.968 micron.

Rotational excitation of isotopologues of HCN and HNC by thermal electron-impact is studied using the molecular {\bf R}-matrix method combined with the adiabatic-nuclei-rotation (ANR) approximation. Rate coefficients are obtained for electron temperatures in the range 5$-$6000 K and for transitions among all levels up to J=8. Hyperfine rates are also derived using the infinite-order-sudden (IOS) scaling method. It is shown that the dominant rotational transitions are dipole allowed, that is those for which $\Delta J=1$. The hyperfine propensity rule $\Delta J=\Delta F$ is found to be stronger than in the case of He$-$HCN collisions. For dipole allowed transitions, electron-impact rates are shown to exceed those for excitation of HCN by He atoms by 6 orders of magnitude. As a result, the present rates should be included in any detailed population model of isotopologues of HCN and HNC in sources where the electron fraction is larger than 10$^{-6}$, for example in interstellar shocks and comets.

We present an analysis of high resolution spectra in the J band of five ultra cool dwarfs from M6 to L0. A new ab initio water vapour line list and existing line lists of FeH and CrH were used for spectra modelling. We find a good fit for the Mn I 12899.76 A line. This feature is one of the few for which we have a reliable oscillator strength. Other atomic features are present but most of the observed features are FeH and water lines. While we are uncertain about the quality of many of the atomic line parameters, the FeH and CrH line lists predict a number of features which are not apparent in our observed spectra. We infer that the main limiting factor in our spectral analysis is the FeH and CrH molecular spectra.

Water is predicted to be among, if not the most abundant molecular species after hydrogen in the atmospheres of close-in extrasolar giant planets (hot-Jupiters) Several attempts have been made to detect water on an exoplanet, but have failed to find compelling evidence for it or led to claims that should be taken with caution. Here we report an analysis of recent observations of the hot-Jupiter HD189733b taken during the transit, where the planet passed in front of its parent star. We find that absorption by water vapour is the most likely cause of the wavelength-dependent variations in the effective radius of the planet at the infrared wavelengths 3.6, 5.8 and 8 microns. The larger effective radius observed at visible wavelengths may be due to either star variability or the presence of clouds/hazes. We explain the most recent thermal infrared observations of the planet during secondary transit behind the star, reporting a non-detection of water on HD189733b, as being a consequence of the nearly isothermal vertical profile of the planet.s atmosphere. Our results show that water is detectable on extrasolar planets using the primary transit technique and that the infrared should be a better wavelength region than the visible, for such searches.

The recently puplished Barber-Tennyson (BT2) synthetic H$_2$O water line list is the most complete and accurate line list in existence. It is finding application in a wide range of astrophysical environments. UKIRT spectra of comet Tempel 1, obtained after the 'Deep Impact' event, revealed several known H$_2$O solar pumped fluorescent (SPF)lines in the 2.8945 to 2.8985 $\mu$m region. In addition, using synthetic spectra produced with BT2, several emission lines were identified that had not previously been recorded in cometary spectra. Unlike the SPF lines, which are transitions from doubly-excited stretch states, these transitions, that we label 'SH', are from states with three or four quanta of vibrational excitation. The SH features were particularly strong during the period 20-40 minutes after impact.

The primordial elements, H, He and Li are included in a low temperature equation of state and monochromatic opacity function. The equation of state and opacity function are incorporated into the stellar evolution code NG-ELMS, which makes use of a non-grey model atmosphere computed at runtime. NG-ELMS is used to compute stellar evolution models for primordial and lithium free element mixtures, for stars in the sub-Solar mass range 0.8--0.15 Msol. We find that lithium has little or no effect upon the structure and observable properties of of stars in this mass range. Furthermore lithium is completely destroyed by fusion before the main sequence in stars of mass less than ~0.7 Mol. We find that on the red giant branch and Hayashi track, the use of a non-grey model atmosphere to provide the upper boundary conditions for the stellar evolution calculation, results in significantly cooler less luminous stars, across the mass range

Magnetic precursors of C-shocks accelerate, compress and heat molecular ions, modifying the kinematics and the physical conditions of the ion fluid with respect to the neutral one. Electron densities are also expected to be significantly enhanced in shock precursors. In this Letter, we present observations of strongly polar ion and neutral molecules such as SiO, H13CO+, HN13C and H13CN, which reveal the electron density enhancements associated with the precursor of the young L1448-mm outflow. While in the ambient gas the excitation of the ions and neutrals is explained by collisional excitation by H2 with a single density of 10E5cmE-3, H13CO+ shows an over excitation in the shock precursor component that requires H2 densities of a factor of >10 larger than those derived from the neutral species. This over excitation in H13CO+ can be explained if we consider an additional excitation by collisions with electrons and an electron density enhancement in the precursor stage by a factor of 500, i.e. a fractional ionization of 5x10E-5. These results show that multiline observations can be used to study the evolution of the ion and electron fluids at the first stages of the C-shock interaction.

Differential cross sections for electron collisions with the O$_2$ molecule in its ground ${X}^{3}\Sigma_g^-$ state, as well as excited ${a}^{1}\Delta_g$ and ${b}^{1}\Sigma_g^+$ states are calculated. As previously, the fixed-bond R-matrix method based on state-averaged complete active space SCF orbitals is employed. In additions to elastic scattering of electron with the O$_2$ ${X}^{3}\Sigma_g^-$, ${a}^{1}\Delta_g$ and ${b}^{1}\Sigma_g^+$ states, electron impact excitation from the ${X}^{3}\Sigma_g^-$ state to the ${a}^{1}\Delta_g$ and ${b}^{1}\Sigma_g^+$ states as well as '6 eV states' of ${c}^{1}\Sigma_u^{-}$, ${A'}^{3}\Delta_u$ and ${A}^{3}\Sigma_u^{+}$ states is studied. Differential cross sections for excitation to the '6 eV states' have not been calculated previously. Electron impact excitation to the ${b}^{1}\Sigma_g^+$ state from the metastable ${a}^{1}\Delta_g$ state is also studied. For electron impact excitation from the O$_2$ ${X}^{3}\Sigma_g^-$ state to the ${b}^{1}\Sigma_g^+$ state, our results agree better with the experimental measurements than previous theoretical calculations. Our cross sections show angular behaviour similar to the experimental ones for transitions from the ${X}^{3}\Sigma_g^-$ state to the '6 eV states', although the calculated cross sections are up to a factor two larger at large scattering angles. For the excitation from the ${a}^{1}\Delta_g$ state to the ${b}^{1}\Sigma_g^+$ state, our results marginally agree with the experimental data except for the forward scattering direction.

Low-energy electron collisions with O$_2$ molecules are studied using the fixed-bond R-matrix method. In addition to the O$_2$ ${X}^3\Sigma_{g}^-$ ground state, integrated cross sections are calculated for elecron collisions with the ${a}^1\Delta_{g}$ and ${b}^1\Sigma_{g}^+$ excited states of O$_2$ molecules. 13 target electronic states of O$_2$ are included in the model within a valence configuration interaction representations of the target states. Elastic cross sections for the ${a}^1\Delta_{g}$ and ${b}^1\Sigma_{g}^+$ excited states are similar to the cross sections for the ${X}^3\Sigma_{g}^-$ ground state. As in case of excitation from the ${X}^3\Sigma_{g}^-$ state, the O$_2^-$ $\Pi_u$ resonance makes the dominant contribution to excitation cross sections from the ${a}^1\Delta_{g}$ and ${b}^1\Sigma_{g}^+$ states. The magnitude of excitation cross sections from the ${a}^1\Delta_{g}$ state to the ${b}^1\Sigma_{g}^+$ state is about 10 time larger than the corresponding cross sections from the ${X}^3\Sigma_{g}^-$ to the ${b}^1\Sigma_{g}^+$ state. For this ${a}^1\Delta_{g}$ $\to$ ${b}^1\Sigma_{g}^+$ transition, our cross section at 4.5 eV agrees well with the available experimental value. These results should be important for models of plasma discharge chemistry which often requires cross sections between the excited electronic states of O$_2$.

A computed list of H$_{2}$$^{16}$O infra-red transition frequencies and intensities is presented. The list, BT2, was produced using a discrete variable representation two-step approach for solving the rotation-vibration nuclear motions. It is the most complete water line list in existence, comprising over 500 million transitions (65% more than any other list) and it is also the most accurate (over 90% of all known experimental energy levels are within 0.3 cm$^{-1}$ of the BT2 values). Its accuracy has been confirmed by extensive testing against astronomical and laboratory data. The line list has been used to identify individual water lines in a variety of objects including: comets, sunspots, a brown dwarf and the nova-like object V838 Mon. Comparison of the observed intensities with those generated by BT2 enables physical values to be derived for these objects. The line list can also be used to provide an opacity for models of the atmospheres of M-dwarf stars and assign previously unknown water lines in laboratory spectra.

We build an accurate database of 5200 HCN and HNC rotation-vibration energy levels, determined from existing laboratory data. 20~000 energy levels in the Harris et al. (2002) linelist are assigned approximate quantum numbers. These assignments, lab determined energy levels and Harris et al (2002) energy levels are incorporated in to a new energy level list. A new linelist is presented, in which frequencies are computed using the lab determined energy levels where available, and the ab initio energy levels otherwise. The new linelist is then used to compute new model atmospheres and synthetic spectra for the carbon star WZ Cas. This results in better fit to the spectrum of WZ Cas in which the absorption feature at 3.56 micron is reproduced to a higher degree of accuracy than has previously been possible. We improve the reproduction of HCN absorption features by reducing the abundance of Si to [Si/H] = --0.5 dex, however, the strengths of the $\Delta v=2$ CS band heads are over-predicted.

We present results of modelling the bulk of the spectral energy distribution (0.35 - 5 micron) for GJ406 (M6V). Synthetic spectra were calculated using the NextGen, Dusty and Cond model atmospheres and incorporate line lists for H2O, TiO, CrH, FeH, CO, MgH molecules as well as the VALD line list of atomic lines. A comparison of synthetic and observed spectra gives Tef = 2800 +/- 100 K. We determine M$_bol = 12.13 +/- 0.10 for which evolutionary models by Baraffe et al. (2003) suggest an age of around 0.1 -- 0.35 Gyr consistent with its high activity. The age and luminosity of GJ406 correspond to a wide range of plausible masses (0.07 -- 0.1 Msun).

V838 Monocerotis had an intriguing, nova-like outburst in January 2002 which has subsequently led to several studies of the object. It is now recognized that the outburst of V838 Mon and its evolution are different from that of a classical nova or other classes of well-known eruptive variables. V838 Mon, along with two other objects that have analogous properties, appears to comprise a new class of eruptive variables. There are limited infrared studies of V838 Mon. Here, we present near-infrared H band (1.5 - 1.75micron) spectra of V838 Mon from late 2002 to the end of 2004. The principal, new result from our work is the detection of several, rotation-vibration lines of water in the H band spectra. The observed water lines have been modeled to first establish that they are indeed due to water. Subsequently the temperature and column densities of the absorbing material, from where the water absorption features originate, are derived. From our analysis, we find that the water features arise from a cool ~750-900 K region around V838 Mon which appears to be gradually cooling with time.

We compare high resolution infrared observations of the CO 3-1 bands in the 2.297-2.310 micron region of M dwarfs and one L dwarf with theoretical expectations. We find a good match between the observational and synthetic spectra throughout the 2000-3500K temperature regime investigated. Nonetheless, for the 2500-3500 K temperature range the temperatures that we derive from synthetic spectral fits are higher than expected from more empirical methods by several hundred K. In order to reconcile our findings with the empirical temperature scale it is necessary to invoke warming of the model atmosphere used to construct the synthetic spectra. We consider that the most likely reason for the back-warming is missing high temperature opacity due to water vapour. We compare the water vapour opacity of the Partridge & Schwenke (1997) line list used for the model atmosphere with the output from a preliminary calculation by Barber & Tennyson (2004). While the Partridge & Schwenke line list is a reasonable spectroscopic match for the new line list at 2000 K, by 4000 K it is missing around 25% of the water vapour opacity. We thus consider that the offset between empirical and synthetic temperature scales is explained by the lack of hot water vapour used for computation of the synthetic spectra. For our coolest objects with temperatures below 2500 K we find best fits when using synthetic spectra which include dust emission. Our spectra also allow us to constrain the rotational velocities of our sources, and these velocities are consistent with the broad trend of rotational velocities increasing from M to L.

HeH$^+$ is found to be the dominant positive ion over a wide range of temperatures and densities relevant to helium rich white dwarfs. The inclusion of HeH$^+$ in ionization equilibrium computations increases the abundance of free electrons by a significant factor. For temperatures below 8000 K, He$^-$ free-free absorption is increased by up to a factor of 5, by the inclusion of HeH$^+$. Illustrative model atmospheres and spectral energy distributions are computed, which show that HeH$^+$ has a strong effect upon the density and pressure structure of helium rich white dwarfs with teff < 8000 K. The inclusion of HeH$^+$ significantly reddens spectral energy distributions and broad band color indices for models with Teff < 5500 K. This has serious implications for existing model atmospheres, synthetic spectra and cooling curves for helium rich white dwarfs.

The wavelength and Einstein A coefficient are calculated for all rotation-vibration transitions of $^4$He$^1$H$^+$, $^3$He$^1$H$^+$, $^4$He$^2$H$^+$ and $^3$He$^2$H$^+$, giving a complete line list and the partition function for $^4$HeH$^+$ and its isotopologues. This opacity is included in the calculation of the total opacity of low-metallicity stars and its effect is analysed for different conditions of temperature, density and hydrogen number fraction. For a low helium number fraction (as in the Sun), it is found that HeH$^+$ has a visible but small effect for very low densities ($\rho\leq 10^{-10}$g cm$^{-3}$), at temperatures around 3500 K. However, for high helium number fraction, the effect of HeH$^+$ becomes important for higher densities ($\rho\leq 10^{-6}$g cm$^{-3}$), its effect being most important for a temperature around 3500 K. Synthetic spectra for a variety of different conditions are presented.

H3+ and the deuterated isotopomers are thought to play an important role in interstellar chemistry. The partition functions of H3+, D2H+ and D3+ are calculated to a temperature of 800 K by explicitly summing the ab initio determined rotation-vibration energy levels of the respective species. These partition functions are used to calculate the equilibrium constants for nine important reactions in the interstellar medium involving H3+ and its deuterated isotopomers. These equilibrium constants are compared to previously determined experimental and theoretical values. The Einstein A coefficients for the strongest dipole transitions are also calculated.

H3+ has been shown (Lenzuni et al. 1991, ApJS, 76, 759) to be the dominant positive ion, in a zero metallicity gas at low temperature and intermediate to high density. It therefore affects both the number of free electrons and the opacity of the gas. The most recent H3+ partition function (Neale and Tennyson, 1995, ApJ, 454, L169) is an order of magnitude larger at 4000 K than all previous partition functions, implying that H3+ is a more important electron donor than previously thought. Here we present new Rosseland mean opacities for a hydrogen-helium gas of 1000 < T(K) < 9000 and -14 < log10(density [g/cc]) < -2. In the calculation of these opacities we have made use of the latest collision induced absorption data as well as the most recent H3+ partition function and line opacity data. It is shown that these updated and new sources of opacity give rise to a Rosseland mean opacity for a hydrogen-helium gas which is in general greater than that calculated in earlier works. The new opacity data are then used to model the evolution of low mass 0.15-0.8 Mo zero metallicity stars, from pre-Main Sequence collapse to Main Sequence turn-off. To investigate the effect of H3+ on the evolution of low mass zero metallicity stars, we repeat our calculations neglecting H3+ as a source of electrons and line opacity. We find that H3+ can have an effect on the structure and evolution of stars of mass ~0.5 Mo or less. A gray atmosphere is used for the calculation, which is sufficient to demonstrate that H3+ affects the evolution of very low mass stars to a greater degree than previously believed.

Model carbon star atmospheres and synthetic spectra have been calculated using the recent HCN/HNC vibration rotation linelist of Harris et al. (2002) ApJ, 578, 657. The calculations are repeated using only HCN lines and show that HNC has a significant effect upon the temperature, density and optical depth of a stellar atmosphere. We fit synthetic spectra in the 2.7-4.0 micron region to observed ISO spectra of the carbon stars WZ Cas and TX Psc obtained by Aoki et al. (1998), A&A, 340, 222. These fits allow us to identify absorption by HNC in the spectrum of WZ Cas at 2.8-2.9 microns, and to determine new independent estimates of effective temperature and log(Nc)/log(No). The findings reported here indicate that absorption by both HCN and HNC is needed to fully explain the observed stellar spectra and represent the first identification of HNC in a star. Q branch absorption by the HCN $\Delta v_2=1$, $\Delta v_3=1$ and $\Delta v_1=1$, $\Delta v_2=-1$ bands at 3.55 and 3.86 microns respectively, are identified in the spectrum of WZ Cas.

M star spectra, at wavelengths beyond 1.35 $\mu$m, are dominated by water vapour yet terrestrial water vapour makes it notoriously difficult to make accurate measurement from ground-based observations. We have used the short wavelength spectrometer on the Infrared Space Observatory at four wavelength settings to cover the 2.5--3.0 $\mu$m region for a range of M stars. The observations show a good match with previous ground-based observations and with synthetic spectra based on the Partridge & Schwenke (1997) line list though not with the SCAN (Jorgensen et al. 2001) line list. We used a least-squared minimisation technique to systematically find best fit parameters for the sample of stars. The temperatures that we find indicate a relatively hot temperature scale for M dwarfs. We consider that this could be a consequence of problems with the Partridge & Schwenke linelist which leads to synthetic spectra predicting water bands which are too strong for a given temperature. Such problems need to be solved in the next generation of water vapour line lists which will extend the calculation of water vapour to higher energy levels with the good convergence necessary for reliable modelling of hot water vapour. Then water bands can assume their natural role as the primary tool for the spectroscopic analysis of M stars.