It is well known that reverberation mapping of active galactic nuclei (AGN) reveals a relationship between AGN luminosity and the size of the broad-line region, and that use of this relationship, combined with the Doppler width of the broad emission line, enables an estimate of the mass of the black hole at the center of the active nucleus based on a single spectrum. This has been discussed in numerous papers over the last two decades. An unresolved key issue is the choice of parameter used to characterize the line width; generally, most researchers use FWHM in favor of line dispersion (the square root of the second moment of the line profile) because the former is easier to measure, less sensitive to blending with other features, and usually can be measured with greater precision. However, use of FWHM introduces a bias, stretching the mass scale such that high masses are overestimated and low masses are underestimated. Here we investigate estimation of black hole masses in AGN based on individual or "single epoch" observations, with a particular emphasis in comparing mass estimates based on line dispersion and FWHM. We confirm the recent findings that, in addition to luminosity and line width, a third parameter is required to obtain accurate masses and that parameter seems to be Eddington ratio. We present simplified empirical formulae for estimating black hole masses from the Hbeta (4861 A) and C IV (1549 A) emission lines.

In 2014 the NGC 5548 Space Telescope and Optical Reverberation Mapping campaign discovered a two-month anomaly when variations in the absorption and emission lines decorrelated from continuum variations. During this time the soft X-ray part of the intrinsic spectrum had been strongly absorbed by a line-of-sight (LOS) obscurer, which was interpreted as the upper part of a disk wind. Our first paper showed that changes in the LOS obscurer produce the decorrelation between the absorption lines and the continuum. A second study showed that the base of the wind shields the BLR, leading to the emission-line decorrelation. In that study, we proposed the wind is normally transparent with no effect on the spectrum. Changes in the wind properties alter its shielding and affect the SED striking the BLR, producing the observed decorrelations. In this work, we investigate the impact of a translucent wind on the emission lines. We simulate the obscuration using XMM-Newton, NuSTAR, and HST observations to determine the physical characteristics of the wind. We find that a translucent wind can contribute a part of the He II and Fe K? emission. It has a modest optical depth to electron scattering, which explains the fainter far-side emission in the observed velocity delay maps. The wind produces the very broad base seen in the UV emission lines and may also be present in the Fe K? line. Our results highlight the importance of accounting for the effects of such winds in the analysis of the physics of the central engine.

We report velocity-delay maps for prominent broad emission lines, Ly_alpha, CIV, HeII and H_beta, in the spectrum of NGC5548. The emission-line responses inhabit the interior of a virial envelope. The velocity-delay maps reveal stratified ionization structure. The HeII response inside 5-10 light-days has a broad single-peaked velocity profile. The Ly_alpha, CIV, and H_beta responses peak inside 10 light-days, extend outside 20 light-days, and exhibit a velocity profile with two peaks separated by 5000 km/s in the 10 to 20 light-day delay range. The velocity-delay maps show that the M-shaped lag vs velocity structure found in previous cross-correlation analysis is the signature of a Keplerian disk with a well-defined outer edge at R=20 light-days. The outer wings of the M arise from the virial envelope, and the U-shaped interior of the M is the lower half of an ellipse in the velocity-delay plane. The far-side response is weaker than that from the near side, so that we see clearly the lower half, but only faintly the upper half, of the velocity--delay ellipse. The delay tau=(R/c)(1-sin(i))=5 light-days at line center is from the near edge of the inclined ring, giving the inclination i=45 deg. A black hole mass of M=7x10^7 Msun is consistent with the velocity-delay structure. A barber-pole pattern with stripes moving from red to blue across the CIV and possibly Ly_alpha line profiles suggests the presence of azimuthal structure rotating around the far side of the broad-line region and may be the signature of precession or orbital motion of structures in the inner disk. Further HST observations of NGC 5548 over a multi-year timespan but with a cadence of perhaps 10 days rather than 1 day could help to clarify the nature of this new AGN phenomenon.

The structure of the broad line region (BLR) is an essential ingredient in the determination of active galactic nuclei (AGN) virial black hole masses, which in turn are important to study the role of black holes in galaxy evolution. Constraints on the BLR geometry and dynamics can be obtained from velocity-resolved studies using reverberation mapping data (i.e. monitoring data). However, monitoring data are observationally expensive and only available for a limited sample of AGN, mostly confined to the local Universe. Here we explore a new version of a Bayesian inference, physical model of the BLR which uses an individual spectrum and prior information on the BLR size from the radius-luminosity relation, to model the AGN BLR geometry and dynamics. We apply our model to a sample of 11 AGN, which have been previously modelled using monitoring data. Our single-epoch BLR model is able to constrain some of the BLR parameters with inferred parameter values that agree within the uncertainties with those determined from the modelling of monitoring data. We find that our model is able to derive stronger constraints on the BLR for AGN with broad emission lines that qualitatively have more substructure and more asymmetry, presumably as they contain more information to constrain the physical model. The performance of this model makes it a practical and cost-effective tool to determine some of the BLR properties of a large sample of low and high redshift AGN, for which monitoring data are not available.

We show that individual (single-epoch) spectra of AGN can constrain some of the geometry and dynamics of the AGN broad line region. Studies of the cosmic influence of supermassive black holes are limited by the current large uncertainties in the determination of black hole masses. One dominant limitation is the unknown geometry, dynamics and line-of-sight inclination of the broad line region, used to probe the central black hole mass. Recent progress has been made to constrain the spatial and kinematic structure of the broad line region using dynamical modelling of AGN monitoring data and an underlying physical model for the broad line region. In this work we test the ability of a modified version of this dynamical modelling code to constrain the broad line region structure using single-epoch spectra. We test our modelling code on single-epoch spectra of nearby Arp 151 by comparing our results with those obtained with monitoring data of this same object. We find that a significant fraction of the broad line region parameters can indeed be adequately constrained, with uncertainties that are comparable to, or at most a factor of ~ a few higher than those obtained from modelling of monitoring data. Considering the wealth of available single-epoch spectroscopic observations, this method is promising for establishing the overall AGN population trends in the geometry and dynamics of the broad line region. This method can be applied to spectra of AGN at low and high redshift making it valuable for studies of cosmological black hole and AGN evolution.

We model the ultraviolet spectra of the Seyfert 1 galaxy NGC~5548 obtained with the Hubble Space Telescope during the 6-month reverberation-mapping campaign in 2014. Our model of the emission from NGC 5548 corrects for overlying absorption and deblends the individual emission lines. Using the modeled spectra, we measure the response to continuum variations for the deblended and absorption-corrected individual broad emission lines, the velocity-dependent profiles of Ly$\alpha$ and C IV, and the narrow and broad intrinsic absorption features. We find that the time lags for the corrected emission lines are comparable to those for the original data. The velocity-binned lag profiles of Ly$\alpha$ and C IV have a double-peaked structure indicative of a truncated Keplerian disk. The narrow absorption lines show delayed response to continuum variations corresponding to recombination in gas with a density of $\sim 10^5~\rm cm^{-3}$. The high-ionization narrow absorption lines decorrelate from continuum variations during the same period as the broad emission lines. Analyzing the response of these absorption lines during this period shows that the ionizing flux is diminished in strength relative to the far-ultraviolet continuum. The broad absorption lines associated with the X-ray obscurer decrease in strength during this same time interval. The appearance of X-ray obscuration in $\sim\,2012$ corresponds with an increase in the luminosity of NGC 5548 following an extended low state. We suggest that the obscurer is a disk wind triggered by the brightening of NGC 5548 following the decrease in size of the broad-line region during the preceding low-luminosity state.

We present ALMA observations of the CO(2-1) and CO(3-2) molecular gas transitions and associated (sub)-mm continua of the nearby Seyfert 1.5 galaxy NGC3227 with angular resolutions 0.085-0.21" (7-15pc). On large scales the cold molecular gas shows circular motions as well as streaming motions on scales of a few hundred parsecs associated with a large scale bar. We fitted the nuclear ALMA 1.3mm emission with an unresolved component and an extended component. The 850$\mu$m emission shows at least two extended components, one along the major axis of the nuclear disk and the other along the axis of the ionization cone. The molecular gas in the central region (1" ~73pc) shows several CO clumps with complex kinematics which appears to be dominated by non-circular motions. While we cannot demonstrate conclusively the presence of a warped nuclear disk, we also detected non-circular motions along the kinematic minor axis. They reach line-of-sight velocities of v-vsys =150-200km/s. Assuming that the radial motions are in the plane of the galaxy, then we interpret them as a nuclear molecular outflow due to molecular gas in the host galaxy being entrained by the AGN wind. We derive molecular outflow rates of $5\,M_\odot\,{\rm yr}^{-1}$ and $0.6\,M_\odot\,{\rm yr}^{-1}$ at projected distances of up to 30pc to the northeast and southwest of the AGN, respectively. At the AGN location we estimate a mass in molecular gas of $5\times 10^{5}\,M_\odot$ and an average column density $N({\rm H}_2) = 2-3\times 10^{23}\,{\rm cm}^{-2}$ in the inner 15pc. The nuclear molecular gas and sub-mm continuum emission of NGC3227 do not resemble the classical compact torus. Rather, these emissions extend for several tens of parsecs and appear connected with the circumnuclear ring in the host galaxy disk, as found in other local AGN. (Abridged)

The flux variations in the emission lines in active galactic nuclei (AGNs) are driven by variations in the ionizing continuum flux --which are usually reflected in the observable UV-optical continuum. The "Reverberation mapping" technique measures the delay between line and continuum variations to determine the size of the line emitting region, this is the basis for measurements of the central black hole mass in AGNs. The Space Telescope and Optical Reverberation Mapping Project (AGN STORM) on NGC 5548 in 2014 is the most intensive multi-wavelength AGN monitoring campaign ever. For most of the campaign, the emission-line variations followed changes in the continuum with a time lag, as expected. However, the lines varied independently of the observed UV-optical continuum during a 60 -- 70 day "holiday." To understand this remarkable phenomenon, we study the intrinsic absorption lines present in NGC 5548. We identify a novel cycle that reproduces the absorption line variability and thus identify the physics that allows the holiday to occur. In our model, variations in this obscurer's line-of-sight covering factor modify the soft X-ray continuum. This leads to changes in the ionization of helium gas in the broad-line region. Ionizing radiation produced by recombining helium then affects the ionization of other species as observed during the AGN STORM holiday. It is likely that any other model which selectively changes the soft X-ray part of the continuum during the holiday can also explain the anomalous emission line behavior observed.

Optical changing-look Active Galactic Nuclei (AGN) are a class of sources that change type within a short timescale of years or decades. This change is characterised by the appearance or disappearance of broad emission lines, often associated with dramatic AGN continuum flux changes that are orders of magnitude larger than those expected from typical AGN variability. In this work we study for the first time the host galaxy of a changing-look AGN, Mrk 590, using high spatial resolution optical and near-infrared observations. We discover that after ~ 10 yr absence, the optical broad emission lines of Mrk 590 have reappeared. The AGN optical continuum flux however, is still ~ 10 times lower than that observed during the most luminous state in the 1990s. The host galaxy shows a 4.5 kpc radius star-forming ring with knots of ionised and cold molecular gas emission. Extended ionised and warm molecular gas emission are detected in the nucleus, indicating that there is a reservoir of gas as close as 60 pc from the black hole. We observe a nuclear gas spiral between radii r ~ 0.5 - 2 kpc, which has been suggested as a dynamical mechanism able to drive the necessary gas to fuel AGN. We also discover blue-shifted and high velocity dispersion [O III] emission out to a radius of 1 kpc, tracing a nuclear gas outflow. The gas dynamics in Mrk 590 suggest a complex balance between gas inflow and outflow in the nucleus of the galaxy.

The discovery of luminous quasars at redshifts up to 7.5 demonstrates the existence of several billion M_sun supermassive black holes (SMBHs) less than a billion years after the Big Bang. They are accompanied by intense star formation in their host galaxies, pinpointing sites of massive galaxy assembly in the early universe, while their absorption spectra reveal an increasing neutral intergalactic medium (IGM) at the epoch of reionization. Extrapolating from the rapid evolution of the quasar density at z=5-7, we expect that there is only one luminous quasar powered by a billion M_sun SMBH in the entire observable universe at z~9. In the next decade, new wide-field, deep near-infrared (NIR) sky surveys will push the redshift frontier to the first luminous quasars at z~9-10; the combination with new deep X-ray surveys will probe fainter quasar populations that trace earlier phases of SMBH growth. The identification of these record-breaking quasars, and the measurements of their BH masses and accretion properties require sensitive spectroscopic observations with next generation of ground-based and space telescopes at NIR wavelengths. High-resolution integral-field spectroscopy at NIR, and observations at millimeter and radio wavelengths, will together provide a panchromatic view of the quasar host galaxies and their galactic environment at cosmic dawn, connecting SMBH growth with the rise of the earliest massive galaxies. Systematic surveys and multiwavelength follow-up observations of the earliest luminous quasars will strongly constrain the seeding and growth of the first SMBHs in the universe, and provide the best lines of sight to study the history of reionization.

We report initial results from a large Gemini program to observe z>~5.7 quasars with GNIRS near-IR spectroscopy. Our sample includes 50 quasars with simultaneous ~0.85-2.5 micron spectra covering the rest-frame ultraviolet and major broad emission lines from Ly-alpha to MgII. We present spectral measurements for these quasars and compare to their lower-redshift counterparts at z=1.5-2.3. We find that when quasar luminosity is matched, there are no significant differences between the rest-UV spectra of z>~5.7 quasars and the low-z comparison sample. High-z quasars have similar continuum and emission line properties and occupy the same region in the black hole mass and luminosity space as the comparison sample, accreting at an average Eddington ratio of ~0.3. There is no evidence for super-Eddington accretion or hypermassive (>10^10 Msun) black holes within our sample. We find a mild excess of quasars with weak CIV lines relative to the control sample. Our results, corroborating earlier studies but with better statistics, demonstrate that these high-z quasars are already mature systems of accreting supermassive black holes operating with the same physical mechanisms as those at lower redshifts.

Swift intensive accretion disk reverberation mapping of four AGN yielded light curves sampled $\sim$200-350 times in 0.3-10 keV X-ray and six UV/optical bands. Uniform reduction and cross-correlation analysis of these datasets yields three main results: 1) The X-ray/UV correlations are much weaker than those within the UV/optical, posing severe problems for the lamp-post reprocessing model in which variations in a central X-ray corona drive and power those in the surrounding accretion disk. 2) The UV/optical interband lags are generally consistent with $ \tau \propto \lambda^{4/3} $ as predicted by the centrally illuminated thin accretion disk model. While the average interband lags are somewhat larger than predicted, these results alone are not inconsistent with the thin disk model given the large systematic uncertainties involved. 3) The one exception is the U band lags, which are on average a factor of $\sim$2.2 larger than predicted from the surrounding band data and fits. This excess appears due to diffuse continuum emission from the broad-line region (BLR). The precise mixing of disk and BLR components cannot be determined from these data alone. The lags in different AGN appear to scale with mass or luminosity. We also find that there are systematic differences between the uncertainties derived by javelin vs. more standard lag measurement techniques, with javelin reporting smaller uncertainties by a factor of 2.5 on average. In order to be conservative only standard techniques were used in the analyses reported herein.

