## Embedding Warm Natural Inflation in f ϕ T ( ) gravity

Sabina Yeasmin ∗ , Atri Deshamukhya †

Department of Physics, Assam University, Silchar, Assam 788011, India

We study warm inflation in the framework of f ( ϕ T ) gravity, where ϕ is the inflaton and T is the trace of the energy-momentum tensor. The inflaton field is assumed to roll on the natural potential and the result is analyzed in light of Planck 2018 and BICEP/Keck 2021 data. We start our work by obtaining the field equations under slow-roll approximations. We then evaluate the scalar and tensor power spectra and their corresponding spectral index and tensor-to-scalar ratio with a temperaturedependent form of the dissipation coefficient during the inflationary era. We find that the warm inflation model in f ( ϕ T ) gravity is compatible with observational bands.

## 1 Introduction

According to the Standard Model of Cosmology and observational data, the universe has undergone two phases of accelerated expansion throughout its evolution, the inflationary phase at its early stage and a second phase of accelerated expansion in recent times. This second phase of accelerated cosmic expansion is believed to be a consequence of the presence of dark energy exerting negative pressure. General Relativity is not ideal for explaining the existence of the dark sector i.e. dark matter and dark energy of the universe [1, 2, 3]. To overcome such shortcomings researchers started studying the different higher-order extensions to General Relativity and exploring alternative theories of gravity [4, 5, 6, 7, 8]. Modified theories of gravity not only address the late-time acceleration but these theories are also applied to explain the early Universe, e.g., the inflationary era.

In modified gravity theories, the Einstein-Hilbert action is modified either by extending the geometric or the matter part or both. Unlike Einstein's theory, whose field equations contain only up to second-order derivatives, the modified theories with higher derivative Ricci/Riemann tensor gravity models include higher derivatives [9]. A significant difference is therefore expected between general relativity and modified gravity predictions.

The f ( R ) gravity model is the first proposed modification of the Einstein gravity where the curvature scalar in the Einstein-Hilbert action has been replaced by a general function of Ricci scalar R . This theory describes both inflation and late-time expansion of the universe. The first viable f ( R ) model to study inflation was proposed by Nojiri and his collaborators [10, 11, 12, 13, 14]. T. Harko et al. [15] introduced another modified gravity theory called f ( R,T ) theory where the gravitational field is assumed to be coupled to the trace T of the energy-momentum tensor of the matter in the action. Later several inflationary models have been studied within the context of the f ( R,T ) gravity theory [16, 17, 18, 19, 20, 21]. There are many other technically allowed modified gravity theories that have been investigated in literature in the context of inflationary cosmology like f ( T ) [22], f ( R,G ) [23], Gauss-Bonnet (GB) theory ( f ( G ) gravity) [24], Brans-Dicke theory [25], etc. where T denotes the Torsion scalar and G is Gauss-Bonnet scalar.

Another recent modification is the f ( ϕ T ) theory of gravity, proposed by Zhang et al. [26], which is an extension of the simplest f ( R,T ) gravity theory. This modified theory of gravity has been investigated by adding a coupling term between the scalar field ϕ and the trace of the energy-momentum tensor into the action in general relativity. They studied cold inflationary dynamics in this framework by considering three well-motivated potentials: chaotic, natural, and Starobinsky, and found that chaotic and natural inflationary

∗ Email:sabyes22@gmail.com

† Email: atri.deshamukhya@aus.ac.in

models are in better agreement with the observational constraints. However, a larger value of the tensor-toscalar ratio was obtained for the Starobinsky model. Although natural inflation in f ( ϕ T ) gravity is consistent with the observational data, a major problem with that is a larger axion decay constant is required to agree with the observational constraints. So, it is difficult to embed this model in fundamental theories like string theory [27]. In this study, we show that it is possible to lower the value of the axion decay constant below the Planck scale if the natural inflationary model in f ( ϕ T ) theory of gravity is studied in the context of the warm inflationary scenario.

In 1995, Berera [28] introduced the idea of warm inflation as an alternative approach to resolving the graceful exit problem of the standard inflation dynamics. The basic concept of warm inflation is that during inflationary period, the inflaton field does interact with other existing fields [29, 30, 31, 32, 33, 34] which leads to vacuum energy decay to the other lighter fields, thereby creating particles simultaneously with the process of inflation. Thus in warm inflation, radiation is produced simultaneously with the expansion of the universe and there can be a smooth transition from the inflationary phase to the radiation-dominated phase without the requirement of invoking any separate reheating phase. In this scenario, dissipative effects play an important role during inflation.

In warm inflation scenario, energy dissipates from the inflaton field into radiation during the inflationary phase and this mechanism is parameterized by a term Γ called dissipation coefficient [31, 35, 36]. Because of this dissipative coupling between inflaton and radiation field inflation can last long even if the potential is not very flat. In this work, we use cubic temperature-dependent dissipation coefficient which affects the prediction of cosmological observables like it suppresses the value of the tensor-to-scalar ratio. The latest Cosmic Microwave Background (CMB) observations of BICEP/Keck [37] combined with the Planck 2018 data [38] sets an bound on the tensor-to-scalar ratio of r &lt; 0 036 . at the 95% confidence level. We show that the warm natural inflation in f ( ϕ T ) gravity model is consistent with this observational bound on r .

The paper is organized as follows: In section 2, we present a brief overview of f ( ϕ T ) gravity in the FLRW background. In section 3, we investigate warm inflation in f ( ϕ T ) gravity and study the cosmological perturbations originating from f ( ϕ T ) gravity during warm inflation and also formulate the corresponding power spectrum, the tilt of spectral index, and tensor-to-scalar ratio. In section 4, we study warm inflation with Natural potential in weak and strong dissipative regimes for variable dissipation coefficient. In section 5, we present our conclusions.

## 2 Overview of f ( ϕ T ) gravity

The action in f ( ϕ T ) gravity theory as proposed in [26] is given by,

$$S = \int \left [ \frac { R } { 2 \kappa } + \lambda f ( \phi ) T + L _ { m } \right ] \sqrt { - g } d ^ { 4 } x$$

