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.

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.

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

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.

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 provide an overview of the science case, instrument configuration and project plan for the next-generation ground-based cosmic microwave background experiment CMB-S4, for consideration by the 2020 Decadal Survey.

We present the science case, reference design, and project plan for the Stage-4 ground-based cosmic microwave background experiment CMB-S4.

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

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

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

The compelling science case for the observation of B-mode polarization in the cosmic microwave background (CMB) is driving the CMB community to expand the observed sky fraction, either by extending survey sizes or by deploying receivers to potential new northern sites. For ground-based CMB instruments, poorly-mixed atmospheric water vapor constitutes the primary source of short-term sky noise. This results in short-timescale brightness fluctuations, which must be rejected by some form of modulation. To maximize the sensitivity of ground-based CMB observations, it is useful to understand the effects of atmospheric water vapor over timescales and angular scales relevant for CMB polarization measurements. To this end, we have undertaken a campaign to perform a coordinated characterization of current and potential future observing sites using scanning 183 GHz water vapor radiometers (WVRs). So far, we have deployed two identical WVR units; one at the South Pole, Antarctica, and the other at Summit Station, Greenland. The former site has a long heritage of ground-based CMB observations and is the current location of the Bicep/Keck Array telescopes as well as the South Pole Telescope. The latter site, though less well characterized, is under consideration as a northern-hemisphere location for future CMB receivers. Data collection from this campaign began in January 2016 at South Pole and July 2016 at Summit Station. Data analysis is ongoing to reduce the data to a single spatial and temporal statistic that can be used for one-to-one site comparison.

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.

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

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

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.

The Atacama Millimeter/submillimeter Array (ALMA) Phasing Project (APP) has developed and deployed the hardware and software necessary to coherently sum the signals of individual ALMA antennas and record the aggregate sum in Very Long Baseline Interferometry (VLBI) Data Exchange Format. These beamforming capabilities allow the ALMA array to collectively function as the equivalent of a single large aperture and participate in global VLBI arrays. The inclusion of phased ALMA in current VLBI networks operating at (sub)millimeter wavelengths provides an order of magnitude improvement in sensitivity, as well as enhancements in u-v coverage and north-south angular resolution. The availability of a phased ALMA enables a wide range of new ultra-high angular resolution science applications, including the resolution of supermassive black holes on event horizon scales and studies of the launch and collimation of astrophysical jets. It also provides a high-sensitivity aperture that may be used for investigations such as pulsar searches at high frequencies. This paper provides an overview of the ALMA Phasing System design, implementation, and performance characteristics.

This paper presents the first detailed investigation of the characteristics of mm/submm phase fluctuation and phase correction methods obtained using ALMA with baseline lengths up to ~15 km. Most of the spatial structure functions (SSFs) show that the phase fluctuation increases as a function of baseline length, with a power-law slope of ~0.6. In many cases, we find that the slope becomes shallower (average of ~0.2-0.3) at baseline lengths longer than ~1 km, namely showing a turn-over in SSF. The phase correction method using water vapor radiometers (WVRs) works well, especially for the cases where PWV >1 mm, which reduces the degree of phase fluctuations by a factor of two in many cases. However, phase fluctuations still remain after the WVR phase correction, suggesting the existence of other turbulent constituent that cause the phase fluctuation. This is supported by occasional SSFs that do not exhibit any turn-over; these are only seen when the PWV is low or after WVR phase correction. This means that the phase fluctuation caused by this turbulent constituent is inherently smaller than that caused by water vapor. Since there is no turn-over in the SSF up to the maximum baseline length of ~15 km, this turbulent constituent must have scale height of 10 km or more, and thus cannot be water vapor, whose scale height is around 1 km. This large scale height turbulent constituent is likely to be water ice or a dry component. Excess path length fluctuation after the WVR phase correction at a baseline length of 10 km is large (>200 micron), which is significant for high frequency (>450 GHz or <700 micron) observations. These results suggest the need for an additional phase correction method, such as fast switching, in addition to the WVR phase correction. We simulated the fast switching, and the result suggests that it works well, with shorter cycle times linearly improving the coherence.

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

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.

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.

