Analytical study of anisotropic compact star models

A simple classification is given of the anisotropic relativistic star models, resembling the one of charged isotropic solutions. On the ground of this database, and taking into account the conditions for physically realistic star models, a method is proposed for generating all such solutions. It is based on the energy density and the radial pressure as seeding functions. Numerous relations between the realistic conditions are found and the need for a graphic proof is reduced just to one pair of inequalities. This general formalism is illustrated with an example of a class of solutions with linear equation of state and simple energy density. It is found that the solutions depend on three free constants and concrete examples are given. Some other popular models are studied with the same method.


Introduction
The study of relativistic stellar structure is now more than 100 years old. It began with the discovery in 1916 by Karl Schwarzschild of a universal vacuum exterior solution [1]. He also gave in the same year the first interior stellar solution [2], which should be matched to the exterior one. It has a constant energy density ρ. For a long time the star interior was considered to be made of a perfect fluid, which has equal radial ( p r ) and tangential ( p t ) pressures. This leads to the isotropic condition p r = p t , imposed on the Einstein equations. However, spherical symmetry demands only the equality of the two tangential pressures. This fact was noticed by Einstein and developed first by Lemaitre in 1933 [3]. He discussed a model sustained solely by p t and with constant ρ. His work remained unnoticed for a long time.
In 1972 Ruderman [4] argued for the first time that nuclear matter at very high densities of the order of 10 15 g/cm 3 may have anisotropic features and its interactions are relativistic. The pioneering work of Bowers and Liang [5] on building a e-mail: boykovi@gmail.com anisotropic models in 1974 gave start to a number of such solutions. Anisotropy may have a lot of sources [6]: a mixture of fluids of different types, presence of a superfluid, existence of a solid core, phase transitions, presence of magnetic field, viscosity, etc. Such models describe compact stellar objects like neutron stars, strange stars, quark stars, boson stars, gravastars, dark stars and others.
The Einstein equations describe the effect of matter upon the metric of spacetime. For static, spherically symmetric fluid solutions the metric may be written in canonical or isotropic coordinates and has two components ν and λ or ν and μ. The energy-momentum tensor is represented by its diagonal components, mentioned above: ρ, p r and p t . There are only three equations for these five characteristics, so that two of them may be chosen freely. They should satisfy, however, a lot of regularity, stability and energy conditions for a realistic model. The situation is analogous to the search of charged isotropic star models [7]. This is not surprising since charge can be looked upon as an effective anisotropy of the model [8]. Different choices of the two given functions have been made.
The simplest one is to propose ansatze for the two metric functions. Then the Einstein equations become expressions for the matter components. One of the first solutions was given in [9], where some of the Tolman isotropic solutions [10] were modified to become anisotropic. Another [11] relies on the metric of the well-known charged Krori and Barua solution [12]. Others were based on different isotropic solutions, found in the past [13,14]. A solution in isotropic coordinates also exists [15].
String theory has inspired embedding of branes like in the Randall-Sundrum model [16]. This rekindled the interest in stellar models embedded in 5-dimensional flat spacetime (embedding class one). They must satisfy the Karmarkar condition [17]. It can be written as a relation between λ and ν where means radial derivative. Consequently, one can choose one of these functions as generating the whole solution. It is interesting that the isotropic condition, which is an equality between the pressures, can be translated into a similar relation, giving different generating functions [18][19][20][21][22]. This is more easily done in isotropic coordinates, but canonical coordinates can be used too [20].
There are just two perfect fluid solutions of the Karmarkar condition -the interior Schwarzschild one [2], which has infinite speed of sound and a cosmological one. When the fluid is anisotropic, a plethora of realistic solutions have been found in the last two years. In some of them the generating metric component is a polynomial [23][24][25][26][27][28], in others it is a rational function [29][30][31]. There are also trigonometric [32,33], hyperbolic [34] and exponential generating metric components [35,36].
There is a group of hybrid solutions with given one metric and one matter function. Such is the model in higherdimensional spacetime with prescribed ν and the anisotropy factor Δ = p t − p r , which measures the diversion from isotropy [37]. An algorithm was given how to find any anisotropic solution provided these potentials are known [38]. It is based on linear differential equations, which are integrable. Different examples were given, but the regularity conditions were not studied. For isotropic solutions it becomes the algorithm of Lake [20]. It was used in [39], where a special function appears. Solutions were also found in isotropic coordinates [40,41] Conformally flat anisotropic spheres have vanishing Weyl tensor. This leads to a linear differential equation for ν , which can be integrated and even a relation between ν and λ is the outcome [42]. Several models were given with p r = 0 or prescribed λ. Recently, a conformally flat model with polytropic equation of state was discussed [43]. More solutions have been given by other authors [44,45].
Closely related are solutions which admit conformal motion. They depend on the conformal factor and a matter component, which can be ρ [46,47], or the mass m [48]. One can add here a model with given λ [49], since the expressions for ρ, m and λ are simply related. There is a model with a linear equation of state (LEOS) between p r and ρ [50], and another one with LEOS between the pressures [51]. Reference [42] has been generalised to non-static and electrically charged solutions [52].
The last group of known solutions are the models with two freely prescribed matter components. The interior Schwarzschild solution [2] is one of them, but it is isotropic with p r = p t and ρ = const. The first anisotropic solution was proposed by Lemaitre [3]. It has vanishing p r and constant ρ and was found independently by Florides [57]. It was studied further in [58]. A very important equation of hydrostatic equilibrium exists, involving only ρ, p r and p t in canonical coordinates. This is the TOV (Tolman, Oppen-heimer and Volkoff) equation [10,59] found initially for isotropic solutions. Its anisotropic version was used by Bowers and Liang [5] to find the first well-known star model, which has constant ρ, while Δ is given in a form suitable to solve the TOV equation. They asserted that anisotropic models may have arbitrarily large surface redshifts, which can explain the big redshifts of quasars. However, when the energy conditions are taken into account, realistic models have bounded redshifts [60]. The surface redshift, the mass and the radius of the star are characteristics that can be measured by astronomers. A model with constant ρ and nonvanishing p r has been given [61]. In [62] the method of [5] was further developed and a model with singular ρ was studied too. Bondi [63] searched for models with large redshifts, including a solution with constant ρ, and another, with constant Q = p r + 2 p t , and found that the latter shows more perspective.
Recently, it was shown that in isotropic coordinates the existence of an EOS gives an expression of λ in terms of μ, which may serve as a generating function. The case with p r = 0 was solved completely [64], giving an example of a new solution. The same was done for solutions with LEOS p r = αρ − β [65] and for the Chaplygin EOS [66].
A number of solutions with given ρ and p r , not linked by an EOS, are known [67][68][69][70]. Recently, solutions with prescribed λ and p r have been given full physical analysis [71][72][73][74]. At first sight, these are hybrid solutions, but since λ is closely related to ρ, we mention them here. Solutions, containing free Δ, exist in several combinations. Thus Δ, ρ solutions either simplify the TOV equation [75,76] or the Einstein equations, which acquire simple solutions [77,78], or even such in hypergeometric functions [79,80].
A subgroup of this last group are models, where p r satisfies some EOS. The simplest one is the so-called γlaw, p r = γρ, linear and without a free term. Anisotropy allows for much more solutions. This resembles the addition of charge to isotropic solutions [7,81]. There are models with prescribed ρ [82,83] or m [84]. For the usual compact stars γ ∈ [0, 1]. For dark energy stars it may be negative [85,86]. Quintessence stars have a second dark energy density ρ q imposed on normal matter [87][88][89][90]. Some earlier models of this type with p r = −ρ can be found in Refs. [91,92].
Another class of models have a linear EOS with a free term, the so-called MIT bag constant, which is suitable for more compact stars, p r = αρ − β. There are models with in addition given m, [93,94] or λ [95][96][97][98][99]. There are also models with quadratic EOS and λ [100][101][102]. The popular polytropic EOS from Newtonian gravity was shown to lead in the anisotropic Newtonian case to the well-known Lane-Emden equations [103] and to their relativistic generalisation in Einstein's gravity [104]. A simple solution for a simple ansatz for λ has been found [105]. Finally a solution with λ and the modified Van der Vaals EOS was presented recently [106].
The above classification of anisotropic star models does not pretend to be exhaustive, especially for the more exotic cases. There are additional recent references. We wanted to draw a global picture, showing where the general method for finding physically realistic solutions, proposed in the present paper, stands. In Sect. 2 the Einstein field equations are given, as well as the definitions of the main characteristics of a static anisotropic star. In Sect. 3 we summarise the conditions for a physically realistic model, amassed during the past decades. In Sect. 4 we argue that the model of type ρ, p r is the easiest one to implement these conditions. In Sect. 5 we derive different relations between the conditions, which reduce their number. In Sect. 6 an example is given-solutions with linear EOS and simple energy density. In this method the main object we study is the tangential pressure, which is done in Sect. 7. In Sect. 8 some other EOS are studied. Section 9 contains our discussion.

