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

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

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

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

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