where a coupling term (coupling between inflaton ( ϕ ) and trace of the energy-momentum tensor ( T ) is added to the Hilbert-Einstein action. Here, κ = 8 πG = 1 /M 2 P , M P is the reduced Planck mass, R stands for Ricci scalar, L m is the matter Lagrangian density, g is the determinant of the metric tensor, λ is the model parameter, and f ( ϕ ) is a dimensionless function of the inflaton ϕ which satisfies the condition f (0) = 0 . This indicates that when inflaton field decays and its energy no longer dominates the universe, Einstein gravity is recovered. Also in the limit λ → 0 , the action in equation (1) returns to the action in the Einstein theory. √

In this work, we have considered the form of f ( ϕ ) as f ( ϕ ) = κϕ . Also, we have used the natural units such that c = ℏ = 1 and used the ( -, + + +) , , sign convention for the metric tensor throughout this manuscript. With the form of f ( ϕ ) = √ κϕ , the action reads

$$S = \int \left [ \frac { R } { 2 \kappa } + \lambda \sqrt { \kappa } \phi T + L _ { m } \right ] \sqrt { - g } d ^ { 4 } x$$

By varying the action (2) with respect to the metric, we get the modified Einstein field equation in the following form,

$$R _ { \alpha \beta } - \frac { 1 } { 2 } g _ { \alpha \beta } R = \kappa T _ { \alpha \beta } ^ { ( e f f ) }$$

where T ( eff ) αβ is the effective stress-energy tensor given by,

$$T _ { \alpha \beta } ^ { ( e f f ) } = T _ { \alpha \beta } - 2 \lambda \sqrt { \kappa } \phi ( T _ { \alpha \beta } - \frac { 1 } { 2 } T g _ { \alpha \beta } + \Theta _ { \alpha \beta } )$$

where T is the trace of the energy-momentum tensor T αβ in Einstein gravity. In flat FLRW background T αβ is given by

$$T _ { \alpha \beta } = \partial _ { \alpha } \phi \partial _ { \beta } \phi + g _ { \alpha \beta } \left ( \frac { 1 } { 2 } \dot { \phi } ^ { 2 } - V ( \phi ) \right )$$

where V ( ϕ ) is the potential of the inflaton. For inflaton field, Θ αβ reads as

$$\Theta _ { \alpha \beta } = - \partial _ { \alpha } \phi \partial _ { \beta } \phi - T _ { \alpha \beta }$$

From equations (4),(5) and (6), the effective energy density and pressure for inflaton field are obtained as

$$T _ { 0 0 } ^ { ( e f f ) } = \rho _ { \phi } ^ { ( e f f ) } = ( 1 + 2 \lambda \sqrt { \kappa } \phi ) \frac { \dot { \phi } ^ { 2 } } { 2 } + ( 1 + 4 \lambda \sqrt { \kappa } \phi ) V$$

$$T _ { i j } ^ { ( e f f ) } = P _ { \phi } ^ { ( e f f ) } g _ { i j } = \left [ ( 1 + 2 \lambda \sqrt { \kappa } \phi ) \frac { \dot { \phi } ^ { 2 } } { 2 } - ( 1 + 4 \lambda \sqrt { \kappa } \phi ) V \right ] g _ { i j }$$

Substituting equation (7) into the 00 component of the modified Einstein field equation (3), we get

$$H ^ { 2 } = \frac { \kappa } { 3 } \left [ \frac { \dot { \phi } ^ { 2 } } { 2 } ( 1 + 2 \lambda \sqrt { \kappa } \phi ) + V ( 1 + 4 \lambda \sqrt { \kappa } \phi ) \right ]$$

which is known as the modified Friedmann equation. Here we have defined the Hubble parameter as H = ˙ a a . Similarly, substituting equation (8) into the ij component of the modified Einstein field equation (3), we get

$$\frac { \ddot { a } } { a } = - \frac { \kappa } { 3 } \left [ \dot { \phi } ^ { 2 } ( 1 + 2 \lambda \sqrt { \kappa } \phi ) - V ( 1 + 4 \lambda \sqrt { \kappa } \phi ) \right ]$$

This is the modified acceleration equation. We can also obtain the expression for ˙ H as,

$$\dot { H } = \frac { \ddot { a } } { a } - H ^ { 2 } = - \frac { \kappa } { 2 } ( p ^ { e f f } + \rho ^ { e f f } ) = - \frac { \kappa } { 2 } \left [ \dot { \phi } ^ { 2 } ( 1 + 2 \lambda \sqrt { \kappa } \phi ) \right ]$$

Equations (10) and (11) are known as the modified Friedmann second equation. In the limit λ = 0 , these equations reduce to those in Einstein gravity.

Also the continuity equation or the modified Klein-Gordon equation can be written as,

$$\ddot { \phi } ( 1 + 2 \lambda \sqrt { \kappa } \phi ) + 3 H \dot { \phi } ( 1 + 2 \lambda \sqrt { \kappa } \phi ) + \lambda \sqrt { \kappa } \dot { \phi } ^ { 2 } + ( 1 + 4 \lambda \sqrt { \kappa } \phi ) V ^ { \prime } + 4 \lambda \sqrt { \kappa } V = 0$$

where the prime denotes derivative with respect to the inflaton field.

## 3 Warm Inflation in f ( ϕ T ) gravity

In the warm inflationary paradigm, the inflaton field interacts with other fields during the inflationary period, leading to the decay of the inflaton field into radiation. To capture this feature, a dissipation co-efficient Γ is added to the equation of motion describing the warm inflationary scenario. The equations that completely describe the dynamics of the warm inflation scenario in f ( ϕ T ) gravity can be described as

$$H ^ { 2 } = \frac { \kappa } { 3 } ( \rho _ { \phi } ^ { ( e f f ) } + \rho _ { \gamma } )$$

$$\dot { \rho } _ { \phi } ^ { ( e f f ) } + 3 H ( \rho _ { \phi } ^ { ( e f f ) } + P _ { \phi } ^ { ( e f f ) } ) = - \Gamma \dot { \phi } ^ { 2 }$$

$$\dot { \rho } _ { \gamma } + 4 H \rho _ { \gamma } = \Gamma \dot { \phi } ^ { 2 }$$

Here ρ ( eff ) ϕ , ρ γ and P ( eff ) ϕ are the energy density of the scalar field, energy density of the radiation field and pressure of the scalar field respectively. Γ is the dissipation coefficient which describes the decay of inflaton into radiation during inflationary phase.

Modified equation of motion of the inflaton field in warm inflation can be derived by substituting equations (7) and (8) into equation (14),

$$\ddot { \phi } ( 1 + 2 \lambda \sqrt { \kappa } \phi ) + 3 H \dot { \phi } ( 1 + 2 \lambda \sqrt { \kappa } \phi + Q ) + \lambda \sqrt { \kappa } \dot { \phi } ^ { 2 } + ( 1 + 4 \lambda \sqrt { \kappa } \phi ) V ^ { \prime } + 4 \lambda \sqrt { \kappa } V = 0 \quad \ ( 1 6 )$$

where Q = Γ 3 H is the dissipation rate which describes the effectiveness at which the inflaton energy converts into radiation. The Q ≫ 1 corresponds to strong dissipation regime, and the Q ≪ 1 corresponds to weak dissipation regime [29]. The presence of extra friction term in the equation of motion of the inflaton suggests that radiation will not be redshifted during inflation, because, inflaton field continuously converts into radiation through dissipation

For inflation to occur and last long enough, the potential energy of the inflaton field must be dominated over both the kinetic energy of the inflaton field and the energy density of the radiation field, however, the radiation energy density still satisfies the condition ρ 1 / 4 γ &gt; H . That means T &gt; H during warm inflation and this condition is considered as the condition for warm inflation. It is also assumed that the production of radiation is quasi-stable. These approximations are called slow-roll approximations which can be quantified as

$$\frac { ( 1 + 2 \lambda \sqrt { \kappa } \phi ) } { 2 } \phi ^ { 2 } + \rho _ { \gamma } < < ( 1 + 4 \lambda \sqrt { \kappa } \phi ) V \\ - \quad \ast \quad \cdot \quad -$$

$$( 1 + 2 \overset { \kappa } { 2 \lambda \sqrt { \kappa } \phi } ) \ddot { \phi } < < 3 H \dot { \phi } ( 1 + 2 \lambda \sqrt { \kappa } \phi + Q )$$

$$\dot { \rho } _ { \gamma } < < 4 H \rho _ { \gamma }$$

$$\stackrel { \cdot \, \cdot \, \cdot } { \sqrt { \kappa } \dot { \phi } ^ { 2 } } < < H \dot { \phi } ^ { 2 }$$

Under these considerations, ρ ϕ ≃ V ( ϕ ) , and the dynamical equations (13), (16) and (15) read as

$$H ^ { 2 } = \frac { \kappa } { 3 } ( 1 + 4 \lambda \sqrt { \kappa } \phi ) V$$

$$3 H \dot { \phi } ( 1 + 2 \lambda \sqrt { \kappa } \phi + Q ) + ( 1 + 4 \lambda \sqrt { \kappa } \phi ) V ^ { \prime } + 4 \lambda \sqrt { \kappa } V = 0$$

$$\rho _ { \gamma } = \frac { \Gamma \dot { \phi } ^ { 2 } } { 4 H } = C \tau ^ { 4 }$$

where C = g ∗ π 2 30 is the Stefan-Boltzmann constant and g ∗ is the number of degrees of freedom for the radiation at temperature τ . (For calculation we will take C = 70 for g ∗ = 200 )

The slow-roll approximation can be parameterized by several slow-roll parameters namely ϵ , η , and β . Using equations (21), (22) and (23), we can derive the modified slow-roll parameters in the framework of f ( ϕ T ) gravity in terms of the inflaton potential V ( ϕ ) as

$$\epsilon \ & = \ \ - \frac { \dot { H } } { H ^ { 2 } } \\ & = \ \frac { 1 } { 2 \kappa ( 1 + 2 \lambda \sqrt { \kappa } \phi + Q ) } \left ( \frac { V ^ { \prime } } { V } + \frac { 4 \lambda \sqrt { \kappa } } { 1 + 4 \lambda \sqrt { \kappa } \phi } \right ) ^ { 2 }$$

$$\eta \ & = \ & - \frac { \ddot { \phi } } { H \dot { \phi } } \\ & = \ & \frac { 1 } { \kappa ( 1 + 2 \lambda \sqrt { \kappa } \phi + Q ) } \left ( \frac { V ^ { \prime \prime } } { V } + \frac { \lambda \sqrt { \kappa } ( 6 + 8 \lambda \sqrt { \kappa } \phi + 8 Q ) } { ( 1 + 4 \lambda \sqrt { \kappa } \phi ) ( 1 + 2 \lambda \sqrt { \kappa } \phi + Q ) } \frac { V ^ { \prime } } { V } \\ & - \ & \frac { 8 \lambda ^ { 2 } \kappa } { ( 1 + 2 \lambda \sqrt { \kappa } \phi + Q ) ( 1 + 4 \lambda \sqrt { \kappa } \phi ) } \right )$$

and

$$\beta \ & = \ - \frac { \dot { \rho } _ { \gamma } } { H \rho _ { \gamma } } \\ & = \ \frac { 1 } { \kappa ( 1 + 2 \lambda \sqrt { \kappa } \phi + Q ) } \frac { \Gamma ^ { \prime } } { \Gamma } \left ( \frac { V ^ { \prime } } { V } + \frac { 4 \lambda \sqrt { \kappa } } { 1 + 4 \lambda \sqrt { \kappa } \phi } \right )$$

The conditions satisfied by the slow-roll parameters for warm inflation to occur are ϵ ≪ 1 , | η | ≪ 1 and | β | ≪ 1 [39]. Violation of the slow-roll conditions marks the end of inflation.

Another significant parameter during inflation is the number of e-folding. This parameter quantifies the extent of expansion of the universe throughout the inflationary epoch. The number of e-foldings when the inflation field ϕ rolls from its value ϕ i to ϕ f is estimated as:

$$N \ & = \ \int _ { t _ { 1 } } ^ { t _ { 2 } } H d t = \int _ { \phi _ { i } } ^ { \phi _ { f } } \frac { H } { \delta } d \phi \\ & = \ - \kappa \int _ { \phi _ { i } } ^ { \phi _ { f } } \frac { ( 1 + 2 \lambda \sqrt { \kappa } \phi + Q ) } { \frac { V ^ { \prime } } { V } + \frac { 4 \lambda \sqrt { \kappa } } { 1 + 4 \lambda \sqrt { \kappa } } } d \phi$$

## 3.1 Perturbation Spectra

In this section, we develop the theory of cosmological perturbations for warm inflation in the f ( ϕ T ) theory of gravity. In warm inflationary models T &gt; H , the fluctuations of the inflaton field will be produced by thermal effects due to radiation. Consequently, the source of density fluctuations in warm inflationary models is the thermal fluctuations in the radiation field rather than quantum fluctuations. These thermal fluctuations in the radiation field produce fluctuations in the inflaton field. These fluctuations lead to perturbations in the metric which in turn perturb the equation of motion of inflaton field. In the following, we will calculate the fluctuations of the inflaton field. To compute the inflaton fluctuation we start with the perturbations of the FRW metric in the spatially flat gauge which is given by

$$d s ^ { 2 } = - ( 1 + 2 A ) d t ^ { 2 } + 2 a ( t ) \partial _ { i } B d t d x ^ { i } + a ^ { 2 } ( t ) \delta _ { i j } d x ^ { i } d x ^ { j }$$

Where A x, t ( ) and B x,t ( ) are spece-time dependent scalar perturbations of the metric. Under this perturbation, we expand the inflaton field ϕ x, t ( ) = ϕ t ( ) + δϕ x, t ( ) , where δϕ x, t ( ) is the linear response due to the thermal stochastic noise. Using the slow-roll conditions, we find the perturbed equation of the inflaton field in the momentum space given by,

$$& ( 1 + 2 \lambda \sqrt { \kappa } \phi ) \ddot { \delta } \phi _ { k } + 3 H ( 1 + 2 \lambda \sqrt { \kappa } \phi + Q ) \dot { \delta } \phi _ { k } + \frac { k ^ { 2 } } { a ^ { 2 } } ( 1 + 2 \lambda \sqrt { \kappa } \phi ) \delta \phi _ { k } = ( 1 + 2 \lambda \sqrt { \kappa } \phi ) \dot { \phi } \dot { A } \\ + & \quad \frac { k ^ { 2 } } { a } ( 1 + 2 \lambda \sqrt { \kappa } \phi ) \dot { \phi } B - ( 2 ( 1 + 4 \lambda \sqrt { \kappa } \phi ) V ^ { \prime } + 8 \lambda \sqrt { \kappa } V + \Gamma \dot { \phi } ) A$$

And the perturbed Einstein equation becomes

$$3 H ^ { 2 } A + \frac { k ^ { 2 } } { a } H B = - 4 \pi G \delta \rho$$

$$H A & = 4 \pi G ( \rho + p ) \delta u$$

In the above equations, ρ and p denote the total energy density and pressure respectively, δρ denotes the perturbations at the linear order of the total energy density, and δu is the scalar part of the linear perturbation of the 4- velocity. Substituting the expressions of A and B from equations (30) and (31) into equation (29), we get the perturbed equation of motion of inflation field as

$$( 1 + 2 \lambda \sqrt { \kappa } \phi ) \ddot { \delta } _ { k } ( t ) + 3 H ( 1 + 2 \lambda \sqrt { \kappa } \phi + Q ) \dot { \delta } \phi _ { k } ( t ) + \frac { k ^ { 2 } } { a ^ { 2 } } ( 1 + 2 \lambda \sqrt { \kappa } \phi ) \delta \phi _ { k } ( t ) = \xi _ { k } ( t ) \quad \quad ( 3 2 )$$

where the thermal stochastic noise source ξ k ( ) t is introduced to describe thermal fluctuations. In the slow-roll regime, the term ¨ δϕ k in equation (32) can be neglected. Thus equation (32) becomes

$$3 H ( 1 + 2 \lambda \sqrt { \kappa } \phi + Q ) \dot { \delta } \phi _ { k } ( t ) + \frac { k ^ { 2 } } { a ^ { 2 } } ( 1 + 2 \lambda \sqrt { \kappa } \phi ) \delta \phi _ { k } ( t ) = \xi _ { k } ( t )$$

The solution of equation (33) is

$$\delta \phi _ { k } ( t ) = \frac { 1 } { 3 H ( 1 + 2 \lambda \sqrt { \kappa } \phi + Q ) } \exp \left [ - \frac { t } { \mathcal { T } } \right ] \int _ { t _ { 0 } } ^ { t } \exp \left [ \frac { t ^ { \prime } } { \mathcal { T } } \right ] \xi _ { k } ( t ^ { \prime } ) d t ^ { \prime } + \delta \phi _ { k } ( t _ { 0 } ) \exp \left [ - \frac { t - t _ { 0 } } { \mathcal { T } } \right ] \quad \quad ( 3 4 )$$

√ √

where, T ( ϕ ) = 3 H (1+2 λ κϕ + Q ) (1+2 λ √ κϕ ) k 2 a 2 = 3 H (1+2 λ κϕ + Q ) (1+2 λ √ κϕ k ) 2 p and k p is the physical wave number. The first term in the right-hand side of the solution (34) is the noise contribution which acts to thermalize δϕ . whereas the second term is the memory term for the initial value of δϕ which is exponentially damping i.e. it becomes negligible with time. If k p of one mode of δϕ k is smaller than the freeze-out physical wave number k F , then the mode will not thermalize during the Hubble time, and for k p ≳ k F , the mode will be thermalized and the memory term will vanish within Hubble time. So, the freeze-out wave number k F is

$$k _ { F } = \sqrt { \frac { 3 H ^ { 2 } ( 1 + 2 \lambda \sqrt { \kappa } \phi + Q ) } { 1 + 2 \lambda \sqrt { \kappa } \phi } }$$

The power spectrum for the scalar fluctuations in warm inflation has the form

$$P _ { R } = \left ( \frac { H } { \dot { \phi } } \right ) ^ { 2 } \delta \phi ^ { 2 }$$

where the fluctuations of scalar field can be obtained through the relation

$$\delta \phi ^ { 2 } \ = \ \frac { k _ { F } \tau } { 2 \pi ^ { 2 } }$$

Combining equations (35),(36) and (37), the expression for scalar power spectrum P R for warm inflation in f ( ϕ T ) gravity leads to

$$P _ { R } \ = \ \frac { H ^ { 3 } \tau } { 2 \pi ^ { 2 } \dot { \phi } ^ { 2 } } \sqrt { \frac { 3 ( 1 + 2 \lambda \sqrt { \kappa } \phi + Q ) } { 1 + 2 \lambda \sqrt { \kappa } \phi } }$$

In the limit λ → 0 , the above expression goes back to the form of P R for warm inflation in Einstein gravity. The power spectrum for the tensor perturbation is

$$P _ { T } \ & = \ \frac { 1 6 H ^ { 2 } } { \pi M _ { P } ^ { 2 } } \\ & = \ \frac { 1 6 ( 1 + 4 \lambda \sqrt { \kappa } \phi ) V } { 3 \pi M _ { P } ^ { 4 } } \\ & \quad. \quad.$$

The spectral index n s and tensor to scalar ratio r are defined as [40, 41]

$$n _ { s } - 1 = \frac { d \ln P _ { R } } { d \ln k }$$

$$r = \frac { P _ { T } } { P _ { R } }$$

where P R is the scalar power spectrum and P T is the tensor power spectrum [39].

## 4 Warm Natural Inflation in f ( ϕ T ) Gravity with Temperature-Dependent Dissipation Coefficient

A successful inflationary model must generate sufficient expansion to solve the horizon problem. To satisfy this constraint, the slope of the inflaton potential must be very flat. This imposes a restriction on the choice of potentials. Moreover, it is expected that such a potential should be motivated by a fundamental theory. The natural inflation model is a well-motivated model where axion plays the role of inflaton field and the shift symmetry present in the axionic theory assures the flat potential required for inflation. We, therefore, attempt to study natural warm inflation in the framework of f ( ϕ, T ) gravity. The form of the natural potential we are considering is given by [42, 43, 44]

$$V ( \phi ) = \mu ^ { 4 } \left ( 1 + \cos \left ( \frac { \phi } { \mathfrak { f } } \right ) \right )$$

where f is the decay constant and µ is the mass scale of the axion. We assume that the inflaton couples with light gauge fields and produces thermal friction which is given by the dissipation coefficient of the form [45, 46, 47]

$$\Gamma ( \tau ) = C _ { \Gamma } \frac { \tau ^ { 3 } } { f ^ { 2 } }$$

where C Γ is a dimensionless factor proportional to the coupling constant between the inflaton field ϕ and gauge field.

With this form of dissipation coefficient, we carry out our further analysis in two dissipative regimes viz weak Q &lt;&lt; 1 and strong Q &gt;&gt; 1 .

## 4.1 Weak Dissipative Regime

In weak dissipative regime, Q &lt;&lt; 1 and the slow-roll parameters ϵ , η and β take the forms,

$$\epsilon = \frac { 1 } { 2 \kappa ( 1 + 2 \lambda \sqrt { \kappa } \phi ) } \left ( \frac { 4 \lambda \sqrt { \kappa } } { 4 \lambda \sqrt { \kappa } \phi + 1 } - \frac { \sin \left ( \frac { \phi } { f } \right ) } { f \left ( \cos \left ( \frac { \phi } { f } \right ) + 1 \right ) } \right ) ^ { 2 }$$

$$\eta \ = \ \frac { 1 } { \kappa ( 1 + 2 \lambda \sqrt { \kappa } \phi ) } \left ( - \frac { 8 \lambda ^ { 2 } \kappa } { ( 2 \lambda \sqrt { \kappa } \phi + 1 ) ( 4 \lambda \sqrt { \kappa } \phi + 1 ) } - \frac { \cos \left ( \frac { \phi } { f } \right ) } { f ^ { 2 } \left ( \cos \left ( \frac { \phi } { f } \right ) + 1 \right ) }$$

$$\begin{matrix} & & \quad \\ & & \quad \lambda \sqrt { \kappa } ( 8 \lambda \sqrt { \kappa } \phi + 6 ) \sin \left ( \frac { \phi } { f } \right ) & & \quad \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \colon \\ & - & \frac { \lambda \sqrt { \kappa } ( 8 \lambda \sqrt { \kappa } \phi + 6 ) \sin \left ( \frac { \phi } { f } \right ) } { f ( 2 \lambda \sqrt { \kappa } \phi + 1 ) ( 4 \lambda \sqrt { \kappa } \phi + 1 ) \left ( \cos \left ( \frac { \phi } { f } \right ) + 1 \right ) \end{matrix}$$

$$\beta = \frac { 1 } { \kappa ( 1 + 2 \lambda \sqrt { \kappa } \phi ) } \frac { \Gamma ^ { \prime } } { \Gamma } \left ( \frac { 4 \lambda \sqrt { \kappa } } { 1 + 4 \lambda \sqrt { \kappa } \phi } - \frac { \sin \left ( \frac { \phi } { f } \right ) } { f \left ( \cos \left ( \frac { \phi } { f } \right ) + 1 \right ) } \right )$$

and the equation of motion of the inflaton field reads

$$3 H \dot { \phi } ( 1 + 2 \lambda \sqrt { \kappa } \phi ) + ( 1 + 4 \lambda \sqrt { \kappa } \phi ) V ^ { \prime } + 4 \lambda \sqrt { \kappa } V = 0$$

Using equations (21) and (42), we obtain

$$\dot { \phi } = - \frac { \sqrt { \mu ^ { 4 } ( 4 \lambda \sqrt { \kappa } \phi + 1 ) \left ( \cos \left ( \frac { \phi } { f } \right ) + 1 \right ) } \left ( \frac { 4 \lambda \sqrt { \kappa } } { 4 \lambda \sqrt { \kappa } \phi + 1 } - \frac { \sin \left ( \frac { \phi } { f } \right ) } { f ( \cos \left ( \frac { \phi } { f } \right ) + 1 \right ) } \right ) } { \sqrt { 3 } \kappa ( 2 \lambda \sqrt { \kappa } \phi + 1 ) }$$

The temperature of the thermal bath can be obtained by using equation (23) as

$$\tau = \left ( \frac { 3 Q \dot { \phi } ^ { 2 } } { 4 C } \right ) ^ { \frac { 4 } { 4 } }$$

Also, the relation (43), reads

$$\tau = \left ( \frac { 3 Q H \mathfrak { f } ^ { 2 } } { C _ { \Gamma } } \right ) ^ { \frac { 1 } { 3 } }$$

Equating equations (49) and (50), we get

$$Q = \frac { \left ( \frac { 1 } { 4 C } \right ) ^ { 3 } \left ( C _ { T } ^ { 4 } \mu ^ { 4 } \right ) \left ( ( 4 \lambda \sqrt { \kappa } \phi + 1 ) \left ( \cos \left ( \frac { \phi } { f } \right ) + 1 \right ) \right ) \left ( \frac { 4 \lambda \sqrt { \kappa } f } { 4 \lambda \sqrt { \kappa } \phi + 1 } - \frac { \sin \left ( \frac { \phi } { f } \right ) } { \cos \left ( \frac { \phi } { f } \right ) + 1 } \right ) ^ { 6 } } { ( 9 \kappa ^ { 5 } f ^ { 1 4 } ) \left ( 2 \lambda \sqrt { \kappa } \phi + 1 \right ) ^ { 6 } }$$

The violation of the slow-roll conditions governs the end of inflation. It is numerically checked that ϵ violates the slow-roll condition first. So, by using the condition ϵ ϕ ( f ) = 1 , we determine the value of the scalar field at the end of inflation. The number of e-foldings (equation (27)) in this regime becomes,

$$N \ = \ - \kappa \int _ { \phi _ { i } } ^ { \phi _ { f } } \frac { ( 1 + 2 \lambda \sqrt { \kappa } \phi ) } { \frac { V ^ { \prime } } { V } + \frac { 4 \lambda \sqrt { \kappa } } { 1 + 4 \lambda \sqrt { \kappa } \phi } } d \phi$$

Using equation (52), the initial field value ( ϕ i ) is calculated for N = 60 . The temperature of the thermal bath can be obtained by inserting the equations (48) and (51) into equations (49),

$$\tau = 0. 1 4 4 \sqrt { \frac { C _ { \Gamma } ^ { 4 } \mu ^ { 8 } ( 4 \lambda \sqrt { \kappa } \phi + 1 ) ^ { 2 } \left ( \cos \left ( \frac { \phi } { \gamma } \right ) + 1 \right ) ^ { 2 } \left ( \frac { 4 \lambda \sqrt { \kappa } \gamma \gamma } { 4 \lambda \sqrt { \kappa } \phi + 1 } - \frac { \sin ( \frac { \phi } { \gamma } ) } { \cos ( \frac { \phi } { \gamma } ) + 1 } \right ) ^ { 6 } \left ( \frac { 4 \lambda \sqrt { \kappa } } { 4 \lambda \sqrt { \kappa } \phi + 1 } - \frac { \sin ( \frac { \phi } { \gamma } ) } { \gamma ( \cos ( \frac { \phi } { \gamma } ) + 1 ) } \right ) ^ { 2 } } } { \kappa ^ { 5 } C ^ { 4 } \gamma ^ { 1 4 } ( 2 \lambda \sqrt { \kappa } \phi + 1 ) ^ { 8 } } \ \ ( 5 3 )$$

From the above relation, we can obtain the temperature of the thermal bath for any value of ϕ . Also, by substituting the relations (38), (39), (48) and (51) into equations (40) and (41), we obtain the spectral index and tensor-to-scalar ratio at initial field value ϕ i .

To evaluate the viability of any inflationary model, it is necessary to check whether the values of cosmological observables are in good agreement with the observational data. So, in order to constrain different parameters used in our model, we consider the latest limits set by Planck 2018 and BICEP/Keck data on n s and r .

The predictions of the spectral tilt and tensor-to-scalar ratio for three representative values of dissipation parameters are shown in Fig. 1. The plots are obtained by varying λ and fixing the values of the mass scale µ = 0 001 . M P and axion decay constant f = 0 4 . M P for 60 e-foldings. It is seen that variation of tensor-toscalar ratio is very slow and with the increasing dissipation parameter ( C Γ ) , its values go on decreasing. The ranges of model parameter λ are found to be [93 6 . , 97 3] . for C Γ = 1 1 . × 10 5 , [93 57 . , 97 27] . for C Γ = 5 × 10 5 and [93 47 . , 97 17] . for C Γ = 7 38 . × 10 5 . The plot shows that the warm natural inflationary scenario in the f ( ϕ T ) gravity framework is consistent with observational constraints on r and n s in weak dissipative regime.

## 4.2 Strong Dissipative Regime

In strong dissipative regime, Q &gt;&gt; 1 . So, the slow-roll parameters ϵ , η and β in this regime become,

$$\epsilon = \frac { 1 } { 2 \kappa ( 2 \lambda \sqrt { \kappa } \phi + Q ) } \left ( \frac { 4 \lambda \sqrt { \kappa } } { 4 \lambda \sqrt { \kappa } \phi + 1 } - \frac { \sin \left ( \frac { \phi } { f } \right ) } { f \left ( \cos \left ( \frac { \phi } { f } \right ) + 1 \right ) } \right ) ^ { 2 }$$

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00595v1-figure-e57d62f0-94af-4cdb-8601-e5fb10038c45.jpg)

r

(a)

C

Γ

= 1 1

.

×

10

5

Figure 1: The spectral index n s , and the tensor-to-scalar ratio r predicted by warm natural inflationary model in the f ( ϕ T ) gravity theory for N = 60 in weak dissipative regime (green solid line). The marginalized joint 95% CL regions for the spectral index n s , and the tensor-to-scalar ratio r from Planck 2018 data alone, and in combinations with BICEP/Keck Array are shown in blue and orange respectively.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00595v1-figure-cba97d07-04df-40d5-8839-883280afb41a.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00595v1-figure-339c8b55-c097-4cca-a752-e3b2247fe91d.jpg)

$$\eta \ = \ \frac { 1 } { \kappa ( 2 \lambda \sqrt { \kappa } \phi + Q ) } \left ( - \frac { 8 \lambda ^ { 2 } \kappa } { ( 2 \lambda \sqrt { \kappa } \phi + Q ) ( 4 \lambda \sqrt { \kappa } \phi + 1 ) } - \frac { \cos \left ( \frac { \phi } { f } \right ) } { f ^ { 2 } \left ( \cos \left ( \frac { \phi } { f } \right ) + 1 \right ) }$$

$$- \ \frac { \lambda \sqrt { \kappa } ( 8 \lambda \sqrt { \kappa } \phi + 6 + 8 Q ) \sin \left ( \frac { \phi } { f } \right ) } { \epsilon \in \Gamma _ { 1 } + \alpha \sqrt { \kappa } \in \Gamma _ { 2 } + 1 \times ( \dots ) \left ( \phi \times \dots \right ) } } \right )$$