The star HR 8799 hosts one of the largest known debris discs and at least four giant planets. Previous observations have found evidence for a warm belt within the orbits of the planets, a cold planetesimal belt beyond their orbits and a halo of small grains. With the infrared data, it is hard to distinguish the planetesimal belt emission from that of the grains in the halo. With this in mind, the system has been observed with ALMA in band 6 (1.34 mm) using a compact array format. These observations allow the inner edge of the planetesimal belt to be resolved for the first time. A radial distribution of dust grains is fitted to the data using an MCMC method. The disc is best fit by a broad ring between $145^{+12}_{-12}$ AU and $429^{+37}_{-32}$ AU at an inclination of $40^{+5}_{-6}${\deg} and a position angle of $51^{+8}_{-8}${\deg}. A disc edge at ~145 AU is too far out to be explained simply by interactions with planet b, requiring either a more complicated dynamical history or an extra planet beyond the orbit of planet b.

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

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.

A major goal of the Atacama Large Millimeter/submillimeter Array (ALMA) is to make accurate images with resolutions of tens of milliarcseconds, which at submillimeter (submm) wavelengths requires baselines up to ~15 km. To develop and test this capability, a Long Baseline Campaign (LBC) was carried out from September to late November 2014, culminating in end-to-end observations, calibrations, and imaging of selected Science Verification (SV) targets. This paper presents an overview of the campaign and its main results, including an investigation of the short-term coherence properties and systematic phase errors over the long baselines at the ALMA site, a summary of the SV targets and observations, and recommendations for science observing strategies at long baselines. Deep ALMA images of the quasar 3C138 at 97 and 241 GHz are also compared to VLA 43 GHz results, demonstrating an agreement at a level of a few percent. As a result of the extensive program of LBC testing, the highly successful SV imaging at long baselines achieved angular resolutions as fine as 19 mas at ~350 GHz. Observing with ALMA on baselines of up to 15 km is now possible, and opens up new parameter space for submm astronomy.

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

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

We present Atacama Large Millimeter/submillimeter Array (ALMA) observations from the 2014 Long Baseline Campaign in dust continuum and spectral line emission from the HL Tau region. The continuum images at wavelengths of 2.9, 1.3, and 0.87 mm have unprecedented angular resolutions of 0.075 arcseconds (10 AU) to 0.025 arcseconds (3.5 AU), revealing an astonishing level of detail in the circumstellar disk surrounding the young solar analogue HL Tau, with a pattern of bright and dark rings observed at all wavelengths. By fitting ellipses to the most distinct rings, we measure precise values for the disk inclination (46.72pm0.05 degrees) and position angle (+138.02pm0.07 degrees). We obtain a high-fidelity image of the 1.0 mm spectral index ($\alpha$), which ranges from $\alpha\sim2.0$ in the optically-thick central peak and two brightest rings, increasing to 2.3-3.0 in the dark rings. The dark rings are not devoid of emission, we estimate a grain emissivity index of 0.8 for the innermost dark ring and lower for subsequent dark rings, consistent with some degree of grain growth and evolution. Additional clues that the rings arise from planet formation include an increase in their central offsets with radius and the presence of numerous orbital resonances. At a resolution of 35 AU, we resolve the molecular component of the disk in HCO+ (1-0) which exhibits a pattern over LSR velocities from 2-12 km/s consistent with Keplerian motion around a ~1.3 solar mass star, although complicated by absorption at low blue-shifted velocities. We also serendipitously detect and resolve the nearby protostars XZ Tau (A/B) and LkHa358 at 2.9 mm.

