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

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 study radial profiles in H$\alpha$ equivalent width and specific star formation rate (sSFR) derived from spatially-resolved SDSS-IV MaNGA spectroscopy to gain insight on the physical mechanisms that suppress star formation and determine a galaxy's location in the SFR-$\rm M_\star$ diagram. Even within the star-forming `main sequence', the measured sSFR decreases with stellar mass, both in an integrated and spatially-resolved sense. Flat sSFR radial profiles are observed for $\rm log(M_\star/ M_\odot) < 10.5$, while star-forming galaxies of higher mass show a significant decrease in sSFR in the central regions, a likely consequence of both larger bulges and an inside-out growth history. Our primary focus is the green valley, constituted by galaxies lying below the star formation main sequence, but not fully passive. In the green valley we find sSFR profiles that are suppressed with respect to star-forming galaxies of the same mass at all galactocentric distances out to 2 effective radii. The responsible quenching mechanism therefore appears to affect the entire galaxy, not simply an expanding central region. The majority of green valley galaxies of $\rm log(M_\star/ M_\odot) > 10.0$ are classified spectroscopically as central low-ionisation emission-line regions (cLIERs). Despite displaying a higher central stellar mass concentration, the sSFR suppression observed in cLIER galaxies is not simply due to the larger mass of the bulge. Drawing a comparison sample of star forming galaxies with the same $\rm M_\star$ and $\rm \Sigma_{1~kpc}$ (the mass surface density within 1 kpc), we show that a high $\rm \Sigma_{1~kpc}$ is not a sufficient condition for determining central quiescence.

We present here the characterization of the main properties of a sample of 98 AGN host galaxies, both type-II and type-I, in comparison with those of about 2700 non-active galaxies observed by the MaNGA survey. We found that AGN hosts are morphologically early-type or early-spirals. For a given morphology AGN hosts are, in average, more massive, more compact, more central peaked and rather pressurethan rotational-supported systems. We confirm previous results indicating that AGN hosts are located in the intermediate/transition region between star-forming and non-star-forming galaxies (i.e., the so-called green valley), both in the ColorMagnitude and the star formation main sequence diagrams. Taking into account their relative distribution in terms of the stellar metallicity and oxygen gas abundance and a rough estimation of their molecular gas content, we consider that these galaxies are in the process of halting/quenching the star formation, in an actual transition between both groups. The analysis of the radial distributions of the starformation rate, specific star-formation rate, and molecular gas density shows that the quenching happens from inside-out involving both a decrease of the efficiency of the star formation and a deficit of molecular gas. All the intermediate data-products used to derive the results of our analysis are distributed in a database including the spatial distribution and average properties of the stellar populations and ionized gas, published as a Sloan Digital Sky Survey Value Added Catalog being part of the 14th Data Release: http://www.sdss.org/dr14/manga/manga-data/manga-pipe3d-value-added-catalog/

Galaxy integrated H{\alpha} star formation rate-stellar mass relation, or SFR(global)-M*(global) relation, is crucial for understanding star formation history and evolution of galaxies. However, many studies have dealt with SFR using unresolved measurements, which makes it difficult to separate out the contamination from other ionizing sources, such as active galactic nuclei and evolved stars. Using the integral field spectroscopic observations from SDSS-IV MaNGA, we spatially disentangle the contribution from different H{\alpha} powering sources for ~1000 galaxies. We find that, when including regions dominated by all ionizing sources in galaxies, the spatially-resolved relation between H{\alpha} surface density ({\Sigma}H{\alpha}(all)) and stellar mass surface density ({\Sigma}*(all)) progressively turns over at high {\Sigma}*(all) end for increasing M*(global) and bulge dominance (bulge-to-total light ratio, B/T). This in turn leads to the flattening of the integrated H{\alpha}(global)-M*(global) relation in the literature. By contrast, there is no noticeable flattening in both integrated H{\alpha}(HII)-M*(HII) and spatially-resolved {\Sigma}H{\alpha}(HII)-{\Sigma}*(HII) relations when only regions where star formation dominates the ionization are considered. In other words, the flattening can be attributed to the increasing regions powered by non-star-formation sources, which generally have lower ionizing ability than star formation. Analysis of the fractional contribution of non-star-formation sources to total H{\alpha} luminosity of a galaxy suggests a decreasing role of star formation as an ionizing source toward high-mass, high-B/T galaxies and bulge regions. This result indicates that the appearance of the galaxy integrated SFR-M* relation critically depends on their global properties (M*(global) and B/T) and relative abundances of various ionizing sources within the galaxies.

