Electromagnétisme

Table des matières

1 Les sources du champ électromagnétique : charges et courants électriques 6

1 Charges immobiles . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
1.1 Charges des particules élémentaires . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
1.2 Force entre deux charges ponctuelles . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
1.3 Distributions étendues de charges . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
2 Charges mobiles : courant électrique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
2.1 Intensité du courant . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
2.2 Densité de courant à 3D . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
2.3 Conducteurs de plus basse dimension . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
3 Loi de conservation de la charge . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
3.1 Cas d’un flot unidimensionnel . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
3.2 Cas du flot tridimensionnel . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
4 Conclusion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13

2 Le champ électrostatique 14
1 Champ électrique créé par une charge ponctuelle unique . . . . . . . . . . . . . . . . . . . . . . 14
1.1 Notion de champ électrique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
1.2 Lignes de champ électrique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
2 Champ électrique total créé par une distribution de charges : principe de superposition . . . . . 15
2.1 Distribution discrète de charges ponctuelles . . . . . . . . . . . . . . . . . . . . . . . . . 16
2.2 Distributions continues de charges . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
3 Mouvement d’une charge dans un champ électrique uniforme . . . . . . . . . . . . . . . . . . . 17
4 Champ électrique et symétries spatiales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18

3 Les propriétés du champ électrostatique 20

1 Circulation et potentiel électrostatique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
1.1 Circulation du champ électrique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
1.2 Potentiel électrique créé par une charge unique . . . . . . . . . . . . . . . . . . . . . . . 20
1.3 Lien avec les notions de travail et d’énergie potentielle . . . . . . . . . . . . . . . . . . . 21
2 Potentiel électrostatique pour les distributions de charges . . . . . . . . . . . . . . . . . . . . . 21
2.1 Distribution discrète de charges ponctuelles . . . . . . . . . . . . . . . . . . . . . . . . . 21
2.2 Distributions continues de charges . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
3 Flux du champ électrique et théorème de Gauss . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
3.1 Cas du champ créé par une seule charge ponctuelle source à travers une sphère centrée . 23
3.2 Cas du champ créé par une seule charge ponctuelle source à travers une surface quelconque 23
3.3 Cas général . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24
4 Exemple d’utilisation du théorème de Gauss : boule chargée . . . . . . . . . . . . . . . . . . . . 24
5 Equations de Maxwell de l’électrostatique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26
6 Exercices . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26

4 Le champ magnétique 28
1 Expériences historiques : Oersted and Ampère . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
1.1 Expérience d’Oersted . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
1.2 Expérience(s) d’Ampère (1820) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
2 Force de Lorentz comme définition du champ magnétique . . . . . . . . . . . . . . . . . . . . . 30
2.1 Force de Lorentz . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30

2
 2.2 Unités et ordres de grandeur . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31
2.3 Considérations de symétrie par rapport à un plan : cas général . . . . . . . . . . . . . . 31
3 Mouvement d’une particule chargée dans un champ magnétique uniforme . . . . . . . . . . . . 32
3.1 Mouvement cyclotron dans un champ magnétique uniforme . . . . . . . . . . . . . . . . 32
3.2 Considérations de symétrie : un exemple sur le mouvement cyclotron . . . . . . . . . . 34
4 Champ créé par des courants sources : loi de Biot et Savart (1820) . . . . . . . . . . . . . . . . 35
4.1 Enoncé pour des courants sources volumiques . . . . . . . . . . . . . . . . . . . . . . . . 35
4.2 Enoncé pour des circuits filiformes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35
4.3 Cas d’un fil rectiligne infini . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35
5 Action d’un champ magnétique sur un élément de circuit électrique . . . . . . . . . . . . . . . . 37
5.1 Force de Laplace . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 37
5.2 Interprétation de l’expérience d’Ampère sur les forces entre fils parallèles . . . . . . . . . 37
6 Exercices . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 37

5 Les propriétés générales du champ magnétique 39

1 Champ magnétique à flux conservatif . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 39
1.1 Cas particulier du champ créé par un fil rectiligne infini . . . . . . . . . . . . . . . . . . 39
1.2 Cas d’un circuit filiforme quelconque . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40
1.3 Démonstration par le calcul de la divergence . . . . . . . . . . . . . . . . . . . . . . . . . 40
2 Circulation du champ magnétique et théorème d’Ampère . . . . . . . . . . . . . . . . . . . . . . 40
3 Les symétries du champ magnétique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 41
4 Exemple d’application du théorème d’Ampère . . . . . . . . . . . . . . . . . . . . . . . . . . . . 42
5 Equations de Maxwell de la magnétostatique . . . . . . . . . . . . . . . . . . . . . . . . . . . . 43
6 Exercices . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 44

6 Le champ électromagnétique et ses sources dans le vide 45

1 Les équations de l’électrodynamique classique . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45
2 Découpler les équations de Maxwell . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 46
3 Les solutions des équations de Maxwell via les potentiels . . . . . . . . . . . . . . . . . . . . . . 47
4 Energie électromagnétique dans le vide . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 47
4.1 Densité de puissance cédée par le champ à la matière chargée . . . . . . . . . . . . . . . 48
4.2 Densité locale d’énergie du champ u(~r, t) . . . . . . . . . . . . . . . . . . . . . . . . . . . 48
4.3 Vecteur de Poynting R(~~ r, t) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 48
4.4 Loi de conservation de l’énergie . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 48
4.5 Expressions de u(~r, t) et du vecteur de Poynting en fonction des champs . . . . . . . . . 48

7 Ondes électromagnétiques dans le vide 50

1 Relation de dispersion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50
2 Structure des ondes et modes de Fourier . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 51
3 Etats de polarisation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53
4 Aspects énergétiques d’une OPPM . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53

8 Champ électrique dans les conducteurs 54

1 Relation constitutive courant-champ électrique . . . . . . . . . . . . . . . . . . . . . . . . . . . 54
2 Origine microscopique de la relation constitutive (Modèle de Drude) . . . . . . . . . . . . . . . 54
2.1 Passage à la loi d’Ohm macroscopique (1826) : conductance et résistance . . . . . . . . . 55
3 Conducteurs en équilibre électrostatique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 55
3.1 L’équilibre électrostatique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 55
3.2 Interface vide/métal . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56
3.3 Boule métallique dans un champ électrique statique . . . . . . . . . . . . . . . . . . . . . 56
4 Condensateurs . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57
4.1 Condensateur à l’équilibre électrostatique . . . . . . . . . . . . . . . . . . . . . . . . . . 57
4.2 Conducteur placé dans un champ extérieur . . . . . . . . . . . . . . . . . . . . . . . . . 58

3
9 Effet Hall : effet d’un champ magnétique statique et uniforme sur un conducteur 60
1 Modèle microscopique de la conduction sous champ magnétique . . . . . . . . . . . . . . . . . . 60
1.1 Modèle de Drude avec champ magnétique . . . . . . . . . . . . . . . . . . . . . . . . . . 60
1.2 Tenseurs de conductivité et de résistivité . . . . . . . . . . . . . . . . . . . . . . . . . . . 60
1.3 Effet Hall classique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 61
2 Origine de la force de Laplace sur un conducteur . . . . . . . . . . . . . . . . . . . . . . . . . . 61

10 Champ électrique dans la matière isolante 62

1 Polarisabilité d’un atome : apparition d’un dipôle induit . . . . . . . . . . . . . . . . . . . . . . 62
2 Vecteur polarisation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 63
3 Charges de polarisation volumiques et surfaciques . . . . . . . . . . . . . . . . . . . . . . . . . . 63
4 Equation de Maxwell pour le champ lissé . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 64
5 Notion de susceptibilité . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 64
5.1 Réponse linéaire et tenseur de susceptibilité électrique . . . . . . . . . . . . . . . . . . . 64
5.2 Milieux linéaires homogènes et isotropes (LHI) . . . . . . . . . . . . . . . . . . . . . . . 64
6 Matière condensée . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 65
6.1 Notion de champ local . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 65
6.2 Lien entre polarisabilité et susceptibilité . . . . . . . . . . . . . . . . . . . . . . . . . . . 65
7 Boule isolante dans un champ électrique statique (équation de Laplace) . . . . . . . . . . . . . . 66

11 Champ magnétique dans la matière 67

1 Dipôles magnétiques et matériaux magnétiques . . . . . . . . . . . . . . . . . . . . . . . . . . . 67
1.1 Matériaux paramagnétiques . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 67
1.2 Matériaux diamagnétiques . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 67
1.3 Matériaux ferromagnétiques . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 67
2 Vecteur aimantation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 67
3 Courants ampériens volumiques et surfaciques . . . . . . . . . . . . . . . . . . . . . . . . . . . . 68
4 Equations de Maxwell pour les champs lissés dans la matière . . . . . . . . . . . . . . . . . . . 68
5 Susceptibilité magnétique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 68
6 Boule aimantée dans un champ magnétique statique . . . . . . . . . . . . . . . . . . . . . . . . 68

12 Ondes dans les milieux matériels 69

1 Equations de Maxwell dans la matière . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 69
1.1 ~ et H
Utilisation des vecteurs D ~ et relations constitutives . . . . . . . . . . . . . . . . . 69
1.2 ~ B
Utilisation des vecteurs E, ~ et P~ et M~ . . . . . . . . . . . . . . . . . . . . . . . . . . . 69
2 Modèle microscopique de l’interaction lumière-matière : isolants . . . . . . . . . . . . . . . . . . 70
2.1 Modèle de l’électron élastiquement lié . . . . . . . . . . . . . . . . . . . . . . . . . . . . 70
3 Modèle microscopique de l’interaction lumière-matière : métaux . . . . . . . . . . . . . . . . . . 71
3.1 Modèle de Drude AC et conductivité optique . . . . . . . . . . . . . . . . . . . . . . . . 71
3.2 Lien avec le modèle de Lorentz (introduit pour les isolants) . . . . . . . . . . . . . . . . 73
4 Fonction de réponse . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 73
4.1 Origine de la dispersion en pulsation ω . . . . . . . . . . . . . . . . . . . . . . . . . . . . 73
4.2 Signification du caractère complexe de la fonction de réponse . . . . . . . . . . . . . . . 74
4.3 Analyse énergétique : parties dissipative et réactive de la fonction de réponse . . . . . . 74

13 Phénomènes aux interfaces : réflexion, transmission et confinement 76

1 Comportement des champs EM à une interface . . . . . . . . . . . . . . . . . . . . . . . . . . . 76
1.1 Analyse des équations de Maxwell "microscopiques" à divergence . . . . . . . . . . . . . 76
1.2 Analyse des équations de Maxwell "microscopiques" à rotationnel . . . . . . . . . . . . . 77
1.3 Analyse des équations de Maxwell pour les champs macroscopiques lissés . . . . . . . . . 77
2 Problème de diffusion (scattering) à une interface . . . . . . . . . . . . . . . . . . . . . . . . . . 77
2.1 Généralités . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 78
2.2 Interface entre deux diélectriques en incidence normale . . . . . . . . . . . . . . . . . . . 78
2.3 Interface vide-métal parfait . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 79

4
14 Electrodynamique des supraconducteurs 80
1 Introduction historique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 80
2 Relation constitutive des frères London . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 80
2.1 Modèle à deux fluides et fonction d’onde macroscopique . . . . . . . . . . . . . . . . . . 80
2.2 Relation phénoménologique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 80
3 Effet Meissner-Ochsenfeld . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 81

15 Annexe sur l’analyse vectorielle 82

1 Gradient . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 82
2 Divergence . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 82
3 Rotationnel . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 82

5
Chapitre 1

Les sources du champ électromagnétique :

charges et courants électriques

Dans ce chapitre, nous énonçons la loi de Coulomb qui donne la force à distance entre deux charges ponctuelles
immobiles. Puis nous passons à la description des charges mobiles qui sont responsables du courant électrique.
Le dernier paragraphe approfondit la notion de conservation locale de la charge.

1 Charges immobiles
Nous nous plaçons dans un référentiel donné, celui du laboratoire (qui est supposé Galiléen), et nous ne
considérons dans ce premier paragraphe que des charges immobiles dans ce référentiel.

1.1 Charges des particules élémentaires

La charge est une quantité que l’on associe à chaque particule élémentaire et qui caractérise son couplage
électromagnétique avec les autres particules. Une particule neutre (neutron, neutrino) a une charge nulle et
ne subit pas de force électromagnétique. L’électron a une charge qe = −e et le proton porte la charge opposée
qp = e. Dans la pratique, on utilise le terme "une charge" pour désigner une particule chargée.
L’unité de charge est le coulomb (symbole C). La charge élémentaire e vaut en Coulomb :

e = 1, 602 . 10−19 C . (1.1)

On voit que le Coulomb est une unité plutôt adaptée aux charges macroscopiques, formées d’un grand nombre
de particules élémentaires.
La charge totale d’un atome ou d’un ion est obtenue en additionnant les charges des électrons et des protons
qui le compose. Dans une réaction chimique ou même nucléaire la charge totale d’un système isolé est conservée.
Par exemple, la réaction nucléaire de désintégration du neutron :

n → p + e− + ν e , (1.2)

en un proton, un électron et un antineutrino, conserve la charge électrique. Ce concept fondamental de conservation

de la charge sera approfondi dans le dernier paragraphe de ce chapitre.
Dans notre entourage, on rencontre essentiellement la matière sous forme d’atomes ou de molécules neutres.
On verra que les charges ont tendance à s’associer par charges opposées afin de former des édifices neutres qui
sont plus stables. Dans les étoiles chaudes ces édifices neutres sont cassés (ionisés) et la matière se retrouve
sous forme d’un plasma formé essentiellement de protons et d’électrons. Sur Terre, les objets macroscopiques
sont principalement neutres mais on peut les électriser en arrachant des charges comme dans l’électrisation par
frottements (triboélectricité) : des électrons sont transférés d’un corps frotté vers l’autre, et on obtient deux
corps chargés. Le première machine électrostatique permettant de créer ainsi des charges a été réalisée par Von
Guericke (1602-1686) en Allemagne.

1.2 Force entre deux charges ponctuelles

Dans les années 1780, Coulomb a réalisé des expériences sur des billes électrisées (par frottements car à
l’époque la pile n’existait pas). En mettant au point une balance à torsion très précise capable de mesurer
de toutes petites forces, Coulomb a réussi à montrer que la force entre deux billes chargées varie comme

6
Figure 1.1 – Loi de Coulomb : Les forces d’interaction entre deux charges ponctuelles de même signe :
q1 q2 > 0. La distance entre les deux charges est notée r. Exercice : faire le dessin similaire pour le cas de deux
charges opposées.

l’inverse du carré de la distance les séparant. C’est une loi fondamentale qui est l’équivalent pour le couplage
électromagnétique de la loi de Newton sur l’attraction gravitationnelle des masses.

Loi de Coulomb
Considérons deux charges ponctuelles et immobiles placées dans le vide (seules dans l’Univers). La force
exercée par la charge 1 sur la charge 2 s’écrit comme (Fig. 1.1) :

1 q1 q2
F~12 = ~u12 , (1.3)
4π0 r2

où r = M1 M2 est la distance entre les deux charges q1 et q2 , et ~u12 est le vecteur unitaire pointant de la charge
1 vers la charge 2. Le facteur
1
= 9.109 (SI), (1.4)
4π0
dépend du système d’unités utilisé. Dans tout ce cours, on utilise le système international (SI) dans lequel les
forces sont en newtons, les charges en coulombs et les distances en mètres.
Remarque 1 : Dans le cadre de cette théorie, une charge seule dans l’univers ne subit pas de force électrique.
Autrement dit, une charge n’agit pas sur elle-même. D’ailleurs, si on essayait d’appliquer la formule de Coulomb
pour deux charges situées au même point (r = 0) on obtiendrait une force infinie.
Remarque 2 : Cette formule reste valable si on inverse les rôles entre les 2 charges. En effet la force exercée
par la charge 2 sur la charge 1 s’écrit comme :
1 q1 q2
F~21 = ~u21 = −F~12 , (1.5)
4π0 r2
ce qui est conforme au principe de l’action-réaction (Fig. 1.1). On a utilisé que ~u21 = −~u12 . Historiquement,
cette loi a été obtenue pour des billes macroscopiques. On verra plus loin dans le cours que tout se passe comme
si la charge était répartie au centre de la bille. De plus on sait maintenant que la force électromagnétique entre
deux particules élémentaires, par exemple entre un proton et un électron au sein d’un atome, est également
donnée par cette formule.

Gravitation vs électrostatique
Considérons deux masses ponctuelles et immobiles placées dans le vide. La force d’attraction gravitationnelle
exercée par la masse 1 sur la masse 2 s’écrit comme :
m1 m2
F~12 = −G ~u12 , (1.6)
r2
où r = M1 M2 est la distance entre les deux masses m1 et m2 , et ~u12 est le vecteur unitaire pointant de 1 vers 2.
Le facteur G est la constante de gravitation universelle :

G = 6, 67.10−11 (SI). (1.7)

7
La force de gravitation est toujours une attraction mutuelle (car les masses sont toutes positives), tandis que la
force électrostatique peut être une répulsion ou une attraction (car il existe des charges positives et des charges
négatives).
Exercice : Montrer qu’entre deux protons (ou entre un proton et un électron), la force électromagnétique
est beaucoup plus intense que la force gravitationnelle.

1.3 Distributions étendues de charges

Afin d’énoncer la loi de Coulomb, nous avons volontairement considéré la situation la plus simple où seulement
une paire de charges est présente. En pratique, les charges macroscopiques sont toujours composées d’un très
grand nombre de charges élémentaires. Nous avons donc besoin d’un formalisme permettant de décrire de grandes
collections de charges étendues sur des régions macroscopiques de l’espace. De nos jours nous savons que la
matière est granulaire à l’échelle atomique : les charges sont des électrons, des ions, des noyaux, etc.... Cependant
une description complète et exhaustive des positions et des mouvements de toutes ces charges élémentaires serait
impossible à réaliser. Il est donc nécessaire de passer à une description continue des charges. Cela est finalement
assez familier puisque pour résoudre des problèmes d’hydrodynamique (du remplissage de votre baignoire aux
écoulements dans des oléoducs) on utilise une description continue des fluides (la mécanique des fluides) au lieu
de suivre le mouvement de chaque molécule d’eau ou de gaz. Pour réaliser cela on utilise l’idée de lissage qui est
transversale dans de nombreux domaines de la physique.

Figure 1.2 – Distribution de charge 3D : Un volume chargé 3D macroscopique, noté V, peut être décomposé
en petits volumes dτ . La taille des cellules élémentaires dτ est petite par rapport au volume total de V, mais
suffisamment grande afin de contenir beaucoup de charges élémentaires (points noirs) comme on le voit sur le
zoom en haut à droite.

a) Densité volumique de charge

La stratégie consiste à réaliser une sorte de découpage ou de maillage du système macroscopique en un
grand nombre de petites cellules (Fig. 1.2). La taille des cellules doit être petite par rapport à l’échelle du
système macroscopique car elle fixe la résolution avec laquelle on pourra décrire des variations de densité au sein
du système. Cependant chacune de ces cellules doit encore contenir un grand nombre de particules sinon on
retombe sur l’aspect granulaire. On considère donc un élément de volume dτ , situé autour du point P , et de
taille intermédiaire : petit à l’échelle macroscopique et grand à l’échelle microscopique. Cet élément de volume,
−−→
repéré par le vecteur position ~r = OP (Figure 1.2), contient une charge dq. On appelle densité volumique de
charge en ~r la quantité :
dq
ρ(~r) = ρ(P ) = , (1.8)
dτ
qui s’exprime donc dans l’unité : C.m−3 . Attention, il vaut mieux voir la formule Eq. (1.8) comme le rapport
des deux nombres dq et dτ plutôt que comme la dérivée d’une fonction q qui dépendrait de la variable τ , ce qui
n’a pas de sens ici. Cette densité est par définition une quantité lissée (ou nivelée) qui décrit les variations de
répartition des charges dans un espace macroscopique avec une résolution de l’ordre de la taille de la cellule

8
élémentaire choisie. Le champ de densité de charge est un champ scalaire : à chaque point de l’espace (repéré
soit par ~r, soit par la lettre P ), on associe un nombre réel noté ρ(~r). Nous avons considéré des charge immobiles :
c’est pourquoi le champ de densité de charge ρ(~r) est statique (c’est à dire indépendant de la date t).
La charge totale Q dans un volume macroscopique V est donnée par l’intégrale volumique de la densité :
ZZZ ZZZ
Q= dq = ρ(~r)dτ (1.9)
V V

Au niveau conceptuel, il est préférable de voir ici l’intégration comme un procédé de sommation sur les petits
éléments de volume dτ plutôt que comme l’opération inverse de la dérivation. On a besoin de savoir sommer des
objects continus : c’est ce que fait l’intégrale triple pour des objets 3D (de même que l’intégrale double pour des
objets 2D, et l’intégrale curviligne à 1D).
Exercice : Comparer cette définition à la façon dont on peut décrire la répartition de masse (généralement
non uniforme) dans un objet macroscopique.
Exercice : Exprimer l’élément de volume dτ en utilisant les coordonnées sphériques (r, θ, ϕ). En déduire la
formule donnant le volume total d’une boule de rayon R.

b) Densité surfacique de charges

Considérons un objet chargé macroscopique modélisé par une nappe ou surface 2D (pas forcément plane).
On peut imaginer un feuillet de graphène ou un film très mince. De manière similaire au cas 3D, on va considérer
ici un élément de surface dS petit à l’échelle macroscopique et situé autour du point P repéré par le vecteur
−−→
position ~r = OP . Par ailleurs, on prend dS grand à l’échelle microscopique (c’est à dire contenant un très grand
nombre de charges élémentaires). Notons qu’ici dS est un nombre strictement positif qui mesure l’aire (en m2 )
de l’élément de surface. Cet élément de surface contient une charge dq proportionnelle à dS. On appelle densité
surfacique de charge en ~r la quantité :
dq
σ(~r) = σ(P ) = . (1.10)
dS
Il s’agit à nouveau d’un champ scalaire représentant une densité lissée, mais l’unité est maintenant le C.m−2 .
La charge totale Q dans une surface macroscopique S est donnée par l’intégrale surfacique de la densité :
ZZ ZZ
Q= dq = σ(~r)dS, (1.11)
S S

qu’il faut à nouveau voir comme une somme continue.

Exercice : Comparer cette définition à la façon dont on peut décrire la répartition de population dans un
espace géographique tel que l’Europe par exemple.
Exercice : Exprimer l’élément de surface dS dans un espace 3D, en utilisant les coordonnées sphériques
(r, θ, ϕ). En déduire la formule donnant la surface d’une sphère de rayon R.
Exercice : Exprimer l’élément de surface dS dans le plan, en utilisant les coordonnées polaires (r, θ). En
déduire la formule donnant la surface d’un disque de rayon R.

c) Densité linéique de charges

Considérons un objet chargé macroscopique modélisé par une ligne 1D (pas forcément droite, ni plane). On
peut imaginer un nanotube de carbone ou une molécule de polymère. On considère un élément de longueur
−−→
dl petit à l’échelle macroscopique et situé autour du point P repéré par le vecteur position ~r = OP . Par
ailleurs, on prend dl très grand à l’échelle microscopique (c’est à dire contenant un très grand nombre de charges
élémentaires). Cet élément de surface contient une charge dq. On appelle densité linéique de charge en ~r la
quantité :
dq
λ(~r) = λ(P ) = (1.12)
dl
Il s’agit à nouveau d’un champ de densité scalaire lissé, mais l’unité est maintenant le C.m−1 . La charge totale
Q porté par un lacet L est donnée par l’intégrale curviligne de la densité :
Z Z
Q= dq = λ(~r)dl (1.13)
L L

9
2 Charges mobiles : courant électrique
Le courant électrique est un mouvement d’ensemble de porteurs de charges, c’est à dire de particules
chargées. Généralement, il s’agit d’électrons libres dans les métaux. Il ne faut pas confondre ce mouvement
d’ensemble ordonné avec le mouvement d’agitation thermique désordonné et incessant des électrons. Le mouvement
d’agitation thermique existe toujours (même en l’absence de pile), il est très rapide mais se moyenne à zéro. En
revanche, le mouvement de dérive est entretenu par une pile. Si la vitesse de dérive est constante dans le temps,
on parle de courant continu (DC direct current en anglais). La première pile a été mise au point par Alessandro
Volta en 1800. Ici nous ne rentrons pas dans la description du fonctionnement de la pile. Nous introduisons
uniquement les notions nécessaires pour décrire le flot de porteurs de charge dans le circuit, à l’extérieur de la
pile.

2.1 Intensité du courant

On quantifie le courant par la quantité de charge qui traverse une section droite (de surface S) du conducteur
par unité de temps : c’est l’intensité I, mesurée en ampères (A). Ainsi si une charge δQ traverse la section droite
d’un conducteur pendant la durée dt, l’intensité du courant est :

δQ
I= (1.14)
dt
L’ampère correspond à une charge de 1 coulomb traversant la section du fil conducteur durant 1 seconde :

1 A = 1 C.s−1 (1.15)

Figure 1.3 – Fil rectiligne cylindrique de section S et d’axe Ox uniformément rempli d’électrons libres (porteurs
de charges) représentés par les points (bleus) avec la densité volumique n (en porteurs par m3 ). A gauche : Le
fil n’est pas connecté a une pile et les porteurs de charges sont immobiles (si on néglige un mouvement d’agitation
thermique désordonné). Le courant est nul : I = 0. A droite : On connecte le fil à une pile. Les porteurs de
charge acquièrent une vitesse de dérive commune ~vd = vd ~ux . Le courant est non nul : I = 6 0. L’image est prise à
la date t. Les porteurs (points en gras) dans le tronçon central de longueur vd dt sont ceux qui, entre les instants
t et t + dt, vont traverser la section coloré en (vert) foncé : leur nombre est donc en moyenne nvd dtS.

2.2 Densité de courant à 3D

On veut relier l’intensité I du courant aux caractéristiques microscopiques du flot de porteurs de charges
responsable de ce courant. On suppose que l’on a un fil rectiligne cylindrique de section S et d’axe Ox
uniformément rempli de porteurs de charges (typiquement des électrons) avec la densité n, en nombre de porteurs
par m3 (Figure 1.3). Ce gaz de porteurs est animé d’un mouvement d’ensemble et de dérive à la vitesse ~vd = vd ~ux .
Il faut compter les porteurs de charge qui traversent une certaine section droite entre les instants t et t + dt : il
s’agit des porteurs qui se situent à la date t dans le cylindre de section S et de longueur vd dt en amont de la
section considérée (Figure 1.3). La charge totale transversant la section S du fil pendant dt est donc donnée par :

δQ = q n (vd dt S), (1.16)

et l’intensité correspondante s’écrit :

I = nqvd S. (1.17)

10
Cette formule exprime que l’intensité est proportionnelle à la section du fil S et à une quantité que l’on va
appeler la densité de courant. On définit donc le vecteur densité de courant volumique :

~j(~r) = nq~vd , (1.18)

dont la norme s’exprime en C.s−1 .m−2 =A.m−2 . Nous avons ici un champ vectoriel ~j3D (~r) : on attache un vecteur
densité de courant à chaque point M (repéré par ~r) de l’espace.
On peut généraliser à une surface quelconque et définir le flux du vecteur densité de courant comme :

~j.−
→
ZZ
I= dS (1.19)
S

−
→
Notons qu’ici dS est bien une quantité vectorielle qui indique non seulement l’aire de l’élément de surface (la
−
→ −→
norme de dS), mais aussi l’orientation de cet élément de surface dans l’espace 3D (direction de dS) et enfin un
−
→
sens de traversée (sens de dS). En particulier, l’équation ci-dessus contient implicitement un choix de convention
−→
pour le signe des courants : on compte positivement les courants allant dans le sens indiqué par dS. C’est la
première fois de ce cours que nous voyons la notion de flux. Cette notion réapparaîtra plus tard dans le cours.
C’est aussi une notion très utile dans de nombreux domaines de la physique et des technologies : échanges
thermiques et flux de chaleur, démographie et flux migratoires, circulation routière et flux de voitures, etc...

