Constraints on inflation revisited: An analysis including the latest local measurement of the Hubble constant

TL;DR

This paper revisits inflation constraints using the latest local Hubble constant measurement, combining multiple cosmological data sets, and finds that certain inflation models are favored or excluded at specific confidence levels.

Contribution

It incorporates the latest local Hubble constant measurement and additional parameters to refine inflation model constraints using comprehensive cosmological data.

Findings

Scale invariance is excluded at 3.3σ level.

Positive ΔN_eff is favored at 1.6σ.

Spontaneously broken SUSY inflation is most favored.

Abstract

We revisit the constraints on inflation models by using the current cosmological observations involving the latest local measurement of the Hubble constant ($H_{0} = 73.00\pm 1.75$ km s $^{-1}$ Mpc$^{-1}$). We constrain the primordial power spectra of both scalar and tensor perturbations with the observational data including the Planck 2015 CMB full data, the BICEP2 and Keck Array CMB B-mode data, the BAO data, and the direct measurement of $H_0$. In order to relieve the tension between the local determination of the Hubble constant and the other astrophysical observations, we consider the additional parameter $N_{\rm eff}$ in the cosmological model. We find that, for the $\Lambda$CDM+$r$+$N_{\rm eff}$ model, the scale invariance is only excluded at the 3.3$\sigma$ level, and $\Delta N_{\rm eff}>0$ is favored at the 1.6$\sigma$ level. Comparing the obtained 1$\sigma$ and 2$\sigma$ contours of $(n_s,r)$ with the theoretical predictions of selected inflation models, we find that both the convex and concave potentials are favored at 2$\sigma$ level, the natural inflation model is excluded at more than 2$\sigma$ level, the Starobinsky $R^2$ inflation model is only favored at around 2$\sigma$ level, and the spontaneously broken SUSY inflation model is now the most favored model.

Figures

Tables

Table 1. Table 1: The fitting results of the cosmological parameters in the Λ Λ \Lambda CDM+ r 𝑟 r , Λ Λ \Lambda CDM+ r 𝑟 r + N eff subscript 𝑁 eff N_{\rm eff} , Λ Λ \Lambda CDM+ r 𝑟 r + d ​ n s / d ​ ln ⁡ k 𝑑 subscript 𝑛 s 𝑑 𝑘 dn_{\rm s}/d\ln k , and Λ Λ \Lambda CDM+ r 𝑟 r + d ​ n s / d ​ ln ⁡ k 𝑑 subscript 𝑛 s 𝑑 𝑘 dn_{\rm s}/d\ln k + N eff subscript 𝑁 eff N_{\rm eff} models using the Planck+BK+BAO+ H 0 subscript 𝐻 0 H_{0} data.

Parameter	$\Lambda$CDM+$r$	$\Lambda$CDM+$r$+$N_{\mathrm{eff}}$	$\Lambda$CDM+$r$+$d{n}_s/d\ln k$	$\Lambda$CDM+$r$+$d{n}_s/d\ln k$+$N_{\mathrm{eff}}$
${\Omega }_b{h}^2$	$0.02238\pm 0.00014$	$0.02253\pm 0.00017$	$0.02242\pm 0.00015$	$0.02253\pm 0.00019$
${\Omega }_c{h}^2$	$0.1177\pm 0.0010$	$0.1216\pm 0.0027$	$0.1174\pm 0.0010$	$0.1218\pm 0.0030$
$100{\theta }_{\mathrm{MC}}$	$1.04106\pm 0.00029$	${1.04062}_{-0.00038}^{+0.00039}$	${1.04109}_{-0.00030}^{+0.00029}$	$1.04061\pm 0.00042$
$\tau$	$0.075\pm 0.012$	$0.074\pm 0.012$	$0.077\pm 0.012$	$0.074\pm 0.014$
$\ln \left({10}^{10}{A}_s\right)$	$3.079\pm 0.023$	$3.088\pm 0.023$	$3.084\pm 0.023$	$3.087\pm 0.027$
$n_s$	${0.9699}_{-0.0039}^{+0.0040}$	${0.9787}_{-0.0065}^{+0.0064}$	${0.9701}_{-0.0042}^{+0.0041}$	$0.9781\pm 0.0080$
$d{n}_s/d\ln k$	…	…	$-{0.0042}_{-0.0066}^{+0.0067}$	${0.0010}_{-0.0073}^{+0.0074}$
$r_{0.002}$ ($2\sigma$)
$N_{\mathrm{eff}}$	…	$3.30\pm 0.16$	…	$3.30\pm 0.18$
${\Omega }_m$	$0.3023\pm 0.0060$	$0.2988\pm 0.0060$	${0.3007}_{-0.0061}^{+0.0060}$	${0.2988}_{-0.0094}^{+0.0085}$
$H_0$	${68.23}_{-0.46}^{+0.47}$	$69.63\pm 0.99$	${68.37}_{-0.50}^{+0.47}$	${69.70}_{-1.30}^{+1.20}$
${\sigma }_8$	${0.8188}_{-0.0085}^{+0.0087}$	$0.8300\pm 0.0110$	${0.8193}_{-0.0086}^{+0.0085}$	$0.8310\pm 0.0120$
${\chi }_{\min }^2$	$13616.988$	$13612.184$	$13615.324$	$13611.122$

