A New Limit on CMB Circular Polarization from SPIDER

TL;DR

This paper reports a new upper limit on the circular polarization of the cosmic microwave background (CMB) using data from SPIDER's 2015 balloon flight, significantly improving previous constraints and demonstrating the potential for future, more sensitive measurements.

Contribution

The study introduces a novel method leveraging the non-zero circular-to-linear polarization coupling of HWP modulators to constrain CMB circular polarization over a wide angular scale range.

Findings

SPIDER sets the first constraints on CMB circular polarization at 95 and 150 GHz.

The new limits improve previous bounds by several orders of magnitude.

Constraints on the Stokes V power spectrum range from 141 to 255 μK² at 150 GHz.

Abstract

We present a new upper limit on CMB circular polarization from the 2015 flight of SPIDER, a balloon-borne telescope designed to search for $B$-mode linear polarization from cosmic inflation. Although the level of circular polarization in the CMB is predicted to be very small, experimental limits provide a valuable test of the underlying models. By exploiting the non-zero circular-to-linear polarization coupling of the HWP polarization modulators, data from SPIDER's 2015 Antarctic flight provide a constraint on Stokes $V$ at 95 and 150 GHz from $33<\ell<307$. No other limits exist over this full range of angular scales, and SPIDER improves upon the previous limit by several orders of magnitude, providing 95% C.L. constraints on $\ell (\ell+1)C_{\ell}^{VV}/(2\pi)$ ranging from 141 $\mu K ^2$ to 255 $\mu K ^2$ at 150 GHz for a thermal CMB spectrum. As linear CMB polarization experiments become increasingly sensitive, the techniques described in this paper can be applied to obtain even stronger constraints on circular polarization.

Tables

Table 1. Table 1: Thicknesses of the 95 GHz HWP components. Spider ’s HWPs are made from single-crystal birefringent sapphires, which are AR coated to maximize transmission (Bryan, 2014 ) . At 95 GHz, the sapphires are AR coated with quartz wafers that are glued at the centers and held by spring clips at the edges. The measured thicknesses of these materials for each HWP are listed in the table. Since the adhesive covers only a small fraction of the total surface area, it is ignored in the calculation of s 𝑠 s . However, the possibility of a narrow gap between the sapphire and the quartz is taken into account.

95 GHz Receivers	X2	X4	X6
Sapphire (mm)	4.97 $\pm$ 0.01	4.94 $\pm$ 0.01	4.97 $\pm$ 0.01
Top Quartz Layer (mm)	0.420 $\pm$ 0.015	0.429 $\pm$ 0.015	0.427 $\pm$ 0.015
Bottom Quartz Layer (mm)	0.419 $\pm$ 0.015	0.419 $\pm$ 0.015	0.422 $\pm$ 0.015
Gap (mm)	0.01 $\pm$ 0.01	0.01 $\pm$ 0.01	0.01 $\pm$ 0.01

Table 2. Table 2: Thicknesses of the 150 GHz HWP components. The 150 GHz HWPs are AR coated with Cirlex, a polyimide film, which is adhered with a melted HDPE bond layer. The measured thicknesses of these materials for each HWP are listed in the table. The uniformity in the thickness of the Cirlex sheets can likely be attributed to a common production batch.

150 GHz Receivers	X1	X3	X5
Sapphire (mm)	3.21 $\pm$ 0.01	3.26 $\pm$ 0.01	3.23 $\pm$ 0.01
Top Cirlex Layer (mm)	0.250 $\pm$ 0.005	0.250 $\pm$ 0.005	0.250 $\pm$ 0.005
Bottom Cirlex Layer (mm)	0.250 $\pm$ 0.005	0.250 $\pm$ 0.005	0.250 $\pm$ 0.005
HDPE Bond Layer (mm)	0.006 $\pm$ 0.001	0.006 $\pm$ 0.001	0.006 $\pm$ 0.001

Table 3. Table 3: Refractive indices of the HWP materials. The listed values assume a temperature of approximately 4 K and observing bands in the range of 50-200 GHz.

Material	Refractive Index ($n$)	Reference
Sapphire (fast axis)	3.019 $\pm$ 0.003	Bryan et al. (2010a)
Sapphire (slow axis)	3.336 $\pm$ 0.003	Bryan et al. (2010a)
Quartz (fused)	1.95 $\pm$ 0.01	Bryan (2014)
Cirlex	1.94 $\pm$ 0.01	Bryan (2014)
HDPE	1.56 $\pm$ 0.01	Lamb (1996)