2.3 Conducteurs de plus basse dimension

Notons pour terminer qu’il existe maintenant des conducteurs 2D ou même 1D. Dans ce cas il faut adapter les
notions de densité de courant. Depuis les années 1970, il existe des conducteurs 2D dans les transistors (au Si ou
AsGa) des ordinateurs et maintenant des téléphones portables. Depuis 2004, la forme 2D du carbone a été isolée :
c’est le graphène (Prix Nobel 2010), voir Fig. 1.4. Depuis les années 1980, il existe des conducteurs réellement
1D au niveau microscopique (conducteurs 1D organiques, nanotubes de carbone, fil quantiques semi-conducteur).

Figure 1.4 – A gauche : Echantillon de graphène suspendu (rectangle violet foncé), et connecté à 4 électrodes de
mesure (violet clair au-dessus de l’échantillon). A droite : Zoom au microscope électronique sur le structure
cristallographique (en nid d’abeille) du feuillet de graphène 2D.

3 Loi de conservation de la charge

Pour conclure ce chapitre nous allons discuter la notion de conservation de la charge. L’idée essentielle est
que si de la charge apparaît (ou disparaît) localement dans un certain petit volume, c’est que des charges sont
entrées (ou sont sorties) de ce volume sous forme de courants : aucune nouvelle charge n’a été créé spontanément
ni aucune charge n’a disparu. Il existe donc une relation entre la densité de charge volumique et les densités de
courants vues séparément aux deux paragraphes précédents.

11
Figure 1.5 – Bilan de charge : Fil rectiligne cylindrique de section S et d’axe Ox à la date t. En entrée et en
sortie du tronçon central de longueur dx, les intensités sont différentes donc la charge totale contenue dans le
tronçon va changer. Exercice : pour ce dessin, est-ce que cette charge va augmenter ou diminuer ? Les porteurs
de charge en-dehors du tronçon central n’ont pas été représentés pour alléger la figure (mais il y en a bien sûr...).

3.1 Cas d’un flot unidimensionnel

Pour simplifier, on commence par traiter le cas d’un courant uniforme dans la section du fil conducteur
(Figure 1.5). On note ρ(~r, t) = ρ(x, t) et ~j(~r, t) = ~j(x, t) en notant x la coordonnée le long du fil. Considérons
un petit segment compris entre les abscisses x et x + dx. A la date t, ce petit segment contient une charge
dq(t) = ρ(x, t)Sdx. A la date t + dt, ce petit segment contient une charge dq(t + dt) = ρ(x, t + dt)Sdx. Le
segment a donc gagné (en valeur algébrique) une charge :

∂ρ(x, t)
ρ(x, t + dt)Sdx − ρ(x, t)Sdx = dtSdx (1.20)
∂t
La conservation de la charge signifie qu’il ne se crée ni se détruit de charges, donc ce gain de charge correspond
à un bilan entre un courant de charges entrantes et un courant de charges sortantes. Entre les dates t et t + dt,
on peut compter que I(x, t)dt charges sont entrées en x, et que I(x + dx, t)dt charges sont sorties en x, donc le
gain (algébrique) :

∂jx (x, t)
I(x, t)dt − I(x + dx, t)dt = jx (x, t)Sdt − jx (x + dx, t)Sdt = − dxSdt, (1.21)
∂x
et en égalisant les deux gains on obtient :

∂ρ(x, t) ∂jx (x, t)

+ =0 (1.22)
∂t ∂x

Il s’agit d’une équation au dérivées partielles qui relie la dérivée temporelle de la charge locale avec la dérivée
spatiale du courant, que l’on appelle l’équation de conservation locale de la charge ou équation de
continuité.

3.2 Cas du flot tridimensionnel

Cette idée de conservation locale de la charge peut se généraliser à des flots réellement 3D dans lesquels
les charges peuvent aussi s’écouler selon les directions y et z, ce qui correspond aux projections du vecteur
densité de courant ~j(~r, t) = ~j(~r, t) suivant les axes y et z. En faisant le bilan sur un petit cube situé entre x et
x + dx, entre y et y + dy et entre z et x + dz, on obtient l’équation de conservation locale de la charge
ou équation de continuité à 3D :

∂ρ(~r, t) ∂jx (~r, t) ∂jy (~r, t) ∂jz (~r, t)

+ + + = 0, (1.23)
∂t ∂x ∂y ∂z

~ = (∂/∂x, ∂/∂y, ∂/∂z) :

qui s’écrit en utilisant l’opérateur différentiel ∇

∂ρ(~r, t) ~ ~
+ ∇.j(~r, t) = 0 . (1.24)
∂t

12
4 Conclusion
Nous avons décrit la force d’interaction électromagnétique entre les charges immobiles. Puis nous avons
expliqué que le courant électrique consiste en un mouvement d’ensemble des charges mobiles. Les notions clefs
pour décrire les charges sont la densité scalaire de charge électrique ρ, et le vecteur densité de courant électrique
~j. Dans les chapitres suivants, nous allons montrer que les charges immobiles génèrent des champs électriques :
c’est le domaine de l’électrostatique. Ensuite, nous reviendrons sur les charges en mouvement et montrons que
les courants génèrent des effets magnétiques, de même que les aimants.

A retenir :
— La loi de Coulomb : force d’interaction en 1/r2 entre deux charges ponctuelles et immobiles.
— Les courants sont constitués de charges en mouvement.
— La valeur de la charge élémentaire e = 1, 6 . 10−19 C. Electron de charge −e et proton de charge +e.
— Savoir utiliser les densités de charge et de courant.

Les unités :
— L’unité de charge électrique : le coulomb (C).
— L’unité de l’intensité d’un courant électrique : l’ampère (A).
— L’unité de densité volumique de courant : A.m−2 .
— Savoir en déduire par raisonnement d’homogénéité les unités des différentes densités 3D, 2D, 1D.

13
Chapitre 2

Le champ électrostatique

Au chapitre précédent, nous avons écrit la loi de Coulomb qui décrit la force à distance entre deux charges
ponctuelles immobiles. Le but de ce chapitre est d’introduire la notion de champ électrique, en se concentrant
d’abord sur la situation la plus simple possible, c’est à dire le champ créé par une charge unique. Nous introduirons
le principe de superposition qui permet de calculer le champ électrique créé par une distribution étendue de charges
à partir de la connaissance du champ créé par une simple charge source ponctuelle. Enfin, nous montrerons
comment l’analyse des symétries d’une distribution de charge permet d’obtenir de précieux renseignements sur la
cartographie des lignes de champ électrique.

1 Champ électrique créé par une charge ponctuelle unique

On peut considérer qu’une charge q1 seule et non nulle dans l’univers modifie les propriétés de tout l’espace
environnant. En effet une seconde charge q2 (également non nulle) passant à proximité subira une force électrique
de Coulomb et sera donc déviée par la première charge au lieu de suivre un mouvement rectiligne uniforme (ce
qui se passerait en vertu du principe d’inertie si q1 était nulle). On dit que la charge q1 créé un "champ de
forces" qui baigne tout l’espace. Le champ électrique correspond à ce champ de force divisé par q2 . La dimension
du champ électrique est donc le newton par coulomb (N.C−1 ), ce qui aussi équivalent au volt par mètre (V.m−1 ).
Dans le (SI), l’unité officielle de champ électrique est le V.m−1 .

1.1 Notion de champ électrique

On revient sur une situation très simple : 2 charges ponctuelles immobiles, situées en M1 et M2 respectivement.
On a vu au chapitre précédent qu’il existe une action mutuelle décrite par une force à distance. La force exercée
par la charge 1 sur la charge 2 est proportionnelle au produit q1 q2 , et s’écrit comme :
→
− 1 q1 q2
F 12 = ~u12 , (2.1)
4π0 r2
où r = M1 M2 est la distance entre les deux charges q1 et q2 et ~u12 est le vecteur unitaire pointant de la charge 1
vers la charge 2.
→
−
Pour obtenir une quantité intrinsèque attachée à la charge q1 , on divise la force F 12 par la charge q2 :
→
−
→
− F 12 1 q1
E 1 (M2 ) = = ~u12 , (2.2)
q2 4π0 r2
→
−
où on dira que E 1 (M2 ) est le champ électrique créé par la charge q1 en M2 .
1. La charge 1 (appelée charge source) modifie les propriétés de l’espace qui l’entoure, en créant un champ,
→
− →
−
noté E 1 (M ), qui s’étend dans tout l’espace vide autour d’elle. Dans la notation E 1 (M ), l’indice 1 nous
rappelle que c’est la charge 1 qui créé le champ, tandis que M représente un point arbitraire qui peut
décrire tout l’espace (en dehors de M1 où ce champ diverge). On voit que le champ "baigne" tout l’espace :
c’est un objet complètement nonlocal bien qu’il soit créé par une charge ponctuelle unique (la source la
plus locale qui soit).
2. Si une charge (appelée charge test) est placée en un point donné M2 de l’espace, elle va se coupler
→
−
localement le champ E 1 (M2 ) créé en M2 , et subir une force égale à sa charge q2 fois le champ.

14
 1.0 1.0

0.5 0.5

0.0 0.0

-0.5 -0.5

-1.0 -1.0

-1.0 -0.5 0.0 0.5 1.0 -1.0 -0.5 0.0 0.5 1.0

Figure 2.1 – Le champ électrostatique créé par une charge ponctuelle placée à l’origine du plan. Il s’agit d’une
coupe, en réalité il y a une symétrie sphérique complète autour du point origine.

On remplace ainsi l’action à distance entre charges discrètes par un couplage local entre une charge et le champ
électrique à l’endroit où se trouve cette charge. Bien sûr la non localité se retrouve dans le fait que la charge
source localisée (et même ponctuelle) crée un champ dans tout l’espace.
Dans le paragraphe suivant, ce point de vue sera généralisé à des situations plus complexes de la manière
suivante : admettons qu’une distribution très compliquée de charges (éventuellement inconnue) créé un certain
champ E(M~ ) en tout point M d’une certaine région de l’espace. Si on place une charge additionnelle (la charge
~
test) q, elle subira une force q E(M ). On peut aussi dire que le champ électrique est le champ de force électrique
par unité de charge électrique test.

1.2 Lignes de champ électrique

Revenons à notre situation fondamentale avec une seule charge q source immobile. Il est judicieux de choisir
l’origine des coordonnées sphériques sur cette charge. On utilisera les coordonnées sphériques (r, θ, ϕ) pour
repérer un point arbitraire M de l’espace. Le champ en M s’écrit alors :

→
− 1 q
E (M ) = ~ur (2.3)
4π0 r2

On peut représenter le champ en traçant le vecteur en plusieurs points de l’espace autour de la charge source
→
−
(Fig. 2.1). Il faut bien noter que le champ E (M ) dépend des trois coordonnées du point M (et pas seulement de
r) car le vecteur unitaire ~ur porte une dépendance en θ et en ϕ. Cependant on doit observer ici les invariances
suivantes concernant d’une part la charge source, et d’autre part le champ qu’elle engendre :
1. La charge source est invariante par toute rotation autour de l’origine. Ces rotations sont quantifiées par
deux angles (qui peuvent être par exemple les angles θ et en ϕ des coordonnées sphériques.)
2. Le champ vectoriel électrostatique Eq. (2.3) est invariant par toute rotation autour de l’origine. Cela
signifie que ce champ ne distingue aucune direction privilégiée de l’espace : il est parfaitement isotrope.
Cela reste vrai pour une charge source consistant en une boule chargée uniformément.

2 Champ électrique total créé par une distribution de charges : principe de

superposition
Dans les paragraphes précédents, nous avons principalement traité du champ électrostatique créé par la source
la plus simple : une charge ponctuelle immobile. Nous passons maintenant au cas de distributions étendues de
charges (toujours immobiles). De telles distributions sont en fait constituées de milliards de charges élémentaires.
On introduit ici l’idée très importante (et justifiée par l’expérience) de principe de superposition. En vertu du
principe de superposition, on peut calculer le potentiel électrostatique créé par une distribution étendue de
charge en décomposant celle-ci en charges élémentaires ponctuelles ou infinitésimales.

15
2.1 Distribution discrète de charges ponctuelles

Soit une collection de particules chargées qi placées en divers points Pi de l’espace (i = 1, 2, 3, ..). Imaginons
que l’on place une charge test q (fictive ou réelle) en un point M . Chaque charge qi exerce sur la charge test q
une force d’interaction Colombienne à distance :
→
− 1 qqi
Fi = ~ui , (2.4)
4π0 ri2
−−→
où ri = Pi M , et ~ui = Pi M /ri est le vecteur unitaire dirigé de Pi vers M .
La force totale subie par la charge test de la part de l’ensemble des charges sources est donc :
→
− X→ − X qi
F = Fi =q ~ui . (2.5)
i i
4π0 ri2
→
−
Pour obtenir le champ électrique créé par les charges sources en M , il suffit de diviser la force totale F par
la charge test q. On obtient alors :
→
− 1 X qi ~ui
E (M ) = , (2.6)
4π0 ri2i
−−→
avec ri = Pi M , et ~ui = Pi M /ri le vecteur unitaire dirigé de Pi vers M . Les charges sources sont immobiles et
situées aux points Pi , tandis que M est un point arbitraire d’observation. Le champ total en M , du aux N
charges ponctuelles, est égal à la somme des champs créés individuellement par chacune des charges sources.

2.2 Distributions continues de charges

Figure 2.2 – Distribution 3D : Un volume chargé 3D, noté V, peut être décomposé en petits volumes dτ . La
taille des cellules élémentaires dτ est petite par rapport au volume total de V, mais suffisamment grande afin de
contenir beaucoup de charges élémentaires (points noirs) comme on le voit sur le zoom en haut à droite.

Distributions 3D :
Dans une description continue en terme de densité volumique de charge, il faut remplacer la somme discrète
du paragraphe précédent par une somme continue, c’est à dire une intégrale triple :

→
−
ZZZ
1 ρ(P )dτ
E (M ) = ~uP M . (2.7)
4π0 P ∈V PM2

Notons que dans le passage au continu Pi devient le point "courant" P sur la distribution source, et que l’on
applique le principe général de passage au continu élaboré au chapitre 2.
X ZZZ
qi (...) → dτ ρ(P )(...) (2.8)
i V

16
Distributions 2D :
Dans une description continue, il faut utiliser une densité surfacique de charge, et remplacer la somme
discrète par une somme continue 2D, c’est à dire une intégrale double :

→
−
ZZ
1 σ(P )dS
E (M ) = ~uP M . (2.9)
4π0 P ∈S PM2

On a appliqué le principe général de passage au continu sous la forme :

X ZZ
qi (...) → dSσ(P )(...). (2.10)
i S

Distributions 1D :
Dans une description continue, il faut utiliser une densité linéique de charge, et remplacer la somme discrète
par une somme continue 1D, c’est à dire une intégrale curviligne :

→
−
Z
1 λ(P )dl
E (M ) = ~uP M . (2.11)
4π0 P ∈L PM2

On a appliqué le principe général de passage au continu sous la forme :

X Z
qi (...) → dlλ(P )(...). (2.12)
i L

Le champ infinitésimal créé par le segment dl (placé en P) est :

−→ 1 λ(P )dl
dE P (M ) = ~uP M . (2.13)
4π0 P M 2

3 Mouvement d’une charge dans un champ électrique uniforme

→
−
Nous avons vu que le champ électrique E (~r) est défini par son action mécanique sur une charge "test".
Pour que la description en terme de champ d’interaction ait tout son sens, il faut que cette force ne dépende
que de la valeur du champ sans qu’il soit nécessaire de préciser quelles sont les dispositifs et sources
qui ont engendré ce champ : l’expression de la force ne doit contenir que le champ lui-même et certaines
caractéristiques de la charge test (ici la charge électrique q). C’est justement le cas pour la force électrique
exercée sur une charge q :

→
− →
−
F = q E (~r) . (2.14)
→
−
Si on donne la forme de E (~r), on peut alors étudier le mouvement de la charge q, sans savoir quelles charges
ont engendré le champ. Cela suppose tout de même que la charge test q introduite pour définir le champ ne
perturbe pas la distribution de charges sources.
→
−
L’exemple le plus simple est le cas du mouvement dans un champ uniforme E , c’est à dire un champ
indépendant du point de l’espace où se situe la particule. Il n’existe pas de champ uniforme à l’échelle de
l’Univers entier bien sûr, mais on peut fabriquer des champs électriques uniformes dans une certaine région de
l’espace. Le dispositif typique est le condensateur plan qui permet de générer des lignes de champ uniforme
entre deux armatures planes portées à une certaine différence de potentiel U .
→
−
Prenons un champ électrique uniforme E = E0 u~y , et une particule de charge q avec une vitesse ~v (t = 0) =
v0 ~ux . On prend l’origine des coordonnées cartésiennes à la position initiale de la particule à t = 0 : ~r(t = 0) = ~0
(Fig. 2.3).
L’équation du Newton est donc :
d~v (t) →
−
m = qE . (2.15)
dt
C’est l’équation d’un mouvement à une accélération constante, donc analogue à l’équation de la chute libre d’un
objet dans le champ de pesanteur terrestre.

17
Figure 2.3 – Mouvement d’un électron accéléré par un champ électrique uniforme :Une charge q,
négative sur ce dessin, se déplace dans un champ électrique uniforme, sous l’influence de la seule force électrique.
La trajectoire est une parabole.

En projection, l’équation de Newton devient :

dvx (t) dvz (t) dvy (t) qE0

= = 0, et = , (2.16)
dt dt dt m
que l’on peut intégrer très simplement en :

qE0
vx (t) = v0 , vz (t) = 0, et vy (t) = t. (2.17)
m
On prend les primitives en tenant compte des conditions initiales :

qE0 t2
x(t) = v0 t, vz (t) = 0, et y(t) = . (2.18)
2m
Ce sont les équations horaires du mouvement. La trajectoire est une parabole d’équation :

qE0 x2
y= . (2.19)
2mv02

On voit que la courbure de la parabole (en son sommet O) est d’autant plus grande que le champ est grand
et la vitesse initiale faible (Fig. 2.3).

4 Champ électrique et symétries spatiales

Considérons une situation physique réelle, notée 1, décrite par la dynamique suivante :

d~v1 (t) ~ 1 (~r1 (t)) ,

m = qE (2.20)
dt
et correspondant à une certaine trajectoire ~r1 (t) bien réelle de la particule.
Maintenant imaginons la même particule (charge q et masse m) qui décrirait la trajectoire obtenue en
appliquant une symétrie miroir par rapport à un plan Π à la trajectoire 1. La trajectoire ainsi obtenue correspond
à une situation hypothétique que nous allons noter 2. A chaque instant, par construction géométrique, la position
sur cette trajectoire 2 serait donnée par :

~r2 (t) = SymΠ [~r1 (t)] , (2.21)

18
et donc la vitesse et l’accélération se transformeraient aussi comme :
 
d~v1 (t) d~v2 (t) d~v1 (t)
~v1 (t) → ~v2 (t) = SymΠ [~v1 (t)], et → = SymΠ . (2.22)
dt dt dt

~ 2 (~r) qui conduirait

On peut se poser la question suivante : existe-il une cartographie de champ magnétique B
exactement à la cinématique (trajectoire 2) hypothétique construite par symétrie miroir de la situation réelle ?
Si un tel champ existe, il doit satisfaire à tout instant à l’équation de Newton :

d~v2 (t) →
−
m = q E 2 (~r2 (t)) . (2.23)
dt
En comparant cette équation avec l’équation initiale Eq. (2.20), on voit qu’il faut qu’à chaque instant :
→
− →
− →
−
E 1 (~r1 (t)) → E 2 (~r2 (t)) = SymΠ [ E 1 (~r1 (t))] . (2.24)
→
−
Et si l’on veut que cela soit vrai pour toute condition initiale, il faut que le champ statique E 2 (~r) satisfasse :
→
− →
−
E 2 (~r2 ) = −SymΠ [ E 1 (~r1 )] . (2.25)

On voit que, sous une opération de symétrie miroir, le vecteur champ magnétique se transforme comme les
vecteurs position, vitesse ou accélération. On dit que ces derniers sont des vecteurs polaires.

19
Chapitre 3

Les propriétés du champ électrostatique

Nous démontrons et discutons les deux propriétés globales du champ électrique. D’une part, la circulation du
champ électrique le long d’un chemin fermé est toujours nulle, ce qui signifie aussi que le champ électrostatique
dérive d’un potentiel. La seconde propriété est exprimée par le théorème de Gauss qui relie le flux du champ
électrique aux charges sources. Ces deux propriétés peuvent être traduites et synthétisées sous la forme de deux
équations aux dérivées partielles (les équations de Maxwell de l’électrostatique).

1 Circulation et potentiel électrostatique

1.1 Circulation du champ électrique
Le champ est un objet étendu. Pour obtenir des informations globales
sur le champ on peut étudier ce que l’on appelle sa circulation le long d’un
chemin 1D abstrait quelconque. Il s’agit de la quantité suivante :
Z Z
CL = ~
E(M ~M =
).dl E. ~
~ dl, (3.1)
M ∈L L

~
dans laquelle on fait la somme des produits scalaires E(M ~ M le long d’un chemin L reliant le point de départ
).dl
A au point d’arrivée B. En général cette circulation dépend de la carte du champ mais aussi du chemin suivi.
Nous allons voir que le champ électrostatique est très particulier car pour lui la circulation de dépend que des
points de départ et d’arrivée.
Une propriété très importante du champ électrostatique est que sa circulation ne dépend que des points de
départ et d’arrivée et non du chemin suivi. Pour le montrer, nous allons revenir au cas le plus simple, celui du
champ créé par une charge ponctuelle source unique (située à l’origine). Puis nous généraliserons au cas d’une
distribution quelconque en utilisant le principe de superposition.

1.2 Potentiel électrique créé par une charge unique

Nous pouvons réécrire la champ électrostatique Eq.(2.3) comme :
 
~ ∂ 1 q ∂V
E(M )=− ~ur = − ~ur , (3.2)
∂r 4π0 r ∂r
avec
1 q
V (M ) = V (r) = + K. (3.3)
4π0 r
La quantité V (M ) est appelée potentiel électrostatique. Elle est définie à une constante additive près (notée K
ci-dessous) car seul le champ est mesurable, ou bien des différences de potentiel. On prendra dans le cas de la
charge unique source : K = 0. Donc :
1 q
V (M ) = V (r) = . (3.4)
4π0 r

Rq 1 : Puisqu’on attache une valeur du potentiel V (M ) à chaque point de l’espace M , le potentiel est un
champ scalaire (au sens mathématique), tout comme le champ de densité de charge.

20
 Rq 2 : On définit les lignes ou surfaces caractérisées par une même valeur du potentiel électrique : les
équipotentielles. Dans le cas d’une charge ponctuelle source, le lieu de tout les points tel que V (M ) = V0
correspond aux points tel que r = cste, il s’agit donc de sphères centrées sur la charge source. Il faut noter que
les lignes de champs (radiales) sont toujours perpendiculaires aux équipotentielles.
Rq 3 sur les unités : volt V et V.m−1 .
On peut maintenant évaluer la circulation du champ électrostatique de la charge isolée le long d’un chemin
quelconque L reliant un point de départ A à un point d’arrivée B :
Z B
→
− →
− →
−
Z Z Z
CL = ~
E. dl = − ∇V · dl = − dV = − dV = V (A) − V (B), (3.5)
L L L A
et on voit que le résultat de dépend que la différence de potentiel électrostatique entre les deux points extrêmes
A et B, et non du chemin suivi entre ces points. Notons aussi que dans l’égalité ci-dessus, on a seulement utilisé
le fait que le champ électrostatique dérivé d’un potentiel (et non l’expression spécifique de ce potentiel).
Sur un chemin fermé (pour lequel A = B sont confondus) la circulation est donc nécessairement nulle :

~ →−
I
E. dl = V (A) − V (A) = 0 , (3.6)
L

pour toute boucle fermée L. Ceci est une propriété typique et remarquable du champ électrostatique (un champ
vectoriel quelconque ne possède pas cette propriété en général).

1.3 Lien avec les notions de travail et d’énergie potentielle

Imaginons que l’on promène une charge test q0 le long du chemin L. Lors de ce déplacement, la charge test
~
subit une force électrostatique q0 E(M ) dont le travail le long de L s’écrit :
Z
WL = ~
q0 E(M ~ M = q0 C L .
).dl (3.7)
M ∈L
On voit que la travail (Work en anglais) est égal à la circulation du champ multiplié par la charge test (celle qui
travaille).
Remarque sur les signes. Le travail hérite donc automatiquement de la propriété d’indépendance par rapport
au chemin suivi :

WL = q0 V (A) − q0 V (B), (3.8)

Lien avec la mécanique : on peut écrire le théorème de l’énergie cinétique entre A et B
Ec (B) − Ec (A) = WL = q0 V (A) − q0 V (B), (3.9)
On peut réécrire cela en regroupant :
Ec (B) + q0 V (B) = Ec (A) + q0 V (A). (3.10)
Cette dernière égalité fait apparaître une nouvelle quantité qui se conserve : la somme de l’énergie potentielle et
de q0 V . Cette quantité conservée s’appelle l’énergie totale de la particule test. La partie q0 V s’appelle l’énergie
potentielle et traduit l’interaction de la charge test q0 avec le champ électrostatique extérieur (représenté par le
potentiel V ).

2 Potentiel électrostatique pour les distributions de charges

2.1 Distribution discrète de charges ponctuelles

Soit une collection de particules chargées qi placées en divers points Pi de l’espace (P1 ,P2 ,P3 ). Imaginons
que l’on place une charge test q (fictive ou réelle) en un point M . Chaque charge qi exerce sur la charge test q
une force d’interaction Colombienne à distance :
Le potentiel électrique créé par les charges sources en M est la somme des potentiels créés par chaque charge
qi :
1 X qi
V (M ) = (3.11)
4π0 ri
i
avec ri = Pi M . Les charges sources sont immobiles et situées aux points Pi , tandis que M est un point arbitraire
d’observation.

21
Figure 3.1 – Flux du champ électrostatique E(M ~ ) à travers une surface fermée Σ. Le flux infinitésimal
~ ~
à travers une surface élémentaire s’écrit E(M ) · dS M par définition, où dS~ M est le vecteur surface orienté vers
l’extérieur de la surface Σ. On obtient le flux sortant total en additionnant tous les flux élémentaires de manière
à recouvrir toute la surface Σ. Les sources (charges) qui créent le champ électrostatique E(M ~ ) ne sont pas
représentées sur ce dessin.

2.2 Distributions continues de charges

Distributions 3D :
Dans une description continue en terme de densité volumique de charge, il faut remplacer la somme discrète
du paragraphe précédent par une somme continue, c’est à dire une intégrale triple :
ZZZ
1 ρ(P )dτ
V (M ) = . (3.12)
4π0 P ∈V PM

où ρ(P ) est la densité de charge volumique en P .

Distributions 2D :
Dans une description continue, il faut utiliser une densité surfacique de charge, et remplacer la somme
discrète par une somme continue 2D, c’est à dire une intégrale double :
ZZ
1 σ(P )dS
V (M ) = , (3.13)
4π0 P ∈S PM2

où σ(P ) est la densité de charge surfacique en P .

Distributions 1D :
Dans une description continue, il faut utiliser une densité linéique de charge, et remplacer la somme discrète
par une somme continue 1D, c’est à dire une intégrale curviligne :
Z
1 λ(P )dl
V (M ) = , (3.14)
4π0 P ∈L PM

où λ(P ) est la densité de charge linéique en P .

3 Flux du champ électrique et théorème de Gauss

