Search

Measurements of B-mode polarization in the CMB sourced from primordial gravitational waves would provide information on the energy scale of inflation and its potential form. To achieve these goals, one must carefully characterize the Galactic foregrounds, which can be distinguished from the CMB by conducting measurements at multiple frequencies. BICEP Array is the latest-generation multi-frequency instrument of the BICEP/Keck program, which specifically targets degree-scale primordial B-modes in the CMB. In its final configuration, this telescope will consist of four small-aperture receivers, spanning frequency bands from 30 to 270 GHz. The 220/270 GHz receiver designed to characterize Galactic dust is currently undergoing commissioning at Stanford University and is scheduled to deploy to the South Pole during the 2024--2025 austral summer. Here, we will provide an overview of this high-frequency receiver and discuss the integration status and test results as it is being commissioned.

SPIDER is a balloon-borne instrument designed to map the cosmic microwave background at degree-angular scales in the presence of Galactic foregrounds. SPIDER has mapped a large sky area in the Southern Hemisphere using more than 2000 transition-edge sensors (TESs) during two NASA Long Duration Balloon flights above the Antarctic continent. During its first flight in January 2015, SPIDER observed in the 95 GHz and 150 GHz frequency bands, setting constraints on the B-mode signature of primordial gravitational waves. Its second flight in the 2022-23 season added new receivers at 280 GHz, each using an array of TESs coupled to the sky through feedhorns formed from stacks of silicon wafers. These receivers are optimized to produce deep maps of polarized Galactic dust emission over a large sky area, providing a unique data set with lasting value to the field. In this work, we describe the instrument's performance during SPIDER's second flight., Comment: Submitted to SPIE Astronomical Telescopes + Instrumentation 2024, JATIS

Using data from the first flight of SPIDER and from Planck HFI, we probe the properties of polarized emission from interstellar dust in the SPIDER observing region. Component separation algorithms operating in both the spatial and harmonic domains are applied to probe their consistency and to quantify modeling errors associated with their assumptions. Analyses spanning the full SPIDER region demonstrate that i) the spectral energy distribution of diffuse Galactic dust emission is broadly consistent with a modified-blackbody (MBB) model with a spectral index of $\beta_\mathrm{d}=1.45\pm0.05$ $(1.47\pm0.06)$ for $E$ ($B$)-mode polarization, slightly lower than that reported by Planck for the full sky; ii) its angular power spectrum is broadly consistent with a power law; and iii) there is no significant detection of line-of-sight decorrelation of the astrophysical polarization. The size of the SPIDER region further allows for a statistically meaningful analysis of the variation in foreground properties within it. Assuming a fixed dust temperature $T_\mathrm{d}=19.6$ K, an analysis of two independent sub-regions of that field results in inferred values of $\beta_\mathrm{d}=1.52\pm0.06$ and $\beta_\mathrm{d}=1.09\pm0.09$, which are inconsistent at the $3.9\,\sigma$ level. Furthermore, a joint analysis of SPIDER and Planck 217 and 353 GHz data within a subset of the SPIDER region is inconsistent with a simple MBB at more than $3\,\sigma$, assuming a common morphology of polarized dust emission over the full range of frequencies. These modeling uncertainties have a small--but non-negligible--impact on limits on the cosmological tensor-to-scalar ratio derived from the \spider dataset. The fidelity of the component separation approaches of future CMB polarization experiments may thus have a significant impact on their constraining power., Comment: 21 pages, 15 figures

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

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

