Precision measurements of cosmic microwave background (CMB) polarization require extreme control of instrumental systematics. In a companion paper we have presented cosmological constraints from observations with the BICEP2 and Keck Array experiments up to and including the 2015 observing season (BK15), resulting in the deepest CMB polarization maps to date and a statistical sensitivity to the tensor-to-scalar ratio of $\sigma(r) = 0.020$. In this work we characterize the beams and constrain potential systematic contamination from main beam shape mismatch at the three BK15 frequencies (95, 150, and 220 GHz). Far-field maps of 7,360 distinct beam patterns taken from 2010-2015 are used to measure differential beam parameters and predict the contribution of temperature-to-polarization leakage to the BK15 B-mode maps. In the multifrequency, multicomponent likelihood analysis that uses BK15, Planck, and WMAP maps to separate sky components, we find that adding this predicted leakage to simulations induces a bias of $\Delta r = 0.0027 \pm 0.0019$. Future results using higher-quality beam maps and improved techniques to detect such leakage in CMB data will substantially reduce this uncertainty, enabling the levels of systematics control needed for BICEP Array and other experiments that plan to definitively probe large-field inflation.

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 self-calibration 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$K-arcmin) 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 Chern-Simons electromagnetic term $g_{a\gamma}\leq 7.2\times 10^{-2}/H_I$ (95% confidence) than the constraint derived from the B-mode 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.

The BICEP and Keck Array experiments are a suite of small-aperture refracting telescopes observing the microwave sky from the South Pole. They target the degree-scale B-mode polarization signal imprinted in the Cosmic Microwave Background (CMB) by primordial gravitational waves. Such a measurement would shed light on the physics of the very early universe. While BICEP2 observed for the first time a B-mode signal at 150 GHz, higher frequencies from the Planck satellite showed that it could be entirely due to the polarized emission from Galactic dust, though uncertainty remained high. Keck Array has been observing the same region of the sky for several years, with an increased detector count, producing the deepest polarized CMB maps to date. New detectors at 95 GHz were installed in 2014, and at 220 GHz in 2015. These observations enable a better constraint of galactic foreground emissions, as presented here. In 2015, BICEP2 was replaced by BICEP3, a 10 times higher throughput telescope observing at 95 GHz, while Keck Array is now focusing on higher frequencies. In the near future, BICEP Array will replace Keck Array, and will allow unprecedented sensitivity to the gravitational wave signal. High resolution observations from the South Pole Telescope (SPT) will also be used to remove the lensing contribution to B-modes.

We present results from an analysis of all data taken by the BICEP2/Keck CMB polarization experiments up to and including the 2015 observing season. This includes the first Keck Array observations at 220 GHz and additional observations at 95 & 150 GHz. The $Q/U$ maps reach depths of 5.2, 2.9 and 26 $\mu$K$_{cmb}$ arcmin at 95, 150 and 220 GHz respectively over an effective area of $\approx 400$ square degrees. The 220 GHz maps achieve a signal-to-noise on polarized dust emission approximately equal to that of Planck at 353 GHz. We take auto- and cross-spectra between these maps and publicly available WMAP and Planck maps at frequencies from 23 to 353 GHz. We evaluate the joint likelihood of the spectra versus a multicomponent model of lensed-$\Lambda$CDM+$r$+dust+synchrotron+noise. The foreground model has seven parameters, and we impose priors on some of these using external information from Planck and WMAP derived from larger regions of sky. The model is shown to be an adequate description of the data at the current noise levels. The likelihood analysis yields the constraint $r_{0.05}<0.07$ at 95% confidence, which tightens to $r_{0.05}<0.06$ in conjunction with Planck temperature measurements and other data. The lensing signal is detected at $8.8 \sigma$ significance. Running maximum likelihood search on simulations we obtain unbiased results and find that $\sigma(r)=0.020$. These are the strongest constraints to date on primordial gravitational waves.

We present the design and characterization of a ground-based absolute polarization angle calibrator accurate to better than 0.1 degrees for use with polarization sensitive cosmic microwave background (CMB) experiments. The calibrator's accuracy requirement is driven by the need to reduce upper limits on cosmic polarization rotation, which is expected to be zero in a large class of cosmological models. Cosmic polarization effects such as cosmic birefringence and primordial magnetic fields can generate spurious B-modes that result in non-zero CMB TB and EB correlations that are degenerate with a misalignment of detector orientation. Common polarized astrophysical sources used for absolute polarization angle calibration have not been characterized to better than 0.5 degrees. Higher accuracy can be achieved through self-calibration methods, however these are subject to astrophysical foreground contamination and inherently assume the absence of effects like cosmic polarization rotation. The deficiencies in these two calibration methods highlight the need for a well characterized polarized source. The calibrator we present utilizes a 76 GHz Gunn oscillator coupled to a frequency doubler, pyramidal horn antenna, and co-rotating wire-grid polarizer. We use an accurate bubble level in combination with four precision-grade aluminum planes located within the enclosure to calibrate the source's linear polarization plane with respect to the local gravity vector to better than the 0.1 degree goal. In 2017 the calibrator was deployed for an engineering test run on the POLARBEAR CMB experiment located in Chile's Atacama Desert and is being upgraded for calibration of the POLARBEAR-2b receiver in 2018. In the following work we present a detailed overview of the calibrator design, systematic control, characterization, deployment, and plans for future CMB experiment absolute polarization calibration.

