Universidade Estadual de Maringá

Centro de Ciências Exatas

Departamento de Física

Cálculo Tensorial e Relatividade Geral

Mateus Reinke Pelicer

Orientadora: Prof.a Dr.a Hatsumi Mukai

Maringá - PR
2016
 Universidade Estadual de Maringá
Centro de Ciências Exatas
Departamento de Física

Cálculo Tensorial e Relatividade Geral

Trabalho de Conclusão de Curso apresentado ao

Departamento de Física da Universidade Esta-
dual de Maringá, realizado sob orientação da
Prof.a Dr.a Hatsumi Mukai como requisito parcial
para a obtenção do título de Bacharel em Física.

Maringá - PR
2016
text

“lembrem de mim
como de um
que ouvia a chuva
como quem assiste missa
como quem hesita, mestiça,
entre a pressa e a preguiça. ”
Paulo Leminski
 Agradecimentos
Agradeço imensamente a meus pais, por todo carinho, paciência e cumpli-
cidade proporcionados durante esta jornada. Agradeço também a minha
namorada, por todo suporte, pela amizade e, claro, pelas bobagens de todo
dia, que facilitam o viver. Gostaria também de agradecer meus professores,
em especial a Hatsumi Mukai, Luis Carlos Malacarne e Breno Ferraz de Oli-
veira, pelo apoio dado e por toda dedicação e disposição ao me auxiliar no
decorrer do curso.

Obrigado.
Resumo

O presente trabalho trata sobre cálculo tensorial e Teoria da Relatividade

Geral (TRG), com maior foco no último. Primeiramente faz-se uma introdu-
ção aos tensores, abrangendo de álgebra a tensores de curvatura. Em seguida
são apresentados os experimentos mentais (gedankenexperiment) que levaram
à gênese da TRG, e são obtidas as equações de campo a partir da formulação
variacional. Por último é derivada a solução de Schwarzschild e analisadas
as consequências decorrentes desta. A importância do estudo da TRG se dá
em várias áreas, dentre estas a cosmologia moderna e a gravitação quântica.

Palavras-Chave: Tensores; Gravitação; Relatividade.

Abstract

The present work deals with tensor calculus and Theory of General Re-
lativity (TRG), focusing in the last one. First, an introduction to tensors is
made, covering from algebra to curvature tensors. Then, the thought expe-
riments (gedankenexperiment) that led to the genesis of TRG are presented,
followed by the field equations, wich are obtained using the variational for-
mulation. The Schwarzschild solution is derived and its consequences are
analyzed. The importance of studying TRG can be seen in many fields,
among wich modern cosmology and quantum gravitation.

Key-Words: Tensors; Gravitation; Relativity.

Sumário

Introdução 9

1 Tensores 11
1.1 A ideia por trás de tensores . . . . . . . . . . . . . . . . . . . 11
1.2 Álgebra Tensorial . . . . . . . . . . . . . . . . . . . . . . . . . 12
1.3 O tensor Métrico . . . . . . . . . . . . . . . . . . . . . . . . . 16
1.4 Cálculo Tensorial . . . . . . . . . . . . . . . . . . . . . . . . . 18
1.5 Geodésicas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
1.6 Curvatura . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
1.7 Tensores Relativos . . . . . . . . . . . . . . . . . . . . . . . . 24

2 Teoria da Relatividade Geral 26

2.1 Breve Histórico . . . . . . . . . . . . . . . . . . . . . . . . . . 26
2.2 O Princípio de Equivalência . . . . . . . . . . . . . . . . . . . 30
2.3 Equações de Campo . . . . . . . . . . . . . . . . . . . . . . . 34
2.4 Limite Newtoniano . . . . . . . . . . . . . . . . . . . . . . . . 40

3 Solução de Schwarzschild 43
3.1 A Métrica . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 43
3.2 Geodésicas e Órbitas . . . . . . . . . . . . . . . . . . . . . . . 47
3.3 O Periélio de Mercúrio . . . . . . . . . . . . . . . . . . . . . . 51
 3.4 Desvio Gravitacional da Luz . . . . . . . . . . . . . . . . . . . 55
3.5 Coordenadas de Lemaître . . . . . . . . . . . . . . . . . . . . 58

Considerações Finais 62

A Derivação da Equação da Geodésica. 67

B Expressão para o Divergente 70

C Simulação das Órbitas 72

Introdução

A Teoria da Relatividade Geral (TRG) foi elaborada no início do século

XX, e publicada na forma como é conhecida hoje em 1915, por Albert Eins-
tein e por David Hilbert, ambos trabalhando independentemente [1]. Nos
anos seguintes ocorreram estudos sistemáticos na área, envolvendo persona-
lidades como Karl Schwarzschild, Hermann Weyl e Georges Lemaître; como
resultado, previsões anteriormente inconcebíveis surgiram, como a existência
de buracos negros, a expansão do Universo e a existência de ondas gravita-
cionais [2].
Com a gênese da relatividade foi estabelecido um novo paradigma na Fí-
sica. Da Teoria da Relatividade Restrita (TRR), surgem dois conceitos de
extrema importância: o primeiro é a ideia de que espaço e tempo não são
entes independentes mas, na realidade, estão correlacionados; o segundo é a
geometrização do espaço-tempo. A gravitação, antes entendida como uma
força devido à interação a distância entre corpos, agora é interpretada como
uma deformação no espaço-tempo, no qual partículas em movimento seguem
geodésicas em espaços curvos, ou seja, “retas locais” [3, 4]. A deformação é
causada pelo tensor Energia-Momento, havendo uma relação entre massa e
energia de um corpo. Apenas destes fatos é possível compreender a impor-
tância da TRG tanto na Física como na epistemologia.
O desenvolvimento científico é incessante, e os avanços tardios que ocor-

9
Cálculo Tensorial e Relatividade Geral 10

reram em áreas como mecânica clássica e eletromagnetismo – que datam dos

séculos XVI e XIX, respectivamente – levam a ponderar o quanto ainda há a
ser desvendado no estudo da relatividade. Seguindo esta linha de raciocínio
o presente trabalho foi escrito a partir de uma densa revisão de literatura,
tendo como objetivo introduzir e analisar os conceitos matemáticos e con-
ceituais da TRG e suas consequências imediatas, visando possibilitar uma
futura pesquisa no campo.
A apresentação deste está exposta da seguinte forma: No Capítulo 1 é
dada uma introdução necessária da matemática necessária para o estudo da
TRG, focando seu uso como ferramenta matemática. Os tópicos aborda-
dos são álgebra, tensor métrico, cálculo tensorial, geodésicas, curvatura e
tensores relativos. No Capítulo 2 são abordados os experimentos mentais
que levam ao princípio de equivalência, seguido pelas equações de campo,
que são deduzidas a partir do princípio variacional. Enfim, é determinada
a constante gravitacional relativística a partir do limite clássico de campos
fracos. No Capítulo 3 é deduzida a solução de Schwarzschild, uma solução
das equações de campo para uma massa esfericamente simétrica, a partir da
qual previsões como o avanço do periélio de Mercúrio, o desvio gravitacional
da luz e a existência de buracos negros são analisados.
Capítulo 1

Tensores

Neste capítulo será desenvolvida uma introdução aos tensores, abrangendo

a sua álgebra, cálculo, geodésicas, tensores de curvatura e tensores relativos,
com o intuito de posteriormente aplicá-los no desenvolvimento da TRG.

1.1 A ideia por trás de tensores

O estudo de tensores é de extrema importância na Física, não somente
na relatividade geral, mas também em áreas como mecânica clássica e eletro-
magnetismo, em que são utilizados para descrever grandezas como potenciais,
stress e momento de inércia (Figura 1.1), e também para escrever algumas
equações de forma mais compacta e elegante, como as equações de Maxwell
e a de continuidade [7, 8]. São inúmeras as interpretações possiveis para os
tensores [9] mas, de acordo com Medrano [10]: “No que se diz respeito ao
significado físico dos tensores, pode-se dizer que estes encontram-se relacio-
nados as propriedades físicas, intrínsecas, do meio”. Tensores de ordem 0, 1
ou 2 podem ser representados por matrizes, já os de ordens maiores precisam
ser descritos a partir das componentes separadamente.

11
Cálculo Tensorial e Relatividade Geral 12

 
σ σ σ
 11 12 13 
σ =  σ21 σ22 σ23
 

 
σ31 σ32 σ33

Figura 1.1: Representação das componentes do tensor stress.

Fonte: Site: RockMechs. 1

Neste trabalho será usada uma abordagem clássica, onde se trabalha com
as componentes dos tensores. Apesar desta abordagem deixar de lado parte
da intuição geométrica, ela é mais elementar e a linguagem empregada no
tratamento remete à geometria [11], o que acaba tornando os conceitos de
mais fácil compreensão.

1.2 Álgebra Tensorial

Seja VN um espaço N-dimensional no qual estão definidos os sistemas
de coordenadas S e S̄, representados pelas coordenadas xµ e x̄µ , respecti-
vamente. As coordenadas do último podem ser escritas como funções das
coordenadas do primeiro:

x̄ν = x̄ν x0 , x1 , ..., xN −1 .

(1.2.1)
1
Disponível em: <[Link]
transformation-rotation/> Acesso em: agosto de 2016.
Cálculo Tensorial e Relatividade Geral 13

Derivando a equação (1.2.1), é obtida a relação da transformação para des-

locamentos infinitesimais

N −1
ν
X ∂ x̄ν µ
dx̄ = dx . (1.2.2)
µ=0
∂xµ

A partir daqui serão adotadas duas convenções:

• Convenção de Einstein: Para não ser necessária a escrita da somatória,

sempre que um índice se repetir em um lado da equação estará implicíto
que há uma soma sobre este índice.

• Intervalo dos índices gregos: Os índices gregos que se repetirem são

somados entre todo o intervalo entre 0 e N − 1.

Assim, a equação (1.2.2) pode ser reescrita como

∂ x̄ν µ
dx̄ν = dx . (1.2.3)
∂xµ

É convencionado que os índices das coordenadas sejam sobrescritos pois na

relação de transformação, do lado direito, a coordenada que está sendo dife-
renciada remete ao mesmo sistema do lado esquerdo. Por esta transformação
são definidas as componentes de um vetor contravariante

∂ x̄µ ν
V̄ µ = V . (1.2.4)
∂xν

Por outro lado, são ditas componentes de um vetor covariante aquelas que
se transformam por
∂xν
V̄µ = Vν . (1.2.5)
∂ x̄µ
Note que os índices subescritos seguem a mesma lógica empregada para os
vetores contravariantes.
Cálculo Tensorial e Relatividade Geral 14

Em um espaço bidimensional é possível dar significado geométrico a es-

tas componentes, como se observa na Figura 1.2. As componentes de um
vetor contravariante são paralelas aos eixos de coordenadas, enquanto as
componentes de um vetor covariante são perpendiculares aos respectivos ei-
xos. Quando se trabalha em espaços euclideanos, não há distinção entre as
componentes.

Figura 1.2: Representação das componentes de vetores.

Fonte: Autor (2016).

Um tensor de ordem maior é definido, de modo formal, a partir do pro-

duto tensorial entre espaços vetoriais [12]. Porém, uma definição equivalente
é aquela que se refere à tranformação de coordenadas de suas componentes
[10]. A ordem de um tensor pode ser representada por (m, n), sendo m a
ordem das componentes contravariantes e n das covariantes. Caso o tensor
tenha somente componentes covariantes (ou contravariantes), refere-se ape-
nas à ordem deste. Tensores contravariantes de segunda ordem são definidos
segundo a lei de transformação