$$- \ \frac { \lambda \sqrt { \kappa } ( 8 \lambda \sqrt { \kappa \phi } + 6 + 8 Q ) \sin \left ( \frac { \phi } { f } \right ) } { f ( 2 \lambda \sqrt { \kappa } \phi + Q ) ( 4 \lambda \sqrt { \kappa } \phi + 1 ) \left ( \cos \left ( \frac { \phi } { f } \right ) + 1 \right ) } \right )$$

$$\beta = \frac { 1 } { \kappa ( 2 \lambda \sqrt { \kappa } \phi + Q ) } \frac { \Gamma ^ { \prime } } { \Gamma } \left ( \frac { 4 \lambda \sqrt { \kappa } } { 1 + 4 \lambda \sqrt { \kappa } \phi } - \frac { \sin \left ( \frac { \phi } { f } \right ) } { f \left ( \cos \left ( \frac { \phi } { f } \right ) + 1 \right ) } \right )$$

and the equation of motion of the inflaton field becomes

$$3 H \dot { \phi } ( 2 \lambda \sqrt { \kappa } \phi + Q ) + ( 1 + 4 \lambda \sqrt { \kappa } \phi ) V ^ { \prime } + 4 \lambda \sqrt { \kappa } V = 0$$

$$\text{From } \Gamma ( \tau ) = C _ { \Gamma } \frac { \tau ^ { 3 } } { \tilde { \jmath } ^ { 2 } }, \, \text{we have } \tau = \left ( \frac { 3 Q \hat { \jmath } ^ { 2 } H } { C _ { \Gamma } } \right ) ^ { \frac { 1 } { 3 } }. \, \text{Equating this with } \tau = \left ( \frac { 3 Q \hat { \jmath } ^ { 2 } } { 4 C } \right ) ^ { \frac { 1 } { 4 } }, \, \text{we get}$$

