Spektrometria gamma w praktyce: kalibracja energii, rozdzielczość i wydajność detektora

Streszczenie

Widmo gamma wygląda jak wykres z pikami, ale praktyczna spektrometria zaczyna się dopiero wtedy, gdy umiemy powiedzieć, co oznacza oś energii, jak szeroki jest pik, ile zliczeń należy do fotopiku, a ile do tła, i jaką część emitowanych fotonów detektor rzeczywiście rejestruje. Bez tych czterech kroków spektrometr jest tylko bardziej efektownym licznikiem impulsów.^1,2

Ten artykuł pokazuje typowy tok pracy: od surowego widma, przez kalibrację energii, ocenę rozdzielczości, wyznaczenie pola piku, korektę tła, krzywą wydajności i ograniczenia czasu martwego, aż do ostrożnej identyfikacji radionuklidu albo oszacowania aktywności.

Rozszerzenie tematu

Widmo nie jest jeszcze wynikiem

Wyobraźmy sobie ekran analizatora wielokanałowego. Na osi poziomej są kanały, na osi pionowej liczba zliczeń. Widać kilka pików i szerokie tło. Początkujący użytkownik chce od razu kliknąć pik i nazwać izotop. Metrolog powinien najpierw zapytać: czy kanały są skalibrowane na energię, czy pik jest rozdzielony od sąsiednich, jak odjęto tło, jaka jest wydajność detektora i czy wynik dotyczy czasu rzeczywistego czy żywego?

Spektrometria gamma jest silna dlatego, że fotony gamma wielu radionuklidów mają charakterystyczne energie. Nie oznacza to jednak, że każdy pik automatycznie identyfikuje izotop. Trzeba uwzględnić rozdzielczość detektora, tło Comptona, linie o podobnych energiach, piki sumacyjne, czas półtrwania, geometrię i kontekst próbki.

Dziunikowski i Kalita opisują spektrometr promieniowania gamma jako układ z detektorem, wzmacniaczem liniowym i analizatorem amplitudy impulsów; ćwiczenia obejmują cechowanie energii, rozdzielczość połówkową i wydajność detekcji.¹ Przewodnik KChRS pokazuje podobną logikę na przykładach dydaktycznych z sondą NaI(Tl) i analizatorem wielokanałowym.²

Co rejestruje analizator wielokanałowy

Analizator wielokanałowy (MCA) tworzy histogram amplitud impulsów. Impuls o małej amplitudzie trafia do niskiego kanału, impuls o większej amplitudzie do wyższego. Jeżeli tor jest liniowy, a detektor zachowuje informację energetyczną, po kalibracji kanał można przeliczyć na energię.

To jest punkt różniący spektrometrię od zwykłego zliczania. Licznik impulsów pyta: ile zdarzeń przekroczyło próg? Spektrometr pyta: jak rozkładają się amplitudy tych zdarzeń? Dzięki temu w widmie gamma można zobaczyć fotopiki, skraje Comptona, obszary tła i czasem cechy aparaturowe.

W dobrym opisie widma trzeba więc podać:

• typ detektora,
• liczbę kanałów,
• zakres energii,
• czas pomiaru,
• czas żywy, jeśli system go raportuje,
• geometrię próbki,
• kalibrację energii,
• sposób wyznaczania pól pików,
• sposób odjęcia tła.

Samo zdjęcie widma jest materiałem poglądowym, ale nie wystarcza jako wynik ilościowy.

Kalibracja energii: kanał staje się keV

Kalibracja energii polega na przypisaniu znanych energii gamma do położeń pików w kanałach analizatora. W najprostszym zakresie przyjmuje się zależność liniową:

E = a * k + b,

gdzie E jest energią, k numerem kanału, a nachyleniem, a b przesunięciem zera analizatora. Dziunikowski i Kalita w ćwiczeniach spektrometrycznych zalecają wyznaczenie zależności między energią fotonów i numerem kanału metodą regresji liniowej oraz określenie zera analizatora.¹

W praktyce potrzebne są co najmniej dwa punkty, ale trzy lub więcej punktów pozwala sprawdzić liniowość i wykryć błąd. Jeżeli jeden punkt kalibracyjny odstaje, nie należy go mechanicznie usuwać bez przyczyny. Trzeba sprawdzić, czy pik nie został źle odczytany, czy nie jest dubletem, czy tło nie przesunęło maksimum i czy układ nie dryfował.

Kalibracja energii powinna obejmować zakres, w którym analizujemy próbkę. Ekstrapolacja daleko poza użyte linie jest słabsza metrologicznie. Jeśli kalibrujemy układ tylko niskimi energiami, a potem interpretujemy pik powyżej 1 MeV, należy wyraźnie zaznaczyć ograniczenie.

Fotopik, skraj Comptona i tło ciągłe

Fotopik, nazywany też pikiem pełnej absorpcji, odpowiada zdarzeniom, w których cała energia fotonu gamma została zdeponowana w detektorze. To właśnie fotopik zwykle wykorzystuje się do identyfikacji energii i liczenia aktywności. Obok niego występuje tło komptonowskie, bo część fotonów rozprasza się i zostawia w detektorze tylko część energii.³

W scyntylatorze NaI(Tl) tło komptonowskie bywa szczególnie wyraźne, a piki są szerokie. W detektorze HPGe piki są węższe, lecz tło i efekty geometrii nadal istnieją. Nie wolno więc sprowadzać spektrometrii do listy energii z tabeli. Trzeba patrzeć na całe widmo.

Praktyczny błąd polega na zliczaniu wszystkiego „pod górką” bez modelu tła. Jeżeli pik leży na pochyłym tle, samo zsumowanie kanałów zawyży pole. Jeżeli pik jest słaby, sposób oszacowania tła może decydować o wyniku bardziej niż statystyka samego piku.

Pole piku i zliczenia netto

