Low-ionisation broad absorption line quasars (LoBALs) mark an important, yet poorly understood, population of quasars showing direct evidence for energetic mass outflows. We outline a sample of 12 luminous (L$_{\rm{bol}}$ $>$10$^{46}$ergs$^{-1}$) LoBALs at 2.0$<$z$<$2.5 - a key epoch in both star formation and black hole accretion, which have been imaged as part of a targeted program with the \textit{Herschel} Spectral and Photometric Imaging REceiver (SPIRE). We present K-band NOTCam spectra for three of these targets, calculating their spectroscopic redshifts, black hole masses and bolometric luminosities, and increasing the total number of LoBAL targets in our sample with spectral information from five to eight. Based on FIR obeservations from \textit{Herschel} SPIRE, we derive prolific SFRs ranging 740 - 2380M$_{\rm{\odot}}$yr$^{-1}$ for the detected targets, consistent with LoBALs existing in an evolutionary phase associated with starburst activity. Furthermore, an upper limit of $<$440M$_{\rm{\odot}}$yr$^{-1}$ is derived for the non-detections, meaning moderate-to-high SFRs cannot be ruled out, even among the undetected targets. Indeed, we detect an enhancement in both the SFRs and FIR fluxes of LoBALs compared to HiBAL and non-BAL quasars, further supporting the evolutionary LoBAL paradigm. Despite this enhancement in SFR however, the environments of LoBALs appear entirely consistent with the general galaxy population at 2.0$<$z$<$2.5.

A series of gravitational instabilities in a circumnuclear gas disk (CND) are required to trigger gas transport to a central supermassive black hole (SMBH) and ignite Active Galactic Nuclei (AGNs). A test of this scenario is to investigate whether an enhanced molecular gas mass surface density ($\Sigma_{\rm mol}$) is found in the CND-scale of quasars relative to a comparison sample of inactive galaxies. Here we performed sub-kpc resolution CO(2-1) observations with ALMA of four low-redshift ($z \sim 0.06$), luminous ($\sim 10^{45}$ erg s$^{-1}$) quasars with each matched to a different star-forming galaxy, having similar redshift, stellar mass, and star-formation rate. We detected CO(2-1) emission from all quasars, which show diverse morphologies. Contrary to expectations, $\Sigma_{\rm mol}$ of the quasar sample, computed from the CO(2-1) luminosity, tends to be smaller than the comparison sample at $r < 500$ pc; there is no systematic enhancement of $\Sigma_{\rm mol}$ in our quasars. We discuss four possible scenarios that would explain the lower molecular gas content (or CO(2-1) luminosity as an actual observable) at the CND-scale of quasars, i.e., AGN-driven outflows, gas-rich minor mergers, time-delay between the onsets of a starburst-phase and a quasar-phase, and X-ray-dominated region (XDR) effects on the gas chemical abundance and excitation. While not extensively discussed in the literature, XDR effects can have an impact on molecular mass measurements particularly in the vicinity of luminous quasar nuclei; therefore higher resolution molecular gas observations, which are now viable using ALMA, need to be considered.

Using a high luminosity ($L_{\rm bol} \sim 10^{47.5} - 10^{48.3}$ erg s$^{-1}$), high redshift ($3.2 < z < 3.8$) quasar sample of 19 quasars with optical and near-infrared spectroscopy, we investigate the reliability of the CIV-based black hole mass estimates ($M_{\rm BH}$). The median logarithm of the CIV- and H$\beta$-based $M_{\rm BH}$ ratios is 0.110 dex with the scatter of 0.647 dex. The CIV-to-H$\beta$ BH mass differences are significantly correlated with the CIV FWHMs, blueshifts and asymmetries. Corrections of the CIV FWHM using the blueshift and asymmetry reduce the scatter of the mass differences by $\sim$ 0.04-0.2 dex. Quasars in our sample accrete at the Eddington ratio $R_{\rm EDD}>0.3$ and cover a considerable range of blueshifts, with 18/19 of the quasars showing CIV blueshifts (with the median value of 1126 km s$^{-1}$) and 14/19 of the quasars showing CIV blueshifts larger than 500 km s$^{-1}$. It suggests that not all quasars with high Eddington ratios show large blueshifts. The Baldwin effect between the CIV rest-frame equivalent width (REW) and the continuum luminosity at 1350 $\overset{\lower.5em\circ}{\mathrm{A}}$ is not seen, likely due to the limited luminosity range of our sample. We find a lack of flux in the red wing of the composite spectrum with larger CIV blueshift, and detect a higher ratio of [OIII] quasars with REW$_{\rm [OIII]}>5$ $\overset{\lower.5em\circ}{\mathrm{A}}$ in the subsample with lower CIV blueshift. It is more likely that they are caused by the combination of the Eddington ratio and the orientation effect.

