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.

Sub-damped Lyman alpha systems (subDLAs; HI column densities of 19.0<=logN(HI)<20.3) are rarely included in the cosmic HI census performed at redshifts z>=1.5, yet are expected to contribute significantly to the overall HI mass budget of the Universe. In this paper, we present a blindly selected sample of 155 subDLAs found along 100 quasar sightlines (with a redshift path length X=475) in the XQ-100 survey to investigate the contribution of subDLAs to the HI mass density of the Universe. The impact of X-Shooter's spectral resolution on sub-DLA identification is evaluated, and found to be sufficient for reliably finding absorbers with logN(HI)>=18.9. We compared the implications of searching for subDLAs solely using HI absorption versus the use of additional metal lines to confirm the identification, and found that metal-selection techniques would have missed 75 subDLAs. Using a bootstrap-Monte Carlo simulation, we computed the column density distribution function (f(N,X)) and the cosmological HI mass density of subDLAs and compared with our previous work based on the XQ-100 damped Lyman alpha systems. We do not find any significant redshift evolution in f(N,X) or cosmological HI mass density for subDLAs. However, subDLAs contribute 10-20 per cent of the total cosmological HI mass density measured at redshifts 2<z<5 (agreeing with previous measurements), and thus have a small but significant contribution to the HI budget of the Universe.

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 reverberation mapping results from the first year of combined spectroscopic and photometric observations of the Sloan Digital Sky Survey Reverberation Mapping Project. We successfully recover reverberation time delays between the $g+i$-band emission and the broad H$\beta$ emission line for a total of 44 quasars, and for the broad H$\alpha$ emission line in 18 quasars. Time delays are computed using the JAVELIN and CREAM software and the traditional interpolated cross-correlation function (ICCF): Using well defined criteria, we report measurements of 32 H$\beta$ and 13 H$\alpha$ lags with JAVELIN, 42 H$\beta$ and 17 H$\alpha$ lags with CREAM, and 16 H$\beta$ and 8 H$\alpha$ lags with the ICCF. Lag values are generally consistent among the three methods, though we typically measure smaller uncertainties with JAVELIN and CREAM than with the ICCF, given the more physically motivated light curve interpolation and more robust statistical modeling of the former two methods. The median redshift of our H$\beta$-detected sample of quasars is 0.53, significantly higher than that of the previous reverberation-mapping sample. We find that in most objects, the time delay of the H$\alpha$ emission is consistent with or slightly longer than that of H$\beta$. We measure black hole masses using our measured time delays and line widths for these quasars. These black hole mass measurements are mostly consistent with expectations based on the local M-sigma relationship, and are also consistent with single-epoch black hole mass measurements. This work increases the current sample size of reverberation-mapped active galaxies by about two-thirds and represents the first large sample of reverberation mapping observations beyond the local universe (z < 0.3).

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.

We investigate the NV absorption signal along the line of sight of background quasars, in order to test the robustness of the use of this ion as criterion to select intrinsic (i.e. physically related to the quasar host galaxy) narrow absorption lines (NALs). We build composite spectra from a sample of $\sim$ 1000 CIV absorbers, covering the redshift range 2.55 < z < 4.73, identified in 100 individual sight lines from the XQ-100 Legacy Survey. We detect a statistical significant NV absorption signal only within 5000 km s$^{-1}$ of the systemic redshift, z$\rm_{em}$. This absorption trough is $\sim$ 15$\sigma$ when only CIV systems with N(CIV) > 10$^{14}$ cm$^{-2}$ are included in the composite spectrum. This result confirms that NV offers an excellent statistical tool to identify intrinsic systems. We exploit the stacks of 11 different ions to show that the gas in proximity to a quasar exhibits a considerably different ionization state with respect to gas in the transverse direction and intervening gas at large velocity separations from the continuum source. Indeed, we find a dearth of cool gas, as traced by low-ionization species and in particular by MgII, in the proximity of the quasar. We compare our findings with the predictions given by a range of Cloudy ionization models and find that they can be naturally explained by ionization effects of the quasar.

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 optical continuum lags for two Seyfert 1 galaxies, MCG+08-11-011 and NGC 2617, using monitoring data from a reverberation mapping campaign carried out in 2014. Our light curves span the ugriz filters over four months, with median cadences of 1.0 and 0.6 days for MCG+08-11-011 and NGC\,2617, respectively, combined with roughly daily X-ray and near-UV data from Swift for NGC 2617. We find lags consistent with geometrically thin accretion-disk models that predict a lag-wavelength relation of $\tau \propto \lambda^{4/3}$. However, the observed lags are larger than predictions based on standard thin-disk theory by factors of 3.3 for MCG+08-11-011 and 2.3 for NGC\,2617. These differences can be explained if the mass accretion rates are larger than inferred from the optical luminosity by a factor of 4.3 in MCG+08-11-011 and a factor of 1.3 in NGC\,2617, although uncertainty in the SMBH masses determines the significance of this result. While the X-ray variability in NGC\,2617 precedes the UV/optical variability, the long 2.6 day lag is problematic for coronal reprocessing models.

