Galactic foregrounds are the main obstacle to observations of the cosmic microwave background (CMB) $B$-mode polarization. In addition to obscuring the inflationary $B$-mode signal by several orders of magnitude, Galactic foregrounds have non-trivial spectral signatures that are partially unknown and distorted by averaging effects along the line-of-sight, within the pixel/beam window, and by various analysis choices (e.g., spherical harmonic transforms and filters). Statistical moment expansion methods provide a powerful tool for modeling the effective Galactic foreground emission resulting from these averaging effects in CMB observations, while blind component separation treatments can handle unknown foregrounds. In this work, we combine these two approaches to develop a new semi-blind component separation method at the intersection of parametric and blind methods, called constrained moment ILC (cMILC). This method adds several constraints to the standard ILC method to de-project the main statistical moments of the Galactic foreground emission. Applications to maps are performed in needlet space and when compared to the NILC method, this helps significantly reducing residual foreground contamination (bias, variance, and skewness) in the reconstructed CMB $B$-mode map, power spectrum, and tensor-to-scalar ratio. We consider sky-simulations for experimental settings similar to those of LiteBIRD and PICO, illustrating which trade-offs between residual foreground biases and degradation of the constraint on $r$ can be expected within the new cMILC framework. We also outline several directions that require more work in preparation for the coming analysis challenges.

Cosmological Birefringence (CB) is a phenomenon, caused by parity violating modifications to electrodynamics, whereby the linear polarisation angle of light changes as photons traverse a vacuum. It is possible to use a number of different analysis techniques to constrain this effect using Cosmic Microwave Background (CMB) polarisation observations. We investigate two different methods of constraining direction dependent birefringence for present and future CMB experiments including BICEP/Keck, Simons Observatory (SO), and LiteBIRD . Specifically we compare the constraints placed on anisotropic CB from a quadratic estimator technique to those derived from estimates of the $B$-mode power-spectrum for the three different experiments. The constraints derived from estimates of the $B$-mode power spectrum are found to be comparable to those derived from quadratic estimator for BICEP/Keck and SO, but not LiteBIRD due to its larger sky coverage. These forecasted upper bounds for CB are converted to constraints on primordial magnetic fields and the coupling between photons and pseudo Nambu-Goldstone bosons. Finally we show that even with the best constraints on CB, for the respective experiments, the potentially induced $B$-mode power can act as a significant contaminant in the prospective measurement of primordial $B$-modes.

The cosmological recombination radiation (CRR) is one of the inevitable $\Lambda$CDM spectral distortions of the cosmic microwave background (CMB). While it shows a rich spectral structure across dm-mm wavelengths, it is also one of the smallest signals to target. Here we carry out a detailed forecast for the expected sensitivity levels required to not only detect but also extract cosmological information from the CRR in the presence of foregrounds. We use ${\tt CosmoSpec}$ to compute the CRR including all important radiative transfer effects and modifications to the recombination dynamics. We confirm that detections of the overall CRR signal are possible with spectrometer concepts like ${\it SuperPIXIE}$. However, for real exploitation of the cosmological information, a $\simeq 50$ times more sensitive spectrometer is required. While extremely futuristic, this could provide independent constraints on the primordial helium abundance, $Y_p$, and probe the presence of extra relativistic degrees of freedom during BBN and recombination. Significantly improving the constraints on other cosmological parameters requires even higher sensitivity (another factor of $\simeq 5$) when considering a combination of a CMB spectrometer with existing CMB data. To a large part this is due to astrophysical foregrounds which interestingly do not degrade the constraints on $Y_p$ and $N_{\rm eff}$ as much. A future CMB spectrometer could thus open a novel way of probing non-standard BBN scenarios, dark radiation and sterile neutrinos. In addition, inflation physics could be indirectly probed using the CRR in combination with existing and forthcoming CMB anisotropy data.

