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.

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.

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.