The BICEP/Keck experiment (BK) is a series of small-aperture refracting telescopes observing degree-scale Cosmic Microwave Background (CMB) polarization from the South Pole in search of a primordial $B$-mode signature. This $B$-mode signal arises from primordial gravitational waves interacting with the CMB, and has amplitude parametrized by the tensor-to-scalar ratio $r$. Since 2016, BICEP3 and the Keck Array have been observing with 4800 total antenna-coupled transition-edge sensor detectors, with frequency bands spanning 95, 150, 220, and 270 GHz. Here we present the optical performance of these receivers from 2016 to 2019, including far-field beams measured in situ with an improved chopped thermal source and instrument spectral response measured with a field-deployable Fourier Transform Spectrometer. As a pair differencing experiment, an important systematic that must be controlled is the differential beam response between the co-located, orthogonally polarized detectors. We generate per-detector far-field beam maps and the corresponding differential beam mismatch that is used to estimate the temperature-to-polarization leakage in our CMB maps and to give feedback on detector and optics fabrication. The differential beam parameters presented here were estimated using improved low-level beam map analysis techniques, including efficient removal of non-Gaussian noise as well as improved spatial masking. These techniques help minimize systematic uncertainty in the beam analysis, with the goal of constraining the bias on $r$ induced by temperature-to-polarization leakage to be subdominant to the statistical uncertainty. This is essential as we progress to higher detector counts in the next generation of CMB experiments.

Branches of cosmic inflationary models, such as slow-roll inflation, predict a background of primordial gravitational waves that imprints a unique odd-parity B-mode pattern in the Cosmic Microwave Background (CMB) at amplitudes that are within experimental reach. The BICEP/Keck (BK) experiment targets this primordial signature, the amplitude of which is parameterized by the tensor-to-scalar ratio r, by observing the polarized microwave sky through the exceptionally clean and stable atmosphere at the South Pole. B-mode measurements require an instrument with exquisite sensitivity, tight control of systematics, and wide frequency coverage to disentangle the primordial signal from the Galactic foregrounds. BICEP Array represents the most recent stage of the BK program, and comprises four BICEP3-class receivers observing at 30/40, 95, 150 and 220/270 GHz. The 30/40 GHz receiver will be deployed at the South Pole during the 2019/2020 austral summer. After 3 full years of observations with 30,000+ detectors, BICEP Array will measure primordial gravitational waves to a precision $\sigma(r)$ between 0.002 and 0.004, depending on foreground complexity and the degree of lensing removal. In this paper we give an overview of the instrument, highlighting the design features in terms of cryogenics, magnetic shielding, detectors and readout architecture as well as reporting on the integration and tests that are ongoing with the first receiver at 30/40 GHz.

The BICEP/Keck (BK) experiment aims to detect the imprint of primordial gravitational waves in the Cosmic Microwave Background polarization, which would be direct evidence of the inflation theory. While the tensor-to-scalar ratio has been constrained to be r_0.05 < 0.06 at 95% c.l., further improvements on this upper limit are hindered by polarized Galactic foreground emissions and removal of gravitational lensing polarization. The 30/40 GHz receiver of the BICEP Array (BA) will deploy at the end of 2019 and will constrain the synchrotron foreground with unprecedented accuracy within the BK sky patch. We will show the design of the 30/40 GHz detectors and test results summarizing its performance. The low optical and atmospheric loading at these frequencies requires our TES detectors to have low saturation power in order to be photon-noise dominated. To realize the low thermal conductivity required from a 250 mK base temperature, we developed new bolometer leg designs. We will present the relevant measured detector parameters: G, Tc, Rn, Psat , and spectral bands, and noise spectra. We achieved a per bolometer NEP including all noise components of 2.07E-17 W/sqrt(Hz), including an anticipated photon noise level 1.54E-17 W/sqrt(Hz).