Equations

$$H^2=\frac{1}{3M_{\mathrm{pl}}^2}\left[\frac{1}{2}{\dot{\varphi }}^2+V(\varphi )\right],$$

$$\ddot{\varphi }+3H\dot{\varphi }+V^{\prime }(\varphi )=0,$$

$$H^2\approx \frac{V(\varphi )}{3M_{\mathrm{pl}}^2},$$

$$3H\dot{\varphi }\approx -V^{\prime }(\varphi ).$$

$$ϵ=\frac{M_{\mathrm{pl}}^2}{2}{\left[\frac{V^{\prime }(\varphi )}{V(\varphi )}\right]}^2,$$

$$\eta =M_{\mathrm{pl}}^2\left[\frac{V^{\prime \prime }(\varphi )}{V(\varphi )}\right],$$

$${\xi }^2=\frac{M_{\mathrm{pl}}^4{V}^{\prime }(\varphi )V^{\prime \prime \prime }(\varphi )}{V^2(\varphi )},$$

$$r=\frac{P_{\mathrm{t}}(k)}{P_{\mathrm{s}}(k)}=16ϵ.$$

$$n_{\mathrm{s}}=1-6ϵ+2\eta ,$$

$$d{n}_{\mathrm{s}}/d\ln k=16ϵ\eta -24{ϵ}^2-2{\xi }^2.$$

$$P_{\mathrm{s}}(k)=A_{\mathrm{s}}{\left(\frac{k}{k_{\ast }}\right)}^{n_{\mathrm{s}}-1+\frac{1}{2}\frac{d{n}_{\mathrm{s}}}{d\ln k}\ln \left(\frac{k}{k_{\ast }}\right)},$$

$$P_{\mathrm{t}}(k)=A_{\mathrm{t}}{\left(\frac{k}{k_{\ast }}\right)}^{n_{\mathrm{t}}+\frac{1}{2}\frac{d{n}_{\mathrm{t}}}{d\ln k}\ln \left(\frac{k}{k_{\ast }}\right)},$$

$$\mathbf{P}=\{{\Omega }_{\mathrm{b}}{h}^2,{\Omega }_{\mathrm{c}}{h}^2,100{\theta }_{\mathrm{MC}},\tau ,\ln (10^{10}{A}_{\mathrm{s}}),n_{\mathrm{s}},r\},$$

$$\small\left.\begin{array}[]{c}r_{0.002}<0.069~{}(2\sigma)\\ \\ n_{\rm s}=0.9699^{+0.0040}_{-0.0039}~{}(1\sigma)\end{array}\right\}\quad\mbox{$\Lambda$CDM+$r$}.$$

$${\rho }_r=\left[1+N_{\mathrm{eff}}\frac{7}{8}{\left(\frac{4}{11}\right)}^{4/3}\right]{\rho }_{\gamma },$$

$$\small\left.\begin{array}[]{c}r_{0.002}<0.071~{}(2\sigma)\\ \\ n_{\rm s}=0.9787^{+0.0064}_{-0.0065}~{}(1\sigma)\end{array}\right\}\quad\mbox{$\Lambda$CDM+$r$+$N_{\rm eff}$}.$$

$$\small\left.\begin{array}[]{c}r_{0.002}<0.077~{}(2\sigma)\\ \\ n_{\rm s}=0.9701^{+0.0041}_{-0.0042}~{}(1\sigma)\\ \\ dn_{\rm s}/d\ln k=-0.0042^{+0.0067}_{-0.0066}~{}(1\sigma)\end{array}\right\}\quad\mbox{$\Lambda$CDM+$r$+$dn_{\rm s}/d\ln k$}.$$

$$\small\left.\begin{array}[]{c}r_{0.002}<0.074~{}(2\sigma)\\ \\ n_{\rm s}=0.9781\pm 0.0080~{}(1\sigma)\\ \\ dn_{\rm s}/d\ln k=0.0010^{+0.0074}_{-0.0073}~{}(1\sigma)\end{array}\right\}\quad\mbox{$\Lambda$CDM+$r$+$dn_{\rm s}/d\ln k$+$N_{\rm eff}$}.$$

$$r=\frac{4n}{N},$$

$$n_{\mathrm{s}}=1-\frac{n+2}{2N},$$

$$r=\frac{8}{(f/M_{\mathrm{pl}}{)}^2}\frac{1+\cos {\theta }_N}{1-\cos {\theta }_N},$$

