Your email was sent successfully. Check your inbox.

An error occurred while sending the email. Please try again.

Proceed reservation?

Export
Filter
  • Fuskeland, U.  (18)
Type of Medium
Publisher
Language
Years
Subjects(RVK)
  • 1
    In: Astronomy & Astrophysics, EDP Sciences, Vol. 675 ( 2023-07), p. A4-
    Abstract: End-to-end simulations play a key role in the analysis of any high-sensitivity cosmic microwave background (CMB) experiment, providing high-fidelity systematic error propagation capabilities that are 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 posterior and prior simulations, respectively. We show that the two options lead to significantly different correlation structures, as prior simulations (contrary to posterior simulations) effectively include cosmic variance, but they exclude realization-specific correlations from non-linear degeneracies. Consequently, they quantify fundamentally different types of uncertainties. We argue that as a result, they also have different and complementary scientific uses, even if this dichotomy is not absolute. In particular, posterior simulations are in general more convenient for parameter estimation studies, while prior simulations are generally more convenient for model testing. Before B EYOND P LANCK , most pipelines used a mix of constrained and random inputs and applied the same hybrid simulations for all applications, even though the statistical justification for this is not always evident. B EYOND P LANCK 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 B EYOND P LANCK posterior simulations and their uses are described extensively in a suite of companion papers. In this work, we consider one important applications of the corresponding prior simulations, namely, code validation. Specifically, we generated a set of one-year LFI 30 GHz prior simulations with known inputs and we used these to validate the core low-level B EYOND P LANCK algorithms dealing with gain estimation, correlated noise estimation, and mapmaking.
    Type of Medium: Online Resource
    ISSN: 0004-6361 , 1432-0746
    RVK:
    RVK:
    Language: English
    Publisher: EDP Sciences
    Publication Date: 2023
    detail.hit.zdb_id: 1458466-9
    SSG: 16,12
    Library Location Call Number Volume/Issue/Year Availability
    BibTip Others were also interested in ...
  • 2
    In: Astronomy & Astrophysics, EDP Sciences, Vol. 675 ( 2023-07), p. A9-
    Abstract: We discuss the treatment of bandpass and beam leakage corrections in the Bayesian B EYOND P LANCK cosmic microwave background (CMB) analysis pipeline as applied to the Planck LFI measurements. As a preparatory step, we first applied three corrections to the nominal LFI bandpass profiles, including the removal of a known systematic effect in the ground measuring equipment at 61 GHz, along with a 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 B EYOND P LANCK , 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 ± 0.03 GHz. The relative corrections were fitted using a spurious map approach that is fundamentally similar to the method pioneered by the WMAP team, but excluding the introduction of many additional degrees of freedom. All the bandpass parameters were sampled using a standard Metropolis sampler within the main B EYOND P LANCK Gibbs chain and the bandpass uncertainties were thus propagated to all other data products in the analysis. In summary, we find that our bandpass model significantly reduces leakage effects. For beam leakage corrections, we adopted the official Planck LFI beam estimates without any additional degrees of freedom and we only marginalized over the underlying sky model. We note that this is the first time that leakage from beam mismatch has been included for Planck LFI maps.
    Type of Medium: Online Resource
    ISSN: 0004-6361 , 1432-0746
    RVK:
    RVK:
    Language: English
    Publisher: EDP Sciences
    Publication Date: 2023
    detail.hit.zdb_id: 1458466-9
    SSG: 16,12
    Library Location Call Number Volume/Issue/Year Availability
    BibTip Others were also interested in ...
  • 3
    In: Astronomy & Astrophysics, EDP Sciences, Vol. 675 ( 2023-07), p. A2-
    Abstract: We present a Gibbs sampling solution to the mapmaking problem for cosmic microwave background (CMB) measurements that builds on existing destriping methodology. Gibbs sampling breaks the computationally heavy destriping problem into two separate steps: noise filtering and map binning. Considered as two separate steps, both are computationally much cheaper than solving the combined problem. This provides a huge performance benefit as compared to traditional methods and it allows us, for the first time, to bring the destriping baseline length to a single sample. Here, we applied the Gibbs procedure to simulated Planck 30 GHz data. We find that gaps in the time-ordered data are handled efficiently by filling them in with simulated noise as part of the Gibbs process. The Gibbs procedure yields a chain of map samples, from which we are able to compute the posterior mean as a best-estimate map. The variation in the chain provides information on the correlated residual noise, without the need to construct a full noise covariance matrix. However, if only a single maximum-likelihood frequency map estimate is required, we find that traditional conjugate gradient solvers converge much faster than a Gibbs sampler in terms of the total number of iterations. The conceptual advantages of the Gibbs sampling approach lies in statistically well-defined error propagation and systematic error correction. This methodology thus forms the conceptual basis for the mapmaking algorithm employed in the B EYOND P LANCK framework, which implements the first end-to-end Bayesian analysis pipeline for CMB observations.
    Type of Medium: Online Resource
    ISSN: 0004-6361 , 1432-0746
    RVK:
    RVK:
    Language: English
    Publisher: EDP Sciences
    Publication Date: 2023
    detail.hit.zdb_id: 1458466-9
    SSG: 16,12
    Library Location Call Number Volume/Issue/Year Availability
    BibTip Others were also interested in ...
  • 4
    In: Astronomy & Astrophysics, EDP Sciences, Vol. 675 ( 2023-07), p. A14-
    Abstract: Using the Planck Low Frequency Instrument (LFI) and WMAP data within the global Bayesian B EYOND P LANCK framework, we constrained the polarized foreground emission between 30 and 70 GHz. We combined, for the first time, full-resolution Planck LFI time-ordered data with low-resolution WMAP sky maps at 33, 40, and 61 GHz. The spectral parameters were 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 the synchrotron emission, we approximated the SED as a power-law in frequency and we find that the low signal-to-noise ratio of the current data strongly limits the number of free parameters that can be robustly constrained. We partitioned 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 β s Spur  = −3.17 ± 0.06 and β s Plane  = −3.03 ± 0.07, which is in good agreement with previous results. For the thermal dust emission, we assumed a modified blackbody model and we fit a single power-law index across the full sky. We find β d  = 1.64 ± 0.03, which is slightly steeper than the value reported in Planck HFI data, but still statistically consistent at the 2 σ confidence level.
    Type of Medium: Online Resource
    ISSN: 0004-6361 , 1432-0746
    RVK:
    RVK:
    Language: English
    Publisher: EDP Sciences
    Publication Date: 2023
    detail.hit.zdb_id: 1458466-9
    SSG: 16,12
    Library Location Call Number Volume/Issue/Year Availability
    BibTip Others were also interested in ...
  • 5
    In: The Open Journal of Astrophysics, The Open Journal, Vol. 6 ( 2023-03-16)
    Type of Medium: Online Resource
    ISSN: 2565-6120
    Language: Unknown
    Publisher: The Open Journal
    Publication Date: 2023
    detail.hit.zdb_id: 3001843-2
    Library Location Call Number Volume/Issue/Year Availability
    BibTip Others were also interested in ...
  • 6
    In: Astronomy & Astrophysics, EDP Sciences, Vol. 675 ( 2023-07), p. A13-
    Abstract: 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.
    Type of Medium: Online Resource
    ISSN: 0004-6361 , 1432-0746
    RVK:
    RVK:
    Language: English
    Publisher: EDP Sciences
    Publication Date: 2023
    detail.hit.zdb_id: 1458466-9
    SSG: 16,12
    Library Location Call Number Volume/Issue/Year Availability
    BibTip Others were also interested in ...
  • 7
    In: Astronomy & Astrophysics, EDP Sciences, Vol. 675 ( 2023-07), p. A7-
    Abstract: We present a Bayesian calibration algorithm for cosmic microwave background (CMB) observations as implemented within the global end-to-end B EYOND P LANCK framework and applied to the Planck Low Frequency Instrument (LFI) data. Following the most recent Planck analysis, we decomposed the full time-dependent gain into a sum of three nearly orthogonal components: one absolute calibration term, common to all detectors, one time-independent term that can vary between detectors, and one time-dependent component that was allowed to vary between one-hour pointing periods. Each term was then sampled conditionally on all other parameters in the global signal model through Gibbs sampling. The absolute calibration is sampled using only the orbital dipole as a reference source, while the two relative gain components were sampled using the full sky signal, including the orbital and Solar CMB dipoles, CMB fluctuations, and foreground contributions. We discuss various aspects of the data that influence gain estimation, including the dipole-polarization quadrupole degeneracy and processing masks. Comparing our solution to previous pipelines, we find good agreement in general, with relative deviations of −0.67% (−0.84%) for 30 GHz, 0.12% (−0.04%) for 44 GHz and −0.03% (−0.64%) for 70 GHz, compared to Planck PR4 and Planck 2018, respectively. We note that the B EYOND P LANCK calibration was performed globally, which results in better inter-frequency consistency than previous estimates. Additionally, WMAP observations were used actively in the B EYOND P LANCK analysis, which both breaks internal degeneracies in the Planck data set and results in an overall better agreement with WMAP. Finally, we used a Wiener filtering approach to smoothing the gain estimates. We show that this method avoids artifacts in the correlated noise maps as a result of oversmoothing the gain solution, which is difficult to avoid with methods like boxcar smoothing, as Wiener filtering by construction maintains a balance between data fidelity and prior knowledge. Although our presentation and algorithm are currently oriented toward LFI processing, the general procedure is fully generalizable to other experiments, as long as the Solar dipole signal is available to be used for calibration.
    Type of Medium: Online Resource
    ISSN: 0004-6361 , 1432-0746
    RVK:
    RVK:
    Language: English
    Publisher: EDP Sciences
    Publication Date: 2023
    detail.hit.zdb_id: 1458466-9
    SSG: 16,12
    Library Location Call Number Volume/Issue/Year Availability
    BibTip Others were also interested in ...
  • 8
    In: Astronomy & Astrophysics, EDP Sciences, Vol. 675 ( 2023-07), p. A16-
    Abstract: 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.
    Type of Medium: Online Resource
    ISSN: 0004-6361 , 1432-0746
    RVK:
    RVK:
    Language: English
    Publisher: EDP Sciences
    Publication Date: 2023
    detail.hit.zdb_id: 1458466-9
    SSG: 16,12
    Library Location Call Number Volume/Issue/Year Availability
    BibTip Others were also interested in ...
  • 9
    In: Astronomy & Astrophysics, EDP Sciences, Vol. 675 ( 2023-07), p. A3-
    Abstract: We describe the computational infrastructure for end-to-end Bayesian cosmic microwave background (CMB) analysis implemented by the BeyondPlanck Collaboration. The code is called Commander3 . It 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 re-implement 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, fast Fourier transform optimization, and parallelization and load-balancing. Thirdly, we quantify the CPU and RAM requirements for the current B EYOND P LANCK analysis, finding that a total of 1.5 TB of RAM is required for efficient analysis and that the total cost of a full Gibbs sample for LFI 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.
    Type of Medium: Online Resource
    ISSN: 0004-6361 , 1432-0746
    RVK:
    RVK:
    Language: English
    Publisher: EDP Sciences
    Publication Date: 2023
    detail.hit.zdb_id: 1458466-9
    SSG: 16,12
    Library Location Call Number Volume/Issue/Year Availability
    BibTip Others were also interested in ...
  • 10
    In: Astronomy & Astrophysics, EDP Sciences, Vol. 675 ( 2023-07), p. A1-
    Abstract: 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.
    Type of Medium: Online Resource
    ISSN: 0004-6361 , 1432-0746
    RVK:
    RVK:
    Language: English
    Publisher: EDP Sciences
    Publication Date: 2023
    detail.hit.zdb_id: 1458466-9
    SSG: 16,12
    Library Location Call Number Volume/Issue/Year Availability
    BibTip Others were also interested in ...
Close ⊗
This website uses cookies and the analysis tool Matomo. Further information can be found on the KOBV privacy pages