Spektrometria neutronowa: kule Bonnera, komory rekombinacyjne i widma neutronów

Streszczenie

Pomiar neutronów przypomina próbę mierzenia deszczu za pomocą powierzchni, która reaguje inaczej na mżawkę, inaczej na grad, a jeszcze inaczej na krople spadające pod kątem. Neutron nie jonizuje bezpośrednio gazu, kryształu ani półprzewodnika; najpierw musi zderzyć się z jądrem albo wywołać reakcję jądrową, a dopiero wtórna cząstka naładowana daje sygnał elektryczny.^1,2

Dlatego spektrometria neutronowa nie jest prostą wersją spektrometrii gamma. Kule Bonnera, wielowarstwowe moderatory i komory rekombinacyjne mierzą nie samo widmo, lecz zestaw odpowiedzi detektorów na widmo; dopiero później, z użyciem funkcji odpowiedzi i metod rekonstrukcji, próbuje się odtworzyć rozkład energii neutronów.²

Rozszerzenie tematu

Dlaczego neutron jest trudny pomiarowo

W liczniku gamma można zacząć od intuicji: foton oddziałuje w detektorze, oddaje energię elektronom, a elektronika mierzy impuls. W przypadku neutronu ta intuicja szybko zawodzi. Neutron jest elektrycznie obojętny, więc nie wytwarza bezpośrednio par jonów ani par elektron-dziura. Sygnał pojawia się dopiero po oddziaływaniu jądrowym: rozproszeniu na jądrze, wychwycie, reakcji (n, alpha), (n, p) albo rozszczepieniu. To wtórne protony, cząstki alfa, jądra litu albo fragmenty rozszczepienia jonizują ośrodek.^1,2

Drugi kłopot polega na tym, że prawdopodobieństwo oddziaływania neutronu bardzo silnie zależy od energii. Dziunikowski i Kalita podkreślają, że wiele przekrojów czynnych dla neutronów termicznych zachowuje się w przybliżeniu zgodnie z regułą 1/v, czyli rośnie, gdy neutron zwalnia.¹ Ten sam neutron po spowolnieniu w moderatorze może być więc znacznie łatwiejszy do zarejestrowania niż neutron szybki. Detektor neutronów jest przez to zwykle układem złożonym: konwerter jądrowy, moderator, geometria, elektronika i algorytm interpretacji tworzą jeden przyrząd pomiarowy.

Trzeci kłopot jest dozymetryczny. W ochronie radiologicznej nie wystarczy powiedzieć, ile neutronów przeleciało przez powierzchnię. Znaczenie ma energia, bo od niej zależą wtórne cząstki, liniowy przekaz energii, współczynniki jakości i przeliczenie fluencji na równoważnik dawki. Maciak zwraca uwagę, że dla pól neutronowych podstawową wielkością metrologiczną jest energetyczny i kierunkowy rozkład fluencji neutronów, a równoważnik dawki otrzymuje się przez ważenie tego rozkładu odpowiednimi współczynnikami fluencja-dawka.²

Detektor termiczny jako punkt wyjścia

Najprostszy dydaktyczny punkt startu to detekcja neutronów termicznych. W licznikach BF3 wykorzystuje się izotop B-10, w licznikach helowych He-3, a w detektorach litowych Li-6. W każdym z tych przypadków neutron termiczny wywołuje reakcję jądrową, której produkty są naładowane i mogą jonizować gaz albo scyntylator.¹

To tłumaczy, dlaczego w praktyce tak często łączy się detektor neutronów z moderatorem. Jeżeli mamy detektor czuły na neutrony termiczne, a pole zawiera neutrony szybkie, można je częściowo spowolnić w materiale bogatym w wodór, na przykład w polietylenie albo parafinie. Neutrony tracą energię w zderzeniach sprężystych z jądrami wodoru, po czym są rejestrowane przez detektor termiczny. Taki układ nie mówi jednak jeszcze, jaka była energia neutronu przed moderacją. Mówi tylko, że po całej drodze przez moderator i detektor powstał rejestrowalny sygnał.

To rozróżnienie jest kluczowe. Licznik z moderatorem może być dobrym monitorem obecności neutronów albo miernikiem dawki po odpowiednim wzorcowaniu, ale nie staje się automatycznie spektrometrem. Spektrometria zaczyna się dopiero wtedy, gdy znamy zależność odpowiedzi układu od energii neutronu i mamy wiele niezależnych odpowiedzi, które da się matematycznie połączyć.

Funkcja odpowiedzi: najważniejsze pojęcie

W spektrometrii gamma często myślimy o piku: foton ma energię, detektor ma kanał, a kalibracja energii przypisuje kanałowi liczbę w keV. W spektrometrii neutronowej trzeba myśleć inaczej. Detektor ma funkcję odpowiedzi R(E), czyli informację, jak silnie reaguje na neutrony o energii E.

Jeżeli w polu neutronowym występuje widmo fluencji Phi_E(E), mierzone zliczenia albo ładunek M_d w danej konfiguracji detektora można zapisać schematycznie jako:

M_d = integral R_d(E) * Phi_E(E) dE.

Po podziale energii na przedziały dostajemy postać macierzową:

M_d = suma_i R_d(E_i) * Phi_E(E_i) * Delta E_i.

Maciak podaje właśnie taki opis dla układów typu Bonner Sphere Spectrometry: odpowiedź pojedynczej konfiguracji detektor-moderator jest splotem funkcji odpowiedzi z widmem neutronów.² To zdanie jest sercem całego tematu. Pojedynczy pomiar nie „pokazuje widma”. Pojedynczy pomiar daje jedną liczbę, w której wymieszane są wszystkie energie zgodnie z czułością układu.

Kule Bonnera