Field equations and definitions
The interior of static spherically symmetric stars is described by the canonical line element where λ and ν are dimensionless and depend only on the radial coordinate r . The Einstein equations read where ρ is the matter density, p r is the radial pressure, p t is the tangential one, means a radial derivative and Here G is the gravitational constant and c is the speed of light.
The gravitational mass in a sphere of radius r is given by Due to kc 2 , its dimension is length. Then Eq. (2) gives The compactness of the star u is defined by and is dimensionless. On the other side, the redshift Z depends on ν: The field equations do not contain ν, but its first and second derivative. One can express ν from Eqs. (2), (3), and (7) as The second derivative ν may be excluded by differentiation of Eq. (3) and combination with the other field equations. The result is where Δ = p t − p r is the anisotropic factor. Combining (10) and (11) one gets the well-known TOV (Tolman, Oppenheimer, Volkoff) equation [10,59] of hydrostatic equilibrium in a relativistic anisotropic star [5] p r = − ρc 2 + p r kr p r + 2m/r 2 The hydrostatic force on the left, F h , is balanced by the gravitational and the anisotropic forces, F g and F a , on the right. This equation is not independent from the field equations, but it is their consequence. It can replace one of them. It is also equivalent to the Bianchi identities T μ ν;μ = 0, which in the static spherically symmetric case have only one non-trivial component [5,75,76,104]. In CGS units G = 6.674 × 10 −8 cm 3 /g s 2 , c = 3 × 10 10 cm/s, k = 2.071 × 10 −48 s 2 /g cm, kc 2 = 1.864 × 10 −27 cm/g. From now on we set G = c = 1, passing to the usual relativistic units. Then k = 8π .
As a whole, we have three field equations for five unknown functions: λ, ν, ρ, p r and p t . We can choose freely two of them, but the model will be physically realistic if a number of regularity, matching and stability conditions are satisfied too.

