100% ont trouvé ce document utile (1 vote)
1K vues8 pages

TD1 PHONON Solution

Télécharger au format pdf ou txt
Télécharger au format pdf ou txt
Télécharger au format pdf ou txt
Vous êtes sur la page 1/ 8

Département de Physique

Master : Physique Appliquée et Ingénierie Physique


Année Universitaire 2019/2020

TD : PHYSIQUE DU SOLIDE : PHONONS


Série N° 1
Exercice 1
On se propose d’étudier les modes de vibrations (phonons) d’une chaîne linéaire infinie dans
l’approximation harmonique. On considère une chaine d’atomes 1D de longueur totale L
transportant N atomes identiques de masse m. Les atomes sont tous séparés à l’état d’équilibre
mécanique par la distance a. On admettra que chaque atome n’est soumis qu’à l’interaction de
ses 2 proches voisins par des forces de rappel élastiques de constante de raideur C. On suppose
les conditions aux limites périodiques de type Born-Von Karman et on veut montrer qu'il existe
des ondes élastiques longitudinales de pulsation 𝛚 et de vecteur d'onde 𝐤 ⃗ = 𝐤𝐮 ⃗⃗⃗⃗⃗𝐱 pouvant se
propager le long de la chaîne infinie d'atomes

1. On appel un le déplacement longitudinal de l’atome d’indice n par apport à sa position


d’équilibre 𝐱 𝟎𝐧 = 𝐧𝐚 . Ecrire l’équation de mouvement de l’atome d’indice n, liant
𝑢𝑛 aux déplacements 𝐮𝐧−𝟏 et 𝐮𝐧+𝟏 de ses proches voisins.

Le PFD appliqué à l'atome numéroté n donne :


𝝏𝟐 𝒖 𝒏
𝒎 = C [𝒖𝒏+𝟏 −𝒖𝒏 ] + C [𝒖𝒏−𝟏 - 𝒖𝒏 ]
𝝏𝒕𝟐
On injecte la solution proposée dans cette équation :

Soit :

2. En notation complexe, le déplacement de la masse numérotée n est:


𝐧
𝒖𝒏 = 𝐀𝐞𝐢(𝛚𝐭−𝐤𝐱𝟎 )
a. Pourquoi A ne dépend-il pas de n ?
L'onde n'est ni stationnaire ni amortie : tous les atomes ont donc même amplitude pour une
onde progressive. Le système est globalement invariant par translation.
Ainsi, A ne dépend pas de l'atome considéré.
b. Quelle contrainte les conditions aux limites imposent-elles à 𝐤?
Les conditions aux limites périodiques impliquent que

𝒖𝒏 = 𝒖𝒏+𝑵 , 𝒏 = 𝟏, … 𝑵
𝒆𝒊𝒌𝒏𝒂 = 𝒆𝒊𝒌(𝒏+𝑵)𝒂 , 𝒏 = 𝟏, … 𝑵
La variable k ne peut donc prendre que les valeurs discrètes :
𝟐𝝅 𝟐𝝅 𝑵 𝑵
𝒌=𝒑 =𝒑 , 𝒑=− … − 𝟏, 𝟎, 𝟏, …
𝑵𝒂 𝑳 𝟐 𝟐

1
Département de Physique
Master : Physique Appliquée et Ingénierie Physique
Année Universitaire 2019/2020

3. Déterminer, la relation de dispersion 𝛚 = 𝐟(𝐤) donnant les valeurs possibles de la


pulsation 𝛚 des ondes susceptibles de se propager sur la chaîne en fonction du vecteur
d’onde 𝐤

On injecte la solution proposée dans l’équation de mouvement:


−𝐦𝝎𝟐 = 𝐂(𝐞𝐢𝐤𝐚 + 𝐞−𝐢𝐤𝐚 − 𝟐)
Soit :
𝑪 𝟏
𝝎(𝒌) = 𝟐√ |𝒔𝒊𝒏 ( 𝒌𝒂)|
𝒎 𝟐

𝛑 𝛑
4. Tracer cette relation dans l’intervalle[− 𝐚 , 𝐚 ]. Justifier le choix de cet intervalle

5. Que déduit-on de la non linéarité de cette relation ?

Cette relation de dispersion est non linéaire : il y a dispersion.


