Tese Brhan

République Algérienne Démocratique et Populaire
‫ــ" ا ـ! ـــ وا ـ ـــ ا ـ ــ ــ‬#‫ـ ـ ـ‬$ ‫وزارة ا‬

Ministère de l’Enseignement Supérieur et de la Recherche Scientifique
Université Mohamed El Bachir El Ibrahimi , Bordj Bou Arréridj
« CD:D:EF; ‫ج‬:; » #‫اھ‬:;<‫ا‬:#= ‫ـ!? @ ? > ا‬A
Faculté des Sciences et de la Technologie !#AF FIJ$ ‫م وا‬F ‫@ ا‬# L
Département des Sciences de la Matière ‫م ا !دة‬F E "NO
Mémoire de fin d’étude

PRESENTÉ EN VUE DE L’OBTENTION
DU DIPLOME DE MASTER
Filière : Physique
Option : Physique des matériaux
THÈME
Etude des propriétés électroniques et optiques
des alliages ternaires AlPxSb1-x
Préparé par :Kolli Bourhan Eddine
Soutenu le : 2Juillet 2018
Devant le jury :
Président Dr Grar Nabila

Encadreur Dr Nour El Houda FARES
Examinateur Dr Kbir Hadda
Remerciements
• Je voudrais d’abord remercier « Allah » que DIEU nous bénisse ! Je tiens

à féliciter ma famille pour son soutien moral, et pour l’encouragement qu’elle m’a apporté
tout au long de mes études.
• Je voudrais témoigner de la reconnaissance à tous ceux qui m’ont permis de mener à
terme ce travail.
• Je tiens particulièrement à exprimer ma gratitude à mon encadreur Dr Fares
Nour El Houda pour son apport substantiel dans la réalisation de ce modeste
travail ainsi que pour son ensemble de consultations, tout au long de cette formation, qui
ont inspiré et nourri mon savoir.
• Je remercie également tous les enseignants dont j’ai eu l’honneur d’être l’étudiant au
département science de la matière de l’université de Bordj Bou Arréridj
• Tous mes remerciements à mes amis pour leurs précieuses suggestions.
• Que toutes les institutions qui ont facilité mes recherches soient remerciées pour leur
dévouement et les conseils qu’elles m’ont prodigués.
Table des matières
Remerciements
Dédicaces
Résumés
Sommaire
Liste des tableaux
Liste des figures
Introduction générale
Références bibliographiques
Chapitre I Généralités sur les semi-conducteurs et leurs alliages

I.1. Introduction 5
I.2. Notions générales sur les semi-conducteurs 5
I.2.1. Semi-conducteur intrinsèque 5
I.2.2. Semi-conducteur extrinsèque 6
I.3. Propriétés structurales 6
I.3.1. La constante de réseau 6
I.3.2. Première zone de Brillouin 7
I.4. Propriétés électroniques 9
I.4.1. La structure de bande 9
I.4.2. Gap énergétique
I.4.3. Densité de charge 9

I.5. Propriétés optiques 10
I.6. Propriétés diélectriques 11
I.7. L'approximation du cristal virtuel 11
I.7.1. L’approximation du cristal virtuel améliorée 12
Chapitre II Méthode de calcule des structures de bandes électroniques

II.1. Introduction 17
II.2. L’équation de Schrödinger d’un solide cristallin 17
II.2.1. L’approximation adiabatique 19
II.2.2. Approximation du champ autocohérent 20
II.2.2.1. Approximation de Hartree 20
II.2.2.2. Approximation de Hartree-Fock 22
II.2.2.3. Approximation de Hartree-Fock-Slater 22
II.3. Les méthodes de calcul de la structure de bandes électroniques 23
II.3.1. Méthodes de calcul des électrons libres 23
II.3.2. Méthodes de calcul des électrons presque libres 24
II.3.3. Méthodes de calcul des ondes planes 24
II.3.4. Méthode de calcul des ondes planes augmentées 25
II.3.5. Méthodes de calcul des ondes planes orthogonales 26
II.3.6. Méthodes des pseudopotentiels 26
II.4. Conclusion 26
Chapitre III Méthode des pseudopotentiels

III.1.Introduction 28
III.2. Principe de la méthode des pseudopotentiels 28
III.3. Théorie de Phillips et Kleinman 29
III.4. Modèle de potentiel 32
III.4.1. Modèle de potentiel de Heine-Abarenkov 33
III.4.2. Modèle de potentiel d’Ashcroft 33
III. 5. Méthode empirique du pseudopotentiel (EPM) 34
III.5.1. Approximation Empirique locale 34
III-6-conclusion 46
Chapitre IV Résultats et Discussions

IV.2. Etude des propriétés électroniques 48
IV.2.1. Structure de bandes électronique 50
IV.2.1.2-Gap d'énergie 52
IV.2.3. Densité de charge électronique 56
IV.3. Etude des propriétés optiques 60
IV.4. Etude des propriétés diélectriques 61
IV.4.1. Constante diélectrique de haute fréquence 61
Références bibliographiques 64
Conclusion générale
Liste des figures
Figure 1.1. Maille élémentaire de la structure Zinc. Blende 6

Figure 1.2 Réseau triclinique primitif de l'espace tridimensionnel 7
Figure I.3 Première zone de Brillouin de la structure zinc blende 8
Figure I.4.a. Semi-conducteur à gap direct.b. Semi-conducteur à gap indirect 9
Figure IV.1 : Structure de bande électronique du composé binaireAlSb 50
Figure IV.2 : Structure de bande électronique du composé binaireAlP 51
Figure IV.3 Structure de bande électronique de l’alliage ternaireAlP0,5Sb0,5 51
Γ 52
Figure IV.4 Variation du gap direct Γ d’énergie en fonction de la
composition x avec et sans désordre de l’alliage ternaire
Figure IV.5 Variation du gap d’énergie indirect en fonction de la composition 54
x avec et sans désordre de l’alliage ternaire AlP Sb
Figure IV.6 Variation des gaps d’énergies direct et indirects et et en 55

fonction de la composition x sans désordre de l’alliage ternaire AlP Sb
Figure IV.7 Variation des gaps d’énergies direct et indirects et et en 55
fonction
de la composition x avec désordre de l’alliage ternaire AlP Sb
Figure IV. 8 Densité de charge électronique au point de la somme des quatre 57
bandes de valence selon la direction [111] du composé binaire AlSb
Figure.I V.9 Densité de charge électronique au point de la somme des quatre 58
bandes de valence selon la direction [111] de l’alliage ternaire Al . In . Sb .
Figure IV.10 Densité totale de charges électroniques de la première bande de
conduction selon la direction [1 1 1 ] du composé binaire AlSb
Figure.V.12 La variation de l’indice de réfraction en fonction de la composition x de
l’alliage
ternaire AlPxSb1-x
Figure IV.13 La variation de la constante diélectrique de haute fréquence en
fonction de la composition x de l’alliage ternaire AlPxSb1-x (VCAA)
Liste des tableaux
Tableau.IV.1 Les facteurs de formes symétriques (Vs) et antisymétriques (VA) 48

ajustés en (R) des composés binaires semi-conducteurs AlPetAlSb et leurs
paramètres de réseaux
Tableau.IV.2.Comparaison entre les gaps d’énergies calculés et observés dans les 49
composés Binaires semi-conducteurs AlP et AlSb
Tableau.IV.3. Paramètres de réseau de l’alliage semiconducteur AlPxSb1-x. 49
Tableau.IV.4.indice de réfraction (n) des semi-conducteurs AlSb ,AlPxSb1-x et AlP 59

