UNIVERSITE IBN ZOHR Numéro d’examen : …………………………….............

FACULTÉ DES SCIENCES

DÉPARTEMENT DE PHYSIQUE Nom et Prénom :………………………………………..
AGADIR
C.N.E :…………………………………………..
ÉPREUVE DE PHYSIQUE DES Note :
MATERIAUX (2023-2024)
–SMP5–Durée : 1 heure 30minutes

I- Structure diamant (4pt):

On considère la structure cubique diamant dans laquelle cristallisent le carbone (sous haute et haute
température) ainsi que les matériaux semi-conducteurs tels que Si, Ge. Cette structure est représentée
ci-dessous. Soit a le paramètre de maille.
1) Donner le nombre d’atomes par maille et indiquer leur position. (1pt)
8 atomes dans une maille.
𝟏 𝟏 𝟏 𝟏 𝟏 𝟏 𝟏 𝟏 𝟏 𝟏 𝟑 𝟑 𝟑 𝟏 𝟑 𝟑 𝟑 𝟏
(𝟎, 𝟎, 𝟎), (𝟎, , ) , ( , 𝟎, ) , ( , , 𝟎) , ( , , ) , ( , , ) , ( , , ) , ( , , )
𝟐 𝟐 𝟐 𝟐 𝟐 𝟐 𝟒 𝟒 𝟒 𝟒 𝟒 𝟒 𝟒 𝟒 𝟒 𝟒 𝟒 𝟒
Le fait de les lister dans cet ordre met en évidence que la structure diamant peut être décrite comme
𝑎1 + ⃗⃗⃗⃗⃗
deux réseaux cubiques faces centrées décalés de 1/4 (⃗⃗⃗⃗⃗ 𝑎2 +𝑎⃗⃗⃗⃗3 ).
2) Identifier le réseau de Bravais et le motif. (1pt)
Réseau de Bravais : cubique à faces centrées.
Motif : Atome du C en (0, 0, 0) et (1/4, 1/4, 1/4).

3) Donner les expressions des vecteurs de bases de la maille primitive ⃗⃗⃗⃗

𝑎1 , ⃗⃗⃗⃗
𝑎2 et ⃗⃗⃗⃗ ⃗)
𝑎3 dans la base (𝑖, 𝑗, 𝑘
de la maille cubique. (0.75pt)
𝑎 𝑎 𝑎
Vecteurs primitifs du cfc: 𝑎1 = (𝑗 + 𝑘⃗ ), 𝑎2 = (𝑖 + 𝑘⃗ ) 𝑒𝑡 𝑎3 = (𝑖 + 𝑗)
2 2 2

4) Etablir le facteur de structure Fhkl. En déduire les conditions d'extinctions. (1.25pt)

Le facteur de structure est égal à : 𝐹ℎ𝑙𝑘 = (𝐹ℎ𝑘𝑙 )𝑐𝑓𝑐 . (𝐹ℎ𝑘𝑙 )𝑚𝑜𝑡𝑖𝑓
ℎ 𝑘 ℎ 𝑙 𝑘 𝑙 ℎ 𝑘 𝑙
𝐹ℎ𝑘𝑙 = ( 𝑒 𝑖2𝜋(ℎ.0+𝑘,0+𝑙.0) + 𝑒 𝑖2𝜋(2 +2+𝑙.0) + 𝑒 𝑖2𝜋(2 +𝑘.0+2) + 𝑒 𝑖2𝜋(ℎ.0+2+2) ) 𝑓 (1 + 𝑒 𝑖2𝜋(4+4+4) )
(ℎ+𝑘+𝑙)
⇒ 𝐹ℎ𝑘𝑙 = 𝑓(1 + (−1)(ℎ+𝑘) + (−1)(𝑘+𝑙) + (−1) (ℎ+𝑙) )(1 + (−1) 2 )
𝐹ℎ𝑘𝑙 = 0 Lorsque h, k et l ne sont pas de même parité (extinctions dues au réseau) ou bien ℎ + 𝑘 + 𝑙 =
4𝑛 + 2 (extin. dues au motif).

II- Energie de cohésion (4pt):