Klasyczny układ kul Bonnera polega na umieszczeniu detektora neutronów termicznych w kulistych moderatorach, zwykle polietylenowych, o różnych średnicach. Oryginalna koncepcja Brambletta, Ewinga i Bonnera wykorzystywała scyntylator LiI(Eu) wzbogacony w Li-6 oraz zestaw kul o średnicach od 2 do 12 cali, dla których wyznaczono funkcje odpowiedzi.²

Sens konstrukcji jest prosty. Dla neutronów termicznych najlepsza może być konfiguracja bez moderatora, bo dodatkowy polietylen tylko osłabia strumień lub zmienia geometrię. Dla neutronów rzędu MeV moderator zaczyna być potrzebny, bo spowalnia neutrony do energii, przy których detektor termiczny reaguje mocniej. Różne grubości moderatora dają różne maksima odpowiedzi w funkcji energii. Zestaw pomiarów wykonanych różnymi kulami tworzy więc kilka rzutów tego samego widma.

W praktyce otrzymuje się wektor zliczeń:

M = [M_0, M_2, M_4, M_6, M_8, M_10, M_12],

gdzie indeksy mogą odpowiadać konfiguracji bez kuli i kulom o różnych średnicach. Do tego potrzebna jest macierz odpowiedzi:

R = detektory x grupy energii.

Zadanie brzmi: znaleźć takie widmo Phi, aby R * Phi możliwie dobrze odtwarzało zmierzone M, a jednocześnie widmo było fizyczne, czyli bez ujemnych fluencji i bez sztucznych oscylacji.

Unfolding, czyli odtwarzanie widma

Problem odwrotny w sferach Bonnera jest źle uwarunkowany. Liczba konfiguracji pomiarowych wynosi zwykle kilka lub kilkanaście, a liczba przedziałów energii w docelowym widmie może być znacznie większa. Maciak opisuje to jako równanie całkowe Fredholma pierwszego rodzaju i wskazuje, że proste odwracanie macierzy może dawać widma niefizyczne: oscylujące albo z ujemnymi wartościami fluencji.²

Dlatego w praktyce stosuje się metody rekonstrukcji: regularyzację, najmniejsze kwadraty z ograniczeniami, metody iteracyjne, maksimum entropii, podejścia Bayesowskie albo algorytmy stochastyczne. Nazwy programów, takie jak GRAVEL, MAXED czy FRUIT, są mniej ważne niż sama zasada: rozwiązanie nie jest prostym wynikiem aparatu. Jest wynikiem aparatu, funkcji odpowiedzi, założeń początkowych i przyjętej metody numerycznej.²

To ma duże konsekwencje dydaktyczne. Dwa różne algorytmy mogą dać podobny równoważnik dawki, ale różniące się szczegóły widma. Ten sam wektor zliczeń może być zgodny z wieloma widmami, jeżeli funkcje odpowiedzi są szerokie i nakładają się na siebie. Dlatego dobry raport ze spektrometrii neutronowej powinien mówić nie tylko „zmierzono widmo”, ale także: jaką macierz odpowiedzi zastosowano, jaki był zakres energii, jakie były ograniczenia, jaki był wektor startowy, jak oszacowano niepewność i jak sprawdzono stabilność rozwiązania.

Co kule Bonnera robią dobrze, a czego nie robią

Kule Bonnera są popularne, bo mogą objąć bardzo szeroki zakres energii neutronów. Dobrze nadają się do pól rozproszonych, mieszanego promieniowania neutronowo-gamma i sytuacji, w których nie znamy z góry kształtu widma. Rozszerzone zestawy ERBSS dodają specjalne wkładki, na przykład metalowe, aby poprawić odpowiedź dla neutronów wysokoenergetycznych powyżej typowego zakresu klasycznych kul.²

Ich ograniczeniem jest rozdzielczość. Moderator nie wybiera jednej energii, tylko tworzy szeroką odpowiedź. Zestaw kul nie działa jak spektrometr germanowy dla gamma, gdzie linia energetyczna może stać się ostrym pikiem. Kule Bonnera dają raczej informację integralną: ile sygnału uzyskano po różnych drogach moderacji. To wystarcza do wielu zastosowań dozymetrycznych, ale wymaga pokory przy interpretacji detali widma.

Drugie ograniczenie jest organizacyjne. Klasyczny pomiar wieloma kulami wymaga powtarzania ustawienia geometrii dla kolejnych moderatorów. Jeżeli pole zmienia się w czasie, trzeba monitorować jego stabilność i wprowadzać korekcje. Maciak opisuje właśnie ten problem jako motywację dla układów z pojedynczą kulą i wieloma detektorami wewnątrz moderatora, czyli koncepcji SSBSS oraz polskiego spektrometru wielowarstwowego SWP-1.²

Spektrometr wielowarstwowy i pomiar pasywny

W układach pasywnych zamiast aktywnego licznika można użyć detektorów aktywacyjnych albo termoluminescencyjnych. Maciak omawia SWP-1, czyli pasywny spektrometr wielowarstwowy oparty na polietylenowej kuli i detektorach MTS-6/MTS-7 rozmieszczonych na różnych głębokościach. Różne głębokości wewnątrz jednej kuli odpowiadają różnym warunkom moderacji, a więc różnym funkcjom odpowiedzi.²

Taka koncepcja rozwiązuje część problemów klasycznego zestawu kul. Geometria zewnętrzna jest jedna, więc nie trzeba przenosić kolejnych moderatorów w ten sam punkt. Detektory pasywne nie mają czasu martwego toru impulsowego ani problemu spiętrzania impulsów. Z drugiej strony wynik nie pojawia się natychmiast: detektory pasywne trzeba odczytać i przeliczyć, a analiza nadal wymaga funkcji odpowiedzi oraz rekonstrukcji widma.

