
The fourth generation of the Sloan Digital Sky Survey (SDSSIV) 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 SDSSIV in its first
two years of operation (July 20142016). 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 (APOGEE2),
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 SDSSIV 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. SDSSIV is
planning to continue to collect astronomical data until 2020, and will be
followed by SDSSV.

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

We present a measurement of the anisotropic and isotropic Baryon Acoustic
Oscillations (BAO) from the extended Baryon Oscillation Spectroscopic Survey
Data Release 14 quasar sample with optimal redshift weights. Applying the
redshift weights improves the constraint on the BAO dilation parameter
$\alpha(z_{\rm eff})$ by 17\%. We reconstruct the evolution history of the BAO
distance indicators in the redshift range of $0.8<z<2.2$. This paper is part of
a set that analyses the eBOSS DR14 quasar sample.

We measure the growth rate and its evolution using the anisotropic clustering
of the extended Baryon Oscillation Spectroscopic Survey (eBOSS) Data Release 14
(DR14) quasar sample, which includes $148\,659$ quasars covering the wide
redshift range of $0.8 < z < 2.2$ and a sky area of $2112.90$ $\rm deg^2$. To
optimise measurements we deploy a redshiftdependent weighting scheme, which
allows us to avoid binning, and perform the data analysis consistently
including the redshift evolution across the sample. We perform the analysis in
Fourier space, and use the redshift evolving power spectrum multipoles to
measure the redshift space distortion parameter $f\sigma_8$ and parameters
controlling the anisotropic projection of the cosmological perturbations. We
measure $f \sigma_8(z=1.52)=0.43 \pm 0.05 $ and $df\sigma_8/dz (z=1.52)=  0.16
\pm 0.08$, consistent with the expectation for a $\Lambda$CDM cosmology as
constrained by the Planck experiment.

We present an anisotropic analysis of Baryon Acoustic Oscillation (BAO)
signal from the SDSSIV 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
distanceredshift relation, relative to a fiducial model, as a quadratic
expansion. The coefficients of this expansion are used to reconstruct the
distanceredshift 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 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 fullshape clustering
measurements with the eBOSS DR14 quasar sample.

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 twopoint
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 bestfitting standard ${\rm \Lambda CDM}$
model to the latest cosmic microwave background data from Planck.

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

The main effect of massive neutrinos on the largescale structure consists in
a few percent suppression of matter perturbations on all scales below their
freestreaming scale. Such effect is of particular importance as it allows to
constraint the value of the sum of neutrino masses from measurements of the
galaxy power spectrum. In this work, we present the first measurements of the
next higherorder correlation function, the bispectrum, from Nbody simulations
that include massive neutrinos as particles. This is the simplest statistics
characterising the nonGaussian properties of the matter and dark matter halos
distributions. We investigate, in the first place, the suppression due to
massive neutrinos on the matter bispectrum, comparing our measurements with the
simplest perturbation theory predictions, finding the approximation of
neutrinos contributing at quadratic order in perturbation theory to provide a
good fit to the measurements in the simulations. On the other hand, as
expected, a linear approximation for neutrino perturbations would lead to
O($f_{\nu}$) errors on the total matter bispectrum at large scales. We then
attempt an extension of previous results on the universality of linear halo
bias in neutrino cosmologies, to nonlinear and nonlocal corrections finding
consistent results with the power spectrum analysis.

We present measurements of the Baryon Acoustic Oscillation (BAO) scale in
redshiftspace 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 sphericallyaveraged 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 bestfit cosmology. All
of our results are consistent with basic largescale 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 distanceladder. 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 (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.

We describe the Sloan Digital Sky Survey IV (SDSSIV), a project encompassing
three major spectroscopic programs. The Apache Point Observatory Galactic
Evolution Experiment 2 (APOGEE2) is observing hundreds of thousands of Milky
Way stars at high resolution and high signaltonoise ratio in the
nearinfrared. The Mapping Nearby Galaxies at Apache Point Observatory (MaNGA)
survey is obtaining spatiallyresolved 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 Xray
AGN and galaxies in Xray clusters, and the Time Domain Spectroscopic Survey
(TDSS), obtaining spectra of variable sources. All programs use the 2.5meter
Sloan Foundation Telescope at Apache Point Observatory; observations there
began in Summer 2014. APOGEE2 also operates a second nearinfrared
spectrograph at the 2.5meter 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,
SDSSIV provides regularly scheduled public data releases; the first one, Data
Release 13, was made available in July 2016.

Galaxy clustering data from current and upcoming large scale structure
surveys can provide strong constraints on primordial nonGaussianity through
the scaledependent halo bias. To fully exploit the information from galaxy
surveys, optimal analysis methods need to be developed and applied to the data.
Since the halo bias is sensitive to local nonGaussianity predominately at
large scales, the volume of a given survey is crucial. Consequently, for such
analyses we do not want to split into redshift bins, which would lead to
information loss due to edge effects, but instead analyse the full sample. We
present an optimal technique to directly constrain local nonGaussianity
parametrised by $f_\mathrm{NL}^\mathrm{loc}$, from galaxy clustering by
applying redshift weights to the galaxies. We derive a set of weights to
optimally measure the amplitude of local nonGaussianity,
$f_\mathrm{NL}^\mathrm{loc}$, discuss the redshift weighted power spectrum
estimators, outline the implementation procedure and test our weighting scheme
against Lognormal catalogs for two different surveys: the quasar sample of the
Extended Baryon Oscillation Spectroscopic Survey (eBOSS) and the emission line
galaxy sample of the Dark Energy Spectroscopic Instrument (DESI) survey. We
find an improvement of 30 percent for eBOSS and 6 percent for DESI compared to
the standard Feldman, Kaiser $\&$ Peacock weights, although these predictions
are sensitive to the bias model assumed.

The low statistical errors on cosmological parameters promised by future
galaxy surveys will only be realised with the development of new, fast,
analysis methods that reduce potential systematic problems to low levels. We
present an efficient method for measuring the evolution of the growth of
structure using Redshift Space Distortions (RSD), that removes the need to make
measurements in redshift shells. We provide sets of galaxyweights that cover a
wide range in redshift, but are optimised to provide differential information
about cosmological evolution. These are derived to optimally measure the
coefficients of a parameterisation of the redshiftdependent matter density,
which provides a framework to measure deviations from the concordance
$\Lambda$CDM cosmology, allowing for deviations in both geometric and/or
growth. We test the robustness of the weights by comparing with alternative
schemes and investigate the impact of galaxy bias. We extend the results to
measure the combined anisotropic Baryon Acoustic Oscillation (BAO) and RSD
signals.

Observed galaxy clustering exhibits local transverse statistical isotropy
around the lineofsight (LOS). The variation of the LOS across a galaxy survey
complicates the measurement of the observed clustering as a function of the
angle to the LOS, as fast Fourier transforms (FFTs) based on Cartesian grids,
cannot individually allow for this. Recent advances in methodology for
calculating LOSdependent clustering in Fourier space include the realization
that power spectrum LOSdependent moments can be constructed from sums over
galaxies, based on approximating the LOS to each pair of galaxies by the LOS to
one of them. We show that we can implement this method using multiple FFTs,
each measuring the LOSweighted clustering along different axes. The N log(N)
nature of FFTs means that the computational speedup is a factor of >1000
compared with summing over galaxies. This development should be beneficial for
future projects such as DESI and Euclid which will provide an order of
magnitude more galaxies than current surveys.