Showing 1–35 of 35 results for author: Jewell, J B

Thermal control of long delay lines in a high-resolution astrophotonic spectrograph

Authors: Gregory P. Sercel, Pradip R. Gatkine, Nemanja Jovanovic, Jeffrey B. Jewell, Luis Pereira da Costa, J. Kent Wallace, Dimitri P. Mawet

Abstract: High-resolution astronomical spectroscopy carried out with a photonic Fourier transform spectrograph (FTS) requires long asymmetrical optical delay lines that can be dynamically tuned. For example, to achieve a spectral resolution of R = 30,000, a delay line as long as 1.5 cm would be required. Such delays are inherently prone to phase errors caused by temperature fluctuations. This is due to the… ▽ More High-resolution astronomical spectroscopy carried out with a photonic Fourier transform spectrograph (FTS) requires long asymmetrical optical delay lines that can be dynamically tuned. For example, to achieve a spectral resolution of R = 30,000, a delay line as long as 1.5 cm would be required. Such delays are inherently prone to phase errors caused by temperature fluctuations. This is due to the relatively large thermo-optic coefficient and long lengths of the waveguides, in this case composed of SiN, resulting in thermally dependent changes to the optical path length. To minimize phase error to the order of 0.05 radians, thermal stability of the order of 0.05° C is necessary. A thermal control system capable of stability such as this would require a fast thermal response and minimal overshoot/undershoot. With a PID temperature control loop driven by a Peltier cooler and thermistor, we minimized interference fringe phase error to +/- 0.025 radians and achieved temperature stability on the order of 0.05° C. We present a practical system for precision temperature control of a foundry-fabricated and packaged FTS device on a SiN platform with delay lines ranging from 0.5 to 1.5 cm in length using inexpensive off-the-shelf components, including design details, control loop optimization, and considerations for thermal control of integrated photonics. △ Less

Submitted 2 April, 2024; v1 submitted 15 March, 2024; originally announced March 2024.

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.

Journal ref: A&A 675, A4 (2023)

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.

Comments: 32 pages, 23 figures, data available from https://www.cosmoglobe.uio.no/

Journal ref: A&A 675, A10 (2023)

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.

Comments: 26 pages, 24 figures. Submitted to A&A. Part of the BeyondPlanck paper suite

Journal ref: A&A 675, A11 (2023)

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.

Comments: 13 pages, 10 figures

MSC Class: 83F05

Journal ref: A&A 675, A12 (2023)

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.

Comments: 9 pages, 8 figures, mild language edits

Journal ref: A&A 675, A5 (2023)

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.

Comments: 11 figures, submitted to A&A. Includes updated instrument model and changes addressing referee comments

Journal ref: A&A 675, A16 (2023)

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.

Comments: 27 pages, 23 figures, 3 tables, part 13 of the BeyondPlanck release. All BeyondPlanck products and software will be released publicly at http://beyondplanck.science. Paper version as accepted by A&A

Journal ref: A&A 675, A13 (2023)

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.

Comments: 10 pages, 9 figures, part 15 of the BeyondPlanck release. All BeyondPlanck products and software will be released publicly at http://beyondplanck.science. Refereed, edited, and accepted by A&A. Added references, added arxiv IDs to BeyondPlanck papers. Figure update. Expanded discussion on AME models. Moved title from BeyondPlanck XV to BeyondPlanck XVI due to internal reordering

Journal ref: A&A 675, A15 (2023)

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.

Comments: 16 Pages, 7 Figures. Part of the BeyondPlanck paper suite

Journal ref: A&A 675, A3 (2023)

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.

Comments: 9 pages, 8 figures. Part of the BeyondPlanck paper suite

Journal ref: A&A 675, A8 (2023)

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.

Comments: 14 pages, 11 figures, this is part 10 in the BeyondPlanck release. All BeyondPlanck products and software will be released publicly at http://beyondplanck.science. Submitted to A&A

Journal ref: A&A 675, A9 (2023)

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.

Comments: 18 pages, 15 figures. All BeyondPlanck products and software will be released publicly at http://beyondplanck.science. This paper has been submitted to A&A. Major revision, matches submitted version

Journal ref: A&A 675, A14 (2023)

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.

Comments: 19 pages, 18 figures. All BeyondPlanck products and software will be released publicly at http://beyondplanck.science

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.

Comments: 26 pages, 29 figures, as accepted in A&A

Journal ref: A&A 675, A6 (2023)

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.

Comments: 11 pages, 10 figures. All BeyondPlanck products and software will be released publicly at http://beyondplanck.science during the online release conference (November 18-20, 2020). Connection details will be made available at the same website. Registration is mandatory for the online tutorial, but optional for the conference

Journal ref: A&A 675, A2 (2023)

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.