$$\dot { \phi } ^ { 6 } = 3 Q ( 4 C ) ^ { 3 } \left ( \frac { f ^ { 2 } H } { C _ { \Gamma } } \right ) ^ { 4 }$$

Expressions of ˙ ϕ and Q can be obtained by solving equations (57) and (58). These expressions are substituted in equations (54), (55) and (56) to obtain the slow-roll parameters ϵ , η and β in terms of model parameters. In this case, we have numerically checked that inflation ends when ϵ reaches unity and found the corresponding value of inflaton field ( ϕ f ) . The number of e-foldings (equation (27)) in this regime becomes,

$$N \ = \ - \kappa \int _ { \phi _ { i } } ^ { \phi _ { f } } \frac { ( 2 \lambda \sqrt { \kappa } \phi + Q ) } { \frac { V ^ { \prime } } { V } + \frac { 4 \lambda \sqrt { \kappa } } { 1 + 4 \lambda \sqrt { \kappa } \phi } } d \phi$$

Using the above definition of N , we obtain the initial field value ( ϕ i ) for 60 e-foldings. From the expressions of the power spectrum P R (equation (38)) and P T (equation (39)), the spectral index n s and tensor-to-scalar ratio r are computed at ϕ i for different values of model parameters.

The depictions of the spectral tilt and tensor-to-scalar ratio are shown in Fig. 2. Here also we consider three representative values of C Γ . In each case, the mass scale and the axion decay constant are set at µ = 0 001 . M P , f = 0 4 . M P . The range of model parameter λ is found to be [88 3 . , 95 1] . for C Γ = 10 6 , [95 75 . , 104 3] . for C Γ = 1 1 . × 10 6 and [103 2 . , 114] for C Γ = 1 2 . × 10 6 . The values of n s and r are in good agreement with Planck 2018 [38] and BICEP [37] data at 95% CL.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00595v1-figure-3bc50775-a884-4c85-b8df-46af4fc06803.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00595v1-figure-dcf93df6-c29f-421a-b6b2-e78b7c82e400.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00595v1-figure-25595cfc-d1dd-4aaa-be3e-bb367a3dcee8.jpg)