We present the first results from an optical reverberation mapping campaign executed in 2014, targeting the active galactic nuclei (AGN) MCG+08-11-011, NGC 2617, NGC 4051, 3C 382, and Mrk 374. Our targets have diverse and interesting observational properties, including a "changing look" AGN and a broad-line radio galaxy. Based on continuum-H$\beta$ lags, we measure black hole masses for all five targets. We also obtain H$\gamma$ and He{\sc ii}\,$\lambda 4686$ lags for all objects except 3C 382. The He{\sc ii}\,$\lambda 4686$ lags indicate radial stratification of the BLR, and the masses derived from different emission lines are in general agreement. The relative responsivities of these lines are also in qualitative agreement with photoionization models. These spectra have extremely high signal-to-noise ratios (100--300 per pixel) and there are excellent prospects for obtaining velocity-resolved reverberation signatures.

Hot, Dust-Obscured Galaxies, or "Hot DOGs", are a rare, dusty, hyperluminous galaxy population discovered by the WISE mission. Predominantly at redshifts 2-3, they include the most luminous known galaxies in the universe. Their high luminosities likely come from accretion onto highly obscured super massive black holes (SMBHs). We have conducted a pilot survey to measure the SMBH masses of five z~2 Hot DOGs via broad H_alpha emission lines, using Keck/MOSFIRE and Gemini/FLAMINGOS-2. We detect broad H_alpha emission in all five Hot DOGs. We find substantial corresponding SMBH masses for these Hot DOGs (~ 10^{9} M_sun), and their derived Eddington ratios are close to unity. These z~2 Hot DOGs are the most luminous AGNs at given BH masses, suggesting they are accreting at the maximum rates for their BHs. A similar property is found for known z~6 quasars. Our results are consistent with scenarios in which Hot DOGs represent a transitional, high-accretion phase between obscured and unobscured quasars. Hot DOGs may mark a special evolutionary stage before the red quasar and optical quasar phases, and they may be present at other cosmic epochs.

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.

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.

Broad (~10,000 km/s), double-peaked emission-line profiles of Balmer lines emitted by active galactic nuclei (AGN) are thought to originate in the outer parts of an accretion disk surrounding a nuclear supermassive black hole (SMBH), at R~1000 gravitational radii and are most frequently observed in the nuclear spectra of low-luminosity AGN (LLAGN) and radio-galaxies. In the present paper we argue that broad double-peaked profiles are present also in the spectra of other Type 1 AGN, such as Seyfert 1 galaxies, suggesting that the inner part of the broad-line region (BLR) is also the outer part of the accretion disk. We use the Palomar spectral survey of nearby galaxies to show that the only difference between Seyfert 1 BLR line profiles and "bona fide" double peakers is that, in most cases, besides a disk component, we need an additional Gaussian component attributed to non-disk clouds. The recognition that the inner and most variable part of the BLR has a disk geometry suggests that the factor 'f ' in the expression to obtain the SMBH mass (M) in Type 1 AGN M(SMBH)=f [R(BLR) (Delta V)^2]/G is f = 1/[sin^2(i)] for the disk dominated sources. Our median inclination i=27 degrees implies f=4.5, very close to the most recent value of f=4.3+/-1.05, obtained from independent studies. We derive a relation between f and the FWHM of the broad profile that may help to reduce the uncertainties in the SMBH mass determinations of AGN.