Thermal Kinetic Inductance Detectors (TKIDs) combine the excellent noise performance of traditional bolometers with a radio frequency (RF) multiplexing architecture that enables the large detector counts needed for the next generation of millimeter-wave instruments. Here we present dark prototype TKID pixels that demonstrate a noise equivalent power NEP = $2 \times 10^{-17} \mathrm{W}/\mathrm{\sqrt{Hz}} $ with a $1/f$ knee at 0.1 Hz, suitable for background-limited noise performance at 150 GHz from a ground-based site. We discuss the optimizations in the device design and fabrication techniques to realize optimal electrical performance and high quality factors at a bath temperature of 250 mK.

We describe an on-sky demonstration of a microwave-multiplexing readout system in one of the receivers of the Keck Array, a polarimetry experiment observing the cosmic microwave background at the South Pole. During the austral summer of 2018-2019, we replaced the time-division multiplexing readout system with microwave-multiplexing components including superconducting microwave resonators coupled to radio-frequency superconducting quantum interference devices at the sub-Kelvin focal plane, coaxial-cable plumbing and amplification between room temperature and the cold stages, and a SLAC Microresonator Radio Frequency system for the warm electronics. In the range 5-6 GHz, a single coaxial cable reads out 528 channels. The readout system is coupled to transition-edge sensors, which are in turn coupled to 150-GHz slot-dipole phased-array antennas. Observations began in April 2019, and we report here on an initial characterization of the system performance.

The Probe of Inflation and Cosmic Origins (PICO) is a proposed probe-scale space mission consisting of an imaging polarimeter operating in frequency bands between 20 and 800 GHz. We describe the science achievable by PICO, which has sensitivity equivalent to more than 3300 Planck missions, the technical implementation, the schedule and cost.

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 Probe of Inflation and Cosmic Origins (PICO) is an imaging polarimeter that will scan the sky for 5 years in 21 frequency bands spread between 21 and 799 GHz. It will produce full-sky surveys of intensity and polarization with a final combined-map noise level of 0.87 $\mu$K arcmin for the required specifications, equivalent to 3300 Planck missions, and with our current best-estimate would have a noise level of 0.61 $\mu$K arcmin (6400 Planck missions). PICO will either determine the energy scale of inflation by detecting the tensor to scalar ratio at a level $r=5\times 10^{-4}~(5\sigma)$, or will rule out with more than $5\sigma$ all inflation models for which the characteristic scale in the potential is the Planck scale. With LSST's data it could rule out all models of slow-roll inflation. PICO will detect the sum of neutrino masses at $>4\sigma$, constrain the effective number of light particle species with $\Delta N_{\rm eff}<0.06~(2\sigma)$, and elucidate processes affecting the evolution of cosmic structures by measuring the optical depth to reionization with errors limited by cosmic variance and by constraining the evolution of the amplitude of linear fluctuations $\sigma_{8}(z)$ with sub-percent accuracy. Cross-correlating PICO's map of the thermal Sunyaev-Zeldovich effect with LSST's gold sample of galaxies will precisely trace the evolution of thermal pressure with $z$. PICO's maps of the Milky Way will be used to determine the make up of galactic dust and the role of magnetic fields in star formation efficiency. With 21 full sky legacy maps in intensity and polarization, which cannot be obtained in any other way, the mission will enrich many areas of astrophysics. PICO is the only single-platform instrument with the combination of sensitivity, angular resolution, frequency bands, and control of systematic effects that can deliver this compelling, timely, and broad science.

We have developed a flexible radio-frequency readout system suitable for a variety of superconducting detectors commonly used in millimeter and submillimeter astrophysics, including Kinetic Inductance detectors (KIDs), Thermal KID bolometers (TKIDs), and Quantum Capacitance Detectors (QCDs). Our system avoids custom FPGA-based readouts and instead uses commercially available software radio hardware for ADC/DAC and a GPU to handle real time signal processing. Because this system is written in common C++/CUDA, the range of different algorithms that can be quickly implemented make it suitable for the readout of many others cryogenic detectors and for the testing of different and possibly more effective data acquisition schemes.

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.

