Anisotropic geodesic fluid in non-comoving spherical coordinates

We start with a recently introduced spherically symmetric geodesic fluid model (arXiv: 1601.07030) whose energy-momentum tensor in the comoving frame is dust-like with nontrivial energy flux. In the non-comoving energy frame (vanishing energy flux) the same EMT contains besides dust only radial pressure. We present Einstein's equations together with the matter equations in static spherically symmetric coordinates. These equations are self-contained (four equations for four unknowns). We solve them analytically except for a resulting nonlinear ordinary differential equation (ODE) for the gravitational potential. This ODE can be rewritten as a Lienard differential equation which, however, may be transformed into a rational Abel differential equation of the first kind. Finally we list some open mathematical problems and outline possible physical applications (galactic halos, dark energy stars) and related open problems.


Introduction
Analytic solutions of the coupled Einstein-matter equations for the stationary anisotropic and spherically symmetric case, without supplying any external input, are rather rare. The only example we know, a conformal flat generalization of the de Sitter space-time, has been published very recently ( [1], section 6). But normally one has to provide some external input. The generic case has been discussed recently by Herrera, Ospino and Di Prisco [2]. The authors of [2] provide only three Einstein equations for five unknowns (energy density, two pressures, two metric functions). So the knowledge of two solution generating functions is required.
In the present paper we describe another model for which we derive an analytic solution except for one remaining ordinary differential equation.
We start with a recently introduced irrotational geodesic fluid model whose energy-momentum tensor (EMT) in the frame comoving with the fluid is dust-like with nontrivial energy flux [3]. Then we pass over to the noncomoving energy frame (vanishing energy flux [4]). Here the same EMT contains besides dust only radial pressure. We consider the resulting Einsteins field equations together with the matter equations in static spherically symmetric coordinates. These equations are self-contained (four equations for four unknowns). We solve them analytically except for a resulting nonlinear ordinary differential equation (ODE) for the gravitational potential.
This ODE turns out to be the general relativistic generalization of a corresponding ODE derived in [5] for the nonrelativistic darkon fluid model. It has been used in [5] as a model for galactic halos.
The paper is organized as follows: We define our model in section 2.
In section 3 we introduce non-comoving coordinates and present the corresponding Einstein equations and the matter equations. In the course of integration of these equations in section 4 we derive a nonlinear ODE for the gravitational potential. We reformulate this ODE in section 5 as a Lienard differential equation and transform it into a rational Abel differential equation of the first kind. In section 6 we list some open mathematical problems.

