-
Cosmological measurements from the CMB and BAO are insensitive to the tail probability in the assumed likelihood
Authors:
Jordan Krywonos,
Simone Paradiso,
Alex Krolewski,
Shahab Joudaki,
Will Percival
Abstract:
When fitting cosmological models to data, a Bayesian framework is commonly used, requiring assumptions on the form of the likelihood and model prior. In light of current tensions between different data, it is interesting to investigate the robustness of cosmological measurements to statistical assumptions about the likelihood distribution from which the data was drawn. We consider the impact of ch…
▽ More
When fitting cosmological models to data, a Bayesian framework is commonly used, requiring assumptions on the form of the likelihood and model prior. In light of current tensions between different data, it is interesting to investigate the robustness of cosmological measurements to statistical assumptions about the likelihood distribution from which the data was drawn. We consider the impact of changes to the likelihood caused by uncertainties due to the finite number of mock catalogs used to estimate the covariance matrix, leading to the replacement of the standard Gaussian likelihood with a multivariate $t$-distribution. These changes to the likelihood have a negligible impact on recent cosmic microwave background (CMB) lensing and baryon acoustic oscillation (BAO) measurements. We then extend our analysis to perform a sensitivity test on the Gaussian likelihoods typically adopted, considering how increasing the size of the tails of the likelihood (again using a $t$-distribution) affects cosmological inferences. For an open $Λ$CDM model constrained by BAO alone, we find that increasing the weight in the tails shifts and broadens the resulting posterior on the parameters, with a $\sim$0.2-0.4$σ$ effect on $Ω_Λ$ and $Ω_{\rm k}$. In contrast, the CMB temperature and polarization constraints in $Λ$CDM showed less than 0.03$σ$ changes in the parameters, except for $\{τ$, ln($10^{10}A_{\rm s})$, $σ_8$, $S_8$, $σ_8Ω_{\rm m}^{0.25}$, $z_{\rm re}$, $10^9A_{\rm s}e^{-2τ}\}$ which shifted by around 0.1-0.2$σ$. If we use solely $\ell < 30$ data, the amplitude $A_{\rm s} e^{-2τ}$ varies in the posterior mean by 0.7$σ$ and the error bars increase by 6%. We conclude, at least for current-generation CMB and BAO measurements, that uncertainties in the shape and tails of the likelihood do not contribute to current tensions.
△ Less
Submitted 6 March, 2024;
originally announced March 2024.
-
Evaluating extensions to LCDM: an application of Bayesian model averaging and selection
Authors:
S. Paradiso,
G. McGee,
W. J. Percival
Abstract:
We employ Bayesian Model Averaging (BMA) as a powerful statistical framework to address key cosmological questions about the universe's fundamental properties. We explore extensions beyond the standard $Λ$CDM model, considering a varying curvature density parameter $Ω_{\rm k}$, a spectral index $\mathrm{n}_{\rm s}=1$ and a varying $n_{\rm run}$, a constant dark energy equation of state (EOS)…
▽ More
We employ Bayesian Model Averaging (BMA) as a powerful statistical framework to address key cosmological questions about the universe's fundamental properties. We explore extensions beyond the standard $Λ$CDM model, considering a varying curvature density parameter $Ω_{\rm k}$, a spectral index $\mathrm{n}_{\rm s}=1$ and a varying $n_{\rm run}$, a constant dark energy equation of state (EOS) $w_0$CDM and a time-dependent one $w_0w_a$CDM. We also test cosmological data against a varying effective number of neutrino species $N_{\rm eff}$. Data from different combinations of cosmic microwave background (CMB) data from the last Planck PR4 analysis, CMB lensing from Planck 2018, baryonic acoustic oscillations (BAO) and the Bicep-KECK 2018 results, are used. We find that the standard $Λ$CDM model is favoured when combining CMB data with CMB lensing, BAO and Bicep-KECK 2018 data against $K-Λ$CDM model $N_{\rm eff}-Λ$CDM with a probability $> 80\%$. When investigating the dark energy EOS, we find that this dataset is not able to express a strong preference between the standard $Λ$CDM model and the constant dark energy EOS model $w_0$CDM, with an approximately split model posterior probability of $\approx 60\%:40\%$ in favour of $Λ$CDM, whereas the time-varying dark energy EOS model is ruled out. Finally, we find that the CMB data alone show a strong preference for a model that includes the running of the spectral index $n_{\rm run}$, with a probability $\approx 90\%$, when compared to the $n_{\rm s}=1$ model and the standard $Λ$CDM. Overall, we find that including the model uncertainty in the considered cases does not significantly impact the Hubble tension.
△ Less
Submitted 5 July, 2024; v1 submitted 4 March, 2024;
originally announced March 2024.
-
LSPE-STRIP on-sky calibration strategy using bright celestial sources
Authors:
R. T. Génova-Santos,
M. Bersanelli,
C. Franceschet,
M. Gervasi,
C. López-Caraballo,
L. Mandelli,
M. Maris,
A. Mennella,
J. A. Rubiño-Martín,
F. Villa,
M. Zannoni,
C. Baccigalupi,
B. Caccianiga,
L. Colombo,
F. Cuttaia,
F. Farsian,
G. Morgante,
S. Paradiso,
G. Polenta,
S. Ricciardi,
M. Sandri,
A. Taylor,
L. Terenzi,
M. Tomasi
Abstract:
In this paper we describe the global on-sky calibration strategy of the LSPE-Strip instrument. Strip is a microwave telescope operating in the Q- and W-bands (central frequencies of 43 and 95 GHz respectively) from the Observatorio del Teide in Tenerife, with the goal to observe and characterise the polarised Galactic foreground emission, and complement the observations of the polarisation of the…
▽ More
In this paper we describe the global on-sky calibration strategy of the LSPE-Strip instrument. Strip is a microwave telescope operating in the Q- and W-bands (central frequencies of 43 and 95 GHz respectively) from the Observatorio del Teide in Tenerife, with the goal to observe and characterise the polarised Galactic foreground emission, and complement the observations of the polarisation of the cosmic microwave background to be performed by the LSPE-SWIPE instrument and other similar experiments operating at higher frequencies to target the detection of the B-mode signal from the inflationary epoch of the Universe. Starting from basic assumptions on some of the instrument parameters (NET, 1/f noise knee frequency, beam properties, observing efficiency) we perform realistic simulations to study the level of accuracy that can be achieved through observations of bright celestial calibrators in the Strip footprint (sky fraction of 30 %) on the determination and characterisation of the main instrument parameters: global and relative gain factors (in intensity and in polarisation), polarisation direction, polarisation efficiency, leakage from intensity to polarisation, beams, window functions and pointing model.
△ Less
Submitted 16 January, 2024; v1 submitted 8 January, 2024;
originally announced January 2024.
-
Measuring the CMB primordial B-modes with Bolometric Interferometry
Authors:
A. Mennella,
P. Ade,
A. Almela,
G. Amico,
L. H. Arnaldi,
J. Aumont,
S. Banfi,
E. S. Battistelli,
B. Bélier,
L. Bergé,
J. -Ph. Bernard,
P. de Bernardis,
M. Bersanelli,
J. Bonaparte,
J. D. Bonilla,
E. Bunn,
D. Buzi,
F. Cacciotti,
D. Camilieri,
F. Cavaliere,
P. Chanial,
C. Chapron,
L. Colombo,
F. Columbro,
A. Coppolecchia
, et al. (89 additional authors not shown)
Abstract:
The Q&U Bolometric Interferometer for Cosmology (QUBIC) is the first bolometric interferometer designed to measure the primordial B-mode polarization of the Cosmic Microwave Background (CMB). Bolometric interferometry is a novel technique that combines the sensitivity of bolometric detectors with the control of systematic effects that is typical of interferometry, both key features in the quest fo…
▽ More
The Q&U Bolometric Interferometer for Cosmology (QUBIC) is the first bolometric interferometer designed to measure the primordial B-mode polarization of the Cosmic Microwave Background (CMB). Bolometric interferometry is a novel technique that combines the sensitivity of bolometric detectors with the control of systematic effects that is typical of interferometry, both key features in the quest for the faint signal of the primordial B-modes. A unique feature is the so-called "spectral imaging", i.e., the ability to recover the sky signal in several sub-bands within the physical band during data analysis. This feature provides an in-band spectral resolution of Δν/ν \sim 0.04 that is unattainable by a traditional imager. This is a key tool for controlling the Galactic foregrounds contamination. In this paper, we describe the principles of bolometric interferometry, the current status of the QUBIC experiment and future prospects.
△ Less
Submitted 5 November, 2023;
originally announced November 2023.
-
The advantage of Bolometric Interferometry for controlling Galactic foreground contamination in CMB primordial B-modes measurements
Authors:
E. Manzan,
M. Regnier,
J-Ch. Hamilton,
A. Mennella,
J. Errard,
L. Zapelli,
S. A. Torchinsky,
S. Paradiso,
E. Battistelli,
M. Bersanelli,
P. De Bernardis,
M. De Petris,
G. D'Alessandro,
M. Gervasi,
S. Masi,
M. Piat,
E. Rasztocky,
G. E Romero,
C. G. Scoccola,
M. Zannoni,
the QUBIC Collaboration
Abstract:
In the quest for the faint primordial B-mode polarization of the Cosmic Microwave Background, three are the key requirements for any present or future experiment: an utmost sensitivity, excellent control over instrumental systematic effects and over Galactic foreground contamination. Bolometric Interferometry (BI) is a novel technique that matches them all by combining the sensitivity of bolometri…
▽ More
In the quest for the faint primordial B-mode polarization of the Cosmic Microwave Background, three are the key requirements for any present or future experiment: an utmost sensitivity, excellent control over instrumental systematic effects and over Galactic foreground contamination. Bolometric Interferometry (BI) is a novel technique that matches them all by combining the sensitivity of bolometric detectors, the control of instrumental systematics from interferometry and a software-based, tunable, in-band spectral resolution due to its ability to perform band-splitting during data analysis (spectral imaging). In this paper, we investigate how the spectral imaging capability of BI can help in detecting residual contamination in case an over-simplified model of foreground emission is assumed in the analysis. To mimic this situation, we focus on the next generation of ground-based CMB experiment, CMB-S4, and compare its anticipated sensitivities, frequency and sky coverage with a hypothetical version of the same experiment based on BI, CMB-S4/BI, assuming that line-of-sight (LOS) frequency decorrelation is present in dust emission but is not accounted for during component separation. We show results from a Monte-Carlo analysis based on a parametric component separation method (FGBuster), highlighting how BI has the potential to diagnose the presence of foreground residuals in estimates of the tensor-to-scalar ratio $r$ in the case of unaccounted Galactic dust LOS frequency decorrelation.
△ Less
Submitted 3 November, 2023;
originally announced November 2023.
-
A convenient approach to characterizing model uncertainty with application to early dark energy solutions of the Hubble tension
Authors:
S. Paradiso,
M. DiMarco,
M. Chen,
G. McGee,
W. J. Percival
Abstract:
Despite increasingly precise observations and sophisticated theoretical models, the discrepancy between measurements of H0 from the cosmic microwave background or from Baryon Acoustic Oscillations combined with Big-Bang Nucleosynthesis versus those from local distance ladder probes -- commonly known as the $H_0$ tension -- continues to perplex the scientific community. To address this tension, Ear…
▽ More
Despite increasingly precise observations and sophisticated theoretical models, the discrepancy between measurements of H0 from the cosmic microwave background or from Baryon Acoustic Oscillations combined with Big-Bang Nucleosynthesis versus those from local distance ladder probes -- commonly known as the $H_0$ tension -- continues to perplex the scientific community. To address this tension, Early Dark Energy (EDE) models have been proposed as alternatives to $Λ$CDM, as they can change the observed sound horizon and the inferred Hubble constant from measurements based on this. In this paper, we investigate the use of Bayesian Model Averaging (BMA) to evaluate EDE as a solution to the H0 tension. BMA consists of assigning a prior to the model and deriving a posterior as for any other unknown parameter in a Bayesian analysis. BMA can be computationally challenging in that one must approximate the joint posterior of both model and parameters. Here we present a computational strategy for BMA that exploits existing MCMC software and combines model-specific posteriors post-hoc. In application to a comprehensive analysis of cosmological datasets, we quantify the impact of EDE on the H0 discrepancy. We find an EDE model probability of $\sim$90% whenever we include the H0 measurement from Type Ia Supernovae in the analysis, whereas the other data show a strong preference for the standard cosmological model. We finally present constraints on common parameters marginalized over both cosmological models. For reasonable priors on models with and without EDE, the H0 tension is reduced by at least 20%.
△ Less
Submitted 8 January, 2024; v1 submitted 10 October, 2023;
originally announced October 2023.
-
Identifying frequency decorrelated dust residuals in B-mode maps by exploiting the spectral capability of bolometric interferometry
Authors:
M. Regnier,
E. Manzan,
J. -Ch Hamilton,
A. Mennella,
J. Errard,
L. Zapelli,
S. A. Torchinsky,
S. Paradiso,
E. Battistelli,
P. De Bernardis,
L. Colombo,
M. De Petris,
G. D'Alessandro,
B. Garcia,
M. Gervasi,
S. Masi,
L. Mousset,
N. Miron Granese,
C. O'Sullivan,
M. Piat,
E. Rasztocky,
G. E. Romero,
C. G. Scoccola,
M. Zannoni
Abstract:
Astrophysical polarized foregrounds represent the most critical challenge in Cosmic Microwave Background (CMB) B-mode experiments. Multi-frequency observations can be used to constrain astrophysical foregrounds to isolate the CMB contribution. However, recent observations indicate that foreground emission may be more complex than anticipated.
We investigate how the increased spectral resolution…
▽ More
Astrophysical polarized foregrounds represent the most critical challenge in Cosmic Microwave Background (CMB) B-mode experiments. Multi-frequency observations can be used to constrain astrophysical foregrounds to isolate the CMB contribution. However, recent observations indicate that foreground emission may be more complex than anticipated.
We investigate how the increased spectral resolution provided by band splitting in Bolometric Interferometry (BI) through a technique called spectral imaging can help control the foreground contamination in the case of unaccounted Galactic dust frequency decorrelation along the line-of-sight.
We focus on the next generation ground-based CMB experiment CMB-S4, and compare its anticipated sensitivities, frequency and sky coverage with a hypothetical version of the same experiment based on BI. We perform a Monte-Carlo analysis based on parametric component separation methods (FGBuster and Commander) and compute the likelihood on the recovered tensor-to-scalar ratio.
The main result of this analysis is that spectral imaging allows us to detect systematic uncertainties on r from frequency decorrelation when this effect is not accounted for in component separation. Conversely, an imager would detect a biased value of r and would be unable to spot the presence of a systematic effect. We find a similar result in the reconstruction of the dust spectral index, where we show that with BI we can measure more precisely the dust spectral index also when frequency decorrelation is present.
The in-band frequency resolution provided by BI allows us to identify dust LOS frequency decorrelation residuals where an imager of similar performance would fail. This opens the prospect to exploit this potential in the context of future CMB polarization experiments that will be challenged by complex foregrounds in their quest for B-modes detection.
△ Less
Submitted 21 February, 2024; v1 submitted 6 September, 2023;
originally announced September 2023.
-
Cosmoglobe: Towards end-to-end CMB cosmological parameter estimation without likelihood approximations
Authors:
J. R. Eskilt,
K. Lee,
D. J. Watts,
V. Anshul,
R. Aurlien,
A. Basyrov,
M. Bersanelli,
L. P. L. Colombo,
H. K. Eriksen,
K. S. F. Fornazier,
U. Fuskeland,
M. Galloway,
E. Gjerløw,
L. T. Hergt,
H. T. Ihle,
J. G. S. Lunde,
A. Marins,
S. K. Nerval,
S. Paradiso,
F. Rahman,
M. San,
N. -O. Stutzer,
I. K. Wehus
Abstract:
We implement support for a cosmological parameter estimation algorithm as proposed by Racine et al. (2016) in Commander, and quantify its computational efficiency and cost. For a semi-realistic simulation similar to Planck LFI 70 GHz, we find that the computational cost of producing one single sample is about 20 CPU-hours and that the typical Markov chain correlation length is $\sim$100 samples. T…
▽ More
We implement support for a cosmological parameter estimation algorithm as proposed by Racine et al. (2016) in Commander, and quantify its computational efficiency and cost. For a semi-realistic simulation similar to Planck LFI 70 GHz, we find that the computational cost of producing one single sample is about 20 CPU-hours and that the typical Markov chain correlation length is $\sim$100 samples. The net effective cost per independent sample is $\sim$2 000 CPU-hours, in comparison with all low-level processing costs of 812 CPU-hours for Planck LFI and WMAP in Cosmoglobe Data Release 1. Thus, although technically possible to run already in its current state, future work should aim to reduce the effective cost per independent sample by one order of magnitude to avoid excessive runtimes, for instance through multi-grid preconditioners and/or derivative-based Markov chain sampling schemes. This work demonstrates the computational feasibility of true Bayesian cosmological parameter estimation with end-to-end error propagation for high-precision CMB experiments without likelihood approximations, but it also highlights the need for additional optimizations before it is ready for full production-level analysis.
△ Less
Submitted 31 October, 2023; v1 submitted 27 June, 2023;
originally announced June 2023.
-
Cosmoglobe DR1 results. II. Constraints on isotropic cosmic birefringence from reprocessed WMAP and Planck LFI data
Authors:
J. R. Eskilt,
D. J. Watts,
R. Aurlien,
A. Basyrov,
M. Bersanelli,
M. Brilenkov,
L. P. L. Colombo,
H. K. Eriksen,
K. S. F. Fornazier,
C. Franceschet,
U. Fuskeland,
M. Galloway,
E. Gjerløw,
B. Hensley,
L. T. Hergt,
D. Herman,
H. T. Ihle,
K. Lee,
J. G. S. Lunde,
S. K. Nerval,
S. Paradiso,
S. K. Patel,
F. Rahman,
M. Regnier,
M. San
, et al. (6 additional authors not shown)
Abstract:
Cosmic birefringence is a parity-violating effect that might have rotated the plane of linearly polarized light of the cosmic microwave background (CMB) by an angle $β$ since its emission. This has recently been measured to be non-zero at a statistical significance of $3.6σ$ in the official Planck PR4 and 9-year WMAP data. In this work, we constrain $β$ using the reprocessed BeyondPlanck LFI and C…
▽ More
Cosmic birefringence is a parity-violating effect that might have rotated the plane of linearly polarized light of the cosmic microwave background (CMB) by an angle $β$ since its emission. This has recently been measured to be non-zero at a statistical significance of $3.6σ$ in the official Planck PR4 and 9-year WMAP data. In this work, we constrain $β$ using the reprocessed BeyondPlanck LFI and Cosmoglobe DR1 WMAP polarization maps. These novel maps have both lower systematic residuals and a more complete error description than the corresponding official products. Foreground $EB$ correlations could bias measurements of $β$, and while thermal dust $EB$ emission has been argued to be statistically non-zero, no evidence for synchrotron $EB$ power has been reported. Unlike the dust-dominated Planck HFI maps, the majority of the LFI and WMAP polarization maps are instead dominated by synchrotron emission. Simultaneously constraining $β$ and the polarization miscalibration angle, $α$, of each channel, we find a best-fit value of $β=0.35^{\circ}\pm0.70^{\circ}$ with LFI and WMAP data only. When including the Planck HFI PR4 maps, but fitting $β$ separately for dust-dominated, $β_{>70\,\mathrm{GHz}}$, and synchrotron-dominated channels, $β_{\leq 70\,\mathrm{GHz}}$, we find $β_{\leq 70\,\mathrm{GHz}}=0.53^{\circ}\pm0.28^\circ$. This differs from zero with a statistical significance of $1.9σ$, and the main contribution to this value comes from the LFI 70 GHz channel. While the statistical significances of these results are low on their own, the measurement derived from the LFI and WMAP synchrotron-dominated maps agrees with the previously reported HFI-dominated constraints, despite the very different astrophysical and instrumental systematics involved in all these experiments.
△ Less
Submitted 3 May, 2023;
originally announced May 2023.
-
Cosmoglobe DR1 results. I. Improved Wilkinson Microwave Anisotropy Probe maps through Bayesian end-to-end analysis
Authors:
D. J. Watts,
A. Basyrov,
J. R. Eskilt,
M. Galloway,
L. T. Hergt,
D. Herman,
H. T. Ihle,
S. Paradiso,
F. Rahman,
H. Thommesen,
R. Aurlien,
M. Bersanelli,
L. A. Bianchi,
M. Brilenkov,
L. P. L. Colombo,
H. K. Eriksen,
C. Franceschet,
U. Fuskeland,
E. Gjerløw,
B. Hensley,
G. A. Hoerning,
K. Lee,
J. G. S. Lunde,
A. Marins,
S. K. Nerval
, et al. (8 additional authors not shown)
Abstract:
We present Cosmoglobe Data Release 1, which implements the first joint analysis of WMAP and Planck LFI time-ordered data, processed within a single Bayesian end-to-end framework. This framework builds directly on a similar analysis of the LFI measurements by the BeyondPlanck collaboration, and approaches the CMB analysis challenge through Gibbs sampling of a global posterior distribution, simultan…
▽ More
We present Cosmoglobe Data Release 1, which implements the first joint analysis of WMAP and Planck LFI time-ordered data, processed within a single Bayesian end-to-end framework. This framework builds directly on a similar analysis of the LFI measurements by the BeyondPlanck collaboration, and approaches the CMB analysis challenge through Gibbs sampling of a global posterior distribution, simultaneously accounting for calibration, mapmaking, and component separation. The computational cost of producing one complete WMAP+LFI Gibbs sample is 812 CPU-hr, of which 603 CPU-hrs are spent on WMAP low-level processing; this demonstrates that end-to-end Bayesian analysis of the WMAP data is computationally feasible. We find that our WMAP posterior mean temperature sky maps and CMB temperature power spectrum are largely consistent with the official WMAP9 results. Perhaps the most notable difference is that our CMB dipole amplitude is $3366.2 \pm 1.4\ \mathrm{μK}$, which is $11\ \mathrm{μK}$ higher than the WMAP9 estimate and $2.5\ σ$ higher than BeyondPlanck; however, it is in perfect agreement with the HFI-dominated Planck PR4 result. In contrast, our WMAP polarization maps differ more notably from the WMAP9 results, and in general exhibit significantly lower large-scale residuals. We attribute this to a better constrained gain and transmission imbalance model. It is particularly noteworthy that the W-band polarization sky map, which was excluded from the official WMAP cosmological analysis, for the first time appears visually consistent with the V-band sky map. Similarly, the long standing discrepancy between the WMAP K-band and LFI 30 GHz maps is finally resolved, and the difference between the two maps appears consistent with instrumental noise at high Galactic latitudes. All maps and the associated code are made publicly available through the Cosmoglobe web page.
△ Less
Submitted 14 March, 2023;
originally announced March 2023.
-
Status of QUBIC, the Q&U Bolometer for Cosmology
Authors:
L. Mousset,
P. Ade,
A. Almela,
G. Amico,
L. H. Arnaldi,
J. Aumont,
S. Banfi,
E. S. Battistelli,
B. Bélier,
L. Bergé,
J. -Ph. Bernard,
P. de Bernardis,
M. Bersanelli,
J. Bonaparte,
J. D. Bonilla,
E. Bunn,
D. Buzi,
D. Camilieri,
F. Cavaliere,
P. Chanial,
C. Chapron,
S. Colombo,
F. Columbro,
A. Coppolecchia,
B. Costanza
, et al. (86 additional authors not shown)
Abstract:
The Q&U Bolometric Interferometer for Cosmology (QUBIC) is a novel kind of polarimeter optimized for the measurement of the B-mode polarization of the Cosmic Microwave Back-ground (CMB), which is one of the major challenges of observational cosmology. The signal is expected to be of the order of a few tens of nK, prone to instrumental systematic effects and polluted by various astrophysical foregr…
▽ More
The Q&U Bolometric Interferometer for Cosmology (QUBIC) is a novel kind of polarimeter optimized for the measurement of the B-mode polarization of the Cosmic Microwave Back-ground (CMB), which is one of the major challenges of observational cosmology. The signal is expected to be of the order of a few tens of nK, prone to instrumental systematic effects and polluted by various astrophysical foregrounds which can only be controlled through multichroic observations. QUBIC is designed to address these observational issues with a novel approach that combines the advantages of interferometry in terms of control of instrumental systematics with those of bolometric detectors in terms of wide-band, background-limited sensitivity.
△ Less
Submitted 6 October, 2022;
originally announced October 2022.
-
BeyondPlanck IV. On end-to-end simulations in CMB analysis -- Bayesian versus frequentist statistics
Authors:
M. Brilenkov,
K. S. F. Fornazier,
L. T. Hergt,
G. A. Hoerning,
A. Marins,
T. Murokoshi,
F. Rahman,
N. -O. Stutzer,
Y. Zhou,
F. B. Abdalla,
K. J. Andersen,
R. Aurlien,
R. Banerji,
A. Basyrov,
A. Battista,
M. Bersanelli,
S. Bertocco,
S. Bollanos,
L. P. L. Colombo,
H. K. Eriksen,
J. R. Eskilt,
M. K. Foss,
C. Franceschet,
U. Fuskeland,
S. Galeotta
, et al. (26 additional authors not shown)
Abstract:
End-to-end simulations play a key role in the analysis of any high-sensitivity CMB experiment, providing high-fidelity systematic error propagation capabilities unmatched by any other means. In this paper, we address an important issue regarding such simulations, namely how to define the inputs in terms of sky model and instrument parameters. These may either be taken as a constrained realization…
▽ More
End-to-end simulations play a key role in the analysis of any high-sensitivity CMB experiment, providing high-fidelity systematic error propagation capabilities unmatched by any other means. In this paper, we address an important issue regarding such simulations, namely how to define the inputs in terms of sky model and instrument parameters. These may either be taken as a constrained realization derived from the data, or as a random realization independent from the data. We refer to these as Bayesian and frequentist simulations, respectively. We show that the two options lead to significantly different correlation structures, as frequentist simulations, contrary to Bayesian simulations, effectively include cosmic variance, but exclude realization-specific correlations from non-linear degeneracies. Consequently, they quantify fundamentally different types of uncertainties, and we argue that they therefore also have different and complementary scientific uses, even if this dichotomy is not absolute. Before BeyondPlanck, most pipelines have used a mix of constrained and random inputs, and used the same hybrid simulations for all applications, even though the statistical justification for this is not always evident. BeyondPlanck represents the first end-to-end CMB simulation framework that is able to generate both types of simulations, and these new capabilities have brought this topic to the forefront. The Bayesian BeyondPlanck simulations and their uses are described extensively in a suite of companion papers. In this paper we consider one important applications of the corresponding frequentist simulations, namely code validation. That is, we generate a set of 1-year LFI 30 GHz frequentist simulations with known inputs, and use these to validate the core low-level BeyondPlanck algorithms; gain estimation, correlated noise estimation, and mapmaking.
△ Less
Submitted 9 September, 2022;
originally announced September 2022.
-
BeyondPlanck X. Planck LFI frequency maps with sample-based error propagation
Authors:
A. Basyrov,
A. -S. Suur-Uski,
L. P. L. Colombo,
J. R. Eskilt,
S. Paradiso,
K. J. Andersen,
R. Aurlien,
R. Banerji,
M. Bersanelli,
S. Bertocco,
M. Brilenkov,
M. Carbone,
H. K. Eriksen,
M. K. Foss,
C. Franceschet,
U. Fuskeland,
S. Galeotta,
M. Galloway,
S. Gerakakis,
E. Gjerløw,
B. Hensley,
D. Herman,
M. Iacobellis,
M. Ieronymaki,
H. T. Ihle
, et al. (15 additional authors not shown)
Abstract:
We present Planck LFI frequency sky maps derived within the BeyondPlanck framework. This framework draws samples from a global posterior distribution that includes instrumental, astrophysical and cosmological parameters, and the main product is an entire ensemble of frequency sky map samples. This ensemble allows for computationally convenient end-to-end propagation of low-level instrumental uncer…
▽ More
We present Planck LFI frequency sky maps derived within the BeyondPlanck framework. This framework draws samples from a global posterior distribution that includes instrumental, astrophysical and cosmological parameters, and the main product is an entire ensemble of frequency sky map samples. This ensemble allows for computationally convenient end-to-end propagation of low-level instrumental uncertainties into higher-level science products. We show that the two dominant sources of LFI instrumental systematic uncertainties are correlated noise and gain fluctuations, and the products presented here support - for the first time - full Bayesian error propagation for these effects at full angular resolution. We compare our posterior mean maps with traditional frequency maps delivered by the Planck collaboration, and find generally good agreement. The most important quality improvement is due to significantly lower calibration uncertainties in the new processing, as we find a fractional absolute calibration uncertainty at 70 GHz of $δg_{0}/g_{0} =5 \cdot 10^{-5}$, which is nominally 40 times smaller than that reported by Planck 2018. However, the original Planck 2018 estimate has a non-trivial statistical interpretation, and this further illustrates the advantage of the new framework in terms of producing self-consistent and well-defined error estimates of all involved quantities without the need of ad hoc uncertainty contributions. We describe how low-resolution data products, including dense pixel-pixel covariance matrices, may be produced directly from the posterior samples without the need for computationally expensive analytic calculations or simulations. We conclude that posterior-based frequency map sampling provides unique capabilities in terms of low-level systematics modelling and error propagation, and may play an important role for future CMB B-mode experiments. (Abridged.)
△ Less
Submitted 30 August, 2022;
originally announced August 2022.
-
BeyondPlanck XI. Bayesian CMB analysis with sample-based end-to-end error propagation
Authors:
L. P. L. Colombo,
J. R. Eskilt,
S. Paradiso,
H. Thommesen,
K. J. Andersen,
R. Aurlien,
R. Banerji,
M. Bersanelli,
S. Bertocco,
M. Brilenkov,
M. Carbone,
H. K. Eriksen,
M. K. Foss,
C. Franceschet,
U. Fuskeland,
S. Galeotta,
M. Galloway,
S. Gerakakis,
E. Gjerløw,
B. Hensley,
D. Herman,
M. Iacobellis,
M. Ieronymaki,
H. T. Ihle,
J. B. Jewell
, et al. (14 additional authors not shown)
Abstract:
We present posterior sample-based cosmic microwave background (CMB) constraints from Planck LFI and WMAP observations derived through global end-to-end Bayesian processing. We use these samples to study correlations between CMB, foreground, and instrumental parameters, and we identify a particularly strong degeneracy between CMB temperature fluctuations and free-free emission on intermediate angul…
▽ More
We present posterior sample-based cosmic microwave background (CMB) constraints from Planck LFI and WMAP observations derived through global end-to-end Bayesian processing. We use these samples to study correlations between CMB, foreground, and instrumental parameters, and we identify a particularly strong degeneracy between CMB temperature fluctuations and free-free emission on intermediate angular scales, which is mitigated through model reduction, masking, and resampling. We compare our posterior-based CMB results with previous Planck products, and find generally good agreement, but with higher noise due to exclusion of HFI data. We find a best-fit CMB dipole amplitude of $3362.7\pm1.4μK$, in excellent agreement with previous Planck results. The quoted uncertainty is derived directly from the sampled posterior distribution, and does not involve any ad hoc contribution for systematic effects. Similarly, we find a temperature quadrupole amplitude of $σ^{TT}_2=229\pm97μK^2$, in good agreement with previous results in terms of the amplitude, but the uncertainty is an order of magnitude larger than the diagonal Fisher uncertainty. Relatedly, we find lower evidence for a possible alignment between $\ell = 2$ and $\ell = 3$ than previously reported due to a much larger scatter in the individual quadrupole coefficients, caused both by marginalizing over a more complete set of systematic effects, and by our more conservative analysis mask. For higher multipoles, we find that the angular temperature power spectrum is generally in good agreement with both Planck and WMAP. This is the first time the sample-based asymptotically exact Blackwell-Rao estimator has been successfully established for multipoles up to $\ell\le600$, and it now accounts for the majority of the cosmologically important information. Cosmological parameter constraints are presented in a companion paper. (Abriged)
△ Less
Submitted 30 August, 2022;
originally announced August 2022.
-
From BeyondPlanck to Cosmoglobe: Open Science, Reproducibility, and Data Longevity
Authors:
S. Gerakakis,
M. Brilenkov,
M. Ieronymaki,
M. San,
D. J. Watts,
K. J. Andersen,
R. Aurlien,
R. Banerji,
A. Basyrov,
M. Bersanelli,
S. Bertocco,
M. Carbone,
L. P. L. Colombo,
H. K. Eriksen,
J. R. Eskilt,
M. K. Foss,
C. Franceschet,
U. Fuskeland,
S. Galeotta,
M. Galloway,
E. Gjerløw,
B. Hensley,
D. Herman,
M. Iacobellis,
H. T. Ihle
, et al. (17 additional authors not shown)
Abstract:
The BeyondPlanck and Cosmoglobe collaborations have implemented the first integrated Bayesian end-to-end analysis pipeline for CMB experiments. The primary long-term motivation for this work is to develop a common analysis platform that supports efficient global joint analysis of complementary radio, microwave, and sub-millimeter experiments. A strict prerequisite for this to succeed is broad part…
▽ More
The BeyondPlanck and Cosmoglobe collaborations have implemented the first integrated Bayesian end-to-end analysis pipeline for CMB experiments. The primary long-term motivation for this work is to develop a common analysis platform that supports efficient global joint analysis of complementary radio, microwave, and sub-millimeter experiments. A strict prerequisite for this to succeed is broad participation from the CMB community, and two foundational aspects of the program are therefore reproducibility and Open Science. In this paper, we discuss our efforts toward this aim. We also discuss measures toward facilitating easy code and data distribution, community-based code documentation, user-friendly compilation procedures, etc. This work represents the first publicly released end-to-end CMB analysis pipeline that includes raw data, source code, parameter files, and documentation. We argue that such a complete pipeline release should be a requirement for all major future and publicly-funded CMB experiments, noting that a full public release significantly increases data longevity by ensuring that the data quality can be improved whenever better processing techniques, complementary datasets, or more computing power become available, and thereby also taxpayers' value for money; providing only raw data and final products is not sufficient to guarantee full reproducibility in the future.
△ Less
Submitted 14 March, 2023; v1 submitted 20 May, 2022;
originally announced May 2022.
-
BeyondPlanck XII. Cosmological parameter constraints with end-to-end error propagation
Authors:
S. Paradiso,
L. P. L. Colombo,
K. J. Andersen,
R. Aurlien,
R. Banerji,
A. Basyrov,
M. Bersanelli,
S. Bertocco,
M. Brilenkov,
M. Carbone,
H. K. Eriksen,
J. R. Eskilt,
M. K. Foss,
C. Franceschet,
U. Fuskeland,
S. Galeotta,
M. Galloway,
S. Gerakakis,
E. Gjerløw,
B. Hensley,
D. Herman,
M. Iacobellis,
M. Ieronymaki,
H. T. Ihle,
J. B. Jewell
, et al. (16 additional authors not shown)
Abstract:
We present cosmological parameter constraints as estimated using the Bayesian BeyondPlanck (BP) analysis framework. This method supports seamless end-to-end error propagation from raw time-ordered data to final cosmological parameters. As a first demonstration of the method, we analyze time-ordered Planck LFI observations, combined with selected external data (WMAP 33-61GHz, Planck HFI DR4 353 and…
▽ More
We present cosmological parameter constraints as estimated using the Bayesian BeyondPlanck (BP) analysis framework. This method supports seamless end-to-end error propagation from raw time-ordered data to final cosmological parameters. As a first demonstration of the method, we analyze time-ordered Planck LFI observations, combined with selected external data (WMAP 33-61GHz, Planck HFI DR4 353 and 857GHz, and Haslam 408MHz) in the form of pixelized maps which are used to break critical astrophysical degeneracies. Overall, all results are generally in good agreement with previously reported values from Planck 2018 and WMAP, with the largest relative difference for any parameter of about 1 sigma when considering only temperature multipoles between 29<l<601. In cases where there are differences, we note that the BP results are generally slightly closer to the high-l HFI-dominated Planck 2018 results than previous analyses, suggesting slightly less tension between low and high multipoles. Using low-l polarization information from LFI and WMAP, we find a best-fit value of tau=0.066 +/- 0.013, which is higher than the low value of tau=0.051 +/- 0.006 derived from Planck 2018 and slightly lower than the value of 0.069 +/- 0.011 derived from joint analysis of official LFI and WMAP products. Most importantly, however, we find that the uncertainty derived in the BP processing is about 30% larger than when analyzing the official products, after taking into account the different sky coverage. We argue that this is due to marginalizing over a more complete model of instrumental and astrophysical parameters, and this results in both more reliable and more rigorously defined uncertainties. We find that about 2000 Monte Carlo samples are required to achieve robust convergence for low-resolution CMB covariance matrix with 225 independent modes.
△ Less
Submitted 20 May, 2022;
originally announced May 2022.
-
BeyondPlanck V. Minimal ADC Corrections for Planck LFI
Authors:
D. Herman,
R. A. Watson,
K. J. Andersen,
R. Aurlien,
R. Banjeri,
M. Bersanelli,
S. Bertocco,
M. Brilenkov,
M. Carbone,
L. P. L. Colombo,
H. K. Eriksen,
M. K. Foss,
C. Franceschet,
U. Fuskeland,
S. Galeotta,
M. Galloway,
S. Gerakakis,
E. Gjerløw,
B. Hensley,
M. Iacobellis,
M. Ieronymaki,
H. T. Ihle,
J. B. Jewell,
A. Karakci,
E. Keihänen
, et al. (14 additional authors not shown)
Abstract:
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 LFI Data Processing Center (DPC) analysis, and relies on the binned rms noise profile of each detector data stream. However, rather than building the corre…
▽ More
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 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 advange of this approach is that only samples which are significantly affected by ADC DNLs are corrected. 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 correctsion are applited 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 HFI, 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.
△ Less
Submitted 6 April, 2022; v1 submitted 25 March, 2022;
originally announced March 2022.
-
From BeyondPlanck to Cosmoglobe: Preliminary $\mathit{WMAP}$ $\mathit Q$-band analysis
Authors:
D. J. Watts,
M. Galloway,
H. T. Ihle,
K. J. Andersen,
R. Aurlien,
R. Banerji,
A. Basyrov,
M. Bersanelli,
S. Bertocco,
M. Brilenkov,
M. Carbone,
L. P. L. Colombo,
H. K. Eriksen,
J. R. Eskilt,
M. K. Foss,
C. Franceschet,
U. Fuskeland,
S. Galeotta,
S. Gerakakis,
E. Gjerløw,
B. Hensley,
D. Herman,
M. Iacobellis,
M. Ieronymaki,
J. B. Jewell
, et al. (18 additional authors not shown)
Abstract:
We present the first application of the Cosmoglobe analysis framework by analyzing 9-year $\mathit{WMAP}$ time-ordered observations using similar machinery as BeyondPlanck utilizes for $\mathit{Planck}$ LFI. We analyze only the $\mathit Q$-band (41 GHz) data and report on the low-level analysis process from uncalibrated time-ordered data to calibrated maps. Most of the existing BeyondPlanck pipeli…
▽ More
We present the first application of the Cosmoglobe analysis framework by analyzing 9-year $\mathit{WMAP}$ time-ordered observations using similar machinery as BeyondPlanck utilizes for $\mathit{Planck}$ LFI. We analyze only the $\mathit Q$-band (41 GHz) data and report on the low-level analysis process from uncalibrated time-ordered data to calibrated maps. Most of the existing BeyondPlanck pipeline may be reused for $\mathit{WMAP}$ analysis with minimal changes to the existing codebase. The main modification is the implementation of the same preconditioned biconjugate gradient mapmaker used by the $\mathit{WMAP}$ team. Producing a single $\mathit{WMAP}$ $\mathit Q$1-band sample requires 22 CPU-hrs, which is slightly more than the cost of a $\mathit{Planck}$ 44 GHz sample of 17 CPU-hrs; this demonstrates that full end-to-end Bayesian processing of the $\mathit{WMAP}$ data is computationally feasible. In general, our recovered maps are very similar to the maps released by the $\mathit{WMAP}$ team, although with two notable differences. In temperature we find a $\sim2\,\mathrm{μK}$ quadrupole difference that most likely is caused by different gain modeling, while in polarization we find a distinct $2.5\,\mathrm{μK}$ signal that has been previously called poorly-measured modes by the $\mathit{WMAP}$ team. In the Cosmoglobe processing, this pattern arises from temperature-to-polarization leakage from the coupling between the CMB Solar dipole, transmission imbalance, and sidelobes. No traces of this pattern are found in either the frequency map or TOD residual map, suggesting that the current processing has succeeded in modelling these poorly measured modes within the assumed parametric model by using $\mathit{Planck}$ information to break the sky-synchronous degeneracies inherent in the $\mathit{WMAP}$ scanning strategy.
△ Less
Submitted 23 May, 2022; v1 submitted 24 February, 2022;
originally announced February 2022.
-
BeyondPlanck XIII. Intensity foreground sampling, degeneracies, and priors
Authors:
K. J. Andersen,
D. Herman,
R. Aurlien,
R. Banerji,
A. Basyrov,
M. Bersanelli,
S. Bertocco,
M. Brilenkov,
M. Carbone,
L. P. L. Colombo,
H. K. Eriksen,
J. R. Eskilt,
M. K. Foss,
C. Franceschet,
U. Fuskeland,
S. Galeotta,
M. Galloway,
S. Gerakakis,
E. Gjerløw,
B. Hensley,
M. Iacobellis,
M. Ieronymaki,
H. T. Ihle,
J. B. Jewell,
A. Karakci
, et al. (19 additional authors not shown)
Abstract:
We present the intensity foreground algorithms and model employed within the BeyondPlanck analysis framework. The BeyondPlanck analysis is aimed at integrating component separation and instrumental parameter sampling within a global framework, leading to complete end-to-end error propagation in the $Planck$ Low Frequency Instrument (LFI) data analysis. Given the scope of the BeyondPlanck analysis,…
▽ More
We present the intensity foreground algorithms and model employed within the BeyondPlanck analysis framework. The BeyondPlanck analysis is aimed at integrating component separation and instrumental parameter sampling within a global framework, leading to complete end-to-end error propagation in the $Planck$ Low Frequency Instrument (LFI) data analysis. Given the scope of the BeyondPlanck analysis, a limited set of data is included in the component separation process, leading to foreground parameter degeneracies. In order to properly constrain the Galactic foreground parameters, we improve upon the previous $\texttt{Commander}$ component separation implementation by adding a suite of algorithmic techniques. These algorithms are designed to improve the stability and computational efficiency for weakly constrained posterior distributions. These are: 1) joint foreground spectral parameter and amplitude sampling, building on ideas from Miramare; 2) component-based monopole determination; 3) joint spectral parameter and monopole sampling; and 4) application of informative spatial priors for component amplitude maps. We find that the only spectral parameter with a significant signal-to-noise ratio using the current BeyondPlanck data set is the peak frequency of the anomalous microwave emission component, for which we find $ν_{\mathrm{p}}=25.3\pm0.5$ GHz; all others must be constrained through external priors. Future works will be aimed at integrating many more data sets into this analysis, both map and time-ordered based, thereby gradually eliminating the currently observed degeneracies in a controlled manner with respect to both instrumental systematic effects and astrophysical degeneracies. When this happens, the simple LFI-oriented data model employed in the current work will need to be generalized to account for both a richer astrophysical model and additional instrumental effects.
△ Less
Submitted 1 October, 2022; v1 submitted 20 January, 2022;
originally announced January 2022.
-
BeyondPlanck XVI. Limits on Large-Scale Polarized Anomalous Microwave Emission from Planck LFI and WMAP
Authors:
D. Herman,
B. Hensley,
K. J. Andersen,
R. Aurlien,
R. Banerji,
M. Bersanelli,
S. Bertocco,
M. Brilenkov,
M. Carbone,
L. P. L. Colombo,
H. K. Eriksen,
M. K. Foss,
C. Franceschet,
U. Fuskeland,
S. Galeotta,
M. Galloway,
S. Gerakakis,
E. Gjerløw,
M. Iacobellis,
M. Ieronymaki,
H. T. Ihle,
J. B. Jewell,
A. Karakci,
E. Keihänen,
R. Keskitalo
, et al. (13 additional authors not shown)
Abstract:
We constrain the level of polarized anomalous microwave emission (AME) on large angular scales using $\textit{Planck}$ LFI and $\textit{WMAP}$ polarization data within a Bayesian CMB analysis framework. We model synchrotron emission with a power-law spectral energy distribution, and the sum of AME and thermal dust emission through linear regression with the $\textit{Planck}$ HFI 353 GHz data. This…
▽ More
We constrain the level of polarized anomalous microwave emission (AME) on large angular scales using $\textit{Planck}$ LFI and $\textit{WMAP}$ polarization data within a Bayesian CMB analysis framework. We model synchrotron emission with a power-law spectral energy distribution, and the sum of AME and thermal dust emission through linear regression with the $\textit{Planck}$ HFI 353 GHz data. This template-based dust emission model allows us to constrain the level of polarized AME while making minimal assumptions on its frequency dependence. We neglect cosmic microwave background fluctuations, but show through simulations that these 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, and for priors steeper than $β_{\mathrm{s}} = -3.1\pm0.1$ we find an upper limit of $p_{\mathrm{AME}}^{\rm max}\lesssim 0.6\,\%$ ($95\,\%$ confidence). In contrast, for $β_{\mathrm{s}}=-3.0\pm0.1$, we find a nominal detection of $p_{\mathrm{AME}}=2.5\pm1.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 $β_\mathrm{s}$. Combining the current $\textit{Planck}$ and $\textit{WMAP}$ observations with measurements from high-sensitivity low-frequency experiments such as C-BASS and QUIJOTE will be critical to improve these limits further.
△ Less
Submitted 26 April, 2022; v1 submitted 10 January, 2022;
originally announced January 2022.
-
BeyondPlanck III. Commander3
Authors:
M. Galloway,
K. J. Andersen,
R. Aurlien,
R. Banerji,
M. Bersanelli,
S. Bertocco,
M. Brilenkov,
M. Carbone,
L. P. L. Colombo,
H. K. Eriksen,
M. K. Foss,
C. Franceschet,
U. Fuskeland,
S. Galeotta,
S. Gerakakis,
E. Gjerløw,
B. Hensley,
D. Herman,
M. Iacobellis,
M. Ieronymaki,
H. T. Ihle,
J. B. Jewell,
A. Karakci,
E. Keihänen,
R. Keskitalo
, et al. (13 additional authors not shown)
Abstract:
We describe the computational infrastructure for end-to-end Bayesian CMB analysis implemented by the BeyondPlanck collaboration. This code is called commander3, and provides a statistically consistent framework for global analysis of CMB and microwave observations, and may be useful for a wide range of legacy, current, and future experiments. The paper has three main goals. Firstly, we provide a h…
▽ More
We describe the computational infrastructure for end-to-end Bayesian CMB analysis implemented by the BeyondPlanck collaboration. This code is called commander3, and provides a statistically consistent framework for global analysis of CMB and microwave observations, and may be useful for a wide range of legacy, current, and future experiments. The paper has three main goals. Firstly, we provide a high-level overview of the existing code base, aiming to guide readers who wish to extend and adapt the code according to their own needs, or to reimplement it from scratch in a different programming language. Secondly, we discuss some critical computational challenges that arise within any global CMB analysis framework, for instance in-memory compression of time-ordered data, FFT optimization, and parallelization and load-balancing. Thirdly, we quantify the CPU and RAM requirements for the current BeyondPlanck analysis, and find that a total of 1.5 TB of RAM is required for efficient analysis, and the total cost of a full Gibbs sample is 170 CPU-hrs, including both low-level processing and high-level component separation, which is well within the capabilities of current low-cost computing facilities. The existing code base is made publicly available under a GNU General Public Library (GPL) license.
△ Less
Submitted 10 January, 2022;
originally announced January 2022.
-
BeyondPlanck VIII. Efficient Sidelobe Convolution and Correction through Spin Harmonics
Authors:
M. Galloway,
M. Reinecke,
K. J. Andersen,
R. Aurlien,
R. Banerji,
M. Bersanelli,
S. Bertocco,
M. Brilenkov,
M. Carbone,
L. P. L. Colombo,
H. K. Eriksen,
M. K. Foss,
C. Franceschet,
U. Fuskeland,
S. Galeotta,
S. Gerakakis,
E. Gjerløw,
B. Hensley,
D. Herman,
M. Iacobellis,
M. Ieronymaki,
H. T. Ihle,
J. B. Jewell,
A. Karakci,
E. Keihänen
, et al. (13 additional authors not shown)
Abstract:
We introduce a new formulation of the Conviqt convolution algorithm in terms of spin harmonics, and apply this to the problem of sidelobe correction for BeyondPlanck, the first end-to-end Bayesian Gibbs sampling framework for CMB analysis. We compare our implementation to the previous Planck LevelS implementation, and find good agreement between the two codes in terms of accuracy, but with a speed…
▽ More
We introduce a new formulation of the Conviqt convolution algorithm in terms of spin harmonics, and apply this to the problem of sidelobe correction for BeyondPlanck, the first end-to-end Bayesian Gibbs sampling framework for CMB analysis. We compare our implementation to the previous Planck LevelS implementation, and find good agreement between the two codes in terms of accuracy, but with a speed-up reaching a factor of 3--10, depending on the frequency bandlimits, $l_{\textrm{max}}$ and $m_{\textrm{max}}$. The new algorithm is significantly simpler to implement and maintain, since all low-level calculations are handled through an external spherical harmonic transform library. We find that our mean sidelobe estimates for Planck LFI agree well with previous efforts. Additionally, we present novel sidelobe rms maps that quantify the uncertainty in the sidelobe corrections due to variations in the sky model.
△ Less
Submitted 10 January, 2022;
originally announced January 2022.
-
BeyondPlanck X. Bandpass and beam leakage corrections
Authors:
T. L. Svalheim,
K. J. Andersen,
R. Aurlien,
R. Banerji,
M. Bersanelli,
S. Bertocco,
M. Brilenkov,
M. Carbone,
L. P. L. Colombo,
H. K. Eriksen,
M. K. Foss,
C. Franceschet,
U. Fuskeland,
S. Galeotta,
M. Galloway,
S. Gerakakis,
E. Gjerløw,
B. Hensley,
D. Herman,
M. Iacobellis,
M. Ieronymaki,
H. T. Ihle,
J. B. Jewell,
A. Karakci,
E. Keihänen
, et al. (14 additional authors not shown)
Abstract:
We discuss the treatment of bandpass and beam leakage corrections in the Bayesian BeyondPlanck CMB analysis pipeline as applied to the Planck LFI measurements. As a preparatory step, we first apply three corrections to the nominal LFI bandpass profiles including removal of a known systematic effect in the ground measuring equipment at 61 GHz; smoothing of standing wave ripples; and edge regulariza…
▽ More
We discuss the treatment of bandpass and beam leakage corrections in the Bayesian BeyondPlanck CMB analysis pipeline as applied to the Planck LFI measurements. As a preparatory step, we first apply three corrections to the nominal LFI bandpass profiles including removal of a known systematic effect in the ground measuring equipment at 61 GHz; smoothing of standing wave ripples; and edge regularization. The main net impact of these modifications is an overall shift in the 70 GHz bandpass of +0.6 GHz; we argue that any analysis of LFI data products, either from Planck or BeyondPlanck, should use these new bandpasses. In addition, we fit a single free bandpass parameter for each radiometer of the form $Δ_i = Δ_0 + δ_i$, where $Δ_0$ represents an absolute frequency shift per frequency band and $δ_i$ is a relative shift per detector. The absolute correction is only fitted at 30 GHz with a full $χ^2$-based likelihood, resulting in a correction of $Δ_{30}=0.24\pm0.03\,$GHz. The relative corrections are fitted using a spurious map approach, fundamentally similar to the method pioneered by the WMAP team, but without introducing many additional degrees of freedom. All bandpass parameters are sampled using a standard Metropolis sampler within the main BeyondPlanck Gibbs chain, and bandpass uncertainties are thus propagated to all other data products in the analysis. In total, we find that our bandpass model significantly reduces leakage effects. For beam leakage corrections, we adopt the official Planck LFI beam estimates without additional degrees of freedom, and only marginalize over the underlying sky model. We note that this is the first time leakage from beam mismatch has been included for Planck LFI maps.
△ Less
Submitted 10 January, 2022;
originally announced January 2022.
-
BeyondPlanck XV. Polarized foreground emission between 30 and 70 GHz
Authors:
T. L. Svalheim,
K. J. Andersen,
R. Aurlien,
R. Banerji,
M. Bersanelli,
S. Bertocco,
M. Brilenkov,
M. Carbone,
L. P. L. Colombo,
H. K. Eriksen,
M. K. Foss,
C. Franceschet,
U. Fuskeland,
S. Galeotta,
M. Galloway,
S. Gerakakis,
E. Gjerløw,
B. Hensley,
D. Herman,
M. Iacobellis,
M. Ieronymaki,
H. T. Ihle,
J. B. Jewell,
A. Karakci,
E. Keihänen
, et al. (13 additional authors not shown)
Abstract:
We constrain polarized foreground emission between 30 and 70 GHz with the Planck Low Frequency Instrument (LFI) and WMAP data within the global Bayesian BeyondPlanck framework. We combine for the first time full-resolution Planck LFI time-ordered data with low-resolution WMAP sky maps at 33, 40 and 61 GHz. Spectral parameters are fit with a likelihood defined at the native resolution of each frequ…
▽ More
We constrain polarized foreground emission between 30 and 70 GHz with the Planck Low Frequency Instrument (LFI) and WMAP data within the global Bayesian BeyondPlanck framework. We combine for the first time full-resolution Planck LFI time-ordered data with low-resolution WMAP sky maps at 33, 40 and 61 GHz. Spectral parameters are fit with a likelihood defined at the native resolution of each frequency channel. This analysis represents the first implementation of true multi-resolution component separation applied to CMB observations for both amplitude and spectral energy distribution (SED) parameters. For synchrotron emission, we approximate the SED as a power-law in frequency and find that the low signal-to-noise ratio of the current data strongly limits the number of free parameters that may be robustly constrained. We partition the sky into four large disjoint regions (High Latitude; Galactic Spur; Galactic Plane; and Galactic Center), each associated with its own power-law index. We find that the High Latitude region is prior-dominated, while the Galactic Center region is contaminated by residual instrumental systematics. The two remaining regions appear to be signal-dominated, and for these we derive spectral indices of $β_{\mathrm s}^{\mathrm{Spur}}=-3.17\pm0.06$ and $β_{\mathrm s}^{\mathrm{Plane}}=-3.03\pm0.07$, in good agreement with previous results. For thermal dust emission we assume a modified blackbody model and we fit a single power-law index across the full sky. We find $β_{\mathrm{d}}=1.64\pm0.03$, which is slightly steeper than reported from Planck HFI data, but still statistically consistent at the 2$σ$ confidence level.
△ Less
Submitted 18 January, 2022; v1 submitted 17 November, 2020;
originally announced November 2020.
-
BeyondPlanck VII. Bayesian estimation of gain and absolute calibration for CMB experiments
Authors:
E. Gjerløw,
H. T. Ihle,
S. Galeotta,
K. J. Andersen,
R. Aurlien,
R. Banerji,
M. Bersanelli,
S. Bertocco,
M. Brilenkov,
M. Carbone,
L. P. L. Colombo,
H. K. Eriksen,
M. K. Foss,
C. Franceschet,
U. Fuskeland,
M. Galloway,
S. Gerakakis,
B. Hensley,
D. Herman,
M. Iacobellis,
M. Ieronymaki,
J. B. Jewell,
A. Karakci,
E. Keihänen,
R. Keskitalo
, et al. (13 additional authors not shown)
Abstract:
We present a Bayesian calibration algorithm for CMB observations as implemented within the global end-to-end BeyondPlanck (BP) framework, and apply this to the Planck Low Frequency Instrument (LFI) data. Following the most recent Planck analysis, we decompose the full time-dependent gain into a sum of three orthogonal components: One absolute calibration term, common to all detectors; one time-ind…
▽ More
We present a Bayesian calibration algorithm for CMB observations as implemented within the global end-to-end BeyondPlanck (BP) framework, and apply this to the Planck Low Frequency Instrument (LFI) data. Following the most recent Planck analysis, we decompose the full time-dependent gain into a sum of three 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 is allowed to vary between one-hour pointing periods. Each term is 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 are 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 anomalous jumps in the instrumental gain. 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 DR4 (Planck 2018). The deviations we find are within expected error bounds, and we attribute them to differences in data usage and general approach between the pipelines. In particular, the BP calibration is performed globally, resulting in better inter-frequency consistency. Additionally, WMAP observations are used actively in the BP analysis, which breaks degeneracies in the Planck data set and results in better agreement with WMAP. Although our presentation and algorithm are currently oriented toward LFI processing, the procedure is fully generalizable to other experiments.
△ Less
Submitted 20 May, 2022; v1 submitted 16 November, 2020;
originally announced November 2020.
-
BeyondPlanck VI. Noise characterization and modelling
Authors:
H. T. Ihle,
M. Bersanelli,
C. Franceschet,
E. Gjerløw,
K. J. Andersen,
R. Aurlien,
R. Banerji,
S. Bertocco,
M. Brilenkov,
M. Carbone,
L. P. L. Colombo,
H. K. Eriksen,
J. R. Eskilt,
M. K. Foss,
U. Fuskeland,
S. Galeotta,
M. Galloway,
S. Gerakakis,
B. Hensley,
D. Herman,
M. Iacobellis,
M. Ieronymaki,
J. B. Jewell,
A. Karakci,
E. Keihänen
, et al. (16 additional authors not shown)
Abstract:
We present a Bayesian method for estimating instrumental noise parameters and propagating noise uncertainties within the global BeyondPlanck Gibbs sampling framework, and apply this to Planck Low Frequency Instrument (LFI) time-ordered data. Following previous literature, we initially adopt a $1/f$ model for the noise power spectral density (PSD), but find the need for an additional lognormal comp…
▽ More
We present a Bayesian method for estimating instrumental noise parameters and propagating noise uncertainties within the global BeyondPlanck Gibbs sampling framework, and apply this to Planck Low Frequency Instrument (LFI) time-ordered data. Following previous literature, we initially adopt a $1/f$ model for the noise power spectral density (PSD), but find the need for an additional lognormal component in the noise model for the 30 and 44\,GHz bands. We implement an optimal Wiener-filter (or constrained realization) gap-filling procedure to account for masked data. We then use this procedure to both estimate the gapless correlated noise in the time-domain, $n_\mathrm{corr}$, and to sample the noise PSD parameters, $ξ^n = \{σ_0, f_\mathrm{knee}, α, A_\mathrm{p}\}$. In contrast to previous \textit{Planck} analyses, we assume piecewise stationary noise only within each pointing period (PID), not throughout the full mission, but we adopt the LFI Data Processing Center (DPC) results as priors on $α$ and $f_\mathrm{knee}$. On average, we find best-fit correlated noise parameters that are mostly consistent with previous results, with a few notable exceptions. However, a detailed inspection of the time-dependent results reveals many important findings. First and foremost, we find strong evidence for statistically significant temporal variations in all noise PSD parameters, many of which are directly correlated with satellite housekeeping data. Second, while the simple $1/f$ model appears to be an excellent fit for the LFI 70 GHz channel, there is evidence for additional correlated noise not described by a $1/f$ model in the 30 and 44 GHz channels, including within the primary science frequency range of 0.1--1 Hz. (Abridged)
△ Less
Submitted 19 July, 2022; v1 submitted 12 November, 2020;
originally announced November 2020.
-
BeyondPlanck II. CMB map-making through Gibbs sampling
Authors:
E. Keihänen,
A. -S. Suur-Uski,
K. J. Andersen,
R. Aurlien,
R. Banerji,
M. Bersanelli,
S. Bertocco,
M. Brilenkov,
M. Carbone,
L. P. L. Colombo,
H. K. Eriksen,
M. K. Foss,
C. Franceschet,
U. Fuskeland,
S. Galeotta,
M. Galloway,
S. Gerakakis,
E. Gjerløw,
B. Hensley,
D. Herman,
M. Iacobellis,
M. Ieronymaki,
H. T. Ihle,
J. B. Jewell,
A. Karakci
, et al. (15 additional authors not shown)
Abstract:
We present a Gibbs sampling solution to the map-making problem for CMB measurements, building 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 benef…
▽ More
We present a Gibbs sampling solution to the map-making problem for CMB measurements, building 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 allows us for the first time to bring the destriping baseline length to a single sample. We apply the Gibbs procedure to simulated Planck 30 GHz data. We find that gaps in the time-ordered data are handled efficiently by filling them with simulated noise as part of the Gibbs process. The Gibbs procedure yields a chain of map samples, from which we may compute the posterior mean as a best-estimate map. The variation in the chain provides information on the correlated residual noise, without 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 total number of iterations. The conceptual advantages of the Gibbs sampling approach lies in statistically well-defined error propagation and systematic error correction, and this methodology forms the conceptual basis for the map-making algorithm employed in the BeyondPlanck framework, which implements the first end-to-end Bayesian analysis pipeline for CMB observations.
△ Less
Submitted 11 November, 2020;
originally announced November 2020.
-
BeyondPlanck I. Global Bayesian analysis of the Planck Low Frequency Instrument data
Authors:
BeyondPlanck Collaboration,
K. J. Andersen,
R. Aurlien,
R. Banerji,
A. Basyrov,
M. Bersanelli,
S. Bertocco,
M. Brilenkov,
M. Carbone,
L. P. L. Colombo,
H. K. Eriksen,
J. R. Eskilt,
M. K. Foss,
C. Franceschet,
U. Fuskeland,
S. Galeotta,
M. Galloway,
S. Gerakakis,
E. Gjerløw,
B. Hensley,
D. Herman,
M. Iacobellis,
M. Ieronymaki,
H. T. Ihle,
J. B. Jewell
, et al. (20 additional authors not shown)
Abstract:
We describe the BeyondPlanck project in terms of motivation, methodology and main products, and provide a guide to a set of companion papers that describe each result in fuller detail. We implement a complete end-to-end Bayesian analysis framework for the Planck LFI observations. The primary product is a full joint posterior distribution $P(ω|d)$, where $ω$ represents the set of all free instrumen…
▽ More
We describe the BeyondPlanck project in terms of motivation, methodology and main products, and provide a guide to a set of companion papers that describe each result in fuller detail. We implement a complete end-to-end Bayesian analysis framework for the Planck LFI observations. The primary product is a full joint posterior distribution $P(ω|d)$, where $ω$ represents the set of all free instrumental, astrophysical, and cosmological parameters. Notable advantages of this approach are seamless end-to-end propagation of uncertainties; accurate modeling of both astrophysical and instrumental effects in the most natural basis for each uncertain quantity; optimized computational costs with little or no need for intermediate human interaction between various analysis steps; and a complete overview of the entire analysis process within one single framework. We focus in particular on low-$\ell$ CMB polarization reconstruction with Planck LFI. We identify several important new effects that have not been accounted for in previous pipelines, including gain over-smoothing and time-variable and non-$1/f$ correlated noise in the 30 and 44 GHz channels. We find that all results are consistent with the $Λ$CDM model, and we constrain the reionization optical depth to $τ=0.066\pm0.013$, with a low-resolution $χ^2$ probability-to-exceed of 32%. This uncertainty is about 30% larger than the official pipelines, arising from taking into account a more complete instrumental model. The marginal CMB Solar dipole amplitude is $3362.7\pm1.4μ\mathrm{K}$, where the error bar is derived directly from the posterior distribution without the need of any ad-hoc instrumental corrections. We are currently not aware of any significant unmodelled systematic effects remaining in the Planck LFI data, and, for the first time, the 44 GHz channel is fully exploited. (Abridged.)
△ Less
Submitted 12 September, 2022; v1 submitted 11 November, 2020;
originally announced November 2020.
-
The large scale polarization explorer (LSPE) for CMB measurements: performance forecast
Authors:
The LSPE collaboration,
G. Addamo,
P. A. R. Ade,
C. Baccigalupi,
A. M. Baldini,
P. M. Battaglia,
E. S. Battistelli,
A. Baù,
P. de Bernardis,
M. Bersanelli,
M. Biasotti,
A. Boscaleri,
B. Caccianiga,
S. Caprioli,
F. Cavaliere,
F. Cei,
K. A. Cleary,
F. Columbro,
G. Coppi,
A. Coppolecchia,
F. Cuttaia,
G. D'Alessandro,
G. De Gasperis,
M. De Petris,
V. Fafone
, et al. (80 additional authors not shown)
Abstract:
[Abridged] The measurement of the polarization of the Cosmic Microwave Background radiation is one of the current frontiers in cosmology. In particular, the detection of the primordial B-modes, could reveal the presence of gravitational waves in the early Universe. The detection of such component is at the moment the most promising technique to probe the inflationary theory describing the very ear…
▽ More
[Abridged] The measurement of the polarization of the Cosmic Microwave Background radiation is one of the current frontiers in cosmology. In particular, the detection of the primordial B-modes, could reveal the presence of gravitational waves in the early Universe. The detection of such component is at the moment the most promising technique to probe the inflationary theory describing the very early evolution of the Universe. We present the updated performance forecast of the Large Scale Polarization Explorer (LSPE), a program dedicated to the measurement of the CMB polarization. LSPE is composed of two instruments: Strip, a radiometer-based telescope on the ground in Tenerife, and SWIPE (Short-Wavelength Instrument for the Polarization Explorer) a bolometer-based instrument designed to fly on a winter arctic stratospheric long-duration balloon. The program is among the few dedicated to observation of the Northern Hemisphere, while most of the international effort is focused into ground-based observation in the Southern Hemisphere. Measurements are currently scheduled in Winter 2021/22 for SWIPE, with a flight duration up to 15 days, and in Summer 2021 with two years observations for Strip. We describe the main features of the two instruments, identifying the most critical aspects of the design, in terms of impact into performance forecast. We estimate the expected sensitivity of each instrument and propagate their combined observing power to the sensitivity to cosmological parameters, including the effect of scanning strategy, component separation, residual foregrounds and partial sky coverage. We also set requirements on the control of the most critical systematic effects and describe techniques to mitigate their impact. LSPE can reach a sensitivity in tensor-to-scalar ratio of $σ_r<0.01$, and improve constrains on other cosmological parameters.
△ Less
Submitted 9 August, 2021; v1 submitted 25 August, 2020;
originally announced August 2020.
-
Exploring Cosmic Origins with CORE: Survey requirements and mission design
Authors:
J. Delabrouille,
P. de Bernardis,
F. R. Bouchet,
A. Achúcarro,
P. A. R. Ade,
R. Allison,
F. Arroja,
E. Artal,
M. Ashdown,
C. Baccigalupi,
M. Ballardini,
A. J. Banday,
R. Banerji,
D. Barbosa,
J. Bartlett,
N. Bartolo,
S. Basak,
J. J. A. Baselmans,
K. Basu,
E. S. Battistelli,
R. Battye,
D. Baumann,
A. Benoît,
M. Bersanelli,
A. Bideaud
, et al. (178 additional authors not shown)
Abstract:
Future observations of cosmic microwave background (CMB) polarisation have the potential to answer some of the most fundamental questions of modern physics and cosmology. In this paper, we list the requirements for a future CMB polarisation survey addressing these scientific objectives, and discuss the design drivers of the CORE space mission proposed to ESA in answer to the "M5" call for a medium…
▽ More
Future observations of cosmic microwave background (CMB) polarisation have the potential to answer some of the most fundamental questions of modern physics and cosmology. In this paper, we list the requirements for a future CMB polarisation survey addressing these scientific objectives, and discuss the design drivers of the CORE space mission proposed to ESA in answer to the "M5" call for a medium-sized mission. The rationale and options, and the methodologies used to assess the mission's performance, are of interest to other future CMB mission design studies. CORE is designed as a near-ultimate CMB polarisation mission which, for optimal complementarity with ground-based observations, will perform the observations that are known to be essential to CMB polarisation scienceand cannot be obtained by any other means than a dedicated space mission.
△ Less
Submitted 14 June, 2017;
originally announced June 2017.
-
Exploring Cosmic Origins with CORE: Cosmological Parameters
Authors:
Eleonora Di Valentino,
Thejs Brinckmann,
Martina Gerbino,
Vivian Poulin,
François R. Bouchet,
Julien Lesgourgues,
Alessandro Melchiorri,
Jens Chluba,
Sebastien Clesse,
Jacques Delabrouille,
Cora Dvorkin,
Francesco Forastieri,
Silvia Galli,
Deanna C. Hooper,
Massimiliano Lattanzi,
Carlos J. A. P. Martins,
Laura Salvati,
Giovanni Cabass,
Andrea Caputo,
Elena Giusarma,
Eric Hivon,
Paolo Natoli,
Luca Pagano,
Simone Paradiso,
Jose Alberto Rubino-Martin
, et al. (103 additional authors not shown)
Abstract:
We forecast the main cosmological parameter constraints achievable with the CORE space mission which is dedicated to mapping the polarisation of the Cosmic Microwave Background (CMB). CORE was recently submitted in response to ESA's fifth call for medium-sized mission proposals (M5). Here we report the results from our pre-submission study of the impact of various instrumental options, in particul…
▽ More
We forecast the main cosmological parameter constraints achievable with the CORE space mission which is dedicated to mapping the polarisation of the Cosmic Microwave Background (CMB). CORE was recently submitted in response to ESA's fifth call for medium-sized mission proposals (M5). Here we report the results from our pre-submission study of the impact of various instrumental options, in particular the telescope size and sensitivity level, and review the great, transformative potential of the mission as proposed. Specifically, we assess the impact on a broad range of fundamental parameters of our Universe as a function of the expected CMB characteristics, with other papers in the series focusing on controlling astrophysical and instrumental residual systematics. In this paper, we assume that only a few central CORE frequency channels are usable for our purpose, all others being devoted to the cleaning of astrophysical contaminants. On the theoretical side, we assume LCDM as our general framework and quantify the improvement provided by CORE over the current constraints from the Planck 2015 release. We also study the joint sensitivity of CORE and of future Baryon Acoustic Oscillation and Large Scale Structure experiments like DESI and Euclid. Specific constraints on the physics of inflation are presented in another paper of the series. In addition to the six parameters of the base LCDM, which describe the matter content of a spatially flat universe with adiabatic and scalar primordial fluctuations from inflation, we derive the precision achievable on parameters like those describing curvature, neutrino physics, extra light relics, primordial helium abundance, dark matter annihilation, recombination physics, variation of fundamental constants, dark energy, modified gravity, reionization and cosmic birefringence. (ABRIDGED)
△ Less
Submitted 5 April, 2017; v1 submitted 30 November, 2016;
originally announced December 2016.