Les phonons ont une structure de paquets d’ondes

6. Donner la pulsation maximale 𝛚𝐦𝐚𝐱 des ondes pouvant se propager le long de la


chaine. Comment vibrent deux atomes voisins.

𝑪
La fréquence maximale est 𝝎𝒎𝒂𝒙 = 𝟐√𝑴 et correspond à K = ∓π / a

𝑢𝑛+1
= 𝑒𝑥𝑝{𝑖𝐾𝑎} = 𝑒𝑥𝑝{𝑖𝜋} = 𝑐𝑜𝑠 𝜋 + 𝑖 𝑠𝑖𝑛 𝜋 = −1
𝑢𝑛

les atomes voisins sont déphasés de 180 degrés

7. Calculer les vitesses de phase 𝐕𝐩 (𝐤) et de groupe 𝐕𝐠 (𝐤).

La vitesse de phase est : (on se place dans la suite dans le cas où )

2
Département de Physique
Master : Physique Appliquée et Ingénierie Physique
Année Universitaire 2019/2020

𝟏
𝝎 𝝀 𝑪 |𝒔𝒊𝒏 (𝟐 𝒌𝒂)|
𝒗𝒑 = = = 𝟐 √
𝒌 𝑻 𝑴 𝒌
Avec T la périodicité et λ la longueur d’onde du mouvement oscillatoire

La vitesse de groupe :
𝒅𝝎 𝝎𝒎𝒂𝒙 𝒂 𝒌𝒂
𝒗𝒈 = = 𝒄𝒐𝒔 ( )
𝒅𝒌 𝟐 𝟐

8. Dans le cas des pulsations faibles 𝛚 ≪ 𝛚𝐦𝐚𝐱 , quelle est l’expression de la vitesse de
propagation du son ? comparer la longueur d’onde 𝛌 et a et indiquer comment vibrent
deux atomes voisins ?

Pour les pulsations faibles 𝝎 ≪ 𝝎𝒎𝒂𝒙 , k → 0(grandes longueurs d’onde)

ω C
vg = ≈ a√ = S = const
k M

S est la vitesse du son. La vitesse du son ne dépend pas de la fréquence.


𝟐𝜋 𝟐𝝅
𝒌= =𝒑 avec p<<N
𝝀 𝑵𝒂

𝑵𝒂
𝝀= >> a
𝒑

𝒖𝒏+𝟏
= 𝒆𝒙𝒑{𝒊𝑲𝒂}
𝒖𝒏

, k→0
𝒖𝒏+𝟏
= 𝒆𝒙𝒑{𝒊𝑲𝒂} = 𝟏
𝒖𝒏

les atomes voisins sont en phase

9. Pour les grandes longueurs d’onde (𝒂 ≪ 𝝀) ), on peut considérer que le cristal est un
milieu continu (modèle macroscopique) et on définit une fonction 𝒖(𝒙, 𝒕) = 𝒖𝒏 (𝒕)
telle que :
𝒖𝒏 (𝒕) = 𝒖(𝒙 = 𝒏𝒂, 𝒕)
1. A partir d’un développement de Taylor, montrer que :
𝝏𝒖 𝝏𝟐 𝒖 𝒂𝟐
𝒖𝒏+𝟏 (𝒕) − 𝒖𝒏 (𝒕) = 𝒖(𝒙 + 𝒂, 𝒕) − 𝒖(𝒙, 𝒕) ≈ ( ) 𝒂 + ( 𝟐 )
𝝏𝒙 𝝏𝒙 𝟐
Regarder que :
𝜕𝑢 𝜕 2 𝑢 𝑎2
𝑢(𝑥 + 𝑎, 𝑡) = 𝑢(𝑥, 𝑡) + ( ) 𝑎 + ( 2 )
𝜕𝑥 𝜕𝑥 2
2. En déduire l’équation aux dérivés partielles :

3
Département de Physique
Master : Physique Appliquée et Ingénierie Physique
Année Universitaire 2019/2020

𝝏𝟐 𝒖 𝑚 𝝏𝟐 𝒖

𝝏𝒙𝟐 𝑐𝑎2 𝝏𝒕𝟐

Remplacer dans l’équation de mouvement :(question 1)


𝜕 2 𝑢𝑛
𝑚 = C [𝑢𝑛+1 −𝑢𝑛 ] + C [𝑢𝑛−1 - 𝑢𝑛 ]
𝜕𝑡 2