∂ x̄α ∂ x̄β µν
T̄ αβ = T , (1.2.6)
∂xµ ∂xν

e tensores covariantes de segunda ordem por:

∂xµ ∂xν
T̄αβ = Tµν ; (1.2.7)
∂ x̄α ∂ x̄β
Cálculo Tensorial e Relatividade Geral 15

já um tensor que possui ambas componentes é denominado tensor misto:

∂xµ ∂ x̄β ν
T̄αβ = T . (1.2.8)
∂ x̄α ∂xν µ

Estas definições se estendem para tensores de ordem maior pela adição

de índices e termos de derivação respectivos aos índices

∂ x̄β1 ∂ x̄βm ∂xν1 ∂xνn µ1 µ2 ···µm

T̄ β1 β2 ···βm α1 α2 ···αn = · · · · · · T ν1 ν2 ...νn . (1.2.9)
∂xµ1 ∂xµm ∂ x̄α1 ∂ x̄αn

Esta é a definição mais geral a partir da lei de transformação para um tensor

de ordem (m, n). Nota-se que as relações de transformações são lineares e
homogêneas, ou seja, uma relação válida em um sistema de coordenadas é
válida em qualquer outro sistema [11]. É evidente pelas definições que vetores
são tensores de ordem 1 e, escalares são os de ordem 0, também chamados
de invariantes. Este últimos são grandezas cujo valor independe do sistema
de coordenadas adotado e são definidos a partir do produto interno

∂ x̄µ ∂xλ ω
V̄ µ Ūµ = V Uλ = δ λω V ω Uω = V λ Uλ , (1.2.10)
∂xω ∂ x̄µ

na qual δ ωλ é a delta de Kronecker 2 . O produto interno também é chamado

de contração visto que, se aplicado a tensores de ordem maior não surgem
invariantes, mas a ordem do tensor é subtraída em 2. A contração de um
tensor de tipo (3, 2) resulta em outro de tipo (2, 1):

T µν ω = T µλνλω = Aµλ B ν λω . (1.2.11)

Outra operação a ser definida é o produto tensorial (ou externo), já men-


2 ω 1 se λ = ω
δ λ =
0 se λ 6= ω
Cálculo Tensorial e Relatividade Geral 16

cionado. O produto tensorial entre um tensor de ordem (m, n) e outro de

ordem (r, s) gera um tensor de ordem (m + r, n + s), por exemplo, o produto
entre um vetor de ordem (1, 1) e outro (0, 2) gera um terceiro de ordem (1, 3),
como segue:
T αβµν = Aαβ Bµν . (1.2.12)

A combinação linear só é definida para tensores dados no mesmo ponto

do espaço e de mesma ordem. Sejam c1 e c2 escalares:

T αµν = c1 Aαµν + c2 B αµν . (1.2.13)

Tensores também podem apresentar propriedades de simetria. Diz-se que

um tensor é simétrico se a troca de dois indíces de mesmo tipo não altera
o valor da componente, esta definição é válida para tensores de qualquer
ordem. Para um tensor de tipo (0, 2), isto significa

Aµν = Aνµ , (1.2.14)

e denomina-se um tensor como antisimétrico se a troca de dois indíces do

mesmo tipo troca o sinal das componentes, novamente para um tensor cova-
riante de segunda ordem:
Aµν = −Aνµ . (1.2.15)

1.3 O tensor Métrico

Um espaço é dito métrico quando é possível definir a distância infinitesi-
mal ds entre dois pontos em função de um tensor covariante simétrico gµν ,
denominado tensor métrico (ou fundamental) que é uma função das coorde-
Cálculo Tensorial e Relatividade Geral 17

nadas [11]:
ds2 = gµν dxµ dxν . (1.3.1)

Este tipo de espaço possui a propriedade de conexão entre as componentes

contravariantes e covariantes do tensor, dada pela relação:

gµν U ν = Uµ . (1.3.2)

Um espaço em que não é possível definir tal conexão é denominado espaço

afim [13]. Pode-se definir um tensor contravariante g µν , denominado tensor
recíproco (ou conjugado) por

gµν g µλ = δ λν , (1.3.3)

tal que os vetores covariantes e contravariantes também estejam relacionados

por
g µν Uν = U µ . (1.3.4)

O tensor métrico determina algumas propriedades do espaço, como ho-

mogeneidade e isotropia, o primeiro ocorre se a métrica independe do ponto
escolhido como referência, já o segundo, quando existe um vetor ~v 6= 0 tal
que ds2 = 0. Uma interpretação física seria, se a métrica é homogênea, dado
um tempo fixo, o espaço é o mesmo de qualquer ponto de referência e, se
isotrópico, o espaço é o mesmo visto de qualquer direção. Na relatividade ge-
ral está diretamente relacionado ao potencial gravitacional. Um espaço que
possui propriedades de homogeneidade e isotropia é o esférico bidimensional
(Figura 1.3), em que a métrica é dada por ds2 = a2 dθ2 + a2 sin2 θdφ2 , onde a
é o raio da esfera, θ a colatitude e φ o azimute.
Cálculo Tensorial e Relatividade Geral 18

θ
φ

Figura 1.3: Representação de um espaço homogêneo e isotrópico.

Fonte: Autor (2016).

1.4 Cálculo Tensorial

Ao se generalizar vetores e tensores para espaços curvos é necessário
abranger alguns conceitos do cálculo em espaços planos. Primeiramente,
a derivada de um tensor não é necessariamente um tensor, visto que a trans-
formação de coordenadas toma a forma

∂T µν ∂ x̄α
 2 µ
∂ x ∂xν λω ∂xµ ∂ 2 xν ∂xµ ∂xν ∂ T̄ λω

λω
= T̄ + λ α ω T̄ + λ ω .
∂xρ ∂xρ ∂ x̄α ∂ x̄λ ∂ x̄ω ∂ x̄ ∂ x̄ ∂ x̄ ∂ x̄ ∂ x̄ ∂ x̄α
(1.4.1)
Para que a diferenciação de um tensor também possua carácter tensorial são
construídos os símbolos de Christoffel de primeira e segunda ordem, respec-
tivamente:  
1 ∂gρµ ∂gρν ∂gµν
Γρµν = ν
+ µ
− (1.4.2)
2 ∂x ∂x ∂xρ
e
1
Γλ µν = g λρ Γρµν . (1.4.3)
2
λ
Em algumas literaturas, estes são representados por [µ ν, ρ] e µν
, respec-
tivamente. São simétricos no índices µ e ν. Os símbolos de Christoffel não
Cálculo Tensorial e Relatividade Geral 19

são tensores, sua transformação de coordenadas é dada pela relação:

ρ ∂ x̄ρ ∂xσ ∂xλ ω ∂ x̄ρ ∂ 2 xω

Γ̄ µν = Γ + . (1.4.4)
∂xω ∂ x̄µ ∂ x̄ν σλ ∂xω ∂ x̄µ ∂ x̄ν

Se generaliza o conceito de diferenciação definindo as derivadas covarian-

tes para tensores contravariantes e covariantes, respectivamente:

∂T µν
∇ρ T µν = ρ
+ Γµ λρ T λν + Γν λρ T µλ , (1.4.5)
∂x

∂Tµν
∇ρ Tµν = ρ
− Γλ µρ Tλν − Γλ νρ Tµλ . (1.4.6)
∂x
Uma propriedade útil é que a derivada do tensor métrico é sempre nula:

∂gµν λ λ ∂gµν
∇ρ gµν = − Γ µρ gλν − Γ νρ gµλ = − Γνρµ − Γµρν
∂xρ  ∂x
ρ
 
∂gµν 1 ∂gµν ∂gρν ∂gµρ 1 ∂gµν ∂gρµ ∂gνρ
= − + − − + −
∂xρ 2 ∂xρ ∂xµ ∂xν 2 ∂xρ ∂xν ∂xµ
(1.4.7)

= 0.

Agora, seja xµ = xµ (u) uma curva parametrizada por u, define-se a deri-

vada absoluta para tensores contravariantes como

δT µν ∂xρ
= ∇ρ T µν (1.4.8)
δu ∂u

e, para covariantes
δTµν ∂xρ
= ∇ρ Tµν . (1.4.9)
δu ∂u
Como os símbolos de Christoffel são nulos em espaços planos, as derivadas
se reduzem ao caso do cálculo diferencial usual.
Cálculo Tensorial e Relatividade Geral 20

1.5 Geodésicas
As geodésicas costumam ser definidas, em espaços planos, como o menor
caminho entre dois pontos. Este conceito pode ser generalizado para espaços
curvos por uma definição mais formal [11]. Na Figura 1.4 representa-se um
conjunto de curvas geodésicas em um espaço curvo, a elipsóide triaxial.

Figura 1.4: Geodésicas em uma elipsóide triaxial.

Fonte: Site: Wikipedia 3 .

Definição: Uma geodésica é uma curva cujo valor é estacionário com

respeito a pequenas variações arbitrárias, mantendo os extremos fixos.
Z B
δ ds = 0. (1.5.1)
A

Supondo que a curva é parametrizada pelo próprio comprimento de arco

xµ = xµ (s) e definindo, por simplificidade, pµ = dxµ /ds como um vetor
unitário tangente à curva

dxµ dxν
gµν = pν pν = 1 (1.5.2)
ds ds
3
Disponível em: <[Link] Acesso em: agosto de
2016.
Cálculo Tensorial e Relatividade Geral 21

e4 L = gµν pµ pν . Portanto, a equação da geodésica é dada por 5

d ∂L ∂L
λ
− λ = 0, (1.5.3)
ds ∂p ∂x

e que pode ser escrita de modo explícito como

d 2 xλ λ dxµ dxν
+ Γ µν = 0. (1.5.4)
ds2 ds ds

As soluções destas equaçõe são determinadas de forma única com as condições

iniciais de xµ e dxµ /ds. Para um parâmetro t qualquer, a equação tem a
seguinte forma

d 2 xλ λ dxµ dxν d2 t/ds2 dxλ

+ Γ µν = − . (1.5.5)
dt2 dt dt (dt/ds)2 dt

Comumente é utilizado um parâmetro tal que dt/ds = 0 e a equação

(1.5.4) pode ser utilizada. Note que a equação (1.4.8) é idêntica a equação
da geodésica, portanto
δ dxλ
= 0. (1.5.6)
δs ds
diz-se que o vetor pλ propaga-se paralelamente ao longo da curva. O conceito
de transporte paralelo pode ser aplicado para definir a curvatura intrínseca
de um espaço. Se um vetor se propaga paralelamente por uma curva fechada
em um espaço plano, o vetor final é paralelo ao inicial, já em um espaço
curvo, isto não é necessariamente válido.
Um tipo de geodésica que deve ser mencionado são as geodésicas nulas,
que ocorrem quando, dado um parâmetro λ, o seguinte conjunto de equações
4
Na física, esta é a lagrangiana.
5
Apêndice A
Cálculo Tensorial e Relatividade Geral 22

deve ser satisfeito:

d2 x ρ ρ dxµ dxν
+ Γ µν = 0, (1.5.7)
dλ2 dλ dλ
gµν pµ pν = 0. (1.5.8)

