Physics · Livre 5 · Bachelor Year 3

Physique universitaire — 3e année

Physique universitaire — 3e année · Bachelor Year 3

17L’ensemble canonique

Les systèmes isolés sont une fiction de théoricien : les échantillons réels siègent dans des thermostats, des pièces, des océans d’air — en contact avec un réservoir qui fixe non leur énergie mais leur température. Le volume de deuxième année a rencontré empiriquement la loi qui en résulte : la probabilité d’un état d’énergie EE porte le facteur eE/kBT\eu^{-E/k_{\text{B}}T}. Ce chapitre obtient ce facteur de Boltzmann en quatre lignes à partir du comptage du Chapitre 16, puis bâtit la machine qui fait de la physique statistique une discipline industrielle : la fonction de partition ZZ, une somme unique d’où l’énergie, l’entropie, la pression, la capacité thermique et les fluctuations tombent toutes par dérivation. Les premières victoires de la machine sont racontées ici : pourquoi les capacités thermiques meurent à basse température (Einstein, 1907 — la première théorie quantique de la matière), pourquoi la capacité thermique de l’hydrogène gravit un escalier, pourquoi le ciel s’amincit exponentiellement — et comment Perrin, en regardant des grains microscopiques sédimenter dans une goutte d’eau, a compté le nombre d’Avogadro et clos le débat sur la réalité des atomes.

17.1 La distribution de Boltzmann, obtenue

Théorème 17.1 (Distribution canonique)

Un système qui échange de l’énergie avec un grand réservoir à la température TT occupe son micro-état ss (d’énergie EsE_s) avec la probabilité

Ps=eEs/kBTZ,Z=seEs/kBT\mathcal P_s = \frac{\eu^{-E_s/k_{\text{B}}T}}{Z} , \qquad Z = \sum_s \eu^{-E_s/k_{\text{B}}T}

(β=1/kBT\beta = 1/k_{\text{B}}T dans toute la suite). ZZ, la fonction de partition, normalise les probabilités — et se trouve contenir toute la thermodynamique du système.

Démonstration. Le système dans l’état ss laisse au réservoir l’énergie EEsE - E_s : par le postulat fondamental appliqué à l’ensemble, PsΩres(EEs)\mathcal P_s \propto \Omega_{\text{res}}(E - E_s). Développons l’entropie du réservoir, énorme et régulière, au premier ordre : Sres(EEs)=Sres(E)EsSres/E=Sres(E)Es/TS_{\text{res}}(E - E_s) = S_{\text{res}}(E) - E_s\,\partial S_{\text{res}}/\partial E = S_{\text{res}}(E) - E_s/T par la définition de la température. Donc ΩreseEs/kBT\Omega_{\text{res}} \propto \eu^{-E_s/k_{\text{B}}T}. (Les ordres supérieurs meurent avec la taille du réservoir.)

Proposition 17.2 (La machine)

À partir de Z(T,V,N)Z(T, V, N) :

E=lnZβ,F=kBTlnZ,S=FT,P=FV,\langle E\rangle = -\frac{\partial\ln Z}{\partial\beta} , \qquad F = -k_{\text{B}}T\ln Z , \qquad S = -\frac{\partial F}{\partial T} , \qquad P = -\frac{\partial F}{\partial V} ,

F=ETSF = \langle E\rangle - TS est l’énergie libre ; et la capacité thermique est une fluctuation :

C=ET=E2E2kBT2.C = \frac{\partial\langle E\rangle}{\partial T} = \frac{\langle E^2\rangle - \langle E\rangle^2} {k_{\text{B}}T^2} .

Pour des sous-systèmes indépendants et discernables Z=zNZ = z^N ; pour NN particules identiques dans une même boîte, Z=zN/N!Z = z^N/N!. À TT et VV fixés, l’équilibre minimise FF : la nature troque énergie contre entropie au taux de change TT — le principe le plus utile de toute la physique de la matière.

Démonstration partielle. βlnZ=EseβEs/Z=E-\partial_\beta\ln Z = \sum E_s\eu^{-\beta E_s}/Z = \langle E\rangle ; une dérivée de plus donne E2E2\langle E^2\rangle - \langle E\rangle^2, et la règle de dérivation composée convertit β\partial_\beta en T\partial_T. Les identifications de FF, SS, PP suivent de la comparaison de  ⁣d(kBTlnZ)\dd(-k_{\text{B}}T\ln Z) avec le  ⁣dF=S ⁣dTP ⁣dV\dd F = -S\dd T - P\dd V thermodynamique ; la factorisation est l’exponentielle d’une somme. Le principe de minimum : que FsystF_{\text{syst}} décroisse, c’est que StotalS_{\text{total}} croisse, puisque ΔSres=ΔEsys/T\Delta S_{\text{res}} = -\Delta E_{\text{sys}}/T.

Le cadre canonique : un petit système qui emprunte de l’énergie à un vaste réservoir. Chaque joule emprunté coûte au réservoir 1/T d’entropie — d’où la remise exponentielle sur les états énergétiques.
Le cadre canonique : un petit système qui emprunte de l’énergie à un vaste réservoir. Chaque joule emprunté coûte au réservoir 1/T1/T d’entropie — d’où la remise exponentielle sur les états énergétiques.

17.2 Premières victoires

Exemple 17.3 (Deux niveaux, en deux lignes)

Pour les niveaux 0,ϵ0, \epsilon : z=1+eβϵz = 1 + \eu^{-\beta\epsilon}, d’où E=Nϵ/(eβϵ+1)\langle E\rangle = N\epsilon/(\eu^{\beta\epsilon} + 1) — le résultat qui coûtait une page de Stirling par la voie microcanonique (Exercice 16.8). Le formalisme canonique est le microcanonique avec la combinatoire prédigérée ; la bosse de Schottky dans C(T)C(T) suit d’une seule dérivation.

Exemple 17.4 (Le solide d’Einstein, et la mort de Dulong–Petit)

Modélisons un cristal par 3N3N oscillateurs quantiques de pulsation ω\omega (Chapitre 9). Par oscillateur,

