Models of Collapsing and Expanding Anisotropic Gravitating Source in f(R,T) Theory of Gravity

Models of Collapsing and Expanding Anisotropic Gravitating Source in f(R,T) Theory of Gravity

G. Abbas and Riaz Ahmed
Department of Mathematics The Islamia University
of Bahawalpur,
Bahawalpur, Pakistan.
abbasg91@yahoo.comahmedoriya@gmail.com

Abstract

In this paper, we have formulated the exact solutions of the non-static
anisotropic gravitating source in f(R,T) gravity which may lead to expansion and collapse. By assuming the no thermal conduction in gravitating source, we have determine parametric solutions in f(R,T) gravity with non-static spherical geometry filled with anisotropic fluid. We have examined the range of parameter for which expansion scalar become negative and positive leading to collapse and expansion, respectively. Further, using the definition of mass function the condition for the trapped surface have been explored and it has been investigated there exists a single horizon in this case. The impact of coupling parameter λ has been discussed in detail in both cases. For the various values of coupling parameter λ, we have plotted energy density, anisotropic pressure and anisotropic parameter in case of collapse and expansion. The physical significance of the graphs has been explained in detail.

1 Introduction

The gravity is the fundamental tool which is so much related to our daily life problems, it still remains the most difficult to understand and interpret interaction from all the others. The gravitational force is the most form of without any sophisticated and deep knowledge, and the first one to be tested experimentally due to the nature and apparatus needed [1].
Although a lot of work has been done [2]-[4] to explore the dynamical aspects of stars without finding the exact solutions in f(R,T) gravity. The f(R,T) theory of gravity has been introduced by Harko et al. [5], for the formulation of this theory they have modified the Lagrangian of general relativity as a general function of R (Ricci scalar) and T (trace of the stress-energy tensor). They have formulated the equations of motion by using the metric approach instead of Platini approach. It has been investigated that the importance of T in the theory may be prominently observed by the exotic form of matter or phenomenological aspects of quantum gravity. The f(R,T) is an explicit generalization of f(R) theory, in which many cosmological and astrophysical results have discussed so far [7]. But, there is a still room to study some cosmological and astronomical processes in f(R,T) theory which have not yet been studied. In our present work the f(R,T) model can be selected in the following form

f(R,T)=f1(R)+f2(T).

(1)

Here, we take f1(R)=R and f2(T)=2λT
where R is Ricci scalar, λ is some positive constant and T is the trace of the stress-energy tensor as already mentioned.

Recently, Zubair et al. [6] analyzed the dynamical stability of cylindrically symmetric collapsing object with locally anisotropic fluid in f(R,T) theory. Alves et al. [8] investigated the existence of spacetime fluctuation in f(R,T) and f(R,Tϕ) theories of gravity. The study of collapse and dynamics of collisionless self-gravitating systems has been carried out by the coupled collisions using the Boltzmann and Poisson equations in f(R,T) gravity [9, 10]. Chakraborty [11] have proved that unknown generalized function f(R,T) can be evaluated in the closed form if this theory obeys the conservation of stress-energy tensor. Sharif and Zubair [12] derived the energy conditions in f(R,T) gravity which corresponds to the results of f(R) gravity. Houndjo et al. [13] investigated some Little Rip model in f(R,T) gravity using the standard reconstruction approach. Also, they remarked that the second law of thermodynamics remains valid for Little Rip model if the temperature inside the horizon is same that of apparent horizon.

Oppenheimer and Snyder [14] investigated the process of collapse in 1939, they observed the contraction of inhomogeneous spherically symmetric dust cloud. This study involves the exterior and interior regions as Schwarzschild and Friedman like solutions respectively.
It has been investigated in [15] that when massive stars collapse by the force of its own gravity, the final fate of such gravitational collapse is a white dwarf, neutron star, a black hole or a naked singularity. Misner and Sharp [16, 17] studied perfect fluid spherically symmetric collapse and also some authors [18]-[32] have discussed the phenomena of gravitational collapse using the dissipative and viscous fluid in general relativity. It has been shown that one can goes beyond the general relativity, the more one has chances of admitting an uncovered singularity [33].
In fact most modified gravity theories go out of their way to avoid making any changes to gravity near massive objects. The reason for this is that there is no evidence of any odd gravitational behavior near massive objects, so all of the modified gravity theories are designed to match the standard gravity at short non-galactic distances. There are also a large number of modified gravity theories that attempt to deal with dark energy, i.e. they explain the observations that lead us to think dark energy exists as a modified gravity in a similar fashion to how modified Newtonian dynamics is supposed to do away with the need for dark matter. As far as it has proposed a modified gravity theory that consistently accounts for all observational data, might work really well for galaxies as well as for cosmology. That being said, the current dark energy and dark matter model for the Universe is not completely without its own warts so modified gravity remains a possibility, albeit less likely on current evidence than the current model.