$$n_{\mathrm{s}}=1-\frac{1}{(f/M_{\mathrm{pl}}{)}^2}\frac{3+\cos {\theta }_N}{1-\cos {\theta }_N},$$

$$\cos \frac{{\theta }_N}{2}=\exp \left(-\frac{N}{2(f/M_{\mathrm{pl}}{)}^2}\right).$$

$$r\simeq 0,$$

$$n_{\mathrm{s}}=1-\frac{1}{N}.$$

$$r\simeq \frac{12}{N^2},$$

$$n_{\mathrm{s}}=1-\frac{2}{N}.$$

Full text

I Introduction

Inflation is the leading paradigm to explain the origin of the primordial density perturbations and the primordial gravitational waves, which is a period of accelerated expansion of the early universe. It can resolve a number of puzzles of the standard cosmology, such as the horizon, flatness, and monopole problems [1, 2, 3, 4], and offer the initial conditions for the standard cosmology. During the epoch, inflation can generate the primordial density perturbations, which seeded the cosmic microwave background (CMB) anisotropies and the large-scale structure (LSS) formation in our universe. Thus, current cosmological observations can be used to explore the nature of inflation. For example, the measurements of CMB anisotropies have confirmed that inflation can provide a nearly scale-invariant primordial power spectrum [5, 6, 7, 8].

Although inflation took place at energy scale as high as $10^{16}$ GeV, where particle physics remains elusive, hundreds of different theoretical scenarios have been proposed. Thus selecting an actual version of inflation has become a major issue in the current study. As mentioned above, the primordial perturbations can lead to the CMB anisotropies and LSS formation, so comparing the predictions of these inflation models with cosmological data can provide the possibility to identify the suitable inflation models.

The astronomical observations measuring the CMB anisotropies have provided an excellent opportunity to explore the physics in the early universe. The Planck collaboration [9] has measured the primordial power spectrum of density perturbations with an unprecedented accuracy. Namely, the spectral index is measured to be $n_{s}=0.968\pm 0.006$ ($1\sigma$), ruling out the scale invariance at more than $5\sigma$, and the running of the spectral index is measured to be $dn_{\rm s}/d\ln k=-0.003\pm 0.007$ ($1\sigma$), from the Planck temperature data combined with the Planck lensing likelihood. The constraint on the tensor-to-scalar ratio is $r_{0.002}<0.11$ at the $2\sigma$ level, also derived by using the Planck temperature data combined with the Planck lensing likelihood. In addition, the Keck Array and BICEP2 collaborations [10] released a highly significant detection of B-mode polarization with inclusion of the first Keck Array B-mode polarization at 95 GHz. These data were taken by the BICEP2 and Keck Array CMB polarization experiments up to and including the 2014 observing season to improve the current constraints on primordial power spectra. The constraint on the tensor-to-scalar ratio is $r_{0.05}<0.09$ at the $2\sigma$ level from the B-mode only data of BICEP2 and Keck Array. The tighter constraint is $r_{0.05}<0.07$ at the $2\sigma$ level when the BICEP2/Keck Array B-mode data are combined with the Planck CMB data plus other astrophysical observations.

The baryon acoustic oscillation (BAO) data can effectively break the degeneracies between cosmological parameters and further improve the constraints on inflation models (see, e.g., Refs. [11, 12, 13, 14, 15]). In this paper, we employ the latest BAO measurements including the Date Release 12 of the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS DR12) [16], the 6dF Galaxy Survey (6dFGS) measurement [17], and the Main Galaxy Sample of Data Release 7 of Sloan Digital Sky Survey (SDSS-MGS) [18].

Recently, Riess et al. [19] reported their new result of direct measurement of the Hubble constant, $H_{0}=73.00\pm 1.75$ km s -1 Mpc*-1*, which is $3.3\sigma$ higher than the fitting result, $H_{0}=66.93\pm 0.62$ km s -1 Mpc*-1*, derived by the Planck collaboration [20] based on the $\Lambda$CDM model assuming $\sum m_{\nu}=0.06$ eV using the Planck TT, TE, EE+lowP data. The strong tension between the new measurement of $H_{0}$ and the Planck data may be from some systematic uncertainties in the measurements or some new physics effects. In order to reconcile the new measurement of $H_{0}$ and the Planck data, one can consider the new physics by adding some extra parameters, such as the parameters describing a dynamical dark energy [21, 22], extra relativistic degrees of freedom [19, 23, 24, 25, 26] and light sterile neutrinos [23, 24, 27, 28, 29, 30, 31].

