## Systematic error cancellation for a four-port interferometric polarimeter

A. Kogut a,* , D. J. Fixsen a,b a Goddard Space Flight Center, Code 665, 8800 Greenbelt Road, Greenbelt, MD, USA 20771

b University of Maryland, College Park, MD, USA 20742

Abstract. The Primordial Inflation Explorer (PIXIE) is an Explorer-class mission concept to measure the gravitationalwave signature of primordial inflation through its distinctive imprint on the linear polarization of the cosmic microwave background (CMB). Its optical system couples a polarizing Fourier transform spectrometer to the sky to measure the differential signal between orthogonal linear polarization states from two co-pointed beams on the sky. The double differential nature of the four-port measurement mitigates beam-related systematic errors common to the two-port systems used in most CMB measurements. Systematic errors coupling unpolarized temperature gradients to a false polarized signal cancel to first order for any individual detector. This common-mode cancellation is performed optically, prior to detection, and does not depend on the instrument calibration. Systematic errors coupling temperature to polarization cancel to second order when comparing signals from independent detectors. We describe the polarized beam patterns for PIXIE and assess the systematic error for measurements of CMB polarization.

Keywords: cosmic microwave background, systematic error, beam patterns, polarimeter, Fourier transform spectrometer.

- * Address all correspondence to: Alan Kogut, alan.j.kogut@nasa.gov

## 1 Introduction

Polarization of the cosmic microwave background (CMB) provides a powerful test of the physics of the early universe. An arbitrary pattern of linear polarization mapped over the sky may be decomposed into a spatially symmetric component (even parity E-modes) and an anti-symmetric component (odd parity B-modes). Scalar sources such as temperature or density perturbations can only generate even-parity E-modes, while gravitational waves created during an inflationary epoch in the early universe can generate either parity. Detection of the B-mode signal in the CMB polarization field is thus recognized as a 'smoking gun' signature of inflation, testing physics at energies inaccessible through any other means. 1-8

The amplitude of the gravitational wave signal depends on the energy scale of inflation as

$$E = 1. 0 6 \times 1 0 ^ { 1 6 } \left ( \frac { r } { 0. 0 1 } \right ) ^ { 1 / 4 } \, \text{GeV}$$

where r is the power ratio of gravitational waves to density fluctuations. 9 In most large-field models, r is predicted to be of order 0.01, corresponding to polarized amplitude 30 nK or energy near the Grand Unified Theory scale, 10 16 GeV. Signals at this amplitude could be detected by a dedicated polarimeter, providing a critical test of a central component of modern cosmology. Detection of a gravitational-wave component in the CMB polarization would have profound implications for both cosmology and high-energy physics. It would provide strong evidence for inflation, provide a direct, model-independent determination of the relevant energy scale, and test physics at energies a trillion times beyond those accessible to particle accelerators. Generation of gravitational waves during inflation is purely a quantum-mechanical process: a detection of the B-mode signal provides direct observational evidence that gravity obeys quantum mechanics.

Characterizing the CMB to measure polarization at the parts-per-billion level requires careful control of systematic errors. A particular concern are systematic errors related to the instrument optics, which can couple the much brighter unpolarized temperature fluctuations into a false polarization signal. All CMB instruments must couple the detectors to the sky, and must therefore account for potential beam-related systematic errors. An extensive literature discusses common effects and mitigation strategies. 10-13

The Primordial Inflation Explorer (PIXIE) is an Explorer-class mission designed to measure the inflationary signature in polarization as well as distortions from the 2.725 K blackbody spectrum induced by energy-releasing processes at more recent cosmological epochs. 14 Its projected sensitivity of a few nK on degree angular scales or larger corresponds to a limit r &lt; 10 -3 at 5 standard deviations. PIXIE differs from most CMB polarimeters in its use of a polarizing Fourier transform spectrometer coupled to the sky through a multi-moded optical system. The double differential nature of the resulting four-port measurement minimizes beam-related systematic errors common to the two-port systems used in most CMB measurements. We describe the polarized beam patterns for PIXIE and assess the systematic error for measurements of CMB polarization.

## 2 PIXIE Optical System

Acommonimplementation for CMB polarimetry images the sky onto a set of polarization-sensitive detectors. Since each detector is sensitive to a single linear polarization from the sky (although two or more detectors may share a physical pixel), the resulting system may be described as a two-port device with any polarization comparison between detectors occurring post-detection. In contrast, the PIXIE optical system forms a four-port device (Fig 1). Reflective optics couple a polarizing Fourier Transform Spectrometer (FTS) to the sky. The FTS introduces an optical phase delay between the two input beams, and routes recombined beams to non-imaging concentrators at each of two output ports. Within each concentrator, a pair of polarization-sensitive detectors measure the power as a function of optical phase delay. Let glyph[vector] E = E x x ˆ+ E y y ˆ represent the electric field incident

Fig 1 Atwo-port system (left) couples a single linear polarization from the sky to a single detector. PIXIE's polarizing Fourier Transform Spectrometer operates as a four-port device (right) with two input ports open to the sky and two output ports terminated by polarization-sensitive detectors. Interfering the two beams cancels the effects of common mode beam ellipticity, as each detector then couples to both linear polarizations from the sky.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00558v1-figure-bb97c63d-5856-451a-9a10-75bfb10994d3.jpg)

from the sky. The power P at the detectors as a function of the phase delay z may be written

$$P _ { L x } \ & = \ 1 / 2 \ \int \{ \ ( E _ { A x } ^ { 2 } + E _ { B y } ^ { 2 } ) + ( E _ { A x } ^ { 2 } - E _ { B y } ^ { 2 } ) \cos ( 4 z \omega / c ) \ \} d \omega \\ P _ { L y } \ & = \ 1 / 2 \ \int \{ \ ( E _ { A y } ^ { 2 } + E _ { B x } ^ { 2 } ) + ( E _ { A y } ^ { 2 } - E _ { B x } ^ { 2 } ) \cos ( 4 z \omega / c ) \ \} d \omega \\ P _ { R x } \ & = \ 1 / 2 \ \int \{ \ ( E _ { A x } ^ { 2 } + E _ { B y } ^ { 2 } ) + ( E _ { B y } ^ { 2 } - E _ { A x } ^ { 2 } ) \cos ( 4 z \omega / c ) \ \} d \omega$$