The method of weighted addition of multi-frequency maps, more commonly referred to as {\it Internal Linear Combination} (ILC), has been extensively employed in the measurement of cosmic microwave background (CMB) anisotropies and its secondaries along with similar application in 21cm data analysis. Here we argue and demonstrate that ILC methods can also be applied to data from absolutely-calibrated CMB experiments to extract average-sky signals in addition to the conventional CMB anisotropies. The performance of the simple ILC method is, however, limited, but can be significantly improved by adding constraints informed by physics and existing empirical information. In recent work, a moment description has been introduced as a technique of carrying out high precision modeling of foregrounds in the presence of inevitable averaging effects. We combine these two approaches to construct a heavily constrained form of the ILC, dubbed \milc, which can be used to recover tiny monopolar spectral distortion signals in the presence of realistic foregrounds and instrumental noise. This is a first demonstration for measurements of the monopolar and anisotropic spectral distortion signals using ILC and extended moment methods. We also show that CMB anisotropy measurements can be improved, reducing foreground biases and signal uncertainties when using the \milc. While here we focus on CMB spectral distortions, the scope extends to the 21cm monopole signal and $B$-mode analysis. We briefly discuss augmentations that need further study to reach the full potential of the method.

The characterization of the spectral energy distribution (SED) of dust emission has become a critical issue in the quest for primordial B-modes. The dust SED is often approximated by a modified blackbody (MBB) emission law but to what extent is it accurate? This paper addresses this question expanding the dust SED, at the power spectrum level, in moments around the MBB law, related to derivatives with respect to the dust spectral index. We present the mathematical formalism and apply it to simulations and Planck total intensity data, from 143 to 857 GHz, because today no polarized data provide the required sensitivity to perform this analysis. With simulations, we demonstrate the ability of high order moments to account for spatial variations in MBB parameters. Neglecting these moment leads to poor fits and a bias in the recovered dust spectral index. We identify the main moments that are required to fit the Planck data. The comparison with simulations help us to disentangle the respective contributions from dust and the Cosmic Infrared Background to the high order moments, but the simulations give an insufficient description of the actual Planck data. Extending our model to CMB B-mode analyses within a simplified framework, we find that ignoring the dust SED distortions, or trying to model them with a single decorrelation parameter, could lead to biases larger than the targeted sensitivity for the next generation of CMB B-mode experiments.

We develop an analytic model for the power spectra of polarized filamentary structures as a way to study the Galactic polarization foreground to the Cosmic Microwave Background. Our approach is akin to the cosmological halo-model framework, and reproduces the main features of the Planck 353 GHz power spectra. We model the foreground as randomly-oriented, three-dimensional, spheroidal filaments, accounting for their projection onto the sky. The main tunable parameters are the distribution of filament sizes, the filament physical aspect ratio, and the dispersion of the filament axis around the local magnetic field direction. The abundance and properties of filaments as a function of size determine the slopes of the foreground power spectra, as we show via scaling arguments. The filament aspect ratio determines the ratio of $B$-mode power to $E$-mode power, and specifically reproduces the Planck-observed dust ratio of one-half when the short axis is roughly one-fourth the length of the long axis. Filament misalignment to the local magnetic field determines the $TE$ cross-correlation, and to reproduce Planck measurements, we need a (three-dimensional) misalignment angle with a root mean squared dispersion of about 50 degrees. These parameters are not sensitive to the particular filament density profile. By artificially skewing the distribution of the misalignment angle, this model can reproduce the Planck-observed (and parity-violating) $TB$ correlation. The skewing of the misalignment angle necessary to explain $TB$ will cause a yet-unobserved, positive $EB$ dust correlation, a possible target for future experiments.

This paper discusses the science case for a sensitive spectro-polarimetric survey of the microwave sky. Such a survey would provide a tomographic and dynamic census of the three-dimensional distribution of hot gas, velocity flows, early metals, dust, and mass distribution in the entire Hubble volume, exploit CMB temperature and polarisation anisotropies down to fundamental limits, and track energy injection and absorption into the radiation background across cosmic times by measuring spectral distortions of the CMB blackbody emission. In addition to its exceptional capability for cosmology and fundamental physics, such a survey would provide an unprecedented view of microwave emissions at sub-arcminute to few-arcminute angular resolution in hundreds of frequency channels, a data set that would be of immense legacy value for many branches of astrophysics. We propose that this survey be carried-out with a large space mission featuring a broad-band polarised imager and a moderate resolution spectro-imager at the focus of a 3.5m aperture telescope actively cooled to about 8K, complemented with absolutely-calibrated Fourier Transform Spectrometer modules observing at degree-scale angular resolution in the 10-2000 GHz frequency range. We propose two observing modes: a survey mode to map the entire sky as well as a few selected wide fields, and an observatory mode for deeper observations of regions of specific interest.