The Probe of Inflation and Cosmic Origins (PICO) is a NASA-funded study of a Probe-class mission concept. The top-level science objectives are to probe the physics of the Big Bang by measuring or constraining the energy scale of inflation, probe fundamental physics by measuring the number of light particles in the Universe and the sum of neutrino masses, to measure the reionization history of the Universe, and to understand the mechanisms driving the cosmic star formation history, and the physics of the galactic magnetic field. PICO would have multiple frequency bands between 21 and 799 GHz, and would survey the entire sky, producing maps of the polarization of the cosmic microwave background radiation, of galactic dust, of synchrotron radiation, and of various populations of point sources. Several instrument configurations, optical systems, cooling architectures, and detector and readout technologies have been and continue to be considered in the development of the mission concept. We will present a snapshot of the baseline mission concept currently under development.

The Probe of Inflation and Cosmic Origins (PICO) is a probe-class mission concept currently under study by NASA. PICO will probe the physics of the Big Bang and the energy scale of inflation, constrain the sum of neutrino masses, measure the growth of structures in the universe, and constrain its reionization history by making full sky maps of the cosmic microwave background with sensitivity 80 times higher than the Planck space mission. With bands at 21-799 GHz and arcmin resolution at the highest frequencies, PICO will make polarization maps of Galactic synchrotron and dust emission to observe the role of magnetic fields in Milky Way's evolution and star formation. We discuss PICO's optical system, focal plane, and give current best case noise estimates. The optical design is a two-reflector optimized open-Dragone design with a cold aperture stop. It gives a diffraction limited field of view (DLFOV) with throughput of 910 square cm sr at 21 GHz. The large 82 square degree DLFOV hosts 12,996 transition edge sensor bolometers distributed in 21 frequency bands and maintained at 0.1 K. We use focal plane technologies that are currently implemented on operating CMB instruments including three-color multi-chroic pixels and multiplexed readouts. To our knowledge, this is the first use of an open-Dragone design for mm-wave astrophysical observations, and the only monolithic CMB instrument to have such a broad frequency coverage. With current best case estimate polarization depth of 0.65 microK(CMB}-arcmin over the entire sky, PICO is the most sensitive CMB instrument designed to date.

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.

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.

We show measurements of thermal kinetic inductance detectors (TKID) intended for millimeter wave cosmology in the 200-300 GHz atmospheric window. The TKID is a type of bolometer which uses the kinetic inductance of a superconducting resonator to measure the temperature of the thermally isolated bolometer island. We measure bolometer thermal conductance, time constant and noise equivalent power. We also measure the quality factor of our resonators as the bath temperature varies to show they are limited by effects consistent with coupling to two level systems.

CMB-S4 is a proposed experiment to map the polarization of the Cosmic Microwave Background (CMB) to nearly the cosmic variance limit for angular scales that are accessible from the ground. The science goals and capabilities of CMB-S4 in illuminating cosmic inflation, measuring the sum of neutrino masses, searching for relativistic relics in the early universe, characterizing dark energy and dark matter, and mapping the matter distribution in the universe have been described in the CMB-S4 Science Book. This Technology Book is a companion volume to the Science Book. The ambitious science goals of CMB-S4, a "Stage-4" experiment, require a step forward in experimental capability from the current Stage=II experiments. To guide this process, we summarize the current state of CMB instrumentation technology, and identify R&D efforts necessary to advance it for use in CMB-S4. The book focuses on technical challenges in four broad areas: Telescope Design; Receiver Optics; Focal-Plane Optical Coupling; and Focal-Plane Sensor and Readout.

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 small-aperture refracting cosmic microwave background (CMB) telescope designed to make sensitive polarization maps in pursuit of a potential B-mode 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 2015-2016 summer deployment season in situ with a thermal chopped source. We generate high-fidelity per-detector beam maps, show the array-averaged beam profile, and characterize the differential beam response between co-located, 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 temperature-to-polarization leakage may be marginalized over or subtracted in analysis - will inform the design of next-generation CMB experiments with many thousands of detectors.