L’énergie de cohésion pour un solide peut être estimée à partir de l’expression
𝑨 𝜶𝒒𝟐
suivante : 𝑼(𝑹) = 𝑵( 𝒏
− ) , où N est le nombre de paires d’ions, R la distance au plus
𝑹 𝟒𝝅𝝐𝟎 𝑹

proche voisin, q la charge ionique, =1,75 la constante de Madelung et A et n des constantes

déterminées expérimentalement.
1) Quel type de liaison ce potentiel permet-t-il de modéliser ? Nommer deux caractéristiques de ce type de
cristaux. (1pt)

1
Ce potentiel permet de modéliser les liaisons 1. Forte énergie de liaison (ou de cohésion)
ioniques. 2. Haute température de fusion
Les cristaux ioniques sont caractérisés par les 3. Faible conductivité électrique
propriétés suivantes : 4. Faible conductivité thermique
2) Expliquer brièvement l’origine des deux termes dans l’expression de l’énergie. (1pt)
Les particules dans les cristaux ioniques sont des ions électriquement chargés.
Le potentiel d’interaction peut être écrit comme la somme d'un potentiel répulsif de champ central de la
𝐴 𝛼𝑞2
forme( ). et un potentiel de Coulomb attractive (− 4𝜋𝜖 𝑅).
𝑅𝑛 0
3) Trouvez la distance interatomique à l’équilibre 𝑅0 . (1pt)
𝑑𝑈0 𝛼𝑞2 𝐴𝑛 𝛼𝑞 2
À l’équilibre : ( ) = 0 ⇒ 𝑑𝑈0
= −𝑁 𝐴𝑛
𝑛+1 + 𝑁 2 ⇒ 𝑅 𝑛+1 =
𝑑𝑅 𝑅=𝑅0 𝑑𝑅 𝑅 4𝜋𝜀0 𝑅 0 4𝜋𝜀0 𝑅02
𝟒𝝅𝜺𝟎 𝑨
⇒ 𝑹𝒏−𝟏
𝟎 = 𝒏
𝜶𝒒𝟐
𝜶𝑵𝒒𝟐 𝟏
4) Montrer que l’énergie de cohésion par atome est donnée par : 𝑼𝟎 =− (𝟏 − ). (1pt)
𝟒𝝅𝝐𝟎𝑹𝟎 𝒏
𝐴 𝛼𝑞 2 𝐴𝛼𝑞 2 𝛼𝑞 2 𝜶𝒒𝟐 𝑵 𝟏
𝑈(𝑅) = 𝑁(𝑅𝑛 − 4𝜋𝜖 𝑅) ⇒ 𝑈(𝑅0 ) = 𝑁 4𝜋𝜀 − 𝑁 4𝜋𝜖 ⇒ 𝑼(𝑹𝟎 ) = − 𝟒𝝅𝝐 (𝟏 − 𝒏)
0 0 𝐴𝑛𝑅0 0 𝑅0 𝟎 𝑹𝟎

III- Dynamique du réseau cristallin (12pt) :

On considère une chaîne de longueur L formée de N atomes identiques de masse m, séparés d’une distance

a (avec L>>a). Les vibrations longitudinales des atomes de cette chaîne seront supposées pouvoir être
décrites dans le cadre de l’approximation harmonique et en ne tenant compte que des interactions entre
premiers voisins, avec une constante de rappel unique  On note en un l’écart de l’atome n par rapport à sa
position d’équilibre. On cherche des solutions de l’équation de mouvement sous la forme d’ondes
harmoniques d’amplitude A, de pulsation  et de vecteur d’onde k : un  Aei ( k na  t )

