🔬 physics

A four parameters modified Gaussian function describing shape of asymmetric alpha lines

Dit artikel presenteert een nieuwe vierparametrische analytische functie voor het nauwkeurig modelleren van asymmetrische alfa-pieken, die is geïntegreerd in de gevestigde ALF-software om de deconvolutie van tot tien overlappende pieken mogelijk te maken.

Oorspronkelijke auteurs: Jerzy Wojciech Mietelski

Gepubliceerd 2026-07-09

📖 4 min leestijd☕ Koffiepauze-leesvoer

CC BY 4.0

Oorspronkelijke auteurs: Jerzy Wojciech Mietelski

Oorspronkelijk artikel gelicentieerd onder CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/). ✨ Dit is een AI-gegenereerde uitleg van het onderstaande artikel. Het is niet geschreven of goedgekeurd door de auteurs. Raadpleeg het oorspronkelijke artikel voor technische nauwkeurigheid. Lees de volledige disclaimer

Stel je voor dat je luistert naar een koor van alfadeeltjes die hun weg zingen naar een detector toe. In een perfecte wereld zouden alle zangers exact dezelfde noot raken, wat een prachtige, symmetrische klokvorm creëert—een klassieke Gaussische vorm. Maar in de echte wereld is het een rommelig proces. Terwijl deze minuscule deeltjes van hun bron naar de detector reizen, botsen ze tegen dingen aan: een "dode laag" op het ingangsvenster van de detector en het bronmateriaal zelf. Deze botsingen stelen een beetje van hun energie, maar niet altijd een gelijke hoeveelheid.

Denk aan hardlopers in een race waarbij sommigen struikelen over een losse steen en anderen niet. De hardlopers die niet struikelen, komen exact op tijd aan, waardoor de hoogenergetische kant van de piek scherp en symmetrisch blijft. Maar de hardlopers die wel struikelen, komen te laat aan, waardoor de laagenergetische kant van de piek uitloopt in een lange, rommelige staart. Dit creëert een "asymmetrische" vorm die lijkt op een klokvorm die aan één kant door de modder is gesleept.

Decennialang worstelden wetenschappers om deze rommelige vorm met wiskunde te beschrijven. Meestal moesten ze een ingewikkelde gereedschapskist van vele verschillende formules gebruiken om de verschillende delen van de curve aan elkaar te plakken. Maar Jerzy Wojciech Mietelski, een onderzoeker van het Instituut voor Kernfysica in Polen, besloot een andere aanpak te proberen. Hij introduceerde een enkele, slimme wiskundige functie die fungeert als een "vormveranderaar".

Deze nieuwe functie is een gemodificeerde Gaussische functie, maar dan met een twist. Het heeft slechts vier instelbare knoppen (parameters) om aan te draaien:

Amplitude: Hoe hoog de piek is.
Positie: Waar de piek zich op de grafiek bevindt.
Vormparameter 1: Controleert hoe de curve zich aan de rechterkant gedraagt (de schone, Gaussische kant).
Vormparameter 2: Controleert de rommelige, lange staart aan de linkerkant.

Hier is de magische truc: Aan de rechterkant van de piek gedraagt de functie zich als een standaard, perfecte klokvorm. Maar aan de linkerkant transformeert het in een complexere formule die van nature die lange, exponentiële staart creëert die wordt veroorzaakt door het energieverlies. Het is alsof je een enkel stuk klei hebt dat aan de ene kant glad en rond kan zijn, maar aan de andere kant uitgerekt en taps toelopend, zonder dat je twee verschillende vormen aan elkaar hoeft te lijmen.

Mietelski heeft dit niet alleen verzonnen; hij heeft het aan het werk gezet in een computercode genaamd ALF, die sinds 1993 in zijn laboratorium draait. Het proces verloopt in twee stappen, zoals het afstemmen van een radio. Eerst kijkt de computer naar een enkele, heldere piek en bepaalt de exacte instellingen voor de vormparameters (de "modder" en de "schone" kanten). Zodra de computer weet hoe de "modder" er specifiek voor dat experiment uitziet, worden die instellingen vastgezet. Vervolgens gebruikt hij deze vaste vorm om een rommel van overlappende pieken—tot wel 10 tegelijk—te ontwarren door alleen de hoogte en positie aan te passen.