Mrk 590 was originally classified as a Seyfert 1 galaxy, but then it underwent dramatic changes: the nuclear luminosity dropped by over two orders of magnitude and the broad emission lines all but disappeared from the optical spectrum. Here we present followup observations to the original discovery and characterization of this "changing look" active galactic nucleus (AGN). The new Chandra and HST observations from 2014 show that Mrk 590 is awakening, changing its appearance again. While the source continues to be in a low state, its soft excess has re-emerged, though not to the previous level. The UV continuum is brighter by more than a factor of two and the broad MgII emission line is present, indicating that the ionizing continuum is also brightening. These observations suggest that the soft excess is not due to reprocessed hard X-ray emission. Instead, it is connected to the UV continuum through warm Comptonization. Variability of the Fe K-alpha emission lines suggests that the reprocessing region is within about 10 light years or 3 pc of the central source. The AGN type change is neither due to obscuration, nor due to one-way evolution from type-1 to type-2, as suggested in literature, but may be related to episodic accretion events.

Disk reverberation mapping of a handful of nearby AGN suggest accretion disk sizes which are a factor few too large for their luminosities, apparently at odds with the standard model. Here, we investigate the likely contribution to the measured delay signature of diffuse continuum emission arising from broad line region gas. We start by constructing spherically symmetric pressure-law BLR models (i.e., $P(r)\propto r^{-s}$) that approximately reproduce the observed emission line fluxes of the strong UV--optical emission-lines in the best-studied source, NGC~5548. We then determine the contribution of the diffuse continuum to the measured continuum flux and inter-band delays, accounting for the observed variability behaviour of the ionizing nuclear continuum. Those pressure-law models that approximately reproduce the observed emission-line luminosities unavoidably produce substantial diffuse continuum emission. This causes a significant contamination of the disk reverberation signature (i.e., wavelength-dependent continuum delays). Qualitatively, the diffuse continuum delay signatures produced by our models resemble that observed for NGC~5548, including the deviation of the lag spectrum above that of a simple power-law in wavelength, short-ward of the Balmer and Paschen jumps. Furthermore, for reasonable estimates of the BLR covering fraction, the delay induced by diffuse continuum emission causes elevated inter-band delays over the entire UV--optical regime; for these pressure-law models, there are no `disk-dominated' wavelength intervals. Thus, the diffuse continuum contribution must be taken into account in order to correctly infer AGN accretion disk sizes based on inter-band continuum delays.

We present the first results from a reverberation-mapping campaign undertaken during the first half of 2012, with additional data on one AGN (NGC 3227) from a 2014 campaign. Our main goals are (1) to determine the black hole masses from continuum-Hbeta reverberation signatures, and (2) to look for velocity-dependent time delays that might be indicators of the gross kinematics of the broad-line region. We successfully measure Hbeta time delays and black hole masses for five AGNs, four of which have previous reverberation mass measurements. The values measured here are in agreement with earlier estimates, though there is some intrinsic scatter beyond the formal measurement errors. We observe velocity dependent Hbeta lags in each case, and find that the patterns have changed in the intervening five years for three AGNs that were also observed in 2007.

We present a Hubble Space Telescope (HST) survey of extended [O III] {\lambda}5007 emission for a sample of 12 nearby (z < 0.12), luminous Type 2 quasars (QSO2s), which we use to measure the extent and kinematics of their AGN-ionized gas. We find the size of the observed [O III] regions scale with luminosity in comparison to nearby, less luminous Seyfert galaxies and radially outflowing kinematics to exist in all targets. We report an average maximum outflow radius of $\sim$600 pc, with gas continuing to be kinematically influenced by the central AGN out to an average radius of $\sim$1130 pc. These findings question the effectiveness of AGN being capable of clearing material from their host bulge in the nearby universe and suggest that disruption of gas by AGN activity may prevent star formation without requiring evacuation. Additionally, we find a dichotomy in our targets when comparing [O III] radial extent and nuclear FWHM, where QSO2s with compact [O III] morphologies typically possess broader nuclear emission-lines.

We report the results of intensive X-ray, UV and optical monitoring of the Seyfert 1 galaxy NGC 4593 with Swift. There is no intrinsic flux-related spectral change in the the variable components in any band with small apparent variations due only to contamination by a second constant component, possibly a (hard) reflection component in the X-rays and the (red) host galaxy in the UV/optical bands. Relative to the shortest wavelength band, UVW2, the lags of the other UV and optical bands are mostly in agreement with the predictions of reprocessing of high energy emission from an accretion disc. The U-band lag is, however, far larger than expected, almost certainly because of reprocessed Balmer continuum emission from the more distant broad line region gas. The UVW2 band is well correlated with the X-rays but lags by ~6x more than expected if the UVW2 results from reprocessing of X-rays on the accretion disc. However, if the lightcurves are filtered to remove variations on timescales >5d, the lag approaches the expectation from disc reprocessing. MEMEcho analysis shows that direct X-rays can be the driver of most of the variations in the UV/optical bands as long as the response functions for those bands all have long tails (up to 10d) in addition to a strong peak (from disc reprocessing) at short lag (<1d). We interpret the tails as due to reprocessing from the surrounding gas. Comparison of X-ray to UVW2 and UVW2 to V-band lags for 4 AGN, including NGC 4593, shows that all have UVW2 to V-band lags which exceed the expectations from disc resprocessing by factor < 2. However the X-ray to UVW2 lags are, mostly, in greater excess from the expectations from disc reprocessing and differ between AGN. The largest excess is in NGC 4151. Absorption and scattering may be affecting X-ray to UV lags.

We study the host galaxies of four Iron Low-Ionization Broad Absorption-line Quasars (FeLoBALs) using Hubble Space Telescope imaging data, investigating the possibility that they represent a transition between an obscured AGN and an ordinary optical quasar. In this scenario, the FeLoBALs represent the early stage of merger-triggered accretion, in which case their host galaxies are expected to show signs of an ongoing or recent merger. Using PSF subtraction techniques, we decompose the images into host galaxy and AGN components at rest-frame ultraviolet and optical wavelengths. The ultraviolet is sensitive to young stars, while the optical probes stellar mass. In the ultraviolet we image at the BAL absorption trough wavelengths so as to decrease the contrast between the quasar and host galaxy emission. We securely detect an extended source for two of the four FeLoBALs in the rest-frame optical; a third host galaxy is marginally detected. In the rest-frame UV we detect no host emission; this constrains the level of unobscured star formation. Thus, the host galaxies have observed properties that are consistent with those of non-BAL quasars with the same nuclear luminosity, i.e., quiescent or moderately starforming elliptical galaxies. However, we cannot exclude starbursting hosts that have the stellar UV emission obscured by modest amounts of dust reddening. Thus, our findings also allow the merger-induced young quasar scenario. For three objects, we identify possible close companion galaxies that may be gravitationally interacting with the quasar hosts.

We present and analyse quasi-simultaneous multi-epoch spectral energy distributions (SEDs) of seven reverberation-mapped AGNs for which accurate black hole mass measurements and suitable archival data are available from the 'AGN Watch' monitoring programs. We explore the potential of optical-UV and X-ray data, obtained within 2 days, to provide more accurate SED-based measurements of individual AGN and quantify the impact of source variability on key measurements typically used to characterise the black hole accretion process plus on bolometric correction factors at 5100 {\AA}, 1350 {\AA} and for the 2-10 keV X-ray band, respectively. The largest SED changes occur on long timescales (>1 year). For our small sample, the 1 micron to 10 keV integrated accretion luminosity typically changes by 10% on short time-scales (over 20 days), by ~30% over a year, but can change by 100% or more for individual AGN. The EUV gap is the most uncertain part of the intrinsic SED, introducing a ~25% uncertainty in the accretion-induced luminosity, relative to the model independent interpolation method that we adopt. That aside, our analysis shows that the uncertainty in the accretion-induced luminosity, the Eddington luminosity ratio and the bolometric correction factors can be reduced (by a factor of two or more) by use of SEDs built from data obtained within 20 days. However, mass accretion rate and mass accretion efficiency are mostly limited by the unknown EUV emission and the unknown details of the central engine and our aspect angle.

Previous studies have shown that the incidence rate of intervening strong MgII absorbers towards GRBs were a factor of 2 - 4 higher than towards quasars. Exploring the similar sized and uniformly selected legacy data sets XQ-100 and XSGRB, each consisting of 100 quasar and 81 GRB afterglow spectra obtained with a single instrument (VLT/X-shooter), we demonstrate that there is no disagreement in the number density of strong MgII absorbers with rest-frame equivalent widths $W_r^{2796} >$ 1 {\AA} towards GRBs and quasars in the redshift range 0.1 < z < 5. With large and similar sample sizes, and path length coverages of $\Delta$z = 57.8 and 254.4 for GRBs and quasars, respectively, the incidences of intervening absorbers are consistent within 1 sigma uncertainty levels at all redshifts. For absorbers at z < 2.3 the incidence towards GRBs is a factor of 1.5$\pm$0.4 higher than the expected number of strong MgII absorbers in SDSS quasar spectra, while for quasar absorbers observed with X-shooter we find an excess factor of 1.4$\pm$0.2 relative to SDSS quasars. Conversely, the incidence rates agree at all redshifts with reported high spectral resolution quasar data, and no excess is found. The only remaining discrepancy in incidences is between SDSS MgII catalogues and high spectral resolution studies. The rest-frame equivalent width distribution also agrees to within 1 sigma uncertainty levels between the GRB and quasar samples. Intervening strong MgII absorbers towards GRBs are therefore neither unusually frequent, nor unusually strong.

During the Space Telescope and Optical Reverberation Mapping Project (STORM) observations of NGC 5548, the continuum and emission-line variability became de-correlated during the second half of the 6-month long observing campaign. Here we present Swift and Chandra X-ray spectra of NGC 5548 obtained as a part of the campaign. The Swift spectra show that excess flux (relative to a power-law continuum) in the soft X-ray band appears before the start of the anomalous emission-line behavior, peaks during the period of the anomaly, and then declines. This is a model-independent result suggesting that the soft excess is related to the anomaly. We divide the Swift data into on- and off-anomaly spectra to characterize the soft excess via spectral fitting. The cause of the spectral differences is likely due to a change in the intrinsic spectrum rather than being due to variable obscuration or partial covering. The Chandra spectra have lower signal-to-noise ratios, but are consistent with Swift data. Our preferred model of the soft excess is emission from an optically thick, warm Comptonizing corona, the effective optical depth of which increases during the anomaly. This model simultaneously explains all the three observations: the UV emission line flux decrease, the soft-excess increase, and the emission line anomaly.

Swift monitoring of NGC 4151 with ~6 hr sampling over a total of 69 days in early 2016 is used to construct light curves covering five bands in the X-rays (0.3-50 keV) and six in the ultraviolet (UV)/optical (1900-5500 A). The three hardest X-ray bands (>2.5 keV) are all strongly correlated with no measurable interband lag while the two softer bands show lower variability and weaker correlations. The UV/optical bands are significantly correlated with the X-rays, lagging ~3-4 days behind the hard X-rays. The variability within the UV/optical bands is also strongly correlated, with the UV appearing to lead the optical by ~0.5-1 day. This combination of >~3 day lags between the X-rays and UV and <~1 day lags within the UV/optical appears to rule out the "lamp-post" reprocessing model in which a hot, X-ray emitting corona directly illuminates the accretion disk, which then reprocesses the energy in the UV/optical. Instead, these results appear consistent with the Gardner & Done picture in which two separate reprocessings occur: first, emission from the corona illuminates an extreme-UV-emitting toroidal component that shields the disk from the corona; this then heats the extreme-UV component which illuminates the disk and drives its variability.

We present a re-calibration of the $M_{BH}-\sigma_{\star}$ relation, based on a sample of 16 reverberation-mapped galaxies with newly determined bulge stellar velocity dispersions ($\sigma_{\star}$) from integral-field spectroscopy (IFS), and a sample of 32 quiescent galaxies with publicly available IFS. For both samples, $\sigma_{\star}$ is determined via two different methods that are popular in the literature, and we provide fits for each sample based on both sets of $\sigma_{\star}$. We find the fit to the AGN sample is shallower than the fit to the quiescent galaxy sample, and that the slopes for each sample are in agreement with previous investigations. However, the intercepts to the quiescent galaxy relations are notably higher than those found in previous studies, due to the systematically lower $\sigma_{\star}$ measurements that we obtain from IFS. We find that this may be driven, in part, by poorly constrained measurements of bulge effective radius ($r_{e}$) for the quiescent galaxy sample, which may bias the $\sigma_{\star}$ measurements low. We use these quiescent galaxy parameterizations, as well as one from the literature, to recalculate the virial scaling factor $f$. We assess the potential biases in each measurement, and suggest $f=4.82\pm1.67$ as the best currently available estimate. However, we caution that the details of how $\sigma_{\star}$ is measured can significantly affect $f$, and there is still much room for improvement.

We present the Swift optical to X-ray spectral energy distributions (SEDs) of 44 quasars at redshifts z~2 observed by Swift, part of a larger program to establish and characterize the optical through X-ray SEDs of moderate-redshift quasars. Here we outline our analysis approach and present preliminary analysis and results for the first third of the full quasar sample. Not all quasars in the sample are detected in X-rays; all of the X-ray detected objects so far are radio loud. As expected for radio loud objects, they are X-ray bright relative to radio-quiet quasars of comparable optical luminosities, with an average alpha_ox = 1.39 +/- 0.03 (where alpha_ox is the power-law slope connecting the monochromatic flux at 2500 Ang and at 2 keV), and display hard X-ray spectra. We find integrated 3000 Ang - 25 keV accretion luminosities of between 0.7*10^(46) erg s^(-1) and 5.2*10^(47) erg s^(-1). Based on single-epoch spectroscopic virial black hole mass estimates, we find that these quasars are accreting at substantial Eddington fractions, 0.1 \le L/LEdd \le 1.

We present the results of an optical spectroscopic monitoring program targeting NGC 5548 as part of a larger multi-wavelength reverberation mapping campaign. The campaign spanned six months and achieved an almost daily cadence with observations from five ground-based telescopes. The H$\beta$ and He II $\lambda$4686 broad emission-line light curves lag that of the 5100 $\AA$ optical continuum by $4.17^{+0.36}_{-0.36}$ days and $0.79^{+0.35}_{-0.34}$ days, respectively. The H$\beta$ lag relative to the 1158 $\AA$ ultraviolet continuum light curve measured by the Hubble Space Telescope is roughly $\sim$50% longer than that measured against the optical continuum, and the lag difference is consistent with the observed lag between the optical and ultraviolet continua. This suggests that the characteristic radius of the broad-line region is $\sim$50% larger than the value inferred from optical data alone. We also measured velocity-resolved emission-line lags for H$\beta$ and found a complex velocity-lag structure with shorter lags in the line wings, indicative of a broad-line region dominated by Keplerian motion. The responses of both the H$\beta$ and He II $\lambda$4686 emission lines to the driving continuum changed significantly halfway through the campaign, a phenomenon also observed for C IV, Ly $\alpha$, He II(+O III]), and Si IV(+O IV]) during the same monitoring period. Finally, given the optical luminosity of NGC 5548 during our campaign, the measured H$\beta$ lag is a factor of five shorter than the expected value implied by the $R_\mathrm{BLR} - L_\mathrm{AGN}$ relation based on the past behavior of NGC 5548.

We conduct a multiwavelength continuum variability study of the Seyfert 1 galaxy NGC 5548 to investigate the temperature structure of its accretion disk. The 19 overlapping continuum light curves (1158 to 9157 angstroms) combine simultaneous HST , Swift , and ground-based observations over a 180 day period from 2014 January to July. Light-curve variability is interpreted as the reverberation response of the accretion disk to irradiation by a central time-varying point source. Our model yields the disk inclination, i, temperature T1 at 1 light day from the black hole, and a temperature-radius slope, alpha. We also infer the driving light curve and find that it correlates poorly with both the hard and soft X-ray light curves, suggesting that the X-rays alone may not drive the ultraviolet and optical variability over the observing period. We also decompose the light curves into bright, faint, and mean accretion-disk spectra. These spectra lie below that expected for a standard blackbody accretion disk accreting at L/LEdd = 0.1

The first generation of E-ELT instruments will include an optical-infrared High Resolution Spectrograph, conventionally indicated as EELT-HIRES, which will be capable of providing unique breakthroughs in the fields of exoplanets, star and planet formation, physics and evolution of stars and galaxies, cosmology and fundamental physics. A 2-year long phase A study for EELT-HIRES has just started and will be performed by a consortium composed of institutes and organisations from Brazil, Chile, Denmark, France, Germany, Italy, Poland, Portugal, Spain, Sweden, Switzerland and United Kingdom. In this paper we describe the science goals and the preliminary technical concept for EELT-HIRES which will be developed during the phase A, as well as its planned development and consortium organisation during the study.

We present the results of a new spectroscopic survey for dusty intervening absorption systems, particularly damped Ly$\alpha$ absorbers (DLAs), towards reddened quasars. The candidate quasars are selected from mid-infrared photometry from the Wide-field Infrared Survey Explorer combined with optical and near-infrared photometry. Out of 1073 candidates, we secure low-resolution spectra for 108 using the Nordic Optical Telescope on La Palma, Spain. Based on the spectra, we are able to classify 100 of the 108 targets as quasars. A large fraction (50 %) is observed to have broad absorption lines (BALs). Moreover, we find 6 quasars with strange breaks in their spectra, which are not consistent with regular dust reddening. Using template fitting we infer the amount of reddening along each line of sight ranging from A(V)$\approx$0.1 mag to 1.2 mag (assuming an SMC extinction curve). In four cases, the reddening is consistent with dust exhibiting the 2175{\AA} feature caused by an intervening absorber, and for two of these, a MgII absorption system is observed at the best-fit absorption redshift. In the rest of the cases, the reddening is most likely intrinsic to the quasar. We observe no evidence for dusty DLAs in this survey. However, the large fraction of BAL quasars hampers the detection of absorption systems. Out of the 50 non-BAL quasars only 28 have sufficiently high redshift to detect Ly$\alpha$ in absorption.

The sub-population of quasars reddened by intrinsic or intervening clouds of dust are known to be underrepresented in optical quasar surveys. By defining a complete parent sample of the brightest and spatially unresolved quasars in the COSMOS field, we quantify to which extent this sub-population is fundamental to our understanding of the true population of quasars. By using the available multiwavelength data of various surveys in the COSMOS field, we built a parent sample of 33 quasars brighter than $J=20$ mag, identified by reliable X-ray to radio wavelength selection techniques. Spectroscopic follow-up with the NOT/ALFOSC was carried out for four candidate quasars that had not been targeted previously to obtain a 100\% redshift completeness of the sample. The population of high $A_V$ quasars (HAQs), a specific sub-population of quasars selected from optical/near-infrared photometry, is found to contribute $21\%^{+9}_{-5}$ of the parent sample. The full population of bright spatially unresolved quasars represented by our parent sample consists of $39\%^{+9}_{-8}$ reddened quasars defined by having $A_V>0.1$, and $21\%^{+9}_{-5}$ of the sample having $E(B-V)>0.1$ assuming the extinction curve of the Small Magellanic Cloud. We show that the HAQ selection works well for selecting reddened quasars, but some are missed because their optical spectra are too blue to pass the $g-r$ color cut in the HAQ selection. This is either due to a low degree of dust reddening or anomalous spectra. We find that the fraction of quasars with contributing light from the host galaxy is most dominant at $z \lesssim 1$. At higher redshifts the population of spatially unresolved quasars selected by our parent sample is found to be representative of the full population at $J<20$ mag. This work quantifies the bias against reddened quasars in studies that are based solely on optical surveys.

We describe the execution and data reduction of the European Southern Observatory Large Programme "Quasars and their absorption lines: a legacy survey of the high-redshift universe with VLT/XSHOOTER" (hereafter `XQ-100'). XQ-100 has produced and made publicly available an homogeneous and high-quality sample of echelle spectra of 100 QSOs at redshifts z~3.5-4.5 observed with full spectral coverage from 315 to 2500 nm at a resolving power ranging from R~4000 to 7000, depending on wavelength. The median signal-to-noise ratios are 33, 25 and 43, as measured at rest-frame wavelengths 1700, 3000 and 3600 Angstrom, respectively. This paper provides future users of XQ-100 data with the basic statistics of the survey, along with details of target selection, data acquisition and data reduction. The paper accompanies the public release of all data products, including 100 reduced spectra. XQ-100 is the largest spectroscopic survey to date of high-redshift QSOs with simultaneous rest-frame UV/optical coverage, and as such enables a wide range of extragalactic research, from cosmology and galaxy evolution to AGN astrophysics.

