Electric-field effects on defect migration energetics in GaN

Resumo Técnico: Efeitos de campo elétrico na energética de migração de defeitos em GaN

Definição do Problema
A confiabilidade de dispositivos baseados em Nitreto de Gálio (GaN) e AlGaN, particularmente em ambientes de alta potência e radiação, é comprometida pela migração e reconfiguração de defeitos pontuais. Embora a formação de defeitos sob irradiação seja bem estudada, o comportamento de defeitos existentes sob a influência combinada de campos elétricos operacionais e temperatura permanece menos compreendido. Em transistores de alta mobilidade eletrônica (HEMTs) e dispositivos de potência, os campos operacionais frequentemente se aproximam do limite de ruptura (~3 MV/cm), podendo impulsionar o transporte de defeitos que leva a instabilidades de tensão de limiar e deriva de desempenho. As abordagens computacionais existentes enfrentam um compromisso: a Teoria do Funcional da Densidade (DFT) oferece alta precisão, mas é restrita a supercélulas pequenas, enquanto a Dinâmica Molecular (MD) clássica com potenciais interatômicos padrão carece da capacidade de modelar a transferência de carga dinâmica e o acoplamento eletrostático necessários para simular a migração de defeitos sob campos eléticos externos.

Metodologia
Para abordar essas limitações, os autores desenvolveram um campo de força reativo ReaxFF especificamente para GaN. A metodologia envolveu três estágios principais:

1. Desenvolvimento do Campo de Força: Um conjunto de dados de treinamento abrangente foi gerado usando cálculos DFT (VASP) dentro do framework do funcional PBE. Este conjunto de dados incluiu dados de equação de estado (EOS) para várias fases de Ga e GaN (wurtzita, sal-gema, zinco, etc.), perfis moleculares de N2 e energias totais de supercélulas contendo defeitos pontuais neutros (vacâncias, antissítios, pares de Frenkel, intersticiais e defeitos de Schottky).
2. Validação: Os parâmetros do ReaxFF foram otimizados contra este conjunto de dados DFT. O campo de força resultante foi validado comparando parâmetros de rede, energias de coesão, módulo de bulk e energias de formação/migração de defeitos contra resultados experimentais e de DFT prévios.
3. Simulações de Campo Elétrico: Usando o modelo ReaxFF validado no LAMMPS, os autores realizaram cálculos de Nudged Elastic Band (CI-NEB) para determinar as barreiras de migração para defeitos pontuais intrínsecos ($V_{Ga}$, $V_N$, $Ga_{i1,oct}$ e $N_{i2,split}$) sob campos elétricos externos variando de 0 a 0,04 V/Å (0–40 MV/cm). As simulações levaram em conta o equilíbrio de carga dinâmica (QEq) e distorções de rede locais.

Principais Contribuições e Resultados

• Validação do Campo de Força: O modelo ReaxFF desenvolvido reproduz com sucesso as propriedades estruturais de DFT (ex: constantes de rede dentro de 0,3% e módulo de bulk dentro de 1,5%) e a energética de defeitos. Embora existam algumas divergências nas energias absolutas de formação de defeitos (até 26% para certas configurações), o modelo captura as tendências relativas e as barreiras de migração com alta fidelidade em relação à literatura.
• Resposta de Campo Não Linear: Contrariando a simples aproximação linear de dipolo-campo ($\Delta E_B = -\mu_{eff} \cdot E$), o estudo revela que os campos elétricos modificam as barreiras de migração de forma anisotrópica e não linear. A resposta depende fortemente do tipo específico de defeito e de seu mecanismo de migração.
  • Vacâncias: A migração de vacâncias de Nitrogênio ($V_N$) exibe uma resposta quase linear e simétrica ao campo, consistente com um modelo de dipolo rígido devido à mínima deformação de rede e variação de carga. As vacâncias de Gálio ($V_{Ga}$) mostram uma assimetria leve; a barreira de avanço aumenta enquanto a barreira de recuo diminui, mas com inclinações desiguais, atribuído ao maior tamanho atômico do Ga induzindo deformação local e redistribuição de carga.
  • Intersticiais: Defeitos intersticiais ($Ga_{i1,oct}$ e $N_{i2,split}$) exibem pronunciada não linearidade e assimetria. Para $Ga_{i1,oct}$, o campo é forte o suficiente para deslocar a configuração do ponto de sela limitante da taxa, causando um cruzamento de qual pico de energia domina a barreira. Esses defeitos envolvem o movimento cooperativo de múltiplos átomos, levando a uma complexa redistribuição de carga parcial e respostas de deformação que divergem significamente das previsões lineares.
• Mecanismo de Desvio: A partida do comportamento linear surge da redistribuição de carga parcial induzida pelo campo e da distorção da rede local ao longo da coordenada de reação. O momento de dipolo efetivo não é fixo, mas evolui continuamente conforme o defeito migra, um fenômeno capturado nativamente pelo campo de força reativo, mas perdido por modelos de dipolo fixo.

Significância e Alegações
O artigo afirma que os campos elétricos não apenas enviesam o movimento de defeitos linearmente, mas podem modificar complexamente o panorama de migração de defeitos através do acoplamento carga-rede.

• Impacto Cinético: O estudo destaca que, embora as reduções de barreira induzidas pelo campo para vacâncias de alta barreira (> 2 eV) possam não alterar significativamente as taxas de difusão em escalas de tempo de dispositivos, os efeitos são cineticamente significativos para defeitos intersticiais com barreiras próximas a 1 eV. Uma redução de ~0,05 eV pode aumentar as taxas de migração em uma ordem de magnitude à temperatura ambiente, potencialmente acelerando a redistribuição de defeitos em regiões de alto campo de dispositivos operacionais.
• Confiabilidade do Dispositivo: Essas descobertas sugerem que, em regiões de alto campo (ex: bordas de gate em HEMTs), os defeitos intersticiais podem tornar-se substancialmente mais móveis, contribuindo para a geração de armadilhas (traps) e deriva da tensão de limiar. A assimetria direcional da resposta implica um desloco líquido de defeitos em direção a regiões específicas do dispositivo, potencialmente impulsionando problemas de confiabilidade a longo prazo.
• Escopo Metodológico: Os autores observam que seus resultados correspondem a defeitos neutros devido à natureza clássica do esquema de equilíbrio de carga QEq. Eles reconhecem que defeitos altamente carregados em condições reais de operação (ex: $V_{Ga}^{3-}$) podem exibir respostas de campo ainda mais fortes, sugerindo que os achados atuais representam um limite inferior para os efeitos de campo. No entanto, o trabalho estabelece que o acoplamento campo-defeito em GaN é anisotrópico e específico para cada defeito, necessitando de simulações reativas para previsões precisas sob condições operacionais realistas.