The Lyman-$\alpha$ (hereafter Ly-$\alpha$) forest is a probe of large-scale matter density fluctuations at high redshift, $z > 2.1$. It consists of HI absorption spectra along individual lines-of-sight. If the line-of-sight density is large enough, 3D maps of HI absorption can be inferred by tomographic reconstruction. In this article, we investigate the Ly-$\alpha$ forest available in the Stripe 82 field ($220\,\mathrm{deg^{2}}$), based on the quasar spectra from SDSS Data Release DR16. The density of observed quasar spectra is $37\,\mathrm{deg^{-2}}$ with a mean pixel signal-to-noise ratio of two per angstrom. This study provides an intermediate case between the average SDSS density and that of the much denser but smaller CLAMATO survey. We derive a 3D map of large-scale matter fluctuations from these data, using a Wiener filter technique. The total volume of the map is $0.94\,\mathrm{h^{-3} Gpc^{3}}$. Its resolution is $13\,\mathrm{h^{-1} Mpc}$, which is related to the mean transverse distance between nearest lines-of-sight. From this map, we provide a catalog of voids and protocluster candidates in the cosmic web. The map-making and void catalog are compared to simulated eBOSS Stripe 82 observations. A stack over quasar positions provides a visualization of the Ly-$\alpha$ quasar cross-correlation. This tomographic reconstruction constitutes the largest-volume high-redshift 3D map of matter fluctuations.

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

We analyze the anisotropic clustering of the Sloan Digital Sky Survey-IV Extended Baryon Oscillation Spectroscopic Survey (eBOSS) Luminous Red Galaxy Data Release 14 (DR14) sample combined with Baryon Oscillation Spectroscopic Survey (BOSS) CMASS sample of galaxies in the redshift range 0.6$<z<$1.0, which consists of 80,118 galaxies from eBOSS and 46,439 galaxies from the BOSS-CMASS sample. The eBOSS-CMASS Luminous Red Galaxy sample has a sky coverage of 1,844 deg$^2$, with an effective volume of 0.9 Gpc$^3$. The analysis was made in configuration space using a Legendre multipole expansion. The Redshift Space Distortion signal is modeled as a combination of the Convolution Lagrangian Perturbation Model and the Gaussian Streaming Model. We constrain the logarithmic growth of structure times the amplitude of dark matter density fluctuations, $f (z_{\rm eff})\sigma_8(z_{\rm eff})=0.454 \pm0.139 $, and the Alcock-Paczynski dilation scales which constraints the angular diameter distance $D_A(z_{eff})=1466.5 \pm 136.6 (r_s/r_s^{\rm fid})$ and $H(z_{\rm eff})=105.8 \pm 16 (r_s^{\rm fid}/r_s) \mathrm{km\,s^{-1}\,Mpc^{-1}}$, 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 at an effective redshift $z_{\rm eff}=0.72$. These results are in full agreement with the current $\Lambda$-Cold Dark Matter ($\Lambda$-CDM) cosmological model inferred from Planck measurements. This study is the first eBOSS LRG full-shape analysis i.e. including Redshift-Space Distortions (RSD) simultaneously with the Alcock-Paczynski (AP) effect and the Baryon Acoustic Oscillation (BAO) scale.

We measure the star formation histories (SFH) of a sample of low-mass galaxies with $M_\ast<10^9M_\odot$ from the SDSS-IV MaNGA survey. The large number of IFU spectra for each galaxy are either combined to reach a high signal to noise ratio or used to investigate spatial variations. We use Bayesian inferences based on full spectrum fitting. Our analysis based on Bayesian evidence ratio indicates a strong preference for a model that allows the presence of an old stellar population, and that an improper model for the SFH can significantly underestimate the old population in these galaxies. The addition of NIR photometry to the constraining data can further distinguish between different SFH model families and significantly tighten the constraints on the mass fraction in the old population. On average more than half of the stellar mass in present-day low-mass galaxies formed 8 Gyrs ago, while about 30\% within the past 4 Gyrs. Satellite galaxies on average have formed their stellar mass earlier than central galaxies, and stars in the outer regions of galaxy are younger than those in the central part. Our results suggest that most of the low-mass galaxies have an early episode of active star formation that produces a large fraction of their present stellar mass.

The Wide Field Infrared Survey Telescope (WFIRST), with its field of view and high sensitivity will make surveys of cosmological large-scale structure possible at high redshifts. We investigate the possibility of detecting Baryon Acoustic Oscillations (BAO) at redshifts $z>7.6$ for use as a standard ruler. We use data from the hydrodynamic simulation \textsc{BlueTides} in conjunction with the gigaparsec-scale Outer Rim simulation and a model for patchy reionization to create mock WFIRST High Latitude Survey grism data for \lya\ emission line selected galaxies at redshifts $z=7.4$ to $z=10$, covering 2280 square degrees. We measure the monopoles of galaxies in the mock catalogues and fit the BAO features. We find that for a line flux of $L = 7\times 10^{-17} \ {\rm erg/s/cm}^{2}$, the $5 \sigma$ detection limit for the current design, the BAO feature is partially detectable (measured in three out of four survey quadrants analysed independently). The resulting root mean square error on the angular diameter distance to $z=7.7$ is 7.9$\%$. If we improve the detection sensitivity by a factor of two (i.e. $L = 3.5\times 10^{-17} \ {\rm erg/s/cm}^{2}$), the distance error reduces to $1.4\%$. We caution that many more factors are yet to be modelled, including dust obscuration, the damping wing due to the intergalactic medium, and low redshift interlopers. If these issues do not strongly affect the results, or different observational techniques (such as use of multiple lines) can mitigate them, WFIRST or similar instruments may be able to constrain the angular diameter distance to the high redshift Universe.

We present catalogues of cosmic voids identified in the distribution of Luminous Red Galaxies (LRGs) and Quasi Stellar Objects (QSOs) in the fourteenth data release (DR14) of the extended Baryon Oscillation Spectroscopic Survey (eBOSS). We perform a multivariate analysis to assess the level of contamination in these catalogues by spurious Poisson underdensities. We find that the LRG void catalogue is largely free from contamination but that the QSO catalogue may be heavily contaminated. We analyse the multipoles of the void-galaxy cross-correlation function in these catalogues to obtain constraints on the growth rate of structure around voids. We find a value of $\beta(z=0.703)=0.58^{+0.33}_{-0.28}$ for the LRG voids and $\beta(z=1.53)=0.15^{+0.13}_{-0.12}$ for the QSO voids.

We present constraints on neutrino masses, the primordial fluctuation spectrum from inflation, and other parameters of the $\Lambda$CDM model, using the one-dimensional Ly$\alpha$-forest power spectrum measured by Palanque-Delabrouille et al. (2013) from SDSS-III/BOSS, complemented by Planck 2015 cosmic microwave background (CMB) data and other cosmological probes. This paper improves on the previous analysis by Palanque-Delabrouille et al. (2015) by using a more powerful set of calibrating hydrodynamical simulations that reduces uncertainties associated with resolution and box size, by adopting a more flexible set of nuisance parameters for describing the evolution of the intergalactic medium, by including additional freedom to account for systematic uncertainties, and by using Planck 2015 constraints in place of Planck 2013. Fitting Ly$\alpha$ data alone leads to cosmological parameters in excellent agreement with the values derived independently from CMB data, except for a weak tension on the scalar index $n_s$. Combining BOSS Ly$\alpha$ with Planck CMB constrains the sum of neutrino masses to $\sum m_\nu < 0.12$ eV (95\% C.L.) including all identified systematic uncertainties, tighter than our previous limit (0.15 eV) and more robust. Adding Ly$\alpha$ data to CMB data reduces the uncertainties on the optical depth to reionization $\tau$, through the correlation of $\tau$ with $\sigma_8$. Similarly, correlations between cosmological parameters help in constraining the tensor-to-scalar ratio of primordial fluctuations $r$. The tension on $n_s$ can be accommodated by allowing for a running ${\mathrm d}n_s/{\mathrm d}\ln k$. Allowing running as a free parameter in the fits does not change the limit on $\sum m_\nu$. We discuss possible interpretations of these results in the context of slow-roll inflation.

We present a measurement of the 1D Ly$\alpha$ forest flux power spectrum, using the complete Baryon Oscillation Spectroscopic Survey (BOSS) and first extended-BOSS (eBOSS) quasars at $z_{\rm qso}>2.1$, corresponding to the fourteenth data release (DR14) of the Sloan Digital Sky Survey (SDSS). Our results cover thirteen bins in redshift from $z_{\rm Ly\alpha}=2.2$ to 4.6, and scales up to $k=0.02\rm \,(km/s)^{-1}$. From a parent sample of 180,413 visually inspected spectra, we selected the 43,751 quasars with the best quality; this data set improves the previous result from the ninth data release (DR9), both in statistical precision (achieving a reduction by a factor of two) and in redshift coverage. We also present a thorough investigation of identified sources of systematic uncertainties that affect the measurement. The resulting 1D power spectrum of this work is in excellent agreement with the one from the BOSS DR9 data.

We present constraints on local primordial non-Gaussianity (PNG), parametrized through $f^{\rm loc}_{\rm NL}$, using the Sloan Digital Sky Survey IV extended Baryon Oscillation Spectroscopic Survey Data Release 14 quasar sample. We measure and analyze the anisotropic clustering of the quasars in Fourier space, testing for the scale-dependent bias introduced by primordial non-Gaussianity on large scales. We derive and employ a power spectrum estimator using optimal weights that account for the redshift evolution of the PNG signal. We find constraints of $-51<f^{\rm loc}_{\rm NL}<21$ at 95% confidence level. These are amont the tightest constraints from Large Scale Structure (LSS) data. Our redshift weighting improves the error bar by 15% in comparison to the unweighted case. If quasars have lower response to PNG, the constraint degrades to $-81<f^{\rm loc}_{\rm NL}<26$, with a 40% improvement over the standard approach. We forecast that the full eBOSS dataset could reach $\sigma_{f^{\rm loc}_{\rm NL}}\simeq 5\text{-}8$ using optimal methods and full range of scales.