We present Atacama Large Millimeter/submillimeter Array (ALMA) 1.3 mm continuum images of the asteroid 3 Juno obtained with an angular resolution of 0.042 arcseconds (60 km at 1.97 AU). The data were obtained over a single 4.4 hr interval, which covers 60% of the 7.2 hr rotation period, approximately centered on local transit. A sequence of ten consecutive images reveals continuous changes in the asteroid's profile and apparent shape, in good agreement with the sky projection of the three-dimensional model of the Database of Asteroid Models from Inversion Techniques. We measure a geometric mean diameter of 259pm4 km, in good agreement with past estimates from a variety of techniques and wavelengths. Due to the viewing angle and inclination of the rotational pole, the southern hemisphere dominates all of the images. The median peak brightness temperature is 215pm13 K, while the median over the whole surface is 197pm15 K. With the unprecedented resolution of ALMA, we find that the brightness temperature varies across the surface with higher values correlated to the subsolar point and afternoon areas, and lower values beyond the evening terminator. The dominance of the subsolar point is accentuated in the final four images, suggesting a reduction in the thermal inertia of the regolith at the corresponding longitudes, which are possibly correlated to the location of the putative large impact crater. These results demonstrate ALMA's potential to resolve thermal emission from the surface of main belt asteroids, and to measure accurately their position, geometric shape, rotational period, and soil characteristics.

We present initial results of very high resolution Atacama Large Millimeter/submillimeter Array (ALMA) observations of the $z$=3.042 gravitationally lensed galaxy HATLAS J090311.6+003906 (SDP.81). These observations were carried out using a very extended configuration as part of Science Verification for the 2014 ALMA Long Baseline Campaign, with baselines of up to 15 km. We present continuum imaging at 151, 236 and 290 GHz, at unprecedented angular resolutions as fine as 23 milliarcseconds (mas), corresponding to an un-magnified spatial scale of ~180 pc at z=3.042. The ALMA images clearly show two main gravitational arc components of an Einstein ring, with emission tracing a radius of ~1.5". We also present imaging of CO(10-9), CO(8-7), CO(5-4) and H2O line emission. The CO emission, at an angular resolution of ~170 mas, is found to broadly trace the gravitational arc structures but with differing morphologies between the CO transitions and compared to the dust continuum. Our detection of H2O line emission, using only the shortest baselines, provides the most resolved detection to date of thermal H2O emission in an extragalactic source. The ALMA continuum and spectral line fluxes are consistent with previous Plateau de Bure Interferometer and Submillimeter Array observations despite the impressive increase in angular resolution. Finally, we detect weak unresolved continuum emission from a position that is spatially coincident with the center of the lens, with a spectral index that is consistent with emission from the core of the foreground lensing galaxy.

Radio and mm-wavelength observations of Sagittarius A* (Sgr A*), the radio source associated with the supermassive black hole at the center of our Galaxy, show that it behaves as a partially self-absorbed synchrotron-emitting source. The measured size of Sgr A* shows that the mm-wavelength emission comes from a small region and consists of the inner accretion flow and a possible collimated outflow. Existing observations of Sgr A* have revealed a time lag between light curves at 43 GHz and 22 GHz, which is consistent with a rapidly expanding plasma flow and supports the presence of a collimated outflow from the environment of an accreting black hole. Here we wish to measure simultaneous frequency-dependent time lags in the light curves of Sgr A* across a broad frequency range to constrain direction and speed of the radio-emitting plasma in the vicinity of the black hole. Light curves of Sgr A* were taken in May 2012 using ALMA at 100 GHz using the VLA at 48, 39, 37, 27, 25.5, and 19 GHz. As a result of elevation limits and the longitude difference between the stations, the usable overlap in the light curves is approximately four hours. Although Sgr A* was in a relatively quiet phase, the high sensitivity of ALMA and the VLA allowed us to detect and fit maxima of an observed minor flare where flux density varied by ~10%. The fitted times of flux density maxima at frequencies from 100 GHz to 19 GHz, as well as a cross-correlation analysis, reveal a simple frequency-dependent time lag relation where maxima at higher frequencies lead those at lower frequencies. Taking the observed size-frequency relation of Sgr A* into account, these time lags suggest a moderately relativistic (lower estimates: 0.5c for two-sided, 0.77c for one-sided) collimated outflow.

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.

