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.
We revisit the impact of finite time responses of bolometric detectors used for deep observations of the cosmic microwave background (CMB). Until now, bolometer transfer functions have been accounted ...for through a two-step procedure by first deconvolving an estimate of their Fourier-space representation from the raw time-ordered data (TOD), and then averaging the deconvolved TOD into pixelized maps. However, for many experiments, including the Planck High Frequency Instrument (HFI), it is necessary to apply an additional low-pass filter to avoid an excessive noise boost, which leads to an asymmetric effective beam. In this paper we demonstrate that this effect can be avoided if the transfer function deconvolution and pixelization operations are performed simultaneously through integrated maximum likelihood mapmaking. The resulting algorithm is structurally identical to the artDeco algorithm introduced by Keih\"anen & Reinecke (2012) for beam deconvolution. We illustrate the relevance of this method with simulated Planck HFI 143 GHz data, and find that the resulting effective beam is both more symmetric than with the two-step procedure, resulting in a sky-averaged ellipticity that is 64 % lower, and an effective beam full-width-at-half-maximum (FWHM) that is 2.3 % smaller. Similar improvements are expected for any other bolometer-based CMB experiments with long time constants.
LiteBIRD is a planned JAXA-led cosmic microwave background (CMB)
B
-mode satellite experiment aiming for launch in the late 2020s, with a primary goal of detecting the imprint of primordial ...inflationary gravitational waves. Its current baseline focal-plane configuration includes 15 frequency bands between 40 and 402 GHz, fulfilling the mission requirements to detect the amplitude of gravitational waves with the total uncertainty on the tensor-to-scalar ratio,
δr
, down to
δr
< 0.001. A key aspect of this performance is accurate astrophysical component separation, and the ability to remove polarized thermal dust emission is particularly important. In this paper we note that the CMB frequency spectrum falls off nearly exponentially above 300 GHz relative to the thermal dust spectral energy distribution, and a relatively minor high frequency extension can therefore result in even lower uncertainties and better model reconstructions. Specifically, we compared the baseline design with five extended configurations, while varying the underlying dust modeling, in each of which the High-Frequency Telescope (HFT) frequency range was shifted logarithmically toward higher frequencies, with an upper cutoff ranging between 400 and 600 GHz. In each case, we measured the tensor-to-scalar ratio
r
uncertainty and bias using both parametric and minimum-variance component-separation algorithms. When the thermal dust sky model includes a spatially varying spectral index and temperature, we find that the statistical uncertainty on
r
after foreground cleaning may be reduced by as much as 30–50% by extending the upper limit of the frequency range from 400 to 600 GHz, with most of the improvement already gained at 500 GHz. We also note that a broader frequency range leads to higher residuals when fitting an incorrect dust model, but also it is easier to discriminate between models through higher
χ
2
sensitivity. Even in the case in which the fitting procedure does not correspond to the underlying dust model in the sky, and when the highest frequency data cannot be modeled with sufficient fidelity and must be excluded from the analysis, the uncertainty on
r
increases by only about 5% for a 500 GHz configuration compared to the baseline.
We present the first application of the C
OSMOGLOBE
analysis framework by analyzing nine-year WMAP time-ordered observations that uses similar machinery to that of B
EYOND
P
LANCK
for the
Planck
Low ...Frequency Instrument (LFI). We analyzed only the
Q
-band (41 GHz) data and report on the low-level analysis process based on uncalibrated time-ordered data to calibrated maps. Most of the existing B
EYOND
P
LANCK
pipeline may be reused for 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 WMAP team. Producing a single WMAP
Q
1-band sample requires 22 CPU-hrs, which is slightly more than the cost of a
Planck
44 GHz sample of 17 CPU-hrs; this demonstrates that a full end-to-end Bayesian processing of the WMAP data is computationally feasible. In general, our recovered maps are very similar to the maps released by the WMAP team, although with two notable differences. In terms of temperature, we find a ∼2 μK quadrupole difference that most likely is caused by different gain modeling, while in polarization we find a distinct 2.5 μK signal that has been previously referred to as poorly measured modes by the WMAP team. In the C
OSMOGLOBE
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 modeling these poorly measured modes within the assumed parametric model by using
Planck
information to break the sky-synchronous degeneracies inherent in the WMAP scanning strategy.
BEYONDPLANCK Andersen, K. J.; Aurlien, R.; Banerji, R. ...
Astronomy and astrophysics (Berlin),
06/2023, Letnik:
675
Journal Article
Recenzirano
Odprti dostop
We describe the B
EYOND
P
LANCK
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 C
OSMOGLOBE
effort, which aims to combine state-of-the-art radio, microwave, and submillimeter data sets into one global astrophysical model.
