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.

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 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 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 isolate a set of quasars that exhibit emergent C iv broad absorption lines (BALs) in their spectra by comparing spectra in the SDSS Data Release 7 and the SDSS/BOSS Data Releases 9 and 10. After visually defining a set of emergent BALs, follow-up observations were obtained with the Gemini Observatory for 105 quasars. We find an emergence rate consistent with the previously reported disappearance rate of BAL quasars given the relative numbers of non-BAL and BAL quasars in the SDSS. We find candidate newly emerged BALs are preferentially drawn from among BALs with smaller balnicity indices, shallower depths, larger velocities, and smaller widths. Within two rest-frame years (average) after a BAL has emerged, we find it equally likely to continue increasing in equivalent width in an observation six months later (average) as it is to start decreasing. From the time separations between our observations, we conclude the coherence time-scale of BALs is less than 100 rest-frame days. We observe coordinated variability among pairs of troughs in the same quasar, likely due to clouds at different velocities responding to the same changes in ionizing flux; and the coordination is stronger if the velocity separation between the two troughs is smaller. We speculate the latter effect may be due to clouds having on average lower densities at higher velocities due to mass conservation in an accelerating flow, causing the absorbing gas in those clouds to respond on different timescales to the same ionizing flux variations.

We study the dramatic decrease in iron absorption strength in the iron low-ionization broad absorption line quasar SDSS J084133.15+200525.8. We report on the continued weakening of absorption in the prototype of this class of variable broad absorption line quasar, FBQS J140806.2+305448. We also report a third example of this class, SDSS J123103.70+392903.6; unlike the other two examples, it has undergone an increase in observed continuum brightness (at 3000~\AA\ rest-frame) as well as a decrease in iron absorption strength. These changes could be caused by absorber transverse motion or by ionization variability. We note that the \mgii\ and UV \feii\ lines in several FeLoBAL quasars are blueshifted by thousands of \kms\ relative to the \Hb\ emission line peak. We suggest that such emission arises in the outflowing winds normally seen only in absorption.

We present the discovery of the highest velocity CIV broad absorption line to date in the z=2.47 quasar SDSS J023011.28+005913.6, hereafter J0230. In comparing the public DR7 and DR9 spectra of J0230, we discovered an emerging broad absorption trough outflowing at ~60,000 km/s, which we refer to as trough A. In pursuing follow up observations of trough A, we discovered a second emergent CIV broad absorption trough outflowing at ~40,000 km/s, namely trough B. In total, we collected seven spectral epochs of J0230 that demonstrate emergent and rapidly (~10 days in the rest-frame) varying broad absorption. We investigate two possible scenarios that could cause these rapid changes: bulk motion and ionization variability. Given our multi-epoch data, we were able to rule out some simple models of bulk motion, but have proposed two more realistic models to explain the variability of both troughs. Trough A is likely an augmented `crossing disk' scenario with the absorber moving at 10,000 < v (km/s) < 18,000. Trough B can be explained by a flow-tube feature travelling across the emitting region at 8,000 < v (km/s) < 56,000. If ionization variability is the cause for the changes observed, trough A's absorber has n_e >= 724 cm^-3 and is at r_{equal} >=2.00 kpc, or is at r < 2.00 kpc with no constraint on the density; trough B's absorber either has n_e >= 1540 cm^-3 and is at r_{equal} >= 1.37 kpc, or is at r < 1.37 kpc with no constraint on the density.

Broad Absorption Line (BAL) trough variability is predominantly due to cloud motion transverse to our line of sight. The rate at which the variability occurs indicates the velocity of the cloud, which can provide constraints on the cloud's distance from the central source. This requires detailed spectroscopy during a variability event. Such spectra have proven elusive, suggesting either the timescale of variability is too short to be caught, or too long to notice until a sufficient amount of time has passed. Photometric monitoring of BAL quasar colours may potentially be used as an early warning system to trigger time resolved spectroscopic monitoring of BAL variability. Towards this end, we are analyzing both BAL and non-BAL colour variability using time series photometry from Stripe 82 in the Sloan Digital Sky Survey.

We test whether the Tinker & Chen model of MgII absorption due to the gaseous halo around a galaxy can reproduce absorption in quasar pairs (both lensed and physical) and lensed triples and quads from the literature. These quasars exhibit absorption from a total of 38 MgII systems spanning z=0.043 - 2.066 with mean redshift <z>=1.099 and weighted mean rest-frame equivalent width of 0.87 Ang. Using the Tinker & Chen model to generate simulated sight-lines, we marginalize the unknown parameters of the absorbing galaxies: dark matter halo mass, impact parameter, and azimuthal angle on the sky. We determine the ability of the model to statistically reproduce the observed variation in MgII absorption strength between paired sight-lines for different values of the gas covering fraction f_c and the characteristic length scale ell_c, within which the variation in absorption equivalent widths between sight-lines exponentially decreases. We find a best-fit f_c=0.60 \pm 0.15 and ell_c<8/h_70 kpc (1\sigma confidence limits), with smaller f_c allowed at larger ell_c. At 99.7% confidence, we are able to rule out f_c>0.87 for all values of ell_c and the region where ell_c<1.0/h_70 kpc and f_c<0.3.

We report sensitive Chandra X-ray non-detections of two unusual, luminous Iron Low-Ionization Broad Absorption Line Quasars (FeLoBALs). The observations do detect a non-BAL, wide-binary companion quasar to one of the FeLoBAL quasars. We combine X-ray-derived column density lower limits (assuming solar metallicity) with column densities measured from ultraviolet spectra and CLOUDY photoionization simulations to explore whether constant density slabs at broad line region densities can match the physical parameters of these two BAL outflows, and find that they cannot. In the "overlapping-trough" object SDSS J0300+0048, we measure the column density of the X-ray absorbing gas to be N_H >= 1.8 x 1024 cm-2. From the presence of Fe II UV78 absorption but lack of Fe II UV195/UV196 absorption, we infer the density in that part of the absorbing region to be n_e ~ 106 cm-3. We do find that a slab of gas at that density might be able to explain this object's absorption. In the Fe III-dominant object SDSS J2215-0045, the X-ray absorbing column density of N_H >= 3.4 x 1024 cm-2 is consistent with the Fe III-derived N_H >= 2 x 1022 cm-2 provided the ionization parameter is log U > 1.0 for both the n_e = 1011 cm-3 and n_e = 1012 cm-3 scenarios considered (such densities are required to produce Fe III absorption without Fe II absorption). However, the velocity width of the absorption rules out its being concentrated in a single slab at these densities. Instead, this object's spectrum can be explained by a low density, high ionization and high temperature disk wind that encounters and ablates higher density, lower ionization Fe III-emitting clumps.