We present the phase characteristics study of the Atacama Large Millimeter/submillimeter Array (ALMA) long (up to 3 km) baseline, which is the longest baseline tested so far using ALMA. The data consist of long time-scale (10 - 20 minutes) measurements on a strong point source (i.e., bright quasar) at various frequency bands (bands 3, 6, and 7, which correspond to the frequencies of about 88 GHz, 232 GHz, and 336 GHz). Water vapor radiometer (WVR) phase correction works well even at long baselines, and the efficiency is better at higher PWV (>1 mm) condition, consistent with the past studies. We calculate the spatial structure function of phase fluctuation, and display that the phase fluctuation (i.e., rms phase) increases as a function of baseline length, and some data sets show turn-over around several hundred meters to 1 km and being almost constant at longer baselines. This is the first millimeter/submillimeter structure function at this long baseline length, and to show the turn-over of the structure function. Furthermore, the observation of the turn-over indicates that even if the ALMA baseline length extends to the planned longest baseline of 15 km, fringes will be detected at a similar rms phase fluctuation as that at a few km baseline lengths. We also calculate the coherence time using the 3 km baseline data, and the results indicate that the coherence time for band 3 is longer than 400 seconds in most of the data (both in the raw and WVR-corrected data). For bands 6 and 7, WVR-corrected data have about twice longer coherence time, but it is better to use fast switching method to avoid the coherence loss.

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

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

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

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

Atacama Large Millimeter/submillimeter Array (ALMA) will be the world largest mm/submm interferometer, and currently the Early Science is ongoing, together with the commissioning and science verification (CSV). Here we present a study of the temporal phase stability of the entire ALMA system from antennas to the correlator. We verified the temporal phase stability of ALMA using data, taken during the last two years of CSV activities. The data consist of integrations on strong point sources (i.e., bright quasars) at various frequency bands, and at various baseline lengths (up to 600 m). From the observations of strong quasars for a long time (from a few tens of minutes, up to an hour), we derived the 2-point Allan Standard Deviation after the atmospheric phase correction using the 183 GHz Water Vapor Radiometer (WVR) installed in each 12 m antenna, and confirmed that the phase stability of all the baselines reached the ALMA specification. Since we applied the WVR phase correction to all the data mentioned above, we also studied the effectiveness of the WVR phase correction at various frequencies, baseline lengths, and weather conditions. The phase stability often improves a factor of 2 - 3 after the correction, and sometimes a factor of 7 improvement can be obtained. However, the corrected data still displays an increasing phase fluctuation as a function of baseline length, suggesting that the dry component (e.g., N2 and O2) in the atmosphere also contributes the phase fluctuation in the data, although the imperfection of the WVR phase correction cannot be ruled out at this moment.

We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of the [CII] 157.7micron fine structure line and thermal dust continuum emission from a pair of gas-rich galaxies at z=4.7, BR1202-0725. This system consists of a luminous quasar host galaxy and a bright submm galaxy (SMG), while a fainter star-forming galaxy is also spatially coincident within a 4" (25 kpc) region. All three galaxies are detected in the submm continuum, indicating FIR luminosities in excess of 10^13 Lsun for the two most luminous objects. The SMG and the quasar host galaxy are both detected in [CII] line emission with luminosities, L([CII]) = (10.0 +/- 1.5)x10^9 Lsun and L([CII]) = (6.5+/-1.0)x10^9 Lsun, respectively. We estimate a luminosity ratio, L([CII])/L(FIR) = (8.3+/-1.2)x10^-4 for the starburst SMG to the North, and L([CII])/L(FIR) = (2.5+/-0.4)x10^-4 for the quasar host galaxy, in agreement with previous high-redshift studies that suggest lower [CII]-to-FIR luminosity ratios in quasars than in starburst galaxies. The third fainter object with a flux density, S(340GHz) = 1.9+/-0.3 mJy, is coincident with a Ly-Alpha emitter and is detected in HST ACS F775W and F814W images but has no clear counterpart in the H-band. Even if this third companion does not lie at a similar redshift to BR1202-0725, the quasar and the SMG represent an overdensity of massive, infrared luminous star-forming galaxies within 1.3 Gyr of the Big Bang.

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

