Poço Quântico de GaAs/AlGaAs: Tutorial de 8 Bandas sobre a Estrutura Cristalina Tipo Zinco-Blenda.

1. Introdução

Neste tutorial, iremos calcular as propriedades eletrónicas e óticas de um poço quântico de GaAs/AlGaAs, utilizando o modelo de oito bandas de tipo zinco-blenda k·p do OghmaNano. Começaremos com um exemplo pré-definido, analisaremos a estrutura das camadas, executaremos a simulação e examinaremos as sub-bandas de condução e de valência resultantes. Por fim, analisaremos os espectros óticos calculados a partir dos estados.

O/A tutorial sobre poços quânticos parabólicos o modelo trata os eletrões e as lacunas utilizando bandas de massa efetiva relativamente simples. Aqui, o Hamiltoniano inclui uma base da banda de condução e as bases da banda de valência correspondentes às lacunas pesadas (HH), lacunas leves (LH) e lacunas separadas (SO), com dois componentes de spin para cada uma: um total de oito estados base. O acoplamento entre estes estados permite que o modelo descreva efeitos que parábolas independentes não conseguem, nomeadamente a mistura da banda de valência e as alterações na curvatura em função do vetor de onda no plano.

2. Abrir o exemplo com oito bandas.

No OghmaNano, clique em «Nova simulação» na faixa de opções «Ficheiro». Clique duas vezes na pasta «Exemplos de GaAs» (??). Nos exemplos de poços quânticos, clique duas vezes em «GaAs/AlGaAs QW 8 band k.p» (??). Isto abre uma simulação pré-configurada; não é necessário introduzir manualmente uma matriz de Hamiltoniana de 8x8.

3. Estrutura de poço quântico e parâmetros dos materiais

A janela principal exibe uma estrutura de três camadas de AlGaAs/GaAs/AlGaAs (??). Abra o editor de camadas usando o botão no lado esquerdo. O exemplo fornecido consiste em duas barreiras de AlGaAs de 40 nm que envolvem uma camada de GaAs de 10 nm:

As extremidades da banda de energia do GaAs situam-se dentro da banda proibida do AlGaAs circundante, proporcionando confinamento ao longo da direção de crescimento da camada. No editor de camadas, a camada central de GaAs é identificada especificamente como uma barreira quântica (??). As propriedades dos materiais óticos e os tipos de camada já estão configurados neste exemplo.

O solucionador de oito bandas requer mais informações sobre o material do que um modelo escalar parabólico de banda: além das extremidades da banda e das massas efetivas, parâmetros que descrevem o acoplamento interbanda, a energia de separação e (quando aplicável) a tensão entram no Hamiltoniano. Pode consultar as funções do GaAs através dos parâmetros elétricos e do editor de código do material (??). A captura de tela mostra os potenciais de deformação da banda de valência b e d; não é necessário modificá-los para esta execução inicial.

Que alterações são introduzidas no cálculo das oito bandas?

Um modelo parabólico atribui uma curvatura distinta a cada banda. A abordagem k·p, por outro lado, constrói uma matriz cujos elementos fora da diagonal principal relacionam os estados básicos subjacentes. De forma esquemática, o problema do estado confinado é:

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

Aqui, \(z\) representa a direção de confinamento; \(k_x\) e \(k_y\) descrevem o movimento paralelo às camadas. A função de onda do envelope \(\Psi_n\) possui oito componentes, que descrevem as contribuições dos diferentes estados de base. Para este tutorial, concentraremos-nos nas energias calculadas, em vez de na álgebra matricial completa. Uma explicação detalhada está disponível no capítulo teoria das oito bandas de zinco-blenda.

4. Executar a simulação.

Volte para a faixa de opções «Ficheiro» e clique no botão azul «Executar simulação». O Hamiltoniano de oito bandas é montado e resolvido numa sequência de vetores de onda no plano. Pode acompanhar o cálculo na aba «Terminal» (??).

O exemplo apresentado mostra a construção da matriz esparsa, seguida pelas energias dos autoestados e mensagens como block converged after 12 iterations. Os estados de condução e de valência são resolvidos separadamente para cada vetor de onda sucessivo. Aguarde até que a execução seja concluída e o botão volte ao estado «Executar simulação» antes de analisar os resultados.

5. Compreensão da estrutura de bandas calculada.

O solucionador calcula um conjunto discreto de autoestados confinados para cada vetor de onda no plano \(k_{\parallel}\). O gráfico mostra, portanto, a dispersão das sub-bandas, em vez do perfil espacial da extremidade da banda ao longo das camadas. Abra a guia «Resultados» e procure na pasta bands ou abra os ficheiros de resumo das bandas exibidos em ??.

5.1 Sub-bandas de condução

O diagrama da banda de condução (??) apresenta dois ramos com curvatura ascendente. Em \(k_{\parallel}=0\), a separação vertical tem origem nos diferentes estados confinados. O aumento da energia em \(k_{\parallel}\) reflete o movimento paralelo ao poço, onde as cargas não estão confinadas.

Próximo do mínimo da sub-banda de condução, a dispersão pode ser aproximada por:

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

No modelo multibanda, isto representa uma aproximação local: a curvatura e, consequentemente, a massa efetiva aparente, não precisam de permanecer constantes em toda a faixa do vetor de onda representada no gráfico. O acoplamento entre as bandas de condução e de valência contribui para este comportamento.

5.2 Subbandas de valência: mistura e curvatura