2 f(R,T) theory of gravity

Harko and his collaborators [5, 33] proposed a generalization of the f(R) theories, as f(R,T) gravity.
It depends on general function of R (Ricci scalar) and T the trace of tensor Tμν but in f(R) theories action depends on just Ricci scalar R. According to the authors, the dependence of theory on T (the trace of tensor Tμν) arises from quantum mechanical aspects which are usually neglected in f(R) or GR theories, for instance. The full action of f(R,T) gravity is in [5] as follows

S=∫d4x√−g(f(R,T)+Lm),

(2)

where Lm is matter Lagrangian and g is determinant of metric gab. Here, we choose Lm=ρ, and above action yields

Gμν=1fR[(fT+1)T(m)μν−ρgμνfT+f−RfR2gμν+(∇μ∇ν−gμν□)fR],

(3)

where T(m)μν is matter stress-energy tensor. The spacetime in this case has the following form

ds2=W2(t,r)dt2−X2(t,r)dr2−Y2(t,r)(dθ2+sin2θdϕ2).

(4)

The stress-energy tensor for anisotropic source is

Tuv(m)=(ρ+P⊥)VuVv−P⊥guv+(Pr−P⊥)XuXv.

(5)

Here ρ, Vu, Xu, Pr,P⊥ are energy density of matter, co-moving four-velocity of the source fluid, radial four vector, radial and tangential pressures, respectively. Also, for the given line element the quantities appearing in Tuv(m), must satisfy

Vu=W−1δu0,VuVu=1,Xu=X−1δu1,XuXu=−1.

(6)

The volumetric rate of expansion Θ is

Θ=1WXY(˙XY+2X˙Y),

(7)

where ⋅=∂t and ′=∂r.
The dimensionless anisotropy is

Δa=Pr−P⊥Pr.

(8)

For the given line element and stress-energy tensor the set of field equation is

3 Generating Solutions

For various possible values of α, we discuss the nature of solutions:

3.1 Collapse solution α=−52

When the value of expansion scalar is negative, we have gravitational collapse, so Eq.(19), implies that Θ<0, if α<−2, for the convenience we take α=−52 and the condition Y′=Y2α, provides Y′=Y−5, which further integrated to

Ytrap=(6r+h1(t))16.

(25)

Here h1(t) is integration function.
Using Eqs.(18), (25) in Eqs.(15)-(17), with some tedious algebra, we obtain the following explicit form of matter variables

Figure 1: This graph shows the variation of ρ with respect to r for the various values of λ and h1(t) = 1Figure 2: This graph shows the variation of Pr with respect to r for the various values of λ and h1(t) = 1Figure 3: This graph shows the variation of P⊥ with respect to r for the various values of λ and h1(t) = 1Figure 4: This graph shows the variation of △a with respect to r for the various values of λ and h1(t) = 1

In the above expressions h1(t) is an arbitrary function of time profile t, by taking h1(t)=1, we have analyzed the results. By choosing α=−52, we get Θ<0 and energy density remains positive and decreasing function of r the graphical behavior of ρ with various values of λ is shown in Fig.1. The radial pressure increases firstly and then decreases continuously with respect to radius at different values of λ as shown in Fig.2, but transverse pressure is decreasing with respect to radius as shown in Fig.3. In account of this we can say pressure is minimum when the value of λ is minimum as shown in 2 and 3. The maximum value of anisotropy occurs near the center of the sphere, so the anisotropic parameter Δa attains maximum value near the center and its values decreases when r increases with various values of λ it is given in Fig.4.

3.2 Expansion with α=32