Pole piku to liczba zliczeń przypisana interesującej linii gamma po odjęciu tła. Najprostszy model używa zakresu zainteresowania (ROI) wokół piku i lokalnego tła oszacowanego z kanałów po lewej i prawej stronie. Bardziej zaawansowane metody dopasowują funkcję piku i funkcję tła.

Uproszczony rachunek wygląda tak:

N_net = N_ROI - N_tlo,

gdzie N_ROI jest sumą zliczeń w oknie piku, a N_tlo oszacowaną liczbą zliczeń tła w tym samym oknie. Jeśli tło mierzono osobno, trzeba przeskalować je do tego samego czasu pomiaru:

N_tlo,sk = N_tlo * t_s / t_b.

Przewodnik KChRS wielokrotnie prowadzi ćwiczenia przez odczyt liczby zliczeń pod fotopikami i odejmowanie tła, na przykład dla K-40, Cs-137 i próbek środowiskowych.² Dydaktycznie to bardzo dobra praktyka, bo student widzi, że „pik” nie jest pojedynczą liczbą, lecz obszarem widma wymagającym decyzji analitycznej.

Rozdzielczość energetyczna i FWHM

Rozdzielczość energetyczna mówi, jak wąski jest pik. W spektrometrii gamma najczęściej używa się szerokości połówkowej FWHM, czyli szerokości piku w połowie jego wysokości. Rozdzielczość można podawać bezwzględnie, na przykład w keV, albo względnie:

R = FWHM / E.

Dziunikowski i Kalita definiują energetyczną zdolność rozdzielczą przez szerokość połówkową piku i omawiają wkłady pochodzące z fluktuacji statystycznych, szumów aparatury i konstrukcji detektora.¹ W praktyce FWHM nie jest tylko parametrem katalogowym. Warto kontrolować go okresowo, bo pogorszenie rozdzielczości może wskazywać na problem z detektorem, elektroniką, zasilaniem, chłodzeniem albo warunkami pomiaru.

Dwa detektory mogą pokazać tę samą energię, ale zupełnie inną rozdzielczość. NaI(Tl) może być wystarczający do prostych ćwiczeń i wielu pomiarów terenowych, ale przy skomplikowanym widmie naturalnych szeregów promieniotwórczych szerokie piki ograniczają jednoznaczność. HPGe rozwiązuje część tych problemów, lecz wymaga bardziej wymagającego stanowiska.⁴

Wydajność detektora: ile fotonów trafia do fotopiku

Jeżeli z próbki emitowanych jest 1000 fotonów gamma, detektor nie zarejestruje wszystkich w fotopiku. Część nie trafi w detektor, część zostanie osłabiona w próbce lub naczyniu, część zdeponuje tylko część energii, a część zdarzeń nie przejdzie przez kryteria analizy. Dlatego potrzebna jest wydajność.

W spektrometrii gamma zwykle interesuje nas wydajność w fotopiku:

epsilon(E, geometria) = liczba zliczeń w fotopiku / liczba fotonów gamma wyemitowanych w kierunku całej przestrzeni.

To wielkość zależna od energii i geometrii. Nie jest jedną stałą detektora. Dla niskich energii istotne jest osłabienie w oknie, powietrzu, osłonie i próbce. Dla wyższych energii rośnie znaczenie ucieczki fotonów rozproszonych i rozmiaru detektora. Dziunikowski i Kalita rozróżniają wydajność całkowitą i wydajność w piku głównym, co jest kluczowe dla interpretacji widm NaI(Tl).¹

Krzywa wydajności powinna być wyznaczona dla tej geometrii, w której mierzymy próbki. Jeśli wzorzec leży na detektorze w małym naczyniu, a próbka jest dużym pojemnikiem, nie wolno bezrefleksyjnie przenosić wydajności.

Od pola piku do aktywności

Dla pojedynczej linii gamma uproszczony wzór na aktywność próbki można zapisać jako:

A = N_net / (t_live * epsilon(E) * I_gamma).

Jeżeli wynik ma być aktywnością właściwą:

a_m = N_net / (t_live * epsilon(E) * I_gamma * m).

Oznaczenia:

• N_net: zliczenia netto w fotopiku,
• t_live: czas żywy pomiaru,
• epsilon(E): wydajność w fotopiku dla energii i geometrii,
• I_gamma: prawdopodobieństwo emisji danej linii gamma na rozpad,
• m: masa próbki.

Ten wzór wygląda prosty, ale każdy składnik wymaga kontroli. Zliczenia netto zależą od tła i ROI. Czas żywy zależy od czasu martwego. Wydajność zależy od geometrii. I_gamma pochodzi z danych jądrowych. Masa próbki zależy od przygotowania i wilgotności. Dlatego wynik bez niepewności jest tylko połową informacji.

Czas martwy, czas żywy i pile-up

W spektrometrii gamma wysoka częstość zliczeń może zniekształcić widmo. Układ potrzebuje czasu na przetworzenie impulsu. Gdy zdarzeń jest za dużo, część z nich nie zostanie zarejestrowana, a część może nałożyć się na siebie. Wtedy rośnie czas martwy i pojawia się pile-up.⁵

W praktyce trzeba rozróżniać:

• real time: czas zegarowy trwania pomiaru,
• live time: czas efektywnie dostępny dla rejestracji impulsów,
• dead time: część czasu, w której tor nie przyjmuje nowych zdarzeń.

Jeżeli spektrometr raportuje czas żywy, do obliczeń aktywności zwykle używa się właśnie jego. Jeżeli system podaje tylko czas rzeczywisty, a czas martwy jest istotny, wynik może być zaniżony. Przy dużym pile-up problem jest głębszy: nie chodzi tylko o utratę zdarzeń, ale także o fałszywe impulsy o zawyżonej amplitudzie.

Identyfikacja radionuklidu: energia to nie wszystko