Table 4. Table 4: Spider ’s polarization efficiency γ 𝛾 \gamma . These values were obtained by combining calculations of the four HWP non-ideality parameters with measurements of the detector cross-polarization response. Since γ 𝛾 \gamma is dominated by the contribution from the HWPs, the same value is used for every detector on a given receiver.

Receiver Name	Frequency	Polarization Efficiency ($\gamma$)
X1	150 GHz	0.959 $\pm$ 0.005
X2	95 GHz	0.965 $\pm$ 0.001
X3	150 GHz	0.950 $\pm$ 0.008
X4	95 GHz	0.964 $\pm$ 0.001
X5	150 GHz	0.956 $\pm$ 0.005
X6	95 GHz	0.964 $\pm$ 0.003

Table 5. Table 5: Spider ’s 95% C.L. limits on ℓ ​ ( ℓ + 1 ) ​ C ℓ V ​ V / ( 2 ​ π ) ℓ ℓ 1 superscript subscript 𝐶 ℓ 𝑉 𝑉 2 𝜋 \ell(\ell+1)C_{\ell}^{VV}/(2\pi) for a CMB source based on Figure 4 .

Bin Center ($\ell$)	95 GHz Limit ($\mu K^2$)	150 GHz Limit ($\mu K^2$)
45	1088	195
70	783	153
95	842	149
120	853	141
145	856	142
170	985	164
195	1032	177
220	1129	197
245	1254	242
270	1455	244
295	1760	255

Equations

$$M_{HWP}=\left[\begin{array}[]{cccc}T&\rho&0&0\\ \rho&T&0&0\\ 0&0&c&-s\\ 0&0&s&c\end{array}\right].$$

$$d=I{M}_{II}+Q{M}_{IQ}+U{M}_{IU}+V{M}_{IV}.$$

$$M_{IV}=s\gamma \sin (2{\theta }_{HWP}-2{\xi }_{det}).$$

Full text

1 Introduction

Anisotropies in the intensity and linear polarization of the Cosmic Microwave Background (CMB) have provided a wealth of information about the history and contents of the universe. Standard cosmological models do not predict a measurable amount of circular polarization, characterized by the Stokes $V$ parameter, in the CMB; as such, any detection of a primordial $V$ signal would be of enormous interest. A variety of secondary physical processes may produce circular polarization in the CMB at very low levels. For instance, Faraday conversion can transform existing linear polarization into circular polarization in both the magnetic fields of galaxy clusters (Cooray et al., 2003) and the relativistic plasma remnants of Population III stars (De & Tashiro, 2015). Magnetic fields in the primordial universe (Giovannini, 2009; Zarei et al., 2010), scattering from the cosmic neutrino background (Mohammadi, 2014), and photon-photon interactions in neutral hydrogen (Sawyer, 2015) have also been shown to potentially produce CMB circular polarization. Additional sources include postulated extensions to QED such as Lorentz-invariance violating operators (Colladay & Kosteleckỳ, 1998; Alexander et al., 2009), axion-like pseudoscalar particles (Finelli & Galaverni, 2009), and non-linear photon interactions (through effective Euler-Heisenberg Lagrangians) (Motie & Xue, 2012). A brief review of some of these generation mechanisms can be found in King & Lubin (2016). Despite the wide range of physical processes they invoke, all of these mechanisms predict levels of circular polarization that are unlikely to be accessible with current technology.

Nevertheless, circular polarization measurements provide a valuable test of the standard cosmological model and the physics behind these generation mechanisms. Yet there are relatively few published limits. MIPOL reported the strongest constraint on large angular scales ($\ell<30$) in 2013, providing 95% C.L. measurements ranging from $\Delta V/T_{\mathrm{CMB}}\leq 2.4\times 10^{-4}$ to $\Delta V/T_{\mathrm{CMB}}\leq 4.3\times 10^{-4}$ at 33 GHz (Mainini et al., 2013). This is roughly an order of magnitude better than the previous 95% C.L. limit of $\Delta V/T_{\mathrm{CMB}}\leq 4\times 10^{-3}$ at 33 GHz at $\ell\approx 10$ (Lubin et al., 1983). On smaller angular scales, the only reported measurement comes from the VLA, which set 95% C.L. limits at 5 GHz between $\Delta V/T_{\mathrm{CMB}}\leq 2.2\times 10^{-4}$ and $\Delta V/T_{\mathrm{CMB}}\leq 0.6\times 10^{-4}$ for a range of angular scales with $\ell>3000$ (Partridge et al., 1988).

