0% ont trouvé ce document utile (0 vote)
507 vues10 pages

Modele Kronig Penney

Ce document décrit les bandes d'énergie pour un puits de potentiel périodique 1D. Il présente des graphiques montrant les solutions de l'équation de Schrödinger pour un ensemble périodique 1D de puits de potentiel pour différentes valeurs du paramètre P caractérisant la forme des puits.

Transféré par

DAMEGARIN SAMBIANI
Copyright
© © All Rights Reserved
Nous prenons très au sérieux les droits relatifs au contenu. Si vous pensez qu’il s’agit de votre contenu, signalez une atteinte au droit d’auteur ici.
Formats disponibles
Téléchargez aux formats PDF, TXT ou lisez en ligne sur Scribd
0% ont trouvé ce document utile (0 vote)
507 vues10 pages

Modele Kronig Penney

Ce document décrit les bandes d'énergie pour un puits de potentiel périodique 1D. Il présente des graphiques montrant les solutions de l'équation de Schrödinger pour un ensemble périodique 1D de puits de potentiel pour différentes valeurs du paramètre P caractérisant la forme des puits.

Transféré par

DAMEGARIN SAMBIANI
Copyright
© © All Rights Reserved
Nous prenons très au sérieux les droits relatifs au contenu. Si vous pensez qu’il s’agit de votre contenu, signalez une atteinte au droit d’auteur ici.
Formats disponibles
Téléchargez aux formats PDF, TXT ou lisez en ligne sur Scribd

�������������������������

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 = Plotv2 , {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 = Plotv2 , 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 = ParametricPlotArcCos[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

��������� ShowPlot1, ParametricPlotArcCos[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] *)

��������� ParametricPlotSolution[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 = ParametricPlotSolution[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π

Vous aimerez peut-être aussi