Estas são as geodésicas que descrevem o movimento da luz na TRG.

1.6 Curvatura
O tensor de Riemann (ou Riemann-Christoffel) descreve a curvatura in-
trínseca do espaço, e é utilizado para construir outros tensores que auxiliam
na descrição da curvatura do espaço-tempo na TRG. Para construir este ten-
sor, assume-se a existência de duas curvas, xλ e xω tal que a junção das duas
seja uma curva fechada num espaço curvo. Se um vetor Uµ é transportado
paralelamente em relação a estas duas curvas, como já mencionado, o vetor
final não é, necessariamente, paralelo ao inicial. Esta afirmação é análoga a
de que as derivadas não são comutativas em espaços curvos

∇λ ∇ω Uµ 6= ∇ω ∇λ Uµ . (1.6.1)

Define-se o tensor de curvatura a partir desta diferença:

∇µ ∇ν Uλ − ∇ν ∇µ Uλ = Rσ λµν Uσ (1.6.2)

com
∂Γσ λν ∂Γσ λν
Rσ λµν = − + Γω λν Γσ ωµ − Γω λµ Γσ ων , (1.6.3)
∂xµ ∂xν
denominado tensor de Riemann ou de curvatura. Ele têm este nome devido
ao teorema que segue [11].
Cálculo Tensorial e Relatividade Geral 23

Teorema 1 É condição suficiente e necessária para que um espaço seja

plano que
Rσ λµν = 0. (1.6.4)

O teorema também é válido pra o tensor de Riemann covariante Rωλµν =

gωσ Rσ λµν . O tensor de Riemann é anti-simétrico nos dois último índices (µ
e ν) e tem simetria cíclica nos subíndices

Rσ λµν + Rσ µνλ + Rσ νλµ = 0. (1.6.5)

A partir deste pode-se construir outros tensores de interesse, como o tensor

de Ricci e o escalar de curvatura. O tensor de Ricci é definido pela contração

Rλν = Rσ λσν . (1.6.6)

O tensor de Ricci é simérico. O escalar de curvatura é dado por

R = g µν Rµν . (1.6.7)

Definidos estes tensores, o tensor de Einstein é definido a partir da iden-

tidade de Bianchi, dada por

∇ρ Rσ λµν + ∇µ Rσ λρν + ∇ν Rσ λρµ = 0 (1.6.8)

tal que, com algumas manipulações obtém-se a identidade

 
R
∇µ Rµν − gµν = 0, (1.6.9)
2

da qual define-se Gµν = Rµν − 12 Rgµν como o tensor de Einstein, que é

simétrico e tem derivada covariante nula.
Cálculo Tensorial e Relatividade Geral 24

1.7 Tensores Relativos

Tensores relativos são uma subclasse dos tensores [14] e se transformam
pela relação
 w
 ∂x  ∂ x̄µ1 ∂xβ1
τ̄ µ1 µ2 ···µm
ν1 ν2 ···νn =   · · · · · · τ α1 α2 ···αmβ1 β2 ···βn (1.7.1)
∂ x̄ ∂xα1 ∂ x̄ν1

na qual |∂x/∂ x̄| é o Jacobiano da transformação de xµ para x̄ν . O expoente

w denota o peso do tensor relativo. Diz-se que τ é um tensor relativo de peso
w.
Esta classe de tensores obedece as mesmas propriedades já definidas:

• A soma de dois tensores de mesma ordem (m, n), mesmo peso w e,

definidos no mesmo ponto, resultam em um tensor de mesma ordem e
mesmo peso que os anteriores.

• O produto de dois tensores relativos definidos no mesmo ponto, um de

ordem (m, n) e peso w e outro de ordem (r, s) e peso u resulta em um
tensor de ordem (m + r, n + s) e peso w + u.

Dois exemplos importantes de tensores relativos são o determinante do

tensor métrico e o elemento infinitesimal de volume. A partir destes se cons-
trói o invariante de volume. Por definição:

∂xµ ∂ x̄ν
ḡαβ = gµν (1.7.2)
∂xα ∂ x̄β

portanto,  2
 ∂x 
ḡ =   g (1.7.3)
∂ x̄
com ḡ e g sendo os determinantes do tensor métrico em relação as respectivas
coordenadas. Já o invariante de volume de um espaço quadridimensional se
Cálculo Tensorial e Relatividade Geral 25

transforma pela relação:  

 ∂ x̄ 
d x̄ =   d4 x
4
(1.7.4)
∂x
Com estes é construído o invariante de volume quadri-dimensional

p p
|ḡ|d4 x̄ = |g|d4 x. (1.7.5)

A partir destes tópicos do cálculo tensorial é possível iniciar uma discussão

da TRG.
Capítulo 2

Teoria da Relatividade Geral

Este capítulo tratará do desenvolvimento da TRG. É enunciado o princí-

pio da equivalência e obtidas as equações de Einstein a partir do princípio
variacional. Em seguida é verificado que estas retornam a equação de Pois-
son para o potencial na mecânica Newtoniana usando uma aproximação de
campos fracos.

2.1 Breve Histórico

Antes de abordar os tópicos de um assunto qualquer é importante discor-
rer sobre sua necessidade, visando uma maior compreensão de sua relevân-
cia e dos motivos históricos que levaram à sua gênese. Consequentemente
propiciando uma visão conceitual mais ampla. Aqui será dada uma breve
introdução histórica.
Segundo a teoria da gravitação de Newton, interações entre corpos ocor-
rem instantaneamente e, mudanças de coordenadas devem obedecer às trans-
formadas de Galileu, fazendo com que as forças sejam independentes do re-
ferencial adotado [15].

26
Cálculo Tensorial e Relatividade Geral 27

Com o advento do eletromagnetismo no século XIX surgiram alguns pro-

blemas teóricos. As equações de Maxwell não eram invariantes sob trans-
formações galileanas e aparentava ser constante a velocidade de propagação
das ondas eletromagnéticas (c) no vácuo [8]. Tais observações trouxeram um
rompimento com a mecânica Newtoniana e resultaram no ressurgimento da
hipótese de um referencial absoluto, o éter luminífero, um meio sutil no qual
as equações de Maxwell seriam válidas. Visando o primeiro problema, no ano
de 1903 o físico Hendrick Lorentz desenvolveu as equações hoje conhecidas
como transformadas de Lorentz, as quais mantinham as equações de Maxwell
invariantes, impondo que a velocidade da luz não pudesse ser alcançada ou
ultrapassada. As transformações de Galileu surgiam como caso especial para
velocidades baixas [16].
Em 1905, Albert Einstein percebeu que os eventos eletromagnéticos que
ocorrem em um sistema de coordenadas devem ocorrer da mesma forma em
qualquer outro sistema se movendo com velocidade constante e postulou o
princípio de Covariância, ditando que as leis da física devem ter a mesma
forma para quaisquer referenciais inerciais. A partir deste, viu-se a necessi-
dade de uma teoria mecânica que obedecesse as transformadas de Lorentz,
e surge a Teoria da Relatividade Restrita (TRR) como generalização da me-
cânica Newtoniana [17].
Em 1908, o físico Hermann Minkowski (Figura 2.1) forneceu uma inter-
pretação geométrica para a TRR, argumentando que espaço e tempo não são
entidades distintas, estando conectadas pelo invariante

ds2 = cdτ 2 = cdt2 − dx2 − dy 2 − dz 2 , (2.1.1)

Cálculo Tensorial e Relatividade Geral 28

Figura 2.1: Físico Hermann Minkowski.

Fonte: Site: Space Time Society. 1

do qual têm-se a métrica de Minkowski, como é conhecida hoje.

 
1 0 0 0
 
 0 −1 0
 
0 
ηµν =  (2.1.2)
 0 0 −1 0
 

 
0 0 0 −1

Minkowski também desenvolveu o conceito de linha de Universo (world line)

de uma partícula, que é a trajetória desta no espaço quadri-dimensional e
chamou de ponto do Universo (world-point) uma coordenada (ct, x, y, z) [3].
Um aspecto interessante da intepretação geométrica é a análise decorrente
dos diagramas de Minkowski (Figura 2.2): na região do cone de luz passado
e futuro (em azul), têm-se que ds2 > 0 e este intervalo é dito do tipo-tempo e
dois eventos ocorrendo nele podem ter relação causal; na região fora do cone,
ds2 < 0, o intervalo é do tipo-espaço e não há relação causal entre eventos;
as retas que separam estas regiões ocorrem para ds2 = 0 e definem intervalos
tipo-luz [7].
A partir das ideias de geometrização do espaço-tempo e do princípio de
1
Disponível em: <Disponível em: <[Link]
Acesso em: setembro de 2016.
Cálculo Tensorial e Relatividade Geral 29

Covariância, Einstein postulou o princípio de Equivalência e iniciou uma

busca por equações que generalizassem a TRR, de modo a incluir campos
gravitacionais. De início, as tentativas se resumiam a englobar os campos
dentro do espaço-tempo de Minkowski porém, após vários esforços em vão,
surgiu a ideia de que a gravitação fosse causada por uma curvatura no espaço-
tempo, causada pela presença de energia ou matéria em determinada região
do Universo. Tal curvatura seria totalmente descrita pelo tensor métrico,
sendo este o análogo relativístico do potencial.

Figura 2.2: Diagrama de Minkowski.

Fonte: Autor (2016).

A base conceitual da TRG estava completa em 1913, sendo necessárias so-

mente as equações de campo. Após diversas tentativas, em novembro de 1915
foram publicadas em sua forma correta [1, 18]. A forma final das equações
foram publicadas tanto por Einstein quanto por David Hilbert, o primeiro
com um atraso de cinco dias em relação ao último, o que causou uma con-
trovérsia sobre a autoria das mesmas. Uma discussão aprofundada se dá na
referência [19].
Cálculo Tensorial e Relatividade Geral 30

Logo após a obtenção das equações de campo, novas predições surgiram

da teoria, como o avanço do periélio dos planetas e a deflexão gravitacional
da luz. A primeira já havia sido observada, porém um planeta em especial
-Mercúrio- sofria um avanço muito maior que o previsto pela teoria Newto-
niana, sendo explicado somente com a TRG [21]. O segundo foi observado
experimentalmente em 1919, mas testes com maior precisão ocorreram ape-
nas na década de 60, corroborando fortemente para a validação da teoria
[37]. Outras predições como a existência de buracos negros e de ondas gravi-
tacionais foram verificadas recentemente e hoje a teoria é amplamente aceita,
constituindo a base da cosmologia moderna [20].

2.2 O Princípio de Equivalência

O Princípio de equivalência é uma generalização do princípio da mecânica
Newtoniana de que não há distinção entre massas gravitacionais e inerciais,
também conhecido como formulação fraca do princípio de equivalência, o
qual é suportado pelos experimentos de Eötvös e Dicke [23, 24]. Para uma
formulação da forma forte toma-se o seguinte experimento de pensamento
(gedankenexperiment): Sejam dois observadores (1 e 2), ambos dentro de
elevadores sem visão do exterior. O primeiro está em repouso e sob ação de
um campo gravitacional constante enquanto o segundo não está sob efeitos
gravitacionais, mas é acelerado por uma força externa na direção oposta à
do campo em 1, como mostra a Figura 2.3.
Cálculo Tensorial e Relatividade Geral 31

Figura 2.3: Observadores 1 e 2 descrevem o movimento de uma bola em queda

