## Neutrinos and dark energy after Planck and BICEP2: data consistency tests and cosmological parameter constraints

Jing-Fei Zhang, a Jia-Jia Geng, a Xin Zhang a,b, 1

a Department of Physics, College of Sciences, Northeastern University, Shenyang 110004, China

b Center for High Energy Physics, Peking University,

Beijing 100080, China

E-mail: jfzhang@mail.neu.edu.cn, gengjiajia163@163.com,

zhangxin@mail.neu.edu.cn

Abstract. The detection of the B-mode polarization of the cosmic microwave background (CMB) by the BICEP2 experiment implies that the tensor-to-scalar ratio r should be involved in the base standard cosmology. In this paper, we extend the ΛCDM+ +neutrino/dark rar diation models by replacing the cosmological constant with the dynamical dark energy with constant w . Four neutrino plus dark energy models are considered, i.e., the w CDM+ + r ∑ m ν , w CDM+ + r N eff , w CDM+ + r ∑ m ν + N eff , and w CDM+ + r N eff + m eff ν, sterile models. The current observational data considered in this paper include the Planck temperature data, the WMAP 9-year polarization data, the baryon acoustic oscillation data, the Hubble constant direct measurement data, the Planck Sunyaev-Zeldovich cluster counts data, the Planck CMB lensing data, the cosmic shear data, and the BICEP2 polarization data. We test the data consistency in the four cosmological models, and then combine the consistent data sets to perform joint constraints on the models. We focus on the constraints on the parameters w , ∑ m ν , N eff , and m eff ν, sterile .

1 Corresponding author.

## Contents

| 1   | Introduction                                              |   1 |
|-----|-----------------------------------------------------------|-----|
| 2   | Methodology                                               |   3 |
| 3   | Results and discussion                                    |   5 |
| 3.1 | Data consistency tests                                    |   5 |
| 3.2 | Constraints on dark energy, dark radiation, and neutrinos |   8 |
| 3.3 | Early-universe parameters and BICEP2                      |  11 |
| 4   | Conclusion                                                |  12 |
| A   | Detailed constraint results                               |  14 |

## 1 Introduction

Discovery of neutrino oscillations has indicated that neutrinos have nonzero masses. Neutrino oscillation experiments can only place limits on the squared mass differences between the neutrino mass eigenstates: the solar and reactor experiments observed ∆ m 2 21 glyph[similarequal] 8 × 10 -5 eV , 2 and the atmospheric and accelerator beam experiments observed ∆ m 2 32 glyph[similarequal] 3 × 10 -3 eV . To measure the absolute masses of neutrinos, different experiments are needed. In fact, 2 cosmological data have been providing tight limits on the total mass of neutrinos.

The cosmic microwave background (CMB) observations have been used to constrain the neutrino mass and possibly the extra relativistic degrees of freedom (sometimes referred to as 'dark radiation') [1-3]. The latest CMB temperature power spectrum measured by the Planck satellite mission provided the tight limits on the total mass of active neutrinos, ∑ m ν , and the effective number of relativistic species, N eff [4]. For example, the Planck+WP+highL data combination (here, WP denotes the WMAP 9-year polarization data, and highL denotes the ACT and SPT temperature data) gives the 95% confidence level (CL) limits: ∑ m &lt; ν 0 66 eV for the case of no extra relics ( . N eff = 3 046) and . N eff = 3 36 . +0 68 . -0 64 . for the case of minimal-mass normal hierarchy for the neutrino masses (only one massive eigenstate with m ν = 0 06 eV). Late-time geometric measurements can be used to help reduce some geometric . degeneracies and thus improve constraints. Therefore, the baryon acoustic oscillation (BAO) data are very useful in the parameter estimation. Note also that the BAO data are proven to be in good agreement with the Planck data, so one can always combine CMB data with BAO data without any question. The Planck+WP+highL+BAO data combination changes the above limits to: ∑ m &lt; . ν 0 23 eV for the case of no extra relics and N eff = 3 30 . +0 54 . -0 51 . for the case of minimal-mass normal hierarchy model.

The measurements of the growth of large-scale structure also play a crucial role in constraining the neutrino mass. The Planck Collaboration reported a result of counts of rich clusters of galaxies by analyzing its sample of thermal Sunyaev-Zeldovich (SZ) clusters, σ 8 (Ω m / . 0 27) 0 3 . = 0 782 . ± 0 010 [5]. . And the weak lensing data can also give another combination of σ 8 and Ω m . The cosmic shear measurement provided by the CFHTLenS survey gave σ 8 (Ω m / . 0 27) 0 46 . = 0 774 . ± 0 040 [6]. . However, based on the six-parameter base ΛCDM cosmology, the Planck data predict more clusters than these astrophysical measurements

observe. The Planck+WP+highL data combination gives σ 8 (Ω m / . 0 27) 0 3 . = 0 87 . ± 0 02 and . σ 8 (Ω m / . 0 27) 0 46 . = 0 89 . ± 0 03, in tensions with the above two measurements at the 4.0 . σ and 2.3 σ levels, respectively. In fact, the tensions might just imply non-zero neutrino masses (see, e.g., Ref. [7]). Increasing neutrino masses could suppress the growth of structure below their free-streaming length, allowing σ 8 to be substantially lower.

On the other hand, the measurement of the Hubble constant H 0 is also very useful in breaking parameter degeneracies in cosmological constraints, but the direct measurement of H 0 is in tension with the Planck result based on the standard cosmology. The direct measurement gives H 0 = 73 8 . ± 2 4 km s . -1 Mpc -1 [8], but the six-parameter ΛCDM model fitting to Planck+WP+highL data gives H 0 = 67 3 . ± 1 2 km s . -1 Mpc -1 ; there is a tension at the 2.4 σ level between the two. One way of reducing the tension is to increase N eff [4]. Increasing the radiation density at fixed θ ∗ (to preserve the angular scales of the acoustic peaks) and fixed z eq (to preserve the early-ISW effect and so first-peak height) leads to the increase of the Hubble expansion rate before recombination and thus the decrease of the age of the universe. The angular scale of the photon diffusion length, θ D , is thus increased, thereby reducing the power in the CMB damping tail at a given multipole. The result is that N eff is positively correlated with H 0 .

Therefore, increasing both ∑ m ν and N eff could reconcile all the mentioned tensions between Planck and other astrophysical observations, and on the other hand the combination of these consistent data sets could determine these parameters more accurately. But in this case one has to assume the model in which the massive active neutrinos coexist with some extra radiation degrees of freedom. Actually, a more natural choice is to consider the model with light sterile neutrinos.

The existence of light massive sterile neutrinos is hinted by the anomalies of shortbaseline neutrino experiments [9-13]. There is evidence for oscillations at a ∆ m 2 of about 1 eV 2 from these experiments. Thus, additional neutrino masses are required. And the additional types of neutrino have to be sterile neutrinos so that they do not interact by the weak interaction and do not affect the width of Z 0 . Through their mixing with active neutrinos and their gravitational interactions, the sterile neutrinos could have a considerable effect on astrophysics and cosmology. Therefore, it is particularly important to search for cosmological evidence of sterile neutrinos through their gravitational effects on the large-scale structure formation and the evolution of the universe. When considering sterile neutrinos in a cosmological model, two extra parameters, N eff and m eff ν, sterile , are added. Hence, in a cosmological model with sterile neutrinos, the tensions of Planck with H 0 direct measurement, counts of Planck SZ clusters, and cosmic shear can be substantially relieved [7, 14, 15].

Recently, the BICEP2 Collaboration reported the detection of the B-mode polarization of the CMB [16]. If the treatment of the foreground model is correct, the BICEP2's result indicates the discovery of the primordial gravitational waves (PGWs). The frontiers of fundamental physics will be pushed forward in an unprecedented way as long as the BICEP2's result is confirmed by upcoming experiments. Adopting the ΛCDM+ r model, the fit to the observed B-mode power spectrum gives an unexpectedly large tensor-to-scalar ratio, r = 0 20 . +0 07 . -0 05 . , with r = 0 disfavored at the 7 σ level [16]. This result is in tension with the upper limit, r &lt; 0 11 (95% CL), given by the fit to the combined Planck+WP+highL data [4]. . To explain and/or reduce this tension, numerous proposals have been put forward; see, e.g., Refs. [17-32]. An intriguing mechanism for relieving this tension is to involve a sterile neutrino species in the cosmological model with PGWs [33-38]. In Ref. [33], the model with PGWs and sterile neutrinos is called ΛCDM+ r + ν s model (or ΛCDM+ + r N eff + m eff ν, sterile model). It

was shown in Refs. [33, 34] that in the ΛCDM+ r + ν s model the tension between Planck and BICEP2 is well relieved, and meanwhile, the other tensions of Planck with other astrophysical observations, such as the H 0 direct measurement, the cluster counts, and the cosmic shear measurement, can all be significantly reduced. In addition, the current consistent cosmological data could provide independent evidence for the existence of sterile neutrino at high statistical significance, i.e., a joint analysis of current data prefer ∆ N eff ≡ N eff -3 046 . &gt; 0 at the 2.7 σ level and a nonzero mass of sterile neutrino at the 3.9 σ level [33]. In Ref. [35], other typical models of neutrino and dark radiation were analyzed with the current observational data in detail.

In this paper, we investigate the effects of dynamical dark energy on the fits of models of neutrinos and dark radiation to the current observations. In Ref. [39], it was demonstrated that the tension between Planck and H 0 direct measurement might hint that dark energy is not the cosmological constant; once a dynamical dark energy is adopted, the tension of Planck with H 0 measurement could be well relieved. Also, it is known that due to the repulsive gravitational force, dark energy could suppress the growth of large-scale structure. Therefore, we wish to see if dark energy together with neutrinos could further reduce the tensions between Planck and other astrophysical observations. On the other hand, we will combine the consistent data sets to constrain the cosmological parameters concerning dark energy and neutrinos.

