Primordial black holes and secondary gravitational waves from generalized power-law non-canonical inflation with quartic potential

Here, generation of PBHs and secondary GWs from non-canonical inflation with quartic potential have been probed. It is illustrated that, quartic potential in non-canonical setup with a generalized power-law Lagrangian density can source a consistent inflationary era with the latest observational data. Besides, we show that our model satisfies the swampland criteria. At the same time, defining a peaked function of inflaton field as non-canonical mass scale parameter $M(\phi)$ of the Lagrangian, gives rise to slow down the inflaton in a while. In this span, namely Ultra-Slow-Roll (USR) stage, the amplitude of the curvature perturbations on small scales enlarges versus CMB scales. It has been illustrated that, further to the peaked aspect of the chosen non-canonical mass scale parameter, the amount of $\alpha$ parameter of the Lagrangian has enlarging impact on the amplitude of the scalar perturbations. As a consequence of adjusting three parameter Cases of this model, three Cases of PBHs in proper mass scopes to explain LIGO-VIRGO events, microlensing events in OGLE data and DM content in its totality, could be produced. In the end, power-law behavior of the current density parameter of gravitational waves $\Omega_{\rm GW_0}$ in terms of frequency has been examined. Also, the logarithmic power index as $n=3-2/\ln(f_c/f)$ in the infrared regime is obtained.


I. INTRODUCTION
It is widely recognized that, quantum fluctuations of scalar field in the inflationary era could produce primal curvature perturbations.The hastened expansion of the universe during the inflation throws the perturbed modes out of the Hubble horizon.These modes come back to the horizon in Radiation Dominated (RD) era.The enlarged perturbed scales can produce the ultra-dense domains in RD era.In the end, gravitational cave-in of these domains causes to born Primordial Black Holes (PBHs) with different masses.The earliest allusion to the context of PBHs was referred to the 1970s [1][2][3][4].Thenceforth, fortunate detection of Gravitational Waves (GWs) ensued from converging black holes with masses of 30M ⊙ (M ⊙ connotes the solar mass) by LIGO-Virgo cooperation [5][6][7][8][9], has returned the PBHs notion to the center of attention.So PBHs could be considered as the origin of traced GWs.Into the bargain, the cryptic essence of Dark Matter (DM) [10], has inspired the researchers to ponder about PBHs as a probable candidate for the totality or a portion of DM content .
In most of the PBHs production models, enhancement of the primordial curvature perturbations takes places in an Ultra-Slow-Roll (USR) phase of the inflationary era.In the USR phase due to a friction dominated process, velocity of the inflaton field decreases in comparison with slow-roll regime.Usually, in this phase the slow-roll condition is broken down by the second slow roll parameter.In this way extra time for enough growth of P R to produce PBHs is provided.Basically, the USR regime in the standard model of inflation could be generated by altering the inflationary potentials.This aim could be reached either by altering the potentials to have inflection points [34][35][36][37], or by setting a bump/dip for the potentials on the scales smaller than CMB scale [33,38].Hence, decreasing the velocity of the inflaton could take place in the vicinity of the inflection point or in the bump/dip passing.Beyond the standard model of inflation, the USR regime could be produced via applying suitable coupling functions to enhance the friction gravitationally [14][15][16][17][18][19][20][21][22][23].
The non-canonical inflationary model with power-law Lagrangian L(X, φ) = X α − V (φ), where X connotes the canonical kinetic term, is a well recognized generalization for the standard model of inflation [12,13,61,62,66].Therein, for α = 1 the canonical standard model of inflation is recovered.It has been shown that, smaller obtained slow-roll parameters in this framework in comparison with canonical case can lead to the longer inflationary era and also larger scalar spectral index (n s ) and smaller tensor-to-scalar ratio (r) [62].Thence, the steep potentials such as quartic and exponential potentials could be resurrected in this setup [61,62,66].On the other hand, it has been proven that in this framework allowable PBHs mass spectra and detectable GWs could be produced on small scales [12,13].
So as to check the theoretical rectitude of the inflationary models, their compatibility with swampland criteria in addition to observational constraints have been verified.The swampland criteria emanate from string theory and comprise two conjectures namely distance conjecture and de Sitter conjecture [67][68][69].The distance conjecture imposes an upper limit on the filed evolution ∆φ < 1.The de Sitter conjecture imposes a lower limit on the potential gradient |V ,φ /V | > 1.The de Sitter conjecture is in conflict with the slow-roll condition in the standard model of inflation due to In the present work, we are interested in studying the production of PBHs and GWs in the framework of generalized power-law non-canonical inflation with quartic potential, in which the mass scale parameter of non-canonical Lagrangian M(φ) depends on the scalar field.We try to employ the field-dependent feature of M(φ) to make the model not only be consistent with observations at CMB scale but also produce PBHs on small scales.To do this, in Sec.II, fundamental relations of generalized power-law non-canonical model are analyzed.
Thence, in Sec.III, observational and theoretical viability of the model thereto enhancing process of the amplitude of the scalar power spectrum at small scales has been expounded.
Thereafter, PBHs mass spectra and current density parameter spectra of secondary GWs are obtained in Sec.IV and Sec.V, respectively.In the end, Sec.VI is allotted to the abridged results of the paper.

