Équation de Boltzmann-Peierls

En physique du solide l'équation de Boltzmann-Peierls décrit l'évolution de la fonction de distribution des phonons dans un solide cristallin. Elle a été établie par Rudolf Peierls en 1929. Diverses approximations de type BGK ont été proposées, la plus précise par Joseph Callaway en 1959. Dans les dernières décennies des efforts ont porté sur le développement de méthodes analogues à celles de la mécanique des fluides, permettant une grande efficacité de la résolution du problème moyennant quelques simplifications acceptables pour les applications pratiques.

Cette équation décrit le transfert d'énergie thermique au niveau microscopique, dans un domaine où l'approximation diffusive (loi de Fourier) n'est pas valide, le libre parcours moyen des phonons étant du même ordre de grandeur ou plus grand que la longueur caractéristique associée à la température $\frac{T}{||\mathrm{\nabla }T||}$. Ce type de problème est présent en microélectronique et la méthode s'applique également aux électrons dans les solides.

L'équation de Boltzmann-Peierls

Les phonons sont des bosons décrits par leur densité numérique $f_{𝐤,s}(𝐱,t,𝐤)$ où $𝐱$ est la variable d'espace, $t$ le temps, $𝐤$ le vecteur d'onde supposé continu et $s=1,2,3$ désigne la polarisation. Cette grandeur est liée à la densité d'états ${\rho }_{𝐤,s}(k)$ caractéristique du matériau et à la distribution statistique de ces états $n_{𝐤,s}(𝐱,t)$ par $f_{𝐤,s}={\rho }_{𝐤,s}{n}_{𝐤,s}$.

Dans ce qui suit on s'intéressera à une polarisation particulière, les calculs étant les mêmes pour toutes : l'indice s disparaît donc.

L'équation de conservation de $f_{𝐤}$ s'écrit :

$\frac{\mathrm{d}{f}_{𝐤}}{\mathrm{d}t}=\frac{\partial f_{𝐤}}{\partial t}+{𝐯}_{𝐤}\cdot \mathrm{\nabla }{f}_{𝐤}=\frac{\mathrm{d}{f}_{𝐤}}{\mathrm{d}T}(\frac{\partial T}{\partial t}+{𝐯}_{𝐤}\cdot \mathrm{\nabla }T)=S_{𝐤}$

où $T$ est la température thermodynamique et ${𝐯}_{𝐤}$ est la vitesse de propagation ou vitesse de groupe liée à la pulsation de l'onde ${\omega }_{𝐤}$ par ${𝐯}_{𝐤}=\frac{\mathrm{d}{\omega }_{𝐤}}{\mathrm{d}k}\frac{𝐤}{k}$. La relation ${\omega }_{𝐤}(k)$ constitue la relation de dispersion.

La connaissance du terme source $S_{𝐤}$ est au centre du problème. Prenons l'exemple simple de la diffusion d'un phonon sur un centre diffuseur comme un défaut du réseau cristallin. La probabilité de transition $𝐤\to {𝐤}^{\prime }$ est proportionnelle à ${\mathrm{\Lambda }}_{𝐤\to {𝐤}^{\prime }}{n}_{𝐤}(1+{n_{𝐤}}^{\prime })$ où ${\mathrm{\Lambda }}_{𝐤\to {𝐤}^{\prime }}$ est la probabilité de transition spontanée.

En prenant en compte la réversibilité ${\mathrm{\Lambda }}_{𝐤\to {𝐤}^{\prime }}={\mathrm{\Lambda }}_{{𝐤}^{\prime }\to 𝐤}$ on obtient le terme source pour la diffusion :

$S_{𝐤}=\frac{1}{2}\sum _{𝐤}\sum _{{𝐤}^{\prime }}{\mathrm{\Lambda }}_{𝐤\to {𝐤}^{\prime }}(n_{𝐤}-n_{{𝐤}^{\prime }})$

le facteur $\frac{1}{2}$ étant destiné à corriger les doubles comptes.

On peut faire de même pour les processus à 3 phonons dont les probabilités de transition sont plus complexes car faisant intervenir une approximation à l'ordre 3 au moins du potentiel interatomique cristallin (il n'y a pas d'interaction phonon-phonon à l'ordre 2).

Le calcul de ces termes est coûteux et la résolution de l'équation de Boltzmann par une méthode stochastique ou déterministe est elle-même coûteuse (calcul d'un grand nombre de fonctions de distribution en tout point et à tout instant, un terme source complexe). Des approximations ont donc été recherchées.

Première approximation

On peut écrire en première approche le terme source sous une forme utilisant la méthode de Bhatnagar-Gross-Krook comme un simple terme de rappel vers la valeur d'équilibre $f_{𝐤}^{(0)}={\rho }_k{n}_{𝐤}^{(0)}$ où $n_{𝐤}^{(0)}$ est donné par la statistique de Bose-Einstein pour un potentiel chimique nul :