Soit un champ électrostatique défini en tout point de l’espace 3D (sauf peut-être sur à l’emplacement des
charges ponctuelles isolées). Pour obtenir des informations globales sur le champ on peut étudier ce que l’on
appelle le flux du champ électrique à travers une surface 2D abstraite quelconque, notée ici Σ. Il
s’agit de la quantité suivante (Fig. 3.1) :
→
− −
→
ZZ
ΦΣ = E (M ) · dS M , (3.15)
M ∈Σ

22
Figure 3.2 – Théorème de Gauss : Sphère de Gauss Σ centrée sur la charge source : a) positive (à gauche)
b) négative (à droite).

−
→
dans laquelle le vecteur dS M est le vecteur surface élémentaire associé à la surface dS autour du point M ,
défini comme il suit :
−
→
— dS M est normal à Σ en M ,
−
→
— dS M est orienté vers l’extérieur à Σ,
−
→
— dS M a pour norme l’aire dS.
La surface Σ est une nappe 2D immergée/plongée dans l’espace 3D habituel (exemples : une sphère, une
surface cylindrique, un parallélépipède...) ; elle peut être ouverte ou fermée. Cependant dans la suite, nous allons
nous intéresser uniquement aux surfaces fermées car elle présente c’est à travers des surfaces fermées que le
flux du champ électrostatique présente une propriété remarquable. Un surface 2D fermée sépare l’espace 3D en
deux zones : l’intérieur et l’extérieur.
Nous présentons ici la seconde propriété générale du champ électrique. Elle est vraie pour tout champ
électrostatique. Avant de la démontrer, nous annonçons cette propriété :
ZZ
~ ~ Qint
E(M ).dS(M )= , (3.16)
M ∈Σ 0

où Qint est la charge à l’intérieur de la surface fictive Σ.

3.1 Cas du champ créé par une seule charge ponctuelle source à travers une sphère
centrée
Nous allons travailler sur le cas particulier le plus simple : une charge source q placée à l’origine. On va
imaginer par la pensée une sphère de rayon R que l’on peut déplacer dans tout l’espace, et on va se demander
quel est le flux du champ à travers cette surface fermée. Pour fixer les choses, on va même la charge q positive
(exercice : refaire le même raisonnement et les dessins avec une charge négative). On voit que :
1. Si la sphère contient la charge : tout le flux est sortant. Le flux sera donc non nul.
2. Si la sphère ne contient pas la charge : le flux est en partie sortant et en partie est entrant. La question
naturelle : est-ce que ces flux antagonistes se compensent exactement pour donner un flux nul ? On peut
essayer de deviner la réponse à partir du dessin.

3.2 Cas du champ créé par une seule charge ponctuelle source à travers une surface
quelconque
Ici on généralise le résultat précédent en autorisant des déformations de la surface fictive ΣG .
~
Le flux élémentaire de E(M ~
) à travers la surface élémentaire dS(M ) est donné par :
−→
~ ~ M = q ~ur · dS M q
dΦ = E(M ).dS 2
= dΩ. (3.17)
4π0 r 4π0

23
Figure 3.3 – Théorème de Gauss : Sphère de Gauss Σ excentrée par rapport à la charge source : a) sphère
(à gauche) b) "sphère" déformée (à droite).

~ M depuis
Dans la dernière égalité, dΩ est l’angle solide (en steradian) sous lequel est vu la surface élémentaire dS
la charge (située en O).
L’angle solide sous lequel est vu une surface depuis un point d’observation O est la projection de cette
surface sur une sphère de rayon unité entourant et centrée sur O. C’est la généralisation 2D de la notion de
taille angulaire (en radian) d’un objet 1D qui est la projection de cet objet sur le cercle trigonométrique de
rayon unité centré sur O.
Si l’observateur est entièrement enveloppé par ΣG , celle-ci se projettera sur la totalité de la sphère unité,
donc l’angle solide sera égal à 4π steradians.
Si l’observateur "regarde" la surface fermée ΣG depuis l’extérieur, il y a un phénomène de compensation
algébrique des contributions positives et négatives dΩ, qui donne un angle solide nul à la fin.

3.3 Cas général

En invoquant le principe de superposition, on peut décomposer toute distribution en N charges ponctuelles,
et le théorème de Gauss est vrai pour chacune d’elle.
ZZ
~ ~ Qint
E(M ).dS(M )= , (3.18)
M ∈Σ 0

où Qint est la charge à l’intérieur de la surface fictive Σ. Pour le cas d’une charge unique source : Qint = q si la
charge source est à l’intérieur et Qint = 0 si la charge est à l’extérieur.

4 Exemple d’utilisation du théorème de Gauss : boule chargée

Le théorème de Gauss est une propriété générale de tout champ électrique. Dans la plupart des cas, ce
théorème ne donne que le flux total qui est donc une sorte de somme des champs électriques perçant une surface
fermée, et il n’est pas possible de "remonter" aux valeurs du champ. En effet, si on donne la valeur d’une somme,
il est en général impossible de remonter aux valeurs individuelles de chaque terme... sauf si tous les termes sont
égaux et qu’on connaît leur nombre. Dans le cas de distributions de charges suffisamment symétriques, il sera
possible d’utiliser les symétries + le théorème de Gauss pour déterminer complètement les valeurs du champ
électrique.
Considérons le cas d’une boule uniformément chargée de rayon R, et portant la charge totale Qtot . La densité
volumique de charge ρ est donc :
Qtot
ρ = 4πr3 , (3.19)
3

1) Analyse des symétries pour obtenir la direction du champ : Tout plan ΠS passant par le centre
O de la boule et le point M d’observation, est un plan de symétrie pour la boule. C’est aussi un plan de symétrie
des charges puisque la densité est uniforme dans la boule. Par conséquent, comme le champ électrique est un
→
−
vecteur polaire, E (M ) doit être contenu dans tous ces plans de symétrie ΠS , et donc dans leur intersection : la
droite (OM ). Le champ électrique est donc radial :
→
−
E (M ) = Er (r, θ, ϕ) ~ur . (3.20)

24
Figure 3.4 – Carl Friedrich Gauss (1777-1855) : Fils unique d’un couple d’ouvriers, Carl Friedrich Gauss est
un génie très précoce dont les études seront financées par le duc de Brunswick. Il étudie par lui-même les oeuvres
de Newton, Euler et Lagrange avant d’entrer à l’université de Göttingen à 18 ans. En mathématiques, Gauss
a des contributions très importantes en théorie des nombres (arithmétique), en algèbre (théorème fondamental
de l’algèbre), en probabilité (méthode des moindres carrés et fonction "en cloche" gaussienne), en géométrie
non-euclidienne (son élève Riemann prolongera ses travaux), en géométrie différentielle,... Nommé professeur
d’astronomie à l’Université de Göttingen, il améliore les techniques de calcul des orbites, fait bâtir l’observatoire
de Göttingen, et se consacre aussi à la géodésie. A partir de 1830, il collabore avec Weber sur des sujets de
physique reliés au magnétisme et à l’électricité : théorème de Gauss et théorie du potentiel.

De plus cette composante radiale ne peut pas dépendre des angles ϕ et θ à cause de l’invariance par les rotations
autour du centre de la boule O. Donc :
→
−
E (M ) = Er (r) ~ur . (3.21)
Ainsi, au lieu d’avoir à déterminer 3 composantes chacune fonction des 3 coordonnées du point M , nous avons
seulement à déterminer une fonction d’une seule variable Eθ (r).
2) Utilisation du théorème de Gauss : comme calculer le flux d’une fonction encore inconnue Eθ (r) ?
L’astuce consiste à choisir une surface d’intégration sur laquelle cette fonction inconnue ne varie pas, donc une
surface avec r constante. Cela dicte le choix de la surface de Gauss (Σ) comme étant une sphère de centre O et
de rayon r (passant par M ). Le flux est donc :
ZZ ZZ ZZ
~ ~
E(M ).dS(M ) = Er (r) ~ur . dS~ur = Er (r) dS = 4πr2 Er (r), (3.22)
M ∈Σ Σ Σ

où on a pu sortir Eθ (r) en dehors de l’intégrale double car r est constante sur Σ. Par ailleurs la somme des dS
donne la surface totale de la sphère de rayon r, soit 4πr2 .
Le théorème d’Ampère donne donc :

Qint (r)
4πr2 Er (r) = , (3.23)
0

où Qint (r) désigne la charge totale emprisonnée dans une sphère fictive de rayon r.
On distingue deux cas :
• A l’extérieur de la boule r > R : La fonction Qint (r) sature à la valeur Qtot indépendante de r, dès
que r dépasse R. Donc :
Qtot
Er (r) = , (3.24)
4π0 r2
qui est un résultat remarquable. Cela signifie qu’à l’extérieur de la boule, le champ est le même que si
toute la charge avait été concentrée en un point particulier : le centre O de la boule.
• A l’intérieur de la boule r < R : le fonction Qint (r) est donnée par :

4πr3 r3
Qint (r) = ρ = 3 Qtot , (3.25)
3 R

25
 et donc le champ est :
ρr
Er (r) =. (3.26)
30
La norme du champ croît linéairement en r. En particulier le champ électrique s’annule au centre de la boule
conformément aux arguments de symétrie.

5 Equations de Maxwell de l’électrostatique

Les propriétés globales (de circulation et de flux) du champ électrique peuvent être exprimées sous la forme
d’équations aux dérivées partielles.

→
− ~ →
− ~ ρ(~r)
∇ ∧ E(~r) = ~0 , ∇ · E(~r) = . (3.27)
0

1) Ces équations relient les dérivées spatiales du champ électrostatique entre elles et aux charges électriques
sources. Il s’agit de situations où toutes les charges sont immobiles : les courants sont donc nuls ainsi que
le champ magnétique. En régime statique, toutes les dérivées temporelles ∂/∂t sont nulles, conduisant à
des équations fermées pour le champ électrique et les charges sources.
2) La première équation est homogène : elle exprime que le champ électrostatique dérive d’un potentiel.
3) La seconde équation relie le champ à ses sources (démontrer le théorème de Gauss à partir de cette
équation locale).
4) Ecrire l’équation de Poisson de l’électrostatique et la solution correspondante.
5) En déduire la loi de Coulomb exprimant le champ électrique en fonction des charges.
6) La loi de conservation de la charge est vérifiée trivialement : ρ(~r) est indépendant du temps et ~j(~r) = ~0.

6 Exercices

~ r) par le théorème de Gauss

Exercices sur le calcul du champ E(~
1) Il faut toujours préalablement trouver la direction des lignes de champ et les invariances du
champ en se basant sur l’analyse des symétries des charges. Ensuite il faut choisir une surface fermée
fictive passant par M afin d’appliquer le théorème de Gauss.

2) Calculer le champ électrique créé par une sphère uniformément chargée.

3) Calculer le champ électrique créé par une boule uniformément chargée.

4) Calculer le champ électrique créé par un cylindre (infini) creux et uniformément chargé.

5) Calculer le champ électrique créé par un cylindre plein (infini) et uniformément chargé.

6) Calculer le champ électrique créé par une nappe plane infinie uniformément chargée. En déduire le
champ entre les armatures d’un condensateur plan.

CQFR
— Formule du champ créé par une charge isolée (autre forme de la loi de Coulomb) : : en 1/r2
— Formule du potentiel créé par une charge isolée : en 1/r
~ →
−
— Le champ électrique dérive d’un potentiel : E(M ) = − ∇V (M ).
— Théorème de Gauss
— Le principe de superposition
— Savoir dessiner les lignes de champ électrique créé par une charge ponctuelle

26
Unités et conventions :
— L’unité de potentiel électrique : le volt (V)
— L’unité de champ électrique : le volt par mètre (V.m−1 )
— Savoir orienter une surface et relier cela au signes possibles du flux à travers cette surface

27
Chapitre 4

Le champ magnétique

Depuis l’Antiquité, on avait déjà remarqué les propriétés fascinantes des aimants. Au Moyen-Age, la boussole
permet d’aider à la navigation et on réalise que la Terre agit comme une sorte de très gros aimant qui interagit
avec la boussole (Gilbert 1600). En 1820, Oersted fait une découverte fondamentale : un courant électrique dévie
la boussole ! C’est la rencontre et bientôt la réunion de deux domaines : l’électricité et le magnétisme. Ampère va
explorer systématiquement ce type d’effets et fonder l’électrodynamique. Ces découvertes expérimentales sont les
signes de l’existence d’un autre champ, appelé champ magnétique, qui se développe en présence d’aimants et/ou
de courants électriques.

1 Expériences historiques : Oersted and Ampère

Avant Oersted, le magnétisme se résumait à constater l’interaction à distance entre des objets aimantés,
par exemple entre la Terre (un gros aimant) et la boussole (un petit aimant), ou bien entre les pôles de petits
aimants. La pile électrique, découverte en 1800 par Volta, est d’abord utilisée par les chimistes pour réaliser
des électrolyses, mais va bientôt devenir l’outil indispensable pour étudier l’électricité de manière approfondie.
En effet, avant la pile de Volta, les physiciens devaient charger des condensateurs (bouteilles de Leyde) en
générant de l’électricité statique (des charges) par frottement. Par conséquent, le courant électrique ne pouvait
être observé que de manière intermittente lors de la décharge de ces condensateurs. Cela ne laissait pas beaucoup
de temps pour observer les effects de ce courant électrique. Grâce à l’amélioration des piles depuis le prototype
de Volta, il devient possible d’établir un courant permanent dans un circuit électrique et d’en observer les effets.
C’est ainsi qu’Hans-Christian Oersted (1777-1851), un professeur danois, constate que le courant électrique
peut produire de la lumière quand il passe dans un fil suffisamment fin. Oersted touche un peu à tout, chimie,
acoustique, optique, électricité, et il est convaincu d’une sorte d’unité des phénomènes physiques. C’est dans le
cadre de cette conception philosophique, qu’il commence à rechercher des effets physiques reliant magnétisme et
électricité.

1.1 Expérience d’Oersted

En 1820, Oersted parvient à démontrer expérimentalement qu’un courant électrique exerce une action à
distance sur une boussole (aiguille aimantée) placée à proximité (Fig. 4.1). Son expérience montre qu’une
aiguille aimantée (montée sur un support qui lui laisse la liberté de ses mouvements) va bouger de manière
à s’orienter perpendiculairement au fil électrique (rectiligne) qui porte le courant. Il observe aussi que l’effet
mécanique change de signe (de direction) si on inverse le sens du courant, comme le montre les figures Fig.
4.1) a) et b). Oersted remarque au passage que ces nouvelles forces, créées par le courant, semblent s’exercer
perpendiculairement au courant.
Cette découverte est fondamentale car elle établit pour la première fois une connexion entre électricité
et magnétisme. Elle est très vite reproduite par Arago devant l’Académie des Sciences à Paris, et devant
Marie-André Ampère.

1.2 Expérience(s) d’Ampère (1820)

Marie-André Ampère (1775-1836) est un mathématicien, physicien, chimiste et philosophe français. Il est
fasciné par la découverte d’Oersted, et se lance dans l’étude systématique des interactions magnétiques entre
courants électriques et aimants, fondant ainsi l’électrodynamique, dans une série de 61 mémoires publiés entre
1820 et 1831. La démarche d’Ampère est fondamentalement expérimentale, mais il réussit à extraire de ses

28
Figure 4.1 – Expérience d’Oersted : Interaction un circuit et une boussole placée au-dessous ce celui-ci
(sans contact). Initialement il n’y a pas de courant dans le circuit, et la boussole est alignée selon l’axe du fil qui
la surplombe (sans contact). On s’arrange aussi pour éliminer l’effet du champ magnétique terrestre. Lorsqu’on
fait passer du courant (I1 6= 0) dans le circuit, la boussole bouge et vient s’orienter perpendiculairement au fil. Si
on inverse la polarité du branchement de la pile (cas b), donc le signe du courant, on observe que la boussole se
place perpendiculairement au fil mais dans le sens opposé à celui du cas précédent a).

expériences des lois mathématiques précises. Ampère sera surnommé par Maxwell "le Newton de l’électricité".
Le physicien James Clerk Maxwell utilisera les lois d’Ampère (et quelques autres dues à Faraday) pour écrire les
équations du champ électromagnétique (années 1860).
Ampère a d’abord l’idée naturelle de remplacer la boussole de l’expérience d’Oersted par un second circuit
électrique (Fig. 4.2). Il réalise l’expérience consistant à tester les actions mécaniques mutuelles entre deux fils
parallèles parcourus par des courants. Il réussit ainsi à démontrer que deux circuits électriques indépendants
(chaque circuit est alimenté par sa propre pile et n’est pas relié à l’autre) exercent l’un sur l’autre des forces
d’interaction à distance. Ces forces n’existent que si les deux conducteurs sont simultanément parcourus par
des courants.
Citons deux conclusions, parmi de nombreuses, extraites du mémoire d’Ampère "De l’action exercée sur un
courant électrique par un autre courant, le globe terrestre ou un aimant (1820)" :
• "Deux courants électriques, s’attirent quand ils se meuvent parallèlement dans le même
sens, ils se repoussent quand ils se meuvent parallèlement en sens contraire".
• "Ces attractions et répulsions sont absolument différentes des attractions et répulsions
électriques ordinaires".
Nous n’allons pas suivre le chemin historique car il est sinueux et il serait plus difficile à exposer que de
présenter le champ magnétique à la lumière des découvertes plus récentes. Cependant on ne peut s’empêcher
de saluer le travail prodigieux d’Ampère qui, à partir d’observations expérimentales, a réussi à écrire les lois
théoriques des interactions magnétiques et fonder l’électrodynamique. Rappelons nous que les physiciens de
cette époque ignoraient la structure atomique de la matière, l’existence des électrons, etc... Ampère a identifié
la notion d’intensité du courant électrique et sa distinction d’avec la tension électrique ; il a formulé les lois
d’attraction répulsion entre les circuits, il a compris que les lois d’interaction entre objets aimantés pouvait être
formulées dans les mêmes termes à condition de "remplacer" un objet magnétique par un système de courants
internes équivalents ; il a inventé le solénoïde (ou bobine) ; il a co-inventé avec Arago l’électro-aimant. Il pensait
qu’un corps aimanté est constitué de courants internes, "les courants ampériens", ce qui n’est pas très loin de la
compréhension moderne du magnétisme.
Pour observer une force mécanique entre les 2 circuits, il faut que I1 6= 0 et I2 6= 0. Avec nos connaissances
sur l’origine microscopique du courant électrique, nous savons que cela implique que des porteurs de charge
libres doivent être mobiles dans les deux circuits. En termes de champ d’action à distance, cela signifie que
les charges mobiles d’un conducteur 1 engendrent un champ (appelé champ magnétique et noté B(~ ~ r)), et que ce
champ n’agit que sur des charges mobiles dans l’autre conducteur 2.
On voit que ce champ magnétique doit être très différent du champ électrique. En effet le champ électrique
peut être créé par des charges immobiles et peut aussi agir aussi sur des charges immobiles. De plus dans toute

29
Figure 4.2 – Expérience d’Ampère : Interaction de deux tronçons de circuits parallèles. Chaque circuit est
indépendant et alimenté par sa propre pile (découverte par Volta en 1800) : a) Intensités de même sens : on
observe une attraction (à gauche) b) Intensités opposées : on observe une répulsion (à droite). Les forces sont
proportionnelles aux intensités I1 et I2 , et décroissent avec la distance d.

l’expérience, les fils restent globalement et localement électriquement neutre puisque de partout le gaz d’électrons
libres de conduction est électriquement compensé par le fond d’ions positifs.
Afin de pouvoir interpréter, cette expérience il nous faut comprendre d’une part comment le conducteur
→
− →
−
1 créé un champ magnétique B 1 (~r) qui baigne tout l’espace, et d’autre part comment ce champ B 1 (~r) agit
localement au niveau du circuit 2 en créant une force (on l’appellera la force de Laplace). Pour cela, dans le
paragraphe suivant, nous allons tout d’abord définir le champ magnétique à partir de la force de Lorentz (1880),
qui est la force exercée sur un champ magnétique sur une charge dans le vide. Puis, nous présenterons la loi
de Biot et Savart qui permet de calculer le champ magnétique créé par un courant électrique donné. En fin de
chapitre, nous présenterons les forces de Laplace exercées sur un circuit par un champ magnétique. Le mécanisme
microscopique (transfert des forces de Lorentz sur les électrons mobiles à une force mécanique sur le réseau
cristallin) à l’origine des forces de Laplace est relié à l’effet Hall classique et sera décrit ultérieurement.

2 Force de Lorentz comme définition du champ magnétique

On sait qu’un aimant ou une bobine (plus généralement un circuit parcouru par des courants électriques)
sont sources d’un champ magnétique qui baigne tout l’espace. Ce champ magnétique va se manifester par son
action mécanique sur une charge "test" q de vitesse ~v (t), dans le référentiel du laboratoire. Comme tout champ
→
−
d’interaction, le champ magnétique B (~r) va même être défini par son action sur une charge "test".
Pour que la description en terme de champ ait tout son sens, il faut que cette force ne dépende que de la valeur
du champ sans qu’il soit nécessaire de préciser quelles sont les dispositifs et sources qui ont engendré ce champ :
l’expression de la force ne doit contenir que le champ lui-même et certaines caractéristiques de la charge test.
Une telle expression n’a été explicitée que très tardivement, en 1885, par le physicien théoricien Hendrik-Antoon
Lorentz (1853-1928).

2.1 Force de Lorentz

On se place dans le référentiel du laboratoire (supposé Galiléen) dans lequel règne un champ magnétique
→
−
statique B (~r) (aussi appelé champ magnétostatique), c’est à dire indépendant du temps. Une charge "test" q en
mouvement, dont la position est donnée par ~r(t) et la vitesse instantanée par ~v (t), subit la force dite de Lorentz :
→
− →
−
F L (t) = q~v (t) ∧ B (~r(t)) . (4.1)

