Source author record

J. R. Mullaney

J. R. Mullaney appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

44works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

44 published item(s)

preprint2022arXiv

Quasar Feedback Survey: Multi-phase outflows, turbulence and evidence for feedback caused by low power radio jets inclined into the galaxy disk

We present a study of a luminous, z=0.15, type-2 quasar (log [L([OIII])/(erg/s)]=42.8) from the Quasar Feedback Survey. It is classified as 'radio-quiet' (log [L(1.4 GHz)/(W/Hz)]=23.8); however, radio imaging reveals ~1 kpc low-power jets (log [Pjet/(erg/s)]=44) inclined into the plane of the galaxy disk. We combine MUSE and ALMA observations to map stellar kinematics and ionised and molecular gas properties. The jets are seen to drive galaxy-wide bi-conical turbulent outflows, reaching W80 = 1000-1300 km/s, in the ionised phase (traced via optical emission-lines), which also have increased electron densities compared to the quiescent gas. The turbulent gas is driven perpendicular to the jet axis and is escaping along the galaxy minor axis, reaching 7.5 kpc on both sides. Traced via CO(3-2) emission, the turbulent material in molecular gas phase is one-third as spatially extended and has 3 times lower velocity-dispersion as compared to ionised gas. The jets are seen to be strongly interacting with the interstellar medium (ISM) through enhanced ionised emission and disturbed/depleted molecular gas at the jet termini. We see further evidence for jet-induced feedback through significantly higher stellar velocity-dispersion aligned, and co-spatial with, the jet axis (<5 deg). We discuss possible negative and positive feedback scenarios arising due to the interaction of the low-power jets with the ISM in the context of recent jet-ISM interaction simulations, which qualitatively agree with our observations. We discuss how jet-induced feedback could be an important feedback mechanism even in bolometrically luminous 'radio-quiet' quasars.

preprint2021arXiv

The post-Herschel view of intrinsic AGN emission: constructing templates for galaxy and AGN emission at IR wavelengths

Measuring the star-forming properties of AGN hosts is key to our understanding of galaxy formation and evolution. However, this topic remains debated, partly due to the difficulties in separating the infrared (i.e. 1--1000 $μ$m) emission into AGN and star-forming components. Taking advantage of archival far-infrared data from Herschel, we present a new set of AGN and galaxy infrared templates, and introduce the spectral energy distribution fitting code IRAGNSEP. Both can be used to measure infrared host galaxy properties, free of AGN contamination. To build these, we used a sample of 100 local ($z$ < 0.3), low-to-high luminosity AGNs (i.e. $L_{\rm bol}~\sim~10^{42--46}~\rm erg~s^{-1}$), selected from the 105-month Swift - BAT X-ray survey, which have archival Spitzer - IRS spectra and Herschel photometry. We first built a set of seven galaxy templates using a sample of 55 star-forming galaxies selected via infrared diagnostics. Using these templates, combined with a flexible model for the AGN contribution, we extracted the intrinsic infrared emission of our AGN sample. We further demonstrate that we can reduce the diversity in the intrinsic shapes of AGN spectral energy distributions down to a set of three AGN templates, of which two represent AGN continuum, and one represents silicate emission. Our results indicate that, on average, the contribution of AGNs to the far-infrared ($λ~\gtrsim$ 50 $μ$m) is not as high as suggested by some recent work. We further show that the need for two infrared AGN continuum templates could be related to nuclear obscuration, where one of our templates appears dominated by the emission of the extended polar dust.

preprint2020arXiv

A binning-free method reveals a continuous relationship between galaxies' AGN power and offset from main sequence

Studies investigating the relationship between AGN power and the star formation rates (SFRs) of their host galaxies often rely on averaging techniques -- such as stacking -- to incorporate information from non-detections. However, averages, and especially means, can be strongly affected by outliers and can therefore give a misleading indication of the "typical" case. Recently, a number of studies have taken a step further by binning their sample in terms of AGN power (approximated by the 2-10keV luminosity of the AGN), and investigating how the SFR distribution differs between these bins. These bin thresholds are often weakly motivated, and binning implicitly assumes that sources within the same bin have similar (or even identical) properties. In this paper, we investigate whether the distribution of host SFRs -- relative to the locus of the star-forming main sequence (i.e., $R_{\rm MS}$) -- changes continuously as a function of AGN power. We achieve this by using a hierarchical Bayesian model that completely removes the need to bin in AGN power. In doing so, we find strong evidence that the $R_{\rm MS}$ distribution changes with 2-10keV X-ray luminosity. The results suggest that higher X_ray luminosity AGNs have a tighter physical connection to the star-forming process than lower X-ray luminosity AGNs, at least within the $0.8 < z < 1.2$ redshift range considered here.

preprint2020arXiv

Processing GOTO data with the Rubin Observatory LSST Science Pipelines I : Production of coadded frames

The past few decades have seen the burgeoning of wide field, high cadence surveys, the most formidable of which will be the Legacy Survey of Space and Time (LSST) to be conducted by the Vera C. Rubin Observatory. So new is the field of systematic time-domain survey astronomy, however, that major scientific insights will continue to be obtained using smaller, more flexible systems than the LSST. One such example is the Gravitational-wave Optical Transient Observer (GOTO), whose primary science objective is the optical follow-up of Gravitational Wave events. The amount and rate of data production by GOTO and other wide-area, high-cadence surveys presents a significant challenge to data processing pipelines which need to operate in near real-time to fully exploit the time-domain. In this study, we adapt the Rubin Observatory LSST Science Pipelines to process GOTO data, thereby exploring the feasibility of using this "off-the-shelf" pipeline to process data from other wide-area, high-cadence surveys. In this paper, we describe how we use the LSST Science Pipelines to process raw GOTO frames to ultimately produce calibrated coadded images and photometric source catalogues. After comparing the measured astrometry and photometry to those of matched sources from PanSTARRS DR1, we find that measured source positions are typically accurate to sub-pixel levels, and that measured L-band photometries are accurate to $\sim50$ mmag at $m_L\sim16$ and $\sim200$ mmag at $m_L\sim18$. These values compare favourably to those obtained using GOTO's primary, in-house pipeline, GOTOPHOTO, in spite of both pipelines having undergone further development and improvement beyond the implementations used in this study. Finally, we release a generic "obs package" that others can build-upon should they wish to use the LSST Science Pipelines to process data from other facilities.

preprint2020arXiv

The evolving AGN duty cycle in galaxies since z$\sim$3 as encoded in the X-ray luminosity function

We present a new modeling of the X-ray luminosity function (XLF) of Active Galactic Nuclei (AGN) out to z$\sim$3, dissecting the contribution of main-sequence (MS) and starburst (SB) galaxies. For each galaxy population, we convolved the observed galaxy stellar mass (M$_{\star}$) function with a grid of M$_{\star}$-independent Eddington ratio ($λ_{\rm EDD}$) distributions, normalised via empirical black hole accretion rate (BHAR) to star formation rate (SFR) relations. Our simple approach yields an excellent agreement with the observed XLF since z$\sim$3. We find that the redshift evolution of the observed XLF can only be reproduced through an intrinsic flattening of the $λ_{\rm EDD}$ distribution, and with a positive shift of the break $λ^{*}$, consistent with an anti-hierarchical behavior. The AGN accretion history is predominantly made by massive (10$^{10}<$M$_{\star}<$10$^{11}$ M$_{\odot}$) MS galaxies, while SB-driven BH accretion, possibly associated with galaxy mergers, becomes dominant only in bright quasars, at $\log$(L$_{\rm X}$/erg s$^{-1}$)$>$44.36 + 1.28$\cdot$(1+z). We infer that the probability of finding highly-accreting ($λ_{\rm EDD}>$ 10%) AGN significantly increases with redshift, from 0.4% (3.0%) at z=0.5 to 6.5% (15.3%) at z=3 for MS (SB) galaxies, implying a longer AGN duty cycle in the early Universe. Our results strongly favor a M$_{\star}$-dependent ratio between BHAR and SFR, as BHAR/SFR $\propto$ M$_{\star}^{0.73[+0.22,-0.29]}$, supporting a non-linear BH buildup relative to the host. Finally, this framework opens potential questions on super-Eddington BH accretion and different $λ_{\rm EDD}$ prescriptions for understanding the cosmic BH mass assembly.

preprint2019arXiv

The galaxy's gas content regulated by the dark matter halo mass results in a super-linear M$_{\rm BH}$-M$_{\star}$ relation

