Issue 
A&A
Volume 627, July 2019



Article Number  A160  
Number of page(s)  16  
Section  Astronomical instrumentation  
DOI  https://doi.org/10.1051/00046361/201834495  
Published online  22 July 2019 
Systematic effects induced by halfwave plate precession into measurements of the cosmic microwave background polarization
^{1}
Physics Department, Università di Roma “Sapienza”, Ple. Aldo Moro 2, 00185 Rome, Italy
^{2}
INFN – Sezione di Roma1, Ple. Aldo Moro 2, 00185 Rome, Italy
email: giuseppe.dalessandro@roma1.infn.it
^{3}
Dipartimento di Fisica e Scienze della Terra, Università degli Studi di Ferrara and INFN – Sezione di Ferrara, Via Saragat 1, 44100 Ferrara, Italy
^{4}
Institut d’Astrophysique Spatiale, CNRS, Univ. ParisSud, Université ParisSaclay, Bât. 121, 91405 Orsay Cedex, France
^{5}
Institut d’Astrophysique de Paris, CNRS, 98 Bis Boulevard Arago, 75014 Paris, France
^{6}
LERMA, Sorbonne Université, Observatoire de Paris, Université PSL, École normale supérieure, CNRS, Paris, France
Received:
23
October
2018
Accepted:
13
June
2019
Context. The primordial Bmode signal in the cosmic microwave background (CMB) represents the smoking gun of cosmic inflation, and measuring it is the main goal of current experimental effort. The most accessible method for measuring polarization features of the CMB radiation is a Stokes polarimeter based on the rotation of a halfwave plate (HWP).
Aims. Current observational cosmology is starting to be limited by systematic effects. A Stokes polarimeter with a rotating HWP has the advantage of mitigating a long list of potential systematics by modulating the linearly polarized component of the radiation, but the rotating HWP itself may introduce new systematic effects that must be under control. This represents one of the most critical parts in the design of a Bmode experiment. It is therefore mandatory to take all the systematic effects into account that the instrumentation can induce. We here present, simulate, and analyze the spurious signal arising from the precession of a rotating HWP.
Methods. We first derived an analytical formula to describe the systematic effect that is induced by the HWP precession on the propagating radiation, using the 3D generalization of the Müller formalism. We then performed several numerical simulations that show the effect induced on the Stokes parameters by this systematic. We also derived and discuss the effect on Bmodes as measured by a satellite experiment.
Results. We derive the analytical formula for the Stokes parameters from a Stokes polarimeter where the HWP follows a precessional motion with an angle θ_{0}. We show the result depending on the HWP inertia tensor, spinning speed, and on θ_{0}. The result of numerical simulations is reported as a simple timeline of the electric fields. Finally, assuming that the entire sky is observed with a satellite mission, we analyze the effect on Bmode measurements.
Conclusions. The effect is not negligible at the sensitivity of current Bmode experiments, therefore this systematic needs to be carefully considered for future experiments.
Key words: instrumentation: polarimeters / cosmic background radiation / cosmology: observations / techniques: polarimetric
© ESO 2019
1. Introduction
The BICEP2 experiment, designed to measure the cosmic microwave background (CMB) polarization, claimed the first detection of primordial Bmodes in 2014 (BICEP2 Collaboration 2014) by measuring the tensortoscalar ratio as . One year later, a joint effort involving the Planck and BICEP2 collaborations (BICEP2/Keck & Planck Collaborations 2015) revised this bound by publishing an upper limit of r < 0.12, which was obtained by removing the residual dust contamination from the BICEP2 maps. More recently, BICEP2 and Keck Array experiments (Keck Array & BICEP2 Collaborations 2016) further reduced this constraint to r < 0.07, which represents the strongest constraint to date on inflationary gravitational waves.
Further improving the constraint on the tensortoscalar ratio represents a hard challenge for current and future CMB experiments, which must accurately take into account gravitational lensing and foreground removal (Errard et al. 2016) in addition to having an excellent control of systematic effects (Wallis et al. 2017). Concerning the control of systematics, some experiments are designed with the capability of selfcalibrating and thus of removing some systematic effects (Piat et al. 2012). For other experiments that lack this feature, an accurate instrumental systematics prevision and laboratory calibrations are mandatory in order to separate systematic errors from scientific data (Natoli et al. 2018; Inoue et al. 2016; D’Alessandro et al. 2015; Johnson et al. 2017).
A standard device for polarization analysis is a Stokes polarimeter, composed of a halfwave plate (HWP) and a polarizer. The HWP (Pisano et al. 2008) induces a phase shift on the linear polarization by birefringence, and the polarizer selects one component, so that by rotating the HWP it is possible to modulate the linearly polarized fraction of the incoming radiation and extract the Stokes components (T, Q, U), given a reference frame. The HWP can spin fast (>1 Hz), see, for example, LSPE Collaboration (2012), Columbro et al. (2019), EBEX Collaboration (2018) and Thornton et al. (2016), in order to modulate the signal at high frequency, and it can also rotate step by step (<1 Hz) (Piat et al. 2012; Bryan et al. 2016) depending on the scanning strategy of the experiment. Systematic effects such as temperaturetopolarization leakage, chromatic HWP behavior, and scanning strategy, have been evaluated in EssingerHileman et al. (2016), Salatino et al. (2011, 2018), and Takakura et al. (2017) and were also measured by Kusaka et al. (2014).
We here analyze the systematic effect that is induced in observation of the CMB polarization by the precession of the HWP, in the specific case of a transmissive plate. We assume that the HWP free from other systematic effects.
In Sect. 2 we review the dynamic of a spinning cylindrical plate and define the precession rate. In Sect. 3 we provide an analytic study of the effect that is induced by the spinning and preceding HWP using the 3D Jones formalism for fully polarized radiation and the 3D generalization of the Müller formalism for partially polarized radiation. In Sect. 4 we use the equations derived before to show some results on electric field components produced with numerical simulations. Finally, we describe in Sect. 5 the effect that is induced on fullsky observation of the CMB, assuming a satellite mission with different scanning strategies.
2. Precessing halfwave plate theory
In this section we present equations to describe a precessing body that rotates along one symmetry axis. We introduce the main variables and show their evolution with time. We report here only the main equations that are essential for the results shown in the following sections. All the details of the computation are provided in Appendix A.
We approximate the HWP as a cylindrically symmetric rigid body, like a coin, and we define a reference system x̂–ŷ–ẑ, with the x̂–ŷ plane coincident with the base of the cylinder, see Fig. 1. We hypothesize that the HWP has a large spin angular momentum L_{s} = I_{s}ω_{s} along the symmetry axis, where I_{s} and ω_{s} are the moment of inertia and the angular velocity, respectively. In the unperturbed case, L_{s} coincides with the ẑ axis.
Fig. 1. Illustration of rigidbody precession where we show the angular momentum dislocations with respect to the x̂ and ŷ axes. 
We now introduce a small angular perturbation (θ_{0}) of the spin momentum L_{s}, expanding it in its components along the x̂ and ŷ axes, η and ξ. By introducing the precession rate γ, defined as the rotation velocity of L_{s} around the ẑ axis, η and ξ are simply defined as
These equations describe the torquefree precession of a spin axis that rotates in space with fixed angle θ_{0}.
The precession rate is given by γ = ω_{s}I_{s}/I_{⊥} (Kleppner & Kolenkow 1973), after I_{s} and I_{⊥} are defined as the inertia momenta parallel and perpendicular to the ẑ axis. Assuming that the cylindrical plate is a thin disk, we can approximate I_{s} = 2I_{⊥} and γ = 2ω_{s}, which means that the disk wobbles twice as fast as it spins.
The apparent rate of a thindisk precession for an observer on the rigid body is given by
In a torquefree precession we can identify two different rotations according to the reference frame we consider. In the fixed laboratory frame, the angular velocity vector rotates around the fixed ẑ axis (where the angular momentum vector lies), tracing the socalled space cone. In the reference frame integral with the rotating body, the angular momentum and angular velocity vector describe a circle around the symmetry axis of the cylinder, tracing the socalled body cone with a precession rate γ′.
3. 3D generalization of a Stokes polarimeter
A description of Stokes polarimetry by using Müller and Jones matrices is available in the literature, see, for instance, Bryan et al. (2010a), O’Dea et al. (2007), and Chuss et al. (2012). In this section and in Appendices B and C, we derive the analytic equation for a traditional Stokes polarimeter and for a Stokes polarimeter where the HWP has a precessional motion.
3.1. Jones formalism
We now use an extended Jones formalism (Jones 1941) to retrieve a formula for electromagnetic (EM) field intensity after the Stokes polarimeter. The Jones matrices used are defined in Appendix B followed by explicit calculations. The traditional Jones formalism can describe the state of fully polarized light with a twodimensional vector i = (E_{x}, E_{y}) and optical elements that change the state of input radiation with 2 × 2 matrices. We extend the Jones vector in three dimensions, i = (E_{x}, E_{y}, 0), where the ẑaxis is the propagation direction of the EM field, and the x̂–ŷ plane is perpendicular to ẑ, thus any optical element is described by 3 × 3 matrices (Sheppard 2011, 2014; OrtegaQuijano & ArceDiego 2013). Combining these matrices, we obtain a general formula for the Jones vector for the onaxis detector of a precessing polarimeter:
For the ideal case with η = ξ = 0°, when no precession occurs, we obtain the outgoing intensity from an ideal polarimeter:
For the general case with η, ξ ≠ 0° we find the intensity at the detector:
where we define the modulating functions g(ξ, θ) and f(η, ξ, θ):
We can write the outgoing intensity through the Stokes parameters as follows:
which is valid only if T^{2} = Q^{2} + U^{2}. The last equation gives the intensity for the onaxis detector of a Stokes polarimeter whose modulating element describes a torquefree precession. It is not merely a function of the HWP orientation about the ẑaxis, but also depends on the displacements from the x̂ and ŷaxes due to the precession.
3.2. Müller formalism
We now use a 3D extended Müller formalism (Sheppard et al. 2016; Samim et al. 2016) to retrieve a formula for an EM field intensity outgoing from the Stokes polarimeter. The Müller formalism (fully defined in Appendix C) is required for the case of CMB experiments because we wish to propagate partially polarized radiation through a polarimeter in order to extract the information about its polarization state. Starting from 3D extended Jones matrices (Appendix B), the 9 × 9 Müller matrix corresponding to each optical element can be easily obtained from Eq. (C.1):
where σ_{n} (n = [0, …, 8]) are the tracenormalized GellMann matrices (Eq. (C.2)). The Stokes polarimeter Müller matrix becomes
The output Stokes vector therefore is
Each component of the Stokes vector is defined in Eq. (C.4), and assuming the field entering the polarimeter has E_{z} = 0,
Based on Eq. (C.5), which describes the relation between the 3D Stokes components and the classic definition of Stokes parameters, we can determine the outgoing intensity:
where m_{ij} are the components of the Stokes polarimeter Müller matrix and are a function of the wobbling angles η, ξ. When η = ξ = 0, the Eq. (12) gives the wellknown formula of the Stokes polarimeter (Eq. (5)).
4. Phenomenology
Starting from Eq. (B.12), we built a simulation to show the effects of the precession for a specific case. We set the spinning frequency at 1 Hz and the ratio I_{⊥}/I_{s} = 0.502. The Jones vector used as input field in the simulation was i_{in} = (1, 0, 0), corresponding to a (1, 1, 0, 0) Stokes vector. Typical outputs of the simulations are reported in Fig. 2. We show the fractional residuals, defined as the difference between the intensity from a precessing HWP and from the ideal nonprecessing case, normalized to the maximum intensity of the ideal polarimeter.
Fig. 2. (a) Fractional residuals between the total intensities in a Stokes polarimeter with and without HWP precession. We fix the ratio I_{⊥}/I_{s} = 0.502 and consider different HWP precession angles θ_{0} = [1°,0.5°,0.25°], normalizing the difference with respect to the maximum intensity of an ideal Stokes polarimeter. (b) Fractional residuals, as before, but in a shorter time interval to highlight the beats induced by precession. (c) Power spectra with peaks at multiples and submultiples of the HWP spinning frequency with an amplitude related to θ_{0}. 
The time lines reported in Fig. 2a show the fractional residual for different amplitudes of the precession angle θ_{0}. These time lines clearly show several beats with an amplitude depending on the precession angle θ_{0} in a nonlinear way.
In Fig. 2c we report the power spectra of time lines for the ideal case (red) and for the precessing cases. The ideal case presents a single line at 4 Hz, as expected from a wave plate that spins at 1 Hz. The beats from a precessing HWP produce spurious peaks at multiples and submultiples of the spin frequency. The most noteworthy peak in the spectra is the one at lower frequencies because it is the result of the very slow modulation imposed by the precession, which in Fig. 2a has a period of about 60 s, implying a peak at 0.016 Hz in the spectrum. In the next section we illustrate how these beats depend on the geometrical parameters of the plate.
4.1. Spinning speed and ratio I_{⊥}/I_{s} effects
For a cylindrical HWP, including its support, with mass m, thickness h, and radius R, the components of the moment of inertia with respect to the principal axes are
where we assume a diagonal inertia tensor:
The frequency for the precessional motion is directly linked to the HWP regular spin frequency f_{s} and to the ratio I_{⊥}/I_{s} (Eq. (A.5)).
We can note that this ratio only depends on the cylinder height and radius as
We therefore explore different configurations as shown in Fig. 3, where we show the fractional residual with respect to the ideal case. We consider input radiation that is horizontally polarized, a precession angle θ_{0} = 1°, spin frequencies f_{s} = [0.1, 0.5, 1.0, 2.0] Hz and I_{⊥}/I_{s} = [0.501, 0.502, 0.506, 0.513]. As an example, these values for the I_{⊥}/I_{s} ratio correspond to an HWP with mass m = 1 kg, radius R = 16 cm, and thickness h ≃ [1.2, 1.7, 3.0, 4.5] cm; then I_{s} is fixed to 0.0128 kg m^{2}. In practice, the ratio of the components of the moment of inertia does not only depend on the plate thickness, but also on the structure of the HWP support.
Fig. 3. Fractional residual for different ratios I_{⊥}/I_{s} and different HWP spinning frequencies. Each column shows a different f_{s} value and each row a different I_{⊥}/I_{s} ratio. A thinner HWP, I_{⊥}/I_{s} → 0.500, has beats in the intensity over an extended period and a thicker HWP has shorter beats. By looking at the figure from left to right, it is clear that the spinning speed compresses the beats. The maximum amplitude remains the same because it only depends on θ_{0}, which is fixed to 1° in this particular simulation. 
Figure 3 shows the dependence of the effect from the HWP spinning speed and the inertia tensor: the simulation shows that a thinner HWP, I_{⊥}/I_{s} → 0.500, has beats in the intensity over an extended period and a thicker HWP has shorter beats. This is clear by looking at Fig. 3 from top to bottom. The value 0.5 is not possible because it corresponds to a null thickness.
The effect of different spinning speeds is to compress the beats. This is clear by looking at Fig. 3 from left to right. The maximum amplitude remains the same because it depends only on θ_{0}, which is fixed to 1° in this particular simulation.
In Fig. 4 we report the power spectra of the time lines for different ratios I_{⊥}/I_{s}. The spectra exhibit the effect discussed above, showing that the beat frequency moves to lower values as the ratio I_{⊥}/I_{s} is reduced toward the minimum value of 0.5 (for I_{⊥}/I_{s} = 0.513, the period of the beats is 10 s, corresponding to a peak in the spectrum at 0.1 Hz).
Fig. 4. Power spectra extracted from time lines showing the shifting of the spurious peak at low frequencies. The simulation considers input light that is horizontally polarized, a precession angle θ_{0} = 1°, and I_{⊥}/I_{s} = [0.501, 0.502, 0.506, 0.513]. 
4.2. Minor effect
As reported for example in Gu et al. (2018), a tilted HWP also changes its transmission properties because path of radiation inside the plate is different from the case of normal incidence. The resulting effect is a variation of retardance, which can be as high as 40% in the case discussed in Gu et al. (2018) for a source at 347 nm with a tilt of 5°. This effect is multiplicative with the ratio wavelength to thickness, which is much higher in the case of plates used in the optical bandwidth than in the case of the millimeter wavelength. In instrumentation devoted to observations of the CMB polarization, the plate thickness is on the same order of magnitude as the wavelength: 3.1 mm POLARBEAR2 (Hill et al. 2016), 3.05 mm SPIDER (Bryan et al. 2010b), 1.62 mm EBEX (EBEX Collaboration 2018), 3.2 mm QUBIC (Aumont et al. 2016), 3 mm LSPESWIPE (LSPE Collaboration 2012), and 3.05 mm ABS (Kusaka et al. 2014). The impact of the effect for millimeter astronomy can be examined through the variation between the input and output Stokes vector, and it is estimated to be orders of magnitude smaller than the effect due to electromagnetic field projection that is analyzed in the rest of this paper.
5. Fullsky simulations
In order to test the effect of the HWP precessional motion on CMB observations, we built an algorithm that can generate a realistic satellite scanning strategy in presence of spinning HWP and produces data timelines. We completed this software with a mapmaking algorithm that collapses data timelines into maps. All simulations are noisefree to better capture the impact of systematic effects.
5.1. Simulation pipeline
The scan simulator takes as inputs the details of a satellite scanning strategy, three spin rates, and two angles (see Das & Souradeep (2015) for a detailed description of the geometrical configuration):

