Alassane Sidibe
email: alassane.sidibe@uwinloc.com

Laura López Mir

Benjamin Dhuiège
email: benjamin.dhuiege@genesink.com

Gaël Depres
email: gael.depres@arjowiggins.com

Alexandru Takacs
email: atakacs@laas.fr

Jan Mennekens

A Thin Paper UHF Antenna on Nanocelloluse Based Substrate for Battery-free Geolocation Tags

This paper presents the design and characterisation of an antenna on a paper-based substrate. The antenna is designed on a nanocellulose foil optimised with a highly conductive ink based on silver nanoparticles. The electrical and dielectric properties of the materials were characterised prior to the antenna simulation. The antenna is designed to operate at the ISM 868 MHz frequency band and has a maximum gain of +1.7 dBi. Its development is part of the first step in the implementation of a flexible battery-free geolocation tag, remotely powered by UHF Energy sources in the European RFID spectrum between 865.6 MHz and 867.6 MHz. An accurate characterisation procedure was carried out and the results show good agreement with the simulations.

I. INTRODUCTION

The trend of new technologies leads to the development and production of several devices on a Printed Circuit Board (PCB), the manufacturing process of which can sometimes be complicated and expensive. There is therefore a drive to reduce costs in mass applications. In addition, with the need to recycle devices, make them smaller and thinner, electronics manufacturers are making many improvements to their processes and adapting them especially to the new generation of Internet of Things (IoT) devices in the Organic and Large Area Electronics (OLAE) branch. For these reasons, Uwinloc [START_REF]UWINLOC | New Indoor Location System for assets tracking[END_REF], a company that develops and commercialises an indoor location solution using a battery-free tag, decided to participate in the European project Madras [START_REF]MADRAS -Advanced materials and processing in organic electronics[END_REF] in order to meet the new challenges of the new markets requirements. Several research activities have been carried out previously on flexible antennas [START_REF] Kirtania | Flexible Antennas: A Review[END_REF]. In [START_REF] Thielens | Fabrication and characterization of flexible spray-coated antennas[END_REF], a Radiofrequency Identification (RFID) antenna on Polyethylene Naphthalate (PEN) substrate with a spray coating method is presented. Our proposal focuses on the design of an antenna with innovative materials. This antenna will cover the European Industrial, Scientific, and Medical (ISM) frequency band (865.6-867.6 MHz) to harvest the radiated electromagnetic (EM) waves by an RFID source. In section II, the materials developed for the antenna are presented and the simulation and characterisation procedure is presented in section III. II. ANTENNA MATERIALS For energy harvesting, the antennas must have a high efficiency and an omnidirectional radiation pattern. The radiation efficiency of the antenna depends on the properties of the conductive ink, the dielectric, and the mismatch losses.

The materials were developed by considering the requirement (printed layer surface resistance lower than 100 mΩ/sq) for a competitive antenna design. Thus, it is necessary to have a high conductive ink.

A. Development of the conductive ink

The formulated ink called Smart Screen B (S-CS91349) is a sustainable conductive ink that has been developed for the printed electronics market and is particularly well suited for applications requiring high conductivity and high resolution on paper substrates. This ink, based on silver nanoparticles (Ag-NP), is perfectly adapted to design conductive tracks on flexible substrates and is suitable to produce antennas for IoT applications (RFID, NFC) using screen-printing technologies. It was developed by Genesink [START_REF]Conductive ink for printed electronics : IOT, display and more[END_REF]. The sheet resistance was measured with a 4-point probe and gives a low resistivity around 20 mΩ/sq for a 2.4 μm thick layer. Therefore, this ink satisfies the specifications required for antenna development. The conductivity of this formulated ink is higher than the standard ink used on flexible or paper-based antennas as in [START_REF] Thielens | Fabrication and characterization of flexible spray-coated antennas[END_REF], [START_REF] Phan | Novel ultra-wideband "robe" antenna on high-loss paper substrate[END_REF]. Its volume resistivity is between 4 and 5.5 μΩ.cm, which corresponds to a skin depth of 3.9 μm at 900 MHz.

B. Paper-based substrate

The advantage of using an innovative nano-cellulose based substrate is to give customers the opportunity to target new applications and develop smart systems for other use-cases. The paper substrate used is the antenna design was produced by Arjowiggins [7] with a certain amount (60%) of cellulose nanofibrils (CNF), which has a transparency of 79 %. This paper is made with pure cellulose, which is enzymatically treated, with endoglucanase, then mechanically treated by refiners. The resulting 60-µm thick Nano-cellulose Foil (NCF) substrate was characterised with cavity resonator method and has a dielectric constant of 4.2 (±0.2%) at 900 MHz and a dissipation factor of 0.135 (±0.015%) at 1 kHz.

III. ANTENNA DESIGN AND CHARACTERISATION

A. Simulation of the proposed antenna

In order to validate the developed materials dedicated for the final battery-free geolocation tag, a modified loop antenna is designed (Fig. 1 (a)). The 3D EM simulation and optimisation of the antenna were performed with Ansys HFSS software. The loop antenna structure was optimised after several simulations steps to improve the return loss and the realized gain. It is simulated with the obtained parameters (dielectric and electrical) of the Ag NP ink and the NCF substrate. The PCB and the connector to be used for the characterisation were clearly modelled in the software to faithfully reproduce the fabricated samples. The simulated return loss of the antenna in Fig. 3 shows a resonant frequency at 855 MHz with a bandwidth of 67 MHz of under -10 dB. From the radiation pattern and the 3D polar plot in Fig. 2, the antenna exhibits an omnidirectional radiation pattern at 867 MHz with a maximum simulated gain of +1.9 dBi.

B. The antenna characterisation

The design was carried out at Eurecat Technology Centre [START_REF]Eurecat Centre Tecnològic -Innovant per a les empreses[END_REF] with screen-printed the Ag-NP ink on the NCF substrate. The nominal thickness of 1-layer ink coated on NCF is 2.7 μm, which is less than the skin depth. However, the power losses caused by the skin effect have been evaluated on the simulation before. They are only 2% of efficiency or 0.25 dBi gain less compared to a 9 μm thick ink. The characterisation procedure consists of using a small PCB with two pads on each side: on one side, a U.FL connector is welded while the other side is connected and fixed to the printed antenna with a conductive and a structural adhesive. Three samples were evaluated in an anechoic chamber at the LAAS-CNRS laboratory, where the antennas are placed under free space conditions to measure the performance in terms of the radiation efficiency as expected from the simulations. They were also placed on a Rohacell foam (εr ≃ 1; tanD ≃ 0.001) to emulate the software simulation condition. As shown in Fig. 3, the measured return loss of the 3 samples fits well, but there is a slight shift of 14 MHz (at the peak resonant frequency) from the simulation. This difference is probably due to the U.FL to SMA cable used. The measured gain of each sample also corresponds to the simulated value with a maximum of +1.6 dBi (±0.1) in E-plane at 867 MHz.

IV. CONCLUSION In this paper, an antenna on NCF foil with a highly conductive Ag-NP operating in the narrow band of the European UHF Spectrum allocated to the RFID application is presented. The antenna achieves a maximum gain of +1.7 dBi in the Eplane at 867 MHz. It has the advantage of being on a thin, flexible, and low-cost substrate. It is also dedicated to being recyclable after the end of use. The manufacturing of the antenna does not require a complicated process (screen-printing) and an accurate characterisation can be performed with a low-cost U.FL connector. This preliminary design has validated the properties of materials developed by the partners. The authors are currently working on an optimised and more compact design that could be used with a small PCB module allowing RF-to-dc conversion, power management and transmission of an Ultra Wide Band (UWB) geolocation signal.

Figure 1 .

 1 Figure 1. (a) Geometries of the designed antenna; (b) Representation of the excitation strategy with a hybridized PCB with U.FL connector.

Figure 2 .

 2 Figure 2. Simulation of the realized gain at 867 MHz: (a) Radiation pattern plot at φ = 0° (red line) and φ =90° (Blue dashed line); (b) 3D polar plot of the gain.

Figure 3 .

 3 Figure 3. Simulated and measured return loss of the 3 samples of the fabricated antenna.

Figure 4 .

 4 Figure 4. (a) Measurement setup on an anechoic chamber, (b) Comparison of the simulated and measured gain in the E-plane (φ=90°) at 867 MHz.