We present an updated version of the mass--metallicity relation (MZR) using integral field spectroscopy data obtained from 734 galaxies observed by the CALIFA survey. These unparalleled spatially resolved spectroscopic data allow us to determine the metallicity at the same physical scale ($\mathrm{R_{e}}$) for different calibrators. We obtain MZ relations with similar shapes for all calibrators, once the scale factors among them are taken into account. We do not find any significant secondary relation of the MZR with either the star formation rate (SFR) or the specific SFR for any of the calibrators used in this study, based on the analysis of the residuals of the best fitted relation. However we do see a hint for a (s)SFR-dependent deviation of the MZ-relation at low masses (M$<$10$^{9.5}$M$_\odot$), where our sample is not complete. We are thus unable to confirm the results by Mannucci et al. (2010), although we cannot exclude that this result is due to the differences in the analysed datasets. In contrast, our results are inconsistent with the results by Lara-Lopez et al. (2010), and we can exclude the presence of a SFR-Mass-Oxygen abundance Fundamental Plane. These results agree with previous findings suggesting that either (1) the secondary relation with the SFR could be induced by an aperture effect in single fiber/aperture spectroscopic surveys, (2) it could be related to a local effect confined to the central regions of galaxies, or (3) it is just restricted to the low-mass regime, or a combination of the three effects.

We present near-infrared spectra, obtained with SINFONI and XShooter observations at ESO VLT, of nine dusty, red QSOs at 1.2<z<2.6. The sources are hard X-ray detected, characterized by cold absorption (N_H>10^{21} - 10^{22} cm^{-2}) and show a broad Ha component in the NIR spectra. We complement this sample with twelve additional sources taken from the literature with similar properties resulting in a total sample of 21 X-ray obscured, intermediate type (1.8-1.9), dusty reddened QSOs. From the broad Ha line we have computed the BH masses through the virial formula and derived Eddington ratios. Moreover, from optical/IR multi-component SED fitting we have derived the stellar mass of their host galaxies and their SFRs. We find that most of the sources in our sample are hosted in starburst and main sequence star-forming galaxies with Eddington ratios lambda>0.1. We find a strong trend with the BH mass i.e. less massive objects are scattered below and above the local relation while the most massive ones are mainly located above it. We also studied the evolution of these sources on the M_BH-M_star plane compared to a sample of optically blue type--1 QSOs and we find that obscured red QSOs show a ratio of M_BH to M_star that increases with redshift which is consistent with or slightly lower than what has been found for blue QSOs. These sources may represent the blow-out phase at the end of the rapid BH growth and immediately preceding the classical blue QSOs typically sampled in optical surveys. They in fact show evidence of outflows in the ionized gas component, but their BH has already fully formed.

Most galaxy evolutionary models require quasar feedback to regulate star formation in their host galaxies. In particular, at high redshift, models expect that feedback associated with quasar-driven outflows is so efficient that the gas in the host galaxy is largely swept away or heated up, hence suppressing star formation in massive galaxies. We observationally investigate this phenomenon by using VLT-SINFONI integral field spectroscopy of the luminous quasar 2QZJ002830.4-281706 at z=2.4. The spectra sample the optical emission lines redshifted into the near-IR. The [OIII]5007 emission-line kinematics map reveals a massive outflow on scales of several kpc. The detection of narrow Halpha emission reveals star formation in the quasar host galaxy, with SFR=100 Msun/yr. However, the star formation is not distributed uniformly, but is strongly suppressed in the region with the highest outflow velocity and highest velocity dispersion. This result indicates that star formation in this region is strongly quenched by the quasar outflow, which is cleaning the galaxy disk of its molecular gas. This is one of the first direct observational proofs of quasar feedback quenching the star formation at high redshift.