Earth revolution velocity ω_{1},

precession velocity ω_{2},

satellite spin ω_{3},

precession angle α, that is, the angle between the satellite spin axis and the SunEarth direction,

boresight angle β, that is, the angle between the focal plane direction and the spin axis.
We only simulated a single detector placed at the center of the focal plane by illuminating a spinning HWP with f_{s} frequency. The systematic affecting the HWP was included in the data at the timeline level, and a simple rebinning mapmaking was used to average all the samples in T,Q,U Stokes parameters maps (Tegmark 1997). In this paper we consider Planck (Tauber et al. 2006), WMAP (Bennett et al. 2003), COrE (Natoli et al. 2018), and LiteBIRD (Sekimoto et al. 2018) scanning strategies. The input parameters we employed for these scanning strategies are listed in Table 1. As sampling rate we used 60 Hz.
Parameters for the scanning strategies adopted in the simulation pipeline.
5.2. Maps and results
5.2.1. Input map
The input sky map we used for fullsky simulations contained the solar dipole and Galactic diffuse foregrounds in temperature and a CMB realization, both in temperature and polarization. We decided to include foregrounds only in temperature in order to highlight the temperaturetopolarization leakage. The input 𝒞_{ℓ} used for the CMB realization is compatible with the best fit of the Planck 2015 release (Ade et al. 2016) with no tensor perturbations.
The foreground field was generated from the Commander solution delivered with the Planck 2015 release (Adam et al. 2016). It includes the primary temperature emissions (Planck Collaboration ES 2015): synchrotron, freefree, spinning dust, CO, and thermal dust emission, without considering their polarization contribution.
This map was modeled in order to highlight the temperaturetopolarization leakage induced by the HWP precession during the observations. We set the resolution parameter of the input map at HEALPix (Gorski et al. 2005) resolution N_{side} = 256, comfortable enough with respect to the Gaussian beam with FWHM = 60 arcmin. In order to evaluate the effect of the parameters chosen for the simulation, we ran a case with N_{side} = 128, finding the same results in terms of angular power spectra residual, except at smaller scales, where the pixel size matters independently of the presence of systematic effects.
5.2.2. Output maps
We performed several simulations with different configurations for the HWP. We varied the spin frequency, precession angle, and I_{s}/I_{⊥} ratio. For each simulation we compared input and output maps and computed the Bmode power spectrum. As a visual example, we show in Fig. 5 (top panel) the output maps for a mission for which we adopted a LiteBIRDlike scanning strategy and solved the Stokes parameters through an HWP with a spin frequency of 1 Hz, a precession angle of θ_{0} = 1°, and I_{⊥}/I_{s} = 0.514.
Fig. 5. Top panel: T, Q, and U maps reconstructed through a Stokes polarimeter where the HWP wobbles. Bottom panel: difference between input and output T, Q, and U maps showing the effect of the waveplate precession on observed maps. 
The residual maps (i.e., difference between output and input maps) in Fig. 5 (bottom panel) show the effect of the HWP wobble that induces variations of a few percent with respect to the input map. The effect is noticeable close to the Galactic plane and close to the maximum and minimum of the solar dipole, where the intensity emission is higher.
The effect on the maps is generated by the coupling between the satellite spin and the polarization modulation, which is affected by the precession, therefore we decided to test several conditions. In particular, when we slowed the HWP spin speed down to 0.1 Hz, the effect of the precession is emphasized. We show this in Fig. 6, where the systematic effect induced in the T, Q, and U maps, reported in histogramequalized color scale, is at the same level as in the input map.
Fig. 6. Temperature and polarization difference maps showing the effect of the waveplate precession for a slowspinning modulator. Modulation parameters: f_{s} = 0.1 Hz, θ_{0} = 1°, I_{⊥}/I_{s} = 0.514. 
5.2.3. Results
The angular power spectra from the output maps shown in Fig. 5 are reported up to l ∼ 200, given the limit imposed by the beam. The relative variations for both E and Bmodes (Fig. 7) show the effect of the precession, combined with the satellite spin, which dominates at small angular scales (l > 150).
Fig. 7. Input minus output relative difference of EE and BB power spectra computed from the maps shown in Fig. 5. 
When a 0.1 Hz spinning HWP is used, the synchronism with the satellite is slightly different and spurious Bmode polarization is detected at different angular scales (l > 100). These two cases show different matching between the systematic effect and the satellite spin (Figs. 8a and 8b).
Fig. 8. Panels a and c: BB and EE power spectra extracted from the input and output maps of Fig. 6; in panel c the two lines overlap. Panels b and d: normalized difference between input and output power spectra. 
5.3. Scan strategy comparison
The few simulations presented so far, assuming a LiteBIRDlike scanning strategy, show the effect of the HWP precession on fullsky maps and angular power spectra. Because the scan strategy can have a role in mitigating this systematic effect that couples temperature and Bmode polarization (Wallis et al. 2017), we implemented simulations, as described in Sect. 5.1, that can reproduce different satellite observational strategies. The results obtained from analyzing these simulations are reported as residual maps (i.e., difference between output and input maps). They are shown in Fig. 9 as root mean square (hereafter RMS) of the residual maps, shown in Table 2, and as Bmodes angular power spectra, shown in Fig. 10.
Fig. 9. Difference between polarization Qmaps reconstructed through a Stokes polarimeter showing the direct effect of the waveplate precession with different scan strategies. Modulation parameters: f_{s} = 0.1 Hz, θ_{0} = 1°, I_{⊥}/I_{s} = 0.508 for Planck, I_{⊥}/I_{s} = 0.514 for COrE and LiteBIRD, and I_{⊥}/I_{s} = 0.510 for WMAP. Scan strategy from the left to the right: Planck, WMAP, COrE, and LiteBIRD. 
RMS values in μK for intensity and polarization for different scan strategies.
Fig. 10. Bmode power spectra (a) recovered from simulations including the systematic and the difference output minus input of BB power spectra (b) normalized with respect to the input spectrum for the Planck, WMAP, COrE, and LiteBIRD scan strategies with different HWP spinning speeds (f_{s} = [0.1, 1] Hz). The dashed line represents the primordial Bmode angular power spectrum in the case r = 0.001, divided by the lensing only Bmodes; the other cosmological parameters are assumed from Planck Collaboration I (2018). 
In Fig. 9 we report the Q residual maps, in histogramequalized color scale, for the case f_{s} = 0.1 Hz, θ_{0} = 1° (the U maps show variations with a similar pattern and similar dynamic range). We made the simulations with several HWP physical parameters. Here we report the results for the following values of the I_{⊥}/I_{s} ratio as representative cases: 0.508 for Planck 0.514 for COrE and LiteBIRD, and 0.510 for WMAP.
The result of this analysis shows that the residual does not only depend on the scanning strategy, but mostly on the combination of scanning strategy, HWP rotation speed, and I_{⊥}/I_{s} ratio.
For example, all scanning strategy simulations show the largest effect in the case of an HWP spinning at 0.1 Hz, and they show the lowest residual in the case of 1 Hz spin frequency. This is also evident from the RMS value, reported in Table 2, and from the angular power spectrum in Fig. 10. Some strategies produce a spurious peak in the angular power spectrum that might be induced by a resonance between satellite spin and HWP wobbling, that is, COrE 0.1 Hz at ℓ ∼ 40, WMAP 0.1 Hz at ℓ ∼ 85, or LiteBIRD 0.1 Hz at ℓ ∼ 120.
The Plancklike scanning strategy (Tauber et al. 2006) does not create particular patterns or structures on larger angular scales, as can be seen in the BB power spectra that we plot in orange and blue in Fig. 10.
On the other hand, the COrElike simulation (with slowly spinning HWP) shows the worst coupling between the satellite spin and waveplate precession in terms of angular structures at large scales, as is visible in terms of spurious Bmodes (plotted in brown in Fig. 10). These simple cases show that the largescale patterns in the residual map are not related to the whole quality of the map, which is better described by the angular power spectrum.
The power spectra and the maps we recovered show the contamination that is generated by the halfwave plate precession systematic for a specified scan strategy. Repeating the analysis with different precession angles, we conclude as expected that the larger the precession angle, the stronger the spurious Bmode signal, and the higher the HWP spin frequency, the greater the mitigation of the systematic effect.
In terms of research for primordial Bmodes, the faster rotation of the wave plate helps to mitigate the systematic effect that is induced by the precession of the modulating element in a Stokes polarimeter by moving the contamination at high ℓ. Fig. 10b illustrates the fractional residual Bmodes due to observation with a wobbling HWP in the case of notensor perturbations, but only lensinginduced Bmodes. The fractional residual power is a good figure of merit of the contamination because nextgeneration CMB polarization experiments are designed to reach a sensitivity that is usually quantified as a fraction of the lensinginduced Bmodes level.
5.4. Temperaturetopolarization leakage
When we include a map with only temperature foregrounds in the simulation input, we can highlight the temperaturetopolarization leakage effect that is induced by the systematic for various scenarios. We verified that polarization foregrounds, removed with the ideal component separation method, leave a residual that is one order of magnitude lower in terms of polarizationtopolarization leakage. The HWP wobble induces Bmodes whose amplitude is proportional to θ_{0}, shown in Fig. 11, where we report the BB power spectra for different precession angles [0.5°, 1°, 1.5°, 2°] extracted from the maps scanned by a COrElike satellite. The induced Bmode signal exceeds the gravitational lensing contribution already for θ_{0} = 1°.
Fig. 11. Bmode power spectra recovered from simulations including the systematic for the COrE scan strategy with fixed HWP spinning speed (f_{s} = 0.1 Hz) and different precession angles. The input map includes only temperature contributions; the polarization contribution in the output map arises from the HWP precession. 
The output polarization components Q and U are shown in Table 3.
RMS values in μK for intensity and polarization residual maps considering a COrElike scan strategy with the HWP spin speed at 0.1 Hz and with different precession angles.
5.5. Neighboring detector to mitigate the systematic
The systematic effect induced by the wobbling can be mitigated by redundancy. Observing the same sky pixel with different phases of the wobbling plate averages the contamination out. This can also be obtained by combining multiple detectors and observing the same pixel at different times. In order to determine this mitigation, we simulated the observation with two different detectors, pointing to different boresight angles β, shifted by 1°, for all the proposed scanning strategies. The modulation parameters were f_{s} = 1 Hz, θ_{0} = 1°, I_{⊥}/I_{s} = 0.508 for Planck, I_{⊥}/I_{s} = 0.514 for COrE and LiteBIRD, and I_{⊥}/I_{s} = 0.510 for WMAP. The combination of the data from the two detectors resulted in a single map with a reduced contamination with respect to the singledetector maps, as reported in Table 4. In the table we report the RMS of the difference between the map with and without the induced systematic effect. This RMS of the residual is very similar for the two singledetector maps, and is reduced in the map we produced when they were combined. The RMS of the residuals scales with a factor , indicating that the contamination is rather uncorrelated among the two detectors.
RMS values, in μK, of the difference between the maps with and without the effect of the HWP wobbling.
6. Conclusions
We described the systematic error induced by the precession of the halfwave plate modulator in a Stokes polarimeter and the effects on CMB polarization measurements. In Sect. 3 we derived the analytical equation (Eq. (12)) for the power on the central detector of the polarimeter when the HWP precesses with an angle θ_{0} and with a frequency γ imposed by the physical properties of the device (the spin frequency f_{s} and the I_{⊥}/I_{s} ratio, including its support). Using Eq. (12), we performed several simulations to assess the level of the systematic effect that is induced by the precession. We found a strong dependence on θ_{0} and γ both for the fractional residual of the signal and for the power spectrum.
In Sect. 5 we developed the simulation of fullsky observation by a satellite mission, including the systematic effect, and quantified its effect on Bmodes that were retrieved from the output map. The HWP precession produces spurious effects on the maps at different angular scales depending on the strategy used to scan the sky. In particular, we implemented four scan strategies, WMAP, Planck, COrE, and LiteBIRDlike. With the specific mechanical properties of the implemented HWP, our simulations show a Bmode induced by leakage from the intensity signal that dominates at different scales for the strategies selected: at ℓ ∼ 40 for COrElike satellite, at ℓ ≃ 85 for a WMAPlike, and ℓ ≃ 120 in the LiteBIRDlike case. The newera CMB experiment will gain some order of magnitude in sensitivity (Abitbol et al. 2017), a few 10^{−4}μK, compared to today. The analytical formula of the systematic effect that is induced by the HWP wobbling enables removing or at least to predicting its effect on data.
In general, the precession effect induces beats in the time lines. In amplitude and frequency, these beats are related to the physical properties of the moving parts. Their effect on the maps depends on specific scanning parameters and on possible coupling with the beat frequency. We recommend that this effect be taken into account in the design of an observation strategy by modeling and measuring the inertia tensor of the moving parts. When the tensor is measured (or modeled), it can be inserted into Eq. (12) to simulate the effect on the timelines. The scanning strategy must avoid any synchrony with the beat frequency. In this case, redundancy helps to cancel most of the contamination, but considering the targeted sensitivity of future CMB experiments (Sekimoto et al. 2018), the precession must still be considered as a possible source of systematic effect.
Acknowledgments
We acknowledge the support of the ASICOSMOS programme, Elia Battistelli, for CPU time on the SPINDUST server. The work has been supported by University of Rome “Sapienza” research funds. LP acknowledges the support of the CNES postdoctoral program.
References
 Abitbol, M. H., Ahmed, Z., Barron, D., et al. 2017, ArXiv eprints [arXiv:1706.02464] [Google Scholar]
 Adam, R., Ade, P. A. R., Aghanim, N., et al. 2016, A&A, 594, A9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Ade, P. A. R., Aghanim, N., Arnaud, M., et al. 2016, A&A, 594, A13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Aumont, J., Banfi, S., Battaglia, P., et al. 2016, ArXiv eprints [arXiv:1609.04372] [Google Scholar]
 Bennett, C. L., Bay, M., Halpern, M., et al. 2003, ApJ, 583, 1 [NASA ADS] [CrossRef] [Google Scholar]
 BICEP2 Collaboration (Ade, P. A. R., et al.) 2014, Phys. Rev. Lett., 112, 241101 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 BICEP2/Keck & Planck Collaborations (Ade, P. A. R., et al.) 2015, Phys. Rev. Lett., 114, 101301 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Bryan, S. A., Montroy, T. E., & Ruhl, J. E. 2010a, Appl. Opt., 49, 6313 [NASA ADS] [CrossRef] [Google Scholar]
 Bryan, S. A., Ade, P. A. R., Amiri, M., et al. 2010b, Proc. SPIE, 7741, 77412B [CrossRef] [Google Scholar]
 Bryan, S., Ade, P., Amiri, M., et al. 2016, Rev. Sci. Inst., 87, 014501 [NASA ADS] [CrossRef] [Google Scholar]
 Chuss, D. T., Wollack, E. J., Pisano, G., et al. 2012, Appl. Opt., 51, 6824 [NASA ADS] [CrossRef] [Google Scholar]
 Columbro, F., Battistelli, E. S., Coppolecchia, A., et al. 2019, Astron. Nachr., 340, 83 [NASA ADS] [CrossRef] [Google Scholar]
 D’Alessandro, G., de Bernardis, P., Masi, S., & Schillaci, A. 2015, Appl. Opt., 54, 9269 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Das, S., & Souradeep, T. 2015, JCAP, 1505, 012 [NASA ADS] [CrossRef] [Google Scholar]
 EBEX Collaboration (Aboobaker, A. M., et al.) 2018, ApJS, 239, 7 [NASA ADS] [CrossRef] [Google Scholar]
 Errard, J., Feeney, S. M., Peiris, H. V., & Jaffe, A. H. 2016, JCAP, 3, 052 [Google Scholar]
 EssingerHileman, T., Kusaka, A., Appel, J. W., et al. 2016, Rev. Sci. Inst., 87, 094503 [Google Scholar]
 GellMann, M. 1962, Phys. Rev., 125, 1067 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Gorski, K. M., Hivon, E., Banday, A. J., et al. 2005, ApJ, 622, 759 [NASA ADS] [CrossRef] [Google Scholar]
 Gu, H., Chen, X., Zhang, C., Jiang, H., & Liu, S. 2018, J. Opt., 20 [Google Scholar]
 Hill, C. A., Beckman, S., Chinone, Y., et al. 2016, Proc. SPIE, 9914, 99142U [CrossRef] [Google Scholar]
 Inoue, Y., Ade, P., Akiba, Y., et al. 2016, Proc. SPIE, 9914, 99141I [CrossRef] [Google Scholar]
 Johnson, B. R., Columbro, F., Araujo, D., et al. 2017, Rev. Sci. Inst., 88, 105102 [NASA ADS] [CrossRef] [Google Scholar]
 Jones, R. C. 1941, J. Opt. Soc. Am., 31, 488 [Google Scholar]
 Keck Array & BICEP2 Collaborations (Ade, P. A. R.) 2016, Phys. Rev. Lett., 116, 031302 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Kleppner, D., & Kolenkow, R. 1973, An Introduction to Mechanics (Cambridge: Cambridge University Press) [Google Scholar]
 Kusaka, A., EssingerHileman, T., Appel, J. W., et al. 2014, Rev. Sci. Inst., 85, 024501 [NASA ADS] [CrossRef] [Google Scholar]
 LSPE Collaboration (Aiola, S., et al.) 2012, Proc. SPIE, 8446, 84467A [NASA ADS] [CrossRef] [Google Scholar]
 Natoli, P., Ashdown, M., Banerji, R., et al. 2018, JCAP, 1804, 022 [NASA ADS] [CrossRef] [Google Scholar]
 O’Dea, D., Challinor, A., & Johnson, B. R. 2007, MNRAS, 376, 1767 [NASA ADS] [CrossRef] [Google Scholar]
 OrtegaQuijano, N., & ArceDiego, J. L. 2013, Opt. Express, 21, 6895 [NASA ADS] [CrossRef] [Google Scholar]
 Piat, M., Battistelli, E., Baù, A., et al. 2012, J. Low Temp. Phys., 167, 872 [NASA ADS] [CrossRef] [Google Scholar]
 Pisano, G., Savini, G., Ade, P. A. R., & Haynes, V. 2008, Appl. Opt., 47, 6251 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Planck Collaboration I. (Akrami, Y., et al.) 2018, A&A, submitted [arXiv:1807.06205] [Google Scholar]
 Planck Collaboration ES. 2015, The Explanatory Supplement to the Planck 2015 Results, http://wiki.cosmos.esa.int/planckpla2015 (ESA) [Google Scholar]
 Salatino, M., de Bernardis, P., & Masi, S. 2011, A&A, 528, A138 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Salatino, M., Lashner, J., Gerbino, M., et al. 2018, Proc. SPIE, 10708, 1070848 [Google Scholar]
 Samim, M., Krouglov, S., & Barzda, V. 2016, Phys. Rev. A, 93, 013847 [NASA ADS] [CrossRef] [Google Scholar]
 Sekimoto, Y., Ade, P., Arnold, K., et al. 2018, Proc. SPIE, 10698, 106981Y [Google Scholar]
 Sheppard, C. J. R. 2011, J. Opt. Soc. Am. A, 28, 2655 [NASA ADS] [CrossRef] [Google Scholar]
 Sheppard, C. J. R. 2014, Phys. Rev. A, 90, 023809 [NASA ADS] [CrossRef] [Google Scholar]
 Sheppard, C. J. R., Castello, M., & Diaspro, A. 2016, J. Opt. Soc. Am. A, 33, 1938 [NASA ADS] [CrossRef] [Google Scholar]
 Takakura, S., Aguilar, M., Akiba, Y., et al. 2017, JCAP, 5, 008 [NASA ADS] [CrossRef] [Google Scholar]
 Tauber, J., Bersanelli, M., Lamarre, J. M., et al. 2006, ArXiv eprints [arXiv:astroph/0604069] [Google Scholar]
 Tegmark, M. 1997, ApJ, 480, L87 [NASA ADS] [CrossRef] [Google Scholar]
 Thornton, R. J., Ade, P. A. R., Aiola, S., et al. 2016, ApJS, 227, 21 [NASA ADS] [CrossRef] [Google Scholar]
 Wallis, C. G. R., Brown, M. L., Battye, R. A., & Delabrouille, J. 2017, MNRAS, 466, 425 [NASA ADS] [CrossRef] [Google Scholar]