ACKNOWLEDGMENT

The authors would like to acknowledge all the partners in the Madras project [2]. This project has received funding from the European Union's Horizon 2020 research and innovation program under grant agreement No. 862492.


Ranajoy Banerji
email: ranajoyb@astro.uio.no

Guillaume Patanchon
email: guillaume.patanchon@apc.univ-paris-diderot.fr

Jacques Delabrouille
email: delabrouille@apc.in2p3.fr

Masashi Hazumi
email: masashi.hazumi@kek.jp

Duc Thuong Hoang
email: hoang@apc.in2p3.fr

Hirokazu Ishino
email: scishino@s.okayama-u.ac.jp

Tomotake Matsumura
email: tomotake.matsumura@ipmu.jp

Bandpass mismatch error for satellite CMB experiments II: correcting for the spurious signal

Keywords: CMBR experiments, CMBR polarisation, inflation

Future Cosmic Microwave Background (CMB) satellite missions aim at using the B-mode polarisation signal to measure the tensor-to-scalar ratio r with a sensitivity σ(r) of the order of ≤ 10 -3 . Small uncertainties in the characterisation of instrument properties such as the spectral filters can lead to a leakage of the intensity signal to polarisation and can possibly bias any measurement of a primordial signal. In this paper we discuss methods for avoiding and correcting for the intensity to polarisation leakage due to bandpass mismatch among detector sets. We develop a template fitting map-maker to obtain an unbiased estimate of the leakage signal and subtract it out of the total signal. Using simulations we show how such a method can reduce the bias on the observed B-mode signal by up to 3 orders of magnitude in power.

Introduction

The Cosmic Microwave Background (CMB) is the relic radiation left over from the epoch when electrons and light nuclei from the primordial plasma first combined to form neutral atoms. In the ΛCDM cosmological scenario, this happened when the Universe was about 380,000 years old. As the properties of the CMB anisotropies depend on the physics that gave rise to initial seeds of structure and also on particle interactions in the primordial plasma, the CMB is a rich source of information about the early Universe. The Planck mission, launched by ESA in 2009 [START_REF] Matsumura | Planck early results. I. The Planck mission[END_REF], delivered to the scientific community full-sky, signal-dominated primordial CMB temperature anisotropy maps on all scales down to about 5 resolution, and provided the state-of-the art measurement of CMB polarization of the E-mode type (of even parity) [START_REF]Planck 2013 results. I. Overview of products and scientific results[END_REF][START_REF]Planck 2015 results. I. Overview of products and scientific results[END_REF]. These observations put strong constraints on the ΛCDM cosmological scenario and on its main parameters [START_REF]Planck 2013 results. XVI. Cosmological parameters[END_REF].

However, the observation of one of the remaining key CMB observables, polarization of the B-mode type (of odd parity), has only just barely started. First detections of such Bmodes due to the conversion of E-modes from the last scattering surface by intervening matter have been reported [START_REF]Measurement of the cosmic microwave background polarization lensing power spectrum with the POLARBEAR experiment[END_REF][START_REF]Measurements of sub-degree B-mode polarization in the cosmic microwave background from 100 square degrees of SPTpol data[END_REF][START_REF]BICEP2 and Keck Array collaborations, Improved constraints on cosmology and foregrounds from BICEP2 and Keck Array cosmic microwave background data with inclusion of 95 GHz band[END_REF][START_REF]The Atacama Cosmology Telescope: two-season ACTPol spectra and parameters[END_REF][START_REF]A measurement of the cosmic microwave background B-mode polarization power spectrum at sub-degree scales from 2 years of POLARBEAR data[END_REF]. However, the very interesting primordial B-modes, predicted to be generated in the early Universe in a phase of rapid expansion known as cosmic inflation, still escape detection. The expected amplitude of such primordial B-modes is set by the value of JCAP07(2019)043 a single parameter, r = ∆ 2 t /∆ 2 s (∆ is the amplitude of the fluctuation), which measures the ratio of tensor perturbations to scalar perturbations in the early universe at some reference scale. As of now, combined observations from the BICEP2 and Keck Array experiments and from the Planck mission yield an upper limit of r < 0.06 at 95% confidence [START_REF]BICEP2 and Keck Array collaborations, Improved constraints on cosmology and foregrounds from BICEP2 and Keck Array cosmic microwave background data with inclusion of 95 GHz band[END_REF][START_REF]2 and Keck Array collaborations, BICEP2/Keck Array x: constraints on primordial gravitational waves using Planck, WMAP and new BICEP2/Keck observations through the 2015 season[END_REF]. Proposed future CMB space missions such as the consecutive versions of CORE [START_REF]COrE (Cosmic Origins Explorer) a white paper[END_REF][START_REF]Exploring cosmic origins with CORE: survey requirements and mission design[END_REF], LiteBIRD [START_REF] Hazumi | LiteBIRD: a small satellite for the study of B-mode polarization and inflation from cosmic background radiation detection[END_REF][START_REF] Matsumura | Mission design of LiteBIRD[END_REF][START_REF] Ishino | LiteBIRD: lite satellite for the study of B-mode polarization and inflation from cosmic microwave background radiation detection[END_REF], PICO [START_REF] Sutin | PICO -the probe of inflation and cosmic origins[END_REF][START_REF] Young | Optical design of PICO, a concept for a space mission to probe inflation and cosmic origins[END_REF], PIXIE [START_REF] Kogut | The primordial inflation explorer (PIXIE): a nulling polarimeter for cosmic microwave background observations[END_REF], PRISM [START_REF]PRISM (Polarized Radiation Imaging and Spectroscopy Mission): an extended white paper[END_REF] or ground-based experiments such as CMB-S4 [START_REF]CMB-S4 collaboration[END_REF][START_REF]CMB-S4 collaboration, CMB-S4 technology book[END_REF] target a sensitivity σ(r) of 10 -3 or better, requiring instrumental sensitivity and control of instrumental systematic effects four orders of magnitude better than for observing CMB temperature anisotropies on large scales.

The measurement of polarization B-modes is a challenging task. Polarization-sensitive detectors used in CMB mapping typically observe a mixture of intensity and polarization, in which the contribution from the sky intensity outshines the primordial polarization B-modes by three to four orders of magnitude, or possibly more. To recover B-mode polarization, intensity signals must be separated from polarization using linear combinations of measurements observing with different polarization angles. Such linear combinations can be obtained either using several observations with the same polarization-sensitive detector re-observing the sky at the exact same location but with different orientations, or using different polarizationsensitive detectors observing with different polarization angles [START_REF]Exploring cosmic origins with CORE: mitigation of systematic effects[END_REF]. In this second case, a mismatch of the spectral response between different detectors, with different polarization angles, leads to a leakage of the dominant intensity signal to polarization.

A companion paper develops the data model for the leakage of signal from intensity to polarization due to bandpass mismatch between detectors, and evaluates the level of the bias on the observed B-mode signal from the CMB sky, assuming a similar level of uncertainty in band properties as in Planck HFI [START_REF] Hoang | Bandpass mismatch error for satellite CMB experiments I: estimating the spurious signal[END_REF]. It shows that if the bandpass mismatch is not taken care of by appropriate means, the bias on the B-modes at low multipoles can be of the same order of magnitude as the primordial signal on large scale for a target sensitivity of σ(r) 10 -3 . In principle, an ideal and achromatic fast rotating Half Wave Plate (HWP) can be used to mitigate this potential systematic effect, and the development of near-ideal broad-band HWPs is actively pursued. Here, we investigate the effectiveness of an alternative approach to mitigate, in the map-making step, the leakage of signal from intensity to polarization due to bandpass mismatch among a set of detectors when no HWP is used in the observation. We extend the preliminary work done in the context of the CORE systematic effects study outlined in section 6 of [START_REF]Exploring cosmic origins with CORE: mitigation of systematic effects[END_REF]. The study of whether it is possible to correct the impact of HWP imperfections is beyond the scope of this paper.

