In order to answer some of the major open questions in the fields of supermassive black hole (SMBH) and galaxy evolution, a complete census of SMBH growth, i.e., active galactic nuclei (AGN), is required. Thanks to deep all-sky surveys, such as those by the Wide-field Infrared Survey Explorer (WISE) and the Spectrum-Roentgen-Gamma (SRG) missions, this task is now becoming feasible in the nearby Universe. Here, we present a new survey, the Local AGN Survey (LASr), with the goal of identifying AGN unbiased against obscuration and determining the intrinsic Compton-thick (CT) fraction. First, we construct the most complete all-sky sample of galaxies within 100 Mpc from astronomical databases (90% completeness for log M_star/M_sun ~9.4), four times deeper than the current local galaxy reference, the Two Micron All-Sky Survey Redshift Survey (2MRS), which turns out to miss 20% of known luminous AGN. These 49k galaxies serve as parent sample for LASr, called LASr-GPS. It contains 4.3k already known AGN, >=82% of these are estimated to have L(12um) < 10^42.3 erg/s, i.e., are low-luminosity AGN. As a first method for identifying Seyfert-like AGN, we use WISE-based infrared colours, finding 221 galaxies at L(12um) >= 10^42.3 erg/s to host an AGN at 90% reliability, This includes 61 new AGN candidates and implies and optical type 2 fraction of 50 to 71%. We quantify the efficiency of this technique and estimate the total number of AGN with L_int(2-10 keV) >= 10^42 erg/s in the volume to be 362^+145_-116 (8.6^+3.5_-2.8 10^-5 Mpc^-3). X-ray brightness estimates indicate the CT fraction to be 40-55% to explain the Swift non-detections of the infrared selected objects. One third of the AGN within 100 Mpc remain to be identified and we discuss the prospects for the eROSITA all-sky survey.

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.

In this work we present a robust quantification of X-ray selected AGN in local ($z \leq 0.25$) dwarf galaxies ($M_\mathrm{*} \leq 3 \times 10^9 \mathrm{M_\odot}$). We define a parent sample of 4,331 dwarf galaxies found within the footprint of both the MPA-JHU galaxy catalogue (based on SDSS DR8) and 3XMM DR7, performed a careful review of the data to remove misidentifications and produced a sample of 61 dwarf galaxies that exhibit nuclear X-ray activity indicative of an AGN. We examine the optical emission line ratios of our X-ray selected sample and find that optical AGN diagnostics fail to identify 85% of the sources. We then calculated the growth rates of the black holes powering our AGN in terms of their specific accretion rates ($\propto L_\mathrm{X}/M_\mathrm{*}$, an approximate tracer of the Eddington ratio). Within our observed sample, we found a wide range of specific accretion rates. After correcting the observed sample for the varying sensitivity of 3XMM, we found further evidence for a wide range of X-ray luminosities and specific accretion rates, described by a power law. Using this corrected AGN sample we also define an AGN fraction describing their relative incidence within the parent sample. We found the AGN fraction increases with host galaxy mass (up to $\approx$ 6%) for galaxies with X-ray luminosities between $10^{39} \mathrm{erg/s}$ and $10^{42} \mathrm{erg/s}$, and by extrapolating the power law to higher luminosities, we found evidence to suggest the fraction of luminous AGN ($L_\mathrm{X} \geq 10^{42.4} \mathrm{erg/s}$) is constant out to $z \approx 0.7$.

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

Using data from the MOSFIRE Deep Evolution Field (MOSDEF) survey, we present a census of AGN-driven ionized outflows in a sample of 159 AGNs at $1.4 \le z \le 3.8$. The sample spans AGN bolometric luminosities of $10^{44-47} \mathrm{~erg~s}^{-1}$ and includes both quiescent and star-forming galaxies extending across three orders of magnitude in stellar mass. We identify and characterize outflows from the \hbeta, [OIII], \halpha ~and [NII] emission line spectra. We detect outflows in $17\%$ of the AGNs, seven times more often than in a mass-matched sample of inactive galaxies in MOSDEF. The outflows are fast and galaxy-wide, with velocities of $\sim 400-3500 ~\mathrm{km~s}^{-1}$ and spatial extents of $0.3-11.0$ kpc. The incidence of outflows among AGNs is independent of the stellar mass of the host galaxy, with outflows detected in both star-forming and quiescent galaxies. This suggests that outflows exist across different phases in galaxy evolution. We investigate relations between outflow kinematic, spatial, and energetic properties and both AGN and host galaxy properties. Our results show that AGN-driven outflows are widespread in galaxies along the star-forming main sequence. The mass-loading factors of the outflows are typically $0.1-1$ and increase with AGN luminosity, capable of exceeding unity at $L_\mathrm{AGN} \gtrsim 10^{46.3} \mathrm{~erg~s}^{-1}$. In these more luminous sources the ionized outflow alone is likely sufficient to regulate star formation, and when combined with outflowing neutral and molecular gas may be able to quench star formation in their host galaxies.

Population synthesis models predict that high-mass X-ray binary (HMXB) populations produced in low metallicity environments should be more X-ray luminous, a trend supported by studies of nearby galaxies. This trend may be responsible for the observed increase of the X-ray luminosity ($L_{\mathrm{X}}$) per star formation rate (SFR) with redshift due to the decrease of metallicity ($Z$) at fixed stellar mass as a function of redshift. To test this hypothesis, we use a sample of 79 $z\sim2$ star-forming galaxies with oxygen abundance measurements from the MOSDEF survey, which obtained rest-frame optical spectra for $\sim1500$ galaxies in the CANDELS fields at $1.37<z<3.80$. Using Chandra data from the AEGIS-X Deep, Deep Field North, and Deep Field South surveys, we stack the X-ray data at the galaxy locations in bins of redshift and $Z$ because the galaxies are too faint to be individually detected. In agreement with previous studies, the average $L_{\mathrm{X}}$/SFR of our $z\sim2$ galaxy sample is enhanced by $\approx0.4-0.8$ dex relative to local HMXB $L_{\mathrm{X}}$-SFR scaling relations. Splitting our sample by $Z$, we find that $L_{\mathrm{X}}$/SFR and $Z$ are anti-correlated with 97% confidence. This observed $Z$ dependence for HMXB-dominated galaxies is consistent both with the local $L_{\mathrm{X}}$-SFR-$Z$ relation and a subset of population synthesis models. Although the statistical significance of the observed trends is weak due to the low X-ray statistics, these results constitute the first direct evidence connecting the redshift evolution of $L_{\mathrm{X}}$/SFR and the $Z$ dependence of HMXBs.

