Simulação de Diodo de Junção PN de Nitreto de Gálio (GaN): Difusão-Deriva, Curvas J–V e Recombinação.
1. Introdução
Um díodo de junção PN de nitreto de gálio (GaN) consiste em regiões de tipo p e tipo n do mesmo semicondutor, unidas para formar uma junção retificadora. Os eletrões e as lacunas redistribuem-se na interface, produzindo uma região de depleção e uma barreira de potencial intrínseca. Sob polarização direta, a barreira é reduzida e a corrente aumenta acentuadamente; sob polarização inversa, a junção suprime a injeção de portadores. A resposta assimétrica resultante entre corrente e tensão é o princípio básico de funcionamento de um retificador semicondutor.
O GaN é um semicondutor III-V com uma banda proibida ampla, apresentando uma banda proibida à temperatura ambiente de aproximadamente 3,4 eV. O seu elevado campo elétrico crítico e a baixa concentração de portadores intrínsecos tornam-no adequado para a conversão de energia, retificação elétrica e eletrónica de alta temperatura. A junção homojunção compacta p-GaN/n-GaN presente neste tutorial oferece uma forma clara de investigar como a barreira interna controla o fluxo de corrente e como as propriedades do material do GaN influenciam uma junção retificadora.
Neste tutorial, o OghmaNano é utilizado para analisar a retificação elétrica, o transporte de portadores por deriva e difusão, a curvatura das bandas de energia e a recombinação. As suas equações acopladas de Poisson e de continuidade de eletrões/lacunas calculam a corrente que flui através do díodo, juntamente com os perfis internos de portadores e eletrostáticos, em função da posição e da tensão aplicada. Isto estabelece uma ligação entre a característica J–V (corrente-tensão) e os processos físicos que ocorrem na junção.
As finas camadas dispostas verticalmente são exibidas na janela principal do dispositivo OghmaNano (??). A visualização é tridimensional, mas a simulação elétrica apresentada aqui é unidimensional, ao longo da direção de crescimento das camadas.
2. Criar uma nova simulação
Na janela principal do OghmaNano, selecione «Nova simulação». A primeira janela que aparece mostra as categorias de simulação disponíveis (??). Abra os exemplos III–V/GaN e selecione a demonstração da diode PN de GaN (??). A ordem exata dos exemplos pode variar entre as diferentes versões do OghmaNano.
Após carregar o exemplo, abra o editor de camadas na barra de ferramentas do dispositivo. Isso exibirá as regiões semicondutoras, suas espessuras e as entradas selecionadas do banco de dados de materiais. O dispositivo é uma junção homopolar PN, com as camadas p-GaN e n-GaN formando a junção eletricamente ativa.
3. Análise da estrutura do dispositivo, dos contactos e do processo de dopagem.
3.1 Estrutura em camadas
A estrutura elétrica é composta por três regiões de GaN. As camadas de p-GaN e n-GaN formam a junção retificadora; a camada de n⁺-GaN proporciona uma região com maior concentração de dopantes, adjacente ao cátodo. A espessura total do semicondutor é de 700 nm, excluindo o ar circundante utilizado para visualização. O editor de camadas é apresentado em ??.
| Camada | Espessura | Dopagem nominal | Função na área elétrica. |
|---|---|---|---|
| p-GaN (dopado com magnésio) | 250 nm | \(N_A=5\times10^{17}\,\mathrm{cm^{-3}}\) | Lado da junção PN com predominância de lacunas. |
| n-GaN (dopado com silício) | 350 nm | \(N_D=5\times10^{17}\,\mathrm{cm^{-3}}\) | Lado da junção PN com predominância de eletrões. |
| GaN dopado com n⁺ | 100 nm | \(N_D=1\times10^{18}\,\mathrm{cm^{-3}}\) | Região do lado do cátodo com maior concentração de impurezas. |
3.2 Contatos
O exemplo utiliza condições de contorno óhmicas ideais no ânodo do lado p e no cátodo do lado n. Essas condições de contorno permitem concentrar-se na resposta elétrica interna da junção. Em diodos de GaN fabricados, a resistência de contato e a ativação dos aceptores de magnésio influenciam a queda de tensão e a injeção de corrente; ambos são parâmetros úteis para estudos mais avançados.
3.3 Perfil de dopagem e região de depleção
Abra o editor de «Dopagem / Íons» na barra de ferramentas «Elétrica» para visualizar as distribuições de aceitadores e doadores (??). Uma única transição abrupta, que vai do material p-GaN dopado com aceitadores para o material n-GaN dopado com doadores, define a junção PN ativa. A região n⁺ introduz uma segunda etapa na concentração de doadores perto do cátodo.
No equilíbrio térmico, a junção desenvolve uma região de depleção à medida que os portadores de carga móveis difundem e deixam para trás os dopantes ionizados. Dentro da aproximação da região de depleção, a largura característica depende da permissividade, da concentração de dopagem em cada lado e do potencial suportado pela junção. Para uma junção abrupta convencional, a relação conhecida é:
\[ W \simeq \sqrt{\frac{2\varepsilon_s}{q}\left(\frac{1}{N_A}+\frac{1}{N_D}\right)(V_{\mathrm{bi}}-V)}. \]Aqui, \(V\) representa a tensão aplicada na junção, \(V_{\mathrm{bi}}\) é o potencial intrínseco e \(\varepsilon_s\) é a permissividade do GaN. Esta expressão fornece uma referência útil para interpretar o campo elétrico autoconsistente calculado pelo OghmaNano, especialmente ao analisar como a dopagem e a tensão aplicada afetam a junção.
4. Parâmetros elétricos e modelos de materiais
4.1 Modelo de material GaN Lua
As três camadas semicondutoras ativas eletricamente fazem referência ao modelo de material Inorganic/GaN/GaN.lua. O banco de dados de materiais elétricos fornece dados como a largura de banda, afinidade eletrônica, densidade de estados, permissividade, mobilidades e parâmetros de recombinação. No exemplo em JSON, os parâmetros db_overrides relevantes estão desativados, portanto, os valores utilizados devem ser aqueles retornados pelo modelo Lua, e não os valores inativos visíveis em outras partes do arquivo salvo. A largura de banda do GaN deve ser de aproximadamente 3,4 eV à temperatura ambiente; verifique isso nos perfis de banda exibidos ou nos parâmetros de material avaliados.
As capturas de tela das interfaces de código de material e editor de parâmetros estão agrupadas abaixo. Alguns dos nomes dos arquivos das capturas de tela originais ainda contêm as etiquetas históricas de GaAs/AlGaAs; a simulação ativa descrita neste tutorial utiliza GaN em todas as etapas. No software, verifique o caminho do material GaN no Lua, em vez de tentar inferir o material a partir do nome do arquivo da imagem.
4.2 Transporte de portadores
A OghmaNano resolve a equação de Poisson para semicondutores, juntamente com as equações de continuidade para eletrões e lacunas, utilizando expressões de corrente de deriva e difusão. As densidades de corrente de eletrões e lacunas dependem do campo local, da mobilidade e do gradiente de concentração de portadores:
\[ J_n=q\mu_n n E+qD_n\frac{dn}{dy},\qquad J_p=q\mu_p p E-qD_p\frac{dp}{dy}. \]Estas expressões descrevem dois mecanismos complementares: a derivação causada pelo campo elétrico interno e a difusão causada pelos gradientes de concentração de portadores de carga. Em estado estacionário, a interação entre estes mecanismos determina a corrente líquida do díodo. O modelo também determina os níveis quase-Fermi de eletrões e lacunas, cuja variação oferece uma perspetiva alternativa sobre as forças que impulsionam os portadores de carga.
4.3 Recombinação
Os eletrões e as lacunas injetados podem recombinar-se no interior do GaN. O exemplo inclui os processos de recombinação configurados no modelo do material, incluindo os termos de Shockley–Read–Hall (SRH) mediados por defeitos, quando ativados pelos dados efetivos do material. A expressão SRH idealizada habitual é:
\[ R_{\mathrm{SRH}}=\frac{np-n_i^2}{\tau_p(n+n_1)+\tau_n(p+p_1)}. \]A vida útil dos centros de recombinação (SRH) depende das densidades de armadilhas, das velocidades térmicas dos portadores de carga e das secções eficazes de captura. Estes parâmetros determinam a eficiência com que os portadores de carga injetados se recombinam e, portanto, influenciam a forma e a magnitude da curva característica corrente-tensão. A análise da taxa de recombinação espacial ajuda a identificar onde a perda de portadores de carga está mais concentrada dentro da junção.
4.4 Malha numérica
A malha elétrica cobre a estrutura de GaN de 700 nm, com 150 pontos ao longo da direção de crescimento, proporcionando um espaçamento médio de aproximadamente 4,7 nm. Este valor representa uma resolução inicial razoável para uma demonstração didática. O editor de malha (??) permite visualizar a distribuição dos nós e ajustá-la ao analisar gradientes acentuados na região de depleção.
5. Executar a simulação J–V.
O exemplo está configurado para uma variação de tensão de 0 V a +5 V, em incrementos de 0,025 V. Inicie a simulação com o comando «Executar simulação». Durante a variação, a janela de terminal exibe o progresso numérico e as informações sobre a convergência (??). Após a conclusão, o separador «Resultados» contém os resultados gerados, incluindo jv.csv e dados relacionados com o dispositivo.
Abra jv.csv para analisar a relação entre a densidade de corrente e a tensão aplicada. A visualização linear (??) e a visualização logarítmica (??) apresentam a mesma resposta elétrica de formas complementares. O gráfico linear mostra como a corrente aumenta em altas tensões de polarização direta; o gráfico logarítmico facilita a identificação do início da injeção e da resposta em baixas correntes.
5.1 Interpretação da curva de polarização direta.
Em baixas tensões de polarização direta, a barreira de potencial limita fortemente a injeção de portadores minoritários, de modo que a corrente resultante permanece baixa. À medida que a tensão aplicada aumenta, a barreira da junção diminui, os eletrões e as lacunas atravessam a junção PN com maior facilidade e a corrente aumenta rapidamente. Num gráfico semilogarítmico, uma região quase exponencial é frequentemente descrita localmente através da equação do díodo.
\[ J\simeq J_0\left[\exp\left(\frac{qV}{n_{\mathrm{id}}k_{\mathrm B}T}\right)-1\right], \]onde \(J_0\) é a corrente de saturação efetiva e \(n_{\mathrm{id}}\) é o fator de idealidade. Esta equação é uma ferramenta útil para a interpretação dos resultados, mas não é um modelo independente de todo o dispositivo simulado. Em particular, uma alteração na inclinação em valores de corrente mais elevados pode indicar resistência em série, efeitos de injeção ou alterações na recombinação. O aumento acentuado da corrente mostrado nas capturas de tela fornecidas é consistente com a retificação em polarização direta, mas um fator de idealidade específico só deve ser indicado após a definição de um intervalo de tensão adequado.
5.2 Análise do comportamento da junção sob polarização reversa.
A curva inicial de variação da corrente em função da tensão (J–V) varia de 0 V a +5 V e ilustra o desenvolvimento da condução direta. Para um exercício complementar, defina uma tensão inicial negativa e examine a alteração na curvatura da banda, na largura da região de depleção e na distribuição do campo elétrico sob polarização reversa. Isso permite uma comparação útil com a injeção direta e ilustra como a junção PN reage à polaridade da tensão aplicada.
6. Análise de imagens da simulação: bandas, densidades de portadores, correntes e recombinação.
O visualizador de imagens instantâneas fornece dados internos em pontos selecionados ao longo da varredura de tensão. As quatro imagens em cada grupo abaixo correspondem a posições selecionadas progressivamente na varredura salva; elas são identificadas como imagens instantâneas 0–3, em vez de apresentarem os valores exatos de tensão, porque a posição do controle deslizante, por si só, não define o valor de polarização. A tensão aplicada selecionada é exibida no visualizador e deve ser lida diretamente para uma interpretação quantitativa.
6.1 Bandas de condução e de valência; níveis de quasi-Fermi
Abra os ficheiros Ec.csv, Ev.csv, Fn.csv e Fp.csv a partir dos dados obtidos. As extremidades da banda de condução e da banda de valência descrevem a barreira de junção eletrostática. Os níveis quase-Fermi dos eletrões e das lacunas mostram como as populações de portadores se desviam do equilíbrio. Em condições de polarização nula, o nível de Fermi de equilíbrio deve ser comum em todo o dispositivo conectado, enquanto uma tensão direta aplicada cria uma separação crescente entre os potenciais eletroquímicos dos portadores.
Numa heterojunção de GaN, ambos os lados utilizam o mesmo semicondutor, pelo que a variação da banda na fronteira PN deve-se principalmente à curvatura eletrostática da banda. Compare a barreira perto da interface p–n em baixa polarização com o perfil após o aumento da injeção direta. O degrau doador n/n⁺ pode afetar a curvatura local sem alterar a banda de energia fundamental do GaN.
6.2 Densidades de portadores de elétrons e lacunas.
Os perfis de densidade de elétrons livres e lacunas livres armazenados (Q_nfree.csv e Q_pfree.csv) fornecem uma representação direta dos portadores majoritários e minoritários. Na primeira imagem, os elétrons predominam no lado do tipo n e as lacunas predominam no lado do tipo p, com mudanças acentuadas perto da região de depleção. Sob polarização direta, a injeção de portadores minoritários aumenta a população de portadores em ambos os lados da junção, alterando, eventualmente, a forma espacial dos perfis de densidade.
A escala logarítmica de densidade é particularmente útil ao comparar portadores minoritários e majoritários, que podem diferir em várias ordens de magnitude. Não confunda a concentração nominal de aceitadores de Mg com a concentração de lacunas livres sem considerar a ionização dos aceitadores; a ativação incompleta de Mg é uma característica importante do GaN tipo p experimental.
6.3 Densidades de corrente de eletrões e lacunas.
Trace os gráficos de Jn.csv e Jp.csv para analisar a distribuição da corrente entre os eletrões e as lacunas. Em estado estacionário, as correntes de eletrões e lacunas podem trocar valores ao longo do dispositivo à medida que os portadores se recombinam. De acordo com as convenções de sinal, a corrente elétrica total deve ser aproximadamente independente da posição quando não houver acumulação líquida de carga.
A captura de tela com baixa corrente apresenta características numéricas bem definidas. Essa estrutura pode surgir quando componentes de corrente pequenos são avaliados a partir de fluxos de deriva e difusão opostos ou quando uma transição abrupta não é resolvida adequadamente. Não se deve interpretá-la automaticamente como um pico de corrente físico. Avalie a convergência do solucionador e compare a corrente total que atravessa o dispositivo antes de fazer uma afirmação quantitativa.
6.4 Mecanismos de recombinação
As imagens que representam a recombinação combinam as taxas de recombinação livre-a-livre e as taxas relacionadas com o mecanismo SRH (R_nfree_to_pfree.csv e R_ss_srh.csv na sequência de capturas de tela guardada). Em níveis de injeção baixos, a recombinação mediada por defeitos pode concentrar-se em torno da junção, onde ambos os tipos de portadores estão presentes; em níveis de polarização direta mais elevados, a distribuição espacial da recombinação alarga-se à medida que as populações de portadores minoritários aumentam. A extensão e o tamanho relativo destes canais dependem dos parâmetros numéricos do material e não podem ser inferidos apenas a partir da altura de um pico.
Ao comparar os perfis de recombinação em diferentes polarizações diretas, é possível identificar os locais onde os eletrões e as lacunas injetadas interagem com maior intensidade. Esta informação espacial complementa a curva J–V e ajuda a relacionar as alterações na corrente com as populações de portadores de carga no interior da junção.
6.5 Perfis das taxas de geração de eletrões e lacunas.
As saídas das taxas de elétrons e lacunas (R_n.csv e R_p.csv) fornecem informações adicionais sobre o equilíbrio local de portadores de carga. Os seus sinais e componentes individuais devem ser interpretados de acordo com as convenções de saída da OghmaNano. Compare as características nítidas observadas em baixas tensões com as distribuições mais amplas observadas em tensões mais elevadas; a interpretação física deve concentrar-se nas tendências consistentes, e não na interpretação de cada pequena variação numérica como um mecanismo separado.
7. Exercícios e interpretação
8. Pressupostos do modelo e análises adicionais
A simulação utiliza um modelo de deriva-difusão unidimensional e em estado estacionário de uma estrutura compacta de semicondutor GaN de 700 nm, com condições de contorno óhmicas nos contactos. Esta formulação é particularmente útil para investigar como a dopagem da junção, os parâmetros do material e a recombinação afetam a estrutura de banda interna e a resposta J–V. Estudos de configurações completas de dispositivos de potência podem ser desenvolvidos com base neste modelo, incorporando a geometria do dispositivo, a resistência de contacto e o comportamento térmico.
A camada de p-GaN é caracterizada por uma densidade nominal de aceitadores de Mg. Como os aceitadores de Mg no GaN são relativamente profundos, a sua ativação elétrica e a densidade resultante de lacunas requerem atenção. Os detalhes dependem do modelo de material utilizado para as estatísticas de portadores e a ionização dos átomos dopantes. Da mesma forma, os defeitos na interface, as armadilhas no material, a mobilidade dos portadores e os coeficientes de recombinação influenciam a resposta elétrica calculada.
A varredura de 0 a 5 V oferece uma visão detalhada da condução direta, desde o início da injeção de portadores até o regime de corrente mais elevada. Estender a varredura para a polarização reversa oferece outra direção útil para analisar o comportamento da região de depleção e do campo elétrico. Estudos mais avançados de dispositivos de potência podem incluir, adicionalmente, os efeitos de campos elevados e o bloqueio de tensão dependente da geometria.
Para obter previsões quantitativas de nível de pesquisa, verifique as propriedades avaliadas do material GaN, analise o resíduo do solucionador e a continuidade da corrente, refine a malha espacial e compare com medições experimentais ou dados da literatura para uma estrutura com correspondência precisa. Essas verificações ajudam a expandir o modelo inicial para uma comparação quantitativa com dispositivos fabricados.
9. Próximos passos
Após concluir a análise inicial da curva J–V, o próximo passo útil é variar a dopagem dos tipos p e n e comparar a deformação resultante da banda e o comportamento da região de depleção. Um estudo complementar consiste em introduzir uma região central com baixa dopagem e comparar as arquiteturas GaN PN e PIN. Essas alterações demonstram como o design das camadas afeta a eletrostática, a injeção direta e o potencial de suporte à tensão reversa.