COSMOGLOBE DR1 results Watts, D. J.; Basyrov, A.; Eskilt, J. R. ...
Astronomy and astrophysics (Berlin),
11/2023, Letnik:
679
Journal Article
Recenzirano
Odprti dostop
We present C
OSMOGLOBE
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 ...directly builds on a similar analysis of the LFI measurements by the B
EYOND
P
LANCK
collaboration, and approaches the cosmic microwave background (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-h, of which 603 CPU-h 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 ± 1.4 μK, which is 11 μK higher than the WMAP9 estimate and 2.5
σ
higher than B
EYOND
P
LANCK
; 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. Relatedly, these updated maps allowed us for the first time to combine WMAP and LFI polarization data into a single coherent model of large-scale polarized synchrotron emission. Still, we identified a few issues that require additional work, including (1) low-level noise modeling; (2) large-scale temperature residuals at the 1–2 μK level; and (3) a strong degeneracy between the absolute
K
-band calibration and the dipole of the anomalous microwave emission component. We conclude that leveraging the complementary strengths of WMAP and LFI has allowed the mitigation of both experiments’ weaknesses, and resulted in new state-of-the-art WMAP sky maps. All maps and the associated code are made publicly available through the C
OSMOGLOBE
web page.
COSMOGLOBE DR1 results Eskilt, J. R.; Watts, D. J.; Aurlien, R. ...
Astronomy and astrophysics (Berlin),
11/2023, Letnik:
679
Journal Article
Recenzirano
Odprti dostop
Cosmic birefringence is a parity-violating effect that might have rotated the plane of the linearly polarized light of the cosmic microwave background (CMB) by an angle
β
since its emission. This ...angle has recently been measured to be nonzero 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 B
EYOND
P
LANCK
LFI and C
OSMOGLOBE
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 nonzero, 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° ±0.70° with LFI and WMAP data only. When including the
Planck
HFI PR4 maps, but fitting
β
separately for dust-dominated,
β
> 70 GHz
, and synchrotron-dominated channels,
β
≤70 GHz
, we find
β
≤70 GHz
= 0.53° ±0.28°. 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.
We implement support for a cosmological parameter estimation algorithm 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 ∼100 samples. The net effective cost per independent sample is ∼2000 CPU-hours, in comparison with all low-level processing costs of 812 CPU-hours for
Planck
LFI and WMAP in C
OSMOGLOBE
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.
BEYONDPLANCK Andersen, K. J.; Herman, D.; Aurlien, R. ...
Astronomy and astrophysics (Berlin),
06/2023, Letnik:
675
Journal Article
Recenzirano
Odprti dostop
We present the intensity foreground algorithms and model employed within the B
EYOND
P
LANCK
analysis framework. The B
EYOND
P
LANCK
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 B
EYOND
P
LANCK
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
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 M
IRAMARE
; (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 B
EYOND
P
LANCK
data set is the peak frequency of the anomalous microwave emission component, for which we find
ν
p
= 25.3 ± 0.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. This work will be organized within the Open Science-based C
OSMOGLOBE
community effort.
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.
COSMOGLOBE DR1 results Watts, D. J.; Fuskeland, U.; Aurlien, R. ...
Astronomy and astrophysics (Berlin),
06/2024, Letnik:
686
Journal Article
Recenzirano
We present the first model of full-sky polarized synchrotron emission that is derived from all WMAP and Planck LFI frequency maps. The basis of this analysis is the set of end-to-end reprocessed C ...OSMOGLOBE Data Release 1 (DR1) sky maps presented in a companion paper, which have significantly lower instrumental systematics than the legacy products from each experiment. We find that the resulting polarized synchrotron amplitude map has an average noise rms per 2° full width at half maximum (FWHM) beam of 3.2 μK at 30 GHz. This is 30% lower than the recently released B EYOND P LANCK model that included only LFI+WMAP Ka – V data, and 29% lower than the WMAP K -band map alone. The mean B -to- E power spectrum ratio is 0.39 ± 0.02, with amplitudes consistent with those measured previously by Planck and QUIJOTE. Assuming a power law model for the synchrotron spectral energy distribution and using the T – T plot method, we find a full-sky inverse noise-variance-weighted mean of the synchrotron polarized spectral index of β s = −3.07 ± 0.07 from the C OSMOGLOBE DR1 K band and 30 GHz, in good agreement with previous estimates. In summary, the novel C OSMOGLOBE DR1 synchrotron model is both more sensitive and systematically cleaner than similar previous models, and it has a more complete error description that is defined by a set of Monte Carlo posterior samples. We believe that these products are preferable over previous Planck and WMAP products for all synchrotron-related scientific applications, including simulations, forecasting, and component separation.
COSMOGLOBE DR1 results Watts, D J; Basyrov, A; Eskilt, J R ...
Astronomy & astrophysics,
11/2023, Letnik:
679
Journal Article
Recenzirano
Odprti dostop
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 ...directly builds on a similar analysis of the LFI measurements by the BEYONDPLANCK collaboration, and approaches the cosmic microwave background (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-h, of which 603 CPU-h 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 ± 1.4 μK, which is 11 μ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. Relatedly, these updated maps allowed us for the first time to combine WMAP and LFI polarization data into a single coherent model of large-scale polarized synchrotron emission. Still, we identified a few issues that require additional work, including (1) low-level noise modeling; (2) large-scale temperature residuals at the 1–2 μK level; and (3) a strong degeneracy between the absolute K-band calibration and the dipole of the anomalous microwave emission component. We conclude that leveraging the complementary strengths of WMAP and LFI has allowed the mitigation of both experiments’ weaknesses, and resulted in new state-of-the-art WMAP sky maps. All maps and the associated code are made publicly available through the COSMOGLOBE web page.
COSMOGLOBE DR1 results Eskilt, J R; Watts, D J; Aurlien, R ...
Astronomy & astrophysics,
11/2023, Letnik:
679
Journal Article
Recenzirano
Odprti dostop
Cosmic birefringence is a parity-violating effect that might have rotated the plane of the linearly polarized light of the cosmic microwave background (CMB) by an angle β since its emission. This ...angle has recently been measured to be nonzero 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 nonzero, 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° ±0.70° with LFI and WMAP data only. When including the Planck HFI PR4 maps, but fitting β separately for dust-dominated, β> 70 GHz, and synchrotron-dominated channels, β≤70 GHz, we find β≤70 GHz = 0.53° ±0.28°. 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.
BEYONDPLANCK Andersen, K J; Aurlien, R; Banerji, R ...
Astronomy & astrophysics,
07/2023, Letnik:
675
Journal Article
Recenzirano
Odprti dostop
We describe the BEYONDPLANCK project in terms of our motivation, methodology, and main products, and provide a guide to a set of companion papers that describe each result in more detail. Building ...directly on experience from ESA’s Planck mission, we implemented a complete end-to-end Bayesian analysis framework for the Planck Low Frequency Instrument (LFI) observations. The primary product is a full joint posterior distribution P(ω ∣ d), where ω represents the set of all free instrumental (gain, correlated noise, bandpass, etc.), astrophysical (synchrotron, free-free, thermal dust emission, etc.), and cosmological (cosmic microwave background – CMB – map, power spectrum, etc.) parameters. Some notable advantages of this approach compared to a traditional pipeline procedure 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. As a practical demonstration of this framework, we focus in particular on low-ℓ CMB polarization reconstruction with Planck LFI. In this process, 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. Modeling and mitigating both previously known and newly discovered systematic effects, we find that all results are consistent with the ΛCDM model, and we constrained the reionization optical depth to τ = 0.066 ± 0.013, with a low-resolution CMB-based χ2 probability to exceed of 32%. This uncertainty is about 30% larger than the official pipelines, arising from taking a more complete instrumental model into account. The marginal CMB solar dipole amplitude is 3362.7 ± 1.4 μK, where the error bar was derived directly from the posterior distribution without the need of any ad hoc instrumental corrections. We are currently not aware of any significant unmodeled systematic effects remaining in the Planck LFI data, and, for the first time, the 44 GHz channel is fully exploited in the current analysis. We argue that this framework can play a central role in the analysis of many current and future high-sensitivity CMB experiments, including LiteBIRD, and it will serve as the computational foundation of the emerging community-wide COSMOGLOBE effort, which aims to combine state-of-the-art radio, microwave, and submillimeter data sets into one global astrophysical model.
BEYONDPLANCK Andersen, K J; Herman, D; Aurlien, R ...
Astronomy & astrophysics,
07/2023, Letnik:
675
Journal Article
Recenzirano
Odprti dostop
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 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 νp = 25.3 ± 0.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. This work will be organized within the Open Science-based COSMOGLOBE community effort.