We present two prescriptions for broadband (~77 - 252 GHz), millimeter-wave antireflection coatings for cryogenic, sintered polycrystalline aluminum oxide optics: one for large-format (700 mm diameter) planar and plano-convex elements, the other for densely packed arrays of quasi-optical elements, in our case 5 mm diameter half-spheres (called "lenslets"). The coatings comprise three layers of commercially-available, polytetrafluoroethylene-based, dielectric sheet material. The lenslet coating is molded to fit the 150 mm diameter arrays directly while the large-diameter lenses are coated using a tiled approach. We review the fabrication processes for both prescriptions then discuss laboratory measurements of their transmittance and reflectance. In addition, we present the inferred refractive indices and loss tangents for the coating materials and the aluminum oxide substrate. We find that at 150 GHz and 300 K the large-format coating sample achieves (97 +/- 2)% transmittance and the lenslet coating sample achieves (94 +/- 3)% transmittance.

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.

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

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.

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.

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.

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

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

As one of the most massive Milky Way satellites, the Sagittarius dwarf galaxy has played an important role in shaping the Galactic disk and stellar halo morphologies. The disruption of Sagittarius over several close-in passages has populated the halo of our Galaxy with large-scale tidal streams and offers a unique diagnostic tool for measuring its gravitational potential. Here we test different progenitor mass models for the Milky Way and Sagittarius by modeling the full infall of the satellite. We constrain the mass of the Galaxy based on the observed orbital parameters and multiple tidal streams of Sagittarius. Our semi-analytic modeling of the orbital dynamics agrees with full $N$-body simulations, and favors low values for the Milky Way mass, $\lesssim 10^{12}M_\odot$. This conclusion eases the tension between $\Lambda$CDM and the observed parameters of the Milky Way satellites.

We report the detection of spatially distinct stellar density features near the apocenters of the Sagittarius (Sgr) stream's main leading and trailing arm. These features are clearly visible in a high-fidelity stellar halo map that is based on RR Lyrae from Pan-STARRS1: there is a plume of stars 10 kpc beyond the apocenter of the leading arm, and there is a "spur" extending to 130 kpc, almost 30 kpc beyond the previously detected apocenter of the trailing arm. Such apocenter substructure is qualitatively expected in any Sgr stream model, as stars stripped from the progenitor at different pericenter passages become spatially separated there. The morphology of these new Sgr stream substructures could provide much-needed new clues and constraints for modeling the Sgr system, including the level of dynamical friction that Sgr has experienced. We also report the discovery of a new, presumably unrelated halo substructure at 80 kpc from the Sun and $10^\circ$ from the Sgr orbital plane, which we dub the Outer Virgo Overdensity.

The extensive span of the Sagittarius (Sgr) stream makes it a promising tool for studying the Milky Way gravitational potential. Characterizing its stellar kinematics can constrain halo properties and provide a benchmark for the Cold Dark Matter galaxy formation paradigm. Accurate models of the disruption dynamics of the Sgr progenitor are necessary to employ this tool. Using a combination of analytic modeling and N-body simulations, we build a new model of the Sgr orbit and resulting stellar stream. In contrast to previous models, we simulate the full infall trajectory of the Sgr progenitor from the time it first crossed the Milky Way virial radius 8 Gyr ago. An exploration of the parameter space of initial phase-space conditions yields tight constraints on the angular momentum of the Sgr progenitor. Our best-fit model is the first to reproduce accurately existing data on the 3D positions and radial velocities of the debris detected 100 kpc away in the MW halo. In addition to replicating the mapped stream, the simulation also predicts the existence of several arms of the Sgr stream extending to hundreds of kiloparsecs. The two most distant stars known in the Milky Way halo coincide with the predicted structure. Additional stars in the newly predicted arms can be found with future data from the Large Synoptic Survey Telescope. Detecting a statistical sample of stars in the most distant Sgr arms would provide an opportunity to constrain the Milky Way potential out to unprecedented Galactocentric radii.