$n_{𝐤}^{(0)}={(e^{\frac{\hslash {\omega }_{𝐤}}{k_{B}T}}-1)}^{-1}$

Ce rappel s'effectue avec le temps caractéristique ${\tau }_{𝐤}$ :

$S_{𝐤}=-\frac{f_{𝐤}-f_{𝐤}^{(0)}}{{\tau }_{𝐤}}$

Cette approximation est bonne lorsque l'on traite des processus dits « résistifs » (car participant à la résistivité thermique) qui ne conservent pas la quantité de mouvement comme la diffusion dite processus umklapp ou pour les phénomènes de diffusion d'un phonon sur un centre diffuseur (un défaut du cristal).

Lorsque plusieurs phénomènes se superposent ce temps caractéristique est obtenu en appliquant la règle de Matthiessen qui exprime l'indépendance des phénomènes :

${\tau }^{-1}=\sum _i{\tau }_i^{-1}$

Deuxième approximation (Callaway)

Si, en plus des processus résistifs, il existe des processus dits « normaux » qui respectent la conservation de la quantité de mouvement une expression différente est proposée par Callaway :

$\frac{\partial f_{𝐤}}{\partial t}+{𝐯}_{𝐤}\cdot \mathrm{\nabla }{f}_{𝐤}=-\frac{f_{𝐤}-f_{𝐤}^{(0)}}{{\tau }_{𝐤}}-\frac{f_{𝐤}-f_{𝐤}^{(1)}}{{\tau }_{𝐤}^{*}}$

où $f_{𝐤}^{(1)}$ est la distribution de Bose-Einstein en repère mobile à la vitesse de dérive $𝐮(𝐱,t)$ :

$f_{𝐤}^{(1)}=\frac{1}{e^{\frac{\hslash }{k_{B}T}({\omega }_{𝐤}+𝐮\cdot 𝐤)}-1}\approx n_{𝐤}^{(0)}+\frac{\hslash 𝐮\cdot 𝐤}{k_{B}T}\frac{e^{\frac{\hslash {\omega }_{𝐤}}{k_{B}T}}}{{(e^{\frac{\hslash {\omega }_{𝐤}}{k_{B}T}}-1)}^2}$

Au contraire de la distribution de Bose-Einstein cette distribution n'est pas isotrope. $f_{𝐤}^{(1)}$ correspond à une distribution de Bose-Einstein observée dans un repère lorentzien se déplaçant à la vitesse $𝐮$ par rapport au solide. Ce terme participe à la propagation et donc indirectement à la résistivité thermique.

Méthodes aux moments

À une époque plus récente les physiciens se sont intéressés aux méthodes utilisées en mécanique des fluides et en transfert radiatif, particulièrement aux méthodes basées sur l'utilisation des moments de la fonction de distribution comme la méthode de Grad ou la méthode MN. On peut remarquer que, dans l'hypothèse du continu, $\frac{𝐤}{k}$ est une distribution angulaire et donc que la quantité $I_{𝐤,s}=\hslash {\omega }_{𝐤}{v}_{𝐤}{f}_{𝐤,s}\frac{𝐤}{k}$ est l'exacte contrepartie de la luminance spectrale du rayonnement de photons.

Moments de l'équation

Prenons l'exemple simple d'un milieu isotrope dans lequel on ne retient que les ondes acoustiques dont la vitesse de propagation $v$ est indépendante du nombre d'onde (approximation $k\to 0$) et connue. La relation de dispersion se réduit à ${\omega }_{𝐤}=vk$. Dans ce qui suit on ignore les polarisations.

On définit les moments suivants ($\otimes$ est le produit tensoriel) :

$E=\underset{𝐤}{\int }\hslash kv{f}_{𝐤}\mathrm{d}𝐤$ énergie volumique ;

$𝐅=\underset{𝐤}{\int }\hslash k{v}^2{f}_{𝐤}\frac{𝐤}{k}\mathrm{d}𝐤$ flux d'énergie, et sa valeur réduite $𝐟=\frac{𝐅}{vE}$ ;