Comments: 61 pages, 31 figures, submitted to A&A. Major revision: Previously reported unresolved effects are now understood and mitigated. This paper version concludes the BeyondPlanck project. All BeyondPlanck products and software are available at https://cosmoglobe.uio.no

Journal ref: A&A 675, A1 (2023)

Cosmological Parameters from CMB Maps without Likelihood Approximation

Authors: Benjamin Racine, Jeffrey B. Jewell, Hans Kristian K. Eriksen, Ingunn K. Wehus

Abstract: We propose an efficient Bayesian MCMC algorithm for estimating cosmological parameters from CMB data without use of likelihood approximations. It builds on a previously developed Gibbs sampling framework that allows for exploration of the joint CMB sky signal and power spectrum posterior, P(s,Cl|d), and addresses a long-standing problem of efficient parameter estimation simultaneously in high and… ▽ More We propose an efficient Bayesian MCMC algorithm for estimating cosmological parameters from CMB data without use of likelihood approximations. It builds on a previously developed Gibbs sampling framework that allows for exploration of the joint CMB sky signal and power spectrum posterior, P(s,Cl|d), and addresses a long-standing problem of efficient parameter estimation simultaneously in high and low signal-to-noise regimes. To achieve this, our new algorithm introduces a joint Markov Chain move in which both the signal map and power spectrum are synchronously modified, by rescaling the map according to the proposed power spectrum before evaluating the Metropolis-Hastings accept probability. Such a move was already introduced by Jewell et al. (2009), who used it to explore low signal-to-noise posteriors. However, they also found that the same algorithm is inefficient in the high signal-to-noise regime, since a brute-force rescaling operation does not account for phase information. This problem is mitigated in the new algorithm by subtracting the Wiener filter mean field from the proposed map prior to rescaling, leaving high signal-to-noise information invariant in the joint step, and effectively only rescaling the low signal-to-noise component. To explore the full posterior, the new joint move is then interleaved with a standard conditional Gibbs sky map move. We apply our new algorithm to simplified simulations for which we can evaluate the exact posterior to study both its accuracy and performance, and find good agreement with the exact posterior; marginal means agree to less than 0.006 sigma, and standard deviations to better than 3%. The Markov Chain correlation length is of the same order of magnitude as those obtained by other standard samplers in the field. △ Less

Submitted 28 March, 2016; v1 submitted 21 December, 2015; originally announced December 2015.

Comments: 9 pages, 3 figures, Published in ApJ

Journal ref: B. Racine et al 2016 ApJ 820 31

Optimized Large-Scale CMB Likelihood And Quadratic Maximum Likelihood Power Spectrum Estimation

Authors: E. Gjerløw, L. P. L. Colombo, H. K. Eriksen, K. M. Górski, A. Gruppuso, J. B. Jewell, S. Plaszczynski, I. K. Wehus

Abstract: We revisit the problem of exact CMB likelihood and power spectrum estimation with the goal of minimizing computational cost through linear compression. This idea was originally proposed for CMB purposes by Tegmark et al.\ (1997), and here we develop it into a fully working computational framework for large-scale polarization analysis, adopting \WMAP\ as a worked example. We compare five different… ▽ More We revisit the problem of exact CMB likelihood and power spectrum estimation with the goal of minimizing computational cost through linear compression. This idea was originally proposed for CMB purposes by Tegmark et al.\ (1997), and here we develop it into a fully working computational framework for large-scale polarization analysis, adopting \WMAP\ as a worked example. We compare five different linear bases (pixel space, harmonic space, noise covariance eigenvectors, signal-to-noise covariance eigenvectors and signal-plus-noise covariance eigenvectors) in terms of compression efficiency, and find that the computationally most efficient basis is the signal-to-noise eigenvector basis, which is closely related to the Karhunen-Loeve and Principal Component transforms, in agreement with previous suggestions. For this basis, the information in 6836 unmasked \WMAP\ sky map pixels can be compressed into a smaller set of 3102 modes, with a maximum error increase of any single multipole of 3.8\% at $\ell\le32$, and a maximum shift in the mean values of a joint distribution of an amplitude--tilt model of 0.006$σ$. This compression reduces the computational cost of a single likelihood evaluation by a factor of 5, from 38 to 7.5 CPU seconds, and it also results in a more robust likelihood by implicitly regularizing nearly degenerate modes. Finally, we use the same compression framework to formulate a numerically stable and computationally efficient variation of the Quadratic Maximum Likelihood implementation that requires less than 3 GB of memory and 2 CPU minutes per iteration for $\ell \le 32$, rendering low-$\ell$ QML CMB power spectrum analysis fully tractable on a standard laptop. △ Less

Submitted 10 September, 2015; v1 submitted 13 June, 2015; originally announced June 2015.

