Search

We use a custom-made calibrator to measure individual detectors' polarization angles of BICEP3, a small aperture telescope observing the cosmic microwave background (CMB) at 95GHz from the South Pole. We describe our calibration strategy and the statistical and systematic uncertainties associated with the measurement. We reach an unprecedented precision for such measurement on a CMB experiment, with a repeatability for each detector pair of $0.02\deg$. We show that the relative angles measured using this method are in excellent agreement with those extracted from CMB data. Because the absolute measurement is currently limited by a systematic uncertainty, we do not derive cosmic birefringence constraints from BICEP3 data in this work. Rather, we forecast the sensitivity of BICEP3 sky maps for such analysis. We investigate the relative contributions of instrument noise, lensing, and dust, as well as astrophysical and instrumental systematics. We also explore the constraining power of different angle estimators, depending on analysis choices. We establish that the BICEP3 2-year dataset (2017--2018) has an on-sky sensitivity to the cosmic birefringence angle of $\sigma = 0.078\deg$, which could be improved to $\sigma = 0.055\deg$ by adding all of the existing BICEP3 data (through 2023). Furthermore, we emphasize the possibility of using the BICEP3 sky patch as a polarization calibration source for CMB experiments, which with the present data could reach a precision of $0.035\deg$. Finally, in the context of inflation searches, we investigate the impact of detector-to-detector variations in polarization angles as they may bias the tensor-to-scalar ratio r. We show that while the effect is expected to remain subdominant to other sources of systematic uncertainty, it can be reliably calibrated using polarization angle measurements such as the ones we present in this paper., Comment: 29 Pages, 17 Figures, 6 Tables, as submitted to PRD

The BICEP/$\textit{Keck}$ (BK) series of cosmic microwave background (CMB) polarization experiments has, over the past decade and a half, produced a series of field-leading constraints on cosmic inflation via measurements of the "B-mode" polarization of the CMB. Primordial B modes are directly tied to the amplitude of primordial gravitational waves (PGW), their strength parameterized by the tensor-to-scalar ratio, $r$, and thus the energy scale of inflation. Having set the most sensitive constraints to-date on $r$, $\sigma(r)=0.009$ ($r_{0.05}<0.036, 95\%$ C.L.) using data through the 2018 observing season ("BK18"), the BICEP/$\textit{Keck}$ program has continued to improve its dataset in the years since. We give a brief overview of the BK program and the "BK18" result before discussing the program's ongoing efforts, including the deployment and performance of the $\textit{Keck Array}$'s successor instrument, BICEP Array, improvements to data processing and internal consistency testing, new techniques such as delensing, and how those will ultimately serve to allow BK reach $\sigma(r) \lesssim 0.003$ using data through the 2027 observing season., Comment: 9 pages, 5 figures. Contribution to the 2024 Cosmology session of the 58th Rencontres de Moriond

Time-Division Multiplexing is the readout architecture of choice for many ground and space experiments, as it is a very mature technology with proven outstanding low-frequency noise stability, which represents a central challenge in multiplexing. Once fully populated, each of the two BICEP Array high frequency receivers, observing at 150GHz and 220/270GHz, will have 7776 TES detectors tiled on the focal plane. The constraints set by these two receivers required a redesign of the warm readout electronics. The new version of the standard Multi Channel Electronics, developed and built at the University of British Columbia, is presented here for the first time. BICEP Array operates Time Division Multiplexing readout technology to the limits of its capabilities in terms of multiplexing rate, noise and crosstalk, and applies them in rigorously demanding scientific application requiring extreme noise performance and systematic error control. Future experiments like CMB-S4 plan to use TES bolometers with Time Division/SQUID-based readout for an even larger number of detectors., Comment: 10 pages, 7 figures, Submitted to Journal of Low Temperature Physics

