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.

Large multi-object spectroscopic surveys require automated algorithms to optimise their observing strategy. One of the most ambitious upcoming spectroscopic surveys is the 4MOST survey. The 4MOST survey facility is a fibre-fed spectroscopic instrument on the VISTA telescope with a large enough field of view to survey a large fraction of the southern sky within a few years. Several Galactic and extragalactic surveys will be carried out simultaneously, so the combined target density will strongly vary. In this paper, we describe a new tiling algorithm that can naturally deal with the large target density variations on the sky and which automatically handles the different exposure times of targets. The tiling pattern is modelled as a marked point process, which is characterised by a probability density that integrates the requirements imposed by the 4MOST survey. The optimal tilling pattern with respect to the defined model is estimated by the tiles configuration that maximises the proposed probability density. In order to achieve this maximisation a simulated annealing algorithm is implemented. The algorithm automatically finds an optimal tiling pattern and assigns a tentative sky brightness condition and exposure time for each tile, while minimising the total execution time that is needed to observe the list of targets in the combined input catalogue of all surveys. Hence, the algorithm maximises the long-term observing efficiency and provides an optimal tiling solution for the survey. While designed for the 4MOST survey, the algorithm is flexible and can with simple modifications be applied to any other multi-object spectroscopic survey.

Understanding the links between the activity of supermassive black holes (SMBH) at the centres of galaxies and their host dark matter haloes is a key question in modern astrophysics. The final data release of the SDSS-IV eBOSS provides the largest contemporary spectroscopic sample of galaxies and QSOs. Using this sample and covering the redshift interval $z=0.7-1.1$, we have measured the clustering properties of the eBOSS QSOs, Emission Line Galaxies (ELGs) and Luminous Red Galaxies (LRGs). We have also measured the fraction of QSOs as a function of the overdensity defined by the galaxy population. Using these measurements, we investigate how QSOs populate and sample the galaxy population, and how the host dark-matter haloes of QSOs sample the underlying halo distribution. We find that the probability of a galaxy hosting a QSO is independent of the host dark matter halo mass of the galaxy. We also find that about 60\% of eBOSS QSOs are hosted by LRGs and about 20-40\% of QSOs are hosted by satellite galaxies. We find a slight preference for QSOs to populate satellite galaxies over central galaxies. This is connected to the host halo mass distribution of different types of galaxies. Based on our analysis, QSOs should be hosted by a very broad distribution of haloes, and their occurrence should be modulated only by the efficiency of galaxy formation processes.

We develop a new Multi-Tracer Halo Occupation Distribution (\texttt{MTHOD}) framework for the galaxy distribution and apply it to the extended Baryon Oscillation Spectroscopic Survey (eBOSS) final data between $z=0.7-1.1$. We obtain a best fit \mthod\, for each tracer and describe the host halo properties of these galaxies. The mean halo masses for LRGs, ELGs and QSOs are found to be $1.9 \times 10^{13} \msolaroh$, $1.1 \times 10^{12} \msolaroh$ and $5 \times 10^{12} \msolaroh$ respectively in the eBOSS data. We use the \texttt{MTHOD} framework to create mock galaxy catalogues and predict auto- and cross-correlation functions for all the tracers. Comparing these results with data, we investigate galactic conformity, the phenomenon whereby the properties of neighbouring galaxies are mutually correlated in a manner that is not captured by the basic halo model. We detect \textsl{1-halo} conformity at more than 3$\sigma$ statistical significance, while obtaining upper limit on \textsl{2-halo} conformity. We also look at the environmental dependence of the galaxy quenching efficiency and find that halo mass driven quenching successfully explains the behaviour in high density regions, but it fails to describe the quenching efficiency in low density regions. In particular, we show that the quenching efficiency in low density filaments is higher in the observed data, as compared to the prediction of the \mthod\ with halo mass driven quenching. The mock galaxy catalogue constructed in this paper is publicly available on https://www.roe.ac.uk/~salam/MTHOD/ .

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 use three semi-analytic models (SAMs) of galaxy formation and evolution, run on the same 1$h^{-1}$Gpc MultiDark Planck2 cosmological simulation, to investigate the properties of [OII] emission line galaxies in the redshift range $0.6<z<1.2$. We compare model predictions with different observational data sets, including DEEP2--Firefly galaxies with absolute magnitudes. We estimate the [OII] luminosity, L[OII], using simple relations derived both from the models and observations and also using a public code. This code ideally uses as input instantaneous star formation rates (SFRs), which are only provided by one of the SAMs under consideration. We use this SAM to study the feasibility of inferring galaxies' L[OII] for models that only provide average SFRs. We find that the post-processing computation of L[OII] from average SFRs is accurate for model galaxies with dust attenuated L[OII]$\lesssim10^{42.2}$erg s$^{-1}$ ($<5\%$ discrepancy). We also explore how to derive the [OII] luminosity from simple relations using global properties usually output by SAMs. Besides the SFR, the model L[OII] is best correlated with the observed-frame $u$ and $g$ broad-band magnitudes. These correlations have coefficients (r-values) above 0.64 and a dispersion that varies with L[OII]. We use these correlations and an observational one based on SFR and metallicity to derive L[OII]. These relations result in [OII] luminosity functions and halo occupation distributions with shapes that vary depending on both the model and the method used. Nevertheless, for all the considered models, the amplitude of the clustering at scales above 1$h^{-1}$Mpc remains unchanged independently of the method used to derive L[OII].

Context. Several new multi-object spectrographs are currently planned or under construction that are capable of observing thousands of Galactic and extragalactic objects simultaneously. Aims. In this paper we present a probabilistic fibre-to-target assignment algorithm that takes spectrograph targeting constraints into account and is capable of dealing with multiple concurrent surveys. We present this algorithm using the 4-metre Multi-Object Spectroscopic Telescope (4MOST) as an example. Methods. The key idea of the proposed algorithm is to assign probabilities to fibre-target pairs. The assignment of probabilities takes the fibre positioner's capabilities and constraints into account. Additionally, these probabilities include requirements from surveys and take the required exposure time, number density variation, and angular clustering of targets across each survey into account. The main advantage of a probabilistic approach is that it allows for accurate and easy computation of the target selection function for the different surveys, which involves determining the probability of observing a target, given an input catalogue. Results. The probabilistic fibre-to-target assignment allows us to achieve maximally uniform completeness within a single field of view. The proposed algorithm maximises the fraction of successfully observed targets whilst minimising the selection bias as a function of exposure time. In the case of several concurrent surveys, the algorithm maximally satisfies the scientific requirements of each survey and no specific survey is penalised or prioritised. Conclusions. The algorithm presented is a proposed solution for the 4MOST project that allows for an unbiased targeting of many simultaneous surveys. With some modifications, the algorithm may also be applied to other multi-object spectroscopic surveys.

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.

Current and future cosmological surveys are targeting star-forming galaxies at z~1 with nebular emission lines. In this work we use a state-of-the-art model of galaxy formation and evolution to explore the large scale environment of star-forming emission line galaxies (ELGs). Model galaxies are selected with cuts such that the samples can be directly compared with the DEEP2, VVDS, eBOSS-SGC and DESI surveys. Their large scale environment have been classified using a velocity-shear-tensor and a tidal-tensor algorithms. Half of model ELGs live in filaments and about a third in sheets. Model ELGs in knots have the largest satellite fractions. We find that the shape of the mean halo occupation distribution of model ELGs varies widely for different large scale environments. To put these results in context, we have compared fixed number density samples of galaxies, generated by either imposing an extra cut in stellar mass, star formation rate (SFR) or L[OII] to the ELGs, or by imposing a single cut in these quantities to the total model galaxy population. For the fixed number density samples, ELGs are close to L[OII] and SFR selected samples for densities above $10^{-4.2}h^{3}{\rm Mpc}^{-3}$. ELGs with an extra cut in stellar mass to fix their number density, present differences in sheets and knots with respect to the other samples. ELGs, SFR and L[OII] selected samples with equal number density, have similar large scale bias but their clustering below separations of $1h^{-1}$Mpc is different.

We analyse publicly available, individual spectra of four, massive ($M>10^{11}M_{\odot}$) early-type galaxies with redshifts in the range 1.4 < z < 2 to determine their stellar content, extending our previous work up to z~2. The wide wavelength range of the VLT/X-Shooter spectroscopic data in the UV-Optical-NIR arms along with the availability of spectro-photometry allows us to explore different techniques to obtain the stellar population properties, namely through age/metallicity sensitive spectral indices, full spectral fitting and broad-band photometric fitting. Moreover, together with the widely used optical Lick indices we consider further indices in the UV rest-frame, and demonstrate that UV indices significantly help the accuracy of the resulting population parameters. We find galaxy ages ranging from 0.2 to 4 Gyr, where the oldest galaxy is found at the lowest redshift, with an excellent agreement between ages determined via indices, full spectral fitting or broad-band colours. These ages are in perfect agreement with ages of local galaxies at the same velocity dispersion when we assume pure passive evolution. Total metallicities derived from indices show some scatter (between less than half-solar to very high values, ([Z/H]~0.6]). We speculate on possible mechanisms explaining these values, but given the sample size and low S/N of the spectra no conclusion can be made. Indices in the UV-rest frame generally lead to similar conclusions as optical indices. For the oldest galaxy (4 Gyr) we show that its UV-indices can only be explained by stellar population models including a UV contribution from old stellar populations, suggesting that old, UV bright populations start to inhabit mature galaxies of a few Gyr of age. This is the highest redshift (z~1.4) detection of the UV-upturn up to date.

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 spectra of emission-line galaxies (ELGs) from the extended Baryon Oscillation Spectroscopic Survey (eBOSS) of the Sloan Digit Sky Survey (SDSS) are used to study the mass-metallicity relation (MZR) at $z\sim0.8$. The selected sample contains about 180,000 massive star-forming galaxies with $0.6 < z < 1.05$ and $9 < {\rm log}(M_{\star}/M_{\odot}) < 12$. The spectra are stacked in bins of different parameters including redshift, stellar mass, star formation rate (SFR), specific star formation rate (sSFR), half-light radius, mass density, and optical color. The average MZR at $z\sim0.83$ has a downward evolution in the MZR from local to high-redshift universe, which is consistent with previous works. At a specified stellar mass, galaxies with higher SFR/sSFR and larger half-light radius have systematically lower metallicity. This behavior is reversed for galaxies with larger mass density and optical color. Among the above physical parameters, the MZR has the most significant dependency on SFR. Our galaxy sample at $0.6<z<1.05$ approximately follows the fundamental metallicity relation (FMR) in the local universe, although the sample inhomogeneity and incompleteness might have effect on our MZR and FMR.

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 present the measurements and modelling of the small-to-intermediate scale (0.1--25 Mpc/h) projected and three-dimensional (3D) redshift-space two-point correlation functions (2PCFs) of local galaxies in the Sloan Digital Sky Survey (SDSS) Data Release 7. We find a clear dependence of galaxy clustering on luminosity in both projected and redshift spaces, generally being stronger for more luminous samples. The measurements are successfully interpreted within the halo occupation distribution (HOD) framework with central and satellite velocity bias parameters to describe galaxy kinematics inside haloes and to model redshift-space distortion (RSD) effects. In agreement with previous studies, we find that more luminous galaxies reside in more massive haloes. Including the redshift-space 2PCFs helps tighten the HOD constraints. Moreover, we find that luminous central galaxies are not at rest at the halo centres, with the velocity dispersion about 30% that of the dark matter. Such a relative motion may reflect the consequence of galaxy and halo mergers, and we find that central galaxies in lower mass haloes tend to be more relaxed with respect to their host haloes. The motion of satellite galaxies in luminous samples is consistent with their following that of the dark matter. For faint samples, satellites tends to have slower motion, with velocity dispersion inside haloes about 85% that of the dark matter. We discuss possible applications of the velocity bias constraints on studying galaxy evolution and cosmology. In the appendix, we characterize the distribution of galaxy redshift measurement errors, which is well described by a Gaussian-convolved double exponential distribution.

