Source author record

C. M. Harrison

C. M. Harrison 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

31works
7topics
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

31 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.

preprint2022arXiv

Sonification and Sound Design for Astronomy Research, Education and Public Engagement

Over the last ten years there has been a large increase in the number of projects using sound to represent astronomical data and concepts. Motivation for these projects includes the potential to enhance scientific discovery within complex datasets, by utilising the inherent multi-dimensionality of sound and the ability of our hearing to filter signals from noise. Other motivations include creating engaging multi-sensory resources, for education and public engagement, and making astronomy more accessible to people who are blind or have low vision, promoting their participation in science and related careers. We describe potential benefits of sound within these contexts and provide an overview of the nearly 100 sound-based astronomy projects that we identified. We discuss current limitations and challenges of the approaches taken. Finally, we suggest future directions to help realise the full potential of sound-based techniques in general and to widen their application within the astronomy community.

preprint2022arXiv

The Quasar Feedback Survey: Revealing the Interplay of Jets, Winds & Emission Line Gas in Type 2 Quasars with Radio Polarization

We present results from a combined radio polarization and emission line study of five type 2 quasars at $z<0.2$ with the Karl G. Jansky Very Large Array (VLA) B-array at 5 GHz and Hubble Space Telescope (HST) [O III] observations. These five sources are known to exhibit close association between radio structures and ionized gas morphology and kinematics. Four sources (J0945+1737, J1000+1242, J1356+1026 and J1430+1339) show polarization in the current data. J1010+1413 is the unpolarized source in our sample. We detect $0.5-1\%$ fractional polarization in the radio cores and a high fractional polarization ($10-30\%$) in the lobes of these sources. The morphological, spectral and polarization properties suggest a jet origin for radio emission in J0945+1737, J1000+1242, J1010+1413 and J1430+1339 whereas the current data cannot fully discern the origin of radio emission (jet or wind) in J1356+1026. An anti-correlation between various polarized knots in the radio and [O III] emission is observed in our sources, similar to that observed in some radio-loud AGN in the literature. This suggests that the radio emission is likely to be depolarized by the emission-line gas. By modeling the depolarization effects, we estimate the size of the emission-line gas clouds to be $\sim(2.8\pm1.7)\times10^{-5}$ parsec and the amount of thermal material mixed with the synchrotron plasma to be $\sim(1.01\pm0.08)\times10^{6}$M$_\odot$ in the lobe of J0945+1737 (which exhibits the most prominent polarization signature in its lobe). The current work demonstrates that the interplay of jets/winds and emission-line gas is most likely responsible for the nature of radio outflows in radio-quiet AGN.

preprint2022arXiv

The resolved chemical abundance properties within the interstellar medium of star-forming galaxies at $\mathbf{ \textit{z} \approx 1.5}$

