Apuntes de Relatividad General

IPN ESFM
Instituto Politécnico Nacional
Escuela Superior de Fı́sica y Matemáticas
Relatividad General
Notas del Curso
Apuntes de:
A. Villa
Escrito por:
A. Villa
Abril 2021
Yo 1 Relatividad General
30 de junio de 2023
Índice general
3
IPN ESFM
La Métrica
El producto escalar asociado a una métrica lo definiremos como
⃗v · ⃗v = v 2
Por el Teorema del Coseno
c2 = a2 + b2 − 2abCos (γ)
Considerando Cos (τ − θ) = Cos (θ)
c2 = a2 + b2 + 2abCos (θ)
c · c = a · a + b · b + 2abCos (θ)
(⃗a + ⃗b) · (⃗a + ⃗b) =
a 2 + a · b + b · a + b2 =
a2 + a · b + b · a + b2 = a2 + b2 + 2abCos (θ)
a2 + a · b + b · a + b2 = a2 + b2 + 2abCos (θ)
2a · b = 2abCos (θ)
a · b = abCos (θ)
De esta manera encontramos el producto interno para dos vectores diferentes, por otro lado
Si recordamos que Cos (θ) = Cos (τ − θ) = Cos (τ ) Cos (θ) + Sen (τ ) Sen (θ)
2abCos (θ) = ab (Cos (τ ) Cos (θ) + Sen (τ ) Sen (θ))

= (aCos (τ ))(bCos (θ)) + (aSen (τ ))(bSen (θ))
= a1 b1 + a2 b2
Considerando que ahı́ están las componentes del vector ⃗a y ⃗b
Ejemplo :
Sean E = {ê1 , ê2 } una base euclidiana ortonormal tal que ê1 · ê2 = δi,j
También sea otra base S = {ê1 + ê2 , ê2 } ≡ {u1 , u2 }
⃗ = 2u1 + 3u2 y B
Sea A ⃗ = u1 − u2
Efectuando el producto tenemos
⃗·B
A ⃗ = (2u1 + 3u2 ) · (u1 − u2 )
= 2u1 · u1 − 2u1 · u2 + 3u2 · u1 − 3u2 · u2
Ahora renombremos los productos de tal manera que
g11 = u1 · u1 , g12 = u1 · u2 + 3u2 , g21 = u2 · u1 , g22 = u2 · u2
A este reordenamiento lo llamaremos ”métrica”
g11 = (ê1 + ê2 ) · (ê1 + ê2 ) = (ê1 · ê1 ) + (ê1 · ê2 ) + (ê2 · ê1 ) + (ê2 · ê2 ) = 1 + 0 + 0 + 1 = 2
g12 = (ê1 + ê2 ) · ê2 = (ê1 · ê2 ) + (ê2 · ê2 ) = 0 + 1 = 1
g12 = ê2 · (ê1 + ê2 ) = (ê2 · ê1 ) + (ê2 · ê2 ) = 0 + 1 = 1
g22 = ê2 · ê2 = (ê2 · ê2 ) = 1
⃗ = a1 u1 + a2 u2 y B
Generalicemos un poco A ⃗ = b1 u1 + b2 u2 luego
⃗·B
A ⃗ = a1 b1 g11 + a1 b2 g12 + a2 b1 g21 + a2 b2 g22
IPN ESFM
Para una base euclidiana E la métrica esta dada por

g g12 1 0
gE = 11 =
g21 g22 0 1
Para la base S, la métrica es

g g12 2 1
gS = 11 =
g21 g22 1 1
Para hablar de transporte paralelo vamos a trabajar los vectores en R3 y utilizaremos las
coordenadas esféricas
x = R Sen (θ) Cos (ϕ)

y = R Sen (θ) Sen (ϕ)
z = R Cos (θ)
vamos a recordar que

∂ i
êi = x êj
∂xj
para la primera variable x1 = R
∂ ∂ ∂
êR = xêx + yêy + zêz
∂R ∂R ∂R
= Sen (θ) Cos (ϕ) êx + Sen (θ) Sen (ϕ) êy + Cos (θ) êz
para la segunda variable x2 = θ

∂ ∂ ∂
êθ = xêx + yêy + zêz
∂θ ∂θ ∂θ
= RCos (θ) Cos (ϕ) êx + RSen (θ) Sen (ϕ) êy − RSen (θ) êz
para la segunda variable x3 = ϕ

∂ ∂ ∂
êϕ = xêx + yêy + zêz
∂ϕ ∂ϕ ∂ϕ
= −RSen (θ) Sen (ϕ) êx + RSen (θ) Cos (ϕ) êy
Como podemos recordar estas bases no son ortonormales.

Ahora planteamos la siguiente cuestión, como varia êR , êθ , êϕ , cuando cambian las coordenadas del
punto, es decir como varı́an los vectores base a medida que va cambiando el punto donde estamos
mirando, esto nos va a ayudar a definir la derivada covariante.
Es fácil ver que, para el cambio en coordenadas cartecianas es 0, es decir
∂ ∂ ∂
êx = 0; êx = 0; êx = 0
∂x ∂y ∂z
Esto nos dice que estos vectores son únicos, que para cualquier componente el vector base es único,
en cambio para los vectores base de las coordenadas esféricas si hay cambios en los vectores base,
es decir, para êR
∂
êR = 0
∂R
∂
êR = Cos (θ) Cos (ϕ) êx + Cos (θ) Sen (ϕ) êy − Sen (θ) êz
∂θ
∂
êR = −Sen (θ) Sen (ϕ) êx + Sen (θ) Cos (ϕ) êy
∂ϕ
IPN ESFM
construyamos la metrica de esta nueva base para conocer los vectores base duales.
 
êR · êR êR · êθ êR · êϕ
g =  êθ · êR êθ · êθ êθ · êϕ 
êϕ · êR êϕ · êθ êϕ · êϕ
Calculemos las componentes de la métrica
gRR = êR · êR = 1

gRθ = êR · êθ = 0
gθR = êθ · êR = 0
gRϕ = êR · êϕ = 0
gϕR = êϕ · êR = 0
gθϕ = êθ · êϕ = 0
gϕθ = êϕ · êθ = 0
Estos cálculos se pueden realizar, ahora
gθθ = θ · θ
= R2 Cos2 (θ)Cos2 (ϕ) + R2 Cos2 (θ)Sen2 (ϕ) + R2 Sen2 (θ)
= R2 Cos2 (θ) + R2 Sen2 (θ) = R2
gϕϕ = ϕ · ϕ
= R2 Sen2 (θ)Sen2 (ϕ) + R2 Sen2 (θ)Cos2 (ϕ) + R2 Sen2 (θ)
= R2 Sen2 (θ)
Longitud de arco
Con los resultados de la seccion anterior tenemos definido lo que es una métrica, luego para
pequeños incrementos ∆S sobre un arco continuo, si tomamos dos vectores R⃗ 1, R
⃗ 2 tal que
⃗1 − R
R ⃗ 2 = (xf − xi )ê1 + (yf − yi )ê2 = ∆xê1 + ∆yê2 = ∆S
Entonces
(∆S)2 = (∆xê1 + ∆yê2 ) · (∆xê1 + ∆yê2 )

= (∆x)2 (ê1 · ê1 ) + ∆x∆y(ê1 · ê2 ) + ∆y∆x(ê2 · ê1 ) + (∆x)2 (ê2 · ê2 )
= (∆x)2 g11 + ∆x∆yg12 + ∆y∆xg21 + (∆x)2 g22
= (∆x)2 g11 + 2∆x∆yg12 + (∆x)2 g22
Para calcular la magnitud de esos pequeños incrementos, si tomamos los incrementos como
diferenciales, entonces
Zb Zb p
dλ
Longitud = dS = g11 (dx)2 + 2g12 dxdy + g22 (dy)2 ·
dλ
a a
s
Z 2 2
dx dx dy dy
= g11 + 2g12 + g22 · dλ
dλ dλ dλ dλ
Ejemplo Sea y = x2 en el intervalo [0, 1], calcular la longitud de arco.