The rest of this paper is organised as follows. In section 2, we model the bandpass mismatch systematic effect for a set of polarization sensitive detectors and discuss how the mismatched signal projects onto the polarization maps. Section 3 is dedicated to the estimation of the noise penalty that occurs when one tries to avoid the bandpass mismatch by making single detector polarisation maps which are combined afterwards, or equivalently by fitting for a pixel-dependent correction term for each detector pair. In section 4, we detail the development of map-making approaches that estimate and correct for the leakage signal to first order. We then demonstrate the effectiveness of the correction method using simulations of data timestreams for a set of detectors and a model of bandpass mismatch leakage due to different foreground contamination of the data stream in section 5 and discuss the results. In section 6 we conclude and give an outlook for the correction method for use in future CMB missions and for including various other systematic contamination effects.

JCAP07(2019)043 2 Model of the bandpass leakage systematic effect

Most of the observations of CMB polarisation are obtained using total-power detectors that integrate radiation in wide frequency bands (with typical bandwidth ∆ν/ν 0.25 -0.35), along a given linear polarization direction. Focal plane pixels can comprise pairs of associated detectors in which two polarimetric sensors (labelled a and b) are sensitive to orthogonal components of the mean squared electric field integrated in the band. Such detector pairs, implemented in practice using polarization sensitive bolometers (PSBs), have been used in the High Frequency Instrument (HFI) of the Planck space mission [START_REF] Lamarre | Planck pre-launch status: the HFI instrument, from specification to actual performance[END_REF]. Multichroic versions using Transition Edge Sensors (TESs), in which each focal plane pixel comprises six detectors, observing each in one of two orthogonal polarisations and one of three independent frequency bands, are being deployed in ground-based experiments [START_REF] Suzuki | Multi-chroic dual-polarization bolometric focal plane for studies of the cosmic microwave background[END_REF][START_REF] Posada | Large arrays of dual-polarized multichroic TES detectors for CMB measurements with the SPT-3G receiver[END_REF], and are foreseen as the baseline for the low frequency telescope in the LiteBIRD space mission [START_REF] Suzuki | The LiteBIRD satellite mission -sub-Kelvin instrument[END_REF].

The frequency band of each detector is defined primarily by a set of stacked filters and on-chip filter banks to the independent radiation, each with a high-pass, low-pass, or band-pass spectral response to the incident radiation. Optical elements such as corrugated horns and lenses also contribute to the selection of photons, as well as does the absorption by the sensors themselves. Overall, each detector i is characterised by a specific spectral response g i (ν). Although all the detectors in a given frequency band are designed to have the same spectral response, small variations from detector to detector are unavoidable in practice, and have been observed in the case of Planck HFI, as illustrated in [START_REF]Planck 2013 results. IX. HFI spectral response[END_REF], Figs 5-10. As discussed in the companion paper [START_REF] Hoang | Bandpass mismatch error for satellite CMB experiments I: estimating the spurious signal[END_REF], and detailed below, this induces a bandpass mismatch between the detector bands, which leads to a mismatch between the contributions of the various astrophysical emission components in different detector timestreams after the various detectors have been calibrated on CMB temperature anisotropies, or the dipole.

Band-integrated observations

For each detector i, the observed band-averaged intensity in the direction of sky pixel p can be written as

I i (p) = c ν dν g i (ν) I c (ν, p), (2.1)
where the sum runs over all components c (such as the CMB, galactic synchrotron, free-free, thermal dust emission, . . . ), g i (ν) is the normalised spectral response of detector i such that dνg i (ν) = 1, and I c (ν, p) is the intensity Stokes parameters of component c at frequency ν and in pixel p. Similar equations hold for Stokes parameters Q and U . This model of band-averaged intensity can be recast as

I i (p) = c γ ic (p) I c (ν 0 , p), (2.2)
where γ ic (p) is a bandpass coefficient that depends on the detector i being considered, the component c and, if the frequency scaling of component c is pixel-dependent, also on the pixel of the sky p. This formalism is particularly relevant when the emission I c (ν, p) of each component can be written as a product of a pixel-dependent template I c (ν 0 , p) at a reference frequency ν 0 , and scaled over frequencies with a single pixel-independent scaling factor A c (ν, ν 0 ):

I c (ν, p) = A c (ν, ν 0 ) I c (ν 0 , p). (2.3)

JCAP07(2019)043

For all components for which this is the case, γ ic does not depend on the pixel p, and we have

γ ic = ν dν g i (ν) A c (ν, ν 0 ). (2.4)
To a very good approximation, this is the case for CMB temperature and polarisation anisotropies, and for kinetic and thermal Sunyaev-Zel'dovich (kSZ and tSZ) effects (neglecting the relativistic corrections that modify the latter when interactions with very hot electron gas are considered) [START_REF] Challinor | Relativistic corrections to the Sunyaev-Zel'dovich effect[END_REF][START_REF] Itoh | Relativistic corrections to the Sunyaev-Zel'dovich effect for clusters of galaxies[END_REF]. The emission laws of other astrophysical components, such as Galactic synchrotron, free-free, and thermal dust emission, are only approximately pixel-independent, with varying degree of approximation. The approximation of pixel-independence will be used in section 4 to develop the method to correct, to first order, polarisation maps from bandpass mismatch leakage. To give a specific example, the emission law for galactic thermal dust is well modelled using a modified blackbody of the form

A dust (ν) ∝ ν β B ν (T dust ), (2.5)
where B ν (T ) is the Planck blackbody function for temperature T dust , and β is a spectral index that depends on the chemical composition and structure of dust grains. Fitting such a model on multi-frequency data from the Planck space mission, both T dust and β are found to be pixel-dependent, with T dust ranging from about 15 K to about 27 K, and β from about 1.2 to about 2.2, with mean values of 1.62 and 19.6 K for the spectral index and the temperature respectively for intensity maps [START_REF]Planck intermediate results. XXII. Frequency dependence of thermal emission from galactic dust in intensity and polarization[END_REF]. As argued in [START_REF] Martínez-Solaeche | A 3D model of polarized dust emission in the milky way[END_REF], this range of temperatures and spectral indices, although compatible with Planck data analysis, may be underestimated by reasons of line-of-sight integration effects.

In the companion paper [START_REF] Hoang | Bandpass mismatch error for satellite CMB experiments I: estimating the spurious signal[END_REF], the bandpass mismatch error is estimated by assuming the pixel-independence of γ ic . Although this is sufficient to evaluate the amplitude of the bandpass mismatch leakage to first order, we will relax this assumption in the present paper. This will allow us to estimate the amplitude of residuals at second order, after first order correction of the bandpass leakage systematic effect, in section 5. We proceed in the rest of the section assuming the pixel-dependence of A c and hence γ ic .

The different detectors can be intercalibrated on a strong signal with pixel-independent emission law, such as the CMB dipole and additionally the CMB anisotropies using a blind independent component analysis (ICA) method such as SMICA [START_REF] Delabrouille | Multi-detector multi-component spectral matching and applications for CMB data analysis[END_REF][START_REF] Cardoso | Component separation with flexible models. Application to the separation of astrophysical emissions[END_REF], as demonstrated on an analysis of WMAP multi-frequency data [START_REF] Patanchon | CMB and foregrounds in WMAP first year data[END_REF]. This amounts to dividing equation (2.2) by the γ i coefficient of CMB temperature anisotropies, and redefining the units of intensity. We can then write

I i (p) = I CMB (p) + c = CMB γ ic (p) I c (ν 0 , p). (2.6)
Note that, by calibrating over the CMB dipole or temperature anisotropies (and thus expressing the I, Q and U Stokes parameters of the observed sky brightness in thermodynamic units of temperature fluctuations of a 2.725 K blackbody), we normalise the bandpass parameter γ i for the CMB to unity. For the other components c, the bandpass parameter γ ic encompasses all of the dependence of the signal upon the band spectral response function g i (ν) for detector i, the spectral emission law of each individual component c and the choice of the reference frequency ν 0 . Since unit conversion coefficients are frequency-dependent, the frequency-dependence of A c (ν, ν 0 , p) and hence the value of γ ic (p) both depend on the choice of units for I c (ν, p).

JCAP07(2019)043

Using the model of equation (2.2) for the various astrophysical components, the signal for detector i at time index t is

X it = c γ ic (p t ) [I c (p t ) + Q c (p t ) cos (2ψ it ) + U c (p t ) sin (2ψ it )] + n it , (2.7)
where p t is the sky pixel pointed to at time t, X is the total band integrated noisy signal, ψ is the orientation of the polarisation angle in the local coordinate frame at point p t on the sky, and n it is the noise of detector i at time t.

Intensity leakage term due to bandpass mismatch

