[1403.3985] BICEP2 I: Detection Of B-mode Polarization at Degree Angular Scales

Authors:  BICEP2 Collaboration: P. A. R Ade (1), R. W. Aikin (2), D. Barkats (3), S. J. Benton (4), C. A. Bischoff (5), J. J. Bock (2,6), J. A. Brevik (2), I. Buder (5), E. Bullock (7), C. D. Dowell (6), L. Duband (8), J. P. Filippini (2), S. Fliescher
Abstract:  We report results from the BICEP2 experiment, a Cosmic Microwave Background (CMB) polarimeter specifically designed to search for the signal of inflationary gravitational waves in the B-mode power spectrum around l=80. The telescope comprised a 26 cm aperture all-cold refracting optical system equipped with a focal plane of 512 antenna coupled transition edge sensor (TES) 150 GHz bolometers each with temperature sensitivity of approx. 300 uk.sqrt(s). BICEP2 observed from the South Pole for three seasons from 2010 to 2012. A low-foreground region of sky with an effective area of 380 square degrees was observed to a depth of 87 nK-degrees in Stokes Q and U. In this paper we describe the observations, data reduction, maps, simulations and results. We find an excess of B-mode power over the base lensed-LCDM expectation in the range 30<l<150, inconsistent with the null hypothesis at a significance of $>5\sigma$. Through jackknife tests and simulations based on detailed calibration measurements we show that systematic contamination is much smaller than the observed excess. We also estimate potential foreground signals and find that available models predict these to be considerably smaller than the observed signal. These foreground models possess no significant cross-correlation with our maps. Additionally, cross-correlating BICEP2 against 100 GHz maps from the BICEP1 experiment, the excess signal is confirmed with $3\sigma$ significance and its spectral index is found to be consistent with that of the CMB, disfavoring synchrotron or dust at $2.3\sigma$ and $2.2\sigma$, respectively. The observed B-mode power spectrum is well-fit by a lensed-LCDM + tensor theoretical model with tensor/scalar ratio $r=0.20^{+0.07}_{-0.05}$, with r=0 disfavored at $7.0\sigma$. Subtracting the best available estimate for foreground dust modifies the likelihood slightly so that r=0 is disfavored at $5.9\sigma$.
[PDF]  [PS]  [BibTex]  [Bookmark]

Discussion related to specific recent arXiv papers
Post Reply
Antony Lewis
Posts: 1936
Joined: September 23 2004
Affiliation: University of Sussex
Contact:

[1403.3985] BICEP2 I: Detection Of B-mode Polarization at D

Post by Antony Lewis » March 19 2014

What an amazing result, which raises a lot of interesting questions!

The high-looking points at [tex]l > 150[/tex] are obviously a concern: without using Keck array data the shape of expected primordial tensor B spectrum is not there. Taken at face value (which the authors warn against) it seems to favour a lot more power on sub-horizon scales than you'd expect from either primordial tensors or lensing. Combined with Planck's preference for a low TT plateau, which disfavours a large tensor amplitude on large scales, a naive application of the Planck+WP+BICEP likelihood to LCDM+[tex]r[/tex]+[tex]n_T[/tex] seems to favour a very non-inflationary [tex]n_T>0[/tex]:

Image

Of course this should not be taken too seriously, since the Keck cross-correlation (their Fig 9) does seem to indicate something more like the expected shape at higher [tex]l[/tex]s, increasing overall confidence in the result and interpretation. The Keck points look systematically lower, which might indicate ultimately something closer to [tex]r\approx 0.16[/tex], and consistent with near scale-invariance. Hopefully more data from Keck will pin this down one way or another in relatively short order.

Even if you assume an inflationary consistency relation for [tex]n_T[/tex] so that it is nearly scale-invariant, the paper points out that the combination with the preference for low power on large scales in Planck then drives a preference for significantly lower scalar power on large scales than you'd expect from a nearly constant [tex]n_s[/tex] vanilla LCDM model. For example allowing running Planck+WP+BICEP gives [tex]n_{run}\approx -0.03\pm 0.02[/tex] (95%):

Image