Sol.
IPN ESFM
Parametrizando la función (
x=λ
y = x2 = λ 2
Después derivamos respecto al parámetro entonces
dx


 =1
dλ
 dy = 2λ

dλ
Sustituyendo en la ecuación de longitud de arco
λ=1s
Z 2 2
dx dx dy dy
Longitd = g11 + 2g12 + g22 · dλ
dλ dλ dλ dλ
λ=0
λ=1
Z
s 2 2
dx dx dy dy
= g11 + 2
g
12 + g22 · dλ
dλ dλ dλ dλ
λ=0
λ=1
Z r
dx dy
= (1)(1)2 + 2
g
12 + (1)(2λ)2 · dλ
dλ dλ
λ=0
λ=1
Z r
dx dy
= (1)(1)2 + 2
g
12 + (1)(2λ)2 · dλ
dλ dλ
λ=0
λ=1
Z p
= 1 + 4λ2 · dλ
λ=0
El parámetro va de 0 a 1, luego integrando

λ=1
Z p
Longitud = 1 + 4λ2 · dλ
λ=0
1
1 p 1 p
= ln 2λ + 4λ2 + 1 + λ 4λ2 + 1
4 2 0
1 √ 1√
= ln 5+2 + 5 = 1,4789...
4 2
La metrica es la manera de medir longitudes, mediante el producto interno entre vectores.
Base Dual
En álgebra lineal, una base dual es un conjunto de vectores que forman una base para el espacio
dual de un espacio vectorial. Para un espacio vectorial V de dimensiones finitas, el espacio dual V ∗
es isomorfo a V y para cualquier conjunto dado de vectores base {ê1 , · · · , ên } de V , hay asociada
una base dual {ê1 , · · · , ên } de V ∗ con la relación
êi · êj = δij
Esta condición va a cumplir independientemente de la métrica que tengamos, y para calcular la
base dual tenemos que calcular ( 1
ê = aê1 + bê2
ê2 = cê1 + dê2
IPN ESFM
Usando la condición que definimos
ê1 · ê1 = (aê1 + bê2 ) · ê1 = a(ê1 · ê1 ) + b(ê2 · ê1 ) = ag11 + bg21 = 1 · · · (I)
ê1 · ê2 = (aê1 + bê2 ) · ê2 = a(ê1 · ê2 ) + b(ê2 · ê2 ) = ag12 + bg22 = 0 · · · (II)
ê2 · ê1 = (cê1 + dê2 ) · ê1 = c(ê1 · ê1 ) + d(ê2 · ê1 ) = cg11 + dg21 = 0 · · · (III)
ê2 · ê2 = (cê1 + dê2 ) · ê2 = c(ê1 · ê2 ) + d(ê2 · ê2 ) = cg12 + dg22 = 1 · · · (IV )
Para resolver, multiplicamos (I) por g12 y a (II) por −g11 y sumamos
ag12 g11 + bg12 g21 = g12

−ag11 g12 − bg11 g22 = 0
b(g12 g21 − g11 g22 ) = g12
Con este resultado podemos notar que es un determinante g12 g21 − g11 g22 = − |g|
De modo que
1
−b |g| = g12 ⇒ b = − g12
|g|
Luego tenemos de la ecuación (II) tenemos
1 1
ag12 − g12 g22 = 0 ⇒ a= g22
|g| |g|
Para c y d multiplicamos (IV ) por g11 y a (III) por −g12 y sumamos
−cg12 g11 − cg12 g21 = 0

cg11 g12 + dg11 g22 = g11
d(g11 g22 − g12 g21 ) = g11
Con este resultado podemos notar que es un determinante g11 g22 − g12 g21 = |g|
De modo que
1
d |g| = g11 ⇒ d = g11
|g|
Luego tenemos de la ecuación (III) tenemos
1 1
cg11 + g11 g21 = 0 ⇒ c=− g21
|g| |g|
Luego si ordenamos las constantes de la manera que

a b 1 g22 −g12 −1
≡ g ij

= = gij
c d |g| −g21 g11
De donde podemos tener que

ij 1 g22 −g12
g =
|g| −g21 g11
Ahora con el resultado de las constantes podemos concluir que
1 1
ê1 = aê1 + bê2 = g22 ê1 − g12 ê2 = g 11 ê1 + g 12 ê2
|g| |g|
1 1
ê2 = cê1 + dê2 = − g21 ê1 + g11 ê2 = g 21 ê1 + g 22 ê2
|g| |g|
Para reducir la escritura de ecuaciones usaremos la notación de Einstein

X 
ê1 = g 11 ê1 + g 12 ê2 = g 1j 


j
 X
X ⇒ êi = g ij êj = g ij êj
ê2 = g 21 ê1 + g 22 ê2 = g 2j 

 j

j
IPN ESFM
Por lo tanto la basa dual esta definida como
êi = g ij êj
êi · êj = δji
Ejemplo Capitulo 2, 27:22

Sean dos bases S1 = {ê1 , ê2 } y S2 = {ê1 , ê2 } Sol.
Vectores covariantes/contravariantes
Sea una base S = {ê1 , ê2 } con una métrica

2 1
gij =
1 2
Como podemos ver que la métrica es diferente a la cartesiana, luego calculemos la matriz inversa
de la métrica gij , entonces
−1 1 g11 −g21
gij =
|g| −g12 g22
Calculando el determinante de gij tenemos
|gij | = (g11 g22 − g12 g21 ) = (2)(2) − (1)(1) = 3
Tenemos que
−1 1 2 −1
gij = = g ij
3 −1 2
Entonces podemos calcular
2 1
ê1 = g 11 ê1 + g 12 ê2 =
ê1 − ê2
3 3
2 21 22 1 2
ê = g ê1 + g ê2 = − ê1 + ê2
3 3
Ahora consideremos un vector
⃗ = 2ê1 + 3ê2
A
y lo expresamos en la base dual, para esto tenemos que invertir los vectores unitarios y expresarlos
en la base dual, para esto tenemos que

1 2
2 − ê1 + ê2 = 2(ê2 )
3 3
2 1
ê1 − ê2 = ê1
3 3
ê2 = ê1 + 2ê2
Ahora recordemos que êi = g ij êj es fácil notar que
gij (êi ) = gij (g ij êj )

−1
= gij (gij êj )
= Iij êj = êj
Entonces, cambiando los indices mudos
êi = gij êj
IPN ESFM
Con esta transformación de vectores para ê2 tendremos que

ê2 = g2j êj
= g21 ê1 + g22 ê2
= (1)ê1 + (2)ê2
⃗ en
lo cual nos deja el mismo resultado para el vector ê2 , ahora calculemos ê1 para encontrar el A
la base dual
ê1 = g1j êj
= g11 ê1 + g12 ê2
= (2)ê1 + (1)ê2
Luego
⃗ = 2ê1 + 3ê2
A
= 2(g11 ê1 + g12 ê2 ) + 3(g21 ê1 + g22 ê2 )
= (2g11 + 3g21 )ê1 + (2g12 + 3g22 )ê2
= (2(2) + 3(1))ê1 + (2(1) + 3(2))ê2
= 7ê1 + 8ê2
Entonces tenemos que el vector, tiene una representasion para la base y su dual.
⃗ = 2ê1 + 3ê2
A
= 7ê1 + 8ê2
En el caso de una base cartesiana hay ninguna diferencia entre la base y la dual. Ahora comohay
para un mismo vector dos representasiones, hay que darle una notacion introduciendo la siguiente
nomeclatura, tomando el vector anterior de ejemplo
⃗ = 2ê1 + 3ê2 = A1 ê1 + A2 ê2
A
= 7ê1 + 8ê2 = A1 ê1 + A2 ê2
⃗
y llamaremos a {A1 , A2 , · · · AN } Componentes Contravariantes de A
⃗
y a {A1 , A2 , · · · AN } Componentes Covariantes de A
con estos a todas las componentes Contravariantes le acompañara una base de vectores Covariantes
o Covectores, y para las componetes Covariantes se asocia una base Contravariantes.
Base Coordenada
Sea una función f (x) = x2 + 3, y calculemos la derivada dx
df , para esto despejamos x
p
x(f ) = f − 3
dx 1 1 1
= √ = √ =
df 2 f −3 2
2 x +3−3 2x
Esto se puede calcular porque x se puede despejar.
Tambien se puede calcular de la forma
f = x2 (f ) + 3
Aquı́ representamos que x depende de f , entonces
d d
(f ) = (x2 (f ) + 3)
df df
dx
1 = 2x · +0
df
dx 1
=
df 2x
IPN ESFM
Esto es derivación implı́cita.
Ahora veamos un caso un poco mas complejo, sean las funciones
f (x, y) = x2 + y 2
g(x, y) = 2x + 3y
∂ ∂ ∂ ∂
Ahora imaginemos que que nos piden que ∂x f, ∂x g, ∂y f, ∂y g, para lo que debemos considerar la
f y g como variables de tal forma que
f = x2 (f, g) + y 2 (f, g) · · · · · · (1)

g(x, y) = 2x(f, g) + 3y(f, g) · · · · · · (2)
Entonces
∂ ∂x2 (f, g) ∂y 2 (f, g)
f= +
∂f ∂f ∂f
∂ ∂
1 = 2x · x + 2y · y
∂f ∂f
∂ ∂2x(f, g) ∂3y(f, g)
g= +
∂f ∂f ∂f
∂ ∂
0=2· x+3· y
∂f ∂f
∂ ∂x2 (f, g) ∂y 2 (f, g)

f= +
∂g ∂g ∂g
∂ ∂
0 = 2x · x + 2y · y
∂g ∂g
∂ ∂2x(f, g) ∂3y(f, g)
g= +
∂f ∂f ∂f
∂ ∂
1=2· x+3· y
∂g ∂g
Reescribamos estas ecuaciones de la siguiente forma particular, considerando los siguientes arreglos.
∂ ! !
∂ ∂ ∂ ∂ ∂
2x 2y ∂f x ∂g x 2x ∂f x + 2y ∂f y 2x ∂g x + 2y ∂g y
∂ ∂ = ∂ ∂ ∂ ∂
2 3 ∂f y ∂g y 2 ∂f x + 3 ∂f y 2 ∂g x + 3 ∂g y

1 0
=
0 1
Si recordamos que
A·B =I
−1
A (A · B) = A−1 (I)
B = A−1
Podemos notar que !

∂ ∂ −1
∂f x ∂g x 2x 2y
∂ ∂ =
∂f y ∂g y
2 3
Calculando esta matriz inversa

−1
2x 2y 1 3 −2y
=
2 3 6x − 4y −2 2x
IPN ESFM
Con estos resultados podemos proponer 2 bases que cumplen que
êa = 4êx + 5êy

êb = 7êx + 8êy
escribamos la matriz asociada a este sistema

−1
4 5 1 8 −5
=−
7 8 3 −7 4
Entonces la transformación se da por

8 5
êx = − êa + êb
3 3
7 4
êy = êa − êb
3 3
Con todo esto regresamos al producto y lo renombramos por sus derivadas
∂ !
∂
2x 2y ∂f x ∂g x 1 0
∂ ∂ =
2 3 ∂f y ∂g y 0 1
! !
∂ ∂ ∂ ∂
∂x f ∂y f ∂f x ∂g x 1 0
∂ ∂ ∂ ∂ =
∂x g ∂y g ∂f y ∂g y 0 1
Como vimos en el ejemplos numéricos, la columnas de las matrices están asociadas las componentes
de las respectivas bases, de manera general tenemos
  z}|{ z}|{
êx êy
z}|{ z}|{ ∂ ∂
 f ∂ f   ∂f x ∂g x
 ∂  
1 0
 ∂x ∂y   ∂ ∂ =
  
 ∂ ∂   y y 0 1
 g g   ∂f ∂g 
∂x ∂y |{z} |{z}
|{z} |{z} êf êg
Ahora ya podemos hablar bien de lo que es una base coordenada

si partimos de un conjunto de coordenadas S1 = {x, y} y tenemos definida una transformacion
lineal de coordenadas a S2 {f, g}, los vectores unitarios de la base S1 por la misma transforma,
generan unos vectores unitarios de la base S2 y estarán dados como
∂ ∂
êf ≡ xêx + yêy
∂f ∂f
∂ ∂
êg ≡ xêx + yêy
∂g ∂g
∂ ∂
êx ≡ f êf + gêg
∂x ∂x
∂ ∂
êy ≡ f êf + gêg
∂y ∂y
Ahora definiremos la NOTACIÓN OFICIAL del calculo tensorial

Tendremos las coordenadas origen
{x1 , x2 , · · · , xN }
con su base asociada
{ê1 , ê2 , · · · , êN }
Las coordenadas generadas serán
′ ′ ′
{x1 , x2 , · · · , xN }
IPN ESFM
con su base asociada

{ê1′ , ê2′ , · · · , êN ′ }
Entonces definiremos los vectores primados
∂ 1 ∂ 2 ∂ N
ê1′ = x ê1 + x ê2 + · · · + x êN
∂x1′ ∂x1′ ∂x1′
∂ i
= x êi
∂x1′
Entonces en general tendremos que
∂ j
êi′ = x êj
∂xi′
∂ j′
êi = x êj ′
∂xi
Tensores
Para empezar a introducir tensores definiremos unas coordenadas de la siguiente manera
a=x+y
b = 4x + y
Con la matriz inversa podemos despejar x y y (el calculo esta en corto)

1 1
x=− a+ b
3 3
4 1
y = a+− b
3 3
Por los resultados de ”BASE COORDENADA”tenemos que
∂ ∂ 1 4
êa = xêx + yêy = − êx + êy
∂a ∂a 3 3
∂ ∂ 1 1
êb = xêx + yêy = êx − êy
∂b ∂b 3 3
Esto nos dice que las bases de coordenadas, no dependen del punto, no depende de la variable.
Por ejemplo, sea un vector

⃗ = 2êx + êy = αêa + β êb
A
⃗ es el mismo en las dos bases, para saber cuando valen las componentes
En magnitud el vector A
α y β sustituiremos
⃗ = αêa + β êb
A

1 4 1 1
= α − êx + êy + β êx − êy
3 3 3 3

1 1 4 1
= − α + β êx + α − β êy = 2êx + êy
3 3 3 3
Esto nos da un sistema de dos ecuaciones, que solucionando tenemos que α = 3 y β = 9, por lo
que
⃗ = 2êx + êy = 3êa + 9êb
A
Ahora veamos su se puede hacer este calculo con la notación que ya conocemos y lo haremos de la
siguiente forma
⃗ = 2êx + êy = 3êa + 9êb
A
= Ax êx + Ay êy = Aa êa + Ab êb
IPN ESFM
Recordando la transformación que tenı́amos para las bases, veamos si funciona para las compo-
nentes tomando a = x + y y b = 4x + y
∂ ? ∂
Aa = aAx + aAy = (1)Ax + (1)Ay = 2 + 1 = 3
∂x ∂y
? ∂ ∂
Ab = bAx + bAy = (4)Ax + (1)Ay = 8 + 1 = 9
∂x ∂y
Con esto logramos mostrar que también nos ayuda a calcular las componentes de un vector y mas
aun no ayuda a comprobar que un vector queda invariante, independiente de la base que trabajemos
y si el elemento queda invariante bajo transformaciones de esta manera entonces tendremos un
”TENSOR la tansformacion bajo la convención de Einstein serán
2
∂ i′ j
 ′
 Ai =
 x A → Tipo(1, 0)
TENSOR ∂xj
A ′ = ∂ xi′ A → Tipo(0, 1)

i j
∂xj
SUBIR Y BAJAR INDICES (vectores covariantes y contravariantes) Para empezar,

subir y bajar indices no es un cambio de coordenadas, es un cambio de base, es cambiar de la base
a la base dual.
Lo que buscamos es, dado un vector A ⃗ que este en componentes contravariantes,convertirlo a com-
ponentes covariantes, y para esto utilizaremos la métrica gij .
Ejemplo: Utilizando las mismas coordenadas {a, b} y {x, y}

Calculemos la métrica de {êa , êb }

gaa gab
gab =
gba gbb
Luego
1 1
gaa = êa · êa = (−êx + 4êy ) · (−êx + 4êy ) = (1 + 16)
9 9
1 1
gbb = êb · êb = (êx − êy ) · (êx − êy ) = (1 + 1)
9 9
1 1
gab = êa · êb = (−êx + 4êy ) · (êx − êy ) = (−1 − 4)
9 9
con esto el la métrica es

1 17 −5
gab =
9 −5 2
De aquı́ la matriz inversa es

2 5
g ab =
5 17
⃗ = 3êa + 9êb = Aa êa + Ab êb ,
Trabajamos en una única base, {a, b} y tomamos de nuevo el vector A
calculemos recordando el cambio de base
êa = g aa êa + g ab êb = 2êa + 5êb

êb = g ba êa + g bb êb = 5êa + 17êb
Con esto
⃗ = Aa (2êa + 5êb ) + Ab (5êa + 17êb )

A
3 9
z }| { z }| {
= (2Aa + 5Ab ) êa + (5Aa + 17Ab ) êb
IPN ESFM
2 1
y con esto tenemos un sistema de ecuaciones que al resolverlo tendremos que Aa = 3 y Ab = 3
Ahora calculemos lo mismo pero utilizando la métrica Ai = gij Aj

a b 17 5 17 15 2
Aa = gaa A + gab A = 3+ − 9= − = ✓
9 9 3 3 3

5 2 15 18 1
Ab = gba Aa + gbb Ab = − 3+ 9=− + = ✓
9 9 9 3 3
⃗ en la base cartesiana y la comparamos en la otra base

Ahora calculemos la longitud del vector A
⃗·A
A ⃗ = (2êx + êy ) · (2êx + êy ) = 4 + 1 = 5
Para la otra base
A⃗·A⃗ = (3êa + 9êb ) · (3êa + 9êb )
= 9(êa · êa ) + 27(êa · êb ) + 27(êb · êa ) + 81((êb · êb ))

17 5 5 2
=9 + 27 − + 27 − + 81
9 9 9 9
= 17 − 15 − 15 + 18 = 5⌟//
Ahora con la notacion del calculo tensorial tenemos
⃗·A⃗ = Ai êi · Aj êj

A
= Ai Aj (êi · êj )
= Ai Aj gij = Ai Ai = Aj Aj
Como
⃗ = 3êa + 9êb = 2 êa + 1 êb
A
3 3
El producto queda
⃗·A
A ⃗ = Ai Ai
= Aa Aa + Ab Ab
2 1
=3· +9· =2+3=5
3 3
Para generalizar mas los tensores vamos a hacer uso de la Fuerza de Lorentz F⃗ = ⃗v × B,
⃗ lo veremos
⃗
de modo muy matemático y para esto consideramos a B un valor fijado, entonces calculamos
F⃗ · ⃗n = ⃗v × B
⃗ · ⃗n = const
Vamos a llamar a esto T, esto es un objeto matematico, que admite dos vectores tales que:
T(⃗v , ⃗n) = const
ahora
v · B = v 1 ê1 + v 2 ê2 + v 3 ê3 × B 1 ê1 + B 2 ê2 + B 3 ê3

= v 1 B 2 ê3 − v 1 B 3 ê2 − v 2 B 1 ê3 + v 2 B 3 ê1 + v 3 B 1 ê2 − v 3 B 2 ê1

= v 2 B 3 − v 3 B 2 ê1 + v 3 B 1 − v 1 B 3 ê2 + v 1 B 2 − v 2 B 1 ê3

Luego
(v · B) · ⃗n = v 2 B 3 − v 3 B 2 n1 + v 3 B 1 − v 1 B 3 n2 + v 1 B 2 − v 2 B 1 n3

organicemos estos numeros de manera interesante

T12 T13 T21 T23 T32 T31
q q q q q q
T(⃗v , ⃗n) = −B 3 v 1 n2 + B 2 v 1 n3 + B 3 v 2 n1 − B 1 v 2 n3 + B 1 v 3 n2 − B 2 v 3 n1
= T12 v 1 n2 + T13 v 1 n3 + T21 v 2 n1 + T23 v 2 n3 + T32 v 3 n2 + T31 v 3 n1
0 0 0
q q q
= T11 v 1 n1 1 2 1 3
+ T12 v n + T13 v n + T21 v n + 2 1
T22 v 2 n2 + T23 v n + T32 v n + T31 v n + T33 v 3 n3
2 3 3 2 3 1
= Tij v i nj
IPN ESFM
Un tensor (al menos de orden 2) es una función de dos vectores que se expresa de esta forma
T(⃗v , ⃗n) = Tij v i nj
ademas si hay un cambio de coordenadas, las componentes del tensor debe cumplir con el cambio
∂ i′ ∂ j ′
Tij = x x Ti′ j ′
∂xi ∂xj
Producto Tensorial de Espacios Vectoriales

Si V es un espacio vectorial, entonces estudiaremos el producto tensorial · ⊗ · entre dos espacios
vectoriales. Para Relatividad nos va ser util estudiar el producto tensorial de V con V , (V ⊗ V ) y
esto va a seguir siendo un nuevo espacio vectorial, que buscaremos definir apartir de su base.
Sea {ê1 , ê2 } una base de V , entonces construiremos la base de V ⊗ V como

{ê1 ⊗ ê1 , ê1 ⊗ ê2 , ê2 ⊗ ê1 , ê2 ⊗ ê2 }
Esta base no ayudara a definir nuevos vectores del espacio V ⊗ V que se escribirán como
⃗ = 5ê1 ⊗ ê2 + 6ê2 ⊗ ê1 + 3ê2 ⊗ ê2
A
Recalco que la base en un espacio vectorial son ê1 , ê2 y sus elementos son combinaciones lineales
con estos elementos, para este nuevo espacio vectorial tendremos que su base son estos nuevos
entes matemáticos ê1 ⊗ ê1 , ê1 ⊗ ê2 , ê2 ⊗ ê1 , ê2 ⊗ ê2 y todos sus elementos deben ser combinaciones
lineales con los mismos, no son términos separados, se convierten en un mismo ente, y los nuevos
elementos de este espacio vectorial se les llama TENSORES.
A este ejemplo anterior se le llama TENSOR de rango 2, porque es el producto de 2 espacios
vectoriales.
Habiamos definido un tensor de manera diferente en las secciones anteriores, y eso se sigue man-
teniendo igual, lo que pasa que ahora estamos dando una nueva notación matematica que nos
permita trabajarlo de manera mas intrı́nseca.
Propiedades
1. (a + b) ⊗ c = a ⊗ c + b ⊗ c multilineal
2. a ⊗ c ̸= c ⊗ a No conmuta
3. ê1 ⊗ ê2 no es una operación, es un elemento base del espacio vectorial
Recordando que tenemos una definición del producto punto êi · êj = gij , definiremos una analogı́a
entre los nuevos entes que definimos.
êi ⊗ êj · êk ⊗ êm ≡ (êi · êk ) (êj · êm ) = gik gjm
Ahora podemos generalizar esto de la siguiente forma

⃗ ⊗ B)
(A ⃗ · (C
⃗ ⊗ D)
⃗ = (Ai êi ⊗ B j êj ) · (C m êm ⊗ Dp êp )
= Ai B j C m Dp (êi ⊗ êj ) · (êm ⊗ êp )
= Ai B j C m Dp gim gjp
= Ai (C m gim )B j (Dp gjp )
= Ai Ci B j Dj = (A ⃗ · C)(
⃗ B ⃗ · D)⌟
⃗ //
Ejemplo (Tensores)
Sea T ∈ y tomando la base dual de Tp M ⊗ Tp M
T = 5ê1 ⊗ ê2 + 2ê2 ⊗ ê2
Si recordamos que definimos

⃗ B)
T(A, ⃗ = Tij Ai B j
IPN ESFM
Ahora calculemos
⃗ ⊗ B)
T · (A ⃗ Como podemos ver ⃗ ⊗ B)
(A ⃗ ∈ Tp M ⊗ Tp M
Entonces estamos calculando un producto punto entre dos elementos del espacio Tp M ⊗ Tp M pero
en este caso T no esta definido como un producto tensorial, por lo que debemos hacer un paso
extra.
⃗ ⊗ B)
T · (A ⃗ = (5ê1 ⊗ ê2 + 2ê2 ⊗ ê2 ) · (A
⃗ ⊗ B)
⃗
⃗ ⊗ B)
= 5(ê1 ⊗ ê2 ) · (A ⃗ + 2(ê2 ⊗ ê2 ) · (A
⃗ ⊗ B)
⃗
= 5(ê1 ⊗ ê2 ) · (Ai êi ⊗ B j êj ) + 2(ê2 ⊗ ê2 ) · (Am êm ⊗ B p êp )
z }| { z }| {
= 5Ai B j (ê1 ⊗ ê2 ) · (êi ⊗ êj ) +2Am B p (ê2 ⊗ ê2 ) · (êm ⊗ êp )
= 5Ai B j (ê1 · êi )(ê2 · êj ) + 2Am B p (ê2 · êm )(ê2 · êp )
= 5Ai B j (δi1 )(δj2 ) + 2Am B p (δm
2
)(δp2 )
= 5(Ai δi1 )(B j δj2 ) + 2(Am δm
2
)(B p δp2 )
= 5A1 B 2 + 2A2 B 2 ⌟//
Pero si esto es un tensor deberia tener la forma de Tij Ai B j entonces
⃗ ⊗ B)
T · (A ⃗ = 5A1 B 2 + 2A2 B 2
T11 T12 T21 T22
q q q q
= 0 A1 B 1 + 5 A1 B 2 + 0 A2 B 1 + 2 A2 B 2
= T11 A1 B 1 + T12 A1 B 2 + T21 A2 B 1 + T22 A2 B 2
= Tij Ai B j ⌟//
Con este resultado podemos decir algunas cosas, por ejemplo, no podemos definir este producto
T·A⃗ porque A⃗ es un elemento de Tp M y no de Tp M ⊗ Tp M
Esto quiere decir que no podemos hacer aperaciones entre elementos de diferentes rangos, no, si
podremos hacer esas operaciones, seran llamadas CONTRACCIONES.
Ahora veremos la condicion que le pedimos al tensor, que debe transformarse de manera espe-
cifica, y lo mostraremos con un ejemplo
Ejemplo Sea T = 5ê1 ⊗ ê2 + 2ê2 ⊗ ê2 ∈ Tp M ⊗ Tp M y recordando las coordenadas que ya
trabajamos
a=x+y

gaa gab 1 17 −5 −1 2 5
g= = ; g =
b = 4x + y gba gbb 9 −5 2 5 17
Ahora ajustaremos a nuestro tensor a las coordenadas que tenemos
T = 5êx ⊗ êy + 2êy ⊗ êy

Txx Txy Tyx Tyy
q q q q
x x x y
= 0 ê ⊗ ê + 5 ê ⊗ ê + 0 ê ⊗ ê + 2 êy ⊗ êy y x
= Txx êx ⊗ êx + Txy êx ⊗ êy + Tyx êy ⊗ êx + Tyy êy ⊗ êy ⌟//
Queremos que este TENSOR sea el mismo, lo veamos desde el sistema de coordenadas que lo
veamos, entonces si lo describimos en las coordenadas {êa , êb } tendremos un elemento de esta
forma
? ? ? ?
T = Taa êa ⊗ êa + Tab êa ⊗ êb + Tba êb ⊗ êa + Tbb êb ⊗ êb
No sabemos que valor tienen los términos de este TENSOR en esta nueva base, pero calcularemos
sus valores para que permanezca invariable la longitud del elemento y sabemos hacerlo porque ya
IPN ESFM
tenemos todas las definiciones, sabemos que

1 1
q q
a ∂ ∂
ê = aêx + aêy = êx + êy
∂x ∂y
4 1
q q
b ∂ x ∂ y
ê = bê + bê = 4êx + êy
∂x ∂y
Luego veamos que sera
êa ⊗ êa = (êx + êy ) ⊗ (êx + êy )

= êx ⊗ êx + êx ⊗ êy + êy ⊗ êx + êy ⊗ êy
êa ⊗ êb = (êx + êy ) ⊗ (4êx + êy )

= 4êx ⊗ êx + êx ⊗ êy + 4êy ⊗ êx + êy ⊗ êy
êb ⊗ êa = (4êx + êy ) ⊗ (êx + êy )

= 4êx ⊗ êx + 4êx ⊗ êy + êy ⊗ êx + êy ⊗ êy
êb ⊗ êb = (4êx + êy ) ⊗ (4êx + êy )

= 16êx ⊗ êx + 4êx ⊗ êy + 4êy ⊗ êx + êy ⊗ êy
Luego el TENSOR lo sera
T = Taa êa ⊗ êa + Tab êa ⊗ êb + Tba êb ⊗ êa + Tbb êb ⊗ êb
= Taa (êx ⊗ êx + êx ⊗ êy + êy ⊗ êx + êy ⊗ êy )

+ Tab (4êx ⊗ êx + êx ⊗ êy + 4êy ⊗ êx + êy ⊗ êy )
+ Tba (4êx ⊗ êx + 4êx ⊗ êy + êy ⊗ êx + êy ⊗ êy )
+ Tbb (16êx ⊗ êx + 4êx ⊗ êy + 4êy ⊗ êx + êy ⊗ êy )
= (Taa + 4Tab + 4Tba + 16Tbb )êx ⊗ êx
+ (Taa + Tab + 4Tba + 4Tbb )êx ⊗ êy
+ (Taa + 4Tab + Tba + 4Tbb )êy ⊗ êx
+ (Taa + Tab + Tba + Tbb )êy ⊗ êy
Esto nos da un sistema de ecuaciones
Taa + 4Tab + 4Tba + 16Tbb = 0

Taa + Tab + 4Tba + 4Tbb = 5
Taa + 4Tab + Tba + 4Tbb = 0
Taa + Tab + Tba + Tbb = 2
Solucionando el sistema tenemos que

4 1 4 1
Taa = ; Tab = − ; Tba = ; Tbb = −
3 3 3 3
Asi pues tenemos que
4 a 1 4 1
T= ê ⊗ êa − êa ⊗ êb + êb ⊗ êa − êb ⊗ êb
3 3 3 3
IPN ESFM
Ahora escribiremos la regla que siguen estas transformaciones
∂ k ∂
Ti′ j ′ = x mTkm
∂xi′ ∂xj ′
Como ya tenemos despejado los valores de x y y
a b
x=− +
3 3
a b
y=4 −
3 3
Con esto, reemplazamos en la formula y tenemos
∂ k ∂ m
Taa = x x Tkm
∂a ∂a
∂ ∂ m ∂ ∂ m
= x x Txm + y x Tym
∂a ∂a ∂a ∂a
∂ ∂ ∂ ∂ ∂ ∂
= x xTxx + yTxy + y xTyx + yTyy
∂a ∂a ∂a ∂a ∂a ∂a
∂ ∂ + ∂ x ∂ yTxy + ∂ y ∂ x Tyx ∂ ∂
= x x Txx + y yTyy
∂a ∂a ∂a ∂a ∂a ∂a ∂a ∂a
∂ ∂ ∂ ∂
= x yTxy + y yTyy Derivando las funciones
∂a
∂a
∂a
∂a
1 4 4 4
= − 5+ 2
3 3 3 3

20 32 4
= − + =
9 9 3
∂ k ∂ m
Tab = x x Tkm
∂a ∂b
∂ ∂ + ∂ x ∂ yTxy + ∂ y ∂ x Tyx ∂ ∂
= x x Txx + y yTyy
∂a ∂b ∂a ∂b ∂a ∂b ∂a ∂b
∂ ∂ ∂ ∂
= x yTxy + y yTyy Derivando las funciones
∂a ∂b
∂a ∂b
1 1 4 1
= − − 5+ − 2
3 3 3 3

5 8 1
= − =−
9 9 3
∂ k ∂ m
Tba = x x Tkm
∂b ∂a
∂ ∂ + ∂ x ∂ yTxy + ∂ y ∂ x Tyx ∂ ∂
= x x Txx + y yTyy
∂b ∂a ∂b ∂a ∂b ∂a ∂b ∂a
∂ ∂ ∂ ∂
= x yTxy + y yTyy
∂b ∂a
∂b ∂a

1 4 1 4
= 5+ − 2
3 3 3 3

20 8 4
= − =
9 9 3
IPN ESFM
∂ k ∂ m
Tbb = x x Tkm
∂b ∂b
∂ ∂ + ∂ x ∂ yTxy + ∂ y ∂ x Tyx ∂ ∂
= x x Txx + y yTyy
∂b ∂b ∂b ∂b ∂b ∂b ∂b ∂b
∂ ∂ ∂ ∂
= x yTxy + y yTyy
∂b ∂b
∂b ∂b
1 1 1 1
= − 5+ − − 2
3 3 3 3

5 2 1
= − + =−
9 9 3
Estos resultados nos garantiza que T es un elemento independiente del sistema de coordenadas.
Contracción
Cuando hablemos o indiquemos que Tij es un tensor o que Ak es un vector, estamos abusando de
⃗ = Ak êk respectivamente.
la notacion y sobreentenderemos que hablamos de T = Tij êi ⊗ êj y de A
k j
Ahora tenemos Tij y A , queremos calcular Tij A y por tanto esto seria
Tij Aj = Ti1 A1 + Ti2 A2
Tenemos
T = T11 ê1 ⊗ ê1 + T12 ê1 ⊗ ê2 + T21 ê2 ⊗ ê1 + T22 ê2 ⊗ ê2
Y
⃗ = A1 ê1 + A2 ê2
A
Si hacemos TA,⃗ (y no estamos haciendo ningún producto), solo estamos juntando estos dos ele-
mentos, podrı́amos crear una nomeclatura para realizar una contracción pero no lo veo necesario,
daremos por hecho que Tij Aj = Bi es un tensor, por lo que
B = B1 ê1 + B2 ê2 = (T11 A1 + T12 A2 )ê1 + (T21 A1 + T22 A2 )ê2
Ahora hablemos un poco de TENSORES DE RANGO SUPERIOR, por ejemplo un tensor de

rango 3 seria
T = Tijk êi ⊗ êj ⊗ êk ∈ TP M ⊗ TP M ⊗ TP M
Entonces tenemos que
êi ⊗ êj ⊗ êk · (ên ⊗ êm ⊗ êp ) = g in g jm g kp

Tambien los tensores mixtos se pueden anotar como
T = T ijk êi ⊗ êj ⊗ êk Tipo (1, 2)
Si quisiéramos bajar ese indice i igual necesitarı́amos la métrica
Tpjk ≡ gpi T ijk
se desarrolla de esta manera
Tpjk = gp1 T 1jk + gp2 T 2jk Para dim = 2
Paremos a notar algo muy especifico en estos tensores que escribimos, y es que deje un espacio
bajo el indice i, este lugar es porque cada indice indica el orden del espacio vectorial del que se
forma la base.
Tomemos el mismo tensor, pero ahora subamos el indice k
T ij p = g kp T ijk
Solo si el tensor es simétrico, no es importante esta separación, (no se si me explico, pero el orden
y lugar del indice, ayuda a identificar el orden de la base, espero poder explicar mejor conforme
avancen mis escritos).
IPN ESFM
Geodésicas
Vamos a empezar con un pequeño ejemplo en R2
y = ax + b, sı́ P (0, 0)
0 = a(0) + b → b=0
sı́ P (2, 1)
1 = a(2) → a = 1/2
1
y= x
2
Figura 1: *
Ejemplo para la curva parabolica proponemos
y = ax2 + bx + c, sı́ P (0, 0)

2
0 = a(0) + b(0) + c → c=0
sı́ P (2, 1)
1 = a(2)2 + b(2) → a = 1/4
1
y = x2 b = 0 Porque nace en el origen
4
Ahora parametrizaremos la de la siguiente manera
x = 2λ + 1
1 1 1
y = x = (2λ + 1) = λ +
2 2 2
Escogimos esta parametrizacion porque realmente no importa el resultado al final sera el que debe
de ser.
x=λ
1 1 1
y = x2 = (λ)2 = λ2
4 4 4
Para este ejemplo vamos a utilizar la métrica euclidiana y derivemos las parametrizaciones
dx dy
=2 =1 Para 1
dλ dλ
dx dy 1
=1 = λ Para 2
dλ dλ 2
Recordando la ecuación de longitud de arco
Zλ2 r
dxi dxj
S= gij · dλ
dλ dλ
λ1
Remplazando tenemos los términos y tomando que gxx = 1, gyy = 1 y gxy = gyx = 0 para el la
parametrizacion 1
s
Zλ2 1 0 0 1
q dx dx q dx dy q dy dx q dy dy
S= gxx + gxy + gyx + gyy · dλ
dλ dλ dλ dλ dλ dλ dλ dλ
λ1
Zλ2 p Zλ2 √ √ Z
λ2
= (2)(2) + (1)(1) · dλ = 5 · dλ = 5 dλ
λ1 λ1 λ1
IPN ESFM
Para que la parametrizacion cumpla y calcular los limites

1
2λ + 1 = 0 → λ=− ; x=0
2
1
2λ + 1 = 2 → λ= ; x=2
2
1
2
√ Z √ 1 √
S= 5 dλ = 5 λ]−
2
1 = 5
2
− 21
para el la parametrizacion 2 tenemos
Zλ2
s 2 2
dx dy
S= gxx + gyy · dλ
dλ dλ
λ1
Zλ2
s 2
λ
= (1)(1)2 + (1) · dλ
2
λ1
Zλ2 r
λ2
= 1+ · dλ
4
λ1
Para que la parametrizacion cumpla y calcular los limites

1
λ=0 → λ=− ; x=0
2
1
λ=2 → λ= ; x=2
2
Z2 r
λ2
S= 1+ · dλ
4
0

1 1p 2 2
−1 x
= λ + 4λ + 2Senh = 2,2956
2 2 2 0
Con estos resultados podemos ver que la longitud de la recta es menor que la de la curva, en el
plano cartesiano podemos asegurar que la distancia mas corta en dos puntos es una recta, entonces
el objetivo sera minimizar esta integral para cualquier métrica y ası́, encontrar la longitud mas
corta entre dos punto.
Como minimizamos una integral?
s 2
Z
dx dx dy
S= g11 (x, y) + g12 (x, y) + · · dλ
dλ dλ dλ
dx dy
Podemos ver que esta integral depende de x, y, dλ y dλ entonces en general tendremos que
Z λ2
dx dx
I(x(λ), y(λ)) = f (x, y, , )dλ
λ1 dλ dλ
Simplifiquemos este problema, para hacerlo mas claro, haciendo lo siguente

Z λ2
dx
I[x(λ)] = f (x, )dλ
λ1 dλ
Minimizar la integrar es encontrar la curva tal que la integral sea la minima.

Supongamos que tenemos la solución, y esta es una función x(λ), tal que
I[x(λ)] = mı́nimo
IPN ESFM
Esta funcion I[x] es un funcional.

Y si tenemos la una solución aproximada una tal que X(λ) = x(λ) + ϵ(λ), pero este pequeño
incremento es tal que (ϵ(λ))2 = 0, ∀λ ademas ϵ(a) = ϵ(b) = 0 Entonces con esta integral (o
funcional) que tenemos, si le sustituimos esta aproximación tendremos
I[X(λ)] = I[x(λ) + ϵ(λ)]
Entonces como estas dos curvas son muy parecidas, podremos hablar de la variacion
Z
dx dϵ dx
δI ≡ I[x + ϵ] − I[x] = f x + ϵ, + − f x, dλ = 0
dλ dλ dλ
Si recordamos el Teorema e Taylor:
∂ ∂
f (x + a, y + b) ≃ f (x, y) + fa + fb
∂x ∂y
Aplicando este teorema

dx dϵ dx ∂ ∂ dϵ
f x + ϵ, + ≃ f x, + f ·ϵ+ dx
f
dλ dλ dλ ∂x ∂ dλ dλ

dx dϵ dx ∂ ∂ dϵ
f x + ϵ, + − f x, ≃ f ·ϵ+ dx
f
dλ dλ dλ ∂x ∂ dλ dλ
y este resultado es lo que tenemos en nuestro funcional por lo que sustituyendo
Zλ2 " #
∂ ∂ dϵ
δI = f ·ϵ+ dx
f dλ = 0
∂x ∂ dλ dλ
λ1
Esto que hicimos es equivalente a lo que se hace en calculo 2 (si mal no recuerdo) cuando tenemos
una función f (x) y calculamos f ′ (x) = 0, ahora tenemos un funcional I y calculamos δI = 0
(variación de un funcional).
Z " #
∂ ∂ dϵ
δI = f ·ϵ+ dx
f dλ
∂x ∂ dλ dλ
Z Z
∂ ∂ dϵ
= f · ϵdλ + dx
f · dλ Integrando por partes
∂x ∂ dλ dλ
Hagamos un cambio de variable y calculemos la segunda integral

!
∂ d ∂
u= dx
f → du = dx
f dλ
∂ dλ dλ ∂ dλ
dϵ
v= dλ → v = ϵ
dλ
Ahora sustituyendo nuestro cambio de variable sobre la segunda integral
Z Z
∂ dϵ
dx
f · dλ = uv − vdu
∂ dλ dλ
" #λ2 Z !
∂ d ∂
= dx
f · ϵ − ϵ· dx
f dλ Por la condicion de ϵ
∂ dλ dλ ∂ dλ
λ1
" #λ2 Z !
∂ d ∂
= dx
f · ϵ

− ϵ· dx
f dλ
∂dλ
dλ ∂ dλ
λ1
IPN ESFM
con este calculo regresemos a la variacion del funcional

Z Z !
∂ d ∂
δI = ϵ · f dλ − ϵ · dx
f dλ
∂x dλ ∂ dλ
Z " !#
∂ d ∂
= f− dx
f ϵdλ
∂x dλ ∂ dλ
Si queremos que esta variacion δI = 0 basta conque

!
∂ d ∂
f− dx
f =0
∂x dλ ∂ dλ
Esta es la Ecuación de EULER - LAGRANGE
Regresemos al problema general que tiene 2 variables

Z λ2
dx dx
I(x(λ), y(λ)) = f (x, y, , )dλ
λ1 dλ dλ
Nuestra funcion con el incremento en 2 variables esta dada

dx dϵx dy dϵy
f x + ϵx , y + ϵy , + , +
dλ dλ dλ dλ
Realizando los mismos calculos tendremos que
 !!    
Z
∂ d ∂ ∂ d  ∂
δI =  f− dx
f ϵx +  f − f  ϵy  dλ
∂x dλ ∂ dλ ∂y dλ ∂ dy
dλ
Para que este funcional de 2 variables de cero, cada una de las Ecuaciones de Euler-Lagrange debe
ser cero por separado.
Esto nos da una pista de la forma fácil de generalizar estos funcionales, usando la notación de
Einstein, si tenemos la integral
Z λ2 i

i i dx
I[x (λ)] = f x, dλ
λ1 dλ
Las ecuaciones de Euler-Lagrange estaran dadas por

!
∂ i d ∂
x = i f
∂f dλ ∂ dx
dλ
Con estos resultados, resolvamos el mismo ejemplo que planteamos al inicio.
Elemplo
Zλ2
s 2 2
dx dy
S= + dλ
dλ dλ
λ1
r 2
dx 2 dy

Para este ejemplo tenemos que f = dλ + dλ Y por la ecuaciones de Euler-Lagrange
tenemos que
 
!
dx
∂ d ∂ d  2 dλ 
f =0 → dx
f =  r =0
∂x dλ ∂ dλ dλ  dx
2

dy
2 
2 dλ + dλ
IPN ESFM
Analogamente para y
 
!
dy
∂ d ∂ d  2 dλ 
f =0 → dy
f =  r
2  = 0

∂y dλ ∂ dλ dλ  dx 2

dy
2 dλ + dλ
De aqui tenemos que

dx
dλ
r 2 = C1
dx 2 dy

dλ + dλ
dy
dλ
r 2 = C2
dx 2 dy

dλ + dλ
Dividiendo estas dos ecuaciones

dy
C2 dλ
= dx
C1 dλ
De donde
dy C2 dx
= → dy = C3 dx
dλ C1 dλ
Resolviendo esta ecuacion diferencial nos da que
y = C3 x + C4 ⇒ y = ax + b
Para calcular las constantes C3 y C4 solo debemos obligar a la ecuación a que pase por los puntos
que nos indican las condiciones iniciales.
Con este método de Euler- Lagrange obtenemos una recta, esto pasa porque estamos en un espacio
euclideo, si estuviéramos en una esfera Nos darı́a las GEODESICAS.
Calcular geodesicas, es encontrar las curvas xi (λ) que minimizan la longitud en cualquier espacio,
plano/curvo.
Por lo tanto, resolver este problema es saber como se va a mover una particula en el Espacio-
Tiempo, cuando esta curvado.
Con las ecuacion de Euler - Lagrange podremos calcular las geodesicas, es decir minimizar la
distancia entre dos puntos, en cualquier espacio. Entonces, simplifiquemos la notación
dxi

∂ d ∂ i
f = f con u =
∂xi dλ ∂ui dλ
Tendremos a f como q
f= gjk uj uk
y su derivada es
∂ 1 ∂ gjk uj uk
f = ·
∂ui ∂ui
p
2 gjk uj uk
realicemos un ejemplo, en un espacio de dimensión 2

∂ gjk uj uk ∂ g11 (u1 )2 + g12 u1 u2 + g21 u2 u1 + g22 (u2 )2
=
∂u1 ∂u1
= g11 2u + g12 u + g21 u2
1 2
= 2(g11 u1 + g12 u2 )
= 2(g1k uk )
IPN ESFM
Entonces podemos deducir que

∂ gjk uj uk
= 2gik uk
∂ui
Para generalizar esta derivada tenemos que

∂ gjk uj uk
∂ui
Tenemos que (ABC)′ = A′ BC + AB ′ C + ABC ′ entonces
0
∂ gjk uj uk ∂ >
∂ ∂
= gjk uj uk + gjk i uj uk + gjk uj i uk
∂ui ∂ui
∂u ∂u
= gjk δij uk + gjk uj δik
= (gjk δij )uk + (gjk δik )uj
= gik uk + gji uj
= gij uj + gik uk = gij uj + gij uj
= 2gij uj
ahora realicemos la derivada con respecto de λ

∂ d ∂
f= f
∂xi dλ ∂ui
d − 1
= gjk uj uk 2 · gip up
dλ
d − 1 − 1 dgip p − 1 dup
= gjk uj uk 2 gip up + gjk uj uk 2 u + gjk uj uk 2 gip
dλ dλ dλ
Ahora anexemos un calculo extra para reducir el calculo anterior
d ∂ dx ∂ dy
f (x, y) = f + f
dλ ∂x dλ ∂y dλ
entonces
d ∂ dxk ∂
f (xk ) = f = f uk
dλ ∂xk dλ ∂xk
por lo que tendremos para el tensor métrico y definiendo la notación.
d ∂ ∂
gip = gip uk = (∂k gip )uk f ≡ ∂k f
dλ ∂xk ∂xk
Entonces
dup

d ∂ d − 1 − 1 d − 1
i
f = gjk uj uk 2 gip up + gjk uj uk 2 gip up + gjk uj uk 2 gip
dλ ∂u dλ dλ dλ
k
dup

d − 1 − ∂
1 dx p 1
j k −2
= gjk uj uk 2 gip up + gjk uj uk 2 gip u + gjk u u gip
dλ ∂xk dλ dλ
p

d − 1 − 1 − 1 du
= gjk uj uk 2 gip up + gjk uj uk 2 ∂k gip uk up + gjk uj uk 2 gip
dλ dλ
ahora calculemos el lado izquierdo de la ecuación de Euler.
p
∂ g ui uk
∂ jk
f=
∂xi ∂xi
− 21
∂ gjk uj uk
q
1 i k
= gjk u u
2 ∂xi
q − 21 " j k
#
1 ∂(gjk ) ∂(u ) ∂(u )
= gjk ui uk uj uk + gjk uk + gjk uj
2 ∂xi ∂xi ∂xi
1 − 1
= gjk ui uk 2 · (∂i gkj )uj uk
2
IPN ESFM
Entonces igualando los dos términos

∂ d ∂
f= f
∂xi dλ ∂ui
dup

1 1
i k −2 j k d 1
j k −2
− 1 − 1
gjk u u · (∂i gkj )u u = gjk u u gip up + gjk uj uk 2 ∂k gip uk up + gjk uj uk 2 gip
2 dλ dλ
p

1 i k − 1
j k d j k − 1
p j k − 1
k p j k − 1 du
()1/2
2
2
2
2
gjk u u · (∂i gkj )u u = gjk u u gip u + gjk u u ∂k gip u u + gjk u u gip
2 dλ dλ
dup

1 j k
1 d
j k 2
1
j k −2
(∂i gjk )u u = gjk u u gjk u u gip up + ∂k gip uk up + gip
2 dλ dλ
p

du 1 1 d − 1
= (∂i gjk )uj uk − gjk uj uk 2 gjk uj uk 2 gip up − ∂k gip uk up

gip
dλ 2 dλ
dup

mi 1 j k
1 d
j k 2
1
j k −2 p k p
g gip = (∂i gjk )u u − gjk u u gjk u u gip u − ∂k gip u u
dλ 2 dλ
p

m du 1 mi j k j k 2
1 d j k −2
1 m p
δp = g (∂i gjk )u u − gjk u u gjk u u δp u − g mi (∂k gip )uk up
dλ 2 dλ
dum

1 mi j k j k 2
1 d j k −2
1 m
= g (∂i gjk )u u − gjk u u gjk u u u − g mi (∂k gip )uk up
dλ 2 dλ
dum

1 d
j k 2
1 m
j k −2 1
− gjk u u gjk u u u = + g mi (∂k gip )uk up − g mi (∂i gjk )uj uk
dλ dλ 2
m

1 d − 1 du 1
− gjk uj uk 2 gjk uj uk 2 um = + g mi 2(∂k gip )uk up − (∂i gjk )uj uk

si j → p
dλ dλ 2
dum

1 d − 1 1
− gjk uj uk 2 gjk uj uk 2 um = + g mi (2(∂k gip ) − (∂i gkp )) uk up
dλ dλ 2
Desarrollemos el siguiente termino y reordenemos
2(∂k gip )uk up = ∂k gip uk up + ∂k gip uk up con k → p
k p k p
= ∂k gip u u + ∂p gik u u
= (∂k gip + ∂p gik ) uk up
Entonces remplazamos
dum

1 d
j k 2
1 m
j k −2 1
− gjk u u gjk u u u = + g mi (∂k gip + ∂p gik − ∂i gkp ) uk up
dλ dλ 2
Ahora para visualizar un resultado necesario, hagamos un ejemplo
dx
x = x(λ), →
dλ
dy dx
y = x2 (λ), → = 2x(λ) ·
dλ dλ
con la métrica EUCLIDIANA
) s 2 2
ux = vy
Z
dx dy
+ dλ
uy = −vx dλ dλ
s 2 2 s 2 2
Z Z
dx dy dx dx
+ dλ = + 4x2 (λ)
dλ
dλ dλ dλ dλ
Z p Z
dx
= 1 + 4x2 (λ) · dλ = 1dλ
dλ
Elijamos la parametrizacion tal que
dx 1
=p
dλ 1 + 4x2 (λ)
IPN ESFM
entonces tendremos que

