Convection and cracking stability of spheres in general relativity

In the present paper we consider convection and cracking instabilities as well as their interplay. We develop a simple criterion to identify equations of state unstable to convection, and explore the influence of buoyancy on cracking (or overturning) for isotropic and anisotropic relativistic spheres. We show that a density profile ρ(r)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\rho (r)$$\end{document}, monotonous, decreasing and concave , i.e. ρ′<0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\rho ' < 0$$\end{document} and ρ′′<0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\rho '' < 0$$\end{document}, will be stable against convection, if the radial sound velocity monotonically decreases outward. We also studied the cracking instability scenarios and found that isotropic models can be unstable, when the reaction of the pressure gradient is neglected, i.e. δRp=0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\delta \mathcal {R}_p = 0$$\end{document}; but if it is considered, the instabilities may vanish and this result is valid, for both isotropic and anisotropic matter distributions.


Introduction
The stability of general relativistic self-gravitating matter distributions has been extensively studied and reported in the literature through several techniques for many years. It is complex multivariate problem which depends on the micro-physics -bulk/shear viscosity, crust on the surface, magnetic field and so on -of the material constituents and their description through a macroscopic equation of state that characterizes the configuration (see standard texts in Relativistic Astrophysics and Neutron Stars [1][2][3][4][5] and references therein).
The studies of stability distinguish two different approaches associated to the global and local scales, where instabilities affect the structure. On the global scale it is examined through the dynamical perturbation scheme which, in the case of spherical symmetry, can be translated into how radial pulsations induce possible disruptions of a stellar body. a e-mail: hector@ula.ve b e-mail: lnunez@uis.edu.co c e-mail: adrianacvr67@gmail.com On the other hand, local stability investigates the effects of convection and/or cracking of the material within the matter distribution. Pulsation is a global phenomenon characterized by the collective motion of the entire body while, convection/cracking occurs locally and is governed by the nearby values of the thermodynamical variables and their gradients (see an interesting review in [6]).
The dynamical instability approach studies the evolution of perturbations on the physical and geometrical variables throughout the matter distribution. In General Relativity, it arose from the seminal works of Chandrasekhar [7,8], Tooper [9,10] and Bardeen [11] and, a decade later, was formalized by Friedman and Schutz [12]. For an anisotropic fluid, this criterion bounds the adiabatic index as where ρ denotes the energy density, P the radial pressure, and v 2 the radial sound speed, respectively [13][14][15][16].
Within the global stability criteria we can also identify the Harrison-Zeldovich-Novikov condition, which implies that dM(ρ c )/dρ c ≥ 0, where M is the total mass of the configuration and ρ c the central density, of the distribution [16]. The stability of a spherical star to convection implies the buoyancy principle which leads to that pressure and energy density must, decrease outwards in any hydrostatic matter configuration [17][18][19]. Finally, the cracking instability approach determines the tidal acceleration profiles generated by perturbations of the energy density and the anisotropy of pressures identifying the changes sign of the total force distribution within the system [20][21][22][23]. This approach has been applied to an anisotropic fluid with barotropic equations of state [24] and, more recently extended to take into account the perturbation of the pressure gradient in both isotropic and anisotropic matter configurations [25,26].
In the present paper we consider convection and cracking instabilities as well as their interplay. We develop a simple criterion to identify equations of state unstable to convection, and also explore the influence of buoyancy on cracking (or overturning) of isotropic and anisotropic relativistic spheres. This paper is organized as follows: Sect. 2 describes the general equations of the theoretical framework of General Relativity. In Section 3 we formulate the concepts of adiabatic stability while cracking for self-gravitating anisotropic matter configurations is discussed in Sect. 4. Next we present in Sect. 5 the acceptability conditions which make any model physically reasonable. The models used and a discussion of our results for isotropic and anisotropic cases are presented in the Sects. 6 and 7. Finally, in Sect. 8 we wrap-up our final remarks.