Correlations between the mass of a supermassive black hole and the properties of its host galaxy (e.g., total stellar mass (M*), luminosity (Lhost)) suggest an evolutionary connection. A powerful test of a co-evolution scenario is to measure the relations MBH-Lhost and MBH-M* at high redshift and compare with local estimates. For this purpose, we acquired HST imaging with WFC3 of 32 X-ray-selected broad-line AGN at 1.2<z<1.7 in deep survey fields. By applying state-of-the-art tools to decompose the HST images including available ACS data, we measured the host galaxy luminosity and stellar mass along with other properties through the 2D model fitting. The black hole mass was determined using the broad Halpha line, detected in the near-infrared with Subaru/FMOS, which potentially minimizes systematic effects using other indicators. We find that the observed ratio of MBH to total M* is 2.7 times larger at z~1.5 than in the local universe, while the scatter is equivalent between the two epochs. A non-evolving mass ratio is consistent with the data at the 2-3 sigma confidence level when accounting for selection effects and their uncertainties. The relationship between MBH-Lhost paints a similar picture. Therefore, our results cannot distinguish whether SMBHs and their total M* and Lhost proceed in lockstep or whether the growth of the former somewhat overshoots the latter, given the uncertainties. Based on a statistical estimate of the bulge-to-total mass fraction, the ratio MBH/M* is offset from the local value by a factor of ~7 which is significant even accounting for selection effects. Taken together, these observations are consistent with a scenario in which stellar mass is subsequently transferred from an angular momentum supported component of the galaxy to the pressure supported one through secular processes or minor mergers at a faster rate than mass accretion onto the SMBH.

The evolution of the galaxy size - stellar mass (Mstellar) relation has been a puzzle for over a decade. High redshift galaxies are significantly more compact than galaxies observed today, at an equivalent mass, but how much of this apparent growth is driven by progenitor bias, minor mergers, secular processes, or feedback from AGN is unclear. To help disentangle the physical mechanisms at work by addressing the latter, we study the galaxy size - Mstellar relation of 32 carefully-selected broad-line AGN hosts at 1.2 < z < 1.7 (7.5 < log M_BH < 8.5; L_bol/L_Edd > 0.1). Using HST with multi-band photometry and state-of-the-art modeling techniques, we measure half-light radii while accounting for uncertainties from subtracting bright central point sources. We find AGN hosts to have sizes ranging from 1 to 6 kpc at Mstellar ~ 0.3 - 1 x 10^11 Msun. Thus, many hosts have intermediate sizes as compared to equal-mass star-forming and quiescent galaxies. While inconsistent with the idea that AGN feedback may induce an increase in galaxy sizes, this finding is consistent with hypotheses in which AGNs preferentially occur in systems with prior concentrated gas reservoirs, or are involved in secular compaction processes perhaps responsible for simultaneously building bulges and shutting down star formation. If driven by minor mergers, which do not grow central black holes as fast as they do bulge-like stellar structures, such a process would explain both the galaxy size - mass relation observed here and the evolution in the black hole, bulge mass relation described in a companion paper.

We report the discovery of 28 quasars and 7 luminous galaxies at 5.7 $\le$ z $\le$ 7.0. This is the tenth in a series of papers from the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, which exploits the deep multi-band imaging data produced by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. The total number of spectroscopically identified objects in SHELLQs has now grown to 93 high-z quasars, 31 high-z luminous galaxies, 16 [O III] emitters at z ~ 0.8, and 65 Galactic cool dwarfs (low-mass stars and brown dwarfs). These objects were found over 900 deg2, surveyed by HSC between 2014 March and 2018 January. The full quasar sample includes 18 objects with very strong and narrow Ly alpha emission, whose stacked spectrum is clearly different from that of other quasars or galaxies. While the stacked spectrum shows N V 1240 emission and resembles that of lower-z narrow-line quasars, the small Ly alpha width may suggest a significant contribution from the host galaxies. Thus these objects may be composites of quasars and star-forming galaxies.

We present ALMA Band~7 observations at $850\mu$m of 20 luminous ($\log\, L_{\rm bol}>46.9$ [erg s$^{-1}$]) unobscured quasars at $z\sim2$. We detect continuum emission for 19/20 quasars. After subtracting an AGN contribution, we measure the total far-IR luminosity for 18 quasars, assuming a modified blackbody model, and attribute the emission as indicative of the star formation rate (SFR). Our sample can be characterized with a log-normal SFR distribution having a mean of 140 $M_\odot$ yr$^{-1}$ and a dispersion of 0.5 dex. Based on an inference of their stellar masses, the SFRs are similar, in both the mean and dispersion, with star-forming main-sequence galaxies at the equivalent epoch. Thus, there is no evidence for a systematic enhancement or suppression (i.e., regulation or quenching) of star formation in the hosts of the most luminous quasars at $z\sim2$. These results are consistent with the Magneticum cosmological simulation, while in disagreement with a widely recognized phenomenological model that predicts higher SFRs than observed here based on the high bolometric luminosities of this sample. Furthermore, there is only a weak relation between SFR and accretion rate onto their supermassive black holes both for average and individual measurements. We interpret these results as indicative of star formation and quasar accretion being fed from the available gas reservoir(s) in their host with a disconnect due to their different physical sizes, temporal scales, and means of gas processing.