Fluid model
Our model is defined by a self-gravitating, irrotational, pressure-less and stress free geodesic fluid whose EMT in the frame comoving with the fluid is dust-like with a nontrivial energy flux.
Therefore our model will be described by the following covariant set of equations (Greek indices run from 0 to 3 and we use the usual summation convention) • Einsteins equations (κ = 8πG, c = 1) with a EMT T µν , decomposed w.r.t. the unit and time-like fluid velocity vector u µ where ρ is the total energy density (comprising baryonic and the so called dark sector contributions) and q µ is the energy flux vector (u µ q µ = 0) in the comoving frame.
• Constraints for u µ Geodesic flow: Irrotational flow: The covariant derivative ▽ µ is given in terms of a torsion-free connection (Christoffel symbols).
• Covariant conservation of the EMT which is a consequence of (1) (Bianchi identities).
3 Einsteins equations and the matter equations in non-comoving, static and spherically symmetric coordinates Our fluid will be assumed to move with radial velocity v relative to the energy frame (EF). Such a choice of relative motion may be related to the observed motion of e.g. a galaxy relative to the microwave background [6].
For static spherically symmetric coordinates in the EF we use Schwarzschild (canonical) coordinates [7] Then u µ and q µ are given by (v = v(r), q = q(r)) with β = (1 − v 2 ) −1/2 . The time-like and space-like unit vectors n µ and s µ are defined by Sometimes it is convenient to use instead of λ(r) the mass function M(r) related to each other by The EMT (2) reads in the energy frame (decomposition of T µν w.r.t. n µ and where the energy density ρ * and the radial pressure p * r in the EF are related to the corresponding kinematic quantities in the comoving frame by In (11) we have used the relation which follows from the requirement of the vanishing energy flux in the EF [8].
With the metric (6) and the EMT (10) we obtain for the Einstein equa- Here and in what follows a prime denotes differentiation of a function w.r.t. to its argument.

Integration of the Einstein and matter equations
To integrate the set of independent equations (14)-(17) we will proceed in three steps: • Equ. (17) can easily be integrated • Insertion of ρ * = −p * r /v 2 from (11) and (18) into the TOV-equation (16) leads to which again can easily be integrated where α is an integration constant.
The result (24) turns out to be equal to the corresponding equation obtained for the stationary solution of the nonrelativistic darkon fluid model in [5]. In this limit we get from the Poisson equation and (24) the following relation between the energy density ρ and the potential φ Then positivity of ρ requires Note that a positive energy density yields a negative radial pressure according to (11).
5 Reformulation of the ODE (23) as a Lienard differential equation or as an Abel differential equation 1 With the transformation the ODE (23) becomes the autonomous ODE The further transformation leads to the mixed Lienard differential equation To transform (30) into an Abel differential equation we proceed as usual [10]: we obtain from (30) with By the further transformation we obtain from (33),(34) a rational Abel differential equation with Unfortunately (36) does not belong to the known integrable cases of rational Abel differential equations. But, as shown in [11](see also [12] and [13]), all integrable rational Abel differential equations consist of classes whose members are related to each other by the equivalence transformation where F, P and Q are arbitrary functions of z satisfying F ′ = 0 and P = 0.
A computer algebra routine has been presented in [11] which allow us to decide whether a given Abel differential equation belongs to one of the known integrable classes.

Open problems
From the results of section 5 follow immediately the following open mathematical problems: • For which values of γ does the Lienard equation (30) has positive so- • Check by means of the computer program presented in [11] whether the Abel equation (36) belongs to one of the known integrable classes.
• If the answer is no, elaborate numerical solutions for equation (36).
• Stability of the stationary solutions.

Physics
Analytic or numerical results for the gravitational potential φ(r) from solutions of either the nonlinear ODE (23) or of any of its equivalent forms given in section 5 will be suitable for the description of either galactic halos or of dark energy stars.

Galactic halos
A star in circular motion in a gravitational potential φ(r) possesses the tangential velocity v tg (r) given by the relationship which holds also in the general relativistic case (see [14]). Keeping in mind that in our model φ is sourced not only by stellar matter but also by the so called dark sector contributions, we may use a solution for φ in (39) for modeling of galactic rotation curves (RCs).
For the description of galactic halos we need gravitational potentials φ which vanish for r → ∞. But in the nonrelativistic case, described by (24), solutions vanish already for a finite but very large distance as shown in [5] by a theorem due to Taliaferro [15] as well as by numerical results. RCs at large radii for "dark matter dominated" galaxies (for a very recent review on the dark matter issues see [16]).
We do not expect any essential modifications of the weak field limit at large radii for the relativistic model presented in this paper.
On the other hand, at small scales, our nonrelativistic model seems to show numerically a vanishing RC already at a very small but nonzero radius [5]. It was not possible to give a definite answer to this point in [5] because the numerical solutions have shown a discontinuous behavior around the critical value b c for b = ϕ ′ (0).

Problem 2:
Behavior of the gravitational potential φ(r) for r → 0 in the weak field limit as well as for the case of strong fields.

Dark energy stars
As has been already stated, our model shows a negative radial pressure for a positive energy density. Therefore it is predestinated for the description of anisotropic dark energy (DE) stars. To do that we have to take a solution of our model as interior solution which have to be matched with the exterior Schwarzschild solution.
Recent treatments of anisotropic DE stars are to be found in [17], [18] and [19]. In all these cases analytic solutions are given by supplying some functions and constants as external input. In [17] proportionalities between energy densities and DE radial pressure as well as analytic expressions for the two metric functions φ and λ are assumed. In [18] two equations of state and an analytic expression for the DE energy density are given. An analytic expression for the mass function and a DE-equation of state for the radial pressure are provided in [19].
In our model no such external inputs are needed. But in order to proceed we have to succeed in finding (approximate) analytic expressions or numerical results for the gravitational potential.