Future detection of high-redshift gamma-ray bursts (GRBs) will be an important tool for studying the early Universe. Fast and accurate redshift estimation for detected GRBs is key for encouraging rapid follow-up observations by ground- and space-based telescopes. Low-redshift dusty interlopers pose the biggest challenge for GRB redshift estimation using broad photometric bands, as their high extinction can mimic a high-redshift GRB. To assess false alarms of high-redshift GRB photometric measurements, we simulate and fit a variety of GRBs using phozzy, a simulation code developed to estimate GRB photometric redshifts, and test the ability to distinguish between high- and low-redshift GRBs when using simultaneously observed photometric bands. We run the code with the wavelength bands and instrument parameters for the Photo-z Infrared Telescope (PIRT), an instrument designed for the Gamow mission concept. We explore various distributions of host galaxy extinction as a function of redshift, and their effect on the completeness and purity of a high-redshift GRB search with the PIRT. We find that for assumptions based on current observations, the completeness and purity range from $\sim 82$ to $88\%$ and from $\sim 84$ to $>99\%$, respectively. For the priors optimized to reduce false positives, only $\sim 0.6\%$ of low-redshift GRBs will be mistaken as a high-redshift one, corresponding to $\sim 1$ false alarm per 500 detected GRBs., Comment: 14 pages, 15 figures, accepted for publication in MNRAS

First identified in 2016 by JAXA's Akatsuki mission, the discontinuity/disruption is a recurrent wave observed to propagate during decades at the deeper clouds of Venus (47--56 km above the surface), while its absence at the clouds' top ($\sim$70 km) suggests that it dissipates at the upper clouds and contributes in the maintenance of the puzzling atmospheric superrotation of Venus through wave-mean flow interaction. Taking advantage of the campaign of ground-based observations undertaken in coordination with the Akatsuki mission since December 2021 until July 2022, we aimed to undertake the longest uninterrupted monitoring of the cloud discontinuity up to date to obtain a pioneering long-term characterization of its main properties and better constrain its recurrence and lifetime. The dayside upper, middle and nightside lower clouds were studied with images with suitable filters acquired by Akatsuki/UVI, amateur observers and NASA's IRTF/SpeX, respectively. Hundreds of images were inspected in search of manifestations of the discontinuity events and to measure key properties like its dimensions, orientation or rotation period. We succeeded in tracking the discontinuity at the middle clouds during 109 days without interruption. The discontinuity exhibited properties nearly identical to measurements in 2016 and 2020, with an orientation of $91^{\circ}\pm 8^{\circ}$, length/width of $4100\pm 800$ / $500\pm 100$ km and a rotation period of $5.11\pm 0.09$ days. Ultraviolet images during 13-14 June 2022 suggest that the discontinuity may have manifested at the top of the clouds during $\sim$21 hours as a result of an altitude change in the critical level for this wave due to slower zonal winds., Comment: 8 pages, 4 figures, 2 animated figures, 1 table

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

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

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

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

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

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

Experiments deploying large arrays of transition-edge sensors (TESs) often require a robust method to monitor gain variations with minimal loss of observing time. We propose a sensitive and non-intrusive method for monitoring variations in TES responsivity using small square waves applied to the TES bias. We construct an estimator for a TES's small-signal power response from its electrical response that is exact in the limit of strong electrothermal feedback. We discuss the application and validation of this method using flight data from SPIDER, a balloon-borne telescope that observes the polarization of the cosmic microwave background with more than 2000 TESs. This method may prove useful for future balloon- and space-based instruments, where observing time and ground control bandwidth are limited., Comment: 7 pages, 3 figures; Proceedings of the 19th International Workshop on Low Temperature Detectors (LTD19); Minor updates to match published version

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

The Photo-z InfraRed Telescope (PIRT) is an instrument on the Gamow Explorer, currently proposed for a NASA Astrophysics Medium Explorer. PIRT works in tandem with a companion wide-field instrument, the Lobster Eye X-ray Telescope (LEXT), that will identify x-ray transients likely to be associated with high redshift gamma-ray bursts (GRBs) or electromagnetic counterparts to gravitational wave (GW) events. After receiving an alert trigger from LEXT, the spacecraft will slew to center the PIRT field of view on the transient source. PIRT will then begin accumulating data simultaneously in five bands spanning 0.5 - 2.5 microns over a 10 arc-minute field of view. Each PIRT field will contain many hundreds of sources, only one of which is associated with the LEXT transient. PIRT will gather the necessary data in order to identify GRB sources with redshift $z > 6$, with an expected source localization better than 1 arcsec. A near real-time link to the ground will allow timely follow-up as a target of opportunity for large ground-based telescopes or the James Webb Space Telescope (JWST). PIRT will also allow localization and characterization of GW event counterparts. We discuss the instrument design, the on-board data processing approach, and the expected performance of the system., Comment: 10 pages, 5 figures