These limits are more than 7 orders of magnitude higher than the best measurements of the linear polarization power spectra, but there are no contemporary experiments with the primary goal of improving them. However, some modern linear polarization experiments, such as Spider, can take advantage of this vast disparity to set stronger limits as a consequence of their polarization modulation techniques.

Spider is a balloon-borne CMB telescope that is searching for a $B$-mode linear polarization signal from cosmic inflation (Fraisse et al., 2013; Rahlin et al., 2014). During its first flight in January 2015, Spider made maps of approximately 10% of the sky with degree-scale angular resolution in 95 and 150 GHz observing bands. The analysis of the linear polarization data from this flight is currently in progress. In this paper, we exploit non-idealities of Spider’s half-wave plate (HWP) polarization modulators to obtain a new upper limit on CMB circular polarization.

The Spider payload features six monochromatic receivers housed in a shared cryostat (Gudmundsson et al., 2015). Each receiver includes a stepped HWP polarization modulator to reduce the potential impact of systematic errors due to beam asymmetries and instrument polarization (Bryan, 2014; Bryan et al., 2016). Although Spider’s antenna-coupled TES bolometers are not sensitive to variations in circular polarization (Ade et al., 2015), non-idealities of the HWPs allow a measurement of the Stokes $V$ parameter after combining maps made at several HWP angles. The calculation of Spider’s circular polarization coupling is described in the next section. Section 3 details how this coupling is used to derive a circular polarization limit. The significance of this result and prospects for future measurements are discussed in Section 4.

2 Coupling to Circular Polarization

A birefringent material forms a half-wave plate when the difference in the optical path length between waves polarized along the fast and slow crystal axes is exactly half of the photon wavelength. An ideal HWP rotates the polarization plane of the light passing through it by $2\theta_{HWP}$, where $\theta_{HWP}$ is the angle between the incoming polarization plane and the slow crystal axis. However, this condition can only be exactly satisfied at a single frequency. Similarly, the single-layer anti-reflection (AR) coatings applied to Spider’s HWPs are not uniformly efficient over the observing bands. When combined, these conditions lead to a frequency-dependent reduction in transmission through the HWPs and induce non-ideal polarization modulation effects as the HWPs are rotated. Since Spider’s observing bands have roughly 20% bandwidths, the magnitude of such effects could be significant.

Following Bryan et al. (2010b), a non-ideal HWP can be modeled with four parameters that can be broadly interpreted as the total transmission $T$, the difference in transmission between the fast and slow axes $\rho$, the linear polarization response $c$, and the coupling to circular polarization $s$. In terms of these parameters, the Mueller matrix of a HWP with its birefringent crystal axes oriented along the horizontal and vertical directions can be written as

[TABLE]

For an ideal HWP, $T=1$, $c=-1$, and $\rho=s=0$. The ideal case captures the effect of the HWP on linear polarization signals and has no coupling between linear and circular polarization ($s=0$). In real HWPs, however, these parameters can deviate significantly from their ideal values. While simulations have shown that these non-idealities are not problematic for detecting a $B$-mode signal at Spider’s anticipated sensitivity level (O’Dea et al., 2011), they allow Spider to measure circular polarization to the extent that $s$ is non-zero.

The sky signal in a detector timestream $d$ is given in terms of the Stokes parameters $I$, $Q$, $U$, and $V$ and the instrument Mueller matrix elements $M_{XY}$ by

[TABLE]

The instrument Mueller matrix is calculated in Bryan et al. (2010b) by multiplying the Mueller matrices of every element in the optical chain, including $M_{HWP}$ from Equation 1. For the purposes of this paper we are interested only in the instrument Mueller matrix element $M_{IV}$. Generalizing the treatment in Bryan et al. (2010b) for arbitrary detector angles, it is straightforward to show that the $V$ parameter couples to a detector timestream as

[TABLE]

Here $\theta_{HWP}$ is the HWP angle and $\xi_{det}$ is the detector angle, both of which are defined relative to the instrument. Note that $M_{IV}$ does not depend on the rotational orientation of the instrument relative to the sky and can be positive, negative, or zero depending on the relative HWP and detector angles. The overall polarization efficiency of the instrument is described by $\gamma$, while $s$ describes the coupling to circular polarization from the HWP non-idealities. Note that $s$ does not appear in the $M_{II}$, $M_{IQ}$, or $M_{IU}$ matrix elements in Equation 2 and therefore is not used in Spider’s linear polarization analysis.

