Partie I : La relativité restreinte

Chapitre 1 : Rappels de mécanique newtonienne

I. Transformation de Galilée
La mécanique classique et la mécanique relativiste se construisent autour du principe de relativité du mouvement.
Considérons deux référentiels R et R' définis par (O; x, y , z, t ) et par (O ′; x ′, y ′, z ′, t ′) . Ce principe déclare que, si R' admet
la vitesse d’entraînement v e par rapport à R, alors R admet la vitesse − ve vis-à-vis de R'. Donc dans le référentiel R, O' est

animé d’une vitesse v e par rapport à O.

Soit un évènement aperçu en M par O et O'. Dans R, les coordonnées de M sont (x, y, z , t ) et dans R', ses coordonnées sont
(x ′, y ′, z ′, t ′) . Quelles relations existe-t-il entre (x, y, z, t ) et (x ′, y ′, z ′, t ′) ?
Pour simplifier mais tout en gardant la généralité du problème, on suppose que les axes Ox et O'x' sont confondus. Les axes
Oy et O'y', Oz et O'z' seront pris respectivement parallèles. Ensuite, la vitesse d’entraînement v e est supposée constante et
ayant la même direction que Ox.

Galilée nous propose la transformation évidente suivante :

x = x′ + ve t ′ x′ = x − ve t
y = y′ y′ = y
ou
z = z′ z′ = z
t = t′ t′ = t

II. Composition des vitesses

Elle s'obtient en différentiant les relations ci–dessus :

dx = dx'+ ve dt '
dt = dt '

dx dx'
vx = = + ve
dt dt '

v x = v′x + ve

On obtient de la même manière :

v y = v′y v z = v′z

1
III. Transformation des accélérations
On trouve, pour chacune des composantes de l'accélération :

a x = a′x a y = a′y a z = a′z

Ce que l'on peut écrire sous forme de vecteur :

r r
a = a′

Donc, la troisième loi de Newton s'écrit :

r r r
F = m a = m a′

F = F′

La loi de force est invariante par changement de référentiel galiléen, ce qui signifie qu’aucune expérience de mécanique ne
permet de différencier deux référentiels galiléens.

2
 Chapitre 2 : Établissement de la transformation de Lorentz

I. Résultat négatif de l'expérience de Michelson–Morley

L'expérience de Michelson–Morley consiste à vérifier l'addition des vitesses donnée par la cinématique newtonienne
lorsque l'une des vitesses est celle de la lumière. Selon cette expérience, il serait possible, en observant des figures
d'interférences, de mesurer une vitesse donnée par :

v =u+c

Le résultat de cette expérience nous dit que cette loi d'addition des vitesses est fausse, ce qui peut paraître contraire à
l'expérience quotidienne. Il faut donc trouver une nouvelle transformation qui tienne compte, qu'à faible vitesse, les vitesses
peuvent s'additionner, ce qui n'est plus vrai à grande vitesse.
Nous pouvons nous apercevoir que, d'après le résultat de cette expérience, la vitesse de la lumière est une vitesse limite,
identique dans tous les référentiels galiléens.
Nous allons établir de manière rigoureuse la transformation de Lorentz.

II. Démonstration de la transformation de Lorentz

Les hypothèses pour déterminer la transformation de Lorentz sont les suivantes :

– invariance de la vitesse de la lumière par changement de repère galiléen ainsi que sa constance et son indépendance vis–
à–vis du temps, de la source qui la crée et de l'observateur.
– homogénéité de l'espace, conduisant à des relations linéaires entre les deux référentiels.
– synchronisation des horloges (t = t ′ = 0 ) de R et R' lorsque les origines O et O' coïncident; les axes ont des directions
parallèles (disparition des termes constants).
– la vitesse relative v e , selon les axes Ox et O'x', n'influence pas les composantes y et z, c'est–à–dire y = y ′ et z = z ′ (ce
que nous vérifierons plus loin).

D'après ces hypothèses, nous pouvons écrire :

x′ = a1x + a2 ct et ct ′ = a3 x + a4 ct

avec a i des constantes pouvant dépendre de v e et de c.

La vitesse du point O' de R' mesurée dans le référentiel R est ve . Ce point admet comme coordonnée x' = 0 , ce qui donne :

a1x + a2 ct = 0

En différentiant, on obtient :

a1 dx = −a2 cdt
dx
a1 = − a2 c
dt

dx
Or est justement la vitesse de O' mesurée dans le repère R, donc :
dt

ve
a2 = −a1
c

On obtient alors :

x ′ = a1 ( x − v e t )

