Cours de MécaStat

COURS DE MÉCANIQUE STATISTIQUE
MASTER PREMIÈRE ANNÉE

UNIVERSITÉ LYON-1
cel-00339827, version 1 - 19 Nov 2008
Monique Combescure
3 janvier 1970
TABLE DES
MATIÈRES
I] QUELQUES NOTIONS ÉLÉMENTAIRES EN STATISTIQUE
I-1 Variables aléatoires
I-2 Distributions de Probabilités à plusieurs Variables
I-3 Sommes de Variables Aléatoires
II] RAPPELS DE MÉCANIQUE CLASSIQUE

II-1 N Particules Classiques de masse m
II-2 L’Espace de Phase
II-3] L’équation de Liouville
II-3-1 Le Rôle de l’Ensemble en Mécanique statistique
II-3-2 Le Flux dans l’Espace de Phase
cel-00339827, version 1 - 19 Nov 2008
II-3-3 Les Ensembles Stationnaires
III] LES ENSEMBLES D’ÉQUILIBRE

III-1 Introduction
III-2 L’Ensemble Microcanonique
III-3 L’Ensemble Canonique
III-4 L’Ensemble Grand-Canonique
IV] L’ÉQUIVALENCE DES ENSEMBLES D’ÉQUILIBRE

IV-1 La Limite Thermodynamique
IV-2 L’Extensivité des Propriétés Thermodynamiques
IV-3 La Dispersion (les “Fluctuations”) de l’énergie dans l’Ensemble canonique
IV-4 Relation entre les Ensembles Microcanonique et Canonique
IV-5 La Méthode de la Phase Stationnaire (parenthèse mathématique)
IV-6 Relation entre les Ensembles Canonique et Grand-Canonique
V] LES PRINCIPES DE LA THERMODYNAMIQUE

V-1 L’Équilibre entre deux Systèmes en Contact
V-2 Le Deuxième Principe de la Thermodynamique
VI] APPLICATION AUX OSCILLATEURS HARMONIQUES
1
2
VI-1 La Fonction de Partition Canonique de l’Oscillateur Harmonique
VII] LES SYSTÈMES MAGNÉTIQUES

VII-1 Systèmes Magnétiques
VII-2 Mécanique Statistique du Modèle d’Ising
VII-3 Calcul de la Fonction de Partition Canonique d’un Système Magnétique Simple :
spins sans interaction
VIII] LES TRANSITIONS DE PHASE

VIII-1 Observations Expérimentales
VIII-2 L’Approximation de Champ Moyen pour un Ferromagnétique
VIII-3 L’Approximation de Champ Moyen pour un gaz
IX] MÉCANIQUE STATISTIQUE DES SYSTÈMES QUANTIQUES

IX-1 Rappels en Mécanique Quantique
IX-2 La Matrice Densité
cel-00339827, version 1 - 19 Nov 2008
IX-3 Les Matrices-Densité Stationnaires
X] LES GAZ PARFAITS QUANTIQUES

X-1 L’Hamiltonien du Gaz Parfait et ses Fonctions Propres
X-2 Caractérisation à l’aide des Nombres d’Occupation
X-3 Calcul de la fonction de Partition d’un Système Quantique de Particules Indépendantes
XI] GAZ DE FERMIONS INDÉPENDANTS

XI-1 Caractérisation de l’État Fondamental
XI-2 Thermodynamique à Très Basse Température
XII] GAZ DE BOSONS INDÉPENDANTS
XIII] GAZ DE PHOTONS

XIII-1 Les Équations de Maxwell
XIII-2 Quantification des modes de Vibration du Champ Électromagnétique
XIII-3 Le Corps Noir : un Système de Photons à l’Équilibre
XIV] SYSTÈMES MOLÉCULAIRES

XIV-1 Gaz monoatomiques
XIV-2 Gaz de molécules diatomiques
XIV-3 Conclusions
Chapitre 1
Quelques Notions Élémentaires en

Statistique
1.1 Variables Aléatoires

cel-00339827, version 1 - 19 Nov 2008
Exemple 1 :
On jette un dé. Le résultat sera 1 ou 2 ou... ou 6, mais ne peut être prédit : c’est une
variable aléatoire (Dé se dit ALEA en latin).
Soit P (n) la probabilité que le résultat soit n (n = 1, 2, ..., 6). Pour un dé “non-pipé”,
on a évidemment
1
P (n) = , ∀n = 1, 2, ...6
6
tandis que pour un “dé pipé” les P (n) peuvent être différents de 1/6. Mais dans tous les
cas on aura
0 ≤ P (n) ≤ 1
6
X
P (n) ≡ P (1) + P (2) + ... + P (6) = 1
1
Cette condition dite “de normalisation” exprime qu’on est sûr que n prendra l’une des
valeurs 1, 2, ..., 6.
Plus généralement une variable aléatoire (discrète) n est spécifiée par :

(i) L’ensemble X des valeurs qu’elle peut prendre
(ii) Une fonction P (n) telle que
0 ≤ P (n) ≤ 1, ∀n ∈ X
3
4
(iii) la condition
X
P (n) = 1
n∈X
On appelle P (n) la loi de probabilité de n.
Exemple 2 :
On mesure la direction du vent au sommet de la Tour Eiffel. Le résultat est un angle entre 0
et 2π (0 est le Nord par convention), mais ne peut être prédit : c’est une variable aléatoire.
Soit P (ϕ0 )∆ϕ0 le probabilité que la direction du vent se situe dans l’intervalle
[ϕ0 − ∆ϕ0 , ϕ0 + ∆ϕ0 ]. Il est sous-entendu que ∆ϕ0 est suffisamment petit pour que
la variation de P (ϕ) dans cet intervalle soit négligeable.
La probabilité P (ϕ0 )∆ϕ0 est égale à l’aire hâchurée sur la figure.
>
cel-00339827, version 1 - 19 Nov 2008
>φ
Elle a une valeur comprise entre 0 et 1. Afin d’obtenir la condition de normalisation, il

faut diviser l’intervalle [0, 2π] en petits sous-intervalles de longueur ∆ϕ0 et sommer sur
les aires hâchurées correspondantes. À la limite ∆ϕ0 → 0 ceci donne
Z 2π
dϕP (ϕ) = 1
0
comme condition de normalisation (c’est l’aire sous la courbe).

Plus généralement, une variable aléatoire continue x est spécifiée par
(i) L’ensemble X des valeurs qu’elle peut prendre
(ii) Une fonction P (x) non-négative telle que
Z
dxP (x) = 1 normalisation
X
Remarquer qu’il n’y a pas de borne supérieure pour P (x). La fonction x 7→ P (x) est
appelée la loi de probabilité ou la densité de probabilité.
5
Souvent on ne fait pas la distinction entre les cas continu et discret et on dit que :
“P (x) est la probabilité d’obtenir la valeur x”
quand on devrait dire :
“P (x)∆x est la probabilité d’ontenir une valeur dans l’intervalle [x − ∆x, x + ∆x]”.
Ce léger abus de langage est permis quand on a compris ce qu’est une distribution de
probabilté.
Les variables aléatoires en physique :

Les variables aléatoires telles qu’on les a définies ci-dessus pourraient être considérées
comme des objets purement mathématiques. On serait alors libre de choisir la fonction
P (x) sans aucune contrainte. Cependant comme le suggèrent les exemples donnés plus
haut, on considérera par la suite certaines grandeurs physiques comme des variables aléatoires,
par exemple :
- Le nombre d’atomes dans un sous-volume délimité de 1 cm3 à l’intérieur d’un récipient
contenant un gaz
cel-00339827, version 1 - 19 Nov 2008
- Le nombre de désintégrations par seconde dans un échantillon donné d’un matériau ra-
dioactif
- La vitesse d’une particule spécifiée dans un gaz contenu dans un récipient
De façon générale il s’agit de grandeurs bien déterminées, mais dont la valeur x varie
rapidement avec le temps, sans qu’on puisse (ou veuille) mesurer cette variation exacte-
ment. On introduit alors une loi de probabilité P (x), censée résumer nos connaissances sur
la grandeur en question à chaque instant.
Il est évident que les lois P (x) décrivant des grandeurs physiques ne peuvent être arbi-
traires. Elles doivent être compatibles avec les lois fondamentales de la physique. Les lois
de probabilité employées dans ce cours seront :
- ou bien déduites de quelques postulats de base plus les lois de la physique
- ou bien postulées sans aucune démonstration, parce qu’elles sont “plausibles”.
Dans tous les cas on les considérera comme justifiées si les résultats de la théorie sont
conformes aux résultats expérimentaux.
Finalement, si la grandeur physique en question est telle qu’on peut l’observer directement
dans une expérience, on peut en déterminer la distribution de probabilité expérimentalement.
P (x) sera connue avec une précision d’autant plus grande que le nombre d’observations
est grand.
EXEMPLES DE LOIS DE PROBABILITÉS :

6
• distributions continues :
1) distribution exponentielle :
P (x) = ae−ax, x ∈ X := [0, ∞), a>0
2) Loi de Lorentz ou de Cauchy
1 a
P (x) = x ∈ R, a 6= 0
π x2 + a2
3) distribution uniforme sur un intervalle :
1
P (x) = , a < b, x ∈ [a, b]
b−a
4) Loi Gaussienne (ou loi normale) :
cel-00339827, version 1 - 19 Nov 2008
(Gauss 1809, Laplace 1790, De Moivre 1733)

1 (x − a)2
P (x) = √ exp − , b 6= 0, x∈R
2πb2 2b2
R
Exercice 1 : Vérifier que dans tous les cas on a bien X
dxP (x) = 1
• Distributions discrètes :
(la variable aléatoire est appelée n au lieu de x)
5) loi de Poisson :
−a
an
P (n) = e , n∈N
n!
6) loi binomiale :

N
P (n) = pn(1 − p)N−n, 0 < p < 1, n ∈ {1, 2, ..., N }
n
P
Exercice 2 : Vérifier que dans ces derniers cas n P (n) = 1
7
Définition 1.1.1 À partir de la loi de probabilité P (x), on peut définir les quantités
suivantes : R
moyenne de x : x̄ ≡ hxi = dxxP (x)R
moyenne d’une fonction x 7→ gR: ḡ = dxg(x)P (x)
mième moment de x : hxmi = dxxmP (x), m ∈ N
l’écart à la valeur moyenne de x : ∆x = x − x̄
l’écart quadratique moyen, ou variance :
σx2 = h(x − x̄)2 i = h(∆x)2 i
la déviation standard ou écart-type : σx
Exercice 3 : Montrer que σx2 = hx2 i − hxi2
Exercice 4 : Les définitions ci-dessus ont été données pour une distribution de proba-
bilité continue. Écrire les définitions analogues pour une loi de probabilité discrète.
cel-00339827, version 1 - 19 Nov 2008
Exercice 5 : Calculer x̄ et σx pour les exemples 1) à 6) de lois de probabilités donnée

plus haut. Quel nouveau phénomène apparaı̂t dans le cas 2) ? Pour les exemples 5) et 6),
adapter les définitions de x̄ et σx.
Définition 1.1.2 La transformée de Fourier de P (x)

Z +∞
G(k) := dx eikxP (x)
−∞
s’appelle la fonction caractéristique de la variable aléatoire . Si X est stricte-

ment contenu dans R, par exemple un intervalle, il est clair qu’on prendra P (x) =
0, si x ∈ /X
Exercice 6 : Montrer les propriétés suivantes :
|G(k)| ≤ 1, ∀k ∈ R
G(0) = 1
(ik)2 ¯2 (ik)3 ¯3
G(k) = 1 + ikx̄ + x + x + ...
2! 3!
à condition que tous les moments soient finis. On dit que G(k) est la fonction génératrice
des moments.
8
Exercice 7 Soit x une variable aléatoire continue décrite par une loi de probabilité
continue P (x). Soit α un paramètre réel. La quantité y ≡ αx est aussi une variable
aléatoire ; montrer que sa loi de probabilité P̃ (y) est donnée par

1 y
P̃ (y) = P
α α
1.1.1 Distributions de probabilités à plusieurs variables

X peut être un sous-ensemble d’un espace vectoriel. Si x1 et x2 sont deux variables
aléatoires, alors le couple (x1 , x2 ) est aussi une variable aléatoire. De même pour un triplet
(x1 , x2 , x3 ), etc...
Exemple 1 : Une particule dans une boite cubique d’arête L. La variable aléatoire est
la position x := (x1 , x2 , x3 ) de la particule.
x3
<
cel-00339827, version 1 - 19 Nov 2008
< x1
< x2
L
X := [0, L]3 ⊂ R3 , P (x) ≡ P (x1 , x2 , x3 )
et la condition de normalisation s’écrit :

Z L Z L Z L
P (x1 , x2 , x3 )dx1 dx2 dx3 = 1
0 0 0 | {z }
Probabilité que les trois coordonnées cartésiennes
de la particule aient des valeurs entre x1 et x1 + dx1 , x2 et x2 + dx2 , x3 et x3 + dx3 .
Exercice 1 : Montrer que si P (x1 , x2 , x3 ) est uniforme, alors P (x1 , x2 , x3 ) = L−3
Exemple 2 : La variable aléatoire est la vitesse d’une particule : X = R3 , x =

(v1 , v2 , v3 ) ≡ ~
v.
9
Une loi de probabilité particulière pour la vitesse d’une particule, qu’on rencontre fréquemment
en physique est la distribution des vitesses de Maxwell, dénotée PM :
3/2
m v2
m~
v) ≡
PM (~ exp −
2πkB T 2kB T
2
v ≡
où ~ v12 + v22 + v32
m est la masse de la particule
T est la température
kB est la constante de Boltzmann
Cette loi sera démontrée plus tard dans ce Cours.
Exercice 3 : Montrer que la distribution de Maxwell est normalisée.
Exemple 4 : La variable aléatoire est la paire (position, vitesse) d’une particule. Donc
x ≡ (~ v ) = (x1 , x2 , x3 , v1 , v2 , v3 ) ∈ R6
r, ~
cel-00339827, version 1 - 19 Nov 2008
Exercice 3 : Montrer que la distribution particulière dans X ≡ [0, L]3 × R3
P (~ v ) = L−3 PM (~
r, ~ v)
est normalisée.
Définitions :
On considérera le cas d’une distribution P (x, y) à deux variables réelles. Les généralisations
à plus que deux variables sont évidentes. Soit g(x, y) une fonction arbitraire de x, y.
Moyenne de g : Z Z
g(x, y) = dx dy g(x, y)P (x, y)
Moments : Z Z
xmy n = dx dy xmy nP (x, y)
Lois marginales :
Il arrive souvent que l’on considère des fonctions g(x, y) qui ne dépendent que de l’une
des variables, soit de x par exemple : g(x, y) = h(x). On a alors :
Z Z
g(x, y) = h(x) = dxh(x) dyP (x, y)
| {z }
≡ P1 (x)
10
On peut donc considérer h(x) comme la moyenne de h(x) par rapport à la loi de
probabilité P1 (x) qui ne dépend que d’une seule variable. Interprétation : P1 (x) est la
probabilité d’obtenir la valeur x, quelle que soit la valeur y que l’on obtient. On appelle
P1 (x) la loi marginale pour x, déduite de la loi P (x, y).
Exercice 5 : Définir également la loi marginale pour y. Vérifier que les distributions
marginales sont positives et normalisées.
Exercice 6 : À partir de la distribution des vitesses de Maxwell, trouver la loi marginale

pour la composante v1 de la vitesse d’une particule.
Exercice 7 : Prendre g(x, y) = x + y ; Montrer que
x + y = x̄ + ȳ
Définition 1.1.3 Soit P (x, y) ≥ 0 une loi de probabilité ; alors

cel-00339827, version 1 - 19 Nov 2008
P (x, y)
P (x|y) := R
dx P (x, y)
est la loi de probabilité de x quand la valeur de y est connue.
R
P2 (y) = dx P (x, y) est la loi marginale pour y. On a
P (x, y) = P2 (y)P (x|y)
Ce concept d’une loi de probabilité (conditionnelle) pour une variable sachant la valeur de
l’autre se généralise de façon évidente aux distributions à plusieurs variables.
11
Définition 1.1.4 Soit P (x, y) une loi de probabilité à deux variables. On définit les
variables aléatoires (de moyenne nulle) :
∆x = x − x̄
∆y = y − ȳ
On introduit la matrice variance µij , i, j ∈ {1, 2} :
µ11 = ∆x2 = σx2
µ22 = ∆y 2 = σy2
µ12 = µ21 = ∆x∆y
µ11 , µ22 sont les variances, et µ12 = µ21 les covariances.
Si ∆x∆y = 0, on dit que x, y sont décorrélées.
Variables aléatoires indépendantes

cel-00339827, version 1 - 19 Nov 2008
Définition 1.1.5 (i) Deux variables aléatoires sont dites indépendantes si leur loi
de probabilité P (x, y) s’écrit comme
P (x, y) = Q(x)R(y)
où Q, R sont des lois de probabilité arbitraires.

Si Q = R, x et y sont appelées variables aléatoires indépendantes et identiques.
(ii) N variables aléatoires x1 , x2 , ..., xN sont indépendantes si
P (x1 , x2 , ..., xN ) = Q1 (x1 )Q2 (x2 )...QN (xN )
Exercice 8 : Montrer que pour la distribution des vitesses de Maxwell les trois composantes
de la vitesse d’une particule sont des variables aléatoires indépendantes (et identiques).
Exercice 9 : Montrer que si deux variables aléatoires sont indépendantes, elles sont non-
corrélées, mais que la réciproque n’est pas vraie. Contre-exemple :
1 1
P (x, y) = δ(x − 1)[δ(y − 1) + δ(y + 1)] + δ(x + 1)δ(y)
4 2
1.1.2 Sommes de variables aléatoires

Soient {x1 , x2 , ..., xN } N variables aléatoires ayant une loi de probabilité connue
P (x1 , x2 , ..., xN ). Ceci implique en particulier que les valeurs xi et ∆x2i sont connues
12
pour tout i = 1, 2, ...N .

La somme X = x1 + x2 + ... + xN est aussi une variable aléatoire. Soit Π(X) sa loi de
probabilité.
On se posera successivement les questions suivantes :

(1) Peut-on calculer la moyenne X̄ ?
(2) Peut-on calculer la variance ∆X 2 ?
(3) Peut-on calculer la loi Π(X) ?
(1) La réponse est la suivante :
X̄ = x1 + x2 + ... + xN
pour une distribution de probabilité P quelconque.

(2) On a :
N
X
cel-00339827, version 1 - 19 Nov 2008
∆X = ∆xi
i=1
N
X X
∆X 2 = ∆x2i + 2 ∆xi∆xj
i=1 i<j
Si les variables xi, i = 1, ...N sont non-corrélées, le second terme s’annule et on a
Proposition 1.1.6 Si les variables xi i = 1, ...N sont non-corrélées, on a

N
X
∆X 2 = ∆x2i
i=1
c’est à dire que la variance d’une somme de variables aléatoires non-corrélées est
la somme de leurs variances.
Un cas particulier est celui où les xi sont des variables aléatoires indépendantes et iden-
tiques, donc P (x1 , x2 , ...xN ) = Q(x1)Q(x2 )...Q(xN ). On a alors :
X = N x1 ∆X 2 = N ∆x21
Ces relations sont vraies quelle que soit la loi de probabilité Q(x).
Remarque Au lieu de considérer la somme X, on peut aussi s’intéresser à la “moyenne
arithmétique”
1
x= (x1 + x2 + ... + xN )
N
13
Pour celle-ci, on a (pour des variables aléatoires indépendantes et identiques) :
x̄ = x1
1
∆x2 = ∆x21
N
ou bien
1
σx = √ σx1
N
La moyenne arithmétique de N variables aléatoires indépendantes et identiques a une
déviation standard qui est réduite d’un facteur N −1/2 par rapport à la déviation-
standard d’une seule variable.
(3) Pour N grand une réponse à la question (3) est donnée par le Théorème de la Limite
Centrale :
Théorème 1.1.7 Supposons que les xi i = 1, ...N soient des variables aléatoires
indépendantes et identiques avec la loi de probabilité Q(x) quelconque. Supposons que
cel-00339827, version 1 - 19 Nov 2008
l’on ait
∆x2i < ∞
Alors, X et ∆X 2 étant données ci-dessus, on a pour N → ∞
!
1 (X − X)2
Π(X) ≃ p exp −
2π∆X 2 2∆X 2
Remarque
(i) Les réponses (1) et (2) sont vraies quel que soit N. Par contre le Théorème de la
Limite Centrale n’est vrai qu’asymptotiquement quand N devient grand.
(ii) Π(X) est une Gaussienne, lorsque N devient très grand, quelle que soit la distri-
bution de probabilité Q(x) qu’on avait au départ (Pourvu néanmoins que les variances
soient finies). Pour cette raison on dit que la Gaussienne est la “loi limite universelle”.
La loi de probabilité asymptotique Π̃(x) pour la moyenne arithmétique des xi se déduit
facilement du Théorème de la Limite Centrale :
√ !
N N (x − x̄)2
Π̃(x) ≃ q exp − N →∞
2π∆x21 2∆x21
Chapitre 2
RAPPELS DE MÉCANIQUE
CLASSIQUE
2.1 N particules classiques de masse m

cel-00339827, version 1 - 19 Nov 2008
Soient N particules identiques de masse m et de positions xi(t), i = 1, 2, ...N

dans un champ extérieur U (x1 , x2 , ...xN ) (par exemple le champ de la pesanteur), et
interagissant deux à deux par un potentiel V (xi − xj ). Un exemple de tel potentiel est
le “potentiel empirique de Lenard-Jones” qui décrit une répulsion à courte portée, suivie
d’une attraction à plus longue portée. Pour un gaz noble constitué par exemple d’atomes
d’Argon, le potentiel est décrit par une courbe qui est la suivante :
V(r)
>
r0
> r
V0
où r0 ≃ 3.8Angström et V0 ≃ −0.01 eV.

Soient pi i = 1, ..., N les impulsions des particules. On a pi = mvi où vi = dx dt
i
est le
vitesse instantanée de la particule i. Écrivant que l’énergie totale est la somme de l’énergie
14
15
cinétique et de l’énergie potentielle pour le système à N -corps ainsi formé, on a :

N
X p2i X
H(x1 , ..., xN , p1 , ..., pN ) = + V (xi − xj ) + U (x1 , ..., xN )
i=1
2m i<j
H considéré comme une fonction des positions et des impulsions de toutes les particules
est appelée “fonction Hamiltonienne”. Elle permet en effet de “résoudre” les équations
du mouvement via les équations dites d’Hamilton qui sont :
dxi ∂H dpi ∂H
ẋi ≡ = ≡ vi ṗi ≡ =− (2.1)
dt ∂pi dt ∂xi
2.2 L’espace de Phase

La variable collective décrivant l’ensemble des positions et impulsions des N particules
est appelée “micro-état” et sera notée Γ. Elle “vit” dans un espace à 6N -dimensions appelé
cel-00339827, version 1 - 19 Nov 2008
espace de phase. La dimension de l’espace de phase peut donc être gigantesque pour un
ensemble macroscopique composé de N = 1023 particules.
Les équations (2.1) permettent d’associer à toute donnée Γ0 à l’instant 0 un point de
l’espace de phase à l’instant t noté Γ(t). L’ensemble des points Γ(t) quand t varie est
appelé “trajectoire dans l’espace de phase”.
Théorème 2.2.1 On a pour tout temps t conservation de l’énergie :
H(Γ(t)) = H(Γ0 ) (2.2)
Preuve : (2.2) résulte de (2.1). En effet

N
dH(Γ(t)) X ∂H ∂H
= · ẋi + · ṗi = 0
dt i=1
∂xi ∂pi
Exercice 1 : Montrer que deux trajectoires différentes dans l’espace de phase ne peuvent
pas se couper.
Exercice 2 : Peut-il y avoir des trajectoires en boucle ?
En vertu du Théorème ci-dessus, les trajectoires Γ(t) dans l’espace de phase “vivent” sur
la surface d’énergie H(Γ(t)) = Cste de dimension 6N − 1.
On doit se demander s’il existe d’autres invariants que l’énergie. Si c’était le cas, alors la
16
surface d’énergie se diviserait en “secteurs”, c’est à dire des hypersurfaces de dimension

plus basse que 6N − 1 sur lesquels les trajectoires seraient confinées.
En général (sauf cas particulier), on admet que ce n’est pas le cas et que H est le seul
invariant du système. C’est l’hypothèse ergodique.
2.3 L’Équation de Liouville

2.3.1 Le rôle de l’Ensemble en Mécanique Statistique
En pratique, il est impossible de résoudre (2.1) dès que N devient grand, donc il est
impossible de connaı̂tre le “micro-état” Γ(t). Pour combler ce manque de connaissance,
on introduit une loi de probabilité P (Γ) telle que P (Γ)dΓ soit la probabilité de trouver
ce système dans un élément de volume dΓ ≡ dx1 ...dxN dp1 ...dpN autour de point Γ de
l’espace de phase. On doit bien évidemment avoir
Z
cel-00339827, version 1 - 19 Nov 2008
dΓ P (Γ) = 1 (2.3)
R6N
Le choix de P (Γ) est motivé par les conditions expérimentales du système. En général,
une distribution P (Γ) évoluera dans le temps. Pour un système à l’équilibre, P (Γ) est
stationnaire (indépendante du temps).
Observables :
Soit A une grandeur physique mesurable (qu’on appellera “observable”), qui a pour valeur
A(Γ) quand le système se trouve dans le micro-état Γ. On peut alors calculer sa moyenne :
Z
hAi := dΓ P (Γ)A(Γ) (2.4)
On considère cette moyenne comme la prédiction théorique pour le résultat d’une me-
sure expérimentale. Pour justifier cette identification, on note d’abord que l’intervalle de
temps ∆t nécessaire pour une mesure expérimentale macroscopique est très long à l’échelle
microscopique. C’est dire que le micro-état Γ(t) varie beaucoup à l’échelle de temps ∆t.
Soit donc R la région de l’espace de phase où P (Γ) diffère sensiblement de zéro.
17
P(Γ) ≈ 0
P(Γ) > 0
Γ(t)
Γ(t + ∆t)
Le vrai résultat expérimental sera la valeur de A(Γ) moyennée le long de la trajectoire

entre les instants t et t + ∆t :
Z t+∆t
1
Aexp = dsA(Γ(s)) (2.5)
∆t t
Mais comme ∆t est long à l’échelle microscopique, ce segment de trajectoire sera nécessairement
long et tortueux, et l’hypothèse essentielle est qu’il passe par une suite de points Γ(s) suf-
cel-00339827, version 1 - 19 Nov 2008
fisamment représentative pour la distribution P (Γ) pour qu’on puisse identifier (2.5) et
(2.4).
On considérera désormais l’identification de (2.4) avec la valeur expérimentale comme

acquise.
On considère souvent une densité ρ(Γ) dans l’espace de phase, avec ρ(Γ) ≥ 0, mais ρ
non normalisée. On peut toujours passer d’une densité ρ à une probabilité P (Γ) par la
relation
ρ(Γ)
P (Γ) = R (2.6)
dΓ ρ(Γ)
On peut représenter une fonction ρ(Γ) ≥ 0 sur l’espace de phase par une collection de
points qui sont localement d’autant plus denses que ρ(Γ) est grande.
.
. .
. · .· .
.
. .
.· .
.
. ·· ·
· . ·
. · · · ·
.
. · .
. . .. · ·
· . ·. . · .
· . .
·
··••·····••···•······•
.
· .
·. . . ·
.
•
· •····••·•·••·•·•·••··•·•
.
. .
· .
. ·· .
··
·· ·· .
. · · ·
··•··•··• ·•· •···•
· .··· .
· .
· ·• ••·• · •
· .
. · . .
·
··
. . .. .. . ·· . .·
.
.
. ·. . ·
. .
. ·
. · ·
··
.
·
· ·
. · ·
·
·
. .
. · . . ·.
. .
espace de phase
18
Ensuite, abandonnant l’interprétation originale d’une loi de probabilité pour un seul

système physique, on peut imaginer que chacun de ces points représente un système phy-
sique différent. D’où le nom “ensemble” pour ρ(Γ).
Remarque 2.3.1 Le formalisme des ensembles que l’on va développer n’exclut pas
d’étudier, le cas échéant,un micro-état précis. Si on sait que le système est dans le
microétat Γ0 , on peut exprimer cette connaissance en posant
ρ(Γ) = δ(Γ − Γ0 )
où δ est la distribution de Dirac.
2.3.2 Le Flux dans l’Espace de Phase

Étant donné un ensemble de points de l’espace de phase à l’instant t0 , on le décrit par
une densité ρ(Γ, t0 ) où l’on a rendu explicite le temps. Maintenant on “met le temps en
cel-00339827, version 1 - 19 Nov 2008
marche”, c’est à dire que chaque point dans le support de ρ(Γ, t0 ) est considéré comme le
point de départ d’une trajectoire dans l’espace de phase :
t 7→ Γ(t) avec Γ(t0 ) = Γ
À tout temps t > t0 les bouts des trajectoires constituent un nouvel ensemble de points,
décrit par une nouvelle densité ρ(Γ, t). Ce sont les équations du mouvement (2.1) qui
déterminent les trajectoires, donc la relation entre ρ(Γ, t0 ) et ρ(Γ, t).
>
pN
p(Γ,t)
>
>
>
p(Γ,t0)
> rN
Cependant avant même d’introduire ces équations du mouvement, on notera que puis-
qu’aucun point à l’instant t0 ne peut disparaı̂tre au cours du temps, de même qu’aucun
nouveau point sans antécédent ne peut apparaı̂tre, la densité ρ est une quantité conservée.
Elle satisfait à l’équation de conservation :
19
∂ρ dΓ
+ ∇ · (ρV ) = 0 où V (Γ) =
∂t dt
Cette équation peut être explicitée de la manière suivante :
∂ρ
+ V · ∇ρ + ρ∇ · V = 0
∂t
où ∇ρ est le vecteur 6N -dimensionnel
∂ρ ∂ρ
∇ρ ≡ ( , ..., ) gradient de ρ
∂x1 ∂pN
et
N X
3
X ∂ ẋik ∂ ṗik
∇·V ≡ + divergence de V (2.7)
i=1 k=1
∂xik ∂pik
Employant maintenant l’équation du mouvement (2.1), on a :

∂H ∂H ∂H ∂H
cel-00339827, version 1 - 19 Nov 2008
V = , ..., ,− , ..., −
∂p1 ∂pN ∂x1 ∂xN
et par suite utilisant la définition (2.7) de ∇ · V ,on a
∇·V =0
c’est à dire que le “flux dans l’espace de phase est incompressible” (un terme employé par
analogie avec l’écoulement d’un fluide réel en trois dimensions).
Il en résulte que
∂ρ dρ
0= + V · ∇ρ = (2.8)
∂t dt
D’où le Théorème de Liouville :
Théorème 2.3.2 La densité ρ(Γ, t) est constante, ou en d’autres termes le volume

dans l’espace de phase est conservé.
L’équation (2.8) peut se réécrire comme :
Équation de Liouville :
N
∂ρ X ∂H ∂ρ ∂H ∂ρ
= − · + ·
∂t i=1
∂pi ∂xi ∂xi ∂pi
20
Cette équation est linéaire en ρ. L’équation de Liouville s’écrit donc aussi comme
∂ρ
= Lρ
∂t
où L, appelé opérateur de Liouville est un opérateur compliqué, mais linéaire, sur
l’espace de toutes les densités ρ(Γ). Cette linéarité permet d’écrire la solution ρ(Γ, t)
comme
ρ(Γ, t) = eLtρ(Γ, 0) (2.9)
Cette solution, cependant, n’est que formelle : l’opérateur eLt, bien que complètement
défini, reste assez insaisissable. Aucune conclusion pratique ne peut être déduite directe-
ment de (9.6).
Exercice 1 : Montrer que l’invariance de l’énergie, exprimée par (2.2) s’écrit aussi comme
V · ∇H = 0
cel-00339827, version 1 - 19 Nov 2008
Exercice 2 : Montrer que :

|V | = |∇H|
Interprétation : le flux de long d’une hypersurface H = cste est d’autant plus rapide
qur les hypersurfaces se trouvent proches les unes des autres (cf. l’effet analogue dans
l’écoulement d’un fluide ordinaire).
2.3.3 Les ensembles stationnaires

En général un ensemble ρ(Γ, t) dépend du temps. Par conséquent toutes les moyennes
hAi calculées via (2.4) et qui correspondent à des observations expérimentales dépendent
du temps.
Or on sait expérimentalement qu’un système physique isolé tend vers un état d’équilibre.
Quand cet équilibre est atteint, les valeurs macroscopiques décrivant le système ne changent
pas avec le temps.
C’est pourquoi on s’intéresse particulièrement à des ensembles ρ(Γ) qui sont station-
naires (indépendants du temps). Il n’est pas difficile d’en trouver : si ρ est une fonction de
H (l’Hamiltonien)
ρ(Γ) = ϕ(H(Γ))
21
on a bien évidemment que

dρ
=0 ∀ϕ
dt
Si le système possède d’autres invariants que H, ϕ peut évidemment aussi dépendre de
ceux-ci. (En pratique ces autres invariants, s’il y en a, ne jouent pas de rôle dans la plupart
des systèmes physiques).
R
Dans la Section 2.3.1, on a déjà considéré la quantité dΓ ρ(Γ) qui a priori intervient
comme facteur de normalisation (pour calculer les moyennes). Pour les ensembles station-
naires, on note cette quantité par le symbole spécifique Z :
Z
Z= dΓ ρ(Γ)
Z est appelée fonction de partition de l’ensemble stationnaire. Elle joue un rôle déterminant
dans la physique statistique de systèmes à l’équilibre.
Exercice : Quelle est la dimension physique de Z ?

cel-00339827, version 1 - 19 Nov 2008
Chapitre 3
LES ENSEMBLES D’ÉQUILIBRE
3.1 Introduction
Comme on l’a vu à la fin du Chapitre précédent, chaque ensemble qui peut s’écrire
comme
cel-00339827, version 1 - 19 Nov 2008
ρ(Γ) = ϕ(H(Γ)) (3.1)

est un ensemble stationnaire. La fonction ϕ est arbitraire, mais en pratique le nombre de
choix est limité. Dans ce Chapitre, on exposera les trois principaux Ensembles Statistiques
à savoir :
• L’Ensemble Microcanonique
• L’Ensemble Canonique
• L’Ensemble Grand-Canonique
et on indiquera les conditions physiques auxquelles chacun d’entre eux correspond.
La physique statistique fournit une base microphysique à la thermodynamique. On verra
en effet qu’on peut retrouver les propriétés thermodynamiques des systèmes à l’aide de
chacun des ensembles ci-dessus. C’est pourquoi on dit que ces ensembles sont équivalents.
Aussi surprenant que cela puisse paraı̂tre, n’importe quel choix de ϕ dans (3.1) donne un
ensemble conduisant à la thermodynamique.
3.2 L’ensemble Microcanonique

L’ensemble microcanonique découle également du postulat de Boltzmann :
Si un système macroscopique est à l’équilibre et possède une énergie entre E et E+δE

alors tous les micro-états de ce système ayant une énergie entre E et E + δE sont a
priori équiprobables
22
23
“A priori” signifie “en l’absence d’informations supplémentaires. Ce postulat nous conduit

à définir l’ensemble microcanonique comme l’ensemble contenant tous les micro-états Γ
ayant une énergie entre E et E + δE avec le même poids :

C si E < H(Γ) < E + δE
ρm(Γ) = (3.2)
0 sinon
Le cas le plus simple est celui où le système est composé de N particules identiques
enfermées dans un volume V . On appelle alors le triplet (N, E, V ) les paramètres fixes,
ou les variables indépendantes de l’ensemble microcanonique.
Remarque 3.2.1 Les deux autres ensembles seront caractérisés par deux autres triplets
de variables indépendantes.
3.2.1 La constante C
cel-00339827, version 1 - 19 Nov 2008
Quelle est sa dimension ?

R
ρm(Γ) est une densité dans l’espace de phase, et dΓ ρm(Γ) est un nombre sans
dimension. Par conséquent :
[dΓ] = (Js)3N donc
[C] = (Js)−3N = (action)−3N
Quelle est son interprétation physique ?

C’est le nombre de micro-états par unité de volume dans l’espace de phase.
La valeur numérique de C peut paraı̂tre sans importance, puisque la probabilité statistique
P (Γ) qui en résulte via (2.6) n’en dépend pas !
Par convention on choisit C tel que ρm(Γ)dΓ est égal au nombre de micro-états quan-
tiques correspondant aux états classiques dans l’élement de volume dΓ considéré ; par
conséquent, on a pour un système classique
Z
dΓ ρm(Γ) = nombre de micro-états quantiques du même
E<H(Γ)<E+δE
système ayant une énergie entre E et E + δE

Pour fixer C, on fait donc appel à la mécanique quantique. On verra que, h étant la
constante de Planck, (h = 6.6 × 10−34 Js) un état quantique à une particule “occupe
un volume de l’espace de phase” R3 égal à h3 . Ici on est dans l’espace de phase à 6N
dimensions. Un micro-état quantique à N particules occupe donc un volume de l’espace
24
de phase égal à N !h3N .

L’apparition du facteur N ! est une conséquence de la mécanique quantique, qui dit que
l’échange des nombres quantiques (des “états”) de deux particules dans un système à N
particules ne donne pas un nouveau micro-état quantique, les particules étant indiscer-
nables.
Par suite on pose :
1
C=
N !h3N
Dans toute la suite de ce Cours, l’indice m fait référence à l’ensemble microcanonique.

Ainsi la fonction de partition microcanonique est définie par
Z Z
1
Zm = dΓ ρm(Γ) = 3N dΓ
h N ! E<H(Γ)<E+δE
3.2.2 Exemples
cel-00339827, version 1 - 19 Nov 2008
1) L’Oscillateur Harmonique Unidimensionnel
p2 1
H(x, p) = + mω 2x2
2m 2
Z
1
Zm = dx dp
h E<H(x,p)<E+δE
L’aire de la région du plan (x, p) contenue entre les deux ellipses H(x, p) = E et
H(x, p) = E + δE est égale à 2πδE/ω d’où
2πδE
Zm =
hω
2) Le Gaz Parfait (enfermé dans un volume V )
N
1 X
H(x1 , ..., xN , p1 , ..., pN ) = p2i
2m i=1
(L’énergie de dépend pas des positions xi des particules)