en utilisant différents modèle
Résumé
Dans ce travail, nous avons fait une étude théorique sur les propriétés électroniques,
optiques et diélectriques de l'alliage semi-conducteur ternaire III-V à base d’Aluminium
cristallisant dans la phase Zinc-blende à savoir : .
La méthode de calcul utilisée pour les états électroniques est celle des pseudos
potentiels. L'effet du désordre compositionnel a été pris en considération par l'introduction
d'un potentiel effectif (VCA amélioré).
Nos résultats sont en très bon accord avec d’autres résultats expérimentaux et
théoriques, ce qui a permis de conclure que la méthode adaptée dans cette étude est valide
dans la description des états électroniques de l’alliage
Abstract
In this work, we have made a theoretical study on the properties of electronic, optical
and dielectric of the III-V ternary semiconductor alloy based on aluminum crystallizing in the
Zinc-blende phase namely:
The calculation method used for the electronic states is the one of the
pseudopotentials. The latter was coupled to the virtual crystal approximation (VCA)in the
treatment of ternary semiconductor alloys. The effect ofcomposition of disorder wastaken into
account byintroducing aneffectivepotential (improved VCA).
Our results are in very good agreement withotherexperimental and theoreticalresults,
whichled to the conclusion that the methodadapted in thisstudyisvalid in the description of the
electronic states of the alloy
!"#$
‫ت‬NOPJL‫ه ا‬RST‫ أ‬VWRSXL YL‫ز‬RKL‫ وا‬Y\WP]L‫ وا‬Y\^‫`_و‬aLb‫اص ا‬PdL‫ ا‬efg Yh_i^ Yj‫_اء درا‬nop RqJr ،IJKL‫ا ا‬D‫ ھ‬FG
F‫ وھ‬،V^sL‫ر ا‬P‫ ط‬FG ‫ر‬PfSù F`L‫م وا‬P\qJLw‫ ا‬efg xJ`Ku F`L‫د ا‬RKpw‫ ا‬Y\yNy
FG FS\{_`L‫\_ ا•~}_اب ا‬y€u D•‫ أ‬xr‫ و‬.YfJ`‚JL‫ت ا‬RaqL‫ ا‬Yzh_‫ ط‬F‫ ھ‬Y\^‫`_و‬aLb‫•ت ا‬R‚fL Yƒxd`XJL‫ب ا‬RX‚L‫ ا‬Yzh_‫ط‬
VCA).shsKu (YLRKG ‫ت‬R^Raƒ‫ل إ‬R•‫† إد‬h_‫‡ ط‬g ‫ن‬RSX‚L‫ا‬
FG R••\\au ‘u F`L‫ ا‬Yzh_}L‫ج أن ا‬R`q`j‫ ا‬eL‫ أدى إ‬RJƒ ، ‫•_ى‬w‫ ا‬Yh_iqL‫ وا‬Y\Sh_‰`L‫ ا‬ŒWR`qL‫• ا‬ƒ x\n IaŽp YzG‫ا‬P`ƒ Rq‰WR`^
VWRSXfL Y\^‫`_و‬aLb‫•ت ا‬R‚L‫’ ا‬O‫ و‬FG Y‚\‚O Yj‫را‬xL‫ه ا‬D‫ھ‬
Introduction
Introduction Générale
L’usage et la connaissance de matériaux aux propriétés particulières a favorisé
l’innovation technologique, qui a introduit dans notre environnement d’innombrables
instruments et outils renfermant des mystères électroniques. Parmi ceux-ci, les
semiconducteurs qui sont des matériaux clés des techniques actuelles, vu la richesse de
leur comportement vis-à-vis des sollicitations auxquelles ils sont soumis [1].
Toutefois ils sont par excellence les matériaux de l’optoélectronique permettant

de contrôler les électrons pour la détection, la modulation et l’émission de la lumière, en
d’autres termes, ils servent à réaliser des composants et systèmes optiques [2].
Parmi les semiconducteurs eux-mêmes, les composés formés des éléments III-V
de la classification de Mendeleïev qui ont des propriétés physiques remarquables entrant
dans la plupart des composants de télécommunication. Par ailleurs les semiconducteurs
IIVI ont longtemps constitué le prototype de matériaux radiatifs, grâce à leur bande
interdite directe qui leur confère une forte efficacité radiative ainsi qu’à leur grande
énergie du gap, conduisant à l’émission de photons situés dans le spectre du visible.
Pour mieux comprendre la structure électronique des semi-conducteurs, il faut

tout d’abord revenir à leur état massif. Leurs propriétés optiques sont liées à leur
structure électronique. La structure électronique d’un semi-conducteur massif présente
une décomposition des énergies accessibles aux électrons sous forme de bandes.
L’étude du comportement des électrons dans un cristal dont le nombre des

électrons est grand es la ésoution de l’équation de Schrödinger est impossible, donc on a
besoin a des méthodes d’approximations.
Les techniques de calcul de la structure électronique mises au point au cours des

dernières décennies sont nombreuses, et en particulier, les méthodes empiriques qui sont
devenues aujourd'hui un outil de base pour le calcul des propriétés électroniques et
structurales des systèmes les plus complexes. Elles sont aussi un outil de choix pour la
prédiction de nouveaux matériaux.
Parmi ces méthodes [3] qui ont été développées pour le calcul de la structure de
bandes on cite la méthode des ondes plans, la méthode des ondes plans augmentées
(A.P.W),la méthode de ondes plans orthogonales (O.P.W) et la méthode pseudo
1
potentiel empiriques (EPM). C'est cette dernière que nous avons utilisé pour étudier
l'alliage ternaire .
L'approche historique le plus simple pour étudier les alliages est l'approximation
du cristal virtuel (VCA), dans lequel l'alliage réel désordonné est représenté par un
solide cristallin d'atomes virtuels qui sont donnés par un moyen d'atomes réels [4], avec
des propriétés qui sont intermédiaires entre celles des constituants purs. Cette théorie
n'est pas " structural" dans le sens qu'elle considère seulement l'occupation moyenne des
sites du réseau, conserve la topologie, mais ignore toutes les informations associées
avec l'arrangement géométrique des atomes autour de chaque site et l'environnement
local des atomes. Cette théorie est donc inadéquate pour l'étude des propriétés "comme
les propriétés électroniques et vibrationnelles " qui dépendent fortement du désordre
substitutionnel, il faut alors développer et améliorer cette approximation. Pour cela on a
introduit l'effet de désordre dans la VCA en ajoutant un potentiel effectif du désordre
dans le potentiel cristallin pour étuder l'alliage ternaire.
Pour aboutir aux objectifs de cette étude, notre thèse acheminera le parcours
suivant :
On débute par une introduction générale, où on expose la problématique de cette

étude, en donnant un bref aperçu sur le choix de la méthode de calcul.
Dans le premier chapitre on expose la théorie des semi-conducteurs et leurs

alliages ainsi que la théorie du cristal virtuel, en tenant compte de l’effet du désordre
compositionnel comme un potentiel effectif.
Tout système électronique est décrit par l’équation de Schrödinger, qui prend en
compte l’énergie cinétique de tous les ions et des électrons et leurs interactions, le
problème est très compliqué vu le grand nombre de particule , donc on est besoin a
des approximations qu’on a exposé au deuxième chapitre.
Une étude bien détaillée, en équations analytiques, où a été mise en évidence de

la méthode des pseudopotentiels couplée à l’approximation du cristal virtuel pour le
calcul des états électroniques de l'alliage ternaire semi-conducteur
cristallisant dans la phase Zinc-blende, est présentée dans le troisième chapitre.
2
Dans le quatrième chapitre on a exposé les principaux résultats obtenus avec une
interprétation de ceux-ci ainsi qu’une comparaison avec d'autres résultats théoriques et
expérimentaux.
Finalement, notre travail est achevé par une conclusion générale.
3
[1] https://crppwww.epfl.ch/avcp/images/AMMAN.PDF.
[2] G. Roosen, « Matériaux semiconducteurs III-V, II-VI et nitrures pour

l’optoélectronique », Ed. Hermes Science (2003).
[3] N. W . Aschroft and N. D. Mermin, "Solid State Physics", Ed Holt, Renehart and
Winston Philadelphia, (1976).
[4]J. Mycielshy and A. M.W. Itowski, Phys.Stat.Sol.(b) 134 (1986) 134
4
Chapitre I
LES SEMICONDUCTEURS ET LEURS ALLIAGES

Chapitre I généralités sur les semi-conducteurs et leurs alliages
I.1. Introduction
Les technologies semi-conducteurs ont connues la plus fort progression industrielle des
vingt dernières années, demeurent un outil essentiel de développement des économies
modernes et cela, grâce à leur propriétés physiques curieuse qui sont connues depuis
longtemps
Parmi ceux- ci, les matériaux semi-conducteurs des groupe III-V et II-VI de la
classification de Mendeleïev, qui ont un champ d’application très élargi.
Les alliages semi-conducteurs sont aussi considérés comme une classe important des
matériaux, grâce à la richesse de leurs structures électroniques et leurs propriétés physiques
qui peuvent être ajustées et contrôlées pour l’obtention d’un composant optoélectronique bien
désiré.
I.2. Notions générales sur les semiconducteurs

Les semi-conducteurs sont des matériaux dont la conductivité électrique intermédiaire
entre celle des conducteurs et des isolants. Cette conductivité varie sous l'effet de la
température, l'éclairement et la présence d'impuretés (dopage, défauts du réseau).
Les semi-conducteurs sont principalement les éléments de la colonne IV du tableau de

classification des éléments (Si, Ge) mais aussi des composés III-V (GaAs, GaN) ou II-VI
(ZnS, CdTe).
Les matériaux semi-conducteurs III-V sont des corps composés formés à partir d'un
élément de la IIIème colonne et d'un élément de la Vème colonne de la classification
périodique de Mendeleïev.
Dans cette thèse nous avons utilisés les alliages binaires des semi-conducteurs III-V:
Antimoniure d'aluminium (AlSb) et Phosphure d'aluminium (AlP).
Il existe deux types de semi-conducteurs, l'un parfait dit intrinsèque et l'autre dopé
appelé extrinsèque.
I.2.1. Semi-conducteur intrinsèque
Un semi-conducteur est dit intrinsèque, si ces propriétés électriques sont entièrement

dues à sa composition chimique et à sa structure cristalline et non pas à la présence
d’impuretés. Ce comportement correspond à son état pur, donc parfait, sans défaut structural
ou impuretés. Toutefois, un semi-conducteurréel n’est jamais parfaitement intrinsèque, ces
5
semi-conducteurs ne conduisent pas ou très peu le courant sauf s’ils sont portés à haute
température [2].
I.2.2. Semi-conducteur extrinsèque
Ces semi-conducteurs sont d’un grand intérêt technique, car leur conductivité
électrique peut être fortement modifiée par l’adjonction dans leurs structures cristallines
d’impuretés spécifiques. Les propriétés électriques étant ainsi déterminées par la présence des
impuretés, on parle de conductivité extrinsèque. L’adjonction délibérée d’impuretés à
l’intérieur d’un semi-conducteur intrinsèque est appelée dopage [1,2].
I-3 Propriétés structurales
L’immense majorité des semi-conducteurs utilisés dans les applications électroniques

ou optoélectroniques sont basés sur la même structure cristalline, celle de diamant (Si , Ge )
ou la structure Zinc blende (GaAs,AlP). Cette structure est présentée sur la figure ci-dessous
[3]
Figure 1.1. Maille élémentaire de la structure Zinc Blende du composé binaire AlSb
I.3.1. La constante de réseau
C’est une grandeur utilisée pour d’écrire la maille d’un cristal, la maille cristalline est
définie par troisparamètres (a, b, c) et par trois angles (α, β, λ) [4].
6
Figure 1.2 Réseau triclinique primitif de l'espace tridimensionnel.
C’est pour cela que la détermination de la constante du réseau est la première étape à
faire pour déterminer la structure d’un cristal. On voit clairement que, la connaissance de cette
constante du réseau permet de calculer la densité des atomes et donc également la densité
électronique.
I.3.2. Première zone de Brillouin
La première zone de Brillouin est définie de manière unique comme la maille

primitive dans l'espace réciproque. La notion de zone de Brillouin est nńécessaire pour
décrire les propriétés vibrationnelles ou électroniques d’un cristal dans lequel la symétrie de
translation joue un rôle essentiel.
La première zone de Brillouin est la cellule de Wigner-Seitz du réseau réciproque, on

peut la construire en traçant les plans bissecteurs des vecteurs joignant à un point quelconque
du réseau réciproque. Cette première zone de Brillouin possède la forme d'un octaèdre
tronqué (Figure.I.3). Cet espace réduit du réseau réciproque est caractérisé par des points et
des lignes d’importance particulière, on les appelle les points et les lignesde haute symétrie
[5].
7
I.3.2. 1. Les points de haute symétrie
: Centre de la première zone de Brillouin avec les coordonnées

: Centre d’une face carrée de l’octaèdre appartenant à l’un des axes , ou avec l’une
des faces carrées. =2 (±1,0,0), =2 (0, ±1,0), =2 (0,0, ±1)
: Centre d’une face hexagonale de l’octaèdre ayant les coordonnées =2 (1,1,1).

: Point W se trouve sur l’un des sommets des faces carrées et de coordonnées =2
(0, 1 2 , 1).
Z : Point est situé sur la ligne qui joint le centre d’une face carrée à l’un des coins de
l’octaèdre et de coordonnées =2 (1, 1 2 , 1).
I.3.2.2. Les lignes de haute symétrie
∆ : Ligne de direction (100), reliant le centre Γ au point X.
Σ : Point appartenant au plan de symétrie kx = ky ou ky = kz ou kx = kz.
Λ : Ligne de direction (100), reliant le centre de la zone (Γ) au centre d’une face hexagonale
qui est le point L de l’octaèdre.
Figure I.3 Première zone de Brillouin de la structure zinc blende avec la représentation des
points et lignes de hautes symétries.
8
I.4. Propriétés électroniques
I.4.1. La structure de bande
La structure de bandes est généralement identique pour les composés (III.V).

C'est.à.dire que le minimum de la bande de conduction et le maximum de la bande de valence
correspondent à la même valeur du vecteur d’onde k dans la structure de bandes. Dans ces
matériaux, la transition entre les bandes de valence et de conduction est le plus souvent
radiative, ce qui leur accorde des propriétés optiques très intéressantes pour la fabrication de
systèmes optoélectroniques.
I.4.2. Gap énergétique
Les courbes Ecv ( ) dites aussi relations de dispersion où est le minimum de la bande
de conduction, le maximum de la bande de valence et le vecteur d’onde associé à un
électron, ces courbes font apparaître deux types de semi-conducteur. Si ce maximum et ce
minimum correspondent à la même valeur de on dit que le semi-conducteur est à gap
direct. Si au contraire, ce maximum et ce minimum correspondent à des valeurs de
différentes, on dit que le semi-conducteur est à gap indirect (Figure I.4)
Figure I.4 a. Semi-conducteur à gap directb. Semi-conducteur à gap indirect
I.4.3. Densité de charge
La densité de charge électronique est un moyen pour l’analyse et l’interprétation d’un

nombre considérable de propriétés des matériaux. En particulier la distribution de la densité
de charge associée à la bande de valence est un aspect important de la structure électronique
puisqu’elle indique la nature des liaisons cristallines. Pour un matériau qui se cristallise dans
9
une structure de type zinc blende la densité décharge sera mieux représentée suivant la
direction <111> qui contient les deux atomes de la maille unitaire. On peut l’obtenir par la
formule suivante :
( )= ∑ , ( )
I.5. Propriétés optiques
La connaissance de l'indice de réfraction (n) des semi-conducteurs est indispensable

pour des dispositifs tels que des cristaux photoniques, des guides d'ondes, des cellules et des
détecteurs d'énergie solaire [7,8]. L’indice ( ) est calculé en utilisant différents modèles et en
tenant compte de l’effet du désordre.
(i) Modèle de Ravindra et Srivastava [9]
n" E$ = 108 I−2
(ii) Modèle de Gupta et Ravindra [10]
n = 4.084 − 0.62E$ I−3
(iii) Modèle de Hervé et Vandamme [11]
45
13.6
n = /1 + 1 2 3 I−4
E$ + 3.4
(iv) Modèle de Reddy et Anjaneyulu [12]
E$ e7 = 36.6 I−5
(v) Modèle de Ravindra [13]
n = 4.16 − 0.85E$ I−6

(vi) Modèle de Reddy et Ahammed [14]
10
n = 4.084 − 0.62E$ I−7
I.6. Propriétés diélectriques
La constante diélectrique appelée aussi permittivité électrique, décrit la réponse d’un

milieu donné à un champ électrique. Elle a des conséquences importantes sur les propriétés
physiques du solide. Elle intervient dans de nombreux domaines, notamment, dans l’étude de
la propagation des ondes électromagnétique, en particulier en optique, dans la détermination
de l’indice de réfraction, la réflexion et l’absorption de la lumière. La constante diélectrique
: (;) et l’indice de réfraction ( ) sont liés par la relation :
:=n I.8
I.7. L'approximation du cristal virtuel
Un semi-conducteur ternaire est supposé uniforme, cependant, l'effet des fluctuations

dû à la composition x est inévitable. Généralement, on considère ces fluctuations comme une
perturbation qui est calculée par le second ordre de l'approche V.C.A [15]. Le cristal ternaire
AB> C4@> est composé des molécules AB avec une fraction molaire x et des molécules AC
avec une fraction (1-x). Ainsi, la propriété physique F(x) peut être représentée comme une
simple interpolation analytique des propriétés de ses composés binaires AC et BC [18].
La valeur moyenne de la propriété physique F étant :
CCCCCC = xFDEF + (1 − x)FDEGH

F(x) I.9
Par exemple la constante du réseau a(x) de l'alliage ternaire AlP> Sb4@> sera donné par la loi
de Vegard [19] :
a(x) ≈ CCCCCC
a(x) = x. aDEN+ (1– x)aDEGH I.10
Où : aDEF est la constante de réseau du composé AlP.
aDEGH est la constante de réseau du composé AlSb.
Par contre, le gap énergétique Eg(x) de l'alliage ternaire est donné par :
11
E$ = xEDEF + (1 − x)EDEGH + bx(1 − x) I.11
Avec :EDEF est le gap du composé AlP.
EDEGH est le gap du composé AlSb.
b étant le paramètre de courbure (bowing) qui est souvent prouvé par l'expérience. L'origine
du bowing est dû à l'aspect structural et au désordre compositionnel qui sont très dominant
dans les fluctuations de l'alliage ternaire. Sachant que la VCA néglige l’effet de désordre
I.8. L’approximation du cristal virtuel améliorée
La méthode du pseudo potentiel empirique EPM couplée à l'approximation du cristal

virtuel VCA pour le calcul des structures de bande des alliages, donne de bon résultat, sauf
que la VCA rend les paramètres du bowing du gap principal aussi faibles que ceux des
résultats trouvés expérimentalement, pour palier à ce problème, on doit prendre en compte
l’effet du désordre dans la VCA. Baldereschi et Maschke en 1975 [19], ont traité l'effet du
désordre compositionnel comme une perturbation, mais leur approche de perturbation ne
pouvait pas expliquer le paramètre de courbure. Dans l’approximation du cristal virtuel
améliorée VCAA, le potentiel dû au désordre compositionnel est ajouté au potentiel
périodique du cristal, le pseudo potentiel du système s’écrit sous la forme [15] :
V(r) = VQRD (r) + VSéU (r) I − 12
VQRD (r): Potentiel périodique du cristal virtuel.
VSéU (r): Potentiel non périodique dû au désordre compositionnel.
Notons P comme étant le paramètre d'ajustement produisant le bowing Enfin en

additionnant ce potentiel du désordre avec le potentiel périodique effectif du cristal, on aura
l'expression finale du pseudopotentiel comme suit :
45
V(r) = VQRD − PVx(1 − x)W X ∆Zr − R \ ] I − 13
\
12
Le potentiel de l'alliage, calculé par la VCA auquel on ajoute le potentiel non périodique
du désordre, est donné par [26]:
`
V(r )= xVDEF (r) + (1 − x)VÊG_ − PVx(1 − x)Wa [VDEF (r) − VDEGH (r)W I − 14
Ce même potentiel peut s’écrire sous forme de séries de Fourrier:

V(r) = ∑d V(G)ecde I − 15
Avec V(G) est le facteur de forme de l'alliage ternaire.
13
[1] E. Mooser, « Introduction à la Physique du solide », 1 Ed. Presses polytechniques et

universitaires romandes (1993).
[2] H. Mathieu, « Physique des semi.conducteurs et des composants électronique », 4

Ed. Masson (1998).
[3] Beiser A. Physique Appliquée : Cours et Problèmes. McGrawHill, série Schaum.
[4] C. Ng –Hélène Ng , physique des semi.conducteurs, 3 Ed. Dunod, (2007).
[5] C. Kittel, « Physique de l’état solide », Ed. Dunod Université (1983).
[6] C. Kittel (trad. Nathalie Bardou, Évelyne Kolb), « Physique de l’état solide »,
Dunod, (1998).
[7] R. Omnès, " Comprendre la mécanique quantique ", édition EDP sciences, France
(1999)
[8] J. P. walter and M. L. Cohen, PhysRev. B4, (1971) 1877.
[9] A. Baldereschi, Phys. Rev. B7, (1973) 5212.
[10] G. Grynberg, A. Aspect, and C. Fabre, Introduction aux lasers et à l’optique

quantique, Ellipses, Paris, (1997).
[11] J.L. Coutaz « Optoélectronique » édition EDP Sciences, France, (2008).
[12] H. Mathieu, H. Fanet «Physique des semi.conducteurs et des composants

électroniques» édition Dunod, Paris, (2009).
[13] M.L. Cohen, J.R. Chelikowsky, «Electronic Structure and Optical properties of
Semiconductors », Springer Series in Solide.State Sciences (1989).
[14] S. Adachi, J, Appl. Phys. 61 (1987) 4869.
15
[15] N. M. Ravindra, P. Ganapathy, J. Choi, Energy gap.refractive index relations in

semiconductors. An overview, Infrared Phys. Technol. 50 (2007) 21.29.
[16] N. M. Ravindra, V. K. Srivastava, Variation of refractive index withenergy gap in

semiconductors, Infrared Phys. 19 (1979) 603.
[17] [V. P. Gupta, N. M. Ravindra, Comments on the Moss formula, Phys. Stat. Sol. B
100 (1980) 715.719.
[18] P. Hervé, L. K. J. Vandamme, General relation betweenrefractive index and

energy gap in semiconductors, Infrared Phys. Technol. 35 (1994) 609.615.
[19] R. R. Reddy, S. Anjaneyulu, Analysis of the Moss and Ravindra relations, Phys.
Stat. Sol. B 174 (1992) k91.k93.
[20] R. R. Reddy, Y. N. Ahammed, A study on the Moss relation, Infrared Phys.

Technol. 36 (1995) 825.830.
[21] A. N. Voloshinski, N. V. Ryzhanova, L. Yu. Vishnekov and A.G. Obukhov, Phys.

Stat. Sol. (b) 133, 735 (1986).
[22] S. J. Lee, T. S. Kown, K.Nahm and C. Koo Kim, J. Phys.Condens.Matter. 2, 3253

(1990).
[23] N. Bouarissa, N. Amrane and H. Aourag, Infrared Phys. Technol. 36, 755 (1995).
[24] N. Bouarissa and H. Aourag, Infrared Phys. Technol. 36, 973 (1995).
[25] J. E. Bernard and A. Zunger, Phys. Rev. B36, 3199, (1987).
[26] A. Baldereschi and K. Maschke, Solid Stat. Commun. (1975) 1699.
16
Chapitre II
METHODES DE CALCUL DES STRUCTURES
DE BANDES ELECTRONIQUES
Chapitre II Méthodes de calcul des structures de bandes électroniques
II.1. Introduction
Parallèlement au progrès effectué dans le domaine expérimental, les méthodes de
calcul de la structure de bandes électroniques, basées sur le développement des concepts de
base de nouveaux algorithmes, ont fait un impact crucial, afin de comprendre les propriétés et
les caractéristiques des matériaux et fournir par la suite, des données complémentaires pour
les expérimentateurs.
II.2. L’équation de Schrödinger d’un solide cristallin
Après les postulats de l’existence d’une onde associée à chaque particule et de

l’expression de sa longueur d’onde. L’équation de Schrödinger constitue le troisième postulat
fondamental de la mécanique quantique [1].
En 1926, guidé par les idées de L. De Broglie, Erwin Schrödinger a proposé une
équation dynamique à l’état stationnaire d’un système quantique appelé équation de
Schrödinger. Pour le cas d’un cristal, cette équation s’écrit sous la forme suivante:
( , … … )( , … , … )= ( , … , … ) II-1
Où :
: L’opérateur hamiltonien du cristal.
: Fonction d’onde du cristal.
E : L’énergie totale du cristal.
(i … Ne) représente les coordonnées des électrons, est le nombre d'électrons dans le
système.
(i … Nα) représente les coordonnées des noyaux, est le nombre d'atomes dans le
système.
Un cristal est une association de particules élémentaires, légères tels que les électrons et
de particules élémentaires lourdes, tels que les noyaux, qui sont en interactions permanentes
entre eux, sous l’action des forces électrostatiques attractives et répulsives.
17
On définit l’Hamiltonien comme étant la somme de toutes les énergies cinétiques et

potentielles présentent dans un cristal [3,4]
= + + + + II-2
Avec :
L’énergie cinétique des électrons :
−ℏ
= = ∆ II − 3
2
Où :
: Nombre total des électrons.
m : Masse de l’électron.
L’énergie cinétique des noyaux.
−ℏ
% %
= = ∆ II − 4
$
2&
$ $
Où :
Nα : Nombre total des noyaux.
M : Masse du noyau.
L’énergie d’interaction entre un électron i et un noyau j :
1 1
= = ) II −5
2 − ) 2
*) *)
18
L’énergie d’interaction entre un noyau α et un noyau β
1 ,$ ,- 1
= = $- II −6
2 | $ − - | 2
$*- $*-
Zαβ : charges des noyaux αet β.
L’énergie d’interaction entre un noyau α et un électron i
$ $
,$
=− =− II − 7
| − $|
$ $
En introduisant les équations II-3, II-4, II-5, II-6 et II-7 dans l’équation II-1,
l’Hamiltonien aura la forme suivante :
1 + + + + 2 ( , … , …) = ( , … … ) II-8
L’équation II-8, contient un grand nombre de particules et de ce fait il est difficile de

résoudre une telle équation à N éléments, des approximations s’imposent pour simplifier le
problème.
II.2.1. L’approximation adiabatique
L’une des premières approximations utilisée, lorsqu’on veut résoudre l’équation de

Schrödinger pour les systèmes complexes contenant un où deux électrons, c’est
l’approximation dite adiabatique de Born Oppenheimer, elle sépare le mouvement de
l’électron de celui du noyau, en se basant sur la notion de différence de masse entre celle du
noyau et celle de l’électron. Par conséquent, on peut négliger en 1ère approximation, l’énergie
cinétique des noyaux devant celle des électrons. Les états propres du système sont alors
caractérisés par des fonctions d’onde qui sont produit d’une fonction d’onde électronique par
une fonction d’onde nucléaire [1].
19
( , )= ( ) ( , ) II-9
L’équation de Schr34 dinger s’écrite alors :
( , ) =E ( ) ( , ) II-10
L’Hamiltonien est découplé en deux Hamiltoniens, l’un électronique et l’autre nucléaire

lié aux noyaux, ces derniers étant supposés fixes à leurs positions d’équilibre et l’équation de
Schrödinger devient :
( + + ) ( , )= ( ) ( , ) II-11
L’étude du mouvement des électrons, en fixant la position des noyaux est une bonne
approximation, elle réduit ainsi un grand nombre de variables, mais la résolution de l’équation
de Schrödinger reste toujours complexe à cause des interactions électrons.
Une deuxième approximation s’impose, d’elle-même, c’est l’approximation du champ

autocohérent.
II.2.2. Approximation du champ autocohérent
II.2.2.1. Approximation de Hartree
Le modèle de Hartree se base sur une expression particulière des solutions, en effet, on
va supposer que les électrons sont indépendants, ceci se traduit par [2] :
( , … )= ( ) ( )… ( ) II-12
5 : Position des électrons i avec i=1…n.
20
Tenant compte de cette hypothèse, et pour contourner la complexité d’un système poly-
électronique à un système mono-électronique, on considère que chaque électrons se déplace
dans le champ moyen des autres électrons et noyaux. L’équation de Schrödinger à une seule
particule (électron) est donnée comme suit :
( )= ( ) II-13
Avec :
ℏ6
=− ∆ + 89 ( )+ ( )
7
II-14
Où :
89 ( ) : Potentiel d’interaction électron-noyau.
( ) : Potentiel de Hartree pour le i ème électron, qui a la forme :
: ; <= (> ? ) B
( )= A
|> : > @ |
II-15
La densité des électrons C5 ( ) étant donnée par :
C( )=− D )( )D II − 16
)
*)
La somme est faite sur les Ne états mono-électroniques.
On remplace II-16 dans II-15 :
| ) ( ′)| B
( )= A II − 17
| − ′ |
)
*)
21
L’équation de Schrödinger aura la forme :
H L
G− ℏ ∆ − D ) ( @ )D B K
G 2 I A + 89 ( )K ( )= ( ) II − 18
G | − @| K
)
F *) J
L’équation II-18 est appelée équation à une particule de Hartree.
L’inconvénient de cette méthode, c’est qu’elle ne tient pas compte de l’antisymétrie de

la fonction d’onde.
II.2.2.2. Approximation de Hartree-Fock
Pour pallier le problème de l’antisymétrie imposant à la fonction d’onde de changer

de signe lors de la permutation de deux électrons, l’approximation Hartree Fock développe la
fonction d’onde sous forme d’un déterminant de Slater qui vérifie le principe de l’exclusion
de Pauli [3]
1 ( N) ( N )
( N … N )= Q R II − 19
O ! ( N) ( N )
Où :
σi : représente le spin
En injectant l’équation II-19dans l’équation de Schrödinger on a :
6
ℏ6 DVW 1> ? 2D DVW ∗ (>)V= (>)D
T− ∆ − ∑) ; A B − ∑) XY Y) ; A B + ( )[ ( )= ( )
7 |> : > ? | |> : > ? |
*) )*
II-20
22
Où :
XY Y) : Symbole de Kronecker
L’équation II-20 représente l’équation de Hartree-Fock à une particule, où on voit

apparaitre dans L’Hamiltonien le terme de Hartree et le terme résultant de l’application du
principe de l’exclusion de Pauli (terme d’échange), ce qui rend les calculs plus complexes.
II.2.2.3. Approximation de Hartree-Fock-Slater
Tenant compte du caractère non local du potentiel d’échange présent dans

l’approximation de Hartree-Fock, en 1951, Slater [4] propose un potentiel d’échange local de
la forme :
B<(>)
( ) = −6 ( ) ^B
\]
II-21
Avec :
α : Paramètre sans dimension.
Malgré la simplicité du potentiel par rapport à celui de Hartree-Fock, la méthode ne

conduit pas à des résultats satisfaisants.
II.3. Les méthodes de calcul de la structure de bandes électroniques
Il existe plusieurs méthodes théoriques pour calculer la structure de bandes

électroniques, elles se basent toutes sur l’approximation de Hartree Fock..
II.3.1. Méthodes de calcul des électrons libres
C’est un modèle développé par Arnold Sommerfeld [8], il considère que le potentiel
dans l’équation de Schrödinger est nul, du fait que les électrons n’interagissent pas entre eux,
les électrons sont libres, donc, tous les niveaux d’énergies leurs sont permis, ces niveaux
forment une bande continue. L’énergie totale du système ne comprend que l’énergie cinétique
des électrons, leur énergie potentielle est négligée.
23
Cette approche permet d’expliquer les propriétés physiques des métaux, mais elle
présente des limites pour la distinction entre isolants et semi-conducteurs.
II.3.2. Méthodes de calcul des électrons presque libres
Dans cette méthode, les électrons sont presque libres, ce qui revient à dire, qu’ils sont
faiblement liés, donc, il suffit d’une faible énergie pour qu’ils se déplacent à la bande de
conduction. Ceci se traduit aussi par, l’existence d’un potentiel périodique du cristal qui est
faible, et qui va être considéré comme perturbation périodique du réseau.
II.3.3. Méthodes de calcul des ondes planes
Pour être conforme avec la périodicité du réseau cristallin, on choisit une forme
particulière de la fonction d’onde établie par Bloch :
_( ) = `_ ( ) II -22
_>
`_ ( ) = ab b>
II − 23
b
Avec :
G : Vecteur du réseau réciproque définit par cd . fd=

]
7
m : Entier naturel
De même le potentiel périodique du cristal V(r) peut être décomposé en série de

Fourier [7] :
( )= (c) b>
II − 24
b
24
En remplaçant les équations II-22 et II-24 dans l’équation de Schrödinger :
ℏ
g− ∆ + ( )h ( ) = ( ) II − 25
2
En multipliant par _?>

et en intégrant sur A B on obtient :
ℏ
− a(i @ )i @ + a(i @ − c) (c) = (i)a(i) II − 26
2
b
En supposant que les V(G) sont connus, on peut calculer les E(k) puis déterminer les C(i′).
II.3.4. Méthode de calcul des ondes planes augmentées
Pour pallier les problèmes des conditions aux limites de la méthode cellulaire aux
frontières de la maille élémentaire, la méthode de calcul des ondes planes augmentées a été
proposée. Sachant que le solide est constitué par un ensemble d’atomes et chaque atome est
constitué par un cœur sphérique de rayon j , où règne un potentiel sphérique et une région
interstitielle où le potentiel est constant. Ce type de potentiel est appelé le potentiel de muffin-
tin (nid d’abeille) [8,9]
II.3.5. Méthodes de calcul des ondes planes orthogonales
Cette méthode a été proposée par Herring en 1940 [10]. En se basant sur le fait
qu’approximer les fonctions d’ondes par des fonctions d’ondes planes, n’est pas une bonne
solution, et du fait que dans la région du cœur la fonction d’onde oscille rapidement, il
propose des ondes planes orthogonales au niveau du cœur et qui ont la forme :
_( )= + k l_m ( ) II − 27
_>
25
Où :
C : Définit le cœur.
k : Constante calculée par l’orthogonalité de _( ) avec chaque niveau du cœur.
l_m ( ): Fonction d’onde atomique.
II.3.6. Méthodes des pseudopotentiels
La méthode du pseudo potentiel, comme la méthode OPW, utilise les propriétés

d'orthogonalité des états de valence et conduction avec les états du cœur. Mais dans le
formalisme du pseudopotentiel l'effet de l'orthogonalité est inclus dans le potentiel sous la
forme d'un potentiel équivalent appelé pseudopotentiel. L'effet d'orthogonalité aux états de
cœur revient à extraire du potentiel cristallin la contribution rapidement variable de la région
du cœur. Le pseudopotentiel est alors lentement variable et se prête bien à une approche du
problème en termes de perturbation [11].
II.4. Conclusion
Les calculs de la structure de bandes électroniques reposent sur un grand nombre

d’hypothèses présentant chacune des difficultés propres dont il faut tester la validité par un
retour à l’expérience. Ces méthodes sont adoptées selon le type de propriétés que l’on veut
étudier , certaines nécessitent un calcul numérique très lourd, comme les méthodes des
premiers principes, d’autres comme l’EPM permet de calculer d’une façon relativement
économique en terme du temps, des paramètres physiques en se basant sur l’utilisation d’un
nombre réduit de facteurs de formes pour le calcul de la structure de bandes et c’est la raison
principale qui rend cette méthode pratique.
26
27
[1] P. Kiréev, « La Physique des semiconducteurs », Ed. Mir, Moscou (1975).
[2] C. Kittel, « Physique de l’état solide », Dunod université (2001).
[3] J.A Pople et R.K. Nesbet, J. Chem. Phys. 22, (1954) 571.
[4] J.C.Slater, Phys. Rev 81, (1951) 385.
[6] S. Bounab, Thèse de Magister, Univ de Constantine, Algéria, (2002).
[7] P. Kiréev, ‘ La physique des semiconducteurs’, Ed Mir. Moscou, (1975).
[8] J. C. Slater Phys. Rev. 51, (1937) 846.
[9] N. W .Aschroft and N. D. Mermin, Solid State Physics, Ed Holt, Renehart and
Winston Philadelphia, (1976).
[10] C. Herring, Phys. Rev. 57, (1940) 1169.
[11] H. Mathieu, H. Fanet « Physique des semi-conducteurs et des composants

électroniques» 6 édition Dunod, Paris, (2009).
Chapitre III
METHODES DES PSEUDOPOTENTIELS

III.1.Introduction :
La technique du pseudopotentiel est une véritable révolution dans la théorie des
structures de bande. Cette méthode s’est avérée, depuis son introduction, très efficace dans
l’investigation et la compréhension des propriétés électroniques des solides.
Le pseudopotentiel a été introduit pour la première fois par Fermi (1934) pour l’étude des
niveaux atomiques observés.[1] en 1940 Convers Herring propose la méthode des ondes
planes orthogonales [2]. La méthode permit de mieux comprendre la nature de structure de
bande des matériaux semi-conducteurs tels que le Silicium et le Germanium et fut la première
à expliquer de manière théorique que le Silicium est un matériau à gap indirect [3]. C’est par
une reformulation de la méthode des ondes planes orthogonales d’Herring, qu’en 1959
Phillips et Kleinman développèrent la première approche formelle (non empirique) de la
notion de pseudopotentiel[4]
III.2. Principe de la méthode des pseudopotentiels :
Dès le XIX ème siècle, les physiciens ont su démontrer qu’une grande partie des
propriétés physique et chimiques des matériaux ne dépendent que du comportement des
électrons de valence. En effet, les orbitales de cœur sont les plus basses en énergie, sont
localisées près du noyau, sont très peu sensible à l’environnement et ne participent pas aux
liaisons chimiques. En outre, elles sont difficiles à représenter sur une base d’ondes planes car
elles possèdent généralement de fortes oscillations autour des noyaux. En revanche, les
orbitales de valence sont peu localisées et s’étendent donc loin du noyau. Ce sont elles qui
déterminent au premier ordre les propriétés physico-chimiques [5]. C’est sur cette propriété
que se base l’approximation dite cœur gelé.
Cette méthode remplace l'effet des électrons du cœur par un pseudopotentiel effectif. Le
système que l'on traite à présent n'est plus le système {noyau nu + électrons} mais {[noyau nu
+ électrons de cœur] + électrons de valence} qui est équivalente à {"ions" + électrons de
valence}[5].
La méthode des ondes planes orthogonales est considérée comme l’ancêtre directe de la
notion du pseudopotentiel. La méthode du pseudopotentiel de PhillipsKleinman est la
première à montrer que la condition d’orthogonalité dans la région decœur entre les états de
cœur et de valence agit comme potentiel répulsif qui tend à s’opposer au potentiel nucléaire
attractif ressenti par les électrons de valence, ainsi, par effet d’annulation on obtient un
28
potentiel faiblement répulsif ou pseudopotentiel[4,6]. Ce dernier peut être traité en utilisant la

méthode des électrons presque libres pour résoudre l’équation de Schrödinger. L’avantage de
cette approximation est que le nombre d’électrons apparaissant de manière explicite dans les
calculs s’en trouve fortement réduit, seuls les électrons de valence sont pris en compte.
III.3. Théorie de Phillips et Kleinman :
Cette théorie consiste donc à remplacer le terme de répulsion cœur-valence (qui est une
interaction non locale) par un potentiel effectif local capable de prendre en compte la
condition d’orthogonalité entre le cœur et la valence afin d’empêcher l’effondrement des
orbitales de valence dans la région spatiale de cœur, mais sans explicitement l’imposer [7,8].
La fonction d’onde orthogonale est la somme d’une onde plane et les états atomiques occupés
du cœur.
= + III − 1
La contrainte d’orthogonalité entre le cœur et la fonction d’onde est défini quant à elle par :
= 0 III − 2
La constante de l’équation III-1 est calculée par la condition d’orthogonalité [9].
=− =− ∗
III − 3
On remplace dans l’équation III-1
= − III − 4
29
La fonction est solution de l’équation de Schrödinger.
!"# =$ III − 5
Où :
!"# : Hamiltonien de Phillips-Kleinman dont l’expression sera explicitée plus loin. En

remplaçant l’équation III-4 dans l’équation III-5 :
!"# − !"# =$ & − ' III − 6
Si on applique l’Hamiltonien ̂)* sur une fonction de valence et une fonction de cœur, on
peut écrire les équations aux valeurs propres suivantes :
!"# + = $, + III − 7
!"# = $, III − 8
Où $ et $ sont les énergies de valence et decoeur respectivement
Ainsi, il en découle que :
!"# −$ =$ −$ III − 9
!"# − $ + $ =$ III − 10
!"# + $ −$ =$ III − 11
30
On pose :
$ −$ = 01 III − 12
Où :
02 : Potentiel répulsif positif
!"# + 01 =$ III − 13
Ainsi l’équation III-12 aura la forme :
3 !"# + 01 4 =$ III − 14
Avec :
)5
!"# = − + 0 III − 15
26
Où :
0 : Potentiel attractif du cœur.
D’où l’équation III-14 s’écrit comme suit :
31
)5
7− +0 + 01 8 =$ III − 16
26
On pose 0"9 = 0 + 01 qui représente le pseudopotentiel, l’équation III-16 sera

reformulée de la manière suivante :
)5
:− + 0"9 ; =$ III − 17
26
En posant,0"9 = 0 + 01 comme étant la somme du potentiel ionique du cœur 0< et du

potentiel répulsif positif02, cette sommation des deux quantités l’une positive et l’autre
négative nous donne un faible potentiel d’où l’appellation du pseudopotentiel, cette approche
réduit au maximum les calculs en outils informatique.
La résolution de l’équation de Schrödinger III-17 sera traitée par la méthode des électrons
presque libres, du fait que le pseudopotentiel0"9 sera considéré comme une perturbation du
système.
III.4. Modèle de potentiel :
D’après l’équation III-17 l’Hamiltonien de Phillips-Kleinman d’un électron dans un cristal est
une sommation de son énergie cinétique et de son pseudopotentiel.
)5
!"# = − + 0"9 III − 18
26
Plusieurs forme de potentiel ont été proposé on a le modèle de Heine Abarenkov, le modèle
d’Ashcroft, le modèle de Shaw et le modèle de Hallers. On va uniquement présenter deux
modèles, celui d’ Heine-Abarenkov et d’Ashcroft.
III.4.1. Modèle de potentiel de Heine-Abarenkov :

32
Heine et Abarenkov ont choisi un rayon de 2= > 2 (rayon du cœur ionique) au delà duquel
le potentiel est purement Coulombien. A l'intérieur de la sphère de rayon 2?, il est supposé
constant. L’expression du modèle de potentiel de Heine et Abarenkov pour l'ion libre
comprend deux termes [10], un terme local de la forme :
0 < 2=
0"9
@
= A−CD
> 2= III − 19
Où :
Z : Nombre d’électrons de valence
Et un terme non local ayant pour forme
− FG , $ IG )G < 2=
0"9
@
E
KL
0 < 2= III − 20
Où :
FM ( , $) : Energie des états profonds.
)M : Opérateur de projection qui donne la composante au nombre quantique l, dans l’espace

des harmoniques sphériques.
( ) : Fonction qui représente l’effet de l’état du cœur possédant plusieurs formes telle que la
forme gaussienne, carré et la forme d’Ashcroft. Les paramètres FM ( , $) doivent être ajustés
de façon à permettre au modèle de potentiel d'avoir le même spectre d'énergie que le vrai
potentiel. Ceci est réalisé en appliquant la méthode du défaut quantique [11].
III.4.2. Modèle de potentiel d’Ashcroft :
33
Ashcroft [12] a repris le modèle de potentiel de Heine-Abarenkov avec une forme plus simple
nulle à l’intérieur, son modèle possède un seul paramètre le rayon du cœur 2

Son modèle, pour lequel l'interaction est purement Coulombienne à I’extérieur du cœur et
0 < 2
0"9 A−CD
> 2 III − 21
toutes constantes FM à l’intérieur du cœur ionique sont nulles, cependant c’est un modèle local
Le modèle de potentiel de Heine-Abarenkov est ainsi simplifié par Ashcroft, il considère que
et ne dépend pas de l’énergie.
III. 5. Méthode empirique du pseudopotentiel (EPM) :
Les pseudopotentiels sont parfois qualifiés d’empiriques, ce qui veut dire qu’ils ne sont pas
obtenus par calcul mais paramétrés pour reproduire aux mieux des résultats expérimentaux de
références [13]. L’approche empirique du pseudo potentiel, est une méthode qui se base, pour
potentiel, appelé les facteurs de forme 0N (O) et0P (O). Les valeurs de (O) des éléments sont
le calcul de la structure de bande, sur l’ajustement théorique des paramètres du pseudo
classifiées et représentées graphiquement dans l’article de Cohen et Heine [14].Ces derniers

seront paramétrés pour reproduire au mieux les résultats expérimentaux.
Le concept de pseudo potentiel empirique fût utilisé pour la première fois pour le calcul des
niveaux énergétiques des métaux alcalins [15,16]. Mais l’approche a englobé avec succès une
douzaine de semi-conducteurs de structure Diamant et Zinc blende, tels que le silicium et le
germanium [3]. Il existe deux approximations Empiriques, locales et non locales.
III.5.1. Approximation Empirique locale :
Cette approximation a été largement utilisée depuis son introduction par Chelikowsky et
Cohen [17,18]. L’Hamiltonien total du cristal peut être écrit comme une sommation des
Hamiltoniens local et non local, en plus de l’Hamiltonien spin-orbit.
34
!"# = !"#
@
+ !"#
Q @
+ !"#
9"
III − 22
Où l’Hamiltonien local !"#

@
vérifie l’équation de Schrödinger à un électron de valence.
L’équation de Schrödinger est donnée par :
!"#
@
=$ III − 23
ℏ5
!"#
@
=− ∆ + 0"9 III − 24
26
L’équation II-23 aura la forme:
ℏ5
T− ∆ + 0"9 U =$ III − 25
26
Avec :
0"9 :Pseudopotentiel périodique, donc il satisfait la relation 0"9 = 0"9 +2
R : Vecteur du réseau direct.
Selon le théorème de Bloch [19] qui donne les solutions de l’équation de Schrödinger pour un
périodique VW ( ) par une onde plane.

potentiel périodique. La fonction d’onde de Bloch est de la forme d’un produit d’un facteur
=V X
III − 26
Où :
35
VW ( ) : Fonction périodique qui a la périodicité du réseau cristallin et qui a la forme
V = V ,X
X
III − 27
X
Où :
k : Vecteur d’onde.
G : Vecteur du réseau réciproque.
En remplaçant l’équation III-27 dans l’équation III-26, la fonction d’onde s’écrit comme suit :
= V ,X
YX
= V YX
YX
III − 28
X X
On substitue l’équation III-28 dans l’équation III-25 :
ℏ5
− ∆& V YX
' + 0"9 V YX
= 0"9 V YX
III − 29
26 YX YX YX
X X X
En développant le 1er terme de l’équation III-29, on aura :
ℏ5 ℏ5
− V YX Z∇
5 YX
\= V YX Z|^ + O|5 YX
\ III − 30
26 26
X X
En remplace l’équation III-30 dans l’équation III-29 :
ℏ5
V YX Z|^ + O|5 YX
\ + 0"9 V YX
= 0"9 V YX
26 YX YX
X X X
36
III-31
_ YX`
En multipliant les deux côtés de l’équation III-31 par
ℏ5
V YX Z|^ + O|5 3X_X a 4
\ + 0"9 V 3X_X a 4
= 0"9 V 3X_X a 4
26 YX YX
X X X
III − 32
En intégrant sur tout le volume du cristal, le premier terme à l’extrême gauche de l’équation
III-32 s’écrit:
YJ
ℏ5 ℏ5 YJ
V YX Z|^ + O|
5 3X_X a 4
\ = V 3X_X a 4
26 26 YX
_J
_J X X
III-33
On passe à la transformée de Fourier :
ℏ5 YJ
ℏ5
|^ + O|5 V 3X_X a 4
= 2b |^ + O|5 V c O − O`
26 YX
_J 26 YX
X X
III − 34
ℏ5 ℏ5
2b |^ + O|5 V YX c O − O = 2b
` |^ + O|5 V#YXa III − 35
26 26
X
Le deuxième terme à gauche de l’équation III-32 après intégration sur tout le volume du
cristal :
37
YJ YJ
0"9 V#YX 3X_X a 4
= V#YX 0"9 X_X`
III − 36
_J X X _J
On passe à la transformée de Fourier :
YJ
V#YX 0"9 3X_X a 4
= 2b V#YX 0"9 O − O ` III − 37
X _J X
On intègre maintenant le second terme de l’équation III-32 sur tout le volume du