Voyage 2050 White Paper highlighting the unique science opportunities using spectral distortions of the cosmic microwave background (CMB). CMB spectral distortions probe many processes throughout the history of the Universe. Precision spectroscopy, possible with existing technology, would provide key tests for processes expected within the cosmological standard model and open an enormous discovery space to new physics. This offers unique scientific opportunities for furthering our understanding of inflation, recombination, reionization and structure formation as well as dark matter and particle physics. A dedicated experimental approach could open this new window to the early Universe in the decades to come, allowing us to turn the long-standing upper distortion limits obtained with COBE/FIRAS some 25 years ago into clear detections of the expected standard distortion signals.

Departures of the energy spectrum of the cosmic microwave background (CMB) from a perfect blackbody probe a fundamental property of the universe -- its thermal history. Current upper limits, dating back some 25 years, limit such spectral distortions to 50 parts per million and provide a foundation for the Hot Big Bang model of the early universe. Modern upgrades to the 1980's-era technology behind these limits enable three orders of magnitude or greater improvement in sensitivity. The standard cosmological model provides compelling targets at this sensitivity, spanning cosmic history from the decay of primordial density perturbations to the role of baryonic feedback in structure formation. Fully utilizing this sensitivity requires concurrent improvements in our understanding of competing astrophysical foregrounds. We outline a program using proven technologies capable of detecting the minimal predicted distortions even for worst-case foreground scenarios.

The Simons Observatory (SO) is a ground-based cosmic microwave background (CMB) experiment sited on Cerro Toco in the Atacama Desert in Chile that promises to provide breakthrough discoveries in fundamental physics, cosmology, and astrophysics. Supported by the Simons Foundation, the Heising-Simons Foundation, and with contributions from collaborating institutions, SO will see first light in 2021 and start a five year survey in 2022. SO has 287 collaborators from 12 countries and 53 institutions, including 85 students and 90 postdocs. The SO experiment in its currently funded form ('SO-Nominal') consists of three 0.4 m Small Aperture Telescopes (SATs) and one 6 m Large Aperture Telescope (LAT). Optimized for minimizing systematic errors in polarization measurements at large angular scales, the SATs will perform a deep, degree-scale survey of 10% of the sky to search for the signature of primordial gravitational waves. The LAT will survey 40% of the sky with arc-minute resolution. These observations will measure (or limit) the sum of neutrino masses, search for light relics, measure the early behavior of Dark Energy, and refine our understanding of the intergalactic medium, clusters and the role of feedback in galaxy formation. With up to ten times the sensitivity and five times the angular resolution of the Planck satellite, and roughly an order of magnitude increase in mapping speed over currently operating ("Stage 3") experiments, SO will measure the CMB temperature and polarization fluctuations to exquisite precision in six frequency bands from 27 to 280 GHz. SO will rapidly advance CMB science while informing the design of future observatories such as CMB-S4.

Following the pioneering observations with COBE in the early 1990s, studies of the cosmic microwave background (CMB) have focused on temperature and polarization anisotropies. CMB spectral distortions - tiny departures of the CMB energy spectrum from that of a perfect blackbody - provide a second, independent probe of fundamental physics, with a reach deep into the primordial Universe. The theoretical foundation of spectral distortions has seen major advances in recent years, which highlight the immense potential of this emerging field. Spectral distortions probe a fundamental property of the Universe - its thermal history - thereby providing additional insight into processes within the cosmological standard model (CSM) as well as new physics beyond. Spectral distortions are an important tool for understanding inflation and the nature of dark matter. They shed new light on the physics of recombination and reionization, both prominent stages in the evolution of our Universe, and furnish critical information on baryonic feedback processes, in addition to probing primordial correlation functions at scales inaccessible to other tracers. In principle the range of signals is vast: many orders of magnitude of discovery space could be explored by detailed observations of the CMB energy spectrum. Several CSM signals are predicted and provide clear experimental targets, some of which are already observable with present-day technology. Confirmation of these signals would extend the reach of the CSM by orders of magnitude in physical scale as the Universe evolves from the initial stages to its present form. The absence of these signals would pose a huge theoretical challenge, immediately pointing to new physics.