1) Écrire l’équation du mouvement d’un atome au voisinage de sa position d’équilibre et en déduire la loi de
dispersion (k). (1.5pt)
L’équation du mouvement : 𝑚 . 𝑢̈ 𝑛 = 𝛼(𝑢𝑛+1 + 𝑢𝑛−1 − 2𝑢𝑛 ) 0.5pt
On cherche des solutions de la forme: 𝑢𝑛 = 𝐴𝑒 𝑖(𝑘𝑛𝑎−𝜔𝑡)
𝑑𝑢𝑛 𝑑 2 𝑢𝑛
= −𝑖𝜔𝑢𝑛 ⇒ = −𝜔2 𝑢𝑛 et 𝑢𝑛±1 = 𝑒 ±𝑖𝐾𝑎 𝑢𝑛
𝑑𝑡 𝑑𝑡 2
⇒ 𝑚 . 𝑢̈ 𝑛 = 𝛼(𝑢𝑛+1 + 𝑢𝑛−1 − 2𝑢𝑛 )
⇒ −𝒎𝝎𝟐 𝒖𝒏 = 𝜶(𝒆𝒊𝒌𝒂 + 𝒆−𝒊𝒌𝒂 − 𝟐)𝒖𝒏 = 𝟐𝜶(𝒄𝒐𝒔𝒌𝒂 − 𝟏)𝒖𝒏 0.5pt
0.5pt
𝟐
𝒌𝒂 𝟐
𝜶 𝒌𝒂
𝒎𝝎 = 𝟒𝜶𝒔𝒊𝒏 ⇒ 𝐋𝐚 𝐥𝐨𝐢 𝐝𝐞 𝐝𝐢𝐬𝐩𝐞𝐫𝐬𝐢𝐨𝐧 𝝎 = 𝟐 √ |𝒔𝒊𝒏 |
𝟐 𝒎 𝟐

2
2) Tracer les courbes de dispersion (k), on précisera les valeurs de (k) au centre et aux bords de la 1ère
zone de Brillouin. Pourquoi peut-on se contenter de tracer ces courbes de dispersion à l’intérieur de la 1 ère
zone de Brillouin pour avoir toutes les valeurs possibles de (k) ? (2pt)
𝟐𝝅
𝜔(𝑘) = 𝜔(𝑘 ± 𝑛 ) 0.5pt
𝒂
de façon générale on peut toujours se ramener à la 1ere zone de Brillouin
en translatant le vecteur d'onde k d'un vecteur du réseau réciproque, sans modifier la
pulsation.
Pour k << /a : 𝝎 = √𝒎 𝒂. 𝒌.
𝛼
0.5pt
C'est le même comportement que pour une onde élastique se
propageant dans un milieu continu (continuum). La vitesse de phase et la vitesse de
groupe sont égales à vs.
𝛼 𝛼
𝝎 = √𝒎 𝒂. 𝒌 = 𝒗𝒔 𝑲 avec 𝒗𝒔 = √𝒎 𝒂.
𝝅 𝛼
Pour 𝑘 = ± 𝒂 on a : 𝝎 = 𝟐√𝒎. 0.5pt

0.5pt

Courbe de dispersion pour une chaine linéaire. La ligne en trait interrompus correspondrait à un
milieu continu 𝝎 = 𝒗𝒔 𝒌.
3) En utilisant les conditions aux limites cycliques de Born et Von Karman, montrer que le vecteur d’onde k ne
peuvent prendre que des valeurs discrètes qu’on précisera. (1pt)
Conditions aux limites périodiques de Born et Von Karman sur une chaîne de longueur
totale 𝑳 = 𝑵𝒂 contenant N cellules : 𝒖𝒏 = 𝒖𝒏+𝑵 ⇒ 𝑨𝒆𝒊(𝒌(𝒏+𝑵)𝒂−𝝎𝒕) = 𝑨𝒆𝒊(𝒌𝒏𝒂−𝝎𝒕) 0.5pt
𝟐𝝅 𝟐𝝅
⇒ 𝒆𝒊𝒌𝑵𝒂 = 𝟏 ⇒ 𝒌=𝒑 =𝒑 𝒂𝒗𝒆𝒄 𝒑 ∈ 𝒁
𝑳 𝑵𝒂 0.5pt
𝝅 𝝅
 avec – 𝒂 ≤ 𝒌 ≤ 𝒂

 le nombre de valeurs possibles de k est N (N modes propres).

4) On se place dans le cadre de l’approximation de Debye.

3
 a) Définir de manière précise l’approximation de Debye en général. (1pt)