livre.
Fonte: Autor (2016).

Supondo que a aceleração aplicada no observador 2 seja igual àquela

propiciada pelo campo gravitacional em 1, e que ambos observadores aban-
donem bolas idênticas, a descrição do movimento que cada observador dará
à respectiva bola será a mesma, não sendo possível distinguir o movimento
devido a forças inerciais (observador 2) e aquele devido a forças gravitaci-
onais (observador 1). O princípio da equivalência em sua forma forte pode
ser escrito como: “Forças inerciais e gravitacionais são equivalentes do ponto
de vista físico, sendo impossível distinguí-las experimentalmente”. Algumas
consequências de grande importância são inferidas diretamente da formula-
ção forte.

• Não existência de referenciais inerciais.

Em um espaço de Minkowski, a linha de Universo de uma partícula livre é

uma linha reta, ou uma geodésica [24], representada por

d 2 xµ
= 0, (2.2.1)
dτ 2
Cálculo Tensorial e Relatividade Geral 32

na qual τ é o tempo próprio. Como a geodésica é definida independentemente

do sistema de coordenadas, a equação de movimento em um referencial qual-
quer é dada por (1.5.4)

d 2 xλ λ dxµ dxν
+ Γ µν = 0, (2.2.2)
dτ 2 dτ dτ

na qual os símbolos de Christoffel são interpretados como a soma das forças

inerciais e gravitacionais agindo sobre a partícula. Portanto, na presença
de um campo gravitacional é impossível tomar Γλ µν = 0 em todos pontos
e o espaço não pode ser descrito pela métrica de Minkowski. Surge aqui a
necessidade das equações de campo da TRG, que devem descrever o tensor
métrico. Determinado gµν é possível descrever a trajetória de um corpo.
Posta a impossibilidade de diferenciar forças gravitacionais e inerciais,
conclui-se que não é possível distinguir referenciais inerciais e não-inerciais.
Porém, ao se tomar uma região suficientemente pequena do espaço, a métrica
pode ser aproximada localmente pela de Minkowski, diz-se que esta região
define um referencial localmente inercial [15] no qual é possível utilizar quais-
quer equações da TRR.

• Desvio gravitacional da luz.

Tomemos dois observadores no elevador, como no experimento mental an-

terior. Suponhamos que, no lugar de bolas em queda livre exista um pequeno
buraco na parede, pelo qual entra um feixe de luz. Um observador em um
elevador livre de campos e de aceleração verá a luz se movimentando em linha
reta, já um observador em um referencial acelerado perceberá que a luz exe-
cuta um movimento curvo (Figura 2.4). Como pelo Princípio de equivalência
não é possível fazer distinção entre um observador acelerado e um observador
inercial em um campo gravitacional, concluí-se que ele só é obedecido se a
Cálculo Tensorial e Relatividade Geral 33

luz sofrer um desvio quando sob ação de um campo gravitacional.

Figura 2.4: Feixe de luz visto de um referencial acelerado.

Fonte: Autor(2016).

Este desvio produz um efeito conhecido como lente gravitacional (Figura

2.5), o qual é causado pela deflexão da luz de uma estrela por um outro
corpo massivo tal que a luminosidade da estrela é ampliada, sendo possível
determinar sua massa[25].

Figura 2.5: Imagem ilustrativa do efeito de lente gravitacional.

Fonte: ALMA (ESO/NAOJ/NRAO)/C. Collao/Japan Meteorological Agency.
Cálculo Tensorial e Relatividade Geral 34

2.3 Equações de Campo

As equações de campo da TRG podem ser derivadas pelo princípio variaci-
onal, como feito por Chandrasekhar (1972) [26] ou, de forma mais elementar
como feito na referência [24]. Aqui será adotada a formulação variacional
que, além de maior elegância, provê também um melhor entedimento físico
da fonte de campo gravitacional e uma maneira mais fácil de englobar campos
eletromagnéticos na teoria [27].
Seja L a lagrangiana escalar do sistema e Ω um volume no espaço-tempo
quadridimensional, a ação é definida como

√
Z
S= L −gd4 x (2.3.1)
Ω

√
sendo g o determinante da métrica e o termo −gd4 x o elemento invariante
√
de volume (o produto L −g é denominado densidade lagrangiana). O sinal
negativo na raíz se dá devido ao determinante da métrica ser sempre negativo
para um espaço-tempo real [28]. Pelo princípio de mínima ação [7]

δS = 0, (2.3.2)

Separando a lagrangiana em dois termos, um para o campo gravitacional e

outro para campos externos (agindo como fontes de campo gravitacional):
L = LG − 2κLM , com κ sendo constante gravitacional relativística e a cons-
tante 2 adotada por conveniência. A equação (2.3.2) pode ser reescrita como

√
Z
δ −g (LG − 2κLM ) d4 x = 0. (2.3.3)
Ω

Para o campo gravitacional, a ação é conhecida como ação de Einstein-

Hilbert, e a lagrangiana é dada pelo invariante de Ricci (LG = R) [29]. Foi
Cálculo Tensorial e Relatividade Geral 35

David Hilbert quem postulou este escalar como a lagrangiana. O argumento

para a escolha é: as equações de campo devem conter o tensor métrico e suas
derivadas de até segunda ordem e as equações surgem de uma variação na
ação, portanto, a lagrangiana deve conter derivadas de primeira ordem em
gµν [28, 30]. Porém, é possível manipular as equações para diminuir a ordem
das derivadas, sendo o escalar de curvatura uma lagrangina aparentemente
aceitável. Retomando a ação

√
Z
δ −g (R − 2κLM ) d4 x = 0 (2.3.4)
Ω

expandindo R, pode-se escrever a equação como

√ √
Z Z
µν
δ Rµν g −gd4 x − δ 2κLM −gd4 x = 0. (2.3.5)
Ω Ω

Primeiramente, analisar-se-á o termo referente ao campo gravitacional:

√ √ √

Z Z Z
µν 4 µν 4
δ Rµν g −gd x = δ (Rµν ) g −gd x + Rµν δ g µν −g d4 x,
Ω Ω Ω
(2.3.6)
usando a definição (1.6.3) para o tensor de Ricci

∂Γρ µν ∂Γρ µρ
 
σ ρ σ ρ
δRµν = δ − + Γ µν Γ σρ − Γ µρ Γ νσ , (2.3.7)
∂xρ ∂xν

é possível simplificar a expressão usando o seguinte lema da geometria Rie-

manniana: “É sempre possível escolher um sistema de coordenadas no qual
todos os símbolos de Christoffel são nulos em um determinado ponto”, tal
sistema de coordenadas é denominado sistema de coordenadas geodésicas
[14]. Adotando tal sistema, os símbolos são nulos no ponto e, portanto as
Cálculo Tensorial e Relatividade Geral 36

derivadas covariantes condizem com as parciais

δRµν = ∇ρ δΓρ µν − ∇ν δΓρ µρ .

(2.3.8)

Note que esta é uma equação tensorial e, portanto deve ser válida em qual-
quer sistema de coordenadas. Como as derivadas covariantes da métrica são
sempre nulas

√ √ 
−gg µν δRµν = −g ∇ρ g µν Γρ µν − ∇ν g µν δΓρ µρ ,



(2.3.9)

que pode ser reescrita fazendo uma troca de índices

√ √
−gg µν δRµν = −g ∇α g µν δΓα µν − g µα δΓρ µρ .

(2.3.10)

Nesta, o termo a ser derivado é um tensor contravariante de ordem 1. Sa-

2
bendo que a expressão do Laplaciano pode ser escrita como

1 ∂ √
∇µ V µ = √ −gV µ


µ
(2.3.11)
−g ∂x

a integral de ação toma a forma

√ ∂ √
Z Z
−gg µν δRµν d4 x = −gg µν δΓα µν − g µα δΓρ µρ d4 x


α
(2.3.12)
Ω Ω ∂x

usando o teorema de Gauss, e lembrando que os extremos devem ser mantidos

fixos, a integral deve ser nula, ou seja,

√
Z
−gg µν δRµν d4 x = 0. (2.3.13)
Ω
2
Apêndice B
Cálculo Tensorial e Relatividade Geral 37

Agora é necessário analisar o termo com a métrica e o determinante

√ √ √
Z Z Z
Rµν δ −gg µν d4 x = µν 4
Rµν g µν δ −gd4 x (2.3.14)


Rµν −gδg d x+
Ω Ω Ω

o segundo termo pode ser escrito como [14]

√ 1 1√
δ −g = − √ δg = − −ggµν δg µν . (2.3.15)
2 −g 2

Quando substituído na expressão (2.3.5), retorna

√ √
Z Z  
µν 4 1
−g Rµν − Rgµν δg µν d4 x.


Rµν δ −gg d x= (2.3.16)
Ω Ω 2

Retomando a ação para o campo gravitacional

√ √
Z Z  
1
δ LG −gd4 x = −g Rµν − Rgµν δg µν d4 x (2.3.17)
Ω Ω 2

no caso dos campos externos serem nulos (LM = 0), como δg µν é livre para
√
variar e −g 6= 0, esta equação traz as equações de campo no vácuo em
forma covariante:
1
Rµν − Rgµν = 0. (2.3.18)
2
Caso haja campos externos é necessário desenvolver o termo com LM .
Supondo que a lagrangiana dependa da métrica e de sua derivada primeira,
usando por simplicidade, a notação g µν ,α = ∂g µν /∂xα

∂ √ √
Z Z  
4

µν ∂
−gLM δg ,α d4 x

µν
δ LM d x = −gLM δg + µν
Ω Ω ∂g µν ∂g ,α
(2.3.19)
Cálculo Tensorial e Relatividade Geral 38

e notando que

√ √
 
∂  µν ∂ ∂
µν
−gLM δg ,α = α −gLM δg − (2.3.20)
∂g µν ,α ∂x ∂g µν ,α
√
 
∂ ∂
−gLM δg µν


∂xα ∂g µν ,α

a integral anterior pode ser reescrita como

∂ √ √
Z Z   
4 ∂ ∂
δg µν d4 x



δ LM d x = µν
−gLM − α µν −gLM
Ω Ω ∂g ∂x ∂g ,α
(2.3.21)
pois o primeiro termo da equação (2.3.20) se anula ao usar o teorema de
Gauss na integral. Definindo o tensor energia-momento como [14]

∂ √ √
  
2
∂ ∂

Tµν =√ −gLM − α −gLM , (2.3.22)
−g ∂g µν ∂x ∂g µν ,α

tal que
√ 1√
Z Z
4
−gLM d x = −gTµν δg µν d4 x. (2.3.23)
Ω Ω 2
Portanto, a equação inicial para a ação é escrita como

√
Z  
1
δ −g Rµν − Rgµν − κTµν δg µν d4 x = 0 (2.3.24)
Ω 2

a partir da qual as equações de campo da TRG são obtidas em forma cova-

riante
1
Rµν − Rgµν = κTµν . (2.3.25)
2
Uma forma alternativa para as equações acima podem ser obtidas multiplicando-
as pelo tensor recíproco g µν
R = −κT (2.3.26)
Cálculo Tensorial e Relatividade Geral 39

que, quando substituindo nas equações iniciais retornam

 
1
Rµν = κ Tµν − T gµν . (2.3.27)
2

A priori não há indícios de que estas equações sejam fisicamente aceitáveis.