cristal et par application de la transformée de Fourier :
YJ
$ V#YX X_X a
= 2b $ V#YX c O − O ` III − 38
X _J X
2b $ V#YX c O − O ` = 2b $ V#YXa III − 39

X
En substituant les équations III-35, III-37 et III-39 dans l’équation III-32
ℏ5
|^ + O|5 V#YX` + V#YX 0"9 O − O ` = c O − O ` $ V#YXa III − 40
26
X
On multiplie par∑X c O − O ` le premier terme à gauche de l’équation III-40 et le terme à

droite de la même équation :
38
ℏ5 c O − O ` $ V#YXa
c O−O |^ + O|5 V#YXa +
`
V#YX 0"9 O−O `
=
26
X X e III − 42
c O − O ` = f 0 gW O ≠ O′ III − 42
Où :
j gW O = O′
e
Avec :
N : Nombre total de cellule dans le volume V
ℏ5
Tk |^ + O|5 − $ l c O − O ` + 0"9 O − O ` U = 0 III − 43
26
X
L’équation III-43 est une équation linéaire homogène, qui admet une solution non triviale si le
déterminant de l’équation séculaire est nul :
ℏ5
m n ok |^ + O|5 − $ l c O − O ` + 0"9 O − O ` o = 0 III − 44
26
La périodicité du cristal permet d’écrire le terme 0"9 O − O ` dans l’équation III-44 sous la
forme de la transformée de Fourier :
1
0 O − O` = 0"9 _ 3X_X a 4
III − 45
0
Dans l’approximation locale, le potentiel 0"9 est représenté par une superposition linéaire
de potentiels atomiques :
r G
0"9 = 0p − 2q − s III − 46
q
39
Avec : m : Nombre de maille élémentaire dans le solide.
l : Nombre d’atomes par maille élémentaire.
2q : Vecteur de translation dans le réseau réel.
s : Vecteur de position liant deux atomes dans la maille élémentaire.
En substituant l’équation III-45 dans l’équation III-44 :
1
0 O − O` = 0p − 2q − s _ 3X_X a 4
III − 47
0
r,G
One pose `
= − 2q − s
1
0 O − O` = 0p ` _ 3X_X a 4 a Y1t Yu `
III − 48
0
r,G
1
0 O − O` = _ 3X_X a 41t _ 3X_X a 4uv
0p _ 3X_X a 4
III − 49
0
r,G
1
0 O − O` = _ 3X_X a 41t _ 3X_X a 4uv
0p _ 3X_X a 4
III − 50
0
r G
Or l’intégrale sur 0 est équivalente à une intégrale sur Ω qui est Volume de la maille
élémentaire de Wigner-Seitz (W.S) et qui contient toute l’information physique sur le système
0p ` _ 3X_X a 4 ` `
= 0p _ 3X_X a 4
III − 51
Sachant que :
_ 3X_X a 41t
= 1 = j DGG III − 52
r r
Ou :
j DGG : Nombre de cellule dans le volume V.
j DGG
0 O − O` = _ X_X a uv
0p ` _ 3X_X a 4 ` `
III − 53
0
G
40
j DGG Ω 1
0 O − O` = _ X_X a uv
0p ` _ 3X_X a 4 ` `
III − 54
0 Ω
G
1 1
0 O − O` = _ X_X a uv
0p ` _ 3X_X a 4 ` `
III − 55
j Ω
G
Finalement on a :
0 O − O ` ≈ Š O − O ` 0‹ O − O ` III-56
Où 0‹ O − O ` :Facteur de forme atomique, qui représente la transformée de Fourier du

potentiel atomique à l’intérieur de la cellule W.S.
Š O − O ` : Facteur de structure qui ne dépend uniquement que de l’emplacement des ions