r

Figure 2: The spectral index n s , and the tensor-to-scalar ratio r predicted by warm natural inflationary model in the f ( ϕ T ) gravity theory for N = 60 in strong dissipative regime (green solid line). The marginalized joint 95% CL regions for the spectral index n s , and the tensor-to-scalar ratio r from Planck 2018 data alone, and in combinations with BICEP/Keck Array are shown in blue and orange respectively.

## 5 Conclusion

Over the past few decades, extensive investigations have been conducted to understand and model the dynamics of the universe. The standard model of cosmology successfully describes the formation and evolution of the universe. However, this model has some problems like horizon and flatness problems which remain as open issues. In order to overcome these shortcomings, the idea of inflationary cosmology appears necessary in the early era of the universe. Even though general relativity provides remarkably accurate predictions for describing cosmological phenomena, it can not explain the dark sector of the universe. For this reason, investigating alternative theories of gravity arises as a promising path.

In this work, we explore natural warm inflation in the framework of f ( ϕ T ) gravity theory. In warm inflation, the interactions between scalar and other fields during inflation are taken into account which provides the dissipation term. Here we consider a cubic temperature-dependent dissipation coefficient. We have investigated this model in two separate dissipative regimes viz weak and strong and in each case, we have evaluated the cosmological observables for certain ranges of model parameters.