BICEP3 is a 520 mm aperture, compact two-lens refractor designed to observe the polarization of the cosmic microwave background (CMB) at 95 GHz. Its focal plane consists of modularized tiles of antenna-coupled transition edge sensors (TESs), similar to those used in BICEP2 and the Keck Array. The increased per-receiver 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 per-receiver throughput while maintaining the advantages of a compact design, BICEP3-class receivers are ideally suited as building blocks for a 3rd-generation 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 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.

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.

BICEP3 is a $550~mm$ aperture telescope with cold, on-axis, 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 degree-angular scales. BICEP3 is designed to house 1280 dual-polarization pixels, which, when fully-populated, 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 2014-2015 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 multi-frequency observation data from Planck, BICEP2, and the Keck Array, BICEP3 is projected to set upper limits on the tensor-to-scalar ratio to $r$ $\lesssim 0.03$ at $95\%$ C.L..

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

SuperSpec is a novel on-chip spectrometer we are developing for multi-object, moderate resolution (R = 100 - 500), large bandwidth (~1.65:1) submillimeter and millimeter survey spectroscopy of high-redshift galaxies. The spectrometer employs a filter bank architecture, and consists of a series of half-wave resonators formed by lithographically-patterned superconducting transmission lines. The signal power admitted by each resonator is detected by a lumped element titanium nitride (TiN) kinetic inductance detector (KID) operating at 100 - 200 MHz. We have tested a new prototype device that achieves the targeted R = 100 resolving power, and has better detector sensitivity and optical efficiency than previous devices. We employ a new method for measuring photon noise using both coherent and thermal sources of radiation to cleanly separate the contributions of shot and wave noise. We report an upper limit to the detector NEP of $1.4\times10^{-17}$ W Hz$^{-1/2}$, within 10% of the photon noise limited NEP for a ground-based R=100 spectrometer.

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.

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.

SuperSpec is a novel on-chip spectrometer we are developing for multi-object, moderate resolution (R = 100 - 500), large bandwidth (~1.65:1) submillimeter and millimeter survey spectroscopy of high-redshift galaxies. The spectrometer employs a filter bank architecture, and consists of a series of half-wave resonators formed by lithographically-patterned superconducting transmission lines. The signal power admitted by each resonator is detected by a lumped element titanium nitride (TiN) kinetic inductance detector (KID) operating at 100-200 MHz. We have tested a new prototype device that is more sensitive than previous devices, and easier to fabricate. We present a characterization of a representative R=282 channel at f = 236 GHz, including measurements of the spectrometer detection efficiency, the detector responsivity over a large range of optical loading, and the full system optical efficiency. We outline future improvements to the current system that we expect will enable construction of a photon-noise-limited R=100 filter bank, appropriate for a line intensity mapping experiment targeting the [CII] 158 micron transition during the Epoch of Reionization

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.

BICEP3 is a 550 mm-aperture refracting telescope for polarimetry of radiation in the cosmic microwave background at 95 GHz. It adopts the methodology of BICEP1, BICEP2 and the Keck Array experiments - it possesses sufficient resolution to search for signatures of the inflation-induced cosmic gravitational-wave background while utilizing a compact design for ease of construction and to facilitate the characterization and mitigation of systematics. However, BICEP3 represents a significant breakthrough in per-receiver sensitivity, with a focal plane area 5$\times$ larger than a BICEP2/Keck Array receiver and faster optics ($f/1.6$ vs. $f/2.4$). Large-aperture infrared-reflective metal-mesh filters and infrared-absorptive cold alumina filters and lenses were developed and implemented for its optics. The camera consists of 1280 dual-polarization pixels; each is a pair of orthogonal antenna arrays coupled to transition-edge sensor bolometers and read out by multiplexed SQUIDs. Upon deployment at the South Pole during the 2014-15 season, BICEP3 will have survey speed comparable to Keck Array 150 GHz (2013), and will significantly enhance spectral separation of primordial B-mode power from that of possible galactic dust contamination in the BICEP2 observation patch.