Although there are strong tensions between the new measurement of $H_{0}$ and other cosmological observations, the result of $H_{0}=73.00\pm 1.75$ km s -1 Mpc*-1* can play an important role in current cosmology due to its reduced uncertainty from 3.3% to 2.4%. In this paper, we combine the new measurement of $H_{0}$ with the Planck data, the BICEP2/Keck Array data and the BAO data to constrain inflation models. The aim of this work is to investigate whether the local determination $H_{0}=73.00\pm 1.75$ km s -1 Mpc*-1* will have a remarkable influence on constraining the primordial power spectra of scalar and tensor perturbations. In order to relieve the tension between the local determination of the Hubble constant and other astrophysical observations, we decide to consider dark radiation, parametrized by $\Delta N_{\rm eff}$ (defined by $N_{\rm eff}-3.046$), in the cosmological model in our analysis. The constraint results of $(n_{s},r)$ will be compared with the theoretical predictions of some typical inflation models to make a model selection analysis.

The structure of the paper is organized as follows. In Sec. II, we briefly introduce the single-field slow-roll inflationary scenario. In Sec. III, we report the results of the constraints on the primordial power spectra with the combination of the Planck data, the BICEP2/Keck Array data, the BAO data and the latest measurement of $H_{0}$. In Sec. IV, we compare the constraint results of $(n_{s},r)$ with the theoretical predictions of some typical inflationary models and show the impacts of the latest measurement of $H_{0}$ on the inflation model selection. Conclusion is given in Sec. V.

II Slow-roll inflationary scenario

In this paper, we only consider the simplest inflationary scenario within the slow-roll paradigm, for which the accelerated expansion of early universe is driven by a homogeneous, slowly rolling scalar field $\phi$. According to the energy density of the inflaton $\rho_{\phi}=\dot{\phi}^{2}/{2}+V(\phi)$, the Friedmann equation becomes

[TABLE]

where $H=\dot{a}/a$ (with $a$ the scale factor of the universe) is the Hubble parameter, $M_{\rm pl}=1/\sqrt{8\pi G}$ is the reduced Planck mass, $V(\phi)$ is the inflaton potential, and the dot denotes the derivative with respect to the cosmic time $t$.

The equation of motion for the inflaton satisfies

[TABLE]

where the prime is the derivative with respect to the inflaton $\phi$. Due to the slow-roll approximation, $\dot{\phi}^{2}\ll 0$ and $\ddot{\phi}\ll 0$, Eqs. (1) and (2) can be reduced to

[TABLE]

[TABLE]

Usually, the inflationary universe can be characterized with the slow-roll parameters, which can be defined as

[TABLE]

[TABLE]

[TABLE]

and so on. The inflaton slowly rolls down its potential $V(\phi)$ as long as $\epsilon\ll 1$ and $|\eta|\ll 1$.

The tensor-to-scalar ratio, which is defined to be the ratio of the tensor spectrum $P_{\rm t}(k)$ to the scalar spectrum $P_{\rm s}(k)$, can be given by the slow-roll approximation as

[TABLE]

Similarly, according to the slow-roll approximation, we can obtain the spectral index

[TABLE]

and the running spectral index

[TABLE]

By constraining these parameters using cosmological observations, we can effectively distinguish between different inflation models.

III Constraints on primordial power spectra

In this section, we make a comprehensive analysis of constraining the primordial power spectra of scalar and tensor perturbations by combining the new measurement of the Hubble constant, $H_{0}=73.00\pm 1.75$ km s -1 Mpc*-1* [19], with the Planck data, the BICEP2/Keck Array data and the BAO data, to investigate how the new measurement of $H_{0}$ affects the constraint results of inflation models. We employ the Planck CMB 2015 data set including the temperature power spectrum (TT), the polarization power spectrum (EE), the cross-correlation power spectrum of temperature and polarization (TE), and the Planck low-$\ell$ ($\ell\leq 30$) likelihood (lowP), as well as the lensing reconstruction, which is abbreviated as “Planck”. We employ all the BICEP2 and Keck Array B-mode data with inclusion of 95 GHz band, abbreviated as “BK”. The BAO data include the CMASS and LOWZ samples from the BOSS DR12 at $z_{\rm eff}=0.57$ and $z_{\rm eff}=0.32$ [16], the 6dFGS measurement at $z_{\rm eff}=0.106$ [17], and the SDSS-MGS measurement at $z_{\rm eff}=0.15$ [18], abbreviated as “BAO”.

The primordial power spectra of scalar and tensor perturbations can be expressed as

[TABLE]

[TABLE]

where $A_{\rm s}$ and $A_{\rm t}$ correspond to the scalar and tensor amplitudes at the pivot scale $k_{*}$, respectively. For the canonical single-field slow-roll inflation model without the inclusion of the running of the spectral index, we have the consistency relation $n_{\rm t}=-r/8$. When the running spectral index is considered, we then have $n_{\rm t}=-r(2-r/8-n_{\rm s})/8$ and $dn_{\rm t}/d\ln k=r(r/8+n_{\rm s}-1)/8$. We uniformly set the pivot scale as $k_{*}=0.002$ Mpc*-1* in this work.