Het artikel meldt dat deze methode is getest over een periode van 30 jaar en meer dan 15.000 spectra succesvol heeft geanalyseerd. Het heeft zelfs strenge "bekwaamheidstests" doorstaan, waarbij het laboratorium moest bewijzen dat zij zaken correct konden meten ten opzien van andere experts. De auteurs merken op dat de resultaten zo betrouwbaar zijn dat de oppervlakte onder de curve (die aangeeft hoeveel radioactief materiaal er is) bijna perfect overeenkomt met het totale aantal deeltjes.

Verrassend genoeg was deze "alfa"-vormveranderaar niet alleen bedoeld voor alfadeeltjes. Het artikel vermeldt dat het team dezelfde wiskunde probeerde te gebruiken om bèta-stralingssignalen van een vloeibaar scintillati spectrometer te scheiden. In één test met 90Sr (Strontium-90) in evenwicht met 90Y (Yttrium-90), slaagde de methode erin de signalen te scheiden met een foutmarge van slechts 3%.

Dus, hoewel het artikel niet beweert dat dit een magische oplossing is voor elk probleem in de fysica, presenteert het een robuust, vier-parametervariabel instrument dat decennialang bewezen heeft te werken. Het verandert een rommelig, asymmetrisch probleem in een oplosbare puzzel, waardoor wetenschappers helder door de ruis van overlappende signalen kunnen kijken, of ze nu kijken naar alfadeeltjes van Antarcticaanse korstmossen of bètapartikels uit dierlijke botten.

Technische Samenvatting: Een vier-parametrische gemodificeerde Gaussische functie voor asymmetrische alfa-lijnen

Probleemstelling
In de alfaspectrometrie vertonen spectrale lijnen vaak asymmetrische vormen in plaats van ideale Gaussische verdelingen. Deze asymmetrie ontstaat primair door energiedegradatie wanneer alfadeeltjes de dode laag (ingangvenster) van de halfgeleiderdetector en het bronmateriaal passeren. De omvang van deze degradatie hangt af van de willekeurige emissiehoek van het deeltje ten opzichte van de detector. Als gevolg hiervan wordt de laagenergetische zijde van de piek breder en nadert deze een exponentieel verval, terwijl de hoogenergetische zijde Gaussisch blijft, slechts gemodificeerd door statistische fluctuaties.

Historisch gezien was er geen enkele, algemeen aanvaarde wiskundige formule om deze vorm te beschrijven. Praktische oplossingen vertrouiden doorgaans op multi-parameter functiesets die worden toegepast op specifieke secties van de spectrale lijn. Hoewel diverse softwarepakketten (bijv. ALFUN-A, MicroSAMPO, ALFA) bestaan om deze uitdagingen aan te pakken, was de specifieke analytische functie die gebruikt werd in de breed geaccepteerde ALF-computer code tot voor dit werk niet gedetailleerd in de literatuur.

Methodologie
Het artikel presenteert een originele analytische functie die is ontworpen om de asymmetrische alfa-piekvorm te beschrijven met slechts vier vrije parameters. Deze functie is geïmplementeerd in de ALF computer code, een op Fortran gebaseerd hulpmiddel ontwikkeld rond 1993 (gemodificeerd van een eerdere gamma-straal code, P.I.M.P.) voor de evaluatie van alfaspectra.

De kern van de methodologie is de volgende gemodificeerde Gaussische functie, $N(x)$ , waarbij $x$ het kanaalnummer vertegenwoordigt:

$N(x) = \alpha \cdot \exp\left( \frac{-(x-M)^2}{\beta \{|x-M| - (x-M)\} + \gamma} \right)$

Parameters:

$\alpha$ : Amplitude (aantal tellingen bij het maximum).
$M$ : Positie van het pieklokaal.
$\beta$ en $\gamma$ : Vormparameters (beide $>0$ ) specifiek voor de bron- en detectorgeometrie.