Z
1 N
Zm = 3N V dp1 ...dpN
h N! E<H(p1 ,..,pN )<E+δE
√
On a donc dans l’intégrale le volume contenu entre deux hypersphères de rayons 2mE
p
et 2m(E + δE). Sachant que le volume d’une hypersphère de rayon R en dimension d
25
est égal à
π d/2
Rd
(d/2)!
on trouve au premier ordre en δE (δE est petit devant E) :

VN 2πmE 3N/2 1 3N δE
Zm(N, V, E) = 3N 2
h N! h (3N/2)! 2 E
N étant grand, on utilise la formule de Stirling :
1
log N ! ≃ N log N − N + log(2πN )
2
ce qui donne :

3 2πmE 3N 3N
log Zm = N log V − N log N + N + N log − log
2 h2 2 2

3N δE
+ + log + O(log N )
2 E
cel-00339827, version 1 - 19 Nov 2008
δE
log N

En regroupant les termes, et en négligeant les termes en O E
et en O N
, il
vient :
V 3 4πmE 5
log Zm = N log + log +
N 2 3h2 N 2
On définit les quantités intensives suivantes :
N E
n= ǫ=
V N
et on obtient à la limite N → ∞ :

5 3 4πmǫ
log Zm =N + log (3.3)
2 2 3(nh3 )2/3
On voit que cette quantité augmente linéairement avec N : elle est dite extensive. Les
quantités V = n−1 N et E = ǫN sont également extensives.
Remarque : Au cours du calcul ci-dessus sont apparus des termes “sur-extensifs” du
type N log N ou N log V ... Ceux-ci ne peuvent figurer dans le résultat final qui représente
une grandeur physique.
Remarque : Dans (3.3) on a trouvé que Zm est exponentielle en N . Ce sera le cas
pour tout système à N corps, et dans chaque ensemble :
Z ≃ eNConstante N →∞
où la constante sera bien sûr fonction des paramètres intensifs.

26
3.2.3 Théorie Générale

À chaque fonction de partition, on associe un potentiel thermodynamique, ainsi
appelé parce qu’il permet de décrire toutes les autres grandeurs thermodynamiques par de
simples dérivations. La fonction de partition microcanonique a un potentiel thermodyna-
mique appelé l’entropie microcanonique :
Sm(N, V, E) = kB log Zm(N, V, E) (3.4)

où kB est la constante de Boltzmann
(kB = 1.38 × 10−23 J K −1 )
Sm dépend des trois paramètres N, V, E. Ses dérivées le long des trois axes principaux
définissent
- La température inverse βm = 1/kB Tm
- la pression pm
- le potentiel chimique µm
cel-00339827, version 1 - 19 Nov 2008
selon les formules suivantes :

1 ∂Sm(N, V, E)
βm(N, V, E) ≡ (3.5)
kB ∂E
1 ∂Sm(N, V, E)
βmpm(N, V, E) ≡ (3.6)
kB ∂V
1 ∂Sm(N, V, E)
βmµm(N, V, E) ≡ − (3.7)
kB ∂N
Remarques :
1) L’indice m rappelle qu’il s’agit de grandeurs définies dans l’ensemble microcanonique.
Il sera supprimé plus tard.
2) Il faudra démontrer que les grandeurs Tm, pm, µm définies ci-dessus coincident avec
celles dont il est question en thermodynamique (et que l’on écrira sans indice).
3) On vérifie aisément que les quantités définies ci-dessus ont la bonne dimension.
4) On trouve les autres grandeurs thermodynamiques à l’aide des formules habituelles, soit
L’énergie libre F = E − T S
L’enthalpie H = E + pV
La Fonction de Gibbs G = N µ
Les formules pour les dérivées partielles permettent d’obtenir la formule suivante pour la
différentielle totale de Sm :
27
1 pm µm
dSm = dE + dV − dN (3.8)
Tm Tm Tm
3.3 L’Ensemble Canonique

L’ensemble microcanonique s’applique, a priori, à des systèmes physiques dont l’énergie
E est strictement constante, c’est à dire des systèmes qui sont parfaitement isolés de leur
environnement.
L’ensemble canonique s’applique, a priori, aux systèmes qui sont en contact thermique
avec leur environnement : ils peuvent absorber ou dégager de l’énergie (de la chaleur). Le
cas idéalisé que l’on considérera ici est celui où l’environnement a une capacité calorifique
infinie, de telle sorte que sa température T est constante. Un tel environnement est appelé
un thermostat.
cel-00339827, version 1 - 19 Nov 2008
Soit un système en contact avec un thermostat. Même à l’équilibre, son énergie peut fluc-
tuer autour d’une moyenne hEi. La valeur de cette moyenne dépend de la température
du système T qui est la même que celle du thermostat. C’est pourquoi dans l’ensemble
canonique, il y a un nouveau paramètre T (sans indice) qui remplace E.
On démontrera plus tard que pour des systèmes macroscopiques (N grand), Tm coincide
avec T .
L’ensemble canonique est défini par le choix suivant de la densité dans l’espace de
phase :
1
ρc(Γ) ≡ 3N e−βH(Γ) (3.9)
h N ! | {z }
Facteur de Boltzmann
On peut voir la distribution canonique comme une somme de distributions microcanoniques
juxtaposées et pondérées avec le poids e−βE .
28
> ƌc(Γ)
e-βH (Γ)
> H (Γ)
E Ε+δΕ Ε+2δΕ Ε+3δΕ
On rappelle que
1
cel-00339827, version 1 - 19 Nov 2008
β=
kB T
Fonction de partition canonique

Z
1
Zc(N, V, T ) = dΓ e−βH(Γ) (3.10)
h3N N ! R6N
Énergie libre canonique
Fc(N, V, T ) = −kB T log Zc(N, V, T ) (3.11)
Entropie canonique
∂Fc(N, V, T )
Sc(N, V, T ) = − (3.12)
∂T
Pression canonique
∂Fc(N, V, T )
pc(N, V, T ) = − (3.13)
∂V
Potentiel chimique canonique
∂Fc(N, V, T )
µc(N, V, T ) = (3.14)
∂N
Différentielle totale
dFc(N, V, T ) = −Sc(N, V, T )dT − pc(N, V, T )dV + µc(N, V, T )dN (3.15)

29
Exercice 1 : Montrer que l’énergie moyenne d’un système décrit par l’ensemble cano-
nique est donnée par
∂ log Zc ∂
hEic = − = (βFc)
∂β ∂β
En partant de ce résultat, et de la formule pour Sc, en déduire que
Fc = hEic − T Sc
3.4 L’Ensemble Grand-Canonique

L’ensemble grand-canonique s’applique, a priori, à des systèmes pouvant échanger de
l’énergie (de la chaleur) et des particules avec leur environnement. Quand le système est
à l’équilibre, l’énergie totale et le nombre total de particules peuvent fluctuer autour de
leurs valeurs moyennes qui sont respectivement hEi et hN i. Les paramètres du système
sont T , µ (le potentiel chimique ) et V (le volume occupé par le système).
cel-00339827, version 1 - 19 Nov 2008
Afin de définir l’ensemble grand-canonique, il convient d’élargir le cadre conceptuel de

l’espace de phase qui nous a suffi jusqu’ici. On considère maintenant en même temps
les espaces de phase pour des systèmes à ... N − 1, N , N + 1,... particules. L’espace de
phase à p particules est noté Γp. Sur cette suite Γ1 , Γ2 , ...ΓN , ... d’espaces de phase, on
considère une densité
1
ρ(p, Γp) = eβµpe−βH(p,Γp )
h3pp!
où H(p, Γp) est l’Hamiltonien à p particules défini sur l’espace de phase Γp.
À cette “densité généralisée” correspondra une définition généralisée des moyennes statis-
tiques :
30
P∞ R
p=0 dΓpA(p, Γp)ρ(p, Γp)
hAig ≡ hA(p, Γp)i = P∞ R (3.16)
p=0 dΓp ρ(p, Γp)
Fonction de partition grand-canonique
∞ Z
X 1
Zg (µ, V, T ) = 3N
eβµN dΓN e−βH(N,ΓN ) (3.17)
N=0
h N! R6N
Potentiel grand-canonique
Ωg (µ, V, T ) = −kB T log Zg (µ, V, T ) (3.18)
Nombre moyen de particules dans l’ensemble grand-canonique
∂Ωg (µ, V, T )
Ng (µ, V, T ) = − (3.19)
∂µ
On définit également la pression et l’entropie grand-canoniques par

cel-00339827, version 1 - 19 Nov 2008
∂Ωg (µ, V, T )
pg (µ, V, T ) = − (3.20)
∂V
∂Ωg (µ, V, T )
Sg (µ, V, T ) = − (3.21)
∂T
Différentielle totale
dΩg(µ, V, T ) = −Sg (µ, V, T )dT − pg (µ, V, T )dV − Ng (µ, V, T )dµ (3.22)
Exercice 2 : À l’aide des définitions précédentes, montrer que
Ng = hN ig
où la moyenne dans l’ensemble grand canonique d’une observable (ici N ) est donnée par
la formule (3.16). Ceci justifie alors que Ng représente bien le nombre moyen de particules
autour duquel fluctue le système dans l’ensemble grand-canonique.
Chapitre 4
L’ÉQUIVALENCE DES
ENSEMBLES D’ÉQUILIBRE
4.1 La Limite Thermodynamique

cel-00339827, version 1 - 19 Nov 2008
En physique statistique, on étudie des systèmes à très grand nombre N de particules. On

pourrait penser que plus N est grand, plus cette étude doit être difficile. Paradoxalement,
le cas limite N → ∞ redevient un cas relativement simple, et cela grâce aux méthodes
statistiques (cela rappelle le calcul des probabilités où la “loi des grands nombres” donne
une limite universelle aux distributions de probabilité : la loi Gaussienne) .
Aussi se sert-on en physique statistique du concept de la limite thermodynamique, qui
n’est rien d’autre que la limite d’un système infini, les grandeurs intensives (voir plus bas)
restant fixes. Plus précisément c’est la limite
 N E
 V , V fixes (microcanonique)
N
V →∞ , T fixes (canonique)
 V
µ, T fixes (grand-canonique)
Bien sûr les systèmes physiques que l’on étudie (hormis l’univers !) ne sont pas infinis ;
néanmoins le nombre de particules d’un système macroscopique (N ≈ 1023 ) est suffisam-
ment grand pour que ses propriétés soient très proches de celles déduites théoriquement,
de la limite thermodynamique. Même un échantillon de gaz d’un volume de 1 mm3 est un
système macroscopique !
31
32
Présentation schématique de la limite thermodynamique (dans un ensemble microca-

nonique) :
8 cm3
1 cm3 6.4 x 1021

32 J
8.1020
0.125 cm3 4J
1020 →∞
0.5 J
4.2 L’extensivité des potentiels thermodynamiques

cel-00339827, version 1 - 19 Nov 2008
Une grandeur physique est appelée extensive si elle est proportionnelle à la taille (au
nombre de particules, ou le volume) du système, tout au moins dans la limite où le système
est très grand.
Plus précisément A est appelé extensif s’il existe une constante a non-nulle telle que
A = aN + [des termes qui croissent moins vite que N ]
ou encore (ce qui est équivalent)
A
= a + [des termes qui tendent vers 0 pour N → ∞]
N
ou encore (équivalent)
A
lim =a
N→∞ N
Un résultat qui est indispensable en thermodynamique est le suivant :
Les potentiels thermodynamiques sont extensifs .
Au chapitre précédent, ce résultat a été démontré pour un cas particulier : l’entropie du

gaz parfait.
Une grandeur physique est dite intensive si dans la limite thermodynamique elle tend vers
une constante.
Il s’ensuit que :
33
- Le rapport entre deux grandeurs extensives est intensif

- Le rapport entre une grandeur extensive et une grandeur intensive est extensif.
Les mêmes règles s’appliquent à la dérivée d’une grandeur par rapport à une autre ; avec
ces règles, on peut véfifier que :
- S, hEi, F... sont extensifs
- p, µ, T... sont intensifs.
Remarque 4.2.1 Si A dépend de plus qu’une grandeur extensive, alors à la limite ther-
modynamique a ne peut dépendre que du rapport entre ces grandeurs. Par exemple pour
l’entropie Sm(N, V, E) on a
Sm(N, V, E) ≃ N s(n, ǫ)
avec n = N/V (densité de particules) ǫ = E/N (énergie par particule) s étant alors
cel-00339827, version 1 - 19 Nov 2008
l’entropie par particule.

Autre exemple :
Fc(N, V, T ) ≃ N f (n, T )
où f est l’énergie libre par particule.

On s’attend donc à ce que la thermodynamique soit complètement déterminée dès qu’on
connaı̂t le potentiel thermodynamique (dans l’un des ensembles) en fonction de deux gran-
deurs intensives.
4.3 La dispersion (les “fluctuations”) de l’énergie dans

l’ensemble canonique
Dans l’Exercice 1 du Chapitre précédent, on a vu que
∂
hHic = (βFc)
∂β
On considère maintenant la dérivée seconde :
∂2 ∂hHic ∂2
−(βF c ) = − = log Zc
∂β 2 ∂β ∂β 2
R R 2
dΓ H2 e−βH dΓ H e−βH
= R − R = hH2 ic − hHi2c
dΓ e−βH dΓ e−βH
34
Le membre de droite est l’écart quadratique moyen (tel que nous l’avons défini au Chapitre
1) dans l’ensemble canonique. La racine carrée de cette expression est “l’écart-type” entre
l’énergie et sa valeur moyenne (aussi défini au Chapitre 1).
Il est instructif de considérer l’écart-type relatif :
q
p ∂
hH ic − hHic
2 2
∂β
hHic 1
= ∝ √
hHic hHic N
Cette propriété résulte du fait que β est intensif, et que hHic est extensif.
Donc pour N → ∞ les systèmes de l’ensemble canonique ont presque
sûrement tous une énergie égale à l’énergie moyenne hHic.
Définissons la capacité calorifique du système (à Volume et Nombre de particules fixés)
comme
∂ 1 ∂
CV ≡ hHic = − 2
hHic > 0
∂T kB T ∂β
cel-00339827, version 1 - 19 Nov 2008
On a donc une relation reliant la capacité calorifique aux fluctuations de l’énergie dans
l’ensemble canonique :
hH2 ic − hHi2c = kB T 2 CV (4.1)
Remarque : Pour N → ∞, les quantités hH2 ic et hHi2c sont ∝ N 2 . Cependant leur

différence est seulement ∝ N .
4.4 Relation entre les ensembles microcanonique et

canonique
La fonction de partition canonique peut se réécrire comme une somme (plus exactement
une intégrale) de fonction de partition microcanonique :
Z Z
1 −βH(Γ)
X
−βE 1
Zc(N, V, T ) = 3N dΓ e = e dΓ
h N! E=δE,2δE...
h3N N ! (E,δE)
| {z }
Zm(N, V, E)
Z
1
= dE e−βE Zm(N, V, E)
δE
Z
1 −βE+ k1 Sm (N,V,E)
= dE e B (4.2)
δE
35
Comme on l’a vu précédemment, on peut s’attendre à ce que Sm(N, V E) soit une

fonction croissante de E. Cependant l’intégrale ci-dessus diverge si Sm augmente plus vite
que kB βE. On admettra que tel n’est pas le cas, et que l’argument de l’exponentielle a un
maximum E = E ∗. La valeur de E ∗ est solution de l’équation

d 1
−βE + Sm(N, V, E) = 0
dE kB
c’est à dire :
1 ∂Sm(N, V, E)
β= ≡ βm(N, V, E)
kB ∂E
d’après (3.5). On a donc
Le maximum de l’intégrand dans (4.2) est atteint pour une valeur E = E ∗ telle que
Tm(N, V, E ∗) = T
où T est la température fixe dans l’ensemble canonique considéré.

cel-00339827, version 1 - 19 Nov 2008
L’intégrand dans (4.2) est très “piqué” autour de E = E ∗ = hEic (comme le suggère le
calcul des fluctuations, quand N devient grand).
On va donc évaluer l’intégrale dans (4.2) en utilisant le méthode de la phase stationnaire (ou
“méthode du col”). Elle est basée sur l’idée que la valeur de l’intégrale sera déterminée,
en bonne approximation quand N devient grand, par un voisinage de l’énergie qui est
d’autant plus petit que l’intégrand est “piqué”.
On utilisera un développement de Taylor autour de E = E ∗, donc on aura besoin de la
dérivée seconde en ce point :