Dans le cadre de l’approximation de Debye :
-On remplace la branche du spectre de vibration par une droite correspondant à un milieu con-
tinu 𝝎 = 𝒗𝒔 𝒌.
-On remplace la 1ere zone de Brillouin par une zone de rayon kD, choisi de telle sorte qu’elle
contient le même nombre N de vecteurs d'onde possibles que la 1ere zone de Brillouin.
b) Exprimer la densité de modes, g(), et vérifier que, dans ce cas de problème unidimensionnel, cette
densité de modes est, en fait, indépendante de .(1pt)
𝒅𝜴 𝟐𝒅𝒌 𝑳
 Densité de modes : 𝒈(𝒌)𝒅𝒌 = 𝒈(𝝎)𝒅𝝎 =
𝜹𝜴
= 𝟐𝝅 = 𝒅𝒌
⁄𝑳 𝝅
0.5pt
 Le chiffre 2 dans la première expression est pour tenir compte de k et -k qui ont le même 𝝎.
𝑳 𝒅𝒌 𝒅𝒌 𝟏 𝑳 𝑵𝒂
⇒ 𝒈(𝝎) = 𝝅 𝒅𝝎 or 𝝎 = 𝒗𝒔 𝒌 ⇒ 𝒅𝝎 = 𝒗 ⇒ 𝒈(𝝎) = 𝝅𝒗 = 𝝅𝒗
𝒔 𝒔 𝒔 0.5pt
5) En déduire les expressions du vecteur d’onde de Debye, kD, la pulsation de DebyeD, et
la température de Debye TD. (1.5pt)
 Le système dans l'espace réel a une longueur L = Na, donc l’espace occupé par chaque état k dans
𝟐𝝅
l'espace réciproque est .
𝑳

 Il y a N valeurs de k dans la première zone de Brillouin.

 kD doit être tel que la « sphère » (à 1D c'est un segment) de Debye contient N valeurs de k on a
𝟐𝒌𝑫 𝑳 𝑵𝒂 𝝅
donc : 𝑵 = 𝟐𝝅 = 𝟐𝒌𝑫 𝟐𝝅 = 𝟐𝒌𝑫 𝟐𝝅 ⇒ 𝒌𝑫 = 𝒂
𝑳
𝝅
 le vecteur d’onde de Debye : 𝒌𝑫 = 𝒂
𝝅
 la pulsation de Debye : 𝝎𝑫 = 𝒗𝒔 𝒌𝑫 = 𝒗𝒔 𝒂 0.5pt
0.5pt
ℏ𝝎𝑫 ℏ𝒗𝒔 𝝅
 la température de Debye : ℏ𝝎𝑫 = 𝒌𝑩 𝑻𝑫 ⇒ 𝑻𝑫 = = 0.5pt
𝒌𝑩 𝒌𝑩 𝒂

6) Donner sous forme d’intégrale définie l’expression de l’énergie interne due aux vibrations du réseau (gaz
de phonons).(1pt)
𝝎𝑫 𝝎𝑫

𝟏 𝒈(𝝎) 1pt
𝑼=∫ ̅ + )ℏ𝝎𝒅𝝎= ∫
𝒈(𝝎)(𝒏 ℏ𝝎
ℏ𝝎 𝒅𝝎 + 𝑼𝟎
𝟐 𝒆𝒙𝒑 (𝒌 𝑻) − 𝟏
𝑩
𝟎 𝟎
𝝎𝑫
𝑵𝒂 ℏ𝝎𝒅𝝎
= ∫ + 𝑼𝟎
𝝅𝒗𝒔 𝟎 𝒆𝒙𝒑(ℏ𝝎/𝒌𝑩 𝑻) − 𝟏
𝟏
̅=
𝒏 ℏ𝝎
représente le nombre d'occupation moyen des phonons dans le mode normal k, pour
𝒆𝒙𝒑( )−𝟏
𝒌𝑩 𝑻
un système en équilibre à la température T. On l’appelle la fonction de distribution de Bose-Einstein.

4
7) Exprimer la capacité calorifique Cv de cette chaîne d’atomes (on notera ħ la constante réduite de Planck,
kB la constante de Boltzmann et T la température). On exprimera Cv en particulier dans les deux cas : basse
 xdx 2  x 2 e x dx 2
