Recent studies show that a universal relation between black-hole (BH) growth and stellar mass ($M_\bigstar$) or star formation rate (SFR) is an oversimplification of BH-galaxy co-evolution, and that morphological and structural properties of host galaxies must also be considered. Particularly, a possible connection between BH growth and host-galaxy compactness was identified among star-forming (SF) galaxies. Utilizing $\approx 6300$ massive galaxies with $I_{\rm 814W}~<~24$ at $z$ $<$ 1.2 in the COSMOS field, we perform systematic partial-correlation analyses to investigate how sample-averaged BH accretion rate ($\rm \overline{BHAR}$) depends on host-galaxy compactness among SF galaxies, when controlling for morphology and $M_\bigstar$ (or SFR). The projected central surface-mass density within 1 kpc, $\Sigma_{1}$, is utilized to represent host-galaxy compactness in our study. We find that the $\rm \overline{BHAR}$-$\Sigma_{1}$ relation is stronger than either the $\rm \overline{BHAR}$-$M_\bigstar$ or $\rm \overline{BHAR}$-SFR relation among SF galaxies, and this $\rm \overline{BHAR}$-$\Sigma_{1}$ relation applies to both bulge-dominated galaxies and galaxies that are not dominated by bulges. This $\rm \overline{BHAR}$-$\Sigma_{1}$ relation among SF galaxies suggests a link between BH growth and the central gas density of host galaxies on the kpc scale, which may further imply a common origin of the gas in the vicinity of the BH and in the central $\sim$ kpc of the galaxy. This $\rm \overline{BHAR}$-$\Sigma_{1}$ relation can also be interpreted as the relation between BH growth and the central velocity dispersion of host galaxies at a given gas content, indicating the role of the host-galaxy potential well in feeding BHs.

Population synthesis models of actively accreting super-massive black holes (or active galactic nuclei -- AGN) predict a large fraction that must grow behind dense, obscuring screens of gas and dust. Deep X-ray surveys are thought to have provided the most complete and unbiased samples of AGN, but there is strong observational evidence that a portion of the population of obscured AGN is being missed. In this paper we use a sample of AGN derived from the deepest X-ray survey to date, the Chandra 7Ms GOODS-South Survey, to investigate the nature of low flux X-ray sources. We make full use of the extensive multi-wavelength coverage of the GOODS-South field, and cross-match our objects with wavelengths from the Radio to the IR. We find the low X-ray flux AGN in our sample have X-ray luminosities that indicate low-luminosity AGN classification, while their radio, infrared and optical counterparts indicate moderate to powerful AGN classification. We find the predicted column densities is on average an order of magnitude higher than the calculated column densities via X-ray detections for X-ray faint sources. We interpret our results as evidence of obscured AGN disguising as low-luminosity AGN via their X-ray luminosities. When we compare the estimation of the obscured AGN space density with and without these objects, we find a difference of 40% in the lowest X-ray luminosity regime probed by our sample.

We present the X-ray source catalog for the 479 ks Chandra exposure of the SDSS J1030+0524 field, that is centered on a region that shows the best evidence to date of an overdensity around a z > 6 quasar, and also includes a galaxy overdensity around a Compton-thick Fanaroff-Riley type II radio galaxy at z = 1.7. Using wavdetect for initial source detection and ACIS Extract for source photometry and significance assessment, we create preliminary catalogs of sources that are detected in the full, soft, and hard bands, respectively. We produce X-ray simulations that mirror our Chandra observation to filter our preliminary catalogs and get a completeness level of > 91% and a reliability level of 95% in each band. The catalogs in the three bands are then matched into a final main catalog of 256 unique sources. Among them, 244, 193, and 208 are detected in the full, soft, and hard bands, respectively. This makes J1030 field the fifth deepest extragalactic X-ray survey to date. The field is part of the Multiwavelength Survey by Yale-Chile (MUSYC), and is also covered by optical imaging data from the Large Binocular Camera (LBC) at the Large Binocular Telescope, near-IR imaging data from the Canada France Hawaii Telescope WIRCam, and Spitzer IRAC. Thanks to its dense multi-wavelength coverage, J1030 represents a legacy field for the study of large-scale structures around distant accreting supermassive black holes. Using a likelihood ratio analysis, we associate multi-band counterparts for 252 (98.4%) of the 256 Chandra sources, with an estimated reliability of 95%. Finally, we compute the cumulative number of sources in each X-ray band, finding that they are in general agreement with the results from the Chandra Deep Fields.

Obscured AGN represent a significant fraction of the entire AGN population, especially at high redshift (~70% at z=3--5). They are often characterized by the presence of large gas and dust reservoirs that are thought to sustain and possibly obscure vigorous star formation processes that make these objects shine at far-IR and sub-mm wavelengths. We exploit ALMA Cycle 4 observations of the continuum (~2.1mm) and high-J CO emission of a sample of six X-ray selected SMGs hosting an obscured AGN at z_spec>2.5 in the 7 Ms CDF-S. We measured the masses and sizes of the dust and molecular gas and we derived the gas density and column density on the basis of a uniform sphere geometry. Finally, we compared the measured column densities with those derived from the Chandra X-ray spectra. We detected both the continuum and line emission for three sources for which we measured both the flux density and size. For the undetected sources, we derived an upper limit on the flux density. We found that the detected galaxies are rich in gas and dust (molecular gas mass in the range <0.5 - 2.7 x 10^10 M_sun for {\alpha}_CO=0.8 and up to ~2 x 10^11~M_sun for {\alpha}_CO=6.5, and dust mass <0.9 - 4.9 x 10^8 M_sun) and compact (gas major axis 2.1-3.0 kpc, dust major axis 1.4-2.7 kpc). The column densities associated with the ISM are on the order of 10^(23-24) cm-2, which is comparable with those derived from the X-ray spectra. For the detected sources we also derived dynamical masses in the range 0.8 - 3.7 x 10^10 M_sun. We conclude that the ISM of high redshift galaxies can substantially contribute to nuclear obscuration up to the Compton-thick (>10^24 cm-2) regime. In addition, we found that all the detected sources show a velocity gradient reminding one rotating system, even though two of them show peculiar features in their morphology that can be associated with a chaotic, possibly merging, structure.

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.

We report the discovery of an extreme X-ray flux rise (by a factor of > 20) of the weak-line quasar SDSS J153913.47+395423.4 (hereafter SDSS J1539+3954) at z = 1.935. SDSS J1539+3954 is the most-luminous object among radio-quiet type 1 AGNs where such dramatic X-ray variability has been observed. Before the X-ray flux rise, SDSS J1539+3954 appeared X-ray weak compared with the expectation from its UV flux; after the rise, the ratio of its X-ray flux and UV flux is consistent with the majority of the AGN population. We also present a contemporaneous HET spectrum of SDSS J1539+3954, which demonstrates that its UV continuum level remains generally unchanged despite the dramatic increase in the X-ray flux, and its C iv emission line remains weak. The dramatic change only observed in the X-ray flux is consistent with a shielding model, where a thick inner accretion disk can block our line of sight to the central X-ray source. This thick inner accretion disk can also block the nuclear ionizing photons from reaching the high-ionization broad emission-line region, so that weak high-ionization emission lines are observed. Under this scenario, the extreme X-ray variability event may be caused by slight variations in the thickness of the disk. This event might also be explained by gravitational light-bending effects in a reflection model.

We report the discovery of a galaxy overdensity around a Compton-thick Fanaroff-Riley type II (FRII) radio galaxy at z=1.7 in the deep multiband survey around the z=6.3 QSO SDSS J1030+0524. Based on a 6hr VLT/MUSE and on a 4hr LBT/LUCI observation, we identify at least eight galaxy members in this structure with spectroscopic redshift z=1.687-1.699, including the FRII galaxy at z=1.699. Most of the identified overdensity members are blue, compact galaxies that are actively forming stars at rates of $\sim$8-60 $M_{\odot}$ yr$^{-1}$. Based on a 500ks Chandra ACIS-I observation we found that the FRII nucleus hosts a luminous QSO ($L_{2-10keV}=1.3\times10^{44}$ erg s$^{-1}$, intrinsic and rest-frame) that is obscured by Compton-thick absorption ($N_H=1.5\pm0.6\times 10^{24}$ cm$^{-2}$). Our Chandra observation, the deepest so far for a distant FRII within a galaxy overdensity, revealed significant diffuse X-ray emission within the region covered by the overdensity. In particular, X-ray emission extending for $\sim$240 kpc is found around the Eastern lobe of the FRII. Four out of the six MUSE star forming galaxies in the overdensity are distributed in an arc-like shape at the edge of this diffuse X-ray emission. The probability of observing by chance four out of the six $z=1.7$ sources at the edge of the diffuse emission is negligible. In addition, these four galaxies have the highest specific star formation rates of the MUSE galaxies in the overdensity and lie above the main sequence of field galaxies of equal stellar mass at z=1.7. We propose that the diffuse X-rays originate from an expanding bubble of gas that is shock-heated by the FRII jet, and that star formation is promoted by the compression of the cold interstellar medium of the galaxies around the bubble, which may be remarkable evidence of positive AGN feedback on cosmological scales. [shortened version]

Possible connections between central black-hole (BH) growth and host-galaxy compactness have been found observationally, which may provide insight into BH-galaxy coevolution: compact galaxies might have large amounts of gas in their centers due to their high mass-to-size ratios, and simulations predict that high central gas density can boost BH accretion. However, it is not yet clear if BH growth is fundamentally related to the compactness of the host galaxy, due to observational degeneracies between compactness, stellar mass ($M_\bigstar$), and star formation rate (SFR). To break these degeneracies, we carry out systematic partial-correlation studies to investigate the dependence of sample-averaged BH accretion rate ($\rm \overline{BHAR}$) on the compactness of host galaxies, represented by the surface-mass density, $\Sigma_\rm e$, or the projected central surface-mass density within 1 kpc, $\Sigma_1$. We utilize 8842 galaxies with H < 24.5 in the five CANDELS fields at z = 0.5-3. We find that $\rm \overline{BHAR}$ does not significantly depend on compactness when controlling for SFR or $M_\bigstar$ among bulge-dominated galaxies and galaxies that are not dominated by bulges, respectively. However, when testing is confined to star-forming galaxies at z = 0.5-1.5, we find that the $\rm \overline{BHAR}$-$\Sigma_1$ relation is not simply a secondary manifestation of a primary $\rm \overline{BHAR}$-$M_\bigstar$ relation, which may indicate a link between BH growth and the gas density within the central 1 kpc of galaxies.

X-ray emission from QSOs has been used to assess SMBH accretion properties up to $z$~6. However, at $z>6$ only ~15 QSOs are covered by sensitive X-ray observations, preventing a statistically significant investigation of the X-ray properties of QSOs in the first Gyr of the Universe. We present new Chandra observations of 10 $z>6$ QSOs, selected to have virial black-hole mass estimates from Mg II line spectroscopy. Adding archival X-ray data for an additional 15 $z>6$ QSOs, we investigate the X-ray properties of the QSO population in the first Gyr of the Universe, focusing in particular on the $L_{UV}-L_{X}$ relation, which is traced by the $\alpha_{ox}$ parameter, and the shape of their X-ray spectra. We performed photometric analyses to derive estimates of the X-ray luminosities, and thus the $\alpha_{ox}$ values and bolometric corrections ($K_{bol}=L_{bol}/L_{X}$). We compared the resulting $\alpha_{ox}$ and $K_{bol}$ distributions with the results found for QSO samples at lower redshift. Finally, we performed a basic X-ray spectral analysis of the brightest $z>6$ QSOs to derive their individual photon indices, and joint spectral analysis of the whole sample to estimate the average photon index. We confirm a lack of significant evolution of $\alpha_{ox}$ with redshift, extending the results from previous works up to $z>6$, and the trend of an increasing bolometric correction with increasing luminosity found for QSOs at lower redshifts. The average power-law photon index of our sample ($\Gamma=2.20_{-0.34}^{+0.39}$ and $\Gamma=2.13_{-0.13}^{+0.13}$ for sources with $<30$ and $>30$ net counts, respectively) is slightly steeper than, but still consistent with, typical QSOs at $z=1-6$. All these results point toward a lack of substantial evolution of the inner accretion-disk/hot-corona structure in QSOs from low redshift to $z>6$. Our data hint at generally high Eddington ratios at $z>6$.

While theoretical arguments predict that most of the early growth of supermassive black holes (SMBHs) happened during heavily obscured phases of accretion, current methods used for selecting $z>6$ quasars (QSOs) are strongly biased against obscured QSOs, thus considerably limiting our understanding of accreting SMBHs during the first Gyr of the Universe from an observational point of view. We report the $Chandra$ discovery of the first heavily obscured QSO candidate in the early universe, hosted by a close ($\approx5$ kpc) galaxy pair at $z=6.515$. One of the members is an optically classified type 1 QSO, PSO167-13. The companion galaxy was first detected as a [C II] emitter by ALMA. An X-ray source is significantly ($P=0.9996$) detected by $Chandra$ in the 2-5 keV band, with $<1.14$ net counts in the 0.5-2 keV band, although the current positional uncertainty does not allow a conclusive association with either PSO167-13 or its companion galaxy. From X-ray photometry and hardness-ratio arguments, we estimated an obscuring column density of $N_H>2\times10^{24}\,\mathrm{cm^{-2}}$ and $N_H>6\times10^{23}\,\mathrm{cm^{-2}}$ at $68\%$ and $90\%$ confidence levels, respectively. Thus, regardless of which of the two galaxies is associated with the X-ray emission, this source is the first heavily obscured QSO candidate at $z>6$.

Neutron star-neutron star mergers are known to be associated with short gamma-ray bursts. If the neutron star equation of state is sufficiently stiff, at least some of such mergers will leave behind a supramassive or even a stable neutron star that spins rapidly with a strong magnetic field (i.e., a magnetar). Such a magnetar signature may have been observed as the X-ray plateau following a good fraction (up to 50%) of short gamma-ray bursts, and it has been expected that one may observe short gamma-ray burst-less X-ray transients powered by double neutron star mergers. A fast X-ray transient (CDF-S XT1) was recently found to be associated with a faint host galaxy whose redshift is unknown. Its X-ray and host-galaxy properties allow several possibleexplanations including a short gamma-ray burst seen off axis, a low-luminosity gamma-ray burst at high redshift, or a tidal disruption event involving an intermediate mass black hole and a white dwarf. Here we report a second X-ray transient, CDF-S XT2, that is associated with a galaxy at redshift z = 0.738. The light curve is fully consistent with being powered by a millisecond magnetar. More intriguingly, CDF-S XT2 lies in the outskirts of its star-forming host galaxy with a moderate offset from the galaxy center, as short bursts often do. The estimated event rate density of similar X-ray transients, when corrected to the local value, is consistent with the double neutron star merger rate density inferred from the detection of GW170817.

We present a detailed X-ray spectral analysis of 1152 AGNs selected in the Chandra Deep Fields (CDFs), in order to identify highly obscured AGNs ($N_{\rm H} > 10^{23}\ \rm cm^{-2}$). By fitting spectra with physical models, 436 (38%) sources with $L_{\rm X} > 10^{42}\ \rm erg\ s^{-1}$ are confirmed to be highly obscured, including 102 Compton-thick (CT) candidates. We propose a new hardness-ratio measure of the obscuration level which can be used to select highly obscured AGN candidates. The completeness and accuracy of applying this method to our AGNs are 88% and 80%, respectively. The observed logN-logS relation favors cosmic X-ray background models that predict moderate (i.e., between optimistic and pessimistic) CT number counts. 19% (6/31) of our highly obscured AGNs that have optical classifications are labeled as broad-line AGNs, suggesting that, at least for part of the AGN population, the heavy X-ray obscuration is largely a line-of-sight effect, i.e., some high-column-density clouds on various scales (but not necessarily a dust-enshrouded torus) along our sightline may obscure the compact X-ray emitter. After correcting for several observational biases, we obtain the intrinsic NH distribution and its evolution. The CT-to-highly-obscured fraction is roughly 52% and is consistent with no evident redshift evolution. We also perform long-term (~17 years in the observed frame) variability analyses for 31 sources with the largest number of counts available. Among them, 17 sources show flux variabilities: 31% (5/17) are caused by the change of NH, 53% (9/17) are caused by the intrinsic luminosity variability, 6% (1/17) are driven by both effects, and 2 are not classified due to large spectral fitting errors.

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.

Massive 10^6-10^10 Msun black holes (BHs) are ubiquitous in local galactic nuclei. They were common by the time the Universe is several Gyr old, and many of them were in place within the first 1~Gyr after the Big Bang. Their quick assembly has been attributed to mechanisms such as the rapid collapse of gas into the nuclei of early protogalaxies, accretion and mergers of stellar-mass BHs accompanying structure formation at early times, and the runaway collapse of early, ultra-dense stellar clusters. The origin of the early massive BHs remains an intriguing and long-standing unsolved puzzle in astrophysics. Here we discuss strategies for discerning between BH seeding models using electromagnetic observations. We argue that the most direct answers will be obtained through detection of BHs with masses M<10^5 Msun at redshifts z>10, where we expect them to first form. Reaching out to these redshifts and down to these masses is crucial, because BHs are expected to lose the memory of their initial assembly by the time they grow well above 10^5 Msun and are incorporated into higher-mass galaxies. The best way to detect 10^4-10^5 Msun BHs at high redshifts is by a sensitive X-ray survey. Critical constraining power is augmented by establishing the properties and the environments of their host galaxies in deep optical/IR imaging surveys. Required OIR data can be obtained with the JWST and WFIRST missions. The required X-ray flux limits (down to 10^{-19} erg/s/cm^2) are accessible only with a next-generation X-ray observatory which has both high (sub-1") angular resolution and high throughput. A combination of deep X-ray and OIR surveys will be capable of probing several generic markers of the BH seed scenarios, and resolving the long-stanging puzzle of their origin. These electromagnetic observations are also highly synergistic with the information from LISA on high-z BH mergers.

Observations in the local universe show a tight correlation between the masses of supermassive black holes (SMBHs; $M_{\rm BH}$) and host-galaxy bulges ($M_{\rm bulge}$), suggesting a strong connection between SMBH and bulge growth. However, direct evidence for such a connection in the distant universe remains elusive. We have studied sample-averaged SMBH accretion rate ($\overline{\rm BHAR}$) for bulge-dominated galaxies at $z=0.5-3$. While previous observations found $\overline{\rm BHAR}$ is strongly related to host-galaxy stellar mass ($M_\star$) for the overall galaxy population, our analyses show that, for the bulge-dominated population, $\overline{\rm BHAR}$ is mainly related to SFR rather than $M_\star$. This ${\overline{\rm BHAR}}$-SFR relation is highly significant, e.g. $9.0\sigma$ (Pearson statistic) at $z=0.5-1.5$. Such a $\overline{\rm BHAR}$-SFR connection does not exist among our comparison sample of galaxies that are not bulge-dominated, for which $M_\star$ appears to be the main determinant of SMBH accretion. This difference between the bulge-dominated and comparison samples indicates that SMBHs only coevolve with bulges rather than the entire galaxies, explaining the tightness of the local $M_{\rm BH}$-$M_{\rm bulge}$ correlation. Our best-fit $\overline{\rm BHAR}$-SFR relation for the bulge-dominated sample is $\log\overline{\rm BHAR} = \log\mathrm{SFR} - (2.48\pm0.05)$ (solar units). The best-fit $\overline{\rm BHAR}/\mathrm{SFR}$ ratio ($10^{-2.48}$) for bulge-dominated galaxies is similar to the observed $M_{\rm BH}/M_{\rm bulge}$ values in the local universe. Our results reveal that SMBH and bulge growth are in lockstep, and thus non-causal scenarios of merger averaging are unlikely the origin of the $M_{\rm BH}$-$M_{\rm bulge}$ correlation. This lockstep growth also predicts that the $M_{\rm BH}$-$M_{\rm bulge}$ relation should not have strong redshift dependence.

We present a multi-wavelength study of seven AGN at spectroscopic redshift >2.5 in the 7 Ms Chandra Deep Field South, selected to have good FIR/sub-mm detections. Our aim is to investigate the possibility that the obscuration observed in the X-rays can be produced by the interstellar medium (ISM) of the host galaxy. Based on the 7 Ms Chandra spectra, we measured obscuring column densities N$_{H, X}$ in excess of 7x10$^{22}$ cm$^{-2}$ and intrinsic X-ray luminosities L$_{X}$>10$^{44}$ erg s$^{-1}$ for our targets, as well as equivalent widths for the Fe K$\alpha$ emission line EW>0.5-1 keV. We built the UV-to-FIR spectral energy distributions by using broad-band photometry from CANDELS and Herschel catalogs. By means of an SED decomposition technique, we derived stellar masses (M$_{*}$~10$^{11}$ Msun), IR luminosities (L$_{IR}$>10$^{12}$ Lsun), star formation rates (SFR~190-1680 Msun yr$^{-1}$) and AGN bolometric luminosities (L$_{bol}$~10$^{46}$ erg s$^{-1}$) for our sample. We used an empirically-calibrated relation between gas masses and FIR/sub-mm luminosities and derived M$_{gas}$~0.8-5.4x10$^{10}$ Msun. High-resolution (0.3-0.7'') ALMA data (when available, CANDELS data otherwise) were used to estimate the galaxy size and hence the volume enclosing most of the ISM under simple geometrical assumptions. These measurements were then combined to derive the column density associated with the ISM of the host, on the order of N$_{H, ISM}$~10$^{23-24}$ cm$^{-2}$. The comparison between the ISM column densities and those measured from the X-ray spectral analysis shows that they are similar. This suggests that, at least at high redshift, significant absorption on kpc scales by the dense ISM in the host likely adds to or substitutes that produced by circumnuclear gas on pc scales (i.e., the torus of unified models). The lack of unobscured AGN among our ISM-rich targets supports this scenario.

We report the discovery of disappearance of Mg ii, Al iii, C iv, and Si iv broad absorption lines (BALs) at the same velocity, accompanied by a new Civ BAL emerging at a higher velocity, in the quasar J0827+4252 at z = 2.038. This is the first report of BAL disappearance (i) over Mg ii, Al iii, C iv, and Si iv ions and (ii) in a weak emission-line quasar (WLQ). The discovery is based on four spectra from the SDSS and one follow-up spectrum from HET/LRS2. The simultaneous Civ BAL disappearance and emergence at different velocities, together with no variations in the CRTS light curve, indicate that ionization changes in the absorbing material are unlikely to cause the observed BAL variability. Our analyses reveal that transverse motion is the most likely dominant driver of the BAL disappearance/emergence. Given the presence of mildly relativistic BAL outflows and an apparently large C iv emission-line blueshift that is likely associated with strong bulk outflows in this WLQ, J0827+4252 provides a notable opportunity to study extreme quasar winds and their potential in expelling material from inner to large-scale regions.

This white paper specifies the footprints, cadence requirements, and total-depth requirements needed to allow the most-successful AGN studies in the four currently selected LSST Deep-Drilling Fields (DDFs): ELAIS-S1, XMM-LSS, CDF-S, and COSMOS. The information provided on cadence and total-depth requirements will also likely be applicable to enabling effective AGN science in any additional DDFs that are chosen.

In deep X-ray surveys, active galactic nuclei (AGNs) with a broad range of luminosities have been identified. However, cosmologically distant low-luminosity AGN (LLAGN, $L_{\mathrm{X}} \lesssim 10^{42}$ erg s$^{-1}$) identification still poses a challenge due to significant contamination from host galaxies. Based on the 7 Ms Chandra Deep Field-South (CDF-S) survey, the longest timescale ($\sim 17$ years) deep X-ray survey to date, we utilize an X-ray variability selection technique to search for LLAGNs that remain unidentified among the CDF-S X-ray sources. We find 13 variable sources from 110 unclassified CDF-S X-ray sources. Except for one source which could be an ultraluminous X-ray source, the variability of the remaining 12 sources is most likely due to accreting supermassive black holes. These 12 AGN candidates have low intrinsic X-ray luminosities, with a median value of $7 \times10^{40}$ erg s$^{-1}$. They are generally not heavily obscured, with an average effective power-law photon index of 1.8. The fraction of variable AGNs in the CDF-S is independent of X-ray luminosity and is only restricted by the total number of observed net counts, confirming previous findings that X-ray variability is a near-ubiquitous property of AGNs over a wide range of luminosities. There is an anti-correlation between X-ray luminosity and variability amplitude for high-luminosity AGNs, but as the luminosity drops to $\lesssim 10^{42}$ erg s$^{-1}$, the variability amplitude no longer appears dependent on the luminosity. The entire observed luminosity-variability trend can be roughly reproduced by an empirical AGN variability model based on a broken power-law power spectral density function.

It is well known that environment affects galaxy evolution, which is broadly related to supermassive black hole (SMBH) growth. We investigate whether SMBH evolution also depends on host-galaxy local (sub-Mpc) and global ($\approx 1-10$ Mpc) environment. We construct the surface-density field (local environment) and cosmic web (global environment) in the COSMOS field at $z=0.3-3.0$. The environments in COSMOS range from the field to clusters ($M_\mathrm{halo} \lesssim 10^{14}\ M_\odot$), covering the environments where ${\approx 99\%}$ of galaxies in the Universe reside. We measure sample-averaged SMBH accretion rate ($\overline{\mathrm{BHAR}}$) from X-ray observations, and study its dependence on overdensity and cosmic-web environment at different redshifts while controlling for galaxy stellar mass ($M_\star$). Our results show that $\overline{\mathrm{BHAR}}$ does not significantly depend on overdensity or cosmic-web environment once $M_\star$ is controlled, indicating that environment-related physical mechanisms (e.g. tidal interaction and ram-pressure stripping) might not significantly affect SMBH growth. We find that $\overline{\mathrm{BHAR}}$ is strongly related to host-galaxy $M_\star$, regardless of environment.

Several works have studied the relation between X-ray, UV, and wind properties in broad absorption line quasars (BALQSOs), generally concluding that the formation of strong winds is tightly connected with the suppression of the ionizing EUV/X-ray emission. The Eddington ratio ($\lambda_{Edd}$), which measures the accretion rate, is also known to be related with outflow and emission-line properties in the general quasar population. Moreover, models describing quasar accretion depend on $\lambda_{Edd}$, which can thus possibly affect the relative production of accelerating UV and ionizing EUV/X-ray radiation. In this work, for the first time, we investigated whether BALQSO X-ray properties are related with the Eddington ratio. We selected a sample of 30 BALQSOs with accurate measurements of black-hole mass and BAL properties from the literature, and we complemented it with 4 additional BALQSOs we observed with \xmm\, to populate the low and high Eddington-ratio regimes. We did not find evidence for a strong relation between $\lambda_{Edd}$ and X-ray suppression, which however shows a significant correlation with the strength of the UV absorption features. These findings are confirmed also by considering a sample of mini-BALQSOs collected from the literature.

We present an X-ray point-source catalog from the XMM-Large Scale Structure survey region (XMM-LSS), one of the XMM-Spitzer Extragalactic Representative Volume Survey (XMM-SERVS) fields. We target the XMM-LSS region with $1.3$ Ms of new XMM-Newton AO-15 observations, transforming the archival X-ray coverage in this region into a 5.3 deg$^2$ contiguous field with uniform X-ray coverage totaling $2.7$ Ms of flare-filtered exposure, with a $46$ ks median PN exposure time. We provide an X-ray catalog of 5242 sources detected in the soft (0.5-2 keV), hard (2-10 keV), and/or full (0.5-10 keV) bands with a 1% expected spurious fraction determined from simulations. A total of 2381 new X-ray sources are detected compared to previous source catalogs in the same area. Our survey has flux limits of $1.7\times10^{-15}$, $1.3\times10^{-14}$, and $6.5\times10^{-15}$ erg cm$^{-2}$ s$^{-1}$ over 90% of its area in the soft, hard, and full bands, respectively, which is comparable to those of the XMM-COSMOS survey. We identify multiwavelength counterpart candidates for 99.9% of the X-ray sources, of which 93% are considered as reliable based on their matching likelihood ratios. The reliabilities of these high-likelihood-ratio counterparts are further confirmed to be $\approx 97$% reliable based on deep Chandra coverage over $\approx 5$% of the XMM-LSS region. Results of multiwavelength identifications are also included in the source catalog, along with basic optical-to-infrared photometry and spectroscopic redshifts from publicly available surveys. We compute photometric redshifts for X-ray sources in 4.5 deg$^2$ of our field where forced-aperture multi-band photometry is available; $>70$% of the X-ray sources in this subfield have either spectroscopic or high-quality photometric redshifts.

We present the results from a $\sim500$ ks Chandra observation of the $z=6.31$ QSO SDSS J1030+0524. This is the deepest X-ray observation to date of a $z\sim6$ QSO. The QSO is detected with a total of 125 net counts in the full ($0.5-7$ keV) band and its spectrum can be modeled by a single power-law model with photon index of $\Gamma = 1.81 \pm 0.18$ and full band flux of $f=3.95\times 10^{-15}$ erg s$^{-1}$ cm$^{-2}$. When compared with the data obtained by XMM-Newton in 2003, our Chandra observation in 2017 shows a harder ($\Delta \Gamma \approx -0.6$) spectrum and a 2.5 times fainter flux. Such a variation, in a timespan of $\sim2$ yrs rest-frame, is unexpected for such a luminous QSO powered by a $> 10^9 \: M_{\odot}$ black hole. The observed source hardening and weakening could be related to an intrinsic variation in the accretion rate. However, the limited photon statistics does not allow us to discriminate between an intrinsic luminosity and spectral change, and an absorption event produced by an intervening gas cloud along the line of sight. We also report the discovery of diffuse X-ray emission that extends for 30"x20" southward the QSO with a signal-to-noise ratio of $\sim$6, hardness ratio of $HR=0.03_{-0.25}^{+0.20}$, and soft band flux of $f_{0.5-2 \: keV}= 1.1_{-0.3}^{+0.3} \times 10^{-15}$ erg s$^{-1}$ cm$^{-2}$, that is not associated to a group or cluster of galaxies. We discuss two possible explanations for the extended emission, which may be either associated with the radio lobe of a nearby, foreground radio galaxy (at $z \approx 1-2$), or ascribed to the feedback from the QSO itself acting on its surrounding environment, as proposed by simulations of early black hole formation.

We imaged with ALMA and ARGOS/LUCI the molecular gas and the dust and stellar continuum in XID2028, an obscured QSO at z=1.593, where the presence of a massive outflow in the ionized gas component traced by the [O III]5007 emission has been resolved up to 10 kpc. This target represents a unique test case to study QSO 'feedback in action' at the peak epoch of AGN-galaxy coevolution. The QSO has been detected in the CO(5-4) transition and in the 1.3mm continuum, at ~30 and ~20 {\sigma} significance respectively, with both emissions confined in the central (<4 kpc) radius area. Our analysis suggests the presence of a fast rotating molecular disc (v~400 km/s) on very compact scales, and well inside the galaxy extent seen in the rest-frame optical light (~10 kpc, as inferred from the LUCI data). Adding available measurements in additional two CO transitions, CO(2-1) and CO(3-2), we could derive a total gas mass of ~10$^{10}$ M$_\odot$, thanks to a critical assessment of CO excitation and the comparison with Rayleigh-Jeans continuum estimate. This translates into a very low gas fraction (<5%) and depletion time scales of 40-75 Myr, reinforcing the result of atypical gas consumption conditions in XID2028, possibly due to feedback effects on the host galaxy. Finally, we also detect at ~5{\sigma} the presence of high velocity CO gas, which we interpret as a signature of galaxy-scale molecular outflow, spatially coincident with the ionised gas outflow. XID2028 represents therefore a unique case where the measurement of total outflowing mass (~500-800 M$_\odot$/yr) including the molecular and atomic components, in both the ionised and neutral phases, has been attempted for a high-z QSO.

Hot Dust-Obscured Galaxies (Hot DOGs) are hyperluminous ($L_{\mathrm{8-1000\,\mu m}}>10^{13}\,\mathrm{L_\odot}$) infrared galaxies with extremely high (up to hundreds of K) dust temperatures. The sources powering both their extremely high luminosities and dust temperatures are thought to be deeply buried and rapidly accreting supermassive black holes (SMBHs). Hot DOGs could therefore represent a key evolutionary phase in which the SMBH growth peaks. X-ray observations can be used to study their obscuration levels and luminosities. In this work, we present the X-ray properties of the 20 most-luminous ($L_{\mathrm{bol}}\gtrsim10^{14}\, L_\odot$) known Hot DOGs at $z=2-4.6$. Five of them are covered by long-exposure ($10-70$ ks) Chandra and XMM-Newton observations, with three being X-ray detected, and we study their individual properties. One of these sources (W0116$-$0505) is a Compton-thick candidate, with column density $N_H=(1.0-1.5)\times10^{24}\,\mathrm{cm^{-2}}$ derived from X-ray spectral fitting. The remaining 15 Hot DOGs have been targeted by a Chandra snapshot (3.1 ks) survey. None of these 15 is individually detected; therefore we applied a stacking analysis to investigate their average emission. From hardness-ratio analysis, we constrained the average obscuring column density and intrinsic luminosity to be log$N_H\,\mathrm{[cm^{-2}]}>23.5$ and $L_X\gtrsim10^{44}\,\mathrm{erg\,cm^{-2}\,s^{-1}}$, which are consistent with results for individually detected sources. We also investigated the $L_X-L_{6\mu\mathrm{m}}$ and $L_X-L_{bol}$ relations, finding hints that Hot DOGs are typically X-ray weaker than expected, although larger samples of luminous obscured QSOs are needed to derive solid conclusions.

Previous studies suggest that the growth of supermassive black holes (SMBHs) may be fundamentally related to host-galaxy stellar mass ($M_\star$). To investigate this SMBH growth-$M_\star$ relation in detail, we calculate long-term SMBH accretion rate as a function of $M_\star$ and redshift [$\overline{\rm BHAR}(M_\star, z)$] over ranges of $\log(M_\star/M_\odot)=\text{9.5--12}$ and $z=\text{0.4--4}$. Our $\overline{\rm BHAR}(M_\star, z)$ is constrained by high-quality survey data (GOODS-South, GOODS-North, and COSMOS), and by the stellar mass function and the X-ray luminosity function. At a given $M_\star$, $\overline{\rm BHAR}$ is higher at high redshift. This redshift dependence is stronger in more massive systems (for $\log(M_\star/M_\odot)\approx 11.5$, $\overline{\rm BHAR}$ is three decades higher at $z=4$ than at $z=0.5$), possibly due to AGN feedback. Our results indicate that the ratio between $\overline{\rm BHAR}$ and average star formation rate ($\overline{\rm SFR}$) rises toward high $M_\star$ at a given redshift. This $\overline{\rm BHAR}/\overline{\rm SFR}$ dependence on $M_\star$ does not support the scenario that SMBH and galaxy growth are in lockstep. We calculate SMBH mass history [$M_{\rm BH}(z)$] based on our $\overline{\rm BHAR}(M_\star, z)$ and the $M_\star(z)$ from the literature, and find that the $M_{\rm BH}$-$M_\star$ relation has weak redshift evolution since $z\approx 2$. The $M_{\rm BH}/M_\star$ ratio is higher toward massive galaxies: it rises from $\approx 1/5000$ at $\log M_\star\lesssim 10.5$ to $\approx 1/500$ at $\log M_\star \gtrsim 11.2$. Our predicted $M_{\rm BH}/M_\star$ ratio at high $M_\star$ is similar to that observed in local giant ellipticals, suggesting that SMBH growth from mergers is unlikely to dominate over growth from accretion.

We systematically analyze X-ray variability of active galactic nuclei (AGNs) in the 7~Ms \textit{Chandra} Deep Field-South survey. On the longest timescale ($\approx~17$ years), we find only weak (if any) dependence of X-ray variability amplitudes on energy bands or obscuration. We use four different power spectral density (PSD) models to fit the anti-correlation between normalized excess variance ($\sigma^2_{\rm nxv}$) and luminosity, and obtain a best-fit power law index $\beta=1.16^{+0.05}_{-0.05}$ for the low-frequency part of AGN PSD. We also divide the whole light curves into 4 epochs in order to inspect the dependence of $\sigma^2_{\rm nxv}$ on these timescales, finding an overall increasing trend. The analysis of these shorter light curves also infers a $\beta$ of $\sim 1.3$ that is consistent with the above-derived $\beta$, which is larger than the frequently-assumed value of $\beta=1$. We then investigate the evolution of $\sigma^2_{\rm nxv}$. No definitive conclusion is reached due to limited source statistics but, if present, the observed trend goes in the direction of decreasing AGN variability at fixed luminosity toward large redshifts. We also search for transient events and find 6 notable candidate events with our considered criteria. Two of them may be a new type of fast transient events, one of which is reported here for the first time. We therefore estimate a rate of fast outbursts $\langle\dot{N}\rangle = 1.0^{+1.1}_{-0.7}\times 10^{-3}~\rm galaxy^{-1}~yr^{-1}$ and a tidal disruption event~(TDE) rate $\langle\dot{N}_{\rm TDE}\rangle=8.6^{+8.5}_{-4.9}\times 10^{-5}~\rm galaxy^{-1}~yr^{-1}$ assuming the other four long outbursts to be TDEs.

We investigate the population of high-redshift ($3\leq z < 6$) AGN selected in the two deepest X-ray surveys, the 7 Ms \textit{Chandra} Deep Field-South and 2 Ms \textit{Chandra} Deep Field-North. Their outstanding sensitivity and spectral characterization of faint sources allow us to focus on the sub-$L_*$ regime (log$L_{\mathrm{X}}\lesssim44$), poorly sampled by previous works using shallower data, and the obscured population. Taking fully into account the individual photometric-redshift probability distribution functions, the final sample consists of $\approx102$ X-ray selected AGN at $3\leq z < 6$. The fraction of AGN obscured by column densities log$N_{\mathrm{H}}>23$ is $\sim0.6-0.8$, once incompleteness effects are taken into account, with no strong dependence on redshift or luminosity. We derived the high-redshift AGN number counts down to $F_{\mathrm{0.5-2\,keV}}=7\times10^{-18}\,\mathrm{erg\,cm^{-2}\,s^{-1}}$, extending previous results to fainter fluxes, especially at $z>4$. We put the tightest constraints to date on the low-luminosity end of AGN luminosity function at high redshift. The space-density, in particular, declines at $z>3$ at all luminosities, with only a marginally steeper slope for low-luminosity AGN. By comparing the evolution of the AGN and galaxy densities, we suggest that such a decline at high luminosities is mainly driven by the underlying galaxy population, while at low luminosities there are hints of an intrinsic evolution of the parameters driving nuclear activity. Also, the black-hole accretion rate density and star-formation rate density, which are usually found to evolve similarly at $z\lesssim3$, appear to diverge at higher redshifts.

We study the X-ray variability properties of distant AGNs in the Chandra Deep Field-South region over 17 years, up to $z\sim 4$, and compare them with those predicted by models based on local samples. We use the results of Monte Carlo simulations to account for the biases introduced by the discontinuous sampling and the low-count regime. We confirm that variability is an ubiquitous property of AGNs, with no clear dependence on the density of the environment. The variability properties of high-z AGNs, over different temporal timescales, are most consistent with a Power Spectral Density (PSD) described by a broken (or bending) power-law, similar to nearby AGNs. We confirm the presence of an anti-correlation between luminosity and variability, resulting from the dependence of variability on BH mass and accretion rate. We explore different models, finding that our acceptable solutions predict that BH mass influences the value of the PSD break frequency, while the Eddington ratio $\lambda_{Edd}$ affects the PSD break frequency and, possibly, the PSD amplitude as well. We derive the evolution of the average $\lambda_{Edd}$ as a function of redshift, finding results in agreement with measurements based on different estimators. The large statistical uncertainties make our results consistent with a constant Eddington ratio, although one of our models suggest a possible increase of $\lambda_{Edd}$ with lookback time up to $z\sim 2-3$. We conclude that variability is a viable mean to trace the accretion history of supermassive BHs, whose usefulness will increase with future, wide-field/large effective area X-ray missions.

Many cosmological studies predict that early supermassive black holes (SMBHs) can only form in the most massive dark matter halos embedded within large scale structures marked by galaxy over-densities that may extend up to 10 physical Mpc. This scenario, however, has not been confirmed observationally, as the search for galaxy over-densities around high-z quasars has returned conflicting results. The field around the z=6.28 quasar SDSSJ1030+0524 (J1030) is unique for multi-band coverage and represents an excellent data legacy for studying the environment around a primordial SMBH. In this paper we present wide-area (25x25 arcmin), Y- and J-band imaging of the J1030 field obtained with the near infrared camera WIRCam at the Canada-France-Hawaii Telescope (CFHT). We built source catalogues in the Y- and J-band, and matched those with our photometric catalogue in the r, z, i bands presented in Morselli et al. (2014). We used these new infrared data together with H and K and Spitzer/IRAC data to refine our selection of Lyman Break Galaxies (LBGs), extending our selection criteria to galaxies in the range 25.2<zAB<25.7. We selected 21 robust high-z candidates in the J1030 field with photometric redshift around 6 and colors i-z>=1.3. We found a significant asymmetry in the distribution of the high-z galaxies in J1030, supporting the existence of a coherent large-scale structure around the quasar. We compared our results with those of Bowler et al. (2015), who adopted similar LBGs selection criteria, and estimated an over-density of galaxies in the field of delta = 2.4, which is significant at >4 sigma. The over-density value and its significance are higher than those found in Morselli et al. (2014), and we interpret this as evidence of an improved LBG selection.

We investigate the dependence of black-hole accretion rate (BHAR) on host-galaxy star formation rate (SFR) and stellar mass ($M_*$) in the CANDELS/GOODS-South field in the redshift range of $0.5\leq z < 2.0$. Our sample consists of $\approx 18000$ galaxies, allowing us to probe galaxies with $0.1 \lesssim \mathrm{SFR} \lesssim 100\ M_\odot\ \mathrm{yr}^{-1}$ and/or $10^8 \lesssim M_* \lesssim 10^{11}\ M_{\odot}$. We use sample-mean BHAR to approximate long-term average BHAR. Our sample-mean BHARs are derived from the $Chandra$ Deep Field-South 7 Ms observations, while the SFRs and $M_*$ have been estimated by the CANDELS team through SED fitting. The average BHAR is correlated positively with both SFR and $M_*$, and the BHAR-SFR and BHAR-$M_*$ relations can both be described acceptably by linear models with a slope of unity. However, BHAR appears to be correlated more strongly with $M_*$ than SFR. This result indicates that $M_*$ is the primary host-galaxy property related to black-hole growth, and the apparent BHAR-SFR relation is largely a secondary effect due to the star-forming main sequence. Among our sources, massive galaxies ($M_* \gtrsim 10^{10} M_{\odot}$) have significantly higher BHAR/SFR ratios than less-massive galaxies, indicating the former have higher black-hole fueling efficiency and/or higher SMBH occupation fraction than the latter. Our results can naturally explain the observed proportionality between $M_{\rm BH}$ and $M_*$ for local giant ellipticals, and suggest their $M_{\rm BH}/M_*$ is higher than that of local star-forming galaxies. Among local star-forming galaxies, massive systems might have higher $M_{\rm BH}/M_*$ compared to dwarfs.

The coeval AGN and galaxy evolution and the observed local relations between SMBHs and galaxy properties suggest some connection or feedback between SMBH growth and galaxy build-up. We looked for correlations between properties of X-ray detected AGN and their FIR detected host galaxies, to find quantitative evidences for this connection, highly debated in the latest years. We exploit the rich multi-wavelength data set available in the COSMOS field for a large sample (692 sources) of AGN and their hosts, in the redshift range $0.1<z<4$. We use X-ray data to select AGN and determine their properties (intrinsic luminosity and nuclear obscuration), and broad-band SED fitting to derive host galaxy properties (stellar mass $M_*$ and star formation rate SFR). We find that the AGN 2-10 keV luminosity ($L_{\rm X}$) and the host $8-1000~\mu m$ star formation luminosity ($L_{\rm IR}^{\rm SF}$) are significantly correlated. However, the average host $L_{\rm IR}^{\rm SF}$ has a flat distribution in bins of AGN $L_{\rm X}$, while the average AGN $L_{\rm X}$ increases in bins of host $L_{\rm IR}^{\rm SF}$, with logarithmic slope of $\sim0.7$, in the redshifts range $0.4<z<1.2$. We also discuss the comparison between the distribution of these two quantities and the predictions from hydro-dynamical simulations. Finally we find that the average column density ($N_H$) shows a positive correlation with the host $M_*$, at all redshifts, but not with the SFR (or $L_{\rm IR}^{\rm SF}$). This translates into a negative correlation with specific SFR. Our results are in agreement with the idea that BH accretion and SF rates are correlated, but occur with different variability time scales. The presence of a positive correlation between $N_H$ and host $M_*$ suggests that the X-ray $N_H$ is not entirely due to the circum-nuclear obscuring torus, but may also include a contribution from the host galaxy.

We present X-ray source catalogs for the $\approx7$ Ms exposure of the Chandra Deep Field-South (CDF-S), which covers a total area of 484.2 arcmin$^2$. Utilizing WAVDETECT for initial source detection and ACIS Extract for photometric extraction and significance assessment, we create a main source catalog containing 1008 sources that are detected in up to three X-ray bands: 0.5-7.0 keV, 0.5-2.0 keV, and 2-7 keV. A supplementary source catalog is also provided including 47 lower-significance sources that have bright ($K_s\le23$) near-infrared counterparts. We identify multiwavelength counterparts for 992 (98.4%) of the main-catalog sources, and we collect redshifts for 986 of these sources, including 653 spectroscopic redshifts and 333 photometric redshifts. Based on the X-ray and multiwavelength properties, we identify 711 active galactic nuclei (AGNs) from the main-catalog sources. Compared to the previous $\approx4$ Ms CDF-S catalogs, 291 of the main-catalog sources are new detections. We have achieved unprecedented X-ray sensitivity with average flux limits over the central $\approx1$ arcmin$^2$ region of $\approx1.9\times10^{-17}$, $6.4\times10^{-18}$, and $2.7\times10^{-17}$ erg cm$^{-2}$ s$^{-1}$ in the three X-ray bands, respectively. We provide cumulative number-count measurements observing, for the first time, that normal galaxies start to dominate the X-ray source population at the faintest 0.5-2.0 keV flux levels. The highest X-ray source density reaches $\approx50\,500$ deg$^{-2}$, and $47\%\pm4\%$ of these sources are AGNs ($\approx23\,900$ deg$^{-2}$).

We briefly review how X-ray observations of high-redshift active galactic nuclei (AGNs) at z = 4-7 have played a critical role in understanding their basic demographics as well as their physical processes; e.g., absorption by nuclear material and winds, accretion rates, and jet emission. We point out some key remaining areas of uncertainty, highlighting where further Chandra and XMM-Newton observations/analyses, combined with new multiwavelength survey data, can advance understanding over the next decade.

We perform long-term ($\approx 15$ yr, observed-frame) X-ray variability analyses of the 68 brightest radio-quiet active galactic nuclei (AGNs) in the 6 Ms $Chandra$ Deep Field-South (CDF-S) survey; the majority are in the redshift range of $0.6-3.1$, providing access to penetrating rest-frame X-rays up to $\approx 10-30$ keV. Twenty-four of the 68 sources are optical spectral type I AGNs, and the rest (44) are type II AGNs. The time scales probed in this work are among the longest for X-ray variability studies of distant AGNs. Photometric analyses reveal widespread photon-flux variability: $90\%$ of AGNs are variable above a 95% confidence level, including many X-ray obscured AGNs and several optically classified type II quasars. We characterize the intrinsic X-ray luminosity ($L_{\rm{X}}$) and absorption ($N_{\rm{H}}$) variability via spectral fitting. Most (74%) sources show $L_{\rm{X}}$ variability; the variability amplitudes are generally smaller for quasars. A Compton-thick candidate AGN shows variability of its high-energy X-ray flux, indicating the size of reflecting material to be $\lesssim 0.3$ pc. $L_{\rm{X}}$ variability is also detected in a broad absorption line (BAL) quasar. The $N_{\rm{H}}$ variability amplitude for our sample appears to rise as time separation increases. About 16% of sources show $N_{\rm{H}}$ variability. One source transitions from an X-ray unobscured to obscured state while its optical classification remains type I; this behavior indicates the X-ray eclipsing material is not large enough to obscure the whole broad-line region.

We exploit the 7 Ms \textit{Chandra} observations in the \chandra\,Deep Field-South (\mbox{CDF-S}), the deepest X-ray survey to date, coupled with CANDELS/GOODS-S data, to measure the total X-ray emission arising from 2076 galaxies at $3.5\leq z < 6.5$. This aim is achieved by stacking the \textit{Chandra} data at the positions of optically selected galaxies, reaching effective exposure times of $\geq10^9\mathrm{s}$. We detect significant ($>3.7\sigma$) X-ray emission from massive galaxies at $z\approx4$. We also report the detection of massive galaxies at $z\approx5$ at a $99.7\%$ confidence level ($2.7\sigma$), the highest significance ever obtained for X-ray emission from galaxies at such high redshifts. No significant signal is detected from galaxies at even higher redshifts. The stacking results place constraints on the BHAD associated with the known high-redshift galaxy samples, as well as on the SFRD at high redshift, assuming a range of prescriptions for X-ray emission due to X- ray binaries. We find that the X-ray emission from our sample is likely dominated by processes related to star formation. Our results show that low-rate mass accretion onto SMBHs in individually X-ray-undetected galaxies is negligible, compared with the BHAD measured for samples of X-ray detected AGN, for cosmic SMBH mass assembly at high redshift. We also place, for the first time, constraints on the faint-end of the AGN X-ray luminosity function ($\mathrm{logL_X\sim42}$) at $z>4$, with evidence for fairly flat slopes. The implications of all of these findings are discussed in the context of the evolution of the AGN population at high redshift.

Improving the capabilities of detecting faint X-ray sources is fundamental to increase the statistics on faint high-z AGN and star-forming galaxies. We performed a simultaneous Maximum Likelihood PSF fit in the [0.5-2] keV and [2-7] keV energy bands of the 4 Ms{\em Chandra} Deep Field South (CDFS) data at the position of the 34930 CANDELS H-band selected galaxies. For each detected source we provide X-ray photometry and optical counterpart validation. We validated this technique by means of a raytracing simulation. We detected a total of 698 X-ray point-sources with a likelihood $\mathcal{L}$$>$4.98 (i.e. $>$2.7$\sigma$). We show that the prior knowledge of a deep sample of Optical-NIR galaxies leads to a significant increase of the detection of faint (i.e. $\sim$10$^{-17}$ cgs in the [0.5-2] keV band) sources with respect to "blind" X-ray detections. By including previous X-ray catalogs, this work increases the total number of X-ray sources detected in the 4 Ms CDFS, CANDELS area to 793, which represents the largest sample of extremely faint X-ray sources assembled to date. Our results suggest that a large fraction of the optical counterparts of our X-ray sources determined by likelihood ratio actually coincides with the priors used for the source detection. Most of the new detected sources are likely star-forming galaxies or faint absorbed AGN. We identified a few sources sources with putative photometric redshift z$>$4. Despite the low number statistics and the uncertainties on the photo-z, this sample significantly increases the number of X--ray selected candidate high-z AGN.

Star-forming galaxies are considered to be the leading candidate sources that dominate the cosmic reionization at z>7, and the search for analogs at moderate redshift showing Lyman continuum (LyC) leakage is currently a active line of research. We have observed a star-forming galaxy at z=3.2 with Hubble/WFC3 in the F336W filter, corresponding to the 730-890A rest-frame, and detect LyC emission. This galaxy is very compact and also has large Oxygen ratio [OIII]5007/[OII]3727 (>=10). No nuclear activity is revealed from optical/near-infrared spectroscopy and deep multi-band photometry (including the 6Ms X-ray, Chandra). The measured escape fraction of ionizing radiation spans the range 50-100\%, depending on the IGM attenuation. The LyC emission is detected at S/N=10 with m(F336W)=27.57+/-0.11 and it is spatially unresolved, with effective radius R_e<200pc. Predictions from photoionization and radiative transfer models are in line with the properties reported here, indicating that stellar winds and supernova explosions in a nucleated star-forming region can blow cavities generating density-bounded conditions compatible with optically thin media. Irrespective to the nature of the ionizing radiation, spectral signatures of these sources over the entire electromagnetic spectrum are of central importance for their identification during the epoch of reionization, when the LyC is unobservable. Intriguingly, the Spitzer/IRAC photometric signature of intense rest-frame optical emissions ([OIII]+Hbeta) observed recently at z~7.5-8.5 is similar to what is observed in this galaxy. Only the James Webb Space Telescope will measure optical line ratios at z>7 allowing a direct comparison with lower redshift LyC emitters, as reported here.

Gas outflows are believed to play a pivotal role in shaping galaxies, as they regulate both star formation and black hole growth. Despite their ubiquitous presence, the origin and the acceleration mechanism of such powerful and extended winds is not yet understood. Direct observations of the cold gas component in objects with detected outflows at other wavelengths are needed to assess the impact of the outflow on the host galaxy interstellar medium (ISM). We observed with the Plateau de Bure Interferometer an obscured quasar at z~1.5, XID2028, for which the presence of an ionised outflow has been unambiguously signalled by NIR spectroscopy. The detection of CO(3-2) emission in this source allows us to infer the molecular gas content and compare it to the ISM mass derived from the dust emission. We then analyze the results in the context of recent insights on scaling relations, which describe the gas content of the overall population of star-forming galaxies at a similar redshifts. The Star formation efficiency (~100) and gas mass (M_gas=2.1-9.5x10^{10} M_sun) inferred from the CO(3-2) line depend on the underlying assumptions on the excitation of the transition and the CO-to-H2 conversion factor. However, the combination of this information and the ISM mass estimated from the dust mass suggests that the ISM/gas content of XID2028 is significantly lower than expected for its observed M$_\star$, sSFR and redshift, based on the most up-to-date calibrations (with gas fraction <20% and depletion time scale <340 Myr). Overall, the constraints we obtain from the far infrared and millimeter data suggest that we are observing QSO feedback able to remove the gas from the host

We present results from the deep XMM-Newton observations of the two brightest X-ray sources in the Chandra Deep Field South (CDFS), PID 203 (z=0.544) and PID 319 (z=0.742). The long exposure of 2.5 Ms over a 10 year period (net 4 yr with a 6 yr gap) makes it possible to obtain high quality X-ray spectra of these two Type I AGN with X-ray luminosity of 10^44 erg/s, typical luminosity for low-redshift PG quasars, track their X-ray variability both in flux and spectral shape. Both sources showed X-ray flux variability of ~10-20 per cent in rms which is similar in the soft (0.5-2 keV) and hard (2-7 keV) bands. PID 203, which has evidence for optical extinction, shows modest amount of absorption (nH~1e21cm^-2) in the X-ray spectrum. Fe K emission is strongly detected in both objects with EW~0.2 keV. The lines in both objects are moderately broad and exhibit marginal evidence for variability in shape and flux, indicating that the bulk of the line emission come from their accretion disks rather than distant tori.

We present the hard-band ($2-10\,\mathrm{keV}$) X-ray luminosity function (HXLF) of $0.5-2\,\mathrm{keV}$ band selected AGN at high redshift. We have assembled a sample of 141 AGN at $3<z\lesssim5$ from X-ray surveys of different size and depth, in order to sample different regions in the $ L_X - z$ plane. The HXLF is fitted in the range $\mathrm{logL_X\sim43-45}$ with standard analytical evolutionary models through a maximum likelihood procedure. The evolution of the HXLF is well described by a pure density evolution, with the AGN space density declining by a factor of $\sim10$ from $z=3$ to 5. A luminosity-dependent density evolution model which, normally, best represents the HXLF evolution at lower redshift, is also consistent with the data, but a larger sample of low-luminosity ($\mathrm{logL_X}<44$), high-redshift AGN is necessary to constrain this model. We also estimated the intrinsic fraction of AGN obscured by a column density $\mathrm{logN_H}\geq23$ to be $0.54\pm0.05$, with no strong dependence on luminosity. This fraction is higher than the value in the Local Universe, suggesting an evolution of the luminous ($\mathrm{L_X>10^{44}\mathrm{erg\,s^{-1}}}$) obscured AGN fraction from $z=0$ to $z>3$.

We have investigated the gas content of a sample of several hundred AGN host galaxies at z$<$1 and compared it with a sample of inactive galaxies, matched in bins of stellar mass and redshift. Gas masses have been inferred from the dust masses, obtained by stacked Herschel far-IR and sub-mm data in the GOODS and COSMOS fields, under reasonable assumptions and metallicity scaling relations for the dust-to-gas ratio. We find that AGNs are on average hosted in galaxies much more gas rich than inactive galaxies. In the vast majority of stellar mass bins, the average gas content of AGN hosts is higher than in inactive galaxies. The difference is up to a factor of ten higher in low stellar mass galaxies, with a significance of 6.5$\sigma$. In almost half of the AGN sample the gas content is three times higher than in the control sample of inactive galaxies. Our result strongly suggests that the probability of having an AGN activated is simply driven by the amount of gas in the host galaxy; this can be explained in simple terms of statistical probability of having a gas cloud falling into the gravitational potential of the black hole. The increased probability of an AGN being hosted by a star-forming galaxy, identified by previous works, may be a consequence of the relationship between gas content and AGN activity, found in this paper, combined with the Schmidt-Kennicutt law for star formation.

We report ALMA Cycle 0 observations at 1.3mm of LESS J033229.4-275619 (XID403), an Ultraluminous Infrared Galaxy at $z=4.75$ in the Chandra Deep Field South hosting a Compton-thick QSO. The source is not resolved in our data at a resolution of $\sim$0.75 arcsec, placing an upper-limit of 2.5 kpc to the half-light radius of the continuum emission from heated-dust. After deconvolving for the beam size, however, we found a $\sim3\sigma$ indication of an intrinsic source size of $0.27\pm0.08$ arcsec (Gaussian FWHM), which would correspond to $r_{half}\sim0.9\pm0.3$ kpc. We build the far-IR SED of XID403 by combining datapoints from both ALMA and Herschel and fit it with a modified blackbody spectrum. For the first time, we measure the dust temperature $T_d=58.5\pm5.3$ K in this system, which is comparable to what has been observed in other high-z submillimeter galaxies. The measured star formation rate is SFR=$1020\pm150$ $M_{\odot}$ yr$^{-1}$, in agreement with previous estimates at lower S/N. Based on the measured SFR and source size, we constrain the SFR surface density to be $\Sigma_{SFR}>26\;M_{\odot}$yr$^{-1}$kpc$^{-2}$ ($\sim200\;M_{\odot}$yr$^{-1}$kpc$^{-2}$ for $r_{half}\sim0.9$ kpc). The compactness of this starburst is comparable to what has been observed in other local and high-z starburst galaxies. If the gas mass measured from previous [CII] and CO(2-1) observations at low resolution is confined within the same dust region, assuming $r_{half}\sim0.9\pm0.3$ kpc, this would produce a column density of $N_H\sim0.3-1.1\times10^{24}$cm$^{-2}$ towards the central SMBH, similar to the column density of $\approx1.4\times10^{24}$cm$^{-2}$ measured from the X-rays. Then, in principle, if both gas and dust were confined on sub-kpc scales, this would be sufficient to produce the observed X-ray column density without any need of a pc-scale absorber [abridged].

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.

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.

X-ray spectroscopy of active galactic nuclei (AGN) offers the opportunity to directly probe the inner regions of the accretion disk. We present the results of our analysis of average AGN XMM-Newton X-ray spectra in the Chandra Deep Field South observation (XMM CDFS). We computed the average spectrum of a sample of 54 AGN with spectroscopic redshifts and signal-to-noise ratio S/N > 15 in the 2-12 keV rest-frame band in at least one EPIC camera. We have taken the effects of combining spectra from sources at different redshifts and both EPIC cameras into account, as well as their spectral resolution; we checked our results using thorough simulations. We explored the Fe line components of distant AGN focusing on the narrow core which arises from regions far from the central engine and on the putative relativistic component (from the accretion disk). The average spectrum shows a highly significant Fe feature. Our model-independent estimates of the equivalent width (EW) suggest a higher EW in a broader range. The line, modelled as an unresolved Gaussian, is significant at 6.8 sigma and has an EW=95 eV (full sample). The current data can be fitted equally well adding a relativistic profile to the narrow component (in the full sample, EW=140 eV and 67 eV respectively for the relativistic and narrow lines). Thanks to the high quality of the XMM CDFS spectra and to the detailed modelling of the continuum and instrumental effects, we have shown that the most distant AGN exhibit a highly significant Fe emission feature. It can be modelled both with narrow and broad lines. We found tantalising evidence for reflection by material both very close and far away from the central engine. The EW of both features are similar to those observed in individual nearby AGN, hence they must be a widespread characteristic of AGN, since otherwise the average values would be smaller than observed.

We present results from a spectral analysis of a sample of high-redshift (z>3) X-ray selected AGN in the 4 Ms Chandra Deep Field South (CDF-S), the deepest X-ray survey to date. The sample is selected using the most recent spectroscopic and photometric information available in this field. It consists of 34 sources with median redshift z=3.7, 80 median net counts in the 0.5-7 keV band and median rest-frame absorption-corrected luminosity $L_{2-10 \rmn{keV}}\approx1.5\times10^{44}\rmn{erg} \rmn{s^{-1}}$. Spectral analysis for the full sample is presented and the intrinsic column density distribution, corrected for observational biases using spectral simulations, is compared with the expectations of X-ray background (XRB) synthesis models. We find that $\approx57$ per cent of the sources are highly obscured ($N_H>10^{23}\rmn{cm^{-2}}$). Source number counts in the $0.5-2\rmn{keV}$ band down to flux $F_{0.5-2 \rmn{keV}}\approx4\times10^{-17}\rmn{erg} \rmn{s^{-1}cm^{-2}}$ are also presented. Our results are consistent with a decline of the AGN space density at z>3 and suggest that, at those redshifts, the AGN obscured fraction is in agreement with the expectations of XRB synthesis models.

We present results on a search of heavily obscured active galaxies z>1.7 using the rest-frame 9-20 keV excess for X-ray sources detected in the deep XMM-CDFS survey. Out of 176 sources selected with the conservative detection criteria (>8 sigma) in the first source catalogue of Ranalli et al., 46 objects lie in the redshift range of interest with the median redshift z~2.5. Their typical rest-frame 10-20 keV luminosity is 1e+44 erg/s, as observed. Among optically faint objects that lack spectroscopic redshift, four were found to be strongly absorbed X-ray sources, and the enhanced Fe K emission or absorption features in their X-ray spectra were used to obtain X-ray spectroscopic redshifts. Using the X-ray colour-colour diagram based on the rest-frame 3-5 keV, 5-9 keV, and 9-20 keV bands, seven objects were selected for their 9-20 keV excess and were found to be strongly absorbed X-ray sources with column density of nH > 0.6e+24 cm-2, including two possible Compton thick sources. While they are emitting at quasar luminosity, ~3/4 of the sample objects are found to be absorbed by nH > 1e+22 cm-2. A comparison with local AGN at the matched luminosity suggests an increasing trend of the absorbed source fraction for high-luminosity AGN towards high redshifts.