Comparison of Silvaco and Synopsys TCAD Predictions Including the Perugia Radiation Damage Model in Silicon Pixel Detectors for the HL-LHC

Este artigo compara as simulações TCAD do Silvaco e do Synopsys que incorporam o modelo de danos por radiação de Perugia para avaliar sua precisão preditiva em métricas de desempenho fundamentais de detectores de pixels de silício sob as condições extremas de radiação esperadas no High Luminosity LHC.

O essencial

Imagine o Grande Colisor de Hádrons (LHC) como uma pista de corrida de partículas massiva e de alta velocidade. Os cientistas estão atualizando-o para a versão "Alta Luminosidade" (HL-LHC), o que significa que eles estarão colidindo partículas com muito mais frequência. O problema? Esse tráfego intenso cria muita "poeira de radiação" que danifica os sensores de silício (as câmeras) que tentam capturar imagens das colisões.

Com o tempo, essa poeira de radiação transforma os sensores de silício em dispositivos "vazados" e "rígidos". Eles começam a perder a capacidade de coletar sinais (como uma câmera perdendo o foco) e exigem tensões muito mais altas para funcionar, o que corre o risco de quebrá-los.

Para corrigir isso antes que aconteça, os cientistas usam simulações computacionais para prever como os sensores se comportarão após anos de radiação. Eles precisam saber: Quanta tensão precisamos? Quanto corrente vazará? O sensor ainda funcionará?

Os Dois "Previsionistas do Tempo"

Neste artigo, os pesquisadores estão testando dois programas de computador diferentes (ferramentas TCAD) que atuam como previsionistas do tempo para esses sensores:

1. Synopsys
2. Silvaco

Ambos os programas usam um conjunto específico de regras chamado "Modelo de Danos por Radiação de Perugia". Pense nesse modelo como um manual de instruções detalhado que diz ao computador exatamente como a "poeira de radiação" danifica o silício, criando pequenas armadilhas e buracos que atrapalham o fluxo elétrico.

O objetivo deste artigo é ver se esses dois " previsionistas" diferentes dão a mesma previsão quando usam o mesmo manual de instruções. Se eles concordarem, significa que o manual é confiável, e os cientistas podem confiar nas previsões, não importa qual software usem.

O Experimento: Um Diodo de Silício Minúsculo

Os pesquisadores construíram um modelo virtual 2D de um sensor de silício minúsculo (um diodo) com 50 micrômetros de espessura (cerca da largura de um fio de cabelo humano). Eles simularam dois cenários:

1. Sensor Novo: Antes de qualquer radiação atingi-lo.
2. Sensor Irradiado: Após ser atingido por uma quantidade massiva de radiação (simulando o ambiente hostil do HL-LHC).

Eles testaram esses sensores em duas temperaturas: uma fria de 248 K (cerca de -25°C) e uma quente de 300 K (temperatura ambiente).

Os Resultados: Os Previsionistas Concordam?

1. O Sensor Novo (Não Irradiado)
Quando o sensor estava novo, ambos os programas de computador concordaram quase perfeitamente sobre quanto de eletricidade fluía através dele e como ele armazenava carga, até cerca de 500 volts.

• A Discrepância: Quando empurraram a tensão muito alta (perto de 700 volts), os programas começaram a discordar ligeiramente sobre exatamente quando o sensor iria "quebrar" (ruptura). Os autores sugerem que isso provavelmente ocorre porque os dois programas usam "grades" digitais (malhas) ligeiramente diferentes para desenhar o sensor, semelhante a como dois aplicativos de mapa diferentes podem desenhar uma estrada ligeiramente diferente.

2. O Sensor Irradiado (O Teste Real)
É aqui que a verdadeira mágica aconteceu. Eles simularam o sensor após ter sido bombardeado por radiação.

• Corrente de Fuga: Ambos os programas previram o "vazamento" (eletricidade indesejada) quase identicamente.
• Tensão de Depleção: Ambos concordaram perfeitamente sobre quanta tensão era necessária para fazer o sensor funcionar novamente.
• Campos Elétricos: Eles mapearam as forças elétricas invisíveis dentro do silício. No meio do sensor (o "bulk"), os dois programas coincidiram quase perfeitamente (dentro de 1% um do outro).
• As "Armadilhas": Eles também observaram as pequenas "armadilhas" criadas pela radiação que capturam elétrons. Os dois programas concordaram sobre o comportamento dessas armadilhas dentro de uma margem muito razoável (cerca de 20%).