Quasar emission lines are often shifted from the systemic velocity due to various dynamical and radiative processes in the line-emitting region. The level of these velocity shifts depends both on the line species and on quasar properties. We study velocity shifts for the line peaks of various narrow and broad quasar emission lines relative to systemic using a sample of 849 quasars from the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project. The coadded (from 32 epochs) spectra of individual quasars have sufficient signal-to-noise ratio (SNR) to measure stellar absorption lines to provide reliable systemic velocity estimates, as well as weak narrow emission lines. The sample also covers a large dynamic range in quasar luminosity (~2 dex), allowing us to explore potential luminosity dependence of the velocity shifts. We derive average line peak velocity shifts as a function of quasar luminosity for different lines, and quantify their intrinsic scatter. We further quantify how well the peak velocity can be measured for various lines as a function of continuum SNR, and demonstrate there is no systematic bias in the line peak measurements when the spectral quality is degraded to as low as SNR~3 per SDSS pixel. Based on the observed line shifts, we provide empirical guidelines on redshift estimation from [OII]3728, [OIII]5008, [NeV]3426, MgII, CIII], HeII1640, broad Hbeta, CIV, and SiIV, which are calibrated to provide unbiased systemic redshifts in the mean, but with increasing intrinsic uncertainties of 46, 56, 119, 205, 233, 242, 400, 415, and 477 km/s, in addition to the measurement uncertainties. These more realistic redshift uncertainties are generally much larger than the formal uncertainties reported by the redshift pipelines for spectroscopic quasar surveys, and demonstrate the infeasibility of measuring quasar redshifts to better than ~200 km/s with only broad lines.

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 XQ-100 survey provides optical and near infrared coverage of 36 blindly selected, intervening damped Lyman alpha systems (DLAs) at 2 < z < 4, simultaneously covering the MgII doublet at 2796A, 2803A, and the Ly-alpha transition. Using the XQ-100 DLA sample, we investigate the completeness of selecting DLA absorbers based on their MgII rest-frame equivalent width (W2796) at these redshifts. Of the 29 DLAs with clean MgII profiles, we find that six (20% of DLAs) have W2796 < 0.6A. The DLA incidence rate of W2796 < 0.6A absorbers is a factor of ~5 higher than what is seen in z~1 samples, indicating a potential evolution in the MgII properties of DLAs with redshift. All of the W2796 < 0.6A DLAs have low metallicities (-2.5 < [M/H] < -1.7), small velocity widths (v90 < 50 km/s), and tend to have relatively low N(HI). We demonstrate that the exclusion of these low W2796 DLAs results in a higher mean N(HI) which in turn leads to a ~7% increase in the cosmological gas density of HI of DLAs at 2 < z < 4; and that this exclusion has a minimal effect on the HI-weighted mean metallicity.

The XQ-100 survey has provided high signal-noise spectra of 100 redshift 3-4.5 quasars with the X-Shooter spectrograph. The metal abundances for 13 elements in the 41 damped Lyman alpha systems (DLAs) identified in the XQ-100 sample are presented, and an investigation into abundances of a variety of DLA classes is conducted. The XQ-100 DLA sample contains five DLAs within 5000 km/s of their host quasar (proximate DLAs; PDLAs) as well as three sightlines which contain two DLAs within 10,000 km/s of each other along the same line-of-sight (multiple DLAs; MDLAs). Combined with previous observations in the literature, we demonstrate that PDLAs with logN(HI)<21.0 show lower [S/H] and [Fe/H] (relative to intervening systems with similar redshift and N(HI)), whilst higher [S/H] and [Si/H] are seen in PDLAs with logN(HI)>21.0. These abundance discrepancies are independent of their line-of-sight velocity separation from the host quasar, and the velocity width of the metal lines (v90). Contrary to previous studies, MDLAs show no difference in [alpha/Fe] relative to single DLAs matched in metallicity and redshift. In addition, we present follow-up UVES data of J0034+1639, a sightline containing three DLAs, including a metal-poor DLA with [Fe/H]=-2.82 (the third lowest [Fe/H] in DLAs identified to date) at z=4.25. Lastly we study the dust-corrected [Zn/Fe], emphasizing that near-IR coverage of X-Shooter provides unprecedented access to MgII, CaII and TiII lines (at redshifts 3-4) to provide additional evidence for subsolar [Zn/Fe] ratio in DLAs.