The expansion of the Universe is understood to have accelerated during two epochs: in its very first moments during a period of Inflation and much more recently, at $z < 1$, when Dark Energy is hypothesized to drive cosmic acceleration. The undiscovered mechanisms behind these two epochs represent some of the most important open problems in fundamental physics. The large cosmological volume at $2 < z < 5$, together with the ability to efficiently target high-$z$ galaxies with known techniques, enables large gains in the study of Inflation and Dark Energy. A future spectroscopic survey can test the Gaussianity of the initial conditions up to a factor of ~50 better than our current bounds, crossing the crucial theoretical threshold of $\sigma(f_{NL}^{\rm local})$ of order unity that separates single field and multi-field models. Simultaneously, it can measure the fraction of Dark Energy at the percent level up to $z = 5$, thus serving as an unprecedented test of the standard model and opening up a tremendous discovery space.

Our current understanding of the Universe is established through the pristine measurements of structure in the cosmic microwave background (CMB) and the distribution and shapes of galaxies tracing the large scale structure (LSS) of the Universe. One key ingredient that underlies cosmological observables is that the field that sources the observed structure is assumed to be initially Gaussian with high precision. Nevertheless, a minimal deviation from Gaussianityis perhaps the most robust theoretical prediction of models that explain the observed Universe; itis necessarily present even in the simplest scenarios. In addition, most inflationary models produce far higher levels of non-Gaussianity. Since non-Gaussianity directly probes the dynamics in the early Universe, a detection would present a monumental discovery in cosmology, providing clues about physics at energy scales as high as the GUT scale.

The hot dense environment of the early universe is known to have produced large numbers of baryons, photons, and neutrinos. These extreme conditions may have also produced other long-lived species, including new light particles (such as axions or sterile neutrinos) or gravitational waves. The gravitational effects of any such light relics can be observed through their unique imprint in the cosmic microwave background (CMB), the large-scale structure, and the primordial light element abundances, and are important in determining the initial conditions of the universe. We argue that future cosmological observations, in particular improved maps of the CMB on small angular scales, can be orders of magnitude more sensitive for probing the thermal history of the early universe than current experiments. These observations offer a unique and broad discovery space for new physics in the dark sector and beyond, even when its effects would not be visible in terrestrial experiments or in astrophysical environments. A detection of an excess light relic abundance would be a clear indication of new physics and would provide the first direct information about the universe between the times of reheating and neutrino decoupling one second later.

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.

Opening up a new window of millimeter-wave observations that span frequency bands in the range of 30 to 500 GHz, survey half the sky, and are both an order of magnitude deeper (about 0.5 uK-arcmin) and of higher-resolution (about 10 arcseconds) than currently funded surveys would yield an enormous gain in understanding of both fundamental physics and astrophysics. In particular, such a survey would allow for major advances in measuring the distribution of dark matter and gas on small-scales, and yield needed insight on 1.) dark matter particle properties, 2.) the evolution of gas and galaxies, 3.) new light particle species, 4.) the epoch of inflation, and 5.) the census of bodies orbiting in the outer Solar System.

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.

We develop a new method, which is based on the optimal redshift weighting scheme, to extract the maximal tomographic information of baryonic acoustic oscillations (BAO) and redshift space distortions (RSD) from the extended Baryon Oscillation Spectroscopic Survey (eBOSS) Data Release 14 quasar (DR14Q) survey. We validate our method using the EZ mocks, and apply our pipeline to the eBOSS DR14Q sample in the redshift range of $0.8<z<2.2$. We report a joint measurement of $f\sigma_8$ and two-dimensional BAO parameters $D_{\rm A}$ and $H$ at four effective redshifts of $z_{\rm eff}=0.98, 1.23, 1.52$ and $1.94$, and provide the full data covariance matrix. Using our measurement combined with BOSS DR12, MGS and 6dFGS BAO measurements, we find that the existence of dark energy is supported by observations at a $7.4\sigma$ significance level. Combining our measurement with BOSS DR12 and Planck observations, we constrain the gravitational growth index to be $\gamma=0.580\pm0.082$, which is fully consistent with the prediction of general relativity. This paper is part of a set that analyses the eBOSS DR14 quasar sample.

The $\beta$-skeleton is a mathematical method to construct graphs from a set of points that has been widely applied in the areas of image analysis, machine learning, visual perception, and pattern recognition. In this work, we apply the $\beta$-skeleton to study the cosmic web. We use this tool on observed and simulated data to identify the filamentary structures and characterize the statistical properties of the skeleton. In particular, we compare the $\beta$-skeletons built from SDSS-III galaxies to those obtained from MD-PATCHY mocks, and also to mocks directly built from the Big MultiDark $N$-body simulation. We find that the $\beta$-skeleton is able to reveal the underlying structures in observed and simulated samples without any parameter fine-tuning. A different degree of sparseness can be obtained by adjusting the value of $\beta$; in addition, the statistical properties of the length and direction of the skeleton connections show a clear dependence on redshift space distortions (RSDs), cosmological effects and galaxy bias. We also find that the $N$-body simulation accurately reproduces the RSD effect in the data, while the MD-PATCHY mocks appear to underestimate its magnitude. Our proof-of-concept study shows that the statistical properties of the $\beta$-skeleton can be used to probe cosmological parameters and galaxy evolution.

We explore the cosmological implications of anisotropic clustering measurements of the quasar sample from Data Release 14 of the Sloan Digital Sky Survey IV Extended Baryon Oscillation Spectroscopic Survey (eBOSS) in configuration space. The $\sim 147,000$ quasar sample observed by eBOSS offers a direct tracer of the density field and bridges the gap of previous BAO measurements between redshift $0.8<z<2.2$. By analysing the two-point correlation function characterized by clustering wedges $\xi_{\rm w_i}(s)$ and multipoles $\xi_{\ell}(s)$, we measure the angular diameter distance, Hubble parameter and cosmic structure growth rate. We define a systematic error budget for our measurements based on the analysis of $N$-body simulations and mock catalogues. Based on the DR14 large scale structure quasar sample at the effective redshift $z_{\rm eff}=1.52$, we find the growth rate of cosmic structure $f\sigma_8(z_{\rm eff})=0.396\pm 0.079$, and the geometric parameters $D_{\rm V}(z)/r_{\rm d}=26.47\pm 1.23$, and $F_{\rm AP}(z)=2.53\pm 0.22$, where the uncertainties include both statistical and systematic errors. These values are in excellent agreement with the best-fitting standard ${\rm \Lambda CDM}$ model to the latest cosmic microwave background data from Planck.

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 investigate the potential sources of theoretical systematics in the anisotropic Baryon Acoustic Oscillation (BAO) distance scale measurements from the clustering of galaxies in configuration space using the final Data Release (DR12) of the Baryon Oscillation Spectroscopic Survey (BOSS). We perform a detailed study of the impact on BAO measurements from choices in the methodology such as fiducial cosmology, clustering estimators, random catalogues, fitting templates, and covariance matrices. The theoretical systematic uncertainties in BAO parameters are found to be 0.002 in the isotropic dilation $\alpha$ and 0.003 in the quadrupolar dilation $\epsilon$. The leading source of systematic uncertainty is related to the reconstruction techniques. Theoretical uncertainties are sub-dominant compared with the statistical uncertainties for BOSS survey, accounting $0.2\sigma_{stat}$ for $\alpha$ and $0.25\sigma_{stat}$ for $\epsilon$ ($\sigma_{\alpha,stat} \sim$0.010 and $\sigma_{\epsilon,stat}\sim$ 0.012 respectively). We also present BAO-only distance scale constraints from the anisotropic analysis of the correlation function. Our constraints on the angular diameter distance $D_A(z)$ and the Hubble parameter $H(z)$, including both statistical and theoretical systematic uncertainties, are 1.5\% and 2.8\% at $z_{\rm eff}=0.38$, 1.4\% and 2.4\% at $z_{\rm eff}=0.51$, and 1.7\% and 2.6\% at $z_{\rm eff}=0.61$. This paper is part of a set that analyzes the final galaxy clustering dataset from BOSS. The measurements and likelihoods presented here are cross-checked with other BAO analysis in \citet{Acacia16}. The systematic error budget concerning the methodology on post-reconstruction BAO analysis presented here is used in \citet{Acacia16} to produce the final cosmological constraints from BOSS.

We present the Data Release 14 Quasar catalog (DR14Q) from the extended Baryon Oscillation Spectroscopic Survey (eBOSS) of the Sloan Digital Sky Survey IV (SDSS-IV). This catalog includes all SDSS-IV/eBOSS objects that were spectroscopically targeted as quasar candidates and that are confirmed as quasars via a new automated procedure combined with a partial visual inspection of spectra, have luminosities $M_{\rm i} \left[ z=2 \right] < -20.5$ (in a $\Lambda$CDM cosmology with $H_0 = 70 \ {\rm km \ s^{-1} \ Mpc ^{-1}}$, $\Omega_{\rm M} = 0.3$, and $\Omega_{\rm \Lambda} = 0.7$), and either display at least one emission line with a full width at half maximum (FWHM) larger than $500 \ {\rm km \ s^{-1}}$ or, if not, have interesting/complex absorption features. The catalog also includes previously spectroscopically-confirmed quasars from SDSS-I, II and III. The catalog contains 526,356 quasars 144,046 are new discoveries since the beginning of SDSS-IV) detected over 9,376 deg$^2$ (2,044 deg$^2$ having new spectroscopic data available) with robust identification and redshift measured by a combination of principal component eigenspectra. The catalog is estimated to have about 0.5% contamination. The catalog identifies 21,877 broad absorption line quasars and lists their characteristics. For each object, the catalog presents SDSS five-band CCD-based photometry with typical accuracy of 0.03 mag. The catalog also contains X-ray, ultraviolet, near-infrared, and radio emission properties of the quasars, when available, from other large-area surveys.

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.