𝜕 2𝑢 𝜕𝑢 𝜕 2 𝑢 𝑎2 𝜕𝑢 𝜕 2 𝑢 𝑎2
𝑚 2 = C [ ( ) 𝑎 + ( 2 ) ] + C [- ( ) 𝑎 + ( 2 ) ]
𝜕𝑡 𝜕𝑥 𝜕𝑥 2 𝜕𝑥 𝜕𝑥 2

𝜕 2𝑢 𝜕 2𝑢 2
𝑚 = C [ ( )𝑎 ]
𝜕𝑡 2 𝜕𝑥 2
D’où
𝜕 2𝑢 𝑚 𝜕 2𝑢

𝜕𝑥 2 𝑐𝑎2 𝜕𝑡 2

3. Quelle est la signification physique de cette équation ? Que représente la


𝑐𝑎2
quantité 𝑚

Cette équation aux dérivées partielles hyperbolique est caractéristique des


𝑐𝑎2
équations de propagation de type D'Alembert. La quantité √ 𝑚 représente la
vitesse de phase.

Exercice 2 : Vibrations d’une chaine atomique


On considère une chaine de 2N atomes ayant une maille élémentaire de longueur a, composé
de deux atomes séparés d’une distance a/2. Le premier atome de la maille est de masse m1 et
le second atome de masse m2 > m1. On note 𝑥𝑛 = 𝑛 𝑎 la position de la maille numéro n, avec
n 2 J0 : N−1K et on suppose des conditions aux limites périodiques. On note 𝑈𝑛 (𝑡) et 𝑉𝑛 (𝑡)
les déplacements longitudinaux des atomes de masse m1 et m2 de la maille numéro n par
rapport à leur position d’équilibre, et on note C la constante élastique entre deux atomes voisins
(voir figure ci-dessous, les positions d’équilibre des atomes sont représentés par des lignes
verticales).

1. Montrer que les équations du mouvement classiques des atomes de la maille numéro n