We exploit the unprecedented depth of integral field data from the KMOS Ultra-deep Rotational Velocity Survey (KURVS) to analyse the strong (H$α$) and forbidden ([NII], [SII]) emission line ratios in 22 main-sequence galaxies at $z\approx1.5$. Using the [NII]/H$α$ emission-line ratio we confirm the presence of the stellar mass $-$ gas-phase metallicity relation at this epoch, with galaxies exhibiting on average 0.13$\pm$0.04 dex lower gas-phase metallicity (12+log(O/H)$_{\rm M13}$=8.40$\pm$0.03) for a given stellar mass ($\log_{10}$($M_{\rm *}$[$M_{\odot}$]=10.1$\pm$0.1) than local main-sequence galaxies. We determine the galaxy-integrated [SII] doublet ratio, with a median value of [SII]$λ$6716/$λ$6731=1.26$\pm$0.14 equivalent to an electron density of log$_{10}$($n_{\rm e}$[cm$^{-3}$])=1.95$\pm$0.12. Utilising CANDELS $HST$ multi-band imaging we define the pixel surface-mass and star-formation rate density in each galaxy and spatially resolve the fundamental metallicity relation at $z\approx1.5$, finding an evolution of 0.05$\pm$0.01 dex compared to the local relation. We quantify the intrinsic gas-phase metallicity gradient within the galaxies using the [NII]/H$α$ calibration, finding a median annuli-based gradient of $Δ$Z/$Δ$R=$-$0.015$\pm$0.005 dex kpc$^{-1}$. Finally we examine the azimuthal variations in gas-phase metallicity, which show a negative correlation with the galaxy integrated star-formation rate surface density ($r_{\rm s}$=$-$0.40, $p_{\rm s}$=0.07) but no connection to the galaxies kinematic or morphological properties nor radial variations in stellar mass surface density or star formation rate surface density. This suggests both the radial and azimuthal variations in interstellar medium properties are connected to the galaxy integrated density of recent star formation.

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

Extended H$α$ over compact far-infrared continuum in dusty submillimeter galaxies -- Insights into dust distributions and star-formation rates at $z\sim2$

Using data from ALMA and near-infrared (NIR) integral field spectrographs including both SINFONI and KMOS on the VLT, we investigate the two-dimensional distributions of H$α$ and rest-frame far-infrared (FIR) continuum in six submillimeter galaxies at $z\sim2$. At a similar spatial resolution ($\sim$0.5" FWHM; $\sim$4.5 kpc at $z=2$), we find that the half-light radius of H$α$ is significantly larger than that of the FIR continuum in half of the sample, and on average H$α$ is a median factor of $2.0\pm0.4$ larger. Having explored various ways to correct for the attenuation, we find that the attenuation-corrected H$α$-based SFRs are systematically lower than the IR-based SFRs by at least a median factor of $3\pm1$, which cannot be explained by the difference in half-light radius alone. In addition, we find that in 40% of cases the total $V$-band attenuation ($A_V$) derived from energy balance modeling of the full ultraviolet(UV)-to-FIR spectral energy distributions (SEDs) is significantly higher than that derived from SED modeling using only the UV-to-NIR part of the SEDs, and the discrepancy appears to increase with increasing total infrared luminosity. Finally, considering all our findings along with the studies in the literature, we postulate that the dust distributions in SMGs, and possibly also in less IR luminous $z\sim2$ massive star-forming galaxies, can be decomposed into three main components; the diffuse dust heated by older stellar populations, the more obscured and extended young star-forming HII regions, and the heavily obscured central regions that have a low filling factor but dominate the infrared luminosity in which the majority of attenuation cannot be probed via UV-to-NIR emissions.

preprint2020arXiv

From Peculiar Morphologies to Hubble-type Spirals: The relation between galaxy dynamics and morphology in star-forming galaxies at z~1.5

We present an analysis of the gas dynamics of star-forming galaxies at z~1.5 using data from the KMOS Galaxy Evolution Survey (KGES). We quantify the morphology of the galaxies using $HST$ CANDELS imaging parametrically and non-parametrically. We combine the H$α$ dynamics from KMOS with the high-resolution imaging to derive the relation between stellar mass (M$_{*}$) and stellar specific angular momentum (j$_{*}$). We show that high-redshift star-forming galaxies at z~1.5 follow a power-law trend in specific stellar angular momentum with stellar mass similar to that of local late-type galaxies of the form j$_*$$\propto$M$_*^{0.53 \pm 0.10}$. The highest specific angular momentum galaxies are mostly disc-like, although generally, both peculiar morphologies and disc-like systems are found across the sequence of specific angular momentum at a fixed stellar mass. We explore the scatter within the j$_{*}$-M$_{*}$ plane and its correlation with both the integrated dynamical properties of a galaxy (e.g. velocity dispersion, Toomre Q$_{\rm g}$, H$α$ star formation rate surface density $Σ_{\rm SFR}$) and its parameterised rest-frame UV/optical morphology (e.g. Sérsic index, bulge to total ratio, Clumpiness, Asymmetry and Concentration). We establish that the position in the j$_{*}$-M$_{*}$ plane is strongly correlated with the star-formation surface density and the Clumpiness of the stellar light distribution. Galaxies with peculiar rest-frame UV/optical morphologies have comparable specific angular momentum to disc-dominated galaxies of the same stellar mass, but are clumpier and have higher star-formation rate surface densities. We propose that the peculiar morphologies in high--redshift systems are driven by higher star formation rate surface densities and higher gas fractions leading to a more clumpy inter-stellar medium.

preprint2020arXiv

KASHz: No evidence for ionised outflows instantaneously suppressing star formation in moderate luminosity AGN at $z$$\sim$$1.4$-$2.6$

As part of our KMOS AGN Survey at High-redshift (KASHz), we present spatially-resolved VLT/KMOS and VLT/SINFONI spectroscopic data and ALMA 870$μ$m continuum imaging of eight $z$=1.4--2.6 moderate AGN ($L_{\rm 2-10 \rm kev}$ = $10^{42} - 10^{45}$ ergs s$^{-1}$). We map [OIII], H$α$ and rest-frame FIR emission to search for any spatial anti-correlation between ionised outflows (traced by the [OIII] line) and star formation (SF; traced by H$α$ and FIR), that has previously been claimed for some high-z AGN and used as evidence for negative and/or positive AGN feedback. Firstly, we conclude that H$α$ is unreliable to map SF inside our AGN host galaxies based on: (i) SF rates inferred from attenuation-corrected H$α$ can lie below those inferred from FIR; (ii) the FIR continuum is more compact than the H$α$ emission by a factor of $\sim 2$ on average; (iii) in half of our sample, we observe significant spatial offsets between the FIR and H$α$ emission, with an average offset of $1.4\pm0.6$ kpc. Secondly, for the five targets with outflows we find no evidence for a spatial anti-correlation between outflows and SF using either H$α$ or FIR as a tracer. This holds for our re-analysis of a famous $z$=1.6 X-ray AGN (`XID 2028') where positive and negative feedback has been previously claimed. Based on our results, any impact on SF by ionised outflows must be subtle, either occurring on scales below our resolution, or on long timescales.

preprint2020arXiv

SUPER III. Broad Line Region properties of AGN at z$\sim$2

The SINFONI survey for Unveiling the Physics and Effect of Radiative feedback (SUPER) was designed to conduct a blind search for AGN-driven outflows on X-ray selected AGN at redshift z$\sim$2 with high ($\sim$2 kpc) spatial resolution, and correlate them to the properties of the host galaxy and central black hole. The main aims of this paper are: a) to derive reliable estimates for the BH mass and accretion rates for the Type-1 AGN in this survey; b) to characterize the properties of the AGN driven winds in the BLR. We analyzed rest-frame optical and UV spectra of 21 Type-1 AGN. We found that the BH masses estimated from H$α$ and H$β$ lines are in agreement. We estimate BH masses in the range Log(M$\rm_{BH}/M_{\odot}$)=8.4-10.8 and Eddington ratios $\rmλ_{Edd}$ =0.04-1.3. We confirm that the CIV line width does not correlate with the Balmer lines and the peak of the line profile is blue-shifted with respect to the [OIII]-based systemic redshift. These findings support the idea that the CIV line is tracing outflowing gas in the BLR, with velocities up to $\sim$4700 km/s. We confirm the strong dependence of the BLR wind velocity with the UV-to-Xray continuum slope, L$\rm_{Bol}$ and $\rmλ_{Edd}$. We inferred BLR mass outflow rates in the range 0.005-3 M$_{\odot}$/yr, showing a correlation with the bolometric luminosity consistent with that observed for ionized winds in the NLR and X-ray winds detected in local AGN, and kinetic power $\sim$10$^{[-7:-4]}\times$ L$\rm_{Bol}$. Finally, we found an anti-correlation between the equivalent width of the [OIII] line with respect to the CIV shift, and a positive correlation with [OIII] outflow velocity. These findings, for the first time in an unbiased sample of AGN at z$\sim$2, support a scenario where BLR winds are connected to galaxy scale detected outflows, and are capable of affecting the gas in the NLR located at kpc scale.

preprint2020arXiv

SUPER IV. CO(J=3-2) properties of active galactic nucleus hosts at cosmic noon revealed by ALMA

Feedback from AGN is thought to be key in shaping the life cycle of their host galaxies by regulating star-formation activity. Therefore, to understand the impact of AGN on star formation, it is essential to trace the molecular gas out of which stars form. In this paper we present the first systematic study of the CO properties of AGN hosts at z~2 for a sample of 27 X-ray selected AGN spanning two orders of magnitude in AGN bolometric luminosity (Lbol= 10^44.7-10^46.9 erg/s) by using ALMA Band 3 observations of the CO(3-2) transition (~1" angular resolution). To search for evidence of AGN feedback on the CO properties of the host galaxies, we compared our AGN with a sample of inactive (i.e., non-AGN) galaxies from the PHIBSS survey with similar redshift, stellar masses, and SFRs. We used the same CO transition as a consistent proxy for the gas mass for the two samples in order to avoid systematics involved when assuming conversion factors. By adopting a Bayesian approach to take upper limits into account, we analyzed CO luminosities as a function of stellar masses and SFRs, as well as the ratio LCO(3-2)/M* (proxy for the gas fraction). The two samples show statistically consistent trends in the LCO(3-2)-Lfir and LCO(3-2)-M* planes. However, there are indications that AGN feature lower CO(3-2) luminosities (0.4-0.7 dex) than inactive galaxies at the 2-3sigma level when we focus on the subset of parameters where the results are better constrained and on the distribution of the mean LCO(3-2)/M*. Therefore, even by conservatively assuming the same excitation factor r31, we would find lower molecular gas masses in AGN, and assuming higher r31 would exacerbate this difference. We interpret our result as a hint of the potential effect of AGN activity (e.g., radiation and outflows), which may be able to heat, excite, dissociate, and/or deplete the gas reservoir of the host galaxies. (abridged)

preprint2019arXiv

Do AGN really suppress star formation?

Active galactic nuclei (AGN) are believed to regulate star formation inside their host galaxies through "AGN feedback". We summarise our on-going study of luminous AGN (z~0.2-3; L_(AGN,bol)>~10^43 erg/s), which is designed to search for observational signatures of feedback by combining observed star-formation rate (SFR) measurements from statistical samples with cosmological model predictions. Using the EAGLE hydrodynamical cosmological simulations, in combination with our Herschel+ALMA surveys, we show that - even in the presence of AGN feedback - we do not necessarily expect to see any relationships between average galaxy-wide SFRs and instantaneous AGN luminosities. We caution that the correlation with stellar mass for both SFR and AGN luminosity can contribute to apparent observed positive trends between these two quantities. On the other hand, the EAGLE simulations, which reproduce our observations, predict that a signature of AGN feedback can be seen in the wide specific SFR distributions of $all$ massive galaxies (not just AGN hosts). Overall, whilst we can not rule out that AGN have an immediate small-scale impact on in-situ star-formation, all of our results are consistent with a feedback model where galaxy-wide in-situ star formation is not rapidly suppressed by AGN, but where the feedback likely acts over a longer timescale than a single AGN episode.

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 observations of Lyman-alpha Blob 1: halo sub-structure illuminated from within

We present new Atacama Large Millimeter/Submillimeter Array (ALMA) 850um continuum observations of the original Lyman-alpha Blob (LAB) in the SSA22 field at z=3.1 (SSA22-LAB01). The ALMA map resolves the previously identified submillimeter source into three components with total flux density S_850 = 1.68+/-0.06 mJy, corresponding to a star formation rate of ~150 M_sun/yr. The submillimeter sources are associated with several faint (m~27 mag) rest-frame ultraviolet sources identified in Hubble Space Telescope Imaging Spectrograph (STIS) clear filter imaging (~5850A). One of these companions is spectroscopically confirmed with Keck MOSFIRE to lie within 20 projected kpc and 250 km/s of one of the ALMA components. We postulate that some of these STIS sources represent a population of low-mass star-forming satellites surrounding the central submillimeter sources, potentially contributing to their growth and activity through accretion. Using a high resolution cosmological zoom simulation of a 10^13 M_sun halo at z=3, including stellar, dust and Ly-alpha radiative transfer, we can model the ALMA+STIS observations and demonstrate that Ly-alpha photons escaping from the central submillimeter sources are expected to resonantly scatter in neutral hydrogen, the majority of which is predicted to be associated with halo substructure. We show how this process gives rise to extended Ly-alpha emission with similar surface brightness and morphology to observed giant 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

The SCUBA-2 Cosmology Legacy Survey: 850um maps, catalogues and number counts

We present a catalogue of nearly 3,000 submillimetre sources detected at 850um over ~5 square degrees surveyed as part of the James Clerk Maxwell Telescope (JCMT) SCUBA-2 Cosmology Legacy Survey (S2CLS). This is the largest survey of its kind at 850um, probing a meaningful cosmic volume at the peak of star formation activity and increasing the sample size of submillimetre galaxies selected at 850um by an order of magnitude. We describe the wide 850um survey component of S2CLS, which covers the key extragalactic survey fields: UKIDSS-UDS, COSMOS, Akari-NEP, Extended Groth Strip, Lockman Hole North, SSA22 and GOODS-North. The average 1-sigma depth of S2CLS is 1.2 mJy/beam, approaching the SCUBA-2 850um confusion limit, which we determine to be ~0.8 mJy/beam. We measure the single dish 850um number counts to unprecedented accuracy, reducing the Poisson errors on the differential counts to approximately 4% at S_850~3mJy. With several independent fields, we investigate field-to-field variance, finding that the number counts on 0.5-1 degree scales are generally within 50% of the S2CLS mean for S_850>3mJy, with scatter consistent with the Poisson and estimated cosmic variance uncertainties, although there is a marginal (2-sigma) density enhancement in the GOODS-North field. The observed number counts are in reasonable agreement with recent phenomenological and semi-analytic models. Finally, the large solid angle of S2CLS allows us to measure the bright-end counts: at S_850>10mJy there are approximately ten sources per square degree, and we detect the distinctive up-turn in the number counts indicative of the detection of local sources of 850um emission and strongly lensed high-redshift galaxies. Here we describe the data collection and reduction procedures and present calibrated maps and a catalogue of sources; these are made publicly available.

preprint2016arXiv

The Spatially Resolved Dynamics of Dusty Starburst Galaxies in a z ~ 0.4 Cluster: Beginning the Transition from Spirals to S0s

To investigate what drives the reversal of the morphology-density relation at intermediate/high redshift, we present a multi-wavelength analysis of 27 dusty starburst galaxies in the massive cluster Cl 0024+17 at z = 0.4. We combine H-alpha dynamical maps from the VLT/FLAMES multi-IFU system with far-infrared imaging using Herschel SPIRE and millimetre spectroscopy from IRAM/NOEMA, in order to measure the dynamics, star formation rates and gas masses of this sample. Most galaxies appear to be rotationally supported, with a median ratio of rotational support to line-of-sight velocity dispersion v/sigma ~ 5 +/- 2, and specific angular momentum lambda_R = 0.83 +/- 0.06 - comparable to field spirals of a similar mass at this redshift. The star formation rates of 3 - 26 M_solar/yr and average 12 CO derived gas mass of 1 x 10^10 M_solar suggest gas depletion timescales of ~ 1Gyr (~ 0.25 of the cluster crossing time). We derive characteristic dust temperatures (mean T_dust = 26 +/- 1 K) consistent with local galaxies of similar far-infrared luminosity, suggesting that the low density gas is yet to be stripped. Taken together, these results suggest that these starbursts have only recently accreted from the field, with star formation rates likely enhanced due to the effects of ram pressure. In order to make the transition to cluster S0s these galaxies must lose ~ 40% of their specific angular momentum. We suggest this must occur > 1 Gyr later, after the molecular gas has been depleted and/or stripped, via multiple tidal interactions with other cluster members.

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

Cold Dust Emission from X-ray AGN in the SCUBA-2 Cosmology Legacy Survey: Dependence on Luminosity, Obscuration & AGN Activity

We study the 850um emission in X-ray selected AGN in the 2 sq-deg COSMOS field using new data from the SCUBA-2 Cosmology Legacy Survey. We find 19 850um bright X-ray AGN in a high-sensitivity region covering 0.89 sq-deg with flux densities of S850=4-10 mJy. The 19 AGN span the full range in redshift and hard X-ray luminosity covered by the sample - 0.7<z<3.5 and 43.2<log10(LX) <45. We report a highly significant stacked 850um detection of a hard X-ray flux-limited population of 699 z>1 X-ray AGN - S850=0.71+/-0.08mJy. We explore trends in the stacked 850um flux densities with redshift, finding no evolution in the average cold dust emission over the redshift range probed. For Type 1 AGN, there is no significant correlation between the stacked 850um flux and hard X-ray luminosity. However, in Type 2 AGN the stacked submm flux is a factor of 2 higher at high luminosities. When averaging over all X-ray luminosities, no significant differences are found in the stacked submm fluxes of Type 1 and Type 2 AGN as well as AGN separated on the basis of X-ray hardness ratios and optical-to-infrared colours. However, at log10(LX) >44.4, dependences in average submm flux on the optical-to-infrared colours become more pronounced. We argue that these high luminosity AGN represent a transition from a secular to a merger-driven evolutionary phase where the star formation rates and accretion luminosities are more tightly coupled. Stacked AGN 850um fluxes are compared to the stacked fluxes of a mass-matched sample of K-band selected non-AGN galaxies. We find that at 10.5<log10(M*/M0)<11.5, the non-AGN 850um fluxes are 1.5-2x higher than in Type 2 AGN of equivalent mass. We suggest these differences are due to the presence of massive dusty, red starburst galaxies in the K-band selected non-AGN sample, which are not present in optically selected catalogues covering a smaller area.

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

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.

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

The SCUBA-2 Cosmology Legacy Survey: the submillimetre properties of Lyman break galaxies at z=3-5

We present statistically significant detections at 850um of the Lyman Break Galaxy (LBG) population at z=3, 4, and 5 using data from the Submillimetre Common User Bolometer Array 2 (SCUBA-2) Cosmology Legacy Survey (S2CLS) in the United Kingdom Infrared Deep Sky Survey Ultra Deep Survey (UKIDSS-UDS) field. We employ a stacking technique to probe beneath the survey limit to measure the average 850um flux density of LBGs at z=3, 4, and 5 with typical ultraviolet luminosities of L(1700A)~10^29 erg/s/Hz. We measure 850um flux densities of (0.25 +/- 0.03, (0.41 +/- 0.06), and (0.88 +/- 0.23) mJy respectively, and find that they contribute at most 20 per cent to the cosmic far-infrared background at 850um. Fitting an appropriate range of spectral energy distributions to the z=3, 4, and 5 LBG stacked 24-850um fluxes, we derive infrared (IR) luminosities of L(8-1000um)~3.2, 5.5, and 11.0x10^11 Lsun (corresponding to star formation rates of ~50-200 Msun/yr) respectively. We find that the evolution in the IR luminosity density of LBGs is broadly consistent with model predictions for the expected contribution of luminous IR galaxy (LIRG) to ultraluminous IR galaxy (ULIRG) type systems at these epochs. We also see a strong positive correlation between stellar mass and IR luminosity. Our data are consistent with the main sequence of star formation showing little or no evolution from z=3 to 5. We have also confirmed that, for a fixed mass, the reddest LBGs (UV slope Beta -> 0) are indeed redder due to dust extinction, with SFR(IR)/SFR(UV) increasing by approximately an order of magnitude over -2<Beta<0 such that SFR(IR)/SFR(UV)~20 for the reddest LBGs. Furthermore, the most massive LBGs also tend to have higher obscured-to-unobscured ratio, hinting at a variation in the obscuration properties across the mass range.

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

The impact of AGN on their host galaxies

In these proceedings I briefly: (1) review the impact (or "feedback") that active galactic nuclei (AGN) are predicted to have on their host galaxies and larger scale environment, (2) review the observational evidence for or against these predictions and (3) present new results on ionised outflows in AGN. The observational support for the "maintenance mode" of feedback is strong (caveat the details); AGN at the centre of massive halos appear to be regulating the cooling of hot gas, which could in turn control the levels of future star formation (SF) and black hole growth. In contrast, direct observational support for more rapid forms of feedback, which dramatically impact on SF (i.e., the "quasar mode"), remains elusive. From a systematic study of the spectra of approx. 24000 z<0.4 AGN we find that extreme ionised gas kinematics are common, and are most prevalent in radio bright AGN (L[1.4GHz] >10^23 W/Hz). Follow-up IFU observations have shown that these extreme gas kinematics are extended over kilo-parsec scales. However, the co-existence of high-levels of SF, luminous AGN activity and radio jets raises interesting questions on the primary drivers and impact of these outflows. Galaxy-wide, high-mass outflows are being observed in an increasing number of AGN and are a plausible mechanism for the depletion of gas; however, there is still much work to be done to determine the physical processes that drive these outflows and to measure the level of impact that they have on their host galaxies.

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

Integral Field Spectroscopy of 2.0<z<2.7 Sub-mm Galaxies; gas morphologies and kinematics

We present two-dimensional, integral field spectroscopy covering the rest-frame wavelengths of strong optical emission lines in nine sub-mm-luminous galaxies (SMGs) at 2.0<z<2.7. The GEMINI-NIFS and VLT-SINFONI imaging spectroscopy allows the mapping of the gas morphologies and dynamics within the sources, and we measure an average Halpha velocity dispersion of sigma=220+-80km/s and an average half light radius of r=3.7+-0.8kpc. The average dynamical measure, V_obs/2sigma=0.9+-0.1 for the SMGs, is higher than in more quiescent star-forming galaxies at the same redshift, highlighting a difference in the dynamics of the two populations. The SMGs' far-infrared SFRs, measured using Herschel-SPIRE far-infrared photometry, are on average 370+-90Mo/yr which is ~2 times higher than the extinction corrected SFRs of the more quiescent star-forming galaxies. Six of the SMGs in our sample show strong evidence for kinematically distinct multiple components with average velocity offsets of 200+-100km/s and average projected spatial offsets of 8+-2kpc, which we attribute to systems in the early stages of major mergers. Indeed all SMGs are classified as mergers from a kinemetry analysis of the velocity and dispersion field asymmetry. We bring together our sample with seven other SMGs with IFU observations to describe the ionized gas morphologies and kinematics in a sample of 16 SMGs. By comparing the velocity and spatial offsets of the SMG Halpha components with sub-halo offsets in the Millennium simulation database we infer an average halo mass for SMGs of 13<log(M[h^-1Mo])<14. Finally we explore the relationship between the velocity dispersion and star formation intensity within the SMGs, finding the gas motions are consistent with the Kennicutt-Schmidt law and a range of extinction corrections, although might also be driven by the tidal torques from merging or even the star formation itself.

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 evolutionary connection between QSOs and SMGs: molecular gas in far-infrared luminous QSOs at z ~ 2.5

We present IRAM Plateau de Bure Interferometer observations of the 12CO(3-2) emission from two far-infrared luminous QSOs at z ~ 2.5 selected from the Herschel-ATLAS survey. These far-infrared bright QSOs were selected to have supermassive black holes (SMBHs) with masses similar to those thought to reside in sub-millimetre galaxies (SMGs) at z ~ 2.5; making them ideal candidates as systems in transition from an ultraluminous infrared galaxy phase to a sub-mm faint, unobscured, QSO. We detect 12CO(3-2) emission from both QSOs and we compare their baryonic, dynamical and SMBH masses to those of SMGs at the same epoch. We find that these far-infrared bright QSOs have similar dynamical but lower gas masses than SMGs. In particular we find that far-infrared bright QSOs have ~50+-23% less warm/dense gas than SMGs, which combined with previous results showing the QSOs lack the extended, cool reservoir of gas seen in SMGs, suggests that they are at a different evolutionary stage. This is consistent with the hypothesis that far-infrared bright QSOs represent a short (~1Myr) but ubiquitous phase in the transformation of dust obscured, gas-rich, starburst-dominated SMGs into unobscured, gas-poor, QSOs.