where ˆ x and ˆ y refer to orthogonal linear polarizations, L and R refer to the detectors in the left and right concentrators, A and B refer to the two input beams, ω is the angular frequency of incident radiation, and the factor of 4 reflects the symmetric folding of the optical path. When both input ports are open to the sky, the power at each detector consists of a dc term proportional to the intensity E 2 x + E 2 y (Stokes I ) plus a term modulated by the phase delay z , proportional to the linear polarization E 2 x -E 2 y (Stokes Q ) in instrument-fixed coordinates. Each detector is thus sensitive to the difference between orthogonal linear polarizations from the two input ports, with the difference now occurring pre-detection. Rotation of the instrument about the beam axis rotates the instrument coordinate system relative to the sky to allow separation of Stokes Q and U parameters on the sky.

Rotation of the instrument relative to the sky can produce systematic errors in the recovered polarization if the instrument beams are not azimuthally symmetric. This effect has been well studied for two-port devices which couple the sky directly to a single polarization-sensitive detector. The dominant systematic error for such a device is temperature-polarization coupling as the beam ellipticity interacts with local gradients in the unpolarized sky intensity, producing a spin-dependent signal degenerate with true polarization. Temperature-polarization coupling can be mitigated in hardware using such techniques as polarization modulation (rapidly switching a single detetor between orthogonal polarization states) or in analysis using a well-measured beam profile.

The beam response of a four-port system is considerably more complicated than the two-port systems commonly used for CMB measurements. Although a four-port system may still employ mitigation strategies such as rotation or beam mapping, the double differential nature of the four-port measurement provides additional mitigation against common spin-dependent polarization errors while simultaneously providing a means to identify and correct residual effects. Figure 2 shows the PIXIE optical path. Consider (in a time-reversed sense) the path through the optics taken by photons leaving the ˆ x detector in the left-side concentrator. Since this detector is sensitive to a single linear polarization, the photons exiting the left-side concentrator are entirely in the ˆ x polarization. A series of polarizing wire grids within the FTS splits the beam and rotates the polarization so that half the initial power exits through port A in the ˆ y polarization while the other half exits through port B in the ˆ x polarization (see, e.g., Appendix A of reference 14). A set of reflective mirrors then couples ports A and B to the sky while preserving the polarization state. Stops at transfer mirror 5 and at the entrance to the concentrator circularize the beam. The entire instrument, including all baffling, remains isothermal at 2.725 K so that rays scattered out of the beam terminate on an isothermal black surface. Such rays contribute to the photon noise budget but do not introduce artifacts in the beams.

Let us define the beam pattern of the concentrator as H x ( θ, φ ) for the ˆ x polarization and H θ,φ y ( ) for the ˆ y polarization, where the angular coordinates θ and φ are referred to the sky. Similarly, we define the beam pattern for the fore-optics (defined as all elements in the optical

chain skyward of the concentrator feed) as F x ( θ, φ ) and F y ( θ, φ ) . Using subscripts L and R to distinguish the two concentrator ports and A and B for the two fore-optic ports, we may re-write Equation 2 as

$$P _ { L x } & \ \infty \ \ H _ { L x } \left [ F _ { A x } E _ { x } ^ { 2 } - F _ { B y } E _ { y } ^ { 2 } \right ] \\ P _ { L y } & \ \infty \ \ H _ { L y } \left [ F _ { A y } E _ { y } ^ { 2 } - F _ { B x } E _ { x } ^ { 2 } \right ] \\ P _ { R x } & \ \infty \ \ H _ { R x } \left [ F _ { B x } E _ { x } ^ { 2 } - F _ { A y } E _ { y } ^ { 2 } \right ] \\ P _ { R y } & \ \infty \ \ H _ { R y } \left [ F _ { B y } E _ { y } ^ { 2 } - F _ { A x } E _ { x } ^ { 2 } \right ],$$

where for clarity we suppress the dependence on angular coordinates ( θ, φ ) as well as the phase delay integral over frequency. Two points are apparent. First, the signal at any single detector depends on the convolution of the concentrator beam profile with the differential beam profile generated by the A- and B-side fore-optics. To the extent that the A- and B-side optics have identical beam patterns, the detectors produce no response from an unpolarized sky, regardless of the intensity gradient on the sky or the ellipticity of the fore-optics. This common mode cancellation is performed optically, prior to detection, and does not depend on the instrument calibration. Second, the beam pattern for the concentrator horn appears only as a common-mode multiplicative factor. Systematic errors coupling temperature anisotropy to polarization thus cancel to second order when comparing signals from independent detectors.

Figure 3 illustrates the multiple levels of common-mode subtraction. An ideal azimuthally symmetric beam would introduce no temperature-polarization coupling. Real beams, however, will have some ellipticity (left column). Rotation of an elliptical beam couples to unpolarized gra-

Fig 2 PIXIE optical signal path. The left panel shows the optical elements within the Fourier Transform Spectrometer while the right panel shows the physical layout.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00558v1-figure-5c3df949-8032-4f6d-9581-22a12d6e4fc2.jpg)

Fig 3 Cartoon illustrating signal cancellation from differential beam profiles. Colored regions indicate the beam shape, while the white lines indicate the polarization state accepted by each beam.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00558v1-figure-2869799f-541b-4e3f-9b6b-537a7f6f6ff4.jpg)