Conditions for a physically realistic model
A comparatively reasonable set of conditions includes C1 The metric potentials are positive and should be finite and free from singularities in the star's interior and at the centre should satisfy e −λ(0) = 1 and e ν(0) = const.
C2 Matching conditions. At the surface of the star r = r s the interior solution should match continuously to the exterior Schwarzschild solution, This determines the metric at the surface In addition, the radial pressure there vanishes, p rs = 0. Neither the energy density nor the tangential pressure are obliged to do so. C3 The interior redshift Z , which, according to Eq. (9), depends only on ν should decrease with the increase of r . The surface redshift and compactness are related, due to Eq. (14): They should be less than the universal bounds, found when different energy conditions hold (see C6). In the isotropic case they are 2 and 8/9, respectively [107]. In the anisotropic case, when DEC holds, they are 5.211 and 0.974. When SEC holds, one has the bounds 3.842 and 0.957 [60]. They are greater than those in the isotropic case, but not arbitrary as asserted in [5]. C4 The density and the pressures should be non-negative inside the star. For ρ this coincides with the null energy condition (NEC). At the centre they should be finite ρ (0) = ρ 0 , p r (0) = p r 0 , p t (0) = p t0 . Moreover, p r 0 = p t0 .
C5 They should reach a maximum at the centre, so ρ (0) = p r (0) = p t (0) = 0 and should decrease monotonously outwards, ρ ≤ 0, p r ≤ 0, p t ≤ 0. The tangential pressure should remain bigger than the radial one, except at the centre, p t ≥ p r . An isotropic model is obtained when this inequality turns into an equality; this is called the isotropic condition.
C6 Energy conditions. The solution should satisfy the dominant energy condition (DEC) ρ ≥ p r , and ρ ≥ p t . When the pressures are positive, DEC is equivalent to the weak energy condition (WEC). It is desirable that even the strong energy condition (SEC) ρ ≥ p r + 2 p t is satisfied. Obviously, the latter encompasses DEC.
C7 Causality condition. It says that the radial and tangential speeds of sound should not surpass the speed of light. The speeds of sound are defined as v 2 r = d p r /dρ and v 2 t = d p t /dρ. Therefore this condition reads C8 The adiabatic index Γ as a criterion of stability. This index is the ratio of two specific heats and should be bigger than 4/3 for stability [6,108,109], C9 Stability against cracking. Cracking was introduced by Herrera [110] as a possibility of breaking of perturbed selfgravitating spheres. Abreu et al. [111] found a simple requirement for avoiding this to happen, namely that the region of stability be C10 The Harrison-Zeldovich-Novikov stability condition [112,113]. It implies that dM (ρ 0 ) /dρ 0 > 0.