z=neβω(n+1/2)=12sinh(βω/2),E=ω2+ωeβω1.z = \sum_n\eu^{-\beta\hbar\omega(n + 1/2)} = \frac{1}{2\sinh(\beta\hbar\omega/2)} , \qquad \langle E\rangle = \frac{\hbar\omega}{2} + \frac{\hbar\omega}{\eu^{\beta\hbar\omega} - 1} .

À haute TT : EkBT\langle E\rangle \to k_{\text{B}}T par oscillateur, C3NkBC \to 3Nk_{\text{B}} — la loi de Dulong–Petit du volume de première année, expliquée. À basse TT le quantum ω\hbar\omega devient inabordable et CC s’effondre exponentiellement — comme mesuré, et inexplicable classiquement. La courbe d’Einstein de 1907, un paramètre par élément, fut la première application des quanta à la matière ordinaire ; le diamant, aux liaisons raides et aux atomes légers (ω/kB1300K\hbar\omega/k_{\text{B}} \approx 1300\,\mathrm{K}), est encore « gelé » à température ambiante — l’anomalie qui intriguait les chimistes depuis quatre-vingts ans (Exercice 17.6).

À gauche : la courbe de capacité thermique d’Einstein — l’équipartition classique retrouvée à haute T, le gel quantique au-dessous de _ E = /k_ B. À droite : le C_V du dihydrogène gazeux gravit un escalier à mesure que la rotation (vers 85\, K) puis la vibration (vers 6000\, K, hors échelle) dégèlent : chaque mouvement ne rejoint l’équipartition que lorsque k_ BT peut payer son quantum.
À gauche : la courbe de capacité thermique d’Einstein — l’équipartition classique retrouvée à haute TT, le gel quantique au-dessous de θE=ω/kB\theta_{\text{E}} = \hbar\omega/k_{\text{B}}. À droite : le CVC_V du dihydrogène gazeux gravit un escalier à mesure que la rotation (vers 85K85\,\mathrm{K}) puis la vibration (vers 6000K6000\,\mathrm{K}, hors échelle) dégèlent : chaque mouvement ne rejoint l’équipartition que lorsque kBTk_{\text{B}}T peut payer son quantum.

Théorème 17.5 (L’équipartition, et sa licence)

Toute coordonnée ou quantité de mouvement qui entre quadratiquement dans l’énergie contribue, dans le régime classique (haute température), pour

ϵ=12kBT\langle\epsilon\rangle = \tfrac12 k_{\text{B}}T

à l’énergie moyenne : 32kBT\tfrac32 k_{\text{B}}T pour un atome de gaz monoatomique, 52\tfrac52 pour une diatomique qui tourne, 3kBT3k_{\text{B}}T pour un oscillateur. La licence expire lorsque kBTk_{\text{B}}T tombe sous l’écart de niveaux du mode : le mode gèle et sa contribution disparaît — la résolution des scandales du dix-neuvième siècle sur les capacités thermiques, dessinée en escalier ci-dessus.

Démonstration partielle. Pour ϵ=ax2\epsilon = ax^2,

ϵ=ax2eβax2 ⁣dxeβax2 ⁣dx=βln ⁣eβax2 ⁣dx=βlnβ1/2=12β.\langle\epsilon\rangle = \frac{\int ax^2\,\eu^{-\beta ax^2}\dd x} {\int\eu^{-\beta ax^2}\dd x} = -\partial_\beta\ln\!\int\eu^{-\beta ax^2}\dd x = -\partial_\beta\ln\beta^{-1/2} = \frac{1}{2\beta} .

Le gel est le calcul de Exemple 17.4, mode par mode.

Exemple 17.6 (Le gaz, canoniquement)

Un atome dans une boîte : z=V/λT3z = V/\lambda_T^3 (le comptage d’états de Proposition 7.5 pondéré par Boltzmann) ; NN atomes identiques : Z=zN/N!Z = z^N/N!. Alors F=NkBT[ln(V/NλT3)+1]F = -Nk_{\text{B}}T[\ln(V/N\lambda_T^3) + 1], et la dérivation livre PV=NkBTPV = Nk_{\text{B}}T, l’entropie de Sackur–Tetrode et E=32NkBT\langle E\rangle = \tfrac32 Nk_{\text{B}}T — tout le gaz parfait à partir d’une intégrale gaussienne. Le facteur de Boltzmann appliqué à l’énergie cinétique d’une molécule donne la distribution des vitesses de Maxwell du volume de première année, désormais obtenue ; appliqué à son énergie potentielle mghmgh, il donne l’atmosphère exponentielle — et, dans une goutte d’eau, l’échelle de grains de Perrin (Problème 17.1).

Méthode 17.7 (L’art du canonique)

(1) Dresser la liste des états et des énergies d’une unité ; calculer zz. (2) Factoriser : Z=zNZ = z^N (ou zN/N!z^N/N! pour des particules identiques qui partagent l’espace) ; prendre le ln\ln tôt. (3) Dériver : en β\beta pour l’énergie, en TT pour l’entropie via FF, en VV pour la pression ; une dérivée seconde pour CC et les fluctuations. (4) Vérifier les deux bouts : la haute TT doit reproduire l’équipartition, la basse TT doit geler avec une queue en eΔ/kBT\eu^{-\Delta/k_{\text{B}}T}. (5) Compétitions (repliement, fixation, alignement) : écrire F=ETSF = E - TS pour chaque possibilité et laisser gagner la plus petite — la transition se situe vers TΔE/ΔST \approx \Delta E/\Delta S.

17.3 Exercices

Exercice 17.1

(a) Écrire le rapport de populations de deux niveaux séparés de ϵ\epsilon à la température TT. (b) Dans une flamme à 2500K2500\,\mathrm{K}, quelle fraction des atomes de sodium siège dans l’état excité à 2.1eV2.1\,\mathrm{eV} de la raie D (dégénérescence 2 pour le fondamental, 6 pour l’excité : rapport de populations 3eβϵ3\eu^{-\beta\epsilon}) ? (c) Pourquoi la flamme flamboie-t-elle malgré tout en jaune (combien d’atomes par cm3^3 suffisent) ? (d) À quelle température la fraction excitée atteindrait-elle 10%10\% ?

