Hétérostructure GaN/AlGaN à puits quantiques : guide pratique sur le modèle à 8 bandes de Wurtzite.
1. Introduction
Dans ce tutoriel, nous allons calculer les états électroniques et la réponse optique d’un puits quantique en GaN, situé entre deux barrières en AlGaN, en utilisant le modèle wurtzite à huit bandes de OghmaNano. L’exemple comprend un puits en GaN de 10 nm et deux couches d’AlGaN de 20 nm. Nous allons d’abord exécuter l’exemple fourni, puis examiner les dispersions des sous-bandes de conduction et de valence, ainsi que les spectres optiques obtenus.
La distinction entre la structure cristalline de la wurtzite et celle de la zincblende est importante ici. Les deux sont des structures cristallines semi-conductrices, mais la wurtzite présente un axe c hexagonal privilégié. La différence de symétrie cristalline modifie la séparation des bandes de valence, le couplage entre les états électroniques et l’intensité de la lumière polarisée dans différentes directions. Par conséquent, un calcul à huit bandes basé sur la structure de la zincblende ne peut pas être simplement réutilisé avec les paramètres des matériaux du GaN.
2. Ouverture de l’exemple de la structure wurtzite à huit bandes.
Dans le ruban « Fichier », cliquez sur « Nouvelle simulation ». Dans la fenêtre qui s’affiche, double-cliquez sur « Démonstrations GaAs » (??). Malgré le nom de la catégorie, elle contient également des exemples de puits quantiques à base de nitrure et de silicium. Dans la fenêtre suivante, double-cliquez sur « GaN Wurtzite 8 band k.p » (??). Cela ouvre la pile de matériaux préparée et sélectionne la démonstration appropriée.
3. Inspection de la puits quantique GaN/AlGaN.
La fenêtre principale affiche la structure à trois couches (??). Dans le calcul réel d’un puits quantique, le confinement se produit le long de l’épaisseur de cette structure : les fonctions d’enveloppe électroniques varient d’une couche à l’autre, tandis que les diagrammes de dispersion décrivent le mouvement parallèle à ces couches.
Cliquez sur « Éditeur de calques » dans la barre d’outils située à gauche. Comme illustré dans ??, l’exemple utilise la géométrie suivante :
| Couche | Épaisseur | Type |
|---|---|---|
| AlGaN | 20 nanomètres | Actif / barrière |
| GaN (puits quantiques) | 10 milles nautiques | puits quantique |
| AlGaN | 20 nanomètres | Actif / barrière |
Pour votre première simulation, ne modifiez pas les paramètres des matériaux fournis. La capture d’écran indique l’épaisseur des couches, mais ne montre pas la fraction molaire d’aluminium, la configuration de contrainte ou les paramètres de polarisation. Ces éléments doivent être vérifiés dans les paramètres des matériaux de l’exemple avant de l’utiliser pour représenter une plaquette expérimentale spécifique.
4. Pourquoi utiliser un hamiltonien de Wurtzite à huit bandes ?
Un modèle parabolique simple attribue une masse effective à chaque porteur de charge et suppose que son énergie varie approximativement en fonction du carré de son vecteur d’onde dans le plan. Ce modèle est utile à proximité d’un point extrême de la bande d’énergie, mais il ne permet pas de décrire la manière dont les différents états de valence interagissent les uns avec les autres lorsque le vecteur d’onde augmente. Dans les matériaux de type wurtzite, l’environnement cristallin hexagonal distingue également le mouvement et la polarisation de la lumière le long de l’axe c de ceux observés dans le plan perpendiculaire.
L’approche à huit bandes conserve une orbitale de type conducteur et trois orbitales de type valence, chacune ayant deux composantes de spin. Conceptuellement, un état propre a la forme suivante :
\[ \Psi_{n,\mathbf{k}_{\parallel}}(\mathbf{r})= e^{i\mathbf{k}_{\parallel}\cdot\mathbf{r}_{\parallel}} \sum_{\alpha=1}^{8}F_{n\alpha}(z)u_{\alpha}(\mathbf{r}), \]
où les fonctions \(F_{n\alpha}(z)\) décrivent le confinement dans le puits quantique, et \(u_\alpha\) représentent les huit états de base de Bloch. Un état confiné donné peut donc être constitué de plusieurs orbitales de la bande de valence ; c’est pourquoi il est trompeur de décrire chaque courbe colorée à l’aide d’une seule étiquette fixe, correspondant à un trou lourd ou un trou léger, lorsqu’on s’éloigne du bord de la bande.
L’ordre des bandes de valence de la wurtzite est influencé par les interactions du champ cristallin et de l’interaction spin-orbite. En utilisant une notation concise, le calcul résout
\[ H_{8\times8}\!\left(k_x,k_y,-i\frac{\partial}{\partial z}\right) \mathbf{F}_n(z)=E_n(\mathbf{k}_{\parallel})\,\mathbf{F}_n(z). \]
Remarquez que le vecteur d’onde dans le plan est conservé explicitement, tandis que le vecteur d’onde correspondant à la direction de croissance est remplacé par un opérateur différentiel. En répétant ce problème de valeurs propres pour différentes valeurs, \(k_{\parallel}\) génère les graphiques de dispersion des sous-bandes présentés ci-dessous. Une analyse plus approfondie des paramètres de base, des termes du champ cristallin et des coefficients des matériaux est fournie dans le Chapitre sur la théorie k·p de la wurtzite..
5. Lancez la simulation.
Retournez à l’onglet « Fichier » et cliquez sur « Exécuter la simulation ». Ouvrez l’onglet « Terminal » pour suivre l’avancement. Dans ??, le terminal identifie le solveur WZ KP8, signale l’assemblage de la matrice creuse, affiche les informations sur la convergence et effectue les calculs pour chaque vecteur d’onde échantillonné. Il résout séparément les états électroniques et de valence à chaque point.
Étant donné qu’il s’agit d’un problème à plusieurs composantes impliquant une fonction d’enveloppe, il est plus complexe que la résolution d’une équation parabolique scalaire indépendante pour chaque bande. Attendez que le solveur d’équations propres et les calculs optiques qui suivent soient terminés avant d’examiner les fichiers résultants.
5.1 Vérification des valeurs des transitions optiques.
Plus tard dans la simulation, le programme affiche les valeurs des éléments de la matrice optique (??). Pour une structure conventionnelle de type wurtzite, orientée selon le plan c, TE désigne le champ électrique optique perpendiculaire à l’axe c, tandis que TM désigne sa composante parallèle à cet axe. Ces polarisations n’ont pas nécessairement la même influence sur une transition électronique donnée.
Une façon utile de concevoir la force optique est de considérer l’élément carré de la matrice de transition :
\[ M_{cv}^{(\mathbf{e})}(\mathbf{k}_{\parallel})= \left\langle\Psi_{c,\mathbf{k}_{\parallel}}\middle| \mathbf{e}\!\cdot\!\mathbf{p}\middle| \Psi_{v,\mathbf{k}_{\parallel}}\right\rangle, \qquad S_{cv}^{(\mathbf{e})}\propto\left|M_{cv}^{(\mathbf{e})}\right|^2. \]
Ici, \(\mathbf{e}\) désigne la polarisation de la lumière et \(\mathbf{p}\) est l’opérateur d’impulsion. Un état peut exister dans la structure de bande électronique, mais contribuer très peu à une polarisation optique particulière si son élément de matrice est faible. Les valeurs relatives TE/TM à la borne fournissent donc des informations physiques supplémentaires, au-delà de la simple position des niveaux d’énergie.
6. Interprétation des diagrammes de bandes électroniques calculés.
L’axe horizontal de chaque graphique représente le vecteur d’onde dans le plan \(k_{\parallel}\), et non la position à l’intérieur du puits quantique. L’axe vertical représente l’énergie d’une sous-bande confinée. Au point \(k_{\parallel}=0\), chaque courbe commence à l’une des énergies de confinement discrètes. En se déplaçant vers la droite, on observe comment le même état électronique se disperse lorsque la charge se déplace parallèlement aux couches.
6.1 États de la bande de conduction
Ouvrez la sortie de la bande de conduction et examinez la dispersion dans ??. Plusieurs niveaux confinés sont visibles, en commençant par différentes énergies proches du bord de la bande de GaN. Pour une faible valeur du vecteur d’onde, une approximation utile pour une sous-bande de conduction est :
\[ E_{c,n}(k_{\parallel})\simeq E_{c,n}(0) +\frac{\hbar^2 k_{\parallel}^2}{2m_{e,\parallel}^{*}}. \]
Le terme \(E_{c,n}(0)\) contient l’énergie de confinement, tandis que \(m_{e,\parallel}^{*}\) décrit la réponse du support aux mouvements dans le plan. Cette relation permet d’expliquer les courbes qui présentent une tendance à la hausse, bien que le calcul complet sur les huit bandes conserve le couplage inter-bandes plutôt que d’imposer une forme parabolique fixe.
6.2 États de la bande de valence
Le diagramme de la bande de valence (??) est plus complexe. Sa branche supérieure s’élève initialement en s’éloignant de \(k_{\parallel}=0\) avant de s’incurver vers le bas ; plusieurs branches inférieures présentent des courbures visiblement différentes. De telles formes sont possibles dans un modèle de bande de valence de type wurtzite couplé, dans lequel la séparation du champ cristallin, le couplage spin-orbite, le confinement et le mélange des bandes agissent de concert.
Il ne faut pas interpréter cette représentation comme une preuve que chaque courbe colorée possède une identité fixe dans les bandes A, B ou C sur l’ensemble de la plage de vecteurs d’onde. Pour établir une attribution définitive, il convient d’examiner les contributions des composantes du vecteur propre sous-jacent. De même, il est intéressant de comparer les valeurs extrêmes inhabituelles et non centrées avec les paramètres de matériau et de déformation choisis, si l’on utilise le modèle de manière quantitative.
7. Transitions optiques, émission stimulée et absorption.
Les bandes électroniques déterminent les énergies des transitions optiques possibles, mais l’intensité de ces transitions dépend également des fonctions d’onde des électrons et des trous, de la polarisation et de l’occupation des porteurs de charge. Pour un vecteur d’onde commun dans le plan, l’énergie du photon pour une transition verticale est de
\[ \hbar\omega=E_{c,n}(\mathbf{k}_{\parallel}) -E_{v,m}(\mathbf{k}_{\parallel}). \]
La règle d’or de Fermi explique pourquoi les éléments de la matrice optique affichés à la fin du processus sont importants. Sous forme schématique, le taux d’une transition est le suivant :
\[ W_{v\rightarrow c}^{(\mathbf{e})}= \frac{2\pi}{\hbar}\left|H_{cv}^{\prime(\mathbf{e})}\right|^2 \delta\!\left(E_c-E_v-\hbar\omega\right). \]
L’élément de la matrice de perturbation \(H'_{cv}\) dépend de la polarisation de la lumière et des états impliqués. La fonction delta impose la conservation de l’énergie ; lors de la construction d’un spectre pratique, les contributions des transitions autorisées sont additionnées pour les états échantillonnés, puis élargies. Les facteurs d’occupation des porteurs déterminent ensuite si l’absorption ou l’émission stimulée prédomine à une longueur d’onde donnée. Un facteur net schématique pratique est \(f_c-f_v\), où \(f_c\) et \(f_v\) représentent respectivement l’occupation des états de conduction et de valence par les électrons.
Étant donné que la couche de GaN présente une large bande interdite, ses spectres calculés apparaissent dans le domaine proche de l’ultraviolet, principalement autour de 330 à 370 nm dans l’exemple fourni. La longueur d’onde et l’énergie du photon sont liées par la formule \(E_{\gamma}=hc/\lambda\).
7.1 Émission stimulée et gain net
Le graphique de l’émission stimulée (??) présente une série de spectres correspondant à différentes conditions de fonctionnement. Le calcul du gain net (??) prend en compte la compétition entre l’émission et l’absorption. Un gain matériel positif indique une amplification au sein de ce modèle électronique, mais il ne prouve pas à lui seul qu’une cavité laser complète a atteint le seuil : il faut également tenir compte du facteur de confinement et des pertes du résonateur.
7.2 Absorption
Les résultats d’absorption sont présentés dans ??. Les contours et les caractéristiques spectrales sont déterminés par les écarts entre les niveaux de valence occupés et les niveaux de conduction disponibles, pondérés par leurs éléments de matrice optique. Dans une hétérostructure de type wurtzite, les spectres TE et TM correspondants peuvent différer, car la composition orbitale d’un niveau détermine l’intensité de son couplage avec chaque polarisation.
Les caractéristiques observées dans ces spectres sont le résultat du modèle configuré et de l’élargissement spectral. Elles ne doivent pas être automatiquement associées à des pics excitoniques spécifiques : un calcul de transition interbande du type k·p pour une seule particule ne constitue pas, en soi, une description complète des excitons liés par l’interaction de Coulomb entre les électrons et les trous.
8. Émission spontanée
Contrairement à l’émission stimulée, l’émission spontanée se produit sans qu’un photon incident ne déclenche la transition. Elle dépend de la présence d’un électron dans l’état de conduction et d’un état vide dans la bande de valence (un trou). Son facteur d’occupation électronique est donc proportionnel à
\[ R_{\mathrm{sp}}(\omega)\ \propto\ \sum_{c,v,\mathbf{k}_{\parallel}} \left|M_{cv}(\mathbf{k}_{\parallel})\right|^2 f_c(1-f_v)\,L\!\left(E_c-E_v-\hbar\omega\right), \]
où \(L\) désigne l’élargissement spectral choisi, et la relation de proportionnalité ne tient pas compte des préfacteurs photoniques. Le résultat est représenté graphiquement dans ??. L’échelle logarithmique permet de comparer des spectres couvrant plusieurs ordres de grandeur. Remarquez les caractéristiques optiques près de la partie du spectre correspondant aux longueurs d’onde les plus courtes, ainsi que la variation de l’amplitude globale dans les différentes conditions représentées.
9. Essayez-le vous-même.
💡 Deux calculs complémentaires utiles :
- Modifiez la largeur de la couche de GaN. Essayez de réduire la couche centrale de 10 nm à 7 nm, en laissant les barrières inchangées. Recalculez les états propres et comparez les niveaux de conduction les plus bas et les énergies de transition optique. Un confinement plus fort augmentera généralement la séparation entre les niveaux confinés, bien que la réponse détaillée de la bande de valence soit multi-bande.
- Comparez les intensités des modes TE et TM. Utilisez les intensités optiques mesurées par le terminal pour identifier les transitions qui favorisent fortement une polarisation particulière. Analysez comment cela se rapporte au caractère de l’état de valence, plutôt que de supposer que toutes les transitions contribuent de manière égale au spectre.
✅ Quels sont les éléments à prendre en compte ?
En général, un puits plus étroit a tendance à augmenter les niveaux d’énergie de confinement des électrons et à modifier les niveaux des trous, ce qui entraîne un décalage des énergies de transition interbande disponibles. Les rapports TE/TM peuvent varier considérablement d’une paire d’états à l’autre, car ils sondent différentes composantes orbitales des fonctions d’onde de la structure cristalline wurtzite. Pour une étude quantitative des matériaux, il est également important d’examiner les effets de la composition de l’alliage, de la contrainte et des champs de polarisation.