Classification of intermediate redshift ($z$ = 0.3--0.8) emission line galaxies as star-forming galaxies, composite galaxies, active galactic nuclei (AGN), or low-ionization nuclear emission regions (LINERs) using optical spectra alone was impossible because the lines used for standard optical diagnostic diagrams: [NII], H$\alpha$, and [SII] are redshifted out of the observed wavelength range. In this work, we address this problem using four supervised machine learning classification algorithms: $k$-nearest neighbors (KNN), support vector classifier (SVC), random forest (RF), and a multi-layer perceptron (MLP) neural network. For input features, we use properties that can be measured from optical galaxy spectra out to $z < 0.8$---[OIII]/H$\beta$, [OII]/H$\beta$, [OIII] line width, and stellar velocity dispersion---and four colors ($u-g$, $g-r$, $r-i$, and $i-z$) corrected to $z=0.1$. The labels for the low redshift emission line galaxy training set are determined using standard optical diagnostic diagrams. RF has the best area under curve (AUC) score for classifying all four galaxy types, meaning highest distinguishing power. Both the AUC scores and accuracies of the other algorithms are ordered as MLP$>$SVC$>$KNN. The classification accuracies with all eight features (and the four spectroscopically-determined features only) are 93.4% (92.3%) for star-forming galaxies, 69.4% (63.7%) for composite galaxies, 71.8% (67.3%) for AGNs, and 65.7% (60.8%) for LINERs. The stacked spectrum of galaxies of the same type as determined by optical diagnostic diagrams at low redshift and RF at intermediate redshift are broadly consistent. Our publicly available code (https://github.com/zkdtc/MLC_ELGs) and trained models will be instrumental for classifying emission line galaxies in upcoming wide-field spectroscopic surveys.

This paper presents stellar mass functions and i-band luminosity functions for Sloan Digital Sky Survey (SDSS) galaxies at $i < 21$ using clustering redshifts, from which we also compute targeting completeness measurements for the Baryon Oscillation Spectroscopic Survey (BOSS). Clustering redshifts is a method of obtaining the redshift distribution of a sample of galaxies with only photometric information by measuring the angular crosscorrelation with a spectroscopic sample in different redshift bins. We construct a spectroscopic sample containing data from the BOSS + eBOSS surveys, allowing us to recover redshift distributions from photometric data out to $z\simeq 2.5$. We produce k-corrected i-band luminosity functions and stellar mass functions by applying clustering redshifts to SDSS DR8 galaxies in small bins of colour and magnitude. There is little evolution in the mass function between $0.2 < z < 0.8$, implying the most massive galaxies form most of their mass before $z = 0.8$. These mass functions are used to produce stellar mass completeness estimates for the Baryon Oscillation Spectroscopic Survey (BOSS), giving a stellar mass completeness of $80\%$ above $M_{\star} > 10^{11.4}$ between $0.2 < z < 0.7$, with completeness falling significantly at redshifts higher than 0.7, and at lower masses. Large photometric datasets will be available in the near future (DECaLS, DES, Euclid), so this, and similar techniques will become increasingly useful in order to fully utilise this data.

In the context of the upcoming SRG/eROSITA survey, we present an N-body simulation-based mock catalogue for X-ray selected AGN samples. The model reproduces the observed hard X-ray AGN luminosity function (XLF) and the soft X-ray logN-logS from redshift 0 to 6. The XLF is reproduced to within $\pm5\%$ and the logN-logS to within $\pm20\%$. We develop a joint X-ray -- optical extinction and classification model. We adopt a set of empirical spectral energy distributions to predict observed magnitudes in the UV, optical and NIR. With the latest eROSITA all sky survey sensitivity model, we create a high-fidelity full-sky mock catalogue of X-ray AGN. It predicts their distributions in right ascension, declination, redshift and fluxes. Using empirical medium resolution optical spectral templates and an exposure time calculator, we find that $1.1\times10^6$ ($4\times10^5$) fiber-hours are needed to follow-up spectroscopically from the ground the detected X-ray AGN with an optical magnitude $21<r<22.8$ ($22.8<r<25$) with a 4-m (8-m) class multi-object spectroscopic facility. We find that future clustering studies will measure the AGN bias to the percent level at redshift $z<1.2$ and should discriminate possible scenarios of galaxy-AGN co-evolution. We predict the accuracy to which the baryon acoustic oscillation standard ruler will be measured using X-ray AGN: better than 3\% for AGN between redshift 0.5 to 3 and better than 1\% using the Ly$\alpha$ forest of X-ray QSOs discovered between redshift 2 and 3. eROSITA will provide an outstanding set of targets for future galaxy evolution and cosmological studies.

This work presents the catalogue of optical spectral properties for all X-ray selected SPIDERS active galactic nuclei (AGN) up to SDSS DR14. SPIDERS (SPectroscopic IDentification of eROSITA Sources) is an SDSS-IV programme that is currently conducting optical spectroscopy of the counterparts to the X-ray selected sources detected in the ROSAT all-sky survey and the XMM-Newton slew survey in the footprint of the Extended Baryon Oscillation Spectroscopic Survey (eBOSS). The SPIDERS DR14 sample is the largest sample of X-ray selected AGN with optical spectroscopic follow-up to date. The catalogue presented here is based on a clean sample of 7344 2RXS ($\rm \bar{z}$ = 0.5) and 1157 XMM-Newton slew survey ($\rm \bar{z}$ = 0.4) type 1 AGN with spectroscopic coverage of the H$\rm \beta$ and/or MgII emission lines. Visual inspection results for each object in this sample are available from a combination of literature sources and the SPIDERS group, which provide both reliable redshifts and source classifications. The spectral regions around the H$\rm \beta$ and MgII emission lines have been fit in order to measure both line and continuum properties, estimate bolometric luminosities, and provide black hole mass estimates using the single-epoch (or photoionisation) method. The use of both H$\rm \beta$ and MgII allows black hole masses to be estimated up to z $\rm \simeq$ 2.5. A comparison is made between the spectral properties and black hole mass estimates derived from H$\rm \beta$ and MgII using the subsample of objects which have coverage of both lines in their spectrum. These results have been made publicly available as an SDSS-IV DR14 value added catalogue.

(Abridged) The Maunakea Spectroscopic Explorer (MSE) is an end-to-end science platform for the design, execution and scientific exploitation of spectroscopic surveys. It will unveil the composition and dynamics of the faint Universe and impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe. Major pillars in the science program for MSE include (i) the ultimate Gaia follow-up facility for understanding the chemistry and dynamics of the distant Milky Way, including the outer disk and faint stellar halo at high spectral resolution (ii) galaxy formation and evolution at cosmic noon, via the type of revolutionary surveys that have occurred in the nearby Universe, but now conducted at the peak of the star formation history of the Universe (iii) derivation of the mass of the neutrino and insights into inflationary physics through a cosmological redshift survey that probes a large volume of the Universe with a high galaxy density. MSE is positioned to become a critical hub in the emerging international network of front-line astronomical facilities, with scientific capabilities that naturally complement and extend the scientific power of Gaia, the Large Synoptic Survey Telescope, the Square Kilometer Array, Euclid, WFIRST, the 30m telescopes and many more.

The 4MOST instrument is a multi-object spectrograph that will address Galactic and extragalactic science cases simultaneously by observing targets from a large number of different surveys within each science exposure. This parallel mode of operation and the survey nature of 4MOST require some distinct 4MOST-specific operational features within the overall operations model of ESO. The main feature is that the 4MOST Consortium will deliver, not only the instrument, but also contractual services to the user community, which is why 4MOST is also described as a facility. This white paper concentrates on information particularly useful to answering the forthcoming Call for Letters of Intent.

The 4MOST Cosmology Redshift Survey (CRS) will perform stringent cosmological tests via spectroscopic clustering measurements that will complement the best lensing, cosmic microwave background and other surveys in the southern hemisphere. The combination of carefully selected samples of bright galaxies, luminous red galaxies, emission-line galaxies and quasars, totalling about 8 million objects over the redshift range $z = 0.15$ to $3.5$, will allow definitive tests of gravitational physics. Many key science questions will be addressed by combining CRS spectra of these targets with data from current or future facilities such as the Large Synoptic Survey Telescope, the Square Kilometre Array and the Euclid mission.

The current status of and motivation for the 4MOST Survey Strategy, as developed by the Consortium science team, are presented here. Key elements of the strategy are described, such as sky coverage, number of visits and total exposure times in different parts of the sky, and how to deal with different observing conditions. The task of organising the strategy is not simple, with many different surveys that have vastly different target brightnesses and densities, sample completeness levels, and signal-to-noise requirements. We introduce here a number of concepts that we will use to ensure all surveys are optimised. Astronomers who are planning to submit a Participating Survey proposal are strongly encouraged to read this article and any relevant 4MOST Survey articles in this issue of The Messenger such that they can optimally complement and benefit from the planned surveys of the 4MOST Consortium.

Groups and clusters of galaxies are a current focus of astronomical research owing to their role in determining the environmental effects on galaxies and the constraints they provide to cosmology. The eROSITA X-ray telescope on board the Spectrum Roentgen Gamma observatory will be launched in 2019 and will have completed eight scans of the full sky when 4MOST starts operating. The experiment will detect groups and clusters of galaxies through X-ray emission from the hot intergalactic medium. The purpose of the 4MOST eROSITA Galaxy Cluster Redshift Survey is to provide spectroscopic redshifts of the optical counterparts to the X-ray emission from 40,000 groups and clusters of galaxies so as to perform dynamical estimates of the total mass and to measure the properties of the member galaxies. The survey aims to obtain precise redshift measurements of the photometrically identified brightest cluster galaxies at redshift $z > 0.7$. At lower redshifts ($z < 0.7$), the programme aims to sample over 15 member galaxies per cluster and enable dynamical mass measurements to calibrate the clusters for cosmological experiments. At $z < 0.2$, eROSITA will also detect X-ray emission from galaxy groups and filaments. 4MOST spectroscopic data from the survey will be used for optical identification of galaxy groups down to eROSITA's mass detection limits of $10^{13} M_\odot$, as well as the detection of the largest filaments for pioneering studies of their X-ray emission.

X-ray and mid-infrared emission are signposts of the accretion of matter onto the supermassive black holes that reside at the centres of most galaxies. As a major step towards understanding accreting supermassive black holes and their role in the evolution of galaxies, we will use the 4MOST multi-object spectrograph to provide a highly complete census of active galactic nuclei over a large fraction of the extragalactic sky observed in X-rays by eROSITA that is visible to 4MOST. We will systematically follow up all eROSITA point-like extragalactic X-ray sources (mostly active galactic nuclei), and complement them with a heavily obscured active galactic nuclei selection approach using mid-infrared data from the Wide-field Infrared Survey Explorer (WISE). The X-ray and mid-infrared flux limits of eROSITA and WISE are well matched to the spectroscopic capabilities of a 4-metre-class telescope, allowing us to reach completeness levels of ~80-90% for all X-ray selected active galactic nuclei with fluxes $f_{0.5-2\ {\rm keV}} > 10^{-14}$ erg s$^{-1}$ cm$^{-2}$; this is about a factor of 30 deeper than the ROSAT all-sky survey. With these data we will determine the physical properties (redshift, luminosity, line emission strength, masses, etc.) of up to one million supermassive black holes, constrain their cosmic evolution and clustering properties, and explore the connection between active galactic nuclei and large-scale structure over redshifts $0 \le z \le 6$.

We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolution spectrographs ($R = \lambda/\Delta\lambda \sim 6500$), and 812 fibres transferring light to the high-resolution spectrograph ($R \sim 20\,000$). After a description of the instrument and its expected performance, a short overview is given of its operational scheme and planned 4MOST Consortium science; these aspects are covered in more detail in other articles in this edition of The Messenger. Finally, the processes, schedules, and policies concerning the selection of ESO Community Surveys are presented, commencing with a singular opportunity to submit Letters of Intent for Public Surveys during the first five years of 4MOST operations.

The combination of Galaxy-Galaxy Lensing (GGL) and Redshift Space Distortion of galaxy clustering (RSD) is a privileged technique to test General Relativity predictions, and break degeneracies between the growth rate of structure parameter $f$ and the amplitude of the linear power-spectrum $\sigma_8$. We perform a joint GGL and RSD analysis on 250 sq. degrees using shape catalogues from CFHTLenS and CFHT-Stripe 82, and spectroscopic redshifts from the BOSS CMASS sample. We adjust a model that includes non-linear biasing, RSD and Alcock-Paczynski effects. We find $f(z=0.57) =0.95\pm0.23$, $\sigma_8(z=0.57)=0.55\pm0.07$ and $\Omega_{\rm m} = 0.31\pm0.08$, in agreement with Planck cosmological results 2018. We also estimate the probe of gravity $E_{\rm G} = 0.43\pm0.10$ in agreement with $\Lambda$CDM-GR predictions of $E_{\rm G} = 0.40$. This analysis reveals that RSD efficiently decreases the GGL uncertainty on $\Omega_{\rm m}$ by a factor of 4, and by 30\% on $\sigma_8$. We use an N-body simulation supplemented by an abundance matching prescription for CMASS to build a set of overlapping lensing and clustering mocks. Together with additional spectroscopic data, this helps us to quantify and correct several systematic errors, such as photometric redshifts. We make our mock catalogues available on the Skies and Universe database.

Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We describe how astrophysical observations will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational interactions with the Standard Model, and compact object abundances. Additionally, we highlight theoretical work and experimental/observational facilities that will complement LSST to strengthen our understanding of the fundamental characteristics of dark matter.

The DESI Legacy Imaging Surveys are a combination of three public projects (the Dark Energy Camera Legacy Survey, the Beijing-Arizona Sky Survey, and the Mayall z-band Legacy Survey) that will jointly image approximately 14,000 deg^2 of the extragalactic sky visible from the northern hemisphere in three optical bands (g, r, and z) using telescopes at the Kitt Peak National Observatory and the Cerro Tololo Inter-American Observatory. The combined survey footprint is split into two contiguous areas by the Galactic plane. The optical imaging is conducted using a unique strategy of dynamically adjusting the exposure times and pointing selection during observing that results in a survey of nearly uniform depth. In addition to calibrated images, the project is delivering a catalog, constructed by using a probabilistic inference-based approach to estimate source shapes and brightnesses. The catalog includes photometry from the grz optical bands and from four mid-infrared bands (at 3.4, 4.6, 12 and 22 micorons) observed by the Wide-field Infrared Survey Explorer (WISE) satellite during its full operational lifetime. The project plans two public data releases each year. All the software used to generate the catalogs is also released with the data. This paper provides an overview of the Legacy Surveys project.

We determine the stellar population properties - age, metallicity, dust reddening, stellar mass and the star formation history - for all spectra classified as galaxies that were published by the Sloan Digital Sky Survey (SDSS data release 14) and by the DEEP2 (data release 4) galaxy surveys. We perform full spectral fitting on individual spectra, making use of high spectral resolution stellar population models. Calculations are carried out for several choices of the model input, including three stellar initial mass functions and three input stellar libraries to the models. We study the accuracy of parameter derivation, in particular the stellar mass, as a function of the signal-to-noise of the galaxy spectra. We find that at low redshift, a signal to noise ratio per pixel around 20 (5) allows a statistical accuracy on $\log_{10}(M^{*}/M_{\odot})$ of 0.2 (0.4) dex, for the Chabrier IMF. For the first time, we study DEEP2 galaxies selected by their \OII luminosity in the redshift range $0.83<z<1.03$, finding that they are consistent with a flat number density in stellar mass in the range $10^9<M/M_{\odot}<10^{11.5}$. We find the resulting stellar mass function based on SDSS or eBOSS in agreement with previous studies (Maraston et al. 2013). We publish all catalogs of properties as well as model spectra of the continuum for these galaxies as a value added catalog of the fourteenth data release of the SDSS. This catalog is about twice as large as its predecessors (DR12) and will aid a variety of studies on galaxy evolution and cosmology.

We study the evolution of star-forming galaxies with $10^{10} M_\odot<M_*<10^{11.6} M_\odot$ over the redshift range of 0.7<z<1.2 using the emission line galaxies (ELGs) in the extended Baryon Oscillation Spectroscopic Survey (eBOSS). By applying the incomplete conditional stellar mass function (ICSMF) model proposed in Guo et al., we simultaneously constrain the sample completeness, the stellar--halo mass relation (SHMR) and the quenched galaxy fraction. We obtain the intrinsic stellar mass functions for star-forming galaxies in the redshift bins of 0.7<z<0.8, 0.8<z<0.9, 0.9<z<1.0 and 1.0<z<1.2, as well as the stellar mass function for all galaxies in the redshift bin of 0.7<z<0.8. We find that the eBOSS ELG sample only selects about 1%-10% of the star-forming galaxy population at the different redshifts, with the lower redshift samples to be more complete. There is only weak evolution in the SHMR of the ELGs from z=1.2 to z=0.7, as well as the intrinsic galaxy stellar mass functions for lower-mass galaxies of $M_*<10^{11} M_\odot$. There is significant decrease of the stellar mass function for star-forming galaxies with redshift at the massive end. Our best-fitting models show that the central ELGs at these redshifts live in halos of mass $M\sim10^{12} M_\odot$ while the satellite ELGs occupy slightly more massive halos of $M\sim10^{12.6} M_\odot$. The average satellite fraction of the observed ELGs varies from 13% to 17%, with the galaxy bias increasing from 1.1 to 1.4 from z=0.7 to 1.2.

A semi-empirical model is presented that describes the distribution of Active Galactic Nuclei (AGN) on the cosmic web. It populates dark-matter halos in N-body simulations (MultiDark) with galaxy stellar masses using empirical relations based on abundance matching techniques, and then paints accretion events on these galaxies using state-of-the-art measurements of the AGN occupation of galaxies. The explicit assumption is that the large-scale distribution of AGN is independent of the physics of black-hole fueling. The model is shown to be consistent with current measurements of the two-point correlation function of AGN samples. It is then used to make inferences on the halo occupation of the AGN population. Mock AGN are found in halos with a broad distribution of masses with a mode of $\approx 10^{12}\,h^{-1} \, M_{\odot}$ and a tail extending to cluster-size halos. The clustering properties of the model AGN depend only weakly on accretion luminosity and redshift. The fraction of satellite AGN in the model increases steeply toward more massive halos, in contrast with some recent observational results. This discrepancy, if confirmed, could point to a dependence of the halo occupation of AGN on the physics of black-hole fueling.

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.

Using the X-ray-selected active galactic nuclei (AGN) from the XMM-XXL north survey and the SDSS Baryon Oscillation Spectroscopic Survey (BOSS) spectroscopic follow-up of them, we compare the properties of X-ray unobscured and obscured broad-line AGN (BLAGN1 and BLAGN2; $N_\textrm{H}$below and above $10^{21.5}$ cm$^{-2}$), including their X-ray luminosity $L_X$, black hole mass, Eddington ratio $\lambda_{\textrm{Edd}}$, optical continuum and line features. We find that BLAGN2 have systematically larger broad line widths and hence apparently higher (lower) $M_{\textrm{BH}}$ ($\lambda_{\textrm{Edd}}$) than BLAGN1. We also find that the X-ray obscuration in BLAGN tends to coincide with optical dust extinction, which is optically thinner than that in narrow-line AGN (NLAGN) and likely partial-covering to the broad line region. All the results can be explained in the framework of a multi-component, clumpy torus model by interpreting BLAGN2 as an intermediate type between BLAGN1 and NLAGN in terms of an intermediate inclination angle.

Studies of galaxy clusters provide stringent constraints on models of structure formation. Provided that selection effects are under control, large X-ray surveys are well suited to derive cosmological parameters, in particular those governing the dark energy equation of state. We forecast the capabilities of the all-sky eROSITA (the extended ROentgen Survey with an Imaging Telescope Array) survey to be achieved by the early 2020s. We bring special attention to modeling the entire chain from photon emission to source detection and cataloguing. The selection function of galaxy clusters for the upcoming eROSITA mission is investigated by means of extensive and dedicated Monte-Carlo simulations. Employing a combination of accurate instrument characterization and of state-of-the-art source detection technique, we determine a cluster detection efficiency based on the cluster fluxes and sizes. Using this eROSITA cluster selection function, we find that eROSITA will detect a total of $\sim 10^5$ clusters in the extra-galactic sky. This number of clusters will allow eROSITA to put stringent constraints on cosmological models. We show that incomplete assumptions on selection effects, such as neglecting the distribution of cluster sizes, induce a bias in the derived value of cosmological parameters. Synthetic simulations of the eROSITA sky capture the essential characteristics impacting the next-generation galaxy cluster surveys and they highlight parameters requiring tight monitoring in order to avoid biases in cosmological analyses.

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 a sample of 329 low to intermediate redshift ($0.05 < z < 0.3$) brightest cluster galaxies (BCGs) in X-ray selected clusters from the SPectroscopic IDentification of eRosita Sources (SPIDERS) survey, a spectroscopic survey within Sloan Digital Sky Survey-IV (SDSS-IV). We define our BCGs by simultaneous consideration of legacy X-ray data from ROSAT, maximum likelihood outputs from an optical cluster-finder algorithm and visual inspection. Using SDSS imaging data, we fit S\'ersic profiles to our BCGs in three bands (\textit{g}, \textit{r}, \textit{i}) with \textsc{SIGMA}, a \textsc{GALFIT}-based software wrapper. We examine the reliability of our fits by running our pipeline on ${\sim}10^{4}$ psf-convolved model profiles injected into 8 random cluster fields, we then use the results of this analysis to create a robust subsample of 198 BCGs. We outline three cluster properties of interest: overall cluster X-ray luminosity ($L_{X}$), cluster richness as estimated by \textsc{redMaPPer} ($ \lambda $) and cluster halo mass ($M_{200}$), which is estimated via velocity dispersion. In general, there are significant correlations with BCG stellar mass between all three environmental properties, but no significant trends arise with either S\'ersic index or effective radius. There is no major environmental dependence on the strength of the relation between effective radius and BCG stellar mass. Stellar mass therefore arises as the most important factor governing BCG morphology. Our results indicate that our sample consists of a large number of relaxed, mature clusters containing broadly homogeneous BCGs up to $z \sim 0.3$, suggesting that there is little evidence for much ongoing structural evolution for BCGs in these systems.

We conduct a comprehensive study of the effects of incorporating galaxy morphology information in photometric redshift estimation. Using machine learning methods, we assess the changes in the scatter and catastrophic outlier fraction of photometric redshifts when galaxy size, ellipticity, S\'{e}rsic index and surface brightness are included in training on galaxy samples from the SDSS and the CFHT Stripe-82 Survey (CS82). We show that by adding galaxy morphological parameters to full $ugriz$ photometry, only mild improvements are obtained, while the gains are substantial in cases where fewer passbands are available. For instance, the combination of $grz$ photometry and morphological parameters almost fully recovers the metrics of $5$-band photometric redshifts. We demonstrate that with morphology it is possible to determine useful redshift distribution $N(z)$ of galaxy samples without any colour information. We also find that the inclusion of quasar redshifts and associated object sizes in training improves the quality of photometric redshift catalogues, compensating for the lack of a good star-galaxy separator. We further show that morphological information can mitigate biases and scatter due to bad photometry. As an application, we derive both point estimates and posterior distributions of redshifts for the official CS82 catalogue, training on morphology and SDSS Stripe-82 $ugriz$ bands when available. Our redshifts yield a 68th percentile error of $0.058(1+z)$, and a catastrophic outlier fraction of $5.2$ per cent. We further include a deep extension trained on morphology and single $i$-band CS82 photometry.

We present the clustering measurements of quasars in configuration space based on the Data Release 14 (DR14) of the Sloan Digital Sky Survey IV extended Baryon Oscillation Spectroscopic Survey. This dataset includes 148,659 quasars spread over the redshift range $0.8\leq z \leq 2.2$ and spanning 2112.9 square degrees. We use the Convolution Lagrangian Perturbation Theory (CLPT) approach with a Gaussian Streaming (GS) model for the redshift space distortions of the correlation function and demonstrate its applicability for dark matter halos hosting eBOSS quasar tracers. At the effective redshift $z_{\rm eff} = 1.52$, we measure the linear growth rate of structure $f\sigma_{8}(z_{\rm eff})= 0.426 \pm 0.077$, the expansion rate $H(z_{\rm eff})= 159^{+12}_{-13}(r_{s}^{\rm fid}/r_s){\rm km.s}^{-1}.{\rm Mpc}^{-1}$, and the angular diameter distance $D_{A}(z_{\rm eff})=1850^{+90}_{-115}\,(r_s/r_{s}^{\rm fid}){\rm Mpc}$, where $r_{s}$ is the sound horizon at the end of the baryon drag epoch and $r_{s}^{\rm fid}$ is its value in the fiducial cosmology. The quoted errors include both systematic and statistical contributions. The results on the evolution of distances are consistent with the predictions of flat $\Lambda$-Cold Dark Matter ($\Lambda$-CDM) cosmology with Planck parameters, and the measurement of $f\sigma_{8}$ extends the validity of General Relativity (GR) to higher redshifts($z>1$) This paper is released with companion papers using the same sample. The results on the cosmological parameters of the studies are found to be in very good agreement, providing clear evidence of the complementarity and of the robustness of the first full-shape clustering measurements with the eBOSS DR14 quasar sample.

The extended Baryon Oscillation Spectroscopic Survey (eBOSS) Data Release 14 sample includes 80,118 Luminous Red Galaxies. By combining these galaxies with the high-redshift tail of the BOSS galaxy sample, we form a sample of LRGs at an effective redshift $z=0.72$, covering an effective volume of 0.9~Gpc$^3$. We introduce new techniques to account for spurious fluctuations caused by targeting and by redshift failures which were validated on a set of mock catalogs. This analysis is sufficient to provide a $2.6$\% measurement of spherically averaged BAO, $D_V(z=0.72) = 2353^{+63}_{-61} (r_d/r_{d,\rm{fid}}) h^{-1}$Mpc, at 2.8$\sigma$ of significance. Together with the recent quasar-based BAO measurement at $z=1.5$, and forthcoming Emission Line Galaxy-based measurements, this measurement demonstrates that eBOSS is fulfilling its remit of extending the range of redshifts covered by such measurements, laying the ground work for forthcoming surveys such as the Dark Energy Spectroscopic Survey and Euclid.

Numerical simulations play a key important role in modern cosmology. Examples are plenty including the cosmic web - large scale structure of the distribution of galaxies in space - which was first observed in N-body simulations and later discovered in observations. The cuspy dark matter halo profiles, the overabundance of satellites, the Too-Big-Too-Fail problem are other examples of theoretical predictions that have a dramatic impact on recent developments in cosmology and galaxy formation. Large observational surveys such as e.g. SDSS, Euclid, and LSST are intimately connected with extensive cosmological simulations that provide statistical errors and tests for systematics. Accurate predictions for baryonic acoustic oscillations and redshift space distortions from high-resolution and large-volume cosmological simulations are required for interpretation of these large-scale galaxy/qso surveys. However, most of the results from extensive computer simulations, that would be greatly beneficial if publicly available, are still in hands of few research groups. Even when the simulation data is available, sharing vast amounts of data can be overwhelming. We argue that there is an effective and simple path to expand the data access and dissemination of numerous results from different cosmological models. Here we demonstrate that public access can be effectively provided with relatively modest resources. Among different results, we release for the astronomical community terabytes of raw data of th popular Bolshoi and MultiDark simulations. We also provide numerous results that are focused on mimicking observational data and galaxy surveys for major projects. Skies and Universes is a community effort: data are produced and shared by many research groups. We offer to other cosmologists and astronomers to host their data products in the skiesanduniverses.org space.

We investigate the mass-metallicity relations for the gaseous (MZRgas) and stellar components (MZRstar) of local star-forming galaxies based on a representative sample from SDSS DR12. The mass-weighted average stellar metallicities are systematically lower than the gas metallicities. This difference in metallicity increases toward galaxies with lower masses and reaches 0.4-0.8 dex at 10^9 Msun (depending on the gas metallicity calibration). As a result, the MZRstar is much steeper than the MZRgas. The much lower metallicities in stars compared to the gas in low mass galaxies implies dramatic metallicity evolution with suppressed metal enrichment at early times. The aim of this paper is to explain the observed large difference in gas and stellar metallicity and to infer the origin of the mass-metallicity relations. To this end we develop a galactic chemical evolution model accounting for star formation, gas inflow and outflow. By combining the observed mass-metallicity relation for both gas and stellar components to constrain the models, we find that only two scenarios are able to reproduce the observations. Either strong metal outflow or a steep IMF slope at early epochs of galaxy evolution is needed. Based on these two scenarios, for the first time we successfully reproduce the observed MZRgas and MZRstar simultaneously, together with other independent observational constraints in the local universe. Our model also naturally reproduces the flattening of the MZRgas at the high mass end leaving the MZRstar intact, as seen in observational data.

Emission line galaxies (ELGs) are used in several ongoing and upcoming surveys (SDSS-IV/eBOSS, DESI) as tracers of the dark matter distribution. Using a new galaxy formation model, we explore the characteristics of [OII] emitters, which dominate optical ELG selections at $z\simeq 1$. Model [OII] emitters at $0.5<z<1.5$ are selected to mimic the DEEP2, VVDS, eBOSS and DESI surveys. The luminosity functions of model [OII] emitters are in reasonable agreement with observations. The selected [OII] emitters are hosted by haloes with $M_{\rm halo}\geq 10^{10.3}h^{-1}{\rm M}_{\odot}$, with ~90% of them being central star-forming galaxies. The predicted mean halo occupation distributions of [OII] emitters has a shape typical of that inferred for star-forming galaxies, with the contribution from central galaxies, $\langle N \rangle_{\left[OII\right]\, cen}$, being far from the canonical step function. The $\langle N \rangle_{\left[OII\right]\, cen}$ can be described as the sum of an asymmetric Gaussian for disks and a step function for spheroids, which plateaus below unity. The model [OII] emitters have a clustering bias close to unity, which is below the expectations for eBOSS and DESI ELGs. At $z\sim 1$, a comparison with observed g-band selected galaxy, which are expected to be dominated by [OII] emitters, indicates that our model produces too few [OII] emitters that are satellite galaxies. This suggests the need to revise our modelling of hot gas stripping in satellite galaxies.

We present measurements of the Baryon Acoustic Oscillation (BAO) scale in redshift-space using the clustering of quasars. We consider a sample of 147,000 quasars from the extended Baryon Oscillation Spectroscopic Survey (eBOSS) distributed over 2044 square degrees with redshifts $0.8 < z < 2.2$ and measure their spherically-averaged clustering in both configuration and Fourier space. Our observational dataset and the 1400 simulated realizations of the dataset allow us to detect a preference for BAO that is greater than 2.8$\sigma$. We determine the spherically averaged BAO distance to $z = 1.52$ to 3.8 per cent precision: $D_V(z=1.52)=3843\pm147 \left(r_{\rm d}/r_{\rm d, fid}\right)\ $Mpc. This is the first time the location of the BAO feature has been measured between redshifts 1 and 2. Our result is fully consistent with the prediction obtained by extrapolating the Planck flat $\Lambda$CDM best-fit cosmology. All of our results are consistent with basic large-scale structure (LSS) theory, confirming quasars to be a reliable tracer of LSS, and provide a starting point for numerous cosmological tests to be performed with eBOSS quasar samples. We combine our result with previous, independent, BAO distance measurements to construct an updated BAO distance-ladder. Using these BAO data alone and marginalizing over the length of the standard ruler, we find $\Omega_{\Lambda} > 0$ at 6.6$\sigma$ significance when testing a $\Lambda$CDM model with free curvature.

We present the first scientific results from the luminous red galaxy sample (LRG) of the extended Baryon Oscillation Spectroscopic Survey (eBOSS). We measure the small and intermediate scale clustering from a sample of more than 61,000 galaxies in the redshift range $0.6 < z < 0.9$. We interpret these measurements in the framework of the Halo Occupation Distribution. The bias of eBOSS LRGs is $2.30 \pm 0.03$, with a satellite fraction of $13\pm3$\% and a mean halo mass of $2.5\times10^{13}h^{-1}M_{\odot}$. These results are consistent with expectations, demonstrating that eBOSS galaxies will be reliable tracers of large scale structure at $z\sim 0.7$. The eBOSS galaxy bias implies a scatter of luminosity at fixed halo mass, $\sigma_{\log L}$, of 0.19 dex. Using the clustering of massive galaxies from BOSS-CMASS, BOSS-LOWZ, and SDSS, we find that $\sigma_{\log L}=0.19$ is consistent with observations over the full redshift range that these samples cover. The addition of eBOSS to previous surveys allows investigation of the evolution of massive galaxies over the past $\sim 7$ Gyr.

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 present a new measurement of the mass-concentration relation and the stellar-to-halo mass ratio over the halo mass range $5\times 10^{12}$ to $2\times 10^{14}M_{\odot}$. To achieve this, we use weak lensing measurements from the CFHT Stripe 82 Survey (CS82), combined with the central galaxies from the redMaPPer cluster catalogue and the LOWZ/CMASS galaxy sample of the Sloan Digital Sky Survey-III Baryon Oscillation Spectroscopic Survey Tenth Data Release. The stacked lensing signals around these samples are modelled as a sum of contributions from the central galaxy, its dark matter halo, and the neighboring halos, as well as a term for possible centering errors. We measure the mass-concentration relation: $c_{200c}(M)=A(\frac{M_{200c}}{M_0})^{B}$ with $A=5.24\pm1.24, B=-0.13\pm0.10$ for $0.2<z<0.4$ and $A=6.61\pm0.75, B=-0.15\pm0.05$ for $0.4<z<0.6$. These amplitudes and slopes are completely consistent with predictions from recent simulations. We also measure the stellar-to-halo mass ratio for our samples, and find results consistent with previous measurements from lensing and other techniques.

We study the galaxy clustering dependence on the $\left[\mathrm{O\scriptsize{II}}\right]$ emission line luminosity in the SDSS DR7 Main galaxy sample at mean redshift $z\sim0.1$. We select volume-limited samples of galaxies with different $\left[\mathrm{O\scriptsize{II}}\right]$ luminosity thresholds and measure their projected, monopole and quadrupole two-point correlation functions. We model these observations using the 1$h^{-1}\rm{Gpc}$ MultiDark Planck cosmological simulation and generate light-cones with the SUrvey GenerAtoR algorithm. To interpret our results, we adopt a modified (Sub)Halo Abundance Matching scheme, accounting for the stellar mass incompleteness of the emission line galaxies. The satellite fraction constitutes an extra parameter in this model and allows to optimize the clustering fit on both small and intermediate scales (i.e. $r_p\lesssim 30h^{-1}\rm{Mpc})$, with no need of any velocity bias correction. We find that, in the local Universe, the $\left[\mathrm{O\scriptsize{II}}\right]$ luminosity correlates with all the clustering statistics explored and with the galaxy bias. This latter quantity correlates more strongly with the SDSS $r$-band magnitude than $\left[\mathrm{O\scriptsize{II}}\right]$ luminosity. In conclusion, we propose a straightforward method to produce reliable clustering models, entirely built on the simulation products, which provides robust predictions of the typical ELG host halo masses and satellite fraction values. The SDSS galaxy data, MultiDark mock catalogues and clustering results are made publicly available.

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.

$N$-body cosmological simulations are an essential tool to understand the observed distribution of galaxies. We use the MultiDark simulation suite, run with the Planck cosmological parameters, to revisit the mass and velocity functions. At redshift $z=0$, the simulations cover four orders of magnitude in halo mass from $\sim10^{11}M_\odot$ with 8,783,874 distinct halos and 532,533 subhalos. The total volume used is $\sim$515 Gpc$^3$, more than 8 times larger than in previous studies. We measure and model the halo mass function, its covariance matrix w.r.t halo mass and the large scale halo bias. With the formalism of the excursion-set mass function, we explicit the tight interconnection between the covariance matrix, bias and halo mass function. We obtain a very accurate ($<2\%$ level) model of the distinct halo mass function. We also model the subhalo mass function and its relation to the distinct halo mass function. The set of models obtained provides a complete and precise framework for the description of halos in the concordance Planck cosmology. Finally, we provide precise analytical fits of the $V_{max}$ maximum velocity function up to redshift $z<2.3$ to push for the development of halo occupation distribution using $V_{max}$. The data and the analysis code are made publicly available in the \textit{Skies and Universes} database.

We present the VISTA-CFHT Stripe 82 (VICS82) survey: a near-infrared (J+Ks) survey covering 150 square degrees of the Sloan Digital Sky Survey (SDSS) equatorial Stripe 82 to an average depth of J=21.9 AB mag and Ks=21.4 AB mag (80% completeness limits; 5-sigma point source depths are approximately 0.5 mag brighter). VICS82 contributes to the growing legacy of multi-wavelength data in the Stripe 82 footprint. The addition of near-infrared photometry to the existing SDSS Stripe 82 coadd ugriz photometry reduces the scatter in stellar mass estimates to delta log(M_stellar)~0.3 dex for galaxies with M_stellar>10^9M_sun at z~0.5, and offers improvement compared to optical-only estimates out to z~1, with stellar masses constrained within a factor of approximately 2.5. When combined with other multi-wavelength imaging of the Stripe, including moderate-to-deep ultraviolet (GALEX), optical and mid-infrared (Spitzer IRAC) coverage, as well as tens of thousands of spectroscopic redshifts, VICS82 gives access to approximately 0.5 Gpc^3 of comoving volume. Some of the main science drivers of VICS82 include (a) measuring the stellar mass function of L^star galaxies out to z~1; (b) detecting intermediate redshift quasars at 2<z<3.5; (c) measuring the stellar mass function and baryon census of clusters of galaxies, and (d) performing optical/near-infrared-cosmic microwave background lensing cross-correlation experiments linking stellar mass to large-scale dark matter structure. Here we define and describe the survey, highlight some early science results and present the first public data release, which includes an SDSS-matched catalogue as well as the calibrated pixel data itself.

We study the first year of the eBOSS quasar sample in the redshift range $0.9<z<2.2$ which includes 68,772 homogeneously selected quasars. We show that the main source of systematics in the evaluation of the correlation function arises from inhomogeneities in the quasar target selection, particularly related to the extinction and depth of the imaging data used for targeting. We propose a weighting scheme that mitigates these systematics. We measure the quasar correlation function and provide the most accurate measurement to date of the quasar bias in this redshift range, $b_Q = 2.45 \pm 0.05$ at $\bar z=1.55$, together with its evolution with redshift. We use this information to determine the minimum mass of the halo hosting the quasars and the characteristic halo mass, which we find to be both independent of redshift within statistical error. Using a recently-measured quasar-luminosity-function we also determine the quasar duty cycle. The size of this first year sample is insufficient to detect any luminosity dependence to quasar clustering and this issue should be further studied with the final $\sim$500,000 eBOSS quasar sample.

We describe the algorithm used to select the Emission Line Galaxy (ELG) sample at $z \sim 0.85$ for the extended Baryon Oscillation Spectroscopic Survey of the Sloan Digital Sky Survey IV, using photometric data from the DECam Legacy Survey. Our selection is based on a selection box in the $g-r$ vs. $r-z$ colour-colour space and a cut on the $g$-band magnitude, to favour galaxies in the desired redshift range with strong [OII] emission. It provides a target density of 200 deg$^{-2}$ on the North Galactic Cap (NGC) and of 240 deg$^{-2}$ on the South Galactic Cap (SGC), where we use a larger selection box because of deeper imaging. We demonstrate that this selection passes the eBOSS requirements in terms of homogeneity. About 50,000 ELGs have been observed since the observations have started in 2016, September. These roughly match the expected redshift distribution, though the measured efficiency is slightly lower than expected. The efficiency can be increased by enlarging the redshift range and with incoming pipeline improvement. The cosmological forecast based on these first data predict $\sigma_{D_V}/D_V = 0.023$, in agreement with previous forecasts. Lastly, we present the stellar population properties of the ELG SGC sample. Once observations are completed, this sample will be suited to provide a cosmological analysis at $z \sim 0.85$, and will pave the way for the next decade of massive spectroscopic cosmological surveys, which heavily rely on ELGs. The target catalogue over the SGC will be released along with DR14.

In current and future surveys, quasars play a key role. The new data will extend our knowledge of the Universe as it will be used to better constrain the cosmological model at redshift $z>1$ via baryon acoustic oscillation and redshift space distortion measurements. Here, we present the first clustering study of quasars observed by the extended Baryon Oscillation Spectroscopic Survey. We measure the clustering of $\sim 70,000$ quasars located in the redshift range $0.9<z<2.2$ that cover 1,168 deg$^2$. We model the clustering and produce high-fidelity quasar mock catalogues based on the BigMultiDark Planck simulation. Thus, we use a modified (Sub)Halo Abundance Matching model to account for the specificities of the halo population hosting quasars. We find that quasars are hosted by halos with masses $\sim10^{12.7}M_\odot$ and their bias evolves from 1.54 ($z=1.06$) to 3.15 ($z=1.98$). Using the current eBOSS data, we cannot distinguish between models with different fractions of satellites. The high-fidelity mock light-cones, including properties of halos hosting quasars, are made publicly available.

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$\alpha$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.

DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.

We present two wide-field catalogs of photometrically-selected emission line galaxies (ELGs) at z=0.8 covering about 2800 deg^2 over the south galactic cap. The catalogs were obtained using a Fisher discriminant technique described in a companion paper. The two catalogs differ by the imaging used to define the Fisher discriminant: the first catalog includes imaging from the Sloan Digital Sky Survey and the Wide-Field Infrared Survey Explorer, the second also includes information from the South Galactic Cap U-band Sky Survey (SCUSS). Containing respectively 560,045 and 615,601 objects, they represent the largest ELG catalogs available today and were designed for the ELG programme of the extended Baryon Oscillation Spectroscopic Survey (eBOSS). We study potential sources of systematic variation in the angular distribution of the selected ELGs due to fluctuations of the observational parameters. We model the influence of the observational parameters using a multivariate regression and implement a weighting scheme that allows effective removal of all of the systematic errors induced by the observational parameters. We show that fluctuations in the imaging zero-points of the photometric bands have minor impact on the angular distribution of objects in our catalogs. We compute the angular clustering of both catalogs and show that our weighting procedure effectively removes spurious clustering on large scales. We fit a model to the small scale angular clustering, showing that the selections have similar biases of 1.35/D_a(z) and 1.28/D_a(z). Both catalogs are publicly available.

We present a direct measurement of the mean halo occupation distribution (HOD) of galaxies taken from the eleventh data release (DR11) of the Sloan Digital Sky Survey-III Baryon Oscillation Spectroscopic Survey (BOSS). The HOD of BOSS low-redshift (LOWZ: $0.2 < z < 0.4$) and Constant-Mass (CMASS: $0.43 <z <0.7$) galaxies is inferred via their association with the dark-matter halos of 174 X-ray-selected galaxy clusters drawn from the XMM Cluster Survey (XCS). Halo masses are determined for each galaxy cluster based on X-ray temperature measurements, and range between ${\rm log_{10}} (M_{180}/M_{\odot}) = 13-15$. Our directly measured HODs are consistent with the HOD-model fits inferred via the galaxy-clustering analyses of Parejko et al. for the BOSS LOWZ sample and White et al. for the BOSS CMASS sample. Under the simplifying assumption that the other parameters that describe the HOD hold the values measured by these authors, we have determined a best-fit alpha-index of 0.91$\pm$0.08 and $1.27^{+0.03}_{-0.04}$ for the CMASS and LOWZ HOD, respectively. These alpha-index values are consistent with those measured by White et al. and Parejko et al. In summary, our study provides independent support for the HOD models assumed during the development of the BOSS mock-galaxy catalogues that have subsequently been used to derive BOSS cosmological constraints.

The Astropy Project (http://astropy.org) is, in its own words, "a community effort to develop a single core package for Astronomy in Python and foster interoperability between Python astronomy packages." For five years this project has been managed, written, and operated as a grassroots, self-organized, almost entirely volunteer effort while the software is used by the majority of the astronomical community. Despite this, the project has always been and remains to this day effectively unfunded. Further, contributors receive little or no formal recognition for creating and supporting what is now critical software. This paper explores the problem in detail, outlines possible solutions to correct this, and presents a few suggestions on how to address the sustainability of general purpose astronomical software.

We present the Data Release 12 Quasar catalog (DR12Q) from the Baryon Oscillation Spectroscopic Survey (BOSS) of the SDSS-III. This catalog includes all SDSS-III/BOSS objects that were spectroscopically targeted as quasar candidates during the full survey and that are confirmed as quasars via visual inspection of the spectra, have luminosities Mi[z=2]<-20.5 (in a LCDM cosmology with H_0 = 70 km/s/Mpc, Omega_M = 0.3, and Omega _L=0.7), and either display at least one emission line with a full width at half maximum (FWHM)larger than 500 km/s or, if not, have interesting/complex absorption features. The catalog also includes previously known quasars (mostly from SDSS-I and II) that were reobserved by BOSS. The catalog contains 297,301 quasars detected over 9,376 square degrees with robust identification and redshift measured by a combination of principal component eigenspectra. The number of quasars with z>2.15 is about an order of magnitude greater than the number of z>2.15 quasars known prior to BOSS. Redshifts and FWHMs are provided for the strongest emission lines (CIV, CIII], MgII). The catalog identifies 29,580 broad absorption line quasars and lists their characteristics. For each object, the catalog presents five-band (u, g, r, i, z) CCD-based photometry together with some information on the optical morphology and the selection criteria. When available, the catalog also provides information on the optical variability of quasars using SDSS and PTF multi-epoch photometry. The catalog also contains X-ray, ultraviolet, near-infrared, and radio emission properties of the quasars, when available, from other large-area surveys. The calibrated digital spectra, covering the wavelength region 3,600-10,500A at a spectral resolution in the range 1,300<R<2,500, can be retrieved from the SDSS Catalog Archive Server.

We present the measurements of the luminosity-dependent redshift-space three-point correlation functions (3PCFs) for the Sloan Digital Sky Survey (SDSS) DR7 Main galaxy sample. We compare the 3PCF measurements to the predictions from three different halo and subhalo models. One is the halo occupation distribution (HOD) model and the other two are extensions of the subhalo abundance matching (SHAM) model by allowing the central and satellite galaxies to have different occupation distributions in the host halos and subhalos. Parameters in all the models are chosen to best describe the projected and redshift-space two-point correlation functions (2PCFs) of the same set of galaxies. All three model predictions agree well with the 3PCF measurements for the most luminous galaxy sample, while the HOD model better performs in matching the 3PCFs of fainter samples (with luminosity threshold below $L^*$), which is similar in trend to the case of fitting the 2PCFs. The decomposition of the model 3PCFs into contributions from different types of galaxy triplets shows that on small scales the dependence of the 3PCFs on triangle shape is driven by nonlinear redshift-space distortion (and not by the intrinsic halo shape) while on large scales it reflects the filamentary structure. The decomposition also reveals more detailed differences in the three models, which are related to the radial distribution, the mean occupation function, and the velocity distribution of satellite galaxies inside halos. The results suggest that galaxy 3PCFs can further help constrain the above galaxy-halo relation and test theoretical models.

Current and future large redshift surveys, as the Sloan Digital Sky Survey IV extended Baryon Oscillation Spectroscopic Survey (SDSS-IV/eBOSS) or the Dark Energy Spectroscopic Instrument (DESI), will use emission-line galaxies (ELG) to probe cosmological models by mapping the large-scale structure of the Universe in the redshift range $0.6 < z < 1.7$. With current data, we explore the halo-galaxy connection by measuring three clustering properties of $g$-selected ELGs as matter tracers in the redshift range $0.6 < z < 1$: (i) the redshift-space two-point correlation function using spectroscopic redshifts from the BOSS ELG sample and VIPERS; (ii) the angular two-point correlation function on the footprint of the CFHT-LS; (iii) the galaxy-galaxy lensing signal around the ELGs using the CFHTLenS. We interpret these observations by mapping them onto the latest high-resolution MultiDark Planck N-body simulation, using a novel (Sub)Halo-Abundance Matching technique that accounts for the ELG incompleteness. ELGs at $z\sim0.8$ live in halos of $(1\pm 0.5)\times10^{12}\,h^{-1}$M$_{\odot}$ and 22.5$\pm2.5$% of them are satellites belonging to a larger halo. The halo occupation distribution of ELGs indicates that we are sampling the galaxies in which stars form in the most efficient way, according to their stellar-to-halo mass ratio.

We present cosmological results from the final galaxy clustering data set of the Baryon Oscillation Spectroscopic Survey, part of the Sloan Digital Sky Survey III. Our combined galaxy sample comprises 1.2 million massive galaxies over an effective area of 9329 deg^2 and volume of 18.7 Gpc^3, divided into three partially overlapping redshift slices centred at effective redshifts 0.38, 0.51, and 0.61. We measure the angular diameter distance DM and Hubble parameter H from the baryon acoustic oscillation (BAO) method after applying reconstruction to reduce non-linear effects on the BAO feature. Using the anisotropic clustering of the pre-reconstruction density field, we measure the product DM*H from the Alcock-Paczynski (AP) effect and the growth of structure, quantified by f{\sigma}8(z), from redshift-space distortions (RSD). We combine measurements presented in seven companion papers into a set of consensus values and likelihoods, obtaining constraints that are tighter and more robust than those from any one method. Combined with Planck 2015 cosmic microwave background measurements, our distance scale measurements simultaneously imply curvature {\Omega}_K =0.0003+/-0.0026 and a dark energy equation of state parameter w = -1.01+/-0.06, in strong affirmation of the spatially flat cold dark matter model with a cosmological constant ({\Lambda}CDM). Our RSD measurements of f{\sigma}_8, at 6 per cent precision, are similarly consistent with this model. When combined with supernova Ia data, we find H0 = 67.3+/-1.0 km/s/Mpc even for our most general dark energy model, in tension with some direct measurements. Adding extra relativistic species as a degree of freedom loosens the constraint only slightly, to H0 = 67.8+/-1.2 km/s/Mpc. Assuming flat {\Lambda}CDM we find {\Omega}_m = 0.310+/-0.005 and H0 = 67.6+/-0.5 km/s/Mpc, and we find a 95% upper limit of 0.16 eV/c^2 on the neutrino mass sum.

Emission line galaxies are one of the main tracers of the large-scale structure to be targeted by the next-generation dark energy surveys. To provide a better understanding of the properties and statistics of these galaxies, we have collected spectroscopic data from the VVDS and DEEP2 deep surveys and estimated the galaxy luminosity functions (LFs) of three distinct emission lines, [OII] ($0.5 < z < 1.3$), H{\beta} ($0.3 < z < 0.8$) and [OIII] ($0.3 < z < 0.8$). Our measurements are based on 35,639 emission line galaxies and cover a volume of $\sim10^7$Mpc$^3$. We present the first measurement of the H{\beta} LF at these redshifts. We have also compiled LFs from the literature that were based on independent data or covered different redshift ranges, and we fit the entire set over the whole redshift range with analytic Schechter and Saunders models, assuming a natural redshift dependence of the parameters. We find that the characteristic luminosity ($L_*$) and density ($\phi_*$) of all LFs increase with redshift. Using the Schechter model over the redshift ranges considered, we find that, for [OII] emitters, the characteristic luminosity $L_*(z=0.5)=3.2\times10^{41}$ erg/s increases by a factor of $2.7 \pm 0.2$ from z=0.5 to 1.3; for H{\beta} emitters $L_*(z=0.3)=1.3\times10^{41}$ erg/s increases by a factor of $2.0 \pm 0.2$ from z=0.3 to 0.8; and for [OIII] emitters $L_*(z=0.3)=7.3\times10^{41}$ erg/s increases by a factor of $3.5 \pm 0.4$ from z=0.3 to 0.8.

The Sloan Digital Sky Survey IV extended Baryonic Oscillation Spectroscopic Survey (SDSS-IV/eBOSS) will observe 195,000 emission-line galaxies (ELGs) to measure the Baryonic Acoustic Oscillation standard ruler (BAO) at redshift 0.9. To test different ELG selection algorithms, 9,000 spectra were observed with the SDSS spectrograph as a pilot survey based on data from several imaging surveys. First, using visual inspection and redshift quality flags, we show that the automated spectroscopic redshifts assigned by the pipeline meet the quality requirements for a reliable BAO measurement. We also show the correlations between sky emission, signal-to-noise ratio in the emission lines, and redshift error. Then we provide a detailed description of each target selection algorithm we tested and compare them with the requirements of the eBOSS experiment. As a result, we provide reliable redshift distributions for the different target selection schemes we tested. Finally, we determine an target selection algorithms that is best suited to be applied on DECam photometry because they fulfill the eBOSS survey efficiency requirements.

MSE is an 11.25m aperture observatory with a 1.5 square degree field of view that will be fully dedicated to multi-object spectroscopy. More than 3200 fibres will feed spectrographs operating at low (R ~ 2000 - 3500) and moderate (R ~ 6000) spectral resolution, and approximately 1000 fibers will feed spectrographs operating at high (R ~ 40000) resolution. MSE is designed to enable transformational science in areas as diverse as tomographic mapping of the interstellar and intergalactic media; the in-situ chemical tagging of thick disk and halo stars; connecting galaxies to their large scale structure; measuring the mass functions of cold dark matter sub-halos in galaxy and cluster-scale hosts; reverberation mapping of supermassive black holes in quasars; next generation cosmological surveys using redshift space distortions and peculiar velocities. MSE is an essential follow-up facility to current and next generations of multi-wavelength imaging surveys, including LSST, Gaia, Euclid, WFIRST, PLATO, and the SKA, and is designed to complement and go beyond the science goals of other planned and current spectroscopic capabilities like VISTA/4MOST, WHT/WEAVE, AAT/HERMES and Subaru/PFS. It is an ideal feeder facility for E-ELT, TMT and GMT, and provides the missing link between wide field imaging and small field precision astronomy. MSE is optimized for high throughput, high signal-to-noise observations of the faintest sources in the Universe with high quality calibration and stability being ensured through the dedicated operational mode of the observatory. (abridged)

We have developed an analytical method based on forward-modeling techniques to characterize the high-mass end of the red sequence (RS) galaxy population at redshift $z\sim0.55$, from the DR10 BOSS CMASS spectroscopic sample, which comprises $\sim600,000$ galaxies. The method, which follows an unbinned maximum likelihood approach, allows the deconvolution of the intrinsic CMASS colour-colour-magnitude distributions from photometric errors and selection effects. This procedure requires modeling the covariance matrix for the i-band magnitude, g-r colour and r-i colour using Stripe 82 multi-epoch data. Our results indicate that the error-deconvolved intrinsic RS distribution is consistent, within the photometric uncertainties, with a single point ($<0.05~{\rm{mag}}$) in the colour-colour plane at fixed magnitude, for a narrow redshift slice. We have computed the high-mass end ($^{0.55}M_i \lesssim -22$) of the $^{0.55}i$-band RS Luminosity Function (RS LF) in several redshift slices within the redshift range $0.52<z<0.63$. In this narrow redshift range, the evolution of the RS LF is consistent, within the uncertainties in the modeling, with a passively-evolving model with $\Phi_* = (7.248 \pm 0.204) \times10^{-4}$ Mpc$^{-3}$ mag$^{-1}$, fading at a rate of $1.5\pm0.4$ mag per unit redshift. We report RS completeness as a function of magnitude and redshift in the CMASS sample, which will facilitate a variety of galaxy-evolution and clustering studies using BOSS. Our forward-modeling method lays the foundations for future studies using other dark-energy surveys like eBOSS or DESI, which are affected by the same type of photometric blurring/selection effects.

We present a study of the clustering and halo occupation distribution of BOSS CMASS galaxies in the redshift range 0.43 < z < 0.7 drawn from the Final SDSS-III Data Release. We compare the BOSS results with the predictions of a Halo Abundance Matching (HAM) clustering model that assigns galaxies to dark matter halos selected from the large BigMultiDark $N$-body simulation of a flat $\Lambda$CDM Planck cosmology. We compare the observational data with the simulated ones on a light-cone constructed from 20 subsequent outputs of the simulation. Observational effects such as incompleteness, geometry, veto masks and fiber collisions are included in the model, which reproduces within 1-$\sigma$ errors the observed monopole of the 2-point correlation function at all relevant scales: from the smallest scales, 0.5 $h^{-1}$ Mpc, up to scales beyond the Baryonic Acoustic Oscillation feature. This model also agrees remarkably well with the BOSS galaxy power spectrum (up to $k\sim1$ $h$ Mpc$^{-1}$), and the Three-point correlation function. The quadrupole of the correlation function presents some tensions with observations. We discuss possible causes that can explain this disagreement, including target selection effects. Overall, the standard HAM model describes remarkably well the clustering statistics of the CMASS sample. We compare the stellar to halo mass relation for the CMASS sample measured using weak lensing in the CFHT Stripe 82 Survey with the prediction of our clustering model, and find a good agreement within 1-$\sigma$. The BigMD-BOSS light-cone including properties of BOSS galaxies and halo properties is made publicly available.

We model the luminosity-dependent projected and redshift-space two-point correlation functions (2PCFs) of the Sloan Digital Sky Survey (SDSS) DR7 Main galaxy sample, using the halo occupation distribution (HOD) model and the subhalo abundance matching (SHAM) model and its extension. All the models are built on the same high-resolution $N$-body simulations. We find that the HOD model generally provides the best performance in reproducing the clustering measurements in both projected and redshift spaces. The SHAM model with the same halo-galaxy relation for central and satellite galaxies (or distinct haloes and subhaloes), when including scatters, has a best-fitting $\chi^2/\rm{dof}$ around $2$--$3$. We therefore extend the SHAM model to the subhalo clustering and abundance matching (SCAM) by allowing the central and satellite galaxies to have different galaxy--halo relations. We infer the corresponding halo/subhalo parameters by jointly fitting the galaxy 2PCFs and abundances and consider subhaloes selected based on three properties, the mass $M_{\rm acc}$ at the time of accretion, the maximum circular velocity $V_{\rm acc}$ at the time of accretion, and the peak maximum circular velocity $V_{\rm peak}$ over the history of the subhaloes. The three subhalo models work well for luminous galaxy samples (with luminosity above $L_*$). For low-luminosity samples, the $V_{\rm acc}$ model stands out in reproducing the data, with the $V_{\rm peak}$ model slightly worse, while the $M_{\rm acc}$ model fails to fit the data. We discuss the implications of the modeling results.

In this paper we present a large database of weak lensing light cones constructed using different snapshots from the Big MultiDark simulation (BigMDPL). The ray-tracing through different multiple planes has been performed with the GLAMER code accounting both for single source redshifts and for sources distributed along the cosmic time. This first paper presents weak lensing forecasts and results according to the geometry of the VIPERS-W1 and VIPERS-W4 field of view. Additional fields will be available on our database and new ones can be run upon request. Our database also contains some tools for lensing analysis. In this paper we present results for convergence power spectra, one point and high order weak lensing statistics useful for forecasts and for cosmological studies. Covariance matrices have also been computed for the different realisations of the W1 and W4 fields. In addition we compute also galaxy-shear and projected density contrasts for different halo masses at two lens redshifts according to the CFHTLS source redshift distribution both using stacking and cross-correlation techniques, finding very good agreement.

The Baryon Acoustic Oscillation (BAO) feature in the power spectrum of galaxies provides a standard ruler to measure the accelerated expansion of the Universe. To extract all available information about dark energy, it is necessary to measure a standard ruler in the local, z<0.2, universe where dark energy dominates most the energy density of the Universe. Though the volume available in the local universe is limited, it is just big enough to measure accurately the long 100 Mpc/h wave-mode of the BAO. Using cosmological N-body simulations and approximate methods based on Lagrangian perturbation theory, we construct a suite of a thousand light-cones to evaluate the precision at which one can measure the BAO standard ruler in the local universe. We find that using the most massive galaxies on the full sky (34,000 sq. deg.), i.e. a K(2MASS)<14 magnitude-limited sample, one can measure the BAO scale up to a precision of 4\% and 1.2\% using reconstruction). We also find that such a survey would help to detect the dynamics of dark energy.Therefore, we propose a 3-year long observational project, named the Low Redshift survey at Calar Alto (LoRCA), to observe spectroscopically about 200,000 galaxies in the northern sky to contribute to the construction of aforementioned galaxy sample. The suite of light-cones is made available to the public.

