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.
Planck intermediate results Aghanim, N; Ashdown, M; Aumont, J ...
Astronomy and astrophysics (Berlin),
12/2016, Letnik:
596
Journal Article
Recenzirano
Odprti dostop
Parity-violating extensions of the standard electromagnetic theory cause in vacuo rotation of the plane of polarization of propagating photons. This effect, also known as cosmic birefringence, has an ...impact on the cosmic microwave background (CMB) anisotropy angular power spectra, producing non-vanishing T-B and E-B correlations that are otherwise null when parity is a symmetry. Here we present new constraints on an isotropic rotation, parametrized by the angle alpha, derived from Planck 2015 CMB polarization data. To increase the robustness of our analyses, we employ two complementary approaches, in harmonic space and in map space, the latter based on a peak stacking technique. The two approaches provide estimates for alpha that are in agreement within statistical uncertainties and are very stable against several consistency tests.Considering the T-B and E-B information jointly, we find alpha=0degrees.31 + or - 0degrees.05 (stat.) + or - 0degrees.28 (syst.) from alpha = 0degrees.35 + or - 0degrees.05 (stat.) + or - 0degrees.28 (syst.) from the stacking approach. These constraints are compatible with no parity violation and are dominated by the systematic uncertainty in the orientation of Planck's polarization-sensitive bolometers.
The two-point correlation function of the galaxy distribution is a key cosmological observable that allows us to constrain the dynamical and geometrical state of our Universe. To measure the ...correlation function we need to know both the galaxy positions and the expected galaxy density field. The expected field is commonly specified using a Monte-Carlo sampling of the volume covered by the survey and, to minimize additional sampling errors, this random catalog has to be much larger than the data catalog. Correlation function estimators compare data–data pair counts to data–random and random–random pair counts, where random–random pairs usually dominate the computational cost. Future redshift surveys will deliver spectroscopic catalogs of tens of millions of galaxies. Given the large number of random objects required to guarantee sub-percent accuracy, it is of paramount importance to improve the efficiency of the algorithm without degrading its precision. We show both analytically and numerically that splitting the random catalog into a number of subcatalogs of the same size as the data catalog when calculating random–random pairs and excluding pairs across different subcatalogs provides the optimal error at fixed computational cost. For a random catalog fifty times larger than the data catalog, this reduces the computation time by a factor of more than ten without affecting estimator variance or bias.
We present a method for beam-deconvolving cosmic microwave background (CMB) anisotropy measurements. The code takes as input the time-ordered data along with the corresponding detector pointings and ...known beam shapes, and produces as output the harmonic aTlm, aElm, and aBlm coefficients of the observed sky. From these one can derive temperature and Q and U polarisation maps. The method is applicable to absolute CMB measurements with wide sky coverage, and is independent of the scanning strategy. We tested the code with extensive simulations, mimicking the resolution and data volume of Planck 30 GHz and 70 GHz channels, but with exaggerated beam asymmetry. We applied it to multipoles up to l = 1700 and examined the results in both pixel space and harmonic space. We also tested the method in presence of white noise.
BEYONDPLANCK Svalheim, T. L.; Andersen, K. J.; Aurlien, R. ...
Astronomy and astrophysics (Berlin),
06/2023, Letnik:
675
Journal Article
Recenzirano
Odprti dostop
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.
We present two novel methods for the estimation of the angular power spectrum of cosmic microwave background (CMB) anisotropies. We assume an absolute CMB experiment with arbitrary asymmetric beams ...and arbitrary sky coverage. The methods differ from the earlier ones in that the power spectrum is estimated directly from the time-ordered data, without first compressing the data into a sky map, and they take into account the effect of asymmetric beams. In particular, they correct the beam-induced leakage from temperature to polarization. The methods are applicable to a case where part of the sky has been masked out to remove foreground contamination, leaving a pure CMB signal, but incomplete sky coverage. The first method (deconvolution quadratic maximum likelihood) is derived as the optimal quadratic estimator, which simultaneously yields an unbiased spectrum estimate and minimizes its variance. We successfully apply it to multipoles up to l = 200. The second method is derived as a weak-signal approximation from the first one. It yields an unbiased estimate for the full multipole range, but relaxes the requirement of minimal variance. We validate the methods with simulations for the 70 GHz channel of Planck surveyor, and demonstrate that we are able to correct the beam effects in the TT, EE, BB and TE spectra up to multipole l = 1500. Together, the two methods cover the complete multipole range with no gap in between.
Beam-deconvolved Planck LFI maps Keihänen, E.; Lindholm, V.; Lopez-Caniego, M. ...
Astronomy and astrophysics (Berlin),
12/2019, Letnik:
632
Journal Article
Recenzirano
Odprti dostop
The Planck Collaboration made its final data release in 2018. In this paper we describe beam-deconvolution map products made from Planck Low Frequency Instrument (LFI) data using the artDeco ...deconvolution code to symmetrize the effective beam. The deconvolution results are auxiliary data products, available through the Planck Legacy Archive. Analysis of these deconvolved survey difference maps reveals signs of residual signal in the 30 GHz and 44 GHz frequency channels. We produce low-resolution maps and corresponding noise covariance matrices (NCVMs). The NCVMs agree reasonably well with the half-ring noise estimates except for 44 GHz, where we observe an asymmetry between EE and BB noise spectra, possibly a sign of further unresolved systematic error. In contrast to the official Planck LFI maps, the beam-deconvolution maps have not been corrected for bandpass mismatch, and the residual noise is not well approximated by white noise.
Late stages of stellar evolution are characterized by copious mass-loss events whose signature is the formation of circumstellar envelopes (CSE). Planck multi-frequency measurements have provided ...relevant information on a sample of Galactic planetary nebulae (PNe) in the important and relatively unexplored observational band between 30 and 857 GHz. Planck enables the assembly of comprehensive PNe spectral energy distributions (SEDs) from radio to far-IR frequencies. Modelling the derived SEDs provides us with information on physical properties of CSEs and the mass content of both main components: ionized gas, traced by the free-free emission at cm–mm waves; and thermal dust, traced by the millimetre and far-IR emission. In particular, the amount of ionized gas and dust has been derived here. Such quantities have also been estimated for the very young PN CRL 618, where the strong variability observed in its radio and millimetre emission has previously prevented constructing its SED. A morphological study of the Helix Nebula was also performed. Planck maps reveal, for the first time, the spatial distribution of the dust inside the envelope, allowing us to identify different components, the most interesting of which is a very extended component (up to 1 pc) that may be related to a region where the slow expanding envelope is interacting with the surrounding interstellar medium.
BEYONDPLANCK Colombo, L. P. L.; Eskilt, J. R.; Paradiso, S. ...
Astronomy and astrophysics (Berlin),
06/2023, Letnik:
675
Journal Article
Recenzirano
Odprti dostop
We present posterior sample-based cosmic microwave background (CMB) constraints from
Planck
LFI and WMAP observations as derived through global end-to-end Bayesian processing within the B
EYOND
P
...LANCK
framework. We first used these samples to study correlations between CMB, foreground, and instrumental parameters. We identified a particularly strong degeneracy between CMB temperature fluctuations and free-free emission on intermediate angular scales (400 ≲
ℓ
≲ 600), mitigated through model reduction, masking, and resampling. We compared our posterior-based CMB results with previous
Planck
products and found a generally good agreement, however, with notably higher noise due to our exclusion of
Planck
HFI data. We found a best-fit CMB dipole amplitude of 3362.7 ± 1.4 μK, which is in excellent agreement with previous
Planck
results. The quoted dipole uncertainty is derived directly from the sampled posterior distribution and does not involve any ad hoc contributions for
Planck
instrumental systematic effects. Similarly, we find a temperature quadrupole amplitude of $ \sigma^{TT}_2=229\pm97\,\muup{\rm K}^2 $ , which is in good agreement with previous results in terms of the amplitude, but the uncertainty is one order of magnitude greater than the naive diagonal Fisher uncertainty. Concurrently, we find less evidence of a possible alignment between the quadrupole and octopole than previously reported, due to a much larger scatter in the individual quadrupole coefficients that is caused both by marginalizing over a more complete set of systematic effects – as well as by requiring a more conservative analysis mask to mitigate the free-free degeneracy. For higher multipoles, we find that the angular temperature power spectrum is generally in good agreement with both
Planck
and WMAP. At the same time, we note that this is the first time that the sample-based, asymptotically exact Blackwell-Rao estimator has been successfully established for multipoles up to
ℓ
≤ 600. It now accounts for the majority of the cosmologically important information. Overall, this analysis demonstrates the unique capabilities of the Bayesian approach with respect to end-to-end systematic uncertainty propagation and we believe it can and should play an important role in the analysis of future CMB experiments. Cosmological parameter constraints are presented in a companion paper.
We present an analysis of the effects of beam deconvolution on noise properties in CMB measurements. The analysis is built around the artDeco beam deconvolver code. We derive a low-resolution noise ...covariance matrix that describes the residual noise in deconvolution products, both in harmonic and pixel space. The matrix models the residual correlated noise that remains in time-ordered data after destriping, and the effect of deconvolution on this noise. To validate the results, we generate noise simulations that mimic the data from the Planck LFI instrument. A χ2 test for the full 70 GHz covariance in multipole range ℓ = 0 − 50 yields a mean reduced χ2 of 1.0037. We compare two destriping options, full and independent destriping, when deconvolving subsets of available data. Full destriping leaves substantially less residual noise, but leaves data sets intercorrelated. We also derive a white noise covariance matrix that provides an approximation of the full noise at high multipoles, and study the properties on high-resolution noise in pixel space through simulations.