Measuring the spectral running from cosmic microwave background and primordial black holes

We constrain the spectral running by combining cosmic microwave background (CMB) data, baryon acoustic oscillation (BAO) data and the constraint from primordial black holes (PBHs). We find that the constraint from PBHs has a significant impact on the running of running of scalar spectral index, and a power-law scalar power spectrum without running is consistent with observational data once the constraint from PBHs is taken into account. In addition, from the constraints on the slow-roll parameters, the derived tensor spectral index in the single-field slow-roll inflation model is quite small, namely $|n_t|\lesssim 9.3\times 10^{-3}$ which will be very difficult to be measured by CMB data only in the future, and the absolute value of derived running of tensor spectral index is not larger than $2.1\times 10^{-4}$ at 95% confidence level.


I. INTRODUCTION
The current cosmological observations, for example the Planck satellite [1,2], provide strong evidence for the standard six-parameter ΛCDM model. However, when the base model is extended and other cosmological parameters are allowed to freely vary, a few anomalies are present. For example, the running of running of the scalar spectral index which is defined by reads β s = 0.025 ± 0.013 (2) from Planck data in [2] and β s = 0.021 ± 0.013 (3) from the combination of Planck [2], BK14 [3] and Baryon Acoustic Oscillation (BAO) [4][5][6] data in [7] at 68% confidence level (CL) on the cosmic scales, where n s is the spectral index of scalar power spectrum.
In 1971, Hawking proposed that the highly overdense region of inhomogeneities would eventually cease expanding and collapse into a black hole in [8]. A large positive value of the running of running may lead to the formation of primordial black holes (PBH) at small scales * lijun@itp.ac.cn † huangqg@itp.ac.cn in the early universe. Different from the astrophysical black holes which should be heavier than a particular mass (around 3 solar mass M ), PBHs can have very small masses. However, the PBHs with masses smaller than 10 −18 M would have completely evaporated by now due to the Hawking radiation. In particular, the quantum emissions from PBHs with mass around 10 −18 M can generate a γ-ray background which should be observed today. Conversely, the non-detection of such a γ-ray background in our Universe provides a stringent constraint on the PBHs with such a mass. Recently various constraints on the abundance of PBHs are present in [9][10][11][12][13][14][15][16][17][18][19][20][21][22]. Ones expect that the absence of PBHs should tightly constrain the running of running of scalar spectral index.
On the other hand, inflation [23][24][25][26] is taken as the leading paradigm for the physics in the very early universe. The initial inhomogenieties and spatial curvature are supposed to be stretched away by the quasiexponential expansion of inflation. Because the Hubble parameter is roughly a constant during inflation, it predicts a nearly scale-invariant power spectrum of the curvature (scalar) perturbations seeded by the quantum fluctuations of inflaton field during inflation. Typically, the scalar spectral index n s is related to the number of e-folds N before the end of inflation by n s − 1 ∼ −1/N , and then β s ∼ −1/N 3 ∼ (n s − 1) 3 ∼ −O(10 −5 ). In this sense, the running of running of scalar spectra index given in Eqs. (2) and (3) seems too large to fit in the typical inflation models.
In this paper we adopt two methods to compare inflationary predictions with current cosmological datasets, in arXiv:1806.01440v1 [astro-ph.CO] 5 Jun 2018 particular including the constraint from the PBHs. The first method consists of a phenomenological parameterization of the primordial spectra of both scalar and tensor perturbations, and the second exploits the analytic slowroll-parameter dependence of primordial perturbations. This paper is organized as follows. In Sec. 2, we explain the datasets used in this paper. In Sec. 3 and Sec. 4, we present the constraints on the phenomenological parameters and the slow-roll parameters with the publicly available codes CosmoMC [27], respectively. A brief summary will be given in Sec. 5.

II. DATA
The full-mission Planck observes the temperature and polarization anisotropies of the cosmic microwave background (CMB) radiation. In this paper we also combine all the data taken by the BICEP2 and Keck Array CMB polarization experiments up to and including the 2014 observing season [3] with the Planck data.
Baryon Acoustic Oscillation (BAO) detections measure the correlation function and power spectrum in the clustering of galaxies. Measuring the position of these oscillations in the matter power spectra at different redshifts removes degeneracies in the interpretation of the CMB anisotropies. The BAO data adopted in this paper include 6dFGS [5], MGS [6], BOSS DR11 Lyα [28], BOSS DR12 with nine anisotropic measurements [29], and eBOSS DR14 [30].
The mass of PBH is roughly given by the horizon mass at the time of formation, namely where γ = 3 −3/2 0.2 [31]. In the comoving units, the horizon scale is R = (aH) −1 ∝ a during radiation domination. Since constant-entropy expansion implies where the number of relativistic degrees of freedom g * is around 3 at equality and 10 2 in the early universe. Therefore where k 10 6 Mpc −1 is the perturbation mode which re-enter the horizon when the one-solar mass PBH can be formed. Therefore the scale correspond to the formation of PBH with a mass 10 −18 M is roughly given by k c 10 15 Mpc −1 . Actually the formation process of the PBHs is still poorly understood. As a conservative estimation, we require that the scalar power spectrum should smaller than the unity at k = k c , namely Otherwise, a lot of PBHs with mass 10 −18 M should be formed in the early Universe and generate an observable γ-ray background. We take Eq. (6) as the constraint on the scalar power spectrum from PBHs.