We present a science forecast for the eBOSS survey, part of the SDSS-IV project, which is a spectroscopic survey using multiple tracers of large-scale structure, including luminous red galaxies (LRGs), emission line galaxies (ELGs) and quasars (both as a direct probe of structure and through the Ly-$\alpha$ forest). Focusing on discrete tracers, we forecast the expected accuracy of the baryonic acoustic oscillation (BAO), the redshift-space distortion (RSD) measurements, the $f_{\rm NL}$ parameter quantifying the primordial non-Gaussianity, the dark energy and modified gravity parameters. We also use the line-of-sight clustering in the Ly-$\alpha$ forest to constrain the total neutrino mass. We find that eBOSS LRGs ($0.6<z<1.0$) (combined with the BOSS LRGs at $z>0.6$), ELGs ($0.6<z<1.2$) and Clustering Quasars (CQs) ($0.6<z<2.2$) can achieve a precision of 1%, 2.2% and 1.6% precisions, respectively, for spherically averaged BAO distance measurements. Using the same samples, the constraint on $f\sigma_8$ is expected to be 2.5%, 3.3% and 2.8% respectively. For primordial non-Gaussianity, eBOSS alone can reach an accuracy of $\sigma(f_{\rm NL})\sim10-15$, depending on the external measurement of the galaxy bias and our ability to model large-scale systematic errors. eBOSS can at most improve the dark energy Figure of Merit (FoM) by a factor of $3$ for the Chevallier-Polarski-Linder (CPL) parametrisation, and can well constrain three eigenmodes for the general equation-of-state parameter (Abridged).

