Simulation of semiconductor detectors in 3D with SolidStateDetectors.jl

Dit artikel introduceert SolidStateDetectors.jl, een open-source Julia-pakket voor efficiënte 3D-simulatie van velden, ladingsdragers en pulsen in halfgeleiderdetectoren, met name germaniumdetectoren.

De kern

Stel je voor dat je een gigantische, ultra-gevoelige valstrik bouwt om de zeldzaamste gebeurtenissen in het universum te vangen, zoals een spook dat probeert door een muur te lopen. Deze "valstrik" is een detector gemaakt van puur germanium (een halfgeleider, net als in je computerchip, maar dan van kristalzuiverheid).

Deze paper introduceert een nieuw, gratis computerprogramma genaamd SolidStateDetectors.jl. Dit programma is als een ultra-realistische virtuele zandbak voor natuurkundigen. Hier is hoe het werkt, vertaald naar alledaagse taal:

1. De Digitale Zandbak (De Detector)

In het echte leven bouwen wetenschappers deze dure kristallen en testen ze met straling. Maar dat kost tijd en geld. Met dit programma kunnen ze eerst in de computer een exacte kopie van hun detector bouwen.

• De Vorm: Je kunt de detector vormen zoals je wilt: cilindrisch, met een klein puntje erin (een "point-contact"), of zelfs met meerdere segmenten (als een taart in plakken).
• De Omgeving: Het programma houdt rekening met de omgeving. Stel je voor dat je de detector in een vacuümkamer zet, of juist onderdompelt in vloeibare argon (een vloeibaar gas). Het programma rekent uit hoe dit de elektrische velden binnenin het kristal verandert, net zoals water de vorm van een zeepbel verandert.

2. Het Onzichtbare Net (Elektrische Velden)

Zodra de detector aan staat, zit er een elektrisch veld in. Denk hierbij aan een onzichtbaar web van krachten dat door het kristal loopt.

• Als een deeltje (zoals een gammastraal) het kristal raakt, maakt het kleine "elektrische vonkjes" (ladingen) los.
• Het programma rekent precies uit hoe deze vonkjes door dat web worden getrokken. Soms worden ze snel naar de ene kant getrokken, soms langzaam naar de andere, afhankelijk van waar ze vandaan komen.

3. De Dans van de Deeltjes (Drift)

Nu begint het echte werk. De losgemaakte vonkjes (elektronen en gaten) beginnen te "dansen" door het kristal.

• De Simulatie: Het programma laat zien hoe deze vonkjes zich verplaatsen, stap voor stap. Het is alsof je een film maakt van een regenbui die door een labyrint valt, waarbij elke druppel een eigen pad heeft.
• Het Signaal: Wanneer deze vonkjes de randen van het kristal raken (de elektroden), geven ze een signaal af. Dit signaal is een elektrische puls (een piek op een grafiek). Het programma voorspelt precies hoe deze puls eruit ziet: hoe snel hij stijgt, hoe hoog hij wordt en hoe hij weer afloopt.

4. De "Spiegel" en de "Valse Spookjes"

Een van de coolste dingen is dat het programma niet alleen kijkt naar het hoofdsignaal, maar ook naar de "echo's" of spiegelbeelden op de andere delen van de detector.

• Waarom is dit belangrijk? Soms denkt een detector dat er één groot deeltje is, terwijl er eigenlijk twee kleine deeltjes zijn die op verschillende plekken vallen (een "multi-site" gebeurtenis). In zoektochten naar donkere materie of zeldzame kernvervalen (zoals dubbel bètaverval) willen we alleen de échte, zeldzame gebeurtenissen zien en de ruis weggooien.
• Door de vorm van de puls te analyseren (de "pulse-shape"), kunnen wetenschappers zeggen: "Ah, dit is een echte gebeurtenis" of "Nee, dit is ruis". Het programma helpt hen om deze regels te perfectioneren voordat ze de dure echte apparatuur bouwen.

5. Vergelijken met de Realiteit

De auteurs hebben dit programma getest met een echte detector. Ze hebben gemeten hoe de echte detector reageerde op straling en vergeleken dit met wat het programma voorspelde.