dans la cellule W.S. On a donc:
1
0‹ O − O ` = 0p ` _ 3X_X a 4 ` `
III − 57
Ω
1
Š O − O` = _ 3X_X a 4uv
III − 58
j
G
Selon Cohen et Bergestresser [19, 20], les facteurs de forme sont déterminés à partir
des gaps optiques expérimentaux, d’où le caractère empirique de l’approche.
Appliquons les équations III-57 et III-58 à la structure Diamant ou Zinc-blende qui

possède deux atomes par motif. Considérant pour cela les notations suivantes :
Ωp : Volume atomique
N : Nombre de maille élémentaire.
Ω = 2jΩp :Volume du cristal .
•L : Paramètre de réseau.
On a :
0 O − O ` = 0XXa ,• _ 3X_X a 4u‘

+ 0XXa ,5 _ 3X_X a 4u’
III − 59
41
s• , s5 :Position des deux atomes présents dans la maille élémentaire pour la structure Diamant
ou Zinc-Blende , les deux atomes sont identiques, d’où s• = −s5 = s
•L •L
s• = s = 1,1,1 n s5 = −s = − 1,1,1
8 8
L’équation III-59 s’écrite sous la forme :
0 O − O ` = 0XXa ,• _ 3X_X a 4u
+ 0XXa ,5 _ 3X_X a 4u
III − 60
0 O − O ` = 0 9 O − O ` <“g[` O − O′ s] + W0 P O − O ` ŠW–[ O − O′ s] III − 61
0 O − O ` = 0 9 O − O ` Š 9 O − O ` + W0 P O − O ` Š P O − O ` III − 62
Avec les facteurs de structure :
Š 9 O − O ` = <“g[` O − O′ s] III-63
Š P O − O ` = ŠW–[ O − O ` s] III − 64
Les facteurs de formes symétrique et antisymétrique :
1
0|X_X`|
9
= 30 a + 0XXa ,5 4 III − 65
2 XX ,•
1
0|X_X`|
P
= 30 a − 0XXa ,5 4 III − 66
2 XX ,•
Le pseudopotentiel des deux ions par maille :
42
2
0XXa , = 0 _ 3X_X a 4
III − 67
Ω
Où :
i=1 et 2 : sont l’anion et le cation pour la structure Zinc-blende ou Diamant. On
substitue l’équation III-62 dans l’équation III-43 :
ℏ5
Tk |^ + O|5 − $# l cXX` + 0XX`
9
ŠXX`
9
+ W0XX`
P
ŠXX`
P
U V#YX` = 0 III − 68
26
X
L’équation III-68 est une équation linéaire homogène écrite en fonction
des facteurs de forme et facteurs de structure, admet des solutions si le déterminant est
nul.
ℏ5
m n Tk |^ + O|5 − $# l cXX` + 0XX`
9
ŠXX`
9
+ W0XX`
P
ŠXX`
P
U V#YX` = 0 III − 69
26
X
Dans l’approche du pseudopotentiel local, on ne peut pas connaitre les
informations sur les états de la région du cœur, cependant, on peut se baser sur une
interpolation des facteurs de forme symétrique et antisymétrique et les comparer avec
les données expérimentales afin d’obtenir une structure de bande des électrons de
valence exacte. Le potentiel atomique possède une symétrie sphérique ce qui rend les
facteurs de forme dépendant de la magnitude (portée) de G [21]. Sachant que le
potentiel s’annule pour des grandes valeurs de G, l’équation séculaire II-69 converge à
partir des valeurs de
5˜ •⁄5
|O| < |O′| = — ™ √11
‹
III-70
43
Pour satisfaire cette condition, on a besoin de trois facteurs de forme symétriques
0 9 O et de trois facteurs de forme antisymétrique
0 9 111 = 0 9
œ 0 9 220 = 0•9
0 9 311 = 0••9
III-71
0 P 111 = 0 P
œ 0 P 220 = 0•P III-72
0 P 311 = 0••P
Les facteurs de forme symétrique , s’obtiennent à partir des structures monoatomiques telles
que le Ge, Si, etc…..d’après la méthode de Heine-Abarenkov[22] ,Tandis que les facteurs de
forme antisymétrique se déterminent expérimentalement d’après Cohen et
Bergestresser[19,20]
La Figure.III.1 ci-dessous décrit les démarches de la procédure de l’approche du