We present an anisotropic analysis of Baryon Acoustic Oscillation (BAO) signal from the SDSS-IV extended Baryon Oscillation Spectroscopic Survey (eBOSS) Data Release 14 (DR14) quasar sample. The sample consists of 147,000 quasars distributed over a redshift range of $0.8 < z < 2.2$. We apply the redshift weights technique to the clustering of quasars in this sample and achieve a 4.6 per cent measurement of the angular distance measurement $D_M$ at $z = 2.2$ and Hubble parameter $H$ at $z=0.8$. We parameterize the distance-redshift relation, relative to a fiducial model, as a quadratic expansion. The coefficients of this expansion are used to reconstruct the distance-redshift relation and obtain distance and Hubble parameter measurements at all redshifts within the redshift range of the sample. Reporting the result at two characteristic redshifts, we determine $D_M(z=1) = 3405\pm305 \ (r_{\rm d} / r_{\rm d, fid}) \ {\rm Mpc}$, $H(z=1) = 120.7\pm 7.3 \ (r_{\rm d,fid} / r_{\rm d}) \ {\rm km} \ {\rm s}^{-1}{\rm Mpc}^{-1}$ and $D_M(z=2) = 5325\pm249 \ (r_{\rm d} / r_{\rm d, fid}) \ {\rm Mpc}$, $H(z=2) = 189.9\pm 32.9 \ (r_{\rm d,fid} / r_{\rm d}) \ {\rm km} \ {\rm s}^{-1}{\rm Mpc}^{-1}$. These measurements are highly correlated. We assess the outlook of BAO analysis from the final quasar sample by testing the method on a set of mocks that mimic the noise level in the final sample. We demonstrate on these mocks that redshift weighting shrinks the measurement error by over 25 per cent on average. We conclude redshift weighting can bring us closer to the cosmological goal of the final quasar sample.

