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

We present Chandra observations of 2106 radio-quiet quasars in the redshift range 1.7<z<2.7 from the Sloan Digital Sky Survey (SDSS), through data release fourteen (DR14), that do not contain broad absorption lines (BAL) in their rest-frame UV spectra. This sample adds over a decade worth of SDSS and Chandra observations to our previously published sample of 139 quasars from SDSS DR5 which is still used to correlate X-ray and optical/UV emission in typical quasars. We fit the SDSS spectra for 753 of the quasars in our sample that have high-quality (exposure time $\gtrapprox$10 ks and off-axis observation angle <10 arcmin) X-ray observations, and analyze their X-ray-to-optical SED properties ($\alpha_{ox}$ and $\Delta\alpha_{ox}$) with respect to the measured CIV and MgII emission-line rest-frame equivalent width (EW) and the CIV emission-line blueshift. We find significant correlations (at the >99.99% level) between $\alpha_{ox}$ and these emission-line parameters, as well as between $\Delta\alpha_{ox}$ and CIV EW. Slight correlations are found between $\Delta\alpha_{ox}$ and CIV blueshift, MgII EW, and the CIV EW to MgII EW ratio. The best-fit trend in each parameter space is used to compare the X-ray weakness ($\Delta\alpha_{ox}$) and optical/UV emission properties of typical quasars and weak-line quasars (WLQs). The WLQs typically exhibit weaker X-ray emission than predicted by the typical quasar relationships. The best-fit relationships for our typical quasars are consistent with predictions from the disk-wind quasar model. The behavior of the WLQs compared to our typical quasars can be explained by an X-ray "shielding" model.

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

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

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

We investigate the connection between the intrinsic C IV absorption line variability and the continuum flux changes of broad absorption line (BAL) quasars using a sample of 78 sources in the Stripe 82 region. The absorption trough variability parameters are measured using the archival multi-epoch spectroscopic data from the Sloan Digital Sky Survey (SDSS), and the continuum flux variability parameters are estimated from the photometric light curves obtained by the SDSS and the Catalina Real-Time Survey (CRTS) surveys. We find evidence for weak correlations (\rho ~ 0.3) between the intrinsic C IV absorption line variability and the quasar continuum variability for the final sample of 78 BAL quasars. The correlation strengths improve (\rho ~0.5) for the "high-SNR" sample sources that have higher spectral signal-to-noise ratio. Using two sub-sets of the high-SNR sample differing on the absorption trough depth, we find that the shallow trough sub-set shows an even stronger correlation (\rho ~ 0.6), whereas the deep trough sub-set does not show any correlation between the absorption line variability and the continuum variability. These results point to the important role of saturation effects in the correlation between the absorption line variability and the continuum variability of BAL quasars. Considering other effects that can also smear the correlation, we conclude that the actual correlation between the absorption line and continuum variability is even stronger.

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.

We report the discovery of disappearance of Mg ii, Al iii, C iv, and Si iv broad absorption lines (BALs) at the same velocity, accompanied by a new Civ BAL emerging at a higher velocity, in the quasar J0827+4252 at z = 2.038. This is the first report of BAL disappearance (i) over Mg ii, Al iii, C iv, and Si iv ions and (ii) in a weak emission-line quasar (WLQ). The discovery is based on four spectra from the SDSS and one follow-up spectrum from HET/LRS2. The simultaneous Civ BAL disappearance and emergence at different velocities, together with no variations in the CRTS light curve, indicate that ionization changes in the absorbing material are unlikely to cause the observed BAL variability. Our analyses reveal that transverse motion is the most likely dominant driver of the BAL disappearance/emergence. Given the presence of mildly relativistic BAL outflows and an apparently large C iv emission-line blueshift that is likely associated with strong bulk outflows in this WLQ, J0827+4252 provides a notable opportunity to study extreme quasar winds and their potential in expelling material from inner to large-scale regions.

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

