
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 consider the power spectrum of a biased tracer observed in a finite volume
in the presence of a largescale overdensity and tidal fields. Expanding both
the observed power spectrum and the source fields (linear power spectrum,
scalar overdensity and tidal field tensor) in spherical harmonics, we
explicitly confirm that each $(\ell,m)$ source generates just the corresponding
$(\ell,m)$ modes of power spectrum in real space. In redshift space, each
$(\ell,m)$ source additionally couples only to $(\ell+2n,m)$ modes of tracer
power spectra. This generalizes the Kaiser formula for monopole, quadrupole and
hexadecapole of the power spectrum to all $(\ell,m)$ modes generated to the
second order in perturbation theory. This formalism can find applications in
constraining the supersample covariance and in the local power spectrum based
bispectrum estimators. As an example application, we forecast the ability to
measure these modes a survey with BOSSlike galaxy number densities.

The existence of the cosmic neutrino background is a fascinating prediction
of the hot big bang model. These neutrinos were a dominant component of the
energy density in the early universe and, therefore, played an important role
in the evolution of cosmological perturbations. In particular, fluctuations in
the neutrino density produced a distinct shift in the temporal phase of sound
waves in the primordial plasma, which has recently been detected in the cosmic
microwave background (CMB). In this paper, we report on the first measurement
of this neutrinoinduced phase shift in the spectrum of baryon acoustic
oscillations (BAO) of the BOSS DR12 data. Constraining the acoustic scale using
Planck data, and marginalizing over the effects of neutrinos in the CMB, we
find evidence for a nonzero phase shift at greater than 95% confidence. We
also demonstrate the robustness of this result in simulations and forecasts.
Besides being a new measurement of the cosmic neutrino background, our work is
the first application of the BAO signal to early universe physics and a
nontrivial confirmation of the standard cosmological history.

Understanding cosmic acceleration is one of the key science drivers for
astrophysics and highenergy physics in the coming decade (2014 P5 Report).
With the Large Synoptic Survey Telescope (LSST) and the Dark Energy
Spectroscopic Instrument (DESI) and other new facilities beginning operations
soon, we are entering an exciting phase during which we expect an order of
magnitude improvement in constraints on dark energy and the physics of the
accelerating Universe. This is a key moment for a matching Small Projects
portfolio that can (1) greatly enhance the science reach of these flagship
projects, (2) have immediate scientific impact, and (3) lay the groundwork for
the next stages of the Cosmic Frontier Dark Energy program. In this White
Paper, we outline a balanced portfolio that can accomplish these goals through
a combination of observational, experimental, and theory and simulation
efforts.

We present measurements of the crosscorrelation of the triplyionized carbon
(CIV) forest with quasars using Sloan Digital Sky Survey Data Release 14. The
study exploits a large sample of new quasars from the first two years of
observations by the Extended Baryon Oscillation Spectroscopic Survey (eBOSS).
The CIV forest is a weaker tracer of largescale structure than the Ly$\alpha$
forest, but benefits from being accessible at redshifts $z<2$ where the quasar
number density from eBOSS is high. Our data sample consists of 287,651 CIV
forest quasars in the redshift range $1.4<z<3.5$ and 387,315 tracer quasars
with $1.2<z<3.5$. We measure largescale correlations from CIV absorption
occuring in three distinct quasar restframe wavelength bands of the spectra
referred to as the CIV forest, the SiIV forest and the Ly$\alpha$ forest. From
the combined fit to the quasarCIV crosscorrelations for the CIV forest and
the SiIV forest, the CIV redshiftspace distortion parameter is $\beta_{\rm
CIV}=0.27_{\ 0.14}^{\ +0.16}$ and its combination with the CIV linear
transmission bias parameter is $b_{\rm CIV}(1+\beta_{\rm CIV})=0.0183_{\
0.0014}^{\ +0.0013}$ ($1\sigma$ statistical error) at the mean redshift
$z=2.00$. Splitting the sample at $z=2.2$ to constrain the bias evolution with
redshift yields the powerlaw exponent $\gamma=0.60\pm0.63$, indicating a
significantly weaker redshiftevolution than for the Ly$\alpha$ forest linear
transmission bias. We demonstrate that CIV absorption has the potential to be
used as a probe of baryon acoustic oscillations (BAO). While the current data
set is insufficient for a detection of the BAO peak feature, the final quasar
samples for redshifts $1.4<z<2.2$ from eBOSS and the Dark Energy Spectroscopic
Instrument (DESI) are expected to provide measurements of the isotropic BAO
scale to $\sim7\%$ and $\sim3\%$ precision, respectively, at $z\simeq1.6$.