d2 1 1 ∂ 2 Sm
−βE + Sm(N, V, E) = (N, V, E ∗)
dE 2 kB
E=E∗ kB ∂E 2
L’argument de l’exponentielle dans (4.2) est (à un facteur β près)
1 ∂ 2 Sm
−E+T Sm(N, V, E) = −E ∗ + T Sm(N, V, E ∗) + (E−E ∗)2 T (N, V, E ∗)+...
| {z } 2 ∂E 2
≡ −Fm(N, V, E ∗)
Il s’ensuit que
Z ∗ !
1 −βF ∗ 1 ∂ 2 Sm
Zc(N, V, T ) = e m dE exp (E − E ∗)2 + ... (4.3)
δE 2kB ∂E 2
36
On démontrera les faits suivants :

(i) Les termes indiqués en ... sont négligeables
(ii) E ∗ = hEic
(iii) La fonction est très piquée
2S
(iv) ∂∂E 2 = − T 2C
m 1
V
<0
∗
(v) Fc = Fm
On va passer aux variables intensives pour rendre la dépendance dans le grand paramètre
N plus explicite :
Sm est extensive, donc, en posant n = N/V, ǫ = E/N on a
Sm(N, V, E) = N sm(n, ǫ)
∂Sm ∂sm
=
∂E ∂ǫ
∂ 2 Sm 1 ∂ 2 sm
=
cel-00339827, version 1 - 19 Nov 2008
∂E 2 N ∂ǫ2
En posant E ∗ ≡ N ǫ∗, (4.3) est alors équivalente à
Z ∗ !
N ∗ N ∂ 2 sm
Zc(N, V, T ) = e−βFm dǫ exp 2
(ǫ − ǫ∗)2 (1 + O(1/N )) (4.4)
δE 2kB ∂ǫ
L’intégrale dans (4.4) est donc celle d’une Gaussienne, centrée autour de ǫ∗ = hǫic, ce
qui montre les points (ii) et (iii) ci-dessus. L’apparition du facteur N dans la Gaussienne fait
que l’intégrand est très “piqué” autour de ǫ∗ et qu’on est dans les conditions d’application
de la méthode de la phase stationnaire (voir Section suivante pour un rappel de ce résultat
mathématique). Les propriétés de la distribution Gaussienne (4.4) font que la dispersion
de l’énergie est donnée par la formule
kB
hǫ2 ic − hǫi2c = 2
−N ∂∂ǫs2m |∗
Par conséquent on a
N kB
hE 2 ic − hEi2c = 2
− ∂∂ǫs2m |∗
En comparant avec la formule (4.1) on déduit le point (iv) ci-dessus, c’est à dire
∗!−1
N ∂ 2 sm
CV = − >0
T2 ∂ǫ2
37
et donc le coefficient de (ǫ − ǫ∗)2 dans l’exponentielle est bien négatif.

On trouve finalement
s
N ∗
−βFm
2πkB 1 ∗
p
Zc(N, V, T ) = e − ∂ 2 sm
= e−βFm 2πkB T 2 CV (4.5)
δE N ∂ǫ2
δE
Par définition, l’énergie libre canonique Fc est donnée par −kB T log Zc(N, V, T ), et on
a vu (Section 3.2, Exercice 1) que
Fc = hEic − T Sc
∗
D’autre part Fm n’a été introduite que comme abréviation de E ∗ − T Sm(N, V, E ∗)
Prenant donc les logarithmes de (4.5) et divisant par β, on trouve :

1 2πkB T 2 CV
hEic − T Sc = hEic − T Sm(N, V, hEic) − kB T log (4.6)
2 (δE)2
cel-00339827, version 1 - 19 Nov 2008
On va maintenant faire tendre N vers l’infini, c’est à dire passer à la limite thermo-
dynamique. Tous les termes figurant dans (4.6) sont extensifs, sauf le dernier qui est en
log N (puisque la capacité calorifique CV est extensive). Il sera négligeable devant les
autres quand N → ∞.
Remarque 4.4.1 Dans la discussion précédente, on n’a pas tenu compte d’une éventuelle
dépendance de δE en N . δE est arbitraire et a seulement été choisie ≪ E (cf. définition
de l’ensemble microcanonique). La valeur de δE est donc sans importance pourvu qu’elle
soit “sous-extensive”.
Le résultat final est donc le suivant :

À la limite thermodynamique on a :
Sc(N, V, T ) = Sm(N, V, hEic) (4.7)
où hEic est l’énergie moyenne calculée dans l’ensemble canonique.
À partir de là, on trouve facilement d’autres paires de grandeurs thermodynamiques

qui sont égales dans les ensembles microcanonique et canonique par des formules analogues
à (4.7) (toujours à la limite thermodynamique).
38
4.5 La Méthode de la Phase Stationnaire (parenthèse

mathématique)
Problème : Soient deux fonctions g(x) et h(x). Trouver la valeur de l’intégrale
Z b
I = dx h(x) eNg(x)
a
sous forme d’un développement asymptotique pour N → ∞.

Il y a deux cas :
(1) g(x) possède un maximum dans (a, b)
(2) g(x) est maximale pour x = a ou b
On ne traitera que le cas (1). Soit x∗ ∈ (a, b) le point x où g(x) est maximale. On a
donc en utilisant le développement de Taylor de g(x) en x = x∗ :
cel-00339827, version 1 - 19 Nov 2008
39
Z b
1 1
I= dx h(x) exp N g(x∗) + N (x − x∗)2 g ′′ (x∗) + N g ′′′ (x∗)(x − x∗)3
a 2 6
+...)
avec g ′′(x∗) < 0 (car x∗ est un maximum).

√
On fait le changement de variables y = N (x − x∗). Il vient :
√
Z N(b−x∗ )
1 ∗ ∗
y 1 1
I = √ eNg(x ) √
dy h(x + √ ) exp g (x )y + √ g ′′′(x∗)y 3
′′ ∗ 2
N N(a−x∗ ) N 2 6 N
où on a négligé les ordres supérieurs à 4 dans le développement de Taylor dans l’exponen-
tielle car ils sont d’ordre 1/N ou plus.
Par suite :
cel-00339827, version 1 - 19 Nov 2008
√
Z N(b−x∗ )
1 ∗ 1 1 ′′ ′′′∗∗
I = √ eNg(x ) √
dy exp g (x )y 1 + √ g (x )y 3
2
(4.8)
N N(a−x∗ ) 2 6 N

∗ y ′ ∗
× h(x ) + √ h (x ) + ...
N
√ √
Noter que N (a − x∗) → −∞ et N (b − x∗) → +∞ quand N → ∞.
Si on remplace les bornes d’intégration dans (4.8) par ±∞, on commet une erreur qui est

de l’ordre de exp − 12 N (a − x∗)2 |g ′′(x∗)| ou

