SPIDERS (The SPectroscopic IDentification of eROSITA Sources) is a large spectroscopic programme for X-ray selected galaxy clusters as part of the Sloan Digital Sky Survey-IV (SDSS-IV). We describe the final dataset in the context of SDSS Data Release 16 (DR16): the survey overall characteristics, final targeting strategies, achieved completeness and spectral quality, with special emphasis on its use as a galaxy cluster sample for cosmology applications. SPIDERS now consists of about 27,000 new optical spectra of galaxies selected within 4,000 photometric red sequences, each associated with an X-ray source. The excellent spectrograph efficiency and a robust analysis pipeline yield a spectroscopic redshift measurement success rate exceeding 98%, with a median velocity accuracy of 20 km s$^{-1}$ (at $z=0.2$). Using the catalogue of 2,740 X-ray galaxy clusters confirmed with DR16 spectroscopy, we reveal the three-dimensional map of the galaxy cluster distribution in the observable Universe up to $z\sim0.6$. We highlight the homogeneity of the member galaxy spectra among distinct regions of the galaxy cluster phase space. Aided by accurate spectroscopic redshifts and by a model of the sample selection effects, we compute the galaxy cluster X-ray luminosity function and we present its lack of evolution up to $z=0.6$. Finally we discuss the prospects of forthcoming large multiplexed spectroscopic programmes dedicated to follow up the next generation of all-sky X-ray source catalogues.

This white paper describes the science case for Very Long Baseline Interferometry (VLBI) and provides suggestions towards upgrade paths for the European VLBI Network (EVN). The EVN is a distributed long-baseline radio interferometric array, that operates at the very forefront of astronomical research. Recent results, together with the new science possibilities outlined in this vision document, demonstrate the EVN's potential to generate new and exciting results that will transform our view of the cosmos. Together with e-MERLIN, the EVN provides a range of baseline lengths that permit unique studies of faint radio sources to be made over a wide range of spatial scales. The science cases are reviewed in six chapters that cover the following broad areas: cosmology, galaxy formation and evolution, innermost regions of active galactic nuclei, explosive phenomena and transients, stars and stellar masers in the Milky Way, celestial reference frames and space applications. The document concludes with identifying the synergies with other radio, as well as multi-band/multi-messenger instruments, and provide the recommendations for future improvements. The appendices briefly describe other radio VLBI arrays, the technological framework for EVN developments, and a selection of spectral lines of astrophysical interest below 100 GHz. The document includes a glossary for non-specialists, and a list of acronyms at the end.

We study the optical linear and circular polarisation in the optically-thin regime of the core and jet of M87. Observations were acquired two days before the Event Horizon Telescope (EHT) campaign in early April 2017. A high degree ($\sim 20$ per cent) of linear polarisation (P$_{\rm lin}$) is detected in the bright jet knots resolved at $\sim 10\, \rm{arcsec}$ to $23\, \rm{arcsec}$ ($0.8$-$1.8\, \rm{kpc}$) from the centre, whereas the nucleus and inner jet show P$_{\rm lin} \lesssim 5$ per cent. The position angle of the linear polarisation shifts by $\sim 90$ degrees from each knot to the adjacent ones, with the core angle perpendicular to the first knot. The nucleus was in a low level of activity (P$_{\rm lin} \sim 2$-$3$ per cent), and no emission was detected from HST-1. No circular polarisation was detected either in the nucleus or the jet above a $3\sigma$ level of P$_{\rm circ} \leq 1.5$ per cent, discarding the conversion of P$_{\rm lin}$ into P$_{\rm circ}$. A disordered magnetic field configuration or a mix of unresolved knots polarised along axes with different orientations could explain the low P$_{\rm lin}$. The latter implies a smaller size of the core knots, in line with current interferometric observations. Polarimetry with EHT can probe this scenario in the future. A steep increase of both P$_{\rm lin}$ and P$_{\rm circ}$ with increasing frequency is expected for the optically-thin domain, above the turnover point. This work describes the methodology to recover the four Stokes parameters using a $\lambda/4$ wave-plate polarimeter.

This paper documents the sixteenth data release (DR16) from the Sloan Digital Sky Surveys; the fourth and penultimate from the fourth phase (SDSS-IV). This is the first release of data from the southern hemisphere survey of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2); new data from APOGEE-2 North are also included. DR16 is also notable as the final data release for the main cosmological program of the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), and all raw and reduced spectra from that project are released here. DR16 also includes all the data from the Time Domain Spectroscopic Survey (TDSS) and new data from the SPectroscopic IDentification of ERosita Survey (SPIDERS) programs, both of which were co-observed on eBOSS plates. DR16 has no new data from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey (or the MaNGA Stellar Library "MaStar"). We also preview future SDSS-V operations (due to start in 2020), and summarize plans for the final SDSS-IV data release (DR17).

Through the years numerous attempts have been made to connect the phenomenology and physics of mass accretion onto stellar-mass and super-massive black holes in a scale-invariant fashion. In this paper, we explore this connection at the radiatively-efficient (and non-jetted) end of accretion modes by comparing the relationship between the luminosity of the accretion disk and corona in the two source classes. We analyse 458 RXTE-PCA archival observations of the X-ray binary (XRB) GX339-4 focusing on the soft and soft-intermediate states, which have been suggested to be analogous to radiatively efficient AGN. The observed scatter in the $\log L_{disk}-\log L_{corona}$ relationship of GX339-4 is high ($\sim0.43\,$dex) and significantly larger than in a representative sample of radiatively-efficient, non- or weakly-jetted AGN ($\sim0.30\,$dex). On the face of it, this would appear contrary to the hypothesis that the systems simply scale with mass. On the other hand we also find that GX339-4 and our AGN sample show different $\dot{m}$ and $\Gamma$ distributions, with the latter being broader in GX339-4 (dispersion of $\sim0.16$ cf. $\sim0.08$ for AGN). GX339-4 also shows an overall softer slope, with mean $\sim2.20$ as opposed to $\sim2.07$ for the AGN sample. Remarkably, once similarly broad $\Gamma$ and $\dot{m}$ distributions are selected, the AGN sample overlaps nicely with GX339-4 observations in the mass-normalised $\log L_{disk}-\log L_{corona}$ plane, with a scatter of $\sim0.30-0.33\,$dex. This indicates that a mass-scaling of properties might hold after all, with our results being consistent with the disk-corona systems in AGN and XRBs exhibiting the same physical processes, albeit under different conditions for instance in terms of temperature, optical depth and/or electron energy distribution in the corona, heating-cooling balance, coronal geometry and/or black hole spin.

Context. Several new multi-object spectrographs are currently planned or under construction that are capable of observing thousands of Galactic and extragalactic objects simultaneously. Aims. In this paper we present a probabilistic fibre-to-target assignment algorithm that takes spectrograph targeting constraints into account and is capable of dealing with multiple concurrent surveys. We present this algorithm using the 4-metre Multi-Object Spectroscopic Telescope (4MOST) as an example. Methods. The key idea of the proposed algorithm is to assign probabilities to fibre-target pairs. The assignment of probabilities takes the fibre positioner's capabilities and constraints into account. Additionally, these probabilities include requirements from surveys and take the required exposure time, number density variation, and angular clustering of targets across each survey into account. The main advantage of a probabilistic approach is that it allows for accurate and easy computation of the target selection function for the different surveys, which involves determining the probability of observing a target, given an input catalogue. Results. The probabilistic fibre-to-target assignment allows us to achieve maximally uniform completeness within a single field of view. The proposed algorithm maximises the fraction of successfully observed targets whilst minimising the selection bias as a function of exposure time. In the case of several concurrent surveys, the algorithm maximally satisfies the scientific requirements of each survey and no specific survey is penalised or prioritised. Conclusions. The algorithm presented is a proposed solution for the 4MOST project that allows for an unbiased targeting of many simultaneous surveys. With some modifications, the algorithm may also be applied to other multi-object spectroscopic surveys.

We look to provide a detailed description of the SPectroscopic IDentification of ERosita Sources (SPIDERS) survey, an SDSS-IV programme aimed at obtaining spectroscopic classification and redshift measurements for complete samples of sufficiently bright X-ray sources. We describe the SPIDERS X-ray Point Source Spectroscopic Catalogue, considering its store of 11,092 observed spectra drawn from a parent sample of 14,759 ROSAT and XMM sources over an area of 5,129 deg$^2$ covered in SDSS-IV by the eBOSS survey. This programme represents the largest systematic spectroscopic observation of an X-ray selected sample. A total of 10,970 (98.9\%) of the observed objects are classified and 10,849 (97.8\%) have secure redshifts. The majority of the spectra (10,070 objects) are active galactic nuclei (AGN), 522 are cluster galaxies, and 294 are stars. The observed AGN redshift distribution is in good agreement with simulations based on empirical models for AGN activation and duty cycle. Forming composite spectra of type 1 AGN as a function of the mass and accretion rate of their black holes reveals systematic differences in the H-beta emission line profiles. This study paves the way for systematic spectroscopic observations of sources that are potentially to be discovered in the upcoming eROSITA survey over a large section of the sky.

We present a statistical analysis of the optical properties of an X-ray selected Type 1 AGN sample, using high signal to-noise ratio ($S/N>20$) spectra of the counterparts of the ROSAT/2RXS sources in the footprint of the SDSS-IV/SPIDERS (Spectroscopic IDentification of eROSITA Sources) programme. The final sample contains 2100 sources. It significantly extends the redshift and luminosity ranges ($\rm z \sim 0.01 - 0.80$ and $\rm L_{0.1-2.4 \,keV}$ ~ $\rm 2.0 \times 10^{41}$ - $\rm 1.0 \times 10^{46} \, erg s^{-1}$ ) used so far in this kind of analysis. By means of a Principal Component Analysis, we derive Eigenvector (EV) 1 and 2 in an eleven dimensional optical and X-ray parameter space, which are consistent with previous results. The validity of the correlations of the Eddington ratio $\rm L/L_{Edd}$ with EV1 and the black hole mass with EV2 are strongly confirmed These results imply that $\rm L/L_{Edd}$ and black hole mass are related to the diversity of the optical properties of Type 1 AGN. Investigating the relation of the width and asymmetry of $\rm H\beta$ and the relative strength of the iron emission $\rm r_{FeII}$ , we show that our analysis supports the presence of a distinct kinematic region: the Very Broad Line Region. Furthermore, comparing sources with a red-asymmetric broad $\rm H\beta$ emission line to sources for which it is blue-asymmetric, we find an intriguing difference in the correlation of the FeII and the continuum emission strengths. We show that this contrasting behaviour is consistent with a flattened, stratified model of the Broad Line Region, in which the FeII emitting region is shielded from the central source.

We present the brightest cluster galaxies (BCGs) catalog for SPectroscoic IDentification of eROSITA Sources (SPIDERS) DR14 cluster program value-added catalog. We list the 416 BCGs identified as part of this process, along with their stellar mass, star formation rates, and morphological properties. We identified the BCGs based on the available spectroscopic data from SPIDERS and photometric data from SDSS. We computed stellar masses and SFRs of the BCGs on the basis of SDSS, WISE, and GALEX photometry using spectral energy distribution fitting. Morphological properties for all BCGs were derived by Sersic profile fitting using the software package SIGMA in different optical bands (g,r,i). We combined this catalog with the BCGs of galaxy groups and clusters extracted from the deeper AEGIS, CDFS, COSMOS, XMM-CFHTLS, and XMM-XXL surveys to study the stellar mass - halo mass relation using the largest sample of X-ray groups and clusters known to date. This result suggests that the mass growth of the central galaxy is controlled by the hierarchical mass growth of the host halo. We find a strong correlation between the stellar mass of BCGs and the mass of their host halos. This relation shows no evolution since z $\sim$ 0.65. We measure a mean scatter of 0.21 and 0.25 for the stellar mass of BCGs in a given halo mass at low ( $0.1<z < 0.3$ ) and high ( $0.3<z<0.65$ ) redshifts, respectively. We further demonstrate that the BCG mass is covariant with the richness of the host halos in the very X-ray luminous systems. We also find evidence that part of the scatter between X-ray luminosity and richness can be reduced by considering stellar mass as an additional variable.

The correlation observed between monochromatic X-ray and UV luminosities in radiatively-efficient active galactic nuclei (AGN) lacks a clear theoretical explanation despite being used for many applications. Such a correlation, with its small intrinsic scatter and its slope that is smaller than unity in log space, represents the compelling evidence that a mechanism regulating the energetic interaction between the accretion disk and the X-ray corona must be in place. This ensures that going from fainter to brighter sources the coronal emission increases less than the disk emission. We discuss here a self-consistently coupled disk-corona model that can identify this regulating mechanism in terms of modified viscosity prescriptions in the accretion disk. The model predicts a lower fraction of accretion power dissipated in the corona for higher accretion states. We then present a quantitative observational test of the model using a reference sample of broad-line AGN and modeling the disk-corona emission for each source in the $L_X-L_{UV}$ plane. We used the slope, normalization, and scatter of the observed relation to constrain the parameters of the theoretical model. For non-spinning black holes and static coronae, we find that the accretion prescriptions that match the observed slope of the $L_X-L_{UV}$ relation produce X-rays that are too weak with respect to the normalization of the observed relation. Instead, considering moderately-outflowing Comptonizing coronae and/or a more realistic high-spinning black hole population significantly relax the tension between the strength of the observed and modeled X-ray emission, while also predicting very low intrinsic scatter in the $L_X-L_{UV}$ relation. In particular, this latter scenario traces a known selection effect of flux-limited samples that preferentially select high-spinning, hence brighter, sources.

We analyze extensive spectroscopic and photometric data of the hypervariable quasar SDSS J131424+530527 (RMID 017) at z=0.456, an optical "changing look" quasar from the Sloan Digital Sky Survey Reverberation Mapping project that increased in optical luminosity by a factor of 10 between 2014 and 2017. The observed broad emission lines all respond in luminosity and width to the changing optical continuum, as expected for photoionization in a stratified, virialized broad emission line region. The luminosity changes therefore result from intrinsic changes in accretion power rather than variable obscuration. The variability is continuous and apparently stochastic, disfavoring an origin as a discrete event such as a tidal disruption flare or microlensing event. It is coordinated on day timescales with blue leading red, consistent with reprocessing powering the entire optical SED. We show that this process cannot work in a standard thin disk geometry on energetic grounds, and would instead require a large covering factor reprocessor. Disk instability models could potentially also explain the data, provided that the instability sets in near the inner radius of a geometrically thick accretion disk.

In the context of the upcoming SRG/eROSITA survey, we present an N-body simulation-based mock catalogue for X-ray selected AGN samples. The model reproduces the observed hard X-ray AGN luminosity function (XLF) and the soft X-ray logN-logS from redshift 0 to 6. The XLF is reproduced to within $\pm5\%$ and the logN-logS to within $\pm20\%$. We develop a joint X-ray -- optical extinction and classification model. We adopt a set of empirical spectral energy distributions to predict observed magnitudes in the UV, optical and NIR. With the latest eROSITA all sky survey sensitivity model, we create a high-fidelity full-sky mock catalogue of X-ray AGN. It predicts their distributions in right ascension, declination, redshift and fluxes. Using empirical medium resolution optical spectral templates and an exposure time calculator, we find that $1.1\times10^6$ ($4\times10^5$) fiber-hours are needed to follow-up spectroscopically from the ground the detected X-ray AGN with an optical magnitude $21<r<22.8$ ($22.8<r<25$) with a 4-m (8-m) class multi-object spectroscopic facility. We find that future clustering studies will measure the AGN bias to the percent level at redshift $z<1.2$ and should discriminate possible scenarios of galaxy-AGN co-evolution. We predict the accuracy to which the baryon acoustic oscillation standard ruler will be measured using X-ray AGN: better than 3\% for AGN between redshift 0.5 to 3 and better than 1\% using the Ly$\alpha$ forest of X-ray QSOs discovered between redshift 2 and 3. eROSITA will provide an outstanding set of targets for future galaxy evolution and cosmological studies.

This work presents the catalogue of optical spectral properties for all X-ray selected SPIDERS active galactic nuclei (AGN) up to SDSS DR14. SPIDERS (SPectroscopic IDentification of eROSITA Sources) is an SDSS-IV programme that is currently conducting optical spectroscopy of the counterparts to the X-ray selected sources detected in the ROSAT all-sky survey and the XMM-Newton slew survey in the footprint of the Extended Baryon Oscillation Spectroscopic Survey (eBOSS). The SPIDERS DR14 sample is the largest sample of X-ray selected AGN with optical spectroscopic follow-up to date. The catalogue presented here is based on a clean sample of 7344 2RXS ($\rm \bar{z}$ = 0.5) and 1157 XMM-Newton slew survey ($\rm \bar{z}$ = 0.4) type 1 AGN with spectroscopic coverage of the H$\rm \beta$ and/or MgII emission lines. Visual inspection results for each object in this sample are available from a combination of literature sources and the SPIDERS group, which provide both reliable redshifts and source classifications. The spectral regions around the H$\rm \beta$ and MgII emission lines have been fit in order to measure both line and continuum properties, estimate bolometric luminosities, and provide black hole mass estimates using the single-epoch (or photoionisation) method. The use of both H$\rm \beta$ and MgII allows black hole masses to be estimated up to z $\rm \simeq$ 2.5. A comparison is made between the spectral properties and black hole mass estimates derived from H$\rm \beta$ and MgII using the subsample of objects which have coverage of both lines in their spectrum. These results have been made publicly available as an SDSS-IV DR14 value added catalogue.

Groups and clusters of galaxies are a current focus of astronomical research owing to their role in determining the environmental effects on galaxies and the constraints they provide to cosmology. The eROSITA X-ray telescope on board the Spectrum Roentgen Gamma observatory will be launched in 2019 and will have completed eight scans of the full sky when 4MOST starts operating. The experiment will detect groups and clusters of galaxies through X-ray emission from the hot intergalactic medium. The purpose of the 4MOST eROSITA Galaxy Cluster Redshift Survey is to provide spectroscopic redshifts of the optical counterparts to the X-ray emission from 40,000 groups and clusters of galaxies so as to perform dynamical estimates of the total mass and to measure the properties of the member galaxies. The survey aims to obtain precise redshift measurements of the photometrically identified brightest cluster galaxies at redshift $z > 0.7$. At lower redshifts ($z < 0.7$), the programme aims to sample over 15 member galaxies per cluster and enable dynamical mass measurements to calibrate the clusters for cosmological experiments. At $z < 0.2$, eROSITA will also detect X-ray emission from galaxy groups and filaments. 4MOST spectroscopic data from the survey will be used for optical identification of galaxy groups down to eROSITA's mass detection limits of $10^{13} M_\odot$, as well as the detection of the largest filaments for pioneering studies of their X-ray emission.

X-ray and mid-infrared emission are signposts of the accretion of matter onto the supermassive black holes that reside at the centres of most galaxies. As a major step towards understanding accreting supermassive black holes and their role in the evolution of galaxies, we will use the 4MOST multi-object spectrograph to provide a highly complete census of active galactic nuclei over a large fraction of the extragalactic sky observed in X-rays by eROSITA that is visible to 4MOST. We will systematically follow up all eROSITA point-like extragalactic X-ray sources (mostly active galactic nuclei), and complement them with a heavily obscured active galactic nuclei selection approach using mid-infrared data from the Wide-field Infrared Survey Explorer (WISE). The X-ray and mid-infrared flux limits of eROSITA and WISE are well matched to the spectroscopic capabilities of a 4-metre-class telescope, allowing us to reach completeness levels of ~80-90% for all X-ray selected active galactic nuclei with fluxes $f_{0.5-2\ {\rm keV}} > 10^{-14}$ erg s$^{-1}$ cm$^{-2}$; this is about a factor of 30 deeper than the ROSAT all-sky survey. With these data we will determine the physical properties (redshift, luminosity, line emission strength, masses, etc.) of up to one million supermassive black holes, constrain their cosmic evolution and clustering properties, and explore the connection between active galactic nuclei and large-scale structure over redshifts $0 \le z \le 6$.

The current status of and motivation for the 4MOST Survey Strategy, as developed by the Consortium science team, are presented here. Key elements of the strategy are described, such as sky coverage, number of visits and total exposure times in different parts of the sky, and how to deal with different observing conditions. The task of organising the strategy is not simple, with many different surveys that have vastly different target brightnesses and densities, sample completeness levels, and signal-to-noise requirements. We introduce here a number of concepts that we will use to ensure all surveys are optimised. Astronomers who are planning to submit a Participating Survey proposal are strongly encouraged to read this article and any relevant 4MOST Survey articles in this issue of The Messenger such that they can optimally complement and benefit from the planned surveys of the 4MOST Consortium.

The 4MOST instrument is a multi-object spectrograph that will address Galactic and extragalactic science cases simultaneously by observing targets from a large number of different surveys within each science exposure. This parallel mode of operation and the survey nature of 4MOST require some distinct 4MOST-specific operational features within the overall operations model of ESO. The main feature is that the 4MOST Consortium will deliver, not only the instrument, but also contractual services to the user community, which is why 4MOST is also described as a facility. This white paper concentrates on information particularly useful to answering the forthcoming Call for Letters of Intent.

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

We determine the stellar population properties - age, metallicity, dust reddening, stellar mass and the star formation history - for all spectra classified as galaxies that were published by the Sloan Digital Sky Survey (SDSS data release 14) and by the DEEP2 (data release 4) galaxy surveys. We perform full spectral fitting on individual spectra, making use of high spectral resolution stellar population models. Calculations are carried out for several choices of the model input, including three stellar initial mass functions and three input stellar libraries to the models. We study the accuracy of parameter derivation, in particular the stellar mass, as a function of the signal-to-noise of the galaxy spectra. We find that at low redshift, a signal to noise ratio per pixel around 20 (5) allows a statistical accuracy on $\log_{10}(M^{*}/M_{\odot})$ of 0.2 (0.4) dex, for the Chabrier IMF. For the first time, we study DEEP2 galaxies selected by their \OII luminosity in the redshift range $0.83<z<1.03$, finding that they are consistent with a flat number density in stellar mass in the range $10^9<M/M_{\odot}<10^{11.5}$. We find the resulting stellar mass function based on SDSS or eBOSS in agreement with previous studies (Maraston et al. 2013). We publish all catalogs of properties as well as model spectra of the continuum for these galaxies as a value added catalog of the fourteenth data release of the SDSS. This catalog is about twice as large as its predecessors (DR12) and will aid a variety of studies on galaxy evolution and cosmology.

