CC BY
10
УДК 530.12:531.51
The Demonstration of Gravitational Phase Transition inside of Fluid Static Ball
Alexandre M. Baranov*
Institute of Engineering Physics & Radio-Electronics, Siberian Federal University, Svobodny, 79, Krasnoyarsk, 660041, Siberian State Technological University, Mira, 82, Krasnoyarsk, 660049, Russia
Received 10.11.2010, received in revised form 10.12.2010, accepted 20.02.2011 It is shown that at the center of gravitating fluid ball there is the phase transition of second kind of gravitational field as a change of field's algebraic type.
Keywords: a gravitational phase transition, Einstein's equations, General Relativity and Gravitation, algebraic classification of Petrov, the theory of catastrophes.
The question about solutions of gravitational equations and researches of symmetric static models in General Relativity and Gravitation (GRG) today is also actual as yesterday. Petrov's algebraic classification of gravitational fields [1] is known as one from high-power tools of gravitational field investigations. Therefore a consideration of this problem from more general positions allows to analyse the fluid ball behaviour from the point of view of the Petrov algebraic classification.
In this paper is considered a spherical model of a gravitating fluid ball from more general base than in the papers [2,3]. Here we suppose that our model has spherical symmetry and is static, i.e. all metric functions don't depend from the time variable.
Metric, Tetrads and Einstein equations
A metric interval is chosen here in Bondi's coordinates and written down as
ds2 = F(r)dt2 + 2L(r)dtdr - r2(dd2 + sind2d<p2), (1)
where functions F = F(r) and L = L(r) are the metric functions, r is a radial variable; t is a time variable; d and f are the angle coordinates. The velocity of light and Newton's gravitational constant GN are chosen equal to the unit. A determinant of a covariant metric tensor gap, which corresponds to (1) equals to det(gap) = g = -L2r4sin20.
We shall introduce the tetrads as four orthonormal basis vectors which base on the given metric (1):
g(0) = ; g(1)M = LSI + 1F Sf; (2)
g(2)M = -^2(sm + isin0Sl); g(3)M = -^2(Sf - isinesl); (3)
gfc) = L-1 Sf; g£) = S£ + 2 FL-1 SI; (4)
*alex_m_bar@mail.ru; AMBaranov@sfu-kras.ru © Siberian Federal University. All rights reserved
= r>? + sib*'; gfc = - («>
where i is an imaginary unit and the greek indexes run over 0,1,2,3.
The Einstein equations in a tetrad form with a source in the form of the energy-momentum tensor (EMT) are written as
G(a)(/3) = R(a)(/3) — 29(a)(i3)R = — &T(a)(íЗ), (6)
where G(a)(p) are the tetrad components of the Einstein tensor; R(a)(p) are the tetrad components of the Ricci tensor; R = R("\a) is the scalar curvature; ^ = 8n is the Einstein constant in the select system of units. The EMT of the perfect pascal fluid is written down in this case as
T(a)(p) = T(aU03) = (M + V)u(a)u(p) — V9(a)(p) = MU(a)U(P) + Vh(a)(p), (7)
where ^(r) is a mass-energy density; p(r) is a pressure of the perfect pascal fluid; U(a) = g(a)?(dx?/ds) is the 4-velocity in the tetrad form; i{a){p) = U(a)U(p) — g(a)(p) is a 3-projector onto the spacelike 3-surface in the 4-dimensional world (3-metric); all functions depend from radial variable only. This 3-metric is orthogonal to the 4-velocity b(a)(p)U(a) = 0.
We will rewrite Einstein's equations (6) in a form allowing to use properties of the energy-momentum tensor after using a connection between the scalar curvature and a trace of the energy-momentum tensor R = ieT :
R(a)03) = —^{T(o)(p) — 2g(a)(p)T^j . (8)
This system of equations are rewritten as the system of three equations after a substitution of the metric functions into the system
X (ln L)' = 2(m + p); (9)
x 2
^(lnL)' — 2 (F + X(lnF)' — (lnF)'(lnL)') = —xp; (10)
x 2 y Fx J
—X2(1—e)+X (ln L) =—x —p), (11)
where all derivatives are taken with respect of variable x = r/R, R is an exterior ball radius; x = ^R2.
The functions: e(x), F(x) and L(x) are connected as
*(*) = ^. (12)
Such connection is appeared from an equality when the relational volume equals zero in a
static case, u?.? = 0.
it1
Now excluding the mass-energy density and the pressure from this system of equations, we get the linear differential equation with variable coefficients for function G(x):
G" + f (x)G' + g(x)G = 0, (13)
where G = %/F, f (x) = (ln y)', y(x) = y/e/x, and a coefficient g(x) is equal to
g(x) = 2(1 -X2 + *. (14)
If we will input a new variable Z = Z(x) according to a relation
dz = (15)
Ve(x)
then the equation (13) goes over in the equation of a nonlinear spatial oscillator with a variable "frequency" Q(Z(x))
G'Zz + Q2(Z (x))G = 0. (16)
As the equation (16) cannot be integrated, generally speaking, in elementary functions therefore the "frequency" Q can be rewritten as derivative
"2 = -d ($) • "7)
where y = x2, and a function $ is an analog of Newton's gravitational potential of an interior part of the fluid ball. The function $ is found out of the gravitational equations through the function e as
$ = 1 - e = X J ^(x)x2dx = 2Qy J ^(y)^ydy. (18)
Besides we can easily find the relation for the pressure out of the system of gravitational equations (9)-(11)
$1
XP = -^ + -(1 - $)(lnF)'. (19)
Now we will write down the function Q2 as an expansion into series with respect to the variable y
Q2(y) = ^ anyn, (20)
n=0
here y < 1.
Further we will choose all an equal to the zero, i.e. Q2 = 0. Then (16) will be transformed into the equation
G'lz = 0 (21)
and function G(Z) = C1 • y + C2, where C1 ,C2 are constants of an integration.
In this case we have a homogeneous model of gravitational ball with a mass density ^0 = const out of the relations (16) and (17), if Q2 = 0. The
using of the connection between functions e, F and L (12) and also the condition of a joint with the exterior solution of Scharzschild [4] leads to the well-known interior solution of Scharzschild [4].
Thereby we must remark that gravitational field of a spheric static distribution of the substance belongs to the algebraic type D or the type 0 of Petrov's classification [1] in an accordance with the theorem [5].
In our case the gravitational field of the homogeneous distribution of a fluid has the algebraic type 0, i.e. the field is described by the conformally-flat solution of the Einstein equations.
We will put n = 0 in the expansion (20) and now we have the constant "frequency" Q2 = = const. If Q0 is positive (Q2 > 0) then a general solution of equation (15) can be written down as a harmonic oscillating function
G(( (x)) = G0cos(Q0C (x) + M (22)
where y0 is a phase displacement.
We must here remark that the function G(() describes spatial oscillations. On the other hand we find the mass density out of (17) and (18)
m(x) = Mo(1 — bx2), (23)
which is the parabolic distribution, and function G takes on the form
G(x) = Go cos(Qo(ln(deM)/(2TC4 + V^x))) + <Po), (24)
dx
where e(x) = 1 — C3x2 + C4x4; C3 = \R2/3; C4 = xR2b/5.
The using of gravitational equations with EMT of the perfect fluid for the given metric (1) simplifies a calculation of the Weyl tensor components. This tensor can be written in the tetrad components in general case as
W(a.m(7)(S) = R(a.m(7)(S) + R(7)[(a.)g(p)](S) — R(S)[(a)g(p)](7) — 3 Rg(l)[(a)g(P)](S), (25)
where the square brackets denote an antisymmetric operation with the indexes. Mapping the Weyl tensor by
Qf = ¿fk § — iekmnSm Sn (26)
onto the 3-dimensional euclidean space with the metric e^ = diag(1,1,1), we will get the traceless 3 x 3 Weyl matrix (where ekmn is an antisymmetric symbol Levi-Civita)
W = (2R2) x2 x diag(2, —1, —1). (27)
This matrix is a canonical form of Weyl's matrix of the algebraic type D according to the algebraic classification of spaces for the researched case and for any x = 0 (n = 2m/R is a compactness, m is Schwarzschild's mass of the ball).
In the point x = 0 (at the center of ball) we have a degeneration of the Weyl matrix of D type into 0 type, i.e. a continual transition x ^ 0 be accompanied by a jumplike change of the algebraic type of space (D ^ 0). This change corresponds to the rank change of Weyl's matrix from r = 3 to r = 0. Herewith the mass density of fluid is constant and equals mo into a neighborhood of point x = 0. In other words the substance is not compressible nearby of the ball center and the algebraic type of gravitational field is the type 0 (the conformally-flat space-time).
Thus the obtained result is in the full correspondence with the statement of the theorem [5]. Moreover, the algebraic classification of gravitational fields is connected to the solution of a cubic characteristic equation
A3 + pA + q = 0, (28)
where the parameters p and q can be found for the given traceless matrix W as
p = — 1 SpW2; q = — 1 SpW3. (29)
The equation (28) can be considered as the extremum condition of some "potential" function
V (X,p,q ) = Vo + 1 A4 +2 pA2 + q\
(30)
where V0 = const.
A deformation (30) describes the cusp catastrophe of Whitney in accordance with the catastrophe theory. Solutions of the equation (28) are, generally speaking, in three areas, on three lines and in one point of plane of the control parameters p and q when we have the static case (see Fig. 1).
Fig. 1. The cusp catastrophe's surface and its projection onto the plane of the control parameters p and q
It is known that the algebraic type of Riemann's space-time can be changed by the infinitesimal perturbations (herewith the rank of Weyl's matrix varies by jump) [6]. The cusp (p = q = 0) accords with the phase transition into type 0 as into a most symmetric "phase" from point of view of the phase transitions of the second kind (Landau's phase theory [7]). Furthermore the parameter p is an analog of a temperature; a derivative dV/dp is the analog of the entropy; d2V/dp2 is the analog of the thermal capacity. Last derivative varies by jump in the cusp (when p = 0 ) also as the thermal capacity.
In our case the equation (28) has two real coincident roots, the sum of which equals the third root with an opposite sign: À2 = À3 = a0 = (n/(2R2))x2; A1 = -2a0. The parameters p and q are expressed over through a0 as p = 3a2; q = 2ajj. A discriminant Q = (p/3)3 + (q/2)2 = 0 describes a semi-cubic parabola (a projection of the fold lines), which corresponds to space-time of D type. The research case corresponds to a curve marked on Fig.1 by a cross because q > 0.
Thus a state of system has not stability in the point of an inflection of function V(A,p,q) with À = a0.
The extremal values of function V(À, p, q) can be written down for all found roots of the equation (28) as
V (A2 ,p) = V (A3 ,p) = 12 P2 + C ;
2
V (Al ,p) = - 3 P2 + C,
(31)
where q and a0 are expressed through the parameter p. There are the jumps of the second derivatives in respect of p from the function V in the point p = q = 0 : A(Vpp) = 1/6 and A(Vpp) = -4/3. These jumps correspond to the rank jump of the Weyl matrix in the cusp [6].
Hence the continual changes of the variable x (x ^ 0) and of parameters p,q (p ^ 0, q ^ 0) lead to the catastrophe: the algebraic type of space-time has the jump at the center of our ball model (D ^ 0).
From the viewpoint of the physical meaning there is the phase transition of second kind. The eigenvalues A are the parameters of the order as in Landau's phase theory [7].
If we will rewrite the components of the Weyl matrix through the mass density then we can write down Weyl's matrix (27) by way of
W = (2R2) (1 - eg) * '»age,-1,-1). (32>
Here one can see the transition of matrix W of the type D into the matrix of type 0 when the mass density e(x) tends toward the constant value e0.
The consideration of case with n =1 leads to the expression Q2 = a0 + a1y according to (20). Then the corresponding distribution of the mass density can be written as (see (17)-(18))
e(x) = ^0(1 - bx2 - cx4), (33)
where c is some constant.
The further analysis carried out similarly in a neighbourhood of the gravitating ball center leads to a previous result of the algebraic types' change by jump. The gravitational field has the algebraic type 0 at the center (and into its small neighborhood) of the ball and in the remaining part of ball the gravitational field belongs to the type D. In the expansion (20) there are not any new results for n > 1.
Summary
In the article the problem of an existence of another algebraic type at the fluid ball center than in the remaining part of ball is considered. The perfect fluid as the source of gravitational field has the arbitrary spherical distribution of the mass density. The eigenvalues task of Weyl's matrix is connected with the algebraic classification of Petrov. This classification is the algebraic classification of gravitational field in the Riemann space-time. At the center of ball the Weyl matrix of spherical distribution of the mass density equals to the zero, i.e. at the center the algebraic type is 0 (the conformally-flat space-time).
On the other hand the theorem [5] says that the spherical gravitational field can be only either the algebraic type D or the type 0 according to Petrov's algebraic classification. In general case the interior gravitational field of our fluid ball has the algebraic type D outside of the center of ball, where we have the type 0.
We can see from Fig.1 the cusp corresponds to the conformally-flat space-time (the algebraic type 0). This correspondence is strong and the algebraic type is changed by jump there. It is the catastrophe or the phase transition of second kind in the physical meaning. The variables Aj (see (28)) play here the role of the order parameter according to the phase theory of Landau [7].And the algebraic types of the Riemann space-time play the roles of "the substance phases" (as in Landau's theory).
In other words into the small neighborhood of the ball center there is the phase transition of the second type in the spheric gravitational field of the fluid ball for any choice of n (n > 0) in (20) as the algebraic type change of Riemann's space-time.
References
[1] A.Z.Petrov, New Methods in General Relativity, Nauka, Moscow, 1966 (in Russian).
[2] A.M.Baranov, Vestnik of Krasnoyarsk State University, Phys. Math. Sci., (2002), no. 1, 5-12 (in Russian).
[3] A.M.Baranov, Z.V.Vlasov, Vestnik of Krasnoyarsk State University, Phys. Math. Sci., (2005), no. 1, 4-11 (in Russian).
[4] J.L.Synge, Relativity: the General Relativity, North-Holland Publishing Company, Amsterdam, 1960.
[5] A.M.Baranov, Vestnik of Krasnoyarsk State University, Phys. Math. Sci., (2006), no. 1, 5-9 (in Russian)
[6] A.M.Baranov, Dep. in VINITI USSR (University of People's Friendship), no. 2632-76, July 13, (1976) (in Russian).
[7] L.D Landau, E.M.Lifshitch, The Statistical Physics, Part 1, Nauka, Moscow, 1976 (in Russian).
Проявление гравитационного фазового перехода в жидком статическом шаре
Александр М. Баранов
Показано, что в центре гравитирующего жидкого шара присутствует фазовый переход 2-го 'рода гравитационного поля как изменение алгебраического типа поля..
Ключевые слова: гравитационный фазовый переход, уравнения Эйнштейна, общая теория относительности и гравитации, алгебраическая классификация Петрова, теория катастроф.
