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

Questo articolo confronta le simulazioni TCAD di Silvaco e Synopsys che incorporano il modello di danno da radiazione Perugia per valutarne l'accuratezza predittiva rispetto alle metriche chiave delle prestazioni dei rivelatori a pixel di silicio nelle condizioni di radiazione estrema previste per l'High Luminosity LHC.

L'essenziale

Immaginate il Large Hadron Collider (LHC) come un'enorme pista di corse per particelle ad alta velocità. Gli scienziati lo stanno aggiornando alla versione "High-Luminosity" (HL-LHC), il che significa che faranno scontrare le particelle con molta più frequenza. Il problema? Questo traffico intenso genera molta "polvere di radiazioni" che danneggia i sensori al silicio (le fotocamere) che cercano di catturare le immagini delle collisioni.

Nel tempo, questa polvere di radiazioni trasforma i sensori al silicio in dispositivi "perdenti" e "rigidi". Iniziano a perdere la capacità di raccogliere segnali (come una fotocamera che perde la messa a fuoco) e richiedono tensioni molto più elevate per funzionare, il che rischia di romperli.

Per risolvere il problema prima che si verifichi, gli scienziati utilizzano simulazioni al computer per prevedere come si comporteranno i sensori dopo anni di radiazioni. Devono sapere: Quanta tensione ci serve? Quanta corrente perderà? Il sensore funzionerà ancora?

I due "meteorologi"

In questo articolo, i ricercatori stanno testando due diversi programmi informatici (strumenti TCAD) che agiscono come meteorologi per questi sensori:

1. Synopsys
2. Silvaco

Entrambi i programmi utilizzano un insieme specifico di regole chiamato "Modello di Danno da Radiazione Perugia". Pensate a questo modello come a un manuale di istruzioni dettagliato che dice al computer esattamente come la "polvere di radiazioni" danneggia il silicio, creando piccole trappole e vuoti che disturbano il flusso elettrico.

L'obiettivo di questo articolo è vedere se questi due diversi "meteorologi" danno la stessa previsione quando utilizzano lo stesso manuale di istruzioni. Se sono d'accordo, significa che il manuale è affidabile e gli scienziati possono fidarsi delle previsioni indipendentemente dal software utilizzato.

L'esperimento: un piccolo diodo al silicio

I ricercatori hanno costruito un modello virtuale bidimensionale di un minuscolo sensore al silicio (un diodo) spesso 50 micrometri (circa la larghezza di un capello umano). Hanno simulato due scenari:

1. Sensore fresco: Prima che lo colpisse qualsiasi radiazione.
2. Sensore irradiato: Dopo essere stato colpito da una massa enorme di radiazioni (simulando l'ambiente ostile dell'HL-LHC).

Hanno testato questi sensori a due temperature: una fresca di 248 K (circa -25°C) e una calda di 300 K (temperatura ambiente).

I risultati: sono d'accordo i meteorologi?

1. Il sensore fresco (non irradiato)
Quando il sensore era nuovo di zecca, entrambi i programmi informatici erano in perfetto accordo su quanta elettricità fluiva attraverso di esso e su come immagazzinava la carica, fino a circa 500 volt.

• La discrepanza: Quando hanno spinto la tensione molto in alto (vicino a 700 volt), i programmi hanno iniziato a non essere d'accordo leggermente su esattamente quando il sensore si sarebbe "rotto" (rottura). Gli autori suggeriscono che ciò è probabilmente dovuto al fatto che i due programmi utilizzano "griglie" (mesh) digitali leggermente diverse per disegnare il sensore, simile a come due diverse app di mappe potrebbero disegnare una strada in modo leggermente diverso.

2. Il sensore irradiato (il vero test)
È qui che è avvenuta la vera magia. Hanno simulato il sensore dopo essere stato bombardato da radiazioni.

• Corrente di dispersione: Entrambi i programmi hanno previsto la "dispersione" (elettricità indesiderata) quasi identicamente.
• Tensione di svuotamento: Entrambi sono stati perfettamente d'accordo su quanta tensione fosse necessaria per far funzionare di nuovo il sensore.
• Campi elettrici: Hanno mappato le forze elettriche invisibili all'interno del silicio. Nel mezzo del sensore (il "bulk"), i due programmi corrispondevano quasi perfettamente (entro l'1% l'uno dall'altro).
• Le "trappole": Hanno anche esaminato le piccole "trappole" create dalle radiazioni che catturano gli elettroni. I due programmi erano d'accordo sul comportamento di queste trappole entro un margine molto ragionevole (circa il 20%).

