Inflation Based on the Tsallis Entropy

We study the inflationary scenario in the Tsallis entropy-based cosmology. The Friedmann equations in this setup can be derived by using the first law of thermodynamics. To derive the relations of the power spectra of the scalar and tensor perturbations in this setup, we reconstruct an $f(R)$ gravity model which is thermodynamically equivalent to our model in the slow-roll approximation. In this way, we find the inflationary observables, including the scalar spectral index and the tensor-to-scalar ratio in our scenario. Then, we investigate two different potentials in our scenario, including the quadratic potential and the potential associated with the natural inflation in which the inflaton is an axion or a pseudo-Nambu-Goldstone boson. We examine their observational viability in light of the Planck 2018 CMB data. We show that although the results of these potentials are in tension with the observations in the standard inflationary setting, their consistency with the observations can be significantly improved within the setup of the Tsallis entropy-based inflation. Moreover, we place constraints on the parameters of the considered inflationary models by using the current observational data.


I. INTRODUCTION
The first inflationary model was proposed by Starobinsky [1] in 1980 and it was based on the addition of the R 2 term in the Einstein-Hilbert action with the motivation to include semi-classical quantum effects to the gravity theory.It is interesting to point out that although this model is the first inflationary model, it is in very good agreement with the current observational data.Then, in 1981, Sato [2,3], suggested a scenario in which the Universe has undergone a rapid acceleration in the early stages of its evolution, and afterward it turned into a fireball with a very high temperature.Subsequently, in the same year, Guth [4] has shown that by including an accelerating phase before the radiation-dominated era, one can resolve the hot Big Bang cosmology problems such as the flatness problem, the horizon problem, and the magnetic monopole problem.Guth's model [4] is known as the old inflation and it is based on a scalar field that goes from a false vacuum towards a true vacuum through a first-order phase transition in the form of a quantum-tunneling process.
The old inflation [4] suffers from some fundamental problems, and to resolve its problems, other inflationary models were suggested later [5][6][7][8][9].For a nice review of the history of the first 30+ years of inflation, see [10].
Along with the evolution of the inflationary models, the theory of cosmological perturbations was developed too [11][12][13][14][15][16].The calculation of the primordial tensor perturbations in the early de Sitter stage was first elaborated by Starobinsky [11] in 1979, in terms of the related quantity which is the spectrum of primordial GW background after the Hubble radius crossing at the radiation-dominated stage.His original motivation was physically sound with the aim to investigate the initial state of the Universe.Then, Mukhanov and Chibisov [12] calculated the spectrum of the perturbations in the range of observable scales in the context of the Starobinsky R 2 inflation [1].In addition, the quantitatively correct expression for the GWs spectrum produced in the Starobinsky model [1] was first presented in [16].The quantum fluctuations of the scalar field during inflation lead to the generation of the perturbations whose growth can seed the Large-Scale Structure (LSS) formation and the anisotropies observed in the Cosmic Microwave Background (CMB) radiation [11][12][13][14][15][16].
Thus, from the observations corresponding to the LSS formation and CMB anisotropies, we can obtain valuable information about the physics of the early Universe.Important observational results about inflation have been presented by Planck collaboration [17], which are obtained from measurements of the CMB anisotropies in both temperature and polarization.Using these observational results, one can discriminate between different inflationary models.
The simplest inflationary model is based on a single scalar field minimally coupled to gravity [18,19].The scalar field responsible for inflation is known as the inflaton.If the potential energy of the inflaton dominates over its kinetic energy, then the inflaton rolls downward its potential slowly.In the chaotic inflationary scenarios [7], the slow-roll motion of the inflaton is provided by the Hubble friction term in the equation of motion of the scalar field, which is predominant in the early stages of the inflaton, and it then suppresses as the inflation goes towards its end.In the new inflation models [5], however, the slow-roll phase of inflation is provided by a plateau-like potential for the inflaton field.
Inflation has occurred in the GUT energy scale which is about two orders of magnitude less than the Planck energy scale.Since inflation has occurred in the regime of high-energy physics, therefore it is expected that quantum gravitational effects had a decisive role in the dynamics of the early Universe.In those energy scales, the gravity theory may be modified due to the quantum gravitational effects.One way to determine the corrections to the gravity theory is the use of the gravity-thermodynamics conjecture implying that there is a deep connection between gravity and thermodynamics [20][21][22][23][24][25][26][27][28].This conjecture implies that using the thermodynamics laws for the Universe as a thermodynamical system, one can derive the gravitational equations governing the evolution of the Universe.In particular, the Friedmann equation can be derived from the first law of thermodynamics [29][30][31][32][33][34].If the entropy is assumed as the Hawking-Bekenstein entropy [35], one can derive the standard Friedmann equation.If we consider some modifications to the entropy of the system, then the gravitational equations will be modified accordingly.In the early Universe, the entropy is expected to be different from the standard Bekenstein-Hawking entropy [35] due to the quantum gravitational effects that appear in the regime of high-energy physics.
One form that can be considered for the entropy of the Universe in its primordial stages is the Tsallis entropy S h ∝ A β [36] in which A is the horizon area and β is a constant parameter known as the Tsallis parameter.This entropy is a generalization of Boltzmann-Gibbs entropy [36,37], and it has been suggested to solve a thermodynamic puzzle.The Boltzmann-Gibbs statistics is not capable of describing the systems with divergent partition functions such as the gravitational systems and requires a non-additive generalization of the entropy definition [38][39][40].Accordingly, Tsallis and Cirto [36,37] introduced an entropy expression leading to non-extensive statistics.The Tsallis entropy has attracted a high level of research interest over the years, and so far, remarkable results have been found based on this entropy in a lot of complex systems such as self-gravitating stellar systems [41,42], black holes [36,37], background radiation [43], neutrinos [44,45], holographic dark energy [46][47][48][49][50] and dark matter [51], thermodynamic gravity [33,52,53], low-dimensional dissipative systems [39], and polymer chains [54].In the special case where the Tsallis parameter is taken as β = 1+∆/2 with 0 ≤ ∆ ≤ 1, the Tsallis entropy reduces to the Barrow entropy [55].
In this paper, we investigate the implications of the Tsallis entropy for the inflationary phase of the early Universe.To provide the accelerated expansion of the Universe in the inflationary phase, we assume the matter-energy content of the Universe follows the form of a canonical scalar field that plays the role of the inflaton.We take the inflaton potential to be in the form of a simple quadratic potential V (φ) = m 2 φ 2 /2 which leads to a chaotic inflationary scenario [7].In the framework of standard inflation, this potential is not favored in light of the current CMB data provided by the Planck 2018 collaboration [17].We show that the results of this potential for the inflationary observables can be improved significantly in the context of Tsallis entropy-based inflation.
One another elegant inflationary scenario that provides a sensible mechanism to generate a flat potential is natural inflation [68].In this class of models, the scalar field φ enjoys a shift symmetry φ → φ + const., which is slightly broken explicitly or due to the nonperturbative quantum effects to a discrete symmetry φ → φ + 2π [68].This feature gives rise to a periodic potential, as appropriate for inflation [68,69].Note that the presence of this symmetry protects the potential from radiative corrections.The scalar field with a flat potential originating from the shift symmetry is known as an axion.In this regard, the inflaton in the natural inflation is an axion or a pseudo-Nambu-Goldstone boson.The original natural inflation is described by a cosine-type periodic potential in the framework of a standard inflationary scenario based on the Einstein gravity where the entropy of the horizon obeys the Bekenstein-Hawking area law [35].In this setup, the natural inflation is not very favored by the latest observation of the Planck 2018 collaboration [17], and its results can satisfy only the 95% CL constraint of Planck 2018 data [17].This point motivates us to study the natural inflation in the setup of the Tsallis inflation and compare its predictions with the Planck 2018 constraints [17].
Although the Friedmann equations of the Friedmann-Robertson-Walker (FRW) Universe in the framework based on the Tsallis entropy have already been derived in [33], the study of inflation in this context also requires the equations of the power spectra of the primordial scalar and tensor perturbations.To this aim, we need the action of the gravity theory.In the absence of a unique action for our Tsallis inflation model, we construct an f (R) gravity model which is thermodynamically equivalent to the Tsallis gravity in the slow-roll regime of inflation.The modified f (R) gravity is a conceivable generalization of general relativity in which f (R) is an arbitrary function of the Ricci scalar R.Over the years, the f (R) gravity models have been extensively studied in the literature (see, e.g., [70][71][72][73][74][75][76][77][78][79][80]).One of the most important motivations that led us to this choice is that it is possible in the f (R) gravity to construct a model whose entropy is proportional to the powers of the horizon area, like the Tsallis entropy.Moreover, the f (R) gravity models have phenomenological effective features that can describe inflation in the regime of high-energy physics.On one side f (R) actions are general enough to cover some basic features of higher-order gravity, on the other side they are sufficiently simple to be easy to work with [71].Besides, there is one other important reason that convinces us that the f (R) theories of gravity are good candidates to help us evaluate the inflationary observables based on the Tsallis entropy.These models can avoid serious problems such as negative energies and related instabilities that are called the Ostrogradski instabilities [71,81,82].This feature made the f (R) theories unique compared to other higher-order gravity theories.
The paper is structured as follows.In Sec.II, we briefly study the thermodynamics of the Tsallis cosmology and the background equations derived from the first law of thermodynamics.In Sec.III, we review the inflationary dynamics and spectra of primordial perturbations in the f (R) gravity.The thermodynamic behavior of the field equations in the f (R) gravity will be discussed in Sec.IV.Then, in Sec.V, we reconstruct an f (R) model which is thermodynamically equivalent to the Tsallis model.Using this equivalence, we derive the equations of the inflationary observables in our Tsallis inflation scenario.In Sec.VI, we examine the observational compatibility of the quadratic and natural potentials with the Planck 2018 data in the framework of the Tsallis inflation.Finally, we summarize our concluding remarks in Sec.VII.

