Streszczenie

Najlepszy sposób zrozumienia aparatury jądrowej to prześledzenie jednego impulsu. W detektorze powstaje krótki sygnał, który wędruje przez kolejne bloki: przedwzmacniacz, wzmacniacz kształtujący, dyskryminator, bramkę, licznik albo analizator.

Tor NIM pokazuje tę drogę w formie fizycznych modułów i kabli. To jego wielka zaleta dydaktyczna: student nie widzi tylko wyniku na ekranie, lecz może wskazać, gdzie powstaje informacja o energii, gdzie zapada decyzja progowa, gdzie liczony jest impuls, a gdzie impuls zostaje odrzucony.1

Rozszerzenie tematu

Punktem startowym jest detektor. W detektorze półprzewodnikowym zdarzenie jądrowe tworzy ładunek proporcjonalny do energii zdeponowanej. W scyntylatorze światło trafia do fotopowielacza i daje impuls prądowy. W liczniku gazowym powstaje lawina jonizacyjna. W każdym przypadku surowy sygnał jest zbyt delikatny albo zbyt krótki, aby od razu traktować go jako wynik.

Przedwzmacniacz jest pierwszym strażnikiem informacji. W spektrometrii półprzewodnikowej często pracuje jako przedwzmacniacz ładunkowy: zamienia ładunek z detektora na napięcie. Jego zadaniem nie jest jeszcze "ładny impuls", lecz zachowanie proporcji między zdarzeniem a sygnałem. Szum, pojemność wejściowa i stabilność zasilania decydują tu o jakości całego widma.

Wzmacniacz kształtujący zmienia impuls w przebieg wygodny do pomiaru amplitudy. Zbyt krótki czas kształtowania może pogorszyć rozdzielczość, zbyt długi zwiększa ryzyko nakładania impulsów. Dlatego dobór kształtowania jest kompromisem między rozdzielczością energetyczną, szybkością zliczeń i odpornością na szum.

Dyskryminator odpowiada na inne pytanie: czy impuls ma zostać uznany za zdarzenie logiczne? Ustawienie progu odcina szum i tło niskoamplitudowe, ale zbyt wysoki próg usuwa także część realnych zdarzeń. Single-channel analyzer idzie krok dalej: przepuszcza tylko impulsy w zadanym oknie amplitudy, czyli tworzy prosty selektor energii.

Analizator wielokanałowy MCA zbiera amplitudy wielu impulsów i buduje widmo. W klasycznym torze analogowym amplituda impulsu po kształtowaniu była przetwarzana na numer kanału. We współczesnym DPP próbkuje się przebieg i filtruje cyfrowo, ale idea pozostaje ta sama: z wielu pojedynczych impulsów powstaje histogram energii.

W dydaktyce NIM pozwala przeprowadzić ćwiczenie bez ukrywania toru wewnątrz oprogramowania. Można odłączyć analizator i obejrzeć impuls na oscyloskopie. Można zmienić próg dyskryminatora i zobaczyć, jak zmienia się liczba zliczeń. Można opóźnić jeden kanał i zrozumieć, czym jest okno koincydencji. To wiedza, której nie daje samo kliknięcie "start acquisition".

Granica bezpieczeństwa jest jasna: ten tor opisujemy jako aparaturę pomiarową. Nie jest to instrukcja budowy układów sterujących detonatorami ani sekwencerów uzbrajania. Te same pojęcia czasu, progów i impulsów mogą występować w wielu dziedzinach, ale w serwisie powinny służyć przede wszystkim metrologii, radiometrii i diagnostyce.

Przedwzmacniacz ładunkowy: fizyka i parametry kluczowe

Przedwzmacniacz ładunkowy (ang. charge-sensitive preamplifier, CSP) jest fundamentalnym elementem toru spektrometrycznego. Jego zadanie jest pozornie proste: zamienić ładunek Q wyemitowany przez detektor na proporcjonalne napięcie wyjściowe V_out = Q/C_f, gdzie C_f to pojemność sprzężenia zwrotnego.

Schemat ideowy: Wzmacniacz transimpedancyjny z rezystorem sprzężenia zwrotnego (opcjonalnie z pojemnością C_f). Gdy detektor emituje ładunek Q, napięcie na wyjściu wzrasta do Q/C_f, po czym powoli opada z stałą czasową τ = R_f · C_f. Typowe wartości: C_f = 1 pF, R_f = 500 MΩ → τ = 500 μs. Powolny zanik (ang. decay tail) jest celowy: pozwala zmierzyć amplitudę przed zanikiem.

Szum wejściowy — czemu jest tak krytyczny: Szum na wyjściu przedwzmacniacza pochodzi z kilku źródeł:

  • Prądowy szum wejściowego tranzystora JFET lub MOSFET (szum białego prądu): $S_i \sim 2qI_G$
  • Napięciowy szum wejściowy (szum termiczny rezystancji i szum 1/f tranzystora): $S_v \sim 4kTR + A_f/f$
  • Szum termiczny rezystorów wejściowych detektora: zależy od temperatury i wartości

Ekvivalentny szum wejściowy wyrażony w elektronach RMS: ENC (Equivalent Noise Charge). Dla HPGe z dobrym przedwzmacniaczem: ENC < 50 e⁻ RMS (co odpowiada FWHM ~1 keV dla Ge). Dla NaI(Tl)+PMT: ENC nie jest limitujące — szum statystyczny scyntylatora jest dominujący (~5-7% FWHM przy 662 keV).

Krioprzewodnictwo (cryogenic operation): Wejściowy tranzystor JFET (np. IF3601, 2SK152) w HPGe często chłodzi się kriogenicznie do 77 K (temperatura ciekłego azotu), co redukuje szum termiczny o czynnik √(77/300) ≈ 0.5. Szum 1/f (flicker noise) jest jeszcze silniej tłumiony przez niską temperaturę. Dlatego kriogeniczne HPGe osiągają FWHM < 2 keV przy 1332 keV (¹³⁷Cs). Cieplejsze (pokojowotempera turowe) HPGe z Peltierem: FWHM 2.5-4 keV. NaI(Tl) w porównaniu: ~7% FWHM = 46 keV przy 662 keV.