An attitude determination system for balloon-borne experiments is presented. The system provides pointing information in azimuth and elevation for instruments flying on stratospheric balloons over Antarctica. In-flight attitude is given by the real-time combination of readings from star cameras, a magnetometer, sun sensors, GPS, gyroscopes, tilt sensors and an elevation encoder. Post-flight attitude reconstruction is determined from star camera solutions, interpolated by the gyroscopes using an extended Kalman Filter. The multi-sensor system was employed by the Balloon-borne Large Aperture Submillimeter Telescope for Polarimetry (BLASTPol), an experiment that measures polarized thermal emission from interstellar dust clouds. A similar system was designed for the upcoming flight of SPIDER, a Cosmic Microwave Background polarization experiment. The pointing requirements for these experiments are discussed, as well as the challenges in designing attitude reconstruction systems for high altitude balloon flights. In the 2010 and 2012 BLASTPol flights from McMurdo Station, Antarctica, the system demonstrated an accuracy of <5' rms in-flight, and <5" rms post-flight.

We present the results of integration and characterization of the SPIDER instrument after the 2013 pre-flight campaign. SPIDER is a balloon-borne polarimeter designed to probe the primordial gravitational wave signal in the degree-scale $B$-mode polarization of the cosmic microwave background. With six independent telescopes housing over 2000 detectors in the 94 GHz and 150 GHz frequency bands, SPIDER will map 7.5% of the sky with a depth of 11 to 14 $\mu$K$\cdot$arcmin at each frequency, which is a factor of $\sim$5 improvement over Planck. We discuss the integration of the pointing, cryogenic, electronics, and power sub-systems, as well as pre-flight characterization of the detectors and optical systems. SPIDER is well prepared for a December 2014 flight from Antarctica, and is expected to be limited by astrophysical foreground emission, and not instrumental sensitivity, over the survey region.

We introduce the light-weight carbon fiber and aluminum gondola designed for the SPIDER balloon-borne telescope. SPIDER is designed to measure the polarization of the Cosmic Microwave Background radiation with unprecedented sensitivity and control of systematics in search of the imprint of inflation: a period of exponential expansion in the early Universe. The requirements of this balloon-borne instrument put tight constrains on the mass budget of the payload. The SPIDER gondola is designed to house the experiment and guarantee its operational and structural integrity during its balloon-borne flight, while using less than 10% of the total mass of the payload. We present a construction method for the gondola based on carbon fiber reinforced polymer tubes with aluminum inserts and aluminum multi-tube joints. We describe the validation of the model through Finite Element Analysis and mechanical tests.

We present the technology and control methods developed for the pointing system of the SPIDER experiment. SPIDER is a balloon-borne polarimeter designed to detect the imprint of primordial gravitational waves in the polarization of the Cosmic Microwave Background radiation. We describe the two main components of the telescope's azimuth drive: the reaction wheel and the motorized pivot. A 13 kHz PI control loop runs on a digital signal processor, with feedback from fibre optic rate gyroscopes. This system can control azimuthal speed with < 0.02 deg/s RMS error. To control elevation, SPIDER uses stepper-motor-driven linear actuators to rotate the cryostat, which houses the optical instruments, relative to the outer frame. With the velocity in each axis controlled in this way, higher-level control loops on the onboard flight computers can implement the pointing and scanning observation modes required for the experiment. We have accomplished the non-trivial task of scanning a 5000 lb payload sinusoidally in azimuth at a peak acceleration of 0.8 deg/s$^2$, and a peak speed of 6 deg/s. We can do so while reliably achieving sub-arcminute pointing control accuracy.