Comments: 13 pages, 13 figures, accepted by ApJS

A CMB Gibbs sampler for localized secondary anisotropies

Authors: Philip Bull, Ingunn K. Wehus, Hans Kristian Eriksen, Pedro G. Ferreira, Unni Fuskeland, Krzysztof M. Gorski, Jeffrey B. Jewell

Abstract: As well as primary fluctuations, CMB temperature maps contain a wealth of additional information in the form of secondary anisotropies. Secondary effects that can be identified with individual objects, such as the thermal and kinetic Sunyaev-Zel'dovich (SZ) effects due to galaxy clusters, are difficult to unambiguously disentangle from foreground contamination and the primary CMB however. We devel… ▽ More As well as primary fluctuations, CMB temperature maps contain a wealth of additional information in the form of secondary anisotropies. Secondary effects that can be identified with individual objects, such as the thermal and kinetic Sunyaev-Zel'dovich (SZ) effects due to galaxy clusters, are difficult to unambiguously disentangle from foreground contamination and the primary CMB however. We develop a Bayesian formalism for rigorously characterising anisotropies that are localised on the sky, taking the TSZ and KSZ effects as an example. Using a Gibbs sampling scheme, we are able to efficiently sample from the joint posterior distribution for a multi-component model of the sky with many thousands of correlated physical parameters. The posterior can then be exactly marginalised to estimate properties of the secondary anisotropies, fully taking into account degeneracies with the other signals in the CMB map. We show that this method is computationally tractable using a simple implementation based on the existing Commander component separation code, and also discuss how other types of secondary anisotropy can be accommodated within our framework. △ Less

Submitted 15 October, 2015; v1 submitted 9 October, 2014; originally announced October 2014.

Comments: Updated to published version. Sections added on instrumental properties and effects of IR sources

Journal ref: ApJS, 219, 10 (2015)

A multi-level solver for Gaussian constrained CMB realizations

Authors: D. S. Seljebotn, K. -A. Mardal, J. B. Jewell, H. K. Eriksen, P. Bull

Abstract: We present a multi-level solver for drawing constrained Gaussian realizations or finding the maximum likelihood estimate of the CMB sky, given noisy sky maps with partial sky coverage. The method converges substantially faster than existing Conjugate Gradient (CG) methods for the same problem. For instance, for the 143 GHz Planck frequency channel, only 3 multi-level W-cycles result in an absolute… ▽ More We present a multi-level solver for drawing constrained Gaussian realizations or finding the maximum likelihood estimate of the CMB sky, given noisy sky maps with partial sky coverage. The method converges substantially faster than existing Conjugate Gradient (CG) methods for the same problem. For instance, for the 143 GHz Planck frequency channel, only 3 multi-level W-cycles result in an absolute error smaller than 1 microKelvin in any pixel. Using 16 CPU cores, this translates to a computational expense of 6 minutes wall time per realization, plus 8 minutes wall time for a power spectrum-dependent precomputation. Each additional W-cycle reduces the error by more than an order of magnitude, at an additional computational cost of 2 minutes. For comparison, we have never been able to achieve similar absolute convergence with conventional CG methods for this high signal-to-noise data set, even after thousands of CG iterations and employing expensive preconditioners. The solver is part of the Commander 2 code, which is available with an open source license at http://commander.bitbucket.org/. △ Less

Submitted 24 February, 2014; v1 submitted 24 August, 2013; originally announced August 2013.

Comments: Updated to match published version (no major changes)

Journal ref: ApJS 210, 24 (2014)

CMB likelihood approximation for banded probability distributions

Authors: E. Gjerløw, K. Mikkelsen, H. K. Eriksen, K. M. Górski, G. Huey, J. B. Jewell, S. K. Næss, G. Rocha, D. S. Seljebotn, I. K. Wehus

Abstract: We investigate sets of random variables that can be arranged sequentially such that a given variable only depends conditionally on its immediate predecessor. For such sets, we show that the full joint probability distribution may be expressed exclusively in terms of uni- and bivariate marginals. Under the assumption that the CMB power spectrum likelihood only exhibits correlations within a banded… ▽ More We investigate sets of random variables that can be arranged sequentially such that a given variable only depends conditionally on its immediate predecessor. For such sets, we show that the full joint probability distribution may be expressed exclusively in terms of uni- and bivariate marginals. Under the assumption that the CMB power spectrum likelihood only exhibits correlations within a banded multipole range, Δl, we apply this expression to two outstanding problems in CMB likelihood analysis. First, we derive a statistically well-defined hybrid likelihood estimator, merging two independent (e.g., low- and high-l) likelihoods into a single expression that properly accounts for correlations between the two. Applying this expression to the WMAP likelihood, we verify that the effect of correlations on cosmological parameters in the transition region is negligible in terms of cosmological parameters for WMAP; the largest relative shift seen for any parameter is 0.06σ. However, because this may not hold for other experimental setups (e.g., for different instrumental noise properties or analysis masks), but must rather be verified on a case-by-case basis, we recommend our new hybridization scheme for future experiments for statistical self-consistency reasons. Second, we use the same expression to improve the convergence rate of the Blackwell-Rao likelihood estimator, reducing the required number of Monte Carlo samples by several orders of magnitude, and thereby extend it to high-l applications. △ Less