Using 16,421 spectra from a sample of 340 quasars ($1.62<z<3.30$) from the SDSS Reverberation Mapping Project, we present an analysis of quasar spectral variability. We confirm the intrinsic Baldwin Effect and brighter-means-bluer trends in which emission line strength and color are associated with changes in luminosity. We then create a composite differential spectrum that captures the wavelength dependence of quasar variability as a function of change in luminosity. When using a bandpass around 1740 \AA\ to describe the luminosity, the differential spectrum follows a power law at longer wavelengths that breaks blueward of 1700 \AA. However, the shape of the differential spectrum, the location of the power law break, and the slope of the intrinsic Baldwin Effect all vary with the choice of bandpass used to define the change in quasar luminosity. We propose that the observed behavior can be explained by inhomogeneous accretion or slim accretion disk models where delays in the reprocessing of incident light on the accretion disk cause local deviations in temperature from the thin disk model. Finally, we quantify the effects on cosmology studies due to the variations in the quasar spectrum in the Lyman-$\alpha$ forest wavelength range. Using the observed spectroscopic signatures to predict the quasar continuum over the interval $1040 < \lambda < 1200$ \AA, we find that the derived spectral templates can reduce the uncertainty of the Lyman-$\alpha$ forest continuum level in individual epochs from 17.2% to 7.7%, near the level where systematic errors in SDSS flux calibration are expected to dominate.

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

We report the discovery of rapid variations of a high-velocity Mg ii broad absorption line (BAL) trough in the quasar SDSS J133356.02+001229.1 (zem ~ 0.9197). Vivek et al. (2012a) revealed the emergence and subsequent near disappearance of a BAL component in this source having an ejection velocity of ~ 28000 km s-1. Our further follow up studies with South African Large Telescope (SALT) reveal the dramatic nature of the absorption line variability in this source. The absorption line emerged again at the same velocity and nearly disappeared within the SALT observations. Our observations allow us to probe variability over timescales of the order of few days to 4.2 years in the QSO rest-frame. The observed velocity stability of BAL absorption does not point to any line of sight acceleration/deceleration of BAL clouds. The ionization parameter of the absorbing cloud is constrained from the column density ratio of Mg II to Fe II ground state absorption. In the absence of strong optical continuum variability, we suggest that photoionization driven BAL variability due to changes in the shielding, multiple streaming clouds across our line of sight in a co-rotating wind or a combination of both as possible explanations for the observed strong equivalent width variations.

We present a detailed study of the quasar-galaxy pair: J1243+4043 - UGC07904. The sight line of the background quasar ( $z_q$= 1.5266) passes through a region of the galaxy ($z_g$=0.0169) at an impact parameter of 6.9 kpc with high metallicity (0.5 Z$_\odot$) and negligible dust extinction. We detect HI 21-cm absorption from the foreground galaxy at arcsecond and milliarcsecond scales. For typical cold neutral medium (CNM) temperatures in the Milky Way, this 21-cm absorber can be classified as a damped Ly$\alpha$ absorber (DLA). We infer the harmonic mean spin temperature of the gas to be $\sim$400 K and for a simple two-phase medium we estimate the CNM-fraction to be $f_{\rm CNM}$ = 0.27. This is remarkably consistent with the CNM fraction observed in the Galaxy and less than that of high-redshift DLAs. The quasar exhibits a core-jet morphology on milliarcsecond scales, corresponding to an overall extent of $\sim$9 pc at $z_g$. We show that the size of CNM absorbing clouds associated with the foreground galaxy is $>$5 pc and they may be part of cold gas structures that extend beyond $\sim$35 pc. Interestingly, the rotation measure of quasar J1243+4043 is higher than any other source in samples of quasars with high-$z$ DLAs. However, we do not find any detectable differences in RMs and polarization fraction of sight lines with or without high-$z$ ($z\ge2$) DLAs or low-$z$ ($z\le0.3$) 21-cm absorbers. Finally, the foreground galaxy UGC07904 is also part of a galaxy group. We serendipitously detect HI 21-cm emission from four members of the group, and a $\sim$80 kpc long HI bridge connecting two of the other members. The latter, together with the properties of the group members, suggests that the group is a highly interactive environment.

