We investigate the properties of a sample of 35 galaxies, detected with ALMA at 1.1 mm in the GOODS-ALMA field (area of 69 arcmin$^2$, resolution = 0.60", RMS $\simeq$ 0.18 mJy beam$^{-1}$). Using the UV-to-radio deep multiwavelength coverage of the GOODS-South field, we fit the spectral energy distributions of these galaxies to derive their key physical properties. The galaxies detected by ALMA are among the most massive at $z$ = 2-4 (M$_{\star,med}$ = 8.5$ \times$ 10$^{10}$ M$_\odot$) and are either starburst or located in the upper part of the galaxy star-forming main sequence. A significant portion of our galaxy population ($\sim$ 40%), located at $z\sim$ 2.5-3, exhibits abnormally low gas fractions. The sizes of these galaxies, measured with ALMA, are compatible with the trend between $H$-band size and stellar mass observed for $z\sim2$ elliptical galaxies suggesting that they are building compact bulges. We show that there is a strong link between star formation surface density (at 1.1 mm) and gas depletion time: the more compact a galaxy's star-forming region is, the shorter its lifetime will be (without gas replenishment). The identified compact sources associated with relatively short depletion timescales ($\sim$100 Myr), are the ideal candidates to be the progenitors of compact elliptical galaxies at $z$ $\sim$ 2.

In this paper, we extend the source detection in the GOODS-ALMA field (69 arcmin$^2$, rms sensitivity $\sigma$ $\simeq$ 0.18 mJy.beam$^{-1}$), to deeper levels than presented in Franco et al. (2018). Using positional information at 3.6 and 4.5 $\mu$m (from Spitzer-IRAC), we explore the presence of galaxies detected at 1.1 mm with ALMA below our original blind detection limit of 4.8-$\sigma$ at which the number of spurious sources starts to dominate over that of real sources. In this Supplementary Catalog, we find a total of 16 galaxies, including 2 galaxies with no counterpart in HST images (also known as optically-dark galaxies) down to a 5$\sigma$ limiting depth of H = 28.2 AB (HST/WFC3 F160W). This brings the total sample of GOODS-ALMA 1.1 mm sources to 35 galaxies. Galaxies in the new sample cover a wider dynamic range in redshift ($z$ = 0.65 - 4.73), are on average twice as large (1.3 vs 0.65 kpc) and and have lower stellar mass (M$_{\star}^{\rm SC}$ = 7.6$\times$10$^{10}$M$_\odot$ vs M$_{\star}^{\rm MC}$ = 1.2$\times$10$^{11}$M$_\odot$). Although exhibiting larger physical sizes, these galaxies have still far-infrared sizes significantly more compact than inferred from their optical emission. We show that the astrometry of the HST image does not only suffer from a global astrometric shift, as already discussed in previous papers, but also from local shifts. These distortions were artificially introduced in the process of building the mosaic of the GOODS-South HST image. By comparing the positions of almost 400 galaxies detected by HST, Pan-STARRS and ALMA, we create a distortion map which can be used to correct for these astrometric issues.

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.

We report the results of optical follow-up observations of 29 gravitational-wave triggers during the first half of the LIGO-Virgo Collaboration (LVC) O3 run with the Gravitational-wave Optical Transient Observer (GOTO) in its prototype 4-telescope configuration (GOTO-4). While no viable electromagnetic counterpart candidate was identified, we estimate our 3D (volumetric) coverage using test light curves of on- and off-axis gamma-ray bursts and kilonovae. In cases where the source region was observable immediately, GOTO-4 was able to respond to a GW alert in less than a minute. The average time of first observation was $8.79$ hours after receiving an alert ($9.90$ hours after trigger). A mean of $732.3$ square degrees were tiled per event, representing on average $45.3$ per cent of the LVC probability map, or $70.3$ per cent of the observable probability. This coverage will further improve as the facility scales up alongside the localisation performance of the evolving gravitational-wave detector network. Even in its 4-telescope prototype configuration, GOTO is capable of detecting AT2017gfo-like kilonovae beyond 200 Mpc in favourable observing conditions. We cannot currently place meaningful electromagnetic limits on the population of distant ($\hat{D}_L = 1.3$ Gpc) binary black hole mergers because our test models are too faint to recover at this distance. However, as GOTO is upgraded towards its full 32-telescope, 2 node (La Palma & Australia) configuration, it is expected to be sufficiently sensitive to cover the predicted O4 binary neutron star merger volume, and will be able to respond to both northern and southern triggers.

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 ($\lambda_{\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 $\lambda_{\rm EDD}$ distribution, and with a positive shift of the break $\lambda^{*}$, 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 ($\lambda_{\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 $\lambda_{\rm EDD}$ prescriptions for understanding the cosmic BH mass assembly.

At a fraction the total cost of an equivalent orbital mission, scientific balloon-borne platforms, operating above 99.7% of the Earth's atmosphere, offer attractive, competitive, and effective observational capabilities -- namely space-like resolution, transmission, and backgrounds -- that are well suited for modern astronomy and cosmology. SuperBIT is a diffraction-limited, wide-field, 0.5 m telescope capable of exploiting these observing conditions in order to provide exquisite imaging throughout the near-IR to near-UV. It utilizes a robust active stabilization system that has consistently demonstrated a 1 sigma sky-fixed pointing stability at 48 milliarcseconds over multiple 1 hour observations at float. This is achieved by actively tracking compound pendulations via a three-axis gimballed platform, which provides sky-fixed telescope stability at < 500 milliarcseconds and corrects for field rotation, while employing high-bandwidth tip/tilt optics to remove residual disturbances across the science imaging focal plane. SuperBIT's performance during the 2019 commissioning flight benefited from a customized high-fidelity science-capable telescope designed with exceptional thermo- and opto-mechanical stability as well as tightly constrained static and dynamic coupling between high-rate sensors and telescope optics. At the currently demonstrated level of flight performance, SuperBIT capabilities now surpass the science requirements for a wide variety of experiments in cosmology, astrophysics and stellar dynamics.

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 $\epsilon \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}$).

We infer and compare the specific X-ray luminosity distributions for a sample of massive (i.e. $\log_{10} (M*/M\odot) > 10.5$) galaxies split according to their far-infrared-derived star-forming properties (i.e., starburst and non-starburst) and redshift. We model each distribution as a power-law with an upper and lower turnover, and adopt a maximum likelihood method to include information from non-detections in the form of upper limits. When we use our inferred distributions to calculate the ratios of high to low sLx AGN (corresponding to above and below $0.1\lambda_{\text{Edd}}$, respectively) we find that starbursts have significantly higher proportions of high sLx AGN compared to their non-starburst counterparts. These findings help explain the increase in average X-ray luminosity in bins of increasing SFR reported by previous studies.

We present 1-7 GHz high-resolution radio imaging (VLA and e-MERLIN) and spatially-resolved ionized gas kinematics for ten z<0.2 type~2 `obscured' quasars (log [L(AGN)/(erg/s)]>~45) with moderate radio luminosities (log [L(1.4GHz)/(W/Hz)]=23.3-24.4). These targets were selected to have known ionized outflows based on broad [OIII] emission-line components (FWHM~800-1800 km/s). Although `radio-quiet' and not `radio AGN' by many traditional criteria, we show that for nine of the targets, star formation likely accounts for <~10 per cent of the radio emission. We find that ~80-90 per cent of these nine targets exhibit extended radio structures on 1-25 kpc scales. The quasars' radio morphologies, spectral indices and position on the radio size-luminosity relationship reveals that these sources are consistent with being low power compact radio galaxies. Therefore, we favour radio jets as dominating the radio emission in the majority of these quasars. The radio jets we observe are associated with morphologically and kinematically distinct features in the ionized gas, such as increased turbulence and outflowing bubbles, revealing jet-gas interaction on galactic scales. Importantly, such conclusions could not have been drawn from current low-resolution radio surveys such as FIRST. Our observations support a scenario where compact radio jets, with modest radio luminosities, are a crucial feedback mechanism for massive galaxies during a quasar phase.

We explore the distribution of RMS=SFR/SFR_MS (where SFR_MS is the star formation rate of "Main Sequence" star-forming galaxies) for AGN hosts at z=1. We split our sample into two bins of X-ray luminosity divided at Lx=2x10^43erg s-1 to investigate whether the RMS distribution changes as a function of AGN power. Our main results suggest that, when the RMS distribution of AGN hosts is modelled as a log-normal distribution (i.e. the same shape as that of MS galaxies), galaxies hosting more powerful X-ray AGNs (i.e. Lx>2x10^43erg s-1) display a narrower RMS distribution that is shifted to higher values compared to their lower Lx counterparts. In addition, we find that more powerful X-ray AGNs have SFRs that are more consistent with that of MS galaxies compared to lower Lx AGNs. Despite this, the mean SFRs (as opposed to RMS) measured from these distributions are consistent with the previously observed flat relationship between SFR and Lx. Our results suggest that the typical star-forming properties of AGN hosts change with Lx , and that more powerful AGNs typically reside in more MS-like star-forming galaxies compared to lower Lx AGNs.

We present the X-UDS survey, a set of wide and deep Chandra observations of the Subaru-XMM Deep/UKIDSS Ultra Deep Survey (SXDS/UDS) field. The survey consists of 25 observations that cover a total area of 0.33 deg$^{2}$. The observations are combined to provide a nominal depth of ~600 ksec in the central 100 arcmin$^{2}$ region of the field that has been imaged with Hubble/WFC3 by the CANDELS survey and $\sim$200 ksec in the remainder of the field. In this paper, we outline the survey's scientific goals, describe our observing strategy, and detail our data reduction and point source detection algorithms. Our analysis has resulted in a total of 868 band-merged point sources detected with a false-positive Poisson probability of $<1\times10^{-4}$. In addition, we present the results of an X-ray spectral analysis and provide best-fitting neutral hydrogen column densities, $N_{\rm H}$, as well as a sample of 51 Compton-thick active galactic nucleus candidates. Using this sample, we find the intrinsic Compton-thick fraction to be 30-35% over a wide range in redshift ($z=0.1-3$), suggesting the obscured fraction does not evolve very strongly with epoch. However, if we assume that the Compton-thick fraction is dependent on luminosity, as is seen for Compton-thin sources, then our results are consistent with a rise in the obscured fraction out to $z\sim3$. Finally, an examination of the host morphologies of our Compton-thick candidates shows a high fraction of morphological disturbances, in agreement with our previous results. All data products described in this paper are made available via a public website.

We present a 69 arcmin$^2$ ALMA survey at 1.1mm, GOODS-ALMA, matching the deepest HST-WFC3 H-band part of the GOODS-South field. We taper the 0"24 original image with a homogeneous and circular synthesized beam of 0"60 to reduce the number of independent beams - thus reducing the number of purely statistical spurious detections - and optimize the sensitivity to point sources. We extract a catalogue of galaxies purely selected by ALMA and identify sources with and without HST counterparts down to a 5$\sigma$ limiting depth of H=28.2 AB (HST/WFC3 F160W). ALMA detects 20 sources brighter than 0.7 mJy in the 0"60 tapered mosaic (rms sensitivity =0.18 mJy/beam) with a purity greater than 80%. Among these detections, we identify three sources with no HST nor Spitzer-IRAC counterpart, consistent with the expected number of spurious galaxies from the analysis of the inverted image; their definitive status will require additional investigation. An additional three sources with HST counterparts are detected either at high significance in the higher resolution map, or with different detection-algorithm parameters ensuring a purity greater than 80%. Hence we identify in total 20 robust detections. Our wide contiguous survey allows us to push further in redshift the blind detection of massive galaxies with ALMA with a median redshift of $z$=2.92 and a median stellar mass of M$_{\star}$ = 1.1 $\times 10^{11}$M$_\odot$. Our sample includes 20% HST-dark galaxies (4 out of 20), all detected in the mid-infrared with IRAC. The near-infrared based photometric redshifts of two of them $z\sim$4.3 and 4.8) suggest that these sources have redshifts $z$>4. At least 40% of the ALMA sources host an X-ray AGN, compared to 14% for other galaxies of similar mass and redshift. The wide area of our ALMA survey provides lower values at the bright end of number counts than single-dish telescopes

We present the star formation rates (SFRs) of a sample of 109 galaxies with X-ray selected active galactic nuclei (AGN) with moderate to high X-ray luminosities (L(2-8keV)= 10^42-10^45 erg/s), at redshifts 1 < z < 4.7, that were selected to be faint or undetected in the Herschel bands. We combine our deep ALMA continuum observations with deblended 8-500{\mu}m photometry from Spitzer and Herschel, and use infrared (IR) SED fitting and AGN - star formation decomposition methods. The addition of the ALMA photometry results in an order of magnitude more X-ray AGN in our sample with a measured SFR (now 37%). The remaining 63% of the sources have SFR upper limits that are typically a factor of 2-10 times lower than the pre-ALMA constraints. With the improved constraints on the IR SEDs, we can now identify a mid-IR (MIR) AGN component in 50% of our sample, compared to only ~1% previously. We further explore the F870{\mu}m/F24{\mu}m-redshift plane as a tool for the identification of MIR emitting AGN, for three different samples representing AGN dominated, star formation dominated, and composite sources. We demonstrate that the F870{\mu}m/F24{\mu}m-redshift plane can successfully split between AGN and star formation dominated sources, and can be used as an AGN identification method.

We present the X-ray properties of the 'Teacup AGN' (SDSS J1430+1339), a $z=0.085$ type 2 quasar which is interacting dramatically with its host galaxy. Spectral modelling of the central quasar reveals a powerful, highly obscured AGN with a column density of $N_{\rm H}=(4.2$-$6.5)\times 10^{23}$ cm$^{-2}$ and an intrinsic luminosity of $L_{\rm 2\mbox{-}10\,keV}=(0.8$-$1.4)\times 10^{44}$ erg s$^{-1}$. The current high bolometric luminosity inferred ($L_{\rm bol}\approx 10^{45}$-$10^{46}$ erg s$^{-1}$) has ramifications for previous interpretations of the Teacup as a fading/dying quasar. High resolution Chandra imaging data reveal a $\approx 10$ kpc loop of X-ray emission, co-spatial with the 'eastern bubble' previously identified in luminous radio and ionised gas (e.g., [OIII] line) emission. The X-ray emission from this structure is in good agreement with a shocked thermal gas, with $T=(4$-$8)\times 10^{6}$ K, and there is evidence for an additional hot component with $T\gtrsim 3\times 10^{7}$ K. Although the Teacup is a radiatively dominated AGN, the estimated ratio between the bubble power and the X-ray luminosity is in remarkable agreement with observations of ellipticals, groups, and clusters of galaxies undergoing AGN feedback.

We explore the kinematics of ionized gas via the [O III] $\lambda$5007 emission lines in active galactic nuclei (AGN) selected on the basis of their mid-infrared (IR) emission, and split into obscured and unobscured populations based on their optical-IR colors. After correcting for differences in redshift distributions, we provide composite spectra of spectroscopically and photometrically defined obscured/Type 2 and unobscured/Type 1 AGN from 3500 to 7000 \AA. The IR-selected obscured sources contain a mixture of narrow-lined Type 2 AGN and intermediate sources that have broad H$\alpha$ emission and significantly narrower H$\beta$. Using both [OIII] luminosities and AGN luminosities derived from optical-IR spectral energy distribution fitting, we find evidence for enhanced large-scale obscuration in the obscured sources. In matched bins of luminosity we find that the obscured population typically has broader, more blueshifted \OIII\ emission than in the unobscured sample, suggestive of more powerful AGN-driven outflows. This trend is not seen in spectroscopically classified samples, and is unlikely to be entirely explained by orientation effects. In addition, outflow velocities increase from small to moderate AGN $E(B-V)$ values, before flattening out (as traced by FWHM) and even decreasing (as traced by blueshift). While difficult to fully interpret in a single physical model, due to both the averaging over populations and the spatially-averaged spectra, these results agree with previous findings that simple geometric unification models are insufficient for the IR-selected AGN population, and may fit into an evolutionary model for obscured and unobscured AGN.

The lack of a strong correlation between AGN X-ray luminosity ($L_X$; a proxy for AGN power) and the star formation rate (SFR) of their host galaxies has recently been attributed to stochastic AGN variability. Studies using population synthesis models have incorporated this by assuming a broad, universal (i.e. does not depend on the host galaxy properties) probability distribution for AGN specific X-ray luminosities (i.e. the ratio of $L_X$ to host stellar mass; a common proxy for Eddington ratio). However, recent studies have demonstrated that this universal Eddington ratio distribution fails to reproduce the observed X-ray luminosity functions beyond z$\sim$1.2. Furthermore, empirical studies have recently shown that the Eddington ratio distribution may instead depend upon host galaxy properties, such as SFR and/or stellar mass. To investigate this further we develop a population synthesis model in which the Eddington ratio distribution is different for star-forming and quiescent host galaxies. We show that, although this model is able to reproduce the observed X-ray luminosity functions out to z$\sim$2, it fails to simultaneously reproduce the observed flat relationship between SFR and X-ray luminosity. We can solve this, however, by incorporating a mass dependency in the AGN Eddington ratio distribution for star-forming host galaxies. Overall, our models indicate that a relative suppression of low Eddington ratios ($\lambda_{\rm Edd}\lesssim$0.1) in lower mass galaxies (M<$10^{10-11}$Msun) is required to reproduce both the observed X-ray luminosity functions and the observed flat SFR/X-ray relationship.

We discuss the spectral analysis of a sample of 63 Active Galactic Nuclei (AGN) detected above a limiting flux of $S(8-24$ keV$)=7\times10^{-14}$ erg/s/cm$^2$ in the multi-tiered NuSTAR Extragalactic Survey program. The sources span a redshift range z=0-2.1 (median $\langle$z$\rangle=$0.58). The spectral analysis is performed over the broad 0.5-24 keV energy range, combining NuSTAR with Chandra and/or XMM-Newton data and employing empirical and physically motivated models. This constitutes the largest sample of AGN selected at $ >10$ keV to be homogeneously spectrally analyzed at these flux levels. We study the distribution of spectral parameters such as photon index, column density ($N_{\rm H}$), reflection parameter ($R$) and 10-40 keV luminosity ($L_{X}$). Heavily obscured ($\log [N_{\rm H}/cm^{-2}]\ge23$) and Compton Thick (CT; $\log [N_{\rm H}/cm^{-2}]\ge24$) AGN constitute $\sim$25% (15-17 sources) and $\sim$2-3% ( 1-2 sources) of the sample, respectively. The observed $N_{\rm H}$ distribution fairly agrees with predictions of Cosmic X-ray Background population synthesis models (CXBPSM). We estimate the intrinsic fraction of AGN as a function of $N_{H}$, accounting for the bias against obscured AGN in a flux-selected sample. The fraction of CT AGN relative to $\log [N_{\rm H}/cm^{-2}]=20-24$ AGN is poorly constrainted, formally in the range 2-56% (90% upper limit of 66%). We derived a fraction ($f_{abs}$) of obscured AGN ($\log [N_{\rm H}/cm^{-2}]=22-24$) as a function of $L_{X}$ in agreement with CXBPSM and previous z < 1 X-ray determinations. Furthermore $f_{abs}$ at z=0.1-0.5 and $\log (L_{X}/erg$ $s^{-1})\approx43.6-44.3$ agrees with observational measurements/trends obtained over larger redshift intervals. We report a significant anti-correlation of $R$ with $L_{X}$ (confirmed by our companion paper on stacked spectra) with considerable scatter around the median $R$ values.

We present sensitive 870$\mu$m continuum measurements from our ALMA programmes of 114 X-ray selected AGN in the CDF-S and COSMOS fields. We use these observations in combination with data from Spitzer and Herschel to construct a sample of 86 X-ray selected AGN, 63 with ALMA constraints at $z=1.5-3.2$ with stellar mass $>2\times10^{10}M_{\odot}$. We constructed broad-band spectral energy distributions in the infrared band (8 - 1000$\mu$m) and constrain star-formation rates (SFRs) uncontaminated by the AGN. Using a hierarchical Bayesian method that takes into account the information from upper limits, we fit SFR and specific SFR (sSFR) distributions. We explore these distributions as a function of both X-ray luminosity and stellar mass. We compare our measurements to two versions of the EAGLE hydrodynamical simulations: the reference model with AGN feedback and the model without AGN. We find good agreement between the observations and that predicted by the EAGLE reference model for the modes and widths of the sSFR distributions as a function of both X-ray luminosity and stellar mass; however, we found that the EAGLE model without AGN feedback predicts a significantly narrower width when compared to the data. Overall, from the combination of the observations with the model predictions, we conclude that (1) even with AGN feedback, we expect no strong relationship between the sSFR distribution parameters and instantaneous AGN luminosity and (2) a signature of AGN feedback is a broad distribution of sSFRs for all galaxies (not just those hosting an AGN) with stellar masses above $\approx 10^{10}$M$_{\odot}$.

We present a study of the average X-ray spectral properties of the sources detected by the NuSTAR extragalactic survey, comprising observations of the E-CDFS, EGS and COSMOS fields. The sample includes 182 NuSTAR sources (64 detected at 8-24 keV), with 3-24 keV fluxes ranging between $f_{\rm 3-24 keV}\approx10^{-14}$ and $6\times10^{-13}$ erg/cm$^2$/s ($f_{\rm 8-24 keV}\approx3\times10^{-14}-3\times10^{-13}$ erg/cm$^2$/s) and redshifts of $z=0.04-3.21$. We produce composite spectra from the Chandra+NuSTAR data ($E\approx2-40$ keV, rest frame) for all the sources with redshift identifications (95%) and investigate the intrinsic, average spectra of the sources, divided into broad-line (BL) and narrow-line (NL) AGN, and also in different bins of X-ray column density and luminosity. The average power-law photon index for the whole sample is $\Gamma=1.65_{-0.03}^{+0.03}$, flatter than $\Gamma\approx1.8$ typically found for AGN. While the spectral slope of BL and X-ray unabsorbed AGN is consistent with typical values ($\Gamma=1.79_{-0.01}^{+0.01}$), a significant flattening is seen in NL AGN and heavily-absorbed sources ($\Gamma=1.60_{-0.05}^{+0.08}$ and $\Gamma=1.38_{-0.12}^{+0.12}$, respectively), likely due to the effect of absorption and to the contribution from Compton reflection to the high-energy flux (E>10 keV). We find that the typical reflection fraction in our spectra is $R\approx0.5$ (for $\Gamma=1.8$), with a tentative indication of an increase of the reflection strength with column density. While there is no significant evidence for a dependence of the photon index with X-ray luminosity in our sample, we find that $R$ decreases with luminosity, with relatively high levels of reflection ($R\approx1.2$) for $L_{\rm 10-40 keV}<10^{44}$ erg/s and $R\approx0.3$ for $L_{\rm 10-40 keV}>10^{44}$ erg/s AGN, assuming $\Gamma=1.8$.

We identify sources with extremely hard X-ray spectra (i.e., with photon indices of Gamma<0.6 in the 13 sq. deg. NuSTAR serendipitous survey, to search for the most highly obscured AGNs detected at >10 keV. Eight extreme NuSTAR sources are identified, and we use the NuSTAR data in combination with lower energy X-ray observations (from Chandra, Swift XRT, and XMM-Newton) to characterize the broad-band (0.5-24 keV) X-ray spectra. We find that all of the extreme sources are highly obscured AGNs, including three robust Compton-thick (CT; N_H > 1.5e24 cm^-2) AGNs at low redshift (z<0.1), and a likely-CT AGN at higher redshift (z=0.16). Most of the extreme sources would not have been identified as highly obscured based on the low energy (<10 keV) X-ray coverage alone. The multiwavelength properties (e.g., optical spectra and X-ray/MIR luminosity ratios) provide further support for the eight sources being significantly obscured. Correcting for absorption, the intrinsic rest-frame 10-40 keV luminosities of the extreme sources cover a broad range, from ~ 5 x 10^42 to 10^45 erg s^-1. The estimated number counts of CT AGNs in the NuSTAR serendipitous survey are in broad agreement with model expectations based on previous X-ray surveys, except for the lowest redshifts (z<0.07) where we measure a high CT fraction of f_CT^obs = 30 (+16 -12) %. For the small sample of CT AGNs, we find a high fraction of galaxy major mergers (50 +/- 33%) compared to control samples of "normal" AGNs.

We present a sample of 10 low-mass active galactic nuclei (AGNs) selected from the 40-month NuSTAR serendipitous survey. The sample is selected to have robust NuSTAR detections at $3 - 24$~keV, to be at $z < 0.3$, and to have optical r-band magnitudes at least 0.5~mag fainter than an $L_\star$ galaxy at its redshift. The median values of absolute magnitude, stellar mass and 2--10 X-ray luminosity of our sample are $\langle M_r\rangle = -20.03$, $\langle M_\star\rangle = 4.6\times10^{9}M_\odot$, and $\langle L_{2-10\mathrm{keV}}\rangle = 3.1\times10^{42}$ erg s$^{-1}$, respectively. Five objects have detectable broad H$\alpha$ emission in their optical spectra, indicating black-hole masses of $(1.1-10.4)\times 10^6 M_\odot$. We find that $30^{+17}_{-10}\%$ of the galaxies in our sample do not show AGN-like optical narrow emission lines, and one of the ten galaxies in our sample, J115851+4243.2, shows evidence for heavy X-ray absorption. This result implies that a non-negligible fraction of low-mass galaxies might harbor accreting massive black holes that are missed by optical spectroscopic surveys and $<10$ keV X-ray surveys. The mid-IR colors of our sample also indicate these optically normal low-mass AGNs cannot be efficiently identified with typical AGN selection criteria based on WISE colors. While the hard ($>10$ keV) X-ray selected low-mass AGN sample size is still limited, our results show that sensitive NuSTAR observations are capable of probing faint hard X-ray emission originating from the nuclei of low-mass galaxies out to moderate redshift ($z<0.3$), thus providing a critical step in understanding AGN demographics in low-mass galaxies.

Tidal disruption events (TDEs), in which stars are gravitationally disrupted as they pass close to the supermassive black holes in the centres of galaxies, are potentially important probes of strong gravity and accretion physics. Most TDEs have been discovered in large-area monitoring surveys of many 1000s of galaxies, and the rate deduced for such events is relatively low: one event every 10$^4$ - 10$^5$ years per galaxy. However, given the selection effects inherent in such surveys, considerable uncertainties remain about the conditions that favour TDEs. Here we report the detection of unusually strong and broad helium emission lines following a luminous optical flare (Mv < -20.1 mag) in the nucleus of the nearby ultra-luminous infrared galaxy F01004-2237. The particular combination of variability and post-flare emission line spectrum observed in F01004-2237 is unlike any known supernova or active galactic nucleus. Therefore, the most plausible explanation for this phenomenon is a TDE -- the first detected in a galaxy with an ongoing massive starburst. The fact that this event has been detected in repeat spectroscopic observations of a sample of 15 ultra-luminous infrared galaxies over a period of just 10 years suggests that the rate of TDEs is much higher in such objects than in the general galaxy population.

We present the first full catalog and science results for the NuSTAR serendipitous survey. The catalog incorporates data taken during the first 40 months of NuSTAR operation, which provide ~20Ms of effective exposure time over 331 fields, with an areal coverage of 13 sq deg, and 497 sources detected in total over the 3-24 keV energy range. There are 276 sources with spectroscopic redshifts and classifications, largely resulting from our extensive campaign of ground-based spectroscopic followup. We characterize the overall sample in terms of the X-ray, optical, and infrared source properties. The sample is primarily comprised of active galactic nuclei (AGNs), detected over a large range in redshift from z = 0.002 - 3.4 (median of <z> = 0.56), but also includes 16 spectroscopically confirmed Galactic sources. There is a large range in X-ray flux, from log( f_3-24keV / erg s^-1 cm^-2 ) ~ -14 to -11, and in rest-frame 10-40 keV luminosity, from log( L_10-40keV / erg s^-1 ) ~ 39 to 46, with a median of 44.1. Approximately 79% of the NuSTAR sources have lower energy (<10 keV) X-ray counterparts from XMM-Newton, Chandra, and Swift/XRT. The mid-infrared (MIR) analysis, using WISE all-sky survey data, shows that MIR AGN color selections miss a large fraction of the NuSTAR-selected AGN population, from ~15% at the highest luminosities (Lx > 10^44 erg s^-1) to ~80% at the lowest luminosities (Lx < 10^43 erg s^-1). Our optical spectroscopic analysis finds that the observed fraction of optically obscured AGNs (i.e., the Type 2 fraction) is F_Type2 = 53(+14-15)%, for a well-defined subset of the 8-24 keV selected sample. This is higher, albeit at a low significance level, than the Type 2 fraction measured for redshift- and luminosity-matched AGNs selected by <10 keV X-ray missions.

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.

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.

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.

We investigate active galactic nuclei (AGN) candidates within the FourStar Galaxy Evolution Survey (ZFOURGE) to determine the impact they have on star-formation in their host galaxies. We first identify a population of radio, X-ray, and infrared-selected AGN by cross-matching the deep $K_{s}$-band imaging of ZFOURGE with overlapping multi-wavelength data. From this, we construct a mass-complete (log(M$_{*}$/M$_{\odot}$) $\ge$ 9.75), AGN luminosity limited sample of 235 AGN hosts over z = 0.2 - 3.2. We compare the rest-frame U - V versus V - J (UVJ) colours and specific star-formation rates (sSFRs) of the AGN hosts to a mass-matched control sample of inactive (non-AGN) galaxies. UVJ diagnostics reveal AGN tend to be hosted in a lower fraction of quiescent galaxies and a higher fraction of dusty galaxies than the control sample. Using 160{\mu}m Herschel PACS data, we find the mean specific star-formation rate of AGN hosts to be elevated by 0.34$\pm$0.07 dex with respect to the control sample across all redshifts. This offset is primarily driven by infrared-selected AGN, where the mean sSFR is found to be elevated by as much as a factor of ~5. The remaining population, comprised predominantly of X-ray AGN hosts, is found mostly consistent with inactive galaxies, exhibiting only a marginal elevation. We discuss scenarios that may explain these findings and postulate that AGN are less likely to be a dominant mechanism for moderating galaxy growth via quenching than has previously been suggested.

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.

We present the infrared (IR) and X-ray properties of a sample of 33 mid-IR luminous quasars ($\nu$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 $\Phi$=(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.

We present the first results from the KMOS AGN Survey at High redshift (KASHz), a VLT/KMOS integral-field spectroscopic survey of z>0.6 AGN. We present galaxy-integrated spectra of 89 X-ray AGN (Lx=10^42-10^45 erg/s), for which we observed [O III] (z=1.1-1.7) or Halpha emission (z=0.6-1.1). The targets have X-ray luminosities representative of the parent AGN population and we explore the emission-line luminosities as a function of X-ray luminosity. For the [O III] targets, ~50 per cent have ionised gas velocities indicative of gas that is dominated by outflows and/or highly turbulent material (i.e., overall line-widths >~600 km/s). The most luminous half (i.e., Lx>6x10^43 erg/s) have a >~2 times higher incidence of such velocities. On the basis of our results, we find no evidence that X-ray obscured AGN are more likely to host extreme kinematics than unobscured AGN. Our KASHz sample has a distribution of gas velocities that is consistent with a luminosity-matched sample of z<0.4 AGN. This implies little evolution in the prevalence of ionised outflows, for a fixed AGN luminosity, despite an order-of-magnitude decrease in average star-formation rates over this redshift range. Furthermore, we compare our Halpha targets to a redshift-matched sample of star-forming galaxies and despite a similar distribution of Halpha luminosities and likely star-formation rates, we find extreme ionised gas velocities are up to ~10x more prevalent in the AGN-host galaxies. Our results reveal a high prevalence of extreme ionised gas velocities in high-luminosity X-ray AGN and imply that the most powerful ionised outflows in high-redshift galaxies are driven by AGN activity.

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.

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.

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.

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.

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.

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.

We use deep panchromatic datasets in the GOODS-N field, from GALEX to the deepest Herschel far-infrared and VLA radio continuum imaging, to explore, using mass-complete samples, the evolution of the star formation activity and dust attenuation of star-forming galaxies to z~4. Our main results can be summarized as follows: i) the slope of the SFR-M correlation is consistent with being constant, and equal to ~0.8 at least up to z~1.5, while its normalization keeps increasing with redshift; ii) for the first time here we are able to explore the FIR-radio correlation for a mass-selected sample of star-forming galaxies: the correlation does not evolve up to z~4; iii) we confirm that galaxy stellar mass is a robust proxy for UV dust attenuation in star-forming galaxies, with more massive galaxies being more dust attenuated, strikingly we find that this attenuation relation evolves very weakly with redshift, the amount of dust attenuation increasing by less than 0.3 magnitudes over the redshift range [0.5-4] for a fixed stellar mass, as opposed to a tenfold increase of star formation rate; iv) the correlation between dust attenuation and the UV spectral slope evolves in redshift, with the median UV spectral slope of star-forming galaxies becoming bluer with redshift. By z~3, typical UV slopes are inconsistent, given the measured dust attenuation, with the predictions of commonly used empirical laws. Finally, building on existing results, we show that gas reddening is marginally larger (by a factor of around 1.3) than stellar reddening at all redshifts probed, and also that the amount of dust attenuation at a fixed ISM metallicity increases with redshift. We speculate that our results support evolving ISM conditions of typical star-forming galaxies such that at z~1.5 Main Sequence galaxies have ISM conditions getting closer to those of local starbursts.

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.

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.

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.

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.

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.

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.

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.

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.

We investigate the effect of active galactic nucleus (AGN) variability on the observed connection between star formation and black hole accretion in extragalactic surveys. Recent studies have reported relatively weak correlations between observed AGN luminosities and the properties of AGN hosts, which has been interpreted to imply that there is no direct connection between AGN activity and star formation. However, AGNs may be expected to vary significantly on a wide range of timescales (from hours to Myr) that are far shorter than the typical timescale for star formation (>~100 Myr). This variability can have important consequences for observed correlations. We present a simple model in which all star-forming galaxies host an AGN when averaged over ~100 Myr timescales, with long-term average AGN accretion rates that are perfectly correlated with the star formation rate (SFR). We show that reasonable prescriptions for AGN variability reproduce the observed weak correlations between SFR and L_AGN in typical AGN host galaxies, as well as the general trends in the observed AGN luminosity functions, merger fractions, and measurements of the average AGN luminosity as a function of SFR. These results imply there may be a tight connection between AGN activity and SFR over galaxy evolution timescales, and that the apparent similarities in rest-frame colors, merger rates, and clustering of AGNs compared to "inactive" galaxies may be due primarily to AGN variability. The results provide motivation for future deep, wide extragalactic surveys that can measure the distribution of AGN accretion rates as a function of SFR.

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)

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 \Gamma=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.

We present a measurement of the average supermassive black hole accretion rate (BHAR) as a function of star formation rate (SFR) for galaxies in the redshift range 0.25<z<0.8. We study a sample of 1,767 far-IR selected star-forming galaxies in the 9 deg^2 Bo\"otes multiwavelength survey field. The SFR is estimated using 250 micron observations from the Herschel Space Observatory, for which the contribution from the AGN is minimal. In this sample, 121 AGNs are directly identified using X-ray or mid-IR selection criteria. We combined these detected AGNs and an X-ray stacking analysis for undetected sources to study the average BHAR for all of the star-forming galaxies in our sample. We find an almost linear relation between the average BHAR (in M_sun/year) and the SFR (in M_sun/year) for galaxies across a wide SFR range 0.85<log SFR<2.56 : log BHAR=(-3.72\pm0.52)+(1.05\pm0.33) log SFR. This global correlation between SFR and average BHAR is consistent with a simple picture in which SFR and AGN activity are tightly linked over galaxy evolution timescales.

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.

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.

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.

We characterize the incidence of active galactic nuclei (AGNs) is 0.3 < z < 1 star-forming galaxies by applying multi-wavelength AGN diagnostics (X-ray, optical, mid-infrared, radio) to a sample of galaxies selected at 70-micron from the Far-Infrared Deep Extragalactic Legacy survey (FIDEL). Given the depth of FIDEL, we detect "normal" galaxies on the specific star formation rate (sSFR) sequence as well as starbursting systems with elevated sSFR. We find an overall high occurrence of AGN of 37+/-3%, more than twice as high as in previous studies of galaxies with comparable infrared luminosities and redshifts but in good agreement with the AGN fraction of nearby (0.05 < z < 0.1) galaxies of similar infrared luminosities. The more complete census of AGNs comes from using the recently developed Mass-Excitation (MEx) diagnostic diagram. This optical diagnostic is also sensitive to X-ray weak AGNs and X-ray absorbed AGNs, and reveals that absorbed active nuclei reside almost exclusively in infrared-luminous hosts. The fraction of galaxies hosting an AGN appears to be independent of sSFR and remains elevated both on the sSFR sequence and above. In contrast, the fraction of AGNs that are X-ray absorbed increases substantially with increasing sSFR, possibly due to an increased gas fraction and/or gas density in the host galaxies.

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.

We present the mid-infrared spectra of seven of the most powerful radio-galaxies known to exist at 1.5 < z < 2.6. The radio emission of these sources is dominated by the AGN with 500 MHz luminosities in the range 10^27.8 - 10^29.1 W/Hz. The AGN signature is clearly evident in the mid-infrared spectra, however, we also detect polycyclic aromatic hydrocarbons emission, indicative of prodigious star formation at a rate of up to ~1000 Msun/yr. Interestingly, we observe no significant correlation between AGN power and star formation in the host galaxy. We also find most of these radio galaxies to have weak 9.7 um silicate absorption features (tau_{9.7 um} < 0.8) which implies that their mid-IR obscuration is predominantly due to the dusty torus that surrounds the central engine, rather than the host galaxy. The tori are likely to have an inhomogeneous distribution with the obscuring structure consisting of individual clouds. We estimate that these radio galaxies have already formed the bulk of their stellar mass and appear to lie at a stage in their evolution where the obscured AGN dominates the energy output of the system but star formation is also prevalent.

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.

We present results from deep X-ray stacking of >4000 high redshift galaxies from z~1 to 8 using the 4 Ms Chandra Deep Field South (CDF-S) data, the deepest X-ray survey of the extragalactic sky to date. The galaxy samples were selected using the Lyman break technique based primarily on recent HST ACS and WFC3 observations. Based on such high specific star formation rates (sSFRs): log SFR/M* > -8.7, we expect that the observed properties of these LBGs are dominated by young stellar populations. The X-ray emission in LBGs, eliminating individually detected X-ray sources (potential AGN), is expected to be powered by X-ray binaries and hot gas. We find, for the first time, evidence of evolution in the X-ray/SFR relation. Based on X-ray stacking analyses for z<4 LBGs (covering ~90% of the Universe's history), we find that the 2-10 keV X-ray luminosity evolves weakly with redshift (z) and SFR as log LX = 0.93 log (1+z) + 0.65 log SFR + 39.80. By comparing our observations with sophisticated X-ray binary population synthesis models, we interpret that the redshift evolution of LX/SFR is driven by metallicity evolution in HMXBs, likely the dominant population in these high sSFR galaxies. We also compare these models with our observations of X-ray luminosity density (total 2-10 keV luminosity per Mpc^3) and find excellent agreement. While there are no significant stacked detections at z>5, we use our upper limits from 5<z<8 LBGs to constrain the SMBH accretion history of the Universe around the epoch of reionization.

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.

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.

We explore the effects of active galactic nuclei (AGN) and star formation activity on the infrared (0.3-1000 microns) spectral energy distributions of luminous infrared galaxies from z = 0.5 to 4.0. We have compiled a large sample of 151 galaxies selected at 24 microns (S24 > 100 uJy) in the GOODS-N and ECDFS fields for which we have deep Spitzer IRS spectroscopy, allowing us to decompose the mid-IR spectrum into contributions from star formation and AGN activity. A significant portion (~25%) of our sample is dominated by an AGN in the mid-IR. Based on the mid-IR classification, we divide our full sample into four sub-samples: z~1 star-forming (SF) sources; z~2 SF sources; AGN with clear 9.7 micron silicate absorption; and AGN with featureless mid-IR spectra. From our large spectroscopic sample and wealth of multi-wavelength data, including deep Herschel imaging at 100, 160, 250, 350, and 500 microns, we use 95 galaxies with complete spectral coverage to create a composite spectral energy distribution (SED) for each sub-sample. We then fit a two-temperature component modified blackbody to the SEDs. We find that the IR SEDs have similar cold dust temperatures, regardless of the mid-IR power source, but display a marked difference in the warmer dust temperatures. We calculate the average effective temperature of the dust in each sub-sample and find a significant (~20 K) difference between the SF and AGN systems. We compare our composite SEDs to local templates and find that local templates do not accurately reproduce the mid-IR features and dust temperatures of our high redshift systems. High redshift IR luminous galaxies contain significantly more cool dust than their local counterparts. We find that a full suite of photometry spanning the IR peak is necessary to accurately account for the dominant dust temperature components in high redshift IR luminous galaxies.

We reproduce the mid-infrared to radio galaxy counts with a new empirical model based on our current understanding of the evolution of main-sequence (MS) and starburst (SB) galaxies. We rely on a simple Spectral Energy Distribution (SED) library based on Herschel observations: a single SED for the MS and another one for SB, getting warmer with redshift. Our model is able to reproduce recent measurements of galaxy counts performed with Herschel, including counts per redshift slice. This agreement demonstrates the power of our 2 Star-Formation Modes (2SFM) decomposition for describing the statistical properties of infrared sources and their evolution with cosmic time. We discuss the relative contribution of MS and SB galaxies to the number counts at various wavelengths and flux densities. We also show that MS galaxies are responsible for a bump in the 1.4 GHz radio counts around 50 {\mu}Jy. Material of the model (predictions, SED library, mock catalogs...) is available online at http://irfu.cea.fr/Sap/Phocea/Page/index.php?id=537.

Using deep 100-160 micron observations in GOODS-S from the GOODS-H survey, combined with HST/WFC3 NIR imaging from CANDELS, we present the first morphological analysis of a complete, FIR selected sample of 52 ULIRGs at z~2. We also make use of a comparison sample of galaxies without Herschel detections but with the same z and magnitude distribution. Our visual classifications of these two samples indicate that the fraction of objects with disk and spheroid morphologies is roughly the same but that there are significantly more mergers, interactions, and irregular galaxies among the ULIRGs. The combination of disk and irregular/interacting morphologies suggests that early stage interactions and minor mergers could play an important role in ULIRGs at z~2. We compare these fractions with those of a z~1 sample across a wide luminosity range and find that the fraction of disks decreases systematically with L_IR while the fraction of mergers and interactions increases, as has been observed locally. At comparable luminosities, the fraction of ULIRGs with various morphological classifications is similar at z~2 and z~1. We investigate the position of the ULIRGs, along with 70 LIRGs, on the specific star formation rate versus redshift plane, and find 52 systems to be starbursts (lie more than a factor of 3 above the main sequence relation). The morphologies of starbursts are dominated by interacting and merging systems (50%). If irregular disks are included as potential minor mergers, then we find that up to 73% of starbursts are involved in a merger or interaction at some level. Although the final coalescence of a major merger may not be required for the high luminosities of ULIRGs at z~2 as is the case locally, the large fraction of interactions at all stages and potential minor mergers suggest that the high star formation rates of ULIRGs are still largely externally triggered at z~2.

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)

We provide evidence for a correlation between the presence of giant clumps and the occurrence of active galactic nuclei (AGN) in disk galaxies. Giant clumps of 10^8-9 Msun arise from violent gravitational instability in gas-rich galaxies, and it has been proposed that this instability could feed supermassive black holes (BH). We use emission line diagnostics to compare a sample of 14 clumpy (unstable) disks and a sample of 13 smoother (stable) disks at redshift z~0.7. The majority of clumpy disks in our sample have a high probability of containing AGN. Their [OIII] emission line is strongly excited, inconsistent with low-metallicity star formation alone. [NeIII] excitation is also higher. Stable disks rarely have such properties. Stacking ultra sensitive Chandra observations (4 Ms) reveals an X-ray excess in clumpy galaxies, which cannot be solely due to star formation and confirms the presence of AGN. The clumpy galaxies in our intermediate-redshift sample have properties typical of gas-rich disk galaxies rather than mergers, being in particular on the Main Sequence of star formation. This suggests that our findings apply to the physically-similar and numerous gas-rich unstable disks at z>1. Using the observed [OIII] and X-ray luminosities, we conservatively estimate that AGN hosted by clumpy disks have typical bolometric luminosities of the order of a few 10^43 erg/s, BH growth rates ~10^-2 Msun/yr, and that these AGN are substantially obscured in X-rays. This moderate-luminosity mode could be sufficient to provide a large fraction of today's BH mass over a couple of Gyr given that our observations suggest a high duty cycle (>10%), accretion bursts with higher luminosities being possible over shorter phases. The observed evolution of disk instabilities with mass and redshift could explain the simultaneous downsizing of star formation and of BH growth.

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.

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.

We present the deepest far-IR observations obtained with Herschel and examine the 3-500um SEDs of galaxies at 0<z<2.5, supplemented by a local reference sample from IRAS, ISO, Spitzer and AKARI data. We find that the ratio of total IR luminosity to rest-frame 8um luminosity, IR8 (=Lir/L8), follows a Gaussian distribution centered on IR8=4 and defines an IR main sequence (MS). A minority population (<20 %) of outliers producing a tail skewed toward higher values of IR8 consist of starbursts (SB) with compact projected star formation densities. IR8 can be used to separate galaxies with normal and extended modes of star formation from compact SBs with high-IR8, high projected IR surface brightness (>3x10^10 Lsun kpc^-2) and a high specific SFR (i.e., SBs). The rest-frame, UV-2700A size of these distant SBs is typically half that of MS galaxies, supporting the correlation between star formation density and SB activity that is measured for the local sample. Locally, (U)LIRGs are systematically in the SB mode, whereas most distant (U)LIRGs form stars in the "normal" MS mode. This confusion between two modes of star formation is the cause of the so-called "mid-IR excess" population of galaxies found at z>1.5 by previous studies. MS galaxies have strong PAH emission line features, a broad far-IR bump resulting from a combination of dust temperatures (Tdust~15-50 K), and an effective Tdust~31 K, as derived from the peak wavelength of their IR SED. Galaxies in the SB regime instead exhibit weak PAH EW and a sharper far-IR bump with an effective Tdust~40 K. Finally, we present evidence that the mid-to-far IR emission of X-ray AGNs is predominantly produced by star formation and that candidate dusty AGNs with a power-law emission in the mid-IR systematically occur in compact, dusty SBs. After correcting for the effect of SBs on IR8, we identify new candidates for extremely obscured AGNs.

We have carried out a survey of long 50ks XMM-Newton observations of a sample of bright, variable AGN. We found a distinctive energy dependence of the variability in RXJ0136.9-3510 where the fractional variability increases from 0.3 to 2 keV, and then remains constant. This is in sharp contrast to other AGN where the X-ray variability is either flat or falling with energy, sometimes with a peak at $\sim$~2 keV superimposed on the overall trend. Intriguingly these unusual characteristics of the variability are shared by one other AGN, namely RE J1034+396, which is so far unique showing a significant X-ray QPO. In addition the broad band spectrum of RXJ0136.9-3510 is also remarkably similar to that of RE J1034+396, being dominated by a huge soft excess in the EUV-soft X-ray bandpass. The bolometric luminosity of RX J0136.9-3510 gives an Eddington ratio of about 2.7 for a black hole mass (from the H beta line width) of $7.9 \times 10^{7}M_{\sun}$. This mass is about a factor of 50 higher than that of RE J1034+396, making any QPO undetectable in this length of observation. Nonetheless, its X-ray spectral and variability similarities suggest that RE J1034+396 is simply the closest representative of a new class of AGN spectra, representing the most extreme mass accretion rates.

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.

We explore the gas-to-dust mass ratio (G/D) and the CO luminosity-to-Mgas conversion factor (a_co) of two well studied galaxies in the GOODS-N field, that are expected to have different star forming modes, the starburst GN20 at z=4.05 and the normal star-forming galaxy BzK-21000 at z=1.52. Detailed sampling is available for their Rayleigh-Jeans emission via ground based mm interferometry (1.1-6.6mm) along with Herschel, PACS and SPIRE data that probe the peak of their infrared emission. Using the physically motivated Draine & Li (2007) models, as well as a modified black body function, we measure the dust mass (Md) of the sources and find 2.0^{+0.7}_{-0.6} x 10^{9} Msun for GN20 and 8.6^{+0.6}_{-0.9} x 10^{8} Msun for BzK-21000. The addition of mm data reduces the uncertainties of the derived Md by a factor of ~2, allowing the use of the local G\D vs metallicity relation to place constraints on the a_co values of the two sources. For GN20 we derive a conversion factor of a_co < 1.0 Msun pc^{-2}(K km s^{-1})^{-1}, consistent with that of local ULIRGs, while for BzK-21000 we find a considerably higher value, a_co ~4.0 Msun pc^{-2}(K km s^{-1})^{-1}, in agreement with an independent kinematic derivation reported previously. The implied star formation efficiency is ~25 Lsun/Msun for BzK-21000, a factor of ~5-10 lower than that of GN20. The findings for these two sources support the existence of different disk-like and starburst star-formation modes in distant galaxies, although a larger sample is required to draw statistically robust results

Using extremely deep PACS 100- and 160um Herschel data from the GOODS-Herschel program, we identify 21 infrared bright galaxies previously missed in the deepest 24um surveys performed by MIPS. These MIPS dropouts are predominantly found in two redshift bins, centred at z ~0.4 and ~1.3. Their S_100/S_24 flux density ratios are similar to those of local LIRGs and ULIRGs, whose silicate absorption features at 18um (at z ~ 0.4) and 9.7um (at z ~ 1.3) are shifted into the 24um MIPS band at these redshifts. The high-z sub-sample consists of 11 infrared luminous sources, accounting for ~2% of the whole GOODS-Herschel sample and putting strong upper limits on the fraction of LIRGs/ULIRGs at 1.0<z<1.7 that are missed by the 24um surveys. We find that a S_100/S_24 > 43 colour cut selects galaxies with a redshift distribution similar to that of the MIPS dropouts and when combined with a second colour cut, S_16/S_8 > 4, isolates sources at 1.0 < z < 1.7. We show that these sources have elevated specific star formation rates (sSFR) compared to main sequence galaxies at these redshifts and are likely to be compact starbursts with moderate/strong 9.7um silicate absorption features in their mid-IR spectra. Herschel data reveal that their infrared luminosities extrapolated from the 24um flux density are underestimated, on average, by a factor of ~3. These silicate break galaxies account for 16% (8%) of the ULIRG (LIRG) population in the GOODS fields, indicating a lower limit in their space density of 2.0 \times 10^(-5) Mpc^(-3). Finally, we provide estimates of the fraction of z < 2 MIPS dropout sources as a function of the 24-, 100-, 160-, 250- and 350um sensitivity limits, and conclude that previous predictions of a population of silicate break galaxies missed by the major 24um extragalactic surveys have been overestimated.

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.

We present a galaxy (SAGE1CJ053634.78-722658.5) at a redshift of 0.14 of which the IR is entirely dominated by emission associated with the AGN. We present the 5-37 um Spitzer/IRS spectrum and broad wavelength SED of SAGE1CJ053634, an IR point-source detected by Spitzer/SAGE (Meixner et al 2006). The source was observed in the SAGE-Spec program (Kemper et al., 2010) and was included to determine the nature of sources with deviant IR colours. The spectrum shows a redshifted (z=0.14+-0.005) silicate emission feature with an exceptionally high feature-to-continuum ratio and weak polycyclic aromatic hydrocarbon (PAH) bands. We compare the source with models of emission from dusty tori around AGNs from Nenkova et al. (2008). We present a diagnostic diagram that will help to identify similar sources based on Spitzer/MIPS and Herschel/PACS photometry. The SED of SAGE1CJ053634 is peculiar because it lacks far-IR emission and a clear stellar counterpart. We find that the SED and the IR spectrum can be understood as emission originating from the inner ~10 pc around an accreting black hole. There is no need to invoke emission from the host galaxy, either from the stars or from the interstellar medium, although a possible early-type host galaxy cannot be excluded based on the SED analysis. The hot dust around the accretion disk gives rise to a continuum, which peaks at 4 um, whereas the strong silicate features may arise from optically thin emission of dusty clouds within ~10 pc around the black hole. The weak PAH emission does not appear to be linked to star formation, as star formation templates strongly over-predict the measured far-IR flux levels. The SED of SAGE1CJ053634 is rare in the local universe but may be more common in the more distant universe. The conspicuous absence of host-galaxy IR emission places limits on the far-IR emission arising from the dusty torus alone.

We present the first measurement of the spatial clustering of mid-infrared selected obscured and unobscured quasars, using a sample in the redshift range 0.7 < z < 1.8 selected from the 9 deg^2 Bootes multiwavelength survey. Recently the Spitzer Space Telescope and X-ray observations have revealed large populations of obscured quasars that have been inferred from models of the X-ray background and supermassive black hole evolution. To date, little is known about obscured quasar clustering, which allows us to measure the masses of their host dark matter halos and explore their role in the cosmic evolution of black holes and galaxies. In this study we use a sample of 806 mid-infrared selected quasars and ~250,000 galaxies to calculate the projected quasar-galaxy cross-correlation function w_p(R). The observed clustering yields characteristic dark matter halo masses of log (M_halo [h^-1 M_sun]) = 12.7^+0.4_-0.6 and 13.3^+0.3_-0.4 for unobscured quasars (QSO-1s) and obscured quasars (Obs-QSOs), respectively. The results for QSO-1s are in excellent agreement with previous measurements for optically-selected quasars, while we conclude that the Obs-QSOs are at least as strongly clustered as the QSO-1s. We test for the effects of photometric redshift errors on the optically-faint Obs-QSOs, and find that our method yields a robust lower limit on the clustering; photo-z errors may cause us to underestimate the clustering amplitude of the Obs-QSOs by at most ~20%. We compare our results to previous studies, and speculate on physical implications of stronger clustering for obscured quasars.

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

Using a suite of X-ray, mid-IR and optical active galactic nuclei (AGN) luminosity indicators, we search for Compton-thick (CT) AGNs with intrinsic L_X>10^42erg/s at z~0.03-0.2, a region of parameter space which is currently poorly constrained by deep narrow-field and high-energy (E>10keV) all-sky X-ray surveys. We have used the widest XMM-Newton survey (the serendipitous source catalogue) to select a representative sub-sample (14; ~10%) of the 147 X-ray undetected candidate CT AGNs in the Sloan Digital Sky Survey (SDSS) with f_X/f_[OIII]<1; the 147 sources account for ~50% of the overall Type-2 AGN population in the SDSS-XMM overlap region. We use mid-IR spectral decomposition analyses and emission-line diagnostics, determined from pointed Spitzer-IRS spectroscopic observations of these candidate CT AGNs, to estimate the intrinsic AGN emission (predicted L_X,2-10keV (0.2-30)x10^42erg/s). On the basis of the optical [OIII], mid-IR [OIV] and 6um AGN continuum luminosities we conservatively find that the X-ray emission in at least 6/14 (>43%) of our sample appear to be obscured by CT material with N_H>1.5x10^24cm^-2. Under the reasonable assumption that our 14 AGNs are representative of the overall X-ray undetected AGN population in the SDSS-XMM parent sample, we find that >20% of the optical Type-2 AGN population are likely to be obscured by CT material. This implies a space-density of log(Phi) >-4.9Mpc^-3 for CT AGNs with L_X>10^42erg/s at z~0.1, which we suggest may be consistent with that predicted by X-ray background synthesis models. Furthermore, using the 6um continuum luminosity to infer the intrinsic AGN luminosity and the stellar velocity dispersion to estimate M_BH, we find that the most conservatively identified CT AGNs in this sample may harbour some of the most rapidly growing black holes (median M_BH~3x10^7M_o) in the nearby Universe, with a median Eddington ratio of ~0.2.

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.

We investigate the local supermassive black hole (SMBH) density function and relative mass accretion rates of all active galactic nuclei (AGNs) identified in a volume-limited sample of infrared (IR) bright galaxies (L_IR > 3 x 10^9 L_sun) to D<15 Mpc (Goulding & Alexander 2009). A database of accurate SMBH mass (M_BH) estimates is compiled from literature sources using physically motivated AGN modeling techniques (reverberation mapping, maser mapping and gas kinematics) and well-established indirect M_BH estimation methods (the M-sigma and M_BH-L_(K,bul) relations). For the three sources without previously published M_BH estimates, we use 2MASS K-band imaging and GALFIT to constrain the bulge luminosities, and hence SMBH masses. In general, we find the AGNs in the sample host SMBHs which are spread over a wide mass range (M_BH ~ (0.1-30) x 10^7 M_sun), but with the majority in the poorly studied M_BH ~ 10^6-10^7 M_sun region. Using sensitive hard X-ray (2-10 keV) and mid-IR constraints we calculate the bolometric luminosities of the AGNs (L_(Bol,AGN)) and use them to estimate relative mass accretion rates. We use these data to calculate the volume-average SMBH growth rate of galaxies in the local Universe and find that the AGNs hosting SMBHs in the mass range M_BH ~ 10^6-10^7 M_sun are dominated by optically unidentified AGNs. These relatively small SMBHs are acquiring a significant proportion of their mass in the present-day, and are amongst the most rapidly growing in the local Universe (SMBH mass doubling times of ~6 Gyrs). Additionally, we find tentative evidence for an increasing volume-weighted AGN fraction with decreasing SMBH mass in the M_BH ~ 10^6-10^8 M_sun range. Overall, we conclude that significant mass accretion onto small SMBHs may be missed in even the most sensitive optical surveys due to absent or weak optical AGN signatures.

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

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

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.

We have defined a sample of 63 AGN with strong forbidden high-ionisation line (FHIL) emission. These lines, with ionisation potentials >~ 100eV, respond to a portion of the spectrum that is often difficult to observe directly, thereby providing constraints on the EUV-soft X-ray continuum. The sources are selected from the Sloan Digital Sky Survey (SDSS) on the basis of their [Fe X]6374A emission, yielding one of the largest and the most homogeneous sample of FHIL-emitting galaxies. We fit a sequence of models to both FHILs ([Fe XI], [Fe X] and [Fe VII]) and lower-ionisation emission lines ([O III], [O I], H-alpha, [N II], [S II]) in the SDSS spectra. These data are combined with X-ray measurements from Rosat, which are available for half of the sample. The correlations between these parameters are discussed for both the overall sample and subsets defined by spectroscopic classifications. The primary results are evidence that: (1) the [Fe X] and [Fe XI] lines are photoionised and their strength is proportional to the continuum flux around 250 eV; (2) the FHIL-emitting clouds form a stratified outflow in which the [Fe X] and [Fe XI] source regions extend sufficiently close to the BLR that they are partially obscured in Seyfert 2s whereas the [Fe VII] source region is more extended and is unaffected by obscuration; (3) narrow-lined Seyfert 1s (NLS1s) tend to have the strongest [Fe X] flux (relative to lower-ionisation lines); and (4) the most extreme [Fe X] ratios (such as [Fe X]/[O III] or [Fe X]/[Fe VII]) are found in the NLS1s with the narrowest broad lines and appear to be an optical-band indication of objects with strong X-ray soft excesses.

We use the photoionisation code Cloudy to determine both the location and the kinematics of the optical forbidden, high ionisation line (hereafter, FHIL) emitting gas in the narrow line Seyfert 1 galaxy Ark 564. The results of our models are compared with the observed properties of these emission lines to produce a physical model that is used to explain both the kinematics and the source of this gas. The main features of this model are that the FHIL emitting gas is launched from the putative dusty torus and is quickly accelerated to its terminal velocity of a few hundred km/s. Iron-carrying grains are destroyed during this initial acceleration. This velocity is maintained by a balance between radiative forces and gravity in this super-Eddington source. Eventually the outflow is slowed at large radii by the gravitational forces of and interactions with the host galaxy. In this model, FHIL emission traces the transition between the AGN and bulge zones of influence.

Based on a new spectroscopic sample observed using the WHT, we examine the kinematic properties of the various emission line regions in narrow line Seyfert 1 galaxies (NLS1s) by modelling their profiles using multiple component fits. We interpret these results by comparison with velocity components observed for different lines species covered in the same spectrum, and equivalent components measured in the spectra of some broad line Seyfert 1s and a representative Seyfert 2 galaxy. We find that the fits to the Halpha and Hbeta line profiles in NLS1s require an additional broad (~3000km/s) component that might correspond to a suppressed broad line region with similar kinematics to those of typical broad line Seyfert 1s. From the profiles of the forbidden high ionisation lines (FHILs) in NLS1s, we find evidence that they appear to trace an `intermediate' velocity region with kinematics between the standard broad and narrow line regions. Weaker evidence of this region is also present in the profiles of the permitted Balmer lines. Finally, we note that despite having similar ionisation potentials, the relative intensities of the highly ionised lines of [Fe X]6374 and [FeXI]7892 show considerable dispersion from one galaxy to another. The interpretation of this requires further modelling, but suggests the possibility of using the ratio as a diagnostic to constrain the physical conditions of the FHIL emitting region and possibly the shape of the spectral energy distribution in the vicinity of 200eV. This spectral region is very difficult to observe directly due to photoelectric absorption both in our Galaxy and intrinsic to the source.