Tomar-se-á o limite de v << c para verificar se estas de fato retornam a
equação de Poisson [24]
∇2 Φ = 4πGρ, (2.3.28)

sendo G a constante Gravitacional, ρ a densidade do corpo e Φ o potencial

gravitacional. Também se busca impor uma limitação na constante κ.
É interessante mencionar que uma das soluções das equações de Eins-
tein, a métrica de Friedmann–Lemaître–Robertson–Walker (FLRW) prevê
um universo em expansão ou em contração [31]. Mas Albert Einstein, como
a grande maioria das pessoas da época, era um defensor do Universo estático
e eterno. Ao surgir o confronto de previsão e crença, a última foi mais forte
e, para que a solução retornasse um universo estático, adicionou outro termo
[32]. Com este, as equações tomam a forma

1
Rµν − Rgµν + Λgµν = κTµν (2.3.29)
2

com Λ da ordem de 10−12 e denominada constante cosmológica. Mas a so-

lução desta equação não era estável e, em 1929 foi observado por Edwin
Hubble que o universo está, de fato, em expansão, e a constante foi retirada.
De acordo com Einstein, este foi o maior erro de sua vida. Apesar disto, após
algumas décadas a constante foi retomada no estudo de teorias inflacionárias
e na teoria quântica de campos [33].
Cálculo Tensorial e Relatividade Geral 40

2.4 Limite Newtoniano

O limite Newtoniano assume uma aproximação de campos fracos, com
o campo tendo como fonte um fluido perfeito [34]. Também assume-se que
a métrica possa ser escrita como a métrica de Minkowski somada a uma
pequena perturbação hµν
gµν = ηµν + hµν (2.4.1)

e que, como v << c, o elemento de linha possa ser aproximado por

ds2 ≈ g00 c2 dt2 . (2.4.2)

Primeiramente, verifiquemos a equação de movimento da partícula em rela-

ção ao parâmetro x0 = ct

d2 x0 /ds2 µ
ẍµ + Γµ ωλ ẋλ ẋω = − x˙ (2.4.3)
(dx0 /ds)2

com ẋν = dxν /dx0 . Para µ = 0, obtém-se a identidade

d2 x0 /ds2
= −Γ0 ωλ ẋλ ẋω . (2.4.4)
(dx0 /ds)2

Esta, se substituída na equação inicial, pode ser omitida ao se notar que

ν
Γν ωλ  Γ0 ωλ 1c dxdt . Convencionando que indíces latinos somam somente
sobre componentes espaciais, isto é, i = 1, 2, 3:

ẍi + Γi ωλ ẋ
λ ω
ẋ ≈ 0 (2.4.5)

e, como ẋk << ẋ0

ẍi ≈ −Γi 00 . (2.4.6)
Cálculo Tensorial e Relatividade Geral 41

No limite quasi-estático, os símbolos de Christoffel são

1 ∂g00
Γi 00 = − g iλ λ (2.4.7)
2 ∂x

notando que g λi ≈ −δ λi , então

1 ∂g00
Γi 00 ≈ . (2.4.8)
2 ∂xi

Agora, usando as equações de campo na forma (2.3.27), e sabendo que o ten-

sor energia-momento para um fluido perfeito é dado por Tµν = (ρ + p/c2 ) uµ uν +
ηµν p, no qual ρ é a densidade, p a pressão e uµ a velocidade [34], no limite
Newtoniano Tµν ≈ T00 = ρc2 :

κ 2
R00 = ρc . (2.4.9)
2

Usando a definição do tensor de Ricci, e usando o resultado obtido para os

símbolos de Christoffel:

∂Γλ 00 1 ∂ 2 g00
R00 ≈ ≈ (2.4.10)
∂xλ 2 ∂xk ∂xk

portanto
∂ 2 g00
= κρc2 . (2.4.11)
∂xk ∂xk
A derivada segunda no índice k nada mais é que o laplaciano. Lembrando
que na mecânica Newtoniana o potencial é descrito pela equação de Poisson
(2.3.28) e que g00 = 1 + h00 , assumindo que

2
h00 = Φ (2.4.12)
c2
Cálculo Tensorial e Relatividade Geral 42

então
∂ 2Φ κ
k k
= c4 ρ (2.4.13)
∂x ∂x 2
e o limite Newtoniano é obtido com a imposição que

8πG
κ= (2.4.14)
c4

e as equações de campo são escritas como

1 8πG
Rµν − Rgµν = 4 Tµν . (2.4.15)
2 c

Portanto, verificou-se que a lagrangiana escolhida era fisicamente aceitável

e que, deduzidas as equações de campo, estas retornavam o limite Newtoni-
ano. Esta verificação é essencial para a aceitação da teoria, dado o sucesso
experimental da mecânica Newtoniana nos séculos que precederam a TRG.
Debruçemo-nos agora sobre as previsões feitas pela solução de Schwarzschild.
Capítulo 3

Solução de Schwarzschild

Este capítulo tratará da solução das equações de campo para um corpo esfe-
ricamente simétrico no vácuo, denominada solução de Schwarzschild. Após
determinada a métrica serão tratadas quantitativamente algumas previsões
da solução.

3.1 A Métrica
As equações de campo são não-lineares, possuindo soluções exatas de
difícil obtenção. Algumas das soluções conhecidas são as já mencionadas de
Schwarzschild e FLRW. Sendo a primeira, o caso mais simples a ser pensado,
representando o campo gravitacional de um ponto material no vácuo. A
Figura 3.1 visa dar um entendimento geométrico desta curvatura em uma
analogia com um espaço bidimensional. A equação a ser resolvida é

Rµν = 0. (3.1.1)

43
Cálculo Tensorial e Relatividade Geral 44

A solução da equação acima foi obtida pela primeira vez por Karl Schwarzs-
child1 em 1916 e é uma boa aproximação para tratar do campo gravitacional
de estrelas [35].

Figura 3.1: Um análogo bidimensional do espaço-tempo de Schwarzschild.

Fonte: Site: What When How. 2

Tomando (t, r, θ, φ) como as coordenadas que descrevem o espaço-tempo

(denominadas coordenadas de Schwarzschild), sendo θ e φ coordenadas po-
lares, t uma coordenada tipo-tempo e r uma coordenada radial. Assume-se
que

• O campo é estático: ∂gµν /∂t = 0;

• A solução é espacialmente simétrica com respeito a origem do sistema

de coordenadas;

• A solução é assintoticamente plana para r → ∞.

Uma métrica que satisfaz os requisitos impostos [34] é do tipo

ds2 = F (r)dt2 − G(r)dr2 − r2 dΩ2 , (3.1.2)

1
Karl Schwarzschild foi um físico alemão que além de TRG, participou de estudos em
Mecânica Quântica e Espectroscopia. Viveu entre 1873 e 1916.
2
Disponível em: <[Link]
space-time-into-a-black-hole-string-theory/> Acesso em: outubro de 2016.
Cálculo Tensorial e Relatividade Geral 45

na qual dΩ2 = dθ2 + sin2 θdφ2 . Se o campo é assintoticamente plano, deve-se

ter como condições de contorno

 F (r → ∞) = c2 ,
(3.1.3)
 G(r → ∞) = 1.

As equações de campo deverão determinar as funções F e G. Primeiramente

é necessário analisar os símbolos de Christoffel. Os componentes não-nulos
são:

F0 F0
Γ0 01 = , Γ1 00 = ,
2F 2G
G0 r
Γ1 11 = , Γ1 22 =− ,
2G G
r sin2 θ 1
Γ1 33 =− , Γ2 12 = ,
G r
1
Γ2 33 = − sin θ cos θ, Γ3 13 = ,
r
Γ3 23 = cot θ.

Partindo da equação (1.6.3), o tensor de Ricci possui os termos não-

diagonais nulos. Já os diagonais são

F 00 F0 F 0 G0 F0
 
R00 =− + + − = 0, (3.1.4)
2G 4G F G rG

F 00 F 0 G0
F0 G0
 
R11 = − + − = 0, (3.1.5)
2F 4F
F G rG
 0 
1 r F G
R22 = −1+ − = 0, (3.1.6)
G 2G F G
R33 = R22 sin2 θ = 0, (3.1.7)

resultando em três equações independentes. Tomando as duas primeiras e,

Cálculo Tensorial e Relatividade Geral 46

lembrando das condições de contorno em (3.1.3):

F 0 G0 d
+ = 0 =⇒ (F G) = 0 =⇒ F G = c2 (3.1.8)
F G dr

que, se substituída em (3.1.7), retorna

dF dr
= . (3.1.9)
c2 − F r

Cuja solução é 

2 k
F (r) = c 1 + , (3.1.10)
r
 −1
k
G(r) = 1 + . (3.1.11)
r
Substituidas as equações (3.1.10) e (3.1.11) na métrica de Schwarzschild
(3.1.2), levam ao resultado
 
2 2 k 1
ds = c 1+ dt2 − dr2 − r2 dΩ2 (3.1.12)
r 1 + kr

sendo necessário determinar a constante k. Tomando r → ∞, como a solução

é assintoticamente plana, a aproximação de campos fracos dita que g00 →
2
1+ c2
Φ com Φ = −GM/r, portanto k = − 2GM
c2
e

   −1
2 2GM 2 2 2GM
ds = 1 − 2 c dt − 1 − 2 dr2 − r2 dΩ2 . (3.1.13)
cr cr

2GM
Nota-se que se r = c2
a solução apresenta uma singularidade, com
g00 → 0 e g11 → ∞. Tal coordenada é determinada raio de Schwarzschild
(rS ) e, para r ≤ rS as coordenadas (t, r, θ, φ) deixam de ser válidas, visto que
a métrica teria um sinal trocado - a coordenada tipo tempo se tornaria tipo
espaço e a coordenada radial, tipo tempo. Porém, o determinante do ten-
Cálculo Tensorial e Relatividade Geral 47

sor métrico em rS é g = −r4 sin2 θ, não apresentando singularidade alguma,

portanto, esta não deve ser uma propriedade intrínseca da solução (como em
r = 0) mas sim das coordenadas utilizadas [14]. Um tipo de transforma-
ção de coordenadas que faz com que a singularidade desapareça são as de
Eddington-Finkelstein, ou também as de Lemaître, como será visto algumas
seções adiante.

3.2 Geodésicas e Órbitas

Definida a métrica de Schwarzschild, é de interesse físico determinar o
movimento de uma partícula teste neste espaço. Seja a equação da geodésica
(1.5.3), considerando o tempo próprio τ como parâmetro e definindo que
ẋ = ∂x/∂τ , a lagrangiana L é dada por

 rS  2  rS −1 2  
L = c2 1 − ṫ − 1 − ṙ − r2 θ̇2 + sin2 θφ̇2 (3.2.1)
r r

da qual se obtém quatro equações geodésicas, das quais duas retornam cons-
tantes de movimento (E e L) relacionadas a energia e momento angular [2].
Assumindo que o movimento ocorra no plano θ = π/2, as equações a serem
resolvidas são dadas por:

 rs −1 rs  rS −2 2 rS
1− r̈ − 2 1 − ṙ + 2 c2 ṫ2 − rφ̇2 = 0 (3.2.2)
r 2r r 2r
 rS 
1− ṫ = E, (3.2.3)
r
r2 φ̇ = L. (3.2.4)

Como o parâmetro assumido é o tempo próprio, notando que ds2 = c2 dτ 2 ,