We present estimates of line-of-sight distortion fields derived from the 95 GHz and 150 GHz data taken by BICEP2, BICEP3, and Keck Array up to the 2018 observing season, leading to cosmological constraints and a study of instrumental and astrophysical systematics. Cosmological constraints are derived from three of the distortion fields concerning gravitational lensing from large-scale structure, polarization rotation from magnetic fields or an axion-like field, and the screening effect of patchy reionization. We measure an amplitude of the lensing power spectrum $A_L^{\phi\phi}=0.95 \pm 0.20$. We constrain polarization rotation, expressed as the coupling constant of a Chern-Simons electromagnetic term $g_{a\gamma} \leq 2.6 \times 10^{-2}/H_I$, where $H_I$ is the inflationary Hubble parameter, and an amplitude of primordial magnetic fields smoothed over 1 Mpc $B_{1\text{Mpc}} \leq 6.6 \;\text{nG}$ at 95 GHz. We constrain the root mean square of optical-depth fluctuations in a simple "crinkly surface" model of patchy reionization, finding $A^\tau<0.19$ ($2\sigma$) for the coherence scale of $L_c=100$. We show that all of the distortion fields of the 95 GHz and 150 GHz polarization maps are consistent with simulations including lensed-$\Lambda$CDM, dust, and noise, with no evidence for instrumental systematics. In some cases, the EB and TB quadratic estimators presented here are more sensitive than our previous map-based null tests at identifying and rejecting spurious B-modes that might arise from instrumental effects. Finally, we verify that the standard deprojection filtering in the BICEP/Keck data processing is effective at removing temperature to polarization leakage., Comment: 34 pages, 19 figures, accepted for publication in The Astrophysical Journal

We characterize Galactic dust filaments by correlating BICEP/Keck and Planck data with polarization templates based on neutral hydrogen (H I) observations. Dust polarization is important for both our understanding of astrophysical processes in the interstellar medium (ISM) and the search for primordial gravitational waves in the cosmic microwave background (CMB). In the diffuse ISM, H I is strongly correlated with the dust and partly organized into filaments that are aligned with the local magnetic field. We analyze the deep BICEP/Keck data at 95, 150, and 220 GHz, over the low-column-density region of sky where BICEP/Keck has set the best limits on primordial gravitational waves. We separate the H I emission into distinct velocity components and detect dust polarization correlated with the local Galactic H I but not with the H I associated with Magellanic Stream I. We present a robust, multifrequency detection of polarized dust emission correlated with the filamentary H I morphology template down to 95 GHz. For assessing its utility for foreground cleaning, we report that the H I morphology template correlates in B modes at a $\sim$10-65$\%$ level over the multipole range $20 < \ell < 200$ with the BICEP/Keck maps, which contain contributions from dust, CMB, and noise components. We measure the spectral index of the filamentary dust component spectral energy distribution to be $\beta = 1.54 \pm 0.13$. We find no evidence for decorrelation in this region between the filaments and the rest of the dust field or from the inclusion of dust associated with the intermediate velocity H I. Finally, we explore the morphological parameter space in the H I-based filamentary model., Comment: 27 pages, 12 figures

Observations of the Cosmic Microwave Background rely on cryogenic instrumentation with cold detectors, readout, and optics providing the low noise performance and instrumental stability required to make more sensitive measurements. It is therefore critical to optimize all aspects of the cryogenic design to achieve the necessary performance, with low temperature components and acceptable system cooling requirements. In particular, we will focus on our use of thermal filters and cold optics, which reduce the thermal load passed along to the cryogenic stages. To test their performance, we have made a series of in situ measurements while integrating the third receiver for the BICEP Array telescope. In addition to characterizing the behavior of this receiver, these measurements continue to refine the models that are being used to inform design choices being made for future instruments., Comment: 9 pages, 8 figures, Proceedings of SPIE 2022

Constraining the Galactic foregrounds with multi-frequency Cosmic Microwave Background (CMB) observations is an essential step towards ultimately reaching the sensitivity to measure primordial gravitational waves (PGWs), the sign of inflation after the Big-Bang that would be imprinted on the CMB. The BICEP Array telescope is a set of multi-frequency cameras designed to constrain the energy scale of inflation through CMB B-mode searches while also controlling the polarized galactic foregrounds. The lowest frequency BICEP Array receiver (BA1) has been observing from the South Pole since 2020 and provides 30 GHz and 40 GHz data to characterize the Galactic synchrotron in our CMB maps. In this paper, we present the design of the BA1 detectors and the full optical characterization of the camera including the on-sky performance at the South Pole. The paper also introduces the design challenges during the first observing season including the effect of out-of-band photons on detectors performance. It also describes the tests done to diagnose that effect and the new upgrade to minimize these photons, as well as installing more dichroic detectors during the 2022 deployment season to improve the BA1 sensitivity. We finally report background noise measurements of the detectors with the goal of having photon noise dominated detectors in both optical channels. BA1 achieves an improvement in mapping speed compared to the previous deployment season., Comment: Proceedings of SPIE Astronomical Telescopes + Instrumentation 2022 (AS22)