pseudopotentiel empirique local [23] .
44
Lire V G
0′O = 0‹9 Cos Os + W0‹9 ŠW–Os
ℏ5
!"9 = − ∆+0 2
26
Résoudre l’équation de Schrödinger
dinger
<•M¨©M $ * n $
<“6ª• M g «•M © g •« ¨ M′ ¬ªé W –¨
FMné 0 O
figure.III.1 :L’approche du pseudopotentiel empirique local
L’approximation du pseudopotentiel empirique local se base sur le bon choix des facteurs de
forme V(G) de la structure Zinc-blende à étudier et dont on connait le paramètre de réseau.
On injecte ces paramètres V(G) dans l’équation de Schrödinger, les potentiels construits et
après résolution de l’équation séculaire par la méthode des électrons presque libres on obtient
les valeurs propres des énergies$ , les coefficients des fonctions d’onde et les fonctions
d’ondes
Ces valeurs sont comparées avec ceux de l’expérience si un grand écart se présente, les
facteurs de forme sont modifiés et le processus est répété jusqu'à la convergence entre les
valeurs expérimentales et les valeurs calculées.
45
III-6-conclusion :
Lorsque la molécule que l'on étudie contient des atomes lourds, le nombre d'électrons peut
devenir très élevé et la résolution de l'équation de Schrödinger s'entrouve alourdie. Par
ailleurs, il est bien connu depuis le XIXe siècle que les propriétés chimiques des atomes et des
molécules proviennent des électrons les plus externes ou électrons de valence, les électrons
les plus internes ou électrons de cœur étant peu impliqués. Dans le but d'alléger les calculs, il
est donc possible d'éviter le traitement explicite des électrons de cœur, en les remplaçant par
un potentiel effectif agissant sur les électrons de valence. Un tel potentiel est appelé
pseudopotentiel. Cette idée fut introduite pour la première fois en chimie quantique en 1959
par Phillips et Kleinmann et s'est surtout développée à partir des années 70. Dans ce chapitre,
le calcul de structures électroniques par la méthode des pseudopotentiels empiriques (EPM) a
été abordé en équations analytiques détaillées.
46
[1] E. Fermi, Nuovo Cimento 11, (1934) 157.
[2] W.C. Herring , Phys. Rev, Vol. 57, (1940) 1169-1177.
[3] F. Bassani, Phys. Rev. , vol. 108, (1957) 263-264.
[4] J.C. Phillips, L. Kleinman, Phys. Rev., (1959) 116.
[5] N. Richard « Actinides et terres rares sous pression : approche pseudopotentiel »
Rapport CEA, France, (2002).
[6] E. Antoncik , J. Phys. Chem. Solids 10, (1959) 314.
[7] Ludovic Castro, thèse de doctorat, univ Toulouse, (2012).
[8] J. C. Phillips, Phys. Rev. 112, 685 (1958).
[9] Jiseok Kim, thèse de doctorat, univ USA, (2011).
[10] N. Gargour, thèse de doctorat, univ Metz, (1990).
[11] F. S. Ham, Solid State Physics 1 (1955) 127. Editeur: Seitz et Turnbull
[12] N. W. Ashcroft, Phys. Lett. 23 (1966) 48.
[13] N. Bouarissa, Philo. Mag. B 80 (2000) 1743.
[14] M. L. Cohen and T. K. Bergstresser, Phys. Rev. 141, (1966) 789.
[15] H. Hellmaan, The Journal of Chemical Physics, vol. 3, 1935, p. 61.
[16] H. Hellmann and W. Kassatotschkin (1925) Acta Physicochim. U. R. S. S. 5, 23.
[17] Chelikowsky, J.R., and Cohen, M.L. Physical Review B 10, 12 (1974), 5095–5107.
[19] E. P. Wigner and F. Seitz, Phys. Rev. 43, 804 , 1933.
[20] M. L. Cohen and T. Bergestresser Phys. Rev. 164, 1069 (1976).
[22] I. V. Abarenkov and V. Heine, Phil. Mag. 13, 529 (1965).
[23] N. Bouarissa, Materials Chemistry and Physics 124 (2010) 336–341.
47
Chapitre IV
Résultats et discussions
IV.1. Etude des propriétés électroniques
La méthode du pseudopotentiel empirique locale combinée avec l’approximation