We present the serendipitous discovery of a projected pair of quasi-stellar objects (QSOs) with an angular separation of $\Delta\theta =4.50$ arcsec. The redshifts of the two QSOs are widely different: one, our programme target, is a QSO with a spectrum consistent with being a narrow line Seyfert 1 AGN at $z=2.05$. For this target we detect Lyman-$\alpha$, \ion{C}{4}, and \ion{C}{3]}. The other QSO, which by chance was included on the spectroscopic slit, is a Type 1 QSO at a redshift of $z=1.68$, for which we detect \ion{C}{4}, \ion{C}{3]} and \ion{Mg}{2}. We compare this system to previously detected projected QSO pairs and find that only about a dozen previously known pairs have smaller angular separation.

We statistically study the physical properties of a sample of narrow absorption line (NAL) systems looking for empirical evidences to distinguish between intrinsic and intervening NALs without taking into account any a priori definition or velocity cut-off. We analyze the spectra of 100 quasars with 3.5 < z$\rm_{em}$ < 4.5, observed with X-shooter/VLT in the context of the XQ-100 Legacy Survey. We detect a $\sim$ 8 $\sigma$ excess in the number density of absorbers within 10,000 km/s of the quasar emission redshift with respect to the random occurrence of NALs. This excess does not show a dependence on the quasar bolometric luminosity and it is not due to the redshift evolution of NALs. It extends far beyond the standard 5000 km/s cut-off traditionally defined for associated absorption lines. We propose to modify this definition, extending the threshold to 10,000 km/s when also weak absorbers (equivalent width < 0.2 \AA) are considered. We infer NV is the ion that better traces the effects of the quasar ionization field, offering the best statistical tool to identify intrinsic systems. Following this criterion we estimate that the fraction of quasars in our sample hosting an intrinsic NAL system is 33 percent. Lastly, we compare the properties of the material along the quasar line of sight, derived from our sample, with results based on close quasar pairs investigating the transverse direction. We find a deficiency of cool gas (traced by CII) along the line of sight associated with the quasar host galaxy, in contrast with what is observed in the transverse direction.

We investigate the origin of the parsec-scale radio emission from the changing-look active galactic nucleus (AGN) of Mrk 590, and examine whether the radio power has faded concurrently with the dramatic decrease in accretion rates observed between the 1990s and the present. We detect a compact core at 1.6 GHz and 8.4 GHz using new Very Long Baseline Array observations, finding no significant extended, jet-like features down to $\sim$1 pc scales. The flat spectral index ($\alpha_{1.6}^{8.4} = 0.03$) and high brightness temperature ($T_{\rm b} \sim 10^{8}\,\rm K$) indicate self-absorbed synchrotron emission from the AGN. The radio to X-ray luminosity ratio of ${\rm log}(L_{\rm R}/L_{\rm X}) \sim -5$, similar to that in coronally active stars, suggests emission from magnetized coronal winds, although unresolved radio jets are also consistent with the data. Comparing new Karl G. Jansky Very Large Array measurements with archival and published radio flux densities, we find $46\%$, $34\%$, and (insignificantly) $13\%$ flux density decreases between the 1990s and the year 2015 at 1.4 GHz, 5 GHz and 8.4 GHz respectively. This trend, possibly due to the expansion and fading of internal shocks within the radio-emitting outflow after a recent outburst, is consistent with the decline of the optical-UV and X-ray luminosities over the same period. Such correlated variability demonstrates the AGN accretion-outflow connection, confirming that the changing-look behaviour in Mrk 590 originates from variable accretion rates rather than dust obscuration. The present radio and X-ray luminosity correlation, consistent with low/hard state accretion, suggests that the black hole may now be accreting in a radiatively inefficient mode.

During an intensive Hubble Space Telescope (HST) Cosmic Origins Spectrograph (COS) UV monitoring campaign of the Seyfert~1 galaxy NGC 5548 performed from 2014 February to July, the normally highly correlated far-UV continuum and broad emission-line variations decorrelated for ~60 to 70 days, starting ~75 days after the first HST/COS observation. Following this anomalous state, the flux and variability of the broad emission lines returned to a more normal state. This transient behavior, characterised by significant deficits in flux and equivalent width of the strong broad UV emission lines, is the first of its kind to be unambiguously identified in an active galactic nucleus reverberation mapping campaign. The largest corresponding emission-line flux deficits occurred for the high-ionization collisionally excited lines, C IV and Si IV(+O IV]), and also He II(+O III]), while the anomaly in Ly-alpha was substantially smaller. This pattern of behavior indicates a depletion in the flux of photons with E_{\rm ph} > 54 eV, relative to those near 13.6 eV. We suggest two plausible mechanisms for the observed behavior: (i) temporary obscuration of the ionizing continuum incident upon BLR clouds by a moving veil of material lying between the inner accretion disk and inner BLR, perhaps resulting from an episodic ejection of material from the disk, or (ii) a temporary change in the intrinsic ionizing continuum spectral energy distribution resulting in a deficit of ionizing photons with energies > 54 eV, possibly due to a transient restructuring of the Comptonizing atmosphere above the disk. Current evidence appears to favor the latter explanation.

We present ground-based optical photometric monitoring data for NGC 5548, part of an extended multi-wavelength reverberation mapping campaign. The light curves have nearly daily cadence from 2014 January to July in nine filters (\emph{BVRI} and \emph{ugriz}). Combined with ultraviolet data from the \emph{Hubble Space Telescope} and \emph{Swift}, we confirm significant time delays between the continuum bands as a function of wavelength, extending the wavelength coverage from 1158\,\AA\ to the $z$ band ($\sim\!9160$\,\AA). We find that the lags at wavelengths longer than the {\it V} band are equal to or greater than the lags of high-ionization-state emission lines (such as He\,{\sc ii}\,$\lambda 1640$ and $\lambda 4686$), suggesting that the continuum-emitting source is of a physical size comparable to the inner broad-line region (BLR). The trend of lag with wavelength is broadly consistent with the prediction for continuum reprocessing by an accretion disk with $\tau \propto \lambda^{4/3}$. However, the lags also imply a disk radius that is 3 times larger than the prediction from standard thin-disk theory, assuming that the bolometric luminosity is 10\% of the Eddington luminosity ($L = 0.1L_{\rm Edd}$). Using optical spectra from the Large Binocular Telescope, we estimate the bias of the interband continuum lags due to BLR emission observed in the filters. We find that the bias for filters with high levels of BLR contamination ($\sim\! 20\%$) can be important for the shortest continuum lags, and likely has a significant impact on the {\it u} and {\it U} bands owing to Balmer continuum emission.

We investigate if the active galactic nucleus (AGN) of Mrk 590, whose supermassive black hole was until recently highly accreting, is turning off due to a lack of central gas to fuel it. We analyse new sub-arcsecond resolution ALMA maps of the $^{12}$CO(3-2) line and 344 GHz continuum emission in Mrk 590. We detect no $^{12}$CO(3-2) emission in the inner 150 pc, constraining the central molecular gas mass to $M({\rm H_2}) \lesssim 1.6 \times 10^5\, {M_{\odot}}$, no more than a typical giant molecular gas cloud, for a CO luminosity to gas mass conversion factor of $\alpha_{\rm CO}\sim 0.8\,{M_{\odot}\,\rm (K \,km\,s^{-1}\,pc^{2}})^{-1}$. However, there is still potentially enough gas to fuel the black hole for another $2.6 \times 10^5$ years assuming Eddington-limited accretion. We therefore cannot rule out that the AGN may just be experiencing a temporary feeding break, and may turn on again in the near future. We discover a ring-like structure at a radius of $\sim 1$ kpc, where a gas clump exhibiting disturbed kinematics and located just $\sim 200$ pc west of the AGN, may be refueling the centre. Mrk 590 does not have significantly less gas than other nearby AGN host galaxies at kpc scales, confirming that gas reservoirs at these scales provide no direct indication of on-going AGN activity and accretion rates. Continuum emission detected in the central 150 pc likely originates from warm AGN-heated dust, although contributions from synchrotron and free-free emission cannot be ruled out.

In previous studies, it has been shown that the long term time average jet power, $\overline{Q}$, is correlated with the spectral index in the extreme ultraviolet (EUV), $\alpha_{EUV}$ (defined by $F_{\nu} \sim \nu^{-\alpha_{EUV}}$ computed between 700\AA\, and 1100\AA\,). Larger $\overline{Q}$ tends to decrease the EUV emission. This is a curious relationship because it connects a long term average over $\sim 10^{6}$ years with an instantaneous measurement of the EUV. The EUV appears to be emitted adjacent to the central supermassive black hole and the most straightforward explanation of the correlation is that the EUV emitting region interacts in real time with the jet launching mechanism. Alternatively stated, the $\overline{Q}$ - $\alpha_{EUV}$ correlation is a manifestation of a contemporaneous (real time) jet power, $Q(t)$, correlation with $\alpha_{EUV}$. In order to explore this possibility, this paper considers the time variability of the strong radio jet of the quasar 1442+101 that is not aberrated by strong Doppler enhancement. This high redshift (z = 3.55) quasar is uniquely suited for this endeavor as the EUV is redshifted into the optical observing window allowing for convenient monitoring. More importantly, it is bright enough to be seen through the Lyman forest and its radio flux is strong enough that it has been monitored frequently. Quasi-simultaneous monitoring (five epochs spanning $\sim 40$ years) show that increases in $Q(t)$ correspond to decreases in the EUV as expected.

As part of the OzDES spectroscopic survey we are carrying out a large scale reverberation mapping study of $\sim$500 quasars over five years in the 30 deg$^2$ area of the Dark Energy Survey (DES) supernova fields. These quasars have redshifts ranging up to 4 and have apparent AB magnitudes between $16.8<r<22.5$ mag. The aim of the survey is to measure time lags between fluctuations in the quasar continuum and broad emission line fluxes of individual objects in order to measure black hole masses for a broad range of AGN and constrain the radius-luminosity ($R-L$) relationship. Here we investigate the expected efficiency of the OzDES reverberation mapping campaign and its possible extensions. We expect to recover lags for $\sim$35-45\% of the quasars. AGN with shorter lags and greater variability are more likely to yield a lag, and objects with lags $\lesssim$6 months or $\sim$1 year are expected be recovered the most accurately. The baseline OzDES reverberation mapping campaign is predicted to produce an unbiased measurement of the $R-L$ relationship parameters for H$\beta$, Mg II $\lambda$2798, and C IV $\lambda$1549. However, extending the baseline survey by either increasing the spectroscopic cadence, extending the survey season, or improving the emission line flux measurement accuracy will significantly improve the $R-L$ parameter constraints for all broad emission lines.