The source-subtracted cosmic infrared background (CIB) fluctuations uncovered in deep Spitzer data cannot be explained by known galaxy populations and appear strongly coherent with unresolved cosmic X-ray background (CXB). This suggests that the source-subtracted CIB contains emissions from significantly abundant accreting black holes (BHs). We show that theoretically such populations would have the angular power spectrum which is largely independent of the epochs occupied by these sources, provided they are at z>~ 4, offering an important test of the origin of the new populations. Using the current measurements we reconstruct the underlying soft X-ray CXB from the new sources and show that its fluctuations, while consistent with a high-z origin, have an amplitude that cannot be reached in direct measurements with the foreseeable X-ray space missions. This necessitates application of the methods developed by the authors to future IR and X-ray datasets, which must cover large areas of the sky in order to measure the signal with high precision. The LIBRAE project within ESA's Euclid mission will probe source-subtracted CIB over ~1/2 the sky at three near-IR bands, and its cross-power with unresolved CXB can be measured then from the concurrent eROSITA mission covering the same areas of the sky. We discuss the required methodology for this measurement and evaluate its projected S/N to show the unique potential of this experimental configuration to accurately probe the CXB from the new BH sources and help identify their epochs.

Twenty years have passed since first light for the Sloan Digital Sky Survey (SDSS). Here, we release data taken by the fourth phase of SDSS (SDSS-IV) across its first three years of operation (July 2014-July 2017). This is the third data release for SDSS-IV, and the fifteenth from SDSS (Data Release Fifteen; DR15). New data come from MaNGA - we release 4824 datacubes, as well as the first stellar spectra in the MaNGA Stellar Library (MaStar), the first set of survey-supported analysis products (e.g. stellar and gas kinematics, emission line, and other maps) from the MaNGA Data Analysis Pipeline (DAP), and a new data visualisation and access tool we call "Marvin". The next data release, DR16, will include new data from both APOGEE-2 and eBOSS; those surveys release no new data here, but we document updates and corrections to their data processing pipelines. The release is cumulative; it also includes the most recent reductions and calibrations of all data taken by SDSS since first light. In this paper we describe the location and format of the data and tools and cite technical references describing how it was obtained and processed. The SDSS website (www.sdss.org) has also been updated, providing links to data downloads, tutorials and examples of data use. While SDSS-IV will continue to collect astronomical data until 2020, and will be followed by SDSS-V (2020-2025), we end this paper by describing plans to ensure the sustainability of the SDSS data archive for many years beyond the collection of data.

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

The LOFAR Two-metre Sky Survey (LoTSS) is an ongoing sensitive, high-resolution 120-168MHz survey of the entire northern sky for which observations are now 20% complete. We present our first full-quality public data release. For this data release 424 square degrees, or 2% of the eventual coverage, in the region of the HETDEX Spring Field (right ascension 10h45m00s to 15h30m00s and declination 45$^\circ$00$'$00$''$ to 57$^\circ$00$'$00$''$) were mapped using a fully automated direction-dependent calibration and imaging pipeline that we developed. A total of 325,694 sources are detected with a signal of at least five times the noise, and the source density is a factor of $\sim 10$ higher than the most sensitive existing very wide-area radio-continuum surveys. The median sensitivity is S$_{\rm 144 MHz} = 71\,\mu$Jy beam$^{-1}$ and the point-source completeness is 90% at an integrated flux density of 0.45mJy. The resolution of the images is 6$''$ and the positional accuracy is within 0.2$''$. This data release consists of a catalogue containing location, flux, and shape estimates together with 58 mosaic images that cover the catalogued area. In this paper we provide an overview of the data release with a focus on the processing of the LOFAR data and the characteristics of the resulting images. In two accompanying papers we provide the radio source associations and deblending and, where possible, the optical identifications of the radio sources together with the photometric redshifts and properties of the host galaxies. These data release papers are published together with a further $\sim$20 articles that highlight the scientific potential of LoTSS.

In this study, we present a modelling of the X-ray emission from the simulated SMBHs within the cosmological hydrodynamical Magneticum Pathfinder Simulation, in order to study the statistical properties of the resulting X-ray Active Galactic Nuclei (AGN) population and their expected contribution to the X-ray flux from galaxy clusters. The simulations reproduce the evolution of the observed unabsorbed AGN bolometric luminosity functions (LFs) up to redshift z~2, consistently with previous works. Furthermore, we study the evolution of the LFs in the soft (SXR) and hard (HXR) X-ray bands by means of synthetic X-ray data generated with the PHOX simulator, that includes an observationally-motivated modelling of an instrinsic absorption component, mimicking the torus around the AGN. The reconstructed SXR and HXR AGN LFs present a remarkable agreement with observational data up to z~2 when an additional obscuration fraction for Compton-thick AGN is assumed, although a discrepancy still exists for the SXR LF at z=2.3. With this approach, we also generate full eROSITA mock observations to predict the level of contamination due to AGN of the intra-cluster medium (ICM) X-ray emission, which can affect cluster detection especially at high redshifts. We find that, at z~1-1.5, for 20-40% of the clusters with M500>3e13 Msun/h, the AGN counts in the observed SXR band exceed by more than a factor of 2 the counts from the whole ICM.

Theoretical models of galaxy formation suggest that the presence of an active galactic nucleus (AGN) is required to regulate the growth of its host galaxy through feedback mechanisms, produced by e.g. AGN-driven outflows. Although such outflows are common both at low and high redshift, a comprehensive picture is still missing. The peak epoch of galaxy assembly (1<z<3) has been poorly explored so far, and current observations in this redshift range are mostly limited to targets with high chances to be in an outflowing phase. This paper introduces SUPER (a SINFONI Survey for Unveiling the Physics and Effect of Radiative feedback), an ongoing ESO's VLT/SINFONI Large Programme. SUPER will perform the first systematic investigation of ionized outflows in a sizeable and blindly-selected sample of 39 X-ray AGN at z~2, which reaches high spatial resolutions (~2 kpc) thanks to the adaptive optics-assisted IFU observations. The outflow morphology and star formation in the host galaxy will be mapped through the broad component of [OIII] and the narrow component of Ha emission lines. The main aim of our survey is to infer the impact of outflows on the on-going star formation and to link the outflow properties to a number of AGN and host galaxy properties. We describe here the survey characteristics and goals, as well as the selection of the target sample. Moreover, we present a full characterization of its multi-wavelength properties: we measure, via spectral energy distribution fitting of UV-to-FIR photometry, stellar masses (4x10^9-2x10^11 Msun), star formation rates (25-680 Msun yr^-1) and AGN bolometric luminosities (2x10^44-8x10^47 erg s^-1), along with obscuring column densities (<2x10^24 cm^-2) and 2-10 keV luminosities (2x10^43-6x10^45 erg s^-1) derived through X-ray spectral analysis. Finally, we classify our AGN as jetted or non-jetted according to their radio and FIR emission.

We report on the discovery of a dramatic X-ray spectral variability event observed in a $z\sim 1$ broad line type-1 QSO. The XMM-Newton spectrum from the year 2000 is characterized by an unobscured power-law spectrum with photon index of $\Gamma\sim 2$, a column density of $N_{\mathrm{H}}\sim 5\times 10^{20}\,\mathrm{cm^{-2}}$, and no prominent reflection component. Five years later, Chandra captured the source in a heavily-obscured, reflection-dominated state. The observed X-ray spectral variability could be caused by a Compton-thick cloud with $N_{\mathrm{H}}\sim 2\times 10^{24}\,\mathrm{cm^{-2}}$ eclipsing the direct emission of the hot corona, implying an extreme $N_{\mathrm{H}}$ variation never before observed in a type-1 QSO. An alternative scenario is a corona that switched off in between the observations. In addition, both explanations require a significant change of the X-ray luminosity prior to the obscuration or fading of the corona and/or a change of the relative geometry of the source/reflector system. Dramatic X-ray spectral variability of this kind could be quite common in type-1 QSOs, considering the relatively few datasets in which such an event could have been identified. Our analysis implies that there may be a population of type-1 QSOs which are Compton-thick in the X-rays when observed at any given time.

The Advanced Telescope for High Energy Astrophysics (Athena) is the X-ray observatory large mission selected by the European Space Agency (ESA), within its Cosmic Vision 2015-2025 programme, to address the "Hot and Energetic Universe" scientific theme (Nandra et al. 2013), and it is provisionally due for launch in the early 2030s. The Square Kilometer Array (SKA) is the next generation radio observatory and consists of two telescopes, one comprised of dishes operating at mid frequencies (SKA1-MID) and located in South Africa, and the other comprised of Log-Periodic antennas operating at low radio frequencies (SKA1-LOW), which will be located in Australia (Braun et al. 2017). The scientific commissioning of the radio telescope is planned to begin in 2021-2022. The SKA-Athena Synergy Team (SAST) has been tasked to single out the potential scientific synergies between Athena and SKA. The astrophysical community was involved in this exercise primarily through a dedicated SKA-Athena Synergy Workshop, which took place on April 24-25, 2017 at SKAO, Jodrell Bank, Manchester. The final result of the synergy exercise, this White Paper, describes in detail a number of scientific opportunities that will be opened up by the combination of Athena and SKA, these include: 1. the Cosmic Dawn; 2. the Evolution of black holes and galaxies; 3. Active galaxy feedback in galaxy clusters; 4. Non-thermal phenomena in galaxy clusters; 5. Detecting the cosmic web; 6. Black-hole accretion physics and astrophysical transients; 7. Galactic astronomy: stars, planets, pulsars and supernovae.

We forecast the potential of the forthcoming X-ray galaxy-cluster survey with eROSITA to constrain dark-energy models. We focus on spatially-flat cosmological scenarios with either constant or time-dependent dark-energy equation-of-state parameters. Fisher information is extracted from the number density and spatial clustering of a photon-count-limited sample of clusters of galaxies up to z~2. We consider different scenarios for the availability of (i) X-ray follow-up observations, (ii) photometric and spectroscopic redshifts, and (iii) accurate knowledge of the observable -- mass relation down to the scale of galaxy groups. With about 125,000 clusters (detected with more than 50 photons and with mass M500c > $10^{13} h^{-1}$ Msun) from an average all-sky exposure of 1.6 ks, eROSITA will give marginalized, one-dimensional, 1-$\sigma$ errors of $\Delta \sigma_8 = \pm~0.008$ ($\sim$1 per cent), $\Delta \Omega_{\rm M} = \pm~0.006$ (2.2 per cent), $\Delta w_0 = \pm~0.07$ (7 per cent), and $\Delta w_a = \pm~0.25$ (optimistic scenario) in combination with (and largely improving upon) current constraints from various cosmological probes (cosmic microwave background, BAOs, Type Ia SNe). Our findings correspond to a dark-energy figure of merit in the range of $116-162$ (after the four years of all-sky survey), making eROSITA one of the first Stage IV experiments to come on line according to the classification of the Dark Energy Task Force. To secure improved mass calibrations and to include high-redshift clusters (z > 0.5) as well as objects at the group-mass scale (M500c < 5 $\times~ 10^{13} h^{-1}$ Msun) will be vital to reach such accuracies.

Using the X-ray-selected active galactic nuclei (AGN) from the XMM-XXL north survey and the SDSS Baryon Oscillation Spectroscopic Survey (BOSS) spectroscopic follow-up of them, we compare the properties of X-ray unobscured and obscured broad-line AGN (BLAGN1 and BLAGN2; $N_\textrm{H}$below and above $10^{21.5}$ cm$^{-2}$), including their X-ray luminosity $L_X$, black hole mass, Eddington ratio $\lambda_{\textrm{Edd}}$, optical continuum and line features. We find that BLAGN2 have systematically larger broad line widths and hence apparently higher (lower) $M_{\textrm{BH}}$ ($\lambda_{\textrm{Edd}}$) than BLAGN1. We also find that the X-ray obscuration in BLAGN tends to coincide with optical dust extinction, which is optically thinner than that in narrow-line AGN (NLAGN) and likely partial-covering to the broad line region. All the results can be explained in the framework of a multi-component, clumpy torus model by interpreting BLAGN2 as an intermediate type between BLAGN1 and NLAGN in terms of an intermediate inclination angle.

Studies of galaxy clusters provide stringent constraints on models of structure formation. Provided that selection effects are under control, large X-ray surveys are well suited to derive cosmological parameters, in particular those governing the dark energy equation of state. We forecast the capabilities of the all-sky eROSITA (the extended ROentgen Survey with an Imaging Telescope Array) survey to be achieved by the early 2020s. We bring special attention to modeling the entire chain from photon emission to source detection and cataloguing. The selection function of galaxy clusters for the upcoming eROSITA mission is investigated by means of extensive and dedicated Monte-Carlo simulations. Employing a combination of accurate instrument characterization and of state-of-the-art source detection technique, we determine a cluster detection efficiency based on the cluster fluxes and sizes. Using this eROSITA cluster selection function, we find that eROSITA will detect a total of $\sim 10^5$ clusters in the extra-galactic sky. This number of clusters will allow eROSITA to put stringent constraints on cosmological models. We show that incomplete assumptions on selection effects, such as neglecting the distribution of cluster sizes, induce a bias in the derived value of cosmological parameters. Synthetic simulations of the eROSITA sky capture the essential characteristics impacting the next-generation galaxy cluster surveys and they highlight parameters requiring tight monitoring in order to avoid biases in cosmological analyses.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) has been in operation since July 2014. This paper describes the second data release from this phase, and the fourteenth from SDSS overall (making this, Data Release Fourteen or DR14). This release makes public data taken by SDSS-IV in its first two years of operation (July 2014-2016). Like all previous SDSS releases, DR14 is cumulative, including the most recent reductions and calibrations of all data taken by SDSS since the first phase began operations in 2000. New in DR14 is the first public release of data from the extended Baryon Oscillation Spectroscopic Survey (eBOSS); the first data from the second phase of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE-2), including stellar parameter estimates from an innovative data driven machine learning algorithm known as "The Cannon"; and almost twice as many data cubes from the Mapping Nearby Galaxies at APO (MaNGA) survey as were in the previous release (N = 2812 in total). This paper describes the location and format of the publicly available data from SDSS-IV surveys. We provide references to the important technical papers describing how these data have been taken (both targeting and observation details) and processed for scientific use. The SDSS website (www.sdss.org) has been updated for this release, and provides links to data downloads, as well as tutorials and examples of data use. SDSS-IV is planning to continue to collect astronomical data until 2020, and will be followed by SDSS-V.

We present the Data Release 14 Quasar catalog (DR14Q) from the extended Baryon Oscillation Spectroscopic Survey (eBOSS) of the Sloan Digital Sky Survey IV (SDSS-IV). This catalog includes all SDSS-IV/eBOSS objects that were spectroscopically targeted as quasar candidates and that are confirmed as quasars via a new automated procedure combined with a partial visual inspection of spectra, have luminosities $M_{\rm i} \left[ z=2 \right] < -20.5$ (in a $\Lambda$CDM cosmology with $H_0 = 70 \ {\rm km \ s^{-1} \ Mpc ^{-1}}$, $\Omega_{\rm M} = 0.3$, and $\Omega_{\rm \Lambda} = 0.7$), and either display at least one emission line with a full width at half maximum (FWHM) larger than $500 \ {\rm km \ s^{-1}}$ or, if not, have interesting/complex absorption features. The catalog also includes previously spectroscopically-confirmed quasars from SDSS-I, II and III. The catalog contains 526,356 quasars 144,046 are new discoveries since the beginning of SDSS-IV) detected over 9,376 deg$^2$ (2,044 deg$^2$ having new spectroscopic data available) with robust identification and redshift measured by a combination of principal component eigenspectra. The catalog is estimated to have about 0.5% contamination. The catalog identifies 21,877 broad absorption line quasars and lists their characteristics. For each object, the catalog presents SDSS five-band CCD-based photometry with typical accuracy of 0.03 mag. The catalog also contains X-ray, ultraviolet, near-infrared, and radio emission properties of the quasars, when available, from other large-area surveys.

SDSS-V will be an all-sky, multi-epoch spectroscopic survey of over six million objects. It is designed to decode the history of the Milky Way, trace the emergence of the chemical elements, reveal the inner workings of stars, and investigate the origin of planets. It will also create an integral-field spectroscopic map of the gas in the Galaxy and the Local Group that is 1,000x larger than the current state of the art and at high enough spatial resolution to reveal the self-regulation mechanisms of galactic ecosystems. SDSS-V will pioneer systematic, spectroscopic monitoring across the whole sky, revealing changes on timescales from 20 minutes to 20 years. The survey will thus track the flickers, flares, and radical transformations of the most luminous persistent objects in the universe: massive black holes growing at the centers of galaxies. The scope and flexibility of SDSS-V will be unique among extant and future spectroscopic surveys: it is all-sky, with matched survey infrastructures in both hemispheres; it provides near-IR and optical multi-object fiber spectroscopy that is rapidly reconfigurable to serve high target densities, targets of opportunity, and time-domain monitoring; and it provides optical, ultra-wide-field integral field spectroscopy. SDSS-V, with its programs anticipated to start in 2020, will be well-timed to multiply the scientific output from major space missions (e.g., TESS, Gaia, eROSITA) and ground-based projects. SDSS-V builds on the 25-year heritage of SDSS's advances in data analysis, collaboration infrastructure, and product deliverables. The project is now refining its science scope, optimizing the survey strategies, and developing new hardware that builds on the SDSS-IV infrastructure. We present here an overview of the current state of these developments as we seek to build our worldwide consortium of institutional and individual members.

In this paper, we investigate 2727 galaxies observed by MaNGA as of June 2016 to develop spatially resolved techniques for identifying signatures of active galactic nuclei (AGN). We identify 303 AGN candidates. The additional spatial dimension imposes challenges in identifying AGN due to contamination from diffuse ionized gas, extra-planar gas and photoionization by hot stars. We show that the combination of spatially-resolved line diagnostic diagrams and additional cuts on H$\alpha$ surface brighness and H$\alpha$ equivalent width can distinguish between AGN-like signatures and high-metallicity galaxies with LINER-like spectra. Low mass galaxies with high specific star formation rates are particularly difficult to diagnose and routinely show diagnostic line ratios outside of the standard star-formation locus. We develop a new diagnostic -- the distance from the standard diagnostic line in the line-ratios space -- to evaluate the significance of the deviation from the star-formation locus. We find 173 galaxies that would not have been selected as AGN candidates based on single-fibre spectral measurements but exhibit photoionization signatures suggestive of AGN activity in the MaNGA resolved observations, underscoring the power of large integral field unit (IFU) surveys. A complete census of these new AGN candidates is necessary to understand their nature and probe the complex co-evolution of supermassive black holes and their hosts.

We release the AllWISE counterparts and Gaia matches to 106,573 and 17,665 X-ray sources detected in the ROSAT 2RXS and XMMSL2 surveys with |b|>15. These are the brightest X-ray sources in the sky, but their position uncertainties and the sparse multi-wavelength coverage until now rendered the identification of their counterparts a demanding task with uncertain results. New all-sky multi-wavelength surveys of sufficient depth, like AllWISE and Gaia, and a new Bayesian statistics based algorithm, NWAY, allow us, for the first time, to provide reliable counterpart associations. NWAY extends previous distance and sky density based association methods and, using one or more priors (e.g., colors, magnitudes), weights the probability that sources from two or more catalogues are simultaneously associated on the basis of their observable characteristics. Here, counterparts have been determined using a WISE color-magnitude prior. A reference sample of 4524 XMM/Chandra and Swift X-ray sources demonstrates a reliability of ~ 94.7% (2RXS) and 97.4% (XMMSL2). Combining our results with Chandra-COSMOS data, we propose a new separation between stars and AGN in the X-ray/WISE flux-magnitude plane, valid over six orders of magnitude. We also release the NWAY code and its user manual. NWAY was extensively tested with XMM-COSMOS data. Using two different sets of priors, we find an agreement of 96% and 99% with published Likelihood Ratio methods. Our results were achieved faster and without any follow-up visual inspection. With the advent of deep and wide area surveys in X-rays (e.g. SRG/eROSITA, Athena/WFI) and radio (ASKAP/EMU, LOFAR, APERTIF, etc.) NWAY will provide a powerful and reliable counterpart identification tool.

Multi-wavelength surveys covering large sky volumes are necessary to obtain an accurate census of rare objects such as high luminosity and/or high redshift active galactic nuclei (AGN). Stripe 82X is a 31.3 deg$^2$ X-ray survey with $Chandra$ and $XMM$-Newton observations overlapping the legacy Sloan Digital Sky Survey (SDSS) Stripe 82 field, which has a rich investment of multi-wavelength coverage from the ultraviolet to the radio. The wide-area nature of this survey presents new challenges for photometric redshifts for AGN compared to previous work on narrow-deep fields because it probes different populations of objects that need to be identified and represented in the library of templates. Here we present an updated X-ray plus multi-wavelength matched catalog, including $Spitzer$ counterparts, and estimated photometric redshifts for 5961 (96% of a total of 6181) X-ray sources, which have a normalized median absolute deviation, $\sigma_{\rm nmad}$ = 0.06 and an outlier fraction, $\eta$ = 13.7%. The populations found in this survey, and the template libraries used for photometric redshifts, provide important guiding principles for upcoming large-area surveys such as $eROSITA$ and 3$XMM$ (in X-ray) and the Large Synoptic Survey Telescope (LSST; optical).

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) began observations in July 2014. It pursues three core programs: APOGEE-2, MaNGA, and eBOSS. In addition, eBOSS contains two major subprograms: TDSS and SPIDERS. This paper describes the first data release from SDSS-IV, Data Release 13 (DR13), which contains new data, reanalysis of existing data sets and, like all SDSS data releases, is inclusive of previously released data. DR13 makes publicly available 1390 spatially resolved integral field unit observations of nearby galaxies from MaNGA, the first data released from this survey. It includes new observations from eBOSS, completing SEQUELS. In addition to targeting galaxies and quasars, SEQUELS also targeted variability-selected objects from TDSS and X-ray selected objects from SPIDERS. DR13 includes new reductions of the SDSS-III BOSS data, improving the spectrophotometric calibration and redshift classification. DR13 releases new reductions of the APOGEE-1 data from SDSS-III, with abundances of elements not previously included and improved stellar parameters for dwarf stars and cooler stars. For the SDSS imaging data, DR13 provides new, more robust and precise photometric calibrations. Several value-added catalogs are being released in tandem with DR13, in particular target catalogs relevant for eBOSS, TDSS, and SPIDERS, and an updated red-clump catalog for APOGEE. This paper describes the location and format of the data now publicly available, as well as providing references to the important technical papers that describe the targeting, observing, and data reduction. The SDSS website, http://www.sdss.org, provides links to the data, tutorials and examples of data access, and extensive documentation of the reduction and analysis procedures. DR13 is the first of a scheduled set that will contain new data and analyses from the planned ~6-year operations of SDSS-IV.