We reproduce the galaxy clustering catalogue from the SDSS-III Baryon Oscillation Spectroscopic Survey Final Data Release (BOSS DR11&DR12) with high fidelity on all relevant scales in order to allow a robust analysis of baryon acoustic oscillations and redshift space distortions. We have generated (6,000) 12,288 MultiDark PATCHY BOSS (DR11) DR12 light-cones corresponding to an effective volume of $\sim192,000\,[h^{-1}\,{\rm Gpc}]^3$ (the largest ever simulated volume), including cosmic evolution in the redshift range from 0.15 to 0.75. The mocks have been calibrated using a reference galaxy catalogue based on the halo abundance matching modelling of the BOSS DR11&DR12 galaxy clustering data and on the data themselves. The production follows three steps. First, we apply the PATCHY code to generate a dark matter field and an object distribution including nonlinear stochastic galaxy bias. Secondly, we run the halo/stellar distribution reconstruction HADRON code to assign masses to the various objects. This step uses the mass distribution as a function of local density and non-local indicators (i.e., tidal field tensor eigenvalues and relative halo exclusion separation for massive objects) from the reference simulation applied to the corresponding PATCHY dark matter and galaxy distribution. Finally, we apply the SUGAR code to build the light cones. The resulting MultiDark PATCHY mock light cones reproduce the number density, selection function, survey geometry, and in general within 1 $\sigma$, for arbitrary stellar mass bins, the power spectrum up to $k=0.3\,h\,{\rm Mpc}^{-1}$, the two-point correlation functions down to a few Mpc scales, and the three-point statistics of the BOSS DR11&DR12 galaxy samples.

