Spherically symmetric isothermal fluids in $f(R,T)$ gravity

We analyze the isothermal property in static fluid spheres within the framework of the modified $f(R, T)$ theory of gravitation. The equation of pressure isotropy of the standard Einstein theory is preserved however, the energy density and pressure are expressed in terms of both gravitational potentials. Invoking the isothermal prescription requires that the isotropy condition assumes the role of a consistency condition and an exact model generalizing that of general relativity is found. Moreover it is found that the Einstein model is unstable and acausal while the $f(R, T)$ counterpart is well behaved on account of the freedom available through an additional coupling constant. The case of a constant spatial gravitational potential is considered and the complete model is determined. This model is markedly different from its Einstein counterpart which is known to be isothermal. Dropping the restriction on the density and imposing a linear barotropic equation of state generates an exact solution and consequently a stellar distribution as the vanishing of the pressure is possible and a boundary hypersurface exists. Finally we comment on the case of relaxing the equation of state but demanding an inverse square fall-off of the density - this case proves intractable.


INTRODUCTION
Phenomenological theories of gravity have been on the increase in recent times. Such ideas purport to resolve the problems which are shortcomings of the standard Einstein's general theory of relativity. Specifically, the late time accelerated expansion of the universe is not a natural consequence of general relativity with the cosmological evolution equations suggesting a decelerating universe. This is inconsistent with experimental data conveyed by the supernovae Type 1a data [1], Baryon Acoustic Oscillations [2] and the WMAP survey involving the cosmic microwave background [3]. In order to resolve the difficulty, proposals of exotic matter fields have emerged. These include dark energy, dark matter, quintessence, phantom fields and the like. These latter fields do not as yet enjoy any experimental support even though their motivations may be sound.
An alternative approach is to reconsider the fundamental geometry prescripts. A modification of the action principle may have the potential to resolve the anomalies with the standard theory. For example, in f (R) theories [5] the action involves a polynomial in the Ricci scalar. It has been demonstrated that such an approach may indeed explain the accelerated expansion of the universe. It has been shown by Goswami et al [6] that the Buchdahl upperbound [7,8] for the mass-radius ratio of general relativity may be improved in f (R) theory with more matter per unit mass being admitted. The results also have implications for our understanding of the dark matter problem. The serious drawback of f (R) theory is the appearance of higher derivative terms which correspond to ghosts. It is usual in gravity theory to have at most second order equations of motion. Moreover it has been demonstrated [9] that f (R) theory is conformally related to the scalar-tensor field theory of Brans and Dicke.
The most general tensor theory of gravity admitting at most second order derivatives is the Lovelock theory [10,11]. The action consists of polynomials in the scalar invariants constructed from the Riemann tensor, Ricci tensor and the Ricci scalar. The drawback in this formalism is that the higher curvature terms are only active in dimensions higher than 4. That is Lovelock theory reduces to standard general relativity in dimensions 3 and 4 and makes a contribution to the dynamics from dimension 5 upwards. A special case of the Lovelock polynomial is the second order term known as the Gauss-Bonnet term that appears in the effective action of heterotic string theory [12]. The exterior field for a spherically symmetric star has been established by Boulware and Deser [13] for the neutral sphere and by Wiltshire [14] for the charged case in the mid 1980s. However, only recently were interior metrics found for perfect fluid astrophysical objects [15][16][17] in Einstein-Gauss-Bonnet gravity that could be matched to the Boulware-Deser [13] exterior metric.
If a scalar tensor action is sought then the most general such theory yielding second order equations of motion is due to Horndeski [18] and consists of the so called Fab Four components of the effective lagrangian. Several studies into its cosmological implications have been undertaken [19] and of late compact objects such as black holes and neutron stars were investigated by Silva et al [20]. Tensor multi-scalar theory of gravity has also recently come into vogue [21].
Harko et al [22] have proposed an action that is a function of the Ricci scalar R and the trace of the energy momentum tensor T which goes by the name f (R, T ) gravity. The equations of motion are indeed second order however the conservation of energy is sacrificed. This is ostensibly a drawback of the theory. However, it was argued by Rastall [23] that spacetime curvature could account for non-compliance with the Newtonian view of energy conservation [23][24][25]. Extensive investigations into the f (R, T ) paradigm have been conducted in recent times [26].
We examine the physically important case of perfect fluids displaying the isothermal property that is an inverse square law fall-off of density as well as a linear equation of state.
In such universes galaxies are considered as pointlike structures. By design such models can only describe cosmological fluids as no hypersurface of vanishing pressure indicating a boundary is present.
The paper is structured as follows: Firstly we review the essential ingredients of the f (R, T ) framework. We then derive the isothermal model in f (R, T ) theory and compare with the solution for Einstein gravity. In the next section we probe the consequences of a constant gravitational potential since it is known in Einstein gravity that a necessary and sufficient condition for isothermal behavior is a constant spatial gravitational potential.
Finally we impose a linear barotropic equation of state on our model but without any restriction on the density profile. Before we conclude with a discussion, we comment on the case of an inverse square fall-off of the density but without imposing an equation of state.

