Search

We present measurements of broad emission lines and virial estimates of supermassive black hole masses ($M_{BH}$) for a large sample of ultra-hard X-ray selected active galactic nuclei (AGNs) as part of the second data release of the BAT AGN Spectroscopic Survey (BASS/DR2). Our catalog includes $M_{BH}$ estimates for a total 689 AGNs, determined from the H$\alpha$, H$\beta$, $MgII\lambda2798$, and/or $CIV\lambda1549$ broad emission lines. The core sample includes a total of 512 AGNs drawn from the 70-month Swift/BAT all-sky catalog. We also provide measurements for 177 additional AGNs that are drawn from deeper Swift/BAT survey data. We study the links between $M_{BH}$ estimates and line-of-sight obscuration measured from X-ray spectral analysis. We find that broad H$\alpha$ emission lines in obscured AGNs ($\log (N_{\rm H}/{\rm cm}^{-2})> 22.0$) are on average a factor of $8.0_{-2.4}^{+4.1}$ weaker, relative to ultra-hard X-ray emission, and about $35_{-12}^{~+7}$\% narrower than in unobscured sources (i.e., $\log (N_{\rm H}/{\rm cm}^{-2}) < 21.5$). This indicates that the innermost part of the broad-line region is preferentially absorbed. Consequently, current single-epoch $M_{BH}$ prescriptions result in severely underestimated ($>$1 dex) masses for Type 1.9 sources (AGNs with broad H$\alpha$ but no broad H$\beta$) and/or sources with $\log (N_{\rm H}/{\rm cm}^{-2}) > 22.0$. We provide simple multiplicative corrections for the observed luminosity and width of the broad H$\alpha$ component ($L[{\rm b}{\rm H}\alpha]$ and FWHM[bH$\alpha$]) in such sources to account for this effect, and to (partially) remedy $M_{BH}$ estimates for Type 1.9 objects. As key ingredient of BASS/DR2, our work provides the community with the data needed to further study powerful AGNs in the low-redshift Universe., Comment: published in ApJS

We present the second catalog and data release of optical spectral line measurements and AGN demographics of the BAT AGN Spectroscopic Survey, which focuses on the of Swift-BAT hard X-ray detected AGNs. We use spectra from dedicated campaigns and publicly available archives to investigate spectral properties of most of the AGNs listed in the 70-month Swift-BAT all-sky catalog; specifically, 743 of the 746 unbeamed and unlensed AGNs (99.6%). We find a good correspondence between the optical emission line widths and the hydrogen column density distributions using the X-ray spectra, with a clear dichotomy of AGN types for NH = 10^22 cm-2. Based on optical emission-line diagnostics, we show that 48%-75% of BAT AGNs are classified as Seyfert, depending on the choice of emission lines used in the diagnostics. The fraction of objects with upper limits on line emission varies from 6% to 20%. Roughly 4% of the BAT AGNs have lines too weak to be placed on the most commonly used diagnostic diagram, [O III]{\lambda}5007/H\b{eta} versus [N II]{\lambda}6584/H{\alpha}, despite the high signal-to-noise ratio (S/N) of their spectra. This value increases to 35% in the [O III]{\lambda}5007/[O II]{\lambda}3727 diagram, owing to difficulties in line detection. Compared to optically-selected narrow-line AGNs in the Sloan Digital Sky Survey, the BAT narrow-line AGNs have a higher rate of reddening/extinction, with H{\alpha}/H\b{eta} > 5 (~ 36%), indicating that hard X-ray selection more effectively detects obscured AGNs from the underlying AGN population. Finally, we present a subpopulation of AGNs that feature complex broad-lines (34%, 250/743) or double-peaked narrow emission lines (2%, 17/743)., Comment: 30 pages, 17 figures, 11 tables, accepted for publication in ApJS; part of BASS DR2 special issue

