Hétérostructure GaAs/AlGaAs à puits quantiques : guide pratique sur le modèle à 8 bandes de type zincblende.

1. Introduction

Dans ce tutoriel, nous allons calculer les propriétés électroniques et optiques d’un puits quantique GaAs/AlGaAs en utilisant le modèle à huit bandes de type zincblende k·p d’OghmaNano. Nous commencerons par un exemple préétabli, examinerons la structure des couches, lancerons la simulation et analyserons les sous-bandes de conduction et de valence obtenues. Enfin, nous examinerons les spectres optiques calculés à partir des états.

Le tutoriel sur les puits quantiques paraboliques il traite les électrons et les trous en utilisant des bandes de masse effective relativement simples. Ici, l’hamiltonien comprend une base de bande de conduction et les bases de bande de valence pour les trous lourds (HH), les trous légers (LH) et les états de séparation (SO), avec deux composantes de spin pour chacun : huit états de base au total. Le couplage entre ces états permet au modèle de décrire des effets que des paraboles indépendantes ne peuvent pas, en particulier le mélange des bandes de valence et les changements de courbure en fonction du vecteur d’onde dans le plan.

2. Ouvrir l’exemple à huit bandes.

Dans OghmaNano, cliquez sur « Nouvelle simulation » dans le ruban « Fichier ». Double-cliquez sur le dossier « Démonstrations GaAs » (??). Parmi les exemples de puits quantiques, double-cliquez sur « GaAs/AlGaAs QW 8 band k.p » (??). Cela ouvre une simulation préconfigurée ; il n’est pas nécessaire de saisir manuellement une matrice hamiltonienne 8 x 8.

3. Structure à puits quantiques et paramètres des matériaux.

La fenêtre principale affiche une structure à trois couches : AlGaAs/GaAs/AlGaAs (??). Ouvrez l’éditeur de couches à l’aide du bouton situé à gauche. L’exemple fourni consiste en deux barrières d’AlGaAs de 40 nm encadrant une couche de GaAs de 10 nm :

Les limites de la bande de GaAs se situent à l’intérieur de l’intervalle de la bande d’AlGaAs environnante, ce qui permet de confiner les électrons dans la direction de croissance de la couche. Dans l’éditeur de couches, la couche centrale de GaAs est spécifiquement désignée comme une hétérostructure quantique (??). Les paramètres relatifs aux matériaux optiques et aux types de couches sont déjà configurés dans cet exemple.

Le solveur à huit bandes nécessite plus d’informations sur les matériaux qu’un modèle scalaire à bande parabolique : en plus des limites de bande et des masses effectives, des paramètres décrivant le couplage inter-bandes, l’énergie de séparation et (le cas échéant) la contrainte sont intégrés à l’hamiltonien. Vous pouvez examiner les fonctions du GaAs via les paramètres électriques et l’éditeur de code du matériau (??). La capture d’écran montre les potentiels de déformation de la bande de valence b et d ; vous n’avez pas besoin de les modifier pour cette première simulation.

Quelles sont les modifications apportées au calcul basé sur les huit bandes ?

Un modèle parabolique attribue une courbure distincte à chaque bande. L’approche k·p, quant à elle, construit une matrice dont les éléments non diagonaux mettent en relation les états de base sous-jacents. Schématiquement, le problème de l’état confiné est le suivant :

\[ H_{8\times8}(k_x,k_y,-i\partial_z)\,\Psi_n(z) =E_n(k_x,k_y)\,\Psi_n(z). \]

Ici. \(z\) quelle est la direction du confinement ? \(k_x\) et \(k_y\) décrire le mouvement parallèle aux couches. La fonction d’onde de l’enveloppe. \(\Psi_n\) comprend huit composantes, qui décrivent les contributions des différents états de base. Pour ce tutoriel, nous nous concentrerons sur les énergies calculées plutôt que sur l’algèbre matricielle complète. Une explication détaillée est disponible dans le chapitre sur la théorie des huit bandes dans le contexte de la structure cristalline de type zincblende.

4. Lancer la simulation.

Retournez à l’onglet « Fichier » et cliquez sur le bouton bleu « Exécuter la simulation ». Le hamiltonien à huit bandes est assemblé et résolu pour une séquence de vecteurs d’onde dans le plan. Vous pouvez suivre le calcul dans l’onglet « Terminal » (??).