Solution

Solution de Exercice 17.1.

(a) P2/P1=(g2/g1)eβϵ\mathcal P_2/\mathcal P_1 = (g_2/g_1)\eu^{-\beta\epsilon}. (b) βϵ=2.1/(8.62×105×2500)=9.7\beta\epsilon = 2.1/(8.62\times10^{-5} \times 2500) = 9.7 : fraction 3e9.72×1043\eu^{-9.7} \approx 2 \times 10^{-4}. (c) Une flamme porte 1015\sim10^{15} atomes de sodium par cm3^3 : même 10410^{-4} d’entre eux, recyclés toutes les quelques nanosecondes, déversent 101910^{19} photons jaunes par seconde — de quoi éblouir. (d) 3eβϵ=0.13\eu^{-\beta\epsilon} = 0.1 : T=ϵ/(kBln30)7200KT = \epsilon/(k_{\text{B}}\ln30) \approx 7200\,\mathrm{K} — une photosphère stellaire, pas une flamme.

Exercice 17.2

Un système à trois niveaux : 0,ϵ,2ϵ0, \epsilon, 2\epsilon. (a) Écrire zz. (b) Calculer E\langle E\rangle et vérifier les deux limites de température. (c) À quelle TT le niveau intermédiaire est-il le plus peuplé en valeur absolue ? (d) Montrer qu’aucune température, si haute soit-elle, ne rend un niveau supérieur plus peuplé qu’un niveau inférieur — quelle notion du chapitre 16 faudrait-il pour cela ?

Solution

Solution de Exercice 17.2.

(a) z=1+eβϵ+e2βϵz = 1 + \eu^{-\beta\epsilon} + \eu^{-2\beta\epsilon}. (b) E=ϵ(eβϵ+2e2βϵ)/z\langle E\rangle = \epsilon(\eu^{-\beta\epsilon} + 2\eu^{-2\beta\epsilon})/z : 0\to 0 à basse TT, ϵ\to \epsilon (le niveau moyen) à haute. (c) P1=1/(eβϵ+1+eβϵ)\mathcal P_1 = 1/(\eu^{\beta \epsilon} + 1 + \eu^{-\beta\epsilon}) croît de façon monotone avec TT, n’approchant sa borne 1/31/3 que pour TT \to \infty. (d) Les poids de Boltzmann ne font que décroître avec l’énergie à T>0T > 0 ; une inversion de population exige les températures négatives de Exemple 16.8, inaccessibles à tout réservoir.

Exercice 17.3

L’atmosphère isotherme. (a) Appliquer le facteur de Boltzmann à l’énergie potentielle mghmgh et obtenir n(h)=n0emgh/kBTn(h) = n_0\eu^{-mgh/k_{\text{B}} T}. (b) Calculer la hauteur caractéristique pour l’air (m=4.8×1026kgm = 4.8 \times 10^{-26}\,\mathrm{kg}) à 288K288\,\mathrm{K}. (c) La pression au sommet de l’Everest en fraction de celle du niveau de la mer. (d) Pourquoi la décroissance de l’atmosphère réelle est-elle proche d’une exponentielle sans l’être exactement (qu’avons-nous tenu pour constant qui ne l’est pas) ?

Solution

Solution de Exercice 17.3.

(a) Le facteur de Boltzmann sur Ep=mghE_p = mgh à TT uniforme. (b) h0=kBT/mg=8.4kmh_0 = k_{\text{B}}T/mg = 8.4\,\mathrm{km}. (c) e8848/84400.35\eu^{-8848/8440} \approx 0.35 : un tiers d’atmosphère — d’où l’oxygène des sommiteux. (d) L’atmosphère réelle n’est pas isotherme : la température baisse avec l’altitude, si bien que le profil s’écarte d’une exponentielle unique.

Exercice 17.4

Comptabilité d’équipartition. Compter les termes quadratiques et prédiser le CVC_V molaire (a) de l’argon ; (b) de N2_2 à température ambiante (rotation active, vibration gelée) ; (c) de N2_2 à 3000K3000\,\mathrm{K} ; (d) d’un solide cristallin (Dulong–Petit). Où les valeurs mesurées 12.512.5, 20.820.8, 26\approx26, 25J/(molK)\approx25\,\mathrm{J}/(\mathrm{mol}\,\mathrm{K}) s’accordent-elles, et qu’enseigne chaque écart ?

Solution

Solution de Exercice 17.4.

(a) Trois translations : 32R=12.5J/(molK)\tfrac32 R = 12.5\,\mathrm{J}/(\mathrm{mol}\,\mathrm{K}), accord exact. (b) On ajoute deux rotations : 52R=20.8J/(molK)\tfrac52 R = 20.8\,\mathrm{J}/(\mathrm{mol}\,\mathrm{K}), comme mesuré : la vibration est gelée. (c) 72R=29.1J/(molK)\tfrac72 R = 29.1\,\mathrm{J}/(\mathrm{mol}\,\mathrm{K}) prédits ; les 26\approx26 mesurés montrent une vibration à peine dégelée (θvib3400K\theta_{\text{vib}} \approx 3400\,\mathrm{K}). (d) 3R=24.9J/(molK)3R = 24.9\,\mathrm{J}/(\mathrm{mol}\,\mathrm{K}) : Dulong–Petit, respectée à température ambiante par la plupart des métaux et violée par le diamant — l’histoire du gel de Exercice 17.6.

Exercice 17.5 ★★

L’oscillateur quantique, canoniquement. (a) Sommer la série géométrique pour zz. (b) Obtenir E\langle E\rangle et identifier n=1/(eβω1)\langle n\rangle = 1/(\eu^{\beta\hbar\omega} - 1) — le résultat de Exercice 9.6, désormais sans effort. (c) Dériver pour C(T)C(T) et vérifier les deux limites. (d) Pourquoi la forme de n\langle n\rangle est-elle sur le point de devenir célèbre (Chapitre 20) ?

Solution

Solution de Exercice 17.5.