Pierwszy krok identyfikacji to porównanie energii piku z tablicami linii gamma. Nie jest to jednak dowód sam w sobie. Trzeba sprawdzić:

• czy widoczne są inne linie tego samego radionuklidu,
• czy intensywności linii są zgodne w granicach wydajności,
• czy energia nie pokrywa się z linią innego radionuklidu,
• czy pik nie jest dubletem,
• czy tło naturalne mogło dać podobny sygnał,
• czy kontekst próbki uzasadnia obecność danego izotopu.

Przykładowo K-40 ma znaną linię około 1461 keV, ale jej obecność w materiale budowlanym albo żywności nie jest sensacją; jest elementem naturalnej promieniotwórczości. Linie z szeregu uranowego i torowego wymagają ostrożności, bo w próbkach środowiskowych pojawia się problem równowagi promieniotwórczej, tła i nakładania pików.²

Dobry opis identyfikacji powinien mówić „przypisano pik o energii...” albo „wynik jest zgodny z obecnością...”, a nie udawać absolutnej pewności tam, gdzie rozdzielczość i tło jej nie dają.

Geometria, samopochłanianie i matryca próbki

Spektrometr mierzy nie tylko radionuklid. Mierzy radionuklid w konkretnej próbce, w konkretnym pojemniku, w konkretnej odległości od detektora. Gęstość, skład chemiczny i grubość próbki mogą zmienić liczbę fotonów, które opuszczą próbkę i dotrą do detektora.

KChRS wskazuje praktycznie, że problem samopochłaniania można ograniczyć, jeśli kalibracja wydajnościowa jest wykonana dla próbek w tej samej geometrii i o podobnych właściwościach osłabiania.² To jest zasada ważniejsza niż wiele późniejszych poprawek: wzorzec powinien możliwie przypominać próbkę.

Jeżeli matryca jest inna, potrzebne są poprawki, symulacja Monte Carlo albo wzorce matrycowe. W artykule ogólnym nie trzeba wchodzić w szczegółowe procedury. Wystarczy podkreślić, że aktywność bezwzględna z widma gamma wymaga geometrii i wydajności, a nie tylko pola piku.

Kontrola jakości widma

Przed interpretacją warto zrobić krótką kontrolę jakości:

• czy piki kalibracyjne są w oczekiwanych kanałach,
• czy FWHM piku kontrolnego nie pogorszyło się,
• czy tło nie wzrosło nietypowo,
• czy czas martwy jest akceptowalny,
• czy nie widać nasycenia albo pile-up,
• czy widmo tła i widmo próbki mają porównywalne ustawienia,
• czy zapisano masę, geometrię, czas i ustawienia toru.

Ta kontrola jest szczególnie ważna w dydaktyce, bo studenci często ufają automatycznemu oprogramowaniu. Program może poprawnie narysować widmo, a jednocześnie wynik analityczny może być słaby z powodu złej geometrii, zbyt krótkiego tła albo przesuniętej kalibracji.

Minimalny opis wyniku spektrometrii gamma

W artykule albo raporcie wynik powinien mieć formę bardziej kompletną niż „wykryto Cs-137”. Minimalny opis powinien zawierać:

• detektor i tor spektrometryczny,
• kalibrację energii,
• kontrolę rozdzielczości,
• czas rzeczywisty i żywy,
• sposób odjęcia tła,
• zakresy ROI lub metodę dopasowania piku,
• wydajność dla energii i geometrii,
• masę lub objętość próbki,
• wynik aktywności z niepewnością,
• datę odniesienia, jeśli aktywność przeliczano w czasie.

Jeżeli artykuł ma tylko charakter jakościowy, część tych elementów można pominąć, ale trzeba wtedy jasno powiedzieć, że widmo pokazuje identyfikację orientacyjną, a nie pełny pomiar aktywności.

Niepewność rozszerzona w spektrometrii gamma — składowe budżetu

ISO/IEC 17025 i GUM (Guide to the Expression of Uncertainty in Measurement) wymagają, by wynik pomiaru był podawany wraz z niepewnością rozszerzoną U = k × u_c (k = 2 dla ~95% poziomu ufności). Budżet niepewności dla aktywności z spektrometrii gamma obejmuje:

Klasa A (statystyczna, z powtórzeń lub Poisson):

• u_stat(N_net): √N_net / N_net = 1/√N_net; dla 10 000 zliczeń: u_rel = 1%;
• u_stat(N_tło): propagacja z tła;

Klasa B (z certyfikatów, modeli, danych fizykalnych):

• u(ε): niepewność kalibracji wydajności (1–5% dla dobrze skalibrowanego HPGe przy Marinelli);
• u(I_gamma): niepewność intensywności gamma z tablic ENSDF/NUDAT (0,1–2%);
• u(A_wzorzec): niepewność certyfikatu wzorca kalibracyjnego (1–3%);
• u(t_live): niepewność czasu żywego (zazwyczaj <0,1%);
• u(m): niepewność masy próbki (waga, 0,1–0,5%);
• u(samopochłanianie): przy kalibracji matrycowo dopasowanej często <2%, bez korekcji — 10–50%;
• u(TCS): przy bliskiej geometrii dla nuklidów kaskadowych (0–5%);
• u(biblioteka): interpretacja niezgodności między bibliotekami nuklidów (<0,5%).

Składowe sumuje się kwadratowo (przy braku korelacji):
u_c/A = √(u_stat² + u_ε² + u_I² + u_A_wzorzec² + u_m² + u_SA² + u_TCS²)

Typowe łączne u_c/A dla pełnie opisanego pomiaru rutynowego: 3–8% (k=1). Przy słabej aktywności (N_net = 100 zliczeń): u_stat = 10% dominuje.

