BICEP2 / Keck Array x: Constraints on Primordial Gravitational Waves using Planck, WMAP, and New BICEP2/Keck Observations through the 2015 Season

TL;DR

Results from an analysis of all data taken by the bicep2/Keck CMB polarization experiments up to and including the 2015 observing season are presented, showing the strongest constraints to date on primordial gravitational waves.

Abstract

We present results from an analysis of all data taken by the BICEP2/Keck CMB polarization experiments up to and including the 2015 observing season. This includes the first Keck Array observations at 220 GHz and additional observations at 95 & 150 GHz. The $Q/U$ maps reach depths of 5.2, 2.9 and 26 $μ$K$_{cmb}$ arcmin at 95, 150 and 220 GHz respectively over an effective area of $\approx 400$ square degrees. The 220 GHz maps achieve a signal-to-noise on polarized dust emission approximately equal to that of Planck at 353 GHz. We take auto- and cross-spectra between these maps and publicly available WMAP and Planck maps at frequencies from 23 to 353 GHz. We evaluate the joint likelihood of the spectra versus a multicomponent model of lensed-$Λ$CDM+$r$+dust+synchrotron+noise. The foreground model has seven parameters, and we impose priors on some of these using external information from Planck and WMAP derived from larger regions of sky. The model is shown to be an adequate description of the data at the current noise levels. The likelihood analysis yields the constraint $r_{0.05}<0.07$ at 95% confidence, which tightens to $r_{0.05}<0.06$ in conjunction with Planck temperature measurements and other data. The lensing signal is detected at $8.8 σ$ significance. Running maximum likelihood search on simulations we obtain unbiased results and find that $σ(r)=0.020$. These are the strongest constraints to date on primordial gravitational waves.

Figures (23)

• Figure 1: Maps of degree angular scale $E$-modes ($50<\ell<120$) at three frequencies made using Keck Array data from the 2015 season only. The similarity of the pattern indicates that $\Lambda$CDM $E$-modes dominate at all three frequencies (and that the signal-to-noise is high). The color scale is in $\mu{\mathrm K}$, and the range is allowed to vary slightly to (partially) compensate for the decrease in beam size with increasing frequency.
• Figure 2: $EE$ and $BB$ auto- and cross-spectra calculated using BICEP2/ Keck 95, 150 & 220 GHz maps and the Planck 353 GHz map. The BICEP2/ Keck maps use all data taken up to and including the 2015 observing season---we refer to these as BK15. The black lines show the model expectation values for lensed-$\Lambda$CDM, while the red lines show the expectation values of the baseline lensed-$\Lambda$CDM+dust model from our previous BK14 analysis ($r=0$, $A_\mathrm{d,353}=4.3$$\mu{\mathrm K^2}$, $\beta_\mathrm{d}=1.6$, $\alpha_\mathrm{d}=-0.4$), and the error bars are scaled to that model. Note that the model shown was fit to $BB$ only and did not use the 220 GHz points (which are entirely new). The agreement with the spectra involving 220 GHz and all the $EE$ spectra (under the assumption that $EE/BB=2$ for dust) is therefore a validation of the model. (The dashed red lines show the expectation values of the lensed-$\Lambda$CDM+dust model when adding strong spectral decorrelation of the dust pattern---see Appendix \ref{['app:decorr']} for further information.)
• Figure 3: Upper: The noise spectra of the BK15 maps for 95 GHz (red), 150 GHz (green) and 220 GHz (blue) after correction for the filtering of signal which occurs due to the beam roll-off and timestream filtering. (Note that no $\ell^2$ scaling is applied.) Lower: The effective sky fraction as calculated from the ratio of the mean noise realization bandpowers to their fluctuation $f_\mathrm{sky}(\ell)=\frac{1}{2\ell \Delta \ell} \left( \frac{\sqrt{2}\bar{N_b}}{\sigma(N_b)} \right)^2$, i.e. the observed number of $B$-mode degrees of freedom divided by the nominal full-sky number. The turn-down at low $\ell$ is due to mode loss to the timestream filtering and matrix purification.
• Figure 4: Results of a multicomponent multi-spectral likelihood analysis of BICEP2/ Keck+WMAP/ Planck data. The red faint curves are the baseline result from the previous BK14 paper (the black curves from Fig. 4 of that paper). The bold black curves are the new baseline BK15 results. Differences between these analyses include adding Keck Array data taken during the 2015 observing season, in particular doubling the 95 GHz sensitivity and adding, for the first time, a 220 GHz channel. (In addition the $\epsilon$ prior is modified.) The upper limit on the tensor-to-scalar ratio tightens to $r_{0.05}<0.072$ at 95% confidence. The parameters $A_\mathrm{d}$ and $A_\mathrm{sync}$ are the amplitudes of the dust and synchrotron $B$-mode power spectra, where $\beta$ and $\alpha$ are the respective frequency and spatial spectral indices. The correlation coefficient between the dust and synchrotron patterns is $\epsilon$. In the $\beta$, $\alpha$ and $\epsilon$ panels the dashed lines show the priors placed on these parameters (either Gaussian or uniform). Broadening or tightening the uniform prior range on $\alpha_s$ and $\alpha_d$ results in very small changes, and negligible changes to the $r$ constraint.
• Figure 5: Constraints in the $r$ vs. $n_s$ plane when using Planck 2015 plus additional data, and when also adding BICEP2/ Keck data through the end of the 2015 season---the constraint on $r$ tightens from $r_{0.05}<0.12$ to $r_{0.05}<0.06$. This figure is adapted from Fig. 21 of Ref. planck2015XIII, with two notable differences: switching lowP to lowT plus a $\tau$ prior of $0.055\pm 0.009$ Ref. planck2016XLVI, and the exclusion of JLA data and the $H_0$ prior.