VJV

UNIVERSIDAD NACIONAL DEL CALLAO

FACULTAD DE CIENCIAS NATURALES Y MATEMÁTICA

ESCUELA PROFESIONAL DE FÍSICA

"SOLUCIÓN A LAS ECUACIONES DE CAMPO GRA VITACIONAL DE

EINSTEIN PARA UN FLUIDO ANISOTRÓPICO CON SIMETRÍA
ESFÉRICA"

TESIS PARA OPTAR EL TÍTULO DE PROFESIONAL

DE LICENCIADO EN FÍSICA

MIGUEL ANGEL DE LA CRUZ CRUZ

CALLAO-PERU

ENER0-2011
 HOJA DE PRESENTACIÓN

"SOLUCIÓN A LAS ECUACIONES DE CAMPO GRA VITACIONAL DE EINSTEIN

PARA UN FLUIDO ANISOTRÓPICO CON SIMETRÍA ESFÉRICA"

Bach. MIGUEL ANGEL DE LA CRUZ CRUZ

Tesis presentada a consideración del cuerpo docente de la Facultad de Ciencias Naturales y

Matemática de la Universidad Nacional del Callao, como parte de los requisitos para
obtener el Titulo Profesional de Licenciado en Física.

Aprobada por:

Lic. Rolando Juan Alva Zavaleta

Mg. Juan Abif!!!:Ldez Velásquez

Lic. Carlos Alberto Lévano Huamaccto

Lic. Rolando Manuel Vega De La Peña

Callao - Perú

Enero- 2011
 FICHA CATALOGRÁFICA

DE LA CRUZ CRUZ, MIGUEL ANGEL

Solución a las Ecuaciones de Campo Gravitacional de Einstein para un Fluido

Anisotrópico con Simetría Esférica, Callao (2011)
X, 105p, 29.7cm. (UNAC, Licenciado en Física, 2011)

Tesis, Universidad Nacional del Callao, Facultad de Ciencias Naturales y

Matemática.
Física
l. UNAC/FCNM 11. Título (serie)
Dedicado a Euclides quien, gracias a la imposibilidad
de la demostración de su V postulado, motivó el
descubrimiento de nuevas geometrías, siendo el
primer beneficiado Albert Einstein.
 AGRADECIMIENTOS

Expreso mi gratitud a todos mis familiares, en especial a mis padres, por su

compresión y estímulo durante mis estudios escolares y, en especial, en esta difícil carrera

de física. Gracias a ellos encuentro motivo para seguir trabajando.

Quiero también agradecer de forma especial a mi asesor de tesis, Lic. Rolando

Manuel Vega De La Peña, por el apoyo incondicional y confianza en el desarrollo de la

tesis. Asimismo, quiero agradecer a. los docentes y amigos Fulgencio Villegas Silva y
Aleja~dro Tuesta Velásquez, como los motivadores y promotores para que mi persona se

dedique a la línea de física teórica y computacional, sin sus consejos y recomendaciones

me hubiese sido. difícil empezar la presente tesis. Hago extensiva mi gratitud al Decano Mg.

Roel Mario Vidal Guzmán, al Director de Escuela Lic. Rolando Juan Alva Zavaleta, al Jefe

de Laboratorio de Física y Química Mg. Juan Abraham Méndez Velásquez, al Director de

Instituto de Investigación Blgo. Edgar Zárate Sarapura, presidentes, miembros de las

comisiones académicas y profesores pertenecientes a la FCNM y a mis ex-alumnos Danny

Vicente Vargas, Pablo Hurtado Quispe, Saúl Sánchez Juárez, Leonardo Briseño Portilla,

Alexander Galván López, Melisa Cárdenas Ochoa, Enrique Urbano Tunque, Daniel Ortiz,

Christian Gómez Luca, Luis Ángel Tovar De La Cruz, Zusana Cecilia Verástegui Muñoz y
a todos mis ex-alumnos de Física I de la Escuela de Matemática y Física II y Física III de la

escuela de Física y Matemática desde el 2005 hasta la fecha.

No olvido en especial a mis mejores amigos, hoy colegas, Fernando Salazar

Espinoza y Roxana Elizabeth Huamán Urquiza con quienes, en momentos buenos y malos,

estuvimos siempre unidos, tratando de hacer las cosas correctas siempre. con la frente en

alto desde que fuimos Tercio Estudiantil (2002-2003), enseñando a nuestros compañeros a

tener mayor identificación hacia nuestra Facultad, y compartiendo numerosas experiencias

en todos estos años. Gracias a ellos pudimos ejecutar proyectos académicos dirigidos a la

comunidad chalaca y hacia nuestros compañeros de la Facultad,· siendo el principal y

primordial beneficiario la Facultad de Ciencias Naturales y Matemática. A todos ellos mi

eterno y sincero agradecimiento.

 RESUMEN

"SOLUCIÓN A LAS ECUACIONES DE CAMPO GRA VITACIONAL DE EINSTEIN

PARA UN FLUIDO ANISOTRÓPICO CON SIMETRÍA ESFÉRICA"

Bach. MIGUEL ANGEL DE LA CRUZ CRUZ

ENERO- 2011

Asesor :Lic. Rolando Manuel Vega De La Peña

Título obtenido : Licenciado en Física

En este trabajo se resuelven las ecuaciones de Einstein para un fluido anisotrópico con

simetría esférica con la métrica estándar de Schwarzschild y, a partir de las funciones

generatriz y anisotrópica, se obtienen los coeficientes métricos que describen el

comportamiento de este fluido en la geometría del espacio-tiempo.

Las soluciones de campo que se obtienen permiten analizar el comportamiento de dicho

fluido a través de sus variables termodinámicas, tales como densidad y presión.

Finalmente, se analizan algunos ejemplos con las funciones generatriz y anisotrópica que

describen el comportamiento local del fluido anisotrópico.

PALABRAS CLAVES

FLUIDO ANISOTRÓPICO

FUNCIÓN GENERATRIZ

FUNCIÓN ANISOTRÓPICA

COEFICIENTES MÉTRICOS
 ABSTRACT

"SOLUTION TO THE EQUATIONS OF GRA VITATIONAL FIELD THE EINSTEIN'S

FOR A FLUID ANISOTROPIC WITH SPHERICAL SYMMETRY"

Bach. MIGUEL ANGEL DE LA CRUZ CRUZ

ENER0-2011

Advisor : Lic. Rolando Manuel Vega De La Peña

Obtained title : Graduate in Physics

In this work was solved Einstein's equations, for a flowing anisotropic with spherical

symmetry with the metric standard of Schwarzschild and, starting from the generatriz and

anisotrópic functions, it was obtained the metric coefficients that it describes the behavior

ofthis fluid on the geometry ofthe space-time.

The field solutions that were obtained allowed to analyze the behavior, through their

thermodynamic variables as density and pressure.

Finally, it was analyzed sorne examples with the generatriz and anisotropic functions that it

describes the local behavior of flowing anisotropic.

KEYWORDS

FLUID ANISOTROPIC

FUNCTION GENERATRIZ

FUNCTION ANISOTROPIC

METRIC COEFFICIENTS
 ÍNDICE

CAPÍTULO!

1 Introducción 1

CAPÍTULOII

11 Fundamentos teóricos 4

11.1 Geometría Riemanniana 4

11.1.1 Transformación de coordenadas 4

. 11.1.2 Vector contravariante 5

11.1.3 Vector covariante 5

11.1.4 Tensores 6

11.1.5 Tensor métrico 7

.. 11.1.6 Álgebra tensorial 11

11.1.7 Símbolos de Christoffel 13

11.1.8 Ley de transformación de los símbolos de Christoffel 14

11.1.9 Derivación covariante de vectores 14

11.1.1 O Derivación covariante de tensores 17

11.1.11 Tensor de curvatura de Riemann 19

II.l.12 Tensor de Ricci 21

II.l.13 Tensor de Einstein 21

 II.l.14 Identidades de Bianchi 22

II.l.15 Geodésicas 22

11.2 Teoría de la relatividad general 24

11.3 Principio de equivalencia 25

11.4 Principio de covarianza general 26

11.5 Ecuaciones de campo gravitacional de Einstein 27

11.6 Soluciones exactas a las ecuaciones de Einstein 29

II.6.1 Solución de Schwarzschild 30

II.6.2 Solución de Robertson-Walker 31

II.6.3 Solución de Reissner-Nordstrom 32

11.6.4 Solución de Kerr 33

II.6.5 Solución de Kerr-Newman 25

CAPÍTULO 111

111 Fluido en relatividad general 35

111.1 Fluido 35

111.2 Fluidos en general 36

111.3 Fluido perfecto 37

111.4 Hidrodinámica clásica 38

111.5 Tensor hidrodinámico relativista 39

111.6 Ley de conservación de la energía-momentum 40

 CAPÍTULO IV

IV Ecuaciones de campo 43

IV.1 Elemento de línea 43

IV.2 Ecuaciones de campo de Einstein 48

· IV.3 Ecuación de equilibrio hidrostático 58

CAPÍTULO V

V Solución de las ecuaciones de campo para una esfera anisotrópica 61

V.1 Cálculo de los coeficientes métricos 61

CAPÍTULO VI

VI Verificación y análisis de la solución obtenida 69

VI.1 Solución global 69

VI.2 Solución local 69

VI.3 Solución de la métrica inhomogénea de Friedmann 81

CAPÍTULO VII

VII Conclusiones 94
 CAPÍTULO VIII

VIII Referencias bibliográficas 98

APÉNDICES

A. Símbolos para las cantidades físicas 101

B. Constantes físicas 103

C. Escalas y unidades naturales 104

CAPÍTULO!

INTRODUCCIÓN

Desde el principio el hombre ha tenido la necesidad de explicar el mundo

que lo rodea, haciéndose las siguientes preguntas: ¿cómo surge el universo?, ¿por

qué los planetas giran alrededor del Sol?, etc. Todas estas preguntas son abordadas

por diferentes áreas tales como la astronomía, cosmología y gravitación. Hasta la

fecha se siguen realizando muchos trabajos en esos campos, llegando inclusive a

reformular, precisar y corregir algunas teorías y/o medidas de los cuerpos celestes,

astros, etc.

Han pasado más de 90 años desde que la teoría de la gravitación de Einstein

fuera concluida y plasmada en el memorable artículo clásico titulado

"Fundamentos de la Teoría de la Relatividad Generar', publicado por la revista

Annalen der Physik, el 20 de marzo de 1916 [1,24,25]. A lo largo de estos años,

numerosas predicciones de la teoría gravitacional de Einstein han sido objeto de

comprobación experimental, verificándose cada vez con mayor exactitud y

obteniéndose así nuevos elementos que la: sustenten. Son numerosas las aplicaciones

de esta teoría en la astronomía, astrofísica, cosmología, etc. [19,20,21,22,23].

Los pnmeros modelos relativistas predijeron que un umverso, que era

considerado como materia gravitacional ordinaria dinámica y contenida, tendería a

la contracción más bien que a la expansión. Para este modelo, Einstein agrega la

1
llamada constante cosmológica con el fin de solucionar un problema donde se

implicaba que el balance era hacia fuera en la contracción y, de esa forma, obtener

una solución estática del universo.