Our velocity relative to the cosmic microwave background (CMB) generates a dipole from the CMB monopole, which was accurately measured by COBE. The relative velocity also modulates and aberrates the CMB fluctuations, generating a small signature of statistical isotropy violation in the covariance matrix. This signature was first measured by Planck 2013. Galaxy surveys are similarly affected by a Doppler boost. The dipole generated from the number count monopole has been extensively discussed, and measured (at very low accuracy) in the NVSS and TGSS radio continuum surveys. For the first time, we present an analysis of the Doppler imprint on the number count fluctuations, using the bipolar spherical harmonic formalism to quantify these effects. Next-generation wide-area surveys with a high redshift range are needed to detect the small Doppler signature in number count fluctuations. We show that radio continuum surveys with the SKA should enable a detection at $\gtrsim 3 \sigma$ in Phase 2, with marginal detection possible in Phase 1.

Measurements of the thermal Sunyaev-Zel'dovich (tSZ) effect have long been recognized as a powerful cosmological probe. Here we assess the importance of relativistic temperature corrections to the tSZ signal on the power spectrum analysis of the Planck Compton-$y$ map, developing a novel formalism to account for the associated effects. The amplitude of the tSZ power spectrum is found to be sensitive to the effective electron temperature, $\bar{T}_e$, of the cluster sample. Omitting the corresponding modifications leads to an underestimation of the $yy$-power spectrum amplitude. Relativistic corrections thus add to the error budget of tSZ power spectrum observables such as $\sigma_8$. This could help alleviate the tension between various cosmological probes, with the correction scaling as $\Delta \sigma_8/\sigma_8 \simeq 0.019\,[k\bar{T}_e\,/\,5\,{\rm keV}]$ for Planck. At the current level of precision, this implies a systematic shift by $\simeq 1\sigma$, which can also be interpreted as an overestimation of the hydrostatic mass bias by $\Delta b \simeq 0.046\,(1-b)\,[k\bar{T}_e\,/\,5\,{\rm keV}]$, bringing it into better agreement with hydrodynamical simulations. It is thus time to consider relativistic temperature corrections in the processing of current and future tSZ data.

We derive full-sky, real-space operators that convert between polarization Stokes Q/U parameters and the coordinate-independent scalar E/B modes that are widely used in Cosmic Microwave Background (CMB) and cosmic shear analysis. We also derive real space operators that decompose the measured Stokes parameters into those corresponding to E-modes and B-modes respectively, without ever evaluating the scalar fields themselves. For all these real space operators we show that the kernels split naturally into angular and radial parts and we show explicitly how the radial extent of these kernels depends on the targeted band-limit. The kernels can be interpreted either as a complex convolving beam or as a Green's function when they are expressed in terms of the forward or inverse rotation Euler angles. We show that an arbitrary radial function can produce E/B-like maps, provided it vanishes at the origin and the antipodal point. These maps are simply filtered versions of the standard E/B maps. We argue that it is possible compute E/B maps in real space with a compactly-supported kernel, an approach that can guarantee the avoidance of known foreground regions and could be employed in a massively-parallel scheme at high-resolution. We show that the spin raising and lowering operators $\eth^2$/$\bar{\eth}^2$ are special cases of these generalized radial functions, and present their band limited versions. The spatial structure of the real space operators provides great intuition for the E/B structure of polarized, filamentary galactic foregrounds. We predict a non-zero B-mode signature that is expected from polarized filaments in the sky. This paper is the first part in a series of papers that explore real-space computation of polarization modes and their applications.

This book lays out the scientific goals to be addressed by the next-generation ground-based cosmic microwave background experiment, CMB-S4, envisioned to consist of dedicated telescopes at the South Pole, the high Chilean Atacama plateau and possibly a northern hemisphere site, all equipped with new superconducting cameras. CMB-S4 will dramatically advance cosmological studies by crossing critical thresholds in the search for the B-mode polarization signature of primordial gravitational waves, in the determination of the number and masses of the neutrinos, in the search for evidence of new light relics, in constraining the nature of dark energy, and in testing general relativity on large scales.

