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

Ce papier compare les simulations TCAD de Silvaco et de Synopsys intégrant le modèle de dommages par rayonnement Perugia afin d'évaluer leur précision prédictive pour les métriques de performance clés des détecteurs à pixels en silicium dans les conditions de rayonnement extrêmes attendues au LHC à haute luminosité.

L'essentiel

Imaginez le Grand collisionneur de hadrons (LHC) comme une piste de course géante et ultra-rapide pour les particules. Les scientifiques le mettent à niveau vers la version « Haute Luminosité » (HL-LHC), ce qui signifie qu'ils vont faire entrer en collision des particules beaucoup plus fréquemment. Le problème ? Ce trafic intense génère beaucoup de « poussière de radiation » qui endommage les capteurs en silicium (les caméras) tentant de photographier les collisions.

Au fil du temps, cette poussière de radiation transforme les capteurs en silicium en dispositifs « fuyards » et « rigides ». Ils commencent à perdre leur capacité à collecter des signaux (comme une caméra qui perd le focus) et nécessitent une tension beaucoup plus élevée pour fonctionner, ce qui risque de les détruire.

Pour corriger cela avant qu'il ne se produise, les scientifiques utilisent des simulations informatiques pour prédire le comportement des capteurs après des années de radiation. Ils doivent savoir : Quelle tension est nécessaire ? Quel courant va fuir ? Le capteur fonctionnera-t-il encore ?

Les deux « prévisionnistes météo »

Dans cet article, les chercheurs testent deux programmes informatiques différents (outils TCAD) qui agissent comme des prévisionnistes météo pour ces capteurs :

1. Synopsys
2. Silvaco

Les deux programmes utilisent un ensemble spécifique de règles appelé le « Modèle de dommages par radiation Perugia ». Imaginez ce modèle comme un manuel d'instructions détaillé qui indique à l'ordinateur exactement comment la « poussière de radiation » endommage le silicium, créant de minuscules pièges et des trous qui perturbent le flux électrique.

L'objectif de cet article est de voir si ces deux « prévisionnistes » différents donnent la même prédiction lorsqu'ils utilisent le même manuel d'instructions. S'ils s'accordent, cela signifie que le manuel est fiable, et les scientifiques peuvent faire confiance aux prédictions, peu importe le logiciel utilisé.

L'expérience : une minuscule diode en silicium

Les chercheurs ont construit un modèle virtuel 2D d'un minuscule capteur en silicium (une diode) épais de 50 micromètres (environ la largeur d'un cheveu humain). Ils ont simulé deux scénarios :