BICEP Array is a degree-scale Cosmic Microwave Background (CMB) experiment that will search for primordial B-mode polarization while constraining Galactic foregrounds. BICEP Array will be comprised of four receivers to cover a broad frequency range with channels at 30/40, 95, 150 and 220/270 GHz. The first low-frequency receiver will map synchrotron emission at 30 and 40 GHz and will deploy to the South Pole at the end of 2019. In this paper, we give an overview of the BICEP Array science and instrument, with a focus on the detector module. We designed corrugations in the metal frame of the module to suppress unwanted interactions with the antenna-coupled detectors that would otherwise deform the beams of edge pixels. This design reduces the residual beam systematics and temperature-to-polarization leakage due to beam steering and shape mismatch between polarized beam pairs. We report on the simulated performance of single- and wide-band corrugations designed to minimize these effects. Our optimized design alleviates beam differential ellipticity caused by the metal frame to about 7% over 57% bandwidth (25 to 45 GHz), which is close to the level due the bare antenna itself without a metal frame. Initial laboratory measurements are also presented.

Targeting faint polarization patterns arising from Primordial Gravitational Waves in the Cosmic Microwave Background requires excellent observational sensitivity. Optical elements in small aperture experiments such as Bicep3 and Keck Array are designed to optimize throughput and minimize losses from transmission, reflection and scattering at millimeter wavelengths. As aperture size increases, cryostat vacuum windows must withstand larger forces from atmospheric pressure and the solution has often led to a thicker window at the expense of larger transmission loss. We have identified a new candidate material for the fabrication of vacuum windows: with a tensile strength two orders of magnitude larger than previously used materials, woven high-modulus polyethylene could allow for dramatically thinner windows, and therefore significantly reduced losses and higher sensitivity. In these proceedings we investigate the suitability of high-modulus polyethylene windows for ground-based CMB experiments, such as current and future receivers in the Bicep/Keck Array program. This includes characterizing their optical transmission as well as their mechanical behavior under atmospheric pressure. We find that such ultra-thin materials are promising candidates to improve the performance of large-aperture instruments at millimeter wavelengths, and outline a plan for further tests ahead of a possible upcoming field deployment of such a science-grade window.

Bicep Array is a cosmic microwave background (CMB) polarization experiment that will begin observing at the South Pole in early 2019. This experiment replaces the five Bicep2 style receivers that compose the Keck Array with four larger Bicep3 style receivers observing at six frequencies from 30 to 270GHz. The 95GHz and 150GHz receivers will continue to push the already deep Bicep/Keck CMB maps while the 30/40GHz and 220/270GHz receivers will constrain the synchrotron and galactic dust foregrounds respectively. Here we report on the design and performance of the Bicep Array instruments focusing on the mount and cryostat systems.

BICEP Array is the newest multi-frequency instrument in the BICEP/Keck Array program. It is comprised of four 550 mm aperture refractive telescopes observing the polarization of the cosmic microwave background (CMB) at 30/40, 95, 150 and 220/270 GHz with over 30,000 detectors. We present an overview of the receiver, detailing the optics, thermal, mechanical, and magnetic shielding design. BICEP Array follows BICEP3's modular focal plane concept, and upgrades to 6" wafer to reduce fabrication with higher detector count per module. The first receiver at 30/40 GHz is expected to start observing at the South Pole during the 2019-20 season. By the end of the planned BICEP Array program, we project $\sigma(r) \sim 0.003$, assuming current modeling of polarized Galactic foreground and depending on the level of delensing that can be achieved with higher resolution maps from the South Pole Telescope.

BICEP3 is a 520mm aperture on-axis refracting telescope observing the polarization of the cosmic microwave background (CMB) at 95GHz in search of the B-mode signal originating from inflationary gravitational waves. BICEP3's focal plane is populated with modularized tiles of antenna-coupled transition edge sensor (TES) bolometers. BICEP3 was deployed to the South Pole during 2014-15 austral summer and has been operational since. During the 2016-17 austral summer, we implemented changes to optical elements that lead to better noise performance. We discuss this upgrade and show the performance of BICEP3 at its full mapping speed from the 2017 and 2018 observing seasons. BICEP3 achieves an order-of-magnitude improvement in mapping speed compared to a Keck 95GHz receiver. We demonstrate $6.6\mu K\sqrt{s}$ noise performance of the BICEP3 receiver.