Red quasars are candidate young objects in an early transition stage of massive galaxy evolution. Our team recently discovered a population of extremely red quasars (ERQs) in the Baryon Oscillation Spectroscopic Survey (BOSS) that has a suite of peculiar emission-line properties including large rest equivalent widths (REWs), unusual "wingless" line profiles, large NV/Lya, NV/CIV, SiIV/CIV and other flux ratios, and very broad and blueshifted [OIII] 5007. Here we present a new catalog of CIV and NV emission-line data for 216,188 BOSS quasars to characterize the ERQ line properties further. We show that they depend sharply on UV-to-mid-IR color, secondarily on REW(CIV), and not at all on luminosity or the Baldwin Effect. We identify a "core" sample of 97 ERQs with nearly uniform peculiar properties selected via i-W3 > 4.6 (AB) and REW(CIV) > 100 A at redshifts 2.0-3.4. A broader search finds 235 more red quasars with similar unusual characteristics. The core ERQs have median luminosity log L (ergs/s) ~ 47.1, sky density 0.010 deg^-2, surprisingly flat/blue UV spectra given their red UV-to-mid-IR colors, and common outflow signatures including BALs or BAL-like features and large CIV emission-line blueshifts. Their SEDs and line properties are inconsistent with normal quasars behind a dust reddening screen. We argue that the core ERQs are a unique obscured quasar population with extreme physical conditions related to powerful outflows across the line-forming regions. Patchy obscuration by small dusty clouds could produce the observed UV extinctions without substantial UV reddening.

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 use the coadded spectra of 32 epochs of Sloan Digital Sky Survey (SDSS) Reverberation Mapping Project observations of 482 quasars with z>1.46 to highlight systematic biases in the SDSS- and BOSS-pipeline redshifts due to the natural diversity of quasar properties. We investigate the characteristics of this bias by comparing the BOSS-pipeline redshifts to an estimate from the centroid of HeII 1640. HeII has a low equivalent width but is often well-defined in high-S/N spectra, does not suffer from self-absorption, and has a narrow component that, when present (the case for about half of our sources), produces a redshift estimate that, on average, is consistent with that determined from [OII] to within 1-sigma of the quadrature sum of the HeII and [OII] centroid measurement uncertainties. The large redshift differences of ~1000 km/s, on average, between the BOSS-pipeline and HeII-centroid redshifts suggest there are significant biases in a portion of BOSS quasar redshift measurements. Adopting the HeII-based redshifts shows that CIV does not exhibit a ubiquitous blueshift for all quasars, given the precision probed by our measurements. Instead, we find a distribution of CIV centroid blueshifts across our sample, with a dynamic range that (i) is wider than that previously reported for this line, and (ii) spans CIV centroids from those consistent with the systemic redshift to those with significant blueshifts of thousands of kilometers per second. These results have significant implications for measurement and use of high-redshift quasar properties and redshifts and studies based thereon.

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.

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 the dependence on data quality of quasar properties measured from the CIV emission line region at high redshifts. Our measurements come from 32 epochs of Sloan Digital Sky Survey (SDSS) Reverberation Mapping Project spectroscopic observations of 482 z>1.46 quasars. We compare the differences between measurements made from the single-epoch and coadded spectra, focusing on the CIV emission line because of its importance for studies of high-redshift quasar demographics and physical properties, including black hole masses. In addition to statistical errors increasing (by factors of ~2-4), we find increasing systematic offsets with decreasing S/N. The systematic difference (measurement uncertainty) in our lowest S/N (<5) subsample between the single-epoch and coadded spectrum (i) CIV equivalent width is 17A (31A), (ii) centroid wavelength is <1A (2A), and fractional velocity widths, Delta V/V, characterized by (iii) the line dispersion, sigma_l, is 0.104 (0.12), and (iv) the mean absolute deviation (MAD) is 0.072 (0.11). These remain smaller than the 1-sigma measurement uncertainties for all subsamples considered. The MAD is found to be the most robust line-width characterization. Offsets in the CIV full-width at half maximum (FWHM) velocity width and the CIV profile characterized by FWHM/sigma_l are only smaller than the statistical uncertainties when S/N>10, although offsets in lower S/N spectra exceed the statistical uncertainties by only a factor of ~1.5. Characterizing the CIV line profile by the kurtosis is the least robust property investigated, as the median systematic coadded--single-epoch measurement differences are larger than the statistical uncertainties for all S/N subsamples.

Reverberation mapping (RM) measurements of broad-line region (BLR) lags in z>0.3 quasars are important for directly measuring black hole masses in these distant objects, but so far there have been limited attempts and success given the practical difficulties of RM in this regime. Here we report preliminary results of 15 BLR lag measurements from the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project, a dedicated RM program with multi-object spectroscopy designed for RM over a wide redshift range. The lags are based on the 2014 spectroscopic light curves alone (32 epochs over 6 months) and focus on the Hbeta and MgII broad lines in the 100 lowest-redshift (z<0.8) quasars included in SDSS-RM; they represent a small subset of the lags that SDSS-RM (including 849 quasars to z~4.5) is expected to deliver. The reported preliminary lag measurements are for intermediate-luminosity quasars at 0.3<~z<0.8, including 9 Hbeta lags and 6 MgII lags, for the first time extending RM results to this redshift-luminosity regime and providing direct quasar black hole mass estimates over ~ half of cosmic time. The MgII lags also increase the number of known MgII lags by several-fold, and start to explore the utility of MgII for RM at high redshift. The location of these new lags at higher redshifts on the observed BLR size-luminosity relationship is statistically consistent with previous Hbeta results at z<0.3. However, an independent constraint on the relationship slope at z>0.3 is not yet possible due to the limitations in our current sample. Our results demonstrate the general feasibility and potential of multi-object RM for z>0.3 quasars.