We derive and numerically implement an algorithm for estimating the 3D power
spectrum of the Lyman$\alpha$ (Ly$\alpha$) forest flux fluctuations. The
algorithm exploits the unique geometry of Ly$\alpha$ forest data to
efficiently measure the crossspectrum between lines of sight as a function of
parallel wavenumber, transverse separation and redshift. The key to fast
evaluation is to approximate the global covariance matrix as blockdiagonal,
where only pixels from the same spectrum are correlated. We then compute the
eigenvectors of the derivative of the signal covariance with respect to
crossspectrum parameters, and project the inversecovarianceweighted spectra
onto them. This acts much like a radial Fourier transform over redshift
windows. The resulting crossspectrum inference is then converted into our
final product, an approximation of the likelihood for the 3D power spectrum
expressed as second order Taylor expansion around a fiducial model. We
demonstrate the accuracy and scalability of the algorithm and comment on
possible extensions. Our algorithm will allow efficient analysis of the
upcoming Dark Energy Spectroscopic Instrument dataset.

We investigate the threepoint correlation between the Lyman$\alpha$ forest
and the CMB weak lensing ($\delta_F \delta_F \kappa$) expressed as the
crosscorrelation between the CMB weak lensing field and local variations in
the forest power spectrum. In addition to the standard gravitational bispectrum
term, we note the existence of a nonstandard systematic term coming from
misestimation of the mean flux over the finite length of Lyman$\alpha$
skewers. We numerically calculate the angular crosspower spectrum and discuss
its features. We integrate it into zerolag correlation function and compare
our predictions with recent results by Doux et al.. We find that our
predictions are statistically consistent with the measurement, and including
the systematic term improves the agreement with the measurement. We comment on
the implication of the response of the Lyman$\alpha$ forest power spectrum to
the longwavelength density perturbations.

We crosscorrelate the positions of damped Lyman$\alpha$ systems (DLAs) and
their parent quasar catalog with a convergence map derived from the Planck
cosmic microwave background (CMB) temperature data. We make consistent
detections of the lensing signal of both samples in both Fourier and
configuration space. By interpreting the excess signal present in the DLA
catalog with respect to the parent quasar catalog as caused by the large scale
structure traced by DLAs, we are able to infer the bias of these objects:
$b_{\rm DLA}=2.6\pm0.9$. These results are consistent with previous
measurements made in crosscorrelation with the Lyman$\alpha$ forest, although
the current noise in the lensing data and the low number density of DLAs limits
the constraining power of this measurement. We discuss the robustness of the
analysis with respect to a number different systematic effects and forecast
prospects of carrying out this measurement with data from future experiments.

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 measurement of baryon acoustic oscillations (BAO) in the
crosscorrelation of quasars with the Ly$\alpha$forest fluxtransmission at a
mean redshift $z=2.40$. The measurement uses the complete SDSSIII data sample:
168,889 forests and 234,367 quasars from the SDSS Data Release DR12. In
addition to the statistical improvement on our previous study using DR11, we
have implemented numerous improvements at the analysis level allowing a more
accurate measurement of this crosscorrelation. We also developed the first
simulations of the crosscorrelation allowing us to test different aspects of
our data analysis and to search for potential systematic errors in the
determination of the BAO peak position. We measure the two ratios
$D_{H}(z=2.40)/r_{d} = 9.01 \pm 0.36$ and $D_{M}(z=2.40)/r_{d} = 35.7 \pm 1.7$,
where the errors include marginalization over the nonlinear velocity of
quasars and the metal  quasar crosscorrelation contribution, among other
effects. These results are within $1.8\sigma$ of the prediction of the
flat$\Lambda$CDM model describing the observed CMB anisotropies. We combine
this study with the Ly$\alpha$forest autocorrelation function
[2017A&A...603A..12B], yielding $D_{H}(z=2.40)/r_{d} = 8.94 \pm 0.22$ and
$D_{M}(z=2.40)/r_{d} = 36.6 \pm 1.2$, within $2.3\sigma$ of the same
flat$\Lambda$CDM model.