A large sub-mm survey with Herschel will enable many exciting science opportunities, especially in an era of wide-field optical and radio surveys and high resolution cosmic microwave background experiments. The Herschel-SPIRE Legacy Survey (HSLS), will lead to imaging data over 4000 sq. degrees at 250, 350, and 500 micron. Major Goals of HSLS are: (a) produce a catalog of 2.5 to 3 million galaxies down to 26, 27 and 33 mJy (50% completeness; 5 sigma confusion noise) at 250, 350 and 500 micron, respectively, in the southern hemisphere (3000 sq. degrees) and in an equatorial strip (1000 sq. degrees), areas which have extensive multi-wavelength coverage and are easily accessible from ALMA. Two thirds of the of the sources are expected to be at z > 1, one third at z > 2 and about a 1000 at z > 5. (b) Remove point source confusion in secondary anisotropy studies with Planck and ground-based CMB data. (c) Find at least 1200 strongly lensed bright sub-mm sources leading to a 2% test of general relativity. (d) Identify 200 proto-cluster regions at z of 2 and perform an unbiased study of the environmental dependence of star formation. (e) Perform an unbiased survey for star formation and dust at high Galactic latitude and make a census of debris disks and dust around AGB stars and white dwarfs.

We present a method of cross-calibrating the polarization angle of a polarimeter using BICEP Galactic observations. \bicep\ was a ground based experiment using an array of 49 pairs of polarization sensitive bolometers observing from the geographic South Pole at 100 and 150 GHz. The BICEP polarimeter is calibrated to +/-0.01 in cross-polarization and less than +/-0.7 degrees in absolute polarization orientation. BICEP observed the temperature and polarization of the Galactic plane (R.A= 100 degrees ~ 270 degrees and Dec. = -67 degrees ~ -48 degrees). We show that the statistical error in the 100 GHz BICEP Galaxy map can constrain the polarization angle offset of WMAP Wband to 0.6 degrees +\- 1.4 degrees. The expected 1 sigma errors on the polarization angle cross-calibration for Planck or EPIC are 1.3 degrees and 0.3 degrees at 100 and 150 GHz, respectively. We also discuss the expected improvement of the BICEP Galactic field observations with forthcoming BICEP2 and Keck observations.

Background Imaging of Cosmic Extragalactic Polarization (BICEP) is a bolometric polarimeter designed to measure the inflationary B-mode polarization of the cosmic microwave background (CMB) at degree angular scales. During three seasons of observing at the South Pole (2006 through 2008), BICEP mapped ~2% of the sky chosen to be uniquely clean of polarized foreground emission. Here we present initial results derived from a subset of the data acquired during the first two years. We present maps of temperature, Stokes Q and U, E and B modes, and associated angular power spectra. We demonstrate that the polarization data are self-consistent by performing a series of jackknife tests. We study potential systematic errors in detail and show that they are sub-dominant to the statistical errors. We measure the E-mode angular power spectrum with high precision at 21 < ell < 335, detecting for the first time the peak expected at ell ~ 140. The measured E-mode spectrum is consistent with expectations from a LCDM model, and the B-mode spectrum is consistent with zero. The tensor-to-scalar ratio derived from the B-mode spectrum is r = 0.03+0.31-0.26, or r < 0.72 at 95% confidence, the first meaningful constraint on the inflationary gravitational wave background to come directly from CMB B-mode polarization.

