In observation of the cosmic microwave background (CMB) polarization, "$EB$~leakage" refers to the artificial $B$-mode signal coming from the leakage of $E$-mode signal when part of the sky is unavailable or excluded. Correction of such leakage is one of the preconditions for detecting primordial gravitational waves via the CMB $B$-mode signal. In this work, we design two independent methods for correcting the $EB$~leakage directly in the pixel domain using standard definitions of the $E$- and $B$-modes. The two methods give consistent results, and both are fast and easy to implement. Tests on a CMB simulation containing zero initial $B$-mode show an efficient suppression of the $EB$ leakage. When combined with the MASTER method to reconstruct the full-sky $B$-mode spectrum in simulations with a relatively simple mask, the error from EB-leakage is suppressed further by more than one order of magnitude at the recombination bump, and up to three orders of magnitude at higher multipoles, compared to a "pure" MASTER scheme under the same conditions. Meanwhile, although the final power spectrum estimation benefits from apodization, the pixel domain correction itself is done without apodization, and thus the methods offer more freedom in choosing an apodization based on specific requirements.

Foreground removal techniques for CMB analyses make specific assumptions about the properties of foregrounds in temperature and in polarization. By investigating the statistics of foreground components more understanding about the degree to which these assumptions are valid can be obtained. In this work we investigate $E$- and $B$-mode maps of the two strongest polarized foregrounds, synchrotron and thermal dust emission, with regards to their similarity with Gaussian processes, their spectral variations and cross-correlations. We perform tests in patches of $\sim3.7^\circ$ size collectively covering the full sky and find most of them to conform with their Gaussian expectation according to the statistics in use. Correlations exhibit distinct differences in $E$- and $B$-mode signals which point towards necessities in foreground removal methods. We discuss potential consequences and possible further directions.

Precise measurement of the foreground synchrotron emission, which contaminates the faint polarized cosmic microwave background radiation (CMB), is a major challenge for the next-generation of CMB experiments. To address this, dedicated foreground measurement experiments are being undertaken at radio frequencies between 2 and 40 GHz. Foreground polarized synchrotron emission measurements are particularly challenging, primarily due to the complicated frequency dependence in the presence of Faraday rotation, and are best recovered through broad fractional-bandwidth polarization measurements at frequencies $\lesssim 5$ GHz. A unique opportunity for measuring the foreground polarized synchrotron emission will be provided by the 15-m SKA-MPG telescope operating in the frequency range 1.7 to 3.5~GHz (S-Band). Here, we present the scope of a Southern sky survey in S-Band at 1 degree angular resolution and explore its added advantage for application of powerful techniques, such as, Stokes $Q$, $U$ fitting and RM-synthesis. A full Southern-sky polarization survey with this telescope, when combined with other on-going efforts at slightly higher frequencies, will provide an excellent frequency coverage for modeling and extrapolating the foreground polarized synchrotron emission to CMB frequencies ($\gtrsim80$~GHz) with rms brightness temperature better than 10~nK per 1 degree$^2$. We find that this survey will be crucial for understanding the effects of Faraday depolarization, especially in low Galactic latitude regions. This will allow better foreground cleaning and thus will contribute significantly in further improving component separation analyses and increase usable sky area for cosmological analysis of the \textit{Planck} data, and the \textit{LiteBIRD} mission in the future.

In CMB science, the simplest idea to remove a contaminated sky region is to multiply the sky map with a mask that is 0 for the contaminated region and 1 elsewhere, which is also called a top-hat masking. Although it is easy to use, such top-hat masking is known to suffer from various leakage problems. Therefore, we want to extend the top-hat masking to a series of semi-analytic functions called the P-filters. Most importantly, the P-filters can seamlessly realize the core idea of masking in CMB science, and, meanwhile, guarantee continuity up to the first derivative everywhere. The P-filters can significantly reduce many leakage problems without additional cost, including the leakages due to low-, high-, and band-pass filtering, and the E-to-E, B-to-B, B-to-E, and E-to-B leakages. The workings of the P-filter are illustrated by using the WMAP and Planck polarization sky maps. By comparison to the corresponding WMAP/Planck masks, we show that the P-filter performs much better than top-hat masking, and meanwhile, has the potential to supersede the principal idea of masking in CMB science. Compared to mask apodization, the P-filter is ``outward'', that tends to make proper use of the region that was marked as 0; whereas apodization is ``inward'', that always kills more signal in the region marked as 1.