$𝖯=\underset{𝐤}{\int }\hslash kv{f}_{𝐤}\frac{𝐤}{k}\otimes \frac{𝐤}{k}\mathrm{d}𝐤$ tenseur de pression et sa valeur normalisée (tenseur d'Eddington) $𝖣=\frac{𝖯}{E}$.

Ces quantités sont les moments de Hausdorff d'une quantité non-négative sur la sphère unité. On montreModèle:, que $𝐅$ est un vecteur propre de $𝖯$ d'où :

$𝖣𝐟=\chi 𝐟$

où $\chi$ est le facteur d'Eddington.

Le tenseur d'Eddington s'exprime en fonction du facteur d'Eddington par :

$𝖣=\frac{1-\chi }{2}𝖨+\frac{3\chi -1}{2}\frac{𝐟}{f}\otimes \frac{𝐟}{f},\frac{1}{3}\le \chi \le 1$

où $𝖨$ est le tenseur unité. $\chi =\frac{1}{3}$ correspond à l'isotropie et $\chi =1$ à un faisceau parallèle. La solution générale possède la symétrie de révolution autour de $𝐟$.

On obtient les équations aux moments de la façon suivante :

• On multiplie l'équation de Boltzmann-Peierls par l'énergie du phonon $\hslash kv$ et on intègre sur k :

$\underset{𝐤}{\int }\frac{\partial f_{𝐤}}{\partial t}\hslash kv\mathrm{d}𝐤+\underset{𝐤}{\int }\hslash k{v}^2\frac{𝐤}{k}\cdot \mathrm{\nabla }{f}_{𝐤}\mathrm{d}𝐤=\underset{𝐤}{\int }{S}_{𝐤}\hslash kv\mathrm{d}𝐤=0$

Le second terme de cette équation est nul du fait de la conservation de l'énergie de tous les processus.

Par identification on obtient :

$\frac{\partial E}{\partial t}+\mathrm{\nabla }\cdot 𝐅=0$

• On multiplie l'équation de Boltzmann-Peierls par $\hslash kv\frac{𝐤}{k}$ et on intègre sur k :

$\underset{𝐤}{\int }\frac{\partial f_{𝐤}}{\partial t}\hslash kv\frac{𝐤}{k}\mathrm{d}𝐤+\underset{𝐤}{\int }\hslash k{v}^2\frac{𝐤}{k}\cdot \mathrm{\nabla }{f}_{𝐤}\frac{𝐤}{k}\mathrm{d}𝐤=\underset{𝐤}{\int }{S}_{𝐤}\hslash kv\frac{𝐤}{k}\mathrm{d}𝐤$

Dans le second membre les processus élastiques conduisent à un bilan nul. On suppose de plus que tous les processus résistifs sont décrits par une approximation BGK à une seule constante de temps $\tau$. Alors par identification on obtient :

$\frac{\partial 𝐅}{\partial t}+v^2\mathrm{\nabla }\cdot 𝖯=-\frac{𝐅}{\tau }$

Comme tout système aux moments il faut à présent faire une hypothèse pour « fermer » (compléter) le système par nature incomplet puisque l'on a deux équations pour les trois inconnues $E,𝐟,\chi$.

Fermeture entropique

On va comme ci-dessus maximiser l'entropie en respectant les moments définis ci-dessus. Cette méthode constitue la méthode M₁. La solution est analytique :

$\chi =\frac{3+4f^2}{5+2\sqrt{4-3f^2}}$

Le système à présent complet est hyperbolique et l'on peut utiliser pour le résoudre les méthodes très efficaces qui ont été développées pour la mécanique des fluides.

On peut reconstruire la fonction de distribution qui n'apparaît pas explicitement dans le problème :

$n_{𝐤}=\frac{1}{e^{\frac{\hslash kv}{k_{B}T}[g_1(f)+g_2(f)𝐟\cdot 𝐤]}-1}$

où g₁ et g₂ sont analytiques. Cette expression comporte une « vitesse de dérive » $𝐮=g_2(f)𝐟$ qui s'interprète comme la manifestation locale du gradient de peuplement de phonons dans chaque direction, avec un maximum dans le sens du flux d'énergie.

L'analyse asymptotique du système aux dérivées partielles permet d'exhiber une limite diffusive :

$𝐅=-\frac{4av{T}^3\tau }{3}\mathrm{\nabla }T$

où a est l'équivalent de la constante radiative modifiée pour prendre en compte le nombre de polarisations possibles (3 au lieu de 2 pour les photons) $a=\frac{12{\pi }^5{k}_B^4}{15h^3{v}^3}$

En identifiant le coefficient avec la conductivité thermique $\lambda$ on obtient le temps caractéristique :

$\tau =\frac{3\lambda }{4av{T}^3}$

Résolution

Transfert de phonons dans une lame de siliciumModèle:, . Comparaison calcul-mesure. « Fourier » désigne le calcul de conduction classique.

Au final on est ramené à un système aux dérivées partielles hyperbolique comportant une équation scalaire sur E et une équation vectorielles sur $𝐟$ pour lequel il existe de nombreux solveurs numériques efficaces développés dans le cadre de la mécanique des fluides.

Il y a cependant une difficulté à écrire les conditions aux limites au bord du matériau. On remarquera d'abord que la notion de phonon repose sur l'hypothèse d'un milieu infini. On supposera toutefois que l'équation de Boltzmann-Peierls reste valide au voisinage de la surface. En l'absence de transmission (on écarte le problème de l'interface entre deux matériaux) on écrit les conditions aux limites comme on le fait dans le cas de la couche de Knudsen :

• la réflexion est soit spéculaire, soit diffuse. On peut également superposer les deux types de réflexion en les pondérant par un coefficient d'« accomodation » ;
• on peut également imposer une distribution d'équilibre à la paroi : cette condition, difficile à justifier physiquement, est utilisée pour étudier des problèmes spécifiques.

Cette approche est également valide pour l'interface entre deux matériaux (résistance de Kapitza).