Efstathiou, Andreas

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.

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 the ultraviolet to submillimetre spectral energy distributions (SEDs) of the HERschel Ultra Luminous Infrared Galaxy Survey (HERUS) sample of 42 local ultraluminous infrared galaxies (ULIRGs) and fit them with a Markov chain Monte Carlo code using the CYprus models for Galaxies and their NUclear Spectra (CYGNUS) radiative transfer models for starbursts, active galactic nucleus (AGN) tori, and host galaxy. The Spitzer IRS spectroscopy data are included in the fitting. Our Bayesian SED fitting method takes comparable time to popular energy balance methods but it is more physically motivated and versatile. All HERUS galaxies harbour high rates of star formation but we also find bolometrically significant AGN in all of the galaxies of the sample. We estimate the correction of the luminosities of the AGN in the ULIRGs due to the anisotropic emission of the torus and find that it could be up to a factor of ∼10 for nearly edge-on tori. We present a comparison of our results with the smooth torus model of Fritz et al. and the two-phase models of Siebenmorgen et al. and SKIRTOR. We find that the CYGNUS AGN torus models fit significantly better the SEDs of our sample compared to all other models. We find no evidence that strong AGN appear either at the beginning or end of a starburst episode or that starbursts and AGN affect each other. IRAS 01003-2238 and Mrk 1014 show evidence of dual AGNs in their SED fits, suggesting a minimum dual AGN fraction in the sample of 5 per cent.

In this work we report the discovery of the hyperluminous galaxy HELP J100156.75 + 022344.7 at a photometric redshift of z 4.3. The galaxy was discovered in the Cosmological Evolution Survey (COSMOS) field, one of the fields studied by the Herschel Extragalactic Legacy Project (HELP). We present the spectral energy distribution (SED) of the galaxy and fit it with the CYprus models for Galaxies and their NUclear Spectra (CYGNUS) multi-component radiative transfer models. We find that its emission is dominated by an obscured quasar with a predicted total 1-1000um luminosity of 3.91+1.69-0.55 × 1013 L and an active galactic nucleus (AGN) fraction of 89 per cent. We also fit HELP J100156.75 + 022344.7 with the CIGALE(CODE INVESTIGATING GALAXY EMISSION) code and find a similar result. This is only the second z 4 hyperluminous obscured quasar discovered to date. The discovery of HELP J100156.75 + 022344.7 in the 2 deg2 COSMOS field implies that a large number of obscured hyperluminous quasars may lie in the HELP fields, which cover 1300 deg2. If this is confirmed, tension between supermassive black hole evolution models and observations will be alleviated. We estimate the space density of objects like HELP J100156.75 + 022344.7 at z 4.5 to be 1.8 × 10-8 Mpc-3. This is slightly higher than the space density of coeval hyperluminous optically selected quasars, suggesting that the obscuring torus in z 4 quasars may have a covering factor 50 per cent.

We present a study of the infrared properties for a sample of seven spectroscopically confirmed submillimeter galaxies (SMGs) at z > 4.0. By combining ground-based near-infrared, Spitzer IRAC and MIPS, Herschel SPIRE, and ground-based submillimeter/millimeter photometry, we construct their spectral energy distributions (SEDs) and a composite model to fit the SEDs. The model includes a stellar emission component at λrest < 3.5 μm, a hot dust component peaking at λrest ∼ 5 μm, and cold dust component which becomes significant for λrest > 50 μm. Six objects in the sample are detected at 250 and 350 μm. The dust temperatures for the sources in this sample are in the range of 40-80 K, and their LFIR ∼ 1013 L⊙ qualifies them as hyper-luminous infrared galaxies. The mean FIR-radio index for this sample is around 〈q〉 = 2.2 indicating no radio excess in their radio emission. Most sources in the sample have 24 μm detections corresponding to a rest-frame 4.5 μm luminosity of Log10(L4.5/L ⊙) = 11 ∼ 11.5. Their L4.5/LFIR ratios are very similar to those of starburst-dominated SMGs at z ∼ 2. The L CO-LFIR relation for this sample is consistent with that determined for local ULIRGs and SMGs at z ∼ 2. We conclude that SMGs at z > 4 are hotter and more luminous in the FIR but otherwise very similar to those at z ∼ 2. None of these sources show any sign of the strong QSO phase being triggered.

We present the Herschel Extragalactic Legacy Project (HELP). This project collates, curates, homogenizes, and creates derived data products for most of the premium multiwavelength extragalactic data sets. The sky boundaries for the first data release cover 1270 deg2 defined by the Herschel SPIRE extragalactic survey fields; notably the Herschel Multi-tiered Extragalactic Survey (HerMES) and the Herschel Atlas survey (H-ATLAS). Here, we describe the motivation and principal elements in the design of the project. Guiding principles are transparent or 'open' methodologies with care for reproducibility and identification of provenance. A key element of the design focuses around the homogenization of calibration, meta data, and the provision of information required to define the selection of the data for statistical analysis. We apply probabilistic methods that extract information directly from the images at long wavelengths, exploiting the prior information available at shorter wavelengths and providing full posterior distributions rather than maximum-likelihood estimates and associated uncertainties as in traditional catalogues. With this project definition paper, we provide full access to the first data release of HELP; Data Release 1 (DR1), including a monolithic map of the largest SPIRE extragalactic field at 385 deg2 and 18 million measurements of PACS and SPIRE fluxes. We also provide tools to access and analyse the full HELP data base. This new data set includes far-infrared photometry, photometric redshifts, and derived physical properties estimated from modelling the spectral energy distributions over the full HELP sky. All the software and data presented is publicly available.

We present a new diagnostic diagram for local ultraluminous infrared galaxies (ULIRGs) and quasars, analysing particularly the Spitzer Space Telescope's infrared spectrograph spectra of 102 local ULIRGs and 37 Palomar Green quasars. Our diagram is based on a special non-linear mapping of these data, employing the kernel principal component analysis method. The novelty of this map lies in the fact that it distributes the galaxies under study on the surface of a well-defined ellipsoid, which, in turn, links basic concepts from geometry to physical properties of the galaxies. Particularly, we have found that the equatorial direction of the ellipsoid corresponds to the evolution of the power source of ULIRGs, starting from the pre-merger phase, moving through the starburst-dominated coalescing stage towards the active galactic nucleus-dominated phase, and finally terminating with the post-merger quasar phase. On the other hand, the meridian directions distinguish deeply obscured power sources of the galaxies from unobscured ones. These observations have also been verified by comparison with simulated ULIRGs and quasars using radiative transfer models. The diagram correctly identifies unique galaxies with extreme features that lie distinctly away from the main distribution of the galaxies. Furthermore, special two-dimensional projections of the ellipsoid recover almost monotonic variations of the two main physical properties of the galaxies, the silicate and polycyclic aromatic hydrocarbon features. This suggests that our diagram naturally extends the well-known Spoon diagram and it can serve as a diagnostic tool for existing and future infrared spectroscopic data, such as those provided by the James Webb Space Telescope.