We present the second generation BLASTbus electronics. The primary purposes of this system are detector readout, attitude control, and cryogenic housekeeping, for balloon-borne telescopes. Readout of neutron transmutation doped germanium (NTD-Ge) bolometers requires low noise and parallel acquisition of hundreds of analog signals. Controlling a telescope's attitude requires the capability to interface to a wide variety of sensors and motors, and to use them together in a fast, closed loop. To achieve these different goals, the BLASTbus system employs a flexible motherboard-daughterboard architecture. The programmable motherboard features a digital signal processor (DSP) and field-programmable gate array (FPGA), as well as slots for three daughterboards. The daughterboards provide the interface to the outside world, with versions for analog to digital conversion, and optoisolated digital input/output. With the versatility afforded by this design, the BLASTbus also finds uses in cryogenic, thermometry, and power systems. For accurate timing control to tie everything together, the system operates in a fully synchronous manner. BLASTbus electronics have been successfully deployed to the South Pole, and flown on stratospheric balloons.

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

SuperSpec is an ultra-compact spectrometer-on-a-chip for millimeter and submillimeter wavelength astronomy. Its very small size, wide spectral bandwidth, and highly multiplexed readout will enable construction of powerful multibeam spectrometers for high-redshift observations. The spectrometer consists of a horn-coupled microstrip feedline, a bank of narrow-band superconducting resonator filters that provide spectral selectivity, and Kinetic Inductance Detectors (KIDs) that detect the power admitted by each filter resonator. The design is realized using thin-film lithographic structures on a silicon wafer. The mm-wave microstrip feedline and spectral filters of the first prototype are designed to operate in the band from 195-310 GHz and are fabricated from niobium with at Tc of 9.2K. The KIDs are designed to operate at hundreds of MHz and are fabricated from titanium nitride with a Tc of 2K. Radiation incident on the horn travels along the mm-wave microstrip, passes through the frequency-selective filter, and is finally absorbed by the corresponding KID where it causes a measurable shift in the resonant frequency. In this proceedings, we present the design of the KIDs employed in SuperSpec and the results of initial laboratory testing of a prototype device. We will also briefly describe the ongoing development of a demonstration instrument that will consist of two 500-channel, R=700 spectrometers, one operating in the 1-mm atmospheric window and the other covering the 650 and 850 micron bands.

SuperSpec is a pathfinder for future lithographic spectrometer cameras, which promise to energize extra-galactic astrophysics at (sub)millimeter wavelengths: delivering 200--500 km/s spectral velocity resolution over an octave bandwidth for every pixel in a telescope's field of view. We present circuit simulations that prove the concept, which enables complete millimeter-band spectrometer devices in just a few square-millimeter footprint. We evaluate both single-stage and two-stage channelizing filter designs, which separate channels into an array of broad-band detectors, such as bolometers or kinetic inductance detector (KID) devices. We discuss to what degree losses (by radiation or by absorption in the dielectric) and fabrication tolerances affect the resolution or performance of such devices, and what steps we can take to mitigate the degradation. Such design studies help us formulate critical requirements on the materials and fabrication process, and help understand what practical limits currently exist to the capabilities these devices can deliver today or over the next few years.

We are developing multi-chroic antenna-coupled TES detectors for CMB polarimetry. Multi-chroic detectors increase the mapping speed per focal plane area and provide greater discrimination of polarized galactic foregrounds with no increase in weight or cryogenic cost. In each pixel, a silicon lens-coupled dual polarized sinuous antenna collects light over a two-octave frequency band. The antenna couples the broadband millimeter wave signal into microstrip transmission lines, and on-chip filter banks split the broadband signal into several frequency bands. Separate TES bolometers detect the power in each frequency band and linear polarization. We will describe the design and performance of these devices and present optical data taken with prototype pixels. Our measurements show beams with percent level ellipticity, percent level cross-polarization leakage, and partitioned bands using banks of 2, 3, and 7 filters. We will also describe the development of broadband anti-reflection coatings for the high dielectric constant lens. The broadband anti-reflection coating has approximately 100 percent bandwidth and no detectable loss at cryogenic temperature. Finally, we will describe an upgrade for the POLARBEAR CMB experiment and installation for the LITEBird CMB satellite experiment both of which have focal planes with kilo-pixel of these detectors to achieve unprecedented mapping speed.