We present an unprecedented view on the morphology and kinematics of the extended narrow-line region (ENLR) and molecular gas around the prototypical hyper-luminous quasar 3C273 ($L\sim10^{47}$ erg/s at z=0.158) based on VLT-MUSE optical 3D spectroscopy and ALMA observations. We find that: 1) The ENLR size of 12.1$\pm$0.2kpc implies a smooth continuation of the size-luminosity relation out to large radii or a much larger break radius as previously proposed. 2) The kinematically disturbed ionized gas with line splits reaching 1000km/s out to 6.1$\pm$1.5kpc is aligned along the jet axis. 3) The extreme line broadening on kpc scales is caused by spatial and spectral blending of many distinct gas clouds separated on sub-arcsecond scales with different line-of-sight velocities. The ENLR velocity field combined with the known jet orientation rule out a simple scenario of a radiatively-driven radial expansion of the outflow. Instead we propose that a pressurized expanding hot gas cocoon created by the radio jet is impacting on an inclined gas disk leading to transverse and/or backflow motion with respect to our line-of-sight. The molecular gas morphology may either be explained by a density wave at the front of the outflow expanding along the jet direction as predicted by positive feedback scenario or the cold gas may be trapped in a stellar over-density caused by a recent merger event. Using 3C273 as a template for observations of high-redshift hyper-luminous AGN reveals that large-scale ENLRs and kpc scale outflows may often be missed due to the brightness of the nuclei and the limited sensitivity of current near-IR instrumentation.

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.

Research over the past decade has shown diminishing empirical evidence for major galaxy mergers being a dominating or even important mechanism for the growth of supermassive black holes in galaxies and the triggering of optically or X-ray selected active galactic nuclei (AGN). We here for the first time test whether such a connection exists at least in the most plausible part of parameter space for this mechanism: the highest specific accretion rate broad-line AGNs at the peak epoch of black hole activity around z = 2. To that end we examine 21 galaxies hosting a high accreting black hole (L/Ledd > 0.7) observed with HST/WFC3 and 92 stellar mass- and redshift- matched inactive galaxies taken from the CANDELS survey. We removed the AGN point sources from their host galaxies and avoided bias in visual classification by adding and then subtracting mock point sources to and from the comparison galaxies, producing matched residual structures for both sets. The resulting samples were joined, randomized, and subsequently visually ranked with respect to perceived strength of structural distortions by 10 experts. The ensuing individual rankings were combined into a consensus sequence and from this we derived merger fractions for both samples. With the merger fractions f$_{m,agn}$ = 0.24 $\pm$ 0.09 for the AGN host galaxy sample and f$_{m,ina}$ = 0.19 $\pm$ 0.04 for the inactive galaxies, we find no significant difference between the AGN host galaxies and inactive galaxies. Also, both samples display comparable fractions of disk-dominated galaxies. These findings are consistent with previous studies for different AGN populations, and we conclude that even black hole growth at the highest specific accretion rates and at the peak of cosmic AGN activity is not predominantly caused by major mergers. (abriged)

We investigate the multi-wavelength properties of host galaxies of 3701 X-ray-selected active galactic nuclei (AGNs) out to z~5 in the Chandra-COSMOS Legacy Survey. Thanks to the extensive multi-wavelength photometry available in the COSMOS field, we derive AGN luminosities, host stellar masses, and star formation rates (SFRs) via a multi-component SED fitting technique. Type 1 and Type 2 AGNs follow the same intrinsic Lx-L6um relation, suggesting that mid-infrared emission is a reasonably good measure of the AGN accretion power regardless of obscuration. We find that there is a strong increase in Type 1 AGN fraction toward higher AGN luminosity, possibly due to the fact that Type 1 AGNs tend to be hosted by more massive galaxies. The AGN luminosity and SFR are consistent with an increase toward high stellar mass, while both the Mstellar-dependence is weaker towards the high-mass end, which could be interpreted as a consequence of quenching both star formation and AGN activity in massive galaxies. AGN host galaxies tend to have SFRs that are consistent with normal star-forming galaxies, independent of AGN luminosities. We confirm that black hole accretion rate and SFR are correlated up to z~5, when forming stars. The majority (~73%) of our AGN sample are faint in the far-infrared, implying that the moderate-luminosity AGNs seem to be still active after the star formation is suppressed. It is not certain whether AGN activity plays a role in quenching the star formation. We conclude that both AGN activity and star formation might be more fundamentally related to host stellar mass.

We report the discovery of a quasar at z = 7.07, which was selected from the deep multi-band imaging data collected by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. This quasar, HSC J124353.93+010038.5, has an order of magnitude lower luminosity than do the other known quasars at z > 7. The rest-frame ultraviolet absolute magnitude is M1450 = -24.13 +/- 0.08 mag and the bolometric luminosity is Lbol = (1.4 +/- 0.1) x 10^{46} erg/s. Its spectrum in the optical to near-infrared shows strong emission lines, and shows evidence for a fast gas outflow, as the C IV line is blueshifted and there is indication of broad absorption lines. The Mg II-based black hole mass is Mbh = (3.3 +/- 2.0) x 10^8 Msun, thus indicating a moderate mass accretion rate with an Eddington ratio 0.34 +/- 0.20. It is the first z > 7 quasar with sub-Eddington accretion, besides being the third most distant quasar, known to date. The luminosity and black hole mass are comparable to, or even lower than, those measured for the majority of low-z quasars discovered by the Sloan Digital Sky Survey, and thus this quasar likely represents a z > 7 counterpart to quasars commonly observed in the low-z universe.

