,

%,

%
,

hhP
XP
X
hH
 Instituto de Fsica Teorica
((
( (

X
h
H
P
%
,

Q
l
eQ
@
l
e
@
Q
IFT Universidade Estadual Paulista
@l

TESE DE DOUTORADO
IFTD.004/12

Topicos em Teoria Quantica de Campos

Eslley Scatena Goncales

Orientador

Antonio Jose Accioly

Agosto de 2012
When I heard the learnd astronomer,
When the proofs, the figures, were ranged in columns before me,
When I was shown the charts and diagrams, to add, divide, and measure them,
When I sitting heard the astronomer where he lectured with much applause in
the lecture-room,
How soon unaccountable I became tired and sick,
Till rising and gliding out I wanderd off by myself,
In the mystical moist night-air, and from time to time,
Lookd up in perfect silence at the stars.

When I Heard the Learnd Astronomer, de Walt Whitman.

i
 Agradecimentos

Quatro anos e muito tempo para que um projeto como esse se desenvolva sem que sofra
imensas modicacoes, sejam elas por vontade propria, sejam por conta das equacoes
que assim quiseram. E por mais que eu queira dizer que este trabalho terminou, sei

que sempre que rele-lo vou querer alterar alguma coisa, um ndice trocado ou uma
discussao que poderia ter sido mais aprofundada. Por outro lado, essa incompletude

se reete tambem no proprio corpo da Ciencia, que esta sempre em desenvolvimento.

Agradeco, entao, aos que me ajudaram a concluir parte dessa empreitada inacaba-
vel comecando pela minha famlia, em especial meu irmao, David, que sempre apoiou

as minhas decisoes. Sou grato aos amigos que nao participaram diretamente deste tra-
balho, mas certamente contribuiram para a minha sanidade mental durante o processo

Fabio, Henrique, Linneu, Marcelo, Italo, Leandro, Nara, Mor, Rone, Claudio, Carol
e Max aos que, de alguma forma, contribuiram diretamente Rodrigo Turcati, Jef-

ferson Morais, Bruno Dias, Carlos Bonilla (pelas discussoes fsicas) e Tamara Prior
(pela revisao do texto na primeira versao) e, tambem, a Tamara West, que tem me
inspirado, incentivado e motivado cada vez mais.

Durante o desenvolvimento deste trabalho tive o prazer de contar com o apoio

dos pesquisadores do CBPF Jose Helayel-Neto, Sebastiao Alves Dias e Sergio Duarte,

sempre dispostos a ajudar e abertos a discussoes. Obrigado.

Sou imensamente grato ao professor Antonio Accioly que, alem de excelente orien-

tador, mostrou-se um grande amigo nesses mais de 6 anos em que trabalhamos juntos.
Por m, agradeco ao IFT pela oportunidade de desenvolver este trabalho e a CA-
PES, pelo apoio nanceiro.

ii
 Resumo

Dois topicos importantes e atuais em teoria quantica de campos, a saber, Teorias com
Derivadas de Ordem mais Alta e Determinacao de Limites para a Massa do Foton,

sao considerados. Como prototipos para as teorias de ordem mais alta selecionamos a
Eletrodinamica de Lee e Wick e a Gravitacao Massiva em 3D. Ambas as teorias sao de

relevancia especial para a fsica contemporanea. A primeira por constituir o alicerce

sobre o qual o Modelo Padrao de Lee e Wick se assenta, enquanto que a outra e a unica

teoria de gravitacao com derivadas de ordem mais alta conhecida ate o momento cuja
versao linearizada e unitaria em nvel de arvore. Novas e interessantes propriedades
relativas as mencionadas teorias sao estudadas em detalhe. Finalmente, limites classicos

e quanticos para a massa do foton sao estimados.

Palavras-chave: Eletrodinamica de Lee e Wick; Gravitacao Massiva em 3D; Foton

Massivo.
Areas de conhecimento: Teoria Quantica de Campos; Gravitacao.

iii
 Abstract

Two important and current topics in quantum eld theory, namely, Higher-Derivative
Theories and Bounds on the Photon Mass, are considered. As prototypes for the

higher-order theories we have singled out Lee-Wick Finite Electrodynamics and Mas-
sive Gravity in 3D. Both theories are of special relevance to contemporary physics. The

former because it is the base on which the Lee-Wick Standard Model rests, while the
other is the only higher-order gravity theory known up to now whose linearized version

is tree-level unitary. Novel and interesting properties related to both theories are analy-
zed in detail. Finally, classical and quantum bounds on the photon mass are estimated.

iv
Sumario

Sumario vii

Introducao 1
Referencias Bibliogracas . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6

I Teorias com Derivadas de Ordem Mais Alta 8

Referencias Bibliogracas . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13

1 Explorando a Eletrodinamica Finita de Lee e Wick 14

1.1 Motivacao . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
1.2 Aspectos gerais do modelo de Lee e Wick . . . . . . . . . . . . . . . . . 16

1.2.1 Analise de um toy model escalar . . . . . . . . . . . . . . . . . . 20

1.3 O tensor energia-momento de Lee e Wick . . . . . . . . . . . . . . . . . 25

1.4 Propagacao ondulatoria na eletrodinamica de Lee e Wick . . . . . . . . 28

1.5 Limites sobre a massa da partcula pesada de Lee e Wick . . . . . . . . 31
1.5.1 Um limite classico sobre a massa da partcula pesada de LW . . 31

1.5.2 Um limite quantico sobre a massa da partcula pesada de LW . 35

1.6 O potencial nao-relativstico de LW . . . . . . . . . . . . . . . . . . . . 38

v
 1.6.1 Mecanica quantica e o limite nao-relativstico da TQC . . . . . 38
1.6.2 Analise comparativa entre os potenciais de LW e Coulomb . . . 44

1.7 Monopolos magneticos . . . . . . . . . . . . . . . . . . . . . . . . . . . 44

1.8 Discussao . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 48

Referencias Bibliogracas . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50

2 Explorando a Gravitacao Massiva em 3D 54

2.1 Motivacao . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 54
2.2 O propagador da gravitacao massiva em (2+1)D . . . . . . . . . . . . . 56

2.2.1 Calculo de g, , R , R e R em primeira ordem em . . 57
2.2.2 Forma gama-gama da lagrangiana de Einstein . . . . . . . . . . 60
2.2.3 A forma bilinear da lagrangiana da 3DHDG . . . . . . . . . . . 61

2.3 Unitariedade em nvel de arvore do modelo 3DHDG . . . . . . . . . . . 65

2.4 Potencial inter-partculas . . . . . . . . . . . . . . . . . . . . . . . . . . 67

2.5 Desvio da Luz . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 72

2.6 Analise dos resultados obtidos no contexto dos modelos BHT e 3DHDG 78

2.6.1 O potencial nao-relativstico para o modelo 3DHDG . . . . . . . 78

2.6.2 Deexao da luz no contexto da 3DHDG . . . . . . . . . . . . . 80
2.6.3 Um estudo comparativo . . . . . . . . . . . . . . . . . . . . . . 81

2.7 Discussao . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 83
Referencias Bibliogracas . . . . . . . . . . . . . . . . . . . . . . . . . . . . 86

Apendice A Prescricao para o calculo do propagador do graviton 89

Apendice B Equacoes de campo e o gauge de Teyssandier 92

B.1 Equacoes de campo do modelo BHT . . . . . . . . . . . . . . . . . . . 92

B.2 Equacoes de campo linearizadas no gauge de Teyssandier . . . . . . . . 95

vi
 B.3 Solucao das equacoes de campo para uma partcula puntual . . . . . . 98

Apendice C Um resultado util para checar a unitariedade da 3DHDG 101

II Determinacao de Limites para a Massa do Foton 103

Referencias Bibliogracas . . . . . . . . . . . . . . . . . . . . . . . . . . . . 106

3 A Massa do Foton e a Interferometria de Linha de Base Extensa 107

3.1 Motivacao . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 107

3.2 Combinando RG, QED Massiva e VLBI para limitar a massa do foton . 111
3.2.1 Calculo da lagrangiana de Interacao . . . . . . . . . . . . . . . . 112

3.2.2 Calculo da funcao de vertice . . . . . . . . . . . . . . . . . . . . 113

3.2.3 Calculo da secao de choque diferencial nao-polarizada . . . . . . 115

3.3 Relacao entre massa, frequencia e parametro de deexao para o foton . 118
3.4 Encontrando limites gravitacionais superiores para a massa do foton . . 123
3.5 Discussao . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 128

Referencias Bibliogracas . . . . . . . . . . . . . . . . . . . . . . . . . . . . 131

Eplogo 135

Referencias Bibliogracas . . . . . . . . . . . . . . . . . . . . . . . . . . . . 138

vii
Lista de Figuras

1.1 Partcula pesada de Lee e Wick decaindo em outras duas mais leves . . . . 23
1.2 Relacoes de dispersao para a eletrodinamica de LW. . . . . . . . . . . . 29

1.3 Geometria do experimento de Plimpton e Lawton. . . . . . . . . . . . . . . 33

1.4 Correcao de vertice para um eletron espalhado por um campo externo. . . . 35

1.5 Contribuicao de menor ordem para o espalhamento elastico e e+ . . . . 39

1.6 Vertice relevante para a QED de LW. . . . . . . . . . . . . . . . . . . . . 40
e2
1.7 O potencial U (em unidades de 4
), como funcao da distancia r . . . . 45

2.1 Vertice relevante para a interacao gravitacional de bosons escalares . . . . . 69

2.2 A forca gravitacional de curto alcance relacionada ao modelo de BHT. . . . 72

2.3 Geometria da deexao da luz. . . . . . . . . . . . . . . . . . . . . . . . . 76

2.4 Possveis regimes gravitacionais preditos pela 3DHDG. . . . . . . . . . . . 79

2.5 Energia potencial gravitacional para os modelos 3DHDG e BHT . . . . 82

3.1 Interacao entre um foton massivo e um campo gravitacional externo. . 113

d
( d )
3.2 d
( d )s=0
como funcao do angulo de espalhamento . . . . . . . . . . . . . 122

viii
Lista de Tabelas

A.1 Tabela multiplicativa para os operadores de Barnes-Rivers . . . . . . . 90

3.1 Secoes de choque para o espalhamento de diferentes partculas. . . . . . 120

3.2 Limites superiores para a massa do foton. . . . . . . . . . . . . . . . . . 124

3.3 Limites para massa do foton obtidos via dispersao da velocidade da luz. 130

ix
Introducao

Podemos dizer, em poucas palavras, que a teoria quantica de campos (QFT) nasceu
da necessidade da fusao da relatividade especial com a mecanica quantica (QM); em

contrapartida, a necessidade de conciliar a relatividade geral com a QM deu origem a

teoria de cordas. Acredita-se que um tratamento correto do campo gravitacional sob o

ponto de vista da QFT, tratamento este por hora ainda desconhecido, devera reproduzir
a teoria da relatividade geral de Einstein no limite de baixas energias ou, de modo mais

geral, o Sistema Einstein-Yang-Mills- Dirac. Em outras palavras, e provavel que a QFT

nada mais seja que o limite de baixas energias relativo a uma teoria de campos efetiva
de uma teoria mais fundamental, tal como a teoria de supercordas.

A QFT e na realidade um conjunto de ideias e ferramentas que combina harmo-

niosamente tres dos maiores temas da fsica moderna: a teoria quantica, o conceito

de campo e o princpio da relatividade. Assim sendo, ela nao so fornece um solido

alicerce para a fsica das partculas elementares contemporanea, mas e tambem um

repositorio de ferramentas que sao essenciais para a fsica nuclear e a fsica atomica,
bem como para a fsica da materia condensada e a astrofsica. Alem disso, esta teoria
tem permitido que se ediquem novas pontes entre a fsica e a matematica. Primitiva-

mente, a motivacao para o estudo da QFT jazia na esperanca de que ela viesse lancar

1
uma nova luz sobre nossa maneira de enxergar as partculas fundamentais da materia
e suas interacoes. De fato, uma estrutura que incorporava, por um lado, a QM (que

tinha sido tao bem sucedida no que tange a resolucao de um grande numero de pro-
blemas da fsica atomica no incio do seculo passado) e, por outro, a teoria de campos

(a linguagem utilizada para descrever o surpreendente quadro da realidade desvelado

por Faraday, Maxwell e Hertz), certamente deveria fornecer alguma indicacao conavel

sobre a natureza fundamental da materia. Na realidade, os resultados superaram em

muito as expectativas. Para nao nos alongarmos muito analisaremos, de passagem,
apenas a QED, o primeiro lho do bem sucedido casamento da relatividade especial

com a QM. Esta teoria prediz, para citar apenas um dos seus sucessos, o momento
magnetico anomalo do eletron corretamente com onze casas decimais. Esta fantastica

concordancia entre os resultados previstos teoricamente e os obtidos experimentalmente

levaram Ryder [1] a perguntar: O que mais se poderia querer de uma teoria fsica?.

A QFT tem tambem sido estudada sob o ponto de vista axiomatico. E crenca
geral que a axiomatizacao e uma especie de polimento, que e aplicado a uma area da
ciencia apos ela ter sido, no que concerne a todos os propositos praticos, completada.

Obviamente, isto nao e verdade nem mesmo no contexto da matematica pura. Em f-

sica teorica, desde o tempo de Newton, o metodo axiomatico tem servido nao somente

para a sistematizacao dos resultados previamente obtidos, mas tambem para auxiliar
na descoberta de novos resultados. O signicado do termo axiomatizacao como em-

pregado na fsica vem gradualmente evolundo em direcao a uma maior generalizacao

e abstracao. Como os resultados experimentalmente vericaveis provenientes da QFT
axiomatica nao sao muitos [2], o metodo axiomatico ocupa um lugar comparativamente

pequeno no contexto da teoria das partculas elementares de nossos dias. E digno de

nota que resultados rigorosos em QFT tem a grande vantagem de possurem um tempo

2
de vida um tanto quanto mais longo quando comparados com a maioria dos trabalhos
da moda em fsica das partculas elementares.

Durante as duas ultimas decadas a QFT passou a ser um topico cada vez mais dis-
cutido no ambito da losoa da fsica. Esta teoria constitui um topico muito atrativo

para os losofos sob o ponto de vista metodologico, semantico bem como ontologico. E
de causar especie, no entanto, que a discussao dos fundamentos conceituais do domnio

quantico tenha sido sempre primazia da QM e nao da QFT. De fato, a QFT e, num
certo sentido, muito mais abrangente que a QM, sendo tambem, ao contrario desta ul-
tima, relativisticamente invariante. Existem pelo menos duas razoes para se favorecer,

no que concerne as reexoes conceituais, a QM em detrimento da QFT. A primeira e

que durante muito tempo acreditou-se que os problemas losocos de maior peso, em

particular o problema da medida, ja tinham aparecido na QM, de modo que uma ana-
lise conceitual da QFT parecia ser redundante. Acreditava-se tambem que um exame

mais acurado da QFT so iria ajudar a obscurecer a visao de seus aspectos centrais,
ja que a QFT e muito mais complexa que a QM padrao, inclusive no que se refere as
diculdades matematicas. Uma segunda razao para se negligenciar a QFT se baseava

no argumento que ela nao tinha ainda atingido o status de teoria completa e consis-
tente. Alem disso, a inexistencia de uma teoria de gravitacao constitua tambem um

fator bastante agravante. Certamente, nao se pode negar que a incorporacao da quarta
forca fundamental (gravitacao) devera levar forcosamente a profundas modicacoes na

QFT como um todo; a versao corrente da QFT nada mais que uma teoria preliminar.
No entanto, se fossemos esperar por uma versao completa da QFT e bastante provavel
que uma analise losoca da mesma nao comecasse nunca. Na verdade, reexoes inter-

pretativas sobre os fundamentos da fsica e suas inconsistencias nao deveriam jamais

ser evitadas, pois elas ajudam na busca da teoria nal. Interessantes discussoes sobre

3
a losoa da QFT podem ser encontradas em [37].
Nosso objetivo neste trabalho e analisar no contexto da QFT dois topicos bastante

relevantes e atuais: Teorias com Derivadas de Ordem Mais Alta e Determinacao de

Limites para a Massa do Foton.

Cabe aqui uma indagacao simples e natural: Qual a razao que nos levou a sele-
cionar entre os inumeros e interessantes topicos que compoem a QFT precisamente

os acima mencionados? A resposta a esta questao e bastante simples: Pretendamos,

por um lado, entender melhor a ideia em voga de que teorias de ordem mais alta se
sao renormalizaveis nao sao unitarias e, se sao unitarias sao assombradas por fantas-

mas. Seria isto uma verdade incontestavel ou uma simples conjectura? Se um mito,
como desmitica-lo?; por outro, queramos encontrar limites para a massa do foton

que nao fossem baseados apenas em calculos de ordem grandeza, como sao em geral os
resultados encontrados na literatura, mas sim em resultados provenientes de calculos

rigorosos realizados no ambito da QFT e onde utilizaramos como dados de entrada as

mais recentes medidas experimentais.
Por razoes didaticas vamos desdobrar o trabalho em duas partes:

Parte I: Teorias com Derivadas de Ordem Mais Alta.

Parte II: Determinacao de Limites para a Massa do Foton.

Na Parte I nos ocupamos com as Teorias Eletromagneticas de Ordem Superior e com as

Teorias Gravitacionais com Derivadas de Ordem Mais Alta; enquanto que na Parte II
nos devotamos a determinacao de limites classicos e quanticos para a massa do foton.

Iniciamos a Parte I explorando a Eletrodinamica Finita de Lee e Wick [8, 9] que,

4
em sua versao atual, e denida pela lagrangiana

1 1
L = F F F 2F ,
4 4m2

onde m e um parametro com dimensao de massa. Limites classicos e quanticos sao

determinados para o parametro m . A propagacao ondulatoria no contexto desta

eletrodinamica de ordem mais alta e entao analisada. Discutimos tambem a interessante

questao relacionada a possibilidade de coexistencia de uma partcula massiva e uma

carga magnetica no ambito da mencionada eletrodinamica [1012]. Em sequencia nos

voltamos para o estudo da Gravitacao Massiva em 3D, que e denida pela lagrangiana

 2R 2  2 3 
L= g 2 + 2 2 R R2 ,
m2 8
onde m2 e um parametro com dimensao de massa. Propriedades interessantes deste
modelo, tais como a presenca de uma forca gravitacional de curto alcance no limite

nao-relativstico e a existencia de um angulo de deexao gravitacional que depende

do parametro de impacto, sao discutidas em detalhe [13, 14]. Limites para a massa
do foton sao obtidos na Parte II por dois caminhos diferentes: via uma combinacao

da Relatividade Geral, QED Massiva e VLBI (Very Long Baseline Interferometry)

e, posteriormente, utilizando o momento magnetico anomalo do eletron calculado no

contexto da QED Massiva [1518]. Estes resultados foram includos no PDG [19].
Encerramos o trabalho analisando a importante e interessante questao, ainda em

aberto, concernente a uma possvel incompatibilidade entre renormalizabilidade e uni-

tariedade em teorias com derivadas de ordem superior. Apresentamos tambem uma
discussao sobre possveis extensoes dos resultados aqui elaborados.

5
Referencias Bibliogracas

[1] L.H. Ryder, Quantum Field Theory, Cambridge University Press (1985).

[2] N. Bogulubov, A. Logunov, and I. Todorov, Introduction to Axiomatic Quantum

Field Theory, Reading (1975).

[3] H. Brown and R. Harre, eds, Philosophical Foundations of Quantum Field Theory,

Clarendon Press (1988).

[4] T. Cao, ed, Conceptual Foundations of Quantum Field Theory, Cambridge Uni-
versity Press (1999).

[5] M. Redhead, Quantum Field Theory for Philosophers, in Asquith and Nickles,
eds., 1983, Proceedings of the Biennial meeting of the Philosophy of Science As-

sociation:PSA 1982, East Lansing, MI: Philosophy of Science Association, vol. 2,

pp. 57-99.

[6] M. Redhead, A philosopher looks at quantum eld theory, in Brown and Harre

1988, pp. 9-23.

[7] M. Redhead, Some philosophical aspects of particle physics, Studies in history and

Philosophy of Science 11 279 (1980).

6
 [8] T. Lee and G. Wick, Nucl. Phys. B 9, 209 (1969).

[9] T. Lee and G. Wick, Phys. Rev. D 2, 1033 (1970).

[10] A. Accioly and E. Scatena, Mod. Phys. Lett. A 25, 269 (2010).

[11] A. Accioly, P. Gaete, J. Helayel-Neto, E. Scatena, and R. Turcati, Mod. Phys.

Lett. A 26 1985 (2011).

[12] A. Accioly, P. Gaete, J. Helayel-Neto, E. Scatena, and R. Turcati, ar-

Xiv:1012.1045v2 [hep-th].

[13] A. Accioly, J. Helayel-Neto, E. Scatena, J. Morais, R. Turcati, and B. Pereira-Dias,

Class. Quantum Grav. 28 225008 (2011).

[14] A. Accioly, J. Helayel-Neto, J. Morais, E. Scatena, and R. Turcati, Phys. Rev. D

83 104005 (2011).

[15] A. Accioly, J. Helayel-Neto and E. Scatena, Phys. Rev. D 82 065026 (2010).

[16] A. Accioly, J. Helayel-Neto, R. Turcati, J. Morais and E. Scatena, Classical Quan-

tum Gravity 27 205010 (2010).

[17] A. Accioly, J. Helayel-Neto and E. Scatena, Phys. Lett. A 374 3806 (2010). This
essay was awarded an honorable mention in the 2010 Essay Competition of the
Gravity Research Foundation.

[18] A. Accioly, J. Helayel-Neto and E. Scatena, Int. J. Mod. Phys. D 19 2393 (2010).

[19] J. Beringer et al. (Particle Data Group), Phys. Rev. D 86 010001 (2012).

7
 Parte I

Teorias com Derivadas de Ordem

Mais Alta

8
 Lagrangianas de ordem superior foram estudadas inicialmente como uma mera cu-
riosidade matematica. No entanto, com o nascimento da Mecanica Quantica e a subse-

quente quantizacao dos campos, estes objetos passaram a serem utilizados como uma
especie de remedio que se acreditava ser ecaz para a cura dos males que aigiam

as teorias fsicas vigentes no nal decada de 40 do seculo passado. Podolsky [1], por
exemplo, deixa isto bem claro, quando apresenta suas razoes para a construcao de uma

eletrodinamica de ordem superior:

Se supusermos que as equacoes da eletrodinamica sao derivaveis a partir de uma
lagrangiana, L, e quisermos preservar o carater linear das equacoes de campo (o prin-

cpio de superposicao) a m de tornarmos a quantizacao facil, entao, a nao ser que

estejamos preparados para introduzir novos campos, a unica maneira de generalizar a

teoria de Maxwell parece ser permitir que a lagrangiana do campo contenha termos
envolvendo derivadas dos campos E e B. Obtemos, entao, como equacoes de campo,

equacoes diferenciais parciais de ordem superior a usual segunda. Longe de ser um

problema, isso parece ser o que era necessario. Pois, os varios metodos propostos de
truncar efeitos de altas frequencias parecem indicar que as derivadas superiores, que se

tornam importantes para as altas frequencias, nao sao tratadas de forma apropriada
pelas equacoes usuais de segunda ordem.

A eletrodinamica de Podolsky e denida pela lagrangiana

1 l2
L = F F + F F ,
4 2

onde l e um parametro com dimensao de comprimento.

Gracas a este trabalho de Podolsky, as teorias de ordem mais alta tornaram-se boas

9
candidatas a sanar problemas devidos aos efeitos de alta frequencia. Estes efeitos nada
mais sao que as divergencias geradas nos calculos de quantidades fsicas mensuraveis

tais como secoes de choque, decaimentos, massas, constantes de acoplamento devidas

ao comportamento extremamente singular dos campos quanticos em distancias muito

pequenas.
No incio da decada de 50, no entanto, Pais e Uhlenbech [2] demonstraram que as

teorias de ordem superior eram assombradas por fantasmas. Apesar de tais teorias nao
terem sido descartadas por completo, o entusiasmo pelas mesmas arrefeceu.
No incio dos anos setenta o interesse pelas teorias de ordem superior recrudesceu

gracas aos trabalhos de Lee e Wick [3, 4] que procuraram contornar o problema da
nao-unitariedade da eletrodinamica por eles proposta via uma serie de procedimentos

engenhosos. E importante notar que a eletrodinamica de Lee e Wick foi construda

anteriormente a QED, com o to de entender a diferenca de massas dos mesons. Se

compararmos as lagrangianas que denem, respectivamente, as teorias de Podolsky e

a de Lee e Wick em sua versao atual, vemos que elas sao equivalentes. Por outro lado,
tudo parece indicar que Lee e Wick nao estavam a par do trabalho de Podolsky, pois a

lagrangiana que eles usaram em seu trabalho nao e exatamente a que hoje denominamos
lagrangiana de Lee e Wick. A lagrangiana proposta originalmente por eles era aquela

resultante da substituicao, na lagrangiana de Maxwell, do campo usual do foton sem

massa A pelo campo complexo

= A + iB ,

onde B e um campo bosonico massivo associado a uma metrica negativa.

Com o advento da QED, a eletrodinamica de Lee e Wick entrou numa fase de

10
ostracismo que durou cerca de quase quatro decadas. Em 2008 esta eletrodinamica
voltou a ser alvo de intensas investigacoes em razao do Modelo Padrao de Lee e Wick

proposto por Grinstein, OConnell e Wise [5] ter utilizado a mesma como paradigma.
O interesse pela eletrodinamica de Lee e Wick, no entanto, ainda nao arrefeceu pois

acredita-se que um melhor conhecimento deste modelo venha lancar novas luzes sobre
o funcionamnto do Modelo Padao de Lee e Wick.

Por outro lado as teorias de gravitacao com derivadas de ordem mais alta foram
sugeridas por Weyl [6] e Eddington [7]. Estas teorias nada mais eram que simples gene-
ralizacoes da relatividade geral obtidas, grosso modo, acrescentando-se a lagrangiana

de Einstein os escalares R2 , R R e R R . Uma interessante discussao sobre

estas teorias pode ser encontrada em Havas [8]. Vericou-se mais tarde que, devido

ao teorema da Gauss-Bonnet, somente era necessario considerar dois dos termos qua-
draticos acima mencionados. Assim sendo, as teorias de gravitacao de ordem superior

passaram, em geral, a serem denidas pela lagrangiana

2 2 2 
L= g R + R + R ,
2 2 2
onde 2 = 32G, sendo G a constante de Newton, e e sao parametros arbitrarios.

Com a constatacao de que a gravitacao nao era renormalizavel dentro do esquema

perturbativo padrao, as teorias de gravitacao de ordem mais alta, consideradas ate

entao como simples extensoes da relatividade geral, passaram a ocupar um lugar de

destaque na busca por uma teoria de gravitacao quantica. Neste sentido e digno de nota
o magistral trabalho de Stelle de 1977 [9], onde e demonstrado que as teorias de ordem

superior sao renormalizaveis juntamente com os seus acoplamentos com a materia.

Infelizmente estas teorias nao sao unitarias devido a presenca de um fantasma massivo

11
de spin 2. Em 1986, Antoniadis e Tomboulis [10] argumentaram que a presenca do
fantasma massivo de spin 2 no propagador nu era inconclusiva, ja que esta excitacao e

instavel. De acordo com eles, a posicao dos polos complexos no propagador vestido e
explicitamente dependente de gauge. Utilizando argumentos padroes de teoria quantica

de campos eles concluram que as teorias de gravitacao de ordem superior sao unitarias.
No ano seguinte Johnston [11] provou que as conjecturas de Antoniadis e Tomboulis

nao sao verdadeiras, uma vez que o par de polos complexos que aparece no propagador
re-somado sao independentes de gauge, implicando na nao-unitariedade das teorias de
gravitacao com derivadas de ordem mais alta.

Assim, acreditava-se ate 2009 que as teorias gravitacionais de ordem superior em

sua versao linearizada nao eram unitarias. Isto explica, em parte, o grande interesse

que a teoria de gravitacao massiva em 3D [12] provocou e, de certo modo continua

provocando, na comunidade cientca ja que a sua versao linearizada e unitaria.

Nos dois captulos que se seguem vamos discutir novas e interessantes propriedades
destas duas intrigantes teorias, a eletrodinamica de Lee e Wick e gravitacao massiva
em 3D, reservando o eplogo para a discussao sobre a possibilidade da unitariedade e

renormalizabilidade poderem coexistir pacicamente no contexto de uma mesma teoria

de ordem mais alta.

12
Referencias Bibliogracas

[1] B. Podolsky, Phys. Rev. 62 66 (1942).

[2] A. Pais and G. E. Uhlenbeck, Phys. Rev. 79 145 (1950).

[3] T. Lee and G. Wick, Nucl. Phys. B 9 209 (1969).

[4] T. Lee and G. Wick, Phys. Rev. D 2 1033 (1970).

[5] B. Grinstein, D. OConnell, and M. Wise, Phys. Rev. D 77 025012 (2008).

[6] H. Weyl, Space-Time- Matter Dover, New York (1952), Chap. IV.

[7] A. Eddington, The Mathematical Theory of Relativity, 2nd ed, Cambridge Univer-

sity Press (1924), Chap. IV.

[8] P. Havas,General Relativity and Gravitation 8 631 (1977).

[9] K. Stelle, Phys. Rev. D 16 953 (1977).

[10] I. Antoniadis and E. Tomboulis, Phys. Rev. D 33 2756 (1986).

[11] D.A. Johnston, Nuclear Physics B 297 721 (1988).

[12] E. Bergshoe, O. Hohm, and P. Townsend, Phys. Rev. Lett. 102 201301 (2009).

13
Captulo 1
Explorando a Eletrodinamica Finita de Lee
e Wick

Estudamos a eletrodinamica de Lee e Wick (LW), i.e., a teoria (invariante) de gauge

U (1) na qual um operador de dimensao 6 contendo derivadas de ordem mais alta e

adicionado a lagrangiana livre do setor U (1). Um limite quantico sobre a partcula

massiva de LW e estimado por meio do computo do momento magnetico anomalo

do eletron no contexto do ja citado modelo. Apresentamos tambem, um estudo

comparativo entre os potenciais de LW e Coulomb.

1.1 Motivacao

Lee e Wick (LW) propuseram uma teoria quantica para a eletrodinamica com o intuito

de entender a diferenca de massas dos mesons antes que a QCD se estabelecesse [1, 2].
Cerca de tres anos atras, as teorias do tipo LW foram redescobertas num outro contexto
quando Grinstein, OConnel e Wise introduziram teorias de gauge nao-Abelianas de LW

[3]. O modelo desenvolvido por eles, conhecido normalmente por Modelo Padrao LW
(MPLW), e naturalmente livre de divergencias quadraticas, fornecendo uma alternativa

para a solucao do problema da hierarquia. Mesmo sendo relativamente recente, ja existe

14
um extenso corpo de literatura relacionada ao MPLW [416], o que claramente mostra o
interesse consideravel que este modelo tem despertado. Estudos fenomenologicos sobre

o MPLW tambem podem ser encontrados em abundancia [1722]. Portanto, e bastante

justicavel explorarmos algo mais simples, como setor eletromagnetico da teoria de LW,

como uma especie de toy model para entendermos a dinamica mais complexa do MPLW.
De certo modo, este tipo de pesquisa e, mutatis mutandis, similar aquelas conduzidas

em teorias com dimensoes menores com o intuito de se ganhar um maior conhecimento

sobre questoes conceituais mais complicadas, as quais sao mais opacas no mundo fsico
quadridimensional.

Nosso objetivo neste captulo e justamente examinar duas caractersticas basicas

da QED de LW: (i) a ordem de magnitude da massa da partcula pesada de LW e,

(ii) o desvio do comportamento do potencial de LW nao-relativstico em comparacao

ao potencial usual de Coulomb. Espera-se que um melhor entendimento destes topicos

nos auxilie a, ao menos em princpio, melhorar nosso conhecimento do MPLW.

Este captulo esta organizado da seguinte maneira: alguns aspectos basicos do mo-
delo de LW serao analisados na Secao 1.2, incluindo uma discussao sobre o seu propa-

gador e renormalizabilidade. Na Secao 1.3 mostramos como e possvel obter o tensor de

energia-momento do modelo atraves de uma simples prescricao. Na Secao 1.4 faremos

uma analise do comportamento ondulatorio do modelo de LW enquanto, na Secao1.5,

calcularemos um limite classico para o cuto da teoria e, posteriormente, um limite

quantico para a massa da partcula pesada de LW sera estimado por meio do calculo do
momento magnetico anomalo do eletron, no contexto do citado modelo. O potencial
nao-relativstico para a interacao de dois fermions no contexto da QED de LW sera

calculado na Secao 1.6 por meio de uma poderosa abordagem baseada na fusao da
mecanica quantica com o limite nao-relativstico da teoria quantica de campos. Uma

15
analise comparativa entre os potenciais de LW e Coulomb e feita logo apos. Em seguida,
na Secao 1.7, fazemos uma breve discussao sobre a possibilidade da existencia de mo-

nopolos magneticos no citado modelo e, nalmente, uma discussao sobre os resultados

obtidos encerra a Secao 1.8.

Ao longo de todo o captulo, temos  = c = 1 e a assinatura da metrica dada por

(+1, -1, -1, -1).

1.2 Aspectos gerais do modelo de Lee e Wick

O modelo abeliano de LW e denido pela seguinte lagrangiana invariante de gauge

1 1
L = F F F 2F . (1.1)
4 4m2
Vamos mostrar que esta lagrangiana descreve (on-shell ) dois campos de spin-1 inde-
pendentes: um sem massa e outro massivo com normas positiva e negativa, respectiva-
mente. Para se alcancar este objetivo, e conveniente adotarmos uma nova formulacao

onde se introduz um campo auxiliar e o termo com derivadas de ordem mais alta
nao comparece. A teoria de campos com os campos vetoriais reais A e Z e com

lagrangiana [23]

1 1 m2 m2 m2
L = A 2Z + A Z

A A + A Z

Z Z , (1.2)
2 2 8 4 8

e equivalente a teoria de campos denida por (1.1). De fato, variando Z , obtemos

2 2
Z = A + 2
2A 2 A , (1.3)
m m

e o par de equacoes acopladas de segunda ordem oriundas de (1.2) e totalmente equiva-

lente as equacoes de quarta ordem geradas a partir de (1.1). O sistema (1.2) separa-se

16
agora em duas lagrangianas, cada uma associada a um campo diferente, quando reali-
zamos a mudanca de variaveis

A = B + C , (1.4)

Z = B C . (1.5)

Em termos de B , C , B B B e C C C , a lagrangiana toma

agora a seguinte forma

1 1 m2
L = B B + C C C C , (1.6)
4 4 2

que nada mais e que a diferenca entre a lagrangiana de Maxwell para B e a lagrangiana

de Proca para C .
Pode-se mostrar que a lagrangiana de LW pode tambem ser escrita como

1 1
L = F F + F F , (1.7)
4 2m2

ou, equivalentemente,

1 1
L = F F + F F . (1.8)
4 4m2

O conteudo de partculas da teoria pode tambem ser obtido diretamente de (1.1).

Para tal, vamos computar os resduos nos polos simples do propagador saturado (con-

tracao do propagador com correntes conservadas). Adicionando a (1.1) o termo de

1
xacao de gauge L = 2 ( A )2 , onde, como e usual, desempenha o papel de

parametro xador de gauge (e notando que devido a estrutura da teoria e a escolha

de um funcional xador de gauge linear, nao sao exigidos neste caso fantasmas de
Faddeev-Popov), podemos escrever a equacao (1.1) como

17
     
1 22 1 2
L = A
2+ 2 +

2 A . (1.9)
2 m m
No espaco dos momenta a equacao anterior toma a seguinte forma

1
L = A O A , (1.10)
2
onde
   
k4 2 1 k2
O
=
k k k

+ 2 . (1.11)
m2 m
k k
Expandindo o operador O em termos dos operadores = k kk2 e = k2
,

camos com
 
k4 2 k 2
O
= k
. (1.12)
m2
Com isso, a tarefa de encontrar o inverso do operador O(= A + B) se torna simples,
pois OO1 = I, onde I = + . Se O1 = C + D,

(A + B)(C + D) = AC + BD = + , (1.13)

1
ou seja, O1 = A
+ B1 .

Utilizando a relacao acima, onde O1 = D (k), obtemos prontamente o propagador

no espaco dos momenta

  2  
m2 k k k
D (k) = 2 2 2 1 + 1 . (1.14)
k (k m2 ) k m2

Contraindo (1.14) com correntes conservadas J (k), resulta

M J D J
J2 J2
= + ,
k 2 k 2 m2

18
o que nos permite concluir, levando em conta que J 2 < 0 [2426], que os sinais dos
resduos de M nos polos k 2 = 0 e k 2 = m2 sao, respectivamente,

ResM(k 2 = 0) > 0, ResM(k 2 = m2 ) < 0.

Vale a pena notar que o sinal errado da partcula pesada indica uma instabilidade

da teoria a nvel classico. Do ponto de vista quantico isto signica que a teoria e
nao-unitaria. Felizmente estas diculdades podem ser contornadas. Realmente, a ins-

tabilidade classica pode ser evitada impondo-se uma condicao de contorno no futuro de
modo a prevenir o crescimento exponencial de certos modos. No entanto, este proce-
dimento leva a violacao de causalidade no modelo [27]; esta acausalidade e suprimida,

por sua vez, em escalas inferiores aquelas associadas as partculas de LW. Por outro
lado, Lee e Wick argumentaram que apesar da presenca dos citados graus de liberdade

relacionados com uma norma nao positiva-denida no espaco de Hilbert, a teoria pode-
ria ser unitaria desde que as novas partculas obtivessem comprimentos de decaimento.

Nao existe uma prova geral da unitariedade em ordem de loop arbitrario para a eletro-
dinamica de LW; a despeito disto, nao existe ate agora nenhum exemplo conhecido de
violacao de unitariedade. Consequentemente, ate ordem em contrario, podemos consi-

derar a eletrodinamica de LW como nita e nao precisamos temer o fantasma massivo

de spin-1 [28].

Em resumo, podemos dizer que o trabalho de LW consiste na introducao de campos

de Pauli-Villars, com propagador com sinal errado, como graus de liberdade, o que

leva a amplitudes que tem melhor comportamento ultravioleta e tornam nita a QED
logaritmicamente divergente.
Uma questao importante concernente a eletrodinamica de LW e se ela e ou nao re-

19
normalizavel. Em princpio, teorias vetoriais massivas nao sao renormalizaveis devido
a seu propagador livre nao se aproximar de zero assintoticamente como no caso sem

massa. Existem, contudo, duas importantes excecoes a esta regra: (i) teorias de gauge
com quebra espontanea de simetria e, (ii) teorias com bosons vetoriais massivos neutros

acoplados a correntes conservadas. Apesar do modelo de LW nao poder ser conside-

rado um sistema puramente massivo, seu propagador livre aproxima-se de zero para

grandes valores de momentum, tornando assim a teoria renormalizavel por contagem

de potencias.

1.2.1 Analise de um toy model escalar

A m de ganharmos um pouco mais de conhecimento sobre como se comporta a part-

cula pesada de Lee e Wick, e interessante fazer a analise de um modelo mais simples.
Para tanto, utilizaremos um campo escalar e admitiremos que a lagrangiana para o

modelo associado a este campo e a seguinte:

1  2 2 1 2 2 1
Lescalar =

2
m g 3 , (1.15)
2M 2 3!
onde e um campo escalar. O correspondente propagador pode ser escrito como

i
D (p) = p4
. (1.16)
p2 M2
m2

Notemos que para M >> m, o propagador acima possui polos em p2 m2 e p2 M 2 .

O resduo no polo p2 = M 2 possui o sinal errado, o que, a princpio, sinaliza uma

instabilidade ou um estado de norma negativa que poderia comprometer a unitarie-

dade do modelo. Contudo, assim como zemos anteriormente, podemos escrever uma
lagrangiana equivalente a anterior utilizando dois campos reais e escalares e tais

20
que

1 1 1 1
Lequiv. = m2 2 2 + M 2 2 g 3 . (1.17)
2 2 2 3!
 
E imediato vericar que a equacao de campo para = M12 2 , quando substituda
em (1.17), recupera a lagrangiana em (1.15). Podemos interpretar melhor os campos

e se diagonalizarmos os termos cineticos. Para tanto, se denirmos = , a

equacao (1.17) pode ser reescrita como

1 1 1 1 1
L = + M 2 2 m2 ( )2 g( )3 . (1.18)
2 2 2 2 3!

Notemos agora que o propagador do campo (se desconsiderarmos a massa m, por

simplicidade) e dado por

i
D (p) = . (1.19)
p2 M2
Contudo, como observado por Lee e Wick, o sinal errado no propagador violaria a
unitariedade da matriz S se estes estados fossem estaveis. Assim, a unitariedade pode

ser preservada se decair em outras partculas mais leves. Vejamos como isso pode ser
obtido reescrevendo a lagrangiana em (1.18).

E possvel diagonalizar os termos massivos atraves de uma transformacao simpletica,

ja que a matriz de massa e simetrica [5]. Deste modo, podemos escrever a parte massiva

da lagrangiana (1.18) na forma LM = T M, com

Se V e uma matriz n n simetrica e de modulo positivo, entao D = S T V S e diagonal, onde S e
uma matriz simpletica.

21

  1 2 1 2
2 (2 + m ) 2
m
T = e M= . (1.20)
1 2 1 2 2
2
m 2
(2 m +M )

Utilizando a transformacao

cosh sinh 0
= , (1.21)
sinh cosh 0

camos com

12 m2 (cosh sinh )2 1 2
2 m (cosh sinh )
2
1 2 2
+ 12 M 2 sinh cosh  
  + 2 M sinh
0
LM = 0 0 .
0
1 2
2 m (cosh sinh )
2 12 m2 (cosh sinh )2
+ 12 M 2 sinh cosh + 12 M 2 cosh2

Se a matriz quadrada acima e diagonal, entao m2 (cosh sinh )2 = M 2 sinh cosh .

Isso implica em
2m2
tanh 2 = . (1.22)
2m2 M 2
Notamos que existe uma solucao para se garantirmos que M > 2m, de modo que
a partcula descrita por 0 decai em duas partculas do campo 0 . A lagrangiana em

termos dos campos redenidos se torna

1 1 1 1 1
L = 0 0 m0 2 20 0 0 + M02 20 (cosh sinh )3 g(0 0 )3 ,
2 2 2 2 3!
(1.23)
onde m0 e M0 sao as massas dos campos diagonalizados e a constante de acoplamento

pode ser redenida como g  = g

3!
(cosh sinh )3 .
A taxa de decaimento para o processo descrito na Figura (1.1) pode ser calculada

22
a partir de [29]

 1 n
d 3 pi
4 (4) 2
d = (2) (p pi ) |Mf i | , (1.24)
i
2Ep i=1
(2)3 (2Ei )

onde Ep e p sao, respectivamente, a energia e momentum da partcula inicial, e Ei e pi ,

a energia e momenta das partculas nais. Alem disso, a amplitude de Feynman, Mf i ,

em nvel de arvore, sera dada pelo vertice ig  .

Figura 1.1: Partcula pesada de Lee e Wick decaindo em outras duas mais leves

No referencial de repouso da partcula 0 , a equacao (1.24) ca

 
(2)4  2 d 3 p1 d 3 p2
= | ig | 3 3
(M E1 E2 ) (3) (p1 + p2 )
2Ep (2) (2E1 ) (2) (2E2 )
2  3
1 g d p1
= 2
(M E1 E2 ). (1.25)
(2) 2Ep 4E1 E2

Como as partculas resultantes possuem a mesma massa e momenta opostos, passando

23
para coordenadas esfericas, temos
 2
1 g 2 p dp
= (M 2E)
2 8Ep E2
g 2 p 
=  , (1.26)
16M E E= M2

onde zemos uso de EdE = pdp e E1 = E2 = E, com |p1 | = p. Agora, lembrando que

p = E 2 m2 , obtemos prontamente

g 2 4m2
= 1 . (1.27)
32M M2

Como foi demonstrado por Cutkosky [30], uma modicacao apropriada do tra-
tamento perturbativo que nos permite escrever o propagador total para a partcula

descrita por , utilizando a auto-energia em um loop, proximo a p2 = M 2 , nos da

i i 2 i
D = + (i(p )) + ...
p2 M 2 p2 M 2 p2 M 2
i
= 2 . (1.28)
p M 2 + (p2 )

Notemos que, ao contrario do que acontece com escalares normais, o sinal de no

denominador e positivo, ou seja,

i
D = , (1.29)
p2 M 2 iM

com dado por (1.27). Este sinal contrario faz com que a ressonancia pareca decair
antes que ela seja propriamente produzida, i.e., possui uma probabilidade dt de decair

no intervalo negativo de tempo dt, o que faz com que a teoria viole a microcausali-
dade. Pode-se pensar que tal comportamento acausal leve a paradoxos, mas se a matriz

S for unitaria, nao havera problemas: estados iniciais com partculas estaveis sao pre-

24
parados em e, como resultado, estados nais sao encontrados em +. A soma
das probabilidades para os estados nais e igual a um e, portanto, nao ha paradoxos.

Macroscopicamente, a causalidade e preservada e a teoria se torna unitaria pois

a parte imaginaria da amplitude de espalhamento ainda e uma quantidade positiva.

Proximo a p2 = M 2 o espalhamento e dominado pelo polo em e parte imaginaria da

amplitude de espalhamento M e, para Im p2 > 0, dada por

M
ImM = g 2 . (1.30)
(p2 M 2 )2 + M 2 2
Assim, o sinal errado do propagador e compensado pelo sinal contrario do comprimento
de decaimento [3].

1.3 O tensor energia-momento de Lee e Wick

Vamos construir o tensor momento-energia simetrico do modelo de Lee e Wick utili-

zando um algoritmo bastante simples baseado nas equacoes de campo da teoria [31].

Esta prescricao consiste basicamente em multiplicar a equacao de campo em ques-

tao por uma derivada apropriada deste campo de modo que a expressao resultante

contenha somente um ndice de espaco-tempo livre e, entao, reescreve-la como uma

quadridivergencia.

Na ausencia de fontes, a equacao de campo para a eletrodinamica de Lee e Wick se

torna

 2
1 + 2 F = 0. (1.31)
m
Multiplicando ambos os lados da equacao por F , resulta

25
 1
F F + F F = 0. (1.32)
m2
Podemos reescrever o primeiro termo da equacao anterior como se segue

1
F F = (F F ) ( F )F = (F F ) ( F F )F
2
1 1
= (F F ) ( F + F )F = (F F ) + ( F )F
2 2
1 1
= (F F ) + (F F ) = (F F + F F ),
4 4

onde zemos uso das identidades de Bianchi ( F + F + F = 0). Ja o segundo

termo da equacao (1.32) pode ser reescrito como

F 2 F = (F 2F ) F 2F
1
= (F 2F ) ( F + F )2F
2
1
= (F 2F ) + F 2F

2  
1 1
= (F 2F ) +

F 2F
F 2 F .
2 2

Porem,

1 1
F 2 F = F 2( F + F ) = F 2 F
2 2
= (F 2F ) F 2F = (F 2F ) + F ( F + F )
 
1
= ( F F ) + F F .
2

26
Com isso, o segundo termo em (1.32) se torna

1 1 1
F 2 F = 2 F 2F + (F 2F + F F ) + F 2F + F F .

m2 m 2

Assim, o tensor de energia-momento pode ser expresso, nalmente, como

1 1  
T = F F + F F + 2
F 2F + F F
4 2m
1  
2 F 2F + F 2F + F F . (1.33)
m

Podemos escrever o momento total do campo eletromagnetico P = (E, p), onde E

e p sao, respectivamente, a energia e o tri-momento totais do campo, como


P =
d3 xT 0 . (1.34)

No caso eletrostatico, as equacoes (1.33) e (1.34) nos dao a seguinte expressao para
a energia



1 1  2 2
 3
E= E 2 ( E) + 2E E d x. (1.35)
2 m
Da identidade vetorial ( E) = ( E) 2 E, encontramos 2 E = ( E),
uma vez que E = 0. Alem disso, [( E)E] = E ( E) + ( E)2 e, entao,
1
se supormos que E E se anule no innito mais rapidamente que r2
, encontramos

1 1
E= [E2 + ( E)2 ]d3 x. (1.36)
2 m2

Conforme a equacao (1.49), o potencial eletrostatico para uma carga puntual pode
ser escrito como V (r) = e
4r
[1 emr ], o qual tende ao valor nito em
4
quando r se

aproxima de zero. Lembrando agora que E = V , obtemos a expressao para o

Em nossas convencoes, F 0i = E i e F ij = ijk Bk .

27
campo eletrostatico devido a carga puntual em questao
  
e 1 mr 1
E(r) = e m r. (1.37)
4r r r

Substituindo a equacao (1.37) em (1.36), apos o calculo da integral, encontramos que a

me2
energia para o campo de uma carga puntual e dada por 2
, o que nos mostra que, ao

contrario da eletrodinamica de Maxwell, a eletrodinamica de LW exibe um valor nito

para a energia em todo o espaco.

1.4 Propagacao ondulatoria na eletrodinamica de

Lee e Wick

Na eletrodinamica de LW, o potencial A obedece as equacoes

 2
1 + 2 (2A A ) = 0, (1.38)
m

que no gauge de Lorentz reduzem-se a

 2
1 + 2 2A = 0, (1.39)
m

A m de encontrarmos as relacoes de dispersao, passamos para o espaco dos mo-

menta via o ansatz A (x) = A (k)eikx . Levando este resultado em (1.39), obtemos

prontamente que
 
k2
1 2 k 2 = 0. (1.40)
m
Como k = (, k), podemos reescrever esta equacao como

2 2 k2
2
( k ) 1 = 0. (1.41)
m2

28
A equacao (1.41) admite duas classes de solucao. A primeira famlia e a das solucoes
usuais da eletrodinamica para as quais vale a relacao de dispersao = |k|. A outra

famlia satisfaz a relacao k2 = 2 m2 , que pode ser colocada na forma mais familiar

|k| = 2 P2 , (1.42)

com P = m. Esta equacao descreve ondas eletromagneticas propagando-se atraves

1 1
de um plasma sem colisoes com um comprimento de atenuacao P
= m
muito

pequeno [3234].
Examinando a equacao (1.42) vemos que |k| e real somente se > P . Se < P ,
|k| e um numero imaginario puro, tais modos nao se propagam e sao chamados de

modos evanescentes (vide Fig. 1.2).

Figura 1.2: Relacoes de dispersao para a eletrodinamica de LW. A linha contnua denota o
modo propagante e a linha tracejada representa a parte imaginaria do modo evanescente.

Vamos estimar na proxima secao um valor para m, o que nos permitira constatar que << 1.

29
 Podemos agora determinar as velocidades de grupo e fase para os modos dispersivos
propagantes. A velocidade de grupo e dada por

|k|
vg = = , (1.43)
|k|

ao passo que a velocidade de fase e determinada pela relacao

vp = . (1.44)
|k|

Analisando (1.43) e (1.44) conclumos que vg < 1 e vp > 1, o que concorda com a ideia

que estamos perseguindo de ser esta propagacao eletromagnetica analoga a propagacao

de ondas eletromagneticas num vacuo tipo-plasma.

Das consideracoes anteriores conclui-se que m e um cuto para o modelo de LW.

Seria fundamental, portanto, estimar um limite para este parametro. Examinemos,
com este intuito, o potencial dado pela equacao (1.49). O fato deste potencial ser

nito quando r 0 e uma indicacao de que a auto-energia e a massa eletromagnetica

de uma partcula carregada puntiforme na teoria de LW sao nitas. Na realidade,

como argumentaremos abaixo, no ambito desta teoria nao so a massa eletromagnetica

e nita, como o celebre problema dos 4/3 da eletrodinamica de Maxwell encontra uma

resolucao natural.
A equacao de movimento de Newton relativa a uma carga estendida no contexto do
modelo de Abraham-Lorentz e dada por [33]

4 2
me v e2 v = Fext , (1.45)
3 3

onde v, me e e sao, respectivamente, a velocidade, a massa eletromagnetica e a carga da

partcula, e Fext e a forca externa aplicada ao eletron. Supomos aqui que, no referencial

30
de repouso instantaneo da partcula, a distribuicao de carga e rgida e esfericamente
simetrica. Alem disso, supomos tambem a ausencia de massa mecanica. Evidentemente

a massa eletromagnetica entra com coeciente errado em (1.45). Por esta razao o fator
4/3 tem sido fonte de acirrados debates. Uma possvel sada para esta diculdade,

entre outras, e calcular a auto-forca que atua numa partcula puntiforme carregada
no contexto da eletrodinamica de LW em pequenas distancias. Neste caso a massa

eletromagnetica entra na equacao de movimento em uma forma consistente com a

relatividade especial; e mais, a equacao de movimento correta nao exibe nem runaway
1
solutions nem comportamento acausal, quando o cuto l( m
) e maior que o raio

classico do eletron [35, 36]. Consequentemente l < 1.4 1015 m, o que nos permite
concluir que m > 141M eV .

1.5 Limites sobre a massa da partcula pesada de

Lee e Wick

Nesta secao apresentaremos dois metodos para se determinar um limite sobre a massa

da partcula de Lee e Wick. O primeiro segue uma abordagem classica e se baseia em

um experimento bem conhecido, enquanto o segundo faz uso das ultimas medidas mais

recentes relativas ao momento magnetico anomalo do eletron.

1.5.1 Um limite classico sobre a massa da partcula pesada de

Lee e Wick

Vale a pena notar que o parametro l introduz uma escala natural de comprimento para

a eletrostatica de LW. Na presenca de uma fonte, Lfonte = A J , as equacoes de

campo concernentes ao modelo de LW sao

31
 (1 + l2 2) F = j , (1.46)

F = 0, (1.47)

onde F = 12 F (0123 = +1). No limite de cargas estaticas, as equacoes anteri-

ores se reduzem a

(1 l2 2 )2 = , (1.48)

com E = .
No experimento de Plimpton e Lawton [37], por sua vez, foi medida uma diferenca
de potencial V entre uma esfera condutora de raio R2 carregada ate um potencial

V0 e uma segunda esfera concentrica descarregada de raio R1 contida na primeira. Da

equacao (1.48) obtemos o potencial eletrostatico


(r) = G(r r )(r )d3 r , (1.49)

|rr |
1 1e l
onde G(r r ) = 4 |rr |
, e (r ) e a densidade de carga. Em coordenadas esfericas,

a equacao anterior se torna

 +1 r 2 +r 2 2rr  cos
Q 1e l
(r) = d(cos ). (1.50)
8 1 r2 + r2 2rr cos

Fazendo u = r2 + r2 2rr cos , obtemos
 r  +r
Q
(1 e l )du,
u
(r) = (1.51)
8rr r  r

e, realizando a integracao, encontramos que o potencial para uma distribuicao esferica

32
de cargas e
Q   r +r
l

r r

(r) = 2r + l e e l . (1.52)
8rr
Agora, de acordo com a Figura (1.3), supondo que a casca esferica externa seja

mantida a um potencial V0 e fazendo r = r = R2 na equacao (1.52), obtemos pronta-

mente

Figura 1.3: Geometria do experimento de Plimpton e Lawton.

8R22 V0
Q=  2R , (1.53)
2
2R2 + l e l 1

o que nos da   R +r 
2 R2 r
V0 R2 2r + l e l e l
(r) =   2R  . (1.54)
l2
2R2 + l e 1 r

33
 Calculando a diferenca de potencial entre as esferas, encontramos
  R +R R R

2l 1 2l 1
V (R2 ) (R1 ) R 2 2R 1 + l e e
=1   2R  . (1.55)
V (R2 ) 2
2R2 + l e l 1 R1

R1 R2
Agora, levando em conta que tanto l
quanto l
sao >> 1, obtemos imediatamente

V l
. (1.56)
V 2R2

No experimento originalmente conduzido por Plimpton e Lawton, um potencial de

3000 volts harmonicamente alternado foi aplicado a esfera exterior. Testes foram feitos

a m de detectar uma mudanca no potencial da esfera interna em relacao a externa. Os

resultados mostraram que nenhuma mudanca no galvanometro foi detectada durante o

experimento, com uma sensibilidade de 1V . Os raios das duas esferas eram de 0.76m e
0.61m, respectivamente. Substituindo os parametros experimentais na equacao (1.56),

encontramos um limite sobre o cuto l < 5.1 1010 m [38].

Muitas melhorias ao experimento de Plimpton e Lawton foram realizadas ao longo
dos anos 60 80 [39]. Em todas as novas experiencias, contudo, a distancia entre os

pontos medidos eram, grosso modo, maiores do que aquela utilizada no experimento
original. Agora, ja que o cuto l introduz uma escala de comprimento natural para a

eletrostatica, como ja comentamos, obviamente para uma dada precisao das medidas,
quanto maior for a distancia entre os pontos medidos, pior sera o limite obtido para o

cuto. Portanto, chegamos a conclusao de que o experimento original de Plimpton e

Lawton e justamente o que fornece o melhor limite superior sobre l.
Em termos da massa, esse cuto nos da um limite inferior igual a m > 0.387

103 M eV .

34
1.5.2 Um limite quantico sobre a massa da partcula pesada
de Lee e Wick

Levando em conta que a QED prediz o momento magnetico anomalo do eletron corre-

tamente ate a decima casa decimal, um limite quantico para a massa m da partcula
pesada de LW pode ser encontrado calculando-se o momento magnetico anomalo do

eletron no contexto da eletrodinamica de LW e comparando o resultado obtido com

aquele da QED [40]. Para tanto, relembramos que o momento magnetico anomalo vem
da correcao de vertice devido ao espalhamento de um eletron por um campo externo,

como mostra a Figura 1.4.

p

p k

q=p p
k

pk
p

Figura 1.4: Correcao de vertice para um eletron espalhado por um campo externo.

Para um eletron espalhado por um campo magnetico estatico externo e, no limite

q 0, a razao giromagnetica e dada por [41]

g = 2[1 + 2F2 (0)].

O fator de forma do eletron, F2 (0), corresponde a um deslocamento no fator-g, geral-

mente escrito como F2 (0) g2
2
, e garante o momento magnetico anomalo do eletron.

Utilizando (1.14) no calculo do diagrama apresentado na Fig. 1.4, pode ser mostrado

35
que

 
1 12 (2 + 3 )
F2 (0) = d1 d2 d3 (1 i ) , (1.57)
0 2 + 3 (2 + 3 )2 + 1

M2
onde m2
, sendo M a massa do eletron. Chamamos atencao ao fato de que o termo
2
  2 
k4m(k2km
k
2) 1+ m k
2 1 que aparece na Eq. (1.14) nao contribui para o fator de
forma F2 (0) porque o propagador sempre aparece acoplado a correntes conservadas.

Integrando a expressao acima, primeiro em relacao a 3 e, em seguida, em relacao

a 2 , temos

 1  11 
1 1 (1 1 )
F2 (0) = d1 d2
0 0 1 1 (1 1 )2 + 1
 1
12
= d1 .
0 1 + (1 1 )2

Agora,


x2 x 1 2  2 
dx = ln x + x(1 2) + 
x2 + x(1 2) + 2
 
1 + 22 4  A B 
+ 2 ln ,
2 1 4  A + B 

onde A 2x + 1 2, e B 1 4; entao


1 1 2 1 + 22 4 1 + 1 4
F2 (0) = + ln + 2 ln . (1.58)
22 2 1 4 1 1 4
Lembrando que 1, o que implica

36
 1 + 22 4 1 2 + 54
,
22 1 4 22


1 + 1 4 2 203 354
ln 2 + 3 + + + ln , (1.59)
1 1 4 3 2
chegamos a conclusao de que

  2     4   
6
2 M 25 M M M
F2 (0) 1 2 + ln +O . (1.60)
2 3 m 12 m m m

O primeiro termo da equacao anterior e igual aquele calculado por Schwinger em 1948
[42]. Desde entao, F2 (0) tem sido calculado ate a ordem de 8 para a QED. O segundo

termo da equacao Eq. (1.60) e a correcao mais importante relacionada ao parametro m

da eletrodinamica de LW. Calculos recentes concernentes a F2 (0) no contexto da QED

dao, para o eletron [43],

F2 (0) = 1 159 652 181.13(.84)) 1012 ,

onde a incerteza vem, em sua maior parte, dos melhores valores para a constante de

estrutura na nao advindos da QED. O atual valor experimental para o momento

magnetico anomalo e, por sua vez [44],

F2 (0) = 1 159 652 180.73(0.28) 1012 .

Comparando o valor teorico predito pela QED com o valor experimental, vemos que
estes resultados concordam em 1 parte em 1011 . Consequentemente,

37
  2
2 M
< 1011 .
3 m
Assim, um limite inferior para a partcula pesada hipotetizada por Lee e Wick e m

132 GeV.
O limite quantico que estimamos e extremamente conavel ja que foi obtido utili-

zando um dos grandes triunfos da QED, a saber, a surpreendente concordancia entre

teoria e experimento ate dez casas decimais para o momento magnetico anomalo do
eletron. Alem disso, este limite e baseado em verdadeiros efeitos quanticos (de loop).

1.6 O potencial nao-relativstico de LW

Vamos agora calcular o potencial nao-relativstico para o espalhamento de dois fermions

no contexto da QED de LW por meio de um metodo baseado na combinacao da meca-

nica quantica com o limite nao-relativstico da teoria quantica de campos. Utilizando
este resultado juntamente com o limite quantico que acabamos de estimar, fazemos em

seguida um estudo comparativo entre os potenciais de LW e Coulomb.

1.6.1 Combinando a mecanica quantica com o limite nao-relativstico

da teoria quantica de campos

Consideremos o espalhamento de dois fermions de massas identicas no contexto do

modelo de Lee e Wick; especicamente, um eletron (e ) e um positron (e+ ), como

mostrado na Fig. 1.5. Agora, de acordo com a aproximacao de Born na teoria de

Schrodinger, a secao de choque diferencial para o espalhamento de duas partculas de

 d  M 
3 2
mesma massa, M , e d =  e Qr
U (r)d r . (A energia potencial U (r) causa
C.M. 4

a troca de momentum Q = p p, onde p e p sao os momenta inicial e nal de uma das

38
partculas). Por outro lado, no limite nao-relativstico, a secao de choque diferencial
 d   M 2
para a interacao entre estes fermions pode ser expressa como d =  N.R.  ,
C.M. 16M

onde MN.R. e o limite nao-relativstico da amplitude de Feynman para o processo em

questao. Como resultado, em termos do momentum trocado k = Q, o potencial

nao-relativstico assume a forma


1 1
U (r) = d3 kMN.R. eikr . (1.61)
4M (2)3
2

Esta formula nos permite construir um potencial tridimensional efetivo (para ser utili-
zado em conexao com a equacao de Schrodinger bem como na mecanica classica) uma

vez que conhecemos o limite nao-relativstico do elemento de matriz covariante.

Figura 1.5: Contribuicao de menor ordem para o espalhamento elastico e e+ no contexto

da eletrodinamica de LW.

A lagrangiana do sistema fermionico em questao pode ser escrita como

1 m2 1
L = F F + F F ( A )2 + (i/ M ) eA . (1.62)
4 2 2

Neste caso, o vertice fundamental descrito pelo termo de interacao na lagrangiana sera

39
dado por
= e . (1.63)

Considerando o caso de partculas com o mesmo spin nos estados iniciais e nais,

podemos escrever (x) como uma expansao em ondas planas e espinores como (x) =
u(p)eipx , tais que ur (p) us (p) = 2p rs . Partindo da lagrangiana para a interacao

das partculas mencionadas, i.e., Lint = e A , obtemos prontamente as regras de

Feynman para o vertice elementar (veja a Fig. 1.6)

Figura 1.6: Vertice relevante para a QED de LW.

M = e2 u(p ) u(p) v(q  ) v(q) (1.64)

 
k2
e m2 2 1 1 m2
= 2 2 u(p 
) u(p) k k v(q  ) v(q). (1.65)
k (k m2 ) k2

Notemos que segundo termo do propagador nao contribuira para a amplitude, pois

encontra-se entre correntes conservadas e e identicamente nulo. Atentemos ao fato

de que, nas nossas convencoes, a carga do eletron e igual a e. Assim, a amplitude

invariante para o processo mostrado na Figura 1.5 e

40
 e2 m2
M= u(p ) u(p)v(q) v(q  ). (1.66)
k 2 (k 2 m2 )

Como estamos interessados no limite nao-relativstico da amplitude, no qual M 2 >> p2 ,

temos

p2 = E 2 p2 = M 2
 1/2
2 2 1/2 p2
E = (M + p ) =M 1+ 2 , mas M 2 >> p2 , entao
M
 
1 p2
E M 1+ M.
2 M2

Assim, nesta aproximacao,

p q = p0 q 0 p q

= M 2 p q M 2.

Alem disso, pela conservacao de momento, temos as seguintes relacoes

p k = p (p k)2 = M 2

q + k = q  (q + k)2 = M 2

(p k)2 = (q + k)2 ,

41
que se reduz a

q k = p k

(q + p) k = 0

(q0 + p0 )k0 = (q + p) k
(q + p) k
k0 = .
(q 0 + p0 )

Ainda considerando o limite nao-relativstico, (q 0 + p0 ) = (Eq0 + Ep0 ) ET  (q + p)

(onde ET e a energia total do sistema) e, portanto

(q + p) k
k0 = 0. (1.67)
ET

Assim, o termo em 1/k 2 na amplitude pode ser escrito como

 
1 1 1 1 k2
= =   2 1 + 02 . (1.68)
k2 k0 k2 k2 1
2
k0 k k
k2

Contudo, o termo mais importante na aproximacao nao-relativstica pode ser obtido

utilizando-se solucoes de momentum nulo para os espinores, ou seja,

ur (p ) us (p) ur (0) us (0),

N.R.
(1.69)

vr (p ) vs (p) vr (0) vs (0),

N.R.
(1.70)

com a normalizacao ur (p)us (p) = vr (p)vs (p) = 2Ep rs [45]. Consequentemente,

ur (0) 0 us (0) = vr (0) 0 vs (0) = 2M rs , (1.71)

ur (0) i us (0) = vr (0) i vs (0) = 0. (1.72)

42
Como resultado,

4e2 m2 M 2
MN.R. = . (1.73)
k2 (k2 + m2 )

Inserindo (1.73) em (1.61) e calculando a integral resultante, temos


1 1
U (r) = MN.R. eikr d3 k
4M 2 (2)3 V

m2 e 2 eikr
= d3 k. (1.74)
(2)3 V k2 (k2 + m2 )

Reescrevendo a integral acima em coordenadas esfericas, com r = |r| e k = |k|, camos

com
 +  1  2
m2 e2 2 eikr cos
U (r) = k dk d(cos ) d. (1.75)
(2)3 0 1 k 2 (k 2 + m2 ) 0

As integrais em d e d(cos ) sao triviais, o que nos da

 +
m2 e2 sin kr dk
U (r) = ,
(2)2 kr k 2 + m2

que pode ser reescrita como

 +  +
e2 sin x sin x
U (r) = dx dx (onde zemos x = kr).
(2)2 r x x2 + m 2 r 2

Resolvendo cada um dos termos como a parte real de uma integral no plano complexo,
obtemos

e2
U (r) = 2
( emr ) e,
(2) r
e2
U (r) = (1 emr ). (1.76)
4r

43
1.6.2 Analise comparativa entre os potenciais de LW e Cou-
lomb

Podemos agora discutir as propriedades do potencial nao-relativstico U (r). Por con-

veniencia, vamos reescreve-lo como

e2  r
U (r) = 1 e l , (1.77)
4r

1
onde l m
, como vimos anteriormente.

Diferentemente do potencial de Coulomb, o qual e singular na origem, U e nito

2
(U (0) = 4l
e
), o que claramente nos mostra que em pequenas distancias l se
comporta como um cuto efetivo.

Quando r
l
 1, U se reduz ao potencial de Coulomb, como era de se esperar.

Apenas para pequenas distancias ( rl 1), U difere signicativamente do poten-

cial coulombiano (veja a Fig. 1.7). Observando a Figura 1.7, notamos que esta
distancia e 1017 m, uma ordem de magnitude maior do que o valor de l.

1.7 Monopolos magneticos

Uma questao importante que diz respeito a QED nita de LW e se a coexistencia de

uma partcula massiva e uma carga magnetica e possvel neste contexto. Para responder

esta questao, introduzimos a corrente magnetica k = (, k) no lado direito da equacao

(1.47). E facil mostrar que o sistema de equacoes de ordem superior modicadas
resultante, a saber,

44
 e2
Figura 1.7: O potencial U (em unidades de como funcao da distancia r. A linha
4 ),
pontilhada representa o potencial de Coulomb nas mesmas unidades.

(1 + l2 2) F = J , (1.78)

F = k , (1.79)

descreve a existencia de uma carga magnetica. De fato, supondo a ausencia de campos

eletricos, cargas e correntes (bem como as correntes magneticas), camos essencial-

mente com duas equacoes para o campo magnetico, as quais possuem a familiar solucao
g
de monopolos de Dirac B = 4r 2
r, onde g e a carga magnetica. Utilizando os metodos
qg
usuais, recuperamos prontamente a condicao de quantizacao de Dirac 4
= n2 , onde q

45
e a carga eletrica e n e um inteiro.
Conseguimos, deste modo, encontrar um sistema consistente de equacoes do tipo

Maxwell + bosons vetoriais massivos + cargas magneticas. Chamamos atencao ao

fato de que o monopolo de Dirac e bosons vetoriais massivos nao podem coexistir no

contexto da eletrodinamica massiva de Proca, pois esta, ao contrario da QED de LW,

nao e invariante de gauge [46]. A propria existencia de um monopolo de Dirac esta,

sem duvidas, ligada a existencia de uma invariancia de gauge da teoria correspondente.

Curiosamente, o sistema formado pelas equacoes (1.78) e (1.79) nao e simetrico
sob a transformacao de dualidade F F , F F , aumentada por j

k , k j . Este fato levanta uma questao interessante: Seria possvel acomodar

simultaneamente cargas magneticas e transformacoes de dualidade no contexto de um

modelo eletromagnetico com derivadas mais altas? Uma boa tentativa nesta direcao
pode ser, por exemplo, o modelo denido pelas equacoes

(1 + l2 2) F = J , (1.80)

(1 + l2 2) F = k , (1.81)

ja que ele e simetrico sob transformacoes de dualidade. Vale a pena notar que (1 +
l2 2) F e identicamente nulo na ausencia de correntes magneticas. Vejamos se este

modelo admite solucoes do tipo monopolo. Para uma carga magnetica estatica de inten-
 r/l 
1e er/l
g
sidade g xa na origem, a solucao das equacoes anteriores e B = 4 r 2 lr
r,

a qual, para largas distancias, se reduz ao resultado de Dirac, como era de se esperar.
Nosso argumento, no entanto, e vericar se esta solucao descreve ou nao um monopolo

magnetico em pequenas distancias.

46
 Para tanto, calculamos o uxo do campo magnetico radial atraves de uma superfcie
esferica S de raio r com um monopolo de intensidade g em seu centro. Realizando os
   r
calculos, encontramos S B dS = g 1 1 + rl e l , o que implica que para r/l 1,

S
B dS 0. Agora, levando em conta que se B = A, S
B dS e nula, chegamos a

conclusao de que A pode existir em toda regiao sob consideracao. Portanto, esta e uma
solucao do tipo de Dirac; nao obstante, a carga magnetica correspondente nao obedece
g 1
a condicao de quantizacao de Dirac. De fato, para r/l 1, B 4l2 r
r, o que implica
que o campo magnetico decresce como 1/r em vez de 1/r3 . Este comportamento bizarro
do campo magnetico certamente nos impede de recuperar a condicao de quantizacao

de Dirac.
Um modo heurstico de ver isto e considerar o movimento de uma partcula de

massa m e carga q no campo de um monopolo magnetico. Da equacao de movimento

da partcula, mr = q r B, encontramos que a variacao do seu momento angular e
qgr 2 d
d
dt
(r mr) = 4l2 dt
(r), um resultado que nos impede de denir um momento angular
total conservado, como no caso dos monopolos de Dirac. Agora, se as distancias nao
sao nem tao grandes nem tao pequenas, o potencial vetor nao pode existir em qualquer

lugar do domnio encerrado por S porque F satisfaz a equacao (1.81) inves de (1.79).
Infelizmente nao podemos superar esta diculdade por meio da introducao do conceito
g er/l
de uma corda assim como fez Dirac, ja que neste caso B (= 4 rl2
) nao se anula
em qualquer lugar do domnio em questao.

Esta analise nos leva a conjectura de que monopolos do tipo de Dirac e trans-
formacoes de dualidade nao podem ser acomodadas no contexto do mesmo modelo
eletromagnetico com derivadas de ordem mais alta [47].

47
1.8 Discussao

Como sabemos, cada campo no Modelo Padrao possui um parceiro de LW com uma

massa de LW associada. Estas massas sao os unicos parametros novos no MPLW

mnimo. Contudo, como observado por Alvarez, Schat, Da Rold e Szynkman, o MPLW
nao fornece nenhuma informacao sobre a origem destas massas, as quais, para resolver

o problema da hierarquia, nao devem ser mais pesadas do que alguns TeV. De fato,
numa analise dos dados eletrofacos (assumindo que todas as massas de LW sao da

mesma ordem) encontramos que a escala de LW deve ser da ordem de 5 TeV [17]. Por
outro lado, o limite para a partcula pesada de LW que estimamos esta na escala de

GeV, o que implica que ela poderia, em princpio, ser determinada no LEP. Contudo,
temos que admitir que o cenario analisado, o qual permite esta estimativa na escala de
GeV, e tal que o termo de derivadas mais altas e adicionado ao foton depois que uma

quebra espontanea de simetria ja aconteceu, de modo que derivadas de ordem mais

alta de Z0 nao estao presentes e a matriz de massa no setor de gauge de bosons neutros

ja foi diagonalizada. Nao obstante, a origem usual das interacoes eletromagneticas e

eletrofracas, como mostra a Teoria Eletrofaca com SU (2) U (1), nos diz que o termo

com dimensao 6, uma vez introduzido na lagrangiana antes da quebra espontanea

da simetria SU (2) U (1), deve ser governado por um parametro com dimensao de
massa2 , o qual deve estar conectado a um parametro originado em uma fsica BSM .

Tal parametro pode ser resultado de uma dimensao extra maior e, assim, acomodaria
de modo mais natural o resultado de um foton massivo na escala de TeV, levando

em conta modelos BSM. Neste contexto, seria oportuno reavaliar a fsica do Z0 no

modelo Eletrofaco de LW. Estamos aqui interessados nos possveis efeitos advindos da

Beyond Standard-Model Alem do Modelo Padrao.

48
associacao dos bosons Z0 e Z0 com o limite sobre o boson de gauge pesado de LW.
Para concluir, alem de reconsiderarmos a ja mencionada fsica envolvendo Z0 no

contexto do modelo de LW, seria razoavel reavaliar tambem a questao da producao de

gravitons massivos no LHC no contexto de uma versao gravitacional do modelo de LW.

Possveis sinais gravitacionais na escala de TeV no LHC, associados a efeitos de dimen-

soes extras maiores ou torcao dinamica num cenario BSM, poderiam ser reproduzidos

numa contrapartida gravitacional do modelo de LW. Isto nos permitiria relacionar di-
retamente o parametro de massa2 tpico do modelo de LW a alguma dimensao extra
maior e, neste contexto, tambem poderia se tornar mais claro a presenca de escalas de

massa da ordem de GeV e TeV para duas geracoes de gravitons pesados do tipo LW.

49
Referencias Bibliogracas

[1] T. Lee and G. Wick, Nucl. Phys. B 9 209 (1969).

[2] T. Lee and G. Wick, Phys. Rev. D 2 1033 (1970).

[3] B. Grinstein, D. OConnell, and M. Wise, Phys. Rev. D 77 025012 (2008).

[4] C. Carone, Phys. Lett. B 677 306 (2009).

[5] C. Carone and R. Lebed, J. Higher Energy Phys. 01 043 (2009).

[6] C. Carone and R. Lebed, Phys. Lett. B 668 221 (2008).

[7] E. Gabrielli, Phys. Rev. D 77 055020 (2008).

[8] F. Wu and M. Zhong, Phys. Lett. B 659 694 (2008).

[9] F. Wu and M. Zhong, Phys. Rev. D 78 085010 (2008).

[10] A. Rodigast and T. Schuster, Phys. Rev. D 79 125017 (2009)

[11] A. van Tonder, arXiv:0810.1928v1.

[12] A. Shalaby, arXiv:0812.3419v2.

50
[13] B. Grinstein, D. OConnell, and M. Wise, Phys. Rev. D 77 065010 (2008).

[14] I. Cho and O-Kab Kwon, arXiv: 1003.2716v1 .

[15] B. Fornal, B. Grinstein, and M. Wise, Phys. Lett. B 674 330 (2009).

[16] Yi-Fu Cai, T. Qiu, R. Brandenberger, and X. Zhang, Phys. Rev. D 80 023511

(2009).

[17] E. Alvarez, C. Schat, L. Da Rold, and A. Szynkman, J. Higher Energy Phys. 04

026 (2008).

[18] T. Underwood and R. Zwicky, Phys. Rev. D 79 035016 (2009).

[19] T. Rizzo, arXiv:0712.1791v2 .

[20] C. Carone and R. Primulando, Phys. Rev. D 80 055020 (2009).

[21] J. Espinosa, B. Grinstein, D. OConnell, and M. Wise, arXiv:0705.1188v1 .

[22] T. Rizzo, arXiv:0704.3458v3 .

[23] A. Accioly, Investigations in the Electromagnetic Sector of the Lee-Wick Standard

Model (submetido ao J. High Energy Phys.).

[24] A. Accioly, Phys. Rev. D 67 127502 (2003).

[25] A. Accioly, Nucl. Phys. B (Proc. Suppl.) 127 100 (2004).

[26] A. Accioly and M. Dias, Int. J. Theor. Phys. 44 1123 (2005).

[27] S. Coleman, in Proceedings of Erice 1969, Ettore Majorana School on Subnuclear

Phenomena (Academic Press, New York, 1970), pp. 282-327.

51
[28] B. Grinstein, D. OConnell, and M. Wise, Phys. Rev. D 79 105019 (2009)

[29] M. Maggiore, A modern introduction to quantum eld theory, Oxford University

Press (2005).

[30] R. E. Cutkosky, P. V. Landsho, D. I. Olive and J. C. Polkinghorne, Nucl. Phys.

B 12, 281 (1969).

[31] A. Accioly, Am. J. Phys. 65, 882 (1997).

[32] M. Lieberman and A. Lichtenberg, Principles of Plasma Discharges and Materials

Processing (Wiley, 1994).

[33] J. Jackson, Classical Electrodynamics, 3rd edition (Wiley, 1998).

[34] R. Santos, Mod Phys. Lett. A 26, 1909 (2011).

[35] J. Frenkel, Phys. Rev. E 54 5859 (1996).

[36] J. Frenkel and R. Santos, Int. J. Mod. Phys. B 13 315 (1999).

[37] S. Plimpton and W. Lawton, Phys. Rev. 50, 1066 (1936).

[38] A. Accioly and E. Scatena, Mod. Phys. Lett. A 25, 269 (2010).

[39] Para uma revisao geral, veja, por exemplo, L.-C. Tu, J. Luo and G. Gilles, Rep.
Prog. Phys. 68, 77 (2005).

[40] A. Accioly, P. Gaete, J. Helayel-Neto, E. Scatena, and R. Turcati, Mod. Phys.

Lett. A 26 1985 (2011).

[41] P. Frampton, Gauge Field Theories (Benjamin/Cummings, California, 1987).

52
[42] J. Schwinger, Phys. Rev. 73 416 (1948).

[43] R. Bouchendira et al Phys. Rev. Lett. 106 080801 (2011).

[44] D. Hanneke, S. Fogwell, and G. Gabrielse, Phys. Rev. Lett. 100 120801 (2008).

[45] Para uma interessante discussao sobre manipulacoes de matrizes de Dirac e espi-

nores independentes da representacao, veja P. Pal, arXiv:physics/0703214v2 .

[46] A. Ignatiev and G. Joshi, Phys. Rev. D 53 984 (1996).

[47] A. Accioly, P. Gaete, J. Helayel-Neto, E. Scatena, and R. Turcati, ar-

Xiv:1012.1045v2 [hep-th].

53
Captulo 2
Explorando a Gravitacao Massiva em 3D

Neste captulo, algumas propriedades interessantes do modelo de gravitacao massiva

de Bergshoe-Hohm-Townsend (BHT), tais como a presenca de uma forca gravita-

cional de curto alcance no limite nao-relativstico e a existencia de um angulo de

deexao dependente do parametro de impacto serao estudadas. Curiosamente, estes
fenomenos nao possuem contrapartida na gravitacao usual de Einstein em 3D. A

m de compreender melhor as duas propriedades gravitacionais citadas, elas tam-

bem serao analisadas no contexto da gravitacao de Einstein-Hilbert com derivadas

de ordem mais alta e com o sinal errado.

2.1 Motivacao

De acordo com uma lenda fsica, teorias de campos baseadas em acoes com derivadas de
quarta ordem (ou maiores) implicam em estados de ghosts nao-fsicos de norma nega-

tiva. Os argumentos apresentados na literatura para justicar esta situacao indesejada

tem, na verdade, evoludo somente para uma classe particular de tais teorias, a saber,
as teorias hbridas nas quais existem tanto derivadas de segunda quanto quarta ordem

[1]; e mais, tal discussao so tem sido formulada dentro do contexto da quantizacao cano-

54
nica [2]. Nao obstante, e notavel que dentre estas teorias hbridas exista uma classe de
modelos livres de ghosts em nvel de arvore. A primeira vista, parece que modelos de

curvatura pura, i.e., sistemas de gravitacao de quarta ordem com lagrangianas do tipo
R + R2 , pertencem a esta elite privilegiada. Na verdade, estes sistemas sao confor-

memente equivalentes a gravitacao de Einstein com um campo escalar massivo [3]. Em

outras palavras, apesar de possuir derivadas de quarta ordem na metrica, estes modelos

sao, em ultima instancia, de segunda ordem nas suas versoes tensoriais-escalares. Ha

apenas dois anos, contudo, Bergshoe, Hohm e Townsend (BHT) [412] propuseram
uma extensao particular da acao de Einstein-Hilbert (EH) com derivadas de ordem

mais alta em tres dimensoes espaco-temporais, a qual e livre de ghosts em nvel de

arvore [4]. A prova da unitariedade dada em [4] se baseia na existencia de uma forma

alternativa da acao de BHT envolvendo um campo auxiliar, cuja versao linearizada

e equivalente a acao tradicional de Fierz-Pauli, a qual sabe-se que nao possui ghosts;

alias, esta prova foi revista em [13]. Podemos tambem nos convencer da ausencia de
ghosts no modelo linearizado de BHT por meio de uma analise canonica [14]. O modelo
de BHT, e claro, viola tal folclore ja que e um exemplo de um sistema de quarta ordem

que nao possui ghosts.

Nossos objetivos neste captulo sao

1. Explorar algumas propriedades interessantes deste notavel modelo que nao possui

equivalente na gravitacao de Einstein usual em espaco-tempo tridimensional.

2. Esclarecer estas propriedades comparando-as com as analogas obtidas no contexto

da 3DHDG .

Um calculo detalhado de como e possvel obter o propagador concernente a gra-

3-Dimensional Higher Derivative Gravity

55
vitacao com derivadas de ordem superior e apresentado na Secao. 2.2, enquanto na
Secao. 2.3 demonstramos como a condicao de preservacao da unitariedade do modelo

recupera a lagrangiana concernente ao modelo BHT. Na Secao 2.4 e mostrado que,

no contexto do modelo de BHT, forcas gravitacionais de curto alcance agem sobre

partculas-teste se movimentando lentamente. Estas forcas atrativas sao uma carac-

terstica distinta da gravitacao massiva em 3D. A expressao analtica para o angulo

de deexao de um raio de luz no contexto do sistema mencionado e obtida na Se-

cao 2.5. Este angulo, ao contrario daquele na gravitacao de Einstein em 3D, depende
do parametro de impacto. Na Secao 2.6 sera feito um estudo comparativo entre os

resultados encontrados nas duas secoes anteriores e os correspondentes encontrados via

3DHDG. Esta analise e de grande ajuda para compreendermos melhor os resultados

das Secoes 2.4 e 2.5. Finalmente, na Secao 2.7, apresentamos alguns comentarios e
observacoes.

Utilizamos unidades naturais, c =  = 1, e a metrica de Minkowski e diag(+1, -1,

-1). O tensor de Ricci e denido por R = R + ....

2.2 O propagador da gravitacao massiva em (2+1)D

A gravitacao massiva em 3D, a qual e tambem conhecida como nova gravitacao mas-
siva e, na verdade, um caso particular da teoria tridimensional mais geral obtida
aumentando-se a gravidade planar com o termo de EH com o sinal errado por meio

de termos quadraticos da curvatura (3DHDG). Levando em conta que em tres dimen-

soes tanto o tensor de curvatura quanto o tensor de Ricci possuem o mesmo numero

de componentes [15], a lagrangiana para esta teoria generica pode ser escrita como

56
  
2 2 2
L= g 2 R + R + R , (2.1)
2 2

onde 2 = 32G, com G sendo o analogo 3D da constante de Newton, e e sao

coecientes livres. Note que as constantes , , e , possuem dimensao de massa

[] = 12 e [] = [] = 1, em unidades naturais.

2.2.1 Calculo de g, , R , R e R em primeira ordem
em

a densidade escalar g

Utilizando a aproximacao de campo fraco, ou seja,

g = + h , (2.2)

onde fazemos << 1 e = = diag(+1, 1, 1), obtemos tambem

g = h + 2 h h + ... . (2.3)

Denindo g det g e utilizando a relacao det A = etr ln A podemos calcular g da
seguinte forma

 1/2
g = (det g )1/2 = exp tr ln g (2.4)

= {exp tr ln [ ( + h )]}1/2 (2.5)

= {(det ) exp tr ln ( + h )}1/2 . (2.6)

57
Como << 1, podemos expandir ln ( + h ), obtendo

% 2 &1/2
g = exp tr(h h h + ...) (2.7)
2
1 2
= exp (h h h + ...) (2.8)
2 2
1 2 1 2
= 1 + ( h h h + ...) + ( h h h + ...)2 + ... . (2.9)
2 4 2 2 4

Deste modo, encontramos

2 2
g = 1 + h + h2 h h + ... . (2.10)
2 8 4

Smbolos de Christoel

Os smbolos de Christoel sao denidos como:

1
= g (g, + g, g, ). (2.11)
2
Substituindo a metrica da aproximacao de campo fraco (2.2), obtemos, em primeira

ordem em ,

= [h + h , h, ]. (2.12)
2 ,
Tensor de Riemman

O tensor de Riemann e denido como:

1
R = (g, + g, g, g, ). (2.13)
2

58
Novamente, utilizando (2.2), camos com:

R = (h, + h, h, h, ). (2.14)
2

Tensor de Ricci

O tensor de Ricci e encontrado pela contracao do tensor de Riemann com a metrica,

ou seja:

R = g R (2.15)

= (h, + h, h, h, ) (2.16)
2
 
= 2h (, + , ) , (2.17)
2

onde zemos
1
= h h, h h . (2.18)
2
Escalar de Curvatura

A contracao do tensor de Ricci nos da o escalar de curvatura R, ou seja

R = g R (2.19)
 
= 2h (, + , ) (2.20)
2
1 
= 2h , , (2.21)
2

com as mesmas denicoes (2.18).

Com isso, agora estamos aptos a calcular a reescrever a lagrangiana em (2.1).

59
2.2.2 Forma gama-gama da lagrangiana de Einstein

Uma forma simples de reescrever o termo livre da lagrangiana e coloca-lo na forma

gama-gama, ou seja, reescreve-lo como produtos simples de smbolos de Christoel.

Sendo assim, consideremos a seguinte acao:

 
 2R   2 
S1 = d x g 2 = d3 x g 2 g R ,
3
(2.22)

onde = 32G. O tensor R pode ser escrito como

R = , + , + . (2.23)

Vamos agora denir a seguinte densidade tensorial:

H gg , (2.24)

de modo que o produto H R se torna

H R = [H ], + H , + [H ], H ,

H ( ). (2.25)

Contudo, temos que a derivada covariante de qualquer densidade tensorial e identica-

mente nula, de modo que

H ; = H , + H + H H = 0, (2.26)

e encontramos

H , = H H H . (2.27)

Com isso, podemos reescrever a acao (2.22) como

60
   
2
S1 = 2 d3 x gg + , (2.28)

a qual so possui produtos de s.

2.2.3 A forma bilinear da lagrangiana da 3DHDG

Para encontrarmos o propagador da teoria em questao, podemos reescrever a lagran-

giana em sua forma bilinear (adicionada de uma lagrangiana xadora de gauge), da

seguinte forma
1
L = h O, h . (2.29)
2
Assim, invertendo o operador O, encontramos o propagador. De modo a facilitar o
trabalho, vamos separar a lagrangiana em quatro termos distintos, L = L1 + L2 + L3 +

Lgf , com

2
L1 = gR (2.30)
2
2
L2 = gR (2.31)
2
2
L3 = gR (2.32)
2
1  , 1 2
Lgf = h h, . (2.33)
2 2

Vamos agora reescrever cada uma das lagrangianas como funcoes explcitas de h .

Calculo de L1

61
Como vimos na secao anterior, podemos escrever L1 na forma gama-gama, ou seja

2 2  
L1 = gR = ( g)g (

+

)
2 2
%
2 1 
= 2 ( g)g g (g, + g, g, )g (g, + g, g, )
4 &
+ g  (g, + g, g, )g (g, + g, g, ) .

Utilizando (2.2), encontramos

1 %
L1 = ( g)g [(h, + h , h, )(h , + h, h, )
2 &

+ (h, + h , h, )(h , + h , h, )] ,

e, rearranjando os termos de maneira conveniente, utilizando (2.3) e (2.10) obtemos

 1 
L1 = h + ( 2 2) h .

(2.34)
2

Como o operador O, em (2.29) e simetrico na troca de por , por e por

, vamos simetrizar o termo entre colchetes de (2.34):

1  1
L1 = h ( + + + ) + ( + )
2 2
1 
+ ( 2 + 2) 2 h .

2

Finalmente, passando para o espaco dos momenta atraves de uma transformacao de

1
Fourier do tipo L(k) = (2)3
eikx d3 x, encontramos:

1 1
L1 = h (k k + k k + k k + k k ) ( k k + k k )
2 2
k 2 
( + ) + k 2 h . (2.35)
2

62
 Calculo de L2
Utilizando (2.21) e (2.10), podemos reescrever L2 como

2 2 1 1 2
L2 = gR = g 2h (h h)
2 2 2 2
2 
= h 2 2 + 22 h .

2

Novamente, simetrizando o termo entre colchetes,

1  
2
L2 = h ( + )2 + 22 h ,

(2.36)
2

que no espaco dos momenta, pode ser escrito como

1  
L2 = h 2 k k k k k 2 ( k k + k k ) + k 4 h . (2.37)
2

Calculo de L3

Utilizando o tensor de Ricci (2.17), L3 se torna

2 %  1 1 &2
L3 = gR = 2 2h (h h) + (h h) ,
2 8 2 2

rearranjando os termos, camos com

2  
L3 = h 2 2 2+ 22+ 22 2 2 h ,
8

que de modo simetrizado, se torna

2  1
L3 = h 2 ( + + + )2
8 2
1 
+ ( + )22 + 22 ( + )2 h .
2

63
Assim, reescrevendo no espaco dos momenta, obtemos

1 2  k 2
L3 = h 2k k k k ( k k + k k + k k + kk )
2 4 2
k 4 
+ ( + ) + k 4 k 2 ( k k + k k ) h . (2.38)
2

Calculo de Lgf
Como vimos, no gauge de De Donder, a lagrangiana xadora de gauge e escrita como

1  , 1 2
Lgf = h h,
2 2
1 1 
= h + 2 h ,

2 4

cuja versao simetrizada e

1  1 1
Lgf = h ( + + + ) + ( + )
2 4 2
1 
2 h .
4

Assim, no espaco dos momenta, Lgf se torna

1 1
Lgf = h (k k + k k + k k + k k )
2 4
1 k2 
(k k + k k ) + h . (2.39)
2 4

Em termos dos operadores de Barnes-Rivers [16, 17] (veja o Apendice A), o operador

64
O na equacao (2.29) se torna

  k 2 % (1) 1
O, = k 2 P (0s) P 2 + P + P (0s) + P (0w)
, 2 2
&
2  (0sw) 
P + P (0ws) + 22 k 4 P (0s)
2 ,
2 k 4  (2) 
+ P + 3P (0s) , . (2.40)
4

Agora podemos utilizar o metodo descrito no Apendice A para encontrar o inverso do

operador de onda O, que nos dara o propagador que estamos procurando. Ele sera
escrito como

M22 (2) 2 (1) M02

P = P + P P (0s) +
k 2 (k 2 M22 ) k2 k 2 (k 2 M02 )

M02 2  (0sw) (0ws)
 2 [2k 2 M02 (2 + 1)] (0w)
2 2 P +P + P ,(2.41)
k (k M02 ) k 2 (k 2 M02 )

4 4
onde M22 = 2
(> 0), M02 = 2 (8+3) (> 0), e e um parametro de gauge. No-
tamos claramente que o sistema 3DHDG e nao-unitario, mesmo em nvel de arvore.
Chamamos atencao ao fato de que o propagador acima e um caso especial daquele

concernente a supergravidade com N = 1 e N = 2 em tres dimensoes, aumentada por

termos do tipo (curvatura)2 e que foi calculado por Nishino e Rajpoot [18]; ele pode

tambem ser obtido por meio do propagador relacionado a gravitacao 3D com o sinal
errado aumentado por um termo de Chern-Simons [1921].

2.3 Unitariedade em nvel de arvore do modelo 3DHDG

Nosso proximo passo sera determinar condicoes sobre os parametros e na lagran-

giana (2.1) de modo que a mesma que unitaria e, deste modo, recuperando o modelo

65
BHT. Para tanto, faremos uso dum procedimento padrao que transforma a tarefa de
checar a unitariedade em nvel de arvore de um dado modelo, o que geralmente e um

trabalho exaustivo, num simples exerccio algebrico [22]. A prescricao consiste basica-
mente em saturar o propagador com correntes externas conservadas, compatveis com

as simetrias do sistema, e examinar posteriormente os resduos do propagador satu-

rado (PS) em cada polo simples. Contraindo o propagador em (2.41) com correntes

conservadas T (k), (k T = k T = 0), obtemos

  
1 1 1 2 1 2 1 1 1 1 2
PS = 2 2
T T + 2+ 2 2
T , (2.42)
k 2 k m2 2 k k m0 2

Supondo que nao existam taquions no modelo, encontramos imediatamente as se-

guintes condicoes

1
> 0, < 0. (2.43)
8 + 3
Por outro lado, os resduos do PS nos polos k 2 = m22 , k 2 = 0, and k 2 = m20 sao,
respectivamente

 1 2
2
Res(PS) |k2 =m2 = T T |k2 =m2 , (2.44)
2
 2  2

Res(PS) |k2 =0 = T 2 T
2
|k2 =0 , (2.45)
1
Res(PS) |k2 =m2 = T 2 |k2 =m2 . (2.46)
0 2 0

Agora, como sabemos, a unitariedade em nvel de arvore de um modelo generico e

garantida se o resduo em cada polo simples de PS e 0. Mantendo em mente que
   
2 1 2 2 2
T 2 T |k2 =m2 > 0 e T T |k2 =0 = 0 (veja o Apendice C), chegamos a conclu-
2

66
sao de que: (i)Res(PS) |k2 =m2 > 0 implica em > 0 e < 0 e, (ii) Res(PS) |k2 =0 = 0.
2

Consequentemente, nao precisamos nos preocupar com estes polos; o problema aparece

em k 2 = m20 ja que Res(PS) |k2 =m2 < 0. Um modo de contornar esta situacao e con-
0

siderar o limite m0 do modelo em discussao, o que nos leva a conclusao de que

= 38 . Portanto, a classe de modelos denidos pela densidade de lagrangiana

  
2R 2 3
L= g 2 + R R2 , (2.47)
2 8
4
e livre de ghosts em nvel de arvore. Por conveniencia, substituindo por 2 m22
, onde

m2 e um parametro de massa, encontramos

  
2R 2 2 3
L= g 2 + 2 2 R R2 , (2.48)
m2 8
que e exatamente a densidade de lagrangiana do modelo BHT.
Este sistema, como ja comentamos, e unitario em nvel de arvore. Nao obstante, nao
esta claro se esta razao entre e sobrevivera a renormalizacao em loops arbitrarios,

ou mesmo em um loop; em outras palavras, unitariedade alem do nvel de arvore

deve ser checada [11]. Na verdade, e mais provavel que o modelo de BHT seja nao-

renormalizavel ja que ele melhora somente o comportamento das projecoes do spin-2

do propagador, mas nao as projecoes do spin-0 [23].

2.4 Potencial inter-partculas

Na mecanica quantica, a nocao de potencial e introduzida para ser utilizada com a

equacao de Schrodinger. Ja na teoria quantica de campos, as interacoes sao visualizadas

como advindas da troca de quanta. E notavel que estes dois conceitos possam ser
mesclados para partculas que se movam lentamente, se nos restringirmos a processos

67
de ordem mais baixa. De fato, de acordo com a aproximacao de Born, a secao de
choque diferencial que concerne, por exemplo, o espalhamento de dois bosons escalares
 d  m 2
=  4 ei(p p) U (r)d2 r , onde p e p sao,

de mesma massa m e dada por d CM

respectivamente, os momenta de entrada e sada das partculas. Por outro lado, a

secao de choque para a interacao dos dois citados bosons via troca de gravitons pode ser
 d   2
expressa como d CM
= M NR 
16m
, onde MNR e o limite nao-relativstico da amplitude

de Feynman para o processo que estamos discutindo. Como resultado, o potencial

inter-partculas assume a forma


1 1
U (r) = d2 kMNR eikr , (2.49)
4m (2)2
2

onde k e o momentum do graviton trocado.

Utilizando a equacao acima, iremos calcular em seguida a expressao para a ener-

gia potencial gravitacional no contexto do modelo de BHT [24]. Consideremos entao

a seguinte lagrangiana descrevendo a interacao gravitacional entre bosons escalares

neutros
1 
Lesc. = m2 2 g (2.50)
2
Novamente, na aproximacao de campo fraco, a lagrangiana acima se torna
 
1 2 2
 h
Lesc. = m

1+
2 2
 
1  2 2
 h
= ( h ) m

1+
2 2
 
1 2 2
 h
= h m

1+ . (2.51)
2 2

A lagrangiana de interacao contera somente os termos em h . Da equacao acima,

68
encontramos

m2 2
Lint. = h +


2 2

1  2 2

Lint. = h m . (2.52)
2 2

Escrevendo o termo entre colchetes no espaco dos momenta e simetrizando os ndices,

encontramos o vertice elementar para o processo que acabamos de descrever

1  
(p, p ) = p p + p p (p p + m2 ) , (2.53)
2

onde supomos que os momenta estao entrando, e multiplicamos por um fator 2, ja que

estamos lidando com bosons.

Figura 2.1: Vertice relevante para a interacao gravitacional de bosons escalares

A amplitude de Feynman e, entao

M = (p, p )P, (q, q  )

m22 1
= 2
A + 2 B, (2.54)
k (k m2 )
2 2 k

69
onde o propagador e dado por

2 (1) m22 (2) 1 (0s) 2  (0sw) 
P = P + P + P + P + P ows
k2 k 2 (k 2 m22 ) k2 k2
 
4 2
+ 2
+ 2 P (0w) , (2.55)
k k

1 2% 1 1
A = (p q)(p q  ) + (p q  )(p q) 2 (q q  )(m2 p p )2 2 (p p )(m2 q q  )2
2 k k
1
+ 2 [(p q) (p q  ) (p q) + (p q  )] (m2 p p )(m2 q q  ) (p p )(q q  )
k
1 &
+ 4 (m2 p p )2 (m2 q q  )2 , (2.56)
k

1 2 1
B = (p p )(q q  ) + (q q  )(m2 p p ) + (p p )(m2 q q  ) + 2 (q q  )(m2 p p )2
2 k
1 1
+ 2 (p p )(m2 q q  )2 + (m2 q q  )(m2 p p ) + 2 (m2 q q  )(m2 p p )2
k k
1 1 
+ 2 (m2 q q  )2 (m2 p p ) + 4 (m2 p p )2 (m2 q q  )2 . (2.57)
k k

No limite nao-relativstico, Eq. (2.54) se reduz a

1 m4
MNR = 2 2 . (2.58)
2 k + m22

Realizando a integracao apropriada utilizando as Eqs. (2.49) e (2.58), chegamos a

conclusao de que a energia potencial e dada por

70
 U (r) = 2m2 GK0 (m2 r), (2.59)

onde K0 e a funcao de Bessel modicada de ordem zero. Note que K0 (x) se comporta
1
como -lnx na origem e como x 2 ex assintoticamente. O potencial nao-relativstico,

por sua vez, tem a forma

V (r) = 2mGK0 (m2 r). (2.60)

Neste ponto, alguns comentarios sao pertinentes:

Assim como o potencial newtoniano, VNewt (r) = 2mG ln rr0 (onde r0 e um regula-

dor infravermelho), o potencial concernente ao modelo BHT possui uma singula-

ridade logartmica na origem.

V (r) 0 conforme r .

No limite em que m2 , V (r) 0, reproduzindo, deste modo, o potencial

usual da gravitacao tridimensional.

Por outro lado, no limite nao-relativstico, uma partcula teste de massa mtest se
movendo num campo gravitacional fraco experiencia uma forca F(r) = mtest V (r).

Agora, levando a Eq. (2.60) em conta, encontramos que F(r) = F (r)r, com

F (r) 2Gmtest mm2 K1 (m2 r), (2.61)

onde K1 e a funcao de Bessel modicada de ordem 1. Relembrando que xK1 (x) e

uma funcao negativa monotonicamente crescente no intervalo 0 x < ja que

71
 d
dx
[xK1 (x)] = xK0 (x) chegamos a conclusao de que F(r) e sempre atrativa (veja a
Fig. 2.2). Esta forca gravitacional de curto alcance nao existe no contexto da relativi-

dade geral tridimensional, sendo peculiar apenas ao modelo de BHT.

A energia potencial dada pela Eq. (2.59) e, como esperado, no nvel linearizado,

equivalente aquela obtida atraves do modelo (unitario) de Pauli-Fierz para um campo

massivo de spin-2.

Figura 2.2: A forca gravitacional de curto alcance relacionada ao modelo de BHT.

2.5 Desvio da Luz

A gravitacao 3D de Einstein e trivial fora das fontes; consequentemente, nao devera-

mos esperar, ao menos em princpio, nenhuma deexao da luz neste contexto. Mas nao

e bem assim na pratica. Na verdade, um raio de luz sofre uma deexao que e propor-
cional a fonte de massa; alem disso, o encurvamento e independente do parametro de
impacto. Isto pode ser mostrado facilmente no caso de uma partcula puntual massiva

m localizada em r = 0. A metrica linearizada correspondente, no gauge de De Donder,

72
e dada por

ds2 = dt2 (1 )(dr2 + r2 d2 ),

onde = 8Gmln rr0 , com r0 sendo um regulador infravermelho, e r e sao as coorde-

nadas polares usuais.

Introduzindo agora uma nova coordenada radial (r ) e angular ( ) atraves da mu-

danca de coordenadas

(1 )r2 = (1 8Gm)r2 ,  = (1 4Gm),

obtemos, em ordem linear em Gm,

ds2 = dt2 dr2 r2 d2 . (2.62)

Curiosamente, a geometria em torno da partcula puntual e localmente identica

aquela de um espaco-tempo plano; contudo, ela nao e globalmente minkowskiniana

ja que o angulo  varia no interalo 0  < 2(1 4Gm). Portanto, a metrica

tridimensional (2.62) descreve um espaco conico com uma cunha de tamanho angular
igual a 8Gm removida e com as duas faces da cunha identicadas. O angulo de

deexao da luz e, portanto, igual a 4Gm e e independente do parametro de impacto.

Este fenomeno e muito similar ao efeito Aharanov-Bohm na eletrodinamica. De fato,

la o padrao de interferencia entre dois feixes coerentes de eletrons se movendo em

lados diferentes de um no solenoide e afetado pelo uxo magnetico dentro do mesmo,
apesar do campo magnetico ser nulo ao longo dos feixes; no nosso caso, por outro lado,

73
o raio de luz se propagando na regiao plana do espaco-tempo em torno da partcula
puntual sente a curvatura do espaco-tempo, a qual e restrita a partcula puntual massiva

[25, 26]. Como veremos em seguida, estes estranhos efeitos geometricos nao acontecem
na gravitacao da 3DHDG. Neste sentido, vamos calcular a deexao gravitacional de

um raio de luz no ambito do modelo de BHT. Para tanto, temos que calcular antes as
equacoes de campos linearizadas relacionadas ao sistema de BHT.

As equacoes de campos concernentes a densidade de lagrangiana

  
2R 2 2 3
L= g 2 + 2 2 R R2 LM , (2.63)
m2 8
onde LM e a densidade de lagrangiana para a materia, sao

 
1 1 2 1 1 3 2 3
G + 2 R g R 2R R g 2R + 2R R g + RR

m2 2 4 4 16 4
2
= T , (2.64)
4

onde T e o tensor de energia-momento, e G R 12 g R e o tensor de Einstein.

As equacoes de campo linearizadas sao dadas por

    
2 1 (lin) 1  T
1+ 2 2h + R + + = T , (2.65)
m2 2 4 2 4 2

1   
2 R(lin)
onde R(lin) = 2
2h
, , 1+ m22
+ 4m22
, h 12 . Note

que aqui os ndices sao levantados (abaixados) utilizando ( ).

Uma exposicao detalhada dos calculos concernentes as equacoes de campo, sua linearizacao e os
resultados no gauge de Teyssandier e apresentada no Apendice B, assim como a solucao das equacoes
para o caso de uma partcula puntual.

74
 Analogamente ao trabalho de Teyssandier em gravitacao com derivadas mais alta
em 4D [27], pode ser demonstrado que e sempre possvel encontrar um sistema de

coordenadas no qual as condicoes de gauge, = 0, na metrica linearizada, sao validas.

Considerando que essas condicoes sao satisfeitas, e imediato mostrar que a solucao geral

da Eq. (2.65) e dada por

h = h(E)
, (2.66)

(E)
onde h e a solucao da equacao de Einstein linearizada no gauge de De Donder, i.e.,

2h(E)
=
(E)
(T n T ), = 0, (2.67)
2

(E) (E)
onde h 12 h(E) , enquanto satisfaz a equacao

1
(2 + m22 ) = (T T ). (2.68)
2 2
(E)
E importante notar que, precisamente neste gauge, as equacoes para e h sao

totalmente desacopladas. Como resultado, a solucao geral da equacao (2.65) se reduz

a uma combinacao linear das solucoes anteriores.

Resolvendo as Eqs. (2.67) e (2.68) para uma partcula puntual de massa m locali-
zada em r = 0, encontramos

m
h00 = K0 (m2 r) (2.69)
8
m  r
h11 = h22 = K0 (m2 r) + 2 ln . (2.70)
8 r0

75
Note que para m2 , a solucao acima reproduz a solucao da equacao de Einstein
em tres dimensoes no gauge de De Donder com o sinal errado, como era de se esperar.

Agora estamos prontos para discutir a deexao de um foton devido ao campo gra-
vitacional gerado por uma massa m. Suponhamos um foton com momentum p vindo

do innito com um parametro de impacto b (veja a Fig. 2.3).

m x

Figura 2.3: Geometria da deexao da luz.

A variacao de momentum em p quando o foton passa pelo campo gravitacional e

dada por



p = p h dx , (2.71)
2

onde a integral e realizada ao longo da linha reta aproximada da trajetoria do foton.

Como consequencia, o deslocamento ao longo da linha reta aproximada e o momentum

sao, respectivamente,

   
dx dx1 , dx1 , 0 , p p1 , p1 , 0 .

Inserindo essas quantidades na Eq. (2.71), obtemos

76
 

p2 = p1 y [h00 + h11 ]dx1 , (2.72)
2

a qual pode ser escrita como


 
d r
p2 = p1 (h00 + h11 ) dx. (2.73)
2 dr y y=b

Com os resultados de (2.69) e (2.70), podemos escrever a Eq. (2.73) simplesmente

como

   

bm2 1 m 2 K 1 m 2 x2 + b 2
p2 = p
4 0 x2 + b2

1
dx.
x2 + b2

Agora, levando em conta que [28]

  '  x2+1
K x2 + y 2 ' dx
(x2 + y 2 )
2 ( + 1)
= K1 (y), > 0, Re() > 1,
+1 y 1

e

z
K+ 1 (z) = e ,
2 2z
encontramos prontamente

2 mp1  
p2 = 1 + em2 b . (2.74)
8

77
   
 py 
Chamando de =  px  o valor absoluto do angulo de deexao, chegamos a

conclusao de que

 
= 4Gm 1 em2 b . (2.75)

Chamamos atencao ao fato de que para m2 , o modulo do angulo de deexao

se reduz a 4Gm, o qual, como ja havamos comentado, e a previsao geometrica da

gravitacao 3D de Einstein [24].

2.6 Analise dos resultados anteriormente obtidos

no contexto dos modelos BHT e 3DHDG

Ja que, como havamos discutido, o modelo de BHT e um caso especial do sistema

3DHDG, uma analise comparativa entre as previsoes destes modelos para o potencial

nao-relativstico certamente nos esclarecera os resultados encontrados nas Secoes 2.4 e

2.5. Para tanto, iremos analisar brevemente estes fenomenos no contexto da 3DHDG

e, em seguida, realizar um estudo comparativo entre os fenomenos correspondentes

obtidos para o modelo de BHT. Vale notar que o estudo da energia potencial no contexto
do modelo 3DHDG tem sua importancia ja que existe uma incrvel semelhanca entre

este potencial e aquele concernente a interacao de dois vortices supercondutores.

2.6.1 O potencial nao-relativstico para o modelo 3DHDG

Um calculo similar aquele da Secao 2.4 leva a seguinte expressao para a energia poten-

cial

78
  
U (r) = 2GM 2 K0 (M0 r) K0 (M2 r) . (2.76)

Ao contrario do potencial newtoniano, UNewt (r) = 2GM 2 ln rr0 , o qual possui uma

singularidade logartmica na origem e nao e limitado no innito, este potencial e ex-

tremamente bem comportado: ele e nito na origem (U (0) = 2GM 2 ln M
M0
2
) e zero no

innito. Curiosamente, ele descreve tres possveis regimes gravitacionais. Assim, de-
pendendo da escolha dos parametros na acao, podemos obter regimes de gravitacao
(M2 < M0 ), anti-gravitacao (M0 < M2 ) ou blindagem gravitacional (M0 = M2 ).

Figura 2.4: Possveis regimes gravitacionais preditos pela 3DHDG.

Na Fig. 2.4 e mostrado um esboco esquematico do potencial gravitacional para as

tres situacoes que acabamos de descrever. E bastante notavel que este esboco quan-

titativo seja muito similar, mutatis mutandis, aquele encontrado por Jacobs e Rebbi
por meio de um estudo computacional aproximado do potencial de dois vortices [29].
De fato, calculando a energia de interacao entre dois vortices, na teoria de Ginzburg-

Landau ou, equivalentemente, no modelo abeliano de Higgs, como funcao da separacao

entre os dois vortices e o parametro , o qual mede a intensidade relativa do auto-

79
acoplamento da materia e o acoplamento eletromagnetico, eles mostraram que o po-
tencial entre os vortices e atrativo para < 1 e repulsivo para > 1. Eles tambem

encontraram que, para = 1, um limite inferior e atingido para todas separacoes e,

portanto, os vortices nao interagem.

2.6.2 Deexao da luz no contexto da 3DHDG

Um procedimento analogo aquele empregado na Secao 2.5 fornece a solucao geral para
as equacoes de campo linearizadas, relacionadas ao modelo 3DHDG, geradas por uma

partcula de massa m localizada em r=0, i.e.,

m  
h00 = K0 (M0 r) K0 (M2 r) , (2.77)
8
m  r
hii = K0 (M0 r) K0 (M2 r) 2 ln . (2.78)
8 r0

Utilizando estes resultados, podemos mostrar que o valor do angulo de deexao, ,

de um foton devido ao campo gravitacional da massa m, e dado por

 
= 4Gm 1 eM2 b . (2.79)

Notem que a excitacao escalar de massa M0 nao contribui em nada para o desvio

da luz. Isso se deve ao fato de que a metrica que diz respeito a gravitacao R +
R2 linearizada com o termo de EH com o sinal errado a teoria obtida atraves da
 
linearizacao das equacoes de campos relacionadas a acao d3 x g 22 R + 2 R2 LM
e conformemente relacionada aquela versao 3D linearizada da relatividade geral com
o sinal errado. Alias, chamando a solucao da gravitacao R + R2 linearizada com o sinal

80
 (R+R2 )
errado de g , obtemos prontamente

(R+R 2) 2)  
g + h(R+R
= (1 ) h(E)

(E)
= (1 )g ,

(E)
onde, e claro, termos de ordem 2 foram descartados. Aqui, h e possuem o mesmo
(E)
sentido que na Secao 2.5, enquanto g e a solucao para a gravitacao 3D de Einstein

linearizada com o sinal errado. Esta e a razao pela qual as expressoes (2.75) e (2.79),
com as devidas mudancas, sao tao similares.

2.6.3 Um estudo comparativo

Agora estamos prontos para realizar uma analise comparativa entre estes dois resultados
obtidos para o sistema 3DHDG e os correspondentes previstos pelo modelo BHT, de

modo a ganharmos um melhor entendimento sobre este ultimo.

Das consideracoes anteriores, chegamos a conclusao de que o nascimento do poten-

cial de BHT vem do potencial da 3DHDG com M2 < M0 (veja a Fig. 2.5). Examinando
o diagrama mostrado na Fig. 2.5, podemos ver claramente que conforme M0 se torna
cada vez maior, o potencial da 3DHDG rapidamente se aproxima do potencial de BHT

ate que, eventualmente, eles se fundem. Portanto, para chegar ao potencial de BHT
partindo da 3DHDG, este ultimo deve necessariamente se tornar singular na origem, o

que ocorre no limite de M0 . E interessante notar que esta e justamente a condicao

necessaria para evitar, em nvel de arvore, o ghost massivo de spin-0 que assombra a

3DHDG. Portanto, a presenca da singularidade no potencial de BHT e, de certo modo,

relacionada a ausencia do ghost em nvel de arvore; em outras palavras, a unitariedade

81
em nvel de arvore e a existencia de uma singularidade no potencial estao, de algum
modo, relacionadas.

Figura 2.5: Energia potencial gravitacional para os modelos 3DHDG com M2 < M0 (linha
contnua) e BHT (linha pontilhada).

Vamos agora melhorar a nossa discussao sobre o desvio da luz no contexto do modelo
BHT. Na Secao 2.5, chegamos a conclusao de que o angulo de deexao gravitacional

previsto para este modelo e independente do parametro de impacto, como era de se

esperar; contudo, este angulo e pequeno? Antes de responder esta questao, lembramos

que uma das razoes para realizar uma pesquisa em gravitacao tridimensional e que a
gravitacao 3D possui uma relevancia fsica direta para modelar fenomenos que estao,
na verdade, dinamicamente connados a uma dimensionalidade menor. De fato, a

fsica gravitacional na presenca de uma corda cosmica retilnea (innitamente longa,

perpendicular a um plano) e adequadamente descrita pela gravitacao 3D [26]. Isto

signica que a resposta para a nossa questao inicial pode ser encontrada analisando a

82
deexao gravitacional em um modelo 4D analogo, de certa maneira, ao modelo BHT,
ou seja, a gravitacao com derivadas de ordem mais alta em 4D com o sinal errado,

denida pela lagrangiana


2R 2 2
L= g 2 + R + R .
2 2

Pode-se mostrar que o modulo do angulo de deexao, , e dado por

 
= 4G 1 eM 2 b ,

2
onde e a densidade de energia da corda, e M 2 24 (> 0). Se levarmos em conta que
|| = 1060 encontramos imediatamente M2 103 cm1 [30, 31]. Como consequencia,
para distancias longe da corda, o angulo de deexao ( G) e pequeno, como era de

se esperar. De fato, obervacoes limitam G 105 , enquanto o cenario de cordas para

a formacao de galaxias requer G 106 , correspondendo a uma escala de grande

unicacao de 1016 GeV [26]. Como resultado, apenas para fotons passando muito
proximos a corda cosmica a deexao da luz ira diferir signicativamente do valor usual

previsto pela gravitacao com derivadas de ordem mais alta em 4D com o sinal errado.
A mesma conclusao, com as mudancas apropriadas, se aplicam ao modelo BHT.

2.7 Discussao

Como se sabe, a gravitacao tridimensional de Einstein sem fontes e sicamente trivial

ja que os tensores de Einstein e Riemann sao equivalentes em D = 3. Alem disso, a
quantizacao do campo gravitacional nao da origem a gravitons propagantes ja que a

83
metrica do espaco tempo e localmente determinada pelas fontes. Consequentemente, a
descricao de fenomenos gravitacionais via gravitacao 3D nos leva a alguns resultados

bizarros, tais quais

Ausencia de forca gravitacional no limite nao-relativstico.

Deexao gravitacional independente do parametro de impacto.

Nao obstante, como ja discutimos, estes estranhos fenomenos nao existem no con-
texto do modelo BHT. Alias, neste contexto, forcas gravitacionais de curto alcance sao

exercidas sobre partculas se movendo lentamente; alem disso, o desvio da luz depende
do parametro de impacto, como deveria. Em outras palavras, as idiossincrasias da

relatividade geral em 3D nao aparecem no sistema BHT.

O grande incentivo para o estudo de modelos gravitacionais em 3D e, na verdade,

tentar encontrar um sistema gravitacional com divergencias ultravioletas menos severas

na teoria de pertubacao. Como a relatividade geral em 3D e dinamicamente trivial,
o modelo BHT, o qual e unitario em nvel de arvore, e um passo importante nesta

direcao. Este tipo de pesquisa conduzida em dimensoes menores certamente nos ajuda
a compreender melhor algumas difceis questoes conceituais, as quais estao presentes e

mais opacas no mundo fsico quadridimensional.

Vale a pena mencionar que a trivialidade da relatividade geral em 3D pode tambem

ser resolvida adicionando a acao de EH em 3D um termo de Chern-Simons, que viola

a paridade. O modelo resultante e conhecido como topological massive gravity (TMG)
[32]. Recentemente, Ahmedov e Aliev [12] mostraram que as equacoes de campo do

modelo BHT consistem de uma equacao massiva (tensorial) do tipo Klein-Gordon com
uma fonte de curvatura quadratica e uma equacao de vnculo; e mais, para espaco-

tempo algebricos do tipo D e N , as equacoes de campos da TMG podem ser pensadas

84
como as razes da equacao massiva do tipo Klein-Gordon, as quais permitiriam ma-
pear todas as solucoes do tipo D e N da TMG naquelas do modelo BHT. Chamamos

atencao ao fato de que, contudo, em contraste a TMG, a gravitacao massiva 3D possui

a grande vantagem de ser uma teoria que preserva a paridade.

85
Referencias Bibliogracas

[1] K. Stelle, Phys. Rev. D 16 953 (1977); Gen. Relativ. Gravit. 9 353 (1978).

[2] P. Manheim, Found. Phys. 37 532 (2007).

[3] B. Whitt, Phys. Lett. B 145 176 (1984).

[4] E. Bergshoe, O. Hohm, and P. Townsend, Phys. Rev. Lett. 102 201301 (2009).

[5] E. Bergshoe, O. Hohm, and P. Townsend, Phys. Rev. D79 124042 (2009).

[6] R. Andringa et. al, Class. Quant. Grav. 27 025010 (2010).

[7] E. Bergshoe, O. Hohm, and P. Townsend, Annals Phys. 325 1118 (2010).

[8] E. Bergshoe et. al, Class. Quant. Grav. 28 015002 (2011).

[9] E. Bergshoe et. al, Class. Quant. Grav. 27 235012 (2010).

[10] M. Nakasone and I. Oda, Prog. Theor. Phys. 121 1389 (2009).

[11] I. Gullu and B. Tekin, Phys. Rev. D 80 064033 (2009).

[12] H. Ahmedov and A. Aliev, Phys. Rev. Lett. 106 021301 (2011).

86
[13] E. Bergshoe, O. Hohm, and P. Townsend, Journal of Physics: Conference Series
229 012005 (2010).

[14] S. Deser, Phys. Rev. Lett. 103 101302 (2009).

[15] A. Staruszkiewicz, Acta Phys. Pol. 24 734 (1963).

[16] P. Nieuwenhuizen, Nucl. Phys. B 60 478 (1973).

[17] I. Antoniadis and E. Tomboulis, Phys. Rev. D 33 2756 (1986).

[18] H. Nishino and S. Rajpoot, Phys. Lett. B 639 110 (2006).

[19] A. Accioly, Phys. Rev. D 67 127502 (2003).

[20] A. Accioly, Nucl. Phys. B (Proc. Suppl.) 127 100 (2004).

[21] A. Accioly and M. Dias, Int. J. Theor. Phys. 44 1123 (2005).

[22] A. Accioly, J. Helayel-Neto, E. Scatena, J. Morais, R. Turcati, and B. Pereira-Dias,

Class. Quantum Grav. 28 225008 (2011).

[23] E. Bergshoe, O. Hohm, and P. Townsend, Cosmology, The Quantum Vacuum

and Zeta Functions: A Workshop with a Celebration of Emilio Elizaldes Sixtieth

Birthday 8-10 Mar 2010, Bellaterra, Barcelona, Spain.

[24] A. Accioly, J. Helayel-Neto, J. Morais, E. Scatena, and R. Turcati, Phys. Rev. D

83 104005 (2011).

[25] A. Vilenkin, Phys. Rev. D 23 852 (1981).

[26] A. Vilenkin and E. Shellard, Cosmic Strings and Other Topological Defects (Cam-

bridge University Press, Cambridge, 2001).

87
[27] P. Teyssandier, Class. Quantum Grav. 6 219 (1989).

[28] I. Gradshteyn and I. Ryzhic, Table of Integrals, and Products (Academic Press,

San Diego, 1994).

[29] L. Jacobs and C. Rebbi, Phys. Rev. B 19 4486 (1979).

[30] A. Accioly and H. Blas, Phys. Rev. D 64 067701 (2001).

[31] A. Accioly and H. Blas, Int. J. Theor. Phys. 43 47 (2004).

[32] S. Deser, R. Jackiw, and S. Templeton, Phys. Rev. Lett. 48 975 (1982); Ann. Phys.

(N.Y.) 140 372 (1982); 185 406 (E) (1988).

[33] A. Accioly, A. Azeredo, H. Mukay, and E. R. Neto, Prog. Theor. Phys. 104 103

(2000).

88
Apendice A
Uma prescricao para o calculo do
propagador do graviton, e uma lista de
algumas identidades que facilitam
enormemente esta tarefa

Para encontrarmos o propagador relacionado a densidade lagrangiana na Eq. (2.1), e

muito conveniente trabalhar em termos dos operadores de Barnes-Rivers no espaco dos
tensores simetricos de posto dois. O conjunto completo de operadores tridimensionais

no espaco dos momenta e [16, 17]

89
 (2) 1
P, = ( + ), (A.1)
2
(1) 1
P, = ( + + + ), (A.2)
2
(0s) 1
P, = , (A.3)
2
(0w)
P, = , (A.4)
(0sw) 1
P, = , (A.5)
2
(0ws) 1
P, = , (A.6)
2
k k k k
onde k2
and k2
sao, respectivamente, os projetores transversais

e longitudinais usuais. A tabela multiplicativa para estes operadores e mostrada na

Tabela I.

Tabela A.1: Tabela multiplicativa para os operadores de Barnes-Rivers

P (2) P (1) P (0s) P (0w) P (0sw) P (0ws)

(2)
P P (2) 0 0 0 0 0
P (1) 0 P (1) 0 0 0 0
P (0s) 0 0 P (0s)
0 P (0sw)
0
P (0w) 0 0 0 P (0w) 0 P (0ws)
P (0sw) 0 0 0 P (0sw) 0 P (0s)
P (0ws) 0 0 P (0ws)
0 P (0w)
0

Para calcular o propagador do graviton nos precisamos da parte bilinear da densi-

1
dade de lagrangiana (2.1). Com o xador de gauge ( )2
2
(gauge de De Donder),
e indo para o espaco dos momenta, reproduzimos a Eq. (2.40). A tarefa de calcular o

operador O e imensamente facilitada se utilizarmos as seguintes identidades

90
   1
P (2) + P (1) + P (0s) + P (0w) , = ( + ),
 2 
(0s) (0w)
(0sw) (0ws)
2P +P + 2(P +P ) = ,
,
  1
2P (1) + 4P (0w)
,
= 2 ( k k + k k + k k + k k ),
k
  1
2P (0w) + 2(P (0sw) + P (0ws) ) = 2 ( k k + k k ),
, k
(0w) 1
P, = 4 (k k k k ).
k

Agora, se escrevermos o operador O na forma generica

O = x1 P (1) + x2 P (2) + xs P (0s) + xw P (0w) + xsw P (0sw) + xws P (0ws) ,

e levarmos em conta que OO1 = I, onde O1 e o propagador, obtemos prontamente

1 (1) 1 (2) 1 
1
O = P + P + xw P (0s) + xs P (0w)
x1 x2 xs xw xsw xws
xsw P (0sw) xws P (0ws) . (A.7)

De (A.7) e (2.40) obtemos (2.41).

91
Apendice B
Equacoes de campo do modelo BHT e o
gauge de Teyssandier

B.1 Equacoes de campo do modelo BHT

As equacoes de campos concernentes a acao

    
3 3 2R 2 2 3 2
S= d xL = dx g 2 + 2 2 R R LM , (B.1)
m2 8
onde LM e a densidade de lagrangiana para a materia, podem ser encontradas por meio

da variacao da lagrangiana em relacao a metrica g . Calculemos a variacao entao de

2
cada um dos termos de (B.1), d3 x gR, d3 x gR2 e d3 x gR .

(i) Variacao de R
Num sistema geodesico, R = , + , , onde a quantidade , i.e.,
o smbolo de Christoel, e denido como = 12 g [g, + g, g, ]. Por-

tanto, temos

1
R = + , = g = [ g + g g ] .
2

92
 Como consequencia, obtemos

1 
R = g ( ) g + g ( g g )
2 
+ g g . (B.2)

(ii) Variacao de R
Claramente,

R = R
1 
= g g + g ( ) g + 2g , (B.3)
2

onde ( ) 12 ( + ).

(iii) Variacao de R

Como R = (g R ), obtemos imediatamente

R = R g g + g 2g .

Agora, levando em conta que g = g g g , chegamos a conclusao de que

R = [R + g 2] g . (B.4)

(iv) Variacao ded3 x gR

Notando que g = 12 gg g , obtemos
   
3 g R
3
d xR g = dx g R
g + d3 x gg 2g
2


d3 x g g . (B.5)

93
 Contudo,
   
3 3
d x gg 2g =

d gg g ,

e d x g g =

d gg g .

Vamos supor que tanto g quanto g se anulem no entorno da regiao de

integracao. Assim, obtemos

  
3 g R
3
d x gR = dx g R
g . (B.6)
2

(v) Variacao de d3 x gR2

Utilizando (B.4) e (B.5), encontramos

 
3 2  g R2
3
d x gR = dx g 2RR 2 R
2

+ 2g 2R g .

(B.7)

2
(vi) Variacao de d3 x gR

Temos
  
3 3 2 g

2
d x gR R
= dx gR g + d3 x g2R R



d3 x g2R R g .

Por outro lado,

 
3  
d x g2R R =
d3 x g 2 R + g R + 2R g

Contudo, no sistema geodesico no qual estamos trabalhando [33],

1 1
R = R + R R R R , R = R
2 2

94
 Consequentemente,
 
3 
d x g2R R =
d3 x g R 2R R
g 
+ 2R R + 2R + 2R g . (B.8)
2

Com isso, nalmente obtemos

 
3 2 3  2 g
dx gR = d x g R R 2R R
2
g 
+ 2R + 2R g . (B.9)
2

Substituindo em (B.1) os valores acima calculados, levando em conta que


T
LM d3 x g g , (B.10)
2

encontramos as equacoes de campo para o modelo BHT como sendo

 
1 1 2 1 1 3 3
G + 2 R g R 2R R g 2R + 2R R2 g + RR
m2 2 4 4 16 4
2
= T , (B.11)
4

onde T e o tensor de energia-momento, e G R 12 g R e o tensor de Einstein.

B.2 Equacoes de campo linearizadas no gauge de

Teyssandier

Agora, de posse das equacoes de campo para o modelo BHT, podemos encontrar as
equacoes linearizadas. Desprezando os termos quadraticos em R e R , alem dos pro-

dutos RR e R R , a equacao (B.11) se torna

95
    
2 22 lin 1 2 2 1 1
+ 2 2 R + 2 2 Rlin R lin
= T , (B.12)
2 m2 2 2 m2 22 m22 2

onde utilizamos as seguintes quantidades linearizadas (em primeira ordem em )

lin  
R = 2 h , + , (B.13)
2 2

Rlin = 2h , e (B.14)
2
1
= h h. (B.15)
2

Ainda podemos reescrever (B.12) como

 
2 22 lin 1 lin 1 1
2
+ 2 2 R R + 2 2 ( 2 ) Rlin = T (B.16)
m2 2 2 m2 2

e, tirando o traco,

1 2 22 1 T
2
+ 2 2 Rlin + 2 2 2Rlin = ,
2 m2 m2 2

podemos reescrever (B.16), com o auxlio de (B.13), como

    
2 1 1 lin 1 2 1 lin
1+ 2 2h + R +
1 + 2 + R +
m2 2 4 2 m2 4m22
   
1 2 1 lin 1
+
1 + 2 + R = T T .

(B.17)
2 m2 4m22 4 2

Se zermos
 
2 1
1+ 2 + Rlin , (B.18)
m2 4m22

96
a equacao (B.17) ca numa forma mais compacta
    
2 1 1 lin 1 1
1+ 2 2h + R + ( + ) =

T T .

m2 2 4 2 4 2
(B.19)
A grande vantagem de se reescrever a equacao (B.17) na forma apresentada em

(B.19) e justamente o fato de que e sempre possvel encontrar um sistema de coorde-

nadas no qual se anule, na metrica linearizada [27]. Assim, a equacao de campo

linearizada pode ser reescrita como

 
  1
2 + m22 = T T , (B.20)
2 2

com
 
1 1 lin
2 2h + R . (B.21)
m2 2
Ou ainda, utilizando a denicao acima,
 
  1 1
2+ m22 = 2 2h + Rlin = T T . (B.22)
2 2 2

Note que e possvel eliminar o termo em Rlin em (B.22) se tomarmos o traco de (B.19)
e o somarmos a derivada de (B.18). Rearranjando os termos e utilizando a denicao
1
(B.14), encontramos 4 Rlin = 8 T . Substituindo em (B.22), camos com
 
1
2 2h T = T T ,
4 2 2

ou ainda

2 (h ) = (T T ) . (B.23)
2

97
 Contudo, a solucao das equacoes de Einstein no gauge de De Donder e dada por

2h(E) (E)
= (T T ) , com = 0, (B.24)
2

(E) (E) (E)

onde h 12 e, portanto, a solucao geral das equacoes de campo, sera dada
pela combinacao dos dois campos

h = h(E)
. (B.25)

B.3 Solucao das equacoes de campo para uma par-

tcula puntual

De posse das equacoes de campo (B.21) e (B.24), podemos encontrar uma solucao geral
para h , dado por (B.25). Para tanto, consideremos uma partcula puntual de massa

M em r = 0, descrita por

T = 0 0 M 2 (0)e (B.26)

T = T = 0 0 M 2 (0) = M 2 (r). (B.27)

No caso estatico, temos que (B.21) sera

 
2 1
( + m22 ) = 00 T T00
2 2

= (2 + m22 ) (r) = M 2 (r) (B.28)
4

que, no espaco dos momenta, se torna

 +  +
d2 k 2 d2 k ikr
2
(k + m22 ) (k)eikr = M e , (B.29)
(2) 4 (2)2

98
nos dando
M
ij (k) = 00 (k) = . (B.30)
4 (k2 + m22 )
Resolvendo para 00 (r), por exemplo,
 +  + 2 
d2 k dk M
00 (r) = (k)e
2 00
ikr
= 2
2 + m2 )
eikr
(2) (2) 4 (k 2
  2
1 + M
= dkk 2 deikr cos
(2)2 4 0 k + m22 0

1 + M M
= dkk 2 J (kr) =
2 0
K0 (rm2 ).
(2) 4 0 k + m2 8

Assim, obtemos nalmente

M
00 (r) = ii (r) = K0 (rm2 ). (B.31)
8

Resolvendo a equacao (B.24) tambem no caso estatico,

h(E)
= (T T )
2
M
h(E)
= 2 (0 0 ) . (B.32)
2

(E) (E)
Notamos que a componente h00 e nula, enquanto que, para hii , temos
 +  +
2 (E) M
d kk2 hii eikr = d2 keikr , (B.33)
2

e entao
(E) M
hii (k) = . (B.34)
2k2
Resolvendo a integral acima do mesmo modo que a anterior, encontramos

(E) M
hii = lim K0 (r), (B.35)
4 0

99
 1
que, nos da, para  ,
r0  
(E) M r
hii = ln , (B.36)
4 r0
onde r0 e um regulador.

Podemos escrever, nalmente, a solucao geral de uma partcula puntual de massa

M localizada em r = 0 para h como sendo

M
h00 = K0 (rm2 ) (B.37)
8   
M r
hii = 2 ln + K0 (rm2 ) . (B.38)
8 r0

100
Apendice C
Um resultado util para checar a
unitariedade de um modelo generico de
gravitacao em 3D

Teorema 1. Se m e a massa de uma partcula fsica generica de spin 2 relacionada a

um dado modelo gravitacional e k e a troca de momentum correspondente, entao

2 1
(T T 2 )|k2 =m2 > 0 and (T 2
T 2 )|k2 =0 = 0.
2
Aqui T (= T ) e a corrente externa conservada.

Notemos que o conjunto de vetores independentes no espaco dos momenta, k

(k 0 , k), k (k 0 , k),  (0, ), onde e um vetor unitario ortogonal a k, e uma

base adequada para expandirmos quaisquer vetores tridimensionais V (k). Utilizando

esta base, podemos escrever o tensor de corrente simetrico da seguinte maneira

T = Ak k + B k k + C  + Dk ( k ) + Ek ( ) + F k ( ) ,

onde a( b) 21 (a b + b a ).
A conservacao das correntes impoe as seguintes restricoes sobre os coecientes A,

101
B, D, E e F :

D 2
Ak 2 + (k + k2 ) = 0 (C.1)
2 0
D
B(k02 + k2 ) + k 2 = 0 (C.2)
2
Ek 2 + F (k02 + k2 ) = 0 (C.3)

Das Eqs. (B1) e (B2), obtemos Ak 4 + B(k02 + k2 )2 , enquanto a Eq. (B3) implica
que E 2 > F 2 . Por outro lado, saturando os ndices de T com momenta k , chegamos

a uma relacao de consistencia para os coecientes A, B e D:

Ak 4 + B(k02 + k2 )2 + Dk 2 (k02 + k2 ) = 0. (C.4)

Apos algumas manipulacoes algebricas utilizando as equacoes anteriores, encontra-

mos

 2
2 1 k 2 (A B) C k2
T T2 = + (E 2 F 2 ),
2 2 2 2
 
2 2 1
T T2 = k 2 2
(E F ) 2C(A B) . (C.5)
2

Portanto,
2 1
(T T 2 )|k2 =m2 > 0 e (T
2
T 2 )|k2 =0 = 0.
2

102
 Parte II

Determinacao de Limites para a

Massa do Foton

103
 As unicas partculas livres conhecidas no momento cujas massas de repouso podem
ser consideradas exatamente iguais a zero sao os fotons e os gravitons; os gluons, ou seja,

as partculas de gauge da QCD, que se supoe terem massa nula, nao sao observaveis
(fenomeno do connamento). Apesar de que a atual expansao acelerada do universo

poderia ser parcialmente atribuda a um efeito tipo massa do graviton, existem debates
acalorados na literatura sobre a possvel existencia de uma massa para esta partcula.

Este nao e caso, porem, no que tange a massa do foton. De fato, a procura por um
valor nao nulo para a massa de repouso do foton tem crescido consideravelmente ao
longo do tempo.

Por outro lado, como e bem conhecido, um foton massivo pode ser acomodado
de modo natural dentro do contexto da QED massiva, que nada mais e que uma

simples modicacao da QED padrao. Alem disso, a QED massiva e invariante sob as
transformacoes da relatividade especial. Como a QED massiva pertence a classe dos

assim chamados modelos que conservam a corrente [1], ela e renormalizavel [2]. Na
realidade a eletrodinamica massiva e, num certo sentido, mais simples teoricamente
que a QED padrao. Consequentemente, ela e uma ferramenta bastante adequada no

trato da questao concernente a uma possvel existencia de uma massa de repouso nao
nula para o foton.

A possvel existencia na natureza de um foton massivo (que suporemos ser des-

crito pela QED massiva) acarretaria alteracoes nos valores de resultados fsicos bem

conhecidos tais como velocidade da luz, campos eletromagneticos estaticos, e assim

por diante. Muitas estimativas tipo ordem de grandeza para a massa do foton
baseadas nestas possveis modicacoes dos valores de alguns resultados fsicos padroes

foram feitas nas ultimas decadas [35]. E importante notar que estas estimativas foram
obtidas empregando-se tanto metodos terrestres quanto extraterrestres.

104
 No captulo que se segue, iremos determinar limites superiores para a massa do
foton por meio de novos enfoques conceituais. Para tanto utilizaremos duas notaveis

predicoes provenientes ,respectivamente, da relatividade geral e da QED, a saber, de-

exao gravitacional da luz e momento magnetico anomalo do eletron, os quais foram

conrmados com grande precisao pelos experimentais. Os valores experimentais mais

recentes concernentes a estas quantidades serao tomados como entrada de modo a nos

fornecer os limites que estamos procurando.

105
Referencias Bibliogracas

[1] A. Ignatiev and G. Joshi, Phys. Rev. D 53 (1996) 984.

[2] D. Boulware, Ann. Phys. 56 (1970) 140.

[3] A. Goldhaber and M. Nieto, Rev. Mod. Phys. 43 (1971) 277.

[4] Liang-Cheng Tu, J. Luo, and G. Gilles, Rep. Prog. Phys. 68 (2005) 77.

[5] A. Goldhaber and M. Nieto, Rev. Mod. Phys. 82 (2010) 939.

106
Captulo 3
A Massa do Foton e a Interferometria de
Linha de Base Extensa

E encontrada uma relacao entre a massa do foton, sua frequencia () e o parametro

de deexao (), determinado pelos experimentais (o qual caracteriza a contribuicao

da curvatura do espaco para a deexao gravitacional). Este surpreendente resultado

nos permite concluir que o conhecimento dos parametros e e tudo que precisamos
para estabelecer limites gravitacionais sobre a massa do foton. Utilizando como

dados de entrada as mais recentes medidas de deexao de ondas de radio pelo Sol
(obtidas via VLBI), sao estimados limites superiores para a massa do foton.

3.1 Motivacao

Com o avanco das tecnicas de interferometria que ocorreu nas ultimas decadas, uma
nova luz foi lancada sobre os fenomenos fsicos nos quais efeitos quanticos e gravita-

cionais estao interligados. Consideremos, por exemplo, o efeito quantico da mudanca

de fase de neutrons causado pela sua interacao com o campo gravitacional da Terra, o

qual foi observado nos experimentos de Colella-Overhauser-Werner (COW) [1] e Bonse-

VLBI (Very Long Baseline Interferometry)

107
Wrobleski [2], realizados com interferometria de neutrons. Estes testes pioneiros tem
ajudado a nos convencer de que existem certos experimentos cujos resultados depen-

dem tanto da constante gravitacional quanto da constante de Planck. O fato de que

a mecanica quantica e a gravidade nao desempenham simultaneamente um papel im-

portante na maioria dos experimentos terrestres certamente e responsavel pela ideia

errada de que ambas nao podem estar intimamente ligadas em certas circunstancias

especiais.
Nas ultimas quatro decadas os experimentos tipo COW se tornaram mais sostica-
dos. O ultimo experimento com interferometria de neutrons [3] relata uma discrepancia

estatisticamente signicativa entre os valores teoricos preditos e os valores experimen-

talmente medidos para a mudanca de fase do neutron devido a gravidade. Supondo a

igualdade da massa inercial e gravitacional para o neutron, os experimentos mostraram

que o fator de fase do neutron e 1% menor do que o previsto [3], o que indica claramente

uma possvel violacao do princpio de equivalencia classico. A primeira vista, parece que
estes ultimos experimentos de interferometria estao em conito com os mais precisos
testes do princpio de equivalencia classico, tanto os realizados por meio de interfero-

metria atomica quanto os baseados em pendulos de torcao. Adunas, Rodriguez-Milla

e Ahluwalia [4, 5] mostraram, contudo, que cada um dos experimentos anteriormente

citados exploram diferentes aspectos da gravidade; e as tecnicas experimentais atuais,

quando associadas aos dados de neutrinos solares, podem ser capazes de explorar vio-

lacoes do princpio de equivalencia classico devidas a mecanica quantica. Eles tambem

predizem uma violacao quantica do princpio de equivalencia classico para a proxima
geracao de experimentos de interferometria atomica. Na verdade, do ponto de vista

operacional, nao se pode dizer que, em princpio, para certos sistemas quanticos, exista
uma igualdade das massas inercial e gravitacional. Portanto, chegamos a conclusao

108
de que a mecanica quantica e o princpio de equivalencia classico nao podem coexistir
pacicamente [6, 7].

Recentemente foi realizado um experimento interessante, nao diretamente relaci-

onado as tecnicas de interferometria, mas que diz respeito ao fato de que neutrons

podem existir em estados quanticos estacionarios formados no campo gravitacional

terrestre. Neste experimento, o estado quantico fundamental de neutrons no campo

gravitacional terrestre foi identicado na transmissao vertical de neutrons entre um

espelho horizontal (embaixo) e um detector/espalhador (em cima) [8, 9]. Este resul-
tado motivou Ernest a fazer uma pesquisa cuidadosa sobre auto-estados gravitacionais

em campos gravitacionais fracos [10, 11]. Na verdade, apesar de todos os estudos da

teoria quantica relacionados a estados atomicos e moleculares, muito pouco tem sido

feito para investigar as propriedades dos hipoteticos estados estacionarios que deveriam
existir em pocos de potenciais gravitacionais centrais similares, em particular, aqueles

com grandes numeros quanticos.

Os experimentos anteriores nos mostram claramente que fenomenos fsicos nos quais
a gravidade (o campo gravitacional terrestre, especicamente) e efeitos quanticos estao

mesclados nao estao mais alem do nosso alcance. Enfatizamos que em todas estas
investigacoes o campo gravitacional da Terra e descrito pela gravitacao newtoniana.

Por outro lado, desde sua publicacao original em 1915, a teoria da relatividade
geral de Einstein continua sendo uma area ativa de pesquisas tanto teoricas quanto

experimentais. Atualmente a teoria tem obtido sucesso em interpretar todos os dados

reunidos ate entao [12]. Alem disso, entre os chamados testes classicos da Relativi-
dade Geral existe um, a deexao gravitacional da luz, que tem sido conrmado com

uma exatidao cada vez maior. Na verdade, espera-se que uma serie de experimentos
melhorados utilizando a VLBI aumentara a precisao atual do parametro de deexao

109
em, pelo menos, um fator de 4 [13]. Alias, o valor atual para e de 0.99980.0003
(68% de nvel de conanca) [13], em acordo com a relatividade geral. Alem disso, e

de conhecimento geral que a deexao gravitacional da luz pelo Sol pode ser medida de
modo mais preciso utilizando-se comprimentos de onda na faixa de radio (por meio de

tecnicas de interferometria) em vez de ondas na frequencia do espectro visvel (com as

tecnicas opticas disponveis) [1417].

De fato, atualmente a VLBI e a tecnica mais precisa que temos a nossa disposicao
para medidas de deexao gravitacional de ondas de radio [13]. Agora, levando em conta
que a busca por limites superiores para a massa do foton tem aumentado ao longo das

ultimas decadas [2022], seria interessante estimar limites gravitacionais para a massa
do foton considerando as medidas mais recentes da deexao gravitacional de ondas de

radio obtidas por meio da VLBI. Este e, precisamente, o objetivo deste captulo.
O captulo e organizado como se segue. Na Secao 3.2 encontramos o angulo de

deexao sofrido por um foton massivo no presente contexto enquanto, na Secao 3.3,
mostramos que existe um vnculo entre a massa do foton, sua frequencia e o parametro
de deexao determinado pelos experimentais (o qual caracteriza a contribuicao da cur-

vatura do espaco para a deexao gravitacional). Partindo deste interessante resultado,

limites superiores para a massa do foton sao encontrados na Secao 3.4, considerando

como dados as medidas mais recentes da deexao gravitacional solar de ondas de ra-
dio obtidas via VLBI [13]. Para concluir, discutimos na Secao 3.5 se os limites que

Em todas estas pesquisas, o foton e descrito pela QED massiva, a qual nada mais e do que a
extensao mais simples da QED padrao. Sua lagrangiana pode ser escrita como
1 2 1
L = F + m2 A2 J A , (3.1)
4 2
onde F (= A A ) e o tensor de intensidade do campo, J e a corrente (eletrica) conservada
e m e a massa do foton. Aqui, ndices sao levantados e abaixados com e , respectivamente.
Na verdade, a QED massiva e teoricamente mais simples do que a teoria padrao [18], alem de ser
renormalizavel [19].

110
estimamos podem ou nao ser melhorados.
Nas nossas convencoes  = c = 1, e a assinatura da metrica e (+ - - -).

3.2 Combinando Relatividade Geral, QED Massiva

e VLBI para limitar gravitacionalmente a massa
do foton

Iniciamos lembrando que a lagrangiana para o campo do foton massivo acoplado mi-
nimamente a gravitacao e dada por


1 m2
L= g g g F F + g A A . (3.2)
4 2

Por outro lado, para pequenas utuacoes em torno da metrica de Minkowski , a

metrica completa pode ser escrita como

g = + h , (3.3)

onde 2 = 32G. Aqui G e a constante de Newton. De modo analogo ao que zemos

no Captulo 2, temos tambem

g = h + 2 h h + ... . (3.4)

e
2 2
g = 1 + h + h2 h h + ... . (3.5)
2 8 4

111
3.2.1 Calculo da lagrangiana de Interacao

A lagrangiana de interacao e encontrada substituindo-se (3.3), (3.4) e (3.5) em (3.2).

Calculando separadamente cada termo de (3.2) via as mencionadas substituicoes, ob-

temos

14 F F

1 1 1
F F = F g g F = F ( h )( h )F
4 4 4
1
= F ( h h )F
4
1 1
= F F + (F h F + F h F ) = F F + (2F F )h
4 4 4 4
1
= F F + h F F . (3.6)
4 2

m2
2
A A

m2 m2 m2
A A = A g A = A ( h )A
2 2 2
m2 m2
= A A A A h . (3.7)
2 2

Consequentemente,

 1 m2 m2  
L = F F + h F F + A A h A A 1 + . (3.8)
4 2 2 2 2

Supondo entao que o foton massivo e espalhado por um campo gravitacional externo
fraco (veja a Figura 3.1), chegamos a lagrangiana para a interacao em questao, em

primeira ordem em ,

112
 
1 m2 2
Lint = hext F F F F +

A A m A A . (3.9)
2 4 2

Note que agora os ndices sao levantados (abaixados) com ( ).

3.2.2 Calculo da funcao de vertice

Partindo da lagrangiana (3.9) e expressando os campos A como A eip

x

 , po-
demos calcular a funcao de vertice da interacao descrita pela Figura 3.1. Calculando
separadamente os termos massivos e os nao-massivos, resulta

 
m2 2 m2 2
2 hext A A m A A = 2 hext A A m A A

2 2 2 2
 
m2 h ext   2 

2
2
= m hext [ 2 ]   ;
2

p

Figura 3.1: Graco de Feynman para a interacao entre um foton massivo e um campo
gravitacional externo.

113
    
2 h
ext F F = h
ext ( A A )( A A )
2
    

hext p p p p p p + p p   ;


1 1
2 hext F F

= h [( A A )( A A )]
2 4 4 ext
1  
hext p p p p p p + p p  
4
1  
= hext p p p p   .
2

Juntando os termos acima, encontramos a funcao de vertice desejada no espaco dos

momenta,



V (p, p ) = hext (k) (m2 p p )( 2 ) + p p
2 
+ 2(p p p p + p p ) . (3.10)

Ainda nos resta determinar a forma de h

ext (k). Supondo que o campo gravitacional

externo fraco e gerado por uma massa puntual M em r=0, encontramos imediatamente

a expressao para o campo gravitacional externo resolvendo as equacoes de Einstein

linearizadas no gauge de De Donder [23]. A expressao resultante e

M
h
ext (r) = ( 2 0 0 ). (3.11)
16r
Portanto, o campo gravitacional no espaco dos momenta, h
ext (k), e dado por

114
 
h
ext (k) = d3 reikr h
ext (r)

M  
= .
0 0
(3.12)
2k2 2

3.2.3 Calculo da secao de choque diferencial nao-polarizada

Agora estamos prontos para calcular a secao de choque diferencial nao-polarizada para
o processo exibido na Fig. 3.1. Para tanto, relembramos que a expressao para a citada

secao de choque e

1 1  2
3 3
d
= M , (3.13)
d (4)2 3 r=1 r =1 rr
com

Mrr = r (p)r (p )V (p, p ), (3.14)

onde r (p) e r (p ) sao os vetores de polarizacao para os bosons vetoriais incidentes e

espalhados, respectivamente. Estes vetores, por sua vez, satisfazem a relacao


3
p p
r (p)r (p) = + . (3.15)
r=1
m2
A equacao (3.13) para a secao de choque se torna entao

1 1 
3 3
d
=  (p)r (p )r (p)r (p )V (p, p )V (p, p )
d (4)2 3 r=1 r =1 r

1 1 2 p p p p p p p p
V V V V V + V V .(3.16)
(4)2 3 m2 m2 m4

115
 A conservacao do momentum e dada por k+p = p e como nao existe troca de
energia com o campo gravitacional externo, p0 = p0 = E, o que nos da tambem |p| =

|p |. Novamente, vamos trabalhar cada termo da expressao acima individualmente a

m de encontrar uma relacao em termos dos momenta envolvidos.

p p p p 1  2
V V = (m p p )(20 0 ) p p + m2
m4 m4

+2(E0 p + E0 p E 2 ) p p

= 4E 4 ;

p p p p
V V
= V V
m2 m2
 2
V p
=
m
 2 2
2m E0
= = 4m2 E 2 ;
m

 
V V = 2p p (p p ) 4E 2 2m2 + 5m4 + 4E 2 m2 + 8E 4 .

Substituindo as quantidades anteriores na equacao (3.16), obtemos a secao de cho-

que diferencial nao-polarizada

 2  
d 2 M 1 2 3
= 3p4 + m4 +2p2 m2 +2p2 (p2 +2m2 ) cos +p4 cos2 , (3.17)
d 4k2 (4)2 3 2

116
onde e o angulo de espalhamento. Como k + p = p e |p| = |p |, temos

k2 = (p p)2

= p2 + p2 2|p ||p| cos

= 2|p|2 (1 cos ).

Utilizando a relacao acima e lembrando que 2 = 32G, (3.17) se torna

 2  
d MG 2 3
= 3p4 + m4 + 2p2 m2 + 2p2 (p2 + 2m2 ) cos + p4 cos2 .
d p2 (1 cos ) 3 2

Como p2 = E 2 m2 , para angulos pequenos camos com

 
d M 2 G2 1  4 
= 2 2
4E 4m2 E 2 + m4
d (2 sin 2 ) (E m )
2 2 2

   2
4 m2
M 2 G2 4E 1 2E 2
= 4  2 2
4 16 E4 1 m E 2
 2
2 2 1 m2
16M G 2E 2
= 2 , (3.18)
4 1 m E 2

onde E e a energia do foton espalhado.

A secao de choque diferencial acima pode ser relacionada a uma trajetoria classica

com um parametro de impacto b atraves da relacao

d
bdb = d. (3.19)
d

Das Eqs. (3.18) e (3.19), chegamos a conclusao de que

117
  
m2
4M G 1 2E 2
= 2 , (3.20)
b 1 mE2

a qual, no limite ultrarelativstico, i.e., E  m, leva a

 
4M G 1 m2
= 1+ ,
b 2 E2
ou seja,

 
m2
= E 1 + (3.21)
2E 2
 
m2
= E 1 + 2 2 , (3.22)
8

4M G
onde e a frequencia do foton massivo e E b
.

3.3 Encontrando uma relacao entre massa, frequen-

cia e parametro de deexao para o foton mas-
sivo

Antes de continuarmos, seguem alguns comentarios pertinentes.

Recentemente foi mostrado que as secoes de choque diferenciais nao-polarizadas

para o espalhamento de diferentes partculas quanticas sao dependentes do spin, o
que esta em desacordo com o princpio da equivalencia classico [7] (veja a Tabela

1). Este resultado levanta uma importante questao: Por que o campo gravitaci-
onal percebe o spin? Porque temos a presenca de uma troca de momentum (k)

nao nula no espalhamento, que e sensvel a estrutura interna (spin) da partcula.

118
 Nao obstante, se escolhermos duas expressoes quaisquer daquelas listadas na Ta-
bela 1, notamos que a diferenca entre elas e sempre extremamente pequena para

angulos de deexao tpicos. Para demonstrar isto para partculas sem massa, por
exemplo, estudamos o comportamento de

d
( d ) d
( d )s ( d
d
)s=0
d
d
, (3.23)
( d )s=0 ( d )s=0

como uma funcao do angulo de espalhamento (veja a Fig. 2). E trivial mostrar
que para pequenos angulos a expressao anterior se reduz a

d
( d ) s2
d
. (3.24)
( d )s=0 2

d
( d )
Para um angulo de deexao tpico, digamos 106 , encontramos d
( d )s=0

1012 . A deteccao de um efeito tao pequeno esta, claramente, alem da tecnologia

atual. Consequentemente, para estes minusculos angulos de deexao, as secoes

de choque nao serao afetadas pelo spin das partculas.

De modo a recuperar os resultados geometricos de Einstein a partir da Tabela 1,

devemos ter k 0; em outras palavras, no limite nao-trivial de pequenas trocas

de momentum, o qual corresponde a um limite de angulo pequeno nao-trivial

ja que |k| = 2 |p| sin 2 , as partculas massivas (nao massivas) se comportam da
mesma maneira, independentemente do spin. De fato, quando o spin e desli-

gado, i.e, para pequenos angulos, obtemos da Tabela 1 que para m = 0,

119
 d
m s d

 2
GM
0 0 sin2 2

 2  2
= 0 0 GM
sin2 2
1+
2

 2
1
0 2
GM
sin2 2
cos2
2

 2   
1
= 0 2
GM
sin2 2
cos2 2 +
4
1 + + 3 cos2 2

 2
0 1 GM
sin2 2
cos4
2

 2   
1
= 0 1 GM
sin2 2 3
+ 23 cos4 2
3
1 3
4
4 cos2
2

 2  
0 2 GM
sin2 2
sin8 2 + cos8
2

Tabela 3.1: Secoes de choque diferenciais nao-polarizadas para o espalhamento de diferentes part-
culas quanticas por um campo gravitacional fraco externo gerado por uma partcula puntual estatica
2
1v2
de massa M . Aqui m e a massa da partcula, s o spin, o angulo de espalhamento e m
p2 = v 2 ,
com v e p sendo a velocidade e o 3-momento, nesta ordem, da partcula incidente.

d 16G2 M 2
, (3.25)
d 4

enquanto para m = 0,

 2
d 16G2 M 2
1+ . (3.26)
d 4 2

Utilizando a Eq. (3.18) conclumos que para m = 0,

120
 4GM
, (3.27)
b

enquanto para m = 0,

 
4GM
1+ . (3.28)
b 2

A primeira equacao nos da o angulo de deexao gravitacional para uma partcula

sem massa um resultado obtido por Einstein ha muito tempo; enquanto a

segunda coincide com a previsao da relatividade geral para a deexao de uma

partcula-teste massiva classica por um campo gravitacional externo fraco [24].
Vale notar que as Eqs. (3.21) e (3.28) sao exatamente a mesma. Resumidamente,

para angulos pequenos os resultados da Tabela 1 nao apenas se reduzem aqueles

previstos pela teoria geometrica de Einstein mas eles tambem concordam com o

princpio de equivalencia classico.

A primeira vista parece que a Eq. (3.22) prediz um angulo de deexao dispersivo
para fotons massivos. Na verdade, esta e uma falsa impressao; ainda mais, e

simples mostrar que a Eq. (3.22) pode ser reescrita como

 
3 v2
= E . (3.29)
2

Apos estas importantes digressoes, retornemos a analise da Eq. (3.22). O primeiro

termo na expressao (3.22) coincide com aquele obtido por Einstein em 1916, resolvendo

121
 d
( d )
Figura 3.2: d
( d )s=0
como funcao do angulo de espalhamento .

a equacao da propagacao da luz no campo gravitacional de um corpo estatico [25],

enquanto o segundo termo e a correcao mais importante para a massa m do foton

massivo.
Por outro lado, o angulo do desvio gravitacional determinado por grupos experi-

mentais e geralmente expresso atraves da relacao [26]

1+
exp = E , (3.30)
2

onde e o parametro de deexao, precisa e inequivocamente determinado pelos expe-

rimentais medindo a deexao de radiacao eletromagnetica pelo campo do Sol.

Das Eqs. (3.22) e (3.30), obtemos

'
m < 2 |1 |. (3.31)

122
 Este surpreendente resultado claramente mostra que existe um vnculo entre a massa
do foton e os parametros e . Alem disso, nos diz que o conhecimento destes para-

metros e tudo que precisamos saber para estabelecer limites superiores para a massa
do foton [27, 28].

3.4 Encontrando limites gravitacionais superiores

para a massa do foton

Agora estamos prontos para encontrar limites gravitacionais superiores para a massa

do foton. Para tanto, faremos uso de medidas recentes da deexao gravitacional solar
de ondas de radio encontradas por Fomalont et al. [13] utilizando a VLBI. As solucoes

para obtidas por meio destas medidas, assim como os limites correspondentes para a
massa do foton que estimamos utilizando a Eq. (3.31), estao organizadas na Tabela 2.
De modo a determinar o melhor valor de ( 1) pelos resultados experimentais, os

autores citados anteriormente minimizaram a seguinte expressao

 2
1 Pd (i) 0.5( 1)Dd (i)
2k = , (3.32)
k d,i d (i)

onde Pd (i) e d (i) sao as medidas de deslocamento de posicao e o erro estimado, nesta

ordem. O termo Dd (i) e a media da predicao do desvio gravitacional previsto pela

relatividade geral, tomada em cada sessao de observacao. Os d s e 2k s diretamente
relacionados as solucoes para sao exibidos na Tabela 2.

Na Tabela 2, os valores relacionados a 43 GHz (sem efeitos de corona) e o mais

preciso, ja que possui o menor valor para 2 . Isto se deve ao atenuamento de alguns

efeitos da corona solar, e o aumento da estimativa dos erros na posicao da fonte. Os

dois valores para 43 GHz sem a correcao dos efeitos de corona mostram as alteracoes

123
 Tipo de Solucao ( 1) 104 d 104 2k m 1012 (M eV )
43 GHz data (sem efeitos de corona) -2.4 3.2 0.9 2.8
43 GHz -1.0 2.6 2.2 1.8
43 GHz - Oct05 -3.2 2.8 1.1 3.2
23 GHz - Oct05 -2.0 2.4 4.7 1.3

Tabela 3.2: Limites superiores para a massa do foton estimados utilizando as solucoes para
encontradas por Fomalont et al. [13].

nas medidas realizadas muito proximas ao Sol na secao de Oct05; a primeira destas
solucoes foi encontrada utilizando dados obtidos durante secoes que duraram varios

dias, enquanto a segunda e baseada em dados obtidos no intervalo de um unico dia de

observacao. Finalmente, os valores para 23 GHz sugere que a refracao devido a corona

(a qual e quatro vezes maior do que a medida em 43 GHz), domina a sensibilidade do

experimento em 23 GHz.

E importante notar que Fomalont et al. [13] utilizaram para os resultados que
obtiveram no artigo uma media das quatro solucoes exibidas na Tabela 2, obtendo =
0.99980.0003. Deste resultado e assumindo que um foton massivo passando proximo

a borda solar possui uma frequencia = 43GHz, o que e perfeitamente justicavel ja

que os dados deles sao, em sua maioria, obtidos em observacoes de 43GHz, onde os

efeitos de refracao da corona solar sao desprezveis alem de 3 graus do Sol nos obtemos
outro limite gravitacional para a massa do foton, a saber, m 2.5 1012 M eV .
Devido a razoes ja discutidas, chegamos a conclusao de que entre os limites gravita-

cionais para a massa do foton que encontramos, o mais conavel e m 2.81012 M eV .

Alem disso, o associado foi determinado, por um lado, utilizando uma unica frequen-

cia; por outro lado, ele possui o melhor valor para 2 . Ademais, os dados utilizados
para o calculo de vem principalmente das observacoes cujos efeitos de refracao da

corona solar sao desprezveis.

124
 Chamamos atencao ao fato de que a Eq. (3.20) pode tambem ser deduzida, a
la Einstein, procurando uma solucao aproximada para equacao de movimento da ge-

odesica de uma partcula teste massiva no campo de Schwarzschild. Adotando esta

abordagem, uma expressao para o angulo de deexao de uma partcula pelo Sol foi
 3
obtido ate a ordem de GM
b
na Ref. [24]. Este tipo de deducao, contudo, e muito
trabalhosa. Por outro lado, Golowich, Gribosky e Pal [29], inves da usual aborda-

gem geometrica, consideraram o fenomeno da curvatura da luz como um problema de

espalhamento quantico. Este tratamento, que nao so e instrutivo como tambem relati-
vamente simples quando o campo gravitacional e fraco, permitiu-lhes facilmente obter
GM
uma expressao para a deexao gravitacional de partculas massivas ate uma ordem b
.
Um resultado identico foi encontrado por Mohanty, Nieves e Pal [30], utilizando um

metodo desenvolvido originalmente por Ohanian e Runi [23].

Seguem alguns comentarios sobre a validade do resultado que obtivemos e como os

mesmos podem ser melhorados num futuro proximo.

De acordo com a relatividade geral, os fotons nao apenas sao deetidos mas
tambem sofrem um atraso devido a curvatura do espaco-tempo produzida por

qualquer massa. Alem disso, o desvio e o atraso da luz sao proporcionais a +

1. Consequentemente, tecnicas utilizadas para medir tal atraso podem tambem

ser empregadas para encontrar limites para a massa do foton. E interessante

notar que ha alguns anos, Bertotti, Iess e Tortora [15] relataram uma medida

do deslocamento de frequencia de fotons (na faixa de radio) trocados com a

espaconave Cassini quando eles passavam proximos ao Sol, levando a um resultado
para que concorda com as previsoes da relatividade geral com uma sensibilidade

que se aproxima a um nvel no qual, teoricamente, desvios sao esperados em

125
 alguns modelos cosmologicos [31, 32].

A Eq. (3.20) foi deduzida assumindo que o campo responsavel pela deexao

do foton e um campo gravitacional estatico. Contudo, como sabemos, tanto o

Sol quanto os planetas nao estao em repouso no Sistema Solar. Na verdade,

eles estao se movendo em relacao ao baricentro do Sistema Solar e o observador.

Este movimento certamente introduz correcoes dependentes da velocidade para a

equacao relativstica da deexao gravitacional da luz. Como consequencia, a cor-

recao a deexao gravitacional induzida pelos movimentos citados anteriormente
devem estar correlacionadas com a correcao da massa do foton exibida na Eq.

(3.22). Este fato nos leva a colocar uma questao importante: Atualmente, a
tecnologia existente e sensvel o suciente para detectar estes pequenos efeitos

da relatividade geral causados pela dependencia temporal do campo gravitaci-

onal? Kopeikin [33] arma que futuros experimentos de deexao de raios de

luz [34] tais como o BepiColombo (com ondas de radio)[35], ASTROD (com la-
ser) [36] e LATOR [37] irao, denitivamente, alcancar a precisao suciente para
medir PPN , PPN e PPN , e compara-los as correcoes pos-newtonianas para o

atraso e deexao causadas pelo movimento de corpos massivos no Sistema So-

lar [38]. Aqui, o desvio da relatividade geral e descrito atraves da comparacao

dos parametros PPN PPN PPN 1, PPN PPN 1 e PPN PPN 1.

Alternativamente, pode-se mostrar que, utilizando a equacao para o pos-PPN

(post-post-Newtonian) que dita o atraso temporal, t, a qual foi obtida por

Kopeikin atraves do acoplamento dos parametros PPN com os termos dependen-
tes da velocidade, para experimentos gravitacionais com a luz se propagando no

post-Newtonian parameters

126
 campo do Sol,
 
  r1 + r2 + r12
t 1 + ln , (3.33)
r1 + r2 r12

com

PPN 2 , (3.34)

onde  (= 5.3 108 ) e a velocidade do Sol (em unidades naturais) em relacao

ao baricentro do Sistema Solar, r12 e a coordenada da distancia entre os pontos

de emissao e observacao, r1 , r2 sao as distancias radiais dos pontos de emissao e

observacao, respectivamente. Agora, notando que as missoes espaciais LATOR e

ASTROD irao medir o parametro PPN com uma precisao aproximada de 109 [36,

37], chegamos a conclusao de que, num futuro proximo, correcoes ao caso estatico
do atraso de tempo explicitamente dependentes da velocidade devido ao campo
gravitacional solar devera, aparentemente, ser levado em conta. Respondamos

a questao levantada acima. Para simplicar vamos nos restringir a medidas de

desvio da luz pelo Sol obtidas por VLBI. Atualmente os grupos experimentais

determinaram os parametros PPN utilizando tecnicas de VLBI com uma precisao

de 104 [13]. Portanto, a citada correcao dependente da velocidade e muito

pequena e pode ser desprezada na determinacao de PPN . Na verdade, a deteccao

de um efeito tao pequeno esta alem da tecnologia atual.

Hoje em dia, como ja dissemos, a VLBI e a tecnica mais precisa que temos a nossa

disposicao para medir a deexao gravitacional de ondas de radio regularmente

[14, 16, 17].

Medidas da deexao da luz com tecnicas opticas podem vir a se tornarem mais

127
 vantajosas para determinar o parametro num futuro proximo [39].

3.5 Discussao

Os limites que encontramos para a massa do foton estao entre aqueles publicados pelo
Particle Data Group, mas certamente sao consideravelmente maiores do que outros

encontrados por la (m < 1 1018 eV ) [40]. Eles sao, todavia, comparaveis a outros
limites existentes [2022]. Vamos, entao, discutir se um limite melhor para a massa do
foton pode ser obtido utilizando a Eq. (3.31). Em primeiro lugar, se as deexoes medi-
'
das utilizando VLBI puderem ser feitas com uma precisao maior, o valor de |1 |
poderia ser reduzido, dando como resultado uma melhor estimativa gravitacional. De

acordo com Fomalont et al. [13], uma serie de experimentos projetados com a VLBI
poderiam aumentar a precisao de futuras medidas em pelo menos um fator de 4. Em

segundo lugar, se as medidas de deexao puderem ser obtidas em frequencias menores,

mantendo o valor do parametro de deexao xo, o limite gravitacional sera melhorado
em proporcao direta a frequencia. Este ponto, contudo, e muito delicado. Na verdade,

como ja mencionamos anteriormente, ate agora os melhores resultados obtidos para a

deexao gravitacional via VLBI sao aqueles advindos principalmente de frequencias em

43GHz, onde os efeitos da corona solar sao desprezveis alem de 3 graus da borda do
Sol. Contudo, a menor frequencia utilizada pelos radio-astronomos e de 2GHz. Infeliz-
mente, as medidas realizadas nesta frequencia sao menos conaveis devido aos efeitos

de refracao na corona, que e inversamente proporcional ao quadrado do comprimento

de onda utilizado. Na verdade, os radio-astronomos utilizam em seus experimentos

uma mistura de diferentes frequencias mas as contribuicoes mais signicativas vem, em

geral, de 43GHz. A possibilidade de melhorar os limites gravitacionais para a massa

128
do foton, neste caso, e bem limitada.
Enfatizamos que, ate entao, apenas duas tentativas de limitar a massa do foton

utilizando medidas de deexao gravitacional foram feitas [41, 42]. Infelizmente, estas
estimativas nao sao muito conaveis pois em ambas os valores do parametro de deexao

utilizados foram superestimados, enquanto os valores da frequencia foram subestima-

dos. Na verdade, as frequencias utilizadas nas estimativas citadas estao na vizinhanca

da menor frequencia medida pelos radio-astronomos ( 2GHz).

E importante notar que as estimativas gravitacionais sobre a massa do foton que
obtivemos, como a maioria das estimativas feitas com o objetivo de limitar a massa do

foton que estao disponveis na literatura, sao essencialmente argumentos de ordem de

magnitude. Todavia, nossos esforcos neste trabalho sao baseados em novas abordagens

conceituais ao assunto; alem disso, no calculo dos limites, utilizamos como dados de
entrada os melhores resultados disponveis atualmente.

Para concluir, observamos que recentemente encontramos um limite quantico para

a massa do foton (m 1.6 1010 M eV ) baseado no calculo do momento magnetico
anomalo do eletron no contexto da eletrodinamica de Proca [43, 44]. Dentre todas

as referencias sobre o assunto que consultamos, esta parece ser a primeira vez que
um limite quantico para a massa do foton e estimado . Vale notar que este limite

quantico e uma ordem de magnitude maior do que o limite obtido atraves da deexao
gravitacional.

Abaixo listamos, para efeito de comparacao, alguns valores de limites para a massa
do foton utilizando metodos de interferometria baseados na dispersao da luz.

Note que o limite para a massa do foton obtido por Boulware e Deser [45] via o efeito Aharonov-
Bohm (o qual esta presente na eletrodinamica massiva) e um limite semiclassico; nosso limite, ao
contrario daquele obtido por Boulware-Deser, se baseia em verdadeiros efeitos (de loop) quanticos.

129
 Autor (ano) Tipo de medida Limites para m (M eV )
Froome (1958) Interferometria de ondas de radio 2.4 1013
Warner et al. (1969) Observacoes do pulsar na Nebulosa do Caranguejo 2.9 1014
Bay et al. (1972) [46] Emissao de pulsar 1.7 1019
Brown et al. (1973) Pulsos de radiacao curtos 7.9 107
Schaefer (1999) Explosoes de raios gama(GRB980703) 2.4 1017
Explosoes de raios gama (GRB930229) 3.4 1012

Tabela 3.3: Alguns limites para massa do foton obtidos utilizando medidas de dispersao da veloci-
dade da luz em diferentes frequencias do espectro eletromagnetico (em ordem cronologica).

130
Referencias Bibliogracas

[1] R. Colella, A. Overhauser and S. Werner, Phys. Rev. Lett. 34 1472 (1975).

[2] U. Bonse and T. Wroblewski, Phys. Rev. Lett. 51 1401 (1983).

[3] K. Littrell, B. Allman and S. Werner, Phys. Rev. A 56 1767 (1997).

[4] G. Adunas, E. Rodriguez-Milla and D. V. Ahluwalia-Khalilova, Gen. Relativ. Gra-

vit. 33 183 (2001).

[5] G. Adunas, E. Rodriguez-Milla and D. V. Ahluwalia-Khalilova, Phys. Lett. B 485

215 (2000).

[6] A. Accioly and R. Paszko, Phys. Rev. D 78 064002 (2008).

[7] A. Accioly and R. Paszko, Adv. Studies Theor. Phys. 3 65 (2009).

[8] V. Nesvizhevsky et al., Nature 415 297 (2002).

[9] V. Nesvizhevsky et al., Phys. Rev. D 67 102002 (2003).

[10] A. Ernest, J. Phys. A: Math. Theor. 42 115207 (2009).

[11] A. Ernest, J. Phys. A: Math. Theor. 42 115208 (2009).

131
[12] S. Turyshev, Usp. Fiz. Nauk 179 3 (2009).

[13] E. Fomalont et al., Astrophys. J. 699 1395 (2009).

[14] D. Robertson, W.Carter and W. Dillenger, Nature 349 768 (1991).

[15] B. Bertotti, L. Iess, and P. Tortora, Nature 425 374 (2003).

[16] S. Shapiro et al., Phys. Rev. Lett. 92 121101 (2005).

[17] S. Lambert and C. Le Poncin-Latte, Astron. Astrophys. 499 331 (2009).

[18] A. Ignatiev and G. Joshi, Phys. Rev. D 53 984 (1996).

[19] D. Boulware, Ann. Phys. 56 140 (1970).

[20] A. Goldhaber and M. Nieto, Rev. Mod. Phys. 43 277 (1971).

[21] Liang-Cheng Tu, J. Luo, and G. Gilles, Rep. Prog. Phys. 68 77 (2005).

[22] A. Goldhaber and M. Nieto, Rev. Mod. Phys. 82 939 (2010).

[23] H. Ohanian and R. Runi, Gravitation and Spacetime (Norton, New York, 1994),

2nd ed.

[24] A. Accioly and S. Ragusa, Class. Quant. Grav. 19 5429 (2002) [Erratum Class.

Quant. Grav. 20 4963 (2003)].

[25] A. Einstein, Ann. Phys. 49 769 (1916).

[26] C. Will, Theory and Experiment in Gravitational Physics (Cambridge: Cambridge

Univ. Press, 1993).

132
[27] A. Accioly, J. Helayel-Neto and E. Scatena, Phys. Lett. A 374 3806 (2010). This
essay was awarded an honorable mention in the 2010 Essay Competition of the

Gravity Research Foundation.

[28] A. Accioly, J. Helayel-Neto and E. Scatena, Int. J. Mod. Phys. D 19 2393 (2010).

[29] E. Golowich, P. Gribosky, and P. Pal, Am. J. Phys. 58(7), 688 (1990).

[30] S. Mohanty, J. Nieves, and P. Pal, Phys. Rev. D 58, 093007 (1998).

[31] T. Damour and A.M. Polyakov, Nucl. Phys. B 423, 532 (1994).

[32] T. Damour, F. Piazza, and G. Veneziano, Phys. Rev. D 66, 046007 (2002).

[33] S. Kopeikin, Mon. Not. R. Astron. Soc. 399, 1539 (2009).

[34] S. Kopeikin and B. Mashhoon, Phys. Rev. D 65, 064025 (2002).

[35] A. Milani et al., Phys. Rev. D 66, 082001 (2002).

[36] W.-T. Ni, Nucl. Phys. B, Proc. Suppl. 166, 153 (2007).

[37] S. Turyshev, M. Shao and K. Nordtvedt, Class. Quantum Grav. 21, 2773 (2004).

[38] J. Plowman and R. Hellings, Class. Quantum Grav. 23, 309 (2006).

[39] S. Klioner, P. Seildemann, and M. Soel, International Astronomical Union Sym-

posium (Cambridge University Press, Cambridge, UK, 2010), Vol. 261.

[40] J. Beringer et al. (Particle Data Group), Phys. Rev. D 86 010001 (2012).

[41] D. Lowenthal, Phys. Rev. D 8 2349 (1973).

[42] A. Accioly and R. Paszko, Phys. Rev. D 69 107501 (2004).

133
[43] A. Accioly, J. Helayel-Neto and E. Scatena, Phys. Rev. D 82 065026 (2010).

[44] A. Accioly, J. Helayel-Neto, R. Turcati, J. Morais and E. Scatena, Classical Quan-

tum Gravity 27 205010 (2010).

[45] D. Boulware and S. Deser, Phys. Rev. Lett. 63 2319 (1989).

[46] Z. Bay and J. White, Phys. Rev. D 5, 796 (1972).

134
Eplogo

Como mencionamos na introducao, este trabalho tinha dois propositos especiais:

Tentar compreender melhor a conjectura em voga sobre a incompatibilidade entre

unitariedade e renormalizabilidade em teorias de ordem mais alta.

Encontrar limites para a massa do foton que na fossem baseados somente em

calculos de ordem de grandeza, mas sim em resultados provenientes de calculos

rigorosos realizados no contexto da QFT e onde fossem utilizados como dados de
entrada as mais recentes medidas experimentais.

Sera que os objetivos foram atingidos? Antes de responder a esta indagacao e necessario
analisar o que os resultados obtidos ate aqui nos dizem sobre os dois topicos previamente

mencionados.
Iniciemos pela conjectura sobre as teorias com derivadas de ordem superior.

No caso especco da versao atual da eletrodinamica de Lee e Wick foi mostrado

em detalhe que esta teria e nao-unitaria, porem renormalizavel. Por sua vez, Lee e
Wick tentaram contornar o problema da nao unitariedade introduzindo modicacoes ad

hoc nas regras usuais da teoria quantica de campos. Entre estas modicacoes engenho-
sas, porem inortodoxas, a assim chamada prescricao de Lee e Wick para a continuacao

135
analtica das amplitudes foi alvo de crticas acerbas [1, 2]. Apesar de altamente insti-
gantes e bem formuladas, as prescricoes de Lee e Wick nao foram ainda incorporadas

ao ferramental da fsica teorica. Assim sendo, e mais prudente, por ora, nos basearmos
na teoria quantica de campos ortodoxa que nos assevera que a eletrodinamica de Lee e

Wick e nao-unitaria. E importante notar, por outro lado, que nossos calculos mostra-
ram claramente que os efeitos devidos aos termos de ordem superior so se manifestam

em distancias muito pequenas, como se esperava. Portanto, fora desta diminuta regiao,
nao ha diferenca sensvel entre as eletrodinamicas de Maxwell e a de Lee e Wick. Em
resumo, dentro do esquema padrao da teoria de campos, a eletrodinamica de Lee e

Wick e renormalizavel mas nao-unitaria.

No que tange a gravitacao massiva em 3D, vericamos que ela e a unica teoria de

ordem superior cuja versao linearizada e unitaria. Sera ela tambem renormalizavel?
Se assim for, a conjectura sobre a incompatibilidade entre unitariedade e renormali-

zabilidade em teorias de ordem superior cai por terra. Esta questao, como todas que
envolvem renormalizacao de terias gravitacionais, e sempre polemica. Em 2009 Oda
[3] mostrou que esta teoria era renormalizavel e unitaria. Bastante recentemente,

no entanto, Muneyuki e Ohta [4, 5] demonstraram que as teorias de gravitacao com

derivadas de ordem superior em 3D sao renormalizaveis, porem aquelas teorias que

demandam uma relacao especial entre seus coecientes, incluindo a nova gravitacao
massiva, nao o sao. Consequentemente, a gravitacao massiva em 3D e unitaria mas

nao-renormalizavel .
Analisemos entao a questao da determinacao de limites classicos e quanticos para a
massa do foton. Como mostramos no texto, os efeitos devidos a uma massa de repouso

nao nula para o foton podem ser incorporados no eletromagnetismo de uma maneira
bastante simples via a QED massiva. Neste esprito, duas interessantes implicacoes

136
decorrentes da possvel existencia de um foton massivo na natureza, a saber, pequenas
alteracoes nos valores conhecidos tanto do momento magnetico do eletron como da

deexao gravitacional da radiacao eletromagnetica, foram utilizados para estabelecer

limites superiores para a massa do foton. Os limites obtidos nao sao tao bons quanto

aqueles recentemente estimados, sendo porem comparaveis a outros limites existentes

e trazem novas contribuicoes para a questao da restricao da massa do foton. Podemos

responder agora, com total seguranca, que os objetivos principais do trabalho foram
plenamente atingidos. Para nalizar esta discussao gostaramos de externar nossa
opiniao sobre a questao da conjectura relativa as teorias de ordem superior, ou seja,

sobre a questao da renormalizabilidade versus unitariedade. Acreditamos, baseado nos

dois exemplos acima, que presentemente e impossvel conciliar numa mesma teoria de

ordem superior os quesitos unitariedade e renormalizabilidade. A discussao sobre a

renormalizabilidade da teoria de gravitacao de ordem superior em 4D apresentada por

Stelle em 1977 [6], que sem duvida continua sendo o modelo ideal para quem deseja se
aventurar neste importante campo de pesquisa, corrobora a nossa crenca.
Para nalizar, vamos comentar rapidamente as investigacoes que seriam, de certa

maneira, uma continuacao natural deste trabalho. Como vimos que nao podemos aco-
modar simultaneamente monopolos de Dirac e transformacoes de dualidade na eletro-

dinamica de Lee e Wick e nem em uma versao estendida da mesma, seria interessante
investigar a possibilidade desta coexistencia em eletrodinamicas com violacao da sime-

tria de Lorentz, tais como a de Meyers-Pospelov [7] e a de Gambini-Pullin [8]. Outro

ponto interessante seria vericar, usando como paradigma o trabalho de Stelle [6], a
questao da renormalizabilidade da teoria de gravitacao massiva em 3D.

137
Referencias Bibliogracas

[1] R. Cutkosky etal., Nucl. Phys. B12, 281 (1969).

[2] D. Boulware and D. Gross , Nucl. Phys. B233, 1 (1984).

[3] I. Oda, J. High Energy Phys. 05 (2009) 064.

[4] K. Muneyuki and N. Ohta, Phys. Rev. D bf 85, 101501 (R) (2012).

[5] N. Ohta, Classical Quantum Gravity bf 29, 015002 (2012).

[6] K. Stelle, Phys. Rev. D 16, 953 (1977).

[7] R. Meyers and M. Pospelov, Phys. Rev. Lett. 90, 211601 (2003).

[8] R. Gambini and J. Pullin, Phys. Rev. D 59, 124021 (1999).

138