Pojemność detektora a szum: ENC rośnie z pojemnością wejściową C_det: $ENC \propto (C_{det} + C_{in}) \cdot \sqrt{S_v}$. Dlatego: detektory małe (niska pojemność) dają lepszą rozdzielczość; kabel między detektorem a przedwzmacniaczem ma pojemność ~100 pF/m → zbyt długi kabel psuje rozdzielczość; przedwzmacniacz montuje się jak najbliżej detektora (często wewnątrz kriostatu HPGe).

Wzmacniacz kształtujący: teoria filtrów i dobór parametrów

Wzmacniacz kształtujący (ang. shaping amplifier, spectroscopy amplifier) ma dwa zadania: wzmocnienie sygnału z przedwzmacniacza (typowo 10-1000×) i kształtowanie go w celu poprawy stosunku sygnał/szum.

Teoria optymalnego filtra: Optymalne filtrowanie szumu białego następuje przez filtr pasmowy (bandpass filter), a optymalny filtr dla 1/f noise to inny kształt. W praktyce używa się filtrów CR-RC^n (Gaussian shaping):

  • CR: różniczkujący, eliminuje długi ogon (decay tail) z przedwzmacniacza
  • RC^n: n-krotne całkowanie (n=4-8), tworzy quasi-gaussowski pik
  • Stała czasowa τ = τ_CR = τ_RC = "shaping time" lub "integration time"

Dobór czasu kształtowania:

  • Krótki τ (0.5-2 μs): lepszy dla dużych rate (mniej pile-up), ale gorszy ENC (szersze pasmo = więcej szumu)
  • Długi τ (4-16 μs): lepszy ENC (węższe pasmo), ale pile-up przy rate >1000 cps
  • Optymalny τ zależy od relacji szumu termicznego do szumu 1/f, i od pojemności detektora
  • Dla HPGe: typowy τ = 4-12 μs. Dla NaI(Tl): τ = 0.5-3 μs (krótszy bo wzmocnienie PMT daje duży sygnał, a szybkość ważniejsza)

Pile-up (nakładanie impulsów): Jeśli dwa zdarzenia w detektorze są bliżej niż ~5τ, ich impulsy nakładają się. Efekt: mierzony impuls ma amplitudę sumy → rejestrowany w złym kanale MCA (suma energii zamiast każdej osobno). Rozwiązanie: pile-up rejector (PUR) — układ logiczny sprawdzający czy impuls jest "czysty" przez pewien czas przed i po.

Wzmacniacz unipolarny vs. bipolarny: Wzmacniacz spectroscopy może generować impuls unipolarny (tylko gaussowskie wzniesienie) lub bipolarny (pik + depresja, sumuje się do zera). Bipolarny: lepsza właściwość baseline (DC level) przy wysokich rate, bo suma ładunku impulsu = 0. Unipolarny: lepsza rozdzielczość energetyczna (cały ładunek w piku, a nie rozłożony po depresji).

Dyskryminator: SCA i jego rola w selekcji energii

Dyskryminator (ang. discriminator) to układ progowy: daje sygnał logiczny gdy amplituda impulsu przekroczy próg. Prosta realizacja: komparator (LM339, AD8611) z regulowanym napięciem odniesienia.

Dyskryminator prędkości przejścia zera (CFD — Constant Fraction Discriminator): Klasyczny dyskryminator progowy cierpi na "walk" — impuls wolny przekracza próg później niż szybki tej samej amplitudy. CFD eliminuje walk przez podział impulsu na dwie kopie: jedną opóźnioną, drugą zanegowaną i osłabioną, a następnie sumuje je. Wynikowy sygnał ma zero-crossing (przejście przez zero) w stałym momencie od szczytu impulsu, niezależnie od amplitudy. Stosowany wszędzie tam, gdzie precyzja czasowa jest ważna: koincydencje, TOF, pomiary czasu życia.

Single Channel Analyzer (SCA): Rozszerzony dyskryminator z dolnym i górnym progiem (Lower Level Discriminator, Upper Level Discriminator — LLD, ULD). Wyjście logiczne tylko gdy: LLD < amplituda < ULD. Pozwala wyselekcjonować impulsy odpowiadające konkretnemu obszarowi energii. Klasyczne zastosowanie: zliczanie fotopieku ¹³⁷Cs (661.7 keV) bez tła Comptona. SCA + scaler = prosty "single-energy counter". Wada SCA: brak pełnego widma — tylko integracja w danym oknie.

Multi Channel Analyzer (MCA): Zbiera całe widmo amplitud. W klasycznym MCA (analog):

  • Analog-to-Digital Converter (ADC): przetwarza amplitudę impulsu na liczbę całkowitą (numer kanału), typowo 1024 lub 4096 kanałów
  • Memory: zlicza impulsy w każdym kanale (histogram)
  • Wyświetlacz lub port komputerowy: transmisja widma do analizy

Czas martwy MCA: czas potrzebny ADC do konwersji (ADC dead time ~ 10-40 μs). Przy rate > 1000 cps istotne straty. Korekcja przez "live time" (rzeczywisty czas zliczania).

Cyfrowy DPP (Digital Pulse Processing): Nowoczesna alternatywa dla całego łańcucha: wzmacniacz → dyskryminator → MCA. Zamiast tego: szybki ADC próbkuje sygnał z przedwzmacniacza (50-250 MS/s, 14-16 bit), a FPGA/DSP wykonuje filtrowanie cyfrowe. Filtr trapezoidalny DPP-PHA: Rise Time (τ_r = czas narastania trapezu), Flat Top (τ_ft = czas plateau dla próbkowania szczytu), Decay Time (τ_d = korekcja stałej zaniku przedwzmacniacza). Zalety DPP: brak pile-up w torze analogowym, lepsza korekcja baseline, konfigurowalny programowo, często lepsza rozdzielczość przy dużych rate.

Typowe ćwiczenie dydaktyczne: charakteryzacja toru NIM krok po kroku

Ćwiczenie w laboratorium fizyki jądrowej z torem NIM jest zazwyczaj podzielone na etapy, z których każdy uczy innego aspektu:

Etap 1: Obserwacja sygnałów na oscyloskopie
Źródło: ²⁴¹Am (60 keV alfa+gamma) lub ¹³⁷Cs (662 keV gamma), detektor NaI(Tl) + PMT. Kolejność obserwacji:

  • Wyjście PMT (przez 50 Ω): krótkie negatywne impulsy, ~100 ns, amplituda 10-100 mV
  • Wyjście przedwzmacniacza (jeśli separowany od PMT): impulsy ~1-10 μs, amplituda 10-100 mV
  • Wyjście wzmacniacza kształtującego: gaussowski pik, ~2-10 μs, amplituda 1-10 V
  • Porównanie: jak zmienia się kształt i amplituda na każdym etapie?

Etap 2: Pomiar plateau dyskryminatora
Ustaw MCA lub skaler. Zmieniaj próg dyskryminatora od 0 do maksymalnego napięcia wyjściowego wzmacniacza. Mierz rate zliczania przy każdym progu. Krzywa "rate vs. próg" pokazuje: obszar szybkiego spadku zliczeń (próg w trakcie ciągłego tła Comptona), ewentualny "half-plateau" przy fotopiku, szybki zanik przy progu wyższym od amplitudy fotopiku.

Etap 3: Widmo energetyczne z MCA
Podłącz wyjście wzmacniacza do wejścia MCA (ADC). Zbierz widmo przez 5-10 minut. Zidentyfikuj: fotopik (ostra linia przy amplitudzie odpowiadającej pełnej energii gamma), край Comptona, Compton continuum, ewentualnie pik backscatter. Oblicz rozdzielczość: FWHM/E_peak × 100% [w %]. Porównaj z danymi producenta.

Etap 4: Wpływ parametrów na widmo
Zmień shaping time (jeśli wzmacniacz ma regulację): obserwuj zmiany FWHM. Zmień HV o ±50 V: obserwuj przesunięcie fotopiku. Dodaj zbyt długi kabel między detektorem a przedwzmacniaczem: obserwuj pogorszenie FWHM (wzrost pojemności kabla). Zmień wzmocnienie wzmacniacza (gain): skaluje całe widmo bez zmiany FWHM w [keV].

Etap 5: Wprowadzenie koincydencji
(dla zaawansowanych) Drugi detektor obok pierwszego. Wyjście dyskryminatora obu detektorów do bramki AND (moduł NIM). Wyjście AND jako gate do MCA. Zbierz widmo w konfiguracji koincydencji. Porównaj z widmem bez koincydencji. Efekt: eliminacja tła kosmicznego i gamm single-channel (komptonowskie przypadkowe).

Kalibracj energetyczna: od kanałów MCA do keV

Po zebraniu widma, MCA daje wynik w kanałach (numerach od 0 do 4095). Przeliczenie na energię wymaga kalibracji przez znane energetycznie źródła:

Procedura kalibracji:

  1. Zmierz widma minimum 3 źródeł o różnych energiach (np. ²⁴¹Am 59.5 keV, ¹³⁷Cs 661.7 keV, ⁶⁰Co 1173+1332 keV)
  2. Zidentyfikuj fotopiki — piki gaussowskie odpowiadające pełnej energii
  3. Wyznacz pozycję każdego fotopiku (dopasowanie gaussa lub centroid): channel_peak
  4. Dopasuj liniową zależność E(keV) = a + b × channel_peak. Współczynniki a i b = parametry kalibracji
  5. Sprawdź liniowość: reszty z dopasowania powinny być < 0.1 kanału (dla ADC 4096 kanałów = < 0.5 keV)

Nieliniowość ADC i MCA:
Nawet dobre ADC mają nieliniowość różniczkową (DNL — Differential Non-Linearity): różne szerokości kanałów. DNL < 0.5% typowo dla ADC wilsona (Wilkinson) używanego w MCA. W widmach gamma, DNL objawia się jako "kanałowy jitter" — pik gaussowski zamiast gładkiego ma podstruktury. Test DNL: widmo pulserowe (generator impulsów monotonnych) powinno dać delta w jednym kanale.

Kalibracja efektywności absolutnej: Nie wystarczy znać energie — potrzebna jest też efektywność absolutna ε(E) do obliczenia aktywności. Procedura: zmierz widmo źródła wzorcowego (certyfikowany ²¹³Bi, ¹⁵²Eu, ⁶⁰Co), znana aktywność A_ref. Oblicz ε(E) = N_detected/(A_ref × I_γ × t_live). Zbuduj krzywą ε(E) przez interpolację (dla detektora NaI(Tl): 3" × 3", geometria 10 cm, typowo ε(662 keV) = 5-15%).

Analizatory wielokanałowe: porównanie rozwiązań

Rynek MCA ewoluował od wielkich stojaków z analogowymi ADC do oprogramowania na PC. Aktualne rozwiązania:

Klasyczny MCA (Standalone): Urządzenie z własnym wyświetlaczem i pamięcią. ORTEC 919 Multichannel Buffer, Canberra S100 MCA. Zalety: niezależny od PC, działa autonomicznie. Wady: ograniczone funkcje analizy, trudna integracja.

USB MCA: Kompaktowe urządzenie podłączane do PC przez USB. Ortec digiDART, Amptek MCA8000A, Amptek dp5. Zalety: tanie, kompaktowe, oprogramowanie PC. Wady: zależny od PC, latencja USB ~ms (nieistotna dla większości spektrometrii).

VME/CAMAC ADC: Moduły do szaf VME lub CAMAC. CAEN V785 (Multi-Event ADC), SILENA 7423. Zalety: integracja z dużymi systemami DAQ, szybkość. Wady: koszt infrastruktury VME.

DPP (Digital Pulse Processor) zintegrowany: CAEN DT5730 (250 MS/s, 14-bit, 8 kanałów, DPP-PHA firmware), XIA Pixie-4, Struck SIS3316. Łączy funkcje wzmacniacza, dyskryminatora i MCA w jednym module cyfrowym. Wejście bezpośrednio z przedwzmacniacza detektora. Całkowity tor sygnałowy: detektor → przedwzmacniacz → DPP → PC. Eliminuje wzmacniacz analogowy, dyskryminator i tradycyjny MCA.

