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

We present a Bayesian phasespace reconstruction of the cosmic largescale
matter density and velocity fields from the SDSSIII 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 610 $h^{1}$ Mpc, and a lognormalPoisson 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 SunyaevZeldovich (kSZ),
integrated SachsWolfe (ISW) measurements, or environmental studies.

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 nonlinear effects on the BAO feature. Using the
anisotropic clustering of the prereconstruction density field, we measure the
product DM*H from the AlcockPaczynski (AP) effect and the growth of structure,
quantified by f{\sigma}8(z), from redshiftspace 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 Bayesian reconstruction method which maps a galaxy distribution
from redshiftspace to realspace inferring the distances of the individual
galaxies. The method is based on sampling density fields assuming a lognormal
prior with a likelihood given by the negative binomial distribution function
modelling stochastic bias. We assume a deterministic bias given by a power law
relating the dark matter density field to the expected halo or galaxy field.
Coherent redshiftspace distortions are corrected in a Gibbssampling procedure
by moving the galaxies from redshiftspace to realspace according to the
peculiar motions derived from the recovered density field using linear theory
with the option to include tidal field corrections from second order Lagrangian
perturbation theory. The virialised distortions are corrected by sampling
candidate realspace positions (being in the neighbourhood of the observations
along the line of sight), which are compatible with the bulk flow corrected
redshiftspace position adding a random dispersion term in high density
collapsed regions. The latter are defined according to the eigenvalues of the
Hessian. This approach presents an alternative method to estimate the distances
to galaxies using the three dimensional spatial information, and assuming
isotropy. Hence the number of applications is very broad. In this work we show
the potential of this method to constrain the growth rate up to $k$ ~ 0.3 $h$
Mpc$^{1}$. Furthermore it could be useful to correct for photometric redshift
errors, and to obtain improved BAO reconstructions.

We present a Bayesian reconstruction algorithm to generate unbiased samples
of the underlying dark matter field from halo catalogues. Our new contribution
consists of implementing a nonPoisson likelihood including a deterministic
nonlinear and scaledependent bias. In particular we present the Hamiltonian
equations of motions for the negative binomial (NB) probability distribution
function. This permits us to efficiently sample the posterior distribution
function of density fields given a sample of galaxies using the Hamiltonian
Monte Carlo technique implemented in the Argo code. We have tested our
algorithm with the Bolshoi $N$body simulation at redshift $z = 0$, inferring
the underlying dark matter density field from subsamples of the halo catalogue
with biases smaller and larger than one. Our method shows that we can draw
closely unbiased samples (compatible within 1$\sigma$) from the posterior
distribution up to scales of about $k$~1 h/Mpc in terms of powerspectra and
celltocell correlations. We find that a Poisson likelihood yields
reconstructions with power spectra deviating more than 10% at $k$=0.2 h/Mpc.
Our reconstruction algorithm is especially suited for emission line galaxy data
for which a complex nonlinear stochastic biasing treatment beyond Poissonity
becomes indispensable.

We study the statistical inference of the cosmological dark matter density
field from nonGaussian, nonlinear and nonPoisson biased distributed tracers.
We have implemented a Bayesian posterior sampling computercode solving this
problem and tested it with mock data based on Nbody simulations.