En 1922, Alexander Friedmann aportó las hoy famosas ecuaciones de

Friedmann [8], verificando que el universo es materia dominante y que la presión

es negligible en comparación con la densidad. Con ello se pudo explicar muchos

comportamientos de algunos cuerpos celestes en expansión o contracción. De otro

lado, el primero en desarrollar la ecuación de Schwarzschild para un fluido

anisotrópico fue R. Tolman [8], como una ampliación a los trabajos de R.

Openheimer. Posteriormente, tuvieron que pasar cerca de 50 años para que S.

Bayin [14] ampliara las soluciones de Tolman y las aplicara para cuerpos

anisotrópicos y cuerpos celestres con radiación.

El objetivo de este trabajo es analizar el comportamiento de un fluido

anisotrópico con simetría esférica utilizando la métrica de Schwarzschild, la cual se

construirá a partir de las funciones generatriz y anisotrópica, para que con ellas se

puedan obtener las coeficientes métricos que describen el comportamiento de este

fluido en la geometría del espacio-tiempo y así, finalmente, encontrar las variables

termodinámicas tales como densidad, presión, etc. para un espacio-tiempo curvo,

como también para un espacio-tiempo plano. Luego se procede a realizar la

comprobación de la solución, señalando con algunos ejemplos las aplicaciones de la

misma.

2
 Como ejemplo adicional se analiza este fluido con la métrica inhomogénea

de Friedmann anulando la función anisotrópica, para así ver el comportamiento de

la densidad, presión y conducción de calor dependiente de las coordenadas espacial

y temporal. Finalmente, se analizan algunas gráficas y soluciones para verificar que

en los extremos el espacio-tiempo se haga plano y encontrar los límites de las

variables termodinámicas con respecto a las coordenadas espacial y temporal, cuyo

análisis detallado se muestra en las conclusiones de esta tesis.

3
 CAPÍTULOII

FUNDAMENTOS TEÓRICOS

11.1 GEOMETRIA RIEMANNIANA

11.1.1 TRANSFORMACIÓN DE COORDENADAS

Para el estudio de campos gravitatorios es .necesario el estudio de coordenadas

curvilíneas cuadridimensionales arbitrarias [2,3]. Una transformación de

coordenadas, es decir, el paso de un sistema de coordenadas x 11 , con f.1 = 0,1,2,3, a

otro x· 11 , está dada por

(2.1)

La transformación se lleva a cabo entre las coordenadas x· 11 y xv, cuando las f 11

son funciones independientes, reales y poseen derivadas parciales, y su jacobiano

sea diferente de cero.

a¡ o a¡3
0
8x 8x 0
,fl ,O ,1 2 ,3 11
J(~) = J(x ,x ,x .,x ) = 8/ = :;t:Ü (2.2)
xv 0 1 2
x x x x 3 axv
' ' '
a¡ o a¡3
8x 3
8x 3

4
 Como el jacobiano es diferente de cero, entonces existe una transformación inversa

de coordenadas:

(2.3)

11.1.2 VECTOR CONTRAVARIANTE

Es todo conjunto de cuatro cantidades Ufl (¡.¡, = 0,1,2,3) que, en una transformación

de coordenadas xfl ~ x·fl , se transforman mediante:

(2.4)

donde, U f1 y U ,fl son las componentes del vector contravariante.

Por ejemplo, para la diferencial de coordenadas, a partir de las reglas de

diferenciación parcial se expresa como:

(2.5)

Entonces, la diferencial de coordenadas es un vector contravariante.

11.1.3 VECTOR COVARIANTE

Es todo conjunto de cuatro cantidades Vfl ( fl = 0,1,2,3) que, en una transformación

de coordenadas xfl ~ x·fl , se transforman mediante:

(2.6)

5
 donde, V.u y V· .u son las componentes del vector covariante.

Así, sea l.fl un escalar, para una transformación de coordenadas el gradiente es:

(2.7)

Bl.fl .
Entonces - - es un vector covanante.
ax.u

11.1.4 TENSORES

Un tensor contravariante de rango 2, r.uv se define como una colección de 16

cantidades de la forma [26,27 ,28]:

Too TOl To2 To3

TlO T¡¡ T12 T13
r.uv = (2.8)
T2o T21 T22 T23
T3o T31 T32 T33

El producto de vectores contravariantes U .u y vv , es definido como:

(2.9)

Para una transformación de coordenadas x ~ x· se obtiene:

(2.10)

6
 De igual forma, el producto de vectores contravariantes UP y Vv , está definido

como:

(2.11)

Para una transformación de coordenadas x ~ x· se obtiene:

(2.12)

Los tensores mixtos (contravariantes y covariantes) de cualquier orden se definen

generalizando las ecuaciones (2.10) y (2.1 1), mediante la siguiente ecuación:

(2.13)

la cual define un tensor mixto de orden m+n, contravariante de m-ésimo orden y

covariante de n-ésimo orden.

11.1.5 TENSOR MÉTRICO

En un espacio n-dimensional En, el cuadrado de la distancia entre dos puntos está

dado por:

(2.14)

7
 11
Haciendo un cambio de coordenadas, de uno local ; a uno general x 11 , la ecuación

(2.14) toma la forma:

(2.15)

De esta ecuación, ds 2 en el sistema de coordenadas local ; , puede escribirse como

la forma bicuadrática

(2.16)

donde,

(2.17)

1lap es llamado tensor de Minkowski.

De acuerdo a la teoría de la relatividad especial [4,24,25], la distancia entre dos

puntos infinitesimales esta dado por:

(2.18)

Estas coordenadas están referidas a un espacio tetradimencional de Minkowski,

donde; x 0 = e t , es la componente temporal y x 1 = x , x 2 =y , x3 = z, son las

componentes espaciales.

Como ds 2 es un invariante, entonces para x1 = x 2 = x3 = O no hay dependencia

espacial, por lo que t es llamada tiempo propio. Si hacemos t = -r , entonces la

coordenada temporal se expresa como:

(2.19)

8
Considerando el sistema de unidades semi-naturales, es decir haciendo e =1, el

tiempo propio se expresa:

(2.20)

El tiempo propio también puede ser expresado en un sistema de coordenadas

arbitrario mediante:

(2.21)

Si g 11 v es llamado tensor métrico contravariante de segundo orden, definido por:

(2.22)

entonces el tiempo propio se expresa como:

(2.23)

Para una transformación de coordenadas x ~ x· , el tensor métrico toma la

siguiente forma:

(2.24)

Aplicando (2.12) a (2.16)

(2.25)

Finalmente:

(2.26)

9
Se observa que g ¡.¡v es un tensor covariante y simétrico . Su inversa es un tensor

contravariante y simétrico. Estos vectores cumplen con la siguiente relación:

(2.27)

donde o; es llamado delta de Kronecker, y cumple

1 Si: a= v
ga = oxa ={ (2.28)
V QXV
o Si a t:-v

Para una transformación de coordenadas x ~ x· , la delta de Kronecker toma la

siguiente forma:

(2.29)

Analizando la última ecuación se observa una contracción de índices (a =v y

f3 = v ), por lo que:

(2.30)

A partir de (2.18) , (2.19) y (2.20) se demuestra también la siguiente relación:

(2.31)

Esta última expresión verifica que gw es un tensor contravariante.

A partir de (2.22), el elemento· de línea se define de la forma:

(2.32)

10
 Con la ayuda de estos dos tensores fundamentales g pv y gpv, se puede subir y bajar

los índices de un tensor. Así, se pueden obtener tensores de diferentes grados de

contravarianza y covarianza, los cuales están dados por las siguientes ecuaciones:

(2.33)
T·······
... Ji...
=g JiV
T .. v ...
········

11.1.6 ÁLGEBRA TENSORIAL

Para formar otro tensor se deben cumplir ciertas reglas, las cuales se indican a

continuación.

(a) La combinación lineal de dos tensores del mismo tipo es un tensor del

mismo tipo. Sean dos tensores mixtos U: y V: ,entonces deben cumplir:

TP
V
=aUPV +bVP
V
(2.34)

donde a y b son escalares. Entonces r: es un tensor, porque en una

transformación de coordenadas se obtiene:

(2.35)

(b) El producto directo de dos tensores consiste en un tensor con todos los

índices superior e inferior de los tensores originales. Sean los tensores U: y

va 'entonces

(2.34)

11
 r:a es un tensor, porque en una transformación de coordenadas se obtiene:

(2.35)

(e) La contracción de un tensor se produce cuando se tienen los índices arriba y

abajo iguales. Sumando sobre sus cuatro valores produce un nuevo tensor

con estos dos índices ausentes. Por ejemplo,

(2.34)

El T¡.¡v es un nuevo tensor ya que cumple

(2.35)

(d) Un índice contravariante puede ser transformado en un índice covariante,

contrayéndolo a través del tensor métrico g w

g ¡.JV yva =ya

V
(2.36)

(e) Un índice covariante puede ser transformado en un índice contravariante,

contrayéndolo a través del tensor métrico gw

(2.37)

12
11.1.7 SÍMBOLOS DE CHRISTOFFEL

Los r son coeficientes denominados símbolos de Christoffel o coeficientes de

conexión afín. A partir del tensor métrico, estos símbolos vienen dados por:

(2.38)

el cual es llamado símbolo de Christoffel de primera clase. Si se quiere subir índices

se multiplica (2.38) con el tensor métrico, obteniéndose:

(2.39)

el cual es llamado símbolo de Christoffel de segunda clase.

Los símbolos de Christoffel no son tensores, debido a que sobre ellos no cumple la

ecuación (2.13). Estos símbolos al ser simétricos en sus índices p y O" , cumplen

las siguientes propiedades:

(2.40)

(2.41)

(2.42)

(2.43)

(2.44)

13
 Los símbolos de Christoffel, son también escritos según las convenciones de la

siguiente forma:

rapa= {a,p,a} (2.45)

(2.46)

11.1.8 LEY DE TRANSFORMACIÓN DE LOS SÍMBOLOS DE CHRISTOFFEL

Sea el símbolo de Christoffel de primera clase

r· 1 agaa
=- · +~+~
aa· aa· } (2.47)
apa 2 { 8x•P 8x·a 8x·a

Utilizando las propiedades del tensor métrico y de (2.4 7) se demuestra también la

siguiente relación:

(2.48)

Esta última ecuación es la ley de transformación de los símbolos de Christoffel de

primer orden.

11.1.9 DERIVACIÓN COVARIANTE DE VECTORES

En las coordenadas cartesianas las diferenciales de un vector dUP de un vector U f.l

forman un vector y las derivadas __

au
f.l forman un tensor. Pero, para el caso de las
axv

14
 BU11
coordenadas curvilíneas, dU11 no es un vector y no es un tensor. Ello es
axv

axv
debido a que los coeficientes de las fórmulas de transformación (2.6) son
8x 11

funciones de las coordenadas.

Para comparar dos vectores separados infinitesimalmente, se debe trasladar uno de

ellos paralelamente a sí mismo hasta el punto donde está localizado el segundo.

Esto es llamado transporte o traslación en paralelo. Consideremos para la

coordenada x 11 el vector contravariante U 11 , para un punto vecino x 11 + dx 11 el

vector es de la forma U 11 + dU 11 • Al ser sometido el vector U 11 a una traslación

hasta el punto x 11 + dx 11 , el cambio que experimenta lo denotamos por &P . La

diferencia DU 11 entre los dos vectores, que están aplicados al mismo punto, es

entonces:

(2.49)

8U 11 , es llamado el desplazamiento paralelo de U y es definido como:

(2.50)

Para un sistema de coordenadas curvilíneas, los coeficientes o símbolos de

Christoffel deben cumplir

(2.51)

En cambio, para un sistema de coordenadas cartesianas

r~ =O (2.52)

15
Sabemos que en una traslación un escalar no varía; entonces, como el producto

vectorial de dos vectores U 11 y V 11 es un escalar, se tiene

(2.53)

Desarrollando la ecuación anterior se obtiene lo siguiente:

Dada la arbitrariedad de V 11 ,

(2.54)

Esta fórmula determina la variación de un vector covariante en una traslación.

Sustituyendo (2.54) y dU 11
au 11
= --dxv en (2.49), se obtiene:
axv

(2.55)

De igual forma para un vector covariante se tiene

(2.56)

16
 Se observa que en (2.55) y (2.56) la expresión entre paréntesis es un tensor ya que

permite la transformación de un vector dxfl a un nuevo vector DU. Estos tensores

son llamados derivadas covariantes de los vectores U f1 y U f1 , representados por

U(:, y Up;v. Entonces DUfl y DUfl se expresan como:

(2.57)

(2.58)

donde las derivadas covariantes para un vector contravariante y un vector covariante

son:

(2.59)

(2.60)

Para r :v = O , la derivada covariante se reduce a la derivada ordinaria.

11.1.10 DERIVACIÓN COVARIANTE DE TENSORES

Sea un escalar de la forma TpvUflVv, entonces el desplazamiento de este escalar es:

(2.61)

Al desarrollar la ecuación anterior se obtiene:

17
 (2.62)

Como U 11 Vv ::t:- O, entonces la expresión entre paréntesis es igual a cero.

Despejando términos se obtiene:

(2.63)

Si se analiza el cambio OI'11 v como un desplazamiento paralelo del tensor, se

tendrían dos tensores T11 v(x 11 ) y T;v (x 11 + dx 11 ) respectivamente. Utilizando (1.63)

y restando estos dos tensores se obtiene:

(2.64)

Además, para un cambio infinitesimal de coordenadas

(2.65)

Reemplazando (2.65) en (2.64) se tiene:

(2.66)

La expresión entre paréntesis es llamada derivada covariante de un tensor

covariante de segundo orden, la cual es denotada en la forma:

18
 (2.66)

Igualmente, aplicando los criterios anteriores, la derivada covariante de un tensor

contravariante de segundo orden es:

(2.67)

En general para un tensor mixto de cualquier orden se tiene:

arpv ...
ypv...
ap ... ;A.
= ap...
axA.
+ 111 yav ... + rv
aA. ap...
Tpa ...
aA. ap...
+ ... - rmaA. yw ... ... - rm yw ... - ...
mp... pA. am .. .
(2.68)

11.1.11 TENSOR DE CURVATURA DE RIEMANN

Para la descripción de un espacio-tiempo curvo es necesario trabajar con un tensor

de curvatura, aplicando para ello los criterios de derivación covariante desarrollados

anteriormente [4,5].

Si se aplica la derivada covariante a Va;p se obtiene:

(2.69)

19
Cambiando el orden de los índices fJ y r en la ecuación (2.69), se obtienen dos
derivadas covariantes va;fJ;r y va;y;fJ ' cuya resta es:

(2.70)

La ecuación anterior puede ser escrita como:

(2.71)

donde el tensor R;p;. está dado por:

(2.72)

El tensor R;p;. es llamado tensor de curvatura de Riemann.

De otro lado, se puede apreciar que este tensor de curvatura satisface las siguientes

relaciones:

(2.73)

(2.74)

(2.75)

(2.76)

20
11.1.12 TENSOR DE RICCI

A partir del tensor de curvatura de Riemann, aplicando contracción, se obtiene:

R aj3 -- R)l - g)lVR

