GaInNAs: Tutorial sobre poços quânticos de nitreto diluído e o modelo k·p de 10 bandas.
1. Introdução
Neste tutorial, simularemos uma estrutura de poço quântico GaInNAs com baixa concentração de nitrogênio, utilizando o modelo de dez bandas k·p do OghmaNano para estruturas de tipo zinco-blenda. Abriremos um exemplo existente, examinaremos a estrutura das camadas, executaremos o cálculo do estado eletrônico e analisaremos os espectros óticos previstos. O exemplo contém um poço quântico GaInNAs de 7 nm dentro de uma estrutura maior de GaAs/AlGaAs, em vez de um poço de três camadas isolado.
Qual a utilidade de um modelo de dez bandas? A teoria k·p de oito bandas convencional descreve uma banda de condução e as bandas de valência de lacunas pesadas, lacunas leves e lacunas separadas, incluindo o spin. A introdução de uma pequena concentração de nitrogénio numa liga III–V pode criar um nível de ressonância relacionado com o nitrogénio, que interage particularmente com a banda de condução do material hospedeiro. Uma descrição de dez bandas adiciona dois componentes de spin a este estado relacionado com o nitrogénio. A interação pode reduzir significativamente a energia de transição fundamental e alterar a curvatura da banda de condução. Isto torna os nitretos diluídos interessantes para dispositivos óticos de comprimentos de onda mais longos.
2. Abrir o exemplo do equalizador de dez bandas.
Selecione «Ficheiro» → «Nova simulação». Na janela «Nova simulação», clique duas vezes em «Demonstrações GaAs» (??). A janela seguinte contém vários exemplos de poços quânticos. Selecione «AlGaAs, nitreto diluído, 10 bandas k.p» (??). Isso carrega a simulação pré-configurada de dez bandas, incluindo as definições de camada e os ficheiros de parâmetros de materiais adequados.
3. Estrutura do dispositivo e parâmetros dos materiais.
A janela principal exibe a estrutura multicamadas (??). Ao contrário das estruturas simples de poços GaAs/AlGaAs utilizadas nos tutoriais anteriores, este exemplo inclui camadas adicionais de revestimento, camadas de transição e camadas de barreira em torno do poço de nitreto diluído. Para verificar a geometria, clique em «Editor de camadas» na barra de ferramentas à esquerda.
O editor de camadas (??) apresenta oito entradas, que vão desde o p-GaAs, passando pelo AlGaAs e pela região ativa, até ao n-GaAs. Em particular, a camada central de GaInNAs (poço quântico) tem 7,0 nm de espessura e o seu tipo de camada está definido como «poço quântico». O exemplo visível também inclui camadas de AlGaAs de 0,16 µm de cada lado da região ativa e uma camada de barreira/SCL de GaAs de 0,02 µm. Não é necessário reconstruir a estrutura manualmente ao utilizar o modelo fornecido.
3.1 Verificação das definições dos materiais
Clique em «Parâmetros elétricos», selecione o separador «GaInNAs (QW)» e abra o editor de código (??). A captura de ecrã ilustra como o ficheiro de material contém grandezas físicas, como a densidade de estados efetiva da banda de valência e a mobilidade elétrica. Por exemplo, a mobilidade dos eletrões em campo baixo, visível, é 0.05 m² V⁻¹ s⁻¹, o que equivale a 500 cm² V⁻¹ s⁻¹. Estes são parâmetros de entrada do modelo do material, e não valores inferidos da dispersão do poço quântico representada no gráfico.
Numa simulação com baixo teor de nitreto, é especialmente importante que a composição, as energias da banda de condução e os parâmetros de acoplamento do nitrogénio sejam mutuamente consistentes. Para esta primeira execução, deixe os valores padrão inalterados; discutiremos um experimento útil sobre a composição no final da página.
4. Que alterações ocorrem quando introduzimos nitrogénio?
Antes de executar o modelo, é importante compreender por que razão o Hamiltoniano apresenta dez bandas em vez de oito. Imagine que a extremidade da banda de condução do semicondutor hospedeiro se desloca para cima no espectro de energia à medida que o vetor de onda do eletrão no plano aumenta. Um estado de ressonância relacionado com o nitrogénio tem a sua própria energia característica. Se os dois estados não interagissem, as suas energias poderiam aproximar-se e cruzar-se. No entanto, num nitreto diluído, o nitrogénio acopla-se aos estados de condução do material hospedeiro, alterando os níveis de energia resultantes.
Um exemplo prático que ilustra este conceito da física é a matriz de cruzamento anti-inverso de duas bandas (BAC):
\[ H_{\mathrm{BAC}}(k)= \begin{pmatrix} E_C(k) & V_N\\ V_N & E_N \end{pmatrix}, \qquad V_N=\beta_N\sqrt{x_N}. \]
Aqui, \(E_C(k)\) representa a dispersão da banda de condução do material hospedeiro, \(E_N\) é a energia do nível relacionado com o nitrogénio, \(x_N\) é a fração de nitrogénio e \(\beta_N\) controla a intensidade do acoplamento. Ao diagonalizar esta pequena matriz, obtêm-se duas ramificações:
\[ E_{\pm}(k)=\frac{E_C(k)+E_N}{2} \pm\frac{1}{2}\sqrt{\big[E_C(k)-E_N\big]^2+4V_N^2}. \]
O termo da raiz quadrada impede que os ramos acoplados se cruzem quando \(V_N\neq0\). Em particular, o ramo inferior, com comportamento semelhante ao de uma banda de condução, pode mover-se para baixo em relação à banda de condução original, reduzindo a banda proibida fundamental. Uma vez que a interação depende da separação entre \(E_C(k)\) e \(E_N\), a dispersão também pode tornar-se fortemente não parabólica.
Este exemplo, que considera um sistema bidimensional, explica a ideia central; não se trata da solução completa do Hamiltoniano do poço quântico apresentada neste tutorial. O modelo de dez bandas k·p mantém a base mais ampla de condução/valência do tipo esfalerita, incluindo o spin, e adiciona os componentes relacionados ao nitrogénio. Em seguida, é incluído o confinamento quântico ao longo da direção de crescimento. Para o tratamento teórico completo, consulte o capítulo sobre o modelo de dez bandas do tipo esfalerita.
5. Executar a simulação.
Volte para a faixa de opções «Ficheiro» e clique em «Executar simulação». Selecione o separador «Terminal» para acompanhar o cálculo (??). O solucionador cria uma matriz hamiltoniana esparsa, obtém os autoestados de condução e valência confinados e repete o cálculo numa sequência de vetores de onda no plano \(k_{\parallel}\).
Na captura de tela, o registo identifica o solucionador KP10, apresenta os níveis de energia calculados e indica a convergência do bloco após um determinado número de iterações. Os pares de níveis de energia próximos podem incluir degenerescências relacionadas ao spin; não devem ser interpretados automaticamente como ramos distintos de anticruzamento do nitrogênio. Aguarde o término da execução antes de analisar os arquivos de saída.
6. Compreensão da estrutura de bandas calculada.
Os gráficos resultantes mostram a energia em função do vetor de onda no plano, e não a energia em função da profundidade no dispositivo. O potencial do poço quântico confina as cargas ao longo da direção de crescimento da camada, produzindo estados discretos; o movimento paralelo às camadas confere a cada estado confinado uma dispersão \(E_n(k_{\parallel})\). Num modelo parabólico simples, essa dispersão é quadrática. Um Hamiltoniano acoplado de dez bandas permite que ele se desvie de uma parábola.
6.1 Estados de condução e a redução da folga.
O gráfico de condução (??) mostra dois ramos calculados que partem de diferentes níveis de energia. O nível de condução mais baixo começa aproximadamente em 0,89 eV na referência de energia escolhida para a simulação. Esse valor não deve ser confundido com a banda proibida do GaInNAs: é uma autoenergia desta estrutura confinada específica. Note que as duas curvas representadas não têm a mesma curvatura. Este é um dos aspetos a investigar ao comparar um cálculo de oito bandas com o modelo de nitreto diluído.
As equações simplificadas do BAC apresentadas acima ajudam a explicar por que a dispersão da condução pode variar em um material que contenha nitrogênio. No entanto, por si só, elas não comprovam a proporção de nitrogênio em nenhuma das curvas coloridas. Para isso, é necessário analisar os autovetores correspondentes ou variar \(x_N\) e observar como os estados evoluem.
6.2 Estados de valência
O diagrama de valência (??) contém vários ramos com formas claramente não parabólicas. As características de «buraco pesado» e «buraco leve» podem combinar-se para um vetor de onda no plano finito; as energias de confinamento, a interação spin-órbita e o acoplamento entre os estados de base afetam todos as curvas. Uma mudança ascendente ou uma alteração na curvatura não indicam necessariamente uma falha numérica. Por outro lado, uma curva traçada de forma suave não é, por si só, prova de que os parâmetros do material estão corretos.
É preferível considerar as cores individuais como ramos de autoestados representados num gráfico, em vez de atribuir permanentemente as etiquetas HH ou LH a todo o intervalo do vetor de onda, sem verificar as suas características.
6.3 Analisando as duas bandas em conjunto.
Ao ativar a dispersão combinada (??), os estados eletrónicos e de valência são colocados na mesma escala de energia. A transição interbanda fundamental está relacionada com a sua separação, e não com a energia absoluta de qualquer um dos estados representados:
\[ \hbar\omega_{nm}(k_{\parallel}) =E_{c,n}(k_{\parallel})-E_{v,m}(k_{\parallel}). \]
Esta expressão descreve uma transição ótica direta, aproximadamente vertical, no vetor de onda no plano: o momento do fotão é pequeno em relação à escala do momento do cristal eletrónico. O espectro observado não é, no entanto, geralmente uma única linha nítida. Muitas transições permitidas contribuem, com intensidades e níveis de ocupação que variam de estado para estado.
7. Dos estados eletrónicos aos espectros óticos
Com os estados calculados disponíveis, o OghmaNano pode construir espectros óticos. A questão física importante não é apenas se um eletrão e uma lacuna têm a diferença de energia correta, mas também a intensidade com que os dois estados se acoplam ao campo eletromagnético e se os seus estados inicial e final estão devidamente ocupados.
7.1 Por que razão algumas transições têm um impacto maior do que outras?
A regra de ouro de Fermi oferece um ponto de partida útil para compreender as taxas de transição. Para uma transição de um estado eletrónico inicial \(i\) para um estado final \(f\), induzida por uma perturbação ótica fraca, a sua forma esquemática é:
\[ W_{i\rightarrow f}=\frac{2\pi}{\hbar} \left|\langle f|H_{\mathrm{int}}|i\rangle\right|^2 \delta(E_f-E_i-\hbar\omega). \]
O elemento da matriz \(\langle f|H_{\mathrm{int}}|i\rangle\) descreve a interação dos estados com a luz. Depende, entre outros fatores, do seu caráter orbital, das funções de envelope e da polarização da luz. A função delta garante a conservação da energia. Num espetro numérico, a transição idealizada e abrupta é representada utilizando uma função apropriada de alargamento da linha. Após a soma sobre as sub-bandas permitidas e os vetores de onda no plano, o resultado é um espetro, e não uma única energia de transição.
Nesse contexto, as ocupações dos níveis de energia distinguem a absorção da emissão. Se \(f_c\) representa a ocupação de um estado de condução e \(f_v\) a ocupação de um estado de elétrons de valência, os fatores característicos são:
\[ \begin{aligned} \text{absorption:}\quad &(1-f_c)f_v,\\ \text{stimulated emission:}\quad &f_c(1-f_v),\\ \text{net stimulated contribution:}\quad &f_c-f_v. \end{aligned} \]
A última expressão representa a diferença entre os processos óticos ascendentes e descendentes. Uma contribuição positiva para o ganho requer uma inversão de população suficiente. Estes fatores explicam a interpretação dos gráficos que se seguem; para mais detalhes, consulte a secção teoria da transição ótica.
7.2 Emissão e absorção estimuladas
Os dois gráficos (?? e ??) ilustram processos óticos concorrentes para a mesma estrutura de bandas de poço quântico. As várias curvas representam a sequência de condições utilizadas no exemplo; as suas magnitudes relativas variam porque a ocupação de portadores determina quais as transições possíveis.
Os espectros estendem-se até o infravermelho próximo, abrangendo, de forma geral, o intervalo de aproximadamente 1,1 a 1,5 µm. Para referência, a energia dos fotões e o comprimento de onda no vácuo estão relacionados pela seguinte equação:
\[ E_{\mathrm{photon}}=\frac{hc}{\lambda} \simeq \frac{1240\ \mathrm{eV\,nm}}{\lambda\,[\mathrm{nm}]}. \]
Dessa forma, 1300 nm corresponde a aproximadamente 0,954 eV. Isso permite uma verificação útil da consistência entre a separação de energia no diagrama de bandas combinado e a escala dos espectros óticos, tendo em conta o confinamento, o intervalo das transições que contribuem e a referência de energia escolhida.
7.3 Ganho líquido: quando a emissão supera a absorção.
A quantidade relevante para um meio amplificador é o equilíbrio entre a emissão estimulada e a absorção. Utilizando a convenção de sinais empregada no gráfico, pode ser representada esquematicamente da seguinte forma:
\[ g_{\mathrm{net}}(\lambda) =g_{\mathrm{stim}}(\lambda)-\alpha(\lambda). \]
O gráfico do ganho líquido calculado (??) inclui tanto curvas negativas quanto positivas. Quando a absorção é predominante, \(g_{\mathrm{net}}\) é negativo. À medida que as populações de eletrões e lacunas se modificam, a emissão estimulada pode tornar-se suficientemente forte para que certas partes do espectro apresentem ganho positivo. A região de transição representada no gráfico apresenta uma estrutura notável, em vez de ser um único pico perfeitamente suave, porque múltiplas transições em estados confinados contribuem para o resultado.
Importante: o simples facto de se obter um ganho material positivo não é, por si só, um indicador de que um laser completo funcionará acima do limiar. As perdas na cavidade, o confinamento ótico, a refletividade dos espelhos e a injeção de portadores também devem ser considerados.
8. Localizar os resultados da simulação
bands e gain.Assim que o cálculo for concluído, abra a guia «Resultados» (??). A captura de ecrã mostra os principais ficheiros de resultados, incluindo: Ec.csv (ramos de condução), Ev.csv (ramos de valência), Ec_Ev.csv (bandas combinadas), gain.csv, absorption.csv, net_gain.csv e spontaneous.csv. Os resultados adicionais estão organizados nas pastas bands e gain. Clique duas vezes num ficheiro CSV para gerar um gráfico ou utilize os dados numéricos guardados para análises adicionais.
Embora este tutorial ilustre a emissão estimulada, a absorção e o ganho líquido, o exemplo também gera uma saída de emissão espontânea. A emissão espontânea é uma recombinação radiativa que não requer a presença de um fotão incidente, e o seu fator de ocupação característico é \(f_c(1-f_v)\). Os seus fatores espectrais absolutos diferem dos da emissão estimulada.
9. Experimente você mesmo: investigue a interação do nitrogênio.
Experimento sugerido: Guarde uma cópia separada deste exemplo antes de o editar. Em seguida, analise a concentração de nitrogénio na definição do material GaInNAs, alterando um parâmetro de cada vez, mantendo a geometria do poço de 7 nm inalterada.
- Registe o valor mais baixo da energia própria dos eletrões em \(k_{\parallel}=0\).
- Altere a fração de nitrogênio, de modo que se mantenha dentro de uma faixa fisicamente adequada para a parametrização selecionada.
- Recalcule a dispersão da condução e compare a forma do ramo inferior.
- Compare a posição da borda de absorção e as principais características de ganho.
A ilustração do modelo BAC prevê que o aumento do acoplamento de nitrogénio pode aumentar a separação entre os ramos acoplados e alterar a energia de condução mais baixa. No entanto, num cálculo completo da liga, a composição pode também afetar outros parâmetros do material. Portanto, interprete quaisquer alterações resultantes utilizando as funções reais do material ativadas na sua simulação, em vez de assumir que a simples fórmula de dois estados representa o modelo completo da liga.
✅ O que procurar
Investigue se a extremidade da banda de condução se desloca e se a sua curvatura se altera à medida que a concentração de nitrogénio varia. Compare os espectros óticos em configurações de ocupação idênticas, de modo a que uma alteração nas energias de transição não seja confundida com uma alteração não relacionada na população de portadores. Se for possível analisar a composição dos autovetores, é particularmente útil observar as alterações nas contribuições da banda de condução do material hospedeiro e as contribuições relacionadas com o nitrogénio para os estados.
Para um segundo experimento, mantenha os parâmetros do material fixos e altere a largura da barreira quântica de 7 nm para 6 nm. O aumento do confinamento quântico geralmente eleva a energia do elétron confinado de menor nível e altera as energias de transição interbanda permitidas. A magnitude desse deslocamento é determinada pelo cálculo acoplado de dez bandas, e não simplesmente pela atribuição de uma única massa efetiva constante a cada estado.