The BICEP3 Polarimeter is a small aperture, refracting telescope, dedicated to the observation of the Cosmic Microwave Background (CMB) at 95GHz. It is designed to target degree angular scale polarization patterns, in particular the very-much-sought-after primordial B-mode signal, which is a unique signature of cosmic inflation. The polarized signal from the sky is reconstructed by differencing co-localized, orthogonally polarized superconducting Transition Edge Sensor (TES) bolometers. In this work, we present absolute measurements of the polarization response of the detectors for more than $\sim 800$ functioning detector pairs of the BICEP3 experiment, out of a total of $\sim 1000$. We use a specifically designed Rotating Polarized Source (RPS) to measure the polarization response at multiple source and telescope boresight rotation angles, to fully map the response over 360 degrees. We present here polarization properties extracted from on-site calibration data taken in January 2022. A similar calibration campaign was performed in 2018, but we found that our constraint was dominated by systematics on the level of $\sim0.5^\circ$. After a number of improvements to the calibration set-up, we are now able to report a significantly lower level of systematic contamination. In the future, such precise measurements will be used to constrain physics beyond the standard cosmological model, namely cosmic birefringence., Comment: Submitted to: SPIE Astronomical Telescopes + Instrumentation (AS22)

For the past decade, the BICEP/Keck collaboration has been operating a series of telescopes at the Amundsen-Scott South Pole Station measuring degree-scale $B$-mode polarization imprinted in the Cosmic Microwave Background (CMB) by primordial gravitational waves (PGWs). These telescopes are compact refracting polarimeters mapping about 2% of the sky, observing at a broad range of frequencies to account for the polarized foreground from Galactic synchrotron and thermal dust emission. Our latest publication "BK18" utilizes the data collected up to the 2018 observing season, in conjunction with the publicly available WMAP and Planck data, to constrain the tensor-to-scalar ratio $r$. It particularly includes (1) the 3-year BICEP3 data which is the current deepest CMB polarization map at the foreground-minimum 95 GHz; and (2) the Keck 220 GHz map with a higher signal-to-noise ratio on the dust foreground than the Planck 353 GHz map. We fit the auto- and cross-spectra of these maps to a multicomponent likelihood model ($\Lambda$CDM+dust+synchrotron+noise+$r$) and find it to be an adequate description of the data at the current noise level. The likelihood analysis yields $\sigma(r)=0.009$. The inference of $r$ from our baseline model is tightened to $r_{0.05}=0.014^{+0.010}_{-0.011}$ and $r_{0.05}<0.036$ at 95% confidence, meaning that the BICEP/Keck $B$-mode data is the most powerful existing dataset for the constraint of PGWs. The up-coming BICEP Array telescope is projected to reach $\sigma(r) \lesssim 0.003$ using data up to 2027., Comment: 8 pages, 6 figures, contribution to the 2022 Cosmology session of the 56th Rencontres de Moriond

The BICEP/Keck Collaboration is currently leading the quest to the highest sensitivity measurements of the polarized CMB anisotropies on degree scale with a series of cryogenic telescopes, of which BICEP Array is the latest Stage-3 upgrade with a total of $\sim32,000$ detectors. The instrument comprises 4 receivers spanning 30 to 270 GHz, with the low-frequency 30/40 GHz deployed to the South Pole Station in late 2019. The full complement of receivers is forecast to set the most stringent constraints on the tensor to scalar ratio $r$. Building on these advances, the overarching small-aperture telescope concept is already being used as the reference for further Stage-4 experiment design. In this paper I will present the development of the BICEP Array 150 GHz detector module and its fabrication requirements, with highlights on the high-density time division multiplexing (TDM) design of the cryogenic circuit boards. The low-impedance wiring required between the detectors and the first-stage SQUID amplifiers is crucial to maintain a stiff voltage bias on the detectors. A novel multi-layer FR4 Printed Circuit Board (PCB) with superconducting traces, capable of reading out up to 648 detectors, is presented along with its validation tests. I will also describe an ultra-high density TDM detector module we developed for a CMB-S4-like experiment that allows up to 1,920 detectors to be read out. TDM has been chosen as the detector readout technology for the Cosmic Microwave Background Stage-4 (CMB-S4) experiment based on its proven low-noise performance, predictable costs and overall maturity of the architecture. The heritage for TDM is rooted in mm- and submm-wave experiments dating back 20 years and has since evolved to support a multiplexing factor of 64x in Stage-3 experiments., Comment: 9 pages, 5 figure, Proceeding of LTD19 submitted to Journal of Low Temperature Physics