We present catalogs of calibrated photometry and spectroscopic redshifts in the Extended Groth Strip, intended for studies of photometric redshifts (photo-z's). The data includes ugriz photometry from CFHTLS and Y-band photometry from the Subaru Suprime camera, as well as spectroscopic redshifts from the DEEP2, DEEP3 and 3D-HST surveys. These catalogs incorporate corrections to produce effectively matched-aperture photometry across all bands, based upon object size information available in the catalog and Moffat profile point spread function fits. We test this catalog with a simple machine learning-based photometric redshift algorithm based upon Random Forest regression, and find that the corrected aperture photometry leads to significant improvement in photo-z accuracy compared to the original SExtractor catalogs from CFHTLS and Subaru. The deep ugrizY photometry and spectroscopic redshifts are well-suited for empirical tests of photometric redshift algorithms for LSST. The resulting catalogs are publicly available. We include a basic summary of the strategy of the DEEP3 Galaxy Redshift Survey to accompany the recent public release of DEEP3 data.

We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolution spectrographs ($R = \lambda/\Delta\lambda \sim 6500$), and 812 fibres transferring light to the high-resolution spectrograph ($R \sim 20\,000$). After a description of the instrument and its expected performance, a short overview is given of its operational scheme and planned 4MOST Consortium science; these aspects are covered in more detail in other articles in this edition of The Messenger. Finally, the processes, schedules, and policies concerning the selection of ESO Community Surveys are presented, commencing with a singular opportunity to submit Letters of Intent for Public Surveys during the first five years of 4MOST operations.

The discoveries made over the past 20 years by Chandra and XMM-Newton surveys in conjunction with multiwavelength imaging and spectroscopic data available in the same fields have significantly changed the view of the supermassive black hole (SMBH) and galaxy connection. These discoveries have opened up several exciting questions that are beyond the capabilities of current X-ray telescopes and will need to be addressed by observatories in the next two decades. As new observatories peer into the early Universe, we will begin to understand the physics and demographics of SMBH infancy (at $z>6$) and investigate the influence of their accretion on the formation of the first galaxies ($\S$ 2.1). We will also be able to understand the accretion and evolution over the cosmic history (at $z\sim$1-6) of the full population of black holes in galaxies, including low accretion rate, heavily obscured AGNs at luminosities beyond the reach of current X-ray surveys ($\S$2.2 and $\S$2.3), enabling us to resolve the connection between SMBH growth and their environment.

Much of the baryonic matter in the Universe, including the most active and luminous sources, are best studied in the X-ray band. Key advances in X-ray optics and detectors have paved the way for the Advanced X-ray Imaging Satellite (AXIS), a Probe-class mission that is a major improvement over Chandra, which has generated a steady stream of important discoveries for the past 2 decades. AXIS can be launched in the late 2020s and will transform our understanding in several major areas of astrophysics, including the growth and fueling of supermassive black holes, galaxy formation and evolution, the microphysics of cosmic plasmas, the time-variable universe, and a wide variety of cutting-edge studies. Relative to Chandra, the AXIS PSF is nearly twice as sharp on-axis; its field of view for subarcsecond imaging 70 times larger by area; its effective area at 1 keV is 10 times larger. The low-Earth orbit ensures a low and stable detector background, resulting in 50 times greater sensitivity than Chandra for extended sources. AXIS has a rapid repointing response with operations similar to Swift, but is 100 times more sensitive for time-domain science. These capabilities open up a vast discovery space and complement the next generation of astronomical observatories. A high-spectral-resolution mission (Athena) operating at the same time as a high-angular-resolution mission (AXIS) greatly increases the range of scientific discovery. AXIS will use lightweight X-ray optics made of thin single-crystal silicon mirrors developed at NASA Goddard. The detector array builds on a long legacy of X-ray CCD and provides improved photon localization, much faster readout time, and broader energy band. The estimated mission costs are consistent with the $1B Probe mission cost guideline.

We map the co-eval growth of galaxies and their central supermassive black holes in detail by measuring the incidence of Active Galactic Nuclei (AGN) in galaxies as a function of star formation rate (SFR) and redshift (to z~4). We combine large galaxy samples with deep Chandra X-ray imaging to measure the probability distribution of specific black hole accretion rates (LX relative to stellar mass) and derive robust AGN fractions and average specific accretion rates. First, we consider galaxies along the main sequence of star formation. We find a linear correlation between the average SFR and both the AGN fraction and average specific accretion rate across a wide range in stellar mass ($M_* \sim 10^{8.5-11.5}M_\odot$) and to at least z~2.5, indicating that AGN in main-sequence galaxies are driven by the stochastic accretion of cold gas. We also consider quiescent galaxies and find significantly higher AGN fractions than predicted, given their low SFRs, indicating that AGN in quiescent galaxies are fuelled by additional mechanisms (e.g. stellar winds). Next, we bin galaxies according to their SFRs relative to the main sequence. We find that the AGN fraction is significantly elevated for galaxies that are still star-forming but with SFRs below the main sequence, indicating further triggering mechanisms enhance AGN activity within these sub-main-sequence galaxies. We also find that the incidence of high-accretion-rate AGN is enhanced in starburst galaxies and evolves more mildly with redshift than within the rest of the galaxy population, suggesting mergers play a role in driving AGN activity in such high-SFR galaxies.

A semi-empirical model is presented that describes the distribution of Active Galactic Nuclei (AGN) on the cosmic web. It populates dark-matter halos in N-body simulations (MultiDark) with galaxy stellar masses using empirical relations based on abundance matching techniques, and then paints accretion events on these galaxies using state-of-the-art measurements of the AGN occupation of galaxies. The explicit assumption is that the large-scale distribution of AGN is independent of the physics of black-hole fueling. The model is shown to be consistent with current measurements of the two-point correlation function of AGN samples. It is then used to make inferences on the halo occupation of the AGN population. Mock AGN are found in halos with a broad distribution of masses with a mode of $\approx 10^{12}\,h^{-1} \, M_{\odot}$ and a tail extending to cluster-size halos. The clustering properties of the model AGN depend only weakly on accretion luminosity and redshift. The fraction of satellite AGN in the model increases steeply toward more massive halos, in contrast with some recent observational results. This discrepancy, if confirmed, could point to a dependence of the halo occupation of AGN on the physics of black-hole fueling.

The NuSTAR serendipitous survey has already uncovered a large number of Active Galactic Nuclei (AGN), providing new information about the composition of the Cosmic X-ray Background. For the AGN off the Galactic plane, it has been possible to use the existing X-ray archival data to improve source localizations, identify optical counterparts, and classify the AGN with optical spectroscopy. However, near the Galactic Plane, better X-ray positions are necessary to achieve optical or near-IR identifications due to the higher levels of source crowding. Thus, we have used observations with the Chandra X-ray Observatory to obtain the best possible X-ray positions. With eight observations, we have obtained coverage for 19 NuSTAR serendips within 12 deg of the plane. One or two Chandra sources are detected within the error circle of 15 of the serendips, and we report on these sources and search for optical counterparts. For one source (NuSTAR J202421+3350.9), we obtained a new optical spectrum and detected the presence of hydrogen emission lines. The source is Galactic, and we argue that it is likely a Cataclysmic Variable. For the other sources, the Chandra positions will enable future classifications in order to place limits on faint Galactic populations, including high-mass X-ray binaries and magnetars.

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 an analysis using the MOSFIRE Deep Evolution Field (MOSDEF) survey on the nature of "MIR-excess" galaxies, which have star formation rates (SFR) inferred from mid-infrared (MIR) data that is substantially elevated relative to that estimated from dust-corrected UV data. We use a sample of $\sim$200 galaxies and AGN at $1.40<z<2.61$ with 24 $\mu$m detections (rest-frame 8$\mu$m) from MIPS/\textit{Spitzer}. We find that the identification of MIR-excess galaxies strongly depends on the methodologies used to estimate IR luminosity ($\rm L_{IR}$) and to correct the UV light for dust attenuation. We find that extrapolations of the SFR from the observed 24 $\mu$m flux, using luminosity-dependent templates based on local galaxies, substantially overestimate $\rm L_{IR}$ in $z\sim2$ galaxies. By including \textit{Herschel} observations and using a stellar mass-dependent, luminosity-independent $\rm L_{IR}$, we obtain more reliable estimates of the SFR and a lower fraction of MIR-excess galaxies. Once stellar mass selection biases are taken into account, we identify $\sim24\%$ of our galaxies as MIR-excess. However, $\rm SFR_{H\alpha}$ is not elevated in MIR-excess galaxies compared to MIR-normal galaxies, indicating that the intrinsic fraction of MIR-excess may be lower. Using X-ray, IR, and optically-selected AGN in MOSDEF, we do not find a higher prevalence for AGN in MIR-excess galaxies relative to MIR-normal galaxies. A stacking analysis of X-ray undetected galaxies does not reveal a harder spectrum in MIR-excess galaxies relative to MIR-normal galaxies. Our analysis indicates that AGN activity does not contribute substantially to the MIR excess and instead implies that it is likely due to the enhanced PAH emission.

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 use deep Chandra X-ray imaging to measure the distribution of specific black hole accretion rates ($L_X$ relative to the stellar mass of the galaxy) and thus trace AGN activity within star-forming and quiescent galaxies, as a function of stellar mass (from $10^{8.5}-10^{11.5} M_\odot$) and redshift (to $z \sim 4$). We adopt near-infrared selected samples of galaxies from the CANDELS and UltraVISTA surveys, extract X-ray data for every galaxy, and use a flexible Bayesian method to combine these data and to measure the probability distribution function of specific black hole accretion rates, $\lambda_{sBHAR}$. We identify a broad distribution of $\lambda_{sBHAR}$ in both star-forming and quiescent galaxies---likely reflecting the stochastic nature of AGN fuelling---with a roughly power-law shape that rises toward lower $\lambda_{sBHAR}$, a steep cutoff at $\lambda_{sBHAR} \gtrsim 0.1-1$ (in Eddington equivalent units), and a turnover or flattening at $\lambda_{sBHAR} \lesssim 10^{-3}-10^{-2}$. We find that the probability of a star-forming galaxy hosting a moderate $\lambda_{sBHAR}$ AGN depends on stellar mass and evolves with redshift, shifting toward higher $\lambda_{sBHAR}$ at higher redshifts. This evolution is truncated at a point corresponding to the Eddington limit, indicating black holes may self-regulate their growth at high redshifts when copious gas is available. The probability of a quiescent galaxy hosting an AGN is generally lower than that of a star-forming galaxy, shows signs of suppression at the highest stellar masses, and evolves strongly with redshift. The AGN duty cycle in high-redshift ($z\gtrsim2$) quiescent galaxies thus reaches $\sim$20 per cent, comparable to the duty cycle in star-forming galaxies of equivalent stellar mass and redshift.

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

Using observations from the first two years of the MOSFIRE Deep Evolution Field (MOSDEF) survey, we study 13 AGN-driven outflows detected from a sample of 67 X-ray, IR and/or optically-selected AGN at $z \sim 2$. The AGN have bolometric luminosities of $\sim10^{44}-10^{46} ~\mathrm{erg~s^{-1}}$, including both quasars and moderate-luminosity AGN. We detect blueshifted, ionized gas outflows in the H$\beta$ , [OIII], H$\alpha$ ~and/or [NII] emission lines of $19\%$ of the AGN, while only 1.8\% of the MOSDEF galaxies have similarly-detected outflows. The outflow velocities span $\sim$300 to 1000 km s$^{-1}$. Eight of the 13 outflows are spatially extended on similar scales as the host galaxies, with spatial extents of 2.5 to 11.0 kpc. Outflows are detected uniformly across the star-forming main sequence, showing little trend with the host galaxy SFR. Line ratio diagnostics indicate that the outflowing gas is photoionized by the AGN. We do not find evidence for positive AGN feedback, in either our small MOSDEF sample or a much larger SDSS sample, using the BPT diagram. Given that a galaxy with an AGN is ten times more likely to have a detected outflow, the outflowing gas is photoionzed by the AGN, and estimates of the mass and energy outflow rates indicate that stellar feedback is insufficient to drive at least some of these outflows, they are very likely to be AGN-driven. The outflows have mass-loading factors of the order of unity, suggesting that they help regulate star formation in their host galaxies, though they may be insufficient to fully quench it.

We present the first spectroscopic measurement of multiple rest-frame optical emission lines at $z>4$. During the MOSFIRE Deep Evolution Field (MOSDEF) survey, we observed the galaxy GOODSN-17940 with the Keck I/MOSFIRE spectrograph. The K-band spectrum of GOODSN-17940 includes significant detections of the [OII]$\lambda\lambda 3726,3729$, [NeIII]$\lambda3869$, and H$\gamma$ emission lines and a tentative detection of H$\delta$, indicating $z_{\rm{spec}}=4.4121$. GOODSN-17940 is an actively star-forming $z>4$ galaxy based on its K-band spectrum and broadband spectral energy distribution. A significant excess relative to the surrounding continuum is present in the Spitzer/IRAC channel 1 photometry of GOODSN-17940, due primarily to strong H$\alpha$ emission with a rest-frame equivalent width of $\mbox{EW(H}\alpha)=1200$ \AA. Based on the assumption of $0.5 Z_{\odot}$ models and the Calzetti attenuation curve, GOODSN-17940 is characterized by $M_*=5.0^{+4.3}_{-0.2}\times 10^9 M_{\odot}$. The Balmer decrement inferred from H$\alpha$/H$\gamma$ is used to dust correct the H$\alpha$ emission, yielding $\mbox{SFR(H}\alpha)=320^{+190}_{-140} M_{\odot}\mbox{ yr}^{-1}$. These $M_*$ and SFR values place GOODSN-17940 an order of magnitude in SFR above the $z\sim 4$ star-forming "main sequence." Finally, we use the observed ratio of [NeIII]/[OII] to estimate the nebular oxygen abundance in GOODSN-17940, finding $\mbox{O/H}\sim 0.2 \mbox{ (O/H)}_{\odot}$. Combining our new [NeIII]/[OII] measurement with those from stacked spectra at $z\sim 0, 2, \mbox{ and } 3$, we show that GOODSN-17940 represents an extension to $z>4$ of the evolution towards higher [NeIII]/[OII] (i.e., lower $\mbox{O/H}$) at fixed stellar mass. It will be possible to perform the measurements presented here out to $z\sim 10$ using the James Webb Space Telescope.

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 investigate the mean star formation rates (SFRs) in the host galaxies of ~3000 optically selected QSOs from the SDSS survey within the Herschel-ATLAS fields, and a radio-luminous sub-sample, covering the redshift range of z = 0.2-2.5. Using WISE & Herschel photometry (12 - 500{\mu}m) we construct composite SEDs in bins of redshift and AGN luminosity. We perform SED fitting to measure the mean infrared luminosity due to star formation, removing the contamination from AGN emission. We find that the mean SFRs show a weak positive trend with increasing AGN luminosity. However, we demonstrate that the observed trend could be due to an increase in black hole (BH) mass (and a consequent increase of inferred stellar mass) with increasing AGN luminosity. We compare to a sample of X-ray selected AGN and find that the two populations have consistent mean SFRs when matched in AGN luminosity and redshift. On the basis of the available virial BH masses, and the evolving BH mass to stellar mass relationship, we find that the mean SFRs of our QSO sample are consistent with those of main sequence star-forming galaxies. Similarly, the radio-luminous QSOs have mean SFRs that are consistent with both the overall QSO sample and with star-forming galaxies on the main sequence. In conclusion, on average QSOs reside on the main sequence of star-forming galaxies, and the observed positive trend between the mean SFRs and AGN luminosity can be attributed to BH mass and redshift dependencies.

The Nuclear Spectroscopic Telescope Array (NuSTAR) provides an improvement in sensitivity at energies above 10 keV by two orders of magnitude over non-focusing satellites, making it possible to probe deeper into the Galaxy and Universe. Lansbury and collaborators recently completed a catalog of 497 sources serendipitously detected in the 3-24 keV band using 13 deg2 of NuSTAR coverage. Here, we report on an optical and X-ray study of 16 Galactic sources in the catalog. We identify eight of them as stars (but some or all could have binary companions), and use information from Gaia to report distances and X-ray luminosities for three of them. There are four CVs or CV candidates, and we argue that NuSTAR J233426-2343.9 is a relatively strong CV candidate based partly on an X-ray spectrum from XMM-Newton. NuSTAR J092418-3142.2, which is the brightest serendipitous source in the Lansbury catalog, and NuSTAR J073959-3147.8 are LMXB candidates, but it is also possible that these two sources are CVs. One of the sources is a known HMXB, and NuSTAR J105008-5958.8 is a new HMXB candidate, which has strong Balmer emission lines in its optical spectrum and a hard X-ray spectrum. We discuss the implications of finding these HMXBs for the surface density (logN-logS) and luminosity function of Galactic HMXBs. We conclude that, with the large fraction of unclassified sources in the Galactic plane detected by NuSTAR in the 8-24 keV band, there could be a significant population of low luminosity HMXBs.

The Advanced Telescope for High ENergy Astrophysics (Athena) is the X-ray observatory mission selected by ESA within its Cosmic Vision 2015-2025 programme to address the Hot and Energetic Universe scientific theme. The ESO-Athena Synergy Team (EAST) has been tasked to single out the potential scientific synergies between Athena and optical/near-infrared (NIR) and sub/mm ground based facilities, in particular those of ESO (i.e., the VLT and ELT, ALMA and APEX), by producing a White Paper to identify and develop the: 1. needs to access ESO ground-based facilities to achieve the formulated Athena science objectives; 2. needs to access Athena to achieve the formulated science objectives of ESO facilities contemporary to Athena; 3. science areas where the synergetic use of Athena and ESO facilities in the late 2020s will result in scientific added value. Community input to the process happened primarily via a dedicated ESO - Athena Synergy Workshop that took place on Sept. 14 - 16, 2016 at ESO, Garching. This White Paper presents the results of the EAST's work, sorted by synergy area, and deals with the following topics: 1. the Hot Universe: Early groups and clusters and their evolution, Physics of the Intracluster medium, Missing baryons in cosmic filaments; 2. the Energetic Universe: Supermassive black hole (SMBH) history, SMBH accretion disks, Active Galactic Nuclei feedback - Molecular outflows, Ultra-fast outflows, Accretion Physics, Transient Science; 3. Observatory Science: Star Formation, Stars. It then discusses the optical-NIR-sub-mm perspective by providing details on VLT/MOONS, the E-ELT instruments, in particular the MOS, VISTA/4MOST, the ESO and ALMA archives, future ALMA and ESO developments, and finally the (likely) ESO - Athena astronomical scene in the 2020s. (abridged)

This paper estimates the specific accretion-rate distribution of AGN using a sample of 4821 X-ray sources from both deep and shallow surveys. The specific accretion-rate distribution is defined as the probability of a galaxy with a given stellar mass and redshift hosting an active nucleus with a certain specific accretion rate. We find that the probability of a galaxy hosting an AGN increases with decreasing specific accretion rate. There is evidence that this trend reverses at low specific accretion rates, $\lambda<10^{-4}-10^{-3}$ (in Eddington units). There is also a break close to the Eddington limit, above which the probability of an accretion event decreases steeply. The specific accretion-rate distribution evolves such that the fraction of AGN among galaxies drops toward lower redshifts. This decrease in the AGN duty cycle is responsible for the strong evolution of the accretion density of the Universe from redshift $z\approx1-1.5$ to the present day. Our analysis also suggests that this evolution is accompanied by a decoupling of accretion events onto black holes from the formation of stars in galaxies. There is also evidence that at earlier times the relative probability of high vs low specific accretion-rate events among galaxies increases. We argue that this differential redshift evolution of the AGN duty cycle with respect to $\lambda$ produces the AGN downsizing trend, whereby luminous sources peak at earlier epochs compared to less luminous ones. Finally, we also find a stellar-mass dependence of the specific accretion-rate distribution, with more massive galaxies avoiding high specific accretion-rate events.

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 results from the MOSFIRE Deep Evolution Field (MOSDEF) survey on the identification, selection biases, and host galaxy properties of 55 X-ray, IR and optically-selected active galactic nuclei (AGN) at $1.4 < z < 3.8$. We obtain rest-frame optical spectra of galaxies and AGN and use the BPT diagram to identify optical AGN. We examine the uniqueness and overlap of the AGN identified at different wavelengths. There is a strong bias against identifying AGN at any wavelength in low mass galaxies, and an additional bias against identifying IR AGN in the most massive galaxies. AGN hosts span a wide range of star formation rate (SFR), similar to inactive galaxies once stellar mass selection effects are accounted for. However, we find (at $\sim 2-3\sigma$ significance) that IR AGN are in less dusty galaxies with relatively higher SFR and optical AGN in dusty galaxies with relatively lower SFR. X-ray AGN selection does not display a bias with host galaxy SFR. These results are consistent with those from larger studies at lower redshifts. Within star-forming galaxies, once selection biases are accounted for, we find AGN in galaxies with similar physical properties as inactive galaxies, with no evidence for AGN activity in particular types of galaxies. This is consistent with AGN being fueled stochastically in any star-forming host galaxy. We do not detect a significant correlation between SFR and AGN luminosity for individual AGN hosts, which may indicate the timescale difference between the growth of galaxies and their supermassive black holes.

We use deep Chandra imaging to measure the distribution of X-ray luminosities (L_X) for samples of star-forming galaxies as a function of stellar mass and redshift, using a Bayesian method to push below the nominal X-ray detection limits. Our luminosity distributions all show narrow peaks at L_X < 10^{42} erg/s that we associate with star formation, as opposed to AGN that are traced by a broad tail to higher L_X. Tracking the luminosity of these peaks as a function of stellar mass reveals an "X-ray main sequence" with a constant slope ~0.63 +/- 0.03 over 8.5 < log M*/Msun < 11.5 and 0.1 < z < 4, with a normalization that increases with redshift as (1+z)^{3.79+/-0.12}. We also compare the peak X-ray luminosities with UV-to-IR tracers of star formation rates (SFRs) to calibrate the scaling between L_X and SFR. We find that L_X \propto SFR^{0.83} x (1+z)^{1.3}, where the redshift evolution and non-linearity likely reflect changes in high-mass X-ray binary populations of star-forming galaxies. Using galaxies with a broader range of SFR, we also constrain a stellar-mass-dependent contribution to L_X, likely related to low-mass X-ray binaries. Using this calibration, we convert our X-ray main sequence to SFRs and measure a star-forming main sequence with a constant slope ~0.76+/-0.06 and a normalization that evolves with redshift as (1+z)^{2.95+/-0.33}. Based on the X-ray emission, there is no evidence for a break in the main sequence at high stellar masses, although we cannot rule out a turnover given the uncertainties in the scaling of L_X to SFR.

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.

The Wide Field Imager (WFI) is one of two instruments for the Advanced Telescope for High-ENergy Astrophysics (Athena). In this paper we summarise three of the many key science objectives for the WFI - the formation and growth of supermassive black holes, non-gravitational heating in clusters of galaxies, and spin measurements of stellar mass black holes - and describe their translation into the science requirements and ultimately instrument requirements. The WFI will be designed to provide excellent point source sensitivity and grasp for performing wide area surveys, surface brightness sensitivity, survey power, and absolute temperature and density calibration for in-depth studies of the outskirts of nearby clusters of galaxies and very good high-count rate capability, throughput, and low pile-up, paired with very good spectral resolution, for detailed explorations of bright Galactic compact objects.

We measure the clustering of X-ray, radio, and mid-IR-selected active galactic nuclei (AGN) at 0.2 < z < 1.2 using multi-wavelength imaging and spectroscopic redshifts from the PRIMUS and DEEP2 redshift surveys, covering 7 separate fields spanning ~10 square degrees. Using the cross-correlation of AGN with dense galaxy samples, we measure the clustering scale length and slope, as well as the bias, of AGN selected at different wavelengths. Similar to previous studies, we find that X-ray and radio AGN are more clustered than mid-IR-selected AGN. We further compare the clustering of each AGN sample with matched galaxy samples designed to have the same stellar mass, star formation rate, and redshift distributions as the AGN host galaxies and find no significant differences between their clustering properties. The observed differences in the clustering of AGN selected at different wavelengths can therefore be explained by the clustering differences of their host populations, which have different distributions in both stellar mass and star formation rate. Selection biases inherent in AGN selection, therefore, determine the clustering of observed AGN samples. We further find no significant difference between the clustering of obscured and unobscured AGN, using IRAC or WISE colors or X-ray hardness ratio.

The Cosmic X-ray Background (CXB) is the total emission from past accretion activity onto supermassive black holes in active galactic nuclei (AGN) and peaks in the hard X-ray band (30 keV). In this paper, we identify a significant selection effect operating on the CXB and flux-limited AGN surveys, and outline how they must depend heavily on the spin distribution of black holes. We show that, due to the higher radiative efficiency of rapidly-spinning black holes, they will be over-represented in the X-ray background, and therefore could be a dominant contributor to the CXB. Using a simple bimodal spin distribution, we demonstrate that only 15 per cent maximally-spinning AGN can produce 50 per cent of the CXB. We also illustrate that invoking a small population of maximally-spinning black holes in CXB synthesis models can reproduce the CXB peak without requiring large numbers of Compton-thick AGN. The spin bias is even more pronounced for flux-limited surveys: 7 per cent of sources with maximally-spinning black holes can produce half of the source counts. The detectability for maximum spin black holes can be further boosted in hard (>10 keV) X-rays by up to ~60 per cent due to pronounced ionised reflection, reducing the percentage of maximally spinning black holes required to produce half of the CXB or survey number counts further. A host of observations are consistent with an over-representation of high-spin black holes. Future NuSTAR and ASTRO-H hard X-ray surveys will provide the best constraints on the role of spin within the AGN population.

The active galactic nuclei X-ray luminosity function traces actively accreting supermassive black holes and is essential for the study of the properties of the active galactic nuclei (AGN) population, black hole evolution, and galaxy-black hole coevolution. Up to now, the AGN luminosity function has been estimated several times in soft (0.5-2 keV) and hard X-rays (2-10 keV). AGN selection in these energy ranges often suffers from identification and redshift incompleteness and, at the same time, photoelectric absorption can obscure a significant amount of the X-ray radiation. We estimate the evolution of the luminosity function in the 5-10 keV band, where we effectively avoid the absorbed part of the spectrum, rendering absorption corrections unnecessary up to NH=10^23 cm^-2. Our dataset is a compilation of six wide, and deep fields: MAXI, HBSS, XMM-COSMOS, Lockman Hole, XMM-CDFS, AEGIS-XD, Chandra-COSMOS, and Chandra-CDFS. This extensive sample of ~1110 AGN (0.01<z<4.0, 41<log L_x<46) is 98% redshift complete with 68% spectroscopic redshifts. We use Bayesian analysis to select the best parametric model from simple pure luminosity and pure density evolution to more complicated luminosity and density evolution and luminosity-dependent density evolution. We estimate the model parameters that describe best our dataset separately for each survey and for the combined sample. We show that, according to Bayesian model selection, the preferred model for our dataset is the luminosity-dependent density evolution (LDDE). Our estimation of the AGN luminosity function does not require any assumption on the AGN absorption and is in good agreement with previous works in the 2-10 keV energy band based on X-ray hardness ratios to model the absorption in AGN up to redshift three. Our sample does not show evidence of a rapid decline of the AGN luminosity function up to redshift four. [abridged]

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

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 examine the host morphologies of heavily obscured active galactic nuclei (AGN) at $z\sim1$ to test whether obscured supermassive black hole growth at this epoch is preferentially linked to galaxy mergers. Our sample consists of 154 obscured AGN with $N_{\rm H}>10^{23.5}$ cm$^{-2}$ and $z<1.5$. Using visual classifications, we compare the morphologies of these AGN to control samples of moderately obscured ($10^{22}$ cm$^{-2}$ $<N_{\rm H}< 10^{23.5}$ cm$^{-2}$) and unobscured ($N_{\rm H}<10^{22}$ cm$^{-2}$) AGN. These control AGN are matched in redshift and intrinsic X-ray luminosity to our heavily obscured AGN. We find that heavily obscured AGN at z~1 are twice as likely to be hosted by late-type galaxies relative to unobscured AGN ($65.3^{+4.1}_{-4.6}\%$ vs $34.5^{+2.9}_{-2.7}\%$) and three times as likely to exhibit merger or interaction signatures ($21.5^{+4.2}_{-3.3}\%$ vs $7.8^{+1.9}_{-1.3}\%$). The increased merger fraction is significant at the 3.8$\sigma$ level. We also find that the incidence of point-like morphologies is inversely proportional to obscuration. If we exclude all point sources and consider only extended hosts, we find the correlation between merger fraction and obscuration is still evident, however at a reduced statistical significance ($2.5\sigma$). The fact that we observe a different disk/spheroid fraction versus obscuration indicates that viewing angle cannot be the only thing differentiating our three AGN samples, as a simple unification model would suggest. The increased fraction of disturbed morphologies with obscuration supports an evolutionary scenario, in which Compton-thick AGN are a distinct phase of obscured SMBH growth following a merger/interaction event. Our findings also suggest that some of the merger-triggered SMBH growth predicted by recent AGN fueling models may be hidden among the heavily obscured, Compton-thick population.

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.

We combine deep X-ray survey data from the Chandra observatory and the wide-area/shallow XMM-XXL field to estimate the AGN X-ray luminosity function in the redshift range z=3-5. The sample consists of nearly 340 sources with either photometric (212) or spectroscopic (128) redshift in the above range. The combination of deep and shallow survey fields provides a luminosity baseline of three orders of magnitude, Lx(2-10keV)~1e43-1e46erg/s at z>3. We follow a Bayesian approach to determine the binned AGN space density and explore their evolution in a model-independent way. Our methodology accounts for Poisson errors in the determination of X-ray fluxes and uncertainties in photometric redshift estimates. We demonstrate that the latter is essential for unbiased measurement of space densities. We find that the AGN X-ray luminosity function evolves strongly between the redshift intervals z=3-4 and z=4-5. There is also suggestive evidence that the amplitude of this evolution is luminosity dependent. The space density of AGN with Lx<1e45erg/s drops by a factor of 5 between the redshift intervals above, while the evolution of brighter AGN appears to be milder. Comparison of our X-ray luminosity function with that of UV/optical selected QSOs at similar redshifts shows broad agreement at bright luminosities, Lx>1e45erg/s. The faint-end slope of UV/optical luminosity functions however, is steeper than for X-ray selected AGN. This implies that the type-I AGN fraction increases with decreasing luminosity at z>3, opposite to trends established at lower redshift. We also assess the significance of AGN in keeping the hydrogen ionised at high redshift. Our X-ray luminosity function yields ionising photon rate densities that are insufficient to keep the Universe ionised at redshift z>4. A source of uncertainty in this calculation is the escape fraction of UV photons for X-ray selected AGN.

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

The intrinsic column density (NH) distribution of quasars is poorly known. At the high obscuration end of the quasar population and for redshifts z<1, the X-ray spectra can only be reliably characterized using broad-band measurements which extend to energies above 10 keV. Using the hard X-ray observatory NuSTAR, along with archival Chandra and XMM-Newton data, we study the broad-band X-ray spectra of nine optically selected (from the SDSS), candidate Compton-thick (NH > 1.5e24 cm^-2) type 2 quasars (CTQSO2s); five new NuSTAR observations are reported herein, and four have been previously published. The candidate CTQSO2s lie at z<0.5, have observed [OIII] luminosities in the range 8.4 < log (L_[OIII]/L_solar) < 9.6, and show evidence for extreme, Compton-thick absorption when indirect absorption diagnostics are considered. Amongst the nine candidate CTQSO2s, five are detected by NuSTAR in the high energy (8-24 keV) band: two are weakly detected at the ~ 3 sigma confidence level and three are strongly detected with sufficient counts for spectral modeling (>~ 90 net source counts at 8-24 keV). For these NuSTAR-detected sources direct (i.e., X-ray spectral) constraints on the intrinsic AGN properties are feasible, and we measure column densities ~2.5-1600 times higher and intrinsic (unabsorbed) X-ray luminosities ~10-70 times higher than pre-NuSTAR constraints from Chandra and XMM-Newton. Assuming the NuSTAR-detected type 2 quasars are representative of other Compton-thick candidates, we make a correction to the NH distribution for optically selected type 2 quasars as measured by Chandra and XMM-Newton for 39 objects. With this approach, we predict a Compton-thick fraction of f_CT = 36^{+14}_{-12} %, although higher fractions (up to 76%) are possible if indirect absorption diagnostics are assumed to be reliable.

We study the evidence for a connection between active galactic nuclei (AGN) fueling and star formation by investigating the relationship between the X-ray luminosities of AGN and the star formation rates (SFRs) of their host galaxies. We identify a sample of 309 AGN with $10^{41}<L_\mathrm{X}<10^{44} $ erg s$^{-1}$ at $0.2 < z < 1.2$ in the PRIMUS redshift survey. We find AGN in galaxies with a wide range of SFR at a given $L_X$. We do not find a significant correlation between SFR and the observed instantaneous $L_X$ for star forming AGN host galaxies. However, there is a weak but significant correlation between the mean $L_\mathrm{X}$ and SFR of detected AGN in star forming galaxies, which likely reflects that $L_\mathrm{X}$ varies on shorter timescales than SFR. We find no correlation between stellar mass and $L_\mathrm{X}$ within the AGN population. Within both populations of star forming and quiescent galaxies, we find a similar power-law distribution in the probability of hosting an AGN as a function of specific accretion rate. Furthermore, at a given stellar mass, we find a star forming galaxy $\sim2-3$ more likely than a quiescent galaxy to host an AGN of a given specific accretion rate. The probability of a galaxy hosting an AGN is constant across the main sequence of star formation. These results indicate that there is an underlying connection between star formation and the presence of AGN, but AGN are often hosted by quiescent galaxies.

We present new measurements of the evolution of the X-ray luminosity functions (XLFs) of unabsorbed and absorbed Active Galactic Nuclei (AGNs) out to z~5. We construct samples containing 2957 sources detected at hard (2-7 keV) X-ray energies and 4351 sources detected at soft (0.5-2 keV) energies from a compilation of Chandra surveys supplemented by wide-area surveys from ASCA} and ROSAT. We consider the hard and soft X-ray samples separately and find that the XLF based on either (initially neglecting absorption effects) is best described by a new flexible model parametrization where the break luminosity, normalization and faint-end slope all evolve with redshift. We then incorporate absorption effects, separately modelling the evolution of the XLFs of unabsorbed ($20<\log N_\mathrm{H}<22$) and absorbed ($22<\log N_\mathrm{H}<24$) AGNs, seeking a model that can reconcile both the hard- and soft-band samples. We find that the absorbed AGN XLF has a lower break luminosity, a higher normalization, and a steeper faint-end slope than the unabsorbed AGN XLF out to z~2. Hence, absorbed AGNs dominate at low luminosities, with the absorbed fraction falling rapidly as luminosity increases. Both XLFs undergo strong luminosity evolution which shifts the transition in the absorbed fraction to higher luminosities at higher redshifts. The evolution in the shape of the total XLF is primarily driven by the changing mix of unabsorbed and absorbed populations.