(a) z=eβω/2/(1eβω)z = \eu^{-\beta\hbar\omega/2}/(1 - \eu^{-\beta\hbar\omega}). (b) E=βlnz\langle E\rangle = -\partial_\beta\ln z donne la forme annoncée avec n=1/(eβω1)\langle n\rangle = 1/(\eu^{\beta\hbar\omega} - 1). (c) CkBC \to k_{\text{B}} à haute TT ; CkB(βω)2eβω0C \approx k_{\text{B}}(\beta \hbar\omega)^2\eu^{-\beta\hbar\omega} \to 0 à basse. (d) Avec ω\hbar\omega pour énergie d’un quantum de lumière, n\langle n\rangle est le nombre thermique de photons par mode : la loi de Planck est à un chapitre d’ici.

Exercice 17.6 ★★

Einstein face aux données. (a) D’après la formule de la figure, à quel T/θET/\theta_{\text{E}} CC est-elle tombée à la moitié de Dulong–Petit ? (b) Le diamant : θE1300K\theta_{\text{E}} \approx 1300\,\mathrm{K} — calculer C/3NkBC/3Nk_{\text{B}} à 300K300\,\mathrm{K} et expliquer l’« anomalie du diamant » du dix-neuvième siècle. (c) Le plomb : θE90K\theta_{\text{E}} \approx 90\,\mathrm{K} — pourquoi le plomb a-t-il toujours été « bien élevé » ? (d) Les mesures à très basse TT montrent CT3C \propto T^3, non une exponentielle : laquelle des hypothèses d’Einstein tombe (tous les oscillateurs à une seule fréquence), et qui l’a réparée (Chapitre 20) ?

Solution

Solution de Exercice 17.6.

(a) Numériquement, C=12×3NkBC = \tfrac12 \times 3Nk_{\text{B}} vers T/θE0.34T/\theta_{\text{E}} \approx 0.34. (b) θE/T=4.3\theta_{\text{E}}/T = 4.3 : C/3NkB=4.32e4.3/(1e4.3)20.26C/3Nk_{\text{B}} = 4.3^2\eu^{-4.3}/(1 - \eu^{-4.3})^2 \approx 0.26 — le diamant à température ambiante n’a guère qu’un quart de la capacité thermique classique : l’« anomalie » est un gel quantique en pleine lumière. (c) θE=90K\theta_{\text{E}} = 90\,\mathrm{K} place le plomb en plein régime classique à 300K300\,\mathrm{K}. (d) L’hypothèse d’une fréquence unique : les solides réels ont un spectre de modes qui descend jusqu’au son de grande longueur d’onde, dont les quanta bon marché donnent la loi en T3T^3 — la réparation de Debye, au Chapitre 20.

Exercice 17.7 ★★

La queue de Maxwell et l’hydrogène manquant. (a) À partir du facteur de Boltzmann sur l’énergie cinétique, écrire la distribution des vitesses et situer la vitesse la plus probable pour N2_2 et H2_2 à 288K288\,\mathrm{K}. (b) La vitesse de libération terrestre est 11.2km/s11.2\,\mathrm{km}/\mathrm{s} : calculer mvesc2/2kBTmv_{\text{esc}}^2/2k_{\text{B}}T pour les deux gaz. (c) La fraction qui s’échappe varie comme emvesc2/2kBT\eu^{-mv_{\text{esc}}^2/2k_{\text{B}}T} : comparer les deux exposants et conclure quel gaz fuit à l’échelle géologique. (d) Relier à la composition observée de l’atmosphère terrestre (sans H2_2 libre) et de celle de Jupiter (surtout H2_2) — quels deux paramètres renversent le verdict ?

Solution

Solution de Exercice 17.7.

(a) f(v)v2emv2/2kBTf(v) \propto v^2\eu^{-mv^2/2k_{\text{B}}T} : vp=2kBT/mv_{\text{p}} = \sqrt{2k_{\text{B}}T/m} : 413m/s413\,\mathrm{m}/\mathrm{s} pour N2_2, 1540m/s1540\,\mathrm{m}/\mathrm{s} pour H2_2. (b) mvesc2/2kBT730mv_{\text{esc}}^2/2k_{\text{B}}T \approx 730 pour N2_2, 5252 pour H2_2. (c) e730\eu^{-730}, c’est jamais ; e521023\eu^{-52} \sim 10^{-23} par temps de séjour est lent — mais la haute atmosphère chaude (1000K\sim1000\,\mathrm{K}) adoucit l’exposant de l’hydrogène jusqu’à 15\sim15 : l’hydrogène s’échappe à l’échelle géologique, l’azote reste. (d) La vitesse de libération et la température de l’exosphère : le puits à 60km/s60\,\mathrm{km}/\mathrm{s} de Jupiter rend astronomique même l’exposant de l’hydrogène — les géantes gazeuses gardent ce que les petits mondes tièdes perdent.

Exercice 17.8 ★★

Les fluctuations rencontrent la réponse. (a) Démontrer que C=(E2E2)/kBT2C = (\langle E^2\rangle - \langle E\rangle^2)/k_{\text{B}}T^2 à partir de deux dérivées de lnZ\ln Z. (b) Le vérifier explicitement sur le système à deux niveaux. (c) Pour NN unités indépendantes, montrer que la fluctuation relative d’énergie décroît comme 1/N1/\sqrt N. (d) Énoncer la morale : ce que la disposition d’un système à absorber de la chaleur (une réponse) a à voir avec le tremblement de son énergie (une fluctuation) — le marché récurrent de la physique statistique.

Solution

Solution de Exercice 17.8.

(a) β2lnZ=E2E2\partial_\beta^2\ln Z = \langle E^2\rangle - \langle E\rangle^2, et C=TE=kBβ2βEC = \partial_T\langle E\rangle = -k_{\text{B}}\beta^2\partial_\beta\langle E\rangle. (b) Les deux membres donnent NkB(βϵ)2eβϵ/(eβϵ+1)2Nk_{\text{B}}(\beta\epsilon)^2\eu^{\beta\epsilon}/( \eu^{\beta\epsilon} + 1)^2. (c) EN\langle E\rangle \propto N, ΔEN\Delta E \propto \sqrt N. (d) L’intensité avec laquelle un système répond au chauffage égale l’ampleur du tremblement spontané de son énergie — réponse et fluctuation sont deux lectures de la même dérivée seconde, un motif (fluctuation–dissipation) qui revient partout en physique.