ELEMENTS OF f (R, T ) THEORY
The f (R, T ) gravity action is given by where f (R, T ) is an arbitrary function of the Ricci scalar R, and T is the trace of the energy momentum tensor T µν . The Lagrangian density L m for the matter field is defined as and its trace by T = g µν T µν . The Lagrangian density L m of matter has the form and is dependent only on the metric tensor components. Variation of the action (1) with respect to the metric g µν generates the field equations where f R (R, T ) = ∂f (R, T )/∂R and f T (R, T ) = ∂f (R, T )/∂T . ∇ µ denotes covariant differentiation and the box operator , is defined via The covariant divergence of Eq. (4) produces the equation which clearly shows that energy is not conserved in this system. With the help of Eq. (3) the quantity Θ µν is expressible as For the purposes of this investigation we consider a perfect fluid source with energymomentum tensor where p is the pressure and ρ the energy density of strange matter, with u µ u µ = 1 and u µ ∇ ν u µ = 0. If we take the matter Lagrangian density to be L m = −p, and the Eq. (6) we obtain the relationship where χ is a coupling constant constant. The field equations are now given by where χ can be positive or negative. Eq. (5) can now be written as and in the case of vanishing χ the law of energy conservation in Einstein gravity is recovered.

FIELD EQUATIONS
In coordinates (t, r, θ, φ) the most general spherically symmetric line element reads as where ν(r) and λ(r) are arbitrary functions of the radial coordinate r only. We consider a comoving fluid 4-velocity field u a = e −ν/2 δ a 0 and a perfect fluid source with energy momentum tensor given in Eq. (4). Additionally we use geometrized units such that the gravitational constant G and the speed of light c are taken as unity. Now Eq. (9) and Eq. (11) generate the field equations where the prime denotes the derivative with respect to the radial coordinate, r. Introducing the transformation e −λ = 1 − 2m(r)/r we obtain where the function m = m(r) represents the gravitational mass. An additional equation may be written from (10) dp dr that reduces to the energy conservation of general relativity when χ = 0. It is possible to rewrite Eqs. (12) and (13) in terms of energy density (ρ) and pressure (p) in the form while the equation of pressure isotropy G r r = G θ θ reduces to Observe that the equation of isotropy is the same for the ordinary Einstein's equations with a perfect fluid source. Therefore any of the well known solutions reported over the past century (for example see Delgaty and Lake [27]) will satisfy (18).

SOLUTION OF THE FIELD EQUATIONS WITH THE ISOTHERMAL PROPERTY
A perfect fluid is said to be isothermal if the density and pressure both obey the inverse square law fall-off and consequently display the equation of state p = γρ for some real number 0 < γ < 1 [28]. Accordingly let us insert where A and B are arbitrary parameters (at this stage) into equations (16) and (17). Observe that the field equations are essentially 3 in number and they contain four unknown functions.
Accordingly, specifying two of the quantities, namely the density and pressure, appears to be over-determining the system. This is true, however, we shall utilise the pressure isotropy equation as a consistency condition and determine the relationship between A and B for the isothermal property to hold. This is a similar route followed by Saslaw et al [28] in dealing with isothermal spheres in standard Einstein gravity.
Introducing (19) into (12) yields the differential equation which is written only in terms of the potential λ. With the help of the substitution e λ = β(r), equation (20) assumes the form where we have set w 1 = 8πB +χ(3B −A). Equation (21) is a Ricatti equation and is solvable in the form where C 1 is a constant of integration. Putting (19) into (13) simplifies it to the form where we have labelled w 2 = 8πA + χ(3A − B). (22) into (23) generates the solution

Now inserting
The isotropy equation (18) becomes and for consistency it is required that the coefficient of r and the constant term simultaneously vanish. This is achieved for which translates to A = χ − 4π 4 (χ 2 + 6πχ + 8π 2 ) and B = 7χ + 20π 4 (χ 2 + 6πχ + 8π 2 ) expressing A and B in terms of the coupling constant χ.
To ensure a subluminal sound speed requires 0 < γ = A B < 1 and this constrains the coupling constant to − 4π < χ < − 20π 7 (28) for the stability of the model. Moreover, this same interval guarantees that both density and pressure remain positive for all radial values. The mass of the infinite sphere as for some constant K. Observe that in the interval of validity above, the mass profile is a monotonically increasing function.
Setting χ = 0 above regains the standard Einstein field equations. Saslaw et al [28] obtained an isothermal cosmological model where the the geometric and dynamical variables are expressed in terms of the parameter α which comes from the equation of state p = αρ and which corresponds to A B in our formalism. We now make some observations about this model. Setting χ = 0 we get A = − 1 8π and B = 5 8π thus giving the specific value α = − 1 5 . Importantly this value is negative showing that the Saslaw model violates causality. In fact the constant value for the metric potential e λ = 1 + 4α (1+α) 2 evaluates to − 1 4 which is not feasible. Accordingly the Saslaw model is not realistic and violates the most basic requirements for physical plausibility.