Reverberation-mapping-based scaling relations are often used to estimate the masses of black holes from single-epoch spectra of AGN. While the radius-luminosity relation that is the basis of these scaling relations is determined using reverberation mapping of the H$\beta$ line in nearby AGN, the scaling relations are often extended to use other broad emission lines, such as MgII, in order to get black hole masses at higher redshifts when H$\beta$ is redshifted out of the optical waveband. However, there is no radius-luminosity relation determined directly from MgII. Here, we present an attempt to perform reverberation mapping using MgII in the well-studied nearby Seyfert 1, NGC 5548. We used Swift to obtain UV grism spectra of NGC 5548 once every two days from April to September 2013. Concurrent photometric UV monitoring with Swift provides a well determined continuum lightcurve that shows strong variability. The MgII emission line, however, is not strongly correlated with the continuum variability, and there is no significant lag between the two. We discuss these results in the context of using MgII scaling relations to estimate high-redshift black hole masses.

We describe the first results from a six-month long reverberation-mapping experiment in the ultraviolet based on 170 observations of the Seyfert 1 galaxy NGC 5548 with the Cosmic Origins Spectrograph on the Hubble Space Telescope. Significant correlated variability is found in the continuum and broad emission lines, with amplitudes ranging from ~30% to a factor of two in the emission lines and a factor of three in the continuum. The variations of all the strong emission lines lag behind those of the continuum, with He II 1640 lagging behind the continuum by ~2.5 days and Lyman alpha 1215, C IV 1550, and Si IV 1400 lagging by ~5-6 days. The relationship between the continuum and emission lines is complex. In particular, during the second half of the campaign, all emission-line lags increased by a factor of 1.3-2 and differences appear in the detailed structure of the continuum and emission-line light curves. Velocity-resolved cross-correlation analysis shows coherent structure in lag versus line-of-sight velocity for the emission lines; the high-velocity wings of C IV respond to continuum variations more rapidly than the line core, probably indicating higher velocity BLR clouds at smaller distances from the central engine. The velocity-dependent response of Lyman alpha, however, is more complex and will require further analysis.

Recent intensive Swift monitoring of the Seyfert 1 galaxy NGC 5548 yielded 282 usable epochs over 125 days across six UV/optical bands and the X-rays. This is the densest extended AGN UV/optical continuum sampling ever obtained, with a mean sampling rate <0.5 day. Approximately daily HST UV sampling was also obtained. The UV/optical light curves show strong correlations (r_max = 0.57 - 0.90) and the clearest measurement to date of interband lags. These lags are well-fit by a \tau propto \lambda^4/3 wavelength dependence, with a normalization that indicates an unexpectedly large disk radius of 0.35 +/- 0.05 lt-day at 1367 A, assuming a simple face-on model. The U-band shows a marginally larger lag than expected from the fit and surrounding bands, which could be due to Balmer continuum emission from the broad-line region as suggested by Korista and Goad. The UV/X-ray correlation is weaker (r_max < 0.45) and less consistent over time. This indicates that while Swift is beginning to measure UV/optical lags in general agreement with accretion disk theory (although the derived size is larger than predicted), the relationship with X-ray variability is less well understood. Combining this accretion disk size estimate with those from quasar microlensing studies suggests that AGN disk sizes scale approximately linearly with central black hole mass over a wide range of masses.

Quasi-stellar objects (QSOs) whose spectral energy distributions (SEDs) are reddened by dust either in their host galaxies or in intervening absorber galaxies are to a large degree missed by optical color selection criteria like the one used by the Sloan Digital Sky Survey (SDSS). To overcome this bias against red QSOs, we employ a combined optical and near-infrared color selection. In this paper, we present a spectroscopic follow-up campaign of a sample of red candidate QSOs which were selected from the SDSS and the UKIRT Infrared Deep Sky Survey (UKIDSS). The spectroscopic data and SDSS/UKIDSS photometry are supplemented by mid-infrared photometry from the Wide-field Infrared Survey Explorer. In our sample of 159 candidates, 154 (97%) are confirmed to be QSOs. We use a statistical algorithm to identify sightlines with plausible intervening absorption systems and identify nine such cases assuming dust in the absorber similar to Large Magellanic Cloud sightlines. We find absorption systems toward 30 QSOs, 2 of which are consistent with the best-fit absorber redshift from the statistical modeling. Furthermore, we observe a broad range in SED properties of the QSOs as probed by the rest-frame 2 {\mu}m flux. We find QSOs with a strong excess as well as QSOs with a large deficit at rest-frame 2 {\mu}m relative to a QSO template. Potential solutions to these discrepancies are discussed. Overall, our study demonstrates the high efficiency of the optical/near-infrared selection of red QSOs.

We investigate and quantify the observed scatter in the empirical relationship between the broad line region size R and the luminosity of the active galactic nucleus (AGN), in order to better understand its origin. This study is motivated by the indispensable role of this relationship in the mass estimation of cosmologically distant black holes, but may also be relevant to the recently proposed application of this relationship for measuring cosmic distances. We study six nearby reverberation-mapped AGN for which simultaneous UV and optical monitoring data exist. We also examine the long-term optical luminosity variations of Seyfert 1 galaxy NGC 5548 and employ Monte Carlo simulations to study the effects of the intrinsic variability of individual objects on the scatter in the global relationship for a sample of ~40 AGN. We find the scatter in this relationship has a correctable dependence on color. For individual AGN, the size of the Hbeta emitting region has a steeper dependence on the nuclear optical luminosity than on the UV luminosity, that can introduce a scatter of ~0.08 dex into the global relationship, due the non-linear relationship between the variations in the ionizing continuum and those in the optical continuum. Also, our analysis highlights the importance of understanding and minimizing the scatter in the relationship traced by the intrinsic variability of individual AGN, since it propagates directly into the global relationship. We find that using the UV luminosity as a substitute for the ionizing luminosity can reduce a sizable fraction of the current observed scatter of ~0.13 dex.

We present a revised measurement of the mass of the central black hole (Mbh) in the Seyfert 1 galaxy NGC 4151. The new stellar dynamical mass measurement is derived by applying an axisymmetric orbit-superposition code to near-infrared integral field data obtained using adaptive optics with the Gemini NIFS spectrograph. When our models attempt to fit both the NIFS kinematics and additional low spatial resolution kinematics, our results depend sensitively on how chi-squared is computed--probably a consequence of complex bar kinematics that manifest immediately outside the nuclear region. The most robust results are obtained when only the high spatial resolution kinematic constraints in the nuclear region are included in the fit. Our best estimates for the BH mass and H-band mass-to-light ratio are Mbh~(3.76+/-1.15)E7 Msun (1-sigma error) and M/L(H-band)~0.34+/-0.03 Msun/Lsun (3-sigma error), respectively (the quoted errors reflect the model uncertainties). Our BH mass measurement is consistent with estimates from both reverberation mapping (3.57[+0.45/-0.37]E7 Msun) and gas kinematics (3.0[+0.75/-2.2]E7 Msun; 1-sigma errors), and our best-fit mass-to-light ratio is consistent with the photometric estimate of M/L(H-band)=0.4+/-0.2 Msun/Lsun. The NIFS kinematics give a central bulge velocity dispersion sigma_c=116+/-3 km/s, bringing this object slightly closer to the M-sigma relation for quiescent galaxies. Although NGC 4151 is one of only a few Seyfert 1 galaxies in which it is possible to obtain a direct dynamical BH mass measurement--and thus, an independent calibration of the reverberation mapping mass scale--the complex bar kinematics makes it less than ideally suited for this purpose.

We investigate whether future measurements of high redshift standard candles (HzSCs) will be a powerful probe of dark energy, when compared to other types of planned dark energy measurements. Active galactic nuclei and gamma ray bursts have both been proposed as potential HzSC candidates. Due to their high luminosity, they can be used to probe unexplored regions in the expansion history of the universe. Information from these regions can help constrain the properties of dark energy, and in particular, whether it varies over time. We consider both linear and piecewise parameterizations of the dark energy equation of state, $w(z)$, and assess the optimal redshift distribution a high-redshift standard-candle survey could take to constrain these models. The more general the form of the dark energy equation of state $w(z)$ being tested, the more useful high-redshift standard candles become. For a linear parameterization of $w(z)$, HzSCs give only small improvements over planned supernova and baryon acoustic oscillation measurements; a wide redshift range with many low redshift points is optimal to constrain this linear model. However to constrain a general, and thus potentially more informative, form of $w(z)$, having many HzSCs can significantly improve limits on the nature of dark energy.

Short duration gamma-ray bursts (SGRBs) are thought to be related to the violent merger of compact objects, such as neutron stars or black holes, which makes them promising sources of gravitational waves. The detection of a 'kilonova'-like signature associated to the Swift-detected GRB 130603B has suggested that this event is the result of a compact object merger. Our knowledge on SGRB has been, until now, mostly based on the absence of supernova signatures and the analysis of the host galaxies to which they cannot always be securely associated. Further progress has been significantly hampered by the faintness and rapid fading of their optical counterparts (afterglows), which has so far precluded spectroscopy of such events. Afterglow spectroscopy is the key tool to firmly determine the distance at which the burst was produced, crucial to understand its physics, and study its local environment. Here we present the first spectra of a prototypical SGRB afterglow in which both absorption and emission features are clearly detected. Together with multiwavelength photometry we study the host and environment of GRB 130603B. From these spectra we determine the redshift of the burst to be z = 0.3565+/-0.0002, measure rich dynamics both in absorption and emission, and a substantial line of sight extinction of A_V = 0.86+/-0.15 mag. The GRB was located at the edge of a disrupted arm of a moderately star forming galaxy with near-solar metallicity. Unlike for most long GRBs (LGRBs), N_HX / A_V is consistent with the Galactic ratio, indicating that the explosion site differs from those found in LGRBs. The merger is not associated with the most star-forming region of the galaxy; however, it did occur in a dense region, implying a rapid merger or a low natal kick velocity for the compact object binary.

The narrow [O III] 4959, 5007 emission-line fluxes in the spectrum of the well-studied Seyfert 1 galaxy NGC 5548 are shown to vary with time. From this we show that the narrow line-emitting region has a radius of only 1-3 pc and is denser (n ~ 10^5 cm^{-3}) than previously supposed. The [O III] line width is consistent with virial motions at this radius given previous determinations of the black hole mass.Since the [O III] emission-line flux is usually assumed to be constant and is therefore used to calibrate spectroscopic monitoring data, the variability has ramifications for the long-term secular variations of continuum and emission-line fluxes, though it has no effect on shorter-term reverberation studies. We present corrected optical continuum and broad Hbeta emission-line light curves for the period 1988 to 2008.

We present a detailed study of the emission from a z = 2.35 galaxy that causes damped Lyman-alpha absorption in the spectrum of the background QSO, SDSS J 2222-0946. We present the results of extensive analyses of the stellar continuum covering the rest frame optical-UV regime based on broad-band HST imaging, and of spectroscopy from VLT/X-Shooter of the strong emission lines: Ly-alpha, [OII], [OIII], [NII], H-alpha and H-beta. We compare the metallicity from the absorption lines in the QSO spectrum with the oxygen abundance inferred from the strong-line methods (R23 and N2). The two emission-line methods yield consistent results: [O/H] = -0.30+/-0.13. Based on the absorption lines in the QSO spectrum a metallicity of -0.49+/-0.05 is inferred at an impact parameter of 6.3 kpc from the centre of the galaxy with a column density of hydrogen of log(N_HI)=20.65+/-0.05. The star formation rates of the galaxy from the UV continuum and H-alpha line can be reconciled assuming an amount of reddening of E(B-V) = 0.06+/-0.01, giving an inferred SFR of 13+/-1 M_sun / yr (assuming a Chabrier IMF). From the HST imaging, the galaxy associated with the absorption is found to be a compact (re=1.12 kpc) object with a disc-like, elongated (axis ratio 0.17) structure indicating that the galaxy is seen close to edge on. Moreover, the absorbing gas is located almost perpendicularly above the disc of the galaxy suggesting that the gas causing the absorption is not co-rotating with the disc. We investigate the stellar and dynamical masses from SED-fitting and emission-line widths, respectively, and find consistent results of 2x10^9 M_sun. We suggest that the galaxy is a young proto-disc with evidence for a galactic outflow of enriched gas. This galaxy hints at how star-forming galaxies may be linked to the elusive population of damped Lyman-alpha absorbers.

Comparison of six high-redshift quasar spectra obtained with the Large Binocular Telescope with previous observations from the Sloan Digital Sky Survey shows that failure to correctly identify absorption and other problems with accurate characterization of the CIV emission line profile in low S/N data can severely limit the reliability of single-epoch mass estimates based on the CIV emission line. We combine the analysis of these new high-quality data with a reanalysis of three other samples based on high S/N spectra of the CIV emission line region. We find that a large scatter between the Hbeta- and CIV-based masses remains even for this high S/N sample when using the FWHM to characterize the BLR velocity dispersion and the standard virial assumption to calculate the mass. However, we demonstrate that using high-quality data and the line dispersion to characterize the CIV line width leads to a high level of consistency between CIV- and Hbeta-based masses, with <0.3 dex of observed scatter, and an estimated ~0.2 dex intrinsic scatter, in the mass residuals.

The quasar Q0918+1636 (z=3.07) has two intervening high-metallicity Damped Lyman-alpha Absorbers (DLAs) along the line of sight, at redshifts of z=2.412 and 2.583. The z=2.583 DLA is located at a large impact parameter of 16.2 kpc, and despite this large impact parameter it has a very high metallicity (consistent with solar), a substantial fraction of H_2 molecules, and it is dusty as inferred from the reddened spectrum of the background QSO. The z=2.412 DLA has a metallicity of [M/H]=-0.6 (based on ZnII and SiII). In this paper we present new observations of this interesting sightline. HST/WFC3 imaging was obtained in the F606W, F105W and F160W bands. This is complemented by ground-based imaging in the u-, g-bands as well as K_s observations in the near-infrared (NIR). In addition, we present further spectroscopy with the ESO/VLT X-Shooter spectrograph. Based on these observations we obtain the following results: By fitting stellar population synthesis models to the photometric SED we constrain the physical properties of the z=2.583 DLA galaxy, and we infer its morphology by fitting a Sersic model to its surface brightness profile. We find it to be a relatively massive (M_star 10^10 M_sun), strongly star-forming (SFR~30 M_sun / yr, dusty (E_(B-V)=0.4) galaxy with a disk-like morphology. We detect most of the strong emission lines from the z=2.583 DLA [OIII],3727, [OIII],4960, [OIII],5007, Hbeta, and Halpha, albeit at low signal-to-noise (SN) ratio except for the [OIII],5007 line. We also detect [OIII],5007 emission from the galaxy counterpart of the z=2.412 DLA at a small impact parameter (<2 kpc). Overall our findings are consistent with the emerging picture that high-metallicity DLAs are associated with relatively (compared to typical DLAs) luminous and massive galaxy counterparts.

We present an updated and revised analysis of the relationship between the Hbeta broad-line region (BLR) radius and the luminosity of the active galactic nucleus (AGN). Specifically, we have carried out two-dimensional surface brightness decompositions of the host galaxies of 9 new AGNs imaged with the Hubble Space Telescope Wide Field Camera 3. The surface brightness decompositions allow us to create "AGN-free" images of the galaxies, from which we measure the starlight contribution to the optical luminosity measured through the ground-based spectroscopic aperture. We also incorporate 20 new reverberation-mapping measurements of the Hbeta time lag, which is assumed to yield the average Hbeta BLR radius. The final sample includes 41 AGNs covering four orders of magnitude in luminosity. The additions and updates incorporated here primarily affect the low-luminosity end of the R-L relationship. The best fit to the relationship using a Bayesian analysis finds a slope of alpha = 0.533 (+0.035/-0.033), consistent with previous work and with simple photoionization arguments. Only two AGNs appear to be outliers from the relationship, but both of them have monitoring light curves that raise doubt regarding the accuracy of their reported time lags. The scatter around the relationship is found to be 0.19(+/-0.02) dex, but would be decreased to 0.13 dex by the removal of these two suspect measurements. A large fraction of the remaining scatter in the relationship is likely due to the inaccurate distances to the AGN host galaxies. Our results help support the possibility that the R-L relationship could potentially be used to turn the BLRs of AGNs into standardizable candles. This would allow the cosmological expansion of the Universe to be probed by a separate population of objects, and over a larger range of redshifts.

Context: A significantly higher incidence of strong (rest equivalent width W_r > 1 {\AA}) intervening Mg II absorption is observed along gamma-ray burst (GRB) sight-lines relative to those of quasar sight-lines. A geometrical explanation for this discrepancy has been suggested: the ratio of the beam size of the source to the characteristic size of a Mg II absorption system can influence the observed Mg II equivalent width, if these two sizes are comparable. Aims: We investigate whether the differing beam sizes of the continuum source and broad-line region of Sloan Digital Sky Survey (SDSS) quasars produce a discrepancy between the incidence of strong Mg II absorbers illuminated by the quasar continuum region and those of absorbers illuminated by both continuum and broad-line region light. Methods: We perform a semi-automated search for strong Mg II absorbers in the SDSS Data Release 7 quasar sample. The resulting strong Mg II absorber catalog is available online. We measure the sight-line number density of strong Mg II absorbers superimposed on and off the quasar C IV 1550 {\AA} and C III] 1909 {\AA} emission lines. Results: We see no difference in the sight-line number density of strong Mg II absorbers superimposed on quasar broad emission lines compared to those superimposed on continuum-dominated spectral regions. This suggests that the Mg II-absorbing clouds typically observed as intervening absorbers in quasar spectra are larger than the beam sizes of both the continuum-emitting regions and broad line-emitting regions in the centers of quasars, corresponding to a lower limit of the order of 10^17} cm for the characteristic size of a Mg II absorbing cloud.