For a set of detectors {i} in a given frequency band, small variations in the actual spectral response g i (ν) between the detectors result in slightly different γ ic parameters, which can be modelled as

γ ic (p) = γ c (p) + δγ ic (p), (2.8)
where γ c is a reference value (e.g., for a reference frequency band) common to all detectors for the component c, and δγ ic (p) is a small detector-dependent and pixel-dependent variation.

We can rewrite the signal (2.7) of detector i, keeping the bandpass variation term coupled to the intensity, as:

X it = [I(p t ) + Q(p t ) cos (2ψ it ) + U (p t ) sin (2ψ it )] + c = CMB δγ ic (p t ) I c (p t ) + n it , (2.9)
where the term in square brackets is the total sky signal in the ideal frequency band, and the second term is the sum of bandpass leakages from intensity due to astrophysical components. The CMB does not appear in the summation as we assume that the data is perfectly calibrated on a CMB anisotropy signal. As the intensity of the various sky emission sources is at least an order of magnitude larger than their polarised amplitude, we neglect in this model the leakage terms from Q and U Stokes parameters. We thus decompose our data model into two parts. One containing the (multicomponent) sky signal we are interested in measuring, and the other, the spurious leakage term which projects foreground intensity onto the polarisation maps when we solve for I, Q and U using a set of different detector measurements in any pixel p. This spurious term is what we wish to eliminate so that it does not project onto the final polarisation maps.

Bandpass mismatch error for a detector pair

In many existing and planned CMB polarisation experiments, detectors are arranged in pairs of orthogonal polarimeters, which we index with i a and i b to represent orthogonal detectors in a pair. Measuring conventionally the angle ψ it as that of detector i a , the half difference of the two measurements from two orthogonal detectors is:

X i [a-b] t = [Q cos (2ψ it ) + U sin (2ψ it )] + c = CMB ∆γ ic I c + n i [a-b] t , (2.10)
where ∆γ ic = (γ iac -γ i b c )/2 is the bandpass mismatch coefficient between detectors i a and i b for component c. We have dropped the explicit pointing notation p t since we are in time domain. The subscript i [a-b] t conventionally indicates that the set of observations X i [a-b] t is obtained from the half difference of X iat and X i b t . The first term on the right hand side

JCAP07(2019)043

of equation (2.10) is the useful polarisation signal, the second term is the spurious bandpass mismatch contamination of polarization by intensity, and the noise is the half difference of the individual detector noise terms. Such data stream differences are useful to immediately eliminate the intensity term from the equations, and working directly with only polarisation and bandpass mismatch leakage terms.

It is interesting to note that the polarisation signal on the sky does not contribute to any leakage due to a bandpass mismatch, for detector pairs. If we had considered the spurious Q and U terms in equation (2.10), they would have appeared with a coefficient proportional to (δγ iac + δγ i b c ) and which is 0. We demonstrate this using simulations in appendix A figure 5.

Bandpass mismatch error in Stokes parameters maps

Assume now that a set of 2N i detectors, arranged in orthogonal pairs {i a , i b }, is used to map sky intensity and polarisation in a given frequency channel. Starting from the signal data model in equation (2.7), we can rewrite it in matrix-vector format, as

X = PS + n, (2.11)
where X is the time-ordered vector of observations of all the 2N i detectors, S = [I, Q, U ] T is the vector of Stokes parameters in pixel space, P is a sparse 'pointing matrix' constituted of as many rows as there are observations and as many columns as pixels, for which each non-zero element in a row is of the form [1, cos(2ψ it ), sin(2ψ it )], and n is the noise. The best least square error estimate of the Stokes parameters S is

S = [P T C -1 n P] -1 P T C -1 n X (2.12)
where C n is the noise covariance matrix. If errors between the measurements are uncorrelated and all have the same variance, C n is proportional to an identity matrix, and the solution becomes S = [P T P] -1 P T X.

(2.13)

If the observations are made with sets of 'optimised configurations' of polarimeters with orientations evenly spread in [0, π], as often is the case, the map-domain noise covariance matrix of all three Stokes parameters, [P T P] -1 , is diagonal [START_REF] Couchot | Optimized polarimeter configurations for measuring the Stokes parameters of the cosmic microwave background radiation[END_REF], and the Q and U Stokes parameters in any pixel p are estimated as

Q(p) ∝ it cos (2ψ it ) X it (2.14) and U (p) ∝ it sin (2ψ it ) X it . (2.15)
where the sum over the time index is done for those corresponding to pixel p. When all detectors are arranged in orthogonal pairs, taking into account that cos (2ψ iat ) = -cos (2ψ i b t ) and sin (2ψ iat ) = -sin (2ψ i b t ), the spurious bandpass leakage term in (2.10) projects on the Q and U maps as errors δ Q and δ U .

δ Q(p) = 1 N hits (p) c I c (p) i∈{ia} ∆γ ic (p) t cos (2ψ it ), (2.16a) δ U (p) = 1 N hits (p) c I c (p) i∈{ia} ∆γ ic (p) t sin (2ψ it ), (2.16b)

JCAP07(2019)043

where the sum over i is over all detector labelled a in the pairs used for making the polarisation maps, the sum over t is for all observations of pixel p with detector pair i at different times, and N hits (p) is the total number of observations by detector pairs in pixel p. The leakage in either map is proportional to the bandpass mismatch parameter ∆γ. If we consider one single pair i of orthogonal detectors, the leakage of intensity in the estimated Q map is proportional to t cos(2ψ it ), and in the U map is proportional to t sin(2ψ it ). As already pointed-out in the companion paper, an observing strategy for which these average terms vanish results in vanishing projection of the intensity leakage, due to bandpass mismatch, onto polarization maps. The bandpass leakage term can significantly contaminate B-mode polarisation maps when these maps are reconstructed from the data streams as described above, using data from a set of detectors with mismatched frequency bands. In the next section, we discuss ways to avoid projecting the bandpass mismatch term onto Stokes parameters maps, and the penalty in terms of noise level in the case where the measurements are not obtained with evenly spread polarimeter angles.

3 Map-making accounting for bandpass mismatch

Single detector map-making

In principle, it is possible to completely avoid bandpass mismatch contamination by making single detector maps of I, Q, and U , and then combining these single detector maps into multidetector maps by simple averaging the polarisation maps (possibly with pixel-dependent coefficients between the maps, for noise-weighting as a function of the pixel-dependent noise level). This is possible as long as the set of angles ψ it allows inverting, for each detector i independently, the linear system that connects Stokes parameters to observations. Observing strategies that make use of an ideal rotating HWP, which would spread the polarimeter angles evenly in [0, π] for each sky pixel without generating HWP-specific systematics, would allow for such single-detector map-making with no loss of sensitivity. At the other extreme, an observing strategy with no HWP, and in which some sky pixels are observed always with the same scanning orientation, does not allow for disentangling I, Q and U for those pixels, and thus does not allow for single-detector map-making. Moreover, depending on the scanning strategy, in the absence of a rotating HWP, it may happen that a large spread of angles can only be achieved after a long period of observation (e.g. one year), which requires, to avoid other possible systematic effects, the stability of the instrument over such long periods. We consider the intermediate case of an observation strategy with no HWP, but in which each pixel is observed with different scanning orientations. Specifically, we consider the case where the spacecraft is located at the Sun-Earth L2 Lagrange point and spins around a spin axis which itself precesses around the anti-solar direction. Such a scan strategy has been considered for CORE [START_REF]Exploring cosmic origins with CORE: survey requirements and mission design[END_REF], LiteBIRD [START_REF] Matsumura | Mission design of LiteBIRD[END_REF] and PICO [START_REF] Sutin | PICO -the probe of inflation and cosmic origins[END_REF]. The choice of the angle between the spin axis and the line of sight (β), between the precession axis and the spin axis (α), and of the periods of precession and of spinning, has an influence on the sky coverage and on the distributions of scanning angles over the sky. A more complete discussion of the trade-offs and of the impact of the parameters of the scanning strategy on the projection of systematic effects can be found in reference [START_REF] Wallis | Optimal scan strategies for future CMB satellite experiments[END_REF].

For a given scanning strategy, single detector maps are obtained by solving, for each detector j and each pixel p, a linear system of equations such as (2.7) for a set of values of t. The linear system of equations linking the observations to the sky Stokes parameters is JCAP07(2019)043 inverted with equation (2.12), or with (2.13) when only uncorrelated white noise is present. This generates no leakage of intensity into polarisation.

