Phonon
Modèle:Voir homonymes En physique, un phonon correspond à une excitation collective dans un arrangement périodique d'atomes constituant une structure cristalline ou amorphe. La déformation est élastique. L'onde qui se propage peut être assimilée à une quasi-particule.
Ils permettent d'expliquer les propriétés physiques des solides :
- la capacité thermique ;
- la conductivité thermique ;
- la capacité à propager le son ;
- la dilatation thermique.
Définition
Le concept de phonon<ref>Modèle:Article</ref> a été créé par Igor TammModèle:Sfn en Modèle:DateModèle:Sfn,<ref>Modèle:Article</ref> et le mot « phonon » (du grec ancien Modèle:Grec ancien / phonê, la voix) a été inventé par Yakov FrenkelModèle:Sfn en Modèle:DateModèle:Sfn,<ref>Modèle:Article</ref>. Le suffixe « -on », qui apparaît dans le nom de nombreuses entités de la physique de la matière condensée (excitons, magnons, etc.) a été calqué sur la fin du mot « électron » (mot inventé par George Stoney en 1891)<ref>Modèle:Article</ref>.
Un phonon est la description quantique du mouvement vibratoire élémentaire dans lequel un réseau d'atomes ou d'ions oscille uniformément à une fréquence donnée<ref>Modèle:Ouvrage</ref>. En mécanique classique cela correspond à un mode normal de vibration. L'analyse de Fourier permet de décomposer toute vibration du réseau comme une superposition de ces modes de vibration élémentaires. Alors que les modes normaux de la mécanique classique sont traités comme des phénomènes ondulatoires, les phonons peuvent être assimilés à une particule de la mécanique quantique.
Dynamique du réseau cristallin
Dans un solide, il existe des forces d'interaction (essentiellement les forces covalentes ou ioniques dans notre cas) qui maintiennent chaque atome près d'une position d'équilibre. L'interaction entre chaque paire d'atomes peut être caractérisée par une fonction d'énergie potentielle <math>\textstyle V_{ij}</math> qui ne dépend que de la distance <math>\textstyle r_{ij}</math> entre chacun des couples de ces <math>\textstyle N</math> atomes et de leur nature. L'énergie potentielle du réseau dans son ensemble est la somme des énergies potentielles d'interaction de chaque paire :
- <math>E = \sum_{i=1}^N \sum_{j=i+1}^N V_{ij}(r_{ij})</math>
Cette expression, caractéristique d'un problème à N corps, ne se prête pas à une résolution que ce soit en mécanique classique ou en mécanique quantique. Il est donc nécessaire de procéder à des approximations pour poursuivre l'analyse. Les trois approximations généralement employées sont :
- restreindre l'étude à un nombre limité d'atomes en utilisant des conditions aux limites périodiques. Cette simplification appelée hypothèse de Born-von Kármán suppose implicitement que l'on saura traiter le cas <math>\textstyle N\to\infty</math>. Une conséquence peu mise en avant est le fait qu'un tel modèle est inapte à traiter les conditions au bord du solide<ref name=Klemens>Modèle:Ouvrage</ref> ;
- restreindre la sommation aux atomes voisins (approximation des voisins les plus proches). En effet, bien que rigoureusement les forces électriques aient une portée infinie dans un solide réel, cette approximation est valide car les forces s'exerçant sur des atomes éloignés sont écrantées et donc négligeables ;
- considérer que le potentiel V est un potentiel harmonique, ce qui est valide lorsque les atomes restent proches de leurs positions d'équilibre. Formellement, cette hypothèse s'applique en effectuant un développement de Taylor du potentiel autour de la valeur d'équilibre et en ne gardant que les termes jusqu'à l'ordre 2. Dans certains cas cette approximation est insuffisante et l'on doit monter en ordre au prix d'une plus grande complexité du problème.
Bien entendu l'analyse prend en compte le fait qu'il s'agit d'un réseau cristallin qui reproduit un certain nombre de cellules de base : le réseau de Bravais.
Pour un potentiel <math>\textstyle V_{ij}</math> harmonique le potentiel global s'écrira en ne retenant que les interactions entre proches voisins (noté pv) :
- <math>E = \frac{1}{2} \sum_{i=1}^N \sum_{pv} \mu_{ij}\omega_{ij}^2R_{ij}^2</math>
où :
- <math>\textstyle \mu_{ij}=\left(m_i^{-1}+m_j^{-1}\right)^{-1}</math> est la masse réduite ;
- <math>\textstyle \omega_{ij}</math> est la pulsation propre au couple (i,j) ;
- <math>\textstyle R_{ij}</math> est la distance <math>\textstyle r_{ij}</math> pour un réseau cristallin au repos.
Réseau unidimensionnel
On traite ci-dessous le réseau à une seule dimension qui permet de comprendre les concepts du domaine.
Approche classique
On étudie donc un tel réseau, homogène (un seul type d'atome), comportant N atomes de masse m positionnés en <math>\textstyle \;x_i\,,\; i=1..N\;</math> liés entre eux par un potentiel harmonique représenté par un ressort de constante élastique C. À l'équilibre la distance entre atomes est a. Avec l'approximation de périodicité et celle de la seule influence des voisins proches, on peut écrire la dynamique du système physique par le système d'équation cinématiques suivant :
- <math>m\frac{d^2 x_i}{dt^2}=C(x_{i+1}-x_{i-1})-2Cx_i\,,\quad i=1..N\,,\quad x_0=x_N\,,\;x_{N+1}=x_1</math>
La solution étant supposée oscillatoire on utilise pour obtenir un découplage des équations du système une transformation de Fourier discrète<ref>Modèle:Ouvrage</ref>. On définit les coordonnées normales :
- <math>Q_k=\sum_{i=1}^N x_i e^{-j\frac{2\pi}{N}ka} \;\Longleftrightarrow\; x_i = \sum_{k=1}^N Q_k e^{jkia}\,, \quad k = {n\pi \over Na}\,,\quad n = 1, ... , N</math>
où j est l'unité imaginaire.
On obtient un système d'équations découplées<ref name=Misra/> :
- <math>\frac{d^2 Q_k}{dt^2}=2\Omega^2(\cos ka-1)Q_k\,,\quad \Omega=\sqrt{\frac{C}{m}}</math>
La solution de chaque équation s'écrit comme une fonction sinusoïdale d'amplitude <math>\textstyle A_k</math> :
- <math>X_k=A_ke^{i\omega_k t}\,,\quad \omega_k=\Omega\sqrt{2(1-\cos ka)}=2\Omega\left|\sin\frac{ka}{2}\right|</math>
Elle constitue un mode normal de nombre d'onde k et de pulsation <math>\textstyle \omega_k</math>.
La relation <math>\textstyle k\to \omega_k</math> constitue la relation de dispersion.
Dans la limite du continu obtenue pour <math>\textstyle N\to\infty\,,a\to 0\,,Na=C^{ste}</math> on obtient un champ d'oscillateurs harmonique de pulsation <math>\textstyle \omega(k)=\Omega ak</math>.
Approche quantique
L'hamiltonien du système s'écrit :
- <math> \mathbf{H} = \sum_{i=1}^N {P_i^2 \over 2m} + {1\over 2} m \omega^2 \sum_{i=1}^N\sum_{pv} (X_i - X_j)^2 </math>
où
- <math>\textstyle X_i</math> est l'opérateur position pour l'atome i ;
- <math>\textstyle P_i</math> est l'opérateur impulsion.
On définit :
- 2N+1 « coordonnées normales » Qk définies comme les transformées de Fourier des opérateurs position Xi ;
- 2N+1 « moments conjugués » Πk définis comme les transformées de Fourier des opérateurs impulsion Pi.
- <math>\begin{matrix}
X_i &=& {1\over\sqrt{N}} \sum_{k} Q_k e^{jkia} &\,, \quad k = {n\pi \over Na}\,, \quad n = 0, \pm 1, \pm 2, ... , \pm N \\ P_i &=& {1\over\sqrt{N}} \sum_{k} \Pi_k e^{-jkia}& \\ \end{matrix}</math>
k est le nombre d'onde du phonon qui prend des valeurs quantifiées résultant du choix de la périodicité du système.
En inversant la transformée de Fourier pour exprimer les Qk en fonction des Xi et les Πk en fonction des Pi, et en utilisant les relations de commutations canoniques entre les Xi et les Pi de la mécanique quantique :
- <math>\left[ Q_k , \Pi_{k'} \right] = i \hbar \delta_{k k'} \,,\quad
\quad \left[ Q_k , Q_{k'} \right] = \left[ \Pi_k , \Pi_{k'} \right] = 0</math>
En d'autres termes, les coordonnées normales Qk et leurs moments conjuguées Πk obéissent aux mêmes relations de commutation que les opérateurs position Xi et impulsion Pi. En fonction de ces grandeurs, l'hamiltonien s'écrit :
- <math> \mathbf{H} = \sum_k \left(
{ \Pi_k\Pi_{-k} \over 2m } + {1\over2} m \omega_k^2 Q_k Q_{-k}\right) </math>
avec
- <math> \omega_k = \Omega\sqrt{2 (1 - \cos ka )}=2\Omega\left|\sin\frac{ka}{2}\right| </math>
Le couplage entre les variables positions a été transformé. Si les Qk et les Πk étaient hermitiens (ce qui n'est pas le cas), l'hamiltonien transformé décrirait N oscillateurs harmoniques non couplées. En fait, cet hamiltonien décrit une théorie quantique des champs de bosons découplés.
Le spectre des énergies propres de cet hamiltonien est obtenu en utilisant les opérateurs adjoints création et annihilation de quantum d'énergie ak† et ak définis comme :
- <math>\begin{matrix}
a_k &=& \sqrt{m\omega_k \over 2\hbar} (Q_k + {i\over m\omega_k} \Pi_{-k}) \\ a_k^\dagger &=& \sqrt{m\omega_k \over 2\hbar} (Q_{-k} - {i\over m\omega_k} \Pi_k) \end{matrix}</math>
Ces opérateurs adjoints vérifient :
- <math>\mathbf{H} = \sum_k \hbar \omega_k \left(a_k^{\dagger}a_k + 1/2\right) </math>
- <math>[a_k , a_{k'}^{\dagger} ] = \delta_{kk'}</math>
- <math>[a_k , a_{k'} ] = [a_k^{\dagger} , a_{k'}^{\dagger} ] = 0.</math>
Comme dans le cas de l'oscillateur harmonique quantique, on peut montrer que les opérateurs ak† et ak correspondent respectivement à la création et l'annihilation d'une excitation d'énergie ℏωk. Cette excitation est un phonon<ref>Définition sur le site Larousse.</ref>, correspondant à l'intervalle d'énergie relative à un phonon entre deux niveaux du mode k :
- <math>E_n=\left(n+\tfrac12\right)\hbar \omega_k \,,\quad n=1,2,3...</math>
On peut immédiatement en déduire deux propriétés importantes des phonons. D'abord, les phonons sont des bosons : tout mode de vibration k peut voir son nombre d'excitation augmenter par l'application répétée de l'opérateur création ak†. En second lieu, chaque phonon est un « mode collectif » correspondant au mouvement de la (quasi) totalité des atomes du réseau. Cette seconde conclusion se voit dans le fait que les opérateurs adjoints contiennent des sommations sur les positions et les impulsions de tous les atomes du réseau.
On peut voir que les excitations générées par les opérateurs adjoints sont des ondes de déplacement d'atomes du réseau en calculant la fonction de corrélation position-position. Soit <math>\textstyle | k \rangle</math> un état pour lequel un seul quantum de mode k est excité, c’est-à-dire :
- <math>\begin{matrix}
| k \rangle = a_k^\dagger | 0 \rangle \end{matrix}</math>
On peut alors montrer que pour deux atomes i et j quelconques :
- <math>\langle k | x_i(t) x_j(0) | k \rangle = \frac{\hbar}{Nm\omega_k} \cos \left[ k(i-j)a - \omega_k t \right] + \langle 0 | x_i(t) x_j(0) |0 \rangle </math>
ce qui décrit une onde du réseau de pulsation ωk et de nombre d'onde k.
Réseau tridimensionnel
La généralisation à trois dimensions du modèle unidimensionnel précédent est aisée (mais assez lourde). Le nombre d'onde k est remplacé par un vecteur à trois dimensions, le vecteur d'onde <math> \vec k </math>. De plus, <math> \vec k </math> est maintenant associé à trois coordonnées normales. L’hamiltonien a la forme :
- <math>
\mathbf{H} = \sum_k \sum_{s=1}^3 \hbar \, \omega_{k,s} \left( a_{k,s}^{\dagger}a_{k,s} + \frac{1}{2} \right) </math>
Le nouvel indice s = 1, 2, 3 correspond à la polarisation des phonons lorsque la maille primitive comporte un seul type d'atome. En effet, dans un modèle unidimensionnel, les atomes ne peuvent vibrer que sur une ligne, et tous les phonons correspondent à une onde longitudinale. En revanche en trois dimensions, la vibration ne se fait plus uniquement dans la direction de propagation, mais peut également lui être perpendiculaire. Elle correspond alors à une onde transverse. Cela donne naissance à des coordonnées normales supplémentaires, qui comme l'expression de l'hamiltonien l'indique, correspondent à des espèces indépendantes de phonons.
Si la maille primitive comporte p types d'atomes, il y aura au total 3 modes acoustiques et 3p-3 modes optiques.
Courbe de dispersion
Dans le cas d'un réseau monoatomique la relation de dispersion <math>\textstyle \omega=f(k)</math> est périodique, de période <math>\textstyle \frac{2\pi}{a}</math>. Elle est nulle en <math>\textstyle k=0</math> et approximativement linéaire à partir de ce point pour les faibles valeurs de k. Dans cette région les vitesses de phase et de groupe sont égales <math>\textstyle \frac{d\omega}{dk}=\frac{\omega}{k}=\Omega a</math>. La valeur <math>\textstyle k=0</math> correspond à une longueur d'onde infinie : le réseau cristallin de déplace en bloc.
Dans le cas général il existe plusieurs atomes différents dans le réseau. On peut calculer la relation de dispersion analytiquement dans les cas d'un réseau unidimensionnel alterné à deux atomes de masses <math>\textstyle m_1</math> et <math>\textstyle m_2</math><ref name=Misra>Modèle:Ouvrage</ref> :
- <math>\omega_\pm^2 = C\left(\frac{1}{m_1} +\frac{1}{m_2}\right) \pm C \sqrt{\left(\frac{1}{m_1} +\frac{1}{m_2}\right)^2-\frac{4\sin^2\frac{ka}{2}}{m_1 m_2}}</math>
Il existe deux branches :
- la branche inférieure, analogue à celle de la chaîne monoatomique et correspondante aux grandes longueurs d'onde est appelée branche acoustique ;
- la branche supérieure est appelée branche optique car elle correspond à des longueurs d'onde dans l'infrarouge. Cette branche correspondant à des oscillations en opposition de phase des voisins, elle est apte à générer un dipôle électrique oscillant dans un réseau ionique, donc à une interaction avec un photon.
Cette configuration est générale et permet d'établir une nomenclature des phonons suivant que ceux-ci sont acoustiques (A) ou optiques (O), polarisés longitudinalement (L) ou transversalement (T). D'où les appellations : LA, LO, TA et TO.
La présence de bandes de fréquences interdites est une caractéristiques de la physique des solides.
Le phonon comme pseudo-particule
Pseudo-quantité de mouvement
Il est tentant de considérer un phonon de vecteur d'onde <math>\textstyle \,\vec k\, </math> comme s'il possédait une quantité de mouvement <math>\textstyle \,\hbar\vec k\,</math>, par analogie avec les photons, ou toutes les ondes correspondant à un corpuscule (dualité onde-corpuscule). Ce n'est pas tout à fait correct, car <math>\textstyle \, \hbar\vec k\, </math> n'est pas réellement une quantité de mouvement physique. Il est nommé pseudo- quantité de mouvement ou « impulsion » de vibration. Ceci est dû au fait que <math>\textstyle \, \vec k \,</math> n'est déterminé qu'à un multiple de vecteur constant près, le vecteur <math>\textstyle \, \vec K \,</math> du réseau réciproque. Par exemple les coordonnées normales Q et Π sont définies de telle manière que :
- <math>\mathrm{Q}_\vec{k} \equiv \mathrm{Q}_{\vec{k} + \vec{K}} \quad;\quad \Pi_\vec{k} \equiv \Pi_{\vec{k} + \vec{K}}</math>
Pour un réseau unidimensionnel <math>\textstyle \;K = 2n\pi/a\,,\;n=1,2,3...</math>. Un phonon de nombre d'onde k est donc équivalent à un nombre infini d'autres phonons de la même famille de nombres d'onde k ± 2π/a, k ± 4π/a, etc.
Généralement, on ne considère que les phonons de vecteurs d'onde <math> \vec k </math> de chaque famille possédant le plus « petit » vecteur <math> \vec k </math>. L'ensemble de ces vecteurs définit la première zone de Brillouin. D'autres zones de Brillouin peuvent être définies comme des copies de la première zone, décalées d'un multiple de vecteurs du réseau réciproque.
Thermodynamique
L'ensemble des phonons présents dans un cristal peut être considéré comme un gaz, au même titre que le gaz de photons. Les phonons étant des bosons leur distribution statistique obéit à la statistique de Bose-Einstein :
- <math>n_{k,s} = \frac{1}{\exp\left(\dfrac{\hbar\omega_{k,s}}{k_\mathrm{B}T}\right) - 1}</math>
L'énergie du milieu est donc :
- <math>E=\sum_k \sum_{s=1}^3 \left(n_{k,s}+\tfrac12\right)\hbar \omega_k </math>
À partir de la fonction de partition <math>\textstyle \,Z=\Pi_k\Pi_{s=1}^3n_{k,s}e^{-\frac{\hbar\omega_{k,s}}{2k_\mathrm{B}T}}\,</math> la physique statistique permet de calculer toutes les propriétés thermodynamiques.
Interactions des phonons
Les phonons peuvent interagir<ref name=Kittel>Modèle:Ouvrage</ref>,<ref>Modèle:Ouvrage</ref>,<ref name=Klemens/> avec une singularité du réseau cristallin (atome étranger ou isotope, lacune, dislocation ou autre défaut). Il en résulte une diffusion qui ne conserve pas la quantité de mouvement.
Ils peuvent aussi réagir avec un autre phonon ; deux cas sont possibles :
- la somme des deux vecteurs d'onde est telle que leur somme reste dans la zone de Brillouin et le vecteur d'onde résultant a cette somme ;
- dans le cas contraire on a affaire à un processus umklapp (umklappprozeß) qui ne conserve pas la quantité de mouvement pour le phonon résultant, l'excédent étant transféré au réseau dans son ensemble. Il faut noter que ce phénomène anharmonique est lié au terme d'ordre 3 dans le développement du potentiel interatomique du milieu.
Une interaction est aussi possible avec les électrons lorsque ceux-ci sont loin de l'équilibre thermodynamique. Le schéma d'interaction est analogue à l'interaction phonon-phonon et peut conduire à une diffusion umklapp.
Enfin on a déjà mentionné plus haut l'interaction photon-phonon pour la branche optique des cristaux ioniques. Une interaction faible peut également se produire avec les branches acoustiques (diffusion Brillouin).
Tous ces phénomènes conservent l'énergie.
Conduction thermique
Il est possible d'écrire une équation cinétique pour la population <math>\textstyle n_{\mathbf{k},s}</math>, c'est l'équation de Boltzmann-Peierls. Plus générale que l'équation de la chaleur, elle permet le calcul des échanges thermiques au niveau microscopique, par exemple en microélectronique.
La capacité calorifique d'un matériau cristallin est accessible par le modèle de Debye, issu des méthodes décrites ci-dessus.
La conductivité thermique λ qui peut être reliée<ref name=Kittel/> à la chaleur spécifique CV et au libre parcours moyen dans le milieu l par l'expression <math>\textstyle \lambda=\frac{1}{3}C_V l</math> mais qui en pratique est issue de l'expérience.
.gif){kind=link}