We determine the low-redshift X-ray luminosity function (XLF), active black hole mass function (BHMF), and Eddington-ratio distribution function (ERDF) for both unobscured (Type 1) and obscured (Type 2) active galactic nuclei (AGN) using the unprecedented spectroscopic completeness of the BAT AGN Spectroscopic Survey (BASS) data release 2. In addition to a straightforward 1/Vmax approach, we also compute the intrinsic distributions, accounting for sample truncation by employing a forward modeling approach to recover the observed BHMF and ERDF. As previous BHMFs and ERDFs have been robustly determined only for samples of bright, broad-line (Type 1) AGNs and/or quasars, ours is the first directly observationally constrained BHMF and ERDF of Type 2 AGN. We find that after accounting for all observational biases, the intrinsic ERDF of Type 2 AGN is significantly skewed towards lower Eddington ratios than the intrinsic ERDF of Type 1 AGN. This result supports the radiation-regulated unification scenario, in which radiation pressure dictates the geometry of the dusty obscuring structure around an AGN. Calculating the ERDFs in two separate mass bins, we verify that the derived shape is consistent, validating the assumption that the ERDF (shape) is mass independent. We report the local AGN duty cycle as a function of mass and Eddington ratio, by comparing the BASS active BHMF with the local mass function for all SMBH. We also present the log N-log S of Swift-BAT 70-month sources., Comment: Accepted by APJS

We examine the contribution of high-redshift (z>6) active galactic nuclei (AGNs) to cosmic hydrogen reionization, by tracing the growth and ionizing output of the first generation of supermassive black holes (SMBHs). Our calculations are anchored to the observed population of z~6 quasars, and trace back the evolving spectral energy distributions (SEDs) of the accretion flows that power these early AGNs and consider a variety of growth histories, including super-Eddington accretion. Compared to a fixed-shape SED, the evolving thin disks can produce ionizing radiation that is higher by up to ~80%. Across a variety of SMBH growth scenarios, the contribution of AGNs to reionization is limited to late epochs (z<7), and remains sub-dominant compared to star-forming galaxies. This conclusion holds irrespective of the (still unknown) space density of low-luminosity z=6 AGNs, and for growth scenarios that allow super-Eddington accretion. The contribution of AGNs to reionization can extend to earlier epochs (z>8) in scenarios with relative slow SMBH mass growth, i.e., for low accretion rates and/or high spins. We finally demonstrate that our framework can reproduce the observed quasar proximity zone sizes, and that compact proximity zones around z=6 quasars can be explained by the late onset of super-Eddington accretion., Comment: Resubmitted to ApJ, after addressing referee report

Precise measurements of black hole masses are essential to understanding the coevolution of these sources and their host galaxies. We develop a novel approach for computing black hole virial masses using measurements of continuum luminosities and emission line widths from partially overlapping, narrow-band observations of quasars; we refer to this technique as single-epoch photometry. This novel method relies on forward-modelling quasar observations for estimating emission line widths, which enables unbiased measurements even for lines coarsely resolved by narrow-band data. We assess the performance of this technique using quasars from the Sloan Digital Sky Survey (SDSS) observed by the miniJPAS survey, a proof-of-concept project of the Javalambre Physics of the Accelerating Universe Astrophysical Survey (J-PAS) collaboration covering $\simeq1\,\mathrm{deg}^2$ of the northern sky using the 56 J-PAS narrow-band filters. We find remarkable agreement between black hole masses from single-epoch SDSS spectra and single-epoch miniJPAS photometry, with no systematic difference between these and a scatter ranging from 0.4 to 0.07 dex for masses from $\log(M_\mathrm{BH})\simeq8$ to 9.75, respectively. Reverberation mapping studies show that single-epoch masses present approximately 0.4 dex precision, letting us conclude that our novel technique delivers black hole masses with only mildly lower precision than single-epoch spectroscopy. The J-PAS survey will soon start observing thousands of square degrees without any source preselection other than the photometric depth in the detection band, and thus single-epoch photometry has the potential to provide details on the physical properties of quasar populations that do not satisfy the preselection criteria of previous spectroscopic surveys., Comment: 16 pages, 10 figures, accepted by A&A