However, the map-making strategy that consists in doing first single detector polarisation maps, and then combining those to produce maps of Q and U (hence, without contamination of polarisation maps by intensity signal from non-CMB components) is in general sub-optimal in terms of noise projection. Indeed, the noise covariance matrix of the I, Q and U Stokes parameters is

P T C -1 n P -1 = σ 2 n N hits    1 cos 2ψ jt sin 2ψ jt cos 2ψ jt cos 2 2ψ jt sin 4ψ jt 2 sin 2ψ jt sin 4ψ jt 2 sin 2 2ψ jt    -1 . (3.1)
where σ 2 n is the variance of the noise in the detector timestreams, and N hits the total number of data points falling in pixel p. Only when each detector i observes each pixel with sets of at least three independent measurements with angles ψ jt evenly spread in [0, π] so that cos 2ψ jt = sin 2ψ jt = 0 and cos 2 2ψ jt = sin 2 2ψ jt = 1/2, is the noise covariance matrix of the I, Q and U Stokes parameters diagonal with minimal volume of the error box [START_REF] Couchot | Optimized polarimeter configurations for measuring the Stokes parameters of the cosmic microwave background radiation[END_REF], so that we have

P T C -1 n P -1 = σ 2 n N hits    1 0 0 0 2 0 0 0 2    . (3.2)
This ideal case can be achieved with an ideal HWP to modulate the polarisation. For a nonideal HWP the bandpass may depend on the HWP angle, and the map-making would require again appropriate measures to avoid projecting bandpass mismatch errors in a similar way as in equation (2.9). Without a HWP a near-ideal spread of scanning angles can be achieved with a spacecraft that rotates around the telescope line of sight, but some unevenness cannot be avoided for off-axis detectors.

In figure 1, we illustrate, for four different scanning strategies, the BB auto power spectrum of the noise projected onto the polarisation maps. For comparison, we also show the theoretical ideal-case noise level that would be achieved if all of the observation time was evenly spread among all pixels, and scanning angles evenly spread in [0, π] for each pixel. Depending on the scan strategy, the noise power spectrum (and hence the white noise variance) after single-detector map-making can be very sub-optimal. This noise increase is due to the large fraction of the sky being sampled poorly leading to a large fraction of the sky having a covariance matrix far from the ideal case, which, for a scan strategy combining rotation and precession, increases with increasing opening angle between the spin axis and the line of sight (spin angle β). In general, a scan strategy with large precession angle, α, and small spin angle, β, (see figure 1 in reference [START_REF] Hoang | Bandpass mismatch error for satellite CMB experiments I: estimating the spurious signal[END_REF]) results in more evenly spread scanning angles in [0, π], and hence smaller noise boosting with single-detector map-making (as compared to map-making using all the detectors with no bandpass mismatch between them).

Leakage subtraction with pair by pair map-making

In WMAP data analysis, timelines consisting of differences between two detectors are used to make maps of I, Q and U . Bandpass mismatch between detectors, if neglected, would result in leakages of intensity signal onto polarisation maps. To avoid this contamination, the WMAP team used a map-making strategy that consists in estimating jointly the three JCAP07(2019)043 Stokes parameters I, Q, U as well as an extra leakage parameter M [START_REF]Three-year Wilkinson Microwave Anisotropy Probe (WMAP) observations: beam profiles, data processing, radiometer characterization and systematic error limits[END_REF][START_REF]Three year Wilkinson Microwave Anisotropy Probe (WMAP) observations: polarization analysis[END_REF]. This yields the same errors on Q and U maps than single-detector map-making followed by map-averaging.

Working directly with timelines of differences between two orthogonal detectors {i a , i b } in a pair, we can write equation (2.8) for the two detectors in the form

γ iac (p) = γ ic (p) + δγ ic (p) (3.3) and γ i b c (p) = γ ic (p) -δγ ic (p), (3.4)
where γ ic (p) is the average of γ iac (p) and γ i b c (p) in pixel p, and δγ ic (p) the half-difference of the two. From equation (2.9), noting that ψ it = ψ iat = (ψ i b t -π/2), and denoting as M i (p) the leakage parameter

M i (p) = c = CMB δγ ic I c (p), (3.5)
for the pair {i a , i b }, we have

X iat (p) [I i (p) + Q(p) cos (2ψ it ) + U (p) sin (2ψ it )] + M i (p) + n iat , (3.6)
and

X i b t (p) [I i (p) -Q(p) cos (2ψ it ) -U (p) sin (2ψ it )] -M i (p) + n i b t . (3.7)
Note that the quantity I i (p) represents the mean intensity for the considered pair, and hence varies from detector pair to detector pair as does the averaged bandpass. This also is the

JCAP07(2019)043

case for M i which represents the half-difference of intensities I ia and I i b . This dependence is made explicit with index i. 1Solving for [I, Q, U, M ], for well-balanced detector noise variance, we get

      I i Q U M i       =       1 0 0 0 0 cos 2 2ψ it sin 4ψ it 2 cos 2ψ it 0 sin 4ψ it 2 sin 2 2ψ it sin 2ψ it 0 cos 2ψ it sin 2ψ it 1       -1       X i [a+b] t (p) X i [a-b] t (p) cos 2ψ it X i [a-b] t (p) sin 2ψ it X i [a-b] t (p)       (3.8)
where, as in equation (2.10), we denote as X i [a-b] t (p) the half difference of the two measurements from the two orthogonal detectors at time t in pixel p, and as X i [a+b] t (p) the half sum (quantities that can be computed directly as half-difference and half-sum data streams when the beams of detectors i a and i b are co-extensive). The noise covariance matrix of the reconstructed I, Q, U and M is

Cov 4 = σ 2 n N hits       1 0 0 0 0 cos 2 2ψ it sin 4ψ it 2 cos 2ψ it 0 sin 4ψ it 2 sin 2 2ψ it sin 2ψ it 0 cos 2ψ it sin 2ψ it 1       -1
.

(3.9)

We see straightforwardly in equations. (3.8) and (3.9) that the half-sum timestream contributes only to the measurement of I, while the half-difference is used to estimate Q, U , and the additional leakage term M . From a direct comparison of the noise covariance matrix for Q, U and M with that of equation (3.1), we also notice that this method produces the same Q and U maps as single detector map-making followed by co-averaging the maps obtained by individual detectors. Hence, the same noise penalty incurs on the reconstruction of Q and U .

In the absence of bandpass leakage that requires estimating M (p) (or equivalently estimating two different intensities for each detector pair, one for detector i a and one for detector i b ), the inverse noise covariance matrix for Q and U would be instead

Cov 3 = σ 2 n N hits    1 0 0 0 cos 2 2ψ it sin 4ψ it 2 0 sin 4ψ it 2 sin 2 2ψ it    -1 . (3.10)
In addition to errors on I being decorrelated from those on Q and U , the average error on polarisation maps is smaller in that case.

Figure 2 shows the histograms of the noise variance for the Q and U components (diagonal elements of P T C -1 n P ), obtained on simulations with Healpix maps with Nside of 256, and for different scanning strategy parameters, after map-making with a pair of detectors, in cases where the leakage component is estimated or neglected (i.e., solving or not for M (p), and hence obtaining noise covariance matrices of either Cov 4 or Cov 3 ). When solving for M , we observe a visible excess of noise for the Q component for all scanning strategies (projected in the Ecliptic coordinate system which allow certain symmetries in the scan strategy). In the cases with α > β and precession period of 96.18 minutes, the excess is of the order of a JCAP07(2019)043 few percent, and the distribution for 3 components shows a steeper tail at the high value end. For the strategy with α < β, the excess is larger, and for some of the pixels the polarization components are poorly constrained. This degradation of the covariance when the leakage component is estimated is due to sub-optimal coverage of polariser angles in each pixel, in particular in a large region near the ecliptic equator.

We note that the U component is not degraded when solving for M , as seen from comparing the two distributions in the figures. This is due to the particular choice of coordinate systems, by reason of symmetries of the scanning strategies in ecliptic coordinates. The study of the global penalty with a full focal plane analysis, and power spectrum estimation with optimal noise weighting of the various map pixels, is postponed to future studies.