We present measurements of the Baryon Acoustic Oscillation (BAO) scale in redshift-space using the clustering of quasars. We consider a sample of 147,000 quasars from the extended Baryon Oscillation Spectroscopic Survey (eBOSS) distributed over 2044 square degrees with redshifts $0.8 < z < 2.2$ and measure their spherically-averaged clustering in both configuration and Fourier space. Our observational dataset and the 1400 simulated realizations of the dataset allow us to detect a preference for BAO that is greater than 2.8$\sigma$. We determine the spherically averaged BAO distance to $z = 1.52$ to 3.8 per cent precision: $D_V(z=1.52)=3843\pm147 \left(r_{\rm d}/r_{\rm d, fid}\right)\ $Mpc. This is the first time the location of the BAO feature has been measured between redshifts 1 and 2. Our result is fully consistent with the prediction obtained by extrapolating the Planck flat $\Lambda$CDM best-fit cosmology. All of our results are consistent with basic large-scale structure (LSS) theory, confirming quasars to be a reliable tracer of LSS, and provide a starting point for numerous cosmological tests to be performed with eBOSS quasar samples. We combine our result with previous, independent, BAO distance measurements to construct an updated BAO distance-ladder. Using these BAO data alone and marginalizing over the length of the standard ruler, we find $\Omega_{\Lambda} > 0$ at 6.6$\sigma$ significance when testing a $\Lambda$CDM model with free curvature.

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.

Deep galaxy surveys have revealed that the global star formation rate (SFR) density in the Universe peaks at 1 < z < 2 and sharply declines towards z = 0. But a clear picture of the underlying processes, in particular the evolution of cold atomic (~100 K) and molecular gas phases, that drive such a strong evolution is yet to emerge. MALS is designed to use MeerKAT's L- and UHF-band receivers to carry out the most sensitive (N(HI)>10$^{19}$ cm$^{-2}$) dust-unbiased search of intervening HI 21-cm and OH 18-cm absorption lines at 0 < z < 2. This will provide reliable measurements of the evolution of cold atomic and molecular gas cross-sections of galaxies, and unravel the processes driving the steep evolution in the SFR density. The large sample of HI and OH absorbers obtained from the survey will (i) lead to tightest constraints on the fundamental constants of physics, and (ii) be ideally suited to probe the evolution of magnetic fields in disks of galaxies via Zeeman Splitting or Rotation Measure synthesis. The survey will also provide an unbiased census of HI and OH absorbers, i.e. cold gas associated with powerful AGNs (>10$^{24}$ W Hz$^{-1}$) at 0 < z < 2, and will simultaneously deliver a blind HI and OH emission line survey, and radio continuum survey. Here, we describe the MALS survey design, observing plan and the science issues to be addressed under various science themes.

