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

Este trabajo compara las simulaciones TCAD de Silvaco y Synopsys que incorporan el modelo de daño por radiación de Perugia para evaluar su precisión predictiva en métricas clave de rendimiento de detectores de píxeles de silicio bajo las condiciones extremas de radiación esperadas en el LHC de Alta Luminosidad.

Lo esencial

Imagina el Gran Colisionador de Hadrones (LHC) como una pista de carreras de partículas masiva y de alta velocidad. Los científicos lo están actualizando a la versión "de Alta Luminosidad" (HL-LHC), lo que significa que harán chocar partículas entre sí con mucha más frecuencia. ¿El problema? Este tráfico intenso genera mucho "polvo radiactivo" que daña los sensores de silicio (las cámaras) que intentan tomar fotografías de las colisiones.

Con el tiempo, este polvo radiactivo convierte los sensores de silicio en dispositivos "con fugas" y "rígidos". Comienzan a perder su capacidad de recolectar señales (como una cámara que pierde el enfoque) y requieren un voltaje mucho más alto para funcionar, lo que corre el riesgo de romperlos.

Para solucionar esto antes de que ocurra, los científicos utilizan simulaciones por computadora para predecir cómo se comportarán los sensores después de años de radiación. Necesitan saber: ¿Cuánto voltaje necesitamos? ¿Cuánta corriente se filtrará? ¿El sensor seguirá funcionando?

Los dos "pronosticadores del tiempo"

En este artículo, los investigadores están probando dos programas informáticos diferentes (herramientas TCAD) que actúan como pronosticadores del tiempo para estos sensores:

1. Synopsys
2. Silvaco

Ambos programas utilizan un conjunto específico de reglas llamado el "Modelo de Daño por Radiación de Perugia". Piensa en este modelo como un manual de instrucciones detallado que le dice a la computadora exactamente cómo el "polvo radiactivo" daña el silicio, creando pequeñas trampas y huecos que alteran el flujo eléctrico.

El objetivo de este artículo es ver si estos dos "pronosticadores" diferentes dan la misma predicción cuando utilizan el mismo manual de instrucciones. Si están de acuerdo, significa que el manual es fiable y los científicos pueden confiar en las predicciones sin importar qué software utilicen.

El Experimento: Un pequeño diodo de silicio

Los investigadores construyeron un modelo virtual 2D de un pequeño sensor de silicio (un diodo) que tiene 50 micrómetros de grosor (aproximadamente el ancho de un cabello humano). Simularon dos escenarios:

1. Sensor nuevo: Antes de que la radiación lo golpee.
2. Sensor irradiado: Después de ser golpeado por una cantidad masiva de radiación (simulando el entorno hostil del HL-LHC).

Probaron estos sensores a dos temperaturas: una fresca de 248 K (aproximadamente -25°C) y una cálida de 300 K (temperatura ambiente).

Los Resultados: ¿Están de acuerdo los pronosticadores?

1. El Sensor Nuevo (No irradiado)
Cuando el sensor estaba totalmente nuevo, ambos programas informáticos estuvieron de acuerdo casi perfectamente sobre cuánta electricidad fluía a través de él y cómo almacenaba la carga, hasta aproximadamente 500 voltios.

• La Discrepancia: Cuando empujaron el voltaje muy alto (cerca de 700 voltios), los programas comenzaron a discrepar ligeramente sobre exactamente cuándo el sensor se "rompería" (ruptura). Los autores sugieren que esto se debe probablemente a que los dos programas utilizan "cuadrículas" (mallas) digitales ligeramente diferentes para dibujar el sensor, similar a cómo dos aplicaciones de mapas diferentes podrían dibujar una carretera de manera ligeramente distinta.

2. El Sensor Irradiado (La prueba real)
Aquí es donde ocurrió la verdadera magia. Simularon el sensor después de haber sido bombardeado con radiación.

• Corriente de Fuga: Ambos programas predijeron la "fuga" (electricidad no deseada) casi idénticamente.
• Voltaje de Agotamiento: Ambos estuvieron de acuerdo perfectamente sobre cuántos voltios se necesitaban para hacer que el sensor funcionara nuevamente.
• Campos Eléctricos: Mapearon las fuerzas eléctricas invisibles dentro del silicio. En el medio del sensor (el "volumen"), los dos programas coincidieron casi perfectamente (dentro del 1% entre sí).
• Las "Trampas": También observaron las pequeñas "trampas" creadas por la radiación que atrapan electrones. Los dos programas estuvieron de acuerdo sobre el comportamiento de estas trampas dentro de un margen muy razonable (aproximadamente el 20%).