du cristal virtuel, qui introduit l’effet du désordre comme un potentiel effectif, est
utilisée pour étudier les propriétés électroniques des alliages semi-conducteurs : à savoir
la structure de bandes, le gap d’énergie, la densité de charge électronique et la masse
effective des électrons, ainsi que leurs variations en fonction de la composition x de
l’alliage semiconducteur.
Le pseudopotentiel est donné en fonction des paramètres du pseudopotentiel V( ),

appelés aussi facteurs de forme, et qui ne sont autre que la Transformée de Fourier du
potentiel atomique à l’intérieur de la cellule Wigner-Seitz.
= IV − 1
Où :
V(G) : Facteur de forme
La nature empirique de la méthode du pseudopotentiel consiste à ajuster les

facteurs de forme, afin de réaliser l'accord le plus étroit des niveaux d’énergies calculés
avec les valeurs théoriques. Les facteurs de forme des composés binaires
semiconducteurs AlSb et AlP sont ajustés et illustrés dans le Tableau. IV.1, ainsi que
leurs paramètres de réseau.[1-2].
Tableau. IV.1 Les facteurs de formes symétriques (Vs) et antisymétriques (VA) ajustés
en (Ryd) des composés binaires semi-conducteurs AlSb et AlP et leurs paramètres de
réseaux.
Composant Les facteurs de forme a (Å)
Binaire VS(3) VS(8) VS(11) VA(3) VA(4) VA(11)
AlSb -0.225597 0.028086 0.062230 0.007109 0.058960 0.004544 6.1355
AlP -0.241305 0.021262 0.249336 0.030061 0.116000 0.232543 5.4635
Connaissant les valeurs des paramètres pseudopotentiels (Facteurs de formes) on

peut calculer les gaps d’énergies qui sont illustrés au Tableau. IV.2. On remarque que
ces valeurs sont en très bon accord avec l’expérience rapportée par la référence [3,4]:
48
Tableau.IV.2. Gaps d’énergies calculés et expérimentaux des composés binaires semi-

conducteurs AlSb et AlP.
Composant binaire EГ-Г (eV) EГ-X (eV) EГ-L (eV)
AlSb 2.30 a,c 1.615 a,c 2.211 a,c
AlP 3.56 b,c 2.52 b,c 3.57 b,c
"
3, 4 et ! .
Le paramètre de réseau de l’alliage semiconducteur AlPxSb1-x , suit la loi de

Vegard [5], dont les valeurs sont présentées au Tableau.V.3.
*+ , -.,
=/ *+ + 1−/ *+ IV − 2
Tableau.IV.3. Paramètres de réseau de l’alliage semiconducteur AlPxSb1-x.

x a (A°)
0 6.1355
0.1 6.0683
0.2 6.0011
0.3 5.9339
0.4 5.8667
0.5 5.7995
0.6 5.7323
0.7 5.6651
0.8 5.5979
0.9 5.5307
1 5.4635
Les facteurs de forme symétriques et antisymétriques des alliages ternaires AlPxSb1-x

sont donnés par la relation ci-dessous : [6]
57
++23 = 1−/ *+ +/ *+ −4 / 1−/ 6 *+ − *+ IV − 3
Avec 4 : Paramètre ajustable qui simule l’effet de désordre.
49
IV.1.1. Structure de bandes électronique :
La structure de bande électronique de l’alliage ternaire AlPxSb1-x est calculée aux

points de hautes symétries dans la zone de Brillouin, en utilisant la VCA avec et sans
tenir compte de l’effet du désordre compositionnel.
Vu la similitude de ces structures pour différentes valeurs de x, seule la structure
de bande électronique de l’alliage AlP0,5Sb0,5 est représentée dans la figure. IV.3. Les
Figures. (V.1-V.2) illustrent les structures de bandes d’énergies des composés binaires
8 9: et AlP respectivement.
Pour l’interprétation de ces structures de bandes, on note que la référence zéro

d’énergie est le maximum de la bande de valence. Dans la Figure. IV.1 on constate que
le maximum de la bande de valence est au point Γ et que le minimum de la bande de
conduction est au point X, d’où le composé binaire AlSb est un semiconducteur à
bande interdite indirecte. Ce résultat est en accord avec les valeurs trouvées à la
référence [3]. La Figure IV.2 illustre la structure de bande du composé binaire AlP qui
est aussi à bande interdite indirect, un résultat qui est aussi en accord avec les valeurs de
la référence [4].
15
AlSb
10
Energie électronique (eV)
-5
-10
-15
Γ X W L Γ Κ X
Vecteur d'onde
Figure IV.1 Structure de bande électronique du composé binaire AlSb
50
15 AlP
10
Energie électronique (ev)

5
-5
-10
-15
Γ X W L Γ Κ X
Vecteur d'onde k
Figure IV.2 Structure de bande électronique du composé binaire AlP
La structure de bande d’énergie de l’alliage ternaire AlP0,5Sb0,5 est

schématisée sur la Figure. IV.3, on voit clairement que l’alliage est un semiconducteur à
gap d’énergie indirect.
15
VCA
AlP0.5Sb0.5
VCAA
10
Energie électronique (eV)
-5
-10
Γ X W L Γ Κ X
Veteur d'onde (k)
Figure IV.3 Structure de bande électronique de l’alliage ternaire AlP0,5Sb0,5
Les trois courbes présentent quatre bandes de valence distinctes, la plus basse des
bandes de valence présente la forme habituelle attendue pour les semiconducteurs
51
Zinc-blende [7]. On note que ces quatre bandes sont moins dispersives que les bandes
de conductions, ceci est dû au fait que les électrons de conduction sont délocalisés.
Un gap d’énergie est remarqué entre la première et la deuxième bande de

valence au point X, ce gap est relié à la différence des potentiels des cations et anions.
Ce gap «antisymétrique » a été proposé comme une mesure qualitative de l'ionicité du
cristal [8].
La largeur totale de la bande de valence semble augmenter avec l'augmentation de
la teneur en P, ceci reflète le changement dans l'ionicité d'AlSb quand plus d'atomes de
P sont ajoutés dedans.
La Figure IV.3 montre un effet de désordre que important sur la structure de

bande électronique de l’alliage ternaire au point x=0.5, cet effet ne peut être négligé.
IV.1.2 Gap d'énergie
Pour connaitre les transitions optiques, l’étudie de la dépendance du gap

énergétique direct =>> et indirect =>? en fonction de la composition x avec et sans tenir
compte de l'effet du désordre s’avère nécessaire.
La variation de l’énergie du gap direct =>> en fonction de x est illustrée à la Figure