Oprogramowanie do analizy widm:

  • Genie 2000 (Mirion): przemysłowy standard, raportowanie ISO 17025 ready
  • Gamma Vision (ORTEC): podobny do Genie 2000, Windows only
  • InterWinner (ITECH): lekki, dobra analiza multipletów
  • FitzPeaks (freeware): stary, ale bardzo dobry
  • GammaVision Express (online): nowe rozwiązanie SaaS
  • W Polsce: CLOR i niektóre uczelnie używają Genie 2000 lub własnych programów napisanych w Pythonie/LabVIEW

Tory alternatywne: NIM vs. USB-DAQ vs. DPP w 2025 roku

Wybór między klasycznym torem NIM a nowoczesnymi rozwiązaniami cyfrowymi zależy od kontekstu:

Kryterium Klasyczny NIM USB MCA DPP (CAEN)
Koszt Średni (moduły używane) Niski Wysoki
Elastyczność Bardzo wysoka Niska Wysoka
Rozdzielczość Dobra (HPGe: <2 keV) Średnia Najlepsza
Rate capability Ograniczona (pile-up) Ograniczona Bardzo dobra
Wartość dydaktyczna Najwyższa (fizyczny podgląd) Niska Średnia
Serwis/utrzymanie Trudny (stare moduły) Łatwy Łatwy
Integracja systemów Trudna Łatwa (USB) Łatwa (Ethernet)

W dydaktyce physyki jądrowej klasyczny tor NIM jest niezastąpiony właśnie ze względu na najwyższą wartość dydaktyczną — możliwość "zobaczenia" każdego etapu toru sygnałowego. Studenci doktoranccy powinni znać oba światy: rozumieć co dzieje się wewnątrz, gdy klikają w interfejs DPP.

Integracja toru NIM z komputerem: od RS-232 do TCP/IP

Historyczny problem toru NIM: brak standardowego interfejsu komputerowego. Rozwiązania w czasie:

RS-232 i GPIB (lat 80.): Wiele MCA i zasilaczy HV miało interfejs RS-232 lub GPIB (IEEE-488). Komunikacja: komenda tekstowa ("SET HV 1200"), odpowiedź ASCII. Program zbierający dane: QuickBasic, Fortran, C na DOS/Windows. Wadą: wolna prędkość (~1200-9600 baud), brak standaryzacji komend.

SCSI/CAMAC (lat 90.): Wiele systemów używało magistrali CAMAC z kontrolerem SCSI lub ISA. Program: biblioteka CAMAC w Fortran lub C, wywołania bezpośrednio do sprzętu.

USB/PCI (lata 2000.): Większość nowych MCA i zasilaczy to USB lub PCI. Biblioteki SDK: C, C++, Python wrappers. Łatwa integracja z LabVIEW (National Instruments), Matlab, Python (PyUSB, ctypes).

Ethernet (od ok. 2010): Nowoczesne systemy (CAEN SY4527, Ortec DSPEC 502) mają interfejs Ethernet z HTTP REST API lub OPC-UA. Możliwa integracja z EPICS, TANGO, LabVIEW. Zdalne sterowanie przez sieć (ważne przy pracy ze źródłami radioaktywnymi za osłoną).

Przykład skryptu Python (zbieranie widma przez USB):
Komunikacja z ORTEC digiDART przez USB HID:

  1. Zainstaluj pakiet ORTEC SDK (Windows/Linux)
  2. Importuj bibliotekę (import ortec_sdk)
  3. Połącz z urządzeniem: device = ortec_sdk.connect()
  4. Ustaw parametry: device.set_hv(1200), device.set_gain(10.0)
  5. Uruchom pomiar: device.start()
  6. Pobierz spektrum po czasie T: spectrum = device.get_spectrum()
  7. Analizuj: szukaj piku w kanałach, oblicz FWHM, kalibruj energię

Współczesny doktorant powinien umieć napisać taki skrypt i zintegrować go z Jupyter Notebook lub matplotlib. Ale rozumienie tego co się dzieje "za warstwą SDK" — to daje klasyczny tor NIM.

Szum i interferencje w torze NIM: skąd pochodzą artefakty w widmach

Widmo gamma rzadko jest idealnie czyste. Najczęstsze artefakty i ich przyczyny:

Podwyższone tło przy niskich energiach: Impuls szumu z elektroniki (wzmacniacz, przedwzmacniacz, PMT dark current) rejestrowany jako zdarzenie o "małej energii". Rozwiązanie: podnieść próg LLD dyskryminatora powyżej obszaru szumów. Test: zmierzyć widmo bez źródła — tło szumowe widoczne przy kanałach 0-50.

"Ciągłe" tło na całym widmie: Fotony Comptona (poniżej krawędzi Comptona), promieniowanie kosmiczne (miony), fotony tła środowiskowego (⁴⁰K, ²³⁸U, ²³²Th z materiałów budowlanych). Nie jest to artefakt elektroniczny — jest to fizyczne tło promieniowania. Redukcja: osłony ołowiane (5-10 cm Pb), wymiana ołowiu "wiekowego" (stary ołów ma mniejszy Pb-210), ewentualnie Compton suppression (koincydencja ze scyntylatorem BGO wokół HPGe).

Piki "sum" (sum peaks): Przy dużym rate, dwa fotony rejestrowane jako jeden. Suma energii: E_sum = E₁ + E₂ → pik przy E₁+E₂. Typowy przykład: ⁶⁰Co, fotony 1173 i 1332 keV, "pik sumacyjny" przy 2505 keV. Intensywność piku sumacyjnego rośnie z aktywnością i z efektywnością geometryczną (bliższe źródło = wyższy kąt bryłowy = więcej sum). Korekcja: zmniejszyć aktywność lub oddalić źródło.

Pik backscatter (280-300 keV): Fotony Comptona rozpraszają się na materiałach wokół detektora (ołów tarczy, stół laboratoryjny) pod kątem 180° i wchodzą do detektora. Energia fotonu backscatter: $E_{bs} = m_e c^2 / (1 + m_e c^2 / E_0)$. Dla ¹³⁷Cs (662 keV): E_bs ≈ 184 keV. Pik backscatter jest zawsze widoczny przy pomiarach ¹³⁷Cs — nie jest artefaktem elektronicznym.

