Ab initio modeling of molecular IR spectra of astrophysical interest: Application to CH4

Abstract

Aims. We describe an ab initio-based numerical method of obtaining infrared spectroscopic data (line list) of polyatomic molecules that allows calculation of complete sets of lines for temperatures up to several thousand Kelvin. While the main focus is on completeness and consistency, not spectroscopic accuracy, the approach is in principle "exact" for line positions and, although not exact for line strengths, of sufficient accuracy to be of value, especially in wavelength regions where there are gaps in reliable experimental data.Methods. Global potential energy and dipole moment hypersurfaces are fitted to the results of ab initio electronic structure calculations. The MULTIMODE software is then used to obtain rovibrational energy levels and dipole transition matrix elements. This information is used to calculate a complete set of Einstein coefficients of spontaneous emission Aij.Results. The method is applied to obtain a spectroscopic database for methane containing over 1.4 million lines up to an upper state energy of 6200 cm-1 (∼9000 K). The emission spectrum of CH4 at 1000 K is calculated with the complete set of Einstein coefficients and compared with the one obtained from the HITRAN database. Gaps in the database are realistically filled in by the calculated spectrum.Conclusions. Consistent and complete databases are important for astrophysical applications. Databases obtained by the method described here fulfill this requirement and are sufficiently accurate for astrophysical applications such as model atmosphere calculations and the corresponding synthetic spectra. © 2009 ESO.