The Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) program monitors 849 active galactic nuclei (AGN) both spectroscopically and photometrically. The photometric observations used in this work span over four years and provide an excellent baseline for variability studies of these objects. We present the photometric light curves from 2014 to 2017 obtained by the Steward Observatory's Bok telescope and the CFHT telescope with MegaCam. We provide details on the data acquisition and processing of the data from each telescope, the difference imaging photometry used to produce the light curves, and the calculation of a variability index to quantify each AGN's variability. We find that the Welch-Stetson J-index provides a useful characterization of AGN variability and can be used to select AGNs for further study.

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.

EarthFinder is a NASA Astrophysics Probe mission concept selected for study as input to the 2020 Astrophysics National Academies Decadal Survey. The EarthFinder concept is based on a dramatic shift in our understanding of how PRV measurements should be made. We propose a new paradigm which brings the high precision, high cadence domain of transit photometry as demonstrated by Kepler and TESS to the challenges of PRV measurements at the cm/s level. This new paradigm takes advantage of: 1) broad wavelength coverage from the UV to NIR which is only possible from space to minimize the effects of stellar activity; 2) extremely compact, highly stable, highly efficient spectrometers (R>150,000) which require the diffraction-limited imaging possible only from space over a broad wavelength range; 3) the revolution in laser-based wavelength standards to ensure cm/s precision over many years; 4) a high cadence observing program which minimizes sampling-induced period aliases; 5) exploiting the absolute flux stability from space for continuum normalization for unprecedented line-by-line analysis not possible from the ground; and 6) focusing on the bright stars which will be the targets of future imaging missions so that EarthFinder can use a ~1.5 m telescope.

Quasar broad emission lines are largely powered by photoionization from the accretion continuum. Increased central luminosity will enhance line emissivity in more distant clouds, leading to increased average distance of the broad-line-emitting clouds and decreased averaged line width, known as the broad-line region (BLR) "breathing". However, different lines breathe differently, and some high-ionization lines, such as C IV, can even show "anti-breathing" where the line broadens when luminosity increases. Using multi-year photometric and spectroscopic monitoring data from the Sloan Digital Sky Survey Reverberation Mapping project, we quantify the breathing effect ($\Delta$log W=$\alpha\Delta$log L) of broad H$\alpha$, H$\beta$, Mg II, C IV, and C III] for statistical quasar samples over $z\approx 0.1-2.5$. We found that H$\beta$ displays the most consistent normal breathing expected from the virial relation ($\alpha\sim-0.25$), Mg II and H$\alpha$ on average show no breathing ($\alpha\sim 0$), and C IV (and similarly C III] and Si IV mostly shows anti-breathing ($\alpha>0$). The anti-breathing of C IV can be well understood by the presence of a non-varying core component in addition to a reverberating broad-base component, consistent with earlier findings. The deviation from canonical breathing introduces extra scatter (a luminosity-dependent bias) in single-epoch virial BH mass estimates due to intrinsic quasar variability, which underlies the long argued caveats of C IV single-epoch masses. Using the line dispersion instead of FWHM leads to less, albeit still substantial, deviations from canonical breathing in most cases. Our results strengthen the need for reverberation mapping to provide reliable quasar BH masses, and quantify the level of variability-induced bias in single-epoch BH masses based on various lines.

We present the results of a pilot Hubble Space Telescope (HST) imaging study of the host galaxies of ten quasars from the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project. Probing more than an order of magnitude in BH and stellar masses, our sample is the first statistical sample to study the BH-host correlations beyond z>0.3 with reliable BH masses from reverberation mapping rather than from single-epoch spectroscopy. We perform image decomposition in two HST bands (UVIS-F606W and IR-F110W) to measure host colors and estimate stellar masses using empirical relations between broad-band colors and the mass-to-light ratio. The stellar masses of our targets are mostly dominated by a bulge component. The BH masses and stellar masses of our sample broadly follow the same correlations found for local RM AGN and quiescent bulge-dominant galaxies, with no strong evidence of evolution in the M_BH-M_*,bulge relation to z~0.6. We further compare the host light fraction from HST imaging decomposition to that estimated from spectral decomposition. We found a good correlation between the host fractions derived with both methods. However, the host fraction derived from spectral decomposition is systematically smaller than that from imaging decomposition by ~30%, indicating different systematics in both approaches. This study paves the way for upcoming more ambitious host galaxy studies of quasars with direct RM-based BH masses at high redshift.

This paper documents the sixteenth data release (DR16) from the Sloan Digital Sky Surveys; the fourth and penultimate from the fourth phase (SDSS-IV). This is the first release of data from the southern hemisphere survey of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2); new data from APOGEE-2 North are also included. DR16 is also notable as the final data release for the main cosmological program of the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), and all raw and reduced spectra from that project are released here. DR16 also includes all the data from the Time Domain Spectroscopic Survey (TDSS) and new data from the SPectroscopic IDentification of ERosita Survey (SPIDERS) programs, both of which were co-observed on eBOSS plates. DR16 has no new data from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey (or the MaNGA Stellar Library "MaStar"). We also preview future SDSS-V operations (due to start in 2020), and summarize plans for the final SDSS-IV data release (DR17).

