The evolution of the primordial curvature perturbation in the ultraslow-roll inflation

A plateau on the background inflaton potential can lead inflation into the ultraslow-roll (USR) stage, during which the primordial curvature perturbation R\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{{\mathscr {R}}}}$$\end{document} can still remarkably evolve outside the horizon. In this paper, the evolution of the Fourier mode Rk\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{{\mathscr {R}}}}_k$$\end{document} and the resulting power spectrum PR(k)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\mathscr {P}}}_{{{\mathscr {R}}}}(k)$$\end{document} are systematically investigated. It is found that, when a scale crosses the horizon well before or after the USR stage, Rk\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{{\mathscr {R}}}}_k$$\end{document} only revolves around the origin in the complex plane. In contrast, if a scale crosses the horizon around the start of the USR stage, the argument of Rk\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{{\mathscr {R}}}}_k$$\end{document} will be nearly frozen and Rk\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{{\mathscr {R}}}}_k$$\end{document} will have a subsequent linear evolution accordingly, in which it may stop rather near or extremely far away from the origin, naturally explaining the entire growing stage of PR(k)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\mathscr {P}}}_{{{\mathscr {R}}}}(k)$$\end{document} on small scales from the minimum to the peak. Moreover, there is a sharp dip in PR(k)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\mathscr {P}}}_{{{\mathscr {R}}}}(k)$$\end{document}, but the lowest point cannot reach 0 exactly. Our work provides a whole picture for Rk\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{{\mathscr {R}}}}_k$$\end{document} and PR(k)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\mathscr {P}}}_{{{\mathscr {R}}}}(k)$$\end{document} in the USR inflation and will help our understanding of the relevant primordial black hole and gravitational wave physics.


