On Stochastic Effects and Primordial Black-Hole Formation

The effect of large quantum fluctuations on primordial black-hole formation during an inflationary phase along a plateau is investigated. By using techniques of stochastic inflation in combination with replica field theory and the Feynman-Jensen variational method, it is non-perturbatively demonstrated that the abundance of primordial black holes is amplified by several orders of magnitude as compared to the classical computation.

Cosmological inflation [1][2][3][4] is a fundamental building block of our current understanding of the Universe. In addition to explaining the flatness and homogeneity of the Cosmos, inflation predicts the generation of perturbations from quantum fluctuations in the early Universe. The most common way for a realization of inflation is via a single scalar field, the so-called inflaton ϕ. Quantum fluctuations of the latter as well as of the metric seed today's structure in the Universe. The predictions of inflation are in remarkable agreement with measurements (cf. Refs. [5][6][7]).
Under certain circumstances, the quantum fluctuations of the inflation field can be dominant over its classical evolution. There are two important cases in which this happens. One is typically at larger values of the inflaton potential V(ϕ), yielding eternally expanding patches of the Universe [8][9][10] (for a review see Ref. [11]). The other case occurs when the inflaton potential possesses a (quasi-)inflection point or one (or a multiple) plateaulike feature. This can be realized in numerous scenarios, such as double inflation [12], radiative plateau inflation [13], or in string theory (e.g. Refs. [14,15]).
Classically, using the slow-roll conditions, where an overdot represents a derivative with regard to cosmic time t, H ≡ . a/a is the Hubble parameter, and a is the scale factor, the Universe inflates as with N being the number of e-folds. On the other hand, the corresponding quantum fluctuations are In turn, noting that the primordial metric perturbation implies that quantum effects can expected to be important, whenever ζ becomes of order one. This is often the case in primordial black hole (PBH) formation [16][17][18][19], which will be discussed further below. Flatness of the inflaton potential, i.e. dV(ϕ)/dϕ = 0, leads to a growth of the primordial power spectrum of comoving curvature perturbations where ϕ ≡ dϕ/dN . The classical equation of motion, becomes which, by virtue of Eq. (5), yields We note that, during the a plateau or inflectionpoint phase, the slow-roll conditions Eq. (1) are violated. arXiv:1906.02744v2 [gr-qc] 6 Nov 2019 Whilst these are tantamount to |ε 1,2 | 1, with the slowroll parameters ε 1 ≡ − . H/H 2 and ε 2 ≡ . ε 1 /(ε 1 H), here, the flatness of the potential implies |ε 2 | 6 (cf. Ref. [20]). Furthermore, as will be discussed below, quantum effects will start to dominate.
In order to proceed and to demonstrate the strength of the quantum effects, we begin by choosing an inflaton potential with a quasi-inflection point as shown in Fig. 1. For definiteness, we will follow Ref. [21] and make the choice whereφ ≡ ϕ/M Pl , with M Pl being the Planck mass, and In turn, the power spectrum can be calculated using the Mukhanov-Sasaki equation [22,23], Of course, as mentioned above, quantum effects cannot be ignored during the inflection-point phase. In order to study these situations, Starobinsky introduced a stochastic framework [9]. Its key idea lies in splitting the inflaton ϕ into long-and short-wavelength modes, and viewing the former as classical objects evolving stochastically in an environment provided by quantum fluctuations of shorter wavelengths. Hence, it constitutes an example of how fundamental properties of quantum fields can be modelled using methods of statistical mechanics; one focusses on the "relevant" degrees of freedom (the long-wavelength modes) and regards the shortwavelength modes as "irrelevant" ones, where "short" and "long" are subject to the Hubble horizon.
In turn, the right-hand side of Eq. (10) acquires a stochastic source term h, which is a Gaussian-distributed random variable with (see Ref. [9]) where ∆ is a known function, depending on derivatives of the mode function (see e.g., Ref. [24] for a detailed presentation). The overbar in Eq. (11) and below denotes the average over the noise due to quantum fluctuations of short wavelengths.
Before we proceed, in order to provide an overview, let us briefly summarize the subsequent computational program (cf. Ref. [25]): 1. Wick-rotate to Euclidean signature, 2. use the replica trick [26], in order to 3. apply the Feynman-Jensen variational principle [27]. This will allow us to go beyond ordinary perturbation theory (cf. Ref. [28]) and to obtain the full power spectrum, and therefore the PBH mass spectrum including stochastic modifications.
After Wick-rotating, we proceed with the replica trick (see the Appendix for a proof), where Z[j] is the generating functional depending on an external current j. The different replica copies are labelled by the indices a, b; their number is equal to m. We define the replicated action S (m) via Besides terms diagonal in replica space (∝ δ ab ), it also contains the non-diagonal part originating from the average over noise.
We apply the Feynman-Jensen variation principle, and therefore define a Gaussian variational action with t ≡ dt and k ≡ d 3 k/(2π) 3 . The variational propagator (G ab ) a,b=1,...,m is a matrix in replica space, which is to be determined. In order to do so, we make the ansatz The self-energy matrix σ mimics the diagonal and the non-diagonal parts in Eq. (13).
Maximizing the right-hand side of the Feynman-Jensen inequality wherein the subscript "var" refers to the variational action (14), allows us to extract the replica structure σ, and in turn the full power spectrum. The latter is obtained via (cf. Ref. [28]) where the matrices 1 and 1ll are in replica space, with 1 ab = δ ab , and 1ll ab = 1 for all a, b. Figure 2 shows the results for the full power spectrum P QD (k) with quantum-diffusion effects (red, dot-dashed curve) in comparison to the power spectrum P (k) without (blue, dashed curve). Both of these spectra are enhanced by the inflection point. As can be seen, the resummed stochastic effects yield an increase of power.
Having described a methodology to effectively incorporate the large quantum effects during the inflection-point phase, we are now in a position to derive the probability β(M ) that an overdense region of mass M has a size exceeding the Jeans length at maximum expansion, so that it can collapse against the pressure and form a black hole. If the horizon-scale fluctuations have a Gaussian distribution with dispersion σ, i.e. root-mean square of the primordial power spectrum P , one expects for the fraction β of horizon patches collapsing to a black hole [19] β ≈ Erfc Here, Erfc is the complementary error function Erfc ≡ 1 − Erf, and δ c is the density-contrast threshold for PBH