Appendix A: Precession theory
We report in this appendix the derivation of the main equations that describe a precessing disk. In particular, we detail the derivation of Eq. (3).
The contribution to L_{x} by the rotation about the xaxis is . We can treat I_{xx} as a constant because moments of inertia about the principal axes are constant for small angular displacements. In addition, the rotation about the ydirection contributes to L_{x} by giving a component L_{s} sin (ξ) on the xaxis. Combining these contributions, we obtain
Because I_{xx} = I_{yy} = I_{⊥}, and by exploring smallangle firstorder approximation, Eq. (A.1) read
Furthermore, based on the same approximation, L_{z} = I_{s}ω_{s}. Because we consider a torquefree system (dL/dt = 0), both L_{s} and ω_{s} are constant, leading to
By introducing and , Eq. (A.3) becomes
In order to solve this coupled system of differential equations, we can differentiate one and substitute the other:
The solution for the harmonic motion is (with A and ϕ arbitrary constants)
while for ω_{y} we obtain
Integrating ω_{x} and ω_{y}, we obtain
In the smallangle approximation we impose that A/γ ≪ 1. These equations reveal that the spin axis rotates around a fixed direction in space. If that direction is along the zaxis, then θ_{x0} = θ_{y0} = 0. Assuming the initial conditions η(t = 0) = θ_{0} and ξ(t = 0)=0, and assuming that A/γ = θ_{0}, we obtain
The last equations describe the torquefree precession of the spin axis that rotates in space at a fixed angle θ_{0} with respect to the zaxis with a frequency of the precession motion given by γ = ω_{s}I_{s}/I_{⊥} (Kleppner & Kolenkow 1973).
Considering a thin disk, we obtain I_{s} = 2I_{⊥} and so γ = 2ω_{s}, which means that the disk wobbles twice as fast as it spins.
Finally, the apparent rate of a thindisk precession to an observer on the rigid body is
Appendix B: 3D Jones matrix calculation
We detail in this appendix the 3D Jones matrix description of a rotating HWP. The Jones matrices used here are those that describe the rotation about the two axes that are orthogonal to the propagation direction (with their inverse transformations):
The angles η and ξ are the precession angles η and ξ derived in Sect. 2 (Eqs. (1) and (2)).
The Jones matrix for an HWP with the fast axis at angle θ with respect to the horizontal axis is
In the end, the matrix for a linear polarizer that transmits the horizontal component of a light beam is
B.1. Nowobble case
By combining these matrices, we obtain a general formula for the Jones vector at the onaxis detector of a polarimeter with modulating element precesses:
These calculations were realized with the Python package SymPy^{1}, which allows symbolic computations with matrices and vectors.
For the ideal case with η = ξ = 0°, when no precession occurs, we obtain the outgoing Jones vector for an ideal polarimeter:
So that the intensity is
where we recognize the modulation terms, at four times the frequency of the HWP rotation, within the expressions cos^{2}(2θ), sin^{2}(2θ) and cos(2θ)sin(2θ). Passing through the Stokes parameters (T, Q, U), we retrieve the intensity at the same detector for an ideal polarimeter:
where the Stokes components T, Q, U are then defined as usual:
B.2. Wobbled case
For a general polarimeter with a wobbling HWP, we obtain a general formula for the Jones vector at the onaxis detector:
For the general case with η, ξ ≠ 0°, we find the intensity at the detector:
where we define the modulating functions g(ξ, θ) and f(η, ξ, θ):
We can write the outgoing intensity through the Stokes parameters as follows:
The last equation gives the intensity at the onaxis detector of a Stokes polarimeter with its modulating element describing a torquefree precession. It is not merely a function of the HWP orientation about the ẑ axis, but also depends on the displacements from the x̂ and ŷaxes due to the precession.
Appendix C: 3D Müller matrix calculation
The Müller formalism is required to propagate partially polarized light (like that of the CBM), therefore we need to calculate Eqs. (B.7) and (B.12) using the Müller formalism. 3D Müller matrices M_{ij} are related to Jones matrices (Samim et al. 2016) by
where J is the associated Jones matrix, and σ_{n} (n = [0, …, 8]) are the GellMann matrices. Eq. (C.1) is valid if and only if we use tracenormalized GellMann matrices (GellMann 1962; Sheppard et al. 2016), which are defined as follows:
Using the common polarization matrix
we can define the Stokes vector in the 3D formalism:
The conventional 2D Stokes parameters are related to the 3D Stokes parameters (optical ordering) by
C.1. Nowobble case
From Eq. (C.1) we can find the analogous 3D Müller matrices for the Jones HWP matrix (Eq. (B.2)) and for the Jones polarizer matrix (Eq. (B.3)),
and by combining the previous matrices, we can determine the 3D polarimeter Stokes vector (η = ξ = 0°):
Through Eq. (C.5), we can derive the intensity:
Equations (C.9) and (B.7) give the same result.
C.2. Wobbled case
For a wobbling HWP we need to calculate the Müller rotation matrices (Eq. (B.1)) from Eq. (C.1):
Finally, we can calculate the 3D Müller matrix for a Stokes polarimeter with a wobbling HWP:
For the sake of clarity, we define the following equations that we use from now on:
By assuming that the field that enters the polarimeter has E_{z} = 0, this is true only for the onaxis detectors, and for all the focal plane assuming a telecentric optic system, we can find the general Stokes vector for a wobbling HWP:
where we used , and through Eq. (C.5) the equivalent intensity
The components of the output Stokes vector that are not null (Eq. (C.24)) are , , and , so that from the definition of the Stokes vector (Eq. (C.4)), it is clearly . This is the case because the polarizing grid does not permit for onaxis rays. If we were to lose this assumption (i.e., for an offaxis detector), it would be easy to verify that the output Stokes vector become a function of Δ_{0}, Δ_{1}, Δ_{2}, Δ_{4}, Δ_{5}, and Δ_{8}, but the components of the output Stokes vector that are not null are always , , and (Eq. (C.25)).
All Tables
RMS values in μK for intensity and polarization residual maps considering a COrElike scan strategy with the HWP spin speed at 0.1 Hz and with different precession angles.
RMS values, in μK, of the difference between the maps with and without the effect of the HWP wobbling.
All Figures
Fig. 1. Illustration of rigidbody precession where we show the angular momentum dislocations with respect to the x̂ and ŷ axes. 

