Simulação de Diodo PIN de GaN: Modelos de Deriva-Difusão, Curvas J–V e Transporte em Materiais com Banda Proibida Ampla.
1. Introdução
Um díodo PIN de nitreto de gálio (GaN) é um dispositivo semicondutor de banda proibida ampla, com uma região do tipo p, uma região comparativamente pouco dopada (nominalmente intrínseca) e uma região do tipo n. Ao contrário de uma junção PN comum, um díodo PIN introduz deliberadamente uma região de derivação substancial entre os terminais dopados. Sob polarização reversa, esta região pode suportar um campo elétrico distribuído; sob polarização direta, os eletrões e as lacunas injetados permitem que a estrutura conduza eletricidade. Estas propriedades tornam a arquitetura PIN relevante para os princípios subjacentes da retificação de potência, bloqueio de tensão reversa e engenharia de dispositivos de banda proibida ampla.
O GaN possui uma largura de banda proibida à temperatura ambiente de aproximadamente 3,4 eV e pode suportar campos elétricos elevados em comparação com o silício. No entanto, os diodos de potência GaN práticos exigem um controlo cuidadoso da espessura da camada de deriva e da dopagem de fundo, da resistência de contacto, dos defeitos cristalinos e da terminação da borda. A região intrínseca de GaN de 800 nm neste exemplo é propositadamente compacta: isso facilita a análise da eletrostática central e da injeção de portadores, mas, por si só, não representa um retificador de potência de alta tensão com as características especificadas.
Este tutorial utiliza as equações acopladas de Poisson e de deriva-difusão do OghmaNano para resolver um diodo PIN de GaN unidimensional. Em vez de tratar a corrente como uma única grandeza de circuito, o modelo analisa o perfil interno da banda, os níveis de quasi-Fermi, as concentrações de eletrões e lacunas e as taxas de recombinação à medida que a polarização é alterada. O objetivo é estabelecer uma ligação entre uma curva J–V intuitiva, que representa a retificação, e os processos físicos que ocorrem no interior do semicondutor.
2. Criar uma nova simulação
Abra o OghmaNano e selecione «Nova simulação». Escolha as demonstrações de semicondutores III–V e abra o exemplo do díodo PIN de GaN. As categorias de simulação disponíveis e a janela de seleção de demonstrações são apresentadas nas Figuras ?? e ??. Se a demonstração já estiver instalada, ao abri-la, deverá ser exibida a estrutura de camadas na janela principal do dispositivo.
A vista principal da simulação (??) apresenta uma representação tridimensional do dispositivo. Embora a visualização seja em 3D, o exemplo numérico apresentado aqui resolve apenas o transporte elétrico ao longo da direção do crescimento da camada. Abra o editor de camadas para visualizar a espessura real e a atribuição de materiais para cada região.
3. Estrutura em camadas e dopagem
O protótipo finalizado utiliza a seguinte estrutura semicondutora eletricamente ativa. A camada intrínseca de 800 nm é a configuração utilizada para a simulação apresentada; versões mais antigas do exemplo podem conter uma camada intrínseca de 400 nm e devem ser atualizadas antes de se tentar reproduzir estes resultados.
| Camada | Espessura | Concentração nominal de impurezas dopantes | Objetivo/Finalidade |
|---|---|---|---|
| p-GaN | 200 nm | Concentração de aceitadores: 5 × 1017 cm−3 | injeção no lado p e região do ânodo |
| GaN intrínseco | 800 nm | Sem recurso intencional a substâncias proibidas. | Região de derivação e região de depleção com baixa concentração de impurezas. |
| n-GaN | 200 nm | Concentração de doadores: 5 × 1017 cm−3 | região de junção do tipo n |
| GaN dopado com n⁺ | 100 nm | Concentração de impurezas doadoras: 1 × 1018 cm−3 | Região de contato traseiro |
A camada ativa de semicondutores tem uma espessura de 1,3 µm. Uma região de ar separada de 200 nm pode aparecer no editor de camadas para a visualização do dispositivo e para a definição da geometria ótica; não faz parte da estrutura PIN eletricamente ativa. As quatro regiões de semicondutores utilizam GaN, pelo que o objetivo das interfaces aqui é alterar a dopagem, e não a composição da liga de semicondutores.
Abra o editor de Contatos (veja a figura abaixo) para verificar o ânodo e o cátodo e, em seguida, abra o editor de Dopagem (veja a figura abaixo). A região p dopada com aceitadores e as regiões n/n⁺ dopadas com doadores devem aparecer em lados opostos de uma ampla região central não dopada. O centro da estrutura não é uma heterojunção de materiais separada: a sua principal característica distintiva é a baixa dopagem intencional.
Na aproximação idealizada de depleção, uma região com baixa carga líquida pode suportar um campo elétrico que se estende por grande parte da sua espessura. A tensão suportada por uma região está relacionada com a integral espacial do seu campo elétrico:
\[V_{\mathrm{region}}=-\int_{y_0}^{y_1} E(y)\,\mathrm{d}y.\]Dessa forma, uma região mais ampla e com uma concentração menor de impurezas pode distribuir um potencial reverso aplicado por uma distância maior. Esta é uma das principais razões pelas quais as estruturas PIN são estudadas para aplicações de bloqueio de tensão. Isso não significa que a tensão de rutura possa ser obtida simplesmente multiplicando a espessura por um campo crítico assumido: os dispositivos reais apresentam campos não uniformes, carga de fundo, contactos e concentração de campo nas bordas.
4. Modelos elétricos, contactos e configuração numérica.
4.1 Parâmetros do material GaN
Abra o editor de parâmetros elétricos para verificar as configurações de transporte e recombinação utilizadas para cada região de GaN. A simulação está configurada para carregar as propriedades do GaN a partir do banco de dados de materiais OghmaNano, com as opções de substituição do banco de dados desativadas. Esta é a opção preferível em vez de atribuir mobilidades, energias de banda ou afinidades eletrónicas arbitrárias no ficheiro JSON do dispositivo. Verifique os dados dos materiais carregados se precisar de uma comparação quantitativamente reproduzível com os resultados experimentais.
No caso do GaN, as principais considerações incluem a grande largura de banda proibida, a mobilidade dos eletrões e das lacunas, as estatísticas dos portadores de carga livres, a recombinação radiativa e a recombinação assistida por armadilhas. A concentração nominal de aceitadores exibida no editor de dopagem não deve ser automaticamente identificada com a densidade de lacunas livres: os aceitadores de magnésio no GaN tipo p não estão completamente ionizados à temperatura ambiente. É importante considerar em que medida o modelo elétrico trata a ionização dos dopantes antes de atribuir uma resistência de contacto física ou determinar uma densidade de portadores medida.
4.2 Contatos e condições de contorno
A demonstração utiliza contactos óhmicos ideais para o ânodo do lado p e o cátodo do lado n. Estes simplificam a interpretação, evitando um segundo elemento retificador na interface metal/semicondutor. Em dispositivos GaN fabricados, especialmente nos contactos p-GaN, obter um comportamento óhmico de baixa resistência pode ser um desafio e envolve a estrutura metálica, a preparação da superfície e a dopagem.
4.3 Rede elétrica
O editor de malha elétrica (mostrado abaixo) apresenta a discretização ao longo da espessura do semicondutor. Ao alterar a região intrínseca de 400 nm para 800 nm, certifique-se de que a malha cobre a pilha ativa atualizada de 1,3 µm. Uma escolha inicial de aproximadamente 250 a 260 pontos proporciona um espaçamento da ordem de 5 nm; o software pode refazer a malha automaticamente, dependendo da configuração. Para um estudo mais quantitativo, repita a simulação com uma malha mais fina e verifique se a curva J–V e os perfis internos permanecem efetivamente inalterados.
4.4 Transporte e recombinação
O OghmaNano associa a equação de Poisson eletrostática às equações de continuidade de eletrões e lacunas. Em uma dimensão, as equações em estado estacionário podem ser escritas de forma esquemática como:
\[\frac{\mathrm{d}}{\mathrm{d}y}\!\left(\varepsilon\frac{\mathrm{d}\phi}{\mathrm{d}y}\right)=-\rho,\qquad \frac{\mathrm{d}J_n}{\mathrm{d}y}=q(R-G),\qquad \frac{\mathrm{d}J_p}{\mathrm{d}y}=-q(R-G),\]onde \(\phi\) é o potencial eletrostático, \(\rho\) a densidade local de carga espacial, \(J_n\) e \(J_p\) as densidades de corrente de eletrões e lacunas, \(G\) a taxa de geração e \(R\) a taxa líquida de recombinação. Com estas convenções de sinal, a densidade de corrente total \(J_n+J_p\) é constante num dispositivo unidimensional, em estado estacionário e com condições de convergência.
O modelo de recombinação pode incluir contribuições diretas de banda para banda (recombinação livre-livre) e de Shockley–Read–Hall (SRH). A importância relativa de cada um depende da qualidade do material, do nível de injeção e dos parâmetros do modelo. A banda de condução direta do GaN torna a recombinação radiativa fisicamente possível, mas a sua presença num modelo não implica que um determinado dispositivo representado seja um LED eficiente.
5. Executar a simulação J–V em condições de baixa luminosidade.
Defina a variação da tensão entre −2,0 V e +5,0 V, em incrementos de 0,025 V. Uma tensão negativa polariza a estrutura PIN em polarização reversa e uma tensão positiva a polariza em polarização direta. A configuração de polarização pretendida é:
"Vstart": -2.0,
"Vstop": 5.0,
"Vstep": 0.025Na janela principal, selecione «Executar simulação». Durante o cálculo, utilize o separador «Terminal» para monitorizar o progresso e a convergência. Após a conclusão, abra o separador «Resultados» e selecione jv.csv para apresentar graficamente a densidade de corrente em função da tensão aplicada. O visualizador de resultados, o terminal do solucionador e a curva J–V resultante são apresentados abaixo.
Em polarização inversa, a corrente simulada é pequena em relação ao seu valor em polarização direta. À medida que a polarização positiva aumenta, a injeção de portadores aumenta rapidamente e o diodo entra em funcionamento. Em tensões de polarização direta mais elevadas, a curva torna-se mais linear num eixo de densidade de corrente linear, indicando que a corrente é cada vez mais influenciada pelo transporte resistivo e outras não idealidades na estrutura simulada. Uma inclinação que se torna mais suave num eixo de corrente logarítmica não deve ser interpretada automaticamente como uma saturação acentuada da corrente: um aumento de corrente aproximadamente linear também se torna progressivamente mais suave num gráfico logarítmico.
Para uma comparação simples com a teoria dos diodos, a característica de polarização direta, que varia de baixa a moderada, pode ser relacionada com:
\[J\approx J_0\left[\exp\!\left(\frac{qV}{nk_{\mathrm B}T}\right)-1\right],\]onde \(J_0\) é uma corrente de saturação efetiva e \(n\) um fator de idealidade efetivo. Um diodo PIN pode apresentar desvios significativos em relação a um regime ideal de Shockley único, pois a recombinação na região de depleção, a injeção na camada intrínseca, a injeção em níveis elevados e a resistência em série podem influenciar a corrente. Portanto, a expressão simples deve ser considerada como uma aproximação diagnóstica, e não como a equação fundamental resolvida pelo OghmaNano.
A capacidade de bloquear a tensão inversa não é equivalente à rutura por avalanche observada. Um teste que varia a tensão até −2 V apenas avalia o comportamento em pequena tensão inversa. A tensão na qual um dispositivo de potência de GaN fabricado entra em rutura depende da física em campos elevados, dos defeitos do material, dos caminhos de fuga, da geometria e da terminação das bordas, e não pode ser determinada apenas com base nesta demonstração.
6. Análise dos resultados da simulação interna.
Após a varredura J–V, o visualizador de imagens da OghmaNano permite que os resultados, resolvidos espacialmente, sejam representados graficamente em pontos de polarização selecionados. As capturas de tela abaixo ilustram quatro posições na varredura armazenada, começando pela mais antiga e terminando na mais recente selecionada. Os seus nomes (-0 a -3) identificam as posições das capturas de tela, e não as tensões. Consulte a tensão selecionada real no visualizador de imagens ao realizar comparações quantitativas.
6.1 Limites das bandas de energia e níveis de quasi-Fermi
Abra o visualizador de imagens e apresente os gráficos de Ec.csv, Ev.csv, Fn.csv e Fp.csv em conjunto. Estes correspondem às extremidades das bandas de condução e de valência, bem como aos níveis quase-Fermi de eletrões e lacunas. As figuras ?? a ?? mostram a sua evolução ao longo da variação da polarização.
Na polarização inversa, espera-se que grande parte da variação do potencial ocorra na região intrínseca com baixa concentração de impurezas, produzindo uma inclinação da banda relativamente extensa. A polarização direta reduz a barreira de injeção efetiva e cria uma separação das quase-níveis de Fermi de elétrons e lacunas, que depende da polarização. Os gradientes espaciais desses níveis são importantes porque determinam o fluxo de corrente; a curvatura da banda e a divisão das quase-níveis de Fermi devem ser consideradas em conjunto, em vez de interpretadas de forma independente.
Como parte de uma verificação física básica, a separação entre \(E_c\) e \(E_v\) nas camadas de GaN deve ser consistente com a banda de energia do GaN presente na base de dados de materiais, e não um valor padrão residual de outro semicondutor. Uma vez que este é um dispositivo PIN de material homogêneo, não se espera que ocorram descontinuidades abruptas na banda de energia induzidas pela composição nas transições p/i/n.
6.2 Densidades de eletrões e lacunas.
Trace os gráficos de Q_nfree.csv e Q_pfree.csv (ou os resultados da densidade de portadores de elétrons e lacunas apresentados pela versão instalada). Em baixas taxas de injeção, as densidades dos portadores majoritários concentram-se nas respetivas regiões dopadas. A região central não dopada contém concentrações de portadores de equilíbrio muito mais baixas. Em forte polarização direta, os portadores injetados ocupam uma fração maior da região intrínseca e a recombinação pode aumentar substancialmente.
A mudança abrupta na densidade nominal de impurezas não implica um perfil idêntico e abrupto de portadores de carga livres. A densidade de portadores é determinada de forma autoconsistente a partir das energias da banda, níveis de Fermi, densidade de estados e neutralidade de carga. Isto é especialmente relevante no lado do p-GaN, caso o modelo do material utilize um tratamento simplificado da ionização incompleta dos aceitadores.
6.3 Recombinação na região intrínseca
Visualize os gráficos de recombinação espacial nos mesmos pontos de polarização. As figuras de recombinação deste exemplo mostram como a região de sobreposição e recombinação significativas entre elétrons e lacunas se desenvolve à medida que a injeção direta aumenta. Interprete a altura e a largura de cada perfil representado em conjunto: um pico local elevado não é necessariamente o que mais contribui para a taxa de recombinação total do dispositivo.
Para um processo direto de transição entre bandas, a taxa de emissão radiativa resultante segue, geralmente, a seguinte forma:
\[R_{\mathrm{rad}}=B\left(np-n_i^2\right),\]onde \(B\) é o coeficiente de recombinação do material. A recombinação SRH depende dos parâmetros de defeito e pode apresentar uma distribuição espacial diferente. Para comparar os mecanismos, represente-os graficamente separadamente, caso ambos estejam disponíveis e ativados no dispositivo carregado. Não presuma que um gráfico de recombinação total sem identificação corresponde a um único mecanismo microscópico.
6.4 Taxas de recombinação/continuidade de eletrões e lacunas.
Os gráficos adicionais que mostram a variação das taxas de elétrons e lacunas são úteis para analisar os pontos onde as concentrações de portadores de carga se alteram mais rapidamente, através da recombinação e da injeção. As quatro capturas de tela abaixo mostram como esses perfis calculados se modificam em função da polarização. Se a taxa exibida for positiva ou negativa, ou se contiver contribuições tanto de geração quanto de recombinação, utilize a legenda do gráfico e o cabeçalho do arquivo de dados para determinar a correta interpretação física antes de comparar as magnitudes.
6.5 A situação atual da conservação como modelo de referência.
Para qualquer ponto de polarização em regime permanente, represente graficamente a densidade de corrente de eletrões \(J_n\) e a densidade de corrente de lacunas \(J_p\), caso esses valores estejam disponíveis. Cada componente pode variar com a posição à medida que os portadores de carga se recombinam, mas o total \(J_n+J_p\) deve ser espacialmente constante num dispositivo 1D convergente, sem termos adicionais de acumulação de carga. Esta é uma verificação particularmente útil quando uma pequena ondulação aparente num traçado de baixa corrente pode, de outra forma, ser confundida com um efeito físico.
7. Exercícios sugeridos
8. Premissas e limitações físicas
Este exemplo tem como objetivo servir de simulação educativa em larga escala do transporte de portadores em materiais GaN. Apresenta uma primeira análise, com significado físico, da eletrostática de dispositivos PIN e da injeção de portadores, mas não se trata de um modelo comercial calibrado. Em particular, os resultados devem ser interpretados tendo em conta as seguintes ressalvas.
- As características ideais dos contactos: as barreiras de contacto metal/GaN e a resistência prática do contacto p-GaN não são explicitamente validadas.
- Dopagem nominal: a ativação de aceitadores do tipo p e a carga residual de fundo em GaN nominalmente intrínseco podem alterar a distribuição física da carga.
- Ruptura em campos magnéticos elevados: a varredura da curva J–V não determina a ruptura em avalanche, a tensão de bloqueio limitada pela terminação ou a fiabilidade sob alta polarização reversa.
- Transporto unidimensional: o modelo não considera o efeito de concentração da corrente lateral, as perdas nas bordas, as placas de campo e a geometria de terminação do dispositivo.
- Defeitos e recombinação: previsões realistas e quantitativas de fugas e recombinação exigem parâmetros calibrados para os defeitos no volume do material e na interface.
- Efeitos térmicos e dinâmicos: o aquecimento interno, as transições de comutação e os efeitos de armazenamento de carga dos portadores minoritários devem ser avaliados separadamente, quando aplicável.
Estas limitações não diminuem o valor do exemplo enquanto ponto de partida. A estrutura em camadas e a curva J–V ilustram o que acontece quando uma região de derivação com baixa concentração de impurezas é inserida num díodo semicondutor de banda proibida ampla. Além disso, fornecem uma base de referência reprodutível para testar modelos mais avançados de transporte e recombinação.
9. Próximos passos
Após analisar o diodo PIN de GaN, compare a sua distribuição de campo com a de uma junção PN de GaN convencional e, em seguida, investigue como a espessura da camada intrínseca e a dopagem residual influenciam a retificação. Trabalhos mais avançados poderiam introduzir estados de armadilha calibrados, contactos não ideais e mecanismos de transporte em campos elevados para examinar os limites do modelo simples de deriva e difusão.