Submitted 1 April, 2013; originally announced April 2013.

Comments: 9 pages, 7 figures, 1 table, submitted to ApJ

Analysis of WMAP 7-year Temperature Data: Astrophysics of the Galactic Haze

Authors: Davide Pietrobon, Krzysztof M. Gorski, James Bartlett, Anthony J. Banday, Gregory Dobler, Loris P. L. Colombo, Luca Pagano, Graca Rocha, Rajib Saha, Jeffrey B. Jewell, Sergi R. Hildebrandt, Hans Kristian Eriksen, Charles R. Lawrence

Abstract: We analyse WMAP 7-year temperature data, jointly modeling the cosmic microwave background (CMB) and Galactic foreground emission. We use the Commander code based on Gibbs sampling. Thus, from the WMAP7 data, we derive simultaneously the CMB and Galactic components on scales larger than 1deg with sensitivity improved relative to previous work. We conduct a detailed study of the low-frequency foregr… ▽ More We analyse WMAP 7-year temperature data, jointly modeling the cosmic microwave background (CMB) and Galactic foreground emission. We use the Commander code based on Gibbs sampling. Thus, from the WMAP7 data, we derive simultaneously the CMB and Galactic components on scales larger than 1deg with sensitivity improved relative to previous work. We conduct a detailed study of the low-frequency foreground with particular focus on the "microwave haze" emission around the Galactic center. We demonstrate improved performance in quantifying the diffuse galactic emission when Haslam 408MHz data are included together with WMAP7, and the spinning and thermal dust emission is modeled jointly. We also address the question of whether the hypothetical galactic haze can be explained by a spatial variation of the synchrotron spectral index. The excess of emission around the Galactic center appears stable with respect to variations of the foreground model that we study. Our results demonstrate that the new galactic foreground component - the microwave haze - is indeed present. △ Less

Submitted 30 July, 2012; v1 submitted 25 October, 2011; originally announced October 2011.

Comments: 16 pages, 16 figures, Published on ApJ

Journal ref: ApJ, 755, 69 (2012)

Bayesian component separation and CMB estimation for the 5-year WMAP temperature data

Authors: C. Dickinson, H. K. Eriksen, A. J. Banday, J. B. Jewell, K. M. Gorski, G. Huey, C. R. Lawrence, I. J. O'Dwyer, B. D. Wandelt

Abstract: A well-tested and validated Gibbs sampling code, that performs component separation and CMB power spectrum estimation, was applied to the {\it WMAP} 5-yr data. Using a simple model consisting of CMB, noise, monopoles and dipoles, a ``per pixel'' low-frequency power-law (fitting for both amplitude and spectral index), and a thermal dust template with fixed spectral index, we found that the low-… ▽ More A well-tested and validated Gibbs sampling code, that performs component separation and CMB power spectrum estimation, was applied to the {\it WMAP} 5-yr data. Using a simple model consisting of CMB, noise, monopoles and dipoles, a ``per pixel'' low-frequency power-law (fitting for both amplitude and spectral index), and a thermal dust template with fixed spectral index, we found that the low-$\ell$ ($\ell < 50$) CMB power spectrum is in good agreement with the published {\it WMAP}5 results. Residual monopoles and dipoles were found to be small ($\lesssim 3 μ$K) or negligible in the 5-yr data. We comprehensively tested the assumptions that were made about the foregrounds (e.g. dust spectral index, power-law spectral index prior, templates), and found that the CMB power spectrum was insensitive to these choices. We confirm the asymmetry of power between the north and south ecliptic hemispheres, which appears to be robust against foreground modeling. The map of low frequency spectral indices indicates a steeper spectrum on average ($β=-2.97\pm0.21$) relative to those found at low ($\sim$GHz) frequencies. △ Less

Submitted 25 March, 2009; originally announced March 2009.

Comments: 18 pages, 19 figures, submitted to ApJ

Journal ref: Astrophys.J.705:1607-1623,2009

CMB likelihood approximation by a Gaussianized Blackwell-Rao estimator

Authors: Ø. Rudjord, N. E. Groeneboom, H. K. Eriksen, Greg Huey, K. M. Górski, J. B. Jewell