In the text 
Fig. 2. (a) Fractional residuals between the total intensities in a Stokes polarimeter with and without HWP precession. We fix the ratio I_{⊥}/I_{s} = 0.502 and consider different HWP precession angles θ_{0} = [1°,0.5°,0.25°], normalizing the difference with respect to the maximum intensity of an ideal Stokes polarimeter. (b) Fractional residuals, as before, but in a shorter time interval to highlight the beats induced by precession. (c) Power spectra with peaks at multiples and submultiples of the HWP spinning frequency with an amplitude related to θ_{0}. 

In the text 
Fig. 3. Fractional residual for different ratios I_{⊥}/I_{s} and different HWP spinning frequencies. Each column shows a different f_{s} value and each row a different I_{⊥}/I_{s} ratio. A thinner HWP, I_{⊥}/I_{s} → 0.500, has beats in the intensity over an extended period and a thicker HWP has shorter beats. By looking at the figure from left to right, it is clear that the spinning speed compresses the beats. The maximum amplitude remains the same because it only depends on θ_{0}, which is fixed to 1° in this particular simulation. 

In the text 
Fig. 4. Power spectra extracted from time lines showing the shifting of the spurious peak at low frequencies. The simulation considers input light that is horizontally polarized, a precession angle θ_{0} = 1°, and I_{⊥}/I_{s} = [0.501, 0.502, 0.506, 0.513]. 

