## Friedmann cosmology with particle creation in modified f ( R, T ) gravity

C. P. Singh 1 and Vijay Singh 2

1 2 , Department of Applied Mathematics, Delhi Technological University (Formerly Delhi College of Engineering) Bawana Road, Delhi-110 042, India. 1 cpsphd@rediffmail.com 2 gtrcosmo@gmail.com

## Abstract

The theoretical and observational consequences of thermodynamics of open systems, which allow particle creation are investigated in modified f ( R,T ) ( R is the Ricci scalar and T is the trace of energy-momentum tensor) theory of gravity within the framework of a flat Friedmann-Robertson-Walker line element. A simplest particular model f ( R,T ) = R +2 ( f T ) and 'gamma-law' equation of state p = ( γ -1) ρ are assumed to explore the role of particle creation in the early and late time evolution of the universe. A power-law expansion model is proposed for 0 ≤ γ ≤ 2 by considering the natural phenomenological particle creation rate ψ = 3 βnH , where β is a pure number of the order of unity. A Big Rip singularity is observed for γ &lt; 0, which describes the phantom cosmology. We observe that the accelerated expansion of the universe is driven by the particle creation without any exotic 'quintessence' component or a cosmological constant. It is also found by density parameter that the model becomes of negative curvature due to particle creation and the entropy increases with the evolution of the universe. Some kinematics tests such as lookback time, luminosity distance, proper distance, angular diameter versus redshift are discussed in detail to observe the role of particle creation in early and late time evolution of the universe.

Keywords: Modified gravity theory; Particle creation.

PACS number(s): 98.80-k, 98.80-cq, 04.20-q.

1 Corresponding author

## 1. Introduction

Many theoretical and observational studies of the universe, such as Type Ia supernovae 1 -3 , cosmic microwave background (CMB) anisotropy , large scale structures (LSS) 4 5 etc., have been shown that a pure Friedmann-Lamaitre-Robertson-Walker (FLRW) cosmology with matter and radiation could not explain all the large scale properties of our universe. The observations suggest that our universe is spatially flat and consists of about 70 percent dark energy with negative pressure, nearly 25 percent of dark matter, which populates the galaxy halos and 5 percent of baryonic matter. The late time cosmic acceleration is usually assumed to be driven by a cosmic fluid/field generically known as dark energy (DE) 6 -8 , which can be observed only by its gravitational effects 9 -11 . Since the time of this discovery in 1998, a large number of probable candidates of DE have been proposed to explain late time acceleration of the universe. The two most accepted DE models are that of a cosmological constant 12 (assumed possibly to be the quantum vacuum energy) and a slowly rolling scalar field (quintessence models) 13 -17 . All of these DE cosmological models are based on the Big Bang cosmology. In the context of early universe, the standard Big Bang model presents theoretical and observational difficulties, such as the singularity problem 18 , flatness and horizon problem, reheating during the inflationary epoch 19 , confliction between the age of the universe and the age of the oldest stars in globular clusters (age problem) 20 , the entropy problem 21 22 , etc.

The flatness and horizon problems together with the entropy problem have been resolved (up to a certain extent) by the inflationary model proposed by Guth 23 and Linde 24 . The age confliction 25 is not an isolated complication, it comes with another serious one trouble that is structure formation through gravitational amplification of small primeval density perturbation. These issues open the door of investigations of different alternative theories of gravity for studying not only unpredicted late time cosmic acceleration but early universe as well 26 -29 . Many pioneer concepts have been introduced by theorists to alleviate the problems of early and late universes. Among the ways to resolve the problems of early universe, Dirac's large number hypothesis inspired a class of new cosmology named particle creation 30 -32 . The steady state model introduced by Bondi and Gold 33 on the foundation of perfect cosmological principle (PCP) also asserts the continuous generation of matter in the universe. Hoyle 34 and Narlikar 35 independently proposed a creation field theory and studied the matter creation during the evolution of the universe. Tryon 36 and Fomin 37 in their individual work proposed a theoretical concept of the creation of the universe as a vacuum fluctuation. Brout et al. 38 -41 putted a strong foundation of simultaneous creation of matter and curvature from a quantum fluctuation of the Minkowskian space-time vacuum.

Later on, Gunzig et al. 42 and Prigogine et al. 43 established the theoretical scenario of matter creation in the framework of cosmology. The authors applied the thermodynamics of an open system to cosmology allowing particle creation and entropy production. They showed that the second law of thermodynamics might be modified to accommodate the flow of energy from gravitational field to the matter field, resulting in the creation of particles and consequently entropy. The work might suggest that at the expense of the gravitational field, particle creation takes place as an irreversible process constrained by the usual requirements of the non-equilibrium thermodynamics, however, the reverse process (matter destruction)

thermodynamically forbidden. Calv´o et al. a 44 extended this new concept of matter creation under adiabatic conditions. The further results were generalized by Lima and Germano 45 through a contravariant formulation allowing specific entropy variation as usually expected for non-equilibrium process in fluids. Lima and Alcaniz 46 and Alcaniz and Lima 47 investigated observational consequences of FRW models driven by adiabatic matter creation through some kinematics tests. After the discovery of accelerating universe the particle creation theory was reconsidered to explain it and got unexpected results. The negative pressure due to particle creation, might play the role of exotic matter component. Zimdahl et al. 48 tested the particle creation with SNe Ia data to get the result of accelerating universe. Qiang et al. 49 studied a model with adiabatic particle creation which is also consistent with SNe Ia data. Singh and Beesham 50 51 , , and Singh 52 53 , studied early universe in FRW cosmology with particle creation through some kinematical tests. This theoretical formulation of continuous creation of matter in the universe may reinterpret several predictions of the standard Big Bang cosmology.

On the other hand, many modifications of general relativity (GR) have been proposed to reconcile the problems plagued in cosmology 54 . Among such modifications the simplest one is replacing the Ricci scalar curvature R with the general function of R called f ( R ) theories of gravity, describe transition from decelerated to accelerated expansion of the universe 55 . Bertolami et al. 56 have generalized f ( R ) theories by introducing an explicit coupling between arbitrary function of the Ricci scalar R and the matter Lagrangian density. Harko 57 has extended this new concept to the arbitrary coupling between matter and geometry. Harko et al. 58 have considered another extension of GR, where the gravitational Lagrangian is given by an arbitrary function of the Ricci scalar R and of the trace T of the stress-energy tensor, so called f ( R,T ) theory of gravity. The justification of choosing T as an argument for the Lagrangian is from exotic imperfect fluids or quantum effects. The authors have argued that due to the coupling between matter and geometry the theory depends on a source term, representing the variation of the matter-stress-energy tensor with respect to the metric. Consequently, the cosmic acceleration in f ( R,T ) gravity results not only from geometrical effect but also from the matter contribution. These interesting features of f ( R,T ) gravity motivate many authors to study it for resolving several issues of current interest in cosmology and astrophysics.