j k
12 d 1 m
j k −2
1 d
− gjk u u gjk u u u = − gjk uj uk 2 (1)um = 0
dλ dλ
dum

1 mi k p
1 d
j k 2
1 m
j k −2
+ g (∂k gip + ∂p gik − ∂i gkp ) u u = − gjk u u gjk u u u
dλ 2 dλ
dum 1
+ g mi (∂k gip + ∂p gik − ∂i gkp ) uk up = 0
dλ 2
Ejemplo: Geodésica
Sea x1 = x, x2 = y y
− 1+1 b
b
! !
−1 − y 0 −1
0
g= 0 1 = gij y su inversa g = y
b
= g ij
1+ yb 0 1+ y
Escribir la ecuación de la geodésica
dum 1
+ g mi (∂k gip + ∂p gik − ∂i gkp ) uk up = 0
dS 2
dx1 dx
Sol.- Para el primer la primera variable u1 = dS = dS
dum 1
+ g mi (∂k gip + ∂p gik − ∂i gkp ) uk up
dS 2
du1 1 1i
= + g (∂k gip + ∂p gik − ∂i gkp ) uk up
dS 2"
0
! #
du1
1 11 0
k p 12 k p
= + g ∂k g1p + ∂p g1k − ∂1
gkp *
u u +
g >
(∂k g2p + ∂p g2k − ∂2 gkp ) u u
dS 2
du1 1
= + g 11 (∂k g1p + ∂p g1k ) uk up
dS 2
du1 *0 + ∂ g