Abstract: We introduce a new CMB temperature likelihood approximation called the Gaussianized Blackwell-Rao (GBR) estimator. This estimator is derived by transforming the observed marginal power spectrum distributions obtained by the CMB Gibbs sampler into standard univariate Gaussians, and then approximate their joint transformed distribution by a multivariate Gaussian. The method is exact for full-sky c… ▽ More We introduce a new CMB temperature likelihood approximation called the Gaussianized Blackwell-Rao (GBR) estimator. This estimator is derived by transforming the observed marginal power spectrum distributions obtained by the CMB Gibbs sampler into standard univariate Gaussians, and then approximate their joint transformed distribution by a multivariate Gaussian. The method is exact for full-sky coverage and uniform noise, and an excellent approximation for sky cuts and scanning patterns relevant for modern satellite experiments such as WMAP and Planck. A single evaluation of this estimator between l=2 and 200 takes ~0.2 CPU milliseconds, while for comparison, a single pixel space likelihood evaluation between l=2 and 30 for a map with ~2500 pixels requires ~20 seconds. We apply this tool to the 5-year WMAP temperature data, and re-estimate the angular temperature power spectrum, $C_{\ell}$, and likelihood, L(C_l), for l<=200, and derive new cosmological parameters for the standard six-parameter LambdaCDM model. Our spectrum is in excellent agreement with the official WMAP spectrum, but we find slight differences in the derived cosmological parameters. Most importantly, the spectral index of scalar perturbations is n_s=0.973 +/- 0.014, 1.9 sigma away from unity and 0.6 sigma higher than the official WMAP result, n_s = 0.965 +/- 0.014. This suggests that an exact likelihood treatment is required to higher l's than previously believed, reinforcing and extending our conclusions from the 3-year WMAP analysis. In that case, we found that the sub-optimal likelihood approximation adopted between l=12 and 30 by the WMAP team biased n_s low by 0.4 sigma, while here we find that the same approximation between l=30 and 200 introduces a bias of 0.6 sigma in n_s. △ Less

Submitted 26 September, 2008; originally announced September 2008.

Comments: 10 pages, 7 figures, submitted to ApJ

Journal ref: Astrophys.J.692:1669-1677,2009

A Markov Chain Monte Carlo Algorithm for analysis of low signal-to-noise CMB data

Authors: J. B. Jewell, H. K. Eriksen, B. D. Wandelt, I. J. O'Dwyer, G. Huey, K. M. Gorski

Abstract: We present a new Monte Carlo Markov Chain algorithm for CMB analysis in the low signal-to-noise regime. This method builds on and complements the previously described CMB Gibbs sampler, and effectively solves the low signal-to-noise inefficiency problem of the direct Gibbs sampler. The new algorithm is a simple Metropolis-Hastings sampler with a general proposal rule for the power spectrum, C_l,… ▽ More We present a new Monte Carlo Markov Chain algorithm for CMB analysis in the low signal-to-noise regime. This method builds on and complements the previously described CMB Gibbs sampler, and effectively solves the low signal-to-noise inefficiency problem of the direct Gibbs sampler. The new algorithm is a simple Metropolis-Hastings sampler with a general proposal rule for the power spectrum, C_l, followed by a particular deterministic rescaling operation of the sky signal. The acceptance probability for this joint move depends on the sky map only through the difference of chi-squared between the original and proposed sky sample, which is close to unity in the low signal-to-noise regime. The algorithm is completed by alternating this move with a standard Gibbs move. Together, these two proposals constitute a computationally efficient algorithm for mapping out the full joint CMB posterior, both in the high and low signal-to-noise regimes. △ Less

Submitted 3 July, 2008; originally announced July 2008.

Comments: Submitted to ApJ

Journal ref: Astrophys.J.697:258-268,2009

Joint Bayesian component separation and CMB power spectrum estimation

Authors: H. K. Eriksen, J. B. Jewell, C. Dickinson, A. J. Banday, K. M. Gorski, C. R. Lawrence

Abstract: We describe and implement an exact, flexible, and computationally efficient algorithm for joint component separation and CMB power spectrum estimation, building on a Gibbs sampling framework. Two essential new features are 1) conditional sampling of foreground spectral parameters, and 2) joint sampling of all amplitude-type degrees of freedom (e.g., CMB, foreground pixel amplitudes, and global t… ▽ More We describe and implement an exact, flexible, and computationally efficient algorithm for joint component separation and CMB power spectrum estimation, building on a Gibbs sampling framework. Two essential new features are 1) conditional sampling of foreground spectral parameters, and 2) joint sampling of all amplitude-type degrees of freedom (e.g., CMB, foreground pixel amplitudes, and global template amplitudes) given spectral parameters. Given a parametric model of the foreground signals, we estimate efficiently and accurately the exact joint foreground-CMB posterior distribution, and therefore all marginal distributions such as the CMB power spectrum or foreground spectral index posteriors. The main limitation of the current implementation is the requirement of identical beam responses at all frequencies, which restricts the analysis to the lowest resolution of a given experiment. We outline a future generalization to multi-resolution observations. To verify the method, we analyse simple models and compare the results to analytical predictions. We then analyze a realistic simulation with properties similar to the 3-yr WMAP data, downgraded to a common resolution of 3 degree FWHM. The results from the actual 3-yr WMAP temperature analysis are presented in a companion Letter. △ Less