Spider’s six receivers are assigned names consisting of the letter ‘X’ followed by the numbers 1 through 6, where the even numbers refer to 95 GHz receivers and the odd numbers to 150 GHz receivers. Each receiver has a dedicated HWP and therefore a unique value of the $s$ non-ideality parameter. It can be calculated as described in Bryan et al. (2010b) from the thicknesses and refractive indices of the HWP materials, the spectrum of the observed source, and the shape of the observing band. Similar HWP modeling techniques have been used for the linear polarization properties of sapphire HWPs by Savini et al. (2006) and found to be in good agreement with experimental measurements (Pisano et al., 2006; Savini et al., 2009; Bryan, 2014).

For the results presented in this paper, uncertainties in the component properties lead to significant uncertainty in $s$ for each HWP, which is quantified with Monte Carlo simulations. We use the temperature derivative of the CMB blackbody spectrum for the source in the baseline case, as well as the thicknesses and uncertainties of the HWP components listed in Tables 1 and 2, and the refractive indices and uncertainties of the materials in Table 3. Since the refractive indices of sapphire should be the same for every HWP, we use the same randomly drawn values of the two indices for all receivers in each iteration of the $s$ calculations.

To take Spider’s observing bands into account, we use Fourier Transform Spectrometer (FTS) measurements made just prior to flight (Gudmundsson, 2014). However, correctly interpreting these measurements relies on knowing the frequency dependence of the coupling to the Rayleigh-Jeans calibration source. Although the intensity of the source has a $\nu^{2}$ frequency dependence, the beam throughput ($A\Omega$) scales as $\nu^{-2}$ in the beam-filling limit. For Spider, the source is not entirely beam filling, and calculations indicate that this coupling should be approximately $\nu^{-1.5}$. This leads to larger absolute values of $s$ than in the beam-filling case. However, due to the relatively large uncertainty in the calculation of this scaling, we adopt a conservative approach in this paper and assume a $\nu^{-2}$ coupling, likely underestimating $|s|$.

The probability distributions of the $s$ parameters for a CMB source for each Spider receiver are shown in Figure 1. These are derived from 10,000 Monte Carlo simulations, which calculate a new value of $s$ for each iteration following the methodology presented by Bryan et al. (2010b), using randomly drawn sets of physical parameters based on the central values and uncertainties listed in Tables 1-3. The distributions of $s$ for the 150 GHz systems exclude zero at roughly 2- to 4-$\sigma$. At 95 GHz, they include $s=0$ within the 1-$\sigma$ range. However, to the extent that these distributions are truly good estimators of the $s$ probability distributions, they can still be used to constrain the amplitude of circular polarization by virtue of the significant probability of non-zero $s$. Note that having three separate HWPs at each frequency greatly improves Spider’s statistical power to constrain $V$. The differences in distributions between receivers at the same frequency are caused by differences in the shapes of the measured observing bands for each receiver and in the measured thicknesses of the actual HWP components.

3 Results

During the 2015 flight, Spider observed approximately 4500 square degrees of sky near the southern Galactic pole, centered around roughly RA=50∘ and Dec=-35∘. For the first 7.5 days used in this analysis, almost the entire region was mapped every 12 sidereal hours following a sinusoidal azimuth scan profile and using a scan width of $\sim$70∘ peak-to-peak. Maps for the remaining 4.5 days covered smaller overlapping sub-regions using narrower sinusoidal azimuth scans with widths of $\sim$35∘ peak-to-peak (Shariff, 2015). The HWPs were held at fixed angles during each of these maps and rotated to new angles between them, following the patterns shown in Figure 2. Over the course of the flight, each receiver observed the sky at 8 discrete HWP angles nominally spaced at 22.5 degree intervals.

The data from individual receivers are combined into 4 independent sets, illustrated by the colored bands in Figure 2, which were optimized for separating the $Q$, $U$, and $V$ signals. Each of these sets is used to construct an independent $V$ map with a binned map-maker (Rahlin, 2016), using the values of polarization efficiency $\gamma$ listed in Table 4. If the $s$ value for each receiver was known exactly, Equations 2 and 3 could be used by the mapmaker to construct $V$ maps directly from the Spider data. Instead, since the values of $s$ actually follow broad probability distributions, and $s$ appears only in the $M_{IV}$ matrix element in Equation 2, we make $V$ maps assuming $s=1$ and later scale the resulting power spectra.

