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

Dieser Artikel vergleicht Silvaco- und Synopsys-TCAD-Simulationen, die das Perugia-Strahlenschädigungsmodell integrieren, um ihre Vorhersagegenauigkeit für wesentliche Leistungsparameter von Silizium-Pixeldetektoren unter den extremen Strahlungsbedingungen zu bewerten, die am High-Luminosity-LHC erwartet werden.

Das Wesentliche

Stellen Sie sich den Large Hadron Collider (LHC) als eine massive, Hochgeschwindigkeits-Teilchenrennstrecke vor. Wissenschaftler rüsten ihn auf die „High-Luminosity"-Version (HL-LHC) auf, was bedeutet, dass sie Teilchen viel häufiger gegeneinander prallen lassen werden. Das Problem? Dieser intensive Verkehr erzeugt viel „Strahlungsstaub", der die Siliziumsensoren (die Kameras) beschädigt, die versuchen, die Kollisionen zu fotografieren.

Im Laufe der Zeit verwandelt dieser Strahlungsstaub die Siliziumsensoren in „undichte" und „steife" Bauteile. Sie beginnen, ihre Fähigkeit zu verlieren, Signale zu sammeln (wie eine Kamera, die den Fokus verliert), und benötigen deutlich höhere Spannungen, um zu funktionieren, was das Risiko ihres Zerbruchs birgt.

Um dies zu beheben, bevor es geschieht, nutzen Wissenschaftler Computersimulationen, um vorherzusagen, wie sich die Sensoren nach Jahren der Strahlung verhalten werden. Sie müssen wissen: Wie viel Spannung benötigen wir? Wie viel Strom wird undicht? Wird der Sensor noch funktionieren?

Die zwei „Wettervorhersager"

In diesem Papier testen die Forscher zwei verschiedene Computerprogramme (TCAD-Tools), die wie Wettervorhersager für diese Sensoren fungieren:

1. Synopsys
2. Silvaco

Beide Programme verwenden einen spezifischen Satz von Regeln namens „Perugia Radiation Damage Model". Betrachten Sie dieses Modell als detaillierte Bedienungsanleitung, die dem Computer genau sagt, wie der „Strahlungsstaub" das Silizium beschädigt, winzige Fallen und Löcher erzeugt, die den elektrischen Fluss durcheinanderbringen.

Das Ziel dieses Papiers ist es, herauszufinden, ob diese beiden verschiedenen „Vorhersager" bei Verwendung derselben Bedienungsanleitung dieselbe Vorhersage treffen. Wenn sie übereinstimmen, bedeutet dies, dass die Anleitung zuverlässig ist, und Wissenschaftler können den Vorhersagen vertrauen, unabhängig davon, welche Software sie verwenden.

Das Experiment: Eine winzige Siliziumdiode

Die Forscher bauten ein virtuelles 2D-Modell eines winzigen Siliziumsensors (einer Diode) mit einer Dicke von 50 Mikrometern (etwa so breit wie ein menschliches Haar). Sie simulierten zwei Szenarien:

1. Frischer Sensor: Bevor ihn Strahlung trifft.
2. Bestrahlter Sensor: Nachdem er von einer massiven Menge Strahlung getroffen wurde (was die harsche Umgebung des HL-LHC simuliert).

Sie testeten diese Sensoren bei zwei Temperaturen: einer kühlen 248 K (etwa -25°C) und einer warmen 300 K (Raumtemperatur).

Die Ergebnisse: Stimmen die Vorhersager überein?

1. Der frische Sensor (Unbestrahlt)
Als der Sensor brandneu war, stimmten beide Computerprogramme bis zu etwa 500 Volt fast perfekt darin überein, wie viel Elektrizität durch ihn floss und wie er Ladung speicherte.

• Die Abweichung: Als sie die Spannung sehr hoch trieben (nahe 700 Volt), begannen die Programme leicht zu divergieren, genau wann der Sensor „zerbrechen" würde (Durchbruch). Die Autoren vermuten, dass dies wahrscheinlich daran liegt, dass die beiden Programme leicht unterschiedliche digitale „Gitter" (Meshes) verwenden, um den Sensor zu zeichnen, ähnlich wie zwei verschiedene Karten-Apps eine Straße leicht unterschiedlich zeichnen könnten.

2. Der bestrahlte Sensor (Der eigentliche Test)
Hier geschah die eigentliche Magie. Sie simulierten den Sensor, nachdem er mit Strahlung bombardiert worden war.