Using data from the All Wavelength Extended Groth Strip International Survey (AEGIS) we statistically detect the extended X-ray emission in the interstellar medium (ISM)/intra-cluster medium (ICM) in both active and normal galaxies at 0.3 <= z <= 1.3. For both active galactic nuclei (AGN) host galaxy and normal galaxy samples that are matched in restframe color, luminosity, and redshift distribution, we tentatively detect excess X-ray emission at scales of 1--10 arcsec at a few sigma significance in the surface brightness profiles. The exact significance of this detection is sensitive to the true characterization of Chandra's point spread function. The observed excess in the surface brightness profiles is suggestive of lower extended emission in AGN hosts compared to normal galaxies. This is qualitatively similar to theoretical predictions of the X-ray surface brightness profile from AGN feedback models, where feedback from AGN is likely to evacuate the gas from the center of the galaxy/cluster. We propose that AGN that are intrinsically under-luminous in X-rays, but have equivalent bolometric luminosities to our sources will be the ideal sample to study more robustly the effect of AGN feedback on diffuse ISM/ICM gas.

Using data from the DEEP2 galaxy redshift survey and the All Wavelength Extended Groth Strip International Survey we obtain stacked X-ray maps of galaxies at 0.7 < z < 1.0 as a function of stellar mass. We compute the total X-ray counts of these galaxies and show that in the soft band (0.5--2,kev) there exists a significant correlation between galaxy X-ray counts and stellar mass at these redshifts. The best-fit relation between X-ray counts and stellar mass can be characterized by a power law with a slope of 0.58 +/- 0.1. We do not find any correlation between stellar mass and X-ray luminosities in the hard (2--7,kev) and ultra-hard (4--7,kev) bands. The derived hardness ratios of our galaxies suggest that the X-ray emission is degenerate between two spectral models, namely point-like power-law emission and extended plasma emission in the interstellar medium. This is similar to what has been observed in low redshift galaxies. Using a simple spectral model where half of the emission comes from power-law sources and the other half from the extended hot halo we derive the X-ray luminosities of our galaxies. The soft X-ray luminosities of our galaxies lie in the range 10^39-8x10^40, ergs/s. Dividing our galaxy sample by the criteria U-B > 1, we find no evidence that our results for X-ray scaling relations depend on optical color.