Of course anything else that can lower the large-scale scalar TT can also help to improve consistency, for example a sterile neutrinos allows for bluer [tex]n_s[/tex] and hence lower large-scale scalar power:

Image

([tex]N_{eff}\approx 4[/tex] gives roughly scale invariant [tex]n_s[/tex] and fits about as well as running, though it is hard to motivate an extra sterile neutrino that thermalized at about the same temperature as the active neutrinos, and it is unlikely to improve joint fits on its own)

I look forward to seeing papers on inflationary transients that only strongly affect large scales (e.g. from bubble nucleation, fast-roll, etc) very shortly...

Please note above plots are just a preliminary quick-look and come with usual health warnings.

Sesh Nadathur
Posts: 3
Joined: March 19 2014
Affiliation: University of Helsinki

[1403.3985] BICEP2 I: Detection Of B-mode Polarization at De

Post by Sesh Nadathur » March 19 2014

When fitting to BICEP data do you use bandpowers 1-5 only (as in their paper), or all bandpowers 1-9?

Subodh Patil
Posts: 4
Joined: March 20 2014
Affiliation: Leiden University
Contact:

[1403.3985] BICEP2 I: Detection Of B-mode Polarization at De

Post by Subodh Patil » March 20 2014

Lots of interesting questions indeed! Although the cross correlation with Keck data does show more of the expected shape at high l's, doesn't it still point to a lack of power at lower l's (c.f. fig 9 of the above reference), suggesting that a preference for a blue tilt would persist in the B2 x Keck cross correlation?

I understand as a naive theorist that there may be lots of remaining foreground issues, uncertainties from polarized dust, EE contamination, sods law etc to be sorted out, but just to put it out there-- there are some models of structure formation that can easily generate a large tensor to scalar ratio, which moreover, have predicted a blue tensor tilt that is complimentary to the scalar tilt--

http://arxiv.org/abs/1403.4927
http://arxiv.org/abs/hep-th/0604126

Antony Lewis
Posts: 1936
Joined: September 23 2004
Affiliation: University of Sussex
Contact:

Re: [1403.3985] BICEP2 I: Detection Of B-mode Polarization a

Post by Antony Lewis » March 20 2014

This was using E and B from all 9 bins. Non-[tex]n_T[/tex] joint constraints don't change very much if you restrict to only five bins or use only B.

If you very crudely replace the BICEP2 points in the likelihood with the BICEP2xKECK points from the figure (thanks Adam Moss), assuming window, covariance, etc, are unchanged, and use just BB, things do change a bit more. For example the [tex]n_T[/tex] plot then looks like this:

Image

Here much of the pull to higher [tex]n_T[/tex] is presumably coming from the tension with Planck lowl-[tex]l[/tex] rather than the BB spectrum directly, though at least it now looks marginally consistent with scale invariant.

Although the cross-likelihood is a fudge, it does seem to be a reasonable handle on the size of non-foreground systematic errors; the value of [tex]r[/tex] shifts down a bit, as also shown here as the contour superimposed on the BICEP2+Planck sample points:

Image

So the Planck+BICEP+running constraint shifts from [tex]r_{0.002}=0.22\pm 0.11[/tex] (95%) to about [tex]r_{0.002}=0.14\pm 0.09[/tex], though still preferring running at over 2\sigma. Note that in running models [tex]r_{0.002}[/tex] gets larger because the large-scale scalar power is decreasing even for the same tensor amplitude; the corresponding numbers at smaller scales are [tex]r_{0.05}=0.2\pm0.08[/tex] and [tex]r_{0.05} = 0.13\pm 0.8[/tex].

Allowing for foregrounds could reduce this by a further 0.04 or so in the BICEP models, taking the signal to a level r\approx 0.1, which sits a bit more comfortably with Planck in simplest inflation scenarios, though tensions in the low-[tex]l[/tex] platau would remain.

Sesh Nadathur
Posts: 3
Joined: March 19 2014
Affiliation: University of Helsinki

[1403.3985] BICEP2 I: Detection Of B-mode Polarization at De

Post by Sesh Nadathur » March 20 2014

