Theory predicts that a supermassive black hole binary (SMBHB) could be observed as a luminous active galactic nucleus (AGN) that periodically varies on the order of its orbital timescale. In X-rays, periodic variations could be caused by mechanisms including relativistic Doppler boosting and shocks. Here we present the first systematic search for periodic AGNs using $941$ hard X-ray light curves (14-195 keV) from the first 105 months of the Swift Burst Alert Telescope (BAT) survey (2004-2013). We do not find evidence for periodic AGNs in Swift-BAT, including the previously reported SMBHB candidate MCG+11$-$11$-$032. We find that the null detection is consistent with the combination of the upper-limit binary population in AGNs in our adopted model, their expected periodic variability amplitudes, and the BAT survey characteristics. We have also investigated the detectability of SMBHBs against normal AGN X-ray variability in the context of the eROSITA survey. Under our assumptions of a binary population and the periodic signals they produce which have long periods of hundreds of days, up to $13$% true periodic binaries can be robustly distinguished from normal variable AGNs with the ideal uniform sampling. However, we demonstrate that realistic eROSITA sampling is likely to be insensitive to long-period binaries because longer observing gaps reduce their detectability. In contrast, large observing gaps do not diminish the prospect of detecting binaries of short, few-day periods, as 19% can be successfully recovered, the vast majority of which can be identified by the first half of the survey.

We present ALMA band-7 data of the [CII] $\lambda157.74\,\mu{\rm m}$ emission line and underlying far-infrared (FIR) continuum for twelve luminous quasars at $z \simeq 4.8$, powered by fast-growing supermassive black holes (SMBHs). Our total sample consists of eighteen quasars, twelve of which are presented here for the first time. The new sources consists of six Herschel/SPIRE detected systems, which we define as "FIR-bright" sources, and six Herschel/SPIRE undetected systems, which we define as "FIR-faint" sources. We determine dust masses for the quasars hosts of $M_{dust} \le 0.2-25.0\times 10^8 M_{\odot}$, implying ISM gas masses comparable to the dynamical masses derived from the [CII] kinematics. It is found that on average the MgII line is blueshifted by $\sim 500\,{\rm km\,s}^{-1}$ with respect to the [CII] emission line, which is also observed when complementing our observations with data from the literature. We find that all of our "FIR-bright" subsample and most of the "FIR-faint" objects lie above the main sequence of star forming galaxies at $z \sim 5$. We detect companion sub-millimeter galaxies (SMGs) for two sources, both FIR-faint, with a range of projected distances of $\sim20-60$ kpc and with typical velocity shifts of $\left|\Delta v\right| \lesssim200\,{\rm km\,s}^{-1}$ from the quasar hosts. Of our total sample of eighteen quasars, 5/18 are found to have dust obscured starforming companions.

I review the current understanding of some key properties of the earliest growing supermassive black holes (SMBHs), as determined from the most up-to-date observations of z>=5 quasars. This includes their accretion rates and growth history, their host galaxies, and the large-scale environments that enabled their emergence less than a billion years after the Big Bang. The available multi-wavelength data show that these SMBHs are consistent with Eddington-limited, radiatively efficient accretion that had to proceed almost continuously since very early epochs. ALMA observations of the hosts' ISM reveal gas-rich, well developed galaxies, with a wide range of SFRs that may exceed ~1000 M_sol/yr. Moreover, ALMA uncovers a high fraction of companion, interacting galaxies, separated by <100 kpc (projected). This supports the idea that the first generation of high-mass, luminous SMBHs grew in over-dense environments, and that major mergers may be important drivers for rapid SMBH and host galaxy growth. Current X-ray surveys cannot access the lower-mass, supposedly more abundant counterparts of these rare z>5 massive quasars, which should be able to elucidate the earliest stages of BH formation and growth. Such lower-mass nuclear BHs will be the prime targets of the deepest surveys planned for the next generation of facilities, such as the upcoming Athena mission and the future Lynx mission concept.