Jamil et al. 59 have reconstructed cosmological models in the framework of f ( R,T ) gravity, showing dust fluid Λ CDM , phantom-non-phantom era and phantom cosmology, which are good agreement with BAO observational data for low-redshift z &lt; 2. Houndjo 60 has reconstructed F R,T ( ) = f 1 ( R ) + f 2 ( T ) using auxiliary scalar field with two known examples of scale factors corresponding to an accelerating universe. Houndjo and Piatella 61 have numerically reconstructed f ( R,T ) theory with dark matter and holographic DE, describing the same expansion history in GR. Houndjo et al. 62 have investigated f ( R,T ) gravity models which reproduce four types of future finite-time singularities. Sharif and Zubair 63 have discussed a non-equilibrium thermodynamics in f ( R,T ) gravity by taking two forms of the energy-momentum tensor of dark components which endorses second law of thermodynamics holds both in phantom and non-phantom phases. Jamil et al. 64 have proved, in general, the first law of black hole thermodynamics is violated in f ( R,T ) gravity. Alvarenga et al. 65 have paid special attention on f ( R,T ) = R +2 ( f T ) assuming special function f ( T ) showing

energy conditions can be satisfied for suitable input parameters. Sharif et al. 66 and Sharif and Zubair 67 have derived energy conditions in the context of f ( R,T ) gravity, which can reduce to the well-known conditions in f ( R ) and GR. Chakraborty 68 has shown that a part of an arbitrary function f ( R,T ) can be determined taking into account of the conservation of stress-energy tensor. Alvarenga et al. 69 have studied the evolution of scalar cosmological perturbations in the background of metric formalism in f ( R,T ) theory. Pasqua et al. 70 have considered modified holographic Ricci DE (MHRDE) model in f ( R,T ) modified gravity, which explained a quintessence-like behavior of the model. Shabani and Farhaudi 71 have studied the cosmological solutions of f ( R,T ) gravity through phase space analysis. Sharif and Zubair 72 have considered f ( R,T ) theory as an effective description for the holographic and new agegraphic DE and reconstructed the corresponding models, which explain phantom or quintessence regimes of the universe. Recently, Singh and Singh 73 have reconstructed the particular form f ( R,T ) = R +2 ( f T ) by assuming de-Sitter and power-law expansion of scale factors and showed the reconstructed forms of f ( R,T ) could explain DE models successfully. However, all these mentioned works have been carried out in FRW spacetime. Some authors have also explored f ( R,T ) theory in anisotropic background 74 -77 and higher dimension 78 . Thus, the newly proposed f ( R,T ) gravity motivates us to study the particle creation phenomena in the early and late time evolution of the universe.

In this paper, we investigate the theoretical and observational significance of particle creation in modified f ( R,T ) gravity theory in the framework of a flat Friedmann-RobertsonWalker model. Exact solutions of the field equations are obtained by assuming the suitable form of f ( R,T ) = R + 2 f ( T ), equation of state, and particle creation rate. We obtain the non-singular power-law solution for the scale factor. A Big Rip singularity is also observed which describes the phantom cosmology. We study some kinematical tests to explain the physical significance of particle creation during early and late time evolution of the universe.

The paper is organized as follows. The thermodynamics of particle creation is presented in the Sec. 2. A brief review of modified f ( R,T ) gravity theory and its field equations are given in Sec. 3. In Sec. 4, the model and field equations are presented with particle creation. The exact solutions of the field equations with suitable assumptions are obtained in Sec. 5. Sec. 6 and its subsections are devoted to study some kinematics tests of the model. Finally, the outcomes are summarized in Sec. 7.

## 2. Thermodynamics of particle creation

If we regard the whole universe as a closed thermodynamical system in which number of the particles in a given volume is constant, then the laws of thermodynamics have the form

$$d ( \rho _ { m } V ) = d Q - p _ { m } d V,$$

and

$$T d S = d ( \rho _ { m } V ) + p _ { m } d V,$$

where ρ m , V , p m , T and S are the energy density, volume, thermodynamical pressure, temperature, and entropy, respectively. Here, dQ is the heat received by the system during time dt . It is observed from (1) and (2) that the entropy production is given by

$$T d S = d Q.$$

Eq. (3) shows that the entropy remains stationary, i.e., dS = 0 for a closed adiabatic system for which dQ = 0. However, if we treat the universe as an open thermodynamic system allowing irreversible matter creation from the energy of the gravitational field, we can account for entropy production right from the beginning, and the second law of thermodynamics is also incorporated into the evolutionary equations in a more meaningful way. Thus, the creation of matter acts as a source of internal energy. In such situation the number of particles N in a given volume V is not to be a constant but is time dependent. Therefore, Eq. (1) modifies as

$$d ( \rho _ { m } V ) = d Q - p _ { m } d V + ( h / n ) \ d ( n V ), \quad \quad \quad$$

where N = nV , n is the particle number density, and h = ( ρ m + p m ) is the enthalpy per unit volume of the system. In case of adiabatic system where dQ = 0, Eq. (4) for an open system reduces to

$$d ( \rho _ { m } V ) + p _ { m } d V = ( h / n ) \ d ( n V ).$$

We see that in such a system the thermal energy is received due to the change of the number of particles. In cosmology, this change may be considered as a transformation of energy from gravitational field to the matter.

In the context of an open system, Eq.(5) can be rewritten as

$$d ( \rho _ { m } V ) = - ( p _ { m } + p _ { c } ) d V,$$

where

$$p _ { c } = - ( h / n ) ( d N / d V ). \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

Equation (6) suggests that the creation of matter in an open thermodynamic system corresponds to a supplementary pressure p c , which must be considered as a part of the cosmological pressure entering into the Einstein field equations (decaying of matter leads to a positive decay pressure) and is equivalent to adding the term p c given by (7) to the thermodynamic pressure p m . It is to be noted that p c is negative or zero depending on the presence or absence of particle creation.

Since the increment in entropy for an adiabatic system is only caused by creation of matter, the entropy is an extensive property of the system. In present scenario, S is proportional to the number of particles included in the system. Therefore, the entropy change dS for open systems from (2) and (5) becomes

$$T d S = \left ( h / n \right ) \, d ( n V ) - \mu \, d ( n V ) = ( T S / N ) d N \Rightarrow \frac { d S } { S } = \frac { d N } { N },$$

where µ is the chemical potential given by µ = ( h -Ts /n ) , s = S/V . Since the second law of thermodynamics is a fundamental law in physics, the presence or absence of particle creation can not affect it. This law basically requires dS ≥ 0, consequently, Eq. (8) gives

$$d N \geq 0.$$

The above inequality implies that the space-time can produce matter, whereas the reverse process is thermodynamically not admissible.

The basic idea of this entire formulation is to modify the usual energy momentum conservation law in an open thermodynamical system, which leads to the explicit use of a balance equation for the number density of the particles created, in addition to Einstein's field equations.

The particle flux vector is given by

$$N ^ { \alpha } = n u ^ { \alpha },$$

and N α is assumed to satisfy the balance equation 44 79 ,

$$N ^ { \alpha } _ { \ \ ; \alpha } = \psi,$$

where the function ψ denotes a source term of particle creation which is positive or negative depending on whether there is production or annihilation of particles. In standard cosmology ψ is usually assumed to be zero.

## 3. A brief review of Modified f ( R,T ) gravity theory

The f ( R,T ) theory is a modified theory of gravity, in which the Einstein-Hilbert Lagrangian is modified by replacing Ricci scalar curvature R by an arbitrary function of R and trace T of energy-momentum tensor, i.e., f ( R,T ). So the gravitational action for f ( R,T ) modified theory of gravity 58 in the units G = 1 = c takes the following form.