Supermassive black holes (SMBHs) are tightly correlated with their hosts but the origin of such connection remains elusive. To explore the cosmic build-up of this scaling relation, we present an empirically-motivated model that tracks galaxy and SMBH growth down to z=0. Starting from a random mass seed distribution at z=10, we assume that each galaxy evolves on the star-forming "main sequence" (MS) and each BH follows the recently-derived stellar mass (M$_{\star}$) dependent ratio between BH accretion rate and star formation rate, going as BHAR/SFR$\propto$M$_{\star}^{0.73[+0.22,-0.29]}$. Our simple recipe naturally describes the BH-galaxy build-up in two stages. At first, the SMBH lags behind the host that evolves along the MS. Later, as the galaxy grows in M$_{\star}$, our M$_{\star}$-dependent BHAR/SFR induces a super-linear BH growth, as M$_{\rm BH}$$\propto$M$_{\star}^{1.7}$. According to this formalism, smaller BH seeds increase their relative mass faster and earlier than bigger BH seeds, at fixed M$_{\star}$, thus setting along a gradually tighter M$_{\rm BH}$-M$_{\star}$ locus towards higher M$_{\star}$. Assuming reasonable values of the radiative efficiency $ε\sim$0.1, our empirical trend agrees with both high-redshift model predictions and intrinsic M$_{\rm BH}$-M$_{\star}$ relations of local BHs. We speculate that the observed non-linear BH-galaxy build-up is reflected in a twofold behavior with dark matter halo mass (M$_{\rm DM}$), displaying a clear turnover at M$_{\rm DM}\sim$2$\times$10$^{12}$M$_{\odot}$. While Supernovae-driven feedback suppresses BH growth in smaller halos (BHAR/SFR$\propto$M$_{\rm DM}^{1.6}$), above the M$_{\rm DM}$ threshold cold gas inflows possibly fuel both BH accretion and star formation in a similar fashion (BHAR/SFR$\propto$M$_{\rm DM}^{0.3}$).

preprint2016arXiv

ALMA observations of a z~3.1 Protocluster: Star Formation from Active Galactic Nuclei and Lyman-Alpha Blobs in an Overdense Environment

We exploit ALMA 870um observations to measure the star-formation rates (SFRs) of eight X-ray detected Active Galactic Nuclei (AGNs) in a z~3.1 protocluster, four of which reside in extended Ly-alpha haloes (often termed Ly-alpha blobs: LABs). Three of the AGNs are detected by ALMA and have implied SFRs of ~220-410~M_sun/yr; the non detection of the other five AGNs places SFR upper limits of <210 M_sun/yr. The mean SFR of the protocluster AGNs (~110-210 M_sun/yr) is consistent (within a factor of ~0.7-2.3) with that found for co-eval AGNs in the field, implying that galaxy growth is not significantly accelerated in these systems. However, when also considering ALMA data from the literature, we find evidence for elevated mean SFRs (up-to a factor of ~5.9 over the field) for AGNs at the protocluster core, indicating that galaxy growth is significantly accelerated in the central regions of the protocluster. We also show that all of the four protocluster LABs are associated with an ALMA counterpart within the extent of their Ly-alpha emission. The SFRs of the ALMA sources within the LABs (~150-410 M_sun/yr) are consistent with those expected for co-eval massive star-forming galaxies in the field. Furthermore, the two giant LABs (with physical extents of >100 kpc) do not host more luminous star formation than the smaller LABs, despite being an order of magnitude brighter in Ly-alpha emission. We use these results to discuss star formation as the power source of LABs.

preprint2016arXiv

ALMA resolves extended star formation in high-z AGN host galaxies