Gedrag van de functie:

Hoogenergetische zijde ( $x > M$ ): De functie convergeert naar een standaard Gaussische verdeling.
Laagenergetische zijde ( $x < M$ ): De functie wordt complexer en gaat bij aanzienlijke afstanden van het pieklokaal ( $x \ll M$ ) over naar een gewone exponentiële functie.
Rol van de parameters: $\gamma$ reguleert de vorm van de rechterzijde van de piek, terwijl $\beta$ en $\gamma$ gezamenlijk de linkerzijde reguleren.

Fittingprocedure (Deconvolutie):
De ALF-code maakt gebruik van een niet-lineaire functie-fittingprocedure om tot 10 overlappende alfa-pieken te deconvuleren. Het proces wordt uitgevoerd in twee afzonderlijke stappen:

Initiële vormbepaling: Een enkele, goed gedefinieerde spectrale lijn (of een fragment dat het maximum bevat) wordt gefit om de vier parameters ( $\alpha, M, \beta, \gamma$ ) te bepalen. Hiervoor wordt een standaard niet-lineaire kleinste-kwadratenmethode gebruikt via 10 iteraties van een lineaire benadering (Taylorreeks-expansie).
Piekdeconvolutie: Zodra $\beta$ en $\gamma$ zijn bepaald, worden deze als vaste constanten behandeld voor het gehele spectrum, aangezien zij afhankelijk zijn van de bronkwaliteit en de bron-detectorgeometrie in plaats van individuele pieken. De tweede stap omvat het fitten van de amplitudes ( $\alpha_i$ ) en posities ( $M_i$ ) voor alle overlappende pieken simultaan. Dit wordt bereikt door de som van functies uit te breiden naar een reeks en de incrementen van de vectorvariabelen op te lossen met behulp van een Cracovian-algoritme.

Het Cracovian-algoritme wordt benadrukt als een specifiek voordeel, omdat het de onzekerheden van de gefitte parameters berekent binnen dezelfde procedure die wordt gebruikt om de parameterwaarden te bepalen.

Belangrijkste bijdragen en resultaten

Analytische eenvoud: Het artikel introduceert een compacte, vier-parametrische functie die effectief de complexe asymmetrie van alfa-lijnen modelleert zonder dat complexe multi-segment fitting vereist is.
Validatie en duurzaamheid: De functie is sinds 1993 in routinegebruik. In de afgelopen 30 jaar is de functie toegepast op meer dan 15.000 spectra. Validatie werd uitgevoerd door de berekende piekoppervlakten te vergelijken met de som van de tellingen, waarbij alle tests positief werden doorstaan.
Deconvolutiecapaciteit: De methode slaagt erin om gedeeltelijk overlappende pieken te scheiden, zoals aangetoond met Th-fractie monsters uit dierlijke botten.
Onverwacht nut: De auteurs merken op dat de functie verrassend nuttig bleek voor het analyseren van bèta-radiatiestadia (specifiek $^{90}\text{Sr}$ / $^{90}\text{Y}$ ) verkregen met een vloeistofscintillatiemeter met een logaritmische kanaalschaal, waarbij resultaten binnen een foutmarge van 3% werden behaald.

Betekenis
Het artikel claimt zijn belang primair door de succesvolle, langdurige validatie van deze specifieke vier-parametrische functie binnen de ALF-software. Door de voorheen niet gedocumenteerde analytische basis van ALF te formaliseren, biedt dit werk een robuust, wiskundig gedefinieerd instrument voor de deconvolutie van overlappende alfa-pieken. De methode biedt een evenwicht tussen analytische strengheid en computationele efficiëntie, waardoor een nauwkeurige bepaling van piekoppervlakken mogelijk is in complexe spectra waar traditionele Gaussische fitting faalt. Het gebruik van het Cracovian-algoritme zorgt er bovendien voor dat de parameteronzekerheden rigoureus worden berekend naast de fit.

Technische Samenvatting: Een vier-parametrische gemodificeerde Gaussische functie voor asymmetrische alfa-lijnen

Meer zoals dit