NIM w pracowni studenckiej: jak czytać schemat blokowy stanowiska spektrometrycznego

Streszczenie

Schemat blokowy stanowiska spektrometrycznego wygląda prosto: detektor, zasilacz, wzmacniacz, analizator, komputer. Dopiero po chwili widać, że każdy blok odpowiada za inną część informacji i każdy może wprowadzić błąd.

NIM jest dobrym narzędziem do nauki czytania takich schematów, bo rozdziela funkcje na moduły. Student może przejść palcem po kablu od detektora do analizatora i zrozumieć, co dzieje się z impulsem na każdym etapie. To umiejętność ważniejsza niż zapamiętanie numeru konkretnego modułu.¹

Krótkie chassis NIM z modułami pomiarowymi: licznikiem/czasomierzem, analizatorem okna, wzmacniaczem i zasilaczem wysokiego napięcia. Fizyczny widok aparatury pomaga czytać schemat blokowy: każdy prostokąt na rysunku odpowiada realnemu modułowi, kablowi, wejściu, wyjściu i decyzji pomiarowej.

Rozszerzenie tematu

Pierwsza zasada czytania schematu brzmi: rozpoznaj, gdzie jest informacja o energii, a gdzie informacja o czasie. W spektrometrii gamma energia jest zakodowana głównie w amplitudzie impulsu po kształtowaniu. Informacja czasowa może być mniej ważna albo używana tylko do odrzucania pile-up i synchronizacji.

Druga zasada: oddziel sygnały analogowe od logicznych. Sygnał analogowy niesie amplitudę i kształt; nie wolno go traktować jak zwykłego "tak/nie". Sygnał logiczny po dyskryminatorze mówi już tylko, że impuls przeszedł próg. To rozróżnienie tłumaczy, dlaczego do widma potrzebujemy analizatora amplitudy, a do zliczania wystarczy skaler.

Trzecia zasada: znajdź zasilanie. Zasilacz wysokiego napięcia dla PMT lub detektora gazowego nie jest tym samym co zasilanie modułów NIM. Napięcie detektora wpływa na czułość i tryb pracy, a napięcia w binie zasilają elektronikę. W schemacie trzeba zaznaczyć oba obiegi.

Czwarta zasada: sprawdź impedancję i terminację. Jeśli moduł wysyła szybki impuls NIM kablem 50 Ω, odbiornik powinien być zgodny z tym standardem. Brak terminacji może podwoić amplitudę widzianą na oscyloskopie albo wygenerować odbicia. To częsty powód, dla którego dwa stanowiska "z tym samym detektorem" dają inne obserwacje.

Piąta zasada: szukaj miejsc decyzji. Próg dyskryminatora, okno SCA, bramka czasowa i ustawienie live time są decyzjami analitycznymi. Nie są neutralnymi elementami aparatury. W sprawozdaniu należy opisać, jak je ustawiono i dlaczego.

W praktyce dobry opis stanowiska powinien zawierać schemat blokowy, typ detektora, napięcie, kolejność modułów, rodzaj sygnałów, terminacje, nastawy wzmacniacza, progi, czas pomiaru, tło i sposób kalibracji. Dopiero taki opis pozwala komuś innemu zrozumieć wynik.

Od konceptu do rzeczywistości: co każdy moduł NIM naprawdę robi

Schemat blokowy to abstrakcja. Zrozumienie stanowiska wymaga wiedzy, co dzieje się fizycznie wewnątrz każdego bloku. W poniższym opisie prześledzimy tor sygnału od momentu, gdy promieniowanie jonizujące wchodzi do detektora, aż do zapisu kanału w histogramie widma.

Detektor — konwersja promieniowania na ładunek elektryczny. W detektorze NaI(Tl) z PMT: promieniowanie gamma powoduje jonizację w krysztale NaI:Tl, która daje fotony scyntylacyjne (emisja w zakresie UV-blue, ~415 nm, wydajność ~38 fotonów/keV). Fotony trafiają na fotokatodę PMT, gdzie część z nich ($\eta \approx 25\%$) wytwarza fotoelektrony. Kaskada dynod wzmacnia ładunek z typowym wzmocnieniem $M = 10^6$. Na anodzie PMT pojawia się impuls prądowy o czasie trwania ~1-2 μs, proporcjonalny do energii pierwotnego fotonu gamma.

W detektorze HPGe: promieniowanie gamma jonizuje atomy kryształu germanowego, tworząc pary elektron-dziura. Pole elektryczne, wymuszone polaryzacją $+2000$ V do $+5000$ V (napięcie biasowe), zbiega nośniki do elektrod. Czas zbioru ładunku $\sim 100{-}300$ ns. Ładunek $Q = E/\epsilon$, gdzie $\epsilon \approx 2{,}96$ eV/para dla Ge, więc dla 1 MeV: $Q = 10^6 / 2{,}96 \approx 3{,}4 \times 10^5$ par, $Q = 3{,}4 \times 10^5 \times 1{,}6 \times 10^{-19} \approx 54$ fC.

Przedwzmacniacz ładunkowy (CSP — Charge Sensitive Preamplifier) — konwertuje ładunek $Q$ na napięcie wyjściowe. Układ całkujący (transimpedancyjny) z pojemnością sprzężenia zwrotnego $C_f$ daje napięcie $V_{out} = Q/C_f$. Dla typowego $C_f = 1$ pF: $V_{out} = 54 \times 10^{-15} / 10^{-12} = 54$ mV. Kształt impulsu: szybkie narastanie (~zbieranie ładunku, 100-300 ns), powolny zanik (stała czasu rozładowania $RC_f \sim 50{-}200$ μs). Przedwzmacniacz jest montowany jak najbliżej detektora — zwykle bezpośrednio przy kriostatie HPGe lub w obudowie detektora scyntylacyjnego. Długi kabel między detektorem a przedwzmacniaczem to błąd projektowy, bo pojemność kabla $C_{cable}$ sumuje się z pojemnością wejściową i bezpośrednio pogarsza szum ($ENC \propto \sqrt{C_{total}}$).