Constraints on inflationary $B$modes using Cosmic Microwave Background
polarization data commonly rely on either template cleaning or crossspectra
between maps at different frequencies to disentangle galactic foregrounds from
the cosmological signal. Assumptions about how the foregrounds scale with
frequency are therefore crucial to interpreting the data. Recent results from
the Planck satellite collaboration claim significant evidence for a
decorrelation in the polarization signal of the spatial pattern of galactic
dust between 353 GHz and 217 GHz. Such a decorrelation would suppress power in
the cross spectrum between high frequency maps, where the dust is strong, and
lower frequency maps, where the sensitivity to cosmological $B$modes is
strongest. Alternatively, it would leave residuals in lower frequency maps
cleaned with a template derived from the higher frequency maps. If not
accounted for, both situations would result in an underestimate of the dust
contribution and thus an upward bias on measurements of the tensortoscalar
ratio, $r$. In this paper we revisit this measurement and find that the
nodecorrelation hypothesis cannot be excluded with the Planck data. There are
three main reasons for this: i) there is significant noise bias in cross
spectra between Planck data splits that needs to be accounted for; ii) there is
strong evidence for unknown instrumental systematics whose amplitude we
estimate using alternative Planck data splits; iii) there are significant
correlations between measurements in different sky patches that need to be
taken into account when assessing the statistical significance. Between
$\ell=5590$ and over $72\%$ of the sky, the dust $BB$ correlation between 217
GHz and 353 GHz is $1.001^{+.004/.021}_{.004/.000}$ ($68\%~stat./syst.$) and
shows no significant trend with sky fraction.

The Lyman$\alpha$ forest is a highly nonlinear field with a lot of
information available in the data beyond the power spectrum. The flux
probability distribution function (PDF) has been used as a successful probe of
smallscale physics. In this paper we argue that measuring coefficients of the
Legendre polyonomial expansion of the PDF offers several advantages over
measuring the binned values as is commonly done. In particular, $n$th
coefficient can be expressed as a linear combination of the first $n$ moments,
allowing these coefficients to be measured in the presence of noise and
allowing a clear route for marginalisation over mean flux. Moreover, in the
presence of noise, our numerical work shows that a finite number of
coefficients are well measured with a very sharp transition into noise
dominance. This compresses the available information into a small number of
wellmeasured quantities. We find that the amount of recoverable information is
a very nonlinear function of spectral noise that strongly favors fewer quasars
measured at better signal to noise.

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 study the covariance properties of real space correlation function
estimators  primarily galaxyshear correlations, or galaxygalaxy lensing 
using SDSS data for both shear catalogs and lenses (specifically the BOSS LOWZ
sample). Using mock catalogs of lenses and sources, we disentangle the various
contributions to the covariance matrix and compare them with a simple
analytical model. We show that not subtracting the lensing measurement around
random points from the measurement around the lens sample is equivalent to
performing the measurement using the lens density field instead of the lens
overdensity field. While the measurement using the lens density field is
unbiased (in the absence of systematics), its error is significantly larger due
to an additional term in the covariance. Therefore, this subtraction should be
performed regardless of its beneficial effects on systematics. Comparing the
error estimates from data and mocks for estimators that involve the
overdensity, we find that the errors are dominated by the shape noise and lens
clustering, that empirically estimated covariances (jackknife and standard
deviation across mocks) are consistent with theoretical estimates, and that
both the connected parts of the 4point function and the supersample
covariance can be neglected for the current levels of noise. While the
tradeoff between different terms in the covariance depends on the survey
configuration (area, source number density), the diagnostics that we use in
this work should be useful for future works to test their
empiricallydetermined covariances.

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.