We present reverberation mapping results for the MgII 2800 A broad emission line in a sample of 193 quasars at 0.35<z<1.7 with photometric and spectroscopic monitoring observations from the Sloan Digital Sky Survey Reverberation Mapping project during 2014 - 2017. We find significant time lags between the MgII and continuum lightcurves for 57 quasars and define a "gold sample" of 24 quasars with the most reliable lag measurements. We estimate false-positive rates for each lag that range from 1-24%, with an average false-positive rate of 11% for the full sample and 8% for the gold sample. There are an additional ~40 quasars with marginal MgII lag detections which may yield reliable lags after additional years of monitoring. The MgII lags follow a radius -- luminosity relation with a best-fit slope that is consistent with alpha=0.5 but with an intrinsic scatter of 0.36dex that is significantly larger than found for the Hb radius -- luminosity relation. For targets with SDSS-RM lag measurements of other emission lines, we find that our MgII lags are similar to the Hb lags and ~2-3 times larger than the CIV lags. This work significantly increases the number of MgII broad-line lags and provides additional reverberation-mapped black hole masses, filling the redshift gap at the peak of supermassive black hole growth between the Hb and CIV emission lines in optical spectroscopy.

Results from a few decades of reverberation mapping (RM) studies have revealed a correlation between the radius of the broad-line emitting region (BLR) and the continuum luminosity of active galactic nuclei. This "radius-luminosity" relation enables survey-scale black-hole mass estimates across cosmic time, using relatively inexpensive single-epoch spectroscopy, rather than intensive RM time monitoring. However, recent results from newer reverberation mapping campaigns challenge this widely used paradigm, reporting quasar BLR sizes that differ significantly from the previously established radius-luminosity relation. Using simulations of the radius--luminosity relation with the observational parameters of SDSS-RM, we find that this difference is not likely due to observational biases. Instead, it appears that previous RM samples were biased to a subset of quasar properties, and the broader parameter space occupied by the SDSS-RM quasar sample has a genuinely wider range of BLR sizes. We examine the correlation between the deviations from the radius-luminosity relation and several quasar parameters; the most significant correlations indicate that the deviations depend on UV/optical SED and the relative amount of ionizing radiation. Our results indicate that single-epoch black-hole mass estimates that do not account for the diversity of quasars in the radius-luminosity relation could be overestimated by an average of ~0.3 dex.

We present an observational study of 22 broad absorption line quasars (BAL QSOs) at $3\lesssim z \lesssim5$ based on optical/near-IR spectroscopy, aiming to investigate quasar winds and their effects. The near-IR spectroscopy covers the \hb\ and/or \mgii\ broad emission lines (BELs) for these quasars, allowing us to estimate their central black hole (BH) masses in a robust way. We found that our BAL QSOs on average do not have a higher Eddington ratio than that from non-BAL QSOs matched in redshift and/or luminosity. In a subset consisting of seven strong BAL QSOs possessing sub-relativistic BAL outflows, we see the prevalence of large \civ-BEL blueshift ($\sim$3100 km s$^{-1}$) and weak \oiii\ emission (particularly the narrow \oiii$\lambda$5007 component), indicative of nuclear outflows affecting the narrow emission-line (NEL) regions. In another subset consisting of thirteen BAL QSOs having simultaneous observations of \mgii\ and \hb, we found a strong correlation between 3000~\AA\ and 5000~\AA\ monochromatic luminosity, consistent with that from non-BAL QSOs matched in redshift and luminosity; however, there is no correlation between \mgii\ and \hb\ in FWHM, likely due to nuclear outflows influencing the BEL regions. Our spectroscopic investigations offer strong evidence that the presence of nuclear outflows plays an important role in shaping the BEL/NEL regions of these quasars and possibly, regulating the growth of central supermassive black holes (SMBHs). We propose that BEL blueshift and BAL could be different manifestations of the same outflow system viewed at different sightlines and/or phases.

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 investigate the performance of different methodologies that measure the time lag between broad-line and continuum variations in reverberation mapping data using simulated light curves that probe a range of cadence, time baseline, and signal-to-noise ratio in the flux measurements. We compare three widely-adopted lag measuring methods: the Interpolated Cross-Correlation Function (ICCF), the z-transformed Discrete Correlation Function (ZDCF) and the MCMC code JAVELIN, for mock data with qualities typical of multi-object spectroscopic reverberation mapping (MOS-RM) surveys that simultaneously monitor hundreds of quasars. We quantify the overall lag detection efficiency, the rate of false detections, and the quality of lag measurements for each of these methods and under different survey designs (e.g., observing cadence and depth) using mock quasar light curves. Overall JAVELIN and ICCF outperform ZDCF in essentially all tests performed. Compared with ICCF, JAVELIN produces higher quality lag measurements, is capable of measuring more lags with timescales shorter than the observing cadence, is less susceptible to seasonal gaps and S/N degradation in the light curves, and produces more accurate lag uncertainties. We measure the Hbeta broad-line region size-luminosity (R-L) relation with each method using the simulated light curves to assess the impact of selection effects of the design of MOS-RM surveys. The slope of the R-L relation measured by JAVELIN is the least biased among the three methods, and is consistent across different survey designs. These results demonstrate a clear preference for JAVELIN over the other two non-parametric methods for MOS-RM programs, particularly in the regime of limited light curve quality as expected from most MOS-RM programs.

We investigate the effects of extended multi-year light curves (9-year photometry and 5-year spectroscopy) on the detection of time lags between the continuum variability and broad-line response of quasars at z>~1.5, and compare with the results using 4-year photometry+spectroscopy presented in a companion paper. We demonstrate the benefits of the extended light curves in three cases: (1) lags that are too long to be detected by the shorter-duration data but can be detected with the extended data; (2) lags that are recovered by the extended light curves but missed in the shorter-duration data due to insufficient light curve quality; and (3) lags for different broad line species in the same object. These examples demonstrate the importance of long-term monitoring for reverberation mapping to detect lags for luminous quasars at high-redshift, and the expected performance of the final dataset from the Sloan Digital Sky Survey Reverberation Mapping project that will have 11-year photometric and 7-year spectroscopic baselines.