Pik anihilacyjny (511 keV): Widoczny gdy źródło emituje pozytony (β⁺) lub gdy foton wysokoenergetyczny inicjuje kreację par. Pary elektron-pozytron anihilują → dwa fotony 511 keV. Typowy: ²²Na (β⁺ + 1274 keV gamma) → pik 511 keV.

Pik rentgenowski: Materiały wokół detektora fluoryzują przy napromieniowaniu gamma — emitują charakterystyczne fotony X. Ołów: Pb Kα ≈ 72-87 keV. Widoczny jako wyraźny pik w widmach przy niskich energiach. Identyfikacja: sprawdzić energię rentgenowską danego materiału.

Interferencje EMC: Falownik (inwerter) klimatyzacji w sali laboratoryjnej może generować impulsy EMC rejestrowane przez detektor. Objaw: piki przy stałych amplitudach, niezwiązanych z energiami gamma. Test: wyłączyć klimatyzację i porównać widma. Rozwiązanie: ekranowanie detektora, filtr sieciowy, uziemienie.

Detektory BGO i NaI w Compton Suppression Spectrometer

Compton Suppression Spectrometer (CSS) jest zaawansowanym torem NIM, w którym HPGe jest otoczony scyntylatorem BGO (Bi₄Ge₃O₁₂). System działa przez antykoincydencję:

Zasada działania: Jeśli foton gamma rozproszy się Comptonem w HPGe i częściowo "ucieknie" do BGO, jest to rejestrowane jednocześnie przez oba detektory w oknie koincydencji (~100 ns). Sygnał logiczny z BGO (przez dyskryminator) używany jest jako VETO (anty-koincydencja) do MCA HPGe — zdarzenie jest odrzucane jeśli BGO jest "aktywny" w tym samym czasie.

Efekt w widmie: Compton continuum dramatycznie zmniejsza się (redukcja 5-50×), bo "uciekające" elektrony Comptona są odrzucane. Szczyt fotopiku pozostaje (foton całkowicie pochłonięty przez HPGe nie trafia do BGO). Stosunek fotopik / tło Comptona może wzrosnąć 5-20×.

Tor logiczny CSS:

  1. HPGe signal → preamplifier → shaping amplifier → MCA input
  2. BGO signal → PMT → NIM discriminator (próg powyżej szumu PMT)
  3. HPGe timing signal → CFD → start Gate Generator (okno 100-200 ns)
  4. BGO logic signal → OR gate (16 segmentów BGO) → ANTI-COINCIDENCE wejście MCA lub Gate
  5. MCA: Accept signal = HPGe signal AND NOT (BGO signal in gate)

Moduły NIM potrzebne do CSS: co najmniej 2 przedwzmacniacze (HPGe + BGO×16), 16 dyskryminatorów BGO, 16 CFD lub dyskryminatorów szybkich, moduł OR (fanout+OR dla sygnałów BGO), generator bramki (gate generator), moduł anty-koincydencji (anti-coincidence unit lub logic AND NOT), wzmacniacz kształtujący HPGe, MCA.

CSS jest używany w spektrometrii środowiskowej i w fizyce jądrowej dla zwiększenia czułości wykrywania słabych linii gamma na tle Comptona.

Tor NIM w badaniach wstrzelań: pomiary pulsacyjne

Specyficznym zastosowaniem toru NIM są badania, gdzie detektor musi rejestrować tylko podczas krótkich "impulsów" sygnału, a nie ciągle. Typowy przykład: doświadczenia z pulsacyjnymi źródłami neutronów lub generatorami neutronów.

Neutronowy generator pulsacyjny (D-T, D-D): Generator neutronów (np. Sodern GENIE 16, ADELPHI DD109) emituje neutronowe "piki" o czasie trwania 1-100 μs, z częstością powtarzania 1-100 Hz. Pomiędzy pulsami: brak sygnału lub bardzo niska aktywność tła.

Bramkowanie z synchronizacją do pulsów: Tor NIM: sygnał TTL "trigger" z generatora neutronów → delay (opóźnienie, żeby zaczekać na "prompt" gammy) → gate generator (okno 50-500 μs) → gate wejście MCA. MCA rejestruje tylko w trakcie bramki. Cel: pomiar opóźnionych fotonów gamma (po rozpadzie produktów aktywacji) bez tła "prompt".

Pomiar czasu lotu (TOF — Time of Flight): Przy precyzyjnych pomiarach sekcji przekroju czynnego, czas między pulsem neutronowym a zdarzeniem w detektorze = informacja o energii neutronów. Tor TOF: trigger z generatora → start TDC (Time-to-Digital Converter). Detektor gamma → stop TDC. Widmo TDC = widmo energetyczne neutronów (E = ½mv², v=L/t). Rozdzielczość TDC: typowo 1-2 ns. Przy drodze lotu L=5 m, rozdzielczość energetyczna neutronów 1 MeV: ΔE/E ~ 2×Δt/t = 2×1ns/500ns = 0.4%.

Bezpieczna praca z aparaturą NIM: kwestia kabli i zasilania

Niekompatybilność BNC/SHV: Złącza BNC i SHV wyglądają podobnie i mogą być fizycznie połączone. Jednak BNC jest znamionowane do 750 V max (i tylko do 500 V dla złączy SHV-mate). Podłączenie HV (1000-3000 V) do BNC grozi:

  1. Przebiciem izolacji kabla BNC (izolacja PTFE w RG-58 wytrzymuje tylko ~1.5 kV, za mała margines przy 3 kV)
  2. Porażeniem prądem — kontakt BNC jest dostępny palcem
  3. Uszkodzeniem wejścia modułu (np. wzmacniacza, który oczekuje sygnału -10 V do +10 V)

Zasada: nigdy nie łącz złączy BNC i SHV, nawet jeśli pasują mechanicznie.