In the text 
Fig. 5. Top panel: T, Q, and U maps reconstructed through a Stokes polarimeter where the HWP wobbles. Bottom panel: difference between input and output T, Q, and U maps showing the effect of the waveplate precession on observed maps. 

In the text 
Fig. 6. Temperature and polarization difference maps showing the effect of the waveplate precession for a slowspinning modulator. Modulation parameters: f_{s} = 0.1 Hz, θ_{0} = 1°, I_{⊥}/I_{s} = 0.514. 

In the text 
Fig. 7. Input minus output relative difference of EE and BB power spectra computed from the maps shown in Fig. 5. 

In the text 
Fig. 8. Panels a and c: BB and EE power spectra extracted from the input and output maps of Fig. 6; in panel c the two lines overlap. Panels b and d: normalized difference between input and output power spectra. 

In the text 
Fig. 9. Difference between polarization Qmaps reconstructed through a Stokes polarimeter showing the direct effect of the waveplate precession with different scan strategies. Modulation parameters: f_{s} = 0.1 Hz, θ_{0} = 1°, I_{⊥}/I_{s} = 0.508 for Planck, I_{⊥}/I_{s} = 0.514 for COrE and LiteBIRD, and I_{⊥}/I_{s} = 0.510 for WMAP. Scan strategy from the left to the right: Planck, WMAP, COrE, and LiteBIRD. 

In the text 
Fig. 10. Bmode power spectra (a) recovered from simulations including the systematic and the difference output minus input of BB power spectra (b) normalized with respect to the input spectrum for the Planck, WMAP, COrE, and LiteBIRD scan strategies with different HWP spinning speeds (f_{s} = [0.1, 1] Hz). The dashed line represents the primordial Bmode angular power spectrum in the case r = 0.001, divided by the lensing only Bmodes; the other cosmological parameters are assumed from Planck Collaboration I (2018). 

In the text 
Fig. 11. Bmode power spectra recovered from simulations including the systematic for the COrE scan strategy with fixed HWP spinning speed (f_{s} = 0.1 Hz) and different precession angles. The input map includes only temperature contributions; the polarization contribution in the output map arises from the HWP precession. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.