II. GENERALIZED POWER-LAW NON-CANONICAL INFLATION
We start by the following common action in which the Lagrangian density L(X, φ) could be defined as varied functions of scalar field φ and kinetic term X ≡ 1 2 ∂ µ φ∂ µ φ [61,62,66,70,71].Here, we introduce a new generalized power-law shape for the Lagrangian density as follows therein, the non-canonical dimensionless α parameter specifies deviation from canonicity.
for α = 1 the canonical version of Lagrangian (2) can be retrieved.The parameter M(φ) with dimension of mass, denotes the non-canonical mass scale of inflation.For the first time we consider M(φ) as a general function of scalar field instead of a constant parameter.V (φ) is the inflationary potential.The homogeneous and isotropic universe is described by the ensuing flat Friedmann-Robertson-Walker (FRW) metric as where t is the cosmic time and a(t) denotes the scale factor.Using the flat FRW metric, the kinetic term transforms into X = φ2 /2 (dot signifies derivative against t).
It is noteworthy that, the following field redefinition results in recovery of the well known form of the non-canonical Lagrangian L( X, ψ) = Xα − U(ψ) from the Lagrangian (2), wherein X = ψ2 /2 and U(ψ) = V (φ(ψ)).In this way, any function M(φ) 4 can be eliminated from the non-canonical kinetic term of the Lagrangian (2) by the field redefinition (4).Nevertheless, some form of M(φ) in the field redefinition (4) may cause complexity in deriving an analytical expression for the potential U(ψ).Hence, it is reasonable to use the Lagrangian (2) in such models depending on the form of M(φ).
In what follows the ensuing form of energy density ρ φ and pressure p φ of the scalar field could be computed from the Lagrangian (2) (see [61,62] to peruse these equations) Utilizing Eqs. ( 5) and ( 6), the first and second Friedmann equations are derived as in which H ≡ ȧ/a designates the Hubble parameter and M p = 1/ √ 8πG is the reduced Planck mass.
Using the Lagrangian (2) and minimizing the action (1) with respect to scalar field φ, the ensuing second order equation of motion can be derived wherein ( ,φ ) signifies derivative against φ.It is worth noting that, for α = 1, all of the preceding equations revert to their canonical versions.
In the following, the first and second Hubble slow-roll parameters are introduced as Under the slow-roll approximation conditions ({ε 1 , ε 2 } ≪ 1), kinetic energy term can be dominated by the potential energy.So, thereunder the first Friedmann equation (7) abbre- Furthermore, under the slow-roll conditions in Eq. ( 9), φ and φ2 could be neglected and using Eq. ( 11), the equation of motion can be recast in where θ = 1 for V ,φ > 0 and θ = −1 for V ,φ < 0.
Applying the definition of the first slow-roll parameter (ε 1 ≡ − Ḣ/H 2 ) and using dN = Hdt to convert the cosmic time to e-folding number N thereto Eqs.(7), the field equation of motion (9) takes the following form wherein ( ,N ) and ( ,N N ) designate the first and 2nd derivative against N, respectively.
So as to analyze the perturbations dynamics in this framework, we pursue the calculations of [72].Therefrom, the curvature power spectrum under the slow-roll approximation at sound horizon passing (c s k = aH) by comoving wavenumber k in non-canonical setup is given by where signifies the square of the scalar perturbations sound speed [62,72].The amplitude of the curvature power spectrum at pivot scale (k * = 0.05 Mpc −1 ) has been constrained by Planck [53].For the power-law Lagrangian (2), the square of the sound speed (15) takes the following form It is obvious that, in order to prevent the model from classical instability we required c 2 s > 0 which results in α > 1/2.
Pursuant to [72], in this setup the scalar spectral index n s is obtained from the curvature power spectrum and it could be written with regard to slow-roll parameter as follows The scalar spectral index is observationally constrained by Planck 2018 as n s = 0.9668 ± 0.0037 (TT,TE,EE+lowE+lensing+BK15+BAO, 68% CL) [53].
Subsequent to [72], the tensor perturbations power spectrum under the slow-roll approximation at k = aH in non-canonical framework is given by It is evident that, the tensor power spectrum in this setup is equal to its canonical counterpart.Because The tensor power spectrum in this setup ia associated to the gravity sector of the action.Employing the scalar ( 14) and tensor (18) power spectra, the tensor-to-scalar ratio r is calculated as There is an upper bound on the tensor-to-scalar ratio established by Planck 2018 as r < 0.063 (TT,TE,EE+lowE+lensing+BK15+BAO, 95% CL) [53].It is worth noting that, recently the mentioned upper bound on r is constricted by BICEP/Keck 2018 to r < 0.036 at 95% CL [73].