General physically realistic solution
As shown in the introduction, one can choose the two free functions out of five in a number of ways. However, only the first three conditions in the previous section are imposed on the metric coefficients. The other concern the components of the energy-momentum tensor ρ, p r , p t . If we choose two of them we can satisfy many of the above conditions, at least partially, beforehand and determine the third one through the TOV equation. It does not contain the metric and replaces Eq. (4) from the original system of Einstein equations. The TOV equation gives a direct expression for p t , The mass m is obtained by integration of the density ρ (see Eq. (6)), so that the r.h.s. of the above equation involves only ρ and p r . If we try to express p r from TOV, we get a Riccati equation, where the coefficients R i can easily be extracted from Eq.
where we have used the opposite of Eq. (6), namely The Abel equation does not possess a general solution too. Thus, the natural step is to choose ρ and p r as free functions and express p t from Eq. (19). Then we find the mass from Eq. (6) and e λ from Eq. (7). It replaces the original Eq.
(2). The coefficient e ν is found up to a constant from Eq. (10), which replaces the original Eq. (3). Hence, instead of the system of Eqs. (2)-(4) we shall use the system of Eqs. (7), (10), and (19), where ρ and p r are chosen to satisfy the parts of conditions C4-C10, referring to them. One can even end with a single generating function ρ for the solution, by imposing an equation of state (EOS) on the radial pressure, Condition C10 concerns only ρ, while C8 concerns only ρ and p r , and they have to be satisfied by properly choosing these two functions.
Another group consists from the conditions upon the metric C1-C3. Equation (6) shows that m (0) = m r (0) = 0, because ρ is regular and finite at the centre. Equation (7) states that e λ(0) = 1. The mass, being an integral of a positive function, increases with the radius (which may be seen also from Eq. (22)) and at the surface Eq. (7) and C2 yield m (r s ) = M, the total mass of the star. Therefore m/r should increase from 0 to M/r s , but should not surpass 1 because a horizon appears then. Thus, e λ is an increasing function. The r.h.s. of Eq. (10) is positive, hence ν > 0 and e ν increases too. The arbitrary constant in it allows one to arrange the fulfilment of Eq. (14). Since it shows that e ν(r s ) < 1, the same is true for e ν(0) and the constant in C1 is smaller than unity. Finally, ν > 0 means that Z decreases going outwards (see Eq. (9)). The surface redshift and the compactness will satisfy the bounds in [60], as long as the DEC or the SEC holds. Thus, conditions C1-C3 are satisfied from general considerations and the pre-arranged requirement (m/r ) > 0.
The third group of requirements consists of the parts of conditions C4-C10, which refer to p t . One of them follows immediately from Eq. (19). At r = 0 we have p t0 = p r 0 because p r (0) = 0, ρ (0) and p r (0) are finite, m r (0) = 0 as shown above. Thus, p t0 is also finite and the solution satisfies C4 completely. There remain parts of C5, C6, C7 and C9 to be satisfied by p t .
The usual method used in the literature is to take two of the five essential characteristics of the model and choose simple expressions for them as polynomials or rational functions. Sometimes trigonometric and even special functions are used. They are supplied with enough parameters to try to satisfy the realistic conditions, starting with C1 and ending with C10. Already at this stage some of the expressions are so complicated that the authors pass to graphical description proofs, by means of figures. Even with one free parameter the graphics become 3D. This method is rather close to numerical simulations and not to an analytical study. Is it possible to reduce the number of graphic proofs? Do any relations exist between the numerous conditions, concerning just a few basic characteristics of the model, so that some of the conditions follow from the others in the general case? These questions will be dealt with in the next section.