II. THERMODYNAMICS OF THE TSALLIS COSMOLOGY
In [33], it has been shown that by starting from the first law of equilibrium thermodynamics, dE = T h dS h + W dV, at the apparent horizon of the FRW Universe and taking the entropy associated with the apparent horizon in the form of the Tsallis entropy [36], one can derive the modified Friedmann equations.Here, E is the total energy content of the Universe, T h is the temperature of the apparent horizon, W is the work density, and V is the volume inside the apparent horizon.For these quantities, we have [33,83,84] where ρ and p represent respectively the energy density and pressure of all matter components in the Universe, and they satisfy the following continuity equation In the above equations, H ≡ ȧ/a denotes the Hubble parameter, rA indicates the radius of the apparent horizon, and the overdot indicates the derivative with respect to the cosmic time t.The horizon entropy is denoted by S h , and it will be taken in the form of the Tsallis entropy which is a non-extensive generalization of Boltzmann-Gibbs entropy [36][37][38]] Here, A = 4πr 2 A is the area of the apparent horizon and β is a real parameter known as the Tsallis parameter that measures the degree of non-extensivity [36].In addition, γ is an unknown constant, and several definitions have been presented for it in some literature for convenience [33,[85][86][87][88].In our investigation, however, this parameter will be treated as an unknown parameter.It is clear that for β = 1 and γ = 1/(4G), where G is the Newton gravitational constant, the Tsallis entropy (6) reduces to the Bekenstein-Hawking entropy [35].
On the flat FRW background, the modified Friedmann equations based on the Tsallis entropy (6) can be derived as follows [33] where the effective gravitational constant is defined as In this paper, we use the unit system in which c = = κ B = 1 where c is the speed of light, is the reduced Planck constant, and κ B is the Boltzmann constant.We define where m P is the Planck mass with a reduced value M P = m P / √ 8π = (8πG) −1/2 which throughout this paper we take it equal to unity, M P = 1.
From Eqs. ( 7), (8), and (9), one can easily show that the standard Friedmann equations based on the Einstein gravity in which the entropy associated with the horizon follows the form of Bekenstein-Hawking entropy [35], are recovered for β = 1 and γ = 1/(4G), as expected.It is worth mentioning that although the author of [33] has defined γ as γ ≡ (2 − β)(4π) 1−β /(4β G), we do not follow this convention in the present work.For a spatially flat Universe, we have rA = 1/H, and so one can rewrite the area of the apparent horizon A in terms of the Hubble parameter H as We assume that the matter-energy content of the Universe to be a scalar field φ in the form of a perfect fluid with the energy-momentum tensor T µ (φ) ν = diag (−ρ φ , p φ , p φ , p φ ).Here, ρ φ and p φ denote the energy density and pressure of the scalar field, respectively, and they are given by In these equations, V (φ) is the potential energy of the scalar field.The energy density ρ φ and pressure p φ of the scalar field fulfill the continuity equation From this equation together with Eqs.(11) and (12), the equation of motion of the scalar field equation will be obtained as where (, φ) indicates the partial derivative with respect to the scalar field φ.As we see, Eq. ( 14) is the same as the one that is valid in the standard scenario based on the Einstein gravity.
To study the inflationary epoch in the Tsallis entropy-based scenario, we need to derive the fundamental relations governing the theory of cosmological perturbations.For this purpose, we need the action of the model.But we don't know the action of the model, and to overcome this problem, we try to reconstruct an f (R) gravity model that is equivalent to our model from the thermodynamic point of view.In the f (R) gravity, the horizon entropy is given by [72,89] where F ≡ df (R)/dR.In the following, we review the basic formulas governing the theory of cosmological perturbations in the f (R) gravity, as well as the thermodynamic behavior of field equations in this theory.Then we apply the obtained results in our Tsallis entropybased model and find the necessary relations to study inflation.

III. INFLATIONARY DYNAMICS AND SPECTRA OF PRIMORDIAL PERTUR-BATIONS IN THE f (R) GRAVITY
The f (R) gravity is described by the following action [72] where g is the determinant of the metric g µν , f (R) is an arbitrary function of the Ricci scalar R, and X ≡ − 1 2 g µν ∂ µ φ∂ ν φ is the canonical kinetic term.For the flat FRW metric, the Ricci scalar R is given by [72] For the f (R) gravity model with action ( 16), the Friedmann equations turn into [72] 3F Following [70,72], we introduce the slow-roll parameters as where In the slow-roll regime, the quantities | ε i | (i = 1, 2, 3, 4) are much smaller than unity.
Under the slow-roll limit, the Ricci scalar R in Eq. ( 17) reduces to The spectrum of the curvature perturbations generated during inflation and in the slowroll limit can be estimated as [72] where It is worth mentioning that P s is computed at the time of horizon exit at which k = aH, where k is the comoving wavenumber.The observational value of the amplitude of scalar perturbations at the CMB pivot scale k * = 0.05 Mpc −1 has been constrained by the Planck 2018 CMB observations to be P s (k * ) ≃ 2.1 × 10 −9 [17].
The scalar spectral index n s during the slow-roll regime, in the framework of the f (R) gravity is given by [72] The observational constraint from the Planck data on the scalar spectral index at the TTTEEE+lowℓ+lowE) [17].
The tensor power spectrum in the framework of the f (R) gravity is given by [72] Also, the tensor spectral index n t in this setup is obtained as This parameter determines the scale dependence of the tensor power spectrum.Currently, there is no precise measurement for this quantity and we hope the future observations can provide some constraints on this observable.
Using Eqs. ( 23) and ( 26), one can find the tensor-to-scalar ratio in the f (R) gravity setting and in the slow-roll regime as [72] r The Planck 2018 data sets an upper bound on the tensor-to-scalar ratio as r < 0.0522 at the CMB pivot scale k * = 0.05 Mpc −1 (68% CL, Planck 2018 TTTEEE+lowℓ+lowE) [17].The most recent upper limit on this parameter is r 0.01 < 0.028 at 95% CL which is obtained at the CMB pivot scale k = 0.01 Mpc −1 using 10 datasets from the BICEP/Keck Array 2015 and 2018, Planck releases 3 and 4, and LIGO-Virgo-KAGRA Collaboration [92].
In the case of f (R) = R, using Eqs.( 24), ( 26), (27), and ( 28), we easily find that As we see, in this case, these results are reduced to the well-known relations in the standard inflation.

