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.

If a significant fraction of the dark matter in the Universe is made of an ultra-light scalar field, named fuzzy dark matter (FDM) with a mass $m_a$ of the order of $10^{-22}-10^{-21}$ eV, then its de Broglie wavelength is large enough to impact the physics of large scale structure formation. In particular, the associated cut-off in the linear matter power spectrum modifies the structure of the intergalactic medium (IGM) at the scales probed by the Lyman-$\alpha$ forest of distant quasars. We study this effect by making use of dedicated cosmological simulations which take into account the hydrodynamics of the IGM. We explore heuristically the amplitude of quantum pressure for the FDM masses considered here and conclude that quantum effects should not modify significantly the non-linear evolution of matter density at the scales relevant to the measured Lyman-$\alpha$ flux power, and for $m_a \geq 10^{-22}$ eV. We derive a scaling law between $m_a$ and the mass of the well-studied thermal warm dark matter (WDM) model that is best adapted to the Lyman-$\alpha$ forest data, and differs significantly from the one infered by a simple linear extrapolation. By comparing FDM simulations with the Lyman-$\alpha$ flux power spectra determined from the BOSS survey, and marginalizing over relevant nuisance parameters, we exclude FDM masses in the range $10^{-22} \leq m_a < 2.3\times 10^{-21}$ eV at 95 % CL. Adding higher-resolution Lyman-$\alpha$ spectra extends the exclusion range up to $2.9\times 10^{-21}$ eV. This provides a significant constraint on FDM models tailored to solve the "small-scale problems" of $\Lambda$CDM.

We use the BOSS DR9 quasar spectra to constrain non-thermal dark matter models: cold-plus-warm dark matter and sterile neutrinos resonantly produced in the presence of a lepton asymmetry. We establish constraints on the warm thermal relic mass $m_{x}$ and its relative abundance $F_{wdm}=\Omega_{wdm}/\Omega_{dm}$ using a suite of cosmological hydrodynamical simulations in 28 C+WDM configurations. The 2D bounds in the $m_{x} - F_{wdm}$ parameter space approximately follow $F_{\rm{wdm}} \sim 0.243 (\rm{keV}/m_x)^{-1.31}$. At 95% C.L., our limits from BOSS data alone imply that $m_{x}> 2.5$ keV if $F_{\rm{wdm}}>80\%$, while masses as low as 0.7 keV are consistent with the data if $F_{\rm{wdm}}<15\%$ of the total dark matter density. We also constrain sterile neutrino mass and mixing angle by further producing the non-linear flux power spectrum of 8 RPSN models, where the input linear power spectrum is computed directly from the particles distribution functions. We find values of lepton asymmetries for which sterile neutrinos as light as 6 keV (resp. 3.5 keV) are consistent with BOSS at $2\sigma$ (resp. $3\sigma$). These limits tighten by close to a factor of 2 for values of lepton asymmetries departing from those yielding the coolest distribution functions. Our bounds can be strengthened if we include higher-resolution data (XQ-100, HIRES MIKE). At these scales, however, the flux power spectrum exhibits a suppression compatible with the one expected from WDM. We do not investigate the mechanism responsible for this suppression, but we show how much our bounds would strengthen under the assumption that it is caused by the evolution of the IGM temperature. A 7 keV neutrino produced in a lepton asymmetry $\mathcal{L} = |n_{\nu_e} - n_{\bar{\nu}_e}|/s = 8 \times 10^{-6}$ is consistent at $1.9\sigma$ with BOSS data, and not necessarily excluded by our higher-resolution bounds.

The free-streaming of keV-scale particles impacts structure growth on scales that are probed by the Lyman-alpha forest of distant quasars. Using an unprecedentedly large sample of medium-resolution QSO spectra from the ninth data release of SDSS, along with a state-of-the-art set of hydrodynamical simulations to model the Lyman-alpha forest in the non-linear regime, we issue one of the tightest bounds to date, from Ly-$\alpha$ data alone, on pure dark matter particles : $m_X > 4.09 \: \rm{keV}$ (95% CL) for early decoupled thermal relics such as a hypothetical gravitino, and correspondingly $m_s > 24.4 \: \rm{keV}$ (95% CL) for a non-resonantly produced right-handed neutrino. This limit depends on the value on $n_s$, and Planck measures a higher value of $n_s$ than SDSS-III/BOSS. Our bounds thus change slightly when Ly-$\alpha$ data are combined with CMB data from Planck 2016. The limits shift to $m_X > 2.96 \: \rm{keV}$ (95% CL) and $m_s > 16.0 \: \rm{keV}$ (95% CL). Thanks to SDSS-III data featuring smaller uncertainties and covering a larger redshift range than SDSS-I data, our bounds confirm the most stringent results established by previous works and are further at odds with a purely non-resonantly produced sterile neutrino as dark matter.

We consider a non-standard cosmological model in which the universe contains as much matter as antimatter on large scales and presents a local baryon asymmetry. A key ingredient in our approach is that the baryon density distribution follows Gaussian fluctuations around a null value $\eta = 0$. Spatial domains featuring a positive (resp. negative) baryonic density value constitute regions dominated by matter (resp. antimatter). At the domains' annihilation interface, the typical density is going smoothly to zero, rather than following an abrupt step as assumed in previous symetric matter-antimatter models. As a consequence, the Cosmic Diffuse Gamma Background produced by annihilation is drastically reduced, allowing to easily pass COMPTEL's measurements limits. Similarly the Compton $y$ distorsion and CMB 'ribbons' are lowered by an appreciable factor. Therefore this model essentially escape previous constrainst on symetric matter-antimatter models. However, we produce an estimation of the CMB temperature fluctuations that would result from this model and confront it to data acquired from the Planck satellite. We construct a angular power spectrum in $\delta T / T_{CMB}$ assuming is can be approximated as an average of $C_\ell$ over a Gaussian distribution of $\Omega_B$ using Lewis & Challinor's CAMB software. The resulting $C_\ell$ are qualitatively satisfying. We quantify the goodness of fit using a simple $\chi^2$ test. We consider two distinct scenarios in which the fluctuations on $\Omega_B$ are compensated by fluctuations on $\Omega_{CDM}$ to assure a spatially flat $\Omega_\kappa = 0$ universe or not. In both cases, out best fit have $\Delta \chi^2 \gtrsim 2400$ (with respect to a fiducial $\Lambda$CDM model), empirically excluding our model by several tens of standard deviations.

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.