Isotropy-violation statistics can highlight polarized galactic foregrounds that contaminate primordial $B$-modes in the Cosmic Microwave Background (CMB). We propose a particular isotropy-violation test and apply it to polarized Planck 353 GHz data, constructing an map that indicates $B$-mode foreground dust power over the sky. We build our main isotropy test in harmonic space via the bipolar spherical harmonic basis, and our method helps us to identify the least-contaminated directions. By this measure, there are regions of low foreground in and around the BICEP field, near the South Galactic Pole, and in the Northern Galactic Hemisphere. There is also a possible foreground feature in the BICEP field. We compare our results to those based on the local power spectrum, which is computed on discs using a version of the method of Planck Int.~XXX (2016). The discs method is closely related to our isotropy-violation diagnostic. We pay special care to the treatment of noise, including chance correlations with the foregrounds. Currently we use our isotropy tool to assess the cleanest portions of the sky, but in the future such methods will allow isotropy-based null tests for foreground contamination in maps purported to measure primordial $B$-modes, particularly in cases of limited frequency coverage.

A stochastic gravitational wave background (SGWB) will affect the CMB anisotropies via weak lensing. Unlike weak lensing due to large scale structure which only deflects photon trajectories, a SGWB has an additional effect of rotating the polarization vector along the trajectory. We study the relative importance of these two effects, deflection \& rotation, specifically in the context of E-mode to B-mode power transfer caused by weak lensing due to SGWB. Using weak lensing distortion of the CMB as a probe, we derive constraints on the spectral energy density ($\Omega_{GW}$) of the SGWB, sourced at different redshifts, without assuming any particular model for its origin. We present these bounds on $\Omega_{GW}$ for different power-law models characterizing the SGWB, indicating the threshold above which observable imprints of SGWB must be present in CMB.

Statistical isotropy (SI) of Cosmic Microwave Background (CMB) fluctuations is a key observational test to validate the cosmological principle underlying the standard model of cosmology. While a detection of SI violation would have immense cosmological ramification, it is important to recognise their possible origin in systematic effects of observations. WMAP seven year (WMAP-7) release claimed significant deviation from SI in the bipolar spherical harmonic (BipoSH) coefficients $A_{ll}^{20}$ and $A_{l-2l}^{20}$. Here we present the first explicit reproduction of the measurements reported in WMAP-7, confirming that beam systematics alone can completely account for the measured SI violation. The possibility of such a systematic origin was alluded to in WMAP-7 paper itself and other authors but not as explicitly so as to account for it accurately. We simulate CMB maps using the actual WMAP non-circular beams and scanning strategy. Our estimated BipoSH spectra from these maps match the WMAP-7 results very well. It is also evident that only a very careful and adequately detailed modelling, as carried out here, can conclusively establish that the entire signal arises from non-circular beam effect. This is important since cosmic SI violation signals are expected to be subtle and dismissing a large SI violation signal as observational artefact based on simplistic plausibility arguments run the serious risk of "throwing the baby out with the bathwater".

We test the statistical isotropy and Gaussianity of the cosmic microwave background (CMB) anisotropies using observations made by the Planck satellite. Our results are based mainly on the full Planck mission for temperature, but also include some polarization measurements. In particular, we consider the CMB anisotropy maps derived from the multi-frequency Planck data by several component-separation methods. For the temperature anisotropies, we find excellent agreement between results based on these sky maps over both a very large fraction of the sky and a broad range of angular scales, establishing that potential foreground residuals do not affect our studies. Tests of skewness, kurtosis, multi-normality, N-point functions, and Minkowski functionals indicate consistency with Gaussianity, while a power deficit at large angular scales is manifested in several ways, for example low map variance. The results of a peak statistics analysis are consistent with the expectations of a Gaussian random field. The "Cold Spot" is detected with several methods, including map kurtosis, peak statistics, and mean temperature profile. We thoroughly probe the large-scale dipolar power asymmetry, detecting it with several independent tests, and address the subject of a posteriori correction. Tests of directionality suggest the presence of angular clustering from large to small scales, but at a significance that is dependent on the details of the approach. We perform the first examination of polarization data, finding the morphology of stacked peaks to be consistent with the expectations of statistically isotropic simulations. Where they overlap, these results are consistent with the Planck 2013 analysis based on the nominal mission data and provide our most thorough view of the statistics of the CMB fluctuations to date.

