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 maps of the ionized gas flux distributions, excitation, star-formation rate SFR, surface mass density $\Sigma_{H+}$, and obtain total values of SFR and ionized gas masses {\it M} for 62 Active Galactic Nuclei (AGN) observed with SDSS-IV MaNGA and compare them with those of a control sample of 112 non-active galaxies. The most luminous AGN -- with $L(\rm{[OIII]}\lambda 5007) \ge 3.8\times 10^{40}\,\mbox{erg}\,\mbox{s}^{-1}$, and those hosted by earlier-type galaxies are dominated by Seyfert excitation within 0.2 effective radius $R_e$ from the nucleus, surrounded by LINER excitation or transition regions, while the less luminous and hosted by later-type galaxies show equally frequent LINER and Seyfert excitation within $0.2\,R_e$. The extent $R$ of the region ionized by the AGN follows the relation $R\propto\,L(\rm{[OIII]})^{0.5}$ -- as in the case of the Broad-Line Region. The SFR distribution over the region ionized by hot stars is similar for AGN and controls, while the integrated SFR -- in the range $10^{-3}-10$\,M$_\odot$\,yr$^{-1}$ is also similar for the late-type sub-sample, but higher in the AGN for 75\% of the early-type sub-sample. We thus conclude that there is no signature of AGN quenching star formation in the body of the galaxy in our sample. We also find that 66\% of the AGN have higher ionized gas masses $M$ than the controls -- in the range 10$^5-3\times10^7$\,M$_\odot$ -- while 75\% of the AGN have higher $\Sigma_{H+}$ within $0.2\,R_e$ than the control galaxies.