III. METHODOLOGY AND VIABILITY OF THE MODEL
It is widely recognized that, creation of the seeds of PBHs during the inflation requires an enhancement around seven order of magnitude in the P R on small scales in proportion to P R (k * ) at CMB scale.In the subsequent stage, the circumstances of occurring such enhancement in the P R through the non-canonical framework with power-law Lagrangian (2) has been expounded.To this end, the ensuing two-parted function is chosen for the non-canonical mass scale parameter M(φ) of the generalized Lagrangian ( 2) wherein The M 0 parameter in (20) with the dimension of mass, is responsible for observational compatibility of the model on CMB scale as well as the α parameter.The second term of (20), ǫ(φ) is a dimensionless peaked function of φ, that is able to generate a localized peak in φ = φ c with height ω and breadth b.The parameters {φ c , b} have dimensions of mass and ω is dimensionless.It can be deduced that, for φ = φ c the ǫ(φ) function melts away (ǫ(φ) ≪ 1) and M(φ) comes back to constant M 0 .Hence, ǫ(φ) is responsible for slowing the inflaton down in the vicinity of φ = φ c (i.e.USR phase) and enhancing the curvature power spectrum on small scales without significant effect on CMB scale.So as to achieve this goal, we need to adjust the peak function parameter {ω, b, φ c }.In this way, not only the compatibility of the model with observational data on CMB scale is insured, but also the USR phase could be generated to enhance the P R on smaller scales to produce PBHs.
It should be noted that, in [12] the authors have investigated production of PBH and GWs in non-canonical framework with power-law Lagrangian through the mechanism of inflection point in a class of steep-deformed exponential potential.Two cases of PBHs are produced in the vicinity of inflection point of the potential (the USR stage) [12].
In [13], production of PBHs and GWs in the same framework as [12], from the quartic potential in the presence of a tiny bump has been investigated.At the moment of bumpassing three cases of PBHs were generated [13].In both mentioned papers [12,13] the power-law Lagrangian contains the constant non-canonical mass scale parameter M, but in the present work we consider a generalized form of this framework with M(φ) as a function of φ defined by Eq. ( 20).
Here, we consider the quartic potential for the model in which, λ ≃ 0.13 is the dimensionless self coupling parameter [74].The quartic potential can source chaotic inflation [75] and it has substantial reheating features [76,77].It is demonstrated that, the quartic potential in the standard model of inflation suffers from some problems.First, because of the large self-coupling constant λ ≃ 0.13, its produced scalar fluctuations on CMB scale cannot be compatible with the latest Planck data [53].
On the other hand, owing to generation of the large tensor fluctuations by this potential, its predicted quantity for r cannot lie in the permitted domain of Planck 2018 data [53].
As regards the quartic potential cannot accommodate to feasible inflationary epoch in the standard model of inflation, in this work we attempt to revive it through the generalized power-law non-canonical setup.Simultaneously, employing a peaked function for M(φ) in this setup, we could examine the generation of PBHs and GWs on smaller scales.
Altogether, this model is comprised of three Cases containing a collection of six parameters {α, M 0 , λ, ω, φ c , b}.Regarding λ ≃ 0.13, the M 0 parameter for each case could be computed via the observational constraint on the scalar power spectrum (P * R ∼ 2.1 × 10 −9 ) at pivot scale (k * = 0.05 Mpc −1 ) [53].Three parameters {ω, φ c , d} do not have any significant effect on CMB scale.So α is the only free parameter which can be used to remedy the large tensor-to-scalar ratio r of the model.So pursuant to [61][62][63][64][65][66], we have chosen α = 17 for all Cases of our model to decrease the r.It is notable that, using the chosen α, Null Energy Condition (NEC) is preserved and the kinetic energy remains positive and dominated by potential energy throughout the inflationary era.Thus, the model is free of the ghost and tachyonic instabilities.Furthermore, the Lagrangian (2) satisfies the Effective Field Theory (EFT) conditions for non-canonical Lagrangian [78].The assigned parameters for all Cases of the model are classified in Table I   for Cases of this model (see the ε 1 plot in Fig. 1).
Concerning the background evolution, the second Friedmann equation ( 8) and the equation of motion (13) should be solved concurrently, via substituting the inflationary potential (22) and the defined M(φ) from ( 20)- (21).The required initial conditions are acquired from the slow-roll equations ( 11)-( 12) with the potential (22).Thereafter, evolutions of scalar field φ, its derivative φ ,N , the first and second slow-roll parameters (ε 1 , ε 2 ) against the efolds number have been plotted in Fig. 1 for Cases A (red lines), B (green lines) and C (blue lines).A small USR region in the proximity of φ = φ c can be seen in each graph of Fig 1.
The short-living USR stage lasts for about two e-folds numbers.In this region, the inflaton speed is decreased and goes toward zero (see φ ,N in Fig. 1).Hence, an approximately flat region in the graph of the field evolution is produced for each Case of the model (notice φ in Fig. 1).The inflationary era of our model is comprised of three consecutive phases namely first SR phase, intermediate USR phase and final SR phase.Following [80], the sharpness of transitions between these phases can be inferable from the evolution of ε 2 .As it can be seen from the ε 2 graph of Fig 1, in the beginning of the USR stage ε 2 takes the negative values then it tends to positive values larger than one, at last it returns to zero at the beginning of the final SR phase.Hence, according to [80] we conclude that, the sharp transitions take place between three phases of the inflationary era.As regards the reduction in the ε 1 during the USR phase (see ε 1 in Fig. 1), according to Eq. ( 14) the amplification of P R can be presumable in this phase.Into the bargain, the value of ε 2 for each Case of the model oversteps its bound (slow-roll condition ε 2 ≪ 1 ) in USR region (see ε 2 in Fig. 1).Ergo, in the USR regime, the slow-roll conditions is violated by ε 2 even though it is held by ε 1 .Albeit, the validity of the slow-roll conditions in the vicinity of CMB scale (N * = 0), is obvious from the ε 1 and ε 2 graphs of Fig. 1.Consequently, it is allowable to compute the scalar spectral index n s and tensor-to-scalar ratio r, applying Eqs. ( 17) and ( 19) under the slow-roll approximation for the model.It is inferable from the listed quantities in Table II that, n s and r apropos of each Case of the model can be compatible with Planck 2018 data (TT,TE,EE+lowE+lensing+BK15+BAO, 95% CL) [53].Furthermore, the values of r apropos of all Cases are compatible with BICEP/Keck 2018 data (r < 0.036 at 95% CL) [73].Hence, the quartic potential through the generalized power-law non-canonical setup is As described erstwhile, in the USR domain, validity of the slow-roll conditions is contravened by ε 2 .Ergo, it is unallowable to compute the curvature power spectrum via Eq.( 14) governed by the slow-roll conditions.Thus, the evolution of the curvature perturbations should be evaluated through the subsequent Mukhanonv-Sasaki (MS) equation [72] all over the inflationary era therein the prime connotes derivative against the conformal time η ≡ a −1 dt, and the CMB constraints [53], PTA constraints [81], the impact on neutron-to-proton ration during the Big Bang Nucleosynthesis (BBN) [82][83][84], and the µ-distortion of CMB [85,86], respectively.
In order to solve the MS equation ( 23) numerically, it is necessary to take the ensuing Bunch-Davies vacuum state [10] into account as the inceptive condition on the scales deep interior the horizon Thenceforth, the curvature perturbations power spectrum is calculated from the subsequent equation The acquired quantities for the curvature power spectra P peak

IV. PBHS MASS SPECTRA
The current section is devoted to evaluate the mass spectra of produced PBHs through the quartic potential in non-canonical framework with the generalized power-law Lagrangian (2).As described in the preceding sections, the seeds of PBHs generate in the inflationary era from the enhanced scalar perturbations.The perturbed scales leave the horizon during the inflation and thereafter when they revert to the horizon in RD era, the ultra-dense districts could be produced.Eventually, PBHs can be born from the gravitational cave-in of the mentioned districts.The mass of PBHs could be defined as a portion of the horizon mass via the subsequent equation 3 connotes the collapse efficiency [4] and g * = 106.75implies the effective number of relativistic degrees of freedom.The fractional abundance f PBH (M) of PBHs specifies that what percentage of the cosmos DM content could be allotted to PBHs and can be obtained from the ensuing equation wherein Ω DM h 2 ≃ 0.12 implies the defined current DM density parameter by Planck data [53].Furthermore, β(M) indicates the generation rate of PBHs and it can be calculated from the Press-Schechter formalism for gaussian distribution of the primal perturbations as follows [87,88] Here "erfc" is the error function complementary and δ c = 0.4 indicates the density threshold [89,90].Moreover, σ 2 (M) denotes the coarse-grained density contrast smoothed on the scale k as in which P R is the curvature power spectrum and W (x) = exp (−x 2 /2) denotes the Gaussian window.
In the next, the PBHs abundance for each Case of Table I can be obtained from the Eqs.( 27)- (30), by substituting the attained P R from the MS equation ( 23) in Eq. ( 30).The consequent quantities for f PBH and M PBH (pertinent to peak location) have been classified in Table II for each Case of Table I.Additionally, the PBHs mass spectra thereto the observational domains have been schemed in Fig. 5 for all Cases of the model.As it can be seen from Table II, there is a huge difference in the PBH masses of the cases A, B and C. It can be inferred from Eq. ( 27) that, PBHs masses are proportional to the scales k wherein the PBHs are produced.Originally these scales are pertinent to the peak scales of the scalar power spectra.So we got the huge differences in the PBHs masses, because of the differences between the wavenumbers of the peaks of the scalar power spectra.The peak position of the scalar power spectrum defines the mass of produced PBH.The key parameter which is responsible to specify the peak position of scalar power spectrum is φ c for each case of the model.So the parameter φ c plays a key role in specifying the mass of produced PBH.
It can be inferable from Fig. 5 that, the consequent PBHs mass spectrum from the parameter Case A could be a promising claimant for 99% of DM content.The anticipated PBHs mass spectrum for the parameter Case B has placed in the permissible zone of the OGLE data [31,52], and it could be proper to be considered as the origin of ultrashorttimescale microlensing events.The foretold PBHs mass spectrum concerning the parameter Case C, has situated in the sensitivity zone of LIGO-VIRGO and its concurrent GWs would be traced by these observatories.The colory districts demonstrate the present observational restriction on the PBHs abundance.

V. SECONDARY GRAVITATIONAL WAVES
Propagated secondary GWs in the cosmos, could be the further upshot of reentry of the perturbed modes to the horizon in the RD era.Whereas, the concurrent GWs with PBHs production could be tracked down by multifarious detectors, so they can be considered as an indirect way to detect PBHs.The present section is earmarked for scrutinizing the secondary GWs coeval with PBHs production in the generalized power-law non-canonical model with the quartic potential.
It can be shown that the present energy density of the secondary GWs is given by [105] Ω therein Θ indicates the Heaviside theta function and η c implies termination time of the growth of Ω GW .The present GWs energy spectrum and its counterpart at η c are associated together through [81] Ω GW 0 h 2 = 0.83 g * 10.75 Thenceforth, the present density parameter spectra of the secondary gravitational waves Ω GW 0 concurrent with PBHs respecting the all Cases of Table I could be attained, through the Eqs.( 31)-( 33) and the obtained P R from the MS equation (23).
The consequent Ω GW 0 spectra have been plotted in Fig. 6 for Cases A (red line), B (green line), C (blue line), thereto sensitivity domains of GWs detectors like SKA (purple domain) [54], EPTA (brown domain) [55][56][57][58], LISA (orange domain) [59,60], BBO (green domain) [106,107] and DECIGO (red domain) [106,108].The predicted Ω GW 0 spectrum of parameter Case A, has lied in the sensitivity domain of LISA.Moreover, the resultant Ω GW 0 spectra for Cases B and C have situated in the sensitivity domain of SKA detector (see Fig. 6).Thus, the rectitude of this model could be verified in light of approaching observations of these detectors.
In the last stage, the slopes of Ω GW 0 spectra at different frequency regions have been estimated for all Cases of the model.It has been exhibited that, the inclination of Ω GW 0 spectrum could be matched with a power-law function of frequency as Ω GW 0 (f ) ∼ f n [23,109,110].The estimated quantities for the power indexes n for all Cases of the model in three frequency regions like f ≪ f c , f < f c and f > f c have been arranged in Table III (f c denotes the frequency of the peak of Ω GW 0 spectrum).It can be inferred that, the consequent quantities for n in the infrared region f ≪ f c could be compatible with the logarithmic equation n = 3 − 2/ ln(f c /f ) acquired in [111][112][113].As it can be seen from the listed results in Table III, there are large differences between the frequencies of GWs.The secondary GWs are produced concurrent with PBHs generation.
The frequencies of GWs are relevant to the wavenubers k through Eq. ( 33).As we mentioned previously, the key parameter to specify the scale k for each case is φ c .On the other word the frequencies of the peaks of the GWs spectra are pertinent to peaks positions of scalar power spectra, which are specified by φ c .

VI. CONCLUSIONS
This study is devoted to exhibit creation of PBHs from the quartic potential in noncanonical inflationary model with a generalized power-law lagrangian (2).The observational anticipations of the quartic potential on CMB scale, through the power-law non-canonical setup have been remedied.The consequent quantities for scalar spectral index n s and tensorto-scalar ratio r for all Cases of the model could place in the allowable domain defined by Planck 2018 (TT,TE,EE+lowE+lensing+BK15+BAO, 95% CL) [53].Additionally, the r quantities respecting all Cases of the model could be compatible with the latest limitation r < 0.036 of BICEP/Keck 2018 data at 95% CL [73] (notice Table II).
Choosing the non-canonical mass scale parameter M(φ) of the Lagrangian (2) as a twoparted peaked function of φ (20)-( 21) leads to make the inflaton slow down in an USR region (in the vicinity of the peak position φ = φ c ) with no effect on CMB scale.Thus, after adjusting the parameters of the model according to Table I, the sufficient enhancement in the scalar power spectrum during the USR region could be produced to born PBHs seeds.In this way not only the quartic potential could originate a feasible inflationary era, but also the PBHs and GWs seeds could be generated.Note that our model consists of 6 free parameters {α, M 0 , λ, ω, φ c , b}.Since λ = 0.13 is fixed, the parameter M 0 can be obtained from the CMB normalization constraint.So α is the only free parameter which can be fine tuned to remedy the results of n s and r for the quartic potential.Three parameters {ω, φ c , d} of the peak function do not have any significant effect on CMB scale.These three parameters should be fine tuned to produce PBHs in different mass range with allowable abundances.The ω parameter is responsible for the height of the peak of P R and subsequent PBH abundance, The φ c parameter is defined the peak position of the P R and the mass of produced PBH and d parameter is responsible for duration of USR stage and breadth of the peak of the P R .The USR stage in this model lasts for about two e-folds numbers.Following [80], from the behavior of ε 2 we inferred that, the sharp transitions could occur between three phases of the inflationary era of this model (first SR phase, intermediate USR phase and final SR phase).Sharp transition between the SR and USR phases may cause quantum corrections like large one-loop corrections in the curvature power spectrum [80,114,115].Moreover, the transition procedure could affect the size of local non-Gaussianity [80].Studying of possible quantum corrections and non-Gaussianity in this framework are beyond the scope of the present work and we leave them for future works.
The background evolution through the concurrent solving of Eqs. ( 8) and ( 13) could be evaluated.Thence, Fig. 1 embodies plots of scalar field φ and its derivative φ ,N , thereto the first ε 1 and second ε 2 slow-roll parameters against the e-fold number N. It is obvious from the mentioned plots that, in the USR region the velocity of the inflaton approaches to zero and ε 1 undergoes a decrease, which yields an increase in the curvature perturbations power spectrum (see Fig. 3).The slow-roll conditions ({ε 1 , ε 2 } ≪ 1) all over the USR region are obeyed by ε 1 , whereas broken by ε 2 momentarily.We also checked that the swampland criteria are satisfied in our model.
In addition, the curvature perturbations power spectra could be computed by way of solving the MS equation (23) for Cases of the model (notice Table II and Fig. 3).The portrayed spectra of P R at CMB scale are confined in the observational constraint of Planck 2018 (P * R ≃ 2.1 × 10 −9 ), while they enhance around seven order of magnitude in the USR region to produce PBHs.Three Cases of PBHs mass spectra according to parameter Cases of the model have been attained via the numerical solutions of MS equation.The consequent PBHs from the parameter Case A could be a promising claimant for 99% of DM content.The anticipated PBHs for the parameter Case B in the permissible zone of the OGLE data [31,52], could be the origin of the ultrashort-timescale microlensing events.The foretold PBHs concerning the parameter Case C, in the sensitivity zone of LIGO-VIRGO could be considered as the source of detected GWs by these detectors.The consequent values for PBHs abundance f PBH and mass M PBH are arranged in Table II and mapped in Fig. 5.
Thereafter, the spectra of Ω GW 0 (the current density parameter of coeval GWs with PBHs) have been obtained and plotted in Fig. 6 for all Cases.The predicted Ω GW 0 spectrum from parameter Case A, has lied in the sensitivity domain of LISA, whereas the resultant Ω GW 0 spectra for Cases B and C have situated in the sensitivity domain of SKA detector (see Fig. 6).Thus, the rectitude of this model could be verified in light of the approaching data of these detectors.
In the last stage, the power-law slopes of Ω GW 0 spectra (Ω GW 0 (f ) ∼ f n ) [23,109,110] at different frequency regions have been estimated for all Cases of the model.Table III embodies the peak frequencies of Ω GW 0 spectra and the estimated values for the power indexes n for all Cases of the model in three frequency regions like f ≪ f c , f < f c and f > f c .It can be inferred that, the consequent quantities for n in the infrared region f ≪ f c could be compatible with the logarithmic equation n = 3−2/ ln(f c /f ) acquired in [111][112][113].

FIG. 1 :
FIG. 1: Evolution of the scalar field φ and its derivative φ ,N as well as the first ε 1 and second ε 2 slow-roll parameters in regard of the e-folds number N for the Cases A (red lines), B (green lines) and C (blue lines).

FIG. 3 :
FIG. 3: The curvature perturbations power spectra computed through the numerical solutions of the MS equation (23) against the comoving wavenumber k respecting the Cases A (red line), B (green line) and C (blue line).The light-green, yellow, cyan, and orange sectors are excluded by

R 4 K/w ¡ - 1 PRFIG. 4 :
FIG.4:The relation between P R and non-canonical α parameter as to Case A of TableI.The amplitude of P R grows up for greater α in USR domain (see the different values of α in the legend of the plot).

FIG. 5 :
FIG. 5: The PBHs mass spectra relevant to Cases A (red line), B (green line) and C (blue line).

in which Ω r 0 h 2 ≃ 4 . 2 ×
10 −5 indicates the present radiation density parameter and g * ≃ 106.75 specifies the effective degrees of freedom in the energy density at η c .The relation between frequency and wavenumber is f = 1.546 × 10 −15 k Mpc −1 Hz.

FIG. 6 :
FIG. 6: The acquired present energy density spectra of secondary gravitational waves Ω GW 0 h 2 against the frequency respecting the parameter Cases A (red line), B (green line) and C (blue line) of Table I, in addition to the observational regions of the GWs detectors such as EPTA (brown domain), SKA (purple domain), LISA (orange domain), DECIGO (red domain) and BBO (green domain).The power-law incline of the Ω GW 0 h 2 spectrum has been schemed via dashed black lines in three frequency regions for Case A.
. Table II embodies the resultant quantities for n s , r and generated PBHs.

TABLE I :
The assigned parameters for Cases A, B, and C considering the fixed λ = 0.13 and

TABLE II :
[75,77,79]ed quantities for Cases of TableIwith regard to the scalar spectral index n s , the tensor-to-scalar ratio r, the peak of the amplitude of the scalar power spectrum P peak R , the wavenumber k peak , the PBHs abundances f peak PBH and masses M peak PBH .Here, the length of the inflationary era is ∆N = N end − N * and the observable quantities n s and r are calculated at CMB horizon passing e-folds number (N * = 0).To the best of our knowledge, deficiencies of the Hot Big Bang (HBB) theory could be rectified through a feasible inflationary epoch with a length of 60-70 e-folds number[75,77,79].The last column of TableIIembodies the length of inflation (∆N = N end − N * ) apropos to Cases of TableI, from the CMB horizon passing point (N * = 0) to the end of inflation (N end ).Whenever the first slow-roll parameter meets one (ε 1 = 1), the inflationary era ends

TABLE III :
The consequent values for frequencies and heights of the peaks of Ω GW 0 h 2 spectra, in addition to the power index n in frequency regions f ≪ f c , f < f c and f > f c respecting Cases A,