First-order bandpass mismatch correction by regression

We now wish to devise a map-making scheme that corrects for bandpass mismatch while paying the minimum penalty in terms of noise boosting. To do so, we adopt the model developed in equation (2.10), but drop the dependence of ∆γ ic on the pixel p, i.e. we approximate γ ic and hence ∆γ ic with its average value on the part of sky of interest for CMB polarisation science. This can be done because (to first order) the frequency scaling of astrophysical components does not vary much from pixel to pixel. For each detector i, the pixel-independent term γ ic is the bandpass coefficient for an average frequency scaling, given by equation (2.4).

Without precise knowledge of the individual detector bands g i (ν), the bandpass mismatch coefficients ∆γ ic between detector pairs cannot be computed accurately a priori. They must then be estimated from the data themselves. To this effect, we assume that after the mission data is analysed to produce intensity maps, and after component separation for intensity has been performed (see [START_REF] Delabrouille | Diffuse source separation in CMB observations[END_REF] for a review), a template of each of the main sky components I c (ν 0 , p) at some suitable reference frequency in the frequency band of interest (or, up to a redefinition of ∆γ ic that encompasses a factor for the average colour correction for the band and the average emission law of component c, a template of the component at any suitable reference frequency ν 0 given by (2.3)) will be available at high signal-to-noise ratio. This is, in fact, already the case for the main astrophysical components that have been seen by the Planck and WMAP space missions. Such templates can be produced, for instance, with the Planck Sky Model (PSM) [START_REF] Delabrouille | The pre-launch Planck sky model: a model of sky emission at submillimetre to centimetre wavelengths[END_REF], or downloaded from the Planck legacy archive at ESA,2 or from NASA's Legacy Archive for Microwave Background Data Analysis. 3 The quality of these component intensity templates will still drastically increase with the next CMB space mission.

Assuming then that time-domain templates I c (derived from the template maps) are available for each of the main components that induce intensity to polarisation leakage by bandpass mismatch, we can recast the data-model adopted from (2.10) with pixel independent bandpass parameters, in a matrix-vector format as,

X = PS + T∆γ + n, (4.1)
where the unknowns are S, the linear polarisation Stokes parameters Q and U in the frequency band of interest for all pixels in the sky, and ∆γ, the bandpass mismatch parameters for each detector pair {i a , i b } and for each component c. As usual n is the noise. The known quantities are the pointing matrix P and the template matrix T which is composed of the individual templates I c . Solving/inverting equation (4.1) is a regression problem where we are trying to estimate S and ∆γ given our knowledge of X, P and T. Minimizing

JCAP07(2019)043

χ 2 = [X -PS -T∆γ] T C -1 n [X -PS -T∆γ] . (4.2)

JCAP07(2019)043

with respect to S and ∆γ yields the maximum likelihood estimators

S = P T C -1 n F T P -1 P T C -1 n F T X, (4.3) ∆γ = T T C -1 n F P T -1 T T C -1 n F P X, (4.4)
where F P = 1 -P P T C -1 n P -1

P T C -1 n , (4.5)
and

F T = 1 -T T T C -1 n T -1 T T C -1 n , (4.6)
with 1 the identity matrix. The F P and F T can be thought of as filtering operators which de-project the component of the signal in the space spanned by P and T respectively. In the case where the noise is well balanced between both detectors and is uncorrelated, C n ∝ 1 and these filtering operators become F P = 1 -P P T P -1 P T and F T = 1 -T T T T -1 T T respectively. Similar approaches to map-making for CMB observations while solving for nuisance terms in this way have been used in other contexts [START_REF] Keihanen | A maximum likelihood approach to the destriping technique[END_REF][START_REF] Poletti | Making maps of cosmic microwave background polarization for B-mode studies: the POLARBEAR example[END_REF][START_REF]Planck intermediate results. XLVI. Reduction of large-scale systematic effects in HFI polarization maps and estimation of the reionization optical depth[END_REF].

The implementation of the solution of (4.3) in the form mentioned above is computationally demanding, consisting of large dense matrices and redundant values. This problem can be bypassed, without requiring sophisticated data distribution and compression algorithms, by performing the map-making in two steps. We first make the unbiased estimation of the set of amplitudes of the templates, ∆γ as in (4.4), and then estimate S as:

S = P T P -1 P T X -T ∆γ . (4.7)
We note that in (4.5) the term P T P -1 P T is simply the map-making equation (written here without inverse noise weighting C -1 n as we consider the white-noise case). The operator P then re-projects this map back into the time domain. The implementation of the estimation of ∆γ therefore does not require large computational resources, and re-utilises a standard maximum likelihood map-maker. This two step method of estimating the Stokes parameters S is equivalent, for our current purpose, to the direct estimation given by (4.3). The algorithm has been tested on toy models in A and B to demonstrate that it gives mathematically expected results.

In summary, the correction technique for the spurious bandpass leakage signal that is proposed here first estimates the spurious bandpass leakage signal, for each of the contributing sky emission components, and for each orthogonal detector pair. It then subtracts out the estimated spurious signal at the timestream level to de-bias the estimation of the polarisation Q and U Stoke's parameter maps, before re-projecting the corrected data streams onto maps of Q and U in a map-making step.

Simulations and results

Simulated data

We evaluate the effectiveness of the correction method described in section 4 using simulated data streams for 8 pairs of polarization sensitive detectors arranged in 4 rows on the focal plane. Two such orthogonal pairs of detectors follow each other on the scan path, are rotated at 45 • relative to each other, and up to a time-translation, they form an optimal configuration.

JCAP07(2019)043

We consider a scan strategy that is similar to those considered for proposed missions such as LiteBIRD, CORE and PICO. We pick a precession angle α = 45 • , and opening angle for the centre of the focal plane of β = 50 • . The spin and precession periods for these particular simulations are 10 and 96.18 minutes respectively, which corresponds closely to the LiteBIRD baseline strategy. The observation time is 1 complete sidereal year. The exact opening angle β + δβ of the detectors themselves depend on their location on the focal plane and the value of δβ for the 4 rows of 4 detectors each are 9.875 , 19.75 , 29.625 and 39.5 .

We perform our simulation and correction separately for two frequency bands centred near 80 GHz (close to the foreground minimum, where dust and synchrotron emission are comparable in amplitude) and 140 GHz (one of the preferred frequency for CMB observations), For the 80 GHz band we consider separately three cases of components, 1) CMB and Thermal Dust, 2) CMB and Synchrotron, 3) CMB, Thermal Dust and Synchrotron. All 80 GHz maps contain the three Stokes parameters I, Q, U and are smoothed by a Gaussian beam of 40 FWHM. For the 140 GHz band we consider thermal dust as the only foreground component since the contribution from synchrotron is small. These maps are smoothed by a Gaussian beam of 20 FWHM.

To simulate the effect of unequal bandpasses, the foreground components are integrated over top-hat bands that are unique for each detector. A mean bandwidth of 25% of the band central frequency is chosen and the band edges of the individual detectors are sampled from a uniform distribution centred on the mean of the band edges and with a width of 1% of the band centre frequency. This generates the same uncertainty in the bandpass parameter as seen in Planck HFI. In µK CMB units, this is equivalent to rescaling the component template maps by the bandpass function γ ic . The configuration of the detectors on the focal plane as well as the frequency bands are shown in table 1. In the simulated maps, we retain the pixel dependence of the spectral parameters (spectral index and dust temperature for thermal dust, and spectral index for synchrotron) of the foreground components, thus making the bandpass function pixel dependent. All input foreground maps are simulated using the PSM tool using Planck post-flight templates. For the Synchrotron emission we assume a power law scaling, the emission template given by remazeilles2014 [START_REF] Remazeilles | An improved source-subtracted and destriped 408 MHz all-sky map[END_REF], and the spectral index template given by mamd2008 [START_REF] Miville-Deschenes | Separation of anomalous and synchrotron emissions using WMAP polarization data[END_REF] with a mean of -3.0 and a standard deviation of 0.06. We note that these values are slightly lower than those quoted by the S-PASS survey [START_REF] Krachmalnicoff | S-PASS view of polarized galactic synchrotron at 2.3 GHz as a contaminant to CMB observations[END_REF] conducted over half the sky. For thermal dust we use the ffp10 model. This is a greybody emission with a mean greybody temperature of 19.4 K and standard deviation 1.24 K, and spectral index with a mean of 1.61 and standard deviation 0.11.