We describe the Sloan Digital Sky Survey IV (SDSS-IV), a project encompassing three major spectroscopic programs. The Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2) is observing hundreds of thousands of Milky Way stars at high resolution and high signal-to-noise ratio in the near-infrared. The Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey is obtaining spatially-resolved spectroscopy for thousands of nearby galaxies (median redshift of z = 0.03). The extended Baryon Oscillation Spectroscopic Survey (eBOSS) is mapping the galaxy, quasar, and neutral gas distributions between redshifts z = 0.6 and 3.5 to constrain cosmology using baryon acoustic oscillations, redshift space distortions, and the shape of the power spectrum. Within eBOSS, we are conducting two major subprograms: the SPectroscopic IDentification of eROSITA Sources (SPIDERS), investigating X-ray AGN and galaxies in X-ray clusters, and the Time Domain Spectroscopic Survey (TDSS), obtaining spectra of variable sources. All programs use the 2.5-meter Sloan Foundation Telescope at Apache Point Observatory; observations there began in Summer 2014. APOGEE-2 also operates a second near-infrared spectrograph at the 2.5-meter du Pont Telescope at Las Campanas Observatory, with observations beginning in early 2017. Observations at both facilities are scheduled to continue through 2020. In keeping with previous SDSS policy, SDSS-IV provides regularly scheduled public data releases; the first one, Data Release 13, was made available in July 2016.

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

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

The wide-area XMM-XXL X-ray survey is used to explore the fraction of obscured AGN at high accretion luminosities, $L_X (\rm 2-10 \, keV) > 10^{44} \, erg \,s ^{-1}$, and out to redshift $z\approx1.5$. The sample covers an area of about $\rm14\,deg^2$ and provides constraints on the space density of powerful AGN over a wide range of neutral hydrogen column densities extending beyond the Compton-thick limit, $\rm N_H\approx10^{24}\,cm^{-2}$. The fraction of obscured Compton-thin ($\rm N_H=10^{22}-10^{24}\,cm^{-2}$) AGN is estimated to be $\approx0.35$ for luminosities $L_X(\rm 2-10\,keV)>10^{44}\,erg\,s^{-1}$ independent of redshift. For less luminous sources the fraction of obscured Compton-thin AGN increases from $0.45\pm0.10$ at $z=0.25$ to $0.75\pm0.05$ at $z=1.25$. Studies that select AGN in the infrared via template fits to the observed Spectral Energy Distribution of extragalactic sources estimate space densities at high accretion luminosities consistent with the XMM-XXL constraints. There is no evidence for a large population of AGN (e.g. heavily obscured) identified in the infrared and missed at X-ray wavelengths. We further explore the mid-infrared colours of XMM-XXL AGN as a function of accretion luminosity, column density and redshift. The fraction of XMM-XXL sources that lie within the mid-infrared colour wedges defined in the literature to select AGN is primarily a function of redshift. This fraction increases from about 20-30% at z=0.25 to about 50-70% at $z=1.5$.

The supermassive black holes in active galactic nuclei (AGN) and stellar-mass black holes in X-ray binaries (XRBs) are believed to work in a similar way. While XRBs evolve rapidly and several sources have undergone a few complete cycles from quiescence to an outburst and back, most AGN remain in the same state over periods of decades, due to their longer characteristic timescale proportional to their size. However, the study of the AGN spectral states is still possible with a large sample of sources. Multi-wavelength observations are needed for this purpose since the AGN thermal disc emission dominates in the ultraviolet energy range, while the up-scattered hot-corona emission is detected in X-rays. We compared simultaneous UV and X-ray measurements of AGN obtained by the XMM-Newton satellite. The non-thermal flux was constrained from the 2-12 keV X-ray luminosity, while the thermal disc component was estimated from the UV flux at 2900A. The hardness (ratio between the X-ray and UV plus X-ray luminosity) and the total luminosity were used to construct the AGN state diagrams. For sources with reliable mass measurements, the Eddington ratio was used instead of the total luminosity. The state diagrams show that the radio-loud sources have on average higher hardness, due to the lack of the thermal disc emission in the UV band, and have flatter intrinsic X-ray spectra. In contrast, the sources with high luminosity and low hardness are radio-quiet AGN. The hardness-Eddington ratio diagram reveals that the average radio-loudness is stronger for low-accreting sources, while it decreases when the accretion rate is close to the Eddington limit. Our results indicate that the general properties of AGN accretion states are similar to those of X-ray binaries. This suggests that the AGN radio dichotomy of radio-loud and radio-quiet sources can be explained by the evolution of the accretion states.

SPIDERS (SPectroscopic IDentification of eROSITA Sources) is an SDSS-IV survey running in parallel to the eBOSS cosmology project. SPIDERS will obtain optical spectroscopy for large numbers of X-ray-selected AGN and galaxy cluster members detected in wide area eROSITA, XMM-Newton and ROSAT surveys. We describe the methods used to choose spectroscopic targets for two sub-programmes of SPIDERS: X-ray selected AGN candidates detected in the ROSAT All Sky and the XMM-Newton Slew surveys. We have exploited a Bayesian cross-matching algorithm, guided by priors based on mid-IR colour-magnitude information from the WISE survey, to select the most probable optical counterpart to each X-ray detection. We empirically demonstrate the high fidelity of our counterpart selection method using a reference sample of bright well-localised X-ray sources collated from XMM-Newton, Chandra and Swift-XRT serendipitous catalogues, and also by examining blank-sky locations. We describe the down-selection steps which resulted in the final set of SPIDERS-AGN targets put forward for spectroscopy within the eBOSS/TDSS/SPIDERS survey, and present catalogues of these targets. We also present catalogues of ~12000 ROSAT and ~1500 XMM-Newton Slew survey sources which have existing optical spectroscopy from SDSS-DR12, including the results of our visual inspections. On completion of the SPIDERS program, we expect to have collected homogeneous spectroscopic redshift information over a footprint of ~7500 deg$^2$ for >85 percent of the ROSAT and XMM-Newton Slew survey sources having optical counterparts in the magnitude range 17<r<22.5, producing a large and highly complete sample of bright X-ray-selected AGN suitable for statistical studies of AGN evolution and clustering.

We present the first fully simultaneous fits to the NIR and X-ray spectral slope (and its evolution) during a very bright flare from Sgr A*, the supermassive black hole at the Milky Way's center. Our study arises from ambitious multi-wavelength monitoring campaigns with XMM-Newton, NuSTAR and SINFONI. The average multi-wavelength spectrum is well reproduced by a broken power-law with $\Gamma_{NIR}=1.7\pm0.1$ and $\Gamma_X=2.27\pm0.12$. The difference in spectral slopes ($\Delta\Gamma=0.57\pm0.09$) strongly supports synchrotron emission with a cooling break. The flare starts first in the NIR with a flat and bright NIR spectrum, while X-ray radiation is detected only after about $10^3$ s, when a very steep X-ray spectrum ($\Delta\Gamma=1.8\pm0.4$) is observed. These measurements are consistent with synchrotron emission with a cooling break and they suggest that the high energy cut-off in the electron distribution ($\gamma_{max}$) induces an initial cut-off in the optical-UV band that evolves slowly into the X-ray band. The temporal and spectral evolution observed in all bright X-ray flares are also in line with a slow evolution of $\gamma_{max}$. We also observe hints for a variation of the cooling break that might be induced by an evolution of the magnetic field (from $B\sim30\pm8$ G to $B\sim4.8\pm1.7$ G at the X-ray peak). Such drop of the magnetic field at the flare peak would be expected if the acceleration mechanism is tapping energy from the magnetic field, such as in magnetic reconnection. We conclude that synchrotron emission with a cooling break is a viable process for Sgr A*'s flaring emission.

We report on the diversity in quasar spectra from the Baryon Oscillation Spectroscopic Survey. After filtering the spectra to mitigate selection effects and Malmquist bias associated with a nearly flux-limited sample, we create high signal-to-noise ratio composite spectra from 58,656 quasars (2.1 \le z \le 3.5), binned by luminosity, spectral index, and redshift. With these composite spectra, we confirm the traditional Baldwin effect (BE, i.e., the anticorrelation of C IV equivalent width (EW) and luminosity) that follows the relation W_\lambda \propto L^{\beta_w} with slope \beta_w = -0.35 \pm 0.004, -0.35 \pm 0.005, and -0.41 \pm 0.005 for z = 2.25, 2.46, and 2.84, respectively. In addition to the redshift evolution in the slope of the BE, we find redshift evolution in average quasar spectral features at fixed luminosity. The spectroscopic signature of the redshift evolution is correlated at 98% with the signature of varying luminosity, indicating that they arise from the same physical mechanism. At a fixed luminosity, the average C IV FWHM decreases with increasing redshift and is anti-correlated with C IV EW. The spectroscopic signature associated with C IV FWHM suggests that the trends in luminosity and redshift are likely caused by a superposition of effects that are related to black hole mass and Eddington ratio. The redshift evolution is the consequence of a changing balance between these two quantities as quasars evolve toward a population with lower typical accretion rates at a given black hole mass.

SPIDERS (The SPectroscopic IDentification of eROSITA Sources) is a program dedicated to the homogeneous and complete spectroscopic follow-up of X-ray AGN and galaxy clusters over a large area ($\sim$7500 deg$^2$) of the extragalactic sky. SPIDERS is part of the SDSS-IV project, together with the Extended Baryon Oscillation Spectroscopic Survey (eBOSS) and the Time-Domain Spectroscopic Survey (TDSS). This paper describes the largest project within SPIDERS before the launch of eROSITA: an optical spectroscopic survey of X-ray selected, massive ($\sim 10^{14}$ to $10^{15}~M_{\odot}$) galaxy clusters discovered in ROSAT and XMM-Newton imaging. The immediate aim is to determine precise ($\Delta_z \sim 0.001$) redshifts for 4,000-5,000 of these systems out to $z \sim 0.6$. The scientific goal of the program is precision cosmology, using clusters as probes of large-scale structure in the expanding Universe. We present the cluster samples, target selection algorithms and observation strategies. We demonstrate the efficiency of selecting targets using a combination of SDSS imaging data, a robust red-sequence finder and a dedicated prioritization scheme. We describe a set of algorithms and work-flow developed to collate spectra and assign cluster membership, and to deliver catalogues of spectroscopically confirmed clusters. We discuss the relevance of line-of-sight velocity dispersion estimators for the richer systems. We illustrate our techniques by constructing a catalogue of 230 spectroscopically validated clusters ($0.031 < z < 0.658$), found in pilot observations. We discuss two potential science applications of the SPIDERS sample: the study of the X-ray luminosity-velocity dispersion ($L_X-\sigma$) relation and the building of stacked phase-space diagrams.

We present the Data Release 12 Quasar catalog (DR12Q) from the Baryon Oscillation Spectroscopic Survey (BOSS) of the SDSS-III. This catalog includes all SDSS-III/BOSS objects that were spectroscopically targeted as quasar candidates during the full survey and that are confirmed as quasars via visual inspection of the spectra, have luminosities Mi[z=2]<-20.5 (in a LCDM cosmology with H_0 = 70 km/s/Mpc, Omega_M = 0.3, and Omega _L=0.7), and either display at least one emission line with a full width at half maximum (FWHM)larger than 500 km/s or, if not, have interesting/complex absorption features. The catalog also includes previously known quasars (mostly from SDSS-I and II) that were reobserved by BOSS. The catalog contains 297,301 quasars detected over 9,376 square degrees with robust identification and redshift measured by a combination of principal component eigenspectra. The number of quasars with z>2.15 is about an order of magnitude greater than the number of z>2.15 quasars known prior to BOSS. Redshifts and FWHMs are provided for the strongest emission lines (CIV, CIII], MgII). The catalog identifies 29,580 broad absorption line quasars and lists their characteristics. For each object, the catalog presents five-band (u, g, r, i, z) CCD-based photometry together with some information on the optical morphology and the selection criteria. When available, the catalog also provides information on the optical variability of quasars using SDSS and PTF multi-epoch photometry. The catalog also contains X-ray, ultraviolet, near-infrared, and radio emission properties of the quasars, when available, from other large-area surveys. The calibrated digital spectra, covering the wavelength region 3,600-10,500A at a spectral resolution in the range 1,300<R<2,500, can be retrieved from the SDSS Catalog Archive Server.

We present the X-ray luminosity (L) versus dynamical mass (M) relation for 63 nearby clusters in the HIFLUGCS. The luminosity measurements are obtained based on ~1.3 Ms of clean XMM data and ROSAT pointed observations. The masses are estimated using optical spectroscopic redshifts of 13647 cluster galaxies in total. Given sufficient numbers of member galaxies in computing the dynamical masses, the L-M relations agree between the disturbed and undisturbed clusters. The cool-core clusters still dominate the scatter in the L-M relation even when a core corrected X-ray luminosity is used, which indicates that the scatter mainly reflects the structure formation history of the clusters. As shown by the clusters with a small number of redshifts, the dynamical masses can be underestimated leading to a biased scaling relation. To investigate the potential of spectroscopic surveys to follow up high-redshift galaxy clusters/groups observed in X-ray surveys for the identifications and mass calibrations, we carried out Monte-Carlo re-sampling of the cluster galaxy redshifts and calibrated the uncertainties of the redshift and dynamical mass estimates when only reduced numbers of galaxy redshifts per cluster are available. The re-sampling considers the SPIDERS and 4MOST configurations, designed for the follow-up of the eROSITA clusters, and was carried out for each cluster at the actual cluster redshift as well as at z=0.2, 0.4, 0.6, and 0.8. For following up very distant cluster/groups, we carried out the mass calibration based on the re-sampling with only 10zs/cluster, and redshift calibration based on the re-sampling with only 5zs/cluster and 10zs/cluster, respectively. Our results demonstrate the power of combining upcoming X-ray and optical spectroscopic surveys for mass calibration. The scatter in the dynamical mass estimates for the clusters with at least ten members is within 50%.

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

The uncertain origin of the recently-discovered `changing-looking' quasar phenomenon -- in which a luminous quasar dims significantly to a quiescent state in repeat spectroscopy over ~10 year timescales -- may present unexpected challenges to our understanding of quasar accretion. To better understand this phenomenon, we take a first step to building a sample of changing-look quasars with a systematic but simple archival search for these objects in the Sloan Digital Sky Survey Data Release 12. By leveraging the >10 year baselines for objects with repeat spectroscopy, we uncover two new changing-look quasars, and a third discovered previously. Decomposition of the multi-epoch spectra and analysis of the broad emission lines suggest that the quasar accretion disk emission dims due to rapidly decreasing accretion rates (by factors of >2.5), while disfavoring changes in intrinsic dust extinction for the two objects where these analyses are possible. Broad emission line energetics also support intrinsic dimming of quasar emission as the origin for this phenomenon rather than transient tidal disruption events or supernovae. Although our search criteria included quasars at all redshifts and transitions from either quasar-like to galaxy-like states or the reverse, all of the clear cases of changing-look quasars discovered were at relatively low-redshift (z ~ 0.2 - 0.3) and only exhibit quasar-like to galaxy-like transitions.

In this paper we describe and publicly release a catalogue consisting of 8445 point-like X-ray sources detected in the XMM-XXL north survey. For the 2512 AGN which have reliable spectroscopy, we present the X-ray spectral fitting which has been computed with a Bayesian approach. We have also applied an X-ray spectral stacking method to different sub-samples. We confirm the well-known Iwasawa-Taniguchi effect in our luminosity-redshift sub-samples, and argue that such an effect is due to a decrease in the covering factor of a distant obscuring torus with increasing X-ray luminosity. By comparing the distribution of the reflection fraction, the ratio of the normalization of the reflected component to the direct radiation, we find that the low-luminosity, low-redshift sub-sample had systematically higher reflection fraction values than the high-redshift, high-luminosity one. On the other hand, no significant difference is found between samples having similar luminosity but different redshift, suggesting that the structure of the torus does not evolve strongly with redshift. Contrary to previous works, we do not find evidence for an increasing photon index at high Eddington ratio. This may be an indication that the structure of the accretion disc changes as the Eddington ratio approaches unity. Comparing our X-ray spectral analysis results with the optical spectral classification, we find that $\sim20$ per cent of optical type-1 AGN show an X-ray absorbing column density higher than $10^{21.5}\,\rm{cm^{-2}}$, and about 50 per cent of type-2 AGN have an X-ray absorbing column density less than $10^{21.5}\,\rm{cm^{-2}}$. We suggest that the excess X-ray absorption shown in the high-luminosity optical type-1 AGN can be due to small-scale dust-free gas within (or close to) the broad line region, while in the low-luminosity ones it can be due to a clumpy torus with a large covering factor.

The northern tile of the wide-area and shallow XMM-XXL X-ray survey field is used to estimate the average dark matter halo mass of relatively luminous X-ray selected AGN [$\rm log\, L_X (\rm 2-10\,keV)= 43.6^{+0.4}_{-0.4}\,erg/s$] in the redshift interval $z=0.5-1.2$. Spectroscopic follow-up observations of X-ray sources in the XMM-XXL field by the Sloan telescope are combined with the VIPERS spectroscopic galaxy survey to determine the cross-correlation signal between X-ray selected AGN (total of 318) and galaxies (about 20,\,000). We model the large scales (2-25\,Mpc) of the correlation function to infer a mean dark matter halo mass of $\log M / (M_{\odot} \, h^{-1}) = 12.50 ^{+0.22} _{-0.30}$ for the X-ray selected AGN sample. This measurement is about 0.5\,dex lower compared to estimates in the literature of the mean dark matter halo masses of moderate luminosity X-ray AGN [$L_X (\rm 2-10\,keV)\approx 10^{42} - 10^{43}\,erg/s$] at similar redshifts. Our analysis also links the mean clustering properties of moderate luminosity AGN with those of powerful UV/optically selected QSOs, which are typically found in halos with masses few times $10^{12}\,M_{\odot}$. There is therefore evidence for a negative luminosity dependence of the AGN clustering. This is consistent with suggestions that AGN have a broad dark matter halo mass distribution with a high mass tail that becomes sub-dominant at high accretion luminosities. We further show that our results are in qualitative agreement with semi-analytic models of galaxy and AGN evolution, which attribute the wide range of dark matter halo masses among the AGN population to different triggering mechanisms and/or black hole fueling modes.

We investigate the role of supermassive black holes in the global context of galaxy evolution by measuring the host galaxy stellar mass function (HGMF) and the specific accretion rate i.e., lambda_SAR, distribution function (SARDF) up to z~2.5 with ~1000 X-ray selected AGN from XMM-COSMOS. Using a maximum likelihood approach, we jointly fit the stellar mass function and specific accretion rate distribution function, with the X-ray luminosity function as an additional constraint. Our best fit model characterizes the SARDF as a double power-law with mass dependent but redshift independent break whose low lambda_SAR slope flattens with increasing redshift while the normalization increases. This implies that, for a given stellar mass, higher lambda_SAR objects have a peak in their space density at earlier epoch compared to the lower lambda_SAR ones, following and mimicking the well known AGN cosmic downsizing as observed in the AGN luminosity function. The mass function of active galaxies is described by a Schechter function with a almost constant Mstar* and a low mass slope alpha that flattens with redshift. Compared to the stellar mass function, we find that the HGMF has a similar shape and that, up to log((Mstar/Msun)~11.5 the ratio of AGN host galaxies to star forming galaxies is basically constant (~10%). Finally, the comparison of the AGN HGMF for different luminosity and specific accretion rate sub-classes with the phenomenological model prediction by Peng et al. (2010) for the "transient" population, i.e. galaxies in the process of being mass-quenched, reveals that low-luminosity AGN do not appear to be able to contribute significantly to the quenching and that at least at high masses, i.e. Mstar>10^(10.7) Msun , feedback from luminous AGN (log(Lbol)>~46 [erg/s]) may be responsible for the quenching of star formation in the host galaxy.

The Extended Baryon Oscillation Spectroscopic Survey (eBOSS) will conduct novel cosmological observations using the BOSS spectrograph at Apache Point Observatory. Observations will be simultaneous with the Time Domain Spectroscopic Survey (TDSS) designed for variability studies and the Spectroscopic Identification of eROSITA Sources (SPIDERS) program designed for studies of X-ray sources. eBOSS will use four different tracers to measure the distance-redshift relation with baryon acoustic oscillations (BAO). Using more than 250,000 new, spectroscopically confirmed luminous red galaxies at a median redshift z=0.72, we project that eBOSS will yield measurements of $d_A(z)$ to an accuracy of 1.2% and measurements of H(z) to 2.1% when combined with the z>0.6 sample of BOSS galaxies. With ~195,000 new emission line galaxy redshifts, we expect BAO measurements of $d_A(z)$ to an accuracy of 3.1% and H(z) to 4.7% at an effective redshift of z= 0.87. A sample of more than 500,000 spectroscopically-confirmed quasars will provide the first BAO distance measurements over the redshift range 0.9<z<2.2, with expected precision of 2.8% and 4.2% on $d_A(z)$ and H(z), respectively. Finally, with 60,000 new quasars and re-observation of 60,000 quasars known from BOSS, we will obtain new Lyman-alpha forest measurements at redshifts z>2.1; these new data will enhance the precision of $d_A(z)$ and H(z) by a factor of 1.44 relative to BOSS. Furthermore, eBOSS will provide improved tests of General Relativity on cosmological scales through redshift-space distortion measurements, improved tests for non-Gaussianity in the primordial density field, and new constraints on the summed mass of all neutrino species. Here, we provide an overview of the cosmological goals, spectroscopic target sample, demonstration of spectral quality from early data, and projected cosmological constraints from eBOSS.

Capitalizing on the all-sky coverage of {\it WISE}, and the 35\% and 50\% sky coverage from SDSS and Pan-STARRS, respectively, we explore the efficacy of $m_{R}$ (optical) - $m_{3.4 \mu m}$ (mid-infrared), hereafter $R-W1$, as a color diagnostic to identify obscured supermassive black hole accretion in wide-area X-ray surveys. We use the $\sim$16.5 deg$^2$ Stripe 82 X-ray survey data as a test-bed to compare $R-W1$ with $R-K$, an oft-used obscured AGN selection criterion, and examine where different classes of objects lie in this parameter space. Most stars follow a well-defined path in $R-K$ vs. $R-W1$ space. We demonstrate that optically normal galaxies hosting X-ray AGN at redshifts $0.5<z<1$ can be recovered with an $R-W1>4$ color-cut, while they typically are not selected as AGN based on their $W1-W2$ colors. Additionally, different observed X-ray luminosity bins favor different regions in $R-W1$ parameter space: moderate luminosity AGN ($10^{43}$ erg s$^{-1} < L_{\rm 0.5-10 keV} < 10^{44}$ erg s$^{-1}$) tend to have red colors while the highest luminosity AGN ($L_{\rm 0.5-10 keV} > 10^{45}$ erg s$^{-1}$) have bluer colors; higher spectroscopic completeness of the Stripe 82X sample is needed to determine whether this is a selection effect or an intrinsic property. Finally, we parameterize X-ray obscuration of Stripe 82X AGN by calculating their hardness ratios (HRs) and find no clear trends between HR and optical reddening. Our results will help inform best-effort practices in following-up obscured AGN candidates in current and future wide-area, shallow X-ray surveys, including the all-sky {\it eROSITA} mission.

This paper presents a survey of X-ray selected active galactic nuclei (AGN) with optical spectroscopic follow-up in a $\sim 18\, \rm{deg^2}$ area of the equatorial XMM-XXL north field. A sample of 8445 point-like X-ray sources detected by XMM-Newton above a limiting flux of $F_{\rm 0.5-10\, keV} > 10^{-15} \rm\,erg\, cm^{-2}\, s^{-1}$ was matched to optical (SDSS) and infrared (WISE) counterparts. We followed up 3042 sources brighter than $r=22.5$ mag with the SDSS BOSS spectrograph. The spectra yielded a reliable redshift measurement for 2578 AGN in the redshift range $z=0.02-5.0$, with $0.5-2\rm\, keV$ luminosities ranging from $10^{39}-10^{46}\rm\,erg\,s^{-1}$. This is currently the largest published spectroscopic sample of X-ray selected AGN in a contiguous area. The BOSS spectra of AGN candidates show a bimodal distribution of optical line widths allowing a separation between broad- and narrow-emission line AGN. The former dominate our sample (70 per cent) due to the relatively bright X-ray flux limit and the optical BOSS magnitude limit. We classify the narrow emission line objects (22 per cent of full sample) using standard BPT diagnostics: the majority have line ratios indicating the dominant source of ionization is the AGN. A small number (8 per cent of full sample) exhibit the typical narrow line ratios of star-forming galaxies, or only have absorption lines in their spectra. We term the latter two classes "elusive'' AGN. We also compare X-ray, optical and infrared color AGN selections in this field. X-ray observations reveal, the largest number of AGN. The overlap between the selections, which is a strong function of the imaging depth in a given band, is also remarkably small. We show using spectral stacking that a large fraction of the X-ray AGN would not be selectable via optical or IR colours due to host galaxy contamination.

The advent of the new generation of X-ray telescopes yielded a significant step forward in our understanding of ionised absorption generated in the accretion discs of X-ray binaries. It has become evident that these relatively weak and narrow absorption features, sporadically present in the X-ray spectra of some systems, are actually the signature of equatorial outflows, which might carry away more matter than that being accreted. Therefore, they play a major role in the accretion phenomenon. These outflows (or ionised atmospheres) are ubiquitous during the softer states but absent during the power-law dominated, hard states, suggesting a strong link with the state of the inner accretion disc, presence of the radio-jet and the properties of the central source. Here, we discuss the current understanding of this field.

We report the discovery of a new "changing-look" quasar, SDSS J101152.98+544206.4, through repeat spectroscopy from the Time Domain Spectroscopic Survey. This is an addition to a small but growing set of quasars whose blue continua and broad optical emission lines have been observed to decline by a large factor on a time scale of approximately a decade. The 5100 Angstrom monochromatic continuum luminosity of this quasar drops by a factor of > 9.8 in a rest-frame time interval of < 9.7 years, while the broad H-alpha luminosity drops by a factor of 55 in the same amount of time. The width of the broad H-alpha line increases in the dim state such that the black hole mass derived from the appropriate single-epoch scaling relation agrees between the two epochs within a factor of 3. The fluxes of the narrow emission lines do not appear to change between epochs. The light curve obtained by the Catalina Sky Survey suggests that the transition occurs within a rest-frame time interval of approximately 500 days. We examine three possible mechanisms for this transition suggested in the recent literature. An abrupt change in the reddening towards the central engine is disfavored by the substantial difference between the timescale to obscure the central engine and the observed timescale of the transition. A decaying tidal disruption flare is consistent with the decay rate of the light curve but not with the prolonged bright state preceding the decay, nor can this scenario provide the power required by the luminosities of the emission lines. An abrupt drop in the accretion rate onto the supermassive black hole appears to be the most plausible explanation for the rapid dimming.

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

We present a study of the X-ray flaring activity of Sgr A* during all the 150 XMM-Newton and Chandra observations pointed at the Milky Way center over the last 15 years. This includes the latest XMM-Newton and Chandra campaigns devoted to monitoring the closest approach of the very red Br-Gamma emitting object called G2. The entire dataset analysed extends from September 1999 through November 2014. We employed a Bayesian block analysis to investigate any possible variations in the characteristics (frequency, energetics, peak intensity, duration) of the flaring events that Sgr A* has exhibited since their discovery in 2001. We observe that the total bright-or-very bright flare luminosity of Sgr A* increased between 2013-2014 by a factor of 2-3 (~3.5 sigma significance). We also observe an increase (~99.9% significance) from 0.27+-0.04 to 2.5+-1.0 day^-1 of the bright-or-very bright flaring rate of Sgr A*, starting in late summer 2014, which happens to be about six months after G2's peri-center passage. This might indicate that clustering is a general property of bright flares and that it is associated with a stationary noise process producing flares not uniformly distributed in time (similar to what is observed in other quiescent black holes). If so, the variation in flaring properties would be revealed only now because of the increased monitoring frequency. Alternatively, this may be the first sign of an excess accretion activity induced by the close passage of G2. More observations are necessary to distinguish between these two hypotheses.

The third generation of the Sloan Digital Sky Survey (SDSS-III) took data from 2008 to 2014 using the original SDSS wide-field imager, the original and an upgraded multi-object fiber-fed optical spectrograph, a new near-infrared high-resolution spectrograph, and a novel optical interferometer. All the data from SDSS-III are now made public. In particular, this paper describes Data Release 11 (DR11) including all data acquired through 2013 July, and Data Release 12 (DR12) adding data acquired through 2014 July (including all data included in previous data releases), marking the end of SDSS-III observing. Relative to our previous public release (DR10), DR12 adds one million new spectra of galaxies and quasars from the Baryon Oscillation Spectroscopic Survey (BOSS) over an additional 3000 sq. deg of sky, more than triples the number of H-band spectra of stars as part of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE), and includes repeated accurate radial velocity measurements of 5500 stars from the Multi-Object APO Radial Velocity Exoplanet Large-area Survey (MARVELS). The APOGEE outputs now include measured abundances of 15 different elements for each star. In total, SDSS-III added 2350 sq. deg of ugriz imaging; 155,520 spectra of 138,099 stars as part of the Sloan Exploration of Galactic Understanding and Evolution 2 (SEGUE-2) survey; 2,497,484 BOSS spectra of 1,372,737 galaxies, 294,512 quasars, and 247,216 stars over 9376 sq. deg; 618,080 APOGEE spectra of 156,593 stars; and 197,040 MARVELS spectra of 5,513 stars. Since its first light in 1998, SDSS has imaged over 1/3 of the Celestial sphere in five bands and obtained over five million astronomical spectra.

In the last decade, a combination of high sensitivity, high spatial resolution observations and of coordinated multi-wavelength surveys has revolutionized our view of extra-galactic black hole (BH) astrophysics. We now know that supermassive black holes reside in the nuclei of almost every galaxy, grow over cosmological times by accreting matter, interact and merge with each other, and in the process liberate enormous amounts of energy that influence dramatically the evolution of the surrounding gas and stars, providing a powerful self-regulatory mechanism for galaxy formation. The different energetic phenomena associated to growing black holes and Active Galactic Nuclei (AGN), their cosmological evolution and the observational techniques used to unveil them, are the subject of this chapter. In particular, I will focus my attention on the connection between the theory of high-energy astrophysical processes giving rise to the observed emission in AGN, the observable imprints they leave at different wavelengths, and the methods used to uncover them in a statistically robust way. I will show how such a combined effort of theorists and observers have led us to unveil most of the SMBH growth over a large fraction of the age of the Universe, but that nagging uncertainties remain, preventing us from fully understating the exact role of black holes in the complex process of galaxy and large-scale structure formation, assembly and evolution.

We argue that the `changing look' AGN recently reported by LaMassa et al. could be a luminous flare produced by the tidal disruption of a super-solar mass star passing just a few gravitational radii outside the event horizon of a $\sim 10^8 M_{\odot}$ nuclear black hole. This flare occurred in a massive, star forming galaxy at redshift $z=0.312$, robustly characterized thanks to repeated late-time photometric and spectroscopic observations. By taking difference-photometry of the well sampled multi-year SDSS Stripe-82 light-curve, we are able to probe the evolution of the nuclear spectrum over the course of the outburst. The tidal disruption event (TDE) interpretation is consistent with the very rapid rise and the decay time of the flare, which displays an evolution consistent with the well-known $t^{-5/3}$ behaviour (with a clear superimposed re-brightening flare). Our analysis places constraints on the physical properties of the TDE, such as the putative disrupted star's mass and orbital parameters, as well as the size and temperature of the emitting material. The properties of the broad and narrow emission lines observed in two epochs of SDSS spectra provide further constraints on the circum-nuclear structure, and could be indicative that the system hosted a moderate-luminosity AGN as recently as a few $10^4$ years ago, and is likely undergoing residual accretion as late as ten years after peak, as seen from the broad H$\alpha$ emission line. We discuss the complex interplay between tidal disruption events and gas accretion episodes in galactic nuclei, highlighting the implications for future TDE searches and for estimates of their intrinsic rates.

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

