UMMTO

Faculté des Sciences

M1 Physique 28 mars 2021
CORRECTION DE L’EXAMEN DE PHYSIQUE STATISTIQUE

Exercice 1 : (13 pts)

Partie 1

1. L’équation de mouvement en 𝑣(𝑥, 𝑡) de chaque particule de masse m soumise à une force de

1 𝜕𝑝(𝑥)
frottement visqueux−𝑚𝛾𝑣 et à la force 𝐹(𝑥) = − s’écrit :
𝑛(𝑥) 𝜕𝑥
𝑑𝑣 1 𝜕𝑝(𝑥) 𝑑𝑣 1 𝜕𝑝(𝑥)
𝑚 𝑑𝑡 = −𝑚𝛾𝑣 − ⟹ + 𝛾𝑣 = − 𝑚 𝑛(𝑥)
𝑛(𝑥) 𝜕𝑥 𝑑𝑡 𝜕𝑥
𝜕𝑝 𝜕𝑛 𝒅𝒗 𝒌𝑻 𝝏𝒏
Or 𝑝 = 𝑛𝐾𝑇 ⇒ 𝜕𝑥 = 𝑘𝑇 𝜕𝑥 ce qui donne au final + 𝜸𝒗 = − 𝒎 𝒏(𝒙) 𝝏𝒙 1pt
𝒅𝒕
𝑑𝑣 𝑑𝑣
2. La solution transitoire est donnée par 𝑑𝑡 + 𝛾𝑣 = 0 ⇒ 𝑣 = −𝛾𝑑𝑡
⇒ 𝑣 = 𝑣0 𝑒 −𝛾𝑡 → 0 lorsque 𝛾𝑡 ≫ 1 autrement dit en régime stationnaire. 0.5 pt
La solution permanente :
𝑘𝑇 𝜕𝑛 𝑘𝑇 𝜕𝑛
⇒ 𝛾𝑣 = − ⇒𝑣=− 1 pt
𝑚 𝑛(𝑥) 𝜕𝑥 𝑚𝛾 𝑛(𝑥) 𝜕𝑥

3. La quantité Γ(𝑥, 𝑡) = 𝑣(𝑥, 𝑡)𝑛(𝑥, 𝑡) est un flux de particules. 0.5 pt

Ce flux est exprimé par la loi phénoménologique de Fick :
𝜕𝑛 𝑘𝑇 𝜕𝑛
Γ(𝑥, 𝑡) = 𝑣(𝑥, 𝑡)𝑛(𝑥, 𝑡) = −𝐷 =−
𝜕𝑥 𝑚𝛾 𝜕𝑥

𝑘𝑇
L’expression du coefficient de diffusion 𝐷 = 𝑚𝛾 = 𝜇𝐾𝑇 1 pt : Equation d’Einstein 0.5 pt

Partie 2
On considère une particule brownienne à 1d, repérée par sa position 𝑥, de masse 𝑚, soumise à une
force de friction −𝑚𝛾𝑣 proportionnelle à sa vitesse 𝑣, ainsi qu’à une force aléatoire 𝐹(𝑡).
1. Les hypothèses sur la force aléatoire de 𝐹(𝑡) de la théorie de Langevin :
< 𝐹(𝑡) > = 0 et < 𝐹(𝑡)𝐹(𝑡 ′ ) > = 𝐶𝛿(𝑡 − 𝑡 ′ ) 1pt
2. Equation de mouvement en 𝑣(𝑡) de la particule.

𝑑𝑣 𝑑𝑣 𝐹(𝑡)
𝑚 𝑑𝑡 + 𝑚𝛾𝑣 = 𝐹(𝑡) ⇒ + 𝛾𝑣 = (1) 1pt
𝑑𝑡 𝑚

Lorsque 𝐹(𝑡) = 0 (solution sans second membre càd solution transitoire) de l’équation (1)
𝑑𝑣
= −𝛾𝑑𝑡 ⇒ 𝑣 = 𝑣(0)𝑒 −𝛾𝑡 = 𝑣0 𝑒 −𝛾𝑡 1pt
𝑣