1
0 2 p

= + g 11 ∂1g1p

p 11 u 1 p
u + ∂ g
2 1p + ∂pg12
*
u u
dS 2
du1 1 : 0u1 u1 + ∂ g · u2 u1 + ∂ g · u1 u2 + ∂ g * 0u2 u2
h i
= + g 11 ∂1g11

2 11 2 11 212
dS 2
du1 1 11
∂2 g11 · u2 u1 + ∂2 g11 · u1 u2

= + g
dS 2
du1
= + g 11 ∂2 g11 · u1 u2
dS
derivando

∂ −1 − yb b
∂2 g11 = =
∂y y2
IPN ESFM
entonces
dum 1
dS 2
du1
= + g 11 ∂2 g11 · u1 u2
dS !
du1 1 b
= + − u1 u2
dS 1 + yb y2
du1 b
= − u1 u2
dS b
1 + y y2
d dx b dx dy
= − ·
dS dS 1 + yb y 2 dS dS
b
=ẍ − · ẋẏ
b
1+ y y2
Por lo que la primera ecuación es
b
ẍ − · ẋẏ = 0
b
1+ y y2
dx2 dy
Ahora para u2 = dS = dS
dum 1
dS 2
du2 1 2i
= + g (∂k gip + ∂p gik − ∂i gkp ) uk up
dS 2"
0
#
2
du 1 21 k p 22
> (∂ g + ∂ g − ∂ g ) u u + g (∂ g + ∂ g − ∂ g ) u u k p
= + g k 1p p 1k 1 kp k 2p p 2k 2 kp
dS 2
du2 1
= + g 22 (∂k g2p + ∂p g2k − ∂2 gkp ) uk up
dS 2
du2 1 22
0 0
1 p 2 p
= + g ∂1g2p +
*
∂pg21 − ∂2 g1p u u + (∂2 g2p + ∂p g22 − ∂2 g2p ) u u
*