The Extended Baryon Oscillation Spectroscopic Survey (eBOSS) will conduct novel cosmological observations using the BOSS spectrograph at Apache Point Observatory. Observations will be simultaneous with the Time Domain Spectroscopic Survey (TDSS) designed for variability studies and the Spectroscopic Identification of eROSITA Sources (SPIDERS) program designed for studies of X-ray sources. eBOSS will use four different tracers to measure the distance-redshift relation with baryon acoustic oscillations (BAO). Using more than 250,000 new, spectroscopically confirmed luminous red galaxies at a median redshift z=0.72, we project that eBOSS will yield measurements of $d_A(z)$ to an accuracy of 1.2% and measurements of H(z) to 2.1% when combined with the z>0.6 sample of BOSS galaxies. With ~195,000 new emission line galaxy redshifts, we expect BAO measurements of $d_A(z)$ to an accuracy of 3.1% and H(z) to 4.7% at an effective redshift of z= 0.87. A sample of more than 500,000 spectroscopically-confirmed quasars will provide the first BAO distance measurements over the redshift range 0.9<z<2.2, with expected precision of 2.8% and 4.2% on $d_A(z)$ and H(z), respectively. Finally, with 60,000 new quasars and re-observation of 60,000 quasars known from BOSS, we will obtain new Lyman-alpha forest measurements at redshifts z>2.1; these new data will enhance the precision of $d_A(z)$ and H(z) by a factor of 1.44 relative to BOSS. Furthermore, eBOSS will provide improved tests of General Relativity on cosmological scales through redshift-space distortion measurements, improved tests for non-Gaussianity in the primordial density field, and new constraints on the summed mass of all neutrino species. Here, we provide an overview of the cosmological goals, spectroscopic target sample, demonstration of spectral quality from early data, and projected cosmological constraints from eBOSS.