We present host stellar velocity dispersion measurements for a sample of 88 broad-line quasars at 0.1<z<1 (46 at z>0.6) from the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project. High signal-to-noise ratio coadded spectra (average S/N~30 per 69 km/s pixel) from SDSS-RM allowed decomposition of the host and quasar spectra, and measurement of the host stellar velocity dispersions and black hole (BH) masses using the single-epoch (SE) virial method. The large sample size and dynamic range in luminosity (L5100=10^(43.2-44.7) erg/s) lead to the first clear detection of a correlation between SE virial BH mass and host stellar velocity dispersion far beyond the local universe. However, the observed correlation is significantly flatter than the local relation, suggesting that there are selection biases in high-z luminosity-threshold quasar samples for such studies. Our uniform sample and analysis enable an investigation of the redshift evolution of the M-sigma relation free of caveats by comparing different samples/analyses at disjoint redshifts. We do not observe evolution of the M-sigma relation in our sample, up to z~1, but there is an indication that the relation flattens towards higher redshifts. Coupled with the increasing threshold luminosity with redshift in our sample, this again suggests certain selection biases are at work, and simple simulations demonstrate that a constant M-sigma relation is favored to z~1. Our results highlight the scientific potential of deep coadded spectroscopy from quasar monitoring programs, and offer a new path to probe the co-evolution of BHs and galaxies at earlier times.

We construct the rest-frame 2--10 keV intrinsic X-ray luminosity function of Active Galactic Nuclei (AGNs) from a combination of X-ray surveys from the all-sky Swift BAT survey to the Chandra Deep Field-South. We use ~3200 AGNs in our analysis, which covers six orders of magnitude in flux. The inclusion of the XMM and Chandra COSMOS data has allowed us to investigate the detailed behavior of the XLF and evolution. In deriving our XLF, we take into account realistic AGN spectrum templates, absorption corrections, and probability density distributions in photometric redshift. We present an analytical expression for the overall behavior of the XLF in terms of the luminosity-dependent density evolution, smoothed two power-law expressions in 11 redshift shells, three-segment power-law expression of the number density evolution in four luminosity classes, and binned XLF. We observe a sudden flattening of the low luminosity end slope of the XLF slope at z>~0.6. Detailed structures of the AGN downsizing have been also revealed, where the number density curves have two clear breaks at all luminosity classes above log LX>43. The two break structure is suggestive of two-phase AGN evolution, consisting of major merger triggering and secular processes.

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

Massive black hole binaries are naturally predicted in the context of the hierarchical model of structure formation. The binaries that manage to lose most of their angular momentum can coalesce to form a single remnant. In the last stages of this process, the holes undergo an extremely loud phase of gravitational wave emission, possibly detectable by current and future probes. The theoretical effort towards obtaining a coherent physical picture of the binary path down to coalescence is still underway. In this paper, for the first time, we take advantage of observational studies of active galactic nuclei evolution to constrain the efficiency of gas-driven binary decay. Under conservative assumptions we find that gas accretion toward the nuclear black holes can efficiently lead binaries of any mass forming at high redshift (> 2) to coalescence within the current time. The observed "downsizing" trend of the accreting black hole luminosity function further implies that the gas inflow is sufficient to drive light black holes down to coalescence, even if they bind in binaries at lower redshifts, down to z~0.5 for binaries of ~10 million solar masses, and z~0.2 for binaries of ~1 million solar masses. This has strong implications for the detection rates of coalescing black hole binaries of future space-based gravitational wave experiments.

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

Heavily obscured, Compton Thick (CT, NH>10^24 cm^-2) AGN may represent an important phase in AGN/galaxy co-evolution and are expected to provide a significant contribution to the cosmic X-ray background (CXB). Through direct X-ray spectra analysis, we selected 39 heavily obscured AGN (NH>3x10^23 cm^-2) in the 2 deg^2 XMM-COSMOS survey. After selecting CT AGN based on the fit of a simple absorbed two power law model to the XMM data, the presence of CT AGN was confirmed in 80% of the sources using deeper Chandra data and more complex models. The final sample of CT AGN comprises 10 sources spanning a large range of redshift and luminosity. We collected the multi-wavelength information available for all these sources, in order to study the distribution of SMBH and host properties, such as BH mass (M_BH), Eddington ratio (\lambda_Edd), stellar mass (M*), specific star formation rate (sSFR) in comparison with a sample of unobscured AGN. We find that highly obscured sources tend to have significantly smaller M_BH and higher \lambda_edd with respect to unobscured ones, while a weaker evolution in M* is observed. The sSFR of highly obscured sources is consistent with the one observed in the main sequence of star forming galaxies, at all redshift. We also present optical spectra, spectral energy distribution (SED) and morphology for the sample of 10 CT AGN: all the available optical spectra are dominated by the stellar component of the host galaxy, and a highly obscured torus component is needed in the SED of the CT sources. Exploiting the high resolution Hubble-ACS images available, we conclude that these highly obscured sources have a significantly larger merger fraction with respect to other X-ray selected samples of AGN. Finally we discuss implications in the context of AGN/galaxy co-evolutionary models, and compare our results with the predictions of CXB synthesis models.

Quasar feedback in the form of powerful outflows is invoked as a key mechanism to quench star formation in galaxies, preventing massive galaxies to over-grow and producing the red colors of ellipticals. On the other hand, some models are also requiring `positive' AGN feedback, inducing star formation in the host galaxy through enhanced gas pressure in the interstellar medium. However, finding observational evidence of the effects of both types of feedback is still one of the main challenges of extragalactic astronomy, as few observations of energetic and extended radiatively-driven winds are available. Here we present SINFONI near infrared integral field spectroscopy of XID2028, an obscured, radio-quiet z=1.59 QSO detected in the XMM-COSMOS survey, in which we clearly resolve a fast (1500 km/s) and extended (up to 13 kpc from the black hole) outflow in the [OIII] lines emitting gas, whose large velocity and outflow rate are not sustainable by star formation only. The narrow component of Ha emission and the rest frame U band flux from HST-ACS imaging enable to map the current star formation in the host galaxy: both tracers independently show that the outflow position lies in the center of an empty cavity surrounded by star forming regions on its edge. The outflow is therefore removing the gas from the host galaxy (`negative feedback'), but also triggering star formation by outflow induced pressure at the edges (`positive feedback'). XID2028 represents the first example of a host galaxy showing both types of feedback simultaneously at work.

We present X-shooter at the Vewry Large Telescope observations of a sample of 10 luminous, X-ray obscured QSOs at z$\sim1.5$ from the XMM-COSMOS survey, expected to be caught in the transitioning phase from starburst to active galactic nucleus (AGN)-dominated systems. The main selection criterion is X-ray detection at bright fluxes (L$>=10^{44}$ erg s$^{-1}$) coupled to red optical-to-NIR-to-MIR colors. Thanks to its large wavelength coverage, X-shooter allowed us to determine accurate redshifts from the presence of multiple emission lines for five out of six targets for which we had only a photometric redshift estimate, with an 80 percent success rate, significantly larger than what is observed in similar programs of spectroscopic follow-up of red QSOs. We report the detection of broad and shifted components in the [OIII]$\lambda\lambda$5007,4959 complexes for six out of eight sources with these lines observable in regions free from strong atmospheric absorptions. The full width half-maximum (FWHM) associated with the broad components are in the range FWHM$\sim900-1600$ km s$^{-1}$, larger than the average value observed in SDSS Type 2 AGN samples at similar observed [OIII] luminosity, but comparable to those observed for QSO/ultraluminous infrared galaxies systems for which the presence of kpc scale outflows have been revealed through integral field unit spectroscopy. Although the total outflow energetics (inferred under reasonable assumptions) may be consistent with winds accelerated by stellar processes, we favour an AGN origin for the outflows given the high outflow velocities oberved (v$>1000$ km s$^{-1}$) and the presence of strong winds also in objects undetected in the far infrared.

We present LOFAR Low Band observations of the Bootes and 3C295 fields. Our images made at 34, 46, and 62 MHz reach noise levels of 12, 8, and 5 mJy beam$^{-1}$, making them the deepest images ever obtained in this frequency range. In total, we detect between 300 and 400 sources in each of these images, covering an area of 17 to 52 deg$^{2}$. From the observations we derive Euclidean-normalized differential source counts. The 62 MHz source counts agree with previous GMRT 153 MHz and VLA 74 MHz differential source counts, scaling with a spectral index of $-0.7$. We find that a spectral index scaling of $-0.5$ is required to match up the LOFAR 34 MHz source counts. This result is also in agreement with source counts from the 38 MHz 8C survey, indicating that the average spectral index of radio sources flattens towards lower frequencies. We also find evidence for spectral flattening using the individual flux measurements of sources between 34 and 1400 MHz and by calculating the spectral index averaged over the source population. To select ultra-steep spectrum ($\alpha < -1.1$) radio sources, that could be associated with massive high redshift radio galaxies, we compute spectral indices between 62 MHz, 153 MHz and 1.4 GHz for sources in the Bo\"otes field. We cross-correlate these radio sources with optical and infrared catalogues and fit the spectral energy distribution to obtain photometric redshifts. We find that most of these ultra-steep spectrum sources are located in the $ 0.7 \lesssim z \lesssim 2.5$ range.

We present near-infrared spectra, obtained with SINFONI and XShooter observations at ESO VLT, of nine dusty, red QSOs at 1.2<z<2.6. The sources are hard X-ray detected, characterized by cold absorption (N_H>10^{21} - 10^{22} cm^{-2}) and show a broad Ha component in the NIR spectra. We complement this sample with twelve additional sources taken from the literature with similar properties resulting in a total sample of 21 X-ray obscured, intermediate type (1.8-1.9), dusty reddened QSOs. From the broad Ha line we have computed the BH masses through the virial formula and derived Eddington ratios. Moreover, from optical/IR multi-component SED fitting we have derived the stellar mass of their host galaxies and their SFRs. We find that most of the sources in our sample are hosted in starburst and main sequence star-forming galaxies with Eddington ratios lambda>0.1. We find a strong trend with the BH mass i.e. less massive objects are scattered below and above the local relation while the most massive ones are mainly located above it. We also studied the evolution of these sources on the M_BH-M_star plane compared to a sample of optically blue type--1 QSOs and we find that obscured red QSOs show a ratio of M_BH to M_star that increases with redshift which is consistent with or slightly lower than what has been found for blue QSOs. These sources may represent the blow-out phase at the end of the rapid BH growth and immediately preceding the classical blue QSOs typically sampled in optical surveys. They in fact show evidence of outflows in the ionized gas component, but their BH has already fully formed.

AGN are known to have complex X-ray spectra that depend on both the properties of the accreting SMBH (e.g. mass, accretion rate) and the distribution of obscuring material in its vicinity ("torus"). Often however, simple and even unphysical models are adopted to represent the X-ray spectra of AGN. In the case of blank field surveys in particular, this should have an impact on e.g. the determination of the AGN luminosity function, the inferred accretion history of the Universe and also on our understanding of the relation between AGN and their host galaxies. We develop a Bayesian framework for model comparison and parameter estimation of X-ray spectra. We take into account uncertainties associated with X-ray data and photometric redshifts. We also demonstrate how Bayesian model comparison can be used to select among ten different physically motivated X-ray spectral models the one that provides a better representation of the observations. Despite the use of low-count spectra, our methodology is able to draw strong inferences on the geometry of the torus. For a sample of 350 AGN in the 4 Ms Chandra Deep Field South field, our analysis identifies four components needed to represent the diversity of the observed X-ray spectra: (abridged). Simpler models are ruled out with decisive evidence in favour of a geometrically extended structure with significant Compton scattering. Regarding the geometry of the obscurer, there is strong evidence against both a completely closed or entirely open toroidal geometry, in favour of an intermediate case. The additional Compton reflection required by data over that predicted by toroidal geometry models, may be a sign of a density gradient in the torus or reflection off the accretion disk. Finally, we release a catalogue with estimated parameters such as the accretion luminosity in the 2-10 keV band and the column density, $N_{H}$, of the obscurer.

We took advantage of the observations carried out by XMM in the COSMOS field during 3.5 years, to study the long term variability of a large sample of AGN (638 sources), in a wide range of redshift (0.1<z<3.5) and X-ray luminosity ($10^{41}<$L(2-10)$<10^{45.5}$). Both a simple statistical method to asses the significance of variability, and the Normalized Excess Variance ($\sigma^{2}_{rms}$) parameter, where used to obtain a quantitative measurement of the variability. Variability is found to be prevalent in most AGN, whenever we have good statistic to measure it, and no significant differences between type-1 and type-2 AGN were found. A flat (slope -0.23+/-0.03) anti-correlation between $\sigma^{2}_{rms}$ and X-ray luminosity is found, when significantly variable sources are considered all together. When divided in three redshift bins, the anti-correlation becomes stronger and evolving with z, with higher redshift AGN being more variable. We prove however that this effect is due to the pre-selection of variable sources: considering all the sources with available $\sigma^{2}_{rms}$ measurement, the evolution in redshift disappears. For the first time we were also able to study the long term X-ray variability as a function of $M_{\rm BH}$ and Eddington ratio, for a large sample of AGN spanning a wide range of redshift. An anti-correlation between $\sigma^{2}_{rms}$ and $M_{\rm BH}$ is found, with the same slope of the anti-correlation between $\sigma^{2}_{rms}$ and X-ray luminosity, suggesting that the latter can be a byproduct of the former one. No clear correlation is found between $\sigma^{2}_{rms}$ and the Eddington ratio in our sample. Finally, no correlation is found between the X-ray $\sigma^{2}_{rms}$ and the optical variability.