dients in the sky to produce a time-dependent signal degenerate with a true polarization signal. If, however, two beams sensitive to opposite polarization states but with the same ellipticity are compared, the common-mode ellipticity cancels for unpolarized emission, leaving no net temperaturepolarization coupling (second column). Only the differential ellipticity produces a net temperaturepolarization coupling, which appears at second order in the beam difference. The PIXIE optics employ such beam cancellation in the A -B comparison for a single detector (third column, with the differential ellipticity greatly exaggerated). Comparisons between different detectors provide an additional level of cancellation. Each concentrator contains two detectors sensitive to orthogonal polarization states (Eq. 2), which view the same sky through the same fore-optics. If the differential ellipticity between the A and B sky beams is the same for the ˆ x polarization as for the ˆ y polarization, the net temperature-polarization coupling for the single-detector output will cancel to second order in the detector-pair difference (right-most column). Alternatively, an orthogonal linear combination of detector pairs can be chosen to cancel the sky signal, thereby isolating any beam effects. Such measurements can be used both as confirmation of the expected amplitude of the beam differences and to correct residual beam effects in the sky data.

The double-differential beam cancellation of PIXIE's four-port optical system reduces the sensitivity to unpolarized gradients on the sky. The following sections use Monte Carlo ray-trace code to evaluate the common-mode and differential beam patterns. We quantify the expected systematic error response for ideal optics and show the minimal degradation in performance after accounting for machining and assembly tolerances.

## 3 Single-Detector Response

Systematic errors in the PIXIE four-port optical system depend on successive differences in the beam patterns. We may write the individual fore-optics beam patterns in terms of the linear combinations

$$F \ & = \ ( F _ { A x } + F _ { A y } + F _ { B x } + F _ { B y } ) / 4 \\ \Delta \ & = \ ( F _ { A x } + F _ { A y } - F _ { B x } - F _ { B y } ) / 4 \\ \delta \ & = \ ( F _ { A x } - F _ { A y } + F _ { B x } - F _ { B y } ) / 4 \\ \epsilon \ & = \ ( F _ { A x } - F _ { A y } - F _ { B x } + F _ { B y } ) / 4$$

Fig 4 Linear combinations of the PIXIE fore-optics showing the common-mode and differential beam patterns. The spatial ( ∆ ) and polarization ( δ ) asymmetries are small compared to the mean beam pattern F . Contours for the common-mode response F are shown at amplitude 0.3, 0.7, and 0.9 to highlight the circular tophat beam structure. Note the change in scale for the three differential beam patterns.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00558v1-figure-13482d5b-359b-4104-bb70-dae92fab8a8f.jpg)

to distinguish the common-mode beam pattern F = F θ, φ ( ) from the differential beam patterns ∆ (A-B spatial asymmetry), δ ( ˆ x -ˆ y polarization asymmetry), and glyph[epsilon1] (spatial/polarization cross term). With these definitions, the individual beam patterns become

$$F _ { A x } \ & = \ \ F + \delta + \Delta + \epsilon \\ F _ { A y } \ & = \ \ F - \delta + \Delta - \epsilon \\ F _ { B x } \ & = \ \ F + \delta - \Delta - \epsilon \\ F _ { B y } \ & = \ \ F - \delta - \Delta + \epsilon \.$$

Note that these four linear combinations represent a complete set, carrying all information for 2 ports in 2 linear polarizations.

Figure 4 shows the common-mode and differential beam patterns, using a Monte Carlo raytrace code to propagate 10 11 rays through the PIXIE fore-optics. As expected, the beams are dominated by the common-mode illumination F . Since F by definition is the average of the beams, it can not generate any differential ellipticity and thus can not generate temperature-polarization coupling regardless of its azimuthal structure.

Differences between the left and right beams are are measured by the A -B spatial asymme-

try ∆( θ, φ ) . Out-of-plane reflections at the secondary mirror and folding flat generate a dipolar modulation in ∆ with rms amplitude 0.015 of the common-mode beam pattern. This is the largest differential mismatch between the beams. The instrument is symmetric about the left-right midplane so that, by design, the A and B beams are mirror images of each other (Figure 5). Structure within one of the beams will thus be reflected left-to-right in the other beam, maximizing the net effect along the left-right direction.

Differences between polarization states are measured by the polarization asymmetry δ θ, φ ( ) . Since both the ˆ x and ˆ y polarizations from the detectors are launched at 45 ◦ relative to the symmetry plane of the instrument ( § 4), this term is small (of order 10 -4 of the common-mode pattern). For completeness, there is also a spatial/polarization cross term glyph[epsilon1] ( θ, φ ) . This term is also small (of order 10 -4 ). As shown below, it does not couple to temperature-polarization mixing, but appears as a small perturbation on the amplitude of the measured polarization signal.

Using these definitions, it is straightforward (if somewhat tedious) to show that

$$P _ { L x } \ & = \ H _ { L x } [ \quad Q \, F \, + \, Q \epsilon \, + \, I ( \ \delta + \Delta ) \ ] \\ P _ { L y } \ & = \ H _ { L y } [ \, - \, Q \, F \, + \, Q \epsilon \, + \, I ( - \delta + \Delta ) \ ] \\ P _ { R x } \ & = \ H _ { R x } [ \quad Q \, F \, - \, Q \epsilon \, + \, I ( \ \delta - \Delta ) \ ] \\ P _ { R y } \ & = \ H _ { R y } [ \, - \, Q \, F \, - \, Q \epsilon \, + \, I ( - \delta - \Delta ) \ ] \.$$

The first term in brackets represents the desired polarized sky signal Q ( θ, φ ) , convolved with the mean fore-optics beam pattern. The second term, Q glyph[epsilon1] , convolves the true sky polarization with the cross beam pattern glyph[epsilon1] ( θ, φ ) . This term is small. The cross beam pattern may be written as the double difference

$$\epsilon = ( F _ { A x } - F _ { A y } ) - ( F _ { B x } - F _ { B y } )$$

and is thus second order in the beam difference. Furthermore, since this term does not mix the Stokes parameter Q with either U or I , it only appears as a scale error in the amplitude of the

