
We present the strongest constraints to date on anisotropies of cosmic
microwave background (CMB) polarization rotation derived from 150 GHz data
taken by the BICEP2/Keck Array CMB experiments up to and including the 2014
observing season (BK14). The definition of the polarization angle in BK14 maps
has gone through selfcalibration in which the overall angle is adjusted to
minimize the observed TB and EB power spectra. After this procedure, the QU
maps lose sensitivity to a uniform polarization rotation but are still
sensitive to anisotropies of polarization rotation. This analysis places
constraints on the anisotropies of polarization rotation, which could be
generated by CMB photons interacting with axionlike pseudoscalar fields or
Faraday rotation induced by primordial magnetic fields. The sensitivity of BK14
maps ($\sim 3\mu$Karcmin) makes it possible to reconstruct anisotropies of the
polarization rotation angle and measure their angular power spectrum much more
precisely than previous attempts. Our data are found to be consistent with no
polarization rotation anisotropies, improving the upper bound on the amplitude
of the rotation angle spectrum by roughly an order of magnitude compared to the
previous best constraints. Our results lead to an order of magnitude better
constraint on the coupling constant of the ChernSimons electromagnetic term
$g_{a\gamma}\leq 7.2\times 10^{2}/H_I$ (95% confidence) than the constraint
derived from the Bmode spectrum, where $H_I$ is the inflationary Hubble scale.
This constraint leads to a limit on the decay constant of $10^{6}\lesssim
f_a/M_{\rm pl}$ at mass range of $10^{33}< m_a< 10^{28}$ eV for $r=0.01$,
assuming $g_{a\gamma}\sim\alpha/(2\pi f_a)$ with $\alpha$ denoting the fine
structure constant. The upper bound on the amplitude of the primordial magnetic
fields is 30nG (95% confidence) from the polarization rotation anisotropies.

We present measurements of polarization lensing using the 150 GHz maps which
include all data taken by the BICEP2 & Keck Array CMB polarization experiments
up to and including the 2014 observing season (BK14). Despite their modest
angular resolution ($\sim 0.5^\circ$), the excellent sensitivity ($\sim
3\mu$Karcmin) of these maps makes it possible to directly reconstruct the
lensing potential using only information at larger angular scales ($\ell\leq
700$). From the autospectrum of the reconstructed potential we measure an
amplitude of the spectrum to be $A^{\phi\phi}_{\rm L}=1.15\pm 0.36$ (Planck
$\Lambda$CDM prediction corresponds to $A^{\phi\phi}_{\rm L}=1$), and reject
the nolensing hypothesis at 5.8$\sigma$, which is the highest significance
achieved to date using an EB lensing estimator. Taking the crossspectrum of
the reconstructed potential with the Planck 2015 lensing map yields
$A^{\phi\phi}_{\rm L}=1.13\pm 0.20$. These direct measurements of
$A^{\phi\phi}_{\rm L}$ are consistent with the $\Lambda$CDM cosmology, and with
that derived from the previously reported BK14 Bmode autospectrum ($A^{\rm
BB}_{\rm L}=1.20\pm 0.17$). We perform a series of null tests and consistency
checks to show that these results are robust against systematics and are
insensitive to analysis choices. These results unambiguously demonstrate that
the Bmodes previously reported by BICEP / Keck at intermediate angular scales
($150\lesssim\ell\lesssim 350$) are dominated by gravitational lensing. The
good agreement between the lensing amplitudes obtained from the lensing
reconstruction and Bmode spectrum starts to place constraints on any
alternative cosmological sources of Bmodes at these angular scales.

We present results from an analysis of all data taken by the BICEP2 & Keck
Array CMB polarization experiments up to and including the 2014 observing
season. This includes the first Keck Array observations at 95 GHz. The maps
reach a depth of 50 nK deg in Stokes $Q$ and $U$ in the 150 GHz band and 127 nK
deg in the 95 GHz band. We take auto and crossspectra between these maps and
publicly available maps from WMAP and Planck at frequencies from 23 GHz to 353
GHz. An excess over lensedLCDM is detected at modest significance in the
95x150 $BB$ spectrum, and is consistent with the dust contribution expected
from our previous work. No significant evidence for synchrotron emission is
found in spectra such as 23x95, or for correlation between the dust and
synchrotron sky patterns in spectra such as 23x353. We take the likelihood of
all the spectra for a multicomponent model including lensedLCDM, dust,
synchrotron and a possible contribution from inflationary gravitational waves
(as parametrized by the tensortoscalar ratio $r$), using priors on the
frequency spectral behaviors of dust and synchrotron emission from previous
analyses of WMAP and Planck data in other regions of the sky. This analysis
yields an upper limit $r_{0.05}<0.09$ at 95% confidence, which is robust to
variations explored in analysis and priors. Combining these $B$mode results
with the (more modeldependent) constraints from Planck analysis of CMB
temperature plus BAO and other data, yields a combined limit $r_{0.05}<0.07$ at
95% confidence. These are the strongest constraints to date on inflationary
gravitational waves.