In the weak dissipative regime, the cosmological observables are functions of model parameter λ , axion decay constant f , mass scale µ , and dissipation coefficient. We have presented the predictions of the spectral tilt and tensor-to-scalar ratio for dissipation parameter C Γ = 1 1 . × 10 5 , C Γ = 5 × 10 5 and C Γ = 7 38 . × 10 5 by varying λ and keeping the other parameters fixed. The obtained value for the spectral index and the tensorto-scalar ratio in the selected range of the model parameter λ falls in the 1 σ confidence level of the Planck 2018 data and the joint Planck and BICEP results. Hence, we claim that this model in weak dissipative regime is capable of making predictions consistent with the observational data.

In the strong dissipative regime also, the cosmological observables are the function of model parameter λ , axion decay constant f , mass scale µ , and dissipation coefficient as expected. Here, we have plotted the predictions of the spectral tilt and tensor-to-scalar ratio for dissipation parameter C Γ = 10 6 , C Γ = 1 1 . × 10 6 and C Γ = 1 2 . × 10 6 by varying λ and keeping the other parameters fixed. The result shows that for some selected ranges of the model parameter λ , values of the spectral index, and the tensor-to-scalar ratio match in the 1 σ confidence level of the Planck 2018 data and their combination with BICEP/Keck data. Hence, this model is consistent with the observational data in strong disspative regime also.