GRAVITY
Let us discuss the relation between the first law of thermodynamics and field equations in the f (R) gravity.In [72,89], it has been shown that the first law of equilibrium thermodynamics, dE = T h dS h + W dV, does not hold at the apparent horizon of the FRW Universe in the f (R) gravity.Consequently, to derive the Friedmann equations, we should apply the non-equilibrium thermodynamics and write the first law of thermodynamics as In this equation, S h is the entropy associated with the apparent horizon and it is still given by Eq. ( 15).Furthermore, S implies the non-equilibrium entropy and involves the non-equilibrium thermodynamic effects of the f (R) gravity [72,89].
From the non-equilibrium relation Taking the differentiation of Eq. ( 1), and then using Eq. ( 4) and also the relation dV = 4πr 2 A dr A , we obtain With the help of Eq. ( 5), we get It is assumed that the entropy of the horizon S h is in the form of Eq. ( 15).Differentiating the entropy (15), it follows that where we have used A = 4πr 2 A in deriving the above equation.Now, using Eqs.( 2), ( 3), ( 4), (31), and (32), in the right side of Eq. ( 29), we find Supposing the matter-energy content of the Universe to be a scalar field φ with the energy density ρ = ρ φ and the pressure p = p φ , we use Eqs.(11) and (12), and then we can rewrite Eq. ( 33) as Note that in deriving Eq. ( 34), we have used the relation rA = 1/H.Finally, with the help of the second and third relations in Eq. ( 20), we find From Eq. ( 29), it is clear that in the absence of the term d S/dS h the first-law of equilibrium thermodynamics dE = T h dS h + W dV on the apparent horizon holds.Therefore, Eq. ( 35) is an important relation in our examination.In the following, we show that in the framework based on the Tsallis entropy and under the slow-roll approximation, the term d S/dS h vanishes, and hence d S = 0. Therefore, this point allows us to use the obtained relations in the f (R) gravity to examine the inflationary models in the Tsallis entropy-based setting.