Cette force est nulle pour une charge immobile. La force est également nulle quand la vitesse est alignée avec la
direction du champ magnétique (le produit vectoriel de deux vecteurs colinéaires est nul). Plus généralement,
~ r) et la vitesse, quand ces deux derniers
cette force est toujours perpendiculaire au plan formé par le champ B(~

30
vecteurs ne sont pas colinéaires. Etant perpendiculaire à la vitesse, la force de Lorentz ne travaille pas et ne peut
donc pas changer l’énergie cinétique de la particule (théorème de l’énergie cinétique). En revanche cette force de
Lorentz magnétique peut courber les trajectoires.
Et si de plus le champ est uniforme (constant par rapport à la variable d’espace), l’expression est encore plus
simple :
→
− →
−
F L (t) = q~v (t) ∧ B . (4.2)
Ici on ne démontre pas l’expression de la force Eq. (4.1) : elle nous sert de définition du champ magnétique et
de nombreuses expériences ont depuis montré sa validité. Notons à nouveau que cette expression de l’interaction
entre une charge et un champ magnétique n’a été écrite pour la première fois qu’en 1885 par Hendrik-Antoon
Lorentz (soit un siècle après la loi de Coulomb). Ceci est du au fait que les particules élémentaires n’ont commencé
à être identifiées qu’avec la théorie atomique et la découverte de l’électron vers la fin du XIXième siècle.

2.2 Unités et ordres de grandeur

L’unité de champ magnétique en (SI) est fixée par la relation Eq. (4.1) puisque q est en coulombs (C), la
vitesse en m.s−1 , et la force en newtons (N). Cette unité (SI) de champ magnétique est appelée le tesla (T).
Nikola Tesla était un physicien et ingénieur américain (d’origine serbe) de la fin du XIXième siècle et début
du XXième qui a notamment travaillé sur la distribution de l’alimentation électrique, en rivalité avec Thomas
Edison.
On les ordres de grandeur suivants pour les champs magnétiques en statique :
• Champ magnétique terrestre en France : 2.10−5 T pour la composante horizontale, et 4.10−5 T
pour la composante verticale.
• Champ dans l’entrefer d’un électro-aimant : 100 mT à 2 T.
• Champ dans une bobine supraconductrice : 5 T à 45 T.

2.3 Considérations de symétrie par rapport à un plan : cas général

Considérons une situation physique réelle, notée 1, décrite par la dynamique suivante :
d~v1 (t) ~ 1 (~r1 (t)) ,
m = q~v1 (t) ∧ B (4.3)
dt
et correspondant à une certaine trajectoire ~r1 (t) bien réelle de la particule.
Maintenant imaginons la même particule (charge q et masse m) qui décrirait la trajectoire obtenue en
appliquant une symétrie miroir par rapport à un plan Π à la trajectoire 1. La trajectoire ainsi obtenue correspond
à une situation hypothétique que nous allons noter 2. A chaque instant, par construction géométrique, la position
sur cette trajectoire 2 serait donnée par :

~r2 (t) = SymΠ [~r1 (t)] , (4.4)

et donc la vitesse et l’accélération se transformeraient aussi comme :

 
d~v1 (t) d~v2 (t) d~v1 (t)
~v1 (t) → ~v2 (t) = SymΠ [~v1 (t)], et → = SymΠ . (4.5)
dt dt dt
~ 2 (~r) qui conduirait
On peut se poser la question suivante : existe-il une cartographie de champ magnétique B
exactement à la cinématique (trajectoire 2) hypothétique construite par symétrie miroir de la situation réelle ?
Si un tel champ existe, il doit satisfaire à tout instant à l’équation de Newton :
d~v2 (t) →
−
m = q~v2 (t) ∧ B 2 (~r2 (t)) . (4.6)
dt
En comparant cette équation avec l’équation initiale Eq. (4.3), on voit qu’il faut qu’à chaque instant :
→
− →
− →
−
B 1 (~r1 (t)) → B 2 (~r2 (t)) = −SymΠ [ B 1 (~r1 (t))] , (4.7)

le signe moins étant du à la présence du produit vectoriel.

→
−
Et si l’on veut que cela soit vrai pour toute condition initiale, il faut que le champ statique B 2 (~r) satisfasse :
→
− →
−
B 2 (~r2 ) = −SymΠ [ B 1 (~r1 )] . (4.8)

31
On voit que, sous une opération de symétrie miroir, le vecteur champ magnétique se transforme différemment
des vecteurs position, vitesse ou accélération. On dit que ces derniers sont des vecteurs polaires, tandis que
le champ magnétique est un vecteur axial. C’est une conclusion très générale et valable quel que soit le
champ magnétique. Nous reviendrons sur ce raisonnement dans un cas particulier plus concret au paragraphe
3.2.
Vous connaissez d’autres vecteurs axiaux en physique : le moment d’une force, le moment cinétique, ... Le
produit vectoriel de deux vecteurs polaires est un vecteur axial, par exemple le moment cinétique est le produit
vectoriel de deux vecteurs polaires : le vecteur position et la quantité de mouvement. On verra que le champ
magnétique s’exprime aussi comme une somme de produits vectoriels via la loi de Biot et Savart dans la section
4. En effet nous aurions pu déduire directement de la forme de l’équation de la dynamique que B ~ est axial, car
force, vitesse et champ sont reliées par un produit vectoriel. La force et la vitesse étant des vecteurs polaires, le
troisième partenaire (ici le champ magnétique) doit être axial.

Conséquence : Imaginons un champ magnétique B(~ ~ r) engendré par des courants sources qui présentent
un plan de symétrie (Π). Cela signifie que les sources sont invariantes par la symétrie de plan (Π), et donc les
champs créés dans les situations 1 et 2 sont identiques ainsi :
→
− →
− →
−
B (~r1 ) → B (~r2 = SymΠ [~r1 ]) = −SymΠ [ B (~r1 )] , (4.9)

Et plus simplement :
→
− →
− →
−
B (~r) → B (SymΠ [~r]) = −SymΠ [ B (~r)] , (4.10)
Il s’agit d’une contrainte sur les lignes de champ, une relation entre le vecteur champ magnétique en un point ~r
d’une part, et en son symétrique SymΠ [~r] d’autre part. Cette contrainte devient encore plus forte lorsque le
point ~r est contenu dans le plan de symétrie car alors ~r = SymΠ [~r], et la relation ci-dessus devient une relation
entre le vecteur champ en ce point et lui même :
→
− →
− →
−
B (~r) → B (~r) = −SymΠ [ B (~r)] , (4.11)
→
−
Pour satisfaire cette dernière contrainte, il faut que le champ magnétique B (~r) soit perpendiculaire au plan de
symétrie Π.
Propriété du champ magnétique dans un plan de symétrie des courants : Soit ΠS un plan de
→
−
symétrie des courants sources. Le champ créé B (M ) en M ∈ ΠS est orthogonal au plan ΠS :

→
−
∀ M ∈ ΠS B (M ) ⊥ ΠS . (4.12)

3 Mouvement d’une particule chargée dans un champ magnétique uniforme

On va s’intéresser à l’action mécanique d’un champ magnétique existant sur le mouvement d’une charge
placée dans le vide. La particule sera soumise directement à la seule force de Lorentz. Pour simplifier, nous
→
−
étudions l’action d’un champ magnétique uniforme et constant B sur le mouvement d’une particule de charge
q dans le vide.

3.1 Mouvement cyclotron dans un champ magnétique uniforme

→
−
On prend l’axe z, selon la direction du champ magnétique de sorte que : B = B0 ~uz . On prend les conditions
initiales suivantes : ~r(0) = ~0 et ~v (0) = v0 ~ux , comme dans la Fig. 4.3.
L’équation du Newton est donc :
d~v (t) →
−
m = q~v (t) ∧ B . (4.13)
dt
• En projection selon l’axe Oz, cela donne : dvz /dt = 0 donc vz (t) = vz (0) = 0.
• Dans le plan Oxy la situation est plus compliquée, car les équations :

dvx (t) qB0 dvy (t) qB0

= vy (t), et =− vx (t) . (4.14)
dt m dt m
sont des équations différentielles couplées.

32
Figure 4.3 – Mouvement cyclotron d’un électron : Une charge q, négative sur ce dessin, se déplace dans
un champ magnétique uniforme, sous l’influence de la seule force de Lorentz. Notez que la force de Lorentz dépend
du temps via la vitesse. La trajectoire est un cercle de rayon Rc car la vitesse initiale était dans le plan 0xy. Si
la vitesse initiale contenait une composante non nulle suivant l’axe du champ (z), alors un mouvement rectiligne
uniforme suivant z se superposerait au mouvement cyclotron circulaire dans xy, et la trajectoire deviendrait une
hélice s’enroulant autour de la direction (0z) du champ magnétique.

Il y a plusieurs méthodes pour découpler ces équations et trouver la solution. L’une d’elles consiste à
considérer le nombre complexe w(t) = vx (t) + ivy (t), dont les partie réelles et imaginaires sont respectivement
les composantes vx (t) et vy (t) de la vitesse instantanée. Graphiquement, dans le plan complexe, le point d’affixe
w(t) représente donc "la pointe" du vecteur vitesse. On peut alors exprimer les deux équations réelles ci-dessus
en une seule équation complexe :
dw(t) dvx (t) dvy (t) qB0 iqB0
= +i = (vy (t) − ivx (t)) = (vx (t) + ivy (t)) = iw(t) . (4.15)
dt dt dt m m
Cette équation différentielle en w(t) est linéaire du premier ordre et s’intègre très facilement en une exponentielle :
 
iqB0
w(t) = w(0) exp t . (4.16)
m

Les conditions initiales (à t = 0) sur la vitesse sont contenues dans w(0) = vx (0) + ivy (0) = v0 . On a donc (en
utilisant la formule d’Euler pour l’exponentielle d’un nombre imaginaire pur) :
   
qB0 qB0
vx (t) = v0 cos t , et vy (t) = v0 sin t . (4.17)
m m

On voit que le vecteur vitesse tourne (on dit parfois effectue un mouvement de précession, ou "précesse") dans
le plan xy, et donc autour de l’axe de rotation 0z, à la vitesse angulaire :
qB0
ωc = , (4.18)
m
appelée pulsation cyclotron. Cette rotation se fait dans le sens trigonométrique si q > 0, et dans le sens
horaire si q < 0.
Pour obtenir la trajectoire, il suffit de calculer les primitives des composantes de la vitesse tout en tenant
compte des conditions initiales. On trouve :
   
mv0 qB0 mv0 qB0
x(t) = sin t , et y(t) = (1 − cos t ). (4.19)
qB0 m qB0 m

Ainsi la trajectoire est un cercle de rayon :

mv0
Rc = . (4.20)
qB0

33
Figure 4.4 – A gauche (situation 1) : Un électron effectue un mouvement cyclotron dans un champ magnétique
uniforme. La trajectoire est parcourue dans le sens trigonométrique. On peut en déduire que le champ magnétique
est orienté dans le sens de (+~uz ) A droite (situation 2) : L’ image de l’électron par la symétrie miroir de plan
(Π) effectue aussi un mouvement circulaire uniforme mais dans le sens anti-trigonométrique. Un observateur en
déduira que l’image de l’électron évolue dans un champ magnétique (non représenté) orienté dans le sens de
−~uz .

parcouru à une vitesse angulaire constante (Fig. 4.3). Ce rayon, appelé "rayon cyclotron", augmente quand la
quantité de mouvement mv0 de la particule augmente, et diminue si le champ magnétique B0 augmente. Dans
un spectromètre de masse, on envoie des particules chargées de vitesse commune v0 connue dans un champ
connu B0 . La mesure du rayon de giration Rc permet d’accéder à la valeur du rapport q/m de la particule.

3.2 Considérations de symétrie : un exemple sur le mouvement cyclotron

Nous revenons ici sur les considérations générales du paragraphe 2.3 selon lesquelles le champ magnétique
est un vecteur axial. Nous allons reproduire le raisonnement mais en utilisant le cas particulier du mouvement
cyclotron d’un électron dans un champ magnétique uniforme.
• 1) Imaginons qu’un cameraman ait filmé le mouvement cyclotron de l’électron représenté sur la Fig. 4.4).
Il s’agit d’un enregistrement d’une situation bien réelle et réalisée dans la Nature (situation 1 évoquée
dans le paragraphe). Le cameraman montre ce film à un physicien qui sait seulement qu’il s’agit d’un
électron. En visualisant ce film, notre physicien peut conclure que l’électron est soumis à une force de
Lorentz correspondant à un champ magnétique orienté dans le sens +~uz .
• 2) Le cameraman place maintenant un grand miroir plan dans le plan (Π) de la Figure Fig. 4.4), puis il
filme l’image de l’électron dans le miroir. Il obtient un film montrant un électron qui tourne dans le sens
opposé (sens horaire dans le plan xy) à celui du "vrai" électron. Ce film correspond la situation 2 du
paragraphe. On peut se demander s’il reflète un mouvement autorisé par les lois de la Nature, après tout
il s’agit d’un montage, "d’un fake" où le cameraman à utilisé un miroir et ne filme qu’une image...
• 3) Le cameraman montre ce second film au physicien qui ne sait pas qu’il s’agit d’un "fake", d’un montage.
En fait, le physicien n’a aucun moyen de dire qu’il s’agit d’un montage. Il conclut que l’électron se déplace
dans un vrai champ magnétique orienté selon les −~uz .

Vous pouvez reprendre le raisonnement de la Fig. 4.4), en partant du même mouvement initial (situation 1)
mais en utilisant d’autres plans miroirs Π, par exemple des plans parallèles à xy ou yz. On obtient toujours la
conclusion suivante :

Conclusion : le champ magnétique dans la situation 2 (transformée par la symétrie miroir) doit être l’opposé
du symétrique du champ magnétique de la situation 1 originale. Pour exprimer cela, on dit que le champ
magnétique est un vecteur axial.

34
4 Champ créé par des courants sources : loi de Biot et Savart (1820)
Jusqu’à présent, nous avons seulement étudié l’action d’un champ magnétique donné sur une charge test
mobile. Nous ne nous sommes pas encore posés la question de l’origine de ce champ magnétique et de sa relation
avec les courants sources qui l’on engendré. En réalisant des expériences en 1820, les physiciens Jean-Baptiste
Biot et Félix Savart ont réussi à exprimer cette relation qui se traduit sous sa forme moderne par une formule
→
−
donnant le champ B (M ) créé en un point M arbitraire en fonction de la distribution de courants sources ~j(P ).

4.1 Enoncé pour des courants sources volumiques

→
−
Le champ magnétique B (M ) créé en M par une distribution volumique de courants étendue ~j(P ) s’écrit :

→
− µ0
ZZZ ~j(P ) ∧ ~u
B (M ) = 2
dτ, (4.21)
4π P ∈V P M

où ~u = ~uP M est le vecteur unitaire pointant de P vers M . On appelle "élément de courant volumique
autour de P " le vecteur infinitésimal ~j(P )dτ autour de P . Par ailleurs la constante de couplage dans le système
international d’unité (SI), s’écrit :
µ0
= 10−7 (SI), (4.22)
4π
où µ0 s’appelle la perméabilité magnétique du vide. En fait son unité exacte est le henry par mètre (H.m−1 ), le
henry étant l’unité d’inductance.
On voit que la contribution provenant d’un élément de courant placé en P , et observée en M , décroît en
1/P M 2 , comme pour le champ newtonien de la gravitation et le champ coulombien électrostatique. En revanche,
les sources sont ici vectorielles et non scalaires, ce qui se manifeste par un produit vectoriel entre le vecteur
courant source et le vecteur unitaire ~uP M . On peut constater l’analogie formelle avec la formule donnant le
champ électrique en fonction de la densité volumique de charge sources ρ(P ) :

→
− →
− 1 µ0 ~j(P )dτ ∧ ~u ↔ ρ(P )dτ ~u,
E (M ) ↔ B (M ), ↔ et (4.23)
4π0 4π

Cette loi exprime aussi un principe de superposition : le champ crée par la somme de distributions de courants
est le somme des champs magnétiques créés par chacune des distributions.

4.2 Enoncé pour des circuits filiformes

Les circuits électriques sont souvent constitués de fils qui sont des tubes de dimensions transverses faibles
par rapport à leur longueur, d’où la possibilité de l’approximation des circuits filiformes qui consiste à remplacer
ces tubes par des lignes de section nulle portant la même intensité I que les fils réels. La correspondance se fait
~ car en effet l’intensité s’écrit I = jS, l’élément de volume dτ = Sdl, et les vecteurs ~j
en écrivant ~j(P )dτ = I dl,
~
et dl sont colinéaires de même sens. On a donc la loi de Biot-Savart pour les circuits filiformes :

→
− µ0
Z ~ ∧ ~u
I dl µ0
Z ~ ∧−
I dl
−→
PM
B (M ) = 2
= , (4.24)
4π P PM 4π P P M 3
−−→
où dans la dernière égalité on a a exprimé ~u = ~uP M = P M /P M le vecteur unitaire pointant de P vers M .
~ = dl
Le vecteur dl ~ P est localisé en un point P du circuit et dirigé le long du circuit dans le sens des I positif.
~
On appelle "élément de courant linéique autour de P " le vecteur infinitésimal I dl.
Que ce soit sous la forme "courants volumiques" Eq. (4.21) ou "circuits filiformes" Eq. (4.21), la formule de
Biot et Savart exprime que le champ magnétique est un produit vectoriel de vecteurs polaires : c’est donc un
vecteur axial (comme nous l’avions anticipé par une analyse de la dynamique dans le paragraphe 2.3).

4.3 Cas d’un fil rectiligne infini

On peut appliquer la formule de Biot-Savart au cas d’un fil rectiligne très long, disons infini pour simplifier
les calculs, et parcouru par un courant d’intensité I. On choisit le système de coordonnées cylindriques M (r, θ, z)

35
Figure 4.5 – Fil rectiligne infini parcouru par un courant : Les lignes de champ s’enroulent autour de la
direction du fil. On a représenté la même situation dans le plan ΠS = (M~ur ~uz ) qui est un plan de symétrie, et
dans le plan ΠA = (M~ur ~uθ ) qui est plan d’antisymétrie (vue de dessus).

d’axe 0z matérialisé par le fil source. On note ~ur , ~uθ , ~uz la base locale en M . On a :

~ = Idz ~uz et −
I dl
−→ ~ ∧−
P M = r~ur − z~uz , donc I dl
−→
P M = Irdz~uθ . (4.25)

D’où :
→
− µ0
Z ~ ∧−
I dl
−→
PM µ0 I
Z ∞
rdz

B (M ) = = ~uθ , (4.26)
4π P ∈V PM3 4π 2 2 3/2
−∞ (r + z )

Comme on intègre sur z (la cote de P ), et que r est fixe (position du point d’observation), on peut faire apparaître
la variable sans dimension u = z/r :
Z ∞  Z ∞ 
→
− µ0 I d(z/r) µ0 I du
B (M ) = 2 3/2
~
u θ = 2 3/2
~uθ , (4.27)
4πr −∞ (1 + (z/r) ) 4πr −∞ (1 + u )

Après intégration on trouve (l’intégrale entre paranthèses vaut 2) :

→
− µ0 I
B (M ) = ~uθ , (4.28)
2πr
Comme on voit sur la Fig. 4.5, le champ magnétique est donc orthoradial. Les lignes de champ magnétique
sont donc des cercles qui "s’enroulent" autour de l’axe du courant. L’intensité du champ décroit comme 1/r. Si
on change le signe de I, tous les vecteurs champs changent de sens. Bien sûr si I = 0, le champ est identiquement
nul de partout.
Remarque de symétrie : Ici on a fait un calcul direct et on a trouvé le champ suivant ~uθ . On remarque
aussi que le plan ΠS = (M~ur ~uz ) contenant le fil et le point M d’observation est un plan de symétrie du courant.
En effet, le courant est contenu dans ce plan, et donc inchangé par la symétrie miroir associée. Donc le champ
→
−
magnétique étant un vecteur axial, B (M ) doit être égal à l’opposé de son image par ce plan miroir, conformément
→
−
à la propriété (5.13). Cela n’est possible que si B (M ) est orthogonal à ΠS , donc effectivement orthoradial.

36
5 Action d’un champ magnétique sur un élément de circuit électrique
5.1 Force de Laplace
Afin d’interpréter son expérience Fig. 4.2, Ampère trouve une formule pour l’action d’un champ magnétique
~ :
~ sur un élément infinitésimal de circuit électrique I dl
B
−→ ~ ∧ B.
~
dF = I dl (4.29)

De son côté Biot, aidé de Savart, avait réussi à établir une expression intégrale pour cette force. Il y aura donc
un débat sur l’attribution de cette loi à Biot-Savart ou à Ampère. On demande à Laplace (qui a alors 71 ans) de
trancher la question en arbitre impartial. Finalement à la demande de Laplace, Ampère parvient à montrer que
sa loi infinitésimale conduit après intégration à la loi globale de Biot-Savart. Ironie de la postérité, on appelle
maintenant force de Laplace la formule Eq. (4.29) qui devrait plutôt s’appeler Biot-Ampère !

5.2 Interprétation de l’expérience d’Ampère sur les forces entre fils parallèles
Reprenons la situation de l’expérience d’Ampère décrite par la figure Fig. (4.2).
Action du circuit 1 sur le circuit 2 : Pour simplifier les calculs, on assimile le circuit 1 au tronçon de fil
1 supposé très long, et on néglige le champ créé par les autres parties du circuit 1 qui servent à boucler le circuit.
On peut donc estimer que le circuit 1 créé un champ magnétique :
→
− µ0 I 1
B 1 (M ) = ~uθ1 , (4.30)
2πr1
où r1 et ~uθ1 correspondent aux coordonnées cylindriques d’axe le fil 1.
~ 2 du circuit 2 :
Ce champ magnétique exerce une force dF~12 sur un élément de circuit I2 dl
−→ ~ 1 = − µ0 I1 I2 dl2 ~ur1 ,
~2∧B
dF 12 = I2 dl (4.31)
2πd
car r1 = d au niveau du fil 2.
On voit que si I1 I2 > 0, alors le fil 1 attire le fil 2 (cas attractif). En revanche si I1 I2 < 0, alors le fil 1
repousse le fil 2 (cas répulsif).
De plus, cette formule a longtemps servi de définition pour l’unité d’intensité : 1 A est l’intensité qui doit
circuler dans deux fils parallèles séparés de 1 m pour que la force soit de 2.10−7 N par mètre de fil.

Action du circuit 2 sur le circuit 1 : On peut faire le raisonnement réciproque dans lequel le circuit 2
agit à distance sur le circuit 1. Le circuit 2 créé un champ :
→
− µ0 I 2
B 2 (M ) = ~uθ2 , (4.32)
2πr2
où r2 et ~uθ2 correspondent aux coordonnées cylindriques d’axe le fil 1. Ce champ magnétique exerce une force
~ 1 du circuit 1 :
dF~21 sur un élément de circuit I1 dl
−→ ~ 2 = − µ0 I1 I2 dl1 ~ur2 ,
~1∧B
dF 21 = I1 dl (4.33)
2πd
−→
car r2 = d au niveau du fil 2. Sur la figure Fig. (4.6), cela correspond à une force dF 21 vers la droite pour
−→
I1 I2 > 0, donc bien opposée à dF 12 .

6 Exercices

~ r) par Biot-Savart
Exercices sur le calcul direct du champ B(~
1) Calculer le champ magnétique créé par une spire de courant sur son axe.

2) Calculer le champ magnétique créé par un fil rectiligne infini parcouru par un courant I.

37
Figure 4.6 – Force dans l’expérience d’Ampère : La force infinitésimale dF ~ 12 = I2 dl~2∧B ~ 1 est la force
~
subie par l’élément de circuit 2 I2 dl2 de la part de la totalité du circuit 1 (assimilé à 1 fil infini). Le champ B ~1
est le champ magnétique créé par l’ensemble du circuit 1 (assimilé à 1 fil infini).

38
Chapitre 5

Les propriétés générales du champ

magnétique

Nous présentons ici les deux grandes propriétés générales valables pour n’importe quel champ magnétique,
quelles que soient les sources (courants et ou aimants) à l’origine du champ. D’une part, le flux du champ
magnétique est conservé le long des tubes de champ et donc nul à travers toute surface fermée. D’autre part,
la circulation du champ magnétique le long d’une boucle fermée est proportionnelle à l’intensité du courant
électrique enlacée par ce lacet.

1 Champ magnétique à flux conservatif

Commençons ici par énoncer (avant de la démontrer) une propriété fondamentale du champ magnétique :
son flux à travers n’importe quelle surface fermée fictive (Σ) est nul :

→
− −
→
ZZ
B (M ).dS ext = 0 . (5.1)
M ∈Σ

On dit que le champ magnétique est à flux conservatif. C’est une propriété globale car elle met en jeu les
valeurs du champ magnétique en des points recouvrant une surface qui peut être d’extension macroscopique
dans l’espace.
Cette propriété globale du flux équivaut mathématiquement à imposer que la divergence du champ magnétique
est nulle en tout point de l’espace :
→
−
div B (M ) = 0 , (5.2)

qui est une propriété locale car elle n’engage que les valeurs du champ au voisinage de chaque point M .
L’équivalence entre Eq. (5.1) et Eq. (5.2) découle du théorème de la divergence (ou Green-Ostrogradsky),
mais nous reviendrons sur cette démonstration dans le chapitre sur les équations de Maxwell. Cette propriété
de divergence nulle de partout (qui est l’une des 4 équations de Maxwell) équivaut à dire que les monopôles
magnétiques n’existent pas. Le physicien théoricien Paul Dirac a cependant étudié les propriétés qu’auraient de
tels monopôles magnétiques s’ils existaient, mais jusqu’ici aucune expérience n’a pu en détecter.

1.1 Cas particulier du champ créé par un fil rectiligne infini

Dans le cas particulier du champ créé par un fil rectiligne infini, on a obtenu une formule explicite Eq. (4.28).
On peut donc calculer directement le flux en injectant Eq. (4.28) dans Eq. (5.1). Comme le montre la Fig. 4.5, les
lignes de champ sont des cercles d’axe Oz, et les tubes de champ sont donc des tores d’axe Oz. Si on considère,
deux sections S1 et S2 d’un tube de champ, on a :

→
− −→ →
− −
→
ZZ ZZ
B (M ).dS = B (M ).dS , (5.3)
M ∈S1 M ∈S2

−
→
avec dS = dS ~uθ . Cette relation implique :

− −
→ →ext − −
→ → − −
→ → − −
→ →
ZZ ZZ ZZ ZZ
B .dS = − B .dS 1 + B .dS 2 B .dS lat = 0 (5.4)
Σ S1 S2 Slat

39
 −
→ −
→ −
→ − −
→ →
car dS ext = −dS 1 = dS 2 , et il n’y a pas de flux à travers la surface latérale d’un tube de champ B .dS lat = 0.
On a donc montré la relation Eq. (5.1) dans le cas d’un tube de champ. On peut généraliser à une surface
fermée quelconque en regroupant les lignes de champ en plusieurs tubes élémentaires.

1.2 Cas d’un circuit filiforme quelconque

Pour généraliser, il faut partir de la formule de Biot-Savart donnant le champ magnétique au point M :

→
− µ0
Z ~ ∧ ~u
I dl
B (M ) = , (5.5)
4π P r2
où on note r = P M et ~u = ~uP M . On ne détaille pas ici la démonstration, mais l’idée est que chaque champ
→
− ~ présente les mêmes symétries de révolution que celles utilisées
élémentaire d B , créé par l’élément de circuit I dl,
pour un fil rectiligne et génère donc un flux infinitésimal nul. Par sommation de ces flux nuls, on obtient un flux
total nul à travers toute surface fermée.

1.3 Démonstration par le calcul de la divergence

Si on n’est pas convaincu par les raisonnements "avec les mains" et de symétrie précédents, on peut aussi
démontrer purement analytiquement que la divergence du champ magnétique est nulle de partout. Pour cela, il
faut calculer la divergence de la formule de Biot-Savart par rapport à M (on dérive par rapport aux coordonnées
de M) : !
→
− µ0 I
Z ~ ∧ ~u
dl
div B (M ) = divM , (5.6)
4π P r2
On utilise la formule d’analyse vectorielle suivante :

~ =− → ~ ~ −→ ~ ~
 
div U ~ ∧V rotU . V − rotV . U , (5.7)

appliquée à U ~ et V
~ = dl ~ = ~u/r2 :
!
~ ∧ ~u
dl −→ ~ ~ −→

~u

divM = rotM (dl) . V − rotM ~ = 0,
. dl (5.8)
r2 r2

car :  
−→ ~ −→ ~u
rotM (dl) = ~0 et rotM = ~0 . (5.9)
r2
~ ne dépend pas des coordonnées de M , et la seconde se vérifie par le
La première égalité découle du fait que dl
calcul direct (formulaire ou internet pour avoir le rotationnel en coordonnées sphériques et on l’applique au cas
d’un champ radial en 1/r2 ).

2 Circulation du champ magnétique et théorème d’Ampère

La circulation du champ magnétique sur une boucle fermée (fictive) C est aussi très particulière car elle
dépend des courants enlacés par le contour fermé C selon la formule :

→
− ~
I
B · dl = µ0 Ienl , (5.10)
P ∈C

où Ienl désigne le courant algébrique enlacé par le contour C. Cette relation est connue sous le nom de théorème
d’Ampère et s’avère très utile pour calculer le champ magnétique dans des situations très symétriques.
On peut vérifier si cette égalité Eq. (5.10) fonctionne pour le cas du champ créé par un fil rectiligne infini,
que nous avons déjà déterminé par la loi de Biot et Savart au chapitre précédent, cf. Eq. (4.28). On allons donc
vérifier si le théorème d’Ampère est est valide dans ce cas particulier. Prenons tout d’abord pour boucle fictive
d’Ampère C un chemin quelconque qui entoure le fil (pas forcément un cercle). L’élément de boucle est noté
~ = dr~ur + rdθ~uθ + dz~uz . Alors la circulation du champ magnétique s’écrit :
dl
→
− ~
I I I
µ0 I dθ
B · dl = ~uθ .rdθ~uθ = µ0 I = µ0 I . (5.11)
P ∈C P ∈C 2πr P ∈C 2π

40
Figure 5.1 – Théorème d’Ampère : a) La boucle d’intégration (C) enlace le fil électrique portant une intensité
I. b) La boucle d’intégration (C) n’enlace pas le courant électrique.

Notez la simplification r/r = 1 dans l’intégrand, qui autorise à ce que r ne soit pas forcément constant et donc
le lacet peut bien avoir une forme quelconque (pas forcément un cercle).

Que vaut l’intégrale de l’angle θ sur un chemin fermé ?

I
dθ
= ?, (5.12)
P ∈C 2π

En fait la réponse n’est pas triviale à cause du fait que l’angle est défini modulo 2π. Cette intégrale est en fait
une quantité topologique qui mesure combien de fois la boucle s’entoure autour de l’axe des z. D’après la Fig.
(5.2), on voit que cette intégrale vaut +1 si la boucle enlace l’axe Oz, et vaut zéro si elle ne l’enlace pas.

3 Les symétries du champ magnétique

En étudiant l’invariance de l’équation de Newton avec une force de Lorentz, nous avons obtenu que le champ
magnétique est un vecteur axial. Et nous avons pu en déduire la propriété suivante :

Propriété du champ magnétique dans un plan de symétrie des courants : Soit ΠS un plan de
→
−
symétrie des courants sources. Le champ créé B (M ) en M ∈ ΠS est orthogonal au plan ΠS :

→
−
∀ M ∈ ΠS B (M ) ⊥ ΠS . (5.13)

Cette propriété peut aussi être démontrée à partir de la formule de Biot et Savart en additionnant deux à deux
les contributions d’éléments de circuits élémentaires symétriques.

Pour un plan d’antisymétrie (ΠA ), on obtient la propriété suivante :

Propriété du champ magnétique dans un plan d’antisymétrie des courants : Soit ΠA un plan de
→
−
symétrie des courants sources. Le champ créé B (M ) en M ∈ ΠA est contenu dans le plan ΠA :

→
−
∀ M ∈ ΠA B (M ) ∈ ΠA . (5.14)

41
Figure 5.2 – Nombre d’enroulement d’un lacet autour d’un point O (notion mathématique) : a)
Le lacet entoure le point O. Lorsqu’on fait le tour du lacet C, l’angle polaire θ augmente de manière monotone de
0 à 2π, l’intégrale sur les dθ résultante vaut 2π. b) Le lacet n’enroule pas le point origine O. L’angle θ commence
par augmenter de θ1 à θ2 avec des dθ > 0, puis θ décroît de θ2 vers θ1 , avec des dθ négatifs qui compensent les
dθ positifs de la première partie du chemin. Cette notion est une caractéristique topologique du lacet et du point
O : sa valeur ne dépend pas de la forme précise de la boucle mais seulement du fait que celle-ci entoure ou pas
l’origine.

~ r)} à partir des symétries des

~ r), (B(~
Méthode pour obtenir des informations sur les champs {E(~
sources
1) Utilisation des plans de symétries pour déterminer la direction des lignes de champ. Pour obtenir
la direction d’un champ en un point M de l’espace, il convient d’utiliser un ou des plan(s) de symétrie
ou d’antisymétrie des distributions sources passant par M .

2) On peut aussi utiliser d’autres plans de symétrie/antisymétrie (ne passant pas M ) pour restreindre
l’étude à un demi-espace par exemple, mais cela ne suffira pas pour obtenir la direction du champ en
tout point.

3) Utilisation des invariances pour déterminer le comportement des composantes non nulles des
champs (ou de la norme des champs) en fonction des différentes coordonnées du point M . Ce travail va
de pair avec le choix d’un système de coordonnées adapté.

4 Exemple d’application du théorème d’Ampère

Reprenons notre exemple du fil rectiligne infini le long de l’axe Oz, et parcouru par un courant I. Nous
avons déjà calculé le champ magnétique créé par ce type de courant en utilisant la formule de Biot et Savart, cf.
Eq. (4.28). Supposons un instant que nous ne connaissions pas ce résultat Eq. (4.28). Nous allons montrer que
le théorème d’Ampère et l’analyse des symétries permettent de calculer très facilement ce champ sans avoir à
calculer d’intégrales.
La méthode est la suivante et consiste en deux étapes :

1) Analyse des symétries pour obtenir la direction du champ magnétique : Le plan ΠS = (OM z),
qui contient le fil et le point M d’observation, est un plan de symétrie pour les courants puisque le fil est invariant

