�������������������������
Description des bandes d' énergie pour un puits de potentiel périodique 1 D
GLL - Avril 2017
(* Notations :
2m
u = α a → variable associée à l'énergie des états; α= E ;
ℏ2
⇒(α×a)2 est proportionnelle à l'énergie, c'est la grandeur recherchée
m×V0
P = ba → grandeur qui caractérise la forme
ℏ2
des puits b: largeur, a: période, V0 :hauteur des puits;
v = k a → variable associée au vecteur d'onde, homogène à une phase rad *)
(*Pour évaluer une cellule graphique, etc.:
sélectionner la cellule par sa bordure latérale à droite ,
puis appuyer sur les touches Shift et Enter,
�����
ou Menu Evaluation >> Evaluate Cells *
���� ������ ��������
����������������������������������������
(*il n'y a de solutions possibles que pour les
énergies = valeurs de u telles que: -1 ≤F[u,P]≤1;
cette condition est matérialisée par les lignes pointillées rouges*)
2 ��� Modèle Kro�[Link]
Sin[u]
������� F[u_, P_] := Cos[u] + P * ;
������� u
Plot[{F[u, 10], F[u, 0], 1, - 1}, {u, 0, 30},
����������������
PlotStyle → {{}, {Dashed, Gray}, {Red, Dotted}, {Red, Dotted}},
�������������� ���� ���� ����� ��������� ����� ���������
PlotLegends → {"cristal P=10", "dans le vide"}, AxesLabel → {"α×a", "F[α×a,P]"}]
����������������� �����������
F[α×a,P]
cristal
������� 1
dans le
α×a
5 10 15 20 25 30
-1
-2
�������� Plot[{F[u, 3], , 1, - 1}, {u, 0, 30}, PlotStyle → {{}, {Red, Dotted}, {Red, Dotted}},
���������������� �������������� ����� ��������� ����� ���������
PlotLegends → {"cristal P=3"}, AxesLabel → {"α×a", "F[α×a,P]"}]
����������������� �����������
F[α×a,P]
2
�������� cristal P=3
1
α×a
5 10 15 20 25 30
-1
������������������������������������������������
ℏ2
Dans le vide, E(k) = 2m
k2
⇒avec les notations choisies, on écrit u2 = (αa)2 =(ka)2
Modèle Kro�[Link] ��� 3
�������������
�������� TicksTable := Table[k * π, {k, - 5, 5, 1}];
�����
(*position des ticks pour l'axe des abscisses*)
Plot0 = Plotv2 , {v, - 5 π, 5 π}, AspectRatio → 1, PlotStyle → {{Black, Dashed}},
���������������� ���������������� �������������� ���� ����
2
Ticks → {TicksTable, All}, AxesLabel → "k×a", "(α×a) "
����������� ���� �����������
2
(α×a)
250
200
150
��������
100
50
k×a
-5 π -4 π -3 π -2 π -π π 2π 3π 4π 5π
4 ��� Modèle Kro�[Link]
�������������������������������������������������������������
�������� Plot1 = Plotv2 , v - 2 π2 , v + 2 π2 , (v - 4 π)2 , (v + 4 π)2 ,
����������������
{v, - π, π}, PlotRange → All, AspectRatio → 1,
������������� ���� ����������������
PlotStyle → {{Black, Dashed}}, AxesLabel → "k×a", "(α×a)2 ",
���� ���� �����������
Ticks → - π, - π 2, 0, π 2, π, Automatic, AxesLabel → "k×a", "(α×a)2 "
����������� �����������
(α×a)2
250
200
150
��������
100
50
π π k×a
-π - π
2 2
����������������������������������������������������
L' équation (1) F[u, P] = Cos[ka] est résolue de manière numérique, en définissant pour chaque
valeur de u le point {k×a,(α×a)2 }, k×a étant solution de l’équation (1) - d’après la résolution de l’éq
de Schrödinger pour un ensemble périodique 1D de puits de potentiel
⇒donc on trace la courbe des points {k×a,u2 } solution de l’équation (1)
�������� PP := 5; (*valeur choisie pour P dans les graphes qui suivent *)
Modèle Kro�[Link] ��� 5
��������������������������������������������
�������� Plot2 = ParametricPlotArcCos[F[u, PP]], u2 , {- ArcCos[F[u, PP]], u ^ 2},
�����������������������⋯ ����������� �����������
{u, - 5 * π, 5 * π}, AspectRatio → 1, PlotStyle -> {Darker[Green]},
���������������� �������������� �������⋯ ����
Ticks -> - π, - π 2, 0, π 2, π, Automatic, AxesLabel → "k×a", "(α×a)2 "
����������� ����������� �����������
(α×a)2
250
200
150
��������
100
50
π π k×a
-π - π
2 2
6 ��� Modèle Kro�[Link]
�������� Show[Plot1, Plot2]
������
(α×a)2
250
200
150
��������
100
50
π π k×a
-π - π
2 2
◼ Comparaison des diagrammes d’énergie pour plusieurs valeurs de P
Modèle Kro�[Link] ��� 7
��������� ShowPlot1, ParametricPlotArcCos[F[u, 2]], u2 ,
������ �����������������������⋯ �����������
{- ArcCos[F[u, 2]], u ^ 2}, ArcCos[F[u, 8]], u2 , {- ArcCos[F[u, 8]], u ^ 2},
����������� ����������� �����������
{u, - 5 * π, 5 * π}, AspectRatio → 1,
����������������
PlotStyle -> {Darker[Green], {Darker[Green]}, Blue, Blue},
�������������� �������⋯ ���� �������⋯ ���� ���� ����
Ticks -> - π, - π 2, 0, π 2, π, Automatic,
����������� �����������
2
AxesLabel → "k×a", "(α×a) ", MaxRecursion → 8
����������� ��������������������������������
(α×a)2
250
200
150
���������
100
50
π π k×a
-π - π
2 2
������������
Contrairement au cas de l’électron dans le vide, le tracé de E(k) dans tout l’espace des k est plus
difficile à expliciter formellement, puisque la fonction ArcCos[x] renvoie toujours des valeurs dans
l’intervalle [-π,π] alors que l’on voudrait sortir de la 1ère zone de Brillouin
⇒il faut décomposer la solution par branche, en repartant du graphe F[u,P] tracé en I
8 ��� Modèle Kro�[Link]
��������� (*Solution [u,P] détermine l'abscisse k×a du graphe,
et permet de sortir de la zone [-π,π] imposée par la fonction ArcCos *)
�����������
Solution[u_, P_] := Which[
����
0 ≤ u ≤ π, ArcCos[F[u, P]],
�����������
π < u ≤ 2 π, - ArcCos[F[u, P]] + 2 π,
�����������
2 π < u ≤ 3 π, ArcCos[F[u, P]] + 2 π,
�����������
3 π < u ≤ 4 π, - ArcCos[F[u, P]] + 4 π,
�����������
4 π < u ≤ 5 π, ArcCos[F[u, P]] + 4 π
�����������
];
(*chaque ligne correspond à une zone
différente des oscillations de la fonction F[u,P] *)
��������� ParametricPlotSolution[u, PP], u2 , - Solution[u, PP], u2 ,
�����������������������������������������������
{u, - 5 * π, 5 * π}, AxesLabel → "k×a", "(α×a)2 ", AspectRatio → 1
����������� ����������������
2
(α×a)
250
200
150
���������
100
50
k×a
-15 -10 -5 5 10 15
Modèle Kro�[Link] ��� 9
��������� Plot3 = ParametricPlotSolution[u, PP], u2 , - Solution[u, PP], u2 ,
�����������������������������������������������
{u, - 5 * π, 5 * π}, AxesLabel → "k×a", "(α×a)2 ", AspectRatio → 1,
����������� ����������������
PlotStyle -> {Darker[Green]}, Ticks → {TicksTable, All}
�������������� �������⋯ ���� ����������� ����
(α×a)2
250
200
150
���������
100
50
k×a
-5 π -4 π -3 π -2 π -π π 2π 3π 4π 5π
10 ��� Modèle Kro�[Link]
��������� Show[Plot0, Plot3]
������
(α×a)2
250
200
150
���������
100
50
k×a
-5 π -4 π -3 π -2 π -π π 2π 3π 4π 5π