In this paper we present the MOSFIRE Deep Evolution Field (MOSDEF) survey. The MOSDEF survey aims to obtain moderate-resolution (R=3000-3650) rest-frame optical spectra (~3700-7000 Angstrom) for ~1500 galaxies at 1.37<z<3.80 in three well-studied CANDELS fields: AEGIS, COSMOS, and GOODS-N. Targets are selected in three redshift intervals: 1.37<z<1.70, 2.09<z<2.61, and 2.95<z<3.80, down to fixed H_AB (F160W) magnitudes of 24.0, 24.5 and 25.0, respectively, using the photometric and spectroscopic catalogs from the 3D-HST survey. We target both strong nebular emission lines (e.g., [OII], Hbeta, [OIII], 5008, Halpha, [NII], and [SII]) and stellar continuum and absorption features (e.g., Balmer lines, Ca-II H and K, Mgb, 4000 Angstrom break). Here we present an overview of our survey, the observational strategy, the data reduction and analysis, and the sample characteristics based on spectra obtained during the first 24 nights. To date, we have completed 21 masks, obtaining spectra for 591 galaxies. For ~80% of the targets we derive a robust redshift from either emission or absorption lines. In addition, we confirm 55 additional galaxies, which were serendipitously detected. The MOSDEF galaxy sample includes unobscured star-forming, dusty star-forming, and quiescent galaxies and spans a wide range in stellar mass (~10^9-10^11.5 Msol) and star formation rate (~10^0-10^3 Msol/yr). The spectroscopically confirmed sample is roughly representative of an H-band limited galaxy sample at these redshifts. With its large sample size, broad diversity in galaxy properties, and wealth of available ancillary data, MOSDEF will transform our understanding of the stellar, gaseous, metal, dust, and black hole content of galaxies during the time when the universe was most active.

