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

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

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

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.

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.