Akredytowane laboratoria PCA muszą mieć formalnie zatwierdzony budżet niepewności i walidować go przez pomiar CRM co roku. Wynik pomiaru bez niepewności nie spełnia wymagań ISO 17025 i nie może być używany jako dowód regulacyjny (np. w postępowaniu o dopuszczenie żywności lub gleby do użytku).

Wzorce kalibracyjne — tabela radionuklidów

Wybór radionuklidów do kalibracji zależy od zakresu energii i charakteru próbek. Poniżej zestawienie często używanych wzorców w polskich laboratoriach środowiskowych:

Radionuklid	E_gamma (keV)	T₁/₂	Intensywność gamma	Uwagi
Am-241	59,54	432,2 lat	35,9%	Niska energia; ostrożnie z absorpcją okienka
Co-57	122,1 / 136,5	272 dni	85,6% / 10,7%	Bliskie linie — test rozdzielczości
Cs-137/Ba-137m	661,7	30,2 roku	85,1%	Najpopularniejszy wzorzec środowiskowy
Co-60	1173,2 / 1332,5	5,27 roku	99,9% / 100%	Dwie linie — test kalibracji liniowości
Eu-152	121,8 / 344,3 / 778,9 / 1408,0	13,5 roku	28,5% / 26,5% / 12,9% / 21,0%	Wielopunktowy wzorzec kalibracyjny
Ba-133	53,2 / 79,6 / 81,0 / 160,6	10,5 roku	2,2% / 2,6% / 34,1% / 0,6%	Niskie energie; dobry do weryfikacji przy <100 keV
Mn-54	834,8	312 dni	99,98%	Prosta linia; aktywacyjna
Zn-65	1115,5	244 dni	50,6%	Wysoka energia; aktywacyjna
Y-88	898,0 / 1836,1	106,6 dni	93,7% / 99,2%	Wysoka energia; aktywacyjna
Bi-207	569,7 / 1063,7	32,9 roku	97,7% / 74,6%	Alternatywa dla Co-60 bez koincydencji sum

Certyfikowane wzorce do kalibracji wydajności dostarcza NIST (SRM), Eckert & Ziegler (Nuclitec), PerkinElmer, Amersham. Kalibracja serwisowa (factory calibration) wzorca jest ważna przez 2–5 lat; po tym czasie aktywność oblicza się przez przeliczenie z daty certyfikacji.

Szybkie sprawdzenie widma przed interpretacją — lista kontrolna

Przed przystąpieniem do analizy ilościowej zaleca się przebiec przez kilka punktów weryfikacyjnych:

1. Czy kalibracja energii jest aktualna? Sprawdź przez pik wzorcowy (np. Cs-137 661,7 keV lub Co-60 1332,5 keV) — czy centrum mieści się w ±1 kanale oczekiwanej pozycji?
2. Czy FWHM jest poprawne? Dla HPGe: FWHM przy 1332 keV (Co-60) powinno wynosić 1,7–2,5 keV; powyżej 3 keV — problem z detektorem lub elektroniką.
3. Czy czas martwy < 5%? Powyżej 10% konieczna ocena pile-up.
4. Czy tło zostało zmierzone w tej samej geometrii? Nie przenoś tła z innego pojemnika lub dnia pomiaru bez normalizacji.
5. Czy odjąłeś tło? Nie porównuj „aktywności" próbki bez odjęcia tła.
6. Czy piki tła (K-40 1461 keV, Bi-214 609 keV, Pb-208 2614 keV) mają oczekiwane natężenie? Anomalie wskazują na problem z geometrią lub na skażenie próbki.
7. Czy wzorzec kalibracyjny miał tę samą matrycę i geometrię? Inaczej: opisz korekcję samopochłaniania.

Ta lista jest krótka i wykonuje się ją w 5 minut — ale eliminuje 80% typowych błędów analitycznych początkujących.

Sumowanie koincydencyjne — efekt pominięty przez amatorów

Sumowanie koincydencyjne (coincidence summing, TCS — true coincidence summing) jest jednym z najbardziej niedocenianych źródeł błędów systematycznych w precyzyjnej spektrometrii gamma. Efekt polega na tym, że radionuklid emituje dwa (lub więcej) fotonów kaskadowo w czasie krótszym niż czas rozdzielczy detektora. Detektor rejestruje je jako jeden impuls o energii będącej sumą obu.

Przykład Co-60: kaskada 1173,2 keV + 1332,5 keV. Jeżeli obydwa fotony zostaną zarejestrowane przez detektor, trafiają do kanału odpowiadającego sumie energii 2505,7 keV. W fotopiku 1173 keV pojawia się deficyt zliczeń, a w pikach sumy pojawia się nadmiar.

Efekt TCS jest silny przy bliskiej geometrii (mały dystans źródło–detektor, wysoka wydajność geometryczna). Przy zwiększeniu odległości maleje prawdopodobieństwo jednoczesnego trafienia obu fotonów. Korekcja TCS jest opisana w normie ISO 11929 i IAEA TECDOC-1361. Oprogramowanie LabSOCS i Efficiency Tracer (Mirion) oblicza korekcje TCS automatycznie.

Nuklidy szczególnie podatne na TCS: Co-60, Eu-152, Eu-154, Ra-226 (przez córki), Tb-160, Sc-46. Dla Cs-137 (jednofotony 661,7 keV) TCS jest negligible. Dlatego w kalibracji objętościowej przy bliskiej geometrii stosuje się Cs-137 lub Ba-133, a nie Co-60, chyba że korekcja TCS jest uwzględniona.

Szeregi promieniotwórcze i równowaga w próbce

W próbkach środowiskowych (gleba, osad, granit, materiały budowlane) aktywność pochodzi głównie z radionuklidów naturalnych: K-40, U-238 i jego pochodne (szereg Ra), Th-232 i jego pochodne. Warunkiem poprawnej interpretacji jest stan równowagi lub nieróżnowagi między nuklidami szeregu.