a)l/3 - )lavf3 (2.77)

donde Ra/3 es denominado tensor de Ricci.

Usando las ecuaciones (2.72) y (2.73), el tensor de Ricci satisface las siguientes

relaciones:

(2.78)

(2.79)

Contrayendo Ra/3 se tiene:

R = Raa = gaf3 R af3 = gaf3g)lV R )la~

(2.80)

donde R es llamada curvatura escalar de Ricci.

11.1.13 TENSOR DE EINSTEIN

El tensor de Eisntein es definido como:

(2.81)

Se verifica también que la derivada covariante de GJlV es:

(2.82)

21

231.1_
11.1.14 IDENTIDADES DE BIANCHI

La derivada covariante de R:Pr;v se reduce a:

(2.83)

A partir de (2.83) se verifica que:

(2.84)

Esta última es conocida como las identidades de Bianchi, siendo sus formas

contraídas las siguientes:

(2.85)

(2.86)

(2.87)

(2.88)

11.1.15 GEODÉSICAS

La ecuación diferencial de una curva que tiene una longitud extremal es llamada la

ecuación geodésica. Para describir esta ecuación, se debe analizar las relaciones que

satisfagan el valor estacionario de la integral:

22
 1= fds (2.89)

donde los límites de integracíón son tomados a partir de dos puntos fijos y

próximos. Por lo tanto, se debe buscar la solución al problema variacional

(2.90)

donde el Lagrangiano L está dado por:

L -( g
-
-
pv ds ds
dxvJ2
dxJ.l - (2.91)

y su valor es igual a la unidad a lo largo de la curva geodésica.

Usando (2.90) y el cálculo variacional, se obtiene la ecuación:

(2.92)

Como IJxP ds ;:¡:. O, se encuentra la ecuación de Lagrange:

!!_ ( aL
ds a(dxJ.l 1ds)
J_(JxJ.l
aL _0
-
(2.93)

A partir de la ecuación anterior, y utilizando el elemento de línea ds 2 = gpvdxpdxv

a lo largo de la línea geodésica, se obtiene:

23
 (2.94)

Esta ecuación es conocida como la ecuación de la geodesica.· Si se analiza

localmente la curvatura, entonces

(2.95)

11.2 TEORÍA DE LA RELATIVIDAD GENERAL

Las ideas físicas que descansan en los fundamentos de la teoría de la relatividad

general fueron formuladas por Einstein inmediatamente después de la creación de

la teoría de la relatividad especial en 1907 [1]. Analizando los procesos físicos en un

sistema no inercial, Einstein estableció que en estos sistemas no es posible dar la

misma definición de la medición del tiempo (simultaneidad) que fue dada para el

caso de sistemas inerciales debido a que, en un sistema no inercial, el tiempo varía

de un punto a otro y depende de la aceleración en un elemento dado del espacio, es

decir de la influencia del potencial gravitacional.

En sus investigaciones físicas, Einstein partió del hecho que la geometría depende

directamente de los procesos físicos y, por ello, no representa una forma de

descripción elegida por conveniencia a la cual la física debía adaptarse. Las leyes de

la geometría, desde este punto de vista, deben considerarse como leyes físicas

integrales; es decir, en la teoría de la relatividad general no existe una geometría y

una cinemática independiente de los procesos físicos, ya que las propiedades de las

escalas y de los relojes se determinan por el campo gravitacional y este campo,

expresado a través de las componentes de los potenciales, es el que representa el

24
 estado físico del espacio, el cual define simultáneamente la gravitación, la inercia y

la métrica.

En la teoría de la relatividad general, de acuerdo a la interpretación de Einstein, el

movimiento de los cuerpos ya no es debido a fuerzas aplicadas sobre él, si no a la

geometría del espacio-tiempo, que está curvo o deformado por la presencia de esta

materia.

11.3 PRINCIPIO DE EQUIVALENCIA

De acuerdo a Einstein, el principio de equivalencia se puede formular con el

siguiente enunciado:

"En cada punto del espacio-tiempo, en un campo gravitacional arbitrario, es

posible elegir un sistema de coordenadas localmente inercial tal que, dentro de una

región suficientemente pequeña del punto en cuestión, las leyes de la naturaleza

toman la misma forma que en los sistemas coordenados cartesianos sin aceleración

en ausencia de la gravitación".

Es decir, localmente el campo gravitacional puede ser considerado como un sistema

acelerado. Este principio permite explicar que la aceleración ya no está asociada a la

fuerza, sino a la gravitación.

Una aplicación del principio de equivalencia planteada por Einstein es el siguiente

experimento ideal:

Consideremos un sistema S que está sometido a un campo gravitacional uniforme

~ ~

g. Todas las partículas que están en este sistema caen con la misma aceleración g,

25
 siendo la característica de los cuerpos su masa gravitatoria mg. Ahora

reemplazamos el sistema por un sistema S· que está acelerando uniformemente con

una aceleración - g . Los cuerpos en este sistema están en movimiento relativo,

siendo la característica de los cuerpos la masa inercial m; y, en este sistema, los

cuerpos están experimentando una aceleración g , lo mismo que el sistema S .

Así, la mecánica dentro de S· es equivalente a la mecánica dentro de S , lo cual es

válido si cumple la siguiente igualdad:

(1.96)

Esta última expresión indica que es imposible para un experimento fisico encontrar

una distinción fisica entre los sistemas S· y S.

11.4 PRINCIPIO DE COVARIANZA GENERAL

El aparato matemático desarrollado por Christoffel, Ricci y otros, mediante el cual

se obtuvo un cálculo diferencial independiente del sistema de coordenadas que

permitía dar una forma invariante a las ecuaciones de la fisica-matemática, fue

utilizado por Einstein llevando sus ecuaciones en forma diferencial; es decir, las

lleva en una representación de covarianza general, siendo éstas invariantes con

respecto a cualquier transformación de las coordenadas espacio-temporales.

Todas las leyes físicas deben estar descritas por ecuaciones covariantes. Este

principio no impone restricciones en el contenido de las leyes físicas, sino en la

forma como ellas se escriben. Una ecuación se dice covariante bajo una

26
 transformación de coordenadas si la forma de la ecuación no cambia por dicha

transformación.

11.5 ECUACIONES DE CAMPO GRAVITACIONAL DE EINSTEIN

Para la descripción de un campo gravitacional se usan las ecuaciones de campo

gravitacional de Einstein [5,6]. Estas ecuaciones son una generalización de las

ecuaciones de campo de Newton, que asume la existencia de solo un potencial que

describe el campo gravitacional.

Los potenciales en relatividad general son identificados con 1O componentes del

tensor métrico g 1-JV de la geometría de un espacio-tiempo curvado Riemanniano.

Así, uno debe tener 1O ecuaciones diferenciales parciales de segundo orden para la

métrica g 1-JV, lo cual se obtiene a partir de la ecuación de Poisson:

V 2 <1>(x) = 4nGp(x) (2.97)

donde G es la constante gravitacional de Newton y p(x) es la densidad de materia

que produce el campo gravitacional.

Para un espacio-tiempo plano, el tensor métrico g w es de la forma:

1 o o o
o -1 o o
g!-JV = (2.98)
o o -1 o
o o o -1

27
Utilizando la identidad (2.27) se verifica

1 o o o
g"v =
o -1 o o (2.99)
o o -1 o
o o o -1

Se observa que los elementos de g w y gw son iguales para un espacio-tiempo

plano, pero para un espacio-tiempo curvo esto no es válido.

La relación entre la geometría del espacio-tiempo, la masa y la energía presente en

él, se establece mediante las ecuaciones de campo de Einstein:

(2.100)

donde:

G¡¡v es el tensor de Einstein definido por la ecuación (2.81)

Tw es el tensor energía-momentum que expresa la distribución de materia y

energía.

Contrayendo los índices 11 y v (2.98) la ecuación se reduce a:

R = -87lGT (2.101)

donde gw g¡¡v = 4 y T = r: = g"vTw es la traza de el tensor energía-momentum.

28
 Reemplazando (2.1 O1) en (2.1 00) se obtiene:

(2.102)

Si se analiza en el vacío, T11 v = O, entonces la ecuación anterior se reduce a:

(2.103)

Esta última ecuación es llamada la ecuación de campo de Einstein para el vacío.

Inicialmente, la ecuación de campo gravitacional fue escrita de la forma:

(2.104)

donde Einstein introdujo una constante adicional A, llamada constante

cosmológica, para que la ecuación de campo gravitacional sea compatible con el

postulado de Mach. Ahora bien, toda teoría gravitacional debe reducirse a la teoría

newtoniana, en el caso de movimientos lentos y campos gravitacionales débiles. Al

parecer la constante A debe ser muy pequeña como para interferir con los éxitos de

la teoría de Newton de la gravitación. En lo que sigue se asumirá que A=O.

11.6 SOLUCIONES EXACTAS A LAS ECUACIONES DE EINSTEIN

La teoría de la relatividad general asocia la curvatura del espacio-tiempo con la

presencia de materia, lo cual se muestra en la ecuación (2.1 00). Debido a la

complejidad de las ecuaciones de campo de Einstein, es difícil encontrar soluciones

exactas excepto en espacios de alta simetría. Las soluciones obtenidas en esta forma

nos dan las características del campo gravitacional, · así como también las

29
 propiedades de la curvatura del espacio-tiempo y del tensor energía-momentum

A continuación se presentan algunas soluciones exactas de interés, cuyas

propiedades locales se han desarrollado anteriormente [5,6,7].

11.6.1 SOLUCIÓN DE SCHWARZSCHILD

Después que Einstein publicara su teoría de la relatividad general, en 1916 Karl

Schwarzschild encontró la solución exacta de la ecuación gravitacional de Einstein

para el vacío, que corresponde al campo gravitacional externo de una distribución

esféricamente simétrica de masa, siendo la métrica en coordenadas esféricas la

siguiente:

(2.105)

donde M es la masa del cuerpo, G la constante gravitacional , e la velocidad de la

luz y r el radio de curvatura.

En la ecuación anterior se observa la existencia de una singularidad para un radio

crítico, llamado radio gravitacional o de Schwarzschild:

2GM
rg =--z- (2.106)
c

30
 Cuando r ) ) rg, la métrica se reduce a un espacio-tiempo plano, donde la

influencia del campo gravitacional es nula.

La singularidad cuando r = rg es ficticia pues se puede eliminar por una

transformación de coordenadas adecuada, siendo la única singularidad real cuando

r =0.

Un cuerpo masivo cuyo radio es menor o igual que su radio gravitacional no puede

emitir luz y es lo que se llama agujero negro. La superficie definida por r = rg es

llamada "horizonte de eventos". Así, r =O se encuentra dentro del horizonte y no

puede verse desde el exterior.

11.6.2 SOLUCIÓN DE ROBERTSON-WALKER

Está basada en el principio cosmológico que dice que el universo es homogéneo e

isotrópico, con máxima simetría espacial. Este espacio de curvatura constante puede

depender del tiempo, siendo la métrica para este universo la siguiente:

(2.107)

donde R 2 (t) es el factor de escala sin dimensión, dependen sólo de t y K. La

curvatura espacial es determinada por el elemento de línea:

(2.108)

Entonces,

31
 K =O ; se tiene una variedad tridimensional plana.

K = +1 ; el espacio es de curvatura positiva constante, cerrado y que tiene un

volumen finito.

K = -1 ; es un espacio de curvatura negativa constante, abierto y que tiene un

volumen infinito.

El factor de escala R(t) opera sobre la parte espacial total, sin considerar la

dirección, inflando estos espacios de manera uniforme, los cuales se expanden o se

dR(t) .. .
contraen ya sea que - - es positivo o negativo.
dt

11.6.3 SOLUCIÓN DE REISSNER-NORDSTROM

Fue encontrada por Reissner (1916) y Nordstrom (1918), y es otra métrica con

simetría esférica que representa el campo gravitacional de un cuerpo de masa M y

carga Q:

Se observa que se reduce a la métrica de Schwarzschild para Q = O. Esta métrica

tiene una singularidad en r = O, y existen dos superficies concéntricas llamadas

horizontes definidas por:

(2.110)

32
 (2.111)

que son las soluciones de la ecuación g 00 = O. La singularidad está dentro del

horizonte interno r = r_. Los dos horizontes se juntan cuando GM 2 =Q y

desaparece si GM 2 < Q 2 • Para r =O queda desnuda, es decir no hay horizonte

que la cubra. Además, se observa que ninguna distribución de materia cargada con

un radio menor que r < ~ puede producir, en su exterior, un espacio-tiempo

2Me

cuya métrica tenga GM 2 < Q2 •

11.6.4 SOLUCIÓN DE KERR

Fue encontrada en 1963 por Roy Kerr, y es una métrica que describe un agujero

negro en rotación. Esta métrica en coordenadas de Boyer y Lindquist es la siguiente:

(2.112)

donde,

(2.113)

(2.114)

M representa la masa y a el momento angular por unidad de masa. Cuando a = O

la métrica se reduce a la de Schwarzschild.

33
 La singularidad de la métrica de Kerr no es un punto, sino un anillo de radio a en el

plano () = ;r . Si M > Jal, la métrica tiene dos horizontes definidos por:

2

(2.115)

(2.116)

Estos horizontes se juntan cuando M = lal y desaparecen para M < lal, dejando
la singularidad desnuda. Las propiedades globales de la métrica de Kerr son muy

similares a las de la métrica de Reissner-Nordstrom.

11.6.5 SOLUCIÓN DE KERR-NEWMAN

Esta solución es una generalización de la métrica de Kerr correspondiente a un

agujero negro con masa M, momento angular aM y carga Q. La métrica de Kerr-

Newman tiene la misma forma que la ecuación (2.112), con la siguiente relación:

(2.117)

Esta métrica se reduce a la de Kerr cuando Q = O, y a la Reissner-Nordstrom

cuando a = O. Sus dos horizontes están en:

(2.118)

(2.119)

34
CAPÍTULO III

FLUIDO EN RELATIVIDAD GENERAL

111.1 FLUIDO

En muchas situaciones en relatividad general, un fluido perfecto puede ser tomado

como el origen del campo gravitacional, como una primera aproximación [8]. En

general, un fluido es una clase especial de medio continuo. Un continuo es una

colección de partículas numerosas cuya dinámica, en forma individual, no puede

ser seguida. Sólo se puede describir. la colección en términos de cantidades

"promedio": número de partículas por unidad de volumen, densidad de energía,

densidad de momentum, presión, temperatura, etc. Por ejemplo, el comportamiento

de un "lago de agua", y el campo gravitacional generado, no depende de cada

molécula de agua en particular; éste depende sólo de las propiedades de la inmensa

colección de moléculas, en promedio. Sin embargo, estas propiedades pueden variar

de punto a punto en este lago: la presión es más grande en la parte inferior que en la

parte alta, y la temperatura puede variar también. La atmósfera, otro fluido, tiene

una densidad que varía con la posición. Así, las partículas en forma individuales no

importan, pero éstas deben ser suficientemente pequeñas para que sean

relativamente homogéneas. Con esto, la velocidad promedio y la energía cinética,

por ejemplo, son las mismas para todos lados en esta colección. Así, una colección

35
 es llamada "elementos". Este término es impreciso para usar en una colección

grande de partículas, que puede ser observado como un simple valor para las

cantidades como densidad, velocidad promedio y temperatura. Si esta colección no

existe (gas rarificado), entonces las aproximaciones continuas no funcionan para

describirlo.

En esta aproximación continua asignamos a cada elemento un valor de densidad,

temperatura, etc. Entonces, estos elementos son llamados "pequeños". Esta

aproximación se expresa matemáticamente asignando a cada punto un valor de

densidad, temperatura, etc. Así, un continuo es definido por vario campos, teniendo

valores en cada punto y en cada tiempo.

De este modo, la notación de continuo se enlaza para sólidos como también para

gases. Como un fluido es un continuo, a veces se consideran "flujos". Esta

definición no es muy precisa, y así la diferencia entre sólidos y fluidos no está muy

bien definida. Así, los sólidos deben fluir bajo altas presiones.

Cada colección de este continuo tiene su tensor energía-momentum, y su suma Tpv

codifica la densidad de energía para el fluido completo.

111.2 FLUIDOS EN GENERAL

Hasta ahora, para describir un fluido, se tiene que partir con la colección más simple

posible de partículas. Para generalizar esto a fluidos reales se debe tomar en cuenta

lo siguiente: (a) aparte del movimiento del fluido, cada partícula tiene la misma

velocidad aleatoria; y (b) pueden haber varias fuerzas entre partículas que

contribuyen en energías potenciales a la energía total.

36
 Para cada elemento del fluido, se debe considerar un marco en el cual el fluido es

como un todo (el momento total espacial es cero). Este es un marco de referencia

fisico, el cual es momentáneamente co-móvil; entonces, los elementos del fluido

pueden ser acelerados, un momento después un marco inercial diferente debe ser el

marco de referencia fisico. Además, dos elementos diferentes del fluido pueden

estar en movimiento relativo uno con respecto al otro, así que ellos no deben tener

los marcos de referencias fisicos. Por lo tanto, el marco de referencia fisico es

específico a un solo elemento del fluido, el cual es una función de posición y

tiempo.

Todas las cantidades escalares asociadas con un elemento del fluido en relatividad

(como la densidad, densidad de energía y temperatura) están definidas en el marco

de referencia físico.

111.3 FLUIDO PERFECTO

Un fluido perfecto en relatividad general está definido como un fluido que no tiene

viscosidad y no tiene conducción de calor en el marco de referencia fisico [5,8,9].

Esto es una generalización del "gas ideal" de la termodinámica ordinaria. Las dos

restricciones en esta definición simplifican enormemente el tensor energía-

momentum rw.
Este fluido, al no tener conducción de calor, implica inmediatamente que en TJLY,

dentro del marco de referencia físico, T 0; =T;o = O ( i = 1,2,3) . La energía puede

circular sólo si las partículas circulan.

37
 Se sabe que la viscosidad es una fuerza paralela en la interfase entre las partículas.

Considerando las coordenadas espaciales, el tensor TiJ debe ser cero cuando i *j
( i = 1,2,3 y j = 1,2,3) . El tensor TiJ debe ser una matriz diagonal simétrica de 3x3.

Además, estos elementos son diagonales en todo marco de referencia físico. Al ser

la viscosidad cero, el número de funciones o variables se reduce a una, siendo ésta

la presión. Para esta matriz 3x3, entonces debe cumplir:

(3.1)

III.4 HIDRODINÁMICA CLÁSICA

Para un fluido perfecto, donde se desprecia la viscosidad y conducción de calor, las

ecuaciones fundamentales que describen su movimiento son:

Conservación de la masa:

ap + v.(pv)= o (3.2)
at

Ecuación de Euler:

pDv - + (v.V )v =-Vp+F

- = pav (3.3)
dt at ex/

donde p =densidad, v = velocidad, p = presión, Fext =fuerza externa aplicada el

fluido.

38
III.5 TENSOR HIDRODINÁMICO RELATIVISTA

Describe la relación entre la materia y el campo, y viene dado por el tensor

simétrico energía-momentum T 11 v [8]. Para el caso de fluido perfecto, donde se

desprecia la transferencia de calor, viscosidad y la anisotropía debido a la presión, el

tensor energía-momentum puede ser escrito como:

(3.4)

donde U 11 es la cuadrivelocidad

1
U 11 = (1, v) (3.5)
-JGv2
p y p son, respectivamente, la densidad de energía y la presión isotrópica (con la

geometrización de unidades e =1 y G = 1), L~s componentes del tensor son:

T 00 ~densidad de energía

T 0 ; ~ densidad de flujo de energía y momentum (3.6)

TiJ ~ tensor de tensiones o de esfuerzo

donde, i = 1,2,3 y j = 1,2,3 . Se verifica también la simetría del tensor

(3.7)

39
 En el sistema de referencia ca-móvil (v =O), el tensor energía-momentum para el

fluido es de la forma [10,11]:

o o o
p

Tll =
o p o o
(3.8)
V
o o p o
o o o p
Entonces, para un fluido anisotrópico su tensor energía-momentum es el siguiente:

o o o
p

Tll =
o p¡ o o
(3.9)
V
o o P2 o
o o o P3
donde, p 1 '* p '* p
2 3, que indica que las presiones son diferentes.

111.6 LEY DE CONSERVACIÓN DE LA ENERGÍA MOMENTUM

La ley de conservación de la energía-momentum en ausencia de fuerzas exteriores

toma una forma tensorial simple en relatividad especial [5,8]. Consideremos un

pequeño cubo de materia de lado a (figura 1).

J.-----
/ "(0,0,0)
/ "

Figura 1

40
Analizando el balance energético, las energías entrantes dentro del cubo para sus

diferentes caras está dado en forma tensorial por:

Respecto a la cara trasera :

Respecto a la cara delantera: - a 2 T 0 x (x =a)

Respecto a la cara derecha :

(3.10)

Respecto a la cara izquierda: - a 2 T 0Y (y= a)

Respecto a la cara inferior : a 2 T9z. (z = O)

Respecto a la cara superior : - a 2T 0 z (z =a)

El signo menos proviene del hecho de que los flujos son opuestos a los ejes

positivos.

La energía del cubo por unidad de tiempo es igual a:

41
Dividiendo entre a 3 y tomando el límite cuando a ~ O, se obtiene:

(3.12)

donde, i = 1,2,3 . En consecuencia:

Tov
,V
=0 (3.13)

Aplicando el mismo razonamiento con las componentes espaciales del cuadri-

momentum, se verifica

Tiv
,V
=0 (3.14)

Así, la conservación de la energía-momentum se escribe en relatividad especial por

cuatro ecuaciones que cumplen:

rw =o
,V
(3.15)

Esta última ecuación es válida en relatividad especial pero, como está escrita en

forma covariante, es válida también en relatividad general debido al principio de

covariancia general [12, 13].

Finalmente, la conservación de la energía-momentum de forma general es escrita

como la derivada covariante del tensor contravariante TJLv de la siguiente manera:

rw =O
;v
(3.16)

42
CAPÍTULO IV

ECUACIONES DE CAMPO

IV.l ELEMENTO DE LÍNEA

Consideremos el elemento de línea más general para representar un campo

gravitacional. Recordando la ecuación (2.32) se tomará la expresión:

(4.1)

Transformando las coordenadas espaciales a coordenadas esféricas

(x,y,z) ~ (r,B,rjJ), la expresión más general para ds 2 con simetría esférica está

dada por [4,5,6]:

con su tensor métrico

A(r) o o o
-C(r) -B(r) o o
= (4.3)
g¡JV
o o -D(r) o
o o o - D(r )sen 2 B
donde,

g 11 =A(r), grr =-B(r), g 88 =-Ji...r), g?? =-D(r)sen 2 B,gr =-C(r)