We present Searches After Gravitational-waves Using ARizona Observatories (SAGUARO), a comprehensive effort dedicated to the discovery and characterization of optical counterparts to gravitational wave (GW) events. SAGUARO utilizes ground-based facilities ranging from 1.5m to 10m in diameter, located primarily in the Northern Hemisphere. We provide an overview of SAGUARO's telescopic resources, pipeline for transient detection, and database for candidate visualization. We describe SAGUARO's discovery component, which utilizes the $5$~deg$^2$ field-of-view optical imager on the Mt. Lemmon 1.5m telescope, reaching limits of $\approx 21.3$~AB mag while rapidly tiling large areas. We also describe the follow-up component of SAGUARO, used for rapid vetting and monitoring of optical candidates. With the onset of Advanced LIGO/Virgo's third observing run, we present results from the first three SAGUARO searches following the GW events S190408an, S190425z and S190426c, which serve as a valuable proof-of-concept of SAGUARO. We triggered and searched 15, 60 and 60 deg$^{2}$ respectively, 17.6, 1.4 and 41.8 hrs after the initial GW alerts. We covered 7.8, 3.0 and 5.1\% of the total probability within the GW event localizations, reaching 3$\sigma$ limits of 19.8, 21.3 and 20.8 AB mag, respectively. Although no viable counterparts associated with these events were found, we recovered 6 known transients and ruled out 5 potential candidates. We also present Large Binocular Telescope spectroscopy of PS19eq/SN2019ebq, a promising kilonova candidate that was later determined to be a supernova. With the ability to tile large areas and conduct detailed follow-up, SAGUARO represents a significant addition to GW counterpart searches.

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 accretion-disk structure measurements from continuum lags in the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project. Lags are measured using the \texttt{JAVELIN} software from the first-year SDSS-RM $g$ and $i$ photometry, resulting in well-defined lags for 95 quasars, 33 of which have lag SNR $>$ 2$\sigma$. We also estimate lags using the \texttt{ICCF} software and find consistent results, though with larger uncertainties. Accretion-disk structure is fit using a Markov Chain Monte Carlo approach, parameterizing the measured continuum lags as a function of disk size normalization, wavelength, black hole mass, and luminosity. In contrast with previous observations, our best-fit disk sizes and color profiles are consistent (within 1.5~$\sigma$) with the \citet{SS73} analytic solution. We also find that more massive quasars have larger accretion disks, similarly consistent with the analytic accretion-disk model. The data are inconclusive on a correlation between disk size and continuum luminosity, with results that are consistent with both no correlation and with the \citet{SS73} expectation. The continuum lag fits have a large excess dispersion, indicating that our measured lag errors are underestimated and/or our best-fit model may be missing the effects of orientation, spin, and/or radiative efficiency. We demonstrate that fitting disk parameters using only the highest-SNR lag measurements biases best-fit disk sizes to be larger than the disk sizes recovered using a Bayesian approach on the full sample of well-defined lags.

The width of the broad emission lines in quasars is commonly characterized either by the full-width-at-half-maximum (FWHM) or the square root of the second moment of the line profile ($\sigma_{\rm line}$), and used as an indicator of the virial velocity of the broad-line region (BLR) in the estimation of black hole (BH) mass. We measure FWHM and $\sigma_{\rm line}$ for H$\alpha$, H$\beta$ and Mg II broad lines in both the mean and root-mean-square (rms) spectra of a large sample of quasars from the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project. We introduce a new quantitative recipe to measure $\sigma_{\rm line}$ that is reproducible, less susceptible to noise and blending in the wings, and scales with the intrinsic width of the line. We compare the four definitions of line width (FWHM and $\sigma_{\rm line}$ in mean and rms spectra, respectively) for each of the three broad lines and among different lines. There are strong correlations among different width definitions for each line, providing justification for using the line width measured in single-epoch spectroscopy as a virial velocity indicator. There are also strong correlations among different lines, suggesting alternative lines to H$\beta$ can be used to estimate virial BH masses. We further investigate the correlations between virial BH masses using different line width definitions and the stellar velocity dispersion of the host galaxies, and the dependence of line shape (characterized by the ratio FWHM/$\sigma_{\rm line}$) on physical properties of the quasar. Our results provide further evidence that FWHM is more sensitive to the orientation of a flattened BLR geometry than $\sigma_{\rm line}$, but the overall comparison between the virial BH mass and host stellar velocity dispersion does not provide conclusive evidence that one particular width definition is significantly better than the others.

We analyze extensive spectroscopic and photometric data of the hypervariable quasar SDSS J131424+530527 (RMID 017) at z=0.456, an optical "changing look" quasar from the Sloan Digital Sky Survey Reverberation Mapping project that increased in optical luminosity by a factor of 10 between 2014 and 2017. The observed broad emission lines all respond in luminosity and width to the changing optical continuum, as expected for photoionization in a stratified, virialized broad emission line region. The luminosity changes therefore result from intrinsic changes in accretion power rather than variable obscuration. The variability is continuous and apparently stochastic, disfavoring an origin as a discrete event such as a tidal disruption flare or microlensing event. It is coordinated on day timescales with blue leading red, consistent with reprocessing powering the entire optical SED. We show that this process cannot work in a standard thin disk geometry on energetic grounds, and would instead require a large covering factor reprocessor. Disk instability models could potentially also explain the data, provided that the instability sets in near the inner radius of a geometrically thick accretion disk.

