We examine the statistical properties of polarization maps from Planck 2018 within the patch of sky observed by the BICEP2/Keck experiment using the one point distribution function (1PDF), skewness, and kurtosis statistics. Our analysis is performed for the Q and U Stokes parameters and for the corresponding E- and B-modes of the CMB signals. We extend our analysis by studying the correlations between CMB polarization maps and residual maps (the difference between the full signal and the CMB map) for the frequency range of 100--217 GHz with both the Q/U and E/B approaches. Although all the CMB maps reveal almost Gaussian statistical properties for Q/U and E/B domains, we have detected very significant anomalies for cross-correlations with residuals at 100 GHz at the level of $3.7\sigma$ for the Commander map and $5.2\sigma$ for NILC, for both the Q and U parameters. Using the NILC--Commander difference, which does not contain a cosmological signal, we find a sub-dominant non-Gaussian component in Q skewness and kurtosis at the level of $4.3\sigma$ and $10\sigma$, respectively. For the B-mode we have found a very high level of cross-correlation (0.63--0.69) between the NILC/Commander maps and the 143 GHz total signal, which cannot be associated with the cosmological component. These strong deviations suggest that remnants of foregrounds, systematic effects, and component separation exist in the 2018 Planck CMB polarization maps in the BICEP2 sky area, which is far away from the Galactic plane. Our analysis also demonstrates the preferability of the Q/U domain over E/B for determination of the statistical properties of the derived CMB signals, due to non-locality of the transition Q/U $\rightarrow$ E/B.

We study the distribution functions of the CMB polarization angle $\psi$, focusing on the Planck 2018 CMB maps. We extend the model of Preece & Battye (2014) of Gaussian correlated $Q$ and $U$ Stokes parameters to allow nonzero means. When the variances of $Q$ and $U$ are equal and their covariance and means are zero, the polarization angle is uniformly distributed. Otherwise the uniform distribution will be modulated by harmonics with $2\psi$ and $4\psi$ phases. These modulations are visible in the Planck 2018 CMB maps. Furthermore, the mean value of $U$ is peculiar compared to the power spectrum.

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.

We study the problem of EB-leakage that is associated with incomplete polarized CMB sky. In the blind case that assumes no additional information about the statistical properties and amplitudes of the signal from the missing sky region, we prove that the recycling method (Liu et al.~2018) gives the unique best estimate of the EB-leakage. Compared to the previous method, this method reduces the uncertainties in the BB power spectrum due to EB-leakage by more than one order of magnitude in the most interesting domain of multipoles, where $\ell$ is between $80$ and $200$. This work also provides a useful guideline for observational design of future CMB experiments.

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.

Better understanding of Galactic foregrounds is one of the main obstacles to detection of primordial gravitational waves through measurement of the B mode in the polarized microwave sky. We generalize the method proposed in \cite{2011PhRvD..83h3003B} and decompose the polarization signals into the E and B families directly in the domain of the Stokes $Q$, $U$ parameters as $(Q,U)\equiv(Q_E, U_E)+(Q_B,U_B)$. This also enables an investigation of the morphology and the frequency dependence of these two families, which has been done in the WMAP K, Ka (tracing synchrotron emission) and Planck 2015 HFI maps (tracing thermal dust). The results reveal significant differences in spectra between the E and B families. The spectral index of the E family fluctuates less across the sky than that of the B family, and the same tendency occurs for the polarization angles of the dust and synchrotron channels. The new insight from WMAP and Planck data on the North Polar Spur and BICEP2 zones through our method clearly indicates that these zones are characterized by very low polarization intensity of the B family compared to the E family. We have detected global structure of the B family polarization angles at high Galactic latitudes which cannot be attributed to the cosmic microwave background or instrumental noise. However, we cannot exclude instrumental systematics as a partial contributor to these anomalies.

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.