The mid-infrared to ultraviolet (0.1 -- 10 $\mu m$) spectral energy distribution (SED) shapes of 407 X-ray-selected radio-quiet type 1 AGN in the wide-field ``Cosmic Evolution Survey" (COSMOS) have been studied for signs of evolution. For a sub-sample of 200 radio-quiet quasars with black hole mass estimates and host galaxy corrections, we studied their mean SEDs as a function of a broad range of redshift, bolometric luminosity, black hole mass and Eddington ratio, and compared them with the Elvis et al. (1994, E94) type 1 AGN mean SED. We found that the mean SEDs in each bin are closely similar to each other, showing no statistical significant evidence of dependence on any of the analyzed parameters. We also measured the SED dispersion as a function of these four parameters, and found no significant dependencies. The dispersion of the XMM-COSMOS SEDs is generally larger than E94 SED dispersion in the ultraviolet, which might be due to the broader ``window function'' for COSMOS quasars, and their X-ray based selection.

Accreting black holes on all mass scales (from stellar to supermassive) appear to follow a nonlinear relation between X-ray luminosity, radio luminosity and BH mass, indicating that similar physical processes drive the central engines in X-ray binaries and active galactic nuclei (AGN). However, in recent years an increasing number of BH systems have been identified that do not fit into this scheme. These outliers may be the key to understand how BH systems are powered by accretion. Here we present results from EVN observations of a sample of low-luminosity AGN (LLAGN) with known mass that have unusually high radio powers when compared with their X-ray luminosity.

We study the incidence of nuclear obscuration on a complete sample of 1310 AGN selected on the basis of their rest-frame 2-10 keV X-ray flux from the XMM-COSMOS survey, in the redshift range 0.3<z<3.5. We classify the AGN as obscured or un-obscured on the basis of either the optical spectral properties and the overall SED or the shape of the X-ray spectrum. The two classifications agree in about 70% of the objects, and the remaining 30% can be further subdivided into two distinct classes: at low luminosities X-ray un-obscured AGN do not always show signs of broad lines or blue/UV continuum emission in their optical spectra, most likely due to galaxy dilution effects; at high luminosities broad line AGN may have absorbed X-ray spectra, which hints at an increased incidence of small-scale (sub-parsec) dust-free obscuration. We confirm that the fraction of obscured AGN is a decreasing function of the intrinsic X-ray luminosity, while the incidence of absorption shows significant evolution only for the most luminous AGN, which appear to be more commonly obscured at higher redshift. We find no significant difference between the mean stellar masses and star formation rates of obscured and un-obscured AGN hosts. We conclude that the physical state of the medium responsible for obscuration in AGN is complex, and mainly determined by the radiation environment (nuclear luminosity) in a small region enclosed within the gravitational sphere of influence of the central black hole, but is largely insensitive to the wider scale galactic conditions.

Quasi-stellar objects (QSOs) occur in galaxies in which supermassive black holes (SMBHs) are growing substantially through rapid accretion of gas. Many popular models of the co-evolutionary growth of galaxies and SMBHs predict that QSOs are also sites of substantial recent star formation, mediated by important processes, such as major mergers, which rapidly transform the nature of galaxies. A detailed study of the star-forming properties of QSOs is a critical test of such models. We present a far-infrared Herschel/PACS study of the mean star formation rate (SFR) of a sample of spectroscopically observed QSOs to z~2 from the COSMOS extragalactic survey. This is the largest sample to date of moderately luminous AGNs studied using uniform, deep far-infrared photometry. We study trends of the mean SFR with redshift, black hole mass, nuclear bolometric luminosity and specific accretion rate (Eddington ratio). To minimize systematics, we have undertaken a uniform determination of SMBH properties, as well as an analysis of important selection effects within spectroscopic QSO samples that influence the interpretation of SFR trends. We find that the mean SFRs of these QSOs are consistent with those of normal massive star-forming galaxies with a fixed scaling between SMBH and galaxy mass at all redshifts. No strong enhancement in SFR is found even among the most rapidly accreting systems, at odds with several co-evolutionary models. Finally, we consider the qualitative effects on mean SFR trends from different assumptions about the star-forming properties of QSO hosts and redshift evolution of the SMBH-galaxy relationship. While limited currently by uncertainties, valuable constraints on AGN-galaxy co-evolution can emerge from our approach.

Massive black holes (MBHs) in contrast to stellar mass black holes are expected to substantially change their properties over their lifetime. MBH masses increase by several order of magnitude over the Hubble time, as illustrated by Soltan's argument. MBH spins also must evolve through the series of accretion and mergers events that grow the MBH's masses. We present a simple model that traces the joint evolution of MBH masses and spins across cosmic time. Our model includes MBH-MBH mergers, merger-driven gas accretion, stochastic fueling of MBHs through molecular cloud capture, and a basic implementation of accretion of recycled gas. This approach aims at improving the modeling of low-redshift MBHs and AGN, whose properties can be more easily estimated observationally. Despite the simplicity of the model, it captures well the global evolution of the MBH population from z\sim6 to today. Under our assumptions, we find that the typical spin and radiative efficiency of MBHs decrease with cosmic time because of the higher incidence of stochastic processes in gas-rich galaxies and MBH-MBH mergers in gas-poor galaxies. At z=0 the spin distribution in gas-poor galaxies peaks at spins 0.4-0.8, and it is not strongly mass dependent. MBHs in gas-rich galaxies have a more complex evolution, with low-mass MBHs at low redshift having low spins, and spins increasing at larger masses and redshifts. We also find that at z>1 MBH spins are on average highest in high luminosity AGN, while at lower redshifts these differences disappear.

We define a quasar-galaxy mixing diagram using the slopes of their spectral energy distributions (SEDs) from 1\mu m to 3000\AA\ and from 1\mu m to 3\mu m in the rest frame. The mixing diagram can easily distinguish among quasar-dominated, galaxy-dominated and reddening-dominated SED shapes. By studying the position of the 413 XMM selected Type 1 AGN in the wide-field "Cosmic Evolution Survey" (COSMOS) in the mixing diagram, we find that a combination of the Elvis et al. (1994, hereafter E94) quasar SED with various contributions from galaxy emission and some dust reddening is remarkably effective in describing the SED shape from 0.3-3\mu m for large ranges of redshift, luminosity, black hole mass and Eddington ratio of type 1 AGN. In particular, the location in the mixing diagram of the highest luminosity AGN is very close (within 1\sigma) to that of the E94 SED. The mixing diagram can also be used to estimate the host galaxy fraction and reddening in quasar. We also show examples of some outliers which might be AGN in different evolutionary stages compared to the majority of AGN in the quasar-host galaxy co-evolution cycle.

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.

We present the X-ray spectral analysis of the 390 brightest extragalactic sources in the Chandra-COSMOS catalog, showing at least 70 net counts in the 0.5-7 keV band. This sample has a 100% completeness in optical-IR identification, with 75% of the sample having a spectroscopic redshift and 25% a photometric redshift. Our analysis allows us to accurately determine the intrinsic absorption, the broad band continuum shape ({\Gamma}) and intrinsic L(2-10) distributions, with an accuracy better than 30% on the spectral parameters for 95% of the sample. The sample is equally divided in type-1 (49.7%) and type-2 AGN (48.7%) plus few passive galaxies at low z. We found a significant difference in the distribution of {\Gamma} of type-1 and type-2, with small intrinsic dispersion, a weak correlation of {\Gamma} with L(2-10) and a large population (15% of the sample) of high luminosity, highly obscured (QSO2) sources. The distribution of the X ray/Optical flux ratio (Log(FX /Fi)) for type-1 is narrow (0 < X/O < 1), while type-2 are spread up to X/O = 2. The X/O correlates well with the amount of X-ray obscuration. Finally, a small sample of Compton thick candidates and peculiar sources is presented. In the appendix we discuss the comparison between Chandra and XMM-Newton spectra for 280 sources in common. We found a small systematic difference, with XMM-Newton spectra that tend to have softer power-laws and lower obscuration.

The Square Kilometre Array (SKA) is the most ambitious radio telescope ever planned. With a collecting area of about a square kilometre, the SKA will be far superior in sensitivity and observing speed to all current radio facilities. The scientific capability promised by the SKA and its technological challenges provide an ideal base for interdisciplinary research, technology transfer, and collaboration between universities, research centres and industry. The SKA in the radio regime and the European Extreme Large Telescope (E-ELT) in the optical band are on the roadmap of the European Strategy Forum for Research Infrastructures (ESFRI) and have been recognised as the essential facilities for European research in astronomy. This "White Paper" outlines the German science and R&D interests in the SKA project and will provide the basis for future funding applications to secure German involvement in the Square Kilometre Array.

M87 is a giant elliptical galaxy located in the centre of the Virgo cluster, which harbours a supermassive black hole of mass 6.4x10^9 M_sun, whose activity is responsible for the extended (80 kpc) radio lobes that surround the galaxy. The energy generated by matter falling onto the central black hole is ejected and transferred to the intra-cluster medium via a relativistic jet and morphologically complex systems of buoyant bubbles, which rise towards the edges of the extended halo. Here we present the first observations made with the new Low-Frequency Array (LOFAR) of M87 at frequencies down to 20 MHz. Images of M87 were produced at low radio frequencies never explored before at these high spatial resolution and dynamic range. To disentangle different synchrotron models and place constraints on source magnetic field, age and energetics, we also performed a detailed spectral analysis of M87 extended radio-halo using these observations together with archival data. We do not find any sign of new extended emissions; on the contrary the source appears well confined by the high pressure of the intra-cluster medium. A continuous injection of relativistic electrons is the model that best fits our data, and provides a scenario in which the lobes are still supplied by fresh relativistic particles from the active galactic nuclei. We suggest that the discrepancy between the low-frequency radio-spectral slope in the core and in the halo implies a strong adiabatic expansion of the plasma as soon as it leaves the core area. The extended halo has an equipartition magnetic field strength of ~10 uG, which increases to ~13 uG in the zones where the particle flows are more active. The continuous injection model for synchrotron ageing provides an age for the halo of ~40 Myr, which in turn provides a jet kinetic power of 6-10x10^44 erg/s.

eROSITA is the primary instrument on the Russian SRG mission. In the first four years of scientific operation after its launch, foreseen for 2014, it will perform a deep survey of the entire X-ray sky. In the soft X-ray band (0.5-2 keV), this will be about 20 times more sensitive than the ROSAT all sky survey, while in the hard band (2-10 keV) it will provide the first ever true imaging survey of the sky at those energies. Such a sensitive all-sky survey will revolutionize our view of the high-energy sky, and calls for major efforts in synergic, multi-wavelength wide area surveys in order to fully exploit the scientific potential of the X-ray data. The design-driving science of eROSITA is the detection of very large samples (~10^5 objects) of galaxy clusters out to redshifts z>1, in order to study the large scale structure in the Universe, test and characterize cosmological models including Dark Energy. eROSITA is also expected to yield a sample of around 3 millions Active Galactic Nuclei, including both obscured and un-obscured objects, providing a unique view of the evolution of supermassive black holes within the emerging cosmic structure. The survey will also provide new insights into a wide range of astrophysical phenomena, including accreting binaries, active stars and diffuse emission within the Galaxy, as well as solar system bodies that emit X-rays via the charge exchange process. Finally, such a deep imaging survey at high spectral resolution, with its scanning strategy sensitive to a range of variability timescales from tens of seconds to years, will undoubtedly open up a vast discovery space for the study of rare, unpredicted, or unpredictable high-energy astrophysical phenomena. In this living document we present a comprehensive description of the main scientific goals of the mission, with strong emphasis on the early survey phases.

The "Cosmic Evolution Survey" (COSMOS) enables the study of the Spectral Energy Distributions (SEDs) of Active Galactic Nuclei (AGN) because of the deep coverage and rich sampling of frequencies from X-ray to radio. Here we present a SED catalog of 413 X-ray (\xmm) selected type 1 (emission line FWHM$>2000$ km s$^{-1}$) AGN with Magellan, SDSS or VLT spectrum. The SEDs are corrected for the Galactic extinction, for broad emission line contributions, constrained variability, and for host galaxy contribution. We present the mean SED and the dispersion SEDs after the above corrections in the rest frame 1.4 GHz to 40 keV, and show examples of the variety of SEDs encountered. In the near-infrared to optical (rest frame $\sim 8\mu m$-- 4000\AA), the photometry is complete for the whole sample and the mean SED is derived from detections only. Reddening and host galaxy contamination could account for a large fraction of the observed SED variety. The SEDs are all available on-line.

Using the wide multi-band photometry available in the COSMOS field we explore the host galaxy properties of a large sample of Active Galactic Nuclei (AGN) obtained by combining X-ray and optical spectroscopic selections. Based on a careful study of their Spectral Energy Distribution (SED), which has been parametrized using a 2-component (AGN+galaxy) model fit, we derived dust-corrected rest-frame magnitudes, colors, stellar masses and star formation rates (SFRs). We find that AGN hosts span a large range of stellar masses and SFRs. No color-bimodality is seen at any redshift in the AGN hosts, which are found to be mainly massive, red galaxies. Once accounting for the color-mass degeneracy in well defined mass-matched samples, we find a residual marginal enhancement of AGN incidence in redder galaxies with lower specific star formation rates, and we argue that this result might emerge because of our ability to properly account for AGN light contamination and dust extinction. Interestingly, we find that the probability for a galaxy to host a black hole growing at any given "specific accretion rate" (i.e. the ratio of X-ray luminosity to the host stellar mass) is almost independent of the host galaxy mass, while it decreases as a power-law with Lx/M. By analyzing the normalization of such probability distribution, we show how the incidence of AGN increases with redshift as rapidly as (1+z)^4, in close resemblance with the overall evolution of the specific star formation rate of the entire galaxy population. Although AGN activity and star formation in galaxies do appear to have a common triggering mechanism, at least in a statistical sense, within the COSMOS sample we do not find strong evidence of any 'smoking gun' signaling powerful AGN influence on the star-forming properties of their hosts galaxies.

Models of galaxy evolution assume some connection between the AGN and star formation activity in galaxies. We use the multi-wavelength information of the CDFS to assess this issue. We select the AGNs from the 3Ms XMM-Newton survey and measure the star-formation rates of their hosts using data that probe rest-frame wavelengths longward of 20 um. Star-formation rates are obtained from spectral energy distribution fits, identifying and subtracting an AGN component. We divide the star-formation rates by the stellar masses of the hosts to derive specific star-formation rates (sSFR) and find evidence for a positive correlation between the AGN activity (proxied by the X-ray luminosity) and the sSFR for the most active systems with X-ray luminosities exceeding Lx=10^43 erg/s and redshifts z~1. We do not find evidence for such a correlation for lower luminosity systems or those at lower redshifts. We do not find any correlation between the SFR (or the sSFR) and the X-ray absorption derived from high-quality XMM-Newton spectra either, showing that the absorption is likely to be linked to the nuclear region rather than the host, while the star-formation is not nuclear. Comparing the sSFR of the hosts to the characteristic sSFR of star-forming galaxies at the same redshift we find that the AGNs reside mostly in main-sequence and starburst hosts, reflecting the AGN - sSFR connection. Limiting our analysis to the highest X-ray luminosity AGNs (X-ray QSOs with Lx>10^44 erg/s), we find that the highest-redshift QSOs (with z>2) reside predominantly in starburst hosts, with an average sSFR more than double that of the "main sequence", and we find a few cases of QSOs at z~1.5 with specific star-formation rates compatible with the main-sequence, or even in the "quiescent" region. (abridged)

The 4MOST consortium is currently halfway through a Conceptual Design study for ESO with the aim to develop a wide-field (>3 square degree, goal >5 square degree), high-multiplex (>1500 fibres, goal 3000 fibres) spectroscopic survey facility for an ESO 4m-class telescope (VISTA). 4MOST will run permanently on the telescope to perform a 5 year public survey yielding more than 20 million spectra at resolution R~5000 ({\lambda}=390-1000 nm) and more than 2 million spectra at R~20,000 (395-456.5 nm & 587-673 nm). The 4MOST design is especially intended to complement three key all-sky, space-based observatories of prime European interest: Gaia, eROSITA and Euclid. Initial design and performance estimates for the wide-field corrector concepts are presented. We consider two fibre positioner concepts, a well-known Phi-Theta system and a new R-Theta concept with a large patrol area. The spectrographs are fixed configuration two-arm spectrographs, with dedicated spectrographs for the high- and low-resolution. A full facility simulator is being developed to guide trade-off decisions regarding the optimal field-of-view, number of fibres needed, and the relative fraction of high-to-low resolution fibres. Mock catalogues with template spectra from seven Design Reference Surveys are simulated to verify the science requirements of 4MOST. The 4MOST consortium aims to deliver the full 4MOST facility by the end of 2018 and start delivering high-level data products for both consortium and ESO community targets a year later with yearly increments.

Bolometric luminosities and Eddington ratios of both X-ray selected broad-line (Type-1) and narrow-line (Type-2) AGN from the XMM-Newton survey in the COSMOS field are presented. The sample is composed by 929 AGN (382 Type-1 AGN and 547 Type-2 AGN) and it covers a wide range of redshifts, X-ray luminosities and absorbing column densities. About 65% of the sources are spectroscopically identified as either Type-1 or Type-2 AGN (83% and 52% respectively), while accurate photometric redshifts are available for the rest of the sample. The study of such a large sample of X-ray selected AGN with a high quality multi-wavelength coverage from the far-infrared (now with the inclusion of Herschel data at 100 micron and 160 micron) to the optical-UV allows us to obtain accurate estimates of bolometric luminosities, bolometric corrections and Eddington ratios. The kbol-Lbol relations derived in the present work are calibrated for the first time against a sizable AGN sample, and rely on observed redshifts, X-ray luminosities and column density distributions. We find that kbol is significantly lower at high Lbol with respect to previous estimates by Marconi et al. (2004) and Hopkins et al. (2007). Black hole masses and Eddington ratios are available for 170 Type-1 AGN, while black hole masses for Type-2 AGN are computed for 481 objects using the black hole mass-stellar mass relation and the morphological information. We confirm a trend between kbol and lambda_Edd, with lower hard X-ray bolometric corrections at lower Eddington ratios for both Type-1 and Type-2 AGN. We find that, on average, Eddington ratio increases with redshift for all Types of AGN at any given Mbh, while no clear evolution with redshift is seen at any given Lbol.

[Abriged] Supermassive black holes (SMBH) lurk in the nuclei of most massive galaxies, perhaps in all of them. The tight observed scaling relations between SMBH masses and structural properties of their host spheroids likely indicate that the processes fostering the growth of both components are physically linked, despite the many orders of magnitude difference in their physical size. This chapter discusses how we constrain the evolution of SMBH, probed by their actively growing phases, when they shine as active galactic nuclei (AGN) with luminosities often in excess of that of the entire stellar population of their host galaxies. Following loosely the chronological developments of the field, we begin by discussing early evolutionary studies, when AGN represented beacons of light probing the most distant reaches of the universe and were used as tracers of the large scale structure. This early study turned into AGN "Demography", once it was realized that the strong evolution (in luminosity, number density) of the AGN population hindered any attempt to derive cosmological parameters from AGN observations directly. Following a discussion of the state of the art in the study of AGN luminosity functions, we move on to discuss the "modern" view of AGN evolution, one in which a bigger emphasis is given to the physical relationships between the population of growing black holes and their environment. This includes observational and theoretical efforts aimed at constraining and understanding the evolution of scaling relations, as well as the resulting limits on the evolution of the SMBH mass function. Physical models of AGN feedback and the ongoing efforts to isolate them observationally are discussed next. Finally, we touch upon the problem of when and how the first black holes formed and the role of black holes in the high-redshift universe.

We present the results of our studies of the aperiodic optical flux variability for SS Cyg, an accreting binary systemwith a white dwarf. The main set of observational data presented here was obtained with the ANDOR/iXon DU-888 photometer mounted on the RTT-150 telescope, which allowed a record(for CCD photometers) time resolution up to 8 ms to be achieved. The power spectra of the source's flux variability have revealed that the aperiodic variability contains information about the inner boundary of the optically thick flow in the binary system. We show that the inner boundary of the optically thick accretion disk comes close to the white dwarf surface at the maximum of the source's bolometric light curve, i.e., at the peak of the instantaneous accretion rate onto the white dwarf, while the optically thick accretion disk is truncated at distances 8.5e9 cm ~10 R_{WD} in the low state. We suggest that the location of the inner boundary of the accretion disk in the binary can be traced by studying the parameters of the power spectra for accreting white dwarfs. In particular, this allows the mass of the accreting object to be estimated.

The black hole mass function of supermassive black holes describes the evolution of the distribution of black hole mass. It is one of the primary empirical tools available for mapping the growth of supermassive black holes and for constraining theoretical models of their evolution. In this review we discuss methods for estimating the black hole mass function, including their advantages and disadvantages. We also review the results of using these methods for estimating the mass function of both active and inactive black holes. In addition, we review current theoretical models for the growth of supermassive black holes that predict the black hole mass function. We conclude with a discussion of directions for future research which will lead to improvement in both empirical and theoretical determinations of the mass function of supermassive black holes.

We present a rest-frame spectral stacking analysis of ~1000 X-ray sources detected in the XMM-COSMOS field in order to investigate the iron K line properties of active galaxies beyond redshift z~1. In Type I AGN that have a typical X-ray luminosity of Lx~1.5e44 erg/s and z~1.6, the cold Fe K at 6.4 keV is weak (EW~0.05keV), in agreement with the known trend. In contrast, high-ionization lines of Fe XXV and Fe XXVI are pronounced. These high-ionization Fe K lines appear to have a connection with high accretion rates. While no broad Fe emission is detected in the total spectrum, it might be present, albeit at low significance, when the X-ray luminosity is restricted to the range below 3e44 erg/s, or when an intermediate range of Eddington ratio around 0.1 is selected. In Type II AGN, both cold and high-ionzation lines become weak with increasing X-ray luminosity. However, strong high-ionization Fe K (EW~0.3 keV) is detected in the spectrum of objects at z>2, while no 6.4 keV line is found. It is then found that the primary source of the high-ionization Fe K emission is those objects detected with Spitzer-MIPS at 24 micron. Given their median redshift of z=2.5, their bolometric luminosity is likely to reach 10^13 Lsun and the MIPS-detected emission most likely originates from hot dust heated by embedded AGN, probably accreting at high Eddington ratio. These properties match those of rapidly growing black holes in ultra-luminous infrared galaxies at the interesting epoch (z=2-3) of galaxy formation.