então L = c2 . Usando as constantes de movimento, a equação inicial é
Cálculo Tensorial e Relatividade Geral 48

reescrita como

2 rS −1 2 2 L2 
 rS −1 2
c = 1− E c − 2 − 1− ṙ (3.2.5)
r r r

que pode ser facilmente reescrita como

 2
1 dr
+ Φef f (r) = ε, (3.2.6)
2 dτ

c2
na qual ε = 2
(E 2 − 1) e

GM L2 GM L2
Φef f = − + 2− 2 3 . (3.2.7)
r 2r cr

Sendo Φef f o potencial efetivo, que difere daquele da mecânica Newtoni-

ana para uma força central apenas pelo termo com 1/r3 , que é uma corre-
ção relativística e implica em algumas conclusões diferentes das Newtonia-
nas. Primeiramente, o efeito assintótico é invertido, na teoria newtoniana
lim Φef f = −∞, já na relativística, lim Φef f = ∞, como pode ser visto na
r→0 r→0
Figura 3.2.
Também, na mecânica Newtoniana os pontos de equilíbrio existem para
qualquer configuração adotada de parâmetros, sendo sempre possível obter
órbitas estáveis. Enquanto que na TRG ao se calcular os pontos de equilíbrio
por 
dΦef f 
= 0, (3.2.8)
dr r=r0
obtém-se " 1/2 #
L2 12G2 M 2

r0 = 1± 1− (3.2.9)
2GM c2 L2

sendo visível que só há órbitas estáveis se L2 ≥ 12G2 M 2 /c2 e também, que

o limite clássico é obtido (Figura 3.2) para L2  12G2 M 2 /c2 .
Cálculo Tensorial e Relatividade Geral 49

0,04 L=8
L=6
0,02 L=4
Φef f

-0,02
Schwarzschild
Newton
-0,04

0 25 50 75 100
r

Figura 3.2: Comparação dos potenciais clássicos e relativísticos variando o

parâmetro L e assumindo GM = 1 e c = 1.
Fonte: Autor (2016).

A diferença nas órbitas clássicas e relativísticas podem ser vistas nas

Figuras 3.3, 3.4, 3.5. Tais imagens foram feitas em simulações usando o
3
método de Euler para diferentes valores de L.

3
Apêndice C
Cálculo Tensorial e Relatividade Geral 50

Clássica
L=4 Relativística
25

-25
-25 0 25 50

Figura 3.3: Comparação das órbitas clássicas e Relativísticas adotando GM = 1,

c = 1 e L = 4.
Fonte: Autor (2016).

Clássica
50 Relativística
L=6

25

-25

-50
-50 -25 0 25 50

Figura 3.4: Comparação das órbitas clássicas e Relativísticas adotando GM = 1,

c = 1 e L = 6.
Fonte: Autor (2016).
Cálculo Tensorial e Relatividade Geral 51

100
Clássica
75 L=8 Relativística

50

25

-25

-50

-75
-75 -50 -25 0 25 50 75 100

Figura 3.5: Comparação das órbitas clássicas e Relativísticas adotando GM = 1,

c = 1 e L = 8.
Fonte: Autor (2016).

Verifica-se que para L muito pequeno a órbita difere radicalmente da

predição clássica, mas conforme este cresce as órbitas se tornam próximas,
como esperado. Também é possível visualizar que não há órbitas ligadas e o
já mencionado avanço no periélio, que será discutido na próxima seção.

3.3 O Periélio de Mercúrio

O periélio é o ponto mais próximo do Sol em uma órbita, e todos planetas
sofrem um desvio neste, antes pensado que devido à influência gravitacional
dos outros planetas. A Figura 3.6 exemplifica como se dá o avanço do periélio.
O de Mercúrio é o mais conhecido, pois apresenta uma anomalia: seu desvio
era aproximadamente o dobro do previsto pelos físicos. Com o advento da
TRG verificou-se que o avanço no periélio ocorre mesmo sem perturbações
externas e, sendo Mercúrio o planeta mais próximo do sol, seu avanço deveria
ser maior. A verificação de que a previsão teórica estava de acordo com os
dados experimentais foi um dos primeiros testes a corroborarem fortemente
Cálculo Tensorial e Relatividade Geral 52

com a teoria. Aqui será determinado qual o avanço do periélio tomando como
base as referências [2, 34, 36].

Figura 3.6: Representação do Periélio de um planeta.

Fonte: Site: Wikipedia 4 .

Na mecânica Clássica, o movimento de um corpo em órbita é descrito por

uma elipse fechada, dada pela equação

1 GM
= 2 (1 + e cos φ) , (3.3.1)
r L

na qual e é a excentricidade, φ é o azimute, G a constante gravitacional, M

a massa do corpo que gera o campo gravitacional e L o momento angular do
corpo em órbita. Para determinar a análoga relativística usa-se as equações
(3.2.4) e (3.2.6), das quais é obtido

2
c2 r 2 c2 E 2
   
dr 2 2GM
+r 1+ 2 1− 2 − 2 r4 = 0. (3.3.2)
dφ L cr L
4
Disponível em: < [Link] > Acesso em: ou-
tubro de 2016.
Cálculo Tensorial e Relatividade Geral 53

A qual, com a mudança de variável u = L2 /(GM r) pode ser escrita como

2
2G2 M 2 3 2L2 ε

du 2
− u + u − 2u = (3.3.3)
dφ L 2 c2 G2 M 2

que, se derivada em φ, é:

d2 u 3G2 M 2 2
− u + u − 1 = 0. (3.3.4)
dφ2 L2 c2

Notando que em u = 1 tem-se o ponto de retorno clássico, a equação (3.3.4)

pode ser tratada como uma perturbação na equação Newtoniana. Se u0
corresponde à solução Newtoniana e up é uma pequena perturbação

u = u0 + up . (3.3.5)

Desprezando os termos quadráticos em up e os cruzados, pode-se separar

(3.3.4) em duas equações. Para o termo u0 tem-se

d2 u0
+ u0 − 1 = 0, (3.3.6)
dφ2

que tem como solução a equação da elipse (3.3.1). Já para o termo pertur-
bativo up :
d2 u p 3G2 M 2 2
+ u p = u, (3.3.7)
dφ2 L 2 c2 0
que, quando substituída a solução para u0 , possui como uma possível solução

3G2 M 2
  
1 2 1 2
up = 1 + e + eφ sin φ − e cos (2φ) . (3.3.8)
L2 c2 2 6

Notando que o primeiro termo é uma constante e o terceiro oscila em torno

da origem, estes podem ser desprezados no cálculo do avanço do periélio.
Cálculo Tensorial e Relatividade Geral 54

O segundo termo, que é cumulativo a cada órbita, deve ser considerado [2].
Têm-se como resultado

G2 M 2
u = 1 + e cos φ + 3e φ sin φ. (3.3.9)
L 2 c2

Portanto, a solução pode ser escrita de forma mais compacta como

1 GM
= 2 {1 + e cos [φ (1 − σ)]} (3.3.10)
r L

com σ = 3G2 M 2 /L2 c2 . A forma anterior é retomada expandindo o coseno

em torno de σ.
O período da órbita pode ser aproximado por uma aproximação binomial
como
3G2 M 2
 
T ≈ 2π 1 + (3.3.11)
L 2 c2
portanto, o avanço do periélio é dado pela expressão

6πG2 M 2
δφ = . (3.3.12)
L2 c2

Para substituir o momento angular L, compara-se (3.3.1) com a equação da

elipse
1 1 + e cos φ
= (3.3.13)
r a (1 − e2 )
na qual a é o semieixo maior. Podemos aproximar L2 ≈ GM a(1 − e2 ), tal
que
6πGM
δφ = . (3.3.14)
c2 a(1
− e2 )
Usando os valores tabulados pela agência Nasa 5 , para o planeta Mercúrio
5
Disponível em: <[Link]
Acesso em outubo de 2016.
Cálculo Tensorial e Relatividade Geral 55

a precessão do periélio é

δφ ≈ 0.10300 /órbita ≈ 4300 /século (3.3.15)

Apesar de ser o maior avanço de periélio do sistema solar, este valor ainda
é consideravelmente baixo, isto representa que a cada século há um avanço
de 43 segundos de arco. Apesar do valor baixo, a medida está de acordo com
o valor determinado experimentalmente [20].

3.4 Desvio Gravitacional da Luz

Como já discutido anteriormente, uma das implicações do princípio de
equivalência é o desvio gravitacional da luz. Nesta seção este desvio será
tratado de forma quantitativa dentro do espaço de Schwarzschild. Supondo
um feixe de luz que passa nos arredores de uma estrela, a deflexão sofrida
por este pode ser determinado fazendo uso das equações geodésicas. Feixes
de luz percorrem geodésicas nulas, que nas teorias relativísticas também são
chamadas geodésicas tipo-luz, implicando que a equação 3.2.1 não possui
mais valor c2 , mas valor nulo. Será usado um parâmetro λ, visto que o
tempo próprio não é válido. Usando as equações (3.2.3) e (3.2.4), pode-se
escrever 2
2GM L2
  
dr
+ 1− 2 − E 2 c2 = 0 (3.4.1)
dλ cr r2
que pode ser colocada em função de φ:
2
E 2 c2 4
  
dr 2GM
+ 1− 2 r2 − r = 0. (3.4.2)
dφ cr L2

Pimeiramente, admitindo que seja possível uma órbita fechada circular para
o feixe (ṙ = r̈ = 0), procura-se sob quais condições isto é satisfeito. Com
Cálculo Tensorial e Relatividade Geral 56

auxílio da equação acima e de (3.2.2)

c2
 
2 2GM 2
φ̇ = 2 1− 2 ṫ , (3.4.3)
r cr

GM 2
φ̇2 = ṫ . (3.4.4)
r3
Igualando (3.4.3) e (3.4.4), verifica-se que a luz pode percorrer trajetórias
circulares se, e somente se, r = 3GM/c2 . Como a condição é extremamente
restritiva, uma perturbação qualquer fará com que o feixe saia de órbita.
Em seguida, vamos analisar a deflexão que sofre um feixe passando pró-
ximo a um corpo massivo. Tomando a equação (3.4.2) escrita em função da
variável u = L2 /(GM r), e derivando-a em φ:

d2 u 3G2 M 2 2
+ u = u. (3.4.5)
dφ2 c2 L2

No caso em que M = 0, a equação acima tem como solução

u = u0 cos (φ − φ0 ) , (3.4.6)

que resulta em uma reta no espaço. Como é esperado que, para M 6= 0, o

desvio seja pequeno [24], vamos aproximar o termo de segunda ordem pela
equação acima, tal que

d2 u 3G2 M 2 2
2
+ u = 2 2
u0 cos2 (φ − φ0 ) . (3.4.7)
dφ cL

Uma das possíveis soluções desta é dada por

G2 M 2 
u = u20 1 + sin2 (φ − φ0 ) .

2 2
(3.4.8)
cL
Cálculo Tensorial e Relatividade Geral 57

Sem perda de generalidade, supondo que φ0 = 0, a solução geral para 1/r é a

soma da última equação com aquela para M = 0, podendo ser escrita como
 
1 1 3GM 1
= cos φ + 2 2 1 + cos (2φ) . (3.4.9)
r r0 2r0 c 3

Supondo que o feixe passa próximo ao corpo gerador de campo gravita-