Wzmacniacz kształtujący (Shaping Amplifier) — filtruje i wzmacnia sygnał z CSP. Usuwa długi ogon impulsu CSP przez różniczkowanie (czas CR), a następnie wygładza sygnał przez całkowanie (filtr RC albo Gaussowski). Wynikiem jest impuls Gaussowski o szerokości proporcjonalnej do stałej czasu kształtowania $\tau_p$. Dla spektrometrii HPGe typowe $\tau_p = 2{-}10$ μs. Amplituda impulsu po kształtowaniu jest proporcjonalna do oryginalnego ładunku $Q$, a więc do energii promieniowania. Zakres wzmocnienia regulowany jest potencjometrem lub przełącznikiem (typowo 1-3000×). Poprawne ustawienie wzmocnienia jest krytyczne: zbyt małe → sygnały nie wypełniają zakresu ADC i gubi się precyzja; zbyt duże → szczyty widma są obcinane.

Dyskryminator (SLD — Single Level Discriminator lub CFD — Constant Fraction Discriminator) — generuje cyfrowy impuls NIM Logic za każdym razem, gdy amplituda impulsu przekracza zadany próg. SLD jest prostszy, ale podatny na błąd czasowy (ang. walk): przy różnych amplitudach próg jest przekraczany w różnym momencie, co wprowadza jitter ~1-5 ns. CFD eliminuje ten efekt przez wewnętrzne porównanie z ułamkowym progiem względnym — jitter spada do ~100-300 ps. W praktyce dydaktycznej do prostej spektrometrii gamma SLD wystarczy. Do pomiaru koincydencji gamma-gamma lub gamma-neutron, gdzie liczy się precyzja czasowa, wymagany jest CFD.

Analizator amplitudy (MCA — Multi-Channel Analyzer lub ADC+histogramizer) — konwertuje amplitudę impulsu na numer kanału (ang. digitization). Wewnętrzny ADC (najczęściej Wilkinsona lub Successive Approximation) próbkuje szczyt impulsu po wstępnym rozciąganiu (peak stretching), przetwarzając amplitudę na liczbę. Wynik (numer kanału) inkrementuje licznik w odpowiedniej komórce histogramu widma. Typowe MCA mają 512, 1024, 2048, 4096 lub 8192 kanałów. Czas martwy ADC (podczas konwersji, ~1-10 μs) jest kompensowany przez korektę live time — mierzony przez wewnętrzny zegar aktywny tylko przy braku konwersji. Analizator może też filtrować sygnały przez okno energetyczne (SCA mode: zlicza tylko sygnały w zadanym przedziale amplitud) lub bramkę logiczną (Gate: akceptuje impulsy tylko przy aktywnym zewnętrznym sygnale bramkowania).

Skaler (Scaler) — zlicza impulsy logiczne NIM w zadanym oknie czasowym. Nie interesuje go amplituda — tylko fakt przekroczenia progu przez sygnał na wejściu. Używany do mierzenia aktywności próbek, weryfikacji martwego czasu aparatury, kalibracji geometrii detektora lub prostego monitoringu poziomu promieniowania. Wiele scalerów ma wejście inhibit (ang. gate), które zatrzymuje zliczanie na czas martwego czasu MCA, zapewniając spójność danych live-time.

Czytanie realnego schematu blokowego: przykład krok po kroku

Wyobraźmy sobie schemat stanowiska widmowego HPGe z koincydencją. Elementy:

[Źródło Gamma] → [HPGe + kriostat LN₂] → [CSP (wbudowany)] → [Wzmacniacz 2016A]
    → [ADC/MCA] (widmo energii)
    → [CFD] → [Moduł koincydencji] → [Gate → ADC]
[Scyntylator NaI] → [PMT] → [CSP] → [Wzmacniacz] → [CFD] → [Moduł koincydencji]
[HV supply] → [PMT]
[HV supply] → [HPGe]
[NIM Bin] → [zasilanie wszystkich modułów]

Czytanie krok po kroku:

Znajdź oba detektory (HPGe i NaI) — to są dwa niezależne "sensory" mierzące to samo zdarzenie z różnych stron lub mierzące koincydentne fotony.
Znajdź oba zasilacze HV — jeden dla HPGe (np. +3000 V), drugi dla PMT NaI (np. +1000 V). Czy są w binie? Czy to dedykowane moduły HV?
Prześledź tor HPGe: CSP → wzmacniacz → do ADC (amplituda, widmo energii). Odpowiedź na pytanie: ile czasu trwa od zdarzenia do wpisu do histogramu? Suma: $t_{collection}$ (~200 ns) + $t_{shaping}$ (~10 μs) + $t_{ADC}$ (~5 μs) = ~15 μs martwego czasu.
Prześledź tor koincydencji: z obydwu wzmacniaczy sygnał idzie do CFD → wyjścia logiczne → moduł koincydencji → sygnał Gate → blokuje (lub odblokowuje) bramkę ADC.
Zidentyfikuj bramkę ADC: przyjmuje sygnał energetyczny TYLKO gdy koincydencja jest aktywna. To oznacza: do histogramu trafiają wyłącznie zdarzenia, w których oba detektory wykryły foton w oknie czasowym ~100 ns.
Oceń, co może pójść źle: czy okno koincydencji jest właściwie ustawione? Czy progi CFD są odpowiednie dla energii fotonu? Czy live time jest poprawnie mierzony przy włączonej koincydencji?

Ten rodzaj analizy — prześledź sygnał, zidentyfikuj decyzje, oceń błędy — jest sednem "czytania schematu blokowego".

Typy schematów blokowych: od prostego do złożonego

Tor prosty (energia): Detektor → CSP → Shaping Amp → MCA. Przykład: pomiar widma Cs-137 z detektorem NaI. Zakres wiedzy: 1 zmienna (energia), 1 histogram.

Tor z dyskryminatorem: Detektor → CSP → Shaping Amp → SCA/Discriminator → Scaler. Przykład: pomiar aktywności próbki w wybranym oknie energetycznym. Tu informacja o energii jest używana do filtrowania, ale nie do budowy widma.

Tor koincydencyjny (energia + czas): 2× (Detektor → CSP → Shaping Amp → CFD) → Moduł koincydencji → Gate → MCA (jeden lub oba). Przykład: spektrometria pozytronów w materii skondensowanej (PALS — Positron Annihilation Lifetime Spectroscopy), koincydencja 1275 keV (narodziny pozytronu) i 511 keV (anihilacja). Wyniki: widmo czasowe (TAC/MCA) zamiast widma energii.

