Mapping the local and distant Universe is key to our understanding of it. For decades, the Sloan Digital Sky Survey (SDSS) has made a concerted effort to map millions of celestial objects to constrain the physical processes that govern our Universe. The most recent and fifth generation of SDSS (SDSS-V) is organized into three scientific ``mappers". Milky Way Mapper (MWM) that aims to chart the various components of the Milky Way and constrain its formation and assembly, Black Hole Mapper (BHM), which focuses on understanding supermassive black holes in distant galaxies across the Universe, and Local Volume Mapper (LVM), which uses integral field spectroscopy to map the ionized interstellar medium in the local group. This paper describes and outlines the scope and content for the nineteenth data release (DR19) of SDSS and the most substantial to date in SDSS-V. DR19 is the first to contain data from all three mappers. Additionally, we also describe nine value added catalogs (VACs) that enhance the science that can be conducted with the SDSS-V data. Finally, we discuss how to access SDSS DR19 and provide illustrative examples and tutorials.

We report the discovery of NGTS-11 c, a transiting warm Neptune ($P \approx 12.8$ d; $M_{p} = 1.2^{+0.3}_{-0.2} M_{\mathrm{Nep}}$; $R_{p} = 1.24 \pm 0.03 R_{\mathrm{Nep}}$), in an orbit interior to the previously reported transiting warm Saturn NGTS-11 b ($P \approx 35.5$ d). We also find evidence of a third outer companion orbiting the K-dwarf NGTS-11. We first detected transits of NGTS-11 c in TESS light curves and confirmed them with follow-up transits from NGTS and many other ground-based facilities. Radial-velocity monitoring with the HARPS and FEROS spectrographs revealed the mass of NGTS-11 c and provides evidence for a long-period companion ($P > 2300$ d; $M_{p} \sin i > 3.6 M_{\mathrm{Jup}}$). Taking into account the two additional bodies in our expanded datasets, we find that the mass of NGTS-11 b ($M_{p} = 0.63 \pm 0.09 M_{\mathrm{Sat}}$; $R_{p} = 0.97 \pm 0.02 R_{\mathrm{Sat}}$) is lower than previously reported ($M_{p} = 1.2 \pm 0.3 M_{\mathrm{Sat}}$). Given their near-circular and compact orbits, NGTS-11 c and b are unlikely to have reached their present locations via high-eccentricity migration. Instead, they probably either formed in situ or formed farther out and then underwent disk migration. A comparison of NGTS-11 with the eight other known systems hosting multiple well-characterized warm giants shows that it is most similar to Kepler-56. Finally, we find that the commonly used 10-day boundary between hot and warm Jupiters is empirically well supported.

As part of the KESPRINT collaboration, we present the discovery and characterization of three exoplanets in the sub-Neptune to super-Neptune regime, spanning key regions of the exo-Neptunian landscape. TOI-1472c and TOI-1648b are newly discovered sub-Neptunes, while TOI-1472b is a previously known super-Neptune for which we provide an improved mass measurement. These planets have orbital periods of 6--15 days and radii of 2.5--4.1 R$_\oplus$, probing regions where planet formation and atmospheric evolution remain poorly understood. We combine TESS transit photometry with ground-based radial velocities to determine precise masses, radii, and orbital properties. TOI-1472b has a mass of $18.0^{+0.84}_{-0.85}$ M$_\oplus$ and a radius of $4.06 \pm 0.10$ R$_\oplus$, TOI-1472c has a mass of $21.1^{+0.96}_{-0.99}$ M$_\oplus$ and a radius of $3.33 \pm 0.08$ R$_\oplus$, and TOI-1648b has a mass of $7.4^{+1.1}_{-1.3}$ M$_\oplus$ and a radius of $2.54^{+0.14}_{-0.12}$ R$_\oplus$. The planets exhibit a range of eccentricities (0.041--0.178), indicating diverse evolutionary histories. TOI-1648b, with a high Transmission Spectroscopy Metric (TSM $\sim$59), is a promising target for atmospheric characterization. Together, these three planets provide precise constraints on the structure, composition, and dynamical evolution of small to intermediate-sized exoplanets, enriching our understanding of the exo-Neptunian landscape.

In the last few decades planet search surveys have been focusing on solar type stars, and only recently the high-mass regimes. This is mostly due to challenges arising from the lack of instrumental precision, and more importantly, the inherent active nature of fast rotating massive stars. Here we report NGTS-33b (TOI-6442b), a super-Jupiter planet with mass, radius and orbital period of 3.6 $\pm$ 0.3 M$_{\rm jup}$, 1.64 $\pm$ 0.07 R$_{\rm jup}$ and $2.827972 \pm 0.000001$ days, respectively. The host is a fast rotating ($0.6654 \pm 0.0006$ day) and hot (T$_{\rm eff}$ = 7437 $\pm$ 72 K) A9V type star, with a mass and radius of 1.60 $\pm$ 0.11 M$_{\odot}$ and 1.47 $\pm$ 0.06 R$_{\odot}$, respectively. Planet structure and Gyrochronology models shows that NGTS-33 is also very young with age limits of 10-50 Myr. In addition, membership analysis points towards the star being part of the Vela OB2 association, which has an age of $\sim$ 20-35 Myr, thus providing further evidences about the young nature of NGTS-33. Its low bulk density of 0.19$\pm$0.03 g cm$^{-3}$ is 13$\%$ smaller than expected when compared to transiting hot Jupiters with similar masses. Such cannot be solely explained by its age, where an up to 15$\%$ inflated atmosphere is expected from planet structure models. Finally, we found that its emission spectroscopy metric is similar to JWST community targets, making the planet an interesting target for atmospheric follow-up. Therefore, NGTS-33b's discovery will not only add to the scarce population of young, massive and hot Jupiters, but will also help place further strong constraints on current formation and evolution models for such planetary systems.

We analyse 829,481 stars from the Next Generation Transit Survey (NGTS) to extract variability periods. We utilise a generalisation of the autocorrelation function (the G-ACF), which applies to irregularly sampled time series data. We extract variability periods for 16,880 stars from late-A through to mid-M spectral types and periods between 0.1 and 130 days with no assumed variability model. We find variable signals associated with a number of astrophysical phenomena, including stellar rotation, pulsations and multiple-star systems. The extracted variability periods are compared with stellar parameters taken from Gaia DR2, which allows us to identify distinct regions of variability in the Hertzsprung-Russell Diagram. We explore a sample of rotational main-sequence objects in period-colour space, in which we observe a dearth of rotation periods between 15 and 25 days. This 'bi-modality' was previously only seen in space-based data. We demonstrate that stars in sub-samples above and below the period gap appear to arise from a stellar population not significantly contaminated by excess multiple systems. We also observe a small population of long-period variable M-dwarfs, which highlight a departure from the predictions made by rotational evolution models fitted to solar-type main-sequence objects. The NGTS data spans a period and spectral type range that links previous rotation studies such as those using data from Kepler, K2 and MEarth.

We re-investigate the HD208487 system to test the reality of the proposed HD208487c world. We also search for additional companions using applied Bayesian statistics and 15+ years of new RV data from the HARPS and the PFS instruments that were taken post-discovery of HD208487b. The RV data was analyzed with GLS Periodograms, followed by Bayesian analysis using the EMPEROR code. We scrutinised various stellar activity indices to search for any corresponding peaks in the power spectra, correlations with the RV measurements, or significant signals from a Bayesian analysis methodology. Finally, photometric data was checked to test for any transits or possible activity manifestations that could lead to possible false RV signals or excess noise. Our analysis points towards a candidate second planet in the system, positioned near the period of a previously proposed and subsequently challenged signal. This signal, HD208487c, would relate to a cool Saturn with an orbital period of 923.06 +2.02 -2.76 d and a minimum mass of Mj sini = 0.32 +/- 0.01Mj. Our analysis also gives rise to a newly discovered candidate planet, HD208487d, which would be the result of a cool super-Neptune/sub-Saturn with a period of 1380.13 +19.20 -8.25 d and a minimum mass of Mj sini = 0.15 +/- 0.01Mj. Neither stellar activity indices nor photometric data show signals statistically matching these periods. We have uncovered a candidate three planet system that would consist of an inner gas giant, a central Saturn and an outer super-Neptune/sub-Saturn. A dynamical analysis suggests that gravitational scattering of an initially ordered, equally-spaced system in a long resonant chain of six Neptunes can explain the current proposed architecture of HD208487. More RVs may also shed light on the reality of a fourth Doppler signal uncovered in the data that sits close to the 2:1 period-ratio with signal of HD208487c.

We statistically validated a sample of hot Neptune candidates applying a two-step vetting technique using DAVE and TRICERATOPS. We performed a systematic validation of 250 transit-like events in the Transiting Exoplanet Survey Satellite (TESS) archive in the parameter region defined by $P\leq 4$ d and $3R_\oplus\leq R\leq 5R_\oplus$. Through our analysis, we identified 18 hot Neptune-sized candidates, with a false positive probability &lt;50\%. Nine of these planet candidates still need to be confirmed. For each of the nine targets we retrieved the stellar parameters using ARIADNE and derived constraints on the planetary parameters by fitting the lightcurves with the juliet package. Within this sample of nine candidates, we statistically validated (i.e, with false positive probability < $0.3\%$) two systems (TOI-277 b and TOI-1288 b) by re-processing the candidates with TRICERATOPS along with follow-up observations. These new validated exoplanets expand the known hot Neptunes population and are high-priority targets for future radial velocities follow-up.

We report the discovery of a new ultra-short period hot Jupiter from the Next Generation Transit Survey. NGTS-6b orbits its star with a period of 21.17~h, and has a mass and radius of $1.330^{+0.024}_{-0.028}$\mjup\, and $1.271^{+0.197}_{-0.188}$\rjup\, respectively, returning a planetary bulk density of 0.711$^{+0.214}_{-0.136}$~g~cm$^{-3}$. Conforming to the currently known small population of ultra-short period hot Jupiters, the planet appears to orbit a metal-rich star ([Fe/H]$=+0.11\pm0.09$~dex). Photoevaporation models suggest the planet should have lost 5\% of its gaseous atmosphere over the course of the 9.6~Gyrs of evolution of the system. NGTS-6b adds to the small, but growing list of ultra-short period gas giant planets, and will help us to understand the dominant formation and evolutionary mechanisms that govern this population.

About one out of 200 Sun-like stars has a planet with an orbital period shorter than one day: an ultra-short-period planet (Sanchis-ojeda et al. 2014; Winn et al. 2018). All of the previously known ultra-short-period planets are either hot Jupiters, with sizes above 10 Earth radii (Re), or apparently rocky planets smaller than 2 Re. Such lack of planets of intermediate size (the "hot Neptune desert") has been interpreted as the inability of low-mass planets to retain any hydrogen/helium (H/He) envelope in the face of strong stellar irradiation. Here, we report the discovery of an ultra-short-period planet with a radius of 4.6 Re and a mass of 29 Me, firmly in the hot Neptune desert. Data from the Transiting Exoplanet Survey Satellite (Ricker et al. 2015) revealed transits of the bright Sun-like star \starname\, every 0.79 days. The planet's mean density is similar to that of Neptune, and according to thermal evolution models, it has a H/He-rich envelope constituting 9.0^(+2.7)_(-2.9)% of the total mass. With an equilibrium temperature around 2000 K, it is unclear how this "ultra-hot Neptune" managed to retain such an envelope. Follow-up observations of the planet's atmosphere to better understand its origin and physical nature will be facilitated by the star's brightness (Vmag=9.8).

We investigate the evolution of the galaxy Star Formation Rate Function (SFRF) and Cosmic Star Formation Rate Density (CSFRD) of $z\sim 0-8$ galaxies in the Evolution and Assembly of GaLaxies and their Environments (EAGLE) simulations. In addition, we present a compilation of UV, IR and H$\alpha$ SFRFs and compare these with the predictions from the EAGLE suite of cosmological hydrodynamic simulations. We find that the constraints implied by different indicators are inconsistent with each other for the highest star-forming objects at z < 2, a problem that is possibly related to selection biases and the uncertainties of dust attenuation effects. EAGLE's feedback parameters were calibrated to reproduce realistic galaxy sizes and stellar masses at z = 0.1. In this work we test if and why those choices yield realistic Star Formation Rates (SFRs) for $z \sim 0-8$ as well. We demonstrate that SNe feedback plays a major role at setting the abundance of galaxies at all star-forming regimes, especially at high redshifts. On the contrary, Active Galactic Nuclei (AGN) feedback becomes more prominent at lower redshifts and is a major mechanism that affects only the highest star-forming systems. Furthermore, we find that galaxies with SFR $\sim 1-10 \, {\rm M_{\odot} \, yr^{-1}}$ dominate the CSFRD at redshifts z < 5, while rare high star-forming galaxies (SFR $\sim 10-100 \,{\rm M_{\odot} \, yr^{-1}}$) contribute significantly only briefly around the peak era ($z \sim 2$) and then are quenched by AGN feedback. In the absence of this prescription objects with SFR $\sim 10-100 \,{\rm M_{\odot} \, yr^{-1}}$ would dominate the CSFRD, while the cosmic budget of star formation would be extremely high. Finally, we demonstrate that the majority of the cosmic star formation occurs in relatively rare high mass halos (${\rm M_{Halo}} \sim 10^{11-13} \, {\rm M_{\odot}}$) even at the earliest epochs.

Observations of circumstellar environments to look for the direct signal of exoplanets and the scattered light from disks has significant instrumental implications. In the past 15 years, major developments in adaptive optics, coronagraphy, optical manufacturing, wavefront sensing and data processing, together with a consistent global system analysis have enabled a new generation of high-contrast imagers and spectrographs on large ground-based telescopes with much better performance. One of the most productive is the Spectro-Polarimetic High contrast imager for Exoplanets REsearch (SPHERE) designed and built for the ESO Very Large Telescope (VLT) in Chile. SPHERE includes an extreme adaptive optics system, a highly stable common path interface, several types of coronagraphs and three science instruments. Two of them, the Integral Field Spectrograph (IFS) and the Infra-Red Dual-band Imager and Spectrograph (IRDIS), are designed to efficiently cover the near-infrared (NIR) range in a single observation for efficient young planet search. The third one, ZIMPOL, is designed for visible (VIR) polarimetric observation to look for the reflected light of exoplanets and the light scattered by debris disks. This suite of three science instruments enables to study circumstellar environments at unprecedented angular resolution both in the visible and the near-infrared. In this work, we present the complete instrument and its on-sky performance after 4 years of operations at the VLT.

Giant stars have been shown to be rich hunting grounds for those aiming to detect giant planets orbiting beyond ~0.5 AU. Here we present two planetary systems around bright giant stars, found by combining the radial-velocity (RV) measurements from the EXPRESS and PPPS projects, and using a Bayesian framework. HIP18606 is a naked-eye (V=5.8 mags) K0III star and is found to host a planet with an orbital period of ~675 days, a minimum mass (Msini) of 0.8 MJ, and a circular orbit. HIP111909 is a bright (V=7.4 mags) K1III star, and hosts two giant planets on circular orbits with minimum masses of Msini=1.2 MJ and Msini=0.8 MJ, and orbital periods of ~490 d and ~890 d, for planets b and c respectively, strikingly close to the 5:3 orbital period ratio. Analysis of 11 known giant star planetary systems arrive at broadly similar parameters to those published, whilst adding a further two new worlds orbiting these stars. With these new discoveries, we have found a total of 13 planetary systems (including three multiple systems) within the 37 giant stars that comprise the EXPRESS and PPPS common sample. Periodogram analyses of stellar activity indicators present possible peaks at frequencies close to proposed Doppler signals in at least two planetary systems, HIP24275 and HIP90988, calling for more long-term activity studies of giant stars. Even disregarding these possible false-positives, extrapolation leads to a fraction of 25-30% of low-luminosity giant stars hosting planets. We find the mass-function exponentially rises towards the lowest planetary masses, however there exists a ~93% probability that a second population of giant planets with minimum masses between 4-5 MJ, is present, worlds that could have formed by the gravitational collapse of fragmenting proto-planetary disks. Finally, our noise modelling reveals a lack of statistical evidence for the presence of correlated noise...(Abridged)

AGN feedback is invoked as one of the most relevant mechanisms that shape the evolution of galaxies. Our goal is to understand the interplay between AGN feedback and the interstellar medium in M51, a nearby spiral galaxy with a modest AGN and a kpc-scale radio jet expanding through the disc of the galaxy. For that purpose, we combine molecular gas observations in the CO(1-0) and HCN(1-0) lines from the Plateau de Bure interferometer with archival radio, X-ray, and optical data. We show that there is a significant scarcity of CO emission in the ionisation cone, while molecular gas emission tends to accumulate towards the edges of the cone. The distribution and kinematics of CO and HCN line emission reveal AGN feedback effects out to r~500pc, covering the whole extent of the radio jet, with complex kinematics in the molecular gas which displays strong local variations. We propose that this is the result of the almost coplanar jet pushing on molecular gas in different directions as it expands; the effects are more pronounced in HCN than in CO emission, probably as the result of radiative shocks. Following previous interpretation of the redshifted molecular line in the central 5" as caused by a molecular outflow, we estimate the outflow rates to be Mdot_H2~0.9Msun/yr and Mdot_dense~0.6Msun/yr, which are comparable to the molecular inflow rates (~1Msun/yr); gas inflow and AGN feedback could be mutually regulated processes. The agreement with findings in other nearby radio galaxies suggests that this is not an isolated case, and probably the paradigm of AGN feedback through radio jets, at least for galaxies hosting low-luminosity active nuclei.