• Het Resultaat: Het klopte bijna perfect! De digitale simulatie zag eruit als een spiegelbeeld van de echte metingen. Dit betekent dat wetenschappers nu vol vertrouwen kunnen zeggen: "Als we deze detector zo bouwen, zal hij zich zo gedragen."

Waarom is dit geweldig?

Vroeger waren zulke simulaties moeilijk, duur of niet openbaar. Dit programma is:

• Gratis en open: Iedereen mag het gebruiken en verbeteren.
• Snel: Het is geschreven in een moderne programmeertaal (Julia) die razendsnel is, zelfs voor complexe 3D-berekeningen.
• Flexibel: Je kunt het aanpassen voor verschillende soorten materialen en vormen.

Kortom: Dit paper introduceert een krachtig nieuw gereedschap dat het mogelijk maakt om in de computer te "proefdraaien" met de bouw van de meest gevoelige stralingsdetectoren ter wereld, zodat we in het echte leven minder fouten maken en meer van de geheimen van het universum kunnen ontdekken.

Technische samenvatting

Probleemstelling

In moderne experimenten voor kern- en deeltjesfysica, zoals zoektochten naar neutrinoloze dubbel-beta-verval ($0\nu\beta\beta$) en donkere materie, spelen hoogzuivere germanium (HPGe) detectoren een cruciale rol. Deze detectoren, vaak met puntcontacten (BEGe) of segmentatie, vereisen een uiterst nauwkeurige analyse van pulsvormen om achtergrondruis te onderdrukken en positiebepaling mogelijk te maken.

Bestaande simulatiepakketten voor HPGe-detectoren hebben echter beperkingen:

1. Ze nemen vaak de omgeving van de detector (zoals steunstructuren, koude vloeistoffen of afscherming) niet mee in de berekening, wat van invloed kan zijn op het elektrische veld.
2. Ze zijn vaak niet open-source of moeilijk uit te breiden voor nieuwe geometrieën en impureitsprofielen.
3. Ze zijn niet altijd geoptimaliseerd voor efficiënte 3D-berekeningen of parallelle verwerking.

Er is behoefte aan een flexibel, open-source hulpmiddel dat volledige 3D-simulaties van elektrische velden, ladingsdrift en geïnduceerde signalen kan uitvoeren, inclusief de invloed van de omgeving, om zo de productie en optimalisatie van detectoren (bijv. voor het LEGEND-experiment) te ondersteunen.

Methodologie

De auteurs introduceren SolidStateDetectors.jl, een open-source softwarepakket geschreven in de programmeertaal Julia. Julia werd gekozen vanwege de combinatie van gebruiksvriendelijke syntax en hoge prestaties, wat essentieel is voor complexe numerieke berekeningen.

De simulatieprocedure omvat de volgende stappen:

1. Definitie van Geometrie en Randvoorwaarden: De gebruiker definieert de detector, de omgeving en elektrische randvoorwaarden via gestructureerde tekstconfiguraties. Constructive Solid Geometry (CSG) wordt gebruikt om vormen te definiëren.
2. Berekening van Elektrische Velden:
   • Het elektrische potentiaalveld $\Phi(\mathbf{r})$ wordt numeriek opgelost via de Poisson-vergelijking (Gauss' wet) met een Successive Over-Relaxation (SOR) algoritme op een adaptief rooster.
   • Het rooster wordt automatisch verfijnd waar potentiaalgradiënten groot zijn.
   • Het pakket kan onontgonnen (undepleted) gebieden modelleren door de ladingsdichtheid daar op nul te stellen.
   • Gewichtspotentiaalvelden ($\Phi_w$) worden berekend voor elke elektrode om geïnduceerde signalen te bepalen (Shockley-Ramo theorema).
3. Ladingsdrift en Signaalvorming:
   • Ladingsdragers (elektronen en gaten) worden gedrift door het berekende elektrische veld volgens specifieke driftsnelheidsmodellen (afhankelijk van temperatuur).
   • De drift wordt stap voor stap berekend (standaard 1 ns stappen).
   • Geïnduceerde signalen op de elektrodes worden berekend als de som van de bijdragen van alle driftende ladingen, gebruikmakend van de gewichtspotentiaalvelden.
4. Omgevingsinvloed: Het pakket kan materialen rondom de detector (zoals koperen schalen of vloeibare argon) meenemen in de berekening, wat van invloed is op het elektrische veld, vooral nabij gepassiveerde oppervlakken.

Belangrijkste Bijdragen

• Volledige 3D-simulatie met symmetrie: Het pakket voert berekeningen uit in volledige 3D, maar kan ook gebruikmaken van symmetrieën (bijv. periodieke randvoorwaarden in de $\phi$-richting) om rekentijd te besparen.
• Omgevingsmodelleren: Het is mogelijk om de invloed van de fysieke omgeving (cryostaten, vloeibare argon, afscherming) op het elektrische veld en de detectie-eigenschappen te simuleren.
• Flexibiliteit en Uitbreidbaarheid: Als open-source pakket in Julia is het eenvoudig aan te passen voor verschillende halfgeleidermaterialen (niet alleen Ge, maar ook Si) en complexe geometrieën.
• Validatie: Het pakket is grondig gevalideerd door vergelijking met analytische oplossingen en experimentele data.

Resultaten

De auteurs hebben SolidStateDetectors.jl gevalideerd aan de hand van een n-type gesegmenteerde BEGe-detectoren:

1. Numerieke Validatie: Vergelijking met analytische oplossingen voor een oneindig lange coaxiale detector toonde een zeer kleine RMS-afwijking (tot $2.1 \cdot 10^{-6}$ V voor cilindrische coördinaten). De berekende capaciteiten kwamen overeen met theoretische waarden.
2. Omgevingsinvloed: Simulaties toonden aan dat de omgeving (bijv. een koperen schaal of onderdompeling in vloeibare argon) significante veranderingen in het elektrische potentiaalveld veroorzaakt, vooral nabij gepassiveerde oppervlakken. Dit is cruciaal voor het begrijpen van achtergronden in zeldzame gebeurtenis-experimenten.
3. Vergelijking met Data:
   • Oppervlakte-gebeurtenissen: Voor 81 keV fotonen (dicht bij het oppervlak) kwamen de gesimuleerde "super-pulsen" goed overeen met gemeten data.
   • Bulk-gebeurtenissen: Voor interacties in het bulk-materiaal (gebruikmakend van Compton-scattering data) toonde de simulatie een goede overeenkomst, inclusief de vorming van spiegelpulsen op niet-verzamelpende segmenten.
   • Pulsvormanalyse (A/E): De simulatie van de $A/E$-parameter (verhouding tussen stroomamplitude en energie) voor $^{208}$Tl en $^{212}$Bi gebeurtenissen bevestigde dat het pakket de scheiding tussen enkel- en meer-site gebeurtenissen nauwkeurig kan voorspellen. De overlevingsfracties voor single-site gebeurtenissen en de onderdrukking van multi-site gebeurtenissen kwamen binnen statistische onzekerheden overeen met de data.

Betekenis

SolidStateDetectors.jl is een essentiële tool voor de volgende generatie experimenten met HPGe-detectoren, zoals LEGEND.

• Optimalisatie: Het stelt onderzoekers in staat om de snijlijnen van kristallen te optimaliseren voor maximale massa en minimale bias-spanning, rekening houdend met onzuiverheidsprofielen.
• Achtergrondreductie: Door nauwkeurige simulaties van pulsvormen en achtergrondmechanismen (zoals oppervlakte-gebeurtenissen beïnvloed door de omgeving), kan de efficiëntie van achtergrondreductietechnieken worden verhoogd.
• Toekomstige Ontwikkeling: Het pakket is ontworpen om te groeien; versie 1.0 (gepland voor eind 2021) zal functies toevoegen zoals simulatie van ladingswolk-interacties (zelfafstoting en diffusie), betere visualisatie en GPU-ondersteuning.

Kortom, het pakket vult een belangrijke lacune in de simulatiewereld door een snel, flexibel en open-source alternatief te bieden dat de complexiteit van moderne halfgeleiderdetectoren in hun volledige omgeving kan modelleren.