1. Capteur neuf : Avant toute exposition à la radiation.
2. Capteur irradié : Après avoir été frappé par une quantité massive de radiation (simulant l'environnement hostile du HL-LHC).

Ils ont testé ces capteurs à deux températures : une fraîcheur de 248 K (environ -25 °C) et une chaleur de 300 K (température ambiante).

Les résultats : les prévisionnistes sont-ils d'accord ?

1. Le capteur neuf (non irradié)
Lorsque le capteur était tout neuf, les deux programmes informatiques étaient en accord presque parfait sur la quantité d'électricité qui y circulait et sur la façon dont il stockait la charge, jusqu'à environ 500 volts.

• La divergence : Lorsqu'ils ont poussé la tension très haut (près de 700 volts), les programmes ont commencé à différer légèrement sur le moment exact où le capteur « casserait » (claquage). Les auteurs suggèrent que cela est probablement dû au fait que les deux programmes utilisent des « grilles » (maillages) numériques légèrement différentes pour dessiner le capteur, de la même manière que deux applications de cartes différentes pourraient tracer une route légèrement différemment.

2. Le capteur irradié (le vrai test)
C'est là que la vraie magie opère. Ils ont simulé le capteur après qu'il ait été bombardé par des radiations.

• Courant de fuite : Les deux programmes ont prédit la « fuite » (électricité indésirable) de manière presque identique.
• Tension de déplétion : Les deux ont été en parfait accord sur la quantité de tension nécessaire pour faire fonctionner à nouveau le capteur.
• Champs électriques : Ils ont cartographié les forces électriques invisibles à l'intérieur du silicium. Au milieu du capteur (le « volume »), les deux programmes correspondaient presque parfaitement (à moins de 1 % l'un de l'autre).
• Les « pièges » : Ils ont également examiné les minuscules « pièges » créés par la radiation qui capturent les électrons. Les deux programmes étaient en accord sur le comportement de ces pièges dans une marge très raisonnable (environ 20 %).

La touche de température :
À température ambiante (300 K), les programmes ont divergé un peu plus aux niveaux de radiation les plus élevés. Cependant, les auteurs notent que ce n'est pas une grande préoccupation car, dans le monde réel, ces capteurs endommagés ne sont presque jamais exploités à température ambiante ; ils sont maintenus très froids pour survivre. Ainsi, c'est l'accord à basse température (248 K) qui compte vraiment, et là, les deux programmes étaient parfaitement synchronisés.

La conclusion

L'article conclut que Synopsys et Silvaco sont comme deux chefs différents suivant exactement la même recette (le modèle Perugia) et aboutissant au même plat délicieux.

Même si les outils logiciels sont différents, lorsqu'ils utilisent le modèle de dommages par radiation Perugia, ils produisent des prédictions presque identiques sur la façon dont les capteurs en silicium survivront aux radiations sévères du futur HL-LHC. Cela donne aux scientifiques la confiance que leurs modèles sont solides et qu'ils peuvent utiliser l'un ou l'autre outil pour concevoir la prochaine génération de détecteurs de particules.

Note : Les auteurs mentionnent qu'ils prévoient d'examiner la « charge collectée » à l'avenir, mais cet article s'est concentré strictement sur la tension, le courant et les champs électriques.

Résumé technique

1. Énoncé du problème

Le Grand collisionneur de hadrons à haute luminosité (HL-LHC) exposera les détecteurs à pixels en silicium à des fluences de radiation 5 à 10 fois supérieures aux niveaux actuels. Cette radiation intense provoque des dommages en volume et en surface dans les capteurs en silicium, entraînant :

• Une augmentation du courant de fuite.
• Des décalages de la tension de fonctionnement (de déplétion).
• Une perte de signal et une dégradation de la collecte de charge.
• Des modifications du comportement de l'électronique de lecture.

Une prédiction précise de ces effets est cruciale pour estimer les tensions de fonctionnement et entraîner les algorithmes de reconstruction. Bien que le modèle de dommages par radiation de l'Université de Pérouse (intégrant une physique spécifique des défauts en volume et en surface) ait été mis en œuvre avec succès dans Synopsys Sentaurus TCAD, il était nécessaire de valider sa portabilité et son pouvoir prédictif dans Silvaco TCAD. Sans une validation croisée entre ces deux outils de simulation standards de l'industrie, la robustesse de la modélisation des dommages par radiation reste non vérifiée pour une adoption plus large.

2. Méthodologie

Les auteurs ont réalisé une étude comparative de simulations de dispositifs utilisant les outils TCAD Synopsys et Silvaco.

• Structure du dispositif : Une simulation 2D d'une diode n-on-p a été utilisée.
  • Épaisseur : 50 µm.
  • Concentration en volume : $5 \times 10^{12} \text{ cm}^{-3}$ (Bore).
  • Pas : 20 µm avec une petite électrode de collecte.
  • Surface : couche d'oxyde de 1 µm.
• Modèles physiques : Les deux outils ont utilisé des modèles physiques identiques pour assurer une comparaison équitable :
  • Statistiques : Fermi-Dirac.
  • Recombinaison : Scharfetter - Shockley Read Hall.
  • Tunneling bande-à-bande : Schenk.
  • Rétrécissement de la bande interdite : Slotboom.
  • Ionisation par impact : van Overstraeten et de Man.
• Mise en œuvre des dommages par radiation : Le modèle de Pérouse a été transféré vers Silvaco. Cela inclut :
  • Modèle en volume : Deux défauts accepteurs et un défaut donneur avec des densités proportionnelles à la fluence ($\Phi$). Il utilise une fonction dépendante de la fluence pour la création d'accepteurs ($N_A$) afin de gérer les fortes fluences ($> 3 \times 10^{15} \text{ neq/cm}^2$).
  • Modèle de surface : Une distribution continue d'énergie de défauts (de type accepteur dans la partie supérieure de la bande interdite, de type donneur près de la bande de valence) et une augmentation linéaire de la charge d'oxyde avec la dose.
• Scénarios de simulation :
  • Fluences : Non irradié, $\Phi = 1.5 \times 10^{15} \text{ neq/cm}^2$, et $\Phi = 1.0 \times 10^{16} \text{ neq/cm}^2$.
  • Températures : $T = 248 \text{ K}$ et $T = 300 \text{ K}$.
• Observables : Caractéristiques Courant-Tension (IV), Capacité-Tension (CV), profils de champ électrique et probabilités d'ionisation des pièges.

3. Contributions clés

• Portabilité du modèle : Transfert réussi des modèles complexes de dommages par radiation en volume et en surface de Pérouse de Synopsys vers Silvaco TCAD.
• Validation croisée : Fourniture de la première comparaison directe de ces deux moteurs TCAD distincts lors de la simulation d'environnements à haute radiation pertinents pour le HL-LHC.
• Analyse des statistiques de défauts : Analyse non seulement des caractéristiques électriques macroscopiques, mais aussi des probabilités d'ionisation microscopiques des pièges, offrant un aperçu plus profond de la cohérence de la mise en œuvre de la physique des défauts.

4. Résultats

Dispositif non irradié

• Caractéristiques IV : Excellent accord entre les outils jusqu'à ~500 V. Des écarts sont apparus à des tensions plus élevées (Synopsys a prédit une rupture à ~700 V contre ~760 V pour Silvaco), probablement dus à des différences dans la génération de maillage.
• Caractéristiques CV : Les deux outils ont prédit avec précision une tension de déplétion d'environ 10 V et ont correspondu à la caractéristique « marche » à ~55 V (attribuée à la déplétion sous l'oxyde).