Fig 5 Schematic of the PIXIE optical system showing the symmetric polarization response at the beam apertures. The Fourier transform spectrometer interferes a single linear polarization from one side of the instrument with the orthogonal polarization from the other side. By construction, the ˆ x polarization on the A side is simply the mirror reflection of the ˆ y polarization on the B side.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00558v1-figure-81f70651-ef37-404e-b925-6705ad7de6e3.jpg)

true sky polarization and may be absorbed by the calibration. The final term represents systematic temperature-polarization coupling.

The left-right symmetry of the PIXIE optics minimizes temperature-polarization coupling. PIXIE's optical design interferes the ˆ x polarization from one beam with the ˆ y polarization from the other beam (Eq. 3). The optical system is symmetric about the central plane, so that the ˆ x polarization from one beam is the mirror reflection of the ˆ y polarization from the other beam (Fig 5). This enforces a reflection symmetry such that

$$F _ { A x } ( \theta, \phi ) \ & = \ F _ { B y } ( \theta, - \phi ) \\ F _ { A y } ( \theta, \phi ) \ & = \ F _ { B x } ( \theta, - \phi )$$

where the azimuthal angle φ is defined from the midline. Note that this left-right symmetry is not equivalent to an ˆ x y -ˆ symmetry since the ˆ x y -ˆ coordinate system is rotated by 45 ◦ with respect to the optical midline. Temperature-polarization mixing thus depends on the linear combinations

$$\delta + \Delta \ & = \ F _ { A x } - F _ { B y } \\ \delta - \Delta \ & = \ F _ { B x } - F _ { A y }$$

proportional to the anti-symmetric component of the difference between the beams.

The spacecraft spin combines with the mirror symmetry of the instrument optics to further minimize temperature-polarization coupling. Each detector is sensitive to a single linear polarization (Stokes Q in a coordinate system fixed with respect to the instrument). The entire spacecraft rotates about the instrument boresight to interchange the roles of ˆ x and ˆ y polarization at the detectors, allowing full characterization of the Stokes Q and U parameters on the sky. True sky polarization is modulated at twice the spacecraft spin frequency,

$$Q _ { \text{inst} } = Q _ { \text{sky} } \cos ( 2 \gamma ) + U _ { \text{sky} } \sin ( 2 \gamma )$$

where γ is the spin angle of the instrument with respect to the sky. Temperature-polarization mixing is dominated by the anti-symmetric component of the differential beam pattern from the instrument fore-optics. Anti-symmetric signals can only appear at odd harmonics of the spacecraft spin, and may readily be distinguished from true sky polarization.

We quantify the suppression of temperature-polarization systematic errors using the spindependent moments of the differential beam patterns. The instantaneous power at each detector depends on the convolution of the beam pattern (in instrument-fixed coordinates) with the sky signal (rotated from sky to instrument coordinates). Azimuthal asymmetry in the beam patterns causes the measured power to vary with the spacecraft spin angle. We thus compute the coefficients

$$a _ { m } \ & = \ \int B ( \Omega ) \cos ( m \phi ) d \Omega \\ b _ { m } \ & = \ \int B ( \Omega ) \sin ( m \phi ) d \Omega$$

where B represents one of the linear combinations of beam patterns (Eq. 4) and d Ω = sin( θ dθdφ ) is computed in instrument coordinates centered on the boresight.

Fig 6 Decomposition of the PIXIE differential beam patterns by spin angle. The common-mode beam F is sensitive only to polarized emission and does not contribute to temperature-polarization systematic errors. The mirror symmetry of the PIXIE optics suppresses temperature-polarization mixing from the A-B spatial asymmetry (beam ∆ ) by a factor of 10 -6 (see text).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00558v1-figure-be0f012b-62db-42b8-803b-93b582303e61.jpg)

Figure 6 shows the power P m = a 2 m + b 2 m as a function of spin moment m . The odd-even asymmetry in spin moment m is superposed atop an overall decrease in power with increasing m . The noise floor at P ≈ 10 -12 reflects shot noise from the discrete ray-trace simulation. Recall that the common-mode beam pattern F is sensitive only to polarized emission on the sky (Eq. 6) and does not create temperature-polarization errors even for the m = 2 case. Systematic errors from temperature-polarization coupling are dominated by the m = 2 mode of the A-B spatial asymmetry ∆ , and are suppressed by a factor 10 -6 relative to the polarization response in the common-mode beam.

## 4 Additional Symmetries

The mirror symmetry of PIXIE's differential 4-port interferometer suppresses systematic errors from temperature-polarization coupling by 6 orders of magnitude for the single-detector response. Additional symmetries between different detectors allow further suppression and identification of beam-related systematic errors. The left and right concentrators are identical, resulting in left-right symmetry

$$H _ { L x } ( \theta, \phi ) \ & = \ H _ { R x } ( \theta, - \phi ) \\ H _ { L y } ( \theta, \phi ) \ & = \ H _ { R y } ( \theta, - \phi )$$

for identical polarization states. This is similar to the left-right symmetry in Eq. 8 except that the symmetry is now between identical polarization states on opposite sides of the instrument.

Differences between the two polarizations ˆ x and ˆ y within a single concentrator can occur, corresponding to the difference between the E-plane and H-plane beam patterns for a single-moded feed. PIXIE's multi-moded operation reduces this effect, which vanishes in the geometric optics limit. We further reduce the effect by rotating the concentrator so that the symmetry axes of the

square aperture lie at ± 45 ◦ relative to the ˆ x and ˆ y polarization vectors (Fig 7). The resulting beams in ˆ x and ˆ y are equivalent linear combinations of the E-plane and H-plane beam patterns, so that

$$H _ { L x } \ & \approx \ H _ { L y } \\ H _ { R x } \ & \approx \ H _ { R y }$$

with residuals resulting from small displacements in the rotation angle. 15 Without loss of generality, we may follow Eq. 5 to decompose the beam pattern from each horn into a component common to all four detectors plus a set of differential beam patterns,

