Efstathiou, Andreas

We present Herschel observations of 6 fine-structure lines in 25 ultraluminous infrared galaxies at z < 0.27. The lines, [O III]52 μm, [N III]57 μm, [O I]63 μm, [N II]122 μm, [O I]145 μm, and [C II]158 μm, are mostly single Gaussians with widths <600 km s-1 and luminosities of 107-109 LO. There are deficits in the [O I]63/L IR, [N II]/L IR, [O I]145/L IR, and [C II]/L IR ratios compared to lower luminosity systems. The majority of the line deficits are consistent with dustier H II regions, but part of the [C II] deficit may arise from an additional mechanism, plausibly charged dust grains. This is consistent with some of the [C II] originating from photodissociation regions or the interstellar medium (ISM). We derive relations between far-IR line luminosities and both the IR luminosity and star formation rate. We find that [N II] and both [O I] lines are good tracers of the IR luminosity and star formation rate. In contrast, [C II] is a poor tracer of the IR luminosity and star formation rate, and does not improve as a tracer of either quantity if the [C II] deficit is accounted for. The continuum luminosity densities also correlate with the IR luminosity and star formation rate. We derive ranges for the gas density and ultraviolet radiation intensity of 101 < n < 102.5 and 102.2 < G 0 < 103.6, respectively. These ranges depend on optical type, the importance of star formation, and merger stage. We do not find relationships between far-IR line properties and several other parameters: active galactic nucleus (AGN) activity, merger stage, mid-IR excitation, and SMBH mass. We conclude that these far-IR lines arise from gas heated by starlight, and that they are not strongly influenced by AGN activity.

We present Herschel photometry and spectroscopy, carried out as part of the Herschel ultraluminous infrared galaxy (ULIRG) survey, and a model for the infrared to submillimetre emission of the ULIRG IRAS 08572+3915. This source shows one of the deepest known silicate absorption features and no polycyclic aromatic hydrocarbon emission. The model suggests that this object is powered by an active galactic nucleus (AGN) with a fairly smooth torus viewed almost edge-on and a very young starburst. According to our model, the AGN contributes about 90 per cent of the total luminosity of 1.1 × 1013 L⊙, which is about a factor of 5 higher than previous estimates. The large correction of the luminosity is due to theanisotropy of the emission of the best-fitting torus. Similar corrections may be necessary for other local and high-z analogues. This correction implies that IRAS 08572+3915 at a redshift of 0.058 35 may be the nearest hyperluminous infrared galaxy and probably the most luminous infrared galaxy in the local (z < 0.2) Universe. IRAS 08572+3915 shows a low ratio of [C II] to IR luminosity (log L[C II]/LIR < -3.8) and a [OI]63 μm to [CII]158 μm line ratio of about 1 that supports the model presented in this Letter.

We examine the power source of 41 local ultraluminous infrared galaxies (ULIRGs) using archival infrared (IR) and optical photometry. We fit the observed spectral energy distributions with starburst and active galactic nucleus (AGN) components, each component being drawn from a family of templates. We find that all of the sample require a starburst, whereas only half require an AGN. In 90 per cent of the sample the starburst provides over half the IR emission, with a mean fractional luminosity of 82 per cent. When combined with other galaxy samples we find that starburst and AGN luminosities correlate over six decades in IR luminosity, suggesting that a common factor governs both luminosities, plausibly the gas masses in the nuclear regions. We find no trend for increasing fractional AGN luminosity with increasing total luminosity, contrary to previous claims. We find that the mid-IR F7.7/C 7.7 line-continuum ratio is no indication of the starburst luminosity, or the fractional AGN luminosity, and therefore that F 7.7/C7.7 is not a reliable diagnostic of the power source in ULIRGs. The radio flux correlates with the starburst luminosity, but shows no correlation with the AGN luminosity, in line with previous results. We propose that the scatter in this correlation is due to a skewed starburst initial mass function and/or relic relativistic electrons from a previous starburst, rather than contamination from an obscured AGN. We show that most ULIRGs undergo multiple starbursts during their lifetime, and by inference that mergers between more than two galaxies may be common amongst ULIRGs. Our results support the evolutionary model for ULIRGs proposed by Farrah et al., where they can follow many different evolutionary paths of starburst and AGN activity in transforming merging spiral galaxies into elliptical galaxies, but that most do not go through an optical quasi-stellar object phase. The lower level of AGN activity in our local sample compared with z ∼ 1 hyperluminous infrared galaxies implies that the two samples are distinct populations. We postulate that different galaxy formation processes at high z are responsible for this difference.