The paper is organized as follows. In Sec. 2, we briefly present the cosmological models and the observational data we use in this work. In Sec. 3, we test the data consistency for the models and present the fit results. Conclusion is given in Sec. 4. More detailed fit results are given in Appendix A.

## 2 Methodology

In this paper, we consider the extension of the ΛCDM cosmology to models in which dark energy has a constant w , i.e., the w CDM cosmology. Of course, since the PGWs are likely to have been detected by BICEP2, the basic framework should be the w CDM+ model. Further r extensions considered in this paper include: (i) active neutrinos with additional parameter ∑ m ν , (ii) extra dark radiation with additional parameter N eff , (iii) active neutrinos coexisting with extra relativistic degrees of freedom, with additional parameters ∑ m ν and N eff , and (iv) massive sterile neutrinos with additional parameters N eff and m eff ν, sterile .

The conventions used in this paper are consistent with those adopted by the Planck team [4], i.e., those used in the camb Boltzmann code. The base parameters for the basic eight-parameter w CDM+ r model are:

$$\{ \omega _ { b }, \, \omega _ { c }, \, 1 0 0 \theta _ { \text{MC} }, \, \tau, \, w, \, n _ { s }, \, \ln ( 1 0 ^ { 1 0 } A _ { s } ), \, r _ { 0. 0 5 } \},$$

where ω b ≡ Ω b h 2 and ω c ≡ Ω c h 2 are the baryon and cold dark matter densities today, respectively, θ MC is the approximation used in CosmoMC to r s ( z ∗ ) /D A ( z ∗ ) (the angular size of the sound horizon at the time of last-scattering), τ is the Thomson scattering optical depth due to reionization, w is the dark energy equation of state parameter, n s and A s are the power-law spectral index and power of the primordial curvature perturbations, respectively, and r 0 05 . is the tensor-to-scalar ratio at the pivot scale k 0 = 0 05 Mpc . -1 . Flat priors for the base parameters are used. Note also that the prior ranges for the base parameters are chosen to be much wider than the posterior in order not to affect the results of parameter

estimation. We use the CosmoMC package [40] to infer the posterior probability distributions of parameters.

Next, we describe the observational data sets used in this paper. Actually, this work is an extension of our previous works [33, 35] to the models of dark energy with constant w , thus we will use the same data sets to Refs. [33, 35] in order to be easier for a direct comparison. The data sets we use include the CMB (Planck+WP), BAO, H 0 , SZ, Lensing (CMB lensing + weak lensing), and BICEP2.

Planck+WP. We use the CMB TT angular power spectrum data from the first release of Planck [4], combined with the CMB large-scale TE and EE polarization power spectrum data form the 9-yr release of WMAP [3]. Note that in some occasions we also abbreviate Planck+WP to CMB for convenience.

BAO. We use the latest BAO measurements from the data release (DR) 11 of the Baryon Oscillation Spectroscopic Survey (BOSS) (part of SDSS-III): D V (0 32)( . r d, fid /r d ) = (1264 ± 25) Mpc and D V (0 57)( . r d, fid /r d ) = (2056 ± 20) Mpc, with r d, fid = 149 28 Mpc [41]. . Note also that there are some other BAO data sets, e.g., 6dFGS ( z = 0 1) [42], SDSS-DR7 . ( z = 0 35) [43], and WiggleZ ( . z = 0 44, 0.60, and 0.73) [44], where the three data from the . WiggleZ survey are correlated (with the inverse covariance matrix given by Ref. [44]). This work is in accordance with our previous papers [33, 35] for the use of the BAO data, i.e., we only use the latest two most accurate BAO data from the BOSS-DR11. This is sufficient for our purpose in breaking the CMB parameter degeneracies.

H . 0 We use the direct measurement of the Hubble constant from the Hubble Space Telescope (HST) observations, H 0 = (73 8 . ± 2 4) km s . -1 Mpc -1 [8].

SZ. We use the result of the counts of clusters of galaxies from the sample of Planck SZ clusters, σ 8 (Ω m / 0 27) . 0 3 . = 0 782 . ± 0 010 [5]. . Note that this result is derived by using the mass function given by Tinker et al. [45]; a different mass function given by Watson et al. [46] leads to a slightly different value, i.e., σ 8 (Ω m / 0 27) . 0 3 . = 0 802 . ± 0 014. . Moreover, the result also depends on the bias (1 -b ) that is assumed to account for all the possible observational biases including departure from hydrostatic equilibrium, absolute instrument calibration, temperature inhomogeneities, residual selection bias, etc. Numerical simulations taking into account of several ingredients of gas physics of clusters give the result of the bias of (1 -b ) = 0 8 . +0 2 . -0 1 . . Adopting the central value, (1 -b ) = 0 8, the Planck team found . that the constraints on Ω m and σ 8 are in good agreement with previous measurements using clusters of galaxies [5]. The result of σ 8 (Ω m / . 0 27) 0 3 . quoted in this paper is derived by fixing (1 -b ) = 0 8. . The result is changed to σ 8 (Ω m / . 0 27) 0 3 . = 0 764 . ± 0 025, if the bias is . allowed to vary in the range of [0 7 . , 1]. In Table 2 of Ref. [5], other values of σ 8 (Ω m / . 0 27) 0 3 . from various data combinations and analysis methods are also given. But in this paper, in accordance with previous works by us and by other authors, we choose to use the result of σ 8 (Ω m / . 0 27) 0 3 . = 0 78 . ± 0 01. .

Lensing. We use two kinds of lensing data, i.e., the CMB lensing power spectrum C φφ glyph[lscript] from the Planck mission [47] and the cosmic shear measurement of weak lensing from the CFHTLenS survey, σ 8 (Ω m / . 0 27) 0 46 . = 0 774 . ± 0 040 [6]. . Of course, they are two absolutely different, physically and observationally independent data sets.

BICEP2. We use the CMB angular power spectra (TT, TE, EE, and BB) data from the BICEP2 [16].

Since it has been proven that the Planck data are in good agreement with the BAO data, we can alway safely combine Planck+WP with BAO. So the Planck+WP+BAO combination is the basic data combination used in this paper. In the six-parameter base ΛCDM model, it

was shown that the Planck data are in tension with the H 0 direct measurement, the cluster counts, and the cosmic shear measurement at the 2-3 σ level [4]. In the seven-parameter ΛCDM+ r model, the Planck temperature data are also in tension with the BICEP2 polarization data [16]. It has been demonstrated that the ingredients such as dark energy, dark radiation, and neutrinos could help relieve these tensions to some extent [33-39]. Therefore, in this paper, we will test the data consistency in the cosmological models with dark energy, dark radiation, and neutrinos.

As mentioned above, we consider four models in this paper, i.e., (i) the w CDM+ + r ∑ m ν model, (ii) the w CDM+ + r N eff model, (iii) the w CDM+ + r ∑ m ν + N eff model, and (iv) the w CDM+ + r N eff + m eff ν, sterile model. We will constrain these models by using the Planck + WP + BAO data, and then compare the derived results with the observations of H 0 , SZ cluster counts, cosmic shear, and tensor-to-scalar ratio. This could test if in these models the Planck+WP+BAO data are consistent with these astrophysical observations. If these data sets are consistent, we can then combine these data together and use them to constrain the parameters concerning dark energy, dark radiation, neutrinos, and so on. But if the data are still in tension, the direct combination of different data sets is not appropriate. The subsequent data fits are based on this data consistency test analysis, and the fit results will be discussed in detail.

## 3 Results and discussion

## 3.1 Data consistency tests

We first constrain the four models using the Planck+WP+BAO data. Figure 1 shows the fit results in the Ω m -H 0 , Ω m -σ 8 , and n s -r 0 002 . planes. To exhibit the effects of dynamical dark energy, we also compare the results of w CDM + neutrino models with those of ΛCDM + neutrino models (purple contours vs. blue contours). The first-row panels are for the ΛCDM/ w CDM+ + r ∑ m ν models, the second-row panels are for the ΛCDM/ w CDM+ + r N eff models, the third-row panels are for the ΛCDM/ w CDM+ + r ∑ m ν + N eff models, and the fourth-row panels are for the ΛCDM/ w CDM+ + r N eff + m eff ν, sterile models. The gray bands are from the observational results, i.e., the HST result of Hubble constant H 0 = 73 8 . ± 2 4 km . s -1 Mpc -1 [8], the Planck result of SZ cluster counts σ 8 (Ω m / . 0 27) 0 3 . = 0 782 . ± 0 010 [5], and . the BICEP2 result of tensor-to-scalar ratio r = 0 20 . +0 07 . -0 05 . [16]. From the comparison of the contours from CMB+BAO with the bands from other astrophysical observations, one can directly test the consistency between different data sets.

First, we discuss the results of the ΛCDM/ w CDM+ + r ∑ m ν models (see the first-row panels). From the Ω m -H 0 plane, we clearly see that the tension on H 0 appearing in the ΛCDM-based model is greatly relieved once a dynamical dark energy is considered. But even in the w CDM-based model, Ω m and H 0 are still in strong degeneracy. Since w is anticorrelated with H 0 , we infer that a larger H 0 may prefer a phantom energy with w &lt; -1. From the Ω m -σ 8 plane, we see that the consideration of dynamical dark energy amplifies the parameter space but nearly does not improve the reduction of the tension in the Ω m -σ 8 plane. Since w is anti-correlated with σ 8 but positively correlated with Ω m , the effect of w on σ 8 (Ω m / . 0 27) 0 3 . would be canceled a lot, leading to the fact that only ∑ m ν plays a significant role in relieving the tension between Planck CMB and SZ cluster counts. From the n s -r 0 002 . plane, we see that w could not affect the value of tensor-to-scalar ratio.