We present the results of a search for red QSOs using a selection based on optical imaging from SDSS and near-infrared imaging from UKIDSS. For a sample of 58 candidates 46 (79%) are confirmed to be QSOs. The QSOs are predominantly dust-reddened except a handul at redshifts z>3.5. The dust is most likely located in the QSO host galaxies. 4 (7%) of the candidates turned out to be late-type stars, and another 4 (7%) are compact galaxies. The remaining 4 objects we could not identify. In terms of their optical spectra the QSOs are similar to the QSOs selected in the FIRST-2MASS red Quasar survey except they are on average fainter, more distant and only two are detected in the FIRST survey. We estimate the amount of extinction using the SDSS QSO template reddened by SMC-like dust. It is possible to get a good match to the observed (restframe ultraviolet) spectra, but for nearly all the reddened QSOs it is not possible to match the near-IR photometry from UKIDSS. The likely reasons are that the SDSS QSO template is too red at optical wavelengths due to contaminating host galaxy light and that the assumed SMC extinction curve is too shallow. Our survey has demonstrated that selection of QSOs based on near-IR photometry is an efficent way to select QSOs, including reddened QSOs, with only small contamination from late-type stars and compact galaxies. This will be useful with ongoing and future wide-field near-IR surveys such as the VISTA and EUCLID surveys. [Abridged]

Accurate distances to celestial objects are key to establishing the age and energy density of the Universe and the nature of dark energy. A distance measure using active galactic nuclei (AGN) has been sought for more than forty years, as they are extremely luminous and can be observed at very large distances. We report here the discovery of an accurate luminosity distance measure using AGN. We use the tight relationship between the luminosity of an AGN and the radius of its broad line region established via reverberation mapping to determine the luminosity distances to a sample of 38 AGN. All reliable distance measures up to now have been limited to moderate redshift -- AGN will, for the first time, allow distances to be estimated to z~4, where variations of dark energy and alternate gravity theories can be probed.

I present mass functions of actively accreting black holes detected in different quasar surveys which in concert cover a wide range of cosmic history. I briefly address what we learn from these mass functions. I summarize the motivation for such a study and the methods by which we determine black hole masses.

We present the final results from a high sampling rate, multi-month, spectrophotometric reverberation mapping campaign undertaken to obtain either new or improved Hbeta reverberation lag measurements for several relatively low-luminosity AGNs. We have reliably measured thetime delay between variations in the continuum and Hbeta emission line in six local Seyfert 1 galaxies. These measurements are used to calculate the mass of the supermassive black hole at the center of each of these AGNs. We place our results in context to the most current calibration of the broad-line region (BLR) R-L relationship, where our results remove outliers and reduce the scatter at the low-luminosity end of this relationship. We also present velocity-resolved Hbeta time delay measurements for our complete sample, though the clearest velocity-resolved kinematic signatures have already been published.