Mild, unavoidable deviations from circular-symmetry of instrumental beams along with scan strategy can give rise to measurable Statistical Isotropy (SI) violation in Cosmic Microwave Background (CMB) experiments. If not accounted properly, this spurious signal can complicate the extraction of other SI violation signals (if any) in the data. However, estimation of this effect through exact numerical simulation is computationally intensive and time consuming. A generalized analytical formalism not only provides a quick way of estimating this signal, but also gives a detailed understanding connecting the leading beam anisotropy components to a measurable BipoSH characterisation of SI violation. In this paper, we provide an approximate generic analytical method for estimating the SI violation generated due to a non-circular (NC) beam and arbitrary scan strategy, in terms of the Bipolar Spherical Harmonic (BipoSH) spectra. Our analytical method can predict almost all the features introduced by a NC beam in a complex scan and thus reduces the need for extensive numerical simulation worth tens of thousands of CPU hours into minutes long calculations. As an illustrative example, we use WMAP beams and scanning strategy to demonstrate the easability, usability and efficiency of our method. We test all our analytical results against that from exact numerical simulations.

Any isotropy violating phenomena on cosmic microwave background (CMB) induces off-diagonal correlations in the two-point function. These correlations themselves can be used to estimate the underlying anisotropic signals. Masking due to residual foregrounds, or availability of partial sky due to survey limitation, are unavoidable circumstances in CMB studies. But, masking induces additional correlations, and thus complicates the recovery of such signals. In this work, we discuss a procedure based on bipolar spherical harmonic (BipoSH) formalism to comprehensively addresses any spurious correlations induced by masking and successfully recover hidden signals of anisotropy in observed CMB maps. This method is generic, and can be applied to recover a variety of isotropy violating phenomena. Here, we illustrate the procedure by recovering the subtle Doppler boost signal from simulated boosted CMB skies, which has become possible with the unprecedented full-sky sensitivity of PLANCK probe.

The European Space Agency's Planck satellite, dedicated to studying the early Universe and its subsequent evolution, was launched 14~May 2009 and scanned the microwave and submillimetre sky continuously between 12~August 2009 and 23~October 2013. In February~2015, ESA and the Planck Collaboration released the second set of cosmology products based on data from the entire Planck mission, including both temperature and polarization, along with a set of scientific and technical papers and a web-based explanatory supplement. This paper gives an overview of the main characteristics of the data and the data products in the release, as well as the associated cosmological and astrophysical science results and papers. The science products include maps of the cosmic microwave background (CMB), the thermal Sunyaev-Zeldovich effect, and diffuse foregrounds in temperature and polarization, catalogues of compact Galactic and extragalactic sources (including separate catalogues of Sunyaev-Zeldovich clusters and Galactic cold clumps), and extensive simulations of signals and noise used in assessing the performance of the analysis methods and assessment of uncertainties. The likelihood code used to assess cosmological models against the Planck data are described, as well as a CMB lensing likelihood. Scientific results include cosmological parameters deriving from CMB power spectra, gravitational lensing, and cluster counts, as well as constraints on inflation, non-Gaussianity, primordial magnetic fields, dark energy, and modified gravity.

The Bipolar Spherical Harmonics (BipoSH) form a natural basis to study the CMB two point correlation function in a non-statistically isotropic (non-SI) universe. The coefficients of expansion in this basis are a generalisation of the well known CMB angular power spectrum and contain complete information of the statistical properties of a non-SI but Gaussian random CMB sky. We use these coefficients to describe the weak lensing of CMB photons in a non-SI universe. Finally we show that the results reduce to the standard weak lensing results in the isotropic limit.

Statistical isotropy (SI) is one of the fundamental assumptions made in cosmological model building. This assumption is now being rigorously tested using the almost full sky measurements of the CMB anisotropies. A major hurdle in any such analysis is to handle the large biases induced due to the process of masking. We have developed a new method of analysis, using the bipolar spherical harmonic basis functions, in which we semi-analytically evaluate the modifications to SI violation induced by the mask. The method developed here is generic and can be potentially used to search for any arbitrary form of SI violation. We specifically demonstrate the working of this method by recovering the Doppler boost signal from a set of simulated, masked CMB skies.

The two point correlation function of the CMB temperature anisotropies is generally assumed to be statistically isotropic (SI). Deviations from this assumption could be traced to physical or observational artefacts and systematic effects. Measurement of non-vanishing power in the BipoSH spectra is a standard statistical technique to search for isotropy violations. Although this is a neat tool allowing a blind search for SI violations in the CMB sky, it is not easy to discern the cause of isotropy violation using this measure. In this article, we propose a novel technique of constructing orthogonal BipoSH estimators, which can be used to discern between models of isotropy violation.

