À la fin du XIXe siècle, on cherche à comprendre la distribution spectrale de l’énergie émise par un corps noir, c’est-à-dire un système en équilibre thermique qui absorbe et réémet parfaitement le rayonnement électromagnétique. Expérimentalement, on mesure une densité spectrale d’énergie $ u(\nu, T) $, c’est-à-dire l’énergie par unité de volume et par unité de fréquence.
Les expériences réalisées avec des cavités chauffées montrent un fait remarquable : lorsque la température augmente, le rayonnement émis ne devient pas seulement plus intense, il se déplace vers des longueurs d’onde de plus en plus courtes. À basse température, l’émission est principalement dans l’infrarouge ; lorsque la température augmente, la cavité rougit, puis jaunit, puis blanchit. Autrement dit, le maximum de la distribution spectrale se déplace vers les hautes fréquences. Le problème théorique consiste donc à expliquer quantitativement cette redistribution de l’énergie en fonction de la température, et à retrouver la forme exacte de la courbe expérimentale.
La tentative classique pour établir cette loi repose sur une modélisation très précise : on enferme le rayonnement dans une cavité métallique maintenue à la température $ T $. Les parois imposent des conditions aux champs électromagnétiques. Le problème devient alors celui des modes propres du champ électromagnétique dans une boîte.
Considérons une cavité cubique de côté $ L $. Les champs électromagnétiques doivent satisfaire les équations de Maxwell avec des conditions aux limites imposées par les parois conductrices : le champ électrique tangentiel s’annule sur les parois.
Ces conditions quantifient les longueurs d’onde possibles. Les composantes du vecteur d’onde $ \vec{k} $ ne peuvent prendre que des valeurs discrètes :
\[ k_x = \frac{n_x \pi}{L}, \quad k_y = \frac{n_y \pi}{L}, \quad k_z = \frac{n_z \pi}{L} \]où $ n_x, n_y, n_z $ sont des entiers positifs.
Chaque triplet $ (n_x, n_y, n_z) $ définit un mode propre du champ. La fréquence correspondante est donnée par la relation de dispersion du vide :
\[ \omega = c |\vec{k}|, \qquad \nu = \frac{c}{2\pi} |\vec{k}| \]La présence des facteurs $ \pi $ dans les composantes du vecteur d’onde n’est pas arbitraire : elle provient directement des conditions aux limites imposées par les parois conductrices de la cavité.
Considérons, pour simplifier, une dimension spatiale. Le champ électrique doit s’annuler sur les parois situées en $ x = 0 $ et $ x = L $. On cherche donc des solutions stationnaires de la forme
\[ E(x) \propto \sin(k_x x) \]La condition $ E(0)=0 $ est automatiquement satisfaite. La condition $ E(L)=0 $ impose
\[ \sin(k_x L) = 0 \]Or le sinus s’annule lorsque son argument vaut un multiple entier de $ \pi $. On doit donc avoir
\[ k_x L = n_x \pi \]d’où
\[ k_x = \frac{n_x \pi}{L} \]Le facteur $ \pi $ apparaît donc parce que les fonctions propres compatibles avec deux nœuds fixes aux extrémités sont des sinus dont la demi-longueur d’onde correspond à la distance entre deux zéros successifs, c’est-à-dire $ \pi/k_x $. Autrement dit, la cavité impose que sa longueur contienne un nombre entier de demi-longueurs d’onde.
En trois dimensions, le raisonnement est identique dans chaque direction.
Les solutions se factorisent en
\[ \sin\!\left(\frac{n_x \pi x}{L}\right) \sin\!\left(\frac{n_y \pi y}{L}\right) \sin\!\left(\frac{n_z \pi z}{L}\right) \]ce qui conduit naturellement à
\[ k_x = \frac{n_x \pi}{L}, \quad k_y = \frac{n_y \pi}{L}, \quad k_z = \frac{n_z \pi}{L} \]Les facteurs $ \pi $ traduisent donc simplement la structure mathématique des solutions d’onde stationnaire dans une boîte à parois imposant des nœuds.
Enfin, lorsqu’on passe de la pulsation $ \omega $ à la fréquence ordinaire $ \nu $, un autre facteur $ 2\pi $ apparaît, mais pour une raison différente : il provient de la relation générale
\[ \omega = 2\pi \nu \]Le $ \pi $ issu des conditions aux limites et le $ 2\pi $ reliant fréquence angulaire et fréquence ordinaire ont donc des origines conceptuellement distinctes, bien qu’ils apparaissent dans les mêmes formules.
Chaque triplet $ (n_x, n_y, n_z) $ définit un mode propre du champ. La fréquence correspondante est donnée par la relation de dispersion du vide :
\[ \omega = c |\vec{k}|, \quad \text{ou} \quad \nu = \frac{c}{2\pi} |\vec{k}| \]Ainsi, le problème du rayonnement thermique devient un problème de comptage des points autorisés dans l’espace des vecteurs d’onde, appelé espace $ k $.
Dans l’espace $ k $, les valeurs autorisées forment un réseau cubique. Chaque point correspond à un mode électromagnétique.
Le nombre de modes dont le vecteur d’onde a une norme comprise entre $ k $ et $ k + dk $ correspond au nombre de points du réseau situés dans une coquille sphérique de rayon $ k $ et d’épaisseur $ dk $.
Le volume d’une coquille sphérique en espace $ k $ vaut :
\[ 4\pi k^2 dk \]Chaque mode occupe un volume élémentaire $ (\pi/L)^3 $ dans cet espace. En effet, rappelons que chaque point du réseau est séparé de ses voisins par une distance \( \pi/L \) selon chaque axe.Le nombre total de modes dans la coquille est donc proportionnel à :
\[ \frac{4\pi k^2 dk}{(\pi/L)^3} \]En tenant compte des deux polarisations indépendantes du champ électromagnétique, on obtient, après simplification et passage à la variable fréquence $ \nu $, une densité de modes proportionnelle à :
\[ \nu^2 d\nu \]C’est ici qu’apparaît le facteur $ \nu^2 $, qui joue un rôle central dans la loi de Rayleigh–Jeans.
Lorsque l’on compte les modes du champ électromagnétique dans la cavité, chaque point du réseau dans l’espace $ k $ correspond à une direction de propagation donnée, c’est-à-dire à un vecteur d’onde $ \vec{k} $.
Mais cela ne suffit pas pour décrire complètement l’onde.
Une onde électromagnétique n’est pas un simple scalaire : c’est un champ vectoriel. Pour une direction de propagation donnée, le champ électrique peut osciller selon différentes directions perpendiculaires à $ \vec{k} $.
Dans le vide, l’une des équations de Maxwell est
\[ \nabla \cdot \vec{E} = \frac{\rho}{\varepsilon_0} \] où $ \rho $ est la densité de charge.Or, dans le vide sans charges libres, $ \rho = 0 $. On obtient donc
\[ \nabla \cdot \vec{E} = 0 \]On cherche maintenant une solution de type onde plane :
\[ \vec{E}(\vec{r},t) = \vec{E}_0\, e^{i(\vec{k}\cdot\vec{r}-\omega t)} \]Ici, $ \vec{E}_0 $ est un vecteur constant (l’amplitude), et toute la dépendance spatiale est contenue dans le facteur exponentiel.
Calculons maintenant la divergence :
\[ \nabla \cdot \vec{E} = \nabla \cdot \left( \vec{E}_0 e^{i(\vec{k}\cdot\vec{r}-\omega t)} \right) \]Le vecteur $ \vec{E}_0 $ est constant, donc le gradient n’agit que sur l’exponentielle.
Or, une dérivée spatiale d’une exponentielle de ce type fait apparaître un facteur $ i\vec{k} $. Plus précisément,
\[ \nabla e^{i\vec{k}\cdot\vec{r}} = i\vec{k}\, e^{i\vec{k}\cdot\vec{r}} \]Donc,
\[ \nabla \cdot \vec{E} = i (\vec{k} \cdot \vec{E}_0) \, e^{i(\vec{k}\cdot\vec{r}-\omega t)} \]Mais Maxwell impose que
\[ \nabla \cdot \vec{E} = 0 \]L’exponentielle n’est jamais nulle. La seule façon pour que l’expression soit nulle est donc que
\[ \vec{k} \cdot \vec{E} = 0 \]Cette équation signifie simplement que le champ électrique est perpendiculaire au vecteur d’onde.
Si l’onde se propage dans la direction de $ \vec{k} $, alors $ \vec{E} $ ne peut pas osciller dans cette direction.
C’est exactement ce qu’on appelle une onde transverse.
Par exemple, si l’onde se propage selon $ z $, alors $ \vec{E} $ ne peut avoir que des composantes selon $ x $ et $ y $, jamais selon $ z $.
C’est cette contrainte géométrique qui explique pourquoi il existe exactement deux polarisations indépendantes pour chaque direction de propagation.
Dans ce plan perpendiculaire à $ \vec{k} $, il existe deux directions indépendantes possibles.
C’est cela qu’on appelle les deux polarisations indépendantes.
Pour un vecteur d’onde donné $ \vec{k} $, il existe un plan perpendiculaire à ce vecteur.
Dans ce plan, on peut choisir deux vecteurs unitaires orthogonaux entre eux :
\[ \vec{e}_1 \quad \text{et} \quad \vec{e}_2 \]Toute polarisation possible est une combinaison linéaire de ces deux directions.
Autrement dit, pour chaque vecteur $ \vec{k} $, il existe deux modes physiques distincts.
Chaque polarisation correspond à un oscillateur harmonique indépendant.
Quand on calcule le nombre de modes dont la fréquence est comprise entre $ \nu $ et $ \nu + d\nu $, le comptage géométrique dans l’espace $ k $ donne le nombre de vecteurs d’onde possibles.
Mais pour chaque vecteur $ \vec{k} $, il faut multiplier par 2 pour tenir compte des deux polarisations.
C’est précisément ce facteur 2 qui conduit, après calcul complet, à la densité de modes :
\[ \frac{8\pi \nu^2}{c^3} \]au lieu d’un facteur $ 4\pi \nu^2 / c^3 $.
En termes modernes, on dirait que le photon possède deux états de polarisation (ou deux degrés de liberté transverses).
Dans le cadre purement classique, cela signifie simplement que pour chaque fréquence autorisée, il existe deux oscillateurs indépendants qui peuvent stocker de l’énergie.
C’est pourquoi la densité d’énergie finale est multipliée par 2.
La physique classique considère maintenant que chaque mode propre du champ est équivalent à un oscillateur harmonique.
Pourquoi ? Parce que les équations de Maxwell, pour chaque mode normal, se réduisent à une équation différentielle du type :
\[ \ddot{q} + \omega^2 q = 0 \]qui est précisément l’équation d’un oscillateur harmonique de pulsation $ \omega $. L’énergie associée à chaque mode prend alors la forme :
\[ E = \frac{1}{2} (p^2 + \omega^2 q^2) \]avec $ q $ et $ p $ jouant le rôle de variables canoniques.
Ainsi, le rayonnement dans la cavité est vu comme un gaz d’oscillateurs harmoniques indépendants.
C’est chaque mode du champ électromagnétique qui est mathématiquement équivalent à un oscillateur harmonique.
Un oscillateur harmonique classique, c’est :
\[ H=\frac{p^2}{2m}+\frac{1}{2}m\omega^2 q^2 \]avec une masse $m$, une position $q$, une impulsion $p$.
Un photon :
Donc un photon n’est pas un oscillateur matériel.
Dans une cavité, le champ électromagnétique peut être décomposé en modes normaux.
Chaque mode se comporte dynamiquement comme une variable $q(t)$ qui satisfait :
\[ \ddot q + \omega^2 q = 0. \]C’est exactement l’équation d’un oscillateur harmonique.
Mais ici :Il n’y a pas de masse matérielle.
Il y a un degré de liberté dynamique.
Le hamiltonien d’un mode du champ s’écrit formellement :
\[ H=\frac{1}{2}(p^2+\omega^2 q^2). \]Cela ressemble à un oscillateur de masse unité.
Ce n’est qu’une analogie mathématique.
Le « $p$ » n’est pas l’impulsion d’une particule massive.
C’est la variable canonique conjuguée à l’amplitude du champ.
Quand on quantifie cet oscillateur, son énergie devient :
\[ E_n = \left(n+\frac12\right)h\nu. \]Chaque quantum d’excitation $ n \to n+1 $ correspond à l’ajout d’un photon.
Donc :
Un photon n’est pas un oscillateur.
C’est une excitation quantifiée d’un oscillateur de champ.
Oui.
Un photon a :
\[ E=h\nu, \qquad p=\frac{E}{c}=\frac{h\nu}{c}. \]Même sans masse au repos, il possède une quantité de mouvement relativiste.
La relation générale relativiste est :
\[ E^2 = p^2 c^2 + m^2 c^4. \]Pour le photon $m=0$, donc :
\[ E=pc. \]Il a donc une quantité de mouvement, mais pas de masse au repos.
| Objet | Ce que c’est |
|---|---|
| Oscillateur harmonique | Modèle mathématique du mode du champ |
| Mode électromagnétique | Degré de liberté du champ dans la cavité |
| Photon | Quantum d’énergie d’un mode |
Donc quand on dit : “Un photon est un oscillateur harmonique”, c’est incorrect.
La formulation correcte est : chaque mode du champ est un oscillateur harmonique. Les photons sont les quanta de ces oscillateurs.
La théorie quantique des champs dit :
La lumière n’est pas faite de “petites billes” au départ. Elle est faite d’un champ électromagnétique quantifié.
Les photons sont les unités minimales d’excitation de ce champ.
La mécanique statistique classique fournit alors un résultat fondamental : le théorème de l’équipartition.
Ce théorème affirme que, à l’équilibre thermique, chaque degré de liberté quadratique dans l’énergie contribue en moyenne pour :
\[ \frac{1}{2} k_B T \]Or, un oscillateur harmonique possède deux termes quadratiques indépendants dans son énergie : un terme cinétique en $ p^2 $ et un terme potentiel en $ q^2 $. Par conséquent, l’énergie moyenne d’un oscillateur harmonique classique vaut :
\[ \langle E \rangle = k_B T \]Ce résultat est indépendant de la pulsation $ \omega $. Tous les modes, qu’ils soient de basse ou de haute fréquence, reçoivent en moyenne la même énergie.
On peut maintenant combiner les deux ingrédients essentiels :
La densité de modes entre $ \nu $ et $ \nu + d\nu $ est proportionnelle à $ \nu^2 d\nu $.
Chaque mode transporte en moyenne une énergie $ k_B T $.
La densité spectrale d’énergie devient donc :
\[ u(\nu, T) \propto \nu^2 k_B T \] En calculant précisément les constantes numériques, on obtient la loi de Rayleigh–Jeans : \[ u(\nu, T) = \frac{8\pi \nu^2}{c^3} k_B T \] Cette loi est parfaitement cohérente avec les principes de la physique classique.Cependant, cette expression croît indéfiniment lorsque la fréquence augmente. L’énergie totale, obtenue en intégrant sur toutes les fréquences, diverge.
Autrement dit, la théorie classique prédit qu’une cavité à température finie devrait contenir une énergie infinie, principalement concentrée dans les hautes fréquences. Cette contradiction flagrante avec l’expérience est connue sous le nom de catastrophe ultraviolette.
C’est précisément ce problème qui conduit Planck, puis Einstein, à introduire une rupture conceptuelle majeure : l’énergie des oscillateurs ne peut pas varier continûment, mais seulement par quanta proportionnels à la fréquence.
Et c’est là que commence véritablement la théorie quantique du rayonnement que vous êtes en train d’analyser.
Dans une cavité cubique de volume $ V=L^3 $, les conditions aux limites imposent
\[ k_x=\frac{n_x\pi}{L},\quad k_y=\frac{n_y\pi}{L},\quad k_z=\frac{n_z\pi}{L}. \] Les points autorisés forment un réseau de pas $ \pi/L $ dans l’espace $ k $.Le nombre de modes dont le vecteur d’onde est compris entre $ k $ et $ k+dk $ correspond au nombre de points du réseau dans une coquille sphérique :
\[ \frac{4\pi k^2 dk}{(\pi/L)^3}. \]On simplifie :
\[ (\pi/L)^3=\frac{\pi^3}{L^3}, \]donc
\[ \frac{4\pi k^2 dk}{(\pi/L)^3} = \frac{4V}{\pi^2} k^2 dk. \]On ne s'intéresse qu'au premier octant de la sphère $ ( n_x,n_y,n_z>0 ) $ On divise donc par 8.
\[ \frac{V}{2\pi^2} k^2 dk. \]On multiplie ensuite par 2 pour les deux polarisations :
\[ \frac{V}{\pi^2} k^2 dk. \]On divise par $ V $ :
\[ \frac{1}{\pi^2} k^2 dk. \]C’est la densité de modes en fonction de $ k $.
Relation de dispersion :
\[ \omega=ck, \qquad \omega=2\pi\nu. \]Donc
\[ k=\frac{2\pi\nu}{c}, \qquad dk=\frac{2\pi}{c} d\nu. \]Alors
\[ k^2 dk = \left(\frac{2\pi\nu}{c}\right)^2 \left(\frac{2\pi}{c} d\nu\right) = \frac{8\pi^3 \nu^2}{c^3} d\nu. \]Théorème de l’équipartition :
\[ \langle E \rangle = k_B T. \]Donc la densité spectrale d’énergie est
\[ u(\nu,T) = \frac{8\pi \nu^2}{c^3} k_B T. \] \[ \boxed{ u(\nu,T) = \frac{8\pi \nu^2}{c^3} k_B T } \]Le « big fail » ne vient pas d’un mauvais comptage géométrique des nœuds en espace $ k $. Le comptage est, en fait, parfaitement correct dans le cadre classique. Le problème vient entièrement de l’hypothèse d’équipartition appliquée à une infinité de modes.
Dans une cavité de volume fini $ V $, les valeurs de $ \vec{k} $ sont discrètes. Mais lorsque $ V $ est macroscopique, l’espacement $ \Delta k = \pi/L $ devient extrêmement petit, et l’ensemble des modes devient quasi continu.
Le nombre de modes de fréquence inférieure à $ \nu $ croît comme
\[ N(\nu) \propto \nu^3. \]Donc la densité spectrale de modes croît comme
\[ g(\nu) \propto \nu^2. \]Il n’y a là aucune erreur mathématique. Cette croissance en $ \nu^2 $ est une conséquence géométrique inévitable de la dimension 3 de l’espace.
Même à très haute fréquence, le nombre de modes continue d’augmenter. Rien, dans l’électromagnétisme classique, n’impose de coupure naturelle.
Le problème surgit lorsqu’on combine cette densité de modes croissante avec le théorème de l’équipartition :
\[ \langle E \rangle = k_B T \]indépendamment de la fréquence.
Autrement dit :
Ce n’est pas une violation de la conservation de l’énergie au sens dynamique.
C’est pire : la théorie prédit qu’un système à température finie possède une énergie totale infinie.
La mécanique statistique classique appliquée à un système possédant une infinité de degrés de liberté quadratiques donne une énergie infinie.
Dans la théorie classique :
En réalité, ce qui se produit expérimentalement n’est pas une diminution du nombre de modes — le nombre de solutions géométriques reste bien proportionnel à $ \nu^2 $.
Ce qui change, c’est l’énergie accessible par mode.
La correction ne consiste pas à modifier le comptage des modes.
Elle consiste à modifier la valeur moyenne de l’énergie par mode.
Planck impose que l’énergie d’un oscillateur ne peut prendre que des valeurs discrètes :
\[ E_n = n h\nu. \]Lorsque $ h\nu \gg k_B T $, l’oscillateur a extrêmement peu de chances d’être excité.
Son énergie moyenne devient alors exponentiellement petite :
\[ \langle E \rangle \sim e^{-h\nu/k_BT}. \]Donc aux hautes fréquences :
La catastrophe ultraviolette n’est donc pas due à un mauvais comptage des modes.
Elle vient de l’application illimitée du théorème d’équipartition à un ensemble infini de degrés de liberté.La théorie classique suppose implicitement :
La physique quantique introduit :
C’est précisément cette suppression exponentielle qui sauve la conservation d’une énergie totale finie.
Le théorème de l’équipartition est un résultat de mécanique statistique classique qui découle de l’intégrale de Boltzmann sur un espace des phases continu. Pour un degré de liberté $x$ qui apparaît de façon quadratique dans l’énergie,
\[ H(x)=\frac{1}{2}a x^2, \]on calcule la moyenne thermique avec la densité $e^{-\beta H}$ (où $\beta = 1/k_B T)$. Le point clé est que $x$ peut prendre toutes les valeurs réelles et que l’on intègre sur un continuum :
\[ \langle H\rangle = \frac{\int_{-\infty}^{+\infty}\frac{1}{2}a x^2,e^{-\beta \frac{1}{2}a x^2},dx}{\int_{-\infty}^{+\infty}e^{-\beta \frac{1}{2}a x^2},dx}. \] Cette intégrale gaussienne donne automatiquement \[ \langle H\rangle=\frac{1}{2}k_B T. \]On n’a pas eu besoin de connaître $a$ : c’est la signature de l’équipartition. Mais ce résultat est vrai parce que l’intégration se fait sur un ensemble continu d’états, avec une densité d’états qui ne “s’éteint” jamais.
Un mode électromagnétique dans une cavité se comporte comme un oscillateur harmonique, avec une coordonnée $q$ et une impulsion $p$ :
\[ H(p,q)=\frac{p^2}{2m}+\frac{1}{2}m\omega^2 q^2. \]Il y a deux termes quadratiques indépendants, donc l’équipartition donne
\[ \left\langle \frac{p^2}{2m}\right\rangle=\frac{1}{2}k_B T, \qquad \left\langle \frac{1}{2}m\omega^2 q^2\right\rangle=\frac{1}{2}k_B T, \] et donc \[ \langle E\rangle=\langle H\rangle=k_B T. \]C’est exactement l’hypothèse qui, combinée au comptage des modes $g(\nu)\propto \nu^2$, mène à Rayleigh–Jeans et à la divergence ultraviolette.
Dans le cadre classique, pour un oscillateur de fréquence très grande ($\omega$ énorme), il n’existe aucun seuil : on peut exciter le mode avec une énergie aussi petite qu’on veut, simplement en prenant des amplitudes $q$ minuscules. Le mode “existe” et l’espace des phases $(p,q)$ autour de l’origine est toujours accessible.
C’est contre-intuitif physiquement parce que, expérimentalement, les modes de haute fréquence ne se remplissent pas à température modérée. Mais dans la statistique classique, rien n’empêche qu’un très grand nombre de modes reçoivent chacun une énergie moyenne $k_BT$. Comme il y a une infinité de modes, l’énergie totale diverge.
Ce n’est pas une violation de la conservation de l’énergie au sens “le système créerait de l’énergie”. C’est plutôt une prédiction absurde : si l’on applique strictement ces hypothèses à un champ ayant une infinité de degrés de liberté, alors l’énergie d’équilibre est infinie, donc le modèle ne peut pas être correct.
Le geste de Planck peut se reformuler ainsi : l’espace des états d’un oscillateur n’est pas continu de la même manière que dans la mécanique classique. Les états accessibles ne forment plus un continuum d’énergies, mais une série discrète :
\[ E_n = n h\nu \quad (\text{ou } E_n=(n+\tfrac12)h\nu \text{ dans la formulation quantique complète}). \]Du point de vue statistique, cela remplace une intégrale continue par une somme sur des niveaux :
\[ \langle E\rangle = \frac{\sum_{n=0}^{\infty} E_n e^{-\beta E_n}} {\sum_{n=0}^{\infty} e^{-\beta E_n}}. \]Avec $E_n=n h\nu$, cette somme se calcule exactement et donne :
\[ \boxed{\;\langle E\rangle=\frac{h\nu}{e^{h\nu/k_BT}-1}\;} \](la version de Planck ; si on inclut le $\tfrac12 h\nu$, cela ajoute simplement une énergie de point zéro qui ne contribue pas au rayonnement thermique mesuré de la même façon).
L’idée physique profonde est qu’il apparaît un coût minimal pour exciter un mode : l’énergie caractéristique est $h\nu$. Quand $\nu$ est très grande, ce coût devient très grand par rapport à $k_BT$, donc la probabilité d’occupation chute.
Regardez le régime ultraviolet $h\nu \gg k_B T$. Alors
\[ e^{h\nu/k_BT} \text{ est gigantesque,} \quad\Rightarrow\quad \langle E\rangle \approx h\nu,e^{-h\nu/k_BT}. \]Donc même si le nombre de modes croît comme $\nu^2$, l’énergie moyenne par mode décroît beaucoup plus vite (exponentiellement). Le produit devient intégrable :
\[ u(\nu,T)=g(\nu)\langle E\rangle = \frac{8\pi \nu^2}{c^3}\cdot \frac{h\nu}{e^{h\nu/k_BT}-1}. \]Et l’énergie totale
\[ \int_0^\infty u(\nu,T),d\nu \]devient finie.
Il y a une façon très parlante de relier cela à l’espace des phases, sans entrer dans toute la mécanique quantique : la quantification revient à dire qu’on ne peut pas découper l’espace ((p,q)) en cellules arbitrairement petites. Il existe une “maille” naturelle de taille de l’ordre de $h$ (plus précisément $h$ ou $2\pi\hbar$ selon les conventions). Au lieu d’une intégrale continue donnant l’équipartition pour tous les $\omega$, on obtient un nombre discret d’états accessibles à une énergie donnée, ce qui modifie radicalement la moyenne thermique.
C’est exactement la raison pour laquelle l’équipartition, qui est un théorème d’intégrale continue, cesse d’être applicable aux modes de haute fréquence : le continuum d’états qu’elle suppose n’existe plus de la même manière.
Imaginons une cavité parfaitement conductrice, initialement vide de tout rayonnement électromagnétique à l’instant $ t=0 $. Les parois sont maintenues à une température constante $ T $. La question n’est plus seulement statistique mais dynamique : comment, selon la théorie classique, le champ électromagnétique apparaît-il dans la cavité et comment atteint-il un état d’équilibre thermique ?
Dans la description classique, rien n’est mystérieux au départ. Les parois métalliques sont constituées de charges — électrons libres et ions du réseau cristallin — en agitation thermique permanente. Or toute charge accélérée rayonne. Les mouvements microscopiques désordonnés des charges produisent donc un rayonnement électromagnétique qui pénètre dans la cavité. Même si le champ est nul initialement, il commence ainsi à se développer.
Le champ électromagnétique dans la cavité peut être décomposé en modes propres, c’est-à-dire en ondes stationnaires compatibles avec les conditions aux limites imposées par les parois. Chaque mode se comporte dynamiquement comme un oscillateur harmonique, dont l’amplitude évolue dans le temps selon une équation du type
\[ \ddot q + \omega^2 q = 0. \]Les fluctuations thermiques des charges dans les parois possèdent un spectre fréquentiel large. Certaines composantes fréquentielles correspondent aux fréquences propres des modes de la cavité. L’énergie est alors progressivement transférée vers ces modes, dont les amplitudes croissent.
Le processus n’est pas unilatéral. Les parois n’émettent pas seulement ; elles absorbent également le rayonnement présent dans la cavité. Le champ électrique des modes exerce une force sur les charges des parois, qui peuvent absorber cette énergie et la redistribuer sous forme thermique par collisions et dissipation.
Il existe donc un échange continu d’énergie entre les parois et les modes du champ. Lorsque l’émission moyenne compense l’absorption moyenne, l’équilibre thermique est atteint.
Dans le cadre de la mécanique statistique classique, chaque mode du champ est un oscillateur harmonique possédant deux termes quadratiques dans son hamiltonien. Le théorème de l’équipartition impose alors que l’énergie moyenne par mode soit
\[ \langle E \rangle = k_B T, \]indépendamment de la fréquence du mode.
Le nombre de modes disponibles dans la cavité croît avec la fréquence comme $ \nu^2 $. Autrement dit, plus la fréquence augmente, plus il existe de modes capables d’absorber de l’énergie.
Si chaque mode reçoit en moyenne une énergie $ k_B T $, alors l’énergie totale contenue dans le champ devient proportionnelle à
\[ \int_0^{\infty} \nu^2 k_B T , d\nu, \]qui diverge.
Du point de vue dynamique, cela signifie que les parois devraient continuellement transférer de l’énergie vers des modes de plus en plus haute fréquence. Aucune limite naturelle n’existe dans la théorie classique pour empêcher l’excitation de ces modes. Il est toujours possible, classiquement, d’exciter un oscillateur de fréquence très élevée avec une amplitude arbitrairement petite, donc avec une énergie arbitrairement faible. L’absence de seuil énergétique rend possible l’alimentation indéfinie d’un nombre infini de degrés de liberté.
La conservation de l’énergie n’est pas violée localement ; c’est la prédiction globale d’équilibre qui devient absurde. La théorie classique implique qu’un système à température finie devrait contenir une énergie infinie.
L’expérience montre pourtant qu’un équilibre stable est atteint et que le spectre d’émission décroît rapidement aux hautes fréquences. Les modes de fréquence élevée ne se remplissent pas indéfiniment ; leur contribution devient négligeable.
Il apparaît donc que le mécanisme d’échange classique est incomplet. La difficulté ne provient ni du comptage des modes ni des équations de Maxwell, mais de l’hypothèse selon laquelle l’énergie d’un oscillateur peut varier continûment sans seuil minimal.
Planck introduit une hypothèse radicale : l’énergie d’un oscillateur de fréquence $ \nu $ ne peut prendre que des valeurs discrètes,
\[ E_n = n h\nu. \]Il existe désormais un coût minimal pour exciter un mode. Lorsque $ h\nu \gg k_B T $, l’énergie thermique typique disponible dans les parois est insuffisante pour peupler le premier niveau excité. Le mode reste alors pratiquement dans son état fondamental.
Ainsi, les hautes fréquences sont thermiquement « gelées ». L’énergie moyenne par mode décroît exponentiellement lorsque la fréquence augmente. La divergence disparaît et l’équilibre devient stable.
Dans la description classique, l’espace des amplitudes possibles est continu et toujours accessible, même à énergie arbitrairement faible. En introduisant un quantum d’action $ h $, la quantification impose une granularité fondamentale : tous les états énergétiques ne sont plus accessibles continûment.
La catastrophe ultraviolette révèle donc une limite intrinsèque de la mécanique classique appliquée à un champ possédant une infinité de degrés de liberté. Ce n’est pas l’électromagnétisme qui est en défaut, mais l’hypothèse de continuité de l’énergie.
C’est précisément cette rupture conceptuelle qui ouvre la voie à la théorie quantique du rayonnement, puis à l’idée plus radicale encore, formulée par Einstein, que la lumière elle-même est constituée de quanta.
Audio
Video
Présentation PDF.
A. Einstein