température 𝑇 ≪ 𝑇𝐷 et haute température 𝑇 ≫ 𝑇𝐷 . On donne : 0 
ex 1 6
et 0 (e x  1) 2

3
. (3pt)

ℏ𝝎 ℏ𝝎
A haute température ℏ𝝎 << 𝒌𝑩 𝑻 on peut écrire : 𝒆𝒙𝒑 𝒌 ≈𝟏+𝒌
𝑩𝑻 𝑩𝑻 1.5pt
𝝎𝑫
𝑵𝒂 𝟏
 𝑙’𝑒𝑛𝑒𝑟𝑔𝑖𝑒 𝑖𝑛𝑡𝑒𝑟𝑛𝑒: U = 𝝅𝒗 ∫ ℏ𝝎
ℏ𝝎 𝒅𝝎 + 𝑼𝟎 =
𝒔 𝒆𝒙𝒑( )−𝟏
𝒌𝑩 𝑻
𝟎
𝝎𝑫
𝑵𝒂 𝟏 𝑵𝒂 𝝎𝑫 𝑵𝒂
𝝅𝒗𝒔
∫ ℏ𝝎 ℏ𝝎 𝒅𝝎 + 𝑼𝟎 ⇒ 𝑼 = 𝝅𝒗 𝒌𝑩 𝑻 ∫𝟎 𝒅𝝎 + 𝑼𝟎 = 𝝅𝒗 𝒌𝑩 𝑻 𝝎𝑫 + 𝑼𝟎 ⇒ 𝑼 = 𝑵 𝒌𝑩 𝑻 +
𝒔 𝒔
𝒌𝑩 𝑻
𝟎

𝑼𝟎
𝒅𝑼
 La capacité calorifique : 𝑪𝑽 = ) = 𝑵𝒌𝑩
𝒅𝑻 𝑽

 Cv ne dépend pas de la température, Même résultat que celui obtenu par la loi de Dulong et
Petit.
A basse température ℏ𝝎 ≫ 𝒌𝑩 𝑻 et dans l’approximation de Debye
1.5pt
ℏ𝝎 +
 On pose 𝒙 = 𝒌 lorsque 𝑻 → 𝟎 ⇒ 𝒙 → +∞
𝑩𝑻
𝝎𝑫 ∞
𝑵𝒂 ℏ𝝎𝒅𝝎 𝑵𝒂 (𝒌𝑩 𝑻)𝟐 𝒙𝒅𝒙
 𝑙’𝑒𝑛𝑒𝑟𝑔𝑖𝑒 𝑖𝑛𝑡𝑒𝑟𝑛𝑒 : U = 𝝅𝒗 ∫ + 𝑼𝟎 = 𝝅𝒗 ∫ + 𝑼𝟎
𝒔 𝒆𝒙𝒑(ℏ𝝎/𝒌𝑩 𝑻)−𝟏 𝒔 ℏ 𝒆𝒙 −𝟏
𝟎 𝟎
∞
𝒙𝒅𝒙 𝝅𝟐 𝑵𝒂 (𝒌𝑩 𝑻)𝟐 𝝅𝟐
Or ∫ = ⇒ 𝑼𝒂𝒄 = 𝝅𝒗 + 𝑼𝟎
𝒆𝒙 −𝟏 𝟔 𝒔 ℏ 𝟔
𝟎

 La capacité calorifique
𝒅𝑼 𝑵𝒂 𝒌𝟐𝑩 𝝅 𝑵 𝟐 𝝅𝟐
𝑪𝑽 = ) = 𝟐 𝑻= 𝒌 𝑻
𝒅𝑻 𝑽 𝒗𝒔 ℏ 𝟔 𝝎𝑫 ℏ 𝑩 𝟑
𝝅𝟐 𝑻
ℏ𝝎𝑫 = 𝒌𝑩 𝑻𝑫 ⇒ 𝑪𝑽 = (𝑵𝒌𝑩 )
𝟑 𝑻𝑫