Introduction
The detection of gravitational waves from the merger of binary black holes marked the new epoch of multimessenger astronomy [1].The characteristics of these black holes, such as their unexpectedly large masses and relatively small spins, suggest that they may not be usual astrophysical ones, but are of primordial origin [2][3][4].Primordial black holes (PBHs) are formed before big bang nucleosynthesis, so they are a powerful probe in the early Universe [5]).At the same time, the first-order scalar perturbations that generate PBHs can also serve as the source of the second-order scalar-induced gravitational waves (SIGWs) [6,7].Today, the Hawking radiation from PBHs can change the background intensities of a e-mail: linan@mail.neu.edu.cn(corresponding author) various cosmic rays [8].Moreover, PBHs are a natural and promising candidate of dark matter (DM) [9,10], and can act as the seeds of supermassive black holes in galactic centers [11], supported by the recent observation of the James Webb Space Telescope on the massive galaxies at high redshifts [12].
Different from astrophysical black holes, the mass range of PBHs can span almost 50 orders of magnitude, from the Planck mass to supermassive scale (10 −38 -10 10 M , with M = 1.989 × 10 30 kg being the solar mass).Because of the Hawking radiation, the PBHs with mass M < 5 × 10 −19 M have already evaporated, but the PBHs with mass M > 5 × 10 −19 M can still exist today, acting as a stable and pressureless candidate of DM.The PBH abundance f PBH is defined as its proportion in DM today.If f PBH 0.1, the PBHs can be considered as an effective candidate of DM; if f PBH 10 −3 , its possibility as DM can be safely excluded in the relevant mass range.Albeit various constraints from evaporation, lensing, dynamics, accretion, and gravitational waves have strictly constrained the upper bounds of f PBH in different mass ranges, there remains an open mass window from 10 −17 M (asteroid mass range) to 10 −13 M (sublunar mass range), where PBHs are possible to compose all DM.For more recent and complete constraints on f PBH , see Refs.[8][9][10] and the references therein.
The PBH abundance f PBH can be calculated from the power spectrum P R (k) of the primordial curvature perturbation R via the Press-Schechter theory [13] or the more general peak theory [14,15].If f PBH ∼ 0.1, R must be large enough, so that P R (k) can be significantly enhanced on small scales (e.g., 10 8 -10 18 Mpc −1 ), with a peak up to O(10 −2 ) at least.However, such a large magnitude cannot be achieved in the usual single-field slow-roll (SR) inflation models, which only generate a nearly scale-invariant power spectrum on large scales.For example, P R (k) ∼ O(10 −9 ) on the scale of cosmic microwave background [16].There-fore, in order to produce abundant PBHs, the SR conditions must be violated on small scales, and this can be realized in the so-called ultraslow-roll (USR) inflation models .
The aim of this paper is to study the evolution of the primordial curvature perturbation R and the resulting power spectrum P R (k) both theoretically and numerically.Although the shape and amplitude of P R (k) have been obtained in various USR inflation models, there remain several important but unclear issues not thoroughly explored before [56,57].For instance, there is always a sharp dip in P R (k) before its steep growth.Of course, the detail of this dip does not affect the PBH abundance and SIGW spectrum, as they are relevant only to the peak of P R (k).In Refs.[58][59][60][61][62][63][64][65][66], the location of the dip was calculated numerically or semi-analytically, but the underlying physical explanation seems still lacking.
In this paper, we will provide a new perspective on the evolution of R in the USR inflation.Our main improvements to previous works are threefold.First, we investigate the Fourier mode R k and its derivative R k,N with respect to the number of e-folds N in the complex plane.In this way, the entire growing stage of P R (k) from the dip to the peak is demonstrated as a natural consequence from an interesting linear evolution of R k .Second, by illustrating the evolutions of θ k and ϕ k (i.e., the arguments of R k and R k,N ), we confirm the linear evolution of R k , indicating that it is the argument θ k rather than the modulus |R k | that is frozen once the relevant scale crosses the horizon in the USR inflation.Third, we carefully study the minimum of P R (k) in the dip and find that it cannot reach 0 exactly, which is also seldom mentioned before.Altogether, we wish to give a whole picture and thorough understanding of R k and P R (k), which will be beneficial to the model building of the USR inflation and help our research on PBH and gravitational wave physics.
Below, we will follow our previous works on the USR inflation [38,39,53], in which we considered an antisymmetric perturbation δV on the background inflaton potential V b in the USR region.There are several advantages in such a construction.First, the USR stage can be separately studied from the SR stage, so it does not spoil the nearly scale-invariant P R (k) on large scales.Second, δV can be connected to V b very smoothly on both sides of the USR region, and the inflaton can certainly surmount the perturbation.Third, there is no modulated oscillation in P R (k), avoiding the overproduction of tiny PBHs.Consequently, the fine-tuning problem frequently met in PBH physics can be largely relieved.In one word, the antisymmetric δV can create a plateau flat enough on V b , inducing a sufficiently long duration of the USR stage, in which P R (k) can be greatly enhanced.
This paper is organized as follows.In Sect.2, we show the basic equation of motion for the primordial curvature perturbation R k and study its power spectrum P R (k) in the USR inflation.In Sect.3, the evolutions of R k for five typical scales with different times of horizon-exit are systematically investigated in the complex plane.The resulting P R (k) is presented in Sect.4, and several notable features in P R (k) are also explained in detail.We conclude in Sect. 5. We work in the natural system of units and set c = h = k B = 1.

Primordial curvature perturbation and power spectrum
In this section, we show the equation of motion for the primordial curvature perturbation R k and discuss the basic property of its power spectrum P R (k).

Basic equations
We start from the single-field inflation model, with the inflaton minimally coupled to gravity, and the relevant action reads where φ is the inflaton field, V (φ) is its potential, R is the Ricci scalar, and m P = 1/ √ 8π G is the reduced Planck mass, respectively.
In the inflationary era, besides the cosmic time t, it is more convenient to utilize the number of e-folds N as the time variable in cosmic dynamics, which is defined as dN = H (t) dt = d ln a, with a = e N being the scale factor, and H = ȧ/a being the Hubble expansion rate.In order to characterize the motion of the inflaton, two important parameters are introduced as In the usual SR inflation, ε and |η| are always much smaller than 1 and are thus named as the SR parameters.However, in the USR stage, both SR conditions can be broken, and their values may be greatly changed, leaving significant influences on cosmic evolution, PBH abundance, and SIGW spectrum.
In terms of the SR parameters, the Friedmann equation for cosmic expansion can be expressed as and the Klein-Gordon equation for the inflaton field φ can be written as

