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.

Understanding the links between the activity of supermassive black holes (SMBH) at the centres of galaxies and their host dark matter haloes is a key question in modern astrophysics. The final data release of the SDSS-IV eBOSS provides the largest contemporary spectroscopic sample of galaxies and QSOs. Using this sample and covering the redshift interval $z=0.7-1.1$, we have measured the clustering properties of the eBOSS QSOs, Emission Line Galaxies (ELGs) and Luminous Red Galaxies (LRGs). We have also measured the fraction of QSOs as a function of the overdensity defined by the galaxy population. Using these measurements, we investigate how QSOs populate and sample the galaxy population, and how the host dark-matter haloes of QSOs sample the underlying halo distribution. We find that the probability of a galaxy hosting a QSO is independent of the host dark matter halo mass of the galaxy. We also find that about 60\% of eBOSS QSOs are hosted by LRGs and about 20-40\% of QSOs are hosted by satellite galaxies. We find a slight preference for QSOs to populate satellite galaxies over central galaxies. This is connected to the host halo mass distribution of different types of galaxies. Based on our analysis, QSOs should be hosted by a very broad distribution of haloes, and their occurrence should be modulated only by the efficiency of galaxy formation processes.

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 combine the most precise small scale ($< 100\, \rm h^{-1}kpc$) quasar clustering constraintsto date with recent measurements at large scales ($> 1\, \rm h^{-1}Mpc$) from the extended Baryon Oscillation Spectroscopic Survey (eBOSS) to better constrain the satellite fraction of quasars at $z\sim 1.5$ in the halo occupation formalism. We build our Halo Occupation Distribution (HOD) framework based on commonly used analytic forms for the one and two-halo terms with two free parameters: the minimum halo mass that hosts a central quasar and the fraction of satellite quasars that are within one halo. Inspired by recent studies that propose a steeper density profile for the dark matter haloes that host quasars, we explore HOD models at kiloparsec scales and best-fit parameters for models with $10\times$ higher concentration parameter. We find that an HOD model with a satellite fraction of $f_{\rm sat} = 0.071_{-0.004}^{+0.009}$ and minimum mass of $\rm M_{m} = 2.31_{-0.38}^{+0.41} \times 10^{12}\, \, \rm h^{-1} M_{\odot}$ for the host dark matter haloes best describes quasar clustering (on all scales) at $z \sim 1.5$. Our results are marginally inconsistent with earlier work that studied brighter quasars, hinting at a luminosity-dependence to the one-halo term.

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 examine recent high-precision measurements of small-scale quasar clustering (at $z\sim0.5-2$ on scales of $\sim25~\mathrm{kpc/h}$) from the SDSS in the context of the MassiveBlackII (MBII) cosmological hydrodynamic simulation and conditional luminosity function (CLF) modeling. At these high luminosities ($g < 20.85$ quasars), the MBII simulation volume ($100~\mathrm{cMpc}/h$ comoving boxsize) has only 3 quasar pairs at distances of $1-4$ Mpc. The black-hole masses for the pairs range between $M_{bh}\sim1-3\times 10^{9}~M_{\odot}/h$ and the quasar hosts are haloes of $M_h\sim1-3\times10^{14}~M_{\odot}/h$. Such pairs show signs of recent major mergers in the MBII simulation. By modeling the central and satellite AGN CLFs as log-normal and Schechter distributions respectively (as seen in MBII AGNs), we arrive at CLF models which fit the simulation predictions and observed luminosity function and the small-scale clustering measured for the SDSS sample. The small-scale clustering of our mock quasars is well-explained by central-satellite quasar pairs that reside in $M_h>10^{14}~M_{\odot}/h$ dark matter haloes. For these pairs, satellite quasar luminosity is similar to that of central quasars. Our CLF models imply a relatively steep increase in the maximum satellite luminosity, $L^*_{\mathrm{sat}}$, in haloes of $M_h>10^{14}~M_{\odot}/h$ with associated larger values of $L^*_{\mathrm{sat}}$ at higher redshift. This leads to increase in the satellite fraction that manifests itself in an enhanced clustering signal at $\lesssim$ 1 Mpc/h. For the ongoing eBOSS-CORE sample, we predict $\sim 200-500$ quasar pairs at $z\sim1.5$ (with $M_h \gtrsim10^{13}~M_{\odot}/h$ and $M_{bh} \gtrsim10^{8}~M_{\odot}/h$) at $\sim25~\mathrm{kpc}$ scales. Such a sample would be $\gtrsim10$ times larger than current pair samples.