We present new measurements of the quasar luminosity function (LF) at $z \sim 6$, over an unprecedentedly wide range of the rest-frame ultraviolet luminosity $M_{1450}$ from $-30$ to $-22$ mag. This is the fifth in a series of publications from the Subaru High-$z$ Exploration of Low-Luminosity Quasars (SHELLQs) project, which exploits the deep multi-band imaging data produced by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. The LF was calculated with a complete sample of 110 quasars at $5.7 \le z \le 6.5$, which includes 48 SHELLQs quasars discovered over 650 deg$^2$, and 63 brighter quasars discovered by the Sloan Digital Sky Survey and the Canada-France-Hawaii Quasar Survey (including one overlapping object). This is the largest sample of $z \sim 6$ quasars with a well-defined selection function constructed to date, and has allowed us to detect significant flattening of the LF at its faint end. A double power-law function fit to the sample yields a faint-end slope $\alpha = -1.23^{+0.44}_{-0.34}$, a bright-end slope $\beta = -2.73^{+0.23}_{-0.31}$, a break magnitude $M_{1450}^* = -24.90^{+0.75}_{-0.90}$, and a characteristic space density $\Phi^* = 10.9^{+10.0}_{-6.8}$ Gpc$^{-3}$ mag$^{-1}$. Integrating this best-fit model over the range $-18 < M_{1450} < -30$ mag, quasars emit ionizing photons at the rate of $\dot{n}_{\rm ion} = 10^{48.8 \pm 0.1}$ s$^{-1}$ Mpc$^{-3}$ at $z = 6.0$. This is less than 10 % of the critical rate necessary to keep the intergalactic medium ionized, which indicates that quasars are not a major contributor to cosmic reionization.

We present near-IR spectroscopy in J- and H-band for a large sample of 243 X-ray selected moderate-luminosity type-1 AGN in the COSMOS, SXDS and E-CDF-S survey fields using the multi-object spectrograph Subaru/FMOS. Our sample covers the redshift range 0.5<z<3.0 and an X-ray luminosity range of $10^{43}\lesssim L_X \lesssim 10^{45}$~erg s$^{-1}$. We provide emission-line properties and derived virial black hole mass estimates, bolometric luminosities and Eddington ratios, based on H$\alpha$ (211), H$\beta$ (63) and MgII (4). We compare line widths, luminosities and black hole mass estimates from H$\alpha$ and H$\beta$ and augment these with commensurate measurements of MgII and CIV detected in optical spectra. We demonstrate the robustness of using H$\alpha$, H$\beta$ and MgII as reliable black hole mass estimators for high-z moderate-luminosity AGN, while the use of CIV is prone to large uncertainties (>0.4 dex). We extend a recently proposed correction based on the CIV blueshift to lower luminosities and black hole masses. While our sample shows an improvement in their CIV black hole mass estimates, the deficit of high blueshift sources reduces its overall importance for moderate-luminosity AGN, compared to the most luminous quasars. In addition, we revisit luminosity correlations between $L_{\rm{bol}}$, $L_X$, $L_{\rm{[OIII]}}$, $L_{5100}$ and $L_{\rm{H}\alpha}$ and find them to be consistent with a simple empirical model, based on a small number of well-established scaling relations. We finally highlight our highest redshift AGN, CID 781, at z=4.6 which shows the lowest black hole mass ($\sim10^8$M$_\odot$) among current near-IR samples at this redshift, and is in a state of fast growth.

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 discovery of 41 new high-z quasars and luminous galaxies, which were spectroscopically identified at 5.7 < z < 6.9. This is the fourth in a series of papers from the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, based on the deep multi-band imaging data collected by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. We selected the photometric candidates by a Bayesian probabilistic algorithm, and then carried out follow-up spectroscopy with the Gran Telescopio Canarias and the Subaru Telescope. Combined with the sample presented in the previous papers, we have now spectroscopically identified 137 extremely-red HSC sources over about 650 deg2, which include 64 high-z quasars, 24 high-z luminous galaxies, 6 [O III] emitters at z ~ 0.8, and 43 Galactic cool dwarfs (low-mass stars and brown dwarfs). The new quasars span the luminosity range from M1450 ~ -26 to -22 mag, and continue to populate a few magnitude lower luminosities than have been probed by previous wide-field surveys. In a companion paper, we derive the quasar luminosity function at z ~ 6 over an unprecedentedly wide range of M1450 ~ -28 to -21 mag, exploiting the SHELLQs and other survey outcomes.

We present the discovery of strong Balmer line absorption in H$\alpha$ to H$\gamma$ in two luminous low-ionization broad absorption line quasars (LoBAL QSOs) at z~1.5, with black hole masses around $10^{10}$ $M_\odot$ from near-IR spectroscopy. There are only two previously known quasars at z>1 showing Balmer line absorption. SDSS J1019+0225 shows blueshifted absorption by ~1400 km/s with an H$\alpha$ rest-frame equivalent width of 13 \AA{}. In SDSS J0859+4239 we find redshifted absorption by ~500 km/s with an H$\alpha$ rest-frame equivalent width of 7 \AA{}. The redshifted absorption could indicate an inflow of high density gas onto the black hole, though we cannot rule out alternative interpretations. The Balmer line absorption in both objects appears to be saturated, indicating partial coverage of the background source by the absorber. We estimate the covering fractions and optical depth of the absorber and derive neutral hydrogen column densities, $N_{\rm{H\,I}}\sim1.3\times 10^{18}$ cm$^{-2}$ for SDSS J1019+0225 and $N_{\rm{H\,I}}\sim9\times 10^{17}$ cm$^{-2}$ for SDSS J0859+4239, respectively. In addition, the optical spectra reveal also absorption troughs in HeI* $\lambda3889$ and $\lambda3189$ in both objects.

