A stochastic simulation of the dislocation-mediated etching of porous GaN distributed Bragg reflectors

Título: Uma simulação estocástica da etching mediada por discordâncias de refletores de Bragg distribuídos (DBRs) de GaN poroso

1. O Problema

Os refletores de Bragg distribuídos (DBRs) de nitreto de gálio (GaN) poroso são estruturas essenciais para dispositivos optoeletrônicos, como LEDs e lasers, pois permitem o controle da extração de luz e o confinamento óptico. A fabricação desses DBRs envolve a porosificação eletroquímica seletiva de camadas dopadas (n-type) em uma estrutura alternada com camadas não intencionalmente dopadas (NID).

O desafio central reside em como o agente etchante (electrolito) acessa as camadas dopadas enterradas, já que as camadas NID não conduzem corrente e, teoricamente, não seriam atacadas.

• Modelo Tradicional ("Kebab"): A interpretação inicial era que cada discordância de rosca (TD - threading dislocation) atuava como um canal vertical contínuo, permitindo que o etchante atravessasse todas as camadas NID e atacasse cada camada dopada sequencialmente, formando uma estrutura semelhante a um "kebab".
• Realidade Experimental: Estudos recentes (usando tomografia FIB e STEM) mostraram que o processo é mais complexo. Nem todas as TDs permanecem ativas em todas as camadas. O caminho de etching pode desviar lateralmente através de uma camada dopada para atingir outra TD, criando estruturas em "cascata". Além disso, a morfologia e a eficiência do etching variam significativamente com a tensão aplicada, e a física exata que controla a ativação/desativação dessas TDs e a formação de poros não é totalmente compreendida.

2. Metodologia

Os autores desenvolveram uma simulação estocástica em Python (utilizando a biblioteca Numba para eficiência) para modelar o processo de etching eletroquímico.

• Modelo Geométrico: A estrutura DBR é representada como uma matriz de pixels (5 nm x 5 nm). O modelo inclui camadas alternadas de GaN dopado e NID, com discordâncias de rosca (TDs) penetrando verticalmente através da estrutura.
• Mecanismo Probabilístico:
  • As pixels de GaN NID têm probabilidade zero de serem etchadas.
  • As pixels de GaN dopado e as TDs têm probabilidades finitas de etching ($P_{doped}$ e $P_{disloc}$, respectivamente).
  • O algoritmo busca pixels vizinhos ao etchante existente e decide, com base nas probabilidades atribuídas, se o pixel é removido (convertido em etchante).
  • Uma regra adicional impede a fusão de poros vizinhos, mantendo domínios de porosidade não interconectados, conforme observado experimentalmente.
• Correlação com Dados Experimentais:
  • Foram sintetizados três tipos de amostras (A, B e C) via MOVPE. A Amostra A é um DBR padrão; as Amostras B e C possuem espessuras de camadas dopadas variáveis (crescentes e decrescentes, respectivamente).
  • Realizou-se etching eletroquímico em diferentes tensões (6V a 12V) e mediu-se a cronoamperometria (corrente vs. tempo).
  • A simulação foi executada variando $P_{doped}$ para simular os efeitos de diferentes tensões aplicadas, mantendo $P_{disloc}$ constante.
  • Os resultados da simulação (curvas de corrente e morfologia) foram comparados com dados experimentais de microscopia eletrônica de varredura (SEM) e tomografia FIB.

3. Contribuições Principais

1. Desenvolvimento de um Modelo Estocástico: Criação de uma ferramenta computacional capaz de reproduzir a evolução temporal do etching e a formação de microestruturas complexas em DBRs de GaN.
2. Validação do Modelo "Cascata": A simulação demonstra que, ao variar a probabilidade de etching, é possível reproduzir a transição entre o modelo "kebab" (TDs ativas em todas as camadas) e o modelo "cascata" (TDs que se ativam e desativam ao longo da profundidade, com propagação lateral).
3. Correlação Tensão-Probabilidade: Estabelecimento de uma correlação qualitativa onde o aumento da tensão de etching experimental corresponde ao aumento da probabilidade de etching das camadas dopadas ($P_{doped}$) na simulação.
4. Generalização da Aplicabilidade: Demonstração de que o modelo pode prever comportamentos de cronoamperometria não apenas para DBRs uniformes, mas também para estruturas com espessuras de camadas variáveis e até para geometrias com trincheiras laterais (onde o etchante acessa as camadas pelos lados).

4. Resultados

• Morfologia e "Kebabs" vs. "Cascades": A simulação reproduz com sucesso a observação experimental de que, em tensões mais baixas, há uma maior proporção de estruturas em cascata e menor eficiência de etching em camadas profundas. Em tensões mais altas (simuladas por $P_{doped}$ mais alto), a proporção de estruturas "kebab" (K5) aumenta, alinhando-se com dados de tomografia FIB.
• Cronoamperometria:
  • As curvas simuladas reproduzem as oscilações de corrente observadas experimentalmente, que correspondem à etching camada por camada.
  • A simulação explica o decaimento das oscilações de corrente: ocorre porque o frente de etching não é uniforme; algumas TDs avançam mais rápido que outras, fazendo com que camadas mais profundas comecem a ser etchadas antes que as camadas superiores sejam completamente removidas, superpondo as correntes e reduzindo a amplitude dos picos.
  • A variação de $P_{doped}$ na simulação consegue mimetizar as mudanças na forma das curvas de corrente ao alterar a tensão experimental (de 6V a 12V).
• Limitações Identificadas: O modelo atual não consegue reproduzir a formação de poros conicos nas camadas NID (observados experimentalmente em altas tensões) nem a variação no tamanho dos poros (pequenos em baixa tensão, grandes em alta tensão), pois o modelo assume probabilidade zero de etching em NID fora da linha da TD.
• Aplicação em Geometrias Alternativas: A simulação adaptada para geometrias com trincheiras laterais sugere que, mesmo quando o etchante entra pelos lados, as TDs podem ser ativadas a certa distância da trincheira, alterando a morfologia dos poros e possivelmente afetando propriedades como a birrefringência.

5. Significância

Este trabalho é significativo porque oferece uma ferramenta de modelagem robusta para entender um processo de fabricação complexo e não totalmente compreendido na tecnologia de semicondutores de nitreto.

• Insights Físicos: Sugere que a tensão aplicada afeta diferentemente a etching da camada dopada e a ativação da TD, indicando mecanismos físicos distintos (possivelmente relacionados ao tunelamento Zener e dopantes não intencionais).
• Otimização de Processos: A capacidade de prever curvas de cronoamperometria e morfologias de poros sem a necessidade de testes experimentais extensivos permite otimizar parâmetros de fabricação para obter DBRs com propriedades ópticas específicas.
• Versatilidade: A demonstração de que o modelo se aplica a estruturas não uniformes e geometrias com trincheiras amplia sua utilidade para o desenvolvimento de uma gama mais ampla de dispositivos fotônicos e eletrônicos baseados em GaN poroso.

Em resumo, a simulação estocástica proposta preenche uma lacuna entre a observação microestrutural e a resposta macroscópica (corrente), fornecendo um quadro teórico para a engenharia de DBRs de GaN poroso mediado por defeitos.