$$S = \frac { 1 } { 8 \pi } \int d ^ { 4 } x \sqrt { - g } \left [ \frac { f ( R, T ) } { 2 } + \mathcal { L } _ { m } \right ] \,,$$

where g is the determinant of the metric tensor g µν and L m corresponds to matter Lagrangian. The energy-momentum tensor T µν , defined by fluid Lagrangian density is given by

$$\ u v o n & \ u \mu \nu, \ u c m u \, \upsilon y \, \ m u \, \text{language} \, \ u \theta \, \gamma \, \upsilon \, \upsilon y \\ & T _ { \mu \nu } = - \frac { 2 } { \sqrt { - g } } \frac { \delta ( \sqrt { - g } \, \mathcal { L } _ { m } ) } { \delta g ^ { \mu \nu } }, \\ \nu T _ { \dots } \, \text{$Rv$$$

and its trace, T = g µν T µν . By assuming that the matter Lagrangian density L m depends only on the metric tensor components g µν , not on its derivatives, we obtain

$$T _ { \mu \nu } = g _ { \mu \nu } \mathcal { L } _ { m } - 2 \frac { \partial \mathcal { L } _ { m } } { \partial g ^ { \mu \nu } }.$$

The equations of motion by varying the action (12) with respect to metric tensor become

$$f _ { R } ( R, T ) R _ { \mu \nu } - \frac { 1 } { 2 } f ( R, T ) g _ { \mu \nu } + ( g _ { \mu \nu } \Box - \nabla _ { \mu } \nabla _ { \nu } ) f _ { R } ( R, T ) & = 8 \pi T _ { \mu \nu } - f _ { T } ( R, T ) T _ { \mu \nu } - f _ { T } ( R, T ) \odot _ { \mu \nu }, \\ \quad \quad \quad$$

where f R and f T denote the derivatives of f ( R,T ) with respect to R and T , respectively. Here, ∇ µ is covariant derivative and /square ≡ ∇ ∇ µ µ is the d'Alembert operator and /circleminus µν is defined by

$$\odot _ { \mu \nu } \equiv g ^ { \alpha \beta } \frac { \delta T _ { \alpha \beta } } { \delta g ^ { \mu \nu } }.$$

Using (14) into (16), we get

$$\odot _ { \mu \nu } = - 2 T _ { \mu \nu } + g _ { \mu \nu } \mathcal { L } _ { m } - 2 g ^ { \alpha \beta } \frac { \partial ^ { 2 } \mathcal { L } _ { m } } { \partial g ^ { \mu \nu } \partial g ^ { \alpha \beta } }.$$

Since the field equations of f ( R,T ) gravity depend on /circleminus µν , i.e., on the physical nature of the matter, several models corresponding to different form of f ( R,T ) would generate depending on the nature of the matter source. The gravitational field equations can be recast in such a form that the higher order corrections coming both from the geometry , and from matter -geometry coupling , provide a stress-energy tensor of geometrical and matter origin, describing an effective source term on right hand side of the standard field equations. The equations of f ( R,T ) gravity are much more complicated with respect to the ones of General Relativity even for FRW metric. For this reason many possible form of f ( R,T ), for example, f ( R,T ) = R + 2 f ( T ), f ( R,T )= µf 1 ( R ) + νf 2 ( T ), where f 1 ( R ) and f 2 ( T ) are arbitrary functions of R and T , and µ and ν are real constants, respectively 58 -62 , and f ( R,T ) = R f T ( ) 68 , etc., may be considered to solve the field equations. We however do not consider the general case, but restrict ourselves to the following form

$$f ( R, T ) = R + 2 f ( T ),$$

where f ( T ) is an arbitrary function of the trace of energy-momentum tensor of matter. The term 2 f ( T ) in the gravitational action modified the gravitational interaction between matter and curvature.

Using (18), one can re-write the gravitational field equations defined in (15) as

$$R _ { \mu \nu } - \frac { 1 } { 2 } R g _ { \mu \nu } = 8 \pi T _ { \mu \nu } - 2 ( T _ { \mu \nu } + \odot _ { \mu \nu } ) f ^ { \prime } ( T ) + f ( T ) g _ { \mu \nu }.$$

Here, a prime stands for derivative of f ( T ) with respect to T .

The main issue now arises on the matter contents of the universe through the energy momentum tensor and consequently on the matter Lagrangian L m and the trace of the energy momentum tensor. We assume the universe is filled with a perfect fluid which is incorporated in the next section.

## 4. Model and field equations

We consider a homogenous and isotropic universe with spatially flat geometry described by flat Friedmann-Robertson-Walker (FRW) metric

$$d s ^ { 2 } = d t ^ { 2 } - a ^ { 2 } ( t ) ( d x ^ { 2 } + d y ^ { 2 } + d z ^ { 2 } ),$$

where a t ( ) is the scale factor, which is a function of cosmic time t only.

In the formalism of particle creation, the second law of thermodynamics leads naturally to a modification of the energy momentum tensor with an additional creation pressure depends on the rate of creation of particles. In the presence of particle creation, the energy momentum tensor of perfect fluid is given by

$$T _ { \mu \nu } & = ( \rho _ { m } + p _ { m } + p _ { c } ) u _ { \mu } u _ { \nu } - ( p _ { m } + p _ { c } ) g _ { \mu \nu },$$

where u µ is the four velocity of the fluid such that u u µ ν = 1, and in comoving coordinates u µ = δ µ 0 .

The trace of energy momentum tensor (21), gives

$$T = \rho _ { m } - 3 ( p _ { m } + p _ { c } ).$$

We treat the scalar invariant L m as the effective pressure of the perfect fluid matter and pressure originated by creation of particles. Therefore, the matter Lagrangian may be assumed as L m = -( p m + p c ), therefore, Eq. (17) becomes

$$\odot _ { \mu \nu } = - 2 T _ { \mu \nu } - ( p _ { m } + p _ { c } ) g _ { \mu \nu }.$$

In view of (23), the field equations (19) can be rewritten as

$$R _ { \mu \nu } - \frac { 1 } { 2 } g _ { \mu \nu } R = 8 \pi T _ { \mu \nu } + 2 \left [ T _ { \mu \nu } + g _ { \mu \nu } ( p _ { m } + p _ { c } ) \right ] f ^ { \prime } ( T ) + g _ { \mu \nu } f ( T ).$$

The field equations (24) for a fluid endowed with matter creation (21) in the background of metric (20), yield

$$3 H ^ { 2 } \ = \ 8 \pi \rho _ { m } + 2 ( \rho _ { m } + p _ { m } + p _ { c } ) f ^ { \prime } ( T ) + f ( T ),$$

$$2 \dot { H } + 3 H ^ { 2 } \ = \ - 8 \pi ( p _ { m } + p _ { c } ) + f ( T ),$$

where H t ( ) = ˙ a/a is the Hubble parameter. A dot denotes derivative with respect to cosmic time t .

For adiabatic matter creation, the pressure p c in Eq. (7) takes the following form 44

$$p _ { c } = - \frac { ( \rho _ { m } + p _ { m } ) } { 3 n H } \ \psi.$$

For a complete formulation of the problem, we consider one more relation between the particle number n and V = a 3 , describing the dynamics of n as a result of matter creation (decay) processes. This relation is given by Eq. (11), which for the metric (20) takes the form

$$\dot { n } + 3 n H = \psi,$$

## 5. Solution of field equations

The field equations (25)-(28) have seven unknowns, namely, a , ρ m , p m , p c , n f T , ( ) and ψ t ( ). We have only three equations plus the constraint (27). Therefore, one needs three more relations in order to construct a definite cosmological scenario.

In first choice, we consider a particular function given by 58

$$f ( T ) = \lambda T, \quad ( \lambda \ \text{ is a constant} ).$$

With this assumption, the field equations (25) and (26) yield

$$3 H ^ { 2 } \ & = \ ( 8 \pi + 3 \lambda ) \rho _ { m } - \lambda ( p _ { m } + p _ { c } ), & ( 3 0 ) \\ \dot { H } \perp \dot { H } ^ { 2 } \ & = \ - \ ( 8 \pi \perp \dot { H } \ \lambda \ / \ n \ \perp \ n \ \cdot \perp \lambda \ \alpha & ( \dot { H } \dot { H } \dot { H } )$$

$$2 \dot { H } + 3 H ^ { 2 } \ = \ - ( 8 \pi + 3 \lambda ) ( p _ { m } + p _ { c } ) + \lambda \rho _ { m }.$$

In order to obtain the exact solution of the field equations, we assume two more additional relations: the equation of state and the matter creation rate. In the cosmological domain, the former is usually given by the gamma-law

$$p _ { m } = ( \gamma - 1 ) \rho _ { m },$$

where γ is a constant lies in the interval [0, 2] and known as the EoS parameter of the perfect fluid.

Using (32) into (30) and (31), and simplifying, we get a single evolution equation for H :

$$2 \dot { H } + ( 8 \pi + 2 \lambda ) ( \gamma \rho _ { m } + p _ { c } ) = 0.$$

We confine our attention to the simple phenomenological expression for the matter creation rate 21

$$\psi ( t ) = 3 \beta n H,$$

where the parameter β lies in the interval (0 , 1) and is assumed to be a constant.

Using (32) and (34) into (27), we have

$$p _ { c } = - \beta \gamma \rho _ { m },$$

Putting (32) and (35) into (30), we obtain

$$\rho _ { m } & = \frac { 3 H ^ { 2 } } { 8 \pi + 4 \lambda - \gamma \lambda ( 1 - \beta ) }. \\ \intertext { t o o \cdot \cdot \cdot \cdot \cdot } _ { t o o \cdot \cdot \cdot \cdot \cdot }.$$

Substituting (35) and (36) into (33) gives

$$\dot { H } + \frac { 3 } { 2 } \frac { \gamma ( 8 \pi + 2 \lambda ) ( 1 - \beta ) } { [ 8 \pi + 4 \lambda - \gamma \lambda ( 1 - \beta ) ] } \, H ^ { 2 } = 0. \\ \intertext { a n n e } \, \frac { \lambda } { \lambda } \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots$$

The solution of (37) for γ = 0 and for all values of λ and β is given by

/negationslash

$$H ( t ) = \left ( C + \frac { 3 } { 2 } \frac { \gamma ( 8 \pi + 2 \lambda ) ( 1 - \beta ) } { [ 8 \pi + 4 \lambda - ( 1 - \beta ) \gamma \lambda ] } \, t \right ) ^ { - 1 }, \\ \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \end{array} \,,$$

where C is an integration constant. For γ = 0, the well known de-Sitter scale factor a t ( ) = a 0 e H t 0 is obtained.

From Eq. (38) we find the following expression for the scale factor

$$a ( t ) = D \left ( C + \frac { 3 } { 2 } A \gamma \, t \right ) ^ { \frac { 2 } { 3 A \gamma } },$$

where D is a new integration constant and A = (8 π +2 )(1 λ -β / π ) 8 +4 λ -γλ (1 -β ).

The above scale factor may be rewritten as

$$a ( t ) = a _ { 0 } \left ( 1 + \frac { 3 } { 2 } A \gamma H _ { 0 } ( t - t _ { 0 } ) \right ) ^ { \frac { 2 } { 3 A \gamma } },$$

where H = H &gt; 0 0 at t = t 0 . The subscript '0' refers to the present values of parameters. Since 0 ≤ γ ≤ 2, we must have A &gt; 0 for expansion of the universe. Also, A &gt; 0 implies λ &gt; 0 as 0 ≤ β &lt; 1.

If γ &lt; 0 then we have a Big Rip singularity at a finite value of cosmic time t br = t -t 0 = -2 3 / H Aγ 0 . Thus, we have a phantom cosmology for γ &lt; 0. If one choose t 0 = 2 H -1 0 / Aγ 3 then (40) takes the familiar form of power-law expansion of the universe, i.e.,

$$a ( t ) = a _ { 0 } \left ( \frac { 3 } { 2 } A \gamma H _ { 0 } \, t \right ) ^ { \frac { 2 } { 3 A \gamma } }.$$

If λ = 0 = β , (39) and (41) reduce to the well-known expressions of power-law expansion of scale factor for a flat FRW model in GR.

By use of (40) one may express the energy density of matter, particle creation pressure and the particle number density as functions of the scale factor a . These parameters respectively have the following forms

$$\rho _ { m } \ = \ \rho _ { 0 } \left ( \frac { a _ { 0 } } { a } \right ) ^ { 3 A \gamma }, \\ \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \ \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \ \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \,$$

$$n \ = \ n _ { 0 } \left ( \frac { a _ { 0 } } { a } \right ) ^ { \dots } \stackrel { \cdot } { 3 ( 1 - \beta ) },$$

$$p _ { c } \ = \ \stackrel { \theta \nu } { - \beta \gamma \rho _ { 0 } \left ( \frac { a _ { 0 } } { a } \right ) ^ { 3 A \gamma } }, \\ \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

where ρ 0 = 3 H / 2 0 [8 π +4 λ -γλ (1 -β )] is the present value of energy density. Here, n 0 is the present value of particle number density for any values of β . The above results show that the transition from one phase to other phase, in the course of expansion, happens exactly as in the standard cosmological model.

The number of particles N in a given volume V is given by

$$N = N _ { 0 } \left ( \frac { a } { a _ { 0 } } \right ) ^ { 3 \beta },$$

which shows that N increases with time. If β = 0, N would remain constant throughout the evolution of the universe and we would recover the standard FRW model of the universe in f ( R,T ) theory. Again, from (8) S = S 0 ( N/N 0 ), the entropy in terms of scale factor is

$$S = S _ { 0 } \left ( \frac { a } { a _ { 0 } } \right ) ^ { 3 \beta }.$$

The deceleration parameter defined as q = -aa/a ¨ ˙ 2 , gives

$$q = - 1 + \frac { 3 \gamma A } { 2 } = \left [ \frac { 3 \gamma } { 2 } \frac { ( 8 \pi + 2 \lambda ) ( 1 - \beta ) } { [ 8 \pi + 4 \lambda - ( 1 - \beta ) \gamma \lambda ] } - 1 \right ] \\. \. \. \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdots$$

which shows that q is independent of cosmic time t . Therefore, q may be positive or negative for a given set of values of β and λ . We know that the model accelerates for q &lt; 0, therefore, the value of A must be 0 &lt; A &lt; 2 3 / γ for an accelerated universe.

As expected, the above solutions reduce to the standard FRW model of GR for β = 0 and λ = 0 and for all values of γ . In what follows, we study the role of f ( R,T ) gravity and particle creation in early and late time evolution of the universe.

$$\text{Case $I$} \colon \gamma = \frac { 2 } { 3 }$$

Fig. 1(a) plots the scale factor as a function of time for γ = 2 3 and some selected values / of λ and β . We observe that if β = 0, q &lt; 0 for all λ &gt; 0, therefore, we find that the universe accelerates in f ( R,T ) gravity without particle creation. Similarly, if λ = 0, i.e., in the absence of f ( T ), q = -β &lt; 0 for any values of β &gt; 0. Thus, the acceleration occurs due to particle creation. The rate of expansion increases more rapidly for non-zero values of β and λ . It is to be noted that if λ = 0 = β then the the marginal inflationary phase of GR is recovered, i.e., a ∼ t and q = 0.

Fig. 1(a). Scale factor as a function of time for γ = 2 3 and some selected values of λ and β .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0633v2-figure-bd60dda6-0395-492c-b877-2457c5918e39.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0633v2-figure-7651810d-e89d-491d-83d7-414ee82f589a.jpg)

Fig.1(b). Scale factor as a function of time for γ = 4 3 and some selected values of λ and β .

Case II: γ = 4 3

/negationslash

In this case, if β = 0 and λ &gt; 0, we have q &gt; 0. This shows that the universe decelerates in the absence of particle creation. If λ = 0 then q ≥ 0 for 0 &lt; β ≤ 1 2, and / q &lt; 0 for 1 / 2 &lt; β &lt; 1. Therefore, the universe decelerates or accelerates depending on the rate of particle creation. However, if λ = 0 and β = 0, the deceleration or acceleration of the universe depend on the following constrains, respectively:

/negationslash

$$0 & < \beta \leq \frac { 1 } { 4 }, \ \lambda > 0 \ \ o r \ \frac { 1 } { 4 } < \beta < \frac { 1 } { 2 }, \ 0 < \lambda < \frac { 6 \pi - 1 2 \pi \beta } { - 1 + 4 \beta }, \\ \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot$$

$$\frac { 1 } { 4 } < \beta \leq \frac { 1 } { 2 }, \, \lambda > \frac { 6 \pi - 1 2 \pi \beta } { - 1 + 4 \beta } \ o r \ \frac { 1 } { 2 } < \beta < 1, \, \lambda > 0. \\ \text{savior of scale factor versus time is shown in Fio} \ \ 1 ( h ) \text{ for some selected values of}$$

The behavior of scale factor versus time is shown in Fig. 1(b) for some selected values of λ and β . The figure shows that the universe accelerates fast due to higher particle creation rate. For λ = 0 = β , a ∼ t 1 2 / and q = 1, which is the radiation phase in GR.

## Case III: γ = 1

/negationslash

In this case, the universe expands with decelerated rate as q &gt; 0 for β = 0 and λ &gt; 0. Fig. 1(c) plots graph between scale factor versus time for some selected values of λ and β . For λ = 0, we have q &gt; 0 for 0 &lt; β &lt; 1 3, and / q &lt; 0 for 1 / 3 &lt; β &lt; 1. The critical case ( β = 1 3, / q = 0), describes a coasting cosmology. For λ = 0 and β = 0, the model decelerates or accelerates under the following constraints:

/negationslash

$$0 < \beta < \frac { 1 } { 3 }, \, 0 < \lambda < \frac { \pi - 3 \pi \beta } { \beta },$$

$$0 < \beta \leq \frac { 1 } { 3 }, \ \lambda > \frac { \pi - 3 \pi \beta } { \beta } \ o r \ \frac { 1 } { 3 } < \beta < 1, \ \lambda > 0,$$

respectively. As expected, for λ = 0 = β , we have a ∼ t 2 3 / and q = 1 2, i.e., / the model reduces to standard matter-dominated era of GR.

Fig. 1(c). Scale factor as a function of time for γ = 1 and some selected values of λ and β .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0633v2-figure-5f177381-9718-448a-a301-a9d74543e114.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0633v2-figure-a9d11d6e-ab37-4d41-b166-683e5a3d36ff.jpg)