The aim of this work is to characterize the average Fe K emission properties of active galactic nuclei (AGNs) in the source rest-frame. We selected a sample of 248 AGNs from the 2XMM catalog, covering a wide redshift range 0 < z < 5 and with the EPIC-PN net 2-10 keV rest-frame counts >=200 and power law photon indices in the range 1.5-2.2. We employed two fully independent rest-frame stacking procedures to compute the mean Fe K profile. The counting statistics for the integrated spectrum is comparable to the one available for the best studied local Seyferts. To identify the artifacts possibly introduced by the stacking procedure, we have carried out simulations. We report that the average Fe K line profile in our sample is best represented by a combination of a narrow and a broad line. The equivalent widths of the narrow and broad (parametrized with a diskline) components are ~30 eV and ~100 eV, respectively. We also discuss the results of more complex fits and the implications of the adopted continuum modeling on the broad line parameters and its detection significance.

We present new constraints on the ratio of black hole (BH) mass to total galaxy stellar mass at 0.3<z<0.9 for a sample of 32 type-1 active galactic nuclei (AGNs) from the XMM-COSMOS survey covering the range M_BH/M_sun~10^(7.2--8.7). Virial M_BH estimates based on H_beta are available from the COSMOS Magellan/IMACS survey. We use high-resolution Hubble Space Telescope (HST) imaging to decompose the light of each type-1 AGN and host galaxy, and employ a specially-built mass-to-light ratio to estimate the stellar masses (M_*). The M_BH-M_* ratio shows a zero offset with respect to the local relation for galactic bulge masses, and we also find no evolution in the mass ratio M_BH/M_*=(1+z)^{0.02+-0.34} up to z~0.9. Interestingly, at the high-M_BH end there is a positive offset from the z=0 relation, which can be fully explained by a mass function bias with a cosmic scatter of 0.3, reaffirming that the intrinsic distribution is consistent with zero evolution. From our results we conclude that since z~0.9 no substantial addition of stellar mass is required: the decline in star formation rates and merger activity at z<1 support this scenario. Nevertheless, given that a significant fraction of these galaxies show a disk component, their bulges are indeed undermassive. This is a direct indication that for the last 7 Gyr the only essential mechanism required in order that these galaxies obey the z=0 relation is a redistribution of stellar mass to the bulge, likely driven by secular processes, i.e., internal instabilities and minor merging.

We explore the connection between black hole growth at the center of obscured quasars selected from the XMM-COSMOS survey and the physical properties of their host galaxies. We study a bolometric regime (<Lbol > 8 x 10^45 erg/s) where several theoretical models invoke major galaxy mergers as the main fueling channel for black hole accretion. We confirm that obscured quasars mainly reside in massive galaxies (Mstar>10^10 Msun) and that the fraction of galaxies hosting such powerful quasars monotonically increases with the stellar mass. We stress the limitation of the use of rest-frame color-magnitude diagrams as a diagnostic tool for studying galaxy evolution and inferring the influence that AGN activity can have on such a process. We instead use the correlation between star-formation rate and stellar mass found for star-forming galaxies to discuss the physical properties of the hosts. We find that at z ~1, ~62% of Type-2 QSOs hosts are actively forming stars and that their rates are comparable to those measured for normal star-forming galaxies. The fraction of star-forming hosts increases with redshift: ~71% at z ~2, and 100% at z ~3. We also find that the the evolution from z ~1 to z ~3 of the specific star-formation rate of the Type-2 QSO hosts is in excellent agreement with that measured for star-forming galaxies. From the morphological analysis, we conclude that most of the objects are bulge-dominated galaxies, and that only a few of them exhibit signs of recent mergers or disks. Finally, bulge-dominated galaxies tend to host Type-2 QSOs with low Eddington ratios (lambda<0.1), while disk-dominated or merging galaxies have at their centers BHs accreting at high Eddington ratios (lambda > 0.1).

We present a study of the multi-wavelength properties, from the mid-infrared to the hard X-rays, of a sample of 255 spectroscopically identified X-ray selected Type-2 AGN from the XMM-COSMOS survey. Most of them are obscured the X-ray absorbing column density is determined by either X-ray spectral analyses (for the 45% of the sample), or from hardness ratios. Spectral Energy Distributions (SEDs) are computed for all sources in the sample. The average SEDs in the optical band is dominated by the host-galaxy light, especially at low X-ray luminosities and redshifts. There is also a trend between X-ray and mid-infrared luminosity: the AGN contribution in the infrared is higher at higher X-ray luminosities. We calculate bolometric luminosities, bolometric corrections, stellar masses and star formation rates (SFRs) for these sources using a multi-component modeling to properly disentangle the emission associated to stellar light from that due to black hole accretion. For 90% of the sample we also have the morphological classifications obtained with an upgraded version of the Zurich Estimator of Structural Types (ZEST+). We find that on average Type-2 AGN have lower bolometric corrections than Type-1 AGN. Moreover, we confirm that the morphologies of AGN host-galaxies indicate that there is a preference for these Type-2 AGN to be hosted in bulge-dominated galaxies with stellar masses greater than 10^10 solar masses.

We present a study of the redshift evolution of the projected correlation function of 593 X-ray selected AGN with I_AB<23 and spectroscopic redshifts z<4, extracted from the 0.5-2 keV X-ray mosaic of the 2.13 deg^2 XMM-COSMOS survey. We introduce a method to estimate the average bias of the AGN sample and the mass of AGN hosting halos, solving the sample variance using the halo model and taking into account the growth of the structure over time. We find evidence of a redshift evolution of the bias factor for the total population of XMM-COSMOS AGN from b(z=0.92)=2.30 +/- 0.11 to b(z=1.94)=4.37 +/- 0.27 with an average mass of the hosting DM halos logM [h^-1 M_sun] ~ 13.12 +/- 0.12 that remains constant at all z < 2. Splitting our sample into broad optical lines AGN (BL), AGN without broad optical lines (NL) and X-ray unobscured and obscured AGN, we observe an increase of the bias with redshift in the range z=0.7-2.25 and z=0.6-1.5 which corresponds to a constant halo mass logM [h^-1 M_sun] ~ 13.28 +/- 0.07 and logM [h^-1 M_sun] ~ 13.00 +/- 0.06 for BL /X-ray unobscured AGN and NL/X-ray obscured AGN, respectively. The theoretical models which assume a quasar phase triggered by major mergers can not reproduce the high bias factors and DM halo masses found for X-ray selected BL AGN with L_BOL ~ 2e45 erg s^-1. Our work extends up to z ~ 2.2 the z <= 1 statement that, for moderate luminosity X-ray selected BL AGN, the contribution from major mergers is outnumbered by other processes, possibly secular such as tidal disruptions or disk instabilities.

We present the results of the analysis of a sample of 17 low-luminosity (L_x < 1e42 erg/s), radio loud AGNs in massive galaxies. The sample is extracted from the SDSS database and it spans uniformly a wide range in optical [OIII] emission line and radio luminosity, but within a narrow redshift range (0.05 < z < 0.11) and a narrow super massive black hole mass range (~ 1e8 M_sun). For these sources we measured core X-ray emission with the Chandra X-ray telescope and radio emission with the VLA. Our main goal is to establish which emission component, if any, can be regarded as the most reliable accretion/jet-power estimator at these regimes. In order to do so, we studied the correlation between emission line properties, radio luminosity, radio spectral slopes and X-ray luminosity, as well as more complex multi-variate relations involving black hole mass, such as the fundamental plane of black hole activity. We find that 15 out of 17 sources of our sample can be classified as Low-Excitation Galaxies (LEG) and their observed properties suggest X-ray and radio emission to originate from the jet basis. We also find that X-ray emission does not appear to be affected by nuclear obscuration and can be used as a reliable jet-power estimator. More generally, X-ray, radio and optical emission appear to be related, although no tight correlation is found. In accordance with a number of recent studies of this class of objects these findings may be explained by a lack of cold (molecular) gaseous structures in the innermost region of these massive galaxies.

We show how accretion rate governs the physical properties of a sample of unobscured broad-line, narrow-line, and lineless active galactic nuclei (AGNs). We avoid the systematic errors plaguing previous studies of AGN accretion rate by using accurate accretion luminosities (L_int) from well-sampled multiwavelength SEDs from the Cosmic Evolution Survey (COSMOS), and accurate black hole masses derived from virial scaling relations (for broad-line AGNs) or host-AGN relations (for narrow-line and lineless AGNs). In general, broad emission lines are present only at the highest accretion rates (L_int/L_Edd > 0.01), and these rapidly accreting AGNs are observed as broad-line AGNs or possibly as obscured narrow-line AGNs. Narrow-line and lineless AGNs at lower specific accretion rates (L_int/L_Edd < 0.01) are unobscured and yet lack a broad line region. The disappearance of the broad emission lines is caused by an expanding radiatively inefficient accretion flow (RIAF) at the inner radius of the accretion disk. The presence of the RIAF also drives L_int/L_Edd < 10^-2 narrow-line and lineless AGNs to 10 times higher ratios of radio to optical/UV emission than L_int/L_Edd > 0.01 broad-line AGNs, since the unbound nature of the RIAF means it is easier to form a radio outflow. The IR torus signature also tends to become weaker or disappear from L_int/L_Edd < 0.01 AGNs, although there may be additional mid-IR synchrotron emission associated with the RIAF. Together these results suggest that specific accretion rate is an important physical "axis" of AGN unification, described by a simple model.

We present the first results of the spectroscopy of distant, obscured AGN as obtained with the ultra-deep (~3.3 Ms) XMM-Newton survey in the Chandra Deep Field South (CDFS). One of the primary goals of the project is to characterize the X-ray spectral properties of obscured and heavily obscured Compton-thick AGN over the range of redhifts and luminosities that are relevant in terms of their contribution to the X-ray background. The ultra-deep exposure, coupled with the XMM detector's spectral throughput, allowed us to accumulate good quality X-ray spectra for a large number of X-ray sources and, in particular, for heavily obscured AGN at cosmological redshifts. Specifically we present the X-ray spectral properties of two high-redshift - z= 1.53 and z=3.70 - sources. The XMM spectra of both are very hard, with a strong iron Kalpha line at a rest-frame energy of 6.4 keV. A reflection-dominated continuum provides the best description of the X-ray spectrum of the z=1.53 source, while the intrinsic continuum of the z=3.70 AGN is obscured by a large column N_H ~ 10^24 cm-2 of cold gas. Compton-thick absorption, or close to it, is unambiguously detected in both sources. Interestingly, these sources would not be selected as candidate Compton thick AGN by some multiwavelength selection criteria based on the mid-infrared to optical and X-ray to optical flux ratios.

It is well accepted that unabsorbed as well as absorbed AGN are needed to explain the nature and the shape of the Cosmic X-ray background, even if the fraction of highly absorbed objects (dubbed Compton-thick sources) substantially still escapes detection. We derive and analyze the absorption distribution using a complete sample of AGN detected by Swift-BAT in the first three years of the survey. The fraction of Compton-thick AGN represents only 4.6% of the total AGN population detected by Swift-BAT. However, we show that once corrected for the bias against the detection of very absorbed sources the real intrinsic fraction of Compton-thick AGN is 20$^{+9}_{-6}$%. We proved for the first time (also in the BAT band) that the anti-correlation of the fraction of absorbed AGN and luminosity it tightly connected to the different behavior of the luminosity functions (XLFs) of absorbed and unabsorbed AGN. This points towards a difference between the two subsamples of objects with absorbed AGN being, on average, intrinsically less luminous than unobscured ones. Moreover the XLFs show that the fraction of obscured AGN might also decrease at very low luminosity. This can be successfully interpreted in the framework of a disk cloud outflow scenario as the disappearance of the obscuring region below a critical luminosity. Our results are discussed in the framework of population synthesis models and the origin of the Cosmic X-ray Background.

We investigate the behaviour of the accretion discs in the outbursts of the low-mass black-hole X-ray binaries (BHXRB), an overview of which we have presented previously. Almost all of the systems in which there are sufficient observations in the most disc dominated states show a variation of the disc luminosity with temperature close to L ~\propto T^4. This in turn implies that in these states, the disc radius, R_in, and the colour correction factor, f_col, are almost constant. Deviations away from the T^4 law are observed at the beginning and end of the most disc dominated states, during the intermediate states. Although these could be explained by an inward motion of the accretion disc, they are more likely to be the result of an increase in the value of f_col as the disc fraction decreases. By comparing the expected and observed disc luminosities, we place approximate limits on the allowed distances and masses of the BHXRB system. In a number of cases, the measured distances and masses of the BHXRB system indicate that it is possible that the black hole may be spinning.

Our analysis is aimed at characterizing the properties of the integrated spectrum of active galactic nuclei (AGNs) such as the ubiquity of the Fe K{\alpha} emission in AGNs and the dependence of the spectral parameters on the X-ray luminosity and redshift. We selected 2646 point sources from the 2XMM catalogue at high galactic latitude (|BII| > 25 degrees) and with the sum of EPIC-PN and EPIC-MOS 0.2-12 keV counts greater than 1000. Redshifts were obtained for 916 sources from the NED. The final sample consists of 507 AGN. Individual source spectra have been summed in the observed frame to compute the integrated spectra in different redshift and luminosity bins over the range 0<z<5. Detailed analysis of these spectra has been performed. We find that the narrow Fe K{\alpha} line at 6.4 keV is significantly detected up to z=1. The line equivalent width decreases with increasing X-ray luminosity in the 2-10 keV band (''IT effect''). The anti-correlation is characterized by the relation log(EWFe) = (1.66 +/- 0.09) + (-0.43 +/- 0.07) log(LX,44), where EWFe is the rest frame equivalent width of the neutral iron K{\alpha} line in eV and LX,44 is the 2-10 keV X-ray luminosity in units of 10^{44} erg s^{-1}. The equivalent width is nearly independent of redshift up to z ~ 0.8 with an average value of 101+/-40 (rms dispersion) eV in the luminosity range 43.5<= logLX <= 44.5. Our analysis also confirmed the hardening of the spectral indices at low luminosities implying a dependence of obscuration on luminosity. We confirm that the neutral narrow Fe K{\alpha} line is an almost ubiquitous feature of AGNs. We find compelling evidence for the ''IT effect'' over a redshift interval larger than probed in any previous study. We detect no evolution of the average rest frame equivalent width of the Fe K{\alpha} line with redshift.

We report the final optical identifications of the medium-depth (~60 ksec), contiguous (2 deg^2) XMM-Newton survey of the COSMOS field. XMM-Newton has detected ~800 X-ray sources down to limiting fluxes of ~5x10^{-16}, ~3x10^{-15}, and ~7x10^{-15} erg/cm2/s in the 0.5-2 keV, 2-10 keV and 5-10 keV bands, respectively. The work is complemented by an extensive collection of multi-wavelength data from 24 micron to UV, available from the COSMOS survey, for each of the X-ray sources, including spectroscopic redshifts for ~50% of the sample, and high-quality photometric redshifts for the rest. The XMM and multiwavelength flux limits are well matched: 1760 (98%) of the X-ray sources have optical counterparts, 1711 (~95%) have IRAC counterparts, and 1394 (~78%) have MIPS 24micron detections. Thanks to the redshift completeness (almost 100%) we were able to constrain the high-luminosity tail of the X-ray luminosity function confirming that the peak of the number density of logL_X>44.5 AGN is at z~2. Spectroscopically-identified obscured and unobscured AGN, as well as normal and starforming galaxies, present well-defined optical and infrared properties. We devised a robust method to identify a sample of ~150 high redshift (z>1), obscured AGN candidates for which optical spectroscopy is not available. We were able to determine that the fraction of the obscured AGN population at the highest (L_X>10^{44} erg s^{-1}) X-ray luminosity is ~15-30% when selection effects are taken into account, providing an important observational constraint for X-ray background synthesis. We studied in detail the optical spectrum and the overall spectral energy distribution of a prototypical Type 2 QSO, caught in a stage transitioning from being starburst dominated to AGN dominated, which was possible to isolate only thanks to the combination of X-ray and infrared observations.

Using Chandra observations we have identified a sample of seven off-nuclear X-ray sources, in the redshift range z=0.072-0.283, located within optically bright galaxies in the COSMOS Survey. Using the multi-wavelength coverage available in the COSMOS field, we study the properties of the host galaxies of these ULXs. In detail, we derived their star formation rate from H_alpha measurements and their stellar masses using SED fitting techniques with the aim to compute the probability to have an off-nuclear source based on the host galaxy properties. We divide the host galaxies in different morphological classes using the available ACS/HST imaging. We find that our ULXs candidates are located in regions of the SFR versus M$_star$ plane where one or more off-nuclear detectable sources are expected. From a morphological analysis of the ACS imaging and the use of rest-frame colours, we find that our ULXs are hosted both in late and early type galaxies. Finally, we find that the fraction of galaxies hosting a ULX ranges from ~0.5% to ~0.2% going from L[0.5-2 keV]=3 x 10^39 erg s^-1 to L[0.5-2 keV]= 2 x 10^40 erg s^-1.

We derive the growth of SMBHs relative to the stellar content of their host galaxy predicted under the assumption of BH accretion triggered by galaxy encounters occurring during their merging histories. We show that, within this framework, the ratio Gamma=(M_BH/M_*)(z)/(M_BH/M_*)(z=0) between the Black Hole mass and the galactic stellar mass (normalized to the local value) depends on both BH mass and redshift. While the average value and the spread of Gamma(z) increase with z, such an effect is larger for massive BHs, reaching values Gamma=5 for massive Black Holes (M>10^9 M_{\odot}) at z>4, in agreement with recent observations of high-redshift QSOs; this is due to the effectiveness of interactions in triggering BH accretion in high-density environments at high redshifts. To test such a model against observations, we worked out specific predictions for sub-samples of the simulated galaxies corresponding to the different observational samples for which measurements of Gamma have been obtained. We found that for Broad Line AGNs at 1<z<2 values of Gamma=2 are expected, with a mild trend toward larger value for increasing BH mass. Instead, when we select from our Monte Carlo simulations only extremely gas rich, rapidly star forming galaxies at 2<z<3, we find low values 0.3<Gamma<1.5, consistent with recent observational findings on samples of sub-mm galaxies; in the framework of our model, these objects end up at z=0 in low-to-intermediate mass BHs (M<10^9 M_{\odot}), and they do not represent typical paths leading to local massive galaxies. The latter have formed preferentially through paths passing above the local M_*-M_BH relation. We discuss how the global picture emerging from the model is consistent with a downsizing scenario, where massive BHs accrete a larger fraction of their final mass at high redshifts z>4.

We present a study of the X-ray to optical properties of a sample of 545 X-ray selected Type 1 AGN, from the XMM-COSMOS survey, over a wide range of redshifts ($0.04<\z<4.25$) and X-ray luminosities ($40.6 \leq \Log \Lhard \leq 45.3$). About 60% of them are spectroscopically identified Type 1 AGN, while the others have a reliable photometric redshift and are classified as Type 1 AGN on the basis of their multi-band Spectral Energy Distributions. We discuss the relationship between UV and X-ray luminosity, as parameterized by the $\alphaox$ spectral slope, and its dependence on redshift and luminosity. We compare our findings with previous investigations of optically selected broad-line AGN (mostly from SDSS). A highly significant correlation between $\alphaox$ and $\lo$ is found, in agreement with previous investigations of optically selected samples. We calculate bolometric corrections, $\kbol$, for the whole sample using hard X-ray luminosities ($\Lhard$), and the Eddington ratios for a subsample of 150 objects for which black hole mass estimates are available. We confirm the trend of increasing bolometric correction with increasing Eddington ratio as proposed in previous works. A tight correlation is found between $\alphaox$ and $\kbol$, which can be used to estimate accurate bolometric corrections using only optical and X-ray data. We find a significant correlation between $\alphaox$ and Eddington ratio, in which $\alphaox$ increases for increasing Eddington ratios.

The co-evolution of host galaxies and the active black holes which reside in their centre is one of the most important topics in modern observational cosmology. Here we present a study of the properties of obscured Active Galactic Nuclei (AGN) detected in the CDFS 1Ms observation and their host galaxies. We limited the analysis to the MUSIC area, for which deep K-band observations obtained with ISAAC@VLT are available, ensuring accurate identifications of the counterparts of the X-ray sources as well as reliable determination of photometric redshifts and galaxy parameters, such as stellar masses and star formation rates. In particular, we: 1) refined the X-ray/infrared/optical association of 179 sources in the MUSIC area detected in the Chandra observation; 2) studied the host galaxies observed and rest frame colors and properties. We found that X-ray selected (L_X>10^{42} erg s^{-1}) AGN show Spitzer colors consistent with both AGN and starburst dominated infrared continuum; the latter would not have been selected as AGN from infrared diagnostics. The host galaxies of X-ray selected obscured AGN are all massive (M_*>10^{10} M_sun) and, in 50% of the cases, are also actively forming stars (1/SSFR<t_{Hubble}) in dusty environments. The median L/LEdd value of the active nucleus is between 2% and 10% depending on the assumed M_{BH}/M_{*} ratio. Finally, we found that the X-ray selected AGN fraction increases with the stellar mass up to a value of ~30% at z>1 and M_*>3x10^{11} M_sun, a fraction significantly higher than in the local Universe for AGN of similar luminosities.

I will describe the constraints available from a study of AGN evolution synthesis models on the growth of the supermassive black holes (SMBH) population in the two main 'modes' observed (kinetic- and radiatively-dominated, respectively). I will show how SMBH mass function evolves anti-hierarchically, i.e. the most massive holes grew earlier and faster than less massive ones, and I will also derive tight constraints on the average radiative efficiency of AGN. An outlook on the redshift evolution of the AGN kinetic luminosity function will also be discussed, thus providing a robust physical framework for phenomenological models of AGN feedback within structure formation. Finally, I will present new constraints on the evolution of the black hole-galaxy scaling relation at 1<z<2 derived by exploiting the full multi-wavelength coverage of the COSMOS survey on a complete sample of 90 type 1 AGN.

We present a catalog of 213 type-2 AGN selected from the zCOSMOS survey. The selected sample covers a wide redshift range (0.15<z<0.92) and is deeper than any other previous study, encompassing the luminosity range 10^{5.5} < Lsun< L[OIII] < 10^{9.1} Lsun. We explore the intrinsic properties of these AGN and the relation to their X-ray emission (derived from the XMM-COSMOS observations). We study their evolution by computing the [OIII]5007A line luminosity function (LF) and we constrain the fraction of obscured AGN as a function of luminosity and redshift. The sample was selected on the basis of the optical emission line ratios, after applying a cut to the signal-to-noise ratio (S/N) of the relevant lines. We used the standard diagnostic diagrams [OIII]/Hbeta versus [NII]/Halpha and ([OIII]/Hbeta versus [SII]/Halpha) to isolate AGN in the redshift range 0.15<z<0.45 and the diagnostic diagram [OIII]/Hbeta versus [OII]/Hbeta to extend the selection to higher redshift (0.5<z<0.92). Combining our sample with one drawn from SDSS, we found that the best description of the evolution of type-2 AGN is a luminosity-dependent density evolution model. Moreover, using the type-1 AGN LF we were able to constrain the fraction of type-2 AGN to the total (type-1 + type-2) AGN population. We found that the type-2 fraction decreases with luminosity, in agreement with the most recent results, and shows signs of a slight increase with redshift. However, the trend with luminosity is visible only after combining the SDSS+zCOSMOS samples. From the COSMOS data points alone, the type-2 fraction seems to be quite constant with luminosity.

