
The fourth generation of the Sloan Digital Sky Survey (SDSSIV) began
observations in July 2014. It pursues three core programs: APOGEE2, MaNGA, and
eBOSS. In addition, eBOSS contains two major subprograms: TDSS and SPIDERS.
This paper describes the first data release from SDSSIV, 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 variabilityselected
objects from TDSS and Xray selected objects from SPIDERS. DR13 includes new
reductions of the SDSSIII BOSS data, improving the spectrophotometric
calibration and redshift classification. DR13 releases new reductions of the
APOGEE1 data from SDSSIII, 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 valueadded catalogs are being released in tandem with
DR13, in particular target catalogs relevant for eBOSS, TDSS, and SPIDERS, and
an updated redclump 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 ~6year operations of
SDSSIV.

If a significant fraction of the dark matter in the Universe is made of an
ultralight 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 cutoff 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 nonlinear 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 wellstudied 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 higherresolution 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
"smallscale problems" of $\Lambda$CDM.

We use the BOSS DR9 quasar spectra to constrain nonthermal dark matter
models: coldpluswarm 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 nonlinear 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 higherresolution data (XQ100,
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
higherresolution bounds.

The freestreaming of keVscale particles impacts structure growth on scales
that are probed by the Lymanalpha forest of distant quasars. Using an
unprecedentedly large sample of mediumresolution QSO spectra from the ninth
data release of SDSS, along with a stateoftheart set of hydrodynamical
simulations to model the Lymanalpha forest in the nonlinear 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 nonresonantly produced righthanded neutrino.
This limit depends on the value on $n_s$, and Planck measures a higher value of
$n_s$ than SDSSIII/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 SDSSIII
data featuring smaller uncertainties and covering a larger redshift range than
SDSSI data, our bounds confirm the most stringent results established by
previous works and are further at odds with a purely nonresonantly produced
sterile neutrino as dark matter.

We consider a nonstandard 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 matterantimatter
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 matterantimatter 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 onedimensional Ly$\alpha$forest power spectrum measured by
PalanqueDelabrouille et al. (2013) from SDSSIII/BOSS, complemented by Planck
2015 cosmic microwave background (CMB) data and other cosmological probes. This
paper improves on the previous analysis by PalanqueDelabrouille 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 tensortoscalar 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 slowroll
inflation.