Submitted 5 November, 2007; v1 submitted 7 September, 2007; originally announced September 2007.

Comments: 23 pages, 16 figures; version accepted for publication in ApJ -- only minor changes, all clarifications. More information about the WMAP3 analysis available at http://www.astro.uio.no/~hke under the Research tab

Journal ref: Astrophys.J.676:10-32,2008

The joint large-scale foreground-CMB posteriors of the 3-year WMAP data

Authors: H. K. Eriksen, C. Dickinson, J. B. Jewell, A. J. Banday, K. M. Gorski, C. R. Lawrence

Abstract: Using a Gibbs sampling algorithm for joint CMB estimation and component separation, we compute the large-scale CMB and foreground posteriors of the 3-yr WMAP temperature data. Our parametric data model includes the cosmological CMB signal and instrumental noise, a single power law foreground component with free amplitude and spectral index for each pixel, a thermal dust template with a single fr… ▽ More Using a Gibbs sampling algorithm for joint CMB estimation and component separation, we compute the large-scale CMB and foreground posteriors of the 3-yr WMAP temperature data. Our parametric data model includes the cosmological CMB signal and instrumental noise, a single power law foreground component with free amplitude and spectral index for each pixel, a thermal dust template with a single free overall amplitude, and free monopoles and dipoles at each frequency. This simple model yields a surprisingly good fit to the data over the full frequency range from 23 to 94 GHz. We obtain a new estimate of the CMB sky signal and power spectrum, and a new foreground model, including a measurement of the effective spectral index over the high-latitude sky. A particularly significant result is the detection of a common spurious offset in all frequency bands of ~ -13muK, as well as a dipole in the V-band data. Correcting for these is essential when determining the effective spectral index of the foregrounds. We find that our new foreground model is in good agreement with template-based model presented by the WMAP team, but not with their MEM reconstruction. We believe the latter may be at least partially compromised by the residual offsets and dipoles in the data. Fortunately, the CMB power spectrum is not significantly affected by these issues, as our new spectrum is in excellent agreement with that published by the WMAP team. The corresponding cosmological parameters are also virtually unchanged. △ Less

Submitted 7 September, 2007; originally announced September 2007.

Comments: 5 pages, 4 figures, submitted to ApJL. Background data are available at http://www.astro.uio.no/~hke under the Research tab

Journal ref: Astrophys.J.672:L87-L90,2008

Bayesian analysis of the low-resolution polarized 3-year WMAP sky maps

Authors: H. K. Eriksen, Greg Huey, A. J. Banday, K. M. Gorski, J. B. Jewell, I. J. O'Dwyer, B. D. Wandelt

Abstract: We apply a previously developed Gibbs sampling framework to the foreground corrected 3-yr WMAP polarization data and compute the power spectrum and residual foreground template amplitude posterior distributions. We first analyze the co-added Q- and V-band data, and compare our results to the likelihood code published by the WMAP team. We find good agreement, and thus verify the numerics and data… ▽ More We apply a previously developed Gibbs sampling framework to the foreground corrected 3-yr WMAP polarization data and compute the power spectrum and residual foreground template amplitude posterior distributions. We first analyze the co-added Q- and V-band data, and compare our results to the likelihood code published by the WMAP team. We find good agreement, and thus verify the numerics and data processing steps of both approaches. However, we also analyze the Q- and V-bands separately, allowing for non-zero EB cross-correlations and including two individual foreground template amplitudes tracing synchrotron and dust emission. In these analyses, we find tentative evidence of systematics: The foreground tracers correlate with each of the Q- and V-band sky maps individually, although not with the co-added QV map; there is a noticeable negative EB cross-correlation at l <~ 16 in the V-band map; and finally, when relaxing the constraints on EB and BB, noticeable differences are observed between the marginalized band powers in the Q- and V-bands. Further studies of these features are imperative, given the importance of the low-l EE spectrum on the optical depth of reionization tau and the spectral index of scalar perturbations n_s. △ Less

Submitted 24 May, 2007; originally announced May 2007.

Comments: 5 pages, 4 figures, submitted to ApJL

Journal ref: Astrophys.J.665:L1-L4,2007

Estimation of Polarized Power Spectra by Gibbs sampling