IV.4, la courbe a été tracé en utilisant les deux méthodes la VCA et la VCA améliorée.
4.0
AIPxSb1-x
3.5
3.0
Gap d'énergie (eV)
2.5
2.0
1.5
1.0
0.5
0.0 0.2 0.4 0.6 0.8 1.0
AlSb Composition x AlP
Figure IV.4 Variation du gap d’énergie direct E>> en fonction de la composition x avec
et sans désordre de l’alliage ternaire AlPA Sb5BA
52
Notons qu’en partant de AlSb (x = 0) à AlP (x = 1), la bande interdite directe (Γ-
Γ) augmente dans les deux cas (avec et sans tenir compte de l'effet de désordre
compositionnel). Les deux courbes présentent le même comportement, qui est non
monotone. Cependant, la non-linéarité semble être différente. En conséquence, on
s'attend bien entendu à des bowings optiques différents.
Une interpolation quadratique par la méthode des moindres carrés sur les deux
courbes nous donne les expressions analytiques de la forme:
EΓΓ x = −2.07x 6 + 3.24x + 2.36 VCA IV-4

6
EΓΓ x = 2.7x − 0.27x + 1.45 VCAA IV-5
En effet, lorsqu’on utilise la VCA seule, c'est-à-dire on a un facteur de désordre

nul (p=0), on a trouvé un bowing optique de -2,07 eV. Cette valeur n'est pas en accord
avec la valeur expérimentale recommandée par Vurgaftman et al. [4], cette valeur est
égale à 2.7 eV. Néanmoins, l'utilisation de la VCA amélioré, c'est-à-dire on tient compte
de l’effet du désordre comme un potentiel effectif dans le potentiel total du cristal 4 ≠
0 , on note un bowing optique de 2,7 eV, ce résultat est obtenu pour un facteur de
désordre égal à P = 0,973. Cette valeur du facteur de désordre indique que l’effet du
désordre compositionnel sur l’énergie de gap direct =>> pour alliage ternaire AlPA Sb5BA
est important comparé à d'autres alliages semiconducteurs ternaires III-V à savoir
l’alliage ternaire AlA In5BA Sb où le facteur de désordre presque négligeable avec une
valeur de 0.075 [9].
La Figure IV.5 présente la variation du gap d’énergie indirect =J

Γ
en fonction de
la composition x avec et sans effet du désordre.
53
2.6 AIPxSb1-x
2.4
2.2
Gap d'énergie (eV) 2.0
1.8
1.6
1.4
1.2
1.0
0.8
0.0 0.2 0.4 0.6 0.8 1.0

Figure IV.5 Variation du gap d’énergie indirect =>J en fonction de la composition x

avec et sans désordre de l’alliage ternaire AlPA Sb5BA
Les deux courbes sont interpolées par les deux équations suivantes :
=Γ? / = −0.64/ 6 + 1.55/ + 1.61 VCA IV-6

=Γ? 6
/ = 4.2/ − 3.24/ + 1.47 VCAA IV-7
La situation semble être différente pour le gap d'énergie indirect =>J de l’alliage
ternaire AlPA Sb5BA . Notons que, lorsqu’on utilise la VCA seule, la courbe varie d’une
façon monotone avec la composition de l'alliage x et avec un bowing optique de -0,64
eV. Cependant, on observe un comportement non monotone de =>J avec un bowing
optique élevé égal à 4,2 eV, et ceci lorsqu’on utilise la VCA amélioré. Encore une fois,
ceci montre la grande contribution de l'effet de désordre compositionnel sur les gaps
énergétiques de l’alliage considéré.
Les Figures (IV.6-IV.7) schématisent les variations des gaps d’énergies direct
=>> et indirect =>J et =>K en fonction de la composition x, avec et sans effet de désordre de
l’alliage ternaire AlPA Sb5BA .
54
4.0
AIPxSb1-x
3.5
Gap d'énergie (eV)

3.0
2.5
2.0 (Γ-Γ)
(Γ-X)
(Γ-L)
1.5
0.0 0.2 0.4 0.6 0.8 1.0
Figure IV.6 Variation des gaps d’énergies direct =>> et indirects et =>K et =>J en
fonction de la composition x sans désordre de l’alliage ternaire AlPA Sb5BA
4.0
AIPxSb1-x
3.5
3.0
Gap d'énergie (eV)
2.5
2.0
1.5
1.0
(Γ-Γ)
(Γ-X)
(Γ-L)
0.5
0.0 0.2 0.4 0.6 0.8 1.0
Figure IV.7 Variation des gaps d’énergies direct =>> et indirects et =>K et =>J en fonction
de la composition x avec désordre de l’alliage ternaire AlPA Sb5BA .
55
On remarque que dans les trois cas, et en passant du composé binaire 8 9: au

composé AlP les énergies de gaps direct=>> , et indirect=>J et =>K ne présentent pas le
même comportement en fonction de la variation de la composition x.
Selon notre étude, on peut conclure que l'absorption optiques des gaps dans AlPxSb1-x
est indirect pour toutes compositions x .
Pour examiner la transition énergétique les courbes du gap direct =>> et indirect
=>J et =>K de l’alliage ternaire AlPxSb1-x sont tracées en fonction de la composition x
sur les Figures (IV.6-IV.7) en utilisant la VCA seule et la VCA améliorée.
Nos résultats sur la Figure IV.6 montrent que lorsque l'effet de désordre
compositionnel est négligé, aucune transition entre les bandes passantes directes et
indirectes n’est possible, donc notre alliage est un semi-conducteur à bande interdite
(Γ-Χ) indirect le long de la composition x.
Néanmoins, la situation semble être complètement différente lorsque le VCA

améliorée est utilisée (Figure IV.7) où on note que le matériau est un semi-conducteur à
bande interdite indirecte (Γ- Χ) dans la gamme de composition d'alliage variant entre 0
et 0.076, puis on remarque qu’il y a une double transition de gap d’énergie. La première
transition vers le gap indirect (Γ- L) se fait à la composition x≈ 0,076 et qui correspond
à une énergie de 1,169 eV. La deuxième transition vers le gap indirect (Γ- X) se produit
à x≈ 0,24 et correspond à énergie de 0,886 eV. En conséquence, alliage ternaire
AlPxSb1-x peut être à gap indirect (Γ- Χ) ou indirect (Γ-L) en fonction de la teneur en P.
IV.1.3. Densité de charge électronique
La densité de charge dans les matériaux semi-conducteurs, est une caractéristique

très importante, elle nous renseigne sur la nature des liaisons cristallines dans le solide.
Son calcul dépend des fonctions d’ondes calculées à partir de la structure de bande. En
utilisant, dans notre cas, la méthode des pseudopotentiels (E.P.M).
La Figure IV.8 représente la densité de charge électronique au point Γ de la
somme des quatre bandes de valence selon la direction [111] du composé binaire AlSb.
56
Densité de charge électronique (units.arb)

25 AlSb
20
15
10
Al Sb
5
-0.50 -0.25 0.00 0.25 0.50
Position atomique (units.atom)
Figure IV. 8 Densité de charge électronique au point Γ de la somme des quatre

bandes de valence selon la direction [111] du composé binaire AlSb
Il existe une distribution asymétrique de charges où une quantité importante se

situe entre les sites atomiques. Ceci est dû aux liaisons partiellement covalente et
ionique qui caractérisent les semiconducteurs III-V.
Le maximum de la densité de charge de valence des électrons semble être
légèrement décalé vers l'anion (Sb), ceci indique que l’anion possède un potentiel plus
fort que le cation. Ainsi, la principale contribution à la formation de la liaison chimique
provient de l'anion. Une petite quantité de charges semble être dans les régions
interstitielles.
La Figure IV.9 représente la densité de charge électronique au point Γ de la

somme des quatre bandes de valence selon la direction [111], avec et sans tenir compte
de l’effet du désordre compositionnel, de l’alliage ternaire AlL.M InL.M Sb.
Quand la teneur des atomes du phosphore augmente dans le composé binaire

AlSb, le profil de la distribution de charge électronique de l’alliage ternaire AlL.M InL.M Sb
, en tenant compte de l’effet du désordre compositionnel, est identique. Néanmoins, le
maximum de la densité de charge de valence des électrons diminue en augmentant la
teneur en P dans AlPxSb1-x, affectant ainsi la contribution à la formation de la liaison
chimique. L'effet du désordre compositionnel sur la densité de charge de valence totale
d'électrons semble être significatif comme on peut l'observer sur la figure. Par
57
conséquent, cet effet devrait être pris dans tous les calculs de la distribution des charges
électroniques de la bande de valence de l’alliage ternaire AlPxSb1-x.
30
Densité de charge électronique (units. arb)

VCA
Al P0.5Sb0.5
VCAA
25
20
15
10
Al
5 (Sb,P)
-0.50 -0.25 0.00 0.25 0.50
Position atomique (units.atom )
Figure.I V.9 Densité de charge électronique au point Γ de la somme des quatre

bandes de valence selon la direction [111] de l’alliage ternaire AlL.M InL.M Sb .
Nous avons également calculé les densités de charge électronique de la première

bande de conduction au point Г selon la direction [111] pour le composé binaire AlSb,
et l’alliage ternaire AlPxSb1-x au point x=0.5. Nos résultats sont représentés sur les
Figure (IV.10-IV.11) respectivement. Notons que pour AlSb la majorité des charges
sont localisées autour du site de l'anion (Sb) où ils atteignent le maximum. Comparé à
l'anion, la quantité de charges autour du cation (Al) est moins prononcée.
Dans la région de liaison, nous observons que le minimum de la densité de charge
de conduction électronique est localisée presque à mi-chemin. La tendance semble être
similaire pour les alliages ternaires semi-conducteurs III-V [9]. Une faible quantité de
charge peut être observée dans les régions interstitielles, mais semble être plus
importante du côté du cation.
Lorsque la teneur du phosphore augmente dans le composé binaire AlSb, la

quantité de charges localisées autour du site de l’anion diminue alors que celles situées
autour du site cationique augmentent. Ainsi, il semble qu'il y ait un transfert de charges
de l'anion vers le cation. Dans la région de liaison, la situation semble être pratiquement
constante. L'effet du désordre compositionnel pour l’alliage ternaire AlPxSb1-x est de
58
nouveau significatif pour la densité de charge électronique de la première bande de

conduction.
16
Densité de charge atomique (units. arb )

AlSb
14
12
10
2
Al Sb
0
-0.50 -0.25 0.00 0.25 0.50
Position atomique (units. atom)

conduction selon la direction [1 1 1 ] du composé binaire AlSb.
Densité de charge électronique (units. arb)
16
VCA
VCAA
14 Al P Sb
0.5 0.5
12
10
2
Al (Sb,P)
0
-0.50 -0.25 0.00 0.25 0.50
Position atomique (units. atom )

conduction selon la direction [1 1 1 ] de l’alliage ternaire AlL.M InL.M Sb .
59
IV.2. Etude des propriétés optiques :
IV.2.1.Indice de réfraction
La connaissance de l'indice de réfraction (N) des semi-conducteurs est

indispensable pour des dispositifs tels que les cristaux photoniques, les guides d'ondes,
les cellules solaires et les détecteurs [10].
En tenant compte de l’effet du désordre, on calcule l’indice (N) selon différents

modèles pour les semiconducteurs AlSb, AlP0,5Sb0,5 et 8 P, les résultats sont dans le
Tableau IV.4
Tableau IV.4 Indice de réfraction n des semiconducteurs AlSb, et AlP0.50Sb0.50

AlP en utilisant différents modèles.
Matériau n calculé à partir de Expt.

Moss Gupta Hervé and Reddy and Ravindra Reddy and [11]
model and Vandamme Anjaneyulu model Ahammed
Ravindra model model model
model
AlSb 2.86 3.08 2.89 3.11 2.78 3.33 3.2327
AlP0.50Sb0.50 3.24 3.48 3.26 3.61 3.33 3.98 -
AlP 2.56 2.52 2.51 2.67 2.02 2.91 2.75
D’après la comparaison des valeurs calculées et des données expérimentales

rapportées dans les références [11], On peut dire que, pour AlSb et AlP, la valeur de (N)
déterminé en utilisant le modèle Reddyet anjaneyulu est en meilleur accord avec
l'expérience [11] que celles calculées à partir des autres modèles. Le modèle de Reddy
et Anjaneyulu a déjà validé ses bons résultats par rapport à l’expérience pour les semi-
conducteurs III-V [12].
La variation de l’indice de réfraction (N) en fonction de la composition x,