Dla dydaktyki jest to bardzo dobry przykład różnicy między pomiarem aktywnym i pasywnym. Pomiar aktywny daje sygnał w czasie rzeczywistym, ale może mieć problemy z przeciążeniem elektroniki. Pomiar pasywny integruje ekspozycję, jest odporny na część problemów impulsowych, lecz traci informację czasową i przenosi ciężar pracy na odczyt oraz model.

Komora wielosygnałowa KW-1

Komora KW-1 pokazuje inną drogę: zamiast zmieniać kulę wokół detektora, buduje się detektor z wieloma aktywnymi objętościami, które widzą neutrony po przejściu przez różne grubości moderatora. Maciak opisuje ją jako równoważną tkance komorę jonizacyjną z polipropylenowymi elektrodami pełniącymi jednocześnie rolę moderatora. Osiem czynnych objętości komory daje osiem sygnałów, a warstwa zawierająca B-10 zwiększa czułość na neutrony po ich spowolnieniu.²

Z punktu widzenia matematyki to nadal jest spektrometria odpowiedzi. Każdy kanał komory ma własną funkcję odpowiedzi. Neutrony o różnych energiach inaczej przenikają przez kolejne warstwy moderatora, inaczej termalizują się w objętościach i inaczej przyczyniają się do jonizacji gazu. Zestaw sygnałów z kanałów tworzy wektor danych podobny do wektora z kul Bonnera, tylko uzyskany równocześnie w jednym detektorze.

Ten przykład jest ważny także dla kategorii elektroniki i aparatury. Spektrometria neutronowa nie kończy się na materiale detektora. Wymaga stabilnego elektrometru, wielu kanałów akwizycji, korekcji składowej gamma, kontroli ciśnienia gazu, modelowania Monte Carlo oraz kalibracji funkcji odpowiedzi. Bez tej warstwy aparaturowej fizycznie elegancki detektor nie daje wiarygodnego widma.

Komory rekombinacyjne

Rekombinacyjna komora jonizacyjna nie jest zwykłym licznikiem impulsów. Jest detektorem dozymetrycznym przeznaczonym do pól mieszanych, na przykład neutrony plus gamma. Maciak opisuje ją jako komorę równoważną tkance, w której materiał elektrod i gaz mają naśladować właściwości dozymetryczne tkanki, a metoda pomiarowa wykorzystuje zjawisko lokalnej rekombinacji jonów w gazie.²

Idea jest następująca: promieniowanie o różnym liniowym przekazie energii tworzy w gazie ślady jonizacji o różnej gęstości. Przy odpowiednich warunkach część jonów rekombinuje zanim zostanie zebrana przez elektrody. Charakterystyka nasycenia komory, czyli zależność zebranego ładunku od napięcia polaryzującego, zawiera informację o jakości promieniowania. W polu mieszanym pozwala to szacować wielkości dozymetryczne, w tym związane z neutronową składową pola.²

To nie jest jednak metoda „rysowania ostrego widma energii neutronów”. Komora rekombinacyjna jest szczególnie użyteczna w dozymetrii pól mieszanych, gdy istotne jest oszacowanie równoważnika dawki, jakości promieniowania albo rozkładu dawki względem liniowego przekazu energii. Może być elementem spektrometrii w szerszym sensie, ale jej naturalnym językiem jest dozymetria i mikrodozymetria, a nie identyfikacja linii energetycznych.

Fluencja, dawka i współczynniki przeliczeniowe

W wielu artykułach o promieniowaniu można rozdzielić pomiar fizyczny od oceny skutku. W neutronach to rozdzielenie jest szczególnie ważne. Fluencja Phi_E(E) mówi, ile neutronów o danej energii przecina powierzchnię. Równoważnik dawki wymaga przeliczenia tej fluencji przez współczynnik zależny od energii:

H = integral h_Phi(E) * Phi_E(E) dE.

Jeżeli widmo jest źle odtworzone, równoważnik dawki może być źle oceniony nawet wtedy, gdy całkowita liczba zliczeń wygląda stabilnie. Energia neutronów decyduje o tym, jakie cząstki wtórne powstaną i jak silny będzie ich efekt biologiczny. Maciak podkreśla, że dla neutronów współczynnik wagowy promieniowania zmienia się silnie z energią, co wymusza znajomość widma w zastosowaniach ochrony radiologicznej.²

W praktyce studenckiej warto więc rozróżniać trzy poziomy:

• zliczenia albo ładunek z detektora,
• odtworzone widmo fluencji,
• wielkości pochodne, takie jak równoważnik dawki, średnia energia albo udział neutronów termicznych, pośrednich i prędkich.

Pominięcie środkowego poziomu jest czasem możliwe w prostych monitorach dawki, ale wtedy wynik zależy od założonego widma wzorcowania. Gdy pole odbiega od tego wzorca, interpretacja może się pogorszyć.

Najważniejsze źródła niepewności

W spektrometrii neutronowej niepewność nie kończy się na sqrt(N). Statystyka zliczeń jest tylko jednym składnikiem. Ważne są także:

• niepewność funkcji odpowiedzi,
• różnice między geometrią rzeczywistą i modelem Monte Carlo,
• stabilność pola w czasie pomiaru,
• rozproszenie od ścian, podłogi, statywów i osłon,
• czułość na promieniowanie gamma,
• czas martwy i spiętrzanie impulsów w układach aktywnych,
• odczyt i powtarzalność detektorów pasywnych,
• założenia algorytmu unfolding,
• wybór zakresu energii i przedziałów grupowych.