$$H \ & = \ ( H _ { L x } + H _ { L y } + H _ { R x } + H _ { R y } ) / 4 \\ \rho \ & = \ ( H _ { L x } - H _ { L y } + H _ { R x } - H _ { R y } ) / 4 \\ \tau \ & = \ ( H _ { L x } + H _ { L y } - H _ { R x } - H _ { R y } ) / 4 \\ \kappa \ & = \ ( H _ { L x } - H _ { L y } - H _ { R x } + H _ { R y } ) / 4$$

so that

$$H _ { L x } & \ = \ H + \rho + \kappa + \tau \\ H _ { L y } & \ = \ H - \rho + \kappa - \tau \\ H _ { R x } & \ = \ H + \rho - \kappa - \tau \\ H _ { R y } & \ = \ H - \rho - \kappa + \tau$$

where the horn parameters are defined analogously to the fore-optics in Eq. 4. Figure 8 shows the common-mode and differential beam patterns from the concentrator horn. Asymmetries from the off-axis orientation appear at the few-percent level. As with the fore-optics, the differential beams are dominated by an anti-symmetric (dipolar) component.

Appendix A shows the full response for each detector. Retaining only terms to first order in the

Fig 7 Schematic showing the orientation of the PIXIE concentrator. The square aperture is rotated 45 ◦ to minimize any differences between the ˆ x and ˆ y polarization.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00558v1-figure-70271e51-1366-4c0d-9fd7-682a3821b913.jpg)

Fig 8 Common-mode and differential beam patterns for the PIXIE feed horn concentrators. The feed horn beam pattern does not directly source T → B systematic errors, but only modulates the effect from the differential foreoptics. The off-axis design creates dipolar modulation in the differential beam patterns ρ and τ , while the square shape is reflected in the quadrupolar modulation for κ . Contours for the common-mode response H are shown at amplitude 0.3, 0.7, and 0.9.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00558v1-figure-0540bdc5-73ba-4160-bfd7-5afe141d34b0.jpg)

beam differences, the signals at each detector become

$$\text{rences, the signals at each} \text{ detector become} \\ P _ { L x } & = \ Q H F & + \ Q \left ( H \epsilon + F \rho + F \kappa + F \tau \right ) \\ & & + \ I H ( \delta + \Delta ) \\ P _ { L y } & = - Q H F & + \ Q \left ( H \epsilon + F \rho - F \kappa + F \tau \right ) \\ & & - \ I H ( \delta - \Delta ) \\ P _ { R x } & = \ Q H F & + \ Q \left ( - H \epsilon + F \rho - F \kappa + F \tau \right ) \\ & & + \ I H ( \delta - \Delta ) \\ P _ { R y } & = - Q H F & + \ Q \left ( - H \epsilon + F \rho + F \kappa - F \tau \right ) \\ & & - \ I H ( \delta + \Delta ) \. \\ \colon \text{errors from the concentrator beam pattern appear in the second (bolarization amplitude)}$$

Systematic errors from the concentrator beam pattern appear in the second (polarization amplitude) term. Although larger in amplitude than the spatial-polarization error Q glyph[epsilon1] , these terms do not couple temperature to polarization and so may be absorbed into the calibration. The final term representing temperature-polarization mixing is dominated at lowest order by the differential error from the fore-optics.

## 5 Combined Detector Response

PIXIE's four detectors share different portions of the optical system (left or right concentrator, ˆ x or ˆ y polarization). Linear combinations of the post-detection signals can either eliminate or isolate specific systematic error signals, providing additional safeguards against temperature-polarization mixing. For example, we may combine all 4 detectors to yield the sum signal

$$\left [ P _ { L x } - P _ { L y } + P _ { R x } - P _ { R y } \right ] / 4 \, & = \, Q \, H F \\ & + \, Q \, \epsilon \tau \\ & + \, I \, H \delta \\ & + \, I \, \Delta \tau$$

where we now retain terms to second order in the differential beam patterns. As before, the first term is the true sky polarization, convolved with the combined common-mode beam pattern from the feed horn and fore-optics. The second term affects only the amplitude of the true sky polarization and may be absorbed into the calibration. The final two terms represent systematic errors coupling temperature anisotropy to polarization.

We use Monte Carlo ray-trace simulations to quantify the expected amplitude of these terms. Table 1 summarizes the common-mode and differential beam patterns for the PIXIE optical system. The differential beam patterns are small compared to the common-mode response. We compare the weighted beam area of the differential beams to the weighted area of the common-mode beam pattern,

$$f = \frac { \int | \Delta ( \theta, \phi ) | d \Omega } { \int | F ( \theta, \phi ) | d \Omega } \,,$$

computed similarly for each of the 6 differential beam patterns. The differential beams have fractional area of a few percent for the concentrator, and 10 -2 to 10 -5 for the more symmetric foreoptics. The differential beams are dominated by a dipolar modulation ( m = 1 ) which does not lead to temperature-polarization mixing. The systematic error response to spin modulation at m = 2 is typically of order 10 -6 or smaller.

We may now quantify the systematic error terms in the post-detection linear combination. The third term I Hδ in Eq. 17 is similar to the temperature-polarization mixing I H ∆ from a single detector (Eq. 9), but reduced in amplitude by a factor of 200 due to replacing the A-B differential beam pattern ∆ with the smaller ˆ x -ˆ y differential beam pattern δ . The lower response to m = 2 modulation from the δ differential beam (compared to the ∆ beam) produces additional systematic error suppression. The final term I ∆ τ also represents temperature-polarization mixing, but now appears at second order in small beam differences and is reduced by a factor 20 in amplitude from the single-detector error. The m = 2 spin modulation of the τ differential beam yields additional suppression.

We may also choose linear combinations of detectors to cancel the polarized sky signal Q HF , thereby isolating specific systematic error signals. Such measurements of the systematic error signals can be used both to correct the sky measurements and as confirmation of the expected effect from beam pattern differences. For example, the orthogonal combination of four detectors