Exercice 17.9 ★★

Boltzmann au laboratoire de chimie. Les vitesses de réaction portent le facteur eEa/kBT\eu^{-E_a/k_{\text{B}}T} (franchir une barrière d’activation EaE_a — Arrhenius). (a) Montrer que la règle empirique « la vitesse double tous les 10K10\,\mathrm{K} près de la température ambiante » correspond à Ea0.55eVE_a \approx 0.55\,\mathrm{eV}. (b) De quel facteur cette réaction ralentit-elle dans un réfrigérateur (5C5\,{}^{\circ}\mathrm{C}) ? (c) Cuire un œuf en altitude : l’eau bout à 93C93\,{}^{\circ}\mathrm{C} sur un col à 2000m2000\,\mathrm{m} — estimer le temps de cuisson supplémentaire. (d) Pourquoi la queue au-delà d’une barrière, et non l’énergie moyenne, gouverne-t-elle la chimie (quelles molécules réagissent) ?

Solution

Solution de Exercice 17.9.

(a) ln2=EaΔT/kBT2\ln2 = E_a\,\Delta T/k_{\text{B}}T^2 avec ΔT=10K\Delta T = 10\,\mathrm{K}, T=298KT = 298\,\mathrm{K} : Ea=0.693kBT2/100.55eVE_a = 0.693\,k_{\text{B}}T^2/10 \approx 0.55\,\mathrm{eV}. (b) De 298298 à 278K278\,\mathrm{K} : facteur eEa(1/2781/298)/kB4.4\eu^{E_a(1/278 - 1/298)/k_{\text{B}}} \approx 4.4 de ralentissement — d’où la conservation des aliments au réfrigérateur. (c) 373366K373 \to 366\,\mathrm{K} : la vitesse chute de 1.4\approx1.4 : l’œuf de onze minutes en demande un quart d’heure. (d) Les réactions se gagnent sur la queue exponentielle des molécules au-dessus de la barrière : décaler un peu la température et la population de cette queue bouge énormément — la moyenne n’y est presque pour rien.

Exercice 17.10 ★★★

Repliement à deux états. Une biomolécule est repliée (énergie 00, une seule configuration) ou dépliée (énergie ΔE>0\Delta E > 0, Ωu=eΔS/kB\Omega_u = \eu^{\Delta S/k_{\text{B}}} configurations). (a) Écrire la fraction repliée en fonction de TT. (b) Montrer que le milieu de la « fusion » est Tm=ΔE/ΔST_{\text{m}} = \Delta E/\Delta S et l’interpréter comme un partage F=ETSF = E - TS. (c) Avec ΔE=3.0eV\Delta E = 3.0\,\mathrm{eV} et ΔS=100kB\Delta S = 100\,k_{\text{B}} (l’unité coopérative d’une petite protéine) : calculer TmT_{\text{m}} et la largeur de la transition. (d) Pourquoi la coopérativité (de nombreux contacts qui lâchent ensemble, ΔE\Delta E et ΔS\Delta S grands) affine-t-elle la fusion — et comment les machines de cyclage thermique de l’ADN (PCR) en tirent-elles exactement parti ?

Solution

Solution de Exercice 17.10.

(a) freplieˊ=1/(1+eβ(ΔETΔS))f_{\text{replié}} = 1/(1 + \eu^{-\beta(\Delta E - T\Delta S)}), l’entropie de l’état déplié étant absorbée dans son énergie libre. (b) En Tm=ΔE/ΔST_{\text{m}} = \Delta E/\Delta S les deux énergies libres s’égalent : moitié-moitié. (c) Tm=4.8×1019/1.38×1021=348KT_{\text{m}} = 4.8 \times 10^{-19}/1.38 \times 10^{-21} = 348\,\mathrm{K} (75C75\,{}^{\circ}\mathrm{C}) ; largeur δTkBTm2/ΔE3.5K\delta T \sim k_{\text{B}}T_{\text{m}}^2/\Delta E \approx 3.5\,\mathrm{K}. (d) La coopérativité multiplie ΔE\Delta E et ΔS\Delta S par le nombre de contacts qui lâchent ensemble, ce qui laisse TmT_{\text{m}} mais rétrécit la largeur 1/ΔE\propto 1/\Delta E : les brins d’ADN se séparent sur deux kelvins, et c’est ce qui permet à une machine PCR de cycler proprement entre « fondu » et « apparié ».

Exercice 17.11 ★★★

Paramagnétisme et refroidissement magnétique. NN spins 12\tfrac12 de moment μ\mu dans un champ BB. (a) Montrer que M=Nμtanh(μB/kBT)M = N\mu\tanh(\mu B/k_{\text{B}}T) et développer en la loi de Curie MNμ2B/kBTM \approx N\mu^2B/k_{\text{B}}T. (b) Évaluer l’alignement μB/kBT\mu B/k_{\text{B}}T pour des moments électroniques à B=1TB = 1\,\mathrm{T}, T=300KT = 300\,\mathrm{K}, puis à 1K1\,\mathrm{K}. (c) Désaimantation adiabatique : aimanter à 1K1\,\mathrm{K}, isoler, réduire BB d’un facteur dix — montrer qu’une entropie constante signifie un μB/kBT\mu B/k_{\text{B}}T constant, donc que TT chute d’un facteur dix. (d) Qu’est-ce qui fixe le plancher de ce réfrigérateur (les interactions entre spins — estimer l’échelle dipolaire μ0μ2/4πa3\mu_0\mu^2/4\pi a^3 pour a=0.5nma = 0.5\,\mathrm{nm}, en unités de température) ?

Solution

Solution de Exercice 17.11.