We aim to constrain the evolution of AGN as a function of obscuration using an X-ray selected sample of $\sim2000$ AGN from a multi-tiered survey including the CDFS, AEGIS-XD, COSMOS and XMM-XXL fields. The spectra of individual X-ray sources are analysed using a Bayesian methodology with a physically realistic model to infer the posterior distribution of the hydrogen column density and intrinsic X-ray luminosity. We develop a novel non-parametric method which allows us to robustly infer the distribution of the AGN population in X-ray luminosity, redshift and obscuring column density, relying only on minimal smoothness assumptions. Our analysis properly incorporates uncertainties from low count spectra, photometric redshift measurements, association incompleteness and the limited sample size. We find that obscured AGN with $N_{H}>{\rm 10^{22}\, cm^{-2}}$ account for ${77}^{+4}_{-5}\%$ of the number density and luminosity density of the accretion SMBH population with $L_{{\rm X}}>10^{43}\text{ erg/s}$, averaged over cosmic time. Compton-thick AGN account for approximately half the number and luminosity density of the obscured population, and ${38}^{+8}_{-7}\%$ of the total. We also find evidence that the evolution is obscuration-dependent, with the strongest evolution around $N_{H}\thickapprox10^{23}\text{ cm}^{-2}$. We highlight this by measuring the obscured fraction in Compton-thin AGN, which increases towards $z\sim3$, where it is $25\%$ higher than the local value. In contrast the fraction of Compton-thick AGN is consistent with being constant at $\approx35\%$, independent of redshift and accretion luminosity. We discuss our findings in the context of existing models and conclude that the observed evolution is to first order a side-effect of anti-hierarchical growth.

We present the results of deep \chandra\ imaging of the central region of the Extended Groth Strip, the AEGIS-X Deep (AEGIS-XD) survey. When combined with previous \chandra\ observations of a wider area of the strip, AEGIS-X Wide (AEGIS-XW; Laird et~al. 2009), these provide data to a nominal exposure depth of 800ks in the three central ACIS-I fields, a region of approximately $0.29$~deg$^{2}$. This is currently the third deepest X-ray survey in existence, a factor $\sim 2-3$ shallower than the Chandra Deep Fields (CDFs) but over an area $\sim 3$ times greater than each CDF. We present a catalogue of 937 point sources detected in the deep \chandra\ observations. We present identifications of our X-ray sources from deep ground-based, Spitzer, GALEX and HST imaging. Using a likelihood ratio analysis, we associate multi band counterparts for 929/937 of our X-ray sources, with an estimated 95~\% reliability, making the identification completeness approximately 94~\% in a statistical sense. Reliable spectroscopic redshifts for 353 of our X-ray sources are provided predominantly from Keck (DEEP2/3) and MMT Hectospec, so the current spectroscopic completeness is $\sim 38$~per cent. For the remainder of the X-ray sources, we compute photometric redshifts based on multi-band photometry in up to 35 bands from the UV to mid-IR. Particular attention is given to the fact that the vast majority the X-ray sources are AGN and require hybrid templates. Our photometric redshifts have mean accuracy of $\sigma=0.04$ and an outlier fraction of approximately 5\%, reaching $\sigma=0.03$ with less than 4\% outliers in the area covered by CANDELS . The X-ray, multi-wavelength photometry and redshift catalogues are made publicly available.

We present results from the MOSFIRE Deep Evolution Field (MOSDEF) survey on rest-frame optical AGN identification and completeness at z~2.3. With our sample of 50 galaxies and 10 X-ray and IR-selected AGN with measured H-beta, [OIII], H-alpha, and [NII] emission lines, we investigate the location of AGN in the BPT, MEx (mass-excitation), and CEx (color-excitation) diagrams. We find that the BPT diagram works well to identify AGN at z~2.3 and that the z~0 AGN/star-forming galaxy classifications do not need to shift substantially at z~2.3 to robustly separate these populations. However, the MEx diagram fails to identify all of the AGN identified in the BPT diagram, and the CEx diagram is substantially contaminated at high redshift. We further show that AGN samples selected using the BPT diagram have selection biases in terms of both host stellar mass and stellar population, in that AGN in low mass and/or high specific star formation rate galaxies are difficult to identify using the BPT diagram. These selection biases become increasingly severe at high redshift, such that optically-selected AGN samples at high redshift will necessarily be incomplete. We also find that the gas in the narrow-line region appears to be more enriched than gas in the host galaxy for at least some MOSDEF AGN. However, AGN at z~2 are generally less enriched than local AGN with the same host stellar mass.