exp − 21 N (b − x∗)2 |g ′′ (x∗)| .
(4.8) devient alors :
∗) Z
eNg(x +∞
1 ′′ ∗ 2 ∗
1 1
I = √ dy exp − |g (x )|y h(x ) + f1 (y) + 2 f2 (y) + ...
N −∞ 2 N N
+ termes exponentiellement petits en N

où fi sont de polynômes en y 2 que l’on peut déduire facilement, et où les termes en
N −1/2 , N −3/2 , ... disparaissent car ils ne contiennent que des termes impairs en y.
Le résultat final est donc :
s
2π Ng(x∗ ) ∗ 1
I = e h(x ) 1 + O( ) (4.9)
N |g ′′ (x∗)| N
40
4.6 Relation entre les ensembles canonique et grand-

canonique
Le raisonnement est en grande partie analogue. On trouve en remplaçant la somme
sur N par une intégrale
Z ∞ Z
βµN −βFC (N) βµN ∗ −βFc (N ∗ )
∞
(N − N ∗)2 ∂ 2 Fc ∗
Zg (µ, V, T ) ≃ dN e e ≃e dN exp (N )
0 0 2 ∂N 2
où N ∗ est solution de
∂Fc
(N ∗)
µ=
∂N
D’une manière analogue à précédemment on a que
N ∗ = hN ig
À la limite thermodynamique, on a :
cel-00339827, version 1 - 19 Nov 2008
Ωg (µ, V, T ) = Fc(hN ig, V, T ) − µhN ig (4.10)
où hN ig est le nombre moyen de particules calculé dans l’ensemble grand-canonique.

La correspondance mathématique
∂Ωg
Ωg (T, µ, V ) = Fc(N, V, T ) − µN avec N = −
∂µ
est appelée “transformation de Legendre”.
On a
∞ Z
∂ log Zg 1 X Nµ
= βµN
e dΓ(N )e−βH(N,Γ(N))
∂β Zg N=0 N !h3N
∞ Z !
X 1
− eβµN dΓ(N )H(N, ΓN )e−βH(N,γ(N)) = −hEig + µhN ig
N=0
N !h3N
d’où on tire
∂ log Zg ∂(βΩg )
hEig = µhN ig − = µhN ig +
∂β ∂β
en remplaçant log Zg par −βΩg .
∂Ωg
hEig = µhN ig + Ωg + β
∂β
Or Sg = − ∂Ω
∂T
g
et
∂Ωg ∂Ωg
β = −T
∂β ∂T
41
donc
∂Ωg ∂Ωg
Fg = hEig −T Sg = Ωg −T +µhN ig +T = Ωg +µhN ig = Fc(hN ig, V, T )
∂T ∂T
(4.11)
à la limite thermodynamique. Des relations analogues à (4.11) s’obtiennent pour les autres
grandeurs équivalentes à la limite thermodynamique.
Remarque 4.6.1 Les Fonctions de partition Zm, Zc, Zg restent différentes à la limite
thermodynamique.
Dans la Section 4.3, on a calculé la fluctuation de l’énergie dans l’ensemble canonique.

On peut de la même manière étudier les fluctuations du nombre de particules dans
l’ensemble grand-canonique.
Le résultat est la relation
hN 2 ig − hN i2g = hN igT kB nκT (4.12)

cel-00339827, version 1 - 19 Nov 2008
où n = hN ig/V est la densité et où

1 ∂V
κT ≡ − (4.13)
V ∂p T
est la compressibilité isothermique.

On montre aisément que
∂hN ig ∂ 2 log Zg
hN 2 ig − hN i2g = kB T = (kB T )2
∂µ ∂µ2
On utilise la relation dite de “Gibbs-Duheim” exprimant que les variables intensives µ, p, T
d’un système ne sont pas indépendantes :
N dµ − V dp + SdT = 0
À température fixe, il reste simplement N dµ − V dp = 0, ce qui implique que

∂ N ∂
=
∂µ T V ∂p T
Il en résulte (identifiant N de l’ensemble canonique avec hN ig à la limite thermodyna-

mique) que
∂hN i hN i ∂hN i
=
∂µ T V ∂p T
42
Par suite
2 hN i ∂hN i
2 ∂n
hN i − hN i = kB T = hN ikB T (4.14)
V ∂p T ∂p T
Or
1 ∂V 1 ∂n
κT = − =
V ∂p T n ∂p T
d’où en reportant dans (4.14), il vient (4.12).
4.7 La Loi Fondamentale de la thermodynamique et

la formule de Gibbs-Duheim
On verra dans le Chapitre suivant le rôle particulier que joue l’entropie d’un système
dans ses propriétés d’évolution vers l’équilibre. Cela se traduit par une Loi qui exprime
les variations de l’entropie en fonction des paramètres extensifs que sont l’énergie E, le
volume V et le nombre de particules N . Soient T la température, p la pression et µ le
cel-00339827, version 1 - 19 Nov 2008
potentiel chimique qui sont les variables intensives. On a :

Pour une transformation thermodynamique réversible
T dS(N, V, E) = dE + pdV − µdN (4.15)
Or l’extensivité de S implique que
S(λN, λV, λE) = λS(N, V, E)
C’est dire que S est une fonction homogène de degré 1 des trois paramètres extensifs.
Plus généralement :
Définition 4.7.1 f (x, y, z) est dite homogène de degré n si on a
f (λx, λy, λz) = λnf (x, y, z) (4.16)
Dérivant (4.16) par rapport à λ, il vient
∂f ∂f ∂f
x +y +z = nλn−1 f (x, y, z) (4.17)
∂λx ∂λy ∂λz
Par suite en faisant n = λ = 1 dans (4.17), il vient :
∂f ∂f ∂f
x +y +z = f (x, y, z) (4.18)
∂x ∂y ∂z
43
Appliquant (4.18) à S(N, V, E) il vient :

∂S ∂S ∂S
S(N, V, E) = N +V +E
∂N E,V ∂V N,E ∂E N,V
et comme (voir formules pour l’ensemble microcanonique) :

∂S µ ∂S p ∂S 1
=− = =
∂N V,E T ∂V N,E T ∂E N,V T
il vient :
T S(N, V, E) = E + pV − µN (4.19)
Différenciant (4.19) on trouve :
T dS + SdT = dE + pdV + V dp − µdN − N dµ (4.20)

cel-00339827, version 1 - 19 Nov 2008
Comparant avec (4.15) il vient :

Formule de Gibbs Duheim :
SdT − V dp + N dµ = 0
Cette formule exprime que les trois variables intensives que sont T, p, µ ne sont pas
indépendantes.
Chapitre 5
LES PRINCIPES DE LA
THERMODYNAMIQUE
Dans ce Chapitre, on cherchera à déduire les principes de la thermodynamique de la

mécanique statistique.
cel-00339827, version 1 - 19 Nov 2008
5.1 Équilibre entre deux systèmes en contact

5.1.1 Contact thermique
Un récipient de volume V est divisé en deux compartiments de volume V1 et V2 respec-
tivement, par une paroi diatherme c’est à dire qui laisse passer de l’énergie. Les nombres
de particules (resp. les énergies) dans les deux compartiments sont N1 et N2 (resp. E1 et
E2 ).
énergie
< >
N1 V1 E1 N2 V2 E2
L’énergie totale E = E1 + E2 est fixée. On s’attend à ce qu’un équilibre s’établisse

entre les deux compartiments, avec un flux d’énergie nul, où E1 et E2 fluctuent dans
des voisinages étroits de leurs valeurs moyennes hE1 i et hE2 i. On cherche à trouver une
équation pour ces valeurs moyennes. On écrit la fonction de partition microcanonique du
44
45
système entier à l’aide de celles des deux compartiments :

X
Zm(N, V, E) = Zm(N1 , V1 , E1 )Zm(N2 , V2 , E − E1 )
E1 =δE,2δE...
Z
1 1
= dE1 exp [Sm(N1 , V1 , E1 ) + Sm(N2 , V2 , E − E1 )] (5.1)
δE kB
Si les deux compartiments sont très grands, la contribution principale à cette intégrale
provient d’une région étroite autour d’une valeur E1∗ qui sera la moyenne recherchée. Pour
trouver E1∗, on écrit que la “phase est stationnaire” :
d
[Sm(N1 , V1 , E1 ) + Sm(N2 , V2 E − E1 )] = 0
dE1
d’où
∂Sm(N1 , V1 , E1 ) ∂Sm(N2 , V2 , E2 )
=
∂E1 ∂E2
c’est à dire que les tempétatures T1 et T2 sont égales.
cel-00339827, version 1 - 19 Nov 2008
Donc il n’y a plus de flux d’énergie entre les deux compartiments (équilibre thermique)
lorsque leurs températures sont égales, un fait bien connu expérimentalement.
5.1.2 Échange d’énergie et de particules

Le raisonnement est analogue. Ici N = N1 + N2 est fixé, mais N1 et N2 séparément
peuvent fluctuer (de même que les énergies E1 et E2 ). On a
N
X X
Zm(N, V, E) = Zm(N1 , V1 , E1 )Zm(N − N1 , V2 , E − E1 )
N1 =0 E1 =0,δE...
Il y a une double sommation dans le plan (E1 , N1 ). La contribution principale au

résultat, quand N et E sont grands proviendra d’une région étroite autour d’un point
(E1∗, N1∗). On trouve E1∗ comme dans la section précédente, donc on doit avoir
T1 = T2
et N1∗ est solution de l’équation
∂Sm(N1 , V1 , E1 ) ∂Sm(N2 , V2 , E2 )
=
∂N1 ∂N2
c’est à dire
µ1 = µ2
46
Donc il n’y a plus de flux d’énergie ni de particules entre les deux compartiments
lorsque leurs températures et leurs potentiels chimiques sont égaux. C’est comme si
la paroi était absente. En particulier les pressions à droite et à gauche à l’équilibre
sont égales selon notre intuition. Donc p est déterminée par T et µ. (On démontrera
cela plus tard)
5.1.3 Paroi mobile mais adiabatique
<
<
N1 E1 N2 E2
V1 V2
cel-00339827, version 1 - 19 Nov 2008
“Adiabatique” signifie qu’il n’y a pas d’échange de particules ni de chaleur entre les
compartiments. Donc N et E sont fixés pour chacun d’eux. Par un raisonnement analogue
à ceux ci-dessus, on démontre que la paroi ne bouge plus (V1 est fixé à des fluctuations
près) quand les pressions des deux compartiments sont égales.
Exercice : Discuter le cas d’une paroi mobile, diatherme mais imperméable.
5.2 Le second principe de la thermodynamique

tot
Considérons le cas de la sous-section 5.1.1 où la paroi est fixe. L’entropie Sm du système
total est donc, lorsque N1 , V1 , E1 sont imposés :
tot
Sm (N, V, E|N1 , V1 , E1 ) = Sm(N1 , V1 , E1 ) + Sm(N2 , V2 , E − E1 )
Imaginons maintenant l’expérience suivante : on rend tout d’un coup la paroi diathermique.
Le système se trouve alors hors d’équilibre, et part à la recherche d’un nouvel équilibre,
qu’il trouvera quand E1 aura atteint la valeur E1∗. Le calcul de l’intégrale (5.1) montre
que dans ce nouvel état d’équilibre, l’entropie du système est plus grande que dans l’ancien
(puisque E1∗ maximise l’intégrand).
Si ensuite on rend la paroi perméable aux particules, le système se trouve à nouveau hors
d’équilibre. Quand le nouvel équilibre (avec N1 = N1∗ ) est atteint, l’entropie a encore
47
augmenté.
On peut généraliser cette expérience et considérer (ou au moins imaginer) que le grand
système est subdivisé en de nombreux sous-volumes (petits mais macroscopiques) numérotés
par i = 1, 2, ...p, chacun caractérisé par des valeurs Ni, Vi, Ei. Si on rend les parois
successivement (ou simultanément) diathermiques (et/ou perméables) il y aura des flux
d’énergie (et/ou de particules) qui feront augmenter l’entropie jusqu’à ce que le nouvel
équilibre soit atteint.
La même chose se passe pour un système initialement soumis à des contraintes de type
arbitraire (mais au niveau macroscopique) quand on les relaxe.
Avant de formuler le deuxième principe de la thermodynamique, il est utile d’introduire

une nouvelle définition :
cel-00339827, version 1 - 19 Nov 2008
On appellera macro-état l’état d’un système macroscopique spécifié seulement par

des grandeurs macroscopiques, sans distinguer les détails microscopiques.
Exemples :
(1) L’état d’un système spécifié par les trois paramètres indépendants de l’un des ensembles
stationnaires : (N, V, E) ou (N, V, T ) ou (µ, V, T ) . Ces macro-états sont appelés des
états thermodynamiques.
(2) Les états considérés ci-dessus, où il y avait une contrainte sous forme de paroi. Dans
ce cas les deux parties se trouvent chacune à l’équilibre, mais dans des états thermodyna-
miques différents.
(3) L’état d’un système qui se trouve hors d’équilibre (même loin) à condition qu’on puisse
le subdiviser (imaginairement) en petits éléments de volume à l’intérieur desquels il y a
équilibre thermodynamique.
∆V
┌•
Volume ∆V
48
Température T (~
r)
Nombre de particules ∆N (~
r)
∆N(~r)
Densité n(~
r ) = ∆V
S’il s’agit d’un gaz parfait, alors on a localement
p(~
r)∆V = ∆N (~
r )kB T (~
r)
Dans toutes les applications n(~ r) et T (~r ) sont des fonctions qui varient lentement avec
~
r = (~
r1 , ~
r2 , ...~
rN ). On parle d’état thermodynamique local.
Dans tous les cas un macro-état correspond à un nombre très grand de micro-états.
Deuxième principe de la thermodynamique
Si un système thermiquement isolé tend vers son état d’équilibre (à travers une suite
de macro-états) , son entropie augmente jusqu’à ce que l’équilibre soit atteint.
Remarque 5.2.1 Pour arriver à cette formulation du deuxième principe de la thermody-

cel-00339827, version 1 - 19 Nov 2008
namique, il a suffi de comparer des états d’équilibre seulement (bien que contraints). Rien
n’est dit sur la vitesse de la convergence (relaxation) vers l’équilibre.
Remarque 5.2.2 Le terme entropie a été forgé en 1865 par le physicien allemand Clausius
à partir de la racine grecque tropi qui évoque l’idée de transformation ou de retour en
arrière.
( Roger Balian, Les états de la matière, p205, Université de tous les savoirs, Odile Jacob,
2002)
Chapitre 6
APPLICATION AUX
OSCILLATEURS HARMONIQUES
6.1 La fonction de partition canonique de l’oscillateur

harmonique
cel-00339827, version 1 - 19 Nov 2008
6.1.1 Un oscillateur unidimensionnel
1 1
H(x, p) = p2 + mω 2 x2
2m 2
Calculons la fonction de partition canonique :
Z Z s s
1 1 2πm 2π
Zc(T ) = dx dp e−βH(x,p) = (6.1)
h h β βmω 2
2πkB T
=
hω
On en déduit successivement l’énergie libre, l’énergie moyenne et l’entropie :
2πkB T
Fc(T ) = −kB T log Zc(T ) = −kB T log (6.2)
hω

∂ log Zc(T ) ∂ 2π 1
hEic = − =− log = = kB T (6.3)
∂β ∂β βhω β

Fc − hEic 2πkB T
Sc(T ) = − = kB 1 + log
T hω
∂hEic
CV = = kB
∂T
49
50
On peut calculer également l’énergie cinétique moyenne et l’énergie potentielle

moyenne :
p2 1
h ic = kB T
2m 2
1 1
h mω 2 x2 ic = kB T
2 2
2
p 1
hEic = h ic + h mω 2 x2 ic
2m 2
Remarque 6.1.1 On a ici un exemple de ce qu’on appelle l’équipartition de l’énergie
(sous-entendu : entre les degrés de liberté du système).
Plus précisément le théorème de l’équipartition de l’énergie dit qu’une énergie 21 kB T
est stockée dans tout degré de liberté (classique) d’un hamiltonien quadratique. D’où la
relation bien connue
1 3
h m~ v 2 ic = kB T
2 2
cel-00339827, version 1 - 19 Nov 2008
pour l’énergie cinétique d’une particule tri-dimensionnelle.
6.1.2 L’oscillateur tridimensionnel
~2
p 1
H(~
r, p
~) = mω 2~
r2
+
2m 2
Se traite comme 3 oscillateurs hamoniques unidimensionnels indépendants :
- La fonction de partition est le cube de (6.1)
- L’énergie libre et ses dérivées acquièrent un facteur 3 par rapport à (6.2) et (6.3)
6.1.3 N oscillateurs harmoniques (tridimensionnels) identiques

et indépendants
r ≡ (~
~ r1 , ~
r2 , ...~ ~ ≡ (~
rN ) p p1 , p
~2 , ...~
pN )
N N
X 1 X 1
H(~
r, p
~) = ~2i +
p mω 2~
ri2
i=1
2m i=1
2
Se traite comme 3N oscillateurs harmoniques indépendants :
3N
2πkB T
Zc(N, T ) = (6.4)
hω
51
2πkB T
Fc(N, T ) = −3N kB T log (6.5)
hω
hEic = 3N kB T (6.6)
CV = 3N kB
Exercice : Vérifier l’équipartition de l’énergie.
6.1.4 N oscillateurs couplés

Cet exemple est d’une très grande importance pratique.
On considère un cristal unidimensionnel de N atomes couplés identiques. Soit ui la coor-

u ≡ (u1 , u2 , ...uN )
donnée et pi l’impulsion de d’atome i (i = 1, 2, ...N ). On notera ~
~ ≡ (p1 , p2 ...pN ). Si le cristal se trouve dans l’état d’énergie la plus basse, deux atomes
p
voisins se trouveront à une distance a l’un de l’autre. Si les atomes s’écartent de cette
cel-00339827, version 1 - 19 Nov 2008
position, le coût en énergie par rapport à l’énergie la plus basse ne dépend que des écarts
relatifs ui+1 − ui et est donné par
N N
1 X 1 X
H(~
u, p
~) = K (ui+1 − ui)2 + p2i (6.7)
2 i=1
2m i=1
Quelques remarques sur cette expression :

(i) On a défini uN+1 ≡ u1 , c’est à dire que le cristal est une chaı̂ne cyclique d’atomes. On
dit aussi qu’on a choisi des conditions aux limites périodiques.
(ii) K est une constante qui dépend de la nature du cristal.
(iii) La force Fi sur le ième atome est due à ses voisins et est donnée par
∂H
Fi = − = −K(2ui − ui−1 − ui+1 )
∂ui
Fi s’annule donc si et seulement si
ui−1 + ui+1
ui =
2
c’est à dire si le ième atome se trouve exactement au milieu de ses voisins. La force de
rappel est linéaire en ui : on a bien affaire à un cristal harmonique.
Remarque 6.1.2 En fait le fait de choisir des conditions au bord périodiques ne cor-
respond en général à aucune réalité physique, mais a des raisons de commodité mathématique.
52
En particulier elles évitent qu’on ait à s’ocuper des effets de bord. On admet alors impli-
citement que, quelles que soient les conditions aux limites, les propriétés volumiques du
système sont les mêmes pour N grand.
Soit à calculer la fonction de partition canonique pour ce système. Dans ce but, on fait
d’abord un changement de variables (que l’on n’explicitera pas ici) :
(~
u, p ~ ~
~) 7→ (ξ, π) (6.8)
qui diagonalise H :
N N
X πk2 X 1
H= + mωk2 ξk2
k=1
2m k=1
2
où les ωk2 sont différents, et sont les valeurs propres de la transformation.
Exercice : Montrer que pour N = 3, on a ω1 = ω2 = 1 ω3 = 0.
cel-00339827, version 1 - 19 Nov 2008
Le problème est donc celui d’oscillateurs harmoniques indépendants (mais non iden-
tiques). La fonction de partition est :
N Z Z
Y βπ 2 β
Zc(N, T ) = h −N
dξk dπk exp − k − mωk2 ξk2
k=1
2m 2
N
Y 2πkB T
=
k=1
hωk
Donc
N
X 2πkB T
Fc(N, T ) = −kB T log
k=1
hωk
hEic = N kB T indépendant des ωk
CV = N kB toujours indépendant de T
Il est évident (ou presque) que pour un cristal tridimensionnel on aurait trouvé
CV = 3N kB = 3nmN kB = 3nmR
où nm est le nombre de moles dans le système, N est le nombre d’Avogadro, R = N kB

la constante des gaz parfaits.
Donc la capacité calorifique d’un corps cristallin est indépendante de la température
et du corps considéré. C’est la loi de Dulong et Petit (1818).
53
Que disent les expériences ?

Cv
— <
3Nkß régime classique
← aT + bT3 <
T
Pour T ≃ 0, on a CV ≃ AT + BT 3 .
La loi de Dulong est Petit est bien satisfaite au dessus d’une centaine (plusieurs centaines
selon les matériaux) de degrés Kelvin. Les déviations aux basses températures peuvent être
expliquées, comme on le verra, par la mécanique quantique.
cel-00339827, version 1 - 19 Nov 2008
À des températures très élevées, il peut y avoir des déviations dues à l’anharmonicité des
forces entre les atomes : il faur alors tenir compte de termes en (ui+1 − ui)4 négligées
dans (6.7).
6.2 L’oscillateur harmonique quantique unidimension-

nel
En mécanique quantique, un système physique ne possède que des états discrèts (les
états propres du hamiltonien), d’énergie Ei.
Règle générale :
Pour le calcul de la fonction de partition quantique on doit remplacer
Z
1 X
dΓ ρ(H(Γ)) par ρ(Ei)
h3N N ! i
L’hamiltonien est

~2 ∂ 2 1 2 2 † 1
H=− + mω x = ~ω a a +
2m ∂x2 2 2
où a, a† sont les opérateurs d’annihilation et de création habituels, et
~ ≡ h/2π.
54
Niveaux d’énergie :

1
Ei = i + ~ω, i = 0, 1, 2, ...
2
donc
∞ 1
X
−βEi e− 2 β~ω
Zc(T ) = e =
i=0
1 − e−β~ω
1
=
2 sinh( 21 β~ω)
Énergie moyenne :
∂ log Zc(T ) 1 1 + e−β~ω
hEic = − = ~ω
∂β 2 1 − e−β~ω
<
〈
〈E c
cel-00339827, version 1 - 19 Nov 2008
quantique ←
← classique
hω
—
2
<
T
À haute température, la courbe quantique s’approche asymptotiquement de la courbe
classique. À basse température hEic s’approche de l’énergie de l’état fondamental 12 ~ω.
Chaleur spécifique :
∂hEic eβ~ω
CV = = kB (β~ω)2
∂T (eβ~ω − 1)2
Comportement asymptotique :
T → ∞ ou β → 0 : CV ≃ kB

~ω 2 − k~ωT
T → 0 ou β → ∞ : CV ≃ kB e B
kB T
On dit que les degrés de liberté de l’oscillateur harmonique sont “gelés”. C’est un effet
quantique.
55
Cv
< —
kß
1
<
kßT
6.3 Cristal harmonique quantique

L’hamiltonien a la même forme que (6.7) mais ui, pi sont des opérateurs
∂
pi = −i~
∂ui
cel-00339827, version 1 - 19 Nov 2008
Par la même transformation que (6.8) on peut mettre le hamiltonien sous la forme
N
X 1
H= ~ωk(a†kak + )
k=1
2
C’est donc un hamiltonien pour N oscillateurs harmoniques quantiques indépendants.
Niveaux d’énergie :
N
X 1
Ei = (nk + )~ωk
k=1
2
où nk ∈ N, ∀k.
On dit que dans le système il y a nk phonons du type k. Le phonon est un quantum de
vibration dans le cristal quantique. k est le nombre d’onde.
Fonction de partition canonique :

N N
Y e−β~ωk /2 Y 1
Zc(N, T ) = =
k=1
1−e −β~ω k
k=1
2 sinh(β~ωk/2)
Énergie moyenne :
N
1X 1 + e−β~ωk
hEic = ~ωk
2 k=1 1 − e−β~ωk
56
Chaleur spécifique :
N
∂hEic X eβ~ωk
CV = = kB (β~ωk)2
∂T k=1
(eβ~ωk − 1)2
Exercice : Prouver que pour T → ∞, CV (N, T ) tend vers la valeur classique qui est
N kB .
Pour T → 0, les choses sont plus compliquées. On s’attend à ce que pour N grand,
les valeurs des pulsations ωk forment quasiment un continuum. On appelle N (ω)dω le
nombre de pulsations entre ω et ω + dω, et on écrit
Z ∞
eβ~ω
CV (N, T ) = kB dω N (ω)(β~ω)2 β~ω (6.9)
0 (e − 1)2
Z ∞
N = dωN (ω)
0
cel-00339827, version 1 - 19 Nov 2008
Considérons maintenant un cristal réel de N atomes, donc N oscillateurs tridimensionnels

couplés. Le nombre de modes de vibration est alors 3N . La formule (6.9) reste valable,
mais on a maintenant Z ∞
dω N (ω) = 3N
0
Pour T → 0 ce sont les basses fréquences qui seules déterminent CV . On utilise que
N (ω) ≃ C1 ω 2 ω→0 (6.10)
En remplaçant cette formule asymptotique dans (6.9) (approximation dite “de Debye”) il
vient Z ∞
eβ~ω
CV ≃ kB dω C1 ω 2 (β~ω)2
0 (eβ~ω − 1)2
Posant x = β~ω il vient
Z ∞
C1 x4 ex
CV ≃ kB dx = csteT 3
(β~)3 0 (e − 1)
x 2
ce qui explique le terme BT 3 trouvé expérimentalement quand T → 0. On montrera

(Chapitre 11) que le terme AT s’explique non pas par les phonons, mais par le mouvement
des électrons.
57
Remarque 6.3.1 En réalité, la forme (6.10) de N (ω) est une bonne approximation pour
ω < ωD, appelée “pulsation de Debye”, et N (ω) ≃ 0 pour ω > ωD. Il est habituel
d’introduire la “température de Debye” :
kB TD = ~ωD
Les formules ci-dessus sont alors transformées en

Z
T 3 TD /T x4 ex
CV = 9N kB dx x
TD 0 (e − 1)2
cel-00339827, version 1 - 19 Nov 2008
Chapitre 7
SYSTÈMES MAGNÉTIQUES
7.1 Introduction aux systèmes magnétiques

Des moments magnétiques sont associés au cortège électronique des atomes sur un
réseau ( le moment orbital+ le spin seront appelés “spin”).
cel-00339827, version 1 - 19 Nov 2008
(Il existe aussi un magnétisme nucléaire mais qui est beaucoup plus faible)
L’interaction Coulombienne entre les électrons plus le principe de Pauli prédit donc une
interaction d’échange entre spins voisins, caractérisée par une constante de couplage J .
Deux cas sont possibles :
- J > 0 : l’état aligné de deux spins possède l’énergie la plus basse
- J < 0 : l’état opposé (antiparallèle) de deux spins possède l’énergie la plus basse.
Soient i = 1, 2, ...N une numérotation des “spins”, et ~

σi la valeur du ième spin. Pour
x y z
des “spins 1/2”, ~
σ = (σ , σ , σ ) où

x 0 1 y 0 i z 1 0
σ = σ = σ =
1 0 −i 0 0 −1
sont appelées matrices de Pauli.

En ne retenant que les degrés de liberté de “spin” du système (c’est à dire en négligeant
ceux associés aux vibrations des atomes autour de leurs position d’équilibre), on obtient
un hamiltonien dit Hamiltonien de Heisenberg :
X N
X
H = −J σi · σ
~ ~ ·
~ j − gµB H σi
~ (7.1)
hi,ji i=1
où hi, ji indique les paires de “spins voisins”

~ est le champ magnétique extérieur
H
58
59
e~
µB ≡ 2m
est le “magnéton de Bohr”.
Noter que
σj = σixσjx + σiy σjy + σiz σjz
σi · ~
~
Beaucoup de cristaux sont anisotropes (les axes cristallins brisent l’isotropie de l’espace).
Cette anisotropie peut se répercuter sur l’hamiltonien des spins, de telle sorte que les trois
termes σiασjα α = x, y, z acquièrent des coefficients différents.
Dans un cas d’extrême anisotropie, deux de ces coefficients peuvent s’annuler (α = x, y)
et l’hamiltonien des spins devient
X N
X
H = −J σiz σjz − gµB Hz σiz (7.2)
hi,ji i=1
C’est l’hamiltonien d’Ising écrit pour la première fois par Lenz (1920), et étudié par
Ising (1925). Historiquement l’hamiltonien de Heisenberg vient un peu plus tard.
L’espace de phase des spins
cel-00339827, version 1 - 19 Nov 2008
On ne considérera que des spins 12 . L’état du spin est alors complètement spécifié par la
valeur propre de σiz que l’on appellera si et qui vaut ±1. Ses états sont donc |sii et on a :
σiz |sii = si|sii σix|sii = | − sii σiy |sii = −isi| − sii
Un système de N spins possède 2N états linéairement indépendants
|s1 s2 ...sN i = |s1 i ⊗ |s2 i ⊗ ... ⊗ |sN i
On a
σiz |s1 s2 ...sN i = si|s1 s2 ...sN i
Les niveaux d’énergie :
Il sont au nombre de 2N . Il faut résoudre
Hψ = ǫψ
où ψ est une combinaison linéaire des |s1 s2 ...sN i. Pour l’hamiltonien de Heisenberg, c’est
un problème difficile. Pour l’hamiltonien d’Ising on a simplement
H|s1 s2 ...sN i = E(s1 , ...sN )|s1 s2 ...sN i (7.3)

où
X N
X
E(s1 , ...sN ) = −J sisj − gµB H si (7.4)
hi,ji i=1
60
7.2 Mécanique statistique du modèle d’Ising

La fonction de partition canonique est
X
Zc(N, H, T ) = e−βE(s1 ,...sN ) (7.5)
états propres
 
X X X X N
X
= ... exp βJ sisj + βgµB H si 
s1 =±1 s2 =±1 sN =±1 hi,ji i=1
Énergie libre canonique du modèle d’Ising :
Fc(N, H, T ) = −kB T log Zc(N, H, T ) (7.6)
Aimantation totale :
∂Fc(N, H, T )
M(N, H, T ) = − (7.7)
cel-00339827, version 1 - 19 Nov 2008
∂H
Susceptibilité magnétique :
1 ∂M 1 ∂ 2 Fc
χ(N, H, T ) = =− (7.8)
N ∂H N ∂H 2
Soit (s1 , s2 , ..., sN ) une configuration de spins d’énergie E(s1 , ..., sN ). Son aimantation
est
N
∂E(s1 , ..., sN ) X
M =− = gµB si
∂H i=1
On a donc :
N
X
gµB h siic = hM ic
i=1
Or
1 ∂ 1 ∂Zc(N, H, T )
M= log Zc(N, H, T ) ==
β ∂H βZc ∂H
!  
N N N
gµB X X X X X X X
= ... si exp βJ sisj + βgµB H si = gµB h si i c
Zc s1 =±1 s2 =±1 sN =±1 i=1 hi,ji i=1 i=1
Ceci montre donc que l’on a :

N
X
M = hM ic = gµB h si i c
i=1
61
De même on a : !2
N
X
M 2 = (gµB )2 si
i=1
de telle sorte que

 !2 
 XN N
X 
hM 2 ic − hM i2c = (gµB)2 h si i c − h sii2c
 
i=1 i=1
Or !  
N N
∂M ∂ gµB X X X X
= si exp βJ sisj + βgµB H si 
∂H ∂H Zc si =±1 i=1 hi,ji i=1
!2  
N N
1 X X X X
= β(gµB )2 si exp βJ sisj + βgµBH si 
Zc si =±1 i=1 hi,ji i=1
P P 2
cel-00339827, version 1 - 19 Nov 2008
N P PN
 s =±1
i
s
i=1 i exp βJ s s
hi,ji i j + βgµ B H s
i=1 i

−β(gµB )2
 Zc 
 !2 
 XN N
X 
= β(gµB)2 h si i c − h sii2c = β hM 2 ic − hM i2c
 
i=1 i=1
Il en résulte que la susceptibilité magnétique satisfait :
N kB T χ = hM 2 ic − hM i2c
c’est à dire que la susceptibilité magnétique mesure l’écart quadratique moyen de l’ai-
mantation dans l’ensemble canonique. On remarquera l’analogie avec la capacité calorifique
CV qui mesure les fluctuations de l’énergie d’un système.
Remarque 7.2.1 si et g sont sans dimension. µB a la dimension d’un moment magnétique.

62
7.3 Calcul de la fonction de partition canonique d’un

système magnétique simple : N spins sans inter-
action
On prend J = 0. On désignera par h la quantité gµB H (à ne pas confondre avec la
constante de Planck h ! ! !). On a alors :
" #N
X X X
Zc(N, H, T ) = ... eβh(s1 +...sN ) = eβsh
s1 =±1 sN =±1 s=±1
D’où
Zc(N, H, T ) = (2 cosh βh)N

Fc(N, H, T ) = −N kB T log(2 cosh βh)
gµB H
M = gµB N tanh
cel-00339827, version 1 - 19 Nov 2008
kB T
2
(gµB) 2 gµB H
χ(H, T ) = 1 − tanh
kB T kB T
Remarques :
• La susceptibilité magnétique χ est indépendante de N
• L’aimantation M est extensive
• Dans les expressions ci-dessus apparaı̂t le rapport entre deux énergies gµB H et
kB T . Il est donc sans dimension. Noter que l’indice B de kB réfère à Boltzmann, tandis
que celui de µB réfère à Bohr !
‹M c /gµ N
‹
<
>H
Dans cette figure, on voit que quand le champ magnétique est très grand, l’aimantation
atteint sa valeur de saturation gµB N qui exprime qu’alors tous les spins sont parallèles
63
au champ. La pente à l’origine est proportionnelle à la susceptibilité magnétique en champ

nul. Celle-ci est donnée par
g 2 µ2B
χ(0, T ) =
kB T
C’est la loi de Curie. Plus T est basse plus il est facile pour le champ d’aligner les spins.
Remarque 7.3.1 On appelle paramagnétique un système dont l’aimantation augmente

avec le champ magnétique, et vaut zéro en champ nul. Un ensemble de spins libres consti-
tue le paramagnétisme le plus simple. Il est à la théorie du magnétisme ce que le gaz parfait
est à la théorie des gaz.
cel-00339827, version 1 - 19 Nov 2008
Chapitre 8
TRANSITIONS DE PHASE
8.1 Observation expérimentale

8.1.1 Ferromagnétisme
Il existe des matériaux magnétiques qui ont une aimantation non nulle même en champ
cel-00339827, version 1 - 19 Nov 2008
nul ; on l’appelle l’aimantation spontanée notée Ms(T ). Le signe de cette aimantation

dépend alors du signe du dernier champ appliqué. Ms(T ) est non nul à des températures
en dessous d’une température critique T < Tc. Tc est appelée température de
Curie.
Exemples :
Fer (Fe) Tc = 1043◦ K
Nickel (Ni) Tc = 627◦ K
Gadolynium (Gd) Tc = 293◦ K = 20◦ C
64
65
<
<
H
En dessous de Tc la valeur de M(T ) augmente quand T baisse. Une courbe type est la
suivante :
>
cel-00339827, version 1 - 19 Nov 2008
>T
0 Tc
On dit qu’à T = Tc le système magnétique subit une transition de phase. Un
système pouvant subir une telle transition de phase est appelé ferromagnétique.
Pour T > Tc : phase paramagnétique
Pour T < Tc : phase ferromagnétique
8.1.2 Transition liquide-gaz

Losqu’on diminue la température d’un gaz, on voit apparaı̂tre la coexistence de liquide
et de gaz. Il existe une température critique (différente pour chaque gaz !) Tc au dessus de
laquelle il n’existe que la phase gaz ; en dessous de Tc les deux phases vapeur et liquide
sont possibles. À T = Tc le système subit une transition de phase.
Le but des sections suivantes est de répondre à la question :

66
La théorie peut-elle expliquer les transitions de phase ?
8.2 L’approximation de champ moyen pour un ferro-

magnétique
(Théorie de Weiss, 1907)
On ne peut calculer Zc et Fc pour un réseau général. On est donc conduit à introduire

des approximations. L’une des approximations possibles est celle dite du champ moyen.
Celle-ci, bien que non rigoureusement justifiée suffit à expliquer le ferromagnétisme.
On remarque que le spin est soumis à un champ magnétique total qui est la somme du
champ extérieur H et des forces d’échange dues à ses spins voisins. On est donc amené à
remplacer l’hamiltonien
X N
X
cel-00339827, version 1 - 19 Nov 2008
H = −J sisj − gµB H si
hi,ji i=1
par un hamiltonien de spins indépendants

N
X
Hm = − h i si
i=1
Ici hi est un champ effectif agissant sur le spin si, que l’on ne sait pas choisir a priori.
Sachant que le champ réel agissant sur si est donné par
X
J sj + gµB H
j proches voisins de i
on prend, en approximation de champ moyen pour hi la moyenne canonique de cette

expression. Or tous les hsj ic sont égaux :
N
1 X M
hsiic = hsii =
N i=1
N gµB
et donc
M
hi ≡ hcm = h + cJ
N gµB
où on désigne par c le nombre de voisins d’un spin quelconque (dépend de la forme du
cristal).
67
On calcule la fonction de partition de l’hamiltonien ainsi approximé :

X X cm
Zccm = ... eβh (s1 +...+sN )
s1 =±1 sN =±1

M
Fccm = −N kB T log (2 cosh βh cm
) = −N kB T log 2 cosh β(h + cJ )
N gµB

∂Fc M
M=− = gµBN tanh β h + cJ (8.1)
∂H N gµB
(8.1) est une équation d’auto-consistance, puisque M apparaı̂t dans les deux membres.
On pose x = hsiic = gµM BN
. x est compris entre −1 et 1. Il satisfait
x = tanh(βh + βcJ x) (8.2)
et la solution de (8.1) sera M = gµB N x En champ nul l’équation (8.2) devient, en

posant
kB T
θ≡
cel-00339827, version 1 - 19 Nov 2008
cJ
x
x = tanh
θ
^y
1
θ petit
θ grand
>
x
xs
On fait une résolution graphique en traçant sur la même figure la droite y = x et la

courbe y = tanh(x/θ). La pente à l’origine de cette dernière est 1/θ. Clairement pour
θ grand il n’existe qu’une solution x = 0 (d’où M = 0). Pour θ petit, il y en a trois :
x = ±xs, x = 0 (il n’y a pas d’expression analytique de xs). Alors on a Ms = gµB N xs
Il n’est pas difficile de démontrer que si on prend la solution générale M(H, T ), alors
lim M(H, T ) = Ms et lim M(H, T ) = −Ms

H↓0 H↑ 0
68
La transition a lieu pour θ = 1 c’est à dire à la température

cJ
Tc =
kB
On voit que Tc augmente si l’énergie d’échange et/ou le nombre de proches voisins
augmente, et dépend donc du type de matériau considéré.
Voisinage de la température critique : Bien qu’il n’existe pas d’expression analytique
pour la solution xs, on peut toutefois montrer comment l’aimantation spontanée se com-
porte au voisinage de la température critique. Soit T = Tc + ∆T avec ∆T ≪ 1 On
a:
cJ cJ ∆T
= ∆T
≃1−
kB T kB Tc(1 + Tc
) Tc
Pour ∆T = 0, il n’y a que la solution nulle. Pour ∆T petit, xs est nécessairement petit.
On doit résoudre :
∆T ∆T 1 ∆T
cel-00339827, version 1 - 19 Nov 2008
x ≃ tanh x(1 − ) ≃ x(1 − )− x3 (1 − )3 + ...

Tc Tc 3 Tc
c’est à dire
∆T 1
x≃x−x − x3 + termes plus petits
Tc 3
La solution est
3∆T
x2 = −
Tc
Il n’y a pas de solution pour ∆T > 0, comme prévu. Pour ∆T < 0 (c’est à dire un
peu en dessous de la température critique) on a donc
s
T T
Ms(T ) ≃ gµBN 3 1 − 0≤1− ≪1
Tc Tc
expérimental
théorique
>T
69
On appelle équation d’état l’expression de M en fonction de T, H. Cette équation

est représentée par la surface ci-dessous ; chaque point de cette surface représente un état
d’équilibre du système.
^
Tc >T
cel-00339827, version 1 - 19 Nov 2008
La théorie du champ moyen est-elle correcte ? Elle l’est dans les grandes lignes, mais pas
dans les détails. En particulier dans les expériences Mc(T ) au voisinage de T = Tc
s’annule bien avec une pente infinie, mais le comportement est plutôt en

T γ
Ms(T ) ∝ 1 − γ ≈ 0.35
Tc
quel que soit le matériau.
Comportement en champ non nul :
On résoud encore graphiquement l’équation (8.2), en posant M∞ = gµBN . Il vient :
M T gµB
= tanh x = x− H
M∞ Tc kB Tc
T /Tc est donc la pente de la droite y = TTc x − b, tandis que son ordonnée à l’origine est
b = −gµB H/kB Tc. On a donc les deux cas de figure ci-dessous, suivant que T < Tc ou
T > Tc :
70
^y ^y
1 1
>x >x
b
b
-1 -1
Comportement au dessus de Tc :
M(N, T, H) est une fonction continue de H, croissant de façon monotone de 0 à M∞.
Quand H → 0 le point d’intersection de la courbe se rapproche de l’origine. Un développement
limité de tanh x est alors possible : tanh x ∼x→0 x et on trouve :
gµBH
x∼ pour T > Tc, H ≪ kB Tc/µB
cel-00339827, version 1 - 19 Nov 2008
kB (T − Tc)
On en déduit la susceptibilité magnétique au dessus de la température critique :
M 1
χ(N, T, H = 0) ∼H→0 = N (gµB)2
H kB (T − Tc)
On retrouve donc la Loi de Curie pour un modèle d’Ising sans interaction avec le rempla-
cement de T par T − Tc au dénominateur. C’est la Loi de Curie-Weiss.
Comportement au-dessous de Tc :
Deux cas apparaissent suivant que H est grand ou petit. Pour H grand il n’y a qu’un
point d’intersection qui correspond à l’aimantation du système à l’équilibre. Pour H petit
il apparaı̂t 3 points d’intersection, mais seul celui où l’aimantation a le même signe que H
correspond à un équilibre stable. Là encore on peut calculer la susceptibilité magnétique
quand H → 0 (ici au dessous de Tc). On trouve :
N 1 Tc − T
χ∼ (gµB)2 pour 0 < ≪1
2 kB (Tc − T ) Tc
Les exposants critiques :
• Pour l’aimantation :
M ∼T →Tc C|T − TC |1/2
• Pour la susceptibilité :
χ ∼T →Tc C|T − Tc|−1
71
Les exposants critiques prédits par la théorie de champ moyen sont donc universels. Ils
ne dépendent ni de la force de couplage J ni du nombre c de proches voisins. Ils sont les
mêmes pour tous les matériaux.
8.3 Solution exacte de modèle d’Ising

Le modèle d’Ising à une et deux dimensions admet une solution analytique. C’est un
résultat non trivial établi en 1928 (1 d) et en 1944 (2 d).
En dimension 1 on considère une chaı̂ne de spins 1/2 fermée sur elle-même. L’atome du
site N + 1 coincide avec l’atome du site 1. Chaque spin interagit avec ses deux voisins.
L’aimantation qui est le moment magnétique par unité de longueur (L est la longueur de
la chaı̂ne) est donnée par
N gµB sinh gµBH

2kB T
M (H, T ) = 1/2
2L
sinh2 gµBH
+ exp(−βJ )
cel-00339827, version 1 - 19 Nov 2008
2kB T
M (H, T ) est une fonction analytique : elle s’annule pour H = 0 quelle que soit la
température T . La susceptibilité magnétique reste finie à toute température non nulle.
Le modèle d’Ising à une dimension est donc paramagnétique. Il ne présente
pas de transition de phase.
En dimension 2, Onsager a publié en 1944 la détermination analytique de l’énergie

libre en champ nul. Il établit alors la nature ferromagnétique du modèle d’Ising à
2 dimensions. C’est en 1952 que Yang donne le calcul de l’aimantation spontanée de ce
système en fonction de la température. Donnons le résultat pour un réseau carré :
L’énergie moyenne magnétique vaut

NJ J N J sinh2 (J/2kB T ) − 1 2
Ē(H, T ) = − tanh − K(x) − 1
2 2kB T 2 sinh(J/kB T ) π
sinh(J/2kB T )
où x := 2 cosh 2
(J/2k T )
et K(x) est l’intégrale elliptique
B
Z π/2
1
K(x) = dφ
0 (1 − x2 sin2 φ)1/2
Existence d’une transition de phase :
Elle se manifeste par un comportement singulier des fonctions thermodynamiques ou de
leurs dérivées. Le paramètre x est toujours positif. Il s’annule pour T = 0 et T → ∞. Il
passe par un maximum égal à 1 pour une température critique Tc telle que
72
J
sinh =1 (8.3)
2kB Tc
Or K(x) est définie pour tout x < 1 mais diverge pour x = 1 : on a pour x → 1

1 2
T − Tc
K(x) ∼ − log(1 − x ) ∝ − log
2 Tc
Cela conduit pour Ē(N, H, T ) à une fonction continue en T = Tc mais dont la dérivée
(la capacité calorifique) diverge comme log |T − Tc|.
De même on trouve qu’il y a une aimantation en champ nul (caractéristique d’un ferro-
magnétique) donnée par : (s est la surface du cristal)
( h i1/8
NgµB 1
M0 (N, T, H = 0) = 2s
1 − sinh4 (J/2kB T )
pour T < Tc
0 pour T > Tc
M0 décroı̂t donc de sa maleur maximale pour T = 0 qui est N gµB /2s à la valeur 0
cel-00339827, version 1 - 19 Nov 2008
obtenue pout T ≥ Tc. Le système passe alors de sa phase ferromagnétique à la phase

paramagnétique à la temérature Tc définie par (8.3).
Comparaison avec la théorie de champ moyen :
• Température de Curie : Tc = cJ/kB en théorie de champ moyen (c est le nombre
de proches voisins, ici 4), tandis que
Tc = 0, 567...J/kB en résolution exacte.
• Comportement critique : χ(H = 0) ∼ C log(T − Tc) en résolution exacte alors
qu’on avait χ(H = 0) ∼ C(T − Tc)−1 en théorie du champ moyen.
Conclusion :
Quoique difficile soit l’obtention d’une résolution exacte, elle permet de contrôler la validité
des solutions approchées comme ici la théorie de champ moyen. La solution exacte prédit
l’existence d’une transition de phase, laquelle est correctement reproduite, qualitativement
(mais non quantitativement) par la théorie de champ moyen.
Les solutions du modèle d’Ising mettent en évidence l’importance des dimensions du
système dans le problème des transitions de phase, puisque les comportements à une et
deux dimensions sont radicalement différents.
8.4 L’approximation du champ moyen pour un gaz

(la théorie de Van der Waals, 1873)
73
8.4.1 Le problème
L’hamiltonien pour les N particules du gaz est le suivant :
N N X
1 X X
H= ~2i
p + ri − ~
U (~ rj ) (8.4)
2m i=1 i=1 j6=i
On note rij ≡ |~
ri − ~
rj |. On a la factorisation suivante :
3N
1 2πm 2
Zc = Zr
N! βh2
où Zr est la fonction de partition “configurationnelle” :
Z Z !
β XX
Zr = d~
r1 ... rN exp −
d~ U (rij )
2 i j6=i
Clairement, si U = 0, on a Zr = V N , ce qui donne la fonction de partition d’un gaz

cel-00339827, version 1 - 19 Nov 2008
parfait.
Forme de U (r) :
U (r) est dû au déplacement des charges électroniques que deux molécules induisent l’une
sur l’autre lorsqu’elles sont proches. Le calcul est quantique et complexe. Pour r petit,
la fonction U (r) augmente très rapidement (cause physique : les nuages électroniques de
deux molécules ne peuvent pas se recouvrir). Quand r → ∞, U (r) ∼ −r −6 (donc la
force F (r) ∼ r −7 est attractive ; c’est la force empirique dite de Van der Waals.) Au
total on est conduit à supposer une forme pour U (r) qui est la suivante :
" 6 #
r0 12 r0
U (r) = u0 −2
r r
(potentiel de Lennard-Jones), avec u0 ∼ 0.01 eV et r0 de l’ordre de quelques Ang-

ström. Le graphe de U (r) est le suivant :
74
^ U(r)
r0
0
r1 >r
-u0
8.4.2 L’approximation de Van der Waals

On imagine que pour la ième particule, les N − 1 autres représentent un fond prati-
cel-00339827, version 1 - 19 Nov 2008
quement continu qui peut être approximé par un potentiel effectif :

X
U (rij ) ≈ Ueff (~
ri )
j6=i
On a donc 3N/2 Z N
1 2πm −βUef f (~
r)
Zc ≈ d~
re
N ! βh2
Choix de Ueff : Pour Ueff = cste, on retomberait sur le gaz parfait. Ici on tient compte
du fait que dans le volume V il existe pour chaque particule un “volume interdit” v de
taille
|v| = (N − 1)b ≃ N b
où b ∝ r13 est le volume de la “sphère d’exclusion” autour d’une particule. On pose donc

U0 pour ~r ∈V −v
Ueff (~r) =
∞ pour ~r∈v
D’où
3N/2
1 2πm N N
Zc ≈ (V − bN )e−βU0 = Zgp (1 − nb)e−βU0 (8.5)
N! βh2
Valeur de U0 :
U0 est l’énergie potentielle d’une particule due à la présence des autres. Ces autres par-
ticules se trouvent (pratiquement uniquement) dans la partie négative du potentiel de la
75
première. Donc U0 est négatif et proportionnel à la densité n = N/V des autres parti-
cules. On pose :
U0 = −an a > 0, [a] = Jm3
Il vient alors
N
Zc ≈ Zgp (1 − bn)eβan
Soit pour l’énergie libre :
Fc = −kB T log Zc ≈ Fgp − kB N T log(1 − bn) − N an
On a :
Fc = termes qui ne dependent pas de n + N kB T (log n − log(1 − bn)) − N an
d’où
∂Fc N kB T b
= + N kB T − Na
cel-00339827, version 1 - 19 Nov 2008
∂n n 1 − bn
et pour la pression, se rappelant que n = N/V :

∂Fc n2 ∂Fc bn2 nkB T
p=− = = kB T n + − an2 = − an2
∂V N ∂n 1 − bn 1 − bn
Équation d’état de Van der Waals :
(p + an2 )(1 − bn) = nkB T (8.6)
est cubique en n
Exercice : Faire a = b = 0, et retrouver la loi des gaz parfaits p = nkB T .
8.5 La transition de phase gaz-liquide

L’équation d’état de Van der Waals conduit à une surface dans l’espace (p, T, v) (où
v = n−1 ) de la forme suivante :
76
p
haute T
isotherme
critique
T = Tc
instable
T
basse T
Tc
cel-00339827, version 1 - 19 Nov 2008
1
v =
n
Le point critique est donné par l’isotherme T = Tc où

 ∂p
 ∂v = 0
(8.7)
 ∂2 p
∂v2
=0
Résolvant ce système de deux équations on trouve :

8a 1 a
kB Tc = nc = pc = (8.8)
27b 3b 27b2
Noter l’on a la loi universelle (indépendante des valeurs de a, b pour chaque fluide) :
pc 3
= = 0.375
nckB Tc 8
Exercice : Retrouver les formules (8.8). On vérifiera que le système (8.7) conduit aux
équations
v3 2a v4 3a
= =
(v − b) 2 kB T (v − b) 3 kB T
Comparaison avec l’expérience :
77
fluide Tc pc/nckB Tc
N e 44 K 0.31
N2 126 K 0.29
CO2 304 K 0.27
H2 O 647 K 0.23
Remarque 8.5.1 La “loi universelle” est donc approchée pour des fluides réels. C’est la
limitation de la théorie de champ moyen qu’on a utilisée pour démontrer l’équation d’état
de Van der Waals.
cel-00339827, version 1 - 19 Nov 2008
Chapitre 9
MÉCANIQUE STATISTIQUE DE
SYSTÈMES QUANTIQUES
Au Chapitre 7, on a déjà étudié la mécanique statistique d’un système quantique : celui

d’une collection de spins. Les spins sont des systèmes quantiques simplifiés, où le nombre
cel-00339827, version 1 - 19 Nov 2008
de degrés de liberté est fini. En général un système quantique a un nombre infini de

degrés de liberté. Ceci est reflété par le fait que les états quantiques ψ sont des éléments
d’un espace d’Hilbert de dimension infinie.
9.1 Rappels des fondements mathématiques de la mécanique

quantique
L’espace des états quantiques est l’espace H = L2 (Rn) où n est la dimension de
l’espace. C’est un espace de Hilbert avec comme produit scalaire
Z
hψ|ϕi := ψ̄(x)ϕ(x)dx
Rn
C’est le “recouvrement” des deux paquets d’onde représentés par ψ et ϕ. Les fonctions
d’onde ϕ, ψ sont dites orthogonales si elles obéissent la relation d’orthogonalité
hϕ|ψi = 0
c’est à dire si le recouvrement des paquets d’onde correspondants est nul.
Proposition 9.1.1 Dans l’espace de Hilbert H il existe une infinité de bases orthonor-
males, c’est à dire d’ensembles {ϕ}n∈N tels que
hϕn|ϕmi = δn,m
78
79
où δn,m est l’indice de Kronecker égal à 1 si n = m et 0 sinon. Si Id est l’opérateur

identité, on a la relation de complétude :
X
Id = |ϕnihϕn|
n∈N
On a : Z
2
hϕ|ϕi = kϕk = |ϕ(x)|2 dx = 1
Rn
car c’est la probabilité totale de trouver le système quelque part dans Rn :
|ϕ(x)|2 dx
est la probabilité de trouver le système dans un voisinage dx de la position x.
Définition 9.1.2 On appelle opérateur A une application de H dans lui-même qui est
linéaire :
A(λψ + µϕ) = λ(Aψ) + µ(Aϕ)
cel-00339827, version 1 - 19 Nov 2008
REMARQUE : Pour un opérateur A donné et un état ϕ normalisé la quantité Aϕ =

hϕ|Aϕi est appelée valeur moyenne de A dans l’état ϕ. En fait le terme anglais
“expectation value” est plus correct puisque il indique, comme en théorie des probabilités,
l”’ espérance” de la quantité A lorsque le système quantique est dans l’état ϕ. Parfois en
physique on utilise le terme d’observable pour un opérateur pour indiquer que c’est une
quantité qu’on peut observer et mesurer.
Exemples d’opérateurs :
• La position X
• L’impulsion P = −i~∇
• Le moment angulaire L = X ∧ P
• L’énergie H = P 2 /2m + V (X)
On voit que la plupart des opérateurs en mécanique quantique sont non-bornés :

dans un état quelconque ϕ on ne sait pas si ψ ′x = xϕ(x) est de carré intégrable autrement
dit si Z
x2 |ϕ(x)|2 dx < ∞
Rn
ceci n’est vrai que pour des états ϕ suffisamment décroissants à l’infini, appelés états du
domaine de X. De même
∂ϕ
ϕ′ (x) = −i~
∂x1
80
n’est dans L2 (Rn) que pour des états ϕ particuliers dits états du domaine de P1 =
∂
−i~ ∂x 1
tels que

∂ϕ(x) 2
Z
∂x dx < ∞

Rn 1
Définition 9.1.3 Un opérateur A est dit borné si ∀ϕ ∈ H on a Aϕ ∈ H. On a alors
kAk = supϕ∈H:kϕk=1kAϕk < ∞
kAk est appelée norme de A.
Définition 9.1.4 Soit A un opérateur borné et {ϕn}n∈N une base orthonormale de H.

On appelle “élements de matrice” de A les nombres
An,m := hϕn|Aϕmi
Exemples d’opérateurs bornés :

• L’opérateur identité est évidemment borné de norme 1.
cel-00339827, version 1 - 19 Nov 2008
1
•A= 1+X2
est borné. En effet :
Z Z
2
|ϕ(x)|2
kAϕk = dx 2 )2
≤ dx|ϕ(x)|2 = 1
R n (1 + x R n
donc on a kAk ≤ 1.
2 2
• Soit H = P +X 2
l’hamiltonien de l’oscillateur harmonique de masse et de fréquence 1.
Il n’est pas borné ; cependant l’opérateur H est inversible. Soit A = H −1 . A est borné
de norme ≤ 2/~ et ses éléments de matrice dans la base des fonctions d’Hermite sont
1
An,m = 1
δn,m
n+ 2
~
Définition 9.1.5 Produit de deux opérateurs :
Soient A, B deux opérateurs dans H. Le produit A · B est l’opérateur défini comme suit :
A · Bϕ = A(Bϕ)
Si A et B sont bornés le produit A · B est aussi borné et on a
kA · Bk ≤ kAkkBk
Commutateur de deux opérateurs :
[A, B] = A · B − B · A
(doit être défini proprement pour des opérateurs non-bornés.)

81
Exemple de commutateur :
[P, X] = −i~Id
En effet (dimension 1) :
(P · X − X · P )ϕ = −i~[(xϕ)′ − xϕ′] = −i~ϕ(x)
Définition 9.1.6 L’adjoint d’un opérateur A sur H est l’opérateur A∗ satisfaisant
hA∗ϕ|ψi = hϕ|Aψi
pour tout ϕ appartenant au domaine de A et tout ψ tel que
|hϕ|Aψi| ≤ Cψ kϕk
Définition 9.1.7 Un opérateur A est dit symétrique si

cel-00339827, version 1 - 19 Nov 2008
hAϕ|ψi = hϕ|Aψi
pour tous ϕ, ψ appartenant au domaine de A.

Un opérateur est dit autoadjoint si A = A∗. Il est nécessairement symétrique. Cependant
c’est une propriété plus forte que d’être symétrique.
On a cependant
Proposition 9.1.8 Si A est borné et symétrique alors il est autoadjoint.
Par la suite on limitera la notion d’observable à celle d’opérateurs autoadjoints.
Exercice : montrer que pour toute observable A, la valeur moyenne de A dans n’importe
quel état (de son domaine) est réelle.
Définition 9.1.9 Opérateurs à trace :

Soit {ϕn} une base orthonormée de H (par exemple la base des fonctions d’Hermite nor-
malisées :)
2 /2
ϕn(x) = Hn(x)e−x
où les Hn(x) sont les polynômes d’Hermite d’ordre n. Alors un opérateur A est dit “à
trace” si ∞
X
hϕn|(A∗A)1/2 ϕni < ∞
n=0
82
On a alors que sa trace trA est donnée par

∞
X
trA = hϕn|Aϕni
n=0
La trace est indépendante du choix de la base.
Proposition 9.1.10 Un opérateur à trace symétrique peut s’écrire sous la forme

X
ρ= cn|ϕnihϕn|
n
où les ϕn sont les vecteurs d’une base orthonormée de H, et

X
cn := hϕn|ρϕni sont tels que |cn| < ∞
n
Exemples d’opérateurs à trace :

• ∀ϕ ∈ H Aϕ = |ϕihϕ| opérateur de projection sur l’état ϕ,
cel-00339827, version 1 - 19 Nov 2008
on a :
trAϕ = kϕk2 = 1
2 2
• soit H = P +X 2
en dimension 1. C’est l’hamiltonien de l’oscillateur harmonique de
fréquence 1. Soit A = (H 2 + 1)−1 . En prenant pour ϕn les fonctions d’Hermite on a
1
Hϕn = (n + )~
2
d’où
1
hϕn|Aϕni =
1 + (n + 1/2)2 ~2
et la série
X 1
n∈N
1 + (n + 1/2)2 ~2
est convergente.
Définition 9.1.11 Spectre d’une observable :
σ(A) = {λ ∈ C : A − λ n′ est pas inversible}
c’est à dire n’a pas d’inverse borné.
Proposition 9.1.12 Une observable ( c’est à dire un opérateur autoadjoint) a toujours

un spectre réel.
83
Physiquement le spectre d’une observable décrit l’ensemble des valeurs possibles de sa

mesure quantique.
Proposition 9.1.13 On dit que a ∈ R est une valeur propre d’une observable A s’il
existe ϕ ∈ H tel que
Aϕ = aϕ
ϕ est alors appelé état propre de A. Toute valeur propre d’une observable appartient à
son spectre, et les états propres ϕ, ψ correspondant à deux valeurs propres différentes sont
orthogonaux c’est à dire obéissent la relation d’orthogonalité
hϕ|ψi = 0
Si le spectre d’une observable A ne comprend que des valeurs propres, on dit qu’il est
purement ponctuel.
Sinon, le complément du spectre ponctuel de A s’il existe est continu.
cel-00339827, version 1 - 19 Nov 2008
Exemples :
• σ(X) = σ(P ) = R. Ce spectre est continu. Il n’y a pas de valeurs propres dans le
spectre. Interprétation : n’importe quel réel peut être le résultat d’une mesure de la posi-
tion ou de l’impulsion dans un état quantique convenable.
• Énergie cinétique H0 = P 2 /2m. Le spectre est continu et coincide avec [0, +∞[.
• σ(H) = {en = (n + 1/2)~ n ∈ N} où H est l’hamiltonien de l’oscillateur har-
monique en dimension 1 défini plus haut. Ce spectre est purement ponctuel (ou encore
discret). Dans ce cas chaque point du spectre est une valeur propre, c’est à dire que l’on
a:
Hϕn = enϕn
où ϕn est la fonction d’Hermite d’ordre n. Interprétation : les “niveaux d’énergie” de

l’oscillateur harmonique sont quantifiés. Ce sont les résultats possibles des mesures de
l’observable H
De plus les fonctions propres ϕn constituent une base orthonormée de H c’est à dire
qu’on a
hϕn|ϕmi = δn,m
• Spectre de l’atome d’hydrogène :
P2 e2
HH = −
2m |x|
84
HH est autoadjoint sur le domaine des états quantiques qui tendent vers 0 convenablement
à l’origine. Son spectre est donc réel .
Le spectre négatif est discret, constitué des valeurs propres :

me4 1
En = − n = 1, 2, ...
2~2 n2
qui correpondent aux états liés de l’électron avec le proton.
Le spectre positif est continu et correspond aux états dits de diffusion où l’électron est
libre (états ionisés).
0
Systèmes quantiques à N particules :

Soit ψ(~r1 , ~
r2 ...~
rN ; t) la fonction d’onde d’un système à N particules. Elle obéit l’équation
d’évolution de Schrödinger
∂ψ
i~ = Hψ (9.1)
cel-00339827, version 1 - 19 Nov 2008
∂t
ψ est une fonction dans l’espace d’Hilbert L2 (R3N ), c’est à dire qu’elle est de carré
intégrable dans l’espace de toutes les coordonnées ~ r1 , ~
r2 , ...~
rN . On appelle les fonc-
tions ψ les états purs du système quantique. On a pour un système de N particules
de masse m interagissant deux à deux :
N
1 X 1X
H= ~2i +
p ri − ~
U (|~ rj |) (9.2)
2m i=1
2 i6=j
où maintenant p
~i est un opérateur :
∂
~i ≡ −i~∇i = −i~
p
∂~
ri
Le produit scalaire de deux états ψ1 (t) et ψ2 (t) représente le recouvrement des fonctions
d’onde dans tout l’espace ; on désigne par r N ∈ R3N le vecteur (sans flèche)
r N = (~
r1 , ~
r2 , ...~
rN )
Z
hψ1(t) | ψ2 (t)i = dr N ψ1∗(r N , t)ψ2 (r N , t)
R3N
et la normalisation de tout état ψ(t) est donnée par

Z
dr N |ψ1 (r N , t)|2 = 1
R3N
85
Soit ϕn(r N ) une base orthonormale, indépendante du temps. On a
hϕn | ϕmi = δn,m
n, m étant des indices quelconque ; il y en a une infinité. Un état quantique ψ(t) peut
se décomposer dans la base des ϕn :
∞
X
N
ψ(r , t) = cn(t)ϕn(r N )
n=0
où les coefficients cn(t) (complexes) sont donnés par
cn(t) = hϕn | ψ(t)i
et la normalisation de ψ(t) implique que

∞
X
|cn(t)|2 = 1 ∀t ∈ R
n=0
cel-00339827, version 1 - 19 Nov 2008
L’évolution dans le temps des cn(t) est déterminée par l’équation de Schrödinger : en
prenant le produit scalaire de (9.1) avec ϕn (se rappeler que les ϕn sont indépendants du
temps), il vient
∞
dcn(t) X
i~ = Hn,mcm(t)
dt m=0
C’est un système infini d’équations différentielles couplées, donné par la matrice infinie
Hn,m = hϕn | Hϕmi
Soit une observable quantique A. Elle représente une quantité physique mesurable, comme
par exemple la position, l’impulsion... A est donc un opérateur. Sa valeur moyenne Ā
est par définition le résultat de la mesure de A dans l’état quantique ψ :
Ā ≡ hψ | Aψi (9.3)
Supposons que l’on choisisse pour ϕn la base des états propres de A (ce qui suppose que A
soit un opérateur autoadjoint à spectre purement discret, on n’entrera pas dans les détails
mathématiques), alors (9.3) se réduit à
∞ X
X ∞ ∞
X
Ā = Amc∗n(t)cm(t)hϕn | ϕmi = An|cn(t)|2
n=0 m=0 n=0
An est la valeur propre de l’observable A correspondant à la fonction propre ϕn : les An

décrivent les résultats possibles d’une mesure. Les coefficients |cn(t)|2 sont une pondération
purement quantique.
86
Remarque 9.1.14 Dans tout ce qui suit on supposera que H donné par (9.2) n’a que des
valeurs propres discrètes (on dit aussi que son “spectre” est discret) . Cela n’est pas un point
trivial, et cela suppose en particulier que le système est confiné dans un sens convenable
(de même que les systèmes classiques que nous avons considérés étaient confinés dans une
grande boite de volume V ).
9.2 L’opérateur densité

Soit ϕn une base orthonormée d’états propres d’un opérateur A (autoadjoint et à
spectre discret), et An les valeurs propres correpondantes.
Soit ∞
X
k
ψ (t) un état pur quelconque = ckn(t)ϕn (9.4)
n=0
B̄ k = hψ k | Bψ ki la moyenne de B dans ψ k (9.5)

cel-00339827, version 1 - 19 Nov 2008
On considère maintenant une collection d’états purs orthonormés ψ k avec une pondération
pk ≥ 0. dire que les pk sont des “pondérations” revient à supposer
X
pk < ∞
k
On définit l’opérateur densité ρ comme

X
ρ= pk|ψ kihψ k| (9.6)
k
Après normalisation on trouve que

pk
Pk = P = la probabilité classique de trouver le système dans ψ k
p
k k
Bien évidemment
X
Pk = 1
k
et on notera
X
Z≡ pk = trρ
k
En effet
XX ′ ′
X X X
trρ = pkhψ k|ψ k ihψ k |ψ ki = pk δk,k′ = pk
k′ k k k′ k
87
La valeur moyenne de B dans cette collection d’états purs est

P k
k pk B̄
X X
hBi ≡ k
PkB̄ = = Z −1 pkhψk , Bψki
k
Z k
P
Proposition 9.2.1 Tout opérateur densité est un opérateur à trace. Sa trace est k pk .
Remarque 9.2.2 Losqu’un opérateur ρ s’écrit

X
pk|ψ kihψ k|
k
P
avec pk ≥ 0, on a trρ = tr(ρ∗ρ)1/2 = k pk si cette dernière somme est finie.
{ϕn}n∈N étant une base orthonormée on a la relation de complétude suivante :

X
Id = |ϕnihϕn| (9.7)
n∈N
cel-00339827, version 1 - 19 Nov 2008
Soient ρn,m les éléments de matrice de ρ dans la base ϕn. Ils sont donc
X
ρn,m ≡ hϕn , ψ kipkhψ k , ϕmi
k
Par suite on a en utilisant deux fois la relation de complétude (9.7) :

∞
1 X X
k
1 X
hBi = Bm,n pkhψk|ϕmihϕn|ψ i = ρn,mBm,n
Z n,m k
Z n,m=0
Puisque Z = Trρ, il vient :

Tr(ρ · B)
hBi = (9.8)
Trρ
En effet
X X X
tr(ρ · B) = hϕn|ρ · B ϕni = hϕn|ρϕmihϕm|Bϕni = ρn,mBm,n
n n,m n,m
Remarque 9.2.3 (i) Parfois on incorpore Z −1 dans la définition de ρ. L’opérateur den-

sité est alors de trace 1. La définition donnée ici fait ressortir davantage l’analogie avec
les formules pour le cas classique.
(ii) ρ est un opérateur dans l’espace d’Hilbert, comme le sont aussi A, B, H... mais qui
dépend a priori du temps.
88
(iii) La “collection d’états purs” dont on est parti est appelée aussi un mélange statis-
tique d’états purs. On pourrait être tenté de parler d’un “ensemble” d’états purs. Cepen-
dant c’est l’opérateur densité qui est le plus proche de la notion d’ensemble en mécanique
statistique classique. C’est lui qui détermine toutes les moyennes de grandeurs physiques
observables.
Comment est-ce que ρ dépend du temps ? On a :

X
ρm,n(t) = pkckn(t)∗ckm(t)
k
et on se rappelle que l’évolution en temps des cm(t) est donnée par

∞
dckm(t) X
i~ = Hm,n′ ckn′ (t)
dt n′ =0
On trouve donc en utilisant l’hermiticité de H :

cel-00339827, version 1 - 19 Nov 2008
∞
d X
i~ ρm,n(t) = (Hm,n′ ρn′ ,n(t) − ρm,n′ (t)Hn′ ,n)
dt n′ =0
c’est à dire, en notant [A, B] le commutateur des opérateurs A, B :

dρ
i~= [H, ρ] (9.9)
dt
qui est le pendant quantique de l’équation de Liouville
∂ρ
= {H, ρ}
∂t
(crochet de Poisson de ρ et H)
Exercice : Écrire un opérateur densité qui décrit un état pur.
Remarque 9.2.4 Trρ est un réel positif indépendant de la base choisie pour l’évaluer,
et de même bien sûr Tr(ρB) est indépendant de la base (mais pas nécessairement positif ).
9.3 Les opérateurs densité stationnaires

Ils sont solutions de
dρ
i~ = [H, ρ] = 0
dt
89
Une condition suffisante est que ρ soit une fonction de H : ρ = f (H). ρ est alors
diagonal dans la base des fonctions propres de H. Notons ψn et En les états propres et
valeurs propres de H :
Hψn = Enψn
En les utilisant pour base, on trouve

∞
X X
Z = Trρ = hψn | ρψni = ρn,n
n=0 n
P
Tr(ρ · B) n ρn,nhψn | Bψni
hBi = = P
Trρ n ρn,n
Interprétation :
En équilibre thermique, la moyenne d’une observable B s’obtient :
(1) En calculant sa valeur moyenne (quantique) dans les états propres du hamiltonien
(2) En moyennant le résultat avec le poids ρn,n sur tous les états propres.
cel-00339827, version 1 - 19 Nov 2008
Un choix différent de la fonction f conduit à différents opérateurs densité, donc à un

ensemble différent de poids ρn,n. Les ρn,n sont les analogues des “ensembles statistiques
classiques”.
Les ensembles stationnaires quantiques

δm,n si Em ∈ (En, En + δE) ensemble microcanonique
• ρn,m =
0 sinon
• ρn,m = δn,me−βEn ensemble canonique

L’ensemble canonique revient donc à choisir pour opérateur densité
ρ = e−βH
On aura donc

Zc = Tr e−βH

Fc = −kB T log Tr e−βH
• ensemble grand-canonique :
Comme en mécanique statistique classique, il faut sommer sur des systèmes à nombre de
particules N variant de 0 à ∞. Pour cela on dénote HN l’hamiltonien à N corps donné
par (9.2).
90
Fonction de partition grand-canonique :

∞
X ∞
X
Zg (µ, V, T ) = βµN
e Zc(N, V, T ) = eβµN Tr e−βHN (9.10)
N=0 N=0
cel-00339827, version 1 - 19 Nov 2008
Chapitre 10
LES GAZ PARFAITS

QUANTIQUES
Un gaz parfait est un système de particules sans interaction.

cel-00339827, version 1 - 19 Nov 2008
Remarque : (sans rapport avec la mécanique quantique)

Un gaz parfait qui n’est pas à l’équilibre ne peut jamais l’atteindre (l’équilibration étant
due à l’échange d’énergie et de quantité de mouvement lors des chocs entre particules).
Les gaz parfaits n’existent pas : il y a toujours des interactions même faibles, qui garan-
tissent qu’il y aura équilibration.
10.1 L’hamiltonien du gaz parfait et ses fonctions propres
N
1 X
H= ~2j
p
2m j=1
N
! N
~2 X ∂2 ∂2 ∂2 X
=− + + = Hj
2m j=1 ∂x2j ∂yj2 ∂zj2 j=1
Hj est (à un facteur près) le Laplacien sur les trois coordonnées de la position ~
rj :
(xj , yj , zj ) de la jème particule.
On suppose que les particules se trouvent dans une boite cubique de volume V = L3 . On
prend des conditions aux limites périodiques, c’est à dire que l’on identifie les positions
~
r et ~
r + L(m1 , m2 , m3 ) pour tous les entiers m1 , m2 , m3 . L > 0 est la “taille de la
boite”.
91
92
H est une somme sur j d’opérateurs à une particule Hj . Les fonctions propres de Hj sont
1 ~
ψ~k(~
rj ) = eik·~rj
L3/2
avec les valeurs propres correspondantes
~2~
k2
E~k =
2m
où le vecteur d’onde ~

k est quantifié c’est à dire prend les valeurs
~ 2π
k= ~
n avec ~
n = (n1 , n2 , n3 ) (10.1)
L
les n1 , n2 , n3 étant des entiers relatifs arbitraires.
Exercice : Vérifier que les états à une particule ψ~k sont orthonormés.
cel-00339827, version 1 - 19 Nov 2008
Remarque 10.1.1 Grâce à la propriété de conditions aux bords périodiques on est

dans une situation de “confinement” qui assure que le spectre de Hj et donc de H est bien
discret. (voir Remarque 9.1.1)
Notation :
2π
k N
= (~
k1 , ~
k2 , ..., ~
kN ) ∈ Z3N
L
r N = (~
r1 , ~ rN ) ∈ R3N
r2 , ..., ~
σ N = (σ1 , σ2 , ..., σN )
où on désigne par σj le spin de la jème particule.

Valeurs propres et états propres de H :
ΨkN (r N ) = ψ~k1 (~
r1 )ψ~k2 (~
r2 )...ψ~kN (~
rN ) (10.2)
EkN = E~k1 + E~k2 + ... + E~kN
les ~
k1 , ~
k2 , ..., ~
kN sont tous de la forme (10.1) mais sont arbitraires. En particulier ils ne
sont pas nécessairement tous différents (ils peuvent même être tous égaux).
Bien que chaque fonction d’onde du type (10.2) soit mathématiquement fonction propre de
H, elles ne sont pas toutes admises pour la description de systèmes physiques. On appelle
93
σ le spin d’une particule, et on note (~

k, σ) son état quantique. Pour une particule dite
“de spin S”, σ peut prendre les 2S + 1 valeurs
−S, −S + 1, ...S − 1, S
Puisque chaque niveau d’énergie à une particule E~k ne dépend que de ~

k et pas de σ,
chaque état à une particule a une dégénerescence égale à 2S + 1.
Deux cas alors se présentent :
• Bosons : Le spin S est entier
Dans ce cas parmi toutes les fonctions propres de H mathématiquement possibles, seules
sont admises les combinaisons linéaires symétrisées. Soit P une permutation des in-
dices j ∈ {1, 2, ..., N }. La fonction propre bosonique est
X
ΨB N
k ,σ N (r N
) = ψ~k1 ,σ1 (~
rP (1) )ψ~k2 ,σ2 (~
rP (2) )...ψ~kN ,σN (~
rP (N) )
P
(à un facteur de normalisation près qui dans cette discussion est sans importance). L’in-
cel-00339827, version 1 - 19 Nov 2008
dice kN , σ N désigne une collection de nombres quantiques (~kj , σj ) parmi lesquels il peut
y avoir des identiques. Une permutation des indices j correspondants ne conduit pas à un
nouvel état. On dit que les particules obéissent à la statistique de Bose-Einstein.
• Fermions : Le spin S est demi-entier

Seules sont admises les combinaisons liméaires antisymétrisées.
X
ΨkN ,σN (r N ) = (−1)Σ(P ) ψ~k1 ,σ1 (~
rP (1) )ψ~k2 ,σ2 (~
rP (2) )...ψ~kN ,σN (~
rP (N) )
P

ψ~k ,σ (~r1 ) . . . ψ~kN ,σN (~
r1 )
1 1
ψ~ (~
r ) . . . ψ~kN ,σn (~
r2 )
= k1 ,σ1 2 (10.3)
. . . . .

ψ~
k1 ,σ1 (~
rN ) . . . ψ~kN ,σN (~
rN )
Σ(P ) est le nombre d’inversions de la permutation P . (−1)Σ(P ) est la “signature de la

permutation P ”.
(10.3) est le déterminant d’une matrice (complexe) N × N , appelé déterminant de
Slater. Il est clair que l’expression ci-dessus est non nulle seulement si dans la collection
(kN , σ N ) n’apparaisent pas de paires identiques (~ kj , σj ) (car alors le déterminant ci-
dessus aurait 2 colonnes identiques). C’est le principe d’exclusion de Pauli : dans une
collection de fermions, deux particules ne peuvent pas être dans le même état quantique.
On dit que les particules obéissent à la statistique de Fermi-Dirac.
94
10.2 Caractérisation à l’aide des nombres d’occupa-

tion
Pour chaque particule quantique, ~
k vit sur un réseau cubique de maille 2π/L. Chaque
point de ce réseau représente 2S + 1 états. On peut caractériser un microétat d’un système
de N particules quantiques
- soit par la donnée de tous les (~
kj , σj ), j = 1, 2, ..., N
- soit en spécifiant pour chaque (~ k, σ) le nombre n~ de particules dans le système se
k,σ
trouvant dans l’état quantique correspondant ; on appelle n~k,σ un nombre d’occupation.
Dans ce dernier cas, il faut donc donner une liste infinie de nombres d’occupation, qui
sera telle que
X
n~k,σ = N
~
k,σ
Clairement pour la statistique de Bose-Einstein, n~k,σ peut prendre n’importe quelle valeur
entière, alors que pour la statistique de Fermi-Dirac, n~k,σ ne peut valoir que 0 ou 1.
cel-00339827, version 1 - 19 Nov 2008
Remarque importante :
On rencontre souvent en physique quantique de systèmes à N corps de particules indépendantes
mais qui ne constituent pas un gaz parfait.
Exemples :
(1) N particules identiques dans un potentiel extérieur. Dans ce cas, on rajoute au hamil-
tonien (9.2) un potentiel V (~
r1 , ~
r2 , ..., ~
rN ) qui “remplace” le confinement introduit par les
conditions aux limites périodiques.
(2) N particules sans interaction mais dans un potentiel différent (par exemple N oscilla-
teurs harmoniques).
(3) N particules couplées, dont l’hamiltonien, après diagonalisation, prend la forme d’ha-
miltoniens de particules indépendantes (appelés quasiparticules).
Qu’est ce qui change dans ce cas ? Les seules choses à changer sont :
(i) l’étiquette (~
k, σ) d’une fonction propre à une particule est remplacée par un nombre
quantique plus général noté i.
(ii) Les ondes planes sont remplacées par des fonctions d’onde plus générales ψi(~
r), ce
qui sera sans répercution sur la thermodynamique.
95
(iii) l’expression de l’énergie propre individuelle ~2~

k2 /2m qui est caractéristique des
particules libres est remplacée par Ei.
10.3 Calcul de la fonction de partition d’un système

quantique de particules indépendantes
Les formules sont les suivantes :
∗
X
Zc(N, V, T ) = exp(−βEkN ,σN )
(kN ,σN )
où l’astérique dans la somme signifie qu’on tient compte de la statistique liée au spin, soit
celle de Bose-Einstein, soit celle de Fermi-Dirac.
N
X X
E kN ,σN = E~kj ,σj = n~k,σE~k,σ
j=1 ~
k,σ
cel-00339827, version 1 - 19 Nov 2008
On a donc
   
∗
X X ∗N
X X
Zc(N, V, T ) = exp −β n~k′ ,σ′ E~k′ ,σ′  = exp −β n~k′ ,σ′ E~k′ ,σ′ 
{n~k,σ : n~k,σ =N } ~
k′ ,σ′ n~k,σ ~
k′ ,σ′
P
où ∗N signifie la double restriction à la statistique de BE ou FD, et au nombre total de

P
particules égal à N . Cette restriction ~k,σ n~kσ = N est mal commode. Pour l’enlever,
on passe à la fonction de partition grand-canonique :
X∞
Zg (µ, V, T ) = eβµN Zc(N, V, T )
N=0
 
∗
X X
= exp β (µ − E~k′ ,σ′ )n~k′ ,σ′ 
n~k,σ ~
k′ ,σ′
 
Y X∗
=  exp β(µ − E~k,σ)n~k,σ 
~
k,σ n~k,σ
Remarque 10.3.1 Un exemple plus simple de l’opération mathématique effectuée est le

suivant :
Soient a < 0 et b < 0
! !
X X X
exp (an + bm) = ean ebm
n,m∈N n∈N m∈N
que l’on vérifiera à titre d’exercice.

96
Il est nécessaire maintenant de préciser si les particules sont des bosons (statistique de
BE), ou des fermions (statistique de FD).
• Bosons :
X 1
exp β(µ − E~k,σ )n~k,σ = β(µ−E~k,σ )
n~kσ ∈N 1−e
d’où
Yh i−1
ZgBE (µ, V, T ) = 1 − eβ(µ−E~k,σ )
~
kσ
• Fermions :
X
exp β(µ − E~k,σ )n~k,σ = 1 + eβ(µ−E~k,σ )
n~k,σ =0,1
d’où
Yh i
ZgF D(µ, V, T ) = 1 + eβ(µ−E~k,σ )
cel-00339827, version 1 - 19 Nov 2008
~
kσ
On prend le logarithme, et on se limite au cas où E~k,σ ≡ E~k ne dépend que de k :

X
log ZgBE,F D(µ, V, T ) = ∓G log 1 ∓ eβ(µ−E~k ) (10.4)
~
k
où G ≡ 2S + 1
Un état quantique individuel “occupe” en moyenne un “volume de l’espace de phase” égal
à (2π~)3 . Comme le volume accessible dans l’espace des positions est égal à V , cela veut
dire que p~ = ~~
k est quantifié à l’intérieur de cellules élémentaires de plus en plus petites
de taille tridimensionnelle
(2π~)3
V
donc pour V grand, ces états constituent “presque” un continuum, et la somme sur ~
k (se
rappeler que ~
k=p~/~) peut être remplacée par une intégrale :
Z
X V
... → 3
d~
k ...
~
(2π) R 3
k
Pour le gaz parfait, on trouve :

Z
V 2~ 2

log ZgBE,F D(µ, V, T ) = ∓G d~
k log 1 ∓ eβµ−β~ k /2m (10.5)
(2π)3
On utilise la formule de Taylor :
97
∞
x
X elx
∓ log(1 ∓ e ) = (±1)l−1 , x<0 (10.6)
l=1
l
Pour µ < 0 il vient

∞
!
βµl Z 2~ 2
X V e βl~ k
log ZgBE,F D(µ, V, T ) = G (±1)l−1 d~
k exp − (10.7)
(2π)3 l=1 l R 3 2m
Utilisant la variable radiale |~

k| = k et notant que lorsque l’intégrand ne dépend que de
|~
k| on peut remplacer d~ k par 4πk2 dk :
∞ Z ∞
GV X eβµl βl~2 k2
log ZgBE,F D (µ, V, T ) = (±1) l−1 2
k dk exp −
2π 2 l=1
l 0 2m
∞ 3/2 Z ∞
GV X eβµl 2m 2
= (±1) l−1
x2 dxe−x
2π 2 l=1
l βl~2 0
cel-00339827, version 1 - 19 Nov 2008
On utilise Z √
∞
2 −x2
π
x dxe =
0 4
Il vient
∞ 3/2 ∞
X (±1)l−1 2mπ V X eβµl
log ZgBE,F D(µ, V, T ) = GV e βµl
=G (±1)l−1
l=1
l βlh2 λ3 l=1
l5/2
où on a posé
h
λ≡ √
2πmkB T
Thermodynamique :
Comme log Zg est linéaire en V , on a :
log Zg = βpV
où p est la pression. On a donc

∞
BE,F D kB T X eβµl
p (µ, V, T ) = G (±1)l−1 (10.8)
λ3 l=1
l5/2
Noter qu’en fait p ne dépend que des variables intensives µ, T . La densité moyenne :

hN ig ∂p
n≡ = (10.9)
V ∂µ V,T
98
satisfait donc ∞
G X eβµl
n= (±1)l−1 (10.10)
λ3 l=1
l3/2
En éliminant µ entre (10.8) et (10.10) on peut trouver l ’équation d’état du gaz parfait
quantique. Rendons cette élimination explicite dans le cas où n est petit ( V → ∞). On
a un développement limité
eβµ = a1 n + a2 n2 + a3 n3 + ...
Substituant dans (10.10), on obtient :

G 2 a21 n2
n ≃ 3 a1 n + a2 n ± 3/2
λ 2
D’où
λ3 a21
a1 = , a2 = ∓ 3/2
G 2
cel-00339827, version 1 - 19 Nov 2008
Par suite

BE,F D GkB T λ3 2 λ6 2 λ6 λ3 n 1 1
p = n∓n ±n = nkB T 1 ∓ − 5/2
λ3 G G2 23/2 G2 25/2 G 23/2 2
3
BE,F D −5/2 λ n
p = nkB T 1 ∓ 2 + ... (10.11)
G

λ3 n
= nkB T 1 ∓ 0.1768 + ...
G
On notera que nλ3 est bien une quantité sans dimension de même que G.
Exercice : Trouver les ordres supérieurs, en particulier a3 , et que la correction d’ordre n2
à (10.11) est (pour FD comme pour BE)
2
λ3 n
−0.0033
G
Solution :
On doit avoir pour n → 0

G 2 3 1 2 2 1 3
n = 3 a1 n + a2 n + a3 n ± 3/2 (a1 n + a2 n ) + 3/2 (a1 n) + ...
λ 2 3
d’où
λ3 a21
a1 = , a2 = ∓
G 23/2
99

2 a31 1 1
a3 = ∓ a1 a2 − = a31 −√
23/2 33/2 4 27
Or

GkB T 2 3 1 2 2 1
p= a1 n + a2 n + a3 n ± (a1 n + a2 n ) + a3 n3
5/2 1
+ ...
λ3 25/2 3
En reportant les valeurs ci-dessus de a1 , a2 , a3 , il vient

−5/2
nλ3 1 2 nλ3 2
p = nkB T 1 ∓ 2 + ( − 5/2 )( ) + ...
G 8 3 G
Remarque 10.3.2 (i) Les effets quantiques sont petits tant que nλ3 ≪ 1. La particule
exhibe un comportement quantique à l’échelle spatiale λ. Si cette dernière est petite devant
la dimension spatiale d ≡ n−1/3 accessible à une particule, alors le gaz quantique est
bien approximé par le gaz parfait classique. Par contre les effets quantiques deviennent
significatifs quand λ > d ou ≃ d.
cel-00339827, version 1 - 19 Nov 2008
(ii) Le signe ± dans (10.11) correspond à BE/FD. Les fermions se “repoussent” à cause
du principe de Pauli, donc la pression quantique est augmentée par rapport à l’expression
classique. par contre pour les bosons, la pression est plus faible que la pression classique.
(iii) Les effets quantiques sont donc significatifs lorsque :
- la température est faible (quelques degrés Kelvin)
- la densité est forte (liquide, solide)
- les particules sont légères (électrons, Hélium)
(iv)Par contre les effets quantiques sont négligeables quand G → ∞
Si les effets quantiques sont très importants, on parle d’un système quantique fortement
dégénéré. Le développement (10.11) est pour des systèmes quantiques faiblement
dégénérés.
10.4 Moyenne des nombres d’occupation

On a :  
X X
Zg (µ, V, T ) = exp β (µ − E~k′ ,σ′ )n~k′ ,σ′ 
{n~k,σ } ~
k′ ,σ′
et par suite, le nombre moyen de particules est :

1 ∂ log Zg
hN ig =
β ∂µ
100
Une question plus fine concerne le nombre moyen de particules dans un niveau
(~
k0 , σ0 ) donné. Il est donné par
 
1 X X
hn~k0 ,σ0 i ≡ n~k0 ,σ0 exp β (µ − E~k′ ,σ′ )n~k′ ,σ′  (10.12)
Zg
{n~k,σ } ~ ′
k ,σ′
∂ log Zg
=−
∂βEk0 ,σ0
On en déduit :
Les moyennes de nombre d’occupation pour les statistiques de Bose-Einstein et Fermi-
Dirac sont
eβ(µ−E~k,σ )
hn~k,σiBE,F D = (10.13)
1 ∓ eβ(µ−E~k,σ )
Limite classique : On rappelle que µ < 0. D’après (10.10), on a

cel-00339827, version 1 - 19 Nov 2008
eβµ = G−1 nλ3 + ...
donc dans la limite classique (nλ3 ≪ 1) on a aussi eβµ ≪ 1. On obtient alors que (10.13)
tend vers la distribution de Maxwell.
Remarque 10.4.1 La densité totale des particules est donnée par

Z
1 X G eβ(µ−E~k )
n= hn~k,σi BE,F D
= d~
k (10.14)
V ~
8π 3 R3 1 ∓ eβ(µ−E~k )
k,σ
Exercice : Montrer que l’on obtient (10.14) en dérivant (10.5) par rapport à µ et en
divisant par βV . On rappelle que E~k ≡ E~k,σ = ~2~
k2 /2m et que
∂ log Zg
Ng = kB T
∂µ
Chapitre 11
GAZ DE FERMIONS
INDÉPENDANTS
11.1 Introduction
cel-00339827, version 1 - 19 Nov 2008
On considère un gaz de fermions fortement dégénéré à basse température.

(1) Quel est l’état fondamental ?
(2) Quand est-ce que la température peut être considérée comme basse ?
Des exemples sont :
- le “gaz” des électrons de conduction dans un métal
- l’ He3 liquide
- les électrons dans un atome
Électrons dans un métal :
Les électrons ont un spin 1/2. Ce sont donc des fermions.
Soit l la maille d’un réseau cristallin. Elle est typiquement de l’ordre de 3 Angström. D’où
n ≈ l−3 ≈ 4 × 1028 m−3
De plus, on a en unités du Système International :
m = 9 × 10−31 kg
kB = 1.4 × 10−23 JK−1
h = 6.6 × 10−34 Js
On suppose maintenant qu’on est à “température ambiante” :
T = 300◦K
101
102
On a donc
h
λ= √ ≈ 5 × 10−9 m
2πmkB T
et par suite
nλ3 = 4 × 125 × 1028−27 = 6000
Pour que nλ3 ≪ 1 il faudrait que λ3 soit 10 000 fois plus petit, donc λ 20 fois plus petit,
donc T ≈ 120 000◦K. La température ambiante est donc une “très basse température”
pour les électrons dans un métal. Les états des N électrons dans un métal qui vont compter
dans un ensemble canonique seront ceux qui sont très près de l’état fondamental.
11.2 Caractérisation de l’état fondamental

On considère un gaz de N électrons libres et indépendants, dans un cube B de volume
V = L3 . Donc
cel-00339827, version 1 - 19 Nov 2008
~2~
k2 2π
E~k,σ = , ~
k= ~
n, n = (nx, ny , nz ) ∈ Z3
~
2m L
11.2.1 Raisonnement direct

Chaque vecteur d’onde ~ k peut “héberger” deux électrons (fermions) : l’un avec un spin
1
σ = 2 , l’autre avec un spin σ = − 12 .
Les surfaces d’énergie E~ = cste dans l’espace des vecteurs d’onde ~
k k sont des sphères. On
rajoute successivement des électrons en commençant par les énergies les plus basses, pour
minimiser l’énergie totale. Dans un volume donné Ω de l’espace de phase (~ r, p
~), on peut
“loger” Vol(Ω)/h3 états quantiques. Ω est ici l’ensemble des points (~
r, p
~) où
r ∈ B,
~ ~∈S
p
S est la sphère de rayon p = ~k. On a donc
4πp3
Vol(Ω) = V
3
Comme chaque état quantique d’energie E~k peut avoir deux configurations de spin pos-
sibles, le nombre total N d’électrons est donc donné par l’équation
V 4 k3 V
N =2 πp3 =
h3 3 3π 2
103
La solution donne donc un vecteur d’onde maximal, appelé vecteur d’onde de Fermi :
kF = (3π 2 n)1/3 (11.1)
La surface (sphérique) |~ k| = kF est appelée la surface de Fermi. L’électron “ajouté

en dernier”, c’est à dire dont l’énergie est la plus grande pour cet état fondamental à N
électrons a une énergie appelée énergie de Fermi donnée par
~2 kF2
EF =
2m
Cette énergie définit par l’équation E = kB T une température caractéristique appelée
température de Fermi
EF
TF =
kB
Exercice : Montrer que la condition nλ3 ≈ 1 correspond à T ≈ TF .
Solution :
cel-00339827, version 1 - 19 Nov 2008
À la température de Fermi on a
h2
kB TF = (3π 2 n)2/3
8mπ 2
d’où √ !3
8π 2 mkB TF (4π)3/2
3π 2 n = =
h λ3
donc
(4π)3/2 8
nλ3 = 2
= √ ≃ 1.5
3π 3 π
On a donc pour l’état fondamental

(
1 pour |~
k| ≤ kF
n~k,σ = (11.2)
0 pour |~
k| > kF
11.2.2 Raisonnement sur la base des formules précédentes

Comme on est dans le cas de fermions, la formule (10.13) donne
1
hn~k,σi = −β(µ−E~k )
1+e
104
Quand T → 0 on trouve

1 si E~k < µ
hn~k,σ iT =0 =
0 si E~k > µ
On se rappelle que E~k = ~2~

k2 /2m. On a donc la caractérisation (11.2) avec
√
2mµ
kF ≡
~
La condition que doit remplir le potentiel chimique µ est que la densité totale donnée par
(10.14) s’obtient quand tous les niveaux tels que |~
k| < kF sont remplis, tandis que tous
les niveaux tels que |~
k| > kF sont vides :
Z Z
1 1 kF
kF3
n= d~
k~k2 = 4π k2 dk =
4π 3 |~
k|<kF 4π 3
0 3π 2
qui est bien la formule (11.1) donnée dans la sous-section précédente.

cel-00339827, version 1 - 19 Nov 2008
Potentiel chimique à température nulle :
~2
µ(T = 0, n) = (3π 2 n)2/3 (11.3)
2m
Pression à température nulle : On utilise la relation
log Zg = βpV
Z
kB T V
p(T → 0, n) ≃ d~k log 1 + eβ(EF −~2 k2 /2m)
V 4π 3 |~k|<kF | {z T →0
}
≃ β(EF − ~2 k2 /2m)
!
~2~
Z
1 k2 1 ~2 kF2 4 ~2 4
= d~
k EF − = πkF3 − πkF5
4π 3 |~
k|<kF 2m 4π 3 2m 3 2m 5
1 ~2 kF5 2 (3π 2 )2/3 ~2

= = n5/3
nEF =
5 3π 2 m 5 5 2m
Donc même pour T = 0, un gaz parfait de fermions exerce une pression, due aux vitesses
non nulles des particules.
Exercice 1 : Montrer que l’énergie du fondamental E0 est donnée par
3
E0 = N EF
5
105
Solution :
X ~2~
k2 V G~2
Z
E(T = 0) = hn~k,σi|T =0 = d~
k~k2
~
2m 8π 3 · 2m |~
k|≤kF
k,σ
Z
V ~2 kF
V
= 4π dkk4 = EF kF3
8π 3 m 0 5m
Or on a
Z Z kF
X VG ~ V 2 V kF3
N (T = 0) = hn~k,σi|T =0 = dk = 2 k dk =
~
8π 3 |~k|≤kF π 0 3π 2
k,σ
d’où en éliminant kF3 , on trouve le résultat.
Exercice 2 : Montrer que l’on a, à température nulle :

2
cel-00339827, version 1 - 19 Nov 2008
pV = E0
3
On verra que cette relation reste vraie aux températures non nulles.
Solution :
On a vu que
2
p= nEF
5
et comme à température nulle on a :
5 E0
EF =
3N
il vient
2N 5E0 2E0
p= =
5V 3N 3V
11.3 Thermodynamique à très basse température

Application : Contribution des électrons à la chaleur spécifique d’un métal ; on se rap-
pelle que la température ambiante est une “très basse température” pour les électrons d’un
métal.
Méthode : Le système est encore proche de son état fondamental, donc on perturbe au
voisinnage de celui-ci. Quel est CV (N, V, T ) à basse température ?
106
dE
On a CV = dT
, donc il faut trouver E(N, V, T ).
3/2 Z ∞
X V 2m η 3/2
E(µ, V, T ) = hn~k,σ iE~k,σ = dη (11.4)
~
2π 2 ~2 0 1 + eβ(η−µ)
k,σ
où on a posé η = E~k,σ = ~2~

k2 /2m comme variable d’intégration. On calcule de même
N (µ, V, T ) :
3/2 Z ∞
X V 2m η 1/2
N (µ, V, T ) = hn~k,σi = dη (11.5)
~
2π 2 ~2 0 1 + eβ(η−µ)
k,σ
La statégie consiste à inverser (11.5) pour obtenir µ(N, V, T ) et à le reporter dans (11.4) ;
cela se fait, à basse température, à l’aide d’un développement limité en kB T /EF = T /TF .
Remarque :
3/2
EF 3√
cel-00339827, version 1 - 19 Nov 2008
= πn0 λ3
kB T 8
donc kB T /EF est petit quand n0 λ3 est grand.
On calcule, pour
η 1/2
f (η) =
η 3/2
l’intégrale
Z ∞ Z µ Z µ Z ∞
f (η) dη f (η) dη f (η) dη
I = = dη f (η) − + (11.6)
0 1+ eβ(η−µ) 0 0 1+ e−β(η−µ) µ 1 + eβ(η−µ)
On suppose que βµ ≫ 1. Cette hypothèse doit être vérifiée a posteriori, puisque, bien
que µ(T = 0) = EF soit finie, µ(T ) est une fonction non encore connue de T .
Dans le second terme du membre de droite de (11.6), on peut remplacer le borne inférieure
de l’intégrale par −∞ à des contributions exponentiellement petites en eβµ.
Deuxième terme : on pose x = µ − η
Troisième terme : on pose x = η − µ.
On a donc Z µ Z ∞
f (µ + x) − f (µ − x)
I= dη f (η) + dx
0 0 eβx + 1
107
^
eβ(η-μ)+1
∆η
< >
cel-00339827, version 1 - 19 Nov 2008
T<<TF
> η
µ
Les contributions principales à l’intégrale sur x proviennent de x ≤ 1/β ≪ µ, donc
au premier ordre de la série de Taylor, on peut remplacer f (µ + x) − f (µ − x) par
2xf ′(µ) dans l’intégrale. Sachant que
Z ∞
z dz π2
=
0 ez + 1 12
on obtient : Z µ
π2
I= dη f (η) + f ′(µ) + ...
0 6β 2
Pour calculer l’expression de E(µ, V, T ) (resp. N (µ, V, T )) quand T → 0, il suffit de
prendre f (η) = η 3/2 (resp. f (η) = η 1/2 ). On trouve :
3/2 Z
V 2m µ
3/2
3 1/2
π2
E(µ, V, T ) = η dη + µ + ...
2π 2 ~2 0 2 6β 2

V (2mµ)5/2 5 π2
= 1+ × + ...
10π 2 m ~3 2 4β 2 µ2
De même,
3/2 Z
V 2m µ
1/2 π2 1
N (µ, V, T ) = η dη + √ + ...
2π 2 ~2 0 6β 2 2 µ
108
3/2
V 2m 2µ3/2 3 π2
= 1+ × + ...
2π 2 ~2 3 4 6β 2 µ2
D’où " 2 #
V (2mµ)5/2 5π 2 kB T
E(µ, V, T ) = 1+ + ... (11.7)
10π 2 m~3 8 µ
" 2 #
V (2mµ)3/2 π2 kB T
N (µ, V, T ) = 1+ + ... (11.8)
3π 2 ~3 8 µ
L’équation (11.8) donne N en fonction de µ, ou inversement µ en fonction de N (ou de
n = N/V ) :
" 2 #−2/3
~2 π2 kB T
µ(N, V, T ) = (3π 2 n)2/3 1 + + ...
|2m {z } 8 µ
= µ0 ≡ µ(T = 0, n) ≡ EF
ou
cel-00339827, version 1 - 19 Nov 2008
" 2 #
π2 kB T
µ(N, V, T ) = µ0 1 − + ... (11.9)
12 µ0
Il en résulte que " 2 #
π2 T
βµ ≃ βµ0 1 −
12 TF
et par suite
" 2 #
TF π2 T TF π2T
βµ ≃ 1− ≃ − ≫1
T 12 TF T 12TF
comme on devait le vérifier.

La substitution de (11.9) dans (11.7) donne
2 !5/2 2 !
V (2mµ0 )5/2 π2 kB T 5π 2 kB T
E= 1− + ... 1+ + ...
10π 2 m~3 12 µ0 8 µ0
" 2 #
3 5π 2 kB T
E(N, V, T ) = N EF 1 + + ...
5 12 EF
dE
où on a utilisé que EF = µ0 . On calcule de même la chaleur spécifique CV = dT
quand
T → 0.
Récapitulation :
On se souvient que l’énergie du fondamental à T = 0 est :
109
3
E0 = N EF
5
(cf Exercice 1)
Pour T → 0 on a :
" 2 #
π2 kB T
µ(N, V, T ) ≃ µ0 1−
12 µ0
" 2 #
5π 2 kB T
E(N, V, T ) ≃ E0 1+ (11.10)
12 EF
π 2 kB T
CV (N, V, T ) ≃ N kB (11.11)
2 EF
Interprétation et commentaires :
(1) CV (N, V, T ) se comporte linéairement en T pour T ≃ 0.( Le résultat classique était
cel-00339827, version 1 - 19 Nov 2008
CV = 23 N kB ). Pour les phonons d’un solide, on a CV ∝ T 3 quand T → 0. Donc
CVtotal = CVel + CVph ≃ AT + BT 3 quand T → 0
points vides
àT=0
occupé }excitation
pour T > 0
vide
surface ↑ points occupés

de Fermi pour T = 0
(2) Interprétation :
Pour T > 0 des “excitations” sont possibles : un électron juste à l’intérieur de la surface de
Fermi va occuper un niveau à l’extérieur, laissant un “trou”. La surface de Fermi devient
légèrement diffuse.
110
〈
^ 〈n k,σ
→
T=0
1
T = 300° k
>
kF k
Exercice : Déterminer l’ordre de grandeur de ∆k/kF pour un métal à 300◦ K.

Solution : Pour eβ(ǫk −µ) = 10−2 , on a
hn~k,σi ≃ 0.99
Ceci a lieu pour

cel-00339827, version 1 - 19 Nov 2008
~2
β ∆kkF ≃ −4.60 = log 10−2
m
On a alors
∆k 2.3 T
≃− ≃ 2.3
kF β~2 kF2 /2m TF
Or on a vu que TF correspond approximativement à la condition nλ3 = 1. On trouve
donc
h2 n2/3
TF ≃ ≃ 6.5 · 104
2πmkB
◦
Donc à T = 300 K, on a
T
≃ 4.6 · 10−3
TF
ce qui donne pour ∆k/kF :
∆k
≃ 10−2
kF
11.4 Thermodynamique à température quelconque

On calcule le grand-potentiel Ω(µ, V, T ) :
Ω(µ, V, T ) = −kB T log Zg

Z
X
β(µ−E~k,σ )
GV kB T ~
2~
k2
β(µ− ~2m )

= −GkB T log 1 + e ≈− d k log 1 + e
~
8π 3
k
111
Z ∞
GV kB T ~2 k2

=− k2 dk log 1 + eβ(µ− 2m )
2π 2 m 0
2 k2

β(µ− ~2m )
2
On fait une intégration par parties, posant dv = k dk et u = log 1 + e . Le
terme tout intégré est nul. Il vient :
Z ~2 k2 Z
∞ ∞
~2 β 4 eβ(µ− 2m
)
~2 1
Ω = −GkB T V dk k 2 k2 = −GV k4 dk ~2 k2
6mπ 2 0 1+e β(µ− ~2m ) 6mπ 2 0 1 + e−β(µ− 2m
)
On compare avec l’énergie :
~2~
Z ∞
X k2 GV ~2 k4 dk
E= hn~k,σi ≈ ~2 k2
~
2m 2π 2 2m 0 1 + eβ( 2m
−µ)
k,σ
Z ∞
GV ~2 k4 dk
= ~2 k2
4mπ 2 0 eβ( 2m
−µ)
+1
Il vient donc :
cel-00339827, version 1 - 19 Nov 2008
2
Ω(µ, V, T ) = − E(µ, V, T ) (11.12)
3
Pour la pression, comme Ω est linéaire en V :

∂Ω Ω 2E
p=− =− =
∂V V 3V
On a donc la relation générale pour n’importe quelle température :
2E
pV =
3
Chapitre 12
GAZ DE BOSONS
INDÉPENDANTS
12.1 Introduction
cel-00339827, version 1 - 19 Nov 2008
Exemple important : He4 liquide. En effet l’hélium4 étant composé de six particules de
spin 1/2 c’est à dire deux proton, deux neutrons et deux électons est un boson. Contrire-
ment à son isotope l’hélium3 qui est composé de deux protons deux électrons et un neutron
qui est par conséquent un fermion, ayant un nombre impair de fermions. L’hélium4 devient
“superfluide” pour T < Tc = 2.17◦K. On tentera d’expliquer ce phénomène par les pro-
priétés particulières des gaz de bosons à basse température.
Pour considérer un gaz de bosons simple, on prendra S = 0 d’où G = 1.
Dans l’état fondamental, toutes les particules se trouvent à ~

k = 0. Donc l’énergie du fon-
damental est nulle.
Question : Quand T → 0 comment est-ce que le système approche cet état fondamental ?
Rappel : pour des bosons de spin 0 on a :

1
hn~ki =
e−β(µ−ǫ~k ) − 1
Il faut donc que µ < 0 pour qu’aucun des hn~ki ne soit négatif. On voit cependant un
mécanisme mathématique qui permet d’avoir des nombres d’occupation très grands dans
le niveau ~
k=0:
Si β(µ − ǫ~k) ≪ 1 on a :
kB T
hn~ki ≃
~2~
k2
2m
−µ
112
113
avec donc pour ~

k=0:
kB T
hn0 i = −
µ
Ce nombre d’occupation moyen peut devenir très grand pour µ ր 0.
12.2 Consensation de Bose-Eintein

P
Reprenons les calculs du Chapitre 10. Le remplacement de ~
k par
Z
V
d~
k
8π 3
conduit à une densité que l’on notera n1 égale à
Z
1 eβ(µ−ǫ~k )
n1 = d~
k
8π 3 1 − eβ(µ−ǫ~k )
cel-00339827, version 1 - 19 Nov 2008
Z ~2 k2 Z ∞
1 ∞
2 eβ(µ− 2m
)
2m 3/2 ∞ 2 X (βµ−x2 )l
1
= k dk ~2 k2
= x dx e
2π 2 0 1 − eβ(µ− 2m ) β~2 2π 20 l=1
q
2
où on a utilisé comme variable d’intégration x = k β~
2
2m
et que pour y = eβµ−x < 1,
P
(ce qui est le cas pour µ < 0), on a 1−y 1
= ∞ l
l=0 y .
Or Z √
∞
2 −lx2 −3/2
π
x dxe =l
0 4
D’où ∞
3/2 X ∞
m eβµl 1 X eβµl
n1 = =
2πβ~2 l=1
l3/2 λ3 l=1
l3/2
D’où en majorant eβµl par 1 (car µ < 0), ce calcul donne
∞
X 1
n1 ≤ λ−3 = ζ(3/2)λ−3 ≈ 2.612λ−3
l=1
l3/2
où ζ(3/2) est la fonction de Riemann. La densité ne peut donc être égale à cette valeur, que
pout T > T0 où T0 est une valeur critique de la température au dessus de laquelle le calcul
est satisfaisant. C’est celle pour laquelle λ = λ0 (λ est une fonction de la température
∝ T −1/2 ), où
nλ30 = 2.612
ce qui donne pour une densité n donnée :

114
2/3
h2 n
T0 =
2πmkB 2.612
Pour T ≤ T0 l’équation n = n1 ne peut être satisfaite car alors
2.612
n> ≥ n1
λ3
Il se passe que les particules additionnelles qui contribuent à la densité n se placent toutes
au niveau ~ k = 0.
Mais dans ces circonstances les calculs du Chapitre 10 sont à revoir. En effet, considérons
P R
l’expression (10.14) de la densité. Le passage de ~k à 8π1 3 d~k suppose que l’intégrand
varie lentement avec ~ k. Comme ce n’est pas le cas ici pour ~ k ≃ 0, il faut donc traiter
séparément la contribution à la densité de ~
k = 0. Au lieu de (10.14) on a alors
1
n= hn~k=0 i +n1
|V {z }
cel-00339827, version 1 - 19 Nov 2008
= n0 densité des particules à ~

k=0
D’où ∞
3 λ3 1 X eβµl
nλ = + (12.1)
V e−βµ − 1 l=1
l3/2
| {z }
≡ n0 λ3
La formule (10.10) obtenue précédemment n’avait pas le premier terme mais était par
ailleurs identique.
Lorsque µ ր 0 le deuxième terme de (12.1) s’approche de
∞
X
l−3/2 = ζ(3/2) = 2.612...
l=1
où ζ est la fonction de Riemann, donc la limite est finie.

Que fait le premier terme de (12.1) ?
Pour µ ր 0 on a
λ3 1 kB T λ3
≃ −
V e−βµ − 1 Vµ
qui tend donc vers +∞.
Remarque 12.2.1 Dans l’expression ci-dessus, on n’a pas le droit de prendre d’abord la
limite thermodynamique V → ∞ ce qui donnerait zéro ; il faut d’abord obtenir tous les
résultats finaux pour un système fini, et seulement après faire tendre V → ∞.
115
Le premier terme de (12.1) n’est pas négligeable devant le second si

kB T λ3
2.612 ≤
V (e−βµ − 1)
ou encore
−βµ
λ3 kB T
e −1<
2.612V
À la limite thermodynamique le membre de droite est très petit, donc µ doit être très
proche de 0. On remplace alors e−βµ − 1 pas −βµ, et on trouve que µ doit être dans
l’intervalle
λ3 kB T
<µ<0
2.612V
donc c’est un intervalle qui tend vers zéro à la limite thermodynamique.
Quand on augmente le potentiel chimique, la densité augmente. Mais quand µ est très
proche de zéro, toutes ces nouvelles particules vont se placer au niveau ~k = 0. On ap-
pele ce phénomène la condensation de Bose-Einstein . Il se produit en dessous de la
température critique définie précédemment comme fonction de le densité, c’est à dire
cel-00339827, version 1 - 19 Nov 2008
2/3
h2 n
T0 = (12.2)
2πmkB 2.612
3
^ nλ
2,612 MB (nλ3= eβμ)

← 2e terme
1
> eβμ
régime régime (µ = 0)
classique quantique
Dans la limite thermodynamique, il ne reste que la courbe représentant le deuxième

terme, plus un “pic delta” à µ = 0.
Considérons l’expression (12.1) divisée par λ3 , pour n. Pour T > T0 il suffit de prendre
en compte le deuxième terme. Pour T < T0 on a µ → 0 à la limite thermodynamique
donc on peut poser µ = 0 dans le deuxième terme, qui s’écrit alors
3/2
2.612 2.612 λ0 3 T
= = n
λ3 λ30 λ T0
| {z }
116
=n
Donc
 P∞
1 eβ µl


 λ3 l=1 l3/2 pour T > T0
n= 3/2 (12.3)

 − 1 + T
 n pour T < T0
V βµ T0
On a que
1
n0 = −
V βµ
est la densité des états à ~
k = 0.
On voit que la densité des particules à ~
k = 0 est donnée par
" 3/2 #
T 2.612
n0 = 1 − n=n− (V → ∞)
T0 λ3
cel-00339827, version 1 - 19 Nov 2008
Cette formule donne la valeur non nulle de la densité du “condensat de bosons” à ~

k =
0 lorsque T < T0 . C’est une condensation non pas dans l’espace de configuration des
particules mais dans l’espace des vecteurs d’onde.
Si T = 0 toutes les particules sont dans le condensat à ~
k = 0. Si 0 < T < T0 une
fraction macroscopique du nombre total de particules est dans le condensat. Par exemple
pour T = 0.63T0 on a N0 ≈ N 2
. Si T > T0 seule une fraction microscopique de N (donc
~
négligeable) est dans l’état k = 0.
n0
^ /n
>T
T0
Le nombre moyen de particules se trouvant dans une “peau” entre k et k + dk est égal
à 4πk2 hn~ki dk et est donné par la figure suivante :
117
^
4̟k2<nk→>
T<T0
← → T>T0
• >k
12.3 Équation d’état pour un gaz de bosons

On calcule d’abord l’Énergie du système :
Il n’y a pas de contribution des particules condensées à ~
k = 0. On a :
cel-00339827, version 1 - 19 Nov 2008
Z
E 1 X 1 eβ(µ−ǫ~k )
= hn~kiǫ~k = d~
kǫ~k
V V ~
8π 3 1 − eβ(µ−ǫ~k )
k6=0
Z ∞ ∞
1 2 ~2 k2 X βµl βl~2 k2
= k dk e exp −
2π 2 0 2m l=1 2m
∞ Z ∞
~2 X βµl 2m 5/2 ∞ 4 −x2
2m 5/2 3~2 X eβµl
= e x dxe =
4mπ 2 l=1 β~2 l 0 β~2 32mπ 3/2 l=1 l5/2
∞ ∞
3kB T 2πmkBT 3/2 X eβµl 3 kB T X eβµl
= = (12.4)
2 h2 l=1
l 5/2 2 λ 3
l=1
l 5/2
Le calcul est analogue à celui fait précédemment pour (12.1) et utilise que
Z ∞
2 3√
dx x4 e−x = π
0 8
Les formules obtenues au Chapitre 10 pour log ZgBE et pBE restent vraies. Le remplace-
P R
ment de ~k par 8π1 3 d~ k est ici inoffensif car le point ~
k = 0 a un “poids” zéro, donc les
particules dans le condensat ne contribuent pas à la pression. En comparant les formules
(12.4) et (10.8) on voit que
2
pV = E
3
118
Pour T < T0 les deux formules prennent une forme particulière puisque µ = 0 :
∞
2E kB T X kB T kB T
p= = l−5/2 = ζ(5/2) = 1.342 ∝ T 5/2 T < T0
3V λ3 l=1
λ3 λ3
Donc la pression ne dépend pas de la densité en dessous de T0 .

Fixons maintenant la température T . Comment se comporte la pression comme fonction
de la densité n ?
La condensation de Bose (donc la propriété µ = 0) se produit lorsque
2/3
h2 n
T <
2πmkB 2.612
donc en dessus d’une densité critique nc(T ) donnée par
−3/2
h2
nc(T ) = 2.612
2πmkB T
cel-00339827, version 1 - 19 Nov 2008
Pour n < nc la pression augmente quand on augmente la densité. En effet on a

∞
kB T X eβµl
p=
λ3 l=1
l5/2
∞
1 X eβµl
n=
λ3 l=1
l3/2
donc si n augmente, µ augmente, mais alors p augmente.
Mais quand n ≥ nc(T ), le potentiel chimique s’annule et la pression devient alors
constante par rapport à n, et égale à
kB T
p = 1.342
λ3
La courbe “séparatrice” qui sépare la région “normale” et celle de la condensation de Bose
s’obtient en éliminant T entre l’équation
2.612
nc =
λ3
et
kB T
p = 1.341
λ3
C’est donc une courbe, dans le plan (1/n, p) de la forme p = f (1/n) avec f (x) ∝ x−5/3 .
119
Exercice : montrer que pour T = T0 on a :

CV
= 1.925
kB N
Solution :
Pout T ≤ T0 on a µ = 0. Donc la formule (12.4) devient
3kB T 3V kB T
E=V ζ(5/2) = 1, 341
2λ3 2λ3
Or 3/2
−3 n T
λ =
2, 612 T0
On en déduit que pour T ≤ T0 on a (se rappelant que nV = N ) :
1, 341 3N kB T 5/2
E= × (12.5)
2, 612 2 T03/2
cel-00339827, version 1 - 19 Nov 2008
Par suite pour T = T0 on a

∂E 1, 341 15
CV := = × N kB = 1, 925N kB
∂T 2, 612 4
d’où le résultat.
La capacité calorifique est continue et T = T0 , mais sa dérivée est discontinue.
La variation de la capacité calorifique en fonction de la température est la suivante :
^ Cv/NkB
1,925
>T
Application à He4 :
En substituant les valeurs de m et n pour He4 liquide dans (12.2), on trouve :
120
T0 = 3.14◦K
En dessous de sa température critique l’Hélium4 est superfluide.

En réalité, cette transition de phase est observée pour T ≈ 2.17◦K ; il s’agit d’une conden-
sation de Bose-Einstein modifiée par les interactions moléculaires.
Remarque 12.3.1 Il a été démontré que le phénomène de condensation de Bose se pro-

duit également pour des gaz de bosons dilués interagissant par des interactions répulsives.
12.4 Entropie d’un gaz de bosons sans interaction

On a que Z ∞
GV 1
Ω = −kB T k2 dk log ~2 k2
2π 2 0 1 − e−β( 2m
−µ)
On fait une intégration par parties en posant

cel-00339827, version 1 - 19 Nov 2008
1
u = log 2 k2 , k2 dk = dv
−β( ~2m −µ)
1−e
Le terme tout intégré est nul. On a donc
Z ∞
GV ~2 k4 dk
Ω = −kB T 2 k2
6π 2 m 0 e β( ~2m −µ)
−1
Or on a vu que l’énergie est donnée par
Z ∞
V ~2 k4 dk
E(µ, V, T ) = G 2 ~2 k2
2π 2m 0 eβ( 2m
−µ)
−1
d’où on déduit (comme pour des fermions)
2
Ω(µ, V, T ) = − E(µ, V, T )
3
On calcule alors l’entropie qui est égale à
∂F
S=−
∂T
où F (N, V, T ) est l’énergie libre du système donnée par
F = Ω + µN
121
En dessous de la température de Bose, on a µ = 0, donc F = Ω(0, V, T ) et

2 ∂E
S=
3 ∂T N,V
et comme pour T < T0 , on a E ∝ T 5/2

∂E 5E
=
∂T 2T
d’où
3/2
5E T
S= ≈ 1.283N kB pour T < T0 (12.6)
3T T0
Où on a utilisé (12.5). En effet
3/2 3/2
5 1.341 T T
S= × N kB ≈ 1.283N kB
2 2.612 T0 T0
cel-00339827, version 1 - 19 Nov 2008
Remarque 12.4.1 Pour T = 0 l’entropie est nulle. Les particules condensées ne contri-
buent pas à l’entropie. Cela s’explique par le fait que les particules dans le condensat sont
dans un état pur bien déterminé, de sorte que leur contribution au nombre d’états acces-
sibles est simplement un facteur 1. Comme log 1 = 0, la contribution du condensat à
l’entropie est donc nulle.
12.5 Comportement du potentiel chimique en fonc-

tion de N, T
On va établir que le potentiel chimique est continu, avec dérivée première continue en
T = T0 . On a (G=0 pour simplifier)
Z ∞ Z ∞
V k2 dk V k2 dk
N = ~2 k2
= ~2 k2
2π 2 0 eβ( 2m −µ)
−1 2π 2 0 eβ 2m −1
Z " #
∞
V 1 1
+ k2 dk 2 k2 −
2π 2 0 e β( ~2m −µ)
−1 eβ~2 k2 /2m −1
La première intégrale vaut
3/2
T
N
T0
122
En effet elle se réécrit comme

∞ Z ∞ Z
V X ∞ 2 2 k2
−βl ~2m
V X 2m 3/2 ∞ 2 2
2
k dke = 2 2
x dxe−x
2π l=1 0 2π l=1 βl~
|0 {z }
√
π
= 4
3/2 3/2
2πmkB T 2.612 N T
= 2.612V = =N
h2 n λ3 T0
Dans le 2ème terme, il y a divergence de l’intégrale à sa borne inférieure quand µ → 0.
On peut le réécrire comme
Z
V 2m 3/2 ∞ 1/2 1 1
dǫ ǫ − βǫ
4π 2 ~2 0 eβ(ǫ−µ) − 1 e −1
Pour contourner cette difficulté, on divise l’intervalle d’intégration en ǫ entre [0, ηkBT ]
et [ηkB T, ∞[, ou η ≪ 1 est un nombre sans dimension. L’intégrale sur l’intervalle
cel-00339827, version 1 - 19 Nov 2008
[ηkB T, ∞[ ne pose pas de difficulté et est O(µ) quand µ → 0. Dans l’intervalle [0, ηkB T ],
on développe les exponentielles au premier ordre en ǫ et µ. Il vient
Z
V 2m 3/2 µ ηkB T ǫ1/2 dǫ
(12.7)
4π 2 ~2 β 0 ǫ(ǫ − µ)
Mais µ < 0 donc en posant ǫ = x2 , on est ramené à l’intégrale

Z √ηkB T s
2 dx 2 ηkB T
= √ arctan = (−µ)−1/2(π + O(µ))
0 x + |µ|
2 −µ −µ
quand µ ր 0, d’où en reportant dans (12.7) il vient que la contribution de 2ème terme
est
3/2 3/2
V 2m p √ V kB T 2πmkB T0 p
− πkB T −µ = −2 π −µ
4π 2 ~2 (kB T0 )3/2 h2
| {z }
= n/2.612
√
2 π N kB T p N kB T p
=− −µ ≈ −1.357 −µ
2.612 (kB T0 )3/2 (kB T0 )3/2
On en déduit que le nombre total de particules pour T & T0 satisfait :
3/2
T kB T p
N =N − 1.357N −µ
T0 (kB T0 )3/2
123
( Notons que T0 est une fonction de N ). Cette relation peut être inversée pour donner
2 " 3/2 #2
kB T0 T0 T
µ=− −1 pour T & T0
1.842 T T0
On observe bien qu’il y a continuité vers 0 de µ quand T ց T0 . Qu’en est-il de sa dérivée

par rapport à T ?
On trouve
µ(T0 + ∆T ) ∝ (∆T )2
d’où
µ(T0 + ∆T )
ր 0 quand ∆T ց 0
∆T
d’où la continuité de la dérivée première de µ au voisinage de la température de Bose. On
a vu que l’énergie était continûment dérivable en T = T0 (CV est continue). Par contre
on peut montrer que la dérivée seconde de E est discontinue en T = T0 .
cel-00339827, version 1 - 19 Nov 2008
La courbe qui représente le potentiel chimique µ en fonction de la température est donc

comme suit :
^ µ/ kBT
T0
ı >T
Chapitre 13
GAZ DE PHOTONS
13.1 Les équations de Maxwell

→
−
Soit E (~ r , t) le champ électrique. À l’intérieur d’une cavité vide, cubique, d’arête L il
satisfait l’équation des ondes :
cel-00339827, version 1 - 19 Nov 2008
→
−
1 ∂2 E →
−
= ∆ E , ∇~ ·→
−
E =0 (13.1)
c2 ∂t2
c est la vitesse de la lumière.
On cherche les solutions de (13.1) sous la forme d’une “onde progressive” :

→
− →
− ~
r, t) = E0 eik·~r−iωt
E (~
On trouve :
(i) ω et ~
k sont reliés par ω = ω~k ≡ c|~
k|.
(ii) Sous des conditions aux limites périodiques on doit avoir
2π
~
k= n
~ n ∈ Z3
avec ~
L
→
− →
−
(iii) E0 · ~
k = 0, c’est à dire que E0 doit être perpendiculaire à ~ k, conformément aux
→
−
équations de Maxwell. E0 doit donc être une combinaison linéaire de deux vecteurs unitaires
orthogonaux que l’on appellera
ê~k,α α = 1, 2 étant l’indice de polarisation
On a donc trouvé que les ondes progressives sont caractérisées par un vecteur d’onde ~
k et
une polarisation α :
→
− ~
E ~k,α = A~k,αê~k,αeik·~r−iω~k t (13.2)
124
125
où A~k,α est une constante arbitraire, réelle. La solution générale de (13.1) est une super-
position de solutions (13.2) pour toutes les valeurs possibles de ~ k et α.
13.2 Quantification des modes de vibration du champ

électromagnétique
(13.2) s’écrit encore sous la forme
→
− ~
E ~k,α = A~k,α(t)ê~k,αeik·~r
où A~k,α(t) doit satisfaire :

d2
A~k,α(t) + c2 k2 A~k,α(t) = 0
dt2
Cette équation est à comparer avec celle de l’oscillateur harmonique d’hamiltonien
p2 1
cel-00339827, version 1 - 19 Nov 2008
H= + mω 2 x2
2m 2
En effet les solutions classiques x(t) pour un tel hamiltonien satisfont
d2
x + ω 2x = 0
dt2
Le mode de vibration (~ k, α) du champ électromagnétique est donc un oscillateur harmo-
nique de fréquence ω = ck (où k = |~ k|). Ceci établi, la quantification s’effectue selon les
règles de la mécanique quantique :
il s’ensuit que l’énergie correspondant à un mode de vibration (~ k, α) ne peut prendre que
des valeurs discrètes :
1
ǫ~k,α = (n~k,α + )~ωk
2
où n~k,α ∈ N.
On appelle n~k,α le nombre de photons du type (~
k, α). Ce nombre n’est pas restreint ;
les photons sont des bosons.
Remarque 13.2.1 (i) L’indice de polarisation α est analogue au spin σ. Les photons ont
un spin S égal à 1. Pourtant, α ne prend que deux valeurs, et non pas 2S + 1 = 3.
(ii) Il n’y a pas de condensation de Bose-Einstein, car on a toujours µ = 0. On ne peut
fixer le nombre total de photons.
(iii) L’impulsion p
~ du photon est donnée par
~ = ~~
p k
126
Or on a que l’énergie ǫ~k d’un photon de nombre d’onde ~

k est quantifiée par
ǫ~k = ~ω~k = ~c|~ p|

k| = c|~
Comme l’énergie d’une particule relativiste de masse m et d’impulsion p

~ est donnée par
p
ǫ= p ~2 c2 + m2 c4
p| indique que le photon a bien une masse nulle.

on voit que la relation ǫ = c|~
13.3 Le potentiel chimique pour un gaz de photons

Nous considérons le système isolé constitué du gaz de photons et du thermostat. On
applique le postulat fondamental que tous les états possibles d’un système isolé à l’équilibre
sont équiprobables. Les paramètres extérieurs du système sont ici l’énergie et le volume.
Le nombre total de photons ne figure pas parmi les paramètres extérieurs car
cel-00339827, version 1 - 19 Nov 2008
il n’est pas fixé : il y a un constant échange de photons entre la cavité et le thermostat

car l’interaction entre la matière et le rayonnement (quantifié) se fait par absorption et
émission de photons.
La probabilité Pl pour que le gaz de photons se trouve dans l’état microscopique l (d’énergie
El, de nombre de photons Nl...) est proportionnelle au nombre d’états accessibles du
thermostat caractérisés par l’énergie Etot − El, mais ayant un nombre quelconque de
photons :
Pl ∝ ΩR(ER = Etot − El, ...)
L’entropie microcanonique du thermostat est
SR(ER, ...) = kB log ΩR(ER, ...)
On développe cette entropie au voisinage de ER = Etot en se limitant au premier ordre

en El. Les termes d’ordre supérieurs sont négligeables par définition du thermostat. En
effet pour qu’un système soit un thermostat il faut que sa température microcanonique
−1
Tm = ∂S ∂E
m
soit pratiquement indépendante de l’énergie E c’est à dire que
∂Sm ∂Sm ∂ 2 Sm ∂Sm

(ER − E) = (ER) − E 2
≈ (ER)
∂E ∂E ∂E ∂E
ce qui conduit à négliger les termes de dérivée seconde de l’entropie devant le terme d’ordre
un :
El
SR(ER = Etot − El) = SR(ER = Etot) −
TR
127
Comme la température T imposée au gaz de photons est par définition
T = TR
on trouve finalement que

El
−k
Pl ∝ e BT
La probabilité d’un état microscopique l du gaz de photons ne dépend que de l’énergie El

et pas du nombre de photons Nl. Par analogie avec l’ensemble grand-canonique, on voit
que le potentiel chimique µP du gaz de photons est toujours nul à l’équilibre.
µP = 0
Les photons sont des bosons identiques, de masse nulle et de spin 1, sans interaction
mutuelle et de potentiel chimique nul.
Remarque 13.3.1 Si un système de particules a un potentiel chimique nul, cela indique

qu’il a liberté totale d’échanger des particules avec son environnement, ce qui est en effet
cel-00339827, version 1 - 19 Nov 2008
le cas pour les photons dans une cavité. Pour des bosons massifs c’est également le cas
en dessous de la température de Bose, où la température, fixée par un thermostat,
fait en sorte que le système atteigne une certaine densité, en empruntant des bosons à
l’environnement et en les condensant dans l’état ~
k = 0 par le phénomène de condensation
de Bose-Einstein.
Si Eλ est l’énergie d’un état individuel de photon noté λ, le nombre moyen d’occupation
Nλ du niveau Eλ est donné par la distribution de Bose-Einstein pour µ = 0, c’est à dire
1
Nλ =
eβEλ −1
13.4 Le corps noir : un système de photons à l’équilibre

On peut aussi intituler cette section Thermodynamique du rayonnement.
Le corps noir est une cavité dont les parois sont maintenues à une température T . Les
molécules des parois peuvent émettre et absorber des photons jusqu’à ce qu’un équilibre
s’établisse. Cet équilibre sera indépendant de la nature du processus d’émission et d’ab-
sorption. Il est entièrement déterminé par
(i) la température T des parois
(ii) les niveaux d’énergie (modes) du champ électromagnétique à l’intérieur de la cavité.
128
On peut étudier ce rayonnement à l’équilibre en faisant un petit trou (de taille négligeable)
dans l’une des parois, et en mesurant le rayonnement qui sort.
Comme on a µ = 0, le nombre moyen de photons à l’équilibre à la température T et pour
l’état caractérisé par (~
k, α) est :
1
hn~k,αi = (13.3)
eβ~ω~k − 1
Le nombre de niveaux ayant un vecteur d’onde de longueur comprise entre k et k + dk
est égal à
V V ω2
2 3
4πk2 dk = dω ≡ ρ(ω) dω
(2π) π 2 c3
ρ(ω) est appelée la densité spectrale des modes propres. Donc en multipliant par le nombre
moyen d’occupation, on obtient que le nombre n(ω, T ) dω de photons ayant une énergie
entre ~ω et ~(ω + dω) est
V ω2
n(ω, T ) dω = dω
cel-00339827, version 1 - 19 Nov 2008
π 2 c3 eβ~ω − 1
L’énergie dE des photons dans cet intervalle de fréquence est donc
~ ω3
dE = ~ωn(ω, T ) dω = V dω
π 2 c3 (eβ~ω − 1)
| {z }
≡ u(ω, T ) densité spectrale de l’énergie
par unité de volume dans la cavité, quantité que l’on peut mesurer.
Loi de Planck pour le rayonnement d’un corps noir :
~ ω3
u(ω, T ) = (13.4)
π 2 c3 (eβ~ω − 1)
^ u(ω,T)
loi classique
60000k
40000k
30000k
>ω
129
On peut considérer deux cas limites :
• ~ω ≪ kB T :
on a
kB T
u(ω, T ) ≃ ω2 (limite classique ; pas de ~)
π 2 c3
C’est la formule de Rayleigh-Jones.
• ~ω ≫ kB T
on a :
~
u(ω, T ) ≃ ω 3 e−~ω/kB T
π 2 c3
C’est la loi de Wien, découverte empiriquement en 1893.
cel-00339827, version 1 - 19 Nov 2008
Exercice : Montrer qu’à T fixée, u(ω, T ) atteint son maximum pour ω = ωm où
kB T
ωm = 2.821
~
Ce maximum se déplace donc proportionellement à T (loi du déplacement de Wien).
On peut ausi vérifier qu’à une fréquence ω donnée, u(ω, T ) est une fonction croissante de
T.
Solution :La dérivée de u(ω) s’annule pour
x
1− = e−x
3
où x = β~ω. Or une résolution graphique ou numérique à la calculette montre que
x = 2, 821.
130
Énergie totale due au rayonnement :

Z ∞ 4 Z ∞
V~ kB T x3
E(V, T ) = V dω u(ω, T ) = dx
0 π 2 c3 ~ 0 ex − 1
| {z }
= π 4 /15
4 3
π2 kB T 8π 5 kB T
E(V, T ) = V ~c = kB T V (13.5)
15 ~c 15 2π~c
Exercice : En utilisant l’équation aux dimensions, vérifier que l’expression du membre de
droite de (13.5) est bien homogène à une énergie.
(13.5) est dite loi de Stephan-Boltzmann. Elle implique que
E(V, T ) ∝ T 4
cel-00339827, version 1 - 19 Nov 2008
Interprétation :
si ω = kc est la fréquence d’un photon, sa longueur d’onde est donnée par
2π 2πc 2π~c
λth = = =
k ω ǫ~k
donc
2π~c
λth =
kB T
est la longueur d’onde d’un photon ayant une énergie thermique. C’est l’analogue de la
longueur d’onde de de Broglie thermique d’une particule massive.
Grand potentiel :
Z
X eβEλ V
Ω(V, T ) = −kB T log =2 k T
3 B
d~
k log 1 − e−~ω~k /kB T
λ
eβEλ − 1 (2π)
On utilise que
Z Z ∞ Z ∞
V V
2 d~
k... = dωρ(ω)... = dωω 2...
(2π)3 0 π 2 c3 0
131
Il vient : 3 Z ∞
V kB T kB T
Ω(V, T ) = dxx2 log(1 − e−x)
π2 ~c
|0 {z }
R ∞ 3
= − 31 0 dx exx−1
(où on a fait une intégration par parties)
Énergie totale : 3
π2 kB T
E(V, T ) = kB T V
15 ~c
Grand potentiel :
3
π2 kB T 1
Ω(V, T ) = − kB T V = − E(V, T )
45 ~c 3
Entropie :
∂Ω π2 4
kB 4 E
S=− = V 4T 3 = ×
∂T 45 (~c)3 3 T
cel-00339827, version 1 - 19 Nov 2008
Pression du rayonnement :
3
∂Ω π2 kB T
p=− = kB T
∂V 45 ~c
Donc
1
pV = E(V, T )
3
(Rappel : Pour un gaz parfait classique, on avait pV = 32 E)
Remarque 13.4.1 (i) Microscopiquement, la pression est due aux chocs des photons sur
les parois.
(ii) L’entropie d’un gaz de photons s’annule lorsque T → 0, comme celle d’un gaz de
fermions, ou d’un gaz de bosons massifs.
Chapitre 14
SYSTÈMES MOLÉCULAIRES
On a vu qu’on peut traiter les interactions dans des systèmes à N corps en utilisant
l’Approximation du Champ Moyen. Que peut-on dire sur des systèmes moléculaires (typi-
quement des gaz) sans faire cette approximation ?
cel-00339827, version 1 - 19 Nov 2008
14.1 Gaz monoatomiques

On suppose que les atomes sont à symétrie sphérique et que leurs interactions deux à
deux ne dépendent que de leur distance. On a donc un Hamiltonien qui s’écrit
N
X ~2i
p 1X
H= + ri − ~
U (|~ rj |)
i=1
2m 2 i<j
Typiquement U (r) sera un potentiel très répulsif à petite distance, modélisant l’impénétrabilité
des atomes, et attractif à longue distance. Un exemple est le “potentiel de Lenard-Jones” :
" 6 #
r0 12 r0
U (r) = U0 −
r r
N
Quand la densité n = V
des particules n’est pas trop importante, on s’attend à ce que
l’équation d’état du gaz parfait
p = kB T n
soit modifiée par des termes en nk avec k = 2, 3, ... qui seront petits devant le terme
linéaire en n quand n ≪ 1. D’où le développement du viriel :
p
= n + B2 (T )n2 + B3 (T )n3 + ...
kB T
132
133
Le coefficient Bq (T ), q = 2, 3, ... s’appelle le qième coefficient du viriel.

La grande-fonction de partition Zg (µV, T ) s’écrit
∞
X ζN
Zg (µ, V, T ) = ZN
N=0
N!
où
YZ
ZN = e−βU(rij ) dr N
i<j V
et
eβµ
ζ≡
λ3
Dans la limite des gaz parfaits, on a n et eβµ ≈ 0, donc on développe en puissances de ζ.
Ordre 1 :
Zg ≈ 1 + ζZ1 = 1 + ζV
pV
cel-00339827, version 1 - 19 Nov 2008
= log Zg ≈ log(1 + ζV ) ≈ ζV
kB T
∂p/kB T
n=ζ ≈ζ
∂ζ
D’où en éliminant ζ, on retrouve l’équation des gaz parfaits
pV = N kB T
Ordre 2 :
ζ2 ζ2
Zg ≈ 1 + ζZ1 + Z2 = 1 + ζV + Z2
2 2
où on a :
Z Z Z Z
−βU(|~
r1 −~
r2 |)
Z2 = d~
r1 d~
r2 e = d~
r1 d~
r2 (1 + f (r12))
V V V V
Z
2
=V +V d~ r|)
r f (|~
où on a posé
~
r1 + ~
r2
f (r) ≡ e−βU(r) − 1 ~ =
R ~ r1 − ~
r=~ r2
2
Soit Z
a≡ d~
r f (r)
134
On a donc
ζ2
(V 2 + aV )
Zg (µV, T ) ≈ 1 + ζV +
2
pV ζ2 2 ζ 2V 2 1
= log Zg ≈ ζV + (V + aV ) − ≈ ζV + ζ 2 aV
kB T 2 2 2
∂pV /kBT
n=ζ ≈ ζ + aζ 2
∂ζ
En éliminant ζ, on a , au deuxième ordre :
Au deuxième ordre du développement du viriel :
p
≈ n + B2 (T )n2 (14.1)
kB T
avec Z
1
B2 (T ) = − d~
r f (r)
2
cel-00339827, version 1 - 19 Nov 2008
Exemple 1 : 
 ∞ si r < σ
U (r) = −ε si σ < r < bσ

0 si bσ < r
b est un nombre sans dimension > 1.
^ U(r)
σ bσ
>r
D’ où
Z ∞ Z σ Z bσ
2 2 βε
B2 (T ) = −2π dr r f (r) = 2π r dr − 2π(e − 1) r 2 dr
0 0 σ
2π
= σ 3 1 − (b3 − 1)(eβε − 1)
3
135
B2 (T ) est donc une fonction monotone croissante de T qui change de signe à

ε/kB
T = TB =
3

log b3b−1
TB Température de Boyle
• T < TB B2 < 0 =⇒ p < pGP
• T > TB B2 > 0 la répulsion entre les atomes augmente la pression.
Exemple 2 :
U (r) = εϕ(r/σ)
(Le potentiel de Lenard-Jones est de ce type avec ε = U0 et σ = r0 ). On définit une
température réduite
kB T
T∗ ≡
ε
On a alors Z ∞
3 2 ϕ(r)
B2 (T ) = −2πσ r dr exp − ∗ −1
cel-00339827, version 1 - 19 Nov 2008
0 T
Comme dans l’exemple 1, il existe une température de Boyle TB , obtenue lorsque le pa-
ramètre sans dimension T ∗ satisfait
Z ∞
ϕ(r)
2π 1 − exp − ∗ r 2 dr = 0
0 T
Pour le potentiel de Lenard-Jones :
1 1
ϕ(r) = −
r 12 r6
on trouve (numériquement)
TB∗ ≈ 3.418
14.2 Gaz de molécules diatomiques

Dans cette section, on négligera les forces d’interactions entre les molécules. Par contre
on s’intéressera à la structure interne des molécules (supposées diatomiques). donc les N
molécules, outre le terme habituel d’énergie cinétique lié aux translations, ont des degrés de
liberté internes qui interviennent dans l’hamiltonien par un terme HI que l’on n’explicitera
pas dans un premier temps.
Les degrés de liberté internes sont indépendants des degrés de liberté de translation.
N
X ~2j
p
H = H0 + HI = + HI
j=1
2m
136
Par conséquent, un niveau d’énergie d’une molécule donnée est caractérisé :

- par le vecteur d’onde ~
k qui repère l’état propre de H0
- par les nombres quantiques de structure interne désignés collectivement par τ , corres-
pondant à des énergies discrétisées ητ .
La fonction de partition d’ une molécule s’écrit donc :
" !#
~2~
Z
V X k2
z= d~k exp −β + ητ
(2π)3 τ
2m
L’indépendance des deux types de degrés de liberté conduit à une factorisation :
z = z0 ζ
où z0 est la fonction de partition de translation et ζ(β) la fonction de partition associée à
la structure interne de la molécule :
3/2
m
z0 = V ≡ V λ−3
2π~2 β
cel-00339827, version 1 - 19 Nov 2008
et
X
ζ(β) ≡ e−βητ (14.2)
τ
Finalement on a l’expression générale pour la fonction de partition d’un gaz parfait
de molécules :
V N −3N
Z= λ [ζ(β)]N
N!
14.2.1 Le phénomène de “gel” des degrés de liberté internes

Soit η0 la plus basse des énergies ητ possibles (niveau fondamental de structure interne),
et η1 l’énergie immédiatement supérieure dans le spectre. Soit g0 la dégénérescence de η0 .
Si kB T ≪ η1 − η0 , tous les termes de la somme dans (14.2) sont négligeables devant les
g0 premiers, de telle sorte que la fonction de partition interne prend la forme
ζ(β) ≃ g0 e−βη0
On dit alors que les degrés de liberté internes sont “gelés”. Ce phénomène est d’origine
purement quantique : c’est l’existence d’une énergie minimale d’excitation η1 − η0 qui en
est la cause. Les températures pour lesquelles se manifeste l’effet de “gel” ne sont pas a priori
impossibles à atteindre. Nous calculerons leur ordre de grandeur à la section suivante pour,
respectivement, les degrés de liberté de rotation et de vibration d’une molécule diatomique.
Un mode d’excitation interne d’une molécule ne se manifeste que si kB T est compa-
rable ou supérieur à l’énergie d’excitation minimale de ce mode.
137
14.2.2 Structure interne des molécules diatomiques

L‘Hamiltonien HI est la somme de trois termes indépendants les uns des autres :
HI = He + Hr + Hv
He est l’Hamiltonien du cortège électronique de la molécule, Hr décrit ses rotations dans

l’espace, tandis que Hv décrit les vibrations des ceux atomes composant la molécule.
• degrés de liberté électroniques :
Les valeurs typiques des énergies d’excitation d’un cortège électronique atomique ou moléculaire
sont de l’ordre de quelques électrons-volts, ce qui correspond à des
températures caractéristiques de
kB T ≃ 1eV ←→ T ≃ 12000◦K
Il faut cependant tenir compte de la “structure fine” du niveau fondamental de la molécule :

cel-00339827, version 1 - 19 Nov 2008
les degrés de liberté électroniques sont bien gelés mais avec un facteur de dégénérescence
égal à g0 = (2L + 1)(2S + 1), où L est le moment cinétique orbital, et S le spin total.
• degrés de liberté de vibration :

Dans l’“approximation de Born-Oppenheimer” (voir Cours de Mécanique Quantique), on
considère que les noyaux (lourds par rapport à l’électron) sont fixes, à une distance ρ l’un
de l’autre, et on calcule l’énergie minimale du système constitué par les deux noyaux et les
électrons comme une fonction u(ρ), présentant un minimum prononcé pour ρ = d où d
est la distance d’équilibre entre les deux noyaux.
u(ρ)
^
d
>ρ
Mais les noyaux vibrent l’une par rapport à l’autre, ce qui signifie que leur distance ρ
oscille autour de d. Cette vibration est équivalente au mouvement d’une particule fictive
138
de masse
m1 m2
mR ≡
m1 + m2
(m1 et m2 sont les masses des deux noyaux)
On a donc affaire à un oscillateur hamonique à une dimension de pulsation
κ
ω=
mR
où la constante de rappel κ est liée à la courbure de u(ρ) au point ρ = d. Les énergies
quantiques sont
1
ητ = (τ + )~ω, τ = 0, 1, 2...
2
La fonction de partition liée à ce mode de vibration est
1
ζ(β) =
2 sinh(β~ω/2)
cel-00339827, version 1 - 19 Nov 2008
Il s’ajoute donc à la capacité calorifique de translation du gaz de molécules diatomiques

un terme dû aux vibrations de la forme
(~ω/2kBT )2
CVv = N kB (14.3)
sinh2 (~ω/2kB T )
Il s’introduit alors une température caractéristique Tv définie par
kB Tv ≃ ~ω
Elle est d’autant plus élevée que les forces interatomiques sont plus intenses et que les
atomes sont plus légers. Elle varie selon les molécules de quelques dizaines à quelques
milliers de ◦K.
Pour T ≪ Tv les vibrations de la molécule sont “gelées” : on constate en effet que (14.3)
tend vers zéro avec T /Tv :
2
Tv
CVv ≃ N kB e−Tv /T pour T ≪ Tv
T
139
Pour T ≫ Tv on retrouve le résultat d’équipartition de l’énergie
CVv ≃ N kB , T ≫ Tr
• Degrés de liberté de rotation pour une molécule diatomique hétéropolaire :

Hétéropolaire signifie qu’elle est constituée de deux atomes différents, comme par exemple
HCl. L’Hamiltonien pour les degrés de liberté de rotation de la molécule s’écrit
1
Hr = ~2
L
2I
où I est le moment d’inertie de la molécule par rapport à son centre de masse :
I = m1 r12 + m2 r22 ≡ mRd2
Les états stationnaires de rotation sont donc les harmoniques sphériques Yj,m où j ∈ N
et m = −j, −j + 1, ..., j − 1, j avec pour valeurs propres
~2
ηj,m = j(j + 1)
cel-00339827, version 1 - 19 Nov 2008
2I
On en déduit une température caractéristique Tr définie par
~2
kB Tr =
2I
Pour la molécule H2 (homopolaire, voir sous-section suivante) Tr ≃ 85◦K. Tr est de
l’ordre de quelques dizaines de ◦K pour des molécules HA où A est un atome autre
que H. Elle est au plus de l’ordre de quelques ◦K pour toutes les autres molécules
diatomiques. On est donc presque toujours dans la situation inverse du “gel” pour la
rotation :
T ≫ Tr
La fonction de partition ζr correspondant à la rotation d’une molécule s’écrit donc :

∞ X +j ∞
X
−βηj,m
X ~2
ζr(β) = e = (2j + 1) exp −β j(j + 1)
j=0 m=−j j=0
2I
√
Si T ≫ Tr (ce qui équivaut à ~ ≪ IkB T ), la somme se ramène à une intégrale
Z ∞
~2 2I
ζr (β) ≃ dj(2j + 1) exp −β j(j + 1) = 2 kB T, T ≫ Tr
0 2I ~
La contribution correspondante à l’énergie individuelle moyenne sera donc
d
η̄r = − log ζr ≃ kB T, T ≫ Tr
dβ
140
d’où pour la capacité calorifique totale une contribution qui s’ajoute à celle du gaz parfait
monoatomique :
CVr ≃ N kB , T ≫ Tr
Ici encore le résultat d’équipartition de l’énergie est satisfait.
• degrés de liberté de rotation pour une molécule homopolaire :
(Par exemple H2 , O2 , N2 ...) On doit alors prendre en compte le postulat de symétrisation

quantique, donc le spin (total) de la molécule qui donne lieu à l’antisymétrisation de la
fonction d’onde si le spin est demi-entier (fermions), et à la symétrisation si le spin est
entier (bosons). Donc les états permis sont
symétriques dans l’échange de deux atomes pour des bosons
antisymétriques dans l’échange des deux atomes pour des fermions.
Le spin s et le moment angulaire l jouent tous deux un rôle dans cette symétrisation.
~ Lz ) satisfont
En effet on a, par échange des deux atomes que les fonctions propres de (L,
cel-00339827, version 1 - 19 Nov 2008
Ylm(π − θ, π + ϕ) = (−1)lYlm(θ, ϕ)
Par suite les états propres de rotation sont

symétriques si l est pair
antisymétriques sinon.
Si les spins sont entiers, on a affaire à des bosons, sinon à des fermions. Donc les états
quantiques pour des
bosons ont un état de spin symétrique et l pair, ou un état de spin antisymétrique et l
impair,
fermions ont un état de spin antisymétrique et l impair, ou un état de spin antisymétrique
et l pair. On aura donc :
ζ B = ♯ {états spin symétriques} × ζpair + ♯ {états spin antisym.} × ζimp.
ζ F = ♯ {états spin symétriques} × ζimp. + ♯ {états spin antisym.} × ζpair

s étant la valeur du spin individuel de l’atome, il y a 2s + 1 états symétriques de la
forme (m, m), et 21 [(2s + 1)2 − (2s + 1)] états symétriques de la forme (m, m′) +
(m′, m), m 6= m′. D’où
♯ { états spin sym.} = (s + 1)(2s + 1)
♯ { états spin antisym.} = s(2s + 1)

141
(Par ♯E on désigne le cardinal de l’ensemble E, c’est à dire le nombre d’éléments dans cet
ensemble). Par contre il n’y a que 21 [(2s + 1)2 − (2s + 1)] états antisymétriques de la
forme (m, m′) − (m′, m), m 6= m′ Notons seulement que sauf dans le cas de H2 ou
D2 ces effets quantiques sont négligeables car on est au-dessus de la température de
gel. En effet, dès que T ≪ Tr , ou a :
T
ζpair = ζimp. ≃
2Tr
et donc
(2s + 1)2 T
ζB = ζF ≃ ×
2 Tr
Pour des molécules hétéropolaires on avait

T
ζ = (2s1 + 1)(2s2 + 1)
Tr
cel-00339827, version 1 - 19 Nov 2008
Le résultat qu’on retiendra dans le cas le plus général est que si s est le spin d’un des
atomes on a
1 T
ζr(β) ≃ (2s1 + 1)(2s2 + 1)
σ Tr
où σ = 1 pour des molécules hétéropolaires et = 2 pour des molécules homopolaires.
14.3 Conclusions
Dans la Section 14-1, nous avons analysé les déviations à la loi des gaz parfaits pour des
gaz monoatomiques interagissant par des potentiels à deux corps. Cela conduit, en basse
densité, à des corrections à l’équation d’état des gaz parfaits appelés “développement du
viriel”.
Dans la Section 14-2, nous avons considéré des gaz parfaits moléculaires (donc sans
interaction) mais où on prend en compte la structure interne des molécules. On n’a
considéré que des molécules diatomiques pour faire plus simple. On obtient des corrections
à loi des gaz parfaits qui dans le cas le plus simple de molécules diatomiques dont
le cortège électronique a un moment cinétique orbital et un spin totals nuls a la forme
suivante :
Z(N, V, T ) = [z0 ζ(β)]N
142
avec s
2π~2
z0 = V λ−3 , λ=
mkB T
eβu0 2IkB T 1
ζ(β) = (2s1 + 1)(2s2 + 1)
σ ~2 sinh(~ω/2kBT )
s1 , s2 sont les spins des deux atomes
σ = 1, 2 suivant que les atomes sont identiques ou non
I est le moment cinétique de l’atome (égal à mRd2 où mR est la masse réduite des deux
noyaux, et d leur distance)
ω est le mode de vibration
−u0 est l’énergie électronique minimale de la molécule.
Plus généralement, si on est à des températures normales T telles que Tr ≪ T il y aura

nv modes de vibration τ contribuant à la fonction de partition, ceux pour lesquels ~ωτ ≤
2kB T , et on a alors que la capacité calorifique totale est égale (molécules diatomiques) à
cel-00339827, version 1 - 19 Nov 2008
nv
5N X (~ωτ /2kB T )2
CV = kB + N kB
2 τ=1
sinh2 (~ωτ /2kB T )
Chacun des termes de la somme correspond à un mode de vibration τ , donc à une

température caractéristique Tv(τ) ; la contribution du mode τ croı̂t de 0 à N kB quand T
croı̂t typiquement de Tv(τ) /10 à Tv(τ) . Il s’ensuit que les divers modes de vibration sont
excités successivement lorsque la température croı̂t, certains pouvant avoir atteint leur
température caractéristique tandis que les autres sont “gelés”. Cependant, il arrive assez
souvent que la molécule se dissocie à cause de la vibration avant que le dernier mode normal
soit excité.

Cours de MécaStat

Transféré par

Droits d'auteur :

Formats disponibles

Cours de MécaStat

Transféré par

Informations du document

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Cours de MécaStat

Transféré par

Droits d'auteur :

Formats disponibles

COURS DE MÉCANIQUE STATISTIQUE

MASTER PREMIÈRE ANNÉE

II] RAPPELS DE MÉCANIQUE CLASSIQUE

II-3-3 Les Ensembles Stationnaires

III] LES ENSEMBLES D’ÉQUILIBRE

IV] L’ÉQUIVALENCE DES ENSEMBLES D’ÉQUILIBRE

V] LES PRINCIPES DE LA THERMODYNAMIQUE

VI] APPLICATION AUX OSCILLATEURS HARMONIQUES

VI-1 La Fonction de Partition Canonique de l’Oscillateur Harmonique

VII] LES SYSTÈMES MAGNÉTIQUES

VIII] LES TRANSITIONS DE PHASE

IX] MÉCANIQUE STATISTIQUE DES SYSTÈMES QUANTIQUES