Uziemienie i pętle mas: Typowy problem w laboratorium z wieloma instrumentami: "pętla masy" (ground loop). Jeśli dwa urządzenia są uziemione w różnych punktach instalacji elektrycznej (różne kontakty gniazdek), prąd 50 Hz może płynąć przez przewód masy kabla sygnałowego, generując szum 50 Hz. W widmie gamma: pik szumowy przy odpowiedniej amplitudzie lub ogólny wzrost tła w całym widmie. Rozwiązanie: wszystkie urządzenia NIM + komputer + oscyloskop uziemione w jednym punkcie lub przez listwa z wspólną masą.

Kable RG-58 a RG-174: Oba to 50 Ω, ale RG-174 jest cieńszy (zewnętrzna średnica 2.8 mm vs 5 mm). RG-174 ma wyższe tłumienie na metr: przy 100 MHz: ~1.2 dB/m vs 0.3 dB/m dla RG-58. Dla sygnałów NIM (pasmo do ~100 MHz): w krótkich kablach (<1 m) różnica nieistotna. Dla kabli >3 m: używaj RG-58 lub lepszego.

Zabezpieczenie przed wyładowaniami elektrostatycznymi (ESD): Moduły NIM mają chronione wejścia, ale niektóre komponenty (wejście JFET przedwzmacniacza) są wrażliwe. Przy podłączaniu/odłączaniu kabli: zakładaj pasek antystatyczny (ESD strap). Przy przechowywaniu modułów: trzymaj złącza zakryte.

Dokumentacja pomiarów i raportowanie wyników toru NIM

Kompletna dokumentacja pomiaru spektrometrycznego powinna zawierać:

Parametry sprzętu:

  • Typ i model detektora, numer seryjny, data ostatniej kalibracji
  • Napięcie HV i polaryzacja (przykład: HPGe ORTEC GEM-C11015-S, S/N 31-P30127A, HV=3500 V ujemne)
  • Wzmacniacz: model, wzmocnienie, shaping time (przykład: ORTEC 572A, Gain=200, τ=4 μs)
  • MCA: model, liczba kanałów, czas konwersji ADC (przykład: ORTEC 919, 4096 kanałów)
  • Oprogramowanie: wersja Genie 2000 lub GammaVision

Parametry pomiaru:

  • Data i godzina startu i stopu
  • Czas rzeczywisty (real time) i czas aktywny (live time)
  • Czas martwy [%] = (real_time - live_time) / real_time × 100%
  • Temperatura otoczenia (wpływa na PMT i SiPM)
  • Geometria pomiaru: odległość źródło-detektor, orientacja

Wyniki i niepewności:

  • Aktywność próbki A [Bq] z niepewnością poszerzoną U (k=2, poziom ufności 95%)
  • Oznaczalność (Minimum Detectable Activity, MDA) — minimalna aktywność, którą można oznaczalnie wykryć przy danych warunkach pomiarowych
  • Krzywa kalibracji efektywności z niepewnością dla każdego punktu

Według normy ISO 17025 i polskich przepisów dotyczących laboratoriów pomiarów jądrowych (PAA), dokumentacja ta musi być przechowywana przez minimum 5 lat i dostępna dla inspekcji.

Tor NIM do pomiarów alfa: specyfika i trudności

Detektory alfa (półprzewodnikowe detector Si-PIN lub PIPS — Passivated Implanted Planar Silicon) wymagają specyficznego toru NIM ze względu na właściwości cząstek alfa:

Fizyczne właściwości cząstek alfa: Cząstka alfa (jądro helu-4) ma krótki zasięg w materii: w powietrzu ~3-7 cm przy energiach 4-9 MeV, w krzemie ~20-50 μm. Pomiary alfa muszą odbywać się w wysokiej próżni lub w gazie helowym — inaczej zasięg cząstki alfa nie pozwoli dotrzeć do detektora. Pomiar aktywności alfa w laboratorium: komora alfa z pumpą próżniową (ciśnienie < 1 Pa).

Tor sygnałowy dla detektora PIPS:

  • Detektor PIPS (np. Canberra A300-15AM): 15 mm², grubość aktywna 300 μm, pojemność ~50 pF, napięcie +50 V
  • Przedwzmacniacz ładunkowy: niski szum (ENC < 500 e⁻ dla Si-PIN), τ_decay = 50-100 μs
  • Wzmacniacz kształtujący: τ_shaping = 4-10 μs (kompromis między szumem a pile-up)
  • MCA: 512-4096 kanałów; dla widm alfa: rozdzielczość FWHM < 30 keV wymagana dla różnicowania izotopów Am, Pu, U

Kalibracja widm alfa:
Standardowe energie alfa: ²⁴¹Am (5485.56 keV, 84.8% + 5442.80 keV, 13%), ²³⁸Pu (5499.03 keV, 70.9% + 5456.30 keV, 28.7%), ²³²U (5263.5 keV, 68.6%). Rozdzielczość typowa dla PIPS w próżni: FWHM ≈ 20-30 keV (0.3-0.5%). Dla porównania HPGe przy gamma 1332 keV: FWHM ≈ 2 keV (0.15%).

Kołowy system kaskadowych pomiarów alfa + gamma (alfa-gamma coincidence): Przy pomiarach Pu i Am w próbkach środowiskowych: pomiar alfa (identyfikacja izotopów) + pomiar gamma w koincydencji (dla ²⁴¹Am: gamma 59.5 keV). Koincydencja alfa-gamma może być użyta do weryfikacji i redukcji tła.

Pomiary beta: od licznika GM do spektrometrii

Elektrony beta mają złożony spektrum energetyczny (rozkład Fermiego od 0 do E_max) i wymagają innego podejścia:

Zliczanie beta z GM: Najprostsze: liczyć impulsy z licznika GM (np. LND 70x, Hamamatsu R9880U lub GMI-16). Energia nie jest mierzona — tylko aktywność. Korekcje: tło naturalne, czas martwy, efektywność geometryczna (kąt bryłowy). Stosowane w: monitoringu zanieczyszczeń beta w środowisku (np. ⁹⁰Sr/⁹⁰Y), dozymetrii.