We analyse the stellar populations in the host galaxies of 53 X-ray selected optically dull active galactic nuclei (AGN) at 0.34<z<1.07 with ultra-deep (m=26.5) optical medium-band (R~50) photometry from the Survey for High-z Absorption Red and Dead Sources (SHARDS). The spectral resolution of SHARDS allows us to consistently measure the strength of the 4000 AA break, Dn(4000), a reliable age indicator for stellar populations. We confirm that most X-ray selected moderate-luminosity AGN (L_X<10^44 erg/s) are hosted by massive galaxies (typically M*>10^10.5 M_sun) and that the observed fraction of galaxies hosting an AGN increases with the stellar mass. A careful selection of random control samples of inactive galaxies allows us to remove the stellar mass and redshift dependencies of the AGN fraction to explore trends with several stellar age indicators. We find no significant differences in the distribution of the rest-frame U-V colour for AGN hosts and inactive galaxies, in agreement with previous results. However, we find significantly shallower 4000 AA breaks in AGN hosts, indicative of younger stellar populations. With the help of a model-independent determination of the extinction, we obtain extinction-corrected U-V colours and light-weighted average stellar ages. We find that AGN hosts have younger stellar populations and higher extinction compared to inactive galaxies with the same stellar mass and at the same redshift. We find a highly significant excess of AGN hosts with Dn(4000)~1.4 and light weighted average stellar ages of 300-500 Myr, as well as a deficit of AGN in intrinsic red galaxies. We interpret failure in recognising these trends in previous studies as a consequence of the balancing effect in observed colours of the age-extinction degeneracy.

We present measurements of the luminosity and color-dependence of galaxy clustering at 0.2<z<1.0 in the PRIsm MUlti-object Survey (PRIMUS). We quantify the clustering with the redshift-space and projected two-point correlation functions, xi(rp,pi) and wp(rp), using volume-limited samples constructed from a parent sample of over 130,000 galaxies with robust redshifts in seven independent fields covering 9 sq. deg. of sky. We quantify how the scale-dependent clustering amplitude increases with increasing luminosity and redder color, with relatively small errors over large volumes. We find that red galaxies have stronger small-scale (0.1<rp<1 Mpc/h) clustering and steeper correlation functions compared to blue galaxies, as well as a strong color dependent clustering within the red sequence alone. We interpret our measured clustering trends in terms of galaxy bias and obtain values between b_gal=0.9-2.5, quantifying how galaxies are biased tracers of dark matter depending on their luminosity and color. We also interpret the color dependence with mock catalogs, and find that the clustering of blue galaxies is nearly constant with color, while redder galaxies have stronger clustering in the one-halo term due to a higher satellite galaxy fraction. In addition, we measure the evolution of the clustering strength and bias, and we do not detect statistically significant departures from passive evolution. We argue that the luminosity- and color-environment (or halo mass) relations of galaxies have not significantly evolved since z=1. Finally, using jackknife subsampling methods, we find that sampling fluctuations are important and that the COSMOS field is generally an outlier, due to having more overdense structures than other fields; we find that 'cosmic variance' can be a significant source of uncertainty for high-redshift clustering measurements.

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.

Chandra data in the COSMOS, AEGIS-XD and 4Ms CDFS are combined with optical/near-IR photometry to determine the rest-frame U-V vs V-J colours of X-ray AGN hosts at mean redshifts 0.40 and 0.85. This combination of colours (UVJ) provides an efficient means of separating quiescent from star-forming, including dust reddened, galaxies. Morphological information emphasises differences between AGN split by their UVJ colours. AGN in quiescent galaxies are dominated by spheroids, while star-forming hosts are split between bulges and disks. The UVJ diagram of AGN hosts is then used to set limits on the accretion density associated with evolved and star-forming systems. Most of the black hole growth since z~1 is associated with star-forming hosts. Nevertheless, ~15-20% of the X-ray luminosity density since z~1, is taking place in the quiescent region of the UVJ diagram. For the z~0.40 subsample, there is tentative evidence (2sigma significance), that AGN split by their UVJ colours differ in Eddington ratio. AGN in star-forming hosts dominate at high Eddington ratios, while AGN in quiescent hosts become increasingly important as a fraction of the total population toward low Eddington ratios. At higher redshift, z~0.8, such differences are significant at the 2sigma level only at Eddington ratios >1e-3. These findings are consistent with scenarios in which diverse accretion modes are responsible for the build-up of SMBHs at the centres of galaxies. We compare our results with the GALFORM semi-analytic model, which postulates two black hole fuelling modes, the first linked to star-formation and the second occuring in passive galaxies. GALFORM predicts a larger fraction of black hole growth in quiescent galaxies at z<1, compared to the data. Relaxing the strong assumption of the model that passive AGN hosts have zero star-formation rate could reconcile this disagreement.

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)

The Advanced Telescope for High-energy Astrophysics (Athena+) is being proposed to ESA as the L2 mission (for a launch in 2028) and is specifically designed to answer two of the most pressing questions for astrophysics in the forthcoming decade: How did ordinary matter assemble into the large scale structures we see today? and how do black holes grow and shape the Universe? For addressing these two issues, Athena+ will provide transformational capabilities in terms of angular resolution, effective area, spectral resolution, grasp, that will make it the most powerful X-ray observatory ever flown. Such an observatory, when opened to the astronomical community, will be used for virtually all classes of astrophysical objects, from high-z gamma-ray bursts to the closest planets in our solar neighborhood. In this paper, we briefly review the core science objectives of Athena+, present the science requirements and the foreseen implementation of the mission, and illustrate its transformational capabilities compared to existing facilities.

We present a study of Spitzer/IRAC and X-ray active galactic nuclei (AGNs) selection techniques in order to quantify the overlap, uniqueness, contamination, and completeness of each. We investigate how the overlap and possible contamination of the samples depends on the IR and X-ray depths. We use Spitzer/IRAC imaging, Chandra and XMM X-ray imaging, and PRism MUlti-object Survey (PRIMUS) spectroscopic redshifts to construct galaxy and AGN samples at 0.2<z<1.2 over 8 deg^2. We construct samples over a wide range of IRAC flux limits (SWIRE to GOODS depth) and X-ray flux limits (10 ks to 2 Ms). We compare IR-AGN samples defined using the IRAC color selection of Stern et al. and Donley et al. with X-ray detected AGN samples. For roughly similar depth IR and X-ray surveys, we find that ~75% of IR-AGN are identified as X-ray AGN. This fraction increases to ~90% when comparing against the deepest X-ray data, indicating that only ~10% of IR-selected AGN may be heavily obscured. The IR-AGN selection proposed by Stern et al. suffers from contamination by star-forming galaxies at various redshifts when using deeper IR data, though the selection technique works well for shallow IR data. While similar overall, the IR-AGN samples preferentially contain more luminous AGN, while the X-ray AGN samples preferentially contain lower specific accretion rate AGN, where the host galaxy light dominates at IR wavelengths. The host galaxy populations of the IR and X-ray AGN samples have similar restframe colors and stellar masses; both selections identify AGN in blue, star-forming and red, quiescent galaxies.

We present an observationally motivated model to connect the AGN and galaxy populations at 0.2<z<1.0 and predict the AGN X-ray luminosity function (XLF). We start with measurements of the stellar mass function of galaxies (from the Prism Multi-object Survey) and populate galaxies with AGNs using models for the probability of a galaxy hosting an AGN as a function of specific accretion rate. Our model is based on measurements indicating that the specific accretion rate distribution is a universal function across a wide range of host stellar mass with slope gamma_1 = -0.65 and an overall normalization that evolves with redshift. We test several simple assumptions to extend this model to high specific accretion rates (beyond the measurements) and compare the predictions for the XLF with the observed data. We find good agreement with a model that allows for a break in the specific accretion rate distribution at a point corresponding to the Eddington limit, a steep power-law tail to super-Eddington ratios with slope gamma_2 = -2.1 +0.3 -0.5, and a scatter of 0.38 dex in the scaling between black hole and host stellar mass. Our results show that samples of low luminosity AGNs are dominated by moderately massive galaxies (M* ~ 10^{10-11} M_sun) growing with a wide range of accretion rates due to the shape of the galaxy stellar mass function rather than a preference for AGN activity at a particular stellar mass. Luminous AGNs may be a severely skewed population with elevated black hole masses relative to their host galaxies and in rare phases of rapid accretion.

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.

A crucial challenge in astrophysics over the coming decades will be to understand the origins of supermassive black holes (SMBHs) that lie at the centres of most, if not all, galaxies. The processes responsible for the initial formation of these SMBHs and their early growth via accretion - when they are seen as Active Galactic Nuclei (AGN) - remain unknown. To address this challenge, we must identify low luminosity and obscured z>6 AGNs, which represent the bulk of early SMBH growth. Sensitive X-ray observations are a unique signpost of accretion activity, uncontaminated by star formation processes, which prevent reliable AGN identification at other wavelengths (e.g. optical, infrared). The Athena+ Wide Field Imager will enable X-ray surveys to be carried out two orders of magnitude faster than with Chandra or XMM-Newton, opening a new discovery space and identifying over 400 z>6 AGN, including obscured sources. Athena+ will also play a fundamental role to enhance the scientific return of future multiwavelength facilities that will probe the physical conditions within the host galaxies of early SMBHs, which is vital for understanding how SMBHs form, what fuels their subsequent growth, and to assess their impact on the early Universe. Follow-up of samples of z>6 galaxies with the Athena+ X-ray Integral Field Unit could also reveal the presence of highly obscured AGNs, thanks to the detection of strong iron lines. Thus, Athena+ will enable the first quantitative measurements of the extent and distribution of SMBH accretion in the early Universe.

Observations in the last decade have provided strong evidence that the growth of supermassive black holes at the centres of galaxies is among the most influential processes in galaxy evolution. Open questions that relate to our current understanding of black hole growth and its relation to the build-up of galaxies at redshifts z=1-4, when most black holes and stars we see in present-day galaxies were put in place, include: what is the nature of AGN feedback and whether it plays a significant role in the evolution of galaxies? what is the dominant population of accreting AGN at that critical epoch? is it dominated by obscured objects as required by many current observations and models? The Athena+ mission concept will provide the technological leap required for a breakthrough in our understanding of AGN and galaxy evolution at the heyday of the Universe. The high throughput of Athena+ will allow the systematic study of the incidence, nature and energetics of AGN feedback processes to z~4 via the identification and measurement of blue-shifted X-ray absorption lines. The excellent survey and spectral capabilities of the Athena+ Wide Field Imager will complete the census of black hole growth by yielding samples of up to 100 times larger than is currently possible of the most heavily obscured, including Compton thick, AGN to redshifts z~4. The demographics of this population relative to their hosts is fundamental for understanding how major black hole growth events relate to the build-up of galaxies.