Relations between the different conditions
Let us take the anti-cracking condition C9 and write it as Since the radial speed of sound from the causality condition is arranged to lie in the interval (0, 1], the l.h.s. is not positive. Combining inequalities (16) and (23) we get This can be written as Multiplying by ρ , which should be negative, results in the pair of inequalities Starting with Eq. (26) we may go backwards to Eq. (24). It can be replaced by the second part of Eq. (16), which gives a weaker pair of inequalities. Equation (24) may be replaced also by the weaker pair given by Eq. (23), which is nothing but Eq. (18). Thus we have proved that Eq. (26) is equivalent to C7 and C9 and if it holds, they hold too. We take Eq. (26) as the basic one that we have to satisfy. It shows that p t and p r decrease monotonously (for p r this is arranged beforehand) and that p t (0) = 0 as long as p r (0) = 0. The latter equality is also pre-arranged. All these enter C5, together with the pre-arranged conditions for ρ. C5 will be satisfied completely when we prove that p t ≥ p r . Equation (26) also shows that Δ increases with r .
In order to do this proof we first study the behaviour of inequalities under differentiation and integration. If g (r ) ≥ 0, g (r ) ≥ 0 does not necessarily follow, because the function g (r ) may oscillate up and down, having regions with negative derivative, but still remaining positive all the time. However, if g (r ) ≥ 0, the definite integral of it r 2 r 1 g (r ) dr ≥ 0 too, provided that r 2 > r 1 , because this is the area of the surface under the graphic of g (r ). The analogous conclusion also holds for a negative g (r ). There are two special radial points in a star model: r = 0 and r = r s . So we shall take integrals between 0 and r , or between r and r s . These lead to corollaries of the initial inequality and not to equivalence.
The easiest function to integrate is a derivative, g = h . This is one of the reasons to start from the end of the list of conditions, where derivatives prevail the inequalities, and work backwards, contrary to what is usually done.
Thus, taking the integral r 0 of Eq. (26) we obtain We have shown earlier that TOV leads to p t0 = p r 0 , hence, Eq. (27) becomes p r 0 ≥ p t ≥ p r . In this way C5 holds in total, provided that Eq. (26) is true. It remains to satisfy C6, that is, the two energy conditions. Let us deal first with DEC. The first part of C7 may be written as Integrating this inequality from r to r s and taking into account that p rs = 0, we find Thus, DEC follows for p r from the causality condition. In fact, we can pre-arrange both of them. Doing the same procedure for p t yields a similar result, Hence, if DEC holds at the surface, it also holds in the interior. Now, at the surface the expression for p t following from Eq. (19) becomes Then ρ s − p ts ≥ 0 turns into The r.h.s. is negative, while the denominator in the l.h.s. is positive. Consequently, a sufficient condition for DEC will be a positive numerator, It is a rather realistic condition. Let us discuss next SEC. A sufficient condition, which depends only on ρ and p r and therefore can be pre-arranged, reads Then SEC follows from the chain of inequalities 2 p t + p r ≤ 2 p t0 + p r = 2 p r 0 + p r ≤ ρ, (35) which are due to the fact that p t decreases and is equal to p r at the origin. Equation (34) is also a rather mild restriction. Summarising, the pair of inequalities (26) are the only conditions that p t must satisfy. Before discussing them we present a relation about C8. A sufficient condition for it to take place is a lower limit on the radial speed of sound which is simpler. Indeed, when SEC (C6) holds and we apply C5 we get ρ ≥ 3 p r . Put in C8 it gives the above inequality. Now let us take Eq. (27) and place in it the expression for p t from Eq. (19). After some manipulations we get −r p r ≤ 2F ≤ −r p r + 2 ( p r 0 − p r ) , Applying the L'Hospital rule to Eq. (6) yields Therefore, F behaves as r 2 when r → 0. The left and the right bounds should behave the same way, otherwise one of the inequalities will break near the origin. This means that p r ∼ p r 0 − c 1 r 2 , p r ∼ −2c 1 r there. Hence if we choose p r = p r 0 − c 1 r n as a simple expression for the radial pressure, only n = 2 may fulfil Eq. (27). If p r satisfies an EOS, we have p r = f (ρ) ρ ρ ∼ r and ρ ∼ ρ 0 −c 1 r 2 . If ρ = ρ 0 −ar n [56], only n = 2 has a chance. Here c 1 , a, n are some positive constants.