We constrain the ratio of black hole (BH) mass to total stellar mass of type-1 AGN in the COSMOS survey at 1<z<2. For 10 AGN at mean redshift z~1.4 with both HST/ACS and HST/NICMOS imaging data we are able to compute total stellar mass M_(*,total), based on restframe UV-to-optical host galaxy colors which constrain mass-to-light ratios. All objects have virial BH mass-estimates available from the COSMOS Magellan/IMACS and zCOSMOS surveys. We find zero difference between the M_BH--M_(*,total)-relation at z~1.4 and the M_BH--M_(*,bulge)-relation in the local Universe. Our interpretation is: (a) If our objects were purely bulge-dominated, the M_BH--M_(*,bulge)-relation has not evolved since z~1.4. However, (b) since we have evidence for substantial disk components, the bulges of massive galaxies (logM_(*,total)=11.1+-0.25 or logM_BH~8.3+-0.2) must have grown over the last 9 Gyrs predominantly by redistribution of disk- into bulge-mass. Since all necessary stellar mass exists in the galaxy at z=1.4, no star-formation or addition of external stellar material is required, only a redistribution e.g. induced by minor and major merging or through disk instabilities. Merging, in addition to redistributing mass in the galaxy, will add both BH and stellar/bulge mass, but does not change the overall final M_BH/M_(*,bulge) ratio. Since the overall cosmic stellar and BH mass buildup trace each other tightly over time, our scenario of bulge-formation in massive galaxies is independent of any strong BH-feedback and means that the mechanism coupling BH and bulge mass until the present is very indirect.

(Abriged) We report on the measurement of the rest frame K-band luminosity and total stellar mass of the hosts of 89 broad line Active Galactic Nuclei detected in the zCOSMOS survey in the redshift range 1<z<2.2. The unprecedented multiwavelength coverage of the survey field allows us to disentangle the emission of the host galaxy from that of the nuclear black hole in their Spectral Energy Distributions. We derive an estimate of black hole masses through the analysis of the broad Mg II emission lines observed in the medium-resolution spectra taken with VIMOS/VLT as part of the zCOSMOS project. We found that, as compared to the local value, the average black hole to host galaxy mass ratio appears to evolve positively with redshift, with a best fit evolution of the form (1+z)^{0.68 \pm0.12 +0.6 -0.3}, where the large asymmetric systematic errors stem from the uncertainties in the choice of IMF, in the calibration of the virial relation used to estimate BH masses and in the mean QSO SED adopted. A thorough analysis of observational biases induced by intrinsic scatter in the scaling relations reinforces the conclusion that an evolution of the MBH-M* relation must ensue for actively growing black holes at early times: either its overall normalization, or its intrinsic scatter (or both) appear to increase with redshift. This can be interpreted as signature of either a more rapid growth of supermassive black holes at high redshift, a change of structural properties of AGN hosts at earlier times, or a significant mismatch between the typical growth times of nuclear black holes and host galaxies.

The physical and evolutionary relation between growing supermassive black holes (AGN) and host galaxies is currently the subject of intense research activity. Nevertheless, a deep theoretical understanding of such a relation is hampered by the unique multi-scale nature of the combined AGN-galaxy system, which defies any purely numerical, or semi-analytic approach. Various physical process active on different scales have signatures in different parts of the electromagnetic spectrum; thus, observations at different wavelengths and theoretical ideas all should contribute towards a "large dynamic range" view of the AGN phenomenon. As an example, I will focus in this review on two major recent observational results on the cosmic evolution of supermassive black holes, focusing on the novel contribution given to the field by the COSMOS survey. First of all, I will discuss the evidence for the so-called "downsizing" in the AGN population as derived from large X-ray surveys. I will then present new constraints on the evolution of the black hole-galaxy scaling relation at 1<z<2 derived by exploiting the full multi-wavelength coverage of the survey on a complete sample of ~90 type 1 AGN.

We present black hole masses and accretion rates for 182 Type 1 AGN in COSMOS. We estimate masses using the scaling relations for the broad Hb, MgII, and CIV emission lines in the redshift ranges 0.16<z<0.88, 1<z<2.4, and 2.7<z<4.9. We estimate the accretion rate using an Eddington ratio L_I/L_Edd estimated from optical and X-ray data. We find that very few Type 1 AGN accrete below L_I/L_Edd ~ 0.01, despite simulations of synthetic spectra which show that the survey is sensitive to such Type 1 AGN. At lower accretion rates the BLR may become obscured, diluted or nonexistent. We find evidence that Type 1 AGN at higher accretion rates have higher optical luminosities, as more of their emission comes from the cool (optical) accretion disk with respect to shorter wavelengths. We measure a larger range in accretion rate than previous works, suggesting that COSMOS is more efficient at finding low accretion rate Type 1 AGN. However the measured range in accretion rate is still comparable to the intrinsic scatter from the scaling relations, suggesting that Type 1 AGN accrete at a narrow range of Eddington ratio, with L_I/L_Edd ~ 0.1.

We present multi-wavelength observations (from optical to sub-millimeter, including Spitzer and SCUBA) of H2XMMJ 003357.2-120038 (also GD158_19), an X-ray selected, luminous narrow-line (Type 2) quasar at z=1.957 selected from the HELLAS2XMM survey. Its broad-band properties can be reasonably well modeled assuming three components: a stellar component to account for the optical and near-IR emission, an AGN component (i.e., dust heated by an accreting active nucleus), dominant in the mid-IR, with an optical depth at 9.7 micron along the line of sight (close to the equatorial plane of the obscuring matter) of tau(9.7)=1 and a full covering angle of the reprocessing matter (torus) of 140 degrees, and a far-IR starburst component (i.e., dust heated by star formation) to reproduce the wide bump observed longward of 70 micron. The derived star-formation rate is about 1500 solar masses per year. The overall modeling indicates that GD158_19 is a high-redshift X-ray luminous, obscured quasar with coeval powerful AGN activity and intense star formation. It is probably caught before the process of expelling the obscuring gas has started, thus quenching the star formation.