The BICEP experiment was designed specifically to search for the signature of inflationary gravitational waves in the polarization of the cosmic microwave background (CMB). Using a novel small-aperture refractor and 49 pairs of polarization-sensitive bolometers, BICEP has completed 3 years of successful observations at the South Pole beginning in 2006 February. To constrain the amplitude of the inflationary B-mode polarization, which is expected to be at least 7 orders of magnitude fainter than the 3 K CMB intensity, precise control of systematic effects is essential. This paper describes the characterization of potential systematic errors for the BICEP experiment, supplementing a companion paper on the initial cosmological results. Using the analysis pipelines for the experiment, we have simulated the impact of systematic errors on the B-mode polarization measurement. Guided by these simulations, we have established benchmarks for the characterization of critical instrumental properties including bolometer relative gains, beam mismatch, polarization orientation, telescope pointing, sidelobes, thermal stability, and timestream noise model. A comparison of the benchmarks with the measured values shows that we have characterized the instrument adequately to ensure that systematic errors do not limit BICEP's 2-year results, and identifies which future refinements are likely necessary to probe inflationary B-mode polarization down to levels below a tensor-to-scalar ratio r = 0.1.

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

We present new measurements of the cosmic microwave background (CMB) polarization from the final season of the Cosmic Anisotropy Polarization MAPper (CAPMAP). The data set was obtained in winter 2004-2005 with the 7 m antenna in Crawford Hill, New Jersey, from 12 W-band (84-100 GHz) and 4 Q-band (36-45 GHz) correlation polarimeters with 3.3' and 6.5' beamsizes, respectively. After selection criteria were applied, 956 (939) hours of data survived for analysis of W-band (Q-band) data. Two independent and complementary pipelines produced results in excellent agreement with each other. A broad suite of null tests as well as extensive simulations showed that systematic errors were minimal, and a comparison of the W-band and Q-band sky maps revealed no contamination from galactic foregrounds. We report the E-mode and B-mode power spectra in 7 bands in the range 200 < l < 3000, extending the range of previous measurements to higher l. The E-mode spectrum, which is detected at 11 sigma significance, is in agreement with cosmological predictions and with previous work at other frequencies and angular resolutions. The BB power spectrum provides one of the best limits to date on B-mode power at 4.8 uK^2 (95% confidence).

South Pole Station offers a unique combination of high, dry, stable conditions and well-developed support facilities. Over the past 20 years, a sequence of increasingly sophisticated CMB experiments at Pole have built on the experience of early pioneering efforts, producing a number of landmark contributions to the field. Telescopes at the South Pole were among the first to make repeated detections of degree-scale CMB temperature anisotropy and to map out the harmonic structure of its acoustic peaks. More recent achievements include the first detection of polarization of the CMB and the most precise measurements of the temperature power spectrum at small angular scales. New CMB telescopes at the South Pole are now making ultra-deep observations of the large-scale polarization of the CMB and of its secondary temperature anisotropies on arcminute scales. These two observing goals represent the current frontiers of CMB research, focused on constraining Inflation and the nature of Dark Energy. The South Pole now hosts an array of CMB observing platforms covering a wide range of angular scales and supporting very long integration times on the cleanest sky available, and thus should play an increasing role in pushing these frontiers of CMB research.