CONSTANT GRAVITATIONAL POTENTIAL
It has been shown that a necessary and sufficient condition for isothermal behaviour, namely an inverse square fall off of the density and pressure, is a constant spatial gravitational potential λ. This is valid in Einstein theory and the more general Lovelock theory [29]. But what are the consequences of a constant potential in f (R, T ) theory? We now examine this question.
Setting Z = k for some constant k in the isotropy equation (18) gives for the remaining temporal potential. Note that k is now restricted through 0 < k < 2.
Introducing the substitutions a 1 = 2(χ + 4π) and a 2 = 8π + 3χ the density and pressure are given by respectively while the sound speed has the remarkably simple constant value dp dρ = 8π χ + 3 (33) and to ensure causal behaviour 0 < dp dρ < 1 it is demanded that χ obeys −4π < χ < − 8π 3 . The expressions governing the energy conditions have the form The active gravitational mass is calculated as where 2 F 1 is the familiar hypergeometric function.
In view of the complexity of the expressions for the dynamical quantities we conduct a qualitative study with the aid of graphical plots. The following parameter values have been appear to be well behaved for the f (R, T ) cases, the weak energy condition is violated for the Einstein case. Finally the plot of the mass profile (Fig 7) is reasonable. The mass increases more rapidly in the case of the f (R, T ) theory than compared to its Einstein counterpart.
In summary, the (f (R, T ) case displays more pleasing physical behavior than the Einstein case. in this approach is that the resulting nonlinear differential equation is difficult to integrate and to date no unique general solution is known. An alternative approach is to specify one of the four variables ν, λ, p or ρ and then to solve the system to reveal the remaining three.
Finally if the density or pressure equation is solvable for r in terms of ρ or p then a linear This equation of state is clearly not the most general one for p = αρ but represents a special case. For example see the seminal work of Tolman [30] wherein some Tolman models do indeed display equations of state. Note that in the Einstein case, specifying the density is tantamount to specifying the potential λ as the G t t = T t t equation only contains λ and ρ and is well known that the left-hand side may be expressed as an entire derivative. This is expressing ν in terms of λ. Substituting (38) into the isotropy equation (18) generates the differential equation governing the behaviour of e λ = β. Obtaining the general solution to (39) has proved elusive in view of the nonlinearity. The method of Lie group analysis was invoked however no symmetries could be detected immediately. However, on careful observation it is seen that in some cases (39) may be solved explicitly.
For the special case α = χ a 2 the isotropy equation becomes where we have redefined c 1 = (α + 1)a 1 and c 2 = −(α + 1)a 1 (αa 1 + a 1 − 4) + 4. Dividing throughout by the first term on the left we may rearrange equation (40) to the form with the help of partial fractions. The solution by quadratures may now be obtained implicitly as where K is a constant of integration. Equation (42) is essentially an algebraic equation in β(r). Clearly for judicious choices of the constants c 1 and c 2 , equation (42) may be solved explicitly to find the gravitational potential function β.
As an example, consider the choice c 1 = −2 and consequently c 2 = −8 follows. Now from c 1 α+1 = 2(χ + 4π) and the original assumption α = χ a 2 we solve simultaneously and obtain the pair We must discard the negative value of α since the causality criterion 0 < α < 1 will be violated. However, note that we are able to obtain the value α = 0.926 which indeed guarantees a subluminal sound speed. For this choice of c 2 equation (42) is solvable and the metric potential evaluates to which corresponds to the Vaidya-Tikekar [31] spheroidal geometry utilised to model superdense relativistic stars. In order to determine the remaining gravitational potential it is prudent to introduce the transformations x = 2K 2 r 2 , Z(x) = e −λ and e ν = y 2 (x) whence the equation of pressure isotropy assumes the form and is now a second order linear differential equation in y. Inserting Z = 1−x 1−2x into (45) generates the potential or in the canonical form The pressure and density are given by where we have made the further simplifications v 3 = 2v 1 + √ 2v 2 , K 1 = χ + 4π and K 2 = χ + 2π. Now we have a complete model with Vaidya-Tikekar [31] geometry and linear barotropic equation of state p = αρ. a defect in this model is that there exists an essential singularity at r = ± 1 2K . While the presence of the singularity is undesirable, it may not be a generic feature of this model. Suitable constants c 1 andc 2 may yet exist that support a well behaved cosmological model. Interestingly, the vanishing of the pressure for a finite r is possible allowing for the interpretation of this model as a bounded astrophysical distribution.

RELAXING THE EQUATION OF STATE
Finally we consider the case where the density displays an inverse square-law fall-off but we refrain from imposing an equation of state. That is the system of field equations is now completely determined and the resulting solution should be inspected for an equation of state. It turns out that equation (16) allows us to write ν in terms of λ and its derivative. When this form is substituted into the isotropy equation (18) the resulting differential equation proves intractable to solve. Note that this situation does not arise in the standard Einstein gravity since on setting χ = 0 for the Einstein case, (16) can be solved explicitly for λ in terms of r. This has been amply demonstrated by Dadhich et al [29,32] for the Einstein case and its generalization pure Lovelock theory.