For all sets of simulations, we consider the same noise properties for each detector, a white noise with a rms of 40 µK √ s, uncorrelated between detectors. The thermal dust templates that are used in equations (4.3) to (4.6) are derived from a single thermal dust intensity map generated at 353 GHz with a delta bandpass. The timeordered template is generated individually for each detector pair using the pair's pointing and orientation data. Similarly, for synchrotron, the individual templates for each detector pair are derived from a single intensity map generated at 30 GHz with a delta bandpass. Both maps are generated using a delta bandpass at their specified frequency. Before generating the time-ordered templates, the foreground template maps are smoothed to the resolution of the frequency band they are correcting for, 20 and 40 FWHM for 140 GHz and 80 GHz bands respectively. For a future mission, these intensity maps for the components are expected to be measured at a very high signal to noise ratio but will still have uncertainties. To simulate these uncertainties we overlay a white noise map of 5 µK.arcmin rms., a reasonable assumption for any future mission of interest.

JCAP07(2019)043

Simulation results

We carry out the correction on the contaminated time-ordered signal using the iterative algorithm described in section 4. We first estimate, using equation (4.4), the amplitude of the foreground templates that fit the leakage signal pair by pair. We then subtract out the estimated leakage signal and project the corrected pair differenced signal onto polarization Q and U maps. We note that the correction algorithm approximates the emission laws for the foreground components as uniform on the sky, and hence the bandpass mismatch parameter, ∆γ, to be pixel independent. The estimated amplitude of the templates, ∆γ, is a global quantity calculated for each detector pair and each foreground component. This allows us a first order correction of the bandpass mismatch signal, with a residual leakage left behind due to the spatial variation of the bandpass mismatch term ∆γ. Since we are interested in the polarisation signal away from the galactic plane, we impose a 25% mask on the sky and the template amplitude is estimated only in the valid region.

To observe and analyse the projection of the leakage alone, we subtract out the projected noise at the map level. The map-making process, being linear, allows us this freedom to independently make noise maps. To produce the map of the leakage signal, we first make the uncorrected polarization map for the multi-detector set. The ideal, leakage less, map is that produced from single detectors individually and then averaged. The leakage projection map is then simply the multi detector map with the ideal map subtracted out of it. Similarly, we compute maps of the residual leakage left after correction by subtracting from the corrected multi-detector polarization map the ideal averaged map.

JCAP07(2019)043

Figure 3. We compare here, for the Stokes Q map, the leakage projected on to the map (left panel) to the residual leakage left after performing the correction (right panel). The first three rows are at 80 GHz and the fourth at 140 GHz. All cases were done for the complete 16 detector set in the band configuration given in 1 and the scanning strategy described in 5.1. Keeping the same colour scale, we see a significant reduction in the leakage to the polarization maps for all cases. The residual after correction can be attributed to the pixel-independent approximation that the algorithm assumes to estimate the template amplitudes.

JCAP07(2019)043 6 Conclusion and discussion

In the absence of an 'ideal' HWP, the bias from the spurious polarisation signal can effectively prevent a detection or reduce the significance of detection for a cosmological signal of r ∼ 10 -3 . This paper has one driving motivation, to find analysis algorithms which can diminish or effectively cancel completely the leakage of signal from intensity to polarisation due to bandpass mismatch systematics in future CMB space missions that will not employ a continuously rotating HWP.

In the absence of a continuously rotating HWP, the effect of bandpass mismatch can be avoided altogether if we produce single detector polarisation maps and combine them afterwards, but this solution is suboptimal in terms of noise projection. The spread of polarisation angles in a sky pixel, and hence the noise covariance matrix, depends upon the scanning strategy employed. For intermediate types of scanning strategies envisaged for future missions, the noise covariance matrix only approaches the optimal case when the precession angle is greater than the opening angle of the instrument. In other cases, a suboptimal spread of polarisation angles boosts the white noise significantly. Moreover, for scans of a shorter period (less than a year), an insufficient spread of angles might prevent inverting the data model altogether. We show that jointly estimating the Stokes parameters I, Q, U along with an additional term for the bandpass leakage, for pair differencing experiments, is equivalent to making single detector maps and averaging them. This incurs the same noise penalty which increases with increasing opening angle of the scan strategy.

To obtain an unbiased estimate of the polarisation Q and U Stokes parameters, without the noise penalty of single detector map-making, we develop a multi-detector map-maker that filters out the intensity to polarisation leakage at the timestream level. We model the leakage term as an additional spurious term proportional to the intensity signal of each foreground component, and the bandpass mismatch parameter for each detector pair. A direct estimation of the unbiased Stokes parameters is computationally intensive and hence we break down the problem in two. First, we estimate the amplitudes of the foreground templates, for each pair of orthogonal detectors. Using this, we subtract out the estimated leakage signal at the timestream level and perform our multi-detector map-making. Using this method we demonstrate that we can reduce the bias on the B mode power by about 3 orders of magnitude. The residual leakage after correction is primarily due to the nonuniformity of the foreground spectral parameters on the sky.

The filtering map-making technique demonstrated above is versatile in correcting for any systematic effect that can be modelled as an additional spurious term derived from a template. Effects such as pointing mismatch and beam ellipticity can be modelled in a similar way. We note that for the analysis of real data, refinements to the present algorithm might be required to take into account low-frequency noise, correlated noise among detectors, as well as possible correlations in the bandpass mismatch between similar detector pairs. These extensions of our algorithm are left for future work.

JCAP07(2019)043

A Spurious signal sources and residuals

We discuss here the validity and impact of some of the hypotheses we have made in designing our data model and our correction algorithm. In particular, we test the effectiveness of the correction algorithm under different conditions of noise and foreground contamination.

An important assumption in designing the correction algorithm is that the foreground emission parameters, such as spectral indices, are constant over the entire sky. In practice, these parameters are spatially variable, and the correction algorithm can only achieve a global fit for each parameter. This allows us only to perform a leading order estimation of the leakage and we expect to be left with residues after correction. We construct a few test cases under ideal conditions of spectral parameters to show that our correction algorithm performs at the theoretical limit it is designed for.

We simulate our toy model skies at the frequencies for detectors 1a and 1b in table 1 and in the same focal plane configuration. All frequency maps including the template skies, required for constructing our time-ordered templates, are obtained using the PSM tool. For the ideal case, we scale the foreground templates using spatially uniform spectral parameters: a modified black-body model for thermal dust emission with a spectral index β dust = 1.62 and dust temperature T dust = 19.6 K, and a power law model for synchrotron emission, with spectral index of β sync = -3.1 (in Rayleigh-Jeans brightness temperature units).

Our first test is designed to evaluate the contribution of the leakage from Intensity to Polarisation as compared to the cross-Polarisation leakage, and the effect of the correction on the leakage from either source. The test is done at the 140 GHz band and the input maps contain only CMB and the dominant thermal dust emission. We perform our test with three setups, the first where the input maps contain the Stokes I alone, the second containing Stokes Q and U and the third with all Stokes I, Q and U parameters. This allows us to determine exactly the source of each leakage component and the residual leakage component. In all three cases we perform our correction using the template constructed from the 353 GHz, uniformly scaled, thermal dust intensity map. Here, we de not add noise to the time stream data nor to the dust template map.

In figure 5 we see a subdominant cross-polar leakage signal in the case where we had just Q and U as input. This can possibly be attributed to the fact that Healpix pixelised polarisation maps suffer from a parallel transport effect of the polarisation angles when the map resolutions are altered or when they are beam smoothed, but they are unlikely to be cross-polar leakage as discussed in section 2.3. Let us call this an 'un-resolved' source of leakage. From the same figure we see that in the case where we have only Stokes I in the input, the residual leakage after correction can be attributed to machine precision error. This shows us that under ideal noiseless conditions when the sky emission scales uniformly over the sky, the algorithm is capable of performing an exact unbiased estimate of the leakage signal. Since the algorithm corrects only for intensity to polarisation leakage, the 'un-resolved' cross-polar leakage is left behind after correction in the case where the input contains all the Stokes parameters.

B Multiple foreground sources and correction

Our second test concerns the correction of the leakage when there are more than one mismatched foreground components. For this purpose we build our test case at 80 GHz where the intensity of thermal dust and synchrotron are comparable to each other. As before, we do not include noise in the timestream or the template map.