Primordial curvature perturbation
Now, we turn to the perturbations on the background spacetime.In Newtonian gauge, the metric with the scalar perturbation Φ reads Here, we neglect the vector perturbation and anisotropic stress, as they are irrelevant to PBH and SIGW physics.Furthermore, a more useful and gauge-invariant primordial curvature perturbation R is defined as and the equation of motion for its Fourier mode R k is the Mukhanov-Sasaki (MS) equation [67,68], The initial conditions of the MS equation can be fixed by the Bunch-Davies vacuum condition [69].In the very early Universe, the asymptotic solution of R k is [70] where z = a φ/H = φ ,N e N . 1 Hence, in the very early stage of inflation when k a H = He N , from Eq. ( 3), we have where N ini is the initial number of e-folds.For convenience of the numerical calculations in Sect.3, it is sufficient to set 1 In Ref. [70], R k is given in terms of the conformal time τ as where dτ = dt/a = dN /(He N ), and where N out is the number of e-folds when the relevant scale k crosses the horizon (i.e., k = a(N out )H (N out ) = H (N out )e N out ).Moreover, we can further decompose R k in the complex plane as where |R k | and θ k are the modulus and argument of R k .Substituting Eq. ( 6) into Eq.( 2), we obtain the equation of motion for |R k | as From Eq. ( 4), the initial conditions of Eq. ( 7) read In general, the inflaton φ rolls down from its potential V (φ), so z = φ ,N e N < 0 when N → N ini .Considering the positivity of |R k |, we have assigned a negative sign in Eq. ( 8).
Next, taking into account the asymptotic solutions in Eqs. ( 4) and ( 5), we obtain the relation between θ k and |R k | as Furthermore, the derivative of R k can also be decomposed as where the argument ϕ k can be obtained from Eq. ( 6) as Substituting Eq. ( 10) into Eq.( 11), after some algebra, we obtain The influences of θ k and ϕ k on the evolution of R k will be explored in Sect.3.3.

