Issue |
A&A
Volume 627, July 2019
|
|
---|---|---|
Article Number | A160 | |
Number of page(s) | 16 | |
Section | Astronomical instrumentation | |
DOI | https://doi.org/10.1051/0004-6361/201834495 | |
Published online | 22 July 2019 |
Systematic effects induced by half-wave 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
e-mail: 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. Paris-Sud, Université Paris-Saclay, 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 B-mode 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 half-wave 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 B-mode 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 B-modes 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 time-line of the electric fields. Finally, assuming that the entire sky is observed with a satellite mission, we analyze the effect on B-mode measurements.
Conclusions. The effect is not negligible at the sensitivity of current B-mode 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 B-modes in 2014 (BICEP2 Collaboration 2014) by measuring the tensor-to-scalar 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 tensor-to-scalar 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 self-calibrating 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 half-wave 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 temperature-to-polarization leakage, chromatic HWP behavior, and scanning strategy, have been evaluated in Essinger-Hileman 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 full-sky observation of the CMB, assuming a satellite mission with different scanning strategies.
2. Precessing half-wave 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 Ls = Isωs along the symmetry axis, where Is and ωs are the moment of inertia and the angular velocity, respectively. In the unperturbed case, Ls coincides with the ẑ axis.
![]() |
Fig. 1. Illustration of rigid-body 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 Ls, expanding it in its components along the x̂ and ŷ axes, η and ξ. By introducing the precession rate γ, defined as the rotation velocity of Ls around the ẑ axis, η and ξ are simply defined as
These equations describe the torque-free precession of a spin axis that rotates in space with fixed angle θ0.
The precession rate is given by γ = ωsIs/I⊥ (Kleppner & Kolenkow 1973), after Is 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 Is = 2I⊥ and γ = 2ωs, which means that the disk wobbles twice as fast as it spins.
The apparent rate of a thin-disk precession for an observer on the rigid body is given by
In a torque-free 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 so-called 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 so-called 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 two-dimensional vector i = (Ex, Ey) and optical elements that change the state of input radiation with 2 × 2 matrices. We extend the Jones vector in three dimensions, i = (Ex, Ey, 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; Ortega-Quijano & Arce-Diego 2013). Combining these matrices, we obtain a general formula for the Jones vector for the on-axis 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 T2 = Q2 + U2. The last equation gives the intensity for the on-axis detector of a Stokes polarimeter whose modulating element describes a torque-free 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 trace-normalized Gell-Mann 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 Ez = 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 mij 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 well-known 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⊥/Is = 0.502. The Jones vector used as input field in the simulation was iin = (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 non-precessing 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⊥/Is = 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⊥/Is 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 fs and to the ratio I⊥/Is (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 fs = [0.1, 0.5, 1.0, 2.0] Hz and I⊥/Is = [0.501, 0.502, 0.506, 0.513]. As an example, these values for the I⊥/Is 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 Is is fixed to 0.0128 kg m2. 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⊥/Is and different HWP spinning frequencies. Each column shows a different fs value and each row a different I⊥/Is ratio. A thinner HWP, I⊥/Is → 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⊥/Is → 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⊥/Is. The spectra exhibit the effect discussed above, showing that the beat frequency moves to lower values as the ratio I⊥/Is is reduced toward the minimum value of 0.5 (for I⊥/Is = 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⊥/Is = [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 POLARBEAR-2 (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 LSPE-SWIPE (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. Full-sky 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 time-lines. We completed this software with a map-making algorithm that collapses data time-lines into maps. All simulations are noise-free 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 Sun-Earth 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 fs frequency. The systematic affecting the HWP was included in the data at the time-line level, and a simple rebinning map-making 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 full-sky 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 temperature-to-polarization 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, free-free, spinning dust, CO, and thermal dust emission, without considering their polarization contribution.
This map was modeled in order to highlight the temperature-to-polarization 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 Nside = 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 Nside = 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 Is/I⊥ ratio. For each simulation we compared input and output maps and computed the B-mode power spectrum. As a visual example, we show in Fig. 5 (top panel) the output maps for a mission for which we adopted a LiteBIRD-like 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⊥/Is = 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 wave-plate 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 histogram-equalized color scale, is at the same level as in the input map.
![]() |
Fig. 6. Temperature and polarization difference maps showing the effect of the wave-plate precession for a slow-spinning modulator. Modulation parameters: fs = 0.1 Hz, θ0 = 1°, I⊥/Is = 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 B-modes (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 B-mode 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 LiteBIRD-like scanning strategy, show the effect of the HWP precession on full-sky maps and angular power spectra. Because the scan strategy can have a role in mitigating this systematic effect that couples temperature and B-mode 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 B-modes angular power spectra, shown in Fig. 10.
![]() |
Fig. 9. Difference between polarization Q-maps reconstructed through a Stokes polarimeter showing the direct effect of the wave-plate precession with different scan strategies. Modulation parameters: fs = 0.1 Hz, θ0 = 1°, I⊥/Is = 0.508 for Planck, I⊥/Is = 0.514 for COrE and LiteBIRD, and I⊥/Is = 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. B-mode 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 (fs = [0.1, 1] Hz). The dashed line represents the primordial B-mode angular power spectrum in the case r = 0.001, divided by the lensing only B-modes; the other cosmological parameters are assumed from Planck Collaboration I (2018). |
In Fig. 9 we report the Q residual maps, in histogram-equalized color scale, for the case fs = 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⊥/Is 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⊥/Is 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 Planck-like 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 COrE-like simulation (with slowly spinning HWP) shows the worst coupling between the satellite spin and wave-plate precession in terms of angular structures at large scales, as is visible in terms of spurious B-modes (plotted in brown in Fig. 10). These simple cases show that the large-scale 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 half-wave 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 B-mode signal, and the higher the HWP spin frequency, the greater the mitigation of the systematic effect.
In terms of research for primordial B-modes, 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 B-modes due to observation with a wobbling HWP in the case of no-tensor perturbations, but only lensing-induced B-modes. The fractional residual power is a good figure of merit of the contamination because next-generation CMB polarization experiments are designed to reach a sensitivity that is usually quantified as a fraction of the lensing-induced B-modes level.
5.4. Temperature-to-polarization leakage
When we include a map with only temperature foregrounds in the simulation input, we can highlight the temperature-to-polarization 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 polarization-to-polarization leakage. The HWP wobble induces B-modes 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 COrE-like satellite. The induced B-mode signal exceeds the gravitational lensing contribution already for θ0 = 1°.
![]() |
Fig. 11. B-mode power spectra recovered from simulations including the systematic for the COrE scan strategy with fixed HWP spinning speed (fs = 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 COrE-like 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 fs = 1 Hz, θ0 = 1°, I⊥/Is = 0.508 for Planck, I⊥/Is = 0.514 for COrE and LiteBIRD, and I⊥/Is = 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 single-detector 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 single-detector 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 half-wave 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 fs and the I⊥/Is 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 full-sky observation by a satellite mission, including the systematic effect, and quantified its effect on B-modes 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 LiteBIRD-like. With the specific mechanical properties of the implemented HWP, our simulations show a B-mode induced by leakage from the intensity signal that dominates at different scales for the strategies selected: at ℓ ∼ 40 for COrE-like satellite, at ℓ ≃ 85 for a WMAP-like, and ℓ ≃ 120 in the LiteBIRD-like case. The new-era 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 ASI-COSMOS 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 e-prints [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 e-prints [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]
- Essinger-Hileman, T., Kusaka, A., Appel, J. W., et al. 2016, Rev. Sci. Inst., 87, 094503 [Google Scholar]
- Gell-Mann, 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., Essinger-Hileman, 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]
- Ortega-Quijano, N., & Arce-Diego, 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 e-prints [arXiv:astro-ph/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 Lx by the rotation about the x-axis is . We can treat Ixx as a constant because moments of inertia about the principal axes are constant for small angular displacements. In addition, the rotation about the y-direction contributes to Lx by giving a component Ls sin (ξ) on the x-axis. Combining these contributions, we obtain
Because Ixx = Iyy = I⊥, and by exploring small-angle first-order approximation, Eq. (A.1) read
Furthermore, based on the same approximation, Lz = Isωs. Because we consider a torque-free system (dL/dt = 0), both Ls 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 small-angle 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 z-axis, 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 torque-free precession of the spin axis that rotates in space at a fixed angle θ0 with respect to the z-axis with a frequency of the precession motion given by γ = ωsIs/I⊥ (Kleppner & Kolenkow 1973).
Considering a thin disk, we obtain Is = 2I⊥ and so γ = 2ωs, which means that the disk wobbles twice as fast as it spins.
Finally, the apparent rate of a thin-disk 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. No-wobble case
By combining these matrices, we obtain a general formula for the Jones vector at the on-axis detector of a polarimeter with modulating element precesses:
These calculations were realized with the Python package SymPy1, 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 cos2(2θ), sin2(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 on-axis 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 on-axis detector of a Stokes polarimeter with its modulating element describing a torque-free 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 Mij are related to Jones matrices (Samim et al. 2016) by
where J is the associated Jones matrix, and σn (n = [0, …, 8]) are the Gell-Mann matrices. Eq. (C.1) is valid if and only if we use trace-normalized Gell-Mann matrices (Gell-Mann 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. No-wobble 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 Ez = 0, this is true only for the on-axis 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 on-axis rays. If we were to lose this assumption (i.e., for an off-axis 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 COrE-like 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 rigid-body 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⊥/Is = 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⊥/Is and different HWP spinning frequencies. Each column shows a different fs value and each row a different I⊥/Is ratio. A thinner HWP, I⊥/Is → 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⊥/Is = [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 wave-plate precession on observed maps. |
In the text |
![]() |
Fig. 6. Temperature and polarization difference maps showing the effect of the wave-plate precession for a slow-spinning modulator. Modulation parameters: fs = 0.1 Hz, θ0 = 1°, I⊥/Is = 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 Q-maps reconstructed through a Stokes polarimeter showing the direct effect of the wave-plate precession with different scan strategies. Modulation parameters: fs = 0.1 Hz, θ0 = 1°, I⊥/Is = 0.508 for Planck, I⊥/Is = 0.514 for COrE and LiteBIRD, and I⊥/Is = 0.510 for WMAP. Scan strategy from the left to the right: Planck, WMAP, COrE, and LiteBIRD. |
In the text |
![]() |
Fig. 10. B-mode 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 (fs = [0.1, 1] Hz). The dashed line represents the primordial B-mode angular power spectrum in the case r = 0.001, divided by the lensing only B-modes; the other cosmological parameters are assumed from Planck Collaboration I (2018). |
In the text |
![]() |
Fig. 11. B-mode power spectra recovered from simulations including the systematic for the COrE scan strategy with fixed HWP spinning speed (fs = 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 (full-text 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 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.