We present high resolution (0.3") Atacama Large Millimeter Array (ALMA) 870um imaging of five z~1.5-4.5 X-ray detected AGN with luminosities of L(2-8keV)>10^42 erg/s. These data provide a >~20x improvement in spatial resolution over single-dish rest-frame FIR measurements. The sub-millimetre emission is extended on scales of FWHM~0.2"-0.5", corresponding to physical sizes of 1-3 kpc (median value of 1.8 kpc). These sizes are comparable to the majority of z=1-5 sub-millimetre galaxies (SMGs) with equivalent ALMA measurements. In combination with spectral energy distribution analyses, we attribute this rest-frame far-infrared (FIR) emission to dust heated by star formation. The implied star-formation rate surface densities are ~20-200 Mo/yr/kpc^2, which are consistent with SMGs of comparable FIR luminosities (i.e., L(IR)~[1-5]x10^(12)Lo). Although limited by a small sample of AGN, which all have high FIR luminosities, our study suggests that the kpc-scale spatial distribution and surface density of star formation in high-redshift star-forming galaxies is the same irrespective of the presence of X-ray detected AGN.

preprint2016arXiv

An enhanced fraction of starbursting galaxies among high Eddington ratio AGNs

We investigate the star-forming properties of 1620 X-ray selected AGN host galaxies as a function of their specific X-ray luminosity (i.e., X-ray luminosity per unit host stellar mass) -- a proxy of the Eddington ratio. Our motivation is to determine whether there is any evidence of a suppression of star-formation at high Eddington ratios, which may hint toward "AGN feedback" effects. Star-formation rates (SFRs) are derived from fits to Herschel-measured far-infrared spectral energy distributions, taking into account any contamination from the AGN. Herschel-undetected AGNs are included via stacking analyses to provide average SFRs in bins of redshift and specific X-ray luminosity (spanning $0.01 \lesssim L_{\rm X}/M_{\ast} \lesssim 100~L_{\odot} ~M_{\odot}^{-1}$). After normalising for the effects of mass and redshift arising from the evolving galaxy main sequence, we find that the SFRs of high specific luminosity AGNs are slightly enhanced compared to their lower specific luminosity counterparts. This suggests that the SFR distribution of AGN hosts changes with specific X-ray luminosity, a result reinforced by our finding of a significantly higher fraction of starbursting hosts among high specific luminosity AGNs compared to that of the general star-forming galaxy population (i.e., 8-10 per cent vs. 3 per cent). Contrary to our original motivation, our findings suggest that high specific luminosity AGNs are more likely to reside in galaxies with enhanced levels of star-formation.

preprint2016arXiv

The NuSTAR Extragalactic Surveys: The Number Counts of Active Galactic Nuclei and the Resolved Fraction of the Cosmic X-ray Background

We present the 3-8 keV and 8-24 keV number counts of active galactic nuclei (AGN) identified in the NuSTAR extragalactic surveys. NuSTAR has now resolved 33-39% of the X-ray background in the 8-24 keV band, directly identifying AGN with obscuring columns up to approximately 1e25 /cm2. In the softer 3-8 keV band the number counts are in general agreement with those measured by XMM-Newton and Chandra over the flux range 5e-15 < S(3 - 8 keV)/(erg/cm2/s) < 1e-12 probed by NuSTAR. In the hard 8-24 keV band NuSTAR probes fluxes over the range 2e-14 < S(8-24 keV)/(erg/cm2/s) < 1e-12, a factor of approximately 100 fainter than previous measurements. The 8-24 keV number counts match predictions from AGN population synthesis models, directly confirming the existence of a population of obscured and/or hard X-ray sources inferred from the shape of the integrated cosmic X-ray background. The measured NuSTAR counts lie significantly above simple extrapolation with a Euclidian slope to low flux of the Swift/BAT 15-55 keV number counts measured at higher fluxes S(15-55 keV) > 1e-11 erg/cm2/s, reflecting the evolution of the AGN population between the Swift/BAT local (z<0.1) sample and NuSTAR's z~1 sample. CXB synthesis models, which account for AGN evolution, lie above the Swift/BAT measurements, suggesting that they do not fully capture the evolution of obscured AGN at low redshifts.

preprint2015arXiv

A remarkably flat relationship between the average star formation rate and AGN luminosity for distant X-ray AGN

In this study we investigate the relationship between the star formation rate, SFR, and AGN luminosity, L(AGN), for ~2000 X-ray detected AGN. The AGN span over three orders of magnitude in X-ray luminosity (10^(42) < L(2-8keV) < 10^(45.5) erg/s) and are in the redshift range z = 0.2 - 2.5. Using infrared (IR) photometry (8 - 500um), including deblended Spitzer and Herschel images and taking into account photometric upper limits, we decompose the IR spectral energy distributions into AGN and star formation components. Using the IR luminosities due to star formation, we investigate the average SFRs as a function of redshift and AGN luminosity. In agreement with previous studies, we find a strong evolution of the average SFR with redshift, tracking the observed evolution of the overall star forming galaxy population. However, we find that the relationship between the average SFR and AGN luminosity is flat at all redshifts and across all the AGN luminosities investigated; in comparison to previous studies, we find less scatter amongst the average SFRs across the wide range of AGN luminosities investigated. By comparing to empirical models, we argue that the observed flat relationship is due to short timescale variations in AGN luminosity, driven by changes in the mass accretion rate, which wash out any underlying correlations between SFR and L(AGN). Furthermore, we show that the exact form of the predicted relationship between SFR and AGN luminosity (and it's normalisation) is highly sensitive to the assumed intrinsic Eddington ratio distribution.

preprint2015arXiv

ALMA and Herschel reveal that X-ray selected AGN and main-sequence galaxies have different star formation rate distributions

Using deep Herschel and ALMA observations, we investigate the star formation rate (SFR) distributions of X-ray selected AGN host galaxies at 0.5<z<1.5 and 1.5<z<4, comparing them to that of normal, star-forming (i.e., "main-sequence", or MS) galaxies. We find 34--55 per cent of AGNs in our sample have SFRs at least a factor of two below that of the average MS galaxy, compared to ~15 per cent of all MS galaxies, suggesting significantly different SFR distributions. Indeed, when both are modelled as log-normal distributions, the mass and redshift-normalised SFR distributions of X-ray AGNs are roughly twice as broad, and peak ~0.4 dex lower, than that of MS galaxies. However, like MS galaxies, the normalised SFR distribution of AGNs in our sample appears not to evolve with redshift. Despite X-ray AGNs and MS galaxies having different SFR distributions, the linear-mean SFR of AGNs derived from our distributions is remarkably consistent with that of MS galaxies, and thus with previous results derived from stacked Herschel data. This apparent contradiction is due to the linear-mean SFR being biased by bright outliers, and thus does not necessarily represent a true characterisation of the typical SFR of X-ray AGNs.

preprint2015arXiv

Mid-infrared luminous quasars in the GOODS-Herschel fields: a large population of heavily-obscured, Compton-thick quasars at z~2

We present the infrared (IR) and X-ray properties of a sample of 33 mid-IR luminous quasars ($ν$L(6 micron)>6x10$^{44}$ erg/s) at redshift z~1-3, identified through detailed spectral energy distribution analyses of distant star-forming galaxies, using the deepest IR data from Spitzer and Herschel in the GOODS-Herschel fields. The aim is to constrain the fraction of obscured, and Compton-thick (CT, N$_H$>1.5x10$^{24}$ cm$^{-2}$) quasars at the peak era of nuclear and star-formation activities. Despite being very bright in the mid-IR band, ~30% of these quasars are not detected in the extremely deep 2 Ms and 4 Ms Chandra X-ray data available in these fields. X-ray spectral analysis of the detected sources reveals that the majority (~67%) are obscured by column densities N$_H$>10$^{22}$ cm$^{-2}$; this fraction reaches ~80% when including the X-ray undetected sources (9 out of 33), which are likely to be the most heavily-obscured, CT quasars. We constrain the fraction of CT quasars in our sample to be ~24-48%, and their space density to be $Φ$=(6.7$\pm$2.2)x10$^{-6}$ Mpc$^{-3}$. From the investigation of the quasar host galaxies in terms of star-formation rates (SFRs) and morphological distortions, as a sign of galaxy mergers/interactions, we do not find any direct relation between SFRs and quasar luminosity or X-ray obscuration. On the other hand, there is tentative evidence that the most heavily-obscured quasars have, on average, more disturbed morphologies than the unobscured/moderately-obscured quasar hosts, which preferentially live in undisturbed systems. However, the fraction of quasars with disturbed morphology amongst the whole sample is ~40%, suggesting that galaxy mergers are not the main fuelling mechanism of quasars at z~2.

preprint2015arXiv

Relationship between star formation rate and black hole accretion at z=2: The different contributions in quiescent, normal and starburst galaxies

We investigate the co-evolution of black-hole-accretion-rate (BHAR) and star-formation-rate (SFR) in $1.5<z<2.5$ galaxies displaying a greater diversity of star-forming properties compared to previous studies. We combine X-ray stacking and far-IR photometry of stellar mass-limited samples of normal star-forming, starburst and quiescent/quenched galaxies in the COSMOS field. We corroborate the existence of a strong correlation between BHAR (i.e. the X-ray luminosity, L_X), and stellar mass (M*) for normal star-forming galaxies, although find a steeper relation than previously reported. We find that starbursts show a factor of 3 enhancement in BHAR compared to normal SF galaxies (against a factor of 6 excess in SFR), while quiescents show a deficit of a factor 5.5 at a given mass. One possible interpretation of this is that the starburst phase does not coincide with cosmologically relevant BH growth, or that starburst-inducing mergers are more efficient at boosting SFR than BHAR. Contrary to studies based on smaller samples, we find the BHAR/SFR ratio of main sequence (MS) galaxies is not mass invariant, but scales weakly as M*^(0.43\pm0.09}, implying faster BH growth in more massive galaxies at $z\sim2$. Furthermore, BHAR/SFR during the starburst is a factor of 2 lower than in MS galaxies, at odds with the predictions of hydrodynamical simulations of merger galaxies that foresee a sudden enhancement of L_X/SFR during the merger. Finally, we estimate that the bulk of the accretion density of the Universe at $z\sim2$ is associated with normal star-forming systems, with only 6(+/-1)% and 11(+/-1)% associated with starburst and quiescent galaxies, respectively.

preprint2015arXiv

Storm in a "Teacup": a radio-quiet quasar with ~10kpc radio-emitting bubbles and extreme gas kinematics

We present multi-frequency (1-8 GHz) VLA data, combined with VIMOS IFU data and HST imaging, of a z=0.085 radio-quiet type 2 quasar (with L(1.4GHz)~5e23 W/Hz and L(AGN)~2e45 erg/s). Due to the morphology of its emission-line region, the target (J1430+1339) has been referred to as the Teacup AGN in the literature. We identify "bubbles" of radio emission that are extended ~10-12 kpc to both the east and west of the nucleus. The edge of the brighter eastern bubble is co-spatial with an arc of luminous ionized gas. We also show that the Teacup AGN hosts a compact radio structure, located ~0.8 kpc from the core position, at the base of the eastern bubble. This radio structure is co-spatial with an ionized outflow with an observed velocity of v=-740 km/s. This is likely to correspond to a jet, or possibly a quasar wind, interacting with the interstellar medium at this position. The large-scale radio bubbles appear to be inflated by the central AGN, which indicates that the AGN can also interact with the gas on >~10 kpc scales. Our study highlights that even when a quasar is formally "radio-quiet" the radio emission can be extremely effective for observing the effects of AGN feedback.

preprint2015arXiv

The composite nature of Dust-Obscured Galaxies (DOGs) at z~2-3 in the COSMOS field: I. A Far-Infrared View

Dust-Obscured galaxies (DOGs) are bright 24 um-selected sources with extreme obscuration at optical wavelengths. They are typically characterized by a rising power-law continuum of hot dust (T_D ~ 200-1000K) in the near-IR indicating that their mid-IR luminosity is dominated by an an active galactic nucleus (AGN). DOGs with a fainter 24 um flux display a stellar bump in the near-IR and their mid-IR luminosity appears to be mainly powered by dusty star formation. Alternatively, it may be that the mid-IR emission arising from AGN activity is dominant but the torus is sufficiently opaque to make the near-IR emission from the AGN negligible with respect to the emission from the host component. In an effort to characterize the astrophysical nature of the processes responsible for the IR emission in DOGs, this paper exploits Herschel data (PACS + SPIRE) on a sample of 95 DOGs within the COSMOS field. We derive a wealth of far-IR properties (e.g., total IR luminosities; mid-to-far IR colors; dust temperatures and masses) based on SED fitting. Of particular interest are the 24 um-bright DOGs (F_24um > 1mJy). They present bluer far-IR/mid-IR colors than the rest of the sample, unveiling the potential presence of an AGN. The AGN contribution to the total 8-1000 um flux increases as a function of the rest-frame 8 um-luminosity irrespective of the redshift. This confirms that faint DOGs (L_8 um< 10^12 L_sun) are dominated by star-formation while brighter DOGs show a larger contribution from an AGN.

preprint2015arXiv

The NuSTAR Extragalactic Surveys: First Direct Measurements of the >10 keV X-Ray Luminosity Function for Active Galactic Nuclei at z>0.1

We present the first direct measurements of the rest-frame 10-40 keV X-ray luminosity function (XLF) of Active Galactic Nuclei (AGNs) based on a sample of 94 sources at 0.1 < z <3, selected at 8-24 keV energies from sources in the NuSTAR extragalactic survey program. Our results are consistent with the strong evolution of the AGN population seen in prior, lower-energy studies of the XLF. However, different models of the intrinsic distribution of absorption, which are used to correct for selection biases, give significantly different predictions for the total number of sources in our sample, leading to small, systematic differences in our binned estimates of the XLF. Adopting a model with a lower intrinsic fraction of Compton-thick sources and a larger population of sources with column densities N_H ~ 10^{23-24} /cm2 or a model with stronger Compton reflection component (with a relative normalization of R ~ 2 at all luminosities) can bring extrapolations of the XLF from 2-10 keV into agreement with our NuSTAR sample. Ultimately, X-ray spectral analysis of the NuSTAR sources is required to break this degeneracy between the distribution of absorbing column densities and the strength of the Compton reflection component and thus refine our measurements of the XLF. Furthermore, the models that successfully describe the high-redshift population seen by NuSTAR tend to over-predict previous, high-energy measurements of the local XLF, indicating that there is evolution of the AGN population that is not fully captured by the current models.

preprint2015arXiv

The NuSTAR Extragalactic Surveys: Initial Results and Catalog from the Extended Chandra Deep Field South

We present initial results and the source catalog from the NuSTAR survey of the Extended Chandra Deep Field South (hereafter, ECDFS) - currently the deepest contiguous component of the NuSTAR extragalactic survey program. The survey covers the full ~30 arcmin x 30 arcmin area of this field to a maximum depth of ~360 ks (~220 ks when corrected for vignetting at 3-24 keV), reaching sensitivity limits of ~1.3 x 10^-14 erg/cm2/s (3-8 keV), ~3.4 x 10^-14 erg/cm2/s (8-24 keV) and ~3.0 x 10^-14 erg/cm2/s (3-24 keV). Fifty four (54) sources are detected over the full field, although five of these are found to lie below our significance threshold once contaminating flux from neighboring (i.e., blended) sources is taken into account. Of the remaining 49 that are significant, 19 are detected in the 8-24 keV band. The 8-24 keV to 3-8 keV band ratios of the twelve sources that are detected in both bands span the range 0.39-1.7, corresponding to a photon index range of Gamma ~ 0.5-2.3, with a median photon index of 1.70 +/- 0.52. The redshifts of the 49 sources in our main sample span the range z = 0.21-2.7, and their rest-frame 10-40 keV luminosities (derived from the observed 8-24 keV fluxes) span the range L(10-40 keV) ~ (0.7-300) x 10^43 erg/s, sampling below the "knee" of the X-ray luminosity function out to z ~ 0.8-1. Finally, we identify one NuSTAR source that has neither a Chandra nor an XMM-Newton counterpart, but that shows evidence of nuclear activity at infrared wavelengths, and thus may represent a genuine, new X-ray source detected by NuSTAR in the ECDFS.

preprint2015arXiv

The NuSTAR Extragalactic Surveys: Overview and Catalog from the COSMOS Field

To provide the census of the sources contributing to the X-ray background peak above 10 keV, NuSTAR is performing extragalactic surveys using a three-tier "wedding cake" approach. We present the NuSTAR survey of the COSMOS field, the medium sensitivity and medium area tier, covering 1.7 deg2 and overlapping with both Chandra and XMM-Newton data. This survey consists of 121 observations for a total exposure of ~3 Ms. To fully exploit these data, we developed a new detection strategy, carefully tested through extensive simulations. The survey sensitivity at 20% completeness is 5.9, 2.9 and 6.4 x 10^-14 erg/cm2/s in the 3-24 keV, 3-8 keV and 8-24 keV bands, respectively. By combining detections in 3 bands, we have a sample of 91 NuSTAR sources with 10^42 -10^45.5 erg/s luminosities and redshift z=0.04-2.5. Thirty two sources are detected in the 8-24 keV band with fluxes ~100 times fainter than sources detected by Swift-BAT. Of the 91 detections, all but four are associated with a Chandra and/or XMM-Newton point-like counterpart. One source is associated with an extended lower energy X-ray source. We present the X-ray (hardness ratio and luminosity) and optical-to-X-ray properties. The observed fraction of candidate Compton-thick AGN measured from the hardness ratio is between 13-20%. We discuss the spectral properties of NuSTAR J100259+0220.6 (ID 330) at z=0.044, with the highest hardness ratio in the entire sample. The measured column density exceeds 10^24 /cm2, implying the source is Compton-thick. This source was not previously recognized as such without the >10 keV data.

preprint2014arXiv

Kiloparsec-scale outflows are prevalent among luminous AGN: outflows and feedback in the context of the overall AGN population

We present integral field unit (IFU) observations covering the [O III]4959,5007 and H-Beta emission lines of sixteen z<0.2 type 2 active galactic nuclei (AGN). Our targets are selected from a well-constrained parent sample of 24,000 AGN so that we can place our observations into the context of the overall AGN population. Our targets are radio-quiet with star formation rates (<~[10-100] Msol/yr) that are consistent with normal star-forming galaxies. We decouple the kinematics of galaxy dynamics and mergers from outflows. We find high-velocity ionised gas (velocity widths of 600-1500 km/s and maximum velocities of <=1700 km/s) with observed spatial extents of >~(6-16) kpc in all targets and observe signatures of spherical outflows and bi-polar superbubbles. We show that our targets are representative of z<0.2, luminous (i.e., L([O III]) > 5x10^41 erg/s) type 2 AGN and that ionised outflows are not only common but also in >=70% (3 sigma confidence) of cases, they are extended over kiloparsec scales. Our study demonstrates that galaxy-wide energetic outflows are not confined to the most extreme star-forming galaxies or radio-luminous AGN; however, there may be a higher incidence of the most extreme outflow velocities in quasars hosted in ultra-luminous infrared galaxies. Both star formation and AGN activity appear to be energetically viable to drive the outflows and we find no definitive evidence that favours one process over the other. Although highly uncertain, we derive mass outflow rates (typically ~10x the SFRs), kinetic energies (~0.5-10% of L[AGN]) and momentum rates (typically >~10-20x L[AGN]/c) consistent with theoretical models that predict AGN-driven outflows play a significant role in shaping the evolution of galaxies.

preprint2014arXiv

Modeling the connection between ultraviolet and infrared galaxy populations across cosmic times

Using a phenomenological approach, we self-consistently model the redshift evolution of the ultraviolet (UV) and infrared (IR) luminosity functions across cosmic time, as well as a range of observed IR properties of UV-selected galaxy population. This model is an extension of the 2SFM (2 star-formation modes) formalism, which is based on the observed "main-sequence" of star-forming galaxies, i.e. a strong correlation between their stellar mass and their star formation rate (SFR), and a secondary population of starbursts with an excess of star formation. The balance between the UV light from young, massive stars and the dust-reprocessed IR emission is modeled following the empirical relation between the attenuation (IRX for IR excess hereafter) and the stellar mass, assuming a scatter of 0.4\,dex around this relation. We obtain a good overall agreement with the measurements of the IR luminosity function up to z~3 and the UV luminosity functions up to z~6, and show that a scatter on the IRX-M relation is mandatory to reproduce these observables. We also naturally reproduce the observed, flat relation between the mean IRX and the UV luminosity at L$_{\rm UV}>$10$^{9.5}$ L$_\odot$. Finally, we perform predictions of the UV properties and detectability of IR-selected samples and the vice versa, and discuss the results in the context of the UV-rest-frame and sub-millimeter surveys of the next decade.

preprint2014arXiv

NuSTAR J033202-2746.8: direct constraints on the Compton reflection in a heavily obscured quasar at z~2

We report NuSTAR observations of NuSTAR J033202-2746.8, a heavily obscured, radio-loud quasar detected in the Extended Chandra Deep Field-South, the deepest layer of the NuSTAR extragalactic survey (~400 ks, at its deepest). NuSTAR J033202-2746.8 is reliably detected by NuSTAR only at E>8 keV and has a very flat spectral slope in the NuSTAR energy band (Gamma=0.55^{+0.62}_{-0.64}; 3-30 keV). Combining the NuSTAR data with extremely deep observations by Chandra and XMM-Newton (4 Ms and 3 Ms, respectively), we constrain the broad-band X-ray spectrum of NuSTAR J033202-2746.8, indicating that this source is a heavily obscured quasar (N_H=5.6^{+0.9}_{-0.8}x10^23 cm^-2) with luminosity L_{10-40 keV}~6.4x10^44 erg s^-1. Although existing optical and near-infrared (near-IR) data, as well as follow-up spectroscopy with the Keck and VLT telescopes, failed to provide a secure redshift identification for NuSTAR J033202-2746.8, we reliably constrain the redshift z=2.00+/-0.04 from the X-ray spectral features (primarily from the iron K edge). The NuSTAR spectrum shows a significant reflection component (R=0.55^{+0.44}_{-0.37}), which was not constrained by previous analyses of Chandra and XMM-Newton data alone. The measured reflection fraction is higher than the R~0 typically observed in bright radio-loud quasars such as NuSTAR J033202-2746.8, which has L_{1.4 GHz}~10^27 W Hz^-1. Constraining the spectral shape of AGN, including bright quasars, is very important for understanding the AGN population, and can have a strong impact on the modeling of the X-ray background. Our results show the importance of NuSTAR in investigating the broad-band spectral properties of quasars out to high redshift.

preprint2014arXiv

NuSTAR Observations of Heavily Obscured Quasars at z ~ 0.5

We present NuSTAR hard X-ray (3-79 keV) observations of three Type 2 quasars at z ~ 0.4-0.5, optically selected from the Sloan Digital Sky Survey (SDSS). Although the quasars show evidence for being heavily obscured Compton-thick systems on the basis of the 2-10 keV to [OIII] luminosity ratio and multiwavelength diagnostics, their X-ray absorbing column densities (N_H) are poorly known. In this analysis: (1) we study X-ray emission at >10 keV, where X-rays from the central black hole are relatively unabsorbed, in order to better constrain N_H; (2) we further characterize the physical properties of the sources through broad-band near-UV to mid-IR spectral energy distribution (SED) analyses. One of the quasars is detected with NuSTAR at >8 keV with a no-source probability of <0.1%, and its X-ray band ratio suggests near Compton-thick absorption with N_H \gtrsim 5 x 10^23 cm^-2. The other two quasars are undetected, and have low X-ray to mid-IR luminosity ratios in both the low energy (2-10 keV) and high energy (10-40 keV) X-ray regimes that are consistent with extreme, Compton-thick absorption (N_H \gtrsim 10^24 cm^-2). We find that for quasars at z ~ 0.5, NuSTAR provides a significant improvement compared to lower energy (<10 keV) Chandra and XMM-Newton observations alone, as higher column densities can now be directly constrained.

preprint2014arXiv

Rapidly growing black holes and host galaxies in the distant Universe from the Herschel Radio Galaxy Evolution Project

We present results from a survey of 70 radio galaxies (RGs) at redshifts 1<z<5.2 using the PACS and SPIRE on-board Herschel. Combined with existing mid-IR photometry from Spitzer and observations obtained with LABOCA, the SEDs of galaxies in our sample are continuously covered across 3.6-870um. The total infrared luminosities of these RGs are such that they almost all are either ultra-or hyper-luminous infrared galaxies. We fit the infrared SEDs with a set of empirical templates which represent dust heated (1) by a variety of SB and (2) by a AGN. We find that the SEDs of RGs require the dust to be heated by both AGN and SB, but the luminosities of these two components are not strongly correlated. Assuming empirical relations and simple physical assumptions, we calculate the SFR, the black hole mass accretion rate (MdotBH), and the black hole mass (MBH) for each RG. We find that the host galaxies and their BHs are growing extremely rapidly, having SFR~100-5000 Msun/yr and MdotBH~1-100 Msun/yr. The mean sSFR of RGs at z>2.5 are higher than the sSFR of typical star-forming galaxies over the same redshift range but are similar or perhaps lower than the galaxy population for RGs at z<2.5. By comparing the sSFR and the specific black hole mass accretion rate, we conclude that BHs in radio loud AGN are already, or soon will be, overly massive compared to their host galaxies in terms of expectations from the local MBH-MGal relation. In order to ``catch up'' with the BH, the galaxies require about an order-of magnitude more time to grow in mass, at the observed SFRs, compared to the time the BH is actively accreting. However, during the current cycle of activity, we argue that this catching-up is likely to be difficult due to the short gas depletion times. Finally, we speculate on how the host galaxies might grow sufficiently in stellar mass to ultimately fall onto the local MBH-MGal relation.

preprint2014arXiv

The host galaxies of X-ray selected Active Galactic Nuclei to z=2.5: Structure, star-formation and their relationships from CANDELS and Herschel/PACS

We study the relationship between the structure and star-formation rate (SFR) of X-ray selected low and moderate luminosity active galactic nuclei (AGNs) in the two Chandra Deep Fields, using Hubble Space Telescope imaging from the Cosmic Assembly Near Infrared Extragalactic Legacy Survey (CANDELS) and deep far-infrared maps from the PEP+GOODS-Herschel survey. We derive detailed distributions of structural parameters and FIR luminosities from carefully constructed control samples of galaxies, which we then compare to those of the AGNs. At z~1, AGNs show slightly diskier light profiles than massive inactive (non-AGN) galaxies, as well as modestly higher levels of gross galaxy disturbance (as measured by visual signatures of interactions and clumpy structure). In contrast, at z~2, AGNs show similar levels of galaxy disturbance as inactive galaxies, but display a red central light enhancement, which may arise due to a more pronounced bulge in AGN hosts or due to extinguished nuclear light. We undertake a number of tests of these alternatives, but our results do not strongly favour one interpretation over the other. The mean SFR and its distribution among AGNs and inactive galaxies are similar at z>1.5. At z<1, however, clear and significant enhancements are seen in the SFRs of AGNs with bulge-dominated light profiles. These trends suggest an evolution in the relation between nuclear activity and host properties with redshift, towards a minor role for mergers and interactions at z>1.5.

preprint2013arXiv

A wide search for obscured Active Galactic Nuclei using XMM-Newton and WISE

We use a combination of the XMM-Newton serendipitous X-ray survey with the optical SDSS, and the infrared WISE all-sky survey in order to check the efficiency of the low X-ray to infrared luminosity selection method in finding heavily obscured AGN. We select sources in the 2-8 keV X-ray band which have a redshift determination in the SDSS catalogue. We match this sample with the WISE catalogue, and fit the SEDs of the 2844 sources which have three, or more, photometric data-points in the infrared. We then select the heavily obscured AGN candidates by comparing their 12 micron AGN luminosity to the observed 2-10 keV X-ray luminosity and their expected intrinsic relation. With this approach we find 20 candidates, and we examine their X-ray and optical spectra. Of the 20 initial candidates, we find nine (64%; out of the 14, for which X-ray spectra could be fit) based on the X-ray spectra, and seven (78%; out of the nine detected spectroscopically in the SDSS) based on the [OIII] line fluxes. Combining all criteria, we determine the final number of heavily obscured AGN to be 12-19, and the number of Compton-thick AGN to be 2-5, showing that the method is reliable in finding obscured AGN, but not Compton-thick. However those numbers are smaller than what would be expected from X-ray background population synthesis models, which demonstrates how the optical-infrared selection and the scatter of the L_x-L_MIR relation introduced by observational constraints limit the efficiency of the method. Finally, we test popular obscured AGN selection methods based on mid-infrared colours, and find that the probability of an AGN to be selected by its mid-infrared colours increases with the X-ray luminosity. However, a selection scheme based on a relatively low X-ray luminosity and mid-infrared colours characteristic of QSOs would not select ~25% of the heavily obscured AGN of our sample. (abridged)

preprint2013arXiv

Concurrent Supermassive Black Hole and Galaxy Growth: Linking Environment and Nuclear Activity in z = 2.23 H-alpha Emitters

We present results from a ~100 ks Chandra observation of the 2QZ Cluster 1004+00 structure at z = 2.23 (hereafter, 2QZ Clus). 2QZ Clus was originally identified as an overdensity of four optically-selected QSOs at z = 2.23 within a 15x15 arcmin^2 region. Narrow-band imaging in the near-IR revealed that the structure contains an additional overdensity of 22 z = 2.23 Halpha-emitting galaxies (HAEs), resulting in 23 unique z = 2.23 HAEs/QSOs. Our Chandra observations reveal that 3 HAEs in addition to the 4 QSOs harbor powerfully accreting supermassive black holes (SMBHs), with 2-10 keV luminosities of ~(8-60) x 10^43 ergs/s and X-ray spectral slopes consistent with unobscured AGN. Using a large comparison sample of 210 HAEs in Chandra-COSMOS (C-COSMOS), we find suggestive evidence that the AGN fraction increases with local HAE galaxy density. The 2QZ Clus HAEs reside in a moderately overdense environment (a factor of ~2 times over the field), and after excluding optically-selected QSOs, we find the AGN fraction is a factor of ~3.5^+3.8_-2.2 times higher than C-COSMOS HAEs in similar environments. Using stacking analyses of the Chandra data and Herschel SPIRE observations at 250 um, we respectively estimate mean SMBH accretion rates (Mdot) and star-formation rates (SFRs) for the 2QZ Clus and C-COSMOS samples. We find that the mean 2QZ Clus HAE stacked 2-10 keV luminosity is QSO-like (~6-10 x 10^43 ergs/s), and the implied Mdot/SFR \sim (1.6-3.2) x10^-3 is broadly consistent with the local MBH/M* relation and z ~ 2 X-ray selected AGN. The C-COSMOS HAEs are on average an order of magnitude less X-ray luminous and have Mdot/SFR \sim (0.2-0.4) x10^-3, comparable to z ~ 1-2 star-forming galaxies with similar mean X-ray luminosities. We estimate that a periodic QSO phase with duty cycle ~2-8% would be sufficient to bring the star-forming galaxies onto the local MBH/M* relation.

preprint2013arXiv

Narrow-line region gas kinematics of 24,264 optically-selected AGN: the radio connection

Using a sample of 24264 optically selected AGNs from the SDSS DR7 database, we characterise how the profile of the [OIII] emission line relates to bolometric luminosity (L_Bol), Eddington ratio, radio loudness, radio luminosity (L_Rad) and optical class (Sy 1/2) to determine what drives the kinematics of this kpc-scale line emitting gas. Spectral stacking is used to characterise how the average [OIII] profile changes as a function of these variables. After accounting for the correlation between L_Bol and L_Rad, we report that L_Rad has the strongest influence on the [OIII] profile, with moderate radio luminosity AGNs (log(L_Rad)=23-25 W/Hz) having the broadest [OIII] profiles. When binned according to Eddington ratio, only AGNs in our highest bin (i.e., R_Edd>0.3) show any signs of having broadened [OIII] profiles, although the small numbers of such extreme AGNs mean we cannot rule out that other processes (e.g., radio jets) are responsible for this broadening. We find no significant difference between the [OIII] profiles of broad and narrow-line Seyfert 1s. The [OIII] profiles of Sy 1 and 2 AGNs show the same trends in terms of line width, but Sy 1 AGNs display a much stronger blue wing, which we interpret as evidence of outflowing ionised gas. Using multi-component line fitting we calculate the proportions of AGNs with broad [OIII] profiles, finding that moderate radio luminosity AGNs are roughly 5 times more likely to have extremely broad [OIII] lines (FWHM>1000km/s) compared to lower L_Rad AGNs. Our results are consistent with the most disturbed gas kinematics being induced by compact radio cores (rather than powerful radio jets), although broadened [OIII] lines are also present, but much rarer, in low L_Rad systems (<10^23 W/Hz). Our catalogue of multi-component fits is freely available as an online resource.

preprint2013arXiv

New insights on the distant AGN population

Current X-ray surveys have proved to be essential tools in order to identify and study AGNs across cosmic time. However, there is evidence that the most heavily obscured AGNs are largely missing even in the deepest surveys. The search for these obscured AGNs is one of the most outstanding issues of extragalactic astronomy, since they are expected to make a major contribution to the high energy peak of the X-ray background (XRB) and might constitute a particularly active and dusty phase of black hole and galaxy evolution. Using a newly developed SED fitting technique to decompose the AGN from the star-formation emission in the infrared band we identify a sample of 17 IR bright quasars at z=1-3. For the majority of these sources the X-ray spectrum is well characterised by an absorbed power law model, revealing that ~25% of the sources are unabsorbed; the remainder of these IR quasars have moderate-to-high absorption (N_H>10^23 cm^-2), and ~25% are not detected in the X-rays and are likely to be Compton thick (N_H>10^24 cm^-2). We therefore find a much higher fraction of obscured than unobscured quasars, indicating that a large fraction of the luminous black hole accretion was very heavily obscured at z~2.

preprint2013arXiv

Nuclear Activity is more prevalent in Star-Forming Galaxies

We explore the question of whether low and moderate luminosity Active Galactic Nuclei (AGNs) are preferentially found in galaxies that are undergoing a transition from active star formation to quiescence. This notion has been suggested by studies of the UV-to-optical colors of AGN hosts, which find them to be common among galaxies in the so-called "Green Valley", a region of galaxy color space believed to be composed mostly of galaxies undergoing star-formation quenching. Combining the deepest current X-ray and Herschel. PACS far-infrared (FIR) observations of the two Chandra Deep Fields (CDFs) with redshifts, stellar masses and rest-frame photometry derived from the extensive and uniform multi-wavelength data in these fields, we compare the rest-frame U-V color distributions and SFR distributions of AGNs and carefully constructed samples of inactive control galaxies. The UV-to-optical colors of AGNs are consistent with equally massive inactive galaxies at redshifts out to z~2, but we show that such colors are poor tracers of star formation. While the FIR distributions of both star-forming AGNs and star-forming inactive galaxies are statistically similar, we show that AGNs are preferentially found in star-forming host galaxies, or, in other words, AGNs are less likely to be found in weakly star-forming or quenched galaxies. We postulate that, among X-ray selected AGNs of low and moderate accretion luminosities, the supply of cold gas primarily determines the accretion rate distribution of the nuclear black holes.

preprint2013arXiv

The NuSTAR Extragalactic Survey: A First Sensitive Look at the High-Energy Cosmic X-ray Background Population

We report on the first ten identifications of sources serendipitously detected by the NuSTAR to provide the first sensitive census of the cosmic X-ray background (CXB) source population at >10 keV. We find that these NuSTAR-detected sources are ~100x fainter than those previously detected at >10 keV and have a broad range in redshift and luminosity (z=0.020-2.923 and L_10-40 keV~4x10^{41}-5x10^{45} erg/s); the median redshift and luminosity are z~0.7 and L_10-40 keV~3x10^{44} erg/s, respectively. We characterize these sources on the basis of broad-band ~0.5-32 keV spectroscopy, optical spectroscopy, and broad-band ultraviolet-to-mid-infrared SED analyzes. We find that the dominant source population is quasars with L_10-40 keV>10^{44} erg/s, of which ~50% are obscured with N_H>10^{22} cm^{-2}. However, none of the ten NuSTAR sources are Compton thick (N_H>10^{24} cm^{-2}) and we place a 90% confidence upper limit on the fraction of Compton-thick quasars (L_10-40 keV>10^{44} erg/s) selected at >10 keV of ~33% over the redshift range z=0.5-1.1. We jointly fitted the rest-frame ~10-40 keV data for all of the non-beamed sources with L_10-40 keV>10^{43} erg/s to constrain the average strength of reflection; we find R<1.4 for Γ=1.8, broadly consistent with that found for local AGNs observed at >10 keV. We also constrain the host galaxy masses and find a median stellar mass of ~10^{11} M_sun, a factor ~5 times higher than the median stellar mass of nearby high-energy selected AGNs, which may be at least partially driven by the order of magnitude higher X-ray luminosities of the NuSTAR sources. Within the low source-statistic limitations of our study, our results suggest that the overall properties of the NuSTAR sources are broadly similar to those of nearby high-energy selected AGNs but scaled up in luminosity and mass.

preprint2012arXiv

Deep Silicate absorption features in Compton-thick AGN predominantly arise due to dust in the host galaxy

We explore the origin of mid-infrared (mid-IR) dust extinction in all 20 nearby (z < 0.05) bona-fide Compton-thick (N_H > 1.5 x 10^24 cm^-2) AGN with hard energy (E > 10 keV) X-ray spectral measurements. We accurately measure the silicate absorption features at lambda~9.7um in archival low-resolution (R~57-127) Spitzer Infrared Spectrograph (IRS) spectroscopy, and show that only a minority (~45%) of nearby Compton-thick AGN have strong Si-absorption features (S_9.7 = ln(f_{int}/f_{obs}) > 0.5) which would indicate significant dust attenuation. The majority (~60%) are star-formation dominated (AGN:SB<0.5) at mid-IR wavelengths and lack the spectral signatures of AGN activity at optical wavelengths, most likely because the AGN emission-lines are optically-extinguished. Those Compton-thick AGN hosted in low-inclination angle galaxies exhibit a narrow-range in Si-absorption (S_9.7 ~ 0-0.3), which is consistent with that predicted by clumpy-torus models. However, on the basis of the IR spectra and additional lines of evidence, we conclude that the dominant contribution to the observed mid-IR dust extinction is dust located in the host galaxy (i.e., due to disturbed morphologies; dust-lanes; galaxy inclination angles) and not necessarily a compact obscuring torus surrounding the central engine.

preprint2012arXiv

Energetic galaxy-wide outflows in high-redshift ultra-luminous infrared galaxies hosting AGN activity

We present integral field spectroscopy observations, covering the [O III]4959,5007 emission-line doublet of eight high-redshift (z=1.4-3.4) ultra-luminous infrared galaxies (ULIRGs) that host Active Galactic Nuclei (AGN) activity, including known sub-millimetre luminous galaxies (SMGs). The targets have moderate radio luminosities that are typical of high-redshift ULIRGs (L(1.4GHz)=10^(24)-10^(25)W/Hz) and therefore are not radio-loud AGN. We de-couple kinematic components due to the galaxy dynamics and mergers from those due to outflows. We find evidence in the four most luminous systems (L([O III])>~10^(43)erg/s) for the signatures of large-scale energetic outflows: extremely broad [O III] emission (FWHM ~ 700-1400km/s) across ~4-15kpc, with high velocity offsets from the systemic redshifts (up to ~850km/s). The four less luminous systems have lower quality data displaying weaker evidence for spatially extended outflows. We estimate that these outflows are potentially depositing energy into their host galaxies at considerable rates (~10^(43)-10^(45)erg/s); however, due to the lack of constraints on the density of the outflowing material and the structure of the outflow, these estimates should be taken as illustrative only. Based on the measured maximum velocities (v(max)~400-1400km/s) the outflows observed are likely to unbind some fraction of the gas from their host galaxies, but are unlikely to completely remove gas from the galaxy haloes. By using a combination of energetic arguments and a comparison to ULIRGs without clear evidence for AGN activity, we show that the AGN activity could be the dominant power source for driving all of the observed outflows, although star formation may also play a significant role in some of the sources.

preprint2012arXiv

GOODS-Herschel: radio-excess signature of hidden AGN activity in distant star-forming galaxies

We present here a new spectral energy distribution (SED) fitting approach that we adopt to select radio-excess sources amongst distant star-forming galaxies in the GOODS-Herschel (North) field and to reveal the presence of hidden, highly obscured AGN. Through extensive SED analysis of 458 galaxies with radio 1.4 GHz and mid-IR 24 um detections using some of the deepest Chandra X-ray, Spitzer and Herschel infrared, and VLA radio data available to date, we have robustly identified a sample of 51 radio-excess AGN (~1300 deg^-2) out to redshift z~3. These radio-excess AGN have a significantly lower far-IR/radio ratio (q<1.68) than the typical relation observed for star-forming galaxies (q~2.2). We find that ~45% of these radio-excess sources have a dominant AGN component in the mid-IR band, while for the remainders the excess radio emission is the only indicator of AGN activity. The fraction of radio-excess AGN increases with X-ray luminosity reaching ~60% at Lx~10^44-10^45 erg/s, making these sources an important part of the total AGN population. However, almost half (24/51) of these radio-excess AGN are not detected in the deep Chandra X-ray data, suggesting that some of these sources might be heavily obscured. We also find that the specific star formation rates (sSFRs) of the radio-excess AGN are on average lower that those observed for X-ray selected AGN hosts, indicating that our sources are forming stars more slowly than typical AGN hosts, and possibly their star formation is progressively quenching.

preprint2012arXiv

GOODS-Herschel: Ultra-deep XMM-Newton observations reveal AGN/star-formation connection

Models of galaxy evolution assume some connection between the AGN and star formation activity in galaxies. We use the multi-wavelength information of the CDFS to assess this issue. We select the AGNs from the 3Ms XMM-Newton survey and measure the star-formation rates of their hosts using data that probe rest-frame wavelengths longward of 20 um. Star-formation rates are obtained from spectral energy distribution fits, identifying and subtracting an AGN component. We divide the star-formation rates by the stellar masses of the hosts to derive specific star-formation rates (sSFR) and find evidence for a positive correlation between the AGN activity (proxied by the X-ray luminosity) and the sSFR for the most active systems with X-ray luminosities exceeding Lx=10^43 erg/s and redshifts z~1. We do not find evidence for such a correlation for lower luminosity systems or those at lower redshifts. We do not find any correlation between the SFR (or the sSFR) and the X-ray absorption derived from high-quality XMM-Newton spectra either, showing that the absorption is likely to be linked to the nuclear region rather than the host, while the star-formation is not nuclear. Comparing the sSFR of the hosts to the characteristic sSFR of star-forming galaxies at the same redshift we find that the AGNs reside mostly in main-sequence and starburst hosts, reflecting the AGN - sSFR connection. Limiting our analysis to the highest X-ray luminosity AGNs (X-ray QSOs with Lx>10^44 erg/s), we find that the highest-redshift QSOs (with z>2) reside predominantly in starburst hosts, with an average sSFR more than double that of the "main sequence", and we find a few cases of QSOs at z~1.5 with specific star-formation rates compatible with the main-sequence, or even in the "quiescent" region. (abridged)

preprint2012arXiv

No clear submillimetre signature of suppressed star formation amongst X-ray luminous AGNs

Many theoretical models require powerful active galactic nuclei (AGNs) to suppress star formation in distant galaxies and reproduce the observed properties of today's massive galaxies. A recent study based on Herschel-SPIRE submillimetre observations claimed to provide direct support for this picture, reporting a significant decrease in the mean star-formation rates (SFRs) of the most luminous AGNs (Lx>10^44 erg/s) at z=1-3 in the Chandra Deep Field-North (CDF-N). In this letter we extend these results using Herschel-SPIRE 250um data in the COSMOS and CDF-S fields to achieve an order of magnitude improvement in the number of sources at Lx>10^44 erg/s. On the basis of our analysis, we find no strong evidence for suppressed star formation in Lx>10^44 erg/s AGNs at z=1-3. The mean SFRs of the AGNs are constant over the broad X-ray luminosity range of Lx~10^43-10^45 erg/s (with mean SFRs consistent with typical star-forming galaxies at z~2; <SFRs>~100-200 Msol/yr). We suggest that the previous CDF-N results were likely due to low number statistics. We discuss our results in the context of current theoretical models.

preprint2012arXiv

The hidden "AGN main sequence": Evidence for a universal black hole accretion to star formation rate ratio since z~2 producing a M_BH-M* relation

Using X-ray stacking analyses we estimate the average amounts of supermassive black hole (SMBH) growth taking place in star-forming galaxies (SFGs) at z~1 and z~2 as a function of galaxy stellar mass (M*). We find the average SMBH growth rate follows remarkably similar trends with M* and redshift as the average star-formation rates (SFRs) of their host galaxies (i.e., (i.e., dM_BH/dt ~ M*^(0.86+/-0.39) for the z~1 sample and dM_BH/dt ~ M*^(1.05+/-0.36) for the z~2 sample). It follows that the ratio of SMBH growth rate to SFR is (a) flat with respect to M* (b) not evolving with redshift and (c) close to the ratio required to maintain/establish a SMBH to M* ratio of ~10^{-3} as also inferred from today's M_BH-M_bulge relationship. We interpret this as evidence that SMBHs have, on average, grown in-step with their host galaxies since at least z~2, irrespective of host galaxy mass and AGN triggering mechanism. As such, we suggest that the same secular processes that drive the bulk of star formation are also responsible for the majority of SMBH growth. From this, we speculate that it is the availability of gas reservoirs that regulate both cosmological SMBH growth and star formation.

preprint2011arXiv

Defining the intrinsic AGN infrared spectral energy distribution and measuring its contribution to the infrared output of composite galaxies

We use infrared spectroscopy and photometry to empirically define the intrinsic, thermal infrared spectral energy distribution (i.e., 6-100 um SED) of typical active galactic nuclei (i.e., 2-10 keV luminosity, Lx=10^{42}-10^{44} ergs/s AGNs). On average, the infrared SED of typical AGNs is best described as a broken power-law at <40 um that falls steeply at >40um (i.e., at far-infrared wavelengths). Despite this fall-off at long wavelengths, at least 3 of the 11 AGNs in our sample have observed SEDs that are AGN-dominated even at 60 um, demonstrating the importance of accounting for possible AGN contribution even at far-infrared wavelengths. Our results also suggest that the average intrinsic AGN 6-100 um SED gets bluer with increasing X-ray luminosity, a trend seen both within our sample and also when we compare against the intrinsic SEDs of more luminous quasars (i.e., Lx>10^{44} ergs/s). We compare our intrinsic AGN SEDs with predictions from dusty torus models and find they are more closely matched by clumpy, rather than continuous, torus models. Next, we use our intrinsic AGN SEDs to define a set of correction factors to convert either monochromatic infrared or X-ray luminosities into total intrinsic AGN infrared (i.e., 8-1000 um) luminosities. Finally, we outline a procedure that uses our newly defined intrinsic AGN infrared SEDs, in conjunction with a selection of host-galaxy templates, to fit the infrared photometry of composite galaxies and measure the AGN contribution to their total infrared output. We verify the accuracy of our SED fitting procedure by comparing our results to two independent measures of AGN contribution. Our SED fitting procedure opens up the possibility of measuring the intrinsic AGN luminosities of large numbers of galaxies with well-sampled infrared data (e.g., IRAS, ISO, Spitzer and Herschel).

preprint2011arXiv

GOODS-Herschel: The far-infrared view of star formation in AGN host galaxies since z~3

We present a study of the infrared properties of X-ray selected, moderate luminosity (Lx=10^{42}-10^{44}ergs/s) active galactic nuclei (AGNs) up to z~3, to explore the links between star formation in galaxies and accretion onto their central black holes. We use 100um and 160um fluxes from GOODS-Herschel -the deepest survey yet undertaken by the Herschel telescope- and show that in >94 per cent of cases these fluxes are dominated by the host. We find no evidence of any correlation between the X-ray and infrared luminosities of moderate AGNs at any redshift, suggesting that star-formation is decoupled from nuclear (AGN) activity. The star formation rates of AGN hosts increase strongly with redshift; by a factor of 43 from z<0.1 to z=2-3 for AGNs with the same X-ray luminosities. This increase is consistent with the factor of 25-50 increase in the specific star formation rates (SSFRs) of normal, star-forming (main-sequence) galaxies. Indeed, the average SSFRs of AGN hosts are only marginally (20 per cent) lower than those of main-sequence galaxies, with this small deficit being due to a fraction of AGNs residing in quiescent (low-SSFR) galaxies. We estimate 79+/-10 per cent of moderate AGNs are hosted in main-sequence galaxies, 15+/-7 per cent in quiescent galaxies and <10 per cent in strongly starbursting galaxies. The fractions of all main sequence galaxies at z<2 experiencing a period of moderate nuclear activity is strongly dependent on galaxy stellar mass (Mstars); rising from a few per cent at Mstars~10^{10}Msun to >20 per cent at Mstars>10^{11}Msun. Our results indicate that it is galaxy stellar mass that is most important in dictating whether a galaxy hosts a moderate luminosity AGN. We argue that the majority of moderate nuclear activity is fuelled by internal mechanisms rather than violent mergers, suggesting that disk instabilities could be an important AGN feeding mechanism.

preprint2011arXiv

X-ray Spectral Constraints for z~2 Massive Galaxies: The Identification of Reflection-Dominated Active Galactic Nuclei

We use the 4Ms CDF-S survey to place direct X-ray constraints on the ubiquity of z~2 heavily obscured AGNs in K<22 BzK galaxies. Forty seven of the 222 BzK galaxies in the central region of the CDF-S are detected at X-ray energies, 11 of which have hard X-ray spectral slopes (Gamma<1) indicating the presence of heavily obscured AGN activity. The other 36 X-ray detected BzK galaxies appear to be relatively unobscured AGNs and starburst galaxies; we use X-ray variability analyses over a rest-frame baseline of ~3 years to further confirm the presence of AGN activity in many of these systems. The majority (7 out of 11) of the heavily obscured AGNs have excess IR emission over that expected from star formation (termed "IR-excess galaxies"). However, we find that X-ray detected heavily obscured AGNs only comprise ~25% of the IR-excess galaxy population, which is otherwise composed of relatively unobscured AGNs and starburst galaxies. We find that the typical X-ray spectrum of the heavily obscured AGNs is better characterized by a pure reflection model than an absorbed power-law model, suggesting extreme Compton-thick absorption in some systems. We verify this result by producing a composite rest-frame 2-20 keV spectrum, which has a similar shape as a reflection-dominated X-ray spectrum and reveals an emission feature at rest-frame energy ~6.4 keV, likely to be due to Fe K. These heavily obscured AGNs are likely to be the distant analogs of the reflection-dominated AGNs recently identified at z~0 with >10 keV observatories. On the basis of these analyses we estimate the space density for typical (intrinsic X-ray luminosities of L_X>1E43 erg/s) heavily obscured and Compton-thick AGNs at z~2. Our space-density constraints are conservative lower limits but they are already consistent with the range of predictions from X-ray background models.

preprint2010arXiv

Highly Ionised Gas as a Diagnostic of the Inner NLR

The spectra of AGN from the ultraviolet to the near infrared, exhibit emission lines covering a wide range of ionisation states, from neutral species such as [O I] 6300A, up to [Fe XIV] 5303A. Here we report on some recent studies of the properties of highly ionised lines (HILs), plus two case studies of individual objects. Future IFU observations at high spatial and good spectral resolution, will probe the excitation and kinematics of the gas in the zone between the extended NLR and unresolved BLR. Multi-component SED fitting can be used to link the source of photoionisation with the strengths and ratios of the HILs.

preprint2010arXiv

The Synergy between Deep X-ray and Infrared Surveys: AGN and Star Formation Activity

We explore the connections between the infrared and X-ray properties of AGNs. Using the well constrained infrared SEDs of a sample of local (i.e., z < 0.1) sample of X-ray AGNs, we develop new diagnostics that exploit 24um and 70um flux densities to identify AGN-dominated systems at z<1.5 and measure their total infrared luminosities. We apply these diagnostics to the X-ray detected AGNs in the Chandra Deep Field South and, in doing so, reveal that the infrared to X-ray luminosity ratio was a factor of ~12 higher in the early Universe compared to today. We explore possible explanations for this dramatic evolution and demonstrate how forthcoming Herschel observations will discriminate between these scenarios and, while doing so, identify potential Compton-thick AGNs. Summary of Mullaney et al., MNRAS, 401, 995 (2010).

preprint2009arXiv

Characterising the Far-infrared Properties of Distant X-ray Detected AGNs: Evidence for Evolution in the Infrared--X-ray Luminosity Ratio

(Abridged) We investigate the far-infrared properties of X-ray sources detected in the Chandra Deep Field-South (CDF-S) survey using the ultra-deep 70um and 24um Spitzer observations taken in this field. We rely on stacking analyses of the 70um data to characterise the average 70um properties of the X-ray sources. Using Spitzer-IRS data of the Swift-BAT sample of z~0 active galactic nuclei (hereafter, AGNs), we show that the 70um/24um flux ratio can distinguish between AGN-dominated and starburst-dominated systems out to z~1.5. From stacking analysis we find that both high redshift and z~0 AGNs follow the same tendency toward warmer 70um/24um colours with increasing X-ray luminosity (L_X). We also show that the 70um flux can be used to determine the infrared (8-1000um) luminosities of high redshift AGNs. We use this information to show that L_X=10^{42-43} erg/s AGNs at high redshifts (z=1-2) have infrared to X-ray luminosity ratios (hereafter, L_IR/L_X) that are, on average, 4.7_{-2.0}^{+10.2} and 12.7^{+7.1}_{-2.6} times higher than AGNs with similar X-ray luminosities at z=0.5-1 and z~0, respectively. By comparison, we find that the L_IR/L_X ratios of L_X=10^{43-44} erg/s AGNs remain largely unchanged across this same redshift interval. We explore the consequences that these results may have on the identification of distant, potentially Compton thick AGNs using L_IR/L_X ratios. We discuss possible scenarios for the observed increase in the L_IR/L_X ratio with redshift, including changes in the dust covering factor of AGNs and/or the star formation rates of their host galaxies. Finally, we show how deep observations to be undertaken by the Herschel Space Observatory will enable us to discriminate between these proposed scenarios and also identify Compton-thick AGNs at high redshifts.

preprint2009arXiv

The Chandra Deep Protocluster Survey: Point-Source Catalogs for a 400 ks Observation of the z = 3.09 Protocluster in SSA22

We present X-ray point-source catalogs for a deep 400 ks Chandra ACIS-I exposure of the SSA22 field. The observations are centred on a z = 3.09 protocluster, which is populated by Lyman break galaxies (LBGs), Lyalpha emitters (LAEs), and extended Lyalpha-emitting blobs (LABs). The survey reaches ultimate (3 count) sensitivity limits of 5.7 X 10^-17 and 3.0 X 10^-16 ergs cm^-2 s^-1 for the 0.5-2 keV and 2-8 keV bands, respectively. These limits make SSA22 the fourth deepest extragalactic Chandra survey yet conducted, and the only one focused on a known high redshift structure. In total, we detect 297 X-ray point sources and identify one obvious bright extended X-ray source over a ~330 arcmin^2 region. In addition to our X-ray catalogs, we provide all available optical spectroscopic redshifts and near-infrared and mid-infrared photometry available for our sources. The basic X-ray and infrared properties of our Chandra sources indicate a variety of source types, although absorbed active galactic nuclei (AGNs) appear to dominate. In total, we have identified 12 X-ray sources (either via optical spectroscopic redshifts or LAE selection) at z =3.06-3.12 that are likely to be associated with the SSA22 protocluster. These sources have X-ray and multiwavelength properties that suggest they are powered by AGN with 0.5-8 keV luminosities in the range of 10^43-10^45 ergs s^-1 We have analysed the AGN fraction of sources in the protocluster as a function of local LAE source density and find suggestive evidence for a correlation between AGN fraction and local LAE source density (at the 96 per cent confidence level), implying that supermassive black hole growth at z = 3 is strongest in the highest density regions.