The field equations
Let us consider a spherically symmetric space-time whose line element is given by with the regularity conditions at r = 0 : e 2ν(0) = const., e −2λ(0) = 1, At the surface of the sphere r = r b , the interior solution should match continuously the exterior Schwarzschild solution, which implies that: e 2ν(r b ) = e −2λ(r b ) = 1 − 2μ b , and the compactness at the surface defined as μ b = 2M/r b .
We shall consider a distribution of matter consisting of a non-Pascalian fluid with the energy-momentum tensor: with energy density ρ(r ), radial pressure P(r ) and tangential pressure P ⊥ (r ) of the fluid being determined by the Einstein field equations as where primes denote differentiation with respect to r . As it is well known, T μ r ;μ = 0 implies the hydrostatic equilibrium equation, or Tolman-Oppenheimer-Volkoff (TOV) equation which, if m(r ) = r 2 1 − e −2λ , can be written for this anisotropic fluid as and together with dm dr = 4πr 2 ρ , (8) constitute the stellar structure equations. In order to obtain the density and pressure profiles we have to provide two equations of state, P = P(ρ) and P ⊥ = P ⊥ (ρ), which for the present work will be assumed barotropic.

Adiabatic convection stability condition
The stability of a spherical star against convection can be easily understood. When a fluid element is displaced downward, if its density, ρ e , increases more rapidly than the surrounding density, ρ s , the element will sink downward and the star will be unstable. On the other hand, if the density of fluid element is less than its surroundings, it will float back and the star will be stable to convection.
Thus we can single out these three cases: 1. If ρ e > ρ s , gravity will tend to push the fluid element downward further and the system will be unstable. 2. If ρ e = ρ s , the system is considered neutral or metastable. 3. If ρ e < ρ s , a restoring force will act on the fluid element and then the system will be stable because it tends to its original state.
Following Bondi [17], let us denote the density, ρ(r p ), of an infinitesimal fluid element at its original position r p and displace this piece of material downward, thus: with δρ(r ) = ρ (r )(−δr ) and r = r p − δr , where r represents the current position of the fluid element, r p its original position and −δr the downward shift. Because ρ (r ) < 0, then δρ(r ) is a positive quantity, and the density of the compressed fluid element at the new displaced position will be greater that the density at its original position r p . On the other hand, expanding the density of the environment at the displaced position we get: The system will be stable against convection if the environment density is greater or equal than the density of the fluid element, we then have: thus ρ (r p ) ≤ ρ (r ) . Now, expanding ρ (r ) around r p we get which becomes the criterion of adiabatic stability against convection. Thus, density profiles with the second derivative less or equal than zero, ρ (r ) ≤ 0, will be stable against adiabatic convective motions. It is clear that parabolic density profiles, ρ = αr 2 + β, with α and β constants, will be stable against this type of convection, because buoyancy condition must also be fulfilled at the center r = 0 of the sphere: ρ < 0 ⇒ ρ c = 0, and ρ c < 0. These profiles have been implemented for the MIT Bag model through a linear equation of state, P = β(ρ −ρ s ), when densities become high enough for a phase transition to quark matter to occur [3,16].

Convection and cracking sources of instability
In this section we shall consider convective instabilities in the framework of cracking induced by perturbation in the density profile. Just for completeness we outline here the main concepts and equations concerning cracking for isotropic and anisotropic matter configurations, for further details, we refer interested readers to [25,26] and references therein.
As in the previous works [25,26], we assume that density fluctuations induce variations into all other physical variables, i.e. m(r ), P(r ), P ⊥ (r ) and their derivatives, generating a non-vanishing total radial force distribution (δR = 0) within the configuration. It is important to stress we are considering local perturbations of density, that can be properly described by any function of compact support, δρ = δρ(r ), defined in a closed interval r r b , where r b is the total radius of the configuration.
Accordingly, local density perturbations, ρ → ρ + δρ, generate fluctuations in mass, radial pressure, tangential pressure and radial pressure gradient, that can be represented up to linear terms in density fluctuation as: where are the radial and tangential sound speeds, respectively. Note that the present perturbation scenario contrasts the original one presented by Herrera and collaborators [20][21][22], where fluctuations in density and anisotropy were considered independent and simultaneous; and it is also different from a previous work [24] because there pressure gradient were not affected by the density perturbation.
Following [26], we formally expand the quantity R emerging from the TOV equation as: as where R 0 (ρ, P, P ⊥ , m, P ) = 0, because initially the configuration is in equilibrium. Next, by using (13)-(16) the above Eq. (19) can be reshaped as: where it is clear that density perturbations δρ(r ) are influencing: the distribution of reacting pressure forces R p , gravity forces R g and anisotropy forces R a . Depending on this effect, each perturbed distribution force can contribute in a different way to the change of sign of δR: each term can be written as with and Notice that if, as in [24], the perturbation δρ is constant and does not affect the pressure gradient, we have: δR p = 0, Thus, only anisotropic matter distribution can present cracking instabilities because δR g > 0 for all r and the possible change of sign for δR should emerge from δR a and the criterion against cracking is written as: more recently this equivalence between the restriction on pressures and velocities was demonstrated in [16] and included as part of the acceptability conditions that have to be considered when building physically reasonable compact object models. We shall discuss these constraints in the next section.