V. INFLATIONARY DYNAMICS IN THE TSALLIS ENTROPY-BASED MODEL
In this section, we assume that the Tsallis entropy-based model and the f (R) gravity are equivalent thermodynamically in the slow-roll regime.With the help of Eq. ( 22), one can easily express the horizon area A in Eq. ( 10) in terms of Ricci scalar R as Since we have supposed that the Tsallis entropy-based model and the f (R) gravity are equivalent in the slow-roll limit, we set Eqs. ( 6) and ( 15) equal to each other.Next, by using Eq. ( 36), we reach the following differential equation Since F (R) ≡ df (R)/dR, solving this differential equation analytically, the function f (R) is obtained as the following form As we see, f (R) is a power-law function of R as f (R) = µR n where µ ≡ (4 −1+2β (3π) −1+β γ G)/(2− β) and n ≡ 2 − β.For β = 1 and γ = 1/(4G), we find f (R) = R, and then Eq. ( 15) reduces to the Bekenstein-Hawking area-law of entropy [35] in the Einstein gravity.
Here, it is worthwhile to point out that the R n Lagrangian was originally regarded in the context of higher derivative theories [93,94], and then applied to inflation [95][96][97], which provides a simple and practical generalization of the Starobinsky R 2 inflation [1].In particular, in [97], this Lagrangian with n ≈ 2 has been investigated to establish a way to measure a deviation from the R 2 inflation [1].However, our methodology in the present paper differs from the approach of [97] in several aspects which are as follows.The action in our model is completely different from the action of [97] because, in our model, the gravitational part of the action consists of only the R n term, but the action of [97] contains the term R + R n which includes also the standard Einstein-Hilbert term.Furthermore, in our model, the contribution of a scalar field has been included in the action beside the R n , but such a contribution is absent in the action of [97].In the analysis of [97], a conformal transformation from the Jordan frame to the Einstein frame has been performed, and the calculations of the inflationary observables are accomplished in the Einstein frame.
In contrast, in our analysis, all calculations are performed in the Jordan frame.Finally, the form of the potentials that we consider in our work differs from the potentials that are regarded in [97] for the Einstein-frame scalar field.
Using Eq. ( 22), one can easily rewrite Eq. ( 38) as the following form From the slow-roll parameters (20), we find that the field equations ( 7) and ( 14) in the slow-roll limit lead to It can be shown that in the slow-roll approximation, the same equation as Eq. ( 40) can also be derived by using the first Friedmann equation ( 18) in the f (R) gravity.This arises from the fact that in the slow-roll regime, the non-equilibrium entropy can be neglected in front of the equilibrium entropy, as we will show explicitly at the end of this section.Applying Eq. ( 40), we can rewrite the function f in Eq. ( 39) in terms of the scalar field φ as Taking the time derivative of Eq. ( 40) and then applying Eq. ( 41), we can rewrite ε 1 in Eq.
Since d/dt = φd/dφ and d/dR = (1/R ,φ )d/dφ, we can also rewrite the slow-roll parameters ε 2 , ε 3 , and ε 4 in Eq. ( 20) in the following forms In deriving Eqs. ( 44), (45), and ( 46), we have also used Eqs.( 21) and (41).Besides, applying Eqs. ( 21), (24), and (41), the scalar power spectrum P s in Eq. ( 23) takes the form Substituting Eqs. ( 43), ( 44), (45), and (46) into Eq.( 25), the scalar spectral index n s reads We can also rewrite the tensor-to-scalar ratio r in Eq. ( 28) in terms of the scalar field φ as where we have used Eqs.( 21), (24), and (41) together with the relations d/dt = φd/dφ and It is convenient to evaluate the inflationary observations in terms of the so-called e-fold number N which measures the growth of the scale factor a during inflation.It is defined as where the subscript "e" refers to end of inflation.The definition (50) leads to Note that the anisotropies observed in the CMB exit the Hubble horizon around N * ≈ 50-60 e-folds before the end of inflation [99,100].The precise value of horizon exit e-fold number N * depends on the energy scale of inflation and also on the details of the reheating process after inflation [99,100].In our model, like most conventional inflationary models, the features of the reheating mechanism after inflation are unknown to us, and therefore it is not possible to determine the precise value of N * .
Using the last equality in Eq. ( 51) and also applying Eqs. ( 40) and ( 41), we reach the following differential equation One can solve Eq. ( 52) to obtain the scalar field φ as a function of the e-fold number N in the slow-roll approximation.In this way, we find the inflationary observables in terms of N.
Using Eqs. ( 40) and ( 41), we can rewrite Eq. ( 35) in the slow-roll limit as Using Eq. ( 22), we obtain R ,φ = 24HH ,φ , and since F = df /dR = f ,φ /R ,φ , we have ). Substituting this relation into Eq.( 53), we get Now substituting Eqs. ( 43) and (45) into Eq.( 54), and then using Eqs.( 9), (40), and (42), we will have which means that in the slow-roll Tsallis entropy-based inflation, the first law of equilibrium thermodynamics dE = T h dS h + W dV on the apparent horizon holds.This point makes us sure to apply the results derived in the f (R) gravity to find the inflationary observables in the Tsallis entropy-based cosmology.