Warm Comptonization models for the soft X-ray excess in AGN do not self-consistently explain the relationship between the Comptonizing medium and the underlying accretion disc. Because of this, they cannot directly connect the fitted Comptonization temperatures and optical depths to accretion disc parameters. Since bulk velocities exceed thermal velocities in highly radiation pressure dominated discs, in these systems bulk Comptonization by turbulence may provide a physical basis in the disc itself for warm Comptonization models. We model the dependence of bulk Comptonization on fundamental accretion disc parameters, such as mass, luminosity, radius, spin, inner boundary condition, and $\alpha$. In addition to constraining warm Comptonization models, our model can help distinguish contributions from bulk Comptonization to the soft X-ray excess from those due to other physical mechanisms, such as absorption and reflection. By linking the time variability of bulk Comptonization to fluctuations in the disc vertical structure due to magnetorotational instability (MRI) turbulence, our results show that observations of the soft X-ray excess can be used to study disc turbulence in the radiation pressure dominated regime. Because our model connects bulk Comptonization to one dimensional vertical structure temperature profiles in a physically intuitive way, it will be useful for understanding this effect in future simulations run in new regimes.

Radiation pressure dominated accretion discs may have turbulent velocities that exceed the electron thermal velocities. Bulk Comptonization by the turbulence may therefore dominate over thermal Comptonization in determining the emergent spectrum. We discuss how to self-consistently resolve and interpret this effect in calculations of spectra of radiation MHD simulations. In particular, we show that this effect is dominated by radiation viscous dissipation and can be treated as thermal Comptonization with an equivalent temperature. We investigate whether bulk Comptonization may provide a physical basis for warm Comptonization models of the soft X-ray excess in AGN. We characterize our results with temperatures and optical depths to make contact with other models of this component. We show that bulk Comptonization shifts the Wien tail to higher energy and lowers the gas temperature, broadening the spectrum. More generally, we model the dependence of this effect on a wide range of fundamental accretion disc parameters, such as mass, luminosity, radius, spin, inner boundary condition, and alpha. Because our model connects bulk Comptonization to one dimensional vertical structure temperature profiles in a physically intuitive way, it will be useful for understanding this effect in future simulations run in new regimes. We also develop a global Monte Carlo code to study this effect in global radiation MHD simulations. This code can be used more broadly to compare global simulations with observed systems, and in particular to investigate whether magnetically dominated discs can explain why observed high Eddington accretion discs appear to be thermally stable.

Bulk velocities exceed thermal velocities for sufficiently radiation pressure dominated accretion flows. We model the contribution of bulk Comptonization to the soft X-ray excess in AGN. Bulk Comptonization is due to both turbulence and the background shear. We calculate spectra both taking into account and not taking into account bulk velocities using scaled data from radiation magnetohydrodynamic (MHD) shearing box simulations. We characterize our results with temperatures and optical depths to make contact with other warm Comptonization models of the soft excess. We chose our fiducial mass, $M = 2 \times 10^6 M_{\odot}$, and accretion rate, $L/L_{\rm Edd} = 2.5$, to correspond to those fit to the super-Eddington narrow line Seyfert 1 (NLS1) RE1034+396. The temperatures, optical depths, and Compton $y$ parameters we find broadly agree with those fit to RE1034+396. The effect of bulk Comptonization is to shift the Wien tail to higher energy and lower the gas temperature, broadening the spectrum. Observations of the soft excess in NLS1s can constrain the properties of disc turbulence if the bulk Comptonization contribution can be separated out from contributions from other physical effects, such as reflection and absorption.