We present the results of the first test plates of the extended Baryon Oscillation Spectroscopic Survey. This paper focuses on the emission line galaxies (ELG) population targetted from the Dark Energy Survey (DES) photometry. We analyse the success rate, efficiency, redshift distribution, and clustering properties of the targets. From the 9000 spectroscopic redshifts targetted, 4600 have been selected from the DES photometry. The total success rate for redshifts between 0.6 and 1.2 is 71\% and 68\% respectively for a bright and faint, on average more distant, samples including redshifts measured from a single strong emission line. We find a mean redshift of 0.8 and 0.87, with 15 and 13\% of unknown redshifts respectively for the bright and faint samples. In the redshift range 0.6<z<1.2, for the most secure spectroscopic redshifts, the mean redshift for the bright and faint sample is 0.85 and 0.9 respectively. Star contamination is lower than 2\%. We measure a galaxy bias averaged on scales of 1 and 10~Mpc/h of 1.72 \pm 0.1 for the bright sample and of 1.78 \pm 0.12 for the faint sample. The error on the galaxy bias have been obtained propagating the errors in the correlation function to the fitted parameters. This redshift evolution for the galaxy bias is in agreement with theoretical expectations for a galaxy population with MB-5\log h < -21.0. We note that biasing is derived from the galaxy clustering relative to a model for the mass fluctuations. We investigate the quality of the DES photometric redshifts and find that the outlier fraction can be reduced using a comparison between template fitting and neural network, or using a random forest algorithm.

As part of the Sloan Digital Sky Survey IV the extended Baryon Oscillation Spectroscopic Survey (eBOSS) will improve measurements of the cosmological distance scale by applying the Baryon Acoustic Oscillation (BAO) method to quasar samples. eBOSS will adopt two approaches to target quasars over 7500 sq. deg. First, a "CORE" quasar sample will combine optical selection in ugriz using a likelihood-based routine called XDQSOz, with a mid-IR-optical color-cut. eBOSS CORE selection (to g < 22 OR r < 22) should return ~ 70 quasars per sq. deg. at redshifts 0.9 < z < 2.2 and ~7 z > 2.1 quasars per sq. deg. Second, a selection based on variability in multi-epoch imaging from the Palomar Transient Factory should recover an additional ~3-4 z > 2.1 quasars per sq. deg. to g < 22.5. A linear model of how imaging systematics affect target density recovers the angular distribution of eBOSS CORE quasars over 96.7% (76.7%) of the SDSS North (South) Galactic Cap area. The eBOSS CORE quasar sample should thus be sufficiently dense and homogeneous over 0.9 < z < 2.2 to yield the first few-percent-level BAO constraint near z~1.5. eBOSS quasars at z > 2.1 will be used to improve BAO measurements in the Lyman-alpha Forest. Beyond its key cosmological goals, eBOSS should be the next-generation quasar survey, comprising > 500,000 new quasars and > 500,000 uniformly selected spectroscopically confirmed 0.9 < z < 2.2 quasars. At the conclusion of eBOSS, the SDSS will have provided unique spectra of over 800,000 quasars.