The Gamow Explorer will use Gamma Ray Bursts (GRBs) to: 1) probe the high redshift universe (z > 6) when the first stars were born, galaxies formed and Hydrogen was reionized; and 2) enable multi-messenger astrophysics by rapidly identifying Electro-Magnetic (IR/Optical/X-ray) counterparts to Gravitational Wave (GW) events. GRBs have been detected out to z ~ 9 and their afterglows are a bright beacon lasting a few days that can be used to observe the spectral fingerprints of the host galaxy and intergalactic medium to map the period of reionization and early metal enrichment. Gamow Explorer is optimized to quickly identify high-z events to trigger follow-up observations with JWST and large ground-based telescopes. A wide field of view Lobster Eye X-ray Telescope (LEXT) will search for GRBs and locate them with arc-minute precision. When a GRB is detected, the rapidly slewing spacecraft will point the 5 photometric channel Photo-z Infra-Red Telescope (PIRT) to identify high redshift (z > 6) long GRBs within 100s and send an alert within 1000s of the GRB trigger. An L2 orbit provides > 95% observing efficiency with pointing optimized for follow up by the James Webb Space Telescope (JWST) and ground observatories. The predicted Gamow Explorer high-z rate is >10 times that of the Neil Gehrels Swift Observatory. The instrument and mission capabilities also enable rapid identification of short GRBs and their afterglows associated with GW events. The Gamow Explorer will be proposed to the 2021 NASA MIDEX call and if approved, launched in 2028., Comment: 14 pages, 8 Figures

