Twenty years have passed since first light for the Sloan Digital Sky Survey (SDSS). Here, we release data taken by the fourth phase of SDSS (SDSS-IV) across its first three years of operation (July 2014-July 2017). This is the third data release for SDSS-IV, and the fifteenth from SDSS (Data Release Fifteen; DR15). New data come from MaNGA - we release 4824 datacubes, as well as the first stellar spectra in the MaNGA Stellar Library (MaStar), the first set of survey-supported analysis products (e.g. stellar and gas kinematics, emission line, and other maps) from the MaNGA Data Analysis Pipeline (DAP), and a new data visualisation and access tool we call "Marvin". The next data release, DR16, will include new data from both APOGEE-2 and eBOSS; those surveys release no new data here, but we document updates and corrections to their data processing pipelines. The release is cumulative; it also includes the most recent reductions and calibrations of all data taken by SDSS since first light. In this paper we describe the location and format of the data and tools and cite technical references describing how it was obtained and processed. The SDSS website (www.sdss.org) has also been updated, providing links to data downloads, tutorials and examples of data use. While SDSS-IV will continue to collect astronomical data until 2020, and will be followed by SDSS-V (2020-2025), we end this paper by describing plans to ensure the sustainability of the SDSS data archive for many years beyond the collection of data.

The design and performance of the LUX-ZEPLIN (LZ) detector is described as of March 2015 in this Conceptual Design Report. LZ is a second-generation dark-matter detector with the potential for unprecedented sensitivity to weakly interacting massive particles (WIMPs) of masses from a few GeV/c2 to hundreds of TeV/c2. With total liquid xenon mass of about 10 tonnes, LZ will be the most sensitive experiment for WIMPs in this mass region by the end of the decade. This report describes in detail the design of the LZ technical systems. Expected backgrounds are quantified and the performance of the experiment is presented. The LZ detector will be located at the Sanford Underground Research Facility in South Dakota. The organization of the LZ Project and a summary of the expected cost and current schedule are given.

We investigate the hydrodynamic evolution of Gamma--Ray Burst (GRB) remnants at late times (>> 1 year), by extending the standard fireball model into the nonrelativistic Sedov--Taylor regime. We calculate the associated synchrotron luminosity as a function of remnant age. The radio luminosity of old remnants depends strongly on the power-law index p= (dlogN/dlogE) of the kinetic energy distribution of shock-accelerated electrons. In addition, the detection probability of remnants in nearby galaxies depends on the GRB rate in the local universe. We predict that if p<2.6 and the BATSE-calibrated GRB rate per comoving volume is independent of redshift (despite the strong evolution in the cosmic star formation rate) then the VLA can discover a 10^3 years-old GRB remnant in the Virgo cluster. Nearby GRB remnants could be searched for and resolved at radio frequencies, and then identified as distinct from supernova remnants based on the detection of the extended emission of high-ionization lines in follow--up optical observations.

We calculate the prompt x-ray emission as a function of viewing angle for beamed Gamma-Ray Burst (GRB) sources. Prompt x-rays are inevitable due to the less highly blueshifted photons emitted at angles greater than 1/gamma relative to the beam symmetry axis, where gamma is the expansion Lorentz factor. The observed flux depends on the combinations (gamma Delta theta) and (gamma theta_v), where (Delta theta) is the beaming angle and theta_v is the viewing angle. We use the observed source counts of gamma-ray-selected GRBs to predict the minimum detection rate of prompt x-ray bursts as a function of limiting sensitivity. We compare our predictions with the results from the Ariel V catalog of fast x-ray transients, and find that Ariel's sensitivity is not great enough to place significant constraints on gamma and (Delta theta). We estimate that a detector with fluence limit ~10^{-7} erg/cm^2 in the 2-10 keV channel will be necessary to distinguish between geometries. Because the x-ray emission is simultaneous with the GRB emission, our predicted constraints do not involve any model assumptions about the emission physics but simply follow from special-relativistic considerations.