We analyse the clustering of the Sloan Digital Sky Survey IV extended Baryon Oscillation Spectroscopic Survey Data Release 14 quasar sample (DR14Q). We measure the redshift space distortions using the power spectrum monopole, quadrupole and hexadecapole inferred from 148,659 quasars between redshifts 0.8 and 2.2 covering a total sky footprint of 2112.9 deg$^2$. We constrain the logarithmic growth of structure times the amplitude of dark matter density fluctuations, $f\sigma_8$, and the Alcock-Paczynski dilation scales which allow constraints to be placed on the angular diameter distance $D_A(z)$ and the Hubble $H(z)$ parameter. At the effective redshift of $z_{\rm eff}=1.52$, $f\sigma_8(z_{\rm eff})=0.420\pm0.076$, $H(z_{\rm eff})=[162\pm 12]\, (r_s^{\rm fid}/r_s)\,{\rm km\, s}^{-1}{\rm Mpc}^{-1}$, and $D_A(z_{\rm eff})=[1.85\pm 0.11]\times10^3\,(r_s/r_s^{\rm fid})\,{\rm Mpc}$, where $r_s$ is the comoving sound horizon at the baryon drag epoch and the superscript `fid' stands for its fiducial value. The errors take into account the full error budget, including systematics and statistical contributions. These results are in full agreement with the current $\Lambda$-Cold Dark Matter ($\Lambda$CDM) cosmological model inferred from Planck measurements. Finally, we compare our measurements with other eBOSS companion papers and find excellent agreement, demonstrating the consistency and complementarity of the different methods used for analysing the data.

With upcoming high-quality data from surveys such as the Extended Baryon Oscillation Spectroscopic Survey (eBOSS) or the Dark Energy Spectroscopic Instrument (DESI), improving the theoretical modeling and gaining a deeper understanding of the effects of neutrinos and dark radiation on structure formation at small scales are necessary, to obtain robust constraints free from systematic biases. Using a novel suite of hydrodynamical simulations that incorporate dark matter, baryons, massive neutrinos, and dark radiation, we present a detailed study of their impact on Lyman-Alpha forest observables. In particular, we accurately measure the tomographic evolution of the shape and amplitude of the small-scale matter and flux power spectra and search for unique signatures along with preferred scales where a neutrino mass detection may be feasible. We then investigate the thermal state of the intergalactic medium (IGM) through the temperature-density relation. Our findings suggest that at k~5h/Mpc the suppression on the matter power spectrum induced by M_nu=0.1 eV neutrinos can reach ~4% at z~3 when compared to a massless neutrino cosmology, and ~10% if a massless sterile neutrino is included; surprisingly, we also find good agreement (~2%) with some analytic predictions. For the 1D flux power spectrum, the highest response to free-streaming effects is achieved at k~0.005 s/km when M_nu=0.1 eV; this k-limit falls in the Lyman-Alpha forest regime, making the small-scale 1D flux power spectrum an excellent probe for detecting neutrino and dark radiation imprints. Our results indicate that the IGM at z~3 provides the best sensitivity to active and sterile neutrinos.

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.

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.

A flat Friedman-Roberson-Walker universe dominated by a cosmological constant ($\Lambda$) and cold dark matter (CDM) has been the working model preferred by cosmologists since the discovery of cosmic acceleration. However, tensions of various degrees of significance are known to be present among existing datasets within the $\Lambda$CDM framework. In particular, the Lyman-$\alpha$ forest measurement of the Baryon Acoustic Oscillations (BAO) by the Baryon Oscillation Spectroscopic Survey (BOSS) prefers a smaller value of the matter density fraction $\Omega_{\rm M}$ compared to the value preferred by cosmic microwave background (CMB). Also, the recently measured value of the Hubble constant, $H_0=73.24\pm1.74 \ {\rm km}\ {\rm s}^{-1} \ {\rm Mpc}^{-1}$, is $3.4\sigma$ higher than $66.93\pm0.62 \ {\rm km}\ {\rm s}^{-1} \ {\rm Mpc}^{-1}$ inferred from the Planck CMB data. In this work, we investigate if these tensions can be interpreted as evidence for a non-constant dynamical dark energy (DE). Using the Kullback-Leibler (KL) divergence to quantify the tension between datasets, we find that the tensions are relieved by an evolving DE, with the dynamical DE model preferred at a $3.5\sigma$ significance level based on the improvement in the fit alone. While, at present, the Bayesian evidence for the dynamical DE is insufficient to favour it over $\Lambda$CDM, we show that, if the current best fit DE happened to be the true model, it would be decisively detected by the upcoming DESI survey.

We present an analysis of the main systematic effects that could impact the measurement of CMB polarization with the proposed CORE space mission. We employ timeline-to-map simulations to verify that the CORE instrumental set-up and scanning strategy allow us to measure sky polarization to a level of accuracy adequate to the mission science goals. We also show how the CORE observations can be processed to mitigate the level of contamination by potentially worrying systematics, including intensity-to-polarization leakage due to bandpass mismatch, asymmetric main beams, pointing errors and correlated noise. We use analysis techniques that are well validated on data from current missions such as Planck to demonstrate how the residual contamination of the measurements by these effects can be brought to a level low enough not to hamper the scientific capability of the mission, nor significantly increase the overall error budget. We also present a prototype of the CORE photometric calibration pipeline, based on that used for Planck, and discuss its robustness to systematics, showing how CORE can achieve its calibration requirements. While a fine-grained assessment of the impact of systematics requires a level of knowledge of the system that can only be achieved in a future study phase, the analysis presented here strongly suggests that the main areas of concern for the CORE mission can be addressed using existing knowledge, techniques and algorithms.

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.

In this study, we probe the transition to cosmic homogeneity in the Large Scale Structure (LSS) of the Universe using the CMASS galaxy sample of BOSS spectroscopic survey which covers the largest effective volume to date, $3\ h^{-3}\ \mathrm{Gpc}^3$ at $0.43 \leq z \leq 0.7$. We study the scaled counts-in-spheres, $\mathcal{N}(<r)$, and the fractal correlation dimension, $\mathcal{D}_2(r)$, to assess the homogeneity scale of the universe using a $Landy\ \&\ Szalay$ inspired estimator. Defining the scale of transition to homogeneity as the scale at which $\mathcal{D}_2(r)$ reaches 3 within $1\%$, i.e. $\mathcal{D}_2(r)>2.97$ for $r>\mathcal{R}_H$, we find $\mathcal{R}_H = (63.3\pm0.7) \ h^{-1}\ \mathrm{Mpc}$, in agreement at the percentage level with the predictions of the $\Lambda$CDM model $\mathcal{R}_H=62.0\ h^{-1}\ \mathrm{Mpc}$. Thanks to the large cosmic depth of the survey, we investigate the redshift evolution of the transition to homogeneity scale and find agreement with the $\Lambda$CDM prediction. Finally, we find that $\mathcal{D}_2$ is compatible with $3$ at scales larger than $300\ h^{-1}\ $Mpc in all redshift bins. These results consolidate the Cosmological Principle and represent a precise consistency test of the $\Lambda CDM$ model.

We discuss the potential of a next generation space-borne Cosmic Microwave Background (CMB) experiment for studies of extragalactic sources. Our analysis has particular bearing on the definition of the future space project, CORE, that has been submitted in response to ESA's call for a Medium-size mission opportunity as the successor of the Planck satellite. Even though the effective telescope size will be somewhat smaller than that of Planck, CORE will have a considerably better angular resolution at its highest frequencies, since, in contrast with Planck, it will be diffraction limited at all frequencies. The improved resolution implies a considerable decrease of the source confusion, i.e. substantially fainter detection limits. In particular, CORE will detect thousands of strongly lensed high-z galaxies distributed over the full sky. The extreme brightness of these galaxies will make it possible to study them, via follow-up observations, in extraordinary detail. Also, the CORE resolution matches the typical sizes of high-z galaxy proto-clusters much better than the Planck resolution, resulting in a much higher detection efficiency; these objects will be caught in an evolutionary phase beyond the reach of surveys in other wavebands. Furthermore, CORE will provide unique information on the evolution of the star formation in virialized groups and clusters of galaxies up to the highest possible redshifts. Finally, thanks to its very high sensitivity, CORE will detect the polarized emission of thousands of radio sources and, for the first time, of dusty galaxies, at mm and sub-mm wavelengths, respectively.

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 perform a tomographic baryon acoustic oscillations analysis using the two-point galaxy correlation function measured from the combined sample of BOSS DR12, which covers the redshift range of $0.2<z<0.75$. Splitting the sample into multiple overlapping redshift slices to extract the redshift information of galaxy clustering, we obtain a measurement of $D_A(z)/r_d$ and $H(z)r_d$ at nine effective redshifts with the full covariance matrix calibrated using MultiDark-Patchy mock catalogues. Using the reconstructed galaxy catalogues, we obtain the precision of $1.3\%-2.2\%$ for $D_A(z)/r_d$ and $2.1\%-6.0\%$ for $H(z)r_d$. To quantify the gain from the tomographic information, we compare the constraints on the cosmological parameters using our 9-bin BAO measurements, the consensus 3-bin BAO and RSD measurements at three effective redshifts in \citet{Alam2016}, and the non-tomographic (1-bin) BAO measurement at a single effective redshift. Comparing the 9-bin with 1-bin constraint result, it can improve the dark energy Figure of Merit by a factor of 1.24 for the Chevallier-Polarski-Linder parametrisation for equation of state parameter $w_{\rm DE}$. The errors of $w_0$ and $w_a$ from 9-bin constraints are slightly improved when compared to the 3-bin constraint result.

We present a Bayesian phase-space reconstruction of the cosmic large-scale matter density and velocity fields from the SDSS-III Baryon Oscillations Spectroscopic Survey Data Release 12 (BOSS DR12) CMASS galaxy clustering catalogue. We rely on a given $\Lambda$CDM cosmology, a mesh resolution in the range of 6-10 $h^{-1}$ Mpc, and a lognormal-Poisson model with a redshift dependent nonlinear bias. The bias parameters are derived from the data and a general renormalised perturbation theory approach. We use combined Gibbs and Hamiltonian sampling, implemented in the \textsc{argo} code, to iteratively reconstruct the dark matter density field and the coherent peculiar velocities of individual galaxies, correcting hereby for coherent redshift space distortions (RSD). Our tests relying on accurate $N$-body based mock galaxy catalogues, show unbiased real space power spectra of the nonlinear density field up to $k\sim0.2\, h$ Mpc$^{-1}$, and vanishing quadrupoles down to $r\sim20\,h^{-1}$ Mpc. We also demonstrate that the nonlinear cosmic web can be obtained from the tidal field tensor based on the Gaussian component of the reconstructed density field. We find that the reconstructed velocities have a statistical correlation coefficient compared to the true velocities of each individual lightcone mock galaxy of $r\sim0.68$ including about 10% of satellite galaxies with virial motions (about $r=0.75$ without satellites). The power spectra of the velocity divergence agree well with theoretical predictions up to $k\sim0.2\,h\,{\rm Mpc}^{-1}$. This work will be especially useful to improve, e.g. BAO reconstructions, kinematic Sunyaev-Zeldovich (kSZ), integrated Sachs-Wolfe (ISW) measurements, or environmental studies.

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.

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.

We perform a tomographic baryon acoustic oscillations (BAO) analysis using the monopole, quadrupole and hexadecapole of the redshift-space galaxy power spectrum measured from the pre-reconstructed combined galaxy sample of the completed Sloan Digital Sky Survey (SDSS-III) Baryon Oscillation Spectroscopic Survey (BOSS) Data Release (DR)12 covering the redshift range of $0.20<z<0.75$. By allowing for overlap between neighbouring redshift slices, we successfully obtained the isotropic and anisotropic BAO distance measurements within nine redshift slices to a precision of $1.5\%-3.4\%$ for $D_V/r_d$, $1.8\% -4.2\%$ for $D_A/r_d$ and $3.7\% - 7.5\%$ for $H \ r_d$, depending on effective redshifts. We provide our BAO measurement of $D_A/r_d$ and $H \ r_d$ with the full covariance matrix, which can be used for cosmological implications. Our measurements are consistent with those presented in \citet{Acacia}, in which the BAO distances are measured at three effective redshifts. We constrain dark energy parameters using our measurements, and find an improvement of the Figure-of-Merit of dark energy in general due to the temporal BAO information resolved. This paper is part of a set that analyses the final galaxy clustering dataset from BOSS.

The BOSS quasar sample is used to study cosmic homogeneity with a 3D survey in the redshift range $2.2<z<2.8$. We measure the count-in-sphere, $N(<\! r)$, i.e. the average number of objects around a given object, and its logarithmic derivative, the fractal correlation dimension, $D_2(r)$. For a homogeneous distribution $N(<\! r) \propto r^3$ and $D_2(r)=3$. Due to the uncertainty on tracer density evolution, 3D surveys can only probe homogeneity up to a redshift dependence, i.e. they probe so-called "spatial isotropy". Our data demonstrate spatial isotropy of the quasar distribution in the redshift range $2.2<z<2.8$ in a model-independent way, independent of any FLRW fiducial cosmology, resulting in $3-\langle D_2 \rangle < 1.7 \times 10^{-3}$ (2 $\sigma$) over the range $250<r<1200 \, h^{-1}$Mpc for the quasar distribution. If we assume that quasars do not have a bias much less than unity, this implies spatial isotropy of the matter distribution on large scales. Then, combining with the Copernican principle, we finally get homogeneity of the matter distribution on large scales. Alternatively, using a flat $\Lambda$CDM fiducial cosmology with CMB-derived parameters, and measuring the quasar bias relative to this $\Lambda$CDM model, our data provide a consistency check of the model, in terms of how homogeneous the Universe is on different scales. $D_2(r)$ is found to be compatible with our $\Lambda$CDM model on the whole $10<r<1200 \, h^{-1}$Mpc range. For the matter distribution we obtain $3-\langle D_2 \rangle < 5 \times 10^{-5}$ (2 $\sigma$) over the range $250<r<1200 \, h^{-1}$Mpc, consistent with homogeneity on large scales.

We present the large-scale 3-point correlation function (3PCF) of the SDSS DR12 CMASS sample of $777,202$ Luminous Red Galaxies, the largest-ever sample used for a 3PCF or bispectrum measurement. We make the first high-significance ($4.5\sigma$) detection of Baryon Acoustic Oscillations (BAO) in the 3PCF. Using these acoustic features in the 3PCF as a standard ruler, we measure the distance to $z=0.57$ to $1.7\%$ precision (statistical plus systematic). We find $D_{\rm V}= 2024\pm29\;{\rm Mpc\;(stat)}\pm20\;{\rm Mpc\;(sys)}$ for our fiducial cosmology (consistent with Planck 2015) and bias model. This measurement extends the use of the BAO technique from the 2-point correlation function (2PCF) and power spectrum to the 3PCF and opens an avenue for deriving additional cosmological distance information from future large-scale structure redshift surveys such as DESI. Our measured distance scale from the 3PCF is fairly independent from that derived from the pre-reconstruction 2PCF and is equivalent to increasing the length of BOSS by roughly 10\%; reconstruction appears to lower the independence of the distance measurements. Fitting a model including tidal tensor bias yields a moderate significance ($2.6\sigma)$ detection of this bias with a value in agreement with the prediction from local Lagrangian biasing.

We search for a galaxy clustering bias due to a modulation of galaxy number with the baryon-dark matter relative velocity resulting from recombination-era physics. We find no detected signal and place the constraint $b_v < 0.01$ on the relative velocity bias for the CMASS galaxies. This bias is an important potential systematic of Baryon Acoustic Oscillation (BAO) method measurements of the cosmic distance scale using the 2-point clustering. Our limit on the relative velocity bias indicates a systematic shift of no more than $0.3\%$ rms in the distance scale inferred from the BAO feature in the BOSS 2-point clustering, well below the $1\%$ statistical error of this measurement. This constraint is the most stringent currently available and has important implications for the ability of upcoming large-scale structure surveys such as DESI to self-protect against the relative velocity as a possible systematic.

We develop a new methodology called double-probe analysis with the aim of minimizing informative priors in the estimation of cosmological parameters. We extract the dark-energy-model-independent cosmological constraints from the joint data sets of Baryon Oscillation Spectroscopic Survey (BOSS) galaxy sample and Planck cosmic microwave background (CMB) measurement. We measure the mean values and covariance matrix of $\{R$, $l_a$, $\Omega_b h^2$, $n_s$, $log(A_s)$, $\Omega_k$, $H(z)$, $D_A(z)$, $f(z)\sigma_8(z)\}$, which give an efficient summary of Planck data and 2-point statistics from BOSS galaxy sample, where $R=\sqrt{\Omega_m H_0^2}\,r(z_*)$, and $l_a=\pi r(z_*)/r_s(z_*)$, $z_*$ is the redshift at the last scattering surface, and $r(z_*)$ and $r_s(z_*)$ denote our comoving distance to $z_*$ and sound horizon at $z_*$ respectively. The advantage of this method is that we do not need to put informative priors on the cosmological parameters that galaxy clustering is not able to constrain well, i.e. $\Omega_b h^2$ and $n_s$. Using our double-probe results, we obtain $\Omega_m=0.304\pm0.009$, $H_0=68.2\pm0.7$, and $\sigma_8=0.806\pm0.014$ assuming $\Lambda$CDM; and $\Omega_k=0.002\pm0.003$ and $w=-1.00\pm0.07$ assuming o$w$CDM. The results show no tension with the flat $\Lambda$CDM cosmological paradigm. By comparing with the full-likelihood analyses with fixed dark energy models, we demonstrate that the double-probe method provides robust cosmological parameter constraints which can be conveniently used to study dark energy models. We extend our study to measure the sum of neutrino mass and obtain $\Sigma m_\nu<0.10/0.22$ (68\%/95\%) assuming $\Lambda$CDM and $\Sigma m_\nu<0.26/0.52$ (68\%/95\%) assuming $w$CDM. This paper is part of a set that analyses the final galaxy clustering dataset from BOSS.

We analyse the broad-range shape of the monopole and quadrupole correlation functions of the BOSS Data Release 12 (DR12) CMASS and LOWZ galaxy sample to obtain constraints on the Hubble expansion rate $H(z)$, the angular-diameter distance $D_A(z)$, the normalised growth rate $f(z)\sigma_8(z)$, and the physical matter density $\Omega_mh^2$. We adopt wide and flat priors on all model parameters in order to ensure the results are those of a `single-probe' galaxy clustering analysis. We also marginalise over three nuisance terms that account for potential observational systematics affecting the measured monopole. However, such Monte Carlo Markov Chain analysis is computationally expensive for advanced theoretical models, thus we develop a new methodology to speed up our analysis. We obtain $\{D_A(z)r_{s,fid}/r_s$Mpc, $H(z)r_s/r_{s,fid}$kms$^{-1}$Mpc$^{-1}$, $f(z)\sigma_8(z)$, $\Omega_m h^2\}$ = $\{956\pm28$ , $75.0\pm4.0$ , $0.397 \pm 0.073$, $0.143\pm0.017\}$ at $z=0.32$ and $\{1421\pm23$, $96.7\pm2.7$ , $0.497 \pm 0.058$, $0.137\pm0.015\}$ at $z=0.59$ where $r_s$ is the comoving sound horizon at the drag epoch and $r_{s,fid}=147.66$Mpc for the fiducial cosmology in this study. In addition, we divide the galaxy sample into four redshift bins to increase the sensitivity of redshift evolution. However, we do not find improvements in terms of constraining dark energy model parameters. Combining our measurements with Planck data, we obtain $\Omega_m=0.306\pm0.009$, $H_0=67.9\pm0.7$kms$^{-1}$Mpc$^{-1}$, and $\sigma_8=0.815\pm0.009$ assuming $\Lambda$CDM; $\Omega_k=0.000\pm0.003$ assuming oCDM; $w=-1.01\pm0.06$ assuming $w$CDM; and $w_0=-0.95\pm0.22$ and $w_a=-0.22\pm0.63$ assuming $w_0w_a$CDM. Our results show no tension with the flat $\Lambda$CDM cosmological paradigm. This paper is part of a set that analyses the final galaxy clustering dataset from BOSS.

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.

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

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.