IX-3 Les Matrices-Densité Stationnaires

X] LES GAZ PARFAITS QUANTIQUES

XI] GAZ DE FERMIONS INDÉPENDANTS

XII] GAZ DE BOSONS INDÉPENDANTS

XIII] GAZ DE PHOTONS

XIV] SYSTÈMES MOLÉCULAIRES

Quelques Notions Élémentaires en

1.1 Variables Aléatoires

Plus généralement une variable aléatoire (discrète) n est spécifiée par :

On appelle P (n) la loi de probabilité de n.

Elle a une valeur comprise entre 0 et 1. Afin d’obtenir la condition de normalisation, il

comme condition de normalisation (c’est l’aire sous la courbe).

Les variables aléatoires en physique :

EXEMPLES DE LOIS DE PROBABILITÉS :

P (x) = ae−ax, x ∈ X := [0, ∞), a>0

2) Loi de Lorentz ou de Cauchy

(Gauss 1809, Laplace 1790, De Moivre 1733)

σx2 = h(x − x̄)2 i = h(∆x)2 i

la déviation standard ou écart-type : σx

Exercice 3 : Montrer que σx2 = hx2 i − hxi2

Exercice 5 : Calculer x̄ et σx pour les exemples 1) à 6) de lois de probabilités donnée