We present reverberation-mapping lags and black-hole mass measurements using the CIV 1549 broad emission line from a sample of 349 quasars monitored as a part of the Sloan Digital Sky Survey Reverberation Mapping Project. Our data span four years of spectroscopic and photometric monitoring for a total baseline of 1300 days. We report significant time delays between the continuum and the CIV 1549 emission line in 52 quasars, with an estimated false-positive detection rate of 10%. Our analysis of marginal lag measurements indicates that there are on the order of 100 additional lags that should be recoverable by adding more years of data from the program. We use our measurements to calculate black-hole masses and fit an updated CIV radius-luminosity relationship. Our results significantly increase the sample of quasars with CIV RM results, with the quasars spanning two orders of magnitude in luminosity toward the high-luminosity end of the CIV radius-luminosity relation. In addition, these quasars are located at among the highest redshifts (z~1.4-2.8) of quasars with black hole masses measured with reverberation mapping. This work constitutes the first large sample of CIV reverberation-mapping measurements in more than a dozen quasars, demonstrating the utility of multi-object reverberation mapping campaigns.

We present absorption variability results for 134 bona fide \mgii\ broad absorption line (BAL) quasars at 0.46~$\lesssim z \lesssim$~2.3 covering days to $\sim$ 10 yr in the rest frame. We use multiple-epoch spectra from the Sloan Digital Sky Survey, which has delivered the largest such BAL-variability sample ever studied. \mgii-BAL identifications and related measurements are compiled and presented in a catalog. We find a remarkable time-dependent asymmetry in EW variation from the sample, such that weakening troughs outnumber strengthening troughs, the first report of such a phenomenon in BAL variability. Our investigations of the sample further reveal that (i) the frequency of BAL variability is significantly lower (typically by a factor of 2) than that from high-ionization BALQSO samples; (ii) \mgii\ BAL absorbers tend to have relatively high optical depths and small covering factors along our line of sight; (iii) there is no significant EW-variability correlation between \mgii\ troughs at different velocities in the same quasar; and (iv) the EW-variability correlation between \mgii\ and \aliii\ BALs is significantly stronger than that between \mgii\ and \civ\ BALs at the same velocities. These observational results can be explained by a combined transverse-motion/ionization-change scenario, where transverse motions likely dominate the strengthening BALs while ionization changes and/or other mechanisms dominate the weakening BALs.

We present a detailed characterization of the 849 broad-line quasars from the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project. Our quasar sample covers a redshift range of 0.1<z<4.5 and is flux-limited to i_PSF<21.7 without any other cuts on quasar properties. The main sample characterization includes: 1) spectral measurements of the continuum and broad emission lines for individual objects from the coadded first-season spectroscopy in 2014; 2) identification of broad and narrow absorption lines in the spectra; 3) optical variability properties for continuum and broad lines from multi-epoch spectroscopy. We provide improved systemic redshift estimates for all quasars, and demonstrate the effects of signal-to-noise ratio on the spectral measurements. We compile measured properties for all 849 quasars along with supplemental multi-wavelength data for subsets of our sample from other surveys. The SDSS-RM sample probes a diverse range in quasar properties, and shows well detected continuum and broad-line variability for many objects from first-season monitoring data. The compiled properties serve as the benchmark for follow-up work based on SDSS-RM data. The spectral fitting tools are made public along with this work.

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

This white paper specifies the footprints, cadence requirements, and total-depth requirements needed to allow the most-successful AGN studies in the four currently selected LSST Deep-Drilling Fields (DDFs): ELAIS-S1, XMM-LSS, CDF-S, and COSMOS. The information provided on cadence and total-depth requirements will also likely be applicable to enabling effective AGN science in any additional DDFs that are chosen.

We systematically investigate short-timescale ($<$10-day rest-frame) CIV broad absorption-line (BAL) variability to constrain quasar-wind properties and provide insights into BAL-variability mechanisms in quasars. We employ data taken by the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project, as the rapid cadence of these observations provides a novel opportunity to probe BAL variability on shorter rest-frame timescales than have previously been explored. In a sample of 27 quasars with a median of 58 spectral epochs per quasar, we have identified 15 quasars ($55^{+18}_{-14}$%), 19 of 37 CIV BAL troughs ($51^{+15}_{-12}$%), and 54 of 1460 epoch pairs ($3.7 \pm 0.5$%) that exhibit significant CIV BAL equivalent-width variability on timescales of less than 10 days in the quasar rest frame. These frequencies indicate that such variability is common among quasars and BALs, though somewhat rare among epoch pairs. Thus, models describing BALs and their behavior must account for variability on timescales down to less than a day in the quasar rest frame. We also examine a variety of spectral characteristics and find that in some cases, BAL variability is best described by ionization-state changes, while other cases are more consistent with changes in covering fraction or column density. We adopt a simple model to constrain the density and radial distance of two outflows appearing to vary by ionization-state changes, yielding outflow density lower limits consistent with previous work.

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

We report on an X-ray and optical/UV study of eight Broad Absorption Line (BAL) to non-BAL transforming quasars at $z\,\approx\,$1.7-2.2 over 0.29-4.95 rest-frame years with at least three spectroscopic epochs for each quasar from the SDSS, BOSS, $Gemini$, and ARC 3.5-m telescopes. New $Chandra$ observations obtained for these objects show their values of $\alpha_{\rm ox}$ and $\Delta{\alpha}_{\rm ox}$, as well as their spectral energy distributions, are consistent with those of non-BAL quasars. Moreover, our targets have X-ray spectral shapes that are, on average, consistent with weakened absorption with an effective power-law photon index of $\Gamma_{\rm eff}\,=\,1.69^{+0.25}_{-0.25}$. The newer $Gemini$ and ARC 3.5-m spectra reveal that the BAL troughs have remained absent since the BOSS observations where the BAL disappearance was discovered. The X-ray and optical/UV results in tandem are consistent with at least the X-ray absorbing material moving out of the line-of-sight, leaving an X-ray unabsorbed non-BAL quasar. The UV absorber might have become more highly ionized (in a shielding-gas scenario) or also moved out of the line-of-sight (in a wind-clumping scenario).

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.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) has been in operation since July 2014. This paper describes the second data release from this phase, and the fourteenth from SDSS overall (making this, Data Release Fourteen or DR14). This release makes public data taken by SDSS-IV in its first two years of operation (July 2014-2016). Like all previous SDSS releases, DR14 is cumulative, including the most recent reductions and calibrations of all data taken by SDSS since the first phase began operations in 2000. New in DR14 is the first public release of data from the extended Baryon Oscillation Spectroscopic Survey (eBOSS); the first data from the second phase of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE-2), including stellar parameter estimates from an innovative data driven machine learning algorithm known as "The Cannon"; and almost twice as many data cubes from the Mapping Nearby Galaxies at APO (MaNGA) survey as were in the previous release (N = 2812 in total). This paper describes the location and format of the publicly available data from SDSS-IV surveys. We provide references to the important technical papers describing how these data have been taken (both targeting and observation details) and processed for scientific use. The SDSS website (www.sdss.org) has been updated for this release, and provides links to data downloads, as well as tutorials and examples of data use. SDSS-IV is planning to continue to collect astronomical data until 2020, and will be followed by SDSS-V.

Broad absorption lines (BALs) in the spectra of quasi-stellar objects (QSOs) originate from outflowing winds along our line of sight; winds are thought to originate from the inner regions of the QSO accretion disk, close to the central supermassive black hole (SMBH). Winds likely play a role in galaxy evolution and aid the accretion mechanism onto the SMBH. BAL equivalent widths can change on typical timescales from months to years; such variability is generally attributed to changes in the covering factor and/or in the ionization level of the gas. We investigate BAL variability, focusing on BAL disappearance. We analyze multi-epoch spectra of more than 1500 QSOs -the largest sample ever used for such a study- observed by different programs from the Sloan Digital Sky Survey-I/II/III (SDSS), and search for disappearing C IV BALs. The spectra rest-frame time baseline ranges from 0.28 to 4.9 yr; the source redshifts range from 1.68 to 4.27. We detect 73 disappearing BALs in the spectra of 67 sources. This corresponds to 3.9% of disappearing BALs, and 5.1% of our BAL QSOs exhibit at least one disappearing BAL. We estimate the average lifetime of a BAL along our line of sight (~ 80-100 yr), which appears consistent with the accretion disk orbital time at distances where winds are thought to originate. We inspect properties of the disappearing BALs and compare them to the properties of our main sample. We also investigate the existence of a correlation in the variability of multiple troughs in the same spectrum, and find it persistent at large velocity offsets between BAL pairs, suggesting that a mechanism extending on a global scale is necessary to explain the phenomenon. We select a more reliable sample of disappearing BALs following Filiz Ak et al. (2012), where a subset of our sample was analyzed, and compare the findings from the two works, obtaining generally consistent results.

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.

We present composite broad-line region (BLR) reverberation-mapping lag measurements for \halpha, \hbeta, \HeII\,$\lambda4686$ and \MgII\ for a sample of 144, $z\lesssim 1$ quasars from the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project. Using only the 32-epoch spectroscopic light curves in the first 6-month season of SDSS-RM observations, we compile correlation-function measurements for individual objects and then coadd them to allow the measurement of the average lags for our sample at mean redshifts of $0.4$ (for \halpha) and $\sim 0.65$ (for the other lines). At similar quasar luminosities and redshifts, the sample-averaged lag decreases in the order of \MgII, \halpha, \hbeta\ and \HeII. This decrease in lags is accompanied by an increase in the mean line width of the four lines, and is roughly consistent with the virialized motion for BLR gas in photoionization equilibrium. These are among the first RM measurements of stratified BLR structure at $z>0.3$. Dividing our sample by luminosity, \halpha\ shows clear evidence of increasing lags with luminosity, consistent with the expectation from the measured BLR size-luminosity relation based on \hbeta. The other three lines do not show a clear luminosity trend in their average lags due to the limited dynamic range of luminosity probed and the poor average correlation signals in the divided samples, a situation that will be improved with the incorporation of additional photometric and spectroscopic data from SDSS-RM. We discuss the utility and caveats of composite-lag measurements for large statistical quasar samples with reverberation-mapping data.

As astronomers increasingly exploit the information available in the time domain, spectroscopic variability in particular opens broad new channels of investigation. Here we describe the selection algorithms for all targets intended for repeat spectroscopy in the Time Domain Spectroscopic Survey (TDSS), part of the extended Baryon Oscillation Spectroscopic Survey within the Sloan Digital Sky Survey-IV. Also discussed are the scientific rationale and technical constraints leading to these target selections. The TDSS includes a large "Repeat Quasar Spectroscopy" (RQS) program delivering ~13,000 repeat spectra of confirmed SDSS quasars, and several smaller "Few-Epoch Spectroscopy" (FES) programs targeting specific classes of quasars as well as stars. The RQS program aims to provide a large and diverse quasar data set for studying variations in quasar spectra on timescales of years, a comparison sample for the FES quasar programs, and opportunity for discovering rare, serendipitous events. The FES programs cover a wide variety of phenomena in both quasars and stars. Quasar FES programs target broad absorption line quasars, high signal-to-noise ratio normal broad line quasars, quasars with double-peaked or very asymmetric broad emission line profiles, binary supermassive black hole candidates, and the most photometrically variable quasars. Strongly variable stars are also targeted for repeat spectroscopy, encompassing many types of eclipsing binary systems, and classical pulsators like RR Lyrae. Other stellar FES programs allow spectroscopic variability studies of active ultracool dwarf stars, dwarf carbon stars, and white dwarf/M dwarf spectroscopic binaries. We present example TDSS spectra and describe anticipated sample sizes and results.

We present inferences on the geometry and kinematics of the broad-Hbeta line-emitting region in four active galactic nuclei monitored as a part of the fall 2010 reverberation mapping campaign at MDM Observatory led by the Ohio State University. From modeling the continuum variability and response in emission-line profile changes as a function of time, we infer the geometry of the Hbeta- emitting broad line regions to be thick disks that are close to face-on to the observer with kinematics that are well-described by either elliptical orbits or inflowing gas. We measure the black hole mass to be log (MBH) = 7.25 (+/-0.10) for Mrk 335, 7.86 (+0.20, -0.17) for Mrk 1501, 7.84 (+0.14, -0.19) for 3C 120, and 6.92 (+0.24, -0.23) for PG 2130+099. These black hole mass measurements are not based on a particular assumed value of the virial scale factor f, allowing us to compute individual f factors for each target. Our results nearly double the number of targets that have been modeled in this manner, and investigate the properties of a more diverse sample by including previously modeled objects. We measure an average scale factor f in the entire sample to be log10(f) = 0.54 +/- 0.17 when the line dispersion is used to characterize the line width, which is consistent with values derived using the normalization of the MBH-sigma relation. We find that the scale factor f for individual targets is likely correlated with the black hole mass, inclination angle, and opening angle of the broad line region but we do not find any correlation with the luminosity.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) began observations in July 2014. It pursues three core programs: APOGEE-2, MaNGA, and eBOSS. In addition, eBOSS contains two major subprograms: TDSS and SPIDERS. This paper describes the first data release from SDSS-IV, Data Release 13 (DR13), which contains new data, reanalysis of existing data sets and, like all SDSS data releases, is inclusive of previously released data. DR13 makes publicly available 1390 spatially resolved integral field unit observations of nearby galaxies from MaNGA, the first data released from this survey. It includes new observations from eBOSS, completing SEQUELS. In addition to targeting galaxies and quasars, SEQUELS also targeted variability-selected objects from TDSS and X-ray selected objects from SPIDERS. DR13 includes new reductions of the SDSS-III BOSS data, improving the spectrophotometric calibration and redshift classification. DR13 releases new reductions of the APOGEE-1 data from SDSS-III, with abundances of elements not previously included and improved stellar parameters for dwarf stars and cooler stars. For the SDSS imaging data, DR13 provides new, more robust and precise photometric calibrations. Several value-added catalogs are being released in tandem with DR13, in particular target catalogs relevant for eBOSS, TDSS, and SPIDERS, and an updated red-clump catalog for APOGEE. This paper describes the location and format of the data now publicly available, as well as providing references to the important technical papers that describe the targeting, observing, and data reduction. The SDSS website, http://www.sdss.org, provides links to the data, tutorials and examples of data access, and extensive documentation of the reduction and analysis procedures. DR13 is the first of a scheduled set that will contain new data and analyses from the planned ~6-year operations of SDSS-IV.

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 describe the Sloan Digital Sky Survey IV (SDSS-IV), a project encompassing three major spectroscopic programs. The Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2) is observing hundreds of thousands of Milky Way stars at high resolution and high signal-to-noise ratio in the near-infrared. The Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey is obtaining spatially-resolved spectroscopy for thousands of nearby galaxies (median redshift of z = 0.03). The extended Baryon Oscillation Spectroscopic Survey (eBOSS) is mapping the galaxy, quasar, and neutral gas distributions between redshifts z = 0.6 and 3.5 to constrain cosmology using baryon acoustic oscillations, redshift space distortions, and the shape of the power spectrum. Within eBOSS, we are conducting two major subprograms: the SPectroscopic IDentification of eROSITA Sources (SPIDERS), investigating X-ray AGN and galaxies in X-ray clusters, and the Time Domain Spectroscopic Survey (TDSS), obtaining spectra of variable sources. All programs use the 2.5-meter Sloan Foundation Telescope at Apache Point Observatory; observations there began in Summer 2014. APOGEE-2 also operates a second near-infrared spectrograph at the 2.5-meter du Pont Telescope at Las Campanas Observatory, with observations beginning in early 2017. Observations at both facilities are scheduled to continue through 2020. In keeping with previous SDSS policy, SDSS-IV provides regularly scheduled public data releases; the first one, Data Release 13, was made available in July 2016.

We investigate broad absorption line (BAL) disappearance and emergence using a 470 BAL-quasar sample over < 0.10-5.25 rest-frame years with at least three spectroscopic epochs for each quasar from the Sloan Digital Sky Survey. We identify 14 disappearing BALs over < 1.73-4.62 rest-frame years and 18 emerging BALs over < 1.46-3.66 rest-frame years associated with the CIV 1548,1550 and/or SiIV 1393,1402 doublets, and report on their variability behavior. BAL quasars in our dataset exhibit disappearing/emerging CIV BALs at a rate of 2.3 and 3.0 per cent, respectively, and the frequency for BAL to non-BAL quasar transitions is 1.7 per cent. We detect four re-emerging BALs over < 3.88 rest-frame years on average and three re-disappearing BALs over < 4.15 rest-frame years on average, the first reported cases of these types. We infer BAL lifetimes along the line of sight to be nominally < 100-1000 yr using disappearing CIV BALs in our sample. Interpretations of (re-)emerging and (re-)disappearing BALs reveal evidence that collectively supports both transverse-motion and ionization-change scenarios to explain BAL variations. We constrain a nominal CIV/SiIV BAL outflow location of < 100 pc from the central source and a radial size of > 1x10^-7 pc (0.02 au) using the ionization-change scenario, and constrain a nominal outflow location of < 0.5 pc and a transverse size of ~0.01 pc using the transverse-motion scenario. Our findings are consistent with previous work, and provide evidence in support of BALs tracing compact flow geometries with small filling factors.

Very few low-ionization broad absorption line (LoBAL) QSOs have been found at high redshifts to date. One high-redshift LoBAL QSO, J0122+1216, was recently discovered at the Lijiang 2.4-m Telescope with an initial redshift determination of 4.76. Aiming to investigate its physical properties, we carried out follow-up observations in the optical and near-IR spectroscopy. Near-IR spectra from UKIRT and P200 confirms that it is a LoBAL, with a new redshift determination of $4.82\pm0.01$ based on the \mgii~ emission-line. The new \mgii~ redshift determination reveals strong blueshifts and asymmetry of the high-ionization emission lines. We estimated a black hole mass of $\sim 2.3\times 10^9 M_\odot$ and Eddington ratio of $\sim 1.0$ according to the empirical \mgii-based single-epoch relation and bolometric correction factor. It is possible that strong outflows are the result of an extreme quasar environment driven by the high Eddington ratio. A lower limit on the outflowing kinetic power ($>0.9\% L_{Edd}$) was derived from both emission and absorption lines, indicating these outflows play a significant role in the feedback process to regulate the growth of its black hole as well as host galaxy evolution.

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 study the time lags between the continuum emission of quasars at different wavelengths, based on more than four years of multi-band ($g$, $r$, $i$, $z$) light-curves in the Pan-STARRS Medium Deep Fields. As photons from different bands emerge from different radial ranges in the accretion disk, the lags constrain the sizes of the accretion disks. We select 240 quasars with redshifts $z \approx 1$ or $z \approx 0.3$ that are relatively emission line free. The light curves are sampled from day to month timescales, which makes it possible to detect lags on the scale of the light crossing time of the accretion disks. With the code JAVELIN, we detect typical lags of several days in the rest frame between the $g$ band and the $riz$ bands. The detected lags are $\sim 2-3$ times larger than the light crossing time estimated from the standard thin disk model, consistent with the recently measured lag in NGC5548 and micro-lensing measurements of quasars. The lags in our sample are found to increase with increasing luminosity. Furthermore, the increase in lags going from $g-r$ to $g-i$ and then to $g-z$ is slower than predicted in the thin disk model, particularly for high luminosity quasars. The radial temperature profile in the disk must be different from what is assumed. We also find evidence that the lags decrease with increasing line ratios between ultraviolet FeII lines and MgII, which may point to changes in the accretion disk structure at higher metallicity.

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

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 present results from the largest systematic investigation of broad absorption line (BAL) acceleration to date. We use spectra of 140 quasars from three Sloan Digital Sky Survey programs to search for global velocity offsets in BALs over timescales of ~2.5-5.5 years in the quasar rest frame. We carefully select acceleration candidates by requiring monolithic velocity shifts over the entire BAL trough, avoiding BALs with velocity shifts that might be caused by profile variability. The CIV BALs of two quasars show velocity shifts consistent with the expected signatures of BAL acceleration, and the BAL of one quasar shows a velocity-shift signature of deceleration. In our two acceleration candidates, we see evidence that the magnitude of the acceleration is not constant over time; the magnitudes of the change in acceleration for both acceleration candidates are difficult to produce with a standard disk-wind model or via geometric projection effects. We measure upper limits to acceleration and deceleration for 76 additional BAL troughs and find that the majority of BALs are stable to within about 3% of their mean velocities. The lack of widespread acceleration/deceleration could indicate that the gas producing most BALs is located at large radii from the central black hole and/or is not currently strongly interacting with ambient material within the host galaxy along our line of sight.

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.

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 use 317,000 emission-line galaxies from the Sloan Digital Sky Survey to investigate line-ratio selection of active galactic nuclei (AGNs). In particular, we demonstrate that "star formation dilution" by HII regions causes a significant bias against AGN selection in low-mass, blue, star-forming, disk-dominated galaxies. This bias is responsible for the observed preference of AGNs among high-mass, green, moderately star-forming, bulge-dominated hosts. We account for the bias and simulate the intrinsic population of emission-line AGNs using a physically-motivated Eddington ratio distribution, intrinsic AGN narrow line region line ratios, a luminosity-dependent Lbol/L[OIII] bolometric correction, and the observed Mbh-sigma relation. These simulations indicate that, in massive (log(M*/Msun) > 10) galaxies, AGN accretion is correlated with specific star formation rate but is otherwise uniform with stellar mass. There is some hint of lower black hole occupation in low-mass (log(M*/Msun) < 10) hosts, although our modeling is limited by uncertainties in measuring and interpreting the velocity dispersions of low-mass galaxies. The presence of star formation dilution means that AGNs contribute little to the observed strong optical emission lines (e.g., [OIII] and Ha) in low-mass and star-forming hosts. However the AGN population recovered by our modeling indicates that feedback by typical (low- to moderate-accretion) low-redshift AGNs has nearly uniform efficiency at all stellar masses, star formation rates, and morphologies. Taken together, our characterization of the observational bias and resultant AGN occupation function suggest that AGNs are unlikely to be the dominant source of star formation quenching in galaxies, but instead are fueled by the same gas which drives star formation activity.

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 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 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 new stellar velocity dispersion measurements for four luminous quasars with the NIFS instrument and the ALTAIR laser guide star adaptive optics system on the Gemini North 8-m telescope. Stellar velocity dispersion measurements and measurements of the supermassive black hole masses in luminous quasars are necessary to investigate the coevolution of black holes and galaxies, trace the details of accretion, and probe the nature of feedback. We find that higher-luminosity quasars with higher-mass black holes are not offset with respect to the MBH-sigma relation exhibited by lower-luminosity AGNs with lower-mass black holes, nor do we see correlations with galaxy morphology. As part of this analysis, we have recalculated the virial products for the entire sample of reverberation-mapped AGNs and used these data to redetermine the mean virial factor hfi that places the reverberation data on the quiescent M_BH-sigma relation. With our updated measurements and new additions to the AGN sample, we obtain <f> = 4.31 +/- 1.05, which is slightly lower than, but consistent with, most previous determinations.