The SDSS-IV/eBOSS has an extensive quasar program that combines several selection methods. Among these, the photometric variability technique provides highly uniform samples, unaffected by the redshift bias of traditional optical-color selections, when $z= 2.7 - 3.5$ quasars cross the stellar locus or when host galaxy light affects quasar colors at $z < 0.9$. Here, we present the variability selection of quasars in eBOSS, focusing on a specific program that led to a sample of 13,876 quasars to $g_{\rm dered}=22.5$ over a 94.5 deg$^2$ region in Stripe 82, an areal density 1.5 times higher than over the rest of the eBOSS footprint. We use these variability-selected data to provide a new measurement of the quasar luminosity function (QLF) in the redshift range $0.68<z<4.0$. Our sample is denser, reaches deeper than those used in previous studies of the QLF, and is among the largest ones. At the faint end, our QLF extends to $M_g(z\!=\!2)=-21.80$ at low redshift and to $M_g(z\!=\!2)=-26.20$ at $z\sim 4$. We fit the QLF using two independent double-power-law models with ten free parameters each. The first model is a pure luminosity-function evolution (PLE) with bright-end and faint-end slopes allowed to be different on either side of $z=2.2$. The other is a simple PLE at $z<2.2$, combined with a model that comprises both luminosity and density evolution (LEDE) at $z>2.2$. Both models are constrained to be continuous at $z=2.2$. They present a flattening of the bright-end slope at large redshift. The LEDE model indicates a reduction of the break density with increasing redshift, but the evolution of the break magnitude depends on the parameterization. The models are in excellent accord, predicting quasar counts that agree within 0.3\% (resp., 1.1\%) to $g<22.5$ (resp., $g<23$). The models are also in good agreement over the entire redshift range with models from previous studies.

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 describe mock data-sets generated to simulate the high-redshift quasar sample in Data Release 11 (DR11) of the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS). The mock spectra contain Ly{\alpha} forest correlations useful for studying the 3D correlation function including Baryon Acoustic Oscillations (BAO). They also include astrophysical effects such as quasar continuum diversity and high-density absorbers, instrumental effects such as noise and spectral resolution, as well as imperfections introduced by the SDSS pipeline treatment of the raw data. The Ly{\alpha} forest BAO analysis of the BOSS collaboration, described in Delubac et al. 2014, has used these mock data-sets to develop and cross-check analysis procedures prior to performing the BAO analysis on real data, and for continued systematic cross checks. Tests presented here show that the simulations reproduce sufficiently well important characteristics of real spectra. These mock data-sets will be made available together with the data at the time of the Data Release 11.

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 present joint constraints on the number of effective neutrino species N_eff and the sum of neutrino masses M_nu, based on a technique which exploits the full information contained in the one-dimensional Lyman-Alpha forest flux power spectrum, complemented by additional cosmological probes. In particular, we obtain N_eff=2.91(+0.21)(-0.22) (95% CL) and M_nu<0.15 eV (95% CL) when we combine BOSS Lyman-Alpha forest data with CMB (Planck+ACT+SPT+WMAP polarization) measurements, and N_eff=2.88(+0.20)(-0.20) (95% CL) and M_nu<0.14 eV (95% CL) when we further add baryon acoustic oscillations. Our results provide evidence for the Cosmic Neutrino Background from N_eff~3 (N_eff=0 is rejected at more than 14 sigma), and rule out the possibility of a sterile neutrino thermalized with active neutrinos (i.e., N_eff=4) - or more generally any decoupled relativistic relic with Delta N_eff ~ 1 - at a significance of over 5 sigma, the strongest bound to date, implying that there is no need for exotic neutrino physics in the concordance LCDM model.

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.

(Abridged) We detect the large-scale structure of Lya emission in the Universe at redshifts z=2-3.5 by measuring the cross-correlation of Lya surface brightness with quasars in SDSS/BOSS. We use a million spectra targeting Luminous Red Galaxies at z<0.8, after subtracting a best fit model galaxy spectrum from each one, as an estimate of the high-redshift Lya surface brightness. The quasar-Lya emission cross-correlation we detect has a shape consistent with a LambdaCDM model with Omega_M =0.30^+0.10-0.07. The predicted amplitude of this cross-correlation is proportional to the product of the mean Lya surface brightness, <mu_alpha>, the amplitude of mass fluctuations, and the quasar and Lya emission bias factors. Using known values, we infer <mu_alpha>(b_alpha/3) = (3.9 +/- 0.9) x 10^-21 erg/s cm^-2 A^-1 arcsec^-2, where b_alpha is the Lya emission bias factor. If the dominant sources of Lya emission are star forming galaxies, we infer rho_SFR = (0.28 +/- 0.07) (3/b_alpha) /yr/Mpc^3 at z=2-3.5. For b_alpha=3, this value is a factor of 21-35 above previous estimates from individually detected Lya emitters, although consistent with the total rho_SFR derived from dust-corrected, continuum UV surveys. 97% of the Lya emission in the Universe at these redshifts is therefore undetected in previous surveys of Lya emitters. Our measurement is much greater than seen from stacking analyses of faint halos surrounding previously detected Lya emitters, but we speculate that it arises from similar Lya halos surrounding all luminous star-forming galaxies. We also detect redshift space anisotropy of the quasar-Lya emission cross-correlation, finding evidence at the 3.0 sigma level that it is radially elongated, consistent with distortions caused by radiative-transfer effects (Zheng et al. (2011)). Our measurements represent the first application of the intensity mapping technique to optical observations.

We present the measurements of the redshift-space three-point correlation function (3PCF) for z~0.5 luminous red galaxies of the CMASS sample in the Sloan Digital Sky Survey-III Baryon Oscillation Spectroscopic Survey Data Release 11. The 3PCF measurements are interpreted within the halo occupation distribution framework using high-resolution N-body simulations, and the model successfully reproduces the 3PCF on scales larger than 1Mpc/h. As with the case for the redshift-space two-point correlation functions, we find that the redshift-space 3PCF measurements also favour the inclusion of galaxy velocity bias in the model. In particular, the central galaxy in a halo is on average in motion with respect to the core of the halo. We discuss the potential of the small-scale 3PCF to tighten the constraints on the relation between galaxies and dark matter haloes and on the phase-space distribution of galaxies.

We present constraints on the parameters of the $\Lambda$CDM cosmological model in the presence of massive neutrinos, using the one-dimensional Ly$\alpha$ forest power spectrum obtained with the Baryon Oscillation Spectroscopic Survey (BOSS) of the Sloan Digital Sky Survey (SDSS) by Palanque-Delabrouille et al. (2013), complemented by additional cosmological probes. The interpretation of the measured Ly$\alpha$ spectrum is done using a second-order Taylor expansion of the simulated power spectrum. BOSS Ly$\alpha$ data alone provide better bounds than previous Ly$\alpha$ results, but are still poorly constraining, especially for the sum of neutrino masses $\sum m_\nu$, for which we obtain an upper bound of 1.1~eV (95\% CL), including systematics for both data and simulations. Ly$\alpha$ constraints on $\Lambda$CDM parameters and neutrino masses are compatible with CMB bounds from the Planck collaboration. Interestingly, the combination of Ly$\alpha$ with CMB data reduces the uncertainties significantly, due to very different directions of degeneracy in parameter space, leading to the strongest cosmological bound to date on the total neutrino mass, $\sum m_\nu < 0.15$~eV at 95\% CL (with a best-fit in zero). Adding recent BAO results further tightens this constraint to $\sum m_\nu < 0.14$~eV at 95\% CL. This bound is nearly independent of the statistical approach used, and of the different combinations of CMB and BAO data sets considered in this paper in addition to Ly$\alpha$. Given the measured values of the two squared mass differences $\Delta m^2$, this result tends to favor the normal hierarchy scenario against the inverted hierarchy scenario for the masses of the active neutrino species.