We present the Herschel SPIRE Fourier Transform Spectroscopy (FTS) atlas for a complete flux-limited sample of local ultraluminous infrared galaxies (ULIRGs) as part of the HERschel Ultra Luminous InfraRed Galaxy Survey (HERUS). The data reduction is described in detail and was optimized for faint FTS sources ,with particular care being taken for the subtraction of the background, which dominates the continuum shape of the spectra. To improve the final spectra, special treatment in the data reduction has been given to any observation suffering from artifacts in the data caused by anomalous instrumental effects. Complete spectra are shown covering 200-671 μm, with photometry in the SPIRE bands at 250, 350, and 500 μm. The spectra include near complete CO ladders for over half of our sample, as well as fine structure lines from [C i] 370 μm, [C i] 609 μm, and [N ii] 205 μm. We also detect H2O lines in several objects. We construct CO spectral line energy distributions (SLEDs) for the sample, and compare their slopes with the far-infrared (FIR) colors and luminosities. We show that the CO SLEDs of ULIRGs can be broadly grouped into three classes based on their excitation. We find that the mid-J (5 < J < 8) lines are better correlated with the total FIR luminosity, suggesting that the warm gas component is closely linked to recent star formation. The higher J transitions do not linearly correlate with the FIR luminosity, consistent with them originating in hotter, denser gas that is unconnected to the current star formation. We conclude that in most cases more than one temperature component is required to model the CO SLEDs.

There is significant scientific value to be gained from combining AKARI fluxes with data at other far-infrared (IR) wavelengths from the Infrared Astronomical Satellite (IRAS) and Herschel missions. To be able to do this we must ensure that there are no systematic differences between the data sets that need to be corrected before the fluxes are compatible with each other. One such systematic effect identified in the Bright Source Catalog version 1 (BSCv1) data is the issue of beam corrections. We determine these for the BSC version 2 (BSCv2) data by correlating ratios of appropriate IRAS and AKARI bands with the difference in 2 Micron All Sky Survey (2MASS) J-band extended and point source magnitudes for sources cross-matched between the IRAS Faint Source Catalog (FSC), AKARI BSCv2 and 2MASS catalogs. We find significant correlations (p < 10 -13) indicating that beam corrections are necessary in the 65 and 90 μm bands. We then use these corrected fluxes to supplement existing data in spectral energy distribution (SED) fits for ultraluminous infrared galaxies (ULIRGs) in the Herschel ULIRG Survey (HERUS). The addition of AKARI fluxes makes little difference to the results of simple (T, β) fits to the SEDs of these sources, though there is a general decrease in reduced Ï ‡ 2 values. The utility of the extra AKARI data, however, is in allowing physically more realistic SED models with more parameters to be fitted to the data. We also extend our analysis of beam correction issues in the AKARI data by examining the Herschel Reference Sample (HRS) galaxies, which have Herschel photometry from 100 to 500 μm and which are more spatially extended than the HERUS ULIRGs. 34 of the HRS sources have good Herschel SEDs and matching data from AKARI. This investigation finds that our simple 2MASS-based beam correction scheme is inadequate for these larger and more complex sources. There are also indications that additional beam corrections at 140 and 160 μm are needed for these sources, extended on scales >1′.

We present optical to far-infrared photometry of 31 reddened QSOs that show evidence for radiatively driven outflows originating from active galactic nuclei (AGNs) in their rest-frame UV spectra. We use these data to study the relationships between the AGN-driven outflows, and the AGN and starburst infrared luminosities. We find that FeLoBAL QSOs are invariably IR-luminous, with IR luminosities exceeding 1012 L ⊙ in all cases. The AGN supplies 76% of the total IR emission, on average, but with a range from 20% to 100%. We find no evidence that the absolute luminosity of obscured star formation is affected by the AGN-driven outflows. Conversely, we find an anticorrelation between the strength of AGN-driven outflows, as measured from the range of outflow velocities over which absorption exceeds a minimal threshold, and the contribution from star formation to the total IR luminosity, with a much higher chance of seeing a starburst contribution in excess of 25% in systems with weak outflows than in systems with strong outflows. Moreover, we find no convincing evidence that this effect is driven by the IR luminosity of the AGN. We conclude that radiatively driven outflows from AGNs can have a dramatic, negative impact on luminous star formation in their host galaxies. We find that such outflows act to curtail star formation such that star formation contributes less than ∼ 25% of the total IR luminosity. We also propose that the degree to which termination of star formation takes place is not deducible from the IR luminosity of the AGN.

We report on our observations of the 79 and 119 μm doublet transitions of OH for 24 local (z < 0.262) ULIRGs observed with Herschel-PACS as part of the Herschel ULIRG Survey (HERUS). Some OH 119 μm profiles display a clear P-Cygni shape and therefore imply outflowing OH gas, while other profiles are predominantly in absorption or are completely in emission. We find that the relative strength of the OH emission component decreases as the silicate absorption increases. This result locates the OH outflows inside the obscured nuclei. The maximum outflow velocities for our sources range from less than 100 to 2000 km s-1, with 15/24 (10/24) sources showing OH absorption at velocities exceeding 700 km s-1 (1000 km s-1). Three sources show maximum OH outflow velocities exceeding that of Mrk231. Since outflow velocities above 500-700 km s-1 are thought to require an active galactic nucleus (AGN) to drive them, about two-thirds of our ULIRG sample may host AGN-driven molecular outflows. This finding is supported by the correlation we find between the maximum OH outflow velocity and the IR-derived bolometric AGN luminosity. No such correlation is found with the IR-derived star formation rate. The highest outflow velocities are found among sources that are still deeply embedded. We speculate that the molecular outflows in these sources may be in an early phase of disrupting the nuclear dust veil before these sources evolve into less-obscured AGNs. Four of our sources show high-velocity wings in their [C II] fine-structure line profiles, implying neutral gas outflow masses of at least (2-4.5) × 108 M.