We present results from an analysis of all data taken by the BICEP2, Keck Array and BICEP3 CMB polarization experiments up to and including the 2018 observing season. We add additional Keck Array observations at 220 GHz and BICEP3 observations at 95 GHz to the previous 95/150/220 GHz data set. The $Q/U$ maps now reach depths of 2.8, 2.8 and 8.8 $\mu{\mathrm K}_{cmb}$ arcmin at 95, 150 and 220 GHz respectively over an effective area of $\approx 600$ square degrees at 95 GHz and $\approx 400$ square degrees at 150 & 220 GHz. The 220 GHz maps now achieve a signal-to-noise on polarized dust emission exceeding 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 and 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 no longer requires a prior on the frequency spectral index of the dust emission taken from measurements on other regions of the sky. This model is an adequate description of the data at the current noise levels. The likelihood analysis yields the constraint $r_{0.05}<0.036$ at 95% confidence. Running maximum likelihood search on simulations we obtain unbiased results and find that $\sigma(r)=0.009$. These are the strongest constraints to date on primordial gravitational waves., Comment: 22 pages, 24 figures, as published in PRL, data and figures available for download at http://bicepkeck.org

We report on the design and performance of the BICEP3 instrument and its first three-year data set collected from 2016 to 2018. BICEP3 is a 52cm aperture, refracting telescope designed to observe the polarization of the cosmic microwave background (CMB) on degree angular scales at 95GHz. It started science observation at the South Pole in 2016 with 2400 antenna-coupled transition-edge sensor (TES) bolometers. The receiver first demonstrated new technologies such as large-diameter alumina optics, Zotefoam infrared filters, and flux-activated SQUIDs, allowing $\sim 10\times$ higher optical throughput compared to the Keck design. BICEP3 achieved instrument noise-equivalent temperatures of 9.2, 6.8 and 7.1$\mu\text{K}_{\text{CMB}}\sqrt{\text{s}}$ and reached Stokes $Q$ and $U$ map depths of 5.9, 4.4 and 4.4$\mu$K-arcmin in 2016, 2017 and 2018, respectively. The combined three-year data set achieved a polarization map depth of 2.8$\mu$K-arcmin over an effective area of 585 square degrees, which is the deepest CMB polarization map made to date at 95GHz., Comment: 35 pages, 35 figures, as submitted to ApJ, data and figures available for download at http://bicepkeck.org

We present an improved search for axion-like polarization oscillations in the cosmic microwave background (CMB) with observations from the Keck Array. An all-sky, temporally sinusoidal rotation of CMB polarization, equivalent to a time-variable cosmic birefringence, is an observable manifestation of a local axion field and potentially allows a CMB polarimeter to detect axion-like dark matter directly. We describe improvements to the method presented in previous work, and we demonstrate the updated method with an expanded dataset consisting of the 2012-2015 observing seasons. We set limits on the axion-photon coupling constant for mass $m$ in the range $10^{-23}$-$10^{-18}~\mathrm{eV}$, which corresponds to oscillation periods on the order of hours to years. Our results are consistent with the background model. For periods between $1$ and $30~\mathrm{d}$ ($1.6 \times 10^{-21} \leq m \leq 4.8 \times 10^{-20}~\mathrm{eV}$), the $95\%$-confidence upper limits on rotation amplitude are approximately constant with a median of $0.27^\circ$, which constrains the axion-photon coupling constant to $g_{\phi\gamma} < (4.5 \times 10^{-12}~\mathrm{GeV}^{-1}) m/(10^{-21}~\mathrm{eV}$), if axion-like particles constitute all of the dark matter. More than half of the collected BICEP dataset has yet to be analyzed, and several current and future CMB polarimetry experiments can apply the methods presented here to achieve comparable or superior constraints. In the coming years, oscillation measurements can achieve the sensitivity to rule out unexplored regions of the axion parameter space., Comment: 14 pages, 3 figures