Next, we discuss the results of the ΛCDM/ w CDM+ + r N eff models (see the second-row panels). In the base ΛCDM model, the Planck prediction of H 0 is in tension with the direct

Figure 1 . The Planck+WP+BAO constraints in the Ω m -H 0 , Ω m -σ 8 , and n s -r 0 002 . planes for the ΛCDM/ w CDM+ + r ∑ m ν models (first-row panels), the ΛCDM/ w CDM+ + r N eff models (second-row panels), the ΛCDM/ w CDM+ + r ∑ m ν + N eff models (third-row panels), and the ΛCDM/ w CDM+ + r N eff + m eff ν, sterile models (fourth-row panels). The gray bands stand for the observational results, i.e., the HST result of Hubble constant H 0 = 73 8 . ± 2 4 km s . -1 Mpc -1 [8], the Planck result of SZ cluster counts σ 8 (Ω m / 0 27) . 0 3 . = 0 782 . ± 0 010 [5], . and the BICEP2 result of tensor-to-scalar ratio r = 0 20 . +0 07 . -0 05 . [16].

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0481v3-figure-3058b3c7-0f38-4086-ad77-6f3dbce980ad.jpg)

measurement of H 0 at about the 2.5 σ level, but once N eff is considered to be free, this tension could be relieved greatly, as shown by the blue contours in the Ω m -H 0 plane for the ΛCDM+ + r N eff model. Further introducing w does not lead to more impact on H 0 (as shown by the purple contours in the Ω m -H 0 plane), from which one can infer that the fit of this model to the CMB+BAO data would give the result of w around -1. In the Ω m -σ 8 plane, we can see that the parameter N eff does not help reduce the tension between Planck data and cluster data, but the parameter w could play a significant role in relieving this tension. Due to the repulsive gravitational force of dark energy, the growth of large-scale structure is suppressed by dark energy. Contrary to our intuition, larger w gives rise to more suppression

Table 1 . Predictions of the values of H 0 , σ 8 (Ω m / 0 27) . 0 3 . , σ 8 (Ω m / 0 27) . 0 46 . , and r 0 002 . given by the ΛCDM+ , ΛCDM+ + r r N eff + m eff ν, sterile , w CDM+ , and r w CDM+ + r N eff + m eff ν, sterile models under the constraints from Planck+WP+BAO. Quantified levels of the tensions with the observational results are also given.

|                           | ΛCDM+ r                    | ΛCDM+ r   | ΛCDM+ r + N eff + m eff ν, sterile   | ΛCDM+ r + N eff + m eff ν, sterile   | w CDM+ r                   | w CDM+ r   | w CDM+ r + N eff + m eff ν, sterile   | w CDM+ r + N eff + m eff ν, sterile   |
|---------------------------|----------------------------|-----------|--------------------------------------|--------------------------------------|----------------------------|------------|---------------------------------------|---------------------------------------|
| Parameter                 | 68% limits                 | tension   | 68% limits                           | tension                              | 68% limits                 | tension    | 68% limits                            | tension                               |
| H 0                       | 67 . 80 +0 . 64 - 0 . 63   | 2 . 4 σ   | 70 . 8 +1 . 7 - 2 . 1                | 1 . 0 σ                              | 69 . 0 +1 . 8 - 2 . 6      | 1 . 6 σ    | 71 . 9 +2 . 3 - 3 . 2                 | 0 . 6 σ                               |
| σ 8 (Ω m / 0 . 27) 0 . 3  | 0 . 857 ± 0 . 015          | 4 . 3 σ   | 0 . 842 +0 . 038 - 0 . 029           | 2 . 0 σ                              | 0 . 868 ± 0 . 026          | 3 . 2 σ    | 0 . 840 +0 . 039 - 0 . 034            | 1 . 7 σ                               |
| σ 8 (Ω m / 0 . 27) 0 . 46 | 0 . 876 +0 . 019 - 0 . 018 | 2 . 3 σ   | 0 . 858 +0 . 038 - 0 . 030           | 1 . 7 σ                              | 0 . 882 +0 . 024 - 0 . 023 | 2 . 3 σ    | 0 . 853 +0 . 041 - 0 . 033            | 1 . 5 σ                               |
| r 0 . 002                 | < 0 . 12 (95%)             |           | < 0 . 20 (95%)                       |                                      | < 0 . 11 (95%)             |            | < 0 . 20 (95%)                        |                                       |

of the structure growth; in other words, larger w leads to smaller σ 8 (or, w and σ 8 are in the anti-correlation). The reason is that for larger w and fixed present dark energy density, dark energy comes to dominate earlier, causing the time of suppressing the growth of linear matter perturbation to be longer. Since w is positively correlated with Ω m , larger w would lead to larger Ω m . The effects of w on σ 8 and Ω m are contrary, and thus some offset happens for the quantity σ 8 (Ω m / . 0 27) 0 3 . , resulting in that increasing w still could not well relieve the tension (as shown by the purple contours in the Ω m -σ 8 plane). Since N eff is positively correlated with n s , and w is also positively correlated with n s , they can both contribute to enhance the value of r , but it seems that even in the w CDM+ + r N eff model the tension of r is still not well relieved.

At last, we discuss the results of both the ΛCDM/ w CDM+ + r ∑ m ν + N eff models and the ΛCDM/ w CDM+ + r N eff + m eff ν, sterile models (see the third-row and the last-row panels). Since the results of the two cases are very similar (see also Ref. [35]), we only focus on the case of sterile neutrino. We can see clearly that in the models involving sterile neutrinos all the tensions could be well relieved; see also Refs. [33-38]. In addition, when w is considered, the tensions will be further reduced, as shown by the Ω m -H 0 and the Ω m -σ 8 panels. In order to quantify how the data consistencies are improved by the sterile neutrino and dynamical dark energy, we give the fit results of H 0 , σ 8 (Ω m / . 0 27) 0 3 . , σ 8 (Ω m / . 0 27) 0 46 . , and r 0 002 . for the ΛCDM+ , the ΛCDM+ + r r N eff + m eff ν, sterile , the w CDM+ , and the r w CDM+ + r N eff + m eff ν, sterile models in Table 1. The levels of tensions between fit results and observational results are also given. Under the constraints from Planck+WP+BAO, considering sterile neutrino improves the tension with H 0 from 2.4 σ to 1.0 σ , the tension with SZ cluster counts from 4.3 σ to 2.0 σ , and the tension with cosmic shear from 2.3 σ to 1.7 σ , respectively; considering dynamical dark energy (with constant w ) improves the above levels of tensions to 1.6 σ , 3.2 σ , and 2.3 σ , respectively; simultaneous consideration of sterile neutrino and dynamical dark energy improves the tension levels to 0.6 σ , 1.7 σ , and 1.5 σ , respectively. But in this case dynamical dark energy nearly does not affect the fit result of r 0 002 . .

From the above analysis, we find that considering dynamical dark energy with constant w could further improve the data consistency. In the w CDM+ + r ∑ m ν model, the CMB+BAO data are basically consistent with H 0 , SZ cluster counts, and cosmic shear data, but still in tension with the BICEP2 data. In the w CDM+ + r N eff model, the CMB+BAO data are consistent with the H 0 measurement, but still in some tension with cluster data and BICEP2 data (though the previous tensions are further reduced to some extent). In the w CDM+ + r ∑ m ν + N eff model and the w CDM+ + r N eff + m eff ν, sterile model, the CMB+BAO data are consistent with all the data sets considered, and thus these data sets can be combined together to perform joint constraints on these two models almost without any question.

CDM+r+

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0481v3-figure-bbc034ab-b225-4a7a-86e2-fe478e1b61d9.jpg)

Σ

mν

Figure 2 . One-dimensional posterior distributions for the w CDM+ + r ∑ m ν model (upper) and the w CDM+ + r N eff model (lower).

## 3.2 Constraints on dark energy, dark radiation, and neutrinos

In this subsection, we give the fit results for the four models. Of course, we focus on the constraint results of dark energy, dark radiation, and neutrinos.

Figure 2 shows the fit results for the w CDM+ + r ∑ m ν model and the w CDM+ + r N eff model. In this figure, the one-dimensional posterior possibility distributions of w and ∑ m ν are plotted for the w CDM+ + r ∑ m ν model, and those of w and N eff are plotted for the w CDM+ + r N eff model. Black curves are from the Planck+WP+BAO constraints, and red curves are from the Planck+WP+BAO+ H 0 +SZ+Lensing constraints.

For the w CDM+ + r ∑ m ν model, the CMB+BAO data combination gives the constraint results:

$$w = - 1. 1 4 _ { - 0. 1 1 } ^ { + 0. 1 7 }, \\ \sum m _ { \nu } < 0. 4 8 \ e V.$$

Here, we show the ± 1 σ errors for w , and show the 95% CL upper limit for ∑ m ν . Note that throughout the paper we quote ± 1 σ errors, but when the parameters cannot be well constrained we only give the 2 σ upper limits for these parameters. The constraints on Ω m and H 0 are: Ω m = 0 293 . +0 021 . -0 018 . and H 0 = 70 1 . +2 3 . -3 2 . km s -1 Mpc -1 . Since the CMB+BAO data are basically consistent with the other astrophysical data (note that in this subsection we do not discuss BICEP2), one can also combine these data together. For convenience, hereafter we also use 'other' to denote the data combination of H 0 +SZ+Lensing. The

CMB+BAO+other data combination leads to:

$$w = - 1. 3 4 ^ { + 0. 1 0 } _ { - 0. 1 2 }, \\ \sum m _ { \nu } = 0. 5 8 ^ { + 0. 1 1 } _ { - 0. 1 0 } \, e V.$$