L’exemple de terminal affiche la matrice éparse en cours de construction, suivi des valeurs d’énergie des états propres et de messages tels que « block converged after 12 iterations ». Les états de conduction et de valence sont résolus séparément pour chaque vecteur d’onde successif. Attendez que l’exécution soit terminée et que le bouton revienne à l’état « Lancer la simulation » avant d’examiner les résultats.

5. Interprétation de la structure de bandes calculée.

Le solveur calcule un ensemble discret d’états propres confinés pour chaque vecteur d’onde dans le plan \(k_{\parallel}\). Le graphique affiche donc la dispersion des sous-bandes plutôt que le profil de la limite de bande dans l’espace à travers les couches. Ouvrez l’onglet « Résultats » et consultez le dossier bands, ou ouvrez les fichiers récapitulatifs des bandes affichés dans ??.

5.1 Sous-bandes de conduction

Le diagramme de la bande de conduction (??) présente deux branches qui s’incurvent vers le haut. Au point \(k_{\parallel}=0\), la séparation verticale est due aux différents états confinés. L’augmentation de l’énergie en \(k_{\parallel}\) correspond à un mouvement parallèle au puits, où les porteurs de charge ne sont pas confinés.

À proximité du minimum d’une sous-bande de conduction, la dispersion peut être approximée par :

\[ E_n(k_{\parallel})\simeq E_n(0) +\frac{\hbar^2 k_{\parallel}^{2}}{2m_{n,\parallel}^{*}}. \]

Dans le modèle à plusieurs bandes, il s’agit d’une approximation locale : la courbure, et donc la masse effective apparente, n’a pas besoin de rester constante sur l’ensemble de la plage de vecteurs d’onde représentée. Le couplage entre la bande de conduction et la bande de valence contribue à ce comportement.

5.2 Sous-bandes de valence : mélange et courbure.

Le diagramme de la bande de valence est plus intéressant (??). Ses deux branches représentées partent de niveaux d’énergie différents, s’incurvent vers le bas et semblent se rejoindre à une valeur finie de \(k_{\parallel}\). Dans la bande de valence de la structure cristalline de type zincblende, les états correspondant aux trous lourds et aux trous légers ne se comportent pas comme deux paraboles complètement indépendantes : pour une valeur finie du vecteur d’onde, les termes de l’hamiltonien peuvent modifier leurs propriétés.

Une façon utile d’exprimer cela consiste à considérer deux états voisins, liés par un élément de matrice \(V\) :

\[ H_{\mathrm{illustrative}}= \begin{pmatrix}E_{\mathrm{HH}}(k)&V(k)\\ V^{*}(k)&E_{\mathrm{LH}}(k)\end{pmatrix}. \]

Si on leur permet de s’apparier, leurs énergies se repoussent à l’approche, et leur caractère HH/LH peut varier d’une branche à l’autre. Un croisement réel est également possible lorsque l’appariement pertinent disparaît par symétrie. La capture d’écran actuelle ne permet pas de distinguer ces cas : pour déterminer le caractère de l’état, il faut examiner les vecteurs propres correspondants et, si nécessaire, calculer davantage de points près du point d’intersection apparent. La matrice 2 x 2 ci-dessus n’est qu’une illustration et ne représente pas l’hamiltonien complet à huit bandes.

5.3 Analyse conjointe des états de conduction et de valence.

La représentation combinée (??) permet de placer les états électroniques et les états des électrons de valence sur une échelle d’énergie commune. Dans cet exemple, les états de valence se situent près de la référence d’énergie, tandis que les états de conduction confinés se trouvent autour de 1,5 eV et au-delà. L’écart entre une paire d’états donnée correspond à l’énergie du photon d’une transition optique verticale possible, en tenant compte de l’élément de matrice optique et de l’occupation des porteurs de charge pertinents.

Étant donné que l’impulsion d’un photon est relativement faible, la première approximation utile consiste à considérer qu’une transition optique conserve le vecteur d’onde dans le plan :

\[ \hbar\omega=E_{c,n}(k_{\parallel})-E_{v,m}(k_{\parallel}), \qquad \lambda=\frac{hc}{\hbar\omega}. \]

Différentes paires d’états électroniques et de trous contribuent donc à différentes longueurs d’onde. Les fonctions d’onde à plusieurs bandes déterminent également l’intensité et la polarisation de chaque transition, au-delà de ce qui peut être déduit des énergies seules.