The BICEP/Keck Array experiment is a series of small-aperture refracting telescopes observing degree-scale Cosmic Microwave Background polarization from the South Pole in search of a primordial $B$-mode signature. 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 use high-fidelity, in-situ measurements of the beam response to estimate the temperature-to-polarization (T $\rightarrow$ P) leakage in our latest data including observations from 2016 through 2018. This includes three years of BICEP3 observing at 95 GHz, and multifrequency data from Keck Array. Here we present band-averaged far-field beam maps, differential beam mismatch, and residual beam power (after filtering out the leading difference modes via deprojection) for these receivers. We show preliminary results of "beam map simulations," which use these beam maps to observe a simulated temperature (no $Q/U$) sky to estimate T $\rightarrow$ P leakage in our real data., Comment: 9 pages, 4 figures

BICEP3 is a 520 mm aperture on-axis refracting telescope at the South Pole, which observes the polarization of the cosmic microwave background (CMB) at 95 GHz to search for the B-mode signal from inflationary gravitational waves. In addition to this main target, we have developed a low-elevation observation strategy to extend coverage of the Southern sky at the South Pole, where BICEP3 can quickly achieve degree-scale E-mode measurements over a large area. An interesting E-mode measurement is probing a potential polarization anomaly around the CMB Cold Spot. During the austral summer seasons of 2018-19 and 2019-20, BICEP3 observed the sky with a flat mirror to redirect the beams to various low elevation ranges. The preliminary data analysis shows degree-scale E-modes measured with high signal-to-noise ratio., Comment: 12 pages, 10 figures; Figure 7 shows the correct file

The BICEP3 CMB Polarimeter is a small-aperture refracting telescope located at the South Pole and is specifically designed to search for the possible signature of inflationary gravitational waves in the Cosmic Microwave Background (CMB). The experiment measures polarization on the sky by differencing the signal of co-located, orthogonally polarized antennas coupled to Transition Edge Sensor (TES) detectors. We present precise measurements of the absolute polarization response angles and polarization efficiencies for nearly all of BICEP3s $\sim800$ functioning polarization-sensitive detector pairs from calibration data taken in January 2018. Using a Rotating Polarized Source (RPS), we mapped polarization response for each detector over a full 360 degrees of source rotation and at multiple telescope boresight rotations from which per-pair polarization properties were estimated. In future work, these results will be used to constrain signals predicted by exotic physical models such as Cosmic Birefringence., Comment: Proceedings submitted to SPIE 2020 (AS111). 12 pages, 5 figures, 2 tables

A detection of curl-type ($B$-mode) polarization of the primary CMB would be direct evidence for the inflationary paradigm of the origin of the Universe. The BICEP/Keck Array (BK) program targets the degree angular scales, where the power from primordial $B$-mode polarization is expected to peak, with ever-increasing sensitivity and has published the most stringent constraints on inflation to date. BICEP Array (BA) is the Stage-3 instrument of the BK program and will comprise four BICEP3-class receivers observing at 30/40, 95, 150 and 220/270 GHz with a combined 32,000+ detectors; such wide frequency coverage is necessary for control of the Galactic foregrounds, which also produce degree-scale $B$-mode signal. The 30/40 GHz receiver is designed to constrain the synchrotron foreground and has begun observing at the South Pole in early 2020. By the end of a 3-year observing campaign, the full BICEP Array instrument is projected to reach $\sigma_r$ between 0.002 and 0.004, depending on foreground complexity and degree of removal of $B$-modes due to gravitational lensing (delensing). This paper presents an overview of the design, measured on-sky performance and calibration of the first BA receiver. We also give a preview of the added complexity in the time-domain multiplexed readout of the 7,776-detector 150 GHz receiver., Comment: Proceedings of SPIE 2020 (AS111). This article supersedes arXiv:1808.00568 and arXiv:2002.05228

We present a constraint on the tensor-to-scalar ratio, $r$, derived from measurements of cosmic microwave background (CMB) polarization $B$-modes with "delensing," whereby the uncertainty on $r$ contributed by the sample variance of the gravitational lensing $B$-modes is reduced by cross-correlating against a lensing $B$-mode template. This template is constructed by combining an estimate of the polarized CMB with a tracer of the projected large-scale structure. The large-scale-structure tracer used is a map of the cosmic infrared background derived from Planck satellite data, while the polarized CMB map comes from a combination of South Pole Telescope, BICEP/Keck, and Planck data. We expand the BICEP/Keck likelihood analysis framework to accept a lensing template and apply it to the BICEP/Keck data set collected through 2014 using the same parametric foreground modelling as in the previous analysis. From simulations, we find that the uncertainty on $r$ is reduced by $\sim10\%$, from $\sigma(r)$= 0.024 to 0.022, which can be compared with a $\sim26\%$ reduction obtained when using a perfect lensing template. Applying the technique to the real data, the constraint on $r$ is improved from $r_{0.05} < 0.090$ to $r_{0.05} < 0.082$ (95\% C.L.). This is the first demonstration of improvement in an $r$ constraint through delensing., Comment: 23 pages, 11 figures; match published version