One can clearly see that once the other data sets are combined, the value of w becomes smaller. In fact, even though the consideration of dynamical dark energy could relieve the tension of H 0 , the strong degeneracy between Ω m and H 0 still exists; recall the discussion in the last subsection. So, when the other data sets are added, the Ω m -H 0 degeneracy is partly broken, resulting in that the Ω m becomes smaller and H 0 becomes larger: Ω m = 0 272 . +0 013 . -0 014 . and H 0 = 73 3 . ± 2 0 km s . -1 Mpc -1 . Since w and H 0 are in anti-correlation, larger H 0 leads to smaller w . This is why w becomes smaller in this case. In addition, the neutrino mass is sensitive to the cluster data, and so the combination with SZ and Lensing data could tighten the constraint on ∑ m ν .

For the w CDM+ + r N eff model, the CMB+BAO data combination gives:

$$w = - 0. 9 6 ^ { + 0. 1 2 } _ { - 0. 0 9 }, \\ N _ { \text{eff} } = 3. 6 3 ^ { + 0. 3 8 } _ { - 0. 4 2 }.$$

In this case, we also have Ω m = 0 306 . +0 017 . -0 015 . and H 0 = 70 1 . +2 0 . -2 4 . km s -1 Mpc -1 . From the discussion in the last subsection, we learn that the CMB+BAO data are not well consistent with the SZ cluster counts data, and so in principle it is not appropriate to combine them together. If we combine these data together (CMB+BAO+other), we get the results:

$$w = - 0. 8 5 ^ { + 0. 0 8 } _ { - 0. 0 6 }, \\ N _ { \text{eff} } = 3. 6 3 ^ { + 0. 3 9 } _ { - 0. 4 3 }.$$

We must keep in mind that the above results are from the tensioned data combination and thus are not reliable. Since the SZ cluster counts data give a lower σ 8 , we obtain a much higher w in this case (recall that w and σ 8 are in anti-correlation; see the relevant discussion in the last subsection). N eff is insensitive to the cluster and cosmic shear data, and so in this case the constraint on N eff is not changed. In this case, we have Ω m = 0 306 . ± 0 012 and . H 0 = 68 7 . ± 1 5 km s . -1 Mpc -1 .

Figure 3 summarizes the main results for the w CDM+ + r ∑ m ν + N eff model and the w CDM+ + r N eff + m eff ν, sterile model. In this figure, the two-dimensional posterior possibility distribution contours in the Ω m -w and ∑ m ν -N eff planes are plotted for the active neutrino model, and those in the Ω m -w and m eff ν, sterile -N eff planes are plotted for the sterile neutrino model. Purple contours are from the CMB+BAO constraints, and blue contours are from the CMB+BAO+other constraints. As discussed in the last subsection, in these two models the CMB+BAO data are consistent with the other observations considered in this paper, and thus the combination of these data sets is viewed as appropriate. In the following we report the fit results for the two models.

For the w CDM+ + r ∑ m ν + N eff model, the CMB+BAO data give the fit results:

$$w = - 1. 0 5 ^ { + 0. 1 7 } _ { - 0. 1 1 }, \\ \sum m _ { \nu } < 0. 6 5 1 \ e V, \quad N _ { \text{eff} } = 3. 6 8 ^ { + 0. 3 8 } _ { - 0. 4 5 }.$$

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0481v3-figure-de93e2d9-7dd2-4411-ba55-a827bf974352.jpg)

0.75

1.00

1.25

1.50

w

0.24

0.27

w

CDM+r+

0.30

mν, eff

sterile

0.33

eff

0.36

5.0

4.5

4.0

3.5

0.00

5.0

4.5

4.0

3.5

0.00

0.25

0.25

w

CDM+r+

0.50

Σ

mν

Σ

mν

0.75

[eV]

CDM+r+

w

0.50

mν,

0.75

+

N

eff

1.00

eff sterile

+

N

1.00

Ωm mν,

eff sterile

[eV]

Figure 3 . Two-dimensional joint, marginalized constraints for the w CDM+ + r ∑ m ν + N eff model (upper) and the w CDM+ + r N eff + m eff ν, sterile model (lower).

w CDM+r+ Σ mν + N eff In this case, we have Ω m = 0 299 . +0 021 . -0 016 . and H 0 = 71 6 . +2 4 . -3 2 . km s -1 Mpc -1 . Furthermore, the CMB+BAO+other data give the fit results:

$$w = - 1. 2 2 ^ { + 0. 1 4 } _ { - 0. 1 2 }, \\ \sum m _ { \nu } = 0. 6 9 ^ { + 0. 1 3 } _ { - 0. 1 5 } \, \text{eV}, \quad N _ { \text{eff} } = 3. 6 6 ^ { + 0. 3 5 } _ { - 0. 4 1 }.$$

In this case, we have Ω m = 0 283 . ± 0 015 and . H 0 = 74 1 . ± 2 0 km s . -1 Mpc -1 . From the upper panels of Fig. 3, one can see that the CMB+BAO data are well consistent with the other astrophysical data for this model.

For the w CDM+ + r N eff + m eff ν, sterile model, the CMB+BAO data combination gives the constraint results:

eff sterile

$$w = - 1. 0 6 ^ { + 0. 1 6 } _ { - 0. 1 1 },$$

mν,

+

N

eff

$$N _ { \text{eff} } = 3. 6 8 ^ { + 0. 3 8 } _ { - 0. 4 5 }, \quad m ^ { \text{eff} } _ { \nu, \text{sterile} } < 0. 7 0 \ e V.$$

In this case, we have Ω m = 0 298 . +0 020 . -0 017 . and H 0 = 71 9 . +2 3 . -3 2 . km s -1 Mpc -1 . Furthermore, the CMB+BAO+other data give the constraint results:

$$w = - 1. 2 0 \pm 0. 1 0, \\ N _ { \text{eff} } = 3. 5 6 ^ { + 0. 1 5 } _ { - 0. 3 3 }, \quad m _ { \nu, \text{steri} } ^ { \text{eff} } = 0. 7 0 \pm 0. 1 6 \ e V.$$

In this case, we have Ω m = 0 283 . +0 013 . -0 014 . and H 0 = 73 7 . ± 2 0 km s . -1 Mpc -1 . The same to the above active neutrino model, we also find from the lower panel of Fig. 3 that the CMB+BAO data are well consistent with the other observations for the sterile neutrino model.

+

N

eff

N

eff

N

eff

r

Figure 4 . Two-dimensional joint, marginalized constraints in the r 0 002 . -n s plane for the four models.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0481v3-figure-9725ee0c-076e-4ce6-8290-abe037645009.jpg)

We find that the w CDM-based neutrino models fit the cosmological data much better than the corresponding ΛCDM-based models. The χ 2 min values of the w CDM-based models are given in Tables 2-5 in Appendix A, and those of the corresponding ΛCDM-based models can be found in Tables 1-4 of Ref. [35]. Note that in these tables actually the values of -ln L max = χ 2 min / 2 are given. Fits to the CMB+BAO+other data lead to an increase of ∆ χ 2 = -10 54 for the . w CDM+ + r ∑ m ν model, ∆ χ 2 = -3 64 for . the w CDM+ + r N eff model, ∆ χ 2 = -4 28 . for the w CDM+ + r ∑ m ν + N eff model, and ∆ χ 2 = -7 5 . for the w CDM+ + r N eff + m eff ν, sterile model, compared to the corresponding ΛCDM-based models. According to the Akaike information criterion, if χ 2 min improves by 2 or more with one additional parameter, its incorporation is justified. In this context, the performance of w CDM-based models (with only one more additional parameter beyond the corresponding ΛCDM-based models) improves χ 2 min significantly.

Furthermore, to see how the additional neutrino/dark radiation parameters improve the fits in the framework of w CDM, we report the χ 2 value of fitting the w CDM+ r model to the CMB+BAO+other data, χ 2 min = 9841 00. . Therefore, compared to this model, we get ∆ χ 2 = -18 12 for the . w CDM+ + r ∑ m ν model, ∆ χ 2 = -1 36 for the . w CDM+ + r N eff model, ∆ χ 2 = -18 42 for . the w CDM+ + r ∑ m ν + N eff model, and ∆ χ 2 = -22 02 for . the w CDM+ + r N eff + m eff ν, sterile model. It is clear to see that the neutrino mass plays a more important role than N eff in the fits.

## 3.3 Early-universe parameters and BICEP2

In this subsection we consider the constraints on n s and r and the BICEP2 data. Figure 4 shows the two-dimensional joint constraints on the four model in the n s -r 0 002 . plane. The

orange contours are for the CMB+BAO data, the purple ones are for the CMB+BAO+other data, and the blue ones are for the CMB+BAO+other+BICEP2 data.

For the w CDM+ + r ∑ m ν model and the w CDM+ + r N eff model, since the tension between Planck and BICEP2 cannot be greatly reduced, we are not interested in these two models for this aspect. We thus only give the 95% CL upper limits of r 0 002 . given by the CMB+BAO data for these two models: r 0 002 . &lt; 0 125 for the . w CDM+ + r ∑ m ν model and r 0 002 . &lt; . 0 165 for the w CDM+ + r N eff model.

In the w CDM+ + r ∑ m ν + N eff model and the w CDM+ + r N eff + m eff ν, sterile model, the tension of r between Planck and BICEP2 can be well relieved, and so we are interested in these two cases. For the w CDM+ + r ∑ m ν + N eff model, the CMB+BAO data give:

$$n _ { s } = 0. 9 8 3 ^ { + 0. 0 1 5 } _ { - 0. 0 1 7 }, \quad r _ { 0. 0 0 2 } < 0. 2 0 0.$$

Furthermore, the CMB+BAO+other+BICEP2 give the results:

$$n _ { s } = 0. 9 8 9 \pm 0. 0 1 3, \quad r _ { 0. 0 0 2 } = 0. 2 0 0 ^ { + 0. 0 3 9 } _ { - 0. 0 4 9 }.$$

For the w CDM+ + r N eff + m eff ν, sterile model, the CMB+BAO data give:

$$n _ { s } = 0. 9 8 4 ^ { + 0. 0 1 4 } _ { - 0. 0 1 7 }, \quad r _ { 0. 0 0 2 } < 0. 1 9 9.$$

Furthermore, the CMB+BAO+other+BICEP2 give the results:

$$n _ { s } = 0. 9 8 7 ^ { + 0. 0 1 4 } _ { - 0. 0 1 3 }, \quad r _ { 0. 0 0 2 } = 0. 1 9 3 ^ { + 0. 0 3 8 } _ { - 0. 0 5 0 }.$$

At last, we discuss the goodness of fits to the full data combination (with BICEP2 involved) for these two models. To see the role the parameter w plays in the fits, we take the ΛCDM-based models as references, and we obtain ∆ χ 2 = -1 78 for the . w CDM+ + r ∑ m ν + N eff model, and ∆ χ 2 = -3 04 . for the w CDM+ + r N eff + m eff ν, sterile model. Furthermore, to see how the two additional parameters relevant to neutrino/dark radiation improve the fits in the framework of w CDM, we choose the w CDM+ r model as a reference. In this case, we have χ 2 min = 9884 78 for the . w CDM+ r model. Therefore, we get ∆ χ 2 = -20 44 for . the w CDM+ + r ∑ m ν + N eff model, and ∆ χ 2 = -20 18 for the . w CDM+ + r N eff + m eff ν, sterile model.

## 4 Conclusion

Based on the standard ΛCDM cosmology, the Planck data are in good agreement with the BAO data, but are in tension with other astrophysical observations, such as the H 0 direct measurement, the SZ cluster counts, the cosmic shear measurement, and so on [4]. Some ingredients such as dark radiation, massive neutrinos, and dynamical dark energy were proposed to be considered in the cosmological model to reduce the tensions. In particular, extra sterile neutrino species, with additional parameters N eff and m eff ν, sterile , was used to resolve the tensions between Planck and other astrophysical observations. Then, after the detection of the B-mode polarization of the CMB by the BICEP2 experiment, which implies that the PGWs were possibly discovered, the base standard cosmology should at least be extended to the ΛCDM+ r model [16]. However, it was found that the BICEP2's result of r is in tension with the Planck's fit result. Again, it was proposed that the sterile neutrino can be used to reconcile the PGW results from BICEP2 and Planck [33, 34]. Other cases concerning neutrinos and dark radiation were also investigated in detail [35].

In this paper, we study the neutrino cosmological models in which the cosmological constant is replaced with the dynamical dark energy with constant w . Four cases are considered, i.e., the w CDM+ + r ∑ m ν model, the w CDM+ + r N eff model, the w CDM+ + r ∑ m ν + N eff model, and the w CDM+ + r N eff + m eff ν, sterile model. The observational data we consider in this paper include the Planck+WP, BAO, H 0 , Planck SZ cluster, Planck CMB lensing, cosmic shear, and BICEP2 data. We first tested the consistency of these data sets in the four cosmological models, and then performed joint constraints on the properties of dark energy and neutrinos.

For the w CDM+ + r ∑ m ν model, the consideration of dynamical dark energy is rather helpful in relieving the tension of H 0 , but even in this model Ω m and H 0 are still in strong degeneracy; the dark energy parameter w does not help much to reduce the tension between Planck and SZ cluster counts, and so only ∑ m ν plays a significant role in this aspect; also, w could not affect the value of r . Thus, the CMB+BAO data are basically consistent with other observations (except for BICEP2). The CMB+BAO constraint gives: w = -1 14 . +0 17 . -0 11 . and ∑ m &lt; ν 0 48 eV. The CMB+BAO+ . H 0 +SZ+Lensing constraint gives: w = -1 34 . +0 10 . -0 12 . and ∑ m ν = 0 58 . +0 11 . -0 10 . eV.

For the w CDM+ + r N eff model, the dark energy parameter w does not lead to more impact on H 0 , but plays a significant role in reducing the tension of SZ cluster counts; larger w leads to smaller σ 8 and larger Ω m , so some offset happens for the quantity σ 8 (Ω m / 0 27) . 0 3 . , resulting in that the consideration of w still could not well relieve the tension between Planck and clusters; w can enhance the value of r in this case (it is positively correlated with n s ), but the tension of r still cannot be well relieved. The CMB+BAO data combination gives the results: w = -0 96 . +0 12 . -0 09 . , N eff = 3 63 . +0 38 . -0 42 . , and r 0 002 . &lt; 0 165. .

The results of the active neutrinos plus dark radiation model and the sterile neutrino model are very similar. So we only took the sterile neutrino model as an example to analyze the data consistency. We found that in this model all the tensions considered are greatly reduced. Under the constraint from CMB+BAO, for the ΛCDM+ r model, the tensions with H 0 , SZ clusters, and cosmic shear are 2.4 σ , 4.3 σ , and 2.3 σ , respectively; for the w CDM+ + r N eff + m eff ν, sterile model, the above tensions are reduced to 0.6 σ , 1.7 σ , and 1.5 σ . The sterile neutrinos could well reconcile the r results from BICEP2 and Planck, and in this case dynamical dark energy does not further improve the r result.

For the w CDM+ + r ∑ m ν + N eff model, the CMB+BAO+ H 0 +SZ+Lensing constraint gives: w = -1 22 . +0 14 . -0 12 . , ∑ m ν = 0 69 . +0 13 . -0 15 . eV, and N eff = 3 66 . +0 35 . -0 41 . ; further adding the BICEP2 data, we have n s = 0 989 . ± 0 013 and . r 0 002 . = 0 200 . +0 039 . -0 049 . . For the w CDM+ + r N eff + m eff ν, sterile model, the CMB+BAO+ H 0 +SZ+Lensing data combination leads to: w = -1 20 . ± 0 10, . N eff = 3 56 . +0 15 . -0 33 . , and m eff ν, sterile = 0 70 . ± 0 16 eV. . Further adding the BICEP2 data gives: n s = 0 987 . +0 014 . -0 013 . and r 0 002 . = 0 193 . +0 038 . -0 050 . .

From our analysis, we found that both the w CDM+ + r ∑ m ν + N eff model and the w CDM+ + r N eff + m eff ν, sterile model are successful in relieving the tensions between Planck and other astrophysical observations. However, in the w CDM+ + r ∑ m ν + N eff model, one has to assume that the massive active neutrinos coexist with dark radiation, and thus we believe that the sterile neutrino model is more natural. To further tightly constrain the property of dark energy, accurate type Ia supernova data are needed. We leave the further analysis on the sterile neutrinos plus dark energy model in future work.

Table 2 . Fitting results for the w CDM+ + r ∑ m ν model. We quote ± 1 σ errors, but for the parameters that cannot be well constrained, we quote the 95% CL upper limits.

|                | Planck+WP+BAO   | Planck+WP+BAO                    | + H 0 +SZ+Lensing   | + H 0 +SZ+Lensing                | +BICEP2   | +BICEP2                       |
|----------------|-----------------|----------------------------------|---------------------|----------------------------------|-----------|-------------------------------|
| Parameter      | Best fit        | 68% limits                       | Best fit            | 68% limits                       | Best fit  | 68% limits                    |
| Ω b h 2        | 0 . 02236       | 0 . 02202 +0 . 00027 - 0 . 00026 | 0 . 02215           | 0 . 02211 ± 0 . 00024            | 0 . 02200 | 0 . 02205 ± 0 . 00023         |
| Ω c h 2        | 0 . 1180        | 0 . 1193 ± 0 . 0023              | 0 . 1172            | 0 . 1171 +0 . 0013 - 0 . 0012    | 0 . 1176  | 0 . 1169 +0 . 0013 - 0 . 0012 |
| 100 θ MC       | 1 . 04130       | 1 . 04124 ± 0 . 00059            | 1 . 04103           | 1 . 04136 +0 . 00055 - 0 . 00056 | 1 . 04119 | 1 . 04136 ± 0 . 00055         |
| τ              | 0 . 098         | 0 . 090 +0 . 012 - 0 . 014       | 0 . 093             | 0 . 093 +0 . 012 - 0 . 014       | 0 . 092   | 0 . 093 +0 . 012 - 0 . 013    |
| Σ m ν          | 0 . 039         | < 0 . 480                        | 0 . 580             | 0 . 584 +0 . 111 - 0 . 097       | 0 . 545   | 0 . 581 +0 . 115 - 0 . 094    |
| w              | - 1 . 01        | - 1 . 14 +0 . 17 - 0 . 11        | - 1 . 30            | - 1 . 34 +0 . 10 - 0 . 12        | - 1 . 35  | - 1 . 34 +0 . 10 - 0 . 12     |
| n s            | 0 . 9651        | 0 . 9605 ± 0 . 0068              | 0 . 9632            | 0 . 9599 ± 0 . 0059              | 0 . 9623  | 0 . 9621 +0 . 0057 - 0 . 0058 |
| ln(10 10 A s ) | 3 . 100         | 3 . 086 +0 . 024 - 0 . 025       | 3 . 088             | 3 . 087 +0 . 023 - 0 . 025       | 3 . 084   | 3 . 086 +0 . 022 - 0 . 025    |
| r 0 . 05       | 0 . 010         | < 0 . 133                        | 0 . 051             | < 0 . 152                        | 0 . 180   | 0 . 181 +0 . 034 - 0 . 039    |
| Ω Λ            | 0 . 701         | 0 . 707 +0 . 018 - 0 . 021       | 0 . 720             | 0 . 728 +0 . 014 - 0 . 013       | 0 . 730   | 0 . 730 ± 0 . 014             |
| Ω m            | 0 . 299         | 0 . 293 +0 . 021 - 0 . 018       | 0 . 280             | 0 . 272 +0 . 013 - 0 . 014       | 0 . 270   | 0 . 270 ± 0 . 014             |
| σ 8            | 0 . 835         | 0 . 830 ± 0 . 037                | 0 . 780             | 0 . 786 ± 0 . 015                | 0 . 797   | 0 . 787 +0 . 015 - 0 . 016    |
| H 0            | 68 . 6          | 70 . 1 +2 . 3 - 3 . 2            | 72 . 0              | 73 . 3 ± 2 . 0                   | 73 . 5    | 73 . 4 ± 2 . 0                |
| r 0 . 002      | 0 . 0094        | < 0 . 125                        | 0 . 046             | < 0 . 143                        | 0 . 171   | 0 . 173 +0 . 034 - 0 . 042    |
| - ln L max     | 4904.82         | 4904.82                          | 4911.44             | 4911.44                          | 4934.33   | 4934.33                       |