We present a sample of 38 intervening Damped Lyman $\alpha$ (DLA) systems identified towards 100 $z>3.5$ quasars, observed during the XQ-100 survey. The XQ-100 DLA sample is combined with major DLA surveys in the literature. The final combined sample consists of 742 DLAs over a redshift range approximately $1.6 < z_{\rm abs} < 5.0$. We develop a novel technique for computing $\Omega_{\rm HI}^{\rm DLA}$ as a continuous function of redshift, and we thoroughly assess and quantify the sources of error therein, including fitting errors and incomplete sampling of the high column density end of the column density distribution function. There is a statistically significant redshift evolution in $\Omega_{\rm HI}^{\rm DLA}$ ($\geq 3 \sigma$) from $z \sim 2$ to $z \sim$ 5. In order to make a complete assessment of the redshift evolution of $\Omega_{\rm HI}$, we combine our high redshift DLA sample with absorption surveys at intermediate redshift and 21cm emission line surveys of the local universe. Although $\Omega_{\rm HI}^{\rm DLA}$, and hence its redshift evolution, remains uncertain in the intermediate redshift regime ($0.1 < z_{\rm abs} < 1.6$), we find that the combination of high redshift data with 21cm surveys of the local universe all yield a statistically significant evolution in $\Omega_{\rm HI}$ from $z \sim 0$ to $z \sim 5$ ($\geq 3 \sigma$). Despite its statistical significance, the magnitude of the evolution is small: a linear regression fit between $\Omega_{\rm HI}$ and $z$ yields a typical slope of $\sim$0.17$\times 10^{-3}$, corresponding to a factor of $\sim$ 4 decrease in $\Omega_{\rm HI}$ between $z=5$ and $z=0$.

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.

We explore the variability of quasars in the MgII and Hbeta broad emission lines and UV/optical continuum emission using the Sloan Digital Sky Survey Reverberation Mapping project (SDSS-RM). This is the largest spectroscopic study of quasar variability to date: our study includes 29 spectroscopic epochs from SDSS-RM over $6$ months, containing 357 quasars with MgII and 41 quasars with Hbeta . On longer timescales, the study is also supplemented with two-epoch data from SDSS-I/II. The SDSS-I/II data include an additional $2854$ quasars with MgII and 572 quasars with Hbeta. The MgII emission line is significantly variable ($\Delta f/f$ 10% on 100-day timescales), a necessary prerequisite for its use for reverberation mapping studies. The data also confirm that continuum variability increases with timescale and decreases with luminosity, and the continuum light curves are consistent with a damped random-walk model on rest-frame timescales of $\gtrsim 5$ days. We compare the emission-line and continuum variability to investigate the structure of the broad-line region. Broad-line variability shows a shallower increase with timescale compared to the continuum emission, demonstrating that the broad-line transfer function is not a $\delta$-function. Hbeta is more variable than MgII (roughly by a factor of $1.5$), suggesting different excitation mechanisms, optical depths and/or geometrical configuration for each emission line. The ensemble spectroscopic variability measurements enabled by the SDSS-RM project have important consequences for future studies of reverberation mapping and black hole mass estimation of $1<z<2$ quasars.

We report the discovery of rapid variations of a high-velocity CIV broad absorption line trough in the quasar SDSS J141007.74+541203.3. This object was intensively observed in 2014 as a part of the Sloan Digital Sky Survey Reverberation Mapping Project, during which 32 epochs of spectroscopy were obtained with the Baryon Oscillation Spectroscopic Survey spectrograph. We observe significant (>4sigma) variability in the equivalent width of the broad (~4000 km/s wide) CIV trough on rest-frame timescales as short as 1.20 days (~29 hours), the shortest broad absorption line variability timescale yet reported. The equivalent width varied by ~10% on these short timescales, and by about a factor of two over the duration of the campaign. We evaluate several potential causes of the variability, concluding that the most likely cause is a rapid response to changes in the incident ionizing continuum. If the outflow is at a radius where the recombination rate is higher than the ionization rate, the timescale of variability places a lower limit on the density of the absorbing gas of n_e > 3.9 x 10^5 cm^-3. The broad absorption line variability characteristics of this quasar are consistent with those observed in previous studies of quasars, indicating that such short-term variability may in fact be common and thus can be used to learn about outflow characteristics and contributions to quasar/host-galaxy feedback scenarios.

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.

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