W glebie nienaruszonej przez procesy chemiczne lub geologiczne przez długi czas zakłada się równowagę wiekową: aktywność Ra-226 = aktywność Pb-210 = aktywność Po-210 itd. To upraszcza obliczenia. Jednak gleba skażona (np. po deszczu kwasowym, po przemyciu przez wodę, po hutnicy) może mieć zaburzoną równowagę.

Typowe linie gamma z szeregu U-238 mierzone przez HPGe:

• Ra-226 przez Pb-214 i Bi-214: 295,2 keV (Pb-214, 19,2%); 352,0 keV (Pb-214, 37,6%); 609,3 keV (Bi-214, 45,5%); 1120,3 keV (Bi-214, 15%); 1764,5 keV (Bi-214, 15,3%);
• Th-232 przez Ac-228 i Tl-208: 338,3 keV (Ac-228, 11,3%); 911,2 keV (Ac-228, 25,8%); 2614,5 keV (Tl-208, 99,8%).

Kluczowa zasada: Ra-226 mierzy się pośrednio przez córki (Pb-214, Bi-214) — po zahermetyzowaniu próbki przez min. 4 tygodnie, aby Rn-222 i córki osiągnęły równowagę wewnątrz próbki. Pominięcie tego kroku prowadzi do zaniżenia aktywności Ra-226 o kilkadziesiąt procent, jeśli Rn-222 zdążył wylecieć z próbki przed pomiarem.

Identyfikacja i unikanie pułapek interpretacyjnych

Kilka typowych pułapek przy interpretacji widm gamma:

Piki od naturalnego tła: laboratoryjne tło zawiera Cs-137 (z opadów atomowych, obecny w ścianach betonowych), K-40 (w cegle, betonie), Ra-226 i Th-232 (granit). Jeśli próbka jest małej aktywności, tło może dominować. Zawsze mierz tło w identycznym pojemniku (pusty pojemnik Marinelli, wypełniony wodą dejonizowaną) przez czas równy co najmniej czasowi pomiaru próbki.

Piki pary elektron-pozyton: powyżej 1022 keV mogą pojawiać się piki ucieczki (single escape peak = E − 511 keV, double escape peak = E − 1022 keV). Bi-214 przy 1764,5 keV daje single escape przy 1253,5 keV. Nieostrożna identyfikacja traktuje taki pik jako niezależną linię gamma.

Piki charakterystyczne ołowiu: osłona ołowiana fluorescuje gdy rejestruje promieniowanie Röntgena — Pb Kα przy 72–75 keV, Kβ przy 84–87 keV. Te linie są obecne w każdym tle z osłoną ołowianą. Nie są sygnałem Pb w próbce.

Suma przypadkowa (accidental summing): przy dużej częstości zliczeń dwa niezwiązane fotony mogą trafić do detektora w tym samym oknie czasowym i zostać zarejestrowane jako jedno zdarzenie o energii sumy. Różni się od TCS tym, że nie wymaga kaskady — jest efektem statystycznym. Redukcja: obniżenie częstości zliczeń (rozcieńczenie próbki, większa odległość od detektora).

Kalibracja wydajności — procedura i wzorce objętościowe

Kalibracja wydajności wyznacza ε(E) dla konkretnej geometrii i detektora. Procedura typowego laboratorium środowiskowego:

1. Wybór geometrii: decydujemy, w jakich pojemnikach mierzymy próbki — Marinelli 1L lub 0,5L, cylinder 100 mL, cylinder 500 mL. Każda geometria wymaga osobnej kalibracji.

2. Przygotowanie wzorca: certyfikowany wzorzec wielopunktowy (np. Eu-152 + Am-241 + Cs-137 zmieszane w matrycy zbliżonej do próbek) w pojemniku identycznym z pomiarowym. NIST SRM 4276 (Gamma-ray spectrometry calibration standard) jest mieszaniną ~15 radionuklidów w objętościowej formie.

3. Pomiar wzorca: minimum 10 000 zliczeń w kluczowych fotopikach (Cs-137, Co-60) dla dobrej statystyki. Czas martwy < 5%.

4. Obliczenie ε(E) dla każdej linii: ε = N_net / (A × I_gamma × t_live), gdzie A jest aktywnością wzorca przeliczoną na datę pomiaru.

5. Dopasowanie krzywej ε(E): logarytmiczno-logarytmiczne dopasowanie (log ε vs. log E) lub wielomian. Oprogramowanie (Genie 2000, GammaVision) wykonuje to automatycznie.

6. Walidacja: pomiar niezależnego CRM (certified reference material) z zewnątrz serii kalibracyjnej. Wynik powinien zgadzać się z wartością certyfikowaną w zakresie 2σ.

Kalibracja musi być powtarzana po: wymianie detektora, zmiana geometrii naczyń, naprawie lub modyfikacji elektroniki, serwisie systemu chłodzenia.

Dopasowanie kształtu piku — metody analityczne

Pole piku N_net wyznacza się przez: (a) metodę ROI z lokalnym tłem, lub (b) dopasowanie funkcji kształtu piku do danych. Metoda (b) jest dokładniejsza przy słabych pikach lub nakładaniu się linii.

Kształt fotopiku w HPGe jest dobrze opisany przez funkcję Gaussa z assymetrycznym prawostronnym eksponentem (step function):

P(E) = G × exp[−(E−E₀)²/(2σ²)] + A × [1 − erf((E−E₀)/(√2 σ))]

Gdzie G to amplituda Gaussa, A to amplituda schodka Comptona (step), E₀ to centrum piku, σ = FWHM/2.355. Parametry E₀, σ, G, A są dopasowywane do danych metodą chi-kwadrat lub maksymalnego prawdopodobieństwa.

W NaI(Tl) kształt piku jest mniej regularny (niesymetryczny) z powodu niekompletnego zbierania ładunku w PM i efektów scyntylacyjnych. Prostszy model z lewostronnym eksponentem (tail parameter) daje lepsze dopasowanie niż sam Gauss.