Solutions with linear EOS and simple energy density
Let us choose ρ as a simple binomial and p r satisfying a linear EOS with a bag constant, This kind of density appears already in the Tolman VII isotropic solution [10]. Anisotropic models with such a density can be found in Refs. [45,56,[67][68][69][70]82,95,99,100,102,105,106]. Models with such LEOS are discussed in [65,[93][94][95][96][97][98][99]. Obviously p r 0 = β (ρ 0 − ρ s ) > 0, p rs = 0. Thus C4 holds for them. ρ and p r have a maximum at the centre and decrease outwards, satisfying their part of C5. We also introduce the notation Now d p r /dρ = β, so the first causality condition in C7 holds for 0 < β ≤ 1. Then DEC holds too. The sufficient condition for SEC (34) becomes This inequality is true, if it is true for ρ = ρ s , which yields Next we can obtain a more refined estimate for the adiabatic index Γ and C8. Inequality (17) leads to It is true when β ≥ 1/3 since the r.h.s. is always negative, because β ≤ 1, while the l.h.s. is non-negative. Hence, C8 holds, in accord with the sufficient condition Eq. (36). When β < 1/3 we get which holds, provided it holds for ρ 0 . This leads to α < 3/ (4 − 3β). Then C8 is also fulfilled.
Then ρ, p r may be written as The mass m is given by Eq. (6), As a function of ρ 0 , the total mass M and its derivative become Thus C10 is satisfied and the specific density and radial pressure satisfy all realistic conditions.
Let us turn to the metric functions next. Putting Eq. (48) into Eq. (7) yields an expression for λ: This is a polynomial of the second degree in x [98]. We have pointed out that m/r should increase with r from 0 to M/r s . Equation (48) gives This will be positive when ρ 0 < 6ρ s , equivalent to α < 5/6. The function e ν is found by integrating Eq. (10) which is a rational function in y. The result is given by Eq. (3.4) from [98]. These functions satisfy the conditions on the metric C1-C3, as was shown before.

Equation (19) may be rewritten as
The inequalities (26) become Using the equations in the previous section we get The function in the middle should lie between two simple constant bounds. It is given by where all of them being positive. Next we find Equation (12) shows that F is closely connected to the gravitational force in the TOV equation, F g = −2F/r . We expect that F y decreases. Then it satisfies Eq. (55) when it does so for r = 0 and r = r s . At r = 0, y = 0 and we have 12γ where γ = αβ. At the surface r s , y = 1 and where we have replaced b with l for simplicity, The final result is four intricate inequalities for the three free constants α, β, b. Equation (41) shows that they hide inside the physical constants ρ 0 , ρ s and r s , from which through Eq. (48) one can find the other main star characteristics M, u s , Z s . In addition, α should satisfy Eqs. (43) and (51), while the surface compactness should satisfy Eq. (33).
A solution with the same functions ρ and p r was given in [98]. It has ρ 0 = 3.98 × 10 15 g/cm 3 , ρ s = 3.29 × 10 15 g/cm 3 so that, finally, is the interval for a realistic b. We can also show graphically that H ≡ F y /γρ 0 decreases for the interval (66), as required, and stays between 1 and 2, because of Eq. (55); see Fig. 1. The middle line is for b = 1, while the lower and the upper ones correspond to the limits in Eq. (66). This is the only time we use a graphic proof.
In [98] b = 1.074, which is in the above region. One obtains also u s = 0.299b = 0.321, a value well below 0.8 and the bounds in [60]. Therefore, such a solution satisfies C1-C10 and is physically realistic. The star is chosen very compact, with radius 3.8 km. The solar mass is Astronomers have observed even more compact, presumably quark stars, such as PSR B0943+10 with radius 2.6 km, described also by an analytical model [101], Table 1 and [102], Table 1. Another observed star mentioned in [102], named RX J1856.5-3754, has radius of 3.5 km, so the above value of 3.8 km is realistic.
Next, let us find another cluster of solutions by choosing ρ s /ρ 0 = 0.82, (α = 0.18), β = 1/3. Then γ = 0.06. Equations (61) and (62)  If we take a constant density solution, ρ = ρ 0 = ρ s , we shall have α = 0 and u s = 0.333b, which is close to the values discussed above. This is not surprising, since the density of the crust is taken high, more than 80% of the central density. However, this solution is unphysical, since the speed of sound in both directions is infinite. The solutions above satisfy C1-C10 and are physically realistic. There is an interesting relation, following from that. Equation (41), written in CGS units, yields when ρ 0 is measured in g/cm 3 , while r s is measured in km. It shows that there are models with different radius, but the same b. Let us choose ρ 0 = ρ 1 × 10 15 g/cm 3 where ρ 1 is a number between 1 and 10. The previous equation becomes Let us take the highest b found in the two examples, b = 1.152. Then the above equality holds for a very compact quark star if r s = 4 km and ρ 1 = 3.863. However, it also holds for a neutron star with r s = 7.861 km and ρ 1 = 1. In both cases the ρ 0 is of order 10 15 g/cm 3 . Equation (48) shows that u s stays constant with b, hence, the total mass M should increase when r s increases. One can call such solutions bisotopic.

Some other EOS
We have already pointed out that models with constant ρ are not physical, because the speed of sound becomes infinite. Conditions C7-C9 become singular and make no sense. The first interior solution was such a model [2] and others followed [3,5,57,58,[61][62][63]76]. They can serve as an approximation to some features of the stars. The same may be said for models with p r = 0 and any ρ [3,57,58,64]. Because of C4 we should have p t0 = 0, and because of C5 p t ≤ 0. Since, according to C4, the pressures are non-negative, p t = 0 throughout the bulk too. What remains is a dust model with only ρ = 0. Such models are unstable and collapse.
They follow from Eq. (39) when ρ s = 0 (α = 1). The energy density vanishes at the surface, so such stars are gaseous and have no crust. The tangential pressure is given again by Eq. (52). Now Then F ∼ y (1 − y) and vanishes at r = 0 (y = 0) and r = r s (y = 1). Furthermore, Thus that is, the tangential pressure changes sign and becomes negative. This is not allowed by C4 and the model is unphysical with such energy density. Our final example are models with the same ρ, but with polytropic EOS [95], where K is the polytropic constant and n is the polytropic index. Once again ρ s = 0 and α = 1. Then Eq. (39) leads to and again F ∼ y (1 − y). On the other side Then and vanishes at the centre and at the surface of the star. Thus Δ x changes sign somewhere. However, C5 necessitates Δ x to be positive so that Δ increases. Hence, this model is unphysical, too. One should change the simple expression for ρ. These examples show that it is not so easy to find physically realistic star models.

Discussion
With five unknown functions and three equations, one has 10 combinations of two free functions, claiming to give the general solution for anisotropic static stars. The panorama given in the Introduction shows that most of them have been applied to find at least one concrete solution. Instead of p t the anisotropic factor Δ has been used. Sometimes m is taken instead of ρ or p r . The general solution was seriously discussed in [38]. The main shortcoming is that the conditions for a realistic model are checked after the choice is made. The expressions for the different characteristics become very involved even for polynomial seeding functions and one has to turn to graphic descriptions, which serve as proofs. Solutions are usually supplied with lots of constants, to negotiate the set C1-C10. Just one constant is enough to turn a 2dimensional plot into a 3-dimensional one, whose projection as a 2-dimensional flow of lines is often used. With two and more constants only partial plots are possible.
We have argued that the combination of free functions ρ, p r seems to be the right choice to reduce the number of graphic proofs. There is a simple chain of relations between ρ, m, and λ, due to Eqs. (6) and (7), so that one may start with some of the last two and p r . The advantage of this solution generation method is that we can arrange beforehand that the part of the condition set, referring to ρ, and p r , is satisfied. This concerns the majority of conditions, namely C4-C10. Besides, C8 and C10 contain only ρ and p r .
The set of conditions C1-C10 was formed over the years. In older papers just the new solution is given and sometimes its matching to the exterior. Later the conditions for monotonous decrease and the energy conditions were added. Then followed the causality and stability conditions, the last being C10, which, however, dates back to 1965. These conditions come from phenomenology, where hydrostatics, nuclear theory and thermodynamics meet with general relativity. In theoretical astrophysics they stand as axioms that must be satisfied. One is tempted to ask whether all these conditions are independent. We have shown that this is not the case and reduced them to the couple of inequalities in Eq. (26). Only Eq. (26) needs a graphic proof in the concrete examples. To emphasise this point we have not given any other graphics, even for illustration purposes. This reduction is possible due to the following factors.
The characteristics of static models, unlike the dynamical ones, depend on one variable, the radius. This is true even for the stability criteria, based on time-dependent perturbations. Thus, there are only ordinary derivatives and all of them can be replaced by r -derivatives, e.g. d p r /dρ = p r /ρ .
Inequalities may be integrated, obtaining other true inequalities, which are corollaries, not equivalent to their parents.
Most of the features of the model have simple monotonous behaviour as functions of the radius. Thus, ρ, p r , p t , Z decrease, while λ, ν, m, Δ, u increase with r . Therefore, one can multiply the inequalities with their derivatives, changing the sign when necessary, without bothering about the region of application-it remains the same. When there is a bump (mainly in p t ), this signals that something is wrong with the causality and/or the anti-cracking conditions. Thus, many of the conditions follow from the others for any solution. When an EOS is specified, some of the conditions are fulfilled, provided only ρ s and ρ 0 are known and not the whole graphic of ρ.
To illustrate this general formalism we have given a solution with simple energy density and simple EOS, namely linear EOS with bag constant. We have found that the general solution depends on three free constants. They implicate the physical constants ρ 0 , ρ s , p r 0 and r s , which lead in turn to other main characteristics like M, u s , and Z s , some of which are measured in astronomy. We have found regions where the solution is realistic and not just a point in the constants' space. In the previous section we have shown why some popular past solutions break certain realistic conditions. This is true even for the γ -law and the polytropic EOS with the same simple density. Hence, other expressions for ρ should be studied. We hope that the application of the method, described in this paper, will lead to many other solutions in the future, obtained in an easier study of genuinely analytic nature.