cional e segue para r → ∞: para o caso em que M = 0 , φ = ±π/2 e o feixe
não sofre deflexão; no caso M 6= 0, se φ = ±π/2 + α, a solução pode ser
escrita como
1 GM
sin α + 2 2 1 + cos2 α = 0


− (3.4.10)
r0 r0 c
aproximando as funções senóides por série de Taylor em primeira ordem,
obtem-se o valor de α:
2GM
α= (3.4.11)
r0 c2
Mas pela equação (3.4.9) vê-se que o desvio é duas vezes este valor, portanto

4GM
δφ = . (3.4.12)
r0 c2

Para se ter noção da ordem de grandeza deste desvio, consideremos o Sol.

O Sol tem raio de aproximadamente 6, 967×108 m, supondo que r0 = 7×108 m,
então
δφ = 8, 423 × 10−6 rad ≈ 1.74” (3.4.13)

O primeiro teste experimental a medir este desvio da luz foi feito em

1919 por Sir Eddington em Sobral, no Ceará, Brasil. No experimento foram
utilizadas dois tipos de lentes, uma detectando 1, 98”, e a outra 0, 93”. Em
1979 testes de maior precisão foram feitos, corroborando fortemente com a
TRG [37,38]. Atualmente, o desvio é a base da técnica de lente gravitacional,
já mencionada. Esta é usada na determinação da distribuição de massa de
Cálculo Tensorial e Relatividade Geral 58

galáxias distantes [39].

3.5 Coordenadas de Lemaître

Dada a necessidade de adotar um novo sistema de coordenadas para des-
crever o espaço-tempo para r 6 rS , nesta seção tais coordenadas serão apre-
sentadas. Antes disto, será apresentado um argumento mais concreto da
necessidade de novas coordenadas. Seja um raio de luz que se move radi-
almente em direção a um corpo - fonte de campo gravitacional. Têm-se
ds2 = dΩ2 = 0 e, portanto

dr  rS 
= ±c 1 − . (3.5.1)
dt r

Supondo que o fóton está em r0 > rS em um tempo inicial t0 :

 
r0 − r rS r0 − rS
t= + ln + t0 . (3.5.2)
c c r − rS

Note que, se r → rS , então t → ∞, sugerindo que a luz se aproxima de

rS , mas nunca o alcança, conclusão claramente absurda. A solução para tal
declaração é que ela reflete a impossibilidade de descrever o espaço-tempo
em r ≤ rS usando as coordenadas de Schwarzschild. Ao mesmo tempo, é
evidente na métrica que para r < rS , a coordenada temporal é tipo espaço e
a radial é tipo tempo, não havendo partículas estacionárias [40].
Quando estas ideias primeiramente vieram a tona, havia uma crença de
que não houvessem corpos cujos raios fossem da ordem do próprio raio de
Schwarzchild, contudo, estudos mostraram que estrelas de Nêutrons podem
colapsar e formar corpos com raio menor que o de Schwarzschild. Tais objetos
são denominados Buracos Negros [41].
Cálculo Tensorial e Relatividade Geral 59

Em r = rS define-se a esfera de Schwarzschild. Para extrair alguma

propriedade do interior da esfera há várias transformações de coordenadas
possíveis, aqui serão usadas as coordenadas de Lemaître (T, R), que formam
um sistema síncrono (g00 = 1) [28]. São definidas como:
Z
F (r)dr
cT = ct +
, (3.5.3)
1 − rrS
Z
dr
R = ct + rS

, (3.5.4)
1− r
F (r)
tal que
1 − rrS
ds2 = c2 dT 2 − F (R)2 dR2 − r2 dΩ2 .


2
(3.5.5)
1 − F (R)
p
Pode-se tomar F (r) = rS /r, tal que a métrica de Schwarzschild, nas coor-
denadas de Lemaître seja

4/3
dR2

2 2 2 3
ds = c dT − h i2/3 − (R − cT ) dΩ2 . (3.5.6)
3 2rS
2rS
(R − cT )

A partir das definições de R e T , obtém-se que

 2/3
3 1/3
r= (R − cT ) rS . (3.5.7)
2

3
Notando que para r = rS , têm-se rS = 2
(R − cT ) e, portanto, a singulari-
dade é eliminada.
Nestas novas coordenadas, supondo novamente um raio de luz se movendo
na direção radial e, usando as relações obtidas para r e rS , têm-se
r
rS
dR = ± cdT. (3.5.8)
r
Cálculo Tensorial e Relatividade Geral 60

que, quando substituída nas definições de R e T :

 r 
rS
dr = ±1 − cdT. (3.5.9)
r

Portanto, nota-se que para raios menores que o de Schwarzchild (r < rS ) é

necessário que dr < 0, e o único movimento possível para o raio de luz é
o decrescente (r → 0). Sendo este o motivo da nomenclatura buraco negro
[41].
Para exemplificar de maneira geométrica, faz-se uso de um diagrama de
cT por R (Figura 3.7): se r é constante, então R − cT também será, produ-
zindo um conjunto de retas. Pelas inequações acima fica evidente que retas
com r > rS estão dentro do cone definido em r = rS (Ponto B), raios de luz
nesta região conseguem escapar para r → ∞; no entando aquelas com r < rS
estão fora do cone (Ponto A), e fadadas a ir em direção à singularidade.
 
 c dT  < 1, r > r
dR S
  (3.5.10)
 c dT  > 1, r < r .
dR S

Dos últimos parágrafos conclui-se que não é possível que uma partícula
dentro do raio de Schwarzschild cruze a superfície em r = rS , que, por este
motivo, também é chamada horizonte de eventos [24]. Esta prisão de partí-
culas impossibilita a verificação do comportamento de partículas na região
interior à superfície de Schwarzschild.
Cálculo Tensorial e Relatividade Geral 61

rS
0 r=
r=
cT

A
B

Figura 3.7: Diagrama de cT por R.

Fonte: Autor(2016).

Os buracos negros ainda são alvo de muitas controvérsias e pesquisa,

sendo importante ressaltar que o buraco negro de Schwarzchild é considerado
o mais simples, visto que não há rotação ou carga. Há aproximadamente 50
anos teorizou-se, a partir do uso de um espaço-tempo curvo dentro da teoria
quântica de campos, que deve haver radiação (radiação Hawking) sendo ex-
pulsa do horizonte de eventos de buracos negros. Apesar de ainda não haver
comprovações experimentais. Este é um dos motivos pelo qual se busca uma
unificação das teorias quânticas e gravitacionais: um maior entendimento do
comportamento de buracos negros. É finalizada esta exposição da TRG com
a expectativa de que, futuramente, áreas como a citada acima possam ser
pesquisadas.
Considerações Finais

Neste trabalho estabeleceu-se uma base matemática e conceitual da Te-

oria da Relatividade Geral, proporcionando o conhecimento necessário para
o início da pesquisa em um leque de áreas. Algumas áreas que podem ser
citadas neste quesito são: a gravitação quântica e a teoria de cordas, que
visam unificar a TRG com teorias quânticas, visto que há problemas na con-
ciliação destas ao se tratar da evolução de buracos negros ou de teoremas
de singularidade [5, 42]; a emissão de radiação por buracos negros, conhecida
por radiação Hawking, que surge ao tomar a teoria quântica de campos em
espaços curvos [2]; teorias f (R), ou teorias de gravitação modificadas, que
visam modificar a ação de Einstein-Hilbert, resultando em termos invarian-
tes de curvatura de ordem maior, que são relevantes no estudo de eventos
próximos a singularidade do buraco negro e do Universo primitivo [22, 43];
entre outras.

62
Referências Bibliográficas

[1] NORTON, J. How Einstein found his field equations: 1912-1915.

Historical studies in the physical sciences, v. 14, n. 2, p. 253–316, 1984.

[2] CARROLL, S. M. Spacetime and Geometry: An introduction to general

relativity. 1. ed. San Francisco: Addison-Wesley Professional, 2003. 513 p.

[3] MINKOWSKI, H. Space and time. In: The Principle of Relativity.

EUA: Dover Publications, Inc, 1952. p. 73–91.

[4] WEATHERALL, J. O. Inertial motion, explanation, and the

foundations of classical spacetime theories. Acesso em: novembro de 2016.
2012, Disponível em: <[Link]

[5] KOX, A. J.; EISENSTAEDT, J. The Universe of General Relativity.

Boston: Birkhäuser, 2006. 383 p.

[6] BERTI, E. et al. Testing general relativity with present and future
astrophysical observations. Classical and Quantum Gravity, v. 32, n. 24, p.
1–179, 2015.

[7] LEMOS, N. A. Mecânica analítica. 2. ed. São Paulo: Editora Livraria

da Física, 2007. 396 p.

[8] GRIFFITHS, D. J.; COLLEGE, R. Introduction to Electrodynamics. 3.

ed. New Jersey: Prentice Hall, 1999. 628 p.

[9] BISHOP, R. L.; GOLDBERG, S. I. Tensor Analysis on Manifolds. New

York: Dover Publications, Inc, 2012. 288 p.

[10] MEDRANO, R. A. Campos vetoriais, espaços lineares e tensores na

física. São José dos Campos, SP: Transtec editorial, 1994. 355 p.

[11] SYNGE, J. L.; SCHILD, A. Tensor calculus. New York: Dover

Publications, Inc, 1969. 336 p.

63
Cálculo Tensorial e Relatividade Geral 64

[12] DAS, A. Tensors: The mathematics of relativity theory and continuum

mechanics. New York: Springer Science & Business Media, 2007. 290 p.

[13] CATTANI, M. Dedução das equações da teoria de gravitação de

Einstein em um curso de graduação. Revista Brasileira de Ensino de Física,
v. 20, n. 1, p. 27–37, Março 1998.

[14] CARMELI, M. Classical Fields: General Relativity and Gauge Theory.

New York: John Wiley & Sons, 1982. 670 p.

[15] NUSSENZVEIG, H. M. Curso de Física Básica: Vol. 4: Ótica,

relatividade e física quântica. 1. ed. São Paulo: Editora Blücher, 1998.
437 p.

[16] LORENTZ, H. A. Electromagnetic phenomena in a system moving

with any velocity smaller than that of light. In: The Principle of Relativity.
EUA: Dover Publications, Inc, 1952. p. 9–34.

[17] EINSTEIN, A. On the electrodynamics of moving bodies. In: The

Principle of Relativity. EUA: Dover Publications, Inc, 1952. p. 35–65.

[18] EINSTEIN, A. The foundations of the general theory of relativity. In:

The Principle of Relativity. EUA: Dover Publications, Inc, 1952. p. 100–164.

[19] EARMAN, J.; GLYMOUR, C. Einstein and Hilbert: Two months in

the history of general relativity. Archive for History of Exact Sciences, v. 19,
n. 3, p. 291–308, 1978.

[20] WILL, C. M. The confrontation between general relativity and

experiment. Living Reviews in Relativity, v. 17, n. 4, 2014. Acesso em:
agosto de 2016. Disponível em: <[Link]

[21] VANKOV, A. A. General Relativity problem of Mercury’s perihelion

advance revisited. Acesso em: dezembro de 2016. 2010, Disponível em:
<[Link]

[22] SOTIRIOU, T. P.; FARAONI, V. f(R) Theories of gravity. Acesso em:

dezembro de 2016. 2010, Disponível em: <[Link]
1726>.

[23] CASTELLANI, O. C. Discussão dos conceitos de massa inercial e massa