Oprogramowanie komercyjne (Genie 2000 SAMPO, GammaVision, Maestro): automatycznie przeszukuje widmo, proponuje identyfikacje z bazy nuklidów, oblicza pola pików i aktywności. Użytkownik musi weryfikować wyniki, szczególnie przy słabych liniach, nakładaniu i linii z tłem środowiskowym.

Naturalne tło i metody jego redukcji

W każdym laboratorium spektrometrii gamma istnieje tło: promieniowanie kosmiczne, naturalne radionuklidy w otoczeniu (K-40, Ra-226, Rn-222 i córki, Th-232), materiały budowlane (ściany, podłogi), elementy detektora i elektroniki. Tło ogranicza MDA (minimal detectable activity).

Redukcja pasywna: osłony ołowiane (5–10 cm stary ołów bez Ra-210, „low background lead"), osłony z OFHC miedzi wewnątrz ołowiu (pochłania fluorescencję Pb Kα przy 72–75 keV), komory betonowe lub granitowe (redukują kosmiczne muony).

Redukcja aktywna: antykoincydencja BGO lub plastikowy scyntylator wokół detektora HPGe. Zdarzenie w zewnętrznym scyntylatorze koincydujące z impulsem HPGe w oknie 100 ns jest odrzucane. Redukuje tło Comptona o czynnik 3–10 i tło kosmiczne. Rozbudowane systemy ultra-niskiego tła (underground labs: Modane, Gran Sasso, Sobieszów) mogą osiągnąć MDA ~0,1 mBq/kg dla Cs-137 przy czasie pomiaru 24–72 h.

Wskaźnik jakości tła F = (A_tło/MDA): im mniejsze tło, tym niższy MDA dla tej samej statystyki. W typowym laboratorium naziemnym MDA dla Cs-137 w 1 kg gleby, pojemnik Marinelli 1L, czas 24 h: ~0,5–2 Bq/kg. W laboratorium ultra-niskiego tła: ~0,05–0,1 Bq/kg.

Oprogramowanie do analizy widm gamma

Kompletny przepływ pracy w laboratorium środowiskowym:

Genie 2000 (Mirion Technologies / Canberra): standard branżowy, zintegrowany z dedykowanymi przedwzmacniaczami DSA i detektorami Mirion. Zawiera SAMPO (analiza piku), biblioteki nuklidów, moduły do analizy bezwzorcowej (LabSOCS) i raporty aktywności.

GammaVision (ORTEC/AMETEK): alternatywny standard, kompatybilny z detektorami ORTEC. Zaawansowane algorytmy identyfikacji, obsługa bibliotek ENSDF, automatyczne korekcje czasu martwego i geometrii.

Maestro (ORTEC): prostszy interfejs do dydaktyki i szybkich pomiarów.

InterSpec / FitzPeaks / Spectrum Manager: narzędzia otwarto-źródłowe lub darmowe, przydatne do nauki.

Krytyczne dla jakości: każde oprogramowanie identyfikacyjne ma bibliotekę nuklidów z energiami i intensywnościami. Różne wydania bibliotek (ENSDF, NUDAT, Table of Isotopes) mogą różnić się wartościami. Zawsze sprawdzaj wersję biblioteki i datę ostatniej aktualizacji.

InterSpec jest darmowym narzędziem open-source (LLNL), dostępnym jako aplikacja desktopowa i web. Obsługuje widma z popularnych formatów (CHN, SPE, SPC, N42), ma silnik identyfikacji nuklidów i narzędzia kalibracji. Przydatny dla laboratoriów akademickich i do celów edukacyjnych. Dokumentacja i kod źródłowy są publicznie dostępne.

Spectrum Manager (dostępny przez IAEA dla krajów członkowskich): narzędzie szkoleniowe IAEA używane w kursach bezpieczeństwa jądrowego. Zawiera bazy danych nuklidów i ćwiczenia identyfikacyjne. Polska uczestniczy w kursach IAEA regularnie przez NCBJ Świerk i IFJ PAN.

Ważna uwaga dotycząca bibliotek: ENSDF (Evaluated Nuclear Structure Data File) jest zarządzany przez BNL (Brookhaven National Laboratory). Aktualizacje są publikowane nieregularnie; wersja 2020+ zawiera znacząco poprawione wartości intensywności niektórych linii (np. Bi-214, Pa-234m). Przy pomiarach do celów regulacyjnych lub publikacji naukowej należy podać wersję i datę źródła danych jądrowych.

Polska praktyka spektrometrii gamma

Polskie laboratoria środowiskowe prowadzą rutynową spektrometrię gamma w ramach systemu monitoringu Państwowej Agencji Atomistyki (sieć DOZ) i CLOR. CLOR ma kilkanaście detektorów HPGe o efektywności względnej 20–40% pracujących przy Marinelli 1L, skalibrowanych na wzorcach CRM z JRC Geel. Wyniki przekazywane są do EURDEP.

Po wypadku w Czarnobylu (1986) Polska uruchomiła stałe stacje monitoringu w kilku lokalizacjach (Belsk, Kraków, Białystok, Wrocław). Po Fukushimie (2011) rozszerzono sieć o automatyczne stacje z detektorami NaI(Tl) raportującymi online.

Laboratorium ultra-niskiego tła IFJ PAN Sobieszów (330 m n.p.m., Karkonosze) oferuje usługi pomiarów środowiskowych o ekstremalnie niskim MDA. Mierzone tu próbki obejmują gleby z Arktyki, osady jeziorne, próbki radiochemiczne Pu i Am, a nawet próbki z meteorytów i trynitytów.

Trzy przykłady rachunkowe

Przykład 1 — kalibracja energii

Dane: w widmie obserwujemy trzy piki na kanałach 298, 662, 1173. Znane linie gamma: Cs-137 (661,7 keV), Co-60 (1173,2 i 1332,5 keV). Przyjmujemy dwa punkty: kanał 662 = 661,7 keV, kanał 1173 = 1173,2 keV.

Nachylenie a = (1173,2 − 661,7) / (1173 − 662) = 511,5 / 511 = 1,00097 keV/kanał. Przesunięcie b = 661,7 − 1,00097 × 662 = 661,7 − 662,6 = −0,9 keV.

Linia kalibracyjna: E = 1,001 × k − 0,9 (keV). Sprawdzamy trzecią linię (kanał 298): E = 1,001 × 298 − 0,9 = 297,4 − 0,9 = 296,5 keV. Z tablic: Ba-133 ma linię przy 302,9 keV (21,6%). Różnica 6,4 keV — możliwy niedopasowany wzorzec lub konieczność nieliniowej kalibracji.

Wniosek: liniowa kalibracja z dwóch punktów jest przybliżeniem. Przy 3 lub więcej punktach i dopasowaniu leastsquares uzyskujemy residua — wskazujące na liniowość toru elektronicznego.

Przykład 2 — wydajność i aktywność Cs-137

Dane: próbka gleby 0,8 kg, pojemnik Marinelli 1L, HPGe. Wzorzec Cs-137 w Marinelli: aktywność certyfikowana A_wzorzec = 150 Bq, zliczenia netto N_wzorzec = 12 000 cps × 300 s = 3 600 000 zliczeń (po korekcie t_live = 298 s). Wydajność: ε(662 keV) = N_wzorzec / (A_wzorzec × I_gamma × t_live) = 3 600 000 / (150 × 0,851 × 298) = 3 600 000 / 38 016 = 94,7 cps/Bq/sonde.

Przeliczamy na bezwymiarowe ε = ε_cps_total / total_emissions = 94,7 × 1 / (total_emissions per second). Prościej: ε = 94,7 zliczeń/Bq·s × 1/s = 94,7 × 10⁻³ / s = 94,7 cps/kBq. Dla próbki: N_net = 1200 zliczeń w 300 s t_live. A = N_net / (ε × I_gamma × t_live) = 1200 / (0,09474 × 0,851 × 298) = 1200 / 24,02 = 50 Bq. Aktywność właściwa: a_m = 50/0,8 = 62,5 Bq/kg.

Niepewność statystyczna: u_stat = √1200/1200 = 2,9%. Łączna (z wydajnością ~3%, I_gamma ~0,1%, masą ~0,5%): u_total ≈ 4,3%.

Wynik: a_m = (62,5 ± 2,7) Bq/kg Cs-137 (k=1).

Przykład 3 — MDA dla K-40

MDA (minimum detectable activity) wyznacza się z wzoru Currie'go dostosowanego do spektrometrii gamma:

MDA = 2,71 + 4,65 × √B / (ε × I_gamma × t_live × m)

Gdzie B jest tłem w ROI fotopiku K-40 (1461 keV). Dane: t_live = 3600 s, ε(1461 keV) = 0,028 (Marinelli 1L, HPGe 30%), I_gamma(K-40) = 10,7%, m = 1,0 kg. Tło w ROI 30 kanałów przy 1461 keV: B = 400 zliczeń.

MDA = (2,71 + 4,65 × √400) / (0,028 × 0,107 × 3600 × 1,0) = (2,71 + 4,65 × 20) / (10,79) = (2,71 + 93) / 10,79 = 95,71 / 10,79 = 8,9 Bq/kg.

Rzeczywiste aktywności K-40 w glebie (~400 Bq/kg) i materiałach budowlanych (~200–600 Bq/kg) są znacznie powyżej MDA. Przy pomiarze wody (K-40 ≈ 50 mBq/L) ta sama geometria jest już niewystarczająca i potrzebny jest długi czas pomiaru lub prekoncentracja.

Pytania otwarte

1. Dlaczego kalibracja energii z dwóch punktów może być niewystarczająca? Jakie źródła nieliniowości mogą występować w torze spektrometrycznym (ADC, wzmacniacz, detektor)?

2. Wzorzec kalibracyjny zawiera Cs-137 (t₁/₂ = 30,2 lat) certyfikowany 3 lata temu z aktywnością 1000 Bq. Oblicz aktywność w dniu kalibracji. Czy 3-letni wiek certyfikacji jest akceptowalny, czy wymaga nowego wzorca?

3. Dlaczego metoda ROI z lokalnym tłem może zawyżać pole piku, gdy pik leży na pochyłym tle Comptona? Jakie modele tła są bardziej adekwatne?

4. Wydajność ε(E) spada monotonicznie dla HPGe w zakresie 100–2000 keV. Dlaczego spada zarówno dla niskich energii (< 200 keV) jak i wyższych? Jakie procesy fizyczne powodują spadek w każdym zakresie?

5. Pole fotopiku piku Co-60 1173 keV i 1332 keV powinno być proporcjonalne do intensywności (99,9% vs. 100%). Jeśli zmierzone proporcje różnią się o 5%, co może to oznaczać: błąd kalibracji, sumowanie, czy inny efekt?

6. W próbce pojemnika Marinelli zamierzona geometria to 1L wody. Próbka jest faktycznie glebą sypką o gęstości 1,4 g/cm³ wypełniającą 85% objętości. Jak to wpływa na aktywność obliczoną z kalibracji wodnej? Zawyżona czy zaniżona?

7. Czas martwy wynosi 15%. Jak wpływa to na obliczoną aktywność, jeśli używamy t_live zamiast t_real? Kiedy wysoki czas martwy może prowadzić do błędów systetycznych (pile-up, efekty dynamiczne)?

8. Laboratorium porównuje aktywność Cs-137 z próbki zmierzonej dwoma metodami: HPGe z Marinelli i scyntylator NaI(Tl) z pojemnikiem cylindrycznym 500 mL. Oba wyniki powinny być zgodne. Jakie czynniki systematyczne mogą powodować rozbieżność i jak je zidentyfikować?

Podsumowanie dydaktyczne

1. Spektrometria gamma jest sekwencją decyzji metrologicznych: kalibracja energii → rozdzielczość → odjęcie tła → pole piku → wydajność → aktywność. Każdy krok wprowadza składową niepewności, którą należy dokumentować.

2. Kalibracja energii musi obejmować zakres energii próbek. Liniowa interpolacja z dwóch punktów jest przybliżeniem — minimum 3 punkty z kontrolą residuów. Dryfty elektroniki (temperatura, napięcie) wymagają kontroli w czasie.

3. FWHM jest miarą rozdzielczości detektora w konkretnych warunkach pracy. Pogorszenie FWHM wskazuje na problem z detektorem lub elektroniką. Regularna kontrola wzorcem (np. Co-57 dla HPGe) jest standardem.

4. Pole piku N_net zależy od wyboru ROI i modelu tła. Lokalny model tła jest akceptowalny dla izolowanych pików; przy nakładaniu linii konieczne jest dopasowanie kształtu funkcji.

5. Wydajność ε(E) jest zależna od energii i geometrii. Kalibracja musi być wykonana dla tej samej geometrii co próbki. Wzorce matrycowe eliminują konieczność osobnej korekty samopochłaniania.

6. Czas żywy t_live musi być używany do obliczania aktywności, nie czas rzeczywisty. Przy czasie martwym > 5% konieczna jest ocena ewentualnego pile-up i deformacji widma.

7. Identyfikacja radionuklidu wymaga potwierdzenia co najmniej 2 linii z poprawnym stosunkiem intensywności, kontekstem próbki i wykluczeniem interferencji. Energia pojedynczego piku nie jest dowodem.

8. Polska sieć spektrometrii gamma (CLOR, IFJ PAN, PAA) opiera się na certyfikowanych wzorcach JRC Geel i akredytacji PCA. Laboratorium ultra-niskiego tła Sobieszów jest polskim centrum dla oznaczeń o najniższych MDA, ważnym dla badań archeometrii, skażeń tła i próbek arktycznych. Polska uczestniczy regularnie w rundach porównań IAEA ALMERA dla matryc środowiskowych (gleba, osad, biota) i uzyskuje wyniki w zakresie dopuszczalnym przez |z-score| ≤ 2 dla radionuklidów Cs-137, K-40, Ra-226, Pb-210 i Am-241. Wyniki są publikowane w raportach IAEA ALMERA i stanowią dowód kompetencji przy audytach PCA. Biorąc udział w PT (proficiency testing) laboratoria mają prawo do stosowania tych wyników jako poparcie akredytacji bez oddzielnej walidacji metody dla tych samych matryc.

Dodatkowe materiały multimedialne

Gotowe uzupełnienia do pracy z widmem: wizualizacja widma gamma, model 3D detektora HPGe i model 3D sondy NaI(Tl). Razem pozwalają porównać wygląd piku z fizyczną geometrią detektora.

Najkrótsze podsumowanie: spektrometria gamma jest sekwencją decyzji metrologicznych. Najpierw trzeba skalibrować oś energii, potem ocenić rozdzielczość, odjąć tło, policzyć pole piku, zastosować wydajność i dopiero wtedy mówić o radionuklidzie lub aktywności.

Powiązane kalkulatory i narzędzia

• Kalkulator: Kalibracja gamma — Dwupunktowa kalibracja energii, przeliczenie kanału na energię i rozdzielczość FWHM.
• Kalkulator: Aktywność gamma — Aktywność z liczby zliczeń w fotopiku po korekcie tła, wydajności, emisji i odzysku.

Ćwiczenia praktyczne

Pierwsze ćwiczenie: student dostaje tabelę kanałów i zliczeń dla syntetycznego widma z dwoma pikami. Ma wyznaczyć maksimum każdego piku, zbudować prostą kalibracyjną E = a * k + b na podstawie jednego znanego piku i dodatkowego punktu kontrolnego, a następnie przeliczyć kanały na energie.

Drugie ćwiczenie: na syntetycznym piku Gaussa obliczyć FWHM w kanałach i w keV, a następnie rozdzielczość względną R = FWHM / E. Potem powtórzyć rachunek dla szerszego piku NaI(Tl) i węższego piku HPGe, aby zobaczyć praktyczny sens rozdzielczości.

Trzecie ćwiczenie: dla wybranego ROI policzyć N_ROI, oszacować lokalne tło z kanałów bocznych i wyznaczyć N_net. Następnie zmienić szerokość ROI i sprawdzić, jak wybór okna wpływa na wynik.

Czwarte ćwiczenie: użyć uproszczonego wzoru A = N_net / (t_live * epsilon * I_gamma) dla danych syntetycznych. Następnie zmienić wydajność o 10% i sprawdzić, czy różnica jest większa niż niepewność statystyczna zliczeń.

Piąte ćwiczenie audytowe: przygotować opis wyniku dla widma materiału budowlanego zawierającego naturalne radionuklidy. Student ma wypisać, które piki traktuje jako jakościowe wskazanie obecności, a które nadają się do obliczenia aktywności, oraz jakie dane kalibracyjne byłyby potrzebne.

Przejdź do ćwiczenia interaktywnego

Powiązane artykuły

• Scyntylatory, fotopowielacze i półprzewodnikowe detektory promieniowania
• NaI(Tl), Ge(Li), HPGe: dlaczego różne detektory dają różne widma gamma
• Kalibracja i niepewność pomiaru w laboratorium jądrowym
• Statystyka zliczeń promieniotwórczych: Poisson, Gauss, chi-kwadrat i błędy aparaturowe
• Naturalna promieniotwórczość żywności i materiałów budowlanych: od widma do Bq/kg