O diagrama da banda de valência é mais interessante (??). As duas ramificações representadas partem de energias diferentes, curvam-se para baixo e parecem convergir num valor finito de \(k_{\parallel}\). Na banda de valência do tipo esfalerita, os estados de lacuna pesada e de lacuna leve não se comportam como duas parábolas completamente independentes: num vetor de onda finito, os termos no Hamiltoniano podem misturar as suas características.

Uma forma útil de expressar isto é considerar dois estados próximos, que estão interligados por um elemento de matriz \(V\):

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

Se for permitido que se combinem, as suas energias repelem-se quando estão prestes a interagir, e o seu caráter HH/LH pode ser trocado entre os ramos. Um cruzamento genuíno também é possível, no qual a combinação relevante desaparece devido à simetria. A captura de ecrã apresentada não é suficiente para distinguir estes casos: para determinar o caráter do estado, examine os vetores próprios correspondentes e, se necessário, calcule mais pontos perto da interseção aparente. A matriz 2x2 acima é apenas uma ilustração, não o Hamiltoniano completo de oito bandas.

5.3 Análise conjunta dos estados de condução e de valência.

O gráfico combinado (??) representa os estados eletrónicos e de eletrões de valência numa escala de energia comum. Neste exemplo, os estados de valência situam-se perto da referência de energia, enquanto os estados de condução confinados estão situados em torno de 1,5 eV e acima. A diferença entre um determinado par de estados corresponde à energia do fotão de uma possível transição ótica vertical, sujeita ao elemento de matriz ótico relevante e à ocupação dos portadores de carga.

Uma vez que o momento do fotão é relativamente pequeno, a primeira aproximação útil é que uma transição ótica preserva o vetor de onda no plano:

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

Diferentes pares de estados de eletrões e lacunas contribuem, portanto, em diferentes comprimentos de onda. As funções de onda multibanda também determinam a intensidade e a polarização de cada transição, indo além do que pode ser inferido apenas a partir das energias.

6. Propriedades óticas do poço quântico

O cálculo de oito bandas fornece energias eletrónicas e funções de onda de estado misto para o modelo ótico. O exemplo gera espectros distintos para emissão estimulada, absorção, ganho líquido e emissão espontânea. As múltiplas curvas em cada gráfico correspondem à sequência de condições de portador utilizadas na simulação fornecida; não se trata de oito bandas hamiltonianas individuais.

6.1 Emissão e absorção estimuladas

A emissão estimulada ocorre quando um fotão incidente induz um eletrão a realizar uma transição radiativa de um estado de condução preenchido para um estado de valência vazio. A absorção é o processo inverso: um eletrão de valência absorve um fotão e move-se para um estado de condução disponível. Os seus fatores de ocupação são, esquematicamente,

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

As curvas de emissão estimulada (??) e as curvas de absorção (??) partilham as energias de transição subjacentes, mas as suas intensidades diferem à medida que as populações dos níveis de energia se alteram. Nestes resultados, grande parte da estrutura espectral situa-se entre 790 e 850 nm; as características exatas dependem das subbandas configuradas, dos elementos de matriz, das populações de portadores e do alargamento das linhas espectrais.

6.2 Ganho líquido: quando a emissão é superior à absorção.

O gráfico do ganho líquido (??) combina as contribuições da emissão estimulada e da absorção, que competem entre si. Subtraindo as suas probabilidades de ocupação, obtém-se um resultado particularmente simples:

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

Dessa forma, a transição passa a gerar ganho quando o estado de condução em que participa está mais ocupado do que o estado correspondente do eletrão de valência. Valores negativos de ganho líquido representam condições em que a absorção é predominante; valores positivos correspondem à emissão estimulada que excede a absorção. O coeficiente de ganho líquido apresentado aqui é o resultado ótico obtido em uma estrutura de poço quântico; não deve ser interpretado automaticamente como o ganho modal de uma cavidade laser completa, que também depende do confinamento e de outras perdas óticas.

6.3 Emissão espontânea

A emissão espontânea não requer um fotão incidente para desencadear a transição. Ocorre quando um estado eletrónico ocupado se recombina de forma radiativa num estado de valência disponível. O fator de ocupação correspondente é, novamente, \(f_c(1-f_v)\), mas a emissão espontânea e a emissão estimulada têm diferentes fatores espectrais e interpretações físicas. Os espectros em ?? utilizam uma escala vertical logarítmica; portanto, um espaçamento vertical igual não significa aumentos absolutos iguais na taxa de fotões emitidos.

Para uma explicação mais completa sobre as características das transições, as profissões relacionadas e o alargamento espectral, consulte as páginas transições óticas e teoria do ganho.

7. Localizar os ficheiros de saída

Após a conclusão da simulação, selecione a opção «Resultados» na janela principal (??). Entre os resultados visíveis, encontram-se Ec.csv, Ev.csv e Ec_Ev.csv para os gráficos de banda, juntamente com gain.csv, absorption.csv, net_gain.csv e spontaneous.csv. As pastas bands e gain contêm resultados adicionais do solucionador. Clique duas vezes num ficheiro CSV para abrir o respetivo gráfico no OghmaNano ou exporte os dados numéricos para análise noutro local.

8. Experimente você mesmo: compare os resultados com o modelo parabólico.

Como segundo experimento, altere a espessura da camada de GaAs de 10 nm para 8 nm, mantendo as duas barreiras de AlGaAs de 40 nm inalteradas. Execute novamente o modelo e compare as energias dos níveis quantizados e o intervalo de comprimento de onda das transições ópticas.