## becomes

$$[ P _ { L x } - P _ { L y } - P _ { R x } + P _ { R y } ] \, / 4 \, & = \, Q \, F \kappa \\ & + \, Q \, \epsilon \tau \\ & + \, I \, \Delta \rho \\ & + \, I \, \delta \kappa \.$$

We may again use Table 1 to estimate the amplitude of each term. Unpolarized CMB signals I have amplitude of order 100 µ K, while the E-mode polarization Q is of order 1 µ K. Multiplying each CMB term by the relative beam area of each beam pattern yields an estimate of the relative amplitude of each term (prior to spin modulation). The difference signal is dominated by the term I ∆ ρ , representing the convolution of the unpolarized CMB anisotropy with the double beam difference ∆ ρ . We may instead choose to compare signals from the two detectors sharing a common concentrator. A similar analysis shows that the detector-pair combination ( P Lx + P Ly ) / 2 is dominated by the term I H ∆ which isolates a single differential beam for measurement and correction. Similar linear combinations can isolate other terms.

## 6 Tolerance

Mirror symmetries within PIXIE's differential optics suppress systematic errors coupling unpolarized structure in the sky to a false polarized signal. Positioning errors in the optical components during assembly can distort the beams from the ideal beam patterns. We quantity the resulting degradation in optical performance using 30 Monte Carlo realizations of the PIXIE optical system. For each realization, we adjust the position of each optical element allowing both translation and rigid-body rotation about its nominal orientation assuming assembly and machining tolerances of ± 0.05 mm drawn from a random Gaussian distribution. After adjusting all optical elements, we follow the paths of 10 9 rays through the adjusted optical system to define the distorted beam patterns.

The PIXIE optical system is robust to typical machining and assembly tolerances. The FTS

Table 1 Spin Modulation of the Common-Mode and Differential Beam Patterns

|                      | Fore-Optics   | Fore-Optics   | Fore-Optics   | Fore-Optics     |
|----------------------|---------------|---------------|---------------|-----------------|
| Parameter            | F             | ∆             | δ             | glyph[epsilon1] |
| Peak Amplitude       | 1             | 3 × 10 - 2    | 2 × 10 - 4    | 1 × 10 - 4      |
| Relative Beam Area f | 1             | 2 × 10 - 2    | 8 × 10 - 5    | 4 × 10 - 5      |
| Power (m=1)          | 7 × 10 - 5    | 2 × 10 - 5    | 9 × 10 - 11   | 2 × 10 - 10     |
| Power (m=2)          | 3 × 10 - 4    | 1 × 10 - 6    | 4 × 10 - 9    | 7 × 10 - 11     |

|                      | Horn Concentrator   | Horn Concentrator   | Horn Concentrator   | Horn Concentrator   |
|----------------------|---------------------|---------------------|---------------------|---------------------|
| Parameter            | H                   | ρ                   | τ                   | κ                   |
| Peak Amplitude       | 1                   | 6 × 10 - 2          | 8 × 10 - 2          | 1 × 10 - 2          |
| Relative Beam Area f | 1                   | 3 × 10 - 2          | 5 × 10 - 2          | 7 × 10 - 3          |
| Power (m=1)          | 5 × 10 - 8          | 2 × 10 - 4          | 9 × 10 - 4          | 4 × 10 - 8          |
| Power (m=2)          | 6 × 10 - 4          | 8 × 10 - 10         | 4 × 10 - 6          | 3 × 10 - 5          |

Fig 9 Differences between the nominal beam patterns from Figure 4 and the distorted beam patterns after allowing for machining and assembly tolerances. Beam patterns are shown from a single Monte Carlo realization in which the position and orientation of each optical element are perturbed about the nominal configuration.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00558v1-figure-43a01d54-1b9b-48c8-a9ce-6d012b80ce2d.jpg)

left and right transfer mirror sets and the mid-plane septum containing the polarizing grids are each machined from a single block of aluminum. The relative position and orientation of the mirrors or grids within each set have the ± 0.02 mm tolerance of computer-aided milling machines. This minimizes relative displacement of these components during assembly (although each set can still be displaced as a rigid body). All optical elements as well as the supporting structure are fabricated from the same material (aluminum) so that self-similar thermal contraction retains optical alignment. Alignment of the primary, secondary, and folding flat mirrors relative to the FTS assumes somewhat looser ± 0.05 mm tolerance typical of pinned construction. Since these components are machined individually, tolerancing errors should be uncorrelated. Figure 9 shows the difference between the nominal beam patterns and the distorted patterns for a single Monte Carlo realization of the distorted optical system. The dominant effect is an angular displacement of order 3 ′ between the A and B beams on the sky, caused by displacements of the primary and secondary mirrors. This in turn creates an anti-symmetric (dipolar) pattern in both the mean beam ( F ) and the A-B spatial asymmetry ( ∆ ). Angular displacement of the beam centroid couples to spin moment m = 1 and does not induce additional temperature-polarization mixing.

Figure 10 compares the spin dependence of the nominal beam patterns to the distorted beams generated from a single Monte Carlo realization of the full optical system. It is similar to the ideal beam patterns shown in Figure 6, except that the position and pointing of each optical element

Fig 10 Effects of machining and assembly tolerances on the differential beam pattern ∆( θ, φ ). We compare the differential beam pattern for the nominal optical configuration to a Monte Carlo realization with all optical elements perturbed from their nominal positions. The distorted patterns are shown as a function of spin moment m for the detectors sensitive to ˆ x sky polarization in both the left-side and right-side concentrators. The distorted optical system still shows suppression of order 10 -5 for temperature-polarization coupling at m = 2 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00558v1-figure-4198c041-fadf-4f30-8298-f10366c759d1.jpg)

