Starburst and AGN activity in ultraluminous infrared galaxies

(Abridged) We examine the power source of 41 local Ultraluminous Infrared Galaxies using archival infrared and optical photometry. We fit the observed Spectral Energy Distributions (SEDs) with starburst and AGN components; each component being drawn from a family of templates. We find all of the sample require a starburst, whereas only half require an AGN. In 90% of the sample the starburst provides over half the IR emission, with a mean fractional luminosity of 82%. When combined with other galaxy samples we find that starburst and AGN luminosities correlate over 6 decades in IR luminosity suggesting that a common factor governs both luminosities, plausibly the gas masses in the nuclear regions. We find that the mid-IR 7.7 micron line-continuum ratio is no indication of the starburst luminosity, or the fractional AGN luminosity, and therefore that this ratio is not a reliable diagnostic of the power source in ULIRGs. We propose that the scatter in the radio-IR correlation in ULIRGs is due to a skewed starburst IMF 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 2001, 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 QSO phase. The lower level of AGN activity in our local sample than in z~1 HLIRGs implies that the two samples are distinct populations. We postulate that different galaxy formation processes at high-z are responsible for this difference.Comment: 24 pages, 8 figures. Accepted for publication in MNRA

The X-ray Spectrum of the z=6.30 QSO SDSS J1030+0524

We present a deep XMM-Newton observation of the z=6.30 QSO SDSS J1030+0524, the second most distant quasar currently known. The data contain sufficient counts for spectral analysis, demonstrating the ability of XMM-Newton to measure X-ray spectral shapes of z~6 QSOs with integration times >100ks. The X-ray spectrum is well fit by a power law with index Gamma=2.12 +/- 0.11, an optical-X-ray spectral slope of a_{ox}=-1.80, and no absorption excess to the Galactic value, though our data are also consistent with a power law index in the range 2.02 < Gamma < 2.5 and excess absorption in the range 0 < N_H(cm^-2) < 8x10^22. There is also a possible detection (2 sigma) of FeKa emission. The X-ray properties of this QSO are, overall, similar to those of lower-redshift radio-quiet QSOs. This is consistent with the statement that the X-ray properties of radio-quiet QSOs show no evolution over 0<z<6.3. Combined with previous results, this QSO appears indistinguishable in any way from lower redshift QSOs, indicating that QSOs comparable to those seen locally existed less than one Gyr after the Big Bang.Comment: ApJ Letters, accepte

Spectral energy distributions of quasars selected in the mid-infrared

We present preliminary results on fitting of SEDs to 142 z>1 quasars selected in the mid-infrared. Our quasar selection finds objects ranging in extinction from highly obscured, type-2 quasars, through more lightly reddened type-1 quasars and normal type-1s. We find a weak tendency for the objects with the highest far-infrared emission to be obscured quasars, but no bulk systematic offset between the far-infrared properties of dusty and normal quasars as might be expected in the most naive evolutionary schemes. The hosts of the type-2 quasars have stellar masses comparable to those of radio galaxies at similar redshifts. Many of the type-1s, and possibly one of the type-2s require a very hot dust component in addition to the normal torus emission.Comment: 4 pages, 2 figures, to appear in the proceedings of The Spectral Energy Distribution of Galaxies, Preston, September 2011, eds R.J. Tuffs & C.C. Popesc

A growth-rate indicator for Compton-thick active galactic nuclei

Due to their heavily obscured central engines, the growth rate of Compton-thick (CT) active galactic nuclei (AGN) is difficult to measure. A statistically significant correlation between the Eddington ratio, {\lambda}$_{Edd}$, and the X-ray power-law index, {\Gamma}, observed in unobscured AGN offers an estimate of their growth rate from X-ray spectroscopy (albeit with large scatter). However, since X-rays undergo reprocessing by Compton scattering and photoelectric absorption when the line-of-sight to the central engine is heavily obscured, the recovery of the intrinsic {\Gamma} is challenging. Here we study a sample of local, predominantly Compton-thick megamaser AGN, where the black hole mass, and thus Eddington luminosity, are well known. We compile results on X-ray spectral fitting of these sources with sensitive high-energy (E> 10 keV) NuSTAR data, where X-ray torus models which take into account the reprocessing effects have been used to recover the intrinsic {\Gamma} values and X-ray luminosities, L$_X$. With a simple bolometric correction to L$_X$ to calculate {\lambda}$_{Edd}$, we find a statistically significant correlation between {\Gamma} and {\lambda}$_{Edd}$ (p = 0.007). A linear fit to the data yields {\Gamma} = (0.41$\pm$0.18)log$_{10}${\lambda}$_{Edd}$+(2.38$\pm$ 0.20), which is statistically consistent with results for unobscured AGN. This result implies that torus modeling successfully recovers the intrinsic AGN parameters. Since the megamasers have low-mass black holes (M$_{BH}\approx10^6-10^7$ M$_{sol}$) and are highly inclined, our results extend the {\Gamma}-{\lambda}$_{Edd}$ relationship to lower masses and argue against strong orientation effects in the corona, in support of AGN unification. Finally this result supports the use of {\Gamma} as a growth-rate indicator for accreting black holes, even for Compton-thick AGN.Comment: Accepted for publication in Ap