There are seven independent free parameters in the base $\Lambda$CDM+$r$ model:

[TABLE]

where $\Omega_{\rm b}h^{2}$ and $\Omega_{\rm c}h^{2}$ denote the present-day densities of baryon and cold dark matter; $\theta_{\rm MC}$ denotes the ratio of the sound horizon $r_{\rm s}$ to the angular diameter distance $D_{\rm A}$ at the last-scattering epoch; $\tau$ denotes the optical depth to reionization; $A_{\rm s}$ and $n_{\rm s}$ denote the amplitude and the spectral index of the primordial power spectra of scalar perturbations, respectively; $r$ denotes the the tensor-to-scalar ratio. When the running is considered, the parameter $dn_{\rm s}/d\ln k$ is added to the cosmological model. In this work, we derive the posterior parameter probabilities by using the Markov Chain Monte Carlo (MCMC) sampler CosmoMC [32].

In Fig. 1, we give one-dimensional marginalized distributions and two-dimensional contours ($1\sigma$ and $2\sigma$) for the parameters $n_{\rm s}$, $r_{\rm 0.002}$ and $H_{0}$ in the $\Lambda$CDM+$r$ model using the Planck+BK+BAO+$H_{0}$ data. The constraint results of the $\Lambda$CDM+$r$ model are summarized in the second column of Table 1. Here we quote $\pm 1\sigma$ limits for every parameter in the $\Lambda$CDM+$r$ model, except for $r$, which is quoted with the $2\sigma$ upper limit. We obtain the constraints on $r$ and $n_{\rm s}$:

[TABLE]

The result of $n_{\rm s}$ for the primordial power spectrum of scalar perturbations excludes the Harrison-Zel’dovich (HZ) scale-invariant spectrum with $n_{\rm s}=1$ at the $7.5\sigma$ level.

In addition, the constraint on the Hubble constant is $H_{0}=68.23^{+0.47}_{-0.46}$ km s*-1* Mpc*-1*, which is $2.6\sigma$ less than the local determination $H_{0}=73.00\pm 1.75$ km s*-1* Mpc*-1*. Namely, the direct measurement of $H_{0}=73.00\pm 1.75$ km s*-1* Mpc*-1* is in tension with the fit result derived by the Planck+BK+BAO+$H_{0}$ data based on the $\Lambda$CDM+$r$ model. As shown in Fig. 2, the green line denotes the one-dimensional posterior distribution for the parameter $H_{0}$ in the $\Lambda$CDM+$r$ model using the Planck+BK+BAO+$H_{0}$ data, and the light red band denotes the new local measurement of $H_{0}$. Obviously, there is a strong tension between the two results.

Next, we consider the extra relativistic degrees of freedom (i.e., the additional parameter $N_{\rm eff}$) in the cosmological model to relieve the tension between the latest measurement of $H_{0}$ and other observational data. The total radiation energy density in the universe is given by

[TABLE]

where $\rho_{\gamma}$ is the energy density of photons. If there are only three-species active neutrinos in the universe, we have the standard value of $N_{\rm eff}=3.046$. Any additional value of $\Delta N_{\rm eff}=N_{\rm eff}-3.046>0$ indicates the existence of some dark radiation in the universe. Now, we follow Planck collaboration [9] to constrain $N_{\rm eff}$ as a free parameter, varying within its prior range of $[0,6]$. Values of $N_{\rm eff}<3.046$ are less well motivated, because such values would require that standard neutrinos are incompletely thermalized or additional photons are produced after the neutrino decoupling, but we still include this range for completeness.

The third column of Table 1 gives the constraint results of the cosmological parameters in the $\Lambda$CDM+$r$+$N_{\rm eff}$ model using the Planck+BK+BAO+$H_{0}$ data. We obtain the constraints on $r$ and $n_{\rm s}$:

[TABLE]

The value of $n_{\rm s}$ becomes larger than that without considering $N_{\rm eff}$. The fit result of $N_{\rm eff}=3.30\pm 0.16$ indicates that $\Delta N_{\rm eff}>0$ is favored at the $1.6\sigma$ level. Due to a positive correlation between $n_{\rm s}$ and $N_{\rm eff}$, as shown in Fig. 3, $\Delta N_{\rm eff}>0$ will lead to a larger $n_{\rm s}$.