We present near-IR spectroscopy of 22 luminous low-ionization broad absorption line quasars (LoBAL QSOs) at redshift 1.3<z<2.5, with 12 objects at z~1.5 and 10 at z~2.3. The spectra cover the rest-frame H$\alpha$ and H$\beta$ line regions, allowing us to obtain robust black hole mass estimates based on the broad H$\alpha$ line. We use these data, augmented by a lower-redshift sample from the Sloan Digital Sky Survey, to test the proposed youth scenario for LoBALs, which suggests that LoBALs constitute an early short-lived evolutionary stage of quasar activity, by probing for any difference in their masses, Eddington ratios, or rest-frame optical spectroscopic properties compared to normal quasars. In addition, we construct the UV to mid-IR spectral energy distributions (SEDs) for the LoBAL sample and a matched non-BAL quasar sample. We do not find any statistically significant difference between LoBAL QSOs and non-BAL QSOs in their black hole mass or Eddington ratio distributions. The mean UV to mid-IR SED of the LoBAL QSOs is consistent with non-BAL QSOs, apart from their stronger reddening. At z>1 there is no clear difference in their optical emission line properties. We do not see particularly weak [OIII] or strong FeII emission. The LoBAL QSOs do not show a stronger prevalence of ionized gas outflows as traced by the [OIII] line, compared to normal QSOs of similar luminosity. We conclude that the optical-MIR properties of LoBAL QSOs are consistent with the general quasar population and do not support them to constitute a special phase of active galactic nucleus evolution.

One of the major unsolved questions on the understanding of the AGN population is the origin of the dichotomy between radio-quiet and radio-loud quasars. The most promising explanation is provided by the spin paradigm, which suggests radio-loud quasars have higher black hole spin. However, the measurement of black hole spin remains extremely challenging. We here aim at comparing the mean radiative efficiencies of carefully matched samples of radio-loud and radio-quiet SDSS quasars at 0.3<z<0.8. We use the [OIII] luminosity as an indirect average tracer of the ionizing continuum in the extreme-UV regime where differences in the SED due to black hole spin are most pronounced. We find that the radio-loud sample shows an enhancement in [OIII] line strength by a factor of at least 1.5 compared to a radio-quiet sample matched in redshift, black hole mass and optical continuum luminosity or accretion rate. We argue that this enhancement is caused by differences in the SED, suggesting higher average bolometric luminosities at fixed accretion rate in the radio-loud population. This suggests that the radio-loud quasar population has on average systematically larger radiative efficiencies and therefore higher black hole spin than the radio-quiet population, providing observational support for the black hole spin paradigm.

We present spectroscopic identification of 32 new quasars and luminous galaxies discovered at 5.7 < z < 6.8. This is the second in a series of papers presenting the results of the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, which exploits the deep multi-band imaging data produced by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. The photometric candidates were selected by a Bayesian probabilistic algorithm, and then observed with spectrographs on the Gran Telescopio Canarias and the Subaru Telescope. Combined with the sample presented in the previous paper, we have now identified 64 HSC sources over about 430 deg2, which include 33 high-z quasars, 14 high-z luminous galaxies, 2 [O III] emitters at z ~ 0.8, and 15 Galactic brown dwarfs. The new quasars have considerably lower luminosity (M1450 ~ -25 to -22 mag) than most of the previously known high-z quasars. Several of these quasars have luminous (> 10^(43) erg/s) and narrow (< 500 km/s) Ly alpha lines, and also a possible mini broad absorption line system of N V 1240 in the composite spectrum, which clearly separate them from typical quasars. On the other hand, the high-z galaxies have extremely high luminosity (M1450 ~ -24 to -22 mag) compared to other galaxies found at similar redshift. With the discovery of these new classes of objects, we are opening up new parameter spaces in the high-z Universe. Further survey observations and follow-up studies of the identified objects, including the construction of the quasar luminosity function at z ~ 6, are ongoing.

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

