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

Through the years numerous attempts have been made to connect the phenomenology and physics of mass accretion onto stellar-mass and super-massive black holes in a scale-invariant fashion. In this paper, we explore this connection at the radiatively-efficient (and non-jetted) end of accretion modes by comparing the relationship between the luminosity of the accretion disk and corona in the two source classes. We analyse 458 RXTE-PCA archival observations of the X-ray binary (XRB) GX339-4 focusing on the soft and soft-intermediate states, which have been suggested to be analogous to radiatively efficient AGN. The observed scatter in the $\log L_{disk}-\log L_{corona}$ relationship of GX339-4 is high ($\sim0.43\,$dex) and significantly larger than in a representative sample of radiatively-efficient, non- or weakly-jetted AGN ($\sim0.30\,$dex). On the face of it, this would appear contrary to the hypothesis that the systems simply scale with mass. On the other hand we also find that GX339-4 and our AGN sample show different $\dot{m}$ and $\Gamma$ distributions, with the latter being broader in GX339-4 (dispersion of $\sim0.16$ cf. $\sim0.08$ for AGN). GX339-4 also shows an overall softer slope, with mean $\sim2.20$ as opposed to $\sim2.07$ for the AGN sample. Remarkably, once similarly broad $\Gamma$ and $\dot{m}$ distributions are selected, the AGN sample overlaps nicely with GX339-4 observations in the mass-normalised $\log L_{disk}-\log L_{corona}$ plane, with a scatter of $\sim0.30-0.33\,$dex. This indicates that a mass-scaling of properties might hold after all, with our results being consistent with the disk-corona systems in AGN and XRBs exhibiting the same physical processes, albeit under different conditions for instance in terms of temperature, optical depth and/or electron energy distribution in the corona, heating-cooling balance, coronal geometry and/or black hole spin.

We look to provide a detailed description of the SPectroscopic IDentification of ERosita Sources (SPIDERS) survey, an SDSS-IV programme aimed at obtaining spectroscopic classification and redshift measurements for complete samples of sufficiently bright X-ray sources. We describe the SPIDERS X-ray Point Source Spectroscopic Catalogue, considering its store of 11,092 observed spectra drawn from a parent sample of 14,759 ROSAT and XMM sources over an area of 5,129 deg$^2$ covered in SDSS-IV by the eBOSS survey. This programme represents the largest systematic spectroscopic observation of an X-ray selected sample. A total of 10,970 (98.9\%) of the observed objects are classified and 10,849 (97.8\%) have secure redshifts. The majority of the spectra (10,070 objects) are active galactic nuclei (AGN), 522 are cluster galaxies, and 294 are stars. The observed AGN redshift distribution is in good agreement with simulations based on empirical models for AGN activation and duty cycle. Forming composite spectra of type 1 AGN as a function of the mass and accretion rate of their black holes reveals systematic differences in the H-beta emission line profiles. This study paves the way for systematic spectroscopic observations of sources that are potentially to be discovered in the upcoming eROSITA survey over a large section of the sky.

We present a statistical analysis of the optical properties of an X-ray selected Type 1 AGN sample, using high signal to-noise ratio ($S/N>20$) spectra of the counterparts of the ROSAT/2RXS sources in the footprint of the SDSS-IV/SPIDERS (Spectroscopic IDentification of eROSITA Sources) programme. The final sample contains 2100 sources. It significantly extends the redshift and luminosity ranges ($\rm z \sim 0.01 - 0.80$ and $\rm L_{0.1-2.4 \,keV}$ ~ $\rm 2.0 \times 10^{41}$ - $\rm 1.0 \times 10^{46} \, erg s^{-1}$ ) used so far in this kind of analysis. By means of a Principal Component Analysis, we derive Eigenvector (EV) 1 and 2 in an eleven dimensional optical and X-ray parameter space, which are consistent with previous results. The validity of the correlations of the Eddington ratio $\rm L/L_{Edd}$ with EV1 and the black hole mass with EV2 are strongly confirmed These results imply that $\rm L/L_{Edd}$ and black hole mass are related to the diversity of the optical properties of Type 1 AGN. Investigating the relation of the width and asymmetry of $\rm H\beta$ and the relative strength of the iron emission $\rm r_{FeII}$ , we show that our analysis supports the presence of a distinct kinematic region: the Very Broad Line Region. Furthermore, comparing sources with a red-asymmetric broad $\rm H\beta$ emission line to sources for which it is blue-asymmetric, we find an intriguing difference in the correlation of the FeII and the continuum emission strengths. We show that this contrasting behaviour is consistent with a flattened, stratified model of the Broad Line Region, in which the FeII emitting region is shielded from the central source.