The DESI Legacy Imaging Surveys are a combination of three public projects (the Dark Energy Camera Legacy Survey, the Beijing-Arizona Sky Survey, and the Mayall z-band Legacy Survey) that will jointly image approximately 14,000 deg^2 of the extragalactic sky visible from the northern hemisphere in three optical bands (g, r, and z) using telescopes at the Kitt Peak National Observatory and the Cerro Tololo Inter-American Observatory. The combined survey footprint is split into two contiguous areas by the Galactic plane. The optical imaging is conducted using a unique strategy of dynamically adjusting the exposure times and pointing selection during observing that results in a survey of nearly uniform depth. In addition to calibrated images, the project is delivering a catalog, constructed by using a probabilistic inference-based approach to estimate source shapes and brightnesses. The catalog includes photometry from the grz optical bands and from four mid-infrared bands (at 3.4, 4.6, 12 and 22 micorons) observed by the Wide-field Infrared Survey Explorer (WISE) satellite during its full operational lifetime. The project plans two public data releases each year. All the software used to generate the catalogs is also released with the data. This paper provides an overview of the Legacy Surveys project.

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

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

We analyse the clustering of the Sloan Digital Sky Survey IV extended Baryon Oscillation Spectroscopic Survey Data Release 14 quasar sample (DR14Q). We measure the redshift space distortions using the power spectrum monopole, quadrupole and hexadecapole inferred from 148,659 quasars between redshifts 0.8 and 2.2 covering a total sky footprint of 2112.9 deg$^2$. We constrain the logarithmic growth of structure times the amplitude of dark matter density fluctuations, $f\sigma_8$, and the Alcock-Paczynski dilation scales which allow constraints to be placed on the angular diameter distance $D_A(z)$ and the Hubble $H(z)$ parameter. At the effective redshift of $z_{\rm eff}=1.52$, $f\sigma_8(z_{\rm eff})=0.420\pm0.076$, $H(z_{\rm eff})=[162\pm 12]\, (r_s^{\rm fid}/r_s)\,{\rm km\, s}^{-1}{\rm Mpc}^{-1}$, and $D_A(z_{\rm eff})=[1.85\pm 0.11]\times10^3\,(r_s/r_s^{\rm fid})\,{\rm Mpc}$, where $r_s$ is the comoving sound horizon at the baryon drag epoch and the superscript `fid' stands for its fiducial value. The errors take into account the full error budget, including systematics and statistical contributions. These results are in full agreement with the current $\Lambda$-Cold Dark Matter ($\Lambda$CDM) cosmological model inferred from Planck measurements. Finally, we compare our measurements with other eBOSS companion papers and find excellent agreement, demonstrating the consistency and complementarity of the different methods used for analysing the data.

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.

Recent work has found evidence for a difference in the bias and dark matter halo masses of WISE-selected obscured and unobscured quasars, implying a distinction between these populations beyond random line-of-sight effects. However, the significance of this difference in the most up-to-date measurements is relatively weak, at $\sim$2$\sigma$ for individual measurements but bolstered by agreement from different techniques, including angular clustering and cross-correlations with cosmic microwave background (CMB) lensing maps. Here, we expand the footprint of previous work, aiming to improve the precision of both methods. In this larger area we correct for position dependent selection effects, in particular fluctuations of the WISE-selected quasar density as a function of Galactic latitude. We also measure the cross-correlation of the obscured and unobscured samples and confirm that they are well-matched in redshift, both centred at $z=1$. Combined with very similar detection fractions and magnitude distributions in the long-wavelength WISE bands, this redshift match strongly supports the fact that IR selection identifies obscured and unobscured quasars of similar bolometric luminosity. Finally, we perform cross-correlations with confirmed spectroscopic quasars, again confirming the results from other methods --- obscured quasars reside in haloes a factor of 3 times more massive than unobscured quasars. This difference is significant at the $\sim$5$\sigma$ level when the measurements are combined, strong support for the idea that obscuration in at least some quasars is tied to the larger environment, and may have an evolutionary component.

We study the first year of the eBOSS quasar sample in the redshift range $0.9<z<2.2$ which includes 68,772 homogeneously selected quasars. We show that the main source of systematics in the evaluation of the correlation function arises from inhomogeneities in the quasar target selection, particularly related to the extinction and depth of the imaging data used for targeting. We propose a weighting scheme that mitigates these systematics. We measure the quasar correlation function and provide the most accurate measurement to date of the quasar bias in this redshift range, $b_Q = 2.45 \pm 0.05$ at $\bar z=1.55$, together with its evolution with redshift. We use this information to determine the minimum mass of the halo hosting the quasars and the characteristic halo mass, which we find to be both independent of redshift within statistical error. Using a recently-measured quasar-luminosity-function we also determine the quasar duty cycle. The size of this first year sample is insufficient to detect any luminosity dependence to quasar clustering and this issue should be further studied with the final $\sim$500,000 eBOSS quasar sample.

We present the most precise estimate to date of the clustering of quasars on very small scales, based on a sample of 47 binary quasars with magnitudes of $g<20.85$ and proper transverse separations of $\sim 25\,h^{-1}$\,kpc. Our sample of binary quasars, which is about 6 times larger than any previous spectroscopically confirmed sample on these scales, is targeted using a Kernel Density Estimation technique (KDE) applied to Sloan Digital Sky Survey (SDSS) imaging over most of the SDSS area. Our sample is "complete" in that all of the KDE target pairs with $17.0 \lesssim R \lesssim 36.2\,h^{-1}$\,kpc in our area of interest have been spectroscopically confirmed from a combination of previous surveys and our own long-slit observational campaign. We catalogue 230 candidate quasar pairs with angular separations of $<8\arcsec$, from which our binary quasars were identified. We determine the projected correlation function of quasars ($\bar W_{\rm p}$) in four bins of proper transverse scale over the range $17.0 \lesssim R \lesssim 36.2\,h^{-1}$\,kpc. The implied small-scale quasar clustering amplitude from the projected correlation function, integrated across our entire redshift range, is $A=24.1\pm3.6$ at $\sim 26.6 ~h^{-1}$\,kpc. Our sample is the first spectroscopically confirmed sample of quasar pairs that is sufficiently large to study how quasar clustering evolves with redshift at $\sim 25 ~h^{-1}$ kpc. We find that empirical descriptions of how quasar clustering evolves with redshift at $\sim 25 ~h^{-1}$ Mpc also adequately describe the evolution of quasar clustering at $\sim 25 ~h^{-1}$ kpc.

DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$\alpha$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.

The BOSS quasar sample is used to study cosmic homogeneity with a 3D survey in the redshift range $2.2<z<2.8$. We measure the count-in-sphere, $N(<\! r)$, i.e. the average number of objects around a given object, and its logarithmic derivative, the fractal correlation dimension, $D_2(r)$. For a homogeneous distribution $N(<\! r) \propto r^3$ and $D_2(r)=3$. Due to the uncertainty on tracer density evolution, 3D surveys can only probe homogeneity up to a redshift dependence, i.e. they probe so-called "spatial isotropy". Our data demonstrate spatial isotropy of the quasar distribution in the redshift range $2.2<z<2.8$ in a model-independent way, independent of any FLRW fiducial cosmology, resulting in $3-\langle D_2 \rangle < 1.7 \times 10^{-3}$ (2 $\sigma$) over the range $250<r<1200 \, h^{-1}$Mpc for the quasar distribution. If we assume that quasars do not have a bias much less than unity, this implies spatial isotropy of the matter distribution on large scales. Then, combining with the Copernican principle, we finally get homogeneity of the matter distribution on large scales. Alternatively, using a flat $\Lambda$CDM fiducial cosmology with CMB-derived parameters, and measuring the quasar bias relative to this $\Lambda$CDM model, our data provide a consistency check of the model, in terms of how homogeneous the Universe is on different scales. $D_2(r)$ is found to be compatible with our $\Lambda$CDM model on the whole $10<r<1200 \, h^{-1}$Mpc range. For the matter distribution we obtain $3-\langle D_2 \rangle < 5 \times 10^{-5}$ (2 $\sigma$) over the range $250<r<1200 \, h^{-1}$Mpc, consistent with homogeneity on large scales.

We report the discovery of a quasar pair at z=5 separated by 21 arcsec. Both objects were identified as quasar candidates using simple color selection techniques applied to photometric catalogs from the CFHT Legacy Survey (CFHTLS). Spectra obtained with the MMT present no discernible offset in redshift between the two objects; on the other hand, there are clear differences in the emission line profiles and in the multiwavelength spectral energy distributions that strongly disfavor the hypothesis that they are gravitationally lensed images of a single quasar. Both quasars are surprisingly bright given their proximity (a projected separation of ~135 kpc), with i=19.4 and i=21.4. Previous measurements of the luminosity function demonstrate that luminous quasars are extremely rare at z=5; the existence of this pair suggests that quasars have strong small-scale clustering at high redshift. Assuming a real-space correlation function of the form $\xi(r) \propto (r/r_0)^{-2}$, this discovery implies a correlation length $r_0 > 20 h^{-1}$ Mpc, consistent with a rapid strengthening of quasar clustering at high redshift as seen in previous observations and predicted by theoretical models where feedback effects are inefficient at shutting down black hole growth at high redshift.

We measure the two-point clustering of spectroscopically confirmed quasars from the final sample of the Baryon Oscillation Spectroscopic Survey (BOSS) on comoving scales of 4 < s < 22 Mpc/h. The sample covers 6950 deg^2 (~ 19 (Gpc/h)^3) and, over the redshift range 2.2 < z < 2.8, contains 55,826 homogeneously selected quasars, which is twice as many as in any similar work. We deduce b_Q = 3.54 +/- 0.10 ; the most precise measurement of quasar bias to date at these redshifts. This corresponds to a host halo mass of ~ 2 x 10^12 ~ M_sun/h with an implied quasar duty cycle of ~1 percent. The real-space projected correlation function is well-fit by a power law of index -2 and correlation length r0 = (8.12 +/- 0.22), Mpc/h over scales of 4 < rp < 25 ~ Mpc/h. To better study the evolution of quasar clustering at moderate redshift, we extend the redshift range of our study to z ~ 3.4 and measure the bias and correlation length of three subsamples over 2.2 < z < 3.4. We find no significant evolution of r0 or bias over this range, implying that the host halo mass of quasars decreases somewhat with increasing redshift. We find quasar clustering remains similar over a decade in luminosity, contradicting a scenario in which quasar luminosity is monotonically related to halo mass at z ~ 2.5. Our results are broadly consistent with previous BOSS measurements, but they yield more precise constraints based upon a larger and more uniform data set.

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

We present a measurement of the quasar luminosity function in the range 0.68<z<4 down to extinction corrected magnitude g_dered=22.5, using a simple and well understood target selection technique based on the time-variability of quasars. The completeness of our sample was derived directly from a control sample of quasars, without requiring complex simulations of quasar light-curves or colors. A total of 1877 quasar spectra were obtained from dedicated programs on the Sloan telescope (as part of the SDSS-III/BOSS survey) and on the Multiple Mirror Telescope. They allowed us to derive the quasar luminosity function. It agrees well with previously published results in the redshift range 0.68<z<2.6. Our deeper data allow us to extend the measurement to z=4. We measured quasar densities to g_dered<22.5, obtaining 30 QSO per deg^2 at z<1, 99 QSO per deg^2 for 1<z<2.15, and 47 QSO per deg^2 at z>2.15. Using pure luminosity evolution models, we fitted our LF measurements, and predicted quasar number counts as a function of redshift and observed magnitude. These predictions are useful inputs for future cosmology surveys such as those relying on the observation of quasars to measure baryon acoustic oscillations.

We have used new deep observations of the Coma cluster from GALEX to identify 13 star-forming galaxies with asymmetric morphologies in the ultraviolet. Aided by optical broad-band and H-alpha imaging, we interpret the asymmetric features as being due to star formation within gas stripped from the galaxies by interaction with the cluster environment. The selected objects display a range of structures from broad fan-shaped systems of filaments and knots (`jellyfish') to narrower and smoother tails extending up to 100 kpc in length. Some of the features have been discussed previously in the literature, while others are newly identified here. As an ensemble, the candidate stripping events are located closer to the cluster centre than other star-forming galaxies; their radial distribution is similar to that of all cluster members, dominated by passive galaxies. The fraction of blue galaxies which are undergoing stripping falls from 40% in the central 500 kpc, to less than 5% beyond 1 Mpc. We find that tails pointing away from (i.e. galaxies moving towards) the cluster centre are strongly favoured (11/13 cases). From the small number of `outgoing' galaxies with stripping signatures we conclude that the stripping events occur primarily on first passage towards the cluster centre, and are short-lived compared to the cluster crossing time. Using infall trajectories from simulations, the observed fraction of blue galaxies undergoing stripping can be reproduced if the events are triggered at a threshold radius of ~1 Mpc and detectable for ~500 Myr. HST images are available for two galaxies from our sample and reveal compact blue knots coincident with UV and H-alpha emission, apparently forming stars within the stripped material. Our results confirm that stripping of gas from infalling galaxies, and associated star formation in the stripped material, is a widespread phenomenon in rich clusters.