Właśnie dlatego artykuły o statystyce zliczeń, czasie martwym oraz kalibracji i niepewności są potrzebnym zapleczem dla spektrometrii neutronowej. Bez nich student widzi tylko nazwę detektora, a nie wie, co naprawdę ogranicza wynik.

Co powinien umieć student po tym artykule

Po pierwsze, powinien umieć wyjaśnić, dlaczego detekcja neutronów jest pośrednia. Po drugie, powinien rozumieć, że moderator nie jest tylko osłoną, ale częścią funkcji odpowiedzi. Po trzecie, powinien umieć zapisać pomiar jako M = R * Phi i powiedzieć, dlaczego odwrócenie tego równania jest trudne.

Po czwarte, powinien odróżniać trzy rodzaje narzędzi: klasyczny zestaw kul Bonnera, pasywny układ wielowarstwowy oraz komorę rekombinacyjną. Każde z nich odpowiada na trochę inne pytanie. Kule Bonnera są uniwersalnym narzędziem widmowym o szerokich funkcjach odpowiedzi. Układy pasywne dobrze integrują ekspozycję i mogą ograniczać problemy elektroniki impulsowej. Komory rekombinacyjne są szczególnie ważne w dozymetrii pól mieszanych, gdzie interesuje nas nie tylko liczba neutronów, lecz także jakość promieniowania.

Najkrótsze podsumowanie brzmi: spektrometr neutronowy nie „widzi energii” wprost. Zbiera kilka odpowiedzi na to samo pole, a widmo jest wynikiem modelu, kalibracji i ostrożnego rozwiązania problemu odwrotnego.

Rozszerzone zestawy kul Bonnera (ERBSS) i detektory długozasięgowe

