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

Dit artikel vergelijkt Silvaco- en Synopsys TCAD-simulaties die het Perugia-stralingschademodel incorporeren om hun voorspellende nauwkeurigheid te evalueren voor belangrijke prestatieparameters van siliciumpixeldetectoren onder de extreme stralingsomstandigheden die op de High Luminosity LHC worden verwacht.

De kern

Stel je de Large Hadron Collider (LHC) voor als een massieve, supersnelle deeltjesrenbaan. Wetenschappers upgraden deze naar de "High-Luminosity"-versie (HL-LHC), wat betekent dat ze deeltjes veel vaker tegen elkaar zullen laten botsen. Het probleem? Deze intense drukte creëert veel "stralingsstof" dat de siliciumsensoren (de camera's) beschadigt die proberen de botsingen vast te leggen.

Na verloop van tijd veranderen deze stralingsdeeltjes de siliciumsensoren in "lekke" en "stijve" apparaten. Ze beginnen hun vermogen om signalen te verzamelen te verliezen (zoals een camera die zijn scherpte verliest) en vereisen een veel hogere spanning om te werken, wat het risico met zich meebrengt dat ze kapotgaan.

Om dit te voorkomen voordat het gebeurt, gebruiken wetenschappers computersimulaties om te voorspellen hoe de sensoren zich na jaren van straling zullen gedragen. Ze moeten weten: Hoeveel spanning hebben we nodig? Hoeveel stroom zal er lekken? Zal de sensor nog steeds werken?

De twee "weervoorspellers"

In dit artikel testen de onderzoekers twee verschillende computerprogramma's (TCAD-tools) die fungeren als weervoorspellers voor deze sensoren:

1. Synopsys
2. Silvaco

Beide programma's gebruiken een specifieke reeks regels genaamd het "Perugia-stralingsbeschadigingsmodel". Denk aan dit model als een gedetailleerde handleiding die de computer precies vertelt hoe de "stralingsstof" het silicium beschadigt, door kleine valkuilen en gaten te creëren die de elektrische stroom verstoren.

Het doel van dit artikel is om te zien of deze twee verschillende "voorspellers" dezelfde voorspelling geven wanneer ze dezelfde handleiding gebruiken. Als ze het eens zijn, betekent dit dat de handleiding betrouwbaar is, en kunnen wetenschappers de voorspellingen vertrouwen, ongeacht welke software ze gebruiken.

Het experiment: Een tiny siliciumdiode

De onderzoekers bouwden een virtueel 2D-model van een tiny siliciumsensor (een diode) die 50 micrometer dik is (ongeveer de breedte van een menselijk haar). Ze simuleerden twee scenario's:

1. Verse sensor: Voordat er straling op valt.
2. Gestraalde sensor: Nadat deze is geraakt door een enorme hoeveelheid straling (wat de harde omgeving van de HL-LHC simuleert).

Ze testten deze sensoren bij twee temperaturen: een koele 248 K (ongeveer -25°C) en een warme 300 K (kamertemperatuur).

De resultaten: Zijn de voorspellers het eens?

1. De verse sensor (niet bestraald)
Toen de sensor gloednieuw was, waren beide computerprogramma's bijna perfect het eens over hoeveel elektriciteit er doorheen stroomde en hoe deze lading opslaat, tot ongeveer 500 volt.

• Het verschil: Toen ze de spanning zeer hoog opdreven (dicht bij 700 volt), begonnen de programma's lichtjes van mening te verschillen over precies wanneer de sensor zou "breken" (doorbraak). De auteurs suggereren dat dit waarschijnlijk komt omdat de twee programma's iets verschillende digitale "roosters" (meshes) gebruiken om de sensor te tekenen, vergelijkbaar met hoe twee verschillende kaart-apps een weg iets anders kunnen tekenen.

2. De gestraalde sensor (de echte test)
Hier gebeurde de echte magie. Ze simuleerden de sensor nadat deze was gebombardeerd met straling.

• Lekkagestroom: Beide programma's voorspelden de "lekstroom" (ongewenste elektriciteit) bijna identiek.
• Uitputtingsvoltage: Beide waren het perfect eens over hoeveel spanning nodig was om de sensor weer werkend te maken.
• Elektrische velden: Ze in kaart gebracht de onzichtbare elektrische krachten binnenin het silicium. In het midden van de sensor (de "bulk") kwamen de twee programma's bijna perfect overeen (binnen 1% van elkaar).
• De "valkuilen": Ze keken ook naar de kleine "valkuilen" die door straling zijn gecreëerd en die elektronen vangen. De twee programma's waren het eens over het gedrag van deze valkuilen binnen een zeer redelijke marge (ongeveer 20%).

De temperatuursdraai:
Bij kamertemperatuur (300 K) waren de programma's het bij de hoogste stralingsniveaus iets minder eens. De auteurs merken echter op dat dit geen grote zorg is, omdat deze beschadigde sensoren in de echte wereld bijna nooit bij kamertemperatuur worden bediend; ze worden zeer koud gehouden om te overleven. De overeenstemming bij de koude temperatuur (248 K) is dus wat echt telt, en daar waren de twee programma's perfect synchroon.

De conclusie

Het artikel concludeert dat Synopsys en Silvaco als twee verschillende koks zijn die precies hetzelfde recept volgen (het Perugia-model) en eindigen met hetzelfde heerlijke gerecht.

Hoewel de softwaretools verschillend zijn, produceren ze, wanneer ze het Perugia-stralingsbeschadigingsmodel gebruiken, bijna identieke voorspellingen over hoe siliciumsensoren de harde straling van de toekomstige HL-LHC zullen overleven. Dit geeft wetenschappers vertrouwen dat hun modellen solide zijn en dat ze elk van de twee tools kunnen gebruiken om de volgende generatie deeltjesdetectoren te ontwerpen.

Opmerking: De auteurs vermelden dat ze van plan zijn om in de toekomst te kijken naar "verzamelde lading", maar dit artikel richtte zich strikt op spanning, stroom en elektrische velden.

Technische samenvatting

1. Probleemstelling

De High-Luminosity Large Hadron Collider (HL-LHC) zal siliciumpixeldetectoren blootstellen aan stralingsfluenzen 5 tot 10 keer hoger dan de huidige niveaus. Deze intense straling veroorzaakt bulk- en oppervlakbeschadiging in siliciumsensoren, wat leidt tot:

• Een toename van de lekstroom.
• Verschuivingen in de operationele (depletie) spanning.
• Signaalverlies en degradatie van ladingscollectie.
• Wijzigingen in het gedrag van de front-end elektronica.

Een nauwkeurige voorspelling van deze effecten is cruciaal voor het schatten van operationele spanningen en het trainen van reconstructiealgoritmen. Hoewel het Perugia-stralingsbeschadigingsmodel van de Universiteit (dat specifieke bulk- en oppervlakdefectfysica omvat) succesvol is geïmplementeerd in Synopsys Sentaurus TCAD, was er behoefte aan validatie van de portabiliteit en voorspellende kracht in Silvaco TCAD. Zonder cross-validatie tussen deze twee belangrijke industriestandaard-simulatiehulpmiddelen blijft de robuustheid van de stralingsbeschadigingsmodellering onbevestigd voor bredere adoptie.

2. Methodologie

De auteurs voerden een vergelijkende studie uit van apparaatsimulaties met behulp van zowel Synopsys- als Silvaco TCAD-hulpmiddelen.

• Apparaatstructuur: Een 2D-simulatie van een n-on-p-diode werd gebruikt.
  • Dikte: 50 µm.
  • Bulkconcentratie: $5 \times 10^{12} \text{ cm}^{-3}$ (Boor).
  • Pitch: 20 µm met een kleine collectie-elektrode.
  • Oppervlak: 1 µm oxide-laag.
• Fysieke modellen: Beide hulpmiddelen gebruikten identieke fysieke modellen om een eerlijke vergelijking te waarborgen:
  • Statistieken: Fermi-Dirac.
  • Recombinatie: Scharfetter - Shockley Read Hall.
  • Band-tot-band tunneling: Schenk.
  • Versmalling van de bandkloof: Slotboom.
  • Impact-ionisatie: van Overstraeten en de Man.
• Implementatie van stralingsbeschadiging: Het Perugia-model werd naar Silvaco geporteerd. Dit omvat:
  • Bulkmodel: Twee acceptor-defecten en één donor-defect met dichtheden evenredig met de fluentie ($\Phi$). Het maakt gebruik van een fluentie-afhankelijke functie voor acceptorcreatie ($N_A$) om hoge fluenzen te hanteren ($> 3 \times 10^{15} \text{ neq/cm}^2$).
  • Oppervlakmodel: Een continue energieverdeling van defecten (acceptor-achtig in het bovenste bandgat, donor-achtig nabij het valentieband) en een lineaire toename van de oxide-lading met de dosis.
• Simulatiescenario's:
  • Fluenzen: Onbestraald, $\Phi = 1.5 \times 10^{15} \text{ neq/cm}^2$, en $\Phi = 1.0 \times 10^{16} \text{ neq/cm}^2$.
  • Temperaturen: $T = 248 \text{ K}$ en $T = 300 \text{ K}$.
• Observabelen: Stroom-Spanning (IV), Capacitantie-Spanning (CV), elektrische veldprofielen en valstroom-ionisatiekansen.

3. Belangrijkste bijdragen

• Modelportabiliteit: Succesvol porteren van de complexe Perugia bulk- en oppervlakstralingsbeschadigingsmodellen van Synopsys naar Silvaco TCAD.
• Cross-validatie: Levering van de eerste directe vergelijking van deze twee verschillende TCAD-engines bij het simuleren van hoge-stralingsomgevingen die relevant zijn voor de HL-LHC.
• Defectstatistieken-analyse: Analyse niet alleen van macroscopische elektrische kenmerken, maar ook van microscopische valstroom-ionisatiekansen, wat dieper inzicht biedt in de consistentie van de implementatie van defectfysica.

4. Resultaten

Onbestraald apparaat

• IV-kenmerken: Uitstekende overeenkomst tussen de hulpmiddelen tot ongeveer 500 V. Afwijkingen verschenen bij hogere spanningen (Synopsys voorspelde doorbraak bij ongeveer 700 V versus Silvaco bij ongeveer 760 V), waarschijnlijk vanwege verschillen in mesh-generatie.
• CV-kenmerken: Beide hulpmiddelen voorspelden nauwkeurig een depletiespanning van ongeveer 10 V en kwamen overeen met het "stap"-kenmerk bij ongeveer 55 V (toegeschreven aan depletie onder het oxide).

Bestraald apparaat

• Lekstroom (IV):
  • Bij 248 K: Opmerkelijk goede overeenkomst over alle fluenzen. Minor afwijkingen traden alleen op bij zeer hoge spanningen en de hoogste fluentie, toegeschreven aan mesh-verschillen.
  • Bij 300 K: Overeenkomst was alleen kwalitatief bij de hoogste fluentie. De auteurs merken echter op dat dit minder kritisch is, aangezien apparaten zelden boven 0°C worden bediend na dergelijke hoge bestraling.
• Depletiespanning (CV): Bij 248 K was de overeenkomst tussen de twee hulpmiddelen "praktisch perfect" voor beide fluenzen.
• Elektrische veldprofielen:
  • In het bulkgebied kwamen de elektrische veldprofielen binnen 1% overeen voor beide fluenzen en alle voorspanningen (100 V, 300 V, 500 V).
  • Lichte afwijkingen werden waargenomen nabij front- en achterzijde-implantaties, opnieuw toegeschreven aan mesh-mismatches en verschillen in implantatiemodellering.
• Valstroom-ionisatiekansen:
  • Voor bulkdefecten (acceptor-valstroom 1 en 2) kwamen de ionisatiekansen overeen binnen een orde van grootte, en binnen 20% in het bulkgebied.
  • Dit niveau van overeenkomst wordt als bevredigend beschouwd gezien de overeenkomst op procent-niveau in elektrische velden.

5. Betekenis

• Validatie van het Perugia-model: De studie bevestigt dat het stralingsbeschadigingsmodel van de Universiteit van Perugia robuust en platformonafhankelijk is, en consistente resultaten oplevert in verschillende TCAD-softwareomgevingen.
• Betrouwbaarheid voor HL-LHC: De hoge mate van overeenkomst in kritieke observabelen (lekstroom, depletiespanning en elektrisch veld) geeft vertrouwen in het gebruik van deze simulaties om sensoreprestaties onder HL-LHC-omstandigheden te voorspellen.
• Operationele richtlijnen: De resultaten ondersteunen het gebruik van deze modellen om veilige operationele spanningen te schatten en signaalverliesmechanismen te begrijpen, wat essentieel is voor het ontwerp en de reconstructiealgoritmen van de opgewaardeerde ATLAS- en CMS-trackers.
• Toekomstig werk: De auteurs plannen om deze validatie uit te breiden naar simulaties van verzamelde lading, wat de ultieme maatstaf is voor trackprestaties.

Conclusie: Het artikel toont aan dat Silvaco- en Synopsys TCAD-hulpmiddelen wederzijds consistente voorspellingen kunnen produceren voor bestraalde siliciumpixeldetectoren wanneer het Perugia-stralingsbeschadigingsmodel wordt gebruikt, waarmee de soundness van het model voor toekomstige toepassingen met hoge-luminositeit-colliders wordt gevalideerd.