
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.

BICEP3, the latest telescope in the BICEP/Keck program, started science
observations in March 2016. It is a 550mm aperture refractive telescope
observing the polarization of the cosmic microwave background at 95 GHz. We
show the focal plane design and detector performance, including spectral
response, optical efficiency and preliminary sensitivity of the upgraded
BICEP3. We demonstrate 9.72$\mu$K$\sqrt{\textrm{s}}$ noise performance of the
BICEP3 receiver.

BICEP3 is a smallaperture refracting cosmic microwave background (CMB)
telescope designed to make sensitive polarization maps in pursuit of a
potential Bmode signal from inflationary gravitational waves. It is the latest
in the BICEP/Keck Array series of CMB experiments at the South Pole, which has
provided the most stringent constraints on inflation to date. For the 2016
observing season, BICEP3 was outfitted with a full suite of 2400 optically
coupled detectors operating at 95 GHz. In these proceedings we report on the
far field beam performance using calibration data taken during the 20152016
summer deployment season in situ with a thermal chopped source. We generate
highfidelity perdetector beam maps, show the arrayaveraged beam profile, and
characterize the differential beam response between colocated, orthogonally
polarized detectors which contributes to the leading instrumental systematic in
pair differencing experiments. We find that the levels of differential
pointing, beamwidth, and ellipticity are similar to or lower than those
measured for BICEP2 and Keck Array. The magnitude and distribution of BICEP3's
differential beam mismatch  and the level to which temperaturetopolarization
leakage may be marginalized over or subtracted in analysis  will inform the
design of nextgeneration CMB experiments with many thousands of detectors.

BICEP3 is a 520 mm aperture, compact twolens refractor designed to observe
the polarization of the cosmic microwave background (CMB) at 95 GHz. Its focal
plane consists of modularized tiles of antennacoupled transition edge sensors
(TESs), similar to those used in BICEP2 and the Keck Array. The increased
perreceiver optical throughput compared to BICEP2/Keck Array, due to both its
faster f/1.7 optics and the larger aperture, more than doubles the combined
mapping speed of the BICEP/Keck program. The BICEP3 receiver was recently
upgraded to a full complement of 20 tiles of detectors (2560 TESs) and is now
beginning its second year of observation (and first science season) at the
South Pole. We report on its current performance and observing plans. Given its
high perreceiver throughput while maintaining the advantages of a compact
design, BICEP3class receivers are ideally suited as building blocks for a
3rdgeneration CMB experiment, consisting of multiple receivers spanning 35 GHz
to 270 GHz with total detector count in the tens of thousands. We present plans
for such an array, the new "BICEP Array" that will replace the Keck Array at
the South Pole, including design optimization, frequency coverage, and
deployment/observing strategies.

A linear polarization field on the sphere can be uniquely decomposed into an
Emode and a Bmode component. These two components are analytically defined in
terms of spin2 spherical harmonics. Maps that contain filtered modes on a
partial sky can also be decomposed into Emode and Bmode components. However,
the lack of full sky information prevents orthogonally separating these
components using spherical harmonics. In this paper, we present a technique for
decomposing an incomplete map into E and Bmode components using E and B
eigenmodes of the pixel covariance in the observed map. This method is found to
orthogonally define E and B in the presence of both partial sky coverage and
spatial filtering. This method has been applied to the BICEP2 and the Keck
Array maps and results in reducing E to B leakage from LCDM Emodes to a level
corresponding to a tensortoscalar ratio of $r<1\times10^{4}$.

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.

BICEP3 is a $550~mm$ aperture telescope with cold, onaxis, refractive optics
designed to observe at the $95~GHz$ band from the South Pole. It is the newest
member of the BICEP/Keck family of inflationary probes specifically designed to
measure the polarization of the cosmic microwave background (CMB) at
degreeangular scales. BICEP3 is designed to house 1280 dualpolarization
pixels, which, when fullypopulated, totals to $\sim$9$\times$ the number of
pixels in a single Keck $95~GHz$ receiver, thus further advancing the
BICEP/Keck program's $95~GHz$ mapping speed. BICEP3 was deployed during the
austral summer of 20142015 with 9 detector tiles, to be increased to its full
capacity of 20 in the second season. After instrument characterization
measurements were taken, CMB observation commenced in April 2015. Together with
multifrequency observation data from Planck, BICEP2, and the Keck Array,
BICEP3 is projected to set upper limits on the tensortoscalar ratio to $r$
$\lesssim 0.03$ at $95\%$ C.L..

We report the results of a joint analysis of data from BICEP2/Keck Array and
Planck. BICEP2 and Keck Array have observed the same approximately 400 deg$^2$
patch of sky centered on RA 0h, Dec. $57.5\deg$. The combined maps reach a
depth of 57 nK deg in Stokes $Q$ and $U$ in a band centered at 150 GHz. Planck
has observed the full sky in polarization at seven frequencies from 30 to 353
GHz, but much less deeply in any given region (1.2 $\mu$K deg in $Q$ and $U$ at
143 GHz). We detect 150$\times$353 crosscorrelation in $B$modes at high
significance. We fit the single and crossfrequency power spectra at
frequencies $\geq 150$ GHz to a lensed$\Lambda$CDM model that includes dust
and a possible contribution from inflationary gravitational waves (as
parameterized by the tensortoscalar ratio $r$), using a prior on the
frequency spectral behavior of polarized dust emission from previous \planck\
analysis of other regions of the sky. We find strong evidence for dust and no
statistically significant evidence for tensor modes. We probe various model
variations and extensions, including adding a synchrotron component in
combination with lower frequency data, and find that these make little
difference to the $r$ constraint. Finally we present an alternative analysis
which is similar to a mapbased cleaning of the dust contribution, and show
that this gives similar constraints. The final result is expressed as a
likelihood curve for $r$, and yields an upper limit $r_{0.05}<0.12$ at 95%
confidence. Marginalizing over dust and $r$, lensing $B$modes are detected at
$7.0\,\sigma$ significance.