Physical acceptability conditions
More of the spherically symmetric perfect fluid "exact solutions" of Einstein field equations found in the literature are of little physical interest because, in addition to solving the structure equations (7) and (8) for a particular set of equations of state -P = P(ρ) and P ⊥ = P ⊥ (ρ)the physical and metric variables have to comply with several acceptability conditions which, over the years were recently compiled in [16] as: C1: Metric potentials, positive, finite and free from singularities; C2: Matching conditions at the surface of the star; C3: Decrease of interior redshift Z with the increase of r ; C4: Positive density and pressures; C5: Density and pressures having a maximum at the center and decreasing monotonically outwards, with P ⊥ ≥ P; C6: Energy conditions. Strong ( The adiabatic index stability criterion as stated in Eq. (1), which is a consequence of the dynamical stability criterion; C9: Stability against cracking as expressed by Eq. (26); C10: Harrison-Zeldovich-Novikov stability condition: Additionally Ivanov [16], demonstrated that these conditions are not independent and can be condensed in five main inequalities: 1. (m/r ) > 0 which fulfills C1, C2 and C3 conditions; 2. 0 ≥ P ⊥ ≥ P accomplishing: C4, C5, C6 (SEC), C7, C9; Clearly if we want the configuration to be stable against convection, we should add a sixth condition, i.e. the adiabatic convection stability condition ρ ≤ 0, to the above mentioned set. In the next section, we shall explore its influence on the stability of isotropic and anisotropic models.

Isotropic and anisotropic models
In this section we select seven exact solutions -describing isotropic or anisotropic fluid spheres -which comply with the Ivanov criteria. With this selection we study the effect of convective instability and the reaction of the pressure gradient to density perturbations.
Four of seven density profiles are among most physically reasonable isotropic solutions (Tolman VII [27], Buchdahl-1 [28], Mehra [29] and Kuchowicz [30]) reported in [31]. The fifth selection corresponds to a one-parameter family of a generalized Tolman IV solution obtained in [32] and allows us to exemplify the correlation of convection stability with a decreasing profile of (v 2 ) . Finally, we study two anisotropic solutions (Gokhroo and Mehra [33] and Sah and Chandra [34]) to illustrate the buoyancy effects in anisotropic matter configurations.

Isotropic solutions
In addition complying with the Ivanov criteria, isotropic model selection has significant physical interest in describing the interior of compact objects. The isotropic solutions shown in Table 1 are: Tolman VII and Mehra are the most frequent parabolic density profiles considered in models of stable neutron stars (see [35][36][37][38] and reference therein) while Buchdahl's solution [28,39] sets limits to the compactness of relativistic spheres. In [31] it is reported that the speed of sound for this solution does not decrease monotonically. We have shown that it can be attained for some particular values of the compactness μ = 2M/r b . Finally, models of charged spheres are frequently based of Kuchowicz solution [40] and in Mehra's solution [29] the density and the speed of sound vanish on the surface.
In reference [32], Lake proposed an algorithm based on the choice of a single monotone seed-function, F, to generate all regular static spherically symmetric perfect fluid solutions of Einstein's equations. Within this scheme we are going to consider the following family of models generated by where C is an arbitrary constant and N is a positive integer that produces an infinite family of analytic solutions. As we can see from the Eq. (27), different values of N recover wellknown solutions: N = 1 corresponds to Tolman IV solution [27]; N = 3 represents Heint IIa [41] solution; N = 4 and N = 5 are Durg IV and Durg V solutions, respectively [42]. The case N = 2 is considered in [43] studying the relationship between the central barionic density and the total mass of observed neutron stars. As can be easily guessed, the main difficulty of the method lies in how to calculate the two integrals that appearing in equation (4) of [32], but fortunately it is possible in the present case, as can seen in the appendix.
Thus, with the help of the following auxiliary function we obtain the density profile for any N , as where: Here, 2 F 1 (a, b; c; d) is the hypergeometric function, that for certain special arguments: (a, b; c; d) automatically evaluates to exact values with K a constant.

Anisotropic solutions
Local anisotropy (unequal stresses: P = P ⊥ ) in compact objects can be associated to different physical scenariossuch as phase transition, density inhomogeneity and electromagnetic field, just to mention a few of them -and has been considered extensively since the work of Bowers and Liang [44]. The unknown physics in the tangential equation of state, P ⊥ = P ⊥ (ρ) is partially compensated by using heuristic criteria: geometric, simplicity or any other assumption relating radial and tangential pressures (see [15,45] and references therein). Our selection of anisotropic solutions shown in the Table 2 are: Florides-Gokhroo-Mehra [33,46,47] and the Sah and Chandra [34].
The Florides-Gokhroo-Mehra profile was due originally to Florides [46], but also corresponds one of the different solutions considered by Stewart [47] and, more recently, by Gokhroo and Mehra [33]. The Florides-Stewart-Gokhroo-Mehra (FSGM) solution represents densities and pressures which, under particular circumstances [48], give rise to an equation of state similar to the Bethe-Börner-Sato newtonian equation of state for nuclear matter [1,2,49].

Modeling performed and discussion of some results
The first effect to be analyzed is the convective instability and its relation to the sign of the gradient of the speed of sound, , which are illustrated in Figs. 1, 2 and 3. We show this effect by plotting the normalized buoyancyρ = ρ /ρ c vŝ r = r/r b and the corresponding gradient (v 2 ) vsr . The first two figures display the convective stability for isotropic  [27] and Heint IIa [41], respectively) are unstable and (v 2 ) change sign. The N = 6 Model is stable to convection having (v 2 ) < 0 for allr . It is interesting to mention, that for this family of solutions, when N increases, we obtain more convective stable models  Tables 1 and 2 -are stable to convection because they have constant ρ < 0. The stability for the other isotropic models having rational density profiles is presented in Fig. 1, where the Buchdahl model becomes unstable because ρ /ρ c changes sign, while Kuchowicz is stable against convective perturbations due to ρ /ρ c > 0 for all r .
In Fig. 3 we illustrate the convective instability for anisotropic matter configurations. Again, models that are stable with the Ivanov criteria are revealed unstable for convection. This is the case of then Sah and Chandra model which is unstable against convective perturbation, but the Florides-Stewart-Gokhroo-Mehra solution is stable because it has a parabolic density profile.
In all the isotropic models analyzed we found an interesting correlation between the stability against convective perturbation and the sign of (v 2 ) . In particular, in the regions where matter configuration has ρ < 0 then we have v 2 as a monotonous decreasing function. And this is more evident in Fig. 2 where we have ploted the buoyancy and the gradient (v 2 )'. We found unstable configurations for N = 1, Tolman IV [27] and for N = 3 Heint IIa [41] with (v 2 ) changing sign within the configuration and for N = 6 a convective stable configuration having (v 2 ) < 0 for allr . This can be easily understood in those regions where P < 0 because Therefore, a distribution having a monotonically decreasing concave pressure profile, i.e. P < 0 and P < 0, will be sta-ble against convection, if the sound velocity monotonically decreases outward from the configuration. It is interesting to see that when it happens we have and this last condition is not only met by the model we have considered here, but also by any relativistic polytropic equation of P = Kρ γ or P = (γ −1)ρ and by other several more realistic numeric equation of state for ultradense matter (see [50], and references therein). The other effect to be discussed is the stability induced by the reaction of the pressure gradient to density perturbations. This is shown in Figs. 4 and 5, where we compare the perturbation on the total force distribution when the pressure gradient is perturbed and when it is not. In these figures we plot δR p = 0, and δR p = 0, respectively. The first case, δR p = 0, represents the recent cracking scheme of Gonzalez-Navarro-Nunez [25,26] while the second one corresponds to a variation of the previous work of Abreu-Hernandez-Nunez [24]. As we have stressed above, in the first case it is assumed a non-constant density perturbation δρ = δρ(r ), which leads to the factor 4πr 2 ρ ρ ≤ 0 in the third term in the perturbation to the gravitational force distribution δR g in Eq. (22), which may cause cracking instabilities.
As it is evident for the models considered, if the pressure gradient is not perturbed, i.e. δR p = 0, δR may change its sign and potential cracking instabilities may appear. On the other hand, if the gradient reacts to the perturbation, δR p = 0, we find that δR does not change sign and the matter configuration becomes stable against cracking. This tendency to make models cracking-stable if the the pressure Fig. 4 Stability against cracking, isotropic N -models: in these figures we show, the stabilizing effect of the reaction of pressure gradient to density perturbation. In these plots we present δR/δρ for two members of the Lake N -Family. As in the previous case this stabilizing effect is only present for N = 1, for N ≥ 2 no instability appears. Again, the more N increases the most stable the configuration is gradient reacts to the density, was previously reported incidentally in reference [26].
This induced stability can be understood if we shape the perturbation of the hydrostatic equation. Clearly, for the isotropic case, R a = δR a = 0 and when δR p = 0, the cracking instability emerges only from the effect of the perturbation on the gravitational force distribution, δR g and particularly from the third term, which is always negative. If δR p = 0 and P > 0, the reaction of the pressure gradient to the density perturbation can neutralize the effect of the negative sign of the above mentioned gravitational term.