Very little is known about the young planet population because the detection of small planets orbiting young stars is obscured by the effects of stellar activity and fast rotation which mask planets within radial velocity and transit data sets. The few planets that have been discovered in young clusters generally orbit stars too faint for any detailed follow-up analysis. Here we present the characterization of a new mini-Neptune planet orbiting the bright (V=9) and nearby K2 dwarf star, HD 18599. The planet candidate was originally detected in TESS light curves from Sectors 2, 3, 29, and 30, with an orbital period of 4.138~days. We then used HARPS and FEROS radial velocities, to find the companion mass to be 25.5$\pm$4.6~M$_\oplus$. When we combine this with the measured radius from TESS, of 2.70$\pm$0.05~R$_\oplus$, we find a high planetary density of 7.1$\pm$1.4~g cm$^{-3}$. The planet exists on the edge of the Neptune Desert and is the first young planet (300 Myr) of its type to inhabit this region. Structure models argue for a bulk composition to consist of 23% H$_2$O and 77% Rock and Iron. Future follow-up with large ground- and space-based telescopes can enable us to begin to understand in detail the characteristics of young Neptunes in the galaxy.

We report the discovery of planetary companions orbiting four low-luminosity giant stars with M$_\star$ between 1.04 and 1.39 M$_\odot$. All four host stars have been independently observed by the EXoPlanets aRound Evolved StarS (EXPRESS) program and the Pan-Pacific Planet Search (PPPS). The companion signals were revealed by multi-epoch precision radial velocities obtained during nearly a decade. The planetary companions exhibit orbital periods between $\sim$ 1.2 and 7.1 years, minimum masses of m$_{\rm p}$sini $\sim$ 1.8-3.7 M$_{jup}$ and eccentricities between 0.08 and 0.42. Including these four new systems, we have detected planetary companions to 11 out of the 37 giant stars that are common targets between the EXPRESS and PPPS. After excluding four compact binaries from the common sample, we obtained a fraction of giant planets (m$_{\rm p} \gtrsim$ 1-2 M$\_{jup}$) orbiting within 5 AU from their parent star of $f = 33.3^{+9.0}_{-7.1} \%$. This fraction is significantly higher than that previously reported in the literature by different radial velocity surveys. Similarly, planet formation models under predict the fraction of gas giant around stars more massive than the Sun.

We observed a transit of WASP-166 b using nine NGTS telescopes simultaneously with TESS observations of the same transit. We achieved a photometric precision of 152 ppm per 30 minutes with the nine NGTS telescopes combined, matching the precision reached by TESS for the transit event around this bright (T=8.87) star. The individual NGTS light curve noise is found to be dominated by scintillation noise and appears free from any time-correlated noise or any correlation between telescope systems. We fit the NGTS data for $T_C$ and $R_p/R_*$. We find $T_C$ to be consistent to within 0.25$\sigma$ of the result from the TESS data, and the difference between the TESS and NGTS measured $R_p/R_*$ values is 0.9$\sigma$. This experiment shows that multi-telescope NGTS photometry can match the precision of TESS for bright stars, and will be a valuable tool in refining the radii and ephemerides for bright TESS candidates and planets. The transit timing achieved will also enable NGTS to measure significant transit timing variations in multi-planet systems.