We predict the rate at which Gamma-Ray Burst (GRB) afterglows should be detected in supernova searches as a function of limiting flux. Although GRB afterglows are rarer than supernovae, they are detectable at greater distances because of their higher intrinsic luminosity. Assuming that GRBs trace the cosmic star formation history and that every GRB gives rise to a bright afterglow, we find that the average detection rate of supernovae and afterglows should be comparable at limiting magnitudes brighter than K=18. The actual rate of afterglows is expected to be somewhat lower since only a fraction of all gamma-ray selected GRBs were observed to have associated afterglows. However, the rate could also be higher if the initial gamma-ray emission from GRB sources is more beamed than their late afterglow emission. Hence, current and future supernova searches can place strong constraints on the afterglow appearance fraction and the initial beaming angle of GRB sources.

We use the existing catalog of Damped Lyman-Alpha (DLA) systems to place constraints on the amount of evolution in the baryonic content of galaxies and on the value of the cosmological constant. The density of cold gas at redshifts z=3+-1 is obtained from the mean HI column density of DLAs per cosmological path length. This path length per unit redshift is in turn a sensitive function of the vacuum density parameter, Omega_v. We compare the total inferred mass of cold gas at high redshifts to that observed in stars today for flat cosmologies. We define "eta" to be net fraction of the baryonic content of local galaxies which was expelled since z=3, and use Bayesian inference to derive confidence regions in the (eta, Omega_v) plane. In all cosmologies we find that eta<0.4 with at least 95% confidence if <25% of the current stellar population formed before z=3. The most likely value of eta is negative, implying a net increase by several tens of percent in the baryonic mass of galaxies since z=3+-1. On the other hand, recent observations of high metal abundances in the intracluster medium of rich clusters (Loewenstein & Mushotzky 1996) require that metal-rich gas be expelled from galaxies in an amount approximately equal to the current mass in stars. Based on our results and the low metallicity observed in DLAs at z>2, we infer that more than half of the baryonic mass processed through galaxies must have been assembled and partly expelled from galaxies after z=2. We expect our constraints to improve considerably as the size of the DLA sample will increase with the forthcoming Sloan Digital Sky Survey.

We discuss several constraints on the properties of gamma-ray bursts (GRB) at cosmological distances. First we use the requirement that burst sources must be optically thin to pair production by photon-photon collisions in order to produce the observed nonthermal spectra. We derive probability distributions for the minimum Lorentz expansion factor, the radiation energy, the maximum baryonic mass, and the maximum surrounding gas density in the events, based on 254 events from the second BATSE catalog. Secondly, the time variabilities of the bursts in the BATSE catalog are used to place an upper bound on the curvature radius of the emitting surface in the events. This is based on the requirement that the emitting region seen by the observer must be sufficiently small to produce the observed variability without violating causality. Finally, we discuss the association of cosmological GRB with galaxies. We consider eight bright and well-localized bursts detected by the Pioneer Venus Orbiter, whose positional error boxes contain no bright galaxies. Using the local luminosity function of galaxies, we calculate the probability for not seeing the GRB host galaxy. This probability tends to increase as the width of the GRB luminosity function increases. However, the allowed width of the GRB luminosity function is restricted by the number-count statistics of the bursts.

We place empirical constraints on the physical properties of $\gamma-$ray burst events at cosmological distances. In particular we derive probability distributions for the radiation energy $E_\gamma$, the minimum Lorentz factor $\gamma_{\rm min}$, the maximum baryonic mass $M_{\rm max}$ and the upper bound on the surrounding gas density $n_{\rm max}$ in the events, based on 169 bursts from the first BATSE catalog. Using peak flux as a distance indicator we probe bursts where the constraints are stronger than average. The resulting variance and skewness of the cosmological probability distributions are calculated in addition to their mean values: $\langle E_\gamma\rangle=4\times 10^{51} h^{-2} {\rm erg}$, $\langle \gamma_{\rm min}\rangle= 5\times 10^2$, $\langle M_{\rm max}\rangle=10^{-5}\xi^{-1} M_\odot$, and $\langle (\gamma/\gamma_{\rm min})^{-5} n_{\rm max}\rangle\lsim 10^{2} \xi^{-1}{\rm cm^{-3}}$, where $\xi$ is the fraction of the total energy which is converted to $\gamma-$rays. The distribution of burst energies ends at about $10^{53} {\rm erg}$, close to the binding energy of a neutron star.