The squeezedlimit bispectrum, which is generated by nonlinear gravitational
evolution as well as inflationary physics, measures the correlation of three
wavenumbers, in the configuration where one wavenumber is much smaller than the
other two. Since the squeezedlimit bispectrum encodes the impact of a
largescale fluctuation on the smallscale power spectrum, it can be understood
as how the smallscale power spectrum "responds" to the largescale
fluctuation. Viewed in this way, the squeezedlimit bispectrum can be
calculated using the response approach even in the cases which do not submit to
perturbative treatment. To illustrate this point, we apply this approach to the
crosscorrelation between the largescale quasar density field and smallscale
Lyman$\alpha$ forest flux power spectrum. In particular, using separate
universe simulations which implement changes in the largescale density,
velocity gradient, and primordial power spectrum amplitude, we measure how the
Lyman$\alpha$ forest flux power spectrum responds to the local,
longwavelength quasar overdensity, and equivalently their squeezedlimit
bispectrum. We perform a Fisher forecast for the ability of future experiments
to constrain local nonGaussianity using the bispectrum of quasars and the
Lyman$\alpha$ forest. Combining with quasar and Lyman$\alpha$ forest power
spectra to constrain the biases, we find that for DESI the expected $1\sigma$
constraint is ${\rm err}[f_{\rm NL}]\sim60$. Ability for DESI to measure
$f_{\rm NL}$ through this channel is limited primarily by the aliasing and
instrumental noise of the Lyman$\alpha$ forest flux power spectrum. The
combination of response approach and separate universe simulations provides a
novel technique to explore the constraints from the squeezedlimit bispectrum
between different observables.

We compare a large suite of theoretical cosmological models to observational
data from the cosmic microwave background, baryon acoustic oscillation
measurements of expansion, Type Ia SNe measurements of expansion, redshift
space distortion measurements of the growth of structure, and the local Hubble
constant. Our theoretical models include parametrizations of dark energy as
well as physical models of dark energy and modified gravity. We determine the
constraints on the model parameters, incorporating the redshift space
distortion data directly in the analysis. To determine whether models can be
ruled out, we evaluate the $p$ value (the probability under the model of
obtaining data as bad or worse than the observed data). In our comparison, we
find the well known tension of H$_0$ with the other data; no model resolves
this tension successfully. Among the models we consider, the large scale growth
of structure data does not affect the modified gravity models as a category
particularly differently than dark energy models; it matters for some modified
gravity models but not others, and the same is true for dark energy models. We
compute predicted observables for each model under current observational
constraints, and identify models for which future observational constraints
will be particularly informative.

We use fluxtransmission correlations in \Lya forests to measure the imprint
of baryon acoustic oscillations (BAO). The study uses spectra of 157,783
quasars in the redshift range $2.1\le z \le 3.5$ from the Sloan Digital Sky
Survey (SDSS) Data Release 12 (DR12). Besides the statistical improvements on
our previous studies using SDSS DR9 and DR11, we have implemented numerous
improvements in the analysis procedure, allowing us to construct a physical
model of the correlation function and to investigate potential systematic
errors in the determination of the BAO peak position. The Hubble distance,
$\DHub=c/H(z)$, relative to the sound horizon is $\DHub(z=2.33)/r_d=9.07 \pm
0.31$. The bestdetermined combination of comoving angulardiameter distance,
$\DM$, and the Hubble distance is found to be
$\DHub^{0.7}\DM^{0.3}/r_d=13.94\pm0.35$. This value is $1.028\pm0.026$ times
the prediction of the flat\lcdm model consistent with the cosmic microwave
background (CMB) anisotropy spectrum. The errors include marginalization over
the effects of unidentified highdensity absorption systems and fluctuations in
ultraviolet ionizing radiation. Independently of the CMB measurements, the
combination of our results and other BAO observations determine the open\lcdm
density parameters to be $\om=0.296 \pm 0.029$, $\ol=0.699 \pm 0.100$ and
$\Omega_k = 0.002 \pm 0.119$.

The motion of the solar system with respect to the cosmic rest frame
modulates the monopole of the Epoch of Reionization 21cm signal into a dipole.
This dipole has a characteristic frequency dependence that is dominated by the
frequency derivative of the monopole signal. We argue that although the signal
is weaker by a factor of $\sim100$, there are significant benefits in measuring
the dipole. Most importantly, the direction of the cosmic velocity vector is
known exquisitely well from the cosmic microwave background and is not aligned
with the galaxy velocity vector that modulates the foreground monopole.
Moreover, an experiment designed to measure a dipole can rely on differencing
patches of the sky rather than making an absolute signal measurement, which
helps with some systematic effects.