This White Paper, submitted to the recent ESA call for science themes to define its future large missions, advocates the need for a transformational leap in our understanding of two key questions in astrophysics: 1) How does ordinary matter assemble into the large scale structures that we see today? 2) How do black holes grow and shape the Universe? Hot gas in clusters, groups and the intergalactic medium dominates the baryonic content of the local Universe. To understand the astrophysical processes responsible for the formation and assembly of these large structures, it is necessary to measure their physical properties and evolution. This requires spatially resolved X-ray spectroscopy with a factor 10 increase in both telescope throughput and spatial resolving power compared to currently planned facilities. Feedback from supermassive black holes is an essential ingredient in this process and in most galaxy evolution models, but it is not well understood. X-ray observations can uniquely reveal the mechanisms launching winds close to black holes and determine the coupling of the energy and matter flows on larger scales. Due to the effects of feedback, a complete understanding of galaxy evolution requires knowledge of the obscured growth of supermassive black holes through cosmic time, out to the redshifts where the first galaxies form. X-ray emission is the most reliable way to reveal accreting black holes, but deep survey speed must improve by a factor ~100 over current facilities to perform a full census into the early Universe. The Advanced Telescope for High Energy Astrophysics (Athena+) mission provides the necessary performance (e.g. angular resolution, spectral resolution, survey grasp) to address these questions and revolutionize our understanding of the Hot and Energetic Universe. These capabilities will also provide a powerful observatory to be used in all areas of astrophysics.

We measure the evolution of the stellar mass function (SMF) from z=0-1 using multi-wavelength imaging and spectroscopic redshifts from the PRism MUlti-object Survey (PRIMUS) and the Sloan Digital Sky Survey (SDSS). From PRIMUS we construct an i<23 flux-limited sample of ~40,000 galaxies at z=0.2-1.0 over five fields totaling ~5.5 deg^2, and from the SDSS we select ~170,000 galaxies at z=0.01-0.2 that we analyze consistently with respect to PRIMUS to minimize systematic errors in our evolutionary measurements. We find that the SMF of all galaxies evolves relatively little since z=1, although we do find evidence for mass assembly downsizing; we measure a ~30% increase in the number density of ~10^10 Msun galaxies since z~0.6, and a <10% change in the number density of all >10^11 Msun galaxies since z~1. Dividing the sample into star-forming and quiescent using an evolving cut in specific star-formation rate, we find that the number density of ~10^10 Msun star-forming galaxies stays relatively constant since z~0.6, whereas the space-density of >10^11 Msun star-forming galaxies decreases by ~50% between z~1 and z~0. Meanwhile, the number density of ~10^10 Msun quiescent galaxies increases steeply towards low redshift, by a factor of ~2-3 since z~0.6, while the number of massive quiescent galaxies remains approximately constant since z~1. These results suggest that the rate at which star-forming galaxies are quenched increases with decreasing stellar mass, but that the bulk of the stellar mass buildup within the quiescent population occurs around ~10^10.8 Msun. In addition, we conclude that mergers do not appear to be a dominant channel for the stellar mass buildup of galaxies at z<1, even among massive (>10^11 Msun) quiescent galaxies.

The PRIsm MUti-object Survey (PRIMUS) is a spectroscopic galaxy redshift survey to z~1 completed with a low-dispersion prism and slitmasks allowing for simultaneous observations of ~2,500 objects over 0.18 square degrees. The final PRIMUS catalog includes ~130,000 robust redshifts over 9.1 sq. deg. In this paper, we summarize the PRIMUS observational strategy and present the data reduction details used to measure redshifts, redshift precision, and survey completeness. The survey motivation, observational techniques, fields, target selection, slitmask design, and observations are presented in Coil et al 2010. Comparisons to existing higher-resolution spectroscopic measurements show a typical precision of sigma_z/(1+z)=0.005. PRIMUS, both in area and number of redshifts, is the largest faint galaxy redshift survey completed to date and is allowing for precise measurements of the relationship between AGNs and their hosts, the effects of environment on galaxy evolution, and the build up of galactic systems over the latter half of cosmic history.

We present evidence that the incidence of active galactic nuclei (AGNs) and the distribution of their accretion rates do not depend on the stellar masses of their host galaxies, contrary to previous studies. We use hard (2-10 keV) X-ray data from three extragalactic fields (XMM-LSS, COSMOS and ELAIS-S1) with redshifts from the Prism Multi-object Survey to identify 242 AGNs with L_{2-10 keV}=10^{42-44} erg /s within a parent sample of ~25,000 galaxies at 0.2<z<1.0 over ~3.4 deg^2 and to i~23. We find that although the fraction of galaxies hosting an AGN at fixed X-ray luminosity rises strongly with stellar mass, the distribution of X-ray luminosities is independent of mass. Furthermore, we show that the probability that a galaxy will host an AGN can be defined by a universal Eddington ratio distribution that is independent of the host galaxy stellar mass and has a power-law shape with slope -0.65. These results demonstrate that AGNs are prevalent at all stellar masses in the range 9.5<log M_*/M_sun<12 and that the same physical processes regulate AGN activity in all galaxies in this stellar mass range. While a higher AGN fraction may be observed in massive galaxies, this is a selection effect related to the underlying Eddington ratio distribution. We also find that the AGN fraction drops rapidly between z~1 and the present day and is moderately enhanced (factor~2) in galaxies with blue or green optical colors. Consequently, while AGN activity and star formation appear to be globally correlated, we do not find evidence that the presence of an AGN is related to the quenching of star formation or the color transformation of galaxies.

We present Keck/LRIS-B spectra for a sample of ten AEGIS X-ray AGN host galaxies and thirteen post-starburst galaxies from SDSS and DEEP2 at 0.2<z<0.8 in order to investigate the presence, properties, and influence of outflowing galactic winds at intermediate redshifts. We focus on galaxies that either host a low-luminosity AGN or have recently had their star formation quenched to test whether these galaxies have winds of sufficient velocity to potentially clear gas from the galaxy. We find, using absorption features of Fe II, Mg II, and Mg I, that six of the ten (60%) X-ray AGN host galaxies and four of the thirteen (31%) post-starburst galaxies have outflowing galactic winds, with typical velocities of ~200 km/s. We additionally find that most of the galaxies in our sample show line emission, possibly from the wind, in either Fe II* or Mg II. A total of 100% of our X-ray AGN host sample (including four red sequence galaxies) and 77% of our post-starburst sample has either blueshifted absorption or line emission. Several K+A galaxies have small amounts of cool gas absorption at the systemic velocity, indicating that not all of the cool gas has been expelled. We conclude that while outflowing galactic winds are common in both X-ray low-luminosity AGN host galaxies and post-starburst galaxies at intermediate redshifts, the winds are likely driven by supernovae (as opposed to AGN) and do not appear to have sufficiently high velocities to quench star formation in these galaxies.

We present the PRIsm MUlti-object Survey (PRIMUS), a spectroscopic faint galaxy redshift survey to z~1. PRIMUS uses a low-dispersion prism and slitmasks to observe ~2,500 objects at once in a 0.18 deg^2 field of view, using the IMACS camera on the Magellan I Baade 6.5m telescope at Las Campanas Observatory. PRIMUS covers a total of 9.1 deg^2 of sky to a depth of i_AB~23.5 in seven different deep, multi-wavelength fields that have coverage from GALEX, Spitzer and either XMM or Chandra, as well as multiple-band optical and near-IR coverage. PRIMUS includes ~130,000 robust redshifts of unique objects with a redshift precision of dz/(1+z)~0.005. The redshift distribution peaks at z=0.6 and extends to z=1.2 for galaxies and z=5 for broad-line AGN. The motivation, observational techniques, fields, target selection, slitmask design, and observations are presented here, with a brief summary of the redshift precision; a companion paper presents the data reduction, redshift fitting, redshift confidence, and survey completeness. PRIMUS is the largest faint galaxy survey undertaken to date. The high targeting fraction (~80%) and large survey size will allow for precise measures of galaxy properties and large-scale structure to z~1.

We present the results of spectroscopic observations in the GOODS-N field completed using DEIMOS on the Keck II telescope as part of the DEEP3 Galaxy Redshift Survey. Observations of 370 unique targets down to a limiting magnitude of R = 24.4 yielded 156 secure redshifts. In addition to redshift information, we provide sky-subtracted one- and two-dimensional spectra of each target. Observations were conducted following the procedures of the Team Keck Redshift Survey (TKRS), thereby producing spectra that augment the TKRS sample while maintaining the uniformity of its spectral database.

We develop a new diagnostic method to classify galaxies into AGN hosts, star-forming galaxies, and absorption-dominated galaxies by combining the [O III]/Hbeta ratio with rest-frame U-B color. This can be used to robustly select AGNs in galaxy samples at intermediate redshifts (z<1). We compare the result of this optical AGN selection with X-ray selection using a sample of 3150 galaxies with 0.3<z<0.8 and I_AB<22, selected from the DEEP2 Galaxy Redshift Survey and the All-wavelength Extended Groth Strip International Survey (AEGIS). Among the 146 X-ray sources in this sample, 58% are classified optically as emission-line AGNs, the rest as star-forming galaxies or absorption-dominated galaxies. The latter are also known as "X-ray bright, optically normal galaxies" (XBONGs). Analysis of the relationship between optical emission lines and X-ray properties shows that the completeness of optical AGN selection suffers from dependence on the star formation rate and the quality of observed spectra. It also shows that XBONGs do not appear to be a physically distinct population from other X-ray detected, emission-line AGNs. On the other hand, X-ray AGN selection also has strong bias. About 2/3 of all emission-line AGNs at L_bol>10^44 erg/s in our sample are not detected in our 200 ks Chandra images, most likely due to moderate or heavy absorption by gas near the AGN. The 2--7 keV detection rate of Seyfert 2s at z~0.6 suggests that their column density distribution and Compton-thick fraction are similar to that of local Seyferts. Multiple sample selection techniques are needed to obtain as complete a sample as possible.

