Simulação de um Diodo de Junção PN em InP: Deriva-Difusão, Retificação e Recombinação de Portadores.
1. Introdução: uma junção retificadora de InP.
O fosfeto de índio (InP) é um semicondutor do grupo III–V com banda de energia direta, utilizado como base para dispositivos eletrónicos de alta velocidade e tecnologias fotónicas integradas. Um díodo de junção PN é uma das suas estruturas de dispositivo mais simples e instrutivas: uma região do tipo p encontra-se com uma região do tipo n, criando um campo elétrico interno e uma barreira de potencial que controla o movimento de eletrões e lacunas. O díodo conduz fortemente sob polarização direta e suprime a injeção de portadores sob polarização inversa, produzindo retificação elétrica.
O InP apresenta uma largura de banda proibida à temperatura ambiente de aproximadamente 1,34 eV e uma mobilidade eletrónica relativamente elevada. Estas propriedades tornam o transporte em junções de InP relevante para a eletrónica de semicondutores compostos e para a compreensão das interfaces de semicondutores utilizadas em plataformas de dispositivos baseadas em InP. O exemplo apresentado neste tutorial demonstra o funcionamento elétrico de uma junção homopolar compacta p-InP/n-InP com uma camada de contacto do lado n⁺.
Em vez de ajustar a curva J–V a uma equação ideal de díodo, a OghmaNano resolve as equações acopladas de Poisson e de continuidade de deriva-difusão. Os dados resultantes revelam tanto a corrente nos terminais quanto a distribuição interna das energias das bandas de condução e de valência, as densidades de eletrões e lacunas, os níveis de quasi-Fermi e a recombinação. Isso permite-nos explicar por que o díodo retifica a corrente, em vez de simplesmente observar que o faz.
A janela principal de simulação, mostrada em ??, apresenta o dispositivo InP com camadas verticais. Embora a geometria seja exibida em três dimensões, o problema elétrico resolvido aqui é unidimensional, ao longo da direção de crescimento.
2. Criar uma nova simulação
Na janela principal do OghmaNano, selecione «Nova simulação». A janela de seleção de simulações será aberta, mostrando as categorias de exemplos disponíveis (??). Abra as demonstrações III–V/InP e escolha o exemplo de díodo PN de InP (??). A localização exata do exemplo pode variar entre as diferentes versões do OghmaNano.
Após abrir o exemplo, selecione o separador «Estrutura do dispositivo» e abra o editor de camadas. Verifique os nomes das camadas, as suas espessuras e os modelos de materiais elétricos associados antes de executar a simulação. As secções que se seguem explicam a contribuição de cada região para a resposta geral do díodo.
3. Análise da estrutura do dispositivo, dos contactos e do processo de dopagem.
3.1 Estrutura da camada semicondutora
O dispositivo elétrico contém três camadas de InP, com uma espessura total do semicondutor de 700 nm. As camadas p-InP e n-InP formam a junção homojunção retificadora, enquanto a região n⁺-InP proporciona uma região do lado do cátodo com uma concentração de dopagem mais elevada. A atribuição das camadas pode ser visualizada no editor (??).
| Camada | Espessura | Dopagem nominal | Função/Papel |
|---|---|---|---|
| p-InP | 250 nm | \(N_A=5\times10^{17}\,\mathrm{cm^{-3}}\) | Região com predominância de lacunas e lado do ânodo. |
| InP dopado com nitrogênio | 350 nm | \(N_D=5\times10^{17}\,\mathrm{cm^{-3}}\) | Região com predominância de eletrões que forma a junção. |
| InP dopado com n⁺ | 100 nm | \(N_D=1\times10^{18}\,\mathrm{cm^{-3}}\) | Camada do lado do cátodo com maior concentração de dopantes. |
A região de ar de 200 nm mostrada no editor de estrutura faz parte da representação visual e da geometria ótica circundante, e não da pilha elétrica do semicondutor. Consequentemente, a sua espessura não é incluída na espessura total do dispositivo elétrico, que é de 700 nm.
3.2 Dopagem e eletrostática intrínseca
Em «Open Electrical», selecione «Doping / Ions» para visualizar as concentrações de impurezas doadoras e aceitadoras (??). A variação líquida da dopagem inverte o sinal na interface entre o material p-InP e o material n-InP, aproximadamente a 250 nm da superfície da pilha de semicondutores eletricamente ativa. A segunda variação, entre o material n-InP e o material n⁺-InP, modifica a densidade de eletrões em direção ao cátodo, mas não cria outra junção PN.
No estado de equilíbrio, os eletrões tendem a difundir-se da região do tipo n para a região do tipo p, enquanto as lacunas difundem-se na direção oposta. Os átomos dopantes ionizados que não estão associados a outros átomos geram um campo que se opõe a uma maior difusão. A região de depleção resultante e a barreira de potencial interna formam a base da resposta de retificação.
\[ \frac{d}{dy}\left(\varepsilon\frac{d\phi}{dy}\right)=-\rho(y), \qquad \rho(y)=q\,[p-n+N_D^+-N_A^-]. \]Aqui, \(\phi\) representa o potencial eletrostático, \(\varepsilon\) a permissividade dielétrica e \(\rho\) a densidade de carga local. Os dopantes ionizados e os portadores de carga móveis determinam o campo elétrico local, enquanto a tensão aplicada modifica a barreira efetiva para a injeção de portadores.
3.3 Contatos elétricos
O modelo utiliza uma barreira ôhmica com predominância de lacunas no ânodo e uma barreira ôhmica com predominância de eletrões no cátodo. Estas condições de contorno são adequadas para estudar a resposta da junção sem introduzir uma barreira de injeção adicional nos contactos. Nos dispositivos fabricados, a metalização e o processamento dos contactos determinam o quão próximos os contactos se aproximam de um comportamento ôhmico ideal.
4. Modelos de materiais e parâmetros de transporte
4.1 Modelos de materiais para simulação da Lua: de onde vêm os parâmetros do InP.
Uma das características importantes deste exemplo de díodo InP é que os parâmetros dos materiais elétricos são definidos utilizando modelos de materiais Lua. Isso estabelece uma ligação clara entre as propriedades físicas de um semicondutor e a solução de deriva-difusão. As regiões p-InP, n-InP e n⁺-InP referenciam todas a mesma definição de material, Inorganic/InP/InP.lua. O seu comportamento elétrico diferente resulta principalmente do perfil de dopagem, e não da atribuição de três conjuntos distintos de constantes de material.
Para analisar o modelo, abra a estrutura do dispositivo → Parâmetros elétricos (o editor de materiais elétricos) e selecione uma das camadas de InP. O separador Editor de código exibe o script Lua do material associado a essa região; o separador Parâmetros exibe as propriedades do material calculadas. Estas duas visualizações devem ser usadas em conjunto. O código define como as propriedades são obtidas, enquanto os valores calculados revelam o que o solucionador elétrico recebe efetivamente.
Inorganic/InP/InP.lua.As imagens em ?? e ?? mostram como se acede ao editor de código. A configuração da camada da simulação é a referência para o modelo de material atribuído: em todos os exemplos, cada região ativa está associada ao ficheiro Lua do InP. Os nomes dos ficheiros das imagens são herdados da sequência de captura de ecrã e não determinam qual o semicondutor que está a ser simulado.
Na simulação fornecida, os interruptores db_overrides estão definidos para false para os parâmetros-chave de transporte e de banda. Assim, o OghmaNano é configurado para obter os valores correspondentes da base de dados de materiais/modelo Lua, em vez de utilizar os valores predefinidos que estão registados noutro local no ficheiro JSON. Isto é importante ao analisar a configuração: os valores efetivos apresentados no separador «Parâmetros» são mais informativos do que os valores predefinidos armazenados num campo de substituição inativo.
4.2 Estrutura de bandas, densidade de estados e mobilidade.
Os primeiros parâmetros a verificar no separador «Parâmetros avaliados» são a largura de banda \(E_g\) e a afinidade eletrónica \(\chi\). A largura de banda define a separação entre as bandas de condução e de valência, enquanto a afinidade eletrónica indica a posição da extremidade da banda de condução em relação a uma referência de energia do vácuo. No caso do InP, a largura de banda à temperatura ambiente é de aproximadamente \(1.34\,\mathrm{eV}\). Uma verificação inicial útil é determinar se a separação entre \(E_c\) e \(E_v\) no diagrama de bandas resultante corresponde à largura de banda avaliada.
As densidades de estados efetivas das bandas de condução e de valência, \(N_c\) e \(N_v\), influenciam a relação entre a concentração de portadores e a posição do nível de quasi-Fermi. Consequentemente, contribuem para as populações de equilíbrio de eletrões e lacunas, bem como para a resposta à injeção sob polarização direta. A permissividade relativa \(\varepsilon_r\) entra na equação de Poisson e influencia o potencial eletrostático e a largura da região de depleção. Estas quantidades devem ser analisadas como um conjunto consistente, e não como parâmetros de ajuste numérico independentes.
As mobilidades do eletrão e do buraco, \(\mu_n\) e \(\mu_p\), determinam a resposta ao desvio dos portadores de carga num campo elétrico; a difusão depende, adicionalmente, dos gradientes de densidade de portadores. No modelo unidimensional, as densidades de corrente são expressas da seguinte forma:
\[ J_n=q\mu_n nE+qD_n\frac{dn}{dy},\qquad J_p=q\mu_p pE-qD_p\frac{dp}{dy}. \]Aqui, \(n\) e \(p\) representam as densidades locais de eletrões e lacunas, \(E\) é o campo elétrico e \(D_n\) e \(D_p\) são os coeficientes de difusão. As duas componentes de corrente podem variar consideravelmente ao longo do díodo, mesmo quando a corrente total em estado estacionário é aproximadamente independente da posição. É por isso que as secções seguintes analisam \(J_n\) e \(J_p\) em conjunto com os perfis de densidade de portadores.
4.3 Como verificar e comparar as configurações dos materiais
No editor de dispositivos elétricos, selecione o material p-InP e examine o seu código Lua e os parâmetros calculados. Em seguida, altere para o material n-InP e, por fim, para o material n⁺-InP. Compare o modelo de material atribuído, a largura da banda proibida, a afinidade eletrónica, as mobilidades, as densidades de estados efetivas, a permissividade e os coeficientes de recombinação. Para esta junção homojunção, o mesmo modelo de InP é utilizado nas três regiões; a principal alteração intencional entre as regiões é a sua concentração de aceitadores ou doadores.
É importante distinguir duas ideias: o modelo de material descreve as propriedades do InP, enquanto o perfil de dopagem estabelece quais os portadores de carga que predominam localmente e onde se encontra a junção. Esta distinção torna-se especialmente útil ao passar de uma junção homopolar PN para uma heterojunção baseada em InP, onde as diferenças entre os modelos de material podem, por si só, gerar descontinuidades na banda de energia. Nesta junção, a curvatura das bandas de condução e de valência, observada nas imagens, reflete principalmente a distribuição eletrostática na transição p/n.
Para garantir a reprodutibilidade do estudo, mantenha o modelo de material Lua fixo ao alterar as concentrações de dopagem. Por outro lado, se estiver a investigar a sensibilidade do resultado em relação a um coeficiente de mobilidade ou recombinação, altere esse parâmetro do material de forma intencional e compare a curva J–V resultante, as densidades de portadores e os perfis de corrente interna. O editor de parâmetros e as visualizações do código Lua estabelecem a ligação entre as premissas de entrada físicas e a saída simulada.
4.4 Parâmetros de recombinação
As regiões de InP incluem a recombinação de Shockley–Read–Hall (SRH) e um canal de recombinação de portadores livres, conforme o modelo de material avaliado. A recombinação de SRH descreve a captura de eletrões e lacunas facilitada por defeitos; a recombinação de portadores livres é particularmente relevante no InP de banda proibida direta, pois as transições eletrão-lacuna podem ocorrer de forma radiativa. No presente exercício, ambos os processos são úteis como mecanismos de perda de portadores que modificam a corrente elétrica.
\[ R_{\mathrm{rad}}=B\left(np-n_i^2\right), \qquad R_{\mathrm{SRH}}=\frac{np-n_i^2}{\tau_p(n+n_1)+\tau_n(p+p_1)}. \]O coeficiente \(B\), os tempos de vida das armadilhas \(\tau_n\) e \(\tau_p\) e as densidades auxiliares \(n_1\) e \(p_1\) determinam a magnitude de cada processo. O principal aspeto a analisar neste tutorial é como a distribuição espacial da recombinação se altera à medida que a polarização faz com que mais eletrões e lacunas se movam para a mesma região.
4.5 Malha e solução numérica
O cálculo elétrico utiliza uma malha unidimensional de 150 pontos ao longo de 700 nm, o que corresponde a um espaçamento médio de aproximadamente 4,7 nm. As equações não lineares acopladas são resolvidas utilizando o método de iteração de Newton e o solucionador linear esparso UMFPACK. Esta malha constitui um ponto de partida útil para analisar a junção e as regiões quase neutras.
A saída inclui informações sobre a convergência e perfis de dispositivos dependentes do desvio. Uma boa primeira verificação consiste em analisar a suavidade da curva J–V e, em seguida, confirmar se as correntes de eletrões e lacunas obedecem à lei geral de conservação da corrente em estado estacionário. O editor de malha é apresentado em ??.
5. Executar a simulação e interpretar a resposta da curva J–V.
Volte para a janela principal e selecione «Executar simulação». O OghmaNano aplica diferentes valores de tensão e regista a solução obtida em cada ponto. O ficheiro de exemplo guardado utiliza uma variação de tensão de 0 a 3 V, com incrementos de 25 mV. O terminal regista o progresso do solucionador durante a variação, enquanto o separador «Saída» permite aceder aos dados resultantes.
Quando a simulação estiver concluída, os resultados estarão disponíveis na aba «Resultados» (??). Abra jv.csv para gerar um gráfico da densidade de corrente em função da tensão aplicada. As capturas de tela incluem tanto uma representação linear da curva J–V (??) quanto uma representação logarítmica da corrente (??).
5.1 Compreensão das características
Em baixas tensões aplicadas, a barreira de junção interna restringe fortemente a injeção de portadores de carga. Ao aumentar a polarização direta, essa barreira diminui, permitindo que os eletrões e as lacunas entrem mais facilmente no lado oposto, o que aumenta a corrente total. No regime intermédio, a corrente é frequentemente aproximadamente exponencial em relação à tensão; em correntes mais elevadas, o transporte e as perdas resistivas podem modificar a inclinação da curva.
\[ J\approx J_0\left[\exp\!\left(\frac{qV}{nk_{\mathrm B}T}\right)-1\right]. \]Esta expressão da «diodo ideal» fornece uma referência útil para analisar um intervalo aproximadamente exponencial, mas não é imposta pelo modelo de deriva-difusão. O fator de idealidade \(n\) e o fator de pré-exponencial \(J_0\) devem ser determinados apenas num intervalo de polarização adequado; podem variar em função da recombinação, da resistência em série e do comportamento do contacto.
O gráfico logarítmico ajuda a identificar o início de uma corrente significativa e as alterações na inclinação. O gráfico linear facilita a interpretação da magnitude e da forma da corrente elevada. Para obter uma imagem mais clara do funcionamento do diodo, analise ambos os gráficos em conjunto com os perfis da banda interna e da densidade de portadores.
6. Funcionamento interno: como o díodo funciona.
O diretório «snapshots» contém perfis espaciais guardados em pontos de polarização selecionados durante a varredura J–V. As capturas de ecrã fornecidas mostram quatro posições ao longo dessa varredura. Elas ilustram uma progressão desde uma baixa injeção até uma polarização direta mais forte. As quatro imagens em cada grupo correspondem às mesmas posições selecionadas no controlo deslizante, permitindo que as diferentes grandezas físicas sejam comparadas diretamente.
6.1 Bandas de condução e de valência e níveis de quasi-Fermi.
Abra o visualizador de imagens e selecione Ec.csv, Ev.csv, Fn.csv e Fp.csv. Estes representam, respetivamente, a extremidade da banda de condução, a extremidade da banda de valência, o nível quase-Fermi dos eletrões e o nível quase-Fermi das lacunas. A sua separação espacial e as suas inclinações revelam a barreira de potencial, a resposta à polarização aplicada e o desvio do equilíbrio.
Em baixas taxas de injeção, as extremidades da banda de condução e de valência curvam-se através da junção, como resultado do campo eletrostático intrínseco. Com o aumento da polarização direta, a barreira efetiva diminui e os níveis quase-Fermi de eletrões e lacunas separam-se à medida que os portadores de carga são impulsionados através do dispositivo. O ponto crucial é que o InP permanece o mesmo semicondutor em ambos os lados da junção: os perfis variáveis refletem a sua eletrostática, definida pela dopagem, e a injeção de portadores.
A comparação entre ?? e ?? ilustra como o desvio terminal altera a estrutura da banda interna. Em particular, os níveis de quasi-Fermi fornecem um indicador direto do estado de não equilíbrio estabelecido durante a condução em regime estacionário.
6.2 Densidades de eletrões e lacunas.
Em seguida, apresente os gráficos Q_nfree.csv e Q_pfree.csv, que mostram a distribuição espacial dos eletrões livres e das lacunas. Em baixas tensões, os eletrões predominam no lado do tipo n e as lacunas no lado do tipo p. As suas densidades diminuem significativamente na região da junção, refletindo a região de depleção.
À medida que a polarização direta aumenta, os portadores minoritários são injetados nos lados opostos. Isso aumenta a sobreposição das populações de eletrões e lacunas e altera a forma de ambos os perfis de densidade. Em níveis de injeção mais elevados, as populações de portadores livres em toda a junção diferem substancialmente dos seus valores em baixa polarização.
Compare a mudança na densidade com a posição da transição de dopagem e, em seguida, com as localizações onde ocorre uma forte curvatura da banda de energia. Isso estabelece uma ligação entre os portadores de carga fornecidos pelo contato e a barreira eletrostática produzida pela junção PN.
6.3 Densidades de corrente de eletrões e lacunas.
Trace os gráficos Jn.csv e Jp.csv. Estes perfis mostram as contribuições dos eletrões e das lacunas para a corrente elétrica que atravessa o dispositivo. Os componentes individuais podem variar em função da posição, uma vez que a recombinação transfere portadores entre as populações de eletrões e lacunas, mas a corrente total em estado estacionário deve permanecer aproximadamente constante quando não houver acumulação de carga líquida.
Em correntes muito baixas, a anulação entre os termos de deriva e difusão pode fazer com que os perfis de corrente individuais pareçam ruidosos ou apresentem uma estrutura abrupta. As capturas de tela obtidas com baixa polarização mostram este comportamento de forma mais evidente do que as capturas obtidas com correntes mais altas. É útil analisar os perfis numa escala adequada e comparar a corrente total com o valor terminal J–V antes de interpretar os detalhes numéricos.
6.4 Perfil de recombinação
A recombinação determina o local onde os eletrões e as lacunas injetadas são removidos das respetivas populações de portadores de carga. Abra o visualizador de imagens com R_nfree_to_pfree.csv e R_ss_srh.csv para comparar os canais de recombinação direta e os canais de recombinação assistida por defeitos, respetivamente. Ambos os processos dependem da disponibilidade de eletrões e lacunas, pelo que as suas distribuições espaciais evoluem à medida que as populações de portadores injetadas se alteram.
No início da condução, a recombinação pode concentrar-se perto da junção, onde as populações de eletrões e lacunas se sobrepõem. Com uma injeção direta mais intensa, a sobreposição estende-se a uma área maior do semicondutor, e os perfis de recombinação podem tornar-se mais amplos. Compare a posição e a largura dos picos de recombinação com as densidades de portadores na mesma polarização.
Ao comparar os dois canais, considere tanto o valor máximo quanto a integral espacial. Um pico local acentuado não precisa dominar a recombinação integrada se outro canal estiver ativo em um volume muito maior. Isso é particularmente importante ao relacionar os perfis de recombinação com a corrente elétrica total.
6.5 Perfis da taxa resultante de geração e recombinação de eletrões e lacunas.
As imagens identificadas como R_n.csv e R_p.csv fornecem informações adicionais sobre as taxas que entram nas equações de continuidade de elétrons e lacunas. Estas podem indicar onde o equilíbrio dos termos locais de geração, captura e recombinação está mais fortemente concentrado. Interprete os seus sinais de acordo com a convenção de representação gráfica utilizada pela OghmaNano.
Considerando todos esses elementos em conjunto, a estrutura de bandas, as densidades de portadores de carga, as correntes e os perfis de taxa fornecem uma explicação completa da característica do diodo: a barreira interna controla a injeção, a tensão aplicada modifica a eletrostática interna e o transporte e a recombinação determinam a corrente de estado estacionário resultante.
7. Exercícios e interpretação
8. Interpretação do modelo e análise complementar.
Este exemplo utiliza um modelo de deriva-difusão de estado estacionário unidimensional de InP, uma estrutura compacta de semicondutor p/n/n⁺ e terminais óhmicos idealizados. É adequado para demonstrar como a retificação elétrica surge da barreira da junção e como esse comportamento está relacionado com a distribuição espacial dos portadores de carga. Estudos mais detalhados do dispositivo podem incluir características de contato medidas, parâmetros de material refinados, variação de temperatura e detalhes geométricos adicionais.
As propriedades do material InP são obtidas a partir do modelo Lua selecionado. Para estudos quantitativos, confirme a largura de banda, a afinidade eletrónica, a mobilidade e os coeficientes de recombinação avaliados no editor de materiais ou nos dados de saída. Verifique a conservação da corrente e os resíduos do solucionador e repita pontos de polarização representativos com uma malha mais fina para estabelecer a estabilidade numérica.
No exemplo fornecido, a variação padrão abrange o intervalo de 0 a 3 V. Um experimento interessante a ser realizado em seguida é incluir uma pequena faixa de polarização reversa, o que permitirá comparar diretamente a condução direta e a corrente de fuga reversa. A faixa de polarização reversa e as características de fuga relevantes do ponto de vista físico dependem do dispositivo e dos mecanismos incluídos na simulação.
Para obter informações mais detalhadas, altere as concentrações de dopagem uma de cada vez. Aumentar a concentração nominal de dopagem geralmente reduz a região de depleção e altera a eletrostática em equilíbrio; modificar os coeficientes de recombinação influencia a corrente dependente da injeção. Estas são comparações úteis, pois os materiais e a geometria da camada podem permanecer fixos.
9. Próximos passos
O díodo PN de InP oferece um ponto de partida útil para estruturas semicondutoras de InP mais avançadas. Uma evolução natural é introduzir uma região com baixa concentração de impurezas e comparar o comportamento eletrostático e de transporte de uma junção PN e uma junção PIN. Outra abordagem é examinar heteroestruturas que envolvem InGaAs/InP, onde as diferenças na composição introduzem descontinuidades de banda adicionais e um novo comportamento na injeção de portadores.