Tor anty-koincydencyjny (veto): Detektor główny → CSP → Shaping Amp → ADC; jednocześnie detektor weta (Compton suppressor lub detektor tła) → Discriminator → Veto → blokada ADC. Każde zdarzenie w detektorze głównym jest odrzucane, jeśli jednocześnie aktywny jest sygnał weta. Efekt: znaczące zmniejszenie tła Comtona w widmie, podwyższenie stosunku szczytu do tła (P/B).

Tor PHA + PSD (discriminacja kształtu impulsu): W detektorach organicznych (ciekłe scyntylatory, BC-501A) kształt impulsu różni się dla neutronu i fotonu gamma. Moduł PSD (Pulse Shape Discriminator) analizuje stosunek obszarów pod wolnym i szybkim składnikiem impulsu, wydając sygnał identyfikacji cząstki. Schemat: Detektor → CSP/FA → Shaping Amp + PSD → 2 gałęzie: [ADC widma dla neutronów] i [ADC widma dla fotonów].

Najczęstsze błędy studentów przy analizie schematów NIM

Praca z aparaturą NIM ujawnia systematyczne błędy koncepcyjne:

Błąd 1: Mylenie "zasilania detektora" z "zasilaniem modułów" — student widzi jeden zasilacz w binie i myśli, że zasila również detektor. Tymczasem HV supply to oddzielny moduł (albo zewnętrzne urządzenie), często o zupełnie innych wymaganiach. Pytanie diagnostyczne: "Jaka jest polaryzacja detektora i od czego ona zależy?"

Błąd 2: Ignorowanie martwego czasu — student wykonuje pomiar, nie notując live time, i interpretuje wynik jako liczbę zliczeń w czasie pomiarowym, choć faktyczny live time był 20% krótszy z powodu dużej intensywności źródła. Pytanie diagnostyczne: "Jaki był wskaźnik wejściowy impulsów i jaki był live time fraction?"

Błąd 3: Nieprawidłowa terminacja — student podłącza wyjście modułu logicznego do wysokoimpedancyjnego wejścia innego modułu bez terminacji 50 Ω, obserwując podwójne lub zwielokrotnione impulsy na oscyloskopie. Pytanie diagnostyczne: "Czy mam terminację na każdym wyjściu szybkiego sygnału NIM?"

Błąd 4: Zbyt mała lub zbyt duża stała kształtowania — przy zbyt małej $\tau$ wzmacniacz nie zbiera pełnego ładunku (incomplete charge collection), szczyt energetyczny w widmie przesuwa się w lewo i się poszerza. Przy zbyt dużej $\tau$ pile-up zdarzeń nakładających się sumuje ładunki dwóch fotonów w jeden "fałszywy" impuls. Pytanie diagnostyczne: "Jaka jest spodziewana aktywność źródła i jaki wynikający z tego wskaźnik pile-up?"

Błąd 5: Pominięcie kalibracji energetycznej — student buduje widmo, ale zapomina o kalibracji: przypisaniu osi X (kanałów) do skali energii (keV) przez pomiar co najmniej dwóch znanych pików. Pytanie diagnostyczne: "Jakich wzorców energetycznych użyłem do kalibracji i w jakim zakresie energii jest ona ważna?"

Błąd 6: Zaniedbanie tła — wynik pomiaru aktywności to zawsze "sygnał + tło". Student zapomina zmierzyć tło (np. ze ślepą próbką lub pustą komorą), przez co raportuje zawyżoną aktywność. Pytanie diagnostyczne: "Jaka jest moja procedura korekcji tła?"

Checklista przed i po pomiarze NIM

Przed pomiarem:

[ ] Czy zasilacz HV jest wyłączony przy podłączaniu detektora?
[ ] Czy napięcie HV jest ustawione prawidłowo (z karty detektora lub instrukcji)?
[ ] Czy wzmacniacz jest rozgrzany (~30 min po włączeniu binu)?
[ ] Czy wzmocnienie i stała kształtowania są właściwe dla planowanego zakresu energii?
[ ] Czy terminacje 50 Ω są założone na wszystkich szybkich kablach BNC?
[ ] Czy MCA jest w trybie PHA (Pulse Height Analysis), a czas pomiaru ustawiony?
[ ] Czy rozdzielczość (FWHM) widma referencyj jest zgodna z danymi katalogowymi detektora?

Po pomiarze:

[ ] Czy zanotowano live time i real time?
[ ] Czy zanotowano warunki pomiaru (napięcie, wzmocnienie, czas)?
[ ] Czy widmo ma piki w oczekiwanych kanałach/energiach?
[ ] Czy wyniki są porównywalne z poprzednimi pomiarami na tym samym stanowisku?
[ ] Czy zostały zmierzone tło i dane do kalibracji energetycznej?
[ ] Czy HV jest wyłączone po zakończeniu (ważne dla żywotności PMT)?

Porównanie typowych wartości parametrów dla NaI i HPGe

Parametr	NaI:Tl + PMT	HPGe
Napięcie detektora	600-1500 V	2000-5000 V
Czas narastania impulsu CSP	~1 μs	~100-300 ns
Optymalne $\tau_p$ kształtowania	0,5-2 μs	2-10 μs
FWHM przy 662 keV	~45-60 keV (7-9%)	<1,8 keV (<0,3%)
Wymaganie chłodzenia	nie (temp. pokojowa)	LN₂ (77 K) lub termoelektryczne
Typowy zakres energii	50-3000 keV	20-10 000 keV
Koszt (orientacyjnie)	1-5 k€	5-50 k€
Zastosowanie dydaktyczne	podstawowe	zaawansowane

Porównanie to wyjaśnia, dlaczego pracownia studencka zwykle zaczyna od NaI:Tl, a eksperyment z HPGe jest traktowany jako zaawansowane ćwiczenie. Nie chodzi tylko o cenę — HPGe wymaga ciągłego uzupełniania ciekłego azotu lub sprawnego systemu termoelektrycznego i jest wrażliwe na przypadkowe odwrócenie polaryzacji lub nadmierny prąd upływu przy nieodpowiednim temperaturze.

Znaczenie dokumentacji stanowiska w metrologii jądrowej

W zaawansowanych pomiarach jądrowych (np. w programach safeguards MAEA, w jądrowej medycynie nuklearnej lub w badaniach środowiskowych) każde stanowisko NIM jest traktowane jako narzędzie pomiarowe wymagające kalibracji, walidacji i dokumentacji zgodnie z normami ISO/IEC 17025 lub IEC 62755.