On the other hand, a larger Hubble constant, $H_{0}=69.63\pm 0.99$ km s*-1* Mpc*-1*, is obtained when the parameter $N_{\rm eff}$ is considered, which is only $1.7\sigma$ less than the local determination $H_{0}=73.00\pm 1.75$ km s*-1* Mpc*-1*. Namely, the tension between $H_{0}=73.00\pm 1.75$ km s*-1* Mpc*-1* and other observational data is greatly alleviated by introducing the parameter $N_{\rm eff}$ in the cosmological model. As showed in Fig. 2, the constraint on $H_{0}$ derived using the Planck+BK+BAO+$H_{0}$ data in the $\Lambda$CDM+$r$+$N_{\rm eff}$ model is much closer to the local measurement of $H_{0}$. In addition, when the free parameter $N_{\rm eff}$ is included in the cosmological model, $\chi^{2}$ decreases from $13616.988$ to $13612.184$. The big $\chi^{2}$ difference, $\Delta\chi^{2}=-4.804$, implies that the $\Lambda$CDM+$r$+$N_{\rm eff}$ model, compared to the $\Lambda$CDM+$r$ model, is more favored by the current Planck+BK+BAO+$H_{0}$ data. Here we note that in this paper we compare models through only a $\chi^{2}_{\rm min}$ comparison, because we constrain these models using the same data combination. In this situation, if one additional parameter can lead to $\chi^{2}_{\rm min}$ decreasing by more than 2, then we say that adding this parameter is reasonable statistically. Thus, we do not employ Bayesian information criterion or Bayesian evidence in this paper, since a $\chi^{2}_{\rm min}$ comparison is sufficient for our task.

Furthermore, we consider the inclusion of the running of the spectral index, $dn_{\rm s}/d\ln k$, in the fit to the Planck+BK+BAO+$H_{0}$ data. Figure 4 gives one-dimensional marginalized distributions and two-dimensional contours ($1\sigma$ and $2\sigma$) for parameters $n_{\rm s}$, $dn_{\rm s}/d\ln k$, $r_{\rm 0.002}$, and $H_{0}$ in the $\Lambda$CDM+$r$+$dn_{\rm s}/d\ln k$ model using the Planck+BK+BAO+$H_{0}$ data. We obtain the constraints on $r$, $n_{\rm s}$ and $dn_{\rm s}/d\ln k$ (see also the fourth column in Table 1):

[TABLE]

We find that $dn_{\rm s}/d\ln k=0$ is well consistent with the Planck+BK+BAO+$H_{0}$ data in this case, and the fit result $H_{0}=68.37^{+0.47}_{-0.50}$ km s*-1* Mpc*-1* is still in tension with the direct $H_{0}$ measurement. The comparison with the $\Lambda$CDM+$r$ model gives $\Delta\chi^{2}=-1.664$, implying that adding the parameter $dn_{\rm s}/d\ln k$ does not effectively improve the fit. The comparison with the $\Lambda$CDM+$r$+$N_{\rm eff}$ model gives $\Delta\chi^{2}=3.14$, explicitly showing that $N_{\rm eff}$ is much more worthy to be added than $dn_{\rm s}/d\ln k$ in the sense of improving the fit.

In Fig. 5, we give one-dimensional marginalized distributions and two-dimensional contours ($1\sigma$ and $2\sigma$) for the parameters $N_{\rm eff}$, $n_{\rm s}$, $dn_{\rm s}/d\ln k$, $r_{\rm 0.002}$, and $H_{0}$ in the $\Lambda$CDM+$r$+$dn_{\rm s}/d\ln k$+$N_{\rm eff}$ model using the Planck+BK+BAO+$H_{0}$ data. We obtain the constraints on $r$, $n_{\rm s}$ and $dn_{\rm s}/d\ln k$ (see also the last column in Table 1):

[TABLE]

We find that the fitting results are almost unchanged comparing to the case of the $\Lambda$CDM+$r$+$N_{\rm eff}$ model (although the parameter space is slightly amplified), as shown in the third and fifth columns of Table 1. The results explicitly show that $dn_{\rm s}/d\ln k=0$ is in good agreement with the current observations. A $\chi^{2}$ comparison shows that, when the additional parameter $dn_{\rm s}/d\ln k$ is included, the $\chi^{2}_{\rm min}$ value decreases only by 1.062 (i.e., $\Delta\chi^{2}=-1.062$), which implies that the running of the spectral index $dn_{\rm s}/d\ln k$ is not deserved to be considered in the cosmological model in the sense of statistical significance.

IV Inflation model selection

In this section, we consider a few simple and representative inflation models and compare them with the constraint results given in the former section. See also Ref. [33] for a preliminary research. In what follows, we give the predictions of these inflation models for $r$ and $n_{\rm s}$. For these inflation models, we uniformly take the number of $e$-folds $N\in[50,60]$. In principle, adding the parameter $N_{\rm eff}$ modifies the radiation density and thereby changes the post-inflationary expansion history, so that the $e$-folding number $N$ becomes dependent on the value of $N_{\rm eff}$. However, practically it is hard to link $N$ to the actual observations. Thus, the usual treatment of considering $N\in[50,60]$ is of course applicable for our analysis.

The simplest class of inflation models has a monomial potential $V(\phi)\propto\phi^{n}$ [34], which is the prototype of the chaotic inflation model. They have the predictions:

[TABLE]

[TABLE]

where $n$ is any positive number. We take $n=2/3$, 1, and 2 as typical examples in this work. See also Refs. [35, 36, 37, 38] for relevant studies of this class of models.

The natural inflation model has the effective one-dimensional potential $V(\phi)=\Lambda^{4}(1+\cos(\phi/f))$ [39, 40], with the predictions:

[TABLE]

[TABLE]

where $\theta_{N}$ is given by

[TABLE]

Note that different values of $n_{\rm s}$ and $r$ result from the different decay constant $f$ when the number of $e$-folds $N$ is set to be a certain value.

The spontaneously broken SUSY (SBS) inflation model has the potential $V(\phi)=V_{0}(1+c\ln(\phi/Q))$ (where $V_{0}$ is dominant and the parameter $c\ll 1$) [41, 42, 43, 44, 45], with the predictions:

[TABLE]

[TABLE]

The Starobinsky $R^{2}$ inflation model is described by the action $S=\frac{M_{\rm pl}^{2}}{2}\int d^{4}x\sqrt{-g}(R+R^{2}/6M^{2})$ (where $M$ denotes an energy scale) [1], with the predictions:

[TABLE]

[TABLE]

In Fig. 6, we plot two-dimensional contours ($1\sigma$ and $2\sigma$) for $n_{\rm s}$ and $r_{\rm 0.002}$ using the Planck+BK+BAO and Planck+BK+BAO+$H_{0}$ data, compared to the theoretical predictions of selected inflation models. The orange contours denote the constraints on the $\Lambda$CDM+$r$ model with the Planck+BK+BAO data, the green contours denote the constraints on the $\Lambda$CDM+$r$ model with the Planck+BK+BAO+$H_{0}$ data, the gray contours denote the constraints on the $\Lambda$CDM+$r$+$N_{\rm eff}$ model with the Planck+BK+BAO data, and the blue contours denote the constraints on the $\Lambda$CDM+$r$+$N_{\rm eff}$ with the Planck+BK+BAO+$H_{0}$ data.

Comparing the orange and green contours, we find that when the direct measurement of $H_{0}$ is included in the data combination, the constraint on the $\Lambda$CDM+$r$ model is only changed a little, i.e., a little right shift of $n_{s}$ is yielded, which does not greatly change the result of inflation model selection (see also Ref. [11] for the case of orange contours). According to both the cases of orange and green contours, the inflation model with a convex potential is not favored; both the inflation model with a monomial potential ($\phi$ and $\phi^{2/3}$ cases) and the natural inflation model are marginally favored at around the 2$\sigma$ level; the SBS inflation model is located at out of the 2$\sigma$ region; the Starobinsky $R^{2}$ inflation model is the most favored model in this case.

When the parameter $N_{\rm eff}$ is considered in the analysis, and if the $H_{0}$ measurement is not used (i.e., using the Planck+BK+BAO data), we find that the parameter space is greatly amplified (mainly for $n_{s}$). Comparing the orange and gray contours, we find that without using the $H_{0}$ measurement the addition of $N_{\rm eff}$ can only amplify the range of $n_{s}$ but cannot lead to an obvious right shift of $n_{s}$.

When the $H_{0}$ measurement is also used, comparing the gray and blue contours, we see that the addition of the $H_{0}$ prior in the combination of data sets for constraining the $\Lambda$CDM+$r$+$N_{\rm eff}$ model leads to a considerable right shift of $n_{s}$ (and also a slight shrink of width for the range of $n_{s}$). In Fig. 3, we explicitly show that $H_{0}$ is positively correlated with $N_{\rm eff}$ and $N_{\rm eff}$ is positively correlated with $n_{s}$, which well explains why the $H_{0}$ prior (with a larger value of $H_{0}$) will lead to a larger value of $n_{s}$ in a cosmological model with $N_{\rm eff}$.

Next, we compare the green and blue contours, which is for the comparison of the $\Lambda$CDM+$r$ and $\Lambda$CDM+$r$+$N_{\rm eff}$ models with the Planck+BK+BAO+$H_{0}$ data, and we see that using the same data sets including the $H_{0}$ measurement, the consideration of $N_{\rm eff}$ yields a tremendous right shift of $n_{s}$ (see also Ref. [12]), which largely changes the result of the inflation model selection. As discussed in the last section, the $\Lambda$CDM+$r$+$N_{\rm eff}$ model is much better than the $\Lambda$CDM+$r$ model for the fit to the current Planck+BK+BAO+$H_{0}$ data, since the inclusion of $N_{\rm eff}$ makes the tension between $H_{0}$ measurement and other observations be greatly relieved and also leads to a much better fit (i.e., the $\chi^{2}_{\rm min}$ value is largely reduced).