Between the BICEP2 and Keck Array experiments, we have deployed over 1500 dual polarized antenna coupled bolometers to map the Cosmic Microwave Background's polarization. We have been able to rapidly deploy these detectors because they are completely planar with an integrated phased-array antenna. Through our experience in these experiments, we have learned of several challenges with this technology- specifically the beam synthesis in the antenna- and in this paper we report on how we have modified our designs to mitigate these challenges. In particular, we discus differential steering errors between the polarization pairs' beam centroids due to microstrip cross talk and gradients of penetration depth in the niobium thin films of our millimeter wave circuits. We also discuss how we have suppressed side lobe response with a Gaussian taper of our antenna illumination pattern. These improvements will be used in Spider, Polar-1, and this season's retrofit of Keck Array.

The Keck Array (SPUD) is a set of microwave polarimeters that observes from the South Pole at degree angular scales in search of a signature of Inflation imprinted as B-mode polarization in the Cosmic Microwave Background (CMB). The first three Keck Array receivers were deployed during the 2010-2011 Austral summer, followed by two new receivers in the 2011-2012 summer season, completing the full five-receiver array. All five receivers are currently observing at 150 GHz. The Keck Array employs the field-proven BICEP/BICEP2 strategy of using small, 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. We describe our efforts to characterize the main beam shape and beam shape mismatch between co-located orthogonally-polarized detector pairs, and discuss the implications of measured differential beam parameters on temperature to polarization leakage in CMB analysis.

The Keck Array (SPUD) began observing the cosmic microwave background's polarization in the winter of 2011 at the South Pole. The Keck Array follows the success of the predecessor experiments Bicep and Bicep2, using five on-axis refracting telescopes. These have a combined imaging array of 2500 antenna-coupled TES bolometers read with a SQUID-based time domain multiplexing system. We will discuss the detector noise and the optimization of the readout. The achieved sensitivity of the Keck Array is 11.5 {\mu}K_(CMB)*sqrt{s} in the 2012 configuration.

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.

Observations of the temperature anisotropy of the Cosmic Microwave Background (CMB) lend support to an inflationary origin of the universe, yet no direct evidence verifying inflation exists. Many current experiments are focussing on the CMB's polarization anisotropy, specifically its curl component (called "B-mode" polarization), which remains undetected. The inflationary paradigm predicts the existence of a primordial gravitational wave background that imprints a unique B-mode signature on the CMB's polarization at large angular scales. The CMB B-mode signal also encodes gravitational lensing information at smaller angular scales, bearing the imprint of cosmological large scale structures (LSS) which in turn may elucidate the properties of cosmological neutrinos. The quest for detection of these signals; each of which is orders of magnitude smaller than the CMB temperature anisotropy signal, has motivated the development of background-limited detectors with precise control of systematic effects. The POLARBEAR experiment is designed to perform a deep search for the signature of gravitational waves from inflation and to characterize lensing of the CMB by LSS. POLARBEAR is a 3.5 meter ground-based telescope with 3.8 arcminute angular resolution at 150 GHz. At the heart of the POLARBEAR receiver is an array featuring 1274 antenna-coupled superconducting transition edge sensor (TES) bolometers cooled to 0.25 Kelvin. POLARBEAR is designed to reach a tensor-to-scalar ratio of 0.025 after two years of observation -- more than an order of magnitude improvement over the current best results, which would test physics at energies near the GUT scale. POLARBEAR had an engineering run in the Inyo Mountains of Eastern California in 2010 and will begin observations in the Atacama Desert in Chile in 2011.

We describe the Cosmic Microwave Background (CMB) polarization experiment called Polarbear. This experiment will use the dedicated Huan Tran Telescope equipped with a powerful 1,200-bolometer array receiver to map the CMB polarization with unprecedented accuracy. We summarize the experiment, its goals, and current status.