1 (4.4)

43
Como la elección del sistema de coordenadas es arbitraria en relatividad general, se

plantea una transformación de coordenadas, que no altere el elemento de línea ds 2 •

Para ello, se define una nueva coordenada de tiempo, de la forma:

t· = t + 'P(r) (4.5)

donde 'P(r) es una función arbitraria de r.

En (4.2), para eliminar el elemento g ,1 que está fuera de la diagonal de la

representación matricial del tensor métrico g pv , se considerará:

dlJ! _ c(r)
(4.6)
dr -- A(r)

A partir de (4.5) y (4.6)

dt· = dt + d'P(r) ~ dt = dt· - dlJ!(r) (4.7)

dlfl(r) =- c(r) dr (4.8)

2A(r)

Remplazando (4.7) y (4.8) en (4.2), se tiene:

44
 - c(r )df'<ir + c(r{- 2~~) dr ]dr- D(r Xde' + sen 1A:W)
2

Haciendo que:

C 2 (r)
E(r)=-+B(r) (4.10)
4A(r)

Reemplazando (4.1 O) en (4.9) se tiene:

Considerando que:

r· = F(r).r (4.12)

y diferenciando esta última expresión,

45
 dr· = F(r )dr + r dF(r) dr = [F(r) + r dF(r )]dr
dr dr

(4.13)

Reemplazando (4.13) en (4.11) se obtiene:

Finalmente, considerando que:

(4.15) en (4.14)

Así, el elemento de línea en su forma estándar es el siguiente:

46
Las funciones G(r), H(r) y K(r), se pueden expresar de forma exponencial de la
forma:

G(r) = ev(r) (4.18)

H(r) = eJ.(r) (4.19)

K(r) = e¡.¡(r) (4.20)

Para que estas funciones sean válidas, las funciones deben tender a cero cuando r

tiende al infinito.

Con estas definiciones el elemento línea toma la forma:

(4.21)

con su tensor métrico

ev o o o
-e;. o o
gf.JV = o (4.22)
o -ef.i o
o o o - e 11 sen
2
(}

De acuerdo a los trabajos de Schwarzschild [15], la función e 11 <rJ es proporcional a

2
r • Considerando esta función como:

(4.23)

47
 Finalmente, el elemento de línea llamada solución de Schwarzschild que se

trabajará de aquí en adelante es:

(4.24)

IV.2 ECUACIONES DE CAMPO DE EINSTEIN

Consideremos un fluido anisotrópico con simétrica esférica, que carece de

viscosidad y conducción de calor. Este fluido es representado mediante el tensor

energía-momentum r:. De acuerdo al principio de covarianza general, se elegirá

un sistema ca-móvil en coordenadas esféricas. Entonces, el tensor r: tendrá una

representación matricial [14], cuyos elementos no nulos son diagonales definidos

de la forma:

p o o o
o - Pr o o
Tll = (4.25)
V
o o -p(J o
o o o -p,p

donde: p = es la densidad, p, =presión respecto al radio r, PB = presión

respecto al ángulo colatitud () y p tP =presión respecto al ángulo azimutal.

De acuerdo a la métrica de Schwarzschild (4.24), la métrica de este elemento línea

toma la representación matricial:

48
 ev o o o
-e;¡ o o
gpv = (4.26)
o o -r2 o
o o o - r 2sen 20

Como g pv es un tensor covariante y simétrico de acuerdo a la condición (2.27),

entonces su inversa g 11 v es:

e-v o o o
-e -A,
o o
gpv = (4.27)
o o -r -2 o
o o o -r- 2sen- 20

De acuerdo a la teoría de la relatividad general, la materia está relacionada

directamente con la curvatura del espacio-tiempo. Las ecuaciones de Einstein están

dadas por la ecuación (2.1 00), para lo cual es necesario calcular previamente los

símbolos de Christoffel. Empleando (2.38) y (2.39) se verifica que:

(4.28)

De la expresión (4.26), el tensor métrico g pv tiene las componentes no nulas:

g 11 =e V , g,, =-e A , g 00 =-r 2 , g~~ =-r 2 sen 2() (4.29)

Igualmente, las componentes no nulas de su inversa gpv son:

(4.30)

49
Calculo de los símbolos de Christoffel.

Con (4.29) ,(4.30) y (4.28) se calculan a continuación los símbolos de Christoffel no

nulos.

_!_ rr _ _!_(-e -A )_g__ (-eA)

1 rrr _
- 2g grr,r - 2 ar

(4.31)

r eer =-re -A (4.32)

r~~r = -21 g rr (- g
~~.r
) = -21 (-e -A {-ar
8 (- r 2 sen 2(}))
-

r ~~r = -rsen 2/lue -A (4.33)

(4.34)

50
 (4.35)

o
r~~ 1 oo (g~~(} ) =-
=-g 1( -21 )( - 2 )~
-r 2 sen(} a(
2 · 2 r aB

r:~ = -senBcosB (4.36)

~ - - ~~ (g ~~.r ) -_ -1 ( -
r,~ _ 1 g 2
1 2
)( - a(- r 2 sen 2 (} ))
2 2 r sen (} ar

1
rr~-
~ -r~-­
~- (4.37)
r

~ _ 1
r(}~ --g ~~(g~~.(} )-- -
1 (- 2 1 2 )( - a (-r sen(}
2 2 ))
2 2 r sen (} aB

(4.38)

rl
rt
1
= 2g
11 (
gtt,r
) 1r
= 2\e
-V f\ ar
a \er V))

(4.39)

51
Cálculo de las componentes del tensor de Ricci.

A partir de (2.78) y (2.79) se calculan a continuación las componentes no nulas del

tensor de Ricci.

aro ar¡>l ar'

- ---
R rr - rO
- - -r¡>l- - -rt+ rrrO rrr¡>l rrr/ roro r¡>lr¡>l r'r'
rr rO+ rr r¡>l + rr r l - rO rO- r¡>l r¡>l- rt rt
8r 8r 8r

(4.40)

- oo
arr ar¡>l
o¡>~ rr00 r'rr + rr00 r¡>lr¡>l + r 000r'r t - rr00 ro(;Ir- r¡>lO¡>l r¡>lO¡>l
R00 ------+
8r ae

8(-re-J..) 1( -;.)11 , _;,.1 _;,.V' _;,.1 2B 8ctg8

= --re Y'L- -re --re -+re --ctg - - -
8r 2 r 2 r ae

(4.41)

52
 = - a (rsen 211ue -A.) - - a (sen (} cos o) -
~ rsen 2(} - k -sen 2n~-A.
e-A.- ~ +sen 211~-A.
~
& ao 2

+ senBcosBctgB- _!_ rsen 2 Be-2 v'

2

(4.42)

art; ['r['r
R 11 ----+ 11 rr+
['r['IJ rrr~P ['r['t
11 riJ+ 11 r¡p- 11 Ir
ar

(4.43)

Cálculo del escalar de Ricci (curvatura escalar).

A partir de (2.80) se obtiene lo siguiente:

53
 =ev ev-A --V"
[ ( 2
1
1 +-A'V' 1 ( )2 -V'
--V' 1 +-A'V'
- )] -e-A [ --V" 1 1 ( )2 +-
--V' A:]
4 4 r 2 4 4 r

(4.44)

Cálculo de las componentes del tensor de Einstein.

A partir de (2.81) se calculan las componentes no nulas del tensor de Einstein.

1
Grr = Rrr- 2grrR

1 +-í!.·v·
=--v" 1 --1(v· )2 +-+-eA
í!.· 1 [ e-Av, --e-AA:v·
1 1
+-e-A (v· )2
2 4 4 r 2 2 2

54
 1 V' 1 ,t
G =-+---e (4.45)
rr r2 r r2

1
Gee = Ree-- geeR
2

=-e" 1(v· -íl: )~ +1+-r

-[1+-r 1 2 [- 1_"íi:v· +-e"
e "v" --e 1_()2
v·
2 2 2 2

G 1 2e-" [ V"+-1 (v· )2 + V' -

= -r íl: --v·íl:
1 ] (4.46)
ee 2 2 r 2

(4.47)

55
=-ev-A[ --V"+-A:V'--V'
1 1 1 ( )2 - - --eve V']
1 [ -A V"--e
1 -A _.1:v·+-e
1 -A (V' )2
2 4 4 r 2 2 2

0 11 =e v [ -e -A ( -
r2
A: +-
1-- 12 ]
r
J r
(4.48)

De (4.25) y (4.26) se obtiene una representación matricial, cuyos elementos no

nulos son diagonales y definidos de la forma:

pev o o o
o p,e
A
o o
= (4.49)
T/lV
o o Por
2
o
o o o pi/Jr 2 sen 2 B

Entonces, las componentes no nulas del tensor Tw son:

(4.50)

(4.51)

56
 (4.52)

T;, = pev (4.53)

Reemplazando las ecuaciones (4.45)-(4.48) y (4.50)-(4.53) en la ecuación de campo

(2.1 00), se obtiene:

1 V' 1 A A
- + - - - e = 8.np e (4.54)
r2 r r2 r

(4.55)

(4.56)

e v[ -e -A(-12- -A,') +-12] = 8;rpe V (4.57)

r r r

Se observa que la ecuaciones (4.55) y (4.56) son iguales, entonces:

Po =pi/J (4.58)

57
 Al ser los valores de la presión p 0 y p ¡p iguales, y definiendo como p .L la presión

tangente, entonces se cumple:

P.1=Pe=P,p (4.59)

y denotando

p=pr (4.60)

A partir de (4.59) y (4.60) las 4 ecuaciones de campo se reducen a 3, siendo las

siguientes:

1 V' 1 A A
-2+ - - -2e =8.npe (4.61)
r r r

(4.62)

V [

r r
J
A: +-
e -e -A ( -12 - - 1
r2
l = 8;rpe V (4.63)

IV.3 ECUACIÓN DE EQUILIBRIO HIDROST ÁTICO

A partir de la ecuación de conservación de energía-momentum (3.16), se obtiene la

ecuación de equilibrio hidrostático. Esta ecuación de conservación para el r: es de

la forma [8,9]:

(4.64)

58
Empleando la ecuación (4.25) y (4.61) se calculan a continuación las componentes

mixtas no nulas de la derivada covariante del tensor T :,li :

o
Tr!J = PJ.- p (4.65)
, r

~ -
Tr;~ - -
r~ T~
r~ ~ + r~r~ yrr

T~ = PJ.- p (4.66)
r;~ r

1 1
=--v·p+-v·(-p)
2 2

59
 1 1
T,.,
'
= --(p + p)v·
2
(4.67)

Sumando (4.65), (4.66) y (4.67), la ecuación (4.64) se reduce a:

r:'
11
=- dp +2 (p.L- p) _ _!_(p+ p)v· =o
dr r 2

(4.68)

60
CAPÍTULO V

SOLUCIÓN DE LAS ECUACIONES DE CAMPO PARA UNA

ESFERA ANISOTRÓPICA

V.l CÁLCULO DE LOS COEFICIENTES MÉTRICOS

De la ecuación (4.63), desarrollando para r E [0, r], se obtiene:

87rp=-e-A - ..1:) +-21

1--
(r2 r r

e--1 = 1- 2g 47rpr2dr (5.1)

r

Se define la masa de un fluido anisotrópico de simetría esférica como [ 16, 17]:

(5.2)

61
Reemplazando (5.2) en (5.1):

(5.3)

El valor del termino -e" es el coeficiente métrico correspondiente a g ,, (con la

geometrización de unidades e = 1 y G = 1). Esto verifica lo definido en (2.1 05).

Ahora, reemplazando (4.68) y (5.3) en (4.61), se reduce a:

87q7 =e--< - 1
( r2
V') - -12
+-
r r

8np = (1-2 r
m(r) )[__!__ + ~[~{p _L - p)- p· )l]- __!__2
r2 r r J r

Se define a continuación una función adimensional g{r) y la llamaremos función

generatriz, definida por:

2m(r)
1---
g = r (5.5)
(r) 8npr 2 + 1

62
Reescribiendo la ecuación (5.4) en términos de la función g(r), se obtiene:

2p·
(5.6)
r(p+ p)

Como la masa está en función del radio r, entonces:

dm(r) 2
dr = m(r) = 4trpr (5.7)

Igualmente, de (5.5) se calcula la derivada de g(r).

(5.8)

Reemplazando (5.7) en (5.8) se obtiene la densidad de energía p del fluido:

63
 (5.9)

Reemplazando la ecuación (5.9) en (5.6) se obtiene:

8np' + (1- g(r) + llJ(r)XI- 3g(r)- g(r)r )np

rg(r)(I + g(r) -m(r))

(5.10)

donde se introdujo una función anisotrópica m(r) definida por:

4(p- p_~_)
ll)(r) = (p+ p) g(r) (5.11)

64
La ecuación (5.10) es una ecuación diferencial de p respecto a r. Los coeficientes

de p y el término independiente son funciones continuas. Esta ecuación se puede

escribir en forma abreviada como:

(5.12)

donde,

(5.13)

(5.14)

Para resolver la ecuación diferencial (5.12) se aplica el método de variación de la

constante.

Si C{r) no está presente, entonces la ecuación (5.12) es una ecuación homogénea de

la forma:

(5.15)

La solución de la ecuación (5.15) es:

(5.16)

65
donde p 1 es la constante de integración no nula de la ecuación diferencial

homogénea.

Ahora, considerando que p 1 = PJ(r) es función de la variable r, la ecuación (5.16)

toma la forma:

(5.17)

Esta última ecuación es una ecuación no-homogénea. Derivando (4.17) se tiene:

(5.18)

Reemplazando (5.18) en (5.15) se obtiene:

, _ e( )ef B(r)dr (5.19)

P1(r)- r

(5.20)

donde, p 0 es una constante de integración.

Reemplazando (5.20) en (5.17) se obtiene:

87rp -_ ( Po + Je (r)e
J
B(r)drd )
r e
-J B(r)dr .
(5.21)

66
Para calcular la densidad de energía p, se procede a dar forma a la ecuación (5.8).

(5.22)

Reemplazado (5.22) en (5.9) se obtiene:

De la ecuación (5.11) se obtiene la presión tangencial:

(5.24)

Finalmente, se calculan los coeficientes métricos en función de g(r) y m(r)·

Reemplazando (5.5) en (5.3) se obtiene:

(5.25)

Ahora, reemplazando (5.25) en (4.54):

67
 8.np =e-A. -1
(r
2 +- 1
V' - -
2
J
r r

V' = (1- g(r)) (5.26)

rg(r)

Resolviendo (5.26) por integración:

Jdv = r(-1- -
Jlrg(r) r r· +
_!_ t,r cte

dr
Irg(r)
v = ---Lnr+Lnv 0 (5.27)

donde, cte = Ln v 0 y v 0 > O• Entonces, para determinar el coeficiente métrico e v :

eV = .:i__ e
J-
2
1
-dr
rg(r)
(5.28)
r

con A= v 2 •

68
CAPÍTULO VI

VERIFICACIÓN Y ANÁLISIS DE LA SOLUCIÓN

OBTENIDA

VI.l SOLUCIÓN GLOBAL

Reemplazando (5.25) y (5.28) en (4.25), el elemento de línea ds 2 con los

coeficientes métricos toma la forma global:

Se observa, que ds 2 no depende de la presión tangencial p .L y la función

anisotrópica, teniendo un carácter radial.

VI.2 SOLUCIÓN LOCAL

La función g(r) debe satisfacer algunos requerimientos generales. Si asociamos una

solución, ésta debe tener significado físico [18]. Asumiendo una presión no

divergente en el origen, entonces para un espacio plano local debe cumplir con las

siguientes condiciones:

69
 r~1

P(r~o)-:¡:. O (6.2)

Con las condiciones (6.2) y la ecuación (5.25) se tiene que:

(6.3)

Se observa que g(r) ~ 1 para r ~O . En realidad la simple elección de g(r) = 1 nos

da el espacio-tiempo plano.

Para el vacío, es decir sin la presencia de presión radial p = p 0 = O, la ecuación

(5.21) se reduce a:

g(r) 1
g¡,)+-=-
r r

efe
g(r) =I-- (6.4)
r

donde efe= 2m, y reemplazando esto en (6.4) se obtiene:

gs(r)= 1-2m (6.5)

r

donde gs(r) corresponde a la solución Schwarzchild para el vacío. Se observa que

esta solución no depende de la función anisotrópica m(r)·

70
 o
-1

-2

-3

Cl -4

-5

-6

-7

-8

-9~--~--~--~----L---~--~--~~--~--~--~
o 10 20 30 40 50 60 70 80 90 100

Gráfica l. Dependencia de la presión respecto al radio de la ecuación (6.4). Se usó el programa Matlab, con
valores de radio de 1 hasta 100 en intervalos de 1 de 0.5, donde se consideró Cte =l.

71
 Ejemplos.

Ahora, para este espacio-tiempo plano, se procede a determinar la métrica de

Schwarzchild g(r) sobre la superficie r = r0 , con la continuidad de g(r)sobre r0 •

(6.6)

• Primer ejemplo.

Considerando un caso trivial definido por:

g(r)=1+ar 2 (6.7)

m(r) = -ar 2 (6.8)

donde a es una constante. Esta selección conduce a un modelo físicamente

razonable, de manera que la función g (r) ~ 1y r ~ O.

Reemplazando (6.7) en (6.6):

g ( r0 ) = 1+ ar02 = 1-2m
-
ro

2m
a=-- (6.9)
r.3
o

Como se usa una métrica de ~chwarzchild en el vacío, la presión radial desaparece.

72
 p(r)= O (6.10)

De la ecuación de la densidad de energía (5.23), se obtiene:

1
8n-p = (1- g(r ))- - g' (r)
r2 r

3a
p=-- (6.11)
8n-

Reemplazando (6.9) en (6.11)

m
p=-- (6.12)
4
-;rr.
3
3 o

La presión tangencial p .L , la cual tiene esta esfera homogénea, es calculada a partir

de la ecuación (5.24). Además, considerando la condición (6.10) se obtiene:

(6.13)

73
Reemplazando (6.7), (6.10) y (6.12) en (5.25), se obtiene el coeficiente métrico

,t 1
-e = - - -
g(r)

(6.14)

Finalmente, reemplazando (6.7) y (6.9) en (5.28) se obtiene el coeficiente métrico

v A2 f r(!+ar1
2)
d
r
e =-e
r

=~e
2 Ln[-r-P
{I+ar
2
1]

74
 (6.15)

La ecuación (6.15) implica que en el límite exterior r0 > 2m .

Su elemento de línea será de la forma:

• Segundo ejemplo.

Considerando el caso:

(6.17)

m(r) = const (6.18)

y reemplazando (6.17) y (6.18) en (5.13) y (5.14), respectivamente, se obtiene:

B = (1- 3b XI-
b +m)
(6.19)
rb(1 +b- m)

75
 e = - (I- b +mXI- b) (6.20)
r 3b(I +b- m)

Reemplazado (6.I9) y (6.20) en (5.2I):

8JqJ = po+ f- (I- b + mXI- b) e J(!-3b)(!-b+m)dr

rb(l+b-OJ) dr e
J -J(!-3b)(!-b+m)dr
rb(!+b-m)
( r3 b(I + b -m)

(l-3b )(1-b+m) L ]
I- b) e
8JqJ =
(
Po - (
b(I+b-m)
J b(!+b-m)

r3
nr
dr e
(!-3b )(1-b+m) Lnr
b(!+b-m)

(I-3b)(l-b+m) I [ (b-IXI-b+m) ]
8JqJ =Por b(!+b-m) +- --:-'----'--'----~- (6.21)
r2 b2 - 6b- bm +m+ I

Estos coeficientes no dependientes de r los denotamos como:

D = (I- 3b XI- b +m) (6.22)

b(I +b -m)

V= (b -I X1- b +m) (6.23)

2
b - 6b- bm +m+ 1

Con (4.50) y (6.23), la ecuación (6.21) se escribe:

-D V
8JqJ =Por +2 (6.24)
r

76
 200

120

a. 100

80

60

40

20

oo 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8

r

Gráfica 2. Dependencia de la presión respecto al radio de la ecuación (6.24). Se usó el programa

Matlab, con valores de D que va de O hasta 2 en intervalos de 0.5, donde se considero C = 1, V=l,
8x=l.

77
Ahora, reemplazando (6.24) en (5.23):

8;rp = (1-b-bV) 21 + p 0 b(D-3)r - D (6.25)

r

Para determinar la presión tangencial se reemplaza (6.25) en (5.24).

m(p+p)
pj_=- +p
4g

8.npJ. = - OJ- ( (1 - b -2 b V) + p 0b(D-3 )r -D + p 0 r -D +-

V)
2
4b r r

8.npj_ = (v- ~(1 +

4b
vX1-b ))-2 + (1- ~(1 + bD- 3b )) Po
r
1
4b rD
(6.26)

Reemplazando (6.17), (6.24) en (5.25), se obtiene el coeficiente métrico -eA:

-eA=- [
2
( Por -D +-2 V]]
b 8nr Sí'Z' r +1

78
 (6.27)

Ahora, remplazando (6.17) en (5.28), se obtiene el coeficiente métrico ev:

(6.28)

Finalmente, su elemento de línea será de la forma:

l-b
ds 2
=A 2
r ¡; dt 2 -(p br
0
2
-D +bV +br 2 }ir 2 -r 2 (d() 2 +sen 2 ()drp2) (6.29)

79
 1

0.9

0.8

o(J 0.7
-E 0.6
·e
• Q)

.....
Q)
e:
Q)
'(3 0.5
'io
(J
0.4

0.3

0.2

0.1
o 50 70 80 90 100
r

Gráfica 3. Dependencia del coeficiente métrico respecto al radio de la ecuación (6.28). Se usó el programa

Matlab, con valores del radio que han variado de O a 100 en intervalos de 1, donde se consideró A= l.

80
VI.3 SOLUCIÓN DE LA MÉTRICA INHOMOGENEA DE FRIEDMANN

Considerado la métrica [4,5]:

ds 2 = dt 2 - B(r,t )ldr 2 + r 2 (d0 2 + sen 2 0d~ 2 )J (6.30)

donde B = variable de expansión o contracción.

Para este ejemplo se considera p =p j_ , entonces la función anisotrópica (5.11)

será nula:

(6.31)

El tensor de energía- momentum para procesos irreversibles es [8,9]:

(6.32)

donde: p = es la densidad, p =presión y q 11 = vector conducción del calor.

De acuerdo a la métrica de inhomogénea de Friedmann (6.30), este elemento de

línea toma la representación matricial:

1 o o o
-Bz o o
g/lV =o o -Bzrz o
(6.33)

o o o - B 2 r 2 sen 20

81
Como g 11 v es un tensor covariante y simétrico de acuerdo a la condición (2.27),

entonces su inversa g 11 v es:

1 o o o
1
B2
o o
1 (6.34)
o ---
B2r2
o
1
o o o 2 2
B r sen 2 0

De la expresión (6.33), el tensor métrico g 11v tiene las componentes no nulas:

gil= 1 ' grr =- B 2 ' goo =- B2 r 2 ' grPrP =-B2 r 2sen 20 (6.35)

Igualmente, las componentes no nulas de su inversa g 11 v son:

g il =1 grr = -1- (6.36)

' B2 '

Calculo de los símbolos de Christoffel.

Con (6.35) ,(6.36) y (4.28) se calculan los símbolos de Christoffel no nulos.

rl _ _!_ _!_(1{-~(-B2)~
11 _
rr - 2 g g rr,l - 2 \ at )

(6.37)

(6.38)

82
 (6.39)

2
,1 B 2 )(~(-B
r'Ir -_ _!_2 g ''(g rr )-_!_(--- at
2 )~
')
1
r rr B·
rlr= ri=B (6.40)

rr,, _ B2 J(~(-
- _!_ g rr ( g,,,, ) -_ _!_(- _1 at B2 )1l)
2 2

r'rr =~
B (6.41)

2 B2 J(~
r 88r -_ _!_ g rr (- g88,r ) -_ _!_ (- _1
2 ar (- B 2 r 2 )~
)

(6.42)

r ,;,;r --
- sen2() (B· r 2 + Br ) (6.43)
B

83
 (6.44)

(6.45)

r:~ = -sen()cos() (6.46)

~
1,~ _
- -1 g ~~ (g #,t ) -_ -1 ( - 1 a(
)( - - B 2 r 2 sen 2 () ))
2 2 B 2 r 2 sen 2 () at
B·
~ -r~ - -
r~~- (6.47)
~- B

~ _-
1,~ - 1g ~~ (g~~.r ) -
_-1( -
2 2
1
2
)(- a( - B 2 r 2 sen 2 () ))
2 2 B r sen () ar

(6.48)

84
 rP _ 1 ,p,p (g,p,pe ) -
re,p--g _-1 (- 1 2 )~
)(-a-(B 2 r 2 sen(}
2 · 2 B 2 r 2 sen 2 (} ae

(6.49)

Cálculo de las componentes del tensor de Ricci.

A partir de (2. 78) y (2. 79) se calculan las componentes no nulas del tensor de Ricci.

~- !( ~ J- !( ~ J- !(~ J-( ~ g J--(~ ~ J-(~ x~ J l

B··
R =-3- (6.50)
tt B

arl arere arrPr,P rl re rl r,p r,pr,p rrre rrr,P re re r,p r,p

R rr -- + Btrr +----¡;;------¡;;- + rr te + rr t,P + t,P t,P + rr re + rr r,P - re re - r,P r,P

2 B" {D. \2 B·
R =2 ( B· ) +BB .. -2-+2~-2-
2
(6.51)
rr B B Br

85
 rB·
--ctgO
B

R = 2rz (B· )z + Brz B-· B·'rz -3 rB· (6.52)

-~ B B

(6.53)

(6.54)

86
Cálculo del escalar de Ricci (curvatura escalar).

A partir de (2.80) se obtiene lo siguiente:

R =R 11 g 11 +R rr g'' +R 88 g 88 +R ¡p¡p g¡p¡p

B'· ( - - 2B
=-3-+
B B2
1J[ ( .)2 +BB . -2-+2~-2-
B"
B
{B}
B
B·]
Br 2

(6.55)

Cálculo de las componentes del tensor de Einstein.

A partir de (2.81) se calculan las componentes no nulas del tensor de Einstein.

1
G,, = R,, - g,,R
2

(6.56)

(
=2B
.)2 +BB . - 2B"- + 2(B·-Y B· 1 ( { B'· (B· )2
- - 2 - - - -B 2 -6--6~+4-
B"
B B2 Br 2 B B2 B3

87
 (6.57)

1
Gee =Roe-- geeR
2

(6.58)

(6.59)

1
G,, =R11 --g 11 R
2

88
 =-3--- B-· {B·
B·· 1 [ -6--6--+4--2
2
B"
3 4
Y
(B· + 8
B·-] Y
B 2 B B B B B3 r

(6.60)

De (6.32) y (6.33) se tiene una representación matricial, cuyos elementos no nulos

son diagonales definidos de la forma:

p -B2q o o
-B2q B2p o o (6.61)
TJlV =
o o B2r2p o
o o o sen 2BB 2r 2p

Entonces, las componentes no nulas del tensor Tpv son:

I;, = p (6.62)

(6.63)

(6.64)

(6.65)

(6.66)

89
Remplazando las ecuaciones (6.56)-(6.60) y (6.62)-(6.66) en la ecuación de campo

(2.1 00), y considerando K = 8nG (con la geometrización de e = 1 ):

B·· 2B·B·- V B2 (6.67)

2 - - ---.n..q
B B2

(6.68)

- (B . )2 r 2 - 2BB·· r 2 +B"r2
- - +B·r Y
- - (B· 2 r2 --.n..p
v B2 2
r (6.69)
B B B

(6.70)

(6.71)

Se observa que la ecuaciones (6.69) y (6.70) son iguales, entonces las 5 ecuaciones

de campo se reducen a 4. Utilizando las ecuaciones de campo (6.68) y (6.69) se

obtiene:

(6.72)

90
Para esta última ecuación diferencial, el B tiene la forma particular de:

(6.73)

donde C y r son constantes.

A partir de (6.44) y las ecuaciones de campo (6.67), (6.68) y (6.71), se obtiene lo

siguiente:

e
1 [ 4r~er
q =K 'J (6.74)

(6.75)

(6.76)

Si consideramos !_ (( 1 , entonces:
T

(6.77)

Las ecuaciones (6.74), (6.75) y (6.76) toman las formas dependientes del radio:

q=4-r
e (6.78)
K-r

91
 (6.79)

~(
2

p = __!_[12e +
K -r
e J
r +e
2 ] (6.80)

Si consideramos !__)) 1 , entonces:

'r

(6.81)

Las ecuaciones (6.74), (6.75) y (6.76) toman las formas dependientes de la variable

temporal:

e ']
1 [ 4r--;:-er
q =K (6.82)

1[
p=- -4eer!_ -3-
K
1]
-r2
(6.83)

1[12eer 3]
p= K
1

+~ (6.84)

92
Se sabe que el parámetro de expansión e es:

B·
e=3- (6.85)
B

De (6.73) en (6.85) Se obtiene:

(6.86)

Se observa de la ecuación (6.85) que e< O, entonces se tiende al colapso.

93
CAPÍTULO VII

CONCLUSIONES

l. A partir de las funciones anisotrópica y generatriz se resolvieron las ecuaciones de

Einstein para un fluido anisotrópico con simetría esférica y, de esa forma, se pudo ·

expresar la densidad y la presión respecto a estas dos funciones.

2. Se verificó analíticamente que la presión p 0 y P,p son numéricamente iguales,

reduciendo de 3 a 2 componentes de la presión para un fluido anisotrópico

(ecuación (4.58)).

3. Se obtuvo también una forma general para la solución de Schwarzschild, cuyos

coeficientes métricos no dependen de la función anisotrópica m(r) ni de la presión

tangencial p 0 = p J., dada por la ecuación (6.1):

4. Al analizar localmente solución de Schwarzschild se encontró que la función

anisotrópica g(r) es la métrica de Schwarzschild (ecuación (2.1 05) ), expresada en

unidades geométricas por la ecuación (6.5):

94
 gs(r)= 1-2m
r

5. Para dos ejemplos propuestos de g(r) y co(r) se encontró que la solución de

Schwarzschild no posee singularidad cuando r =O (ecuaciones (6.16) y (6.29)).

6. Se encontró en el primer ejemplo de la solución local (ecuación (6.15)) que en el

límite exterior, es decir fuera de la distribución, el radio gravitacional r0 debe

cumplir r 0 > 2M, muy similar con la solución de Schwarzschild para el vacío que

corresponde a un campo gravitatorio producido por un cuerpo de masa m con

simetría esférica.

7. Con las condiciones de la función generatriz (6.17) y función anisotrópica (6.18),

para que la presión se anule, es necesario que el exponente -D sea menor que la

unidad (incluido el cero). Con ello se garantiza que la presión radial se anule para

un radio grande.

8. Para la ecuación (6.27), para que el coeficiente métrico o potencial decaiga hasta

tender al cero, el exponente 2-D y la variable b deben ser menor que la unidad

(incluido el cero). Con ello se garantiza que para un radio grande se tiene un

espacio-tiempo plano.

95
9. En la gráfica para la ecuación (6.28), para que el coeficiente métrico o potencial

decaiga hasta tender al cero, el exponente t~h tiene que disminuir de la unidad

(incluyendo el cero). Con ello se garantiza que para un radio grande pero de

exponente pequeño se tiene un espacio-tiempo plano.

10. A partir de la métrica inhomogénea de Friedmann (6.30) para el fluido

anisotrópico, se analizaron las variables termodinámicas de densidad, presión y

conducción de calor, encontrándose dependencia exponencial del radio y tiempo.

Sin importar la dependencia de radio y/o tiempo para el fluido anisotrópico, éste

tendría un colapso gravitacional.

11. Al considerar la métrica inhomogénea de Friedmann, se anuló la función

anisotrópica (6.31) sólo trabajando con presiones con simetría en las 3 direcciones

espaciales. Al hacer eso no se pierde el carácter general, debido a que las variables

termodinámicas dependen de la variable radial y temporal.

t
12. Para -(( 1 se obtuvo que las variables densidad (6.80) y presión (6.79) tienen
r

carácter radial y, cuando el radio tienda al finito, en los extremos se tendría un

espacio-tiempo plano. En cambio, para !_)) 1 las variables densidad (6.84) y

r

presión (6.83) tienen carácter temporal.

96
13. Al hacer la consideración de !_ (( 1 ó !_)) 1, se observa que en las ecuaciones
'r 'r

(6.74) y (6.82) la variable conducción de calor q se mantiene invariable,

conservando su forma en comparación a las demás variables termodinámicas.

97
CAPÍTULO VIII

REFERENCIAS BIBLIOGRÁFICAS

[1] EINSTEIN, A., Die Grundlage der allgemeinen Relativitatstheorie (on the general

theory ofrelativity), Annalen der Physik, 49, pp.769-822, 1916.

[2] SOKOLNIKOFF, l. S., Análisis Tensorial, Editorial Index-Prial, 1971.

[3] CARMELI, M., Group Theory and General Relativity, McGraw-Hill, London,

1977.

[4] LANDAU, L.D. y LIFSHITZ, E.M., Teoría Clásica de los Campos, Volumen 2 del

Curso de Física Teórica, Editorial Reverté, Barcelona, 2008.

[5] WEINBERG, S., Gravitation and Cosmology: Principies and Applications of the

General Theory of Relativity, John Wiley & Sons, New York, 1972.

[6] MISNER, C., THORNE, K. and WHEELER, J.A., Gravitation, W.H. Freeman and

Company, New York, 1973.

[7] HA WKING, S.W. and ELLIS, G.F.R., The Large Scale Structure of Space-Time,

Cambridge University Press, London, 1973.

[8] TOLMAN R., Relativity Thermodynamics and Cosmology, Oxford University

Press, London, 1934.

[9] PHAM MAU QUAN, Thermodynamique d'un fluide relativiste, Séminare L. de

Broglie, Théories physiques, tome 24 (1954-1955), exp. no 3, p.l-19, 1954.

98
[10] HERRERA, L., DI PRISCO A. and OSPINO, J., Conformally flat anisotropic

spheres in general relativity, arXiv:gr-qc/0 10205 8v 1, 200 l.

[ 11] HERRERA, L., Anisotropic geodesic fluid spheres in general relativity, arXiv:gr-

qc/0207040vl, 2002.

[12] MANSOURI, R, On the non-existence of time-dependen! fluid spheres in general

relativity obeying an equation ofstate, Annales de l'I.H.P., section A, tome 27, n°

2 (1977), p.175-183, 1977.

[13] CULETU, H., Anisotropic fluid with time dependent viscosity coefficients,

arXiv:0711.0062v3 [hep-th], 2008.

[14] BAYIN, S., Anisotropic fluid in general relativity, Physics Review D, Vol 26, pp.

1262-1274, 1982.

[15] SCHWARZSCHILD, K., On the Gravitational Field of a Mass Point According to

Einstein 's Theory, Sitzungsberichte der Koniglich Preussischen Akademie der

Wissenschaften zu Berlin, Phys-Math, 1916, 189-196. (Transl. by S. Antoci, A.

Loinger (1999) from the original paper. physics/9905030).

[16] KRISHNA RAO, J., ANNAPURNA, M. and TRIVEDI, M., Static charged

spheres with anisotropic pressure in general relativity, PRANAMA Joumal of

Physics, Vol. 54, N° 2, pp. 215-225, 2000.

[17] MAK, M.K. and HARKO, T., Anisotropic stars in general relativity, Proceedings

ofthe Royal Society A, 459, pp. 393-408, 2003.

[18] WOODHOUSE, N.M., General Relativity, Springer Inc, Berlin, 2007.

[19] JOSHI, P., Global Aspects in Gravitation and Cosmology, Clarendon Press,

Oxford, 1996.

99
[20] BERGER,S., HOJMAN,R. and SANTAMARINA,J., Joumal of Mathematical,

Vol. 28, No. 2949, 1987.

[21] GAETE, P. and HOJMAN, R., Joumal ofMathematical Physics, Vol. 31, No 140,

1990.

[22] CANUTO, V., Neutron Stars: General Review, Solvay Conference on

Astrophysics and Gravitation, 1973.

[23] HERRERA, L. and PONCE DE LEÓN, J., Jouma1 of Mathematica1 Physics, Vol.

26, No 1114, 1987.

[24] EINSTEIN, A., El Significado de la Relatividad, Planeta-De Agostini, México,

1985.

[25] HACYAN, S., Relatividad Especial para Estudiantes de Física, Fondo de Cultura

Económica, México, 1995.

[26] MATTHEW, S., Principies & Applications of Tensor Analysis, W. Sams & Co.,

1963.

[27] NELSON, E., Tensor Analysis, Princeton University Press, Princeton, 1967.

[28] SOKOLNIKOFF, l. S., Análisis Tensorial, Editorial Index-Prial, 1971.

100
 APÉNDICE A

SÍMBOLOS PARA LAS CANTIDADES FÍSICAS

A Energía libre. Número de moléculas por mol.

e Valor de la velocidad de la luz.

E Energía.

F Potencial termodinámico.

g(t), g Función dependiente del tiempo de un elemento de línea para

modelos cosmológicos homogéneos.

m Masa. Magnitud de un objeto celeste.

N Número de moléculas.

h Constante de Planck.

p Presión.

Q Calor.

101
r Coordenada radial.

R Constante de gas.

S Entropía.

t Tiempo.

T Temperatura.

u Energía.

V Velocidad.

x,y,z Coordenadas espaciales.

p Densidad

102
 APÉNDICEB

CONSTANTES FÍSICAS

Velocidad de la luz en el vacío e= 2.99792458x10 8 rnls

Constante de Planck reducida 1i = b 12tr = 1.054 571596(82)x1 0"34 Js

Constante de los gases R =8.314J 1K

Volumen estándar de un gas V0 = 22.4/1 mol

Carga del electrón (magnitud) e= 1.602176 462(63)xl0" 19 e

= 4.803 204 20 (19)xl0- 10 esu

Masa del electrón me= 0.510 998 902(21)MeV/c 2

= 9.109 38188(72)x1 0"31 kg

Permitividad del espacio libre &0 = 8.854187 817x10- 12 C 2 1N.m 2

Constante de estructura fina 2

a= e 14tr&0 1iC = 11137.035 999 76(50)

Constante de Stefan-Boltzmann K= 1.380x10-6 J 1K

Constante gravitacional G = 6.673(10)xl0- 11 m 3 1kg.s 2

103
 APÉNDICE C

ESCALAS Y UNIDADES NATURALES

La escala de Planck.

Usando combinaciones apropiadas de 1i (li = h/27!), e y G es posible definir

unidades naturales de masa, longitud y tiempo, conocidas como masa de Planck

( mP ), longitud de Plan ek ( zp ) y tiempo de Planek ( Tp ) respectivamente.

{hc" 8
mP = ~G = 2.1767(16)xl0- kg (C.1)

(C.2)

TP = /PfMf =5.390547712x10-
c~7
44
s (C.3)

Unidades seminaturales.

Se define e = 1, de tal modo que:

m = {h (CA)
p ~a·

l =T =-JhG
p p
(C.5)

Equivalentemente:

104
 (C.6)

Con estas unidades se trabaja comúnmente en física de partículas, pero no son

recomendables en relatividad clásica.

Unidades de Planck o unidades naturales.

Se define:

(C.7)

de tal modo que:

m p =1 p =Tp =1 (C.8)

Este sistema de unidades sirve para trabajar tanto en relatividad clásica como en

física de partículas, pero resulta confuso ya que es sencillamente imposible usar el

análisis dimensional para comprobar la consistencia de los resultados que se

obtengan.

105