We present a census of the active black hole population at 1 < z < 2, by constructing the bivariate distribution function of black hole mass and Eddington ratio, employing a maximum likelihood ...fitting technique. The study of the active black hole mass function (BHMF) and the Eddington ratio distribution function (ERDF) allows us to clearly disentangle the active galactic nuclei (AGN) downsizing phenomenon, present in the AGN luminosity function, into its physical processes of black hole mass downsizing and accretion rate evolution. We are utilizing type-1 AGN samples from three optical surveys (VVDS, zCOSMOS and SDSS), that cover a wide range of 3 dex in luminosity over our redshift interval of interest. We investigate the cosmic evolution of the AGN population as a function of AGN luminosity, black hole mass and accretion rate. Compared to z = 0, we find a distinct change in the shape of the BHMF and the ERDF, consistent with downsizing in black hole mass. The active fraction or duty cycle of type-1 AGN at z ~ 1.5 is almost flat as a function of black hole mass, while it shows a strong decrease with increasing mass at z = 0. We are witnessing a phase of intense black hole growth, which is largely driven by the onset of AGN activity in massive SMBHs (supermassive black holes) towards z = 2. We finally compare our results to numerical simulations and semi-empirical models and while we find reasonable agreement over certain parameter ranges, we highlight the need to refine these models in order to match our observations.
Context.
Obscured active galactic nuclei (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 FIR and submillimeter wavelengths. Studying the physical properties of obscured AGN and their host galaxies is crucial to shedding light on the early stages of a massive system lifetime.
Aims.
We aim to investigate the contribution of the interstellar medium (ISM) to the obscuration of quasars in a sample of distant highly star forming galaxies and to unveil their morphological and kinematics properties.
Methods.
We exploit Atacama Large Millimeter/submillimeter Array Cycle 4 observations of the continuum (∼2.1 mm) and high-
J
CO emission of a sample of six X-ray selected, FIR detected galaxies hosting an obscured AGN at
z
spec
> 2.5 in the 7 Ms
Chandra
Deep Field-South. We measured the masses and sizes of the dust and molecular gas by fitting the images, visibilities, and spectra, 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.
Results.
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 from the root mean square of the images. We found that the detected galaxies are rich in gas and dust (molecular gas mass in the range < 0.5–2.7 × 10
10
M
⊙
for
α
CO
= 0.8 and up to ∼2 × 10
11
M
⊙
for
α
CO
= 6.5, and dust mass < 0.9–4.9 × 10
8
M
⊙
) and generally 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 × 10
10
M
⊙
.
Conclusions.
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 compare various star formation rate (SFR) indicators for star-forming galaxies at 1.4 < z < 2.5 in the COSMOS field. The main focus is on the SFRs from the far-IR (PACS-Herschel data) with those ...from the ultraviolet, for galaxies selected according to the BzK criterion. FIR-selected samples lead to a vastly different slope of the SFR–stellar mass (M
*) relation, compared to that of the dominant main-sequence population as measured from the UV, since the FIR selection picks predominantly only a minority of outliers. However, there is overall agreement between the main sequences derived with the two SFR indicators, when stacking on the PACS maps the BzK-selected galaxies. The resulting logarithmic slope of the SFR–M
* relation is ∼0.8–0.9, in agreement with that derived from the dust-corrected UV luminosity. Exploiting deeper 24 μm Spitzer data, we have characterized a subsample of galaxies with reddening and SFRs poorly constrained, as they are very faint in the B band. The combination of Herschel with Spitzer data has allowed us to largely break the age/reddening degeneracy for these intriguing sources, by distinguishing whether a galaxy is very red in B-z because of being heavily dust reddened, or whether because star formation has been (or is being) quenched. Finally, we have compared our SFR(UV) to the SFRs derived by stacking the radio data and to those derived from the Hα luminosity of a sample of star-forming galaxies at 1.4 < z < 1.7. The two sets of SFRs are broadly consistent as they are with the SFRs derived from the UV and by stacking the corresponding PACS data in various mass bins.
Using spectroscopic observations taken for the Visible Multi-Object Spectrograph (VIMOS) Ultra-Deep Survey (VUDS) we report here on the discovery of PCl J1001+0220, a massive proto-cluster of ...galaxies located at zspec ~ 4.57 in the COSMOS field. With nine spectroscopic members, the proto-cluster was initially detected as a ~12σ spectroscopic overdensity of typical star-forming galaxies in the blind spectroscopic survey of the early universe (2 < z ≲ 6) performed by VUDS. It was further mapped using a new technique developed which statistically combines spectroscopic and photometric redshifts, the latter derived from a recent compilation of incredibly deep multi-band imaging performed on the COSMOS field. Through various methods, the descendant mass of PCl J1001+0220 is estimated to be log (Mh/M⊙)z=0 $\log(\mathcal{M}_{h}/\mathcal{M}_{\odot})_{z=0}$log(Mh/M⊙)z=0 ~ 14.5–15 with a large amount of mass apparently already in place at z ~ 4.57. An exhaustive comparison was made between the properties of various spectroscopic and photometric member samples and matched samples of galaxies inhabiting less dense environments at the same redshifts. Tentative evidence is found for a fractional excess of older galaxies more massive in their stellar content amongst the member samples relative to the coeval field, an observation which suggests the pervasive early onset of vigorous star formation for proto-cluster galaxies. No evidence is found for the differences in the star formation rates (SFRs) of member and coeval field galaxies either through estimating by means of the rest-frame ultraviolet or through separately stacking extremely deep Very Large Array 3 GHz imaging for both samples. Additionally, no evidence for pervasive strong active galactic nuclei (AGN) activity is observed in either environment. Analysis of Hubble Space Telescope images of both sets of galaxies as well as their immediate surroundings provides weak evidence for an elevated incidence of galaxy–galaxy interaction within the bounds of the proto-cluster. The stacked and individual spectral properties of the two samples are compared, with a definite suppression of Lyα seen in the average member galaxy relative to the coeval field ( fesc, Lyα = 1.8−1.7+0.3 $f_{esc, \, \textrm{Ly}\alpha} = 1.8^{+0.3}_{-1.7}$fesc, Lyα=1.8−1.7+0.3% and 4.0−0.8+1.0 $4.0^{+1.0}_{-0.8}$4.0−0.8+1.0%, respectively). This observation along with other lines of evidence leads us to infer the possible presence of a large, cool, diffuse medium within the proto-cluster environment evocative of a nascent intracluster medium forming in the early universe.
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 gigayear 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 (≈5 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 Atacama large millimeter array (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 NH > 2 × 1024 cm−2 and NH > 6 × 1023 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.
Based on a sample of over 1800 radio AGN at redshifts out to z ~ 5, which have typical stellar masses within ~3 × (1010 − 1011)M⊙, and 3 GHz radio data in the COSMOS field, we derived the 1.4 GHz ...radio luminosity functions for radio AGN (L1.4 GHz~ 1022 − 1027 W Hz-1) out to z ~ 5. We constrained the evolution of this population via continuous models of pure density and pure luminosity evolutions, and we found best-fit parametrizations of Φ∗ ∝ (1 + z)(2.00 ± 0.18) − (0.60 ± 0.14)z, and L∗ ∝ (1 + z)(2.88 ± 0.82) − (0.84 ± 0.34)z, respectively, with a turnover in number and luminosity densities of the population at z ≈ 1.5. We converted 1.4 GHz luminosity to kinetic luminosity taking uncertainties of the scaling relation used into account. We thereby derived the cosmic evolution of the kinetic luminosity density provided by the AGN and compared this luminosity density to the radio-mode AGN feedback assumed in the Semi-Analytic Galaxy Evolution (SAGE) model, i.e., to the redshift evolution of the central supermassive black hole accretion luminosity taken in the model as the source of heating that offsets the energy losses of the cooling, hot halo gas, and thereby limits further stellar mass growth of massive galaxies. We find that the kinetic luminosity exerted by our radio AGN may be high enough to balance the radiative cooling of the hot gas at each cosmic epoch since z ~ 5. However, although our findings support the idea of radio-mode AGN feedback as a cosmologically relevant process in massive galaxy formation, many simplifications in both the observational and semi-analytic approaches still remain and need to be resolved before robust conclusions can be reached.
We study the relation of AGN accretion, star formation rate (SFR) and stellar mass (M
*) using a sample of ≈8600 star-forming galaxies up to z = 2.5 selected with Herschel imaging in the GOODS and ...COSMOS fields. For each of them we derive SFR and M
*, both corrected, when necessary, for emission from an active galactic nucleus (AGN), through the decomposition of their spectral energy distributions (SEDs). About 10 per cent of the sample are detected individually in Chandra observations of the fields. For the rest of the sample, we stack the X-ray maps to get average X-ray properties. After subtracting the X-ray luminosity expected from star formation and correcting for nuclear obscuration, we derive the average AGN accretion rate for both detected sources and stacks, as a function of M
*, SFR and redshift. The average accretion rate correlates with SFR and with M
*. The dependence on SFR becomes progressively more significant at z > 0.8. This may suggest that SFR is the original driver of these correlations. We find that average AGN accretion and star formation increase in a similar fashion with offset from the star-forming ‘main-sequence’. Our interpretation is that accretion on to the central black hole and star formation broadly trace each other, irrespective of whether the galaxy is evolving steadily on the main-sequence or bursting.
ABSTRACT We present the catalog of optical and infrared counterparts of the Chandra COSMOS-Legacy Survey, a 4.6 Ms Chandra program on the 2.2 deg2 of the COSMOS field, combination of 56 new ...overlapping observations obtained in Cycle 14 with the previous C-COSMOS survey. In this Paper we report the i, K, and 3.6 m identifications of the 2273 X-ray point sources detected in the new Cycle 14 observations. We use the likelihood ratio technique to derive the association of optical/infrared (IR) counterparts for 97% of the X-ray sources. We also update the information for the 1743 sources detected in C-COSMOS, using new K and 3.6 m information not available when the C-COSMOS analysis was performed. The final catalog contains 4016 X-ray sources, 97% of which have an optical/IR counterpart and a photometric redshift, while 54% of the sources have a spectroscopic redshift. The full catalog, including spectroscopic and photometric redshifts and optical and X-ray properties described here in detail, is available online. We study several X-ray to optical (X/O) properties: with our large statistics we put better constraints on the X/O flux ratio locus, finding a shift toward faint optical magnitudes in both soft and hard X-ray band. We confirm the existence of a correlation between X/O and the the 2-10 keV luminosity for Type 2 sources. We extend to low luminosities the analysis of the correlation between the fraction of obscured AGNs and the hard band luminosity, finding a different behavior between the optically and X-ray classified obscured fraction.
ABSTRACT
We present a study designed to measure the average Lyman-continuum escape fraction (〈fesc〉) of star-forming galaxies at z ≃ 3.5. We assemble a sample of 148 galaxies from the VANDELS ...spectroscopic survey at 3.35 ≤ zspec ≤ 3.95, selected to minimize line-of-sight contamination of their photometry. For this sample, we use ultra-deep, ground-based, U-band imaging and Hubble Space Telescope V-band imaging to robustly measure the distribution of $\mathcal {R_{\rm obs}}\, =(L_{\rm LyC}/L_{\rm UV})_{\rm obs}$. We then model the $\mathcal {R_{\rm obs}}$ distribution as a function of 〈fesc〉, carefully accounting for attenuation by dust, the intergalactic medium and the circumgalactic medium. A maximum likelihood fit to the $\mathcal {R_{\rm obs}}$ distribution returns a best-fitting value of $\langle f_{\rm esc}\rangle =0.07^{+0.02}_{-0.02}$, a result confirmed using an alternative Bayesian inference technique (both techniques exclude 〈fesc〉 = 0.0 at >3σ). By splitting our sample in two, we find evidence that 〈fesc〉 is positively correlated with Ly α equivalent width (Wλ(Ly α)), with high and low Wλ(Lyα) subsamples returning values of $\langle f_{\rm esc}\rangle =0.12^{+0.06}_{-0.04}$ and $\langle f_{\rm esc} \rangle =0.02^{+0.02}_{-0.01}$, respectively. In contrast, we find evidence that 〈fesc〉 is anticorrelated with intrinsic UV luminosity and UV dust attenuation; with low UV luminosity and dust attenuation subsamples both returning best fits in the range 0.10 ≤ 〈fesc〉 ≤ 0.22. We do not find a clear correlation between fesc and galaxy stellar mass, suggesting stellar mass is not a primary indicator of fesc. Although larger samples are needed to further explore these trends, our results suggest that it is entirely plausible that the low dust, low-metallicity galaxies found at z ≥ 6 will display the 〈fesc〉 ≥ 0.1 required to drive reionization.
We present the results from a ~ 500 ks Chandra observation of the z = 6.31 QSO SDSS J1030 + 0524. This is the deepest X-ray observation to date of a z ~ 6 QSO. The QSO is detected with a total of 125 ...net counts in the full (0.500A0–7 keV) band and its spectrum can be modeled by a single power-law model with photon index of Γ = 1.81 ± 0.18 and full band flux of f = 3.95 × 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 (ΔΓ ≈ −0.6) spectrum and a 2.5 times fainter flux. Such a variation, in a timespan of ~ 2 yr rest-frame, is unexpected for such a luminous QSO powered by a > 109M⨀ 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″ × 20″ southward of the QSO with a signal-to-noise ratio (S/N) of approximately six, hardness ratio of HR = 0.03+0.20−0.25 HR = 0.03 - 0.25 + 0.20 $\text{HR} = {0.03}_{-0.25}^{+0.20}$ , and soft band flux of f0.5– keV = 1.1+0.3−0.3 × 10−15 erg s−1 cm−2 f 0.5 - 2 keV = 1.1 - 0.3 + 0.3 × 10 - 15 erg s - 1 cm - 2 $ f_{0.5-2\,\text{keV}}={1.1}_{-0.3}^{+0.3}\enspace \times \enspace 1{0}^{-15\enspace }\mathrm{erg}\enspace {\mathrm{s}}^{-1}\;\rm{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 ≈ 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.