has been perturbed by an amount randomly chosen from a Gaussian normal distribution of width 0.05 mm. We now also include the illumination of the (perturbed) fore-optics by the (perturbed) concentrator. For clarity, we compare the spin decomposition for the nominal and distorted configurations for a single choice of differential beam. Temperature-polarization coupling for a single detector is dominated by the A-B differential beam ∆( θ, φ ) (Eq. 16). PIXIE has 4 detectors; we show the distorted beam decomposition for detectors observing the same ( ˆ x ) sky polarization from either the left or right concentrator. Compared to the ideal system, the distorted optical system has a larger response to systematic error coupling at m = 2 , but the response is still suppressed by a factor of 10 5 compared to the true sky polarization.

## 7 Discussion

Systematic errors coupling unpolarized anisotropy to a false polarized signal are a common concern to CMB polarimeters. PIXIE's optical design provides several layers of mitigation compared to instruments imaging the CMB across a large (kilo-pixel) focal plane. Missions employing kilopixel arrays across large fields of view must account for the systematic degradation in beam shape from coma and shear for detectors farther from the center of the focal plane. All four PIXIE detectors, in contrast, lie at the center of the focal plane, allowing eam deformation to be minimized (Fig 4).

More importantly, PIXIE's four-port optical system provides three distinct levels of differential measurement. The Fourier transform spectrometer produces a signal that depends on the difference between two nearly-identical beams on the sky. This differential measurement is performed optically, prior to detection, and is independent of detector calibration. We use ray-trace simulations to evaluate the differential beam patterns after removing the common-mode response. The differential beams can be described in terms of the spatial asymmetry between the A- and B-sides

of the instrument, the polarization asymmetry between the ˆ x and ˆ y response, plus a cross term for the mixed spatial-polarization difference. All of the differential beams are small compared to the common-mode response. The largest effect is the spatial (A-B) asymmetry, which has only 1.5% of the common-mode response. The other differential beams have response below 0.01%.

PIXIE's symmetric design further reduces systematic error response from the differential beam cancellation. The FTS interferes the ˆ x polarization from the A-side beam with the ˆ y polarization from the B-side beam (Fig 5). The optical system is symmetric about the mid-plane between the two sides, which forces the ˆ x polarization from one beam to be the mirror reflection of the ˆ x polarization from the other beam. The A-B mirror reflection combines with the A-B beam subtraction to produce an anti-symmetric (dipole) response in the differential beam patterns. The anti-symmetric part of the differential beam pattern does not contribute to the systematic error from temperature-polarization coupling. Each detector samples a single polarization state; the entire instrument spins about the boresight to allow full sampling of the sky polarization. True polarized signals appear at twice the spin frequency, while anti-symmetric signals can only appear at odd harmonics of the spin. Systematic errors from temperature-polarization coupling thus depend only on the m = 2 component of the differential beam patterns, which are dominated by a dipole ( m = 1 ) response. Ray-trace models of the PIXIE beams show that the response at m = 2 is reduced by an additional factor of 10 6 or more. In principle, the optical system could further be optimized to suppress the m = 2 differential beam response, moving power to other m values that do not participate in temperature-polarization mixing. This has not yet been done but is planned for future development.

Finally, we may follow the common practice for CMB measurements and combine the postdetection signals from individual detectors. The four detectors are mounted in identical concentrators and view the same sky direction through the same optical path. Combining all four detectors cancels the leading effects from differential beams in the single-detector signal, reducing the systematic error response by a factor of 1000 or more compared to the individual detectors. Conversely, orthogonal linear combinations of 2 or 4 detectors can cancel the polarized sky signal to isolate, identify, and model specific systematic effects from the individual differential beam patterns.

Systematic error suppression in the differential PIXIE optics is robust against typical machining and assembly tolerances. We combine ray-trace optical simulations with Monte Carlo realizations of distorted PIXIE optics to evaluate both the individual beam patterns and the resulting systematic error response. Each Monte Carlo realization of then PIXIE optics perturbs each optical element (mirrors, folding flats, polarizing grids, etc) in both position and orientation by an amount drawn from a Gaussian distribution whose width is set by typical machining/assembly tolerances of 0.05 mm. The dominant effect of such tolerance errors is an angular displacement of the A-side beams relative to the B-side beams. The two beams are normally co-pointed on the sky; after accounting for tolerances the beams are typically mis-aligned by 3 ′ . This is small compared to the 2.6 ◦ width of the common-mode beams; the resulting dipolar beam asymmetries predominantly effect the m = 1 spin moment and do not couple efficiently to polarization. The distorted optical system still provides suppression of the m = 2 temperature-polarization systematic error by factor of order 10 5 .

## Appendix A: Full Single-Detector Systematic Error

Expanding Eqs. 6 and 15 yields individual detector signals

$$P _ { L x } = \ Q \, H F & \quad \ \ + \ Q \, H \epsilon \, + \, I \, H \delta \, + \, I \, H \Delta \\ & \quad + \ Q \, F \rho \, + \, Q \, \rho \epsilon \, + \, I \, \rho \delta \, + \, I \, \rho \Delta \\ & \quad + \ Q \, F \kappa \, + \, Q \, \kappa \epsilon \, + \, I \, \kappa \delta \, + \, I \, \kappa \Delta \\ & \quad + \ Q \, F \tau \, + \, Q \, \tau \epsilon \, + \, I \, \tau \delta \, + \, I \, \tau \Delta$$

$$P _ { L y } = - Q \, H F & \quad + \, Q \, H \epsilon \, - \, I \, H \delta \, + \, I \, H \Delta \\ & + \, Q \, F \rho \, - \, Q \, \rho \epsilon \, + \, I \, \rho \delta \, - \, I \, \rho \Delta \\ & - \, Q \, F \kappa \, + \, Q \, \kappa \epsilon \, - \, I \, \kappa \delta \, + \, I \, \kappa \Delta \\ & + \, Q \, F \tau \, - \, Q \, \tau \epsilon \, + \, I \, \tau \delta \, - \, I \, \tau \Delta$$

