Voir la source de Théorie k·p

En [[physique du solide]], la '''théorie k·p''' est basée sur la [[théorie de la perturbation (mécanique quantique)|théorie de perturbation]] de la mécanique quantique et est une méthode empirique utilisée pour calculer la [[structure de bande]] et les propriétés optiques des solides cristallins<ref name=Yu2.6/>{{,}}<ref name=Kittel/>{{,}}<ref name=Harrison>
{{Ouvrage|langue=en|auteur1=W.P. Harrison|titre=Electronic Structure and the Properties of Solids|sous-titre=the physics of the chemical bond|lieu=New York|éditeur=[[Dover Publications]]|année=1989|numéro d'édition=Reprint|année première édition=1980|pages totales=158 ''ff''|isbn=0-486-66021-4}}</ref>.
Cette théorie a été appliquée notamment dans le contexte du modèle de Lüttinger-Kohn (d'après [[Joaquin Luttinger]] et [[Walter Kohn]]), et du modèle de Kane (d'après [[Evan Kane]]). Elle a également donné naissance au modèle de [[masse effective]] très fréquemment utilisé en physique du solide.

==Dérivation==

===Théorème de Bloch et vecteurs d'onde===
{{voir|Onde de Bloch}}

Selon la [[mécanique quantique]] (dans l'[[méthode de Hartree–Fock|approximation des électrons indépendants]]), la [[fonction d'onde]] des [[électrons]] dans un matériau peut être décrite par l'[[équation de Schrödinger]] :
:<math>\left(\frac{p^2}{2m}+V\right)\psi = E\psi</math>
où '''p''' est l'[[quantité de mouvement|opérateur de quantité de mouvement]], ''V'' est le [[potentiel électrostatique]], et ''m'' est la masse d'un électron (et où l'effet du [[interaction spin-orbite|couplage spin-orbite]] a été négligé).

Dans un [[cristal]], ''V'' est une [[fonction périodique]], qui a la périodicité du [[réseau cristallin]]. Le [[théorème de Bloch]] prouve que les solutions de cette [[équation différentielle]] peuvent être écrites sous la forme :
:<math>\psi_{n,\mathbf{k}}(\mathbf{r}) = e^{i\mathbf{k}\cdot\mathbf{r}} u_{n,\mathbf{k}}(\mathbf{r})</math>
où '''k''' est le [[vecteur d'onde]], ''n'' est l'indice discret de la [[structure de bande|bande]], et ''u''<sub>''n'','''k'''</sub> est une fonction qui a la même périodicité que le cristal, c'est-à-dire qui satisfait :
:<math>u_{n,\mathbf{k}}(\mathbf{r} + \mathbf{R}) = u_{n,\mathbf{k}}(\mathbf{r})</math>
où '''R''' est un vecteur du [[réseau de Bravais]].

Pour chaque bande (indice discret ''n''), il y a une relation entre le vecteur d'onde '''k''' et l'énergie de l'état ''E''<sub>''n'','''k'''</sub>, appelée [[relation de dispersion]]. Le calcul de cette relation de dispersion (ou structure de bande) est l'application principale de la théorie '''k·p'''.

===Théorie de la perturbation===
{{voir|Théorie de la perturbation (mécanique quantique)}}

La fonction périodique ''u''<sub>''n'','''k'''</sub> satisfait l'équation de Schrödinger réécrite de la manière suivante<ref name=Yu2.6>
{{Ouvrage|langue=en|auteur1=P. Yu|auteur2=M. Cardona|titre=Fundamentals of Semiconductors|sous-titre=Physics and Materials Properties|éditeur=[[Springer (publisher)|Springer]]|année=2005|numéro d'édition=3|passage=Section 2.6, pp. 68 ''ff'|isbn=3-540-25470-6|lire en ligne=https://books.google.com/books?id=W9pdJZoAeyEC&pg=PA244&dq=isbn%3D3540254706}}</ref> :
:<math>H_{\mathbf{k}} u_{n,\mathbf{k}}=E_{n,\mathbf{k}}u_{n,\mathbf{k}}</math>
où l'[[Opérateur hamiltonien|Hamiltonien]] est donné par : 
:<math>H_{\mathbf{k}} = \frac{p^2}{2m} + \frac{\hbar \mathbf{k}\cdot\mathbf{p}}{m} + \frac{\hbar^2 k^2}{2m}  + V </math>
Ici, '''k''' est un vecteur constitué de trois nombres réels dont l'unité est l'inverse d'une distance, alors que '''p''' est un vecteur d'opérateurs ; plus précisément,
:<math>\mathbf{k}\cdot\mathbf{p} = k_x (-i\hbar \frac{\partial}{\partial x}) + k_y (-i\hbar \frac{\partial}{\partial y}) + k_z (-i\hbar \frac{\partial}{\partial z})</math>
On peut écrire cet hamiltonien comme une somme de deux termes :
:<math>H=H_0+H_{\mathbf{k}}', \;\; H_0 = \frac{p^2}{2m}+V, \;\; H_{\mathbf{k}}' = \frac{\hbar^2 k^2}{2m} + \frac{\hbar \mathbf{k}\cdot\mathbf{p}}{m}</math>
Cette expression est à la base de la [[Théorie de la perturbation (mécanique quantique)|perturbation]]. L'hamiltonien "non-perturbé" est ''H''<sub>0</sub>, qui est en fait l'hamiltonien exact à '''k'''=0 (c'est-à-dire en Γ). La "perturbation" dans ce cas est donc le terme <math>H_{\mathbf{k}}'</math>. La méthode résultante est finalement appelée "théorie '''k·p'''" (ou "méthode '''k·p'''" ou "théorie de perturbation '''k·p'''"), en référence du terme proportionnel à ''k''·''p'' qui apparaît. Le résultat de cette analyse est une expression pour ''E''<sub>''n'','''k'''</sub> et ''u''<sub>''n'','''k'''</sub> en termes des énergies et fonctions d'ondes en '''k'''=0.

Le terme de "perturbation" <math>H_{\mathbf{k}}'</math> devient progressivement plus petit lorsque '''k''' s'approche de zéro. Pour cette raison, la théorie '''k·p''' est général adaptée et plus précise pour des valeurs de '''k''' petites (ou proches de point de référence). Cependant, si suffisamment de termes sont inclus dans l'expansion de la perturbation, alors la précision du calcul peut être étendue à toute la [[zone de Brillouin]].

===Expression pour une bande non-dégénérée===

Pour une bande non-dégénérée (c'est-à-dire une bande qui a une énergie en '''k'''=0 différente de toute autre bande), avec un [[extremum]] à '''k'''=0, et sans [[interaction spin-orbite]], le résultat de la théorie ''k''·''p'' est (pour le  [[Théorie de la perturbation (mécanique quantique)|premier ordre non-trivial]])<ref name=Yu2.6/> :
:<math>u_{n,\mathbf{k}} = u_{n,0}+\frac{\hbar}{m}\sum_{n' \neq n}\frac{\langle u_{n,0} | \mathbf{k}\cdot\mathbf{p} | u_{n',0} \rangle}{E_{n,0}-E_{n',0}} u_{n',0}</math>
:<math>E_{n,\mathbf{k}} = E_{n,0}+\frac{\hbar^2 k^2}{2m} + \frac{\hbar^2}{m^2} \sum_{n'\neq n} \frac{|\langle u_{n,0} | \mathbf{k}\cdot\mathbf{p} | u_{n',0} \rangle |^2}{E_{n,0}-E_{n',0}}</math>

Comme '''k''' est un vecteur de nombres réels (plutôt qu'un vecteur d'opérateurs linéaires plus compliqués), les éléments de matrice de ces expressions peuvent être réécrites :
:<math>\langle u_{n,0} | \mathbf{k}\cdot\mathbf{p} | u_{n',0} \rangle = \mathbf{k} \cdot \langle u_{n,0} | \mathbf{p} | u_{n',0} \rangle</math>
C'est-à-dire qu'on peut calculer l'énergie de ''tout'' point '''k''' au moyen d'un petit nombre de paramètres empiriques, ''E''<sub>''n'',0</sub> et <math>\langle u_{n,0} | \mathbf{p} | u_{n',0} \rangle</math>. Ces derniers étant appelés "éléments de matrice optique", intimement reliés aux transitions du [[Dipôle électrostatique#Moment dipolaire|moments dipolaire électrique]]. Ces paramètres sont typiquement extraits de données expérimentales.

En pratique, la somme sur les bandes '''n''' n'inclut bien souvent que les bandes les plus proches, puisque ce sont celles qui contribuent le plus au calcul (étant donné le dénominateur). Cependant, pour étendre la précision du calcul, particulièrement pour des '''k''' éloignés, les bandes additionnelles doivent être incluses dans l'expansion ci-dessus.

====Masse effective====
{{Voir|Masse effective}}

D'après l'expression dérivée ci-dessus pour la relation de dispersion, une expression simplifiée pour la masse effective de la bande de conduction d'un semi-conducteur peut être dérivée<ref name=Harrison/>. Pour cela, l'énergie ''E<sub>n0</sub>'' est fixée au minimum de la bande de conduction ''E<sub>c0</sub>'' et seules les énergies proches du maximum de la bande de valence sont incluses dans la somme, pour lesquelles la différence d'énergie au dénominateur est la plus petite. Le dénominateur est alors équivalent à la valeur de la bande interdite ''E<sub>g</sub>'', donnant la relation de dispersion :
:<math>E_c(\boldsymbol k ) \approx E_{c0} +\frac{(\hbar k)^2}{2m} +\frac{\hbar ^2}{{E_g}m^2}\sum_n  {|\langle u_{c,0}|\mathbf{k}\cdot\mathbf{p}| u_{n,0} \rangle |^2}  </math>
La masse effective dans la direction ℓ est alors :
:<math> \frac{1} {m}_{\ell} = {{1} \over {\hbar^2}} \sum_{ m} \cdot {{\partial^{\ 2} E_{c} (\boldsymbol{k})} \over {\partial k_{\ell} \partial k_ m}} \approx \frac{1}{m}+\frac{2}{E_gm^2}\sum_{m,\ n} {\langle u_{c,0}|p_{\ell}| u_{n,0} \rangle }{\langle u_{n,0}|p_{m}| u_{c,0} \rangle }  </math>
En ignorant les détails des éléments de matrice, les conséquences clés de cette dérivation sont que la masse effective varie comme la bande interdite du semi-conducteur et tend vers zéro lorsque la bande interdite tend vers zéro<ref name=Harrison/>.
Une approximation utile des éléments de matrice dans un semi-conducteur à [[Théorie des bandes|bande interdite directe]] est<ref name=DirectGap> Un semi-conducteur à bande interdite ''directe'' est un semi-conducteur où le maximum de la bande de valence et le minimum de la bande de conduction sont à la même position dans l'espace des '''k''', la plupart du temps au point Γ ('''k''' = 0), tel que l'[[arséniure de gallium]] par exemple.</ref> :
:<math>\frac{2}{E_gm^2}\sum_{m,\ n} {|\langle u_{c,0}|p_{\ell}| u_{n,0} \rangle |}{|\langle u_{c,0}|p_{m}| u_{n,0} \rangle |} \approx 20\mathrm{eV} \frac{1}{mE_{g}} \ , </math>

Au contraire de la bande de conduction, les interactions ''spin-orbite'' doivent être introduites afin de décrire les énergies de la bande de valence et plusieurs bandes doivent être considérées individuellement.
Le calcul est donné dans le livre de Yu et Cardona<ref name=valence>See Yu & Cardona, ''op. cit.'' pp. 75-82</ref>. Dans ce cas, les porteurs libres sont des ''[[trou d'électron|trous]]'' et sont divisés en deux catégories : les trous ''légers'' et les trous ''lourds'', dont les masses sont anisotropes.

===Théorie k·p avec interaction spin-orbite===

Lorsque l'[[interaction spin-orbite]] est incluse, l'équation de Schrödinger pour ''u'' est<ref name=Kittel>
{{Ouvrage|langue=en|auteur1=C. Kittel|titre=Quantum Theory of Solids|lieu=New York|éditeur=[[John Wiley & Sons|Wiley]]|année=1987|numéro d'édition=Second Revised Printing|pages totales=186–190|isbn=0-471-62412-8}}</ref>:
:<math>H_{\mathbf{k}} u_{n,\mathbf{k}}=E_{n,\mathbf{k}}u_{n,\mathbf{k}}</math>
où<ref name="a">
{{article|langue=en
 |auteur=Evan O. Kane
 |année=1957
 |titre=Band Structure of Indium Antimonide
 |journal=[[Journal of Physics and Chemistry of Solids]]
 |volume=1 |pages=249
 |doi=10.1016/0022-3697(57)90013-6
|bibcode = 1957JPCS....1..249K }}</ref>
:<math>H_{\mathbf{k}} = \frac{p^2}{2m} + \frac{\hbar}{m}\mathbf{k}\cdot\mathbf{p} + \frac{\hbar^2 k^2}{2m}  + V + \frac{\hbar}{4 m^2 c^2} (\nabla V \times (\mathbf{p}+\hbar\mathbf{k}))\cdot \vec \sigma </math>
où <math>\vec \sigma=(\sigma_x,\sigma_y,\sigma_z)</math> est un vecteur constitué de trois [[matrices de Pauli]]. Cet hamiltonien peut être sujet à la même analyse issue de la théorie de la perturbation que celle présentée plus haut.

===Calcul dans le cas dégénéré===

Dans le cas de bandes dégénérée ou quasi-dégénérée, en particulier la [[bande de valence]] de la plupart des semi-conducteurs, les équations peuvent être analysées par les méthodes de la [[théorie de la perturbation (mécanique quantique)|théorie de perturbation dégénérée]]<ref name=Yu2.6/>{{,}}<ref name=Kittel/>. Les modèles de ce type sont notamment le modèle à 3 bandes de "[[Luttinger-Kohn]]"<ref>
{{article|langue=en
 |auteur=J. M. Luttinger, W. Kohn
 |année=1955
 |titre=Motion of Electrons and Holes in Perturbed Periodic Fields
 |journal=[[Physical Review]]
 |volume=97 |numéro= |pages=869
 |doi=10.1103/PhysRev.97.869
|bibcode = 1955PhRv...97..869L }}</ref> et de "[[Dresselhaus-Kip-Kittel (DKK)]]"<ref>
{{article|langue=en
 |auteur=G. Dresselhaus, A. F. Kip, C. Kittel
 |année=1955
 |titre=Cyclotron Resonance of Electrons and Holes in Silicon and Germanium Crystals
 |journal=[[Physical Review]]
 |volume=98 |pages=368
 |doi=10.1103/PhysRev.98.368
|bibcode = 1955PhRv....98..368K }}</ref>. Le modèle de [[Evan Kane|Kane]] inclut en plus la bande de conduction dans sa dérivation, ainsi que l'interaction entre bande de valence et bande de conduction <ref name="a" />.

==Voir aussi==
{{Col-début}}
'''Structure de bande électronique'''
* [[Structure de bande]]
* [[Théorie des bandes]]
* [[Pseudo-potentiel]]
'''Propriétés des semi-conducteurs'''
* [[Bande interdite]]
* [[Masse effective]]
* [[Densité d'états électroniques]]
* [[Surface de Fermi]]
'''Fonctions d'onde'''
* [[Onde de Bloch]]
* [[Fonction de Wannier]]
'''Théorie fondamentale'''
* [[Théorie de la fonctionnelle de la densité]]
* [[Théorie des groupes]]
{{col-fin}}

==Notes et références==
{{Traduction/Référence|en|k·p perturbation theory|746596248}}
{{Références}}

{{Palette|Calcul de la structure électronique}}
{{Portail|physique}}

[[Catégorie:Physique du solide]]