Spektrometria beta z detektorem plastikowym:
Scyntylator plastikowy (BC-400, EJ-200) ma niskie Z (głównie C, H) i minimalny rozpraszanie Comptona → prawie liniowy odpowiedź na elektrony. Tor: detektor plastik + PMT → preamplifier → shaping → MCA. Rozdzielczość: ~30% FWHM przy 1 MeV (słaba, bo beta jest ciągłe i rozpraszanie X duże). Stosowane w: pomiarach ⁹⁰Y (E_max = 2.28 MeV), kalibracyjnych pomiarach spektrów beta.

Licznik proporcjonalny 2π i 4π dla beta: Dla absolutnych pomiarów aktywności: geometria 2π (detektor nad płaską próbką) lub 4π (próbka w centrum sfery gazowej). Wysokie pole elektryczne w komorze proporcjonalnej, gaz napełniający P-10 (Ar+CH₄). Efektywność geometryczna: ~50% (2π) lub ~100% (4π) dla emisji izotropowej. Stosowane w laboratoriach kalibracyjnych (CLOR, BIPM).

Rola toru NIM w historii fizyki cząstek

Tor NIM nie był tylko narzędziem pomiarowym — był infrastrukturą, na której zbudowano odkrycia fizyki XX wieku:

CERN 1960-1980: Odkrycie cząstki W i Z (1983, Rubbia i van der Meer, Nobel 1984) wymagało toru NIM do identyfikacji elektronów i mionów z kolajdera SppS. Detektory UA1 i UA2 używały kombinacji komór drucikowych (MWPC), TPC i scyntylatorów, wszystkie ze wzmacniaczami kształtującymi i MCA (kaskady CAMAC z ADC).

SLAC 1969-1979: Doświadczenia DIS (głęboko nieelastyczne rozpraszanie e-p), które odkryły kwarki (Nobel 1990, Friedman, Kendall, Taylor), używały precyzyjnych spektrometrów magnetycznych z wieloma detektorami scyntylacyjnymi. Tor NIM zapewniał precyzyjne timing i selekcję zdarzeń.

Odkrycia w Polsce: W INP PAN Kraków (Instytut Fizyki Jądrowej) tor NIM był używany w badaniach struktury jąder atomowych, przede wszystkim w spektrometrii gamma izotopów neutronowo-bogatych. Lata 1970-1990: eksperymenty z reaktorem Krakus dostarczały neutrony, a tory NIM zbierały widma gamma z detektorów Ge(Li) (poprzedniki HPGe).

Eksperymenty historyczne z CAMAC: LEP (Large Electron-Positron Collider, 1989-2000) używał hybrydy NIM+CAMAC. Sygnały z detektorów wchodziły do modułów NIM (dyskryminatory, logika), a cyfrowizacja i DAQ przez CAMAC. Eksperymenty OPAL, DELPHI, ALEPH, L3 — każdy z tysiącami kanałów NIM.

Przyszłość toru NIM: czy modułowość przeżyje

Modułowość NIM jest wartościową ideą architektoniczną, która wykracza poza lata 1960:

Analogowe NIM (klasyczne): Wciąż produkowane przez CAEN, ORTEC/Ametek, Mesytec, Wiener. Nowe moduły NIM pojawiają się co roku. Powód: wiele laboratoriów ma sprawne NIM biny i stare moduły. Rynek serwisowy i nowych modułów jest wciąż aktywny.

NIM cyfrowy: Niektóre firmy produkują moduły NIM z cyfrowym wnętrzem i analogowym interfejsem NIM. CAEN DT5780 (dual DPP w formacie pulpitu), Struck SIS3302 (flash ADC w formacie NIM). Łączą elastyczność NIM z możliwościami cyfrowego przetwarzania.

μTCA i ATCA jako następcy VME: W fizyce cząstek i systemach LLRF (Low Level RF) akceleratorów, μTCA (micro-Telecommunications Computing Architecture) i ATCA (Advanced Telecommunications Computing Architecture) zastępują VME dla szybkiego DAQ. Nie zastępują NIM — NIM i μTCA koegzystują różnych warstwach systemu.

FPGA jako "moduł universalny": Moduły NIM z FPGA (np. Struck SIS3316 flash ADC 250 MS/s + Xilinx Kintex-7) są de facto programowalnymi modułami — jeden moduł fizyczny, wiele różnych funkcji przez firmware. Eliminuje potrzebę wymiany modułów przy zmianie eksperymentu.

Rynek chińskich zamienników NIM: Od ok. 2015 roku pojawiły się chińskie producenci modułów NIM (Jiuan, Zhuhai Ritec, inne). Tańsze, ale zazwyczaj bez certyfikacji (CE, UL). Używane głównie w edukacji i podstawowych laboratoriach. Jakość bywa zmienna — warto sprawdzić specyfikacje kluczowych parametrów (szum wzmacniacza, liniowość ADC).

Wnioski: dlaczego doktorant fizyki jądrowej powinien znać tor NIM

Znajomość toru NIM jest fundamentem doktorantury w fizyce eksperymentalnej:

  1. Myślenie blokowe: Każdy system pomiarowy (nawet zintegrowany DPP) ma bloki funkcjonalne. Zrozumienie toru NIM pozwala myśleć o systemie strukturalnie — gdzie jest wzmocnienie, gdzie dyskryminacja, gdzie zegar.

  2. Diagnoza problemów: Kiedy widmo wychodzi złe, trzeba wiedzieć, który blok jest winny. Doktorant znający tor NIM może systematycznie testować: "odłącz wzmacniacz, podłącz oscyloskop, sprawdź impuls na wyjściu".

  3. Rozumienie ograniczeń: Każdy system ma ograniczenia — czas martwy, pile-up, szum, nieliniowość. Tor NIM "odsłania" te ograniczenia — nie można ich ukryć za warstwą oprogramowania jak w zintegrowanym systemie DPP.

  4. Komunikacja z inżynierami i technikami: W każdym laboratorium są technicy, którzy przez dekady pracowali z NIM. Znajomość NIM umożliwia prawdziwą komunikację — "potrzebuję dyskryminatora z progiem 20 mV i wyjściem NIM fast" to pełna specyfikacja, zrozumiała przez każdego.

  5. Historia nauki jako kontekst: Tor NIM jest częścią historii fizyki XX wieku. Eksperymenty, które odkryły strukturę jądra, W i Z bozony, kwarki, koincydencje beta-gamma — wszystkie były zbudowane na NIM. Zrozumienie narzędzia to zrozumienie kontekstu, w którym były wykonane.

  6. Polska tradycja eksperymentalna: W polskich ośrodkach (IFJ PAN, NCBJ, wydziały fizyki UJ, UW, PWr) tor NIM był i jest obecny. Wiele wyników naukowych w polskich publikacjach jądrowych z lat 70.-2000 było generowane przez tory NIM. Doktorant znający NIM ma dostęp do rozumienia tej literatury.