• Leckstrom: Beide Programme sagten den „Leckstrom" (unerwünschte Elektrizität) fast identisch voraus.
• Verarmungsspannung: Beide stimmten perfekt darin überein, wie viel Spannung benötigt wurde, um den Sensor wieder funktionsfähig zu machen.
• Elektrische Felder: Sie kartierten die unsichtbaren elektrischen Kräfte im Inneren des Siliziums. In der Mitte des Sensors (dem „Bulk") stimmten die beiden Programme fast perfekt überein (innerhalb von 1 % zueinander).
• Die „Fallen": Sie untersuchten auch die winzigen „Fallen", die durch Strahlung erzeugt werden und Elektronen einfangen. Die beiden Programme stimmten im Verhalten dieser Fallen innerhalb eines sehr vernünftigen Margins überein (etwa 20 %).

Die Temperatur-Wende:
Bei Raumtemperatur (300 K) stimmten die Programme bei den höchsten Strahlungsniveaus etwas weniger überein. Die Autoren stellen jedoch fest, dass dies keine große Sorge ist, da diese beschädigten Sensoren in der realen Welt fast nie bei Raumtemperatur betrieben werden; sie werden sehr kalt gehalten, um zu überleben. Daher ist die Übereinstimmung bei der kalten Temperatur (248 K) das, was wirklich zählt, und dort waren die beiden Programme perfekt synchron.

Das Fazit

Das Papier kommt zu dem Schluss, dass Synopsys und Silvaco wie zwei verschiedene Köche sind, die exakt dasselbe Rezept (das Perugia-Modell) befolgen und am Ende dasselbe köstliche Gericht servieren.

Obwohl die Software-Tools unterschiedlich sind, produzieren sie bei Verwendung des Perugia-Strahlungsschadensmodells nahezu identische Vorhersagen darüber, wie Siliziumsensoren der harschen Strahlung des zukünftigen HL-LHC standhalten werden. Dies gibt Wissenschaftlern das Vertrauen, dass ihre Modelle solide sind und dass sie entweder Tool verwenden können, um die nächste Generation von Teilchendetektoren zu entwerfen.

Hinweis: Die Autoren erwähnen, dass sie in Zukunft die „gesammelte Ladung" untersuchen wollen, aber dieses Papier konzentrierte sich streng auf Spannung, Strom und elektrische Felder.

Technische Zusammenfassung

1. Problemstellung

Der High-Luminosity Large Hadron Collider (HL-LHC) wird Silizium-Pixeldetektoren einer Strahlungsflussdichte aussetzen, die das Fünf- bis Zehnfache der aktuellen Werte beträgt. Diese intensive Strahlung verursacht Volumen- und Oberflächenschäden in Siliziumsensoren, was zu folgenden Effekten führt:

• Erhöhter Leckstrom.
• Verschiebungen der Betriebs- (Verarmungs-)Spannung.
• Signalverlust und Verschlechterung der Ladungssammlung.
• Änderungen im Verhalten der Frontend-Elektronik.

Eine genaue Vorhersage dieser Effekte ist entscheidend für die Abschätzung der Betriebsspannungen und das Training von Rekonstruktionsalgorithmen. Während das Strahlenschadensmodell der Universität Perugia (das spezifische Volumen- und Oberflächendefektphysik integriert) erfolgreich in Synopsys Sentaurus TCAD implementiert wurde, bestand die Notwendigkeit, seine Portierbarkeit und Vorhersagekraft in Silvaco TCAD zu validieren. Ohne eine Kreuzvalidierung zwischen diesen beiden branchenstandardkonformen Simulationstools bleibt die Robustheit der Strahlenschadensmodellierung für eine breitere Anwendung unbestätigt.

2. Methodik

Die Autoren führten eine vergleichende Studie von Bauelementesimulationen mit beiden Synopsys- und Silvaco-TCAD-Tools durch.

• Bauelementstruktur: Eine 2D-Simulation einer n-on-p-Diode wurde verwendet.
  • Dicke: 50 µm.
  • Volumenkonzentration: $5 \times 10^{12} \text{ cm}^{-3}$ (Bor).
  • Pitch: 20 µm mit einer kleinen Sammelelektrode.
  • Oberfläche: 1 µm Oxidschicht.
• Physikmodelle: Beide Tools verwendeten identische physikalische Modelle, um einen fairen Vergleich zu gewährleisten:
  • Statistik: Fermi-Dirac.
  • Rekombination: Scharfetter - Shockley Read Hall.
  • Band-zu-Band-Tunneln: Schenk.
  • Bandlückenverkleinerung: Slotboom.
  • Stoßionisation: van Overstraeten und de Man.
• Implementierung der Strahlenschäden: Das Perugia-Modell wurde auf Silvaco portiert. Dies umfasst:
  • Volumenmodell: Zwei Akzeptordefekte und ein Donordefekt mit Dichten, die proportional zur Flussdichte ($\Phi$) sind. Es wird eine flussdichteabhängige Funktion für die Akzeptorerzeugung ($N_A$) verwendet, um hohe Flussdichten ($> 3 \times 10^{15} \text{ neq/cm}^2$) zu behandeln.
  • Oberflächenmodell: Eine kontinuierliche Energieverteilung von Defekten (akzeptorähnlich im oberen Bandlückenbereich, donorähnlich in der Nähe des Valenzbands) und ein linearer Anstieg der Oxidladung mit der Dosis.
• Simulationsszenarien:
  • Flussdichten: Unbestrahlt, $\Phi = 1.5 \times 10^{15} \text{ neq/cm}^2$ und $\Phi = 1.0 \times 10^{16} \text{ neq/cm}^2$.
  • Temperaturen: $T = 248 \text{ K}$ und $T = 300 \text{ K}$.
• Beobachtbare Größen: Strom-Spannungs-Kennlinien (IV), Kapazitäts-Spannungs-Kennlinien (CV), elektrische Feldprofile und Ionisationswahrscheinlichkeiten von Fallen.

3. Hauptbeiträge

• Modellportierbarkeit: Erfolgreiche Portierung der komplexen Perugia-Volumen- und Oberflächenstrahlenschadensmodelle von Synopsys zu Silvaco TCAD.
• Kreuzvalidierung: Bereitstellung des ersten direkten Vergleichs dieser beiden unterschiedlichen TCAD-Engines bei der Simulation hochstrahlungsintensiver Umgebungen, die für den HL-LHC relevant sind.
• Analyse der Defektstatistik: Analyse nicht nur makroskopischer elektrischer Kenngrößen, sondern auch mikroskopischer Ionisationswahrscheinlichkeiten von Fallen, was tiefere Einblicke in die Konsistenz der Implementierung der Defektphysik bietet.

4. Ergebnisse

Unbestrahltes Bauelement

• IV-Kennlinien: Hervorragende Übereinstimmung zwischen den Tools bis zu ca. 500 V. Abweichungen traten bei höheren Spannungen auf (Synopsys sagte einen Durchbruch bei ca. 700 V voraus, Silvaco bei ca. 760 V), wahrscheinlich aufgrund von Unterschieden bei der Netzgenerierung.
• CV-Kennlinien: Beide Tools sagten präzise eine Verarmungsspannung von ca. 10 V voraus und stimmten mit dem „Schritt"-Merkmals bei ca. 55 V überein (zugeschrieben der Verarmung unter dem Oxid).

Bestrahltes Bauelement

• Leckstrom (IV):
  • Bei 248 K: Bemerkenswert gute Übereinstimmung über alle Flussdichten hinweg. Geringe Abweichungen traten nur bei sehr hohen Spannungen und höchster Flussdichte auf, was auf Netzunterschiede zurückgeführt wird.
  • Bei 300 K: Die Übereinstimmung war bei höchster Flussdichte nur qualitativ. Die Autoren stellen jedoch fest, dass dies weniger kritisch ist, da Bauelemente nach einer derart hohen Bestrahlung selten oberhalb von 0°C betrieben werden.
• Verarmungsspannung (CV): Bei 248 K war die Übereinstimmung zwischen den beiden Tools für beide Flussdichten „praktisch perfekt".
• Elektrische Feldprofile:
  • Im Volumenbereich stimmten die elektrischen Feldprofile für beide Flussdichten und alle Vorspannungen (100 V, 300 V, 500 V) innerhalb von 1% überein.
  • Leichte Abweichungen wurden in der Nähe von Vorder- und Rückseitenimplantationen beobachtet, wiederum zugeschrieben Netzfehlanpassungen und Unterschieden bei der Implantationsmodellierung.
• Ionisationswahrscheinlichkeiten von Fallen:
  • Für Volumenfehler (Akzeptorfallen 1 und 2) stimmten die Ionisationswahrscheinlichkeiten innerhalb einer Größenordnung überein und innerhalb von 20% im Volumenbereich.
  • Dieses Niveau der Übereinstimmung wird angesichts der prozentualen Übereinstimmung der elektrischen Felder als zufriedenstellend betrachtet.

5. Bedeutung

• Validierung des Perugia-Modells: Die Studie bestätigt, dass das Strahlenschadensmodell der Universität Perugia robust und plattformunabhängig ist und konsistente Ergebnisse über verschiedene TCAD-Softwareumgebungen hinweg liefert.
• Zuverlässigkeit für den HL-LHC: Der hohe Grad der Übereinstimmung bei kritischen beobachtbaren Größen (Leckstrom, Verarmungsspannung und elektrisches Feld) gibt Vertrauen in die Verwendung dieser Simulationen zur Vorhersage der Sensorleistung unter HL-LHC-Bedingungen.
• Betriebliche Leitlinien: Die Ergebnisse unterstützen die Verwendung dieser Modelle zur Abschätzung sicherer Betriebsspannungen und zum Verständnis von Signalverlustmechanismen, was für das Design und die Rekonstruktionsalgorithmen der upgegradeten ATLAS- und CMS-Spurdetektoren essenziell ist.
• Zukünftige Arbeiten: Die Autoren planen, diese Validierung auf Simulationen der gesammelten Ladung auszudehnen, was das ultimative Maß für die Spurenerfassungsleistung darstellt.

Fazit: Das Papier zeigt, dass Silvaco- und Synopsys-TCAD-Tools bei Verwendung des Perugia-Strahlenschadensmodells gegenseitig konsistente Vorhersagen für bestrahlte Silizium-Pixeldetektoren liefern können und damit die Solidität des Modells für zukünftige Anwendungen bei Hoch-Luminositäts-Collidern validieren.