We present an estimate of the black hole mass function (BHMF) of broad line quasars (BLQSOs) that self-consistently corrects for incompleteness and the statistical uncertainty in the mass estimates, based on a sample of 9886 quasars at 1 < z < 4.5 drawn from the Sloan Digital Sky Survey. We find evidence for `cosmic downsizing' of black holes in BLQSOs, where the peak in their number density shifts to higher redshift with increasing black hole mass. The cosmic mass density for black holes seen as BLQSOs peaks at z ~ 2. We estimate the completeness of the SDSS as a function of black hole mass and Eddington ratio, and find that at z > 1 it is highly incomplete at M_BH < 10^9 M_Sun and L / L_Edd < 0.5. We also estimate a lower limit on the lifetime of a single BLQSO phase and we place constraints on the maximum mass of a black hole in a BLQSO. Our estimated distribution of BLQSO Eddington ratios peaks at L / L_Edd ~ 0.05 and has a dispersion of ~ 0.4 dex, implying that most BLQSOs are not radiating at or near the Eddington limit; however the location of the peak is subject to considerable uncertainty. The steep increase in number density of BLQSOs toward lower Eddington ratios is expected if the BLQSO accretion rate monotonically decays with time. Furthermore, our estimated lifetime and Eddington ratio distributions imply that the majority of the most massive black holes spend a significant amount of time growing in an earlier obscured phase, a conclusion which is independent of the unknown obscured fraction. These results are consistent with models for self-regulated black hole growth, at least for massive systems at z > 1, where the BLQSO phase occurs at the end of a fueling event when black hole feedback unbinds the accreting gas, halting the accretion flow.

The most distant quasars known, at redshifts z=6, generally have properties indistinguishable from those of lower-redshift quasars in the rest-frame ultraviolet/optical and X-ray bands. This puzzling result suggests that these distant quasars are evolved objects even though the Universe was only seven per cent of its current age at these redshifts. Recently one z=6 quasar was shown not to have any detectable emission from hot dust, but it was unclear whether that indicated different hot-dust properties at high redshift or if it is simply an outlier. Here we report the discovery of a second quasar without hot-dust emission in a sample of 21 z=6 quasars. Such apparently hot-dust-free quasars have no counterparts at low redshift. Moreover, we demonstrate that the hot-dust abundance in the 21 quasars builds up in tandem with the growth of the central black hole, whereas at low redshift it is almost independent of the black hole mass. Thus z=6 quasars are indeed at an early evolutionary stage, with rapid mass accretion and dust formation. The two hot-dust-free quasars are likely to be first-generation quasars born in dust-free environments and are too young to have formed a detectable amount of hot dust around them.

We present results from a multi-month reverberation mapping campaign undertaken primarily at MDM Observatory with supporting observations from around the world. We measure broad line region (BLR) radii and black hole masses for six objects. A velocity-resolved analysis of the H_beta response shows the presence of diverse kinematic signatures in the BLR.

A detailed analysis of the data from a high sampling rate, multi-month reverberation mapping campaign, undertaken primarily at MDM Observatory with supporting observations from telescopes around the world, reveals that the Hbeta emission region within the broad line regions (BLRs) of several nearby AGNs exhibit a variety of kinematic behaviors. While the primary goal of this campaign was to obtain either new or improved Hbeta reverberation lag measurements for several relatively low luminosity AGNs (presented in a separate work), we were also able to unambiguously reconstruct velocity-resolved reverberation signals from a subset of our targets. Through high cadence spectroscopic monitoring of the optical continuum and broad Hbeta emission line variations observed in the nuclear regions of NGC 3227, NGC 3516, and NGC 5548, we clearly see evidence for outflowing, infalling, and virialized BLR gas motions, respectively.

We cross-correlate the SDSS DR3 quasar sample with FIRST and the Vestergaard et al. black hole (BH) mass sample to compare the mean accretion histories of optical and radio quasars. We find significant statistical evidence that radio quasars have a higher mean Eddington ratio Lambda at z > 2 with respect to optical quasars, while the situation is clearly reverse at z < 1. At z > 2 radio quasars happen to be less massive than optical quasars; however, as redshift decreases radio quasars appear in increasingly more massive BHs with respect to optical quasars. These two trends imply that radio sources are not a mere random subsample of optical quasars. No clear correlation between radio activity and BH mass and/or accretion rate is evident from our data, pointing to other BH properties, possibly the spin, as the driver of radio activity. We have checked that our main results do not depend on any evident bias. We perform detailed modelling of reasonable accretion histories for optical and radio quasars, finding that radio quasars grow by a factor of a few, at the most, since z ~ 4. The comparison between the predicted mass function of active radio quasars and the observed optical luminosity function of radio quasars, implies a significantly lower probability for lower mass BHs to be radio loud at all epochs, in agreement with what is observed in the local universe.

We present the first results from a high sampling rate, multi-month reverberation mapping campaign undertaken primarily at MDM Observatory with supporting observations from telescopes around the world. The primary goal of this campaign was to obtain either new or improved Hbeta reverberation lag measurements for several relatively low luminosity AGNs. We feature results for NGC 4051 here because, until now, this object has been a significant outlier from AGN scaling relationships, e.g., it was previously a ~2-3sigma outlier on the relationship between the broad-line region (BLR) radius and the optical continuum luminosity - the R_BLR-L relationship. Our new measurements of the lag time between variations in the continuum and Hbeta emission line made from spectroscopic monitoring of NGC 4051 lead to a measured BLR radius of R_BLR = 1.87 (+0.54 -0.50) light days and black hole mass of M_BH = 1.73 (+0.55 -0.52) x 10^6 M_sun. This radius is consistent with that expected from the R_BLR-L relationship, based on the present luminosity of NGC 4051 and the most current calibration of the relation by Bentz et al. (2009a). We also present a preliminary look at velocity-resolved Hbeta light curves and time delay measurements, although we are unable to reconstruct an unambiguous velocity-resolved reverberation signal.

We identify a sample of 74 high-redshift quasars (z>3) with weak emission lines from the Fifth Data Release of the Sloan Digital Sky Survey and present infrared, optical, and radio observations of a subsample of four objects at z>4. These weak emission-line quasars (WLQs) constitute a prominent tail of the Lya+NV equivalent width distribution, and we compare them to quasars with more typical emission-line properties and to low-redshift active galactic nuclei with weak/absent emission lines, namely BL Lac objects. We find that WLQs exhibit hot (T~1000 K) thermal dust emission and have rest-frame 0.1-5 micron spectral energy distributions that are quite similar to those of normal quasars. The variability, polarization, and radio properties of WLQs are also different from those of BL Lacs, making continuum boosting by a relativistic jet an unlikely physical interpretation. The most probable scenario for WLQs involves broad-line region properties that are physically distinct from those of normal quasars.

We present mass functions of distant actively accreting supermassive black holes residing in luminous quasars discovered in the Large Bright Quasar Survey, the Bright Quasar Survey, and the Fall Equatorial Stripe of the Sloan Digital Sky Survey (SDSS). The quasars cover a wide range of redshifts (0 <~ z <= 5) and were subject to different selection criteria and flux density limits. These samples are thus complementary and can help us gain additional insight on the true underlying black hole mass distribution, free from selection effects and mass estimation errors through future studies. We see evidence that the active z~4 black hole population is somewhat different than that at lower z. In particular, there is a sharp increase in the space density of the detected active black holes (M_BH >~ 10^8 Msun) between redshifts ~4 and ~2.5. Also, the z~4 SDSS quasar mass function has a somewhat flatter high mass-end slope, beta = -1.75 +- 0.56, compared to the mass functions based on quasars below z of 3, which display typical slopes of beta =~ -3.3; the latter are consistent with the mass functions at similar redshifts based on the SDSS Data Release 3 quasar catalog presented by Vestergaard et al. We see clear evidence of cosmic downsizing in the comoving space density distribution of active black holes in the LBQS sample alone. In forthcoming papers, further analysis, comparison, and discussion of these mass functions will be made with other existing black hole mass functions, notably that based on the SDSS DR3 quasar catalog. We present the relationships used to estimate the black hole mass based on the MgII emission line; the relations are calibrated to the Hbeta and CIV relations by means of several thousand high quality SDSS spectra. Mass estimates of the individual black holes of these samples are also presented.

A brief overview of the methods commonly used to determine or estimate the black hole mass in quiescent or active galaxies is presented and it is argued that the use of mass-scaling relations is both a reliable and the preferred method to apply to large samples of distant quasars. The method uses spectroscopic measurements of a broad emission-line width and continuum luminosity and currently has a statistical 1 sigma uncertainty in the absolute mass values of about a factor of 4. Potentially, this accuracy can be improved in the future. When applied to large samples of distant quasars it is evident that the black hole masses are very large, of order 1 to 10 billion solar masses, even at the highest redshifts of 4 to 6. The black holes must build up their mass very fast in the early universe. Yet they do not grow much larger than that: a maximum mass of about 10 billion solar masses is also observed. Preliminary mass functions of active black holes are presented for several quasar samples, including the Sloan Digital Sky Survey. Finally, common concerns related to the application of the mass scaling relations, especially for high redshift quasars, are briefly discussed.

We investigate the relationship between black hole mass and bulge luminosity for AGNs with reverberation-based black hole mass measurements and bulge luminosities from two-dimensional decompositions of Hubble Space Telescope host galaxy images. We find that the slope of the relationship for AGNs is 0.76-0.85 with an uncertainty of ~0.1, somewhat shallower than the M_BH \propto L^{1.0+/-0.1} relationship that has been fit to nearby quiescent galaxies with dynamical black hole mass measurements. This is somewhat perplexing, as the AGN black hole masses include an overall scaling factor that brings the AGN M_BH-sigma relationship into agreement with that of quiescent galaxies. We discuss biases that may be inherent to the AGN and quiescent galaxy samples and could cause the apparent inconsistency in the forms of their M_BH-L_bulge relationships.

We present high-resolution HST images of all 35 AGNs with optical reverberation-mapping results, which we have modeled to create a nucleus-free image of each AGN host galaxy. From the nucleus-free images, we determine the host-galaxy contribution to ground-based spectroscopic luminosity measurements at 5100A. After correcting the luminosities of the AGNs for the contribution from starlight, we re-examine the Hbeta R-L relationship. Our best fit for the relationship gives a powerlaw slope of 0.52 with a range of 0.45 - 0.59 allowed by the uncertainties. This is consistent with our previous findings, and thus still consistent with the naive assumption that all AGNs are simply luminosity-scaled versions of each other. We discuss various consistency checks relating to the galaxy modeling and starlight contributions, as well as possible systematic errors in the current set of reverberation measurements from which we determine the form of the R-L relationship.

We describe a Bayesian approach to estimating quasar black hole mass functions (BHMF) when using the broad emission lines to estimate black hole mass. We show how using the broad line mass estimates in combination with statistical techniques developed for luminosity function estimation leads to statistically biased results. We derive the likelihood function for the BHMF based on the broad line mass estimates, and derive the posterior distribution for the BHMF, given the observed data. We develop our statistical approach for a flexible model where the BHMF is modelled as a mixture of Gaussian functions. Statistical inference is performed using markov chain monte carlo (MCMC) methods. Our method has the advantage that it is able to constrain the BHMF even beyond the survey detection limits at the adopted confidence level, accounts for measurement errors and the intrinsic uncertainty in broad line mass estimates, and provides a natural way of estimating the probability distribution of any quantities derived from the BHMF. We conclude by using our method to estimate the local active BHMF using the z < 0.5 Bright Quasar Survey sources. At z = 0.2, the quasar BHMF falls off approximately as a power law with slope ~ 2 for M_{BH} > 10^8. Our analysis implies that z < 0.5 broad line quasars have a typical Eddington ratio of ~ 0.4 and a dispersion in Eddington ratio of < 0.5 dex (abridged).

We explore the nature of systematic errors that can arise in measurement of black hole masses from single-epoch spectra of active galactic nuclei (AGNs) by utilizing the many epochs available for NGC 5548 and PG1229+204 from reverberation mapping databases. In particular, we examine systematics due to AGN variability, contamination due to constant spectral components (i.e., narrow lines and host galaxy flux), data quality (i.e., signal-to-noise ratio, S/N), and blending of spectral features by comparing the precision and accuracy of single-epoch mass measurements to those of recent reverberation mapping studies. We calculate masses by characterizing the broad Hbeta emission line by both the full width at half maximum and the line dispersion and demonstrate the importance of removing narrow emission-line components and host starlight. We find that the reliability of line width measurements rapidly decreases for S/N lower than ~10 to 20 (per pixel) and that fitting the line profiles instead of direct measurement of the data does not mitigate this problem but can, in fact, introduce systematic errors. We also conclude that a full spectral decomposition to deblend the AGN and galaxy spectral features is unnecessary except to judge the contribution of the host galaxy to the luminosity and to deblend any emission lines that may inhibit accurate line width measurements. Finally, we present an error budget which summarizes the minimum observable uncertainties as well as the amount of additional scatter and/or systematic offset that can be expected from the individual sources of error investigated. In particular, we find that the minimum observable uncertainty in single-epoch mass estimates due to variability is <~ 0.1 dex for high S/N (>~ 20 per pixel) spectra.

We describe a Bayesian approach to estimating luminosity functions. We derive the likelihood function and posterior probability distribution for the luminosity function, given the observed data, and we compare the Bayesian approach with maximum-likelihood by simulating sources from a Schechter function. For our simulations confidence intervals derived from bootstrapping the maximum-likelihood estimate can be too narrow, while confidence intervals derived from the Bayesian approach are valid. We develop our statistical approach for a flexible model where the luminosity function is modeled as a mixture of Gaussian functions. Statistical inference is performed using Markov chain Monte Carlo (MCMC) methods, and we describe a Metropolis-Hastings algorithm to perform the MCMC. The MCMC simulates random draws from the probability distribution of the luminosity function parameters, given the data, and we use a simulated data set to show how these random draws may be used to estimate the probability distribution for the luminosity function. In addition, we show how the MCMC output may be used to estimate the probability distribution of any quantities derived from the luminosity function, such as the peak in the space density of quasars. The Bayesian method we develop has the advantage that it is able to place accurate constraints on the luminosity function even beyond the survey detection limits, and that it provides a natural way of estimating the probability distribution of any quantities derived from the luminosity function, including those that rely on information beyond the survey detection limits.

We report on 293 quasars with strong NIV] lambda 1486 or NIII] lambda 1750 emission lines (rest-frame equivalent width > 3 \AA) at 1.7 < z < 4.0 selected from the Sloan Digital Sky Survey (SDSS) Fifth Data Release. These nitrogen-rich (N-rich) objects comprise ~1.1% of the SDSS quasars. The comparison between the N-rich quasars and other quasars shows that the two quasar subsets share many common properties. We also confirm previous results that N-rich quasars have much stronger Lya and NV lambda 1240 emission lines. Strong nitrogen emission in all ionization states indicates high overall nitrogen abundances in these objects. We find evidence that the nitrogen abundance is closely related to quasar radio properties. The radio-loud fraction in the NIII]-rich quasars is 26% and in the NIV]-rich quasars is 69%, significantly higher than ~8% measured in other quasars with similar redshift and luminosity. Therefore, the high nitrogen abundance in N-rich quasars could be an indicator of a special quasar evolution stage, in which the radio activity is also strong.

We present the mass functions of actively accreting supermassive black holes over the redshift range 0.3 <= z <= 5 for a well-defined, homogeneous sample of 15,180 quasars from the Sloan Digital Sky Survey Data Release 3 (SDSS DR3) within an effective area of 1644 square degrees. This sample is the most uniform statistically significant subset available for the DR3 quasar sample. It was used for the DR3 quasar luminosity function, presented by Richards et al., and is the only sample suitable for the determination of the SDSS quasar black hole mass function. The sample extends from i = 15 to i = 19.1 at z less than about 3 and to i = 20.2 for z greater than about 3. The mass functions display a rise and fall in the space density distribution of active black holes at all epochs. Within the uncertainties the high-mass decline is consistent with a constant slope of beta of about -3.3 at all epochs. This slope is similar to the bright end slope of the luminosity function for epochs below z = 4. Our tests suggest that the down-turn toward lower mass values is due to incompleteness of the quasar sample with respect to black hole mass. Further details and analysis of these mass functions will be presented in forthcoming papers.

In this work we use a sample of 318 radio-quiet quasars (RQQ) to investigate the dependence of the ratio of optical/UV flux to X-ray flux, alpha_ox, and the X-ray photon index, Gamma_X, on black hole mass, UV luminosity relative to Eddington, and X-ray luminosity relative to Eddington. Our sample is drawn from the SDSS, with X-ray data from ROSAT and Chandra, and optical data mostly from the SDSS; 153 of these sources have estimates of Gamma_X from Chandra. We estimate M_BH using standard estimates derived from the Hbeta, Mg II, and C IV broad emission lines. Our sample spans a broad range in black hole mass (10^6 < M_BH / M_Sun < 10^10) and redshift (z < 4.8). We find that alpha_ox increases with increasing M_BH and L_UV / L_Edd, and decreases with increasing L_X / L_Edd. In addition, we confirm the correlation seen in previous studies between Gamma_X and M_BH and both L_UV / L_Edd and L_X / L_Edd; however, we also find evidence that the dependence of Gamma_X of these quantities is not monotonic, changing sign at M_BH ~ 3 x 10^8 M_Sun. We argue that the alpha_ox correlations imply that the fraction of bolometric luminosity emitted by the accretion disk, as compared to the corona, increases with increasing accretion rate relative to Eddington. In addition, we argue that the Gamma_X trends are caused by a dependence of X-ray spectral index on accretion rate. We discuss our results within the context of accretion models with comptonizing corona, and discuss the implications of the alpha_ox correlations for quasar feedback. To date, this is the largest study of the dependence of RQQ X-ray parameters on black hole mass and related quantities, and the first to attempt to correct for the large statistical uncertainty in the broad line mass estimates.

We present a stellar dynamical estimate of the black hole (BH) mass in the Seyfert 1 galaxy, NGC 4151. We analyze ground-based spectroscopy as well as imaging data from the ground and space, and we construct 3-integral axisymmetric models in order to constrain the BH mass and mass-to-light ratio. The dynamical models depend on the assumed inclination of the kinematic symmetry axis of the stellar bulge. In the case where the bulge is assumed to be viewed edge-on, the kinematical data give only an upper limit to the mass of the BH of ~4e7 M_sun (1 sigma). If the bulge kinematic axis is assumed to have the same inclination as the symmetry axis of the large-scale galaxy disk (i.e., 23 degrees relative to the line of sight), a best-fit dynamical mass between 4-5e7 M_sun is obtained. However, because of the poor quality of the fit when the bulge is assumed to be inclined (as determined by the noisiness of the chi^2 surface and its minimum value), and because we lack spectroscopic data that clearly resolves the BH sphere of influence, we consider our measurements to be tentative estimates of the dynamical BH mass. With this preliminary result, NGC 4151 is now among the small sample of galaxies in which the BH mass has been constrained from two independent techniques, and the mass values we find for both bulge inclinations are in reasonable agreement with the recent estimate from reverberation mapping (4.57[+0.57/-0.47]e7 M_sun) published by Bentz et al.

We present Gemini near-infrared spectroscopic observations of six luminous quasars at z=5.8$\sim$6.3. Five of them were observed using Gemini-South/GNIRS, which provides a simultaneous wavelength coverage of 0.9--2.5 $\mu$m in cross dispersion mode. The other source was observed in K band with Gemini-North/NIRI. We calculate line strengths for all detected emission lines and use their ratios to estimate gas metallicity in the broad-line regions of the quasars. The metallicity is found to be supersolar with a typical value of $\sim$4 Z_{\sun}, and a comparison with low-redshift observations shows no strong evolution in metallicity up to z$\sim$6. The FeII/MgII ratio of the quasars is 4.9+/-1.4, consistent with low-redshift measurements. We estimate central BH masses of 10^9 to 10^{10} M_{\sun} and Eddington luminosity ratios of order unity. We identify two MgII $\lambda\lambda$2796,2803 absorbers with rest equivalent width W_0^{\lambda2796}>1 \AA at 2.2<z<3 and three MgII absorbers with W_0^{\lambda2796}>1.5 \AA at z>3 in the spectra, with the two most distant absorbers at z=4.8668 and 4.8823, respectively. The redshift number densities (dN/dz) of MgII absorbers with W_0^{\lambda2796}>1.5 \AA are consistent with no cosmic evolution up to z>4.

We have investigated a sample of 5088 quasars from the Sloan Digital Sky Survey Second Data Release in order to determine how the frequency and properties of broad absorptions lines (BALs) depend on black hole mass, bolometric luminosity, Eddington fraction (L/L_Edd), and spectral slope. We focus only on high-ionization BALs and find a number of significant results. While quasars accreting near the Eddington limit are more likely to show BALs than lower $L/L_{Edd}$ systems, BALs are present in quasars accreting at only a few percent Eddington. We find a stronger effect with bolometric luminosity, such that the most luminous quasars are more likely to show BALs. There is an additional effect, previously known, that BAL quasars are redder on average than unabsorbed quasars. The strongest effects involving the quasar physical properties and BAL properties are related to terminal outflow velocity. Maximum observed outflow velocities increase with both the bolometric luminosity and the blueness of the spectral slope, suggesting that the ultraviolet luminosity to a great extent determines the acceleration. These results support the idea of outflow acceleration via ultraviolet line scattering.

We describe results from a new ground-based monitoring campaign on NGC 5548, the best studied reverberation-mapped AGN. We find that it was in the lowest luminosity state yet recorded during a monitoring program, namely L(5100) = 4.7 x 10^42 ergs s^-1. We determine a rest-frame time lag between flux variations in the continuum and the Hbeta line of 6.3 (+2.6/-2.3) days. Combining our measurements with those of previous campaigns, we determine a weighted black hole mass of M_BH = 6.54 (+0.26/-0.25) x 10^7 M_sun based on all broad emission lines with suitable variability data. We confirm the previously-discovered virial relationship between the time lag of emission lines relative to the continuum and the width of the emission lines in NGC 5548, which is the expected signature of a gravity-dominated broad-line region. Using this lowest luminosity state, we extend the range of the relationship between the luminosity and the time lag in NGC 5548 and measure a slope that is consistent with alpha = 0.5, the naive expectation for the broad line region for an assumed form of r ~ L^alpha. This value is also consistent with the slope recently determined by Bentz et al. for the population of reverberation-mapped AGNs as a whole.

We present low-resolution ultraviolet spectra of 14 low redshift (z<0.8) quasars observed with HST/STIS as part of a Snap project to understand the relationship between quasar outflows and luminosity. By design, all observations cover the CIV emission line. Nine of the quasars are from the Hamburg-ESO catalog, three are from the Palomar-Green catalog, and one is from the Parkes catalog. The sample contains a few interesting quasars including two broad absorption line (BAL) quasars (HE0143-3535, HE0436-2614), one quasar with a mini-BAL (HE1105-0746), and one quasar with associated narrow absorption (HE0409-5004). These BAL quasars are among the brightest known (though not the most luminous) since they lie at z<0.8. We compare the properties of these BAL quasars to the z<0.5 Palomar-Green and z>1.4 Large Bright Quasar samples. By design, our objects sample luminosities in between these two surveys, and our four absorbed objects are consistent with the v ~ L^0.62 relation derived by Laor & Brandt (2002). Another quasar, HE0441-2826, contains extremely weak emission lines and our spectrum is consistent with a simple power-law continuum. The quasar is radio-loud, but has a steep spectral index and a lobe-dominated morphology, which argues against it being a blazar. The unusual spectrum of this quasar resembles the spectra of the quasars PG1407+265, SDSSJ1136+0242, and PKS1004+13 for which several possible explanations have been entertained.

We present new observations leading to an improved black hole mass estimate for the Seyfert 1 galaxy NGC 4593 as part of a reverberation-mapping campaign conducted at the MDM Observatory. Cross-correlation analysis of the H_beta emission-line light curve with the optical continuum light curve reveals an emission-line time delay of 3.73 (+-0.75) days. By combining this time delay with the H_beta line width, we derive a central black hole mass of M_BH = 9.8(+-2.1)x10^6 M_sun, an improvement in precision of a factor of several over past results.

We present Spitzer observations of thirteen z~6 quasars using the Infrared Array Camera (IRAC) and Multiband Imaging Photometer for Spitzer (MIPS). All the quasars except SDSS J000552.34-000655.8 (SDSS J0005-0006) were detected with high S/N in the four IRAC channels and the MIPS 24um band, while SDSS J0005-0006 was marginally detected in the IRAC 8.0um band, and not detected in the MIPS 24um band. We find that most of these quasars have prominent emission from hot dust as evidenced by the observed 24um fluxes. Their spectral energy distributions (SEDs) are similar to those of low-redshift quasars at rest-frame 0.15-3.5 um, suggesting that accretion disks and hot-dust structures for these sources already have reached maturity. However, SDSS J0005-0006 has an unusual SED that lies significantly below low-redshift SED templates at rest-frame 1 and 3.5 um, and thus shows a strong near-IR (NIR) deficit and no hot-dust emission. Type I quasars with extremely small NIR-to-optical flux ratios like SDSS J0005-0006 are not found in low-redshift quasar samples, indicating that SDSS J0005-0006 has different dust properties at high redshift. We combine the Spitzer observations with X-ray, UV/optical, mm/submm and radio observations to determine bolometric luminosities for all the quasars. We find that the four quasars with central black-hole mass measurements have Eddington ratios of order unity.

We have undertaken a new ground-based monitoring campaign to improve the estimates of the mass of the central black hole in NGC 4151. We measure the lag time of the broad H beta line response compared to the optical continuum at 5100 A and find a lag of 6.6 (+1.1/-0.8) days. We combine our data with the recent reanalysis of UV emission lines by Metzroth et al. to calculate a weighted mean of the black hole mass, M_BH = 4.57 (+0.57/-0.47) x 10^7 M_sun. The absolute calibration of the black hole mass is based on normalization of the AGN black hole mass - stellar velocity dispersion (M_BH - sigma_*) relationship to that of quiescent galaxies by Onken et al. The scatter in the M_BH - sigma_* relationship suggests that reverberation-mapping based mass measurements are typically uncertain by a factor of 3-4.

We present an analysis of the mid-infrared (MIR) and optical properties of type 1 (broad-line) quasars detected by the Spitzer Space Telescope. The MIR color-redshift relation is characterized to z=3, with predictions to z=7. We demonstrate how combining MIR and optical colors can yield even more efficient selection of active galactic nuclei (AGN) than MIR or optical colors alone. Composite spectral energy distributions (SEDs) are constructed for 259 quasars with both Sloan Digital Sky Survey and Spitzer photometry, supplemented by near-IR, GALEX, VLA and ROSAT data where available. We discuss how the spectral diversity of quasars influences the determination of bolometric luminosities and accretion rates; assuming the mean SED can lead to errors as large as a factor of 2 for individual quasars. Finally, we show that careful consideration of the shape of the mean quasar SED and its redshift dependence leads to a lower estimate of the fraction of reddened/obscured AGNs missed by optical surveys as compared to estimates derived from a single mean MIR to optical flux ratio.

We try to identify the sources and systematics of the uncertainty of the masses of the central black holes in active galactic nuclei (AGNs) which have been obtained by reverberation-mapping methods. We characterize the broad H$\beta$ emission-line profiles by the ratio of their full-width at half maximum (FWHM) to their line dispersion, i.e., the second moment of the line profile. We use this parameter to separate the reverberation-mapped AGNs into two populations, the first with narrower H$\beta$ lines with relatively extended wings, and the second with broader lines that are relatively flat-topped. The first population is characterized by higher Eddington ratios than the second. Within each population, we calibrate the black-hole mass scale by comparison of the reverberation-based mass with that predicted by the bulge velocity dispersion. We also use the distribution of ratios of the reverberation-based mass to the velocity-dispersion mass prediction in a comparison with a ``generalized thick disk'' model in order to see if inclination can plausibly account for the observed distribution. We find that the line dispersion is a less biased parameter in general than FWHM for black hole mass estimation, although we show that it is possible to empirically correct for the bias introduced by using FWHM to characterize the emission-line width. We also argue that inclination effects are apparent in a small subset of objects with the narrowest emission lines. Our principal conclusion is that the H$\beta$ profile is sensitive primarily to Eddington ratio, but that inclination effects play a role in some cases. (abridged)

We present a spectroscopic survey using the MMT/Hectospec fiber spectrograph of 24 micron sources selected with the Spitzer Space Telescope in the Spitzer First Look Survey. We report 1296 new redshifts for 24 micron sources, including 599 with f(24micron) > 1 mJy. Combined with 291 additional redshifts for sources from the Sloan Digital Sky Survey (SDSS), our observing program was highly efficient and is ~90% complete for i' < 21 mag and f(24micron) > 1 mJy, and is 35% complete for i' < 20.5 mag and 0.3 mJy < f(24micron) < 1 mJy. Our Hectospec survey includes 1078 and 168 objects spectroscopically classified as galaxies and QSOs, respectively. Combining the Hectospec and SDSS samples, we find 24 micron-selected galaxies to z < 0.98 and QSOs to z < 3.6, with mean redshifts of <z(gal)> = 0.27 and <z(QSO)> =1.1. As part of this publication, we include the redshift catalogs and the reduced spectra; these are also available online (http://mips.as.arizona.edu/~papovich/fls) and through the NASA/IPAC Infrared Science Archive (http://irsa.ipac.caltech.edu).

We have obtained high resolution images of the central regions of 14 reverberation-mapped active galactic nuclei (AGN) using the Hubble Space Telescope Advanced Camera for Surveys High Resolution Camera to account for host-galaxy starlight contamination of measured AGN luminosities. We measure the host-galaxy starlight contribution to the continuum luminosity at 5100 A through the typical ground-based slit position and geometry used in the reverberation-mapping campaigns. We find that removing the starlight contribution results in a significant correction to the luminosity of each AGN, both for lower luminosity sources, as expected, but also for the higher luminosity sources such as the PG quasars. After accounting for the host galaxy starlight, we revisit the well-known broad-line region radius--luminosity relationship for nearby AGN. We find the power-law slope of the relationship for the H beta line to be 0.518 +/- 0.039, shallower than previously reported and consistent with the slope of 0.5 expected from the naive theoretical assumption that all AGN have, on average, the same ionizing spectrum and the same ionization parameter and gas density in the H beta line-emitting region.

We present four improved empirical relationships useful for estimating the central black hole mass in nearby AGNs and distant luminous quasars alike using either optical or UV single-epoch spectroscopy. These mass-scaling relationships between line widths and luminosity are based on recently improved empirical relationships between the broad-line region size and luminosities in various energy bands and are calibrated to the improved mass measurements of nearby AGNs based on emission-line reverberation mapping. The mass-scaling relationship based on the Hbeta line luminosity allows mass estimates for low-redshift sources with strong contamination of the optical continuum luminosity by stellar or non-thermal emission, while that based on the C IV lambda 1549 line dispersion allows mass estimates in cases where only the line dispersion (as opposed to the FWHM) can be reliably determined. We estimate that the absolute uncertainties in masses given by these mass-scaling relationships are typically around a factor of 4. We include in an Appendix mass estimates for all the Bright Quasar Survey (PG) quasars for which direct reverberation-based mass measurements are not available.

We present the discovery of seven quasars at z>5.7, selected from ~2000 deg^2 of multicolor imaging data of the Sloan Digital Sky Survey (SDSS). The new quasars have redshifts z from 5.79 to 6.13. Five are selected as part of a complete flux-limited sample in the SDSS Northern Galactic Cap; two have larger photometric errors and are not part of the complete sample. One of the new quasars, SDSS J1335+3533 (z=5.93), exhibits no emission lines; the 3-sigma limit on the rest-frame equivalent width of Ly alpha+NV line is 5 A. It is the highest redshift lineless quasar known, and could be a gravitational lensed galaxy, a BL Lac object or a new type of quasar. Two new z>6 quasars, SDSS 1250+3130 (z=6.13) and SDSS J1137+3549 (z=6.01), show deep Gunn-Peterson absorption gaps in Ly alpha. These gaps are narrower the complete Gunn-Peterson absorption troughs observed among quasars at z>6.2 and do not have complete Ly beta absorption.

We reinvestigate the relationship between the characteristic broad-line region size (R_blr) and the Balmer emission-line, X-ray, UV, and optical continuum luminosities. Our study makes use of the best available determinations of R_blr for a large number of active galactic nuclei (AGNs) from Peterson et al. Using their determinations of R_blr for a large sample of AGNs and two different regression methods, we investigate the robustness of our correlation results as a function of data sub-sample and regression technique. Though small systematic differences were found depending on the method of analysis, our results are generally consistent. Assuming a power-law relation R_blr \propto L^\alpha, we find the mean best-fitting \alpha is about 0.67+/-0.05 for the optical continuum and the broad H\beta luminosity, about 0.56+/-0.05 for the UV continuum luminosity, and about 0.70+/-0.14 for the X-ray luminosity. We also find an intrinsic scatter of about 40% in these relations. The disagreement of our results with the theoretical expected slope of 0.5 indicates that the simple assumption of all AGNs having on average same ionization parameter, BLR density, column density, and ionizing spectral energy distribution, is not valid and there is likely some evolution of a few of these characteristics along the luminosity scale.

We study the low-contrast Fe II emission blends in the ultraviolet (1250--2200A) and optical (4000--6000A) spectra of the Seyfert 1 galaxy NGC 5548 and show that these features vary in flux and that these variations are correlated with those of the optical continuum. The amplitude of variability of the optical Fe II emission is 50% - 75% that of Hbeta and the ultraviolet Fe II emission varies with an even larger amplitude than Hbeta. However, accurate measurement of the flux in these blends proves to be very difficult even using excellent Fe II templates to fit the spectra. We are able to constrain only weakly the optical Fe II emission-line response timescale to a value less than several weeks; this upper limit exceeds all the reliably measured emission-line lags in this source so it is not particularly meaningful. Nevertheless, the fact that the optical Fe II and continuum flux variations are correlated indicates that line fluorescence in a photoionized plasma, rather than collisional excitation, is responsible for the Fe II emission. The iron emission templates are available upon request.

Black-hole masses of distant quasars cannot be measured directly, but can be estimated to within a factor 3 to 5 using scaling relationships involving the quasar luminosity and broad-line width. Why such relationships are reasonable is summarized. The results of applying scaling relationships to data of quasars at a range of redshifts (z <= 6.3) are presented. Luminous quasars typically have masses of order 1 billion solar masses even at the highest redshifts. The fact that such massive black holes appear as early as at z ~ 6 indicate that black holes form very early or build up mass very fast.

We present improved black hole masses for 35 active galactic nuclei (AGNs) based on a complete and consistent reanalysis of broad emission-line reverberation-mapping data. From objects with multiple line measurements, we find that the highest precision measure of the virial product is obtained by using the cross-correlation function centroid (as opposed to the cross-correlation function peak) for the time delay and the line dispersion (as opposed to full width half maximum) for the line width and by measuring the line width in the variable part of the spectrum. Accurate line-width measurement depends critically on avoiding contaminating features, in particular the narrow components of the emission lines. We find that the precision (or random component of the error) of reverberation-based black hole mass measurements is typically around 30%, comparable to the precision attained in measurement of black hole masses in quiescent galaxies by gas or stellar dynamical methods. Based on results presented in a companion paper by Onken et al., we provide a zero-point calibration for the reverberation-based black hole mass scale by using the relationship between black hole mass and host-galaxy bulge velocity dispersion. The scatter around this relationship implies that the typical systematic uncertainties in reverberation-based black hole masses are smaller than a factor of three. We present a preliminary version of a mass-luminosity relationship that is much better defined than any previous attempt. Scatter about the mass-luminosity relationship for these AGNs appears to be real and could be correlated with either Eddington ratio or object inclination.

We calibrate reverberation-based black hole masses in active galactic nuclei (AGNs) by using the correlation between black hole mass, M, and bulge/spheroid stellar velocity dispersion, sigma. We use new measurements of sigma for 6 AGNs and published velocity dispersions for 10 others, in conjunction with improved reverberation mapping results, to determine the scaling factor required to bring reverberation-based black hole masses into agreement with the quiescent galaxy M-sigma relationship. The scatter in the AGN black hole masses is found to be less than a factor of 3. The current observational uncertainties preclude use of the scaling factor to discriminate between broad-line region models.

Mass estimates, based on scaling relationships, are presented of central black holes in luminous quasars at a range of redshifts (z < 0.5, 1.2 ~< z ~< 6.3). The data show that very massive (about 10^9 solar masses or larger) black holes appear already at z ~ 6, indicating that they form very early or very fast.

The applicability and apparent uncertainties of the techniques currently available for measuring or estimating black-hole masses in AGNs are briefly summarized.

This paper joins a series compiling consistent emission line measurements of large AGN spectral databases, useful for reliable statistical studies of emission line properties. It is preceded by emission line measurements of 993 spectra from the Large Bright Quasar Survey (Forster et al. 2001) and 174 spectra of AGN obtained from the Faint Object Spectrograph (FOS) on HST prior to the installation of COSTAR (Kuraszkiewicz et al. 2002). This time we concentrate on 220 spectra obtained with the FOS after the installation of COSTAR, completing the emission line analysis of all FOS archival spectra. We use the same automated technique as in previous papers, which accounts for Galactic extinction, models blended optical and UV iron emission, includes Galactic and intrinsic absorption lines and models emission lines using multiple Gaussians. We present UV and optical emission line parameters (equivalent widths, fluxes, FWHM, line positions) for a large number (28) of emission lines including upper limits for undetected lines. Further scientific analyses will be presented in subsequent papers.

A statistical study is presented of (a) the frequency of narrow CIV lambda 1549 absorption lines in 1.5 <~ z <~ 3.6 radio-quiet and radio-loud quasars, and of (b) the UV and radio properties of the absorbed quasars. The quasar sample is unbiased with respect to absorption properties and the radio-quiet and radio-loud subsamples are well matched in redshift and luminosity. A similarly high incidence (>~50%) of narrow CIV absorbers is detected for the radio-quiet and radio-loud quasars, and a constant ~25% of all the quasars, irrespective of radio type display associated CIV absorbers stronger than a rest equivalent width of 0.5A. Both radio-quiet and radio-loud quasars with narrow absorption lines have systematically redder continua, especially strongly absorbed objects. There is evidence of inclination dependent dust reddening and absorption for the radio quasars. An additional key result is that the most strongly absorbed radio quasars have the largest radio source extent. This result is in stark contrast to a recent study of the low-frequency selected Molonglo survey in which a connection between the strength of the narrow absorbers and the (young) age of the radio source has been proposed. The possible origin of these discrepant results is discussed and may be related to the higher source luminosity for the quasars studied here.

Black-hole masses of the highest redshift quasars (4 <~ z <~ 6) are estimated using a previously presented scaling relationship, derived from reverberation mapping of nearby quasars, and compared to quasars at lower redshift. It is shown that the central black holes in luminous z >~ 4 quasars are very massive (>~ 10^9 solar masses). It is argued that the mass estimates of the high-z quasars are not subject to larger uncertainties than those for nearby quasars. Specifically, the large masses are not overestimates and the lack of similarly large black-hole masses in the nearby Universe does not rule out their existence at high-z. However, AGN host galaxies do not typically appear fully formed or evolved at these early epochs. This supports scenarios in which black holes build up mass very fast in a radiatively inefficient (or obscured) phase relative to the stars in their galaxies. Additionally, upper envelopes of black-hole mass of approximately 10^{10} solar masses and bolometric luminosity of ~ 10^{48} erg/s are observed at all redshifts.

We present results of the analysis of near infrared spectroscopic observations of 6 high-redshift quasars (z > 4), emphasizing the measurement of the ultraviolet FeII/MgII emission line strength in order to estimate the beginning of intense star formation in the early universe. To investigate the evolution of the FeII/MgII ratio over a wider range in cosmic time, we measured this ratio for composite quasar spectra which cover a redshift range of 0 < z < 5 with nearly constant luminosity, as well as for those which span ~6 orders of magnitude in luminosity. A detailed comparison of the high-redshift quasar spectra with those of low-redshift quasars with comparable luminosity shows essentially the same FeII/MgII emission ratios and very similar continuum and line spectral properties, i.e. a lack of evolution of the relative iron to magnesium abundance of the gas in bright quasars since z=5. Current nucleosynthesis and stellar evolution models predict that alpha-elements like magnesium are produced in massive stars ending in type II SNe, while iron is formed predominantly in SNe of type Ia with intermediate mass progenitors. This results in an iron enrichment delay of 0.2 to 0.6 Gyr. We conclude that intense star formation activity in the host galaxies of z>4 quasars must have started already at an epoch corresponding to z_f = 6 to 9, when the age of the universe was ~0.5 Gyr (H_o = 72 km/s/Mpc, Omega_M = 0.3, Omega_Lambda = 0.7). This epoch corresponds well to the re-ionization era of the universe.

We examine rest-frame ultraviolet spectra of 70 high redshift quasars (z>3.5) to study the chemical enrichment history of the gas closely related to the quasars, and thereby estimate the epoch of first star formation. The fluxes of several ultraviolet emission lines were investigated within the framework of the most recent photoionization models to estimate the metallicity of the gas associated with the high-z quasars. Standard photoionization parameters and the assumption of secondary nitrogen enrichment indicate an average abundance of Z/Z_sol = 4 to 5 in the line emitting gas. Assuming a time scale of t_evol = 0.5 - 0.8 Gyrs for the chemical enrichment of the gas, the first major star formation for quasars with z>=4 should have started at a redshift of z_f = 6 - 8, corresponding to an age of the universe of several 10^8 yrs (H_o = 65 km/s/Mpc, Omega_M = 0.3, Omega_Lambda = 0.7). We note that this also appears to be the era of re-ionization of the universe. Finally, there is some evidence for a positive luminosity - metallicity relation in this high redshift quasar sample.

We present observations of 11 high redshift quasars ($3.9 \la z \la 5.0$) observed with low spectral resolution in the restframe ultraviolet using FORS 1 at the VLT UT 1. The emission-line fluxes of strong permitted and intercombination ultraviolet emission lines are measured to estimate the chemical composition of the line emitting gas. Comparisons to photoionization calculations indicate gas metallicities in the broad emission line region in the range of solar to several times solar. The average of the mean metallicity of each high-z quasar in this sample is $Z/Z_\odot = 4.3 \pm 0.3$. Assuming a chemical evolution time scale of $\tau_{evol} \simeq 0.5 - 0.8$ Gyrs, we derive a redshift of $z_f \simeq 6 {\rm to} 8$ for the onset of the first major star formation episode (H$_o = 65$ km s$^{-1}$ Mpc$^{-1}$, $\Omega_M = 0.3$, $\Omega_\Lambda = 0.7$), corresponding to an age of the universe of several $10^8$ yrs at this epoch. We note that this epoch is also supposed to be the era of re-ionization of the universe.

We report on the analysis of a large sample of 744 type 1 Active Galactic Nuclei, including quasars and Seyfert 1 galaxies across the redshift range from 0 < z < 5 and spanning nearly 6 orders of magnitude in continuum luminosity. We discuss correlations of continuum and emission line properties in the rest-frame ultraviolet and optical spectral ranges. The well established Baldwin Effect is detected for almost all emission lines from OVI1034 to [OIII]5007. Their equivalent widths are significantly anti-correlated with the continuum strength, while they are nearly independent of redshift. This is the well known Baldwin Effect. Its slope beta, measured as log W_lambda ~ beta * log lambda * L_lambda (1450A), shows a tendency to become steeper towards higher luminosity. The slope of the Baldwin Effect also increases with the ionization energy needed to create the individual lines. In contrast to this general trend, the NV1240 equivalent width is nearly independent of continuum luminosity and remains nearly constant. The overall line behaviors are consistent with softer UV continuum shapes and perhaps increasing gas metallicity in more luminous Active Galactic Nuclei.

An empirical relationship, of particular interest for studies of high redshift active galactic nuclei (AGNs) and quasars, between the masses of their central black-holes and rest-frame ultraviolet (UV) parameters measured in single-epoch AGN spectra is presented. This relationship is calibrated to recently measured reverberation masses of low-redshift AGNs and quasars. An empirical relationship between single-epoch rest-frame optical spectrophotometric measurements and the central masses is also presented. The UV relationship allows reasonable estimates of the central masses to be made of high-redshift AGNs and quasars for which these masses cannot be directly or easily measured by the techniques applicable to the lower luminosity, nearby AGNs. The central mass obtained by this method can be estimated to within a factor of ~3 for most objects. This is reasonable given the intrinsic uncertainty of a factor less than 2 in the primary methods used to measure the central masses of nearby inactive and active galaxies, namely resolved gas and stellar kinematics in the underlying host galaxy and reverberation-mapping techniques. The UV relationship holds good potential for being a powerful tool to study black-hole demographics at high redshift as well as to statistically study the fundamental properties of AGNs. The broad line region size - luminosity relationship is key to the calibrations presented here. The fact that its intrinsic scatter is also the main source of uncertainty in the calibrations stresses the need for better observational constraints to be placed on this relationship. The empirically calibrated relationships presented here will be applied to quasar samples in forthcoming work.

A large, high-quality spectral data base of well-selected, moderate-redshift radio-loud and radio-quiet quasars is used to characterize the incidence of narrow associated CIV lambda 1549 absorption, and how this may depend on some quasar properties, including radio-type. Preliminary results of this study are presented.

In order to derive information on the star formation history in the early universe we observed 6 high-redshift (z=3.4) quasars in the near-infrared to measure the relative iron and \mgii emission strengths. A detailed comparison of the resulting spectra with those of low-redshift quasars show essentially the same FeII/MgII emission ratios and very similar continuum and line spectral properties, indicating a lack of evolution of the relative iron to magnesium abundance of the gas since z=3.4 in bright quasars. On the basis of current chemical evolution scenarios of galaxies, where magnesium is produced in massive stars ending in type II SNe, while iron is formed predominantly in SNe of type Ia with a delay of ~1 Gyr and assuming as cosmological parameters H_o = 72 km/s Mpc, Omega_M = 0.3, and Omega_Lambda = 0.7$, we conclude that major star formation activity in the host galaxies of our z=3.4 quasars must have started already at an epoch corresponding to z_f ~= 10, when the age of the universe was less than 0.5 Gyrs.

We present the results of a two-year long optical monitoring program of the narrow-line Seyfert 1 galaxy Akn 564. The majority of this monitoring project was also covered by X-ray observations (RXTE) and for a period of ~50 days, we observed the galaxy in UV (HST) and X-rays (RXTE & ASCA) simultaneously with the ground-based observations. Rapid and large-amplitude variations seen in the X-ray band, on a daily and hourly time-scale, were not detected at optical and UV wavelengths, which in turn exhibited much lower variability either on short (one day) or long (several months) time-scales. The only significant optical variations can be described as two 2--4 day events with ~10% flux variations. We detect no significant optical line variations and thus cannot infer a reverberation size for the broad-line region. Similarly, the large X-ray variations seem to vanish when the light curve is smoothed over a period of 30 days. The UV continuum follows the X-rays with a lag of ~0.4 days, and the optical band lags the UV band by ~2 days. No significant correlation was found between the entire X-ray dataset and the optical band. Focusing on a 20-day interval around the strongest optical event we detect a significant X-ray--optical correlation with similar events seen in the UV and X-rays. Our data are consistent with reprocessing models on the grounds of the energy emitted in this single event. However, several large X-ray flares produced no corresponding optical emission.

Iron emission is often a severe contaminant in optical-ultraviolet spectra of active galaxies. Its presence complicates emission line studies. A viable solution, already successfully applied at optical wavelengths, is to use an empirical iron emission template. We have generated FeII and FeIII templates for ultraviolet active galaxy spectra based on HST archival 1100 - 3100 A spectra of IZw1. Their application allows fitting and subtraction of the iron emission in active galaxy spectra. This work has shown that in particular CIII] lambda 1909 can be heavily contaminated by other line emission, including iron transitions. Details of the data processing, generation, and use of the templates, are given by Vestergaard & Wilkes (2001).

We present an empirical template spectrum suitable for fitting/subtracting and studying the FeII and FeIII line emission in the restframe UV spectra of active galatic nuclei (AGNs), the first empirical UV iron template to cover the full 1250 - 3090 A range. Iron emission is often a severe contaminant in optical--UV spectra of AGNs. Its presence complicates and limits the accuracy of measurements of both strong and weak emission lines and the continuum emission, affecting studies of line and continuum interrelations, the ionization structure, and elemental abundances in AGNs. Despite the wealth of work on modeling the AGN FeII emission and the need to account for it in observed AGN spectra, there is no UV template electronically available to aid this process. The iron template we present is based on HST spectra of the Narrow Line Seyfert 1, IZw1. Its intrinsic narrow lines (~900 km/s) and rich iron spectrum make the template particularly suitable for use with most AGN spectra. The iron emission spectrum and the line identifications and measurements are presented and compared with the work of Laor et al. We illustrate the application of the derived FeII and FeIII templates by fitting and subtracting the iron emission from UV spectra of four high-z quasars and of the nearby quasar, 3C273. We briefly discuss the small discrepancies between this observed iron emission and the UV template, and compare the template with previously published ones. We discuss the advantages and limitations of the templates and of the template fitting method. We conclude that the templates work sufficiently well to be a valuable and important tool for eliminating and studying the iron emission in AGNs, at least until accurate theoretical iron emission models are developed. (Abridged)

We present measurements of the optical/UV emission lines for a large homogeneous sample of 993 quasars from the Large Bright Quasar Survey. Our largely automated technique accounts for continuum breaks and galactic reddening, and we perform multicomponent fits to emission line profiles, including the effects of blended iron emission, and of absorption lines both galactic and intrinsic. Here we describe the fitting algorithm and present the results of line fits to the LBQS sample, including upper limits to line equivalent widths when warranted. The distribution of measured line parameters, principally equivalent width and FWHM, are detailed for a variety of lines, including upper limits. We thus initiate a large-scale investigation of correlations between the high energy continuum and emission lines in quasars, to be extended to complementary samples using similar techniques. High quality, reproducible measurements of emission lines for uniformly selected samples will advance our understanding of active galaxies, especially in a new era of large surveys selected by a variety of complementary methods.

We report the discovery of the gamma-ray burst GRB 000131 and its optical afterglow. The optical identification was made with the VLT 84 hours after the burst following a BATSE detection and an Inter Planetary Network localization. GRB 000131 was a bright, long-duration GRB, with an apparent precursor signal 62 s prior to trigger. The afterglow was detected in ESO VLT, NTT, and DK1.54m follow-up observations. Broad-band and spectroscopic observations of the spectral energy distribution reveals a sharp break at optical wavelengths which is interpreted as a Ly-alpha absorption edge at 6700 A. This places GRB 000131 at a redshift of 4.500 +/- 0.015. The inferred isotropic energy release in gamma rays alone was approximately 10^54 erg (depending on the assumed cosmology). The rapid power-law decay of the afterglow (index alpha=2.25, similar to bursts with a prior break in the lightcurve), however, indicates collimated outflow, which relaxes the energy requirements by a factor of < 200. The afterglow of GRB 000131 is the first to be identified with an 8-m class telescope.

We present evidence that the high-velocity CIV lambda 1549 emission line gas of radio-loud quasars may originate in a disk-like configuration, in close proximity to the accretion disk often assumed to emit the low-ionization lines. For a sample of 36 radio-loud z~2 quasars we find the 20--30% peak width to show significant inverse correlations with the fractional radio core-flux density, R, the radio axis inclination indicator. Highly inclined systems have broader line wings, consistent with a high-velocity field perpendicular to the radio axis. By contrast, the narrow line-core shows no such relation with R, so the lowest velocity CIV-emitting gas has an inclination independent velocity field. We propose that this low-velocity gas is located at higher disk-altitudes than the high-velocity gas. A planar origin of the high-velocity CIV-emission is consistent with the current results and with an accretion disk-wind emitting the broad lines. A spherical distribution of randomly orbiting broad-line clouds and a polar high-ionization outflow are ruled out.

We report on the results of a three-year program of coordinated X-ray and optical monitoring of the narrow-line Seyfert 1 galaxy NGC 4051. The rapid continuum variations observed in the X-ray spectra are not detected in the optical, although the time-averaged X-ray and optical continuum fluxes are well-correlated. Variations in the flux of the broad Hbeta line are found to lag behind the optical continuum variations by 6 days (with an uncertainty of 2-3 days), and combining this with the line width yields a virial mass estimate of about 1.1 million solar masses, at the very low end of the distribution of AGN masses measured by line reverberation. Strong variability of He II 4686 is also detected, and the response time measured is similar to that of Hbeta, but with a much larger uncertainty. The He II 4686 line is almost five times broader than Hbeta, and it is strongly blueward asymmetric, as are the high-ionization UV lines recorded in archival spectra of NGC 4051. The data are consistent with the Balmer lines arising in a low to moderate inclination disk-like configuration, and the high-ionization lines arising in an outflowing wind, of which we observe preferentially the near side. Previous observations of the narrow-line region morphology of this source suggest that the system is inclined by about 50 degrees, and if this is applicable to the broad Hbeta-emitting region, a central mass of about 1.4 million solar masses can be inferred. During the third year of monitoring, both the X-ray continuum and the He II 4686 line went into extremely low states, although the optical continuum and the Hbeta broad line were both still present and variable. We suggest that the inner part of the accretion disk may have gone into an advection-dominated state, yielding little radiation from the hotter inner disk.

We present supportive evidence that the Boroson and Green eigenvector 1 is not driven by source orientation. Until recently it was generally accepted that eigenvector 1 does not depend on orientation as it strongly correlates with [OIII]5007 emission, thought to be an isotropic property. However, recent studies of radio-loud AGN have questioned the isotropy of [OIII] emission and concluded that [OII]3727 emission is isotropic. In this paper we investigate the relation between eigenvector 1 and [OII] emission in radio-quiet BQS (Bright Quasar Survey) quasars, and readdress the issue of orientation as the driver of eigenvector 1. We find significant correlations between eigenvector 1 and orientation independent [OII] emission, which implies that orientation does not drive eigenvector 1. The luminosities and equivalent widths of [OIII] and [OII] correlate with one another, and the range in luminosities and equivalent widths is similar. This suggests that the radio-quiet BQS quasars are largely free of orientation dependent dust effects and ionization dependent effects in the narrow-line region. We also conclude that neither the [OIII] emission nor the [OII]/[OIII] ratio are dependent on orientation in our radio-quiet BQS quasar sample, contrary to recent results found for radio-loud quasars.

We study the UV spectra of NLS1 galaxies and compare them with typical Seyfert 1 galaxies and quasars. The NLS1 spectra show narrower UV lines as well as weaker CIV lambda 1549 and CIII] lambda 1909 emission. We show that these line properties are due to a lower ionization parameter and somewhat higher BLR cloud densities. These modified conditions can be explained by the hotter big blue bumps observed in NLS1s, which are in turn due to higher L/L_Edd ratios, as shown by our accretion disk and corona modeling of the NLS1 continua. We also present evidence that the Boroson & Green eigenvector 1, which is correlated with the optical and UV emission-line properties, is not driven by orientation and hence NLS1s, which have extreme eigenvector 1 values, are not viewed from an extreme viewing angle.

VLA imaging at kiloparsec-scale resolution of sixteen core-dominated radio-loud QSOs is presented. Many objects appear to display variable radio emission and their radio morphologies are significantly smaller than those of steep-spectrum quasars, consistent with these objects being observed at sight lines close to their (relativistic, $\gamma \approx$ 4-7) jet axes. The usefulness of the radio source orientation indicator R_V, being defined as ratio of radio core and rest frame optical V-band luminosity, is confirmed.

We report the discovery of X-ray emission from a radio knot at a projected distance of 25" from the nucleus of the Seyfert galaxy, 3C 120. The data were obtained with the ROSAT High Resolution Imager (HRI). Optical upper limits for the knot preclude a simple power law extension of the radio spectrum and we calculate some of the physical parameters for thermal bremsstrahlung and synchrotron self-Compton models. We conclude that no simple model is consistent with the data but if the knot contains small regions with flat spectra, these could produce the observed X-rays (via synchrotron emission) without being detected at other wavebands.

We report results from five years of monitoring of the Einstein Cross (Q 2237+0305) with the Nordic Optical Telescope. The photometry, mainly in the R and I bands, has been performed by a PSF fitting and `cleaning' procedure, in which the four image components as well as the lens galaxy and its nucleus are iteratively removed. The resulting lightcurves exhibit several microlensing features; one event may have a timescale as short as 14 days. Variations of timescales of several years are found in all four images. This becomes even more convincing when our data are combined with data published for 1986-89. No clear high amplification event was observed during the period. A brightening of all four components during 1994 is interpreted as intrinsic variation.

We report the detection of an optical jet in the nearby Seyfert 1 radio galaxy 3C 120. The optical jet coincides with the well-known radio jet and emits continuum radiation (B,V',I) with a radio-to-optical spectral index of 0.65. There are no clear optical counterparts to the radio knots, although the optical condensation A of the galaxy, which includes the bright 4" radio knot, is found to be 12 % polarized with the electric field vectors perpendicular to the jet. These findings indicate that 3C 120 contains the 6th known extragalactic optical synchrotron jet, quite similar in its properties to the jet of PKS 0521-36. The outer parts of the jet is the faintest known optical jet and was discovered as the result of a dedicated effort to detect it. It is therefore possible that more optical jets can be discovered in systematic searches by combining deep imaging in the optical or near-IR with careful galaxy subtraction methods.