The broad band spectra of two Swift/BAT AGNs obtained from Suzaku follow-up observations are studied: NGC 612 and NGC 3081. Fitting with standard models, we find that both sources show similar spectra characterized by a heavy absorption with $N_{\rm{H}} \simeq 10^{24} \ \rm{cm}^{-2}$, the fraction of scattered light is $f_{\rm{scat}} = 0.5-0.8%$, and the solid angle of the reflection component is $\Omega/2\pi = 0.4-1.1$. To investigate the geometry of the torus, we apply numerical spectral models utilizing Monte Carlo simulations by Ikeda et al. (2009) to the Suzaku spectra. We find our data are well explained by this torus model, which has four geometrical parameters. The fit results suggest that NGC 612 has the torus half opening-angle of $\simeq 60^{\circ}-70^{\circ}$ and is observed from a nearly edge-on angle with a small amount of scattering gas, while NGC 3081 has a very small opening angle $\simeq 15^\circ$ and is observed on a face-on geometry, more like the deeply buried "new type" AGNs found by Ueda et al. (2007). We demonstrate the potential power of direct application of such numerical simulations to the high quality broad band spectra to unveil the inner structure of AGNs.

We quantify the fraction of galaxies at moderate redshifts (0.1<z<0.5) that appear red-and-dead in the optical, but in fact contain obscured star formation detectable in the infrared (IR), with the PRIsm MUlti-object Survey (PRIMUS). PRIMUS has measured ~120,000 robust redshifts with a precision of sigma_z/(1+z)~0.5% over 9.1 square degrees of the sky to the depth of i~23 (AB), up to redshift z~1. We specifically targeted 6.7 square degree fields with existing deep IR imaging from the Spitzer Space Telescope from the SWIRE and S-COSMOS surveys. We select in these fields an i band flux-limited sample (i<20 mag in the SWIRE fields and i<21 mag in the S-COSMOS field) of 3310 red-sequence galaxies at 0.1<z<0.5 for which we can reliably classify obscured star-forming and quiescent galaxies using IR color. Our sample constitutes the largest galaxy sample at intermediate redshift to study obscured star formation on the red sequence, and we present the first quantitative analysis of the fraction of obscured star-forming galaxies as a function of luminosity. We find that on average, at L ~ L*, about 15% of red-sequence galaxies have IR colors consistent with star-forming galaxies. The percentage of obscured star-forming galaxies increases by ~8% per mag with decreasing luminosity from the highest luminosities to L~0.2L*. Our results suggest that a significant fraction of red-sequence galaxies have ongoing star formation and that galaxy evolution studies based on optical color therefore need to account for this complication.

We present the results of spectroscopic, narrow-band and X-ray observations of a z=2.30 protocluster in the field of the QSO HS 1700+643. Using a sample of BX/MD galaxies, which are selected to be at z~2.2-2.7 by their rest-frame ultraviolet colours, we find that there are 5 protocluster AGN which have been identified by characteristic emission-lines in their optical/near-IR spectra; this represents an enhancement over the field significant at ~98.5 per cent confidence. Using a ~200 ks Chandra/ACIS-I observation of this field we detect a total of 161 X-ray point sources to a Poissonian false-probability limit of 4x10^{-6} and identify 8 of these with BX/MD galaxies. Two of these are spectroscopically confirmed protocluster members and are also classified as emission-line AGN. When compared to a similarly selected field sample the analysis indicates this is also evidence for an enhancement of X-ray selected BX/MD AGN over the field, significant at ~99 per cent confidence. Deep Lya narrow-band imaging reveals that a total of 4/123 Lya emitters (LAEs) are found to be associated with X-ray sources, with two of these confirmed protocluster members and one highly likely member. We do not find a significant enhancement of AGN activity in this LAE sample over that of the field (result significant at only 87 per cent confidence). The X-ray emitting AGN fractions for the BX/MD and LAE samples are found to be 6.9_{-4.4}^{+9.2} and 2.9_{-1.6}^{+2.9} per cent, respectively, for protocluster AGN with L_{2-10 keV}>4.6x10^{43} erg s^{-1} at z=2.30. These findings are similar to results from the z=3.09 protocluster in the SSA 22 field found by Lehmer et al. (2009), in that both suggest AGN activity is favoured in dense environments at z>2.

We present new observational determinations of the evolution of the 2-10keV X-ray luminosity function (XLF) of AGN. We utilise data from a number of surveys including both the 2Ms Chandra Deep Fields and the AEGIS-X 200ks survey, enabling accurate measurements of the evolution of the faint end of the XLF. We combine direct, hard X-ray selection and spectroscopic follow-up or photometric redshift estimates at z<1.2 with a rest-frame UV colour pre-selection approach at higher redshifts to avoid biases associated with catastrophic failure of the photometric redshifts. Only robust optical counterparts to X-ray sources are considered using a likelihood ratio matching technique. A Bayesian methodology is developed that considers redshift probability distributions, incorporates selection functions for our high redshift samples, and allows robust comparison of different evolutionary models. We find that the XLF retains the same shape at all redshifts, but undergoes strong luminosity evolution out to z~1, and an overall negative density evolution with increasing redshift, which thus dominates the evolution at earlier times. We do not find evidence that a Luminosity-Dependent Density Evolution, and the associated flattening of the faint-end slope, is required to describe the evolution of the XLF. We find significantly higher space densities of low-luminosity, high-redshift AGN than in prior studies, and a smaller shift in the peak of the number density to lower redshifts with decreasing luminosity. The total luminosity density of AGN peaks at z=1.2+/-0.1, but there is a mild decline to higher redshifts. We find >50% of black hole growth takes place at z>1, with around half in Lx<10^44 erg/s AGN.

One of the main themes in extragalactic astronomy for the next decade will be the evolution of galaxies over cosmic time. Many future observatories, including JWST, ALMA, GMT, TMT and E-ELT will intensively observe starlight over a broad redshift range, out to the dawn of the modern Universe when the first galaxies formed. It has, however, become clear that the properties and evolution of galaxies are intimately linked to the growth of their central black holes. Understanding the formation of galaxies, and their subsequent evolution, will therefore be incomplete without similarly intensive observations of the accretion light from supermassive black holes (SMBH) in galactic nuclei. To make further progress, we need to chart the formation of typical SMBH at z>6, and their subsequent growth over cosmic time, which is most effectively achieved with X-ray observations. Recent technological developments in X-ray optics and instrumentation now bring this within our grasp, enabling capabilities fully matched to those expected from flagship observatories at longer wavelengths.

We present the AEGIS-X survey, a series of deep Chandra ACIS-I observations of the Extended Groth Strip. The survey comprises pointings at 8 separate positions, each with nominal exposure 200ks, covering a total area of approximately 0.67 deg2 in a strip of length 2 degrees. We describe in detail an updated version of our data reduction and point source detection algorithms used to analyze these data. A total of 1325 band-merged sources have been found to a Poisson probability limit of 4e-6, with limiting fluxes of 5.3e-17 erg/cm2/s in the soft (0.5-2 keV) band and 3.8e-16 erg/cm2/s in the hard (2-10 keV) band. We present simulations verifying the validity of our source detection procedure and showing a very small, <1.5%, contamination rate from spurious sources. Optical/NIR counterparts have been identified from the DEEP2, CFHTLS, and Spitzer/IRAC surveys of the same region. Using a likelihood ratio method, we find optical counterparts for 76% of our sources, complete to R(AB)=24.1, and, of the 66% of the sources that have IRAC coverage, 94% have a counterpart to a limit of 0.9 microJy at 3.6 microns (m(AB)=23.8). After accounting for (small) positional offsets in the 8 Chandra fields, the astrometric accuracy of the Chandra positions is found to be 0.8 arcsec RMS, however this number depends both on the off-axis angle and the number of detected counts for a given source. All the data products described in this paper are made available via a public website.

We present a new method for determining the sensitivity of X-ray imaging observations, which correctly accounts for the observational biases that affect the probability of detecting a source of a given X-ray flux, without the need to perform a large number of time consuming simulations. We use this new technique to estimate the X-ray source counts in different spectral bands (0.5-2, 0.5-10, 2-10 and 5-10keV) by combining deep pencil-beam and shallow wide-area Chandra observations. The sample has a total of 6295 unique sources over an area of $\rm 11.8deg^2$ and is the largest used to date to determine the X-ray number counts. We determine, for the first time, the break flux in the 5-10 keV band, in the case of a double power-law source count distribution. We also find an upturn in the 0.5-2keV counts at fluxes below about 6e-17erg/s/cm2. We show that this can be explained by the emergence of normal star-forming galaxies which dominate the X-ray population at faint fluxes. The fraction of the diffuse X-ray background resolved into point sources at different spectral bands is also estimated. It is argued that a single population of Compton thick AGN cannot be responsible for the entire unresolved X-ray background in the energy range 2-10keV.

We combine Lyman-break colour selection with ultradeep (> 200 ks) Chandra X-ray imaging over a survey area of ~0.35 deg^2 to select high redshift AGN. Applying careful corrections for both the optical and X-ray selection functions, the data allow us to make the most accurate determination to date of the faint end of the X-ray luminosity function (XLF) at z~3. Our methodology recovers a number density of X-ray sources at this redshift which is at least as high as previous surveys, demonstrating that it is an effective way of selecting high z AGN. Comparing to results at z=1, we find no evidence that the faint slope of the XLF flattens at high z, but we do find significant (factor ~3.6) negative evolution of the space density of low luminosity AGN. Combining with bright end data from very wide surveys we also see marginal evidence for continued positive evolution of the characteristic break luminosity L*. Our data therefore support models of luminosity-dependent density evolution between z=1 and z=3. A sharp upturn in the the XLF is seen at the very lowest luminosities (Lx < 10^42.5 erg s^-1), most likely due to the contribution of pure X-ray starburst galaxies at very faint fluxes.