The MaNGA Survey (Mapping Nearby Galaxies at Apache Point Observatory) is one of three core programs in the Sloan Digital Sky Survey IV. It is obtaining integral field spectroscopy (IFS) for 10K nearby galaxies at a spectral resolution of R~2000 from 3,622-10,354A. The design of the survey is driven by a set of science requirements on the precision of estimates of the following properties: star formation rate surface density, gas metallicity, stellar population age, metallicity, and abundance ratio, and their gradients; stellar and gas kinematics; and enclosed gravitational mass as a function of radius. We describe how these science requirements set the depth of the observations and dictate sample selection. The majority of targeted galaxies are selected to ensure uniform spatial coverage in units of effective radius (Re) while maximizing spatial resolution. About 2/3 of the sample is covered out to 1.5Re (Primary sample), and 1/3 of the sample is covered to 2.5Re (Secondary sample). We describe the survey execution with details that would be useful in the design of similar future surveys. We also present statistics on the achieved data quality, specifically, the point spread function, sampling uniformity, spectral resolution, sky subtraction, and flux calibration. For our Primary sample, the median r-band signal-to-noise ratio is ~73 per 1.4A pixel for spectra stacked between 1-1.5 Re. Measurements of various galaxy properties from the first year data show that we are meeting or exceeding the defined requirements for the majority of our science goals.

We present post-cryogenic Spitzer imaging at 3.6 and 4.5 micron with the Infrared Array Camera (IRAC) of the Spitzer/HETDEX Exploratory Large-Area (SHELA) survey. SHELA covers $\sim$deg$^2$ of the Sloan Digital Sky Survey "Stripe 82" region, and falls within the footprints of the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX) and the Dark Energy Survey. The HETDEX blind R $\sim$ 800 spectroscopy will produce $\sim$ 200,000 redshifts from the Lyman-$\alpha$ emission for galaxies in the range 1.9 < z < 3.5, and an additional $\sim$200,000 redshifts from the [OII] emission for galaxies at z < 0.5. When combined with deep ugriz images from the Dark Energy Camera, K-band images from NEWFIRM, and other ancillary data, the IRAC photometry from Spitzer will enable a broad range of scientific studies of the relationship between structure formation, galaxy stellar mass, halo mass, AGN, and environment over a co-moving volume of $\sim$0.5 Gpc$^3$ at 1.9 < z < 3.5. Here, we discuss the properties of the SHELA IRAC dataset, including the data acquisition, reduction, validation, and source catalogs. Our tests show the images and catalogs are 80% (50%) complete to limiting magnitudes of 22.0 (22.6) AB mag in the detection image, which is constructed from the weighted sum of the IRAC 3.6 and 4.5 micron images. The catalogs reach limiting sensitivities of 1.1 $\mu$Jy at both 3.6 and 4.5 micron (1$\sigma$, for R=2 arcsec circular apertures). As a demonstration of science, we present IRAC number counts, examples of highly temporally variable sources, and galaxy surface density profiles of rich galaxy clusters. In the spirit of Spitzer Exploratory programs we provide all images and catalogs as part of the publication.

Mapping Nearby Galaxies at Apache Point Observatory (MaNGA), one of three core programs in the Sloan Digital Sky Survey-IV (SDSS-IV), is an integral-field spectroscopic (IFS) survey of roughly 10,000 nearby galaxies. It employs dithered observations using 17 hexagonal bundles of 2 arcsec fibers to obtain resolved spectroscopy over a wide wavelength range of 3,600-10,300A. To map the internal variations within each galaxy, we need to perform accurate {\it spectral surface photometry}, which is to calibrate the specific intensity at every spatial location sampled by each individual aperture element of the integral field unit. The calibration must correct only for the flux loss due to atmospheric throughput and the instrument response, but not for losses due to the finite geometry of the fiber aperture. This requires the use of standard star measurements to strictly separate these two flux loss factors (throughput versus geometry), a difficult challenge with standard single-fiber spectroscopy techniques due to various practical limitations. Therefore, we developed a technique for spectral surface photometry using multiple small fiber-bundles targeting standard stars simultaneously with galaxy observations. We discuss the principles of our approach and how they compare to previous efforts, and we demonstrate the precision and accuracy achieved. MaNGA's relative calibration between the wavelengths of H$\alpha$ and H$\beta$ has a root-mean-square (RMS) of 1.7%, while that between [NII] $\lambda$6583A and [OII] $\lambda$3727A has an RMS of 4.7%. Using extinction-corrected star formation rates and gas-phase metallicities as an illustration, this level of precision guarantees that flux calibration errors will be sub-dominant when estimating these quantities. The absolute calibration is better than 5% for more than 89% of MaNGA's wavelength range.

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