Fig. 1(d). Scale factor as a function of time for γ = 1 2 and some selected values of λ and β .

$$\text{Case IV} \colon \gamma = \frac { 1 } { 2 }$$

In this case, if λ = 0 = β , a ∼ t 4 3 / and q = -1 4, / which corresponds to the present accelerated phase of the universe of the standard FRW universe in GR. Since the universe accelerates even in absence of both f ( T ) and particle creation, therefore, the contribution of f ( R,T ) gravity or particle creation just enhance the rate of acceleration of the universe. Fig. 1(d) plots the dynamics of scale factor versus t , which is similar to case I.

## 6. Kinematic tests

Now, we derive some kinematical relations of the model as proposed in the preceding sections.

## 6.1 The density parameter

The density parameter, defined as Ω m = ρ m /ρ c , where ρ c = 3 H / π 2 8 , is given by

$$\Omega _ { m } & = \frac { 8 \pi } { 8 \pi + 4 \lambda - ( 1 - \beta ) \gamma \lambda }. \\ \text{that to } \mathcal { O } \ \ < 1 \text{ for all values of } \mathcal { O } < \sim \ < 2 \text{ } \mathcal { O } \times \mathcal { O } \times 1 \text{ and } \mathcal { U } \searrow \mathcal { H } \text{one}.$$