We study the degeneracy of theoretical gravitational waveforms for binary black hole mergers using an aligned-spin effective-one-body model. After appropriate truncation, bandpassing, and matching, we identify regions in the mass--spin parameter space containing waveforms similar to the template proposed for GW150914, with masses $m_1 = 36^{+5}_{-4} M_\odot$ and $m_2 = 29^{+4}_{-4} M_\odot$, using the cross-correlation coefficient as a measure of the similarity between waveforms. Remarkably high cross-correlations are found across broad regions of parameter space. The associated uncertanties exceed these from LIGO's Bayesian analysis considerably. We have shown that waveforms with greatly increased masses, such as $m_1 = 70 M_\odot$ and $m_2 = 35 M_\odot$, and strong anti-aligned spins ($\chi_1=0.95$ and $\chi_2=-0.95$) yield almost the same signal-to-noise ratio in the strain data for GW150914.

We propose a blind, template-free method for the extraction of a common signal between the Hanford and Livingston detectors and apply it especially to the GW150914 event. We construct a log-likelihood method that maximizes the cross-correlation between each detector and the common signal and minimizes the cross-correlation between the residuals. The reliability of this method is tested using simulations with an injected common signal. Finally, our method is used to assess the quality of theoretical gravitational wave templates for GW150914.

To date, the LIGO collaboration has detected three gravitational wave (GW) events appearing in both its Hanford and Livingston detectors. In this article we reexamine the LIGO data with regard to correlations between the two detectors. With special focus on GW150914, we report correlations in the detector noise which, at the time of the event, happen to be maximized for the same time lag as that found for the event itself. Specifically, we analyze correlations in the calibration lines in the vicinity of 35\,Hz as well as the residual noise in the data after subtraction of the best-fit theoretical templates. The residual noise for the other two events, GW151226 and GW170104, exhibits similar behavior. A clear distinction between signal and noise therefore remains to be established in order to determine the contribution of gravitational waves to the detected signals.

Understanding the properties of dust emission in the microwave domain is an important premise for the next generation of cosmic microwave background (CMB) experiments, devoted to the measurement of the primordial $B$-modes of polarization. In this paper, we compare three solutions to thermal dust emission by the Planck Collaboration \cite{PlanckDust03,planck_16,planck_com} to point out significant differences between their respective parameters (the spectral index $\beta$, the optical depth $\tau$ and the dust temperature $T_d$). These differences originate from e.g. the priors on the parameters or the contribution of the Cosmic infrared background (CIB).

Cosmology has made enormous progress through studies of the cosmic microwave background, however the subtle signals being now sought such as B-mode polarisation due to primordial gravitational waves are increasingly hard to disentangle from residual Galactic foregrounds in the derived CMB maps. We revisit our finding that on large angular scales there are traces of the nearby old supernova remnant Loop I in the WMAP 9-year map of the CMB and confirm this with the new SMICA map from the Planck satellite.

Observational cosmology is entering an era in which high precision will be required in both measurement and data analysis. Accuracy, however, can only be achieved with a thorough understanding of potential sources of contamination from foreground effects. Our primary focus will be on non- Gaussian effects in foregrounds. This issue will be crucial for coming experiments to determine B-mode polarization. We propose a novel method for investigating a data set in terms of skewness and kurtosis in locally defined regions that collectively cover the entire sky. The method is demonstrated on two sky maps: (i) the SMICA map of Cosmic Microwave Background fluctuations provided by the Planck Collaboration and (ii) a version of the Haslam map at 408 MHz that describes synchrotron radiation. We find that skewness and kurtosis can be evaluated in combination to reveal local physical information. In the present case, we demonstrate that the local properties of both maps are predominantly Gaussian. This result was expected for the SMICA map; that it also applies for the Haslam map is surprising. The approach described here has a generality and flexibility that should make it useful in a variety of astrophysical and cosmological contexts.

We calculate weighted mosaic correlations between the recently published Planck 2015 foreground maps - both anomalous microwave emission (AME) maps, free-free emission, synchrotron radiation and thermal dust emission. The weighting coefficients are constructed taking account of the signal-to-error ratio given by the data product. Positive correlation is found for AME compared with thermal dust emission as well as synchrotron radiation. We find AME and free-free emission tending to be anti-correlated, however, when investigating different scales, their relationship appears to be more complex. We argue that dust particles responsible for AME are pushed out of hot zones in the interstellar medium (ISM).