The two fundamental assumptions of the standard cosmological model - that the initial fluctuations are statistically isotropic and Gaussian - are rigorously tested using maps of the cosmic microwave background (CMB) anisotropy from the Planck satellite. Deviations from isotropy have been found and demonstrated to be robust against component separation algorithm, mask choice and frequency dependence. Many of these anomalies were previously observed in the WMAP data, and are now confirmed at similar levels of significance (about 3 sigma). However, we find little evidence for non-Gaussianity, with the exception of a few statistical signatures that seem to be associated with specific anomalies. In particular, we find that the quadrupole-octopole alignment is also connected to a low observed variance of the CMB signal. A power asymmetry is now found to persist to scales corresponding to about l=600, and can be described in the low-l regime by a phenomenological dipole modulation model. However, any primordial power asymmetry is strongly scale-dependent and does not extend to arbitrarily small angular scales. Finally, it is plausible that some of these features may be reflected in the angular power spectrum of the data, which shows a deficit of power on similar scales. Indeed, when the power spectra of two hemispheres defined by a preferred direction are considered separately, one shows evidence for a deficit in power, while its opposite contains oscillations between odd and even modes that may be related to the parity violation and phase correlations also detected in the data. Although these analyses represent a step forward in building an understanding of the anomalies, a satisfactory explanation based on physically motivated models is still lacking.

A stochastic gravitational wave background (SGWB) would gravitationally lens the cosmic microwave background (CMB) photons. We find that the lensing due to gravitational waves(GW) is more efficient as compared to lensing due to scalar density perturbations. Though the effect of lensing due to GW is found to effect all the four CMB power spectra, its effect is most prominently seen in the CMB polarization power spectra. This suggests that the measurements of the CMB angular power spectra could be used to constraining the energy density ($\Omega_{GW}$) of the SGBW. In our analysis we find that the most stringent constraints on $\Omega_{GW}$ are due to measurements of the angular power spectra of CMB temperature anisotropies. We show that in the future it will be possible to place more stringent bounds on $\Omega_{GW}$ using improved upper limits or detections of the angular power spectra of the B-modes of CMB polarization at large multipoles.

We provide a detailed treatment and comparison of the weak lensing effects due to large-scale structure (LSS), or scalar density perturbations and those due to gravitational waves(GW) or tensor perturbations, on the temperature and polarization power spectra of the Cosmic Microwave Background (CMB). We carry out the analysis both in real space by using the correlation function method, as well as in the spherical harmonic space. We find an intriguing similarity between the lensing kernels associated with LSS lensing and GW lensing. It is found that the lensing kernels only differ in relative negative signs and their form is very reminiscent of even and odd parity bipolar spherical harmonic coefficients. Through a numerical study of these lensing kernels, we establish that lensing due to GW is more efficient at distorting the CMB spectra as compared to LSS lensing, particularly for the polarization power spectra. Finally we argue that the CMB B-mode power spectra measurements can be used to place interesting constraints on GW energy densities.

Mild, unavoidable deviations from circular-symmetry of instrumental beams in current Cosmic Microwave Background (CMB) experiments now pose a significant challenge to deriving high precision inferences from the high sensitivity and resolution of CMB measurements. We present analytic results, verified by numerical simulations, that CMB maps of cosmological signal that respect underlying statistical isotropy (SI) symmetry, measured with an instrument that has mildly non-circular (NC) beam would, nevertheless, exhibit SI violation. Further, we show that appropriate observable measures constructed within the Bipolar spherical harmonic (BipoSH) representation of SI violation capture subtle NC-beam effects coupled with the scan strategy of the instrument. Accompanying their latest 7-year data release, the WMAP team published very high significance measurements of non-zero BipoSH spectra, A^{20}_{l l} and A^{20}_{l-2 l}, in the "W" and "V" band of the experiment. We present a strong case that the BipoSH measurement are primarily explained by the quadrupolar (m=2) component of the NC-beams, b_{l2}, of the respective channels. The fact that subtle levels of non-circularity, e.g., WMAP beams |b_{l2}|/b_{l0} < 0.01, lead to measurable BipoSH spectra points to the immense promise and potential of the BipoSH representation. The key result of this work is that using BipoSH measurements it is possible to estimate an equivalent single hit,'parallel transported' effective NC-beam that matches both the angular power spectrum and the non-zero BipoSH measurements of the observed maps. Hence, BipoSH analysis provides a very simple and effective characterization of CMB maps made with NC-beam. (Abridged)