Accretion onto black holes at rates above the Eddington limit has long been discussed in the context of supermassive black hole (SMBH) formation and evolution, providing a possible explanation for the presence of massive quasars at high redshifts (z$\gtrsim$7), as well as having implications for SMBH growth at later epochs. However, it is currently unclear whether such `super-Eddington' accretion occurs in SMBHs at all, how common it is, or whether every SMBH may experience it. In this work, we investigate the observational consequences of a simplistic model for super-Eddington accretion flows -- an optically thick, geometrically thin accretion disc (AD) where the inner-most parts experience severe photon-trapping, which is enhanced with increased accretion rate. The resulting spectral energy distributions (SEDs) show a dramatic lack of rest-frame UV, or even optical, photons. Using a grid of model SEDs spanning a wide range in parameter space (including SMBH mass and accretion rate), we find that large optical quasar surveys (such as SDSS) may be missing most of these luminous systems. We then propose a set of colour selection criteria across optical and infra-red colour spaces designed to select super-Eddington SEDs in both wide-field surveys (e.g., using SDSS, 2MASS and WISE) and deep & narrow-field surveys (e.g., COSMOS). The proposed selection criteria are a necessary first step in establishing the relevance of advection-affected super-Eddington accretion onto SMBHs at early cosmic epochs.

We have conducted 22 GHz radio imaging at 1" resolution of 100 low-redshift AGN selected at 14-195 keV by the Swift-BAT. We find a radio core detection fraction of 96%, much higher than lower-frequency radio surveys. Of the 96 radio-detected AGN, 55 have compact morphologies, 30 have morphologies consistent with nuclear star formation, and 11 have sub-kpc to kpc-scale jets. We find that the total radio power does not distinguish between nuclear star formation and jets as the origin of the radio emission. For 87 objects, we use optical spectroscopy to test whether AGN physical parameters are distinct between radio morphological types. We find that X-ray luminosities tend to be higher if the 22 GHz morphology is jet-like, but find no significant difference in other physical parameters. We find that the relationship between the X-ray and core radio luminosities is consistent with the $L_R/L_X \sim 10^{-5}$ of coronally active stars. We further find that the canonical fundamental planes of black hole activity systematically over-predict our radio luminosities, particularly for objects with star formation morphologies.

We investigate the cosmic evolution of the ratio between black hole mass (MBH) and host galaxy total stellar mass (Mstellar) out to z~2.5 for a sample of 100 X-ray-selected moderate-luminosity, broad-line active galactic nuclei (AGNs) in the Chandra-COSMOS Legacy Survey. By taking advantage of the deep multi-wavelength photometry and spectroscopy in the COSMOS field, we measure in a uniform way the galaxy total stellar mass using a SED decomposition technique and the black hole mass based on broad emission line measurements and single-epoch virial estimates. Our sample of AGN host galaxies has total stellar masses of 10^10-12Msun, and black hole masses of 10^7.0-9.5Msun. Combining our sample with the relatively bright AGN samples from the literature, we find no significant evolution of the MBH-Mstellar relation with black hole-to-host total stellar mass ratio of MBH/Mstellar~0.3% at all redshifts probed. We conclude that the average black hole-to-host stellar mass ratio appears to be consistent with the local value within the uncertainties, suggesting a lack of evolution of the MBH-Mstellar relation up to z~2.5.

We present a detailed study of ionized outflows in a large sample of ~650 hard X-ray detected AGN. Using optical spectroscopy from the BAT AGN Spectroscopic Survey (BASS) we are able to reveal the faint wings of the [OIII] emission lines associated with outflows covering, for the first time, an unexplored range of low AGN bolometric luminosity at low redshift (z~0.05). We test if and how the incidence and velocity of ionized outflow is related to AGN physical parameters: black hole mass, gas column density, Eddington Ratio, [OIII], X-ray, and bolometric luminosities. We find a higher occurrence of ionized outflows in type 1.9 (55%) and type 1 AGN (46%) with respect to type 2 AGN (24%). While outflows in type 2 AGN are evenly balanced between blue and red velocity offsets with respect to the [OIII] narrow component, they are almost exclusively blueshifted in type 1 and type 1.9 AGN. We observe a significant dependence between the outflow occurrence and accretion rate, which becomes relevant at high Eddington ratios (> -1.7). We interpret such behaviour in the framework of covering factor-Eddington ratio dependence. We don't find strong trends of the outflow maximum velocity with AGN physical parameters, as an increase with bolometric luminosity can be only identified when including samples of AGN at high luminosity and high redshift taken from literature.

We present a numerical framework for the variability of active galactic nuclei (AGN), which links the variability of AGN over a broad range of timescales and luminosities to the observed properties of the AGN population as a whole, and particularly the Eddington ratio distribution function (ERDF). We have implemented our framework on GPU architecture, relying on previously published time series generating algorithms. After extensive tests that characterise several intrinsic and numerical aspects of the simulations, we describe some applications used for current and future time domain surveys and for the study of extremely variable sources (e.g., "changing look" or flaring AGN). Specifically, we define a simulation setup which reproduces the AGN variability observed in the PTF/iPTF survey, and use it to forward model longer light curves of the kind that may be observed within the LSST main survey. Thanks to our effcient implementations, these simulations are able to cover for example over 1 Myr with a roughly weekly cadence. We envision that this framework will become highly valuable to prepare for, and best exploit, data from upcoming time domain surveys, such as for example LSST.

We present a multi wavelength analysis of 28 of the most luminous low-redshift narrow-line, ultra-hard X-ray selected active galactic nuclei (AGN) drawn from the 70 month Swift/BAT all-sky survey, with bolometric luminosities of log(L_bol/erg/s) > 45.25. The broad goal of our study is to determine whether these objects have any distinctive properties, potentially setting them aside from lower-luminosity obscured AGN in the local Universe. Our analysis relies on the first data release of the BAT AGN Spectroscopic Survey (BASS/DR1) and on dedicated observations with the VLT, Palomar, and Keck observatories. We find that the vast majority of our sources agree with commonly used AGN selection criteria which are based on emission line ratios and on mid-infrared colours. Our AGN are predominantly hosted in massive galaxies (9.8 < log(M_*/M_sun) < 11.7); based on visual inspection of archival optical images, they appear to be mostly ellipticals. Otherwise, they do not have distinctive properties. Their radio luminosities, determined from publicly available survey data, show a large spread of almost 4 orders of magnitude - much broader than what is found for lower X-ray luminosity obscured AGN in BASS. Moreover, our sample shows no preferred combination of black hole masses (M_BH) and/or Eddington ratio (lambda_Edd), covering 7.5 < log(M_BH/M_sun) < 10.3 and 0.01 < lambda_Edd < 1. Based on the distribution of our sources in the lambda_Edd-N_H plane, we conclude that our sample is consistent with a scenario where the amount of obscuring material along the line of sight is determined by radiation pressure exerted by the AGN on the dusty circumnuclear gas.

We study the sudden optical and ultraviolet (UV) brightening of 1ES 1927+654, which until now was known as a narrow-line active galactic nucleus (AGN). 1ES 1927+654 was part of the small and peculiar class of "true Type-2" AGN, which lack broad emission lines and line-of-sight obscuration. Our high-cadence spectroscopic monitoring captures the appearance of a blue, featureless continuum, followed several weeks later by the appearance of broad Balmer emission lines. This timescale is generally consistent with the expected light travel time between the central engine and the broad-line emission region in (persistent) broad-line AGN. Hubble Space Telescope spectroscopy reveals no evidence for broad UV emission lines (e.g., CIV1549, CIII]1909, MgII2798), probably owing to dust in the broad-line emission region. To the best of our knowledge, this is the first case where the lag between the change in continuum and in broad-line emission of a "changing-look" AGN has been temporally resolved. The nature and timescales of the photometric and spectral evolution disfavor both a change in line-of-sight obscuration and a change of the overall rate of gas inflow as driving the drastic spectral transformations seen in this AGN. Although the peak luminosity and timescales are consistent with those of tidal disruption events seen in inactive galaxies, the spectral properties are not. The X-ray emission displays a markedly different behavior, with frequent flares on timescales of hours to days, and will be presented in a companion publication.

We have performed a very long baseline interferometry (VLBI) survey of local (z < 0.05) ultra hard X-ray (14-195 keV) selected active galactic nuclei (AGN) from the Swift Burst Alert Telescope (BAT) using KVN, KaVA, and VLBA. We first executed fringe surveys of 142 BAT-detected AGN at 15 or 22 GHz. Based on the fringe surveys and archival data, we find 10/279 nearby AGN (~4%) VLBI have 22 GHz flux above 30 mJy. This implies that the X-ray AGN with a bright nuclear jet are not common. Among these 10 radio-bright AGN, we obtained 22 GHz VLBI imaging data of our own for four targets and reprocessed archival data for six targets. We find that, although our 10 AGN observed with VLBI span a wide range of pc-scale morphological types, they lie on a tight linear relation between accretion luminosity and nuclear jet luminosity. Our result suggests that a powerful nuclear radio jet correlates with the accretion disc luminosity. We also probed the fundamental plane of black hole activity at VLBI scales (e.g., few milli-arcsecond). The jet luminosity and size distribution among our sample roughly fit into the proposed AGN evolutionary scenario, finding powerful jets after the blow-out phase based on the Eddington ratio (\lambda_{Edd})-hydrogen column density (N_{H}) relation. In addition, we find some hints of gas inflow or galaxy-galaxy merger in the majority of our sample. This implies that gas supply via tidal interactions in galactic scale may help the central AGN to launch a powerful parsec-scale jet.

The recently released 105-month {\it Swift}-Burst Alert Telescope (BAT) all-sky hard X-ray survey catalog presents an opportunity to study astrophysical objects detected in the deepest look at the entire hard X-ray (14$-$195 keV) sky. Here we report the results of a multifrequency study of 146 blazars from this catalog, quadrupling the number compared to past studies, by utilizing recent data from the {\it Fermi}-Large Area Telescope (LAT), Swift-BAT, and archival measurements. In our $\gamma$-ray analysis of $\sim$10 years of the LAT data, 101 are found as $\gamma$-ray emitters, whereas, 45 remains LAT undetected. We model the broadband spectral energy distributions with a synchrotron-inverse Compton radiative model. On average, BAT detected sources host massive black holes ($M_{\rm bh}\sim10^9$ M$_{\odot}$) and luminous accretion disks ($L_{\rm d}\sim10^{46}$ erg s$^{-1}$). At high-redshifts ($z>2$), BAT blazars host more powerful jets with luminous accretion disks compared to those detected only with the {\it Fermi}-LAT. We find good agreement in the black hole masses derived from the single-epoch optical spectroscopic measurements and standard accretion disk modeling approaches. Other physical properties of BAT blazars are similar to those known for {\it Fermi}-LAT detected objects.

LISA will open the mHz band of gravitational waves (GWs) to the astronomy community. The strong gravity which powers the variety of GW sources in this band is also crucial in a number of important astrophysical processes at the current frontiers of astronomy. These range from the beginning of structure formation in the early universe, through the origin and cosmic evolution of massive black holes in concert with their galactic environments, to the evolution of stellar remnant binaries in the Milky Way and in nearby galaxies. These processes and their associated populations also drive current and future observations across the electromagnetic (EM) spectrum. We review opportunities for science breakthroughs, involving either direct coincident EM+GW observations, or indirect multimessenger studies. We argue that for the US community to fully capitalize on the opportunities from the LISA mission, the US efforts should be accompanied by a coordinated and sustained program of multi-disciplinary science investment, following the GW data through to its impact on broad areas of astrophysics. Support for LISA-related multimessenger observers and theorists should be sized appropriately for a flagship observatory and may be coordinated through a dedicated mHz GW research center.

The ubiquitous variability of quasars across a wide range of wavelengths and timescales encodes critical information about the structure and dynamics of the circumnuclear emitting regions that are too small to be directly resolved, as well as detailed underlying physics of accretion and feedback processes in these active supermassive black holes. We emphasize the importance of studying quasar variability with time-domain spectroscopy, focusing on two science cases: (1) reverberation mapping (RM) to measure the broad-line region sizes and black hole masses in distant quasars; (2) spectroscopic follow-up of extreme variability quasars that dramatically change their continuum and broad-line flux within several years. We highlight the need for dedicated optical-infrared spectroscopic survey facilities in the coming decades to accompany wide-area time-domain imaging surveys, including: (1) the next phase of the Sloan Digital Sky Survey (SDSS-V; ~2020-2025), an all-sky, time-domain multi-object spectroscopic survey with 2.5m-class telescopes; (2) the planned Maunakea Spectroscopic Explorer, a dedicated 10m-class spectroscopic survey telescope with a 1.5 sq. deg field-of-view and multiplex of thousands of fibers in both optical and near-IR (J+H) to begin operations in 2029; (3) the Time-domain Spectroscopic Observatory (TSO), a proposed Probe-class ~1.3m telescope at L2, with imaging and spectroscopy (R=200, 1800) in 4 bands (0.3 - 5 micron) and rapid slew capability to 90% of sky, which will extend the coverage of Hbeta to z=8.

We report on the discovery and follow-up of a peculiar transient, OGLE17aaj, which occurred in the nucleus of a weakly active galaxy. We investigate whether it can be interpreted as a new candidate for a tidal disruption event (TDE). We present the OGLE-IV light curve that covers the slow 60-day-long rise to maximum along with photometric, spectroscopic, and X-ray follow-up during the first year. OGLE17aaj is a nuclear transient exhibiting some properties similar to previously found TDEs, including a long rise time, lack of colour-temperature evolution, and high black-body temperature. On the other hand, its narrow emission lines and slow post-peak evolution are different from previously observed TDEs. Its spectrum and light-curve evolution is similar to F01004-2237 and AT 2017bgt. Signatures of historical low-level nuclear variability suggest that OGLE17aaj may instead be related to a new type of accretion event in active super-massive black holes.

Accreting supermassive black holes (SMBHs) can exhibit variable emission across the electromagnetic spectrum and over a broad range of timescales. The variability of active galactic nuclei (AGNs) in the ultraviolet and optical is usually at the few tens of per cent level over timescales of hours to weeks. Recently, rare, more dramatic changes to the emission from accreting SMBHs have been observed, including tidal disruption events, 'changing look' AGNs and other extreme variability objects. The physics behind the 're-ignition', enhancement and 'shut-down' of accretion onto SMBHs is not entirely understood. Here we present a rapid increase in ultraviolet-optical emission in the centre of a nearby galaxy, marking the onset of sudden increased accretion onto a SMBH. The optical spectrum of this flare, dubbed AT 2017bgt, exhibits a mix of emission features. Some are typical of luminous, unobscured AGNs, but others are likely driven by Bowen fluorescence - robustly linked here with high-velocity gas in the vicinity of the accreting SMBH. The spectral features and increased ultraviolet flux show little evolution over a period of at least 14 months. This disfavours the tidal disruption of a star as their origin, and instead suggests a longer-term event of intensified accretion. Together with two other recently reported events with similar properties, we define a new class of SMBH-related flares. This has important implications for the classification of different types of enhanced accretion onto SMBHs.

Major galaxy mergers are thought to play an important part in fuelling the growth of supermassive black holes. However, observational support for this hypothesis is mixed, with some studies showing a correlation between merging galaxies and luminous quasars and others showing no such association. Recent observations have shown that a black hole is likely to become heavily obscured behind merger-driven gas and dust, even in the early stages of the merger, when the galaxies are well separated (5 to 40 kiloparsecs). Merger simulations further suggest that such obscuration and black-hole accretion peaks in the final merger stage, when the two galactic nuclei are closely separated (less than 3 kiloparsecs). Resolving this final stage requires a combination of high-spatial-resolution infrared imaging and high-sensitivity hard-X-ray observations to detect highly obscured sources. However, large numbers of obscured luminous accreting supermassive black holes have been recently detected nearby (distances below 250 megaparsecs) in X-ray observations. Here we report high-resolution infrared observations of hard-X-ray-selected black holes and the discovery of obscured nuclear mergers, the parent populations of supermassive-black-hole mergers. We find that obscured luminous black holes (bolometric luminosity higher than 2x10^44 ergs per second) show a significant (P<0.001) excess of late-stage nuclear mergers (17.6 per cent) compared to a sample of inactive galaxies with matching stellar masses and star formation rates (1.1 per cent), in agreement with theoretical predictions. Using hydrodynamic simulations, we confirm that the excess of nuclear mergers is indeed strongest for gas-rich major-merger hosts of obscured luminous black holes in this final stage.

We quantify the luminosity contribution of active galactic nuclei (AGN) to the 12 $\mu$m, mid-infrared (MIR; 5-38 $\mu$m), and the total IR (5-1000 $\mu$m) emission in the local AGN detected in the all-sky 70-month Swift/Burst Alert Telescope (BAT) ultra hard X-ray survey. We decompose the IR spectral energy distributions (SEDs) of 587 objects into AGN and starburst components using AGN torus and star-forming galaxy templates. This enables us to recover the AGN torus emission also for low-luminosity end, down to $\log (L_{14-150}/{\rm erg}~{\rm s}^{-1}) \simeq 41$, which typically have significant host galaxy contamination. We find that the luminosity contribution of the AGN to the 12 $\mu$m, the MIR, and the total IR band is an increasing function of the 14-150 keV luminosity. We also find that for the most extreme cases, the IR pure-AGN emission from the torus can extend up to 90 $\mu$m. The obtained total IR AGN luminosity through the IR SED decomposition enables us to estimate the fraction of the sky obscured by dust, i.e., the dust covering factor. We demonstrate that the median of the dust covering factor is always smaller than that of the X-ray obscuration fraction above the AGN bolometric luminosity of $\log (L_{\rm bol}/{\rm erg}~{\rm s}^{-1}) \simeq 42.5$. Considering that X-ray obscuration fraction is equivalent to the covering factor coming from both the dust and gas, it indicates that an additional neutral gas component, along with the dusty torus, is responsible for the absorption of X-ray emission.

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.

In order to relate the observed evolution of the galaxy stellar mass function and the luminosity function of active galactic nuclei (AGN), we explore a co-evolution scenario in which AGN are associated only with the very last phases of the star-forming life of a galaxy. We derive analytically the connections between the parameters of the observed quasar luminosity functions and galaxy mass functions. The $(m_{\rm bh}/m_{*})_{Qing}$ associated with quenching is given by the ratio of the global black hole accretion rate density (BHARD) and star-formation rate density (SFRD) at the epoch in question. Observational data on the SFRD and BHARD suggests $(m_{\rm bh}/m_{*})_{Qing} \propto (1+z)^{1.5}$ below redshift 2. This evolution reproduces the observed mass-luminosity plane of SDSS quasars, and also reproduces the local $m_{\rm bh}/m_{*}$ relation in passive galaxies. The characteristic Eddington ratio, $\lambda^*$, is derived from both the BHARD/SFRD ratio and the evolving $L^*$ of the AGN population. This increases up to $z \sim 2$ as $\lambda^* \propto (1+z)^{2.5}$ but at higher redshifts, $\lambda^*$ stabilizes at the physically interesting Eddington limit, $\lambda^* \sim 1$. The new model may be thought of as an opposite extreme to our earlier co-evolution scenario in Caplar et al. 2015. The main observable difference between the two co-evolution scenarios, presented here and in Caplar et al. 2015, is in the active fraction of low mass star-forming galaxies. We compare the predictions with the data from deep multi-wavelength surveys and find that the "quenching" scenario developed in the current paper is much to be preferred.

The bulk of the X-ray emission in Active Galactic Nuclei (AGN) is produced very close to the accreting supermassive black hole (SMBH), in a corona of hot electrons which up scatters optical and ultraviolet photons from the accretion flow. The cutoff energy ($E_{\rm C}$) of the primary X-ray continuum emission carries important information on the physical characteristics of the X-ray emitting plasma, but little is currently known about its potential relation with the properties of accreting SMBHs. Using the largest broad-band (0.3-150 keV) X-ray spectroscopic study available to date, we investigate how the corona is related to the AGN luminosity, black hole mass and Eddington ratio ($\lambda_{\rm Edd}$). Assuming a slab corona the median values of the temperature and optical depth of the Comptonizing plasma are $kT_{\rm e}=105 \pm 18$ keV and $\tau=0.25\pm0.06$, respectively. When we properly account for the large number of $E_{\rm C}$ lower limits, we find a statistically significant dependence of the cutoff energy on the Eddington ratio. In particular, objects with $ \lambda_{\rm Edd}>0.1$ have a significantly lower median cutoff energy ($E_{\rm C}=160\pm41$ keV) than those with $\lambda_{\rm Edd}\leq 0.1$ ($E_{\rm C}=370\pm51$ keV). This is consistent with the idea that radiatively compact coronae are also cooler, because they tend to avoid the region in the temperature-compactness parameter space where runaway pair production would dominate. We show that this behaviour could also straightforwardly explain the suggested positive correlation between the photon index ($\Gamma$) and the Eddington ratio, being able to reproduce the observed slope of the $\Gamma-\lambda_{\rm Edd}$ trend.

In large optical survey at high redshifts ($z>2$), the CIV broad emission line is the most practical alternative to estimate the mass ($M_{\text{BH}}$) of active super-massive black holes (SMBHs). However, mass determinations obtained with this line are known to be highly uncertain. In this work we use the Sloan Digital Sky Survey Data Release 7 and 12 quasar catalogues to statistically test three alternative methods put forward in the literature to improve CIV-based $M_{\text{BH}}$ estimations. These methods are constructed from correlations between the ratio of the CIV line-width to the low ionization line-widths (H$\alpha$, H$\beta$ and MgII) and several other properties of rest-frame UV emission lines. Our analysis suggests that these correction methods are of limited applicability, mostly because all of them depend on correlations that are driven by the linewidth of the CIV profile itself and not by an interconnection between the linewidth of the CIV line with the linewidth of the low ionization lines. Our results show that optical CIV-based mass estimates at high redshift cannot be a proper replacement for estimates based on IR spectroscopy of low ionization lines like H$\alpha$, H$\beta$ and MgII

We characterize the environments of local accreting supermassive black holes by measuring the clustering of AGN in the Swift/BAT Spectroscopic Survey (BASS). With 548 AGN in the redshift range 0.01<z<0.1 over the full sky from the DR1 catalog, BASS provides the largest, least biased sample of local AGN to date due to its hard X-ray selection (14-195 keV) and rich multiwavelength/ancillary data. By measuring the projected cross-correlation function between the AGN and 2MASS galaxies, and interpreting it via halo occupation distribution (HOD) and subhalo-based models, we constrain the occupation statistics of the full sample, as well as in bins of absorbing column density and black hole mass. We find that AGN tend to reside in galaxy group environments, in agreement with previous studies of AGN throughout a large range of luminosity and redshift, and that on average they occupy their dark matter halos similar to inactive galaxies of comparable stellar mass. We also find evidence that obscured AGN tend to reside in denser environments than unobscured AGN, even when samples were matched in luminosity, redshift, stellar mass, and Eddington ratio. We show that this can be explained either by significantly different halo occupation distributions or statistically different host halo assembly histories. Lastly, we see that massive black holes are slightly more likely to reside in central galaxies than black holes of smaller mass.

We present a framework to link and describe AGN variability on a wide range of timescales, from days to billions of years. In particular, we concentrate on the AGN variability features related to changes in black hole fuelling and accretion rate. In our framework, the variability features observed in different AGN at different timescales may be explained as realisations of the same underlying statistical properties. In this context, we propose a model to simulate the evolution of AGN light curves with time based on the probability density function (PDF) and power spectral density (PSD) of the Eddington ratio ($L/L_{\rm Edd}$) distribution. Motivated by general galaxy population properties, we propose that the PDF may be inspired by the $L/L_{\rm Edd}$ distribution function (ERDF), and that a single (or limited number of) ERDF+PSD set may explain all observed variability features. After outlining the framework and the model, we compile a set of variability measurements in terms of structure function (SF) and magnitude difference. We then combine the variability measurements on a SF plot ranging from days to Gyr. The proposed framework enables constraints on the underlying PSD and the ability to link AGN variability on different timescales, therefore providing new insights into AGN variability and black hole growth phenomena.

We use a phenomenological model which connects the galaxy and AGN populations to investigate the process of AGN triggering through major galaxy mergers at z~0. The model uses stellar mass functions as input and allows the prediction of AGN luminosity functions based on assumed Eddington ratio distribution functions (ERDFs). We show that the number of AGN hosted by merger galaxies relative to the total number of AGN increases as a function of AGN luminosity. This is due to more massive galaxies being more likely to undergo a merger and does not require the assumption that mergers lead to higher Eddington ratios than secular processes. Our qualitative analysis also shows that to match the observations, the probability of a merger galaxy hosting an AGN and accreting at a given Eddington value has to be increased by a factor ~10 relative to the general AGN population. An additional significant increase of the fraction of high Eddington ratio AGN among merger host galaxies leads to inconsistency with the observed X-ray luminosity function. Physically our results imply that, compared to the general galaxy population, the AGN fraction among merger galaxies is ~10 times higher. On average, merger triggering does however not lead to significantly higher Eddington ratios.

We present a catalog of hard X-ray sources detected in the first 105 months of observations with the Burst Alert Telescope (BAT) coded mask imager on board the Swift observatory. The 105 month Swift-BAT survey is a uniform hard X-ray all-sky survey with a sensitivity of $8.40\times 10^{-12}\ {\rm erg\ s^{-1}\ cm^{-2}}$ over 90% of the sky and $7.24\times 10^{-12}\ {\rm erg\ s^{-1}\ cm^{-2}}$ over 50% of the sky in the 14-195 keV band. The Swift-BAT 105 month catalog provides 1632 (422 new detections) hard X-ray sources in the 14-195 keV band above the 4.8{\sigma} significance level. Adding to the previously known hard X-ray sources, 34% (144/422) of the new detections are identified as Seyfert AGN in nearby galaxies (z<0.2). The majority of the remaining identified sources are X-ray binaries (7%, 31) and blazars/BL Lac objects (10%, 43). As part of this new edition of the Swift-BAT catalog, we release eight-channel spectra and monthly sampled light curves for each object in the online journal and at the Swift-BAT 105 month Web site.

We present a long-term, multi-wavelength project to understand the epoch of fastest growth of the most massive black holes by using a sample of 40 luminous quasars at z~4.8. These quasars have rather uniform properties, with typical accretion rates and black hole masses of L/L_Edd~0.7 and M_BH~10^9 M_sol. The sample consists of "FIR-bright" sources with a previous Herschel/SPIRE detection, suggesting SFR > 1000 M_sol/yr, as well as of "FIR-faint" sources for which Herschel stacking analysis implies a typical SFR of ~400 M_sol/yr. Six of the quasars have been observed by ALMA in [CII]{\lambda}157.74 micron line emission and adjacent rest-frame 150 micron continuum, to study the dusty cold ISM. ALMA detected companion, spectroscopically confirmed sub-mm galaxies (SMGs) for three sources - one FIR-bright and two FIR-faint. The companions are separated by ~14-45 kpc from the quasar hosts, and we interpret them as major galaxy interactions. Our ALMA data therefore clearly support the idea that major mergers may be important drivers for rapid, early SMBH growth. However, the fact that not all high-SFR quasar hosts are accompanied by interacting SMGs, and their ordered gas kinematics observed by ALMA, suggest that other processes may be fueling these systems. Our analysis thus demonstrates the diversity of host galaxy properties and gas accretion mechanisms associated with early and rapid SMBH growth.

Hard X-ray ($\geq 10$ keV) observations of Active Galactic Nuclei (AGN) can shed light on some of the most obscured episodes of accretion onto supermassive black holes. The 70-month Swift/BAT all-sky survey, which probes the 14-195 keV energy range, has currently detected 838 AGN. We report here on the broad-band X-ray (0.3-150 keV) characteristics of these AGN, obtained by combining XMM-Newton, Swift/XRT, ASCA, Chandra, and Suzaku observations in the soft X-ray band ($\leq 10$ keV) with 70-month averaged Swift/BAT data. The non-blazar AGN of our sample are almost equally divided into unobscured ($N_{\rm H}< 10^{22}\rm cm^{-2}$) and obscured ($N_{\rm H}\geq 10^{22}\rm cm^{-2}$) AGN, and their Swift/BAT continuum is systematically steeper than the 0.3-10 keV emission, which suggests that the presence of a high-energy cutoff is almost ubiquitous. We discuss the main X-ray spectral parameters obtained, such as the photon index, the reflection parameter, the energy of the cutoff, neutral and ionized absorbers, and the soft excess for both obscured and unobscured AGN.

Supermassive black holes reside in the nuclei of most galaxies. Accurately determining their mass is key to understand how the population evolves over time and how the black holes relate to their host galaxies. Beyond the local universe, the mass is commonly estimated assuming virialized motion of gas in the close vicinity to the active black holes, traced through broad emission lines. However, this procedure has uncertainties associated with the unknown distribution of the gas clouds. Here we show that the comparison of black hole masses derived from the properties of the central accretion disc with the virial mass estimate provides a correcting factor, for the virial mass estimations, that is inversely proportional to the observed width of the broad emission lines. Our results suggest that line-of-sight inclination of gas in a planar distribution can account for this effect. However, radiation pressure effects on the distribution of gas can also reproduce our findings. Regardless of the physical origin, our findings contribute to mitigate the uncertainties in current black hole mass estimations and, in turn, will help to further understand the evolution of distant supermassive black holes and their host galaxies.

We explore the relationship between X-ray absorption and optical obscuration within the BAT AGN Spectroscopic Survey (BASS) which has been collecting and analyzing the optical and X-ray spectra for 641 hard X-ray selected ($E>14$ keV) active galactic nuclei (AGN). We use the deviation from a linear broad H$\alpha$-to-X-ray relationship as an estimate of the maximum optical obscuration towards the broad line region and compare the $A_{\rm V}$ to the hydrogen column densities ($N_{\rm H}$) found through systematic modeling of their X-ray spectra. We find that the inferred columns implied by $A_{\rm V}$ towards the broad line region (BLR) are often orders of magnitude less than the columns measured towards the X-ray emitting region indicating a small scale origin for the X-ray absorbing gas. After removing 30\% of Sy 1.9s that potentially have been misclassified due to outflows, we find that 86\% (164/190) of the Type 1 population (Sy 1--1.9) are X-ray unabsorbed as expected based on a single obscuring structure. However, 14\% (26/190), of which 70\% (18/26) are classified as Sy 1.9, are X-ray absorbed, suggesting the broad line region itself is providing extra obscuration towards the X-ray corona. The fraction of X-ray absorbed Type 1 AGN remains relatively constant with AGN luminosity and Eddington ratio, indicating a stable broad line region covering fraction.

The large majority of the accreting supermassive black holes in the Universe are obscured by large columns of gas and dust. The location and evolution of this obscuring material have been the subject of intense research in the past decades, and are still highly debated. A decrease in the covering factor of the circumnuclear material with increasing accretion rates has been found by studies carried out across the electromagnetic spectrum. The origin of this trend has been suggested to be driven either by the increase in the inner radius of the obscuring material with incident luminosity due to the sublimation of dust; by the gravitational potential of the black hole; by radiative feedback; or by the interplay between outflows and inflows. However, the lack of a large, unbiased and complete sample of accreting black holes, with reliable information on gas column density, luminosity and mass, has left the main physical mechanism regulating obscuration unclear. Using a systematic multi-wavelength survey of hard X-ray-selected black holes, here we show that radiation pressure on dusty gas is indeed the main physical mechanism regulating the distribution of the circumnuclear material. Our results imply that the bulk of the obscuring dust and gas in these objects is located within the sphere of influence of the black hole (i.e., a few to tens of parsecs), and that it can be swept away even at low radiative output rates. The main physical driver of the differences between obscured and unobscured accreting black holes is therefore their mass-normalized accretion rate.

We present results of a ground-based near-infrared campaign with Palomar TripleSpec, Keck NIRSPEC, and Gemini GNIRS to target two samples of reddened active galactic nucleus (AGN) candidates from the 31 deg$^2$ Stripe 82 X-ray survey. One sample, which is $\sim$89\% complete to $K<16$ (Vega), consists of eight confirmed AGNs, four of which were identified with our follow-up program, and is selected to have red $R-K$ colors ($>4$, Vega). The fainter sample ($K>17$, Vega) represents a pilot program to follow-up four sources from a parent sample of 34 that are not detected in the single-epoch SDSS catalog and have {\it WISE} quasar colors. All twelve sources are broad-line AGNs (at least one permitted emission line has a FWHM exceeding 1300 km s$^{-1}$) and span a redshift range $0.59 < z < 2.5$. Half the ($R-K$)-selected AGNs have features in their spectra suggestive of outflows. When comparing these sources to a matched sample of blue Type 1 AGNs, we find the reddened AGNs are more distant ($z > 0.5$) and a greater percentage have high X-ray luminosities ($L_{\rm X,full} > 10^{44}$ erg s$^{-1}$). Such outflows and high luminosities may be consistent with the paradigm that reddened broad-line AGNs represent a transitory phase in AGN evolution as described by the major merger model for black hole growth. Results from our pilot program demonstrate proof-of-concept that our selection technique is successful in discovering reddened quasars at $z > 1$ missed by optical surveys.

We present the first catalog and data release of the Swift-BAT AGN Spectroscopic Survey (BASS). We analyze optical spectra of the majority of AGN (77%, 641/836) detected based on their 14-195 keV emission in the 70-month Swift BAT all-sky catalog. This includes redshift determination, absorption and emission line measurements, and black hole mass and accretion rate estimates for the majority of obscured and un-obscured AGN (74%, 473/641) with 340 measured for the first time. With ~90% of sources at z<0.2, the survey represents a significant census of hard-X-ray selected AGN in the local universe. In this first catalog paper, we describe the spectroscopic observations and datasets, and our initial spectral analysis. The FWHM of the emission lines show broad agreement with the X-ray obscuration (~94%), such that Sy 1-1.8 have NH<10^21.9 cm^-2, and Seyfert 2, have NH>10^21.9 cm^-2. Seyfert 1.9 show a range of column densities. Compared to narrow line AGN in the SDSS, the X-ray selected AGN have a larger fraction of dusty host galaxies suggesting these types of AGN are missed in optical surveys. Using the most sensitive [OIII]/Hbeta and [NII]/Halpha emission line diagnostic, about half of the sources are classified as Seyferts, ~15% reside in dusty galaxies that lack an Hbeta detection, but for which the line upper limits imply either a Seyfert or LINER, ~15% are in galaxies with weak or no emission lines despite high quality spectra, and a few percent each are LINERS, composite galaxies, HII regions, or in known beamed AGN.

We use a phenomenological model to show that black hole growth in the local Universe (z < 0.1) can be described by two separate, mass independent Eddington ratio distribution functions (ERDFs). We assume that black holes can be divided into two independent groups: those with radiatively efficient accretion, primarily hosted by optically blue and green galaxies, and those with radiatively inefficient accretion, which are mainly found in red galaxies. With observed galaxy stellar mass functions as input, we show that the observed AGN luminosity functions can be reproduced by using mass independent, broken power law shaped ERDFs. We use the observed hard X-ray and 1.4 GHz radio luminosity functions to constrain the ERDF for radiatively efficient and inefficient AGN, respectively. We also test alternative ERDF shapes and mass dependent models. Our results are consistent with a mass independent AGN fraction and AGN hosts being randomly drawn from the galaxy population. We argue that the ERDF is not shaped by galaxy-scale effects, but by how efficiently material can be transported from the inner few parsecs to the accretion disc. Our results are incompatible with the simplest form of mass quenching where massive galaxies host higher accretion rate AGN. Furthermore, if reaching a certain Eddington ratio is a sufficient condition for maintenance mode, it can occur in all red galaxies, not just the most massive ones.

To fully understand cosmic black hole growth we need to constrain the population of heavily obscured active galactic nuclei (AGN) at the peak of cosmic black hole growth ($z\sim$1-3). Sources with obscuring column densities higher than $\mathrm{10^{24}}$ atoms $\mathrm{cm^{-2}}$, called Compton-thick (CT) AGN, can be identified by excess X-ray emission at $\sim$20-30 keV, called the "Compton hump". We apply the recently developed Spectral Curvature (SC) method to high-redshift AGN (2<z<5) detected with Chandra. This method parametrizes the characteristic "Compton hump" feature cosmologically redshifted into the X-ray band at observed energies <10 keV. We find good agreement in CT AGN found using the SC method and bright sources fit using their full spectrum with X-ray spectroscopy. In the Chandra deep field south, we measure a CT fraction of $\mathrm{17^{+19}_{-11}\%}$ (3/17) for sources with observed luminosity $\mathrm{>5\times 10^{43}}$ erg $\mathrm{s^{-1}}$. In the Cosmological evolution survey (COSMOS), we find an observed CT fraction of $\mathrm{15^{+4}_{-3}\%}$ (40/272) or $\mathrm{32\pm11 \%}$ when corrected for the survey sensitivity. When comparing to low redshift AGN with similar X-ray luminosities, our results imply the CT AGN fraction is consistent with having no redshift evolution. Finally, we provide SC equations that can be used to find high-redshift CT AGN (z>1) for current (XMM-Newton) and future (eROSITA and ATHENA) X-ray missions.

We investigate the star formation properties of a large sample of ~2300 X-ray-selected Type 2 Active Galactic Nuclei (AGNs) host galaxies out to z~3 in the Chandra COSMOS Legacy Survey in order to understand the connection between the star formation and nuclear activity. Making use of the existing multi-wavelength photometric data available in the COSMOS field, we perform a multi-component modeling from far-infrared to near-ultraviolet using a nuclear dust torus model, a stellar population model and a starburst model of the spectral energy distributions (SEDs). Through detailed analysis of SEDs, we derive the stellar masses and the star formation rates (SFRs) of Type 2 AGN host galaxies. The stellar mass of our sample is in the range 9 < log M_{stellar}/M_{\odot} < 12 with uncertainties of ~0.19 dex. We find that Type 2 AGN host galaxies have, on average, similar SFRs compared to the normal star-forming galaxies with similar M_{stellar} and redshift ranges, suggesting no significant evidence for enhancement or quenching of star formation. This could be interpreted in a scenario, where the relative massive galaxies have already experienced substantial growth at higher redshift (z>3), and grow slowly through secular fueling processes hosting moderate-luminosity AGNs.

We study the observed relation between accretion rate (in terms of L/L_Edd) and shape of the hard X-ray spectral energy distribution (namely the photon index Gamma_X) for a large sample of 228 hard X-ray selected, low-redshift active galactic nuclei (AGN), drawn from the Swift/BAT AGN Spectroscopic Survey (BASS). This includes 30 AGN for which black hole mass (and therefore L/L_Edd) is measured directly through masers, spatially resolved gas or stellar dynamics, or reverberation mapping. The high quality and broad energy coverage of the data provided through BASS allow us to examine several alternative determinations of both Gamma_X and L/L_Edd. For the BASS sample as a whole, we find a statistically significant, albeit very weak correlation between Gamma_X and L/L_Edd. The best-fitting relations we find, Gamma_X=0.15 log(L/L_Edd)+const., are considerably shallower than those reported in previous studies. Moreover, we find no corresponding correlations among the subsets of AGN with different M_BH determination methodology. In particular, we find no robust evidence for a correlation when considering only those AGN with direct or single-epoch M_BH estimates. This latter finding is in contrast to several previous studies which focused on z>0.5 broad-line AGN. We discuss this tension and conclude that it can be partially accounted for if one adopts a simplified, power-law X-ray spectral model, combined with L/L_Edd estimates that are based on the continuum emission and on single-epoch broad line spectroscopy in the optical regime. We finally highlight the limitations on using Gamma_X as a probe of supermassive black hole evolution in deep extragalactic X-ray surveys.

We present new ALMA band-7 data for a sample of six luminous quasars at z~4.8, powered by fast-growing supermassive black holes (SMBHs) with rather uniform properties: the typical accretion rates and black hole masses are L/L_Edd~0.7 and M_BH~10^9 M_sol. Our sample consists of three "FIR-bright" sources, which were individually detected in previous Herschel/SPIRE observations, with star formation rates of SFR>1000 M_sol/yr, and three "FIR-faint" sources for which Herschel stacking analysis implies a typical SFR of ~400 M_sol/yr. The dusty interstellar medium in the hosts of all six quasars is clearly detected in the ALMA data and resolved on scales of ~2 kpc, in both continuum (\lambda_rest ~150 um) and [C II]157.74 um line emission. The continuum emission is in good agreement with the expectations from the Herschel data, confirming the intense SF activity in the quasar hosts. Importantly, we detect companion sub-millimeter galaxies (SMGs) for three sources -- one FIR-bright and two FIR-faint, separated by ~14-45 kpc and <450 km/s from the quasar hosts. The [C II]-based dynamical mass estimates for the interacting SMGs are within a factor of ~3 of the quasar hosts' masses, while the continuum emission implies SFR(quasar)~(2-11)xSFR(SMG). Our ALMA data therefore clearly support the idea that major mergers are important drivers for rapid early SMBH growth. However, the fact that not all high-SFR quasar hosts are accompanied by interacting SMGs and the gas kinematics as observed by ALMA suggest that other processes may be fueling these systems. Our analysis thus demonstrates the diversity of host galaxy properties and gas accretion mechanisms associated with early and rapid SMBH growth.

We estimate the accretion rates onto the supermassive black holes that power 20 of the highest-redshift quasars, at z>5.8, including the quasar with the highest redshift known to date -- ULAS J1120 at z=7.09. The analysis is based on the observed (rest-frame) optical luminosities and reliable "virial" estimates of the BH masses of the quasars, and utilizes scaling relations derived from thin accretion disk theory. The mass accretion rates through the postulated disks cover a wide range, dM_disk/dt~4-190 Msol/yr, with most of the objects (80\%) having dM_disk/dt~10-65 Msol/yr, confirming the Eddington-limited nature of the accretion flows. By combining our estimates of dM_disk/dt with conservative, lower limits on the bolometric luminosities of the quasars, we investigate which alternative values of \eta best account for all the available data. We find that the vast majority of quasars (~85\%) can be explained with radiative efficiencies in the range \eta~0.03-0.3, with a median value close to the commonly assumed \eta=0.1. Within this range, we obtain conservative estimates of \eta>0.14 for ULAS J1120 and SDSS J0100 (at z=6.3), and of >0.19 for SDSS J1148 (at z=6.41; assuming their BH masses are accurate). The implied accretion timescales are generally in the range t_acc = M_BH / dM_BH/dt ~0.1-1 Gyr, suggesting that most quasars could have had ~1-10 mass e-foldings since BH seed formation. Our analysis therefore demonstrates that the available luminosities and masses for the highest-redshift quasars can be explained self-consistently within the thin, radiatively efficient accretion disk paradigm. Episodes of radiatively inefficient, "super-critical" accretion may have occurred at significantly earlier epochs (i.e., z~10 or above).

We provide a comprehensive census of the near-Infrared (NIR, 0.8-2.4 $\mu$m) spectroscopic properties of 102 nearby (z < 0.075) active galactic nuclei (AGN), selected in the hard X-ray band (14-195 keV) from the Swift-Burst Alert Telescope (BAT) survey. With the launch of the James Webb Space Telescope this regime is of increasing importance for dusty and obscured AGN surveys. We measure black hole masses in 68% (69/102) of the sample using broad emission lines (34/102) and/or the velocity dispersion of the Ca II triplet or the CO band-heads (46/102). We find that emission line diagnostics in the NIR are ineffective at identifying bright, nearby AGN galaxies because ([Fe II] 1.257$\mu$m/Pa$\beta$ and H$_2$ 2.12$\mu$m/Br$\gamma$) identify only 25% (25/102) as AGN with significant overlap with star forming galaxies and only 20% of Seyfert 2 have detected coronal lines (6/30). We measure the coronal line emission in Seyfert 2 to be weaker than in Seyfert 1 of the same bolometric luminosity suggesting obscuration by the nuclear torus. We find that the correlation between the hard X-ray and the [Si VI] coronal line luminosity is significantly better than with the [O III] luminosity. Finally, we find 3/29 galaxies (10%) that are optically classified as Seyfert 2 show broad emission lines in the NIR. These AGN have the lowest levels of obscuration among the Seyfert 2s in our sample ($\log N_{\rm H} < 22.43$ cm$^{-2}$), and all show signs of galaxy-scale interactions or mergers suggesting that the optical broad emission lines are obscured by host galaxy dust.

We systematically investigate the near- (NIR) to far-infrared (FIR) photometric properties of a nearly complete sample of local active galactic nuclei (AGN) detected in the Swift/Burst Alert Telescope (BAT) all-sky ultra hard X-ray (14-195 keV) survey. Out of 606 non-blazar AGN in the Swift/BAT 70-month catalog at high galactic latitude of $|b|>10^{\circ}$, we obtain IR photometric data of 604 objects by cross-matching the AGN positions with catalogs from the WISE, AKARI, IRAS, and Herschel infrared observatories. We find a good correlation between the ultra-hard X-ray and mid-IR (MIR) luminosities over five orders of magnitude ($41 < \log (L_{14-195}/{\rm erg}~{\rm s}^{-1})< 46$). Informed by previous measures of the intrinsic spectral energy distribution of AGN, we find FIR pure-AGN candidates whose FIR emission is thought to be AGN-dominated with low starformation activity. We demonstrate that the dust covering factor decreases with the bolometric AGN luminosity, confirming the luminosity-dependent unified scheme. We also show that the completeness of the WISE color-color cut in selecting Swift/BAT AGN increases strongly with 14-195 keV luminosity.

We characterize the optical variability of quasars in the intermediate Palomar Transient Factory (iPTF) and Palomar Transient Factory (PTF) surveys. We re-calibrate the $r$-band light curves for $\sim$28,000 luminous, broad-line AGNs from the SDSS, producing a total of $\sim$2.4 million photometric data points. We utilize both the structure function (SF) and power spectrum density (PSD) formalisms to search for links between the optical variability and the physical parameters of the accreting supermassive black holes that power the quasars. The excess variance (SF$^{2}$) of the quasar sample tends to zero at very short time separations, validating our re-calibration of the time-series data. We find that the the amplitude of variability at a given time-interval, or equivalently the time-scale of variability to reach a certain amplitude, is most strongly correlated with luminosity with weak or no dependence on black hole mass and redshift. For a variability level of SF($\tau$)=0.07 mag, the time-scale has a dependency of $\tau \propto L^{0.4}$. This is broadly consistent with the expectation from a simple Keplerian accretion disk model, which provides $\tau \propto L^{0.5}$. The PSD analysis also reveals that many quasar light curves are steeper than a damped random walk. We find a correlation between the steepness of the PSD slopes, specifically the fraction of slopes steeper than 2.5, and black hole mass, although we cannot exclude the possibility that luminosity or Eddington ratio are the drivers of this effect. This effect is also seen in the SF analysis of the (i)PTF data, and in a PSD analysis of quasars in the SDSS Stripe 82.

We investigate the observed relationship between black hole mass ($M_{\rm BH}$), bolometric luminosity ($L_{\rm bol}$), and Eddington ratio (${\lambda}_{\rm Edd}$) with optical emission line ratios ([NII] {\lambda}6583/H{\alpha}, [SII] {\lambda}{\lambda}6716,6731/H{\alpha}, [OI] {\lambda}6300/H{\alpha}, [OIII] {\lambda}5007/H{\beta}, [NeIII] {\lambda}3869/H{\beta}, and HeII {\lambda}4686/H{\beta}) of hard X-ray-selected AGN from the BAT AGN Spectroscopic Survey (BASS). We show that the [NII] {\lambda}6583/H{\alpha} ratio exhibits a significant correlation with ${\lambda}_{\rm Edd}$ ($R_{\rm Pear}$ = -0.44, $p$-value=$3\times10^{-13}$, {\sigma} = 0.28 dex), and the correlation is not solely driven by $M_{\rm BH}$ or $L_{\rm bol}$. The observed correlation between [NII] {\lambda}6583/H{\alpha} ratio and $M_{\rm BH}$ is stronger than the correlation with $L_{\rm bol}$, but both are weaker than the ${\lambda}_{\rm Edd}$ correlation. This implies that the large-scale narrow lines of AGN host galaxies carry information about the accretion state of the AGN central engine. We propose that the [NII] {\lambda}6583/H{\alpha} is a useful indicator of Eddington ratio with 0.6 dex of rms scatter, and that it can be used to measure ${\lambda}_{\rm Edd}$ and thus $M_{\rm BH}$ from the measured $L_{\rm bol}$, even for high redshift obscured AGN. We briefly discuss possible physical mechanisms behind this correlation, such as the mass-metallicity relation, X-ray heating, and radiatively driven outflows.

We present the measurement of the projected and redshift space 2-point correlation function (2pcf) of the new catalog of Chandra COSMOS-Legacy AGN at 2.9$\leq$z$\leq$5.5 ($\langle L_{bol} \rangle \sim$10$^{46}$ erg/s) using the generalized clustering estimator based on phot-z probability distribution functions (Pdfs) in addition to any available spec-z. We model the projected 2pcf estimated using $\pi_{max}$ = 200 h$^{-1}$ Mpc with the 2-halo term and we derive a bias at z$\sim$3.4 equal to b = 6.6$^{+0.60}_{-0.55}$, which corresponds to a typical mass of the hosting halos of log M$_h$ = 12.83$^{+0.12}_{-0.11}$ h$^{-1}$ M$_{\odot}$. A similar bias is derived using the redshift-space 2pcf, modelled including the typical phot-z error $\sigma_z$ = 0.052 of our sample at z$\geq$2.9. Once we integrate the projected 2pcf up to $\pi_{max}$ = 200 h$^{-1}$ Mpc, the bias of XMM and \textit{Chandra} COSMOS at z=2.8 used in Allevato et al. (2014) is consistent with our results at higher redshift. The results suggest only a slight increase of the bias factor of COSMOS AGN at z$\gtrsim$3 with the typical hosting halo mass of moderate luminosity AGN almost constant with redshift and equal to logM$_h$ = 12.92$^{+0.13}_{-0.18}$ at z=2.8 and log M$_h$ = 12.83$^{+0.12}_{-0.11}$ at z$\sim$3.4, respectively. The observed redshift evolution of the bias of COSMOS AGN implies that moderate luminosity AGN still inhabit group-sized halos at z$\gtrsim$3, but slightly less massive than observed in different independent studies using X-ray AGN at z$\leq2$.

We investigate the sensitivity of the colour-based quasar selection algorithm of the Sloan Digital Sky Survey to several key physical parameters of supermassive black holes (SMBHs), focusing on BH spin ($a_{\star}$) at the high BH-mass regime ($M_{BH} \geqslant10^9\, M_{\odot}$). We use a large grid of model spectral energy distribution, assuming geometrically-thin, optically-thick accretion discs, and spanning a wide range of five physical parameters: BH mass $M_{BH}$, BH spin $a_{\star}$, Eddington ratio $L / L_{Edd}$ , redshift $z$, and inclination angle $inc$. Based on the expected fluxes in the SDSS imaging ugriz bands, we find that $\sim 99.8\%$ of our models with $M_{BH} \leqslant 10^{9.5}\, M_{\odot}$ are selected as quasar candidates and thus would have been targeted for spectroscopic follow-up. However, in the extremely high-mass regime, $\geqslant 10^{10} M_{\odot}$, we identify a bias against slowly/retrograde spinning SMBHs. The fraction of SEDs that would have been selected as quasar candidates drops below $\sim50\%$ for $a_{\star} <0$ across $0.5<z<2$. For particularly massive BHs, with $M_{BH} \simeq 3\times10^{10}\, M_{\odot}$, this rate drops below $\sim20\%$, and can be yet lower for specific redshifts. We further find that the chances of identifying any hypothetical sources with $M_{BH} = 10^{11}\, M_{\odot}$ by colour selection would be extremely low at the level of $\sim 3\%$. Our findings, along with several recent theoretical arguments and empirical findings, demonstrate that the current understanding of the SMBH population at the high-$M_{BH}$, and particularly the low- or retrograde-spinning regime, is highly incomplete.

We present the largest high-redshift (3<z<6.85) sample of X-ray-selected active galactic nuclei (AGN) on a contiguous field, using sources detected in the Chandra COSMOS Legacy survey. The sample contains 174 sources, 87 with spectroscopic redshift, the other 87 with photometric redshift (z_phot). In this work we treat z_phot as a probability weighted sum of contributions, adding to our sample the contribution of sources with z_phot<3 but z_phot probability distribution >0 at z>3. We compute the number counts in the observed 0.5-2 keV band, finding a decline in the number of sources at z>3 and constraining phenomenological models of X-ray background. We compute the AGN space density at z>3 in two different luminosity bins. At higher luminosities (logL(2-10 keV) > 44.1 erg/s) the space density declines exponentially, dropping by a factor ~20 from z~3 to z~6. The observed decline is ~80% steeper at lower luminosities (43.55 erg/s < logL(2-10 keV) < 44.1 erg/s), from z~3 to z~4.5. We study the space density evolution dividing our sample in optically classified Type 1 and Type 2 AGN. At logL(2-10 keV) > 44.1 erg/s, unobscured and obscured objects may have different evolution with redshift, the obscured component being three times higher at z~5. Finally, we compare our space density with predictions of quasar activation merger models, whose calibration is based on optically luminous AGN. These models significantly overpredict the number of expected AGN at logL(2-10 keV) > 44.1 erg/s with respect to our data.

We present new Keck/MOSFIRE K-band spectroscopy for a sample of 14 faint, X-ray-selected active galactic nuclei (AGNs) in the COSMOS field. The data cover the spectral region surrounding the broad Balmer emission lines, which enables the estimation of black hole masses (M_BH) and accretion rates (in terms of L/L_Edd). We focus on 10 AGN at z~3.3, where we observe the Hbeta spectral region, while for the other four z~2.4 sources we use the Halpha broad emission line. Compared with previous detailed studies of unobscured AGNs at these high redshifts, our sources are fainter by an order of magnitude, corresponding to number densities of order ~10^-6--10^-5 Mpc^-3. The lower AGN luminosities also allow for a robust identification of the host galaxy emission, necessary to obtain reliable intrinsic AGN luminosities, BH masses and accretion rates. We find the AGNs in our sample to be powered by supermassive black holes (SMBHs) with a typical mass of M_BH~5*10^8 M_Sol -- significantly lower than the higher-luminosity, rarer quasars reported in previous studies. The accretion rates are in the range L/L_Edd~0.1-0.5, with an evident lack of sources with lower L/L_Edd (and higher M_BH), as found in several studies of faint AGNs at intermediate redshifts. Based on the early growth expected for the SMBHs in our sample, we argue that a significant population of faint z~6 AGNs, with M_BH~10^6 M_Sol, should be detectable in the deepest X-ray surveys available, but this is not observed. We discuss several possible explanations for the apparent absence of such a population, concluding that the most probable scenario involves an evolution in source obscuration and/or radiative efficiencies.

This is the third paper in a series describing the spectroscopic properties of a sample of 39 AGN at $z \sim 1.5$, selected to cover a large range in black hole mass ($M_{BH}$) and Eddington ratio ($L/L_{Edd}$). In this paper, we continue the analysis of the VLT/X-shooter observations of our sample with the addition of 9 new sources. We use an improved Bayesian procedure, which takes into account intrinsic reddening, and improved $M_{BH}$ estimates, to fit thin accretion disc (AD) models to the observed spectra and constrain the spin parameter ($a_*$) of the central black holes. We can fit 37 out of 39 AGN with the thin AD model, and for those with satisfactory fits, we obtain constraints on the spin parameter of the BHs, with the constraints becoming generally less well defined with decreasing BH mass. Our spin parameter estimates range from $\sim$$-$0.6 to maximum spin for our sample, and our results are consistent with the "spin-up" scenario of BH spin evolution. We also discuss how the results of our analysis vary with the inclusion of non-simultaneous GALEX photometry in our thin AD fitting. Simultaneous spectra covering the rest-frame optical through far-UV are necessary to definitively test the thin AD theory and obtain the best constraints on the spin parameter.

Heavily obscured accretion is believed to represent an important stage in the growth of supermassive black holes, and to play an important role in shaping the observed spectrum of the Cosmic X-ray Background (CXB). Hard X-ray (E$>$10 keV) selected samples are less affected by absorption than samples selected at lower energies, and are therefore one of the best ways to detect and identify Compton-thick (CT, $\log N_{\rm\,H}\geq 24$) Active Galactic Nuclei (AGN). In this letter we present the first results of the largest broad-band (0.3-150 keV) X-ray spectral study of hard X-ray selected AGN to date, focusing on the properties of heavily obscured sources. Our sample includes the 834 AGN (728 non-blazar, average redshift $z\simeq 0.055$) reported in the 70-months catalog of the all-sky hard X-ray Swift/BAT survey. We find 55 CT AGN, which represent $7.6^{+1.1}_{-2.1}\%$ of our non-blazar sample. Of these, 26 are reported as candidate CT AGN for the first time. We correct for selection bias and derive the intrinsic column density distribution of AGN in the local Universe in two different luminosity ranges. We find a significant decrease in the fraction of obscured Compton-thin AGN for increasing luminosity, from $46\pm3\%$ (for $\log L_{\rm\,14-195} = 40-43.7$) to $39\pm3\%$ (for $\log L_{\rm\,14-195} = 43.7-46$). A similar trend is also found for CT AGN. The intrinsic fraction of CT AGN with $\log N_{\rm\,H}=24-25$ normalised to unity in the $\log N_{\rm H} = 20-25$ range is $27\pm4\%$, and is consistent with the observed value obtained for AGN located within 20 Mpc.

This is the second in a series of papers aiming to test how the mass ($M_{\rm BH}$), accretion rate ($\dot{M}$) and spin ($a_{*}$) of super massive black holes (SMBHs) determine the observed properties of type-I active galactic nuclei (AGN). Our project utilizes a sample of 39 unobscured AGN at $z\simeq1.55$ observed by VLT/X-shooter, selected to map a large range in $M_{\rm BH}$ and $L/L_{\rm edd}$ and covers the most prominent UV-optical (broad) emission lines, including H$\alpha$, H$\beta$, MgII, and CIV. This paper focuses on single-epoch, "virial" $M_{\rm BH}$ determinations from broad emission lines and examines the implications of different continuum modeling approaches in line width measurements. We find that using a "local" power-law continuum instead of a physically-motivated thin disk continuum leads to only slight underestimation of the FWHM of the lines and the associated $M_{\rm BH}\left({\rm FWHM}\right)$. However, the line dispersion $\sigma_{\rm line}$ and associated $M_{\rm BH}\left(\sigma_{\rm line}\right)$ are strongly affected by the continuum placement and provides less reliable mass estimates than FWHM-based methods. Our analysis shows that H$\alpha$, H$\beta$ and MgII can be safely used for virial $M_{\rm BH}$ estimation. The CIV line, on the other hand, is not reliable in the majority of the cases, this may indicate that the gas emitting this line is not virialized. While H$\alpha$ and H$\beta$ show very similar line widths, the mean ${\rm FWHM\left(MgII\right)}$ is about 30% narrower than ${\rm FWHM\left(H\beta\right)}$. We confirm several recent suggestions to improve the accuracy in CIV-based mass estimates, relying on other UV emission lines. Such improvements do not reduce the scatter between CIV-based and Balmer-line-based mass estimates.

We report herschel observations of 100 very luminous, optically selected AGNs at z=2-3.5 with log(LUV)(erg/sec)> 46.5, where LUV=L1350A. The distribution in LUV is similar to the general distribution of SDSS AGNs in this redshift and luminosity interval. We measured SF luminosity, LSF, and SFR in 34 detected sources by fitting combined SF and WISE-based torus templates. We also obtained statistically significant stacks for the undetected sources in two luminosity groups. The sample properties are compared with those of very luminous AGNs at z>4.5. The main findings are: 1) The mean and the median SFRs of the detected sources are 1176 and 1010 Msun/yr, respectively. The mean SFR of the undetected sources is 148 Msun/yr. The ratio of SFR to BH accretion rate is approximately 80 for the detected sources and less than 10 for the undetected sources. There is no difference in LAGN and only a very small difference in L(torus) between detected and undetected sources. 2) The redshift distribution of LSF and LAGN for the most luminous, redshift 2-7 AGNs are different. The highest LAGN are found at z=~3. However, LSF of such sources peaks at z=~5. Assuming the objects in our sample are hosted by the most massive galaxies at those redshifts, we find many of them are below the main-sequence of SF galaxies at z=2-3.5. 3) The SEDs of dusty tori at high redshift are similar to those found in low redshift, low luminosity AGNs. Herschel upper limits put strong constraints on the long wavelength SED ruling out several earlier suggested torus templates. 4) We find no evidence for a luminosity dependence of the torus covering factor in sources with log(LAGN)=44-47.5. This conclusion is based on the highly uncertain and non-uniformally treated LAGN in many earlier studies. The median covering factors over this range are 0.68 for isotropic dust emission and 0.4 for anisotropic emission.

We present the catalog of optical and infrared counterparts of the Chandra COSMOS-Legacy Survey, a 4.6 Ms Chandra program on the 2.2 square degrees of the COSMOS field, combination of 56 new overlapping observations obtained in Cycle 14 with the previous C-COSMOS survey. In this Paper we report the i, K, and 3.6 micron identifications of the 2273 X-ray point sources detected in the new Cycle 14 observations. We use the likelihood ratio technique to derive the association of optical/infrared (IR) counterparts for 97% of the X-ray sources. We also update the information for the 1743 sources detected in C-COSMOS, using new K and 3.6 micron information not available when the C-COSMOS analysis was performed. The final catalog contains 4016 X-ray sources, 97% of which have an optical/IR counterpart and a photometric redshift, while 54% of the sources have a spectroscopic redshift. The full catalog, including spectroscopic and photometric redshifts and optical and X-ray properties described here in detail, is available online. We study several X-ray to optical (X/O) properties: with our large statistics we put better constraints on the X/O flux ratio locus, finding a shift towards faint optical magnitudes in both soft and hard X-ray band. We confirm the existence of a correlation between X/O and the the 2-10 keV luminosity for Type 2 sources. We extend to low luminosities the analysis of the correlation between the fraction of obscured AGN and the hard band luminosity, finding a different behavior between the optically and X-ray classified obscured fraction.

We study the correlation between the [O III]$\lambda 5007$ and X-ray luminosities of local Active Galactic Nuclei (AGNs), using a complete, hard X-ray ($>10$ keV) selected sample in the Swift/BAT 9-month catalog. From our optical spectroscopic observations at the South African Astronomical Observatory and the literature, a catalog of [O III]$\lambda 5007$ line flux for all 103 AGNs at Galactic latitudes of $|b|>15^\circ$ is complied. Significant correlations with intrinsic X-ray luminosity ($L_{\rm X}$) are found both for observed ($L_{\rm [O~III]}$) and extinction-corrected ($L_{\rm [O~III]}^{\rm cor}$) luminosities, separately for X-ray unabsorbed and absorbed AGNs. We obtain the regression form of $L_{\rm [O~III]}$ $\propto L_{\rm 2-10\; keV}^{1.18\pm0.07}$ and $L_{\rm [O~III]}^{\rm cor}$ $\propto L_{\rm 2-10\; keV}^{1.16\pm0.09}$ from the whole sample. The absorbed AGNs with low ($<$0.5\%) scattering fractions in soft X-rays show on average smaller $L_{\rm [O~III]}/L_{\rm X}$ and $L_{\rm [O~III]}^{\rm cor}/L_{\rm X}$ ratios than the other absorbed AGNs, while those in edge-on host galaxies do not. These results suggest that a significant fraction of this population are buried in tori with small opening angles. By using these $L_{\rm [O~III]}$ vs. $L_{\rm X}$ correlations, the X-ray luminosity function of local AGNs (including Compton thick AGNs) in a standard population synthesis model gives much better agreement with the [O III]$\lambda 5007$ luminosity function derived from the Sloan Digital Sky Survey than previously reported. This confirms that hard X-ray observations are a very powerful tool to find AGNs with high completeness.

We investigated AGN activity in low-mass galaxies, an important regime that can shed light onto BH formation and evolution, and their interaction with their host galaxies. We identified 336 AGN candidates from a parent sample of $\sim 48,000$ nearby low-mass galaxies ($M_{\rm \star} \leq 10^{9.5}M_\odot$, $z < 0.1$) in the SDSS. We selected the AGN using the classical BPT diagram, a similar optical emission line diagnostic based on the HeII$\lambda$4686 line, and mid-IR color cuts. Different criteria select host galaxies with different physical properties such as stellar mass and optical color, and only 3 out of 336 sources fulfill all three criteria. This could be in part due to selection biases. The resulting AGN fraction of $\sim 0.7 \%$ is at least one order of magnitude below the one estimated for more massive galaxies. At optical wavelengths, the HeII-based AGN selection appears to be more sensitive to AGN hosted in star-forming galaxies than the classical BPT diagram, at least in the low-mass regime. The archival X-ray and radio data available for some of the optically selected AGN candidates seem to confirm their AGN nature, but follow-up observations are needed to confirm the AGN nature of the rest of the sample, especially in the case of mid-IR selection. Our sample will be important for future follow-up studies aiming to understand the relation between BHs and host galaxies in the low-mass regime.

We investigate the relationship between X-ray and optical line emission in 340 nearby AGN selected above 10 keV using Swift BAT. We find a weak correlation between the extinction corrected [O III] and hard X-ray luminosity (14-195 keV) with a [OIII] large scatter (R_Pear = 0.64, sigma = 0.62 dex) and a similarly large scatter with the intrinsic 2-10 keV to [O III] luminosities (RPear=0.63, sigma = 0.63 dex). Correlations of the hard X-ray fluxes with the fluxes of high-ionization narrow lines ([O III], He II, [Ne III] and [Ne V]) are not significantly better than with the low ionization lines (Halpha, [SII]). Factors like obscuration or physical slit size are not found to be a significant part of the large scatter. In contrast, the optical emission lines show much better correlations with each other (sigma = 0.3 dex) than with the X-ray flux. The inherent large scatter questions the common usage of narrow emission lines as AGN bolometric luminosity indicators and suggests that other issues such as geometrical differences in the scattering of the ionized gas or long term AGN variability are important.

We explore the connections between the evolving galaxy and AGN populations. We present a simple phenomenological model that links the evolving galaxy mass function and the evolving quasar luminosity function, which makes specific and testable predictions for the distribution of host galaxy masses for AGN of different luminosities. We show that the $\phi^{*}$ normalisations of the galaxy mass function and of the AGN luminosity function closely track each other over a wide range of redshifts, implying a constant "duty cycle" of AGN activity. The strong redshift evolution in the AGN $L^*$ can be produced by either an evolution in the distribution of Eddington ratios, or in the $m_{bh}/m_{*}$ mass ratio, or both. To try to break this degeneracy we look at the distribution of AGN in the SDSS ($m_{bh},L$) plane, showing that an evolving ratio $m_{bh}/m_{*} \propto (1+z)^2$ reproduces the observed data and also reproduces the local relations which connect the black hole population with the host galaxies for both quenched and star-forming populations. We stress that observational studies that compare the masses of black holes in active galaxies at high redshift with those in quiescent galaxies locally will always see much weaker evolution. Evolution of this form would produce, or could be produced by, a redshift-independent $m_{bh} - \sigma$ relation and could explain why the local $m_{bh} - \sigma$ relation is tighter than $m_{bh} - m_{*}$ even if $\sigma$ is not directly linked to black hole growth. Irrespective of the evolution of $m_{bh}/m_{*}$, the model reproduces both the appearance of "downsizing" and the so-called "sub-Eddington boundary" without any mass-dependence in the evolution of black hole growth rates.

Supermassive black holes (SMBHs) and their host galaxies are generally thought to coevolve, so that the SMBH achieves up to about 0.2 to 0.5% of the host galaxy mass in the present day. The radiation emitted from the growing SMBH is expected to affect star formation throughout the host galaxy. The relevance of this scenario at early cosmic epochs is not yet established. We present spectroscopic observations of a galaxy at redshift z = 3.328, which hosts an actively accreting, extremely massive BH, in its final stages of growth. The SMBH mass is roughly one-tenth the mass of the entire host galaxy, suggesting that it has grown much more efficiently than the host, contrary to models of synchronized coevolution. The host galaxy is forming stars at an intense rate, despite the presence of a SMBH-driven gas outflow.

We have newly identified a substantial number of type 1 active galactic nuclei (AGN) featuring weak broad-line regions (BLRs) at z < 0.2 from detailed analysis of galaxy spectra in the Sloan Digital Sky Survey Data Release 7. These objects predominantly show a stellar continuum but also a broad H-alpha emission line, indicating the presence of a low-luminosity AGN oriented so that we are viewing the central engine directly without significant obscuration. These accreting black holes have previously eluded detection due to their weak nature. The new BLR AGNs we found increased the number of known type 1 AGNs by 49%. Some of these new BLR AGNs were detected at the Chandra X-ray Observatory, and their X-ray properties confirm that they are indeed type 1 AGN. Based on our new and more complete catalogue of type 1 AGNs, we derived the type 1 fraction of AGNs as a function of [OIII] 5007 emission luminosity and explored the possible dilution effect on the obscured AGN due to star-formation. The new type 1 AGN fraction shows much more complex behavior with respect to black hole mass and bolometric luminosity than suggested by the existing receding torus model. The type 1 AGN fraction is sensitive to both of these factors, and there seems to be a sweet spot (ridge) in the diagram of black hole mass and bolometric luminosity. Furthermore, we present a hint that the Eddington ratio plays a role in determining the opening angles.

Over the past 15 years, examples of exotic radio-quiet quasars with intrinsically weak or absent broad emission line regions (BELRs) have emerged from large-scale spectroscopic sky surveys. Here, we present spectroscopy of seven such weak emission line quasars (WLQs) at moderate redshifts (z=1.4-1.7) using the X-shooter spectrograph, which provides simultaneous optical and near-infrared spectroscopy covering the rest-frame ultraviolet through optical. These new observations effectively double the number of WLQs with spectroscopy in the optical rest-frame, and they allow us to compare the strengths of (weak) high-ionization emission lines (e.g., CIV) to low-ionization lines (e.g., MgII, Hb, Ha) in individual objects. We detect broad Hb and Ha emission in all objects, and these lines are generally toward the weaker end of the distribution expected for typical quasars (e.g., Hb has rest-frame equivalent widths ranging from 15-40 Ang.). However, these low-ionization lines are not exceptionally weak, as is the case for high-ionization lines in WLQs. The X-shooter spectra also display relatively strong optical FeII emission, Hb FWHM <4000 km/s, and significant CIV blueshifts (1000-5500 km/s) relative to the systemic redshift; two spectra also show elevated ultraviolet FeII emission, and an outflowing component to their (weak) MgII emission lines. These properties suggest that WLQs are exotic versions of "wind-dominated" quasars. Their BELRs either have unusual high-ionization components, or their BELRs are in an atypical photoionization state because of an unusually soft continuum.

The radius and surface composition of an exploding massive star,as well as the explosion energy per unit mass, can be measured using early UV observations of core collapse supernovae (SNe). We present the first results from a simultaneous GALEX/PTF search for early UV emission from SNe. Six Type II SNe and one Type II superluminous SN (SLSN-II) are clearly detected in the GALEX NUV data. We compare our detection rate with theoretical estimates based on early, shock-cooling UV light curves calculated from models that fit existing Swift and GALEX observations well, combined with volumetric SN rates. We find that our observations are in good agreement with calculated rates assuming that red supergiants (RSGs) explode with fiducial radii of 500 solar, explosion energies of 10^51 erg, and ejecta masses of 10 solar masses. Exploding blue supergiants and Wolf-Rayet stars are poorly constrained. We describe how such observations can be used to derive the progenitor radius, surface composition and explosion energy per unit mass of such SN events, and we demonstrate why UV observations are critical for such measurements. We use the fiducial RSG parameters to estimate the detection rate of SNe during the shock-cooling phase (<1d after explosion) for several ground-based surveys (PTF, ZTF, and LSST). We show that the proposed wide-field UV explorer ULTRASAT mission, is expected to find >100 SNe per year (~0.5 SN per deg^2), independent of host galaxy extinction, down to an NUV detection limit of 21.5 mag AB. Our pilot GALEX/PTF project thus convincingly demonstrates that a dedicated, systematic SN survey at the NUV band is a compelling method to study how massive stars end their life.

We investigate early black hole (BH) growth through the methodical search for $z\gtrsim5$ AGN in the $Chandra$ Deep Field South. We base our search on the $Chandra$ 4-Ms data with flux limits of $9.1\times\ 10^{-18}$ (soft, 0.5 - 2 keV) and $5.5\times\ 10^{-17}\ \mathrm{erg}\ \mathrm{s}^{-1}\ \mathrm{cm}^{-2}$ (hard, 2 - 8 keV). At $z\sim5$ this corresponds to luminosities as low as $\sim10^{42}$ ($\sim10^{43}$) $\mathrm{erg}\ \mathrm{s}^{-1}$ in the soft (hard) band and should allow us to detect Compton-thin AGN with $M_\mathrm{BH}>10^7 M_{\odot}$ and Eddington ratios > 0.1. Our field ($0.03~\mathrm{deg}^2$) contains over 600 $z\sim5$ Lyman Break Galaxies. Based on lower redshift relations we would expect $\sim20$ of them to host AGN. After combining the $Chandra$ data with GOODS/ACS, CANDELS/WFC3 and $Spitzer$/IRAC data, the sample consists of 58 high-redshift candidates. We run a photometric redshift code, stack the GOODS/ACS data, apply colour criteria and the Lyman Break Technique and use the X-ray Hardness Ratio. We combine our tests and using additional data find that all sources are most likely at low redshift. We also find five X-ray sources without a counterpart in the optical or infrared which might be spurious detections. We conclude that our field does not contain any convincing $z\gtrsim5$ AGN. Explanations for this result include a low BH occupation fraction, a low AGN fraction, short, super-Eddington growth modes, BH growth through BH-BH mergers or in optically faint galaxies. By searching for $z\gtrsim5$ AGN we are setting the foundation for constraining early BH growth and seed formation scenarios.

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.

While performing a survey to detect recoiling supermassive black holes, we have identified an unusual source having a projected offset of 800 pc from a nearby dwarf galaxy. The object, SDSS J113323.97+550415.8, exhibits broad emission lines and strong variability. While originally classified as a supernova (SN) because of its nondetection in 2005, we detect it in recent and past observations over 63 yr and find over a magnitude of rebrightening in the last 2 years. Using high-resolution adaptive optics observations, we constrain the source emission region to be <12 pc and find a disturbed host-galaxy morphology indicative of recent merger activity. Observations taken over more than a decade show narrow [O III] lines, constant ultraviolet emission, broad Balmer lines, a constant putative black hole mass over a decade of observations despite changes in the continuum, and optical emission-line diagnostics consistent with an active galactic nucleus (AGN). However, the optical spectra exhibit blueshifted absorption, and eventually narrow Fe II and [Ca II] emission, each of which is rarely found in AGN spectra. While this peculiar source displays many of the observational properties expected of a potential black hole recoil candidate, some of the properties could also be explained by a luminous blue variable star (LBV) erupting for decades since 1950, followed by a Type IIn SN in 2001. Interpreted as an LBV followed by a SN analogous to SN 2009ip, the multi-decade LBV eruptions would be the longest ever observed, and the broad Halpha emission would be the most luminous ever observed at late times (>10 yr), larger than that of unusually luminous supernovae such as SN 1988Z, suggesting one of the most extreme episodes of pre-SN mass loss ever discovered.

The physics of active super massive black holes (BHs) is governed by their mass (M_BH), spin (a*) and accretion rate ($\dot{M}$). This work is the first in a series of papers with the aim of testing how these parameters determine the observable attributes of active galactic nuclei (AGN). We have selected a sample in a narrow redshift range, centered on z~1.55, that covers a wide range in M_BH and $\dot{M}$, and are observing them with X-shooter, covering rest wavelengths ~1200-9800 \AA. The current work covers 30 such objects and focuses on the origin of the AGN spectral energy distribution (SED). After estimating M_BH and $\dot{M}$ based on each observed SED, we use thin AD models and a Bayesian analysis to fit the observed SEDs in our sample. We are able to fit 22/30 of the SEDs. Out of the remaining 8 SEDs, 3 can be fit by the thin AD model by correcting the observed SED for reddening within the host galaxy and 4 can be fit by adding a disc wind to the model. In four of these 8 sources, Milky Way-type extinction, with the strong 2175\AA\ feature, provides the best reddening correction. The distribution in spin parameter covers the entire range, from -1 to 0.998, and the most massive BHs have spin parameters greater than 0.7. This is consistent with the "spin-up" model of BH evolution. Altogether, these results indicate that thin ADs are indeed the main power houses of AGN, and earlier claims to the contrary are likely affected by variability and a limited observed wavelength range.

The radiative efficiencies ($\eta$) of 72 luminous unobscured Active Galactic Nuclei (AGNs) at $z\sim1.5-3.5$, powered by some of the most massive black holes (BHs), are constrained. The analysis is based on accretion disk (AD) models, which link the continuum luminosity at rest-frame optical wavelengths and the black hole mass ($M_{\rm BH}$) to the accretion rate through the AD, $\dot{M}_{\rm AD}$. The data are gathered from several literature samples with detailed measurements of the ${\rm H}\beta$ emission line complex, observed at near-IR bands. When coupled with standard estimates of bolometric luminosities ($L_{\rm bol}$), the analysis suggests high radiative efficiencies, with most of the sources showing $\eta>0.2$ - that is, higher than the commonly assumed value of 0.1, and the expected value for non-spinning BHs ($\eta=0.057$). Even under more conservative assumptions regarding $L_{\rm bol}$ (i.e., $L_{\rm bol}=3\times \lambda L_{\lambda}$[5100$\AA$]), most of the extremely massive BHs in the sample (i.e., $M_{\rm BH} > 3\times10^9\,M_{\rm \odot}$) show radiative efficiencies which correspond to very high BH spins ($a_{\rm *}$), with typical values well above $a_{\rm *}\simeq0.7$. These results stand in contrast to the predictions of a "spin-down"scenario, in which a series of randomly-oriented accretion episodes lead to $a_{\rm *}\simeq0$. Instead, the analysis presented here strongly supports a "spin-up" scenario, which is driven by either prolonged accretion or a series of anisotropically-oriented accretion episodes. Considering the fact that these extreme BHs require long-duration continuous accretion to account for their high masses, it is argued that the most probable scenario for the SMBHs under study is that of an almost continuous sequence of randomly- yet not isotropically-oriented accretion episodes.

We report Herschel Spitzer and Wise observations of 44 z=4.8 optically selected active galactic nuclei (AGNs). This sample contains the highest mass black holes (BHs) at this redshift. Ten of the objects were detected by Herschel and five show emission that is not associated with the AGNs. The star formation (SF) luminosity (LSF) obtained by fitting the spectral energy distribution (SED) with standard SF templates is in the range 10^(46.62-47.21) erg/s corresponding to SF rates of 1090-4240 Msun/yr. Fitting with luminous submillimeter galaxy SEDs gives SF rates that are smaller by 0.05 dex when using all bands and 0.1 dex when ignoring the 250-mic band. A 40K gray-body fits to only the 500-mic fluxes reduce LSF by about a factor two. Stacking analysis of 29 undetected sources give significant signals in all three bands. A SF template fit indicates LSF=10^(46.19-46.23) erg/s depending on the assume AGN contribution. A 40K fit to the stacked 500-mic flux gives LSF=10^(45.95) erg/s. The mean BH mass (mbh) and AGN luminosity (LAGN) of the detected sources are significantly higher than those of the undetected ones. The spectral differences are seen all the way from UV to far infrared wavelengths. The mean optical-UV spectra are similar to the ones predicted for thin accretion disks around BHs with the measured masses and accretion rates. We suggest two alternative explanations to the correlation of LSF, LAGN and mbh, one involving no AGN feedback and the second a moderate feedback that affects, but not totally quench SF in 3/4 of the sources. We compare our LSF and LAGN to lower redshift samples and show a new correlation between LSF and mbh. We also examine several rather speculative ideas about the host galaxy properties including the possibility that the detected sources are above the SF mass sequence (MS) at z=4.8, perhaps in mergers, and most of the undetected sources are on the MS.

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.

We investigate the accretion rate, bolometric luminosity, black hole (BH) growth time and BH spin in a large AGN sample under the assumption that all such objects are powered via thin or slim accretion discs (ADs). We use direct estimates of the mass accretion rate, Mdot, to show that many currently used values of Lbol and Ledd are either under estimated or over estimated because they are based on bolometric correction factors that are adjusted to the properties of moderately accreting active galactic nuclei (AGN) and do not take into account the correct combination of BH mass, spin and accretion rate. The consistent application of AD physics to our sample of Sloan Digital Sky Survey (SDSS) AGN leads to the following findings: 1. Even the most conservative assumption about the radiative efficiency of fast accreting BHs shows that many of these sources must contain slim ADs. We illustrate this by estimating the fraction of such objects at various redshifts. 2. Many previously estimated BH growth times are inconsistent with the AD theory. In particular, the growth times of the fastest accreting BHs were over estimated in the past by large factors with important consequences to AGN evolution. 3. Currently used bolometric correction factors for low accretion rate very massive SDSS BHs, are inconsistent with the AD theory. Applying the AD set of assumptions to such objects, combined with standard photoionization calculations of broad emission lines, leads to the conclusion that many such objects must contain fast spinning BHs.

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.

The time-variable electromagnetic sky has been well-explored at a wide range of wavelengths. Numerous high-energy space missions take advantage of the dark Gamma-ray and X-ray sky and utilize very wide field detectors to provide almost continuous monitoring of the entire celestial sphere. In visible light, new wide-field ground-based surveys cover wide patches of sky with ever decreasing cadence, progressing from monthly-weekly time scale surveys to sub-night sampling. In the radio, new powerful instrumentation offers unprecedented sensitivity over wide fields of view, with pathfinder experiments for even more ambitious programs underway. In contrast, the ultra-violet (UV) variable sky is relatively poorly explored, even though it offers exciting scientific prospects. Here, we review the potential scientific impact of a wide-field UV survey on the study of explosive and other transient events, as well as known classes of variable objects, such as active galactic nuclei and variable stars. We quantify our predictions using a fiducial set of observational parameters which are similar to those envisaged for the proposed ULTRASAT mission. We show that such a mission would be able to revolutionize our knowledge about massive star explosions by measuring the early UV emission from hundreds of events, revealing key physical parameters of the exploding progenitor stars. Such a mission would also detect the UV emission from many tens of tidal-disruption events of stars by super massive black holes at galactic nuclei and enable a measurement of the rate of such events. The overlap of such a wide-field UV mission with existing and planned gravitational-wave and high-energy neutrino telescopes makes it especially timely.

We analyze several large samples of AGN in order to establish the best tools required to study the evolution of black hole mass (M_BH) and normalized accretion rate (L/L_Edd). The data include spectra from the SDSS, 2QZ and 2SLAQ surveys at z<2, and a compilation of smaller samples with 0<z<5. We critically evaluate the usage of the MgII(2798A) and CIV(1549A) lines, and adjacent continuum bands, as estimators of M_BH and L/L_Edd, by focusing on sources where one of these lines is observed together with Hbeta. We present a new, luminosity-dependent bolometric correction for the monochromatic luminosity at 3000A, L_3000, which is lower by a factor of 1.75 than those used in previous studies. We also re-calibrate the use of L_3000 as an indicator for the size of the broad emission line region (R_BLR) and find that R_BLR is proportional to L_3000^0.62. We find that FWHM(MgII)\simeqFWHM(Hb) for all sources with FWHM(MgII)<6000 km/s. Beyond this FWHM, the MgII line width seems to saturate. The spectral region of the MgII line can thus be used to reproduce Hb-based estimates of M_BH and L/L_Edd, with negligible systematic differences and a scatter of 0.3 dex. The width of the CIV line, on the other hand, shows no correlation with either that of the Hb or the MgII lines and we could not identify the reason for this discrepancy. The scatter of M_BH(CIV), relative to M_BH(Hb) is of almost 0.5 dex. Moreover, 46% of the sources have FWHM(CIV)<FWHM(Hb), in contrast with the basic premise of the virial method, which predicts FWHM(CIV)/FWHM(Hbeta)~1.9. Thus, the CIV line cannot be used to obtain precise estimates of M_BH. We conclude by presenting the observed evolution of M_BH and L/L_Edd with cosmic epoch. The steep rise of L/L_Edd with redshift up to z~1 flattens towards the expected maximal value of L/L_Edd~1, with lower-M_BH sources showing higher values of L/L_Edd at all redshifts. [Abridged]

We report new Herschel observations of 25 z=4.8 extremely luminous optically selected active galactic nuclei (AGNs). Five of the sources have extremely large star forming (SF) luminosities, L_SF, corresponding to SF rates (SFRs) of 2800-5600 M_sol/yr assuming a Salpeter IMF. The remaining sources have only upper limits on their SFRs but stacking their Herschel images results in a mean SFR of 700 +/- 150 M_sol/yr. The higher SFRs in our sample are comparable to the highest observed values so far, at any redshift. Our sample does not contain obscured AGNs, which enables us to investigate several evolutionary scenarios connecting super-massive black holes and SF activity in the early universe. The most probable scenario is that we are witnessing the peak of SF activity in some sources and the beginning of the post-starburst decline in others. We suggest that all 25 sources, which are at their peak AGN activity, are in large mergers. AGN feedback may be responsible for diminishing the SF activity in 20 of them but is not operating efficiently in 5 others.

We queried the Spitzer archive for high-resolution observations with the Infrared Spectrograph of optically selected active galactic nuclei (AGN) for the purpose of identifying sources with resolved fine-structure lines that would enable studies of the narrow-line region (NLR) at mid-infrared wavelengths. By combining 298 Spitzer spectra with 6 Infrared Space Observatory spectra, we present kinematic information of the NLR for 81 z<=0.3 AGN. We used the [NeV], [OIV], [NeIII], and [SIV] lines, whose fluxes correlate well with each other, to probe gas photoionized by the AGN. We found that the widths of the lines are, on average, increasing with the ionization potential of the species that emit them. No correlation of the line width with the critical density of the corresponding transition was found. The velocity dispersion of the gas, sigma, is systematically higher than that of the stars, sigma_*, in the AGN host galaxy, and it scales with the mass of the central black hole, M_BH. Further correlations between the line widths and luminosities L, and between L and M_BH, are suggestive of a three dimensional plane connecting log(M_BH) to a linear combination of log(sigma) and log(L). Such a plane can be understood within the context of gas motions that are driven by AGN feedback mechanisms, or virialized gas motions with a power-law dependence of the NLR radius on the AGN luminosity. The M_BH estimates obtained for 35 type 2 AGN from this plane are consistent with those obtained from the M_BH-sigma_* relation.

We fitted the optical to mid-infrared (MIR) spectral energy distributions (SEDs) of ~15000 type-I, 0.75<z<2, active galactic nuclei (AGNs) in an attempt to constrain the properties of the physical component responsible for the rest-frame near-infrared (NIR) emission. We combine optical spectra from the Sloan Digital Sky Survey (SDSS) and MIR photometry from the preliminary data release of the Wide Infrared Survey Explorer (WISE). The sample spans a large range of AGN properties: luminosity, black hole mass, and accretion rate. Our model has two components: a UV-optical continuum source and very hot, pure-graphite dust clouds. We present the luminosity of the hot-dust component and its covering factor, for all sources, and compare it with the intrinsic AGN properties. We find that the hot-dust component is essential to explain the (rest) NIR emission in almost all AGNs in our sample, and that it is consistent with clouds containing pure-graphite grains and located between the dust-free broad line region (BLR) and the "standard" torus. The covering factor of this component has a relatively narrow distribution around a peak value of ~0.13, and it correlates with the AGN bolometric luminosity. We suggest that there is no significant correlation with either black hole mass or normalized accretion rate. The fraction of hot-dust-poor AGNs in our sample is ~15-20%, consistent with previous studies. We do not find a dependence of this fraction on redshift or source luminosity.

We present new Gemini-North/NIRI and VLT/SINFONI H-band spectroscopy for a flux limited sample of 40 z~4.8 active galactic nuclei, selected from the Sloan Digital Sky Survey. The sample probably contains the most massive active black holes (BHs) at this redshift and spans a broad range in bolometric luminosity, 2.7x10^46< L_bol < 2.4x10^47 erg/sec. The high-quality observations and the accurate fitting of the MgII(2800A) line, enable us to study, systematically, the distribution of BH mass (M_BH) and normalized accretion rate (L/L_Edd) at z~4.8. We find that 10^8 < M_BH < 6.6x10^9 M_sun, with a median of ~8.4x10^8 M_sun. We also find that 0.2 < L/L_Edd < 3.9 with a median of ~0.6. Most of these sources had enough time to grow to their observed mass at z~4.8 from z=20, assuming a range of seed BH masses, with ~40% that are small enough to be stellar remnants. Compared to previously studied samples at z~2.4 and 3.3, the masses of the z~4.8 BHs are typically lower by ~0.5 dex. and their L/L_Edd is higher by a similar factor. The new z~4.8 sample can be considered as the progenitor population of the most massive BHs at z~2.4 and 3.3. Such an evolutionary interpretation requires that the growth of the BHs from z~4.8 to z~3.3 and z~2.4 proceeds with short duty cycles, of about 10-20%, depending on the particular growth scenario.

We present Gemini-North K-band spectra of two representative members of the class of high-redshift quasars with exceptionally weak rest-frame ultraviolet emission lines (WLQs), SDSS J114153.34+021924.3 at z=3.55 and SDSS J123743.08+630144.9 at z=3.49. In both sources we detect an unusually weak broad H_beta line and we place tight upper limits on the strengths of their [O III] lines. Virial, H_beta-based black-hole mass determinations indicate normalized accretion rates of L/L_Edd=0.4 for these sources, which is well within the range observed for typical quasars with similar luminosities and redshifts. We also present high-quality XMM-Newton imaging spectroscopy of SDSS J114153.34+021924.3 and find a hard-X-ray photon index of Gamma=1.91^{+0.24}_{-0.22} which supports the virial L/L_Edd determination in this source. Our results suggest that the weakness of the broad-emission lines in WLQs is not a consequence of an extreme continuum-emission source but instead due to abnormal broad-emission line region properties.

We propose a novel method to estimate M_*/M_BH, the ratio of stellar mass (M_*) to black hole mass (M_BH), at various redshifts using two recent observational results: the correlation between the bolometric luminosity of active galactic nuclei (AGN) and the star formation rate (SFR) in their host galaxies, and the correlation between SFR and M_* in star-forming (SF) galaxies. Our analysis is based on M_BH and L_bol measurements in two large samples of type-I AGN at z~1 and z~2, and the measurements of M_*/M_BH in 0.05<z<0.2 red galaxies. We find that M_*/M_BH depends on M_BH at all redshifts. At z~2, M_*/M_BH 280 and ~40 for M_BH=10^8 and M_BH=10^9 M_sol, respectively. M_*/M_BH grows by a factor of ~4-8 from z~2 to z~0 with extreme cases that are as large as 10-20. The evolution is steeper than reported in other studies, probably because we treat only AGN in SF hosts. We caution that estimates of M_*/M_BH evolution which ignore the dependence of this ratio on M_BH can lead to erroneous conclusions.

We present new low resolution Spitzer mid-infrared spectroscopy of a sample of 20 ROSAT selected local Narrow Line Seyfert 1 galaxies (NLS1s). We detect strong AGN continuum in all and clear PAH emission in 70% of the sources. The 6.2 micron PAH luminosity spans three orders of magnitudes, from ~10^(39) erg/s to ~10^(42) erg/s providing strong evidence for intense ongoing star formation in the circumnuclear regions of these sources. Using the IRS/Spitzer archive we gather a large number of additional NLS1s and their broad line counterparts (BLS1s) and constructed NLS1 and BLS1 sub-samples to compare them in various ways. The comparison shows a clear separation according to FWHM(H_beta) such that objects with narrower broad H_beta lines are the strongest PAH emitters. We test this division in various ways trying to remove biases due to luminosity and aperture size. Specifically, we find that star formation activity around NLS1 AGN is larger than around BLS1 of the same AGN luminosity. The above result seems to hold over the entire range of distance and luminosity. Moreover the star formation rate is higher in low black hole mass and high L/L_Edd systems indicating that black hole growth and star formation are occurring simultaneously.

We present new H and K bands spectroscopy of 15 high luminosity active galactic nuclei (AGNs) at redshifts 2.3-3.4 obtained on Gemini South. We combined the data with spectra of additional 29 high-luminosity sources to obtain a sample with 10^{45.2}<\lambda L_{\lambda}(5100A)<10^{47.3} ergs/sec and black hole (BH) mass range, using reverberation mapping relationships based on the H_beta method, of 10^{8.8}-10^{10.7} M_sun. We do not find a correlation of L/L_Edd with M_BH but find a correlation with \lambda L_{\lambda}(5100A) which might be due to selection effects. The L/L_Edd distribution is broad and covers the range ~0.07-1.6, similar to what is observed in lower redshift, lower luminosity AGNs. We suggest that this consistently measured and calibrated sample gives the best representation of L/L_Edd at those redshifts and note potential discrepancies with recent theoretical and observational studies. The lower accretion rates are not in accord with growth scenarios for BHs at such redshifts and the growth times of many of the sources are longer than the age of the universe at the corresponding epochs. This suggests earlier episodes of faster growth at z>~3 for those sources. The use of the C IV method gives considerably different results and a larger scatter; this method seems to be a poor M_BH and L/L_Edd estimator at very high luminosity.

We present line and continuum measurements for 9818 SDSS type-I active galactic nuclei (AGNs) with z le 0.75. The data are used to study the four dimensional space of black hole mass, normalized accretion rate (Ledd), metalicity and redshift. The main results are: 1. Ledd is smaller for larger mass black holes at all redshifts. 2. For a given black hole mass Ledd propto z^gamma or (1+z)^delta where the slope gamma increases with black hole mass. The mean slope is similar to the star formation rate slope over the same redshift interval. 3. The FeII/Hb line ratio is significantly correlated with Ledd. It also shows a weaker negative dependence on redshift. Combined with the known dependence of metalicity on accretion rate, we suggest that the FeII/Hb line ratio is a metalicity indicator. 4. Given the measured accretion rates, the growth times of most AGNs exceed the age of the universe. This suggests past episodes of faster growth for all those sources. Combined with the FeII/Hb result, we conclude that the broad emission lines metalicity goes through cycles and is not a monotonously decreasing function of redshift. 5. FWHM(OIII) is a poor proxy of sigma_* especially for high Ledd. 6. We define a group of narrow line type-I AGNs (NLAGN1s) by their luminosity (or mass) dependent Hb line width. Such objects have Ledd>0.25 and they comprise 8% of the type-I population. Other interesting results include negative Baldwin relationships for EW(Hb) and EW(FeII) and a relative increase of the red part of the Hb line with luminosity.

We combine emission line and X-ray luminosities for 45 sources from the Chandra Deep Field South (CDF-S), and seven HELLAS sources, to obtain a new sample of 52 X-ray selected type-II active galactic nuclei (AGNs). Eighteen of our sources are very luminous with a typical, absorption-corrected 2-10 keV luminosity of few 10^{44} erg/sec (type-II QSOs). We compare the emission line properties of the new sources with emission line and X-ray luminosities of known low redshift, mostly lower luminosity AGNs by using a composite spectrum. We find that ${\rm L_{[OIII]}}$/${\rm L_{2-10}}$ and ${\rm L_{[OII]}}$/${\rm L_{2-10}}$ decrease with L(2-10 keV) such that ${\rm L_{[OIII]}}/{\rm L_{2-10}} \propto {\rm L_{2-10}}^{-0.42}$. The trend was already evident, yet neglected in past low redshift samples. This lead to erroneous calibration of the line-to-X-ray luminosity in earlier AGN samples. The analysis of several type-I samples shows the same trend with a similar slope but a median ${\rm L_{[OIII]}}$/${\rm L_{2-10}}$ which is larger by a factor of about two compared with optically selected type-II samples. We interpret this shift as due to additional reddening in type-II sources and comment in general on the very large extinction in many type-II objects and the significantly smaller average reddening of the SDSS type-II AGNs. The decrease of ${\rm L_{[OIII]}}$/${\rm L_{2-10}}$ with L(2--10 keV) is large enough to suggest that a significant fraction of high luminosity high redshift type-II AGNs have very weak emission lines that may have escaped detection in large samples. A related decrease of EW([OIII]) with optical continuum luminosity is demonstrated by an analysis of 12,000 type-I SDSS AGNs.