42
 →
−
par la symétrie plane correspondante. Par conséquent, comme le champ magnétique est un vecteur axial, B (M )
est orthogonal à ΠS , donc orthoradial :
→
−
B (M ) = Bθ (r, θ, z) ~uθ . (5.15)
De plus cette composante orthoradiale ne peut pas dépendre de z à cause de l’invariance par translation selon
Oz, et ne peut pas dépendre de θ à cause de l’invariance par rotation autour de l’axe Oz. Donc :
→
−
B (M ) = Bθ (r) ~uθ . (5.16)
Ainsi, au lieu d’avoir à déterminer 3 composantes chacune fonction des 3 coordonnées du point M , nous avons
seulement à déterminer une fonction d’une seule variable Bθ (r).
2) Utilisation du théorème d’Ampère : Mais comment calculer la circulation Eq. (5.10) d’une fonction
encore inconnue Bθ (r) ? Trouver la primitive d’une fonction que l’on ne connaît pas, cela semble perdu d’avance...
L’astuce consiste à choisir un contour d’intégration sur lequel cette fonction inconnue ne varie pas, donc un
contour avec r constante. Cela dicte le choix du chemin fermé d’Ampère (C) comme étant un cercle d’axe Oz et
de rayon r (passant par M ). La circulation se calcule alors très simplement, en sortant Bθ (r) de l’intégrale, ce
qui donne :

→
− ~
I I I
B · dl = Bθ (r) ~uθ . rdθ~uθ = rBθ (r) dθ = 2πrBθ (r), (5.17)
C C C
c’est à dire simplement Bθ (r) multiplié par la circonférence du cercle C.
Comme le lacet d’intégration (C) enlace le fil portant une intensité électrique I, la théorème d’Ampère
donne :
2πrBθ (r) = µ0 I, (5.18)

Conclusion : le champ magnétique créé par un fil rectiligne infini parcouru par un courant I est donné
par :
→
− µ0 I
B (M ) = ~uθ . (5.19)
2πr
Ceci est en accord avec le calcul direct utilisant la loi de Biot et Savart.

5 Equations de Maxwell de la magnétostatique

En régime statique, toutes les dérivées temporelles ∂/∂t sont nulles, conduisant à des équations fermées pour
le champ magnétique et les courants sources :
→
− ~ →
− ~
∇ · B(~r) = 0 , ∇ ∧ B(~r) = µ0~j(~r) (5.20)

1) Ces équations relient les dérivées spatiales du champ magnétostatique entre elles et aux courants sources.
Elles sont locales, mais leurs solutions sont non locales.
2) La première équation est homogène : elle exprime que le champ magnétostatique est à flux conservatif
(absence de monopôles magnétiques), i.e. il s’écrit comme le rotationnel d’un potentiel vecteur A(~~ r).
3) La seconde équation relie le champ à ses sources (démontrer le théorème d’Ampère à partir de cette
équation).
→
−
4) La loi de conservation de la charge est vérifiée car ρ = 0 et ∇.~j(~r) = 0.
5) Ecrire l’équation de Poisson pour le potentiel vecteur A(~ ~ r) et la solution correspondante. Ré-exprimer
~
cette solution dans la cas de circuits filiformes en utilisant : ~jdτ ∼ I dl.
6) En déduire la loi de Biot-Savart exprimant le champ magnétique en fonction des courants :

→
− µ0
ZZZ ~j(P )dτ ∧ ~u
B (M ) = , (5.21)
4π P ∈V PM2
où ~u = ~uP M est le vecteur unitaire pointant de P vers M . Pour cela, on utilisera la formule d’analyse
vectorielle :
!
~j(P )
   
→
− 1 →
− →
− 1 →
− 1
∇M ∧ = ∇ M ∧ ~j(P ) + ∇ M ∧ ~j(P ) = ∇ M ∧ ~j(P ). (5.22)
PM PM PM PM

43
6 Exercices

~ r) par le théorème d’Ampère

Exercices sur le calcul du champ B(~
1) Il faut toujours préalablement trouver la direction des lignes de champ et les invariances du
champ en se basant sur l’analyse des symétries des charges. Ensuite il faut choisir une surface fermée
fictive passant par M afin d’appliquer le théorème de Gauss.

2) Calculer le champ magnétique créé par un fil rectiligne infini parcouru par un courant I.

3) Calculer le champ magnétique créé par une bobine (infinie) sur son axe.

4) Calculer le champ magnétique créé par une nappe plane de courant surfacique ~jS (vecteur uniforme).

44
Chapitre 6

Le champ électromagnétique et ses sources

dans le vide

1 Les équations de l’électrodynamique classique

Electrodynamique classique : on considère un ensemble de N charges ponctuelles qi (généralement
en mouvement) qui interagissent entre elles via le champ électromagnétique. Les charges créent le champ
électromagnétique, et le champ agit sur les charges par les forces de Lorentz. Par conséquent, toutes les charges
interagissent entre elles via le champ : c’est un problème à N corps.

Etat des charges : à chaque instant t, l’état d’une charge ponctuelle qi est déterminé par sa position et sa
vitesse, i.e. le couple {~ri (t), ~vi (t)}, comme en mécanique classique. C’est une description discrète avec un nombre
fini de degrés de liberté.

Etat du champ électromagnétique : à chaque instant t, l’état du champ électromagnétique est déterminé
~ r, t)}, attachés à chaque point de l’espace ~r. Le champ est donc un objet
~ r, t), (B(~
par une infinité de vecteurs {E(~
extrêmement non local, car ~r est une variable continue qui décrit tout l’espace (à la différence de ~ri (t)). Cette
description continue implique un nombre infini de degrés de liberté.

Equation de la dynamique pour le champ électromagnétique : les équations de Maxwell (publiées

dans les années 1864-66) sont la synthèse des lois de l’électrodynamique découvertes par Coulomb, Laplace,
Gauss, Ampère, Faraday et quelques autres. Elles gouvernent l’évolution temporelle du champ électromagnétique
en présence de densité de charge et de courant :

→
− ~
∇ · B(~r, t) = 0, (6.1)
→
− ~ ~ r, t)
∂ B(~
∇ ∧ E(~r, t) = − , (6.2)
∂t
→
− ~ 1
∇ · E(~r, t) = ρ(~r, t), (6.3)
0
→
− ~ ~ r, t)
∂ E(~
∇ ∧ B(~r, t) = µ0~j(~r, t) + µ0 0 , (6.4)
∂t

avec 0 = 8, 85.10−12 F.m−1 et µ0 = 4π.10−7 H.m−1 , dans le système international.

1) Ces équations sont des équations aux dérivées partielles (EDP) qui relient entre elles les dérivées spatiales
et temporelles des composantes des champs, ainsi que les densités de charge ρ(~r, t) et de courant électrique
~j(~r, t).
2) On distingue les équations homogènes (dites équations de structure du champ) et les équations qui
contiennent les sources ρ(~r, t) et ~j(~r, t).
3) Ces EDP sont linéaires dans le sens qu’elles renferment le principe de superposition suivant : si l’on
additionne deux distributions sources, 1 et 2, selon ρ(~r, t) = ρ1 (~r, t) + ρ2 (~r, t) et ~j(~r, t) = ~j1 (~r, t) + ~j2 (~r, t),
~ r, t) = E
alors les champs correspondants sont E(~ ~ 1 (~r, t) + E
~ 2 (~r, t) , et B(~
~ r, t) = B~ 1 (~r, t) + B ~ 2 (~r, t).

45
 4) Attention, cela ne signifie pas que les équations de Maxwell ne décrivent que les phénomènes linéaires.
En effet pour décrire complètement un phénomène de couplage du champ électromagnétique et de la
matière, il faut compléter les équations de Maxwell par une relation constitutive qui relie les sources
charges et courants aux champs. Si la relation constitutive est non linéaire en champs, alors les équations
obtenues pour les champs seuls deviennent non linéaires en champ.
5) Les champs électriques et magnétiques sont couplés via les dérivées temporelles.

Equation de la dynamique des charges ponctuelles : le mouvement des charges ponctuelles obéit aux
équations de Newton(-Lorentz) qui s’écrivent ici :

d~vi (t) ~ ri (t), t) + ~vi (t) ∧ B(~

~ ri (t), t)),
m = qi (E(~ (6.5)
dt

en l’absence de toute autre force que les forces de Lorentz. Il s’agit d’équations différentielles ordinaires (c’est à
dire impliquant des dérivées totales par rapport au temps).

Problème complet auto-cohérent : Toutes les charges (mobiles ou immobiles) contribuent et donnent
naissance au champ électromagnétique, qui "à son tour" va agir sur le mouvement de ces charges. On a donc un
problème self-consistent (auto-cohérent) très complexe à résoudre. Les charges interagissent entre elles via le
champ : c’est le problème à N corps qui n’a pas de solution exacte. On doit donc comme toujours en physique
utiliser des approximations pour le résoudre.

Relations constitutives : Il est important de comprendre qu’il n’existe pas une seule, ni une meilleure
approximation que les autres. En fait, il existe plusieurs types d’approximations, chacune adaptée à une grande
classe de matériaux, telles que : les conducteurs, les isolants, les supraconducteurs , les ferromagnétiques, etc...
C’est en faisant certaines approximations sur ces équations que l’on obtient une relation constitutive entre les
champs et les sources. On verra que l’on peut classer les matériaux selon leurs propriétés électromagnétiques, et
qu’il existe un type de relation constitutive pour chaque catégorie.

Exercice 1 : Densités volumiques de charge et de courant

Les charges électriques sont discrètes mais doivent apparaître dans les équations de Maxwell sous la forme de
champs sources ρ(~r, t) et ~j(~r, t).
1) Rappeler les définitions de ρ(~r, t) et ~j(~r, t), et leurs unités.
2) Exprimer ρ(~r, t) et ~j(~r, t) en fonction des positions {~ri (t)} et des vitesses {~vi (t)} des charges ponctuelles.

Exercice 2 : Conservation de la charge et équations de Maxwell

1) Généraliser l’équation de conservation de la charge ci-dessus pour des densités de charges ρ(~r, t), et de
courants ~j(~r, t) quelconques à 3D.
2) Montrer que les équations de Maxwell impliquent (automatiquement) la conservation de la charge.
Commenter sur le rôle du courant de déplacement.

2 Découpler les équations de Maxwell

Les équations de Maxwell relient les dérivées premières des composantes des champs électrique et magnétique.
En manipulant ces équations, on peut obtenir d’autres EDP s’appliquant uniquement au champ électrique (ou

46
uniquement au champ magnétique).

Exercice 3 : Des équations de Maxwell vers les équations de d’Alembert

1) En utilisant la formule d’analyse vectorielle du rotationnel du rotationnel, obtenir une EDP pour le
champ électrique seul.
2) Obtenir l’équation similaire pour le champ magnétique.
3) Que deviennent ces équations découplées dans une région de l’espace vide de charges et de courants ?
Quelle est la conséquence très importante de ces équations ? Quelle application numérique s’impose ?
4) A-t-on réellement ainsi réussi à découpler les phénomènes électriques et magnétiques ?

3 Les solutions des équations de Maxwell via les potentiels

On fait apparaître des intermédiaires de calcul que l’on appelle des potentiels (scalaire et vecteur), et qui
obéissent aussi à des EDP. Les solutions de ces EDP s’écrivent en fonction des distributions de charge et courant
sources.

~ r, t)} vers les potentiels {V (~r, t), (A(~

~ r, t), (B(~
Exercice 4 : Des champs {E(~ ~ r, t)}

1) A partir des équations de Maxwell dites de structure (équations homogènes), rappeler comment on
~ r, t).
introduit le potentiel scalaire V (~r, t) et le potentiel vecteur A(~
2) A partir des équations de Maxwell restantes, établir les équations que satisfont les potentiels V (~r, t) et
~ r, t).
le potentiel vecteur A(~
3) En quoi la description en termes de potentiels peut sembler plus simple que la description en termes de
champs ? (compter le nombre de fonctions scalaires et composantes nécessaires dans chaque description).
4) Qu’appelle-t-on l’invariance de jauge en électrodynamique classique ?
5) Dans la jauge de Lorenz, comment s’écrivent les équations 2). Est-ce que ces équations sont équivalentes
~ r, t) sont réellement
à celles de Maxwell ? Est-ce que les potentiels V (~r, t) et le potentiel vecteur A(~
découplés ?

Les solutions de ces équations de Poisson sont connues en mathématiques et s’écrivent :

ρ(P, t − P M/c)
ZZZ
1
V (M, t) = dτ , (6.6)
4π0 P ∈V PM

µ0
ZZZ ~j(P, t − P M/c)
~
A(M, t) = dτ . (6.7)
4π P ∈V PM
Pourquoi appelle-t-on ces équations les solutions des potentiels retardés ?

Remarque : en fait, en physique relativiste et en physique quantique ces potentiels sont plus que de
simples intermédiaires de calcul et jouent un rôle encore plus fondamental que les champs eux-mêmes. Quel effet
quantique démontre la "réalité" du potentiel vecteur (en dehors de son effet via la force de Lorentz) ?

4 Energie électromagnétique dans le vide

L’énergie totale du système {champ+particules chargées} est conservée. Cependant le champ peut céder de
l’énergie aux particules chargées via le travail des forces de Lorentz, et les particules chargées accélérées peuvent
rayonner et donc fournir de l’énergie au champ.

47
4.1 Densité de puissance cédée par le champ à la matière chargée

On considère un élément de volume dτ qui contient des charges mobiles de densité ρm se déplaçant à la
vitesse moyenne ~v . La charge mobile dq = ρm dτ reçoit du champ la puissance :

~ v = ~j.Edτ,
dP = (ρm dτ E).~ ~ (6.8)

car la force magnétique de Lorentz ne travaille pas (force orthogonale à la vitesse). Remarquez aussi que la force
électrique ne travaille que sur les charges mobiles.

4.2 Densité locale d’énergie du champ u(~r, t)

L’énergie électromagnétique est localisée. Un volume dτ , centré autour de la position ~r, contient l’énergie
électromagnétique :

dU = u(~r, t)dτ. (6.9)

Cette relation définit la densité volumique d’énergie électromagnétique u(~r, t), qui s’exprime en J.m−3 .

4.3 ~ r, t)
Vecteur de Poynting R(~

Le vecteur de Poynting R ~ est le vecteur dont le flux à travers la surface infinitésimale dS est égal à la
puissance rayonnée par le champ à travers dS :

~ r, t).−
dPray = R(~
→
dS, (6.10)

et la norme du vecteur de Poynting s’exprime donc en W.m−2 .

4.4 Loi de conservation de l’énergie

Exercice 5 : Loi de conservation de l’énergie électromagnétique

On recherche une expression locale de la conservation de l’énergie.
1) Expliquer qualitativement l’idée de conservation de l’énergie électromagnétique.
2) On considère un volume macroscopique V. Exprimer l’énergie U (t) contenue dans ce volume à la date
t, ainsi que U (t + dt).
3) Ecrire l’énergie électromagnétique rayonnée vers l’extérieur de V entre t et t + dt.
4) Ecrire l’énergie cédée par le champ à la matière entre t et t + dt.
5) En déduire un bilan exprimant la conservation de l’énergie.

4.5 Expressions de u(~r, t) et du vecteur de Poynting en fonction des champs

En mécanique du point, on peut établir le théorème de l’énergie cinétique ou totale à partir de l’équation de
la dynamique de Newton. Pour le champ, les équations de la dynamique sont les équations de Maxwell et elles

48
impliquent aussi une loi de conservation de l’énergie.

Exercice 6 : Expression de la densité d’énergie électromagnétique et du vecteur de Poynting

On veut montrer que les équations de Maxwell impliquent directement l’équation de la conservation de l’énergie
établie à l’exercice 5.
~ r, t) ∧ B(~
1) Calculer la divergence de (E(~ ~ r, t)) en utilisant la formule d’analyse vectorielle :

~ ∧ B)
∇.(E ~ = −E.(∇
~ ~ + B.(∇
∧ B) ~ ~
∧ E). (6.11)

~
2) Exprimer E.(∇ ~ et B.(∇
∧ B) ~ ~ en utilisant les équations de Maxwell.
∧ E)
3) En déduire une expression de type loi de conservation locale.
~ r, t).
4) Par identification, donner les expressions de u(~r, t) et R(~

49
Chapitre 7

Ondes électromagnétiques dans le vide

Il est remarquable que les ondes électromagnétiques puissent exister sans support matériel contrairement
aux autres types d’ondes. Le vide est le "milieu" le plus simple. Ses relations constitutives s’expriment par
l’absence de charges et de courants : ρ = 0 et ~j = ~0. Cependant, il existe des solutions non triviales dans un tel
vide : les ondes électromagnétiques. Evidemment cela suppose que loin de la région vide des sources, des charges
accélérées, ont donné naissance à ces ondes.

1 Relation de dispersion

La relation de dispersion relie les périodes spatiales et temporelles d’une onde. La périodicité spatiale est
caractérisée par la longueur d’onde λ, ou de manière équivalente par le module du vecteur d’onde k = 2π/λ. La
périodicité temporelle est caractérisée par la période T , ou de manière équivalente par la pulsation ω = 2π/T
(ou encore la fréquence ν = 1/T ).
La relation de dispersion est une notion très générale de la physique ondulatoire qui n’est ni restreinte aux
ondes électromagnétiques, ni à la physique classique. En effet, il existe des relations de dispersion pour les
différents types d’ondes à la surface d’un fluide, pour les ondes acoustiques, pour les phonons dans les solides,
etc... De plus, via la dualité onde-corpuscule, cette notion s’étend aussi à la physique des particules et à la
physique des solides. Par exemple, on peut parler de la relation de dispersion d’un électron dans le vide, de la
relation de dispersion d’un électron dans un semiconducteur, etc...
De manière générale, la relation de dispersion indique comment se propage une onde dans un milieu donné
(qui peut être le vide dans le cas des ondes électromagnétiques).

Phénomène de dispersion (milieu dispersif ) : on dit qu’un milieu est dispersif lorsque les ondes de
différentes pulsations (ou longueurs d’ondes) s’y propagent à des vitesses différentes. Un milieu sera non
dispersif lorsque les ondes de différentes pulsations s’y propagent à la même vitesse vϕ , ce qui correspond à
une relation de dispersion linéaire du type ω = vϕ k.

Méthode pour obtenir la relation de dispersion :

~
On injecte une solution de type onde plane (du type ei(k.~r−ωt) en notation complexe) dans l’équation d’onde.
Ceci implique une relation entre ω et le module du vecteur d’onde k.
Note : le préfacteur peut être une amplitude scalaire, vectorielle,... on peut même parfois dans les problèmes
couplés avoir un système d’équations linéaires homogènes dont la relation de dispersion s’obtient en écrivant
la nullité du déterminant.

Propagation des EMs dans le vide : il est remarquable que les ondes EM puissent se propager dans
le vide. Dans le vide, la relation de dispersion est très simple. Qualitativement, on peux comprendre que le
vide (continu, isotrope et homogène sans échelle de longueur caractéristique) ne puisse pas privilégier certaines

50
longueurs d’onde.

Exercice 1 : relation de dispersion d’une onde dans le vide

1) Rappeler l’équation de d’Alembert vérifiée par le champ magnétique dans le vide.
2) En déduire la relation de dispersion d’une onde électromagnétique dans le vide.
3) Le vide est-t-il un milieu dispersif ou non dispersif pour les ondes EM ?

2 Structure des ondes et modes de Fourier

L’onde électromagnétique est un champ vectoriel (il y a même deux vecteurs en chaque point). Pour obtenir
la structure d’une OM il faut revenir aux équations de Maxwell couplées pour savoir comment les vecteurs
~ r, t), B(~
E(~ ~ r, t) sont orientés les uns par rapport aux autres.
Toute onde peut se décomposer en modes de Fourier selon les sommes continues :

→
− d~k dω →
− ~ →
− d~k dω →
− ~
Z Z
~ ~
E (~r, t) = 3
E (k, ω)ei(k.~r−ωt) , B (~r, t) = 3
B (k, ω)ei(k.~r−ωt) , (7.1)
(2π) 2π (2π) 2π

→
−
où E (~k, ω) désigne l’amplitude du mode de Fourier caractérisé par le vecteur d’onde ~k et la pulsation ω. Un
tel mode de Fourier est aussi appelé OPPM pour onde plane progressive et monochromatique, de "profil" en
~
ei(k.~r−ωt) . Le qualificatif "onde plane progressive" correspond au fait que le vecteur d’onde ~k est unique et bien
défini. Le qualificatif "monochromatique" signifie que la pulsation ω est bien définie. Bien sûr en toute rigueur
une OPPM est un objet idéal puisqu’elle devrait être infiniment étendue dans l’espace et le temps à la fois.
Cependant, on peut donner deux types d’arguments pour répondre à l’objection de la phrase précédente. Primo,
dans certaines conditions expérimentales, on peut se rapprocher dans une région de l’espace des conditions
d’une OPPM par exemple en réalisant un faisceau de lumière bien collimaté et constitué de rayons parallèles
et en restreignant l’étude à une zone où les surfaces d’ondes sont planes localement. Secondo, toute onde de
profil spacial et temporel arbitraire peut toujours être décomposer en OPPM. L’OPPM est donc une "brique
élémentaire" très intéressante : si on connaît son comportement dans un milieu donné, on peut en déduire ce
qu’il va se passer pour un paquet d’OPPM quelconque. Cela est du à la linéarité des équations de Maxwell dans
le vide. Même en optique non-linéaire, on peut obtenir beaucoup d’information en analysant la structure d’un
mode de Fourier particulier puis comment ces différents modes se combinent à cause des processus non linéaires
→
− →
−
(voir chapitre sur l’optique non linéaire). Notons que E (~r, t) et B (~r, t) sont des champs à composantes réelles
→
− →
−
donc il faut respecter certaines conditions sur les coefficients E (~k, ω) et B (~k, ω) afin que l’intégrale sur tous les
modes (terme de droite) soit réelle.
Pour écrire les équations de Maxwell dans l’espace de Fourier, il faut aussi exprimer toutes les grandeurs qui
y interviennent comme des superpositions de modes propres. C’est le cas des sources notamment. Ainsi il faut
écrire :

d~k dω ~ →
− d~k dω →− ~
Z Z
i(~k.~
r−ωt) ~
ρ(~r, t) = 3
ρ( k, ω)e , j (~
r , t) = 3
j (k, ω)ei(k.~r−ωt) , (7.2)
(2π) 2π (2π) 2π

→
−
où ρ(~k, ω) et j (~k, ω) désignent respectivement les amplitudes du mode de Fourier (~k, ω) de la densité de charge
électrique, et de la densité de courant électrique.

Equations de Maxwell en représentation de Fourier pour l’espace et le temps : mathématiquement,

les OPPMs ont un comportement remarquablement simple par rapport à la dérivation spatiale et temporelle.
En effet, la dérivation par rapport à l’espace correspond à la simple multiplication par le vecteur d’onde, et la
dérivation par rapport au temps correspond à la simple multiplication par la pulsation. Ainsi les équations aux
dérivées partielles (très compliqué) deviennent de simples équations algébriques pour un mode de Fourier donné :

51
 →
− ~ ~
i k · B( k, ω) = 0, (7.3)
→
− ~ ~ ~ ~k, ω),
i k ∧ E(k, ω) = iω B( (7.4)
→
− ~ ~ 1 ~
i k · E( k, ω) = ρ(k, ω), (7.5)
0
→
− ~ ~
i k ∧ B( ~ ~k, ω),
k, ω) = µ0~j(~k, ω) − iωµ0 0 E( (7.6)
(7.7)

Cette représentation est adaptée pour analyser la structure d’une onde dans le vide où elle devient encore
plus simple car : ρ(~k, ω) = 0 et ~j(~k, ω) = ~0. On voit que les champs sont transverses dans le vide et que
{~k, E(
~ ~k, ω), B(
~ ~k, ω)} forment un trièdre trirectangle. En introduisant le vecteur unitaire ~u défini par ~k = k~u (où
k > 0 est le module du vecteur d’onde), on peut exprimer le champ magnétique en fonction du champ électrique
ou inversement :

→
− →
−
~ ~k, ω) =k ~ ~k, ω) = u ∧ E(
~ ~k, ω),
B( ∧ E( (7.8)
ω c
→
−
~ ~ c2 k ~ ~k, ω) = −c→
− ~ ~k, ω),
E(k, ω) = − ∧ B( u ∧ B( (7.9)
ω
(7.10)

Equations de Maxwell en représentation de Fourier pour l’espace (seulement) : dans les problèmes
avec une condition initiate temporelle, on peut avoir intérêt à utiliser une représentation mixte qui détermine
l’évolution temporelle d’un mode de Fourier spatial de vecteur d’onde bien défini :

→
− ~ ~
i k · B(k, t) = 0, (7.11)
→
− ~ ~ ~ ~k, t)
∂ B(
i k ∧ E( k, t) = − , (7.12)
∂t
→
− ~ ~ 1 ~
i k · E(k, t) = ρ(k, t), (7.13)
0
→
− ~ ~ ~ ~k, t)
∂ E(
i k ∧ B( k, t) = µ0~j(~k, t) + µ0 0 , (7.14)
∂t
(7.15)

Equations de Maxwell en représentation de Fourier pour le temps (seulement) : dans les pro-
blèmes avec des conditions de confinement spatial, on peut avoir intérêt à utiliser une représentation mixte qui
détermine l’évolution spatiale d’un mode de Fourier harmonique de pulsation bien définie ω :

→
− ~
∇ · B(~r, ω) = 0, (7.16)
→
− ~ ~ r, ω),
∇ ∧ E(~r, ω) = iω B(~ (7.17)
→
− ~ 1
∇ · E(~r, ω) = ρ(~r, ω), (7.18)
0
→
− ~ ~ r, ω),
∇ ∧ B(~r, ω) = µ0~j(~r, ω) − iωµ0 0 E(~ (7.19)
(7.20)

52
3 Etats de polarisation
Dans le vide, l’onde est transverse. Elle a donc deux composantes perpendiculaires à la direction de
propagation. La polarisation de l’onde concerne l’étude de la différence de phase φ entre ces deux composantes.

Méthode pour obtenir l’état de polarisation d’une OPPM

1) On se place dans un plan d’onde, par exemple x = 0 et on regarde la courbe décrite par l’extrémité du
~ = 0, t) au cours du temps dans ce plan.
vecteur E(x

2) Suivant les valeurs de φ, on obtient divers états de polarisation. Retrouver les différents cas.

3) On peut toujours décomposer un état de polarisation dans la base de deux états de polarisations
indépendants. Par exemple comment exprimeriez vous un état de polarisation rectiligne en fonction
des états de polarisation circulaires ?

4) Pour modifier l’état de polarisation, on fait passer l’OPPM dans une lame d’un matériau "actif" qui
introduit une différence de phase supplémentaire entre les composantes. Si cette lame est comprise
entre x = 0 et x = e, on regarde quelle est la nouvelle différence de phase en x = e pour obtenir l’état
de polarisation en sortie.

4 Aspects énergétiques d’une OPPM

Il faut savoir calculer la densité d’énergie et le vecteur de Poynting pour le cas particulier d’une onde plane
progressive monochromatique se propageant dans la direction ~uz . Attention les grandeurs énergétiques sont
quadratiques donc on ne peut pas utiliser la représentation complexe. Il faut revenir aux parties réelles pour
exprimer les grandeurs énergétiques instantanées, et calculer ensuite leurs valeurs moyennes sur une période de
l’onde.
Densité volumique d’énergie instantanée : dans une OPPM, les parties électrique et magnétique portent
la même énergie donc on peut écrire la densité d’énergie instantanée :
0 ~ 2 1 ~2 ~ 2 (~r, t) = 1 B~ 2 (~r, t).
uem (~r, t) = E (~r, t) + B (~r, t) = 0 E (7.21)
2 2µ0 µ0