Therefore, it is clear that Ω m &lt; 1 for all values of 0 ≤ γ ≤ 2, 0 &lt; β &lt; 1 and λ &gt; 0. Hence the universe is negatively curved. In absence of both λ and β , we have Ω m = 1 for all γ , i.e., the flat model of GR is recovered.

## 6.2 Lookback time-redshift

The lookback time ∆ t = t 0 -t ( z ), is the difference between the age of the universe at the present time z = 0 and the age of the universe when a particular right ray at redshift z was emitted.

For a given redshift z , the scale factor a t ( z ) is related to a 0 by

$$a ( t _ { z } ) = a _ { 0 } ( 1 + z ) ^ { - 1 }.$$

From (41) and (53), the cosmic time in terms of redshift is given by

$$t ( z ) = \frac { 2 H _ { 0 } ^ { - 1 } } { 3 \gamma A } ( 1 + z ) ^ { - \frac { 3 \gamma A } { 2 } },$$

where H 0 is expressed in the usual practical observational units of km s -1 Mpc -1 and its value is believing to be somewhere between 50-80 km s -1 Mpc -1 . However, H 0 is dimensionally similar to the reciprocal of time. The reciprocal of Hubble constant is called the Hubble time t H : t H = H -1 0 , where t H is expressed in s and H 0 in s -1 . If H 0 is expressed in km s -1 Mpc -1 and t H in gigayears (1 gr = 1 milion years = 10 9 years ) then t H = 977 8 . /H 0 .

Therefore, from (54) we have

$$t _ { 0 } - t ( z ) = \frac { 2 H _ { 0 } ^ { - 1 } } { 3 \gamma A } \left [ 1 - ( 1 + z ) ^ { - \frac { 3 \gamma A } { 2 } } \right ]. \quad \quad \quad$$

Fig. 2 plots lookback time versus redshift for γ = 1 and some selected values of λ and β . All models coincide for lower redshift since they follow the same behavior. The graph shows that the lookback time increases for higher values of β . Thus, the universes with larger matter creation rate are older.

For small values of redshift, (55) becomes

$$H _ { 0 } \left ( t _ { 0 } - t ( z ) \right ) = z - \left ( 1 + \frac { 3 \gamma A } { 2 } \right ) z ^ { 2 } + \cdots.$$

Taking lim z →∞ in (55), the present age of the universe is

$$t _ { 0 } = \frac { 2 H _ { 0 } ^ { - 1 } } { 3 \gamma A } = \frac { H _ { 0 } ^ { - 1 } } { 1 + q }.$$

Thus, the age of the universe depends on both parameters β and λ .

Fig. 2. Lookback time versus redshift for γ = 1, H 0 = 60 and some selected values of λ and β .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0633v2-figure-1eedcd36-670e-46af-b01c-92479e8382a0.jpg)

Fig. 3. Proper distance versus redshift for γ = 1, H 0 = 60 and some selected values of λ and β .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0633v2-figure-aa646556-35a7-45b4-ab38-002d5cbf5a5e.jpg)

## 6.3 Proper distance-redshift

The proper distance between the source and observer is defined as d z ( ) = a r z 0 ( ), where r z ( ) is the radial distance of the object at light emission in term of redshift given by

