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

Using spatially resolved spectroscopy from SDSS-IV MaNGA we have demonstrated that low ionisation emission line regions (LIERs) in local galaxies result from photoionisation by hot evolved stars, not active galactic nuclei. LIERs are ubiquitous in both quiescent galaxies and in the central regions of galaxies where star formation takes place at larger radii. We refer to these two classes of galaxies as extended LIER (eLIER) and central LIER (cLIER) galaxies respectively. cLIERs are late type galaxies located around the green valley, in the transition region between the star formation main sequence and quiescent galaxies. These galaxies display regular disc rotation in both stars and gas, although featuring a higher central stellar velocity dispersion than star forming galaxies of the same mass. cLIERs are consistent with being slowly quenched inside-out; the transformation is associated with massive bulges, pointing towards the importance of bulge growth via secular evolution. eLIERs are morphologically early types and are indistinguishable from passive galaxies devoid of line emission in terms of their stellar populations, morphology and central stellar velocity dispersion. Ionised gas in eLIERs shows both disturbed and disc-like kinematics. When a large-scale flow/rotation is observed in the gas, it is often misaligned relative to the stellar component. These features indicate that eLIERs are passive galaxies harbouring a residual cold gas component, acquired mostly via external accretion. Importantly, quiescent galaxies devoid of line emission reside in denser environments and have significantly higher satellite fraction than eLIERs. Environmental effects thus represent the likely cause for the existence of line-less galaxies on the red sequence.

We present a study on the stellar age and metallicity distributions for 1105 galaxies using the STARLIGHT software on MaNGA integral field spectra. We derive age and metallicity gradients by fitting straight lines to the radial profiles, and explore their correlations with total stellar mass M*, NUV-r colour and environments, as identified by both the large scale structure (LSS) type and the local density. We find that the mean age and metallicity gradients are close to zero but slightly negative, which is consistent with the inside-out formation scenario. Within our sample, we find that both the age and metallicity gradients show weak or no correlation with either the LSS type or local density environment. In addition, we also study the environmental dependence of age and metallicity values at the effective radii. The age and metallicity values are highly correlated with M* and NUV-r and are also dependent on LSS type as well as local density. Low-mass galaxies tend to be younger and have lower metallicity in low-density environments while high-mass galaxies are less affected by environment.

The nature of warm, ionized gas outside of galaxies may illuminate several key galaxy evolutionary processes. A serendipitous observation by the MaNGA survey has revealed a large, asymmetric H$\alpha$ complex with no optical counterpart that extends $\approx8"$ ($\approx6.3$ kpc) beyond the effective radius of a dusty, starbursting galaxy. This H$\alpha$ extension is approximately three times the effective radius of the host galaxy and displays a tail-like morphology. We analyze its gas-phase metallicities, gaseous kinematics, and emission-line ratios, and discuss whether this H$\alpha$ extension could be diffuse ionized gas, a gas accretion event, or something else. We find that this warm, ionized gas structure is most consistent with gas accretion through recycled wind material, which could be an important process that regulates the low-mass end of the galaxy stellar mass function.

We present results on the clustering of 282,068 galaxies in the Baryon Oscillation Spectroscopic Survey (BOSS) sample of massive galaxies with redshifts 0.4<z<0.7 which is part of the Sloan Digital Sky Survey III project. Our results cover a large range of scales from ~0.5 to ~90 Mpc/h. We compare these estimates with the expectations of the flat LCDM cosmological model with parameters compatible with WMAP7 data. We use the MultiDark cosmological simulation together with a simple halo abundance matching technique, to estimate galaxy correlation functions, power spectra, abundance of subhaloes and galaxy biases. We find that the LCDM model gives a reasonable description to the observed correlation functions at z~0.5, which is a remarkably good agreement considering that the model, once matched to the observed abundance of BOSS galaxies, does not have any free parameters. However, we find a deviation (>~10%) in the correlation functions for scales less than ~1 Mpc/h and ~10-40 Mpc/h. A more realistic abundance matching model and better statistics from upcoming observations are needed to clarify the situation. We also estimate that about 12% of the "galaxies" in the abundance-matched sample are satellites inhabiting central haloes with mass M>~1e14 M_sun/h. Using the MultiDark simulation we also study the real space halo bias b(r) of the matched catalogue finding that b=2.00+/-0.07 at large scales, consistent with the one obtained using the measured BOSS projected correlation function. Furthermore, the linear large-scale bias depends on the number density n of the abundance-matched sample as b=-0.048-(0.594+/-0.02)*log(n/(h/Mpc)^3). Extrapolating these results to BAO scales we measure a scale-dependent damping of the acoustic signal produced by non-linear evolution that leads to ~2-4% dips at ~3 sigma level for wavenumbers k>~0.1 h/Mpc in the linear large-scale bias.

We report the discovery of a low-mass companion orbiting the metal-rich, main sequence F star TYC 2949-00557-1 during the MARVELS (Multi-object APO Radial Velocity Exoplanet Large-area Survey) Pilot Project. The host star has an effective temperature T_eff = 6135 +/- 40 K, log(g) = 4.4 +/- 0.1 and [Fe/H] = 0.32 +/- 0.01, indicating a mass of M = 1.25 +/- 0.09 M_\odot and R = 1.15 +/- 0.15 R_\odot. The companion has an orbital period of 5.69449 +/- 0.00023 days and straddles the hydrogen burning limit with a minimum mass of 64 M_J, and may thus be an example of the rare class of brown dwarfs orbiting at distances comparable to those of "Hot Jupiters." We present relative photometry that demonstrates the host star is photometrically stable at the few millimagnitude level on time scales of hours to years, and rules out transits for a companion of radius greater than 0.8 R_J at the 95% confidence level. Tidal analysis of the system suggests that the star and companion are likely in a double synchronous state where both rotational and orbital synchronization have been achieved. This is the first low-mass companion detected with a multi-object, dispersed, fixed-delay interferometer.

We present evidence for a substantial overdensity of stars in the direction of the constellations of Hercules and Aquila. The Cloud is centered at a Galactic longitude of about 40 degrees and extends above and below the Galactic plane by at least 50 degrees. Given its off-centeredness and height, it is unlikely that the Hercules-Aquila Cloud is related to the bulge or thick disk. More likely, this is a new structural component of the Galaxy that passes through the disk. The Cloud stretches about 80 degrees in longitude. Its heliocentric distance lies between 10 and 20 kpc so that the extent of the Cloud in projection is roughly 20 kpc by 15 kpc. It has an absolute magnitude of -13 and its stellar population appears to be comparable to, but somewhat more metal-rich than, M92.