Poço Quântico de GaN/AlGaN: Tutorial de 8 Bandas sobre o Modelo k·p da Estrutura Wurtzita.
1. Introdução
Neste tutorial, iremos calcular os estados eletrónicos e a resposta ótica de um poço quântico de GaN entre duas barreiras de AlGaN, utilizando o modelo de oito bandas de wurtzite k·p do OghmaNano. O exemplo apresenta um poço de GaN com 10 nm e duas camadas de AlGaN com 20 nm. Primeiro, executaremos o exemplo fornecido e, em seguida, analisaremos as dispersões das sub-bandas de condução e de valência, bem como os espectros óticos resultantes.
A distinção entre a estrutura cristalina da wurtzita e da esfalerita é importante neste caso. Ambas são estruturas cristalinas semicondutoras, mas a wurtzita apresenta uma preferência pelo eixo c hexagonal. A diferente simetria cristalina altera a divisão da banda de valência, o acoplamento entre os estados eletrónicos e a intensidade da luz polarizada em diferentes direções. Consequentemente, um cálculo de oito bandas baseado na estrutura da esfalerita não pode ser simplesmente reutilizado com os parâmetros do material GaN.
2. Abrir o exemplo de wurtzite com oito bandas.
Na faixa de opções «Ficheiro», clique em «Nova simulação». Na janela que aparece, clique duas vezes em «Exemplos de GaAs» (??). Apesar do nome da categoria, também contém exemplos de poços quânticos de nitreto e silício. Na janela seguinte, clique duas vezes em «GaN Wurtzite 8 band k.p» (??). Isso abre a pilha de materiais preparada e seleciona a demonstração apropriada.
3. Inspeção da camada quântica de GaN/AlGaN
A janela principal apresenta a estrutura de três camadas (??). No cálculo real do poço quântico, o confinamento ocorre ao longo da espessura desta estrutura: as funções de onda eletrónicas variam entre as camadas, enquanto os gráficos de dispersão descrevem o movimento paralelo a elas.
Clique em «Editor de camadas» na barra de ferramentas à esquerda. Como mostrado em ??, o exemplo utiliza a seguinte geometria:
| Camada | Espessura | Tipo |
|---|---|---|
| AlGaN | 20 nanómetros | Ativo / barreira |
| GaN (poço quântico) | 10 nm | Poço quântico |
| AlGaN | 20 nanómetros | Ativo / barreira |
Na sua primeira tentativa, mantenha as configurações dos materiais fornecidos inalteradas. A captura de tela mostra as espessuras das camadas, mas não exibe a fração molar de Al, a configuração de deformação ou os parâmetros de polarização; estes devem ser verificados nas configurações dos materiais do exemplo antes de utilizá-lo para representar uma amostra experimental específica.
4. Por que usar um Hamiltoniano de wurtzite com oito bandas?
Um modelo parabólico simples atribui uma massa efetiva a cada portador de carga e assume que a sua energia varia aproximadamente como o quadrado do seu vetor de onda no plano. Isto é útil perto de um extremo da banda, mas não consegue descrever como os vários estados de valência se misturam entre si à medida que o vetor de onda aumenta. Nos materiais de wurtzita, o ambiente cristalino hexagonal também distingue o movimento e a polarização da luz ao longo do eixo c daqueles no plano perpendicular.
A abordagem de oito bandas mantém uma orbital semelhante à de condução e três orbitais semelhantes às de valência, cada uma com dois componentes de spin. Em termos conceituais, um autoestado tem a seguinte forma:
\[ \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}), \]
onde as funções \(F_{n\alpha}(z)\) descrevem o confinamento ao longo do poço quântico e \(u_\alpha\) representam os oito estados básicos de Bloch. Assim, um determinado estado confinado pode ter contribuições de mais de uma orbital da banda de valência; é por isso que descrever cada curva colorida usando um único rótulo fixo para buracos pesados ou buracos leves pode ser enganoso, especialmente quando nos afastamos da borda da banda.
A organização da banda de valência da wurtzita é afetada pelas interações do campo cristalino e do momento angular. Em uma notação concisa, o cálculo resolve:
\[ 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). \]
Note que o vetor de onda no plano é mantido explicitamente, enquanto o vetor de onda na direção do crescimento é substituído por um operador diferencial. Ao repetir este problema de autovalores para diferentes valores de \(k_{\parallel}\), obtêm-se os gráficos de dispersão das sub-bandas apresentados abaixo. Uma discussão mais aprofundada da base, dos termos do campo cristalino e dos coeficientes dos materiais é apresentada no capítulo teoria de Wurtzite k·p.
5. Execute a simulação.
Volte para a faixa de opções «Ficheiro» e clique em «Executar simulação». Abra o separador «Terminal» para acompanhar o progresso. Em ??, o terminal identifica o solucionador WZ KP8, apresenta informações sobre a montagem da matriz esparsa, imprime dados sobre a convergência e avança através dos vetores de onda amostrados. Resolve os estados eletrónicos e de valência separadamente em cada ponto.
Uma vez que este é um problema de função de envoltória com múltiplos componentes, a sua resolução é mais complexa do que a resolução de um potencial parabólico escalar independente para cada banda. Aguarde até que o solucionador de autovalores e os cálculos óticos subsequentes estejam concluídos antes de verificar os ficheiros resultantes.
5.1. Inspeção das intensidades das transições óticas
Mais tarde, durante a execução do programa, o terminal apresenta os valores da intensidade dos elementos da matriz ótica (??). Para uma estrutura convencional de wurtzita no plano c, TE refere-se ao campo elétrico ótico perpendicular ao eixo c, enquanto TM refere-se à sua componente paralela a esse eixo. Estas polarizações não precisam de se acoplar igualmente a uma determinada transição eletrónica.
Uma forma útil de analisar a intensidade óptica é através do elemento da matriz de transição ao quadrado:
\[ 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. \]
Aqui, \(\mathbf{e}\) especifica a polarização da luz e \(\mathbf{p}\) é o operador do momento. Um estado pode estar presente na estrutura de bandas eletrónicas, mas contribuir muito pouco para uma polarização ótica específica se o seu elemento de matriz for pequeno. Os valores relativos de TE/TM no terminal fornecem, portanto, informações físicas adicionais para além das posições dos níveis de energia.
6. Compreensão das bandas eletrónicas calculadas
O eixo horizontal de cada gráfico de bandas representa o vetor de onda no plano \(k_{\parallel}\), e não a posição dentro do poço quântico. O eixo vertical representa a energia de uma sub-banda confinada. Em \(k_{\parallel}=0\), cada curva começa em um dos níveis de energia de confinamento discretos. Ao mover-se para a direita, observa-se como o mesmo estado eletrónico se dispersa quando a carga se move paralelamente às camadas.
6.1 Estados da banda de condução
Abra a saída da banda de condução e examine a dispersão em ??. Vários níveis confinados são visíveis, começando em diferentes energias próximas à borda da banda do GaN. Para um vetor de onda pequeno, uma aproximação útil para uma sub-banda de condução é:
\[ E_{c,n}(k_{\parallel})\simeq E_{c,n}(0) +\frac{\hbar^2 k_{\parallel}^2}{2m_{e,\parallel}^{*}}. \]
O termo \(E_{c,n}(0)\) contém a energia de confinamento, enquanto \(m_{e,\parallel}^{*}\) descreve a resposta do portador ao movimento no plano. Esta relação ajuda a explicar as curvas que apresentam um aumento gradual, embora o cálculo completo com oito bandas mantenha o acoplamento interbanda em vez de impor uma forma parabólica fixa.
6.2 Estados da banda de valência
O diagrama da banda de valência (??) é mais complexo. O seu ramo superior sobe inicialmente, afastando-se de \(k_{\parallel}=0\), antes de se curvar para baixo; vários ramos inferiores apresentam curvaturas visivelmente diferentes. Formas como estas são possíveis num modelo de banda de valência de wurtzite acoplado, no qual a divisão do campo cristalino, o acoplamento spin-órbita, o confinamento e a mistura de bandas atuam em conjunto.
Este gráfico não deve ser interpretado como prova de que cada curva colorida possui uma identidade fixa nas bandas A, B ou C ao longo de todo o intervalo do vetor de onda. Para realizar uma atribuição definitiva, analise as contribuições dos componentes do vetor próprio subjacente. Da mesma forma, é útil comparar os extremos incomuns e descentralizados com o material e os parâmetros de deformação escolhidos, caso esteja a utilizar o modelo de forma quantitativa.
7. Transições ópticas, emissão estimulada e absorção.
As bandas eletrónicas determinam as energias das possíveis transições óticas, mas as intensidades das transições também dependem das funções de onda dos eletrões e das lacunas, da polarização e da ocupação dos portadores de carga. Para um vetor de onda comum no plano, a energia do fotão para uma transição vertical é:
\[ \hbar\omega=E_{c,n}(\mathbf{k}_{\parallel}) -E_{v,m}(\mathbf{k}_{\parallel}). \]
A regra de ouro de Fermi explica por que os elementos da matriz ótica impressos no terminal são importantes. De forma esquemática, a taxa para uma transição é:
\[ 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). \]
O elemento da matriz de perturbação \(H'_{cv}\) depende da polarização da luz e dos estados envolvidos. A função delta impõe a conservação da energia; ao construir um espectro prático, as contribuições das transições permitidas são somadas sobre os estados amostrados e alargadas. Os fatores de ocupação das portadoras determinam então se a absorção ou a emissão estimulada predominam em um determinado comprimento de onda. Um fator líquido esquemático conveniente é \(f_c-f_v\), onde \(f_c\) e \(f_v\) são, respetivamente, as ocupações eletrónicas dos estados de condução e de valência.
Uma vez que a camada de GaN apresenta uma grande diferença de energia eletrónica, os seus espectros calculados surgem na região do ultravioleta próximo, predominantemente entre 330 e 370 nm, como se pode observar no exemplo fornecido. O comprimento de onda e a energia do fotão estão relacionados pela equação \(E_{\gamma}=hc/\lambda\).
7.1 Emissão estimulada e ganho líquido
O gráfico da emissão estimulada (??) contém uma família de espectros para diferentes condições de portadora. O cálculo do ganho líquido (??) considera a competição entre a emissão e a absorção. Um ganho material positivo indica amplificação dentro deste modelo eletrónico, mas, por si só, não demonstra que uma cavidade laser completa atingiu o limiar: o fator de confinamento e as perdas do ressonador também devem ser considerados.
7.2 Absorção
Os resultados da absorção são apresentados em ??. As suas bordas e características espectrais são determinadas pelas separações entre os níveis de valência ocupados e os níveis de condução disponíveis, ponderadas pelos seus elementos de matriz ótica. Numa estrutura do tipo wurtzite, os respetivos espectros TE e TM podem diferir, uma vez que o conteúdo orbital de um estado determina a intensidade com que este se acopla a cada polarização.
As características presentes nestes espectros são o resultado do modelo configurado e do alargamento espectral. Não devem ser automaticamente identificadas com picos excitónicos específicos: um cálculo de transição interbanda do tipo k·p para uma única partícula não constitui, por si só, uma descrição completa dos excitões ligados Coulombianamente entre eletrões e lacunas.
8. Emissão espontânea
Ao contrário da emissão estimulada, a emissão espontânea ocorre sem que um fotão incidente inicie a transição. Depende da disponibilidade de um eletrão no estado de condução e de um estado vazio na banda de valência (uma lacuna). O seu fator de ocupação eletrónica é, portanto, proporcional a
\[ 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), \]
onde \(L\) representa o alargamento espectral escolhido e a relação de proporcionalidade exclui os pré-fatores fotónicos. O resultado é apresentado em ??. A escala logarítmica permite comparar espectros que abrangem várias ordens de magnitude. Observe as características óticas perto da parte do espectro com comprimentos de onda mais curtos e a mudança na magnitude geral nas diferentes condições apresentadas.
9. Experimente você mesmo.
💡 Dois cálculos úteis para complementar a análise:
- Altere a largura da camada de GaN. Tente reduzir a camada central de 10 nm para 7 nm, mantendo as barreiras inalteradas. Recalcule os autoestados e compare os níveis de condução mais baixos e as energias de transição ótica. Um confinamento mais forte geralmente aumentará a separação entre os níveis confinados, embora a resposta detalhada da banda de valência seja multibanda.
- Compare as intensidades de TE e TM. Utilize as intensidades óticas reportadas pelo terminal para identificar as transições que favorecem fortemente uma polarização específica. Analise como isso se relaciona com o caráter do estado de valência, em vez de presumir que todas as transições contribuem igualmente para o espectro.
✅ Em que deve prestar atenção?
Um poço mais estreito geralmente aumenta as energias de confinamento dos eletrões e altera os níveis de lacunas, modificando as energias de transição interbanda disponíveis. As relações TE/TM podem variar significativamente de um par de estados para outro, pois investigam diferentes componentes orbitais das funções de onda da estrutura cristalina wurtzita. Para um estudo quantitativo do material, também é importante analisar os efeitos da composição da liga, da tensão e dos campos de polarização.