When expansion scalar attains positive values, we have expanding solution, so Eq.(19), implies that Θ>0, if α>−2, for the convenience we take α=32 and assume that

Y=(r2+r20)−1+h2(t),

(31)

where h2(t) is integration function and r0>0. For simplicity, we take F(t,r)=1+h2(t)(r2+r20) and Y=F(r2+r20), then Eqs.(18), (31) with Eqs.(15)-(17), with simplification, we obtain the following explicit form of matter variables

Figure 5: This graph shows the variation of ρ with respect to r for the various values of λ and h2(t) = 1Figure 6: This graph shows the variation of Pr with respect to r for the various values of λ and h2(t) = 1Figure 7: This graph shows the variation of P⊥ with respect to r for the various values of λ and h2(t) = 1Figure 8: This graph shows the variation of △a with respect to r for the various values of λ and h2(t) = 1

.
In this case, when α=32 we get Θ>0 and energy density remains finitely positive for selected arbitrary time profile t, by fixing the values of α, the variation is given to λ under such restrictions the density ρ can be observed decreasing function shown in Fig.5. The radial pressure decreases continuously with respect to radius at different values of λ, but transverse pressure is increasing with respect to radius as shown in Figs.6 and 7 respectively. In account of this we can say pressure has reverse effects as observed in previous case. It is observed that the anisotropy is increasing function of r but with respect to the value of λ anisotropy is going to be maximum as the value of λ is being increased as shown in Fig.8.

4 Conclusion

Motivated by f(R,T) theory of gravity formulated by Harko et al. [5], a lot of work related to cosmology and stability of dynamics of collapsing stellar system has been done in recent years [34]-[49]. This theory has a vide range of cosmological and astrophysical applications in the era of modern physics. According to the available observation data, our universe is in the phase of accelerating expansion, to explain the physical significance of this phenomena a number of modifications to GR have been proposed. The f(R,T) theory of gravity is one has been the center of attention of researchers in current era this type of theories seems to provide a capacity of working successfully for dark matter. The conformal relation of f(R,T) to GR with a self-interacting scalar field has been discussed by Zubair et. al.[50].
Here, we have developed the generating solutions to collapse and expansion of fluid sphere in f(R,T) theory of gravity. For the interior matter distribution, the collapse solution yields a unique trapped surface. It has been investigated that during gravitational contraction, the phase transition would occur for the massive stellar system, for instance most of condensed matter configuration transit to a π-meson condensed state. The gravitational collapse is a highly dissipative process and a lot amount of heat energy is released during gravitational collapse Herrera at. el. [51]. In order to model the inhomogeneous cosmological solutions Collins [52] have explored the non-static expanding solutions. Also, the inclusion of anisotropic stress in the fluid source are much important, the influence of non-zero anisotropic parameter δa on the late late time evaluation of universe with non-homogeneous background has been explored by Barrow et. al. [53]. Due to the valid selection of R(t,r) and α, one can obtain anisotropy interconnection which collapses or expands by Glass [54]. In this paper, we extend the work of Glass [54] to the f(R,T) theory of gravity.

In this paper interior solution for anisotropic fluids have been discussed in detail, which are being used in modeling of anisotropic stars in the context of modified theory of gravity f(R,T). Using the auxiliary form of of the metric functions, we have determined the trapping conditions for fluid sphere in f(R,T) gravity. The resulting solutions have been classified as collapsing and expanding depending on the the nature of scalar expansion. The matter density, radial and transverse pressures, anisotropic parameter and mass function have been calculated in the context of f(R,T) theory of gravity. For the collapse solution when α=−52) the density decreases as the value of λ increase as shown in Fig. 1. The radial pressure Pr and matter density ρ have maximum values at center and it decreases from center to the surface of star. It has been observed that the anisotropy will be directed outward when Pr>P⊥ this gives that △a>0 as observed graphically in Figs.4 and 8. It is found in fig 4 that anisotropy decreases with the increase in radius.
Further, the expansion of gravitating source would occur when α=32, and Θ the expansion scalar is positive. In this case matter density decreases as shown in Fig.5. The radial/transverse pressures and anisotropic parameter with various values of λ have reverse behavior as compared to the case of gravitational collapse.