$$P _ { R x } = \ Q \, H F & \quad \ \ - \ Q \, H \epsilon \, + \, I \, H \delta \, - \, I \, H \Delta \\ & \quad + \, Q \, F \rho \, - \, Q \, \rho \epsilon \, + \, I \, \rho \delta \, - \, I \, \rho \Delta \\ & \quad - \, Q \, F \kappa \, + \, Q \, \kappa \epsilon \, - \, I \, \kappa \delta \, + \, I \, \kappa \Delta \\ & \quad - \, Q \, F \tau \, + \, Q \, \tau \epsilon \, - \, I \, \tau \delta \, + \, I \, \tau \Delta$$

$$P _ { R y } = - Q \, H F & \quad \ \ \ - \, Q \, H \epsilon \, - \, I \, H \delta \, - \, I \, H \Delta \\ & + \, Q \, F \rho \, + \, Q \, \rho \epsilon \, + \, I \, \rho \delta \, + \, I \, \rho \Delta \\ & + \, Q \, F \kappa \, + \, Q \, \kappa \epsilon \, + \, I \, \kappa \delta \, + \, I \, \kappa \Delta \\ & - \, Q \, F \tau \, - \, Q \, \tau \epsilon \, - \, I \, \tau \delta \, - \, I \, \tau \Delta$$

where now we retain all terms to second order.

## References

- 1 V. A. Rubakov, M. V. Sazhin, and A. V. Veryaskin, 'Graviton creation in the inflationary universe and the grand unification scale,' Physics Letters B 115 , 189-192 (1982).
- 2 R. Fabbri and M. D. Pollock, 'The effect of primordially produced gravitons upon the anisotropy of the cosmological microwave background radiation,' Physics Letters B 125 , 445-448 (1983).
- 3 L. F. Abbott and M. B. Wise, 'Constraints on generalized inflationary cosmologies,' Nuclear Physics B 244 , 541-548 (1984).
- 4 A. G. Polnarev, 'Polarization and anisotropy induced in the microwave background by cosmological gravitational waves,' Astronomicheskii Zhurnal 62 , 1041-1052 (1985).
- 5 R. L. Davis, H. M. Hodges, G. F. Smoot, et al. , 'Cosmic microwave background probes models of inflation,' Physical Review Letters 69 , 1856-1859 (1992).
- 6 L. P. Grishchuk, 'Cosmological perturbations of quantum-mechanical origin and anisotropy of the microwave background,' Physical Review Letters 70 , 2371-2374 (1993).

- 7 M. Kamionkowski, A. Kosowsky, and A. Stebbins, 'Statistics of cosmic microwave background polarization,' Physical Review D 55 , 7368-7388 (1997).
- 8 U. Seljak and M. Zaldarriaga, 'Signature of Gravity Waves in the Polarization of the Microwave Background,' Physical Review Letters 78 , 2054-2057 (1997).
- 9 D. H. D. H. Lyth and A. A. Riotto, 'Particle physics models of inflation and the cosmological density perturbation,' Physics Reports 314 , 1-146 (1999).
- 10 W. Hu, M. M. Hedman, and M. Zaldarriaga, 'Benchmark parameters for CMB polarization experiments,' Physical Review D 67 , 043004 (2003).
- 11 D. O'Dea, A. Challinor, and B. R. Johnson, 'Systematic errors in cosmic microwave background polarization measurements,' Monthly Notices of the Royal Astronomical Society 376 , 1767-1783 (2007).
- 12 C. Rosset, V. B. Yurchenko, J. Delabrouille, et al. , 'Beam mismatch effects in cosmic microwave background polarization measurements,' Astronomy and Astrophysics 464 , 405-415 (2007).
- 13 M. Shimon, B. Keating, N. Ponthieu, et al. , 'CMB polarization systematics due to beam asymmetry: Impact on inflationary science,' Physical Review D 77 , 083003 (2008).
- 14 A. Kogut, D. J. Fixsen, D. T. Chuss, et al. , 'The Primordial Inflation Explorer (PIXIE): a nulling polarimeter for cosmic microwave background observations,' Journal of Cosmology and Astroparticle Physics 7 , 25 (2011).
- 15 A. J. Kogut and D. J. Fixsen, 'Polarized beam patterns from a multimoded feed for observations of the cosmic microwave background,' Journal of Astronomical Telescopes, Instruments, and Systems 4 , 014006 (2018).

Alan Kogut is an astrophysicist at NASA's Goddard Space Flight Center. He received his A.B. from Princeton University and his PhD from the University of California at Berkeley. His research focuses on observations of the frequency spectrum and linear polarization of the cosmic microwave background and diffuse astrophysical foregrounds at millimeter and sub-mm wavelengths.

Dale Fixsen is an astrophysicist at the University of Maryland. He received his B.S. from Pacific Lutheran University, Tacoma Washington and his PhD from Princeton University. His research focuses on the spectrum, temperature, and polarization of the cosmic microwave background and the radio and infrared cosmic backgrounds.
A. Kogut D. J. Fixsen
Systematic error cancellation for a four-port interferometric polarimeter
The Primordial Inflation Explorer (PIXIE) is an Explorer-class mission
concept to measure the gravitational-wave signature of primordial inflation
through its distinctive imprint on the linear polarization of the cosmic
microwave background (CMB). Its optical system couples a polarizing Fourier
transform spectrometer to the sky to measure the differential signal between
orthogonal linear polarization states from two co-pointed beams on the sky. The
double differential nature of the four-port measurement mitigates beam-related
systematic errors common to the two-port systems used in most CMB measurements.
Systematic errors coupling unpolarized temperature gradients to a false
polarized signal cancel to first order for any individual detector. This
common-mode cancellation is performed optically, prior to detection, and does
not depend on the instrument calibration. Systematic errors coupling
temperature to polarization cancel to second order when comparing signals from
independent detectors. We describe the polarized beam patterns for PIXIE and
assess the systematic error for measurements of CMB polarization.