We investigate the responsiveness of the 2798AA Mg II broad emission line in AGN on timescales of several years. Our study is based on a sample of extremely variable AGN as well as a broad population sample. The observed response of the line in previous studies has been mixed. By focussing on extreme variability ($|\Delta g|>$ 1) we find that Mg II clearly does respond to the continuum. However, the degree of responsiveness varies strikingly from one object to another; we see cases of Mg II changing by as much as the continuum, more than the continuum, or very little at all. In 74% of the highly variable sample the behaviour of Mg II corresponds with that of H$\beta$, with 30% of the objects showing large variations in both lines. We do not detect any change in the line width that would correspond to Broad Line Region `breathing', in accordance with results from literature. Some of the objects in our highly variable sample show a clear asymmetry in the Mg II profile. This skewness can be both to the blue and the red of the line centre. Results from our broad population sample show that highly variable quasars have lower Eddington ratios. This result holds for the variability of the continuum, but the correlation is significantly reduced for the variability of the Mg II line. For the first time, we present an overview of the value of the intrinsic Baldwin Effect for Mg II in a large sample.

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

Despite their factor of ~10^8 difference in black hole mass, several lines of evidence suggest possible similarities between black hole accretion flows in active galactic nuclei (AGN) and Galactic X-ray binaries. However, it is still unclear whether the geometry of the disk-corona system in X-ray binaries directly scale up to AGN, and whether this analogy still holds in different accretion states. We test this AGN/X-ray binary analogy, by comparing the observed correlations between the UV-to-X-ray spectral index (alpha_OX) and Eddington ratio in AGN to those predicted from observations of X-ray binary outbursts. This approach probes the geometry of their disk-corona systems as they transition between different accretion states. We use new Chandra X-ray and ground-based rest-UV observations of faded 'changing-look' quasars to extend this comparison to lower Eddington ratios of <10^-2, where observations of X-ray binaries predict a softening of alpha_OX in AGN. We find that the observed correlations between alpha_OX and Eddington ratio of AGN displays a remarkable similarity to accretion state transitions in prototypical X-ray binary outbursts, including an inversion of this correlation at a critical Eddington ratio of ~10^-2. Our results suggest that the structures of black hole accretion flows directly scale across a factor of ~10^8 in black hole mass and across different accretion states, enabling us to apply theoretical models of X-ray binaries to explain AGN phenomenology.

Physical models of X-ray binary outbursts can aid in understanding the origin of 'changing-look' active galactic nuclei (AGN), if we can establish that these two black hole accretion phenomena are analogous. Previously, studies of the correlation between the UV-to-X-ray spectral index alpha_OX and Eddington ratio using single-epoch observations of changing-look AGN samples have revealed possible similarities to the spectral evolution of outbursting X-ray binaries. However, direct comparisons using multi-epoch UV/X-ray light curves of individual changing-look AGN undergoing dramatic changes in Eddington ratio have been scarce. Here, we use published Swift UV/X-ray light curves of two changing-look AGN (NGC 2617 and ZTF18aajupnt) to examine the evolution of their alpha_OX values during outburst. We show that the combination of these two changing-look AGN can trace out the predicted spectral evolution from X-ray binary outbursts, including the inversion in the evolution of alpha_OX as a function of Eddington ratio. We suggest that the spectral softening that is observed to occur below a critical Eddington ratio in both AGN and X-ray binaries is due to reprocessing of Comptonized X-ray emission by the accretion disk, based on the X-ray to UV reverberation lags previously observed in NGC 2617. Our results suggest that the physical processes causing the changing-look AGN phenomenon are similar to those in X-ray binary outbursts.

We study the sudden optical and ultraviolet (UV) brightening of 1ES 1927+654, which until now was known as a narrow-line active galactic nucleus (AGN). 1ES 1927+654 was part of the small and peculiar class of "true Type-2" AGN, which lack broad emission lines and line-of-sight obscuration. Our high-cadence spectroscopic monitoring captures the appearance of a blue, featureless continuum, followed several weeks later by the appearance of broad Balmer emission lines. This timescale is generally consistent with the expected light travel time between the central engine and the broad-line emission region in (persistent) broad-line AGN. Hubble Space Telescope spectroscopy reveals no evidence for broad UV emission lines (e.g., CIV1549, CIII]1909, MgII2798), probably owing to dust in the broad-line emission region. To the best of our knowledge, this is the first case where the lag between the change in continuum and in broad-line emission of a "changing-look" AGN has been temporally resolved. The nature and timescales of the photometric and spectral evolution disfavor both a change in line-of-sight obscuration and a change of the overall rate of gas inflow as driving the drastic spectral transformations seen in this AGN. Although the peak luminosity and timescales are consistent with those of tidal disruption events seen in inactive galaxies, the spectral properties are not. The X-ray emission displays a markedly different behavior, with frequent flares on timescales of hours to days, and will be presented in a companion publication.

Microlensing by stars within distant galaxies acting as strong gravitational lenses of multiply-imaged quasars, provides a unique and direct measurement of the internal structure of the lensed quasar on nano-arcsecond scales. The measurement relies on the temporal variation of high-magnification caustic crossings which vary on timescales of days to years. Multiwavelength observations provide information from distinct emission regions in the quasar. Through monitoring of these strong gravitational lenses, a full tomographic view can emerge with Astronomical-Unit scale resolution. Work to date has demonstrated the potential of this technique in about a dozen systems. In the 2020s there will be orders of magnitude more systems to work with. Monitoring of lens systems for caustic-crossing events to enable triggering of multi-platform, multi-wavelength observations in the 2020s will fulfill the potential of quasar microlensing as a unique and comprehensive probe of active black hole structure and dynamics.

Dwarf carbon (dC) stars, main sequence stars showing carbon molecular bands, were initially thought to be an oxymoron since only AGB stars dredge carbon into their atmospheres. Mass transfer from a former AGB companion that has since faded to a white dwarf seems the most likely explanation. Indeed, a few types of giants known to show anomalous abundances --- notably, the CH, Ba and CEMP-s stars --- are known to have a high binary frequency. The dC stars may be the enhanced-abundance progenitors of most, if not all, of these systems, but this requires demonstrating a high binary frequency for dCs. Here, for a sample of 240 dC stars targeted for repeat spectroscopy by the SDSS-IV's Time Domain Spectroscopic Survey, we analyze radial velocity variability to constrain the binary frequency and orbital properties. A handful of dC systems show large velocity variability ($>$100 km s$^{-1}$). We compare the dCs to a control sample with a similar distribution of magnitude, color, proper motion, and parallax. Using MCMC methods, we use the measured $\Delta$RV distribution to estimate the binary fraction and the separation distribution assuming both a unimodal and bimodal distribution. We find the dC stars have an enhanced binary fraction of 95\%, consistent with them being products of mass transfer. These models result in mean separations of less than 1 AU corresponding to periods on the order of 1 year. Our results support the conclusion that dC stars form from close binary systems via mass transfer.

Changing-look quasars are a newly-discovered class of luminous active galactic nuclei that undergo rapid ($\lesssim$10 year) transitions between Type 1 and Type 1.9/2, with an associated change in their continuum emission. We characterize the host galaxies of four faded changing-look quasars using broadband optical imaging. We use \textit{gri} images obtained with the Gemini Multi Object Spectrograph (GMOS) on Gemini North to characterize the surface brightness profiles of the quasar hosts and search for [O III] $\lambda4959,\lambda5007$ emission from spatially extended regions, or voorwerpjes, with the goal of using them to examine past luminosity history. Although we do not detect, voorwerpjes surrounding the four quasar host galaxies, we take advantage of the dim nuclear emission to characterize the colors and morphologies of the host galaxies. Three of the four galaxies show morphological evidence of merger activity or tidal features in their residuals. The three galaxies which are not highly distorted are fit with a single S\'ersic profile to characterize their overall surface brightness profiles. The single-S\'ersic fits give intermediate S\'ersic indices between the $n=1$ of disk galaxies and the $n=4$ of ellipticals. On a color-magnitude diagram, our changing-look quasar host galaxies reside in the blue cloud, with other AGN host galaxies and star-forming galaxies. On a color-S\'ersic index diagram the changing-look quasar hosts reside with other AGN hosts in the "green valley". Our analysis suggests that the hosts of changing-look quasars are predominantly disrupted or merging galaxies that resemble AGN hosts, rather than inactive galaxies.

Time-domain Astrophysics (TDA), a foundation of Astronomy, has become a major part of current and projected (2020s) astrophysics. While much has been derived from temporal measures of flux and color, the real physics comes from spectroscopy. With LSST coming on line in 2022, with TDA as one of its original drivers, the deluge of Transients and new types of variables will be truly astronomical. With multi-wavelength targeted EM surveys and multi-messenger (e.g. LIGO-international and advanced ICECUBE), and the possibility of full-sky and full-time X-ray imagers, the discovery of new Transients and Variables will flood telescopes on the ground and in space, and this just for multi-band imaging without spectroscopy. In this White Paper we briefly summarize several long-standing major science objectives that can be realized with TDA imaging and spectroscopy (near UV to mid IR) from space. We provide a brief description of how these can be achieved with the Time-domain Spectroscopic Observatory (TSO), a Probe-Class mission concept that ELTs on the ground and Flagship missions in space can not achieve on their own.

Long gamma-ray bursts (GRBs) are the most luminous known electromagnetic radiation sources in the Universe for the 3 to 300 sec of their prompt flashes (isotropic X/ gamma-ray luminosities up to 10^53 ergs/sec). Their afterglows have first day rest-frame UV/optical absolute magnitudes AB = -30 to -23. This luminous continuum nUV-nIR back-light provides the ultimate probe of the SFR(z) back to the first Pop III to II.5 stars, expected to be massive and GRB progenitors. GRB afterglow spectra in the first 1 to 3 hours will directly measure their host galaxy ionization fraction x_i vs. z in the Epoch of Reionization (EOR), tracing the growth of structure. Only 28% of Swift GRBs have measured redshifts due to limited followup at R, J >21. Some ~25% of GRBs are optically dark due to dust absorption in their host galaxies, but those with low NH in their X-ray spectra are likely at z >7. Current 8-10m telescopes and coming ELTs cannot pursue optically dark GRBs promptly, nor can JWST or WFIRST slew within 0.5 to 1 days of a GRB. The Time-domain Spectroscopic Observatory (TSO) is a proposed Probe-class 1.3m telescope at L2, with imaging and spectroscopy (R = 200, 1800) in 4 bands (0.3 to 5um) and rapid slew capability to 90% of sky. TSO would finally utilize z > 6 to 12 GRBs as the most direct probe of the SFR(z), EOR(z), and possibly the first direct detection of the core collapse of the very first (Pop III) stars.

With the epochal first detection of gravitational waves from a binary neutron star (NS) merger with the GW170817 event, and its direct confirmation that NS-NS mergers are significant sources of the of the r-process nucleosynthesis of heavy elements, an immense new arena for prompt EM (X-rays through IR and radio) studies of fundamental physics has been opened. Over the next decade, GW observatories will expand in scale and sensitivity so the need for facilities that can provide prompt, high sensitivity, broad-band EM followup becomes more urgent. NS-NS or NS-black hole (BH) mergers will be instantly recognized (and announced) by the LIGO-international collaboration. LSST will be a prime resource for rapid tiling of what will usually be large (~10-100 degree squared) error boxes. X-ray through IR Telescopes in space with (nearly) full-sky access that can rapidly image and tile are crucial for providing the earliest imaging and spectroscopic studies of the kilonova emission immediately following NS-NS mergers. The Time-domain Spectroscopic Observatory (TSO) is a proposed Probe-class 1.3 m 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. TSO nUV-mid-IR spectra will enable new constraints on NS structure and nucleosynthesis.

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.

We report the results of a Sloan Digital Sky Survey-IV eBOSS program to target X-ray sources and mid-infrared-selected WISE AGN candidates in a 36.8 deg$^2$ region of Stripe 82. About half this survey (15.6 deg$^2$) covers the largest contiguous portion of the Stripe 82 X-ray survey. This program represents the largest spectroscopic survey of AGN candidates selected solely by their WISE colors. We combine this sample with X-ray and WISE AGN in the field identified via other sources of spectroscopy, producing a catalog of 4847 sources that is 82% complete to $r\sim22$. Based on X-ray luminosities or WISE colors, 4730 of these sources are AGN, with a median sample redshift of $z\sim1$. About 30% of the AGN are optically obscured (i.e., lack broad lines in their optical spectra). BPT analysis, however, indicates that 50% of the WISE AGN at $z<0.5$ have emission line ratios consistent with star-forming galaxies, so whether they are buried AGN or star-forming galaxy contaminants is currently unclear. We find that 61% of X-ray AGN are not selected as MIR AGN, with 22% of X-ray AGN undetected by WISE. Most of these latter AGN have high X-ray luminosities ($L_{\rm x} > 10^{44}$ erg s$^{-1}$), indicating that MIR selection misses a sizable fraction of the highest luminosity AGN, as well as lower luminosity sources where AGN heated dust is not dominating the MIR emission. Conversely, $\sim$58% of WISE AGN are undetected by X-rays, though we do not find that they are preferentially redder than the X-ray detected WISE AGN.

SDSS J2232-0806 (the 'Big Dipper') has been identified as a 'slow-blue nuclear hypervariable': a galaxy with no previously known active nucleus, blue colours and large-amplitude brightness evolution occurring on a timescale of years. Subsequent observations have shown that this source does indeed contain an active galactic nucleus (AGN). Our optical photometric and spectroscopic monitoring campaign has recorded one major dimming event (and subsequent rise) over a period of around four years; there is also evidence of previous events consistent with this in archival data recorded over the last twenty years. Here we report an analysis of the eleven optical spectra obtained to date and we assemble a multiwavelength data set including infrared, ultraviolet and X-ray observations. We find that an intrinsic change in the luminosity is the most favoured explanation of the observations, based on a comparison of continuum and line variability and the apparent lagged response of the hot dust. This source, along with several other recently-discovered 'changing-look' objects, demonstrate that AGN can exhibit large-amplitude luminosity changes on timescales much shorter than those predicted by standard thin accretion disc models.

Active galactic nuclei (AGN) that show strong rest-frame optical/UV variability in their blue continuum and broad line emission are classified as "changing-look" AGN, or at higher luminosities changing look quasars (CLQs). These surprisingly large and sometimes rapid transitions challenge accepted models of quasar physics and duty cycles, offer several new avenues for study of quasar host galaxies, and open a wider interpretation of the cause of differences between broad and narrow line AGN. To better characterize extreme quasar variability, we present follow-up spectroscopy as part of a comprehensive search for CLQs across the full SDSS footprint using spectroscopically confirmed quasars from the SDSS DR7 catalog. Our primary selection requires large-amplitude (|\Delta g|>1 mag, |\Delta r|>0.5 mag) variability over any of the available time baselines probed by the SDSS and Pan-STARRS 1 surveys. We employ photometry from the Catalina Sky Survey to verify variability behavior in CLQ candidates where available, and confirm CLQs using optical spectroscopy from the William Herschel, MMT, Magellan, and Palomar telescopes. For our adopted S/N threshold on variability of broad H\beta emission, we find 17 new CLQs, yielding a confirmation rate of >~ 20%. These candidates are at lower Eddington ratio relative to the overall quasar population which supports a disk-wind model for the broad line region. Based on our sample, the CLQ fraction increases from 10% to roughly half as the continuum flux ratio between repeat spectra at 3420 Angstroms increases from 1.5 to 6. We release a catalog of over 200 highly variable candidates to facilitate future CLQ searches.

[Abridged] We present UV/optical observations of PS16aqv (SN 2016ard), a Type I superluminous supernova (SLSN-I) classified as part of our search for low-$z$ SLSNe. PS16aqv is a fast evolving SLSNe-I that reached a peak absolute magnitude of $M_{r} \approx -22.1$. The lightcurves exhibit a significant undulation at 30 rest-frame days after peak, with a behavior similar to undulations seen in the slowly fading SLSN-I SN 2015bn. This similarity strengthens the case that fast and slow SLSNe-I form a continuum with a common origin. At $\approx\!80$ days after peak, the lightcurves exhibit a transition to a slow decline, followed by significant subsequent steepening, indicative of a plateau phase or a second significant undulation. Deep limits at $\approx280$ days after peak imply a tight constraint on the nickel mass, $M_{\rm Ni} \lesssim 0.35$ M$_{\odot}$ (lower than for previous SLSNe-I), and indicate that some SLSNe-I do not produce significantly more nickel than normal Type Ic SNe. Using MOSFiT, we model the lightcurve with a magnetar central engine model and find $P_{\rm spin} \approx 0.9$ ms, $B \approx 1.5 \times 10^{14}$ G, and $M_{\rm ej} \approx 16$ M$_{\odot}$. The implied rapid spin-down time and large reservoir of available energy coupled with the high ejecta mass may account for the fast evolving lightcurve and slow spectroscopic evolution. We also study the location of PS16aqv in its host galaxy and find that it occurred at an offset of $2.46 \pm 0.21$ kpc from the central star-forming region. We find the host galaxy exhibits low metallicity and spatially varying extinction and star formation rate, with the explosion site exhibiting lower values than the central region. The complexity seen in the lightcurves of PS16aqv and other events highlights the importance of obtaining well-sampled lightcurves for exploring deviations from a uniform decline.

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.

We present the Data Release 14 Quasar catalog (DR14Q) from the extended Baryon Oscillation Spectroscopic Survey (eBOSS) of the Sloan Digital Sky Survey IV (SDSS-IV). This catalog includes all SDSS-IV/eBOSS objects that were spectroscopically targeted as quasar candidates and that are confirmed as quasars via a new automated procedure combined with a partial visual inspection of spectra, have luminosities $M_{\rm i} \left[ z=2 \right] < -20.5$ (in a $\Lambda$CDM cosmology with $H_0 = 70 \ {\rm km \ s^{-1} \ Mpc ^{-1}}$, $\Omega_{\rm M} = 0.3$, and $\Omega_{\rm \Lambda} = 0.7$), and either display at least one emission line with a full width at half maximum (FWHM) larger than $500 \ {\rm km \ s^{-1}}$ or, if not, have interesting/complex absorption features. The catalog also includes previously spectroscopically-confirmed quasars from SDSS-I, II and III. The catalog contains 526,356 quasars 144,046 are new discoveries since the beginning of SDSS-IV) detected over 9,376 deg$^2$ (2,044 deg$^2$ having new spectroscopic data available) with robust identification and redshift measured by a combination of principal component eigenspectra. The catalog is estimated to have about 0.5% contamination. The catalog identifies 21,877 broad absorption line quasars and lists their characteristics. For each object, the catalog presents SDSS five-band CCD-based photometry with typical accuracy of 0.03 mag. The catalog also contains X-ray, ultraviolet, near-infrared, and radio emission properties of the quasars, when available, from other large-area surveys.

We present an improved photometric error analysis for the 7,100 CRTS (Catalina Real-Time Transient Survey) optical light curves for quasars from the SDSS (Sloan Digital Sky Survey) Stripe 82 catalogue. The SDSS imaging survey has provided a time-resolved photometric data set which greatly improved our understanding of the quasar optical continuum variability: Data for monthly and longer time-scales are consistent with a damped random walk (DRW). Recently, newer data obtained by CRTS provided puzzling evidence for enhanced variability, compared to SDSS results, on monthly time-scales. Quantitatively, SDSS results predict about 0.06 mag root-mean-square (rms) variability for monthly time-scales, while CRTS data show about a factor of 2 larger rms, for spectroscopically confirmed SDSS quasars. Our analysis has successfully resolved this discrepancy as due to slightly underestimated photometric uncertainties from the CRTS image processing pipelines. As a result, the correction for observational noise is too small and the implied quasar variability is too large. The CRTS photometric error correction factors, derived from detailed analysis of non-variable SDSS standard stars that were re-observed by CRTS, are about 20-30%, and result in reconciling quasar variability behaviour implied by the CRTS data with earlier SDSS results. An additional analysis based on independent light curve data for the same objects obtained by the Palomar Transient Factory provides further support for this conclusion. In summary, the quasar variability constraints on weekly and monthly time-scales from SDSS, CRTS and PTF surveys are mutually compatible, as well as consistent with DRW model.

We present an optical to X-ray spectral analysis of four 'hypervariable' AGN (HVAs) discovered by comparing Pan-STARRS data to that from SDSS over a 10 year baseline (Lawrence et al 2016). There is some evidence that these objects are X-ray loud for their corresponding UV luminosities, but given that we measured them in a historic high state, it is not clear whether to take the high-state or low-state as typical of the properties of these HVAs. We estimate black hole masses based on MgII and H$\alpha$ emission line profiles, and either the high and low state luminosities, finding mass ranges $\log(M_{\rm BH}/M_{\odot}) = 8.2-8.8$ and $\log(M_{\rm BH}/M_{\odot}) = 7.9-8.3$ respectively. We then fit energy conserving models to the SEDs, obtaining strong constraints on the bolometric luminosity and $\alpha_{\rm OX}$. We compare the SED properties with a larger, X-ray selected AGN sample for both of these scenarios, and observe distinct groupings in spectral shape versus luminosity parameter space. In general, the SED properties are closer to normal if we assume that the low-state is representative. This supports the idea that the large slow outbursts may be due to extrinsic effects (for example microlensing) as opposed to accretion rate changes, but a larger sample of HVAs is needed to be confident of this conclusion.

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.

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.

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.

We analyze four extreme active galactic nuclei (AGN) transients to explore the possibility that they are caused by rare, high-amplitude microlensing events. These previously unknown type-I AGN are located in the redshift range 0.6-1.1 and show changes of > 1.5 mag in the g-band on a time-scale of ~years. Multi-epoch optical spectroscopy, from the William Herschel Telescope, shows clear differential variability in the broad line fluxes with respect to the continuum changes and also evolution in the line profiles. In two cases, a simple point-source, point-lens microlensing model provides an excellent match to the long-term variability seen in these objects. For both models, the parameter constraints are consistent with the microlensing being due to an intervening stellar mass object but as yet there is no confirmation of the presence of an intervening galaxy. The models predict a peak amplification of 10.3/13.5 and an Einstein time-scale of 7.5/10.8 yr, respectively. In one case, the data also allow constraints on the size of the CIII] emitting region, with some simplifying assumptions, to to be ~1.0-6.5 light-days and a lower limit on the size of the MgII emitting region to be > 9 light-days (half-light radii). This CIII] radius is perhaps surprisingly small. In the remaining two objects, there is spectroscopic evidence for an intervening absorber but the extra structure seen in the light curves requires a more complex lensing scenario to adequately explain.

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.

We discuss 76 large amplitude transients (Delta-m>1.5) occurring in the nuclei of galaxies, nearly all with no previously known Active Galactic Nucleus (AGN). They have been discovered as part of the Pan-STARRS1 (PS1) 3pi survey, by comparison with SDSS photometry a decade earlier, and then monitored with the Liverpool Telescope, and studied spectroscopically with the William Herschel Telescope (WHT). Based on colours, light curve shape, and spectra, these transients fall into four groups. A few are misclassified stars or objects of unknown type. Some are red/fast transients and are known or likely nuclear supernovae. A few are either radio sources or erratic variables and so likely blazars. However the majority (~66%) are blue and evolve slowly, on a timescale of years. Spectroscopy shows them to be AGN at z~ 0.3 - 1.4, which must have brightened since the SDSS photometry by around an order of magnitude. It is likely that these objects were in fact AGN a decade ago, but too weak to be recognised by SDSS; they could then be classed as "hypervariable" AGN. By searching the SDSS Stripe 82 quasar database, we find 15 similar objects. We discuss several possible explanations for these slow blue hypervariables - (i) unusually luminous tidal disruption events; (ii) extinction events; (iii) changes in accretion state; and (iv) large amplitude microlensing by stars in foreground galaxies. A mixture of explanations (iii) and (iv) seems most likely. Both hold promise of considerable new insight into the AGN phenomenon.

We present a study of the cataclysmic variable star PT Per based on archival XMM-Newton X-ray data and new optical spectroscopy from the WHT with ISIS. The X-ray data show deep minima which recur at a period of 82 minutes and a hard, unabsorbed X-ray spectrum. The optical spectra of PT Per show a relatively featureless blue continuum. From an analysis of the X-ray and optical data we conclude that PT Per is likely to be a magnetic cataclysmic variable of the polar class in which the minima correspond to those phase intervals when the accretion column rotates out of the field of view of the observer. We suggest that the optical spectrum, obtained around 4 years after the X-ray coverage, is dominated by the white dwarf in the system, implying that PT Per was in a low accretion state at the time of the observations. An analysis of the likely system parameters for PT Per suggests a distance of $\approx90$ pc and a very low-mass secondary, consistent with the idea that PT Per is a "period-bounce" binary.

The Time-Domain Spectroscopic Survey (TDSS) is an SDSS-IV eBOSS subproject primarily aimed at obtaining identification spectra of ~220,000 optically-variable objects systematically selected from SDSS/Pan-STARRS1 multi-epoch imaging. We present a preview of the science enabled by TDSS, based on TDSS spectra taken over ~320 deg^2 of sky as part of the SEQUELS survey in SDSS-III, which is in part a pilot survey for eBOSS in SDSS-IV. Using the 15,746 TDSS-selected single-epoch spectra of photometrically variable objects in SEQUELS, we determine the demographics of our variability-selected sample, and investigate the unique spectral characteristics inherent in samples selected by variability. We show that variability-based selection of quasars complements color-based selection by selecting additional redder quasars, and mitigates redshift biases to produce a smooth quasar redshift distribution over a wide range of redshifts. The resulting quasar sample contains systematically higher fractions of blazars and broad absorption line quasars than from color-selected samples. Similarly, we show that M-dwarfs in the TDSS-selected stellar sample have systematically higher chromospheric active fractions than the underlying M-dwarf population, based on their H-alpha emission. TDSS also contains a large number of RR Lyrae and eclipsing binary stars with main-sequence colors, including a few composite-spectrum binaries. Finally, our visual inspection of TDSS spectra uncovers a significant number of peculiar spectra, and we highlight a few cases of these interesting objects. With a factor of ~15 more spectra, the main TDSS survey in SDSS-IV will leverage the lessons learned from these early results for a variety of time-domain science applications.

We present a systematic search for changing-look quasars based on repeat photometry from SDSS and Pan-STARRS1, along with repeat spectra from SDSS and SDSS-III BOSS. Objects with large, |\Delta g|>1 mag photometric variations in their light curves are selected as candidates to look for changes in broad emission line (BEL) features. Out of a sample of 1011 objects that satisfy our selection criteria and have more than one epoch of spectroscopy, we find 10 examples of quasars that have variable and/or "changing-look" BEL features. Four of our objects have emerging BELs; five have disappearing BELs, and one object shows tentative evidence for having both emerging and disappearing BELs. With redshifts in the range 0.20 < z < 0.63, this sample includes the highest-redshift changing-look quasars discovered to date. We highlight the quasar J102152.34+464515.6 at z = 0.204. Here, not only have the Balmer emission lines strongly diminished in prominence, including H$\beta$ all but disappearing, but the blue continuum $f_{\nu} \propto \nu^{1/3}$ typical of an AGN is also significantly diminished in the second epoch of spectroscopy. Using our selection criteria, we estimate that >15% of strongly variable luminous quasars display changing-look BEL features on rest-frame timescales of 8 to 10 years. Plausible timescales for variable dust extinction are factors of 2-10 too long to explain the dimming and brightening in these sources, and simple dust reddening models cannot reproduce the BEL changes. On the other hand, an advancement such as disk reprocessing is needed if the observed variations are due to accretion rate changes.

We review results from our monitoring observations of several lensed quasars performed in the optical, UV, and X-ray bands. Modeling of the multi-wavelength light curves provides constraints on the extent of the optical, UV, and X-ray emission regions. One of the important results of our analysis is that the optical sizes as inferred from the microlensing analysis are significantly larger than those predicted by the theoretical-thin-disk estimate. In a few cases we also constrain the slope of the size-wavelength relation. Our size constraints of the soft and hard X-ray emission regions of quasars indicate that in some objects of our sample the hard X-ray emission region is more compact than the soft and in others the soft emission region is smaller. This difference may be the result of the relative strengths of the disk-reflected (harder and extended) versus corona-direct (softer and compact) components in the quasars of our sample. Finally, we present the analysis of several strong microlensing events where we detect an evolution of the relativistic Fe line profile as the magnification caustic traverses the accretion disk. These caustic crossings are used to provide constraints on the innermost stable circular orbit (ISCO) radius and the accretion disk inclination angle of the black hole in quasar RX J1131-1231.

We conduct a pilot investigation to determine the optimal combination of color and variability information to identify quasars in current and future multi-epoch optical surveys. We use a Bayesian quasar selection algorithm (Richards et al. 2004) to identify 35,820 type 1 quasar candidates in a 239 square degree field of the Sloan Digital Sky Survey (SDSS) Stripe 82, using a combination of optical photometry and variability. Color analysis is performed on 5-band single- and multi-epoch SDSS optical photometry to a depth of r ~22.4. From these data, variability parameters are calculated by fitting the structure function of each object in each band with a power law model using 10 to >100 observations over timescales from ~1 day to ~8 years. Selection was based on a training sample of 13,221 spectroscopically-confirmed type-1 quasars, largely from the SDSS. Using variability alone, colors alone, and combining variability and colors we achieve 91%, 93%, and 97% quasar completeness and 98%, 98%, and 97% efficiency respectively, with particular improvement in the selection of quasars at 2.7<z<3.5 where quasars and stars have similar optical colors. The 22,867 quasar candidates that are not spectroscopically confirmed reach a depth of i ~22.0; 21,876 (95.7%) are dimmer than coadded i-band magnitude of 19.9, the cut off for spectroscopic follow-up for SDSS on Stripe 82. Brighter than 19.9, we find 5.7% more quasar candidates without confirming spectra in sky regions otherwise considered complete. The resulting quasar sample has sufficient purity (and statistically correctable incompleteness) to produce a luminosity function comparable to those determined by spectroscopic investigations. We discuss improvements that can be made to the process in preparation for performing similar photometric selection and science on data from post-SDSS sky surveys.

We analyze the optical, UV, and X-ray microlensing variability of the lensed quasar SDSS J0924+0219 using six epochs of Chandra data in two energy bands (spanning 0.4-8.0 keV, or 1-20 keV in the quasar rest frame), 10 epochs of F275W (rest-frame 1089A) Hubble Space Telescope data, and high-cadence R-band (rest-frame 2770A) monitoring spanning eleven years. Our joint analysis provides robust constraints on the extent of the X-ray continuum emission region and the projected area of the accretion disk. The best-fit half-light radius of the soft X-ray continuum emission region is between 5x10^13 and 10^15 cm, and we find an upper limit of 10^15 cm for the hard X-rays. The best-fit soft-band size is about 13 times smaller than the optical size, and roughly 7 GM_BH/c^2 for a 2.8x10^8 M_sol black hole, similar to the results for other systems. We find that the UV emitting region falls in between the optical and X-ray emitting regions at 10^14 cm < r_1/2,UV < 3x10^15 cm. Finally, the optical size is significantly larger, by 1.5*sigma, than the theoretical thin-disk estimate based on the observed, magnification-corrected I-band flux, suggesting a shallower temperature profile than expected for a standard disk.

Active galactic nuclei (AGNs) are well-known to exhibit flux variability across a wide range of wavelength regimes, but the precise origin of the variability at different wavelengths remains unclear. To investigate the relatively unexplored near-IR variability of the most luminous AGNs, we conduct a search for variability using well sampled JHKs-band light curves from the 2MASS survey calibration fields. Our sample includes 27 known quasars with an average of 924 epochs of observation over three years, as well as one spectroscopically confirmed blazar (SDSSJ14584479+3720215) with 1972 epochs of data. This is the best-sampled NIR photometric blazar light curve to date, and it exhibits correlated, stochastic variability that we characterize with continuous auto-regressive moving average (CARMA) models. None of the other 26 known quasars had detectable variability in the 2MASS bands above the photometric uncertainty. A blind search of the 2MASS calibration field light curves for AGN candidates based on fitting CARMA(1,0) models (damped-random walk) uncovered only 7 candidates. All 7 were young stellar objects within the {\rho} Ophiuchus star forming region, five with previous X-ray detections. A significant {\gamma}-ray detection (5{\sigma}) for the known blazar using 4.5 years of Fermi photon data is also found. We suggest that strong NIR variability of blazars, such as seen for SDSSJ14584479+3720215, can be used as an efficient method of identifying previously-unidentified {\gamma}-ray blazars, with low contamination from other AGN.

Context. The variations of solar activity over long time intervals using a solar activity reconstruction based on the cosmogenic radionuclide 10Be measured in polar ice cores are studied. Methods. By applying methods of nonlinear dynamics, the solar activity cycle is studied using solar activity proxies that have been reaching into the past for over 9300 years. The complexity of the system is expressed by several parameters of nonlinear dynamics, such as embedding dimension or false nearest neighbors, and the method of delay coordinates is applied to the time series. We also fit a damped random walk model, which accurately describes the variability of quasars, to the solar 10Be data and investigate the corresponding power spectral distribution. The periods in the data series were searched by the Fourier and wavelet analyses. The solar activity on the long-term scale is found to be on the edge of chaotic behavior. This can explain the observed intermittent period of longer lasting solar activity minima. Filtering the data by eliminating variations below a certain period (the periods of 380 yr and 57 yr were used) yields a far more regular behavior of solar activity. A comparison between the results for the 10Be data with the 14C data shows many similarities. Both cosmogenic isotopes are strongly correlated mutually and with solar activity. Finally, we find that a series of damped random walk models provides a good fit to the 10Be data with a fixed characteristic time scale of 1000 years, which is roughly consistent with the quasi-periods found by the Fourier and wavelet analyses.

With upcoming all sky surveys such as LSST poised to generate a deep digital movie of the optical sky, variability-based AGN selection will enable the construction of highly-complete catalogs with minimum contamination. In this study, we generate $g$-band difference images and construct light curves for QSO/AGN candidates listed in SDSS Stripe 82 public catalogs compiled from different methods, including spectroscopy, optical colors, variability, and X-ray detection. Image differencing excels at identifying variable sources embedded in complex or blended emission regions such as Type II AGNs and other low-luminosity AGNs that may be omitted from traditional photometric or spectroscopic catalogs. To separate QSOs/AGNs from other sources using our difference image light curves, we explore several light curve statistics and parameterize optical variability by the characteristic damping timescale ($\tau$) and variability amplitude. By virtue of distinguishable variability parameters of AGNs, we are able to select them with high completeness of 93.4% and efficiency (i.e., purity) of 71.3%. Based on optical variability, we also select highly variable blazar candidates, whose infrared colors are consistent with known blazars. One third of them are also radio detected. With the X-ray selected AGN candidates, we probe the optical variability of X-ray detected optically-extended sources using their difference image light curves for the first time. A combination of optical variability and X-ray detection enables us to select various types of host-dominated AGNs. Contrary to the AGN unification model prediction, two Type II AGN candidates (out of 6) show detectable variability on long-term timescales like typical Type I AGNs. This study will provide a baseline for future optical variability studies of extended sources.

Over the last decade, quasar sample sizes have increased from several thousand to several hundred thousand, thanks mostly to SDSS imaging and spectroscopic surveys. LSST, the next-generation optical imaging survey, will provide hundreds of detections per object for a sample of more than ten million quasars with redshifts of up to about seven. We briefly review optical quasar selection techniques, with emphasis on methods based on colors, variability properties and astrometric behavior.

A damped random walk is a stochastic process, defined by an exponential covariance matrix that behaves as a random walk for short time scales and asymptotically achieves a finite variability amplitude at long time scales. Over the last few years, it has been demonstrated, mostly but not exclusively using SDSS data, that a damped random walk model provides a satisfactory statistical description of observed quasar variability in the optical wavelength range, for rest-frame timescales from 5 days to 2000 days. The best-fit characteristic timescale and asymptotic variability amplitude scale with the luminosity, black hole mass, and rest wavelength, and appear independent of redshift. In addition to providing insights into the physics of quasar variability, the best-fit model parameters can be used to efficiently separate quasars from stars in imaging surveys with adequate long-term multi-epoch data, such as expected from LSST.

We present three complete seasons and two half-seasons of SDSS r-band photometry of the gravitationally lensed quasar SBS 0909+532 from the U.S. Naval Observatory, as well as two seasons each of SDSS g-band and r-band monitoring from the Liverpool Robotic Telescope. Using Monte Carlo simulations to simultaneously measure the system's time delay and model the r-band microlensing variability, we confirm and significantly refine the precision of the system's time delay to \Delta t_{AB} = 50^{+2}_{-4} days, where the stated uncertainties represent the bounds of the formal 1\sigma\ confidence interval. There may be a conflict between the time delay measurement and a lens consisting of a single galaxy. While models based on the Hubble Space Telescope astrometry and a relatively compact stellar distribution can reproduce the observed delay, the models have somewhat less dark matter than we would typically expect. We also carry out a joint analysis of the microlensing variability in the r- and g-bands to constrain the size of the quasar's continuum source at these wavelengths, obtaining log[(r_{s,r}/cm) [cos{i}/0.5]^{1/2}] = 15.3 \pm 0.3 and log[(r_{s,g}/cm) [cos{i}/0.5]^{1/2}] = 14.8 \pm 0.9, respectively. Our current results do not formally constrain the temperature profile of the accretion disk but are consistent with the expectations of standard thin disk theory.

We present 11.2 micron observations of the gravitationally lensed, radio-loud z_s=2.64 quasar MG0414+0534, obtained using the Michelle camera on Gemini North. We find a flux ratio anomaly of A2/A1= 0.93 +/- 0.02 for the quasar images A1 and A2. When combined with the 11.7 micron measurements from Minezaki et al. (2009), the A2/A1 flux ratio is nearly 5-sigma from the expected ratio for a model based on the two visible lens galaxies. The mid-IR flux ratio anomaly can be explained by a satellite (substructure), 0.3" Northeast of image A2, as can the detailed VLBI structures of the jet produced by the quasar. When we combine the mid-IR flux ratios with high-resolution VLBI measurements, we find a best-fit mass between 10^(6.2) and 10^(7.5) M_sol inside the Einstein radius for a satellite substructure modeled as a singular isothermal sphere at the redshift of the main lens (z_l=0.96). We are unable to set an interesting limit on the mass to light ratio due to its proximity to the quasar image A2. While the observations used here were technically difficult, surveys of flux anomalies in gravitational lenses with the James Webb Space Telescope will be simple, fast, and should well constrain the abundance of substructure in dark matter haloes.

We report the discovery by the Palomar Transient Factory (PTF) of the transient source PTF11agg, which is distinguished by three primary characteristics: (1) bright, rapidly fading optical transient emission; (2) a faint, blue quiescent optical counterpart; and (3) an associated year-long, scintillating radio transient. We argue that these observed properties are inconsistent with any known class of Galactic transients, and instead suggest a cosmological origin. The detection of incoherent radio emission at such distances implies a large emitting region, from which we infer the presence of relativistic ejecta. The observed properties are all consistent with the population of long-duration gamma-ray bursts (GRBs), marking the first time such an outburst has been discovered in the distant universe independent of a high-energy trigger. We searched for possible high-energy counterparts to PTF11agg, but found no evidence for associated prompt emission. We therefore consider three possible scenarios to account for a GRB-like afterglow without a high-energy counterpart: an "untriggered" GRB (lack of satellite coverage), an "orphan" afterglow (viewing-angle effects), and a "dirty fireball" (suppressed high-energy emission). The observed optical and radio light curves appear inconsistent with even the most basic predictions for off-axis afterglow models. The simplest explanation, then, is that PTF11agg is a normal, on-axis long-duration GRB for which the associated high-energy emission was simply missed. However, we have calculated the likelihood of such a serendipitous discovery by PTF and find that it is quite small (~ 2.6%). While not definitive, we nonetheless speculate that PTF11agg may represent a new, more common (> 4 times the on-axis GRB rate at 90% confidence) class of relativistic outbursts lacking associated high-energy emission.

We present the Data Release 9 Quasar (DR9Q) catalog from the Baryon Oscillation Spectroscopic Survey (BOSS) of the Sloan Digital Sky Survey III. The catalog includes all BOSS objects that were targeted as quasar candidates during the survey, are spectrocopically confirmed as quasars via visual inspection, have luminosities Mi[z=2]<-20.5 (in a $\Lambda$CDM cosmology with H0 = 70 km/s/Mpc, $\Omega_{\rm M}$ = 0.3, and $\Omega_{\Lambda}$ = 0.7) and either display at least one emission line with full width at half maximum (FWHM) larger than 500 km/s or, if not, have interesting/complex absorption features. It includes as well, known quasars (mostly from SDSS-I and II) that were reobserved by BOSS. This catalog contains 87,822 quasars (78,086 are new discoveries) detected over 3,275 deg$^{2}$ with robust identification and redshift measured by a combination of principal component eigenspectra newly derived from a training set of 8,632 spectra from SDSS-DR7. The number of quasars with $z>2.15$ (61,931) is ~2.8 times larger than the number of z>2.15 quasars previously known. Redshifts and FWHMs are provided for the strongest emission lines (CIV, CIII], MgII). The catalog identifies 7,533 broad absorption line quasars and gives their characteristics. For each object the catalog presents five-band (u,g,r,i,z) CCD-based photometry with typical accuracy of 0.03 mag, and information on the morphology and selection method. The catalog also contains X-ray, ultraviolet, near-infrared, and radio emission properties of the quasars, when available, from other large-area surveys.

We investigate the use of optical photometric variability to select and identify blazars in large-scale time-domain surveys, in part to aid in the identification of blazar counterparts to the ~30% of gamma-ray sources in the Fermi 2FGL catalog still lacking reliable associations. Using data from the optical LINEAR asteroid survey, we characterize the optical variability of blazars by fitting a damped random walk model to individual light curves with two main model parameters, the characteristic timescales of variability (tau), and driving amplitudes on short timescales (sigma). Imposing cuts on minimum tau and sigma allows for blazar selection with high efficiency E and completeness C. To test the efficacy of this approach, we apply this method to optically variable LINEAR objects that fall within the several-arcminute error ellipses of gamma-ray sources in the Fermi 2FGL catalog. Despite the extreme stellar contamination at the shallow depth of the LINEAR survey, we are able to recover previously-associated optical counterparts to Fermi AGN with E > 88% and C = 88% in Fermi 95% confidence error ellipses having semimajor axis r < 8'. We find that the suggested radio counterpart to Fermi source 2FGL J1649.6+5238 has optical variability consistent with other gamma-ray blazars, and is likely to be the gamma-ray source. Our results suggest that the variability of the non-thermal jet emission in blazars is stochastic in nature, with unique variability properties due to the effects of relativistic beaming. After correcting for beaming, we estimate that the characteristic timescale of blazar variability is ~3 years in the rest-frame of the jet, in contrast with the ~320 day disk flux timescale observed in quasars. The variability-based selection method presented will be useful for blazar identification in time-domain optical surveys, and is also a probe of jet physics.

We provide a quantitative description and statistical interpretation of the optical continuum variability of quasars. The Sloan Digital Sky Survey (SDSS) has obtained repeated imaging in five UV-to-IR photometric bands for 33,881 spectroscopically confirmed quasars. About 10,000 quasars have an average of 60 observations in each band obtained over a decade along Stripe 82 (S82), whereas the remaining ~25,000 have 2-3 observations due to scan overlaps. The observed time lags span the range from a day to almost 10 years, and constrain quasar variability at rest-frame time lags of up to 4 years, and at rest-frame wavelengths from 1000A to 6000A. We publicly release a user-friendly catalog of quasars from the SDSS Data Release 7 that have been observed at least twice in SDSS or once in both SDSS and the Palomar Observatory Sky Survey, and we use it to analyze the ensemble properties of quasar variability. Based on a damped random walk (DRW) model defined by a characteristic time scale and an asymptotic variability amplitude that scale with the luminosity, black hole mass, and rest wavelength for individual quasars calibrated in S82, we can fully explain the ensemble variability statistics of the non-S82 quasars such as the exponential distribution of large magnitude changes. All available data are consistent with the DRW model as a viable description of the optical continuum variability of quasars on time scales of ~5-2000 days in the rest frame. We use these models to predict the incidence of quasar contamination in transient surveys such as those from PTF and LSST.

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.

Periodic light curve behavior predicted for some binary black hole systems might be detected in large samples, such as the multi-million quasar sample expected from the Large Synoptic Survey Telescope (LSST). We investigate the false-alarm probability for the discovery of a periodic signal in light curves simulated using damped random walk (DRW) model. This model provides a good description of observed light curves, and does not include periodic behavior. We used the Lomb-Scargle periodogram to search for a periodic signal in a million simulated light curves that properly sample the DRW parameter space, and the LSST cadence space. We find that even a very conservative threshold for the false-alarm probability still yields thousands of "good" binary black hole candidates. We conclude that the future claims for binary black holes based on Lomb-Scargle analysis of LSST light curves will have to be interpreted with caution.

We develop a method for separating quasars from other variable point sources using SDSS Stripe 82 light curve data for ~10,000 variable objects. To statistically describe quasar variability, we use a damped random walk model parametrized by a damping time scale, tau, and an asymptotic amplitude (structure function), SF_inf. With the aid of an SDSS spectroscopically confirmed quasar sample, we demonstrate that variability selection in typical extragalactic fields with low stellar density can deliver complete samples with reasonable purity (or efficiency, E). Compared to a selection method based solely on the slope of the structure function, the inclusion of the tau information boosts E from 60% to 75% while maintaining a highly complete sample (98%) even in the absence of color information. For a completeness of C=90%, E is boosted from 80% to 85%. Conversely, C improves from 90% to 97% while maintaining E=80% when imposing a lower limit on tau. With the aid of color selection, the purity can be further boosted to 96%, with C= 93%. Hence, selection methods based on variability will play an important role in the selection of quasars with data provided by upcoming large sky surveys, such as Pan-STARRS and the Large Synoptic Survey Telescope (LSST). For a typical (simulated) LSST cadence over 10 years and a photometric accuracy of 0.03 mag (achieved at i~22), C is expected to be 88% for a simple sample selection criterion of tau>100 days. In summary, given an adequate survey cadence, photometric variability provides an even better method than color selection for separating quasars from stars.

We model the time variability of ~9,000 spectroscopically confirmed quasars in SDSS Stripe 82 as a damped random walk. Using 2.7 million photometric measurements collected over 10 years, we confirm the results of Kelly et al. (2009) and Koz{\l}owski et al. (2010) that this model can explain quasar light curves at an impressive fidelity level (0.01-0.02 mag). The damped random walk model provides a simple, fast [O(N) for N data points], and powerful statistical description of quasar light curves by a characteristic time scale (tau) and an asymptotic rms variability on long time scales (SF_inf). We searched for correlations between these two variability parameters and physical parameters such as luminosity and black hole mass, and rest-frame wavelength. We find that tau increases with increasing wavelength with a power law index of 0.17, remains nearly constant with redshift and luminosity, and increases with increasing black hole mass with power law index of 0.21+/-0.07. The amplitude of variability is anti-correlated with the Eddington ratio, which suggests a scenario where optical fluctuations are tied to variations in the accretion rate. The radio-loudest quasars have systematically larger variability amplitudes by about 30%, when corrected for the other observed trends, while the distribution of their characteristic time scale is indistinguishable from that of the full sample. We do not detect any statistically robust differences in the characteristic time scale and variability amplitude between the full sample and the small subsample of quasars detected by ROSAT. Our results provide a simple quantitative framework for generating mock quasar light curves, such as currently used in LSST image simulations. (abridged)

We present radio and X-ray observations, as well as optical light curves, for a subset of 26 BL Lac candidates from the Sloan Digital Sky Survey (SDSS) lacking strong radio emission and with z<2.2. Half of these 26 objects are shown to be stars, galaxies, or absorbed quasars. We conclude that the other 13 objects are Active Galactic Nuclei (AGN) with abnormally weak emission features; ten of those 13 are definitively radio-quiet, and, for those with available optical light curves, their level of optical flux variability is consistent with radio-quiet quasars. We cannot exclude the possibility that some of these 13 AGN lie on the extremely radio-faint tail of the BL Lac distribution, but our study generally supports the notion that all BL Lac objects are radio-loud. These radio-quiet AGN appear to have intrinsically weak or absent broad emission line regions, and, based on their X-ray properties, we argue that some are low-redshift analogs to weak line quasars (WLQs). SDSS BL Lac searches are so far the only systematic surveys of the SDSS database capable of recovering such exotic low-redshift WLQs. There are 71 more z<2.2 radio-quiet BL Lac candidates already identified in the SDSS not considered here, and many of those might be best unified with WLQs as well. Future studies combining low- and high-redshift WLQ samples will yield new insight on our understanding of the structure and formation of AGN broad emission line regions.

We present the first resolved mid-IR (11 micron) observations of the four-image quasar lens H1413+117 using the Michelle camera on Gemini North. All previous observations (optical, near-IR, and radio) of this lens show a "flux anomaly," where the image flux ratios cannot be explained by a simple, central lens galaxy. We attempt to reproduce the mid-IR flux ratios, which are insensitive to extinction and microlensing, by modeling the main lens as a singular isothermal ellipsoid. This model fails to reproduce the flux ratios. However, we can explain the flux ratios simply by adding to the model a nearby galaxy detected in the H-band by HST/NICMOS-NIC2. This perturbing galaxy lies 4.0" from the main lens and it has a critical radius of 0.63" +/- 0.02" which is similar to that of the main lens, as expected from their similar H-band fluxes. More remarkably, this galaxy is not required to obtain a good fit to the system astrometry, so this represents the first clear detection of an object through its effect on the image fluxes of a gravitational lens. This is a parallel to the detections of visible satellites from astrometric anomalies, and provides a proof of the concept of searching for substructure in galaxies using anomalous flux ratios.

Significant progress in the description of quasar variability has been recently made by employing SDSS and POSS data. Common to most studies is a fundamental assumption that photometric observations at two epochs for a large number of quasars will reveal the same statistical properties as well-sampled light curves for individual objects. We critically test this assumption using light curves for a sample of $\sim$2,600 spectroscopically confirmed quasars observed about 50 times on average over 8 years by the SDSS stripe 82 survey. We find that the dependence of the mean structure function computed for individual quasars on luminosity, rest-frame wavelength and time is qualitatively and quantitatively similar to the behavior of the structure function derived from two-epoch observations of a much larger sample. We also reproduce the result that the variability properties of radio and X-ray selected subsamples are different. However, the scatter of the variability structure function for fixed values of luminosity, rest-frame wavelength and time is similar to the scatter induced by the variance of these quantities in the analyzed sample. Hence, our results suggest that, although the statistical properties of quasar variability inferred using two-epoch data capture some underlying physics, there is significant additional information that can be extracted from well-sampled light curves for individual objects.

Measured gravitational waveforms from black hole binary inspiral events directly determine absolute luminosity distances. To use these data for cosmology, it is necessary to independently obtain redshifts for the events, which may be difficult for those without electromagnetic counterparts. Here it is demonstrated that certainly in principle, and possibly in practice, clustering of galaxies allows extraction of the redshift information from a sample statistically for the purpose of estimating mean cosmological parameters, without identification of host galaxies for individual events. We extract mock galaxy samples from the 6th Data Release of the Sloan Digital Sky Survey resembling those that would be associated with inspiral events of stellar mass black holes falling into massive black holes at redshift z ~ 0.1 to 0.5. A simple statistical procedure is described to estimate a likelihood function for the Hubble constant H_0: each galaxy in a LISA error volume contributes linearly to the log likelihood for the source redshift, and the log likelihood for each source contributes linearly to that of H_0. This procedure is shown to provide an accurate and unbiased estimator of H_0. It is estimated that a precision better than one percent in H_0 may be possible if the rate of such events is sufficiently high, on the order of 20 to z = 0.5.

In 2000, Lamb and Reichart predicted that gamma-ray bursts (GRBs) and their afterglows occur in sufficient numbers and at sufficient brightnesses at very high redshifts (z > 5) to eventually replace quasars as the preferred probe of element formation and reionization in the early universe and to be used to characterize the star-formation history of the early universe, perhaps back to when the first stars formed. Here we report the discovery of the afterglow of GRB 050904 and the identification of GRB 050904 as the first very high redshift GRB. We measure its redshift to be 6.39(+0.11,-0.12), which is consistent with the reported spectroscopic redshift (6.29 +/- 0.01). Furthermore, just redward of Ly-alpha the flux is suppressed by a factor of three on the first night, but returns to expected levels by the fourth night. We propose that this is due to absorption by molecular hydrogen that was excited to rovibrational states by the GRB's prompt emission, but was then overtaken by the jet. Now that very high redshift GRBs have been shown to exist, and at least in this case the afterglow was very bright, observing programs that are designed to capitalize on this science will likely drive a new era of study of the early universe, using GRBs as probes.