Before making these maps, glitches such as cosmic ray hits, payload transmitter signals, and thermal transients are identified and removed from the detector timestreams. This pipeline is shared with the linear polarization analysis and will be described more extensively in a future publication. Some detectors have been excluded from this analysis due to undesirable remaining timestream features, but a number of them may be recovered for future results. Here we use 681 detectors at 95 GHz and 1117 detectors at 150 GHz, rejecting an average of approximately 30% of the data from these timestreams. For this result we subtract a fifth-order polynomial fit in azimuth from each scan (approximately 30 seconds of data) to remove scan-synchronous noise.

Only part of Spider’s observing region is used for the circular polarization analysis, masking data outside of $30^{\circ}\leq\mathrm{RA}\leq 70^{\circ}$ and $-55^{\circ}\leq\mathrm{Dec}\leq-15^{\circ}$. The leakage from other signals to $V$ is subtracted in map-space from full timestream signal simulations based on Planck 100 and 143 GHz temperature-only input maps (Planck Collaboration, 2016a). This is dominated by $T$-to-$V$ leakage from the polynomial timestream filter, which is at the level of $\sim 10\;\mu K$ in the original $s=1$ maps and roughly 2 orders of magnitude smaller than our $V$ sensitivity. The $E$-to-$V$ leakage is about 5 orders of magnitude lower. The $V$ pipeline was verified through simulations in which an input signal-only $V$ map was observed following Spider’s scan strategy and then recovered after applying the same flagging and filtering to the re-observed timestreams.

The cross-spectra of the $s=1$ maps are estimated with PolSpice (Chon et al., 2004), which takes the sky mask into account. We apply a transfer function to account for the effects of timestream filtering and beam smoothing, where the beam correction is derived from map-domain fits to Planck temperature maps. The transfer function is obtained by comparing $TT$ spectra from smoothed Planck maps of the Spider region to spectra from simulated re-observations of those Planck maps that include Spider’s pointing and timestream filtering.

The $s=1$ cross-spectra for pairs of maps at a given frequency are then combined with Monte Carlo simulations. In each iteration, values of $s$ for receivers $i$ and $j$ are drawn from the distributions shown in Figure 1, and the cross-spectra are then scaled by $1/(s_{i}s_{j})$. Note that the selected $s$ values are slightly correlated due to the common sapphire indices. We calculate the weighted mean of the resulting values in each $\ell$ bin, weighting by the variance of the $s=1$ map cross-spectra in that bin. This process is repeated 10,000 times, and the mean and error in each $\ell$ bin are derived from the resulting distribution.

At 150 GHz we cross every pair of maps from each of the three receivers and four independent sets, excluding the auto-spectra, for a total of 66 cross-spectra. This includes crossing maps made simultaneously on different receivers because the noise has been shown to be no more correlated than in any of the other map pairs. At 95 GHz we only cross maps from the same receiver because the $s$ distributions allow both positive and negative $s$ values with significant probabilities. A sign error on either $s_{i}$ or $s_{j}$ (but not both) relative to the true value would flip the sign of the cross-spectrum, potentially suppressing real $V$ signals upon averaging. By restricting ourselves to the 18 cross-spectra that can be constructed from single-receiver maps, we ensure that $s_{i}s_{j}=s_{i}^{2}$ is always positive, thus avoiding this problem at the price of a small noise penalty.

Figure 3 shows Spider’s $VV$ CMB spectra at 95 and 150 GHz, neither of which indicate a significant detection of circular polarization. The mean values and errors are derived from the distributions of the $s$-scaled cross-spectra, and the spread in each of those distributions has contributions from both the distribution of the various cross-spectra and the distributions of $s$ values. Figure 4 shows the 95% C.L. limits on CMB circular polarization derived from these spectra, and the numerical values are provided in Table 5. Although the measurements are made at different frequencies, they are expressed in units of CMB temperature, which are the equivalent fluctuations of a 2.73 K blackbody required to produce the measured intensity variations. This result represents an improvement of several orders of magnitude over the previous best upper limit (Mainini et al., 2013) at a complementary range of angular scales.