We derive constraints on cosmological parameters and tests of dark energy models from the combination of baryon acoustic oscillation (BAO) measurements with cosmic microwave background (CMB) and Type Ia supernova (SN) data. We take advantage of high-precision BAO measurements from galaxy clustering and the Ly-alpha forest (LyaF) in the BOSS survey of SDSS-III. BAO data alone yield a high confidence detection of dark energy, and in combination with the CMB angular acoustic scale they further imply a nearly flat universe. Combining BAO and SN data into an "inverse distance ladder" yields a 1.7% measurement of $H_0=67.3 \pm1.1$ km/s/Mpc. This measurement assumes standard pre-recombination physics but is insensitive to assumptions about dark energy or space curvature, so agreement with CMB-based estimates that assume a flat LCDM cosmology is an important corroboration of this minimal cosmological model. For open LCDM, our BAO+SN+CMB combination yields $\Omega_m=0.301 \pm 0.008$ and curvature $\Omega_k=-0.003 \pm 0.003$. When we allow more general forms of evolving dark energy, the BAO+SN+CMB parameter constraints remain consistent with flat LCDM. While the overall $\chi^2$ of model fits is satisfactory, the LyaF BAO measurements are in moderate (2-2.5 sigma) tension with model predictions. Models with early dark energy that tracks the dominant energy component at high redshifts remain consistent with our constraints. Expansion history alone yields an upper limit of 0.56 eV on the summed mass of neutrino species, improving to 0.26 eV if we include Planck CMB lensing. Standard dark energy models constrained by our data predict a level of matter clustering that is high compared to most, but not all, observational estimates. (Abridged)

We describe the algorithm used to select the Luminous Red Galaxy (LRG) sample for the extended Baryon Oscillation Spectroscopic Survey (eBOSS) of the Sloan Digital Sky Survey IV (SDSS-IV) using photometric data from both the SDSS and the Wide-Field Infrared Survey Explorer (WISE). LRG targets are required to meet a set of color selection criteria and have z-band and i-band MODEL magnitudes z < 19.95 and 19.9 < i < 21.8, respectively. Our algorithm selects roughly 50 LRG targets per square degree, the great majority of which lie in the redshift range 0.6 < z < 1.0 (median redshift 0.71). We demonstrate that our methods are highly effective at eliminating stellar contamination and lower-redshift galaxies. We perform a number of tests using spectroscopic data from SDSS-III/BOSS to determine the redshift reliability of our target selection and its ability to meet the science requirements of eBOSS. The SDSS spectra are of high enough signal-to-noise ratio that at least 89% of the target sample yields secure redshift measurements. We also present tests of the uniformity and homogeneity of the sample, demonstrating that it should be clean enough for studies of the large-scale structure of the universe at higher redshifts than SDSS-III/BOSS LRGs reached.

From the Sloan Digital Sky Survey (SDSS) Data Release 12, which covers the full Baryonic Oscillation Spectroscopic Survey (BOSS) footprint, we investigate the possible variation of the fine-structure constant over cosmological time-scales. We analyse the largest quasar sample considered so far in the literature, which contains 13175 spectra (10363 from SDSS-III/BOSS DR12 + 2812 from SDSS-II DR7) with redshift $z<\,$1. We apply the emission-line method on the [O III] doublet (4960, 5008 A) and obtain $\Delta\alpha/\alpha= \left(0.9 \pm 1.8\right)\times10^{-5}$ for the relative variation of the fine-structure constant. We also investigate the possible sources of systematics: misidentification of the lines, sky OH lines, H$\,\beta$ and broad line contamination, Gaussian and Voigt fitting profiles, optimal wavelength range for the Gaussian fits, chosen polynomial order for the continuum spectrum, signal-to-noise ratio and good quality of the fits. The uncertainty of the measurement is dominated by the sky subtraction. The results presented in this work, being systematics limited, have sufficient statistics to constrain robustly the variation of the fine-structure constant in redshift bins ($\Delta z\approx$ 0.06) over the last 7.9 Gyr. In addition, we study the [Ne III] doublet (3870, 3969 A) present in 462 quasar spectra and discuss the systematic effects on using these emission lines to constrain the fine-structure constant variation. Better constraints on $\Delta\alpha/\alpha\ $ ($<$10$^{-6}$) using the emission-line method would be possible with high-resolution spectroscopy and large galaxy/qso surveys.

We present the rest-frame near-ultraviolet (NUV) spectroscopy of star-forming galaxies (SFGs) at 0.6<z<1.2 from the Extended Baryon Oscillation Spectroscopic Survey (eBOSS) in SDSS-IV. One of the eBOSS programs is to obtain 2 arcsec (about 15 kpc) fiber spectra of about 200,000 emission-line galaxies (ELGs) at redshift z>0.6. We use the data from the pilot observations of this program, including 8620 spectra of SFGs at 0.6<z<1.2. The median composite spectra of these SFGs at 2200 Ang < \lambda < 4000 Ang feature asymmetric, preferentially blueshifted non-resonant emission, Fe II*, and blueshifted resonant absorption, e.g., Fe II and Mg II, indicating ubiquitous outflows driven by star formation at these redshifts. For the absorption lines, we find a variety of velocity profiles with different degrees of blueshift. Comparing our new observations with the literature, we do not observe the non-resonant emission in the small-aperture (<40 pc) spectra of local star-forming regions with the Hubble Space Telescope, and find the observed line ratios in the SFG spectra to be different from those in the spectra of local star-forming regions, as well as those of quasar absorption-line systems in the same redshift range. We introduce an outflow model that can simultaneously explain the multiple observed properties and suggest that the variety of absorption velocity profiles and the line ratio differences are caused by scattered fluorescent emission filling in on top of the absorption in the large-aperture eBOSS spectra. We develop an observation-driven, model-independent method to correct the emission-infill to reveal the true absorption profiles. Our results show that eBOSS and future dark-energy surveys (e.g., DESI and PFS) will provide rich datasets of NUV spectroscopy for astrophysical applications.

The third generation of the Sloan Digital Sky Survey (SDSS-III) took data from 2008 to 2014 using the original SDSS wide-field imager, the original and an upgraded multi-object fiber-fed optical spectrograph, a new near-infrared high-resolution spectrograph, and a novel optical interferometer. All the data from SDSS-III are now made public. In particular, this paper describes Data Release 11 (DR11) including all data acquired through 2013 July, and Data Release 12 (DR12) adding data acquired through 2014 July (including all data included in previous data releases), marking the end of SDSS-III observing. Relative to our previous public release (DR10), DR12 adds one million new spectra of galaxies and quasars from the Baryon Oscillation Spectroscopic Survey (BOSS) over an additional 3000 sq. deg of sky, more than triples the number of H-band spectra of stars as part of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE), and includes repeated accurate radial velocity measurements of 5500 stars from the Multi-Object APO Radial Velocity Exoplanet Large-area Survey (MARVELS). The APOGEE outputs now include measured abundances of 15 different elements for each star. In total, SDSS-III added 2350 sq. deg of ugriz imaging; 155,520 spectra of 138,099 stars as part of the Sloan Exploration of Galactic Understanding and Evolution 2 (SEGUE-2) survey; 2,497,484 BOSS spectra of 1,372,737 galaxies, 294,512 quasars, and 247,216 stars over 9376 sq. deg; 618,080 APOGEE spectra of 156,593 stars; and 197,040 MARVELS spectra of 5,513 stars. Since its first light in 1998, SDSS has imaged over 1/3 of the Celestial sphere in five bands and obtained over five million astronomical spectra.

We present the [OII] luminosity function measured in the redshift range 0.1<z<1.65 with unprecedented depth and accuracy. Our measurements are based on medium resolution flux-calibrated spectra of emission line galaxies with the FORS2 instrument at VLT and with the SDSS-III/BOSS spectrograph. The FORS2 spectra and the corresponding catalog containing redshifts and line fluxes are released along with this paper. In this work we use a novel method to combine the aforementioned surveys with GAMA, zCOSMOS and VVDS, which have different target selection, producing a consistent weighting scheme to derive the [OII] luminosity function. The measured luminosity function is in good agreement with previous independent estimates. The comparison with two state-of-the-art semi-analytical models is good, which is encouraging for the production of mock catalogs of [OII] flux limited surveys. We observe the bright end evolution over 8.5 Gyr: we measure the decrease of log L* from 42.4 erg/s at redshift 1.44 to 41.2 at redshift 0.165 and we find that the faint end slope flattens when redshift decreases. This measurement confirms the feasibility of the target selection of future baryonic acoustic oscillation surveys aiming at observing [OII] flux limited samples.

We present a weak lensing mass map covering ~124 square degrees of the Canada-France-Hawaii Telescope Stripe 82 Survey (CS82). We study the statistics of rare peaks in the map, including peak abundance, the peak-peak correlation functions and the tangential-shear profiles around peaks. We find that the abundance of peaks detected in CS82 is consistent with predictions from a Lambda-CDM cosmological model, once noise effects are properly included. The correlation functions of peaks with different signal-to-noise ratio (SNR) are well described by power laws, and there is a clear cross-correlation between the Sloan Digital Sky Survey III/Constant Mass galaxies and high SNR peaks. The tangential-shear profiles around peaks increase with peak SNR. We fit analytical models to the tangential-shear profiles, including a projected singular isothermal sphere (SIS) model and a projected Navarro, Frenk & White (NFW) model, plus a two-halo term. For the high SNR peaks, the SIS model is rejected at ~3-sigma. The NFW model plus a two-halo term gives more acceptable fits to the data. Some peaks match the positions of optically detected clusters, while others are relatively dark. Comparing dark and matched peaks, we find a difference in lensing signal of a factor of 2, suggesting that about half of the dark peaks are false detections.

We investigate the luminosity and colour dependence of clustering of CMASS galaxies in the Sloan Digital Sky Survey-III Baryon Oscillation Spectroscopic Survey Tenth Data Release. The halo occupation distribution framework is adopted to model the projected two-point correlation function measurements on small and intermediate scales (from $0.02$ to $60\,h^{-1}{\rm {Mpc}}$) and to interpret the observed trends and infer the connection of galaxies to dark matter halos. We find that luminous red galaxies reside in massive halos of mass $M{\sim}10^{13}$--$10^{14}\,h^{-1}{\rm M_\odot}$ and more luminous galaxies are more clustered and hosted by more massive halos. The strong small-scale clustering requires a fraction of these galaxies to be satellites in massive halos, with the fraction at the level of 5--8 per cent and decreasing with luminosity. The characteristic mass of a halo hosting on average one satellite galaxy above a luminosity threshold is about a factor $8.7$ larger than that of a halo hosting a central galaxy above the same threshold. At a fixed luminosity, progressively redder galaxies are more strongly clustered on small scales, which can be explained by having a larger fraction of these galaxies in the form of satellites in massive halos. Our clustering measurements on scales below $0.4\,h^{-1}{\rm {Mpc}}$ allow us to study the small-scale spatial distribution of satellites inside halos. While the clustering of luminosity-threshold samples can be well described by a Navarro-Frenk-White (NFW) profile, that of the reddest galaxies prefers a steeper or more concentrated profile. Finally, we also use galaxy samples of constant number density at different redshifts to study the evolution of luminous galaxies, and find the clustering to be consistent with passive evolution in the redshift range of $0.5 \lesssim z \lesssim 0.6$.

The Sloan Digital Sky Survey (SDSS) has been in operation since 2000 April. This paper presents the tenth public data release (DR10) from its current incarnation, SDSS-III. This data release includes the first spectroscopic data from the Apache Point Observatory Galaxy Evolution Experiment (APOGEE), along with spectroscopic data from the Baryon Oscillation Spectroscopic Survey (BOSS) taken through 2012 July. The APOGEE instrument is a near-infrared R~22,500 300-fiber spectrograph covering 1.514--1.696 microns. The APOGEE survey is studying the chemical abundances and radial velocities of roughly 100,000 red giant star candidates in the bulge, bar, disk, and halo of the Milky Way. DR10 includes 178,397 spectra of 57,454 stars, each typically observed three or more times, from APOGEE. Derived quantities from these spectra (radial velocities, effective temperatures, surface gravities, and metallicities) are also included.DR10 also roughly doubles the number of BOSS spectra over those included in the ninth data release. DR10 includes a total of 1,507,954 BOSS spectra, comprising 927,844 galaxy spectra; 182,009 quasar spectra; and 159,327 stellar spectra, selected over 6373.2 square degrees.

The next generation of wide-field spectroscopic redshift surveys will map the large-scale galaxy distribution in the redshift range 0.7< z<2 to measure baryonic acoustic oscillations (BAO). The primary optical signature used in this redshift range comes from the [OII] emission line doublet, which provides a unique redshift identification that can minimize confusion with other single emission lines. To derive the required spectrograph resolution for these redshift surveys, we simulate observations of the [OII] (3727,3729) doublet for various instrument resolutions, and line velocities. We foresee two strategies about the choice of the resolution for future spectrographs for BAO surveys. For bright [OII] emitter surveys ([OII] flux ~30.10^{-17} erg /cm2/s like SDSS-IV/eBOSS), a resolution of R~3300 allows the separation of 90 percent of the doublets. The impact of the sky lines on the completeness in redshift is less than 6 percent. For faint [OII] emitter surveys ([OII] flux ~10.10^{-17} erg /cm2/s like DESi), the detection improves continuously with resolution, so we recommend the highest possible resolution, the limit being given by the number of pixels (4k by 4k) on the detector and the number of spectroscopic channels (2 or 3).