## Acknowledgments

We acknowledge the use of CosmoMC . We thank Yun-He Li for helpful discussion. JFZ is supported by the Provincial Department of Education of Liaoning under Grant No. L2012087. XZ is supported by the National Natural Science Foundation of China under Grant No. 11175042 and the Fundamental Research Funds for the Central Universities under Grant No. N120505003.

## A Detailed constraint results

In this appendix, we give the detailed constraint results for the four cosmological models considered in this paper. The one- and two-dimensional joint, marginalized posterior probability distributions of the parameters for the models are shown in Figs. 5-8. Detailed fit values for the cosmological parameters are given in Tables 2-5. In the tables, the ± 1 σ errors are quoted, but for the parameters that cannot be well constrained, only the 2 σ upper limits are given.

Figure 5 and Table 2 summarize the fit results for the w CDM+ + r ∑ m ν model. Figure 6 and Table 3 summarize the fit results for the w CDM+ + r N eff model. Figure 7 and Table 4 summarize the fit results for the w CDM+ + r ∑ m ν + N eff model. Figure 8 and Table 5 summarize the fit results for the w CDM+ + r N eff + m eff ν, sterile model.

## References

- [1] E. Komatsu et al. [WMAP Collaboration], 'Five-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Cosmological Interpretation,' Astrophys. J. Suppl. 180 , 330 (2009) [arXiv:0803.0547 [astro-ph]].
- [2] E. Komatsu et al. [WMAP Collaboration], 'Seven-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Cosmological Interpretation,' Astrophys. J. Suppl. 192 , 18 (2011) [arXiv:1001.4538 [astro-ph.CO]].
- [3] G. Hinshaw et al. [WMAP Collaboration], 'Nine-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Cosmological Parameter Results,' Astrophys. J. Suppl. 208 , 19 (2013) [arXiv:1212.5226 [astro-ph.CO]].

Figure 5 . Cosmological constraints on the w CDM+ + r ∑ m ν model.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0481v3-figure-2821b031-c92a-474b-8e25-8eac642390f5.jpg)

- [4] P. A. R. Ade et al. [Planck Collaboration], 'Planck 2013 results. XVI. Cosmological parameters,' arXiv:1303.5076 [astro-ph.CO].
- [5] P. A. R. Ade et al. [Planck Collaboration], arXiv:1303.5080 [astro-ph.CO].
- [6] J. Benjamin, L. Van Waerbeke, C. Heymans, M. Kilbinger, T. Erben, H. Hildebrandt, H. Hoekstra and T. D. Kitching et al. , arXiv:1212.3327 [astro-ph.CO].
- [7] R. A. Battye and A. Moss, 'Evidence for massive neutrinos from CMB and lensing observations,' Phys. Rev. Lett. 112 , 051303 (2014) [arXiv:1308.5870 [astro-ph.CO]].
- [8] A. G. Riess, L. Macri, S. Casertano, H. Lampeitl, H. C. Ferguson, A. V. Filippenko, S. W. Jha and W. Li et al. , 'A 3% Solution: Determination of the Hubble Constant with the Hubble

Figure 6 . Cosmological constraints on the w CDM+ + r N eff model.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0481v3-figure-2c24508b-4624-471d-b8ee-8110150c9998.jpg)

Space Telescope and Wide Field Camera 3,' Astrophys. J. 730 , 119 (2011) [Erratum-ibid. 732 , 129 (2011)] [arXiv:1103.2976 [astro-ph.CO]].

- [9] A. Aguilar-Arevalo et al. [LSND Collaboration], 'Evidence for neutrino oscillations from the observation of anti-neutrino(electron) appearance in a anti-neutrino(muon) beam,' Phys. Rev. D 64 , 112007 (2001).
- [10] A. A. Aguilar-Arevalo et al. [MiniBooNE Collaboration], 'Improved Search for ¯ ν µ → ¯ ν e Oscillations in the MiniBooNE Experiment,' Phys. Rev. Lett. 110 , 161801 (2013).
- [11] G. Mention, M. Fechner, T. .Lasserre, T. .A. Mueller, D. Lhuillier, M. Cribier and A. Letourneau, 'The Reactor Antineutrino Anomaly,' Phys. Rev. D 83 , 073006 (2011).

|                | Planck+WP+BAO   | Planck+WP+BAO                    | + H 0 +SZ+Lensing   | + H 0 +SZ+Lensing                | +BICEP2   | +BICEP2                          |
|----------------|-----------------|----------------------------------|---------------------|----------------------------------|-----------|----------------------------------|
| Parameter      | Best fit        | 68% limits                       | Best fit            | 68% limits                       | Best fit  | 68% limits                       |
| Ω b h 2        | 0 . 02230       | 0 . 02245 +0 . 00037 - 0 . 0004  | 0 . 02278           | 0 . 02291 +0 . 00038 - 0 . 00037 | 0 . 02287 | 0 . 02302 ± 0 . 00036            |
| Ω c h 2        | 0 . 1238        | 0 . 127 +0 . 0054 - 0 . 0061     | 0 . 1179            | 0 . 1207 +0 . 0049 - 0 . 0053    | 0 . 1214  | 0 . 1234 ± 0 . 0048              |
| 100 θ MC       | 1 . 04106       | 1 . 04057 +0 . 00073 - 0 . 00072 | 1 . 04095           | 1 . 04114 +0 . 00071 - 0 . 00072 | 1 . 04112 | 1 . 04092 +0 . 00068 - 0 . 00069 |
| τ              | 0 . 098         | 0 . 095 +0 . 013 - 0 . 016       | 0 . 087             | 0 . 091 +0 . 014 - 0 . 016       | 0 . 081   | 0 . 093 +0 . 013 - 0 . 015       |
| w              | - 0 . 974       | - 0 . 962 +0 . 120 - 0 . 087     | - 0 . 886           | - 0 . 847 +0 . 075 - 0 . 059     | - 0 . 832 | - 0 . 819 +0 . 067 - 0 . 053     |
| N eff          | 3 . 4           | 3 . 63 +0 . 38 - 0 . 42          | 3 . 41              | 3 . 63 +0 . 39 - 0 . 43          | 3 . 66    | 3 . 86 +0 . 38 - 0 . 41          |
| n s            | 0 . 975         | 0 . 983 +0 . 016 - 0 . 017       | 0 . 985             | 0 . 995 ± 0 . 016                | 0 . 996   | 1 . 005 ± 0 . 015                |
| ln(10 10 A s ) | 3 . 116         | 3 . 115 ± 0 . 032                | 3 . 082             | 3 . 094 +0 . 032 - 0 . 035       | 3 . 078   | 3 . 104 +0 . 03 - 0 . 034        |
| r 0 . 05       | 0 . 002         | < 0 . 158                        | 0 . 024             | < 0 . 203                        | 0 . 192   | 0 . 188 +0 . 034 - 0 . 042       |
| Ω Λ            | 0 . 694         | 0 . 694 +0 . 015 - 0 . 017       | 0 . 701             | 0 . 694 ± 0 . 012                | 0 . 69    | 0 . 692 ± 0 . 011                |
| Ω m            | 0 . 306         | 0 . 306 +0 . 017 - 0 . 015       | 0 . 299             | 0 . 306 ± 0 . 012                | 0 . 31    | 0 . 308 ± 0 . 011                |
| σ 8            | 0 . 838         | 0 . 838 +0 . 033 - 0 . 037       | 0 . 776             | 0 . 772 ± 0 . 013                | 0 . 763   | 0 . 77 ± 0 . 012                 |
| H 0            | 69 . 3          | 70 . 1 +2 . 0 - 2 . 4            | 68 . 7              | 68 . 7 ± 1 . 5                   | 68 . 4    | 69 . 2 ± 1 . 4                   |
| r 0 . 002      | 0 . 002         | < 0 . 165                        | 0 . 023             | < 0 . 224                        | 0 . 205   | 0 . 207 +0 . 040 - 0 . 054       |
| - ln L max     | 4903.64         | 4903.64                          | 4919.82             | 4919.82                          | 4940.82   | 4940.82                          |

Table 3 . Fitting results for the w CDM+ + r N eff model. We quote ± 1 σ errors, but for the parameters that cannot be well constrained, we quote the 95% CL upper limits.

Table 4 . Fitting results for the w CDM+ + r ∑ m ν + N eff model. We quote ± 1 σ errors, but for the parameters that cannot be well constrained, we quote the 95% CL upper limits.