Dispositif irradié

• Courant de fuite (IV) :
  • À 248 K : Accord remarquablement bon pour toutes les fluences. De légères déviations se sont produites uniquement à des tensions très élevées et à la fluence la plus élevée, attribuées à des différences de maillage.
  • À 300 K : L'accord était uniquement qualitatif à la fluence la plus élevée. Cependant, les auteurs notent que cela est moins critique car les dispositifs sont rarement exploités au-dessus de 0 °C après une telle irradiation élevée.
• Tension de déplétion (CV) : À 248 K, l'accord entre les deux outils était « pratiquement parfait » pour les deux fluences.
• Profils de champ électrique :
  • Dans la région en volume, les profils de champ électrique étaient en accord à moins de 1 % pour les deux fluences et toutes les tensions de polarisation (100 V, 300 V, 500 V).
  • De légères divergences ont été observées près des implants avant et arrière, à nouveau attribuées à des incohérences de maillage et à des différences de modélisation des implants.
• Probabilités d'ionisation des pièges :
  • Pour les défauts en volume (pièges accepteurs 1 et 2), les probabilités d'ionisation étaient en accord à un ordre de grandeur, et à moins de 20 % dans la région en volume.
  • Ce niveau d'accord est considéré comme satisfaisant compte tenu de l'accord au niveau du pourcentage dans les champs électriques.

5. Importance

• Validation du modèle de Pérouse : L'étude confirme que le modèle de dommages par radiation de l'Université de Pérouse est robuste et indépendant de la plateforme, produisant des résultats cohérents dans différents environnements logiciels TCAD.
• Fiabilité pour le HL-LHC : Le haut degré d'accord dans les observables critiques (courant de fuite, tension de déplétion et champ électrique) apporte confiance dans l'utilisation de ces simulations pour prédire les performances des capteurs dans les conditions du HL-LHC.
• Conduite opérationnelle : Les résultats soutiennent l'utilisation de ces modèles pour estimer les tensions de fonctionnement sûres et comprendre les mécanismes de perte de signal, ce qui est essentiel pour la conception et les algorithmes de reconstruction des trajectographes ATLAS et CMS mis à niveau.
• Travaux futurs : Les auteurs prévoient d'étendre cette validation aux simulations de charge collectée, qui est la métrique ultime pour les performances de trajectographie.

Conclusion : Le papier démontre que les outils TCAD Silvaco et Synopsys peuvent produire des prédictions mutuellement cohérentes pour les détecteurs à pixels en silicium irradiés lorsqu'ils utilisent le modèle de dommages par radiation de Pérouse, validant la solidité du modèle pour les futures applications de collisionneurs à haute luminosité.