It is found that the model parameter λ has a minimal impact on the tensor-to-scalar ratio, as its values exhibit slow variation. So, it can be concluded that the tensor-to-scalar ratio is less sensitive to the change in the model parameter. The model produces smaller values of tensor-to-scalar ratio in strong dissipative regime. Because, the scalar power spectrum is completely dominated by the dissipation, which strongly suppresses the tensor-to-scalar ratio. Additionally, the value of the axion decay constant can be lowered below the Planck scale which is 0 4 . M P in this model. We can realize the natural warm inflation scenario in f ( ϕ T ) gravity for much smaller values of f (sub-Planckian scale) compared to the corresponding cold inflation case, which makes it theoretically more sound.

## Acknowledgement

We are thankful to Prof. Sudhakar Panda for fruitful insights.

## References

| [1]   | Alam, U., Sahni, V., and Starobinsky, A. A. (2007). Exploring the properties of dark energy using type-Ia supernovae and other datasets. Journal of Cosmology and Astroparticle Physics, 2007(02), 011.            |
|-------|--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| [2]   | Ishak, M. (2019). Testing general relativity in cosmology. Living Reviews in Relativity, 22, 1-204.                                                                                                                |
| [3]   | Ferreira, P. G. (2019). Cosmological tests of gravity. Annual Review of Astronomy and Astrophysics, 57(1), 335-374.                                                                                                |
| [4]   | Sotiriou, T. P., and Faraoni, V. (2010). f ( R ) theories of gravity. Reviews of Modern Physics, 82(1), 451-497.                                                                                                   |
| [5]   | Nojiri, S. I., and Odintsov, S. D. (2007). Introduction to modified gravity and gravitational alternative for dark energy. International Journal of Geometric Methods in Modern Physics, 4(01), 115-145.           |
| [6]   | De Felice, A., and Tsujikawa, S. (2010). f ( R ) theories. Living Reviews in Relativity, 13(1), 1-161.                                                                                                             |
| [7]   | Zaregonbadi, R., Farhoudi, M., and Riazi, N. (2016). Dark matter from f ( R,T ) gravity. Physical Review D, 94(8), 084052.                                                                                         |
| [8]   | Haghani, Z., Shiravand, M., and Shahidi, S. (2018). Energy conditions in mimetic- f ( R ) gravity. Inter- national Journal of Modern Physics D, 27(05), 1850049.                                                   |
| [9]   | Shankaranarayanan, S., and Johnson, J. P. (2022). Modified theories of gravity: Why, how and what?. General Relativity and Gravitation, 54(5), 44.                                                                 |
| [10]  | Nojiri, S. I., and Odintsov, S. D. (2003). Modified gravity with negative and positive powers of curvature: Unification of inflation and cosmic acceleration. physical Review D, 68(12), 123512.                   |
| [11]  | Nojiri, S. I., and Odintsov, S. D. (2007). Introduction to modified gravity and gravitational alternative for dark energy. International Journal of Geometric Methods in Modern Physics, 4(01), 115-145.           |
| [12]  | Nojiri, S., Odintsov, S. D., and Oikonomou, V. (2017). Modified gravity theories on a nutshell: Inflation, bounce and late-time evolution. Physics Reports, 692, 1-104.                                            |
| [13]  | Nojiri, S., Odintsov, S. D., and Oikonomou, V. K. (2019). k-essence f ( R ) gravity inflation. Nuclear Physics B, 941, 11-27.                                                                                      |
| [14]  | Nojiri, S. I., and Odintsov, S. D. (2007, May). Modified gravity and its reconstruction from the universe expansion history. In Journal of Physics: Conference Series (Vol. 66, No. 1, p. 012005). IOP Publishing. |
| [15]  | Harko, T., Lobo, F. S., Nojiri, S. I., and Odintsov, S. D. (2011). f ( R,T ) gravity. Physical Review D-Particles, Fields, Gravitation, and Cosmology, 84(2), 024020.                                              |
| [16]  | Bhattacharjee, S., Santos, J. R. L., Moraes, P. H. R. S., and Sahoo, P. K. (2020). Inflation in f ( R,T gravity. The European Physical Journal Plus, 135(7), 576.                                                  |
| [17]  | Yeasmin, S., Deb, B., and Deshamukhya, A. (2023). Warm inflation in f ( R,T ) gravity. Chinese Journal of Physics, 85, 359-374.                                                                                    |
| [18]  | Deb, B., and Deshamukhya, A. (2022). Inflation in f ( R,T ) gravity with double-well potential. Interna- tional Journal of Modern Physics A, 37(18), 2250127.                                                      |
| [19]  | Ashmita, Sarkar, P., and Das, P. K. (2022). Inflationary cosmology in the modified f ( R,T ) gravity. International Journal of Modern Physics D, 31(16), 2250120.                                                  |

| [20]   | Gamonal, M. (2021). Slow-roll inflation in f ( R,T ) gravity and a modified Starobinsky-like inflationary model. Physics of the Dark Universe, 31, 100768.                                                                                 |
|--------|--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| [21]   | Maurya, S. K., Banerjee, A., and Tello-Ortiz, F. (2020). Buchdahl model in f ( R,T ) gravity: A com- parative study with standard Einstein's gravity. Physics of the Dark Universe, 27, 100438.                                            |
| [22]   | Jamil, M., Momeni, D., and Myrzakulov, R. (2015). Warm intermediate inflation in F ( T ) gravity. International Journal of Theoretical Physics, 54, 1098-1112.                                                                             |
| [23]   | De Laurentis, M., Paolella, M., and Capozziello, S. (2015). Cosmological inflation in F ( R,G ) gravity. Physical Review D, 91(8), 083531.                                                                                                 |
| [24]   | Sharif, M., and Ikram, A. (2016). Warm inflation in f ( G ) theory of gravity. Journal of Experimental and Theoretical Physics, 123, 40-50.                                                                                                |
| [25]   | Tirandari, M., and Saaidi, K. (2017). Anisotropic inflation in Brans-Dicke gravity. Nuclear Physics B, 925, 403-414.                                                                                                                       |
| [26]   | Zhang, X., Chen, C. Y., and Reyimuaji, Y. (2022). Modified gravity models for inflation: In conformity with observations. Physical Review D, 105(4), 043514.                                                                               |
| [27]   | Banks, T., Dine, M., Fox, P. J., and Gorbatov, E. (2003). On the possibility of large axion decay constants. Journal of cosmology and astroparticle physics, 2003(06), 001.                                                                |
| [28]   | Berera, A., and Fang, L. Z. (1995). Thermally induced density perturbations in the inflation era. Physical Review Letters, 74(11), 1912.                                                                                                   |
| [29]   | Berera, A. (1995). Warm inflation. Physical Review Letters, 75(18), 3218.                                                                                                                                                                  |
| [30]   | Berera, A. (1997). Interpolating the stage of exponential expansion in the early universe: Possible alternative with no reheating. Physical Review D, 55(6), 3346.                                                                         |
| [31]   | Berera, A. (2000). Warm inflation in the adiabatic regime-a model, an existence proof for inflationary dynamics in quantum field theory. Nuclear Physics B, 585(3), 666-714.                                                               |
| [32]   | Berera, A., Gleiser, M., and Ramos, R. O. (1998). Strong dissipative behavior in quantum field theory. Physical Review D, 58(12), 123508.                                                                                                  |
| [33]   | Deshamukhya, A., and Panda, S. (2009). Warm tachyonic inflation in a warped background. Interna- tional Journal of Modern Physics D, 18(14), 2093-2106.                                                                                    |
| [34]   | Yeasmin, S., and Deshamukhya, A. (2022). Effect of Dissipation on Warm Chromo-Natural Inflation. arXiv preprint arXiv:2203.12213.                                                                                                          |
| [35]   | De Oliveira, H. P., and Ramos, R. O. (1998). Dynamical system analysis for inflation with dissipation. Physical Review D, 57(2), 741.                                                                                                      |
| [36]   | Hall, L. M., Moss, I. G., and Berera, A. (2004). Scalar perturbation spectra from warm inflation. Physical Review D, 69(8), 083525.                                                                                                        |
| [37]   | Ade, P. A., et al. (BICEP/Keck Collaboration). (2021). Improved constraints on primordial gravitational waves using Planck, WMAP, and BICEP/Keck observations through the 2018 observing season. Physical review letters, 127(15), 151301. |
| [38]   | Akrami, Y., et al. (2020). Planck 2018 results-X. Constraints on inflation. Astronomy and Astrophysics, 641, A10.                                                                                                                          |
| [39]   | Berera, A., Moss, I. G., and Ramos, R. O. (2009). Warm inflation and its microphysical basis. Reports on Progress in Physics, 72(2), 026901.                                                                                               |
| [40]   | Taylor, A. N., and Berera, A. (2000). Perturbation spectra in the warm inflationary scenario. Physical Review D, 62(8), 083517.                                                                                                            |

| [41]   | Peng, Z. P., Yu, J. N., Zhang, X. M., and Zhu, J. Y. (2018). Perturbation spectra in the warm k-inflation. Physical Review D, 97(6), 063523.                                            |
|--------|-----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| [42]   | Freese, K., Frieman, J. A., and Olinto, A. V. (1990). Natural inflation with pseudo Nambu-Goldstone bosons. Physical Review Letters, 65(26), 3233.                                      |
| [43]   | Mishra, H., Mohanty, S., and Nautiyal, A. (2012). Warm natural inflation. Physics Letters B, 710(2), 245-250.                                                                           |
| [44]   | Visinelli, L. (2011). Natural warm inflation. Journal of Cosmology and Astroparticle Physics, 2011(09), 013.                                                                            |
| [45]   | Berghaus, K. V., Graham, P. W., and Kaplan, D. E. (2020). Minimal warm inflation. Journal of Cos- mology and Astroparticle Physics, 2020(03), 034.                                      |
| [46]   | Montefalcone, G., Aragam, V., Visinelli, L., and Freese, K. (2023). Observational constraints on warm natural inflation. Journal of Cosmology and Astroparticle Physics, 2023(03), 002. |
| [47]   | Reyimuaji, Y., and Zhang, X. (2021). Warm-assisted natural inflation. Journal of Cosmology and As- troparticle Physics, 2021(04), 077.                                                  |
Sabina Yeasmin Atri Deshamukhya
Embedding Warm Natural Inflation in $f(φ)T$ gravity
We study warm inflation in the framework of $f(\phi)T$ gravity, where $\phi$
is the inflaton and $T$ is the trace of the energy-momentum tensor. The
inflaton field is assumed to roll on the natural potential and the result is
analyzed in light of Planck 2018 and BICEP/Keck 2021 data. We start our work by
obtaining the field equations under slow-roll approximations. We then evaluate
the scalar and tensor power spectra and their corresponding spectral index and
tensor-to-scalar ratio with a temperature-dependent form of the dissipation
coefficient during the inflationary era. We find that the warm inflation model
in $f(\phi)T$ gravity is compatible with observational bands.