|                | Planck+WP+BAO   | Planck+WP+BAO                    | + H 0 +SZ+Lensing   | + H 0 +SZ+Lensing                | +BICEP2   | +BICEP2                          |
|----------------|-----------------|----------------------------------|---------------------|----------------------------------|-----------|----------------------------------|
| Parameter      | Best fit        | 68% limits                       | Best fit            | 68% limits                       | Best fit  | 68% limits                       |
| Ω b h 2        | 0 . 02217       | 0 . 02245 +0 . 00038 - 0 . 00041 | 0 . 02244           | 0 . 02253 +0 . 00034 - 0 . 00037 | 0 . 02272 | 0 . 02261 +0 . 00034 - 0 . 00035 |
| Ω c h 2        | 0 . 1253        | 0 . 1274 +0 . 0054 - 0 . 0062    | 0 . 1224            | 0 . 1250 +0 . 0047 - 0 . 0055    | 0 . 1256  | 0 . 1276 +0 . 0048 - 0 . 0047    |
| 100 θ MC       | 1 . 04043       | 1 . 04047 +0 . 00076 - 0 . 00075 | 1 . 04095           | 1 . 04060 +0 . 00071 - 0 . 00072 | 1 . 04089 | 1 . 04036 ± 0 . 00067            |
| τ              | 0 . 096         | 0 . 095 +0 . 014 - 0 . 016       | 0 . 093             | 0 . 099 +0 . 013 - 0 . 016       | 0 . 104   | 0 . 102 +0 . 013 - 0 . 015       |
| Σ m ν          | 0 . 052         | < 0 . 651                        | 0 . 61              | 0 . 69 +0 . 13 - 0 . 15          | 0 . 60    | 0 . 72 +0 . 15 - 0 . 14          |
| w              | - 0 . 99        | - 1 . 05 +0 . 17 - 0 . 11        | - 1 . 19            | - 1 . 22 +0 . 14 - 0 . 12        | - 1 . 11  | - 1 . 18 +0 . 14 - 0 . 11        |
| N eff          | 3 . 51          | 3 . 68 +0 . 38 - 0 . 45          | 3 . 52              | 3 . 66 +0 . 35 - 0 . 41          | 3 . 73    | 3 . 87 ± 0 . 36                  |
| n s            | 0 . 979         | 0 . 983 +0 . 015 - 0 . 017       | 0 . 978             | 0 . 981 +0 . 014 - 0 . 015       | 0 . 988   | 0 . 989 ± 0 . 013                |
| ln(10 10 A s ) | 3 . 115         | 3 . 116 +0 . 031 - 0 . 036       | 3 . 102             | 3 . 119 +0 . 031 - 0 . 035       | 3 . 131   | 3 . 129 +0 . 03 - 0 . 032        |
| r 0 . 05       | 0 . 009         | < 0 . 188                        | 0 . 047             | < 0 . 202                        | 0 . 172   | 0 . 190 +0 . 035 - 0 . 04        |
| Ω Λ            | 0 . 700         | 0 . 701 +0 . 016 - 0 . 021       | 0 . 720             | 0 . 717 ± 0 . 015                | 0 . 709   | 0 . 714 +0 . 015 - 0 . 016       |
| Ω m            | 0 . 300         | 0 . 299 +0 . 021 - 0 . 016       | 0 . 280             | 0 . 283 ± 0 . 015                | 0 . 291   | 0 . 286 +0 . 016 - 0 . 015       |
| σ 8            | 0 . 847         | 0 . 820 ± 0 . 037                | 0 . 771             | 0 . 776 +0 . 015 - 0 . 017       | 0 . 774   | 0 . 773 +0 . 015 - 0 . 016       |
| H 0            | 70 . 3          | 71 . 6 +2 . 4 - 3 . 2            | 73 . 5              | 74 . 1 ± 2 . 0                   | 72 . 9    | 74 . 3 +1 . 9 - 2 . 1            |
| r 0 . 002      | 0 . 009         | < 0 . 200                        | 0 . 044             | < 0 . 213                        | 0 . 177   | 0 . 200 +0 . 039 - 0 . 049       |
| - ln L max     | 4904.45         | 4904.45                          | 4911.29             | 4911.29                          | 4932.17   | 4932.17                          |

Table 5 . Fitting results for the w CDM+ + r N eff + m eff ν, sterile model. We quote ± 1 σ errors, but for the parameters that cannot be well constrained, we quote the 95% CL upper limits.

|                  | Planck+WP+BAO   | Planck+WP+BAO                    | + H 0 +SZ+Lensing   | + H 0 +SZ+Lensing                | +BICEP2   | +BICEP2                       |
|------------------|-----------------|----------------------------------|---------------------|----------------------------------|-----------|-------------------------------|
| Parameter        | Best fit        | 68% limits                       | Best fit            | 68% limits                       | Best fit  | 68% limits                    |
| Ω b h 2          | 0 . 02242       | 0 . 02248 +0 . 00033 - 0 . 00038 | 0 . 02232           | 0 . 02253 +0 . 00028 - 0 . 0003  | 0 . 02267 | 0 . 02261 ± 0 . 00029         |
| Ω c h 2          | 0 . 1242        | 0 . 1275 +0 . 0056 - 0 . 0061    | 0 . 1192            | 0 . 1228 +0 . 0035 - 0 . 005     | 0 . 1267  | 0 . 1265 +0 . 0047 - 0 . 0048 |
| 100 θ MC         | 1 . 04065       | 1 . 04045 +0 . 00075 - 0 . 00074 | 1 . 04134           | 1 . 04075 +0 . 00066 - 0 . 00065 | 1 . 04035 | 1 . 04044 ± 0 . 00068         |
| τ                | 0 . 095         | 0 . 096 +0 . 013 - 0 . 016       | 0 . 087             | 0 . 095 +0 . 013 - 0 . 016       | 0 . 099   | 0 . 100 +0 . 013 - 0 . 016    |
| m eff ν, sterile | 0 . 005         | < 0 . 701                        | 0 . 67              | 0 . 70 ± 0 . 16                  | 0 . 62    | 0 . 68 ± 0 . 15               |
| w                | - 1 . 030       | - 1 . 060 +0 . 160 - 0 . 110     | - 1 . 210           | - 1 . 200 ± 0 . 100              | - 1 . 167 | - 1 . 170 +0 . 105 - 0 . 096  |
| N eff            | 3 . 44          | 3 . 72 +0 . 31 - 0 . 45          | 3 . 33              | 3 . 56 +0 . 15 - 0 . 33          | 3 . 80    | 3 . 81 +0 . 28 - 0 . 35       |
| n s              | 0 . 975         | 0 . 984 +0 . 014 - 0 . 017       | 0 . 962             | 0 . 974 +0 . 011 - 0 . 017       | 0 . 986   | 0 . 987 +0 . 014 - 0 . 013    |
| ln(10 10 A s )   | 3 . 111         | 3 . 119 +0 . 03 - 0 . 034        | 3 . 087             | 3 . 110 +0 . 028 - 0 . 034       | 3 . 122   | 3 . 126 +0 . 029 - 0 . 034    |
| r 0 . 05         | 0 . 007         | < 0 . 188                        | 0 . 006             | < 0 . 179                        | 0 . 203   | 0 . 186 +0 . 034 - 0 . 041    |
| Ω Λ              | 0 . 709         | 0 . 702 +0 . 017 - 0 . 020       | 0 . 720             | 0 . 717 +0 . 014 - 0 . 013       | 0 . 716   | 0 . 713 +0 . 014 - 0 . 013    |
| Ω m              | 0 . 291         | 0 . 298 +0 . 020 - 0 . 017       | 0 . 280             | 0 . 283 +0 . 013 - 0 . 014       | 0 . 284   | 0 . 287 +0 . 013 - 0 . 014    |
| σ 8              | 0 . 851         | 0 . 816 +0 . 038 - 0 . 041       | 0 . 773             | 0 . 773 ± 0 . 014                | 0 . 776   | 0 . 772 +0 . 014 - 0 . 015    |
| H 0              | 71 . 2          | 71 . 9 +2 . 3 - 3 . 2            | 73 . 0              | 73 . 7 ± 2 . 0                   | 74 . 3    | 74 . 1 ± 2 . 0                |
| r 0 . 002        | 0 . 006         | < 0 . 199                        | 0 . 006             | < 0 . 183                        | 0 . 210   | 0 . 193 +0 . 038 - 0 . 050    |
| - ln L max       | 4904.85         | 4904.85                          | 4909.49             | 4909.49                          | 4932.30   | 4932.30                       |

Figure 7 . Cosmological constraints on the w CDM+ + r ∑ m ν + N eff model.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0481v3-figure-2096792e-9a8b-41c1-84b3-ea9b7263acf0.jpg)

- [12] C. Giunti, M. Laveder, Y. F. Li, Q. Y. Liu and H. W. Long, 'Update of Short-Baseline Electron Neutrino and Antineutrino Disappearance,' Phys. Rev. D 86 , 113014 (2012) [arXiv:1210.5715 [hep-ph]].
- [13] C. Giunti, M. Laveder, Y. F. Li and H. W. Long, 'Pragmatic View of Short-Baseline Neutrino Oscillations,' Phys. Rev. D 88 , 073008 (2013) [arXiv:1308.5288 [hep-ph]].
- [14] M. Wyman, D. H. Rudd, R. A. Vanderveld and W. Hu, ' ν ΛCDM: Neutrinos help reconcile Planck with the Local Universe,' Phys. Rev. Lett. 112 , 051302 (2014) [arXiv:1307.7715 [astro-ph.CO]].
- [15] J. Hamann and J. Hasenkamp, 'A new life for sterile neutrinos: resolving inconsistencies using hot dark matter,' JCAP 1310 , 044 (2013) [arXiv:1308.3255 [astro-ph.CO]].

Figure 8 . Cosmological constraints on the w CDM+ + r N eff + m eff ν, sterile model.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0481v3-figure-fa26a920-1b1d-41aa-bbc8-33a35a22b958.jpg)