We 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 multi-wavelength observations that trace more than 40 years in the life of the active galactic nucleus (AGN) in Mrk 590, traditionally known as a classic Seyfert 1 galaxy. From spectra recently obtained from HST, Chandra, and the Large Binocular Telescope, we find that the activity in the nucleus of Mrk 590 has diminished so significantly that the continuum luminosity is a factor of 100 lower than the peak luminosity probed by our long baseline observations. Furthermore, the broad emission lines, once prominent in the UV/optical spectrum, have all but disappeared. Since AGN type is defined by the presence of broad emission lines in the optical spectrum, our observations demonstrate that Mrk 590 has now become a "changing look" AGN. If classified by recent optical spectra, Mrk 590 would be a Seyfert ~1.9-2, where the only broad emission line still visible in the optical spectrum is a weak component of Halpha. As an additional consequence of this change, we have definitively detected UV narrow-line components in a Type 1 AGN, allowing an analysis of these emission-line components with high-resolution COS spectra. These observations challenge the historical paradigm that AGN type is only a consequence of the line of sight viewing angle toward the nucleus in the presence of a geometrically-flattened, obscuring medium (i.e., the torus). Our data instead suggest that the current state of Mrk 590 is a consequence of the change in luminosity, which implies the black hole accretion rate has significantly decreased.

A large reverberation mapping study of the Seyfert 1 galaxy NGC 7469 has yielded emission-line lags for Hbeta 4861 and He II 4686 and a central black hole mass measurement of about 10 million solar masses, consistent with previous measurements. A very low level of variability during the monitoring campaign precluded meeting our original goal of recovering velocity-delay maps from the data, but with the new Hbeta measurement, NGC 7469 is no longer an outlier in the relationship between the size of the Hbeta-emitting broad-line region and the AGN luminosity. It was necessary to detrend the continuum and Hbeta and He II 4686 line light curves and those from archival UV data for different time-series analysis methods to yield consistent results.

The Sloan Digital Sky Survey Reverberation Mapping project (SDSS-RM) is a dedicated multi-object RM experiment that has spectroscopically monitored a sample of 849 broad-line quasars in a single 7 deg$^2$ field with the SDSS-III BOSS spectrograph. The RM quasar sample is flux-limited to i_psf=21.7 mag, and covers a redshift range of 0.1<z<4.5. Optical spectroscopy was performed during 2014 Jan-Jul dark/grey time, with an average cadence of ~4 days, totaling more than 30 epochs. Supporting photometric monitoring in the g and i bands was conducted at multiple facilities including the CFHT and the Steward Observatory Bok telescopes in 2014, with a cadence of ~2 days and covering all lunar phases. The RM field (RA, DEC=14:14:49.00, +53:05:00.0) lies within the CFHT-LS W3 field, and coincides with the Pan-STARRS 1 (PS1) Medium Deep Field MD07, with three prior years of multi-band PS1 light curves. The SDSS-RM 6-month baseline program aims to detect time lags between the quasar continuum and broad line region (BLR) variability on timescales of up to several months (in the observed frame) for ~10% of the sample, and to anchor the time baseline for continued monitoring in the future to detect lags on longer timescales and at higher redshift. SDSS-RM is the first major program to systematically explore the potential of RM for broad-line quasars at z>0.3, and will investigate the prospects of RM with all major broad lines covered in optical spectroscopy. SDSS-RM will provide guidance on future multi-object RM campaigns on larger scales, and is aiming to deliver more than tens of BLR lag detections for a homogeneous sample of quasars. We describe the motivation, design and implementation of this program, and outline the science impact expected from the resulting data for RM and general quasar science.

This document presents the results from the Distances subgroup of the Cosmic Frontier Community Planning Study (Snowmass 2013). We summarize the current state of the field as well as future prospects and challenges. In addition to the established probes using Type IA supernovae and baryon acoustic oscillations, we also consider prospective methods based on clusters, active galactic nuclei, gravitational wave sirens and strong lensing time delays.

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.

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.

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.