Dokumentacja obejmuje:

Świadectwo kalibracji detektora — wyznaczona przez producenta lub laboratorium wzorcujące efektywność bezwzględna $\varepsilon(E)$ (liczba wykrytych impulsów / liczba fotonów emitowanych przez źródło w daną stronę) jako funkcja energii
Wykaz modułów i wersji oprogramowania — każda zmiana modułu lub wersji firmware może zmienić charakterystykę toru
Logi kalibracji energetycznej — data, źródło wzorcowe, wyznaczone kanały pików, równanie kalibracji
Wyniki testów wewnętrznych — test rozdzielczości FWHM, test liniowości energetycznej, test martwego czasu

Nawet w laboratorium dydaktycznym wdrożenie prostej dokumentacji (zeszyt laboratoryjny z powyższymi informacjami) jest cenną praktyką, uczącą kultury metrologicznej, która jest bezpośrednio przenoszalna do środowiska profesjonalnego.

Ewolucja: od papierowych schematów blokowych do cyfrowych systemów akwizycji

Klasyczne schematy blokowe NIM, rysowane ołówkiem na papierze milimetrowym lub tablicy, są dziś zastępowane przez narzędzia CAD (Altium Designer, KiCad) lub specjalistyczne oprogramowanie dokumentacji eksperymentów (np. ELN — Electronic Lab Notebooks). Jednak strukturalna logika schematu blokowego pozostaje niezmieniona: sensor → kondycjonowanie sygnału → konwersja (A/D) → analiza.

W nowoczesnych systemach cyfrowego przetwarzania sygnałów (DPP — Digital Pulse Processing) moduły FPGA zastępują wiele analogowych bloków NIM. Zamiast osobnego wzmacniacza kształtującego, SCA i ADC, FPGA realizuje wszystkie te funkcje cyfrowo: próbkuje bezpośrednio wyjście CSP z częstotliwością 100-500 MHz i wykonuje kształtowanie numeryczne (filtr trapezoidal, Gaussian lub CUSP). Schemat blokowy staje się wtedy: Detektor → CSP → ADC 14-bit → FPGA (filtry + discriminatory + histogramizer) → komputer.

Dla studenta to zmiana techniczna, ale nie koncepcyjna: pytania o impedancję wejściową ADC, o cyfrowy odpowiednik terminacji (anti-aliasing filter), o "programowe okna" energetyczne i o mierzenie martwego czasu w FPGA — wszystko to są bezpośrednie analogie problemów klasycznego toru NIM. Zrozumienie analogowego NIM jest więc koniecznym warunkiem rozumienia systemów cyfrowych.

Kalibracja energetyczna: od kanałów do keV

Widmo w MCA jest zbiorem par (kanał, liczba zliczeń). Oś X to numer kanału (0 do N-1, gdzie N = liczba kanałów MCA). Aby przypisać energię fotonu do kanału, potrzebna jest kalibracja energetyczna — czyli funkcja mapująca numer kanału na energię w keV lub MeV.

Liniowa kalibracja energetyczna zakłada, że zależność kanał–energia jest liniowa:

$$E = a \cdot K + b$$

gdzie $K$ jest numerem kanału, $a$ jest nachyleniem (keV/kanał, zakres np. 0,5-2 keV/kanał), a $b$ jest przesunięciem (keV, idealne 0, w praktyce kilka keV). Do wyznaczenia dwóch parametrów $a$ i $b$ potrzebne są co najmniej dwa punkty wzorcowe — dwa dobrze znane szczyty w widmie źródła wzorcowego.

Typowe źródła wzorcowe w laboratorium dydaktycznym:

Izotop	Energia fotonu [keV]	T₁/₂
Cs-137	661,66	30,17 lat
Co-60	1173,2 + 1332,5	5,27 lat
Na-22	511,0 (anihilacja) + 1274,5	2,60 lat
Eu-152	121,8; 344,3; 1408,0; wiele	13,5 lat
Am-241	59,5	432 lat
Ba-133	81,0; 356,0	10,5 lat

Źródło Eu-152 jest szczególnie cenne, bo dostarcza kilkanaście pików wzorcowych w szerokim zakresie energii 40-1500 keV, pozwalając na weryfikację liniowości kalibracji. Nieliniowości ADC Wilkinsona rzędu 0,05-0,1% kanału widoczne są dopiero przy precyzyjnej wielopunktowej kalibracji.

Procedura kalibracji krok po kroku:

Ustaw wzmocnienie wzmacniacza tak, by najwyższy spodziewany szczyt mieścił się w ~80% zakresu kanałów MCA (np. dla Co-60: 1332,5 keV w kanale ~3500 z 4096).
Zmierz widmo źródła wzorcowego przez czas wystarczający do zebrania statystyki ~10 000 zliczeń w szczytowym kanale każdego piku (czas zależny od aktywności źródła i geometrii).
Dopasuj krzywą Gaussa do każdego piku, wyznaczając centroid kanału $K_i$ z precyzją poniżej 1 kanału (oprogramowanie GammaVision, Genie 2000, Maestro, lub własny skrypt Python).
Wykonaj regresję liniową $(K_i, E_i)$ metodą najmniejszych kwadratów. Wynik: $a$ i $b$ z niepewnościami $\sigma_a$ i $\sigma_b$.
Sprawdź residua: dla dobrego ADC Wilkinsona residua powinny być <0,1% zakresu energii.
Zanotuj kalibrację w logu laboratoryjnym z datą i warunkami pomiaru.

Kalibracja musi być powtarzana po każdej zmianie napięcia HV, wzmocnienia wzmacniacza, temperatury detektora lub po długiej przerwie w pracy (wzmacniacze analogowe mają dryft termiczny wzmocnienia ~0,01-0,1%/°C).

Obliczanie efektywności detektora

Samo widmo energetyczne mówi, ile zliczeń zarejestrowano w danym piku. By odpowiedzieć na pytanie "ile atomów rozpada się na sekundę w próbce?" (czyli jaka jest aktywność), potrzebna jest efektywność detektora $\varepsilon(E)$.

Efektywność jest funkcją energii i zależy od:

Geometrii (odległość próbka–detektor, kąt bryłowy $\Omega$)
Własności materiału detektora (gęstość, Z, grubość)
Pochłaniania promieniowania między próbką a detektorem (obudowa, powietrze)
Martwego czasu i efektów pile-up

Definicja:

$$\varepsilon(E) = \frac{N_{detected}(E)}{N_{emitted}(E)}$$

Dla wyznaczenia efektywności używa się źródła wzorcowego o dokładnie znanych aktywnościach i rozgałęzieniach gamma (ang. branching ratios). Zmierzone zliczenia $N_{peak}$ w piku energii $E$:

$$N_{peak} = A \cdot t \cdot I_\gamma \cdot \varepsilon(E)$$

gdzie $A$ jest aktywnością źródła w chwili pomiaru [Bq], $t$ jest live time [s], $I_\gamma$ jest intensywnością gamma (frakcją rozpadów emitujących foton energii $E$). Stąd:

$$\varepsilon(E) = \frac{N_{peak}}{A \cdot t \cdot I_\gamma}$$

Wyznaczona empirycznie krzywa efektywności $\varepsilon(E)$ jest charakterystyczna dla danego stanowiska (geometria + detektor + warunki). Zmiana odległości detektora o 1 mm zmienia efektywność o kilka procent (dla geometrii bliskiej) — dlatego precyzyjne pomiary wymagają mechanicznego stojaka zapewniającego powtarzalną geometrię.

Analiza piku i wyznaczanie aktywności

Szczyt (pik) w widmie gamma ma kształt zbliżony do Gaussa (bo szum detektora ma rozkład normalny). Procedura analizy piku:

Wycięcie okna — wybierz zakres kanałów zawierający pik i tło lokalne (zwykle kilka szerokości $\sigma$ po obu stronach centroidu).
Odejmowanie tła — tło pod pikiem wyznacza się przez interpolację liniową (uproszczone) lub dopasowanie wielomianowe w pasmach bocznych. Metoda całkowania brutto: $N_{net} = N_{gross} - B$, gdzie $B$ jest estymowanym tłem pod pikiem.
Dopasowanie Gaussa — model: $y(K) = A_0 \cdot \exp\left[-(K-K_0)^2/(2\sigma^2)\right] + B(K)$, gdzie $A_0$ jest amplitudą piku, $K_0$ centroidem kanałowym, $\sigma$ odchyleniem standardowym. FWHM = $2\sqrt{2\ln 2} \cdot \sigma \approx 2{,}355 \sigma$.
Całkowanie — aktywność netto $N_{net} = A_0 \cdot \sigma \cdot \sqrt{2\pi}$ (całka Gaussa). Metoda numeryczna (sumowanie kanałów po odjęciu tła) jest dokładniejsza, gdy pik ma nieliniowe ogony (np. wkład Comptona, odpowiedź niepełna).
Propagacja niepewności — statystyczna niepewność $\sigma_{N_{net}} = \sqrt{N_{gross} + B_{est}}$ (Poisson). Do niepewności musimy dodać: niepewność kalibracji efektywności (~1-5%), niepewność aktywności źródła wzorcowego (~2-5% z certyfikatu), niepewność live time (~0,1%), niepewność geometrii (~0,5-2%).

Wynik: aktywność próbki $A_{sample} = N_{net} / (t \cdot I_\gamma \cdot \varepsilon(E))$ z całkowitą niepewnością złożoną $\sigma_{rel,total} \approx \sqrt{\sigma_{stat}^2 + \sigma_{\varepsilon}^2 + \sigma_{A_{ref}}^2 + ...}$.

Pomiar czasu połowicznego zaniku

Jednym z klasycznych ćwiczeń laboratoryjnych z NIM jest pomiar czasu połowicznego zaniku krótkożyciowego izotopu (np. Ag-108 z T₁/₂ = 2,37 min, Ba-137m z T₁/₂ = 2,55 min, lub Rn-222 z T₁/₂ = 3,82 dnia).

Procedura dla Ag-108 (aktywacja neutronowa folii srebrnej):

Aktywuj foilę srebrną w źródle neutronowym (np. Cf-252 lub w reaktorze) przez kilka minut.
Wyjmij foilę, umieść w pobliżu detektora GM lub scyntylatora i uruchom MCA lub skaler z bramką czasową 10 s.
Rejestruj zliczenia co 10 s przez 20-30 minut.
Narysuj wykres $\ln(N(t))$ vs $t$. Prawo zaniku: $N(t) = N_0 \cdot e^{-\lambda t}$, więc $\ln N(t) = \ln N_0 - \lambda t$.
Dopasuj linię prostą metodą najmniejszych kwadratów — nachylenie daje $\lambda$, a $T_{1/2} = \ln 2 / \lambda$.

Dokładność wyznaczenia T₁/₂ zależy od:

Liczby punktów czasowych i statystyki w każdym punkcie ($\sigma_{stat} \propto 1/\sqrt{N_i}$)
Korekcji na czas martwy (szczególnie ważne na początku, gdy aktywność jest duża)
Korekcji na tło (odejmowanie zliczen przy pustej komorze)
Efektu pile-up w pierwszych pomiarach (można zmniejszyć przez atenuację próbki lub zwiększenie odległości)

To ćwiczenie uczy jednocześnie: metody naukowej pomiarów fizycznych, analizy statystycznej, korekcji systematycznych i znaczenia opisu niepewności.

Sprzęt analogowy vs cyfrowy: co student traci i co zyskuje

W ostatnich latach wiele laboratoriów migruje z klasycznych systemów analogowych NIM do cyfrowych systemów DPP (Digital Pulse Processing), takich jak CAEN DT5720 (250 MSPS, 12-bit), Mesytec MDPP-16, lub Red Pitaya z własnym firmware jądrowym.

Co student zyskuje przy cyfrowym DPP:

Brak oddzielnych modułów wzmacniacza, SCA, MCA — wszystko w jednym module
Bezpośredni zapis surowych danych (listmode) ze znacznikami czasu — możliwość późniejszej analizy off-line
Łatwiejsze ustawianie parametrów przez GUI komputera zamiast potencjometrów
Brak kalibracji ruchomych elementów mechanicznych

Co student traci przy cyfrowym DPP:

Widoczność sygnału na każdym etapie — klasyczny NIM pozwala wpinać oscyloskop między każdy moduł i "zobaczyć" impuls przed i po kształtowaniu. W systemie DPP sygnał analogowy istnieje tylko przez kilkadziesiąt ns przed ADC i jest niewidoczny bez specjalnego narzędzia diagnostycznego.
Intuicja dla rozkładu na komponenty — klasyczny schemat blokowy wymaga zrozumienia, co każdy blok robi z sygnałem. W "czarnej skrzynce" DPP ta wiedza jest ukryta w parametrach filtru trapezoidalnego.
Doświadczenie diagnostyczne — gdy coś nie działa w klasycznym NIM, student systematycznie sprawdza każdy moduł oscyloskopem. W DPP debugowanie wymaga wiedzy o wewnętrznym oprogramowaniu FPGA.

Dlatego wiele uczelni zachowuje klasyczne ćwiczenia z NIM analogowym jako wstęp do nauki aparatury jądrowej, nawet jeśli właściwe pomiary badawcze są wykonywane na systemach cyfrowych.

Statystyka pomiarów radiometrycznych i minimalna wykrywalna aktywność

Pomiary jądrowe są statystyczne z natury: każdy rozpad jest zdarzeniem losowym. Dla $N$ zliczeń w czasie $t$:

$$\bar{N} = A \cdot \varepsilon \cdot t, \quad \sigma_N = \sqrt{N}$$

Względna precyzja pomiaru wynosi $\sigma_N/N = 1/\sqrt{N}$. Aby osiągnąć precyzję 1%, potrzeba $N \geq 10 000$ zliczeń. Przy aktywności próbki $A = 100$ Bq i efektywności $\varepsilon = 0{,}1$, wskaźnik zliczeń $\dot{N} = 10$ s⁻¹, więc czas dla 10 000 zliczeń: $t = 1000$ s ≈ 17 min. Dla aktywności 1 Bq: $t = 100 000$ s ≈ 28 godzin — co wyjaśnia, dlaczego nisko-tłowe liczniki HPGe na potrzeby środowiskowe mierzą próbki przez wiele godzin lub dni.

Minimalna Wykrywalna Aktywność (MDA, ang. Minimum Detectable Activity) jest standardową miarą czułości stanowiska. Zgodnie z kryterium Currie (1968) dla jednoczesnego uwzględnienia błędu I i II rodzaju:

$$L_D = 2{,}71 + 4{,}65 \cdot \sqrt{B}$$

gdzie $L_D$ jest minimalną liczbą netto zliczeń sygnału wykrywalną z 95% pewnością (przy $\alpha = 0{,}05$ dla błędu I rodzaju i $\beta = 0{,}05$ dla II rodzaju), a $B$ jest mierzoną liczbą zliczeń tła. Przeliczona na aktywność:

$$MDA = \frac{L_D}{t \cdot I_\gamma \cdot \varepsilon(E)}$$

Dla HPGe z tłem $B = 100$ zliczeń w oknie piku po 1000 s pomiaru:

$$L_D = 2{,}71 + 4{,}65 \cdot \sqrt{100} = 2{,}71 + 46{,}5 = 49{,}2 \text{ zliczeń}$$

Przy $I_\gamma = 1$, $\varepsilon = 0{,}1$, $t = 1000$ s:

$$MDA = \frac{49{,}2}{1000 \times 1 \times 0{,}1} = 0{,}49 \text{ Bq}$$

Zmniejszenie tła przez ekranowanie ołowiowe zwiększa czułość proporcjonalnie do $1/\sqrt{B}$. To wyjaśnia, dlaczego spektrometry HPGe do pomiarów środowiskowych działają w masywnych komorach ołowianych grubości 10-15 cm.

Laboratorium NIM jako środowisko nauki eksperymentu

Pracownia NIM jest wyjątkowym środowiskiem dydaktycznym, bo łączy trzy poziomy wiedzy naraz:

Poziom fizyki — student obserwuje zjawiska jądrowe: rozpad beta, gamma, neutronowe, koincydencje anihilacyjne. Każde zdarzenie jest losowe, każde widmo ma charakter probabilistyczny, każdy pik jest rozkładem statystycznym.

Poziom elektroniki — student wchodzi w interakcję z układami analogowymi i cyfrowymi, musi rozumieć impedancję, wzmocnienie, kształtowanie, terminację, szum. Może wpiąć oscyloskop i "zobaczyć" fizyczność sygnałów.

Poziom metrologii — student musi świadomie oceniać niepewności, kalibrować aparaturę, kontrolować warunki pomiaru, porównywać wyniki z oczekiwaniami i dokumentować procedurę.

Połączenie tych trzech poziomów sprawia, że klasyczna pracownia NIM — nawet jeśli nie dostarcza najnowocześniejszego sprzętu — jest wyjątkowo efektywnym środowiskiem rozwijania myślenia eksperymentalnego. Umiejętność "czytania schematu blokowego" jest w istocie umiejętnością rozumowania przyczynowo-skutkowego: co wchodzi do każdego bloku, jak ten blok przetwarza informację, co wychodzi i co może pójść źle.

To myślenie przenosi się bezpośrednio na rozumienie nowoczesnych systemów DPP-FPGA, a dalej na projektowanie i weryfikację skomplikowanych systemów instrumentacyjnych w reaktorach, akceleratorach, satelitach czy systemach dozymetrycznych. Kultura laboratoryjna NIM — precyzja, dokumentacja, reprodukowalność, diagnoza — jest fundamentem pracy eksperymentalnej w fizyce.

Komputerowe oprogramowanie do analizy widm gamma

Oprócz fizycznej obsługi aparatury NIM, student musi opanować oprogramowanie do akwizycji i analizy danych spektrometrycznych. Dostępne są trzy kategorie narzędzi:

Komercyjne oprogramowanie producenta:

Maestro (ORTEC/AMETEK) — klasyczne środowisko MCA dla systemów ORTEC, format danych .Spe. Intuicyjny interfejs, dedykowane biblioteki nukleonowe (GammaVision), automatyczne dopasowanie pików. Standardowe oprogramowanie laboratoryjne w wielu uczelniach.
Genie 2000 (Canberra/Mirion) — zaawansowany pakiet obejmujący akwizycję, kalibrację, analizę pików, raportowanie i zarządzanie bazą danych próbek. Stosowane w profesjonalnych laboratoriach kontroli środowiska i safeguards.
MCDWIN (Amptek) — proste oprogramowanie do modułów MCA Amptek (PX5, DP5). Darmowe, lekkie, przydatne w prostych ćwiczeniach.