We measure the gravitational lensing shear signal around dark matter halos
hosting CMASS galaxies using light sources at $z\sim 1$ (background galaxies)
and at the surface of last scattering at $z\sim 1100$ (the cosmic microwave
background). The galaxy shear measurement uses data from the CFHTLenS survey,
and the microwave background shear measurement uses data from the {\it Planck}
satellite. The ratio of shears from these crosscorrelations provides a purely
geometric distance measurement across the longest possible cosmological lever
arm. This is because the matter distribution around the halos, including
uncertainties in galaxy bias and systematic errors such as miscentering,
cancels in the ratio for halos in thin redshift slices. We measure this
distance ratio in three different redshift slices of the CMASS sample, and
combine them to obtain a $17\%$ measurement of the distance ratio,
$r=0.390^{+0.070}_{0.062}$ at an effective redshift of $z=0.53$. This is
consistent with the predicted ratio from the {\it Planck} bestfit $\Lambda$CDM
cosmology of $r=0.419$.

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV groundbased dark
energy experiment that will study baryon acoustic oscillations (BAO) and the
growth of structure through redshiftspace distortions with a widearea 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
magnitudelimited 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.

DESI (Dark Energy Spectropic Instrument) is a Stage IV groundbased dark
energy experiment that will study baryon acoustic oscillations and the growth
of structure through redshiftspace distortions with a widearea galaxy and
quasar redshift survey. The DESI instrument is a roboticallyactuated,
fiberfed spectrograph capable of taking up to 5,000 simultaneous spectra over
a wavelength range from 360 nm to 980 nm. The fibers feed ten threearm
spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500,
depending on wavelength. The DESI instrument will be used to conduct a
fiveyear survey designed to cover 14,000 deg$^2$. This powerful instrument
will be installed at prime focus on the 4m Mayall telescope in Kitt Peak,
Arizona, along with a new optical corrector, which will provide a threedegree
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.

We present the largescale 3point correlation function (3PCF) of the SDSS
DR12 CMASS sample of $777,202$ Luminous Red Galaxies, the largestever sample
used for a 3PCF or bispectrum measurement. We make the first highsignificance
($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 2point correlation
function (2PCF) and power spectrum to the 3PCF and opens an avenue for deriving
additional cosmological distance information from future largescale structure
redshift surveys such as DESI. Our measured distance scale from the 3PCF is
fairly independent from that derived from the prereconstruction 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 baryondark matter relative velocity resulting from recombinationera
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 2point 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 2point
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 largescale structure surveys such as
DESI to selfprotect against the relative velocity as a possible systematic.

We analyse the Baryon Acoustic Oscillation (BAO) signal of the final Baryon
Oscillation Spectroscopic Survey (BOSS) data release (DR12). Our analysis is
performed in Fourierspace, using the power spectrum monopole and quadrupole.
The dataset includes $1\,198\,006$ galaxies over the redshift range $0.2 < z <
0.75$. We divide this dataset into three (overlapping) redshift bins with the
effective redshifts $\zeff = 0.38$, $0.51$ and $0.61$. We demonstrate the
reliability of our analysis pipeline using Nbody simulations as well as $\sim
1000$ MultiDarkPatchy mock catalogues, which mimic the BOSSDR12 target
selection. We apply density field reconstruction to enhance the BAO
signaltonoise ratio. By including the power spectrum quadrupole we can
separate the lineofsight and angular modes, which allows us to constrain the
angular diameter distance $D_A(z)$ and the Hubble parameter $H(z)$ separately.
We obtain two independent $1.6\%$ and $1.5\%$ constraints on $D_A(z)$ and
$2.9\%$ and $2.3\%$ constraints on $H(z)$ for the low ($\zeff=0.38$) and high
($\zeff=0.61$) redshift bin, respectively. We obtain two independent $1\%$ and
$0.9\%$ constraints on the angular averaged distance $D_V(z)$, when ignoring
the AlcockPaczynski effect. The detection significance of the BAO signal is of
the order of $8\sigma$ (postreconstruction) for each of the three redshift
bins. Our results are in good agreement with the Planck prediction within
$\Lambda$CDM. This paper is part of a set that analyses the final galaxy
clustering dataset from BOSS. The measurements and likelihoods presented here
are combined with others in~\citet{Alam2016} to produce the final cosmological
constraints from BOSS.