We present a search for axion-like polarization oscillations in the cosmic microwave background (CMB) with observations from the Keck Array. A local axion field induces an all-sky, temporally sinusoidal rotation of CMB polarization. A CMB polarimeter can thus function as a direct-detection experiment for axion-like dark matter. We develop techniques to extract an oscillation signal. Many elements of the method are generic to CMB polarimetry experiments and can be adapted for other datasets. As a first demonstration, we process data from the 2012 observing season to set upper limits on the axion-photon coupling constant in the mass range $10^{-21}$-$10^{-18}~\mathrm{eV}$, which corresponds to oscillation periods on the order of hours to months. We find no statistically significant deviations from the background model. For periods larger than $24~\mathrm{hr}$ (mass $m < 4.8 \times 10^{-20}~\mathrm{eV}$), the median 95%-confidence upper limit is equivalent to a rotation amplitude of $0.68^\circ$, which constrains the axion-photon coupling constant to $g_{\phi\gamma} < \left ( 1.1 \times 10^{-11}~\mathrm{GeV}^{-1} \right ) m/\left (10^{-21}~\mathrm{eV} \right )$, if axion-like particles constitute all of the dark matter. The constraints can be improved substantially with data already collected by the BICEP series of experiments. Current and future CMB polarimetry experiments are expected to achieve sufficient sensitivity to rule out unexplored regions of the axion parameter space., Comment: 25 pages, 6 figures, 2 tables

A microwave SQUID multiplexer ($\mu$MUX) has been optimized for coupling to large arrays of superconducting transition-edge sensor (TES) bolometers. We present the scalable cryogenic multiplexer chip design in a 1820-channel multiplexer configuration for the 4-8 GHz rf band. The key metrics of yield, sensitivity, and crosstalk are determined through measurements of 455 readout channels, which span 4-5 GHz. The median white-noise level is 45 pA/$\sqrt{\textrm{Hz}}$, evaluated at 2 Hz, with a 1/f knee $\leq$ 20 mHz after common-mode subtraction. The white-noise level decreases the sensitivity of a TES bolometer optimized for detection of the cosmic microwave background at 150 GHz by only 3%. The measured crosstalk between any channel pair is $\leq$ 0.3%., Comment: 6 pages, 5 figures

Families of cosmic inflation models predict a primordial gravitational-wave background that imprints B-mode polarization pattern in the Cosmic Microwave Background (CMB). High sensitivity instruments with wide frequency coverage and well-controlled systematic errors are needed to constrain the faint B-mode amplitude. We have developed antenna-coupled Transition Edge Sensor (TES) arrays for high-sensitivity polarized CMB observations over a wide range of millimeter-wave bands. BICEP Array, the latest phase of the BICEP/Keck experiment series, is a multi-receiver experiment designed to search for inflationary B-mode polarization to a precision $\sigma$(r) between 0.002 and 0.004 after 3 full years of observations, depending on foreground complexity and the degree of lensing removal. We describe the electromagnetic design and measured performance of BICEP Array low-frequency 40-GHz detector, their packaging in focal plane modules, and optical characterization including efficiency and beam matching between polarization pairs. We summarize the design and simulated optical performance, including an approach to improve the optical efficiency due to mismatch losses. We report the measured beam maps for a new broad-band corrugation design to minimize beam differential ellipticity between polarization pairs caused by interactions with the module housing frame, which helps minimize polarized beam mismatch that converts CMB temperature to polarization ($T \rightarrow P$) anisotropy in CMB maps., Comment: 8 pages, 7 figures, Accepted by the Journal of Low Temperature Physics (Proceedings of the 18th International Workshop on Low Temperature Detectors)

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., Comment: 9 pages, 5 figures, presented at LTD18 in Milan (July 2019), accepted on JLTP (February 2020)

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)., Comment: Accepted for publication in Journal of Low Temperature Physics

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., Comment: 8 pages, 3 figures. Accepted by the Journal of Low Temperature Physics (Proceedings of the 18th International Workshop on Low Temperature Detectors)