FIG. 3:
Fractions β (blue, dashed) and βQD (red, dotdashed) of horizon patches collapsing to PBH as functions of e-folds N . As in Fig. 2, the upper curve depicts the results of incorporating quantum diffusion. formation. A simple analytic argument [19] suggests a value of δ c ≈ 1/3, but more precise arguments -both numerical [29] and analytical [30] -suggest a somewhat larger value: δ c = 0.45. We will use this value in remainder of our work. The result for β are given in Fig. 3, both with and without the effect of quantum diffusion. Due to the exponential dependence of β on the power spectrum [see Eq. (19)], already a moderat increase of the latter by less than an O(1) factor, yields an increase of the former by more than 30 orders of magnitude! Clearly, this has consequences for the formation of PBHs as possible dark-matter candidates (see Ref. [31] for a review). In particular, the extendedness of the mass function generally allows for a larger PBH dark-matter fraction than monochromatic cases (see Refs. [29,[32][33][34][35]). We should note that there is room for where the region of enhancement exactly lies. The closer towards the end of inflation it occurs, the lighter the holes are; but also, at the same time, the lower β has to be. This is because during radiation domination its growth is approximately proportional to the scale factor. Hence, in model building, special care needs to be taken in order for quantum diffusion not to over-produce PBHs, and therefore rule out, the model.
In other scenarios and using different approaches, other recent works [36][37][38] also reach the conclusion that quantum diffusion leads to an increase of the power spectrum, and, in turn, of β. However, firstly, the scenario investigated in this work, can be viewed as a general prototype class of inflationary models in which quantum diffusion becomes relevant. Secondly, the utilized methods -the replica trick and the Feynman-Jensen variational principle -are non-perturbative in nature; the stochastic effects are effectively resummed into a generalised selfenergy contribution.
We should note that these technique are not limited to linear noise potentials and could include arbitrary self-interactions (cf. the original Ref. [28] for the context of spin glasses, and Refs. [25,[39][40][41] for its application to stochastic inflation). Furthermore, it can easily be applied to multi-field scenarios. Let us finally point out that the replica method used in this work has a close connection to the functional renormalisation group (cf. Refs. [42,43]  The replica trick is much older than one might think. According to Giogio Parisi (see Ref. [44]) it can be dated back to at least the fourteenth century when the bishop of Lisieux applied a similar trick to define non-integer powers! The first non-trivial application of the replica trick to modern physics has been done by Edwards and Anderson in 1975 [45].