VI. OBSERVATIONAL CONSTRAINTS
In this section, we apply the obtained results in the previous section to investigate the observational consistency of two different inflationary potentials in the framework of the Tsallis entropy-based inflation.These potentials are the quadratic and natural potentials that are not in good agreement with the current CMB data in the setting of standard inflation.

A. Quadratic potential
Let us continue studying the inflationary scenario in the Tsallis entropy-based setting by considering the quadratic potential [7] V where m is the inflaton mass.Applying the potential (56) and also using Eqs.( 9), ( 22), (40), and ( 42) , the slow-roll parameters in Eqs. ( 43), ( 44), (45), and ( 46) take the form The parameter Θ in these equations is defined as It is easy to show that for β = 1, we have Θ = γ.Therefore, for β = 1 and Θ = 1/(4G), we find the same relations in the Einstein gravity.
With the help of Eq. ( 69), we can rewrite Eq. ( 64) in terms of the e-fold number N and the parameters β, Θ, and m as Applying Eq. ( 69) in Eq. ( 65), we take Now, we can substitute Eq. ( 69) into Eqs.( 57), (58), and ( 59), and obtain the slow-roll parameters as the following forms Using Eq. ( 69) into Eq.( 61), we find the scalar power spectrum as As we see from Eq. ( 75), the power spectrum of the curvature perturbation is a function of the e-fold number N and three free parameters β, Θ, and m.One can use Eq. ( 75) and then impose the CMB normalization at the observable scale to find a constraint on the parameter m.
Applying Eq. ( 69), the scalar spectral index (62) and the tensor-to-scalar ratio (63) take the forms In the limit of standard inflation where β = 1, Eqs. ( 76) and ( 77) reduces to n s = 1 − 4/(2N + 1) and r = 16/(2N + 1), respectively.These results are the same equations that we find in the setup of the standard inflation.Let us now consider the limits where β ≪ 1 and N ≫ 1.From Eqs. ( 76) and ( 77), the leading contributions to n s and r become This means that in the regime β ≪ 1 and N ≫ 1, our theory reduces to the R 2 inflationary model proposed by Starobinsky [1].In this regime, the power spectrum (75) reduces to the form Using Eq. ( 80) and then impose the CMB normalization at the pivot scale k * = 0.05 Mpc −1 [17], we find In Fig. 1, the parameter β has been taken as a varying parameter in the range 0 < β < 2.
It should be noted that each value of β in this range is related to a special case of our powerlaw f (R) scenario, and in the case with β = 1 and γ = 1/(4G), our f (R) model reduces to the Einstein general relativity (GR).This does not mean that a transition from f (R) to GR has occurred during inflation in our scenario.In other words, to investigate each case of our scenario, we should fix the value of β at the first step, and it is not the case that this parameter varies during inflation and causes a transition from f (R) to GR.Since we cannot determine the parameter γ in our investigation, therefore we cannot determine the time at which the f (R) gravity in Eq. ( 38) transits to GR.
Fig. 1 shows that the prediction of the potential m 2 φ 2 /2 in the standard setting is not in good consistency with the Planck 2018 observations [17], while in the framework of the Tsallis entropy-based inflationary scenario, it can be in very good agreement with these data, and its results can lie inside the 68% CL region of Planck 2018 data [17].From the figure, we see that for small β, the model shows better consistency with the observations and the prediction of the model can enter the 68% CL region of these data.As we have proved, in the limit β ≪ 1, the inflationary observables n s and r approach to the same values in the Starobinsky R 2 inflation [1].
With the help of Eqs. ( 76) and ( 77), and also the Planck observational constraints on the r-n s plane, we can estimate the ranges of the parameter β for which the results of the model in the r − n s plane are consistent with the 68% CL region of the Planck 2018 data [17].In the case N * = 50, if 0 < β 0.045 the result of the model is in agreement with the 68% CL constraint of these data, and for N * = 60 the prediction of our model can enter the 68% CL region of the Planck 2018 data [17], provided that 0 < β 0.011.
Here, we are interested in applying the recent constraint of [92] on r 0.01 to present some observational constraint on the model parameter β.For this purpose, we should determine the e-fold number at which the comoving wavenumber k = 0.01 Mpc −1 exits the Hubble horizon during inflation.To do so, we examine the behavior of the comoving wavenumber k as a function of the e-fold number N, at which the mode with comoving wavenumber k leaves the Hubble horizon, k = aH.With the help of this relation, we can easily find which a(N * ) and H(N * ) are the scale factor and the Hubble parameter at the time of horizon exit of the mode k * = 0.05 Mpc −1 , respectively.Using Eq. ( 50), we can obtain the scale factor a in terms of the e-fold number N as where we have normalized the scale factor to its value at the epoch of horizon crossing of the mode k * = 0.05 Mpc −1 .Finally, applying Eqs. ( 71) and ( 84) in Eq. ( 83), the comoving number k can be found as a function of N as Using Eq. ( 85) and setting N * = 60, we plot in Fig. 2 the variation of the comoving wavenumber k against the e-fold number N. In the figure, we have also specified the comoving wavenumber k = 0.01 Mpc −1 and its corresponding e-fold number which is N 0.01 ≃ 61.62.
In Fig. 3, with the help of Eq. ( 77), the variation of the tensor-to-scalar ratio r 0.01 is plotted as a function of the parameter β by taking N 0.01 ≃ 61.62.The gray shaded region is excluded by the constraint reported by Galloni et al. [92] at the CMB pivot scale Here, it should be noted that throughout this paper, we have worked in the framework of the f (R) gravity given by Eq. ( 38), and we have treated GR as a special case of our f (R) gravity model that is realized by taking β = 1 and γ = 1/(4G).To clarify this point, we note that for a given set of model parameters (β, γ), one can find Θ from Eq. ( 60 function and causes the deviation from GR.Since γ is a free parameter in our model, we may take its value such that the regime of the f (R) gravity dominates over the GR regime throughout inflation, and consequently any transition would not occur during inflation at all.This means that in the effective action which may contain the f (R) term together with the Einstein-Hilbert term R, the contribution of the latter will be negligible compared to the former contribution, and therefore the effective action will be reduced to the action (16) which is used in the present work.Besides, since our investigation is not able to determine the precise value of γ, it is not possible to specify the time of the transition from f (R) to GR that may happen during inflation or in the post-inflationary Universe, in the present work.However, if future studies provide some understanding for us about the nature of the inflaton field, we may determine its effective mass, and accordingly estimate value of the γ.This enables us to estimate the time of such a transition.In addition, we may provide some observational constraints of our model parameter β in the setup of the Tsallis inflation, like the analysis performed in [63] for the case of the Barrow cosmology.

B. Natural inflation
The next model that we consider in our investigation is the natural inflation model in which the inflaton field is presumed to be an axion or pseudo-Nambu-Goldstone boson with a cosine-type periodic potential [68,101] where Λ is some non-perturbatively generated scale and φ 0 is the scale that determines the curvature of the potential.Both of these constants have dimensions of mass.It seems that the super-Planckian value of the scale φ 0 , i.e. φ 0 M P , is impossible in the context of string theory because all known controlled string theory constructions are restricted to φ 0 < M P [69,102].In the framework of standard slow-roll inflation, the prediction of the natural potential in the r − n s plane is in tension with the latest observations [17], in the sense that its prediction can lie within the 95% CL region of these data only for N * = 60.
This issue motivates us to study natural inflation in the Tsallis entropy-based scenario to see whether the novel framework can improve the consistency of the model with the Planck 2018 observations [17].
With the help of Eqs. ( 9) and ( 86), we can rewrite Eq. ( 42) as Substituting Eqs. ( 9) and ( 86) in Eq. ( 40), the Hubble parameter H is found as the following form We can apply the above equation together with Eq. ( 84) in Eq. ( 83), to find the comoving number k as To find the scalar field χ in terms of the e-fold number N, we perform the same steps as the previous subsection.In the first step, using Eq. ( 87) for a given set of the model parameters (β, η), we solve the equation ε 1 (χ e ) = 1 numerically to find χ e .After that, with the help of the obtained value of χ e , we solve differential equation ( 94), numerically, and consequently we find χ = χ(N).Substituting χ(N) into Eqs.( 87), ( 88), ( 89), ( 90), ( 95), (96), and ( 99), the slow-roll parameters ε 1 , ε 2 , ε 3 , ε 4 , the inflationary observables n s and r, and also the comoving number k are derived in terms of N, respectively.
In Fig. 6 Using the Planck observational constraints on r-n s plane and also Eqs. ( 95) and ( 96), we also can estimate the ranges of parameter η for which the prediction of our model for To present some constraint on the model parameter η by using the recent constraint of [92] on r 0.01 , we determine the e-fold number at which the comoving wavenumber k = 0.01 Mpc −1 The marginalized joint 68% and 95% CL regions of the Planck 2018 TTTEEE+lowℓ+lowE data [17] are specified by dark and light blue, respectively.
exits the Hubble horizon during inflation, with the help of Eq. ( 99) by taking β = 1/2 and N * = 60.It is found that N 0.01 ≃ 61.6.
In Fig. 7, using Eq. ( 96), we plot the variation of the tensor-to-scalar ratio r 0.01 against the parameter η by setting β = 1/2 and N 0.01 ≃ 61.6.The gray-shaded region is excluded by the constraint on the upper bound on r 0.01 reported by Galloni et al. [92].The prediction of the model is in agreement with this constraint if η 0.593.
The evolution of the slow-roll parameters ε 1 , ε 2 , ε 3 , and ε 4 versus the e-fold number N for the natural potential ( 86) is plotted in Fig. 8, by using Eqs.( 87 In Fig. 9, with the help of Eq. ( 97), the evolution of the function f (R) as a function of

VII. CONCLUSIONS
Inflation has occurred in the regime of high energy physics at which the gravitational theory is expected to be modified.Therefore, the entropy-area relation may undergo some modifications in those energy scales.This motivated us to regard the entropy of the early Universe to be in the form of the Tsallis entropy which is a generalization for the Bekenstein- Hawking entropy [35] and possesses the non-additivity and non-extensivity property.This form of entropy has a relation with the horizon area as S h = γA β , in which A is the area of the horizon and β and γ are unknown constants.We have studied the inflationary era in the context of the Tsallis entropy-based cosmology.Since there is no definite action for this setup, it was not possible to derive the power spectra of the primordial scalar and tensor perturbations in our Tsallis inflation scenario.To resolve this issue, we reconstructed an f (R) model which is thermodynamically equivalent to our setting in the slow-roll approximation.
This equivalence allows us to use the equations of the scalar and tensor power spectra  Moreover, we have proved that in the limits β ≪ 1 and N ≫ 1, the behavior of the quadratic potential in the r−n s plane coincides exactly with the prediction of the Starobinsky R 2 inflation [1].Furthermore, we have estimated that for N * = 50, the prediction of the model is compatible with the 68% CL constraint of the Planck 2018 observations, if 0 < β 0.045.In the case of N * = 60, the model is consistent with the 68% CL region of the observational data provided that 0 < β 0.011.The recent constraint r 0.01 < 0.028 (95% CL) [92] on the tensor-to-scalar ratio at the scale k = 0.01 Mpc −1 constrains this parameter to 0 < β 0.202.
We further examined the viability of the natural potential (86) in the inflationary setting based on the Tsallis entropy.In the framework of the standard inflation, the prediction of this potential in the r − n s plane is not very preferred according to the current CMB observations, regarding the fact that its results can satisfy only the 95% CL constraints of the Planck 2018 data [17].In this case, we have found the inflationary observables n s and r in terms of the model parameters β and η numerically.We have focused on the cases β = 1/2 and drawn the r − n s diagram for N * = 50, 60, with varying η in the range of η > 0. We have demonstrated that in the framework of the Tsallis entropy-based inflation, the natural potential provides a great fit to the Planck 2018 data [17].The prediction of the model lies inside the 68% CL region of these data.Furthermore, we presented some observational constraints on the model parameters by using the Planck 2018 data [17].Our results imply that for N * = 50, 60, the results of our model can lie inside the 68% CL region of the Planck 2018 data [17], provided that 0.172 η 0.433 and 0.281 η 0.613, respectively.The observational bound on r 0.01 provided in [92], also gives rise to the condition η 0.593 for the model parameter.

m ≃ 6 .
614 × 10 −7 Θ for N * = 60 .(82) As we see the value of the parameter m depends on what value the parameter Θ takes.With the help of Eqs.(76) and (77), we can plot the r − n s diagram and compare the prediction of our model with the Planck 2018 CMB data [17].Fig. 1 shows the prediction of our model in the r − n s plane for two typical values of N * and varying β in the range of 0 < β < 2. The dashed and solid black curves illustrate the results of the model for N * = 50 and N * = 60, respectively.Besides, The red solid line between the dashed and solid black curves shows the prediction of the potential in the standard inflation which corresponds to β = 1 in our model, in the range of 50 ≤ N * ≤ 60.Moreover, the prediction of the Starobinsky R 2 inflationary model [1] has been shown by the orange solid line with 50 ≤ N * ≤ 60.

FIG. 1 . 2 inflation [ 1 ]
FIG. 1.The r − n s diagram for the quadratic potential (56) in the slow-roll inflation based on the Tsallis entropy for two different values of N * with varying β in the range of 0 < β < 2. The results for N * = 50 and N * = 60 are shown by the dashed and solid black curves, respectively.The red solid line between the dashed and solid black curves shows the results of the potential in the standard inflation, in the range of 50 ≤ N * ≤ 60.Moreover, the result of the Starobinsky R 2 inflation [1] is specified by the orange solid line with 50 ≤ N * ≤ 60.The marginalized joint 68% and 95% CL regions of the Planck 2018 TTTEEE+lowℓ+lowE data [17] are specified by dark and light blue, respectively.

1 ]FIG. 2 .
FIG. 2. Evolution of the comoving wavenumber k versus the e-fold number N for the quadratic potential (56) in the Tsallis inflationary setting.The pink horizontal and vertical dashed lines specify k = 0.01 Mpc −1 and its corresponding e-fold number N 0.01 ≃ 61.62, respectively.

FIG. 3 .
FIG.3.Variation of the tensor-to-scalar ratio r 0.01 versus the parameter β by setting N 0.01 ≃ 61.62, for the quadratic potential(56) in the framework of the Tsallis inflation.The gray-shaded region is excluded by the constraint on the upper bound on r 0.01 , reported by Galloni et al.[92].The pink vertical dashed line specifies β ≃ 0.202 which is the maximum value of the parameter β for which the model satisfies the constraint r 0.01 < 0.028.

FIG. 5 .
FIG. 5. Evolution of the function f (R) versus the e-fold number N for the quadratic potential (56) in the Tsallis inflationary scenario with β = 0.01 and some typical values of γ.Also, the red curve illustrates the Ricci scalar R predicted by general relativity (GR).
Fig.6clears that the result of the potential(86) in the standard inflation for N * = 50 lies completely outside the allowed regions of the Planck 2018 observations[17] and only for N * = 60 enters the 95% CL region of these data.In the Tsallis entropy-based inflationary setting, however, its prediction can lie well within the 68% CL region allowed by the Planck 2018 data[17].

FIG. 6 .
FIG. 6.The r−n s diagram of the natural potential (86) in the Tsallis inflationary setting by taking β = 1/2 for two different values of N * with varying η in the range of η > 0. The results of our model for N * = 50 and N * = 60 are shown by the dashed and solid black curves, respectively.The red dashed (N * = 50) and solid (N * = 60) curves show the prediction of the standard natural inflation.
)-(90).In this figure, we have considered β = 1/2 and η = 0.65.The figure shows that the slow-roll approximation holds during inflation.This point verifies the viability of our investigation which is based on the slow-roll approximation.Moreover, from the figure, we see that the first slow-roll parameter ε 1 reaches unity at the end of inflation with N = 0.