We report the results of our radio, optical and infra-red studies of a peculiar radio source 4C~35.06, an extended radio-loud AGN at the center of galaxy cluster Abell 407 ($z=0.047$). The central region of this cluster hosts a remarkably tight ensemble of nine galaxies, the spectra of which resemble those of passive red ellipticals, embedded within a diffuse stellar halo of $\sim$1~arcmin size. This system (named the `Zwicky's Nonet') provides unique and compelling evidence for a multiple-nucleus cD galaxy precursor. Multifrequency radio observations of 4C~35.06 with the Giant Meterwave Radio Telescope (GMRT) at 610, 235 and 150 MHz reveal a system of 400~kpc scale helically twisted and kinked radio jets and outer diffuse lobes. The outer extremities of jets contain extremely steep spectrum (spectral index -1.7 to -2.5) relic/fossil radio plasma with a spectral age of a few$\,\times (10^7 - 10^8)$ yr. Such ultra-steep spectrum relic radio lobes without definitive hot-spots are rare, and they provide an opportunity to understand the life-cycle of relativistic jets and physics of black hole mergers in dense environments. We interpret our observations of this radio source in the context of the growth of its central black hole, triggering of its AGN activity and jet precession, all possibly caused by galaxy mergers in this dense galactic system. A slow conical precession of the jet axis due to gravitational perturbation between interacting black holes is invoked to explain the unusual jet morphology.

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 report on the diversity in quasar spectra from the Baryon Oscillation Spectroscopic Survey. After filtering the spectra to mitigate selection effects and Malmquist bias associated with a nearly flux-limited sample, we create high signal-to-noise ratio composite spectra from 58,656 quasars (2.1 \le z \le 3.5), binned by luminosity, spectral index, and redshift. With these composite spectra, we confirm the traditional Baldwin effect (BE, i.e., the anticorrelation of C IV equivalent width (EW) and luminosity) that follows the relation W_\lambda \propto L^{\beta_w} with slope \beta_w = -0.35 \pm 0.004, -0.35 \pm 0.005, and -0.41 \pm 0.005 for z = 2.25, 2.46, and 2.84, respectively. In addition to the redshift evolution in the slope of the BE, we find redshift evolution in average quasar spectral features at fixed luminosity. The spectroscopic signature of the redshift evolution is correlated at 98% with the signature of varying luminosity, indicating that they arise from the same physical mechanism. At a fixed luminosity, the average C IV FWHM decreases with increasing redshift and is anti-correlated with C IV EW. The spectroscopic signature associated with C IV FWHM suggests that the trends in luminosity and redshift are likely caused by a superposition of effects that are related to black hole mass and Eddington ratio. The redshift evolution is the consequence of a changing balance between these two quantities as quasars evolve toward a population with lower typical accretion rates at a given black hole mass.

We describe the redmonster automated redshift measurement and spectral classification software designed for the extended Baryon Oscillation Spectroscopic Survey (eBOSS) of the Sloan Digital Sky Survey IV (SDSS-IV). We describe the algorithms, the template standard and requirements, and the newly developed galaxy templates to be used on eBOSS spectra. We present results from testing on early data from eBOSS, where we have found a 90.5% automated redshift and spectral classification success rate for the luminous red galaxy sample (redshifts 0.6$\lesssim z \lesssim$1.0). The redmonster performance meets the eBOSS cosmology requirements for redshift classification and catastrophic failures, and represents a significant improvement over the previous pipeline. We describe the empirical processes used to determine the optimum number of additive polynomial terms in our models and an acceptable $\Delta\chi_r^2$ threshold for declaring statistical confidence. Statistical errors on redshift measurement due to photon shot noise are assessed, and we find typical values of a few tens of km s$^{-1}$. An investigation of redshift differences in repeat observations scaled by error estimates yields a distribution with a Gaussian mean and standard deviation of $\mu\sim$0.01 and $\sigma\sim$0.65, respectively, suggesting the reported statistical redshift uncertainties are over-estimated by $\sim$54%. We assess the effects of object magnitude, signal-to-noise ratio, fiber number, and fiber head location on the pipeline's redshift success rate. Finally, we describe directions of ongoing development.

The Baryon Oscillation Spectroscopic Survey (BOSS) has collected more than 150,000 $2.1 \leq z \leq 3.5$ quasar spectra since 2009. Using this unprecedented sample, we create a composite spectrum in the rest-frame of 102,150 quasar spectra from 800 \AA\ to 3300 \AA\ at a signal-to-noise ratio close to 1000 per pixel ($\Delta v$ of 69 km~s$^{-1}$). Included in this analysis is a correction to account for flux calibration residuals in the BOSS spectrophotometry. We determine the spectral index as a function of redshift of the full sample, warp the composite spectrum to match the median spectral index, and compare the resulting spectrum to SDSS photometry used in target selection. The quasar composite matches the color of the quasar population to within 0.02 magnitudes in $g-r$, 0.03 magnitudes in $r-i$, and 0.01 magnitudes in $i-z$ over the redshift range $2.2<z<2.6$. The composite spectrum deviates from the imaging photometry by 0.05 magnitudes around $z = 2.7$, likely due to differences in target selection as the quasar colors become similar to the stellar locus at this redshift. Finally, we characterize the line features in the high signal-to-noise composite and identify nine faint lines not found in the previous composite spectrum from SDSS.

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

We study the transient (i.e. emerging or disappearing) C IV broad absorption line (BAL) components in 50 radio detected QSOs using multi-epoch spectra available in Sloan Digital Sky Survey DR10. We report the detectionof 6 BALQSOs having at least one distinct transient C IV absorption component. Based on the structure function analysis of optical light curves, we suggest that the transient absorption is unlikely to be triggered by continuum variations. Transient absorption components usually have low C IV equivalent widths (< 8 \AA), high ejection velocities (> 10000 \kms) and typically occur over rest-frame timescales > 800 days. The detection rate of transient C IV absorption seen in our sample is higher than that reported in the literature. Using a control sample of QSOs, we show that this difference is most likely due to the longer monitoring time-scale of sources in our sample while the effect of small number statistics cannot be ignored. Thus, in order to establish the role played by radio jets in driving the BAL outflows, we need a larger sample of radio detected BALs monitored over more than 3 years in the QSO's rest frame. We also find that the transient phenomenon in radio detected and radio quiet BALs does not depend on any of the QSO properties i.e. the Eddington ratio, black hole mass, bolometric luminosity and optical-to-IR colours. All this suggests that transient BAL phenomenon is simply the extreme case of BAL variability.

Using long-slit optical spectra obtained with the 2-m telescope at IUCAA Girawali Observatory, we show that the radio source J094221.98+062335.2 (z=0.123) is associated with a galaxy pair undergoing a major merger. Its companion is a normal star-forming galaxy infalling with a velocity of 185 km/s at a projected separation of 4.8 kpc. Using the Westerbork Synthesis Radio Telescope (WSRT) and Giant Metrewave Radio Telescope (GMRT) we detect a strong HI 21-cm absorption at the systemic redshift of the radio galaxy with N(HI)~9x10^21 cm^{-2} for an assumed spin-temperature of 100 K. Such a strong HI 21-cm absorption is rare and has been seen only in a few compact radio sources associated with similar merging galaxy pairs. Milliarcsecond resolution Very Long Baseline Array (VLBA) observations resolve the radio source into a compact symmetric object with the hotspot separation of 89 pc. The 21-cm absorption is detected in the VLBA spectra towards both the radio lobes albeit with a strong optical depth gradient. We show that the strong 21-cm absorption is consistent with it being arising from a clumpy circumnuclear disk/torus. We also detect two weaker absorption lines redshifted with respect to the radio source in the WSRT/GMRT spectrum. They probably represent cold (i.e. T< 10^4 K) HI gas falling into the radio source. The presence of high concentration of HI gas in the circumnuclear regions and signature of infalling cold gas allows us to conjecture that the young radio source may have been triggered by the gas infall due to the ongoing merger.

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

We present results of our time variability studies of Mg II and Al III absorption lines in a sample of 22 Low Ionization Broad Absorption Line QSOs (LoBAL QSOs) at 0.2 <= zem <= 2.1 using the 2m telescope at IUCAA Girawali Observatory over a time-scale of 10 days to 7.69 years in the QSO's rest frame. Spectra are analysed in conjunction with photometric light curves from Catalina Real-Time Transient Survey. Long time-scale (i.e >= 1 year) absorption line variability is seen in 8 cases (36% systems) while only 4 of them (i.e 18% systems) show variability over short time-scales (i.e < 1 year). We notice a tendency of highly variable LoBAL QSOs to have high ejection velocity, low equivalent width and low redshift. The detection rate of variability in LoBAL QSOs showing Fe fine-structure lines (FeLoBAL QSOs) is less than that seen in non-Fe LoBAL QSOs. Absorption line variability is more frequently detected in QSOs having continuum dominated by Fe emission lines compared to rest of the QSOs. Confirming these trends with a bigger sample will give vital clues for understanding the physical distinction between different BAL QSO sub-classes. We correlate the absorption line variability with various parameters derived from continuum light curves and find no clear correlation between continuum flux and absorption line variabilities. However, sources with large absorption line variability also show large variability in their light curves. We also see appearance/disappearance of absorption components in 2 cases and clear indications for profile variations in 4 cases. The observed variability can be best explained by a combination of process driven by continuum variations and clouds transiting across the line of sight.

We study the time variability of five Fe Low ionization Broad Absorption Line (FeLoBAL) QSOs using repeated spectroscopic observations with the 2m telescope at IUCAA Girawali observatory (IGO) spanning an interval of upto 10 years. We report a dramatic variation in Al III and Fe III fine-structure lines in the spectra of SDSS J221511.93-004549.9 (z_em ~ 1.478). However, there is no such strong variability shown by the C IV absorption. This source is known to be unusual with (i) the continuum emission dominated by Fe emission lines, (ii) Fe III absorption being stronger than Fe II and (iii) the apparent ratio of Fe III UV 48 to Fe III UV 34 absorption suggesting an inverted population ratio. This is the first reported detection of time variability in the Fe III fine-structure lines in QSO spectra. There is a strong reduction in the absorption strength of these lines between year 2000 and 2008. Using the template fitting techniques, we show that the apparent inversion of strength of UV lines could be related to the complex spectral energy distribution of this QSO. The observed variability can be related to change in the ionization state of the gas or due to transverse motion of this absorbing gas. The shortest variability timescale of Al III line gives a lower limit on the electron density of the absorbing gas as n_e >= 1.1 x 10^4 cm^-3. The remaining 4 FeLoBALs do not show any changes beyond the measurement uncertainties either in optical depth or in the velocity structure. We present the long-term photometric light curve for all of our sources. Among them only SDSS J221511.93-004549.9 shows significant (>= 0.2 mag) variability.

We report a dynamically evolving low ionization broad absorption line flow in the QSO SDSS J133356.02+001229.1 (at z_em = 0.9197). These observations are part of our ongoing monitoring of low ionization broad absorption line (BAL) QSOs with the 2m telescope at IUCAA Girawali observatory (IGO). The broad Mg II absorption with an ejection velocity of 1.7x10^4 km/s, found in the Sloan Digital Sky Survey (SDSS) spectra, has disappeared completely in our IGO spectra. We found an emerging new component at an ejection velocity of 2.8 x 10^4 km/s. During our monitoring period this component has shown strong evolution both in its velocity width and optical depth and nearly disappeared in our latest observations. Acceleration of a low velocity component seen in SDSS spectrum to a higher velocity is unlikely as the Mg II column densities are always observed to be higher for the new component. We argue that the observed variations may not be related to ionization changes and are consistent with absorption produced by multi-streaming flow transiting across our line of sight. We find a possible connection between flux variation of the QSO and N(Mg II) of the newly emerged component. This could mean the ejection being triggered by changes in the accretion disk or dust reddening due to the outflowing gas.

We report long-slit spectroscopic observations of SDSS J092712+294344 carried-out at the recently commissioned 2m telescope in IUCAA Girawali Observatory, India. This AGN-like source is known to feature three sets of emission lines at zem = 0.6972, 0.7020 and 0.7128. Different scenarios such as a recoiling black hole after asymmetric emission of gravitational waves, binary black holes and possible merging systems are proposed for this object. We test these scenarios by comparing our spectra with that from the Sloan Digital Sky Survey (SDSS), obtained 4 years prior to our observations. Comparing the redshifts of [OIII]4960,5008 we put a 3 sigma limit on the relative acceleration to be less than 32 km s^-1 yr^-1 between different emitting regions. Using the 2D spectra obtained at different position angles we show that the [OIII]5008 line from the zem = 0.7128 component is extended beyond the spectral point spread function. We infer the linear extent of this line emitting region is ~8 kpc. We also find a tentative evidence for an offset between the centroid of the [OIII]5008 line at zem = 0.7128 and the QSO trace when the slit is aligned at a position angle of 299 degrees. This corresponds to the zem = 0.7128 system being at an impact parameter of ~1 kpc with respect to the zem = 0.6972 in the north west direction. Based on our observations we conclude that the binary black hole model is most unlikely. The spatial extent and the sizes are consistent with both black hole recoil and merging scenarios.