1/4
A 𝑡 = 0, 𝑣(0) = 𝑣0 .
Lorsque𝐹(𝑡) ≠ 0, la solution particulière de l’équation Langevin
𝑣(𝑡) = 𝐴(𝑡)𝑒 −𝛾𝑡
𝑑𝑣 𝐹(𝑡) 𝑑(𝐴(𝑡)𝑒 −𝛾𝑡 ) 𝐹(𝑡)
Vérifie l’équation (1) : 𝑑𝑡 + 𝛾𝑣 = ⇒ + 𝛾𝐴(𝑡)𝑒 −𝛾𝑡 =
𝑚 𝑑𝑡 𝑚

𝑑𝐴(𝑡) 𝐹(𝑡) 𝑑𝐴(𝑡) 𝐹(𝑡)

⇒ −𝛾𝐴(𝑡)𝑒 −𝛾𝑡 + 𝑒 −𝛾𝑡 + 𝛾𝐴(𝑡)𝑒 −𝛾𝑡 = ⇒ = 𝑒 𝛾𝑡
𝑑𝑡 𝑚 𝑑𝑡 𝑚
1 𝑡 𝛾𝑡 ′ 𝑒 −𝛾𝑡 𝑡 𝛾𝑡 ′
⇒ 𝐴(𝑡) = ∫ 𝑒 𝐹(𝑡 ′ )𝑑𝑡 ′ ⇒ 𝑣(𝑡) = ∫ 𝑒 𝐹(𝑡 ′ )𝑑𝑡 ′
𝑚 0 𝑚 0
3. La solution générale de l’équation de Langevin s’écrit : 𝑆𝑔 = 𝑆𝑡 + 𝑆𝑝 =
1 𝑡 ′
𝑣(𝑡) = 𝑣0 𝑒 −𝛾𝑡 + 𝑚 𝑒 −𝛾𝑡 ∫0 𝐹(𝑡 ′ ) 𝑒 𝛾𝑡 𝑑𝑡 ′ 1 pt

4. La vitesse moyenne < 𝑣(𝑡) >de la particule s’écrit :