Otwarte narzędzia naukowe:

ROOT (CERN) — zaawansowany framework do analizy danych jądrowych, język C++ z interfejsem Python (PyROOT). Standardowe narzędzie w fizyce cząstek i jądrowej, krzywa uczenia się bardzo stroma. Ćwiczenie: histogramowanie, dopasowanie funkcji, analiza statystyczna widm.
Radware (ORNL, David Radford) — pakiet do analizy widm gamma, szczególnie popularny w jądrowej spektroskopii gamma (widma koincydencyjne). Format danych .spk.
Python + biblioteki — skrypty własne z bibliotekami NumPy, SciPy (dopasowanie krzywych, curve_fit), Matplotlib (wizualizacja). Najgibkszy i najdydaktyczniejszy wybór: student widzi cały algorytm analizy.

Przykładowy skrypt Python do analizy piku w widmie:

import numpy as np
from scipy.optimize import curve_fit
import matplotlib.pyplot as plt

# Wczytanie danych .Spe (format ORTEC Maestro)
with open('widmo.Spe') as f:
    lines = f.readlines()
# Znalezienie sekcji z danymi
data_start = next(i for i, l in enumerate(lines) if '$DATA:' in l) + 2
n_channels = int(lines[data_start-1].split()[-1]) + 1
counts = np.array([int(lines[data_start+i]) for i in range(n_channels)])
channels = np.arange(n_channels)

# Dopasowanie Gaussa do piku Cs-137 (661.66 keV)
# Zakładamy wstępnie, że pik jest w okolicach kanału 1324
win = slice(1280, 1370)
x, y = channels[win], counts[win]

def gauss_plus_bg(x, A, x0, sigma, b0, b1):
    return A * np.exp(-(x-x0)**2/(2*sigma**2)) + b0 + b1*x

p0 = [max(y), 1324, 5, min(y), 0]
popt, pcov = curve_fit(gauss_plus_bg, x, y, p0=p0)
A, x0, sigma, b0, b1 = popt
perr = np.sqrt(np.diag(pcov))

fwhm_ch = 2.355 * sigma
n_net = A * sigma * np.sqrt(2*np.pi)

print(f"Centroid: kanał {x0:.2f} ± {perr[1]:.2f}")
print(f"FWHM: {fwhm_ch:.2f} kanałów")
print(f"Zliczenia netto: {n_net:.0f} ± {np.sqrt(n_net):.0f}")

Tego rodzaju skrypt jest przykładem, jak fizyk jądrowy analizuje dane poza środowiskiem graficznym — i dlaczego rozumienie algorytmu za dopasowaniem Gaussa jest ważne dla interpretacji wyników.

Projekt ćwiczenia: kompletne stanowisko spektrometryczne od zera

Doświadczenie pedagogiczne wskazuje, że najbardziej efektywne ćwiczenie z NIM polega na tym, że student sam buduje stanowisko od zera — nie tylko uruchamia gotowe. Etapy budowy stanowiska spektrometrycznego NaI(Tl) z PMT:

Krok 1: Sprawdzenie modułów. Identyfikacja modułów w binie: zasilacz HV, wzmacniacz widmowy, MCA, ew. dyskryminator. Sprawdzenie dokumentacji każdego modułu (katalog producenta). Identyfikacja wejść/wyjść przez oznaczenia frontpanela.

Krok 2: Zasilanie i masa. Włączenie zasilacza binu (główny wyłącznik). Sprawdzenie napięć szyn (+6V, +12V, +24V, -6V, -12V, -24V) multimetrem na wyjściu 44-pin złącza (przy wysuniętym module pomiarowym). Weryfikacja wartości względem specyfikacji standardu NIM.

Krok 3: Podłączenie detektora bez napięcia. Podłączenie kabla BNC od wyjścia detektora NaI+PMT do wejścia wzmacniacza, bez włączania HV. Wpiąć oscyloskop na wyjście wzmacniacza. Włączyć wzmacniacz i obserwować, czy na wyjściu nie widać samodzielnych sygnałów (nie powinno być — tylko szum ~mV).

Krok 4: Włączenie HV. Podłączenie kabla HV od modułu zasilacza HV do złącza HV w obudowie PMT. Ustawienie napięcia HV na 0 V. Włączenie modułu HV. Powolne zwiększanie napięcia do wartości nominalnej (z karty PMT lub poprzednich nastaw) — np. 900 V. Obserwacja na oscyloskopie: pojawienie się przypadkowych impulsów (to sygnały od naturalnego tła promieniowania).

Krok 5: Wzmocnienie i kształtowanie. Regulacja wzmocnienia wzmacniacza tak, by impulsy z tła były widoczne jako sygnały o amplitudzie 1-8 V (zamiast przycięte lub zbyt małe). Regulacja stałej kształtowania (jeśli jest dostępna) na $\tau_p = 1$ μs dla NaI.

Krok 6: Podłączenie MCA. Kabel BNC od wyjścia wzmacniacza do wejścia MCA. Uruchomienie oprogramowania (np. Maestro). Ustawienie czasu pomiaru 60 s. Start akwizycji — obserwacja narastającego histogramu.

Krok 7: Kalibracja. Po zebraniu widma tła, umieszczenie źródła Cs-137 przy detektorze i pomiar 120 s. Identyfikacja piku 661,66 keV. Jeśli widoczny jest też pik anihilacyjny 511 keV (z naturalnej aktywności K-40 lub innych źródeł) — użycie obu do kalibracji energetycznej. Wpisanie kalibracji w Maestro.

Krok 8: Diagnoza. Celowo odłączyć terminację 50 Ω lub zmienić wzmocnienie. Obserwacja jak zmienia się jakość widma. Przywrócenie poprawnych nastaw.

Całe ćwiczenie zajmuje 2-3 godziny laboratoryjne i jest wysoce efektywne: student przez własne ręce przechodzi przez każdy etap, który na gotowym stanowisku byłby niewidoczny. Błędy popełniane przez studenta (np. brak terminacji, złe wzmocnienie, zbyt wysokie HV powodujące wyładowania w PMT) są korygowane na bieżąco przez prowadzącego i stają się trwałą nauką eksperymentalną.