(a) z=2cosh(βμB)z = 2\cosh(\beta\mu B) ; M=Nμtanh(βμB)Nμ2B/kBTM = N\mu\tanh(\beta\mu B) \approx N\mu^2B/k_{\text{B}}T à petit argument : le 1/T1/T de Curie. (b) μBB/kBT\mu_{\text{B}}B/k_{\text{B}}T : 2.2×1032.2 \times 10^{-3} à 300K300\,\mathrm{K} ; 0.670.67 à 1K1\,\mathrm{K} — de l’indifférence à l’alignement marqué. (c) SS ne dépend que de μB/kBT\mu B/k_{\text{B}}T ; baisser BB à entropie fixée entraîne TT proportionnellement : 1K0.1K1\,\mathrm{K} \to 0.1\,\mathrm{K}. (d) Quand kBTk_{\text{B}}T atteint l’énergie spin–spin, l’entropie n’est plus commandée par le champ : μ0μB2/4πa37×1026J5mK\mu_0\mu_{\text{B}}^2/4\pi a^3 \approx 7 \times 10^{-26}\,\mathrm{J} \sim 5\,\mathrm{mK} — le plancher classique (les moments nucléaires, mille fois plus faibles, le repoussent au microkelvin).

Exercice 17.12 ★★★

L’escalier de rotation, quantitativement. Les niveaux de rotation d’une diatomique sont EJ=BJ(J+1)E_J = BJ(J+1), de dégénérescence 2J+12J + 1 (Chapitre 10). (a) Écrire zrotz_{\text{rot}} et montrer que pour kBTBk_{\text{B}}T \gg B la somme devient l’intégrale kBT/Bk_{\text{B}}T/B : le kBk_{\text{B}} d’équipartition retrouvé dans CVC_V. (b) Définir θrot=B/kB\theta_{\text{rot}} = B/k_{\text{B}} et l’évaluer pour H2_2 (B=7.5meVB = 7.5\,\mathrm{meV}) et N2_2 (B=0.25meVB = 0.25\,\mathrm{meV}) : quel gaz montre la marche de rotation à des températures accessibles ? (c) Représenter tout l’escalier CV(T)C_V(T) de l’hydrogène avec ses trois paliers et ses deux contremarches, en chiffrant les deux. (d) Le comportement mesuré à basse TT de H2_2 est encore compliqué par le mélange ortho–para 3:1 de Exercice 14.8 : dites en une phrase comment la statistique des spins nucléaires atteint la capacité thermique d’un gaz.

Solution

Solution de Exercice 17.12.

(a) zrot=J(2J+1)eβBJ(J+1)(2J+1)eβBJ(J+1) ⁣dJ=kBT/Bz_{\text{rot}} = \sum_J(2J+1)\eu^{-\beta BJ(J+1)} \to \int(2J+1)\eu^{-\beta BJ(J+1)}\dd J = k_{\text{B}}T/B : alors E=kBT\langle E\rangle = k_{\text{B}}T et Crot=kBC_{\text{rot}} = k_{\text{B}}. (b) H2_2 : θrot=87K\theta_{\text{rot}} = 87\,\mathrm{K} — la marche se situe dans la plage du laboratoire ; N2_2 : 2.9K2.9\,\mathrm{K}, gelée seulement au voisinage de l’hélium liquide, si bien que l’azote montre toujours 52R\tfrac52 R. (c) Paliers 32R\tfrac32 R (au-dessous de 50K\sim50\,\mathrm{K}), 52R\tfrac52 R (de 200\sim200 à 1000K\sim1000\,\mathrm{K}), puis montée vers 72R\tfrac72 R près de θvib6000K\theta_{\text{vib}} \approx 6000\,\mathrm{K}. (d) Les JJ pairs et impairs appartiennent à des espèces de spin nucléaire différentes qui s’interconvertissent lentement, si bien que le CVC_V de l’hydrogène froid dépend de son histoire ortho–para — la statistique nucléaire auditée au calorimètre.

17.4 Problème : compter Avogadro dans une goutte d’eau

Problème 17.1

Problème du week-end — les grains de Perrin et la réalité des atomes

En 1908, Jean Perrin mit en suspension dans l’eau des grains de résine microscopiques, les laissa sédimenter et les compta couche par couche au microscope. L’échelle exponentielle qu’il trouva était le facteur de Boltzmann rendu visible — et de sa hauteur caractéristique il tira le nombre d’Avogadro, convainquant les derniers sceptiques que les atomes existent. Prix Nobel, 1926. Données : grains de gomme-gutte de rayon r=0.212µmr = 0.212\,\text{µ}\mathrm{m}, masse volumique ρg=1207kg/m3\rho_{\text{g}} = 1207\,\mathrm{kg}/\mathrm{m}^{3} ; eau ρw=999kg/m3\rho_{\text{w}} = 999\,\mathrm{kg}/\mathrm{m}^{3} ; T=293KT = 293\,\mathrm{K} ; g=9.81m/s2g = 9.81\,\mathrm{m}/\mathrm{s}^{2}.

Partie I — Le facteur de Boltzmann visible.

  1. Un grain dans l’eau subit la pesanteur diminuée de la poussée d’Archimède : calculer sa masse effective m=43πr3(ρgρw)m' = \tfrac43\pi r^3(\rho_{\text{g}} - \rho_{\text{w}}) et son poids.
  2. Écrire le profil de concentration d’équilibre n(h)n(h) à partir du facteur de Boltzmann.
  3. Calculer la hauteur caractéristique h0=kBT/mgh_0 = k_{\text{B}}T/m'g avec le kBk_{\text{B}} moderne.
  4. Pourquoi les grains doivent-ils être si précisément calibrés (comment h0h_0 dépend-elle de rr) ?
  5. Comparer h0h_0 à la même formule pour les molécules d’air : pourquoi l’atmosphère des grains fait-elle quelques micromètres quand celle de l’air fait des kilomètres ?
  6. Perrin compta (dans une série) les concentrations relatives 100:55:30:17100 : 55 : 30 : 17 à quatre profondeurs équidistantes de 30µm30\,\text{µ}\mathrm{m} : vérifier que c’est bien une échelle exponentielle et en extraire h0h_0.