However, Spider’s limits depend on the chosen source spectrum through the calculations of the HWP coupling parameters $s$. Many of the methods for generating CMB circular polarization described in Section 1 predict polarization signals with spectra of the form $\nu^{-1}$ or $\nu^{-3}$. We therefore recalculate Spider’s $s$ distributions for such source spectra and find that the $VV$ limits in Figure 4 typically become lower. Still expressed in CMB temperature units, they scale by factors of 0.39 and 0.10 respectively at 95 GHz and 1.02 and 0.30 respectively at 150 GHz. In all of these cases, Spider’s circular polarization limits are still many orders of magnitude above the predicted cosmological signals.

Similarly, the $VV$ limits presented in this paper can also be extended to upper limits on foreground circular polarization by recomputing $s$ with the appropriate source spectra. King & Lubin (2016) suggest that $\nu^{-3.4}$ is a reasonable model for synchrotron circular polarization. With this source spectrum, Spider’s circular polarization limits in Figure 4 scale by 0.08 at 95 GHz and 0.11 at 150 GHz, still using CMB temperature units. To obtain an estimate of Spider’s limit on the circular polarization of thermal dust, we use the linear polarization model of $\nu^{3.5}$ for the source spectrum (Planck Collaboration, 2016b) since we are not aware of any circularly polarized dust models. This leads to circular polarization limits that scale from Figure 4 by 0.27 at 95 GHz and 0.56 at 150 GHz. Although these models of the source spectra are relatively uncertain, the predicted $V$ foreground signals are many orders of magnitude below Spider’s sensitivity level.

4 Conclusion

This paper presents a new upper limit on CMB circular polarization from $33<\ell<307$ at 95 and 150 GHz. It was obtained by exploiting a non-ideality of the HWP polarization modulators used by Spider to measure linear polarization during a 2015 Antarctic flight. This represents an improvement of several orders of magnitude over the previous limit, providing 95% C.L. constraints on $\ell(\ell+1)C_{\ell}^{VV}/(2\pi)$ ranging from 141 $\mu K^{2}$ to 255 $\mu K^{2}$ at 150 GHz for a thermal CMB spectrum. When recalculated for $\nu^{-1}$ and $\nu^{-3}$ source spectra, this limit scales by 1.02 and 0.30 respectively. Data from Spider’s second flight, planned for December 2018, could provide increased sensitivity at 95 and 150 GHz as well as a new measurement at 280 GHz over the same range of angular scales.

As linear polarization experiments become increasingly sensitive, the techniques described in this paper can be applied to provide stronger constraints on CMB circular polarization. Several current and planned experiments use either HWPs or Variable-delay Polarization Modulators (VPMs) (Miller et al., 2016), both of which can be used to measure $V$. Although the current limit is many orders of magnitude larger than the most optimistic signal predictions, these measurements provide an observational test of the standard cosmological model and a wide range of physical processes. Since this limit is still about four orders of magnitude above modern linear polarization measurements, a dedicated experiment with better $V$-coupling could make significant improvements using existing technology.

Spider is supported in the U.S. by the National Aeronautics and Space Administration under grants NNX07AL64G, NNX12AE95G, and NNX17AC55G issued through the Science Mission Directorate and by the National Science Foundation through PLR-1043515. Logistical support for the Antarctic deployment and operations was provided by the NSF through the U.S. Antarctic Program. Support in Canada is provided by the Natural Sciences and Engineering Research Council and the Canadian Space Agency. Support in Norway is provided by the Research Council of Norway. Support in Sweden is provided by the Swedish Research Council through the Oskar Klein Centre (Contract No. 638-2013-8993). K.F. acknowledges support from DoE grant DE-SC0007859 at the University of Michigan. We also wish to acknowledge the generous support of the David and Lucile Packard Foundation, which has been crucial to the success of the project. The collaboration is grateful to the British Antarctic Survey, particularly Sam Burrell, for invaluable assistance with data and payload recovery after the 2015 flight. We thank Brendan Crill and Tom Montroy for significant contributions to Spider’s development. JMN wishes to thank Glenn Starkman for useful discussions about methods of generating CMB circular polarization. The computations described in this paper were performed on the GPC supercomputer at the SciNet HPC Consortium (Loken et al., 2010). SciNet is funded by the Canada Foundation for Innovation under the auspices of Compute Canada, the Government of Ontario, Ontario Research Fund - Research Excellence, and the University of Toronto.