We report a detection of the baryon acoustic oscillation (BAO) feature in the flux-correlation function of the Ly{\alpha} forest of high-redshift quasars with a statistical significance of five standard deviations. The study uses 137,562 quasars in the redshift range $2.1\le z \le 3.5$ from the Data Release 11 (DR11) of the Baryon Oscillation Spectroscopic Survey (BOSS) of SDSS-III. This sample contains three times the number of quasars used in previous studies. The measured position of the BAO peak determines the angular distance, $D_A(z=2.34)$ and expansion rate, $H(z=2.34)$, both on a scale set by the sound horizon at the drag epoch, $r_d$. We find $D_A/r_d=11.28\pm0.65(1\sigma)^{+2.8}_{-1.2}(2\sigma)$ and $D_H/r_d=9.18\pm0.28(1\sigma)\pm0.6(2\sigma)$ where $D_H=c/H$. The optimal combination, $\sim D_H^{0.7}D_A^{0.3}/r_d$ is determined with a precision of $\sim2\%$. For the value $r_d=147.4~{\rm Mpc}$, consistent with the CMB power spectrum measured by Planck, we find $D_A(z=2.34)=1662\pm96(1\sigma)~{\rm Mpc}$ and $H(z=2.34)=222\pm7(1\sigma)~{\rm km\,s^{-1}Mpc^{-1}}$. Tests with mock catalogs and variations of our analysis procedure have revealed no systematic uncertainties comparable to our statistical errors. Our results agree with the previously reported BAO measurement at the same redshift using the quasar-Ly{\alpha} forest cross-correlation. The auto-correlation and cross-correlation approaches are complementary because of the quite different impact of redshift-space distortion on the two measurements. The combined constraints from the two correlation functions imply values of $D_A/r_d$ and $D_H/r_d$ that are, respectively, 7% low and 7% high compared to the predictions of a flat $\Lambda$CDM cosmological model with the best-fit Planck parameters. With our estimated statistical errors, the significance of this discrepancy is $\approx 2.5\sigma$.

The initial shear field plays a central role in the formation of large-scale structures, and in shaping the geometry, morphology, and topology of the cosmic web. We discuss a recent theoretical framework for the shear tensor, termed the `peak/dip picture', which accounts for the fact that halos/voids may form from local extrema of the density field - rather than from random spatial positions; the standard Doroshkevich's formalism is generalized, to include correlations between the density Hessian and shear field at special points in space around which halos/voids may form. We then present the `peak/dip excursion-set-based' algorithm, along with its most recent applications - merging peaks theory with the standard excursion set approach.

The quest for the neutrino mass is a central goal in contemporary cosmology, subject to intense scrutiny, and among different large-scale structure tracers the Lyman-Alpha forest is re-emerging as a unique tool to probe the neutrino mass at high-redshift - through characteristic imprints on the transmitted Lyman-Alpha flux. A detailed modeling of the low-density regions of the intergalactic medium in presence of massive neutrinos on scale ranging from a few to hundreds of megaparsecs is required, if one wants to interpret state-of-the-art data from observations of quasar spectra. To this end, we provide a suite of hydrodynamical simulations made with Gagdet-3, spanning different volumes and having a range of resolutions and neutrino masses (from M_nu=0.1 to 0.8 eV, assuming 3 degenerate species), specifically designed to meet the requirements of the Baryon Acoustic Spectroscopic Survey (BOSS). We adopted a particle-type implementation of massive neutrinos, and chose cosmological parameters compatible with the latest Planck (2013) results. While the resolution requirements match the quality of the SDSS-III/BOSS data, our numerical simulations will also establish a useful theoretical ground for upcoming surveys such as SDSS-IV/eBOSS and DESI. In the very near future, data from leading spectroscopic surveys will allow measuring the absolute mass scale of neutrinos, and determining the exact nature of the neutrino mass hierarchy; hence, we expect that this modeling will become increasingly useful.

The signature left in quasar spectra by the presence of neutral hydrogen in the Universe allows one to constrain the sum of the neutrino masses with improved sensitivity, with respect to laboratory experiments, and may shed a new light on the neutrino mass hierarchy and on the absolute mass scale of neutrinos. Constraints on cosmological parameters and on the dark energy equation of state can also be derived, from a joint parameter estimation procedure. However, this requires a detailed modeling of the line-of-sight power spectrum of the transmitted flux in the Lyman-Alpha (LyA) forest on scales ranging from a few to hundreds of Mpcs, which in turns demands the inclusion and careful treatment of cosmological neutrinos. To this end, we present here a suite of state-of-the-art hydrodynamical simulations with cold dark matter, baryons and massive neutrinos, specifically targeted for modeling the low-density regions of the IGM as probed by the LyA forest at high-redshift. The simulations span volumes ranging from (25 Mpc/h)^3 to (100 Mpc/h)^3, and are made using either 3 X 192^3~21 millions or 3 X 768^3~1.4 billion particles. The resolution of the various runs can be further enhanced, so that we can reach the equivalent of 3 X 3072^3~87 billion particles in a (100 Mpc/h)^3 box size. The chosen cosmological parameters are compatible with the latest Planck (2013) results, although we also explore the effect of slight variations in the main cosmological and astrophysical parameters. We adopt a particle-type implementation of massive neutrinos, and consider three degenerate species having masses M_nu =0.1, 0.2, 0.3, 0.4 and 0.8 eV, respectively. We improve on previous studies in several ways, in particular with updated routines for IGM radiative cooling and heating processes, and initial conditions based on 2LPT rather than the Zeldovich approximation.

We present a suite of cosmological N-body simulations with cold dark matter and baryons aiming at modeling the low-density regions of the IGM as probed by the Lyman-$\alpha$ forests at high redshift. The simulations are designed to match the requirements imposed by the quality of BOSS and eBOSS data. They are made using either 2x768$^3$ or 2x192$^3$ particles, spanning volumes ranging from (25 Mpc.h$^{-1})^3$ for high-resolution simulations to (100 Mpc.h$^{-1})^3$ for large-volume ones. Using a splicing technique, the resolution is further enhanced to reach the equivalent of simulations with 2x3072$^3$= 58 billion particles in a (100 Mpc.h$^{-1}$)^3 box size, i.e. a mean mass per gas particle of 1.2x10$^5$M_sun.h$^{-1}$. We show that the resulting power spectrum is accurate at the 2% level over the full range from a few Mpc to several tens of Mpc. We explore the effect on the one-dimensional transmitted-flux power spectrum of 4 cosmological parameters ($n_s, \sigma_8, \Omega_m, H_0$) and 2 astrophysical parameters ($T_0, \gamma$) related to the heating rate of the IGM. By varying the input parameters around a central model chosen to be in agreement with the latest Planck results, we built a grid of simulations that allows the study of the impact on the flux power spectrum of these six relevant parameters. We improve upon previous studies by not only measuring the effect of each parameter individually, but also probing the impact of the simultaneous variation of each pair of parameters. We thus provide a full second-order expansion, including cross-terms, around our central model. We check the validity of the second-order expansion with independent simulations obtained either with different cosmological parameters or different seeds. Finally, a comparison to the one-dimensional Ly-$\alpha$ forest power spectrum obtained with BOSS shows an excellent agreement.

We measure the large-scale cross-correlation of quasars with the Lyman alpha forest absorption, using over 164,000 quasars from Data Release 11 of the SDSS-III Baryon Oscillation Spectroscopic Survey. We extend the previous study of roughly 60,000 quasars from Data Release 9 to larger separations, allowing a measurement of the Baryonic Acoustic Oscillation (BAO) scale along the line of sight $c/(H(z=2.36) ~ r_s) = 9.0 \pm 0.3$ and across the line of sight $D_A(z=2.36) / ~ r_s = 10.8 \pm 0.4$, consistent with CMB and other BAO data. Using the best fit value of the sound horizon from Planck data ($r_s=147.49 Mpc$), we can translate these results to a measurement of the Hubble parameter of $H(z=2.36) = 226 \pm 8 km/s / Mpc$ and of the angular diameter distance of $D_A(z=2.36) = 1590 \pm 60 Mpc$. The measured cross-correlation function and an update of the code to fit the BAO scale (baofit) are made publicly available.

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.

In the context of the study of the Integrated Sachs Wolfe effect (ISW), we construct a template of the projected density distribution up to $z\simeq 0.7$ by using the Luminous Galaxies (LGs) from the Sloan Digital Sky Survey DR8. We use a photo-z catalogue trained with more than a hundred thousand galaxies from BOSS in the SDSS DR8 imaging area. We consider two different LG samples whose selection matches that of SDSS-III/BOSS: the LOWZ sample ($z\in [0.15,0.5]$) and the CMASS sample ($z\in[0.4,0.7]$). When building the LG density maps we use the information from star density, survey footprint, seeing conditions, sky emission, dust extinction and airmass to explore the impact of these artifacts on the two LG samples. In agreement with previous studies, we find that the CMASS sample is particularly sensitive to Galactic stars, which dominate the contribution to the auto-angular power spectrum below $\ell=7$. Other potential systematics affect mostly the low multipole range ($\ell\in[2,7]$), but leave fluctuations on smaller scales practically unchanged. The resulting power spectra in the multipole range $\ell\in[2,100]$ for the LOWZ, CMASS and LOWZ+CMASS samples are compatible with linear $\Lambda$CDM expectations and constant bias values of $b=1.98 \pm 0.11$, $2.08\pm0.14$ and $1.88\pm 0.11$, respectively, with no traces of non-Gaussianity: $f_{\rm NL}^{\rm local}=59\pm 75$ at 95% confidence level for the full LOWZ+CMASS sample in the range $\ell\in[4,100]$. After cross-correlating WMAP-9yr data with the LOWZ+CMASS LG density field, the ISW signal is detected at the level of 1.62--1.69$\,\sigma$. While this result is in close agreement with predictions from Monte Carlo simulations in the concordance $\Lambda$CDM model, it cannot rule out by itself an Einstein-de Sitter scenario, and has a moderately low signal compared to previous studies conducted on subsets of this LG sample.

We have developed two independent methods to measure the one-dimensional power spectrum of the transmitted flux in the Lyman-$\alpha$ forest. The first method is based on a Fourier transform, and the second on a maximum likelihood estimator. The two methods are independent and have different systematic uncertainties. The determination of the noise level in the data spectra was subject to a novel treatment, because of its significant impact on the derived power spectrum. We applied the two methods to 13,821 quasar spectra from SDSS-III/BOSS DR9 selected from a larger sample of over 60,000 spectra on the basis of their high quality, large signal-to-noise ratio, and good spectral resolution. The power spectra measured using either approach are in good agreement over all twelve redshift bins from $<z> = 2.2$ to $<z> = 4.4$, and scales from 0.001 $\rm(km/s)^{-1}$ to $0.02 \rm(km/s)^{-1}$. We determine the methodological and instrumental systematic uncertainties of our measurements. We provide a preliminary cosmological interpretation of our measurements using available hydrodynamical simulations. The improvement in precision over previously published results from SDSS is a factor 2--3 for constraints on relevant cosmological parameters. For a $\Lambda$CDM model and using a constraint on $H_0$ that encompasses measurements based on the local distance ladder and on CMB anisotropies, we infer $\sigma_8 =0.83\pm0.03$ and $n_s= 0.97\pm0.02$ based on \ion{H}{i} absorption in the range $2.1<z<3.7$.