Standardowe kule Bonnera (do 12") tracą czułość na neutronach powyżej ~20 MeV, co stanowi problem przy akceleratorach protonowych i spallacyjnych źródłach neutronów (SNS, ESS). Rozszerzone zestawy (Extended Range BSS, ERBSS) stosują modyfikacje:

Kule z ołowiem i polietylenem: Na zewnątrz moderatora polietylenowego umieszcza się warstwę Pb lub Cu. Neutrony prędkie wywołują w nich reakcje (n,xn), generując kilka neutronów wtórnych o niższej energii — które moderator efektywnie zwalnia. Zwiększa to czułość powyżej 20 MeV i umożliwia pomiar spektrów kosmicznych (E > 100 MeV).

Sb-Be i D₂O: Dodatkowe moderatory (deuterowa woda, beryl z antymonem) zmieniają kształt funkcji odpowiedzi, wypełniając luki w widmie.

Detector BSS na kosmiczne neutronowe: Na dużych wysokościach (loty lotnicze, stacje orbitalne) widmo neutronowe sięga GeV. Zestawy ERBSS z detektorem H*(10) stosuje się do oceny dawki personelu lotniczego — coroczna dawka pilota może wynosić 2–5 mSv, głównie od promieniowania kosmicznego.

Pasywne detektory aktywacyjne — metoda foliowa

Gdy pomiar aktywny (elektronika wysokiego napięcia) nie jest możliwy — np. w rdzeniu reaktora lub komorze próżniowej akceleratora — stosuje się detektory aktywacyjne: metaliczne folie lub druty eksponowane na neutronowe pole, a następnie mierzone na spektrometrze gamma po ekspozycji.

Zasada: Neutron wywołuje reakcję (n,γ), (n,p) lub (n,α) w folii. Aktywność produktu A(t) = Φ × σ × N × (1 − e^(−λt_exp)), gdzie t_exp — czas ekspozycji. Po napromieniowaniu mierzy się A przez spektrometrię gamma, skąd wyznacza Φ.

Nuklidy stosowane w foliach aktywacyjnych:

• Au-197(n,γ)Au-198: T₁/₂ = 2,69 d, linia gamma 411 keV; czuły na neutronowe termiczne i rezonansowe (rezonans 4,9 eV). Najszerzej stosowany — złoto nie koroduje, jest jednoizotopowe.
• In-115(n,γ)In-116m: T₁/₂ = 54 min; silny rezonans przy 1,46 eV. Para Au+In pozwala rozróżnić składową termiczną od epitermicznej metodą kadmową.
• Co-59(n,γ)Co-60: T₁/₂ = 5,27 lat; akumulacja Co-60 daje informację o fluencji całkowitej w długich ekspozycjach. Stosowany w dozymetrii reaktorowej (CRP — capsule radiation pressure).
• Dy-164(n,γ)Dy-165: T₁/₂ = 2,33 h; duży σ dla termicznych (2700 barn). Stosowany gdy potrzeba szybkiej odpowiedzi (krótki T₁/₂ pozwala na wielokrotne ekspozycje).
• Ni-58(n,p)Co-58: Reakcja progowa E_th ≈ 0,9 MeV. T₁/₂ Co-58 = 70,8 d. Pozwala selektywnie mierzyć składową prędkich neutronów.
• Al-27(n,α)Na-24: Progowa E_th ≈ 3,3 MeV. T₁/₂ Na-24 = 14,96 h. Stosowany do pomiaru strumienia szybkich neutronów w elektronice lotniczej i kosmicznej (test SEU).

Metoda kadmowa: Otoczenie folii kadmem (Cd pochłania neutronowe < 0,5 eV) pozwala mierzyć osobno składową termiczną i kadmową (epitermiczna+prędka). Różnica aktywności folie bez Cd − folie z Cd = składowa termiczna.

Ograniczenia: Brak informacji spektralnej w czasie pomiaru (rekonstrukcja widma wymaga zestawu wielu folie); nie nadają się do szybkiej kontroli dawki (konieczność aktywacji i pomiaru). Stosowane głównie do wzorcowania i dozymetrii reaktorowej.

Rem-metry i przenośne monitory pola neutronowego

Kule Bonnera wymagają kilkunastu detektorów i zewnętrznego systemu unfolding — zbyt skomplikowane do codziennego monitoringu bezpieczeństwa. W praktyce ochrony radiologicznej stosuje się rem-metry: urządzenia przenośne odpowiadające w sposób zadowalający na pole neutronowe proporcjonalnie do H*(10).

WENDI-II (Wide Energy Neutron Detection Instrument): Moderator polietylenowy w kształcie cylindra z wewnętrzną warstwą wolframu i He-3 licznikiem proporcjonalnym. Funkcja odpowiedzi zbliżona do h(E) z ICRP 74 w zakresie 0,025 eV – 5 GeV. Czułość: 0,5 zliczeń/nSv. Ważący moderator zapewnia „naturalną" odpowiedź bez potrzeby unfolding. Stosowany w ochronie przy cyklotronach, akceleratorach i reaktorach.

Berthold LB 6411 / LB 123: Miernik oparty na moderatorze polietylenowym z licznikiem BF₃ lub He-3. Zakres energii: termiczne–10 MeV. Stosowany powszechnie w europejskich elektrowniach jądrowych.

TEPC (Tissue Equivalent Proportional Counter): Komora proporcjonalna wypełniona gazem tkankopodobnym (propan lub izobutan z domieszką CO₂). Mierzy rozproszenie energii liniowe (LET) bezpośrednio, co pozwala wyznaczać zarówno dawkę pochłoniętą, jak i współczynnik jakości bez znajomości widma. Podstawowe narzędzie w dozymetrii lotów kosmicznych (ISS, ESA). Wysoka cena i złożoność obsługi ograniczają użycie naziemne.

Porównanie instrumentów polowych:

Instrument	Zakres E neutronów	Odczyt	Dokładność	Przenośność
Kule Bonnera (zestaw)	Termiczne – 1 GeV	Widmo + H*(10)	Najlepsza (~15%)	Niska (>20 kg)
WENDI-II	0,025 eV – 5 GeV	H*(10)	~20–30%	Dobra (~5 kg)
Rem-metr BF₃	Termiczne – 10 MeV	H*(10)	~30%	Dobra (~3 kg)
TEPC	Termiczne – GeV	D + Q → H	~20%	Umiarkowana (~8 kg)
Komora rekombinacyjna	Prędkie + wysokoenergetyczne	Widmo	Umiarkowana	Umiarkowana

Reakcje jądrowe w detekcji neutronów i przekroje czynne

Wykrywanie neutronów opiera się na reakcjach jądrowych z konwersją energii kinetycznej neutronu (i energii Q reakcji) na naładowane produkty. Kluczowe reakcje:

¹⁰B(n,α)⁷Li: Q = 2,31 MeV (stan podstawowy ⁷Li, 6%) lub Q = 2,79 MeV (stan wzbudzony ⁷Li, 94%). Cząstka alfa i jądro litu dzielą energię. Przy 2,79 MeV: α ma ~1,47 MeV, ⁷Li* ~0,84 MeV. Stosowany w BF₃ (gaz), licznikach borowych (warstwami B₄C), scyntylatorach borowych (szkło litowo-borowe). Przekrój czynny dla neutronów termicznych: σ = 3835 barn (ogromny). Reguła 1/v: σ ∝ 1/v → σ ∝ E^(−½).

³He(n,p)³H: Q = 0,764 MeV. Proton i tryt dzielą energię: p ~0,573 MeV, ³H ~0,191 MeV. Stosowany w licznikach He-3 pod ciśnieniem 2–10 atm. Przekrój dla neutronów termicznych: σ = 5327 barn — najwyższy z dostępnych. Excellent resolution, lowest detection threshold. Niedobór He-3 po 2008 roku ogranicza dostępność.

⁶Li(n,α)³H: Q = 4,78 MeV. Duża energia Q → łatwo odróżnić od gamma. Stosowany w scyntylatorach LiI(Eu) i szkle litowym (Li-6 wzbogacone). Przekrój σ = 940 barn dla neutronów termicznych.

¹H(n,n)¹H (elastyczne rozproszenie na protonach): Q = 0. Proton odrzutu niesie maksymalnie energię neutronu (przy zderzeniu czołowym). Dla neutronów szybkich (>keV) to dominujący mechanizm w organicznych scyntylatorach i licznikach protonowych. Brak pragu energetycznego — działa od najniższych energii (z mniejszą czułością).

Rozszczepienie ²³⁵U, ²³³U, ²³⁹Pu: fragmenty rozszczepienia mają masę ~100 u i energię ~100 MeV — łatwo detekcja. Komory szczelinowe (fission chambers) z cienkimi warstwami U-235. Stosowane do monitoringu strumienia w rdzeniu reaktora.

Typy detektorów neutronowych — tabela

Detektor	Reakcja	Energia termiczna σ	Zakres energii	FWHM	Zastos.
BF₃ (gaz)	¹⁰B(n,α)⁷Li	3835 b	Termiczne–pośrednie	n/d	Monitor reaktora
He-3 (gaz)	³He(n,p)³H	5327 b	Termiczne–pośrednie	~2% przy 0,764 MeV	Spektr. termicznych
LiI(Eu) scynty	⁶Li(n,α)³H	940 b	Termiczne	~5–8%	Kule Bonnera, portale
Szkło Li-6	⁶Li(n,α)³H	940 b	Termiczne	~10%	Kule Bonnera, monitoring
Plastikowy scynty	¹H(n,n)¹H (proton)	niskie dla szybkich	Prędkie (>0,5 MeV)	Słabe	Monitoring szybkich
Scynty stilben	¹H(n,n)¹H + n-γ discrim.	—	Prędkie	Umiarkowane	Spektr. n-γ mixed field
Detektor ⁶Li-ZnS	⁶Li(n,α)³H	940 b	Termiczne	~8%	Portale, monitoring
Komora szczelinowa	²³⁵U(n,f)	587 b	Termiczne–epitermiczne	Słabe	Monitoring rdzenia
Detektory aktywacyjne	Au, In, Co, Dy	Różne	Szerokie (różne rezonanse)	Brak w czasie pomiarów	Pasywna dosymetria

Neutronowe źródła kalibracyjne

Wzorcowanie spektrometrów neutronowych wymaga źródeł o dobrze znanych widmach:

Am-241/Be (AmBe): Reakcja ²⁴¹Am → α → ⁹Be(α,n)¹²C. Widmo ciągłe 0–11 MeV z garbem przy ~4 MeV. Aktywność Am-241: typowo 1–37 GBq, co daje strumień ~10⁶–10⁷ n/s·GBq (ok. 60 n/s na MBq Am-241). Zalety: długi T₁/₂ Am-241 = 432 lat, źródło trwałe. Wady: złożone widmo, nie monoenergetyczne.

⁲⁵²Cf (Cf-252): Rozszczepienie spontaniczne, prawdopodobieństwo rozszczepienia = 3,1%. Widmo Maxwellowskie o temperaturze T = 1,42 MeV, E_mean ≈ 2,1 MeV. Strumień: 2,3 × 10⁶ n/s na µg Cf-252. T₁/₂ = 2,645 lat — wymaga regularnego wzorcowania aktywności. Stosowany jako standard do wzorcowania detektorów neutronowych i jako źródło inicjujące start-up reaktorów.

Pu-Be (PuBe): Starsza wersja AmBe, stosowana w historycznych instalacjach. Widmo podobne do AmBe. Pu-239 ma T₁/₂ = 24 100 lat — długotrwały, ale problem radiologiczny i regulacyjny (materiał jądrowy). Nowe instalacje preferują AmBe.

Reakcje (p,n) i (d,n) z akceleratora: T(p,n)³He przy 14 MeV, D(d,n)³He przy 2,5 MeV, D+T(d,n)⁴He przy 2,5 MeV → 14,1 MeV (fuzja). Monoenergetyczne lub quasi-monoenergetyczne neutrony. Stosowane do wzorcowania funkcji odpowiedzi detektorów. Dostępne w laboratoriach z akceleratorami Van de Graaffa, cyklotronami lub generatorami neutronów DD/DT.

Zastosowania spektrometrii neutronowej

Ochrona radiologiczna przy reaktorach: Reaktor MARIA w Świerku i projektowane reaktory SMR generują pola neutronów o złożonym widmie. Pracownicy mogą otrzymywać dawki od neutronów termicznych, pośrednich i prędkich. Dokładne widmo jest konieczne do poprawnego wyznaczenia dawki skutecznej (współczynnik wagowy neutronów wᴿ zmienia się od 2,5 do 20 w zależności od energii).

Terapia hadronowa (Kraków, Gliwice): Leczenie nowotworów protonami i jonami węgla generuje wtórne neutrony, które mogą dawać dawkę poza obszarem leczonym. Spektrometria neutronowa (kule Bonnera) pozwala ocenić tę dawkę i zoptymalizować plan leczenia. IFJ PAN Kraków prowadzi prace w tym zakresie w ramach Centrum Cyklotronowego Bronowice (CCB).

Safeguards MAEA: Pasywna koincydencyjna zliczalność neutronów spontanicznie rozszczepialnych nuklidów (Pu-240, Cf-252) jest podstawową metodą nieinwazyjnej weryfikacji ilości plutonu w pojemnikach. IAEA stosuje detektory He-3 w geometrii HLNCC (high level neutron coincidence counter). Dryfowanie od reguły 1/v przy He-3 i zanik tej czułości jest jednym z powodów poszukiwania alternatyw dla He-3.

BNCT (Bor Neutron Capture Therapy): Leczenie nowotworów przez napromienianie termicznymi neutronami po podaniu boru-10 do tkanki rakowej. Reakcja ¹⁰B(n,α) zabija komórkę rakową. Charakteryzacja widma neutronów przy głowicy napromieniania jest krytyczna dla dawkowania. Prace nad BNCT prowadzono m.in. w BNC (Petten, Holandia) i aktualnie w reaktorach badawczych Japonii.

Geofizyka (neutronowe mierniki wodoru): Opisane osobno w artykule o spowolnieniu neutronów i miernikach wodoru.

Polska infrastruktura i badania

W Polsce spektrometria neutronowa prowadzona jest przez:

• IFJ PAN Kraków (Centrum Cyklotronowe Bronowice): pomiary neutronów w terapii hadronowej, kalibracja kul Bonnera na wiązce protonowej;
• NCBJ Świerk (reaktor MARIA): monitoring neutronów w pobliżu rdzenia, badania aktywacyjne, dozymetria personelu;
• Politechnika Warszawska (Zakład Dozymetrii): komory rekombinacyjne KW-1, SWP-1, praca badawcza M. Maciaka;
• AGH Kraków (Wydział Fizyki i Informatyki): pomiary pola neutronowego przy akceleratorach medycznych, kalibracje kul Bonnera;
• Centrum Fizyki i Techniki Jądrowej CLOR: monitoring spektrometry środowiskowy.

Polska jest uczestnikiem programu IAEA dla harmonizacji metod spektrometrii neutronowej i uczestniczyła w międzynarodowych porównaniach kul Bonnera organizowanych przez PTB (Niemcy) i NPL (UK). Wyniki polskich laboratoriów są porównywalne z europejskim standardem, co potwierdzono w raportach EURADOS (European Radiation Dosimetry Group).

Trzy przykłady rachunkowe

Przykład 1 — przeliczenie zliczeń na fluencję przez przekrój czynny He-3

Detektor He-3: cylindryczny, objętość czynna V = 500 cm³, ciśnienie p = 4 atm, temperatura T = 293 K. Gaz doskonały: n = pV/RT = (4×101325 × 500×10⁻⁶) / (8,314 × 293) = 202 650 × 5×10⁻⁴ / 2436 = 0,0416 mol. Liczba atomów He-3: N = 0,0416 × 6,022×10²³ = 2,50×10²² atomów.

Przekrój czynny dla neutronów termicznych (2200 m/s): σ = 5327 barn = 5327×10⁻²⁴ cm². Szybkość reakcji przy fluencji neutronów termicznych Φ = 1 n/cm²·s:
R = N × σ × Φ = 2,50×10²² × 5327×10⁻²⁴ × 1 = 0,133 zliczeń/s.

Czułość detektora: S = 0,133 zliczenia/s na 1 n/cm²·s. Przy pomiarze 200 zliczeń/s: Φ_thermal = 200 / 0,133 = 1504 n/cm²·s.

Uwaga: to jest czułość dla neutronów termicznych. Dla neutronów prędkich σ jest kilka rzędów niższy, więc ta sama fluencja daje o wiele mniej zliczeń — stąd konieczność moderatora.

Przykład 2 — widmo z kul Bonnera i problem unfolding

Kule Bonnera mają 7 konfiguracji (bez kuli, kule 2", 4", 6", 8", 10", 12"). Zmierzone zliczenia (syntetyczne, n/s):
M = [150, 400, 600, 750, 680, 500, 320]

Uproszczona macierz odpowiedzi R (7 konfiguracji × 5 grup energii: termiczne, epitermiczne, eV-keV, prędkie 0,1–1 MeV, prędkie >1 MeV):
R ≈ macierz 7×5 z wartościami wzrastającymi dla większych kul przy wyższych energiach.

Liczba równań (7) > liczba grup energii (5) → układ nadokreślony — można użyć MNK (metoda najmniejszych kwadratów). Rozwiązanie daje wektor fluencji Φ dla każdej grupy, skąd liczy się równoważnik dawki H = Σ h_i × Φ_i × dE_i.

Wniosek: problem jest dobrze postawiony tylko gdy R jest dobrze znana i kondycja macierzy jest dobra. Przy złym uwarunkowaniu nawet małe błędy w M mogą powodować duże oscylacje Φ.

*Przykład 3 — równoważnik dawki ambientowy H(10) z widma**

Widmo neutronów (po rekonstrukcji z kul Bonnera) ma składowe:

• Termiczne (< 0,5 eV): Φ₁ = 1×10⁴ n/cm²·s;
• Epitermiczne (0,5 eV–100 keV): Φ₂ = 2×10³ n/cm²·s;
• Prędkie (>100 keV–1 MeV): Φ₃ = 5×10² n/cm²·s;
• Wysokoenergetyczne (>1 MeV): Φ₄ = 1×10² n/cm²·s.

Współczynniki fluencja → H*(10) z ICRP 74 (tabele): h₁ = 9 pSv·cm², h₂ = 40 pSv·cm², h₃ = 300 pSv·cm², h₄ = 400 pSv·cm².

H*(10) = Σ hᵢ × Φᵢ = 9×10⁴ + 8×10⁴ + 15×10⁴ + 4×10⁴ = 36×10⁴ pSv/s = 360 nSv/s = 1296 µSv/h.

To przekracza działający limit 20 µSv/h dla pracowników narażonych (ALARA) — konieczne środki ochrony. Gdybyśmy znali tylko całkowitą fluencję (suma Φ₁+Φ₂+Φ₃+Φ₄ = 12 600 n/cm²·s) i użyli jednego środkowego h = 50 pSv·cm²: H*(10) = 630 000 pSv/s = 2268 µSv/h — błąd 75% z powodu nieznajomości widma.

Pytania otwarte

1. Dlaczego neutron o energii 1 eV ma dużo wyższy przekrój czynny na reakcję z He-3 niż neutron o energii 1 MeV? Jaka jest zależność σ(E) dla reakcji ³He(n,p)³H i jak wpływa to na projekt detektora?

2. W kalibracji kul Bonnera stosuje się źródło AmBe o złożonym widmie. Jak to wpływa na dokładność wyznaczenia funkcji odpowiedzi dla poszczególnych kul? Czy lepiej byłoby używać monoenergetycznych neutronów z akceleratora?

3. Algorytm MAXED (maksimum entropii) i algorytm GRAVEL (iteracyjny MNK z regularyzacją) mogą dać różne widma dla tych samych danych pomiarowych. Jakie kryteria pozwalają wybrać lepsze rozwiązanie? Czy wystarczy porównanie reszt ||R·Φ − M||?

4. Komora rekombinacyjna KW-1 daje osiem sygnałów. Czy to wystarczy do dokładnej rekonstrukcji widma neutronów o energiach od termicznej do 100 MeV? Jakie są główne ograniczenia rozdzielczości energetycznej?

5. Neutronowe mierniki dawki osobistej (PADC — poliwęglanowe detektory śladowe) są tanie i pasywne. Jakie jest ich ograniczenie w zakresie energii i w jaki sposób wpływa to na wiarygodność dozymetrii osobistej przy reaktorze?

6. Przy BNCT charakteryzacja wiązki neutronów jest krytyczna dla dawkowania. Jakie parametry widma neutronów są najważniejsze dla obliczenia dawki biologicznej w tkance nowotworowej? Dlaczego same kule Bonnera mogą być niewystarczające?

7. Safeguards MAEA używa pasywnej koincydencji neutronów do ilościowego oznaczania Pu-240 bez otwierania pojemnika. Pu-240 ma szybkie rozszczepienie spontaniczne emitujące neutrony. Jakie inne nuklidy interferują i jak je odróżnić od sygnału Pu-240?

8. Polska sieć monitoringu posiada detektory NaI(Tl) i HPGe, ale brakuje stałej sieci kul Bonnera. Jakie korzyści przyniosłoby dodanie spektrometrów neutronowych do sieci PAA i gdzie (geograficznie) byłoby to najważniejsze z perspektywy ochrony radiologicznej?

Podsumowanie dydaktyczne

1. Neutrony wykrywa się pośrednio przez reakcje jądrowe (¹⁰B, ³He, ⁶Li, ²³⁵U, H). Sygnał pochodzi od wtórnych cząstek naładowanych, nie od neutronu samego. Przekroje czynne dla neutronów termicznych (∝ 1/v) są o wiele większe niż dla szybkich.

2. Detekcja neutronów prędkich wymaga moderatora (polietylen, woda) do spowolnienia do zakresu termicznego, gdzie czułość jest wysoka. Moderator jest częścią funkcji odpowiedzi i musi być uwzględniony w kalibracji.

3. Kule Bonnera (BSS) są standardem spektrometrii neutronowej o szerokim zakresie energii. Różne kule mają różne funkcje odpowiedzi — zestaw pomiarów tworzy układ równań dla widma fluencji.

4. Unfolding (rekonstrukcja widma) jest problemem odwrotnym, źle uwarunkowanym. Rozwiązanie (widmo) zależy od funkcji odpowiedzi, algorytmu (GRAVEL, MAXED), założeń startowych i regularyzacji. Dwa algorytmy mogą dać różne widma — raport musi opisać metodę.

5. Komory rekombinacyjne (KW-1) i pasywne układy wielowarstwowe (SWP-1) są alternatywami dla klasycznych kul Bonnera, oferując równoległą akwizycję w jednym detektorze lub eliminację problemów z elektroniką impulsową.

6. Z widma fluencji oblicza się równoważnik dawki przez ważenie współczynnikami h(E) z ICRP 74. Nieznajomość widma może prowadzić do błędów dawki o 50–100%, bo wᴿ neutronów zmienia się od 2,5 do 20 w zależności od energii.

7. Polskie ośrodki (IFJ PAN, NCBJ, PW, AGH) mają kompetencje w spektrometrii neutronowej. Najważniejsze zastosowania krajowe: terapia hadronowa (CCB Kraków), monitoring MARII (Świerk) i dozymetria akceleratorów medycznych.

8. Kluczowa metareguła: spektrometr neutronowy nie „widzi" energii neutronów bezpośrednio. Zbiera odpowiedzi kilku konfiguracji na pole i rekonstruuje widmo przez matematykę. Wiarygodność wyniku zależy od jakości funkcji odpowiedzi, liczby konfiguracji i metody rekonstrukcji.

Dodatkowe materiały multimedialne

Ćwiczenia praktyczne

Pierwsze ćwiczenie: przygotować syntetyczną macierz odpowiedzi dla pięciu konfiguracji detektora i ośmiu grup energii. Wartości nie muszą pochodzić z rzeczywistego detektora; mają tylko odzwierciedlać szerokie, przesunięte maksima odpowiedzi dla różnych grubości moderatora.

Drugie ćwiczenie: zdefiniować sztuczne widmo Phi, na przykład z dominującą składową termiczną, słabszą składową pośrednią i garbem w zakresie MeV. Obliczyć M = R * Phi, a następnie dodać szum Poissona odpowiadający skończonej liczbie zliczeń.

Trzecie ćwiczenie: spróbować odtworzyć Phi metodą najmniejszych kwadratów bez ograniczeń. Sprawdzić, czy pojawiają się ujemne wartości fluencji albo oscylacje między sąsiednimi grupami energii. Następnie powtórzyć obliczenie z ograniczeniem Phi_i >= 0.

Czwarte ćwiczenie: dodać prostą regularyzację gładkości, na przykład karę za duże różnice między sąsiednimi grupami energii. Porównać trzy widma: prawdziwe, odtworzone bez regularyzacji i odtworzone z regularyzacją.

Piąte ćwiczenie: obliczyć umowny równoważnik dawki z odtworzonego widma, używając syntetycznych współczynników h_Phi(E). Student ma pokazać, że dwa widma o podobnej liczbie całkowitych neutronów mogą dawać różne wartości dawki, jeżeli różnią się udziałem neutronów szybkich.

Szóste ćwiczenie: zaprojektować specyfikację przyszłego kalkulatora do serwisu. Minimalny zestaw funkcji powinien obejmować edycję macierzy odpowiedzi, wybór widma testowego, generowanie zliczeń, rekonstrukcję metodą nieujemnych najmniejszych kwadratów, wykres reszt R * Phi - M oraz raport niepewności. Ćwiczenie ma być całkowicie numeryczne i nie wymaga żadnego źródła neutronów.

Przejdź do ćwiczenia interaktywnego

Powiązane artykuły

• Od licznika Geigera do spektrometru gamma: tor pomiarowy promieniowania jonizującego
• Scyntylatory, fotopowielacze i półprzewodnikowe detektory promieniowania
• Statystyka zliczeń promieniotwórczych: Poisson, Gauss, chi-kwadrat i błędy aparaturowe
• Czas martwy, pile-up i gubienie impulsów w torach zliczających
• Kalibracja i niepewność pomiaru w laboratorium jądrowym