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 explaining the late time accelerated expansion of the universe remains an open question despite the successes of general relativity in satisfying fundamental experiments such as solar system tests. The seminal review article by Debono and Smoot [1] considers these anomalies and examines the question 'why consider alternative theories'. The accelerated expansion of the unia e-mail: hansrajs@ukzn.ac.za verse has been confirmed by a number of programmes such as supernovae Type 1a data [2][3][4], Baryon Acoustic Oscillations [5] and the WMAP survey involving the cosmic microwave background [6]. 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 [7] 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. [8] that the Buchdahl upperbound [9,10] 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 [11] 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 [12,13]. 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 provided that the field equations are generated from the metric tensor only. If the Lagrangian includes coupling with a scalar field such as a dilatonic field then the dynamics and geometry in 4 dimensions are influenced by the higher curvature effects [14]. Lovelock theory reduces to standard general relativity in dimensions 3 and 4 if no scalar fields are involved in the action and makes an active 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 [15]. The exterior field for a spherically symmetric star has been established by Boulware and Deser [16] for the neutral sphere and by Wiltshire [17] for the charged case in the mid 1980s. However, only recently were interior metrics found for perfect fluid astrophysical objects [18][19][20] in Einstein-Gauss-Bonnet gravity that could be matched to the Boulware-Deser [16] 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 [21] and consists of the so called Fab Four components of the effective lagrangian. Several studies into its cosmological implications have been undertaken [22] and of late compact objects such as black holes and neutron stars were investigated by Silva et al. [23]. Tensor multi-scalar theory of gravity has also recently come into vogue [24].
Harko et al. [25] 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 [26] that spacetime curvature could account for non-compliance with the Newtonian view of energy conservation [26][27][28]. Extensive investigations into the f (R, T ) paradigm have been conducted in recent times [29]. The effect of the f (R, T ) modification on radiating stars was discussed by Yousaf et al. [30]. Moreover Yousaf and other collaborators [31][32][33][34][35][36][37] have considered compact structures within this framework but with anisotropic pressures.
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 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 energy-momentum 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 Following Harko et al we consider the simplest version of The field equations are now given by where χ can be positive or negative. Equation (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 Eqs. (9) and (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 [38]) 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 [39]. Accordingly let us insert where A and B are arbitrary parameters (at this stage) into Eqs. (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. [39] 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 ), Eq. (20) assumes the form where we have set w 1 = 8π B + χ(3B − A). Equation (21) is a Riccati 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). Now inserting (22) into (23) generates the solution 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 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π or χ > 4π for the stability of the model. A positive pressure exists for −4π < χ < − 20 7 π or χ > −2π while the positivity of the energy density is satisfied for −4π < χ < −2π or χ > 4π . Reconciling all of these boils down to χ > 4π for a stable causal configuration of perfect fluid with a linear barotropic equation of state and with the isothermal property. The mass of the infinite sphere behaves as for some constant K . Observe that in the interval of validity above, the mass profile is a monotonically increasing function.

The Einstein isothermal model
The Saslaw et al. isothermal solution in Einstein gravity is given by the geometric variables and dynamical quantities where H is an integration constant and α = p ρ is a proportionality constant obeying 0 < α < 1. The positivity of e λ requires α < −3 − 2 √ 2 or α > −3 + 2 √ 2. Moreover, a positive density demands the window −3 − 2 √ 2 < α < −3 + 2 √ 2 or α > 0 while the interval α < −3 − 2 √ 2 or −3 + 2 √ 2 < α < 0 guarantees a positive pressure. Finally to ensure a causal stable fluid the constraint 0 < α < 1 must be enforced. Now routine checks show that there exists no values of α such that all these constraints may be simultaneously satisfied in some interval on the real line. Accordingly the Saslaw model violates one or more physical requirements and is therefore not realistic. The modified gravity model we have presented in the previous section does not suffer this defect provided that the coupling constant satisfies χ > 4π .

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 [40]. 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 ( 3 4 ) and to ensure causal behaviour 0 < dp dρ < 1 it is demanded that χ obeys −4π < χ < − 8π 3 or approximately −12.566 < χ < −8.378. 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 used to generate the plots c 1 = 1, c 2 = 2 and k = 1.5. These choices have been made following empirical fine tuning. Additionally we consider three different values for the coupling χ namely χ = −12 (thick curve), χ = −10 (thin curve) and χ = 0 (dashed curve). These values of χ are so chosen to coincide with the interval of validity ensuring a subluminal sound speed determined by Eq. (34). We provide two independent choices of χ to indicate that the graphical profile is generic in the causally valid region. The third plot χ = 0 corresponds to the model in the standard Einstein framework.   (Figs. 4, 5, 6) 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 In summary, the f (R, T ) case displays more pleasing physical behavior than the Einstein case.

Equation of state
In the Einstein framework, imposing the equation of state p = αρ determines a relationship between the metric poten- Fig. 7 Mass versus radial value r tials ν and λ. It is possible to isolate ν and substitute this into the equation of pressure isotropy -also an equation connecting ν and λ. The caveat 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 barotropic equation can easily be determined albeit that the expressions are lengthy. For example, if the density profile is prescribed in such a way that the resulting equation can be arranged as a polynomial equation up to quartic order in r , then the equation can be solved for r in terms of ρ. Substituting r in the expression for p gives the equation of state. 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 [41] 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 not the case in the f (R, T ) scenario where both ν and λ appear in the same equation with ρ. For this reason the incompressible fluid (constant density) solution is still unknown in f (R, T ) gravity. However, there is some extra latitude present through the constant χ and an equation of state may be determined as shall be demonstrated below.
Imposing the equation of state p = αρ results in the relationship expressing ν in terms of λ. Substituting (39) into the isotropy equation (18) generates the differential equation governing the behaviour of e λ = β. Obtaining the general solution to (40) 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 (40) 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 Eq. (41) 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 (43) is essentially an algebraic equation in β(r ). Clearly for judicious choices of the constants c 1 and c 2 , Eq. (43) 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 or given approximately numerically as {(χ , α)} = {(−13. 0854, 0.926495), (−6.51411, −1.16523)}. 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 Eq. (43) is solvable and the metric potential evaluates to which corresponds to the Vaidya-Tikekar [42] 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 (46) generates the potential or in the canonical form where we have put v 1 = √ 1 − 2K 2 r 2 and v 2 = √ 1 − 4K 2 r 2 . 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 [42] 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 Eq. (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. [40,43] for the Einstein case and its generalization pure Lovelock theory.

Conclusion
We have analysed the isothermal property in the framework of f (R, T ) theory. Demanding an inverse square fall-off of the density and the equation of state p = αρ yielded an exact model where the proportionality constant α is expressed in terms of the coupling constant χ . For stability and to prevent super-luminal behavior of the fluid the value of χ was constrained to a certain negative window. On setting χ = 0 we regain the Saslaw et al model for standard Einstein gravity and we discover that it is not physically reasonable. In contrast, the f (R, T ) model displayed the necessary features corresponding to expectations, namely a positive definite density and pressure and a sound speed obeying causality. While it is known that a constant spatial potential guarantees isothermal behaviour in the Einstein case and its generalization Lovelock gravity, such a prescription behaves completely differently in the ( f (R, T ) gravity framework. Dropping the inverse square law requirement and requiring an equation of state, the f (R, T ) model is indeed solvable in at least one special case. We have given a prescription to determine other models which satisfy the field equations and the equation of state. The case of an inverse square fall-off of the density without an equation of state did not yield an exact solution.