We present measurements of the genus topology of luminous red galaxies (LRGs) from the Sloan Digital Sky Survey (SDSS) Data Release 7 catalog, with unprecedented statistical significance. To estimate the uncertainties in the measured genus, we construct 81 mock SDSS LRG surveys along the past light cone from the Horizon Run 3, one of the largest N-body simulations to date that evolved 7210^3 particles in a 10815 Mpc/h size box. After carefully modeling and removing all known systematic effects due to finite pixel size, survey boundary, radial and angular selection functions, shot noise and galaxy biasing, we find the observed genus amplitude to reach 272 at 22 Mpc/h smoothing scale with an uncertainty of 4.2%; the estimated error fully incorporates cosmic variance. This is the most accurate constraint of the genus amplitude to date, which significantly improves on our previous results. In particular, the shape of the genus curve agrees very well with the mean topology of the SDSS LRG mock surveys in the LCDM universe. However, comparison with simulations also shows small deviations of the observed genus curve from the theoretical expectation for Gaussian initial conditions. While these discrepancies are mainly driven by known systematic effects such as those of shot noise and redshift-space distortions, they do contain important cosmological information on the physical effects connected with galaxy formation, gravitational evolution and primordial non-Gaussianity. We address here the key role played by systematics on the genus curve, and show how to accurately correct for their effects to recover the topology of the underlying matter. In a forthcoming paper, we provide an interpretation of those deviations in the context of the local model of non-Gaussianity.

The Lyman-$\beta$ forest refers to a region in the spectra of distant quasars that lies between the rest-frame Lyman-$\beta$ and Lyman-$\gamma$ emissions. The forest in this region is dominated by a combination of absorption due to resonant Ly$\alpha$ and Ly$\beta$ scattering. When considering the 1D Ly$\beta$ forest in addition to the 1D Ly$\alpha$ forest, the full statistical description of the data requires four 1D power spectra: Ly$\alpha$ and Ly$\beta$ auto-power spectra and the Ly$\alpha$-Ly$\beta$ real and imaginary cross-power spectra. We describe how these can be measured using an optimal quadratic estimator that naturally disentangles Ly$\alpha$ and Ly$\beta$ contributions. Using a sample of approximately 60,000 quasar sight-lines from the BOSS Data Release 9, we make the measurement of the one-dimensional power spectrum of fluctuations due to the Ly$\beta$ resonant scattering. While we have not corrected our measurements for resolution damping of the power and other systematic effects carefully enough to use them for cosmological constraints, we can robustly conclude the following: i) Ly$\beta$ power spectrum and Ly$\alpha$-Ly$\beta$ cross spectra are detected with high statistical significance; ii) the cross-correlation coefficient is $\approx 1$ on large scales; iii) the Ly$\beta$ measurements are contaminated by the associated OVI absorption, which is analogous to the SiIII contamination of the Ly$\alpha$ forest. Measurements of the Ly$\beta$ forest will allow extension of the usable path-length for the Ly$\alpha$ measurements while allowing a better understanding of the physics of intergalactic medium and thus more robust cosmological constraints.

We measure the luminosity and color dependence and the redshift evolution of galaxy clustering in the Sloan Digital Sky Survey-III Baryon Oscillation Spectroscopic Survey Ninth Data Release. We focus on the projected two-point correlation function (2PCF) of subsets of its CMASS sample, which includes about 260,000 galaxies over ~3,300 sq. deg in the redshift range 0.43<z<0.7. To minimize the selection effect on galaxy clustering, we construct well-defined luminosity and color subsamples by carefully accounting for the CMASS galaxy selection cuts. The 2PCF of the whole CMASS sample, if approximated by a power-law, has a correlation length of r_0=7.93\pm0.06Mpc/h and an index of \gamma=1.85\pm0.01. Clear dependences on galaxy luminosity and color are found for the projected 2PCF in all redshift bins, with more luminous and redder galaxies generally exhibiting stronger clustering and steeper 2PCF. The color dependence is also clearly seen for galaxies within the red sequence, consistent with the behavior of SDSS-II main sample galaxies at lower redshifts. At a given luminosity (k+e corrected), no significant evolution of the projected 2PCFs with redshift is detected for red sequence galaxies. We also construct galaxy samples of fixed number density at different redshifts, using redshift-dependent magnitude thresholds. The clustering of these galaxies in the CMASS redshift range is found to be consistent with that predicted by passive evolution. Our measurements of the luminosity and color dependence and redshift evolution of galaxy clustering will allow for detailed modeling of the relation between galaxies and dark matter halos and new constraints on galaxy formation and evolution.

We use the Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 9 (DR9) to detect and measure the position of the Baryonic Acoustic Oscillation (BAO) feature in the three-dimensional correlation function in the Lyman-alpha forest flux fluctuations at a redshift z=2.4. The feature is clearly detected at significance between 3 and 5 sigma (depending on the broadband model and method of error covariance matrix estimation) and is consistent with predictions of the standard LCDM model. We assess the biases in our method, stability of the error covariance matrix and possible systematic effects. We fit the resulting correlation function with several models that decouple the broadband and acoustic scale information. For an isotropic dilation factor, we measure 100x(alpha_iso-1) = -1.6 ^{+2.0+4.3+7.4}_{-2.0-4.1-6.8} (stat.) +/- 1.0 (syst.) (multiple statistical errors denote 1,2 and 3 sigma confidence limits) with respect to the acoustic scale in the fiducial cosmological model (flat LCDM with Omega_m=0.27, h=0.7). When fitting separately for the radial and transversal dilation factors we find marginalised constraints 100x(alpha_par-1) = -1.3 ^{+3.5+7.6 +12.3}_{-3.3-6.7-10.2} (stat.) +/- 2.0 (syst.) and 100x(alpha_perp-1) = -2.2 ^{+7.4+17}_{-7.1-15} +/- 3.0 (syst.). The dilation factor measurements are significantly correlated with cross-correlation coefficient of ~ -0.55. Errors become significantly non-Gaussian for deviations over 3 standard deviations from best fit value. Because of the data cuts and analysis method, these measurements give tighter constraints than a previous BAO analysis of the BOSS DR9 Lyman-alpha forest sample, providing an important consistency test of the standard cosmological model in a new redshift regime.