Vecteur de Poynting instantané :

!
~ ~ ~ ~ r, t) ~ 2 (~r, t)
~ r, t) = E(~r, t) ∧ B(~r, t) = E(~r, t) ∧
R(~
~uz ∧ E(~
=
E ~ 2 (~r, t)~uz ,
~uz = 0 cE (7.22)
µ0 µ0 c µ0 c

Dans une OPPM, on a une relation à tout instant entre la vecteur de Poynting et la densité d’énergie
électromagnétique :
~ r, t) = c uem (~r, t)~uz ,
R(~ (7.23)
qui indique que de l’énergie est transportée dans la direction de propagation à la vitesse c.
Valeurs moyennes temporelles : Prenons une OPPM polarisée suivant ~ux qui s’écrit en notation valeurs
réelles :
~ r, t) = E0 cos(kz − ωt)~ux .
E(~ (7.24)
La valeur moyenne temporelle du vecteur de Poynting et de la densité d’énergie s’écrivent :

~ = 0 c E02 E02
hRi ~uz , huem i = 0 , (7.25)
2 2
et sont indépendantes de la position (ceci est due à la structure "idéale" d’OPPM).

53
Chapitre 8

Champ électrique dans les conducteurs

Dans ce chapitre, on insiste sur la notion de relation constitutive qui va caractériser la façon dont le matériau
(ici un conducteur) interagit/répond à un champ électrique. L’idée est d’introduire une relation entre les sources
(charges et courants) et les champs électromagnétiques afin de rendre les équations fermées. On obtient ainsi le
comportement effectif du champ électromagnétique dans un milieu. La relation constitutive peut être devinée
de manière phénoménologique en s’inspirant des résultats expérimentaux, ou bien justifiée de manière plus
microscopique. Ici notre modélisation microscopique restera très très élémentaire, mais elle sera suffisante pour
aboutir à la forme correcte de la relation constitutive.
Nous analyserons la physique résultante dans deux cas très différents : avec courant électrique continu et
sans courant. Dans les deux cas, la nature conductrice du matériau est exprimée par une relation linéaire entre le
courant et le champ électrique. Nous verrons dans un chapitre suivant la situation plus complexe d’un conducteur
parcouru par des courants alternatifs qui sera forcément le siège de champs électrique et magnétique couplés.

1 Relation constitutive courant-champ électrique

2 Origine microscopique de la relation constitutive (Modèle de Drude)

Le champ électrique est maintenu par la pile, et cause le mouvement de dérive des électrons. On peut faire
un petit modèle très simple :
d~v ~ − m ~v ,
m = qE (8.1)
dt τ
où on reconnaît la force exercée par le champ électrique sur le porteur de charge, et où une force dissipative a
été introduite pour traduire les collisions dans le matériau qui limite la vitesse de dérive des porteurs. La force
qE~ i électrique tend à accélérer les électrons de manière uniforme, tandis que la force dissipative −m~v /τ tend à
réduire leur vitesse. Qualitativement, on peut en déduire la possibilité d’un état stationnaire dans lequel la force
dissipative compense exactement la force électrique.
Dans ce régime stationnaire, les deux forces se compensent et on obtient une vitesse de dérive du porteur :

~
q Eτ
~v = , (8.2)
m
En réalité il s’agit d’une vitesse moyenne sur un grand nombre de porteurs de charges présents avec une densité
n dans le conducteur. Cela correspond ainsi à une densité de courant volumique :

2 ~
~j = nq~v = nq τ E , (8.3)
m
Ainsi le courant est proportionnel au champ électrique qui créé le mouvement, ce qui est logique. L’effet est
linéaire par rapport à la cause :
~ .
~j = σ0 E (8.4)

Le coefficient de proportionnalité s’appelle la conductivité σ. Elle caractérise la nature intrinsèque conductrice

du matériau : elle est grande si la densité de porteurs est grande (q et m sont en général fixés, charge et masse de
l’électron). On a ainsi un modèle microscopique qui fournit la forme locale de la loi d’Ohm. Cette loi aurait aussi

54
pu être postulée de manière phénoménologique, comme traduisant la réponse linéaire en courant des porteurs de
charge à un champ électrique (entretenu par un générateur) :

→
− 1 ~
E = j = ρ0~j , (8.5)
σ0

où ρ est la résistivité du matériau. Cependant dans ce cadre phénoménologique, ρ0 (ou σ0 ) devrait être déterminée
par les expériences seules, et nous n’aurions pas la relation avec les paramètres microscopiques que sont m, e, τ ,
et n.

2.1 Passage à la loi d’Ohm macroscopique (1826) : conductance et résistance

On obtient la loi d’Ohm macroscopique en intégrant le long des lignes/tubes de courant de A vers B :
→
− ~ ~ = ρ0 I dl,
E · dl = ρ0~j · dl (8.6)
S
ce qui donne :
V (A) − V (B) = RiAB . (8.7)
La constante R donnée par
L
R = ρ0 , (8.8)
S
est la résistance du fil de longueur L et de section S. On peut enfin se demander pourquoi ce fluide/gaz d’électrons
libres ne s’écoule pas en dehors du fil de cuivre. Cela peut s’expliquer qualitativement comme il suit. Au coeur
du métal, une charge négative est toujours compensée par le fond des ions positifs immobiles du réseau cristallin.
Si un électron s’échappe et sort du métal, cela va créer un excédent de charges positives dans le métal : cette
charge positive exerce une force de rappel sur l’électron ce qui l’empêche de s’évader. On dit que les électrons
sont confinés dans fil. Pour arracher un électron hors d’un métal il faut fournir a cet électron une énergie de
l’ordre de quelques eV (travail d’extraction).

3 Conducteurs en équilibre électrostatique

Dans ce paragraphe, nous allons nous intéresser à la situation d’un conducteur/métal à l’équilibre électrosta-
tique, c’est à dire en l’absence de courant et de tout mouvement d’ensemble des charges. Cela paraît simple
à priori : il suffit de couper l’alimentation (enlever la pile ou ouvrir un interrupteur). En réalité il existe des
situations plus complexes : par exemple on peut approcher un matériau isolant chargé à proximité de la surface
d’un métal. Dans ce cas les charges de l’isolant vont attirer ou repousser (selon leur signe) les électrons libres du
métal provoquant ainsi un mouvement, un réarrangement des électrons au sein du métal jusqu’à obtention d’un
état d’équilibre où toutes les charges sont à nouveau immobiles.

3.1 L’équilibre électrostatique

L’équilibre électrostatique est établi lorsque plus aucune charge ne bouge, donc en particulier en l’absence de
tout courant électrique : I = 0 et ~j = ~0. Cela implique l’absence de circuit bouclé alimenté par une pile. On sait
qu’un métal est un composé dans lequel il existe des porteurs de charge (les électrons libres) qui se mettent
en mouvement dès qu’apparaît le moindre petit champ électrique. Donc pour qu’il n’y ait pas de courant, il
faut que le champ électrique soit identiquement nul dans tout le volume du conducteur. A l’intérieur d’un
conducteur en équilibre électrostatique, le champ électrique est nul :
→
−
E i = ~0 . (8.9)

Conséquence 1 : Le potentiel électrostatique est le même en tout point d’un conducteur en équilibre
électrostatique. Un conducteur en équilibre électrostatique est un volume équipotentiel.

V (M ) = const . (8.10)

Conséquence 2 : En utilisant le théorème de Gauss, on doit conclure à l’absence de charges en volume dans
le métal. Les charges ne peuvent exister que bloquées/plaquées contre les bords du métal : interface métal/vide
extérieur.

55
3.2 Interface vide/métal
On peut appliquer le théorème de Gauss au voisinage d’une surface vide/métal. Le champ est identiquement
nul du coté métal. En revanche il peut exister un champ électrique non nul à l’extérieur du métal. Prenons une
surface plane définie par z = 0. Métal z < 0 et vide z > 0 (ici on aura donc ~uz = ~next ).
Comme le métal est une équipotentielle, on peut montrer que le champ électrique doit être perpendiculaire à
la surface du métal donc selon ~uz = ~next .
Prenons maintenant comme surface/boîte de Gauss un parallélépipède à cheval sur l’interface métal/vide :
→
− →
−
E (z = 0+ )Stop ~uz + E (z = 0− )(−Sbottom ~uz ) + 0 + 0 = Qint /0 , (8.11)

D’où :
→
− σ
E (z = 0+ ) = ~next , (8.12)
0
→
−
car E (z = 0− ) = ~0, avec σ la densité surfacique (en C.m−2 ) de charges accumulée sur l’interface métal/vide. On
appelle ce résultat le théorème de Coulomb.

3.3 Boule métallique dans un champ électrique statique

On considère donc une boule conductrice de centre O et de rayon R. On applique un champ uniforme
~
E0 = E0 ~uz , par exemple réalisé entre les armatures d’un condensateur plan (prendre E0 > 0 pour fixer les idées
et les dessins). On veut déterminer comment les lignes de champ électrique et les équipotentielles sont modifiées
par la présence de la boule conductrice. Il est commode de calculer le potentiel électrique plutôt que le champ.
Le potentiel électrique obéit à l’équation de Poisson ou Laplace. On va donc chercher le potentiel électrique
correspondant à une boule conductrice immergée dans un champ extérieur uniforme E ~ 0 en résolvant l’équation
~
de Laplace avec les bonnes conditions aux limites. L’effet de E0 est aussi de créer une distribution de charges
à la surface de la boule, car les charges libres du métal sont libres de se déplacer jusqu’à être bloquer par les
limites de la boule. En revanche on suppose que les armatures du condensateur qui créé E ~ 0 = E0 ~uz sont très
éloignées de sorte que les charges sur ces armatures ne sont pas affectées.

1) Démontrer l’équation de Poisson en partant des équations de Maxwell de l’électrostatique :

ρ
∆V + = 0, (8.13)
0
où ρ est la densité volumique de charges.
2) Que devient l’équation de Poisson dans l’espace r > R vide entourant la boule ? Faire alors la séparation
des variables V (r, θ) = f (r)g(θ) en utilisant l’expression suivante du Laplacien :
   
1 ∂ ∂ 1 ∂ ∂
∆= 2 r2 + 2 sin θ , (8.14)
r ∂r ∂r r sin θ ∂θ ∂θ

Justifier l’omission (volontaire) de la coordonnée ϕ.

3) Justifier que g(θ) = A cos θ.
4) Ecrire l’équation différentielle pour la partie radiale f (r). Chercher les solutions en loi du puissance du
type f (r) = rp . Quelles sont les valeurs de p possibles ?
5) Justifier, d’après ce qui précède, que le potentiel à l’extérieur de la boule métallique (r > R) peut s’écrire
de la manière suivante.

V ext (r, θ) = Cr−2 + Dr cos θ .

(8.15)
Il faut maintenant déterminer les coefficients inconnus C et D.
6) Quand r → ∞, l’influence de la sphère doit disparaître et on doit retrouver le potentiel correspondant au
~ 0 = E0 ~uz . En déduire la valeur de D.
champ uniforme E
7) Calculer le le champ électrique à l’extérieur de la boule.
8) Que vaut le champ électrique à l’intérieur de la boule ? Exprimer la condition de passage en r = R sur
le champ électrique. Justifier que la composante tangentielle du champ électrique est toujours continue
même en présence de charges surfaciques. En déduire la valeur de C.

56
 9) Analyser la discontinuité du champ électrique en r = R, et en déduire la densité surfacique de charge qui
apparaît à la surface de la boule. Peut-il y avoir une densité de charge volumique non nulle à l’intérieur
de la boule ? Pourquoi ?
10) Dessiner l’allure des lignes de champ et des équipotentielles autour de la boule.
11) Que vaut le potentiel sur la boule conductrice ?

4 Condensateurs
Un condensateur est un ensemble de deux armatures métalliques séparées par un isolant (ou du vide). Il
permet de stocker du champ électrostatique et l’énergie électrostatique associée. Historiquement, avant l’invention
de la pile par Volta, les pionniers de l’électricité utilisaient des bouteilles de Leyde (les premiers condensateurs)
pour stocker de l’électricité en vue d’une décharge ultérieure. Dans un condensateur plan, il est possible de
réaliser un champ électrique uniforme. On peut voir ce condensateur comme une boîte à champ électrique.

4.1 Condensateur à l’équilibre électrostatique

On obtient un condensateur quand on présente deux surfaces métalliques en regard et portant des charges
opposées. En connectant deux surfaces métalliques séparées/isolées (entre elles) à une machine électrostatique,
Von Kleist (Allemagne, 1745) et Musschenbroek (Pays Bas à Leyde, 1746) réalisèrent (indépendamment) les
premiers condensateurs, qui permirent de stocker de l’électricité (en fait des charges) pour une utilisation
ultérieure (lors de la décharge). On appela ce dispositif la bouteille de Leyde. De nos jours (et depuis l’invention
de la pile en 1800), on peut charger un condensateur en connectant simplement les deux plaques aux deux bornes
d’un générateur continu. On peut décharger le condensateur en reliant les deux armatures par un conducteur
(court-circuit).
Prenons la géométrie la plus simple : deux armatures planes (de surface S) et parallèles séparées par du vide
sur une épaisseur d. L’armature dite positive (notée P ) porte la charge Q > 0, tandis que l’armature négative
(N ) porte la charge −Q, de manière à ce que le système soit globalement neutre. On s’intéresse au champ créé
par ces charges dans l’espace vide entre les armatures décrites comme des plans infinis portant des densités
surfaciques uniformes σ = Q/S et −σ. C’est une bonne approximation à condition d’être suffisamment loin des
bords des plaques. En appliquant le théorème de Gauss, on obtient le champ électrique :

~ = σ u~z = Q/S u~z .

E (8.16)
0 0
On veut maintenant savoir quelle est la différence de potentiel électrique qui se développe entre les plaques
de ce condensateur plan. Pour cela, on évalue la circulation du champ entre les deux armatures (avec u~z orienté
de P vers N ). La circulation infinitésimale s’exprime comme :

~ = Q/S u~z .dl.

~ dl
E. ~ (8.17)
0

On circule de la plaque P vers la plaque N sur le chemin qui nous arrange (cela est possible car la circulation
du champ électrostatique est indépendante du chemin suivi, et de plus chaque plaque métallique est une
équipotentielle à l’équilibre) :
Z N Z N
E. ~ = V (P ) − V (N ) = Q/S
~ dl ~ = Qd .
u~z .dl (8.18)
P 0 P S0

Le résultat est très simple : un champ uniforme dont la norme est égale a la différence de potentiel U =
V (P ) − V (N ) appliquée entre les deux armatures divisée par la distance entre ces armatures. Il est très pratique
d’avoir un champ uniforme si l’on veut faire des expériences de physique fondamentale dans lesquelles on sonde
un échantillon en lui appliquant un champ extérieur. On définit la capacité comme le rapport entre la charge
stockée et la différence de potentielle nécessaire à la stocker :

Q Q 0 S
C= = = . (8.19)
U V (P ) − V (N ) d

L’unité de capacité est le farad (F) en hommage à Michael Faraday.

57
 Figure 8.1 – Introduction d’un échantillon métallique dans un condensateur.

4.2 Conducteur placé dans un champ extérieur

Imaginons que l’on prenne un échantillon conducteur (non connecté à une pile), et que l’on veuille y faire
pénétrer un champ électrique extérieur E ~ e . On va appliquer sur cet échantillon conducteur (passif) un champ
électrique extérieur, ce dernier étant créé par un condensateur. En pratique, on introduit l’échantillon métallique
dans le condensateur. Que se passe-t-il dans l’échantillon métallique ? Tout d’abord, les porteurs libres sentent
le champ E ~ e via la force électrique q E
~ e . Cette force pousse les électrons vers les bords de l’échantillon : il y
a accumulation des électrons sur un bord et déficit sur le bord opposé. Ces charges accumulées créent une
sorte de "contre-condensateur" qui produit son propre champ de réponse de la matière E ~ rep de sens opposé au
champ extérieur initialement appliqué. La norme de ce champ de réponse E ~ rep augmente avec l’accumulation
progressive des charges sur les bords de l’échantillon. Le phénomène se stabilise lorsque E ~ tot = E
~e + E~ rep = ~0 :
on a alors atteint l’équilibre électrostatique. En réalité, l’établissement de l’équilibre est très très rapide (temps
caractéristique de l’ordre de la picoseconde ou moins).

Protocole 1 : introduction à charge fixée

On réalise l’expérience en maintenant constantes les charges (+Q, −Q) sur les armatures du condensateur :
Q0 = Q. Pour cela, on isole le condensateur (on débranche l’alimentation) après l’avoir chargé et avant
l’introduction de l’échantillon métallique. La circulation du champ électrique d’une armature à l’autre peut
varier car les armatures sont "dans le vide" (non reliées à des références de potentiel).
Comme la charge sur les armatures est constante, le champ créé dans le vide n’a pas changé non plus après
introduction de l’échantillon de métal :
σ Q/S
E= = (8.20)
0 0
Mais attention la circulation totale du champ entre les armatures va changer car dans le métal le champ est nul,
on a donc
Q(d − a)/S
U 0 = V (P ) − V (N ) = E(d − a) + 0.a = . (8.21)
0
D’où la nouvelle capacité :
Q 0 S
C0 = 0
= . (8.22)
U d−a
La capacité augmente car pour une même charge, un voltage plus petit est nécessaire. Notons que la capacité
obtenue est la même dans les deux modes opératoires bien que les charges et les voltages aient évolués de
manières bien différentes.
Notons que malgré les notations P et N ce type de condensateur n’a pas de polarité et peut être branché
dans un sens ou dans l’autre au générateur. On a simplement ici appelé pour nos calculs, P l’armature relié au
plus du générateur DC. En revanche dans les condensateurs chimiques, il faut faire attention au branchement et
il y a un symbole + ou P sur le composant indiquant le sens de branchement.

58
 Figure 8.2 – Introduction d’un échantillon métallique dans un condensateur à tension U fixée

Protocole 2 : introduction à voltage inter-armatures U fixé :

On réalise l’expérience en maintenant constante la différence de potentiel entre les armatures du condensateur.
Cela signifie que pendant l’introduction de l’échantillon, la circulation du champ électrique d’une armature à
l’autre doit rester constante. Par contre comme les deux armatures sont reliées à une alimentation, la charge des
armatures (+Q et −Q initialement) peut changer et devenir (+Q0 ,Q0 ) par un transfert d’électrons via le circuit
électrique.
Soit a l’épaisseur (suivant z de la plaque métallique). On a donc a < d (sinon ça rentre pas !). Le champ est
nul dans le métal donc il ne contribuera pas à la circulation sur une longueur a. Il faut donc que le champ dans
les espaces vides séparant l’échantillon des armatures soit renforcé, notons le E 0 .

U = V (P ) − V (N ) = E 0 (d − a). (8.23)

Donc la charge surfacique devient :

0 U
σ 0 = E 0 0 = . (8.24)
(d − a)
Et donc la nouvelle charge sur les armatures :
U 0 S
Q0 = σ 0 S = . (8.25)
(d − a)

D’où la nouvelle capacité :

Q0 0 S
C0 == . (8.26)
U d−a
La capacité augmente car pour un même voltage, il y a davantage de charges stockées sur les armatures.

59
Chapitre 9

Effet Hall : effet d’un champ magnétique

statique et uniforme sur un conducteur

Motivation : comprendre la conduction électrique en présence d’un champ magnétique uniforme qui tend à dé-
vier les électrons. Une conséquence spectaculaire est l’effet Hall classique, surtout visible dans les semiconducteurs
(il est plus faible dans les métaux).

1 Modèle microscopique de la conduction sous champ magnétique

On cherche à comprendre l’influence d’un champ magnétique extérieur sur la conduction électrique.

1.1 Modèle de Drude avec champ magnétique

Reprenons le modèle de Drude du chapitre sur les conducteurs, et ajoutons un ingrédient : la force de Lorentz
exercée par un champ magnétique uniforme sur un électron libre (ou plus généralement un porteur de charge)
de charge q :
d~v ~ − m ~v + q~v ∧ B
~.
m = qE (9.1)
dt τ
Qualitativement, on voit que l’état stationnaire, dans lequel la force électrique et la force dissipative (de type
frottements fluides) se compensaient, va être modifié par la présence du champ magnétique. Cela va se traduire
par un non alignement entre le champ électrique et la vitesse, donc le courant local.
En écrivant que le vecteur vitesse est constant, on obtient :

~0 = q E~ − 1 ~v + q ~v ∧ B
~, (9.2)
m τ m
d’où la relation entre le champ électrique et le courant stationnaire :

~ = m ~v − ~v ∧ B
E ~. (9.3)
qτ

Or la vitesse de dérive stationnaire qui apparaît dans cette équation, est proportionnelle au courant local ~j = nq~v ,
donc :
~ = m ~j − 1 ~j ∧ B
E ~, (9.4)
nq 2 τ nq
où on voit apparaître la résistivité de Drude en champ nul ρ0 = m/(nq 2 τ ).

1.2 Tenseurs de conductivité et de résistivité

B
Ex = ρ0 jx − jy , (9.5)
nq
B
Ey = ρ0 jy + jx , (9.6)
nq
Ez = ρ0 jz (9.7)

60
La troisième relation montre que le champ magnétique suivant z ne modifie pas la dynamique selon l’axe Oz. On
a avantage à exprimer les deux premières relations sous une forme matricielle qui montre bien que la linéarité
de la relation entre le courant et le champ électrique est préservée, mais que ces deux vecteurs ne
sont plus parallèles (ils ne sont plus colinéaires).
! 
B
ρ0 − nq
 
Ex jx
= B , (9.8)
Ey nq ρ0 jy

Une autre façon d’écrire cette relation de réponse linéaire consiste à factoriser ρ0 :

− qB
    
Ex 1 m τ jx
= ρ0 qB , (9.9)
Ey mτ 1 jy

car B/(nqρ0 ) = (qB/m)τ . On voit apparaître naturellement la pulsation cyclotron qB/m.

1.3 Effet Hall classique

Dans la géométrie de l’effet Hall, on impose que jy = 0 à cause des bords latéraux de l’échantillon et/ou du
fait que l’on connecte un voltmètre (d’impédance "infinie") entre les bords latéraux. Ainsi :
! 
B
ρ0 − nq
 
Ex jx
= B , (9.10)
Ey nq ρ0 0

donc la seconde ligne donne :

B
Ey = jx , (9.11)
nq
c’est à dire un champ électrique transverse (selon y) pour un courant injecté longitudinal (selon x). On peut
calculer le voltage transverse (dit de Hall) :
Z
B
UH = V (Ly /2) − V (−Ly /2) = dV = −Ey Ly = − Ly jx , (9.12)
nq

tandis que l’intensité électrique (en convention positive pour un courant vers les x croissants) :

I = jx Ly Lz , (9.13)

donc :
B
UH = − I, (9.14)
nqLz
ce qui permet de définir la résistance de Hall RH = UH /I. Dans les conventions adoptées (courant positif selon
x croissants et tension positive selon y croissants), la résistance de Hall est positive pour des porteurs de charges
négatifs (comme les électrons) :
B
RH = − . (9.15)
nqLz
Elle serait négative pour des porteurs de charges positifs, comme les trous dans les semiconducteurs dopés p. La
mesure de la résistance de Hall permet de déterminer le signe des porteurs, ce qui est impossible avec une simple
mesure de conductance en champ nul (car ρ0 est en q 2 donc indifférent au signe de q).

2 Origine de la force de Laplace sur un conducteur

61
Chapitre 10

Champ électrique dans la matière isolante

Dans ce chapitre, nous nous intéressons au champ électrique dans la matière, c’est à dire dans des milieux
dilués ou denses constitués d’un très grand nombre d’atomes ou de molécules. Dans les deux chapitres précédents,
nous avons considérés des milieux conducteurs. Dans ce chapitre, nous allons présenter un formalisme adapté au
comportement électromagnétique des matériaux isolants. Un isolant est un matériau dans lequel les électrons
restent liés aux atomes ou molécules. Si on introduit un échantillon isolant entre les armatures d’un condensateur
chargé, l’effet du champ électrique extérieur appliqué consiste seulement à déformer les nuages électroniques et
donc à induire des dipôles : les électrons ne peuvent pas s’échapper des atomes/molécules. Il peut s’agir d’un
dipôle atomique ou moléculaire induit ou d’un dipôle moléculaire permanent qui s’oriente dans le champ externe
(molécule hétéropolaire). Selon le matériau, la polarisation électrique induite peut être plus ou moins grande
pour une même sollicitation électrique (i.e. pour un même champ externe appliqué) : les matériaux isolants ont
des susceptibilités différentes (à cause de mécanismes de polarisation microscopiques différents et aussi parce
que des atomes d’éléments chimiques différents vont se polariser différemment).
Comme pour les conducteurs, on ne peut pas décrire ou suivre les positions de toutes les charges, et les
champs microscopiques {~e, ~b} seraient très irréguliers avec des forces divergences au niveau des électrons et des
protons. On va donc adopter une approche de champ lissé à une échelle mésoscopique. Cela signifie que l’on
moyenne le champ microscopique sur une échelle spatiale grande par rapport à la taille des atomes mais petite
par rapport aux échelles macroscopiques. Le champ électrique macroscopique est défini comme
ZZZ
~ r) = h~eiv = 1
E(~ d~ρ ~e(~
ρ), (10.1)
v
où l’intégrale porte sur un volume v mésoscopique autour de ~r.

1 Polarisabilité d’un atome : apparition d’un dipôle induit

Considérons un atome isolé (par exemple au sein d’un gaz dilué pour bien pouvoir négliger l’influence des
atomes voisins). Nous devrions en principe décrire les électrons de l’atome par la mécanique quantique, mais
que peut dire la physique classique sur l’interaction d’un atome avec un champ électrique ? On peut utiliser
un modèle classique de l’atome. Considérons que le noyau est une charge +Q ponctuelle et que le "nuage
électronique" est décrit par une boule chargée négativement de manière uniforme. Appelons a le rayon de cette
boule, et −Q sa charge totale.
4πa3
−Q = ρ . (10.2)
3
La rigidité de l’atome est assurée par
Quand le champ électrique extérieur appliqué est nul E ~ = ~0, le centre des charge négatives N (centre de
la boule électronique) coïncide avec le centre des charges positives, le noyau P . Quand on applique un champ
non nul, le noyau positif et le nuage électronique sont attirés dans des directions opposées, ce qui engendre un
moment dipolaire électrique p~. Pour calculer ce dipôle, il faut évaluer le décalage N P . Pour cela, on va faire
l’hypothèse que le nuage reste sphérique ce qui permet de calculer le champ électrique créé par ce nuage en
utilisant le théorème de Gauss : −−→
~ N (M ) = ρN M .
E (10.3)
30
La force électrique attractive exercée par le nuage sur le noyau P est donc une sorte de force de rappel "élastique"
qui assure la cohésion de l’atome (dans ce modèle). Et par action-réaction, le nuage exerce la force opposée

62
sur le noyau. On impose alors l’équilibre mécanique du noyau P en écrivant que la force électrique extérieure
compense cette force de rappel (du au champ électrique interne du nuage) :

QE ~ ext = →
~ N (P ) + QE −
0 . (10.4)

Ce qui donne :
−−→
ρN P ~ ext = →
+E
−
0 , (10.5)
30
et finalement :
p~ = QE ~ ext .
~ N (P ) = 4πa3 0 E (10.6)

2 Vecteur polarisation
A l’échelle mésoscopique dV, on a beaucoup de dipôles atomiques dans la cellule et le moment dipolaire
électrique total est bien sûr proportionnel au volume de la cellule dV. On définit donc le vecteur polarisation,
comme la densité vectorielle de moment dipolaire électrique par unité de volume :
P
ip
~i
P~ (~r) = , (10.7)
dV

où la somme porte sur les dipôles p~i inclus dans dV, centré sur la position ~r.
Le champ macroscopique total E(~ ~ r) est la somme du champ extérieur appliqué E
~ a (~r) et du champ E
~ d (~r)
créé par la matière polarisée :

~ r) = E
E(~ ~ a (~r) + E
~ d (~r), (10.8)

où E~ d (~r) est aussi appelé le champ de dépolarisation (car il est souvent de sens opposé au champ appliqué, au
sein de l’échantillon isolant).

3 Charges de polarisation volumiques et surfaciques

On considère un échantillon isolant décrit par un champ de polarisation P~ (~r).

Exercice de cours : charges de polarisation équivalentes

1) Calculer le potentiel électrostatique créé par les dipôles contenus dans l’échantillon V.
2) Montrer que ce potentiel est équivalent au potentiel créé par des charges ponctuelles de densités
volumique ρ(~r) et surfacique σ(~r) données par :
→
−
ρ(~r) = − ∇.P~ (~r) et σ(~r) = P~ .~next , (10.9)

