Issue 
A&A
Volume 647, March 2021



Article Number  A128  
Number of page(s)  18  
Section  Cosmology (including clusters of galaxies)  
DOI  https://doi.org/10.1051/00046361/202039585  
Published online  19 March 2021 
Planck constraints on the tensortoscalar ratio
^{1}
Université ParisSaclay, CNRS/IN2P3, IJCLab, 91405 Orsay, France
email: tristram@lal.in2p3.fr
^{2}
Université de Toulouse, UPSOMP, IRAP, 31028 Toulouse cedex 4, France
^{3}
CNRS, IRAP, 9 Av. colonel Roche, BP 44346, 31028 Toulouse cedex 4, France
^{4}
Jet Propulsion Laboratory, California Institute of Technology, 4800 Oak Grove Drive, Pasadena, California, USA
^{5}
Warsaw University Observatory, Aleje Ujazdowskie 4, 00478 Warszawa, Poland
^{6}
Computational Cosmology Center, Lawrence Berkeley National Laboratory, Berkeley, California, USA
^{7}
Space Sciences Laboratory, University of California, Berkeley, California, USA
^{8}
Institute of Theoretical Astrophysics, University of Oslo, Blindern, Oslo, Norway
^{9}
Instituto de Física de Cantabria (CSICUniversidad de Cantabria), Avda. de los Castros s/n, Santander, Spain
^{10}
Haverford College Astronomy Department, 370 Lancaster Avenue, Haverford, Pennsylvania, USA
^{11}
Department of Physics & Astronomy, University of British Columbia, 6224 Agricultural Road, Vancouver, British Columbia, Canada
Received:
2
October
2020
Accepted:
21
December
2020
We present constraints on the tensortoscalar ratio r using Planck data. We use the latest release of Planck maps, processed with the NPIPE code, which produces calibrated frequency maps in temperature and polarisation for all Planck channels from 30 GHz to 857 GHz using the same pipeline. We computed constraints on r using the BB angular power spectrum, and we also discuss constraints coming from the TT spectrum. Given Planck’s noise level, the TT spectrum gives constraints on r that are cosmicvariance limited (with σ_{r} = 0.093), but we show that the marginalised posterior peaks towards negative values of r at about the 1.2σ level. We derived Planck constraints using the BB power spectrum at both large angular scales (the ‘reionisation bump’) and intermediate angular scales (the ‘recombination bump’) from ℓ = 2 to 150 and find a stronger constraint than that from TT, with σ_{r} = 0.069. The Planck BB spectrum shows no systematic bias and is compatible with zero, given both the statistical noise and the systematic uncertainties. The likelihood analysis using B modes yields the constraint r < 0.158 at 95% confidence using more than 50% of the sky. This upper limit tightens to r < 0.069 when Planck EE, BB, and EB power spectra are combined consistently, and it tightens further to r < 0.056 when the Planck TT power spectrum is included in the combination. Finally, combining Planck with BICEP2/Keck 2015 data yields an upper limit of r < 0.044.
Key words: cosmology: observations / cosmic background radiation / cosmological parameters / gravitational waves / methods: data analysis
© M. Tristram et al. 2021
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1. Introduction
While the idea of cosmic inflation was introduced about 40 years ago to solve inherent problems with the canonical hot bigbang model (Brout et al. 1978; Starobinsky 1980; Kazanas 1980; Sato 1981; Guth 1981; Linde 1982, 1983; Albrecht & Steinhardt 1982), attention quickly focused on using it as a means to generate cosmological perturbations from quantum fluctuations (Mukhanov & Chibisov 1981, 1982; Hawking 1982; Guth & Pi 1982; Starobinsky 1982; Bardeen et al. 1983; Mukhanov 1985). These perturbations include a tensor component (i.e. gravitational waves) as well as the scalar component (i.e. density variations). Inflationary gravitational waves entering the horizon between the epoch of recombination and the present day generate a tensor contribution to the largescale cosmic microwave background (CMB) anisotropy. Hence, primordial tensor fluctuations contribute to the CMB anisotropies, both in temperature (T) and in polarisation (E and B modes; Seljak 1997; Kamionkowski et al. 1997; Seljak & Zaldarriaga 1997).
As described in Planck Collaboration VI (2020) and Planck Collaboration X (2020), the comoving wavenumbers of tensor modes probed by the CMB temperature anisotropy power spectrum have k ≲ 0.008 Mpc^{−1}, with very little sensitivity to higher wavenumbers because gravitational waves decay on subhorizon scales. The corresponding multipoles in the harmonic domain are ℓ ≲ 100, for which the scalar perturbations dominate with respect to tensor modes in temperature. The tensor component can be fitted together with the scalar one, and the precision of the Planck constraint is limited by the cosmic variance of the largescale anisotropies.
In polarisation, the EE and TE spectra also contain a tensor signal coming from the lastscattering and reionisation epochs. The BB power spectrum, however, is treated differently when determining the tensor contribution, since the model does not predict any primordial scalar fluctuations in BB. As a consequence, a primordial Bmode signal would be a direct signature of tensor modes. However, depending on the amplitude of the tensortoscalar ratio, such a signal may be masked by Emode power that is transformed to Bmode power through lensing by gravitational potentials along the line of sight (socalled ‘BB lensing’, Zaldarriaga & Seljak 1998). BB lensing has been measured with high accuracy by Planck in both harmonic (Planck Collaboration VIII 2020) and map (Planck Collaboration Int. XLI 2016) domains, as well as by groundbased observatories POLARBEAR (POLARBEAR Collaboration 2017), SPTpol (Sayre et al. 2020), and ACTPol (Choi et al. 2020). But a primordial BB tensor signal has not been detected yet.
The scalar and tensor CMB angular power spectra are plotted in Fig. 1 for the Planck 2018 cosmology and for two values of the tensortoscalar ratio, namely r = 0.1 and r = 0.01. For a further discussion of the tensortoscalar ratio and its implications for inflationary models, see Planck Collaboration XXII (2014), Planck Collaboration XX (2016), and Planck Collaboration X (2020). We note that the signal from tensor modes in EE is similar to that in BB modes, which makes EE (in particular at low multipoles) an important data set for tensor constraints. Indeed, limits set by cosmic variance alone for fullsky spectra are σ_{r}(TT) = 0.072, σ_{r}(EE) = 0.023, and σ_{r}(BB) = 0.000057 for r = 0. In this paper, we make use of a polarised EB likelihood, which consistently includes the correlated polarisation fields E and B, and covers the range of multipoles where tensor modes can be constrained using Planck data (i.e. from ℓ = 2 to ℓ = 150).
Fig. 1. Scalar (thick solid lines) versus tensor spectra for r = 0.1 (dashed lines) and r = 0.01 (dotted lines). Spectra for TT are in black, EE in blue, and BB in red. The red solid line corresponds to the signal from BB lensing. 
At present the tightest Bmode constraints on r come from the BICEP/Keck measurements (BK15; BICEP2 Collaboration 2018), which cover approximately 400 deg^{2} centred on RA = 0^{h}, Dec = −57° .5. These measurements probe the peak of the Bmode power spectrum at around ℓ = 100, corresponding to gravitational waves with k ≈ 0.01 Mpc^{−1} that enter the horizon during recombination (i.e. somewhat smaller than the scales contributing to the Planck temperature constraints on r). The results of BK15 give a limit of r < 0.07 at 95% confidence, which tightens to r < 0.06 in combination with Planck temperature and other data sets.
Planck Collaboration V (2020) presented Planck Bmode constraints from the 100 and 143GHz HFI channels with a 95% upper limit of r < 0.41 (at a somewhat larger pivot scale, as described in the next section), using only a limited number of multipoles around the socalled ‘reionisation bump’ (2 ≤ ℓ ≤ 29). Using Planck NPIPE maps (Planck Collaboration Int. LVII 2020), called Planck Release 4 (PR4), we are now able to constrain the BB power spectrum for a much larger number of modes, including both the reionisation bump at large angular scales (ℓ ≲ 30) and the socalled ‘recombination bump’ at intermediate scales (50 ≲ ℓ ≲ 150). In this paper, we first describe, in Sect. 2, the cosmological model used throughout the analysis. We then detail the data and the likelihoods in Sect. 3. Section 4 focuses on constraints from TT and in particular the impact of the lowℓ data in temperature. Section 5 gives constraints from the BB angular power spectrum using Planck data, while results from the full set of polarisation power spectra are given in Sect. 6. In Sect. 7, we combine all data sets to provide the most robust constraints on BB coming from Planck and in combination with other CMB data sets, such as the results from the BICEP/Keck Collaboration. Finally, we provide details of several parts of our analysis in a set of appendices, specifically describing the transfer function for BB, the HiLLiPoP likelihood, largescale polarised power spectra, the crossspectrum correlation matrix, comparison between PR3 and PR4, robustness tests, triangle plots for ΛCDM+r parameters, and comparison with other BB spectrum measurements.
2. Cosmological model
We use the baseΛCDM model, which has been established over the last couple of decades to be the simplest viable cosmological model, in particular with the Planck results (e.g. Planck Collaboration VI 2020). In this model, we assume purely adiabatic, nearly scaleinvariant perturbations at very early times, with curvaturemode (scalar) and tensormode power spectra parameterised by
where A_{s} and A_{t} are the initial superhorizon amplitudes for curvature and tensor perturbations, respectively. The primordial spectral indexes for scalar (n_{s}) and tensor (n_{t}) perturbations are taken to be constant. This means that we assume no ‘running’, i.e. a pure powerlaw spectrum with dn_{s}/dlnk = 0. We set the pivot scale at k_{0} = 0.05 Mpc^{−1}, which roughly corresponds to approximately the middle of the logarithmic range of scales probed by Planck; with this choice, n_{s} is not strongly degenerate with the amplitude parameter A_{s}. Note that for historical reasons, the definitions of n_{s} and n_{t} differ, so that a scaleinvariant scalar spectrum corresponds to n_{s} = 1, while a scaleinvariant tensor spectrum corresponds to n_{t} = 0.
The latetime parameters, on the other hand, determine the linear evolution of perturbations after they reenter the Hubble radius. We use the basis (Ω_{b}h^{2}, Ω_{c}h^{2}, θ_{*}, τ) following the approach in Planck cosmological studies (Planck Collaboration VI 2020), where Ω_{b}h^{2} is the baryon density today, Ω_{c}h^{2} is the cold dark matter (CDM) density today, θ_{*} is the observed angular size of the sound horizon at recombination, and τ is the reionisation optical depth.
The amplitude of the smallscale linear CMB power spectrum is proportional to A_{s}e^{−2τ}. Because Planck measures this amplitude very accurately, there is a tight linear constraint between τ and lnA_{s}. For this reason, we usually adopt lnA_{s} as a base parameter with a flat prior; lnA_{s} has a significantly more Gaussian posterior than A_{s}. A linear parameter redefinition then allows the degeneracy between τ and A_{s} to be explored efficiently. Note that the degeneracy between τ and A_{s} is broken by the relative amplitudes of largescale temperature and polarisation CMB anisotropies and by the effect of CMB lensing.
We define r ≡ A_{t}/A_{s}, the primordial tensortoscalar ratio defined explicitly at the scale k_{0} = 0.05, Mpc^{−1}. Our constraints are only weakly sensitive to the tensor spectral index, n_{t}. We adopt the singlefieldinflation consistency relation n_{t} = −r/8. Note that the Planck Collaboration also discussed r constraints for k_{0} = 0.002, Mpc^{−1} (Planck Collaboration V 2020). Given the definitions in Eqs. (1) and (2), the tensortoscalar ration scales with (0.05/0.002)^{−r/8}, which means that r_{0.002} is lower by 4% at r ≃ 0.1 compared to r_{0.05} and less than 0.4% lower for r < 0.01.
In this work, we use an effective tensortoscalar ratio r_{eff}, which we extend into the negative domain by modifying the Boltzmannsolver code CLASS (Blas et al. 2011). While negative tensor amplitudes are unphysical, this approach will allow us to derive posteriors without boundaries, facilitating detection of potential biases, and enabling us to determine a more accurate statistical definition of the constraints on r. With r_{eff} we are able to independently discuss both the uncertainty of r (σ_{r}) and corresponding upper limits (depending on the maximum a posteriori probability). In the rest of this paper, we simply write r as the effective tensortoscalar ratio, and report upper limits for positive tensor amplitudes, for which r_{eff} = r. We use 95% confidence levels when reporting upper limits, and a 68% confidence interval with the maximum a posteriori probability.
3. Data and likelihoods
3.1. Data and simulations
The sky measurements used in this analysis are the PR4 maps available from the Planck Legacy Archive^{1} (PLA) and from the National Energy Research Scientific Computing Center (NERSC)^{2}. They have been produced with the NPIPE processing pipeline, which creates calibrated frequency maps in temperature and polarisation from the Planck Low Frequency Instrument (LFI) and High Frequency Instrument (HFI) data. As described in Planck Collaboration Int. LVII (2020), NPIPE processing includes several improvements, resulting in lower levels of noise and systematics in both frequency and componentseparated maps at essentially all angular scales, as well as notably improved internal consistency between the various frequencies.
NPIPE achieves an overall lower noise level in part by incorporating the data acquired during the 4minute spacecraft repointing manoeuvres that take place between the 30to70min stable science scans. Residual systematics are suppressed using a mitigation strategy that combines aspects of both LFI and HFI processing pipelines. Most importantly, gain fluctuations, bandpass mismatch, and other systematics are formulated into timedomain templates that are fitted and subtracted as a part of the mapmaking process. Degeneracies between sky polarisation and systematic templates are broken by constructing a prior of the polarised foreground sky using the extreme polarisationsensitive frequencies (30, 217, and 353 GHz).
Moreover, the PR4 release comes with 400 simulations of signal, noise, and systematics, componentseparated into CMB maps, which allow for an accurate characterisation of the noise and systematic residuals in the Planck maps. This is important because Planck polarisation data are cosmicvariancedominated only for a few multipoles at very large scales in EE (ℓ < 8, as shown in Fig. 2). These simulations, even though limited in number, represent a huge effort in terms of CPU time. They are essential in order to compute the following two additional quantities.
Fig. 2. Variances for crossspectra in EE and BB based on PR4 simulations, including: cosmic (sample) variance (black); analytic statistical noise (red); and PR4 noise from Monte Carlo simulations (green), including noise and systematics with (solid line) or without (dashed line) correction for the transfer function. The sky fraction used here is 80%, as it illustrates well the effect of both systematics and transferfunction corrections (Planck Collaboration Int. LVII 2020). 
First, the endtoend transfer function from the data reduction (including TOI processing, mapmaking and template fitting for mitigation of systematics, component separation, and powerspectrum estimation). The transfer function is defined as the ratio between the estimated output power spectrum and the input one, averaged over all the simulations (see Sect. 4.3 of Planck Collaboration Int. LVII 2020, for the details).
Second, the covariance of the data (here we use the crosspower spectra), which is the only way to propagate uncertainties when those are dominated by systematics (from the instrument or from foregrounds).
Note that these two quantities estimated from the simulations are directly related to two different characteristics of the final parameter posteriors: the bias of the mean (the transfer function); and the width of the posterior (as propagated into parameter constraints by the covariance matrix in the likelihood). They can be separated from each other, meaning that one systematic effect can easily produce a significant bias without any strong impact on the variance, while another effect can produce a large increase of the variance with no associated bias.
The NPIPE simulations include the systematic effects relevant for polarisation studies, specifically analoguetodigitalconverter nonlinearities, gain fluctuations, bandpass mismatch between detectors, correlated noise (including 4K line residuals), and fullbeam convolutions for each detector.
The use of a polarisation prior in NPIPE processing causes a suppression of largescale (ℓ < 20) CMB polarisation, which needs to be corrected. As explained in Planck Collaboration Int. LVII (2020), allowing for a nontrivial transfer function is a compromise between measuring very noisy but unbiased largescale polarisation from all lowℓ modes, and filtering out the modes that are most affected by the calibration uncertainties left in the data by the Planck scan strategy. As detailed in Planck Collaboration Int. LVII (2020), the transfer function to correct for this bias is determined from simulations. It is then used to correct the power spectrum estimates, just as instrumental beam and pixel effects must be deconvolved. Due to the fact that E modes dominate the CMB polarisation, the simulations do not yield a definitive measurement of the Bmode transfer function. We have chosen to conservatively deconvolve the Emode transfer function from the Bmode spectrum in order to provide a robust upper limit on the true Bmode spectrum. Indeed, when regressing the templates fitted during the mapmaking process with pure E and B CMB maps, we found a similar impact on the EE and BB power spectra (see Appendix A). Moreover, in the situation where primordial Bmode power is not detected, the transfer function correction essentially increases the variance estimate at low multipoles, which propagates the uncertainty induced by the degeneracy between the sky and the systematic templates used in NPIPE. Note that this uncertainty is small compared to the impact of systematics in the error budget (see Fig. 2).
To compute unbiased estimates of the angular power spectra, we perform crosscorrelations of two independent splits of the data. As shown in Planck Collaboration Int. LVII (2020), the most appropriate split for the Planck data is represented by the detectorset (hereafter ‘detset’) maps, comprising two subsets of maps at each frequency, with nearly independent noise characteristics, made by combining half of the detectors. This was obtained by processing each split independently, in contrast to the detset maps produced in the previous Planck releases. Note that timesplit maps (made from, e.g. ‘oddeven rings’ or ‘halfmission data’) share the same instrumental detectors, and therefore exhibit noise correlations due to identical spectral bandpasses and optical responses. The use of timesplit maps is subject to systematic biases in the crosspower spectra (see Sect. 3.3.3 in Planck Collaboration V 2020), as well as underestimation of the noise properties in computing the halfdifferences (which must be compensated by a rescaling of the noise in the PR3 as described in Appendix A.7 of Planck Collaboration III 2020). Hence we use detset splits here.
Uncertainties at the powerspectrum level are dominated by noise and systematics, as illustrated in Fig. 2. Thanks to the NPIPE processing, we are now able to show the impact of the systematics at low ℓ. This is illustrated by comparing the PR4 endtoend noise (based on the Monte Carlo simulations, including instrumental noise, systematics, and foreground uncertainties, and corrected for the transfer function both in EE and BB) with the propagation of the statistical noise coming from the analytic pixelpixel covariance matrix. The systematic uncertainties dominate at ℓ ≲ 15, then slowly decrease so that the effective uncertainties converge towards the analytic estimate at higher multipoles.
3.2. polarised sky masks
Foreground residuals in the foregroundcleaned maps dominate the polarised CMB signal near the Galactic plane. To avoid contamination from these residuals in the cosmological analysis, we mask the Galactic plane. We use a series of different retained sky fractions (from 30% to 70%) to check the consistency of our results with respect to foreground residuals (Fig. 3).
Fig. 3. Galactic masks used for the Planck likelihoods. The mask shown in dark blue indicates the sky rejected in order to retain a 70% sky fraction for analysis. The masks shown in light blue, green, orange and red incrementally omit further parts of the sky, corresponding in turn to 60, 50, 40 and 30% retained sky fractions, the latter shown in white. 
The masks used in this analysis are a combination of a mask for polarisation intensity (to avoid polarised foreground residuals), a mask for total intensity (to avoid potential temperaturetopolarisation leakage residuals), and the confidence mask for component separation provided by the Planck Collaboration. The intensity mask is obtained by thresholding the combination of the 353GHz intensity map (which traces dust) scaled to 143 GHz, and the 30GHz intensity map (which traces synchrotron) scaled to 100 GHz. The polarisation map is constructed similarly. Both foreground tracers are smoothed beforehand with a 10° Gaussian window function.
The impact of the emission of extragalactic polarised sources on the power spectra is negligible, given the Planck resolution and noise level. The confidence mask for component separation ensures the masking of the strongest sources, which could also produce residuals through temperaturetopolarisation leakage.
3.3. Likelihoods
Table 1 summarises the likelihoods used in this analysis, which are described below.
Summary of the likelihoods used in this paper.
3.3.1. Lowℓ temperature likelihood
We use the Planck public lowℓ temperatureonly likelihood based on the PR3 CMB map recovered from the componentseparation procedure (specifically Commander) described in detail in Planck Collaboration V (2020). At large angular scales, Planck temperature maps are strongly signaldominated, and there is no expected gain in updating this likelihood with the PR4 data.
As discussed in Planck Collaboration XX (2016), the lowℓ temperature data from Planck have a strong impact on the r posterior and the derivation of the corresponding constraints. This is because the deficit of power in the measured C_{ℓ}s at lowℓ in temperature (see the discussions in Planck Collaboration XVI 2014 and Planck Collaboration Int. LI 2017) lowers the probability of tensor models, which ‘add’ power at low multipoles. This shifts the maximum in the posterior of r towards low values (or even negative values when using r_{eff}, as we show in Sect. 4).
3.3.2. Highℓ likelihood
At small angular scales (ℓ > 30), we use the HiLLiPoP likelihood, which can include the TT, TE, and/or EE power spectra. HiLLiPoP has been used as an alternative to the public Planck likelihood in the 2013 and 2015 Planck releases (Planck Collaboration XV 2014; Planck Collaboration XI 2016), and is described in detail in Couchot et al. (2017a). In this paper, the HiLLiPoP likelihood is applied to the PR4 detset maps at 100, 143, and 217 GHz. We focus on the TT spectra, since there is marginal additional information at small scales in TE or EE for tensor modes, due to Planck noise. We only make use of TE in Sect. 7 in order to help constrain the spectral index n_{s}. The likelihood is a spectrumbased Gaussian approximation, with semianalytic estimates of the C_{ℓ} covariance matrix based on the data. The crossspectra are debiased from the effects of the mask and the beam leakage using Xpol (a generalisation to polarisation of the algorithm presented in Tristram et al. 2005^{3}) before being compared to the model, which includes CMB and foreground residuals. The beam window functions are evaluated using QUICKPOL (Hivon et al. 2017), adapted to the PR4 data. These adaptations include an evaluation of the beamleakage effect, which couples temperature and polarisation modes due to the beam mismatch between individual detectors.
The model consists of a linear combination of the CMB power spectrum and several foregrounds residuals. These are:
– Galactic dust (estimated directly from the 353GHz channel);
– The cosmic infrared background (as measured in Planck Collaboration XXX 2014);
– Tthermal SunyaevZeldovich emission (based on the Planck measurement reported in Planck Collaboration XXI 2014);
– Kinetic SunyaevZeldovich emission, including homogeneous and patchy reionisation components from Shaw et al. (2012) and Battaglia et al. (2013);
– A tSZCIB correlation consistent with both models above; and
– Unresolved point sources as a Poissonlike power spectrum with two components (extragalactic radio galaxies and infrared dusty galaxies).
On top of the cosmological parameters associated with the computation of the CMB spectrum, with HiLLiPoP we sample seven foreground amplitudes (one per emission source, the spectral energy density rescaling the amplitude for each crossfrequency being fixed) and six nuisance parameters (one overall calibration factor plus intercalibrations for each map). See Appendix B for more details.
3.3.3. Largescale polarised likelihood
We construct a polarised EB likelihood based on power spectra, focusing on the large scales where the tensor signal is dominant. Because it carries very little information about the tensor modes, we do not include the TE spectrum in this analysis.
In polarisation, especially at large angular scales, foregrounds are stronger relative to the CMB than in temperature, and cleaning the Planck frequencies using C_{ℓ} templates in the likelihood (as done in temperature) is not accurate enough. In order to clean sky maps of polarised foregrounds, we use the Commander componentseparation code (Eriksen et al. 2008), with a model that includes three polarised components, namely the CMB, synchrotron, and thermal dust emission. Commander was run on each detset map independently, as well as on each realisation from the PR4 Monte Carlo simulations. Maps are available on the PLA in HEALPix^{4} format (Górski et al. 2005) at a resolution N_{side} = 2048.
To compute unbiased estimates of the angular power spectra, we calculate the crosscorrelation of the two detset maps. We make use of two different angular crosspower spectra estimators (described below), which are then concatenated to produce a fullmultipolerange power spectrum. There is no information loss in this process, since the covariances are deduced using Monte Carlo simulations including the correlations over the entire multipole range.
– For multipoles 2 ≤ ℓ ≤ 35, we compute power spectra using an extension of the quadratic maximum likelihood estimator (Tegmark & de OliveiraCosta 2001) adapted for crossspectra in Vanneste et al. (2018)^{5}. At multipoles below 40, it has been shown to produce unbiased polarised power spectra with almost optimal errors. We use downgraded N_{side} = 16 maps after convolution with a cosine apodizing kernel . The signal is then corrected with the PR4 transfer function, to compensate for the filtering induced by the degeneracies between the signal and the templates for systematics in the mapmaking procedure (see Sect. 3.1).
– For multipoles 35 < ℓ < 300, we compute power spectra with a classical pseudoC_{ℓ} estimator Xpol (Sect. 3.3.2). We used N_{side} = 1024 maps and the native beam of Commander maps (i.e. 5′). In this case, we apodize the mask (see Sect. 3.2) with a 1° Gaussian taper. Given the low signaltonoise ratio in polarisation, we bin the spectra with Δℓ = 10.
The EE, BB, and EB power spectra estimates are presented in Fig. 4 for 50% of the sky, which provides the best combination of sensitivity and freedom from foreground residuals. Power spectra computed on different sky fractions (using masks from Sect. 3.2) are compared in Fig. C.1. A simple χ^{2} test on the first 34 multipoles shows no significant departure from the Planck 2018 ΛCDM model for any of these spectra. The ‘probability to exceed’ values (PTE) for the EE, BB, and EB spectra on the first 34 multipoles are 0.27, 0.21, and 0.26, respectively. The most extreme multipole in the BB spectrum is ℓ = 6, which, for a Gaussian distribution, would correspond conditionally to a 3.4σ outlier (reducing to 2.3σ after taking into account the lookelsewhere effect, including the first 34 multipoles). However, at such low multipoles, the distribution is not Gaussian and the PTE are certainly higher than the numbers of σ would suggest. In EE, the largest deviation from the model is for ℓ = 17 at 3.1σ and in EB it is ℓ = 19 at 2.7σ.
Fig. 4. EE, BB, and EB power spectra of the CMB computed on 50% of the sky with the PR4 maps at low (left panels) and intermediate multipoles (right panels). The Planck 2018 ΛCDM model is plotted in black. Grey bands represent the associated cosmic variance. Error bars are deduced from the PR4 Monte Carlo simulations. Correlations between data points are given in Appendix D. A simple χ^{2} test shows no significant departure from the model for any of these spectra. 
The C_{ℓ} covariance matrix is computed from the PR4 Monte Carlos. For each simulation, we compute the power spectra using both estimators. The statistical distribution of the recovered C_{ℓ} then naturally includes the effect of the components included in the Monte Carlo, namely the CMB signal, instrumental noise, Planck systematic effects incorporated in the PR4 simulations (see Sect. 3.1), componentseparation uncertainties, and foreground residuals. The residual power spectra (both for the simulations and the data) are shown in Fig. C.2.
Given the Planck noise level in polarisation, we focus on multipoles below ℓ = 150, which contain essentially all the information on tensor modes in the Planck CMB angular power spectra. At those scales, and given Planck noise levels, the likelihood function needs to consistently take into account the two polarisation fields E and B, as well as all correlations between multipoles and modes (EE, BB, and EB).
LoLLiPoP (LOwℓ LIkelihood on POlarised Powerspectra) is a Planck lowℓ polarisation likelihood based on crossspectra, and was previously applied to Planck EE data for investigating the reionisation history in Planck Collaboration Int. XLVII (2016). The version used here is updated to use crossspectra calculated on componentseparated CMB detset maps processed by Commander from the PR4 frequency maps. Systematic effects are considerably reduced in crosscorrelation compared to autocorrelation, and LoLLiPoP is based on crosspower spectra for which the bias is zero when the noise is uncorrelated between maps. It uses the approximation presented in Hamimeche & Lewis (2008), modified as described in Mangilli et al. (2015) to apply to crosspower spectra. The idea is to apply a change of variable C_{ℓ} → X_{ℓ} so that the new variable X_{ℓ} is nearly Gaussiandistributed. Similarly to Hamimeche & Lewis (2008), we define
where , are the measured crosspower spectra, C_{ℓ} are the power spectra of the model to be evaluated, is a fiducial model, and O_{ℓ} are the offsets needed in the case of crossspectra. For multidimensional CMB modes (here we restrict ourselves to E and B fields only), the C_{ℓ} generalise to C_{ℓ}, a 2 × 2 matrix of power spectra,
and the g function is applied to the eigenvalues of (with C^{−1/2} the square root of the positivedefinite matrix C). In the case of autospectra, the offsets O_{ℓ} are given by the noise bias effectively present in the measured power spectra. For crosspower spectra, the noise bias is zero, and we use effective offsets defined from the C_{ℓ} noise variance:
The distribution of the new variable X_{ℓ} ≡ vecp(X_{ℓ}), the vector of distinct elements of X_{ℓ}, can be approximated as Gaussian, with a covariance given by the covariance of the C_{ℓ}s. The likelihood function of the C_{ℓ} given the data is then
Uncertainties are incorporated into the C_{ℓ}covariance matrix M_{ℓℓ′}, which is evaluated after applying the same pipeline (including Commander component separation and crossspectrum estimation on each simulation) to the Monte Carlo simulations provided in PR4. While foreground emission and the cleaning procedure are kept fixed in the simulations (so that we cannot include uncertainties arising from an imperfect foreground model), the resulting C_{ℓ} covariance consistently includes CMB sample variance, statistical noise, and systematic residuals, as well as foregroundcleaning uncertainties, together with the correlations induced by masking. These uncertainties are then propagated through the likelihood up to the level of cosmological parameters. Figures of the correlation matrices are given in Appendix D.
Using this approach, we are able to derive three different likelihoods, one using only information from E modes (lowlE), one using only information from B modes (lowlB), and one using EE+BB+EB spectra (lowlEB). We have used these likelihoods from ℓ = 2 up to ℓ = 300 with a nominal range of ℓ = [2, 150], since multipoles above ℓ ≃ 150 do not contribute to the result due to the Planck noise (see Sect. 5).
The approach used in this paper is different from the one used for the Planck 2018 results. Indeed, in Planck Collaboration V (2020), the probability density of the polarised spectra at low multipoles was modelled with a polynomial function adjusted on simulations in which only τ is varied, with all other cosmological parameters in a ΛCDM model fixed to the Planck 2018 bestfit values. As a consequence, the probability density is not proportional to the likelihood when the model is not ΛCDM (and in particular for our case ΛCDM+r), and even in the ΛCDM case it neglects correlations with other parameters that affect the posterior on τ. In addition, the simulations used in Planck Collaboration V (2020) were generated with the same CMB realisation for the mapmaking solution. Cosmic variance was included afterwards by adding CMB realisations on top of noiseonly maps, neglecting correlations between foregrounds or systematic templates and the CMB. The information in polarisation at lowℓ was then extracted using a polynomial function fitted to the distribution from simulations. While this is supposed to empirically take into account the effects of systematics on the likelihood shape, it does not include ℓbyℓ correlations, and is limited in the C_{ℓ} power that one can test (for example imposing a strong prior on the EE power at ℓ = 3). As a consequence, the combination of those two effects reduces the covariance, especially at low multipoles, leading to error bars (especially on τ) that are underestimated.
4. Constraints from TT
To derive constraints on the tensortoscalar ratio from the temperature power spectrum, we use the highℓ HiLLiPoP likelihood for 30 ≤ ℓ ≤ 2500, and the Commander likelihood (lowT) in temperature for ℓ < 30, with a prior on the reionisation optical depth to break the degeneracy with the scalar amplitude A_{s}. We use a Gaussian prior τ = 0.055 ± 0.009. For the baseΛCDM model, using PR4 data, we obtain the same results as presented in Planck Collaboration VI (2020).
We now describe the results obtained when fitting the tensortoscalar ratio r in addition to the six ΛCDM parameters (Ω_{b}h^{2}, Ω_{c}h^{2}, θ_{*}, A_{s}, n_{s}, τ). In Planck Collaboration X (2020), the constraint from TT is reported as r_{0.002} < 0.10 (95% CL) using PR3 data. This is much lower than the expected 2σ upper bound on r. Indeed, when we calculate r_{eff} as proposed in Sect. 2, we find that the maximum of the posterior is in the negative region by about 1.7σ. That the maximum happens to fall at negative values is the major reason for the apparently strong constraint on r.
With PR4 data, after marginalising over the other cosmological parameters and the nuisance parameters, we find that the maximum of the posterior is negative by less than 1.2σ when using HiLLiPoP in temperature (hlpTT) along with lowT. As discussed in Planck Collaboration X (2020), this result is related to the lowℓ deficit in the temperature power spectrum. Indeed, removing lowT from the likelihood moves the maximum of the posterior closer to zero, as illustrated in Fig. 5. The corresponding posterior maximum and 68% confidence interval are
Fig. 5. Constraints on the tensortoscalar ratio r_{0.05} based on highℓ temperature data from Planck PR4 (hlpTT), in combination with lowT, and with a prior on τ. 
Using the temperature power spectrum from PR4, we recover the same constraints on other parameters, in particular the scalar spectral tilt n_{s}, as found using PR3 data (see Appendix E). With the full posterior distribution on r, we are able to accurately derive the maximum probability and the uncertainty σ_{r}. The width of the posterior is consistent with the PR3 results. Using only highℓ data, with a prior on the reionisation optical depth τ, we find σ_{r} = 0.12 for TT (consistent with the cosmic variance limit). Note that we find σ_{r} = 0.43 for TE, indicating that TE is much less constraining for r than TT. When adding information from low multipoles in temperature, σ_{r} reduces to 0.094, but at the price of pushing the maximum distribution towards negative values. The posterior maximum is slightly shifted towards zero thanks to the small differences in HiLLiPoP compared to the public Planck likelihood (see Appendix B). The fact that the distribution peaks in the nonphysical domain can be considered as a statistical fluctuation (with a significance between 1 and 2σ, depending on the data set used), which on its own is not a serious problem. However, the fact that this behaviour is strongly related to the deficit of power at lowℓ in temperature is worth noting.
After integrating the positive part of the rposterior, the final upper limits from the Planck temperature power spectrum using PR4 are
5. Constraints from BB
To derive constraints on the tensortoscalar ratio from BB using the PR4 maps, we sample the likelihood with a fixed ΛCDM model based on the Planck 2018 best fit, to which we add tensor fluctuations with a free amplitude parametrised by the tensortoscalar ratio r. We use the LoLLiPoP likelihood described in Sect. 3.3.3, restricted to BB only (referred as ‘lowlB’). As discussed in Sect. 3.3.3, we construct the C_{ℓ} covariance matrix using the PR4 Monte Carlo simulations, which include CMB signal, foreground emission, realistic noise, and systematic effects.
Before giving the final constraints coming from the Planck BB spectra, we should distinguish between the two different regimes, corresponding to large scales (the reionisation bump) and intermediate scales (the recombination bump). Across the reionisation bump, uncertainties are dominated by systematic residuals, as discussed in Sect. 3.1, while foreground residuals may bias the results. Across the recombination bump, uncertainties are dominated by statistical noise; however, systematic effects, as well as foreground residuals, can still bias constraints on r. In order to test the effects of potential foreground residuals, we calculate the posterior distributions of r using various Galactic masks, as described in Sect. 3.2. While large sky fractions (f_{sky} > 60%) show deviations from r = 0, the posteriors for 40, 50, and 60% of the sky are consistent with zero (Fig. F.1). As a robustness test, we also calculate the posterior distribution when changing the range of multipoles (Fig. F.2) and find consistent results, with posteriors compatible with r = 0. Multipoles above ℓ ≃ 150 do not contribute to the result, since the noise in BB is too high. For the rest of this paper, unless otherwise noted, we use a sky fraction of 50%, and compute the likelihood over the range of multipoles from ℓ = 2 to ℓ = 150.
For the reionisation and recombination bumps we find
Both results are obtained over 50% of the sky, with multipoles in the range ℓ = [2, 35] for the former and ℓ = [50, 150] for the latter. With these ranges of multipoles, and given the statistics of the PR4 maps, we can see that the reionisation bump (σ_{r} = 0.110) and the recombination bump (σ_{r} = 0.113) contribute equally to the overall Planck sensitivity to the tensortoscalar ratio.
We can combine the results from the two bumps in order to give the overall constraints on the tensortoscalar ratio from the PlanckBB spectrum (Fig. 6). The full constraint on r from the PR4 BB spectrum over 50% of the sky, including correlations between all multipoles between ℓ = 2 and ℓ = 150, is
Fig. 6. Posterior distribution of r from PR4 data, using LoLLiPoP and the BB spectrum on 50% of the sky (black). Constraints from the reionisation bump and the recombination bump are plotted in red and blue, respectively. Constraints from Planck BB with the full multipole range ℓ = [2, 150] are in black. 
This is fully compatible with no tensor signal, and we can derive an upper limit by integrating the posterior distribution out to 95%, after applying the physical prior r > 0, which yields
This result can be compared with the BICEP2/Keck Array constraints (BICEP2 Collaboration 2018) of
with σ_{r} = 0.02 compared to σ_{r} = 0.069 for the Planck result presented in this analysis
6. Additional constraints from polarisation
As shown in Fig. 1, the EE tensor spectrum is similar in amplitude to the BB tensor spectrum, even though the scalar mode in EE is stronger. Given that noise dominates the tensor signal at all multipoles in both EE and BB, we expect the likelihood for EE to give useful constraints on r. We thus present the constraints from polarised lowℓ data (ℓ < 150) using different combinations of the LoLLiPoP likelihood (specifically EE, BB, and EE+BB+EB) in Fig. 7. We emphasise that EE+BB+EB is a likelihood of the correlated polarisation fields E and B and not the combination of individual likelihoods (see Sect. 3.3.3).
Fig. 7. Posterior distributions for r from Planck polarised lowℓ data (ℓ < 150) using LoLLiPoP and the EE, BB, and EE+BB+EB spectra. The dashed black line is obtained from EB data by fitting a BB tensor model. The sky fraction used here is f_{sky} = 50%. 
The first thing to notice is that the posterior distribution for EE peaks at r = 0.098 ± 0.097, while the other modes give results compatible with zero within 1σ. Given the lower sensitivity of lowlE to r (σ_{r} ≃ 0.10) compared to that of lowlB (σ_{r} ≃ 0.07), this is mitigated when adding the information from other modes. The posterior distributions for r give
As a consistency check, Fig. 7 also shows the constraints when fitting the BB tensor model on the EB data power spectrum, which is compatible with zero (r = −0.012 ± 0.068) as expected.
Using polarisation data, Planck’s sensitivity to the tensortoscalar ratio reaches σ_{r} = 0.046. Combining all Planck polarisation modes (EE, BB, and EB) out to ℓ = 150 leads to the following upper limit:
Note that this constraint is almost independent of the other ΛCDM parameters, and in particular the reionisation optical depth τ. To demonstrate this, using the same data set (lowlB and lowlEB), we derive 2dimensional constraints for τ and r and plot them in Fig. 8. The constraint is stable when sampling for τ. Indeed, in this case, we obtain
Fig. 8. LoLLiPoP posterior distribution in the τ–r plane using lowlE (blue), lowlB (red), and lowlEB (black). The sky fraction here is f_{sky} = 50%. 
and for the reionisation optical depth
compatible with lowlE results, while lowlB shows no detection of τ, since BB is dominated by noise.
7. Combined results
Up to this point, the constraints on r have been derived relative to a fixed fiducial ΛCDM spectrum based on the Planck 2018 results. Including the Planck temperature likelihoods (both lowT and hlpTT) in a combined analysis of the Planck CMB spectra allows us to properly propagate uncertainties from other cosmological parameters to r, as well as to selfconsistently derive constraints in the n_{s}–r plane. In this section, we combine the lowT and hlpTT with the lowℓ polarised likelihood lowlEB to sample the parameter space of the ΛCDM+r model. The comparison of contours at 68% and 95% confidence levels between PR3 and PR4 data is presented in Fig. G.1.
We also include the BK15 constraints from BICEP2 Collaboration (2018). When combining Planck and BK15, we neglect the correlation between the two data sets and simply multiply the likelihood distributions. This is justified because the BK15 spectra are estimated on 1% of the sky, while the Planck analysis is derived from 50% of the sky.
Figure 9 gives posteriors on r after marginalisation over the nuisance and the other ΛCDM cosmological parameters. We obtain the following 95% CL upper limits:
Fig. 9. Posterior distributions for r after marginalisation over the nuisance parameters and the other ΛCDM parameters, for the Planck temperature data (hlpTT+lowT) in combination with BK15 and the largescale polarised Planck likelihood (lowlEB). 
Figure 10 shows the constraints in the r–n_{s} plane for Planck data in combination with BK15. The constraints from the full combination of Planck data are comparable to those from BK15. The addition of the highℓ TE likelihood produces tighter constraints on the spectral index n_{s} (as already reported in Planck Collaboration VI 2020).
Fig. 10. Marginalised joint 68% and 95% CL regions for n_{s} and r_{0.05} from Planck alone (hlp+lowT+lowlEB) and in combination with BK15. The solid lines correspond to hlpTT+lowT+lowlEB, while the filled regions include TE and correspond to hlpTTTE+lowT+lowlEB. 
There have been several other attempts to constrain the value of r, particularly through measurements of the BB power spectrum. As we have already stressed, there is a weak limit from the TT spectrum and at the current sensitivity level for r, the constraints from EE are about as powerful as those from BB; hence the tensor constraints in this paper are derived from a combination of BB limits with those coming from TT and EE. We show in Appendix H a comparison of our BB limits with those of other experiments.
8. Conclusions
In this paper, we have derived constraints on the amplitude of tensor perturbations using Planck PR4 data. We investigated the intrinsic sensitivity of the TT spectrum, which is cosmicvariance limited, and found σ_{r} = 0.094 using the full range of multipoles. We noted the impact of the lowℓ anomaly, which pushes the maximum posterior distribution towards negative values of r_{eff} at roughly the 1σ level.
For the first time, we analysed the PlanckBB spectrum for r and obtained σ_{r} = 0.069, which is lower than in temperature. The PlanckBmode spectrum, being dominated by noise, gives a constraint on r that is fully compatible with zero from both low and intermediate multipoles, in other words from both the reionisation and recombination peaks. Multipoles above ℓ ≃ 150 do not contribute to the result, since the noise in BB is too high.
Using an appropriate likelihood in polarisation, we showed that the PlanckEE spectrum is also sensitive to the amplitude of the tensortoscalar ratio r. The combined constraints from PlanckEE and BB, including EB correlations, lead to a sensitivity on r of σ_{r} = 0.046, two times better than in temperature. We also investigated the impact of foreground residuals using different Galactic cuts and by varying the range of multipoles used in the polarised likelihood. Finally, by combining temperature and polarisation constraints, we derived the posterior distribution on r marginalised over the ΛCDM cosmological parameters and nuisance parameters, including uncertainties from systematics (both instrumental and astrophysical). The result gives an upper limit of r < 0.056 at the 95% confidence level using Planck data only. In combination with the BICEP/Keck measurements from 2015, this constraint is further reduced to r < 0.044 (95% CL), the tightest limit on r to date.
Acknowledgments
Planck is a project of the European Space Agency (ESA) with instruments provided by two scientific consortia funded by ESA member states and led by Principal Investigators from France and Italy, telescope reflectors provided through a collaboration between ESA and a scientific consortium led and funded by Denmark, and additional contributions from NASA (USA). Some of the results in this paper have been derived using the HEALPix package. This research used resources of the National Energy Research Scientific Computing Center (NERSC), a U.S. Department of Energy Office of Science User Facility operated under Contract No. DEAC0205CH11231. We gratefully acknowledge support from the CNRS/IN2P3 Computing Center for providing computing and dataprocessing resources needed for this work.
References
 Albrecht, A., & Steinhardt, P. J. 1982, Phys. Rev. Lett., 48, 1220 [NASA ADS] [CrossRef] [Google Scholar]
 Bardeen, J. M., Steinhardt, P. J., & Turner, M. S. 1983, Phys. Rev. D, 28, 679 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Battaglia, N., Natarajan, A., Trac, H., Cen, R., & Loeb, A. 2013, ApJ, 776, 83 [Google Scholar]
 Bennett, C. L., Larson, D., Weiland, J. L., et al. 2013, ApJS, 208, 20 [Google Scholar]
 Béthermin, M., Daddi, E., Magdis, G., et al. 2012, ApJ, 757, L23 [NASA ADS] [CrossRef] [Google Scholar]
 BICEP2 Collaboration 2018, Phys. Rev. Lett., 121, 221301 [CrossRef] [PubMed] [Google Scholar]
 Blas, D., Lesgourgues, J., & Tram, T. 2011, JCAP, 2011, 034 [NASA ADS] [CrossRef] [Google Scholar]
 Brout, R., Englert, F., & Gunzig, E. 1978, Ann. Phys., 115, 78 [NASA ADS] [CrossRef] [Google Scholar]
 Choi, S. K., Hasselfield, M., Ho, S. P. P., et al. 2020, JCAP, 12, 045 [Google Scholar]
 Couchot, F., HenrotVersillé, S., Perdereau, O., et al. 2017a, A&A, 602, A41 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Couchot, F., HenrotVersillé, S., Perdereau, O., et al. 2017b, A&A, 597, A126 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Eriksen, H. K., Jewell, J. B., Dickinson, C., et al. 2008, ApJ, 676, 10 [NASA ADS] [CrossRef] [Google Scholar]
 Górski, K. M., Hivon, E., Banday, A. J., et al. 2005, ApJ, 622, 759 [NASA ADS] [CrossRef] [Google Scholar]
 Guth, A. H. 1981, Phys. Rev. D, 23, 347 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Guth, A. H., & Pi, S. 1982, Phys. Rev. Lett., 49, 1110 [NASA ADS] [CrossRef] [Google Scholar]
 Hamimeche, S., & Lewis, A. 2008, Phys. Rev. D, 77, 103013 [Google Scholar]
 Hawking, S. 1982, Phys. Lett. B, 115, 295 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Hivon, E., Mottet, S., & Ponthieu, N. 2017, A&A, 598, A25 [CrossRef] [EDP Sciences] [Google Scholar]
 Kamionkowski, M., Kosowsky, A., & Stebbins, A. 1997, Phys. Rev. Lett., 78, 2058 [Google Scholar]
 Kazanas, D. 1980, ApJ, 241, L59 [NASA ADS] [CrossRef] [Google Scholar]
 Linde, A. D. 1982, Phys. Lett. B, 108, 389 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Linde, A. D. 1983, Phys. Lett. B, 129, 177 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Mangilli, A., Plaszczynski, S., & Tristram, M. 2015, MNRAS, 453, 3174 [NASA ADS] [CrossRef] [Google Scholar]
 Mukhanov, V. F. 1985, JETP Lett., 41, 493 [NASA ADS] [Google Scholar]
 Mukhanov, V. F., & Chibisov, G. 1981, JETP Lett., 33, 532 [Google Scholar]
 Mukhanov, V. F., & Chibisov, G. 1982, Sov. Phys. JETP, 56, 258 [Google Scholar]
 Planck Collaboration XV. 2014, A&A, 571, A15 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XVI. 2014, A&A, 571, A16 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XXI. 2014, A&A, 571, A21 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XXII. 2014, A&A, 571, A22 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XXX. 2014, A&A, 571, A30 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XI. 2016, A&A, 594, A11 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XX. 2016, A&A, 594, A20 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XXII. 2016, A&A, 594, A22 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration III. 2020, A&A, 641, A3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration V. 2020, A&A, 641, A5 [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration VI. 2020, A&A, 641, A6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration VIII. 2020, A&A, 641, A8 [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration X. 2020, A&A, 641, A10 [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration Int. XXX. 2016, A&A, 586, A133 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration Int. XLI. 2016. A&A, 596, A102 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration Int. XLVII. 2016, A&A, 596, A108 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration Int. LI. 2017, A&A, 607, A95 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration Int. LVII. 2020, A&A, 643, A42 [CrossRef] [EDP Sciences] [Google Scholar]
 POLARBEAR Collaboration (Ade, P. A. R., et al.) 2017, ApJ, 848, 121 [NASA ADS] [CrossRef] [Google Scholar]
 POLARBEAR Collaboration (Adachi, S., et al.) 2020, ApJ, 897, 55 [Google Scholar]
 Sato, K. 1981, MNRAS, 195, 467 [NASA ADS] [CrossRef] [Google Scholar]
 Sayre, J. T., Reichardt, C. L., Henning, J. W., et al. 2020, Phys. Rev. D, 101, 122003 [Google Scholar]
 Seljak, U. 1997, ApJ, 482, 6 [Google Scholar]
 Seljak, U., & Zaldarriaga, M. 1997, Phys. Rev. Lett., 78, 2054 [Google Scholar]
 Shaw, L. D., Rudd, D. H., & Nagai, D. 2012, ApJ, 756, 15 [Google Scholar]
 Starobinsky, A. A. 1980, Phys. Lett. B, 91, 99 [NASA ADS] [CrossRef] [Google Scholar]
 Starobinsky, A. A. 1982, Phys. Lett. B, 117, 175 [NASA ADS] [CrossRef] [Google Scholar]
 Tegmark, M., & de OliveiraCosta, A. 2001, Phys. Rev. D, 64, 063001 [NASA ADS] [CrossRef] [Google Scholar]
 Tristram, M., MacíasPérez, J. F., Renault, C., & Santos, D. 2005, MNRAS, 358, 833 [NASA ADS] [CrossRef] [Google Scholar]
 Tucci, M., Toffolatti, L., de Zotti, G., & MartínezGonzález, E. 2011, A&A, 533, A57 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Vanneste, S., HenrotVersillé, S., Louis, T., & Tristram, M. 2018, Phys. Rev. D, 98, 103526 [NASA ADS] [CrossRef] [Google Scholar]
 Zaldarriaga, M., & Seljak, U. 1998, Phys. Rev., D, 58 [Google Scholar]
Appendix A: The NPIPEBB transfer function
The total signal simulations in the PR4 release have insufficient S/N to determine the CMB BB transfer function (Planck Collaboration Int. LVII 2020) directly. We therefore study the use of the wellmeasured CMB EE transfer function in place of the unknown BB function. Such an approach is obviously approximate, but should be sufficient for an analysis providing an upper limit to the BB amplitude, provided that the suppression of the EE power in NPIPE is at least as strong as the BB suppression.
NPIPE filtering occurs during the destriping process, when timedomain templates are fitted against and subtracted from the TOD. The filtering happens because the signal model in the calibration step does not include CMB polarisation. The processing approach mimics simple linear regression, except that it is performed in a subspace that does not include skysynchronous degrees of freedom.
Here we set up a simplified test where we generate one survey’s worth of singledetector, timeordered data corresponding to pure CMB E and pure CMB B modes, and regress out all of the destriping templates, as visualised in Appendix E in Planck Collaboration Int. LVII (2020). We then project both the unfiltered and filtered CMB signals into maps, and compute the ‘filtering function’ for this test, which measures the amount of signal suppression between them as
For the purpose of this test, singledetector, singlesurvey data can be projected onto an intensity map, since the vast majority of the sky pixels are only observed in one orientation. In Eq. (A.1), is the pseudoautospectrum of a map binned from the unfiltered signal outside a ±15° cut in Galactic latitude, while is the pseudocrossspectrum between the filtered and unfiltered maps. This test produces a CMB EE filtering function that resembles the magnitude and angular extent of the transfer function measured from the total signal simulations. This indicates that the missing elements of the data processing (full sky coverage, multiple surveys, multiple detectors, and projection operator) do not prevent us from drawing conclusions about the similarities between the E and Bmode transfer functions.
Results of 300 Monte Carlo realisations are shown in Fig. A.1. They demonstrate good agreement between the EE and BB filtering functions and broad compatibility between the full simulation transfer function and our approximate test. There is a statistically significant dip in the EE/BB filtering function ratio in the range ℓ ≃ 10–20, suggesting that correcting the BB spectra with the EE transfer function might result in an overcorrection by a few percent, marginally inflating the upper limit. Given the approximate nature of this test and the modest size of the correction, we have not sought to tighten our upper limits by this difference.
Fig. A.1. Linear regression test to compare coupling between CMB E and B modes and the NPIPE destriping templates. Top: solid lines show the median bandpowers in individual multipole bins. The filled bands represent the 68% confidence region around the median, as measured over 300 realisations. Bottom left: we compare measured EE and BB filtering functions derived from our simplified model with the EE transfer function taken from PR4. Bottom right: ratio of the EE and BB filtering functions. 
Appendix B: The hillipop likelihood
HiLLiPoP (Highℓ Likelihood on Polarised Powerspectra) is one of the highℓ likelihoods developed for analysis of the Planck data. It was used as part of the 2013 (Planck Collaboration XV 2014) and 2015 (Planck Collaboration XI 2016) releases and also described in Couchot et al. (2017b) and Couchot et al. (2017a).
HiLLiPoP is very similar to the Planck public likelihood (plik). HiLLiPoP is a Gaussian likelihood based on crossspectra from the HFI 100, 143, and 217GHz detset splitmaps. The crossspectra are estimated using a pseudoC_{ℓ} algorithm with a mask adapted to each frequency to reduce the contamination from Galactic emission and point sources. The C_{ℓ} model includes foreground residuals on top of the CMB signal. These foreground residuals are both Galactic (dust emission) and extragalactic (CIB, tSZ, kSZ, SZ×CMB, and point sources). HiLLiPoP also introduced nuisance parameters to take into account map calibration uncertainties.
The most significant differences compared with plik are that HiLLiPoP uses:
– detsetsplit maps instead of time splits (so that the crossspectra do not need to account for a noisecorrelation correction as in plik);
– pointsource masks that were obtained from a procedure that extracts compact Galactic structures;
– a Galactic dust C_{ℓ} model that (as a result of the mask difference) follows closely and is parametrised by the power law discussed in Planck Collaboration Int. XXX (2016), while plik uses an adhoc effective function in ℓ;
– foreground templates derived from Planck measurements (Planck Collaboration Int. XXX 2016 for the dust emission, Planck Collaboration XXX 2014 for the CIB, and Planck Collaboration XXII 2016 for the SZ) with a free amplitude for each emission mechanism, but spectral energy distributions fixed by Planck measurements;
– a twocomponent model for the signal from unresolved point sources, which incorporates the contribution from extragalactic radio (Tucci et al. 2011) and infrared dusty (Béthermin et al. 2012) galaxies, as well as taking into account the variation of the flux cut across the sky and in ‘incompleteness’ of the source catalogue at each frequency;
– all the 15 crossspectra built from the 100, 143, and 217GHz detset maps (while plik keeps only five of them);
– all multipole values (while plik bins the power spectra).
In the end, we have a total of six instrumental parameters (for map calibration) and nine astrophysical parameters (seven for TT, one for TE, and one for EE) in addition to the cosmological parameters (see Table B.1).
Nuisance parameters for the HiLLiPoP likelihood.
Using HiLLiPoP on the PR3 data, we recover essentially identical constraints on the 6parameter baseΛCDM model as the public Planck likelihood. The results for the ΛCDM+r model with PR4 are discussed in Appendix E.
Appendix C: Largescale polarised angular power spectra
We estimate CMB largescale polarised power spectra by crosscorrelating the two independent detset splits, A and B, after Commander component separation. As detailed in Sect. 3.3.3, we use two different powerspectrum estimators for the low (2 ≤ ℓ ≤ 35) and the intermediate (35 < ℓ < 300) multipole ranges. For the lower multipole range we use a quasiQML estimator (Vanneste et al. 2018), while for the higher multipoles we use a classic pseudoC_{ℓ} approach (a generalisation to polarisation of the method presented in Tristram et al. 2005) with a binning of Δℓ = 10.
Figure C.1 shows the reconstructed EE, BB, and EB power spectra for various sky fractions, from 30 to 70%. The 6parameter ΛCDM model based on the best fit to the Planck 2018 data is plotted in black, together with the cosmic variance (computed for the full sky) in grey. The spectra show a remarkable consistency with the model, and are largely insensitive to sky fraction. The BB and EB spectra are dominated by noise. At low multipoles, the BB spectrum computed on the largest sky fraction (70%) exhibits an excess at very low multipoles (ℓ ≤ 5), attributable to Galactic residuals. For both BB and EB, a reduction of the sky fraction corresponds to a larger dispersion.
Fig. C.1. EE, BB, and EB power spectra computed from the PR4 maps for sky fractions from 30 to 70%. The black lines represent the ΛCDM model and the grey bands show its associated fullsky cosmic variance. 
Figure C.2 shows the residuals for the EE, BB, and EB power spectra compared to the bestfit baseΛCDM model from the Planck 2018 results. The plot on the left shows the residuals from the Monte Carlo simulations computed as the average spectra over the simulations divided by the uncertainty in the mean, . The average here is computed independently for each multipole. The plot on the right shows the residuals of the PR4 data compared to the ΛCDM model divided by the spectrum uncertainties (i.e. the dispersion over the simulations).
Fig. C.2. Residuals of the EE, BB, and EB power spectra. Left: mean value of the Monte Carlo simulations compared to the error on the mean, . Right: PR4 data compared to the standard deviation σ_{ℓ}. Grey bands show the 1, 2, and 3σ levels. 
Appendix D: Crossspectrum correlation matrix
Uncertainties are propagated to the likelihood function through the C_{ℓ} covariance matrices (see Sect. 3.3.3). We use Monte Carlo simulations to estimate the covariance over the multipoles considered in this analysis (2 ≤ ℓ ≤ 150). Four hundred simulations are available, allowing us to invert the estimated covariance matrices of rank 45 for BB and 135 for the full EE+BB+EB. While the uncertainties would be better determined with a larger sample of simulations, we have seen no impact on the results when estimating the covariance matrix on a reduced number of simulations, and we conclude that our sample of 400 is enough to propagate the uncertainties of the C_{ℓ}s to r, as well as to τ.
Figure D.1 shows the correlations of the computed from the covariance matrix. Correlations from bin to bin are below 15%, except for the nexttoneighbour bins at ℓ < 40. Figure D.2 shows the correlations computed from the full covariance matrix, including EE, BB, and EB spectra for all multipole bins from 2 to 150.
Fig. D.1. Left: correlation matrix for . Right: enlargement of the upperleft corner. Covariances are estimated from 400 endtoend simulations, including signal, noise, systematics, and foreground residuals. 
Fig. D.2. Correlation matrix for the , , and . Covariances are estimated from 400 endtoend simulations, including signal, noise, systematics, and foreground residuals. 
Appendix E: ΛCDM+r parameters for PR3 and PR4
In this paper, we use a new data set, PR4, and an alternative likelihood, HiLLiPoP. Figure E.1 shows that we obtain essentially the same parameter values for the ΛCDM+r model based on the temperature power spectrum as are obtained from PR3 and plik (available on the Planck Legacy Archive). The differences between PR3 and PR4 are primarily seen in polarisation.
Fig. E.1. Constraint contours (at 68 and 95% confidence) on parameters of a ΛCDM+r model using PR4 and HiLLiPoP (red), compared to those obtained with PR3 and plik (blue). Both sets of results use the Commander likelihood for temperature multipoles ℓ ≤ 30, and a Gaussian prior to constrain τ. 
Appendix F: Robustness tests
As a test of robustness, we computed the posterior distributions for r_{eff} considering spectra estimated on different sky fractions from 30 to 70% (Sect. 3.2). Figure F.1 shows results for the BB (left) and EE+BB+EB (right) likelihoods. A sky fraction of 50% provides the best combination of sensitivity (which increases with the sky fraction) and freedom from foreground residuals (which decreases with sky fraction). Figure F.2 shows that the BB posterior is robust with respect to the choice of ℓ_{min} up to ℓ_{min} ≈ 20, and that the width of the posterior is stable for ℓ_{max} ≳ 150, a consequence of the increase of signaltonoise ratio with multipole.
Fig. F.1. Posterior distributions for r using LoLLiPoP and BB (panel a) and EE + BB + EB (panel b) for sky fractions from 30% to 70%. The multipole range is the same in all cases, ℓ = [2, 145]. The best combination of sensitivity and freedom from foreground residuals is achieved at 50%. 
Fig. F.2. Posterior distributions for r using LoLLiPoP and the BB power spectrum for various multipole ranges. The sky fraction is fixed at 50%. 
Appendix G: Triangle plot for ΛCDM+r parameters
Figure G.1 shows the differences in parameter values obtained for PR3 (plikTT+lowT+lowE, Planck Legacy Archive) and PR4 (hlpTT+lowT+lowlEB, this analysis) when both temperature and polarisation data are included. There are no significant differences. In general, uncertainties are slightly smaller with PR4.
Fig. G.1. Contour constraints for ΛCDM+r parameters for PR3 (plikTT+lowT+lowE) and for PR4 (hlpTT+lowT+lowlEB). 
Appendix H: Comparison with other BB measurements
PlanckTT, EE, and BB power spectra all contribute to the constraints on r derived in this paper. Most previous limits on r have been determined only from BB measurements. For the convenience of readers who might be interested in how PlanckBB data alone compare to those of other experiments, Fig. H.1 shows the PlanckBB bandpowers as upper limits at 95% CL, along with other measurements from the BICEP2/Keck Array (BICEP2 Collaboration 2018), SPTpol (Sayre et al. 2020), ACTPol (Choi et al. 2020), POLARBEAR (POLARBEAR Collaboration 2017, 2020), and WMAP (Bennett et al. 2013).
Fig. H.1. Bmode polarisation power spectrum measurements from Planck PR4 (this work), SPTpol (Sayre et al. 2020), ACTPol (Choi et al. 2020), POLARBEAR (POLARBEAR Collaboration 2017, 2020), the BICEP2/Keck Array (BICEP2 Collaboration 2018), and WMAP (Bennett et al. 2013). Uncertainties are 68% confidence levels, while bandpowers compatible with zero at 95% CL are plotted as upper limits. The solid red curve shows the BB lensing spectrum based on the Planck 2018 ΛCDM bestfit model (Planck Collaboration VI 2020). The dashed red curve corresponds to the BB power spectrum for r = 0.044, the upper limit obtained in this work when combining Planck (TT, EE, and BB) with BK15. 
All Tables
All Figures
Fig. 1. Scalar (thick solid lines) versus tensor spectra for r = 0.1 (dashed lines) and r = 0.01 (dotted lines). Spectra for TT are in black, EE in blue, and BB in red. The red solid line corresponds to the signal from BB lensing. 

In the text 
Fig. 2. Variances for crossspectra in EE and BB based on PR4 simulations, including: cosmic (sample) variance (black); analytic statistical noise (red); and PR4 noise from Monte Carlo simulations (green), including noise and systematics with (solid line) or without (dashed line) correction for the transfer function. The sky fraction used here is 80%, as it illustrates well the effect of both systematics and transferfunction corrections (Planck Collaboration Int. LVII 2020). 

In the text 
Fig. 3. Galactic masks used for the Planck likelihoods. The mask shown in dark blue indicates the sky rejected in order to retain a 70% sky fraction for analysis. The masks shown in light blue, green, orange and red incrementally omit further parts of the sky, corresponding in turn to 60, 50, 40 and 30% retained sky fractions, the latter shown in white. 

In the text 
Fig. 4. EE, BB, and EB power spectra of the CMB computed on 50% of the sky with the PR4 maps at low (left panels) and intermediate multipoles (right panels). The Planck 2018 ΛCDM model is plotted in black. Grey bands represent the associated cosmic variance. Error bars are deduced from the PR4 Monte Carlo simulations. Correlations between data points are given in Appendix D. A simple χ^{2} test shows no significant departure from the model for any of these spectra. 

In the text 
Fig. 5. Constraints on the tensortoscalar ratio r_{0.05} based on highℓ temperature data from Planck PR4 (hlpTT), in combination with lowT, and with a prior on τ. 

In the text 
Fig. 6. Posterior distribution of r from PR4 data, using LoLLiPoP and the BB spectrum on 50% of the sky (black). Constraints from the reionisation bump and the recombination bump are plotted in red and blue, respectively. Constraints from Planck BB with the full multipole range ℓ = [2, 150] are in black. 

In the text 
Fig. 7. Posterior distributions for r from Planck polarised lowℓ data (ℓ < 150) using LoLLiPoP and the EE, BB, and EE+BB+EB spectra. The dashed black line is obtained from EB data by fitting a BB tensor model. The sky fraction used here is f_{sky} = 50%. 

In the text 
Fig. 8. LoLLiPoP posterior distribution in the τ–r plane using lowlE (blue), lowlB (red), and lowlEB (black). The sky fraction here is f_{sky} = 50%. 

In the text 
Fig. 9. Posterior distributions for r after marginalisation over the nuisance parameters and the other ΛCDM parameters, for the Planck temperature data (hlpTT+lowT) in combination with BK15 and the largescale polarised Planck likelihood (lowlEB). 

In the text 
Fig. 10. Marginalised joint 68% and 95% CL regions for n_{s} and r_{0.05} from Planck alone (hlp+lowT+lowlEB) and in combination with BK15. The solid lines correspond to hlpTT+lowT+lowlEB, while the filled regions include TE and correspond to hlpTTTE+lowT+lowlEB. 

In the text 
Fig. A.1. Linear regression test to compare coupling between CMB E and B modes and the NPIPE destriping templates. Top: solid lines show the median bandpowers in individual multipole bins. The filled bands represent the 68% confidence region around the median, as measured over 300 realisations. Bottom left: we compare measured EE and BB filtering functions derived from our simplified model with the EE transfer function taken from PR4. Bottom right: ratio of the EE and BB filtering functions. 

In the text 
Fig. C.1. EE, BB, and EB power spectra computed from the PR4 maps for sky fractions from 30 to 70%. The black lines represent the ΛCDM model and the grey bands show its associated fullsky cosmic variance. 

In the text 
Fig. C.2. Residuals of the EE, BB, and EB power spectra. Left: mean value of the Monte Carlo simulations compared to the error on the mean, . Right: PR4 data compared to the standard deviation σ_{ℓ}. Grey bands show the 1, 2, and 3σ levels. 

In the text 
Fig. D.1. Left: correlation matrix for . Right: enlargement of the upperleft corner. Covariances are estimated from 400 endtoend simulations, including signal, noise, systematics, and foreground residuals. 

In the text 
Fig. D.2. Correlation matrix for the , , and . Covariances are estimated from 400 endtoend simulations, including signal, noise, systematics, and foreground residuals. 

In the text 
Fig. E.1. Constraint contours (at 68 and 95% confidence) on parameters of a ΛCDM+r model using PR4 and HiLLiPoP (red), compared to those obtained with PR3 and plik (blue). Both sets of results use the Commander likelihood for temperature multipoles ℓ ≤ 30, and a Gaussian prior to constrain τ. 

In the text 
Fig. F.1. Posterior distributions for r using LoLLiPoP and BB (panel a) and EE + BB + EB (panel b) for sky fractions from 30% to 70%. The multipole range is the same in all cases, ℓ = [2, 145]. The best combination of sensitivity and freedom from foreground residuals is achieved at 50%. 

In the text 
Fig. F.2. Posterior distributions for r using LoLLiPoP and the BB power spectrum for various multipole ranges. The sky fraction is fixed at 50%. 

In the text 
Fig. G.1. Contour constraints for ΛCDM+r parameters for PR3 (plikTT+lowT+lowE) and for PR4 (hlpTT+lowT+lowlEB). 

In the text 
Fig. H.1. Bmode polarisation power spectrum measurements from Planck PR4 (this work), SPTpol (Sayre et al. 2020), ACTPol (Choi et al. 2020), POLARBEAR (POLARBEAR Collaboration 2017, 2020), the BICEP2/Keck Array (BICEP2 Collaboration 2018), and WMAP (Bennett et al. 2013). Uncertainties are 68% confidence levels, while bandpowers compatible with zero at 95% CL are plotted as upper limits. The solid red curve shows the BB lensing spectrum based on the Planck 2018 ΛCDM bestfit model (Planck Collaboration VI 2020). The dashed red curve corresponds to the BB power spectrum for r = 0.044, the upper limit obtained in this work when combining Planck (TT, EE, and BB) with BK15. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.