La svolta della temperatura:
A temperatura ambiente (300 K), i programmi non erano d'accordo un po' di più ai livelli di radiazione più alti. Tuttavia, gli autori notano che questo non è una grande preoccupazione perché, nel mondo reale, questi sensori danneggiati sono quasi mai operati a temperatura ambiente; vengono mantenuti molto freddi per sopravvivere. Quindi, l'accordo alla temperatura fredda (248 K) è ciò che conta davvero, e lì i due programmi erano perfettamente sincronizzati.

La conclusione

L'articolo conclude che Synopsys e Silvaco sono come due chef diversi che seguono esattamente la stessa ricetta (il Modello Perugia) e finiscono con lo stesso piatto delizioso.

Anche se gli strumenti software sono diversi, quando utilizzano il modello di danno da radiazione Perugia, producono previsioni quasi identiche su come i sensori al silicio sopravviveranno alle radiazioni ostili del futuro HL-LHC. Questo dà agli scienziati la fiducia che i loro modelli siano solidi e che possano utilizzare uno qualsiasi dei due strumenti per progettare la prossima generazione di rivelatori di particelle.

Nota: Gli autori menzionano che intendono esaminare la "carica raccolta" in futuro, ma questo articolo si è concentrato strettamente su tensione, corrente e campi elettrici.

Sommario tecnico

1. Enunciato del Problema

Il Large Hadron Collider ad Alta Luminosità (HL-LHC) sottoporrà i rivelatori a pixel di silicio a flussi di radiazione da 5 a 10 volte superiori ai livelli attuali. Questa intensa radiazione provoca danni nel volume e sulla superficie dei sensori di silicio, portando a:

• Aumento della corrente di dispersione.
• Spostamenti della tensione operativa (di svuotamento).
• Perdita di segnale e degradazione della raccolta della carica.
• Modifiche al comportamento dell'elettronica di front-end.

Una previsione accurata di questi effetti è fondamentale per stimare le tensioni operative e addestrare gli algoritmi di ricostruzione. Sebbene il Modello di Danno da Radiazione dell'Università di Perugia (che incorpora fisica specifica dei difetti volumetrici e superficiali) sia stato implementato con successo in Synopsys Sentaurus TCAD, era necessario convalidare la sua portabilità e capacità predittiva in Silvaco TCAD. Senza una convalida incrociata tra questi due principali strumenti di simulazione standard del settore, la robustezza della modellazione del danno da radiazione rimaneva non verificata per un'adozione più ampia.

2. Metodologia

Gli autori hanno eseguito uno studio comparativo delle simulazioni di dispositivo utilizzando entrambi gli strumenti TCAD Synopsys e Silvaco.

• Struttura del Dispositivo: È stata utilizzata una simulazione 2D di un diodo n-on-p.
  • Spessore: 50 µm.
  • Concentrazione di volume: $5 \times 10^{12} \text{ cm}^{-3}$ (Boro).
  • Passo: 20 µm con un piccolo elettrodo di raccolta.
  • Superficie: strato di ossido di 1 µm.
• Modelli Fisici: Entrambi gli strumenti hanno utilizzato modelli fisici identici per garantire un confronto equo:
  • Statistiche: Fermi-Dirac.
  • Ricombinazione: Scharfetter - Shockley Read Hall.
  • Tunneling banda-banda: Schenk.
  • Restringimento del bandgap: Slotboom.
  • Ionizzazione per impatto: van Overstraeten e de Man.
• Implementazione del Danno da Radiazione: Il modello di Perugia è stato trasferito su Silvaco. Ciò include:
  • Modello di Volume: Due difetti accettori e un difetto donatore con densità proporzionali al flusso ($\Phi$). Utilizza una funzione dipendente dal flusso per la creazione di accettori ($N_A$) per gestire flussi elevati ($> 3 \times 10^{15} \text{ neq/cm}^2$).
  • Modello di Superficie: Una distribuzione continua di energia dei difetti (simili ad accettori nella parte superiore del bandgap, simili a donatori vicino alla banda di valenza) e aumento lineare della carica nell'ossido con la dose.