We now compare the predictions of the above typical inflation models with the fit results of $(n_{s},r)$ corresponding to the blue contours. We see that, in this case, neither the concave potential nor the convex potential is excluded by the current data. But, it seems that, when comparing the two, the inflation model with the concave potential is more favored by the data. The natural inflation model is now excluded by the data at more than the 2$\sigma$ level. For the inflation models with a monomial potential, we find that the $\phi^{2}$ model is entirely excluded, the $\phi$ model is only marginally favored (at the edge of the 2$\sigma$ region), and the $\phi^{2/3}$ model is still well consistent with the current data (located in the 1$\sigma$ region). Now, the Starobinsky $R^{2}$ inflation model is not well favored, because it is located at the edge of the 2$\sigma$ region and actually the $N=50$ point even lies out of the 2$\sigma$ region. We find that in this case the most favored model is the SBS inflation model, which locates near the center of the contours.

Actually, the brane inflation model is also well consistent with the current data in this case (for previous analyses of brane inflation, see, e.g., Refs. [46, 47]). We leave a comprehensive analysis for the brane inflation model to a future work.

From the analysis in this paper, we have found that the inclusion of the latest local measurement of the Hubble constant can exert significant influence on the model selection of inflationary models, but one must be aware of that the result is dependent on the assumption of dark radiation in the cosmological model. Without the addition of the parameter $N_{\rm eff}$, the $H_{0}$ measurement is in tension with the Planck observation, and the $H_{0}$ prior actually does not greatly influence the fit result of the primordial power spectra (see the comparison of the orange and green contours in Fig. 6). The $H_{0}$ tension can be largely relieved provided that the parameter $N_{\rm eff}$ is considered in the model (the tension is reduced from 2.6$\sigma$ to 1.7$\sigma$). The inclusion of the $H_{0}$ measurement in the combination of data sets, together with the consideration of $N_{\rm eff}$ in the cosmological model, leads to a tremendous right shift of $n_{s}$ (see the comparison of the green and blue contours in Fig. 6), which greatly changes the situation of the inflation model selection. Future experiments on accurately measuring the Hubble constant and searching for light relics (dark radiation) would further test the robustness of our result in this paper.

V Conclusion

In this paper, we investigate how the constraints on the inflation models are affected by considering the latest local measurement of the Hubble constant in the cosmological global fit. We constrain the primordial power spectra of both scalar and tensor perturbations by using the current cosmological observations including the Planck 2015 CMB full data, the BICEP2 and Keck Array CMB B-mode data, the BAO data, and the direct measurement of $H_{0}$. In order to relieve the tension between the local determination of the Hubble constant and the other astrophysical observations, we consider the additional parameter $N_{\rm eff}$ in the cosmological model. We make comparison for the $\Lambda$CDM+$r$, $\Lambda$CDM+$r$+$N_{\rm eff}$, $\Lambda$CDM+$r$+$dn_{\rm s}/d\ln k$, and $\Lambda$CDM+$r$+$dn_{\rm s}/d\ln k$+$N_{\rm eff}$ models.

We find that the inclusion of $N_{\rm eff}$ indeed effectively relieves the tension. Comparing the $\Lambda$CDM+$r$ and $\Lambda$CDM+$r$+$N_{\rm eff}$ models, the tension is reduced from 2.6$\sigma$ to 1.7$\sigma$. The comparison also shows that the addition of one parameter, $N_{\rm eff}$, leads to the decrease of $\chi^{2}$ by 4.804. When the running of the spectral index $dn_{\rm s}/d\ln k$ is considered, we find that the fit results are basically not changed and $dn_{\rm s}/d\ln k=0$ is well consistent with the current data. Therefore, it is meaningful to consider the $\Lambda$CDM+$r$+$N_{\rm eff}$ model when the latest measurement of the Hubble constant is included in the analysis.

We constrain the $\Lambda$CDM+$r$+$N_{\rm eff}$ model using the current Planck+BK+BAO+$H_{0}$ data. We find that, in this case, the scale invariance is only excluded at the 3.3$\sigma$ level and $\Delta N_{\rm eff}>0$ is favored at the 1.6$\sigma$ level. We then compare the obtained 1$\sigma$ and 2$\sigma$ contours of $(n_{s},r)$ with the theoretical predictions of some selected typical inflation models. We find that, in this case, both the convex and concave potentials are favored at the 2$\sigma$ level, although the concave potential is more favored. The natural inflation model is now excluded at more than 2$\sigma$ level, the Starobinsky $R^{2}$ inflation model becomes only favored at around 2$\sigma$ level, and the most favored model becomes the SBS inflation model.

Acknowledgements.

This work was supported by the National Natural Science Foundation of China (Grants No. 11522540 and No. 11690021), the National Program for Support of Top-notch Young Professionals, and the Provincial Department of Education of Liaoning (Grant No.  L2012087).