$$r ( z ) = \int _ { t } ^ { t _ { 0 } } \frac { d t } { a ( t ) } = \frac { H _ { 0 } ^ { - 1 } } { a _ { 0 } \left ( \frac { 3 \gamma A } { 2 } - 1 \right$$

Consequently, the proper distance becomes

$$d ( z ) & = \frac { H _ { 0 } ^ { - 1 } } { \left ( \frac { 3 \gamma A } { 2 } - 1 \right ) } \left [ 1 - ( 1 + z ) ^ { - ( \frac { 3 \gamma A } { 2 } - 1 ) } \right ]. \\ \cdot. & \quad \, \quad \, \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \,.$$

The proper distance as a function of redshift for some selected values of β and λ are displayed in Fig. 3. We observe that the f ( T ) contribution in f ( R,T ) and particle creation gives rise to proper distance.

Equation (59) can be rewritten as

$$H _ { 0 } d ( z ) = z - \frac { 3 \gamma A } { 4 } z ^ { 2 } + \cdots.$$

From (59), it is observed that the distance dz is maximum at z →∞ . Hence,

$$H _ { 0 } d ( z \rightarrow \infty ) = \frac { 1 } { \frac { 3 \gamma A } { 2 } - 1 } = \frac { 1 } { q }.$$

## 6.4 Luminosity distance-redshift

The best-known way to trace the evolution of the universe observationally is to look into the redshift-luminosity distance relation. The luminosity distance d l is defined by the relation

d 2 l = l 4 πL , where l is the luminosity of the object and L is the measured flux from the object. In standard FRW cosmology it is defined in terms of redshift as

$$d _ { l } = a _ { 0 } ( 1 + z ) r ( z ) = ( 1 + z ) d ( z ).$$

From (59) and (62), we get

$$d _ { l } H _ { 0 } = \frac { 1 } { \left ( \frac { 3 \gamma A } { 2 } - 1 \right ) } \left [ ( 1 + z ) - ( 1 + z ) ^ { - ( \frac { 3 \gamma A } { 2 } - 2 ) } \right ]. \quad \quad \quad ( 6 3 ) \\ \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \,.$$

The graph between Luminosity distance and redshift for some selected values of β and λ is plotted in Fig. 4. One may observe that the luminosity distance corresponding to any specific value of redshift rises due to f ( R,T ) gravity and particle creation.

Expanding (63) for small z, we find

$$H _ { 0 } d _ { l } = z - \frac { 1 } { 2 } \left ( \frac { 3 \gamma A } { 2 } - 2 \right ) z ^ { 2 } + \cdots.$$

As expected, we find the same behavior for different models at z /lessmuch 1 and the possible difference in behaviors for different models come at large redshift ( z /greatermuch 1). In Fig. 4 we observe that all curves start off with the linear Hubble law ( z = d H l 0 ) for small z , but then, only the curve for q = 1, i.e., β = 0 = λ stays linear all the way. We also note that for the small redshift the luminosity distance is larger for lower values of q.

Fig. 4. Luminosity distance versus redshift for γ = 1, H 0 = 60 and some selected values of λ and β .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0633v2-figure-487a16e3-d0bd-433c-b908-0e427d81fb43.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0633v2-figure-ec7fb51d-2d83-4a7c-b736-390f84febec2.jpg)

Fig. 5. Angular diameter distance versus redshift for γ = 1, H 0 = 60 and some selected values of λ and β .

## 6.5 Angular diameter distance-redshift

The angular diameter distance d A is the ratio of physical transverse size of an object to its angular size (in radians). In terms of z , it is given by

$$d _ { A } = \frac { d ( z ) } { 1 + z } = \frac { d _ { l } } { ( 1 + z ) ^ { 2 } }.$$

Using (59), we have

$$H _ { 0 } d _ { A } & = \frac { 1 } { \left ( \frac { 3 \gamma A } { 2 } - 1 \right ) } \left [ \left ( 1 + z \right ) ^ { - 1 } - \left ( 1 + z \right ) ^ { - \frac { 3 \gamma A } { 2 } } \right ]. \\ \intertext { m a t t h a m a t t h a m a t t h a m a t t h a m a t t h a m a t t h a m a t t h a m a t t h a m a t t h a m a t t h a m a t t h a m a t t h a m a t t h a m a t t h a m a t t h a m a t t h a m a t t h a m a t t h a m a t t h a m a t t h a m a t t h a m a t t h a m a t t h a m a t t h a m a t t h a m a t t h a m a t t h a m a t t h a m a t t h a m a t t h a m a t t h a m a t t h a m a t t h a m a t t h a m b a t t h a m a t t h a m a t t h a m b a t t h a m a t t h a m b a t t h a m b a t t h a m b a t t h a m b a t t h a m b a t t h a m b a t t h a m b a t t h a m b a t t h a m b a t t h a m b a t t h a m b a t t h a m b a t t h a m b a t t h a m b a t t h a m b a t t h a m b a t t h a m b a t t h a m b a t t h a m b a t t h a m b a t t h a m b a t t h a m b a t t h a m b a t t h a m b a t t h a m b a t t h a m b a t t h a m b a t t h a m b a t t h a m b a t t a m b a t t h a m b a t t h a m b a t t h a m b a t t h a m b a t t h a m b a t t h a m b a t t h a m b a t t h a m b a t t h a m b a t t h a m b a t t h a m b a t t h a m b a t t h a m b a t t h a m b a t t h a m b a t t h a m b a t t h a m b a t t h a m b a t t h a m b a t t h a m b a t t h a m b a t t h a m b a t t h a m b a t t h a m b a t t h a m b a t t h a m b a t t h a m b a t t h a m b b a t t h a m b a t t h a m b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b a t t h a m b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b b a t t h a m b b a t t h a m b b a t t h a a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b b a t t h a m b b a t t h a m b b a t t h a mb b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a mb b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a a t t h a m b b a t t h a m b b a t t h a m b b a a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a mb b a t t h a m b b a t t h a m b b a t t h a m b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b a t t h a m b b a t t h a m b b a t t h a m b a t t h a m b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b b a t t h a m b a$$

In Fig. 5 we plot the angular diameter distance versus redshift for some selected values of β and λ . The graph shows that the f ( R,T ) and particle creation enhance the angular distance. The angular diameter distance initially increases with increasing z and eventually begins to decrease.

On expanding (66), we get

$$H _ { 0 } d _ { A } = z + \left [ 1 - \frac { \left ( \frac { 3 \gamma A } { 2 } + 1 \right ) \left ( \frac { 3 \gamma A } { 2 } + 2 \right ) } { 2 \left ( \frac { 3 \gamma A } { 2 } - 1 \right ) } \right ] z ^ { 2 } + \cdots. \\ \intertext { h i o n }$$

## 7 Conclusion

In this paper, we have studied a flat FRW cosmological model described by an open thermodynamic system including particle creation at the expense of gravitational field in f ( R,T ) theory of gravity. We have obtained exact solutions for the scale factor and various physical quantities by assuming a suitable form of f ( R,T ) = R +2 ( f T ) and 'gamma-law' equation of state. The model exhibits non-singular power-law expansion of the universe for 0 ≤ γ ≤ 2. The model shows Big Rip singularity for γ &lt; 0. The dynamics of the scale factor and other physical quantities have been examined through some graphical representations in various phases of evolution of the universe with some selected values of λ and β .

It has been observed that the scale factor always increases with decelerated and accelerated rate depending upon the contribution of particle creation and the parameter λ . It has been noticed that the universe expands with slow rate in early time but fast during late time. The energy density and effective pressure always decrease with time and both tend to zero in late time for 0 ≤ γ ≤ 2. The number of particles increase with time in all phases. The number of particles in the absence of particle production remain constant throughout the evolution of the universe, which is quit obvious. The deceleration parameter has been found as a constant, which exhibits both decelerated and accelerated universe under some constraints on different parameters. The density parameter shows that the model becomes open in the presence of particle creation in f ( R,T ) theory.

We have also discussed some observational consequences of the model through some kinematics tests such as lookback time, proper distance, luminosity distance and angular diameter distance with respect to redshift. The results for the cosmological tests are found to be compatible with the present observations. These tests are found to be depending on λ and β . The universe with particle creation is always older than the usual FRW model in absence of particle creation. The model of Lima et al. 21 may be recovered for λ = 0.

In summary, we have studied a cosmological model with particle creation in f ( R,T ) gravity theory to understand the current acceleration expansion of the universe. We have found that the negative pressure due to the matter creation may play the role of dark energy and derive the accelerating expansion of the universe in f ( R,T ) theory. We may expect that the

process of particle creation is also an ingredient which accounts this unexpected observational results. The changes introduced by the particle creation process, which is quantified by the parameter β , provide reasonable observational results. The new fact justifying the present work is that we have considered the thermodynamics approach for which particle creation is at the expense of the gravitational field. A general expression relating the energy densities and particle number density as function of scale factor have been established. One may find that the particle creation changes the predictions of standard cosmology, thereby alleviating the problem of reconciling observations with the inflationary scenario. In future work, we plan to constraint cosmological model with matter creation using complimentary astronomical observations.

Acknowledgement VS expresses his sincere thank to University Grants Commission (UGC), India, for providing Senior Research Fellowship (SRF).

## References

- 1. A.G. Riess et al., Astrophys. J. 116 , 1009 (1998).
- 2. S. Perlmutter et al., Astrophys. J. 517 , 565 (1999).
- 3. R.G. Vishwakarma, Mon. Not. Roy. Astron. Soc. 331 , 776 (2002).
- 4. C.B. Netterfield et al., Astrophys. J. 571 , 604 (2002).
- 5. D.N. Spergel et al., Astro J. Suppl. 148 , 175 (2003) [astro-ph/0302209].
- 6. V. Sahni, Lec. Notes Phys. 653 , 141 (2004) [astro-ph/0403324].
- 7. J. Frieman et al., Ann. Rev. Astron. Astrophys. 46 , 385 (2008) [astro-ph/0803.0982].
- 8. R.R. Caldwell and M. Kamionkowski, Ann. Rev. Nucl. Part. Sci. 59 , 397 (2008) [astroph/0903.0866].
- 9. L.M. Krauss and M. Turner, Gen. Rel. Grav. 27 , 1137 (1995).
- 10. M.S. Turner and M. White, Phys. Rev. D 56 , R4439 (1997).
- 11. T. Chiba, N. Sugiyama and T. Nakamura, Mon. Not. Roy. Astron. Soc. 289 , L5 (1997).
- 12. T. Padmanabhan Phys. Rept. 380 , 235 (2003) [hep-th/0212290].
- 13. J. Martin, Mod. Phys. Lett. A 23 , 1252 (2008) [astro-ph/0803.4076].
- 14. U. Alam et al., Mon. Not. Roy. Astron. Soc. 354 , 275 (2004)
- 15. T. Chiba et al., Phys. Rev. D 62 , 023511 (2000).
- 16. T. Padmanabhan and T.R. Chaudhury, Phys. Rev. D 66 , 081301 (2002).
- 17. M.C. Bento et al., Phys. Rev. D 66 , 043507 (2002).
- 18. L.R.W. Abramo and J.A.S. Lima, Class. Quant. Grav. 13 , 2953 (1996).
- 19. W. Zimdhal and P. Pav´n, o Mon. Not. Roy. Astron. Soc. 266 , 872 (1994).
- 20. J.A.S. Lima and L.R.W. Abramo, Phys. Lett. A 257 , 123 (1999).
- 21. J.A.S. Lima, A.S. Germano and L.R.W. Abramo, Phys. Rev. D 53 , 4287 (1996).

- 22. I. Brevik and G. Stokkan, Astrophys. Space. Sci. 239 , 89 (1996).
- 23. A.H. Guth, Phys. Rev. D 23 , 347 (1981)
- 24. A.D. Linde, Phys. Lett. 108B , 389 (1981)
- 25. W.L. Freedman, A. Olinto, J. Friemaan and D. Schramm, Proceedings of the 18th texas Symposium on Relativistic Astrophys. World Scietific (1998) .
- 26. C. Brans and R.H. Dicke, Phys. Rev. D 124 , 925 (1961).
- 27. P.G. Bergmann, Int. J. Theor. Phys. 1 , 25 (1968).
- 28. K.Nortvedt, Astrophys. J. 161 , 1059 (1970).
- 29. R.V. Wagoner, Phys. Rev. D 1 , 3209 (1970).
- 30. P.A.M. Dirac, Nature 139 , 323 (1937).
- 31. V. Canuto, P.J. Adoms, S.H. Hsieh and E. Tsiang, Nature 261 , 438 (1976).
- 32. V. Canuto, P.J. Adoms, S.H. Hsieh and E. Tsiang, Phys. Rev. D 16 , 1643 (1977).
- 33. H. Bondi and T. Gold, Mon. Not. Roy. Astron. Soc. 108 , 252 (1948).
- 34. F. Hoyle, Mon. Not. Roy. Astron. Soc. 109 , 365 (1949).
- 35. J.V. Narlikar, Nonstandard cosmologies , in Vth Brazialian school of cosmology and Gravitation, ed. M. Novello (World Scientific, Singapore) 152 (1987).
- 36. E.P. Tryon, Nature 246 , 396 (1973).
- 37. P.I. Fomin, Dokl. Akad. Nank SSSR A 9 , 831 (1975).
- 38. R. Brout, F. Englert and E. Gunzig, Ann. Phys. (NY) 115 , 78 (1978).
- 39. R. Brout, F. Englert and E. Gunzig, Gen. Rel. Grav. 1 , 1 (1979).
- 40. R. Brout, F. Englert, F. Fr´re, J.M. Gunzig, P. Nardone, P. Spindel, C. Tuffin and E. Gunzig, e Nucl. Phys. B 170 , 228 (1980).
- 41. R. Brout, F. Englert and P. Spindel, Phys. Rev. Lett. 43 , 417 (1979).
- 42. E. Gunzig, J. Gehenian and I. Prigogine, Nature 330 , 621 (1983).
- 43. I. Prigogine, J. Gohenian, E. Gunzig and P. Nardone, Proc. Natl. Acad. Sci. USA 85 , 7428 (1988).
- 44. M.O. Calv˜o, J.A.S. Lima and I. Waga, a Phys. Lett. A 162 , 223 (1992).
- 45. J.A.S. Lima and A.S. Germano, Phys. Lett. A 170 , 373 (1992).
- 46. J.A.S. Lima and J.S. Alcaniz, Astron. Astrophys. 348 , 1 (1999).
- 47. J.S. Alcaniz and J.A.S. Lima, Astron. Astrophys. 349 , 72 (1999).
- 48. W. Zimdahl et al., Phys. Rev. D 64 , 063501 (2001).
- 49. Y. Qiang, T. Jiezhang and Y. Ze-Long, Astrophys. Space Sci. 311 , 407 (2007) [arXiv,astro-ph/0503123].
- 50. C.P. Singh and A. Beesham, Astrophys. Space Sci. 336 , 469 (2011).
- 51. C.P. Singh and A. Beesham, Int. J. Theor. Phys. 51 , 3951 (2012).
- 52. C.P. Singh, Astrophys. Space Sci. 338 , 411 (2012).

- 53. C.P. Singh, Mod. Phys. Lett. A 27 , 1250070 (2012).
- 54. T. Clifton, P.G. Ferreira, A. Padilla and C. Skordis, Phys. Rep. 513 , 1 (2012) [arXiv,astroph/1106.2476].
- 55. S. Nojiri andS.D. Odintsov, Phys. Rev. D 74 , 086005 (2006).
- 56. O. Bertolami, C.G. Bochmer, T. Harko and F.S.N. Lobo, Phys. Rev. D 75 , 104016 (2007).
- 57. T. Harko, Phys. Lett. B 669 , 376 (2008).
- 58. T. Harko, F.S.N. Lobo, S. Nojiri and S.D. Odintsov, Phys. Rev. D 84 , 024020 (2011) [arXiv,grqc/1104.2669].
- 59. M. Jamil, D. Momeni, M. Raza and R. Myrzakulov, Eur. Phys. J. C 72 , 1999 (2012).
- 60. M.J.S. Houndjo, Int. J. Mod. Phys. D 21 , 1250003 (2012) [arXiv,astro-ph/1107.3887].
- 61. M.J.S. Houndjo and O.F. Piattella, Int. J. Mod. Phys. D 2 , 1250024 (2012) [arXiv,grqc/1111.4275].
- 62. M.J.S. Houndjo, C.E.M. Batista, J.P. Campos and O.F. Piattella, Canadian. J. Phys. 91 , 548 (2013) [arXiv,gr-qc/1203.6084].
- 63. M. Sharif and M. Zubair, J. Cosmo. Astropart. Phys. 21 , 28 (2012) [arXiv,gr-qc/1204.0848].
- 64. M. Jamil, D. Momeni and R. Myrzakulov, Chin. Phys. Lett. 29 , 109801 (2012).
- 65. F.G. Alvarenga, M.J. Houndjo, A.V. Monwanou and J.B. Chabi Oron, J. Mod. Phys. 4 , 130 (2013) [arXiv,gr-qc/1205.4678].
- 66. M. Sharif, S. Rani and R. Myrzakulov, Eur. Phys. J. Plus 128 , 123 (2013) [arXiv,grqc/1210.2714].
- 67. M. Sharif and M. Zubair, J. Phys. Soc. Jpn. 82 , 014002 (2013) [arXiv,gr-qc/1210.3878].
- 68. S. Chakraborty, Gen. Rel. Grav. 45 , 2039 (2013) [arXiv,gen-ph/1212.3050].
- 69. F.G. Alvarenga, A. de la Cruz-Dombriz, M.J.S. Houndjo, M. E. Rodrigues and D. S´ez-G´mez, a o Phys. Rev. D 87 , 103526 (2013) [arXiv,gr-qc/1302.1866].
- 70. A. Pasqua, S. Chattopadhyay and I. Khomenkoc, Canadian J. Phys. 91 , 632 (2013).
- 71. H. Shabani and M. Farhoudi, Phys. Rev. D 88 , 044048 (2013).
- 72. M. Sharif and M. Zubair, J. Phys. Soc. Jpn. 82 , 064001 (2013) [arXiv,gr-qc/1310.1067].
- 73. C.P. Singh and V. Singh, Gen. Rel. Grav. 46 , 1696 (2014).
- 74. K.S. Adhav, Astrphys. Space. Sci. 339 , 365 (2013).
- 75. S.D. Katore and A.Y. Shaikh, Prespacetime J. 3 , 1087 (2012).
- 76. M. Sharif and M. Zubair, J. Phys. Soc. Jpn. 81 , 114005 (2012) [arXiv,gr-qc/1301.2251].
- 77. V.U.M. Rao and D. Neelima, Astrphys. Space. Sci. 345 , 427 (2013).
- 78. N. Ahmed and A. Pradhan, Int. J. Theor. Phys. 53 , 289 (2014).
- 79. I. Prigogine, J. Gohenian, E. Gunzig and P. Nardone, Gen. Rel. Grav. 21 , 767 (1989).
C. P. Singh Vijay Singh
Friedmann cosmology with particle creation in modified $f(R,T)$ gravity
The theoretical and observational consequences of thermodynamics of open
systems, which allow particle creation are investigated in modified $f(R,T)$
($R$ is the Ricci scalar and $T$ is the trace of energy-momentum tensor) theory
of gravity within the framework of a flat Friedmann-Robertson-Walker line
element. A simplest particular model $f(R,T)=R+2f(T)$ and "gamma-law" equation
of state $p=(\gamma-1)\rho$ are assumed to explore the role of particle
creation in the early and late time evolution of the universe. A power-law
expansion model is proposed for $0\leq \gamma \leq 2$ by considering the
natural phenomenological particle creation rate $\psi=3\beta nH$, where $\beta$
is a pure number of the order of unity. A Big Rip singularity is observed for
$\gamma<0$, which describes the phantom cosmology. We observe that the
accelerated expansion of the universe is driven by the particle creation
without any exotic "quintessence" component or a cosmological constant. It is
also found by density parameter that the model becomes of negative curvature
due to particle creation and the entropy increases with the evolution of the
universe. Some kinematics tests such as lookback time, luminosity distance,
proper distance, angular diameter versus redshift are discussed in detail to
observe the role of particle creation in early and late time evolution of the
universe.