Learning the fundamental mid-infrared spectral components of galaxies with non-negative matrix factorization

The mid-infrared (MIR) spectra observed with the Spitzer Infrared Spectrograph (IRS) provide a valuable data set for untangling the physical processes and conditions within galaxies. This paper presents the first attempt to blindly learn fundamental spectral components of MIR galaxy spectra, using non-negative matrix factorization (NMF). NMF is a recently developed multivariate technique shown to be successful in blind source separation problems. Unlike the more popular multivariate analysis technique, principal component analysis, NMF imposes the condition that weights and spectral components are non-negative. This more closely resembles the physical process of emission in the MIR, resulting in physically intuitive components. By applying NMF to galaxy spectra in the Cornell Atlas of Spitzer/IRS sources, we find similar components amongst different NMF sets. These similar components include two for active galactic nucleus (AGN) emission and one for star formation. The first AGN component is dominated by fine structure emission lines and hot dust, the second by broad silicate emission at 10 and 18 μm. The star formation component contains all the polycyclic aromatic hydrocarbon features and molecular hydrogen lines. Other components include rising continuums at longer wavelengths, indicative of colder grey-body dust emission. We show an NMF set with seven components can reconstruct the general spectral shape of a wide variety of objects, though struggle to fit the varying strength of emission lines. We also show that the seven components can be used to separate out different types of objects. We model this separation with Gaussian mixtures modelling and use the result to provide a classification tool. We also show that the NMF components can be used to separate out the emission from AGN and star formation regions and define a new star formation/AGN diagnostic which is consistent with all MIR diagnostics already in use but has the advantage that it can be applied to MIR spectra with low signal-to-noise ratio or with limited spectral range. The seven NMF components and code for classification are available at https://github.com/pdh21/NMF_software/

Keck spectroscopy of z=1-3 ULIRGs from the Spitzer SWIRE survey

(Abridged) High-redshift ultra luminous infrared galaxies contribute the bulk of the cosmic IR background and are the best candidates for very massive galaxies in formation at z>1.5. We present Keck/LRIS optical spectroscopy of 35 z>1.4 luminous IR galaxies in the Spitzer Wide-area Infra-Red Extragalactic survey (SWIRE) northern fields (Lockman Hole, ELAIS-N1, ELAIS-N2). The primary targets belong to the ``IR-peak'' class of galaxies, having the 1.6 micron (restframe) stellar feature detected in the IRAC Spitzer channels.The spectral energy distributions of the main targets are thoroughly analyzed, by means of spectro-photometric synthesis and multi-component fits (stars + starburst dust + AGN torus). The IR-peak selection technique is confirmed to successfully select objects above z=1.4, though some of the observed sources lie at lower redshift than expected. Among the 16 galaxies with spectroscopic redshift, 62% host an AGN component, two thirds being type-1 and one third type-2 objects. The selection, limited to r'<24.5, is likely biased to optically-bright AGNs. The SEDs of non-AGN IR-peakers resemble those of starbursts (SFR=20-500 Msun/yr) hosted in massive (M>1e11 Msun) galaxies. The presence of an AGN component provides a plausible explanation for the spectroscopic/photometric redshift discrepancies, as the torus produces an apparent shift of the peak to longer wavelengths. These sources are analyzed in IRAC and optical-IR color spaces. In addition to the IR-peak galaxies, we present redshifts and spectral properties for 150 objects, out of a total of 301 sources on slits.Comment: Accepted for publications on Astronomy and Astrophysics (acceprance date March 8th, 2007). 33 pages. The quality of some figures have been degrade

The extraordinary mid-infrared spectral properties of FeLoBAL Quasars

We present mid-infrared spectra of six FeLoBAL QSOs at 1<z<1.8, taken with the Spitzer space telescope. The spectra span a range of shapes, from hot dust dominated AGN with silicate emission at 9.7 microns, to moderately obscured starbursts with strong Polycyclic Aromatic Hydrocarbon (PAH) emission. The spectrum of one object, SDSS 1214-0001, shows the most prominent PAHs yet seen in any QSO at any redshift, implying that the starburst dominates the mid-IR emission with an associated star formation rate of order 2700 solar masses per year. With the caveats that our sample is small and not robustly selected, we combine our mid-IR spectral diagnostics with previous observations to propose that FeLoBAL QSOs are at least largely comprised of systems in which (a) a merger driven starburst is ending, (b) a luminous AGN is in the last stages of burning through its surrounding dust, and (c) which we may be viewing over a restricted line of sight range.Comment: ApJ, accepte