- [16] P. A. R. Ade et al. [BICEP2 Collaboration], 'Detection of B-Mode Polarization at Degree Angular Scales by BICEP2,' Phys. Rev. Lett. 112 , 241101 (2014) [arXiv:1403.3985 [astro-ph.CO]].
- [17] H. Liu, P. Mertsch and S. Sarkar, 'Fingerprints of Galactic Loop I on the Cosmic Microwave Background,' arXiv:1404.1899 [astro-ph.CO].
- [18] K. Harigaya and T. T. Yanagida, 'Discovery of Large Scale Tensor Mode and Chaotic Inflation in Supergravity,' arXiv:1403.4729 [hep-ph].
- [19] K. Nakayama and F. Takahashi, 'Higgs Chaotic Inflation and the Primordial B-mode Polarization Discovered by BICEP2,' arXiv:1403.4132 [hep-ph].

| [20]   | R. H. Brandenberger, A. Nayeri and S. P. Patil, 'Closed String Thermodynamics and a Blue Tensor Spectrum,' arXiv:1403.4927 [astro-ph.CO].                                                                     |
|--------|---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| [21]   | C. R. Contaldi, M. Peloso and L. Sorbo, 'Suppressing the impact of a high tensor-to-scalar ratio on the temperature anisotropies,' arXiv:1403.4596 [astro-ph.CO].                                             |
| [22]   | V. c. Miranda, W. Hu and P. Adshead, 'Steps to Reconcile Inflationary Tensor and Scalar Spectra,' arXiv:1403.5231 [astro-ph.CO].                                                                              |
| [23]   | M. Gerbino, A. Marchini, L. Pagano, L. Salvati, E. Di Valentino and A. Melchiorri, 'Blue Gravity Waves from BICEP2 ?,' arXiv:1403.5732 [astro-ph.CO].                                                         |
| [24]   | J. McDonald, 'Negative Running of the Spectral Index, Hemispherical Asymmetry and Consistency of Planck with BICEP2,' arXiv:1403.6650 [astro-ph.CO].                                                          |
| [25]   | D. K. Hazra, A. Shafieloo, G. F. Smoot and A. A. Starobinsky, 'Ruling out the power-law form of the scalar primordial spectrum,' arXiv:1403.7786 [astro-ph.CO].                                               |
| [26]   | D. K. Hazra, A. Shafieloo, G. F. Smoot and A. A. Starobinsky, 'Whipped inflation,' arXiv:1404.0360 [astro-ph.CO].                                                                                             |
| [27]   | A. Kehagias and A. Riotto, 'Remarks about the Tensor Mode Detection by the BICEP2 Collaboration and the Super-Planckian Excursions of the Inflaton Field,' arXiv:1403.4811 [astro-ph.CO].                     |
| [28]   | D. H. Lyth, 'BICEP2, the curvature perturbation and supersymmetry,' arXiv:1403.7323 [hep-ph].                                                                                                                 |
| [29]   | C. Bonvin, R. Durrer and R. Maartens, 'Can primordial magnetic fields be the origin of the BICEP2 data?,' arXiv:1403.6768 [astro-ph.CO].                                                                      |
| [30]   | J. Lizarraga, J. Urrestilla, D. Daverio, M. Hindmarsh, M. Kunz and A. R. Liddle, 'Can topological defects mimic the BICEP2 B-mode signal?,' arXiv:1403.4924 [astro-ph.CO].                                    |
| [31]   | A. Moss and L. Pogosian, 'Did BICEP2 see vector modes? First B-mode constraints on cosmic defects,' arXiv:1403.6105 [astro-ph.CO].                                                                            |
| [32]   | J. Chluba, L. Dai, D. Jeong, M. Kamionkowski and A. Yoho, 'Linking the BICEP2 result and the hemispherical power asymmetry through spatial variation of r ,' arXiv:1404.2798 [astro-ph.CO].                   |
| [33]   | J. -F. Zhang, Y. -H. Li and X. Zhang, 'Sterile neutrinos help reconcile the observational results of primordial gravitational waves from Planck and BICEP2,' arXiv:1403.7028 [astro-ph.CO].                   |
| [34]   | C. Dvorkin, M. Wyman, D. H. Rudd and W. Hu, 'Neutrinos help reconcile Planck measurements with both Early and Local Universe,' arXiv:1403.8049 [astro-ph.CO].                                                 |
| [35]   | J. -F. Zhang, Y. -H. Li and X. Zhang, 'Cosmological constraints on neutrinos after BICEP2,' Eur. Phys. J. C 74 , 2954 (2014) [arXiv:1404.3598 [astro-ph.CO]].                                                 |
| [36]   | Y. -H. Li, J. -F. Zhang and X. Zhang, 'Tilt of primordial gravitational wave spectrum in a universe with sterile neutrinos,' Sci. China Phys. Mech. Astron. 57 , 1455 (2014) [arXiv:1405.0570 [astro-ph.CO]]. |
| [37]   | M. Archidiacono, N. Fornengo, S. Gariazzo, C. Giunti, S. Hannestad and M. Laveder, 'Light sterile neutrinos after BICEP-2,' JCAP 1406 , 031 (2014) [arXiv:1404.1794 [astro-ph.CO]].                           |
| [38]   | J. Bergstrm, M. C. Gonzalez-Garcia, V. Niro and J. Salvado, 'Statistical tests of sterile neutrinos using cosmology and short-baseline data,' arXiv:1407.3806 [hep-ph].                                       |
| [39]   | M. Li, X. -D. Li, Y. -Z. Ma, X. Zhang and Z. Zhang, 'Planck Constraints on Holographic Dark Energy,' JCAP 1309 , 021 (2013) [arXiv:1305.5302 [astro-ph.CO]].                                                  |
| [40]   | A. Lewis and S. Bridle, 'Cosmological parameters from CMB and other data: A Monte Carlo approach,' Phys. Rev. D 66 , 103511 (2002) [astro-ph/0205436].                                                        |

| [41]   | L. Anderson et al. [BOSS Collaboration], 'The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: Baryon Acoustic Oscillations in the Data Release 10 and 11 galaxy samples,' arXiv:1312.4877 [astro-ph.CO].                                                           |
|--------|----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| [42]   | F. Beutler, C. Blake, M. Colless, D. H. Jones, L. Staveley-Smith, L. Campbell, Q. Parker and W. Saunders et al. , 'The 6dF Galaxy Survey: Baryon Acoustic Oscillations and the Local Hubble Constant,' Mon. Not. Roy. Astron. Soc. 416 , 3017 (2011) [arXiv:1106.3366 [astro-ph.CO]].              |
| [43]   | W. J. Percival et al. [SDSS Collaboration], 'Baryon Acoustic Oscillations in the Sloan Digital Sky Survey Data Release 7 Galaxy Sample,' Mon. Not. Roy. Astron. Soc. 401 , 2148 (2010) [arXiv:0907.1660 [astro-ph.CO]].                                                                            |
| [44]   | C. Blake, E. Kazin, F. Beutler, T. Davis, D. Parkinson, S. Brough, M. Colless and C. Contreras et al. , 'The WiggleZ Dark Energy Survey: mapping the distance-redshift relation with baryon acoustic oscillations,' Mon. Not. Roy. Astron. Soc. 418 , 1707 (2011) [arXiv:1108.2635 [astro-ph.CO]]. |
| [45]   | J. L. Tinker, A. V. Kravtsov, A. Klypin, K. Abazajian, M. S. Warren, G. Yepes, S. Gottlober and D. E. Holz, 'Toward a halo mass function for precision cosmology: The Limits of universality,' Astrophys. J. 688 , 709 (2008) [arXiv:0803.2706 [astro-ph]].                                        |
| [46]   | W. A. Watson, I. T. Iliev, A. D'Aloisio, A. Knebe, P. R. Shapiro and G. Yepes, 'The halo mass function through the cosmic ages,' Mon. Not. Roy. Astron. Soc. 433 , 1230 (2013) [arXiv:1212.0095 [astro-ph.CO]].                                                                                    |
| [47]   | P. A. R. Ade et al. [Planck Collaboration], 'Planck 2013 results. XVII. Gravitational lensing by large-scale structure,' arXiv:1303.5077 [astro-ph.CO].                                                                                                                                            |
Jing-Fei Zhang Jia-Jia Geng Xin Zhang
Neutrinos and dark energy after Planck and BICEP2: data consistency tests and cosmological parameter constraints
The detection of the B-mode polarization of the cosmic microwave background
(CMB) by the BICEP2 experiment implies that the tensor-to-scalar ratio $r$
should be involved in the base standard cosmology. In this paper, we extend the
$\Lambda$CDM+$r$+neutrino/dark radiation models by replacing the cosmological
constant with the dynamical dark energy with constant $w$. Four neutrino plus
dark energy models are considered, i.e., the $w$CDM+$r+\sum m_\nu$, $w$CDM+r +
$N_{\rm eff}$, $w$CDM+r + $\sum m_\nu$ + $N_{\rm eff}$, and $w$CDM+r + $N_{\rm
eff}$ + $m_{\nu,{\rm sterile}}^{\rm eff}$ models. The current observational
data considered in this paper include the Planck temperature data, the WMAP
9-year polarization data, the baryon acoustic oscillation data, the Hubble
constant direct measurement data, the Planck Sunyaev-Zeldovich cluster counts
data, the Planck CMB lensing data, the cosmic shear data, and the BICEP2
polarization data. We test the data consistency in the four cosmological
models, and then combine the consistent data sets to perform joint constraints
on the models. We focus on the constraints on the parameters $w$, $\sum m_\nu$,
$N_{\rm eff}$, and $m_{\nu,{\rm sterile}}^{\rm eff}$.