Thermo
Thermo
LU3PY103
Cours de Thermodynamique
Éric Brunet et Lionel Forêt
avec d’importants emprunts à un autre poly rédigé par
Éric Brunet, Thierry Hocquet et Xavier Leyronas
Une théorie est d’autant plus admirable que ses fondements sont simples,
qu’elle se rapporte à des domaines variés et que son domaine d’application est étendu.
En ce sens, la thermodynamique classique m’impressionne beaucoup.
C’est la seule théorie physique de portée universelle dont je suis convaincu que,
tant qu’on s’en tient à son champ d’application, elle ne sera jamais remise en cause.
Albert Einstein, 1949, Notes for an Autobiography
eric.brunet@sorbonne-universite.fr
axel.buguin@sorbonne-universite.fr
malik.chami@sorbonne-universite.fr
matthieu.delbecq@sorbonne-universite.fr
lionel.foret@sorbonne-universite.fr
thierry.fouchet@sorbonne-universite.fr
fabien.gerbal@sorbonne-universite.fr
jules.grucker@sorbonne-universite.fr
thomas.risler@sorbonne-universite.fr
2
Table des matières
1 Rappels du L1 et du L2 9
1.1 Systèmes de la thermodynamique et variables d’état . . . . . . . . . . . . . . 9
1.1.1 Grandeurs . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
1.1.2 Définir de quoi on parle . . . . . . . . . . . . . . . . . . . . . . . . . 10
1.1.3 Quelques quantités . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
1.1.4 Modèles de gaz . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
1.2 Transformations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
1.2.1 Types de transformation . . . . . . . . . . . . . . . . . . . . . . . . . 12
1.2.2 Le premier principe . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
1.2.3 Travail des forces extérieures . . . . . . . . . . . . . . . . . . . . . . . 14
1.2.4 Échanges thermiques . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
1.2.5 Diagrammes pression-volume . . . . . . . . . . . . . . . . . . . . . . 15
1.3 Différents types d’échanges thermiques . . . . . . . . . . . . . . . . . . . . . 16
1.3.1 La conduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
1.3.2 La convection . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
1.3.3 Le rayonnement . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
1.4 Flèche du temps et entropie statistique . . . . . . . . . . . . . . . . . . . . . 17
1.5 Entropie et thermodynamique . . . . . . . . . . . . . . . . . . . . . . . . . . 17
1.6 Machines thermiques . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
1.7 Transitions de phases d’un corps pur . . . . . . . . . . . . . . . . . . . . . . 18
1.7.1 Le diagramme de phases . . . . . . . . . . . . . . . . . . . . . . . . . 18
1.7.2 Diagrammes isothermes . . . . . . . . . . . . . . . . . . . . . . . . . 19
1.7.3 Enthalpie de changement d’état . . . . . . . . . . . . . . . . . . . . . 21
1.7.4 Formule de Clapeyron . . . . . . . . . . . . . . . . . . . . . . . . . . 22
3
2.3.3 Probabilité du microétat d’une partie du système . . . . . . . . . . . 28
2.3.4 Probabilité d’une variable interne macroscopique . . . . . . . . . . . . 29
2.3.5 Comportement moyen et fluctuations d’une variable interne macrosco-
pique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30
2.3.6 Évolution irréversible d’une variable interne macroscopique . . . . . . 31
2.4 Fluctuations et irréversibilité . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
2.5 Entropie . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34
2.5.1 Définition et exemple . . . . . . . . . . . . . . . . . . . . . . . . . . . 34
2.5.2 Propriété de l’entropie . . . . . . . . . . . . . . . . . . . . . . . . . . 35
2.5.3 Équilibre thermique et température . . . . . . . . . . . . . . . . . . . 36
2.5.4 Équilibre mécanique et pression . . . . . . . . . . . . . . . . . . . . . 38
2.5.5 Potentiel chimique . . . . . . . . . . . . . . . . . . . . . . . . . . . . 39
2.5.6 Résumé des conditions d’équilibre . . . . . . . . . . . . . . . . . . . . 40
2.6 Identité thermodynamique . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40
2.6.1 Cas d’un fluide simple . . . . . . . . . . . . . . . . . . . . . . . . . . 40
2.6.2 Identité thermodynamique généralisée . . . . . . . . . . . . . . . . . . 42
2.6.3 Entropie et chaleur . . . . . . . . . . . . . . . . . . . . . . . . . . . . 42
2.7 Utilisation pratique de la statistique microcanonique . . . . . . . . . . . . . 43
2.7.1 Le système à deux niveaux . . . . . . . . . . . . . . . . . . . . . . . . 44
3 Fonctions thermodynamiques 45
3.1 Réservoirs : thermostat, pressostat . . . . . . . . . . . . . . . . . . . . . . . 45
3.2 Équilibre d’un système en contact avec un réservoir . . . . . . . . . . . . . . 46
3.2.1 Variation d’entropie du réservoir . . . . . . . . . . . . . . . . . . . . . 46
3.2.2 Inégalité de Clausius . . . . . . . . . . . . . . . . . . . . . . . . . . . 48
3.2.3 Potentiels thermodynamiques . . . . . . . . . . . . . . . . . . . . . . 48
3.2.4 Exemple d’utilisation du potentiel thermodynamique . . . . . . . . . 50
3.2.5 Énergie libre et enthalpie libre . . . . . . . . . . . . . . . . . . . . . . 51
3.3 Les fonctions thermodynamiques . . . . . . . . . . . . . . . . . . . . . . . . 52
3.3.1 Présentation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 52
3.3.2 Variables naturelles et dérivées premières des fonctions thermodyna-
miques . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53
3.3.3 Relations de Maxwell et dérivées secondes des fonctions thermodyna-
miques . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 54
3.3.4 Extensivité et additivité . . . . . . . . . . . . . . . . . . . . . . . . . 54
3.3.5 Généralisations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 55
3.3.6 Exemple : diagramme de phases et formule de Clapeyron . . . . . . . 56
3.4 Coefficients calorimétriques et thermoélastiques . . . . . . . . . . . . . . . . 57
3.4.1 Définitions des coefficients calorimétriques . . . . . . . . . . . . . . . 57
3.4.2 Relations de Clapeyron . . . . . . . . . . . . . . . . . . . . . . . . . . 58
3.4.3 Définitions des coefficients thermoélastiques . . . . . . . . . . . . . . 59
3.4.4 Relations entre les coefficients . . . . . . . . . . . . . . . . . . . . . . 60
4 Mélanges 63
4.1 Mélange de gaz parfaits, mélange idéal . . . . . . . . . . . . . . . . . . . . . 64
4.1.1 Expression du potentiel chimique . . . . . . . . . . . . . . . . . . . . 64
4.1.2 Application : taux d’humidité, évaporation, ébullition . . . . . . . . . 66
4
4.2 Mélanges dilués . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 68
4.2.1 Expression du potentiel chimique . . . . . . . . . . . . . . . . . . . . 68
4.2.2 Application : la pression osmotique . . . . . . . . . . . . . . . . . . . 70
4.2.3 Application : Déplacement de la température de transition . . . . . . 71
4.3 Mélange binaire non dilué sous plusieurs phases . . . . . . . . . . . . . . . . 72
4.3.1 Variance du système . . . . . . . . . . . . . . . . . . . . . . . . . . . 72
4.3.2 Diagramme de phases d’un mélange binaire . . . . . . . . . . . . . . 73
4.3.3 Application : la distillation . . . . . . . . . . . . . . . . . . . . . . . . 75
6 La fonction de partition 87
6.1 Distribution de l’énergie d’un système . . . . . . . . . . . . . . . . . . . . . . 88
6.1.1 Énergie moyenne . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 88
6.1.2 Variance et écart-type . . . . . . . . . . . . . . . . . . . . . . . . . . 89
6.1.3 Fluctuations de l’énergie . . . . . . . . . . . . . . . . . . . . . . . . . 89
6.1.4 Distribution de l’énergie . . . . . . . . . . . . . . . . . . . . . . . . . 90
6.2 Liens avec les fonctions thermodynamiques . . . . . . . . . . . . . . . . . . . 92
6.2.1 L’énergie libre . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 92
6.2.2 L’entropie, la pression, le potentiel chimique . . . . . . . . . . . . . . 93
6.3 Factorisation de Z en l’absence d’interactions . . . . . . . . . . . . . . . . . 94
6.3.1 Un premier exemple dans le cas discernable . . . . . . . . . . . . . . 94
6.3.2 Le cas général pour des particules discernables . . . . . . . . . . . . . 95
6.3.3 Le cas des particules indiscernables . . . . . . . . . . . . . . . . . . . 96
6.3.4 La factorisation de Z . . . . . . . . . . . . . . . . . . . . . . . . . . . 97
6.3.5 Généralisation de la factorisation . . . . . . . . . . . . . . . . . . . . 98
6.4 Exemple : assemblée de N oscillateurs harmoniques . . . . . . . . . . . . . . 98
6.5 La limite classique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 99
6.6 Résumé . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 100
5
7.3.1 Contribution des termes de vibration . . . . . . . . . . . . . . . . . . 106
7.3.2 Contribution des termes de rotation . . . . . . . . . . . . . . . . . . . 107
7.3.3 Les autres termes, et un bilan . . . . . . . . . . . . . . . . . . . . . . 109
7.4 La capacité thermique du gaz parfait . . . . . . . . . . . . . . . . . . . . . . 110
7.4.1 Le gaz monoatomique . . . . . . . . . . . . . . . . . . . . . . . . . . 110
7.4.2 Le gaz diatomique . . . . . . . . . . . . . . . . . . . . . . . . . . . . 110
7.5 Retour sur la limite classique (section optionnelle) . . . . . . . . . . . . . . . 111
6
Avant propos et bibliographie
La thermodynamique est une branche de la physique qui étudie les systèmes macrosco-
piques (c’est-à-dire constitués d’un grand nombre de particules) dans le cadre d’une ap-
proche énergétique. C’est une théorie qui s’applique à de nombreux systèmes et qui permet
d’établir des relations très générales entre les coefficients décrivant les états de la matière.
La physique statistique donne une interprétation microscopique aux quantités étu-
diées en thermodynamique. Dans les cas simples, les postulats de la physique statistique
permettent de retrouver, comprendre et interpréter les lois de la thermodynamique.
Ces notes de cours s’adressent à des étudiants ayant déjà de solides notions de ther-
modynamique, et plus spécifiquement qui ont suivi en L2 le cours de thermodynamique de
Sorbonne Université. Le premier chapitre de ce document rappelle rapidement toutes les
notions supposées connues. Le chapitre 2 expose les bases de la physique statistique pour un
système isolé, c’est-à-dire la physique statistique microcanonique, et donne une interpréta-
tion microscopique à l’entropie. Les propriétés de l’entropie sont revues en détail. Dans le
chapitre 3, on introduit les potentiels thermodynamiques, et on montre comment les utiliser
pour établir des relations entre les différentes propriétés des systèmes. Les chapitres 5 et 6
traitent de la physique statistique canonique, en contact avec un thermostat. Les chapitres
7 et 8 sont des applications de la physique statistique aux gaz diatomiques et aux gaz de
photons. Enfin, le chapitre 4 clôt ce cours par un retour à la thermodynamique, où l’on étudie
les propriétés des mélanges sous une ou plusieurs phases. Quelques points mathématiques
sont expliqué dans l’appendice A.
Bibliographie sommaire
• Lhuillier et Rous Introduction à la thermodynamique, Dunod (1996)
Livre concis, idéal pour revoir les notions de base de la thermodynamique
• Jancovici, Thermodynamique et physique statistique, Nathan (1999)
Un autre livre concis, qui aborde également la physique statistique
• Simon, Énergie et entropie, Armand Colin (1979)
Une référence pour approfondir les concepts et les bases de la description microscopique
• Diu, Guthmann, Lederer et Roulet, Thermodynamique, Hermann (2007)
Ouvrage de référence, très épais, très complet. Contient de nombreux compléments et suppléments.
• Landau et Lifshitz Physique Statistique, Ellipses (1994)
Plus difficile (en deuxième lecture) : Un classique, indémodable, extrêmement concis !
• Callen Thermodynamics and an introduction to thermostatistics, Wiley (1985)
Livre en anglais, dont la lecture est hautement recommandée.
• Reif Cours de Physique de Berkeley. Tome 5 : Physique Statistique, Dunod (2000)
Un autre classique, détaillé, avec des illustrations, dont vous apprécierez la lecture !
• Reif Fundamentals of statistical and thermal physics, Mc Graw-Hill (1964)
Livre en anglais du même auteur que le précédent, plus avancé et rapide dans sa présentation.
7
8
Chapitre 1
Rappels du L1 et du L2
Ce chapitre ne sera pas traité en amphi, puisqu’il n’y a aucune nouvelle information. Il
s’agit juste d’un rappel synthétique de toutes les notions déjà vues en L1 et en L2 qui ne
seront pas revues dans le cours de L3.
Les notions déjà vues en L2 et revues en L3 sont la notion d’entropie, en physique
statistique et en thermodynamique, chapitre 2.
9
• L’unité de pression est le pascal (Pa). On utilise aussi 1 bar = 105 Pa et 1 atm =
1,013 bar.
• Un angström (Å) vaut 10−10 m. Un atome commun ou une molécule simple mesure
quelques angströms.
• Dans un solide ou un liquide, les distances typiques entre atomes ou molécules voisines
sont de quelques angströms. Dans un gaz aux conditions de température et de pression
usuelle, les distances typiques entre atomes ou molécules voisines sont de quelques
nanomètres. Ces distances typiques peuvent être évaluées en prenant la racine cubique
du volume par particule.
10
1.1.3 Quelques quantités
• Bien qu’un système à l’équilibre thermodynamique soit macroscopiquement au re-
pos, les particules ont un mouvement microscopique, et donc une énergie cinétique.
L’énergie interne U du système est la somme de cette énergie cinétique microsco-
pique et de l’énergie d’interaction entre les particules.
• La température T d’un système est une fonction d’état telle que a) à l’équilibre
thermodynamique, deux systèmes séparés par une paroi diatherme sont à la même
température (c’est le « principe zéro » de la thermodynamique) et b) quand deux
systèmes de températures différentes sont en contact thermique, l’énergie a tendance
à s’écouler du corps le plus chaud vers le corps le plus froid.
• La température peut être définie par T = pV /(nR) pour un gaz très dilué (limite du
gaz parfait). En calculant la pression cinétique du gaz parfait, c’est-à-dire la pression
causée par tous les chocs des particules sur les parois, on montre que, pour un gaz par-
fait, l’énergie cinétique microscopique de translation 3 Ec est donnée par Ec = 23 nRT .
• En fait, pour n’importe quel système qui peut être décrit par la mécanique classique
(pas seulement le gaz parfait), l’énergie cinétique de translation est Ec = 23 nRT . Cette
propriété, dont on reparlera dans la section 5.2.5, est une définition possible de la tem-
pérature (comme une mesure de l’agitation thermique), mais elle ne fonctionne plus
pour des systèmes fortement quantiques comme des électrons dans un métal.
• La bonne définition de la température, celle qui marche dans tous les cas, est celle
donnée à partir de l’entropie par T1 = ∂U∂S
n,V
, voir le chapitre 2.
• La vitesse quadratique moyenne vq est la vitesse telle que le système aurait l’éner-
gie cinétique de translation Ec = 32 nRT si on supposait que les particules avaient
toutes cette vitesse vq . On trouve vq = 3kB T /m = 3RT /M avec m la masse d’une
p p
11
• L’entropie S est une fonction d’état. Pour un système isolé, l’entropie augmente si la
transformation est irréversible et reste constante si la transformation est réversible.
• D’un point de vue microscopique, l’entropie d’un système isolé est donnée par S =
kB ln Ω, où Ω est le nombre de micro-états compatibles avec l’état thermodynamique
(macroscopique) du système.
1.2 Transformations
1.2.1 Types de transformation
• Une transformation d’un système thermodynamique est le passage (dû à des mo-
difications de paramètres extérieurs) d’un état d’équilibre interne vers un autre état
d’équilibre interne. La variation d’une variable d’état X lors d’une transformation est
notée ∆X = Xfinal − Xinitial . Pour une variation infinitésimale, cette variation sera
plutôt notée dX.
• Une transformation est dite quasistatique si le système est à tout moment dans un
état d’équilibre interne. (Par exemple, le café à 50°C qui refroidit subit une transfor-
mation quasistatique.) On peut décomposer une transformation R quasistatique en une
succession de transformations infinitésimales et écrire ∆X = dX.
• Souvent, on considère des transformations quasistatiques où l’on suppose en plus que
toutes les parties du systèmes sont à tout moment à l’équilibre mécanique. En parti-
culier, quand le système est un fluide fermé par un piston, on doit avoir p = pext à tout
moment pour une transformation de ce type.
• Une transformation réversible est une transformation telle que le système est à tout
instant à l’équilibre, de manière interne et avec l’extérieur. Par exemple, la fusion de la
12
glace dans un environnement à 0°C est une transformation réversible. Toute transfor-
mation réversible est quasistatique, la réciproque n’est pas vraie. Si une transformation
est réversible, alors la transformation inverse (en échangeant les états initial et final)
est également une transformation réversible possible. Ainsi, la solidification de l’eau
liquide dans un environnement à 0°C est une transformation réversible possible.
• Une transformation est
• monobare si la pression extérieure est constante.
• isobare si la pression du système est constante et égale à la pression extérieure.
• monotherme si la température extérieure est constante.
• isotherme si la température du système est constante et égale à la température
extérieure.
• isochore si le volume du système est constant.
• adiabatique si le système n’échange pas de chaleur avec l’extérieur.
• Une transformation telle que le volume diminue est une compression. Une transfor-
mation telle que le volume augmente est une détente.
• Lors d’une transformation quasistatique, adiabatique, avec p = pext , pour un gaz parfait,
on a les lois de Laplace
pV γ = cste, p1−γ T γ = cste, T V γ−1 = cste, (1.3)
où γ = cp /cV est l’indice adiabatique.
• Un thermostat est un dispositif dont la température est constante, quelle que soit la
chaleur échangée. Un système qui n’est en contact thermique qu’avec un thermostat
subit des transformations monothermes.
• Un pressostat est un dispositif dont la pression est constante. Un système dont toutes
les parois mobiles donnent sur un pressostat donné subit des transformations mono-
bares.
• On considère souvent que l’atmosphère joue le rôle de thermostat et de pressostat.
13
• Lors d’une transformation, la variation d’énergie totale du système est égale à l’énergie
reçue ; c’est le premier principe :
∆Etotale = W + Q. (1.4)
Dans la plupart des cas, on considère des systèmes immobiles : Ec-macro et Ep-macro sont
nuls ou constants, et le premier principe s’écrit
∆U = W + Q, (1.6)
δW = −pext dV.
14
• On introduit l’enthalpie H d’un système par
H = U + pV. (1.9)
15
se lit sur le diagramme en considérant l’aire comprise entre la courbe de la trans-
formation et l’axe des abscisses (l’aire « sous la courbe ») : pour une compression,
W = (aire sous la courbe), et pour une détente, W = −(aire sous la courbe).
• Un cycle est une transformation dont l’état final est égal à l’état initial. Sur un dia-
gramme pression-volume, un cycle est représenté par une courbe fermée.
• Toujours sous l’hypothèse p = pext , le travail reçu par le système lors d’un cycle est :
W = (aire dans le cycle) si la représentation de la transformation tourne dans le sens
anti-horaire (on dit alors que le cycle est récepteur), et W = −(aire dans le cycle) si
la représentation de la transformation tourne dans le sens horaire (on dit alors que le
cycle est moteur).
• Au cours d’un cycle, puisque l’état final du système est égal à son état initial, aucune
variable d’état ne varie : ∆U = 0, ∆S = 0, etc.
λS
ϕ= (T1 − T2 ). (1.11)
L
1.3.2 La convection
• Dans un fluide, l’énergie est transportée par le mouvement des masses de fluide plus
ou moins chaudes. C’est la convection. On distingue la convection forcée, quand le
fluide est mis en mouvement de manière extérieure, et la convection libre, quand le
mouvement se met en place spontanément.
• Quand un solide est présent dans un fluide à une température différente, on modélise
le courant de chaleur ϕ à travers une surface S du corps 1 vers le corps 2 comme pour
la conduction : ϕ = hS(T1 − T2 ) ; c’est la loi de Newton.
16
1.3.3 Le rayonnement
• Un corps peut absorber un rayonnement électromagnétique, et donc gagner de
l’énergie. (Il peut aussi réfléchir le rayonnement et ne pas gagner d’énergie. En général,
il absorbe une partie du rayonnement et réfléchit l’autre, dans des proportions qui
dépendent du corps et de la longueur d’onde.)
• Spontanément, tout corps émet du rayonnement électromagnétique. Pour un
corps noir (parfaitement absorbant), la distribution en fréquence et en intensité de
ce rayonnement ne dépend que de la température du corps : il est infrarouge aux
températures ambiantes et devient visible vers 1000°C. C’est le rayonnement du
corps noir.
17
• Le second principe de la thermodynamique implique que les machines thermiques les
plus efficaces sont les machines réversibles, et que toutes les machines réversibles ont
la même efficacité.
• La machine thermique réversible la plus simple est le cycle de Carnot pour un gaz
parfait : ce cycle est composé de deux isothermes et deux adiabatiques et on trouve
par un calcul direct que, pour ce cycle réversible, Qchaud /Tchaud + Qfroid /Tfroid = 0.
• À partir de l’égalité précédente et du premier principe W + Qchaud + Qfroid = 0, on peut
obtenir les efficacités des machines réversibles. Par exemple, un moteur thermique
réversible a une efficacité η = 1 − Tfroid /Tchaud . C’est le rendement de Carnot.
• Pour une machine thermique quelconque, le second principe implique
Qchaud Qfroid
+ ≤ 0. (1.12)
Tchaud Tfroid
C’est l’inégalité de Clausius. À partir de cette inégalité, on peut montrer par
exemple que le rendement d’un moteur est toujours inférieur au rendement de Car-
not.
18
la pression pS (T ). Réciproquement, à une pression p donnée, le changement de phase
ne peut avoir lieu qu’à la température T telle que p = pS (T ).
• À titre d’exemple, les diagrammes de phase de l’eau et du CO2 :
106
105
104 Point critique:
3 p = 220 atm
10
T = 374◦ C
Fluide supercritique
100
p Eau liquide
(atm)
10 Diagramme de
1 Glace phases de l’eau
0,1
Vapeur
0,01 Point triple:
10−3 p = 0,006 atm
T = 0,01◦ C
10−4
−100 −50 0 50 100 150 200 250 300 350 400 450
T (◦ C)
104
Solide
3
10
Point critique:
Liquide T = 31◦ C
Fluide supercritique
100 p = 73 atm
p
(atm) Point triple: Diagramme de
10 T = −57◦ C phases du CO2
p = 5 atm
Gaz
0,1
−100 −80 −60 −40 −20 0 20 40 60
T (◦ C)
Notez l’échelle logarithmique pour la pression, et la présence de plusieurs phases solides
pour l’eau. Notez aussi que le CO2 ne peut pas être liquide en-dessous de 5 atm : à
pression ambiante, la neige carbonique (CO2 solide) se sublime sans fondre. Si on
comprime du CO2 gazeux à 20°C, les premières gouttes apparaissent vers 60 atm. À
40°C, on ne voit jamais apparaître de goutte car le fluide est supercritique.
19
p
Liquide seul
C
ML M MG
pS (T2 ) T1
Liquide + Gaz
Gaz seu TC
l
T2
v = V /n
vL (T2 ) v vG (T2 )
Figure 1.1 – Un diagramme isotherme pour un gaz de van der Waals, avec trois courbes
correspondant respectivement à la température critique TC , une température T1 > TC et une
température T2 < TC . Chaque courbe s’appelle une isotherme.
20
l’isotherme critique (à température Tc , en rouge sur la figure), les différences entre
liquide et gaz disparaissent et le fluide est dit supercritique.
p
Fluide supercritique
C
Liquide
Gaz
TC
courbe courbe
d’ébullition de rosée
Liquide + Gaz
v = V /n
Figure 1.2 – Dans le diagramme (v, p), les régions correspondant au différents états du
système. On a aussi aussi tracé les isothermes de la figure 1.1.
21
1.7.4 Formule de Clapeyron
• L’enthalpie molaire lα β à la température T est reliée à la pression de coexistence pα β (T )
entre ces deux phases par la formule de Clapeyron.
dpα β
lα β = T (vβ − vα ) , (1.13)
dT
où vα et vβ sont les volumes molaires des phases α et β.
dp
• On voit donc que si lα β > 0, alors vβ − vα et dTα β ont le même signe. Pour la plupart
des corps purs, on a vsolide < vliquide < vgaz et, comme lfusion , lvaporisation et lsublimation sont
positives, on en déduit que les pentes des lignes de coexistence du diagramme de phase
sont positive, voir par exemple le diagramme du CO2 dans la section 1.7.1. Cependant,
pour l’eau, on a vglace > veau liquide ce qui implique que la ligne de coexistence solide-
liquide du diagramme de phase de l’eau a une pente négative.
22
Chapitre 2
23
Figure 2.1 – Trois microétats d’un système de 12 particules classiques. Les flèches repré-
sentent les vitesses des particules
24
2.1.3 Statistiques microcanonique et canonique
Comme nous allons le voir, la distribution de probabilité Pσ suit une loi très générale,
indépendante de la nature du système et de ces constituants. Cette loi ne dépend en fait que
de la situation dans laquelle se trouve le système : isolé (son énergie est fixée) ou alors en
contact avec un réservoir d’énergie, aussi appelé thermostat (sa température est fixée). On
parle de statistique microcanonique pour un système isolé et de statistique canonique pour
un système en contact avec un thermostat. Si la statistique canonique est la plus utilisée
en pratique, ses principes se déduisent de la statistique microcanonique qui constitue le
fondement de la physique statistique. Nous allons donc commencer par étudier la statistique
microcanonique qui permet en outre de mieux comprendre l’origine de la seconde loi de la
thermodynamique et la nature de l’entropie.
On pose
Ω= Le nombre de microétats accessibles .
au système étant donné son macroétat
alors
1
Pσ =
Ω
Remarques :
• « Compatibles avec son macroétat » signifie en particulier que, parmi tous les microé-
tats possibles du système, il ne faut considérer que ceux pour lesquels l’énergie du
système vaut U .
• De façon générale, dans la description quantique d’un système, le nombre Ω de microé-
tats correspond au nombre d’états propres de l’hamiltonien ayant l’énergie U . C’est un
nombre fini.
• Pour les systèmes décrits en bonne approximation par la mécanique classique, le mi-
croétat est caractérisé par des grandeurs continues (vitesses, coordonnées de positions).
Déterminer le nombre de microétats, Ω, est impossible, celui-ci est infini. On s’en sort
en supposant que positions et vitesses ne sont connues qu’avec une précision donnée,
25
petite mais non nulle ; il y a alors un nombre fini de positions et de vitesses différentes,
et Ω est fini. Bien sûr, l’expression de Ω dépend alors des précisions choisies. Il s’avert
qu’en pratique, ce choix n’a guère d’importance.
26
Pour chacune de ces neuf possibilités, il y a huit manières de placer la particule 2 . Par
exemple, en partant du troisième choix ci-dessus, on a :
2 1 2 1 1 1 1 1 1 1
2 2 2
2 2 2
L’effet de numéroter les boules a donc été de compter chaque microétat six fois. On
corrige ce surcomptage en divisant notre résultat par 6, où 6 = 3! = 3 × 2 × 1 est le nombre
de permutations de { 1 , 2 , 3 }, c’est-à-dire le nombre d’ordres dans lequel on peut lister ces
trois éléments. On trouve finalement :
9×8×7
Ω= = 84.
6
En généralisant maintenant ce résultat à un N et un λV quelconque, on trouve
N termes
z }| {
λV × (λV − 1) × (λV − 2) × · · · × (λV − N + 1)
Ω(N, V ) = . (2.2)
N!
Le numérateur est le nombre de manières de disposer N particules numérotées dans λV
cases. Les particules étant, en fait, indiscernables, on corrige le surcomptage en divisant par
N ! = N × (N − 1) × · · · × 2 × 1, qui est le nombre de permutations à N éléments.
En multipliant par (λV − N )! le numérateur et le dénominateur, on forme au numérateur
(λV )! et on obtient
(λV )! λV
Ω(N, V ) = = . (2.3)
N !(λV − N )! N
Le notation λV (coefficient binomial ) se lit « N parmi λV » et représente le nombre de
N
manières de choisir N objets parmi λV .
Tous les microétats du système sont équiprobables. La probabilité de l’un d’entre eux est
donc 1/Ω(N, V ).
27
On note Ncase , le nombre de particules dans une case donnée. Ce nombre peut valoir 1 ou
0 (case occupée ou non) selon le microétat du système. Si on regarde la figure représentant
l’ensemble des microétats de 3 particules sur 9 cases, on voit que pour une case donnée
(la première case en haut, à gauche par exemple), il y a 28 microétats où Ncase = 1, et 56
microétats où Ncase = 0. On applique la formule (2.1) à l’observable A = Ncase ; Pσ = 1/84
indépendamment de σ, et on somme sur les 84 microétats. Pour 28 d’entre-eux, Aσ = (Ncase )σ
vaut 1, et pour les 56 autres Aσ vaut 0. On obtient donc ⟨Ncase ⟩ = 28×(1× 84
1 1
)+56×(0× 84 )=
28
84
= 1
3
.
Si on généralise, il y a λV −1
microétats où la case est occupée (nombre de façons de
N −1
placer les N − 1 autre particules sur les λV − 1 cases restantes) et λVN−1 microétats où la
case est vide (nombre de façons de placer les N particules sur les λV − 1 cases restantes).
On a alors,
" # " #
λV − 1 1 λV − 1 1 N
⟨Ncase ⟩ = × 1 × λV + × 0 × λV = . (2.4)
N −1 N
N N
λV
28
α β α β α β α β
l’état σα donné ; il s’agit en effet du nombre de manières de ranger les particules restantes
(ici, il y en a quatre) dans les quinze cases de la partie β. Avec le postulat d’équiprobabilité
on en déduit que la probabilité d’observer l’état σα donné est 1365/3 268 760 ≈ 0,04 %.
De manière générale, on a
Ω(N − Nα , Vβ )
Proba(σα ) = , (2.6)
Ω(N, V )
Ω(Nα , Vα )Ω(N − Nα , Vβ )
Proba(Nα ) = . (2.7)
Ω(N, V )
Avec l’expression (2.3) de Ω(N, V ), on peut calculer les probabilités (2.7) pour les valeurs
29
de l’exemple (λVα = 10, λVβ = 15, N = 10). On trouve
30
On voit clairement que plus ⟨Nα ⟩ est grand, plus la distribution de Nα est piquée autour
de sa valeur moyenne. Pour les plus grandes valeurs, on pourra négliger les fluctuations et
écrire
Nα ≈ ⟨Nα ⟩ si ⟨Nα ⟩ est grand. (2.8)
Application numérique : pour l’atmosphère standard, on a ρ = N V
= kBpT ≈ 2,4 1025 m−3 . À
l’intérieur d’un récipient de 1 ℓ ouvert sur l’atmosphère, on trouve ⟨Nα ⟩ = ρVα ≈ 2,4 1022 .
Les fluctuations de Nα sont négligeables et peuvent être ignorées.
On peut être plus précis. Un calcul un peu long permet de trouver l’écart type ∆ du
nombre de particules Nα . On trouve, toujours dans la limite N → ∞, V → ∞, ρ = N V
:
q q
ρ
1 − λρ ⟨Nα ⟩. (2.9)
∆ = ρ 1 − λ Vα =
31
partie β (à droite). On retire cette cloison (on dit qu’on a « levé une contrainte ») ; les
particules ont maintenant accès à tout le volume. On attend d’être à nouveau à l’équilibre
thermodynamique.
α β α β
Figure 2.3 – Deux micro-états du système : à gauche, avant qu’on ne lève la contrainte et,
à droite, après qu’on a levé la contrainte.
Posons Nα et Nβ les nombres respectifs de particules dans les compartiments α et β.
Dans le macroétat initial, on a par construction Nα = N et Nβ = 0. Dans le macroétat final,
les valeurs de Nα et Nβ sont aléatoires (elles dépendent du microétat). On a cependant vu
que les valeurs moyennes sont données par
N
⟨Nα ⟩ = ρVα , ⟨Nβ ⟩ = ρVβ , avec ρ = . (2.12)
V
Ω(N, Vα )
pretour = Proba(Nα = N ) = . (2.13)
Ω(N, V )
(Nombre de cas favorables sur nombre total de cas . C’est un cas particulier de (2.7).) On
peut vérifier (exercice) que l’on a pretour ≤ (Vα /V )N . Si on prend par exemple Vα = 12 V et
N = 1022 , on trouve donc pretour ≤ 2−(10 ) , qui est effectivement extrêmement proche de
22
0, suffisamment proche pour qu’on puisse être certain que l’événement ne se produira pas
pendant la durée de vie de l’univers.
32
On arrive à la conclusion que lorsqu’on a levé la contrainte (en permettant aux
particules de circuler entre les deux parties), le système a évolué de manière irréversible
vers un nouvel état d’équilibre.
Si on prend un système à l’équilibre avec N = 1020 particules, on voit donc que la valeur
aléatoire de X coïncide avec la valeur moyenne ⟨X⟩ sur environ 10 chiffres significatifs. On
peut donc ignorer les fluctuations et écrire X = ⟨X⟩. Cette valeur d’équilibre ⟨X⟩ est celle
qui maximise la fonction (2.14).
Comme dans la section 2.3.6, on considère maintenant la transformation suivante :
1. Initialement, la valeur de la variable X est fixée par une contrainte à X = Xinitial .
2. Cette contrainte est levée de sorte que X devient une variable libre.
3. On attend que le système retourne à l’équilibre thermodynamique.
Dans l’état initial, X = Xinitial et le nombre de microétats est Ωcontraint (U, V, N, Xinitial ).
Dans l’état final, X est égal à la valeur d’équilibre thermodynamique sans contrainte ⟨X⟩ (à
33
des fluctuations indétectables près), et le nombre de microétats est Ω(U, V, N ). Si Xinitial ̸=
⟨X⟩, la probabilité d’observer le système revenir à un certain instant dans son état initial
avec X = Xinitial est quasi-nul car on a alors Ωcontraint (U, V, N, Xinitial ) ≪ Ω(U, V, N ).
Encadré 2.3 : Évolution irréversible
Dans la limite macroscopique, une fois relâchée la contrainte qui fixait sa valeur, une va-
riable interne X dans un système isolé évolue irréversiblement vers sa valeur d’équilibre
⟨X⟩ qui est aussi sa valeur la plus probable.
Une évolution inverse a une probabilité quasi-nulle de se produire .
2.5 Entropie
2.5.1 Définition et exemple
Encadré 2.4 : Définition de l’entropie
S = kB ln Ω
Pour illustrer cette définition, reprenons l’exemple du gaz sur réseau. Connaissant son
nombre de microétats (2.3), on peut calculer l’entropie du gaz en application la définition
de l’encadré 2.4 :
(λV )! h i
S = kB ln Ω = kB ln = kB ln[(λV )!] − ln[N !] − ln[(λV − N )!] . (2.15)
N !(λV − N )!
Pour un grand système, on a λV ≫ 1, N ≫ 1 et λV − N ≫ 1. Les factorielles sont alors
données en très bonne approximation par la formule de Stirling.
Encadré 2.5 : La formule de Stirling
ln(n!) ∼ n ln n − n.
Les termes sans logarithme se simplifient, et l’on obtient finalement l’entropie du gaz fonction
des variables qui caractérisent le macroétat du système :
h i
S(N, V ) = kB λV ln(λV ) − N ln N − (λV − N ) ln(λV − N ) . (2.16)
34
2.5.2 Propriété de l’entropie
On peut déduire de sa définition (encadré 2.4) des propriétés importantes de l’entropie.
Encadré 2.6 : Propriétés de l’entropie
S = S(U, V, n).
S(U, V, n) = n × s(u, v)
qui est l’entropie qu’aurait le système si X était fixée par une contrainte.
Un système isolé évolue donc vers l’état d’équilibre compatible avec les contraintes
pour lequel S est maximum. En particulier, lorsqu’on lève une contrainte dans
un système isolé, le système cherche un nouvel état d’équilibre et son entropie ne
peut qu’augmenter.
L’entropie est croissante lors d’une transformation adiabatique
Pour un système isolé ou n’échangeant que du travail avec l’extérieur,
∆S ≥ 0.
35
Remarques :
• Les variables de l’entropie (U , V , n, . . .) sont des grandeurs physiques conservées dans
un système isolé et qui sont définis même en dehors de l’équilibre thermodynamique
(contrairement à S). Elles sont de nature géométrique (comme V ), mécanique (comme
U ), chimique, électromagnétique, . . .. On considèrera toujours que ces variables sont
des grandeurs continues, même si ce n’est qu’une approximation. La fonction S sera
supposée continue et plusieurs fois dérivable en toutes ses variables.
• L’additivité se comprend facilement pour un système formé de deux sous-systèmes α
et β isolés l’un de l’autre : le nombre total de microétats du système entier est alors
Ω = Ωα × Ωβ où Ωα et Ωβ sont les nombres de microétats des sous-systèmes, et on
obtient pour l’entropie du système entier :
La propriété d’additivité reste vraie même si les deux sous-systèmes ne sont pas isolés
l’un de l’autre, à condition qu’ils soient macroscopiques, qu’il ne s’agisse pas d’un
mélange de deux fluides, et que les interactions soient de courte portée. (La portée
des forces intermoléculaires de deux systèmes non chargés est de l’ordre de quelques
dizaines de nanomètres.) En effet, les interactions mutuelles entre deux sous-systèmes
macroscopiques sont très faibles sauf dans une région de l’espace de faible volume (égal
à la surface de l’interface fois la portée de l’interaction), négligeable devant le volume
total du système. Tout se passe alors comme si les deux sous-systèmes étaient isolés
l’un de l’autre.
• L’extensivité découle de l’additivité : si on juxtapose λ copies d’un même système de
macroétat (U, V, n), on obtient (par additivité) un système d’entropie λS(U, V, n) ayant
une énergie λU , un volume λV , et une quantité λn de matière.
• Le principe potentiel de l’entropie (qui dit qu’elle est maximale à l’équilibre) est une
conséquence de l’encadré 2.3 qui dit qu’une variable X évolue vers sa valeur la plus
probable (les fluctuations étant négligeables), et de l’équation (2.14) qui implique que
la valeur la plus probable est celle qui maximise Ωcontraint (U, V, n, X).
• Au cours d’une transformation, un système isolé ne reçoit ni travail, ni chaleur, ni
matière. La transformation est induite par relâchement d’une contrainte. D’après le
principe précédent, l’entropie ne peut donc que croître au cours d’une telle transfor-
mation : ∆S ≥ 0.
Considérons maintenant un système fermé qui ne peut échanger de l’énergie que sous
forme de travail avec un système purement mécanique (sans entropie). L’ensemble
formé par les deux systèmes est isolé et, en appliquant le principe précédent, on obtient
∆S +∆Ssystème mécanique ≥ 0. Un système mécanique n’ayant pas d’entropie, on retrouve
∆S ≥ 0.
On peut généraliser à toute transformation adiabatique, au cours de laquelle le système
ne reçoit que du travail (mais pas de chaleur ni de matière), car n’importe quel travail
pourrait en principe être fourni à l’identique par un système purement mécanique.
36
niveau microscopique dans la zone de contact, voir figure 2.4. La répartition de l’énergie U
entre les deux sous-systèmes n’est pas fixée a priori ; Uα l’énergie du sous-système α est une
variable libre de fluctuer. (L’énergie de β fluctue aussi mais sa valeur Uβ = U − Uα se déduit
directement de celle de Uα .)
α β
Uα énergie Uβ = U − Uα
Figure 2.4 – Un système isolé constitué de deux sous-systèmes pouvant échanger de l’éner-
gie.
Afin d’alléger les notations, les dépendances en V et N des entropies, non-pertinentes
ici, seront omises. D’après l’encadré 2.6, la valeur d’équilibre de Uα est celle qui maximise
l’entropie S(U, Uα ) qu’aurait le système si Uα était fixée, c’est-à-dire si les deux sous-systèmes
ne pouvaient pas échanger d’énergie. Par la propriété d’additivité, l’entropie à Uα fixée est
où Sα (Uα ) et Sβ (Uβ ) sont les entropies qu’auraient les parties α et β si elles étaient isolées
avec l’énergie Uα ou Uβ . La valeur de Uα à l’équilibre maximise cette fonction. On se rappelle
que U est fixée, et on obtient donc,
∂Uβ
Or ∂Uα U
= −1 car Uβ = U − Uα . Finalement, la valeur de Uα à l’équilibre est telle que
∂Sα ∂Sβ
= . (2.20)
∂Uα ∂Uβ
Tα = Tβ . (2.21)
Les grandeurs ∂U
∂S
et T sont donc nécessairement liées de façon univoque, c’est-à-dire que ∂U
∂S
doit être une fonction de T . On peut montrer que pour retrouver les dépendances correctes
en T dans des équations d’état connues, il faut définir la température comme
37
Encadré 2.7 : Définition de la température
1 ∂S
=
T ∂U V,n
piston de surface A
⃗f = − ∂Ep ⃗ex
∂x
⃗ex
0 x
Figure 2.5 – On considère un système à gauche sur la surface duquel un système mécanique
schématisé par un ressort exerce une force.
La position de la cloison x n’est pas fixée ; c’est une variable libre macroscopique dont la
valeur à l’équilibre doit maximiser l’entropie qu’aurait le système si x était fixé. On trouve
donc la valeur d’équilibre de x en écrivant que la dérivée par rapport à x de l’entropie du
système contraint (à x fixé) est nulle, sachant que le volume V = Ax et l’énergie interne
U = E − Ep (x) dépendent de x.
∂ ∂S ∂U ∂S ∂V
S(U = E − Ep (x), V = Ax) = 0 = × + × . (2.22)
∂x ∂U V ∂x ∂V U ∂x
38
Or ∂S
∂U V
= T1 , ∂V
∂x
= A et ∂U
∂x
= − ∂E
∂x
p
(on rappelle que E est fixé). On a donc
∂Ep ∂S
− + AT =0 (2.23)
∂x ∂V U
∂S
p=T
∂V U,n
Pour une application pratique de cette formule, reprenons l’exemple du gaz sur réseau.
À partir de l’entropie du gaz (2.16), on peut calculer la pression du gaz :
p ∂S h i λV − N h N i
= = kB λ ln(λV ) + λ − λ ln(λV − N ) − λ = −λkB ln = −λkB ln 1 − .
T ∂V λV λV
On remarque en particulier que si on prend la limite λ → ∞, on obtient l’équation des gaz
parfaits :
p N
= kB . (2.24)
T V
Ce n’est pas étonnant de retrouver un gaz parfait : si λ → ∞, c’est-à-dire si le nombre de
cases disponibles par unité de volume diverge, le système devient semblable à un système
où les particules ont des positions continues et n’interagissent pas, c’est-à-dire à un gaz
parfait. Ce qui est remarquable, c’est qu’il est possible de retrouver aussi facilement cette
loi empirique de la thermodynamique à partir du seul postulat de Boltzmann, S = kB ln Ω.
(En fait, il est délicat de parler de température dans un système sans énergie. Pour compléter le calcul,
il faudrait prendre en compte l’énergie cinétique des particules ; on trouverait que Ω est multiplié par une
fonction de E et de N (pas de V !) et l’expression de p/T serait inchangée.)
∂S
µ = −T
∂n U,V
39
À l’équilibre, la température est la même dans les deux sous-systèmes, et la répartition
des particules est telle que
µ(α) = µ(β) , (2.26)
où µ(α) et µ(β) sont les potentiels chimiques dans les deux sous-systèmes.
Notons que pour un système formé d’un mélange de plusieurs espèces, l’entropie dépend
des quantités de matière de chacune des espèces :
S = S(U, V, n1 , n2 , . . .) (2.27)
où n1 , n2 , etc. sont les quantités de matière des espèces 1, 2, etc. On doit alors définir le
potentiel chimique de chaque constituant comme
∂S
µi = −T . (2.28)
∂ni U,V,nj̸=i
Dans un système constitué d’un mélange des espèces 1 et 2, séparé en deux sous-systèmes α
1 , n2 , n1 et n2 . Il y a aussi quatre
et β, il y a quatre quantités de matière intéressantes : n(α) (α) (β) (β)
membrane laissant passer l’espèce 1 mais pas l’espèce 2, il y a trois quantités de matières
fixées et constantes qui sont n(α)2 , n2 et n1 + n1 . Par contre, n1 est une variable libre
(β) (α) (β) (α)
Remarque : deux systèmes qui peuvent échanger du volume ou des particules peuvent aussi
échanger de l’énergie : en plus de l’égalité de pression ou de potentiel chimique, il y aura
également égalité de température à l’équilibre.
40
Encadré 2.11 : Identité thermodynamique
Remarques :
• Cette relation donne la différence d’entropie dS entre deux états infiniment voisins
(U, V, n) et (U + dU, V + dV, n + dn). Elle peut également être vue comme une manière
compacte d’écrire que T1 = ∂U
∂S
, p = ∂V
V,n T
∂S
U,n
et Tµ = − ∂S
∂n U,V
.
• Pour un système contenant plusieurs espèces 1, 2, 3, . . ., l’entropie dépend de plusieurs
quantités de matières S = S(U, V, n1 , n2 , n3 , . . .) et a pour différentielle
1 p µ1 µ2 µ3
dS = dU + dV − dn1 − dn2 − dn3 − · · · (2.29)
T T T T T
La température étant une grandeur positive, on voit que S est une fonction croissante de
U . Il est donc possible d’inverser S(U, V, n) vis à vis de la variable U qui peut ainsi s’exprimer
comme une fonction de S, V, n. La différentielle de U (S, V, n) est obtenue directement à partir
de l’identité thermodynamique :
Encadré 2.12 : Différentielle de U
Pour un système dont les seules variables sont n, V et U , on a
dU = T dS − p dV + µ dn. (2.30)
∂U ∂U ∂U
= T, = −p, = µ. (2.31)
∂S V,n ∂V S,n ∂n S,V
Le travail reçu par un système fermé (dn = 0) au cours d’une transformation réversible
faisant varier le volume V d’une quantité infinitésimale dV est donc
δWrév = −p dV.
Remarque : comme il est bien connu que les variables « naturelles » de l’entropie sont U , V ,
n, il nous arrivera, pour alléger, d’écrire ∂U
∂S
à la place de l’écriture précise ∂U
∂S
V,n
.
Attention toutefois lors de calcul de dérivée partielle au jeu de variables utilisé. Par
exemple, pour un gaz parfait monoatomique, on a U = 32 nRT . On ne peut cependant
! pas déduire de l’encadré précédent que p = 0 ou que µ = 32 RT . La raison est que (2.31)
donne les dérivées de la fonction U (S, V, n), à S constant, et non les dérivées de la
fonction U (T, V, n), à T constant.
41
de modifier le volume V en fournissant au système un travail δWrév , et que l’ensemble
système + dispositif mécanique est isolé et subit une transformation réversible. Alors, on
a dS = 0 et dU = δWrév , ce qui donne directement le résultat indiqué.
1 fX
dS = dU − dX + · · · ou dU = T dS + fX dX + · · ·
T T
C’est l’identité thermodynamique généralisée. En terme de dérivées partielles :
∂S ∂U
fX = −T = .
∂X U,... ∂X S,...
Le travail reçu par le système au cours d’une transformation réversible faisant varier
X d’une quantité infinitésimale dX est
δWrév = fX dX.
dU = δQ + δW, (2.32)
où δQ et δW sont les quantités d’énergie reçues par le système sous forme de chaleur et
sous forme de travail. Pour une transformation quasistatique, la valeur de dU est également
42
donnée par l’identité thermodynamique ; d’après l’encadré 2.13,
dU = T dS + δWrév . (2.33)
Dans le cas d’un fluide, δWrév = −p dV . Pour un fil, δWrév = f dL, etc. Dans des cas
complexes, δWrév peut être une somme de plusieurs contributions.
Si la transformation est de plus réversible, on a par définition δW = δWrév . On en déduit,
en comparant les deux équations précédentes, que la chaleur reçue lors d’une transformation
réversible est δQrév = T dS.
Cette expression de la chaleur est en fait plus générale :
Encadré 2.14 : Variation d’entropie et chaleur reçue
Pour une transformation infinitésimale réversible ou, plus généralement, une transfor-
mation infinitésimale quasistatique telle que δW = δWrév , on a
δQ = T dS.
Pour mettre en pratique cette démarche, nous l’appliquons à l’étude de deux systèmes
très simples. L’un, le modèle de Langmuir est présenté dans une vidéo. L’autre est le système
à deux niveaux.
43
2.7.1 Le système à deux niveaux
On considère un système isolé constitué de N particules fixes. (Par exemple, ce sont les
atomes d’un solide.) On suppose que ces particules peuvent être dans deux états :
• le fondamental, d’énergie 0,
• l’état excité, d’énergie ϵ > 0,
et on néglige l’énergie d’interaction entre particules. Par exemple, un métal paramagnétique
plongé dans un champ magnétique B ⃗ peut être vu comme un système à deux niveaux.
Chaque atome possède un moment magnétique ; dans le fondamental, ce moment est aligné
avec le champ B ⃗ et dans l’état excité il est opposé au champ B.
⃗
On cherche à déterminer le nombre Ω(N, U ) de microétats pour un macroétat caractérisé
par son nombre N de particules et son énergie U .
Les valeurs possibles de l’énergie U du système sont 0, ϵ, 2ϵ, 3ϵ, . . ., N ϵ : le système a
l’énergie U = M ϵ lorsque M particules sont dans l’état excité et N −M dans le fondamental.
Il y a un seul microétat d’énergie nulle (toutes les particules sont dans le fondamental). Il
y a N microétats d’énergie ϵ (il n’y a qu’une particule excitée, mais ce peut-être n’importe
laquelle dans le système).
De manière générale Ω(N, U ) est donné par le nombre de manière de choisir les U/ϵ
particules excitées dans le système. Cet Ω est donc donné par un binomial, comme dans la
section précédente :
N
Ω(N, U ) = , (2.34)
U/ϵ
à comparer à (2.3).
À partir de là, tous les calculs sont similaires à ceux de la section précédente, et on se
contente ici de donner quelques résultats. On calcule l’entropie S = kB ln Ω comme avant,
avec l’aide de la formule de Stirling, et on trouve
h U U U U i
S = kB N ln N − ln − N − ln N − , (2.35)
ϵ ϵ ϵ ϵ
à comparer à (2.16).
À partir de l’entropie, on peut obtenir la température T (N, U ) du système
1 ∂S kB N ϵ − U
= = ln . (2.36)
T ∂U ϵ U
On peut alors inverser cette relation pour avoir l’énergie interne U (N, T ) du système, puis
calculer la capacité thermique C = ∂U/∂T , etc. Dans le cas où le système à deux niveaux
est un système magnétique, on peut aussi déterminer l’aimantation du système en fonction
de la température et retrouver la formule de Curie : dans la limite des faibles champs B ⃗ et
des hautes températures T , l’aimantation est proportionnelle B/T .
⃗
44
Chapitre 3
Fonctions thermodynamiques
Le critère d’évolution vers l’équilibre d’un système isolé est connu : d’après le second
principe, l’entropie S du système devient maximale.
En pratique, on considère souvent des systèmes en contact avec un réservoir de chaleur,
c’est-à-dire un thermostat, ou un réservoir de chaleur et de volume, c’est-à-dire un ther-
mostat + pressostat. (D’autres types de réservoirs sont évidemment possibles.) On note, de
manière générique, le réservoir par la lettre R. Dans ces cas-là, l’entropie du système n’est
pas maximale, et on a besoin d’un nouveau critère pour déterminer l’équilibre.
La manière de faire consiste à considérer l’entropie totale de l’univers, c’est-à-dire du
système avec le réservoir :
Sunivers = S + SR . (3.1)
L’univers peut alors être considéré comme isolé, et on peut appliquer le second principe :
∆Sunivers ≥ 0, Sunivers est maximale à l’équilibre. (3.2)
Dans ce chapitre, nous allons mettre l’inégalité (3.2) sous une forme plus pratique. Ceci
nous conduira à introduire les fonctions thermodynamiques F et G, qui se révéleront être
extrêmement commodes. En particulier, les relations de Maxwell avec les fonctions F et G
et autres conduisent à des relations non-intuitives sur les propriétés thermodynamiques des
systèmes, caractérisés par les coefficients calorimétriques et thermoélastiques.
45
• Un réservoir de particules ou réservoir de matière est un système qui peut recevoir
ou fournir de la chaleur et de la matière sans que sa température ou son potentiel
chimique ne change.
dTR = 0, dµR = 0 quels que soient δQR et dnR µR : potentiel chimique du réservoir
dnR : matière reçue par le réservoir
On suppose toujours qu’un réservoir est en équilibre interne (toutes ses variables ther-
modynamiques sont donc bien définies). Les transformations subies par un réservoir
sont donc toujours supposées quasistatiques.
Remarques :
• Les réservoirs sont des systèmes idéaux, on ne peut en fabriquer que des approxima-
tions.
• En général, un réservoir physique est un « gros » système. Par exemple, l’océan fait
un excellent thermostat (si on réalise une réaction exothermique dans une éprouvette
plongée dans l’océan, l’océan va recevoir de la chaleur mais sa température ne va pas
changer). De même, l’atmosphère fait un excellent pressostat.
• Un exemple de « petit » thermostat est un mélange d’eau et de glace, qui est donc
forcément à 0◦ C : si on fournit de la chaleur à ce système, un peu de glace fond, mais
la température ne change pas. De même, quand on fait la cuisine, un « bain Marie »
est un thermostat à 100◦ C.
46
évolutions quasistatiques, on peut écrire l’identité thermodynamique appliquée au réservoir :
dSR = T1R dUR + TpRR
dVR . Comme pR et TR sont constants, on peut intégrer et obtenir
∆UR pR
∆SR = + ∆VR . (3.3)
TR TR
∆U = W − pR ∆V + Q, (3.4)
Remarques :
• Cette relation est très générale ; en particulier, on n’a nulle part supposé que la trans-
formation était quasistatique, ou réversible.
• Cette relation reste valable si le système échange une chaleur Q′ avec un autre dispositif,
en plus de la chaleur Q échangée avec le réservoir. Pour le voir, il suffit de remplacer
partout W par W + Q′ dans le calcul ci-dessus.
• Le calcul précédent se généralise facilement au cas où le système échange avec le ré-
servoir une autre quantité conservée X que le volume, par exemple des molécules, ou
des charges, ou de la surface. Dans ce cas, on remplace dans les équations précédentes
tous les termes en −pR ∆V par fX,R ∆X où fX,R est la quantité conjuguée de X (voir
encadré 2.13), constante pour un réservoir, par exemple le potentiel chimique si X est
une quantité de molécules, un potentiel électrostatique si X est une charge, une tension
de surface si X est une aire. On obtient que quelle soit la grandeur X échangée, la
variation d’entropie du réservoir est donnée par la relation de l’encadré.
• Une manière de retenir cette formule est de se dire que tout se passe pour le réservoir
comme si la transformation était réversible. Donc, d’après l’encadré 2.14, δQR = TR dS
qui s’intègre en QR = TR ∆SR car TR est constant. Comme QR (chaleur reçue par le
réservoir) vaut −Q (avec Q la chaleur reçue par le système depuis ce réservoir), on
retrouve bien l’encadré 3.2.
47
3.2.2 Inégalité de Clausius
En utilisant que l’entropie de l’ensemble isolé {système + réservoir+ dispositif mécanique}
ne peut que croître, et que l’on peut toujours supposer le dispositif mécanique comme étant
purement mécanique et donc sans entropie, on a ∆Sunivers = ∆S + ∆SR ≥ 0, et on obtient
directement la célèbre inégalité de Clausius :
Encadré 3.3 : augmentation minimale d’entropie
Pour tout système pouvant échanger de la chaleur avec un unique réservoir à la tem-
pérature TR , on a
Q
∆S ≥ ,
TR
l’égalité n’ayant lieu que pour une transformation réversible.
48
Pour une transformation sans échange de travail (W = 0), Φ ne peut que diminuer, et :
Φ = U − TR S.
Φ = U − TR S + pR V.
Remarques :
• Le potentiel thermodynamique Φ n’est pas une fonction d’état, puisqu’il dépend de
l’état du système (U , S, V ) et de celui du réservoir (TR , et pR si le réservoir est un
pressostat).
• On peut généraliser et définir un potentiel thermodynamique pour un système qui
échange avec un réservoir de l’énergie et une quantité X autre que le volume, fX,R
étant la quantité conjuguée de X constante pour le réservoir. On montre alors que
−TR ∆Sunivers = ∆ (U − TR S − fX,R X) − W ≤ 0 et on peut définir un potentiel ther-
modynamique Φ = U − TR S − fX,R X, tel que ∆Φ ≤ W . On voit donc que la définition
du potentiel thermodynamique dépend du type d’échange existant entre le système et
le réservoir.
Pour un système dont les seuls échanges sont avec un réservoir R, on a, à l’équilibre
thermodynamique,
(Ceci se généralise à d’autres quantités que le système pourrait échanger avec le réservoir.)
Montrons que T = TR à l’équilibre thermodynamique si le système échange de l’énergie
avec le réservoir. Pour cela, trouvons la valeur de U qui minimise Φ, le volume (et les autres
paramètres) étant fixés ; à l’équilibre thermodynamique, le minimum est atteint et on doit
avoir ∂U
∂Φ
V
= 0. Si le réservoir est un thermostat, Φ = U − TR S. Si c’est aussi un pressostat,
il faut ajouter le terme pR V , mais ce terme est constant puisqu’on ne fait varier que U . On
a donc, dans tous les cas,
∂Φ ∂S TR
= 1 − TR =1− , (3.8)
∂U V ∂U V T
et on retrouve bien T = TR . On peut aller plus loin en utilisant le fait que, puisque Φ est
∂2Φ
minimale à l’équilibre, on doit également avoir ∂U 2 V ≥ 0. On obtient donc
∂ 2Φ TR ∂T
= × ≥ 0. (3.9)
∂U 2 V T2 ∂U V
On a donc ∂U ∂T
V
≥ 0, et donc CV = ∂U ∂T V
≥ 0. Pour tout système thermodynamique à
l’équilibre, la capacité thermique à volume constant est positive.
49
On montre de même que pour un système qui échange chaleur et volume avec un réservoir,
on retrouve T = TR et p = pR à l’équilibre en minimisant Φ. Par ailleurs, en utilisant le
fait que Φ est minimale et en écrivant que certaines dérivées secondes sont positives, on
arrive également à montrer que ∂V ∂p T
≤ 0 : quand on augmente la pression de n’importe quel
système thermodynamique à l’équilibre à température fixée, son volume diminue. Autrement
dit, le coefficient de compressibilité isotherme χT introduit plus loin (voir encadré 3.11) est
toujours positif.
Ces trois relations avec la loi de conservation n(α) + n(β) = ntot permettent en principe de
déterminer l’équilibre.
Supposons maintenant que T et p sont fixés et égaux à TR et pR . On part d’une situation
hors équilibre thermodynamique avec µ(α) > µ(β) . Lors du retour vers l’équilibre, Φ ne peut
50
que diminuer, ce qui implique avec (3.11) que dn(α) < 0 : les particules se déplacent de la
région où le potentiel chimique est le plus grand vers la région où le potentiel chimique est
le plus faible.
F = U − T S, (3.13)
appelée énergie libre. Pour un système en contact uniquement avec un thermostat et subissant
une transformation telle que les températures initiale Ti et finale Tf du système sont égales
à TR , on a
∆F = ∆(U − T S) = ∆(U − TR S) = ∆Φ si Ti = Tf = TR .
(Cette égalité est vraie même si la température T du système varie au cours de la transfor-
mation.) Un cas particulier (mais important) de l’encadré 3.4 est donc :
Encadré 3.5 : Énergie libre
F = U − T S.
Lorsqu’un système subit une transformation au cours de laquelle les seuls échanges de
chaleur sont avec un thermostat à la température TR , et telle que les températures
initiale Ti et finale Tf sont égales à TR , on a
∆F ≤ W.
Remarque : L’inégalité de l’encadré est souvent écrite dans l’autre sens : −W ≤ −∆F . La
quantité −W est le travail fourni par le système, et représente donc le travail utile que l’on
récupère d’une machine. Ce travail est inférieur à la diminution −∆F de l’énergie libre de
la machine. L’énergie libre est donc bien l’énergie qu’il est envisageable de récupérer pour
fournir un travail utile, celle qui est « libre ».
De même, on introduit la fonction thermodynamique (fonction d’état)
G = U + pV − T S, (3.14)
appelée enthalpie libre. Pour une transformation telle que les températures initiale et finale
sont égales à TR , et les pressions initiale et finale égales à pR , on a ∆G = ∆(U −TR S +pR V ).
On obtient donc :
51
Encadré 3.6 : Système en contact avec un thermostat et un pressostat
G = U + pV − T S.
∆G ≤ W.
Remarque : Comme dans le cas précédent, l’inégalité de l’encadré est souvent écrite dans
l’autre sens, −W ≤ −∆G, pour montrer que le travail utile récupérable est limité par la
diminution de l’enthalpie libre.
dU = T dS − p dV + µ dn.
52
Les différentielles de ces fonctions se déduisent de la différentielles de U :
dF = −S dT − p dV + µ dn,
dH = T dS + V dp + µ dn,
dG = −S dT + V dp + µ dn.
53
où Uéquilibre = Uéquilibre (T, V, n) et Séquilibre = Séquilibre (T, V, n) sont les valeurs de l’énergie
interne et de l’entropie du système quand ce dernier est à l’équilibre thermodynamique à
température T , les valeurs de n et V étant fixées. Pour déterminer l’expression de F (T, V, n),
il faut d’abord calculer Uéquilibre et Séquilibre comme on l’a fait dans l’exemple ci-dessus.
La transformation permettant d’obtenir la fonction F s’appelle une transformation de
Legendre.
∂ 2f ∂ 2f
∂ ∂f ∂ ∂f
:= = =: . (3.22)
∂y ∂x ∂y ∂x y x ∂x ∂y x y ∂x ∂y
∂S ∂p
= . (3.23)
∂V T,n ∂T V,n
dF = −S d T −p d V +µ dn
Il y a une dérivée de Maxwell pour chaque paire de variables de chaque fonction ther-
modynamique. Par exemple, les relations de Maxwell associées à F sont, en plus de (3.23),
∂S
∂n T,V
∂µ
= − ∂T V,n
et ∂n
∂p
T,V
∂µ
= − ∂V T,n
.
54
Pour une fonction extensive quelconque, la propriété d’extensivité s’exprime mathéma-
tiquement de la façon suivante : si on multiplie par λ la taille du système, c’est-à-dire si on
multiplie par λ chacune des variables extensives, alors la valeur de la fonction est multipliée
par λ. En prenant les exemples de F et G, on a donc
F (T, λV, λn) = λ F (T, V, n), G(T, p, λn) = λ G(T, p, n). (3.24)
V
F (T, V, n) = n f (T, v) avec v = et f (T, v) = F (T, v × 1 mol, 1 mol). (3.25)
n
Ici, v est le volume molaire et f est l’énergie libre molaire.
La même opération pour G donne G(T, p, n) = n g(T, p), avec g l’enthalpie libre molaire,
mais cette relation n’est presque jamais écrite sous cette forme. En effet, si on la dérive par
rapport à n, on voit que g = ∂G ∂n T,p
= µ, et on écrit plutôt G(T, p, n) = n µ(T, p). En écrivant
la différentielle de cette dernière relation, on obtient
c’est-à-dire que le potentiel chimique est égal à l’enthalpie libre molaire. C’est la relation
d’Euler appliquée à G.
Toujours en supposant un corps pur sous une seule phase, le potentiel chimique ne
dépend que de T et p et sa différentielle est donnée par
dµ = −s dT + v dp,
avec s = S
n
l’entropie molaire et v = V
n
le volume molaire. C’est la relation de Gibbs-
Duhem.
3.3.5 Généralisations
La relation fondamentale de la thermodynamique change selon le système étudié ; celle
donnée dans l’encadré 3.7 est pour un corps pur sous une seule phase. Dans un cas plus
compliqué, il peut y avoir des termes supplémentaires ; par exemple, pour un mélange de
deux espèces, il faut remplacer µ dn par µ1 dn1 + µ2 dn2 dans U , et donc dans les fonctions
F , G et H. Dans ce cas particulier important qui est au cœur du chapitre 4, la propriété
d’extensivité pour G s’écrit
55
la relation d’Euler associée devient (exercice)
G = µ1 n1 + µ2 n2 . (3.28)
dU = T dS − p dV + µ dn + fX dX .
56
Cette égalité et les suivantes ne sont évidemment valables qu’en un point où p = psol,liq (T ).
La relation de Gibbs-Duhem donne les dérivées de µ et permet d’écrire
dpsol,liq dpsol,liq
vsolide × − ssolide = vliquide × − sliquide . (3.31)
dT dT
En regroupant,
dpsol,liq sliquide − ssolide lfusion
= = , (3.32)
dT vliquide − vsolide T (vliquide − vsolide )
où lfusion = T (sliquide − ssolide ) est l’enthalpie molaire de fusion, voir section 1.7.3. C’est la
formule de Clapeyron.
Application numérique à faire en exercice : pour de l’eau près de 0°C en utilisant les
masse volumique de la glace Msolide = 0,92 g · cm−3 et de l’eau liquide Mliquide = 1 g · cm−3 ,
et en utilisant l’enthalpie massique de changement de fusion de l’eau Lfusion = 334 kJ · kg−1 ,
dpsol,liq
on trouve dT ≈ −1,4 107 Pa · K−1 , ce qui veut dire qu’il faut augmenter la pression de
140 bar pour que la température de fusion baisse de 1°C.
dU = T dS − p dV, dF = −S dT − p dV,
(3.33)
dH = T dS + V dp, dG = −S dT + V dp.
57
Encadré 3.10 : Les coefficients calorimétriques
CV ℓ
dS = dT + dV ⇐⇒ δQrév = T dS = CV dT + ℓ dV. (3.34)
T T
Avec l’identité thermodynamique (3.33), on obtient
dU = CV dT + (ℓ − p) dV, (3.35)
et donc
∂S ∂U ∂S ∂U
CV = T = et ℓ=T =p+ . (3.36)
∂T V ∂T V ∂V T ∂V T
Cp h
dS = dT + dp ⇐⇒ δQrév = T dS = Cp dT + h dp. (3.37)
T T
Avec l’identité thermodynamique (3.33), on obtient
dH = Cp dT + (h + V ) dp, (3.38)
et donc
∂S ∂H ∂S ∂H
Cp = T = et h=T = −V + . (3.39)
∂T p ∂T p ∂p T ∂p T
Attention : dS est une différentielle, mais T dS = δQrév n’en n’est pas une.
Pour une transformation isochore (dV = 0) du fluide, on a δW = 0 et (3.35) donne
Z Tf
δQ = dU = CV dT, ⇔ Q = ∆U = CV dT. (3.40)
Ti
Pour une transformation isobare (dp = 0) d’un fluide en permanence à l’équilibre méca-
nique, on a δW = −p dV , dH = dU + p dV = δQ et (3.38) donne
Z Tf
δQ = dH = Cp dT, ⇔ Q = ∆H = Cp dT. (3.41)
Ti
58
ne présentons ici que la plus efficace, qui repose sur la remarque suivante. Faisant le choix
d’un jeu de variables, il est judicieux de trouver la « bonne » fonction d’état, celle dont la
différentielle exprime l’identité thermodynamique dans les variables retenues. Les relations
de Maxwell écrites avec cette fonction d’état aboutissent alors aux calculs les plus simples.
CV et ℓ en fonction de (T, V )
On veut exprimer ℓ = T ∂V ∂S
T
, voir (3.36). D’après (3.33), il apparaît clairement que
la « bonne » fonction d’état est l’énergie libre. La relation de Maxwell associée à dF =
−S dT − p dV donne
∂S ∂p
= . (3.42)
∂V T ∂T V
D’où immédiatement
∂p
ℓ=T , (3.43)
∂T V
qui est la première relation de Clapeyron. Ensuite, écrivant la relation de Maxwell sur la
différentielle dS donnée en (3.34) ;
∂(CV /T ) ∂(ℓ/T ) ∂CV ∂ 2p
= =⇒ =T , (3.44)
∂V T ∂T V ∂V T ∂T 2 V
Cp et h en fonction de (T, p)
De même, mais en partant cette fois de la fonction G, on obtient (exercice) les deux
relations de Clapeyron suivantes :
∂V ∂Cp ∂ 2V
h = −T , = −T , (3.45)
∂T p ∂p T ∂T 2 p
qui permettent en particulier d’obtenir la dépendance de Cp en pression à partir de l’équation
d’état. Par exemple, on obtient que le Cp du gaz parfait ne dépend pas de p, mais seulement
de T (et n).
1 ∂V
χT = − en Pa−1 ,
V ∂p T
59
qui indique à quel point il est facile ou difficile de comprimer le fluide en aug-
mentant la pression, tout en maintenant la température constante. χT est une
quantité intensive, et positive, comme nous l’avons déjà mentionné section 3.2.3.
• Le coefficient de compressibilité isentropique (ou compressibilité adiabatique)
1 ∂V
χS = − en Pa−1 ,,
V ∂p S
qui indique à quel point il est facile ou difficile de comprimer le fluide en augmen-
tant la pression, tout en maintenant l’entropie constante, c’est-à-dire sans échange
de chaleur. χS est une quantité intensive et nous allons voir que χT ≥ χS ≥ 0.
• Le coefficient de dilatation thermique
1 ∂V
α= en K−1 ,
V ∂T p
qui indique comment le volume évolue avec la température. α est une quantité
intensive qui est généralement positive, mais pas toujours. (Par exemple, si on
étend la définition de α aux solides isotropes, alors le α du caoutchouc est négatif.)
Par exemple, pour le gaz parfait, on obtient facilement χT = 1/p et α = 1/T à par-
tir de l’équation d’état. Pour de l’eau liquide aux conditions usuelles, on a χT ≈ χS ≈
4,5 10−10 Pa−1 ; l’eau est donc environ 20 000 fois moins compressible qu’un gaz. Le coefficient
de dilatation thermique de l’eau liquide à pression ambiante varie de manière significative
avec la température : il vaut 2 10−4 K−1 à 20°C (environ 15 fois moins que pour un gaz à la
même température), 7 10−4 K−1 à 100°C, il s’annule à 4°C et est négatif entre 0 et 4°C ! Un
kilogramme d’eau occupe 1 ℓ à 4°C, mais 1,04 ℓ à 100°C, soit une augmentation relative de
4% = α ∆T avec ∆T ≈ 100°C et α ≈ 4 10−4 K−1 .
Pour bien comprendre la différence entre χT et χS :
• Lorsqu’on comprime un gaz, on fournit un travail et donc une énergie. Si on com-
prime suffisamment rapidement le gaz, il n’a pas le temps d’évacuer de la chaleur vers
l’extérieur et sa température augmente. La compression doit alors être traitée comme
adiabatique et, en supposant qu’elle reste quasiment réversible, elle sera également
isentropique puisque T dS = δQ = 0. Le coefficient pertinent sera χS . Comme expliqué
dans une vidéo, ce coefficient apparaît dans l’expression de la vitesse du son dans un
fluide.
• Si on comprime le gaz suffisamment lentement pour que l’équilibre thermique soit à
tout instant réalisé, la compression est isotherme et il faut considérer le coefficient χT .
Il existe beaucoup de relations entre tous les coefficients que nous venons de définir ; dans
ce paragraphe, nous en présentons les principales.
60
Relation de Mayer
La relation de Mayer est la relation entre les deux capacités calorifiques qui ont été
définies, à volume constant CV et pression constante Cp . On part de la différentielle du
volume
∂V ∂V
dV = dT + dp = αV dT − χT V dp, (3.46)
∂T p ∂p T
où l’on a utilisé les définitions de α et χT de l’encadré 3.11. On injecte dans (3.34) pour
arriver à
(3.47)
T dS = CV dT + ℓ αV dT − χT V dp = CV + ℓ αV dT − ℓ χT V dp.
La première équation est l’une des formes de la relation de Mayer. Pour un gaz parfait,
α = T1 , ℓ = p et on obtient Cp − CV = nR.
On définit l’indice adiabatique γ comme le rapport des capacités calorifiques
Cp
γ= . (3.49)
CV
Avec la relation de Mayer (3.48), le rapport γ permet d’exprimer CV et Cp :
1 γ
CV = ℓαV et Cp = ℓαV. (3.50)
γ−1 γ−1
La relation de Mayer a plusieurs formes. Par exemple, on peut éliminer ℓ dans (3.48). On
part de la relation
∂p ∂T ∂V
× × = −1. (3.51)
∂T V ∂V p ∂p T
(Propriété mathématique générale valable dès que trois quantités sont reliées entre elles, voir
section A.1.3.) Le premier facteur est ℓ/T , voir (3.43). Les deuxième et troisième facteurs
sont respectivement αV
1
et −V χT , voir encadré 3.11. On en déduit que ℓ = αT
χT
et, dans (3.48),
que 2
α2 T V
T ∂V
Cp − CV = = . (3.52)
χT V χT ∂T p
Comme χT ≥ 0, cette expression implique que Cp ≥ CV .
Formule de Reech
On peut calculer le rapport des deux compressibilités χT et χS . Pour une transformation
quasistatique et isentropique, on a
T dS = CV dT + ℓ dV = Cp dT + h dp = 0 (3.53)
et donc
ℓ h
dT = − dV = − dp. (3.54)
CV Cp
61
On s’intéresse à la relation entre dV et dp, établie en supposant dS = 0, et on obtient
∂V CV h CV
= =− V χT , (3.55)
∂p S Cp ℓ Cp
χT Cp
= = γ, (3.56)
χS CV
62
Chapitre 4
Mélanges
Lorsqu’un corps est pur et sous une seule phase, on rappelle que, par extensivité, son
potentiel chimique peut s’écrire comme une fonction de la température T et de la pression
p, et qu’il vérifie la relation de Gibbs-Duhem, voir l’encadré 3.9 :
Encadré 4.2 : Relation de Gibbs-Duhem
Pour un corps pur sous une seule phase
µ = µ(T, p) et dµ = −s dT + v dp,
63
Bien évidemment,
x1 + x2 + · · · = 1. (4.3)
L’étude thermodynamique d’un mélange, même sous une seule phase, est plus compliquée
que celle d’un corps simple :
• Une fonction thermodynamique dépend de plus que trois variables ; par exemple, l’en-
thalpie libre G dépend de T , de p, et de toutes les quantités de matières n1 , n2 , etc.
• Il y a autant de potentiels chimiques qu’il y a d’espèces. On notera
∂G ∂G
µ1 = , µ2 = , etc. (4.4)
∂n1 T,p,les ni autres que n1 ∂n2 T,p,les ni autres que n2
64
Encadré 4.3 : Pressions partielles
∂F ∂F1
µ1 = = (4.9)
∂n1 V,T,n2 ∂n1 V,T
avec µ01 le potentiel chimique d’un gaz parfait constitué de l’espèce 1 seule.
Pour aller plus loin, nous avons besoin d’en savoir plus sur le potentiel chimique µ0
d’un gaz parfait constitué d’une seule espèce. La relation de Gibbs-Duhem (3.9) donne
∂µ0
∂p T
= V /n = RT /p pour un gaz parfait. En intégrant par rapport à p en gardant T fixé,
nous trouvons
µ0 (T, p) = RT ln(p) + f (T ), (pour un gaz parfait) (4.11)
où la constante d’intégration f (T ) est indépendante de p, mais dépend de T et de la nature
du gaz. Cette expression permet d’obtenir µ01 (T, x1 p) = µ01 (T, p) + RT ln(x1 ) et, finalement :
Encadré 4.4 : Mélange idéal
65
Remarques :
• Pour un mélange idéal de liquides, par exemple, le µ0i (T, p) dans (4.12) est bien sûr le
potentiel chimique du liquide no i, pas le potentiel chimique d’un gaz parfait.
• Dans un mélange de gaz réels, si les interactions sont négligeables, on pourra bien
sûr souvent considérer que le mélange est idéal (limite des gaz parfaits). Plus surpre-
nant, beaucoup de mélanges de liquides dont les molécules ont des tailles et structures
similaires sont bien décrits par l’équation du mélange idéal.
• Exercice : à partir de (4.12), montrer que si on mélange (à p et T fixés) un volume V1
de l’espèce 1 pure et un volume V2 de l’espèce 2 pure et que le mélange est idéal, alors
le volume total est V1 + V2 .
• Si on mélange 1 ℓ d’eau et 1 ℓ d’éthanol, le mélange fait 1,92 ℓ. Ce mélange n’est donc
pas idéal.
66
à comparer avec (3.29). En combinant ces dernières équations, on obtient µ0G (T, ps (T )) =
µ0G (T, pvapeur ) et donc
ps (T )
pvapeur = ps (T ) =⇒ x= = xs (T, patm ). (4.17)
patm
Par exemple, ps (20°C) = 0,023 bar pour l’eau. À l’équilibre, sous pression atmosphérique
p = 1 bar et à cette température de 20°C, l’eau s’évapore donc avec une fraction molaire de
2,3 % de vapeur. À T = 0 °C, ps (0°C) = 0,006 bar et la fraction molaire tombe à 0,6 %. Ces
valeurs montrent que, si de l’air chargé en humidité à l’équilibre se refroidit, la vapeur en
excès se liquéfie sous forme de fines gouttelettes, formant le brouillard ou la buée.
Lorsque l’évaporation a lieu en milieu ouvert, la vapeur produite par le liquide est diluée
ou évacuée, de sorte que l’équilibre n’est jamais atteint, et l’évaporation se poursuit jusqu’à
l’assèchement du liquide. En dehors de l’équilibre, on définit le degré hygrométrique ϕ de la
phase gazeuse (supposée homogène), ou taux d’humidité relative, comme étant le rapport
du taux observé x = pvapeur /patm sur le taux à l’équilibre xs (T, patm ) = ps (T )/patm :
x pvapeur
ϕ= = ⩽ 1. (4.18)
xs ps (T )
67
exacte, et la pression du système doit être égale à la pression de vapeur saturante ps (T ) du
corps. La température d’ébullition s’obtient donc par
ps (T ) = patm
68
Le potentiel chimique du soluté vérifie
69
où le volume molaire du solvant V /n1 était noté v dans (4.27). En mettant les trois dernières
équations dans (4.23), on trouve
0 V n2 ′ n1
µ1 (T, p, x2 ) ≈ µ1 (T, p) + (p1 − p) + fint ,T . (4.24)
n1 V V
On cherche maintenant à estimer p1 . Pour cela, il suffit d’utiliser p = −∂F/∂V ; en dérivant (4.22)
par rapport à V et en multipliant par -1, on obtient
n1 n2 ′ n1
p = p1 + p2 + f ,T ,
V 2 int V
où p et p1 ont déjà été définis et où p2 est la pression qu’aurait le soluté seul dans le volume V .
Cette relation entre les pressions permet de simplifier (4.24) ; on obtient
V
µ1 (T, p, x2 ) ≈ µ01 (T, p) − p2 .
n1
Mais p2 est la pression du soluté seul, c’est-à-dire la pression d’un gaz parfait : p2 V = n2 RT . En
écrivant n2 /n1 ≈ x2 , on retrouve (4.20).
Passons maintenant au calcul de µ2 . On dérive (4.22) par rapport à n2 et on trouve
n n
1 1
µ2 (T, p, x2 ) = µgaz parfait (T, p2 ) + fint , T = RT ln(p2 ) + f (T ) + fint ,T , (4.25)
V V
où l’on a utilisé (4.11) pour exprimer le potentiel chimique du gaz parfait. Le volume du mé-
lange est une fonction de p, T , n1 et n2 . En fait, comme n2 ≪ n1 , on a au premier ordre
V = V (p, T, n1 , n2 ) ≈ V (p, T, n1 , 0) = n1 v(p, T ) par extensivité. On a alors p2 = n2 RT /V ≈
n2 RT /[n1 v(p, T )] ≈ x2 RT /v(p, T ). En reportant ces expressions de p2 et de V dans (4.25), on
1
retrouve bien le résultat attendu (4.21) avec Ψ (T, p) = f (T ) + fint v(p,T ) , T + RT ln[RT /v(p, T )].
soluté 2 dans les deux compartiments. En supposant les deux systèmes dilués, (4.20) donne
µ01 (T, pA ) − RT xA 0 B
2 = µ1 (T, pB ) − RT x2 . (4.26)
On peut relier les deux potentiels chimiques en utilisant le fait que dans un liquide ordinaire,
le volume molaire v dépend peu de la pression (c’est l’incompressibilité). On part de la
∂µ0
relation de Gibbs-Duhem (3.9) : ∂p1 T = v. En intégrant cette relation entre pA et pB en
supposant v constant, on obtient
µ01 (T, pB ) − µ01 (T, pA ) = (pB − pA )v. (4.27)
Dans (4.26), ceci permet d’obtenir :
70
Encadré 4.6 : Pression osmotique — loi de van’t Hoff
trouve alors une pression osmotique de 1,4 bar, ce qui est loin d’être négligeable. Si la mem-
brane de A n’est pas très solide, il y a risque de rupture.
Par définition de T , les potentiels chimiques pour le solvant 1 dans les phases liquide et
∗
µL1 (T ∗ , p, xL2 ) = µG ∗ G
1 (T , p, x2 ).
µL,0 ∗ ∗ L G,0 ∗ ∗ G
1 (T , p) − RT x2 = µ1 (T , p) − RT x2 ,
∂µα,0
µα,0 ∗
1 (T , p) ≈ µα,0
1 (T0 , p) + 1
(T ∗ − T0 ) = µα,0 α ∗
1 (T0 , p) − s1 (T0 , p)(T − T0 ),
∂T p
où sα1 (T0 , p) est l’entropie molaire du solvant pur dans la phase α. Puis, en utilisant que
1 (T0 , p) = µ1 (T0 , p) (parce qu’on a supposé que T0 était la température d’équilibre
µL,0 G,0
liquide/gaz du solvant pur à la pression p) et que sG 1 (T0 , p) − s1 (T0 , p) = lvap /T0 , avec lvap
L
l’enthalpie molaire (la chaleur latente molaire) de vaporisation, voir la section 1.7.3, on trouve
2 L G
RT 0 x 2 − x 2
T ∗ − T0 = . (4.29)
lvap
71
Si xG2 ≪ x2 , nous trouvons que la température d’équilibre liquide/vapeur à pression donnée
L
augmente. En particulier, la température d’ébullition (voir section 4.1.2) augmente elle aussi.
C’est la loi de Raoult de l’ébullioscopie.
Nous pouvons reprendre le même type de raisonnement pour la fusion du solvant (tran-
sition de phase entre un solide et un liquide). Si on suppose que xS2 ≪ xL2 , la température
de fusion diminue selon la loi de Raoult de la crysoscopie :
RT02 xL2
T ∗ − T0 = − . (4.30)
lfusion
Prenons l’exemple de l’eau à pression atmosphérique. Si on ajoute 16 g de sel (NaCl) dans
1 ℓ d’eau, on trouve que la fraction molaire de soluté dans le liquide est xL2 = 10−2 (un pour
cent de soluté). Pour l’abaissement de la température de congélation, nous trouvons 1◦ C.
Pour l’augmentation de la température de vaporisation, nous trouvons 0,3◦ C.
de l’espèce i ∈ {1, 2} dans la phase α ∈ {L, G}. Là encore, ce choix n’est pas très judicieux.
Il vaut mieux commencer par choisir deux variables extensives qui décrivent les tailles des
deux phases : la quantité totale de liquide nLtot = nL1 + nL2 et la quantité totale de gaz
nGtot = n1 + n2 ; les propriétés physiques du système ne dépendent pas de ces deux variables.
G G
Il reste quatre variables à choisir ; assez naturellement on peut choisir T , p et les fractions
molaires de l’espèce 2 dans les deux phases : xL2 = nL2 /nLtot et xG 2 = n2 /ntot . Cependant, ces
G G
72
Pour résumer :
Encadré 4.7 : Variance d’un mélange binaire
• Un mélange binaire sous une seule phase peut être décrit par une variable
extensive ntot et trois variables intensives T , p et x2 . Ces trois variables sont
indépendantes ; on dit que la variance du système est v = 3.
• Un mélange binaire en équilibre sous deux phases α et β peut être décrit
par deux variables extensives nphase
tot
α
et nphase
tot
β
et quatre variables intensives T ,
p et x2
phase α
et x2
phase β
. Cependant, seules deux de ces variables intensives sont
indépendantes ; on dit que la variance du système est v = 2.
Remarques :
• Si par exemple on impose la température et la pression dans un mélange eau-éthanol
sous phases liquide et gazeuse, alors les compositions des deux phases sont déterminées
par l’équilibre thermodynamique. Si on décide plutôt d’imposer la température et la
composition du liquide, alors la pression et la composition du gaz sont déterminés par
l’équilibre thermodynamique.
• La règle générale, dite règle de Gibbs, est que la variance d’un système de c composantes
sous ϕ phases est v = 2 + c − ϕ. Ainsi, pour un corps pur (c = 1) sous une seule phase
(ϕ = 1), la variance est v = 2 et l’on peut librement choisir les deux variables intensives
T et p. Si le corps pur est sous deux phases, la variances est v = 1 : on peut encore
choisir la température T , mais la pression est alors déterminée par le système (c’est la
pression de vapeur saturante pour un équilibre liquide-gaz). Si le corps pur est présent
simultanément sous trois phases, la variance est v = 0 : le système impose la valeur de
la température et de la pression ; c’est le point triple.
Cela signifie que si on est à un équilibre liquide-vapeur d’un mélange eau+éthanol à une
pression et à une température donnée, alors on n’a pas le choix de la proportion d’éthanol
dans chacune des phases. Ces proportions sont déterminés par l’équilibre thermodynamique.
On fixe la pression p une bonne fois pour toute (typiquement, p = 1 atm) et on s’intéresse
au comportement du système en fonction de la température. On trace les deux fonctions
xL2 = xL2 (T, p) et xG
2 = x2 (T, p) mais en mettant x2 en abscisse et T en ordonnée ; voir un
G
73
T
région I: gaz
T2
G
x2
L
x2
T1
région II: liquide
x2
0 1
Figure 4.1 – Un diagramme de phase isobare d’un mélange binaire typique. Pour la pression
considérée, T1 est la température de coexistence liquide/gaz de l’espèce 1 seule et T2 est la
température de coexistence liquide/gaz de l’espèce 2 seule.
reste la dernière goutte de liquide, dont la fraction molaire de l’espèce 2 est xL2 , l’abscisse du
point CL . À plus haute température, le système est totalement en phase gazeuse (point D).
Remarque : contrairement à ce qui se passe pour un corps pur, la température du mélange
binaire évolue pendant le changement de phase.
Pour un système de n moles dans la région de coexistence liquide/gaz, par exemple au
point M de la figure 4.2, on peut facilement calculer le nombre nL de moles dans la phase
74
T
D gaz
C
CL
MG M
ML
BG
B
liquide
A
x2
0 1
liquide. Il suffit d’écrire que le nombre total de moles de l’espèce 2 est d’une part donnée par
x2 n et d’autre part donné par xL2 nL + xG 2 (n − nL ), où n − nL est bien sûr le nombre nG de
moles dans la phase gazeuse. On obtient la règle des moments :
Encadré 4.8 : Règle des moments
nL x 2 − xG distance entre M et MG
= L 2
= , (4.33)
n x2 − xG2 distance entre ML et MG
2.
où les points ML , M et MG indiqués sur la figure 4.2 ont pour abscisses xL2 , x2 et xG
75
T (en °C)
100
p = 1 atm
98
96 gaz
94
92
90
88
86
84 L1 G point
• • 1
82 azéotrope
liquide
80 L2
• •
78
A1 • A2 •
76 xéthanol
0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1
eau éthanol
pure pur
76
Chapitre 5
Système Thermostat
Rappelons qu’un thermostat est un système idéal, supposé tellement grand que sa tempéra-
ture ne change pas, même lorsqu’il reçoit (ou donne) de l’énergie sous forme de chaleur (voir
encadré 3.1). En particulier, le thermostat impose sa température au système. On consi-
dère que le système ne peut échanger de la chaleur qu’avec le thermostat, et que l’univers,
c’est-à-dire l’ensemble {système + thermostat}, est isolé.
L’objet d’étude est le système, et on s’intéresse aux microétats du système (et non ceux
du thermostat). Comme le système n’est pas isolé, son énergie peut fluctuer et tous les
microétats, quelle que soit leur énergie, sont accessibles. (Par comparaison, en microcano-
nique, l’énergie du système a une valeur constante fixée par l’expérimentateur, et seuls les
microétats ayant cette énergie sont accessibles.)
On note les microétats du système par la lettre σ, et on utilise les notations suivantes :
(
Eσ l’énergie du système lorqu’il est dans le microétat σ,
Pour un microétat σ donné :
Pσ la probabilité que le système soit dans le microétat σ.
77
5.1.2 Le facteur de Boltzmann
Le résultat fondamental, qui sera démontré à la fin de ce chapitre, est le suivant :
Encadré 5.1 : Le facteur de Boltzmann
Pour un système à l’équilibre thermodynamique en contact avec un thermostat à la
température T ,
Eσ
Pσ est proportionnel à exp − . (5.1)
kB T
Remarques :
• L’expression (5.1) est bien sûr valable pour un système thermodynamique constitué
d’un grand nombre de particules. De manière peut-être surprenante, elle reste valable
pour un système, composé d’un petit nombre de particules, voire d’une seule particule !
La raison est que même si le système est « petit », l’ensemble {système + thermostat}
est « grand ».
• Dans l’expression (5.1), le thermostat n’intervient que par l’intermédiaire de sa tem-
pérature T ; la manière dont il est construit n’a aucune influence sur les propriétés du
système.
• Pour un système quantique, les microétats que l’on considère sont, comme en microca-
nonique, les états propres de l’hamiltonien. C’est un ensemble discret et on peut définir
la probabilité (5.1) de chacun de ces états.
• Pour un système classique, l’ensemble des microétats est continu. Il faut alors remplacer
la probabilité Pσ par une densité de probabilité ; voir la section 5.2.5 pour plus de
détails.
• En physique statistique, on introduit souvent, pour simplifier les notations, la tempé-
rature inverse
1
β= . (5.2)
kB T
• Le facteur de proportionnalité est souvent noté 1
Z
, où Z s’appelle la fonction de parti-
tion. L’équation (5.1) s’écrit alors
1 −βEσ
Pσ = e . (5.3)
Z
La fonction de partition Z dépend de la température T du thermostat, des propriétés
du système mais, bien évidemment, elle ne dépend pas de σ. Nous aborderons en grand
détail les propriétés de la fonction Z dans le chapitre suivant.
5.2 Applications
5.2.1 Le système à deux niveaux
Le système le plus simple que l’on puisse imaginer en physique statistique est le système
à deux niveaux (ou deux états), c’est-à-dire un système où il n’y a que deux microétats dif-
férents possibles. Ces deux niveaux peuvent par exemple être les deux orientations possibles
du moment cinétique selon l’axe des z pour une particule de spin 1/2, ou les états symé-
trique et antisymétrique d’une molécule d’ammoniac (NH3 ), etc. Selon le système considéré,
78
on peut noter ces deux niveaux {↑, ↓}, ou {+, −}, ou {A, S}, etc. Ici, de manière générique,
on décide de les noter tout simplement {1, 2} :
σ ∈ {1, 2}.
et donc
e−βE1 e−βE2
P1 = , P 2 = . (5.6)
e−βE1 + e−βE2 e−βE1 + e−βE2
Une fois qu’on a les probabilités des états, on peut calculer des propriétés statistiques comme
l’énergie moyenne dans le système :
E1 e−βE1 + E2 e−βE2 E2 − E1
⟨E⟩ = E1 P1 + E2 P2 = = E1 + . (5.7)
e−βE1 + e−βE2 1 + eβ(E2 −E1 )
L’énergie moyenne ⟨E⟩, pour une particule unique, n’est jamais l’énergie du système : celle-
ci, à un moment donné, est soit E1 , soit E2 . Cependant, si on fait de nombreuses mesures
de cette énergie et qu’on moyenne ces mesures, on obtient une quantité qui tend vers ⟨E⟩.
La moyenne d’une somme est égale à la somme des moyenne. La moyenne de n’importe quel
terme du membre de droite est donnée par (5.7). On en déduit que l’énergie moyenne totale
est
E2 − E1
énergie totale = N énergie d’une particule = N E1 + . (5.8)
1 + eβ(E2 −E1 )
L’expression (5.8) est une moyenne, et l’énergie du système à un instant n’est jamais exac-
tement égale à cette valeur. Cependant, si N est très grand, la loi des grands nombres nous
79
apprend que la différence entre l’énergie du système et sa valeur moyenne est négligeable.
En thermodynamique, avec N ∼ 1020 , la différence entre l’énergie du système et sa valeur
moyenne est quasiment impossible à mesurer, et on identifie (5.8) avec l’énergie interne du
système :
E2 − E1
U = N E1 + . (5.9)
1 + eβ(E2 −E1 )
Plus de détails sur les fluctuations d’énergie sont donnés dans le chapitre suivant.
L’énergie interne obtenue est une fonction croissante de la température (rappelons que
β = kB1T ) et, en dérivant par rapport à T , on obtient la capacité thermique C du système.
Nous vous invitons à faire le calcul et à vérifier que
C
N kB
0,4
C X 2 eX |E2 − E1 |
= avec X = . (5.10)
N kB (1 + eX )2 kB T 0,2
0
0 0,5 1 1,5 kB T
|E2 − E1 |
• En physique statistique, on peut considérer des systèmes « petits » (avec une seule
particule, ou quelques-unes). On peut alors calculer des probabilités et des valeurs
moyennes, comme dans (5.7).
• Pour faire de la thermodynamique, il faut que le système considéré soit « grand »
! (avec 1020 particules, par exemple). On peut alors calculer l’énergie interne ou la
capacité thermique, comme dans (5.9) ou (5.10).
Calculer des quantités thermodynamiques (énergie interne, entropie, capacité ther-
mique, etc.) pour un système « petit » n’aurait pas de sens.
80
la limite basse température, on dit parfois que le système est « gelé » : toutes les particules
sont « coincées » dans le fondamental, et très peu d’entre elles arrivent à s’exciter.
Comme expliqué dans une vidéo, les notions que nous venons d’introduire se généralisent
à n’importe quel système quantique. On dira qu’un système quelconque est dans sa limite
basse température, ou qu’il est gelé, si kB T est très petit devant le gap, c’est-à-dire la différence
d’énergie entre le fondamental et le premier état excité.
Les notions de « basse température » et « haute température » dépendent du système
considéré. Par exemple :
• Dans un modèle simple de l’atome, le nuage électronique est dans sa limite basse
température si (typiquement) T ≪ 105 K, c’est-à-dire : toujours.
• L’aimantation d’un métal paramagnétique est dans sa limite basse température
si (typiquement) T ≪ 1 K, même pour un champ très fort, c’est-à-dire : jamais.
! Ces notions dépendent même de quelle propriété du système on regarde. Par exemple,
pour un gaz diatomique, en plus de son déplacement du centre de masse, une particule
peut tourner et vibrer. On trouve que
• Les propriétés liées à la rotation sont dans la limite basse température si (typi-
quement) T ≪ 100 K, c’est-à-dire : presque jamais.
• Les propriétés liées aux vibrations sont dans la limite basse température si (typi-
quement) T ≪ 2 000 K, c’est-à-dire : toujours dans les conditions habituelles.
La difficulté ici est que les variables utilisées pour décrire la particule varient continûment.
Il faut alors renoncer à donner la probabilité de chaque microétat, et plutôt raisonner en
terme de densité de probabilité. Pour une seule variable q, on définit la densité f (q) de manière
suivante : quels que soient a < b, on a
Z b
f (q) dq = probabilité que la variable aléatoire q soit entre a et b.
a
f (q) dq = probabilité que la variable aléatoire q soit entre les nombres q et q + dq.
81
Ici, ces deux intégrales portent sur le domaine où R ∞peut varier la variable aléatoire
RL q. Par
exemple, selon les cas, ces intégrales peuvent R 2πêtre −∞ si q est une vitesse, ou 0 si q est une
coordonnée dans une boîte de taille L, ou 0 si q est un angle.
Dans le problème qui nous intéresse, il y a six variables. On introduit une fonction densité
fσ = f (⃗r,⃗p) = f (x, y, z, px , py , pz ) telle que
Sans surprise, le résultat fondamental de l’encadré 5.1 s’adapte directement aux systèmes
classiques avec des variables continues :
Encadré 5.2 : Le facteur de Boltzmann pour des variables continues
Remarques :
• En notant 1/C la constante de proportionnalité et β = 1/(kB T ), on peut donc écrire
1 −βEσ
fσ = e , (5.13)
C
à comparer à (5.3) dans le cas discret : Pσ = Z1 exp(−βEσ ), avec Z la fonction de
partition du système.
• Le chapitre suivant est consacré à l’étude des fonctions de partitions ; en particulier,
nous introduirons la fonction de partition Z d’un système classique et nous verrons
qu’elle est proportionnelle, mais pas égale à la quantité C dans (5.13).
• La quantité C est dimensionnée. En effet, dans l’exemple d’une particule, fσ d⃗r d⃗p est
une probabilité, et donc sans dimension. Cela implique que, dans cet exemple, C a la
dimension de d⃗r d⃗p, à savoir L3 × (M · L · T −1 )3 . Si on avait choisi comme variables la
position et la vitesse (plutôt que l’impulsion), la dimension de C aurait été L3 (L·T −1 )3 .
Pour un système de deux particules (par exemple), la dimension de C aurait encore
été différente.
• La fonction de partition d’un système, quantique ou classique, est une quantité sans
dimension, indépendante du choix des variables utilisées. On voit donc bien que C ne
peut pas être égal à Z.
Prenons un exemple simple pour illustrer l’encadré 5.2. Dans une boîte cubique de côté
L, on considère une particule classique non-relativiste de masse m soumise au champ de
gravité. Son énergie est
⃗p 2
Eσ = E(x, y, z, px , py , pz ) = mgz + , (5.14)
2m
82
où x, y et z sont dans l’intervalle [0, L].
On applique (5.12), avec β = kB1T et en utilisant « ∝ » pour « est proportionnel à » :
p2
⃗
fσ = f (x, y, z, vx , px , py , pz ) ∝ e−βmgz e−β 2m . (5.15)
Du point de vue des probabilités, cela signifie que la position et l’impulsion de la particule
sont des variables aléatoires indépendantes : on peut déterminer d’abord la position grâce à
la densité fpos puis, de manière indépendante de cette position, déterminer l’impulsion grâce
à la densité fimp .
La densité fimp est examinée en détail dans la section suivante. Ici, intéressons-nous à la
densité fpos .
La première chose qu’on remarque est que fpos ne dépend que de z. Cela signifie que les
variables x et y sont choisies de manière uniforme dans l’intervalle [0, L] : aucune valeur
n’est favorisée. Par contre, la densité de probabilité fpos décroît avec l’altitude z : a priori,
les particules ont plus tendance à être en bas plutôt qu’en haut. Voyons à quel point :
• Pour une molécule de diazote (masse molaire 28 g/mol) à température ambiante, on
trouve βmg = 1/(8 700 m). À moins que la boîte ne fasse plusieurs kilomètres de haut,
on obtient donc que βmgz ≪ 1 et e−βmgz ≈ 1 pour toutes les valeurs de z. La fonction
fpos est donc, avec une très bonne approximation, également indépendant de z, toutes
les valeurs de z sont équiprobables et on conclut que l’on peut négliger la gravité pour
étudier un gaz dans une boîte.
• Pour un petit grain de sable fin ayant une masse m = 1 μg, à température ambiante, on
trouve βmg = 1/(4 10−13 m). On voit que fpos ∝ e−βmgz est quasiment nul sauf pour
z = 0 ; les grains de sable, même très légers, restent au fond de la boîte ; l’agitation
thermique n’est pas suffisante pour les faire décoller.
83
Ce résultat est expliqué en détail dans une vidéo, et le calcul de l’énergie cinétique de
translation moyenne est également fait en TD. Une conséquence immédiate de ce dernier
résultat est que l’énergie cinétique de translation totale d’un système classique est Ec =
3
2
nRT . Cette relation peut servir de définition de la température dans un système classique.
Dans la section précédente, nous avions obtenu la distribution fimp dans le cas d’une
particule seule dans une boîte ; cette distribution reste valable pour n’importe quel système
classique constitué d’un nombre quelconque de particules, quels que soient le potentiel
d’interaction ou le potentiel extérieur auxquels sont soumises ces particules. En particulier, il
s’applique au mouvement des atomes dans les solides, les liquides ou les gaz. Il ne s’applique
cependant pas au mouvement des électrons dans un métal conducteur, car ces électrons
constituent un système fortement quantique.
⃗p 2
Eσ = + V (x, y, z), (5.17)
2m
avec V le potentiel extérieur agissant sur la particule. À l’équilibre thermodynamique, l’éner-
gie moyenne de la particule est
D ⃗p 2 E 3
⟨E⟩ = + V (x, y, z) = kB T + V (x, y, z) , (5.18)
2m 2
où l’on a utilisé l’encadré 5.3 pour la moyenne de l’énergie cinétique. (Exercice : sauriez-vous
calculer V (x, y, z) dans le cas simple V = mgz pour une particule dans une boîte ?)
Il est remarquable que la moyenne de l’énergie cinétique est indépendante de la masse m
de la particule. En fait, si l’on suppose que l’énergie du système s’écrit sous la forme
Eσ = · · · + bq 2 + · · · , (5.19)
où q est l’une des variables permettant de décrire le système (q peut être une position, une
vitesse, un angle, un moment cinétique, ...), alors on vérifie que l’énergie moyenne s’écrit
1
⟨E⟩ = ⟨· · · ⟩ + kB T + ⟨· · · ⟩. (5.20)
2
C’est le théorème d’équipartition de l’énergie.
Encadré 5.4 : Théorème d’équipartition de l’énergie
Remarques :
⃗p 2 p2x p2y p2z
• 2m = 2m + 2m + 2m . Il y a trois termes quadratiques dans l’énergie cinétique d’une
particule, d’où le résultat de l’encadré 5.3.
84
• Pour un système classique quelconque constitué de N particules, l’énergie cinétique de
translation moyenne est donc ⟨Etranslation ⟩ = 23 N kB T .
• Pour un gaz parfait monoatomique, U = ⟨Etranslation ⟩ = 32 N kB T .
• De manière similaire, on peut retrouver l’énergie interne d’un gaz parfait diatomique
à l’aide du théorème d’équipartition de l’énergie ; voir la vidéo.
85
86
Chapitre 6
La fonction de partition
et donc
Encadré 6.1 : La fonction de partition
où la somme est effectuée sur tous les microétats σ accessibles du système, Eσ est
l’énergie du système dans l’état σ et β = kB1T .
Z dépend de la température T du thermostat (par l’intermédiaire de β) et des
propriétés caractérisant le système. Par exemple, on a Z = Z(T, V, N ) pour un fluide
simple, avec V le volume et N le nombre de particules du système.
87
6.1 Distribution de l’énergie d’un système
Dans un système en contact avec un thermostat, l’énergie fluctue puisque le système
passe d’un microétat à un autre en suivant les probabilités (6.1), et que tous les microétats
ont des énergies différentes.
Cependant, en thermodynamique classique, on considère que l’énergie d’un système ma-
croscopique en contact avec un thermostat est fixée : l’énergie du système est celle qui
minimise ϕ = U − T0 S, voir encadré 3.4.
Comment l’énergie peut-elle fluctuer en physique statistique, et être constante en ther-
modynamique ? La réponse est que les fluctuations existent toujours, mais que pour un grand
système (un système tel qu’on étudie en thermodynamique), ces fluctuations sont extraordi-
nairement petites, et essentiellement impossibles à mesurer.
Nous allons commencer par calculer l’énergie moyenne, puis caractériser ces fluctuations.
∂ h −βEσ i
e −βEσ
Eσ = − e , (6.4)
∂β
puis d’utiliser que la somme d’une dérivée est égale à la dérivée de la somme
où on a reconnu Z (voir encadré 6.1) dans la dernière somme. La dernière expression peut
s’écrire comme la dérivée du logarithme, et on obtient finalement
Encadré 6.2 : Énergie moyenne
Remarque : pour effectuer cette dérivée, on doit écrire Z comme une fonction de β plutôt
que de T . On peut passer d’une dérivée par rapport à β à une dérivée par rapport à T , et
réciproquement ; puisque β = 1/(kB T ) et donc T = 1/(kB β), on a
∂ ∂T ∂ 1 ∂ ∂ ∂ ∂
= =− 2
= −kB T 2 et = −kB β 2 (6.6)
∂β ∂β ∂T kB β ∂T ∂T ∂T ∂β
88
6.1.2 Variance et écart-type
On rappelle que l’amplitude typique des fluctuations d’une variable aléatoire est souvent
caractérisée par son écart-type, qui est égal à la racine carrée de sa variance. En notant σE
l’écart-type et Var(E) la variance de E, on a
D 2 E 2
(6.7)
p
Var(E) = E − ⟨E⟩ = E2 − E , σE = Var(E).
Pour calculer Var(E) en utilisant la deuxième expression de (6.7), il nous reste à calculer
⟨E 2 ⟩ puisque ⟨E⟩ est déjà donné par l’encadré 6.2 ou par (6.5). En utilisant la même astuce
que pour le calcul de ⟨E⟩, on trouve
X 1 X −βEσ 2 1 ∂ 2 X −βEσ 1 ∂ 2Z
⟨E 2 ⟩ = Pσ Eσ2 = e Eσ = e = , (6.8)
σ
Z σ Z ∂β 2 σ Z ∂β 2
∂ 2 ln Z ∂⟨E⟩
Var(E) = 2
=− .
∂β ∂β
89
Encadré 6.4 : Fluctuations de l’énergie
C’est une forme de la loi des grands nombres. En thermodynamique, pour un petit système
de N = 1020 particules, l’écart-type σE est typiquement 1010 fois plus petite que ⟨E⟩. Cela
veut dire que, typiquement, la variable aléatoire E et sa moyenne ⟨E⟩ sont égales sur environ
dix chiffres significatifs, et les fluctuations sont presque impossibles à mesurer.
Pour un système de petite taille (quelques particules), les fluctuations ne sont pas négli-
geables, il ne faut pas confondre E et ⟨E⟩, et la notion d’énergie interne n’a pas vraiment
de sens parce que le système n’est pas thermodynamique.
Revenons au système de grande taille. En identifiant ⟨E⟩ et U et en transformant la
dérivée par rapport à β en dérivée par rapport à T , voir (6.6), on obtient
∂⟨E⟩ ∂U
= −kB T 2 = −kB T 2 C (6.10)
∂β ∂T
Var(E) = kB T 2 C, (6.11)
qui relie les fluctuations de E (membre de gauche) à la capacité thermique C. Par exemple,
pour un gaz parfait monoatomique, on a U = ⟨E⟩ = 23 nRT = 32 N kB T , et donc C = 23 N kB .
On en déduit que
r
3 3
2 2
Var(E) = N kB T , et donc σE = N kB T. (6.12)
2 2
√
On retrouve bien que σE est environ N fois plus petit que ⟨E⟩.
Pour avoir l’énergie E, il faut que le système soit dans un microétat σ tel que Eσ = E. La
probabilité d’un de ces microétats est Pσ = Z1 e−βEσ = Z1 e−βE , mais il ne faut pas oublier
qu’il y a a priori plusieurs microétats ayant l’énergie E. Posons donc :
90
On avait déjà introduit Ω(E) en microcanonique. On a alors
1
P (E) = Ω(E) e−βE . (6.13)
Z
(6.14)
ln P (E) = −β0 E − T0 S(E) − ln Z,
où β0 = 1/(kB T0 ) et où l’on a utilisé la définition de l’entropie, voir encadré 2.4, pour écrire
ln Ω(E) = k1B S(E) = β0 T0 S(E). On voit donc que :
∂S ∂ 2 S ϵ2
S(U + ϵ) = S(U ) + ϵ+ + ··· (6.16)
∂U ∂U 2 2
On utilise
∂ 2S
∂S 1 ∂ 1 1 ∂T 1 ∂U
= , 2
= =− 2 =− avec C = . (6.17)
∂U T ∂U ∂U T T ∂U CT 2 ∂T
ϵ ϵ2
S(U + ϵ) = S(U ) + − + ··· (6.18)
T 2T 2 C
On sait que la température T du système est égale à la température T0 du thermostat et on
trouve alors dans (6.14)
ϵ2
(6.19)
ln P E = U + ϵ = −β U − T0 S(U ) − ln Z − + ···
2kB T02 C
91
Encadré 6.5 : Fluctuations de l’énergie
(E − U )2
P (E) est proportionnel à exp − .
2kB T 2 C
Remarques :
• On sait que pour une variable gaussienne, la variance se lit directement comme la moitié
du dénominateur dans l’exponentielle. On retrouve donc le théorème de fluctuation-
dissipation (6.11) : Var(E) = kB T 2 C, et en particulier on retrouve bien que les fluc-
tuations d’énergie sont négligeables pour un système thermodynamique.
• En mécanique quantique, l’énergie E est une variable à valeurs discrètes, mais on
l’approche par une variable continue pour écrire la distribution gaussienne.
92
on avait ln P (U ) ∼ λN , cela voudrait dire que P (U ) ∼ eλN , et donc que la probabilité de
l’énergie la plus probable serait extraordinairement faible (forcément, on aurait λ < 0). En
fait, on a vu que E pouvait être décrit par une distribution gaussienne de moyenne U et de
variance kB T 2 C. Cela implique (propriété des gaussiennes) que P (U ) = (2πkB T 2 C)−1/2 ∝
N −1/2 (puisque C est extensif). On a donc ln P (U ) ∝ ln N .
L’expression précise de ln P (U ) n’est pas importante ; ce qui est important c’est que
ln P (U ) ≪ F . On en déduit donc que βF = − ln Z et donc
Encadré 6.7 : L’énergie libre
F = −kB T ln Z.
Une manière de se convaincre de cette expression est de vérifier qu’elle est compatible
avec la relation de Gibbs-Helmholtz U = −T 2 ∂(F/T
∂T
)
, qui a été vue en TD. Partant de
l’encadré 6.7, on trouve
∂(F/T ) ∂ ln Z ∂ ln Z
−T 2 = kB T 2 =− (6.21)
∂T ∂T ∂β
où l’on a utilisé (6.6). On reconnaît l’expression de U = ⟨E⟩, encadré 6.2, comme attendu.
93
L’expression (6.25) de l’entropie est en fait générique. Par exemple, en microcanonique,
P 1 par-
tant du postulat de Boltzman Pσ = 1/Ω (voir encadré 2.1), on obtient S = −kB σ Ω ln Ω1 =
kB ln Ω, qui est la définition de l’entropie (voir encadré 2.4). L’expression (6.25) est également
utilisée (sans le préfacteur kB , et avec un logarithme en base 2 plutôt qu’un logarithme népé-
rien) pour définir l’entropie en théorie de l’information (Shannon, 1949). Dans ce contexte,
l’entropie mesure l’information manquante sur un système.
Par exemple, les atomes d’un solide sont discernables : chaque atome vibre autour d’une
position spécifique dont il ne peut s’éloigner. Les particules composant un fluide sont indis-
cernables : chaque particule peut potentiellement se retrouver n’importe où dans le système.
94
(On a écrit que E11 = E1 + E1 , que E12 = E1 + E2 , etc.) Notez que l’état « 12 » est différent
de l’état « 21 », parce qu’on a supposé les particules discernables : il y a la particule de
gauche, et la particule de droite, et on ne peut pas les confondre. On vérifie rapidement que
2
Pour deux particules : Z = e−βE1 + e−βE2 . (6.28)
De même,
Pour trois particules : σ ∈ {111, 112, 121, 122, 211, 212, 221, 222},
(6.29)
Z = e−β(E1 +E1 +E1 ) + e−β(E1 +E1 +E2 ) + e−β(E1 +E2 +E1 ) + . . . + e−β(E2 +E2 +E2 ) .
et, après avoir écrit les huit termes, on vérifie que
3
Pour trois particules : Z = e−βE1 + e−βE2 . (6.30)
Z = zN (6.31)
Autrement dit, pour choisir σ, on doit donner les états θ1 , θ2 , . . ., θN de chacune des parti-
cules. Dans l’exemple du système à deux niveaux, on a θi ∈ {1, 2}.
Avec cette notation, il est clair que sommer sur tous les microétats σ revient à sommer
sur tous les θ1 , et sur tous les θ2 , et sur tous les θ3 , . . ., et sur tous les θN :
X XX X
= ··· . (6.33)
σ θ1 θ2 θN
Comme on suppose que les particules sont sans interactions, l’énergie du système est la
somme des énergies des particules :
E = ϵ1 + ϵ2 + · · · + ϵN , (6.34)
où ϵi est l’énergie de la particule i. Par hypothèse, cette énergie ne dépend que de l’état θi
de la particule i : pour un microétat σ = [θ1 , θ2 , . . . , θN ], une version plus précise de (6.34)
est
Eσ = ϵ1 (θ1 ) + ϵ2 (θ2 ) + · · · + ϵN (θN ), (6.35)
95
où ϵi (θi ) est l’énergie de la particule numéro i quand elle est dans l’état θi . Par exemple,
pour le système à deux niveaux, toutes les particules sont identiques et on a ϵi (1) = E1 et
ϵi (2) = E2 , indépendamment de i.
En utilisant les équations (6.33) et (6.35) dans la définition de la fonction de partition Z
(encadré 6.1), on trouve
X XX X
Z= e−βEσ = ··· e−βϵ1 (θ1 ) e−βϵ2 (θ2 ) · · · e−βϵN (θN ) ,
σ θ θ θ
1X2 N
X X (6.36)
−βϵ1 (θ1 ) −βϵ2 (θ2 ) −βϵN (θN )
= e e ··· e .
θ1 θ2 θN
on obtient
Z = z1 × z2 × · · · × zN . (6.38)
Remarquez que zi ne dépend que de la i-ème particule. D’après sa définition (6.37), zi est
la fonction de partition qu’on aurait obtenue en choisissant pour système la particule i toute
seule. La fonction de partition à N particules discernables sans interactions se factorise donc
en un produit de N fonctions de partition à une particule. Dans l’exemple du système à deux
niveaux, toutes les particules sont identiques et ont la même fonction de partition que l’on
a noté tout simplement z, sans mettre l’indice i.
De même, on vérifie aisément que la loi (6.1) donnant la probabilité d’un microétat σ se
factorise elle aussi :
1 −βEσ 1 1 1
Pσ = e = e−βϵ1 (θ1 ) e−βϵ2 (θ2 ) · · · e−βϵN (θN ) , (6.39)
Z z1 z2 zN
ce qui implique que N particules sans interactions ont des comportements indépendants : la
probabilité d’observer la i-ème particule dans le microétat θi est donnée par
1 −βϵi (θi )
P (θi ) = e , (6.40)
zi
indépendamment des microétats des autres particules. L’énergie moyenne de la particule i
est alors
∂ ln zi
⟨ϵi ⟩ = − , (6.41)
∂β
Un cas particulier important est celui où toutes les particules sont identiques, comme
dans l’exemple du système à deux niveaux ; on a alors z1 = z2 = · · · = zN et l’on note z
(sans indice) la fonction de partition d’une quelconque de ces particules.
96
discernables, c’est-à-dire si on peut effectivement étiqueter toutes les particules : il y a la
particule 1, la particule 2, etc. C’est le cas dans un solide. Si les particules sont indiscernables
(comme dans un fluide), un tel étiquetage n’a pas de sens réel et le microétat de la liste (6.32)
est identique au microétat
Particule 1 : État θ2 ,
Particule 2 : État θ1 ,
Le même microétat σ que (6.32) = (6.42)
...
Particule N : État θN ,
ou à tout autre microétat obtenu par permutation des θi . Autrement dit, le nombre de
microétats n’est pas égal au nombre de listes comme dans le membre de droite de (6.32) ou
(6.42). Dans la somme (6.36), si on suppose que tous les états θi sont différents, alors chaque
microétat est compté N ! fois, c’est-à-dire autant de fois qu’il y a de manières de permuter
les θi . On peut alors corriger ce surcomptage en divisant le résultat par N !.
6.3.4 La factorisation de Z
Encadré 6.9 : Factorisation de Z
Dans un système composé de N particules sans interaction, Z se factorise.
• Si les particules sont différentes, elles sont nécessairement discernables et
Z = z1 × z2 × · · · × zN
avec zi la fonction de partition de la particule i.
• Si les particules sont identiques mais discernables (comme dans un solide),
Z = zN
où z est la fonction de partition d’une particule
• Si les particules sont identiques et indiscernables (comme dans un fluide),
zN
Z≈
N!
Il faut être conscient que le facteur N ! dans l’encadré 6.9 n’est qu’une approximation
basée sur l’hypothèse que toutes les particules sont dans des états différents. Cette approxi-
mation n’est bonne qu’à haute température (à définir), lorsque l’on peut ignorer la mécanique
quantique et traiter le système classiquement. Pour traiter certains problèmes comme le gaz
d’électrons dans un métal ou la condensation de Bose-Einstein, il faudra corriger cette ap-
proximation et faire un dénombrement correct des microétats.
Pour un système thermodynamique constitué de N particules identiques sans interac-
tions, en utilisant les expressions de l’encadré 6.9 et la formule de Stirling, encadré 2.5, les
expressions de U = ⟨E⟩, encadré 6.2, et F , encadré 6.7, deviennent
∂ ln z
U = −N , F = −N kB T ln z (particules discernables),
∂β
(6.43)
∂ ln z h z i
U = −N , F = −N kB T ln + 1 . (particules indiscernables).
∂β N
En particulier, puisque F est extensif, on remarque que z doit être intensif dans le cas
discernable et extensif dans le cas indiscernable.
97
6.3.5 Généralisation de la factorisation
La propriété suivante est une généralisation de la factorisation :
Encadré 6.10 : Factorisation de Z
Soit un système dont l’énergie s’écrit comme une somme de deux termes, et où chaque
terme fait intervenir des variables qui n’apparaissent pas dans l’autre terme,
Le premier terme ne dépend que des positions et vitesses des atomes constituant le réseau
cristallin. Le deuxième terme ne dépend que des positions et vitesses des électrons de conduc-
tion, et le troisième ne dépend que des orientations des spins des différentes particules. En
appliquant le résultat de factorisation, on obtient donc
98
La fonction de partition z pour un oscillateur est
∞ ∞
X 1
X 1 1
z= e−βϵn = e−β 2 ℏω e−βℏωn = e−β 2 ℏω . (6.47)
n=0 n=0
1 − e−βℏω
On considère maintenant une assemblée de N ≫ 1 oscillateurs indépendants discernables.
On a Z = z N et, en utilisant (6.43), on trouve,
1 1 N ℏω
F = N ℏω + N kB T ln 1 − e−βℏω , (6.48)
U = N ℏω + βℏω .
2 2 e −1
La capacité thermique molaire s’obtient par c = 1 ∂U
n ∂T
. On trouve
c
R
1
0,8
0,6
(βℏω)2 eβℏω
c=R 2 0,4 (6.49)
eβℏω − 1 0,2
0
0 0,5 1 1,5 2 2,5 3 kB T
ℏω
La capacité thermique molaire est une fonction du rapport entre l’échelle d’énergie ther-
mique kB T et l’échelle d’énergie microscopique ℏω.
Dans ce système le gap, c’est-à-dire la différence entre le fondamental et le premier état
excité, vaut ℏω. Dans la limite basse température, c’est-à-dire pour kB T ≪ ℏω (ou βℏω ≫ 1),
on trouve que U tend vers l’énergie du fondamental (N ℏω/2) et c tend brutalement vers zéro
à la manière de la loi d’Arrhénius :
U ℏω
≈ + ℏωe−βℏω , c ≈ R(βℏω)2 e−βℏω pour βℏω ≫ 1. (6.50)
N 2
C’est le phénomène décrit dans la section 5.2.3.
À haute température, c’est-à-dire pour kB T ≫ ℏω (ou βℏω ≪ 1), on trouve
U ≈ N kB T = nRT, c≈R pour βℏω ≪ 1. (6.51)
C’est le résultat classique ; en effet, en mécanique classique, un oscillateur harmonique
est tout simplement une masse m accrochée à un ressort de raideur k et de longueur à
vide ℓ0 en une dimension. Son énergie est ϵ = 21 mv 2 + 12 k(ℓ − ℓ0 )2 (avec v la vitesse de la
masse, ℓ la longueur du ressort). Il y a deux termes quadratiques et donc ⟨ϵ⟩ = kB T d’après
le théorème d’équipartition de l’énergie, encadré 5.4. Pour N oscillateurs, on retrouve bien
U = N kB T = nRT et c = R.
99
En mécanique classique dans les cas les plus simples, on doit donner toutes les impulsions
et positions des particules pour décrire un microétat
(6.52)
σ = ⃗r1 ,⃗p1 ,⃗r2 ,⃗p2 , . . . ,⃗rN ,⃗pN ,
et sommer sur tous les microétats revient à « sommer », c’est-à-dire intégrer sur ces variables :
Z Z
1 1
Z= · · · d⃗r1 d⃗p1 d⃗r2 d⃗p2 . . . d⃗rN d⃗pN e−βE(⃗r1 ,⃗p1 ,⃗r2 ,⃗p2 ,...,⃗rN ,⃗pN ) , (6.53)
N ! h3N
où E(. . .) est l’énergie du système dans le microétat correspondant. Le préfacteur 1/h3N est
impossible à deviner dans le seul cadre de la mécanique classique ; nous le justifierons dans
le chapitre suivant. Admettons d’ici là qu’il faut mettre autant de préfacteurs 1/h qu’il faut
pour que Z soit une quantité sans dimension : N particules en trois dimensions font 3N
facteurs « dp » et 3N facteurs « dx », d’où la puissance 3N . Le facteur 1/N ! a été mis entre
crochets parce qu’il n’est présent que pour un système de particules indiscernables (fluide),
pour les mêmes raisons que plus haut. Il faut l’omettre pour des particules discernables
(solide).
Remarquons que si les particules sont sans interactions, c’est-à-dire si
E(⃗r1 ,⃗p1 ,⃗r2 ,⃗p2 , . . . ,⃗rN ,⃗pN ) = ϵ(⃗r1 ,⃗p1 ) + ϵ(⃗r2 ,⃗p2 ) + . . . + ϵ(⃗rN ,⃗pN ), (6.54)
alors on retrouve la propriété de partition se factorisation :
ZZZZZZ
1 1
Z= N
z , avec z = 3 d⃗r d⃗p e−βϵ(⃗r,⃗p) . (6.55)
N! h
6.6 Résumé
Encadré 6.11 : Résumé
La fonction de partition Z apparaît dans la relation de Boltzmann
1 X
Pσ = e−βEσ , Z= e−βEσ .
Z σ
Z dépend de T (qui est fixé) mais pas de E (qui fluctue). Z ne dépend pas de σ.
L’énergie moyenne d’un système en contact avec un thermostat est
∂ ln Z
⟨E⟩ = −
∂β
Dans un petit système, l’énergie fluctue beaucoup.
les fluctuations d’énergie sont négligeables,
Dans un grand système, E = ⟨E⟩ = Ele plus probable = U,
F = −kB T ln Z.
Dans un système de N particules sans interactions, on a
1 N
Z = z N (particules discernables) ou Z ≈ z (particules indiscernables),
N!
où z est la fonction de partition d’une particule seule.
100
Chapitre 7
On considère un gaz formé de molécules identiques (un corps pur, donc). Pour étudier ce
système, nous avons vu que la quantité essentielle à calculer était la fonction de partition Z
définie par l’encadré 6.1. Dans le cas d’un gaz, en première approximation, les interactions
entre molécules peuvent être négligées et l’énergie microscopique totale du gaz est égale
à la somme des énergies individuelles des molécules qui le composent ; c’est le modèle du
gaz parfait. De ce fait, la fonction de partition totale du système Z se factorise (voir la
section 6.3),
zN
Z= , (7.1)
N!
où z est la fonction de partition d’une seule molécule de gaz. On aura alors, voir (6.43),
∂ ln z h z i
U = −N , F = −N kB T ln + 1 . (7.2)
∂β N
Il reste donc à calculer ln z.
Avant de s’intéresser au calcul de ln z, on remarque que l’on peut déjà écrire une expression
pour le potentiel chimique µ : comme z est la fonction de partition à une particule, c’est-à-
dire la fonction de partition obtenue en ne mettant qu’une seule particule dans le système, il
est clair que z ne peut pas dépendre de N . On peut alors dériver F par rapport à n = N/NA
et obtenir, après un petit calcul,
∂F ∂F z z
µ= = NA = NA − kB T ln = −RT ln . (7.3)
∂n V,T ∂N V,T N N
∂U ∂F
µ= = .
! ∂n S,V ∂n T,V
101
On peut également calculer la pression à partir de (7.2). On a déjà remarqué après (6.43)
que, puisque F était extensif, z devait nécessairement l’être également. Or z dépend de V
(extensif) et T (intensif), mais pas de N , et donc la seule écriture extensive possible est
z(V, T ) = V × y(T ) pour une certaine fonction y. On trouve alors
∂F ∂F N kB T
p=− =− = , (7.4)
∂V n,T ∂V N,T V
c’est-à-dire l’équation des gaz parfaits.
7.1 Factorisation de z
On cherche maintenant à calculer z pour obtenir l’énergie interne du gaz, voir (7.2), et
les autres fonctions thermodynamiques. On considère donc une seule molécule du gaz dans
une enceinte de volume V à la température T .
Pour décrire l’état de cette particule, il faut certainement donner sa position et son
impulsion, mais ce n’est en général pas suffisant. Si la molécule est polyatomique, il faut
décrire son mouvement dans le référentiel barycentrique, mouvement qui peut être analysé
en termes de rotation de la molécule et de vibration des distances entre les atomes. On
pourrait aller plus loin et considérer également les excitations du nuage électronique et,
pourquoi pas, du noyau atomique.
En première approximation, toutes ces contributions à l’énergie de la molécule se dé-
couplent et l’on peut écrire l’énergie ϵ de la molécule comme une somme de termes :
L’énergie moyenne est alors la somme des énergies moyennes des différentes contributions :
∂ ln zx
⟨ϵ⟩ = ⟨ϵtranslation ⟩ + ⟨ϵrotation ⟩ + ⟨ϵvibration ⟩ + ⟨ϵélectrons ⟩ + ⟨ϵnoyaux ⟩ + · · · avec ⟨ϵx ⟩ = − ,
∂β
et l’énergie interne est donnée par U = N ⟨ϵ⟩. Enfin, en dérivant U par rapport à T et en
divisant par la quantité de matière n = N/NA , on obtient la capacité thermique molaire cV
à volume constant, qui se décompose en
∂⟨ϵx ⟩
cV = ctranslation + crotation + cvibration + célectrons + cnoyaux + · · · avec cx = NA . (7.7)
∂T
La capacité thermique du gaz est donc la somme des contributions des différents degrés de
liberté de la molécule, et on va maintenant les étudier séparément.
102
7.2 Contribution des termes de translation
Pour calculer ztranslation , on se place dans le cadre le plus général de la mécanique quan-
tique. Les états de translation (sur lesquels on somme) sont les états propres du hamiltonien
de translation,
⃗2
P
Htranslation = , (7.8)
2m
où P⃗ est l’opérateur impulsion et m la masse de la molécule. Les énergies sont les valeurs
propres associées.
h2 2 2 2
(7.9)
ϵtranslation (nx , ny , nz ) = n + n + n z .
8mL2 x y
De manière équivalente, un état propre peut être caractérisé par une impulsion ⃗p dont les
composantes prennent des valeurs discrètes strictement positives 1 , px = 2L
h
nx , py = 2L
h
ny ,
h
pz = 2L nz , de sorte que l’énergie s’écrive
⃗p 2
ϵtranslation (⃗p) = . (7.10)
2m
La fonction de partition de translation est la somme sur tous les états, c’est à dire la
somme sur toutes les impulsions ⃗p possibles :
X
ztranslation = e−βϵtranslation (⃗p) . (7.11)
⃗p autorisés
Nous allons montrer que cette somme peut être remplacée par une intégrale et ainsi
être calculée. Pour cela, commençons par considérer le cas d’une boite à une dimension.
L’impulsion a une seule composante p dont les valeurs possibles sont des multiples entiers
positifs de h/(2L). Si on note
103
Figure 7.1 – La courbe rouge correspond à la fonction f (p). Les rectangles jaunes ont une hauteur
de f (p) aux valeurs autorisées de p et une largeur correspondant à l’écart entre les valeurs de p soit
h/(2L). Le graphe de droite montre la même chose sur un intervalle plus étroit.
correspond donc à l’aire d’un des rectangles sur la figure 7.1. La somme sur tous les rectangles
est alors approximativement l’aire sous la courbe, soit
Z ∞
X h
f (p) × ≃ f (p) dp. (7.13)
p autorisés
2L 0
R∞ R∞
En utilisant la parité de la fonction f , qui implique que 0 f (p) dp = 1
2 −∞
f (p) dp, on
obtient Z ∞
X L
f (p) ≃ f (p) dp. (7.14)
p autorisés −∞ h
où on remarque que L3 est le volume V de l’enceinte. Le résultat obtenu pour une enceinte
cubique est en réalité général et s’applique à n’importe quelle géométrie :
104
Encadré 7.1 : densité d’impulsions pour les états de translation
Pour sommer sur les états décrivant la translation d’une particule dans une boîte de
volume V , on somme sur les impulsions autorisées. Il y a
V d3⃗p
impulsions autorisées dans d3⃗p.
h3
(C’est ce qu’on appelle la densité d’états en impulsion.) Les impulsions autorisées sont
tellement nombreuses que l’on peut remplacer la somme par une intégrale
X X Z Z Z V d3⃗p
... = ... = ... .
états de
h3
⃗p autorisés
translation
V d3⃗p −β ⃗p 2
ZZZ
X p2
⃗
ztranslation = e −β 2m
= e 2m . (7.16)
h3
⃗p autorisés
En mécanique classique, on dirait qu’il faut intégrer sur toutes les impulsions et po-
sitions de la particule. Mais en mécanique quantique, on fait une somme sur les états
! propres d’impulsion ; pour ces états la particule est entièrement délocalisée dans tout
le volume, et il n’y a aucune somme à faire sur le volume.
105
Encadré 7.2 : Longueur d’onde thermique de de Broglie
106
L’oscillateur harmonique a déjà été décrit dans la section 6.4. En reformulant légèrement
les résultats obtenus :
• L’échelle d’énergie caractéristique (qui sépare les domaines basse et haute tempéra-
tures) est
∆vibration = ℏω. (7.20)
C’est la différence d’énergie entre le fondamental et le premier état excité.
• À haute température, on retrouve les résultats classiques
ℏω
⟨ϵvibration ⟩ ≈ , cvibration ≪ R si kB T ≪ ∆vibration (7.22)
2
L’expression exacte de cvibration et sa représentation graphique sont données dans (6.49).
Pour fixer les idées, on a ∆vibration /kB ≈ 3 390 K pour la molécule de diazote. (Cette quan-
tité se mesure par spectroscopie, mais le résultat que l’on trouve est plutôt σ = 2 359 cm−1 .
Il s’agit du nombre d’onde correspondant : σ = λ1 = νc , avec λ la longueur d’onde et ν la
fréquence, et donc ℏω = hν = hcσ.)
À température ambiante T ≈ 300 K, on voit que pour le diazote, kB T ≪ ∆vibration et
cvibration est négligeable ; les modes de vibrations sont gelés à température ambiante. En
utilisant la formule exacte (6.49), on trouve cvibration ≈ 0,0016 R pour le diazote.
−ℓ ≤ m ≤ ℓ. (7.24)
107
Les énergies propres sont donc de la forme
ℓ(ℓ + 1) ℏ2
ϵrotation (ℓ, m) = , ℓ = 0, 1, 2, . . . et m = −ℓ, −ℓ + 1, . . . , ℓ − 1, ℓ. (7.26)
2I
Cette énergie ne dépend pas de m. Par contre, le niveau d’énergie correspondant à une
certaine valeur de ℓ est 2ℓ + 1 fois dégénéré, où 2ℓ + 1 fois est le nombre de valeurs de m
possibles.
Le fondamental a une énergie nulle, le premier état excité une énergie ℏ2 /I. On définit
l’échelle d’énergie caractéristique par
ℏ2
∆rotation = . (7.27)
2I
(On suit la convention habituelle de définir ∆rotation comme la moitié du gap.)
On écrit la fonction de partition en sommant sur tous les états quantiques :
∞ X
X +ℓ ∞
X
zrotation = e −βϵrotation (ℓ,m)
= (2ℓ + 1)e−β∆rotation ℓ(ℓ+1) . (7.28)
ℓ=0 m=−ℓ ℓ=0
108
Calculons un ordre de grandeur. Le crotation
moment d’inertie de deux masses m sé- R
parées par une distance L par rapport 1,2
au milieu des deux masses est I = 2 × 1
m( 2 ) = 2 mL . Si on considère la molécule
L 2 1 2 0,8
de diazote, on a m = (14 g/mol)/NA ≈ 0,6
2 10−26 kg et L ≈ 1 Å = 10−10 m. Cela 0,4
donne I ≈ 10−46 kg m2 . On obtient alors 0,2
kB T
avec (7.27) que ∆rotation /kB ≈ ℏ2 /(2IkB ) ≈ 0
3 K. La valeur exacte est 2,9 K. Cela signi- 0 0,5 1 1,5 2 2,5 ∆rotation
fie qu’à température ambiante, on est dans
la limite haute température pour la rotation, et que crotation ≈ R avec une très grande
précision pour une molécule diatomique.
Pour un gaz usuel à température ambiante, on est dans la limite haute température pour la
rotation et la limite basse température pour la vibration. On dit que la rotation est totalement
excitée, et que la vibration est gelée.
Pour chaque degré de liberté à part la translation, on pourra introduire de la même
manière une échelle d’énergie ∆degré de liberté donnée par le gap telle que cdegré de liberté ≈ 0 si
kB T ≪ ∆degré de liberté . (Pour les degrés de translation, l’expression (7.32) devient fausse à
basse température parce que l’approximation Z = z N /N ! devient incorrecte, et pas parce
kB T devient de l’ordre du gap ∆translation , qui est extrêmement faible.)
On peut estimer l’échelle d’énergie ∆électrons liée à l’excitation du nuage électronique en
considérant le cas simple de l’atome d’hydrogène. On sait d’après la loi de Rydberg que les
niveaux d’énergie de l’électron d’un atome d’hydrogène sont donnés par En = −(13,6 eV)/n2 .
Le gap est donc ∆électron = E2 − E1 = 34 × 13,6 eV = 1,6 10−18 J, ce qui correspond à une
température ∆électron /kB ≈ 120 000 K. On voit que l’électron de l’atome d’hydrogène va
rester coincé dans son fondamental à toute température raisonnable, et il en sera également
de même pour les nuages électroniques des molécules usuelles.
Les noyaux des atomes restent également coincés dans leur fondamental ; les énergies
nécessaires pour les exciter sont de l’ordre du MeV, bien au delà de ce qui est accessible
thermiquement.
109
Quelques valeurs expérimentales (copiées de- rotation vibration électrons
puis l’ouvrage Physique Statistique, L. Couture
et R. Zitoun, ellipses 1992), sont données dans la HD 65,7 K 5 492 K 132 000 K
table ci-contre. On a indiqué, pour chaque degré CO 2,8 K 3 122 K 93 600 K
de liberté, la valeur de ∆/kB , en kelvin. HCl 15,2 K 4 302 K 63 300 K
Encadré 7.3 : Degrés de libertés H2 87,5 K 6 324 K 132 000 K
gelés et excités N2 2,9 K 3 395 K 99 700 K
O2 2,1 K 2 274 K 11 400 K
• Les degrés de translation sont tou-
jours excité.
• Les degrés de rotation sont presque
toujours excités, sauf pour quelques
molécules très légères à basse tempé-
rature.
• Les degrés de vibration sont gelés à
température ambiante, mais peuvent
être excités si on chauffe le gaz.
• On peut ignorer les autres degrés de
liberté (électrons et noyaux), qui sont
toujours gelés.
kB T ℏ2
zrotation ≈ avec ∆rotation = . (7.36)
∆rotation 2I
Voyons si on peut retrouver ce résultat par une approche purement classique.
On considère une molécule constituée de deux atomes de masses m1 et m2 séparés par une
distance ℓ0 fixe, tournant autour de leur centre de gravité G. On sait, par la réduction du problème à
deux corps, que ce système est équivalent à une seule particule virtuelle de masse réduite m = mm11+m
m2
2
à une distance ℓ0 du centre, voir le cours de mécanique du L2. On repère la position de la particule
virtuelle par ses coordonnées sphériques θ et ϕ. La vitesse de cette particule (dans le référentiel du
centre de masse) est ℓ0 (θ̇ ⃗eθ + sin θ ϕ̇⃗eϕ ), et son énergie cinétique de rotation est
1 1
ϵrotation = I(θ̇2 + sin2 θ ϕ̇2 ) = I(ωθ2 + sin2 θ ωϕ2 ), (7.37)
2 2
avec I = mℓ20 le moment d’inertie de la molécule, et ωθ = θ̇ et ωϕ = ϕ̇ les vitesses angulaires.
L’état instantané de la molécule peut donc être donné par les quatre variables θ, ϕ, ωθ et ωϕ .
On voit bien que la fonction de partition doit essentiellement être une intégrale de exp(−βϵrotation )
sur ces quatre variables, mais il y a peut-être des préfacteurs à écrire, qui ne sont pas évident à
deviner.
Voici la règle générale : pour chaque variable que l’on choisit pour décrire un degré de liberté,
on doit intégrer sur cette variable et sur l’impulsion généralisée correspondante (au sens de la
mécanique lagrangienne), puis diviser par h pour chaque paire coordonnée-impulsion. (On rappelle
que l’impulsion généralisée associée à une variable q est p = ∂L ∂ q̇ , où L est le lagrangien.)
C’est en vertu de cette règle que, dans l’approche classique présentée plus haut, on a choisi
comme variables les positions et les impulsions (et non les positions et les vitesses).
Dans le cas de la rotation qui nous intéresse ici, le lagrangien est confondu avec l’énergie cinétique
ϵrotation , et on trouve les impulsions généralisées suivantes :
111
on peut alors réécrire l’énergie cinétique en terme de ces impulsions (c’est l’hamiltonien classique)
p2θ p2ϕ
ϵrotation = + . (7.39)
2I 2I sin2 θ
La fonction de partition classique est alors
π 2π p2 p2
Z Z Z Z
1 ϕ
θ+
zrotation = dθ dϕ dpθ dpϕ e−β 2I 2I sin2 θ (7.40)
h2 0 0
Exercice assez difficile pour les plus motivés En changeant un peu le calcul ci-dessus, on
peut obtenir d’un seul coup les termes de vibration et de rotation dans l’approche classique. Il
suffit de dire que la distance ℓ entre les particules évolue au cours du temps autour de la distance
idéale ℓ0 . On prend en compte l’énergie potentielle harmonique associée 21 mω 2 (ℓ − ℓ0 )2 , et l’énergie
cinétique associée p2ℓ /(2m). On n’oublie pas que I est défini en terme de ℓ0 , et que l’énergie cinétique
de rotation doit donc être corrigée puisque la distance ℓ n’est plus égale à ℓ0 . On intègre donc
finalement sur six variables, et on doit trouver (exercice) :
kB T
zrotation + vibration = zrotation × zvibration × 1 + . (7.41)
Iω 2
Une variante de ce calcul (qui conduit évidemment au même résultat) consiste à décrire la position
de la particule en coordonnées cartésiennes, et d’écrire zrotation + vibration comme une intégrale sur
x, y, z, px , py et pz , comme dans (6.53). Les intégrales sur les impulsions sont immédiates, et les
intégrales restantes sur x, y, z se traitent après changement de variable pour passer en sphérique.
ℏ2
L’expression (7.41) n’est valable (approche classique) que si kB T ≫ ∆rotation = 2I et kB T ≫
∆vibration = ℏω. On voit que pour que l’approximation zrotation + vibration ≈ zrotation × zvibration soit
valable, il faut de plus que kB T ≪ Iω 2 = 12 ∆2vibration /∆rotation .
112
Chapitre 8
puisque les photons sont absorbés et émis par les parois. Cela implique que l’on ne pourra
pas suivre le même programme que pour un gaz parfait de molécules (à savoir calculer z,
la fonction de partition à une seule particule, en déduire Z ≈ z N /N ! (encadré 6.9), puis
calculer E, F , etc.) parce qu’on ne connait pas a priori la valeur de N . On ne peut donc
pas factoriser Z sur les particules mais, comme on va le voir, on pourra factoriser Z sur les
113
états quantiques. Cette factorisation permettra (en Master) de traiter les gaz quantiques et
la condensation de Bose-Einstein.
On a numéroté les états quantiques θ (caractérisés par leur impulsion et leur polarisation)
de manière arbitraire. La liste qui précède est infiniment longue, puisqu’il y a un nombre
114
infini d’états quantiques possibles, mais ce n’est pas très grave. Comme les photons sont des
bosons, il peut y en avoir un nombre arbitraire dans chaque état quantique. Les nombres
nθ peuvent donc prendre toutes les valeurs entières entre zéro et l’infini. Cependant, il n’y a
qu’un nombre fini d’états quantiques θ tels que nθ ̸= 0.
115
Remarque : cette factorisation est très différente de celle de l’encadré 6.9, obtenue pour des
particules massives, qui était Z = z1 × z2 × · · · × zN ou Z = z N /N !. Pour les particules
massives, on avait un produit de N termes, et chaque terme correspondait à une particule.
Dans le cas du gaz de photon, Z est le produit d’un nombre infini de termes, et chaque terme
correspond à un état quantique θ.
e−βnθ ϵθ
Proba(nθ ) = , (8.10)
ζθ
L’énergie moyenne présente dans cet état quantique est alors ⟨nθ ⟩ ϵθ .
On montre facilement à partir de (8.10) (exercice) que l’on a, comme pour (6.41),
∂ ln ζθ
énergie moyenne dans l’état θ = ⟨nθ ⟩ ϵθ = − (8.12)
.
∂β
Insistons encore une fois : l’équation (8.12) donne l’énergie moyenne dans un état θ, où la
moyenne est faite sur le nombre de particules pouvant occuper cet état, alors que l’équa-
tion (6.41) donne l’énergie moyenne portée par une particule, où la moyenne est faite sur
tous les états quantiques dans lesquels on peut trouver cette particule.
Il est temps de calculer ζθ . L’expression de l’encadré 8.1 est une somme géométrique de
raison e−βϵθ et on obtient immédiatement
1
ζθ = . (8.13)
1 − e−βϵθ
L’équation (8.12) permet alors de déterminer ⟨nθ ⟩. On trouve, après simplification,
116
Encadré 8.2 : Le taux d’occupation
On voit donc que le nombre nθ de photons dans un état quantique θ est une fonction
croissante de la température et décroissante de l’énergie ϵθ de chaque photon dans cet état
quantique.
région où l’impulsion a une norme entre p et p + dp. On a donc, après division par V ,
Z ∞
4πp2 dp pc 8πc ∞ p3 dp
Z
U
= 2 = 3 . (8.19)
V 0 h3 eβpc − 1 h 0 eβpc − 1
117
La signification de cette expression est que, dans le four, les photons dont l’impulsion a une
norme comprise entre p et p + dp transportent une énergie par unité de volume égale à
U 8π 5 (kB T )4 π 2 (kB T )4
= = . (8.22)
V 15(hc)3 15(ℏc)3
118
longueur d’onde par unité de temps qu’il en absorberait si cet élément de surface faisait
partie de la paroi d’un four à la température T .
Ce rayonnement est le plus important pour un « corps noir », objet idéal qui absorbe
parfaitement toute lumière, sans aucune réflexion.
Remarque : Savoir quels photons sont absorbés et réfléchis par un morceau de matière, c’est
en connaître la couleur. Un objet qui réfléchit 50% des photons entre 500 nm et 510 nm et
absorbe tout le reste apparaît vert. Cet objet émet spontanément comme un corps noir, sauf
entre 500 nm et 510 nm où l’émission est deux fois plus faible.
cd
considérer que la région d’où peuvent venir les photons
t
est une demi-boule centrée en dS et de rayon c dt. r
On se place en coordonnées sphériques (r, θ, ϕ) avec dS
un axe des z perpendiculaire à dS et rentrant dans le
four. Considérons un petit volume dV = r2 sin θ dr dθ dϕ au point (r, θ, ϕ). La surface dS
est vue depuis ce point avec un angle solide cos θ dS/r2 , et comme toutes les directions sont
équiprobables (isotropie), on obtient que
cos θ dS
un photon donné dans dV a une probabilité d’arriver dans dS. (8.24)
4πr2
Si on pose e l’énergie par unité de volume dans le four, on trouve alors que
On intègre maintenant sur toute la demi-boule, c’est-à-dire pour 0 < r < c dt, 0 < θ < π
2
et
0 < ϕ < 2π et on trouve après un petit calcul
ec
(énergie arrivant sur dS pendant dt) = dt dS, (8.26)
4
et donc la puissance reçue par unité de surface de la paroi du four est ec/4.
119
Encadré 8.4 : Loi du rayonnement du corps noir
2π 5 kB4 π 2 kB4
σ= 3 2
= 3 2
≈ 5,67 10−8 W/m2 /K4 .
15h c 60ℏ c
Remarques : la valeur de σ est très simple à retenir : cinq six sept huit. La quantité σT 4
est parfois appelée l’exitance, ou émittance, du corps. Quand on fait un bilan énergétique
prenant en compte le rayonnement du corps noir, il ne faut pas oublier de comptabiliser
l’énergie électromagnétique reçue venant de l’environnement : Un corps noir de surface S à
la température T dans une enceinte à la température Text émet avec une puissance SσT 4 ,
mais reçoit avec une puissance SσText
4
le rayonnement émis par l’enceinte.
Exitance spectrale (kW/m2 /nm)
visible
infrarouge
ultraviolet
violet
vert
rouge
infrarouge
70 4 10−4
T
=
60 50
0K
50 3 10−4
T
40
=
2 10−4
55
30
00
T= T = 400
K
K
20 45 10−4
00
10 K
T = 3 000 K T = 300 K
0 0
0 500 1 000 1 500 2 000 0 2 000 4 000 6 000 8 000 10 000
λ (nm) λ (nm)
Figure 8.1 – Intensité lumineuse (c’est-à-dire « exitance spectrale ») émise par un corps
noir en fonction de la longueur d’onde pour différentes températures
La longueur d’onde d’un photon, c’est sa couleur, et la puissance émise, c’est l’intensité
lumineuse. Pour déterminer quelle est la puissance émise par le corps noir en fonction de la
longueur d’onde, on utilise encore (8.26), mais en prenant pour e l’expression donnée par
(8.23). On obtient
Encadré 8.5 : Loi de Planck
La « couleur » du rayonnement émis par un corps noir est décrit par la quantité :
Remarque : cette quantité (sans le dλ) est parfois appelée l’exitance (ou émittance) spectrale.
À faible λ, une approximation de cette expression est 2πhc2 e−βhc/λ /λ5 , qui est la loi de
Wien (1896). À grand λ, une approximation est 2πckB T /λ4 , qui est la loi de Rayleigh-Jeans
120
(1900) ; cette dernière expression est le résultat de la mécanique classique (h n’apparaît pas)
et conduit à la « catastrophe ultraviolette », une divergence de la puissance émise à faible
longueur d’onde.
La figure 8.1 donne cette intensité lumineuse en fonction de la longueur d’onde pour
différentes températures.
2,9 106 nm · K
λmax = .
T
En reportant dans l’expression de l’encadré 8.5, on voit que ce maximum d’intensité
varie en T 5 .
Par exemple, pour une température de 5 800 K (la température à la surface du soleil), on
trouve λmax = 500 nm, bien dans le visible. Pour un corps humain, la température de 310 K
donne λmax = 9 400 nm, dans l’infrarouge.
121
122
Appendice A
Compléments mathématiques
∂f f (x + ϵ, y) − f (x, y)
= lim . (A.1)
∂x y ϵ→0 ϵ
On note souvent la (ou les) variables maintenue(s) constante(s) en indice, pour se sou-
venir de quelle fonction on dérive. Cette notation est importante en thermodynamique,
où l’on définit souvent plusieurs fonctions avec le même nom, par exemple U (T, V, n)
ou U (S, V, n). La notation ∂V∂U
T,n
correspond à une dérivée de la première fonction et
vaut 0 pour un gaz parfait. La notation ∂V ∂U
S,n
correspond à une dérivée de la deuxième
fonction et vaut −p pour un gaz parfait. Cependant, si on est certain qu’il n’y a pas
d’ambiguïté, on peut écrire ∂f ∂x
plutôt que ∂f
∂x y
.
• La différentielle de f s’écrit
∂f ∂f
df = dx + dy (A.2)
∂x y ∂y x
La différentielle peut se comprendre comme une manière compacte de noter toutes les
dérivées partielles. Elle correspond également à un accroissement de f : si x augmente
de dx et y de dy, alors f augmente de df .
123
(Dans le membre de droite, il s’agit de la fonction f (u, v). Le fait de noter avec le même
symbole f deux fonctions différentes est une spécificité des physiciens ; un mathématicien ne
ferait jamais ça.)
Écrivons les différentielles des quatre fonctions.
∂u ∂u ∂v ∂v
du = dx + dy, dv = dx + dy (A.4)
∂x y ∂y x ∂x y ∂y x
∂f ∂f ∂f ∂f
df = du + dv, df = dx + dy. (A.5)
∂u v ∂v u ∂x y ∂y x
∂f ∂f ∂u ∂f ∂v
= + (A.7)
∂x y ∂u v ∂x y ∂v u ∂x y
(A.8)
f (x, y) = f u(x, y), x
∂f ∂f ∂u ∂f ∂f ∂f ∂u
= + et = (A.9)
∂x y ∂u x ∂x y ∂x u ∂y x ∂u x ∂y x
Par exemple, si
(A.10)
H(T, p, n) = H S(T, p, n), p, n
alors
∂H ∂H ∂S ∂H
= + . (A.11)
∂p T,n ∂S p,n ∂p T,n ∂p S,n
124
En développant et en comparant avec la différentielle de x(y, z),
∂x ∂x
dx = dy + dz, (A.13)
∂y z ∂z y
on obtient −1
∂x ∂y ∂x ∂x ∂y
= , =− . (A.14)
∂y z ∂x z ∂z y ∂y z ∂z x
En multipliant la deuxième expression par − ∂x
∂z
y
, on trouve
Encadré A.1 : Relations entre les dérivées
−1
∂x ∂y ∂x ∂y ∂z
= , = −1. (A.15)
∂y z ∂x z ∂y z ∂z x ∂x y
∂p ∂T ∂V
× × = −1. (A.16)
∂T V,n ∂V p,n ∂p T,n
et donc
ln(N !) ∼ N ln N − N. (A.20)
Remarques :
• Cette formule est due à James Stirling, mathématicien écossais du XVIIIe siècle, qui
n’a rien à voir avec Robert Stirling, écossais lui aussi, et inventeur du « moteur de
Stirling » au XIXe siècle.)
125
√
• On obtient (A.20) en prenant le logarithme de (A.19) et en négligeant le terme ln 2πN .
Ce√terme est effectivement négligeable : pour N = 1020 , on a N ln N ≈ 46 1020 et
ln 2πN ≈ 24.
• Les approximations (A.19) et (A.20) sont très bonnes :
√
N N! 2πN N N e−N ln(N !) N ln N − N
10 3,629 10 6
3,599 106 15,10 13,03
100 9,333 10157
9,325 10157
363,7 360,5
1000 4,0239 10 2567
4,0235 10 2567
5912 5908
1 000 000 8,263932 105 565 708 8,263931 105 565 708 1,281552 107 1,281551 107
Voici un argument rapide permettant d’obtenir (A.20). (Obtenir (A.19) demanderait un peu
plus de travail...)
ln(N !) = ln N × (N − 1) × (N − 2) × . . . × (N − [N − 1])
= ln N + ln(N − 1) + ln(N − 2) + · · · + ln(N − [N − 1]).
En utilisant ln(N − k) = ln N + ln 1 − Nk , on a
h 1 2 N − 1 i
ln(N !) = N ln N + ln 1 − + ln 1 − + · · · + ln 1 −
Z 1 N N N
≈ N ln N + N dx ln(1 − x)
0
où l’on a reconnu dans la somme entre crochets l’approximation habituelle d’une intégrale
par une somme discrète prise sur un grand nombre de valeurs très rapprochées. Le calcul de
l’intégrale donne le résultat.
126
où l’on a fait le changement de variable x = ⟨X⟩ + u∆X pour obtenir la dernière expression.
Numériquement, on obtient
Proba |X − ⟨X⟩| < ∆X = 68,3% Proba |X − ⟨X⟩| < 2∆X = 95,4%
(A.24)
Proba |X − ⟨X⟩| < 3∆X = 99,7% Proba |X − ⟨X⟩| < 6∆X = 99,9999998%
127