1 −𝛾𝑡 𝑡 ′
< 𝑣(𝑡) >=< 𝑣0 𝑒 −𝛾𝑡 > +< 𝑒 ∫ 𝐹(𝑡 ′ ) 𝑒 𝛾𝑡 𝑑𝑡 ′ >
𝑚 0
Seule la force F(t) est régit par une loi de distribution :
𝑡
1 ′
⇒< 𝑣(𝑡) >= 𝑣0 𝑒 −𝛾𝑡 + 𝑒 −𝛾𝑡 ∫ < 𝐹(𝑡 ′ ) > 𝑒 𝛾𝑡 𝑑𝑡 ′
𝑚 0
−𝛾𝑡 ′)
⇒< 𝑣(𝑡) >= 𝑣0 𝑒 car < 𝐹(𝑡 > = 0 1pt
représentation graphiquement : 𝑣(𝑡)
𝑣0

𝑣0 = 0 ⇒< 𝑣 > = 0.
𝑡
5. Calcul de < 𝑣 2 (𝑡) >
1 −𝛾𝑡 𝑡 ′
𝑣(𝑡) = 𝑒 ∫ < 𝐹(𝑡 ′ ) > 𝑒 𝛾𝑡 𝑑𝑡 ′
𝑚 0
Ou
1 −𝛾𝑡 𝑡
𝑣(𝑡) = 𝑒 ∫ < 𝐹(𝑠) > 𝑒 𝛾𝑠 𝑑𝑠
𝑚 0

2 (𝑡)
1 −2𝛾𝑡 𝑡 𝑡
′
𝑣 = 2𝑒 ∫ 𝑑𝑡′ ∫ 𝐹(𝑡 ′ )𝐹(𝑠)𝑒 𝛾(𝑡 +𝑠) 𝑑𝑠
𝑚 0 0
1 𝑡 𝑡
′
⇒< 𝑣2 (𝑡) >= 2 𝑒−2𝛾𝑡 ∫ 𝑑𝑡′ ∫ < 𝐹(𝑡′ )𝐹(𝑠) > 𝑒𝛾(𝑡 +𝑠) 𝑑𝑠
𝑚 0 0
1 𝑡 𝑡 𝐶 𝑡
′
⇒< 𝑣2 (𝑡) >= 2 𝑒−2𝛾𝑡 ∫ 𝑑𝑡′ ∫ 𝐶𝛿(𝑡′ − 𝑠)𝑒𝛾(𝑡 +𝑠) 𝑑𝑠 = 2 𝑒−2𝛾𝑡 ∫ 𝑑𝑡′ 𝑒𝛾2𝑡′
𝑚 0 0 𝑚 0
2( ) 𝐶 −2𝛾𝑡 +2𝛾𝑡 𝐶 −2𝛾𝑡
⇒< 𝑣 𝑡 > = 2𝑒 [𝑒 − 1] = 2 [1 − 𝑒 ] 1.5 pt
2𝛾𝑚 2𝛾𝑚

2/4
 6. On considère que l’équilibre thermodynamique est atteint⇒ 𝛾𝑡 ≫ 1
𝐶
⇒< 𝑣 2 (𝑡) > → 1.5 pt
2𝛾𝑚2

7. Le théorème d’équipartition de l’énergie donne :

1 1 𝐾𝑇 𝐶 𝐾𝑇
𝑚 < 𝑣 2 > = 2 𝐾𝑇 ⇔< 𝑣 2 > = = 2𝛾𝑚2 ⇒ 𝐶 = 2𝑚𝛾𝐾𝑇=2𝑚2 𝛾 2 𝑚𝛾=2𝑚2 𝛾 2 𝐷
2 𝑚

𝑪
⇒ 𝑫 = 𝟐𝒎𝟐 𝜸𝟐 1.5 pt

Exercice 2 (7 pts)

1. Fonction de distribution de Fermi Dirac 𝑓0 à l’équilibre thermodynamique.

1
𝑓0 = 𝜀−𝜀𝑓 1pt
1+𝑒 𝑘𝑇
2. Equation de Boltzmann à l’approximation du temps de relaxation 𝜏(𝜀) correspondant au problème
donné.
𝑞𝐸 𝜕𝑓 𝑓−𝑓0
=− 1pt
𝑚 𝜕𝑣𝑧 𝜏(𝜀)
3. Résolution par approximations successives et en se limitant en termes linéaires dans le champ
électrique E.
𝑞𝐸 𝜕𝑓 𝑓 − 𝑓0 𝑞𝐸 𝜕𝑓
=− ⇒ 𝑓 = 𝑓0 − 𝜏(𝜀) 𝑜𝑟 𝑞 = −𝑒
𝑚 𝜕𝑣𝑧 𝜏(𝜀) 𝑚 𝜕𝑣𝑧
𝑒𝐸 𝜕𝑓0 𝑒𝐸 2 𝜕 2 𝑓
𝑓 = 𝑓0 + 𝜏(𝜀) + [𝜏(𝜀) ]
𝑚 𝜕𝑣𝑧 𝑚 𝜕𝑣𝑧 2
𝑒𝐸 𝜕𝑓
𝑓 = 𝑓0 + 𝜏(𝜀) 𝑚 𝜕𝑣0 en se limitant en termes linéaires dans le champ électrique E. 1pt
𝑧

1
4. 𝜀 = 2 𝑚𝑣𝑧 2 (le mvt est suivant z)

𝜕𝑓0 𝜕𝑓0 𝜕𝜀 𝜕𝑓0

= = 𝑚𝑣𝑧
𝜕𝑣𝑧 𝜕𝜀 𝜕𝑣𝑧 𝜕𝜀
La loi de distribution 𝑓 hors équilibre s’écrit alors
𝜕𝑓0
𝑓 = 𝑓0 + 𝑒𝐸𝜏𝑣𝑧 1pt
𝜕𝜀

5. Le développement limité au premier ordre de𝑓0 donne .

𝜕𝑓0
𝑓 = 𝑓0 (𝜀 + 𝑒𝐸𝜏𝑣𝑧 ) = 𝑓0 + 𝑒𝐸𝜏𝑣𝑧 1pt
𝜕𝜀

3/4
 Cette expression n’est autre que la fonction 𝑓0 dans laquelle l’énergie a été décalée de la
quantité 𝑒𝐸𝜏𝑣𝑧

6. Les électrons sont considérés ici comme des particules quantiques. La dégénérescence est égale à
2𝑒 2𝑒
2 et la densité de courant s’écrit : 𝐽𝑧 = − ℎ3 ∫ 𝑣𝑧 𝑓𝑑 3 𝑝 = − ℎ3 ∫ 𝑣𝑧 𝑓4𝜋𝑝2 𝑑𝑝
2𝑒 ∞ 𝜕𝑓0 ∞ 𝜕𝑓 2𝑒 2 𝐸𝜏
= −4𝜋 ℎ3 ∫0 [𝑓0 + 𝑒𝐸𝜏𝑣𝑧 ] 𝑣𝑧 𝑝2 𝑑𝑝 = −4𝜋
∫0 [𝑣𝑧 2 𝜕𝜀0 ] 𝑝2 𝑑𝑝
𝜕𝜀 ℎ3
3 1
2𝑒 2 𝐸𝜏 ∞ 𝜕𝑓 2𝑒 2 𝐸𝜏 ∞ 𝑣 2 𝜕𝑓 𝑒 2 𝐸𝜏 ∞ 2𝜀 𝜕𝑓
= −4𝜋 ℎ3 ∫0 [𝑣𝑧 2 𝜕𝜀0 ] 𝑝2 𝑑𝑝=−4𝜋 ℎ3 ∫0 [ 3 𝜕𝜀0 ] 𝑝2 𝑑𝑝=−8𝜋 ℎ3 ∫0 [3𝑚 𝜕𝜀0 ] 𝑚2 √2𝜀 2 𝑑𝜀
car ∫ 𝑓0 𝑣𝑧 𝑝2 𝑑𝑝=0
3
3
2 𝑒 2 𝐸 8𝜋√2𝑚2 ∞ 𝜕𝑓0
𝐽𝑧 = − ∫0 𝜏𝜀 2 𝑑𝜀 1pt
3 𝑚 ℎ3 𝜕𝜀

𝜀−𝜀𝑓
𝜕𝑓0 1 𝑒 𝑘𝑇
7. = − 𝑘𝑇 𝜀−𝜀𝑓 2
présente un pic (négatif) très prononcé au voisinage de 𝜀 = 𝜀𝐹
𝜕𝜀
(1+𝑒 𝑘𝑇 )

eF 𝜕𝑓0
Donc peut être remplacée par la -𝛿(𝜀 − 𝜀𝐹 ) ⇒ 𝐽𝑧 =
0 𝜕𝜀
3 3
3 3
2 𝑒2𝐸 8𝜋√2𝑚2 ∞ 2 𝑒2𝐸 8𝜋√2𝑚2
3 𝑚 ℎ3
∫0 𝜏𝜀 2 (𝜀 − 𝜀𝐹 )𝑑𝜀 = 3 𝑚 ℎ3
𝜀𝐹 2 𝜏(𝜀𝐹 ) 0.5 pt

8. Le niveau de Fermi 𝜀𝐹 est très voisin de sa valeur pour T=0 K, 𝜀𝐹0 ⇒

𝟑 𝟑
0 𝟐 𝑒2𝐸 √𝟐𝒎𝟐 𝑒 2𝐸
𝜏(𝜀𝐹 ) = 𝜏(𝜀𝐹 ) et 𝟖𝝅 𝜺𝑭 𝟎 𝟐 𝜏(𝜀𝐹 ) = 𝑛 𝜏(𝜀𝐹 ) = 𝐽𝑧 = 𝜎𝐸 ⇒
𝟑 𝒎 𝒉𝟑 𝒎
𝑛𝑒 2 𝜏(𝜀 𝐹)
𝜎= 0.5 pt
𝒎
f'(e)

-0.3

4/4