-
A Foreground Model Independent Bayesian CMB Temperature and Polarization Signal Reconstruction and Cosmological Parameter Estimation over Large Angular Scales
Authors:
Albin Joseph,
Ujjal Purkayastha,
Rajib Saha
Abstract:
Recent CMB observations have resulted in very precise observational data. A robust and reliable CMB reconstruction technique can lead to efficient estimation of the cosmological parameters. We demonstrate the performance of our methodology using simulated temperature and polarization observations using cosmic variance limited future generation PRISM satellite mission. We generate samples from the…
▽ More
Recent CMB observations have resulted in very precise observational data. A robust and reliable CMB reconstruction technique can lead to efficient estimation of the cosmological parameters. We demonstrate the performance of our methodology using simulated temperature and polarization observations using cosmic variance limited future generation PRISM satellite mission. We generate samples from the joint distribution by implementing the CMB inverse covariance weighted internal-linear-combination (ILC) with the Gibbs sampling technique. We use the Python Sky Model (PySM), d4f1s1 to generate the realistic foreground templates. The synchrotron emission is parametrized by a spatially varying spectral index, whereas the thermal dust emission is described as a two-component dust model. We estimate the marginalized densities of CMB signal and theoretical angular power spectrum utilizing the samples from the entire posterior distribution. The best-fit cleaned CMB map and the corresponding angular power spectrum are consistent with the CMB realization and the sky angular power spectrum, implying an efficient foreground minimized reconstruction. The likelihood function estimated by making use of the Blackwell-Rao estimator is used for the estimation of the cosmological parameters. Our methodology can estimate the tensor to scalar ratio $r\ge 0.0075$ for the chosen foreground models and the instrumental noise levels. Our current work demonstrates an analysis pipeline starting from the reliable estimation of CMB signal and its angular power spectrum to the case of cosmological parameter estimation using the foreground model independent Gibbs-ILC method.
△ Less
Submitted 30 January, 2023; v1 submitted 15 September, 2022;
originally announced September 2022.
-
A Partial-Sky Gibbs ILC Approach for the Estimation of CMB Posterior over Large Angular Scales of the Sky
Authors:
Vipin Sudevan,
Ujjal Purkayastha,
Rajib Saha
Abstract:
In this article we present a formalism to incorporate the partial-sky maps to the Gibbs ILC algorithm to estimate the joint posterior density of the Cosmic Microwave Background (CMB) signal and the theoretical CMB angular power spectrum given the observed CMB maps. In order to generate the partial-sky maps we mask all the observed CMB maps provided by WMAP and Planck satellite mission using a Gaus…
▽ More
In this article we present a formalism to incorporate the partial-sky maps to the Gibbs ILC algorithm to estimate the joint posterior density of the Cosmic Microwave Background (CMB) signal and the theoretical CMB angular power spectrum given the observed CMB maps. In order to generate the partial-sky maps we mask all the observed CMB maps provided by WMAP and Planck satellite mission using a Gaussian smoothed mask formed on the basis of thermal dust emissions in Planck 353 GHz map. The central galactic region from all the input maps is removed after the application of the smoothed mask. While implementing the Gibbs ILC method on the partial-sky maps, we convert the partial-sky cleaned angular power spectrum to the full-sky angular power spectrum using the mode-mode coupling matrix estimated from the smoothed mask. The main products of our analysis are partial-sky cleaned best-fit CMB map and an estimate of the underlying full-sky theoretical CMB angular power spectrum along with their error estimates. We validate the methodology by performing detailed Monte Carlo simulations after using realistic models of foregrounds and detector noise consistent with the WMAP and Planck frequency channels used in our analysis. We can estimate the posterior density and full-sky theoretical CMB angular power spectrum, without any need to explicitly model the foreground components, from partial-sky maps using our method. Another important feature of this method is that the power spectrum results along with the error estimates can be directly used for cosmological parameter estimations.
△ Less
Submitted 15 October, 2021;
originally announced October 2021.
-
A Foreground Model Independent Estimation of Joint Posterior of CMB E mode Polarization over Large Angular Scales
Authors:
Ujjal Purkayastha,
Vipin Sudevan,
Rajib Saha
Abstract:
Ever since Cosmic Microwave Background (CMB) signal is being measured by various satellites based observations with increasing experimental accuracies there has been a parallel increase in the demand for a CMB reconstruction technique which can provide accurate estimates of CMB signal and the theoretical angular power spectrum along with reliable statistical error estimates associated with them. I…
▽ More
Ever since Cosmic Microwave Background (CMB) signal is being measured by various satellites based observations with increasing experimental accuracies there has been a parallel increase in the demand for a CMB reconstruction technique which can provide accurate estimates of CMB signal and the theoretical angular power spectrum along with reliable statistical error estimates associated with them. In this work, we estimate the joint posterior of CMB E mode signal (S) and corresponding theoretical angular power spectrum (C^E_l) over large angular scales given the simulated polarization observations of future generation COrE satellite mission. To generate samples from the joint distribution we employ the ILC technique with prior information of CMB E mode covariance matrix augmented by a Gibbs sampling technique. We estimate the marginalized densities of S and C^E_l using the samples from full-posterior. The best fit cleaned E mode map and the corresponding angular power spectrum agree well with the input E mode map and the sky power spectrum implying accurate reconstruction using COrE like observations. Using the samples C^E_l of all Gibbs chains we estimate the likelihood function P(C^E_l|D) of any arbitrary C^E_l given simulated observed maps (D) of COrE mission following Blackwell-Rao estimator. The likelihood function can be seamlessly integrated to the cosmological parameter estimation method. Apart from producing an accurate estimate of E mode signal over large angular scales our method also builds a connection between the component reconstruction and reliable cosmological parameter estimation using CMB E mode observations over large angular scales. The entire method does not assume any explicit models for E mode foreground components in order to remove them, which is an attractive property since foreground modelling uncertainty does not pose as a challenge in this case.
△ Less
Submitted 18 April, 2020;
originally announced April 2020.
-
An Accurate Reconstruction of CMB E Mode Signal over Large Angular Scales using Prior Information of CMB Covariance Matrix in ILC Algorithm
Authors:
Ujjal Purkayastha,
Vipin Sudevan,
Rajib Saha
Abstract:
In the recent years, the internal-linear-combination (ILC) method was investigated extensively in the context of reconstruction of Cosmic Microwave Background (CMB) temperature anisotropy signal using observations obtained by WMAP and Planck satellite missions. In this article, we, for the first time, apply the ILC method to reconstruct the large scale CMB E mode polarization signal, which serves…
▽ More
In the recent years, the internal-linear-combination (ILC) method was investigated extensively in the context of reconstruction of Cosmic Microwave Background (CMB) temperature anisotropy signal using observations obtained by WMAP and Planck satellite missions. In this article, we, for the first time, apply the ILC method to reconstruct the large scale CMB E mode polarization signal, which serves as the unique probe of ionization history of the Universe, using simulated observations of 15 frequency CMB polarization maps of future generation COrE satellite mission. We find that usual ILC cleaned E mode map is highly erroneous due to presence of a chance-correlation between CMB and astrophysical foreground components in the empirical covariance matrix which is used to estimate the weight factors. The cleaned angular power spectrum for E mode is strongly biased and erroneous due to these chance correlation factors. In order to address the issues of bias and errors we extend and improve the usual ILC method for CMB E mode reconstruction by incorporating prior information of theoretical E mode angular power spectrum while estimating the weights for linear combination of input maps. Using the E mode covariance matrix effectively suppresses the CMB-foreground chance correlation power leading to an accurate reconstruction of cleaned CMB E mode map and its angular power spectrum. We provide a comparative study of the performance of the usual ILC and the new method over large angular scales of the sky and show that the later produces significantly statistically improved results than the former. The new E mode CMB angular power spectrum contains neither any significant negative bias at the low multipoles nor any positive foreground bias at relatively higher mutlipoles. The error estimates of the cleaned spectrum agree very well with the cosmic variance induced error.
△ Less
Submitted 30 March, 2020;
originally announced March 2020.
-
Separating CMB Stokes Q and U polarization signals from Non-Gaussian Emissions
Authors:
Ujjal Purkayastha,
Rajib Saha
Abstract:
In this paper we estimate diffuse foreground minimized Cosmic Microwave Background (CMB) Stokes Q and U polarization maps based upon the fundamental concept of Gaussian nature of CMB and strong non-Gaussian nature of astrophysical polarized foregrounds using WMAP nine year published polarization maps. We excise regions of the sky that define position of the known point sources, regions that are st…
▽ More
In this paper we estimate diffuse foreground minimized Cosmic Microwave Background (CMB) Stokes Q and U polarization maps based upon the fundamental concept of Gaussian nature of CMB and strong non-Gaussian nature of astrophysical polarized foregrounds using WMAP nine year published polarization maps. We excise regions of the sky that define position of the known point sources, regions that are strongly contaminated by either the detector noise or by the diffuse foregrounds or both, and then perform foreground minimizations over the surviving sky regions that constitute approximately $50\%$ of the full sky area. We critically evaluate performance of foreground minimizations in several ways and show that our foreground minimization method removes significant foregrounds from input maps. The cleaned Stokes \{Q, U\} polarization maps have less EE and BB power from relevant sky region compared to WMAP foreground-reduced Stokes \{Q, U\} polarization maps at different multipole ranges. We validate our methodology by performing detailed Monte Carlo simulations. The main driving machinery of our method is an internal-linear-combination (ILC) approach, however, unlike simple variance minimization performed in usual ILC method, the core of the method being dependent on the theoretically well motivated concept of Gaussianity of CMB polarization a direct connection is established between observations and models of inflation. Additionally, the method, like the usual ILC method, is independent on modeling uncertainties of polarized foregrounds. It will be useful to apply our method in future generation low-noise CMB polarization experiments.
△ Less
Submitted 6 July, 2017;
originally announced July 2017.