III. CONSTRAINTS ON THE SPECTRAL RUNNING FROM CMB AND PBHS
In this section, the power spectra of the scalar and tensor perturbations are parameterized by where A s (A t ) is the scalar (tensor) amplitude at the pivot scale k * = 0.05 Mpc −1 , n s is the scalar spectral index, α s ≡ dn s /d ln k is the running of scalar spectral index, β s ≡ d 2 n s /d ln k 2 is the running of running of scalar spectral index, n t is the tensor spectral index, and α t ≡ dn t /d ln k is the running of tensor spectral index. Usually we introduce a new parameter, namely the tensor-to-scalar ratio r, to quantify the tensor amplitude compared to the scalar amplitude at the pivot scale: It is known that the single parameter Harrison-Zeldovich spectrum (n s = 1) does not fit the data and at least the parameters A s and n s in the expansion of the primordial scalar power spectrum are needed. Here the spectral index of tensor power spectrum is set as n t = −r/8 which is nothing but the consistency relation to lowest order in the single-field slow-roll inflation model, and α t = 0. We consider the six parameters in the standard ΛCDM model, i.e. the baryon density parameter Ω b h 2 , the cold dark matter density Ω c h 2 , the angular size of the horizon at the last scattering surface θ MC , the optical depth τ , the scalar amplitude A s and the scalar spectral index n s . We extend this scenario by adding the running of the scalar spectral index α s , the running of running β s and the tensor amplitude, or equivalently the tensor-to-scalar ratio r. We constrain all of these 9 parameters in the ΛCDM+r + α s + β s model by adopting two different data combinations, namely CMB+BAO and CMB+BAO+PBH, respectively. The results are given in Table. I and Fig. 1. The constraint on the running of running of scalar spectral index is β s = 0.020 ± 0.013 (10) at 68% CL from CMB+BAO, and β s = −0.0038 +0.0069 −0.0025 (11) at 68% CL from CMB+BAO+PBH. We see that the constraint on the running of running of scalar spectral index is significantly affected by adding the constraint of PBHs, and a power-law scalar power spectrum without running is consistent with the data once the constraint from PBHs is taken into account.

IV. CONSTRAINTS ON THE SLOW-ROLL PARAMETERS
In this section, we focus on the cannonical single-field slow-roll inflation model in which the inflation is driven by the inflaton potential V (φ). The dynamics of inflation is govern by where M p = 1/ √ 8πG is the reduced Planck energy scale and the dot and prime denote the derivative with respective to the cosmic time t and the inflaton field φ, respec-tively. The inflaton field slowly rolls down its potential if 1 and |η| 1, where The amplitude of scalar and tensor perturbation power spectra are given by in [32], where c 0.08145. See [33,34] as well. And hence, we have where Here the nine parameters sampled in the CosmoMC are {Ω b h 2 , Ω c h 2 , τ, θ MC , A s , , η, ξ, σ}, and then n s , α s ,   β s , r, n t and α t are all taken as the derived parameters. The constraints on the slow-roll parameters { , η, ξ, σ} and the contour plots of these slow-roll parameters are illustrated in Table. II and Fig. 2. We see that the constraint on the fourth slow-roll parameter σ is improved once the constraint from PBHs is included.
Here we also work out the predictions of the slowroll inflation model constrained by the observational data. The parameters {r, n t , α t } characterizing the tensor power spectrum can be obtained by adopting the Eqs. (18), (19) and (20). The results are showed in Table. II and Fig. 3. In particular, we notice that both the derived tensor spectral index n t and its running α t are negative in the constrained single-field inflation model, and the constraint on them at 95% C.L. are from CMB+BAO, and −n t < 9.3 × 10 −3 , from CMB+BAO+PBH. However, in [35], the optimistic estimation indicates that the uncertainty at 68% C.L. on the tensor spectral index is σ nt 1.1 × 10 −2 due to the cosmic variance only for CMB multipoles less than 300. Because the small-scale CMB B-modes are dominated by CMB lensing, it implies that it is very difficult to measure the tensor spectral index for the single-field slowroll inflation model by only using the CMB data.

V. SUMMARY
In this paper we use two methods to constrain the slowroll inflation models by combining CMB, BAO and the constraint from PBHs. Even though a positive running of running of scalar spectral index is slightly preferred by the data combination of CMB and BAO datasets, a power-law scalar power spectrum without running is consistent with the data once the constraint from PBHs is taken into account.
We can also directly constrain the slow-roll parameters from the observational data. An advantage of this method is that we can work out the predictions of singlefield slow-roll inflation model by using these constrained slow-roll parameters. For example, we illustrate the predictions of the parameters characterizing the tensor power spectrum, and find that both the tensor spectral index and its running are negative and their absolute values are not larger than 9.3 × 10 −3 and 2.1 × 10 −4 at 95% C.L., respectively. Our results imply that it is very difficult to measure these two parameters in the future.