Modelling the spectral energy distribution of ULIRGs^*

II. The energetic environment and the dense interstellar medium

¹ INAOE, Luis Enrique Erro 1, 72840 Tonantzintla, Puebla, Mexico
² INAF - Osservatorio Astronomico di Padova, Vicolo dell'Osservatorio, 5, 35122 Padova, Italy e-mail: olga.vega@oapd.inaf.it
³ SISSA, Strada Costiera, 34131 Trieste, Italy
⁴ INAF - Osservatorio Astronomico di Trieste, via Tiepolo 11, 34131 Trieste, Italy
⁵ Laboratoire AIM, CEA/DSM – CNRS – Université Paris Diderot, DAPNIA/Service d'Astrophysique, Bât. 709, CEA-Saclay, 91191 Gif-sur-Yvette Cedex, France

Received: 19 October 2007
Accepted: 23 February 2008

Abstract

Aims. By using the spectral energy distribution (SED) from the near-infrared to the radio of a statistically significant number of luminous infrared galaxies we determine important physical parameters for this population of objects. In particular we constrain the optical depth towards the luminosity source, the star formation rate, the star formation efficiency and the AGN fraction.

Methods. We fit the near-infrared to radio spectral energy distributions of a sample of 30 luminous and ultra-luminous infrared galaxies with pure starburst models or models that include both starburst and AGN components.

Results. We find that although about half of our sample have best-fit models that include an AGN component, only 30% (9/30) have an AGN that accounts for more than 10% of the infrared luminosity from 8 to 1000 μm, whereas all have an energetically dominant starburst. Our derived AGN fractions are generally in good agreement with measurements of the mid-infrared line ratios, Ne[V]/Ne[II] and O[IV]/Ne[II] by Spitzer IRS, but much lower than those derived from PAH equivalent widths or the mid-infrared spectral slope. Our models determine the mass of dense molecular gas within which active star formation takes place via the extinction required to reproduce the infrared part of the SED. Assuming that this mass is that traced by the HCN molecule, we reproduce the observed linear relation between HCN flux and infrared luminosity found previously. We also find that the star formation efficiency, as defined by the current star formation rate per unit molecular gas mass, falls as the starburst ages.

Conclusions. If the evolution of ULIRGs includes a phase in which an AGN contributes an important fraction to the infrared luminosity, this phase should last an order of magnitude less time than the starburst phase. However, we find no convincing evidence that an energetically important AGN is associated with a particular phase of the starburst. Because the mass of dense molecular gas that we derive is consistent with observations of the HCN molecule, it should be possible to estimate the mass of dense, star-forming molecular gas in such objects when molecular line data are not available.