A Reviravolta da Temperatura:
Na temperatura ambiente (300 K), os programas discordaram um pouco mais nos níveis mais altos de radiação. No entanto, os autores observam que isso não é uma grande preocupação porque, no mundo real, esses sensores danificados quase nunca são operados na temperatura ambiente; eles são mantidos muito frios para sobreviver. Portanto, o acordo na temperatura fria (248 K) é o que realmente importa, e lá, os dois programas estavam em perfeita sincronia.

A Conclusão

O artigo conclui que Synopsys e Silvaco são como dois chefs diferentes seguindo exatamente a mesma receita (o Modelo de Perugia) e terminando com o mesmo prato delicioso.

Embora as ferramentas de software sejam diferentes, quando usam o modelo de danos por radiação de Perugia, elas produzem previsões quase idênticas sobre como os sensores de silício sobreviverão à radiação hostil do futuro HL-LHC. Isso dá aos cientistas confiança de que seus modelos são sólidos e de que podem usar qualquer uma das ferramentas para projetar a próxima geração de detectores de partículas.

Nota: Os autores mencionam que planejam examinar a "carga coletada" no futuro, mas este artigo focou estritamente em tensão, corrente e campos elétricos.

Resumo técnico

1. Declaração do Problema

O Grande Colisor de Hádrons de Alta Luminosidade (HL-LHC) submeterá detectores de pixels de silício a fluências de radiação 5 a 10 vezes superiores aos níveis atuais. Essa radiação intensa causa danos no volume e na superfície dos sensores de silício, levando a:

• Aumento da corrente de fuga.
• Deslocamentos na tensão operacional (depleção).
• Perda de sinal e degradação da coleta de carga.
• Modificações no comportamento da eletrônica de front-end.

A previsão precisa desses efeitos é crítica para estimar tensões operacionais e treinar algoritmos de reconstrução. Embora o Modelo de Danos por Radiação da Universidade de Perugia (incorporando física específica de defeitos volumétricos e superficiais) tenha sido implementado com sucesso no Synopsys Sentaurus TCAD, havia a necessidade de validar sua portabilidade e poder preditivo no Silvaco TCAD. Sem a validação cruzada entre essas duas ferramentas de simulação padrão da indústria, a robustez da modelagem de danos por radiação permanecia não verificada para adoção mais ampla.

2. Metodologia

Os autores realizaram um estudo comparativo de simulações de dispositivos utilizando ambas as ferramentas TCAD, Synopsys e Silvaco.

• Estrutura do Dispositivo: Uma simulação 2D de um diodo n-on-p foi utilizada.
  • Espessura: 50 µm.
  • Concentração no Volume: $5 \times 10^{12} \text{ cm}^{-3}$ (Boro).
  • Passo (Pitch): 20 µm com um pequeno eletrodo de coleta.
  • Superfície: Camada de óxido de 1 µm.
• Modelos Físicos: Ambas as ferramentas utilizaram modelos físicos idênticos para garantir uma comparação justa:
  • Estatísticas: Fermi-Dirac.
  • Recombinação: Scharfetter - Shockley Read Hall.
  • Tunelamento banda-a-banda: Schenk.
  • Estreitamento da banda proibida: Slotboom.
  • Ionização por impacto: van Overstraeten e de Man.
• Implementação de Danos por Radiação: O modelo de Perugia foi portado para o Silvaco. Isso inclui:
  • Modelo Volumétrico: Dois defeitos aceitadores e um defeito doador com densidades proporcionais à fluência ($\Phi$). Utiliza uma função dependente da fluência para a criação de aceitadores ($N_A$) para lidar com altas fluências ($> 3 \times 10^{15} \text{ neq/cm}^2$).
  • Modelo de Superfície: Uma distribuição contínua de energia de defeitos (semelhantes a aceitadores na banda proibida superior, semelhantes a doadores próximos à banda de valência) e aumento linear da carga do óxido com a dose.