We present the first high-angular resolution study of the MonR2 star-forming complex carried out with the Submillimeter Array at (sub-)millimeter wavelengths. We image the continuum and molecular line emission toward the young stellar objects in MonR2 at 0.85mm and 1.3mm, with resolutions ranging from 0.5" to ~3". While free-free emission dominates the IRS1 and IRS2 continuum, dust thermal emission prevails for IRS3 and IRS5, giving envelope masses of ~0.1-0.3 M_Sun. IRS5 splits into at least two sub-arcsecond scale sources, IRS5B and the more massive IRS5A. Our 12CO(2-1) images reveal 11 previously unknown molecular outflows in the MonR2 clump. Comparing these outflows with known IR sources in the IRS5 and IRS3 subclusters allows for tentative identification of driving stars. Line images of molecular species such as CH3CN or CH3OH show that, besides IRS3 (a well-known hot molecular core), IRS5 is also a chemically active source in the region. The gas excitation temperature derived from CH3CN lines toward IRS5 is 144 \pm 15 K, indicating a deeply embedded protostar at the hot-core evolutionary stage. SED fitting of IRS5 gives a mass of ~7 M_Sun and a luminosity of 300 L_Sun for the central source. The derived physical properties of the CO outflows suggest that they contribute to the turbulent support of the MonR2 complex and to the gas velocity dispersion in the clump's center. The detection of a large number of CO outflows widespread across the region supports the competitive accretion scenario as origin of the MonR2 star cluster.

The unusual morphologies of the Andromeda spiral galaxy (M31) and its dwarf companion M32 have been characterized observationally in great detail. The two galaxies' apparent proximity suggests that Andromeda's prominent star-forming ring as well as M32's compact elliptical structure may result from a recent collision. Here we present the first self-consistent model of the M31-M32 interaction that simultaneously reproduces observed positions, velocities, and morphologies for both galaxies. Andromeda's spiral structure is resolved in unprecedented detail, showing that a rare head-on orbit is not necessary to match Andromeda's ring-like morphology. The passage of M32 through Andromeda's disk perturbs the disk velocity structure. We find tidal stripping of M32's stars to be inefficient during the interaction, suggesting that some cEs are intrinsically compact. Additionally, the orbital solution implies that M32 is currently closer to the Milky Way than models have typically assumed, a prediction that may be testable with upcoming observations.

We report the detection of a new radio recombination line (RRL) maser object toward the IRS2 source in the MonR2 ultracompact HII region. The continuum emission at 1.3mm and 0.85mm and the H30a and H26a lines were observed with the Submillimeter Array (SMA) at angular resolutions of about 0.5"-3". The SMA observations show that the MonR2-IRS2 source is very compact and remains unresolved at spatial scales <=400AU. Its continuum power spectrum at millimeter wavelengths is almost flat (alpha=-0.16, with S_nu proportional to nu^alpha), indicating that this source is dominated by optically thin free-free emission. The H30a and H26a RRL emission is also compact and peaks toward the position of the MonR2-IRS2 source. The measured RRL profiles are double-peaked with the H26a line showing a clear asymmetry in its spectrum. Since the derived line-to-continuum flux ratios (80 and 180kms-1 for H30a and H26a, respectively) exceed the LTE predictions, the RRLs toward MonR2-IRS2 are affected by maser amplification. The amplification factors are however smaller than those found toward the emission line star MWC349A, indicating that MonR2-IRS2 is a weakly amplified maser. Radiative transfer modelling of the RRL emission toward this source shows that the RRL masers arise from a dense and collimated jet embedded in a cylindrical ionized wind, oriented nearly along the direction of the line-of-sight. High-angular resolution observations at sub-millimeter wavelengths are needed to unveil weakly amplified RRL masers in very young massive stars.

We test competing models that aim at explaining the nature of stars in the Milky Way that are well away (|z|$\gtrsim$ 1kpc) from the midplane, the so-called thick disk: the stars may have gotten there through orbital migration, through satellite mergers and accretion, or through heating of pre-existing thin disk stars. Sales et al. (2009) proposed the eccentricity distribution of thick disk stars as a diagnostic to differentiate between these mechanisms. Drawing on SDSS DR7, we have assembled a sample of 34,223 G-dwarfs with 6-D phase-space information and metallicities, and have derived orbital eccentricities for them. Comparing the resulting eccentricity distributions, p(e|z), with the models, we find that: a) the observed p(e|z) is inconsistent with that predicted by orbital migration only, as there are more observed stars of high and of very low eccentricity; b) scenarios where the thick disk is made predominantly through abrupt heating of a pre-existing thin disk are also inconsistent, as they predict more high-eccentricity stars than observed; c) the observed p(e|z) fits well with a "gas-rich merger" scenario, where most thick disk stars were born from unsettled gas in situ.