We present the single-epoch black hole mass estimators based on the CIV (1549 A) broad emission line, using the updated sample of the reverberation-mapped AGNs and high-quality UV spectra. By performing multi-component spectral fitting analysis, we measure the CIV line widths (FWHM_CIV) and line dispersion (sigma_CIV) and the continuum luminosity at 1350 A (L_1350) to calibrate the CIV-based mass estimators. By comparing with the Hbeta reverberation-based masses, we provide new mass estimators with the best-fit relationships, i.e., M_BH \propto L_1350 ^ (0.50+-0.07) sigma_CIV ^2 and M_BH \propto L_1350 ^ (0.52+-0.09) FWHM_CIV ^ (0.56+-0.48). The new CIV-based mass estimators show significant mass-dependent systematic difference compared to the estimators commonly used in the literature. Using the published Sloan Digital Sky Survey QSO catalog, we show that the black hole mass of high-redshift QSOs decreases on average by ~0.25 dex if our recipe is adopted.

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.

We present velocity-resolved reverberation results for five active galactic nuclei. We recovered velocity-delay maps using the maximum-entropy method for four objects: Mrk 335, Mrk 1501, 3C120, and PG2130+099. For the fifth, Mrk 6, we were only able to measure mean time delays in different velocity bins of the Hbeta emission line. The four velocity-delay maps show unique dynamical signatures for each object. For 3C120, the Balmer lines show kinematic signatures consistent with both an inclined disk and infalling gas, but the HeII 4686 emission line is suggestive only of inflow. The Balmer lines in Mrk 335, Mrk 1501, and PG 2130+099 show signs of infalling gas, but the HeII emission in Mrk 335 is consistent with an inclined disk. We also see tentative evidence of combined virial motion and infalling gas from the velocity-binned analysis of Mrk 6. The maps for 3C120 and Mrk 335 are two of the most clearly defined velocity-delay maps to date. These maps constitute a large increase in the number of objects for which we have resolved velocity-delay maps and provide evidence supporting the reliability of reverberation-based black hole mass measurements.

We use a combination of reverberation mapping data and single-epoch spectra of the CIV emission line in a sample of both low and high-redshift active galactic nuclei (AGNs) to investigate sources of the discrepancies between CIV- and Hbeta-based single-epoch black hole mass estimates. We find that for all reverberation mapped sources, there is a component of the line profile that does not reverberate, and the velocity characteristics of this component vary from object-to-object. The differing strength and properties of this non-variable component are responsible for much of the scatter in CIV-based black hole masses compared to Hbeta masses. The CIV mass bias introduced by this non-variable component is correlated with the shape of the CIV line, allowing us to make an empirical correction to the black hole mass estimates. Using this correction and accounting for other sources of scatter such as poor data quality and data inhomogeneity reduces the scatter between the CIV and Hbeta masses in our sample by a factor of ~2, to only ~0.2 dex. We discuss the possibility that this non-variable CIV component originates in an orientation-dependent outflow from either the proposed broad line region (BLR) disk-wind or the intermediate line region (ILR), a high-velocity inner extension of the narrow line region (NLR).

We present the results from a detailed analysis of photometric and spectrophotometric data on five Seyfert 1 galaxies observed as a part of a recent reverberation mapping program. The data were collected at several observatories over a 140-day span beginning in 2010 August and ending in 2011 January. We obtained high sampling-rate light curves for Mrk 335, Mrk 1501, 3C120, Mrk 6, and PG2130+099, from which we have measured the time lag between variations in the 5100 Angstrom continuum and the H-beta broad emission line. We then used these measurements to calculate the mass of the supermassive black hole at the center of each of these galaxies. Our new measurements substantially improve previous measurements of MBH and the size of the broad line-emitting region for four sources and add a measurement for one new object. Our new measurements are consistent with photoionization physics regulating the location of the broad line region in active galactic nuclei.

In 2008 AGILE and Fermi detected gamma-ray flaring activity from the unidentified EGRET source 3EG J1236+0457, recently associated with a flat spectrum radio quasar GB6 J1239+0443 at z=1.762. The optical counterpart of the gamma-ray source underwent a flux enhancement of a factor 15-30 in 6 years, and of ~10 in six months. We interpret this flare-up in terms of a transition from an accretion-disk dominated emission to a synchrotron-jet dominated one. We analysed a Sloan Digital Sky Survey (SDSS) archival optical spectrum taken during a period of low radio and optical activity of the source. We estimated the mass of the central black hole using the width of the CIV emission line. In our work, we have also investigated SDSS archival optical photometric data and UV GALEX observations to estimate the thermal-disk emission contribution of GB6 J1239+0443. Our analysis of the gamma-ray data taken during the flaring episodes indicates a flat gamma-ray spectrum, with an extension of up to 15 GeV, with no statistically-relevant sign of absorption from the broad line region, suggesting that the blazar-zone is located beyond the broad line region. This result is confirmed by the modeling of the broad-band spectral energy distribution (well constrained by the available multiwavelength data) of the flaring activity periods and by the accretion disk luminosity and black hole mass estimated by us using archival data.