El Giro de la Temperatura:
A temperatura ambiente (300 K), los programas discreparon un poco más a los niveles de radiación más altos. Sin embargo, los autores señalan que esto no es una gran preocupación porque, en el mundo real, estos sensores dañados casi nunca se operan a temperatura ambiente; se mantienen muy fríos para sobrevivir. Por lo tanto, el acuerdo a la temperatura fría (248 K) es lo que realmente importa, y allí los dos programas estuvieron perfectamente sincronizados.

La Conclusión

El artículo concluye que Synopsys y Silvaco son como dos chefs diferentes siguiendo exactamente la misma receta (el Modelo de Perugia) y terminando con el mismo plato delicioso.

Aunque las herramientas de software son diferentes, cuando utilizan el modelo de daño por radiación de Perugia, producen predicciones casi idénticas sobre cómo los sensores de silicio sobrevivirán a la radiación hostil del futuro HL-LHC. Esto da a los científicos confianza de que sus modelos son sólidos y de que pueden utilizar cualquiera de las dos herramientas para diseñar la próxima generación de detectores de partículas.

Nota: Los autores mencionan que planean examinar la "carga recolectada" en el futuro, pero este artículo se centró estrictamente en voltaje, corriente y campos eléctricos.

Resumen técnico

Aquí se presenta un resumen técnico detallado del artículo "Comparación de las predicciones de Silvaco y Synopsys TCAD incluyendo el Modelo de Daño por Radiación de Perugia en Detectores de Píxeles de Silicio para el HL-LHC".

1. Planteamiento del Problema

El Gran Colisionador de Hadrones de Alta Luminosidad (HL-LHC) someterá a los detectores de píxeles de silicio a flujos de radiación de 5 a 10 veces superiores a los niveles actuales. Esta intensa radiación causa daños volumétricos y superficiales en los sensores de silicio, lo que conduce a:

• Aumento de la corriente de fuga.
• Desplazamientos en el voltaje operativo (de agotamiento).
• Pérdida de señal y degradación de la recolección de carga.
• Modificaciones en el comportamiento de la electrónica de entrada.

La predicción precisa de estos efectos es crítica para estimar los voltajes operativos y entrenar algoritmos de reconstrucción. Si bien el Modelo de Daño por Radiación de la Universidad de Perugia (que incorpora física específica de defectos volumétricos y superficiales) se ha implementado con éxito en Synopsys Sentaurus TCAD, existía la necesidad de validar su portabilidad y poder predictivo en Silvaco TCAD. Sin una validación cruzada entre estas dos herramientas de simulación estándar en la industria, la robustez de la modelización del daño por radiación permanece sin verificar para una adopción más amplia.

2. Metodología

Los autores realizaron un estudio comparativo de simulaciones de dispositivos utilizando ambas herramientas TCAD, Synopsys y Silvaco.

• Estructura del Dispositivo: Se utilizó una simulación 2D de un diodo n-on-p.
  • Espesor: 50 µm.
  • Concentración de volumen: $5 \times 10^{12} \text{ cm}^{-3}$ (Boro).
  • Paso: 20 µm con un pequeño electrodo de recolección.
  • Superficie: Capa de óxido de 1 µm.
• Modelos Físicos: Ambas herramientas utilizaron modelos físicos idénticos para garantizar una comparación justa:
  • Estadísticas: Fermi-Dirac.
  • Recombinación: Scharfetter - Shockley Read Hall.
  • Túnel banda a banda: Schenk.
  • Reducción de la banda prohibida: Slotboom.
  • Ionización por impacto: van Overstraeten y de Man.
• Implementación del Daño por Radiación: El modelo Perugia se portó a Silvaco. Esto incluye:
  • Modelo Volumétrico: Dos defectos aceptores y un defecto donador con densidades proporcionales al flujo ($\Phi$). Utiliza una función dependiente del flujo para la creación de aceptores ($N_A$) para manejar flujos altos ($> 3 \times 10^{15} \text{ neq/cm}^2$).
  • Modelo Superficial: Una distribución continua de energía de defectos (tipo aceptor en la parte superior de la banda prohibida, tipo donador cerca de la banda de valencia) y aumento lineal de la carga del óxido con la dosis.