Trzy klasyczne zadania toru NIM: zliczanie, spektrometria i timing

Dla porządkowego usystematyzowania — tor NIM może realizować trzy fundamentalnie różne zadania pomiarowe, i każde wymaga innego zestawu modułów:

Zadanie 1: Zliczanie (rate pomiarowy)
Cel: zmierzenie aktywności próbki. Informacja: liczba zdarzeń na sekundę (cps, cpm, Bq).
Tor minimalny: Detektor → Preamplifier → Dyskryminator → Skaler. Ewentualnie wzmacniacz kształtujący między preamplifier a dyskryminatorem (jeśli sygnał za słaby lub zaszumiony). MCA nie jest potrzebne — tylko liczba zliczeń w wybranym czasie. Typowe zastosowania: kontrola zanieczyszczenia radioaktywnego, monitoring promieniowania, weryfikacja aktywności próbki przed spektrometrią.

Kluczowy parametr: próg dyskryminatora. Za niski → szum zliczany. Za wysoki → realne zdarzenia odrzucane. Metoda wyznaczenia optymalnego progu: plateau curve (patrz wyżej).

Zadanie 2: Spektrometria energetyczna
Cel: zmierzenie widma energetycznego emitowanego promieniowania. Informacja: histogramy amplitud (energia depozytowana w detektorze).
Tor: Detektor → Preamplifier → Shaping Amplifier → MCA (ADC + histogramer). Opcjonalne: PUR (pile-up rejector) między wzmacniaczem a MCA, Dyskryminator SCA jako bramka do MCA (pomiar w oknie energetycznym).

Kluczowe parametry: shaping time τ, gain wzmacniacza, liczba kanałów MCA, czas pomiaru. FWHM w keV jest miarą jakości toru. Zastosowania: identyfikacja izotopów promieniotwórczych, analiza aktywacyjna, pomiary rozpadów, spektrometria środowiskowa.

Zadanie 3: Timing i koincydencje
Cel: zmierzenie relacji czasowych między zdarzeniami lub selekcja zdarzeń przez koincydencję. Informacja: czas między zdarzeniami lub liczba zdarzeń spełniających warunek koincydencji.
Tor: 2 lub więcej detektorów → 2 lub więcej ścieżek (każda: Detektor → Preamplifier → CFD) → Bramka koincydencji (AND) → Skaler lub MCA (w gałęzi energetycznej jednego detektora).

Kluczowe parametry: rozdzielczość czasowa (CRT), okno koincydencji Δt, źródła jitteru (CFD, scyntylator, PMT TTS). Zastosowania: PET, koincydencje beta-gamma, pomiary czasu życia pozytonów (PALS), Rossi-alpha dla reaktorów.

Specyfika pracy z detektorem HPGe: kriogenia i procedury

Detektor HPGe jest najdroższym i najdelikatniejszym elementem toru NIM w laboratorium spektrometrycznym. Procedury operacyjne:

Chłodzenie kryostatem azotowym:

  • Pojemność kriostatu: 1.5-7 litrów ciekłego azotu (LN₂)
  • Czas chłodzenia od RT do 77 K: 24-48 godzin przy pierwszym uruchomieniu
  • Uzupełnianie LN₂: co 5-10 dni (czas autonomii zależy od izolacji)
  • Alarm poziomu LN₂: elektroniczny sygnał gdy poziom spada do 30% → alarm w systemie EPICS/TANGO
  • NIGDY nie podawaj HV na ciepły detektor HPGe (uszkodzenie kryształu)

Procedura startu HPGe po wygrzewaniu kriostatu:

  1. Napełnij LN₂ (przez zawór wlewowy na górze kriostatu)
  2. Poczekaj 24h na stabilizację temperatury (temperatura 77 K ± 1 K)
  3. Włącz przedwzmacniacz (napięcie ±12V z NIM bin)
  4. Sprawdź prąd upływu detektora: powinien być <10 pA (przy zimnym HPGe)
  5. Rampa HV do nominalnej wartości (100 V/min)
  6. Sprawdź widmo testowe (źródło ¹³⁷Cs): pozycja piku, FWHM, tło
  7. Jeśli FWHM > 2× nominalna: prawdopodobnie kryształ nie jest w pełni schłodzony lub jest uszkodzony

Elektrody HPGe — typy:

  • n-type HPGe (BEGe, GLOLEX): elektroda wewnętrzna (+), elektroda zewnętrzna (-). Napięcie 3-5 kV
  • p-type HPGe (GEM, GMX): elektroda wewnętrzna (-), elektroda zewnętrzna (+). Napięcie 1-3 kV
  • HPGe "Low Energy" (LEGE): cienka wejściowa elektroda (0.3 μm Lithium surface), dla energii <100 keV
    Błąd polaryzacji: brak lub słaby sygnał (zaporowa praca złącza), przy dużym prądzie upływu możliwe uszkodzenie.

Powiązane kalkulatory i narzędzia

Powiązane artykuły

Ćwiczenia kontrolne

  1. Dla każdego bloku toru NIM wpisz, jaka informacja jest zachowywana: energia, czas, liczba zdarzeń albo decyzja logiczna.
  2. Wyjaśnij, dlaczego zmiana progu dyskryminatora może zmienić wynik aktywności próbki.
  3. Narysuj dwa warianty toru: do prostego zliczania oraz do spektrometrii gamma.

Przejdź do ćwiczenia interaktywnego