Conclusions and final remarks
Stability is a key concept when considering self-gravitating stellar models, because only those in stable equilibrium are of astrophysical interest. A we have stated above, in addition to solving the structure equations (7) and (8) for a particular set of equations of state: P = P(ρ) and P ⊥ = P ⊥ (ρ), the emerging physical variables have to comply with the several acceptability conditions stated in Sect. 5.
The stability of a spherical star against convection has almost been forgotten in most of the stability analysis. It is a very simple criterion which implements the Archimedes principle in any hydrostatic matter configuration [17][18][19].  (plates a, c, e). It clear that these models are unstable to cracking when the reaction of the force distribution is not considered, i.e. δR p = 0. Notice that, contrary to what was shown in [26] Mehra models are stable when δR p is taken into account. This may happen because a possible mismatch in pressure and density for this solution in this reference This criterion has proven to be interesting because several reasonable models become unstable against convection and should be included in the acceptability criteria to guarantee physically interesting models of compact objects.
In this work we have shown that: 1. A density profiles with its second derivative with respect to the radial marker less or equal than zero, ρ < 0, will be stable against convective motions. This is a very simple criterion to identify potential convection instabilities within spherical matter configurations and it should be added to the above mentioned acceptability Ivanov criteria. 2. A decreasing concave pressure profile P(r ), i.e. P < 0 and P < 0, will be stable against convection if the radial sound velocity decreases monotonically outward in spherically matter configurations. Equivalently, if ∂ 2 P ∂ρ 2 > 0, then (v 2 ) < 0 implies stability against convective motions. It should be pointed out that to illustrate this effect we have implemented a new family of exact solution by using the method propose by in Ref. [32]. 3. From (19) we obtained the possible sources of cracking: the reaction of the pressure gradient δR p ; the perturbation to the gravitational force distribution δR g and the perturbation of the anisotropy δR a . Convection may cause cracking instability only when the perturbation of the pressure gradient is considered, i.e. δR p = 0, in this case cracking density perturbations affect the pressure gradient, depending on the values of (v 2 ) and v 2 ρ ρ . 4. Isotropic and anisotropic models considered can be unstable to cracking when the reaction of the pressure gradient is neglected, i.e. δR p = 0, but if taken into account, the instabilities may vanish. Thus, there is a stabilizing effect against cracking, when the perturbation gradient is affected by density perturbation, δR p = 0.
Local perturbed schemes are based on the reaction of the fluid variables to a density fluctuation that drives the system out of its equilibrium, i.e. we are exclusively considering perturbations under which the system is dynamically unstable. One way to achieve this is to assume for a barotropic fluid, that pressure gradients are not affected and we have shown that this occurs in several of the models considered. Convection contributes to the pressure gradient reaction but instead of developing or increasing cracking, it stabilizes the configuration.
Finally, it should be stressed that the results we have presented of possible instabilities for local perturbation schemes (convection and/or cracking) have to be considered as tendencies that could lead potential evolution of fluids within a relativistic matter configuration, but this should emerge from the full integration of the Einstein Equations. , then the density can be written as The resulting pressure is then Because of the boundary condition P b ≡ P(r b ) = 0, it follows that the constant K is: Where: