Madam - a map-making method for CMB experiments Keihänen, E.; Kurki-Suonio, H.; Poutanen, T.
Monthly notices of the Royal Astronomical Society,
06/2005, Letnik:
360, Številka:
1
Journal Article
Recenzirano
Odprti dostop
We present a new map-making method for cosmic microwave background (CMB) measurements. The method is based on the destriping technique, but it also utilizes information about the noise spectrum. The ...low-frequency component of the instrument noise stream is modelled as a superposition of a set of simple base functions, whose amplitudes are determined by means of maximum-likelihood analysis, involving the covariance matrix of the amplitudes. We present simulation results with 1/f noise and show a reduction in the residual noise with respect to ordinary destriping. This study is related to Planck Low Frequency Instrument (LFI) activities.
BEYONDPLANCK Herman, D.; Hensley, B.; Andersen, K. J. ...
Astronomy and astrophysics (Berlin),
07/2023, Letnik:
675
Journal Article
Recenzirano
Odprti dostop
We constrained the level of polarized anomalous microwave emission (AME) on large angular scales using
Planck
Low-Frequency Instrument (LFI) and WMAP polarization data within a Bayesian cosmic ...microwave background (CMB) analysis framework. We modeled synchrotron emission with a power-law spectral energy distribution, as well as the sum of AME and thermal dust emission through linear regression with the
Planck
High-Frequency Instrument (HFI) 353 GHz data. This template-based dust emission model allowed us to constrain the level of polarized AME while making minimal assumptions on its frequency dependence. We neglected CMB fluctuations, but show through simulations that these fluctuations have a minor impact on the results. We find that the resulting AME polarization fraction confidence limit is sensitive to the polarized synchrotron spectral index prior. In addition, for prior means
β
s
< −3.1 we find an upper limit of
p
AME
max
≲ 0.6% (95% confidence). In contrast, for means
β
s
= −3.0, we find a nominal detection of
p
AME
= 2.5 ± 1.0% (95% confidence). These data are thus not strong enough to simultaneously and robustly constrain both polarized synchrotron emission and AME, and our main result is therefore a constraint on the AME polarization fraction explicitly as a function of
β
s
. Combining the current
Planck
and WMAP observations with measurements from high-sensitivity low-frequency experiments such as C-BASS and QUIJOTE will be critical to improve these limits further.
We study destriping as a map-making method for temperature-and-polarization data for cosmic microwave background observations. We present a particular implementation of destriping and study the ...residual error in output maps, using simulated data corresponding to the 70 GHz channel of the Planck satellite, but assuming idealized detector and beam properties. The relevant residual map is the difference between the output map and a binned map obtained from the signal + white noise part of the data stream. For destriping it can be divided into six components: unmodeled correlated noise, white noise reference baselines, reference baselines of the pixelization noise from the signal, and baseline errors from correlated noise, white noise, and signal. These six components contribute differently to the different angular scales in the maps. We derive analytical results for the first three components. This study is related to Planck LFI activities.
BEYONDPLANCK Herman, D; Hensley, B; Andersen, K J ...
Astronomy and astrophysics (Berlin),
07/2023, Letnik:
675
Journal Article
Recenzirano
Odprti dostop
We constrained the level of polarized anomalous microwave emission (AME) on large angular scales using Planck Low-Frequency Instrument (LFI) and WMAP polarization data within a Bayesian cosmic ...microwave background (CMB) analysis framework. We modeled synchrotron emission with a power-law spectral energy distribution, as well as the sum of AME and thermal dust emission through linear regression with the Planck High-Frequency Instrument (HFI) 353 GHz data. This template-based dust emission model allowed us to constrain the level of polarized AME while making minimal assumptions on its frequency dependence. We neglected CMB fluctuations, but show through simulations that these fluctuations have a minor impact on the results. We find that the resulting AME polarization fraction confidence limit is sensitive to the polarized synchrotron spectral index prior. In addition, for prior means βs < −3.1 we find an upper limit of pAMEmax ≲ 0.6% (95% confidence). In contrast, for means βs = −3.0, we find a nominal detection of pAME = 2.5 ± 1.0% (95% confidence). These data are thus not strong enough to simultaneously and robustly constrain both polarized synchrotron emission and AME, and our main result is therefore a constraint on the AME polarization fraction explicitly as a function of βs. Combining the current Planck and WMAP observations with measurements from high-sensitivity low-frequency experiments such as C-BASS and QUIJOTE will be critical to improve these limits further.
Madam- a map-making method for CMB experiments Keihänen, E.; Kurki-Suonio, H.; Poutanen, T.
Monthly notices of the Royal Astronomical Society,
06/2005, Letnik:
360, Številka:
1
Journal Article
We present a system-level description of the Low Frequency Instrument (LFI) considered as a differencing polarimeter, and evaluate its expected performance. The LFI is one of the two instruments on ...board the ESA Planck mission to study the cosmic microwave background. It consists of a set of 22 radiometers sensitive to linear polarisation, arranged in orthogonally-oriented pairs connected to 11 feed horns operating at 30, 44 and 70 GHz. In our analysis, the generic Jones and Mueller-matrix formulations for polarimetry are adapted to the special case of the LFI. Laboratory measurements of flight components are combined with optical simulations of the telescope to investigate the values and uncertainties in the system parameters affecting polarisation response. Methods of correcting residual systematic errors are also briefly discussed. The LFI has beam-integrated polarisation efficiency >99% for all detectors, with uncertainties below 0.1%. Indirect assessment of polarisation position angles suggests that uncertainties are generally less than 0$\fdg$5, and this will be checked in flight using observations of the Crab nebula. Leakage of total intensity into the polarisation signal is generally well below the thermal noise level except for bright Galactic emission, where the dominant effect is likely to be spectral-dependent terms due to bandpass mismatch between the two detectors behind each feed, contributing typically 1–3% leakage of foreground total intensity. Comparable leakage from compact features occurs due to beam mismatch, but this averages to < 5 × 10-4 for large-scale emission. An inevitable feature of the LFI design is that the two components of the linear polarisation are recovered from elliptical beams which differ substantially in orientation. This distorts the recovered polarisation and its angular power spectrum, and several methods are being developed to correct the effect, both in the power spectrum and in the sky maps. The LFI will return a high-quality measurement of the CMB polarisation, limited mainly by thermal noise. To meet our aspiration of measuring polarisation at the 1% level, further analysis of flight and ground data is required. We are still researching the most effective techniques for correcting subtle artefacts in polarisation; in particular the correction of bandpass mismatch effects is a formidable challenge, as it requires multi-band analysis to estimate the spectral indices that control the leakage.
ABSTRACT
Extraction of the cosmic microwave background (CMB) angular power spectrum is a challenging task for current and future CMB experiments due to the large data sets involved. Here we describe ...an implementation of Monte Carlo apodized spherical transform estimator (MASTER) described in Hivon et al., which exploits the destriping technique as a map‐making method. In this method a noise estimate based on destriped noise‐only Monte Carlo (MC) simulations is subtracted from the pseudo‐angular power spectrum. As a working case we use realistic simulations of the Planck low‐frequency instrument (LFI). We found that the effect of destriping on a pure sky signal is minimal and requires no correction. Instead we found an effect related to the distribution of detector pointings, which affects the high‐ℓ part of the power spectrum. We correct for this by subtracting a ‘signal bias’ estimated by MC simulations. We also give analytical estimates for this signal bias. Our method is fast and accurate enough (the estimator is unbiased and errors are close to theoretical expectations for maximal accuracy) to estimate the CMB angular power spectra for current and future CMB space missions. This study is related to Planck LFI activities.
BEYONDPLANCK Keihänen, E; A.-S. Suur-Uski; Andersen, K J ...
Astronomy and astrophysics (Berlin),
07/2023, Letnik:
675
Journal Article
Recenzirano
Odprti dostop
We present a Gibbs sampling solution to the mapmaking problem for cosmic microwave background (CMB) measurements that builds on existing destriping methodology. Gibbs sampling breaks the ...computationally heavy destriping problem into two separate steps: noise filtering and map binning. Considered as two separate steps, both are computationally much cheaper than solving the combined problem. This provides a huge performance benefit as compared to traditional methods and it allows us, for the first time, to bring the destriping baseline length to a single sample. Here, we applied the Gibbs procedure to simulated Planck 30 GHz data. We find that gaps in the time-ordered data are handled efficiently by filling them in with simulated noise as part of the Gibbs process. The Gibbs procedure yields a chain of map samples, from which we are able to compute the posterior mean as a best-estimate map. The variation in the chain provides information on the correlated residual noise, without the need to construct a full noise covariance matrix. However, if only a single maximum-likelihood frequency map estimate is required, we find that traditional conjugate gradient solvers converge much faster than a Gibbs sampler in terms of the total number of iterations. The conceptual advantages of the Gibbs sampling approach lies in statistically well-defined error propagation and systematic error correction. This methodology thus forms the conceptual basis for the mapmaking algorithm employed in the BEYONDPLANCK framework, which implements the first end-to-end Bayesian analysis pipeline for CMB observations.
BEYONDPLANCK Gjerløw, E; Ihle, H T; Galeotta, S ...
Astronomy and astrophysics (Berlin),
07/2023, Letnik:
675
Journal Article
Recenzirano
Odprti dostop
We present a Bayesian calibration algorithm for cosmic microwave background (CMB) observations as implemented within the global end-to-end BEYONDPLANCK framework and applied to the Planck Low ...Frequency Instrument (LFI) data. Following the most recent Planck analysis, we decomposed the full time-dependent gain into a sum of three nearly orthogonal components: one absolute calibration term, common to all detectors, one time-independent term that can vary between detectors, and one time-dependent component that was allowed to vary between one-hour pointing periods. Each term was then sampled conditionally on all other parameters in the global signal model through Gibbs sampling. The absolute calibration is sampled using only the orbital dipole as a reference source, while the two relative gain components were sampled using the full sky signal, including the orbital and Solar CMB dipoles, CMB fluctuations, and foreground contributions. We discuss various aspects of the data that influence gain estimation, including the dipole-polarization quadrupole degeneracy and processing masks. Comparing our solution to previous pipelines, we find good agreement in general, with relative deviations of −0.67% (−0.84%) for 30 GHz, 0.12% (−0.04%) for 44 GHz and −0.03% (−0.64%) for 70 GHz, compared to Planck PR4 and Planck 2018, respectively. We note that the BEYONDPLANCK calibration was performed globally, which results in better inter-frequency consistency than previous estimates. Additionally, WMAP observations were used actively in the BEYONDPLANCK analysis, which both breaks internal degeneracies in the Planck data set and results in an overall better agreement with WMAP. Finally, we used a Wiener filtering approach to smoothing the gain estimates. We show that this method avoids artifacts in the correlated noise maps as a result of oversmoothing the gain solution, which is difficult to avoid with methods like boxcar smoothing, as Wiener filtering by construction maintains a balance between data fidelity and prior knowledge. Although our presentation and algorithm are currently oriented toward LFI processing, the general procedure is fully generalizable to other experiments, as long as the Solar dipole signal is available to be used for calibration.
BEYONDPLANCK Herman, D; Watson, R A; Andersen, K J ...
Astronomy and astrophysics (Berlin),
07/2023, Letnik:
675
Journal Article
Recenzirano
Odprti dostop
We describe the correction procedure for Analog-to-Digital Converter (ADC) differential non-linearities (DNL) adopted in the Bayesian end-to-end BEYONDPLANCK analysis framework. This method is nearly ...identical to that developed for the official Planck Low Frequency Instrument (LFI) Data Processing Center (DPC) analysis, and relies on the binned rms noise profile of each detector data stream. However, rather than building the correction profile directly from the raw rms profile, we first fit a Gaussian to each significant ADC-induced rms decrement, and then derive the corresponding correction model from this smooth model. The main advantage of this approach is that only samples which are significantly affected by ADC DNLs are corrected, as opposed to the DPC approach in which the correction is applied to all samples, filtering out signals not associated with ADC DNLs. The new corrections are only applied to data for which there is a clear detection of the non-linearities, and for which they perform at least comparably with the DPC corrections. Out of a total of 88 LFI data streams (sky and reference load for each of the 44 detectors) we apply the new minimal ADC corrections in 25 cases, and maintain the DPC corrections in 8 cases. All these corrections are applied to 44 or 70 GHz channels, while, as in previous analyses, none of the 30 GHz ADCs show significant evidence of non-linearity. By comparing the BEYONDPLANCK and DPC ADC correction methods, we estimate that the residual ADC uncertainty is about two orders of magnitude below the total noise of both the 44 and 70 GHz channels, and their impact on current cosmological parameter estimation is small. However, we also show that non-idealities in the ADC corrections can generate sharp stripes in the final frequency maps, and these could be important for future joint analyses with the Planck High Frequency Instrument (HFI), Wilkinson Microwave Anisotropy Probe (WMAP), or other datasets. We therefore conclude that, although the existing corrections are adequate for LFI-based cosmological parameter analysis, further work on LFI ADC corrections is still warranted.