Typowe problemy i ich diagnoza podczas pierwszego kontaktu z NIM

Poniższy katalog objawów i diagnoz jest oparty na wieloletnich doświadczeniach laboratoryjnych:

Objaw	Możliwa przyczyna	Diagnoza i naprawa
Brak sygnałów na wyjściu wzmacniacza	Brak zasilania modułu, brak HV, zły kabel	Sprawdź świecenie LED na module, zmierz HV, zamień kabel
Wszystkie kanały MCA w pierwszym binie	Wzmocnienie za małe lub HV za małe	Zwiększ wzmocnienie lub napięcie HV o 5-10%
Wszystkie kanały MCA w ostatnim binie	Wzmocnienie za duże lub HV za duże	Zmniejsz wzmocnienie lub napięcie HV
Szeroki, niegaussowski pik	Pile-up, złe kształtowanie, duże szumy	Zmniejsz aktywność źródła, zwiększ odległość, sprawdź terminację
Pik się przesuwa w czasie	Dryft wzmocnienia (rozgrzewanie), niestabilność HV	Odczekaj 30 min po włączeniu, wymień moduł HV
Widoczne impulsy przy wyłączonym HV	Oscylacja wzmacniacza lub zakłócenia RF	Sprawdź terminację, zamknij obudowę, odsuń kable zasilania
Linia bazowa faluje z 50 Hz	Pętla masy	Sprawdź uziemienie, użyj jednego punktu masy
Brak koincydencji przy bliskim źródle	Złe ustawienie progu CFD lub okna koincydencji	Sprawdź progi dyskryminatora oscyloskopem, ustaw okno 100-200 ns

Katalog ten jest wartościowym narzędziem samodiagnozy i wprowadza studenta w kulturę systematycznego rozwiązywania problemów technicznych — co jest umiejętnością cenioną zarówno w badaniach, jak i w przemyśle.

Uwagi końcowe dla prowadzącego i studenta

Dokumentacja stanowiska to nie formalność — to dowód powtarzalności. Wynik pomiaru, który nie może być powtórzony przez kogoś innego na podstawie opisu, nie jest wynikiem naukowym. Standard NIM i kultura pomiarów jądrowych od początku kładły nacisk na opis konfiguracji jako integralną część wyniku. Doktorant, który to rozumie już na pierwszym roku, unika dziesiątek godzin straconych na niepowtarzalne eksperymenty. To jest prawdziwy cel ćwiczenia w pracowni NIM.

Typowe błędy studentów przy pierwszym czytaniu schematu NIM

Schematy blokowe stanowisk NIM są z pozoru proste, ale kilka typowych nieporozumień pojawia się regularnie:

Błąd 1: Mylenie złącz BNC i SHV. Kabel BNC i kabel SHV (Safe High Voltage) wyglądają podobnie, ale są różne elektrycznie i mechanicznie — SHV ma cofnięty kontakt chroniący przed porażeniem napięciem HV. Użycie kabla BNC tam, gdzie powinien być SHV (czyli do zasilania HV detektora), grozi porażeniem i uszkodzeniem modułów.

Błąd 2: Podłączanie wyjścia sygnałowego do wejścia HV. W schemacie wyróżnia się wyraźnie wejście HV (zasilacz → detektor, napięcie 500–3000 V) i wejście sygnałowe (detektor → przedwzmacniacz, napięcie mV). Zamiana podłączeń niszczy elektronikę sygnałową.

Błąd 3: Pominięcie terminacji. Wyjście szybkiego dyskryminatora to sygnał NIM fast (50 Ω). Podłączenie go do wejścia o impedancji 1 MΩ (np. wejście oscyloskopu bez terminatora) podwaja widzianą amplitudę i powoduje odbicia. Terminacja 50 Ω na wyjściu lub wejściu jest konieczna.

Błąd 4: Ignorowanie live time. Pomiar przez 300 s real time przy czasie martwym 20% to tylko 240 s live time. Aktywność obliczona bez korekcji jest zaniżona o 20%. Każdy system MCA pokazuje oba czasy — czytaj live time, nie real time.

Błąd 5: Kalibracja "raz na zawsze". Zmiana napięcia HV, wymiana kabla, przeniesienie detektora — każda z tych operacji może przesunąć pozycję piku w widmie. Kalibracja energetyczna jest ważna tylko dla warunków, w których była wykonana. Przed każdą serią pomiarów sprawdź pozycję znanych pików.

Minimalny opis stanowiska do sprawozdania

Według dobrych praktyk laboratoryjnych i wymagań ISO 17025, opis stanowiska pomiarowego powinien zawierać:

typ i model detektora, numer seryjny, datę kalibracji
napięcie HV i polaryzację (np. "HPGe, +3500 V, n-type")
wzmacniacz kształtujący: model, wzmocnienie (gain), czas kształtowania (shaping time)
MCA lub DPP: model, liczba kanałów, czas konwersji ADC
czas pomiaru: real time i live time (czas martwy w %)
geometrię: odległość źródło-detektor, orientacja
temperaturę otoczenia (wpływa na SiPM i PMT)
numer kalibracji energetycznej i efektywności

Brak któregokolwiek z tych elementów uniemożliwia odtworzenie pomiaru przez inny zespół lub w przyszłości.

Ćwiczenia kontrolne

Dla schematu NaI(Tl) -> PMT -> wzmacniacz -> MCA dopisz, gdzie powstaje światło, ładunek, napięcie, amplituda i kanał widma.
Zaznacz na schemacie miejsca, w których błędna nastawa może zmienić liczbę zliczeń.
Przygotuj krótką checklistę do odbioru stanowiska NIM przed ćwiczeniem.

Przejdź do ćwiczenia interaktywnego

Powiązane materiały

Przypisy

U.S. DOE, Standard NIM Instrumentation System, 1990, https://inis.iaea.org/records/5sq70-myv38/files/21080298.pdf?download=1. Zakres wykorzystania: definicje modułu, binu, złączy, sygnałów i zasilania.
↩
WIENER, Introduction to NIM and CAMAC instrumentation, https://www.wiener-d.com/wp-content/uploads/2021/05/WIENER_NIM_CAMAC_introduction.pdf. Zakres wykorzystania: praktyczne przykłady modułów: wzmacniacze, ADC, dyskryminatory i zasilacze.
↩