Partie II — Peser l’invisible.

  1. Inverser : à partir du h0h_0 mesuré à la question 6 et du mgm'g connu, extraire kBk_{\text{B}}.
  2. La constante des gaz R=8.314J/(molK)R = 8.314\,\mathrm{J}/(\mathrm{mol}\,\mathrm{K}) était connue de la chimie macroscopique : la combiner à votre kBk_{\text{B}} pour obtenir le nombre d’Avogadro NA=R/kBN_{\text{A}} = R/k_{\text{B}}.
  3. Les séries de Perrin donnèrent NAN_{\text{A}} entre 5.55.5 et 7.2×10237.2 \times 10^{23} : comparer à la valeur moderne 6.022×10236.022 \times 10^{23} et commenter l’exploit, vu son microscope et son chronomètre.
  4. À partir de NAN_{\text{A}}, calculer la masse d’un unique atome d’hydrogène — le nombre que les atomistes voulaient depuis un siècle.
  5. Expliquer la structure logique : quelles mesures macroscopiques (RR ; taille et masse volumique des grains ; l’échelle) se combinent pour peser un atome, sans qu’aucun atome n’ait jamais été vu ?
  6. Pourquoi les grains — 101010^{10} masses atomiques chacun — obéissent-ils à la même statistique de Boltzmann que les molécules (qu’est-ce qui, dans la démonstration de Théorème 17.1, se soucie de la taille) ?

Partie III — Le tremblement qui scelle l’affaire.

  1. Ces mêmes grains tremblent : la formule d’Einstein de 1905 pour le mouvement brownien donne x2=2Dt\langle x^2\rangle = 2Dt avec D=kBT/6πηrD = k_{\text{B}}T/6\pi\eta r (eau : η=103Pas\eta = 10^{-3}\,\mathrm{Pa}\,\mathrm{s}). Calculer DD pour les grains de Perrin.
  2. Jusqu’où un grain vagabonde-t-il en une minute ? Perrin pouvait-il le mesurer avec un oculaire micrométrique et un chronomètre ?
  3. Perrin vérifia x2t\langle x^2\rangle \propto t et en extrait kBk_{\text{B}} de nouveau, indépendamment : pourquoi deux voies sans rapport (une échelle statique ; un tremblement dynamique) vers un même nombre ont-elles porté un tel poids de preuve ?
  4. Le tremblement est l’équipartition appliquée au grain : chaque composante de vitesse porte 12kBT\tfrac12 k_{\text{B}}T. Estimer la vitesse thermique quadratique moyenne du grain (m4.8×1017kgm \approx 4.8 \times 10^{-17}\,\mathrm{kg}).
  5. Pourquoi cette vitesse n’est-elle jamais vue directement (qu’est-ce qui interrompt le vol libre après quelques nanomètres), et que voit-on à la place ?
  6. Dites quels deux chapitres de ce livre se rencontrent dans cette observation : la mécanique de la traînée (volume de deuxième année, viscosité) et la statistique de ce chapitre.