In the last two decades we have seen important mutual stimulations between the community working on electrodynamics of continuous media and the community working on spacetime structure. This is highlighted by the publication of two important monographs from two communities: Foundations of Classical Electrodynamics by F. W. Hehl and Yu. N. Obukhov (Birkh\"auser, Boston 2003) and Differential Forms in Electromagnetics by I. V. Lindell (IEEE Press-Wiley, Piscataway, NJ 2004; see also a new book "Multiforms, Dyadics, and Electromagnetic Media" by the same author in 2015). Starting around 1960, magnetoelectric effects and magnetoelectric media have been a focus of study. Somewhat later, the constitutive tensor density framework was used to construct spacetime structure theoretically and empirically. Earlier, in putting Maxwell equations into a form compatible with general relativity, Einstein noticed that the Maxwell equations can be formulated in a form independent of the metric gravitational potential in 1916; only the constitutive tensor density depends on the metric. This was followed by further stimulating and clarifying works of Weyl in 1918, Murnaghan in 1921, Kottler in 1922 and Cartan in 1923. Year 2016 is the centennial of the start of this historical development. International Journal of Modern Physics D of World Scientific is celebrating the occasion by publishing this special issue on Spacetime Structure and Electrodynamics. We have invited researchers to submit original research articles that would contribute to a better understanding of the phenomena inherent to spacetime structure and electrodynamic media.

One of the most powerful probes of new physics is the polarized Cosmic Microwave Background (CMB). The detection of a nonzero polarization angle rotation between the CMB surface of last scattering and today could provide evidence of Lorentz-violating physics. The purpose of this paper is twofold. First we review one popular mechanism for polarization rotation of CMB photons: the pseudo-Nambu-Goldstone boson. Second, we propose a method to use the Polarbear experiment to constrain Lorentz-violating physics in the context of the Standard-Model Extension, a framework to standardize a large class of potential Lorentz-violating terms in particle physics.

POLARBEAR-2 (PB-2) is a cosmic microwave background (CMB) polarization experiment that will be located in the Atacama highland in Chile at an altitude of 5200 m. Its science goals are to measure the CMB polarization signals originating from both primordial gravitational waves and weak lensing. PB-2 is designed to measure the tensor to scalar ratio, r, with precision {\sigma}(r) < 0.01, and the sum of neutrino masses, {\Sigma}m{\nu}, with {\sigma}({\Sigma}m{\nu}) < 90 meV. To achieve these goals, PB-2 will employ 7588 transition-edge sensor bolometers at 95 GHz and 150 GHz, which will be operated at the base temperature of 250 mK. Science observations will begin in 2017.

A linear polarization field on the sphere can be uniquely decomposed into an E-mode and a B-mode component. These two components are analytically defined in terms of spin-2 spherical harmonics. Maps that contain filtered modes on a partial sky can also be decomposed into E-mode and B-mode 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 B-mode 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 E-modes to a level corresponding to a tensor-to-scalar 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$K-arcmin) of these maps makes it possible to directly reconstruct the lensing potential using only information at larger angular scales ($\ell\leq 700$). From the auto-spectrum 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 no-lensing hypothesis at 5.8$\sigma$, which is the highest significance achieved to date using an EB lensing estimator. Taking the cross-spectrum 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 B-mode auto-spectrum ($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 B-modes 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 B-mode spectrum starts to place constraints on any alternative cosmological sources of B-modes at these angular scales.

Radiation pressure dominated accretion discs around compact objects may have turbulent velocities that greatly exceed the electron thermal velocities within the disc. Bulk Comptonization by the turbulence may therefore dominate over thermal Comptonization in determining the emergent spectrum. Bulk Comptonization by divergenceless turbulence is due to radiation viscous dissipation only. It can be treated as thermal Comptonization by solving the Kompaneets equation with an equivalent "wave" temperature, which is a weighted sum over the power present at each scale in the turbulent cascade. Bulk Comptonization by turbulence with non-zero divergence is due to both pressure work and radiation viscous dissipation. Pressure work has negligible effect on photon spectra in the limit of optically thin turbulence, and in this limit radiation viscous dissipation alone can be treated as thermal Comptonization with a temperature equivalent to the full turbulent power. In the limit of extremely optically thick turbulence, radiation viscous dissipation is suppressed, and the evolution of local photon spectra can be understood in terms of compression and expansion of the strongly coupled photon and gas fluids. We discuss the consequences of these effects for self-consistently resolving and interpreting turbulent Comptonization in spectral calculations in radiation MHD simulations of high luminosity accretion flows.

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 cross-spectra between these maps and publicly available maps from WMAP and Planck at frequencies from 23 GHz to 353 GHz. An excess over lensed-LCDM 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 multi-component model including lensed-LCDM, dust, synchrotron and a possible contribution from inflationary gravitational waves (as parametrized by the tensor-to-scalar 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 model-dependent) 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.

The polarization of the Cosmic Microwave Background (CMB) provides a plethora of information about the early universe. Most notably, gravitational waves from the Inflationary epoch (the leading explanation of the origin of the universe) create a unique CMB polarization $B$-mode signal. An unambiguous detection of the inflationary $B$-mode signal would be a window into the physics of the universe as it was $10^{-36}$ seconds after the Big Bang, at energy scales many orders of magnitude larger than what the LHC can produce. However, there are several instrumental and astrophysical sources that can obfuscate the inflationary $B$-mode signal. One of the most difficult parameters to calibrate for CMB telescopes is the absolute orientation of the antenna's polarization sensitive axis. A miscalibration of the polarization orientation rotates the much brighter $E$-mode signal, producing a false $B$-mode signal. The current best uncertainty on polarization orientation in the CMB community is $0.5^\circ$, set from extrapolating IRAM measurements of the Crab Nebula supernova remnant at 90 GHz to 150 GHz, where the CMB signals peak. This accuracy is not sufficient to convincingly detect $B$-modes predicted by currently allowable models of Inflation. We suggest to precisely measure the Crab Nebula's polarization, which can be calibrated absolutely to $0.1^\circ$ from measurements of the polarized emission of Mars, and use these data to calibrate current and upcoming CMB experiments. In addition to inflationary $B$-modes, more precise calibration will allow us to better constrain the sum of the neutrino masses and set limits on exotic physics such as parity violation through cosmic polarization rotation.

We constrain anisotropic cosmic birefringence using four-point correlations of even-parity $E$-mode and odd-parity $B$-mode polarization in the cosmic microwave background measurements made by the POLARization of the Background Radiation (POLARBEAR) experiment in its first season of observations. We find that the anisotropic cosmic birefringence signal from any parity-violating processes is consistent with zero. The Faraday rotation from anisotropic cosmic birefringence can be compared with the equivalent quantity generated by primordial magnetic fields if they existed. The POLARBEAR nondetection translates into a 95% confidence level (C.L.) upper limit of 93 nanogauss (nG) on the amplitude of an equivalent primordial magnetic field inclusive of systematic uncertainties. This four-point correlation constraint on Faraday rotation is about 15 times tighter than the upper limit of 1380 nG inferred from constraining the contribution of Faraday rotation to two-point correlations of $B$-modes measured by Planck in 2015. Metric perturbations sourced by primordial magnetic fields would also contribute to the $B$-mode power spectrum. Using the POLARBEAR measurements of the $B$-mode power spectrum (two-point correlation), we set a 95% C.L. upper limit of 3.9 nG on primordial magnetic fields assuming a flat prior on the field amplitude. This limit is comparable to what was found in the Planck 2015 two-point correlation analysis with both temperature and polarization. We perform a set of systematic error tests and find no evidence for contamination. This work marks the first time that anisotropic cosmic birefringence or primordial magnetic fields have been constrained from the ground at subdegree scales.

We present an overview of the design and status of the \Pb-2 and the Simons Array experiments. \Pb-2 is a Cosmic Microwave Background polarimetry experiment which aims to characterize the arc-minute angular scale B-mode signal from weak gravitational lensing and search for the degree angular scale B-mode signal from inflationary gravitational waves. The receiver has a 365~mm diameter focal plane cooled to 270~milli-Kelvin. The focal plane is filled with 7,588 dichroic lenslet-antenna coupled polarization sensitive Transition Edge Sensor (TES) bolometric pixels that are sensitive to 95~GHz and 150~GHz bands simultaneously. The TES bolometers are read-out by SQUIDs with 40 channel frequency domain multiplexing. Refractive optical elements are made with high purity alumina to achieve high optical throughput. The receiver is designed to achieve noise equivalent temperature of 5.8~$\mu$K$_{CMB}\sqrt{s}$ in each frequency band. \Pb-2 will deploy in 2016 in the Atacama desert in Chile. The Simons Array is a project to further increase sensitivity by deploying three \Pb-2 type receivers. The Simons Array will cover 95~GHz, 150~GHz and 220~GHz frequency bands for foreground control. The Simons Array will be able to constrain tensor-to-scalar ratio and sum of neutrino masses to $\sigma(r) = 6\times 10^{-3}$ at $r = 0.1$ and $\sum m_\nu (\sigma =1)$ to 40 meV.

In a companion paper we have reported a $>5\sigma$ detection of degree scale $B $-mode polarization at 150 GHz by the BICEP2 experiment. Here we provide a detailed study of potential instrumental systematic contamination to that measurement. We focus extensively on spurious polarization that can potentially arise from beam imperfections. We present a heuristic classification of beam imperfections according to their symmetries and uniformities, and discuss how resulting contamination adds or cancels in maps that combine observations made at multiple orientations of the telescope about its boresight axis. We introduce a technique, which we call "deprojection", for filtering the leading order beam-induced contamination from time ordered data, and show that it removes power from BICEP2's $BB$ spectrum consistent with predictions using high signal-to-noise beam shape measurements. We detail the simulation pipeline that we use to directly simulate instrumental systematics and the calibration data used as input to that pipeline. Finally, we present the constraints on $BB$ contamination from individual sources of potential systematics. We find that systematics contribute $BB$ power that is a factor $\sim10\times$ below BICEP2's 3-year statistical uncertainty, and negligible compared to the observed $BB$ signal. The contribution to the best-fit tensor/scalar ratio is at a level equivalent to $r=(3-6)\times10^{-3}$.

Atmosphere is one of the most important noise sources for ground-based cosmic microwave background (CMB) experiments. By increasing optical loading on the detectors, it amplifies their effective noise, while its fluctuations introduce spatial and temporal correlations between detected signals. We present a physically motivated 3d-model of the atmosphere total intensity emission in the millimeter and sub-millimeter wavelengths. We derive a new analytical estimate for the correlation between detectors time-ordered data as a function of the instrument and survey design, as well as several atmospheric parameters such as wind, relative humidity, temperature and turbulence characteristics. Using an original numerical computation, we examine the effect of each physical parameter on the correlations in the time series of a given experiment. We then use a parametric-likelihood approach to validate the modeling and estimate atmosphere parameters from the POLARBEAR-I project first season data set. We derive a new 1.0% upper limit on the linear polarization fraction of atmospheric emission. We also compare our results to previous studies and weather station measurements. The proposed model can be used for realistic simulations of future ground-based CMB observations.

We describe a low-cost, open-access, CubeSat-based calibration instrument that is designed to support ground-based and sub-orbital experiments searching for various polarization signals in the cosmic microwave background (CMB). All modern CMB polarization experiments require a robust calibration program that will allow the effects of instrument-induced signals to be mitigated during data analysis. A bright, compact, and linearly polarized astrophysical source with polarization properties known to adequate precision does not exist. Therefore, we designed a space-based millimeter-wave calibration instrument, called CalSat, to serve as an open-access calibrator, and this paper describes the results of our design study. The calibration source on board CalSat is composed of five "tones" with one each at 47.1, 80.0, 140, 249 and 309 GHz. The five tones we chose are well matched to (i) the observation windows in the atmospheric transmittance spectra, (ii) the spectral bands commonly used in polarimeters by the CMB community, and (iii) The Amateur Satellite Service bands in the Table of Frequency Allocations used by the Federal Communications Commission. CalSat would be placed in a polar orbit allowing visibility from observatories in the Northern Hemisphere, such as Mauna Kea in Hawaii and Summit Station in Greenland, and the Southern Hemisphere, such as the Atacama Desert in Chile and the South Pole. CalSat also would be observable by balloon-borne instruments launched from a range of locations around the world. This global visibility makes CalSat the only source that can be observed by all terrestrial and sub-orbital observatories, thereby providing a universal standard that permits comparison between experiments using appreciably different measurement approaches.

The Keck Array is a system of cosmic microwave background (CMB) polarimeters, each similar to the BICEP2 experiment. In this paper we report results from the 2012 and 2013 observing seasons, during which the Keck Array consisted of five receivers all operating in the same (150 GHz) frequency band and observing field as BICEP2. We again find an excess of B-mode power over the lensed-$\Lambda$CDM expectation of $> 5 \sigma$ in the range $30 < \ell < 150$ and confirm that this is not due to systematics using jackknife tests and simulations based on detailed calibration measurements. In map difference and spectral difference tests these new data are shown to be consistent with BICEP2. Finally, we combine the maps from the two experiments to produce final Q and U maps which have a depth of 57 nK deg (3.4 $\mu$K arcmin) over an effective area of 400 deg$^2$ for an equivalent survey weight of 250,000 $\mu$K$^{-2}$. The final BB band powers have noise uncertainty a factor of 2.3 times better than the previous results, and a significance of detection of excess power of $> 6\sigma$.

BICEP2 and the Keck Array are polarization-sensitive microwave telescopes that observe the cosmic microwave background (CMB) from the South Pole at degree angular scales in search of a signature of inflation imprinted as B-mode polarization in the CMB. BICEP2 was deployed in late 2009, observed for three years until the end of 2012 at 150 GHz with 512 antenna-coupled transition edge sensor bolometers, and has reported a detection of B-mode polarization on degree angular scales. The Keck Array was first deployed in late 2010 and will observe through 2016 with five receivers at several frequencies (95, 150, and 220 GHz). BICEP2 and the Keck Array share a common optical design and employ the field-proven BICEP1 strategy of using small-aperture, cold, on-axis refractive optics, providing excellent control of systematics while maintaining a large field of view. This design allows for full characterization of far-field optical performance using microwave sources on the ground. Here we describe the optical design of both instruments and report a full characterization of the optical performance and beams of BICEP2 and the Keck Array at 150 GHz.

BICEP1 is a millimeter-wavelength telescope designed specifically to measure the inflationary B-mode polarization of the Cosmic Microwave Background (CMB) at degree angular scales. We present results from an analysis of the data acquired during three seasons of observations at the South Pole (2006 to 2008). This work extends the two-year result published in Chiang et al. (2010), with additional data from the third season and relaxed detector-selection criteria. This analysis also introduces a more comprehensive estimation of band-power window functions, improved likelihood estimation methods and a new technique for deprojecting monopole temperature-to-polarization leakage which reduces this class of systematic uncertainty to a negligible level. We present maps of temperature, E- and B-mode polarization, and their associated angular power spectra. The improvement in the map noise level and polarization spectra error bars are consistent with the 52% increase in integration time relative to Chiang et al. (2010). We confirm both self-consistency of the polarization data and consistency with the two-year results. We measure the angular power spectra at 21 <= l <= 335 and find that the EE spectrum is consistent with Lambda Cold Dark Matter (LCDM) cosmology, with the first acoustic peak of the EE spectrum now detected at 15sigma. The BB spectrum remains consistent with zero. From B-modes only, we constrain the tensor-to-scalar ratio to r = 0.03+0.27-0.23, or r < 0.70 at 95% confidence level.

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 cross-correlation in $B$-modes at high significance. We fit the single- and cross-frequency 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 tensor-to-scalar 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 map-based 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.

We have developed antenna-coupled transition-edge sensor (TES) bolometers for a wide range of cosmic microwave background (CMB) polarimetry experiments, including BICEP2, Keck Array, and the balloon borne SPIDER. These detectors have reached maturity and this paper reports on their design principles, overall performance, and key challenges associated with design and production. Our detector arrays repeatedly produce spectral bands with 20%-30% bandwidth at 95, 150, or 220~GHz. The integrated antenna arrays synthesize symmetric co-aligned beams with controlled side-lobe levels. Cross-polarized response on boresight is typically ~0.5%, consistent with cross-talk in our multiplexed readout system. End-to-end optical efficiencies in our cameras are routinely 35% or higher, with per detector sensitivities of NET~300 uKrts. Thanks to the scalability of this design, we have deployed 2560 detectors as 1280 matched pairs in Keck Array with a combined instantaneous sensitivity of ~9 uKrts, as measured directly from CMB maps in the 2013 season. Similar arrays have recently flown in the SPIDER instrument, and development of this technology is ongoing.

Recent measurements of the Cosmic Microwave Background $B$-mode polarization power spectrum by the BICEP2 and POLARBEAR experiments have demonstrated new precision tools for probing fundamental physics. Regardless of origin, the fact that we can detect sub-$\mu$K CMB polarization represents a tremendous technological breakthrough. Yet more information may be latent in the CMB's polarization pattern. Because of its tensorial nature, CMB polarization may also reveal parity-violating physics via a detection of cosmic polarization rotation. Although current CMB polarimeters are sensitive enough to measure one degree-level polarization rotation with $>5\sigma$ statistical significance, they lack the ability to differentiate this effect from a systematic instrumental polarization rotation. Here, we motivate the search for cosmic polarization rotation from current CMB data as well as independent radio galaxy and quasar polarization measurements. We argue that an improvement in calibration accuracy would allow the precise measurement of parity- and Lorentz-violating effects. We describe the CalSat space-based polarization calibrator that will provide stringent control of systematic polarization angle calibration uncertainties to $0.05^\circ$ -- an order of magnitude improvement over current CMB polarization calibrators. CalSat-based calibration could be used with current CMB polarimeters searching for $B$-mode polarization, effectively turning them into probes of cosmic parity violation, i.e. without the need to build dedicated instruments.

We report on the design and performance of the BICEP2 instrument and on its three-year data set. BICEP2 was designed to measure the polarization of the cosmic microwave background (CMB) on angular scales of 1 to 5 degrees ($\ell$=40-200), near the expected peak of the B-mode polarization signature of primordial gravitational waves from cosmic inflation. Measuring B-modes requires dramatic improvements in sensitivity combined with exquisite control of systematics. The BICEP2 telescope observed from the South Pole with a 26~cm aperture and cold, on-axis, refractive optics. BICEP2 also adopted a new detector design in which beam-defining slot antenna arrays couple to transition-edge sensor (TES) bolometers, all fabricated on a common substrate. The antenna-coupled TES detectors supported scalable fabrication and multiplexed readout that allowed BICEP2 to achieve a high detector count of 500 bolometers at 150 GHz, giving unprecedented sensitivity to B-modes at degree angular scales. After optimization of detector and readout parameters, BICEP2 achieved an instrument noise-equivalent temperature of 15.8 $\mu$K sqrt(s). The full data set reached Stokes Q and U map depths of 87.2 nK in square-degree pixels (5.2 $\mu$K arcmin) over an effective area of 384 square degrees within a 1000 square degree field. These are the deepest CMB polarization maps at degree angular scales to date. The power spectrum analysis presented in a companion paper has resulted in a significant detection of B-mode polarization at degree scales.

(abridged for arXiv) We report results from the BICEP2 experiment, a cosmic microwave background (CMB) polarimeter specifically designed to search for the signal of inflationary gravitational waves in the B-mode power spectrum around $\ell\sim80$. The telescope comprised a 26 cm aperture all-cold refracting optical system equipped with a focal plane of 512 antenna coupled transition edge sensor 150 GHz bolometers each with temperature sensitivity of $\approx300\mu\mathrm{K}_\mathrm{CMB}\sqrt{s}$. BICEP2 observed from the South Pole for three seasons from 2010 to 2012. A low-foreground region of sky with an effective area of 380 square deg was observed to a depth of 87 nK deg in Stokes $Q$ and $U$. We find an excess of $B$-mode power over the base lensed-LCDM expectation in the range $30< \ell< 150$, inconsistent with the null hypothesis at a significance of $> 5\sigma$. Through jackknife tests and simulations we show that systematic contamination is much smaller than the observed excess. We also examine a number of available models of polarized dust emission and find that at their default parameter values they predict power $\sim(5-10)\times$ smaller than the observed excess signal. However, these models are not sufficiently constrained to exclude the possibility of dust emission bright enough to explain the entire excess signal. Cross correlating BICEP2 against 100 GHz maps from the BICEP1 experiment, the excess signal is confirmed and its spectral index is found to be consistent with that of the CMB, disfavoring dust at $1.7\sigma$. The observed $B$-mode power spectrum is well fit by a lensed-LCDM + tensor theoretical model with tensor-to-scalar ratio $r=0.20^{+0.07}_{-0.05}$, with $r=0$ disfavored at $7.0\sigma$. Accounting for the contribution of foreground dust will shift this value downward by an amount which will be better constrained with upcoming data sets.

Cosmic Microwave Background (CMB) polarimeters aspire to measure the faint $B$-mode signature predicted to arise from inflationary gravitational waves. They also have the potential to constrain cosmic birefringence which would produce non-zero expectation values for the CMB's $TB$ and $EB$ spectra. However, instrumental systematic effects can also cause these $TB$ and $EB$ correlations to be non-zero. In particular, an overall miscalibration of the polarization orientation of the detectors produces $TB$ and $EB$ spectra which are degenerate with isotropic cosmological birefringence, while also introducing a small but predictable bias on the $BB$ spectrum. The \bicep three-year spectra, which use our standard calibration of detector polarization angles from a dielectric sheet, are consistent with a polarization rotation of $\alpha = -2.77^\circ \pm 0.86^\circ \text{(statistical)} \pm 1.3^\circ \text{(systematic)}$. We revise the estimate of systematic error on the polarization rotation angle from the two-year analysis by comparing multiple calibration methods. We investigate the polarization rotation for the \bicep 100 GHz and 150 GHz bands separately to investigate theoretical models that produce frequency-dependent cosmic birefringence. We find no evidence in the data supporting either these models or Faraday rotation of the CMB polarization by the Milky Way galaxy's magnetic field. If we assume that there is no cosmic birefringence, we can use the $TB$ and $EB$ spectra to calibrate detector polarization orientations, thus reducing bias of the cosmological $B$-mode spectrum from leaked $E$-modes due to possible polarization orientation miscalibration. After applying this "self-calibration" process, we find that the upper limit on the tensor-to-scalar ratio decreases slightly, from $r<0.70$ to $r<0.65$ at $95\%$ confidence.

The design and performance of a wide bandwidth linear polarization modulator based on the Faraday effect is described. Faraday Rotation Modulators (FRMs) are solid-state polarization switches that are capable of modulation up to ~10 kHz. Six FRMs were utilized during the 2006 observing season in the Background Imaging of Cosmic Extragalactic Polarization (BICEP) experiment; three FRMs were used at each of BICEP's 100 and 150 GHz frequency bands. The technology was verified through high signal-to-noise detection of Galactic polarization using two of the six FRMs during four observing runs in 2006. The features exhibit strong agreement with BICEP's measurements of the Galaxy using non-FRM pixels and with the Galactic polarization models. This marks the first detection of high signal-to-noise mm-wave celestial polarization using fast, active optical modulation. The performance of the FRMs during periods when they were not modulated was also analyzed and compared to results from BICEP's 43 pixels without FRMs.

The BICEP2 and Keck Array experiments are designed to measure the polarization of the cosmic microwave background (CMB) on angular scales of 2-4 degrees (l=50-100). This is the region in which the B-mode signal, a signature prediction of cosmic inflation, is expected to peak. BICEP2 was deployed to the South Pole at the end of 2009 and is in the middle of its third year of observing with 500 polarization-sensitive detectors at 150 GHz. The Keck Array was deployed to the South Pole at the end of 2010, initially with three receivers--each similar to BICEP2. An additional two receivers have been added during the 2011-12 summer. We give an overview of the two experiments, report on substantial gains in the sensitivity of the two experiments after post-deployment optimization, and show preliminary maps of CMB polarization from BICEP2.

The Keck Array is a cosmic microwave background (CMB) polarimeter that will begin observing from the South Pole in late 2010. The initial deployment will consist of three telescopes similar to BICEP2 housed in ultra-compact, pulse tube cooled cryostats. Two more receivers will be added the following year. In these proceedings we report on the design and performance of the Keck cryostat. We also report some initial results on the performance of antenna-coupled TES detectors operating in the presence of a pulse tube. We find that the performance of the detectors is not seriously impacted by the replacement of BICEP2's liquid helium cryostat with a pulse tube cooled cryostat.

In addition to its potential to probe the Inflationary cosmological paradigm, millimeter-wave polarimetry is a powerful tool for studying the Milky Way galaxy's composition and magnetic field structure. Towards this end, presented here are Stokes I, Q, and U maps of the Galactic plane from the millimeter-wave polarimeter BICEP covering the Galactic longitude range 260 - 340 degrees in three atmospheric transmission windows centered on 100, 150, and 220 GHz. The maps sample an optical depth 1 < AV < 30, and are consistent with previous characterizations of the Galactic millimeter-wave frequency spectrum and the large-scale magnetic field structure permeating the interstellar medium. Polarized emission is detected over the entire region within two degrees of the Galactic plane and indicates that the large-scale magnetic field is oriented parallel to the plane of the Galaxy. An observed trend of decreasing polarization fraction with increasing total intensity rules out the simplest model of a constant Galactic magnetic field throughout the Galaxy. Including WMAP data in the analysis, the degree-scale frequency spectrum of Galactic polarization fraction is plotted between 23 and 220 GHz for the first time. A generally increasing trend of polarization fraction with electromagnetic frequency is found, which varies from 0.5%-1.5%at frequencies below 50 GHz to 2.5%-3.5%above 90 GHz. The BICEP and WMAP data are fit to a two-component (synchrotron and dust) model showing that the higher frequency BICEP data are necessary to tightly constrain the amplitude and spectral index of Galactic dust. Furthermore, the dust amplitude predicted by this two-component fit is consistent with model predictions of dust emission in the BICEP bands.

BICEP is a ground-based millimeter-wave bolometric array designed to target the primordial gravity wave signature on the polarization of the cosmic microwave background (CMB) at degree angular scales. Currently in its third year of operation at the South Pole, BICEP is measuring the CMB polarization with unprecedented sensitivity at 100 and 150 GHz in the cleanest available 2% of the sky, as well as deriving independent constraints on the diffuse polarized foregrounds with select observations on and off the Galactic plane. Instrument calibrations are discussed in the context of rigorous control of systematic errors, and the performance during the first two years of the experiment is reviewed.