Estimates of the mass of super-massive black holes (BHs) in distant active galactic nuclei (AGNs) can be obtained efficiently only through single-epoch spectra, using a combination of their broad emission-line widths and continuum luminosities. Yet the reliability and accuracy of the method, and the resulting mass estimates, M_BH, remain uncertain. A recent study by Croom using a sample of SDSS, 2QZ and 2SLAQ quasars suggests that line widths contribute little information about the BH mass in these single-epoch estimates and can be replaced by a constant value without significant loss of accuracy. In this Letter, we use a sample of nearby reverberation-mapped AGNs to show that this conclusion is not universally applicable. We use the bulge luminosity (L_Bulge) of these local objects to test how well the known M_BH - L_Bulge correlation is recovered when using randomly assigned line widths instead of the measured ones to estimate M_BH. We find that line widths provide significant information about M_BH, and that for this sample, the line width information is just as significant as that provided by the continuum luminosities. We discuss the effects of observational biases upon the analysis of Croom and suggest that the results can probably be explained as a bias of flux-limited, shallow quasar samples.

We present the first results from a detailed analysis of photometric and spectrophotometric data on the narrow-line Seyfert 1 galaxy Mrk 335, collected over a 120-day span in the fall of 2010. From these data we measure the lag in the He II 4686 broad emission line relative to the optical continuum to be 2.7 \pm 0.6 days and the lag in the H\beta 4861 broad emission line to be 13.9 \pm 0.9 days. Combined with the line width, the He II lag yields a black hole mass, MBH = (2.6 \pm 0.8)\times 10^7 Msun. This measurement is consistent with measurements made using the H\beta 4861 line, suggesting that the He II emission originates in the same region as H\beta, but at a much smaller radius. This constitutes the first robust lag measurement for a high-ionization line in a narrow-line Seyfert 1 galaxy.

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.

Using a sample of high-redshift lensed quasars from the CASTLES project with observed-frame ultraviolet or optical and near-infrared spectra, we have searched for possible biases between supermassive black hole (BH) mass estimates based on the CIV, Halpha and Hbeta broad emission lines. Our sample is based upon that of Greene, Peng & Ludwig, expanded with new near-IR spectroscopic observations, consistently analyzed high S/N optical spectra, and consistent continuum luminosity estimates at 5100A. We find that BH mass estimates based on the FWHM of CIV show a systematic offset with respect to those obtained from the line dispersion, sigma_l, of the same emission line, but not with those obtained from the FWHM of Halpha and Hbeta. The magnitude of the offset depends on the treatment of the HeII and FeII emission blended with CIV, but there is little scatter for any fixed measurement prescription. While we otherwise find no systematic offsets between CIV and Balmer line mass estimates, we do find that the residuals between them are strongly correlated with the ratio of the UV and optical continuum luminosities. Removing this dependency reduces the scatter between the UV- and optical-based BH mass estimates by a factor of approximately 2, from roughly 0.35 to 0.18 dex. The dispersion is smallest when comparing the CIV sigma_l mass estimate, after removing the offset from the FWHM estimates, and either Balmer line mass estimate. The correlation with the continuum slope is likely due to a combination of reddening, host contamination and object-dependent SED shapes. When we add additional heterogeneous measurements from the literature, the results are unchanged.

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 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 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 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 present the results of a recent reverberation-mapping campaign undertaken to improve measurements of the radius of the broad line region and the central black hole mass of the quasar PG 2130+099. Cross correlation of the 5100 angstrom continuum and H-beta emission-line light curves yields a time lag of 22.9 (+4.4 - 4.3) days, corresponding to a central black hole mass MBH= 3.8 (+/- 1.5) x 10^7 Msun. This value supports the notion that previous measurements yielded an incorrect lag. We re-analyzed previous datasets to investigate the possible sources of the discrepancy and conclude that previous measurement errors were apparently caused by a combination of undersampling of the light curves and long-term secular changes in the H-beta emission-line equivalent width. With our new measurements, PG 2130+099 is no longer an outlier in either the R-L or the MBH-Sigma relationships.

We have measured the host-galaxy starlight contribution to four lower-luminosity AGNs (NGC 3516, NGC 4593, IC 4329A, and NGC 7469). We include these objects with new broad line region measurements for NGC 4151 and NGC 4593 to present a revised version of the radius-luminosity relationship for AGNs.

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