• Escenarios de Simulación:
  • Flujos: Sin irradiar, $\Phi = 1.5 \times 10^{15} \text{ neq/cm}^2$, y $\Phi = 1.0 \times 10^{16} \text{ neq/cm}^2$.
  • Temperaturas: $T = 248 \text{ K}$ y $T = 300 \text{ K}$.
• Observables: Corriente-Voltaje (IV), Capacitancia-Voltaje (CV), perfiles de campo eléctrico y probabilidades de ionización de trampas.

3. Contribuciones Clave

• Portabilidad del Modelo: Portado con éxito el complejo modelo de daño por radiación volumétrico y superficial de Perugia desde Synopsys a Silvaco TCAD.
• Validación Cruzada: Proporcionó la primera comparación directa de estos dos motores TCAD distintos al simular entornos de alta radiación relevantes para el HL-LHC.
• Análisis de Estadísticas de Defectos: Analizó no solo las características eléctricas macroscópicas, sino también las probabilidades de ionización de trampas microscópicas, ofreciendo una visión más profunda de la consistencia en la implementación de la física de defectos.

4. Resultados

Dispositivo Sin Irradiar

• Características IV: Excelente acuerdo entre las herramientas hasta ~500 V. Aparecieron discrepancias a voltajes más altos (Synopsys predijo ruptura a ~700 V frente a Silvaco a ~760 V), probablemente debido a diferencias en la generación de mallas.
• Características CV: Ambas herramientas predijeron con precisión un voltaje de agotamiento de ~10 V y coincidieron con la característica de "escalón" a ~55 V (atribuida al agotamiento bajo el óxido).

Dispositivo Irradiado

• Corriente de Fuga (IV):
  • A 248 K: Acuerdo notablemente bueno en todos los flujos. Se produjeron desviaciones menores solo a voltajes muy altos y al flujo más alto, atribuidas a diferencias en la malla.
  • A 300 K: El acuerdo fue cualitativo solo en el flujo más alto. Sin embargo, los autores señalan que esto es menos crítico, ya que los dispositivos rara vez se operan por encima de 0°C después de una irradiación tan alta.
• Voltaje de Agotamiento (CV): A 248 K, el acuerdo entre las dos herramientas fue "prácticamente perfecto" para ambos flujos.
• Perfiles de Campo Eléctrico:
  • En la región volumétrica, los perfiles de campo eléctrico coincidieron dentro del 1% para ambos flujos y todos los voltajes de polarización (100 V, 300 V, 500 V).
  • Se observaron ligeras discrepancias cerca de las implantaciones frontal y trasera, nuevamente atribuidas a desajustes de malla y diferencias en la modelización de implantaciones.
• Probabilidades de Ionización de Trampas:
  • Para defectos volumétricos (trampas aceptoras 1 y 2), las probabilidades de ionización coincidieron dentro de un orden de magnitud, y dentro del 20% en la región volumétrica.
  • Este nivel de acuerdo se considera satisfactorio dado el acuerdo a nivel porcentual en los campos eléctricos.

5. Significado

• Validación del Modelo de Perugia: El estudio confirma que el modelo de daño por radiación de la Universidad de Perugia es robusto e independiente de la plataforma, produciendo resultados consistentes en diferentes entornos de software TCAD.
• Fiabilidad para el HL-LHC: El alto grado de acuerdo en observables críticos (corriente de fuga, voltaje de agotamiento y campo eléctrico) proporciona confianza en el uso de estas simulaciones para predecir el rendimiento de los sensores bajo condiciones del HL-LHC.
• Guía Operativa: Los resultados apoyan el uso de estos modelos para estimar voltajes operativos seguros y comprender los mecanismos de pérdida de señal, lo cual es esencial para el diseño y los algoritmos de reconstrucción de los rastreadores actualizados de ATLAS y CMS.
• Trabajo Futuro: Los autores planean extender esta validación a simulaciones de carga recolectada, que es la métrica definitiva para el rendimiento de rastreo.

Conclusión: El artículo demuestra que las herramientas TCAD de Silvaco y Synopsys pueden producir predicciones mutuamente consistentes para detectores de píxeles de silicio irradiados cuando se utiliza el modelo de daño por radiación de Perugia, validando la solidez del modelo para futuras aplicaciones de colisionadores de alta luminosidad.