We investigate the effects of Active Galactic Nuclei (AGN) on the gas kinematics of their host galaxies, using MaNGA data for a sample of 62 AGN hosts and 109 control galaxies (inactive galaxies). We compare orientation of the line of nodes (kinematic Position Angle - PA) measured from the gas and stellar velocity fields for the two samples. We found that AGN hosts and control galaxies display similar kinematic PA offsets between gas and stars. However, we note that AGN have larger fractional velocity dispersion $\sigma$ differences between gas and stars [$\sigma_{frac}=(\sigma_{\rm gas}-\sigma_{stars})/\sigma_{\rm stars}$] when compared to their controls, as obtained from the velocity dispersion values of the central (nuclear) pixel (2.5" diameter). The AGN have a median value of $\sigma_{\rm frac}$ of $<\sigma_{frac}>_{\rm AGN}=0.04$, while the the median value for the control galaxies is $<\sigma_{frac}>_{\rm CTR}=-0.23$. 75% of the AGN show $\sigma_{frac}>-0.13$, while 75% of the normal galaxies show $\sigma_{\rm frac}<-0.04$, thus we suggest that the parameter $\sigma_{\rm frac}$ can be used as an indicative of AGN activity. We find a correlation between the [OIII]$\lambda$5007 luminosity and $\sigma_{frac}$ for our sample. Our main conclusion is that the AGN already observed with MaNGA are not powerful enough to produce important outflows at galactic scales, but at 1-2 kpc scales, AGN feedback signatures are always present on their host galaxies.

Super-massive black holes (SMBH) are present at the center of most galaxies, with the related mass accretion processes giving origin to outflows in Active Galactic Nuclei (AGN). It has been presumed that only intense winds from luminous AGN were able to suppress star formation until the discovery of a new class of galaxies with no recent star formation and with the nucleus in a quiescent state showing kpc scale outflows. We used SDSS MaNGA and Gemini Integral Field Spectroscopy of the prototype Red Geyser Akira and found that the orientation of the outflow changes by about 50$^\circ$ from its nucleus to kpc scales. A possible interpretation is that the outflow is produced by a precessing accretion disk due to a misalignment between the orientation of the disk and the spin of the SMBH. The precession of the central source is also supported by a similar change in the orientation of the ionization pattern. Although similar behavior has commonly being reported for collimated relativistic jets, the precession of an AGN wide wind is reported here for the first time, implying on a larger work surface of the wind, which in turn increases the star formation suppression efficiency of the outflow.

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.

We present spatially resolved stellar population age maps, average radial profiles and gradients for the first 62 Active Galactic Nuclei (AGN) observed with SDSS-IV MaNGA to study the effects of the active nuclei on the star formation history of the host galaxies. These results, derived using the STARLIGHT code, are compared with a control sample of non-active galaxies matching the properties of the AGN hosts. We find that the fraction of young stellar populations (SP) in high-luminosity AGN is higher in the inner ($R \leq 0.5\,R_e$) regions when compared with the control sample; low-luminosity AGN, on the other hand, present very similar fractions of young stars to the control sample hosts for the entire studied range ($1\,R_e$). The fraction of intermediate age SP of the AGN hosts increases outwards, with a clear enhancement when compared with the control sample. The inner region of the galaxies (AGN and control galaxies) presents a dominant old SP, whose fraction decreases outwards. We also compare our results (differences between AGN and control galaxies) for the early and late-type hosts and find no significant differences. In summary, our results suggest that the most luminous AGN seems to have been triggered by a recent supply of gas that has also triggered recent star formation ($t\,\leq\,40\,Myrs$) in the central region.

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 new evidence for AGN feedback in a subset of 69 quenched low-mass galaxies ($M_{\star} \lesssim 5\times10^{9}$ M$_{\odot}$, $M_{\rm{r}} > -19$) selected from the first two years of the SDSS-IV MaNGA survey. The majority (85 per cent) of these quenched galaxies appear to reside in a group environment. We find 6 galaxies in our sample that appear to have an active AGN that is preventing on-going star-formation; this is the first time such a feedback mechanism has been observed in this mass range. Interestingly, five of these six galaxies have an ionised gas component that is kinematically offset from their stellar component, suggesting the gas is either recently accreted or outflowing. We hypothesise these six galaxies are low-mass equivalents to the "red geysers" observed in more massive galaxies. Of the other 63 galaxies in the sample, we find 8 do appear for have some low-level, residual star formation, or emission from hot, evolved stars. The remaining galaxies in our sample have no detectable ionised gas emission throughout their structures, consistent with them being quenched. This work shows the potential for understanding the detailed physical properties of dwarf galaxies through spatially resolved spectroscopy.

We report the characterization of the first $62$ MaNGA Active Galactic Nuclei (AGN) hosts in the Fifth Product Launch (MPL-5) and the definition of a control sample of non-active galaxies. This control sample - comprising two galaxies for each AGN - was selected in order to match the AGN hosts in terms of stellar mass, redshift, visual morphology and inclination. The stellar masses are in the range $9.4<\log(M/M_\odot)<11.5$, and most objects have redshifts $\leq 0.08$. The AGN sample is mostly comprised of low-luminosity AGN, with only 17 nuclei with $L([OIII]\lambda 5007)\geq 3.8\times 10^{40}$ erg s$^{-1}$ (that we call "strong AGN"). The stellar population of the control sample galaxies within the inner $1$-$3$ kpc is dominated by the old ($\sim$ $4$ - $13$ Gyr) age component, with a small contribution of intermediate age ($\sim 640$-$940$ Myr) and young stars ($\leq 40$ Myr) to the total light at $5700\AA$. While the weaker AGN show a similar age distribution to that of the control galaxies, the strong AGN show an increased contribution of younger stars and a decreased contribution of older stars. Examining the relationship between the AGN stellar population properties and $L([OIII])$, we find that with increasing $L([OIII])$, the AGN exhibit a decreasing contribution from the oldest ($>4$ Gyr) stellar population relative to control galaxies, but have an increasing contribution from the younger components with ages $\sim 40$ Myr. We also find a correlation of the mean age differences (AGN - control) with $L([OIII])$, in the sense that more luminous AGN are younger than the control objects, while the low-luminosity AGN seem to be older. These results support a connection between the growth of the galaxy bulge via formation of new stars and the growth of the Supermassive Black Hole via matter accretion in the AGN phase.

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.

The understating of the origin of the H2 line emission from the central regions of galaxies represent an important key to improve our knowledge about the excitation and ionization conditions of the gas in these locations. Usually these lines can be produced by Starburts, shocks and/or radiation from an active galactic nucleus (AGN). Luminous Infrared Galaxies (LIRG) represent ideal and challenging objects to investigate the origin of the H2 emission, as all processes above can be observed in a single object. In this work, we use K-band integral field spectroscopy to map the emission line flux distributions and kinematics and investigate the origin of the molecular and ionized gas line emission from inner 1.4x2.4 kpc2 of the LIRG NGC6240, known to be the galaxy with strongest H2 line emission. The emission lines show complex profiles at locations between both nuclei and surrounding the northern nucleus, while at locations near the southern nucleus and at 1" west of the northern nucleus, they can be reproduced by a single gaussian component. We found that the H2 emission is originated mainly by thermal processes, possible being dominated by heating of the gas by X-rays from the AGN at locations near both nuclei. For the region between the northern and southern nuclei shocks due to the interacting process may be the main excitation mechanism, as indicated by the high values of the H2l2.12um/Brg line ratio. A contribution of fluorescent excitation may also be important at locations near 1" west of the northern nucleus, which show the lowest line ratios. The [FeII]l2.072um/Brg ratio show a similar trend as observed for H2l2.12um/Brg, suggesting that [FeII] and H2 line emission have similar origins. Finally, the [CaVIII]l2.32um coronal line emission is observed mainly in regions next to the nuclei, suggesting it is originated gas ionized by the radiation from the AGN.