Difference between revisions of "LFI systematic effect uncertainties"
|Line 1:||Line 1:|
== Overview ==
== Overview ==
|Line 103:||Line 104:|
Revision as of 21:05, 17 March 2013
- 1 Overview
- 2 Summary of uncertainties due to systematic effects
- 3 Effects independent of sky signal
- 4 Effects dependent on sky signal
- 5 References
Known systematic effects in the Planck-LFI data can be divided into two broad categories: effects independent of the sky signal, which can be considered as additive or multiplicative spurious contributions to the measured timelines, and effects which are dependent on the sky and that cannot be considered independently from the observation strategy.
Here we report a brief summary of these effects, all the details can be found in #planck2013-p02a Planck-2013-III .
Summary of uncertainties due to systematic effects
In this section we provide a top-level overview of the uncertainties due to systematic effects in the Planck-LFI CMB temperature maps and power spectra. Table 1 provides a list of these effects with short indications of their cause, strategies for removal and references to sections and/or papers where more information is found.
Table 1. List of known instrumental systematic effects in Planck LFI.
The impact of 1/f noise has been assessed using half-ring noise maps normalised to the white noise estimate at each pixel obtained from the white noise covariance matrix, so that a perfectly white noise map would be Gaussian and isotropic with unit variance. Deviations from unity trace the contribution of residual 1/f noise in the final maps, which ranges from 0.06% at 70 GHz to 2% at 30 GHz. Pixel uncertainties due to other systematic effects have been calculated on simulated maps degraded at Nside = 128 at 30 and 44 GHz and Nside = 256 at 70 GHz in order to approximate the optical beam size. This downgrading has been applied in all cases a systematic effect has been evaluated at map level.
In Table 2 we list the r.m.s. and the difference between the 99% and the 1% quantities in the pixel value distributions. For simplicity we refer to this difference as peak-to-peak (p-p) difference although it neglects outliers but effectively approximates the peak-to-peak variation of the effect on the map. The last row in Table 2 reports the r.m.s. value of survey difference maps (first − second survey, after a correction which takes into account the white noise. These values thus represent the level of r.m.s. uncertainty per pixel due to known and unknown systematic effects.
Table 2. Summary of systematic effects uncertainties on maps in μKCMB.
Angular power spectra have been obtained from full resolution (Nside = 1024) systematic effect maps at each frequency using HEALPix Anafast #gorski2005. We have then evaluated the propagation of the various effects in the final CMB map by assuming a simple ILC component separation. In Fig. 1 we show how the power spectrum of the various effects compared with the Planck temperature spectrum, with the noise level coming from the half-ring difference maps and with the residual map obtained from a difference map between survey 1 and survey 2. The large plot in the top panel shows the power spectra obtained from frequency-independent maps resulting from the weighted-average of frequency maps. Spectra in the three small plots in the lower panel, instead, show contributions of systematic effects from individual frequency maps.
Our assessment shows that the global impact of systematic effects uncertainties is at least two order of magnitudes less than the CMB power spectrum and demonstrates the robustness of Planck-LFI temperature anisotropy measurements. Comparison between the total simulated systematic effects and residual signal obtained by differencing survey 1 and survey 2 maps highlights an excess signal in the multipole range that is not completely accounted for in our simulations. This excess comes mainly from the 30 GHz channel and is likely to be caused by galactic emissions picked up by beam sidelobes. Understanding this excess and further reducing the level of residual systematic uncertainties is the primary goal of our current analysis to obtain polarization measurements with a level of purity comparable to what has been achieved with temperature anisotropies.
Figure 1. The angular power spectrum of the various systematic effects compared to the Planck temperature spectrum. The black curve, representing the total contribution, has been derived from a map where all the systematic effects have been summed. Top panel: power spectra obtained from frequency independent maps resulting from the weighted-average of individual systematic effects frequency maps. Bottom panel: contributions of systematic effects from individual frequency maps. The CMB curve corresponds to the Planck best-fit model presented in #planck2013-p08 Planck-2013-XV .
Effects independent of sky signal
Noise correlations and 1/f noise
As described in #seiffert2002, planck2011-1-4 imperfect matching of components generates isolation between the complementary diodes of a receiver between −10 and −15 dB. This imperfect isolation leads to a small anti correlated component in the white noise that is cancelled by a weighted average of the time ordered data from the two diodes of each receiver as the first step of analysis. This avoids the complication of tracking the anti correlated white noise throughout the analysis. We treat the combined diode data as the raw data, and calibration, noise estimation, map making etc. are performed on these combined data. The weights were determined from some initial estimates of the calibrated noise for each detector, and are kept fixed for the entire mission.
We estimate the signal subtracted noise power spectrum of each receiver on 5-day time periods. Except for specific, mostly well understood events, shorter time scale noise estimation does not produce any evident trends. For nearly all the radiometers our noise model is a very good approximation of the power spectrum.
Over the course of the nominal mission, the noise is well fit by the model, with the exception of the early parts of sky survey 3. During this time, thermal instabilities brought on by the switchover from the nominal to the redundant sorption cooler cause poor fits and some changes in the parameters.
The LFI is susceptible to temperature fluctuations in the 300 K back-end modules, in the 4 K reference loads and in the 20 K focal plane.
Temperature of 70 GHz reference loads is actively con- trolled by a proportional-integral-derivative (PID) system and is very stable (δTrms ∼ 0.13 mK). Reference loads of 30 and 44 GHz channels, instead, do not benefit from active thermal control. Their temperature is consequently more unstable and susceptible to major system-level events like, for example, the switchover to the redundant sorption cooler.
The 20 K LFI focal plane temperature is measured by a sensor placed on the feed horn flange of the LFI28 receiver. The temperature during the first sky survey was very stable. Towards the end of the first year of operations the sorption cooler performance started to degrade and its stability was maintained with a series of controlled temperature changes. The switchover to the redundant cooler left a clear signature on all the main LFI temperatures. After this operation the level of temperature fluctuations on the focal plane increased unexpectedly, and this was later understood to be the effect of liquid hydrogen that was still present in the cold-end of the nominal cooler because the degraded compressor system was not able to absorb all the hydrogen that was present in the cooler line. Although this effect was later mitigated by a series of dedicated operations, most of the third sky survey suffered from a higher-than-nominal level of temperature variations.
The temperature of the 300 K electronics box is measured by one of its temperature sensors. During the first sky survey the back-end temperature suffered from a 24-hour fluctuation caused by the satellite transponder that was switched on daily during contact with the ground station. After day 258 the system was left continuously on and the 24-hour modulation disappeared. This operation caused an increase of the absolute temperature level. The second temperature change occurred in correspondence to the sorption cooler switch over operation. A yearly temperature modulation due to the satellite rotation around the sun and a temperature spike need also to be considered. This was caused by an operational anomaly that caused the satellite to fail repointing for an entire day with a corresponding temperature increase of the warm units.
Details about the thermal stability performance of Planck can be found in #planck2011-1-3, while the susceptibility of the LFI to temperature variations is discussed in #terenzi2009b.
The signal detected by the radiometers can vary because of fluctuations in the front-end and back-end amplifiers bias voltages. In the LFI these fluctuations occurred according to two time scales:
- slow electric drifts, due to thermal changes in the power supply, in the RF amplifiers and in the detector diodes;
- fast and sudden electric instabilities, arising in the warm electronics or from electromagnetic interference effects, and affecting both the cold amplifiers and the warm detector diodes.
The effect of slow drifts is suppressed by the pseudo-correlation architecture of the differential radiometers. Fast electric changes produce quasi-random fluctuations and abrupt steep drops or jumps in the signal. If jumps are caused by instabilities in the front-end bias voltage then the effect involves the output voltage of both diodes in the radiometer. When the jumps occur in the back-end detector diodes (so-called “pop-corn noise”) they impact only the output voltage of the corresponding diode and affect sky and reference load samples. In both cases the differenced signal is largely immune from these effects.
This effect is caused by pickup from the housekeeping electronics clock that occurs after the detector diodes and before the ADC converters #meinhold2009, mennella2010, mennella2011. This spurious signal is detected in the radiometer time-domain outputs as a one second square wave with a rising edge near 0.5 s and a falling edge near 0.75 s in on-board time. In the frequency domain it appears at multiples of 1 Hz. Frequency spikes are present at some level in the output from all detectors, but affect the 44 GHz data most strongly because of the low voltage output and high post-detection gain values in that channel. For this reason spikes are removed from the 44 GHz time-ordered data via template fitting, as described in #planck2013-p02 Planck-2013-II .
Effects dependent on sky signal
Straylight contamination arise from the spurious signal pickup from the telescope far sidelobes. Main sources of straylight contamination are the Galaxy, especially at 30 GHz, and the cosmological dipole, mainly detected in the directions of the main and sub reflectors spillover. In principle we should also include the straylight contribution from the orbital dipole, but its effect is of a factor ten lower than the cosmic dipole so that it can safely be neglected in this framework (but it has been considered in the calibration pipeline). Further details about the Planck optical system are reported in #tauber2010 and the LFI beam properties are provided in #sandri2010.
Straylight impacts the measured signal essentially in two ways: (i) through direct contamination and coupling with the main beam sky signal and (ii) in the photometric calibration of the radiometer detected signal. The impact on calibration and the adopted mitigation strategies are described in #planck2013-p02b Planck-2013-V . Because of the beam orientation, the straylight fingerprint is different in odd surveys compared to even surveys. The Galaxy, for example, is detected by the sub spillover in the odd surveys and by the main spillover in the even surveys. Because the sub spillover points approximately in the main beam direction, the Galaxy straylight pattern is close to the galactic plane. The main spillover, instead, points at about 85 degrees from the main beam so that the Galaxy is re-imaged onto a ring.
ADC non linearity
The linearity of analogue-to-digital converters (ADC) requires that the voltage step sizes between successive binary outputs are constant over the entire input dynamic range. If these steps are not constant we have a non linearity in the ADC response that leads to calibration errors.
The typical fingerprint of ADC non linearity is a variation of the detector voltage output white noise not paired by a detectable variation in the voltage level. This effect was observed in the LFI radiometer data for the first time in flight, where drops of few percent were observed in the voltage white noise but not in the output level over periods of few weeks. The typical amplitude of the region where the non linearity occurs is of the order of 1 mV, corresponding to about three bits in the ADC. The ADC effect is strongest (3 to 6%) in the 44 GHz channels, because of their lower detector voltages. The ADC non linearity effect has been characterised from flight data and removed from the TOIs according to the procedure described in the main LFI data processing paper #planck2013-p02 Planck-2013-II .
Imperfect photometric calibration
An important set of systematic effects are those related with the photometric calibration of the radiometers. Such effects are discussed at length by #planck2013-p05 Planck-2013-XXVIII . There are three different kinds of systematic effects that can affect the calibration.
- Incorrect assumptions regarding the calibration signal.In the case of LFI, the signal used for the calibration is the dipolar field caused by the motion of the Solar System with respect to the CMB rest frame and by the motion of the spacecraft around the Sun. We model the former using the values quoted by #hinshaw2009 and the latter using the spacecraft’s attitude information. Any error in the numbers would directly lead to an error in the calibration of Planck-LFI data.
- Incorrect treatment of the calibration signal. To actually use any previous knowledge of the CMB dipole, we need to convolve the signal with the beam response of LFI radiometer. Any error in this step would produce a systematic effect in the map not only because of the wrong shape expected for the calibration signal, but also because of the removal of the (wrong) dipole from the calibrated maps done by the Planck-LFI pipeline #planck2013-p02 Planck-2013-II . Possible types of errors include: systematic errors in WMAP’s dipole estimate, wrong convolution of the expected dipole with the radiometer beams, uncorrect masking of the Galaxy when fitting the observed signal with the dipole, etc.
- Incorrect reconstruction of gain fluctuations. Some of the algorithms we used in calibrating LFI data for this release use the radiometer equation and the recorded variations of the radiometers total-power output to track gain changes. In principle, any deviation in the behaviour of the radiometer from the ideal case (e.g., ADC non linearities) can therefore induce systematic effects in the gain curves.
Pointing uncertainties are translated into uncertainties in pixel temperature measurements. If pointing uncertainties are not constant in time then the statistics of the sky anisotropy measurements is not preserved with a consequent impact on power spectrum and cosmological parameters. In Planck-LFI pointing uncertainties arise from:
- Satellite pointing determination. The Planck Attitude Control Movement System guarantees a pointing accuracy of ∼ 2′′ #planck2013-p01 Planck-2013-I , which is well within scientific requirements. However, small non idealities in the system and errors in the attitude reconstruction caused, for example, by thermoelastic effects, can affect the data.
- Uncertainties in the focal plane geometry reconstruction. The measurement of Planck-LFI focal plane geometry is based on the determination of the beam pointing with respect to the nominal line of sight exploiting Jupiter observation. The peak of each beam has been determined by fitting data with a bivariate gaussian function which may be not representative of the real beam centre.
(Planck) Low Frequency Instrument
Cosmic Microwave background
(Hierarchical Equal Area isoLatitude Pixelation of a sphere, <ref name="Template:Gorski2005">HEALPix: A Framework for High-Resolution Discretization and Fast Analysis of Data Distributed on the Sphere, K. M. Górski, E. Hivon, A. J. Banday, B. D. Wandelt, F. K. Hansen, M. Reinecke, M. Bartelmann, ApJ, 622, 759-771, (2005).
[LFI meaning]: absolute calibration refers to the 0th order calibration for each channel, 1 single number, while the relative calibration refers to the component of the calibration that varies pointing period by pointing period.
JFET, Bellow, PAU et REU
JFET, Bellow, PAU et REU
analog to digital converter