Power spectrum
The statistical property of cosmological perturbations can be encoded in the two-point correlation function or equivalently the power spectrum in the Fourier space.The dimensionless power spectrum P R (k) of the primordial curvature perturbation R k is defined as In the SR inflation, R k becomes almost frozen once the scale k crosses the horizon, so P R (k) can be effectively calculated at k = a H for simplicity.However, in the USR inflation, R k can still significantly evolve after the horizon-exit, so P R (k) must be evaluated at the end of inflation when k a H.In the SR inflation, on large scales, P R (k) can be usually formulated in a power-law form as where k * is the pivot scale, A s is the amplitude of P R (k) at k * , and n s is the scalar spectral index, respectively.Moreover, n s and the tensor-to-scalar ratio r can be expressed in terms of the SR parameters as n s = 1 − 4ε + 2η and r = 16ε.In the Planck satellite experiment, on the pivot scale k * = 0.05 Mpc −1 , we have the central values A s = 2.10 × 10 −9 and n s = 0.965±0.008,and also r < 0.06 at 95% confidence level [16].
However, on small scales where PBHs and SIGWs are expected, the constraints on P R (k) are still rather loose [5,[8][9][10] (the observation of cosmic microwave background only covers about 10 e-folds on large scales).In the USR inflation, if there exist PBHs with considerable abundance f PBH ∼ 0.1, P R (k) must possess a peak up to O(10 −2 ) on small scales.In addition, n s − 1 is approximately equal to 4 in the growing stage of P R (k) [31,58,61], and there is even a sharp dip before the steep growth of P R (k), all to be discussed in Sect. 4.

Evolution of R k in the complex plane
In this section, we first construct an USR inflation model, with an antisymmetric perturbation on the background inflaton potential in the USR region.Then, the evolutions of R k for five typical scales are systematically studied in the complex plane.Furthermore, these behaviors are reinterpreted by exploring the evolutions of the arguments θ k and ϕ k .
By numerically solving the Friedmann and Klein-Gordon equations, we obtain the two parameters ε and η throughout inflation.In the SR stage, ε and |η| are both much smaller than 1 and vary smoothly.In contrast, during the USR stage, their values change dramatically and may even approach O(1), as shown in Fig. 1 The reason for choosing these numbers is that they correspond to five typical positions in the power spectrum P R (k) (i.e., the nearly scale-invariant region, sharp dip, steep growth, peak, and falling stage, respectively, to be shown in Fig. 9).Our basic results are summarized in Fig. 2, and the detailed analysis of each case can be found in Figs.
Below, we analyze these five typical cases in order.First, we consider the scale with N out = 50, which crosses the horizon much earlier than the start of the USR stage.In Fig. 3, the evolution of R k from N = 47 to 52 is shown in the complex plane, and the difference of the number of e-folds between two adjacent dots is set to be ΔN = 0.01.It can be observed that R k revolves clockwise from outside to the origin and eventually stops near it.The interval between the dots gradually decreases from outside to inside, meaning Fig. 3 The evolution of R k from N = 47 to 52 for the scale with N out = 50.The horizontal and vertical axes are the real and imaginary parts of R k , respectively.The number of e-folds N is illustrated in color, and the difference ΔN between two adjacent dots is set to be 0.01.In the complex plane, R k revolves clockwise around the origin and eventually stagnates near it.Also, it decelerates monotonically, as the interval between the dots decreases from outside to inside.This behavior indicates that, when N out is much smaller than N s , |R k | is almost constant once the relevant scale crosses the horizon, and the influence from the USR stage is negligible that R k decelerates monotonically, consistent with the trend of the black line before N s in Fig. 2b.This can also be seen in Fig. 2a, in which |R k | remains almost constant after N out , indicating that the USR stage has negligible effect on the large-scale curvature perturbation.
Second, we discuss an important and complicated case with N out = 52.89.In the present situation, the USR stage remarkably influences the evolution of R k , as shown in Fig. 4, with N ∈ (50, 62.9).From Fig. 4a, in the early phase before   4 The evolution of R k from N = 50 to 62.9 for the scale with N out = 52.89.In Fig. 4a, before the USR stage, the evolution of R k still resembles that in Fig. 3 (ΔN = 0.01).In Fig. 4b, due to the influence of the USR stage, the horizontal and vertical coordinates of the dots are highly linearly correlated (from N = 56.8 to 62.9), with the Pearson correlation coefficient being 0.999999999999984 (ΔN = 0.1).In this linear evolution, R k eventually stops at a position extremely near the origin.This characteristic behavior causes a remarkable drop of |R k | and thus explains the notable dip in P R (k) N s (from N = 50 to 56.8), the behavior of R k basically resembles that in Fig. 3.However, once the USR stage begins, R k first stagnates for a fairly long time and then evolves along a perfectly straight line in the complex plane, as shown in Fig. 4b.During this unusual linear evolution from N = 56.8 to 62.9, R k moves directly towards the origin and finally stops extremely near it, inducing a significant decrease of |R k | (roughly between N = 60 and 62, as can also be seen in Fig. 2a).It is this decrease of |R k | that naturally explains the notable dip in the power spectrum P R (k) in Fig. 9.
Third, we study the scale with N out = 55.This scale also crosses the horizon before N s , but there is a distinction from the former cases.The evolution of R k from N = 52 to 65 is shown in Fig. 5.We clearly observe that the linear evolution at present is rather different from that in Fig. 4b.Instead of stopping near the origin, R k first approaches the origin, then departs from it, and finally converges to a faraway position.This seemingly strange phenomenon can be understood as follows.From Fig. 2b, the speed |R k,N | for the N out = 55 case (red line) is always larger than that for the N out = 52.89case (blue line), and this is the essential reason why R k moves farther now.Moreover, the behavior of R k causes |R k | to first fall and then rise again to an even higher value, as shown in the red line in Fig. 2a.This process occurs roughly from N = 59 to 61 and corresponds to the growing stage of P R (k) in Fig. 9.
Fourth, we discuss another special scale with N out = 58.84,which lies within the USR stage.In Fig. 6, we show the evolution of R k from N = 55.5 to 62. Different from the three previous cases, R k no longer approaches the origin.In fact, after the USR stage starts, R k almost immediately enters the linear evolution but with a relatively large distance to the origin (see Fig. 6a) and eventually stops at an extremely distant point (see Fig. 6b), consistent with the trend of the purple line in Fig. 2. Consequently, the final value of is the largest among all five cases, corresponding to the peak of P R (k) in Fig. 9.
Last, we set N out = 62, which is a typical case with the relevant scale crossing the horizon after the USR stage.We show the evolution of R k in Fig. 7. Compared with the above four cases, the most distinguishing character here is that R k first revolves away from the origin and then turn backwards, but it neither stagnates near the origin nor possesses a further linear evolution.As a result, |R k | first increases, then decreases, and finally tends to a relatively large value.This can also be seen from the orange line in Fig. 2a and thus explains the falling stage of P R (k) in Fig. 9.

Evolutions of θ k and ϕ k
Below, we adopt a new point of view to reinterpret the evolution of R k , especially for its linear evolution process.In First, θ k,N can be solved from Eqs. ( 7) and (10).In Fig. 8a, we observe that all |θ k,N | drops before N s , meaning that the angular velocity of R k around the origin decreases monotonically.When N > N s , the drop of |θ k,N | continues for the small-scale cases with N out = 58.84(purple line) and 62 (orange line).This means that the angular velocity of R k keeps decreasing, and R k will never be close to the origin, as can be seen in Figs. 6 and 7.However, for the large-scale cases with N out = 50 (black line), 52.89 (blue line), and 55 (red line), |θ k,N | increases in the USR stage, and R k can approach the origin, as shown in Figs. 3, 4 and 5.More interestingly, there is a cusp in the red line with N out = 55.This spiky behavior of |θ k,N | induces a jump of nearly π for the argument θ k , naturally indicating the linear evolution of R k by the origin (but still with a tiny distance), totally in agreement with Fig. 5.
Similarly, ϕ k,N can be solved from Eqs. ( 7) and ( 12).In Fig. 8b, we clearly find that, except for the case with N out = 62 (orange line), all other |ϕ k,N | greatly decreases for a period of time during the USR stage, meaning that ϕ k almost keeps constant in these cases.As ϕ k describes the

Power spectrum in the USR inflation
Now, in Fig. 9, we summarize all our results of the primordial curvature perturbation R k in Sect. 3 in the power spectrum P R (k).Here, we should mention that, besides R k , k 3 is also an essential factor that determines the shape of P R (k), as shown in Fig. 2a.On the one hand, on large scales (e.g., the pivot scale k * = 0.05 Mpc −1 in the Planck satellite experiment), P R (k) is nearly scale-invariant with an amplitude of 2.10 × 10 −9 [16], because the relevant scales cross the horizon so early that the USR stage has little effect on them.On the other hand, on small scales (e.g., k > 10 10 Mpc −1 ), the five typical scales are marked as dots with different colors, and several notable features can be found in P R (k).
First, there is a high peak up to 0.0222 in P R (k) for the scale with N out = 58.84.As a result, there can be PBHs of 10 −17 M with the abundance f PBH ∼ 0.1.This can be understood from Fig. 6, in which the linear evolution of R k during the USR stage makes it converge to a position rather far away from the origin and thus enhances P R (k) significantly.
Second, the steep growth of P R (k) is located around k ∼ 10 13 Mpc −1 with N out ∼ 55, where the scalar spectral index n s − 1 is around 4, in agreement with its theoretical value [31,58,61].
Third, a very noteworthy aspect of P R (k) is that there is a sharp dip on the scale with N out = 52.89.We attribute the formation of this dip to the linear evolution of R k in Fig. 4b.In this situation, R k stops at a point extremely near Fig. 9 The power spectrum P R (k) in the USR inflation.On large scales, P R (k) is nearly scale-invariant with an amplitude of 2.10 × 10 −9 (the asterisk denotes the pivot scale k * = 0.05 Mpc −1 in the Planck satellite experiment [16]).The five typical scales with different N out = 50, 52.89, 55, 58.84, and 62 are marked as dots with different colors.The height of the peak of P R (k) reaches 0.0222, corresponding to N out = 58.84.The slope of P R (k) around N out ∼ 55 is approximately 4, consistent with its theoretical value [31,58,61].Moreover, there is a sharp dip in P R (k) before its steep growth, with the minimum being 1.61 × 10 −14 , corresponding to N out = 52.89 the origin, so it greatly decreases P R (k) accordingly.Here, we should also point out that, albeit the values of |R k | at the end of inflation significantly differ for the scales with N out = 52.89,55, and 58.84, there is no essential difference between them.The unique distinction is the linear evolution of R k that converges at quite different final positions.
Last, an interesting issue worth exploring is whether the dip can reach 0 exactly.The answer is actually no.If it could, there must be a special scale, with R k located just at the origin at the end of inflation.Therefore, during the linear evolution of R k towards the origin, its argument θ k must be a constant, so θ k,N = 0.However, from Eq. ( 10), θ k,N can never be 0. Another easy way to clarify this point is to perform a linear fit of the real and imaginary parts of R k in its linear evolution and then check whether the intercept b vanishes.If not, |R k | is always nonzero, so the dip can never reach 0 exactly.(Even if b vanishes for a special N out , we still need to examine whether the linear evolution of R k really terminates at the origin or merely passes through it.)In our model, the lowest point of the dip corresponds to N out = 52.89.Hence, we plot the intercept b as a function of N out ∈ (52.4,53.4) in Fig. 10.We clearly see that b is positive-definite, so it is sufficient to allow us to conclude that the lowest point of the dip is definitely not 0. In fact, the minimum of P R (k) is 1.61 × 10 −14 , around 10 −5 that of its pivot value with our model parameters, as shown in Fig. 9.

Conclusion
The USR inflation is receiving increasing research interest in recent years.The fundamental reason is that the primordial curvature perturbation R k is no longer frozen after the horizon-exit but can still increase significantly, so a notable peak emerges in the power spectrum P R (k) on small scales.As a result, there can be abundant PBHs at certain mass windows as an effective candidate of DM and also intensive SIGWs at certain frequencies to be hopefully discovered by the next-generation ground-and space-based detectors.
In this paper, we systematically investigate the evolution of R k during the USR inflation and its influence on P R (k).Five typical scales with different numbers of e-folds N out when they cross the horizon are investigated in order, which correspond to five typical regions in P R (k) (i.e., the nearly scale-invariant region, sharp dip, steep growth, peak, and falling stage, respectively).Our basic conclusions can be summarized as follows.
First, different from previous studies mainly focused on |R k | at the end of inflation, we provide an alternative perspective on the entire evolution of R k by solving the MS equation in the complex plane.Only in this way can we thoroughly understand the interesting behaviors of R k in the revolving process around the origin and the possible linear evolution during the USR stage.
Second, for the very large scale with N out = 50 smaller than N s = 56.81 or the very small scale with N out = 62 larger than N e = 60.93,there is only a revolving process around the origin.Therefore, R k eventually stops at a position near the origin.For the former, this corresponds to the nearly scaleinvariant region of P R (k), and for the latter, this corresponds to the falling stage of P R (k).
Third, for the moderate scales with N out = 52.89,55, and 58.84 around N s , besides the early revolving process, there appears a subsequent linear evolution of R k during the USR stage.This means that, different from the fact that the modulus |R k | is frozen after the horizon-exit in the usual SR inflation, now it is the argument θ k that becomes frozen (or with a jump of nearly π ), once the scale crosses the horizon in the USR inflation.Nevertheless, R k eventually terminates at quite different positions.When N out = 52.89,R k converges at a position extremely near the origin; when N out = 58.84,R k converges at a position extremely far away from the origin; the case with N out = 55 lies in the middle.In all, these three situations actually stem from the same reason of linear evolution, so they can naturally explain the entire growing stage of P R (k) from the dip to the peak.
Fourth, we explore the evolutions of θ k and ϕ k (i.e., the arguments of R k and R k,N ), which were seldom discussed before.All |θ k,N | drops before N s , but during the USR stage, it decreases for the small-scale cases but increases for the large-scales ones (even with a cusp when N out = 55).These behaviors indicate the constancy or a jump of π for θ k , resulting in the linear evolution of R k by the origin but not through it.Moreover, |ϕ k,N | monotonically decreases during the USR stage, and an almost constant ϕ k strengthens the conclusion above.
Last, we also discuss whether the lowest point of P R (k) can reach 0 exactly.We linearly fit the real and imaginary parts of R k and illustrate the intercept b as a function of N out around the dip.It is found that b is positive-definite, so the minimum of P R (k) cannot be 0. In fact, this is a general result and can be further verified by looking at the phase diagram in the R k -R k,N plane, to be shown elsewhere in our future work.
Altogether, by studying the evolution of the primordial curvature perturbation R k in the complex plane, we wish to provide a whole picture and thorough understanding of the power spectrum P R (k) in the USR inflation.This work will be helpful to the model building of the USR inflation and contribute our knowledge of PBH and gravitational wave physics.
. The starting and ending points of the USR stage are determined by η(N s ) = η(N e ) = 0.With the model parameters listed above, we have N s = 56.81,N e = 60.93.3.2Evolution of R kNow, we start to investigate the evolution of the primordial curvature perturbation R k .Since R k is complex, we need to numerically solve the MS equation in the complex plane, in which |R k | indicates the distance from the origin to the position of R k , and |R k,N | is the speed of motion of R k .In the following, the evolutions of |R k | and |R k,N | are systematically studied for five typical scales, for which the numbers of e-folds when they cross the horizon are N out = 50, 52.89,55, 58.84, 62.

Fig. 1
Fig. 1 The evolutions of ε and |η| as a function of the number of e-folds N .In the early SR stage, ε ∼ O(10 −4 ), and |η| ∼ O(10 −2 ).The starting and ending points of the USR stage are determined by η(N s ) = η(N e ) = 0, with N s = 56.81 and N e = 60.93.In the USR stage, ε drops almost exponentially, and |η| may even reach O(1)

Fig. 2
Fig. 2 The evolutions of k 3 |R k | 2 /(2π 2 ) and |R k,N | from N = 49 to 63 for five typical scales, with the numbers of e-folds when they cross the horizon being N out = 50 (black line), 52.89 (blue line), 55 (red line), 58.84 (purple line), and 62 (orange line), respectively.The vertical dashed lines indicate the different N out for each case, with

Fig.
Fig.4 The evolution of R k from N = 50 to 62.9 for the scale with N out = 52.89.In Fig.4a, before the USR stage, the evolution of R k still resembles that in Fig.3(ΔN = 0.01).In Fig.4b, due to the influence of the USR stage, the horizontal and vertical coordinates of the dots are highly linearly correlated (from N = 56.8 to 62.9), with the Pearson

Fig. 5
Fig. 5 The evolution of R k from N = 52 to 65 for the scale with N out = 55 (ΔN = 0.01).The Pearson correlation coefficient for the linear evolution process with N ∈ (56, 65) is 0.99999999793965.Different from Fig. 4b, R k no longer stops near the origin but at a faraway position, leading to the steep growth of P R (k).(It should be noted that the linear evolution does not intersect with the origin.)

Fig. 6 999999996557677 Fig. 7
Fig. 6 The evolution of R k from N = 55.5 to 62 for the case with N out = 58.84,with the relevant scale crossing the horizon within the USR stage.In Fig. 6a, R k first revolves around the origin, but does not come very close to it any longer.With the beginning of the USR stage, R k almost immediately starts the linear evolution (ΔN = 0.01).The subsequent linear evolution from N = 58.8 to 62 is shown in

Fig. 8 ,
we show |θ k,N | and |ϕ k,N | (i.e., the time derivatives of the arguments of R k and R k,N ) as a function of N , which characterize the changes in the direction of motion of R k and R k,N , respectively.

Fig. 8
Fig. 8 The evolutions of |θ k,N | and |ϕ k,N | from N = 49 to 63 for five typical scales with different N out .In Fig. 8a, when N < N s , all |θ k,N | drops, so the angular velocity of R k around the origin decreases.When N > N s , this tendency remains for the small-scale cases with N out = 58.84(purple line) and 62 (orange line), so the angular velocity of R k keeps decreasing, and R k can never be close to the origin.In

Fig. 10
Fig. 10 The intercept b of the linear fit of the real and imaginary parts of R k in its linear evolution as a function of N out .The positivity of b indicates that the lowest point of the dip in P R (k) cannot reach 0 exactly.As the vicinity of the sharp dip is very narrow around N out = 52.89,we plot the figure with N out ∈ (52.4,53.4) here.(Actually, b > 0 even in a large domain from N out = 51 to 60.)