The EXperiment for Cryogenic Large-Aperture Intensity Mapping (EXCLAIM) is a cryogenic balloon-borne instrument that will survey galaxy and star formation history over cosmological time scales. Rather than identifying individual objects, EXCLAIM will be a pathfinder to demonstrate an intensity mapping approach, which measures the cumulative redshifted line emission. EXCLAIM will operate at 420-540 GHz with a spectral resolution R=512 to measure the integrated CO and [CII] in redshift windows spanning 0 < z < 3.5. CO and [CII] line emissions are key tracers of the gas phases in the interstellar medium involved in star-formation processes. EXCLAIM will shed light on questions such as why the star formation rate declines at z < 2, despite continued clustering of the dark matter. The instrument will employ an array of six superconducting integrated grating-analog spectrometers (micro-spec) coupled to microwave kinetic inductance detectors (MKIDs). Here we present an overview of the EXCLAIM instrument design and status., Comment: 10 pages, 4 figures

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., Comment: 9 pages, 11 figures, Accepted by the Journal of Low Temperature Physics (Proceedings of the 18th International Workshop on Low Temperature Detectors)

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., Comment: 18 pages, 9 figures, 4 tables. Minor correction to Figure 2 and caption

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., Comment: 23 pages, 23 figures, as accepted by PRL, data and figures available for download at http://bicepkeck.org/

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., Comment: 12 pages and 14 figures

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., Comment: 15 pages, 13 figures

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., Comment: 10 pages, 7 figures. Published in Proc. SPIE. Presented at SPIE Astronomical Telescopes and Instrumentation Conference 10708: Millimeter, Submillimeter, and Far-Infrared Detectors and Instrumentation for Astronomy XI, June 2018

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., Comment: Published in Proc. SPIE. Presented at SPIE Astronomical Telescopes and Instrumentation Conference 10708: Millimeter, Submillimeter, and Far-Infrared Detectors and Instrumentation for Astronomy XI, June 2018. 14 pages, 7 figures

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., Comment: 12 pages, 9 figures, Published in Proc. SPIE. Presented at SPIE Astronomical Telescopes and Instrumentation Conference 10708: Millimeter, Submillimeter, and Far-Infrared Detectors and Instrumentation for Astronomy XI, June 2018

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., Comment: 9 pages, 8 figures: conference proceedings, 53rd Rencontres de Moriond, Cosmology 2018

Silicon optics with wide bandwidth anti-reflection (AR) coatings, made of multi-layer textured silicon surfaces, are developed for millimeter and submillimeter wavelengths. Single and double layer AR coatings were designed for an optimal transmission centered on 250 GHz, and fabricated using the DRIE (Deep Reaction Ion Etching) technique. Tests of high resistivity silicon wafers with single-layer coatings between 75 GHz and 330 GHz are presented and compared with the simulations., Comment: A more detailed article has been written and is now available in arXiv (arXiv:1803.05168) and published in Applied Optics (doi: 10.1364/AO.57.005196)

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., Comment: 6 pages, 2 figures

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., Comment: 11 pages, 10 figures. To be published in Proc. SPIE. Presented at SPIE Astronomical Telescopes and Instrumentation Conference 9914: Millimeter, Submillimeter, and Far-Infrared Detectors and Instrumentation for Astronomy VIII, June 2016

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., Comment: 17 pages, 6 figures. To be published in Proc. SPIE. Presented at SPIE Astronomical Telescopes and Instrumentation Conference 9914: Millimeter, Submillimeter, and Far-Infrared Detectors and Instrumentation for Astronomy VIII, June 2016

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., Comment: 17 pages, 9 figures. Presented at SPIE Astronomical Telescopes and Instrumentation Conference 9914: Millimeter, Submillimeter, and Far-Infrared Detectors and Instrumentation for Astronomy VIII, June 2016

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., Comment: 12 pages, 8 figures

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}$., Comment: 21 pages, 23 figures, minor updates to reflect accepted version

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.., Comment: 7 pages, LTD-16 proceedings

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., Comment: 17 pages, 18 figures, minor updates to reflect accepted version, data and figures available for download at http://bicepkeck.org/