Authors: D. L. Larson, H. K. Eriksen, B. D. Wandelt, K. M. Gorski, Greg Huey, J. B. Jewell, I. J. O'Dwyer

Abstract: Earlier papers introduced a method of accurately estimating the angular cosmic microwave background (CMB) temperature power spectrum based on Gibbs sampling. Here we extend this framework to polarized data. All advantages of the Gibbs sampler still apply, and exact analysis of mega-pixel polarized data sets is thus feasible. These advantages may be even more important for polarization measuremen… ▽ More Earlier papers introduced a method of accurately estimating the angular cosmic microwave background (CMB) temperature power spectrum based on Gibbs sampling. Here we extend this framework to polarized data. All advantages of the Gibbs sampler still apply, and exact analysis of mega-pixel polarized data sets is thus feasible. These advantages may be even more important for polarization measurements than for temperature measurements. While approximate methods can alias power from the larger E-mode spectrum into the weaker B-mode spectrum, the Gibbs sampler (or equivalently, exact likelihood evaluations) allows for a statistically optimal separation of these modes in terms of power spectra. To demonstrate the method, we analyze two simulated data sets: 1) a hypothetical future CMBPol mission, with the focus on B-mode estimation; and 2) a Planck-like mission, to highlight the computational feasibility of the method. △ Less

Submitted 2 October, 2006; v1 submitted 31 July, 2006; originally announced August 2006.

Comments: 8 pages, 5 figures. High-resolution version available from http://www.astro.uio.no/~hke/docs/larson_et_al_2006.ps.gz; accepted for publication in ApJ

Journal ref: Astrophys.J.656:653-660,2007

A re-analysis of the three-year WMAP temperature power spectrum and likelihood

Authors: H. K. Eriksen, Greg Huey, R. Saha, F. K. Hansen, J. Dick, A. J. Banday, K. M. Gorski, P. Jain, J. B. Jewell, L. Knox, D. L. Larson, I. J. O'Dwyer, T. Souradeep, B. D. Wandelt

Abstract: We analyze the three-year WMAP temperature anisotropy data seeking to confirm the power spectrum and likelihoods published by the WMAP team. We apply five independent implementations of four algorithms to the power spectrum estimation and two implementations to the parameter estimation. Our single most important result is that we broadly confirm the WMAP power spectrum and analysis. Still, we do… ▽ More We analyze the three-year WMAP temperature anisotropy data seeking to confirm the power spectrum and likelihoods published by the WMAP team. We apply five independent implementations of four algorithms to the power spectrum estimation and two implementations to the parameter estimation. Our single most important result is that we broadly confirm the WMAP power spectrum and analysis. Still, we do find two small but potentially important discrepancies: On large angular scales there is a small power excess in the WMAP spectrum (5-10% at l<~30) primarily due to likelihood approximation issues between 13 <= l <~30. On small angular scales there is a systematic difference between the V- and W-band spectra (few percent at l>~300). Recently, the latter discrepancy was explained by Huffenberger et al. (2006) in terms of over-subtraction of unresolved point sources. As far as the low-l bias is concerned, most parameters are affected by a few tenths of a sigma. The most important effect is seen in n_s. For the combination of WMAP, Acbar and BOOMERanG, the significance of n_s =/ 1 drops from ~2.7 sigma to ~2.3 sigma when correcting for this bias. We propose a few simple improvements to the low-l WMAP likelihood code, and introduce two important extensions to the Gibbs sampling method that allows for proper sampling of the low signal-to-noise regime. Finally, we make the products from the Gibbs sampling analysis publically available, thereby providing a fast and simple route to the exact likelihood without the need of expensive matrix inversions. △ Less

Submitted 2 October, 2006; v1 submitted 5 June, 2006; originally announced June 2006.

Comments: 14 pages, 7 figures. Accepted for publication in ApJ. Numerical results unchanged, but interpretation sharpened: Likelihood approximation issues at l=13-30 far more important than potential foreground issues at l <= 12. Gibbs products (spectrum and sky samples, and "easy-to-use" likelihood module) available from http://www.astro.uio.no/~hke/ under "Research"

Journal ref: Astrophys.J.656:641-652,2007

Cosmological Parameter Constraints as Derived from the Wilkinson Microwave Anisotropy Probe Data via Gibbs Sampling and the Blackwell-Rao Estimator

Authors: M. Chu, H. K. Eriksen, L. Knox, K. M. Gorski, J. B. Jewell, D. L. Larson, I. J. O'Dwyer, B. D. Wandelt

Abstract: We study the Blackwell-Rao (BR) estimator of the probability distribution of the angular power spectrum, P(C_l|d), by applying it to samples of full-sky no-noise CMB maps generated via Gibbs sampling. We find the estimator, given a set of samples, to be very fast and also highly accurate, as determined by tests with simulated data. We also find that the number of samples required for convergence… ▽ More We study the Blackwell-Rao (BR) estimator of the probability distribution of the angular power spectrum, P(C_l|d), by applying it to samples of full-sky no-noise CMB maps generated via Gibbs sampling. We find the estimator, given a set of samples, to be very fast and also highly accurate, as determined by tests with simulated data. We also find that the number of samples required for convergence of the BR estimate rises rapidly with increasing l, at least at low l. Our existing sample chains are only long enough to achieve convergence at l less than about 40. In comparison with P(C_l|d) as reported by the WMAP team we find significant differences at these low l values. These differences lead to up to approximately 0.5 sigma shifts in the estimates of parameters in a 7-parameter Lambda CDM model with non-zero dn_s/dln(k). Fixing dn_s/dln(k)= 0 makes these shifts much less significant. Unlike existing analytic approximations, the BR estimator can be straightforwardly extended for the case of power spectra from correlated fields, such as temperature and polarization. We discuss challenges to extending the procedure to higher l and provide some solutions. △ Less

Submitted 26 November, 2004; originally announced November 2004.

Comments: 11 pages, 8 figures, to be submitted to Phys. Rev. D

Journal ref: Phys.Rev.D71:103002,2005

Power spectrum estimation from high-resolution maps by Gibbs sampling

Authors: H. K. Eriksen, I. J. O'Dwyer, J. B. Jewell, B. D. Wandelt, D. L. Larson, K. M. Gorski, S. Levin, A. J. Banday, P. B. Lilje

Abstract: We revisit a recently introduced power spectrum estimation technique based on Gibbs sampling, with the goal of applying it to the high-resolution WMAP data. In order to facilitate this analysis, a number of sophistications have to be introduced, each of which is discussed in detail. We have implemented two independent versions of the algorithm to cross-check the computer codes, and to verify tha… ▽ More We revisit a recently introduced power spectrum estimation technique based on Gibbs sampling, with the goal of applying it to the high-resolution WMAP data. In order to facilitate this analysis, a number of sophistications have to be introduced, each of which is discussed in detail. We have implemented two independent versions of the algorithm to cross-check the computer codes, and to verify that a particular solution to any given problem does not affect the scientific results. We then apply these programs to simulated data with known properties at intermediate (N_side = 128) and high (N_side = 512) resolutions, to study effects such as incomplete sky coverage and white vs. correlated noise. From these simulations we also establish the Markov chain correlation length as a function of signal-to-noise ratio, and give a few comments on the properties of the correlation matrices involved. Parallelization issues are also discussed, with emphasis on real-world limitations imposed by current super-computer facilities. The scientific results from the analysis of the first-year WMAP data are presented in a companion letter. △ Less

Submitted 23 August, 2004; v1 submitted 1 July, 2004; originally announced July 2004.

Comments: 15 pages, 11 figures, accepted for publication in ApJS

Journal ref: Astrophys.J.Suppl. 155 (2004) 227-241

Bayesian Power Spectrum Analysis of the First-Year WMAP data

Authors: I. J. O'Dwyer, H. K. Eriksen, B. D. Wandelt, J. B. Jewell, D. L. Larson, K. M. Gorski, A. J. Banday, S. Levin, P. B. Lilje

Abstract: We present the first results from a Bayesian analysis of the WMAP first year data using a Gibbs sampling technique. Using two independent, parallel supercomputer codes we analyze the WMAP Q, V and W bands. The analysis results in a full probabilistic description of the information the WMAP data set contains about the power spectrum and the all-sky map of the cosmic microwave background anisotrop… ▽ More We present the first results from a Bayesian analysis of the WMAP first year data using a Gibbs sampling technique. Using two independent, parallel supercomputer codes we analyze the WMAP Q, V and W bands. The analysis results in a full probabilistic description of the information the WMAP data set contains about the power spectrum and the all-sky map of the cosmic microwave background anisotropies. We present the complete probability distributions for each C_l including any non-Gaussianities of the power spectrum likelihood. While we find good overall agreement with the previously published WMAP spectrum, our analysis uncovers discrepancies in the power spectrum estimates at low l multipoles. For example we claim the best-fit Lambda-CDM model is consistent with the C_2 inferred from our combined Q+V+W analysis with a 10% probability of an even larger theoretical C_2. Based on our exact analysis we can therefore attribute the "low quadrupole issue" to a statistical fluctuation. △ Less

Submitted 7 December, 2004; v1 submitted 1 July, 2004; originally announced July 2004.

Comments: 5 pages. 4 figures. For additional information and data see http://www.astro.uiuc.edu/~iodwyer/research#wmap

Journal ref: Astrophys.J. 617 (2004) L99-102