ELECTRODYNAMIQUE
ELECTRODYNAMIQUE
ELECTRODYNAMIQUE
ÉLECTRODYNAMIQUE
ET RELATIVITÉ
E D
ε1 > ε0 ε1 < ε0
NOTES DE COURS
(PHQ-420)
Université de Sherbrooke
Faculté des Sciences
Avril 1996
ÉLECTRODYNAMIQUE
ET RELATIVITÉ
NOTES DE COURS
(PHQ-420)
par
David Sénéchal
Professeur
Département de physique
c 1996, David Sénéchal, Faculté des Sciences, Université de Sherbrooke.
Tous droits réservés.
Table des Matières
Préface . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
1. Électrostatique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
1.1. Lois de Coulomb et de Gauss . . . . . . . . . . . . . . . . . . . . . . . . . . 4
1. Charge électrique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
2. Loi de Coulomb . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
3. Loi de Gauss . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
1.2. Potentiel électrique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
1. Définition . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
2. Énergie potentielle d’une charge dans un champ . . . . . . . . . . . . . . . . 9
3. Énergie potentielle d’une distribution de charge . . . . . . . . . . . . . . . . 11
1.3. Électrostatique et conducteurs . . . . . . . . . . . . . . . . . . . . . . . . . . 12
1. Équation de continuité . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
2. Loi d’Ohm . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
3. Distribution des charges dans un conducteur . . . . . . . . . . . . . . . . . 14
4. Capacité . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
1.4. Multipôles . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
1. Moment Dipolaire Électrique . . . . . . . . . . . . . . . . . . . . . . . . 17
2. Moment Quadrupolaire Électrique . . . . . . . . . . . . . . . . . . . . . . 18
3. Théorème d’addition et forme générale . . . . . . . . . . . . . . . . . . . . 19
4. Énergie d’une distribution de charge dans un potentiel externe . . . . . . . . . 20
1.5. Méthodes de solution de l’équation de Laplace . . . . . . . . . . . . . . . . . . . 22
1. Propriétés des fonctions harmoniques . . . . . . . . . . . . . . . . . . . . . 22
2. Séparation des variables (général) . . . . . . . . . . . . . . . . . . . . . . 23
3. Séparation des variables en coordonnées cartésiennes . . . . . . . . . . . . . . 23
4. Séparation des variables en coordonnées sphériques . . . . . . . . . . . . . . 25
5. Séparation des variables en coordonnées cylindriques . . . . . . . . . . . . . . 28
1.6. Diélectriques . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29
1. Origine microscopique de la polarisation . . . . . . . . . . . . . . . . . . . 29
2. Densité de polarisation . . . . . . . . . . . . . . . . . . . . . . . . . . . 32
3. Induction électrique . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
4. Constante diélectrique . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
5. Conditions de continuité pour D et E . . . . . . . . . . . . . . . . . . . . 34
6. Énergie électrostatique et calcul des forces électriques . . . . . . . . . . . . . 36
1.A. Notation indicielle . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 39
1.B. Polynômes de Legendre . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40
Problèmes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 43
iii
2. Magnétostatique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
2.1. Loi d’Ampère . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
1. Loi de Biot et Savart . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
2. Potentiel vecteur . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 51
3. Loi d’Ampère . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 51
2.2. Mouvement dans des champs électrique et magnétique . . . . . . . . . . . . . . . 53
1. Force de Lorentz . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53
2. Formulations Lagrangienne et Hamiltonienne . . . . . . . . . . . . . . . . . 53
2.3. Dipôles magnétiques . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 54
1. Moment dipolaire . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 54
2. Moment magnétique et moment cinétique . . . . . . . . . . . . . . . . . . . 56
3. Force et couple sur un moment dipolaire . . . . . . . . . . . . . . . . . . . 58
4. Précession de Larmor . . . . . . . . . . . . . . . . . . . . . . . . . . . 59
2.4. Magnétisme macroscopique . . . . . . . . . . . . . . . . . . . . . . . . . . . 61
1. Aimantation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 61
2. Diamagnétisme . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 62
3. Paramagnétisme . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 63
4. Ferromagnétisme . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 64
5. Conditions aux limites des champs B et H . . . . . . . . . . . . . . . . . . 65
6. Potentiel scalaire magnétique . . . . . . . . . . . . . . . . . . . . . . . . 66
Problèmes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 68
3. Électrodynamique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 70
3.1. Induction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 70
1. Loi de Faraday . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 70
2. Inductance . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 71
3.2. Équations de Maxwell . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 72
1. Courant de Déplacement . . . . . . . . . . . . . . . . . . . . . . . . . . 72
2. Équations de Maxwell . . . . . . . . . . . . . . . . . . . . . . . . . . . 73
3.3. Potentiels électromagnétiques . . . . . . . . . . . . . . . . . . . . . . . . . . 74
1. Transformations de jauge . . . . . . . . . . . . . . . . . . . . . . . . . . 74
2. Couplage minimal . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 75
3. Équation du mouvement pour les potentiels . . . . . . . . . . . . . . . . . . 76
3.4. Énergie du champ électromagnétique . . . . . . . . . . . . . . . . . . . . . . . 77
3.5. Ondes Électromagnétiques . . . . . . . . . . . . . . . . . . . . . . . . . . . . 79
1. Équation d’Onde et Solution de d’Alembert . . . . . . . . . . . . . . . . . . 79
2. Ondes monochromatiques . . . . . . . . . . . . . . . . . . . . . . . . . . 80
3. Densité et flux d’énergie . . . . . . . . . . . . . . . . . . . . . . . . . . 81
4. Rayonnement dipolaire . . . . . . . . . . . . . . . . . . . . . . . . . . . 82
5. Formule de Larmor . . . . . . . . . . . . . . . . . . . . . . . . . . . . 85
3.A. Théorème de Helmholtz . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 85
Problèmes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 86
iv
4. Relativité Restreinte . . . . . . . . . . . . . . . . . . . . . . . . . . . 88
4.1. Transformations de Lorentz . . . . . . . . . . . . . . . . . . . . . . . . . . . 88
1. Principe de la relativité restreinte . . . . . . . . . . . . . . . . . . . . . . 88
2. Transformations de Lorentz . . . . . . . . . . . . . . . . . . . . . . . . . 89
3. Intervalle . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 90
4. Contraction des Longueurs et Dilatation du Temps . . . . . . . . . . . . . . . 91
4.2. Formalisme covariant . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 91
1. Quadri-vecteurs et quadri-tenseurs . . . . . . . . . . . . . . . . . . . . . . 91
2. Tenseur métrique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 94
3. Exemples . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 95
4.3. Formulation relativiste de l’électrodynamique . . . . . . . . . . . . . . . . . . . 97
1. quadri-vecteur du potentiel électromagnétique . . . . . . . . . . . . . . . . . 97
2. Tenseur électromagnétique . . . . . . . . . . . . . . . . . . . . . . . . . 98
3. Invariants du champ . . . . . . . . . . . . . . . . . . . . . . . . . . . . 100
4.4. Applications . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 101
1. Champ d’une charge ponctuelle en mouvement uniforme . . . . . . . . . . . . 101
2. Rayonnement d’une charge accélérée . . . . . . . . . . . . . . . . . . . . . 102
4.5. Formulation Lagrangienne . . . . . . . . . . . . . . . . . . . . . . . . . . . . 103
1. Principe de la moindre action . . . . . . . . . . . . . . . . . . . . . . . . 103
2. Action invariante d’une particule chargée dans un champ . . . . . . . . . . . . 104
3. Action du champ . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 105
4. Tenseur d’énergie-impulsion . . . . . . . . . . . . . . . . . . . . . . . . . 106
Problèmes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 109
A. Facteurs de conversion SI-gaussien . . . . . . . . . . . . . . . . . . . 110
v
Préface
Électrostatique
La charge électrique est une quantité strictement conservée lors des interactions fondamen-
tales : les particules chargées peuvent changer de nature, mais la charge totale n’est pas modifiée.
La conservation de la charge constitue l’une des symétries de base sur lesquelles se fonde notre
compréhension de l’Univers.
q1 q2
F12 = 2 r̂12 (1.1)
r12
Explications :
F12 est la force exercée par la charge q2 au point r2 sur q1 au point r1 .
r12 = |r1 − r2 | est la distance entre les deux charges.
q1 , q2 sont les charges électriques.
r̂12 = (r1 − r2 )/r12 est le vecteur unité allant de q2 vers q1
Remarques
N
X qi
Ftot. = q1 r̂ (1.2)
r2 1i
i=2 1i
2
Les unités de base du système gaussien sont le centimètre, le gramme et la seconde, d’où l’acronyme cgs.
Le dyne est l’unité de la force (1 g.cm/s2 ) et l’erg l’unité d’énergie (1 g.cm2 /s2 ).
6 1. Électrostatique
Champ électrique
Le champ électrique E(r) est la force par unité de charge exercée au point r par une distribution
quelconque de charge sur une charge test que l’on prend infiniment petite, pour ne pas perturber
la distribution :
F
E = lim (1.3)
q→0 q
Bien sûr, il s’agit ici d’une définition opérationnelle du champ électrique, à partir d’une mesure
de la force sur une charge test infinitésimale. Il ne peut s’agir d’une définition fondamentale
puisque la limite q → 0 est en pratique impossible, la charge étant quantifiée. Il est préférable de
considérer le champ électrique comme une quantité fondamentale créée par des sources et dont
l’influence sur un charge ponctuelle q se traduit par une force électrique :
F = qE (1.4)
Il est bien entendu que le champ électrique ressenti par une charge ponctuelle est causé par toutes
les charges de l’espace, à l’exception de cette charge elle-même. D’après (1.2), le champ créé par
un ensemble de charges ponctuelles est
X qi
E(r) = (r − ri ) (1.5)
i
|r − ri |3
Le champ électrique est donc un concept qui sert d’intermédiaire dans l’influence mutuelle des
charges électriques : on considère en premier lieu comment des charges peuvent créer un champ
électrique en tout point pour ensuite considérer l’influence de ce champ sur une charge donnée.
Distribution de charge
On considère souvent des distributions continues de charge qu’on note ρ(r). La charge contenue
dans un élément de volume d3 r est alors égale à dq = ρ(r)d3 r. La densité ρ correspondant à une
charge ponctuelle q située au point r0 s’écrit alors ρ(r) = δ(r − r0 ) où δ(. . .) est la fonction de
Dirac. Notons qu’en coordonnées cartésiennes, cette notation signifie
Si on désire, dans le contexte de la physique classique, que la densité ρ soit une quantité continue,
il faut évidemment se limiter à des échelles de distance suffisamment grandes pour que les fluctuations
de la charge dans un élément de volume mésoscopique soient négligeables. Par mésoscopique, on entend
une échelle de grandeur petite en comparaison des grandeurs macroscopiques, mais grande par rapport
aux échelles atomiques. La fonction delta n’est évidemment pas continue, mais son utilisation permet
d’appliquer le concept de densité de charge à des charges ponctuelles. L’utilisation d’une densité de charge
ρ est donc parfaitement générale. Dans un contexte quantique, la densité de charge associée à un électron
est e|ψ|2 , où ψ est la fonction d’onde associée : le concept de densité de charge est alors incontournable.
En vertu du principe de superposition, le champ créé par une distribution statique arbitraire
ρ est
ρ(r0 )
Z
E(r) = d3 r0 0 3
(r − r0 ) (1.7)
|r − r |
1.1. Lois de Coulomb et de Gauss 7
Les crochets [ ] indiquent que ΦE est uneH fonctionnelle de la surface S. Si la surface est fermée,
le signe d’intégrale comporte un cercle : . La loi de Gauss stipule que le flux électrique sortant
d’une surface fermée est égal à la charge totale incluse à l’intérieur de la surface :
I
da · E = 4πqtot. (1.9)
S
La loi de Gauss est donc valide pour une sphère de rayon arbitraire.
2. Pour une surface S quelconque enveloppant une charge ponctuelle, on utilise le théorème de
la divergence : I Z
da · E = d3 r ∇·E (1.11)
S V
1 ∂ 2 1 ∂ 1 ∂Eϕ
∇·E = 2
(r Er ) + (Eθ sin θ) + (1.12)
r ∂r r sin θ ∂θ r sin θ ∂ϕ
ce qui implique, pour le champ produit par une charge ponctuelle centrée à l’origine, que
∇·E = 0 sauf à la position même de la charge. Donc la forme précise de S n’a pas
d’importance car seule l’origine contribue à l’intégrale. On en conlcut que la loi de Gauss
est valable pour une surface quelconque, dans le cas d’une charge ponctuelle.
3. Dans le cas d’un ensemble de charges ponctuelles (ou d’une distribution continue) la loi de
Gauss est valide en vertu du principe de superposition.
Dans le cas d’une distribution continue de charge, la loi de Gauss s’écrit comme suit (on a
utilisé le théorème de la divergence):
Z I Z
3
4π d r ρ(r) = da · E = d3 r ∇·E (1.13)
V S V
Dans ce qui précède, nous sommes partis de la loi de Coulomb et du principe de superposition pour
arriver à la loi de Gauss (1.14). Nous aurions pu faire l’inverse, c’est-à-dire adopter comme point de départ
8 1. Électrostatique
la loi de Gauss et démontrer la loi de Coulomb. En fait, la solution générale de l’équation différentielle
(1.14) est donnée par l’Éq. (1.7). Il faut cependant noter qu’à toute solution de (1.14) on peut ajouter
le rotationnel d’un champ vectoriel quelconque A, puisque que ∇·(∇∧A) = 0. Qu’une telle addition soit
en réalité impossible ici provient d’une propriété additionnelle du champ électrique en électrostatique : le
champ de force F = qE associé au champ électrique est conservatif, c’est-à-dire que ∇∧E = 0. Ces idées
seront précisées dans l’annexe 3.A.
∇2 Φ = −4πρ (1.18)
La solution de cette équation nous est connue en principe : c’est l’Éq. (1.17). Cependant, la
solution (1.17) n’est utile en pratique que si on connait précisément la densité de charge ρ, ce qui
n’est pas souvent le cas. Il découle de notre solution de l’équation de Poisson l’identité suivante,
pour le potentiel créé par une charge unité située à r0 :
1
∇2 = 4πδ(r − r0 ) (1.19)
|r − r0 ]
Le potentiel n’est pas une fonction unique a priori, puisqu’on peut y ajouter une constante
sans changer le champ électrique. Par convention on fixe à zéro le potentiel à l’infini (si possible),
pour lever cette ambiguı̈té.
3
L’expression de l’opérateur gradient ∇ en coordonnées cartésiennes, cylindriques et sphériques est donnée
en annexe de ce manuel.
1.2. Potentiel électrique 9
Ce travail est égal au gain en énergie cinétique de la particule réalisé pendant sont trajet sur la
courbe, en vertu de l’équation du mouvement F = mv̇:
Z
W = m dl · v̇
ZC
= m dt v · v̇
ZC (1.21)
d 1 2
= m dt ( v )
C dt 2
1 1
= mv2 |final − mv2 |initial
2 2
Dans le cas de la force électrique, ce travail est l’intégrale d’un gradient, par conséquent égal à la
différence des valeurs du potentiel aux points initial et final de la trajectoire ri et rf :
Z
W =q dl · E = −q(Φ(rf ) − Φ(ri )) (1.22)
C
Le fait que la force électrique soit ‘conservative’ est évidemment lié au fait qu’elle est irrota-
tionnelle. On a en effet l’équation
qui provient directement du fait que le rotationnel d’un gradient est nul (∇∧∇ = 0). Cela signifie
que le travail de la force électrique sur une charge qui se déplace d’un point ri à un point rf ne
dépend pas du chemin emprunté pour aller de ri à rf , mais seulement de la différence de potentiel
entre ces deux points.
1
mv 2 = eΦ(x) (1.26)
2
La densité de courant est négative (= −J), puisque les électrons se dirigent vers la droite et la charge vers
la gauche. Si ρ est la densité de charge, on a J = −ρv. Cette densité de courant ne dépend pas de x, sinon
il y aurait accumulation continue de charge entre les deux plaques. Puisque la vitesse doit augmenter avec
x, on s’attend à ce que la densité de charge diminue (en valeur absolue) quand x augmente. Le potentiel
peut être déterminé par l’équation de Poisson :
r
4πJ m
Φ00 = −4πρ = = 4πJ (1.27)
v 2eΦ
d √
r
0 001 d m
ΦΦ = (Φ0 )2 = 8πJ Φ (1.28)
2 dx 2e dx
En intégrant, on obtient
m√
r
(Φ0 )2 = 16πJ Φ (1.29)
2e
La constante d’intégration est nulle : en effet, cette constante égale à E02 , où E0 est le champ électrique
à x = 0. Si E0 > 0, les électrons subissent une force qui les repousse vers la plaque, ce qui empêche tout
courant. Si, au contraire, E0 < 0, les électrons quittent trop rapidement la plaque, et le nuage électronique
crée un champ compensatoire qui ramène E0 à zéro dans un régime permanent. On a donc
Φ0 √
= 4 πJ(m/2e)1/4 (1.30)
Φ1/4
ou encore
h i2/3
x4/3
p
Φ(x) = 9πJ m/2e (1.32)
Par conséquent, le champ électrique et la densité de charge ont les comportements suivants :
Prenons maintenant la limite continue de cette expression. On suppose que tout élément de
volume d3 r, aussi petit soit il, contient une charge ρ(r)d3 r et que la charge est divisible à volonté.
L’expression (1.36) pour l’énergie potentielle de la configuration devient
Z Z
1 3 1
U= d r ρ(r)Φ(r) = − d3 r Φ∇2 Φ (1.37)
2 8π
où nous avons utilisé l’équation de Poisson. La différence subtile entre cette expression et
l’Éq. (1.36) est qu’un élément de charge ρ(r)d3 r, dans la limite ou l’élément de volume tend
vers zéro, est trop petit pour contribuer appréciablement au potentiel Φ(r) à ce point précis et
donc on néglige la différence entre Φ(i) (ri ) et Φ(ri ). Par conséquent, si on applique l’expression
(1.37) à une distribution de charges ponctuelles, on trouve plus que l’énergie nécessaire à mettre
en place les charges ponctuelles : on trouve aussi l’énergie nécessaire à construire chacune des
charges ponctuelles à partir d’un ‘fluide’ électrique infiniment divisible.
Intégrons par parties l’expression (1.37), en utilisant la décomposition
Φ∇2 Φ = ∇·(Φ∇Φ) − ∇Φ · ∇Φ (1.38)
Comme E = −∇Φ, on obtient
Z I
1 3 1
2
U= drE + da · E Φ (1.39)
8π 8π ∞
où la seconde intégrale a été convertie en intégrale de surface par le théorème de la divergence.
L’indice ∞ signifie que cette intégrale est calculée sur une surface (une sphère, par exemple) située
à l’infini, c’est-à-dire entourant tout l’espace. Cette intégrale s’annule pour une distribution de
charge concentrée dans une région finie de l’espace. En effet, du point de vue d’une sphère quasi-
infinie de rayon R au centre de laquelle se situe une distribution de charge d’étendue minime en
comparaison de R, le champ électrique et le potentiel se comportent respectivement comme Q/R2
et Q/R, Q étant la charge totale de la distribution. Comme l’aire de la sphère est proportionnelle
à R2 , l’intégrale de surface se comporte au total comme ∼ 1/R et tend vers zéro quand R → ∞.
Il nous reste donc que le premier terme de (1.39). Il est alors légitime d’identifier E2 /8π à une
densité d’énergie électrique E:
1 2
E = E (1.40)
8π
12 1. Électrostatique
L’intégrale de cette densité est l’énergie électrostatique totale de l’électron, qu’on suppose égale à son
énergie au repos :
Z ∞ 2
2 1 e 1
mc = 4πr2 dr = e2 /rc (1.42)
8π rc r4 2
e2
r0 = = 2, 82 × 10−13 cm (1.43)
mc2
Ce rayon n’est malheureusement pas le rayon réel de l’électron : les expériences de la physique des particules
n’ont pas réussi à détecter une structure interne pour l’électron jusqu’à des distances de l’ordre de 10−16 cm.
Le rayon classique de l’électron est plutôt une grandeur caractéristique associée à l’interaction des électrons
avec la lumière, comme nous le verrons plus tard.
Si la surface S est fermée, le courant qui en sort est relié à la dérivée de la charge totale
incluse dans S, en vertu de la conservation de la charge électrique :
I Z
∂
da · J = − d3 r ρ(r) (1.45)
S ∂t V
1.3. Électrostatique et conducteurs 13
J = ρv (1.48)
la mesure d’intégration d3 r a les unités L3 , la fonction delta doit compenser pour que le résultat soit sans
unités.
14 1. Électrostatique
n /
E=0
vide
E ≈ 4πρs n (1.53)
E=0
conducteur
Cette expression se démontre en appliquant la loi de Gauss sur une surface fermée infiniment
mince qui chevauche la surface du conducteur, comme illustré.
La densité volumique d’énergie associée à la surface d’un conducteur est
1 2
E = 2πρ2s (1.54)
8π
Cette quantité est aussi égale à la pression exercée par le conducteur sur sa surface : l’énergie
gagnée par une expansion de la surface du conducteur (déplacement virtuel) sur une distance δx
5
En pratique, la conductivité dépend de la fréquence du champ électrique appliqué. On donne ici les valeurs
de σ à fréquence nulle. La loi d’Ohm (1.50) doit être comprise comme reliant les transformées de Fourier en
fréquence : J(ω) = σ(ω)E(ω).
1.3. Électrostatique et conducteurs 15
pour un élément de surface δA est égale à 2πρ2s δxδA, ce qui mène à la pression citée. C’est aussi
la force exercée sur l’élément de surface par tous les autres éléments de surface du conducteur :
le champ E peut être décomposé en deux parties E1 et E2 , dues respectivement à l’élément
de surface à r et à la somme de toutes les autres sources. Tout juste à l’intérieur (resp. à
l’extérieur) du conducteur, on a E1 = −2πnρs (resp. +2πnρs ), en vertu du théorème de Gauss.
Donc E2 = 2πnρs tout juste à l’intérieur, puisque E1 + E2 = 0 à cet endroit. Or la force par unité
de surface exercée sur l’élément est ρs E2 .
1.3.4 Capacité
En vertu du principe de superposition, le potentiel d’un conducteur est proportionnel aux
charges qui sont placées dans le système. Soit Qi les charges placées dans un système de N
conducteurs (i = 1, . . . , N ) et Φi le potentiel qui en résulte sur le conducteur i. On a la relation
linéaire
X
Φi = Bij Qj (1.55)
j
où les coefficients Bij forment une matrice carrée qui ne dépend que de la configuration spatiale
des conducteurs en question. On préfère généralement considérer la relation inverse :
X
Qi = Cij Φj (1.56)
j
où les coefficients Cii sont appelés coefficients de capacité car ils sont une mesure de la capacité
d’un conducteur à absorber de la charge électrique à un potentiel donné. Les coefficients Cij
(i 6= j) sont appelés coefficients d’influence électrostatique. Pour un conducteur simple, on a
Q = CΦ où C est la capacité du conducteur.
L’énergie potentielle d’un système de conducteurs est donc
1X 1X
U= Qi Φi = C ΦΦ (1.57)
2 i 2 ij ij i j
pour un conducteur simple, cela revient à U = 12 CΦ2 . On peut donc introduire une définition
différente des coefficients Cij , comme étant les dérivées partielles secondes de l’énergie électrosta-
tique par rapport aux potentiels :
∂2U
Cij = (1.58)
∂Φi ∂Φj
D’après cette définition, il est évident que la matrice C est symétrique (Cij = Cji ) car l’ordre de
différentiation est sans importance.
Remarques :
1. L’énergie U doit être une forme quadratique définie positive pour que l’état de moindre
énergie soit Qi = 0.
2. On peut démontrer en général que Cii > 0 et Cij < 0 si i 6= j (le signe de la charge induite est
l’opposée du signe du potentiel inducteur).
16 1. Électrostatique
On a donc
ab ab b2
C11 = , C12 = C21 = − , C22 = (1.62)
b−a b−a b−a
1.4 Multipôles
Dans cette section nous appliquons la formule (1.17) au cas où la distribution ρ(r0 ) est
confinée dans un espace petit comparé à la distance au point d’observation. Il est intuitivement
clair que le potentiel créé par une distribution de charge doit ressembler au potentiel d’une charge
ponctuelle q (q est la charge totale de la distribution) si on se place suffisamment loin des sources.
Cependant, des corrections existent à cette approximation, qui sont capitales si la charge totale
est nulle. Si a est la taille typique de la distribution de charge (plus précisément le rayon de la plus
petite sphère qui englobe toute la distribution) et si r est la distance entre le point d’observation
et la distribution, un développement systématique du potentiel en une série de puissances de a/r
est possible : c’est le développement multipolaire.
Une connaissance des polynômes de Legendre est essentielle pour la compréhension de cette
section; voir l’annexe à la fin de ce chapitre. Le résultat principal dont nous avons besoin est le
développement en série suivant [cf. Éq. (1.209)]:
X r0l ∞
1
= P (t)
|r − r0 | rl+1 l
l=0
1.4. Multipôles 17
où r ≡ |r|, r0 ≡ |r0 |, t = cos 6 rr0 et Pl est le polynôme de Legendre d’ordre l. Pour que cette
série converge, il faut que r0 < r. Si nous substituons ce développement dans la formule (1.17),
l’expression du potentiel devient à son tour une série infinie :
∞
l
X
Φ(r) = Φ(2 ) (r) (1.63)
l=0
l
où Φ(2 ) (r) est le potentiel créé par le multipôle d’ordre l:
Z
l 1
Φ(2 ) (r) = l+1 d3 r0 ρ(r0 )r0l Pl (cos γ) (1.64)
r
l
Le potentiel multipolaire Φ(2 ) (r) diminue avec la distance comme 1/rl+1 . De plus, il dépend de
l
la direction, sauf pour l = 0. Si le rayon caractéristique de la distribution est a, Φ(2 ) sera aussi
proportionnel à al fois une quantité ayant les unités de la charge électrique.
Le premier terme de la série (1.64) est le potentiel monopolaire (l = 0):
Z
1 q
Φ(1) (r) = d3 r0 ρ(r0 ) = (1.65)
r r
Il s’agit du potentiel obtenu en négligeant la structure interne de la distribution de charge, c.-à-d.,
en ne considérant que la charge totale, ou monopôle. Ce terme domine les autres aux grandes
distances, sauf s’il s’annule.
1
Φ(2) (r) = d·r (1.68)
r3
Dans le cas où la charge totale de la distribution est nulle, ce potentiel est le terme est le plus
important de la série (1.64).
Le champ électrique correspondant est
1
E = −∇(1/r3 )d · r − (1/r3 )d = 3(d · r)r − dr2
5
(1.69)
r
Remarquons que le champ dipolaire varie comme 1/r3 et qu’il est deux fois plus fort le long de
l’axe du dipôle que perpendiculairement et a une direction opposée.
18 1. Électrostatique
Comme exemple d’une distribution ponctuelle de charge nulle possédant un moment dipolaire
non nul, considérons une charge −q placée à l’origine (r0 = 0) et une charge +q située à r0 = a.
La distribution de charge est alors
La définition (1.67) donne d = qa. En d’autres termes, le moment dipolaire est selon l’axe qui
relie les deux charges, vers la charge positive. Sa valeur absolue est alors qa, où a est la distance
entre les charges.
La charge totale (monopôle) ne dépend pas de l’origine des coordonnées. Cependant, le
moment dipolaire dépend en général du choix de l’origine. Si on calcule le moment dipolaire par
rapport à une nouvelle origine O0 située à une position a par rapport à l’ancienne origine O, on
trouve Z
d0 = d3 r ρ(r)(r − a) = d − qa (1.71)
où d est le moment dipolaire par rapport à l’origine O. On remarque que si q = 0, alors le moment
dipolaire est indépendant de l’origine. Il s’agit d’une propriété générale : le moment multipolaire
non nul le plus bas ne dépend pas de l’origine des coordonnées, mais tous les autres en dépendent.
Ici nous avons utilisé la convention de sommation sur les indices répétés et nous avons défini le
tenseur quadrupolaire:
Z
Qij = d3 r0 ρ(r0 )(3x0i x0j − δij r02 ) (1.73)
Remarques
1. La trace du tenseur quadrupolaire s’annule : Qii = 0 (indice sommé). Donc on peut écrire
1
Φ(4) (r) = Q xx (1.74)
2r5 ij i j
2. Ce tenseur est symétrique (Qij = Qji ). Par conséquent, il est possible de trouver un système
d’axes orthogonaux (dits axes principaux du tenseur) dans lequel le tenseur est diagonal.
Celui-ci comporte alors trois composantes (Q11 , Q22 et Q33 ) reliées par la contrainte de trace
nulle : Q11 + Q22 + Q33 = 0.
1.4. Multipôles 19
3. Si la distribution de charge possède une symétrie de rotation autour d’un axe (ex. z) alors
Qxx = Qyy . En appelant Qzz = Q le moment quadrupolaire de la distribution, il s’ensuit que
Qxx = Qyy = − 21 Q. Dans ce cas, le potentiel quadrupolaire devient
Q
Φ(4) (r) = (2z 2 − x2 − y 2 ) (1.75)
4r5
ou, en coordonnées polaires,
Q
Φ(4) (r) = (3 cos2 θ − 1) (1.76)
4r3
4. Si ρ > 0 et que la distribution est allongée, alors Q > 0. Si, au contraire, la distribution est
aplatie, alors Q < 0. L’inverse est vrai si ρ < 0. Le moment quadrupolaire d’une distribution
de charge comportant une densité à peu près uniforme nous permet de caractériser la forme
de cette distribution.
Figure 1.1. Équipotentielles d’un dipôle (à gauche) et d’un quadrupôle à symétrie
azimutale (à droite)
5. On peut obtenir un quadrupôle en plaçant côte à côte deux dipôles de signes opposés (ce
procédé est général pour obtenir un l-pôle à l’aide de deux (l − 1)-pôles).
où γ est l’angle entre les point (θ0 , ϕ0 ) et (θ, ϕ) sur la sphère (on utilise les coordonnées sphériques
pour r et r0 ). On obtient alors le développement
0 l
1 1 X 4π r ∗
= Ylm (θ0 , ϕ0 )Ylm (θ, ϕ) (1.78)
|r − r0 | r 2l + 1 r
l,m
20 1. Électrostatique
L’avantage de cette formulation en fonction des harmoniques sphériques est que l’on connait bien
la dépendance angulaire de ces derniers, ainsi que la façon dont ils se transforment lorsqu’on
effectue une rotation ou une inversion de l’espace.
Cette expression est une extension simple de l’énergie potentielle qΦ(r) d’une charge ponctuelle
dans un potentiel. Nous devons comparer cette expression avec les formules (1.37) et (1.39) pour
l’énergie interne de la distribution, qui est l’énergie nécessaire pour assembler les charges qui
la forment. Dans (1.37) le potentiel Φ est celui-là même qui est produit par la distribution ρ,
alors qu’ici il s’agit d’un potentiel externe produit par une autre source et qui satisfait donc
∇2 Φ = 0 dans la région de l’espace considérée ici. On s’intéresse ici à l’énergie d’interaction de la
distribution avec le potentiel électrostatique externe.
Pour évaluer approximativement cette énergie, on peut développer Φ(r) en série de Taylor
autour de l’origine des coordonnées, qu’on place convenablement quelque part au centre de notre
distribution :
1
Φ(r) = Φ(0) + xi ∂i Φ(0) + xi xj ∂i ∂j Φ(0) + . . . (1.82)
2
En substituant ce développement dans l’expression pour W , on trouve
Z Z
1
W = qΦ(0) + ∂i Φ(0) d r ρ(r)xi + ∂i ∂j Φ(0) d3 r ρ(r)xi xj + . . .
3
(1.83)
2
où q est la charge totale de la distribution. Le deuxième terme est simplement −E · d et le
troisième peut être écrit comme
Z
1 3 2
∂ ∂ Φ(0) Qij + 3δij d r ρ(r)r (1.84)
6 i j
Remarques
On constate que la force ressentie par un dipôle dépend non pas du champ lui-même, mais de son
gradient :
Fdip. = pi ∂i E0
ou Fdip. = (d · ∇)E0 (1.89)
De même, le couple qui s’exerce sur un dipôle électrique est, à l’ordre zéro,
Z
N= d3 r ρ(r)r ∧ E(0) = d ∧ E(0) (1.90)
Ces résultats sont tout-a-fait compatibles avec l’énergie d’interaction −E · d du dipôle dans un
champ électrique.
On peut facilement observer dans la vie courante les effets de la force (d · ∇)E agissant sur un dipôle.
Lorsqu’un objet devient chargé d’électricité – par exemple, un peigne qu’on vient de frotter à ses cheveux –
il peut attirer des objets légers, comme des bouts de papier. On peut expliquer le phénomène comme suit :
la charge déposée sur le peigne crée un champ électrique diminuant avec la distance qui polarise le papier.
Les dipôles microscopiques contenus dans le papier sont parallèles au champ et attirés dans la direction où
le champ augmente, c’est-à-dire vers le peigne. Le même peigne peut aussi attirer vers lui un mince filet
d’eau coulant d’un robinet, car l’eau est une substance très polarisable (voir la section sur les diélectriques).
La notion de moment quadrupolaire peut être appliquée entre autres à la description des noyaux
atomiques. L’état quantique d’un noyau peut être caractérisé notamment par son spin I et la projection M
de ce spin selon l’axe z. Dans un état quantique donné |I, M i, la distribution de charge peut être telle que
le moment quadrupolaire Q du noyau est non nul (en raison de la symétrie azimutale de l’état quantique
considéré, le tenseur Qij n’a qu’une seule composante indépendante non nulle). Si le noyau est placé à un
endroit où le champ électrique possède un gradient non nul, par exemple dans une position donnée d’une
structure cristalline, alors le terme quadrupolaire de l’énergie d’interaction (1.85) s’ajoute à l’hamiltonien
H du noyau et lève la dégénérescence entre les états |I, M i ayant des valeurs de M différentes. On peut
alors observer des sauts quantiques entre ces états, à l’aide d’ondes radio : c’est ainsi que la forme des
noyaux est observée.
22 1. Électrostatique
Cette constante est évidemment zéro dans le cas des conditions aux limites de Dirichlet
3. La solution Φ n’admet aucun maximum ou minimum dans le domaine où ∇2 Φ = 0. Dans le
cas contraire, il y aurait un flux électrique non nul sur une sphère infinitésimale entourant
l’extrémum, ce qui impliquerait une charge non-nulle à l’intérieur, ce qui n’est pas le cas
puisque ∇2 Φ = 0. On en tire la conclusion qu’il ne peut y avoir de configuration statique de
charges ponctuelles qui soit stable, car aucune de ces charges ne peut trouver un minimum
local de son énergie potentielle.
Le plan z = 0 est alors une équipotentielle à Φ = 0 et a le comportement correct autour de (0, 0, a). Il
s’agit donc de la solution unique. Au loin, le système ressemble à un dipôle de moment d = 2aqẑ. L’énergie
1.5. Méthodes de solution de l’équation de Laplace 23
potentielle de la charge q en présence du plan est U = −q 2 /(2a). La force exercée sur la charge est donc
donnée par la dérivée par rapport à a de l’énergie potentielle :
∂U q2
F=− = − 2 ẑ (1.94)
∂a 2a
Cette force n’est pas égale à la force de Coulomb entre l’image et la charge (q 2 /4a2 ): la raison en est
qu’il est impossible de déplacer la charge sans que l’image soit également déplacée; la notion de source
indépendante de la charge-test ne s’applique alors pas et F 6= qE. Dans ce genre de problème, il faut
calculer la force à partir du gradient de l’énergie potentielle.
Nous cherchons des solutions particulières (les fonctions de base mentionnées ci-haut) sous la
forme d’un produit de fonctions à une variable :
X 00 Y 00 Y 00
+ + =0 (1.97)
X Y Y
Cette équation étant valide pour toutes valeurs de x, y et z, chacun des termes doit être constant :
X 00 Y0 Z 00
= a2 = b2 = −a2 − b2 (1.98)
X Y Z
24 1. Électrostatique
Nous avons écrit des constantes au carré par convenance, sans supposer nécessairement que a2 et
b2 sont des nombres réels positifs : en réalité, a et b peuvent être complexes. Ce sont les conditions
aux limites qui imposeront des restrictions sur a et b.
La solution correspondant à des valeurs particulières de a et b est alors
√
a2 +b2
Φ(r) = C(a, b)e±ax e±by e±iz (1.99)
où C(a, b) est une constante quelconque. Dans chaque facteur, le signe + ou − peut être choisi :
il y a donc 8 solutions indépendantes pour chaque valeur de a et b. Les fonctions cosh(ax) et
sinh(ax) peuvent être utilisées au lieu de e±ax . De même, cos(αx) et sin(αx) peuvent être utilisés
si a = iα est imaginaire.
et
2
Zm,n (z) = sinh(γmn (c − z)) γmn = π 2 (m2 /a2 + n2 /b2 ) (1.102)
Ces fonctions ne sont pas normalisées. Toutes ces fonctions satisfont aux conditions aux limites sur toutes
les faces sauf sur z = 0. De plus, elles sont orthogonales. On cherche donc une combinaison linéaire telle
que
X
cmn Φmn (x, y, 0) = V ∀x, y
m,n
1
hΦmn |Φrs i = ab sinh2 (γmn c)δmr δns (1.105)
4
1.5. Méthodes de solution de l’équation de Laplace 25
Si on prend le produit scalaire de Φrs avec la condition aux limites (1.103), on obtient
16V
crs = (1.107)
rsπ 2 sinh(γmn c)
Nous avons donc trouvé les coefficients de la combinaison linéaire unique qui satisfait à toutes les conditions
aux limites. Le problème est résolu.
∂2Φ
1 ∂ 2 ∂Φ 1 ∂ ∂Φ 1
r + sin θ + =0 (1.108)
r2 ∂r ∂r r2 sin θ ∂θ ∂θ r2 sin2 θ ∂ϕ2
Profitons de notre expérience en mécanique quantique pour accélérer notre analyse : l’opérateur
du carré du moment cinétique L2 s’écrit, à un facteur multiplicatif h̄2 près, comme suit
1 ∂2
1 ∂ ∂
L2 = − sin θ + (1.110)
sin θ ∂θ ∂θ sin2 θ ∂ϕ2
Les deux membres de cette équation doivent être constants séparément. Dans le cas du deuxième
terme, on connaı̂t la solution : il s’agit des harmoniques sphériques Ylm (θ, ϕ):
s
2l + 1 (l − m)! m
Ylm (θ, ϕ) = P (cos θ)eimϕ (1.113)
4π (l + m)! l
X Blm
Φ(r) = Am
l r
l
+ l+1 Ylm (θ, ϕ) (1.117)
r
l,m
Dans le cas d’un problème à symétrie azimutale, seuls les termes avec m = 0 contribuent :
X Bl
l
Φ(r) = Al r + l+1 Pl (cos θ) (1.118)
r
l
Notez que les constantes Al et Bl ne sont pas les mêmes ici que dans l’équation précédente.
Si on désire faire abstraction d’une connaissance quelconque de la mécanique quantique, alors
il faut démontrer de manière indépendante la solution générale (1.117). Indiquons plutôt comment
démontrer la solution générale (1.118) dans le cas d’une symétrie azimutale, ce qui est plus simple.
On suppose alors que la solution est de la forme Φ(r, θ) = R(r)Θ(θ). La séparation des variables
nous donne l’équation (1.115) pour la fonction radiale, où la constante l est à priori quelconque.
Les solutions R(r) sont les mêmes que ci-haut. La fonction Θ, elle, satisfait à l’équation de
Legendre (1.207) lorsqu’on l’exprime en fonction de t ≡ cos θ. Des solutions à cette équation
existent pour toutes les valeurs réelles de l, mais c’est uniquement quand l est un entier non
négatif que ces solutions sont finies à cos θ = ±1, ce qui est requis ici. Dans ce cas, les solutions
sont les polynômes de Legendre et on trouve en fin de compte la solution générale (1.118).
Exemple 1.6: Potentiel autour d’une sphère métallique dans un champ uniforme
Une sphère conductrice est placée dans un champ électrique uniforme. Trouver comment l’introduction
de cet objet perturbe le potentiel électrique.
Solution : Dans le but de minimiser son énergie, la sphère conductrice va accumuler des charges
négatives dans la direction −z et des charges positives dans la direction +z, de façon à créer une
distribution dipolaire induite. Voyons comment ceci fonctionne quantitativement. On sait que le potentiel
à l’infini a la forme Φ = −E0 z = −E0 r cos θ ou encore Φ = −E0 rP1 (cos θ). En raison de la symétrie
azimutale, on écrit la solution sous la forme (1.118). Quand r → ∞, seul P1 doit contribuer et donc Al = 0
pour l 6= 1 et A1 = −E0 . D’autre part, Φ doit être constant sur la sphère (r = a):
X Bl
A1 a cos θ + Pl (cos θ) = const. (1.119)
al+1
l
Pour que cette identité soit valide ∀θ, il faut que les coefficients des puissances de cos θ soient identiques de
part et d’autre :
Bl = 0 si l > 1 − E0 a + B1 /a2 = 0 B0 /a = const. (1.120)
1.5. Méthodes de solution de l’équation de Laplace 27
On voit que B0 = q, la charge nette déposée sur la sphère, c.-à-d. une donnée du problème. La dépendance
supplémentaire en 1/r2 est celle d’un dipôle électrique induit par le champ externe. Le moment dipolaire
correspondant a comme valeur absolue p = E0 a3 . Nous aurions pu conclure dès le départ que la charge
induite sur la sphère ne comporterait aucun moment multipolaire supérieur, pour la raison suivante :
comme le champ extérieur dans lequel la sphère est plongée ne possède aucun gradient ou dérivée d’ordre
supérieur, l’énergie électrostatique de la sphère ne peut être abaissée en créant un quadrupôle induit (ou
aucun moment d’ordre supérieur). De tels moments multipolaires sont alors inutiles et ne sont pas induits.
(int)
X (ext)
X Bl
Φs (r) = Al rl Pl (cos θ) Φs (r) = Pl (cos θ) (1.122)
rl+1
l l
(int)
Nous avons utilisé la symétrie azimutale du problème et la condition que Φs doit être fini quand r = 0
(ext)
et Φs fini quand r → ∞. D’autre part, d’après le développement (1.62) on peut écrire
X l
r
q Pl (cos θ) (r < a)
al+1
l
Φq (r) = (1.123)
X al
q Pl (cos θ) (r > a)
rl+1
l
(int) (ext)
Sur la sphère (r = b) le potentiel total peut être écrit à la fois à l’aide de Φs et de Φs :
X al
l
Φ(b, θ) = Al b + q Pl (cos θ)
bl+1
l
(1.124)
X Bl + qal
= Pl (cos θ)
bl+1
l
Comme cette expression doit être une constante indépendante de θ, chaque coefficient de Pl doit être nul,
sauf peut-être celui de P0 . On a donc Al = −qal /b2l+1 et Bl = −qal , pour l > 0. Le potentiel de la sphère
28 1. Électrostatique
est alors A0 + q/b = B0 + q/b. Comme la sphère est neutre, le potentiel à l’infini doit être q/r: on en conclut
que A0 = B0 = 0. La solution complète est donc
q X qrl 2l+1
Φ(r, θ) = + 1 − (a/b) Pl (cos θ) (r < a)
a al+1
l=1
q X qal 2l+1
(1.125)
= + 1 − (r/b) Pl (cos θ) (a < r < b)
r rl+1
l=1
q
= (r > b)
r
Remarquez que le potentiel à l’extérieur de la sphère est isotrope, en dépit du fait que la distribution
de charge sur la sphère et à l’intérieur de celle-ci n’est pas isotrope. C’est une propriété générale : le
potentiel à l’extérieur d’une sphère conductrice sera toujours q/r si aucune charge n’est située à l’extérieur,
qu’importe la distribution spatiale de la charge q à l’intérieur de la sphère. C’est une conséquence du
théorème d’unicité où les conditions aux limites sont ici fixées sur la sphère et à l’infini.
1 ∂2Φ ∂2Φ
1 ∂ ∂Φ
r + 2 + =0 (1.126)
r ∂r ∂r r ∂ϕ2 ∂z 2
1 ∂ 2 Q r2 ∂ 2 Z
r ∂ ∂R
r + + =0 (1.128)
R ∂r ∂r Q ∂ϕ2 Z ∂z 2
∂2Q
= −n2 Q ⇒ Q(ϕ) = e±inϕ (1.129)
∂ϕ2
n2 1 ∂2Z
1 ∂ ∂R
r − 2 + =0 (1.130)
rR ∂r ∂r r Z ∂z 2
Notez que l’ordre de la fonction de Bessel (ou de Neumann) est l’exposant angulaire n. La solution
générale est donc
X
Φ(r, ϕ, z) = (Amn Jn (km r) + Bmn Nn (km r))
m,n (1.135)
(cos(nϕ) + Cmn sin(nϕ))(cosh(km z) + Dmn sinh(km z))
On a supposé que k peut prendre un ensemble discret de valeurs {km }. Si k prend des valeurs
imaginaires, alors on doit remplacer les fonctions de Bessel et de Neumann par les fonctions de
Bessel modifiées In et Kn .
Le cas où k = 0 doit être considéré séparément. Le potentiel ne dépend alors pas de s et le
problème est réduit à celui de l’équation de Laplace sur le plan xy. L’équation radiale est dans ce
cas
d dR
r r − n2 R = 0 (1.136)
dr dr
dont la solution est R(r) = r±n , sauf si n = 0, dans lequel cas elle est R = const. ou R = ln r. On
peut donc écrire la solution générale :
∞
X B
Φ(r, ϕ, z) = A0 + B0 ln r + An r + nn
n
(cos(nϕ) + Cn sin(nϕ)) (1.137)
n=1
r
1.6 Diélectriques
Dans le but de comprendre l’effet d’un champ électrique sur un matériau, il est souhaitable
de comprendre en premier l’effet d’un champ électrique sur une molécule ou un atome de ce
matériau. Considérons à cet effet un atome isolé, en présence d’un champ électrique E. Sous
l’influence de ce champ, le nuage électronique de l’atome est déplacé par rapport au noyau et
l’atome neutre développe un moment dipolaire électrique non nul d. Si le champ appliqué est
faible en comparaison du champ électrique typique que ressent un électron dans l’atome, on est
en droit de supposer que le moment dipolaire induit d sera proportionnel au champ appliqué:
d = αE (1.138)
eE
eE = mω02 x =⇒ x= (1.140)
mω02
Eh̄2 a20
d = ex ∼ (1.141)
me2
Signalons que Ea ≡ e/a20 est le champ électrique typique qui doit exister dans un atome. On peut
donc finalement écrire
E h̄2 E
d∼ = ea (1.142)
Ea me Ea 0
où ea0 serait le moment dipolaire typique d’un électron dans un atome si l’électron était à un
endroit précis dans l’atome, plutôt que réparti dans tout le volume de l’atome. On voit que le
moment dipolaire induit est faible devant ea0 si le champ électrique appliqué E est petit devant
le champ atomique Ea . On constate aussi que le moment dipolaire induit est proportionnel au
champ appliqué , avec une polarisabilité d’ordre α ∼ ea0 /Ea = a30 (la polarisabilité a les unités
d’un volume).
La polarisabilité d’un atome ou d’une molécule attribuable au réaménagement du nuage
électronique – tel que décrit très approximativement ci-dessus – est qualifiée de polarisabilité
électronique et donne lieu à un moment dipolaire somme toute assez faible. La réponse d’un
matériau à un champ externe est beaucoup plus importante si celui-ci est composé de molécules
1.6. Diélectriques 31
polaires, c’est-à-dire possédant un dipole électrique permanent (la molécule d’eau est l’exemple
le plus connu de molécule polaire). L’existence d’un moment dipolaire permanent d0 dans un
molécule est liée à la forme précise des orbitales moléculaires qui effectuent la liaison entre les
différents atomes de la molécule : ces orbitales ne sont pas centrées géométriquement sur l’atome
central de la molécule (l’atome d’oxygène, dans le cas du H2 O) et il peut en résulter un moment
dipolaire électrique permanent. En l’absence de champ électrique externe, les molécules sont
cependant orientées au hasard, de sorte que le moment dipolaire moyen dans un petit volume est
nul. Ce type de polarisabilité d’un milieu est qualifié de paraélectrique.
Étudions l’effet d’un champ électrique uniforme E appliqué à une molécule polaire placée dans
un environnement à température absolue T . Si θ est l’angle entre le moment dipolaire d0 de la
molécule et le champ électrique et si ϕ est l’angle azimutal associé, l’énergie de la molécule dans
ce champ est W = −d · E = −d0 E cos θ. La probabilité n(θ, ϕ) que la molécule soit orientée dans
la direction (θ, ϕ) est donnée par la distribution de Boltzmann :
1 1 d0 E cos θ
n(θ, ϕ) = ed0 E cos θ/kB T ≈ 1+ (1.143)
Z Z kB T
où nous avons effectué un développement de Taylor au premier ordre, en supposant que le champ
appliqué est faible, c’est-à-dire que l’énergie du dipôle permanent dans le champ est petite en
comparaison de l’énergie thermique typique : d0 E kB T (on dit qu’il s’agit d’un ‘développement
à température élevée’). La constante de normalisation Z est déterminée en sommant sur toutes
les orientations possibles :
Z
1
1= sin θdθdϕ n(θ, ϕ) = 4π =⇒ Z = 4π (1.144)
Z
(seul le premier terme de (1.143) contribue). On peut ensuite calculer la composante moyenne du
moment dipolaire dans la direction du champ électrique :
Z
d0
hd0 cos θi = sin θdθdϕ n(θ, ϕ) cos θ
4π
Z
d0 d0 E
= sin θdθdϕ cos2 θ (1.145)
4π kB T
d2 E
= 0
3kB T
(ici, seul le deuxième terme de (1.143) contribue). La relation entre le moment dipolaire moyen
et le champ électrique est encore linéaire et permet d’introduire encore un fois la notion de
polarisabilité:
d20
hdi = αE α= (1.146)
3kB T
Cette fois-ci, la polarisabilité résulte d’une réorientation moyenne des molécules en présence d’un
champ, et varie en raison inverse de la température, pour des températures suffisamment élevées.6
Cette dépendance en température est appelée loi de Curie et est aussi observée dans le domaine
magnétique.
6
Cette dépendance en 1/T n’est vraie qu’à volume constant et non à pression constante. Elle est assez bien
vérifiée expérimentalement.
32 1. Électrostatique
P · (r − r0 )
Z
Φ(r) = d3 r0
|r − r0 |3
Z
1
= d3 r0 P · ∇0 (1.147)
|r − r0 |
0
3 0 ∇ ·P
Z I
P
=− d r 0
+ da ·
|r − r | |r − r0 |
Dans la dernière équation, nous avons intégré par parties et l’intégrale de surface est prise sur
la surface de l’échantillon de matériau diélectrique. L’interprétation est la suivante : le potentiel
créé par la distribution de dipôles équivaut à celui créé par une densité de charge volumique
ρ0 = −∇·P plus une densité superficielle de charge ρ0s = P · n, où n est la normale qui sort de
l’échantillon :
ρ0 = −∇·P ρ0s = P · n (1.148)
Cette distribution de charge – volumique ou superficielle – est bien sûr liée au matiériau et mérite
de ce fait le surnom de charge liée. Pour la distinguer de la charge externe au matériau qu’on
pourrait y introduire artificiellement, on affuble d’un prime (0 ) les symboles relatifs à la charge
liée.
La définition de P donnée plus haut (et dans les manuels les plus connus) est ambigüe. En effet, si on
imagine un cristal de type ionique dans lequel les charges positives sont déplacées en moyenne par rapport
aux charges négatives, comment doit-on définir le dipôle microscopique créé à chaque position du cristal?
Pour une charge positive donnée, entourée de charges négatives, le dipôle microscopique dépend du voisin
négatif qu’on choisit pour le former. On peut ainsi obtenir diverses valeurs de P, physiquement équivalentes.
La solution de ce paradoxe est que deux valeurs différentes de P obtenues de cette façon diffèrent par un
1.6. Diélectriques 33
vecteur constant P0 . La charge liée au sein du diélectrique est indépendante de cette constante, alors que
la charge liée à la surface du diélectrique provient des ions non appariés laissés pour compte par les deux
définitions de P. Ces ions non appariés ne sont pas les mêmes dans les deux définitions et la différence de
densité superficielle ρs qui en résulte est précisément −P0 · n. Ainsi, la valeur précise de P au sein d’un
diélectrique n’a pas de sens physique absolu : c’est la variation de P (dans l’espace et dans le temps) qui
est importante. Bien sûr, si le milieu est linéaire, cette ambiguité dans la définition de P est levée par la
relation P ∝ E.
D = E + 4πP (1.149)
L’utilité du champ D provient de ce que ses sources sont les charges libres, en principe
contrôlées par l’observateur, et non les charges liées au diélectrique.
Remarques
3. Il existe aussi des milieux non linéaires, dans lesquels un terme quadratique relativement
important s’ajoute au terme linéaire. Les conséquences sont importantes dans certaines tech-
niques optiques (fibres optiques).
34 1. Électrostatique
4. Dans un milieu linéaire et isotrope, la relation entre les charges libre et liée est aussi linéaire :
χe 1−ε
ρ0 = −χe ∇·E = − ∇·D = ρ (1.153)
ε ε
En appelant ρtot. la densité de charge totale (libre + liée) on obtient
1
ρtot. = ρ + ρ0 = ρ (1.154)
ε
où ρs est la densité superficielle de charge libre. En d’autres mots, la composante de D per-
pendiculaire à l’interface est discontinue par la densité superficielle de charge libre, alors que la
composante de E parallèle à l’interface est continue. Si aucune charge libre n’est présente dans le
système, alors la composante normale de D est continue. Dans des milieux linéaires et isotropes,
cela revient à dire que la composante normale de E est discontinue et sa composante parallèle
continue.
Justifions maintenant la première de ces conditions aux limites. Pour ce faire, considérons
encore une fois une surface abstraite de type ‘boı̂te à pilule’ à cheval sur l’interface, et intégrons
la relation ∇·D = 4πρ dans cette surface. Le théorème de la divergence donne
Z
∇·D = (D2 − D1 ) · nA = 4πρs A (1.156)
V
(ici A est l’aire de la surface d’intégration) ce qui est bien la condition ci-haut. Pour justifier la
deuxième condition, on considère plutôt une boucle d’intégration rectangulaire, infiniment étroite
et qui longe l’interface de part et d’autre. Soit m̂ la normale à cette boucle. Le vecteur qui décrit
le côté long de la boucle est de longueur l et parallèle à m̂ ∧ n. Intégrons à travers cette boucle le
flux de ∇∧E = 0. On trouve, avec le théorème de Stokes :
Comme cette relation est vraie pour n’importe laquelle orientation du contour (pourvu qu’il soit
dans le plan) ou encore pour n’importe laquelle valeur de m̂, on en déduit que n ∧ (E2 − E1 ) = 0,
la deuxième des relations ci-haut.
E D
ε1 > ε0 ε1 < ε0
On sait que le potentiel à l’infini a la forme Φ0 = −E0 z = −E0 r cos θ ou encore Φ0 = −E0 rP1 (cos θ).
La forme générale du potentiel extérieur Φ(ext.) est
∞
X Bl
Φ(ext.) (r) = Φ0 + Pl (cos θ) (1.158)
rl+1
l=0
On a séparé la contribution Φ0 du champ externe de celle de la sphère. Cette dernière doit être finie quand
r → 0 et nulle quand r → ∞, ce qui justifie la forme employée. Il suffit dès lors de respecter les conditions
aux limites.
1. Le potentiel doit être partout continu. Donc
Al = Bl a−(2l+1) (1.160)
∞
(int.)
X
Er = E0 P1 (cos θ) − lAl al−1 Pl (cos θ) (1.162)
l=0
(ext.) (int.)
La condition à respecter est ε0 Er = ε 1 Er . En égalisant les coefficients de Pl (cos θ) on obtient
B0 = 0
−3
ε0 (E0 + 2B1 a ) = ε1 (E0 − A1 ) (1.163)
−l−2 l−1
ε0 ((l + 1)Bl a ) = ε1 (−lAl a ) (l > 1)
36 1. Électrostatique
ε1 − ε0 ε1 − ε0 3
A1 = E0 B1 = a E0 (autres=0) (1.164)
ε1 + 2ε0 ε1 + 2ε0
ε1 − ε0 a3 E0 cos θ
Φ(ext.) (r) = Φ0 + (1.165)
ε1 + 2ε0 r2
−3ε0
Φ(int.) (r) = E0 z (1.166)
ε1 + 2ε0
Le champ à l’intérieur est donc uniforme :
3ε0
E(int.) = E0 ẑ (1.167)
ε1 + 2ε0
Remarques
1. Si ε1 > ε0 alors E (int.) < E0 : la polarisation induite à l’intérieur tend à compenser le champ. Par contre
D est plus intense à l’intérieur qu’à l’extérieur : les lignes de D sont continues.
2. Si ε1 < ε0 alors E (int.) > E0 .
3. le potentiel à l’extérieur est Φ0 plus le potentiel d’un dipôle d’intensité (si ε0 = 1)
ε1 − 1 3
|d| = a E0 (1.168)
ε1 + 2
Une connaissance de l’énergie électrostatique d’un système nous permet de calculer la force
exercée sur un conducteur ou un diélectrique en utilisant le principe du travail virtuel. Supposons
qu’un ensemble de conducteurs, séparés par un ou plusieurs milieux diélectriques, comporte une
partie mobile libre de se déplacer et dont la position est donnée par une coordonnée généralisée
x. Par exemple, x pourrait être la position d’un conducteur ou la position de l’interface entre
deux milieux diélectriques (on pense ici à la position d’un bloc de matière diélectrique libre de se
déplacer à l’intérieur d’un condensateur). L’énergie électrostatique U (x) dépend en général de ce
paramètre. On s’intéresse ici à la force F que ressent l’objet libre de se déplacer. Pour déterminer
cette force, on fait subir à l’objet un déplacement virtuel δx, compatible avec la conservation de
l’énergie. Le travail effectué par la force F lors de ce déplacement virtuel est le travail virtuel δW :
δW = F δx (1.172)
Si le système est isolé, ce travail ne peut être fourni que par l’énergie potentielle du système, qui
doit diminuer si δW est positif :
∂U
δW = −δU = − δx (1.173)
∂x
Il s’ensuit que la force généralisée est
∂U
F =− (1.174)
∂x q
où l’indice q signifie que la force est calculée pour un système isolé, dans lequel la charge sur
chacun des conducteurs est par conséquent constante.
La situation est différente si les conducteurs sont maintenus à un potentiel constant. Dans ce
cas, le système n’est pas isolé puisqu’il est relié à une source quelconque de potentiel (une pile, ou
un réservoir de charge à un potentiel fixe). Les sources de potentiel peuvent donc accomplir un
travail δWs sur le système et l’expression (1.173) doit être modifiée :
P
La variation δ Ũ à potentiel constant est précisément δU − a δqa Φa et la force ressentie par l’ojet
est alors !
∂ Ũ
F =− (1.178)
∂x
V
où l’indice V signifie que la force est calculée pour un système ouvert, dans lequel le potentiel sur
chacun des conducteurs est constant.
38 1. Électrostatique
La distinction qui existe entre U et Ũ est la même qui existe en thermodynamique entre
l’énergie interne U et l’énergie libre de Helmholtz F = U − T S: U est minimisé dans un système
isolé, alors que c’est F qui est minimisé à l’équilibre dans un système non isolé, à température
constante (T et S sont respectivement la température absolue et l’entropie). De même, à pression
constante, c’est l’enthalpie G = F − P V qui est minimisée (P et V sont la pression et le volume
respectivement).
Écrivons maintenant l’énergie électrostatique (1.171) sous une autre forme, en inversant le
calcul qui a été fait pour aboutir à (1.171):
Z
1
U =− d3 r ∇Φ · D
8π
Z Z (1.179)
1 1
=− da · DΦ + d3 r Φ∇·D
8π 8π
Dans la première équation, on intègre sur tout le volume des diélectriques, qui est limité par les
surfaces conductrices. En intégrant par parties, on obtient la deuxième équation, où la normale
de l’intégrale de surface est extérieure au diélectrique, c’est-à-dire intérieure aux conducteurs. Le
deuxième terme du membre de droite d’annule car, par hypothèse, il n’y a pas de charge libre à
l’intérieur du diélectrique. Le premier terme est intégré sur les surfaces de tous les conducteurs
en présence et, en inversant le signe de l’élément de surface, donne simplement la charge qa sur le
conducteur a fois le potentiel Φa , qui est constant sur chaque conducteur. Plus précisément,
1X
U= q Φ (1.180)
2 a a a
Par conséquent,
X 1X
Ũ = U − q a Φa = − q Φ = −U (1.181)
a
2 a a a
On peut donc exprimer la force ressentie par un objet comme
∂U
F = (1.182)
∂x V
∂U 1 ∂C
F = = V2 (V constant)
∂x 2 ∂x
(1.183)
∂U Q2 ∂C
F =− = (Q constant)
∂x 2C 2 ∂x
Comme Q = CV , on voit que la force est la même dans les deux cas. Il reste à calculer la capacité C. Pour
cela, intégrons la densité d’énergie électrostatique E · D/8π dans l’espace entre les plaques, en négligeant les
1.A. Notation indicielle 39
effets de bord. Comme E = V /d et que D = εE dans la partie du condensateur remplie par le diélectrique,
on trouve
1 V 2h
i
U= εLxd + L(` − x)d (1.184)
8π d
L h i
C= ` + (ε − 1)x (1.185)
4πd
et donc
∂C L V 2L
= (ε − 1) F = (ε − 1) > 0 (1.186)
∂x 4πd 8πd
On constate que le diélectrique est attiré par le condensateur, qu’on se place à charge constante ou à
potentiel constant.
Appendice 1.A
Notation indicielle
Cette section constitue une introduction à la notation indicielle souvent employée dans les
calculs ou démonstrations de ces notes.
Un vecteur A peut être représenté par ses composantes : A = (A1 , A2 , A3 ). On écrit la ie
composante Ai . En particulier, on emploie l’abréviation ∂i pour le gradient :
∂
∂i ≡ (1.187)
∂xi
On emploie couramment la convention de sommation, qui stipule qu’on doit sommer sur les
indices qui apparaissent en paires. Par exemple, on écrit
3
X
Ai Bi pour Ai Bi (1.188)
i=1
ou on écrit
3
X ∂Φ
Ai Ai ∂k ΦBk Bj pour Bj Ai Ai B (1.189)
∂xk k
i,k=1
Dans cette notation, l’indice sur lequel on somme (l’indice de sommation) est l’analogue d’une
variable d’intégration : le symbole précis utilisé importe peu. Ainsi, on a Ai Bi = Ak Bk .
Le symbole de Kronecker δij est défini comme suit :
(
1 si i=j
δij = (1.190)
0 si i 6= j
Il s’agit en fait des éléments de matrice de la matrice unité. On a donc δij Aj = Ai ou encore
∂i xj = δij .
40 1. Électrostatique
Le symbole (ou tenseur) de Levi-Civita εijk est défini comme suit en dimension 3: εijk = 0
sauf
ε123 = ε231 = ε312 = 1 ε213 = ε132 = ε321 = −1 (1.191)
Ce symbole est donc complètement antisymétrique dans ses indices (ex. εijk = −εjik ) et est
invariant par permutations cycliques de ses indices (εijk = εjki ). Naturellement, il s’annule si
deux des indices sont identiques. En dimension d, ce symbole comporte d indices et est défini de
manière analogue. On a la propriété suivante :
εijk εmnk = δim δjn − δin δjm (1.192)
Appendice 1.B
Polynômes de Legendre
Les polynômes de Legendre sont essentiels à la formulation du développement multipolaire.
Cette annexe est une introduction à ces fonctions spéciales particulièrement simples.
On peut définir les polynômes de Legendre par la fonction génératrice
∞
1 X
G(x, t) = √ = xl Pl (t) (1.196)
1 − 2xt + x2 l=0
Les polynômes Pl (t) sont obtenus en effectuant correctement le développement binomial. Les
premiers polynômes sont
P0 (t) = 1
P1 (t) = t
1 (1.197)
P2 (t) = (3t2 − 1)
2
1
P3 (t) = (5t3 − 3t)
2
1.B. Polynômes de Legendre 41
P0
P1
P3
P4
—1 +1 t
P2
—1
En d’autres termes, les polynômes d’ordre pair sont pairs et ceux d’ordre impair sont impairs.
Deuxièmement, on trouve facilement la valeur des polynômes à t = 1. Pour cette valeur, la
fonction génératrice est
1
G(x, 1) =
1−x
X
= xl (série géométrique)
l=0
(1.201)
∞
X
= xl Pl (1)
l=0
La série géométrique ne converge que si −1 < x < 1. La série (1.196) possède le même domaine
de convergence.
En pratique, le calcul des polynômes de Legendre s’effectue à l’aide de relations de récurrence,
comme celle-ci :
(2l + 1)tPl (t) = (l + 1)Pl+1 (t) + lPl−1 (t) (1.203)
Cette relation permet de déterminer tous les polynômes de Legendre en connaissant les expres-
sions pour P0 et P1 . Elle se démontre de la manière suivante : on calcule la dérivée de la fonction
génératrice par rapport à x. On trouve que
∂G(x, t)
(1 − 2xt + x2 ) + (x − t)G(x, t) = 0 (1.204)
∂x
Dans cette relation on substitue le développement en série (1.196) et on obtient une identité sous
la forme d’une développement en série de puissances, puisque les facteurs (1 − 2xt + x2 ) et x − t
sont eux-mêmes des polynômes en x et t. En appliquant l’égalité pour chaque puissance de x, on
trouve la relation (1.203).
On démontre aussi la formule suivante, dite formule de Rodrigues:
1 dl 2
Pl (t) = (t − 1)l (1.205)
2l l! dtl
Cette formule se démontre par calcul direct, en comparant avec l’Éq. (1.198). Elle nous permet de
démontrer la relation d’orthogonalité suivante :
Z 1
2
dt Pm (t)Pn (t) = δmn (1.206)
−1 2l + 1
Il suffit pour cela d’appliquer la formule de Rodrigues dans l’intégrale et d’intégrer par parties
complètement l’un des deux facteurs. La normalisation du cas m = n s’obtient plutôt en intégrant
le carré de la formule (1.196) et en utilisant l’orthogonalité.
La dernière propriété des polynômes de Legendre que nous allons démontrer est que ce sont
des solutions de l’équation différentielle de Legendre :
d 2 dPl
(1 − t ) + l(l + 1)Pl = 0 (1.207)
dt dt
Pour démontrer ceci, considérons une charge unité située à la position r0 et considérons le potentiel
électrique au point r créé par cette charge :
1
Φ(r) = (1.208)
|r − r0 |
Supposons que |r0 | < |r| et choisissons nos axes de telle sorte que r0 soit sur l’axe des z. Définissons
r = |r| et r0 = |r0 |. L’angle entre les vecteurs r et r0 est la coordonnée angulaire θ des coordonnées
sphériques. Remplaçons cette coordonnée par t = cos θ. On peut alors utiliser la série (1.196) pour
décrire le potentiel :
1 1 1
0
= p
|r − r | r 1 − 2tr0 /r + (r0 /r)2
∞ (1.209)
X r0l
= P (t)
rl+1 l
l=0
1.B. Polynômes de Legendre 43
∂2Φ
2 1 ∂ 2 ∂Φ 1 ∂ ∂Φ 1
∇ Φ= 2 r + 2 sin θ + (1.210)
r ∂r ∂r r sin θ ∂θ ∂θ r2 sin2 θ ∂ϕ2
X r0l d
2 dPl
l(l + 1)Pl (t) + (1 − t ) =0 (1.212)
rl+1 dt dt
l
Comme il s’agit d’une identité, valable pour tout r > r0 et pour toute valeur de t, on conclut
que les polynômes de Legendre satisfont à l’équation différentielle (1.207). Les solutions de cette
équation linéaire du deuxième ordre sont uniques si on spécifie deux conditions initiales ou
deux conditions aux limites. Par exemple, les polynômes de Legendre sont les solutions uniques
satisfaisant aux conditions Pl (1) = 1 et Pl (−1) = (−1)l .
Problèmes
Problème 1.1
Un disque de rayon a porte une densité superficielle de charge constante ρs . L’axe des z est perpendiculaire
à ce disque et passe par son centre. Calculez le potentiel électrique Φ(z) le long de l’axe des z. Calculez le
champ électrique E dans les deux limites z a et z a et montrez que les résultats obtenus sont ceux
auxquels on doit s’attendre respectivement d’une charge ponctuelle et d’un plan infini.
Problème 1.2
Une sphère isolante de rayon a porte une charge q répartie uniformément dans son volume. Cette sphère est
entourée d’une enveloppe conductrice sphérique de rayon intérieur a et de rayon extérieur b > a. Trouvez le
champ électrique et le potentiel Φ partout et tracez leur profil en fonction de r.
Problème 1.3
Considérez deux charges q1 et q2 situées respectivement à l’origine et à r = (0, 0, a). Soit E le champ
produit par ces deux charges. À partir de l’expression (1.40) pour la densité d’énergie E, montrez que
l’énergie d’interaction entre les deux charges est U = q1 q2 /a. Indice : écrivez E = E1 + E2 , où E1 et E2
sont respectivement produits par les charges 1 et 2, et ne considérez dans la densité E que le terme qui
provient de la présence simultanée des deux charges.
Problème 1.4
Soit une particule ponctuelle de charge q ayant une trajectoire r(t). La densité de courant et la densité de
charge associées sont J(r, t) = qvδ(r − r(t)) et ρ(r, t) = qδ(r − r(t)). Démontrez explicitement que l’équation
de continuité ∇ · J + ∂ρ/∂t = 0 est satisfaite.
44 1. Électrostatique
Problème 1.5
Considérez l’expression J = ρv pour la densité de courant. Démontrez explicitement que l’équation de
continuité est satisfaite pour cette expression. Considérez pour cela une région V , délimité par une surface
S. Si on suit le flot déterminé par la vitesse v(r), cette région est transportée en un temps δt en une
autre région V 0 (de volume généralement différent). La charge totale à l’intérieur de cette région doit être
inchangée. Le concept d’élément de volume en co-mouvement avec v est ici utile.
Problème 1.6
Un filament chauffant émet des électrons en son extrémité, située à l’origine. On suppose que les électrons
s’échappent de cette extrémité également dans toutes les directions et qu’ils forment un nuage électronique
de densité de charge −ρ(r) autour de l’origine (ρ est positif). De plus, les électrons situés à une distance
r de l’origine s’éloignent de celle-ci, en raison de la répulsion électrostatique des autres électrons, avec une
vitesse v(r). Plaçons-nous en situation de régime permanent, c’est-à-dire supposons que la densité ρ(r) et
la vitesse v(r) ne dépendent pas du temps.
a) Montrez que ρ(r)v(r) = I/(4πr2 ), où I est le courant émis par le filament, en valeur absolue.
b) Montrez que le potentiel électrique Φ(r) satisfait à l’équation différentielle suivante :
r
0 2 00 m
2rΦ + r Φ = I
2eΦ
Indice : comme dans l’exemple 1.1 des notes, on suppose que l’Éq. (1.26) est valable (ceci fixe la valeur du
potentiel à zéro à l’origine, en supposant que les électrons sont émis avec une vitesse négligeable).
Problème 1.7
Supposons que la loi de Coulomb est légèrement modifiée par un facteur exponentiel provenant d’une masse
non nulle du photon, de sorte que le potentiel électrique produit par une particule de charge q située à
l’origine est Φ(r) = qe−µr /r, où µ est une constante représentant l’inverse de la portée de l’interaction
électromagnétique.
a) Montrez que l’équation de Poisson pour le potentiel serait modifiée pour devenir l’équation de
Helmholtz :
(∇2 − µ2 )Φ = −4πρ
Commencez par le démontrer pour une charge ponctuelle située à l’origine et généralisez le résultat en vous
servant du principe de superposition.
b) La masse du photon associée à la constante µ est m = hµ/c, où h est la constante de Planck et c la
vitesse de la lumière. Sachant que la limite supérieure pour la masse du photon est de 4.10−48 g, à quelle
portée minimale µ−1 cela correspond-il? Quelle est la signification pratique de cette distance?
Problème 1.8
Étant donné un ensemble de N conducteurs portés à des valeurs Φi (i = 1, 2 . . . N ) du potentiel
électrostatique et supportant des charges Qi , montrez, à partir de l’expression E2 /8π pour la densité
d’énergie électrostatique, que l’énergie U du système est minimisée quand le potentiel satisfait à l’équation
de Laplace ∇2 Φ = 0 dans les régions situées entre les conducteurs, et que cette valeur minimum est donnée
par
1X
U= Qi Φi (1.213)
2
i
Problème 1.9
Étant donnés deux conducteurs portant des charges égales et opposées ±Q, on définit la capacité comme
le rapport Q/V où V est la différence de potentiel entre les deux conducteurs. Exprimer C en fonction des
coefficients de capacité Cij (i, j = 1, 2).
1.B. Polynômes de Legendre 45
Problème 1.10
Soit deux plaques conductrices carrées de côté l, parallèles et séparées par une distance d telle que l d.
a) Montrez que la capacité est C = l2 /4πd.
b) Si les plaques ne sont pas parfaitement parallèles, de telle façon que la distance entre les plaques soit
d + δ sur un bord et d − δ sur l’autre bord, avec d δ, montrez que la capacité est modifiée de la manière
suivante, au premier ordre en δ:
l2 δ2
C= 1+ 2 (1.214)
4πd 3d
Problème 1.11
Trois sphères conductrices identiques sont placées aux sommets d’un triangle équilatéral. Si la première
sphère est maintenue à un potentiel V et les deux autres à un potentiel nul, les charges sur les sphères sont
respectivement q, q 0 et q 0 . Montrez que si les trois sphères étaient à un potentiel commun W , chaque sphère
porterait une charge (2q 0 + q)W/V . Calculez les potentiels des trois sphères quand les charges portées par
les sphères sont respectivement Q, 0 et 0. Cette question est tirée de l’examen de l’A.C.P. (1993).
Problème 1.12
Calculer le moment quadrupolaire Q des distributions de charge suivantes (prenez l’origine au centre de
l’objet):
a) Une tige très mince de longueur l, chargée uniformément avec charge totale q.
b) Un disque très mince de rayon r, chargé uniformément avec charge totale q.
Problème 1.13
Exprimez en fonction des moments Qij la contribution quadrupolaire au couple N ressenti par une
distribution de charge placée dans un champ électrique non uniforme.
Problème 1.14
Montrez que, dans le vide, la moyenne du potentiel Φ sur la surface d’une sphère est égale au potentiel au
centre de la sphère.
Problème 1.15
Étant donné un plan conducteur infini situé à z = 0 et un dipôle ponctuel d = pẑ situé à r = (0, 0, a),
calculez la force ressentie par ce dipôle.
Problème 1.16
Une tige conductrice cylindrique de rayon a et de longueur quasi-infinie est placée parallèlement à l’axe des
z dans un champ électrique uniforme qui lui est perpendiculaire (E = E0 r̂). Calculez le potentiel électrique
en dehors de la tige. Exprimez la réponse en coordonnées cylindriques. Calculez aussi la densité superficielle
de charge induite sur la surface de la tige.
Problème 1.17
Considérez maintenant deux tiges identiques, comme dans le problème précédent, parallèles l’une à l’autre
et séparées par une distance d, telle que d a. Les deux tiges sont encore placées dans un champ électrique
uniforme qui leur est perpendiculaire, comme ci-haut. Calculez l’énergie d’interaction par unité de longueur
des deux tiges (Supposez, en première approximation, que la présence de la première tige n’affecte pas la
distribution de charge sur la deuxième tige et vice-versa. Limitez vos calculs au premier ordre non-trivial
en a/d). Montrez que la force par unité de longueur entre les deux tiges varie comme ∼ 1/d3 .
46 1. Électrostatique
Problème 1.18
Une sphère conductrice neutre est placée dans un champ électrique uniforme. Calculez la densité superfi-
cielle de charge ρs induite, partout sur la surface de la sphère.
Problème 1.19
Une charge ponctuelle q est placée à l’intérieur d’une sphère conductrice neutre de rayon b, à une distance
a du centre (voir l’exercice résolu ci-haut). Quelle est la force exercée sur cette charge en fonction de b?
Peut-on remplacer la charge induite sur la sphère par un ensemble d’images?
Problème 1.20
Une charge ponctuelle q est placée à l’extérieur d’une sphère conductrice neutre de rayon b, à une distance
a du centre (a > b). Trouver le potentiel Φ partout à l’extérieur de la sphère. Quel est le moment dipolaire
induit sur la sphère?
Problème 1.21
Une coquille conductrice de rayon a est sciée en deux hémisphères qui sont ensuite réunis par un joint
isolant d’épaisseur négligeable à l’équateur. Les hémisphères nord et sud sont maintenus respectivement à
des potentiels +V et −V . Utilisez les coordonnées sphériques où l’axe des z coı̈ncide avec l’axe polaire
de la sphère. Exprimez le potentiel Φ(r) partout à l’intérieur et à l’extérieur de la coquille, comme une
série en polynômes de Legendre. Indiquez lesquels parmi les coefficients sont nuls et calculez les 2 premiers
coefficients non nuls.
z
—V
Problème 1.22
En démontrant l’équation (1.86) pour la force exercée sur une distribution de charge, pourquoi néglige-t-on
l’effet du champ électrique créé par la distribution elle-même?
Problème 1.23
On veut ici démontrer rigoureusement que le potentiel à l’extérieur d’une sphère conductrice creuse de rayon
a contenant une charge q est Φ = q/r, quelle que soit la façon dont cette charge est répartie à l’intérieur
de la sphère. À cette fin, nous allons procéder à une inversion de l’espace, c.-à-d. à une transformation
conforme qui transforme tout point à l’intérieur de la sphère en un point à l’extérieur et vice-versa.
a) Considérez la transformation r → (a2 /r2 )r ou, en coordonnées sphériques, (r, θ, ϕ) → (a2 /r, θ, ϕ).
Définissons un potentiel transformé
Φ0 (r, θ, ϕ) = rα Φ(a2 /r, θ, ϕ)
Démontrez que si Φ satisfait à l’équation de Laplace, alors Φ0 aussi, pourvu que l’exposant α soit bien
choisi. Quelle est la valeur de α appropriée?
1.B. Polynômes de Legendre 47
b) À l’aide de cette transformation, montrez comment la solution Φ(r, θ, ϕ) = q/r à l’extérieur de la sphère
conductrice se démontre à partir de la solution Φ0 (r, θ, ϕ) = const. à l’intérieur d’une sphère conductrice ne
contenant aucune charge.
c) Refaites la partie (a), cette fois en dimension 2, en coordonnées polaires planes. Quelle est la valeur de
α cette fois-ci?
Problème 1.24
Montrez que le champ électrique à l’intérieur d’une sphère uniformément polarisée est uniforme et donné
par E = −(4π/3)P. Montrez ensuite que le champ électrique Es existant dans une cavité sphérique au sein
d’un milieu diélectrique uniformément polarisé est relié au champ E dans le milieu par
4π
Es = E + P (1.215)
3
Indice : considérez deux sphères chargées uniformément (avec des charges égales mais opposées) superposées
avec un déplacement relatif.
Problème 1.25
Considérez une coquille sphérique de rayon intérieur a et de rayon extérieur b (b > a). Cette coquille est
faite d’un matériau diélectrique linéaire de constante ε et est plongée dans un champ électrique externe
uniforme E = E0 ẑ. Le milieu à l’extérieur et à l’intérieur de la coquille est de constante diélectrique ε = 1.
Calculez E, D et P partout.
Problème 1.26
Un condensateur est formé de deux cylindres métalliques concentriques de rayons a et b respectivement
(b > a). On peut utiliser un tel condensateur pour mesurer la susceptibilité électrique d’un liquide, comme
suit. On applique une différence de potentiel V constante aux bornes de ce condensateur (c.-à-d. entre
les deux cylindres) et on le plonge dans le liquide. Dans le but de minimiser l’énergie totale du système,
sous la contrainte que le potentiel est constant, une certaine quantité de liquide diélectrique monte entre
les cylindres et forme une colonne de hauteur h (voir figure). La hauteur de la colonne est déterminée en
équilibrant les forces électriques et gravitationnelles.
a) À l’aide de la loi de Gauss, montrez que le champ électrique entre les cylindres est relié à la différence
de potentiel de la manière suivante :
V
E=
r ln(b/a)
où χe est la susceptibilité électrique. À partir de cette expression, trouvez la force verticale exercée sur la
colonne de liquide, en n’oubliant pas que la différence de potentiel V est constante.
48 1. Électrostatique
c) Quelle est la force gravitationnelle s’exerçant sur la colonne de liquide? Obtenez une expression pour la
susceptibilité χe en équilibrant les forces gravitationnelle et électrique.
h b
b
a a
Magnétostatique
1. Étant donné les unités utilisées pour la charge et fixées par la forme que nous avons choisi
pour la loi de Coulomb, la constante de proportionnalité c a les dimensions (unités) de la
vitesse. D’après les mesures expérimentales elle coı̈ncide avec la vitesse de la lumière dans le
vide2 :
c = 2.99792 . . . × 1010 cm/s
2. La ressemblance avec la loi de Coulomb est frappante. Cependant les éléments de circuits
n’existent pas de façon isolée! Pour deux circuits fermés, la force totale diminue plus rapide-
ment avec la distance. Cette force est alors une double intégrale le long des circuits :
dl1 ∧ (dl2 ∧ r̂12 )
I I
I1 I2
F12 = 2 2
c 1 2 r12
1
la porte n’est pas irrémédiablement fermée, mais ces monopôles seraient alors de rares accidents issus de
l’univers primordial et très massifs.
2
Sans même avoir touché à l’électrodynamique, il est dès lors manifeste que la lumière est probablement un
phénomène électromagnétique. Ce fait échappa cependant longtemps aux pionniers de l’électromagnétisme,
en raison sans doute des systèmes d’unités boiteux alors en usage!
50 2. Magnétostatique
Cependant,
dl1 ∧ (dl2 ∧ r̂12 ) dl (dl · r̂ ) r̂ (dl · dl2 )
2 = 2 12 12 − 12 12
r12 r12 r12
r̂ (dl · dl2 )
= −dl2 (dl1 · ∇(1/r12 )) − 12 12
r12
L’intégrale du premier terme est nulle, car la circulation d’un gradient sur un contour fermé
s’annule. On peut donc écrire
dl1 · dl2
I I
II
F12 = − 122 2 r̂12
c 1 2 r12
J1 (r0 ) ∧ (r − r0 )
Z
1
B(r) = d3 r0 (2.4)
c |r − r0 |3
Z
1
F= d3 r J(r) ∧ B(r) (2.5)
c
où F est la force totale exercée sur la distribution de courant J2 par le champ magnétique B créé
par une autre distribution de courant J1 . Il est entendu que les distributions de courant J1 et
J2 obéissent à la version stationnaire de l’équation de continuité: ∇·J = 0. Cette généralisation
provient du fait qu’on peut toujours considérer une distribution de courant sans divergence comme
une assemblée continue de circuits : considérons les lignes de champ d’une distribution J, c’est-à-
dire l’ensemble des courbes auxquelles J est tangent en chaque point. Ces courbes sont fermées
du fait que ∇·J = 0. En prenant un ensemble de telles courbes proches les unes des autres (un
faisceau de section da) on obtient l’équivalent d’un circuit fermé dans lequel circule un courant
dI = |J|da. On peut donc remplacer l’intégrale de chemin sur une circuit par une intégrale de
volume : Z Z
I dl(· · ·) −→ d3 r J(· · ·) (2.6)
C
3
La convention arbitraire pour le sens du produit vectoriel affecte le signe de B, mais non le résultat physique
(la force), car cette convention est appliquée deux fois en succession.
2.1. Loi d’Ampère 51
B = ∇∧A (2.7)
∇·B = 0 (2.10)
4π
∇∧B = J (2.12)
c
Comme pour la loi de Gauss avec la loi de Coulomb, on peut démontrer la loi d’Ampère à
partir de celle de Biot et Savart. On écrit tout d’abord
∇∧B = ∇(∇·A) − ∇2 A (2.13)
1. La quantité ∇·A s’annule puisque la densité de courant J est stationnaire (par hypothèse):
Z
1 1
∇·A(r) = d3 r0 J(r0 ) · ∇
c |r − r0 |
Z
1 1
=− d3 r0 J(r0 ) · ∇0 (2.14)
c |r − r0 |
∇0 · J(r0 )
Z
1
= d3 r0
c |r − r0 |
52 2. Magnétostatique
(∇0 est le gradient par rapport à r0 ). En passant de la deuxième à la troisième équation, nous
avons procédé à une intégration par partie dont le terme de surface s’annule, puisque par
hypothèse la distribution de courant tend vers zéro suffisamment rapidement à l’infini. Cette
expression est nulle puisque ∇·J = 0.
2. La similitude entre (2.8) (composante par composante) et l’expression correspondante pour
le potentiel électrostatique [éq. (1.17)] nous mène à l’équation de Poisson pour le potentiel
vecteur :
4π 4π
∇2 Ai = − Ji ou ∇2 A = − J (2.15)
c c
La loi d’Ampère est donc démontrée par ces deux remarques.
4πI 2I
2πrB(r) = =⇒ B(r) = (r > a) (2.16)
c cr
Si, au contraire, r < a, le courant traversant la boucle est I(r) = Ir2 /a2 et l’induction magnétique est alors
2Ir
B(r) = (r < a) (2.17)
ca2
On peut aussi résoudre ce problème à l’aide du potentiel vecteur, par analogie avec une situation
électrostatique. Si nous avions un cylindre infini uniformément chargé (charge électrique par unité de
longueur ρl ), alors le théorème de Gauss nous permettrait d’obtenir l’expression suivante pour le potentiel
électrique :
−2ρl ln(r/a) (r > a)
Φ(r) = 2 2 (2.18)
−ρl r − a (r < a)
a2
où nous avons ajusté les constantes afin que Φ soit continu et nul à r = a. La similitude entre les éq.
(2.15) et (1.18) nous permet de traduire ce résultat en magnétostatique. L’équivalent de ρl est ici I/c et le
potentiel vecteur créé par ce courant est alors
2I
− ln(r/a)ẑ (r > a)
c
A(r) = (2.19)
I r2 − a2
−
ẑ (r < a)
c a2
En substituant cette expression dans la relation (2.5) pour la force ressentie par la distribution,
on trouve Xq
i
F= vi ∧ B(ri ) (2.22)
i
c
Il s’agit en fait d’un résultat qui va au-delà de l’équation (2.5), car il est valide même quand le
champ magnétique est non-stationnaire. Il est en effet difficile d’imaginer que la force exercée sur
une seule charge par un champ magnétique puisse dépendre de la présence des autres charges
dans la distribution de courant. On conclut donc que la force électromagnétique totale exercée sur
une charge ponctuelle en mouvement est
qv
F = qE + ∧B (2.23)
c
1 q
L= mv2 − qΦ(r) + A(r) · v (2.24)
2 c
Pour établir que ce lagrangien est bien adéquat, il suffit de démontrer que les équations de
Lagrange coincident avec la formule (2.23) ci-haut. Les équations de Lagrange sont
d ∂L ∂L
− =0
dt ∂v ∂r
d h q i q
= mv + A(r) + q∇Φ − ∇(Ai (r)vi ) (2.25)
dt c c
q dxi q
= ma + ∂i A − qE − − ∇(Ai (r)vi )
c dt c
Remarquons que
∂i Avi − ∇(Ai (r)vi ) = −v ∧ (∇∧A) = −v ∧ B (2.26)
La force de Lorentz (2.23) s’ensuit immédiatement.
Passons maintenant à l’hamiltonien. Le moment conjugué p est
∂L q
p= = mv + A(r) (2.27)
∂v c
54 2. Magnétostatique
puisque ∇·J = 0.
Le deuxième terme (le dipôle) est
Z Z
(2) 1 3 0 0 0 1
Ak = 2 d r Jk (r )r cos γ = 3 d3 r0 Jk (r0 )x0i xi (2.33)
cr cr
Notons en passant que
∂m (xi xk Jm ) = xk Ji + xi Jk (2.34)
Le membre de gauche est une divergence et donc on obtient, en intégrant,
Z
d3 r0 (x0i Jk + x0k Ji ) = 0 (2.35)
m∧r
A(2) = (2.40)
r3
3n(n · m) − m
B(2) = (2.43)
r3
Remarques :
1. La forme du champ est la même que pour un dipôle électrique.
2. Les unités du moment dipolaire magnétique sont les mêmes que celles du moment dipolaire
électrique. Le facteur de 1/c a comme conséquence que le moment dipolaire magnétique
typique sera beaucoup plus petit que son analogue électrique, à moins que les vitesses typiques
de la distribution de courant soient voisines de c.
56 2. Magnétostatique
3. Cette forme est valide loin des sources et constitue le terme dominant dans le développement
multipolaire magnétique. Une approximation supplémentaire consiste R 3 à ajouter un terme
à l’origine des coordonnées qui assure que l’intégrale du champ d r B soit reproduite
correctement :
3n(n · m) − m 8π
B= + mδ(r) (2.44)
r3 3
Le dernier terme ne contribue qu’à l’origine (il n’est donc pas partie du développement
multipolaire) et joue un rôle important dans le traitement de la structure hyperfine de
l’hydrogène.
Selon la correspondence (2.6) et la définition (2.39), le moment dipolaire d’un circuit de
courant est Z
I
m= r ∧ dl (2.45)
2c
Pour une boucle plane avec origine à l’intérieur de la boucle, il est clair que 12 r ∧ dl est l’élément
d’aire de la surface entourée par la boucle. On peut alors écrire
I
|m| = × Aire de la boucle (2.46)
c
1 X X q
i
m= qi (ri ∧ vi ) = Li (2.48)
2c i i
2M ic
où Li et Mi sont le moment cinétique et la masse de la charge i. Pour un rapport q/M commun,
on peut alors écrire
q
m= L (2.49)
2M c
en fonction du moment cinétique total L. En général, la constante de proportionnalité entre le
moment magnétique d’un système et son moment cinétique est appelée rapport gyromagnétique et
noté γ. Dans le cas d’un système de particules en mouvement ayant le même rapport q/M , on a
q
γ= (2.50)
2M c
Si le moment cinétique d’un atome ne provenait que des moments cinétique orbitaux des
électrons, la formule (2.49) pourrait s’appliquer à un atome dans son entier et, comme le moment
cinétique orbital est quantifié en multiples entiers de h̄, le moment magnétique d’un atome serait
un multiple entier du magnéton de Bohr µB :
eh̄
µB = (2.53)
2me c
−e
m= S (2.54)
me c
où S est l’opérateur du spin. L’addition de tous les moments cinétiques orbitaux et de tous les
spins dans un atome à plusieurs électrons conduit à la forme générale suivante pour le moment
magnétique de l’atome :
ge
m = γJ γ≡ (2.55)
2me c
où J est le moment cinétique total de l’atome, γ est le facteur gyromagnétique de l’atome et g
un facteur numérique appelé facteur de Landé, généralement différent de l’unité.4 Le facteur de
Landé ge d’un électron est égal à -2, selon la théorie de Dirac.
4
Il ne faut pas voir la relation (2.55) comme une relation entre des opérateurs, car le facteur de Landé
change d’un niveau d’énergie à un autre : il dépend des nombres quantiques principal et orbital. Le ‘facteur
de Landé’ d’un atome est en fait le facteur g associé à l’état fondamental à plusieurs électrons de l’atome.
58 2. Magnétostatique
La mesure du facteur ge d’un électron isolé peut être effectuée par une expérience de résonance
magnétique sur un faisceau d’électrons. La technique de résonance est très précise et produit le
résultat suivant :
ge = −2(1 + a) aexp. = 0, 001 596 521 884(43) (2.56)
(les chiffres entre parenthèses sont l’incertitude sur les deux derniers chiffres significatifs). La
théorie quantique de l’électrodynamique permet de calculer des corrections à g par rapport à sa
valeur dans la théorie de Dirac. On doit pour cela utiliser la théorie des perturbations. Le résultat
au quatrième ordre est
ath. = 0, 001 596 521 92(17)(73)(108) (2.57)
L’accord entre théorie et expérience est le meilleur de toute la physique – 9 chiffres après la
virgule!
Le proton, tout en étant aussi une particule de spin 12 , possède un facteur g bien différent :
gp = 2 × 2, 79, alors que le neutron, particule neutre de spin 21 , possède néanmoins un moment
magnétique, avec gn = −2 × 1, 93. Ces valeurs étranges proviennent du fait que le proton et le
neutron sont des particules composites formées de quarks, alors que l’électron est une particule
ponctuelle.
On conclut que l’énergie potentielle d’un dipôle dans un champ magnétique est
W = −m · B (2.63)
Ceci signifie aussi que même un champ uniforme peut exercer une influence sur un dipôle,
en exerçant un couple sur la distribution de courant de façon à aligner le moment dipolaire
avec le champ pour minimiser W . Si θ est l’angle que fait m avec B, l’énergie potentielle est
W = Bm cos θ; le couple est alors N = −∂W/∂θ ou N = Bm sin θ; sa direction est perpendiculaire
au plan qui contient m et B. En tenant compte de son signe conventionnel, on peut donc écrire
N=m∧B (2.64)
2.3. Dipôles magnétiques 59
On peut aussi démontrer cette relation directement : la composante Ni du couple exercé sur la
distribution de courant est
Z
1
Ni = imn njk d3 r xm Jj (r)Bk (r) (2.65)
c
1 8π
Udip. = (m1 · m2 − 3(n · m1 )(n · m2 )) − m · m2 δ(r1 − r2 ) (2.67)
r 3 3 1
Cette interaction dipolaire joue un rôle souvent important dans les solides magnétiques, en par-
ticulier dans les ferroaimants. Elle est aussi la source de la structure hyperfine de l’hydrogène et
d’autres atomes.
dm
= γm ∧ B (2.68)
dt
De cette équation on déduit que la composante de m dans la direction de B est constante
dans le temps, alors que les deux autres composantes forment un vecteur tournant dans le plan
perpendiculaire à B. En effet, rappelons que la dérivée temporelle d’un vecteur A tournant à une
vitesse angulaire constante ω autour d’un axe n est
dA
=ω∧A ω = ωn (2.69)
dt
60 2. Magnétostatique
Cette équation reste inchangée si on ajoute à A une composante constante dans la direction de
ω. On conclut que le moment magnétique précesse autour de B avec un vecteur vitesse angulaire
En pratique, B0 est créé par un électroaimant, souvent supraconducteur, alors que B1 est créé par
une modeste solénoı̈de ou par une cavité résonante (B1 B0 ). Comme il est difficile d’analyser le
mouvement d’un moment magnétique dans un champ tournant, plaçons-nous dans un référentiel
tournant à la même vitesse que B1 . On sait que la dérivée temporelle d’un vecteur A dans un
référentiel tournant, en fonction de sa dérivée dans un référentiel inertiel, est
dA dA
= −ω∧A (2.72)
dt R dt I
où B = B0 + B1 . Dans le référentiel tournant, on constate que la vitesse angulaire ω joue le rôle
d’un champ magnétique ‘fictif’, contribuant à un champ magnétique effectif
J(r0 ) cM(r0 ) ∧ (r − r0 )
Z
1 3 0
A(r) = dr + (2.75)
c V |r − r0 | |r − r0 |3
∇0 ∧ M(r0 ) M(r0 )
Z Z I
3 0 0 0 1 3 0
d r M(r ) ∧ ∇ = dr + da ∧ (2.76)
V |r − r0 | V |r − r0 | S |r − r0 |
après une intégration par partie (attention au signe!). Cette dernière équation a l’interprétation
suivante : L’aimantation agit comme une densité de courant J0 = c∇∧M liée au matériau, plus
une densité superficielle de courant K0 cM ∧ n apparaissant à la surface du matériau :
J0 = c∇∧M K0 = cM ∧ n (2.77)
4π
∇∧B = J + 4π∇∧M (2.78)
c
62 2. Magnétostatique
Le rôle de ce nouveau vecteur H est le même que celui de D en électrostatique, c’est-à-dire d’être
insensible aux courants liés (aux charges liées, dans le cas de D) sauf par l’intermédiaire des
conditions aux limites. Il ne faut pas chercher dans H une signification à caractère plus intuitif
que ce qui est fourni par sa définition et par la loi d’Ampère.*
Profitons de ce moment pour écrire toutes ensemble les lois de l’électromagnétisme statique :
2.4.2 Diamagnétisme
Les différents matériaux se distinguent par la relation qui existe entre M et B ou entre H
et B. Cette relation s’obtient par des considérations microscopiques impliquant évidemment la
mécanique quantique et ne peut être comprise par des arguments classiques seulement.
Conidérons premièrement les matériaux dits diamagnétiques. Dans ceux-ci les atomes ne
possèdent pas de moment magnétique permanent, car les couches électroniques sont complètement
remplies : le moment cinétique de l’atome dû aux électrons d’annule, ainsi que le moment
magnétique associé. Cependant, ces atomes peuvent développer un moment magnétique induit
si un champ magnétique externe B0 est appliqué. Contrairement au phénomène diélectrique
correspondant, le moment induit sera ici opposé au champ appliqué. Ceci peut se comprendre
intuitivement à l’aide de la loi de Lenz pour l’induction magnétique (voir le chapitre suivant):
l’application du champ magnétique externe induit un courant dans l’atome qui en retour crée un
champ induit qui s’oppose au changement du champ, c’est-à-dire qui s’oppose au champ externe.
Donc le moment magnétique associé à cette boucle de courant atomique est opposé au champ.
l’intensité de ce phénomène est très faible et proportionnelle au champ appliqué: on a donc une
relation linéaire entre l’aimantation induite M et le champ magnétique B ou H. On écrit
B = µH ou M = χm H (2.82)
* Notons que la terminologie utilisée peut varier d’un auteur à l’autre. On appelle parfois H
l’intensité du champ magnétique, alors que B est appelé tout bonnement champ magnétique. C’est
le cas de Marion et de Landau et Lifchitz.
2.4. Magnétisme macroscopique 63
2.4.3 Paramagnétisme
Si un matériau est composé d’atomes dont les couches électroniques ne sont que partiellement
remplies, ce qui est le cas notamment des éléments dits de transition, alors ces atomes portent un
moment magnétique permanent. Cependant, ces moments permanents sont généralement orientés
de manière aléatoire en raison de l’agitation thermique, de sorte que l’aimantation M est en
moyenne zéro en l’absence de champ appliqué. De tels matériaux sont dits paramagnétiques.
Si un champ externe B0 est appliqué au matériau, les moments magnétiques microscopiques
auront tendance à s’aligner avec le champ B dans le but de minimiser leur énergie potentielle
−m·B. L’aimantation sera alors proportionnelle au champ B, pourvu que cette énergie potentielle
magnétique soit petite par rapport à l’énergie thermique kB T . On peut alors écrire une relation
linéaire comme ci-haut, sauf que maintenant χat. est positif et µ > 1. Le calcul simple précédant
l’éq. (1.146) peut être repris ici sans modification et donne
m20
hmi = χat. B χat. = (2.87)
3kB T
où m0 est la grandeur du moment dipolaire magnétique permanent. On constate que la suscep-
tibilité magnétique χm varie en fonction inverse de la température (loi de Curie). L’ordre de
grandeur typique de m0 est eh̄/me c. Donc
e2 h̄2 1
2
2 2 e /a0
χat. ∼ 2 2 = a0 α (2.88)
me c kB T kB T
où e2 /a0 est l’énergie d’ionisation Eion. typique d’un atome. La susceptibilité paramagnétique est
donc plus grande que la susceptibilité diamagnétique par un facteur Eion. /kB T . À température
ambiente, ce facteur est de l’ordre de 102 – 103 . La même remarque s’applique au rapport de la
polarisabilité électronique à la polarisabilité paraélectrique. Les valeurs typiques de χm varient
entre 10−5 et 10−2 .
64 2. Magnétostatique
2.4.4 Ferromagnétisme
Lorsque la température d’un matériau paramagnétique est suffisamment faible, il se pro-
duit généralement un changement de phase qui affecte beaucoup les propriétés magnétiques du
matériau. Dans les matériaux dits ferromagnétiques, les moments magnétiques permanents des
atomes s’alignent les uns sur les autres de manière spontanée, en l’absence de champ externe. Par
exemple, le fer, le nickel et le cobalt sont ferromagnétiques à température ambiente. Le change-
ment de phase ferromagnétique-paramagnétique se produit à la température dite de Curie, notée
TC et qui dépend des propriétés microscopiques de chaque substance.
L’origine du phénomène est liée à l’interaction électrostatique entre les électrons et au principe
d’exclusion de Pauli. Rappelons que ce principe stipule que la fonction d’onde d’un système de plusieurs
électrons doit être antisymétrique quand deux électrons sont échangés. Considérons par exemple un cristal
de Nickel, dont chaque atome comporte un électron dans une couche supérieure. La fonction d’onde de deux
électrons appartenant à des atomes voisins de Nickel doit être antisymétrique. Cette antisymétrie peut être
réalisée soit dans la partie spatiale de la fonction d’onde, soit dans la partie ‘spin’. Étant donné la forte
répulsion électrique entre les deux électrons, une fonction d’onde spatialement antisymétrique garantit que
les deux électrons ne seront pas au même endroit en même temps et donc minimise l’énergie. Il s’ensuit
que la partie ‘spin’ de la fonction d’onde doit être symétrique et que les spins des deux électrons sont
parallèles. Donc il y a alignement des moments magnétiques dans ce matériau et l’énergie associée à cet
alignement est d’origine électrique, donc beaucoup plus grande qu’une énergie magnétique typique, comme
l’interaction dipolaire par exemple. L’argument donné ici est simpliste et ne s’applique pas de manière
universelle, sinon tous les matériaux paramagnétiques deviendraient ferromagnétiques à des températures
suffisamment basses, ce qui n’est pas le cas. Certains matériaux ont la tendance opposée : les spins des
atomes magnétiques voisins préfèrent être antiparallèles. On parle alors d’antiferromagnétisme. Dans ce
cas les deux atomes magnétiques sont souvent séparés spatialement par un atome non-magnétique. Ces
considérations mènent à une interaction effective entre deux spins S1 et S2 , appelée interaction d’échange:
où J est une constante de l’ordre de l’énergie électrostatique entre les deux électrons correspondants. Dans
le cas ferromagnétique, J est positif, alors qu’il est négatif dans le cas antiferromagnétique.
Hystérèse
La relation existant entre B et H dans un matériau ferromagnétique est non-linéaire. Non seule-
ment est-elle non linéaire, mais elle dépend aussi de l’histoire du matériau, c’est-à-dire des champs
appliqués aux instants précédents. On dit qu’il y a hystérèse. Une relation typique est indiquée à la
figure 2.1. La raison de cette irréversibilité dans le comportement de l’aimantation est la suivante :
les moments magnétiques microscopiques ne sont pas globalement alignés dans le matériau, mais
plutôt regroupés en domaines mésoscopiques. Dans chacun des domaines il y a alignement des
dipôles atomiques, mais l’orientation de M varie d’un domaine à l’autre. La séparation spontanée
du milieu en différents domaines survient dans le but de minimiser l’énergie magnétique totale de
l’échantillon, c’est-à-dire l’énergie d’un domaine dans le champ créé par les autres domaines. Au
fond, il s’agit ici de l’interaction dipolaire qui prend sa revanche sur l’interaction d’échange. Cette
dernière est beaucoup plus forte, mais a une portée beaucoup moindre que la première. Il est donc
avantageux d’avoir un alignement des moments à courte distance (à l’intérieur des domaines)
mais une orientation presque antiparallèle aux longues distances. Lorsqu’un champ externe est
appliqué, certains domaines, ayant une aimantation parallèle au champ appliqué, sont avantagés
par rapport aux autres et croissent à leur dépend. Les murs séparant les différents domaines sont
donc mobiles. Cependant, le mouvement des murs dissipe de l’énergie et il s’agit là de l’origine
2.4. Magnétisme macroscopique 65
B
Br
H
Hc
de l’hystérèse. Dans les matériaux dits durs, les murs se déplacent difficilement, de sorte que
l’aimantation conserve longtemps sa valeur. Ces matériaux ont une valeur élevée de Hc (voir la
figure) et sont utilisés comme aimants permanents. Au contraire, dans les matériaux mous, les
murs se déplacent facilement et l’aimantation s’aligne sur le champ appliqué sans trop dissiper
d’énergie, ce qui est utile dans les transformateurs, etc.
4π
(H2 − H1 ) ∧ n = − K (2.90)
c
En effet, soit n1 le vecteur normal à une boucle rectangulaire infinitésimale traversant l’interface
et n2 = n1 ∧ n. Alors n2 · (H2 − H1 ) = (4π/c)n1 · K ou
4π
(n1 ∧ n) · (H2 − H1 ) = n1 · (n ∧ (H2 − H1 )) = n ·K (2.91)
c 1
d’où on peut conclure la formule ci-haut, en raison de l’arbitraire de n1 .
Solution : Utilisons les coordonnées cylindriques avec origine au centre de la barre, laquelle s’étend de
z = − 12 l à z = 12 l. La normale à la barre est r̂. L’aimantation est M = Mẑ. M étant uniforme à l’intérieur,
∇∧M = 0 et la seule trace que laisse cette aimantation est une densité superficielle
Commençons par calculer la contribution à l’induction magnétique à l’origine d’une mince boucle de courant
située à z: la loi de Biot-Savart donne
adϕ ϕ̂ ∧ r0
Z
I
dB = −
c r3
(2.93)
2πI a2
= ẑ
c r3
Z l/2 dz
Bz (0) = 2πa2 M
−l/2 (a2 + z 2 )3/2
l/2 (2.94)
= 4πM p
a2 + (l/2)2
= 4πM sin θmax
où θmax est l’angle que fait le bord extrême du barreau avec le plan z = 0. Dans le cas du barreau très
mince, on trouve B = 4πM au centre et donc H(0) = 0. À l’extrémité, il faut intégrer différemment :
Z l dz
Bz (l/2) = 2πa2 M
0 (a2 + z 2 )3/2
l (2.95)
= 2πM √
a2 + l2
0
= 2πM sin θmax
Dans la limite du barreau étroit, on trouve Bz (l/2) = 2πM et donc H(l/2) = −2πM: le champ est opposé
à l’induction magnétique. On voit que dans ce cas la relation entre B et H n’est pas linéaire. En dehors du
barreau, H = B puisque M = 0.
Dans la limite du barreau étroit et près de sa surface, on trouve que B ∼ 0 en dehors et B ∼ 4πM à
l’intérieur, en utilisant la loi d’Ampère sous sa forme intégrale.
Très loin du barreau, on doit retrouver le champ d’un dipôle magnétique m = πa2 lM.
H = −∇Φm (2.96)
2.4. Magnétisme macroscopique 67
Il y a alors une parfaite analogie entre les situations électrostatiques sans charge libre (ρ = 0,
ρs = 0) et les situations magnétostatiques sans courant libre (J = 0, K = 0). Cette analogie se
manifeste avec la correspondance suivante :
En effet, les lois différentielles et les conditions aux limites sont les mêmes dans ce cas :
∇·D = 0 ∇·B = 0
E = −∇Φ H = −∇Φm
⇔
Ek continu Hk continu
D⊥ continu B⊥ continu
Une différence essentielle demeure, cependant : la relation entre B et H peut très souvent être non
linéaire, ce qui est plus rare entre E et D. D’autre part, les situations électrostatiques comportent
souvent des corps conducteurs sur lesquels le potentiel électrique Φ est constant, ce qui ne se
produit pas pour le potentiel scalaire magnétique.
Si la relation entre B et H est linéaire (B = µH) alors la relation ∇·B = 0 implique
∇2 Φm = 0: le potentiel magnétique satisfait à l’équation de Laplace. La méthode de séparation
des variables peut alors être utilisée, comme dans le cas électrostatique.
Dans le cas plus général où la relation entre B et H n’est pas linéaire, comme par exemple en
présence d’aimants permanents, on ne peut qu’affirmer ceci :
On définit alors une densité de ‘charge magnétique’ ξ = −∇·M et une densité superficielle de
charge magnétique ξs = n · M à une interface, comme pour la charge liée à un diélectrique. La
solution à l’équation de Poisson nous est heureusement connue :
Z 0 Z
3 0 ξ(r ) M
Φm (r) = d r 0
+ da · (2.99)
|r − r | |r − r0 |
où dΩ0 = sin θ0 dθ0 dϕ0 , (θ0 , ϕ0 ) étant les coordonnées d’un point de source sur la sphère. Si γ est l’angle entre
r et r0 , on peut écrire
Z X rl
Φm (r) = a2 M dΩ0 cos θ0 <
P (cos γ)
l+1 l
(2.101)
r
l=0 >
où r< = r, r> = a si r < a et vice-versa si r > a. On utilise ensuite la formule d’addition des harmoniques
sphériques :
l
4π X ∗ 0 0
Pl (cos γ) = Ylm (θ , ϕ )Ylm (θ, ϕ) (2.102)
2l + 1
m=−l
Substituons cette expression dans Φm , en remarquant que cos θ0 = CY10 (θ0 , ϕ0 ), avec C =
p
4π/3:
4π X rl Z
Φm (r) = a2 M C <
l+1
dΩ0 Y10 (θ0 , ϕ0 )Ylm
∗
(θ0 , ϕ0 )Ylm (θ, ϕ) (2.103)
2l + 1 r
l,m >
4π cos θ
Φm (r) = M a3 2 (r > a) (2.106)
3 r
Ce qui implique que H (= B) a la forme dipolaire (2.43). On remarque que le moment dipolaire associé
est alors égal à M fois le volume de la sphère. À l’intérieur de la sphère, on a plutôt Φm = 43 πM z, ce qui
implique
4 8
H = − πM B = πM (2.107)
3 3
Le champ H s’oppose donc à l’aimantation, contrairement à B .
Problèmes
Problème 2.1
Une boucle circulaire de rayon a porte un courant I. Calculez le champ magnétique partour le long de l’axe
de la boucle.
Problème 2.2
Un solénoı̈de cylindrique de rayon a et de longueur L comporte n tours par unité de longueur et porte un
courant I. On suppose que L a, c.-à-d. que le solénoı̈de est très étroit. Calculez B et A à l’intérieur et à
l’extérieur du solénoı̈de, en restant toutefois proche de celui-ci (de sorte que le solénoı̈de apparait infini).
2.4. Magnétisme macroscopique 69
Problème 2.3
Exprimez la contribution quadrupolaire A(4) d’une densité de courant arbitraire au vecteur potentiel en
fonction de J.
Problème 2.4
Une sphère de rayon a porte une charge Q uniformément distribuée sur sa surface et est en rotation sur
elle-même avec une vitesse angulaire ω. Calculez le moment dipolaire magnétique associé.
Problème 2.5
Une coquille sphérique de rayon intérieur a et de rayon extérieur b est faite d’un matériau magnétique
linéaire de perméabilité µ. Cette coquille est placée dans un champ d’induction magnétique externe
uniforme B0 = B0 ẑ. Montrez que l’induction magnétique à l’intérieur de la coquille est
9µ
Bint = B0 ẑ (2.108)
(2µ + 1)(µ + 2) − 2(a/b)3 (µ − 1)2
Dans ce problème, il est indiqué d’utiliser le potentiel scalaire magnétique et la solution à l’équation de
Laplace en coordonnées sphériques avec symétrie azimutale. Les conditions aux limites aux deux interfaces
(à r = a et r = b) sont importantes. Un tel objet permet d’écranter un champ magnétique externe : il y a
un facteur de réduction R défini par Bint. = RB0 . Discutez de la dépendance de ce facteur en fonction de
µ quand µ 1 et en fonction de l’épaisseur de la coquille.
Problème 2.6
Considérez une disque aimanté de rayon a et d’épaisseur h, tel que h a (disque mince). Ce disque porte
une aimantation uniforme M = Mẑ (on prend ẑ le long de l’axe de symétrie du disque). Calculez B et
H au centre du disque, le long de l’axe de symétrie [considérez par exemple l’origine (0, 0, 0) au milieu
du disque et un point sur la surface : (0, 0, h/2)]. Faites le calcul de deux façons différentes : (a) à l’aide
de la loi de Biot-Savart et du courant lié. (b) à l’aide de la ‘charge magnétique’ et du potentiel scalaire
magnétique.
Problème 2.7
Un matériau supraconducteur a la propriété d’exclure tout champ magnétique de son intérieur. Ainsi, si
on plonge un corps supraconducteur dans un champ magnétique uniforme, une densité superficielle de
courant K sera induite à la surface du supraconducteur pour annuler le champ magnétique à l’intérieur du
supraconducteur. Considérons une sphère supraconductrice de rayon a, plongée dans un champ magnétique
uniforme B0 = B0 ẑ. L’objet de ce problème est de calculer le champ magnétique partout à l’extérieur de la
sphère et de calculer la densité superficielle de courant induite.
a) Sachant que la seule source de courant se trouve sur la surface de la sphère, écrivez le champ B
à l’extérieur comme le gradient d’une fonction (B = −∇Φm ) et trouvez B partout à l’extérieur en
solutionnant l’équation de Laplace pour ΦM , avec les conditions aux limites appropriées, à savoir que la
composante normale de B doit s’annuler à la surface de la sphère. Quel est le dipôle magnétique induit de
la sphère? Pourquoi peut-on utiliser l’équation de Laplace dans ce problème?
b) Trouvez la densité superficielle de courant K0 sur la sphère, à partir de la solution trouvée plus haut.
c) Décrivez les ressemblances et les différences entre ce problème et celui d’une sphère conductrice plongée
dans un champ électrique uniforme en ce qui a trait aux conditions aux limites et aux lignes de champ.
d) Trouvez une expression possible pour le potentiel vecteur A à l’extérieur et à l’intérieur de la sphère.
CHAPITRE 3
Électrodynamique
3.1 Induction
3.1.1 Loi de Faraday
On définit la force électromotrice (f.é.m.) V d’un circuit fermé C comme étant la circulation
du champ électrique autour de ce circuit :
I
V = dl · E (3.1)
C
En électrostatique et magnétostatique la force électromotrice le long d’un circuit fermé est nulle
puisque le champ électrique est alors irrotationnel.* Dans les situations où les courants dépendent
du temps le champ électrique n’est plus irrotationnel et une f.é.m. non nulle est possible, comme
l’ont montré les travaux de M. Faraday à partir des années 1830. La loi de Faraday stipule que la
force électromotrice d’un circuit fermé est proportionnelle à la dérivée par rapport au temps du
flux magnétique passant à travers la boucle du circuit.
Le flux magnétique à travers une surface S est évidemment défini comme
Z
ΦB [S] = da · B (3.2)
S
Comme ∇·B = 0, ce flux ne dépend pas de la forme précise de la surface mais uniquement du
contour C qui lui sert de frontière. En effet, en utilisant le théorème de Stokes pour B, le flux
magnétique peut être exprimé comme la circulation du potentiel vecteur le long du circuit :
Z
ΦB [C] = dl · A (3.3)
C
* On parle souvent de la force électromotrice d’une source de courant, par exemple d’une pile.
Il s’agit alors du travail par unité de charge accompli par le champ électrique le long d’un circuit
allant du pôle positif au pôle négatif de la pile. Un tel circuit n’est pas fermé; pour le fermer il faut
parcourir l’intérieur de la pile, du pôle négatif au positif. Le travail total est alors nul, dans une
situation stationnaire.
3.1. Induction 71
En utilisant le théorème de Stokes pour E, on peut écrire cette loi sous forme différentielle :
1 ∂B
∇∧E = − (3.5)
c ∂t
Nous supposons ici que le circuit C est immobile, c’est-à-dire qu’il ne se déplace pas dans le
temps. Dans le cas contraire une force électromotrice supplémentaire peut venir du fait que les
charges du fil conducteur subissent la force de Lorentz lorsque le fil se déplace.
Remarquons que le champ électrique induit par la variation du champ magnétique tend à
produire un courant dont l’effet est de s’opposer, par le champ magnétique qu’il produit, à la
variation du flux magnétique. C’est la loi de Lenz. Remarquons aussi que la force de Lorentz
agissant sur une particule chargée est toujours perpendiculaire au champ magnétique; ce dernier
ne peut donc produire un travail sur un ensemble de charges en mouvement. Cependant, une
variation de ce champ magnétique dans le temps induit un champ électrique qui, lui, peut fournir
un travail mécanique.
3.1.2 Inductance
Considérons une boucle C dans laquelle circule un courant I. Le champ magnétique produit
par cette boucle dépend linéairement de I, ainsi bien sûr que le flux magnétique ΦB traversant la
boucle et dû à son courant. On peut donc écrire la relation
ΦB = cLI (3.6)
où L est l’inductance de la boucle, une quantité qui ne dépend que de la forme de la boucle (bien
entendu, c est la vitesse de la lumière). Il s’ensuit que la f.é.m. produite par une variation du
courant I de la boucle est
dI
V = −L (3.7)
dt
(i)
Considérons ensuite un système de N boucles Ci , portant des courants Ii . Le flux ΦB
traversant la boucle Ci est une fonction linéaire des courants circulant dans toutes les boucles :
(i)
X
ΦB = c Mij Ij (3.8)
j
Les coefficients Mij sont les coefficients d’inductance mutuelle. Évidemment, Mii = Li . Donnons
une expression générale pour ces coefficients, à partir de la loi de Biot et Savart sous la forme
(2.8). Supposons que ri est la coordonnée d’un point sur la boucle Ci . Le potentiel vecteur A(ri )
dû au courant circulant dans Cj est alors
Ij dlj
I
A(ri ) = (3.9)
c Cj |ri − rj |
où dlj est l’élément de longueur de la boucle Cj au point rj . D’après l’expression (3.3), le
coefficient d’inductance mutuelle est alors
dli · dlj
I I
1
Mij = 2 (3.10)
c Ci j Cj |ri − rj |
72 3. Électrodynamique
Cette expression est également valide dans le cas i = j. On voit que l’inductance a les unités
d’une accélération inverse.
Enfin, nous pouvons facilement démontrer que le travail nécessaire à la mise en place d’un
courant I dans un circuit d’inductance L est égal à W = 21 LI 2 . En effet, la puissance (travail par
˙ ou encore P = 1 Ld(I 2 )/dt.
unité de temps) exercée par la f.é.m. induite est égal à P = V I = LI I, 2
En intégrant par rapport au temps de I = 0 jusqu’au courant final on obtient W = 12 LI 2 . Ce
travail a permis d’emmagasiner de l’énergie dans le champ magnétique produit par le circuit.
2N I
Z b dr
Z h 2N I b
dz = h ln (3.11)
c a r 0 c a
Le flux magnétique traversant le circuit est N fois ce résultat, en raison des N enroulements. L’inductance
est donc
2N 2 b
L = 2 h ln (3.12)
c a
Comme illustration, considérons un circuit comprenant une capacité qui se charge. Il n’y
a aucun courant de charge entre les plaques de la capacité, même si la circulation de H est
3.2. Équations de Maxwell 73
visiblement non nulle. C’est le courant de déplacement (lié à l’augmentation de E entre les
plaques) qui contribue, par la loi d’Ampère, à la circulation de H.
Dans un conducteur soumis à un champ sinusoı̈dal (fréquence ω) le courant de déplacement
vaut environ ω/σ fois J. C’est une fraction appréciable du courant de charge aux fréquences
élevées (ω ∼ σ).
Remarques
1. Il s’agit d’un ensemble de 8 équations du premier ordre, couplées, pour les 12 quantités E, D,
B et H. Pour qu’une solution soit possible, il faut également fournir des équations d’état:
4. dans le vide (ρtot. = 0 et Jtot. = 0), il y a symétrie de dualité par rapport à l’échange
E→B B → −E (3.20)
Toute solution des équations de Maxwell dans le vide devra refléter cette symétrie.
74 3. Électrodynamique
1 ∂ξ
Φ→Φ− A → A + ∇ξ (3.23)
c ∂t
Ceci se vérifie aisément puisque ∇∧∇ = 0 et que les deux dérivées de ξ s’annulent mutuellement
dans E.
La non unicité des potentiels nous permet d’imposer des conditions sur ces derniers sans que
cela ne restreigne les champs. On désigne ces conditions sous le nom de jauges:
1. Jauge de Lorentz: on impose la condition
1 ∂Φ
∇·A + =0 (3.24)
c ∂t
Si les potentiels ne satisfont pas à cette condition, alors il est toujours possible de trouver une
transformation de jauge (c.-à-d. une fonction ξ) telle que les nouveaux potentiels A0 et Φ0 y
satisfassent. En effet, on a
1 ∂Φ0 1 ∂Φ 1 ∂2ξ
∇·A0 + = ∇·A + + ∇2 ξ − 2 2 (3.25)
c ∂t c ∂t c ∂t
Il suffit alors de résoudre l’équation d’onde avec source pour ξ pour trouver une solution à
la condition de Lorentz. La condition de Lorentz ne fixe pas les potentiels de façon unique
puisqu’une transformation de jauge telle que
1 ∂2ξ
∇2 ξ − =0 (3.26)
c2 ∂t2
n’affecte pas la condition de Lorentz.
3.3. Potentiels électromagnétiques 75
∇·A = 0 (3.27)
Pour satisfaire cette condition à l’aide d’une transformation de jauge et à partir d’un potentiel
A quelconque, il suffit de résoudre l’équation de Poisson :
∇·A(r0 )
Z
0 1
2
∇·A = ∇·A + ∇ ξ ⇒ ξ(r) = d3 r0 (3.28)
4π |r − r0 |
Dans le cas statique, ces deux jauges sont bien sûr équivalentes.
En mécanique quantique, cet hamiltonien peut être utilisé tel quel, sauf que p devient
l’opérateur différentiel −ih̄∇. L’équation de Schrödinger pour la fonction d’onde ψ(r, t) est alors
∂ψ h̄2
ih̄ =− [∇ − (iq/h̄c)A(r, t)]2 ψ + qΦ(r, t)ψ (3.30)
∂t 2m
Une question se pose immédiatement : étant donné que ce sont les potentiels qui apparaissent
dans cet hamiltonien et non les champs E et B, est-ce que l’équation de Schrödinger est affectée
par les transformations de jauge? La réponse à cette question est qu’on peut étendre la définition
d’une transformation de jauge en incluant une multiplication de la fonction d’onde ψ par une
phase eiqξ/h̄c et que l’équation de Schrödinger est alors invariante lors d’une transformation de
jauge :
1 ∂ξ
Φ→Φ−
c ∂t
A → A + ∇ξ (3.31)
ψ → eiqξ/h̄c ψ
Pour démontrer ce fait plus aisément, il est commode de définir les opérateurs différentiels
suivants, qu’on appelle dérivées covariantes:
∂ iq
Dt ≡ + Φ
∂t h̄ (3.32)
iq
D ≡∇− A
h̄c
En fonction de ces opérateurs différentiels, l’équation de Schrödinger prend la forme suivante :
h̄2 2
ih̄Dt ψ = − D ψ (3.33)
2m
76 3. Électrodynamique
La propriété principale des dérivées covariantes est la covariance sous transformation de jauge,
c’est-à-dire que, lors d’une transformation de jauge, les dérivées covariantes sont modifiées par un
simple facteur de phase local :
Il est alors manifeste que l’équation de Schrödinger est invariante lors des transformations de
jauge. Le fait crucial ici est que la phase qui multiplie ψ dépend de r et de t. On sait qu’en
mécanique quantique la fonction d’onde peut être multipliée par une phase globale sans en
modifier le sens physique. On apprend ici que la fonction d’onde d’une particule chargée peut
être multipliée par une phase locale sans non plus modifier son sens physique, pourvu que
cette multiplication s’accompagne d’un changement des potentiels électromagnétiques , tel qu’en
l’Éq. (3.31). Le fait de pouvoir transformer une invariance globale en invariance local grâce à
la présence de potentiels scalaire et vecteur est à la base des théories de jauge qui décrivent les
interactions forte (QCD), faible et électromagnétique .
1 ∂
∇2 Φ + ∇·A = −4πρ (3.35)
c ∂t
1 ∂ 1 ∂2A 4π
∇(∇·A) − ∇2 A + ∇Φ + 2 2 = J (3.36)
c ∂t c ∂t c
On s’est servi de l’identité
∇∧(∇∧A) = ∇(∇·A) − ∇2 A (3.37)
1 ∂2Φ
∇2 Φ − = −4πρ
c2 ∂t2 (3.38)
1 ∂2A 4π
∇2 A − 2 2 = − J
c ∂t c
∇2 Φ = −4πρ (3.39)
ρ(r0 , t)
Z
Φ(r, t) = d3 r0 (3.40)
|r − r0 |
(en l’absence de sources, on peut donc prendre Φ = 0). L’équation (3.36) devient alors
1 ∂2A 4π 1 ∂Φ
∇2 A − =− J+ ∇ (3.41)
c2 ∂t2 c c ∂t
3.4. Énergie du champ électromagnétique 77
∇0 · J
Z
4πJl = −∇ d3 r0 (3.42)
|r − r0 |
∇0 · J
Z
∂Φ
∇ = −∇ d3 r0 (3.43)
∂t |r − r0 |
4π
∇2 Φ = −4πρ et ∇2 A = − J (3.45)
c
(le courant est toujours transverse dans ce cas) Les solutions sont bien sûr celles que l’on connaı̂t.
1
E= (E · D + B · H) (3.46)
8π
alors que le flux d’énergie (le courant associé à la densité E) est donné par le vecteur de Poynting:
c
S= E∧H (3.47)
4π
Le vecteur S représente la quantité d’énergie passant par unité de temps et de surface à travers
un élément de surface perpendiculaire à S. Le théorème de Poynting tel que formulé ici n’est
valable que pour un matériau linéaire.
Pour démontrer ce théorème, considérons premièrement le travail effectué sur les charges et
les courants d’un système fermé. La dérivée par rapport au temps de l’énergie cinétique des
particules matérielles chargées (c’est-à-dire la puissance) est
Z
dK
= d3 r J · E (3.48)
dt V
Cette relation se démontre en considérant la forme explicite de la densité de courant pour des
particules ponctuelles. On trouve
Z X X
d3 r J · E = qi vi · E = vi · Fi (3.49)
V i i
78 3. Électrodynamique
où Fi est la force de Lorentz agissant sur la particule i. La force magnétique est orthogonale à la
vitesse et donc n’effectue aucun travail sur les charges en mouvement.
Substituons pour J son expression tirée de la loi d’Ampère :
Z
dK c 3 1 ∂D
= dr ∇∧H − ·E (3.50)
dt 4π V c ∂t
1 ∂B
∇·(E ∧ H) = H · ∇∧E − E · ∇∧H = − H · − E · ∇∧H (3.51)
c ∂t
où nous avons utilisé la loi de Faraday dans la dernière équation. En substituant pour E · ∇∧H,
on trouve
Z Z
dK 1 3 ∂B ∂D
=− d r H· +E· − d3 r ∇·S (3.52)
dt 4π V ∂t ∂t V
Le dernier terme peut être transformé en intégrale de surface. Supposons maintenant que le
milieu considéré est linéaire, de sorte que H est proportionnel à B en tout temps, comme D est
proportionnel à E. Le premier terme peut alors s’écrire comme une dérivée. On trouve finalement
Z I
∂ 1 3
K+ d r (H · B + E · D) = − da · S (3.53)
∂t 8π V S
∂E
+ ∇·S = −J · E (3.54)
∂t
∂E
+ ∇·S = −σE2 (3.55)
∂t
Le terme de droite, avec son signe négatif, représente les pertes d’énergie causées par la résistance
du conducteur.
Le vecteur de Poynting n’est pas seulement égal au flux d’énergie du champ : il est aussi
proportionnel à la densité d’impulsion associée au champ électromagnétique . On montre que
cette densité d’impulsion, dans le vide, est précisément donnée par S/c2 .
3.5. Ondes Électromagnétiques 79
1 ∂2E
∇2 E − =0 (3.56)
c2 ∂t2
À cette fin, écrivons les équations de Maxwell dans le vide, en l’absence de sources :
∇·E = 0 ∇·B = 0
1 ∂B 1 ∂E (3.57)
∇∧E + =0 ∇∧B − =0
c ∂t c ∂t
Prenons le rotationnel de la loi de Faraday :
1 ∂ 1 ∂2E
∇∧(∇∧E) + ∇∧B = ∇(∇·E) − ∇2 E + 2 2 (3.58)
c ∂t c ∂t
Dans le dernier terme nous avons utilisé la loi d’Ampère. Puisque ∇·E = 0, on obtient l’équation
d’onde (3.56) pour E. En prenant plutôt le rotationnel de la loi d’Ampère et en utilisant la loi
de Faraday dans le dernier terme, on trouverait la même équation d’onde, mais en fonction du
champ magnétique B.
Une solution de l’équation d’onde a été proposée par d’Alembert dans le cas des ondes planes,
c’est-à-dire quand les champs ne varient que dans une direction, disons l’axe des x. L’équation
d’onde devient alors (pour E)
∂2E 1 ∂2E
− =0 (3.59)
∂x2 c2 ∂t2
Définissons les variables u1 = x − ct et u2 = x + ct dont les opérateurs de dérivées sont
∂ 1 ∂ 1 ∂ ∂ 1 ∂ 1 ∂
= − = + (3.60)
∂u1 2 ∂x c ∂t ∂u2 2 ∂x c ∂t
80 3. Électrodynamique
E = E1 (u1 ) + E2 (u2 )
(3.62)
= E1 (x − ct) + E2 (x + ct)
soit la somme de deux ondes se propageant à la vitesse c, mais dans des directions opposées :
positive pour E1 et négative pour E2 . Les fonctions E1 et E2 ne sont pas contraintes davantage par
l’équation d’onde. Autrement dit, l’onde progressive a une forme quelconque, pas nécessairement
sinusoı̈dale. Cette onde conservera sa forme au cours de la propagation (pas de dispersion), ce qui
est caractéristique d’une solution de l’équation d’onde.
Comme l’équation d’onde est linéaire, la solution générale de l’équation d’onde en trois
dimensions est une superposition d’ondes planes quelconques se propageant dans des directions
arbitraires. Il s’agit bien sûr d’une solution décrite par un nombre infini de paramètres. Notons
cependant que l’équation d’onde ne suffit pas à déterminer complètement le champ E, puisqu’il y
a plus d’information dans l’ensemble des équations de Maxwell que dans la seule équation d’onde.
Une onde plane se déplaçant dans une direction quelconque s’écrit plutôt comme ceci :
Le vecteur k est le vecteur d’onde, sa grandeur k = |k| est le nombre d’onde et son inverse est
proportionnel à la longueur d’onde: λ = 2π/k.
En appliquant la loi de Gauss ∇·E = 0 à l’onde (3.65), on trouve k · E = 0, c’est-à-dire
que le champ électrique est perpendiculaire au vecteur d’onde (c’est-à-dire à la direction de
propagation). D’autre part, la loi de Faraday nous permet de retrouver le champ magnétique
associé à cette onde :
1 ∂B ω
∇∧E + = 0 =⇒ ik ∧ E − i B = 0 (3.66)
c ∂t c
(le champ magnétique associé oscille bien sûr à la même fréquence ω). On trouve donc
k
B = k̂ ∧ E k̂ = (3.67)
|k|
3.5. Ondes Électromagnétiques 81
En conclusion, dans une onde électromagnétique plane se propageant dans le vide, les vecteurs E,
B et k sont mutuellement perpendiculaires et les champs E et B ont la même amplitude et la
même phase. Les ondes électromagnétiques sont donc transversales et non longitudinales, car les
quantités oscillantes sont perpendiculaires à la direction de propagation. Les ondes électromagné-
tiques de vecteur d’onde k existent donc en deux polarisations possibles.
Signalons que toute onde électromagnétique peut être considérée comme une superposition
linéaire d’ondes planes de vecteurs d’ondes différents. La solution générale à l’équation d’onde
dans le vide est donc de la forme suivante :
d3 k
Z
E(r, t) = Re E (k) ei(k·r−ωt) (3.68)
(2π)3 0
où E0 (k) est une amplitude qui dépend du vecteur d’onde k.
Notons que tout ceci peut changer lorsqu’on considère des ondes électromagnétiques se
propageant dans des milieux matériels. Les amplitudes des champs E et B ne sont alors plus
les mêmes et leurs phases peuvent être différentes si le milieu absorbe progressivement l’onde.
D’autre part, la vitesse de phase ω/k dépend alors de la fréquence et est généralement plus petite
que c. Ceci mène au phénomène de dispersion, par lequel une onde non sinusoı̈dale (un paquet
d’onde) change de forme au cours de sa propagation.
1 ∂2A 4π
∇2 A − 2 2
=− J (3.76)
c ∂t c
Résoudre cette équation requiert une analyse un peu détaillée que nous ne pouvons faire ici.
Indiquons simplement la solution :
J(r0 , t − |r − r0 |/c)
Z
1
A(r, t) = d3 r0 (3.77)
c |r − r0 |
La caractéristique principale de cette solution est qu’elle ressemble fort à celle qu’on aurait si la
densité de courant J était indépendante du temps [Éq. (2.8)], sauf que la densité de courant qui
figure dans l’intégrale n’est pas prise au temps d’observation t, mais au temps retardé t0 :
|r − r0 |
t0 = t − (3.78)
c
En fait, t0 précède t par le temps que prendrait un signal lumineux pour se rendre du point de
source r0 au point d’observation r. Le potentiel (3.77) porte par conséquent le nom de potentiel
retardé. Mise à part cette différence, d’ailleurs prévisible et très plausible, le potentiel (3.77) a la
même allure que son vis-à-vis statique (2.8).
Supposons maintenant que la densité de courant oscille dans le temps à une fréquence ω:
et l’intégration est prise sur la région qui contient les charges en mouvement.
Il y a ici deux longueurs caractéristiques : la dimension d du système radiant (longueur de
l’antenne, par exemple) et la longueur d’onde λ = 2π/k du rayonnement. Nous allons supposer
que le point d’observation est très éloigné: r d et r λ. On dit alors que l’observateur se
situe dans la zone de rayonnement. Prenons l’origine des coordonnées au centre de la source des
courants. Si n est le vecteur unité dans la direction du point d’observation et r la distance à
l’origine, on peut faire l’approximation suivante suffisamment loin des sources :
|r − r0 | ≈ r − n · r0 (3.81)
3.5. Ondes Électromagnétiques 83
On écrit donc
eikr
Z
0
A(r) = d3 r0 J(r0 ) e−ikn·r (3.82)
cr
Nous allons aussi faire l’approximation que la dimension d du système radiant est beaucoup
plus petite que la longueur d’onde du rayonnement, de sorte que kn · r0 est toujours très petit
dans la région où J est appréciable. On peut alors remplacer l’exponentielle par l’unité. Cette
approximation porte le nom d’approximation dipolaire, pour une raison qui deviendra bientôt
claire. On écrit donc
eikr
Z
Adip. (r) = d3 r0 J(r0 ) (3.83)
cr
Dans la section 2.3.1 nous avons montré que cette intégrale s’annule dans une situation statique.
Dans une situation dynamique, ce n’est plus le cas. En effet,
Z Z Z
3 3 ∂
d r Jk (r) = − d r ∂i Ji xk = d3 r ρ(r)xk = −iωdk (3.84)
∂t
d = d0 e−iωt (3.85)
eikr
Adip. (r) = −ik d (k = ω/c) (3.86)
r
Cependant,
eikr
∇(eikr /r) = (ikr − 1) n (3.88)
r2
Dans la zone de rayonnement, le premier terme de la parenthèse domine et on peut négliger le
deuxième. Donc
eikr
B = k2 n∧d (3.89)
r
Le champ électrique s’obtient facilement du champ magnétique en observant que, dans la zone de
rayonnement, on a localement affaire à une onde plane de vecteur d’onde kn et par conséquent
eikr
E = B ∧ n = k2 (d − (d · n)n) (3.90)
r
Observons que le champ électrique est bien sûr perpendiculaire à n et orienté dans le plan
commun à n et d.
L’énergie transportée par cette onde peut être calculée à l’aide du vecteur de Poynting
moyenné dans le temps :
c
hSi = |E |2 n (3.91)
8π 0
84 3. Électrodynamique
Figure 3.1. Patron de rayonnement dipolaire : coupe à angle ϕ fixe et vue tridimen-
sionnelle
ck 4 2
Z
c 4 2
P = k d0 dΩ sin2 θ = d (3.93)
8π 3 0
de la lumière incidente et est proportionnelle à ω 4 . Par conséquent, la lumière bleue est plus
efficacement diffusée que la lumière rouge (environ 5 fois plus) et le ciel apparaı̂t bleu, puisque
la lumière nous provenant du ciel est essentiellement la lumière solaire diffusée par l’atmosphère.
En contrepartie, plus la distance parcourue dans l’atmosphère par la lumière qui nous parvient
directement du Soleil est grande, plus sa composante bleue est faible (elle a été diffusée) et plus le
Soleil apparaı̂t rougeâtre. Ceci explique la couleur du Soleil couchant et son apparence rougeâtre
encore plus prononcée lorsque l’atmosphère est humide ou polluée (les particules en suspension et
les gouttelettes d’eau sont d’excellent diffuseurs).
(il est sous-entendu que nous prenons la partie réelle des quantités complexes). Cette particule
ponctuelle constitue à elle seule une distribution de charge, avec moments dipolaire, quadrupo-
laire, etc. En particulier, le moment dipolaire de cette distribution oscille dans le temps à une
fréquence ω:
d = d0 e−iωt d0 = er0 (3.95)
On peut donc supposer que cette particule émet un rayonnement dipolaire. En fait, la forme
précise de son rayonnement ne sera pas exactement donnée par la formule (3.92): celle-ci ne sera
correcte que dans la limite où |r0 |ω c. D’après la formule (3.93), la puissance totale rayonnée
par cette particule est
ck 4 2 2e2 a2
P = d0 = (3.96)
3 3c3
√
où a = ω 2 r0 / 2 est l’accélération quadratique moyenne de la particule dans son mouvement
d’oscillation. Il se trouve que cette formule est beaucoup plus générale que ce qui est impliqué ici.
Elle donne la puissance rayonnée par une particule chargée d’accélération instantanée a et porte
le nom de formule de Larmor :
2e2 a2
P = (3.97)
3c3
Sous cette forme, cependant, elle n’est valable que pour une particule non relativiste.
Appendice 3.A
Théorème de Helmholtz
Un champ vectoriel F fini et continu qui s’annule à l’infini peut être exprimé comme la somme
d’un gradient et d’un rotationnel :
F = ∇φ + ∇∧A (3.98)
Pour prouver ce théorème, il faut obtenir une expression pour φ et A en fonction de F. Pour
ce faire, définissons le champ W suivant :
F(r0 )
Z
1
W(r) = d3 r0 (3.99)
4π |r − r0 |
86 3. Électrodynamique
∇2 W = −F (3.100)
Alors
F = −∇2 W = −∇(∇·W) + ∇∧(∇∧W) = ∇φ + ∇∧A (3.102)
∇0 · F(r0 )
Z
1
φ(r) = − d3 r0 (3.103)
4π |r − r0 |
∇0 ∧ F(r0 )
Z
1
A(r) = d3 r0 (3.104)
4π |r − r0 |
Le champ F peut donc être décomposé en une somme de deux termes : F = Fl + Ft , où
Fl = ∇φ est la partie longitudinale (irrotationnelle) et Ft = ∇∧A la partie transverse (sans
divergence). Si un champ F est sans divergence, on conclut donc que seul Ft est non nul, alors
que si F est irrotationnel, seul Fl est non nul.
Problèmes
Problème 3.1
Quelles raisons ont motivé les corporations distribuant de l’énergie électriques à adopter le système de
courant alternatif, par opposition à celui de courant continu défendu par Edison? Pourquoi le courant est-il
triphasé?
Problème 3.2
Considérez une boucle C dans laquelle circule un courant I, qu’on déplace de manière rigide avec une
vitesse v dans un champ magnétique B(r). En considérant la force de Lorentz agissant sur la charge en
mouvement dans chaque élément de circuit, démontrez que la f.é.m. causée par ce mouvement est donnée
par I
1
V = dl · (v ∧ B) (3.105)
c C
Cette force électromotrice est également présente dans le référentiel qui se déplace avec la boucle. À quoi
serait alors due cette f.é.m., alors que la vitesse de la boucle est nulle dans ce référentiel?
Problème 3.3
Une coquille sphérique mince de rayon a fait d’un matériau conducteur est en rotation autour de l’axe
z avec une vitesse angulaire ω. Cette coquille est placée dans un champ magnétique uniforme B = B0 ẑ.
Calculez la force électromotrice induite entre le pôle nord et l’équateur de la coquille.
3.A. Théorème de Helmholtz 87
Problème 3.4
Un transformateur primitif est fait d’une pièce de fer doux (un matériau ferromagnétique de très grande
perméabilité et de faible hystérèse) ayant la forme d’un tore ou d’une barre repliée sur elle-même. Deux
circuits (un primaire et un secondaire) s’enroulent autour de ce tore, avec respectivement N1 et N2 tours.
Expliquez comment transformer les tensions avec un tel dispositif et expliquez le rôle du fer doux.
Problème 3.5
Soit un long câble coaxial formé de deux cyclindres conducteurs concentriques de rayons a et b respective-
ment (a < b). Il s’agit d’un circuit dont l’une des composantes entoure l’autre complètement (on suppose
que les deux cylindres sont reliés entre eux à l’extrémité du câble).
a) En supposant que l’épaisseur de ces cylindres est négligeable, calculez l’inductance par unité de longueur
de ce câble.
b) Répétez le calcul, mais cette fois en supposant que le cylindre intérieur est plein et que le courant s’y
distribue de façon uniforme.
c) Avec les données de (a), calculez l’énergie magnétique de deux façons : par la formule U = 21 LI 2 et par
l’intégration directe de la densité d’énergie magnétique, et montrez que les résultats coincident.
CHAPITRE 4
Relativité Restreinte
compatible avec le principe de relativité, i.e., telle que la vitesse de propagation de la lumière soit
la même dans les deux référentiels.
1. Cette transformation doit être linéaire, sinon un corps au repos dans S ne serait pas en
mouvement rectiligne uniforme dans S 0 . Les fonctions fi et g ci-haut sont donc linéaires en xi
et t.
2. Les coordonnées perpendiculaires à la vitesse relative v peuvent être dénotées r⊥ . La seule
transformation acceptable de ces coordonnées est r0⊥ = K(v)r⊥ , où K(v) est un facteur
numérique qui dépend de v, tel que K(0) = 1. Par symétrie, K(v) ne doit pas dépendre de la
direction de v: K(−v) = K(v). Puisque la transformation inverse est
On doit avoir K(v)2 = 1 ou K(v) = 1. Donc les coordonnées transverses à la vitesse relative
ne sont pas affectées par la transformation : r0⊥ = r⊥ .
3. On cherche donc une transformation linéaire (x, t) → (x0 , t0 ) qu’on peut écrire sous la forme
suivante :
x0 = ax + bt
(4.3)
t0 = ex + f t
Les quatre constantes a, b, e, f ne dépendent que de v, la vitesse relative.
La vitesse u = dx/dt d’un objet dans S se transforme alors de la manière suivante :
où a(v) est une fonction de v encore à déterminer et A(v) est la matrice de transformation. On
sait que a(0) = 1. De plus, a(v) ne devrait pas dépendre du signe de v, mais seulement de son
carré, en raison de l’invariance par rapport à l’inversion de l’espace : les changements simultanés
x → −x x0 → −x0 v → −v (4.7)
90 4. Relativité Restreinte
ne devraient pas affecter la transformation. Or il se trouve que cette inversion change a(v) en
a(−v) et c’est tout. On en conclut que a(−v) = a(v). Finalement, on sait que la transformation
inverse (i.e. de S 0 à S) s’obtient en renversant le signe de v: A(v)A(−v) = 1 ou encore
1 −v 1 v 1 0
a(v)a(−v) = (4.8)
−v/c2 1 v/c2 1 0 1
De cette contrainte on trouve que
1
a(v) = p ≡γ (4.9)
1 − v 2 /c2
On peut donc écrire la transformation finale comme
x0 = γ(x − vt) y0 = y
(4.10)
t0 = γ(t − vx/c2 ) z0 = z
Elle porte le nom de transformation de Lorentz. Il est manifeste que deux référentiels inertiels ne
peuvent se déplacer l’un par rapport à l’autre à une vitesse supérieure à c.
Mentionnons que la loi inverse de transformation des vitesses s’écrit
u−v u0 + v
u0 = ou u = (4.11)
1 − uv/c2 1 + u0 v/c2
Il est alors manifeste que même si u0 et v sont arbitrairement proche de c, la vitesse u observée
dans S ne peut pas dépasser c. La vitesse de la lumière est donc une vitesse limite.
4.1.3 Intervalle
On appelle événement un point (ct, r) de l’espace-temps. On introduit un facteur de c devant
le temps t pour que tous les membres du multiplet aient les mêmes unités; en fait, ct peut servir
de coordonnée temporelle aussi bien que t. Étant donné deux événements (ct1 , r1 ) et (ct2 , r2 ) on
définit l’intervalle entre ces événements comme
s2 ≡ c2 t2 − r2 (r ≡ r1 − r2 t ≡ t1 − t2 ) (4.12)
L’intervalle est un invariant, c’est-à-dire une quantité qui reste la même dans tous les référentiels
inertiels. Pour s’en convaincre il suffit de choisir l’origine de deux référentiels S et S 0 au point r2 ,
au temps t2 et de substituer les transformations de Lorentz dans la définition de s02 :
s02 = (ct0 )2 − x02 − y 02 − z 02
= γ 2 (c2 t2 + (v/c)2 x2 − 2xtv) − γ 2 (x2 + v 2 t2 − 2xtv) − y 2 − z 2 (4.13)
= (ct)2 − x2 − y 2 − z 2
L’intervalle est l’analogue dans l’espace-temps de la distance (au carré) en l’espace de dimen-
sion 3: l’intervalle est invariant par rapport aux changements de référentiel inertiel tout comme
la distance est invariante par rapport aux rotations des axes x̂, ŷ et ẑ. La différence importante
est que l’intervalle peut être positif ou négatif alors que la distance carrée est toujours positive ou
nulle. Deux événements séparés par un intervalle nul sont dits sur le même cône de lumière, car
un signal lumineux émis à r1 au temps t1 atteindra r2 exactement à t2 (on suppose que t1 < t2 ).
Deux événements séparés par un intervalle négatif ne peuvent pas être reliés de manière causale,
i.e. par un signal se propageant à une vitesse égale ou inférieure à c. Par contre, deux événements
séparés par un intervalle positif peuvent être reliés par des signaux se propageant plus lentement
que c.
4.2. Formalisme covariant 91
D’autre part, imaginons une horloge qui produit un tic avec une période T dans son référentiel
propre, toujours à la même position x. Désignons par t1 et t2 les deux instants auxquels se
produisent deux tics consécutifs. Comme ct0 = γ(ct − βx), on a, pour une même position x, la
relation t02 − t01 = γ(t2 − t1 ), ou encore T 0 = γT :
T
T0 = p (4.15)
1 − v 2 /c2
Ceci signifie que la période de l’horloge en mouvement telle que mesurée par l’observateur est plus
grande que sa période au repos.
L’élément de volume d3 r dans un référentiel en mouvement est contracté comme les longueurs
le sont, car une seule des trois différentielles perpendiculaires qui le forment est contractée, dans
la direction du mouvement. On écrit donc d3 r0 = d3 r/γ. L’élément de temps dt, lui, est dilaté:
dt0 = γdt. Cependant, l’élément d’espace-temps est invariant : d3 rdt = d3 r0 dt0 , puisque les deux
effets se compensent.
où Λµ ν est une matrice de transformation. L’indice de gauche numérote les rangées, l’indice de
droite les colonnes et les indices répétés sont sommés. Pour la transformation (4.10) cette matrice
est
γ −βγ 0 0
−βγ γ 0 0
Λ= (4.17)
0 0 1 0
0 0 0 1
En notation matricielle, la transformation de la coordonnée s’écrit simplement x0 = Λx, où x est le
vecteur-colonne abritant les composantes xµ .
Examinons maintenant comment se transforme l’opérateur gradient augmenté de la dérivée
par rapport au temps, noté
∂ 1 ∂
∂µ = = ,∇ (4.18)
∂xµ c ∂t
Exprimé en fonction des nouvelles coordonnées obtenues par transformation de Lorentz, ceci
devient
∂x0ν ∂
∂µ = = Λν µ ∂ν0 (4.19)
∂xµ ∂x0ν
En fonction de Λ−1 , la matrice inverse de Λ, on écrit
∂µ0 = (Λ−1 )ν µ ∂ν ou d0 = (Λ−1 )t d (4.20)
où d est le vecteur-colonne abritant les composantes ∂µ . Bien entendu, Λ−1 est la matrice de la
transformation inverse obtenue simplement en changeant le signe de la vitesse :
γ βγ 0 0
βγ γ 0 0
Λ−1 = (4.21)
0 0 1 0
0 0 0 1
1
On écrit parfois 4-vecteur pour abréger quadri-vecteur.
4.2. Formalisme covariant 93
Contraction
La contraction d’un indice covariant avec un indice contravariant se définit comme une opération
par laquelle les deux indices sont identifiés et sommés. Spécifiquement, la contraction d’un vecteur
covariant Aµ avec un vecteur contravariant B µ est le produit Aµ B µ , dans lequel µ est sommé de
0 à 4. La quantité résultante est un invariant, comme on peut le constater d’après les lois de
transformation ci-haut :
A0µ B 0µ = (Λ−1 )α µ Λµ β Aα B β
= δβα Aα B β (4.24)
α
= Aα B
Par exemple, l’opérateur différentiel xµ ∂µ est un invariant, c’est-à-dire que sa forme est la même
dans tous les référentiels inertiels.
En notation matricielle, l’opération de contraction ci-haut s’écrit
Au contraire, la contraction de deux indices covariants (ou deux indices contravariants) n’est
pas invariante. Autrement dit, l’expression Aµ Aµ (on somme sur µ) n’est pas numériquement la
même dans tous les référentiels. L’apparition d’une telle expression dans un calcul signifie qu’on
y effectue des opérations qui ne sont pas valables dans tous les référentiels et est probablement
l’effet d’une erreur!
Tenseurs
Ce sont les propriétés de transformation qui définissent la notion de vecteur. On généralise cette
notion à des quantité comportant plus d’un indice, qu’on appelle généralement tenseurs. Un
quadri-tenseur de rang n est une quantité comportant 4n composantes spécifiées par n indices
grecs allant de 0 à 3. Spécifiquement, un tenseur contravariant de rang 2 est un ensemble de
quantités B µν se transformant comme suit lors d’un changement de référentiel :
B 0µν = Λµ α Λν β B αβ (4.26)
Un tenseur deux fois covariant se noterait plutôt Aµν et se transformerait comme suit :
C 0µ ν = Λµ α (Λ−1 )β ν C α β (4.29)
1 0 0 0
0 −1 0 0
g= (4.31)
0 0 −1 0
0 0 0 −1
Par définition, une transformation de Lorentz est une transformation linéaire spécifiée par une
matrice Λ telle que l’intervalle s2 est invariant. Autrement dit,
ce qui signifie, en comparant avec (4.28), que gµν se comporte comme un tenseur covariant de rang
deux. On l’appelle le tenseur métrique car il définit l’intervalle en fonction de xµ et l’intervalle est
la mesure naturelle des ‘distances’ dans l’espace-temps. Le tenseur métrique a ceci de particulier
qu’il a exactement la même forme (les mêmes composantes) dans tous les référentiels, même s’il
se transforme comme un tenseur.
En prenant l’inverse de la relation matricielle ci-haut, on trouve
Ce qui signifie, en comparant avec (4.27), que la matrice g−1 définit un tenseur contravariant de
rang deux. On note ses composantes g µν . D’après la forme explicite de g, on voit que g−1 = g, ou
g µν = gµν .
Le tenseur métrique permet de transformer un quadri-vecteur contravariant en quadri-vecteur
covariant, par contraction des indices. Ainsi, la quantité gµα B α est un quadri-vecteur covariant,
4.2. Formalisme covariant 95
(gA)0 = gA0
= g(Λ−1 )t A
(4.38)
= Λ(Λ−1 )g(Λ−1 )t A
= ΛgA
Notons que la contraction gµα g αν est précisément le delta de Kronecker δµν . En notation ma-
tricielle, ceci signifie simplement que g2 = 1. Dans cette notation matricielle, le tenseur covariant
Bµν obtenu de B µν par abaissement des indices est gBg, alors que le tenseur mixte B µ ν est Bg et
que l’autre tenseur mixte Bµ ν est gB.
4.2.3 Exemples
Temps propre
Le temps propre τ d’un objet en mouvement est simplement le temps tel qu’il s’écoule dans le
référentiel de l’objet. Dans le référentiel inertiel S 0 qui se déplace avec l’objet à un instant donné,
la différentielle de temps propre coı̈ncide avec dt0 . Cependant, l’objet peut en général être accéléré
de sorte que le référentiel S 0 n’est pas le même à tous les instants. Quel que soit le référentiel dans
lequel on observe l’objet, la différentielle de temps propre peut toujours être exprimée comme
suit :
1q µ p dt
dτ ≡ dx dxµ = dt 1 − u2 /c2 = (4.40)
c γ
où u est la vitesse de l’objet, qui peut dépendre du temps. En effet, dans le référentiel S 0 ,
dxµ dxµ = c2 dτ 2 et de plus l’expression ci-haut est invariante, puisque dxµ dxµ est un invariant.
Pour l’obtenir le temps propre total écoulé sur la trajectoire de l’objet, on intègre dτ :
Z Z Z
dt p
τ= dτ = = dt 1 − u2 /c2 (4.41)
γ
96 4. Relativité Restreinte
quadri-vitesse
Considérons un objet se déplaçant à une vitesse instantanée u. On définit sa quadri-vitesse uµ
comme
dxµ
uµ = (4.42)
dτ
La différentielle dτ est invariante, alors que dxµ se transforme comme les coordonnées; donc uµ
est une quadri-vecteur. En fonction de la vitesse u, les composantes explicites sont
uµ : (cγ, γu)
Par contraction on obtient l’invariant uµ uµ = c2 . Le temps propre écoulé le long d’une trajectoire
peut alors s’écrire comme Z Z
1
τ = dτ = 2 uµ dxµ (4.43)
c
Quadri-accélération
On définit aussi la quadri-accélération
duµ duµ
aµ ≡ =γ (4.44)
dτ dt
duµ
uµ aµ = uµ
dτ
1 d (4.46)
= (u uµ )
2 dτ µ
=0
quadri-courant
Un autre quadri-vecteur est formé par la densité de courant et la densité de charge :
J µ : (cρ, J) (4.48)
∂ρ
∇·J + = ∂µ J µ = 0 (4.49)
∂t
4.3. Formulation relativiste de l’électrodynamique 97
Comme la conservation de la charge doit être valide dans tous les référentiels et que ∂µ est
un quadri-vecteur, le quadri-courant J µ en est forcément un lui aussi. Une autre façon de se
convaincre que J µ est un quadri-vecteur est de considérer l’expression de la densité de courant
associée à un ensemble de charges ponctuelles en mouvement :
X
J(r) = q (i) v(i) δ(r − r(i) ) (4.50)
i
Ici q (i) , r(i) et v(i) sont respectivement la charge, la position et la vitesse de la ie particule. La
fonction delta dans la formule ci-haut n’est pas un invariant de Lorentz, puisqu’elle est définie par
la relation Z
d3 r(i) δ(r − r(i) ) = 1 (4.51)
Comme q (i) et la fonction delta modifiée sont des invariants, on conclut que J µ est bel et bien un
quadri-vecteur.
1 ∂2Φ
∇2 Φ − = −4πρ
1 ∂Φ
c2 ∂t2 ∇·A + =0 (4.55)
1 ∂2A 4π c ∂t
∇2 A − 2 2 = − J
c ∂t c
98 4. Relativité Restreinte
où ξ est une fonction scalaire (i.e. invariante). Évidemment, le passage de Aµ vers Aµ se fait à
l’aide du tenseur métrique et nous utiliserons à la fois les versions covariante et contravariante des
quadri-vecteurs ou quadri-tenseurs dans ce qui suit.
Lors d’un changement de référentiel, les composantes temporelle et spatiales de J µ (ou de
µ
A ) se combinent. Par exemple, considérons un fil rectiligne selon l’axe des x dans lequel circule
un courant I et qui ne porte pas de charge nette. On sait qu’une telle configuration de charge
donne lieu à un champ magnétique qui circule autour du fil. Si on procède à un changement
de référentiel S → S 0 le long de cet axe, la densité de courant dans S 0 est différente de J et
une densité de charge ρ0 = −vγJ/c2 apparait. Cette densité de charge fait apparaı̂tre un champ
électrique radial dans S 0 . Il s’agit d’un fait capital : non seulement les phénomènes magnétiques
et électriques sont-ils unifiés par les équations de Maxwell, mais la distinction entre les deux n’est
pas invariante : elle est fonction du référentiel.
0 Ex Ey Ez
−Ex 0 −Bz By
Fµν : (4.61)
−Ey Bz 0 −Bx
−Ez −By Bx 0
4.3. Formulation relativiste de l’électrodynamique 99
Ce tenseur renferme donc toutes les composantes des champs magnétique et électrique. Les
équations de Maxwell s’expriment de façon particulièrement concise en fonction de Fµν . Pour le
démontrer il suffit de constater que
∂ ν Fµν = ∂ ν (∂µ Aν − ∂ν Aµ )
4π (4.63)
= −∂ ν ∂ν Aµ = − Jµ
c
où nous avons utilisé l’Éq. (4.56). Si on développe cette équation composante par composante on
retrouve les lois de Gauss et d’Ampère-Maxwell.
Quant aux deux autres équations de Maxwell, on les retrouve de la façon suivante : définissons
le tenseur de Levi-Civita en 4 dimensions εµνρλ et considérons l’expression
Comme εµνρλ est antisymétrique dans tous ses indices et que ∂ν ∂ρ et ∂ν ∂λ sont symétriques, cette
expression s’annule. En développant cette expression composante par composante, on constate
qu’il s’agit de la loi de Faraday et de l’équation ∇·B = 0. Notons en passant qu’on définit le
tenseur électromagnétique dual
1
F̃ µν = εµνρλ Fρλ (4.65)
2
L’expression explicite de ce tenseur dual est la même que pour Fµν sauf pour le fait que les rôles
de E et de B sont échangés. Plus précisément, F̃ µν s’obtient à partir de F µν par les substitutions
E → B et B → −E:
0 −Bx −By −Bz
B 0 Ez −Ey
F̃ µν : x (4.66)
By −Ez 0 Ex
Bz Ey −Ex 0
On peut donc écrire les équations de Maxwell comme
4π µ
∂ν F νµ = J ∂ν F̃ µν = 0 (4.67)
c
On obtient explicitement
Ex0 = Ex Bx0 = Bx
Ey0 = γ(Ey − βBz ) By0 = γ(By + βEz ) (4.69)
Ez0 = γ(Ez + βBy ) Bz0 = γ(Bz − βEy )
Exprimons ces transformations différemment en appelant respectivement E⊥ et Ek les com-
posantes de E perpendiculaire et parallèle à la direction de v, la vitesse relative des référentiels.
On écrit alors
v
E0⊥ = γ E⊥ + ∧ B⊥ E0k = Ek
c (4.70)
v
B0k = Bk
B0⊥ = γ B⊥ − ∧ E⊥
c
Remarques
1 1
Tr (FF) = − F µν Fµν = E2 − B2 (4.72)
2 2
Cette quantité est donc indépendante du référentiel. Ceci était a priori évident puisque F µν Fµν
est une contraction complète des indices.2
Un autre invariant s’obtient en considérant plutôt la contraction F̃ µν Fµν . On calcule que
1 µν
F̃ Fµν = E · B (4.73)
4
Donc le produit scalaire E · B est invariant. Cependant, il faut remarquer que cette quantité
change de signe lors d’une transformation de parité, car E est un vecteur polaire et B un vecteur
axial. Une autre façon d’obtenir cet invariant est de considérer l’expression det F, en notant que
det Λ = det Λt = 1.
2
Le même calcul exercé sur Tr F donne des résultats triviaux puisque Tr F = 0 dans tous les référentiels.
4.4. Applications 101
Remarques
1. Si E et B sont de même grandeur dans un référentiel, ils le sont dans tous les référentiels
(E2 − B2 = 0).
2. Si B = 0 et E 6= 0 dans un référentiel, alors E ne pourra être plus petit dans un autre
référentiel.
3. Si E et B sont perpendiculaires dans un référentiel, ils le seront dans tous les référentiels
(E · B = 0).
4.4 Applications
4.4.1 Champ d’une charge ponctuelle en mouvement uniforme
Nous nous intéressons ici au champ électrique produit par une charge ponctuelle q se
déplaçant à une vitesse constante. Nous allons calculer ce champ au moyen d’une transforma-
tion de Lorentz effectuée à partir du référentiel S dans lequel la charge est au repos. Dans ce
référentiel, le champ électrique qu’elle produit est
r q
E = q 3 ou Ei = 3 xi (4.74)
r r
Appliquons les transformations (4.70) pour calculer les champs produits par une charge ponctuelle
q, cette fois se déplaçant avec une vitesse v = vẑ. Nous nous plaçons donc dans un référentiel S 0
se déplaçant à une vitesse −v par rapport à S. On trouve
Ex0 = γEx Ey0 = γEy Ez0 = Ez (4.75)
0
Nous nous intéressons au champ lorsque la charge est à l’origine dans S , c’est-à-dire au temps
t0 = 0. La relation entre les coordonnées est alors
x = x0 y = y0 z = γz 0 (4.76)
Donc la distance à l’origine dans S est
p
r= x02 + y 02 + γ 2 z 02 (4.77)
On peut alors écrire
x0i
Ei0 = qγ (4.78)
(x02 + y 02 + γ 2 z 02 )3/2
Soit θ l’angle que fait le rayon vecteur avec v. Alors z 0 = r0 cos θ et x02 + y 02 = r02 sin2 θ et
x02 + y 02 + γ 2 z 02 = r02 sin2 θ + γ 2 (1 − sin2 θ)
(4.79)
= γ 2 r02 (1 − β 2 sin2 θ)
On peut donc écrire
(1 − β 2 )r0
E0 = q (4.80)
r03 (1 − β 2 sin2 θ)3/2
Remarques
Remarques
où uµ est la quadri-vitesse de la particule. La constante de proportionalité −mc2 à été ajoutée par
convention (m est la masse de la particule). Puisque dτ = dt/γ, on peut écrire
Z p
2
S0 = −mc dt 1 − v 2 /c2 (4.89)
L’hamiltonien est dans ce cas égal à l’énergie cinétique de la particule. Dans l’approximation non
relativiste, H devient
1
H0 ≈ mc2 + mv 2 (4.94)
2
On reconnaı̂t l’énergie cinétique usuelle, plus une constante : l’énergie au repos mc2 .
Le premier terme du membre de gauche est par définition la force F (la dérivée de la quantité de
mouvement par rapport au temps). On peut donc écrire
e e
F = −e∇Φ − Ȧ + vi ∇Ai (4.102)
c c
Notons que la dérivée temporelle Ȧ est une dérivée totale par rapport au temps et non une
dérivée partielle. Cette dérivée est non seulement due à la variation temporelle de A à une
position donnée, mais aussi au déplacement de la particule dans le temps d’un point à un autre :
∂A
Ȧ = + (v · ∇)A (4.103)
∂t
mais
−vj ∂j Ai + vj ∂i Aj = εijk εkrs vj ∂r As = (v ∧ (∇∧A))i (4.104)
Donc on obtient finalement
e
F = eE + v ∧ B (4.105)
c
Ce qui est effectivement la force de Lorentz.
où L est la densité lagrangienne et d4 x = d3 rdt. Le lagrangien est alors l’intégrale de la densité
lagrangienne sur l’espace : Z
L= d3 r L (4.107)
De plus, cette action doit consister de deux parties : une décrivant la dynamique propre du champ
et une décrivant l’interaction du champ avec la matière chargée. Cette dernière partie nous est
déjà connue : il s’agit de S1 ci haut. En exprimant cette dernière en fonction du quadri-vecteur de
densité de courant au lieu de la quadri-vitesse, on a
Z
1
S1 = − d4 x J µ Aµ (4.108)
c
Il reste à trouver S0 , l’action propre du champ ou la densité lagrangienne correspondante L0
Pour que S0 soit invariant, il faut que L0 aussi soit invariant; mais attention! le lagrangien,
lui, n’est pas invariant. La raison en est que l’élément de volume d’espace-temps d4 x est invariant
par rapport aux transformations de Lorentz, alors que d3 r ne l’est pas. En effet,
0
4 0 ∂x
d x = det d4 x = det Λ d4 x (4.109)
∂x
106 4. Relativité Restreinte
où Λ est la matrice de la transformée de Lorentz. On vérifie sans peine que det Λ = 1. Donc
d4 x0 = d4 x ou d3 r0 dt0 = d3 rdt. Mais dt0 6= dt, et donc d3 r0 =
6 d3 r.
Par analogie avec la mécanique d’une particule, on veut que L0 dépende de Aµ et de ses
premières dérivées. Comme les équations de Maxwell sont linéaires, il fait que L0 soit au plus
quadratique dans les champs. La seule possibilité invariante est alors
1 µν 1
L0 = − F Fµν = − F µν ∂µ Aν (4.110)
16π 8π
Le facteur de −1/16π a été ajouté dans le but de retrouver les équations de Maxwell dans leur
normalisation habituelle par rapport au courant. La densité lagrangienne totale est donc
1 1 µν
L = − J µ Aµ − F ∂µ Aν (4.111)
c 8π
Utilisons maintenant le principe de la moindre action pour obtenir les équations du mouvement.
Il faut pour cela ajouter une variation δAµ au quadri-vecteur du potentiel et s’assurer que la
variation correspondante de l’action est nulle au premier ordre. Notons que
Si cette variation s’annule pour n’importe quelle variation δAν , il faut que le tenseur électroma-
gnétique satisfasse aux équations suivantes :
4π ν
∂µ F µν = J (4.115)
c
Cette équation constitue une partie des équations de Maxwell. L’autre partie provient de la
définition même de Fµν en fonction de Aµ et peut être exprimée comme ∂µ F̃ νµ = 0 (voir ci-haut).
Tr T = T µ µ = 0 (4.117)
puisque Tr 1 = 4.
4.5. Formulation Lagrangienne 107
Un calcul explicite donne les expressions suivantes pour les différentes composantes de T:
1 1
T ij = − Ei Ej + Bi Bj − δij (E2 + B2 )
4π 2
1 (4.118)
T 0k = S
c k
1
T 00 = (E2 + B2 ) = E
8π
La partie spatiale de T, en fait la matrice −Tij , porte le nom de tenseur des contraintes de
Maxwell.
Le tenseur d’énergie-impulsion satisfait à la condition suivante :
1 λν
∂µ T µν = F Jλ (4.119)
c
Dans la deuxième équation, nous avons utilisé les équations de Maxwell et nous avons anti-
symétrisé le second facteur du deuxième terme. Les trois derniers termes peuvent aussi s’écrire
1 να ρλ
g F (∂ρ Fλα + ∂λ Fαρ + ∂α Fρλ ) (4.121)
2
après avoir changé les indices de sommation et transposé au besoin. On reconnaı̂t entre par-
enthèses le deuxième groupe d’équations de Maxwell. Ce terme s’annule donc et il reste l’équation
(4.119).
L’interprétation physique de la relation (4.119) est la conservation de l’énergie et de
l’impulsion. Pour le démontrer, considérons premièrement la composante ν = 0 ce cette rela-
tion :
1
∂µ T µ0 = F λ0 Jλ
c (4.122)
1
=− E·J
c
Intégrons cette dernière relation dans un volume V , en utilisant l’expression (4.50) pour le courant
associé à un ensemble de charges ponctuelles :
Z
1 X (i) (i)
d3 r ∂µ T µ0 = − q v · E(r(i) )
V c i
1 dW (4.123)
=−
c dt
Z Z
1 ∂ 3 00 ∂ k0
= drT + d3 r T
c ∂t V V ∂xk
108 4. Relativité Restreinte
où W est l’énergie cinétique des sources du champ contenue dans le volume V . Dans la dernière
équation, le dernier terme est l’intégrale de la divergence d’un vecteur et peut se ramener à une
intégrale de surface. On écrit donc
Z I
∂ 3 00
W+ drT = −c dak T k0 (4.124)
∂t V S
Ceci signifie que le vecteur dont les composantes sont cT k0 est le flux de la quantité apparaissant
dans le membre de gauche, c’est-à-dire l’énergie. Ce vecteur n’est rien d’autre que S, le vecteur de
Poynting, et la densité d’énergie du champ électromagnétique est T 00 = (E2 + B2 )/8π.
Passons maintenant aux composantes spatiales (ν = i) de l’équation (4.119). Reprenant le
calcul précédent, on trouve
1
∂µ T µi = F λi Jλ
c
1
= −Ei ρ − εikl Bl Jk (4.125)
c
1
= − ρE + J ∧ B
c i
car F ik = εikl Bl . En utilisant l’expression du courant et de la densité de charge pour des particules
ponctuelles et en intégrant sur un volume V , le membre de droite devient la force de Lorentz totale
exercée sur les particules comprises dans ce volume :
Z
d3 r ∂µ T µi = −(Ftot. )i
V
Z Z (4.126)
1 ∂ 3 0i ∂ ki
= drT + d3 r T
c ∂t V V ∂xk
Le deuxième terme est l’intégrale d’une divergence et se transforme en intégrale de surface : cette
relation devient donc
Z I
∂ 1 3 0i
(ptot. )i + drT = − dak T ki (4.127)
∂t c V S
où p est l’impulsion totale des sources du champ dans le volume V . L’interprétation de
cette équation est que T 0i /c est la composante i de la densité de l’impulsion du champ
électromagnétique et que le vecteur dont les composantes sont T ki (k = 1, 2, 3) est le flux de
la composante i de l’impulsion du champ. La densité d’impulsion est donc reliée au vecteur de
Poynting :
1
densité d’impulsion : S (4.128)
c2
Dans le cas statique (aucune dépendance en temps), la force totale exercée sur les sources par le
champ est donc
I
(Ftot. )i = − dak T ki (4.129)
4.5. Formulation Lagrangienne 109
Problèmes
Problème 4.1
Un peu de science-fiction. Supposons qu’il soit possible de construire un vaisseau spatial capable de
transporter des êtres vivants d’un système planétaire à un autre. Dans le but d’offrir une gravitation
‘artificielle’ aux passagers, l’astronef évolue avec une accélération constante g dans son référentiel propre
(c.-à-d. le référentiel inertiel dont la vitesse coı̈ncide avec celle de l’astronef à un moment donné). On
suppose que cette accélération est maintenue jusqu’à mi-chemin, pour ensuite être renversée (changement
de signe) jusqu’à destination. Si d est la distance parcourue (référentiel terrestre), quel est le temps T que
prend le voyage, tel que mesuré par les passagers ? Si g = 10 m/s2 et d = 4, 3 années-lumière (distance de
l’étoile la plus proche) quel est le temps requis ? Quel serait T pour d = 2 × 106 années-lumière (distance
de la galaxie d’Andromède) ?
Problème 4.2
Un long fil rectiligne, orienté selon x̂, au repos dans le référentiel S, porte un courant continu I. Un
référentiel S 0 se déplace à une vitesse v = vx̂ par rapport à S. Quels sont les champs électrique et
magnétique dans S 0 (en fonction de la distance par rapport au fil) et à quelles densité de charge et quel
courant sont-ils attribuables?
Problème 4.3
Une charge ponctuelle e évolue dans des champs électrique (E = Eŷ) et magnétique (B = Bẑ) uniformes et
perpendiculaires, tels que E B. Donnez une expression pour la trajectoire de la particule en fonction du
temps (c’est-à-dire (x(t), y(t), z(t))) en supposant que r(0) = 0 et ṙ(0) = 0. Pour ce faire, transportez-vous
dans le référentiel où E = 0, trouvez la trajectoire dans ce référentiel et ensuite transformez le tout dans le
référentiel original. Supposez toujours que E B et montrez que la trajectoire est une cycloı̈de (tracez-la!).
Problème 4.4
Une onde plane a la forme suivante :
En partant du principe que la phase de l’onde est un invariant, montrez que la quantité (ω/c, k) est un
quadri-vecteur contravariant. Dans un référentiel S une onde électromagnétique plane se propage dans le
vide à une fréquence ω. Soit S 0 un référentiel se déplaçant à une vitesse v par rapport à S. La fréquence ω 0
de l’onde dans S 0 est différente (effet Doppler). Si v est perpendiculaire à k, montrez que
ω
ω0 = p (effet Doppler transverse)
1 − v 2 /c2
Quantité SI gaussien
longueur 1 mètre 102 centimètre
temps 1 seconde 1 seconde
masse 1 kilogramme 103 gramme
force 1 newton 105 dyne
7
énergie 1 joule 10 erg
puissance 1 watt 107 erg/s
charge électrique 1 coulomb 3.109 statcoulomb
potentiel électrique 1 volt 1/300 statvolt
champ électrique E 1 volt/m (1/3).104 statvolt/m
polarisation 1 coulomb/m2 3 × 4π.105 statcoulomb/cm2
résistance 1 ohm (1/3)2 .10−11 statohm
conductivité 1 mho/m (3)2 .109 1/sec
capacité 1 farad (3)2 .1011 statfarad
induction magnétique B 1 weber/m2 104 gauss
champ magnétique H 1 ampere/m (4π).10−3 oersted
aimantation 1 weber/m2 (1/4π).104 gauss
inductance 1 henry (1/3)2 .10−11 stathenry
Note : le chiffre 3 apparaissant dans la troisième colonne signifie en fait la valeur numérique de la
vitesse de la lumière (divisé par 108 ), à savoir 2,99792.
Annexe B
Variables Électromagnétiques, SI vs
gaussien
Quantité SI gaussien
D ε0 E + P E + 4πP
1
H B−M B − 4πM
µ0
∇·D ρ 4πρ
∇·B 0 0
∂B 1 ∂B
∇×E − −
∂t c ∂t
∂D 1 ∂D
∇×H J+ 4πJ +
∂t c ∂t
1
F q(E + v × B) q E+ v×B
c
c
S E×H E×H
4π
Annexe C
Formulaire
analyse vectorielle :
εijk εmnk = δim δjn − δin δjm
A × (B × C) = B(A · C) − C(A · B)
(A × B) · (C × D) = (A · C)(B · D) − (A · D)(B · C)
∇ × (∇ × A) = ∇(∇ · A) − ∇2 A
coordonnées cartésiennes :
∂ ∂ ∂
∇ = x̂
+ ŷ + ẑ
∂x ∂y ∂z
∂Ax ∂Ay ∂Az
∇·A= + +
∂x ∂y ∂z
∂A ∂Ay
∂Az y ∂Ax ∂Az ∂Ax
∇×A= − x̂ + − ŷ + − ẑ
∂y ∂z ∂z ∂x ∂x ∂y
coordonnées cylindriques :
∂ 1 ∂ ∂
∇ = r̂ + ϕ̂ + ẑ
∂r r ∂ϕ ∂z
1 ∂ 1 ∂Aϕ ∂Az
∇·A= (rAr ) + +
r ∂r r ∂ϕ ∂z
∂A
1 ∂Az ϕ ∂Ar ∂Az 1 ∂ ∂Ar
∇×A= − r̂ + − ϕ̂ + (rAϕ ) − ẑ
r ∂ϕ ∂z ∂z ∂r r ∂r ∂ϕ
1 ∂2Φ ∂2Φ
2 1 ∂ ∂Φ
∇ Φ= r + 2 +
r ∂r ∂r r ∂ϕ2 ∂z 2
coordonnées sphériques :
∂ 1 ∂ 1 ∂
∇ = r̂ + θ̂ + ϕ̂
∂r r ∂θ r sin θ ∂ϕ
1 ∂ 1 ∂ 1 ∂Aϕ
∇ · A = 2 (r2 Ar ) + (sin θAθ ) +
r ∂r r sin θ ∂θ r sin θ ∂ϕ
1 ∂ ∂Aθ 1 ∂Ar 1 ∂ 1 ∂ ∂Ar
∇×A= (sin θAϕ ) − r̂ + − (rAϕ ) θ̂ + (rAθ ) − ϕ̂
r sin θ ∂θ ∂ϕ r sin θ ∂ϕ r ∂r r ∂r ∂θ
∂2Φ
2 1 ∂ 2 ∂Φ 1 ∂ ∂Φ 1
∇ Φ= 2 r + 2 sin θ + 2 2
r ∂r ∂r r sin θ ∂θ ∂θ r sin θ ∂ϕ2
Bibliographie