We examine the clustering of quasars over a wide luminosity range, by utilizing 901 quasars at $\overline{z}_{\rm phot}\sim3.8$ with $-24.73<M_{\rm 1450}<-22.23$ photometrically selected from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) S16A Wide2 date release and 342 more luminous quasars at $3.4<z_{\rm spec}<4.6$ having $-28.0<M_{\rm 1450}<-23.95$ from the Sloan Digital Sky Survey (SDSS) that fall in the HSC survey fields. We measure the bias factors of two quasar samples by evaluating the cross-correlation functions (CCFs) between the quasar samples and 25790 bright $z\sim4$ Lyman Break Galaxies (LBGs) in $M_{\rm 1450}<-21.25$ photometrically selected from the HSC dataset. Over an angular scale of \timeform{10.0"} to \timeform{1000.0"}, the bias factors are $5.93^{+1.34}_{-1.43}$ and $2.73^{+2.44}_{-2.55}$ for the low and high luminosity quasars, respectively, indicating no luminosity dependence of quasar clustering at $z\sim4$. It is noted that the bias factor of the luminous quasars estimated by the CCF is smaller than that estimated by the auto-correlation function (ACF) over a similar redshift range, especially on scales below \timeform{40.0"}. Moreover, the bias factor of the less-luminous quasars implies the minimal mass of their host dark matter halos (DMHs) is $0.3$-$2\times10^{12}h^{-1}M_{\odot}$, corresponding to a quasar duty cycle of $0.001$-$0.06$.

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

We present the luminosity function of z=4 quasars based on the Hyper Suprime-Cam Subaru Strategic Program Wide layer imaging data in the g, r, i, z, and y bands covering 339.8 deg^2. From stellar objects, 1666 z~4 quasar candidates are selected by the g-dropout selection down to i=24.0 mag. Their photometric redshifts cover the redshift range between 3.6 and 4.3 with an average of 3.9. In combination with the quasar sample from the Sloan Digital Sky Survey in the same redshift range, the quasar luminosity function covering the wide luminosity range of M1450=-22 to -29 mag is constructed. It is well described by a double power-law model with a knee at M1450=-25.36+-0.13 mag and a flat faint-end slope with a power-law index of -1.30+-0.05. The knee and faint-end slope show no clear evidence of redshift evolution from those at z~2. The flat slope implies that the UV luminosity density of the quasar population is dominated by the quasars around the knee, and does not support the steeper faint-end slope at higher redshifts reported at z>5. If we convert the M1450 luminosity function to the hard X-ray 2-10keV luminosity function using the relation between UV and X-ray luminosity of quasars and its scatter, the number density of UV-selected quasars matches well with that of the X-ray-selected AGNs above the knee of the luminosity function. Below the knee, the UV-selected quasars show a deficiency compared to the hard X-ray luminosity function. The deficiency can be explained by the lack of obscured AGNs among the UV-selected quasars.

We report the serendipitous discovery of a quadruply lensed source at $z_{\rm s}=3.76$, HSC~J115252+004733, from the Hyper Suprime-Cam (HSC) Survey. The source is lensed by an early-type galaxy at $z_{\rm l}=0.466$ and a satellite galaxy. Here, we investigate the properties of the source by studying its size and luminosity from the imaging and the luminosity and velocity width of the Ly-$\alpha$ line from the spectrum. Our analyses suggest that the source is most probably a low-luminosity active galactic nucleus (LLAGN) but the possibility of it being a compact bright galaxy (e.g., a Lyman-$\alpha$ emitter or Lyman Break Galaxy) cannot be excluded. The brighter pair of lensed images appears point-like except in the HSC $i$-band (with a seeing $\sim0.5"$). The extended emission in the $i$-band image could be due to the host galaxy underneath the AGN, or alternatively, due to a highly compact lensed galaxy (without AGN) which appears point-like in all bands except in $i$-band. We also find that the flux ratio of the brighter pair of images is different in the Ks-band compared to optical wavelengths. Phenomena such as differential extinction and intrinsic variability cannot explain this chromatic variation. While microlensing from stars in the foreground galaxy is less likely to be the cause, it cannot be ruled out completely. If the galaxy hosts an AGN, then this represents the highest redshift quadruply imaged AGN known to date, enabling study of a distant LLAGN. Discovery of this unusually compact and faint source demonstrates the potential of the HSC survey.

The most frequently proposed model for the origin of quasars holds that the high accretion rates seen in luminous active galactic nuclei are primarily triggered during major mergers between gas-rich galaxies. While plausible for decades, this model has only begun to be tested with statistical rigor in the past few years. Here we report on a Hubble Space Telescope study to test this hypothesis for $z=2$ quasars with high super-massive black hole masses ($M_\mathrm{BH}=10^9-10^{10}~M_\odot{}$), which dominate cosmic black hole growth at this redshift. We compare Wide Field Camera 3 $F160W$ (rest-frame $V$-band) imaging of 19 point source-subtracted quasar hosts to a matched sample of 84 inactive galaxies, testing whether the quasar hosts have greater evidence for strong gravitational interactions. Using an expert ranking procedure, we find that the quasar hosts are uniformly distributed within the merger sequence of inactive galaxies, with no preference for quasars in high-distortion hosts. Using a merger/non-merger cutoff approach, we recover distortion fractions of $f_\mathrm{m,qso}=0.39\pm{}0.11$ for quasar hosts and $f_\mathrm{m,gal}=0.30\pm{}0.05$ for inactive galaxies (distribution modes, 68% confidence intervals), with both measurements subjected to the same observational conditions and limitations. The slight enhancement in distorted fraction for quasar hosts over inactive galaxies is not significant, with a probability that the quasar fraction is higher of $P(f_\mathrm{m,qso}>f_\mathrm{m,gal}) = 0.78$ ($0.78\,\sigma{}$), in line with results for lower mass and lower $z$ AGN. We find no evidence that major mergers are the primary triggering mechanism for the massive quasars that dominate accretion at the peak of cosmic quasar activity.

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.

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

Understanding the relationship between galaxies hosting active galactic nuclei (AGN) and the dark matter halos in which they reside is key to constraining how black-hole fueling is triggered and regulated. Previous efforts have relied on simple halo mass estimates inferred from clustering, weak gravitational lensing, or halo occupation distribution modeling. In practice, these approaches remain uncertain because AGN, no matter how they are identified, potentially live a wide range of halo masses with an occupation function whose general shape and normalization are poorly known. In this work, we show that better constraints can be achieved through a rigorous comparison of the clustering, lensing, and cross-correlation signals of AGN hosts to a fiducial stellar-to-halo mass relation (SHMR) derived for all galaxies. Our technique exploits the fact that the global SHMR can be measured with much higher accuracy than any statistic derived from AGN samples alone. Using 382 moderate luminosity X-ray AGN at z<1 from the COSMOS field, we report the first measurements of weak gravitational lensing from an X-ray selected sample. Comparing this signal to predictions from the global SHMR, we find that, contrary to previous results, most X-ray AGN do not live in medium size groups ---nearly half reside in relatively low mass halos with Mh~10^12.5 Msun. The AGN occupation function is well described by the same form derived for all galaxies but with a lower normalization---the fraction of halos with AGN in our sample is a few percent. By highlighting the relatively "normal" way in which moderate luminosity X-ray AGN hosts occupy halos, our results suggest that the environmental signature of distinct fueling modes for luminous QSOs compared to moderate luminosity X-ray AGN is less obvious than previously claimed.

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.

The redshift evolution of the black hole - bulge relations is an essential observational constraint for models of black hole - galaxy coevolution. In addition to the observational challenges for these studies, conclusions are complicated by the influence of selection effects. We demonstrate that there is presently no statistical significant evidence for cosmological evolution in the black hole-bulge relations, once these selection effects are taken into account and corrected for. We present a fitting method, based on the bivariate distribution of black hole mass and galaxy property, that accounts for the selection function in the fitting and is therefore able to recover the intrinsic black hole - bulge relation unbiased. While prior knowledge is restricted to a minimum, we at least require knowledge of either the sample selection function and the mass dependence of the active fraction, or the spheroid distribution function and the intrinsic scatter in the black hole - bulge relation. We employed our fitting routine to existing studies of the black hole-bulge relation at z~1.5 and z~6, using our current best knowledge of the distribution functions. There is no statistical significant evidence for positive evolution in the MBH-M* ratio out to z~2. At z~6 the current constraints are less strong, but we demonstrate that the large observed apparent offset from the local black hole-bulge relation at z~6 is fully consistent with no intrinsic offset. The method outlined here provides a tool to obtain more reliable constraints on black hole - galaxy co-evolution in the future.

To investigate the validity of the assumption that quiescent galaxies and active galaxies follow the same black hole mass - stellar velocity dispersion relation, as required for the calibration of M_BH estimators for broad line AGNs, we determine and compare the M-sigma relations, respectively, for quiescent and active galaxies. For the quiescent galaxy sample, composed of 72 dynamical M_BH measurements, we update sigma for 28 galaxies using homogeneous H-band measurements that are corrected for galaxy rotation. For active galaxies, we collect 25 reverberation-mapped AGNs and improve sigma measurement for two objects . Combining the two samples, we determine the virial factor $f$, first by scaling the active galaxy sample to the M-sigma relation of quiescent galaxies, and second by simultaneously fitting the quiescent and active galaxy samples, as f=5.1_{-1.1}^{+1.5} and f=5.9_{-1.5}^{+2.1}, respectively. The M-sigma relation of active galaxies appears to be shallower than that of quiescent galaxies. However, the discrepancy is caused by a difference in the accessible M_BH distribution at given sigma, primarily due to the difficulty of measuring reliable stellar velocity dispersion for the host galaxies of luminous AGNs. Accounting for the selection effects, we find that active and quiescent galaxies are consistent with following intrinsically the same M-sigma relation.

To calibrate stellar velocity dispersion measurements from optical and near-IR stellar lines, and to improve the black hole mass (MBH)- stellar velocity dispersion (sigma) relation, we measure stellar velocity dispersions based on high quality H-band spectra for a sample of 31 nearby galaxies, for which dynamical MBH is available in the literature. By comparing velocity dispersions measured from stellar lines in the H-band with those measured from optical stellar lines, we find no significant difference, suggesting that optical and near-IR stellar lines represent the same kinematics and that dust effect is negligible for early-type galaxies. Based on the spatially-resolved rotation and velocity dispersion measurements along the major axis of each galaxy, we find that a rotating stellar disk is present for 80% of galaxies in the sample. For galaxies with a rotation component, velocity dispersions measured from a single aperture spectrum can vary by up to ~20%, depending on the size of the adopted extraction aperture. To correct for the rotational broadening, we derive luminosity-weighted velocity dispersion within the effective radius of each galaxy, providing uniformly measured velocity dispersions to improve the MBH-sigma relation.

We present an investigation of sample selection effects that influence the observed black hole - bulge relations and its evolution with redshift. We provide a common framework in which all kinds of selection effects on the BH-bulge relations can be investigated, but our main emphasis is on the consequences of using broad-line AGN and their host galaxies to search for evolution in the BH-bulge relation. We identified relevant sources of bias that were not discussed in the literature so far. A particularly important effect is caused by the fact that the active fraction among SMBHs varies considerably with BH mass, in the sense that high-mass BHs are less likely to be active than lower mass ones. In the connection with intrinsic scatter of the BH-bulge relation this effect implies a bias towards a low BH mass at given bulge property. This effect adds to the bias caused by working with luminosity or flux limited samples that were already discussed by others. A quantitative prediction of these biases requires (i) a realistic model of the sample selection function, and (ii) knowledge of relevant underlying distribution functions. For low-redshift AGN samples we can naturally reproduce the flattening of the relation observed in some studies. When extending our analysis to higher redshift samples we are clearly hampered by limited empirical constraints on the various relevant distribution functions. Using a best-guess approach for these distributions we estimate the expected magnitude of sample selection biases for a number of recent observational attempts to study the BH-bulge evolution. In no case do we find statistically significant evidence for an evolving BH-bulge relation. We suggest a possible practical approach to circumvent several of the most problematic issues connected with AGN selection; this could become a powerful diagnostic in future investigations (abridged).

Stellar dynamical modeling is a powerful method to determine the mass of black holes in quiescent galaxies. However, in previous work the presence of a dark matter halo has been ignored in the modeling. Gebhardt & Thomas (2009) showed that accounting for a dark matter halo increased the black hole mass of the massive galaxy M87 by a factor of two. We used a sample of 12 galaxies to investigate the effect of accounting for a dark matter halo in the dynamical modeling in more detail, and also updated the masses using improved modeling. The sample of galaxies possesses Hubble Space Telescope and ground based observations of stellar kinematics. Their black hole masses have been presented before, but without including a dark matter halo in the models. Without a dark halo, we find a mean increase in the estimated mass of 1.5 for the whole sample compared to previous results. We attribute this change to using a more complete orbit library. When we include a dark matter halo, along with the updated models, we find an additional increase in black hole mass by a factor of 1.2 in the mean, much less than for M87. We attribute the smaller discrepancy in black hole mass to using data that better resolves the black hole's sphere of influence. We redetermined the M-sigma and M-L relationships using our updated black hole masses and found a slight increase in both normalization and intrinsic scatter.

We estimated black hole masses and Eddington ratios for a well defined sample of local (z<0.3) broad line AGN from the Hamburg/ESO Survey (HES), based on the Hbeta line and standard recipes assuming virial equilibrium for the broad line region. The sample represents the low-redshift AGN population over a wide range of luminosities, from Seyfert 1 galaxies to luminous quasars. From the distribution of black hole masses we derived the active black hole mass function (BHMF) and the Eddington ratio distribution function (ERDF) in the local universe, exploiting the fact that the HES has a well-defined selection function. While the directly determined ERDF turns over around L/L_Edd ~ 0.1, similar to what has been seen in previous analyses, we argue that this is an artefact of the sample selection. We employed a maximum likelihood approach to estimate the intrinsic distribution functions of black hole masses and Eddington ratios simultaneously in an unbiased way, taking the sample selection function fully into account. The resulting ERDF is well described by a Schechter function, with evidence for a steady increase towards lower Eddington ratios, qualitatively similar to what has been found for type~2 AGN from the SDSS. Comparing our best-fit active BHMF with the mass function of inactive black holes we obtained an estimate of the fraction of active black holes, i.e. an estimate of the AGN duty cycle. The active fraction decreases strongly with increasing black hole mass. A comparison with the BHMF at higher redshifts also indicates that, at the high mass end, black holes are now in a less active stage than at earlier cosmic epochs. Our results support the notion of anti-hierarchical growth of black holes, and are consistent with a picture where the most massive black holes grew at early cosmic times, whereas at present mainly smaller mass black holes accrete at a significant rate.

We present a determination of the local (z=0) luminosity function of optically selected type 1 (broad-line) Active Galactic Nuclei. Our primary resource is the Hamburg/ESO Survey (HES), which provides a well-defined sample of more than 300 optically bright AGN with redshifts z<0.3 and blue magnitudes B<17.5. AGN luminosities were estimated in two ways, always taking care to minimise photometric biases due to host galaxy light contamination. Firstly, we measured broad-band B_J (blue) magnitudes of the objects over small apertures of the size of the seeing disk. Secondly, we extracted H alpha and H beta broad emission line luminosities from the spectra which should be entirely free of any starlight contribution. The resulting AGN luminosity function (AGNLF) is consistent with a single power law, also when considering the effects of number density evolution within the narrow redshift range. We compared our AGNLF with the H alpha luminosity function of lower luminosity Seyfert 1 galaxies by Hao et al. (2005) and found a smooth transition between both, with excellent agreement in the overlapping region. From the combination of HES and SDSS samples we constructed a single local AGNLF spanning more than 4 orders of magnitude in luminosity. It shows only mild curvature which can be well described as a double power law with slope indices of -2.0 for the faint end and -2.8 for the bright end. We predicted the local AGNLF in the soft X-ray domain and compared this to recent literature data. We also compared the local AGNLF with results obtained at higher redshifts and find strong evidence for luminosity-dependent evolution, in the sense that AGN with luminosities around M_B~-19 are as common in the local universe as they were at z=1.5, supporting the 'AGN downsizing' picture (abridged).