6. Propriétés optiques du puits quantique.

Le calcul à huit bandes fournit les énergies électroniques et les fonctions d’onde d’état mixte pour le modèle optique. L’exemple produit des spectres distincts pour l’émission stimulée, l’absorption, le gain net et l’émission spontanée. Les différentes courbes de chaque graphique correspondent à la séquence des conditions de porteur utilisées par la simulation fournie ; il ne s’agit pas de huit bandes hamiltoniennes distinctes.

6.1 Émission et absorption stimulées

L’émission stimulée se produit lorsqu’un photon incident induit un électron à effectuer une transition radiative d’un état de conduction occupé vers un état de valence inoccupé. L’absorption est le processus inverse : un électron de valence absorbe un photon et passe à un état de conduction disponible. Leurs facteurs d’occupation sont, schématiquement,

\[ \begin{aligned} \text{stimulated emission:}&\quad f_c(1-f_v),\\ \text{absorption:}&\quad f_v(1-f_c). \end{aligned} \]

Les courbes d’émission stimulée (??) et les courbes d’absorption (??) partagent les mêmes énergies de transition sous-jacentes, mais leurs amplitudes diffèrent en fonction des niveaux d’occupation. Dans les résultats présentés, une grande partie de la structure spectrale se situe autour de 790 à 850 nm ; les caractéristiques précises dépendent des sous-bandes configurées, des éléments de matrice, des populations de porteurs et de l’élargissement spectral.

6.2 Gain net : lorsque les émissions dépassent l’absorption.

La courbe du gain net (??) combine les contributions concurrentes de l’émission stimulée et de l’absorption. En soustrayant leurs probabilités d’occupation, on obtient un résultat particulièrement simple :

\[ f_c(1-f_v)-f_v(1-f_c)=f_c-f_v. \]

Ainsi, la transition devient productive lorsqu’il y a plus de porteurs dans l’état de conduction participant que dans l’état correspondant de l’électron de valence. Les valeurs nettes d’amplification négatives représentent des conditions où l’absorption est prédominante, tandis que les valeurs positives correspondent à une émission stimulée qui dépasse l’absorption. Le coefficient d’amplification net présenté ici est le résultat obtenu dans une structure à puits quantiques ; il ne doit pas être interprété automatiquement comme l’amplification modale d’une cavité laser complète, qui dépend également du confinement et d’autres pertes optiques.

6.3 Émission spontanée

L’émission spontanée ne nécessite pas l’intervention d’un photon incident pour déclencher la transition. Elle se produit lorsqu’un électron, se trouvant dans un état énergétique excité, se recombine de manière radiative pour occuper un état énergétique disponible. Le facteur d’occupation correspondant est à nouveau \(f_c(1-f_v)\), mais l’émission spontanée et l’émission stimulée présentent des préfacteurs spectraux et des interprétations physiques différents. Les spectres présentés dans ?? utilisent une échelle verticale logarithmique ; par conséquent, un espacement vertical égal n’implique pas une augmentation absolue égale du taux d’émission des photons.

Pour une explication plus détaillée des intensités de transition, des professions liées au transport et de l’élargissement spectral, veuillez consulter le document suivant : transitions optiques et théorie des jeux pages.

7. Localisation des fichiers de sortie

Une fois la simulation terminée, sélectionnez « Sortie » dans la fenêtre principale (??). Parmi les résultats affichés, vous trouverez Ec.csv, Ev.csv et Ec_Ev.csv pour les graphiques de bandes, ainsi que gain.csv, absorption.csv, net_gain.csv et spontaneous.csv. Les dossiers bands et gain contiennent d’autres résultats du solveur. Double-cliquez sur un fichier CSV pour ouvrir son graphique dans OghmaNano, ou exportez les données numériques pour les analyser ailleurs.

8. Testez-le vous-même : comparez les résultats avec le modèle parabolique.

Dans le cadre d’une deuxième expérience, modifiez l’épaisseur de la couche de GaAs, en la faisant passer de 10 nm à 8 nm, tout en conservant les deux barrières d’AlGaAs à une épaisseur de 40 nm. Exécutez à nouveau la simulation et comparez les énergies propres et la plage de longueurs d’onde des transitions optiques.