Partie IV — Ce qui fut réglé.

  1. Ostwald et Mach tenaient les atomes pour des fictions comptables : dites en une phrase pourquoi un NAN_{\text{A}} compté à partir d’échelles de grains a mis fin à cette position.
  2. Citer trois autres voies des années 1900 qui ont convergé vers le même NAN_{\text{A}} (le bleu du ciel, Exercice 16.12 ; l’électrolyse jointe à la charge de l’électron ; la production d’hélium par radioactivité) — pourquoi la convergence a-t-elle compté plus qu’aucune valeur isolée ?
  3. L’échelle de Perrin est un équilibre entre lesquelles deux monnaies de ce chapitre (l’énergie qui tire vers le bas, l’entropie qui étale vers le haut), tarifé à quel taux ?
  4. Usages modernes de la même physique : les ultracentrifugeuses analytiques font tourner des protéines à 105g10^5g pour comprimer leurs « atmosphères » en échelles mesurables — montrer que multiplier gg par 10510^5 divise h0h_0 d’autant, et estimer h0h_0 pour une protéine de masse effective 1022kg10^{-22}\,\mathrm{kg} à 105g10^5g, 293K293\,\mathrm{K}.
  5. Doubler le rayon du grain divise h0h_0 par huit : encadrer la fenêtre pratique de tailles de grains entre « échelle trop haute pour voir un gradient » et « échelle plus mince qu’un grain » pour un champ de microscope profond de 100µm100\,\text{µ}\mathrm{m}.
  6. Depuis 2019 le SI définit exactement kBk_{\text{B}} et NAN_{\text{A}} : dites ce que mesure aujourd’hui une expérience de Perrin exacte (un contrôle de cohérence, ou un étalonnage de l’appareil et des grains) — et pourquoi la physique est inchangée.
  7. Résumer le résultat nommé : la concentration d’un grain de 0.2µm0.2\,\text{µ}\mathrm{m} est divisée par deux tous les 35µm\sim35\,\text{µ}\mathrm{m} de hauteur ; lue avec n(h)=n0emgh/kBTn(h) = n_0\eu^{-m'gh/k_{\text{B}}T}, cette échelle a livré NA6×1023N_{\text{A}} \approx 6 \times 10^{23} — des atomes comptés, non conjecturés, dans une goutte d’eau sur une platine de microscope.
Solution

Solution de Problème 17.1.

1. V=43πr3=4.0×1020m3V = \tfrac43\pi r^3 = 4.0 \times 10^{-20}\,\mathrm{m}^{3} : m=VΔρ=8.3×1018kgm' = V \Delta\rho = 8.3 \times 10^{-18}\,\mathrm{kg}, poids mg=8.1×1017Nm'g = 8.1 \times 10^{-17}\,\mathrm{N}. 2. n(h)=n0emgh/kBTn(h) = n_0\,\eu^{-m'gh/k_{\text{B}}T} : la loi barométrique, réduite à une lame de microscope. 3. h0=kBT/mg=4.04×1021/8.1×101750µmh_0 = k_{\text{B}}T/m'g = 4.04 \times 10^{-21}/8.1 \times 10^{-17} \approx 50\,\text{µ}\mathrm{m}. 4. h01/r3h_0 \propto 1/r^3 : une dispersion de 10%10\% sur le rayon fait 30%30\% sur la hauteur caractéristique — des grains polydisperses brouillent l’échelle en bouillie. Perrin les a fractionnés pendant des mois par centrifugations répétées. 5. Même formule, des masses distantes de 101010^{10} : l’« atmosphère » des grains est 101010^{10} fois moins haute — les kilomètres se réduisent à quelques dizaines de micromètres, ce qui est précisément ce qui la rend observable en entier au microscope. 6. Rapports successifs 0.550.55, 0.550.55, 0.570.57 : constants à l’erreur de comptage près — exponentiel. h0=30µm/ln(100/55)50µmh_0 = 30\,\text{µ}\mathrm{m}/\ln(100/55) \approx 50\,\text{µ}\mathrm{m}. 7. kB=mgh0/T=8.1×1017×5.0×105/2931.4×1023J/Kk_{\text{B}} = m'g\,h_0/T = 8.1 \times 10^{-17} \times 5.0 \times 10^{-5}/293 \approx 1.4 \times 10^{-23}\,\mathrm{J}/\mathrm{K}. 8. NA=R/kB6.0×1023mol1N_{\text{A}} = R/k_{\text{B}} \approx 6.0 \times 10^{23}\,\mathrm{mol}^{-1}. 9. À quelques pour cent près ici (avec des données idéalisées) ; les séries réelles de Perrin se dispersaient de ±15%\pm15\% autour de la valeur moderne — stupéfiant pour des grains comptés à la main, et absolument décisif pour l’ordre de grandeur. 10. mH=103kg/mol/NA=1.7×1027kgm_{\text{H}} = 10^{-3}\,\mathrm{kg}/\mathrm{mol}/N_{\text{A}} = 1.7 \times 10^{-27}\,\mathrm{kg}. 11. La chimie macroscopique fournit RR ; la microscopie optique et la pesée fournissent mm' ; le comptage fournit l’échelle : trois mesures de paillasse triangulent la masse d’un atome que personne ne peut voir. 12. Rien : la démonstration n’a utilisé que « système qui échange de l’énergie avec un réservoir » — le facteur de Boltzmann est aveugle à la taille, et c’est exactement ce que Perrin a vérifié. 13. D=kBT/6πηr=4.04×1021/4.0×1091012m2/sD = k_{\text{B}}T/6\pi\eta r = 4.04 \times 10^{-21}/4.0 \times 10^{-9} \approx 10^{-12}\,\mathrm{m}^{2}/\mathrm{s}. 14. 2Dt=1.2×101011µm\sqrt{2Dt} = \sqrt{1.2 \times 10^{-10}} \approx 11\,\text{µ}\mathrm{m} par minute : confortablement mesurable avec un réticule d’oculaire et de la patience. 15. Deux phénomènes indépendants, deux formules indépendantes, un seul nombre : l’accord entre l’échelle statique et le tremblement dynamique ne laissait aucune niche au hasard — l’hypothèse moléculaire prédisait les deux. 16. kBT/m=4.04×1021/4.8×10179mm/s\sqrt{k_{\text{B}}T/m} = \sqrt{4.04 \times 10^{-21}/ 4.8 \times 10^{-17}} \approx 9\,\mathrm{mm}/\mathrm{s} par composante. 17. Le grain est heurté 101910^{19} fois par seconde et oublie sa vitesse en quelques nanomètres : le vol balistique est inobservable, et ce que l’œil voit est son intégrale — la marche aléatoire diffusive. 18. La traînée de Stokes (la viscosité des fluides du volume de deuxième année) fournit le 6πηr6\pi\eta r ; l’ensemble canonique fournit le kBTk_{\text{B}}T : le DD d’Einstein est leur quotient, la mécanique et la statistique dans une seule fraction. 19. Une fiction ne se compte pas : dès lors que NAN_{\text{A}} est le rapport de deux nombres mesurés, avec des barres d’erreur, les atomes deviennent des objets d’expérience — Ostwald a capitulé par écrit en 1909. 20. La diffusion qui bleuit le ciel, l’électrolyse jointe à la charge mesurée de l’électron, l’hélium accumulé à partir du radium : quatre voies physiques sans rapport convergeant vers un même 6×10236\times10^{23} ont fait de ce nombre une propriété de la nature plutôt que d’une théorie. 21. L’énergie gravitationnelle qui tire les grains vers le bas, l’entropie configurationnelle qui les étale vers le haut, échangées au taux TT : l’échelle est le minimum de F=ETSF = E - TS. 22. h01/gh_0 \propto 1/g : à 105g10^5g, h0=kBT/(m×105g)=4.04×1021/9.8×101740µmh_0 = k_{\text{B}}T/(m'\times10^5g) = 4.04 \times 10^{-21}/ 9.8 \times 10^{-17} \approx 40\,\text{µ}\mathrm{m} pour la protéine — l’équilibre de sédimentation, l’expérience de Perrin refaite chaque jour dans les départements de biochimie. 23. De « trop haute » (h0h_0 \gg profondeur du champ : aucun gradient visible) à « trop mince » (h0rh_0 \lesssim r) : les rayons utilisables s’étendent grossièrement de 0.10.1 à 0.5µm0.5\,\text{µ}\mathrm{m} — le choix de Perrin n’était pas de la chance mais de la conception. 24. kBk_{\text{B}} étant désormais exact par définition, la même expérience étalonne les grains (leur taille ou leur masse volumique) ou audite le montage : la physique — l’échelle de Boltzmann — est intacte ; seule a bougé la grandeur qui compte comme inconnue. 25. La population d’un grain de 0.2µm0.2\,\text{µ}\mathrm{m} est divisée par deux tous les 35µm\sim35\,\text{µ}\mathrm{m} ; lue à travers n0emgh/kBTn_0\eu^{-m'gh/k_{\text{B}} T}, l’échelle a rendu NA6×1023N_{\text{A}} \approx 6 \times 10^{23} : le nombre d’Avogadro compté grain par grain — et le débat atomique clos sur une platine de microscope.