JCAP07(2019)043

Figure 5. We compare here the intensity to polarisation versus the cross-polar bandpass leakage to the polarisation spectrum (BB), for a single orthogonal detector pair at 140 GHz, and the only contaminating source being thermal dust. The spectral index and temperature for the thermal dust component is uniform throughout the sky. To demonstrate the leakage we provide two input cases, one with only the intensity map scanned (blue), and one with only the polarisation maps scanned (green), the latter showing an 'un-resolved' leakage component. The correction for the intensity to polarisation leakage component is demonstrated using two input cases, one with only the intensity map scanned and one with all intensity and polarisation maps scanned. In the former case, the residual leakage is down at the level of machine precision (orange), and the latter showing a residue at the level of the 'un-resolved' polarisation leakage signal since it was unaccounted for in the correction procedure. Figure 6. We demonstrate the correction algorithm for multiple input foreground sources (thermal dust and synchrotron) that scale uniformly over the sky and with no input noise in the timestream or on the template maps. In our input we provide only the intensity maps to keep our results unambiguous due to any 'un-resolved' polarisation leakage signal. The spectrum of the residue (orange) shows that the multiple contaminants were corrected for upto machine precision.

The leakage signal in figure 6 contains contribution from both thermal dust and synchrotron. The residual leakage is limited by the spurious cross-polar leakage we have seen in the first test, suggesting that our estimator is unbiased even when more than one leakage JCAP07(2019)043

Figure 7. We demonstrate the limit of the correction algorithm due to intrinsic noise in the timestream and the template maps, and also the limit due to spatially varying spectral parameters of the foreground component. This is performed at 140 GHz and with thermal dust as the only contaminant. We construct two input cases, one where the thermal dust has a variable spectral index (scaling) and greybody temperature, and one where these are uniform. The BB spectrum of the bandpass leakage for these cases are shown in orange and cyan respectively. Both cases have noise in the timestream and in the template maps. In the latter case, with uniform scaling of the thermal dust, the residual spectra (blue) varies significantly more than the former (red) and is due to the uncertainty on estimating the amplitude of the spurious signal. For the former case, higher order leakage terms not accounted for by the correction algorithm provide a floor for the residual spectra. component is estimated simultaneously, and there is no cross-talk between components that may arise due to possible degeneracies.

C Limit due to spatial variation of spectral parameters

Our third test case is designed to check whether the correction algorithm is noise limited or limited by the spatial variation of the spectral parameters. We construct two sets of input maps for the same detector pair as used in the previous tests. For simplicity we consider only the CMB and thermal dust component in the 140 GHz band. In addition to the uniformly scaled foreground emission maps we construct our second set of sky maps with a spatially varying spectral index and dust temperature, as we have done in section 5. The thermal dust template map used in both cases contains a 5 µK.arcmin white noise. For both sets of input maps we perform five sets of simulations and correction, each with a different white noise realisation of r.m.s 40 µK √ s in the timestream. In figure 7 we see that when our input and template was uniformly scaled, the residual has a very wide scatter. This is attributed to the uncertainty in estimating the amplitude of the leakage signal. However, when our input has spatially varying spectral parameters, the residue from the map alone is dominant and masks the effect of the uncertainty in measuring the leakage amplitude.

Figure 1 .

 1 Figure 1. Single detector, full-sky noise spectra in polarisation for four scan strategies and for a scan duration of 1 year are shown for an input noise rms of 40 µK √ s.The scan strategies are chosen to demonstrate a wide variety of precession angle (α) and spin angle (β) combinations, keeping α + β = 95 • . The scan strategies have a common spin period of 0.1 rpm, a precession period of 96.18 minutes and a sampling rate of 40 Hz. An ideal and uniform scan would have projected the noise on the theoretical line (solid red). The noise for the different configurations are however more than the theoretical value, increasing with a decrease in α. This demonstrates the boosting of the noise spectra as a result of inadequate angular variation and is exacerbated for lower values of precession angle α.

Figure 2 .

 2 Figure 2. Histogram of the noise variance of the Q (left column) and U (right column) polarization component maps, projected in the Ecliptic coordinate system, in the cases where the temperature leakage component is estimated as described by equation (3.9) (blue), and if we assume no bandpass mismatch error and noise described by equation (3.10) (red). The corresponding scan strategies for each row is labelled and each of them has a common spin period of 0.1 rpm, a precession period of 96.18 minutes and a sampling rate of 40 Hz. The solid vertical lines correspond to the mean of the histograms of the same colour and the dotted vertical lines their median.

Table 1 .

 1 The table shows the 16 detector used in our simulations, their position on the focal plane and their band properties. A positive value of position means it is further away from the spin axis.

	Pair Detector Pol angle Position	ν min	ν max	ν mid	ν min	ν max	ν mid
			Degrees Arc-mins		80 GHz		140 GHz
		1a	0	9.875	69.66 89.68 79.67 123.16 158.01 140.58
	1	1b	90	9.875	70.39 90.36 80.38 122.19 157.17 139.68
		2a	45	9.875	69.97 90.28 80.02 122.09 158.15 140.12
	2	2b	135	9.875	70.21 90.03 80.12 121.96 157.08 139.52
		3a	0	19.75	69.7 89.97 79.84 122.62 158.04 140.33
	3	3b	90	19.75	70.11 89.84 79.97 122.62 157.02 139.82
		4a	45	19.75	69.86 89.88 79.87 122.97 157.02 140.00
	4	4b	135	19.75	70.32 90.16 80.24 122.03 157.19 139.61
		5a	0	29.625	69.82 90.01 79.92 121.85 157.56 139.71
	5	5b	90	29.625	69.72 89.74 79.73 123.16 157.58 140.37
		6a	45	29.625	70.01 89.81 79.91 122.84 156.91 139.88
	6	6b	135	29.625	70.36 89.70 80.03 122.37 157.54 139.95
		7a	0	39.5	70.15 90.06 80.01 121.83 156.84 139.34
	7	7b	90	39.5	70.34 90.25 80.29 122.52 157.66 140.09
		8a	45	39.5	69.74 89.66 79.70 122.69 158.09 140.39
	8	8b	135	39.5	69.85 90.26 80.05 122.20 158.02 140.11

In principle, the other Stokes parameters, Q and U also depend on the bandpass of each detector, but we neglect this dependence as second order in our map-making model.

https://pla.esac.esa.int/pla/.

https://lambda.gsfc.nasa.gov/.

Acknowledgments RB thanks the funding from CNRS-IN2P3 and from the Research Council of Norway through grant 263011 during the development and writing of this paper, and also NERSC for the computing resources. MH is supported in part by JSPS KAKENHI No. JP15H05891. TM was partly supported by World Premier International Research Center Initiative (WPI Initiative), MEXT, Japan, and in part by JSPS KAKENHI No. 18KK0083.

JCAP07(2019)043

To recapitulate the relevant results from [START_REF] Hoang | Bandpass mismatch error for satellite CMB experiments I: estimating the spurious signal[END_REF], the spurious signal for each detector pair and for each foreground component is proportional to the foreground intensity signal, the bandpass mismatch parameter and the second order crossing moments cos 2ψ and sin 2ψ . The signal also scales in power inversely with the number of detector pairs. It is also important to mention that the shape of the leakage signal is determined by the foreground intensity signal and the scan strategy.

In figure 3, we illustrate the magnitude of the correction at the polarization map level by comparing the leakage map to the residual leakage left after the correction has been performed. On the full set of 16 detectors, we see a significant reduction in the leakage on to the polarization maps after correction. The residual leakage is dominated by the effect of the spatial variation of the emission laws of the foreground components. We demonstrate this explicitly in appendix C where we show the effect of noise, spatial variation of spectral parameters and multiple foregrounds. We see that under ideal conditions of no noise and no spatial variations, we are capable of achieving a correction at the level of machine precision even when we use multiple templates.

We subsequently quantify the leakage to the BB polarization power spectrum in figure 4. For this we simply estimate the power spectrum of the masked leakage and the residual leakage maps using Healpix anafast. In all cases, we see a significant improvement, of at least 2 orders of magnitude in power, to the bias in the BB spectrum. The residual leakage scales in a way that is similar to that of the leakage itself.

JCAP07(2019)043