où la seconde formule est définie sur la surface extérieure du volume V.

3) Montrer que l’on peut déduire la forme de la densité surfacique à partir de la formule de la densité
→
−
volumique ρ(~r) = − ∇.P~ (~r).
4) Donner une interprétation simple de ces formules en termes de compensation des charges entre dipôles.

Dans les situations où la polarisation de l’isolant est uniforme, il existe une astuce très commode pour
~ d (~r). Elle consiste à passer par le calcul du potentiel créé par les dipôles :
calculer le champ dépolarisant E

P~ (r~0 )dr~0 .(~r − r~0 ) (~r − r~0 ) ~0

ZZZ ZZZ
1 1
Vd (~r) = 3 = P~ . 3 dr = P~ .E
~ aux , (10.10)
4π0 4π0

V ~0 V ~0
r − ~r r − ~r

63
où E~ aux est formellement le "champ électrique" qui correspondrait à une densité de charge volumique égale à 1
(attention il n’a donc pas la dimension d’un champ électrique).

Exercice : champ de dépolarisation dans une boule isolante uniformément polarisée

On applique un champ E~ a uniforme sur une boule isolante neutre. On suppose qu’elle se polarise uniformément
en volume avec un vecteur polarisation P~ .
1) Calculer le potentiel électrostatique créé par les dipôles contenus dans une boule de rayon R uniformé-
ment polarisée. On distinguera le potentiel à l’extérieur et à l’intérieur de la boule.
2) En déduire le champ créé par les dipôles contenus dans la boule (à l’intérieur et à l’extérieur de la
boule).
3) Ecrire le champ total.

4 Equation de Maxwell pour le champ lissé

L’idée est de séparer deux types de charges (les charges libres et les charges liées) en écrivant :
→
−
ρ(~r) = − ∇.P~ (~r) + ρlibre (~r), (10.11)

puis de ré-injecter cette expression dans l’équation de Maxwell-Gauss :

→
− ~ ~ r) = 0 E(~
~ r) + P~ (~r).
∇ · D(~r) = ρlibre (~r), avec D(~ (10.12)

~ r) qui ne "voit" que

A l’échelle mésoscopique, on peut travailler avec ce champ déplacement électrique D(~
les charges libres. C’est le champ que peut contrôler l’expérimentateur en ajustant les charges libres sur les
armatures des condensateurs par exemple. L’équation de Maxwell-structure reste inchangée pour les champs
lissés.
~ r).
Exercice : écrire l’équivalent du théorème de Gauss pour D(~

5 Notion de susceptibilité
On a "caché" les charges de polarisation dans un nouveau vecteur D(~ ~ r). Cependant pour décrire l’isolant
par un système fermé d’équations, il faut ajouter une relation entre la réponse de l’isolant (sa polarisation) et
le champ électrique. Cette relation peut être postulée de manière phénoménologique, ou dérivée à partir d’un
modèle microscopique.

5.1 Réponse linéaire et tenseur de susceptibilité électrique

Pour des champs assez faibles, cette relation est linéaire, et "le coefficient de proportionnalité" est en général
une matrice (tenseur de susceptibilité). Si le matériau est homogène, cette matrice ne dépend pas du point
considéré et on peut écrire :

P~ (~r) = 0 χE(~
~ r), (10.13)

où χ est une matrice (ou tenseur) s’appliquant sur les trois composantes du champ électrique.

5.2 Milieux linéaires homogènes et isotropes (LHI)

Dans les milieux isotropes, cette matrice est proportionnelle à l’identité et on a donc simplement :

P~ (~r) = 0 χE(~
~ r), (10.14)

où χ est un nombre (en général positif). De manière équivalente, on peut travailler avec le vecteur déplacement
électrique :

~ r) = 0 (1 + χ)E(~
D(~ ~ r) = 0 r E(~
~ r). (10.15)

64
Dans ces équations, on appelle susceptibilité le coefficient χ et constante diélectrique (ou permittivité relative)
le coefficient r . Ces coefficients sont sans unités (car on a mis en facteur 0 ). Ils nous renseignent sur la façon
dont le milieu matériel répond à une excitation électrique statique. On verra plus tard, qu’il existe une fonction
de réponse χ(ω) qui décrit comment la matière réagit à une excitation électrique de pulsation ω (onde EM par
exemple, ou encore un échantillon placé dans un condensateur soumis à une tension alternative).

6 Matière condensée
Dans la matière condensée, la densité des atomes ou molécules est très élevée ce qui rend subtil le lien entre
champ macroscopique lissé et champ réellement "ressenti" par un atome donnée (le champ local).

6.1 Notion de champ local

Lorsque l’on veut modéliser l’interaction de la matière avec un champ électrique, il faut évaluer la force
électrique subie par les électrons. Celle-ci est proportionnelle au champ microscopique local et non au champ
moyenné dans une cellule mésoscopique. Dans les isolants, il y a une différence entre le champ local et le champ
~
macroscopique lissé. Dans certaines conditions, Lorentz a établi la relation entre le champ lissé macroscopique E
et le champ ressenti microscopiquement par un dipôle atomique ou un électron E ~ local :

~
E ~+ P .
~ local = E (10.16)
30
Cette relation n’est pas universelle, car elle dépend du matériau considéré et de l’environnement microscopique
local des atomes. Pour obtenir cette relation, l’idée est de séparer le champ en deux contributions : le champ
créé par les charges lointaines pour lesquelles il est possible de moyenner, et une partie proche pour laquelle il
faut sommer individuellement les champs créés par chaque atome/molécule (dans un cristal à symétrie cubique
cette contribution est nulle ce qui conduit à la relation ci-dessus).

6.2 Lien entre polarisabilité et susceptibilité

La polarisabilité correspond au coefficient de proportionnalité α entre le dipôle microscopique et le champ
local excitateur :
~ local
p~ = 0 αE (10.17)

Il ne faut donc pas confondre la polarisabilité avec la susceptibilité qui exprime la proportionnalité entre le
vecteur polarisation et le champ électrique macroscopique. Notons que α est typiquement de l’ordre de grandeur
d’un volume microscopique, tandis que la susceptibilité est sans unités. Dans un milieu dilué (un gaz), il suffit
de multiplier par la densité d’atomes pour passer de χ à α. En revanche, dans un milieu dense (solide ou
liquide), il faut utiliser une formule de type Lorentz pour relier susceptibilité (macroscopique) et polarisabilité
(microscopique) :
~
P~ = N 0 αE ~ + NαP
~ local = N 0 αE ⇒ P~ (1 −
Nα ~
) = N 0 α E (10.18)
3 3
D’où l’on déduit la susceptibilité en fonction de la polarisabilité :

Nα
χ = r − 1 = . (10.19)
1 − N α/3

On peut inverser cette relation sous la forme suivante (relation de Clausius-Mossoti) :

Nα r − 1
= , (10.20)
3 r + 2

et on remarque que polarisabilité nulle (α = 0) correspond à une constante diélectrique relative r = 1. Notons
aussi que pour les milieux dilués, on retrouve bien par passage à la limite : χ = r − 1 = N α  1.

65
7 Boule isolante dans un champ électrique statique (équation de Laplace)
On considère donc une boule de centre O et de rayon R remplie d’un isolant polarisé uniformément. L’isolant
est supposé linéaire, homogène et isotrope (LHI) de susceptibilité χ. On applique un champ uniforme E ~ 0 = E0 ~uz
pour polariser l’isolant (prendre E0 > 0 pour fixer les idées et les dessins). Il est souvent très commode de
calculer les potentiels électriques plutôt que les champs. Les potentiels électriques obéissent à l’équation de
Poisson ou Laplace. On va donc chercher le potentiel électrique correspondant à une boule isolante immergée
dans un champ extérieur uniforme E ~ 0 en résolvant l’équation de Laplace avec les bonnes conditions aux limites.

1) Démontrer l’équation de Poisson dans le diélectrique LHI :

ρl
∆V + = 0, (10.21)

où ρl est la densité volumique de charges libres.
2) Faire la séparation des variables V (r, θ) = f (r)g(θ) en utilisant l’expression suivante du Laplacien :
   
1 ∂ 2 ∂ 1 ∂ ∂
∆= 2 r + 2 sin θ , (10.22)
r ∂r ∂r r sin θ ∂θ ∂θ

Justifier l’omission (volontaire) de la coordonnée ϕ.

3) Justifier que g(θ) = A cos θ.
4) Ecrire l’équation différentielle pour la partie radiale f (r). Chercher les solutions en loi du puissance du
type f (r) = rp . Quelles sont les valeurs de p possibles ?
5) Justifier, d’après ce qui précède, que le potentiel peut s’écrire de la manière suivante.
A l’intérieur de la sphère (r < R) :

V int (r, θ) = Ar−2 + Br cos θ .

(10.23)

A l’extérieur de la sphère (r > R) :

V ext (r, θ) = Cr−2 + Dr cos θ .

(10.24)

Il faut maintenant déterminer les coefficients inconnus A, B, C et D.

6) Il ne peut pas y avoir de singularité du potentiel en r = 0. En déduire une condition sur un des coefficients
inconnus.
7) Quand r → ∞, l’influence de la sphère doit disparaître et on doit retrouver le potentiel correspondant au
~ 0 = E0 ~uz . En déduire la valeur de D.
champ uniforme E
8) Exprimer les conditions de passage en r = R qui sont d’une part la continuité du potentiel V , et d’autre
~ (au passage justifier pourquoi cette composante
part la continuité de la composante normale du vecteur D
normale est continue alors que la composante normale du champ électrique ne l’est pas). Résoudre le
système de deux équations obtenu pour les coefficients B et C.
9) En déduire le champ électrique à l’intérieur de la boule.

66
Chapitre 11

Champ magnétique dans la matière

On pense à un aimant, mais en fait la plupart des matériaux ont des propriétés magnétiques plus ou
moins marquées. Certains matériaux sont paramagnétiques : l’aimantation induite est orientée dans le même
sens que le champ magnétique extérieur appliqué. D’autres matériaux sont diamagnétiques : une très faible
aimantation induite se développe dans le sens opposé au champ magnétique appliqué. Il existe aussi des matériaux
ferromagnétiques qui possèdent une aimantation spontanée en l’absence de champ magnétique appliqué (en
dessous d’une certaine température Tc ). Les supraconducteurs ont aussi des propriétés magnétiques spectaculaires
comme l’effet Meissner : écrantage complet d’un champ magnétique extérieur appliqué (diamagnétisme parfait).

1 Dipôles magnétiques et matériaux magnétiques

Rappeler le dipôle magnétostatique. Présenter brièvement l’origine du moment magnétique atomique. Décrire
la phénoménologie des matériaux paramagnétiques, diamagnétiques, ferromagnétiques. Pour étudier la réponse
magnétique d’un matériau, l’expérience typique consiste à introduire un échantillon dans une bobine (ou un
électroaimant) qui créé le champ appliqué.

1.1 Matériaux paramagnétiques

Moments magnétiques majoritairement orientés dans la même direction que le champ magnétique appliqué.

1.2 Matériaux diamagnétiques

Moments magnétiques majoritairement orientés dans la direction opposée que le champ magnétique appliqué.

1.3 Matériaux ferromagnétiques

Les moments magnétiques s’alignent spontanément dans une certaine direction en l’absence de champ
magnétique appliqué.

2 Vecteur aimantation
A l’échelle mésoscopique dV, on a beaucoup de dipôles atomiques dans la cellule et le moment dipolaire
magnétique total est bien sûr proportionnel au volume de la cellule dV. On définit donc le vecteur aimantation,
comme la densité vectorielle de moment dipolaire magnétique par unité de volume :
P
~ (~r) = im
~i
M , (11.1)
dV

où la somme porte sur les m

~ i inclus dans dV, centré sur la position ~r.

67
3 Courants ampériens volumiques et surfaciques
~ (~r).
On considère un échantillon décrit par un champ d’aimantation M

Exercice de cours : courants d’aimantation équivalents

1) Calculer le potentiel vecteur créé par les dipôles magnétiques contenus dans l’échantillon V.
2) Montrer que ce potentiel vecteur est équivalent au potentiel que créeraient des courants volumiques :

~j(~r) = →
− ~ (~r),
∇ ∧M (11.2)

et des courants surfaciques sur la frontière de V dont on explicitera la densité surfacique.

3) Montrer que l’on peut déduire la forme de la densité surfacique à partir de la formule de la densité
→
− ~ (~r).
volumique de courants ∇ ∧ M
3) Donner une interprétation simple de ces formules en termes de compensation entre les "boucles de
courants microscopiques".

4 Equations de Maxwell pour les champs lissés dans la matière

Comme pour les diélectriques, l’idée est de séparer deux types de sources. Ici ce sont des courants : les
courants libres et les courants ampériens correspondant à la réponse de la matière. On écrit ainsi :

~j(~r) = ~jlibre (~r) + →

− ~ (~r),
∇ ∧M (11.3)

puis on ré-injecte cette expression dans l’équation de Maxwell-Ampère microscopique :

→
− ~ r) = ~jlibre (~r), ~ r) = B(~
~ r ) − µ0 M
~ (~r).
∇ ∧ H(~ avec µ0 H(~ (11.4)

~ r) qui ne "voit" que

A l’échelle mésoscopique, on peut travailler avec ce champ "excitation magnétique" H(~
les courants libres. C’est le champ que peut contrôler l’expérimentateur en réglant les courants libres dans les
bobines par exemple. En exercice, écrire le théorème d’Ampère pour H(~~ r).

5 Susceptibilité magnétique
Mêmes idées que pour la susceptibilité électrique, sauf qu’ici on regarde la réponse magnétique en terme
~ r) :
d’aimantation à une sollicitation magnétique décrite par le vecteur excitation magnétique H(~

~ (~r) = χm H(~
M ~ r), (11.5)

où χm est une matrice (ou tenseur) s’appliquant sur les trois composantes du champ excitation magnétique.
Dans un milieu LHI (linéaire homogène isotrope), la susceptibilité magnétique est simplement un nombre qui
peut être soit positif (paramagnétisme) soit négatif (diamagnétisme).
~ H}
On peut décrire la réponse magnétique de la matière en utilisant le couple {B, ~ :

~ r) = µ0 (H(~
B(~ ~ r) + M
~ (~r)) = µ0 (1 + χm )H(~
~ r) = µ0 µr H(~
~ r), (11.6)

ce qui permet de définir au passage la perméabilité magnétique relative du matériau comme : µr = 1 + χm .

6 Boule aimantée dans un champ magnétique statique

68
Chapitre 12

Ondes dans les milieux matériels

Cette partie est consacrée au comportement des ondes électromagnétiques dans les milieux matériels dilués
(vapeurs atomiques) ou denses (isolants ou conducteurs). L’objectif est de relier les phénomènes optiques
(dispersion, réflection, absorption) et les grandeurs qui les caractérisent (indice optique, coefficient de réflexion
ou d’absorption) à des modèles microscopiques de l’interaction entre la matière condensée et la lumière.

1 Equations de Maxwell dans la matière

On peut écrire les équations de Maxwell complètes en régime variable pour les champs macroscopiques.

1.1 ~ et H
Utilisation des vecteurs D ~ et relations constitutives

L’idée est d’utiliser les vecteurs déplacement électrique et excitation magnétique pour ne faire apparaître que
les charges et les courants libres éventuels. Ainsi seules les 2 équations de Maxwell présentant des sources sont à
reformuler :

→
− ~
∇ · B(~r, t) = 0, (12.1)
→
− ~ ~ r, t)
∂ B(~
∇ ∧ E(~r, t) = − , (12.2)
∂t
→
− ~
∇ · D(~r, t) = ρlibre (~r, t), (12.3)
→
− ~
~ r, t) = ~jlibre (~r, t) + ∂ D(~r, t) .
∇ ∧ H(~ (12.4)
∂t

On remarque que 0 = 8, 85.10−12 F.m−1 et µ0 = 4π.10−7 H.m−1 ont disparu, mais c’est juste le résultat du jeu
de réécriture des champs.

Des équations de Maxwell vers les équations d’onde dans la matière

1) En utilisant les équations de Maxwell ci-dessus et les relations constitutives D(~ ~ r) = 0 r E(~
~ r) et
~ ~ ~
B(~r) = µ0 µr H(~r), obtenir des équations découplées pour le champ D(~r, t) seul.
~ r, t).
2) Obtenir l’équation similaire pour le champ excitation magnétique H(~
3) Que deviennent ces équations découplées dans une région de l’espace vide de charges et de courants
libres ? Quelle est la conséquence très importante de ces équations ? Quelle notion importante apparaît
ici ?

1.2 ~ B
Utilisation des vecteurs E, ~ et P~ et M
~
~ et B
Dans certains contextes, il peut être intéressant de privilégier les vecteurs E ~ en faisant apparaître
→
− ~
explicitement les charges et les courants liées ou ampériens sous leur forme ρ(~r) = − ∇.P (~r), ~jd = ∂ P~ /∂t, et
~jA (~r) = →
−
∇ ∧M~ (~r). On obtient les équations de Maxwell sous la forme :

69
 →