The BICEP team are explicitly cautioning against drawing any conclusions from the Keck data, which are still preliminary. For instance, although one could use the B2xKeck spectrum to speculate that with KeckxKeck the high \ell points might come down, at low \ell B2xKeck appears to be consistent with or lower than B2xB1 - and of course B1xB1 shows no detection of \r.

But if the preference for blue \n_T comes as much from Planck low-\ell TT as from any of the high-\ell BB points, then perhaps any detection of sizable \r would also favour a blue tilt.

Sesh Nadathur
Posts: 3
Joined: March 19 2014
Affiliation: University of Helsinki

[1403.3985] BICEP2 I: Detection Of B-mode Polarization at De

Post by Sesh Nadathur » March 20 2014

Actually, my naive speculations were probably completely off-base because the Keck error bars will be so much smaller than B1 ...

Subodh Patil
Posts: 4
Joined: March 20 2014
Affiliation: Leiden University
Contact:

[1403.3985] BICEP2 I: Detection Of B-mode Polarization at De

Post by Subodh Patil » March 20 2014

Thanks, Antony. This is incredibly informative!

Early days yet, but just to comment specifically on the `early transient' models of inflation that exploit initial time fast roll-- you'd need a trajectory that stays off attractor for long enough, or a very tuned (and transiently non-flat) potential to get the requisite suppression of power at in large scales. Both of these are in tension with the shift symmetry seemingly required by a the super-Planckian field excursion over the scales we're seeing in the CMB that r = 0.2 implies.

Another naive question to anyone who may know: my understanding was that if r turned out to be of the order of 0.01 - 0.001, we would have needed a PRISM like observation to have meaningfully measured a tilt to anywhere close the accuracy with which we measure the scalar tilt. Now that the indications are that r is at least ten times larger than than, could ground based observations get ever get to the stage where their measurements of the tensor tilt can be comparable to the measurement of the scalar tilt? i.e. if actually there, could one hope even in principle to infer the tensor spectrum to be > 0 with greater than at least 3 sigma confidence, say, from the ground? what are ACTPol and SPTPol's sensitivities in light of r possibly being ~ 0.2?

Antony Lewis
Posts: 1936
Joined: September 23 2004
Affiliation: University of Sussex
Contact:

Re: [1403.3985] BICEP2 I: Detection Of B-mode Polarization a

Post by Antony Lewis » March 20 2014

The tensor tilt is limited by cosmic variance for large [tex]r[/tex] values, and the BB spectrum falls rapidly at [tex]l>150[/tex]. This gives a significantly smaller leverarm than with TT and larger cosmic variance on the smallest scales you can see: the ideal measurement of [tex]n_T[/tex] from the CMB alone would be like measuring [tex]n_s[/tex] from Planck restricted roughly to [tex]2 \le l \le 200[/tex] (without parameter degeneracies). So the best you could do is not good enough to distinguish inflation models that differ only by [tex]\Delta n_T = O(|n_s-1|)[/tex]. See astro-ph/0305411.

Antony Lewis
Posts: 1936
Joined: September 23 2004
Affiliation: University of Sussex
Contact:

Re: [1403.3985] BICEP2 I: Detection Of B-mode Polarization a

Post by Antony Lewis » March 20 2014

For people interested in how parameters move around with different parameter and data combinations, and a more detailed comparison of shifts moving from BICEP2 to fudged BICEP2xKECK BB, here are provisional parameter constraints from a grid of models (including breakdowns of best-fit chi-squareds) as output from CosmoMC March 2014 version:
bicep_planck.pdf

Anze Slosar
Posts: 183
Joined: September 24 2004
Affiliation: Brookhaven National Laboratory
Contact:

[1403.3985] BICEP2 I: Detection Of B-mode Polarization at De

Post by Anze Slosar » March 21 2014

Ok, before I die by tables, is it right that the bottom line is that surprisingly little changes if you add bicep, with or without running? The usual low-z suspects like Om, sigma8 and h remain pretty much unchanged... In other words, low redshift probes will have hard time telling these scenarios apart.

Subodh Patil
Posts: 4
Joined: March 20 2014
Affiliation: Leiden University
Contact:

[1403.3985] BICEP2 I: Detection Of B-mode Polarization at De

Post by Subodh Patil » March 21 2014

Thanks for that reference, Antony. I wasn't aware of it. If I understand their detectability limits correctly (specifically their fig 3), then it appears that if r > 0.1, one can expect |\Delta n_T| \gtrsim 0.01 to 0.03 as one ranges the angular resolution of the observation from 1' to 5' and \Delta_p = \omega^{-1/2} (where \omega is the weight per solid angle of the stokes parameters) ranges from 15\sqrt{2} to \sqrt{2} \mu K...

This is clearly a very interesting range as far as testing various universality classes of inflationary realizations/ alternatives to inflation are concerned. I'm curious how the various planned/ in progress ground/ space based polarization observations stack up wrt those specs...

Antony Lewis
Posts: 1936
Joined: September 23 2004
Affiliation: University of Sussex
Contact:

Re: [1403.3985] BICEP2 I: Detection Of B-mode Polarization a

Post by Antony Lewis » March 21 2014

Thanks for pointing that out, I think I mis-remembered the details. The error on [tex]n_T[/tex] will always be much worse than the error on [tex]n_s[/tex] for sure if you measure to small scales, but as you say it looks as though you can still in principle just get something potentially useful for large [tex]r\approx 0.1[/tex].

Antony Lewis
Posts: 1936
Joined: September 23 2004
Affiliation: University of Sussex
Contact:

Re: [1403.3985] BICEP2 I: Detection Of B-mode Polarization a

Post by Antony Lewis » March 22 2014

There's another interesting point regarding tensor tilt (pointed out to me by Antony Challinor), in that the shape of the large-scale BB spectrum is close white. This comes from the blue spectrum of the sources, so that the low-[tex]L[/tex] (but above reionization) scales in BB are dominated by signiciantly smaller horizon-size modes. The close-to-white shape tells you very little about the spectrum of the underlying sources except that they are quite blue (the same thing happens for the [tex]l<1000[/tex] white lensing BB, and also large-scale 21cm/number count angular power spectra with sharp window functions). The large values of [tex]n_T[/tex] seen in the above plots therefore work better than you might expect, since the actual range of scales being probed by the BICEP data points is significantly narrowed that you you might expect.

This is illustrated in the following figure of the (normalized) tensor-mode transfer functions

Image

This shows where in [tex]k[/tex] each observable [tex]C_l[/tex] is coming from. Vertical lines denote [tex]l = k\chi_*[/tex], and the bottom panel is the middle paanel (for exact L values L=45,73,109,144) multiplied by the BICEP window functions. As you can see at the same [tex]L<100[/tex], the TT tensor constributions come from a significantly larger scale than the BB contributions, which bunch up around 0.01/Mpc. You can therefore significantly lower the tensor contributions to Planck TT while keeping the BICEP data points roughly the same, by putting most of the tensor power around 0.01/Mpc.

Gansukh Tumurtushaa
Posts: 34
Joined: October 05 2013
Affiliation: Sogang univ.

[1403.3985] BICEP2 I: Detection Of B-mode Polarization at De

Post by Gansukh Tumurtushaa » August 15 2014

Dear all,

I have little bit confused with a difference of the following two data,
first one is a result of BICEP2(lets say A) for r_{0.002} vs n_{s} which shown here, http://bicepkeck.org/B2_2014_i_figs/running_rvsns.pdf, and
second one is the one posted above, r_{0.002} vs n_{s} plot (lets say B) with combinations of green,+grey, +red, and +blue).

The problem is that on both blue contours of A and the grey contours of B shows the same r_{0.002} vs n_{s} plot of Planck+WP+highL+BICEP2 data with running. But if you look carefully at the contours, then you will see that the r_{0.002} values are different from one another. Therefore my question is simply why they are different while data sets are same for both? Does anyone can explain me? that would be highly appreciated.

Another simple question is that what does the 0.002 or 0.5 the subscript of r imply? I think they are the value of k_0 which appear in the denominator of the expression of the power spectrum like P_{s}=A_{s}(k/k_{0})^{n_{s}-1}. Is that right?

Thank you in advance.

Post Reply