The correlation observed between monochromatic X-ray and UV luminosities in radiatively-efficient active galactic nuclei (AGN) lacks a clear theoretical explanation despite being used for many applications. Such a correlation, with its small intrinsic scatter and its slope that is smaller than unity in log space, represents the compelling evidence that a mechanism regulating the energetic interaction between the accretion disk and the X-ray corona must be in place. This ensures that going from fainter to brighter sources the coronal emission increases less than the disk emission. We discuss here a self-consistently coupled disk-corona model that can identify this regulating mechanism in terms of modified viscosity prescriptions in the accretion disk. The model predicts a lower fraction of accretion power dissipated in the corona for higher accretion states. We then present a quantitative observational test of the model using a reference sample of broad-line AGN and modeling the disk-corona emission for each source in the $L_X-L_{UV}$ plane. We used the slope, normalization, and scatter of the observed relation to constrain the parameters of the theoretical model. For non-spinning black holes and static coronae, we find that the accretion prescriptions that match the observed slope of the $L_X-L_{UV}$ relation produce X-rays that are too weak with respect to the normalization of the observed relation. Instead, considering moderately-outflowing Comptonizing coronae and/or a more realistic high-spinning black hole population significantly relax the tension between the strength of the observed and modeled X-ray emission, while also predicting very low intrinsic scatter in the $L_X-L_{UV}$ relation. In particular, this latter scenario traces a known selection effect of flux-limited samples that preferentially select high-spinning, hence brighter, sources.

The role played by the intergalactic medium (IGM) in the X-ray absorption towards high-redshift sources has recently drawn more attention in spectral analysis studies. Here, we study the X-ray absorption towards 15 flat-spectrum radio quasars at $z>2$, relying on high counting statistic ($\gtrsim10\,000$ photons) provided by XMM-Newton, with additional NuSTAR (and simultaneous Swift-XRT) observations when available. Blazars can be confidently considered to have negligible X-ray absorption along the line of sight within the host galaxy, likely swept by the kpc-scale relativistic jet. This makes our sources ideal for testing the absorption component along the IGM. Our new approach is to revisit the origin of the soft X-ray spectral hardening observed in high-z blazars in terms of X-ray absorption occurring along the IGM, with the help of a low-z sample used as comparison. We verify that the presence of absorption in excess of the Galactic value is the preferred explanation to explain the observed hardening, while intrinsic energy breaks, predicted by blazars' emission models, can easily occur out of the observing energy band in most sources. First, we perform an indirect analysis comparing the inferred amount of absorption in excess of the Galactic value with a simulated IGM absorption contribution, that increases with redshift and includes both a minimum component from diffuse IGM metals, and the additional contribution of discrete denser intervening regions. Then, we directly investigate the warm-hot IGM with a spectral model on the best candidates of our sample, obtaining an average IGM density of $n_0=1.01^{+0.53}_{-0.72}\times10^{-7}$ cm$^{-3}$ and temperature of $\log(T/\text{K})=6.45^{+0.51}_{-2.12}$. A more dedicated study is currently beyond reach, but our results can be used as a stepping stone for future more accurate analysis, involving Athena.

We study the X-ray absorption of a complete sample of 99 bright Swift gamma-ray bursts. Over the last few years, a strong correlation between the intrinsic X-ray absorbing column density (N_H(z)) and the redshift was found. This absorption excess in high-z GRBs is now thought to be due to the overlooked contribution of the absorption along the intergalactic medium, by means of both intervening objects and the diffuse warm-hot intergalactic medium along the line of sight. In this work we neglect the absorption along the IGM, because our purpose is to study the eventual effect of a radical change in the Galactic absorption model on the N_H(z) distribution. Therefore, we derive the intrinsic absorbing column densities using two different Galactic absorption models, the Leiden Argentine Bonn HI survey and the more recent model including molecular hydrogen. We find that, if on the one hand the new Galactic model considerably affects the single column density values, on the other hand there is no drastic change in the distribution as a whole. It becomes clear that the contribution of Galactic column densities alone, no matter how improved, is not sufficient to change the observed general trend and it has to be considered as a second-order correction. The cosmological increase of N_H(z) as a function of redshift persists and, in order to explain the observed distribution, it is necessary to include the contribution of both the diffuse intergalactic medium and the intervening systems along the line of sight of the GRBs.