SPIDERS (The SPectroscopic IDentification of eROSITA Sources) is a large spectroscopic programme for X-ray selected galaxy clusters as part of the Sloan Digital Sky Survey-IV (SDSS-IV). We describe the final dataset in the context of SDSS Data Release 16 (DR16): the survey overall characteristics, final targeting strategies, achieved completeness and spectral quality, with special emphasis on its use as a galaxy cluster sample for cosmology applications. SPIDERS now consists of about 27,000 new optical spectra of galaxies selected within 4,000 photometric red sequences, each associated with an X-ray source. The excellent spectrograph efficiency and a robust analysis pipeline yield a spectroscopic redshift measurement success rate exceeding 98%, with a median velocity accuracy of 20 km s$^{-1}$ (at $z=0.2$). Using the catalogue of 2,740 X-ray galaxy clusters confirmed with DR16 spectroscopy, we reveal the three-dimensional map of the galaxy cluster distribution in the observable Universe up to $z\sim0.6$. We highlight the homogeneity of the member galaxy spectra among distinct regions of the galaxy cluster phase space. Aided by accurate spectroscopic redshifts and by a model of the sample selection effects, we compute the galaxy cluster X-ray luminosity function and we present its lack of evolution up to $z=0.6$. Finally we discuss the prospects of forthcoming large multiplexed spectroscopic programmes dedicated to follow up the next generation of all-sky X-ray source catalogues.

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

Large area catalogs of galaxy clusters constructed from ROSAT All Sky Survey provide the base for our knowledge on the population of clusters thanks to the long-term multiwavelength efforts on their follow-up. Advent of large area photometric surveys superseding in depth previous all-sky data allows us to revisit the construction of X-ray cluster catalogs, extending the study to lower cluster masses and to higher redshifts and to provide the modelling of the selection function. We perform a wavelet detection of X-ray sources and make extensive simulations of the detection of clusters in the RASS data. We assign an optical richness to each of the 24,788 detected X-ray sources in the 10,382 square degrees of SDSS BOSS area, using redMaPPer version 5.2. We name this survey COnstrain Dark Energy with X-ray (CODEX) clusters. We show that there is no obvious separation of sources on galaxy clusters and AGN, based on distribution of systems on their richness. This is a combination of increasing number of galaxy groups and their selection as identification of an X-ray sources either by chance or due to groups hosting an AGN. To clean the sample, we use a cut on the optical richness at the level corresponding to the 10\% completeness of the survey and include it into the modelling of cluster selection function. We present the X-ray catalog extending to a redshift of 0.6 down to X-ray fluxes of $10^{-13}$ ergs s$^{-1}$ cm$^{-2}$. We provide the modelling of the sample selection and discuss the redshift evolution of the high end of the X-ray luminosity function (XLF). Our results on $z<0.3$ XLF are in agreement with previous studies, while we provide new constraints on the $0.3<z<0.6$ XLF. We find a lack of strong redshift evolution of the XLF, provide exact modeling of the effect of low number statistics and AGN contamination and present the resulting constraints on the flat $\Lambda$CDM.

We perform the calibration of the X-ray luminosity--mass scaling relation on a sample of 344 CODEX clusters with $z<0.66$ using the dynamics of their member galaxies. Spectroscopic follow-up measurements have been obtained from the SPIDERS survey, leading to a sample of 6,658 red member galaxies. We use the Jeans equation to calculate halo masses, assuming an NFW mass profile and analyzing a broad range of anisotropy profiles. With a scaling relation of the form $L_{\rm{X}} \propto \text{A}_{\rm{X}}M_{\text{200c}}^{\text{B}_{\rm{X}}} E(z)^2 (1+z)^{\gamma_{\rm{X}}}$, we find best fit parameters $ \text{A}_{\rm{X}}=5.7^{+0.4}_{-0.5} (\pm0.6)\times10^{43}\,\mathrm{erg\,s^{-1}}$, $\text{B}_{\rm{X}}=2.5\pm0.2(\pm0.07)$, $\gamma_{\rm{X}}=-2.6^{+1.1}_{-1.2}(\pm0.74)$, where we include systematic uncertainties in parentheses and for a pivot mass and redshift of $3\times10^{14}M_\odot$ and 0.16, respectively. We compare our constraints with previous results, and we combine our sample with the SPT SZE--selected cluster subsample observed with \XMM\, extending the validity of our results to a wider range of redshifts and cluster masses.