Large sets of objects with spectroscopic redshift measurements will be needed for imaging dark energy experiments to achieve their full potential, serving two goals:_training_, i.e., the use of objects with known redshift to develop and optimize photometric redshift algorithms; and_calibration_, i.e., the characterization of moments of redshift (or photo-z error) distributions. Better training makes cosmological constraints from a given experiment stronger, while highly-accurate calibration is needed for photo-z systematics not to dominate errors. In this white paper, we investigate the required scope of spectroscopic datasets which can serve both these purposes for ongoing and next-generation dark energy experiments, as well as the time required to obtain such data with instruments available in the next decade. Large time allocations on kilo-object spectrographs will be necessary, ideally augmented by infrared spectroscopy from space. Alternatively, precision calibrations could be obtained by measuring cross-correlation statistics using samples of bright objects from a large baryon acoustic oscillation experiment such as DESI. We also summarize the additional work on photometric redshift methods needed to prepare for ongoing and future dark energy experiments.

We estimate the galaxy stellar mass function and stellar mass density for star-forming and quiescent galaxies with 0.2<z<4. We construct a deep K<24 sample of 220000 galaxies selected using the UltraVISTA DR1 data release. Our analysis is based on precise 30-band photometric redshifts. By comparing these photometric redshifts with 10800 spectroscopic redshifts from the zCOSMOS bright and faint surveys, we find a precision of sigma(dz/(1+z))=0.008 at i<22.5 and sigma(dz/(1+zs))=0.03 at 1.5<z<4. We derive the stellar mass function and correct for the Eddington bias. We find a mass-dependent evolution of the global and star-forming populations. This mass-dependent evolution is a direct consequence of the star formation being quenched in galaxies more massive than M>10^10.7Msun. For the mass function of the quiescent galaxies, we do not find any significant evolution of the high-mass end at z<1; however we observe a clear flattening of the faint-end slope. From z~3 to z~1, the density of quiescent galaxies increases over the entire mass range. Their comoving stellar mass density increases by 1.6 dex between z~3 and z~1 and by less than 0.2dex at z<1. We infer the star formation history from the mass density evolution and we find an excellent agreement with instantaneous star formation rate measurements at z<1.5, while we find differences of 0.2dex at z>1.5 consistent with the expected uncertainties. We also develop a new method to infer the specific star formation rate from the mass function of star-forming galaxies. We find that the specific star formation rate of 10^10Msun galaxies increases continuously in the redshift range 1<z<4. Finally, we compare our results with a semi-analytical model and find that these models overestimate the density of low mass quiescent galaxies by an order of magnitude, while the density of low-mass star-forming galaxies is successfully reproduced.

The baryon acoustic oscillation (BAO) feature in the two-point correlation function of galaxies supplies a standard ruler to probe the expansion history of the Universe. We study here several galaxy selection schemes, aiming at building an emission-line galaxy (ELG) sample in the redshift range $0.6<z<1.7$, that would be suitable for future BAO studies, providing a highly biased galaxy sample. We analyse the angular galaxy clustering of galaxy selections at the redshifts 0.5, 0.7, 0.8, 1 and 1.2 and we combine this analysis with a halo occupation distribution (HOD) model to derive the properties of the haloes these galaxies inhabit, in particular the galaxy bias on large scales. We also perform a weak lensing analysis (aperture statistics) to extract the galaxy bias and the cross-correlation coefficient and compare to the HOD prediction. We apply this analysis on a data set composed of the photometry of the deep co-addition on Sloan Digital Sky Survey (SDSS) Stripe 82 (225 deg$^2$), of Canda-France-Hawai Telescope/Stripe 82 deep \emph{i}-band weak lensing survey and of the {\it Wide-Field Infrared Survey Explorer}infrared photometric band W1. The analysis on the SDSS-III/constant mass galaxies selection at $z=0.5$ is in agreement with previous studies on the tracer, moreover we measure its cross-correlation coefficient $r=1.16\pm0.35$. For the higher redshift bins, we confirm the trends that the brightest galaxy populations selected are strongly biased ($b>1.5$), but we are limited by current data sets depth to derive precise values of the galaxy bias. A survey using such tracers of the mass field will guarantee a high significance detection of the BAO.

The Baryon Acoustic Oscillation (BAO) feature in the power spectrum of galaxies can be used as a standard ruler to probe the accelerated expansion of the Universe. In this paper, we study several galaxy selection schemes aiming at building an emission-line galaxy (ELG) sample in the redshift range $0.6 < z < 1.7$, that would be suitable for future BAO studies using the Baryonic Oscillation Spectroscopic Survey (BOSS) spectrograph on the Sloan Digital Sky Survey (SDSS) telescope. We explore two different color selections using both the SDSS and the Canada-France-Hawaii Telescope Legacy Survey (CFHTLS) photometry in the $u, g, r, i$ bands and evaluate their performance for selecting bright ELG. This study confirms the feasibility of massive ELG surveys using the BOSS spectrographs on the SDSS telescope for a BAO detection at redshift $z\sim1$, in particular for the proposed eBOSS experiment.

We perform a spectroscopic analysis of 492,450 galaxy spectra from the first two years of observations of the Sloan Digital Sky Survey-III/Baryonic Oscillation Spectroscopic Survey (BOSS) collaboration. This data set has been released in the ninth SDSS data release, the first public data release of BOSS spectra. We show that the typical signal-to-noise ratio of BOSS spectra is sufficient to measure stellar velocity dispersion and emission line fluxes for individual objects. The typical velocity dispersion of a BOSS galaxy is 240 km/s, with an accuracy of better than 30 per cent for 93 per cent of BOSS galaxies. The distribution in velocity dispersion is redshift independent between redshifts 0.15 and 0.7, which reflects the survey design targeting massive galaxies with an approximately uniform mass distribution in this redshift interval. The majority of BOSS galaxies lack detectable emission lines. We analyse the emission line properties and present diagnostic diagrams using the emission lines [OII], Hbeta, [OIII], Halpha, and [NII] (detected in about 4 per cent of the galaxies). We show that the emission line properties are strongly redshift dependent and that there is a clear correlation between observed frame colours and emission line properties. Within in the low-z sample around 0.15<z<0.3, half of the emission-line galaxies have LINER-like emission line ratios, followed by Seyfert-AGN dominated spectra, and only a small fraction of a few per cent are purely star forming galaxies. AGN and LINER-like objects, instead, are less prevalent in the high-z sample around 0.4<z<0.7, where more than half of the emission line objects are star forming. This is a pure selection effect caused by the non-detection of weak Hbeta emission lines in the BOSS spectra. Finally, we show that star forming, AGN and emission line free galaxies are well separated in the g-r vs r-i target selection diagram.

The Baryon Oscillation Spectroscopic Survey (BOSS) is designed to measure the scale of baryon acoustic oscillations (BAO) in the clustering of matter over a larger volume than the combined efforts of all previous spectroscopic surveys of large scale structure. BOSS uses 1.5 million luminous galaxies as faint as i=19.9 over 10,000 square degrees to measure BAO to redshifts z<0.7. Observations of neutral hydrogen in the Lyman alpha forest in more than 150,000 quasar spectra (g<22) will constrain BAO over the redshift range 2.15<z<3.5. Early results from BOSS include the first detection of the large-scale three-dimensional clustering of the Lyman alpha forest and a strong detection from the Data Release 9 data set of the BAO in the clustering of massive galaxies at an effective redshift z = 0.57. We project that BOSS will yield measurements of the angular diameter distance D_A to an accuracy of 1.0% at redshifts z=0.3 and z=0.57 and measurements of H(z) to 1.8% and 1.7% at the same redshifts. Forecasts for Lyman alpha forest constraints predict a measurement of an overall dilation factor that scales the highly degenerate D_A(z) and H^{-1}(z) parameters to an accuracy of 1.9% at z~2.5 when the survey is complete. Here, we provide an overview of the selection of spectroscopic targets, planning of observations, and analysis of data and data quality of BOSS.

The Baryon Acoustic Oscillation (BAO) feature in the power spectrum of galaxies provides a standard ruler to probe the accelerated expansion of the Universe. The current surveys covering a comoving volume sufficient to unveil the BAO scale are limited to redshift $z \lesssim 0.7$. In this paper, we study several galaxy selection schemes aiming at building an emission-line-galaxy (ELG) sample in the redshift range $0.6<z<1.7$, that would be suitable for future BAO studies using the Baryonic Oscillation Spectroscopic Survey (BOSS) spectrograph on the Sloan Digital Sky Survey (SDSS) telescope. We explore two different colour selections using both the SDSS and the Canada France Hawai Telescope Legacy Survey (CFHT-LS) photometry in the u, g, r, and i bands and evaluate their performance selecting luminous ELG. From about 2,000 ELG, we identified a selection scheme that has a 75 percent redshift measurement efficiency. This result confirms the feasibility of massive ELG surveys using the BOSS spectrograph on the SDSS telescope for a BAO detection at redshift $z\sim1$, in particular the proposed eBOSS experiment, which plans to use the SDSS telescope to combine the use of the BAO ruler with redshift space distortions using emission line galaxies and quasars in the redshift $0.6<z<2.2$.

The Sloan Digital Sky Survey III (SDSS-III) presents the first spectroscopic data from the Baryon Oscillation Spectroscopic Survey (BOSS). This ninth data release (DR9) of the SDSS project includes 535,995 new galaxy spectra (median z=0.52), 102,100 new quasar spectra (median z=2.32), and 90,897 new stellar spectra, along with the data presented in previous data releases. These spectra were obtained with the new BOSS spectrograph and were taken between 2009 December and 2011 July. In addition, the stellar parameters pipeline, which determines radial velocities, surface temperatures, surface gravities, and metallicities of stars, has been updated and refined with improvements in temperature estimates for stars with T_eff<5000 K and in metallicity estimates for stars with [Fe/H]>-0.5. DR9 includes new stellar parameters for all stars presented in DR8, including stars from SDSS-I and II, as well as those observed as part of the SDSS-III Sloan Extension for Galactic Understanding and Exploration-2 (SEGUE-2). The astrometry error introduced in the DR8 imaging catalogs has been corrected in the DR9 data products. The next data release for SDSS-III will be in Summer 2013, which will present the first data from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) along with another year of data from BOSS, followed by the final SDSS-III data release in December 2014.

Building on the legacy of the Sloan Digital Sky Survey (SDSS-I and II), SDSS-III is a program of four spectroscopic surveys on three scientific themes: dark energy and cosmological parameters, the history and structure of the Milky Way, and the population of giant planets around other stars. In keeping with SDSS tradition, SDSS-III will provide regular public releases of all its data, beginning with SDSS DR8 (which occurred in Jan 2011). This paper presents an overview of the four SDSS-III surveys. BOSS will measure redshifts of 1.5 million massive galaxies and Lya forest spectra of 150,000 quasars, using the BAO feature of large scale structure to obtain percent-level determinations of the distance scale and Hubble expansion rate at z<0.7 and at z~2.5. SEGUE-2, which is now completed, measured medium-resolution (R=1800) optical spectra of 118,000 stars in a variety of target categories, probing chemical evolution, stellar kinematics and substructure, and the mass profile of the dark matter halo from the solar neighborhood to distances of 100 kpc. APOGEE will obtain high-resolution (R~30,000), high signal-to-noise (S/N>100 per resolution element), H-band (1.51-1.70 micron) spectra of 10^5 evolved, late-type stars, measuring separate abundances for ~15 elements per star and creating the first high-precision spectroscopic survey of all Galactic stellar populations (bulge, bar, disks, halo) with a uniform set of stellar tracers and spectral diagnostics. MARVELS will monitor radial velocities of more than 8000 FGK stars with the sensitivity and cadence (10-40 m/s, ~24 visits per star) needed to detect giant planets with periods up to two years, providing an unprecedented data set for understanding the formation and dynamical evolution of giant planet systems. (Abridged)

The Sloan Digital Sky Survey (SDSS) started a new phase in August 2008, with new instrumentation and new surveys focused on Galactic structure and chemical evolution, measurements of the baryon oscillation feature in the clustering of galaxies and the quasar Ly alpha forest, and a radial velocity search for planets around ~8000 stars. This paper describes the first data release of SDSS-III (and the eighth counting from the beginning of the SDSS). The release includes five-band imaging of roughly 5200 deg^2 in the Southern Galactic Cap, bringing the total footprint of the SDSS imaging to 14,555 deg^2, or over a third of the Celestial Sphere. All the imaging data have been reprocessed with an improved sky-subtraction algorithm and a final, self-consistent photometric recalibration and flat-field determination. This release also includes all data from the second phase of the Sloan Extension for Galactic Understanding and Evolution (SEGUE-2), consisting of spectroscopy of approximately 118,000 stars at both high and low Galactic latitudes. All the more than half a million stellar spectra obtained with the SDSS spectrograph have been reprocessed through an improved stellar parameters pipeline, which has better determination of metallicity for high metallicity stars.