Définition 1.1.2 La transformée de Fourier de P (x)

s’appelle la fonction caractéristique de la variable aléatoire . Si X est stricte-

Exercice 6 : Montrer les propriétés suivantes :

1.1.1 Distributions de probabilités à plusieurs variables

X := [0, L]3 ⊂ R3 , P (x) ≡ P (x1 , x2 , x3 )

et la condition de normalisation s’écrit :

Exemple 2 : La variable aléatoire est la vitesse d’une particule : X = R3 , x =

Exercice 3 : Montrer que la distribution particulière dans X ≡ [0, L]3 × R3

Exercice 6 : À partir de la distribution des vitesses de Maxwell, trouver la loi marginale

Définition 1.1.3 Soit P (x, y) ≥ 0 une loi de probabilité ; alors

est la loi de probabilité de x quand la valeur de y est connue.

P (x, y) = P2 (y)P (x|y)

µ11 = ∆x2 = σx2

Variables aléatoires indépendantes

où Q, R sont des lois de probabilité arbitraires.

P (x1 , x2 , ..., xN ) = Q1 (x1 )Q2 (x2 )...QN (xN )

1.1.2 Sommes de variables aléatoires

pour tout i = 1, 2, ...N .

On se posera successivement les questions suivantes :

pour une distribution de probabilité P quelconque.

Si les variables xi, i = 1, ...N sont non-corrélées, le second terme s’annule et on a

Proposition 1.1.6 Si les variables xi i = 1, ...N sont non-corrélées, on a

Pour celle-ci, on a (pour des variables aléatoires indépendantes et identiques) :

2.1 N particules classiques de masse m

Soient N particules identiques de masse m et de positions xi(t), i = 1, 2, ...N

où r0 ≃ 3.8Angström et V0 ≃ −0.01 eV.

cinétique et de l’énergie potentielle pour le système à N -corps ainsi formé, on a :