We present a new algorithm to sample the constrained eigenvalues of the initial shear field associated with Gaussian statistics, called the `peak/dip excursion-set-based' algorithm, at positions which correspond to peaks or dips of the correlated density field. The computational procedure is based on a new formula which extends Doroshkevich's unconditional distribution for the eigenvalues of the linear tidal field, to account for the fact that halos and voids may correspond to maxima or minima of the density field. The ability to differentiate between random positions and special points in space around which halos or voids may form (peaks/dips), encoded in the new formula and reflected in the algorithm, naturally leads to a straightforward implementation of an excursion set model for peaks and dips in Gaussian random fields - one of the key advantages of this sampling procedure. In addition, it offers novel insights into the statistical description of the cosmic web. As a first physical application, we show how the standard distributions of shear ellipticity and prolateness in triaxial models of structure formation are modified by the constraint. In particular, we provide a new expression for the conditional distribution of shape parameters given the density peak constraint, which generalizes some previous literature work. The formula has important implications for the modeling of non-spherical dark matter halo shapes, in relation to their initial shape distribution. We also test and confirm our theoretical predictions for the individual distributions of eigenvalues subjected to the extremum constraint, along with other directly related conditional probabilities. Finally, we indicate how the proposed sampling procedure naturally integrates into the standard excursion set model, potentially solving some of its well-known problems, and into the ellipsoidal collapse framework. (abridged)

We report a detection of the baryon acoustic oscillation (BAO) feature in the three-dimensional correlation function of the transmitted flux fraction in the \Lya forest of high-redshift quasars. The study uses 48,640 quasars in the redshift range $2.1\le z \le 3.5$ from the Baryon Oscillation Spectroscopic Survey (BOSS) of the third generation of the Sloan Digital Sky Survey (SDSS-III). At a mean redshift $z=2.3$, we measure the monopole and quadrupole components of the correlation function for separations in the range $20\hMpc<r<200\hMpc$. A peak in the correlation function is seen at a separation equal to $(1.01\pm0.03)$ times the distance expected for the BAO peak within a concordance $\Lambda$CDM cosmology. This first detection of the BAO peak at high redshift, when the universe was strongly matter dominated, results in constraints on the angular diameter distance $\da$ and the expansion rate $H$ at $z=2.3$ that, combined with priors on $H_0$ and the baryon density, require the existence of dark energy. Combined with constraints derived from Cosmic Microwave Background (CMB) observations, this result implies $H(z=2.3)=(224\pm8){\rm km\,s^{-1}Mpc^{-1}}$, indicating that the time derivative of the cosmological scale parameter $\dot{a}=H(z=2.3)/(1+z)$ is significantly greater than that measured with BAO at $z\sim0.5$. This demonstrates that the expansion was decelerating in the range $0.7<z<2.3$, as expected from the matter domination during this epoch. Combined with measurements of $H_0$, one sees the pattern of deceleration followed by acceleration characteristic of a dark-energy dominated universe.

We describe fitting methods developed to analyze fluctuations in the Lyman-{\alpha} forest and measure the parameters of baryon acoustic oscillations (BAO). We apply our methods to BOSS Data Release 9. Our method is based on models of the three-dimensional correlation function in physical coordinate space, and includes the effects of redshift-space distortions, anisotropic non-linear broadening, and broadband distortions. We allow for independent scale factors along and perpendicular to the line of sight to minimize the dependence on our assumed fiducial cosmology and to obtain separate measurements of the BAO angular and relative velocity scales. Our fitting software and the input files needed to reproduce our main BOSS Data Release 9 results are publicly available.

We present the BOSS Lyman-alpha (Lya) Forest Sample from SDSS Data Release 9, comprising 54,468 quasar spectra with zqso > 2.15 suitable for Lya forest analysis. This data set probes the intergalactic medium with absorption redshifts 2.0 < z_alpha < 5.7 over an area of 3275 square degrees, and encompasses an approximate comoving volume of 20 h^-3 Gpc^3. With each spectrum, we have included several products designed to aid in Lya forest analysis: improved sky masks that flag pixels where data may be unreliable, corrections for known biases in the pipeline estimated noise, masks for the cores of damped Lya systems and corrections for their wings, and estimates of the unabsorbed continua so that the observed flux can be converted to a fractional transmission. The continua are derived using a principal component fit to the quasar spectrum redwards of restframe Lya (lambda > 1216 Ang), extrapolated into the forest region and normalized by a linear function to fit the expected evolution of the Lya forest mean-flux. The estimated continuum errors are ~5% rms. We also discuss possible systematics arising from uncertain spectrophotometry and artifacts in the flux calibration; global corrections for the latter are provided. Our sample provides a convenient starting point for users to analyze clustering in BOSS Lya forest data, and it provides a fiducial data set that can be used to compare results from different analyses of baryon acoustic oscillations in the Lya forest. The full data set is available from the SDSS-III DR9 web site.

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.

A budget neutral strategy is proposed for NSF to lead the implementation of multimessenger astronomy and astrophysics, as outlined in the Astro2010 Decadal Survey. The emerging capabilities for simultaneous measurements of physical and astronomical data through the different windows of electromagnetic, hadronic and gravitational radiation processes call for a vigorous pursuit of new synergies. The proposed approach is aimed at the formation of new collaborations and multimessenger data-analysis, to transcend the scientific inquiries made within a single window of observations. In view of budgetary constraints, we propose to include the multimessenger dimension in the ranking of proposals submitted under existing NSF programs.

The initial shear field, characterized by a primordial perturbation potential, plays a crucial role in the formation of large scale structures. Hence, considerable analytic work has been based on the joint distribution of its eigenvalues, associated with Gaussian statistics. In addition, directly related morphological quantities such as ellipticity or prolateness are essential tools in understanding the formation and structural properties of halos, voids, sheets and filaments, their relation with the local environment, and the geometrical and dynamical classification of the cosmic web. To date, most analytic work has been focused on Doroshkevich's unconditional formulae for the eigenvalues of the linear tidal field, which neglect the fact that halos (voids) may correspond to maxima (minima) of the density field. I present here new formulae for the constrained eigenvalues of the initial shear field associated with Gaussian statistics, which include the fact that those eigenvalues are related to regions where the source of the displacement is positive (negative): this is achieved by requiring the Hessian matrix of the displacement field to be positive (negative) definite. The new conditional formulae naturally reduce to Doroshkevich's unconditional relations, in the limit of no correlation between the potential and the density fields. As a direct application, I derive the individual conditional distributions of eigenvalues and point out the connection with previous literature. Finally, I outline other possible theoretically- or observationally-oriented uses, ranging from studies of halo and void triaxial formation, development of structure-finding algorithms for the morphology and topology of the cosmic web, till an accurate mapping of the gravitational potential environment of galaxies from current and future generation galaxy redshift surveys.

We present two large cosmological N-body simulations, called Horizon Run 2 (HR2) and Horizon Run 3 (HR3), made using 6000^3 = 216 billions and 7210^3 = 374 billion particles, spanning a volume of (7.200 Gpc/h)^3 and (10.815 Gpc/h)^3, respectively. These simulations improve on our previous Horizon Run 1 (HR1) up to a factor of 4.4 in volume, and range from 2600 to over 8800 times the volume of the Millennium Run. In addition, they achieve a considerably finer mass resolution, down to 1.25x10^11 M_sun/h, allowing to resolve galaxy-size halos with mean particle separations of 1.2 Mpc/h and 1.5 Mpc/h, respectively. We have measured the power spectrum, correlation function, mass function and basic halo properties with percent level accuracy, and verified that they correctly reproduce the LCDM theoretical expectations, in excellent agreement with linear perturbation theory. Our unprecedentedly large-volume N-body simulations can be used for a variety of studies in cosmology and astrophysics, ranging from large-scale structure topology, baryon acoustic oscillations, dark energy and the characterization of the expansion history of the Universe, till galaxy formation science - in connection with the new SDSS-III. To this end, we made a total of 35 all-sky mock surveys along the past light cone out to z=0.7 (8 from the HR2 and 27 from the HR3), to simulate the BOSS geometry. The simulations and mock surveys are already publicly available at http://astro.kias.re.kr/Horizon-Run23/.

We combine the physics of the ellipsoidal collapse model with the excursion set theory to study the shapes of dark matter halos. In particular, we develop an analytic approximation to the nonlinear evolution that is more accurate than the Zeldovich approximation; we introduce a planar representation of halo axis ratios, which allows a concise and intuitive description of the dynamics of collapsing regions and allows one to relate the final shape of a halo to its initial shape; we provide simple physical explanations for some empirical fitting formulae obtained from numerical studies. Comparison with simulations is challenging, as there is no agreement about how to define a non-spherical gravitationally bound object. Nevertheless, we find that our model matches the conditional minor-to-intermediate axis ratio distribution rather well, although it disagrees with the numerical results in reproducing the minor-to-major axis ratio distribution. In particular, the mass dependence of the minor-to-major axis distribution appears to be the opposite to what is found in many previous numerical studies, where low-mass halos are preferentially more spherical than high-mass halos. In our model, the high-mass halos are predicted to be more spherical, consistent with results based on a more recent and elaborate halo finding algorithm, and with observations of the mass dependence of the shapes of early-type galaxies. We suggest that some of the disagreement with some previous numerical studies may be alleviated if we consider only isolated halos.

We use the statistics of regions above or below a temperature threshold (excursion sets) to study the cosmic microwave background (CMB) anisotropy in models with primordial non-Gaussianity of the local type. By computing the full-sky spatial distribution and clustering of pixels above/below threshold from a large set of simulated maps with different levels of non-Gaussianity, we find that a positive value of the dimensionless non-linearity parameter f_NL enhances the number density of the cold CMB excursion sets along with their clustering strength, and reduces that of the hot ones. We quantify the robustness of this effect, which may be important to discriminate between the simpler Gaussian hypothesis and non-Gaussian scenarios, arising either from non-standard inflation or alternative early-universe models. The clustering of hot and cold pixels exhibits distinct non-Gaussian signatures, particularly at angular scales of about 75 arcmin (i.e. around the Doppler peak), which increase linearly with f_NL. Moreover, the clustering changes strongly as a function of the smoothing angle. We propose several statistical tests to maximize the detection of a local primordial non-Gaussian signal, and provide some theoretical insights within this framework, including an optimal selection of the threshold level. We also describe a procedure which aims at minimizing the cosmic variance effect, the main limit within this statistical framework.

In addition to the maximum likelihood approach, there are two other methods which are commonly used to reconstruct the true redshift distribution from photometric redshift datasets: one uses a deconvolution method, and the other a convolution. We show how these two techniques are related, and how this relationship can be extended to include the study of galaxy scaling relations in photometric datasets. We then show what additional information photometric redshift algorithms must output so that they too can be used to study galaxy scaling relations, rather than just redshift distributions. We also argue that the convolution based approach may permit a more efficient selection of the objects for which calibration spectra are required.

Noisy distance estimates associated with photometric rather than spectroscopic redshifts lead to a mis-estimate of the luminosities, and produce a correlated mis-estimate of the sizes. We consider a sample of early-type galaxies from the SDSS DR6 for which both spectroscopic and photometric information is available, and apply the generalization of the V_max method to correct for these biases. We show that our technique recovers the true redshift, magnitude and size distributions, as well as the true size-luminosity relation. We find that using only 10% of the spectroscopic information randomly spaced in our catalog is sufficient for the reconstructions to be accurate within about 3%, when the photometric redshift error is dz = 0.038. We then address the problem of extending our method to deep redshift catalogs, where only photometric information is available. In addition to the specific applications outlined here, our technique impacts a broader range of studies, when at least one distance-dependent quantity is involved. It is particularly relevant for the next generation of surveys, some of which will only have photometric information.

We present measurements of the clustering of hot and cold patches in the microwave background sky as measured from the Wilkinson Microwave Anisotropy Probe (WMAP) five-year data. These measurements are compared with theoretical predictions which assume that the cosmological signal obeys Gaussian statistics. We find significant differences from the simplest Gaussian-based prediction. However, the measurements are sensitive to the fact that the noise is spatially inhomogeneous (e.g., because different parts of the sky were observed for different lengths of time). We show how to account for this spatial inhomogeneity when making predictions. Differences from the Gaussian-based expectation remain even after this more careful accounting of the noise. In particular, we note that hot and cold pixels cluster differently within the same temperature thresholds at few-degree scales. While these findings may indicate primordial non-Gaussianity, we discuss other plausible explanations for these discrepancies. In addition, we find some deviations from Gaussianity at sub-degree scales, especially in the W band, whose origin may be associated with extragalactic dust emission.

Noisy distance estimates associated with photometric rather than spectroscopic redshifts lead to a mis-estimate of the luminosities, and produce a correlated mis-estimate of the sizes. We consider a sample of early-type galaxies from the SDSS DR6 and apply the generalization of the V_max method to correct for these biases. We show that our technique recovers the true redshift, magnitude and size distributions, as well as the true size-luminosity relation. Regardless the specific application outlined here, our method impacts a broader range of studies, when at least one distance-dependent quantity is involved.

Many physical properties of galaxies correlate with one another, and these correlations are often used to constrain galaxy formation models. Such correlations include the color-magnitude relation, the luminosity-size relation, the Fundamental Plane, etc. However, the transformation from observable (e.g. angular size, apparent brightness) to physical quantity (physical size, luminosity), is often distance-dependent. Noise in the distance estimate will lead to biased estimates of these correlations, thus compromising the ability of photometric redshift surveys to constrain galaxy formation models. We describe two methods which can remove this bias. One is a generalization of the V_max method, and the other is a maximum likelihood approach. We illustrate their effectiveness by studying the size-luminosity relation in a mock catalog, although both methods can be applied to other scaling relations as well. We show that if one simply uses photometric redshifts one obtains a biased relation; our methods correct for this bias and recover the true relation.