calculée à partir des différents modèles, est tracée sur la Figure IV.12. On constate que
pour tous les modèles étudiés, (N) décroit avec la concentration de phosphore. Le
comportement des courbes est monotone. Généralement c’est un trait commun pour la
plupart des semiconducteurs III-V [13],
60
où l'énergie minimale correspond à la plus grande valeur de n. La diminution ou

l’augmentation de l’indice de réfraction est directement lié à la valeur du gap en
fonction de x. Il diminue quand le gap augmente et vice versa.
4.5
AIPxSb1-x
4.0
Indice de réfraction
3.5
3.0
Moss Model
Gupta and Ravinda model
2.5 Hervé and Vandamme model
Reddy and Anjaneyulu model
Ravinda model
2.0 Reddy and Ahammed model
0.0 0.2 0.4 0.6 0.8 1.0

Figure.V.12 La variation de l’indice de réfraction en fonction de la composition x de l’alliage

ternaire AlPxSb1-x
IV.3. Etude des propriétés diélectriques :
IV.3.1. Constante diélectrique de haute fréquence :
La connaissance de (N) nous permet de calculer la constante diélectrique de

haute fréquence O∞ qui est considérée comme un paramètre important dans la
conception des dispositifs optiques. Les constantes diélectriques de haute fréquence des
composés binaires AlSb et AlP calculées sont 9,67 et 7,12 respectivement. Ces valeurs
sont comparées par rapport aux valeurs expérimentales qui sont de 9,88 (pour AlSb) et
de 7,4 (pour AlP) citées dans la Réf. [14], un accord à 2-4% est atteint.
Tableau IV.5 Constante diélectrique de haute fréquence (ε∞) des semiconducteurs

AlSb, AlP0.50Sb0.50, AlP.
Material ε∞
AlSb 9.67a
9.88b
61
AlP0.50Sb0.50 13
AlP 7.12a
7.4b
a
Present work
b
[ 14 ]
La variation de la constante diélectrique de haute fréquence en fonction de la

concentration du phosphore de l’alliage ternaire AlPxSb1-x est portée sur la Figure IV. 13.
14
Constante diélectrique de haute fréquence
AIPxSb1-x
13
12
11
10
0.0 0.2 0.4 0.6 0.8 1.0

Figure IV.13 La variation de la constante diélectrique de haute fréquence en fonction de

la composition x de l’alliage ternaire AlPxSb1-x (VCAA)
L’augmentation de la composition x en passant de 0 (AlSb) à 1 (AlP) conduit à

la décroissance non monotone. Ceci suggère que AlPxSb1-x devient progressivement
un bon isolant en augmentant la teneur en P. La tendance de O∞ est opposée à celle x et
qualitativement semblable à celle de N .
En raison de l'absence de données concernant la constante diélectrique de haute

fréquence OP pour AlP0.50Sb0.50 dans la littérature, notre résultat servira comme référence
pour des investigations futures.
62
IV.4. Conclusion
En résumé, notre étude théorique a été faite sur les propriétés électroniques et
optiques de l'alliage semi-conducteur ternaire III-V à base d’Aluminium cristallisant
dans la phase Zinc-blende à savoir : AlPxSb1-x . La méthode de calcul est la méthode des
pseudopotentiels empiriques (l'EPM) couplée avec la VCA améliorée qui prend en
compte l'effet du désordre de la composition x. les caractéristiques à savoir : la
structure de bande, la densité de charge électronique, l’indice de réfraction (n), la
constante diélectrique de haute fréquence OP ont été étudiés en fonction de la
concentration de Phosphore.
Nos résultats sont en bon accord avec d’autres résultats expérimentaux et

théoriques, sauf que par manque de références théoriques ou expérimentales, d’autres
résultats restent prédictifs ou peuvent servir de référence pour les recherches futures.
63
Références bibliographiques- chapitre IV
[1] C. Kittel, « Physique de l’état solide, Ed Dunod, Paris, (2001).
[2] M.L. Cohen, J.R. Chelikowsky, Electronic Structure and Optical Properties of
Semiconductors, Springer, Berlin, 1988.
[3] C. Alibert, A. Joullié, A. M. Joullié, C. Ance, Phys. Rev. B 27 (1983) 4946.
[4] I. Vurgaftman, J. R. Meyer, L. R. Ram-Mohan, J. Appl. Phys. 89 (2001) 5815.
[5] L. Vegard, Z. Phys. 5 (1921) 17.
[6] S.J.Lee, T.S. Kwon, K. Nahm, C.K. Kim, J. Phys. Condens. Matter, 2 (1990) 3253.
[7] N. Bouarissa, Phys. Lett. A 245 (1998) 285.
[8] J.R. Chelikowsky, T. J. Wagner A. Tin, J.H. Weaver, Phys. Rev. B40 (1989) 9644.
[9] N.E.-H. Fares, N. Bouarissa, Infrared Physics & Technology 71 (2015) 396.
[10] S. Adachi, J. Appl. Phys. 61 (1987) 4869.
[11] M. J. Weber, Handbook of Optical Materials, Springer, Berlin, 2003.
[12] M.S. Al-Assiri, N. Bouarissa, Superlatt. Microstruct. 59 (2013) 144.
[13] N. Bouarissa, Mater. Lett. 60 (2006) 2974.
[ 14 ]S. Adachi, Properties of Group IV, III–V and II–VI Semiconductors, Wiley, Chichester,
2005.
Conclusion Générale
Les composants électroniques mettent à profit les propriétés des électrons dans
les semiconducteurs, il est par conséquent nécessaire de préciser les propriétés
physiques de ces matériaux semiconducteurs, à savoir les propriétés électroniques
optiques et diélectriques qui conditionnent les caractéristiques de ces composants
électroniques.
Les calculs sont principalement effectués par l’emploi de la méthode du

peuseudopotentiel empirique EPM, couplée à l’approximation du cristal virtuel VCA,
où l'effet du désordre compositionnel est ajouté comme un potentiel effectif, cette
méthode a fait preuve d’une efficacité particulière pour l’estimation des gaps.
L‘étude a été faite sur les composés binaires ,etAlP ainsi que sur l’alliage ternaire
AlPxSb1-x . Les énergies de gap direct Γ Γ et indirect X Γ calculées par la méthode
des pseudopotentiels pour les différents semiconducteurs cité, coïncident parfaitement
avec les valeurs expérimentales et théoriques de la littérature .
D’après les tracés des structures de bandes électroniques des composés binaires
parents et AlPet de l’alliageAlPXSb1-X, on remarque clairement que l’effet du
désordre est très fort et qu’on peu pas le négliger .On remarque aussi que les quatre
bandes de valence sont moins dispersives que les quatre bandes de conduction et cela
est dûe à délocalisation des électrons de conductions dans leurs orbitales.
D’autres propriétés électroniques ont été calculées à partir de la structure de

bandes électroniques telle que la masse effective des électrons des composés binaires
AlPet AlSb . Ces valeurs sont en bon accord avec les valeurs expéromentale
.enaugmentant la concentration de phosphore on déduit que la mobilité des électrons
diminue.
La densité de charge et l’ionicité sont deux grandeurs fondamentales dans la

compréhension des liaisons chimiques dans les semiconducteurs, la densité a été
tracée au point Γ de la zone de Brillouin, on a remarqué que la distribution des
charges n’est pas équitablement répartie dans la région inter-atomique, on peut
conclure que les liaisons sont partiellement covalentes et ioniques.
La variation de l’indice de réfraction (n) en fonction de la concentration de
phosphor pour différents modèles a été clairement présentée. En comparant nos
résultats à ceux de la littérature, on remarque que le modèle de Reddy
etAnjaneyuludonne un très bon accord, et ceci nous a permet de conclure que la
diminution ou de l’indice de réfraction est étroitement lié à la valeur du gap
énergétique.
La connaissance de l’indice de réfraction nous a permis de déterminer la constante

diélectrique de haute fréquence ( ∞), qui est une valeur cruciale dans la conception
des dispositifs optoélectroniques. La variation de ( ∞) en fonction de la concentration
de phosphor nous a donne un comportement décroissant et monotone. Néanmoins en
absence de la référence expérimentale ou théorique, nos résultats restent prédictifs.
Certains résultats ont été parfaitement vérifiés par une comparaison avec les
valeurs expérimentales et théoriques et qui sont en bon accord. Par contre certains
résultats restent une prédiction par manque de références théoriques ou
expérimentales, et peuvent servir comme référence pour les recherches futures.
Finalement, quoique la méthode utilisée est simple par rapport aux méthodes du
premier principe dites (ab-initio), et qu’elle est plus économique vis à vis du temps de
calcul, elle donne des résultats en très bon accord avec d’autres résultats théoriques et
expérimentaux.

Tese Brhan

Transféré par

Droits d'auteur :

Formats disponibles

Tese Brhan

Transféré par

Informations du document

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Tese Brhan

Transféré par

Droits d'auteur :

Formats disponibles

République Algérienne Démocratique et Populaire

‫ــ" ا ـ! ـــ وا ـ ـــ ا ـ ــ ــ‬#‫ـ ـ ـ‬$ ‫وزارة ا‬

Mémoire de fin d’étude

Préparé par :Kolli Bourhan Eddine

Soutenu le : 2Juillet 2018

Président Dr Grar Nabila

• Je voudrais d’abord remercier « Allah » que DIEU nous bénisse ! Je tiens

Chapitre I Généralités sur les semi-conducteurs et leurs alliages

I.4.2. Gap énergétique

I.4.3. Densité de charge 9

Chapitre II Méthode de calcule des structures de bandes électroniques

Chapitre III Méthode des pseudopotentiels

Chapitre IV Résultats et Discussions

Figure 1.1. Maille élémentaire de la structure Zinc. Blende 6

Figure IV.6 Variation des gaps d’énergies direct et indirects et et en 55

Tableau.IV.1 Les facteurs de formes symétriques (Vs) et antisymétriques (VA) 48

Tableau.IV.4.indice de réfraction (n) des semi-conducteurs AlSb ,AlPxSb1-x et AlP 59

Toutefois ils sont par excellence les matériaux de l’optoélectronique permettant

Pour mieux comprendre la structure électronique des semi-conducteurs, il faut

L’étude du comportement des électrons dans un cristal dont le nombre des

Les techniques de calcul de la structure électronique mises au point au cours des

On débute par une introduction générale, où on expose la problématique de cette

Dans le premier chapitre on expose la théorie des semi-conducteurs et leurs

Une étude bien détaillée, en équations analytiques, où a été mise en évidence de

Finalement, notre travail est achevé par une conclusion générale.

[2] G. Roosen, « Matériaux semiconducteurs III-V, II-VI et nitrures pour

[4]J. Mycielshy and A. M.W. Itowski, Phys.Stat.Sol.(b) 134 (1986) 134

LES SEMICONDUCTEURS ET LEURS ALLIAGES

I.2. Notions générales sur les semiconducteurs

Les semi-conducteurs sont principalement les éléments de la colonne IV du tableau de

I.2.1. Semi-conducteur intrinsèque

Un semi-conducteur est dit intrinsèque, si ces propriétés électriques sont entièrement

I.2.2. Semi-conducteur extrinsèque

I-3 Propriétés structurales

L’immense majorité des semi-conducteurs utilisés dans les applications électroniques

I.3.1. La constante de réseau

Figure 1.2 Réseau triclinique primitif de l'espace tridimensionnel.

I.3.2. Première zone de Brillouin

La première zone de Brillouin est définie de manière unique comme la maille

La première zone de Brillouin est la cellule de Wigner-Seitz du réseau réciproque, on

I.3.2. 1. Les points de haute symétrie

: Centre de la première zone de Brillouin avec les coordonnées

: Centre d’une face hexagonale de l’octaèdre ayant les coordonnées =2 (1,1,1).

I.3.2.2. Les lignes de haute symétrie

∆ : Ligne de direction (100), reliant le centre Γ au point X.

Σ : Point appartenant au plan de symétrie kx = ky ou ky = kz ou kx = kz.

I.4. Propriétés électroniques

I.4.1. La structure de bande

La structure de bandes est généralement identique pour les composés (III.V).

I.4.2. Gap énergétique

Figure I.4 a. Semi-conducteur à gap directb. Semi-conducteur à gap indirect

I.4.3. Densité de charge

La densité de charge électronique est un moyen pour l’analyse et l’interprétation d’un

I.5. Propriétés optiques

La connaissance de l'indice de réfraction (n) des semi-conducteurs est indispensable

(i) Modèle de Ravindra et Srivastava [9]

n" E$ = 108 I−2

(ii) Modèle de Gupta et Ravindra [10]

n = 4.084 − 0.62E$ I−3