We have used Chandra observations of eight bulge-less disk galaxies to identify new ultraluminous X-ray source (ULX) candidates, study their high mass X-ray binary (HMXB) population, and search for low-luminosity active galactic nuclei (AGN). We report the discovery of 16 new ULX candidates in our sample of galaxies. Eight of these are found in the star forming galaxy NGC 337, none of which are expected to be background contaminants. The HMXB luminosity function of NGC 337 implies a star formation rate (SFR) of 6.8$^{+4.4}_{-3.5}$ \msun\ yr$^{-1}$, consistent at 1.5$\sigma$ with a recent state of the art SFR determination. We also report the discovery of a bright ULX candidate (X-1) in ESO 501-23. X-1's spectrum is well fit by an absorbed power law with $\Gamma = 1.18^{+0.19}_{-0.11}$ and N$\rm{_H}$ = 1.13$^{+7.07}_{-1.13} \times 10^{20}$ cm$^{-2}$, implying a 0.3-8 keV flux of $1.08^{+0.05}_{-0.07} \times 10^{-12}$ \esc. Its X-ray luminosity (L$_X$) is poorly constrained due to uncertainties in the host galaxy's distance, but we argue that its spectrum implies L$_X > 10^{40}$ \es. An optical counterpart to this object may be present in HST imaging. We also identify ULX candidates in IC 1291, PGC 3853, NGC 5964 and NGC 2805. We find no evidence of nuclear activity in the galaxies in our sample, placing a flux upper limit of $4 \times 10^{-15}$ \esc\ on putative AGN. Additionally, the type II-P supernova SN 2011DQ in NGC 337, which exploded 2 months before our X-ray observation, is undetected.

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 have undertaken a new ground-based monitoring campaign on the BLRG 3C390.3 to improve the measurement of the size of the BLR and to estimate the black hole mass. Optical spectra and g-band images were observed in 2005 using the 2.4m telescope at MDM Observatory. Integrated emission-line flux variations were measured for Ha, Hb, Hg, and for HeII4686, as well as g-band fluxes and the optical AGN continuum at 5100A. The g-band fluxes and the AGN continuum vary simultaneously within the uncertainties, tau=(0.2+-1.1)days. We find that the emission-line variations are delayed with respect to the variable g-band continuum by tau(Ha)=56.3(+2.4-6.6)days, tau(Hb)=44.3(+3.0_-3.3)days, tau(Hg)=58.1(+4.3-6.1)days, and tau(HeII4686)=22.3(+6.5-3.8)days. The blue and red peak in the double peaked line profiles, as well as the blue and red outer profile wings, vary simultaneously within +-3 days. This provides strong support for gravitationally bound orbital motion of the dominant part of the line emitting gas. Combining the time delay of Ha and Hb and the separation of the blue and red peak in the broad double-peaked profiles in their rms spectra, we determine Mbh_vir=1.77(+0.29-0.31)x10^8Msol and using sigma_line of the rms spectra Mbh_vir=2.60(+0.23-0.31)x10^8Msol for the central black hole of 3C390.3, respectively. Using the inclination angle of the line emitting region the mass of the black hole amounts to Mbh=0.86(+0.19-0.18)x10^9 Msol (peak-separation) and Mbh=1.26(+0.21-0.16)x10^9 Msol (sigma_line), respectively. This result is consistent with the black hole masses indicated by simple accretion disk models to describe the observed double-peaked profiles, derived from the stellar dynamics of 3C390.3, and with the AGN radius-luminosity relation. Thus, 3C390.3 as a radio-loud AGN with a low Eddington ratio, Ledd/Lbol=0.02, follows the same AGN radius-luminosity relation as radio-quiet AGN.

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.

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.

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 results of a search for nuclear X-ray activity in nearby galaxies using Chandra archival data in a sample of 62 galaxies from the Spitzer Infrared Nearby Galaxy Survey (SINGS). We detect 37 nuclear X-ray sources; seven of these are new detections. Most of the nuclear X-ray sources are likely to be AGNs. The fraction of galaxies hosting AGNs is thus about 60%, much higher than that found with optical searches, and demonstrates the efficacy of X-ray observations to find hidden AGNs in optically normal galaxies. We find that the nuclear X-ray sources are preferentially present in earlier type galaxies. Unlike what is observed at high redshift, we do not find a strong correlation between the AGN luminosity and the 24 micron luminosity of the host galaxy; we find a strong correlation with the 3.6 micron luminosity instead. This suggests that at the present epoch the accretion rate depends on the total mass of the galaxy, as perhaps does the black hole mass.

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.

Repeat scans by the Sloan Digital Sky Survey (SDSS) of a 278 square degree stripe along the Celestial equator have yielded an average of over 10 observations each for nearly 8,000 spectroscopically confirmed quasars. Over 2500 of these quasars are in the redshift range such that the CIV emission line is visible in the SDSS spectrum. Utilising the width of these CIV lines and the luminosity of the nearby continuum, we estimate black hole masses for these objects. In an effort to isolate the effects of black hole mass and luminosity on the photometric variability of our dataset, we create several subsamples by binning in these two physical parameters. By comparing the ensemble structure functions of the quasars in these bins, we are able to reproduce the well-known anticorrelation between luminosity and variability, now showing that this anticorrelation is independent of the black hole mass. In addition, we find a correlation between variability and the mass of the central black hole. By combining these two relations, we identify the Eddington ratio as a possible driver of quasar variability, most likely due to differences in accretion efficiency.