• Scenari di Simulazione:
  • Flussi: Non irradiato, $\Phi = 1.5 \times 10^{15} \text{ neq/cm}^2$, e $\Phi = 1.0 \times 10^{16} \text{ neq/cm}^2$.
  • Temperature: $T = 248 \text{ K}$ e $T = 300 \text{ K}$.
• Osservabili: Corrente-Tensione (IV), Capacità-Tensione (CV), profili del campo elettrico e probabilità di ionizzazione delle trappole.

3. Contributi Chiave

• Portabilità del Modello: Trasferimento riuscito dei complessi modelli di danno da radiazione volumetrico e superficiale di Perugia da Synopsys a Silvaco TCAD.
• Convalida Incrociata: Fornita la prima comparazione diretta di questi due motori TCAD distinti nella simulazione di ambienti ad alta radiazione rilevanti per l'HL-LHC.
• Analisi delle Statistiche dei Difetti: Analizzate non solo le caratteristiche elettriche macroscopiche, ma anche le probabilità di ionizzazione delle trappole microscopiche, offrendo una visione più approfondita della coerenza nell'implementazione della fisica dei difetti.

4. Risultati

Dispositivo Non Irradiato

• Caratteristiche IV: Ottimo accordo tra gli strumenti fino a ~500 V. Sono comparse discrepanze a tensioni più elevate (Synopsys ha previsto la rottura a ~700 V contro ~760 V di Silvaco), probabilmente dovute a differenze nella generazione della mesh.
• Caratteristiche CV: Entrambi gli strumenti hanno previsto accuratamente una tensione di svuotamento di ~10 V e hanno corrisposto alla caratteristica "gradino" a ~55 V (attribuita allo svuotamento sotto l'ossido).

Dispositivo Irradiato

• Corrente di Dispersione (IV):
  • A 248 K: Accordi notevolmente buoni su tutti i flussi. Minime deviazioni si sono verificate solo a tensioni molto elevate e al flusso più alto, attribuite a differenze nella mesh.
  • A 300 K: L'accordo è stato solo qualitativo al flusso più alto. Tuttavia, gli autori notano che questo è meno critico poiché i dispositivi vengono raramente operati sopra 0°C dopo tale alta irradiazione.
• Tensione di Svuotamento (CV): A 248 K, l'accordo tra i due strumenti è stato "praticamente perfetto" per entrambi i flussi.
• Profili del Campo Elettrico:
  • Nella regione di volume, i profili del campo elettrico hanno concordato entro 1% per entrambi i flussi e per tutte le tensioni di polarizzazione (100 V, 300 V, 500 V).
  • Sono state osservate lievi discrepanze vicino agli impianti frontali e posteriori, nuovamente attribuite a disallineamenti della mesh e differenze nella modellazione degli impianti.
• Probabilità di Ionizzazione delle Trappole:
  • Per i difetti di volume (trappole accettrici 1 e 2), le probabilità di ionizzazione hanno concordato entro un ordine di grandezza e entro 20% nella regione di volume.
  • Questo livello di accordo è considerato soddisfacente data l'accordo a livello percentuale nei campi elettrici.

5. Significato

• Convalida del Modello di Perugia: Lo studio conferma che il modello di danno da radiazione dell'Università di Perugia è robusto e indipendente dalla piattaforma, producendo risultati coerenti attraverso diversi ambienti software TCAD.
• Affidabilità per l'HL-LHC: L'alto grado di accordo negli osservabili critici (corrente di dispersione, tensione di svuotamento e campo elettrico) fornisce fiducia nell'uso di queste simulazioni per prevedere le prestazioni dei sensori nelle condizioni dell'HL-LHC.
• Linee Guida Operative: I risultati supportano l'uso di questi modelli per stimare le tensioni operative sicure e comprendere i meccanismi di perdita di segnale, essenziale per la progettazione e gli algoritmi di ricostruzione dei tracciatori ATLAS e CMS aggiornati.
• Lavori Futuri: Gli autori pianificano di estendere questa convalida alle simulazioni di carica raccolta, che è la metrica definitiva per le prestazioni di tracciamento.

Conclusione: Il documento dimostra che gli strumenti TCAD Silvaco e Synopsys possono produrre previsioni reciprocamente coerenti per rivelatori a pixel di silicio irradiati quando si utilizza il modello di danno da radiazione di Perugia, convalidando la solidità del modello per future applicazioni di collider ad alta luminosità.