The Robinson Telescope (BICEP) is a ground-based millimeter-wave bolometric array designed to study the polarization of the cosmic microwave background radiation (CMB) and galactic foreground emission. Such measurements probe the energy scale of the inflationary epoch, tighten constraints on cosmological parameters, and verify our current understanding of CMB physics. Robinson consists of a 250-mm aperture refractive telescope that provides an instantaneous field-of-view of 17 degrees with angular resolution of 55 and 37 arcminutes at 100 GHz and 150 GHz, respectively. Forty-nine pair of polarization-sensitive bolometers are cooled to 250 mK using a 4He/3He/3He sorption fridge system, and coupled to incoming radiation via corrugated feed horns. The all-refractive optics is cooled to 4 K to minimize polarization systematics and instrument loading. The fully steerable 3-axis mount is capable of continuous boresight rotation or azimuth scanning at speeds up to 5 deg/s. Robinson has begun its first season of observation at the South Pole. Given the measured performance of the instrument along with the excellent observing environment, Robinson will measure the E-mode polarization with high sensitivity, and probe for the B-modes to unprecedented depths. In this paper we discuss aspects of the instrument design and their scientific motivations, scanning and operational strategies, and the results of initial testing and observations.

The Princeton IQU Experiment (PIQUE) and the Cosmic Anisotropy Polarization MAPper (CAPMAP) are experiments designed to measure the polarization of the Cosmic Microwave Background (CMB) on sub-degree scales in an area within 1 degree of the North Celestial Pole using heterodyne correlation polarimeters and off-axis telescopes located in central New Jersey. PIQUE produced the tightest limit on the CMB polarization prior to its detection by DASI, while CAPMAP has recently detected polarization at l~1000. The experimental methods and instrumentation for these two projects are described in detail with emphasis on the particular challenges involved in measuring the tiny polarized component of the CMB.

Polarization results from the Cosmic Anisotropy Polarization MAPper (CAPMAP) experiment are reported. These are based upon 433 hours, after cuts, observing a 2 square degree patch around the North Celestial Pole (NCP) with four 90 GHz correlation polarimeters coupled to optics defining $4\arcmin$ beams. The E-mode flat bandpower anisotropy within $\ell=940^{+330}_{-300}$ is measured as 66$^{+69}_{-29} \mu$K$^2$; the 95% Confidence level upper limit for B-mode power within $\ell=1050^{+590}_{-520}$ is measured as 38 $\mu$K$^2$.

I describe here the new CAPMAP (Cosmic Anisotropy Polarization MAPer) instrument, which performed its first season of observing between February and April 2003 from the Crawford Hill 7-meter antenna in NJ. CAPMAP is based on the design for the PIQUE instrument, but has better sensitivity and superior angular resolution.

We update the limit from the 90 GHz PIQUE ground-based polarimeter on the magnitude of any polarized anisotropy of the cosmic microwave radiation. With a second year of data, we have now limited both Q and U on a ring of 1 degree radius. The window functions are broad: for E-mode polarization, the effective l is <l> = 191 +143 -132. We find that the E-mode signal can be no greater than 8.4 microK (95% CL), assuming no B-mode polarization. Limits on a possible B-mode signal are also presented.

We compute upper limits on CMB cross-polarization by cross-correlating the PIQUE and Saskatoon experiments. We also discuss theoretical and practical issues relevant to measuring cross-polarization and illustrate them with simulations of the upcoming BOOMERanG 2002 experiment. We present a method that separates all six polarization power spectra (TT, EE, BB, TE, TB, EB) without any other "leakage" than the familiar EE-BB mixing caused by incomplete sky coverage. Since E and B get mixed, one might expect leakage between TE and TB, between EE and EB and between BB and EB - our method eliminates this by preserving the parity symmetry under which TB and EB are odd and the other four power spectra are even.

A ground-based polarimeter, PIQUE, operating at 90 GHz has set a new limit on the magnitude of any polarized anisotropy in the cosmic microwave background. The combination of the scan strategy and full width half maximum beam of 0.235 degrees gives broad window functions with average multipoles, l = 211+294-146 and l = 212+229-135 for the E- and B-mode window functions, respectively. A joint likelihood analysis yields simultaneous 95% confidence level flat band power limits of 14 and 13 microkelvin on the amplitudes of the E- and B-mode angular power spectra, respectively. Assuming no B-modes, a 95% confidence limit of 10 microkelvin is placed on the amplitude of the E-mode angular power spectrum alone.