• Cenários de Simulação:
  • Fluências: Não irradiado, $\Phi = 1.5 \times 10^{15} \text{ neq/cm}^2$, e $\Phi = 1.0 \times 10^{16} \text{ neq/cm}^2$.
  • Temperaturas: $T = 248 \text{ K}$ e $T = 300 \text{ K}$.
• Observáveis: Corrente-Tensão (IV), Capacitância-Tensão (CV), perfis de campo elétrico e probabilidades de ionização de armadilhas.

3. Principais Contribuições

• Portabilidade do Modelo: Portagem bem-sucedida dos complexos modelos de danos por radiação volumétricos e superficiais de Perugia do Synopsys para o Silvaco TCAD.
• Validação Cruzada: Fornecimento da primeira comparação direta entre essas duas engines TCAD distintas ao simular ambientes de alta radiação relevantes para o HL-LHC.
• Análise de Estatísticas de Defeitos: Análise não apenas das características elétricas macroscópicas, mas também das probabilidades de ionização de armadilhas microscópicas, oferecendo uma visão mais profunda da consistência da implementação da física de defeitos.

4. Resultados

Dispositivo Não Irradiado

• Características IV: Excelente concordância entre as ferramentas até ~500 V. Discrepâncias apareceram em tensões mais altas (o Synopsys previu ruptura em ~700 V versus ~760 V no Silvaco), provavelmente devido a diferenças na geração de malha.
• Características CV: Ambas as ferramentas previram com precisão uma tensão de depleção de ~10 V e corresponderam à característica de "degrau" em ~55 V (atribuída à depleção sob o óxido).

Dispositivo Irradiado

• Corrente de Fuga (IV):
  • A 248 K: Concordância notavelmente boa em todas as fluências. Desvios menores ocorreram apenas em tensões muito altas e na fluência mais elevada, atribuídos a diferenças de malha.
  • A 300 K: A concordância foi apenas qualitativa na fluência mais alta. No entanto, os autores observam que isso é menos crítico, pois os dispositivos raramente são operados acima de 0°C após tal alta irradiação.
• Tensão de Depleção (CV): A 248 K, a concordância entre as duas ferramentas foi "praticamente perfeita" para ambas as fluências.
• Perfis de Campo Elétrico:
  • Na região volumétrica, os perfis de campo elétrico concordaram dentro de 1% para ambas as fluências e todas as tensões de polarização (100 V, 300 V, 500 V).
  • Pequenas discrepâncias foram observadas perto das implantações frontal e traseira, novamente atribuídas a incompatibilidades de malha e diferenças na modelagem de implantes.
• Probabilidades de Ionização de Armadilhas:
  • Para defeitos volumétricos (armadilhas aceitadoras 1 e 2), as probabilidades de ionização concordaram dentro de uma ordem de magnitude, e dentro de 20% na região volumétrica.
  • Este nível de concordância é considerado satisfatório, dado o acordo em nível percentual nos campos elétricos.

5. Significado

• Validação do Modelo de Perugia: O estudo confirma que o modelo de danos por radiação da Universidade de Perugia é robusto e independente de plataforma, produzindo resultados consistentes em diferentes ambientes de software TCAD.
• Confiabilidade para o HL-LHC: O alto grau de concordância em observáveis críticos (corrente de fuga, tensão de depleção e campo elétrico) fornece confiança no uso dessas simulações para prever o desempenho dos sensores sob condições do HL-LHC.
• Orientação Operacional: Os resultados apoiam o uso desses modelos para estimar tensões operacionais seguras e para entender mecanismos de perda de sinal, o que é essencial para o projeto e algoritmos de reconstrução dos rastreadores ATLAS e CMS atualizados.
• Trabalho Futuro: Os autores planejam estender essa validação para simulações de carga coletada, que é a métrica definitiva para o desempenho de rastreamento.

Conclusão: O artigo demonstra que as ferramentas TCAD Silvaco e Synopsys podem produzir previsões mutuamente consistentes para detectores de pixels de silício irradiados ao utilizar o modelo de danos por radiação de Perugia, validando a solidez do modelo para futuras aplicações em colisores de alta luminosidade.