gravitacional. Revista Brasileira de Ensino de Física, v. 23, p. 356–359,
Setembro 2001.
Cálculo Tensorial e Relatividade Geral 65

[24] PAPAPETROU, A. Lectures on General Relativity. Dordrecht-Holanda:

D. Reidel Publishing Company, 2012. 218 p.
[25] REFSDAL, S.; BONDI, H. The gravitational lens effect. Monthly
Notices of the Royal Astronomical Society, v. 128, n. 4, p. 295–306, Agosto
1964.
[26] CHANDRASEKHAR, S. On the “derivation” of Einstein’s field
equations. American Journal of Physics, v. 40, n. 2, p. 224–234, 1972.
[27] OLIVEIRA, T. B. R. d. F. Teorias f (R) de Gravidade na Formulação
de Palatini. 103 p. Dissertação (Mestrado) — Universidade Federal do Rio
Grande do Norte, Natal - RN, Julho 2010.
[28] LANDAU, L. D.; LIFSHITZ, E. M. The Classical Theory of Fields:
Course of theoretical physics Volume 2. 4. ed. Oxford: Butterworth-
Heinenann, 2013. 402 p.
[29] CORRY, L.; RENN, J.; STACHEL, J. Belated decision in the
Hilbert-Einstein priority dispute. Science, American Association for the
Advancement of Science, v. 278, n. 5341, p. 1270–1273, Novembro 1997.
[30] BRADING, K. A.; RYCKMAN, T. A. Hilbert’s ‘Foundations of
Physics’: Gravitation and electromagnetism within the axiomatic method.
Studies in History and Philosophy of Science Part B: Studies in History and
Philosophy of Modern Physics, v. 39, n. 1, p. 102–153, Janeiro 2008.
[31] FRIEDMANN, A. 125. On the curvature of space. Zeitschrift für
Physik, v. 10, p. 377–386, 1922.
[32] EINSTEIN, A. Cosmological considerations on the general theory of
relativity. In: The Principle of Relativity. EUA: Dover Publications, Inc,
1952. p. 175–188.
[33] STRAUMANN, N. The history of the cosmological constant
problem. Acesso em: setembro de 2016. 2008, Disponível em: <https:
//[Link]/abs/gr-qc/0208027>.
[34] FOSTER, J.; NIGHTINGALE, J. D. A short course in General
Relativity. 3. ed. New York: Springer Science & Business Media, 2010.
292 p.
[35] SCHWARZSCHILD, K. On the gravitational field of a mass point
according to einstein’s theory. Sitzungsberichte der Königlich Preussischen
Akademie der Wissenschaften, Berlin, v. 1, p. 189–196, 1916.
Cálculo Tensorial e Relatividade Geral 66

[36] BIESEL, O. The precession of mercury’s perihelion. Acesso em:

setembro de 2016. 2008, Disponível em: <[Link]
edu/~morrow/papers/>.

[37] DYSON, F. W.; EDDINGTON, A. S.; DAVIDSON, C. A determination

of the deflection of light by the sun’s gravitational field, from observations
made at the total eclipse of may 29, 1919. Philosophical Transactions of
the Royal Society of London A: Mathematical, Physical and Engineering
Sciences, The Royal Society, v. 220, n. 571-581, p. 291–333, Janeiro 1920.

[38] KENNEFICK, D. Testing relativity from the 1919 eclipse—a question

of bias. Physics Today, v. 62, n. 3, p. 37–42, Março 2009.

[39] REFSDAL, S. On the possibility of determining Hubble’s parameter

and the masses of galaxies from the gravitational lens effect. Monthly
Notices of the Royal Astronomical Society, v. 128, n. 4, p. 307–310, 1964.

[40] CHOW, T. L. Gravity, Black Holes, and the Very Early Universe:
An introduction to general relativity and cosmology. New York: Springer
Science & Business Media, 2008. 300 p.

[41] FROLOV, V.; NOVIKOV, I. Black Hole Physics: Basic concepts and
new developments. Dordrecht-Holanda: Klwer Academid Publishers, 2012.
795 p.

[42] KIEFER, C. Quantum gravity (general) and applications. In:

Compendium of Quantum Physics. [S.l.]: Springer, 2009. p. 565–572.

[43] FELICE, A. D.; TSUJIKAWA, S. f(R) Theories. Living Reviews in

Relativity, v. 13, n. 3, 2010. Acesso em: novembro de 2016. Disponível em:
<[Link]
Apêndice A

Derivação da Equação da
Geodésica.

Sejam xν = xν (u, v) uma família de infinitas curvas, com o parâmetro v

fixo sobre cada curva, e u variando (Figura A.1).

Figura A.1: Família de curvas xµ

Fonte: Autor(2016).

O elemento de linha pode ser escrito como

ds2 = gµν pµ pν du2 , (A.0.1)

com pµ = dxµ /du. Seja L o arco de curva, tal que

Z B
L= ds (A.0.2)
A

67
Cálculo Tensorial e Relatividade Geral 68

e, L(xµ , pµ ) = gµν pµ pν . Então:

Z u2 √
L= Ldu, (A.0.3)
u1

derivando em relação a v:
" √ √  λ #
u2 λ
∂ L ∂x ∂ L ∂ ∂x
Z
dL
= λ
+ du. (A.0.4)
dv u1 ∂x ∂v ∂pλ ∂v ∂u

Reescrevendo o último termo na integral, obtém-se

( √ √ ! √ ! )
u2
∂ L ∂xλ ∂ L ∂xλ ∂ L ∂xλ
Z
dL ∂ ∂
= λ
+ λ
− du
dv u1 ∂x ∂v ∂u ∂p ∂v ∂u ∂pλ ∂v
(A.0.5)
que com uma integração por partes resulta em

" √ #u2 Z 2 " √ √ !#

u
dL ∂ L λ ∂ L ∂ ∂ L
δv = δL = λ
δx + λ
− λ
δxλ du, (A.0.6)
dv ∂p 1 u1 ∂x ∂u ∂p
u

o primeiro termo é nulo, pois os extremos são mantidos fixos e δL deve ser
nulo. Portanto,
√ ! √
d ∂ L ∂ L
− =0 (A.0.7)
du ∂pλ ∂xλ

que é a equação contraída da geodésica. Para a obtenção do resultado ex-

plícito, primeiramente é necessário expandir a equação acima, note que v se
torna uma constante da curva e não mais um parâmetro:
 
d ∂L ∂L 1 dL ∂L
− = , (A.0.8)
du ∂pλ ∂p λ 2L du ∂pλ
Cálculo Tensorial e Relatividade Geral 69

considerando u = s e dL/ds = 0:
 
d ∂L ∂L
− = 0. (A.0.9)
ds ∂pλ ∂xλ

Substituindo a expressão para L e simplificando alguns termo, reescrevemos

a equação (A.0.9) como

d ∂gµν µ ν
(2gλµ pµ ) − p p = 0, (A.0.10)
ds ∂xλ

que, ao se expandir o primeiro termo, é

∂gλµ µ ν dpµ 1 ∂gµν µ ν

p p + gλµ − p p = 0. (A.0.11)
∂xν ds 2 ∂xλ

Escrevendo  
∂gλµ µ ν 1 ∂gλµ ∂gλν
ν
p p = + p µ pν , (A.0.12)
∂x 2 ∂xν ∂xµ
a equação da geodésica fica

d 2 xµ dxµ dxν
gλµ + Γλµν = 0, (A.0.13)
ds ds ds

que, multiplicada pelo tensor recíproco, retorna a equação da geodésica em

forma explícita
d 2 xλ λ dxµ dxν
+ Γ µν = 0. (A.0.14)
ds2 ds ds
Apêndice B

Expressão para o Divergente

Pela definição de derivada covariante, o divergente de um tensor contra-

variante de ordem 1 é

∂V µ
∇µ V µ = + Γµ µν V ν (B.0.1)
∂xµ

∂gµλ
com Γµ µν = 21 g µλ ∂xν
. O determinante do tensor métrico é dado por

g = αβγ···N g1α g2β · · · gN n , (B.0.2)

na qual αβ···N é o tensor de Levi-Civita, que possui a seguinte propriedade:

Λαβ···N = (−1)Λ+1 αβ···N . Portanto, a derivada do determinante pode ser
escrita como

∂g ∂g1α ∂gnN
µ
= µ
(−1)1+α βγ···N g2β · · · gnN +· · ·+ µ
(−1)N +n αβ···N −1 g1α · · · gm(N −1) .
∂x ∂x ∂x
(B.0.3)

70
Cálculo Tensorial e Relatividade Geral 71

O cofator do elemento gµν é definido como ∆gΛµ = (−1)Λ+µ g1α g2β · · · gnN ,
mas pode ser escrito como ∆gµν = g µν g, portanto

∂g ∂gλω
µ
= g µν µ g. (B.0.4)
∂x ∂x

Voltando para a equação do divergente

∂V µ 1 V ν ∂g
∇µ V µ = + (B.0.5)
∂xµ 2 g ∂xν

que, claramente, pode ser escrita como

1 ∂ √
∇µ V µ = √ µ
(V µ g) , (B.0.6)
g ∂x

como esperado.
 Apêndice C

Simulação das Órbitas

Para a simulação das órbitas foram escritos dois programas em linguagem

C. Como há mudança em apenas duas linhas do código (linhas 19 e 20),
abaixo segue o código completo para o caso relativístico e, em seguida, a
alteração feita para o caso clássico.

Simulação das órbitas para o potencial Relativístico

1 # include < stdio .h >
2 # include < stdlib .h >
3 # include < math .h >
4

5 # define G 1
6 # define M 1
7 # define L 8
8 # define c 1
9 # define r0 55
10 # define p0 M_PI
11 # define vr0 0.0
12 # define R 1e -5

72
 Cálculo Tensorial e Relatividade Geral 73

13 # define t0 0.0
14 # define tf 12000
15 # define dt 1e -5
16 # define mult ( int ) ( tf /( dt *10000) )
17

18 double f ( double r ) {
19 return ( -1.0) * G * M /( r * r ) + L * L /( r * r * r ) -
3* G * M * L * L /( c * c * r * r * r * r ) ;
20 }
21

22 double g ( double vr ) {
23 return vr ;
24 }
25

26 double h ( double r ) {
27 return L /( r * r ) ;
28 }
29

30 int main ( int argc , char ** argv ) {

31

32 int counter = 0;
33 double r = r0 , p = p0 , t = t0 , vr = vr0 ;
34 double f1 , g1 , h1 ;
35 FILE * outp = fopen ( " rel . dat " , " w " ) ;
36

37 while (t < tf ) {
38

39 if ( counter % mult == 0) {
40 fprintf ( outp , " % e ␣ % e ␣ % e \ n " , r * cos ( p ) , r * sin ( p ) , t ) ;
 Cálculo Tensorial e Relatividade Geral 74

41 }
42 f1 = f ( r ) ;
43 g1 = g ( vr ) ;
44 h1 = h ( r ) ;
45

46 r += dt * g1 ;
47 vr += dt * f1 ;
48 p += dt * h1 ;
49 t += dt ;
50 counter ++;
51 }
52

53 fclose ( outp ) ;
54 return 0;
55 }

Mudança feita para o caso clássico

1 double f ( double r ) {
2 return ( -1.0) * G * M /( r * r ) + L * L /( r * r * r ) ;
3 }