The Integrated Sachs-Wolfe (ISW) effect produces a secondary temperature anisotropy of CMB. The main contribution comes from z<2, where dark energy leads to a decay of potentials. As the same photons are gravitationally lensed by these decaying potentials, there exists a high degree of correlation between the ISW effect and CMB lensing, leading to a non-zero three-point correlation (bispectrum) of the observed temperature anisotropy. This ISW-lensing bispectrum, whose shape resembles that of the so-called "local-form" primordial bispectrum parametrized by fNL, is known to be the largest contamination of fNL. In order to avoid a spurious detection of primordial non-Gaussianity, we need to remove the ISW-lensing bias. In this work, we investigate three debiasing methods: (I) subtraction of an expected, ensemble average of the ISW-lensing bispectrum; (II) subtraction of a measured ISW-lensing bispectrum; and (III) direct subtraction of an estimated ISW signal from an observed temperature map. One may use an estimation of the ISW map from external non-CMB data or that from the CMB data themselves. As the methods II and III are based on fewer assumptions about the nature of dark energy, they are preferred over the method I. While the methods I and II yield unbiased estimates of fNL with comparable error bars, the method III yields a biased result when the underlying primordial fNL is non-zero and the ISW map is estimated from a lensing potential reconstructed from the observed temperature map. One of the sources of the bias is a lensing reconstruction noise bias which is independent of fNL and can be calculated precisely, but other fNL-dependent terms are difficult to compute reliably. We thus conclude that the method II is the best, model-independent way to remove the ISW-lensing bias of fNL, enabling us to test the physics of inflation with smaller systematic errors.

Statistical isotropy (SI) has been one of the simplifying assumptions in cosmological model building. Experiments like WMAP and PLANCK are attempting to test this assumption by searching for specific signals in the Cosmic Microwave Background (CMB) two point correlation function. Modifications to this correlation function due to gravitational lensing by the large scale structure (LSS) surrounding us have been ignored in this context. Gravitational lensing will induce signals which mimic isotropy violation even in an isotropic universe. The signal detected in the Bipolar Spherical Harmonic (BipoSH) coefficients $A^{20}_{ll}$ by the WMAP team may be explained by accounting for the lensing modifications to these coefficients. Further the difference in the amplitude of the signal detected in the V-band and W-band maps can be explained by accounting for the differences in the designed angular sensitivity of the instrumental beams. The arguments presented in this article have crucial implications for SI violation studies. Constraining SI violation will only be possible by complementing CMB data sets with all sky measurements of the large scale dark matter distribution. Till that time, the signal detected in the BipoSH coefficients from WMAP-7 could also be yet another suggested evidence of strong deviations from the standard $\Lambda$CDM cosmology based on homogeneous and isotropic FRW models.

Gaussianity of temperature fluctuations in the Cosmic Microwave Background(CMB) implies that the statistical properties of the temperature field can be completely characterized by its two point correlation function. The two point correlation function can be expanded in full generality in the bipolar spherical harmonic(BipoSH) basis. Looking for significant deviations from zero for Bipolar Spherical Harmonic(BipoSH) Coefficients derived from observed CMB maps forms the basis of the strategy used to detect isotropy violation. In order to quantify "significant deviation" we need to understand the distributions of these coefficients. We analytically evaluate the moments and the distribution of the coefficients of expansion($A^{LM}_{l_1 l_2}$), using characteristic function approach. We show that for BipoSH coefficients with M=0 an analytical form for the moments up to any arbitrary order can be derived. For the remaining BipoSH coefficients with $M\neq0$, the moments derived using the characteristic function approach need to be supplemented with a correction term. The correction term is found to be important particularly at low multipoles. We provide a general prescription for calculating these corrections, however we restrict the explicit calculations only up to kurtosis. We confirm our results with measurements of BipoSH coefficients on numerically simulated statistically isotropic CMB maps.