Modern CMB analysis pipelines regularly employ complex time-domain filters, beam models, masking, and other techniques during the production of sky maps and their corresponding angular power spectra. However, these processes can generate couplings between multipoles from the same spectrum and from different spectra, in addition to the typical power attenuation. Within the context of pseudo-$C_\ell$ based, MASTER-style analyses, the net effect of the time-domain filtering is commonly approximated by a multiplicative transfer function, $F_{\ell}$, that can fail to capture mode mixing and is dependent on the spectrum of the signal. To address these shortcomings, we have developed a simulation-based spectral correction approach that constructs a two-dimensional transfer matrix, $J_{\ell\ell'}$, which contains information about mode mixing in addition to mode attenuation. We demonstrate the application of this approach on data from the first flight of the SPIDER balloon-borne CMB experiment., Comment: 14 pages, 7 figures; updated to match published version

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

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

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

We present the XFaster analysis package. XFaster is a fast, iterative angular power spectrum estimator based on a diagonal approximation to the quadratic Fisher matrix estimator. XFaster uses Monte Carlo simulations to compute noise biases and filter transfer functions and is thus a hybrid of both Monte Carlo and quadratic estimator methods. In contrast to conventional pseudo-$C_\ell$ based methods, the algorithm described here requires a minimal number of simulations, and does not require them to be precisely representative of the data to estimate accurate covariance matrices for the bandpowers. The formalism works with polarization-sensitive observations and also data sets with identical, partially overlapping, or independent survey regions. The method was first implemented for the analysis of BOOMERanG data, and also used as part of the Planck analysis. Here, we describe the full, publicly available analysis package, written in Python, as developed for the analysis of data from the 2015 flight of the SPIDER instrument. The package includes extensions for self-consistently estimating null spectra and for estimating fits for Galactic foreground contributions. We show results from the extensive validation of XFaster using simulations, and its application to the SPIDER data set., Comment: 18 pages, 11 figures

We present the first linear polarization measurements from the 2015 long-duration balloon flight of SPIDER, an experiment designed to map the polarization of the cosmic microwave background (CMB) on degree angular scales. Results from these measurements include maps and angular power spectra from observations of 4.8% of the sky at 95 and 150 GHz, along with the results of internal consistency tests on these data. While the polarized CMB anisotropy from primordial density perturbations is the dominant signal in this region of sky, Galactic dust emission is also detected with high significance; Galactic synchrotron emission is found to be negligible in the SPIDER bands. We employ two independent foreground-removal techniques in order to explore the sensitivity of the cosmological result to the assumptions made by each. The primary method uses a dust template derived from Planck data to subtract the Galactic dust signal. A second approach, employing a joint analysis of SPIDER and Planck data in the harmonic domain, assumes a modified-blackbody model for the spectral energy distribution of the dust with no constraint on its spatial morphology. Using a likelihood that jointly samples the template amplitude and $r$ parameter space, we derive 95% upper limits on the primordial tensor-to-scalar ratio from Feldman-Cousins and Bayesian constructions, finding $r<0.11$ and $r<0.19$, respectively. Roughly half the uncertainty in $r$ derives from noise associated with the template subtraction. New data at 280 GHz from SPIDER's second flight will complement the Planck polarization maps, providing powerful measurements of the polarized Galactic dust emission., Comment: 29 pages, 13 figures

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

In this work we describe upgrades to the Spider balloon-borne telescope in preparation for its second flight, currently planned for December 2021. The Spider instrument is optimized to search for a primordial B-mode polarization signature in the cosmic microwave background at degree angular scales. During its first flight in 2015, Spider mapped ~10% of the sky at 95 and 150 GHz. The payload for the second Antarctic flight will incorporate three new 280 GHz receivers alongside three refurbished 95- and 150 GHz receivers from Spider's first flight. In this work we discuss the design and characterization of these new receivers, which employ over 1500 feedhorn-coupled transition-edge sensors. We describe pre-flight laboratory measurements of detector properties, and the optical performance of completed receivers. These receivers will map a wide area of the sky at 280 GHz, providing new information on polarized Galactic dust emission that will help to separate it from the cosmological signal., Comment: 13 pages, 8 figures; as published in the conference proceedings for SPIE Millimeter, Submillimeter, and Far-Infrared Detectors and Instrumentation for Astronomy X (2020)

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

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

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

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

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

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

Planetary-scale waves are thought to play a role in powering the yet-unexplained atmospheric superrotation of Venus. Puzzlingly, while Kelvin, Rossby and stationary waves manifest at the upper clouds (65--70 km), no planetary-scale waves or stationary patterns have been reported in the intervening level of the lower clouds (48--55 km), although the latter are probably Lee waves. Using observations by the Akatsuki orbiter and ground-based telescopes, we show that the lower clouds follow a regular cycle punctuated between 30$^{\circ}$N--40$^{\circ}$S by a sharp discontinuity or disruption with potential implications to Venus's general circulation and thermal structure. This disruption exhibits a westward rotation period of $\sim$4.9 days faster than winds at this level ($\sim$6-day period), alters clouds' properties and aerosols, and remains coherent during weeks. Past observations reveal its recurrent nature since at least 1983, and numerical simulations show that a nonlinear Kelvin wave reproduces many of its properties., Comment: 21 pages, 10 figures, 2 animated figures and 2 tables

Future mm-wave and sub-mm space missions will employ large arrays of multiplexed Transition Edge Sensor (TES) bolometers. Such instruments must contend with the high flux of cosmic rays beyond our atmosphere that induce "glitches" in bolometer data, which posed a challenge to data analysis from the Planck bolometers. Future instruments will face the additional challenges of shared substrate wafers and multiplexed readout wiring. In this work we explore the susceptibility of modern TES arrays to the cosmic ray environment of space using two data sets: the 2015 long-duration balloon flight of the SPIDER cosmic microwave background polarimeter, and a laboratory exposure of SPIDER flight hardware to radioactive sources. We find manageable glitch rates and short glitch durations, leading to minimal effect on SPIDER analysis. We constrain energy propagation within the substrate through a study of multi-detector coincidences, and give a preliminary look at pulse shapes in laboratory data., Comment: 9 pages, 6 figures, Proceedings of the 18th International Workshop on Low Temperature Detectors

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

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

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

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

$\gamma$-CoV$_{2}$O$_{6}$ is a quasi one-dimensional spin-$\frac{3}{2}$ magnet that possesses two distinct magnetic orders in the ground state with modulation vectors $k_\mathrm{1}$ = ($\frac{1}{2}$, 0, 0) and $k_\mathrm{2}$ = ($\frac{1}{4}$, 0, -$\frac{1}{4}$), respectively. Here, we use muon spin relaxation and rotation to reveal the thermodynamics of the magnetic phase separation in this compound. In the paramagnetic (PM) region, short-range correlated spin clusters emerge at $T_\mathrm{m}$ $\simeq$ 26 K at the $\it{partial}$ expense of the PM volume. Upon further cooling, we show that these emergent clusters become spatially coherent at $T_\mathrm{{N2}}$ = 7.5 K and eventually form the $k_\mathrm{2}$ order at $T^{\star}$ = 5.6 K, while the remaining PM spins are driven into the $k_\mathrm{1}$ state at $T_\mathrm{{N1}}$ = 6.6 K. These results stress magnetic microphase inhomogeneity as a thermodynamic precursor for the ground state phase separation in weakly coupled spin-$\frac{3}{2}$ chains.

Silicon and Mg in differentiated rocky bodies exhibit heavy isotope enrichments that have been attributed to evaporation of partially or entirely molten planetesimals. We evaluate the mechanisms of planetesimal evaporation in the early solar system and the conditions that controlled attendant isotope fractionations. Energy balance at the surface of a body accreted within ~1 Myr of CAI formation and heated from within by 26Al decay results in internal temperatures exceeding the silicate solidus, producing a transient magma ocean with a thin surface boundary layer of order < 1 meter that would be subject to foundering. Bodies that are massive enough to form magma oceans by radioisotope decay (ge 0.1%) can retain hot rock vapor even in the absence of ambient nebular gas. We find that a steady-state rock vapor forms within minutes to hours and results from a balance between rates of magma evaporation and atmospheric escape. Vapor pressure buildup adjacent to the surfaces of the evaporating magmas would have inevitably led to an approach to equilibrium isotope partitioning between the vapor phase and the silicate melt. Numerical simulations of this near-equilibrium evaporation process for a body with a radius of ~ 700 km yield a steady-state far-field vapor pressure corresponding to 95% saturation. Approaches to equilibrium isotope fractionation between vapor and melt should have been the norm during planet formation due to the formation of steady-state rock vapor atmospheres and/or the presence of protostellar gas. We model the Si and Mg isotopic composition of bulk Earth and show that the best fit is for a carbonaceous chondrite-like source material with about 12% loss of Mg and 15% loss of Si resulting from near-equilibrium evaporation into the solar protostellar disk of hydrogen gas on timescales of 10,000 to 100,000 years., Comment: 35 pages, 15 figures

During the years 2016 to 2018, the instruments Akatsuki/IR2 (JAXA) and IRTF/SpeX (NASA) acquired a large set of images at 1.74, 2.26 and 2.32 {\mu}m to study the nightside mid-to-lower clouds (48-60 km) of Venus. Here we summarize the rich variety of cloud morphologies apparent in these images: from frequent wave packets and billows caused by shear instabilities, to features reported decades ago like the circum-equatorial belts, bright blotches and equatorial troughs, and previously unseen features like dark spots, sharp dark streaks at mid latitudes and fully-developed vortices., Comment: 17 pages, 2 figures, 1 table

Pluto's atmospheric haze settles out rapidly compared with geological timescales. It needs to be accounted for as a surface material, distinct from Pluto's icy bedrock and from the volatile ices that migrate via sublimation and condensation on seasonal timescales. This paper explores how a steady supply of atmospheric haze might affect three distinct provinces on Pluto. We pose the question of why they each look so different from one another if the same haze material is settling out onto all of them. Cthulhu is a more ancient region with comparatively little present-day geological activity, where the haze appears to simply accumulate over time. Sputnik Planitia is a very active region where glacial convection, as well as sublimation and condensation rapidly refresh the surface, hiding recently deposited haze from view. Lowell Regio is a region of intermediate age featuring very distinct coloration from the rest of Pluto. Using a simple model haze particle as a colorant, we are not able to match the colors in both Lowell Regio and Cthulhu. To account for their distinct colors, we propose that after arrival at Pluto's surface, haze particles may be less inert than might be supposed from the low surface temperatures. They must either interact with local materials and environments to produce distinct products in different regions, or else the supply of haze must be non-uniform in time and/or location, such that different products are delivered to different places.

Context. Pluto's tenuous nitrogen (N2) atmosphere undergoes strong seasonal effects due to high obliquity and orbital eccentricity, and has been recently (July 2015) observed by the New Horizons spacecraft. Goals are (i) construct a well calibrated record of the seasonal evolution of surface pressure on Pluto and (ii) constrain the structure of the lower atmosphere using a central flash observed in 2015. Method: eleven stellar occultations by Pluto observed between 2002 and 2016 are used to retrieve atmospheric profiles (density, pressure, temperature) between $\sim$5 km and $\sim$380 km altitude levels (i.e. pressures from about 10 microbar to 10 nanobar). Results: (i) Pressure has suffered a monotonic increase from 1988 to 2016, that is compared to a seasonal volatile transport model, from which tight constraints on a combination of albedo and emissivity of N2 ice are derived; (ii) A central flash observed on 2015 June 29 is consistent with New Horizons REX profiles, provided that (a) large diurnal temperature variations (not expected by current models) occur over Sputnik Planitia and/or (b) hazes with tangential optical depth of about 0.3 are present at 4-7 km altitude levels and/or (c) the nominal REX density values are overestimated by an implausibly large factor of about 20% and/or (d) higher terrains block part of the flash in the Charon facing hemisphere., Comment: 21 pages, 11 figures

We describe the design and performance of the near-infrared (1.51--1.70 micron), fiber-fed, multi-object (300 fibers), high resolution (R = lambda/delta lambda ~ 22,500) spectrograph built for the Apache Point Observatory Galactic Evolution Experiment (APOGEE). APOGEE is a survey of ~ 10^5 red giant stars that systematically sampled all Milky Way populations (bulge, disk, and halo) to study the Galaxy's chemical and kinematical history. It was part of the Sloan Digital Sky Survey III (SDSS-III) from 2011 -- 2014 using the 2.5 m Sloan Foundation Telescope at Apache Point Observatory, New Mexico. The APOGEE-2 survey is now using the spectrograph as part of SDSS-IV, as well as a second spectrograph, a close copy of the first, operating at the 2.5 m du Pont Telescope at Las Campanas Observatory in Chile. Although several fiber-fed, multi-object, high resolution spectrographs have been built for visual wavelength spectroscopy, the APOGEE spectrograph is one of the first such instruments built for observations in the near-infrared. The instrument's successful development was enabled by several key innovations, including a "gang connector" to allow simultaneous connections of 300 fibers; hermetically sealed feedthroughs to allow fibers to pass through the cryostat wall continuously; the first cryogenically deployed mosaic volume phase holographic grating; and a large refractive camera that includes mono-crystalline silicon and fused silica elements with diameters as large as ~ 400 mm. This paper contains a comprehensive description of all aspects of the instrument including the fiber system, optics and opto-mechanics, detector arrays, mechanics and cryogenics, instrument control, calibration system, optical performance and stability, lessons learned, and design changes for the second instrument., Comment: 81 pages, 67 figures, PASP, accepted