3
Considérons maintenant un rayon lumineux émis au point O à un instant pris comme origine (t = t ' = 0 ) , les deux origines
étant confondues à cet instant. Nous supposons que ce rayon lumineux se propage suivant l'axe Ox dans le sens des x
croissants. L'équation de propagation s'écrit :

– dans R : x = ct
– dans R' : x ′ = ct ′

Nous obtenons l'identité suivante, en croisant les membres des deux égalités :

x ct′ = x′ ct

ce qui entraîne, d'après l'expression de ct' et de ce que l'on vient d'obtenir :

x (a3 x + a4 ct ) = a1 (x − ve t )ct
⇔ a3 x 2 + a4 x ct = a1x ct − a1ve ct 2

Nous identifions terme à terme :

a1 = a 4

et en utilisant le fait que x = ct ,

a1ve
a3 = −
c

Nous avons donc obtenu la transformation suivante :

x′ = a1 (x − ve t )
 v 
ct ′ = a1  ct − e x 
 c 

Nous pouvons donc écrire :

 v 
x′ − ct′ = a1  x − vet − ct + e x 
 c 
 ve 
x′ − ct′ = a1 (x − ct ) 1 + 
 c 

 v 
D'après la relativité du mouvement, nous pouvons dire que : x − ct = a1 ( x ′ − ct ′ ) 1 − e 
 c 

En multipliant membre à membre, nous obtenons :

 v  v 
1 = a12 1 + e  1 − e 
 c  c 
soit :

1
a1 = ±
v e2
1−
c2

donc :

1
x′ = ± (x − v e t )
v e2
1−
c2

4
Pour déterminer le signe, nous effectuons l'approximation v e << c : x ′ = ± ( x − v e t ) . Le signe + doit être conservé pour que
la relation s'identifie à la transformation de Galilée lors des vitesses faibles.

En résumé, nous avons :

x ′ = γ e (x − v e t )
 v 
ct ′ = γ e  ct − e x 
 c 

−1 / 2
 v2 
avec γ e = 1 − e2  .
 c 
 

Démontrons maintenant le fait que y = y ′ et z = z ′ .

Supposons que la vitesse v e entraîne une diminution (ou une augmentation) des longueurs transversales. Pour cela,
imaginons deux cylindres creux dont les axes principaux coïncident avec l'axe Ox. Le premier cylindre s'étend à l'infini du
côté des x négatifs et est solidaire du référentiel R. Le second cylindre s'étend à l'infini du côté des x positifs et est solidaire
du référentiel R'.
Plaçons–nous dans le référentiel R. Le cylindre n°2 se déplace en direction du premier cylindre, et, par notre supposition, il
subit une diminution de son rayon et peut passer à l'intérieur du premier cylindre de manière permanente dans le temps.
Plaçons–nous maintenant dans le référentiel R'. Le cylindre n°1 semble se déplacer en direction du cylindre n°2 et subit une
diminution de son rayon. Le cylindre n°1 passe à l'intérieur du deuxième cylindre de manière permanente dans le temps.

Il y a contradiction selon le référentiel où nous nous plaçons et la seule possibilité qui puisse exister est : y = y ′ et z = z ′ .

Nous pouvons ainsi résumer la loi de transformation de Lorentz sous la forme usuelle suivante :

x ′ = γ e ( x − β e ct )
y′ = y ve
z′ = z
en posant β e =
c
(
et γ e = 1 − β e2 )−1 / 2
ct ′ = γ e (ct − β e x )

III. Autres écritures de la transformation de Lorentz

1. Écriture vectorielle

Nous posons :

r = x i + y j + zk r ′ = x ′ i + y ′ j + z ′k

La transformation s'écrit alors :

(
r = r ′ + (γ e − 1) β e ⋅ r ′ )β e
+ γ e ct ′ β e
βe
(
ct = γ e ct ′ + β e ⋅ r ′ )

5
2. Écriture hyperbolique

Nous introduisons la notion de rapidité re définie par :

ve
th re = = βe
c

Nous en déduisons :

−1 / 2
  2
 exp re − exp ( − re )  
γe = 1−   = ch re
  exp r + exp ( − r )  
  e e  


et γ e β e = ch re ⋅ th re = sh re

D'où :

x = x ′ ch re + ct ′ sh re
ct = ct ′ ch re + x ′ sh re

6
 Chapitre 3 : Bases de la relativité restreinte

I. Établissement de la composition des vitesses

La loi de composition des vitesses s'obtient en différentiant la transformation de Lorentz :

Nous avons :

dx = γ e ( dx′ + vedt ′)
 v 
c dt = γ e  c dt′ + e dx′ 
 c 

En divisant membre à membre et en multipliant par c, nous obtenons :

dx dx ′ + v e dt ′
=
dt v
dt ′ + e dx ′
c2

En divisant numérateur et dénominateur par dt ′ :

v ′x + v e
vx =
v v′
1+ e x
c2

La formule fournit bien à faible vitesse la loi d'addition galiléenne des vitesses : v x = v ′x + v e .
Calculons par exemple la vitesse de la lumière émise par une source animée d'une vitesse v :

v+c c 1+ v
c =c
( )
=
vc 1+ v
1+ c
c2

ce qui signifie que la vitesse de la lumière est invariante par changement de référentiel galiléen et indépendante de la vitesse
de la source.

Nous obtenons les transformations pour v y et v z :

v ′y v ′z
vy = vz =
 v v′   v e v ′x 
γ e 1 + e x  γ e 1 + 
 c2   c2 

II. Composition des rapidités

ve vx v′x
Soit th re = , th rx = et th rx′ = alors :
c c c

th re + th rx′
th rx = = th (rx′ + re )
1 + th re ⋅ th rx′

Nous avons donc :

rx = rx′ + re

Ainsi, contrairement à la mécanique classique où ce sont les vitesses qui s'ajoutent, ce sont les rapidités qui s'additionnent en
relativité restreinte.

7
III. Intervalle entre deux évènements
Nous allons montrer l'invariance de la quantité c 2 t 2 − x 2 − y 2 − z 2 :

c 2 t 2 − x 2 − y 2 − z 2 = γ e2 (ct ′ + β e x ′ ) 2 − γ e2 ( x ′ + β e ct ′ ) 2 − y ′ 2 − z ′ 2
( )
= γ e2 c 2 t ′ 2 + β e2 x ′ 2 − x ′ 2 − β e2 c 2 t ′ 2 − y ′ 2 − z ′ 2
= ( )
γ e2 1 − β e2 c 2 t ′ 2 − γ e2 1 − β e2 ( )x ′2 − y′2 − z′2
2 ′2 ′2 ′2 ′2
=c t −x −y −z

On définit l'intervalle entre deux évènements par :

2
s12 = c 2 (t 2 − t1 ) 2 − ( x 2 − x1 )2 − ( y 2 − y1 ) 2 − ( z 2 − z1 ) 2

On définit l'intervalle élémentaire par :

ds 2 = c 2 dt 2 − dx 2 − dy 2 − dz 2

c'est ce qu'on appelle la métrique de l'espace–temps de la relativité restreinte. Elle peut également s'écrire sous la forme :

ds 2 = c 2 dt 2 − dr 2 − r 2 dθ 2 − r 2 sin 2 θ dϕ 2
Il existe trois sortes d'intervalles :

– si ds 2 < 0 , l'intervalle est du genre espace.

– si ds 2 = 0 , l'intervalle est du genre lumière.
– si ds 2 > 0 , l'intervalle est du genre temps.

IV. Formalisme quadrivectoriel et introduction aux tenseurs

Le temps étant une composante à part entière, il est plus pratique d'utiliser dès lors des vecteurs particuliers de
dimension 4, appelés quadrivecteurs. Nous allons définir les quadrivecteurs position, vitesse, accélération et quantité de
mouvement. Les quadrivecteurs seront notés avec un tilde (~) pour ne pas les confondre avec d'autres quantités.
Un quadrivecteur possède une pseudo–norme définie de la manière suivante. Soit un quadrivecteur a~ = (a1 , a 2 , a 3 , a 4 ) . Le
carré de sa pseudo–norme est :
a~ = a 2 − a 2 − a 2 − a 24 1 2 3

1. Quadrivecteur position

Le quadrivecteur position a pour composantes les trois composantes d'espace et une de temps. Il est noté :

~
x = ( x, y , z , ct ) ou ~ ( )
x = r , ct

2. Quadrivecteur vitesse

L'intervalle élémentaire ds 2 peut se mettre sous la forme :

 dx 2 + dy 2 + dz 2   2 
ds 2 = c 2 dt 2 1 −  = c 2 dt 2 1 − v 
 c 2 dt 2   c2 
   

soit :
c dt
ds =
γ

On définit la durée propre élémentaire par :

ds dt
dτ = =
c γ

8
Le quadrivecteur vitesse est donné par :

d~
x  d r d (ct ) 
v~ = = ,
dτ  dτ dτ 
= γ v, γ c ( )
Le carré de sa pseudo–norme est :

 v2 
v~ = γ 2c 2 − γ 2v 2 = γ 2c 2 1 − 2  = c 2
 c 
 

3. Quadrivecteur accélération

Il est donné par la relation :

dv~ d
a~ = = (γ vr , γ c )
dτ dτ

Nous avons donc :

r r r
d (γ v ) d (γ v ) dv dγ r dγ dγ
=γ =γ2 +γ v et c =γc
dτ dt dt dt dτ dt

dγ
Il reste à calculer :
dt
−3 / 2
dγ 1  v 2   v  dv v dv
= − 1 −  − 2  =γ3
dt 2  c 2   c2  dt c 2 dt

Nous obtenons finalement :

 r dβ r 4 dβ 
a~ =  γ 2 a + γ 4 β v , γ cβ 
 dt dt 

4. Quadrivecteur quantité de mouvement et énergie de la particule libre

En relativité restreinte, le quadrivecteur quantité de mouvement est le produit de la masse (au repos) du corps considéré
par sa quadrivitesse.
~ r
p = mv~ = (γ mv , γ m c )

Le carré de sa pseudo–norme est :

~
p
2
= m 2 v~
2
= m 2c 2

La partie spatiale de la force se définit comme la dérivée par rapport au temps de la quantité de mouvement :

F=
dp
=
d γ mv ( )
dt dt

Calculons maintenant l'énergie de la particule libre. Par définition, nous avons :

r r
δE = F ⋅ dl
r
d (γ mv ) r
= ⋅ v dt
dt
r r
= v ⋅ d ( γ mv )

9
En développant le dernier terme, nous obtenons :

r r r r mc 2 c2
v ⋅ d (γ mv ) = γ mv dv + mv 2 dγ = dγ + mv 2 dγ car dv = dγ
γ2 γ 3v

r r  1 v 2 
v ⋅ d (γ mv ) = mc 2  + dγ = mc 2 dγ
 γ 2 c2 
 

Nous avons finalement :

δE = d γ m c 2 ( )
soit :

E = γ mc 2

C'est l'énergie de la particule libre.

Nous aurions pu trouver ce résultat de manière plus rapide en utilisant l'invariance de la pseudo–norme. Nous avons vu que
~ 2
p = m 2 c 2 donc :

d  ~  = 0
2
p
 
ce qui entraîne :

~
p ⋅ d~
p=0
r r
γ m c ⋅ d (γ m c ) − γ m v ⋅ d (γ m v ) = 0
[( ) r
γ m d γ m c 2 − v ⋅ d (γ m v ) = 0
r
]
c'est–à–dire :
r r
v ⋅ d (γ m v ) = d γ m c 2 ( )
D'où le résultat écrit ci–dessus :

E = γ mc 2
Faisons maintenant le développement limité à l'ordre 1 de cette expression pour vérifier que l'on retrouve l'énergétique
classique.

−1 / 2
 v 2   v2  2
γ m c 2 = m c 2 1 − ≈ m c 2 1 +  = m c 2 + mv
 c 2   2c 2  2
  

Nous retrouvons bien l'énergie cinétique de la mécanique classique mais avec un terme constant supplémentaire ne dépendant
que de la masse de la particule. Ainsi, même au repos, une particule possède une énergie, appelée énergie de masse donnée
par la formule la plus célèbre de la physique :

E0 = m c 2

C'est l'équivalence entre la masse et l'énergie. On dit que la matière est une forme organisée d'énergie.

Par extension, l'énergie cinétique d'une particule libre dans le cadre de la relativité restreinte est donnée par :

E k = (γ − 1)m c 2

10
Établissons maintenant deux relations entre la quantité de mouvement et l'énergie.
r r r
– Calculons la quantité E 2 − p 2 c 2 , avec p = γ mv :

r  v2 
E 2 − p 2 c 2 = γ 2 m 2 c 4 − γ 2 m 2 v 2 c 2 = γ 2 m 2 c 4 1 − 
 c2 
 

r
E 2 − p 2c 2 = m 2c 4

Ainsi, pour un photon :

E = pc

– Établissons le lien entre E k et pc :

p 2c 2 = E 2 − m 2c 4
[ ]
= γ 2 m 2 c 4 + (γ − 1) 2 m 2 c 4 − (γ − 1) 2 m 2 c 4 − m 2 c 4
[ 2
= m 2 c 4 γ 2 − (γ − 1) − 1 + E k2]
= 2 E k m c 2 + E k2

Nous avons donc :

(
pc = E k E k + 2 m c 2 )
5. Le tenseur métrique de la relativité restreinte

Le carré de l'intervalle élémentaire ds 2 possède une forme particulière dans le cas de la relativité restreinte, qui est :

ds 2 = c 2 dt 2 − dx 2 − dy 2 − dz 2

mais une forme plus générale est donnée par :

ds 2 = η11dx 2 + η12dx dy + η13dx dz + η14 dx cdt

+ η 22dy 2 + η21dy dx + η 23dy dz + η 24dy cdt
+ η33dz 2 + η31dz dx + η32 dz dy + η34 dz cdt
+ η 44c 2 dt 2 + η 41cdt dx + η 42cdt dy + η43cdt dz

ou, en posant x µ = ( x, y , z , ct ) , où µ est un indice :

ds 2 =
∑η
µν
µν dx µ dxν

Le produit de dx par dy étant commutatif, il en résulte que η12 = η 21 , et de même pour les autres indices, donc :

ds 2 = η11dx 2 + 2η12 dx dy + 2η13 dx dz + 2η14 dx cdt

+ η 22 dy 2 + 2η 23dy dz + 2η 24 dy cdt
+ η33dz 2 + 2η34 dz cdt
+ η 44 c 2 dt 2

Les axes étant choisis orthogonaux, les termes croisés disparaissent et il ne reste plus que la métrique de la relativité
restreinte, écrite sous cette forme :
ds 2 = η 44c 2 dt 2 + η11dx 2 + η22 dy 2 + η33dz 2

11
La quantité η µν est appelée tenseur métrique de la relativité restreinte et peut s'écrire sous forme matricielle :

 −1 0 0 0
 
0 −1 0 0
η µν = 
 0 0 −1 0
 
 0 0 0 1 

V. Avance du périhélie de Mercure en relativité restreinte

1. Équation de mouvement en mécanique classique

Une particule soumise à la force gravitationnelle du Soleil obéit à l'équation différentielle suivante :

d 2u GM
2
+u =
dϕ C2

1
Dans cette équation, on a posé u = , où r est la distance de la particule au Soleil, ϕ l'angle qui définit avec r la position de
r
la particule, G, M et C sont respectivement la constante de la gravitation universelle, la masse du Soleil et la constante des
aires.
C2
La solution de cette équation est une ellipse dont le paramètre p est .
GM

L'application de cette équation à l'orbite de Mercure et la prise en compte des perturbations des autres planètes aboutit à une
précession du périhélie d'une valeur de 531" par siècle, au lieu des 574" observés. Calculons la précession de Mercure dans le
cadre de la relativité restreinte.

2. Équation de mouvement en relativité restreinte

L'énergie totale s'écrit :

E = γ mc 2 + E p

Soit :

E − Ep
γ =
mc 2
r
La position de la particule s'écrit : r = r er , la vitesse vaut alors :

r
v = r&er + rθ& eθ

Le carré de la vitesse est donc :

v 2 = r& 2 + r 2θ& 2

On en déduit :

1 v2 1
γ 2
= 1−
c2
= 1−
c2
(r& 2
+ r 2θ&2 )

1 dr dr du & θ& du
En posant comme précédemment u = , on a : r& = = θ =− 2
r dt du dθ u dθ

12
Soit :

1  θ&2  du  θ& 2  1
2
1− 2  4 dθ
  + 2 = 2
c u   u  γ

Pour simplifier cette expression, on fait appel à la conservation du moment cinétique. Celui–ci s'écrit :
r r r r r
L = r∧ p = r∧ (γ mv ) = γ mr 2θ& ez

Ainsi, on peut extraire θ& de l'expression de L :

Lu 2
θ& =
γm

On obtient alors :

L2   du   1
2
1
1− 
2  dθ
 + u2  = 2
2
c γ m 2
  γ

ce qui donne :
2
L2  du  
2
2   E − E p 
γ 2 = 1+   + u =
m 2c 2  dθ    mc 2 
2  2 
 du  m 2c 2  E − E p 
  + u 2 = 2   − 1
 dθ  L  mc 2  
 

Dérivons cette dernière expression par rapport à θ :

du  d 2u  2m  E − E p  dE p dr
2 + u  = − 2  
 dr dθ

dθ  dθ 2  L  mc 2
 

du  d 2u  2m  E − E p  dE p du
2 + u  = 2 2  
 dr dθ
dθ  dθ 2  Lu
  mc
2


du
Il apparaît la solution triviale = 0 , qui correspond à une trajectoire circulaire. Intéressons–nous à l'autre solution :
dθ

d 2u m  E − Ep  dE p
+u =  
dθ 2
L u  mc 2
2 2  dr


GMm
L'énergie gravitationnelle s'écrit : E p = − , donc
r

dE p
= GMmu 2
dr

d 2u GMm 2  E + GMmu 
+u =  
dθ 2 L2  mc 2 

d 2u GMm 2 E
+ (1 − ε )u =
dθ 2 L2 mc 2

2
 GMm 
avec ε =   .
 Lc 

13
 GMm2
Le second membre de l'équation différentielle est constant et s'identifie à pour de petites vitesses. En effet,
L2
E v2 GM
2
≈ 1+ − .
mc 2c 2 rc 2

GMm 2 E 1−ε
On peut poser 2 2
= :
L mc p

d 2u 1− ε
2
+ (1 − ε )u =
dθ p

La solution de cette équation est de la forme :

u (θ ) =
1 + e cos ( )
1 − ε (θ − θ 0 )
p

qui est une ellipse dont le demi–grand axe tourne lentement.

On pose θ0 = 0 . Au bout d'une révolution, on obtient : 1 − ε θ = 2π , soit :

θ = 2π (1 − ε )−1 / 2 ≈ 2π + επ

L'avance supplémentaire vaut donc :

2
 GMm 
∆θ = επ = π  
 Lc 

GMm2 1
Or 2
= , donc :
L p

πGM
∆θ =
c2 p

Pour Mercure, p = 5,546.1010 m . Avec M = 1,989.1030 kg , nous trouvons :

∆θ = 8,364.10−8 rad

La période de révolution de Mercure étant de 87.9 jours, nous avons pour un siècle :

∆θ × 100 × 365.25
δ = ≈ 3,476.10−5 rad
87.9

Soit :

δ = 7,17"

On remarque que la précession prévue par la relativité restreinte ne correspond pas à l'avance observée. L'avance de 43" par
siècle sera expliquée par la relativité générale.

14
 Partie II : L'électromagnétisme

Chapitre 1 : Outils mathématiques

I. Analyse vectorielle
1. L'opérateur gradient

Soit une fonction scalaire g (x, y , z ) = g (M ) indépendant du temps. Le gradient est défini de la manière suivante :

r r
( )
dg = g M + dM − g (M ) = gradg ⋅ dM

Le gradient est donc un vecteur dont les composantes sont, dans un repère cartésien :

 ∂g ∂g ∂g 
gradg =  , , 
 ∂x ∂y ∂z 
 
r
Il est aussi noté ∇g .

r
Si G = gradg , alors :

B r r B r B

∫ A
G ⋅ dM =
∫ A
gradg ⋅ dM =
∫ A
dg = g (B ) − g ( A)

L'intégrale est alors indépendante du chemin choisi.

2. L'opérateur divergence
r r
Soit un vecteur A de composantes (A x , A y , A z ) . La divergence de A est, dans un repère cartésien :

r ∂A x ∂A y ∂A z
div A = + +
∂x ∂y ∂z

r
L'opérateur divergence est noté aussi ∇ ⋅ , mais nous n'emploierons pas cette notation.

3. L'opérateur rotationnel
r
Le rotationnel d'un vecteur A = (A x , A y , Az ) s'écrit :

r r
B = rot A

r
Les composantes de B sont, dans un repère cartésien :

r  ∂A ∂A y ∂A x ∂A z ∂A y ∂A x 
B= z − , − , − 
 ∂y ∂z ∂z ∂x ∂x ∂y 


15
4. Le laplacien et le dalembertien :

Le laplacien scalaire dans un repère cartésien est défini par :

∂2g ∂2g ∂2g

∆g = + +
∂x 2 ∂y 2 ∂z 2

Le laplacien vectoriel est défini par :

r r r
r ∂2 A ∂2 A ∂2 A
∆A = + +
∂x 2 ∂y 2 ∂z 2

Le dalembertien s'écrit :
r
r ∂2A r

A= − ∆A
c 2 ∂t 2

avec ici c la vitesse de la lumière.

5. Relations entre les opérateurs

Les propriétés des opérateurs de l'analyse vectorielle sont très utilisées en électromagnétisme. Nous donnons ici
les principales relations entre les différents opérateurs définis ci–dessus.

(
∆f = div grad f )
( ) r
rot grad f = 0
( ) r
div rot A = 0
( ) r r
(r
rot rot A = grad div A − ∆A )
r r r
E = − gradV ⇔ rot E = 0
r r r
div B = 0 ⇔ B = rot A

II. Flux et circulation d'un vecteur– Théorèmes de l'analyse vectorielle

1. Flux d'un champ de vecteurs

r r
Le flux élémentaire dφ d'un vecteur A à travers une surface orientée dS est défini par :
r r
dφ = A ⋅ dS

Nous pouvons comparer le flux au débit d'un fleuve, qui est la quantité d'eau qui passe à travers une surface.

2. Circulation d'un vecteur

r
La circulation V d'un vecteur E le long d'une courbe quelconque est :

r r
V =
∫ E ⋅ dl
C

16
3. Théorème de Green–Ostrogradski
r
Soit un élément de volume dτ = dx dy dz . Le flux d'un champ de vecteurs A à travers la surface inférieure
 0 
r  
dS =  0  est :
 − dx dy 
 

dφ inf = − A z dx dy

r
où A z est la projection du vecteur A selon l'axe z.

De même, sur la surface supérieure :

 ∂A 
dφ sup =  A z + z dz  dx dy
 ∂z 

Le flux à travers les surfaces inférieures et supérieures est donc :

∂A z ∂A z
dφ inf + sup = dx dy dz = dτ
∂z ∂z

En effectuant un raisonnement similaire sur les quatre autres surfaces, nous en déduisons que le flux à travers
toutes les surfaces élémentaires est :

 ∂A ∂A y ∂A z 
dφ =  x + +  dτ
 ∂x ∂y ∂z 
 

r r r
dφ = A ⋅ dS = div A⋅ dτ

Le flux total est donc :

r r r
φ=
∫∫
S
A ⋅ dS =
∫∫∫
V
div A ⋅ dτ

Ceci constitue le théorème de Green–Ostrogradski.

4. Théorème de Stokes
r
Soit un contour fermé PQRS. La circulation d'un vecteur A = (A x , A y , Az ) le long de PQ est A x dx . Le long de
 ∂A  ∂A
RS, elle vaut : −  A x + x dy  dx . La somme est donc : − x dx dy .
 ∂y  ∂y
∂A y
De même, le long de QR et SP, nous avons : dx dy .
∂x

Le long du contour fermé PQRS, la circulation élémentaire est :

 ∂A y ∂Ax 
 −  dx dy
 ∂x ∂y 
 

17
Nous pouvons étendre ce résultat pour un contour fini quelconque et à la surface limitée par ce contour. La
r
circulation d'un vecteur A est égale au flux du rotationnel de ce vecteur :

r r r r
∫
C
A ⋅ dl =
∫∫
S
rot A ⋅ dS

C'est le théorème de Stokes.

III. Introduction aux tenseurs

1. Tenseur de Levi–Civita

Dans le cas de l'électromagnétisme, nous utilisons un tenseur à 4 indices qui ressemble au symbole de
Kronecker.

Nous posons :

1 si deux indices subissent une permutation paire

µνρσ 
ε = − 1 si deux indices subissent une permutation impaire
0 si deux indices sont identiques


Nous avons par exemple :

ε 1234 = ε 1423 = 1
ε 1243 = −1
ε 1224 = 0

2. Tenseurs d'ordre deux symétriques ou antisymétriques

Un tenseur quelconque d'ordre deux est décomposable en une partie symétrique et une partie antisymétrique :

1
T µν = (Tµν + Tνµ ) + 1 (Tµν − Tνµ )
2 2

1
Le tenseur sera dit symétrique si T µν = Tνµ donc T µν = (Tµν + Tνµ ) .
2

1
Le tenseur sera dit antisymétrique si T µν = −Tνµ donc T µν = (Tµν − Tνµ ) .
2

Remarque : Un tenseur antisymétrique a la propriété suivante : T µµ = 0 .

18
 Chapitre 2 : Électromagnétisme de Maxwell

I. Électrostatique–Magnétostatique

1. Équation de la conservation de la charge

r
Le vecteur densité de courant j est proportionnel à la charge volumique et à la vitesse de ces charges.
r r
j = −N q v

où N représente le nombre de charges par unité de volume et q la charge portée par une particule.

L'intensité du courant est définie comme le flux du vecteur densité de courant. Nous pouvons donc écrire :
r
r r dl r dτ
dφ = j ⋅ dS = − N q dS = − N q
dt dt

en posant N q dτ = ρ dτ . Nous obtenons alors :

r r d
∫∫
S
j ⋅ dS = −
dt ∫∫∫ ρ dτ
V

Appliquons le théorème de Green :

r r r ∂ρ
∫∫ j ⋅ dS = ∫∫∫ div j ⋅ dτ = −∫∫∫ ∂t dτ
S V V

d'où :

r ∂ρ
div j = −
∂t

Cette équation est appelée équation de la conservation de la charge ou équation de continuité.

2. Électrostatique–Théorème de Gauss

La circulation d'un champ électrostatique le long d'un contour fermé est nulle, c'est–à–dire que le champ
électrostatique est conservatif.

r r
∫
C
E ⋅ dl = 0

Cette relation équivaut à :

r r
rot E = 0

c'est une première équation de Maxwell, pour le champ électrostatique.

19
 r
Cette relation entraîne que E dérive d'un potentiel (scalaire) :
r
E = −grad V

r
Théorème de Gauss : Le flux du champ électrostatique (valable nous le verrons pour un champ électrique) E à
travers une surface fermée S est proportionnel au nombre de charges contenues dans le volume délimité par la
surface S.

r r 1
φ=
∫∫ E ⋅ dS = ε
S
0
Qint

Avec Q int =
∫∫∫ ρ dτ , nous pouvons exprimer la forme locale du théorème de Gauss de la manière suivante :
V

r r r 1
∫∫
S
E ⋅ dS =
∫∫∫
V
div E ⋅ dτ =
ε0 ∫∫∫ ρ dτ
V

Ainsi, nous avons :

r ρ
div E =
ε0

ce qui constitue une deuxième équation de Maxwell.

r
Or E = − gradV donc :

ρ
∆V = −
ε0

Cette dernière relation est appelée équation de Poisson.

Le théorème de Gauss nous permet de retrouver la loi de Coulomb :

À une distance r d'une source ponctuelle Q, nous avons :

r r Q
∫∫
S
E ⋅ dS = E ⋅ 4π r 2 =
ε0

soit :
Q
E=
4π ε 0 r 2

3. Magnétostatique–Théorème d'Ampère

La présence de monopôles magnétiques n'ayant pas (encore) été détecté, le flux du champ magnétique est nul.

r r
∫∫ S
B ⋅ dS = 0

ce qui entraîne :
r
div B = 0

ce qui constitue une troisième équation de Maxwell.

20
Ainsi, nous avons :
r r
B = rot A

r
où A est appelé potentiel vecteur.

Théorème d'Ampère : La circulation du champ magnétique le long d'un contour fermé est proportionnelle à
l'intensité du courant dans ce contour.

r r
∫
C
B ⋅ dl = µ 0 I

r r
où I =
∫∫ j ⋅ dS .
S

Nous pouvons alors écrire :

r r r
∫∫ rot B = µ ∫∫ j ⋅ dS
S
0
S

c'est–à–dire :
r r
rot B = µ 0 j

Il s'agit d'une dernière équation de Maxwell, pour le champ magnétostatique.

Cette dernière relation peut s'écrire :

( ) r r
rot rot A = µ 0 j
r r r
( )
grad div A − ∆A = µ 0 j

r r
Nous pouvons imposer à A une jauge, dite jauge de Coulomb : div A = 0 , de manière à avoir :
r r
∆A = − µ 0 j

Résumons les équations de Maxwell des champs électrostatique et magnétostatique :

r ρ r
div E = div B = 0
ε0

r r r r
rot B = µ 0 j rot E = 0

II. Électromagnétisme

1. Équations de Maxwell

L'expérience montre que la variation temporelle du flux d'un champ magnétique crée un courant électrique qui
s'oppose au flux, nommé courant induit :

dφ d r r
V =−
dt
=−
dt ∫∫
S
B ⋅ dS

21
Nous avons, par définition :

r r r r
V =
∫ E ⋅ dl = ∫∫ rot E ⋅ dS
C S

Ainsi,
r
r r ∂B r
∫∫
S
rot E ⋅ dS = −
∫∫
S
∂t
⋅ dS

c'est–à–dire,
r
r ∂B
rot E = −
∂t

Nous en déduisons de ces mêmes relations une définition du champ électrique :

r r d r r d r r d r r
∫
C
E ⋅ dl = −
dt ∫∫S
B ⋅ dS = −
dt ∫∫
S
rot A ⋅ dS = −
dt ∫
C
A ⋅ dl

c'est–à–dire en comparant le premier membre et le dernier membre de l'égalité,

r
 r ∂A  r
∫
C
E +


 ⋅ dl = 0
∂t 

d'où
r
r ∂A
E = −gradV −
∂t

r r
En magnétostatique, nous avions : rot B = µ 0 j , ce qui entraîne, d'après div rot B = 0 : ( r)
r
div j = 0

Mais cette relation, qui était vraie en régime statique ou lentement variable, ne l'est plus maintenant. Cherchons à
remplacer les équations de Maxwell de la magnétostatique et de l'électrostatique afin qu'elles soient valables en
électromagnétisme.
r
L'équation div B = 0 n'est pas à modifier, puisqu'elle ne dépend que de la présence éventuelle de monopôles. Les
seules équations qui peuvent subir une transformation sont donc :

r ρ r r
div E = rot B = µ 0 j
ε0

Nous allons les rendre dépendantes du temps en les écrivant :

r
r ρ ∂f r r ∂F
div E = + rot B = µ 0 j +
ε 0 ∂t ∂t

r
où f est une fonction scalaire et F une fonction vectorielle dépendante du temps.

22
Calculons tout d'abord la divergence de la seconde relation (qui est nul par définition) :

( )r r ∂
div rot B = µ0 div j +
∂t
r
div F ( )
∂ r
=
∂t
(
− µ0 ρ + div F = 0 )
Ainsi :
r
div F = µ 0 ρ

r r ρ ∂f
Appliquons la définition de E dans la relation div E = + :
ε 0 ∂t

r
 ∂A  ρ ∂f
div  − gradV − = +
 ∂t  ε 0 ∂t

r ρ
⇔ − ∆V −
∂
∂t
(
div A + f = )
ε0

Nous pouvons alors appliquer la jauge suivante :

r
div A = − f

r
Nous pouvons écrire f = 0 sans que A ne soit changé, nous obtenons alors :

r ρ
div E =
ε0

r ρ
Ainsi, en électromagnétisme, l'équation div E = est conservée. À partir de cette relation, nous pouvons
ε0
r
établir l'expression de F :

r 1 r
ρ = ε 0 div E = div F
µ0

r r
F = µ 0ε 0 E

Ainsi :
r
r r ∂E 

rot B = µ 0 j + ε 0
 ∂t 


Nous avons les quatre équations de Maxwell de l'électromagnétisme :

r ρ r
div E = div B = 0
ε0

r r
r r ∂E  r ∂B

rot B = µ 0 j + ε 0 rot E = −
 ∂t  ∂t


r
Remarquons que ces quatre équations ne sont pas indépendantes. D'une part, nous avons construit rot B à partir
r
de div E . D'autre part, nous pouvons montrer que les deux autres équations sont intimement liées.

23
D'après l'analyse vectorielle, nous savons que div rot E = 0 . ( r)
r
 ∂B 
( r)
Or div rot E = −div  
( r
)
 = − ∂ div B , ce qui entraîne :
 ∂t  ∂t

r
div B = 0

2. Propagation des champs électriques et magnétiques dans le vide

Les équations de Maxwell dans le vide s'écrivent :

r r
div E = 0 div B = 0

r r
r ∂E r ∂B
rot B = µ 0 ε 0 rot E = −
∂t ∂t

( r)
Calculons rot rot E :
r
∂2E
( )
r ∂ r
( )
rot rot E = − rot B = − µ 0 ε 0
∂t ∂t 2

( r) ( r
) r r
or rot rot E = grad div E − ∆E . Comme div E = 0 , nous obtenons :

r
r ∂2E
∆E + µ 0 ε 0 =0
∂t 2

Il s'agit de l'équation de d'Alembert qui décrit la propagation des ondes progressives. La vitesse à laquelle se
1
déplace le champ électrique est , qui n'est rien d'autre que la vitesse de la lumière c.
µ 0ε 0

r
Nous pouvons faire le même calcul pour B et obtenir :
r
r ∂2B
∆B + µ 0 ε 0 =0
∂t 2

3. Jauge de Lorenz

Nous avons :
r r
r ∂A r r r r ∂E 
E = − gradV − B = rot A 
rot B = µ 0 j + ε 0
∂t  ∂t 


De cette dernière relation, nous en déduisons :

r
r ∂  ∂A  
( )
r

rot rot A = µ 0 j + ε 0
 ∂t 
− gradV −
∂t  

r
r r r 1  ∂V  1 ∂ 2 A
( )
⇔ grad div A − ∆A = µ 0 j − grad  −
c2  ∂t  c 2 ∂t 2
r
∂2 A r r  r ∂V 
⇔ − ∆A = µ 0 j − grad  div A + 
2 2
c ∂t  c 2 ∂t 

24
c'est–à–dire :

r r  r ∂V 

A = µ 0 j − grad div A + 
 c 2 ∂t 

r r ∂V
Comme A et V ne sont pas totalement définis, nous pouvons choisir comme jauge div A + = 0 . C'est la
c 2 ∂t
jauge de Lorenz.

Nous avons ainsi, sous la jauge de Lorenz :

r r

A = µ0 j

Cherchons une expression similaire pour le potentiel scalaire.

r ρ
Nous avons div E = donc :
ε0

r
 ∂A  ρ
div  − gradV − =
 ∂t  ε 0

r ρ
⇔ ∆V +
∂
∂t
(
div A = −
ε0
)
r ∂V
Donc, d'après la jauge de Lorenz div A = − , nous obtenons :
c 2 ∂t

∂ 2V ρ
∆V − =−
c 2 ∂t 2 ε0

c'est–à–dire :

ρ

V =
ε0

25
 Chapitre 3 : Électromagnétisme et relativité

I. Apport du formalisme quadridimensionnel

Le but est d'utiliser le formalisme quadridimensionnel pour simplifier les équations de l'électromagnétisme.

1. Quadrivecteur densité de charges

~
Nous pouvons construire un quadrivecteur j , caractérisant la densité de courant et la densité volumique de
charges de la manière suivante :

~ r
j = jµ = j, ρ c ( )
La constante c (vitesse de la lumière) est ici pour garder toutes les composantes du quadrivecteur homogènes.

2. Quadrivecteur potentiel

Dans le formalisme du précédent chapitre, nous avons croisé le potentiel scalaire et le potentiel vecteur. Nous
pouvons former le quadrivecteur potentiel :

~ r V
A = A µ =  A, 
 c

~
Le quadrivecteur potentiel A peut s'exprimer à l'aide de ses composantes covariantes :

~  r V
A = Aµ =  − A, 
 c

3. L'opérateur dérivée partielle

Nous construisons un opérateur, qui généralise l'opérateur gradient, défini de la manière suivante :

r ∂ 
∂ µ =  ∇, 
 c∂t 

On a ainsi :
 r ∂ 
∂ µ =  − ∇, 
 c ∂t 

L'opérateur dérivée partielle se transforme, selon la relativité restreinte :

(
∂1′ = γ e ∂1 − β e ∂ 4 )
2′ 2
∂ =∂
3′
∂ = ∂3
(
∂ 4 ′ = γ e ∂ 4 − β e ∂1 )

26
4. Transformation du quadrivecteur potentiel

Entre deux référentiels galiléens, la transformation du quadrivecteur potentiel s'écrit :

 V′
Ax = γ e  A′x + β e 
 c 
Ay = A′y
Az = Az′
V V′ 
= γ e  + β e A′x 
c  c 

5. Tenseur de Faraday

Le tenseur de Faraday est un tenseur de deuxième ordre, défini par :

Fµν = ∂ µ Aν − ∂ν Aµ

~
c'est le commutateur des dérivées premières de A . Par sa définition, le tenseur de Faraday est antisymétrique.

Ses formes covariante et contravariante s'écrivent sous forme matricielle :

 0 −Bz By −E x / c   0 −Bz By E x / c
   
B 0 − B − E y / c B 0 − B E y / c
Fµν =  F µν = 
z x z x
 − By Bx 0 − Ez / c  − By Bx 0 Ez / c
   
 E x / c E y / c E z / c 0   − E x / c − E y / c − E z / c 0 

6. Transformation einsteinienne des champs électrique et magnétique

(
Ex′ = cF 1′ 4′ = c ∂1′ A4′ − ∂ 4′ A1′ )
V   V
( ) (
Ex′ = c γ e2 ∂1 − βe∂ 4  − β e A1  − cγ e2 ∂ 4 − β e∂1  A1 − β e  )
c   c
  V   V V  
( )
Ex′ = c γ e2  1 − β e2  ∂1 − ∂ 4 A1  + βe  ∂ 4 − ∂ 4  + βe ∂1 A1 − ∂1 A1  ( )
  c   c c 
 1V 
E′x = c  ∂ − ∂ 4 A1  = cF 14
 c 

d'où finalement :
Ex′ = E x

Nous avons ensuite :

(
E′y = c F 2′4′ = c ∂ 2′ A4′ − ∂ 4′ A2′ )
V 
 c 
(
E′y = c γ e∂ 2  − β e A1  − cγ e ∂ 4 − β e∂1 A2 )
 2V 
E′y = c γ e  ∂
 c
2
(
− ∂ A + β e ∂1 A2 − ∂ 2 A1 
4

)
(
E′y = c γ e F 24 + β e F 12 )
Soit :
(
E′y = γ e E y − β ecBz )
De la même manière, nous aurions trouvé :

(
E′z = γ e E z + βecB y )

27
Nous calculons ensuite :

Bx′ = F 3′2′ = ∂ 3′ A2′ − ∂ 2′ A3′

Bx′ = ∂3 A2 − ∂ 2 A3 = F 32

Bx′ = Bx

Nous avons :
B′y = F 1′ 3′ = ∂1′ A3′ − ∂ 3′ A1′
 V
( )
B′y = γ e ∂1 − β e∂ 4 A3 − γ e∂3  A1 − β e 
 c
  V 
B′y = γ e  ∂1 A3 − ∂3 A1 + β e  ∂3 − ∂ 4 A3  
  c 
(
B′y = γ e F 13 + β e F 34 )
 E 
B′y = γ e  B y + βe z 
 c 

Nous aurions de même :

 Ey 
B′z = γ e  Bz − β e 
 c 

En résumé, nous avons :

E′x = Ex cBx′ = cBx

(
E′y = γ e E y − β ecBz ) (
cB′y = γ e cB y + βe E z )
(
E′z = γ e E z + β ecB y ) (
cB′z = γ e cBz − β e E y )
Nous avons ci-dessous les expressions vectorielles de la transformation du champ électromagnétique :

E//′ = E// c B//′ = c B//

(
E⊥′ = γ e E⊥ + c β e ∧ B⊥ ) (
c B⊥′ = γ e c B⊥ − βe ∧ E⊥ )
Nous en déduisons les invariants suivants :

( )( )
E ⋅ c B = E// + E⊥ ⋅ c B// + B⊥ = E// ⋅ c B// + E⊥ ⋅ c B⊥

or on a : E// ⋅ c B// = E//′ ⋅ c B//′ . On peut montrer aussi E⊥ ⋅ c B⊥ = E⊥′ ⋅ c B⊥′ , donc :

E ⋅ c B = E ′ ⋅ c B′

Un autre invariant important est E 2 − c 2 B 2 .

(
E 2 − c 2 B 2 = E// + E⊥ ) − c (B
2
2
// + B⊥ )
2

E 2 − c 2 B 2 = E//2 − c 2 B//2 + E⊥2 − c 2 B⊥2

E 2 − c 2 B 2 = E//′2 − c 2 B//′2 + E⊥′2 − c 2 B⊥′2
E 2 − c 2 B 2 = E ′ 2 − c 2 B′ 2

28
II. Équations de Maxwell sous forme tensorielle
Comme les composantes du tenseur de Faraday sont celles des champs magnétique et électrique, nous devinons
que les équations de Maxwell, écrites avec ce formalisme, vont faire intervenir les dérivées premières de ce
tenseur.

Les deux équations, indépendantes l'une de l'autre, qui sont identiques aux quatre équations de Maxwell, sont les
suivantes :

∂ µ F µν = µ 0 j ν

ε µνρσ ∂ν Fρσ = 0

r
avec j µ = ( j , ρ c ) .

Vérifions que nous retrouvons bien les équations de Maxwell "classiques" avec ces deux formules.

Dans la première de ces deux relations, faisons ν = 4 :

∂ µ F µ4 = µ0 j 4 = µ0 ρ c
⇔ ∂ 1 F 14 + ∂ 2 F 24 + ∂ 3 F 34 = µ 0 ρ c
E  E  E 
⇔ ∂x x +∂y y +∂z z  = µ0 ρ c
     
 c   c   c 
⇔ ∂ x E x + ∂ y E y + ∂ z E z = µ0c 2 ρ

r ρ
⇔ div E =
ε0

De même, en faisant ν = 1 , nous trouvons :

∂ µ F µ1 = µ 0 j 1
⇔ ∂ 2 F 21 + ∂ 3 F 31 + ∂ 4 F 41 = µ 0 j 1
1
⇔ ∂ y Bz − ∂ z B y − ∂t Ex = µ0 jx
c2
[ ]r
⇔ rot B x = µ 0 j x + µ 0ε 0 ∂ t E x

En faisant pour ν = 2 et ν = 3 , nous trouvons finalement :

r
r r ∂E 

rot B = µ 0 j + ε 0
 ∂t 


Le deuxième groupe d'équations se trouve en procédant de la même manière.

Sachant que ∑ε
ρ ≠σ
µνρσ =2
∑ε
ρ <σ
µνρσ , nous avons en posant µ = 4 :

(
2 ε 4123 ∂ 1 F23 + ε 4213 ∂ 2 F13 + ε 4312 ∂ 3 F12 = 0)
⇔ ∂ x Bx + ∂ y B y + ∂ z Bz = 0

r
⇔ div B = 0

29
En prenant µ = 1 , nous avons :

( )
2 ε 1234 ∂ 2 F34 + ε 1324 ∂ 3 F24 + ε 1423 ∂ 4 F23 = 0
⇔ −∂ y E z + ∂ z E y − ∂ t B x = 0

[ ] r
⇔ rot E x = −∂ t B x

Ainsi,
r
r ∂B
rot E = −
∂t

r
Nous avons donc les équivalences suivantes, avec j µ = ( j , ρ c ) :

 r ρ
div E =
 ε 0
∂ µ F µν = µ 0 j ν ⇔ r
rot B = µ  rj + ε ∂E 
r
0 0

  ∂t 

r
div B = 0
 r
ε µνρσ ∂ν Fρσ = 0 ⇔  r ∂B
rot E = −
 ∂t

Notons également que le tenseur dual s'écrit par définition :

1 µνρσ
∗
F µν = ε Fρσ
2

ce qui revient à écrire les équations de Maxwell sous la forme :

∂ µ F µν = µ 0 j ν

∂ µ ∗ F µν = 0

Jauge de Lorenz :

Nous avions précédemment, sous la jauge de Lorenz :

r r

A = µ0 j
ρ

V =
ε0
V  ρc r r
Nous pouvons écrire maintenant, avec
 =
2
= µ 0 ρ c et
A = µ 0 j :
 c  c ε0

Aµ = µ0 j µ

30
 Partie III : La relativité générale

Chapitre 1 : Analyse tensorielle

I. Vecteurs contravariants, vecteurs covariants

r
Un vecteur contravariant, noté x = x i e i , est le vecteur polaire classique. L'usage veut que l'on mette ses
composantes avec un indice "en haut".Il y a sommation sur l'indice répété 2 fois, mais le signe ∑ a été omis
volontairement (convention d'Einstein).

Lorsque la base à laquelle est rapporté le vecteur est quelconque, il est possible de décomposer le vecteur selon
d'autres composantes que les composantes contravariantes, ce sont les composantes covariantes. Celles–ci sont
définies de la manière suivante, avec un indice "en bas" :

r
x j = x ⋅e j

Le lien entre les deux types de composantes est établi par :

r
( ) ( )
x ⋅ e j = x i e i ⋅ e j = x i e i ⋅ e j = x i g ij = x j

x j = gij xi

Le produit scalaire étant commutatif, nous avons g ij = g ji . De même, nous pourrions montrer que x j = g ij x i et
g ij = g ji
.

Lorsque la base est orthonormée, nous avons g ij = g ij .

Nous pouvons définir le produit scalaire :

r r
( )
x ⋅ y = x i y j ei ⋅ e j = gij xi y j
= xjy j

On a aussi :
r r
x ⋅ y = x i yi

L'objet g ij est appelé tenseur fondamental ou tenseur métrique.

II. Symboles de Cristoffel

Lorsqu'on définit une base associée à un point M, on dit que celle–ci est une base naturelle. Les vecteurs
unitaires linéairement indépendants sont définis de la manière suivante :

∂OM
ek =
∂u k

où u k sont les coordonnées curvilignes du point M par rapport à un repère cartésien fixe d'origine O.

31
Définissons la différentielle du vecteur OM :

∂ OM
d OM = du k = e k du k
∂u k

Lorsque l'espace dans lequel on travaille n'est pas euclidien, le repère naturel varie de point en point. Quand on
passe du vecteur OM au vecteur OM ′ = OM + d OM , la base du repère naturel est changée en e i + d e i .
Nous trouvons que les d e i peuvent se mettre sous la forme :

d e i = Γi jk du k e j

où Γi j k est un coefficient, appelé symbole de Cristoffel de deuxième espèce.

Exemple en coordonnées sphériques :

r r r
Soit un vecteur OM = r sin θ cos ϕ i + r sin θ sin ϕ j + r cos θ k

Les vecteurs (non–unitaires) formant une base associés au point M sont :

r r r r
e1 = ∂ 1 M = sin θ cos ϕ i + sin θ sin ϕ j + cos θ k
r r r r
e 2 = ∂ 2 M = r cos θ cos ϕ i + r cos θ sin ϕ j − r sin θ k
r r r
e 3 = ∂ 3 M = − r sin θ sin ϕ i + r sin θ cos ϕ j

Nous avons donc :

g 11 = e1 ⋅ e1 = 1
g 22 = e 2 ⋅ e 2 = r 2
g 33 = e 3 ⋅ e 3 = r 2 sin 2 θ

donc le carré de l'intervalle entre deux points séparés d'une distance infinitésimale est :

ds 2 = dr 2 + r 2 dθ 2 + r 2 sin 2 θ dϕ 2

Nous trouvons ensuite :

1 1
d e1 = dθ e 2 + dϕ e 3
r r
1
d e 2 = − rdθ e1 + dr e 2 + cot an θ dϕ e 3
r
1 
d e 3 = − r sin 2 θ dϕ e1 − sin θ cos θ dϕ e 2 +  dr + cot an θ dθ  e 3
r 

Ainsi, par exemple :

3 1 3
Γ1 3 = Γ2 3 = cot an θ
r

On définit aussi un autre symbole de Cristoffel de première espèce par la relation :

Γkji = g jl Γk l i

32
Trouvons maintenant le lien qui existe entre les g ij et les symboles de Cristoffel. Nous avons :

g ij = e i ⋅ e j

donc :
dg ij = e i ⋅ d e j + e j ⋅ d e i

et d'après la formule de d e i :

dgij = ei ⋅ Γ j l k du k el + e j ⋅ Γi l k du k el

= g il Γ j l k du k + g jl Γi lk du k

(
= Γ jik + Γijk du k )
D'autre part, nous avons aussi :
∂g ij
dg ij = du k
∂u k

donc, par identification :

Γ jik + Γijk = ∂ k g ij

Par permutation circulaire, nous obtenons :

Γkji + Γ jki = ∂ i g jk
Γihj + Γkij = ∂ j g ki

Nous en déduisons donc que :

1
Γijk = (∂ k g ij + ∂ i g jk − ∂ j g ki )
2

d'où :
1 il
Γj i
k = g il Γ jlk =
2
(
g ∂ k g jl + ∂ j g lk − ∂ l g jk )

III. Dérivée covariante

r
Soit un vecteur v = v i e i . Cherchons comment varie ce vecteur. Sa différentielle totale s'écrit :

r
dv = e i ⋅ dv i + v i ⋅ d e i

Nous avons vu plus haut que :

d e i = Γi jk du k e j

Donc :

r
dv = ei ⋅ dvi + vi ⋅ Γi jk du k e j

(
= dv j + vi ⋅ Γi jk du k e j )
= (∂ v k
j
+ vi ⋅ Γi jk ) du k
ej

33
La quantité entre parenthèse de la dernière égalité s'appelle dérivée covariante. Dans les coordonnées
cartésiennes habituelles, les symboles de Cristoffel sont nuls et la dérivée covariante s'identifie à la dérivée
ordinaire ∂ k . La dérivée covariante sera notée maintenant ∇ k .

Pour un vecteur covariant, la dérivée covariante s'écrit :

∇ k v j = ∂ k v j − vi ⋅ Γ j i k

Nous pouvons étendre cette définition pour les tenseurs. Prenons par exemple le tenseur t ijk . Sa dérivée
covariante est :

∇ l tijk = ∂ l tijk + tijm Γl k m − tmj

k
Γl mi − tim
k
Γl mj

Dérivée covariante du tenseur métrique :

∇ k gij = ∂ k gij − g mj Γkmi − g imΓkmj

∇ k gij = ∂ k gij − Γkji − Γkij
1 1
∇ k gij = ∂ k gij −
2
( ) (
∂ i g kj + ∂ k g ji − ∂ j g ik − ∂ j g ki + ∂ k g ij − ∂ i g jk
2
)
∇ k gij = 0

Ceci constitue le théorème de Ricci.

Nous pouvons étendre toutes les définitions précédentes dans un espace comportant plus de trois dimensions.
Nous noterons alors les objets mathématiques avec des indices grecs.

IV. Tenseur de Riemann–Cristoffel

Nous définissons un tenseur d'ordre 4, décrivant la courbure d'un espace à quatre dimensions (indices grecs), une
σ
fois contravariant et trois fois covariant, Rµνρ . C'est le commutateur des dérivées covariantes d'un vecteur :

( ) σ
∇ ρ (∇ν ) vµ − ∇ν ∇ ρ vµ = − Rµνρ vσ

En passant les calculs, nous trouvons :

σ λ σ λ σ
Rµνρ = ∂σ Γµνσ − ∂ν Γµρσ + Γµν Γλρ − Γµρ Γλν

Nous voyons que :

σ
R µνρ σ
= − R µρν

Nous pouvons construire un tenseur de Riemann–Cristoffel complètement covariant :

β
R µνρσ = gνβ R µρσ

Nous avons ainsi :

1
Rµνρσ =
2
( )
∂ ρ ∂ µ gσν + ∂σ ∂ν g µρ − ∂ ρ ∂ν g µσ − ∂σ ∂ µ g ρν − Γρλν Γµ λ σ + Γσ λν Γµ λρ

De cette relation, nous en déduisons la première identité de Bianchi :

R µνρσ + R µρσν + R µσνρ = 0

34
 σ .
Nous obtenons la deuxième identité de Bianchi en considérant la dérivée covariante de R µνρ

λ λ λ
∇ α R µρσ = ∂ ρ ∂ α Γσ µ − ∂ σ ∂ α Γρ µ

Par permutation circulaire sur les indices ρ , σ , α et en sommant, nous trouvons :

λ
∇ α R µρσ λ
+ ∇ ρ R µσα λ
+ ∇ σ R µαρ =0

V. Tenseur de Ricci, tenseur d'Einstein et courbure scalaire

Le tenseur de Ricci est obtenu en contractant deux indices du tenseur de Riemann–Cristoffel :

λ λ λ ρ λ ρ λ
Rµν = Rµλν = ∂ν Γµλ − ∂ λ Γµν + Γµλ Γρν − Γµν Γρλ

La courbure scalaire R s'écrit :

R = R µµ = g µν R µν

Exemple dans le cas de la sphère de rayon a :

ds 2 = a 2 dθ 2 + a 2 sin 2 θ dϕ 2

où l'on a fixé r = a .

Les symboles de Cristoffel non nuls sont, avec x 1 = θ et x 2 = ϕ :

1
Γ2 2 = − sin θ cos ϕ
1 1
Γ1 2 = Γ2 1 = cotan θ

Les composantes non nulles du tenseur de Ricci sont :

2 2 2
R11 = ∂ 1 Γ1 2 + Γ1 2 Γ1 2 = −1
1 1 2
R 22 = −∂ 1 Γ2 2 + Γ2 2 Γ1 2 = − sin 2 θ

La courbure scalaire est donc :

2
R = g 11 R11 + g 22 R 22 = −
a2

Le tenseur d'Einstein se déduit de la deuxième identité de Bianchi. Cette dernière s'écrit :

λ
∇ α R µρσ λ
+ ∇ ρ R µσα λ
+ ∇ σ R µαρ =0

Faisons α = λ :

λ
∇ λ R µρσ λ
+ ∇ ρ R µσλ λ
+ ∇ σ R µλρ =0

donc

λ
∇ λ R µρσ − ∇ ρ R µσ + ∇ σ R µρ = 0

35
 ( )
or ∇ σ R µρ = ∇ σ g µν Rνρ = g µν ∇ σ Rνρ

donc, en multipliant par g µν :

νλ
∇ λ R ρσ − ∇ ρ Rσν + ∇ σ Rνρ = 0

Nous faisons ν = σ :
∇ λ Rρλ − ∇ ρ Rσσ + ∇σ Rσρ = 0

Nous changeons ensuite l'indice de sommation λ par σ :

2∇ σ R ρσ − ∇ ρ Rσσ = 0
⇔ 2∇ σ R ρσ − δ ρσ ∇ σ R = 0
 1 
⇔ ∇ σ  R ρσ − δ ρσ R  = 0
 2 

Nous obtenons donc un tenseur, appelé tenseur d'Einstein, dont la dérivée covariante est nulle :

1
G ρσ = R ρσ − δ ρσ R
2

c'est–à–dire, écrit sous sa forme entièrement covariante :

1
G µν = R µν − g µν R
2

VI. Équation des géodésiques dans un espace courbe

La métrique s'écrit :
(
ds 2 = g µν dx µ dxν ⇔ ds = g µν dx µ dxν )1 / 2
c'est–à–dire :
1/ 2
 dx µ dxν 
ds =  g µν ds
 ds ds 


B
La courbe géodésique est une courbe extrémale qui vérifie δ
∫ A
ds = 0 .

Cette relation équivaut à :

B
δ
∫ L(x , x& ) ds = 0
A
α α

1/ 2
 dx µ dxν  dxα
avec L (
xα )
, x& α =  g µν
 ds ds 
et x&α = .
ds


∂L d ∂L
Les équations de Lagrange s'écrivent : = , nous avons donc :
∂x α ds ∂x& α

∂L 1 −1 / 2 ∂g µν
= (g µν x& µ x& ν ) x& µ x&ν
∂x α 2 α ∂x

36
 ∂L 1
∂x&α
(
= g µν x& µ x&ν
2
∂
∂x&α
)
−1 / 2
(
g µν x& µ x&ν )
∂L 1 −1 / 2  µ ν 
(
= g µν x& µ x&ν  g µν ∂x& x&ν + g µν ∂x& x& µ 
)
∂x&α 2  ∂x&α
∂x&α 
 
∂L 1
∂x&α
(µ
= g µν x& x&
2
ν −1 / 2
) (
µ ν ν
g µν δα x& + g µν δα x& µ
)
∂L 1
∂x&α
(
= g µν x& µ x&ν
2
−1 / 2
) (
gαν x&ν + gαµ x& µ )
Par symétrie sur les indices µ et ν dans la première parenthèse, nous pouvons poser :

∂L
(
= g µν x& µ x& ν )−1 / 2 g αµ x& µ
∂x& α

donc :
d ∂L
(
= g µν x& µ x& ν )−1 / 2 dsd (g αµ x& µ )
ds ∂x& α

et ainsi, en combinant les deux relations :

d 1
ds
( )
gαµ x& µ = ∂α g µν x& µ x&ν
2
d  µ 1
 gαµ  x& + gαµ &x&µ = ∂α g µν x& µ x&ν
 ds  2
 ∂  1
 λ gαµ  x& µ x& λ + gαµ &x&µ = ∂α g µν x& µ x&ν
 ∂x  2

λ est un indice muet. Faisons λ = ν :

1 1
gαµ &x&µ =
2
(
∂α g µν x& µ x&ν − ∂ν gαµ x& µ x&ν = ∂α g µν − 2∂ν gαµ x& µ x&ν
2
)

Par symétrie sur les indices µ et ν dans le membre x& µ x&ν , nous avons :

(
2∂ν g αµ x& µ x& ν = ∂ν g αµ + ∂ µ g αν x& µ x& ν )
Nous obtenons donc :
1
gαµ &x&µ =
2
(
∂α g µν − ∂ν gαµ − ∂ µ gαν x& µ x&ν )

ou :
1 βα
g αβ &x&α +
2
(
g ∂ν gαµ + ∂ µ gαν − ∂α g µν x& µ x&ν = 0 )
1 βα
&x&β +
2
(
g ∂ µ gαν + ∂ν gαµ − ∂α g µν x& µ x&ν = 0 )
&x&β + Γµβν x& µ x&ν = 0

d 2 xβ dx µ dxν
2
+ Γµβν =0
ds ds ds

C'est l'équation des géodésiques.

37
VII. Tenseur d'énergie–impulsion
Dans un repère lié à une particule, l'énergie de cette particule est E = m c 2 . La densité d'énergie correspondante
est ρ 0 c 2 , que nous définirons égale à la composante T 44 du tenseur T µν .

Cherchons maintenant l'expression de ce tenseur dans n'importe quel repère. La loi de transformation tensorielle
s'écrit :

T µ 'ν ' = Λµµ'Λνν' T µν

∂x µ ' ∂xν ' µν
T µ 'ν ' = T
∂x µ ∂xν
∂x µ ' ∂xν '
T µ 'ν ' = T 44
∂x 4 ∂x 4
dx µ ' dxν '
T µ 'ν ' = ρ 0c 2
dx 4 dx 4

dx µ  γ r 
d'où, en posant u µ = =  v , γ  la quadrivitesse en unités c,
cdt  c 

T µ 'ν ' = u µ 'uν ' ρ 0c 2

Ainsi, le tenseur d'énergie–impulsion T µν a le caractère tensoriel et est symétrique.

En posant ρ = γ 2 ρ 0 , nous avons :

 ρ u1u1 ρ u1u 2 ρ u1u 3 ρ c u1 

 2 1 
ρ u u ρ u 2u 2 ρ u 2u 3 ρ c u 2 
T µν =  3 1
ρ u u ρ u 3u 2 ρ u 3u 3 ρ c u 3 
 
 ρ c u1 ρ c u 2 ρ c u 3 ρ c 2 

Interprétation physique de Tij :

Tij est homogène à une force par unité de surface. Lorsque i = j , il s'agit de la pression p. Lorsque i ≠ j , il s'agit
de la viscosité.

 p ρ u1u 2 ρ u1u 3 ρ cu1 

 2 1 
ρu u p ρ u 2u 3 ρ cu 2 
T µν =  3 1
ρ u u ρ u 3u 2 p ρ cu 3 
 
 ρ cu1 ρ cu 2 ρ cu 3
ρ c 2 

dans le cas des fluides parfaits, il n'y a pas de transferts d'énergie donc T4i = Ti 4 = 0 , et pas de viscosité. Le
tenseur d'énergie impulsion se résume à :

p 0 0 0 
 
0 p 0 0 
T µν = 
0 0 p 0 
 
 0 0 0 ρ c 2 

38
Établissons sa formule générale :

∂x µ ' ∂xν ' µν

T µ 'ν ' = T
∂x µ ∂xν
∂x µ ' ∂xν '
T µ 'ν ' = u′ µ u′ν ρ c 2 + p
∂xi ∂xi
or

∂x µ ' ∂xν ' µν ∂x µ ' ∂xν ' ∂x µ ' ∂xν '

η µ 'ν ' = η = 4 − i
∂x µ ∂xν ∂x ∂x 4 ∂x ∂xi

d'où

∂x µ ' ∂xν '

= u′µ u′ν − η µ 'ν '
∂xi ∂xi

Donc,

(
T µ 'ν ' = u′ µ u′ν ρ c 2 + u′µ u′ν − η µ 'ν ' p )
T µ 'ν '
( 2
= ρ c + p u u − pη) µ' ν ' µ 'ν '

Dans un espace décrit par la métrique g µν , le tenseur d'énergie–impulsion s'écrit :

( )
T µν = ρ c 2 + p u µ uν − pg µν

Dérivée covariante dans le cas général :

Nous avons : T µν = ρ u µ uν donc :

∂T µν
=0
∂x µ

Ainsi, dans un système où les coordonnées curvilignes sont utilisées :

∇ µ T µν = 0

Nous voyons donc que T µν a les mêmes propriétés tensorielles que le tenseur d'Einstein G µν .

39
 Chapitre 2 : La gravitation newtonienne

I. La loi de la gravitation universelle

Formulée dans les Principia Mathematica en 1687, la loi de la gravitation universelle s'écrit :

r GMm
F=− er
r2

Elle est une généralisation des lois de Képler du mouvement des planètes. Elle permit de prévoir le mouvement
parabolique ou hyperbolique de certaines comètes et aussi d'établir le mouvement de mobiles soumis à la
pesanteur à la surface de la Terre.
Elle a également permis de découvrir la planète Neptune en observant le mouvement d'Uranus, mais elle a été
incapable d'expliquer l'avance supplémentaire du périhélie de Mercure de 43" par siècle.

II. L'équation de Poisson

r GM r GM
En posant g = − e r , nous voyons que g dérive du potentiel U = − :
2 r
r

r
g = −grad U

r
Nous en déduisons : div g = − ∆U où ∆ désigne le laplacien.
r
Calculons le flux de g :

r r GM
dφ = g ⋅ dS = − 2 dS cosθ = −GM dΩ
r

cos θ
dΩ représente l'angle solide élémentaire défini par : dΩ = dS .
r2

Donc :

φ = −4π GM

Par ailleurs, le théorème de Gauss–Ostrogradski nous permet de dire :

r
∫∫ g ⋅ dS = ∫∫∫ div g ⋅ dV
r r
φ=
S V

De plus, nous pouvons écrire :

M =
∫∫∫ ρ dV
V

Donc :

∫∫∫ div g ⋅ dV = −4π G ∫∫∫ ρ dV

r
φ=
V V

40
c'est–à–dire :
r
div g = −∆U = −4π Gρ

∆U = 4π Gρ

C'est l'équation de Poisson. Nous pouvons en déduire l'équation de Laplace ∆U = 0 , valable par exemple pour le
mouvement des planètes.

III. Équation de mouvement d'une particule dans un champ gravitationnel

newtonien
Une particule de masse m est soumise à la loi de la gravitation universelle :

GMm
F =− er
r2

En appliquant la troisième loi de Newton, nous obtenons :

r r
F = ma

avec :

r d 2 OM GM
a= =− er
2
dt r2

or OM = r e r donc :
r d OM
v= = r& e r + rθ& eθ
dt

r
r dv
a= ( ) (
= &r& − rθ& 2 e r + rθ&& + 2 r&θ& eθ )
dt

ce qui, en vertu de la loi de Newton, revient à écrire :

rθ&& + 2 r&θ& = 0
1 2 &&
r
(r θ + 2 rr&θ& = 0 )
1 d 2&
r dt
( )
r θ =0

c'est–à–dire :

r 2θ& = C

La constante C est appelée constante des aires.

Nous avons ensuite :

GM
r
a=− (
e r = &r& − rθ& 2 e r )
r2

GM C2
− = &r& −
2
r r3

41
 1
En posant u = ,
r

du 1 dr
=−
dt u 2 dt

Nous en déduisons :

d 2r d  1 du  d  1 du dθ 
2
= &r& =  −  =  − 
dt dt  u 2 dt  dt  u 2 dθ dt 

dθ d  du  dθ d 2 u d 2u
&r& = −C  = −C = −C 2 u 2
dt dθ  dθ  dt dθ 2 dθ 2

Nous avons ainsi :

d 2u
− GM u 2 = −C 2 u 2 − C 2u 3
dθ 2

En divisant par − C 2 u 2 , nous obtenons :

d 2u GM
+u =
2
dθ C2

C'est l'équation de mouvement d'une particule dans un champ gravitationnel créé par la masse M.

La solution de cette équation est :

1 1 + e cos (θ − θ 0 )
u= =
r p

C2
p= est appelé paramètre de la conique et e est son excentricité.
GM

Cas de la déviation d'un photon :

Si on considère en mécanique classique qu'un photon possède une masse mγ , alors sa trajectoire est solution de
L
l'équation de mouvement établie ci–dessus. On peut montrer que, dans ce cas : C = , où L est le moment
mγ
r r
cinétique associé au photon. Par définition, L = AB ∧ mγ v , donc L2 = (b m γ c )2 = constante. Ainsi, l'équation de
mouvement se réécrit :

d 2u GM
+u =
2
dθ b2c2

1 GM mγ2
dont la solution est : u (θ ) = = (1 + e cos θ ) .
r (θ ) L2

Pour un photon venant de l'infini, on a :

r → ∞, u → 0

42
c'est–à–dire :

1 + e cos θ = 0

On a ainsi :

1 δ
cos θ = − = sin
e 2

où δ est la déviation recherchée.

Lorsque θ → 0, r → b donc :

L2 1
b=
GM mγ2 1 + e

L2 bc2 bc2
⇒e= −1 = −1 ≈
GM mγ2 b GM GM

Donc

δ 1 GM δ
sin = = ≈
2
2 e bc 2

c'est–à–dire :

2GM
δ =
bc2

La théorie newtonienne de la gravitation prévoit une déviation du photon dans le cas où il posséderait une masse.
Dans le cas où un photon frôlerait le bord du Soleil ( b = 6.96 . 10 8 m , M = 1.989 . 10 30 kg ), nous obtenons une
déviation de 0.875".

43
 Chapitre 3 : Bases de la relativité générale

I. Principes de la relativité générale

1. Principe d'équivalence

En relativité restreinte, les transformations de Lorentz ne sont valables que pour des référentiels galiléens. La
question qui se pose est : pourquoi existe–t–il des repères privilégiés? Einstein a répondu qu'il n'existe pas de
référentiels privilégiés et a proposé le principe d'équivalence comme base d'une nouvelle théorie. Le principe
d'équivalence est un principe selon lequel, localement, il est impossible de distinguer un système soumis à un
champ gravitationnel de celui soumis à une accélération.

Considérons un laboratoire en chute libre dans un champ gravitationnel. Les corps lâchés sans vitesse initiale se
trouvant dans le laboratoire flottent librement s'ils n'ont pas de vitesse initiale, et ont une trajectoire rectiligne
uniforme s'ils ont une vitesse initiale constante. Le laboratoire apparaît pour ceux qui sont à l'intérieur comme un
laboratoire libre de toute force. C'est ce que l'on observe dans les stations spatiales. Pour un observateur
extérieur, le laboratoire est soumis à l'attraction gravitationnelle. Si maintenant le laboratoire est infiniment loin
de toute masse, les expériences menées à l'intérieur mènent aux mêmes résultats que ci-dessus.

Si ce même laboratoire est à la surface d'un corps (comme la Terre), les objets lâchés sans vitesse initiale chutent
avec un mouvement uniformément accéléré vers le sol, tandis que les objets ayant une vitesse initiale suivent une
trajectoire dont la forme est un arc de parabole. Si le laboratoire, loin de toute masse, est muni d'un moteur sous
le plancher permettant de l'accélérer de manière uniforme, les objets auront un mouvement identique que s'ils
étaient soumis à la pesanteur.

Ces expériences permettent de conclure que le champ gravitationnel peut être interprété comme un champ
accéléré. Cela est équivalent à dire que la masse inertielle (participant au mouvement) est identique à la masse
grave (participant aux interactions gravitationnelles).

2. La géométrie de l'espace–temps

L'idée d'Einstein est de dire que l'espace–temps dans lequel nous vivons n'est pas celui de la relativité restreinte
mais un espace–temps déformé par la présence de matière. Nous pouvons alors formuler que la courbure de
l'espace–temps est proportionnelle à la matière et l'énergie présentes, ce qui s'écrit :

courbure ∝ matière

Pour formuler mathématiquement cette relation, nous utilisons le formalisme tensoriel. Nous connaissons le
tenseur d'énergie–impulsion T µν , dont la dérivée covariante est nulle, qui représente la matière. Nous devons lui
associer un tenseur caractérisant la géométrie et ayant les mêmes propriétés que lui. Nous avons vu que ce
tenseur est le tenseur d'Einstein G µν . Nous pouvons donc écrire la relation suivante :

1
G µν = R µν − g µν R = χT µν
2

où χ est un coefficient de proportionnalité (constante de couplage).

Nous pouvons cependant remarquer qu'il existe une équation plus générale :

1
R µν − g µν R + Λg µν = χT µν
2

où Λ est appelée constante cosmologique.

Cette nouvelle relation est également valable car la dérivée covariante du premier membre est nulle.

44
3. Détermination de la constante χ

Le but est de rapprocher l'équation tensorielle d'Einstein à l'équation de Poisson : ∆U = 4π Gρ .

L'accélération dans le champ de pesanteur décrit par la mécanique newtonienne est :

r
r d 2x
g = 2 = −grad U
dt

GM
avec U = − , lorsque la source du champ gravitationnel possède une masse M.
r

Dans le cas d'un champ gravitationnel faible, on pose g µν = η µν + h µν , où η µν est la métrique de la relativité
restreinte. On a h µν << 1 et h µν (∞ ) → 0 (on retrouve la métrique de la relativité restreinte à une distance infinie
de la source).

D'après l'équation des géodésiques établie plus haut, la seule composante significative est celle contenant
µ = ν = 4 (car le tenseur d'énergie–impulsion se résume à sa composante T44 = ρ c 2 ) :

2
d 2xβ β  dx 
4
+ Γ4 4   =0
ds 2  ds 
 

1 1
( )
avec Γ4 β4 = − g βν ∂ν g 44 = − η βν + h βν ∂ν (η44 + h44 ) .
2 2

1 1
Donc Γ4 β4 ≈ − η βν ∂ν h44 = − ∂ β h44 .
2 2

Nous avons donc :

2
d 2 xβ 1 β  dx 4 
− ∂ h44 =0
ds 2 2  ds 
 

Faisons β = 4 :

2
d 2 x4 1  dx 4 
− ∂ 4 h44   =0
ds 2 2  ds 
 

d 2x4 d 2t
or le champ gravitationnel newtonien est stationnaire donc ∂ 4 h44 = 0 , ce qui signifie que =c =0
ds 2 ds 2
dt
c'est–à–dire = cst .
ds

2 2
d 2 xβ 1 4  dx 4   dt 
En divisant − ∂ h44 = 0 par   , nous trouvons :
ds 2 2  ds   ds 
 

d 2xβ 1 4
2
− ∂ h44 c 2 = 0
dt 2

45
soit vectoriellement, pour β ≤ 3 :

r
d 2x c2 r
2
+ grad h44 = 0
dt 2
r
d 2x c2
⇔ 2
=− grad h44 = −grad U
dt 2

c2
Donc U = h44 + cste , or h µν (∞ ) → 0 et U (∞ ) → 0 donc cste = 0, ce qui revient à écrire :
2

2U 2GM
h44 = =
c2 c2r

Écrivons l'équation d'Einstein sous une autre forme ( Λ = 0 ). Nous avons :

1
R µν − g µν R = χT µν
2

Multiplions membre à membre par g µν , sachant que g µν g µν = 4 :

R − 2 R = χg µν T µν
R = − χT

En injectant ce résultat dans l'équation d'Einstein, nous obtenons alors :

 1 
R µν = χ  T µν − g µν T 
 2 

Nous avons alors :

 1  1
R44 = χ  T44 − g 44T  = χρ c 2
 2  2

avec R44 = ∂ 4Γ4β β − ∂ β Γ4β 4 + Γ4 λ4 Γ4β λ − Γ4 λ4 Γββ 4 . En faisant l'hypothèse g µν = η µν + h µν , il ne reste que :

1 β
R44 = −∂ β Γ4β 4 = ∂ ∂ β h44
2

On obtient ainsi :

∂ β ∂ β h44 = χρ c 2

h44 = χρ c 2

2U
Or on a vu que h44 = , donc :
c2
2

h44 =
U
c2

À la limite statique,
U = − ∆U , c'est–à–dire :

2

h44 = − ∆U = χρ c 2
c2

46
Or l'équation de Poisson s'écrit ∆U = 4π Gρ , on a finalement :

8π G
χ=−
c4

Ainsi, l'équation d'Einstein s'écrit :

1 8π G
G µν = R µν − g µν R + Λg µν = − T µν
2 c4

Remarquons la petitesse du coefficient de proportionnalité χ ( ≈ −2.10 −43 ) qui met bien en évidence que la
gravitation est une interaction très faible.

II. Métrique de Schwarzschild

1. Établissement de la métrique de Schwarzschild

Schwarzschild a proposé dès 1916 une solution des équations d'Einstein en faisant les hypothèses suivantes :

– la constante cosmologique Λ est nulle.

– on considère le cas extérieur, c'est–à–dire en dehors de la source du champ gravitationnel ( ρ = 0 ).
– le champ gravitationnel est à symétrie sphérique, faible et stationnaire.
– à une distance infinie de la source, la métrique de l'espace–temps est celle de la relativité restreinte.

D'après ce qui précède, la métrique s'écrit sous la forme :

ds 2 = c 2eν dt 2 − eλ dr 2 − r 2 dθ 2 − r 2 sin 2 θ dϕ 2

avec λ (∞ ) → 0 et ν (∞ ) → 0 pour satisfaire la dernière hypothèse.

dν
En posant r = x 1 , θ = x 2 , ϕ = x 3 , ct = x 4 , les symboles de Cristoffel non nuls sont, avec ν ' = :
dr

1 ν' 1 ν' 1 λ' 1 3 1

Γ1 4 = Γ4 4 = eν -λ Γ1 1 = Γ2 2 = − r e - λ Γ1 3 =
2 2 2 r

1 2 1 2 3
Γ3 3 = − r sin 2 θ e -λ Γ1 2 = Γ3 3 = − sin θ cos θ Γ2 3 = cotan θ
r

Dans l'hypothèse du cas extérieur, nous avons T µν = 0 (donc T = 0 ), c'est–à–dire :

R µν = 0
soit, pour chaque composante :

ν ′′ λ ′ν ′ ν ′ 2 λ′  r 
R11 = − + − R 22 = e -λ 1 + (ν ′ − λ ′ ) − 1
2 4 4 r  2 
 ν ′′ λ ′ν ′ ν ′ 2 ν ′ 
R33 = R 22 sin 2 θ R 44 = − eν -λ  − + + 
 2 4 4 r 


Nous avons notamment R11 = 0 et R 44 = 0 donc :

47
 ν ′′ λ ′ν ′ ν ′ 2 λ′ 
− + − = 0
2 4 4 r 

 ⇒ ν ′ = − λ ′ soit ν = − λ + cste .

ν ′′ λ ′ν ′ ν ′ 2 ν′
− + + = 0
2 4 4 r 

Mais la condition d'avoir la métrique de la relativité restreinte à l'infini impose cste = 0.

Comme ν = −λ , nous trouvons pour R 22 :

R 22 = eν (1 + ν ′ r ) − 1 = 0

Nous posons y = eν pour obtenir :

r y′ + y = 1

y est de la forme y = a + b r m ⇒ y ′ = b m r m −1 .

D'où :

ry ′ + y = a + b (1 + m ) r m = 1

Comme nous devons avoir asymptotiquement la métrique de la relativité restreinte, nous avons a = 1 , donc
m = −1 car ni b ni r ne sont nuls.

Nous avons finalement :

b
eν = 1 +
r

La métrique de Schwarzschild s'écrit donc :

 b dr 2
ds 2 = c 2 1 +  dt 2 − − r 2 dθ 2 − r 2 sin 2 θ dϕ 2
 r b
1+
r

2. Détermination de b

Le terme temporel de la métrique de Schwarzschild s'écrit en g 44 = η 44 + h44 = 1 + h44 . Or, nous avons déjà déterminé pour
un champ faible la valeur de h44 . Nous obtenons ainsi :

2U 2GM rS
h 44 = =− =−
2 2 r
c c r

où rS est appelé rayon de Schwarzschild.

La métrique de Schwarzschild s'écrit finalement :

 r  dr 2
ds 2 = c 2 1 − S  dt 2 − − r 2 dθ 2 − r 2 sin 2 θ dϕ 2
 r  rS
1−
r

48
Il existe une métrique de Schwarzschild plus générale tenant compte de la constante cosmologique :

 r Λr 2  2 dr 2
ds 2 = c 2 1 − S +  dt − − r 2 dθ 2 − r 2 sin 2 θ dϕ 2
 r 3  rS Λr 2
  1− +
r 3

Comme nous le verrons par la suite, la valeur numérique de la constante cosmologique est très faible, surtout pour les cas
dans le système solaire.

3. Mouvement d'une particule dans un espace décrit par la métrique de Schwarzschild

Nous rappelons tout d'abord que l'équation de mouvement dans un champ gravitationnel décrit par la mécanique
de Newton est :

d 2u GM
+u =
2
dθ C2

Pour établir l'équation de mouvement dans le cadre de la relativité générale, nous avons besoin de :

d 2xβ β dx µ dxν
+ Γµ ν =0
2 ds ds
ds
et

 r  dr 2
ds 2 = c 2 1 − S  dt 2 − − r 2 dθ 2 − r 2 sin 2 θ dϕ 2
 r  r S
1−
r

rS
Cette dernière relation s'écrit, en posant eν = e −λ = 1 − :
r

2 2 2 2
 dt   dr   dθ   dϕ 
1 − c 2 eν   + e λ   + r 2   + r 2 sin 2 θ   =0
 ds   ds   ds   ds 

Les équations des géodésiques s'écrivent, par exemple pour la variable r :

2 2 2 2
d 2r 1  dr  1  dθ  1  dϕ  1  dt 
+ Γ1 1   + Γ2 2   + Γ3 3   + Γ4 4   = 0
ds 2  ds   ds   ds   ds 

dν
soit, en posant ν ′ = :
dr

2 2 2 2
d 2r λ ′  dr   dθ   dϕ  ν′  dt 
+   − r e -λ   − r sin θ e -λ   + eν -λ   = 0
ds 2 2  ds   ds   ds  2  ds 

2
d 2θ 2 dr dθ  dϕ 
+ − sin θ cos θ   =0
ds 2 r ds ds  ds 

d 2ϕ 2 dr dϕ dθ dϕ
+ + 2 cotan θ =0
2 r ds ds ds ds
ds

d 2t ν ′ dr dt
+2 =0
ds 2 2 ds ds

49
 π
Le champ est à symétrie sphérique donc nous pouvons choisir θ = ⇒ dθ = 0 :
2
2 2 2
d 2r λ ′  dr   dϕ  ν′  dt 
+   − r e -λ   + eν -λ   = 0 (1)
ds 2 2  ds   ds  2  ds 

d 2ϕ 2 dr dϕ
+ =0 ( 2)
2 r ds ds
ds

d 2t ν ′ dr dt
+2 =0 (3)
2 2 ds ds
ds

Ainsi, la forme modifiée de la métrique de Schwarzschild se réécrit :

2 2 2
 dt   dr   dϕ 
1 − c 2 eν   + e λ   + r 2   =0 ( 4)
 ds   ds   ds 

dϕ dϕ
D'après cette dernière relation, nous voyons que est une fonction de r : = u (r ) donc (2) s'écrit :
ds ds

d u (r )
2 dr
+ u (r ) = 0
ds r ds
d u (r ) dr 2 dr
⇒ + u (r ) = 0
dr ds r ds
d u (r ) dr
⇒ = −2
u (r ) r

ce qui entraîne :

A
u (r ) =
r2

où A est une constante.

dt r
De même, nous voyons que = u (ν ) = k e −ν . Donc en multipliant (4) par eν = 1 − S , nous avons :
ds r

2 2
 rS   dr   dϕ   rS 
1 −  +   + r 2   1 −  − c 2 k 2 = 0
 r   ds   ds   r 

r   dϕ  
2 2 2
 dr   dϕ 
⇔   + r2  − S 1 + r 2   = c 2k 2 −1
 ds   ds  r   ds  

dϕ A
En utilisant le fait que = , nous obtenons :
ds r2

2
A 2  dr  A 2 rS  A 2 
  + − 1+ = c 2k 2 −1
4  dϕ  
r   r 2 r  r 2 

50
 1
En posant u = et en divisant par A 2 :
r

2 2 2
 du  r
  + u 2 − S u − rS u 3 = c k − 1
 dϕ  A2 A2
 

Nous dérivons par rapport à ϕ :

du d 2 u du rS du du
2 + 2u − − 3r S u 2 =0
dϕ dϕ 2 dϕ A 2 dϕ dϕ
du  d 2 u r 3r 
⇒ +u − S − S u2  = 0

dϕ  dϕ 2
2A 2 2 


du
La première solution = 0 correspond à une trajectoire circulaire, résultant de la symétrie sphérique.
dϕ

La deuxième solution s'écrit :

d 2u GM 3GM
+u = + u2
2 2 2
dϕ c A c2

3GM
En remarquant que u 2 est faible devant l'unité et en comparant avec l'équation obtenue en champ
c2
newtonien, nous avons finalement :

d 2u GM 3GM
+u = + u2
2 2
dϕ C c2

avec C = cA la constante des aires (obtenue en mécanique classique).

Nous avons donc l'équation de mouvement d'une particule dans un champ gravitationnel décrit par la métrique
de Schwarzschild :

d 2u 1 3GM 2
+u = + 2 u
dϕ 2 p c

4. Application à l'avance du périhélie d'une planète

Dans le système solaire, les perturbations planétaires font avancer le périhélie des planètes. Pour Mercure,
l'avance totale observée est de 574" par siècle. Les perturbations des planètes calculées avec l'équation de la
d 2u GM
mécanique classique 2
+u = donnent une avance du périhélie de 531" par siècle. Il y a un écart résiduel
dϕ C2
de 43". Voyons si le terme supplémentaire de la relativité générale permet d'expliquer cette avance
supplémentaire.

1 + e cos (ϕ − ϕ 0 )
La solution u 0 = de l'équation classique est une solution approchée de l'équation relativiste
p
d 2u 1 3GM 2
+u = + 2 u .
dϕ 2 p c

51
Nous obtenons donc :

d 2u 1 3GM
+u = + (1 + e cos (ϕ − ϕ 0 ))2
dϕ 2 p c2 p2

3GM 1
En négligeant le terme en e 2 (qui n'a pas d'effet cumulatif) et remarquant que << , nous avons :
c2 p p

d 2u 1 6GM
+u = + e cos (ϕ − ϕ 0 )
2
dϕ p c2 p2

La solution de cette équation est celle d'une ellipse perturbée :

u = A cos (α (ϕ − ϕ 0 )) + B α ≈1
du
u′ = = − A α sin (α (ϕ − ϕ 0 ))
dϕ
u ′′ = − Aα 2 cos (α (ϕ − ϕ 0 ))

Ainsi,

( )
u ′′ + u = A 1 − α 2 cos (α (ϕ − ϕ 0 )) + B

On identifie les constantes A et B en faisant α = 1 dans u = A cos (α (ϕ − ϕ 0 )) + B et en utilisant la solution

classique :

e 1
A= , B=
p p

On pose α = 1 + ε avec ε << 1 (ε peut être positif ou négatif). En négligeant le terme en ε 2 , on a :

1 − 2eε cos ((1 + ε )(ϕ − ϕ0 ))

u′′ + u =
p

Nous avons donc :

−2eε 6GM
= 2 2e
p p c
3GM
⇔ε =−
p c2

Le périhélie se produit lorsque r est minimal, donc quand u est maximal, c'est–à–dire quand (1 + ε )ϕ = 2π n , où n
est un entier positif. Ainsi,

ϕ = 2π n (1 + ε ) −1 ≈ 2π n (1 − ε )

Pour deux périhélies successifs, n = 1 ,

ϕ = 2π (1 − ε ) = 2π − 2π ε

La différence, ∆ϕ = −2π ε , correspond à une avance supplémentaire (car ε < 0 ) du périhélie d'une valeur :

6π GM 6π GM
∆ϕ = =
pc 2
(
c a 1− e2
2
)

52
Dans le cas de Mercure, p = 5.546.1010 m . Avec M = 1.989.10 30 kg , nous trouvons :

∆ϕ = 5.019.10 −7 rad

La période de révolution de Mercure étant de 87.9 jours, nous avons pour un siècle :

∆ϕ × 100 × 365.25
δ = ≈ 2.086.10 − 4 rad
87.9

c'est–à–dire δ = 43.02" par siècle, ce qui est conforme avec ce qui est observé.

5. Déviation des rayons lumineux en relativité générale

2GM
Rappelons que si on attribue au photon une masse mγ , on trouve une déviation δ = par la gravitation
bc2
newtonienne.

Dans un espace décrit par la métrique de Schwarzschild, nous avions trouvé :

d 2u GM 3GM
2
+u = 2 2
+ u2
dϕ A c c2

dϕ
avec r 2 = A . La géodésique d'un photon est ds = 0 , donc A devient infini, ce qui revient à écrire :
ds

d 2u 3GM
+u = u2
2
dϕ c2

cos ϕ
On résout tout d'abord u ′′ + u = 0 dont la solution est : u 0 = , solution que l'on porte dans l'équation de
b
mouvement ci–dessus.

3GM
u ′′ + u = cos 2 ϕ
b 2c 2

Une solution u1 est de la forme :

u1 = A cos 2 ϕ + B

donc

u1′ = −2 A sin ϕ cos ϕ

u1′′ = −2 A cos 2ϕ = −4 A cos 2 ϕ + 2 A

u1′′ + u1 = −3 A cos 2 ϕ + 2 A + B

En identifiant, on a :

 GM

− 3 A =
3GM A = − 2 2
  b c
 b2c2 ⇔ 
2 A + B = 0  B = −2 A = 2GM
  b2c2

53
La solution générale est u = u 0 + u1 :

cos ϕ GM 2GM
u= − cos 2 ϕ +
2 2
b b c b2c2

Le rayon lumineux d'une étoile semble provenir d'une distance infinie, donc, quand r → +∞, u → 0 :

cos ϕ GM 2GM
− cos 2 ϕ + =0
b b 2c 2 b 2c 2
GM 2GM
⇒ cos ϕ − cos 2 ϕ + =0
2
bc bc 2

α α α α  α2
La déviation α s'exprime par sin = cos ϕ . Comme α est petit, sin ≈ et sin 2  ≈
 ≈ 0 . Il reste donc :
2 2 2 2 4

α 2GM
=−
2 bc 2

c'est–à–dire, en valeur absolue :

4GM
α=
bc 2

Nous voyons donc que la relativité générale prédit une déviation double de celle prévue par la mécanique
classique ( α = 2δ ). Cette déviation, mesurée pour la première fois en 1919 par Sir Arthur Eddington lors d'une
éclipse totale de Soleil, fournit une déviation de 1.98", comparable à la valeur issue de la relativité générale. Des
mesures plus récentes sont en accord avec la Relativité Générale.

6. Trous noirs

Le concept de trous noirs fut imaginé bien avant la Relativité Générale. Il fut envisagé au XVIIIe siècle, pour des
corps dont la vitesse de libération dépasse celle de la lumière. Considérons un photon se déplaçant radialement
dans la métrique de Schwarzschild :

 r  dr 2
ds 2 = c 2 1 − S dt 2 − =0
 r  r
1− S
r
dr
c dt =
r
1− S
r
soit :

dr  r 
= c 1 − S 
dt  r 

Si le rayon de Schwarzschild devient supérieur à r, la lumière ne peut s'extraire de l'attraction du corps, ce

dernier est un trou noir. Si le rayon du Soleil avait une valeur de 2950m, celui-ci serait un trou noir. La Terre
serait un trou noir si elle avait un rayon de 8.8mm.

54