Recent claimed detections of tidal disruption events (TDEs) in multi-wavelength data have opened potential new windows into the evolution and properties of otherwise dormant supermassive black holes (SMBHs) in the centres of galaxies. At present, there are several dozen TDE candidates, which share some properties and differ in others. The range in properties is broad enough to overlap other transient types, such as active galactic nuclei (AGN) and supernovae (SNe), which can make TDE classification ambiguous. A further complication is that "TDE signatures" have not been uniformly observed to similar sensitivities or even targeted across all candidates. This chapter reviews those events that are unusual relative to other TDEs, including the possibility of TDEs in pre-existing AGN, and summarises those characteristics thought to best distinguish TDEs from continuously accreting AGN, strongly flaring AGN, SNe, and Gamma-Ray Bursts (GRBs), as well as other potential impostors like stellar collisions, "micro-TDEs," and circumbinary accretion flows. We conclude that multiple observables should be used to classify any one event as a TDE. We also consider the TDE candidate population as a whole, which, for certain host galaxy or SMBH characteristics, is distinguishable statistically from non-TDEs, suggesting that at least some TDE candidates do in fact arise from SMBH-disrupted stars., Comment: 57 pages, 17 figures. Accepted for publication in Springer Space Science Reviews. Chapter in ISSI review "The Tidal Disruption of Stars by Massive Black Holes" vol. 79. Corrected several typos from previous submission

In this study, we use the SWIFT/BAT AGN sample, which has received extensive multiwavelength follow-up analysis as a result of the BAT AGN Spectroscopic Survey (BASS), to develop a diagnostic for nuclear obscuration by examining the relationship between the line-of-sight column densities ($N_{\rm{H}}$), the 2-10 keV-to-$12\,\rm{\mu m}$ luminosity ratio, and WISE mid-infrared colors. We demonstrate that heavily obscured AGNs tend to exhibit both preferentially ''redder'' mid-infrared colors and lower values of $L_{\rm{X,\,Obs.}}$/$L_{12\,\rm{\mu m}}$ than less obscured AGNs, and we derive expressions relating $N_{\rm{H}}$ to the $L_{\rm{X,\,Obs.}}$/$L_{12\,\rm{\mu m}}$ and $L_{22\,\rm{\mu m}}$/$L_{4.6\,\rm{\mu m}}$ luminosity ratios as well as develop diagnostic criteria using these ratios. Our diagnostic regions yield samples that are $\gtrsim80$% complete and $\gtrsim60$% pure for AGNs with log($N_{\rm{H}})\geq24$, as well as $\gtrsim85$% pure for AGNs with $\rm{log}(N_{\rm{H}})\gtrsim23.5$. We find that these diagnostics cannot be used to differentiate between optically star forming galaxies and active galaxies. Further, mid-IR contributions from host galaxies that dominate the observed $12~\rm{\mu m}$ emission can lead to larger apparent X-ray deficits and redder mid-IR colors than the AGNs would intrinsically exhibit, though this effect helps to better separate less obscured and more obscured AGNs. Finally, we test our diagnostics on two catalogs of AGNs and infrared galaxies, including the XMM-Newton XXL-N field, and we identify several known Compton-thick AGNs as well as a handful of candidate heavily obscured AGNs based upon our proposed obscuration diagnostics., Comment: 28 pages: 1-19 main text, 20-28 appendix and references. 15 figures. Accepted for publication in ApJS; part of BASS DR2 special issue