− ~
∇ · B(~r, t) = 0, (12.5)
→
− ~ ~ r, t)
∂ B(~
∇ ∧ E(~r, t) = − , (12.6)
∂t
→
− ~ 1 1 →− ~
∇ · E(~r, t) = ρlibre (~r, t) − ∇ · P (~r, t), (12.7)
0 0
→
− ~ ∂ P~ (~r, t) →
− ~
~ (~r, t) + µ0 0 ∂ E(~r, t) .
∇ ∧ B(~r, t) = µ0~jlibre (~r, t) + µ0 + µ0 ∇ ∧ M (12.8)
∂t ∂t

2 Modèle microscopique de l’interaction lumière-matière : isolants

Lorentz a formulé une théorie de l’interaction onde électromagnétique-matière entièrement basée sur la
physique classique : la matière y est décrite par des oscillateurs harmoniques et le rayonnement électromagnétique
par les équations de Maxwell. La matière est représentée comme un milieu contenant des oscillateurs pouvant
avoir des fréquences propres très différentes. A l’époque de Lorentz (fin du 19-ième siècle) on ne connaissait pas
l’origine de ces dipôles. Maintenant on peut les associer à différents types de degrés de libertés microscopiques :
les dipôles atomiques correspondant aux électrons liés autour des noyaux (fréquences typiques dans l’UV), les
dipôles correspondant aux vibrations des atomes (typiquement dans IR), et les dipôles associés aux électrons
libres (pas de fréquence correspondant à une force de rappel mais une fréquence plasma qui dépend de la densité
électronique). Nous supposons pour commencer que l’on peut isoler un type d’oscillateurs caractérisé par une
fréquence propre bien définie ω0 (et donc ignorer les oscillateurs ayant d’autres fréquences).

2.1 Modèle de l’électron élastiquement lié

Pour simplifier, on commence par considérer de la matière diluée, comme un gaz atomique par exemple, et
que ω0 correspond à une raie de transition bien définie de cet atome. On considère un électron lié à un noyau
atomique par une force de rappel de type élastique et soumis à la force électrique due à l’onde électromagnétique.
La dynamique est décrite par l’équation suivante :

d2~r d~r
m 2
+ mγ + mω02~r = −eE0 cos(ωt)~u, (12.9)
dt dt
où −mω02~r représente la force de rappel élastique, −eE0 cos(ωt)~u la force due au champ électrique de l’onde (à
polarisation rectiligne), et −mγd~r/dt une force d’amortissement. Le coefficient γ = 1/τ est un taux de relaxation
décrivant le retour à l’équilibre si on débranche le champ électrique. C’est le champ macroscopique qui apparaît
car dans un milieu dilué il n’y a pas de différence importante entre ce champ et le champ local ressenti par
l’électron.
En injectant une solution forcée ~r(t) = r0 e−iωt ~u, et on obtient pour le dipôle microscopique p~(t) = −e~r(t) :

e2 E0 e−iωt ~
p~ = ~u = 0 α(ω)E. (12.10)
m(−ω 2 − iωγ + ω02 )

La polarisabilité atomique α(ω) (définie ci-dessus) s’exprime comme :

e2
α(ω) = , (12.11)
m0 (−ω 2 − iωγ + ω02 )

et son unité est celle d’un volume, souvent de l’ordre du volume microscopique du dipôle (typiquement en nm3 ).
Dans un milieu homogène, on passe au vecteur polarisation macroscopique (densité volumique de moment
dipolaire électrique) en multipliant le dipôle élémentaire par la densité volumique d’atomes N dans le milieu :

N e2 E0 e−iωt
P~ = N p~ = ~u (12.12)
m(−ω 2 − iωγ + ω02 )

On en déduit la susceptibilité diélectrique complexe (sans unités) :

N e2 1
χ(ω) = N α(ω) = 2 , (12.13)
m0 (ω0 − ω 2 − iωγ)

70
puis immédiatement la constante diélectrique complexe (permittivité diélectrique relative) :
N e2 1
r (ω) = 1 + χ(ω) = 1 + . (12.14)
m0 (ω02 − ω 2 − iωγ)
On voit apparaître une fréquence typique : s
N e2
ωP = , (12.15)
m0
qui est caractéristique de la densité de dipôles N et que l’on appelle la fréquence plasma. Cette fréquence
jouera un rôle très important pour les métaux (pour lesquels il n’y a pas d’autre fréquence caractéristique :
ω0 = 0). On peut alors écrire la constante diélectrique sous la forme :
ωP2
r (ω) = 1 + , (12.16)
ω02 − ω 2 − iωγ
mettant en jeu les trois fréquences caractéristiques : le fréquence de résonance ω0 , la fréquence plasma ωP , et le
taux d’amortissement γ.

3 Modèle microscopique de l’interaction lumière-matière : métaux

On s’intéresse maintenant à la réponse électromagnétique d’un milieu conducteur (métal, plasma) qui contient
une certaine densité de porteurs de charge libres de se déplacer sous l’influence des forces électromagnétiques.
On se limite à une seule espèce de porteurs, par exemple le gaz d’électrons d’un métal. Dans un plasma, il y a
en général au moins deux types de porteurs, typiquement les électrons et les ions positivement chargés.

3.1 Modèle de Drude AC et conductivité optique

On reprend le modèle de Drude :
d~v m
m = qE0 cos(ωt)~u − ~v , (12.17)
dt τ
et on va l’analyser cette fois-ci en régime haute-fréquence. La force exercée par l’onde électromagnétique sur
l’électron, qE0 cos(ωt)~u, se limite à la force électrique car la force de Lorentz magnétique est négligeable, le
champ magnétique de l’onde étant en B = E/c. Ce modèle est formellement le même que le modèle de Lorentz
dans lequel on a enlevé le terme de force de rappel élastique : ω0 = 0. C’est normal car on veut décrire des
électrons libres d’un métal, et non les électrons liés d’un isolant. Notre analyse du modèle de Lorentz reposait sur
la notion de dipôle électrique et de vecteur polarisation. Pour un métal, les électrons sont fortement délocalisés,
et la réponse en terme de courants électriques semble plus naturelle que celle de dipôle.
On va donc écrire les équations pour le courant ~j(t) = N q~v (t),

d~j 1 N q 2 E0
+ ~j = cos(ωt)~u, (12.18)
dt τ m
et on passe en notation complexe :
N q2τ ~
(−iωτ + 1) ~jω = Eω , (12.19)
m
D’où finalement une relation constitutive :
~jω = σ(ω)E ~ω, (12.20)
de type réponse linéaire avec une conductivité (AC ou optique) dépendant de la fréquence comme :
N q 2 τ /m σ0
σ(ω) = = . (12.21)
(1 − iωτ ) (1 − iωτ )
Cette équation constitutive doit être complété par le fait que dans un métal, l’électroneutralité est rétablie
très rapidement devant toutes les autres échelles temporelles, ce qui permet de poser ρ = 0. Pour savoir comment
une onde peut se propager dans un métal il suffit d’écrire les équations de Maxwell complétées par ces relations
constitutives. En particulier en évaluant le rotationnel du rotationnel du champ électrique, nous obtenons une
équation d’évolution pour le champ électrique :
 
2 2~ 2 iωτ 2 ~ω ,
c k Eω = ω + ω E (12.22)
1 − iωτ P

71
dans laquelle apparaît la pulsation plasma :

N q2 σ0
ωP2 = = . (12.23)
m0 0 τ

Ainsi on a vu apparaître deux grandeurs, la conductivité statique σ0 et la pulsation plasma ωP , qui sont deux
paramètres permettant de quantifier la densité N de porteurs de charges libres dans le système, ainsi que les
caractéristiques de ces porteurs (via leur charge q et leur masse m). On peut utiliser l’un ou l’autre selon le
contexte.
On doit analyser les différents régimes possibles.

Regime AC basse fréquence : ω  1/τ

On peut négliger ωτ  1 au dénominateur mais il faut le conserver au numérateur :

σ0
c2 k 2 = ω 2 + iωτ ωP2 ' iωτ ωP2 = iω


, (12.24)
0

car le terme en ω domine le terme en ω 2 à basse fréquence. On donc une relation de dispersion qui conduit à un
vecteur d’onde k ayant une partie imaginaire :
r
σ0 ω
k = (1 + i) . (12.25)
2c2 0

La longueur de pénétration :
s r
2c2 0 2
δ= = , (12.26)
σ0 ω σ0 µ0 ω

tend vers zéro pour un conducteur idéal σ0 → ∞.

Regime de relaxation : 1/τ  ω < ωP

Dans la limite ωτ  1, on a la relation de dispersion :

c2 k 2 = ω 2 − ωP2 < 0 ,


(12.27)

d’où un vecteur d’onde imaginaire pur.

Regime haute fréquence (transparent) : 1/τ  ωP < ω

Au-dessus de la pulsation plasma, on a

c2 k 2 = ω 2 − ωP2 > 0 ,


(12.28)

et donc le vecteur d’onde redevient réel et il y a propagation : le gaz d’électrons de conduction devient transparent
à l’onde électromagnétique.

Conclusion, remarques :

On voit que la dispersion résulte d’une compétition entre le terme en ω 2 qui vient du courant de déplacement
toujours présent (y compris dans le vide) et qui favorise la propagation, et entre le terme en −ωP2 qui est lié à la
matière (proportionnel à N ). On remarque que le terme en −ωP2 porte un signe opposé au terme de courant de
déplacement : il a donc tendance à empêcher la propagation et peut même la supprimer complètement lorsque
la densité de matière est assez forte : ωP > ω.
On a envisagé uniquement ce que l’on appelle les modes transverses électriques puisque l’on a fait dès
~ = 0 et donc ~k ⊥ E.
le départ l’hypothèse ρ = 0, ce qui revient à div E ~ Il existe aussi un mode longitudinal dans
~
lequel le champ électrique est parallèle au vecteur k et qui correspond à des oscillations de la densité de charge.

72
3.2 Lien avec le modèle de Lorentz (introduit pour les isolants)
Nous avons vu que les métaux contiennent des électrons libres de se déplacer dans tout l’échantillon. Ces
électrons libres sont responsables du phénomène de conduction électrique et dominent la réponse du métal à une
onde électromagnétique. Nous avons écrit une relation constitutive exprimant la réponse linéaire en courant à
un champ électrique. Le coefficient de proportionnalité étant la conductivité (ac ou optique) σ(ω). Cependant,
ces électrons de conduction peuvent aussi être décrits par le modèle de Lorentz à condition d’éliminer la force
de rappel, c’est à dire de prendre ω0 = 0. Dans la théorie classique, il y a un caractère un peu arbitraire de la
séparation entre électrons liées et électrons de conductions dans la théorie classique.
Un métal solide (ou liquide) est aussi un milieu dense, cependant les électrons de conduction forment un
fluide impliquant que le champ local ressenti par un électron est égal au champ lissé macroscopique. Donc
contrairement au cas des isolants, nous ne devons pas faire la correction de type Clausius-Mossoti et on obtient
simplement :

N e2 1 ωP2
r (ω) = n(ω)2 = 1 − = 1 − (12.29)
m0 (ω 2 + iωγ) (ω 2 + iωγ)
Pour les fréquence au-dessous de la la fréquence plasma, l’indice est imaginaire pur et le métal sera parfaitement
réfléchissant. Pour les fréquences au-dessus de la fréquence plasma, l’indice devient réel et le métal (ou plutôt le
gaz d’électrons de conduction) devient parfaitement transparent.
Si on néglige l’amortissement, on obtient : Pour un amortissement nul γ = 0 (ou très faible), on obtient :

N e2 1 ωP2
r (ω) = n(ω)2 = 1 − = 1 − (12.30)
m0 ω 2 ω2
En fait, il existe une relation entre la constante diélectrique (ou la susceptibilité χ(ω)) et la conductivité
σ(ω). En effet on peut écrire l’équation de Maxwell-Ampère en considérant que les courants sont de la forme σ E ~
ou bien ∂ P~ /∂t, on a donc la relation :
iσ(ω)
χ(ω) = , (12.31)
0 ω
qui relie susceptibilité et conductivité d’un même matériau.
Vérifions cette relation sur l’exemple du gaz électronique d’un métal pour lequel nous avons obtenu d’une
part (modèle de Drude) :
N q 2 τ /m σ0
σ(ω) = = . (12.32)
(1 − iωτ ) (1 − iωτ )
et d’autre part (modèle de Lorentz avec ω0 = 0) :

N e2 1 ωP2
χ(ω) = r (ω) − 1 = − = (12.33)
m0 (ω 2 + iωγ) (−ω 2 − iωγ)

4 Fonction de réponse
On a vu apparaître la notion de coefficient de réponse de type susceptibilité reliant la réponse d’un milieu
(sa polarisation ou son aimantation) à l’excitation. Dans ce paragraphe, on discute plus en détail cette notion
très générale dans le cadre de la théorie de la réponse linéaire.

4.1 Origine de la dispersion en pulsation ω

En présence d’une sollicitation variable dans le temps, la matière répond avec un certain retard temporel.
On peut donc écrire la polarisation à la date t comme une intégrale sur les valeurs du champ excitateur aux
instants antérieurs à t : Z ∞
~
P (t) = ~ − τ ),
dτ χ̃(τ )E(t (12.34)
0
où χ(τ ) est une fonction décroissante de τ (attention cette fonction porte des unités contrairement à la
susceptibilité que nous avons utilisée jusqu’ici). La variable d’intégration τ représente le délai temporel entre
l’excitation et la réponse. Cette relation intégrale est tout à fait linéaire mais elle n’est plus instantanée. Elle se
simplifie si on passe dans l’espace de Fourier en écrivant les décompositions en modes de Fourier harmoniques :
Z ∞ Z ∞
dω ~ −iωt dω ~
~
E(t) = E(ω)e , et ~
P (t) = P (ω)e−iωt . (12.35)
−∞ 2π −∞ 2π

73
On obtient :
Z ∞ Z ∞ Z ∞ Z ∞ 
dω ~ dω
P~ (t) = dτ χ̃(τ ) E(ω)e−iω(t−τ ) = iωτ
dτ χ̃(τ )e ~
E(ω)e−iωt
, (12.36)
0 −∞ 2π −∞ 2π 0

et on peut ainsi exprimer l’amplitude du mode ω de la polarisation comme :

Z ∞
~ ~
P (ω) = χ(ω)E(ω) avec χ(ω) = dτ χ̃(τ )eiωτ . (12.37)
0

Cette dépendance en ω de la susceptibilité est appelée la dispersion temporelle. On vient de voir qu’elle est
intimement liée au retard de réponse de la matière, c’est à dire aux processus microscopiques de réponse et
de relaxation. Notons aussi que χ(ω) est la transformée de Fourier de χ̃(τ ) : ce sont des fonctions différentes
(d’ailleurs l’une est sans unités et l’autre est homogène à l’inverse d’un temps). Si la réponse de la matière était
instantanée χ̃(τ ) = 2χ0 δ(τ ), on aurait alors χ(ω) = χ0 indépendante de ω (cela peut arriver dans les situations
où on peut négliger les effets de retard). Enfin, une conséquence importante de la dispersion temporelle de χ(ω)
va être le fait que l’indice optique n(ω) va lui même dépendre de la fréquence.

4.2 Signification du caractère complexe de la fonction de réponse

On peut décomposer la susceptibilité électrique χ(ω) = χ1 (ω) + iχ2 (ω) en parties réelle χ1 et imaginaire χ2 .
La fonction de réponse χ(ω) est "le coefficient de proportionnalité" entre les représentations complexes de deux
grandeurs physiques réelles, donc les parties réelle χ1 et imaginaire χ2 ont une signification physique importante
(contrairement à la représentation complexe ou de Fresnel d’une grandeur oscillante dont celle la partie réelle est
physique, et la partie imaginaire est redondante).
~
Prenons le champ électrique E(t) = E~ 0 cos(ωt) avec E~ 0 réel, et écrivons la réponse, i.e. la polarisation
instantanée du milieu P~ (t) = Re[P~ (ω)e−iωt ], comme :

~ 0 e−iωt ] = 0 Re[χ(ω)e−iωt ]E
P~ (t) = 0 Re[χ(ω)E ~ 0 = 0 (χ1 (ω) cos(ωt) + χ2 (ω) sin(ωt)) E
~0 (12.38)

Il est clair que la partie réelle χ1 (ω) quantifie la réponse en phase avec l’excitation électromagnétique. Dans
le modèle de l’électron élastiquement lié, la fonction χ1 (ω) présente une forte discontinuité à la fréquence propre.
En revanche, la partie imaginaire χ2 (ω) décrit la réponse en quadrature de phase avec l’excitation. Dans
le modèle de Lorentz, χ2 (ω) présente un pic à la fréquence propre, dont la largeur est donnée par le taux
d’amortissement γ = 1/τ . Ceci est généralisable à n’importe quelle fonction de réponse linéaire caractérisant
le rapport d’amplitude et de déphasage de deux grandeurs oscillantes à la même fréquence : impédance d’un
élément de circuit, fonction de transfert d’une portion de circuit électrique, polarisabilité, constante diélectrique,
conductivité optique, etc...

4.3 Analyse énergétique : parties dissipative et réactive de la fonction de réponse

Selon le type de fonction de réponse considérée....conductivité susceptibilité....Dans ce paragraphe, on
sous-entend la dépendance en pulsation des fonctions de réponse mais on écrira χ1 = χ1 (ω) et χ2 = χ2 (ω) pour
ne pas surcharger.
Dans un milieu isolant, ou plus généralement si l’on décide de décrire les courants par le terme ∂P/∂t, la
densité volumique d’énergie transférée par unité de temps du champ électromagnétique vers la matière s’écrit en
instantané à la date t :

∂ P~ (t) ~ ~
.E(t) = 0 ω −χ1 sin(ωt) cos(ωt) + χ2 cos2 (ωt) E

2
0 (12.39)
∂t
Lorsqu’on prend la moyenne temporelle sur une période T = 2π/ω, seule la partie en cos2 (ωt) contribue et :

∂ P~ (t) ~ ~2
0 E0
h .E(t)i = ωχ2 , (12.40)
∂t 2
qui a bien la dimension d’une densité volumique de puissante. On voit que c’est la partie imaginaire de la
susceptibilité qui correspond à l’échange d’énergie de l’onde vers la matière. Dans le modèle de Lorentz, χ2 (ω)
présente un pic à la pulsation propre qui traduit donc une absorption efficace de l’onde par l’oscillateur lorsqu’il
y a résonance.

74
 Dans un conducteur, ou plus généralement, si l’on décide de décrire les courants par une relation d’Ohm
~ la densité volumique de puissance cédée par le champ électromagnétique à la matière s’écrit en
locale ~j = σ E,
instantané à la date t :
~
~j(t).E(t) ~2
= (χ1 cos(ωt) + χ2 sin(ωt)) E (12.41)
0

Lorsqu’on prend la moyenne temporelle sur une période T = 2π/ω, seule la partie en cos2 (ωt) contribue et on a
donc :
~2
σ1 E
~
h~j(t).E(t)i = 0
, (12.42)
2
On rappelle qu’en théorie classique, il existe une sorte d’ambiguité ou de liberté sur la représentation des
courants conduisant à deux types de description pour un même système : soit en terme de susceptibilité soit en
termes de conductivité. La relation générale étant :

iσ(ω)
χ(ω) = , (12.43)
0 ω
on voit que les deux formules obtenues sont parfaitement compatibles et donnent la même puissance volumique
moyenne, puisque χ2 = σ1 /0 ω.

75
Chapitre 13

Phénomènes aux interfaces : réflexion,

transmission et confinement

1 Comportement des champs EM à une interface

On va établir les relations de passage aux interfaces des différents champs. La méthode est très générale et
pourra être utilisée pour d’autres EDP dans d’autres contextes (équation de la chaleur, équation de Schrödinger,...)
Puisque les équations de Maxwell sont des EDP d’ordre 1 par rapport aux variables d’espace, les composantes
des champs électrique et magnétique doivent être des fonctions continues de l’espace. Cependant à l’interface
entre deux milieux, ils peuvent varier très brutalement et donc "apparaître discontinus" à l’échelle d’observation
macroscopique.
Dans toute cette section on prendra une interface plane, d’équation z = 0, séparant deux milieux différents.
La méthode présentée est décrite en détail dans Feynman (tome 2 d’électromagnétisme, chapitre 33), et présente
l’avantage de pouvoir être utilisée dans d’autres contextes.

1.1 Analyse des équations de Maxwell "microscopiques" à divergence

Prenons l’équation de Maxwell-Gauss que l’on écrit de manière explicite :

∂Ex (x, y, z) ∂Ey (x, y, z) ∂Ez (x, y, z) ρ(x, y, z)

+ + = . (13.1)
∂x ∂y ∂z 0

Intégrons cette équation dans une fine couche autour de l’interface, entre z = −δ et z = δ :
Z z=δ Z z=δ Z z=δ Z z=δ
∂Ex (x, y, z) ∂Ey (x, y, z) ∂Ez (x, y, z) ρ(x, y, z)dz
dz + dz + dz = . (13.2)
z=−δ ∂x z=−δ ∂y z=−δ ∂z z=−δ 0

Dans le passage à la limite δ → 0, les deux premiers termes à gauche tendent vers zero, le troisième terme
s’intègre comme le "saut de Ez " à l’interface et le terme de droite ne peut être non nul qu’en présence d’une
accumulation surfacique de charges en z = 0 :

σ(x, y)
Ez (x, y, 0+ ) − Ez (x, y, 0− ) = , (13.3)
0

où σ désigne la densité surfacique de charges accumulée sur l’interface (en z = 0)en C.m−2 . La composante
normale à l’interface du champ électrique est donc continue en z = 0 en l’absence d’accumulation de charges.
~ est nulle) s’analyse de manière tout à fait
L’équation de "Maxwell-flux conservatif" (la divergence de B
similaire et conduit à la conclusion que la composante normale à l’interface du champ magnétique est
toujours continue à cause de l’absence de terme source dans l’équation de "Maxwell-flux conservatif".

Bz (x, y, 0+ ) − Bz (x, y, 0− ) = 0 . (13.4)

76
1.2 Analyse des équations de Maxwell "microscopiques" à rotationnel
Prenons l’équation de Maxwell-Ampère. S’agissant d’une équation vectorielle, nous avons trois équations
scalaires à analyser simultanément. Seules deux équations contiennent les dérivées "intéressantes" ∂/∂z, c’est à
dire celles qui peuvent prendre des valeurs très importantes au franchissement de la frontière z = 0 et engendrer
des discontinuités à l’échelle macroscopique :
∂Bx (x, y, z) ∂Bz (x, y, z) ∂Ey (x, y, z)
− = µ0 jy + µ0 0 (13.5)
∂z ∂x ∂t
∂Bz (x, y, z) ∂By (x, y, z) ∂Ex (x, y, z)
− = µ0 jx + µ0 0 (13.6)
∂y ∂z ∂t
En présence d’une interface en z = 0, seuls les termes en ∂/∂z peuvent conduirent à des valeurs très grandes et
on peut donc écrire "symboliquement" au sens du passage à la limite δ → 0 :
Z z=δ Z z=δ
∂Bx (x, y, z)
dz = µ0 jy dz , (13.7)
z=−δ ∂z z=−δ
Z z=δ Z z=δ
∂By (x, y, z)
dz = − µ0 jx dz . (13.8)
z=−δ ∂z z=−δ

D’où :

Bx (x, y, 0+ ) − Bx (x, y, 0− ) = µ0 jSy , (13.9)

+ −
By (x, y, 0 ) − By (x, y, 0 ) = −µ0 jSx . (13.10)

La composante tangentielle (dans le plan de l’interface Oxy) du champ magnétique peut donc présenter une
discontinuité en présence d’une accumulation de courants surfaciques s’écoulant dans le plan Oxy.
L’équation de Maxwell-Faraday s’analyse de manière tout à fait similaire et conduit à la conclusion que la
composante tangentielle à l’interface du champ électrique est toujours continue à cause de l’absence
de terme source :

Ex (x, y, 0+ ) − Ex (x, y, 0− ) = 0 , (13.11)

+ −
Ey (x, y, 0 ) − Ey (x, y, 0 ) = 0 . (13.12)

1.3 Analyse des équations de Maxwell pour les champs macroscopiques lissés
Les équations de Maxwell dans la matière ont la même structure que les équations analysées ci-dessus à la
différence du fait que les termes sources impliquent seulement les charges libres et les courants libres, et que
certains coefficients numériques différent (les µ0 et 0 ). On n’a donc pas à refaire tout le travail et on trouve
tout de suite que :

Dz (x, y, 0+ ) − Dz (x, y, 0− ) = σl (x, y) , (13.13)

où σl est la densité surfacique de charges libres sur l’interface (en z = 0)en C.m−2 .

Hx (x, y, 0+ ) − Hx (x, y, 0− ) = jSly , (13.14)

+ −
Hy (x, y, 0 ) − Hy (x, y, 0 ) = −jSlx , (13.15)

où jSlx et jSyx sont les composantes des courants surfaciques libres en A.m−1 .

2 Problème de diffusion (scattering) à une interface

Dans un problème de diffusion (au sens de scattering en anglais) on a typiquement une onde (ou une ou
plusieurs particules) qui arrive de l’infini et qui s’approche d’une zone d’interaction pour donner naissance à une
ou plusieurs ondes (ou particules) diffusées qui s’éloignent de la zone d’interaction. Ici l’objet diffuseur est une
interface entre deux milieux.

77
2.1 Généralités
On considère deux milieux, notés 1 et 2, occupant respectivement les demi-espaces z > 0 et z < 0. On veut
savoir comment une onde incidente de type OPPM donc de vecteur d’onde ~ki et de pulsation ω bien définis
venant du milieu 1 va réagir à la présence d’une discontinuité de milieu en z = 0. L’interface va donner lieu à
une onde réfléchie dans le milieu 1, caractérisée par {~kr , ω} et à une onde transmise dans le milieu 2, caractérisée
par {~kt , ω}.
Conservation de la pulsation. Tout d’abord, il faut remarquer que toutes les ondes oscillent temporellement
à la même pulsation ω. Ceci est du à la linéarité des équations de Maxwell et au fait que l’interface est statique.
Au niveau plus microscopique, on peut dire que l’onde incidente provoque un mouvement forcé (accéléré) des
charges à la pulsation ω. Ces dernières rayonnent donc des ondes à la même pulsation. Mathématiquement, les
équations de passage des champs en z = 0 doivent être vérifiées pour toute date t ce qui implique ce synchronisme
des différentes ondes incidente, réfléchie et transmise.
Conservation du vecteur d’onde tangentiel. Par ailleurs, l’objet diffuseur a une géométrie très particu-
lière : c’est un plan (z = 0) et il y a donc invariance par translation dans ce plan. Ceci implique que les ondes
incidente, transmise et réfléchie ont la même composante tangentielle du vecteur d’onde. Cette conservation
de ~k tangentiel en conjonction avec les relations de dispersions dans les milieux 1 et 2, fixe la direction des
vecteurs d’onde réfléchi et transmis par rapport à la direction du vecteur d’onde incident : ce sont les lois de
Snell-Descartes.
Calcul des coefficients de transmission et réflexion. On peut calculer les amplitudes des ondes réfléchie
et transmise. Celles-ci sont proportionnelles à l’amplitude de l’onde incidente. On définit donc des coefficients de
réflexion r et de transmission t. Pour cela on écrit la structure des ondes incidente, réfléchie et transmise. Pour
les champs électriques on a

E ~ i0 ei(~ki .~r−ωt) ,
~ i (~r, t) = E (13.16)
E ~ i0 ei(~kr .~r−ωt) ,
~ r (~r, t) = rE (13.17)
E ~ i0 ei(~kt .~r−ωt) ,
~ t (~r, t) = tE (13.18)
et pour les champs magnétiques, on utilise la relation issue de Maxwell-Faraday :
→
−
B~ i (~r, t) = ki ∧ E
~ i0 ei(~ki .~r−ωt) , (13.19)
ω
→
−
~ kr ~ i(~kr .~r−ωt)
Br (~r, t) = r ∧ E i0 e , (13.20)
ω
→
−
~ kt ~ i(~kt .~r−ωt)
Bt (~r, t) = t ∧ E i0 e , (13.21)
ω
Ensuite on exprime les relations de passage pour les composantes tangentielles et normales des divers champs.

2.2 Interface entre deux diélectriques en incidence normale

On considère deux isolants d’indices réels n1 et n2 occupant respectivement les demi-espaces z > 0 et z < 0.
On veut savoir comment une onde incidente de type OPPM en incidence normale (venant de z = ∞ par exemple)
va être réfléchie et transmise par l’interface z = 0. Le vecteur d’onde incident est donc ~ki = −k1 ~uz (avec k1 > 0),
le vecteur d’onde de l’onde réfléchie est ~kr = k1 ~uz , et celui de l’onde transmise dans le milieu 2 s’écrit ~kt = −k2 ~uz
(avec k2 > 0). On considère une polarisation rectiligne, par exemple selon ~ux . Donc :
~ i (~r, t) = Ei0 ei(−k1 z−ωt) ~ux ,
E (13.22)
~ r (~r, t) = rEi0 ei(k1 .z−ωt) ~ux ,
E (13.23)
~ t (~r, t) = tEi0 ei(−k2 z−ωt) ~ux ,
E (13.24)
et les champs magnétiques sont :

~ i (~r, t) = −n1 Ei0 ei(−k1 z−ωt) ~uy ,

B (13.25)
c
~ r (~r, t) = n1 rE i0
B ei(k1 .z−ωt) ~uy , (13.26)
c
~ t (~r, t) = −n2 tEi0 ei(−k2 z−ωt) ~uy ,
B (13.27)
c

78
où on a utilisé les relations de dispersion dans les milieux 1 et 2.
A incidence normale, tous les champs sont purement tangentiels à l’interface ce qui simplifie les calculs.
De plus dans un isolant les courants surfaciques sont nuls et on a donc continuité des champs électriques et
magnétiques tangentiels (voir problème 4 pour les cas plus compliqués d’incidence oblique). Cela donne :
~ i (z = 0, t) + E
E ~ r (z = 0, t) = E
~ t (z = 0, t), (13.28)
~ i (z = 0, t) + B
B ~ r (z = 0, t) = B
~ t (z = 0, t), (13.29)

ce qui donne deux équations permettant de déterminer les deux inconnues r et t :

1 + r = t, (13.30)
n1 (1 − r) = n2 t, (13.31)

d’où les coefficients en amplitude :

n1 − n2 2n1
r= , avec t= . (13.32)
n1 + n2 n1 + n2

2.3 Interface vide-métal parfait

Ici l’onde transmise sera évanescente. A rédiger.

79
Chapitre 14

Electrodynamique des supraconducteurs

Le but de ce chapitre est de montrer qu’un supraconducteur diffère d’un conducteur normal de conductivité
infinie (le conducteur parfait). En effet, un supraconducteur possède en plus des propriétés magnétiques remar-
quables dont le conducteur parfait est dépourvu. Typiquement un champ magnétique pénètre facilement dans un
bon conducteur, tandis qu’il a tendance à être complètement expulsé (effet Meissner) dans un supraconducteur.
Tout cela est codé dans la relation constitutive des frères London.

1 Introduction historique
Certains métaux sont supraconducteurs. Cela signifie qu’ils perdent toute résistivité électrique (et toute
dissipation par effet Joule) en dessous d’une certaine température critique (notée Tc ). C’est le cas du niobium
(Tc = K), de l’aluminium, du plomb et de l’étain. En revanche, d’autres très bons conducteurs comme le
cuivre ou l’or ne présentent pas de phase supraconductrices jusqu’aux températures les plus basses accessibles.
Historiquement, le mercure (Hg) fut le premier supraconducteur à être découvert en 1911 par Heike Kamerlingh
Onnes. En 1986, famille des supraconducteurs à haute Tc .
En plus de l’annulation de résistivité, les supraconducteurs ont tendance à expulser le champ magnétique
appliqué : effet Meisner-Ochsenfeld (1933).

2 Relation constitutive des frères London

Dans le but d’expliquer la supraconductivité, les frères London postulent un modèle à deux fluides.

2.1 Modèle à deux fluides et fonction d’onde macroscopique

Une partie du gaz électronique se comporte comme dans un métal normal, et l’autre partie acquiert un
caractère quantique macroscopique.

2.2 Relation phénoménologique

On ne va pas démontrer la relation constitutive de London. On va juste motiver cette expression.
Partons du modèle de Drude AC qui s’écrit :

ne2 τ 1
σ(ω) = , (14.1)
m 1 − iωτ
où n est la densité électronique, m la masse de l’électron, τ le temps de relaxation et ω la fréquence.
Dans la limite τ → ∞, la conductivité AC devient purement imaginaire et la conduction est purement non
dissipative car le courant est en quadrature de phase avec le champ électrique.
Cependant la règle de somme sur la partie réelle de la conductivité doit être respectée, ce qui signifie qu’il
existe un pic de Drude à ω = 0 qui porte tout le poids de la partie réelle.
Ainsi, on peut supposer que le fluide supraconducteur à pour conductivité électrique :

πns e2 ns e2 1
σ(ω) = δ(ω) − , (14.2)
me me iω

80
 Ensuite on peut écrire l’équation de Maxwell-Faraday qui donne deux expressions pour le rotationnel du
champ électrique :
!
~j
   
−→ ~ ~ −
→ ~ −→ ~ −
→ −→ ~ me −→ ~ me
rotE = iω B = iω rotA et rotE = rot = rot j iω = iω rot j , (14.3)
σ ns e2 ns e 2

D’où par identification le relation constitutive plausible :

2
~j = − ns e A
~, (14.4)
me
qui est l’équation phénoménologique de London. Cette relation est surprenante car elle relie une grandeur
physique (un courant) supposée être indépendante du choix de jauge, à un potentiel vecteur magnétique A ~ qui
dépend de la jauge. Cette relation ne peut avoir de sens que dans une jauge donnée. Comment savoir laquelle ?
Il faut que l’équation de conservation de la charge soit vérifiée. Cela conduit à la condition de jauge dite de
London div A ~ = 0 (parfois aussi appelé dans le contexte de l’électrostatique la jauge de Coulomb).
On peut remarquer qu’il s’agit d’une relation de réponse linéaire entre courant et potentiel vecteur. Au
niveau des unités, le coefficient entre j et A doit être homogène à 1/L2 /µ0 . Ainsi le préfacteur de la relation
de London ci-dessus introduit une longueur caractéristique importante. il s’agit de la longueur de London λL
définie par la relation :
1 ns e 2
= , (14.5)
λ2L µ0 me
d’où : r
me
λL = . (14.6)
µ 0 ns e 2
La longueur de London diminue (écrantage plus efficace) avec la densité supraconductrice.

3 Effet Meissner-Ochsenfeld
Ici on va écrire l’électrodynamique d’un matériau décrit par la relation constitutive de London + les équations
de Maxwell. On obtient une équation découplée pour A ~ seul, ou de manière équivalente pour B ~ seul. Cette
équation montre qu’un champ magéntique doit décroitre dans un tel matériau de manière exponentielle sur une
longueur λL . La longueur de London est donc la longueur caractérisant l’écrantage du champ magnétique dans
un supraconducteur.

81
Chapitre 15

Annexe sur l’analyse vectorielle

Voici un bref rappel des définitions des opérateurs différentiels utiles pour le reste de ce cours. Ce rappel
est donné en coordonnées cartésiennes (il existe des formulaires pour ces opérateurs exprimés en coordonnées
sphériques ou cylindriques).

1 Gradient
~ =
Le gradient d’un champ scalaire f (~r) s’obtient en appliquant les trois composantes de l’opérateur ∇
(∂/∂x, ∂/∂y, ∂/∂z) sur cette fonction scalaire :

→
− ∂f (~r) ∂f (~r) ∂f (~r)
∇f (~r) = ~ux + ~uy + ~uz (15.1)
∂x ∂y ∂z

Exercice : Montrer que l’accroissement infinitésimal df de la fonction f lorsque ~r → ~r + d~r

2 Divergence
~ r) = (Ax (~r), Ay (~r), Az (~r)) s’obtient en faisant le produit scalaire de
La divergence d’un champ de vecteur A(~
→
−
l’opérateur différentiel nabla ∇ = (∂/∂x, ∂/∂y, ∂/∂z) avec le champ vectoriel :

→
− ~ ∂Ax (~r) ∂Ay (~r) ∂Az (~r)
∇.A(~r) = + + (15.2)
∂x ∂y ∂z

Exercice : Montrer que la divergence du gradient d’un champ scalaire f (~r) s’exprime comme :
→
− →
−
∇. ∇(f (~r)) = ∆f (~r), (15.3)

où le symbole ∆ dénote un opérateur différentiel (appelé le Laplacien) et défini par :

 2
∂2 ∂2

∂
∆f (~r) = + + f (~r). (15.4)
∂x2 ∂y 2 ∂z 2

3 Rotationnel
Le rotation d’un champ de vecteur A(~~ r) = (Ax (~r), Ay (~r), Az (~r)) s’obtient en faisant le produit vectoriel de
→
−
l’opérateur différentiel nabla ∇ = (∂/∂x, ∂/∂y, ∂/∂z) avec le champ vectoriel :
     
→
− ~ ∂Az (~r) ∂Ay (~r) ∂Ax (~r) ∂Az (~r) ∂Ay (~r) ∂Ax (~r)
∇ ∧ A(~r) = − ~ux + − ~uy + − ~uz (15.5)
∂y ∂z ∂z ∂x ∂x ∂y

Exercice 1 : Montrer que la divergence du rotationnel d’un champ vectoriel est identiquement nulle.
Exercice 2 : Montrer que le rotationnel d’un gradient est identiquement nul.

82