We present the brightest cluster galaxies (BCGs) catalog for SPectroscoic IDentification of eROSITA Sources (SPIDERS) DR14 cluster program value-added catalog. We list the 416 BCGs identified as part of this process, along with their stellar mass, star formation rates, and morphological properties. We identified the BCGs based on the available spectroscopic data from SPIDERS and photometric data from SDSS. We computed stellar masses and SFRs of the BCGs on the basis of SDSS, WISE, and GALEX photometry using spectral energy distribution fitting. Morphological properties for all BCGs were derived by Sersic profile fitting using the software package SIGMA in different optical bands (g,r,i). We combined this catalog with the BCGs of galaxy groups and clusters extracted from the deeper AEGIS, CDFS, COSMOS, XMM-CFHTLS, and XMM-XXL surveys to study the stellar mass - halo mass relation using the largest sample of X-ray groups and clusters known to date. This result suggests that the mass growth of the central galaxy is controlled by the hierarchical mass growth of the host halo. We find a strong correlation between the stellar mass of BCGs and the mass of their host halos. This relation shows no evolution since z $\sim$ 0.65. We measure a mean scatter of 0.21 and 0.25 for the stellar mass of BCGs in a given halo mass at low ( $0.1<z < 0.3$ ) and high ( $0.3<z<0.65$ ) redshifts, respectively. We further demonstrate that the BCG mass is covariant with the richness of the host halos in the very X-ray luminous systems. We also find evidence that part of the scatter between X-ray luminosity and richness can be reduced by considering stellar mass as an additional variable.

We use galaxy dynamical information to calibrate the richness-mass scaling relation of a sample of 428 galaxy clusters that are members of the CODEX sample with redshifts up to z~0.7. These clusters were X-ray selected using the ROSAT All-Sky Survey (RASS), cross-matched to associated systems in the redMaPPer catalog from the Sloan Digital Sky Survey. The spectroscopic sample we analyze was obtained in the SPIDERS program and contains ~7800 red member galaxies. Adopting NFW mass and galaxy density profiles and a broad range of orbital anisotropy profiles, we use the Jeans equation to calculate halo masses. Modeling the scaling relation as $\lambda \propto \text{A}_{\lambda} {M_{\text{200c}}}^{\text{B}_{\lambda}} ({1+z})^{\gamma_{\lambda}}$, we find the parameter constraints $\text{A}_{\lambda}=38.6^{+3.1}_{-4.1}\pm3.9$, $\text{B}_{\lambda}=0.99^{+0.06}_{-0.07}\pm0.04$, and $\gamma_{\lambda}=-1.13^{+0.32}_{-0.34}\pm0.49$. We find good agreement with previously published mass trends with the exception of those from stacked weak lensing analyses. We note that although the lensing analyses failed to account for the Eddington bias, this is not enough to explain the differences. We suggest that differences in the levels of contamination between pure redMaPPer and RASS+redMaPPer samples could well contribute to these differences. The redshift trend we measure is more negative than but statistically consistent with previous results. We suggest that our measured redshift trend reflects a change in the cluster galaxy red sequence fraction with redshift, noting that the trend we measure is consistent with but somewhat stronger than an independently measured redshift trend in the red sequence fraction. We also examine the impact of a plausible model of correlated scatter in X-ray luminosity and optical richness, showing it has negligible impact on our results.