We present the host galaxy molecular gas properties of a sample of 213 nearby (0.0110^44 erg/s) increases by ~10-100 between a molecular gas mass of 10^8.7 Msun and 10^10.2 Msun. Higher Eddington ratio AGN galaxies tend to have higher molecular gas masses and gas fractions. Higher column density AGN galaxies (Log NH>23.4) are associated with lower depletion timescales and may prefer hosts with more gas centrally concentrated in the bulge that may be more prone to quenching than galaxy wide molecular gas. The significant average link of host galaxy molecular gas supply to SMBH growth may naturally lead to the general correlations found between SMBHs and their host galaxies, such as the correlations between SMBH mass and bulge properties and the redshift evolution of star formation and SMBH growth., Comment: 53 pages, 37 figures, accepted in ApJ

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., Comment: 26 pages, 10 figures, 6 tables. Submitted after Referee acceptance

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., Comment: 13 pages of text, 6 figures. Invited review talk at IAU Symposium 356, "Nuclear Activity in Galaxies across Cosmic Time", Oct 2019

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., Comment: 22 pages, 10 figures, 2 tables. Accepted for publication in MNRAS

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., Comment: 17 pages, including 8 figures and 4 tables. Accepted for publication in ApJ

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., Comment: 16 pages, 7 figures, accepted for publication in ApJ

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., Comment: 24 pages, 13 figures, accepted by MNRAS 07/01/2020

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., Comment: Accepter for publication in ApJ. Main figure Fig. 14. Simulations code available at https://github.com/nevencaplar/AGN-Variability-Simulations

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., Comment: Accepted for publication in MNRAS (ID: stz2309)

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., Comment: 11+2 pages, 6 figures, Accepted for publication in MNRAS, Comments are welcome

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., Comment: accepted for publication in ApJ

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., Comment: Science White Paper for the US Astro 2020 Decadal Survey

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., Comment: Submitted to the Astro2020 Decadal Survey call for science white papers

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., Comment: Published in Nature Astronomy

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., Comment: To appear in the 8 November 2018 issue of Nature. This is the authors' version of the work

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., Comment: 21 pages, 15 figures, accepted for publication in ApJ. The full list of Table 1 is available at http://www.kusastro.kyoto-u.ac.jp/~ichikawa/Table1_MR_20181107.txt

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., Comment: 28 pages, 14 figures, accepted by ApJ

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., Comment: 5 pages, 2 figures, letter accepted for publication in MNRAS

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., Comment: 11 pages, 3 figues, accepted for publication in MNRAS

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., Comment: Accepted for publication in ApJS. The Swift-BAT 105-month Survey public website can be found at this URL: https://swift.gsfc.nasa.gov/results/bs105mon/

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., Comment: 5 pages of text, 4 figures. Proceeding of the meeting "Quasars at All Cosmic Epochs" (Padova, April 2017). Published in "Frontiers in Astronomy and Space Science" (available at http://www.frontiersin.org/articles/10.3389/fspas.2017.00049)

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., Comment: 14 pages, 10 figures, 1 Appendix, submitted to ApJ

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., Comment: To appear in the 28 September 2017 issue of Nature. This is the authors' version of the work

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., Comment: ApJS in press. The full tables will be available on ApJS and the on the BASS website [www.bass-survey.com] once the paper is published

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., Comment: 22 pages, 14 figures, 9 tables. Accepted for publication in ApJ

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., Comment: Accepted ApJ, see www.bass-survey.com for data

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., Comment: 33 pages, 15 figures, accepted for publication in ApJ, Fig. 6 shows the main result

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., Comment: 12 pages, 6 figures, Accepted for publication in ApJ

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., Comment: Resubmitted to MNRAS after addressing referee's comments

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., Comment: 38 pages, 19 figures. Accepted for publication in MNRAS

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., Comment: This version matches the one published in ApJ. 21 pages, 9 figures

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., Comment: 24 pages, 11 figures, accepted for publication in ApJ. The full list of Table 1 is available at http://www.kusastro.kyoto-u.ac.jp/~ichikawa/table1_MRT.txt

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., Comment: 21 pages, 20 figures, accepted for publication in ApJ

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)., Comment: This version matches the one published in ApJL. 5 text pages, 4 figures

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., Comment: Accepted for publication in MNRAS, 9 pages, 5 figures, 1 table

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

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., Comment: 17 pages, 10 figures, accepted for publication in MNRAS

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., Comment: 26 Pages, 15 figures, Accepted for publication in MNRAS

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., Comment: ApJ, accepted. 16 pages and 6 figures

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., Comment: 53 pages, 10 diagrams, accepted by ApJ

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., Comment: 24 pages, 14 figures, accepted for publication in MNRAS