We explore the cosmic evolution of radio AGN with low radio powers (L_1.4GHz < 5\times10^25 W/Hz) out to z=1.3 using to-date the largest sample of ~600 low luminosity radio AGN at intermediate redshift drawn from the VLA-COSMOS survey. We derive the radio luminosity function for these AGN, and its evolution with cosmic time assuming two extreme cases: i) pure luminosity and ii) pure density evolution. The former and latter yield L_*\propto(1+z)^(0.8+/-0.1), and Phi*\propto (1+z)^(1.1+/-0.1), respectively, both implying a fairly modest change in properties of low radio-power AGN since z=1.3. We show that this is in stark contrast with the evolution of powerful (L_1.4GHz > 5\times10^25 W/Hz) radio AGN over the same cosmic time interval, constrained using the 3CRR, 6CE, and 7CRS radio surveys by Willott et al. (2001). We demonstrate that this can be explained through differences in black hole fueling and triggering mechanisms, and a dichotomy in host galaxy properties of weak and powerful AGN. Our findings suggest that high and low radio-power AGN activity is triggered in different stages during the formation of massive red galaxies. We show that weak radio AGN occur in the most massive galaxies already at z~1, and they may significantly contribute to the heating of their surrounding medium and thus inhibit gas accretion onto their host galaxies, as recently suggested for the `radio mode' in cosmological models.

We present a study of the host galaxies of AGN selected from the zCOSMOS survey to establish if accretion onto supermassive black holes and star formation are explicitly linked up to z~1. We identify 152 galaxies that harbor AGN, based on XMM observations of 7543 galaxies (i<22.5). Star formation rates (SFRs), including those weighted by stellar mass, are determined using the [OII]3727 emission-line, corrected for an AGN contribution. We find that the majority of AGN hosts have significant levels of star formation with a distribution spanning ~1-100 Msun yr^-1. The close association between AGN activity and star formation is further substantiated by an increase in the AGN fraction with the youthfulness of their stars as indicated by the rest-frame color (U-V) and spectral index Dn(4000); we demonstrate that mass-selection alleviates an artifical peak falling in the transition region when using luminosity-limited samples. We also find that the SFRs of AGN hosts evolve with cosmic time in a manner that closely mirrors the overall galaxy population and explains the low SFRs in AGNs (z<0.3) from the SDSS. We conclude that the conditions most conducive for AGN activity are a massive host galaxy and a large reservoir of gas. Furthermore, a direct correlation between mass accretion rate onto SMBHs and SFR is shown to be weak although the average ratio is constant with redshift, effectively shifting the evidence for a co-evolution scenario in a statistical manner to smaller physical scales. Our findings illustrate an intermittent scenario with an AGN lifetime substantially shorter than that of star formation and underlying complexities regarding fueling over vastly different physical scales yet to be determined [Abridged].

(abridged) We take advantage of the deep Chandra and Spitzer coverage of a large area (more than 10 times the area covered by the Chandra deep fields, CDFs in the COSMOS field, to extend the search of highly obscured, Compton-thick active nuclei to higher luminosity. These sources have low surface density and large samples can be provided only through large area surveys, like the COSMOS survey. We analyze the X-ray properties of COSMOS MIPS sources with 24$\mu$m fluxes higher than 550$\mu$Jy. For the MIPS sources not directly detected in the Chandra images we produce stacked images in soft and hard X-rays bands. To estimate the fraction of Compton-thick AGN in the MIPS source population we compare the observed stacked count rates and hardness ratios to those predicted by detailed Monte Carlo simulations including both obscured AGN and star-forming galaxies. The density of lower luminosity Compton-thick AGN (logL(2-10keV)=43.5-44) at z=0.7--1.2 is $(3.7\pm1.1) \times10^{-5}$ Mpc$^{-3}$, corresponding to $\sim67%$ of that of X-ray selected AGN. The comparison between the fraction of infrared selected, Compton thick AGN to the X-ray selected, unobscured and moderately obscured AGN at high and low luminosity suggests that Compton-thick AGN follow a luminosity dependence similar to that discovered for Compton-thin AGN, becoming relatively rarer at high luminosities. We estimate that the fraction of AGN (unobscured, moderately obscured and Compton thick) to the total MIPS source population is $49\pm10%$, a value significantly higher than that previously estimated at similar 24$\mu$m fluxes. We discuss how our findings can constrain AGN feedback models.

We present a new measurement of the space density of high redshift (3.0<z<4.5), X-ray selected QSOs obtained by exploiting the deep and uniform multiwavelength coverage of the COSMOS survey. We have assembled a statistically large (40 objects), X-ray selected (F_{0.5-2 keV} >10^{-15} cgs), homogeneous sample of z>3 QSOs for which spectroscopic (22) or photometric (18) redshifts are available. We present the optical (color-color diagrams) and X-ray properties, the number counts and space densities of the z>3 X-ray selected quasars population and compare our findings with previous works and model predictions. We find that the optical properties of X-ray selected quasars are not significantly different from those of optically selected samples. There is evidence for substantial X-ray absorption (logN_H>23 cm^{-2}) in about 20% of the sources in the sample. The comoving space density of luminous (L_X >10^{44} erg s^-1) QSOs declines exponentially (by an e--folding per unit redshift) in the z=3.0-4.5 range, with a behavior similar to that observed for optically bright unobscured QSOs selected in large area optical surveys. Prospects for future, large and deep X-ray surveys are also discussed.

The growth of supermassive black holes (SMBH) through accretion is accompanied by the release of enormous amounts of energy which can either be radiated away, as happens in quasars, advected into the black hole, or disposed of in kinetic form through powerful jets, as is observed, for example, in radio galaxies. Here, I will present new constraints on the evolution of the SMBH mass function and Eddington ratio distribution, obtained from a study of AGN luminosity functions aimed at accounting for both radiative and kinetic energy output of AGN in a systematic way. First, I discuss how a refined Soltan argument leads to joint constraints on the mass-weighted average spin of SMBH and of the total mass density of high redshift (z~5) and "wandering" black holes. Then, I will show how to describe the downsizing trend observed in the AGN population in terms of cosmological evolution of physical quantities (black hole mass, accretion rate, radiative and kinetic energy output). Finally, the redshift evolution of the AGN kinetic feedback will be briefly discussed and compared with the radiative output of the evolving SMBH population, thus providing a robust physical framework for phenomenological models of AGN feedback within structure formation.

We present a comprehensive synthesis model for the AGN evolution and the growth of supermassive black holes in the Universe. We solve the continuity equation for SMBH mass function using the locally determined one as a boundary condition, and the hard X-ray luminosity function as tracer of the AGN growth rate distribution, supplemented with a luminosity-dependent bolometric correction and an absorbing column distribution. Differently from most previous semi-analytic and numerical models, we do not assume any specific distribution of Eddington ratios, rather we determine it empirically by coupling the mass and luminosity functions. SMBH show a very broad accretion rate distribution, and we discuss the consequences of this fact for our understanding of observed AGN fractions in galaxies. We confirm previous results and demonstrate that, at least for z<1.5, SMBH mass function evolves anti-hierarchically, i.e. the most massive holes grew earlier and faster than less massive ones. For the first time, we find hints of a reversal of such a downsizing behaviour at redshifts above the peak of the black hole accretion rate density (z~2). We also derive tight constraints on the (mass weighted) average radiative efficiency of AGN: we find that 0.065<\xi_0 \epsilon_{rad}< 0.07$, where \xi_0 is the local SMBH mass density in units of 4.3x10^5 M_sun Mpc^{-3}. We trace the cosmological evolution of the kinetic luminosity function of AGN, and find that the overall efficiency of SMBH in converting accreted rest mass energy into kinetic power, ranges between 3 and 5 \times 10^{-3}. Such a ``kinetic efficiency'' varies however strongly with SMBH mass and redshift, being maximal for very massive holes at late times, as required for the AGN feedback by many galaxy formation models in Cosmological contexts. (Abriged)

(Abridged) We have studied the relationship among nuclear radio and X-ray power, Bondi rate and the kinetic luminosity of sub-Eddington active galactic nuclear (AGN) jets. Besides the recently discovered correlation between jet kinetic and Bondi power, we show that a clear correlation exists also between Eddington-scaled kinetic power and bolometric luminosity, given by: Log(L_kin/L_Edd)=0.49*Log(L_bol/L_Edd)-0.78. The measured slope suggests that these objects are in a radiatively inefficient accretion mode, and has been used to put stringent constraints on the properties of the accretion flow. We found no statistically significant correlations between Bondi power and bolometric AGN luminosity, apart from that induced by their common dependence on L_kin. Analyzing the relation between kinetic power and radio core luminosity, we are then able to determine, statistically, both the probability distribution of the mean jets Lorentz factor, peaking at \Gamma~7, and the intrinsic relation between kinetic and radio core luminosity, that we estimate as: Log(L_kin)=0.81*Log(L_R)+11.9, in good agreement with theoretical predictions of synchrotron jet models. With the aid of these findings, quantitative assessments of kinetic feedback from supermassive black holes in the radio mode will be possible based on accurate determinations of the central engine properties alone. As an example, Sgr A* may follow the correlations of radio mode AGN, based on its observed radiative output and on estimates of the accretion rate both at the Bondi radius and in the inner flow. If this is the case, the SMBH in the Galactic center is the source of ~ 5 times 10^38 ergs/s of mechanical power, equivalent to about 1.5 supernovae every 10^5 years.

We derive the kinetic luminosity function for flat spectrum radio jets, using the empirical and theoretical scaling relation between jet power and radio core luminosity. The normalization for this relation is derived from a sample of flat spectrum cores in galaxy clusters with jet-driven X-ray cavities. The total integrated jet power at z=0 is W_{tot} ~ 3x10^40 ergs/s/Mpc^{3}. By integrating W_{tot} over red-shift, we determine the total energy density deposited by jets as e_{tot} ~ 2x10^{58} ergs/Mpc^{3}. Both W_{tot} and e_{tot} are dominated by low luminosity sources. Comparing e_{tot} to the local black hole mass density rho_{BH} gives an average jet production efficiency of epsilon_{jet} = e_{jet}/rho_{BH}c^2 ~ 3%. Since black hole mass is accreted mainly during high luminosity states, epsilon_{jet} is likely much higher during low luminosity states.

We present a first attempt to derive the cosmological evolution of the kinetic luminosity function of AGN based on the joint evolution of the flat spectrum radio and hard X-ray selected AGN luminosity functions. An empirical correlation between jet power and radio core luminosity is found, which is consistent with the theoretical assumption that, below a certain Eddington ratio, SMBH accrete in a radiatively inefficient way, while most of the energy output is in the form of kinetic energy. We show how the redshift evolution of the kinetic power density from such a low-mdot mode of accretion makes it a good candidate to explain the so-called ``radio mode'' of AGN feedback as postulated in many galaxy formation schemes.

The X-ray spectra of luminous Seyfert 1 galaxies often appear to be reflection dominated. In a number of Narrow Line Seyfert 1 (NLS1) galaxies and galactic black holes in the very high state, the variability of the continuum and of the iron line are decoupled, the reflected component being often much less variable than the continuum. These properties have been interpreted as effects of gravitational light bending. In this framework, we present detailed Monte-Carlo simulations of the reflection continuum in the Kerr metric. These calculations confirm that the spectra and variability behaviour of these sources can be reproduced by the light bending model. As an alternative to the light bending model, we show that similar observational properties are expected from radiation pressure dominated discs subject to violent clumping instabilities and, as a result, have a highly inhomogeneous two-phase structure. In this model, most of the observed spectral and variability features originate from the complex geometrical structure of the inner regions of near-Eddington accretion flows and are therefore a signature of accretion physics rather than general relativity.

Observational evidence accumulated over the past decade indicates that accretion discs in X-ray binaries are viscously stable unless they accrete very close to the Eddington limit. This is at odds with the most basic standard accretion disc theory, but could be explained by either having the discs to be much cooler whereby they are not radiation pressure dominated, or by a more sophisticated viscosity law. Here we argue that the latter is taking place in practice, on the basis of a stability analysis that assumes that the magneto-rotational-instability (MRI) responsible for generating the turbulent stresses inside the discs is also the source for a magnetically dominated corona. We show that observations of stable discs in the high/soft states of black hole binaries, on the one hand, and of the strongly variable microquasar GRS 1915+105 on the other, can all be explained if the magnetic turbulent stresses inside the disc scale proportionally to the geometric mean of gas and total pressure with a constant of proportionality (viscosity parameter) having a value of a few times 10^{-2}. Implications for bright AGN are also briefly discussed

We discuss the expected X-ray spectral and variability properties of black hole accretion discs at high luminosity, under the hypothesis that radiation pressure dominated discs are subject to violent clumping instabilities and, as a result, have a highly inhomogeneous two-phase structure. After deriving the full accretion disc solutions explicitly in terms of the parameters of the model, we study their radiative properties both with a simple two-zones model, treatable analytically, and with radiative transfer simulations which account simultaneously for energy balance and Comptonisation in the hot phase, together with reflection, reprocessing, ionization and thermal balance in the cold phase. We show that, if not only the density, but also the heating rate within these flows is inhomogeneous, then complex reflection-dominated spectra can be obtained for a high enough covering fraction of the cold phase. In general, large reflection components in the observed X-ray spectra should be associated with strong soft excesses, resulting from the combined emission of ionized atomic emission lines. The variability properties of such systems are such that, even when contributing to a large fraction of the hard X-ray spectrum, the reflection component is less variable than the power-law like emission originating from the hot Comptonising phase, in agreement with what is observed in many Narrow Line Seyfert 1 galaxies and bright Seyfert 1. Our model falls within the family of those trying to explain the complex X-ray spectra of bright AGN with ionized reflection, but presents an alternative, specific, physically motivated, geometrical setup for the complex multi-phase structure of the inner regions of near-Eddington accretion flows.

Deep observations with the Very Large Array of A0620-00, performed in 2005 August, resulted in the first detection of radio emission from a black hole binary at X-ray luminosities as low as 10^{-8.5} times the Eddington limit. The measured radio flux density, of 51 \pm 7 microJy at 8.5 GHz, is the lowest reported for an X-ray binary system so far, and is interpreted in terms of partially self-absorbed synchrotron emission from outflowing plasma. Making use of the estimated outer accretion rate of A0620-00 in quiescence, we demonstrate that the outflow kinetic power must be energetically comparable to the total accretion power associated with such rate, if it was to reach the black hole with the standard radiative efficiency of 10 per cent. This favours a model for quiescence in which a radiatively inefficient outflow accounts for a sizable fraction of the missing energy, and, in turn, substantially affects the overall dynamics of the accretion flow. Simultaneous observations in the X-ray band, with Chandra, confirm the validity of a non-linear radio/X-ray correlation for hard state black hole binaries down to low quiescent luminosities, thereby contradicting some theoretical expectations. Taking the mass term into account, the A0620-00 data lie on the extrapolation of the so called Fundamental Plane of black hole activity, which has thus been extended by more than 2 orders of magnitude in radio and X-ray luminosity. With the addition of the A0620-00 point, the plane relation provides an empirical proof for the scale-invariance of the jet-accretion coupling in accreting black holes over the entire parameter space observable with current instrumentation.

We compare all the available observational data on the redshift evolution of the total stellar mass and star formation rate density in the Universe with the mass and accretion rate density evolution of supermassive black holes, estimated from the hard X-ray selected luminosity function of quasars and active galactic nuclei. We find that on average black hole mass must have been higher at higher redshift for given spheroid stellar mass. Moreover, we find negative redshift evolution of the disk/irregulars to spheroid mass ratio. The total accretion efficiency is constrained to be between 0.06 and 0.12, depending on the exact value of the local SMBH mass density, and on the critical accretion rate below which radiatively inefficient accretion may take place.

The fundamental plane of black hole activity is a non-linear correlation among radio core luminosity, X-ray luminosity and mass of all accreting black holes, both of stellar mass and supermassive, found by Merloni, Heinz and Di Matteo (2003) and, independently, by Falcke, Koerding and Markoff (2004). Here we further examine a number of statistical issues related to this correlation. In particular, we discuss the issue of sample selection and quantify the bias introduced by the effect of distance in two of the correlated quantities. We demonstrate that the fundamental plane relation cannot be a distance artifact, and that its non-linearity must represent an intrinsic characteristic of accreting black holes. We also discuss possible future observational strategies to improve our understanding of this correlation.

We interpret the rapid correlated UV/optical/ X-ray variability of XTEJ1118+480 as a signature of the coupling between the X-ray corona and a jet emitting synchrotron radiation in the optical band.We propose a scenario in which the jet and the X-ray corona are fed by the same energy reservoir where large amounts of accretion power are stored before being channelled into either the jet or the high energy radiation. This time dependent model reproduces the main features of the rapid multi-wavelength variability of XTEJ1118+480. A strong requirement of the model is that the total jet power should be at least a few times larger than the observed X-ray luminosity. This would be consistent with the overall low radiative efficiency of the source. We present independent arguments showing that the jet probably dominates the energetic output of all accreting black holes in the low-hard state.

We examine the disc--jet connection in stellar mass and supermassive black holes by investigating the properties of their compact emission in the hard X-ray and radio bands. We compile a sample of ~100 active galactic nuclei with measured mass, 5 GHz core emission, and 2-10 keV luminosity, together with 8 galactic black holes with a total of about 50 simultaneous observations in the radio and X-ray bands. Using this sample, we study the correlations between the radio (L_R) and the X-ray (L_X) luminosity and the black hole mass (M). We find that the radio luminosity is correlated with both M and L_X, at a highly significant level. We show how this result can be used to extend the standard unification by orientation scheme to encompass unification by mass and accretion rate.

We review the statistical properties of relativistic Doppler boosting relevant for studies of relativistic jets from compact objects based on radio--X-ray(--mass) correlations, such as that found in black-hole X-ray binaries in the low/hard state, or the ``fundamental plane'' of Merloni, Heinz, & DiMatteo. We show that the presence of only moderate scatter in such relations does not necessarily imply low Lorentz factors of the jets producing the radio emission in the samples under consideration. Applying Doppler beaming statistics to a large sample of XRBs and AGN, we derive a limit on the width of the Lorentz factor distribution of black holes with relativistic jets: If the X-rays are unbeamed (e.g., if they originate in the accretion disk or in the slower, innermost part of the jet), the width of the \beta\Gamma distribution should be about one order of magnitude or less. If the scatter about the ``fundamental plane'' is entirely dominated by relativistic beaming, a lower limit on the mean Lorentz factor <\beta\Gamma> > 5 can be derived. On the other hand, if the X-rays are boosted by the same factor as the radio emission, we show that the observed scatter cannot be reasonably explained by Doppler boosting alone.

We examine possible phenomenological constraints for the joint evolution of supermassive black holes (SMBH) and their host spheroids. We compare all the available observational data on the redshift evolution of the total stellar mass and star formation rate density in the Universe with the mass and accretion rate density evolution of supermassive black holes, estimated from the hard X-ray selected luminosity function of quasars and active galactic nuclei (AGN) for a given radiative efficiency, \epsilon. We assume that the ratio of the stellar mass in spheroids to the black hole mass density evolves as (1+z)^{-\alpha}, while the ratio of the stellar mass in disks+irregulars to that in spheroids evolves as (1+z)^{-\beta}, and we derive constraints on \alpha, \beta and \epsilon. We find that \alpha>0 at the more than 4-sigma level, implying a larger black hole mass at higher redshift for a given spheroid stellar mass. The favored values for \beta are typically negative, suggesting that the fraction of stellar mass in spheroids decreases with increasing redshift. This is consistent with recent determinations that show that the mass density at high redshift is dominated by galaxies with irregular morphology. In agreement with earlier work, we constrain \epsilon to be between 0.04 and 0.11, depending on the exact value of the local SMBH mass density, but almost independently of \alpha and \beta.

I present a new method to unveil the history of cosmic accretion and the build-up of Supermassive Black Holes (SMBH) in the nuclei of galaxies, based on observations of the evolving radio and (hard) X-ray luminosity functions of AGN. The fundamental plane of black hole activity discovered by Merloni, Heinz & Di Matteo (2003) is used as a mass and accretion rate estimator. I adopt the local BH mass function as a boundary condition to integrate backwards in time the continuity equation for the SMBH evolution, neglecting the role of mergers. Under the most general assumption that accretion proceeds in a radiatively efficient way above a certain rate, and in a radiatively inefficient way below, the redshift evolution of the mass and accretion rate functions are calculated self-consistently. The only tunable parameters are the accretion efficiency and the critical ratio of the X-ray to Eddington luminosity at which the transition between accretion modes takes place. The evolution of the BH mass function between z=0 and z ~ 3 shows clear signs of an anti-hierarchical behaviour: while the majority of the most massive objects (M > 10^9) were already in place at z ~ 3, lower mass ones mainly grew at progressively lower redshift. As an example, I will discuss the consequences of these results for the lifetimes of accreting black holes.

I present a new method to unveil the history of cosmic accretion and the build-up of SMBH in the nuclei of galaxies, based on observations of the evolving radio and (hard) X-ray luminosity functions of AGN. The fundamental plane of black hole activity discovered by Merloni, Heinz & Di Matteo (2003) is used as a mass and accretion rate estimator. I adopt the local BH mass function as a boundary condition to integrate backwards in time the continuity equation for the SMBH evolution, neglecting the role of mergers. Under the most general assumption that accretion proceeds in a radiatively efficient way above a certain rate, and in a radiatively inefficient way below, the redshift evolution of the mass and accretion rate functions are calculated self-consistently. The only tunable parameters are the accretion efficiency and the critical ratio of the X-ray to Eddington luminosity at which the transition between accretion modes takes place. For fiducial values of these parameters, I found that half (85%) of the local BH mass density was accumulated at redshift z<1 (z<3), mostly in radiatively efficient episodes of accretion. The evolution of the BH mass function between z=0 and z~3 shows clear signs of an anti-hierarchical behaviour: while the majority of the most massive objects (M > 10^9) were already in place at z~3, lower mass ones mainly grew at progressively lower redshift. Also, the average accretion rate decreases with time. Consequently, sources in the radiatively inefficient regime of accretion only begin to dominate the comoving accretion energy density in the universe at z<1. I discuss the implications of these results for the efficiency of accretion onto SMBH, the quasars lifetimes and duty cycles and the history of AGN feedback in the form of mechanical energy output (abriged).

We interpret the rapid correlated UV/optical/ X-ray variability of XTE J1118+480 as a signature of the coupling between the X-ray corona and a jet emitting synchrotron radiation in the optical band. We propose a scenario in which the jet and the X-ray corona are fed by the same energy reservoir where large amounts of accretion power are stored before being channelled into either the jet or the high energy radiation. This time dependent model reproduces the main features of the rapid multi-wavelength variability of XTE J1118+480. Assuming that the energy is stored in the form of magnetic field, we find that the required values of the model parameters are compatible with both a patchy corona atop a cold accretion disc and a hot thick inner disc geometry. The range of variability timescales for the X-ray emitting plasma are consistent with the dynamical times of an accretion flow between 10 and 100 Schwarzschild radii. On the other hand, the derived range of timescales associated with the dissipation in the jet extends to timescales more than 10 times larger, confirming the suggestion that the generation of a powerful outflow requires large scale coherent poloidal field structures. A strong requirement of the model is that the total jet power should be at least a few times larger than the observed X-ray luminosity. This would be consistent with the overall low radiative efficiency of the source. We present independent arguments showing that the jet probably dominates the energetic output of all accreting black holes in the low-hard state.

Studying a sample of both strongly and weakly active galactic nuclei with measured masses and 5 GHz and 2-10 keV core luminosities, together with a few galactic black holes simultaneously observed in the radio and X-ray bands, Merloni, Heinz, & Di Matteo (2003) showed that the sources are correlated through a ``fundamental plane'' relationship in the three-dimensional (log L_R, log L_X, log M) space. Here I elaborate on how such a relationship can be used to infer directly mass and accretion rate of any black hole given its radio and X-ray fluxes, complementing the information obtained from optical/UV surveys. As an example, I show how the local X-ray and radio luminosity functions, coupled with the black hole mass function derived from the M-\sigma relation, provide us with an accretion rate function. It is found that the typical X-ray Eddington ratio of an active black hole at redshift zero is about 5 \times 10^{-4}.

We examine the disc--jet connection in stellar mass and supermassive black holes by investigating the properties of their compact emission in the X-ray and radio bands. We compile a sample of ~100 active galactic nuclei with measured mass, 5 GHz core emission, and 2-10 keV luminosity, together with 8 galactic black holes with a total of ~50 simultaneous observations in the radio and X-ray bands. Using this sample, we study the correlations between the radio (L_{R}) and the X-ray (L_{X}) luminosity and the black hole mass (M). We find that the radio luminosity is correlated with {\em both} M and L_{X}, at a highly significant level. In particular, we show that the sources define a ``fundamental plane'' in the three-dimensional (log L_{R},log L_{X},log M) space, given by log L_{R}=(0.60^{+0.11}_{-0.11}) log L_{X} +(0.78^{+0.11}_{-0.09}) log M + 7.33^{+4.05}_{-4.07}, with a substantial scatter of \sigma_{R}=0.88. We compare our results to the theoretical relations between radio flux, black hole mass, and accretion rate derived by Heinz and Sunyaev (2003). Such relations depend only on the assumed accretion model and on the observed radio spectral index. Therefore, we are able to show that the X-ray emission from black holes accreting at less than a few per cent of the Eddington rate is unlikely to be produced by radiatively efficient accretion, and is marginally consistent with optically thin synchrotron emission from the jet. On the other hand, models for radiatively inefficient accretion flows seem to agree well with the data.

Recent X-ray observations have shown evidence for exceptionally broad and skewed iron K alpha emission lines from several accreting black hole systems. The lines are assumed to be due to fluorescence of the accretion disk illuminated by a surrounding corona and require a steep emissivity profile increasing in to the innermost radius. This appears to question both standard accretion disc theory and the zero torque assumption for the inner boundary condition, both of which predict a much less extreme profile. Instead it argues that a torque may be present due to magnetic coupling with matter in the plunging region or even to the spinning black hole itself. Discussion so far has centered on the torque acting on the disc. However the crucial determinant of the iron line profile is the radial variation of the power radiated in the corona. Here we study the effects of different inner boundary conditions on the coronal emissivity and on the profiles of the observable Fe K alpha lines. We argue that in the extreme case where a prominent highly redshifted component of the iron line is detected, requiring a steep emissivity profile in the innermost part and a flatter one outside, energy from the gas plunging into the black hole is being fed directly to the corona.

We present a systematic, analytical study of geometrically thin, optically thick accretion disc solutions for magnetized turbulent flows, with an alpha-like viscosity prescription. Under the only assumptions that (1) Magneto-Rotational instability (MRI) generates the turbulence that produces the anomalous viscosity needed for accretion to proceed, and that (2) the magnetic field amplified by the instability saturates due to buoyant vertical escape, we are able to self-consistently solve the disc structure equations including the fraction of power, f, that is carried off by vertical Poynting flux (and likely dissipated outside the optically thick disc). For low-viscosity discs, we obtain stable high-f solutions at low accretion rates, when gas pressure dominates, and unstable, low-f, radiation pressure dominated solutions at high accretion rates. For high viscosity discs, instead, a new thermally and viscously stable, radiation pressure dominated solution is found, characterized by f~1 and appearing only above a critical accretion rate (of the order of few tenths of the Eddington one). We discuss the regimes of validity of our assumptions, and the astrophysical relevance of our solutions. We conclude that our newly discovered thin disc solutions, possibly accompanied by powerful, magnetically dominated coronae and outflows, should be seriously considered as models for black holes accreting at super-Eddington rates.

I review our current theoretical understanding of the different possible rotating modes of accretion onto a black hole. I discuss both thick adiabatic flows and radiatively efficient thin disc solutions. I present a new self-consistent unified analytical description of two-phase thin disc--corona systems in which magneto-rotational instability is responsible for angular momentum transport in the disc and for the corona generation. Finally, I briefly discuss the role of magnetic fields in bridging the gap between accretion discs theory and jet production mechanisms

We apply a model of magnetically dominated coronae above standard accretion discs to the low/hard state of galactic black holes. When the disc-corona coupling is accounted for self-consistently assuming that magneto-rotational instability is at work in the disc, and that the corona is generated by buoyant escape of disc magnetic structures, then the model predicts powerful, X-ray emitting coronae at low accretion rates. A main consequence is discussed: the possibility that the corona itself is the launching site of powerful, MHD driven jets/outflows. This depends crucially of the coronal scaleheight. Finally, we present the first radial profiles of a corona a different accretion rates, and discuss their implications for high frequency variability.

A model for the inner regions of accretion flows is presented where, due to disc instabilities, cold and dense material is clumped into deep sheets or rings. Surrounding these density enhancements is hot, tenuous gas where coronal dissipation processes occur. We expect this situation to be most relevant when the accretion rate is close to Eddington and the disc is radiation-pressure dominated, and so may apply to Narrow-Line Seyfert~1 (NLS1) galaxies. In this scenario, the hard X-ray source is obscured for most observers, and so the detected X-ray emission would be dominated by reflection off the walls of the sheets. A simple Comptonization calculation shows that the large photon-indices characteristic of NLS1s would be a natural outcome of two reprocessors closely surrounding the hard X-ray source. We test this model by fitting the XMM-Newton spectrum of the NLS1 1H 0707-495 between 0.5 and 11 keV with reflection dominated ionized disc models. A very good fit is found with three different reflectors each subject to the same \Gamma=2.35 power-law. An iron overabundance is still required to fit the sharp drop in the spectrum at around 7 keV. We note that even a small corrugation of the accretion disc may result in \Gamma > 2 and a strong reflection component in the observed spectrum. Therefore, this model may also explain the strength and the variability characteristics of the MCG-6-30-15 Fe K line. The idea needs to be tested with further broadband XMM-Newton observations of NLS1s.

We present a model for magnetic structured coronae above accretion discs. On the shortest timescales, spatially and temporally correlated coronal flares can explain X-ray temporal and spectral variability observed in Seyfert galaxies. In particular, power density spectra, flux-spectral index and flux-variance correlations are naturally accounted for by the model. More dramatic spectral variations (i.e. state transitions in GBHC) are associated with parameters varying on longer timescales, such as accretion rate, coronal strength or geometry of the inner disc. In the framework of the standard Shakura--Sunyaev accretion disc theory, here we discuss why energetically dominant coronae at low accretion rates are ideal sites for launching powerful MHD driven outflows. Then, if the outflow is radiatively inefficient, then so is the source overall, even without advection being relevant for the dynamics of the accretion flow. This could be an alternative scenario for LLAGN and GBHC in their low/hard state, and may have consequences for our understanding of the accretion history of the universe.

The spectral energy distributions of galactic black holes in the low/hard state and of low-luminosity AGN possess many common features, the most prominent being: compact, flat (or inverted) spectrum radio cores with high brightness temperatures; excess red and infrared emission, often correlated with the radio flux; an extremely weak (or absent) quasi-thermal hump and a hard X-ray power-law with high energy cut-off. These sources are thought to be accreting at low rates and advection (or convection) dominated accretion flows are usually considered the best candidates to explain them. Here we present an alternative possibility, involving strong, unbound, magnetic coronae generated by geometrically thin, optically thick accretion discs at low accretion rates. First we show that, if angular momentum transport in the disc is due to magnetic turbulent stresses, the magnetic energy density and effective viscous stresses inside the disc are proportional to the geometric mean of the total (gas plus radiation) and gas pressure. Therefore the corona is less powerful in a radiation pressure dominated disc, and the relative fraction of the power liberated in the corona increases as the accretion rate decreases. Furthermore, we discuss reasons why energetically dominant coronae are ideal sites for launching powerful jets/outflows, both MHD and thermally driven. In analysing the spectral properties of such coronal outflow dominated accretion discs, we reach the important conclusion that if the jet/outflow is, as is likely, radiatively inefficient, then so is the source overall, even without advection of energy into the black hole being relevant for the dynamics of the accretion flow.

X-ray observations of Seyfert 1 galaxies offer the unique possibility of observing spectral variability on timescales comparable to the dynamical time of the inner accretion flow. They typically show highly variable lightcurves, with Power Density Spectra characterized by `red noise' and a break at low frequencies. Time resolved spectral analysis have established that spectral variability on the shortest timescales is important in all these sources, with the spectra getting softer at high fluxes. Here we present a model that is able to explain a number of the above mentioned properties in terms of magnetic flares shining above a standard accretion disc and producing the X-ray spectrum via inverse Compton scattering soft photons (both intrinsic and reprocessed thermal emission from the accretion disc and locally produced synchrotron radiation). We show that the fundamental heating event, likely caused by magnetic reconnection, must be compact, with typical size comparable to the accretion disc thickness and must be triggered at a height at least an order of magnitude larger than its size; the spatial and temporal distribution of flares are not random: the heating of the corona proceeds in correlated trains of events in an avalanche fashion. The amplitude of the avalanches obeys a power-law distribution and determines the size of the active regions where the spectrum is produced. With our model we simulate X-ray lightcurves that reproduce the main observational properties of the Power Density Spectra and of the X-ray continuum short-term variability of Seyfert 1 galaxies. By comparing them with observations of MGC--6-30-15, we are able to infer that the corona in this source must have a large optical depth (tau >1.5) and small average covering fraction.

We discuss the effects of a comptonizing corona on the appearance of the reflection components, and in particular of the reflection hump, in the X-rays spectra of accreting black holes. Indeed, in the framework of a thermal corona model, we expect that part (or even all, depending on the coronal covering factor) of the reflection features should cross the hot plasma, and thus suffer Compton scattering, before being observed. We have studied in detail the dependence of these effects on the physical (i.e. temperature and optical depth) and geometrical (i.e. inclination angle) parameters of the corona, concentrating on the slab geometry . Due to the smoothing and shifting towards high energies of the comptonized reflection hump, the main effects on the emerging spectra appear above 100 keV. We have also investigated the importance of such effects on the interpretation of the results obtained with the standard fitting procedures. We found that fitting Comptonization models, taking into account comptonized reflection, by the usual cut-off power law + uncomptonized reflection model, may lead to an underestimation of the reflection normalization and an overestimation of the high energy cut-off. We discuss and illustrate the importance of these effects by analysing recent observational results as those of the galaxy NGC 4258. We also find that the comptonizing corona can produce and/or emphasize correlations between the reflection features characteristics (like the iron line equivalent width or the covering fraction) and the X-ray spectral index similar to those recently reported in the literature. We also underline the importance of these effects when dealing with accurate spectral fitting of the X-ray background.

Although more than 2000 astronomical gamma-ray bursts (GRBs) have been detected, the precise progenitor responsible for these events is unknown. The temporal phenomenology observed in GRBs can significantly constrain the different models. Here we analyse the time histories of a sample of bright, long GRBs, searching for the ones exhibiting relatively long (more than 5 per cent of the total burst duration) quiescent times, defined as the intervals between adjacent episodes of emission during which the gamma-rays count rate drops to the background level. We find a quantitative relation between the duration of an emission episode and the quiescent time elapsed since the previous episode. We suggest here that the mechanism responsible for the extraction and the dissipation of energy has to take place in a meta-stable configuration, such that the longer the accumulation period, the higher is the stored energy available for the next emission episode.

Within the framework of the internal-external shocks model for gamma-ray bursts, we study the various mechanisms that can give rise to quiescent times in the observed gamma-ray light-curves. In particular, we look for the signatures that can provide us with evidence as to whether or not the central engine goes dormant for a period of time comparable to the duration of the gaps. We show that the properties of the prompt gamma-ray and X-ray emission can in principle determine whether the quiescent episodes are due to a modulated relativistic wind or a switching off of the central engine. We suggest that detailed observations of the prompt afterglow emission from the reverse shock will strongly constrain the possible mechanisms for the production of quiescent times in gamma-ray bursts.

We show how the simultaneous presence of a strong quasi periodic oscillation (QPO) of period of about 10 seconds in the optical and X-ray lightcurves of the X-ray transient XTE J1118+480 can be used to obtain information about the nature of the accretion flow around the source. The unusually high optical-to-X-ray flux ratio and the QPO observed simultaneously in both energy bands suggest that a significant fraction of the optical flux might originate close to the central source, where most of the X-rays are produced and be indicative of a magnetically dominated corona. We also show how the temporal evolution of the QPO can provide us with information on both the inner radius and the viscous properties of the optically thick accretion disc.

Most astrophysical sources powered by accretion onto a black hole, either of stellar mass or supermassive, when observed with hard X-rays show signs of a hot Comptonizing component in the flow, the so-called corona, with observed temperatures and optical depths lying in a narrow range (0.1 < \tau < 1 and 1x10^9 K < T < 3x10^9 K). Here we argue that these facts constitute strong supporting evidence for a magnetically-dominated corona. We show that the inferred thermal energy content of the corona, in all black hole systems, is far too low to explain their observed hard X-ray luminosities, unless either the size of the corona is at least of the order of 10^3 Schwarzschild radii, or the corona itself is in fact a reservoir, where the energy is mainly stored in the form of a magnetic field generated by a sheared rotator (probably the accretion disc). We briefly outline the main reasons why the former possibility is to be discarded, and the latter preferred.

The simultaneous presence of a strong quasi periodic oscillation of period of about 10 seconds in the optical and X-ray lightcurves of the X-ray transient XTE J1118+480 suggests that a significant fraction of the optical flux originates from the inner part of the accretion flow, where most of the X-rays are produced. We present a model of magnetic flares in an accretion disc corona where thermal cyclo-synchrotron emission contributes significantly to the optical emission, while the X-rays are produced by inverse Compton scattering of the soft photons produced by dissipation in the underlying disc and by the synchrotron process itself. Given the observational constraints, we estimate the values for the coronal temperature, optical depth and magnetic field intensity, as well as the accretion rate for the source. Within our model we predict a correlation between optical and hard X-ray variability and an anticorrelation between optical and soft X-rays. We also expect optical variability on flaring timescales (about tens of milliseconds), with a power density spectrum similar to the one observed in the X-ray band. Finally we use both the available optical/EUV/X-ray spectral energy distribution and the low frequency time variability to discuss limits on the inner radius of the optically thick disc.

We present a critical analysis of the usual interpretation of the multicolour disc model parameters for black hole candidates in terms of the inner radius and temperature of the accretion disc. Using a self-consistent model for the radiative transfer and the vertical temperature structure in a Shakura-Sunyaev disc, we simulate the observed disc spectra, taking into account doppler blurring and gravitational redshift, and fit them with multicolour models. We show not only that such a model systematically underestimates the value of the inner disc radius, but that when the accretion rate and/or the energy dissipated in the corona are allowed to change the inner edge of the disc, as inferred from the multicolour model, appears to move even when it is in fact fixed at the innermost stable orbit.

We compute exactly the Lense-Thirring precession frequency for point masses in the Kerr metric, for arbitrary black hole mass and specific angular momentum. We show that this frequency, for point masses at or close to the innermost stable orbit, and for holes with moderate to extreme rotation, is less than, but comparable to the rotation frequency. Thus, if the quasi periodic oscillations (QPOs) observed in the modulation of the X-ray flux from some black holes candidates are due to Lense-Thirring precession of orbiting material, we predict that a separate, distinct QPO ought to be observed in each object.