dS 2
du2 1
+ g 22 −∂2 g1p · u1 up + (∂2 g2p + ∂p g22 − ∂2 g2p ) · u2 up

=
dS 2
du2 1 22 h ((((
( i
= + g −∂2 g11 · u1 u1 + (∂( 2 g21 +(
( ∂1 g( 2 1
22 − ∂2 g21 ) · u u + −∂2
g 1 2 2 2
12 · u u + (∂2 g22 + ∂2 g22 − ∂2 g22 ) · u u
dS 2 ( (
du2 1 22
−∂2 g11 · u1 u1 + ∂2 g22 · u2 u2

= + g
dS 2
derivando

1
∂ 1+ yb b
∂2 g22 = = 2
∂y

y 2 1 + yb
IPN ESFM
entonces
dum 1
dS 2
du2 1 22
−∂2 g11 · u1 u1 + ∂2 g22 · u2 u2

= + g
dS 2
du2 1 b b 1 b b
= − (1 + ) 2 · (u1 )2 + (1 + ) 2 2
2 · (u )
dS 2 y y 2 y 2 b
y 1+ y
du2 b b b
= − 2 (1 + ) · (u1 )2 + · (u2 )2
dS 2y y 2 1 + yb y 2
du2 b b b
= − 2 (1 + ) · (u1 )2 + · (u2 )2
dS 2y y 2 1 + yb y 2
b b b
=ÿ − (1 + ) · (ẋ)2 + · (ẏ)2
2y 2 y b
2 1 + y y2
Por lo que la segunda ecuación es

b b b
ÿ − 2 1 + · (ẋ)2 + · (ẏ)2 = 0
2y y b
2 1 + y y2
Estas dos ecuaciones son las geodésicas de la métrica que propusimos.

b
ẍ − · ẋẏ = 0
1 + yb y 2

b b b
ÿ − 2 1 + · (ẋ)2 + · (ẏ)2 = 0
2y y b
2 1 + y y2
Transporte paralelo
Para hablar de transporte paralelo vamos a trabajar los vectores en R3 y utilizaremos las
coordenadas esféricas
x = R Sen (θ) Cos (ϕ)

y = R Sen (θ) Sen (ϕ)
z = R Cos (θ)
vamos a recordar que

∂ i
êi = x êj
∂xj
para la primera variable x1 = R
∂ ∂ ∂
êR = xêx + yêy + zêz
∂R ∂R ∂R
= Sen (θ) Cos (ϕ) êx + Sen (θ) Sen (ϕ) êy + Cos (θ) êz
para la segunda variable x2 = θ

∂ ∂ ∂
êθ = xêx + yêy + zêz
∂θ ∂θ ∂θ
= RCos (θ) Cos (ϕ) êx + RSen (θ) Sen (ϕ) êy − RSen (θ) êz
IPN ESFM
para la segunda variable x3 = ϕ

∂ ∂ ∂
êϕ = xêx + yêy + zêz
∂ϕ ∂ϕ ∂ϕ
= −RSen (θ) Sen (ϕ) êx + RSen (θ) Cos (ϕ) êy
Como podemos recordar estas bases no son ortonormales.

Ahora planteamos la siguiente cuestión, como varia êR , êθ , êϕ , cuando cambian las coordenadas del
punto, es decir como varı́an los vectores base a medida que va cambiando el punto donde estamos
mirando, esto nos va a ayudar a definir la derivada covariante.
Es fácil ver que, para el cambio en coordenadas cartecianas es 0, es decir
∂ ∂ ∂
êx = 0; êx = 0; êx = 0
∂x ∂y ∂z
Esto nos dice que estos vectores son únicos, que para cualquier componente el vector base es único,
en cambio para los vectores base de las coordenadas esféricas si hay cambios en los vectores base,
es decir, para êR
∂
êR = 0
∂R
∂
êR = Cos (θ) Cos (ϕ) êx + Cos (θ) Sen (ϕ) êy − Sen (θ) êz
∂θ
∂
êR = −Sen (θ) Sen (ϕ) êx + Sen (θ) Cos (ϕ) êy
∂ϕ
construyamos la metrica de esta nueva base para conocer los vectores base duales.
 
êR · êR êR · êθ êR · êϕ
g =  êθ · êR êθ · êθ êθ · êϕ 
êϕ · êR êϕ · êθ êϕ · êϕ
Calculemos las componentes de la métrica
gRR = êR · êR = 1
gRθ = êR · êθ = 0
gθR = êθ · êR = 0
gRϕ = êR · êϕ = 0
gϕR = êϕ · êR = 0
gθϕ = êθ · êϕ = 0
gϕθ = êϕ · êθ = 0
Estos cálculos se pueden realizar, ahora

gθθ = θ · θ
= R2 Cos2 (θ)Cos2 (ϕ) + R2 Cos2 (θ)Sen2 (ϕ) + R2 Sen2 (θ)
= R2 Cos2 (θ) + R2 Sen2 (θ) = R2
gϕϕ = ϕ · ϕ
= R2 Sen2 (θ)Sen2 (ϕ) + R2 Sen2 (θ)Cos2 (ϕ) + R2 Sen2 (θ)
= R2 Sen2 (θ)
Entonces tenemos que la métrica es
 
1 0 0
gi′ j ′ = 0 R2 0 
0 0 R2 Sen2 (θ)
IPN ESFM
Ahora calculemos êR , recordando que êi = g ij êj de las secciones anteriores y tendremos que
êR = g RR êR + g Rθ êθ + g Rϕ êϕ
Necesitaremos calcular la inversa de la métrica, que es un calculo relativamente fácil recordando

las propiedades de las matrices diagonales, y esta es
 
1 0 0
′ ′ 1
g i j = 0 R2 0 
1
0 0 R2 Sen 2 (θ)
Luego completando el calculo podremos tener el resultado
êR = g RR êR + g Rθ êθ + g Rϕ êϕ

1 0 0
=g RR g Rθ
> ê + g Rϕ
>

êR +
θ
> ê

ϕ
= êR
Análogamente para los otros vectores
êθ = g θR êR + g θθ êθ + g θϕ êϕ

0 0
=g θR
> ê + g θθ ê + g θϕ

R θ
> ê

ϕ
1
= 2 êθ
R
Y para el siguiente
êϕ = g ϕR êR + g ϕθ êθ + g ϕϕ êϕ

0 0
=g ϕR
> ê + g ϕθ

R êθ + g êϕ
>
ϕϕ
1
= 2 êϕ
R Sen2 (θ)
Ahora tenemos que

∂ êR
=0
∂R
Luego respecto de θ
∂ êR
= Cos (θ) Cos (ϕ) êx + Cos (θ) Sen (ϕ) êy − Sen (θ) êz
∂θ
1
= R (Cos (θ) Cos (ϕ) êx + Cos (θ) Sen (ϕ) êy − Sen (θ) êz )
R
1
= êθ
R
Luego respecto a ϕ
∂ êR
= −Sen (θ) Sen (ϕ) êx + Sen (θ) Cos (ϕ) êy
∂ϕ
1
= R (−Sen (θ) Sen (ϕ) êx + Sen (θ) Cos (ϕ) êy )
R
1
= êϕ
R
Vamos a introducir la nomenclatura que se utiliza en Relatividad
∂ êR ϕ
≡ ΓR θ
RR êR + ΓRR êθ + ΓRR êϕ
∂R
IPN ESFM
Esto lo podrı́amos tomar como una definición, que nos dice que cuando derivamos un vector unitario
con respecto a las coordenadas de la base, esto en general va a ser una combinación lineal de los
vectores de la base. Ahora recordando que
∂ êR
=0
∂R
0 0 0
ϕ
ΓR θ
>
>
>

=
RR êR + Γ
RR êθ + RR êϕ
Γ
Como todos los sı́mbolos son cero podemos deducir que
ΓiRR = 0
Luego para las derivadas respecto a las otras variables tendremos

∂ êR 1
= êθ
∂θ R
ϕ
= ΓR θ
Rθ êR + ΓRθ êθ + ΓRθ êϕ
0 0
=ΓR>
θ
Rθ êR + ΓRθ êθ + ΓϕRθ
> ê

ϕ

1
ΓθRθ =
R
y
∂ êR 1
= êϕ
∂ϕ R
ϕ
= ΓR θ
Rϕ êR + ΓRϕ êθ + ΓRϕ êϕ
0 0
=ΓR > ê + Γϕ ê
ΓθRϕ
Rϕ êR +
>

θ Rϕ ϕ

1
ΓϕRϕ =
R
Estos son los resultados de derivar êR respecto a las componentes de la base, entonces para el
vector êθ , calculando.
∂ êθ ϕ
= ΓR θ
θR êR + ΓθR êθ + ΓθR êϕ
∂R
Como
êθ = R Cos (θ) Cos (ϕ) êx + R Cos (θ) Sen (ϕ) êy − R Sen (θ) êz
Derivando respecto de R
∂ êθ
= Cos (θ) Cos (ϕ) êx + Cos (θ) Sen (ϕ) êy − Sen (θ) êz
∂R
ϕ
= ΓR θ
θR êR + ΓθR êθ + ΓθR êϕ
⃗ con este resultado

Para calcular los Gamma utilizaremos un resultado que ya conocemos A1 = ê1 · A
podemos tener que
∂ êθ
ΓR R
θR = ê ·
∂R
∂ êθ
= (êR ) ·
∂R
= (êR ) · (Cos (θ) Cos (ϕ) êx + Cos (θ) Sen (ϕ) êy − Sen (θ) êz )
= Sen (θ) Cos (θ) Cos2 (ϕ) + Cos (θ) Sen (θ) Sen2 (ϕ) − Cos (θ) Sen (θ)

= Cos (θ) Sen (θ) Cos2 (ϕ) + Sen2 (ϕ) − 1 = 0

IPN ESFM
Luego
∂ êθ
ΓθθR = êθ ·
∂R
1 ∂ êθ
= êθ ·
R2 ∂R

1 ∂ êθ
= 2 (R Cos (θ) Cos (ϕ) êx + R Cos (θ) Sen (ϕ) êy − R Sen (θ) êz ) ·
R ∂R
1
= 2 R Cos2 (θ) Cos2 (ϕ) + R Cos2 (θ) Sen2 (ϕ) + R Sen2 (θ)

R
:1
 
R 
= 2 Cos2 (θ) (Cos 2
Sen2 (ϕ)) + Sen2 (θ)

(ϕ)
+
R
1 1
Cos2 (θ) + Sen2 (θ) =

=
R R
Y por ultimo
∂ êθ
ΓϕθR = êϕ ·
∂R
1 ∂ êθ
= êϕ ·
R2 Sen2 (θ) ∂R

1 ∂ êθ
= 2 (−R Sen (θ) Sen (ϕ) êx + R Sen (θ) Cos (ϕ) êy ) ·
R Sen2 (θ) ∂R
1
= 2 ( −R Cos (θ) Sen (θ) Cos (ϕ) Sen (ϕ) + R Cos (θ) Sen (θ) Cos (ϕ) Sen (ϕ)) = 0
R Sen2 (θ)
Seguiremos calculando la variacion de los vectores unitarios respecto a sus componentes, como
el calculo es un poco tedioso, solo realizaremos pocos y al final solo escribiremos los resultados
dejando los cálculos como ejercicios del capitulo. entonces
∂ êϕ ϕ
= ΓR θ
ϕR êR + ΓϕR êθ + ΓϕR êϕ
∂R
= −Sen (θ) Sen (ϕ) êx + Sen (θ) Cos (ϕ) êy
De forma análoga al calculo anterior
∂ êϕ
ΓR R
ϕR = ê ·
∂R
∂ êϕ
= (êR ) ·
∂R
= (êR ) · (−Sen (θ) Sen (ϕ) êx + Sen (θ) Cos (ϕ) êy )
= −Sen2 (θ) Cos (ϕ) Sen (ϕ) + Sen2 (θ) Cos (ϕ) Sen (ϕ) = 0

Luego
∂ êϕ
ΓθϕR = êθ ·
∂R
1 ∂ êϕ
= ê θ ·
R2 ∂R

1 ∂ êϕ
= 2 (R Cos (θ) Cos (ϕ) êx + R Cos (θ) Sen (ϕ) êy − R Sen (θ) êz ) ·
R ∂R
1
= 2 (−R Cos (θ) Cos (ϕ) Sen (θ) Sen (ϕ) + R Cos (θ) Cos (ϕ) Sen (θ) Sen (ϕ)) = 0
R
IPN ESFM
Y por ultimo
∂ êϕ
ΓϕϕR = êϕ ·
∂R
1 ∂ êϕ
= êϕ ·
R2 Sen2 (θ) ∂R

1 ∂ êϕ
= 2 (−R Sen (θ) Sen (ϕ) êx + R Sen (θ) Cos (ϕ) êy ) ·
R Sen2 (θ) ∂R
1 1
R Sen2 (θ) Sen2 (ϕ) + R Sen2 (θ) Cos2 (ϕ) =

= 2
R Sen2 (θ) R
Con esto tenemos todos los resultados de la variación de los vectores base respecto a R, ası́ como
sus respectivos valores de Gamma, ası́ que para los siguientes cálculos vamos a reducir el calculo
al mı́nimo. Ahora vamos a calcular la variación de los vectores base respecto de θ por lo que
∂ êR ϕ
= ΓR θ
Rθ êR + ΓRθ êθ + ΓRθ êϕ
∂θ
1
= Cos (θ) Cos (ϕ) êx + Cos (θ) Sen (ϕ) êy − Sen (θ) êz = êθ
R
Luego cada sı́mbolo se calcula como
∂ êR 1
ΓR R
Rθ = ê · = êR · êθ = 0
∂θ R
∂ êR 1 1
ΓθRθ = êθ · = êθ · êθ =
∂θ R R
∂ êR 1
ΓϕRθ = êϕ · = êϕ · êθ = 0
∂θ R
Para la segunda parcial
∂ êθ ϕ
= ΓR θ
θθ êR + Γθθ êθ + Γθθ êϕ
∂θ
= −R Sen (θ) Cos (ϕ) êx + R Sen (θ) Sen (ϕ) êy − R Cos (θ) êz = −R êR
Entonces los valores de las constantes Gamma son
∂ êθ
ΓR R
θθ = ê · = êR · −R êR = −R
∂θ
∂ êθ
Γθθθ = êθ · = êθ · −R êR = 0
∂θ
∂ êθ
Γϕθθ = êϕ · = êϕ · −R êR = 0
∂θ
Para la ultima parcial
∂ êϕ ϕ
= ΓR θ
ϕθ êR + Γϕθ êθ + Γϕθ êϕ
∂θ
= −R Cos (θ) Sen (ϕ) êx + R Cos (θ) Cos (ϕ) êy
IPN ESFM
Y sus constantes Gamma quedan como

∂ êϕ ∂ êϕ
ΓR R
ϕθ = ê · = êR ·
∂θ ∂θ
= −R Cos (θ) Sen (θ) Cos (ϕ) Sen (ϕ) + R Cos (θ) Sen (θ) Cos (ϕ) Sen (ϕ) = 0
∂ êϕ 1 ∂ êϕ
Γθϕθ = êθ · = 2 êθ ·
∂θ R ∂θ
1
= 2 −R Cos (θ) Cos (ϕ) Sen (ϕ) + R Cos2 (θ) Cos (ϕ) Sen (ϕ) = 0
2

R
∂ êϕ 1 ∂ êϕ
Γϕϕθ = êϕ · = 2 êϕ ·
∂θ R Sen2 (θ) ∂θ
1 Sen (θ)
R2 Cos (θ) Sen (θ) Sen2 (ϕ) + R2 Cos (θ) Sen (θ) Cos2 (ϕ) =

= 2 2
R Sen (θ) Cos (θ)
Ahora derivando respecto de ϕ derivamos todos los vectores base, análogamente como los anteriores
∂ êR ϕ
= ΓR θ
Rϕ êR + ΓRϕ êθ + ΓRϕ êϕ
∂ϕ
1
= −Sen (θ) Sen (ϕ) êx + Sen (θ) Cos (ϕ) êy = êϕ
R
Luego cada sı́mbolo se calcula como
∂ êR 1
ΓR R
Rϕ = ê · = êR · êϕ = 0
∂ϕ R
∂ êR 1
ΓθRϕ = êθ · = êθ · êϕ = 0
∂ϕ R
∂ ê R 1 1
ΓϕRϕ = êϕ · = êϕ · êϕ =
∂ϕ R R
Para la segunda parcial
∂ êθ ϕ
= ΓR θ
θϕ êR + Γθϕ êθ + Γθϕ êϕ
∂ϕ
Entonces los valores de las constantes Gamma son

∂ êθ
ΓR R
θϕ = ê ·
∂ϕ
= −R Cos (θ) Sen (θ) Cos (ϕ) Sen (ϕ) + R Cos (θ) Sen (θ) Cos (ϕ) Sen (ϕ) = 0
∂ êθ
Γθθϕ = êθ ·
∂ϕ
1
= 2 −R2 Cos2 (θ) Cos (ϕ) Sen (ϕ) + R2 Cos2 (θ) Cos (ϕ) Sen (ϕ) = 0

R
∂ êθ
Γϕθϕ = êϕ ·
∂ϕ
1 Cos (θ)
R Cos (θ) Sen (θ) Sen2 (ϕ) + R2 Cos (θ) Sen (θ) Cos2 (ϕ) =
2
= 2 2
R Sen (θ) Sen (θ)
Para la ultima parcial

∂ êϕ ϕ
= ΓR θ
ϕϕ êR + Γϕϕ êθ + Γϕϕ êϕ
∂ϕ
IPN ESFM
Y sus constantes Gamma quedan como

∂ êϕ
ΓR R
ϕϕ = ê ·
∂ϕ
= −RSen2 (θ) Cos2 (ϕ) − RSen2 (θ) Sen2 (ϕ) = −R Sen2 (θ)
∂ êϕ
Γθϕϕ = êθ ·
∂ϕ
1
= 2 −R2 Cos (θ) Sen (θ) Cos2 (ϕ) − R2 Cos (θ) Sen (θ) Sen2 (ϕ) = −Cos (θ) Sen (θ)

R
∂ êϕ
Γϕϕϕ = êϕ ·
∂ϕ
1
R Sen2 (θ) Sen (ϕ) Cos (ϕ) − R2 Sen2 (θ) Sen (ϕ) Cos (ϕ) = 0
2
= 2 2
R Sen (θ)
En resumen tenemos que
Cos (θ)
ΓR
θθ = −R, Γθϕϕ = −Sen (θ) Cos (θ) , Γϕϕθ = Γϕθϕ =
Sen (θ)
1
ΓR 2
ϕϕ = −RSen (θ) , ΓθθR = ΓθRθ = ΓϕϕR = ΓϕRϕ =
R
Con esto podemos notar que los indices de abajo nos indica que los Gamma son simétricos siempre
que el indice de arriba sea el mismo.
Ahora hablemos de paralelismo, imaginemos que tenemos un campo A ⃗ = 4êx , este sabemos que
mide 4 y siempre apunta en la dirección de z, pero supongamos que queremos trabajar en las
⃗ a la nueva base.
coordenadas (r, θ, ϕ), entonces vamos a convertir A
⃗ = A1 êR + A2 êθ + A3 êϕ = 4êz

A
Entonces sus componentes serán

⃗ = êR · 4êz
A1 = êR · A
= (Sen (θ) Cos (ϕ) êx + Sen (θ) Sen (ϕ) êy + Cos (θ) êz ) · 4êz = 4 Cos (θ)
A2 = êθ · 4êx
1 4
= 2 (R Cos (θ) Cos (ϕ) êx + R Cos (θ) Sen (ϕ) êy − R Sen (θ) êz ) · 4êz = − Sen (θ)
R R
A3 = êϕ · 4êx
1
= 2 (−R Sen (θ) Sen (ϕ) êx + R Sen (θ) Cos (ϕ) êy ) · 4êz = 0
R Sen2 (θ)
Entonces el vector en la nueva basa queda determinado como
⃗ = 4 Cos (θ) êR − 4 Sen (θ) êθ = AR êR + Aθ êθ + Aϕ êϕ = Ai êi
A
R
Ahora nos preguntamos cuanto es la derivada de este vector respecto a las componentes y utiliza-
remos la notación xi , entonces
⃗
∂A ∂ ∂Aj ∂ êj
j
= A êj = êj + Aj
∂xi ∂xi ∂xi ∂xi
∂Aj
= êj + Aj Γkji êk
∂xi
De este resultado podemos reordenar los indices de tal forma que
⃗
∂A ∂Ak
k
∂A

j k j k
= êk + A Γji êk = + A Γji êk
∂xi ∂xi ∂xi
IPN ESFM
Con esto tenemos la derivada de un vector respecto a las coordenadas i-esima de la base.
Cual es la interpretación de esto en un mundo euclidio, donde se ha trabajado toda la carrera?.
⃗ =
En el sistema euclidiano tenemos la base {êx , êy , êz }, luego para un vector cualquiera A
3êx + 2êy + 4êz y lo derivamos tendremos algo conocido.
⃗
∂A

∂3 ∂ êx

∂2 ∂ êy

∂4 ∂ êz

= êx + 3 + êy + 2 + êz + 4 =0
∂x ∂x ∂x ∂y ∂y ∂z ∂z
Esto lo tenemos claro porque los vectores unitarios de la base son constantes, y no tienes variación,
entonces cuando derivamos vectores en la base euclidiana, A ⃗ = Ax êx + Ay êy + Az êz tenemos que
⃗
∂A ∂Ax ∂Ay ∂Az
= êx + êy + êz
∂x ∂x ∂x ∂x
No tenemos que preocuparnos de derivar los vectores unitarios, porque sabemos que siempre van a
ser cero, pero para vectores que si cambian, entonces la derivada de esos vectores se expresa con las
constantes Gamma que completan la combinación lineal. Y podemos denotar esta derivada como
⃗
∂A
k
∂A

j k k

= i
+ A Γ ji êk ≡ ∇i A êk
∂x ∂x
Esto lo llamaremos derivada covariante y se difiere de la derivada parcial en que la esta es un caso
particular en la derivación.
Para el campo vectorial A ⃗ que tenemos de ejemplo, podemos notar que este deberı́a ser 0 porque
no cambia, aun que este en otra basa, calculemos la derivada en para comprobarlo como ejemplo.
Entonces tomemos el primer vector unitario con k = 1 y cambiando respecto a las coordenadas
esféricas tenemos
∂AR
∇ 1 A1 = ∇ R AR = + Aj Γ R
jR
∂R
∂AR R R
+ A ΓRR + Aθ ΓR ϕ R

= θR + A ΓϕR = 0
∂R
Como la componente AR = 4 Cos (θ) no depende de R la parcial es cero y revisando los Gamma
que ya calculamos tenemos que todos son cero, por lo que toda la derivada es cero.
∂Aθ
∇ 1 A2 = ∇ R Aθ = + Aj ΓθjR
∂R
∂Aθ R θ
+ A ΓRR + Aθ ΓθθR + Aϕ ΓθϕR

=
∂R
" 0 0
#
4 R θ 4 1 ϕ θ
ΓRR + − Sen (θ)
>

= 2 Sen (θ) + A +A
ΓϕR
>
=0
R R R
Para el siguiente caso

0
ϕ
∂A
∇ 1 A3 = ∇ R Aϕ = + Aj ΓϕjR
∂R
" 0 0 0
#
∂Aϕ R ϕ>
θ ϕ ϕ ϕ
= + A ΓRR + A ΓθR +
> A>
ΓϕR = 0
∂R
Completar que esto es cero

Si tenemos que la derivada covariente de un vector para cualquier base es cero, esto nos dice
que el vector es constante en el espacio.

Apuntes de Relatividad General

Cargado por

Copyright:

Formatos disponibles

Apuntes de Relatividad General

Cargado por

Información del documento

Descripción original:

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Apuntes de Relatividad General

Cargado por

Copyright:

Formatos disponibles

IPN ESFM

Instituto Politécnico Nacional

Escuela Superior de Fı́sica y Matemáticas

Notas del Curso

2abCos (θ) = ab (Cos (τ ) Cos (θ) + Sen (τ ) Sen (θ))

Considerando que ahı́ están las componentes del vector ⃗a y ⃗b

Ahora renombremos los productos de tal manera que

g11 = u1 · u1 , g12 = u1 · u2 + 3u2 , g21 = u2 · u1 , g22 = u2 · u2

A este reordenamiento lo llamaremos ”métrica”

Para una base euclidiana E la métrica esta dada por

Para la base S, la métrica es    

x = R Sen (θ) Cos (ϕ)

vamos a recordar que

para la segunda variable x2 = θ

para la segunda variable x3 = ϕ

Como podemos recordar estas bases no son ortonormales.

Calculemos las componentes de la métrica

gRR = êR · êR = 1

Estos cálculos se pueden realizar, ahora

(∆S)2 = (∆xê1 + ∆yê2 ) · (∆xê1 + ∆yê2 )

Ejemplo Sea y = x2 en el intervalo [0, 1], calcular la longitud de arco.

El parámetro va de 0 a 1, luego integrando

Usando la condición que definimos

ag12 g11 + bg12 g21 = g12

Para c y d multiplicamos (IV ) por g11 y a (III) por −g12 y sumamos

−cg12 g11 − cg12 g21 = 0

Luego si ordenamos las constantes de la manera que

De donde podemos tener que  

Para reducir la escritura de ecuaciones usaremos la notación de Einstein

Por lo tanto la basa dual esta definida como

Ejemplo Capitulo 2, 27:22

|gij | = (g11 g22 − g12 g21 ) = (2)(2) − (1)(1) = 3

ê2 = ê1 + 2ê2

Ahora recordemos que êi = g ij êj es fácil notar que

gij (êi ) = gij (g ij êj )

Entonces, cambiando los indices mudos

êi = gij êj

Con esta transformación de vectores para ê2 tendremos que

Esto es derivación implı́cita.

Ahora veamos un caso un poco mas complejo, sean las funciones

f = x2 (f, g) + y 2 (f, g) · · · · · · (1)

∂ ∂x2 (f, g) ∂y 2 (f, g)

Podemos notar que !

Calculando esta matriz inversa

Con estos resultados podemos proponer 2 bases que cumplen que

êa = 4êx + 5êy

escribamos la matriz asociada a este sistema

Entonces la transformación se da por

Ahora ya podemos hablar bien de lo que es una base coordenada

Ahora definiremos la NOTACIÓN OFICIAL del calculo tensorial

con su base asociada

Con la matriz inversa podemos despejar x y y (el calculo esta en corto)

Por ejemplo, sea un vector

SUBIR Y BAJAR INDICES (vectores covariantes y contravariantes) Para empezar,

Ejemplo: Utilizando las mismas coordenadas {a, b} y {x, y}

con esto el la métrica es  

êa = g aa êa + g ab êb = 2êa + 5êb

⃗ = Aa (2êa + 5êb ) + Ab (5êa + 17êb )

⃗ en la base cartesiana y la comparamos en la otra base

Para la base S, la métrica es

De donde podemos tener que

con esto el la métrica es