FIG. 7 .
FIG.7.Variation of the tensor-to-scalar ratio r 0.01 against the parameter η for the natural potential(86) in the Tsallis inflationary scenario with β = 1/2 and N 0.01 ≃ 61.6.The gray-shaded region is excluded by the constraint on the upper bound on r 0.01 , reported by Galloni et al.[92].The pink vertical dashed line specifies η ≃ 0.593 which is the minimum value of the parameter η for which the model agrees with the constraint r 0.01 < 0.028.

FIG. 9 .
FIG. 9. Evolution of the function f (R) against of the e-fold number N for the natural potential (86) in the Tsallis inflation with β = 1/2 and η = 0.65, for some typical values of γ.The prediction of general relativity (GR) for the Ricci scalar R is also shown by the red color.
obtained in the f (R) gravity for our Tsallis inflation model.We have considered two different inflationary potentials in our scenario and checked their viability with the Planck 2018 observations[17].First, we studied the observational consistency of the quadratic potential (56) which provides a chaotic inflation model.In the standard inflationary setting based on the Einstein gravity, the consistency of this potential is not favored by the Planck 2018 observational data[17].This motivated us to investigate whether this potential can be resurrected in light of the Planck 2018 results in the setting of the Tsallis entropy-based inflation.We have derived analytic formulas for n s and r in terms of the parameter β and the e-fold number N, and then plotted the r − n s diagram for N * = 50, 60, with varying the parameter β in the range of 0 < β < 2. Our results imply that this potential can be in excellent consistency with the Planck 2018 data in the Tsallis entropy-based scenario, such that its results can lie inside the 68% CL region of the observational data.
[17]d different values of N * is compatible with 68% CL constraint of the Planck 2018 data[17].In the case N * = 50, the model is compatible with 68% CL constraint of the Planck 2018 data for 0.172 η 0.433.We also find that for N * = 60, the results of our model can lie inside the 68% CL region of the Planck 2018 data, if 0.281 η 0.613.