𝑑2 𝑈𝑛 (𝑡)
𝑚1 = −2𝐶𝑈𝑛 (𝑡) + 𝐶(𝑉𝑛 (𝑡) + 𝑉𝑛−1 (𝑡))
𝑑𝑡 2
𝑑2 𝑉𝑛 (𝑡)
(𝑡) + 𝑈𝑛+1 (𝑡)) − 2𝐶𝑉𝑛 (𝑡)
{ 𝑚2 𝑑𝑡 2 = 𝐶(𝑈𝑛

4
Département de Physique
Master : Physique Appliquée et Ingénierie Physique
Année Universitaire 2019/2020

D’après la loi de Hooke la force appliquée est proportionnelle au déplacements relatif s des
atomes en interaction

𝑑2 𝑈𝑛 (𝑡)
𝑚1 = 𝐶(𝑉𝑛 (𝑡) − 𝑈𝑛 (𝑡)) − 𝐶(𝑈𝑛 (𝑡) − 𝑉𝑛−1 (𝑡))
𝑑𝑡 2
𝑑 2 𝑉𝑛 (𝑡)
(𝑡) − 𝑉𝑛 (𝑡)) − 𝐶(𝑉𝑛 (𝑡) − 𝑈𝑛+1 (𝑡))
{ 2 𝑑𝑡 2 = 𝐶(𝑈𝑛
𝑚

D’où
𝑑2 𝑈𝑛 (𝑡)
𝑚1 = −2𝐶𝑈𝑛 (𝑡) + 𝐶(𝑉𝑛 (𝑡) + 𝑉𝑛−1 (𝑡))
𝑑𝑡 2
𝑑2 𝑉𝑛 (𝑡)
(𝑡) + 𝑈𝑛+1 (𝑡)) − 2𝐶𝑉𝑛 (𝑡)
{ 2 𝑑𝑡 2 = 𝐶(𝑈𝑛
𝑚

2. On cherche des solutions de la forme 𝑈𝑛 (𝑡) = 𝑈0 𝑒 𝑖(𝑘𝑥𝑛−𝜔𝑡) (t) et 𝑉𝑛 (𝑡) = 𝑉0 𝑒 𝑖(𝑘𝑥𝑛−𝜔𝑡) ,


𝑈0 et 𝑉0 étant des constantes a priori complexes. Comment s’interprètent les constantes
k et ω. Rappeler leurs unités.
K : vecteur d’onde en m-1
ω : pulsation en rad.s-1
D’autres unités sont usuelles

En insérant ces solutions dans les équations du mouvement, montrer que les équations du
mouvement peuvent se réécrire sous la forme :
𝐴 𝐵 𝑈0 0
( )( ) = ( )
𝐶 𝐷 𝑉0 0
−𝐦𝟏𝝎 𝑈0 = − 𝑪( 𝟐𝑈0 − 𝑉0 (𝟏 + 𝒆−𝒊𝒌𝒂 ) )
𝟐

−𝐦𝟐𝝎𝟐 𝑉0 = − 𝑪( − 𝑈0 (𝟏 + 𝒆𝒊𝒌𝒂 ) + 𝟐𝑉0 )

−𝐦𝟏𝝎𝟐 + 𝟐𝑪 𝑪(𝟏 + 𝒆−𝒊𝒌𝒂 ) 𝑈0 0


( ) ( ) = ( )
𝑪(𝟏 + 𝒆𝒊𝒌𝒂 ) −𝐦𝟏𝝎𝟐 + 𝟐𝑪 𝑉0 0
Où on calculera les coefficients A, B, C et D en fonction de ω, de k et des données du problème.
4. Pour avoir une solution non nulle, il s’ensuit que le déterminant de cette matrice doit être
égal à 0. Montrer que cette contrainte donne la relation de dispersion du système suivante

𝜔02 𝑘𝑎
𝜔2 (𝑘) = (1 ± √1 − 𝛾 2 𝑠𝑖𝑛2 ( ))
2 2

𝑚1 +𝑚2 𝑚1 𝑚2
Avec : 𝜔02 = 2𝐶 et 𝛾 2 = 4 (𝑚 2
𝑚1 𝑚2 1 +𝑚2 )

5
Département de Physique
Master : Physique Appliquée et Ingénierie Physique
Année Universitaire 2019/2020

On obtiendra une solution pourvu que le déterminant de la matrice précédente est nul :
−𝐦𝟏𝝎𝟐 + 𝟐𝑪 𝑪(𝟏 + 𝒆−𝒊𝒌𝒂 )
| |=0
𝑪(𝟏 + 𝒆𝒊𝒌𝒂 ) −𝐦𝟐𝝎𝟐 + 𝟐𝑪

(−𝐦𝟏𝝎𝟐 + 𝟐𝑪)(−𝐦𝟐𝝎𝟐 + 𝟐𝑪) − 𝑪𝟐 (1+ 𝒆−𝒊𝒌𝒂 )( 𝟏 + 𝒆𝒊𝒌𝒂 ) = 𝟎


(𝐦𝟏𝐦𝟐𝝎𝟒 −𝟐𝑪𝝎𝟐 (𝒎𝟏 + 𝒎𝟐) + 𝟒𝑪𝟐 − 𝟐𝑪𝟐 (𝟏 + 𝒄𝒐𝒔(𝒌𝒂)) = 𝟎
𝒌𝒂
(𝐦𝟏𝐦𝟐𝝎𝟒 −𝟐𝑪𝝎𝟐 (𝒎𝟏 + 𝒎𝟐) + 𝟒𝑪𝟐 𝒔𝒊𝒏𝟐 ( ) = 𝟎
𝟐
D’où
𝒌𝒂
𝐶(𝑚1 + 𝑚2) ± √𝑪𝟐(𝑚1 + 𝑚2)2 − 4𝑚1𝑚1𝑪𝟐 𝒔𝒊𝒏𝟐 ( 𝟐 )
𝜔2 (𝑘) =
𝑚1𝑚2

2𝐶(𝑚1 + 𝑚2) 4𝑚1𝑚1 𝒌𝒂


𝜔2 (𝑘) = (1 ± √1 − 𝒔𝒊𝒏𝟐 ( )
𝑚1𝑚2 (𝑚1 + 𝑚2)2 𝟐

D’où
𝜔02 𝑘𝑎
𝜔2 (𝑘) = (1 ± √1 − 𝛾 2 𝑠𝑖𝑛2 ( ))
2 2
𝑚1 +𝑚2 𝑚1 𝑚2
Avec : 𝜔02 = 2𝐶 et 𝛾 2 = 4 (𝑚 2
𝑚1 𝑚2 1 +𝑚2 )

Comment appelle-t-on la solution correspondant au signe +?


Brache optique
Et celle correspondant au signe -.
Branche acoustique
Représenter schématiquement l’allure des deux branches dans la première zone de Brillouin.

Exercice 3: Anomalie de Kohn

6
Département de Physique
Master : Physique Appliquée et Ingénierie Physique
Année Universitaire 2019/2020

Soit un matériau unidimensionnel constitué d’atomes identiques équidistants de a. Il s’agit d’un


métal et l’on suppose que, via les électrons de conduction, la constante de rappel qu’exerce
l’atome n+p sur l’atome en position n est de la forme :
𝐬𝐢𝐧 (𝒑𝒌𝟎 𝒂)
𝑪𝒑 = 𝑨
𝒑𝒂
Où A et k0 sont des constantes et p prend toutes les valeurs entières.
Etablir sous forme d’une série illimitée, la relation de dispersion des phonons 𝝎 = 𝒇(𝒌). Sans
𝝏(𝝎𝟐 )⁄
chercher à effectuer cette sommation, montrer que la dérivée 𝝏𝒌 devient infinie pour
une valeur particulière de k que l’on déterminera. Donner l’allure de la courbe de dispersion.

Solution :
D’après le cours on a :
𝑁

−𝑀𝜔 = ∑ 𝛷|𝑛′ −𝑛| (𝑒 𝑖𝑘(𝑛 −𝑛)𝑎 − 1)
2

𝑛′ =1
2
𝜔2 = − ∑ 𝐶𝑝 (1 − 𝑐𝑜𝑠 (𝑘𝑝𝑎))
𝑀
𝑃>0
Substituons Cp par son expression :
2𝐴 𝑠𝑖𝑛 (𝑝𝑘0 𝑎)
𝜔2 = ∑ (1 − 𝑐𝑜𝑠 (𝑘𝑝𝑎))
𝑀 𝑝𝑎
𝑃>0
Soit :
𝜕𝜔2 2𝐴
= ∑ 𝑠𝑖𝑛(𝑝𝑘0 𝑎) 𝑠𝑖𝑛(𝑘𝑝𝑎)
𝜕𝑘 𝑀
𝑃>0
𝐴
= ∑[𝑐𝑜𝑠(𝑝(𝑘0 − 𝑘)𝑎) − 𝑐𝑜𝑠(𝑝(𝑘0 + 𝑘)𝑎)]
𝑀
𝑃>0

Quand k=+-k0 l’un des deux cosinus devient alors unitaire :

𝜕𝜔2 𝐴
𝑙𝑖𝑚 = 𝑙𝑖𝑚 ∑[1 − 𝑐𝑜𝑠(2𝑝𝑘0 𝑎)] = + ∞
𝑘→𝑘0 𝜕𝑘 𝑘→𝑘0 𝑀
𝑃>0
2
𝜕𝜔 𝐴
𝑙𝑖𝑚 = 𝑙𝑖𝑚 ∑[𝑐𝑜𝑠(2𝑝𝑘0 𝑎) − 1] = − ∞
𝑘→−𝑘0 𝜕𝑘 𝑘→𝑘0 𝑀
𝑃>0
La courbe de dispersion présente un point de rebroussement pour 𝑘 ⃗⃗ = ⃗⃗⃗⃗⃗
𝑘0 ± 𝐺⃗
Son allure est représenté sur la figure suivante :
En réalité la constante k0 est telle que : ⃗⃗⃗⃗⃗ ⃗⃗⃗⃗⃗
𝑘0 = 2𝑘 ⃗⃗⃗⃗⃗
𝐹 . (𝑘𝐹 est le vecteur de Fermi)

7
Département de Physique
Master : Physique Appliquée et Ingénierie Physique
Année Universitaire 2019/2020

Anomalie de Kohn dans la dispersion des phonons dans chaque dimension

Cela peut être compris en notant que la capacité des électrons à filtrer les ions change
soudainement en raison du changement du nombre de paires électron-trou qui peuvent être
générées en dessous et au-dessus k = 2kF

Vous aimerez peut-être aussi