Streszczenie

Analiza aktywacyjna neutronami jest pomiarem po fakcie: najpierw stabilne jądra w próbce pochłaniają neutrony i zamieniają się w radionuklidy, a dopiero później mierzy się promieniowanie tych produktów. Z pozoru jest to prosta metoda aktywacji i odczytu, ale jej sens fizyczny leży w równaniu łączącym przekrój czynny, strumień neutronów, liczbę jąder tarczy, czas aktywacji, czas chłodzenia i wydajność detekcji.1

Ten artykuł opisuje analizę aktywacyjną jako model metrologiczny i obliczeniowy. Nie jest instrukcją pracy ze źródłem neutronów; ćwiczenia są oparte na danych syntetycznych, tak aby można było zrozumieć metodę bez opisywania operacyjnej geometrii aktywacji ani procedur laboratoryjnych z materiałami promieniotwórczymi.

Rozszerzenie tematu

Co właściwie mierzy analiza aktywacyjna

W zwykłej spektrometrii gamma próbka jest już promieniotwórcza, a detektor rejestruje jej fotony. W analizie aktywacyjnej próbka może być początkowo stabilna. Informację tworzy się przez reakcję jądrową: neutron zostaje pochłonięty albo wywołuje reakcję, powstaje radionuklid, a jego późniejszy rozpad daje mierzalne promieniowanie beta lub gamma.1

Najprostszy przykład myślowy wygląda tak. Mamy pierwiastek X, którego izotop X-A pochłania neutron:

X-A + n -> X-(A+1)*.

Jeżeli powstały izotop jest promieniotwórczy, jego aktywność po napromienieniu mówi coś o liczbie jąder X-A, przekroju czynnym reakcji i fluencji neutronów. Metoda może więc służyć w dwóch kierunkach. Jeżeli znamy strumień neutronów, możemy oznaczać skład próbki. Jeżeli znamy próbkę i przekrój czynny, możemy wyznaczać strumień neutronów.

Dziunikowski i Kalita opisują metodę aktywacyjną właśnie w tym duchu: aktywność wytworzonego izotopu zależy głównie od przekroju czynnego reakcji i gęstości strumienia neutronów padających na tarczę.1 To zdanie jest ważniejsze niż konkretne przykłady materiałów. Analiza aktywacyjna jest pomostem między mikroskopową wielkością jądrową, czyli sigma, a makroskopowym wynikiem radiometrycznym, czyli zliczeniami w detektorze.

Przekrój czynny i strumień

Przekrój czynny sigma nie jest powierzchnią geometryczną jądra. Jest miarą prawdopodobieństwa zajścia określonej reakcji dla neutronu o danej energii. W analizie aktywacyjnej interesuje nas zwykle przekrój czynny na aktywację, na przykład wychwyt radiacyjny (n, gamma) albo inna reakcja prowadząca do radionuklidu możliwego do pomiaru.1

Strumień neutronów Phi opisuje, ile neutronów przechodzi przez jednostkową powierzchnię w jednostce czasu. W prostym, jednogrupowym modelu tempo produkcji aktywowanego nuklidu można zapisać:

P = N * sigma * Phi,

gdzie N oznacza liczbę jąder tarczy. W jednostkach trzeba pamiętać, że 1 barn = 10^-24 cm2, więc strumień podawany w n/(cm2 s) dobrze pasuje do przekroju czynnego w barnach po przeliczeniu.

Liczba jąder tarczy wynika z masy, składu izotopowego i masy molowej:

N = m * f * N_A / M,

gdzie m to masa pierwiastka lub związku po odpowiednim przeliczeniu stechiometrii, f to udział izotopu aktywowanego, N_A to stała Avogadra, a M to masa molowa. W praktyce błąd w m, f albo składzie próbki wchodzi do wyniku tak samo realnie jak błąd zliczeń.

Narastanie aktywności podczas aktywacji

Podczas napromieniania produkt aktywacji nie tylko powstaje, ale równocześnie się rozpada. Jeżeli N* oznacza liczbę jąder produktu promieniotwórczego, to najprostsze równanie ma postać:

dN*/dt = P - lambda * N*.

Po czasie aktywacji t_a aktywność produktu na końcu napromieniania wynosi:

A_0 = P * (1 - exp(-lambda * t_a)).

Po czasie chłodzenia t_s aktywność spada do:

A(t_s) = P * (1 - exp(-lambda * t_a)) * exp(-lambda * t_s).

Po podstawieniu P = N sigma Phi dostajemy klasyczny wzór:

A(t_s) = N * sigma * Phi * (1 - exp(-lambda * t_a)) * exp(-lambda * t_s).

Dziunikowski i Kalita podają równoważną postać dla aktywowanej próbki, używając liczby jąder w jednostce objętości, objętości próbki, przekroju czynnego, gęstości strumienia neutronów, czasu aktywacji i czasu studzenia.1 Ten wzór warto rozumieć jakościowo. Krótkie napromienianie daje aktywność prawie proporcjonalną do czasu. Bardzo długie napromienianie nie zwiększa aktywności bez końca, bo układ zbliża się do nasycenia, gdzie tempo produkcji równoważy tempo rozpadu.

Nasycenie i wybór półokresu

Czynnik 1 - exp(-lambda t_a) jest czynnikiem nasycenia. Dla t_a << T_1/2 można przybliżyć go przez lambda t_a, więc aktywność rośnie prawie liniowo z czasem aktywacji. Dla t_a rzędu kilku półokresów czynnik zbliża się do jedności i dalsze wydłużanie czasu daje coraz mniejszy zysk.

To prowadzi do praktycznej, ale nadal nieoperacyjnej reguły metrologicznej: dobry produkt aktywacji musi mieć półokres pasujący do pytania pomiarowego. Jeżeli półokres jest zbyt krótki, sygnał może zaniknąć zanim zostanie dobrze zmierzony. Jeżeli jest zbyt długi, aktywność narasta wolno, a metoda może być mało czuła w rozsądnym czasie. Jeżeli w próbce powstaje kilka radionuklidów, czasy aktywacji i chłodzenia decydują o tym, który sygnał dominuje w widmie.1

W źródłowych ćwiczeniach Dziunikowskiego i Kality pojawia się właśnie ten problem: różne produkty aktywacji mają różne półokresy i różne promieniowanie, więc warunki pomiaru dobiera się tak, aby sygnał interesującego nuklidu był możliwie czytelny.1 W artykule nie potrzebujemy jednak instrukcji wykonania takiego ćwiczenia. Wystarczy model: aktywność produktu jest filtrowana przez czas aktywacji, czas chłodzenia i selekcję linii pomiarowej.

Jak z aktywności otrzymać strumień neutronów

Jeżeli mierzymy aktywność produktu i znamy N, sigma, t_a, t_s oraz lambda, to z prostego modelu można wyznaczyć strumień:

Phi = A(t_s) / [N * sigma * (1 - exp(-lambda t_a)) * exp(-lambda t_s)].

To równanie jest dydaktycznie bardzo mocne, bo pokazuje sens każdego składnika:

  • większa liczba jąder tarczy daje większą aktywność,
  • większy przekrój czynny daje większą aktywność,
  • większy strumień neutronów daje większą aktywność,
  • dłuższa aktywacja pomaga tylko do nasycenia,
  • dłuższe chłodzenie zmniejsza sygnał,
  • krótszy półokres szybciej daje sygnał, ale też szybciej go traci.

Nie wolno jednak mylić aktywności z surową liczbą zliczeń. Detektor mierzy impulsy w określonej geometrii, w określonym czasie i w określonym zakresie energii. Jeżeli produkt emituje foton gamma o intensywności I_gamma, detektor ma wydajność epsilon(E), a pomiar trwa t_m, to liczba zliczeń netto w fotopiku jest w przybliżeniu:

N_peak = A_mid * I_gamma * epsilon(E) * t_m * C,

gdzie A_mid oznacza aktywność reprezentatywną dla czasu pomiaru, a C zbiera poprawki: geometrię, czas martwy, samopochłanianie, koincydencje, tło i inne składniki zależne od konkretnej aparatury.2,3

Strumień nie zawsze jest monoenergetyczny

Wzór N sigma Phi jest wygodny, ale zakłada, że można użyć jednego przekroju czynnego i jednego strumienia. W rzeczywistości neutrony mają rozkład energii, a przekrój czynny zależy od energii. Wtedy tempo produkcji trzeba zapisać jako całkę:

P = N * integral sigma(E) * phi(E) dE.

To jest ten sam typ problemu, który pojawił się w spektrometrii neutronowej. Detektor albo próbka aktywacyjna nie mierzy całego widma bezpośrednio. Mierzy splot widma neutronów z funkcją reakcji, czyli przekrojem czynnym. Jeżeli przekrój ma silny pik rezonansowy w pewnym zakresie energii, próbka może być bardzo czuła na ten zakres, a mało czuła na inny.

Dziunikowski i Kalita zwracają uwagę, że dla indu aktywacja zachodzi nie tylko od neutronów termicznych, ale także w zakresie nadtermicznym, gdzie występuje rezonans.1 To jest dobra lekcja ogólna: nazwa „pomiar strumienia termicznego” może ukrywać składowe spoza termicznego przedziału energii, jeżeli reakcja ma istotny przekrój poza tym zakresem.

Metoda kadmowa jako idea korekty

W dydaktyce często wprowadza się pojęcie różnicowego pomiaru termicznej i nadtermicznej składowej aktywacji. Kadm ma bardzo duży przekrój czynny na pochłanianie neutronów termicznych, więc pomiar z filtrem kadmowym i bez niego może rozdzielać wkład termiczny od pozostałych wkładów. Dziunikowski i Kalita opisują ten sens przez różnicę aktywności całkowitej oraz aktywności skorygowanej o wkład nadtermiczny.1

W artykule wystarczy wersja modelowa:

A_th = A_total - F_Cd * A_Cd,

gdzie F_Cd jest poprawką uwzględniającą to, że filtr nie jest idealnym przełącznikiem „termiczne/nie-termiczne”. Z punktu widzenia metrologii ważne są trzy wnioski. Po pierwsze, materiał filtrujący zmienia widmo neutronów, a nie tylko „włącza” lub „wyłącza” sygnał. Po drugie, wynik zależy od przekrojów czynnych w funkcji energii. Po trzecie, poprawka F_Cd jest częścią budżetu niepewności, a nie kosmetycznym dodatkiem.

Nie podajemy tu geometrii, sekwencji czynności ani parametrów pracy ze źródłem neutronów. Dla studentów serwisu istotna jest logika korekty widmowej i rachunek niepewności, które można w pełni przećwiczyć na danych syntetycznych.

Analiza aktywacyjna jako oznaczanie składu

Jeżeli strumień neutronów jest znany lub kontrolowany przez laboratorium wzorcujące, analiza aktywacyjna staje się metodą oznaczania pierwiastków. Po aktywacji próbki mierzy się charakterystyczne linie gamma produktów aktywacji. Energia linii pomaga rozpoznać radionuklid, a pole fotopiku po poprawkach pozwala oszacować ilość pierwiastka.

W takim zastosowaniu równanie rozwiązuje się nie względem Phi, lecz względem N albo masy pierwiastka:

N = A(t_s) / [sigma * Phi * (1 - exp(-lambda t_a)) * exp(-lambda t_s)].

Jeżeli reakcja dotyczy tylko jednego izotopu naturalnego, trzeba uwzględnić jego abundancję izotopową. Jeżeli próbka ma złożoną matrycę, dochodzą korekty samopochłaniania gamma, interferencje linii, aktywacja innych pierwiastków i tło. Z tego powodu realna NAA jest metodą bardzo czułą, ale nie jest wolna od klasycznej metrologii widmowej.

Analiza aktywacyjna jako pomiar pola neutronowego

Drugi kierunek jest odwrotny: próbka aktywacyjna jest znana, a szukamy strumienia neutronów. Takie próbki nazywa się czasem monitorami aktywacyjnymi albo foliami aktywacyjnymi. Ich zaletą jest to, że mogą być małe i nie wymagają aktywnej elektroniki w samym miejscu ekspozycji. Po ekspozycji odczytuje się ich aktywność w kontrolowanym torze radiometrycznym.

Metoda jest szczególnie pouczająca dla studentów, bo mierzy pole neutronowe inną drogą niż detektory aktywne. Zamiast impulsów w czasie rzeczywistym dostajemy zintegrowany efekt reakcji jądrowych w próbce. To dobrze pokazuje różnicę między strumieniem chwilowym, fluencją zintegrowaną i aktywnością produktu.

Ograniczeniem jest czas. Jeżeli pole neutronowe zmieniało się w trakcie ekspozycji, wynik aktywacyjny integruje te zmiany z wagą zależną od rozpadu produktu. Dla produktu o krótkim półokresie późniejsza część ekspozycji ma większy wpływ na aktywność końcową niż wcześniejsza, bo wcześniejszy produkt zdążył częściowo zaniknąć. Dla produktu długowiecznego ważniejsza jest całkowita fluencja.

Widmo gamma produktu aktywacji

W praktyce wiele analiz aktywacyjnych kończy się spektrometrią gamma. Trzeba wtedy wykonać te same kroki, które pojawiają się w artykule o spektrometrii gamma w praktyce: kalibracja energii, wybór fotopiku, odejmowanie tła, obliczenie pola netto, korekta wydajności, ocena niepewności i zapis warunków pomiaru.2

Produkty aktywacji mogą emitować kilka linii gamma. To bywa zaletą, bo kilka linii pozwala sprawdzić zgodność wyniku. Bywa też problemem, bo linie różnych radionuklidów mogą się nakładać. Detektor HPGe daje tu przewagę nad NaI(Tl), ponieważ węższe piki ułatwiają rozdzielanie interferencji. W prostych ćwiczeniach modelowych można jednak używać syntetycznych widm, w których pozycje pików i tło są znane.

Ważny szczegół: aktywność w czasie pomiaru też się zmienia. Jeżeli półokres jest porównywalny z długością pomiaru, nie wystarczy wziąć N_peak / t_m i traktować aktywności jako stałej. Trzeba uwzględnić zanik w trakcie pomiaru albo przeliczyć wynik do określonego czasu odniesienia.

Budżet niepewności

Analiza aktywacyjna jest dobrym przykładem metody, w której jeden ładny wzór ukrywa długi budżet niepewności. Składniki typowe to:

  • statystyka zliczeń w fotopiku i tle,
  • kalibracja wydajności detektora dla energii linii gamma,
  • intensywność emisji I_gamma,
  • masa próbki i udział izotopu tarczowego,
  • przekrój czynny albo przekrój uśredniony po widmie,
  • czas aktywacji, chłodzenia i pomiaru,
  • korekta zaniku w trakcie pomiaru,
  • samopochłanianie w próbce,
  • geometria źródło-detektor,
  • czas martwy i pile-up,
  • interferencje linii gamma,
  • poprawki widmowe, na przykład kadmowe lub rezonansowe.

Jeżeli wynik ma być użyty do wyznaczenia strumienia neutronów, niepewność przekroju czynnego i widma neutronów może dominować nad samą statystyką zliczeń. Jeżeli wynik ma oznaczać pierwiastek w próbce, duże znaczenie może mieć skład matrycy i zgodność wzorca z próbką.

Co metoda mówi o broni jądrowej i testach

Analiza aktywacyjna jest ważna także dla zrozumienia skutków wybuchów jądrowych, ale trzeba zachować właściwy poziom opisu. Neutrony z reakcji jądrowych mogą aktywować materiały konstrukcyjne, grunt, wodę, powietrze albo elementy aparatury. W efekcie radioaktywność po zdarzeniu nie pochodzi wyłącznie z produktów rozszczepienia, lecz także z radionuklidów powstałych przez wychwyt neutronów i inne reakcje.4

Dla serwisu najcenniejszy jest tu język metrologiczny. Aktywacja neutronowa tłumaczy, dlaczego skład otoczenia i widmo neutronów wpływają na późniejszy ślad radiologiczny. Nie trzeba do tego żadnej instrukcji operacyjnej. Wystarczy rozumieć, że aktywność końcowa jest ważoną sumą reakcji sigma(E) * phi(E) po wszystkich istotnych nuklidach i po czasie.

Minimalny model obliczeniowy

Najprostszy kalkulator dydaktyczny powinien przyjmować:

  • masę pierwiastka albo liczbę jąder tarczy,
  • udział izotopu aktywowanego,
  • przekrój czynny sigma,
  • strumień neutronów Phi,
  • półokres produktu,
  • czas aktywacji,
  • czas chłodzenia,
  • czas pomiaru,
  • energię i intensywność linii gamma,
  • wydajność detektora.

Wynikiem powinny być: aktywność na końcu aktywacji, aktywność po chłodzeniu, oczekiwana liczba zliczeń w fotopiku, niepewność statystyczna i wykres narastania oraz zaniku aktywności. Drugi tryb kalkulatora powinien odwracać równanie: z aktywności albo zliczeń netto wyznaczać strumień neutronów.

Taki kalkulator byłby naturalnym rozszerzeniem kalkulatora aktywacji oraz kalkulatora łańcuchów rozpadu. W wersji zaawansowanej można dodać widmo neutronów w grupach energii i przekrój sigma(E), aby student zobaczył różnicę między przekrojem jednogrupowym a całką po widmie.

Najkrótsze podsumowanie: analiza aktywacyjna nie jest tylko techniką „wywołania radioaktywności”. Jest precyzyjnym rachunkiem produkcji i zaniku, w którym każdy etap pomiaru zostawia ślad w wyniku.

Metoda kadmowa: separacja składowej termicznej i nadtermicznej

Neutron „termiczny” w sensie metrologicznym to neutron o energii porównywalnej z energią cieplną otoczenia (około 0,025 eV w temperaturze 20°C). W praktyce w reaktorze lub przy pracy ze spowalniaczem istnieje szerokie widmo energii, od frakcji meV do kilkunastu MeV dla neutronów szybkich. Dla wielu nukleonów przekrój czynny na wychwyt radiacyjny (n,γ) spada z energią proporcjonalnie do 1/v poza obszarami rezonansowymi. To sprawia, że termiczne neutrony dominują w aktywacji dla tych nukleidów. Ale w obecności rezonansów (silnych wzrostów sigma w określonym zakresie energii) wkład nadtermiczny może być znaczny, a nawet dominujący.

Kadm naturalny absorbuje neutrony termiczne z przekrojem czynnym około 2500 barnów dla neutronów 0,025 eV, natomiast jest praktycznie przezroczysty dla neutronów o energiach powyżej granicznej energii kadmu (około 0,5 eV, zwanej kadm-cutoff). Folię lub pojemnik kadmowy umieszcza się wokół próbki aktywacyjnej i mierzy aktywność z filtrem oraz bez filtra. Różnica (po poprawce na nieidealne tłumienie) daje składową termiczną. Aktywność z filtrem odpowiada składowej nadtermicznej i rezonansowej.

Współczynnik kadmowy CDR (ang. cadmium ratio) definiuje się jako stosunek aktywności bez filtru do aktywności z filtrem:

CDR = A_total / A_Cd.

Dla materiałów z przewagą aktywacji termicznej CDR jest duży (>10). Dla nukleidów z silnym rezonansem w zakresie nadtermicznym CDR może wynosić kilka lub nawet zbliżać się do jedności. Indium-115 ma CDR w typowych warunkach reaktorowych na poziomie 2–5, co wskazuje na istotny wkład rezonansowy (główny rezonans 115In(n,γ) leży przy 1,46 eV z przekrojem czynnym około 30 000 barnów).

Wyznaczanie składowej termicznej z CDR wymaga poprawki Westcotta, uwzględniającej fakt, że rozkład Maxwella termicznych neutronów i odchylenia sigma od prawa 1/v mogą zaburzyć prostą zależność CDR od stosunku składowych strumienia. Czynnik g Westcotta koryguje efektywny przekrój czynny dla konkretnej temperatury moderatora:

sigma_eff = g(T) * sigma_0,

gdzie sigma_0 jest przekrojem czynnym tabelarycznym dla energii 0,025 eV. Dla większości nukleidów g bliskie jest jedności w zakresie 20–100°C, ale dla pierwiastków z silnymi rezonansami bliskimi progu termicznego (np. 197Au) odchylenie może sięgać kilku procent.

Całka rezonansowa i jej znaczenie dla analizy aktywacyjnej

Całka rezonansowa I₀ jest miarą aktywacji przez nadtermiczną część widma neutronów. Definiuje się ją jako całkę przekroju czynnego po nadtermicznym widmie 1/E:

I0 = integral(E_Cd, inf) sigma(E) / E dE.

Dla złota (197Au) I₀ wynosi około 1550 barnów przy sigma_0 około 98,7 barnów, co oznacza, że aktywacja rezonansowa w typowym reaktorze jest porównywalna lub wyższa niż termiczna, gdy gęstość strumienia nadtermicznego jest duża.

W pomiarach absolutnych NAA lub przy wyznaczaniu strumienia za pomocą monitorów aktywacyjnych aktywność monitora opisuje się wzorem dwuskładnikowym:

A = N * Phi_th * [sigma_0 * g(T) + I0 * alpha * f] * (1 - exp(-lambda * t_a)) * exp(-lambda * t_s),

gdzie f = Phi_epi / Phi_th jest stosunkiem strumieni nadtermicznego do termicznego, a alpha opisuje kształt widma nadtermicznego (odchylenie od idealnej 1/E, typowo alpha od −0,05 do +0,1 w reaktorach powszechnych).

Nieznajomość f i alpha może prowadzić do systematycznego błędu w wyznaczonym strumieniu. Dlatego użycie dwóch monitorów o różnych stosunkach I₀/sigma₀ pozwala wyznaczyć oba parametry jednocześnie — jest to metoda podwójnego monitora, stosowana w standardowej dozymetrii reaktorowej.

Kluczowe materiały stosowane jako monitory aktywacyjne

Wybór monitora aktywacyjnego zależy od zakresu energii neutronów, wymaganej czułości, dostępności emiterów gamma o dobrze znanych parametrach jądrowych i łatwości obsługi chemicznej. Poniższa tabela zestawia najczęściej stosowane materiały:

Nukleid Reakcja T₁/₂ produktu sigma_th [barn] I₀ [barn] Linia gamma [keV] Główne zastosowanie
197Au (n,γ) 198Au 2,69 dni 98,7 1550 411,8 Strumień termiczny i epitermiczny
115In (n,γ) 116mIn 54,29 min 202 2630 1293 Strumień epitermalny, CDR
55Mn (n,γ) 56Mn 2,578 h 13,4 14,0 846,8 Strumień termiczny
59Co (n,γ) 60Co 5,272 lat 37,2 74 1173/1332 Aktywacja długoterminowa
23Na (n,γ) 24Na 14,96 h 0,533 0,31 1368 Strumień termiczny, czysta 1/v
51V (n,γ) 52V 3,75 min 4,92 2,7 1434 Szybki odczyt
186W (n,γ) 187W 23,9 h 37,9 480 685,7 Wychwyt epitermiczny

Złoto jest historycznie najczęściej stosowanym monitorem ze względu na dobre właściwości chemiczne (bierne, nie utlenia się), relatywnie długi półokres produktu (2,7 dnia), jedną dominującą linię gamma i dobrze zmierzony przekrój czynny. Jednak wysoki stosunek I₀/sigma₀ (~15,7) sprawia, że wynik pomiaru złotym foliami jest silnie zależny od widma nadtermicznego.

Mangan ma znacznie niższy stosunek I₀/sigma₀ (~1,04), co czyni go lepszym monitorem „czystego” strumienia termicznego. Sód (Na-23) jest jeszcze lepszym monitorem dla aktywacji 1/v, ale jego krótki półokres produktu (24Na: 14,96 h) wymusza szybki pomiar.

Zastosowania NAA w różnych dziedzinach nauki

Analiza aktywacyjna neutronami jest jedną z najczulszych metod wielopierwiastkowych, zdolną wykryć stężenia rzędu ng/g (ppb) lub nawet pg/g (ppt) dla sprzyjających pierwiastków. Nie wymaga roztwarzania próbki chemicznej przed aktywacją (w wariancie bez destrukcji, INAA — instrumental NAA), co eliminuje ryzyko kontaminacji z odczynników.

Geochemia i petrologia: NAA jest stosowana do oznaczania pierwiastków ziem rzadkich (REE) w skałach i minerałach. Lantanowce mają przekroje czynne na aktywację (n,γ) rzędu 1–60 000 barnów i tworzą produkty aktywacji z dobrymi liniami gamma. Wzorzec frakcjonowania REE (diagram pająka) pozwala identyfikować procesy magmowe, metamorficzne i hydrotermalne. Klasyczny przykład: oznaczanie Sm, Eu, Gd, Tb, Dy w meteorytach pozwala identyfikować pierwotny skład Ziemi i porównywać z chondrytem.

Archeologia i ochrona dziedzictwa: NAA ceramiki identyfikuje surowce gliniane przez wzorzec pierwiastków śladowych. Różne złoża gliny mają charakterystyczne odciski geochemiczne. Metoda pozwala przyporządkować fragmenty ceramiki do konkretnych centrów produkcji, co jest dowodem na szlaki handlowe w prehistorii. Obsydian i szkło wulkaniczne oznacza się tą metodą w archeologii Bliskiego Wschodu (wzorzec Bi, Rb, Zr, Ba).

Medycyna sądowa i kryminalistyka: NAA włosów i paznokci może ujawniać historię ekspozycji na metale ciężkie (As, Sb, Hg, Se). Wzrost włosów około 1 cm/miesiąc pozwala rekonstruować oś czasu ekspozycji. Metoda była stosowana do analizy włosów historycznych postaci.

Ochrona środowiska: Monitorowanie emisji przemysłowych, analiza pyłów zawieszonych PM2.5 i PM10, oznaczanie pierwiastków śladowych w wodzie i osadach. NAA pyłu atmosferycznego pozwala odróżniać źródła naturalne (gleba, morze) od antropogenicznych (transport, energetyka, przemysł ciężki) na podstawie wzorców pierwiastkowych.

Przemysł jądrowy: NAA jest stosowana do kontroli czystości materiałów paliwowych i konstrukcyjnych reaktorów. Wymagania na zawartość boru w graficie reaktorowym (poniżej 0,4 ppm) lub absorbentów neutronowych w cyrkonowych koszulkach paliwowych są weryfikowane m.in. metodą aktywacyjną.

Badania żywności: Oznaczanie Se, Zn, Cu, Fe, Mn, Co w produktach żywnościowych i suplementach. NAA jest metodą referencyjną przy walidacji innych technik (ICP-MS, AAS) ze względu na brak destrukcji matrycy.

Reaktory jako źródła neutronów i polska perspektywa

Polska ma kilkudziesięcioletnie doświadczenie w analitycznym NAA, głównie dzięki reaktorowi badawczemu MARIA w Świerku, zarządzanemu przez Narodowe Centrum Badań Jądrowych (NCBJ). Reaktor MARIA (typ kanałowy, moc 30 MW) dostarcza gęstości strumienia termicznego do 3 × 10¹⁴ n/(cm²·s) w kanałach pionowych, co odpowiada możliwościom wiodących europejskich reaktorów badawczych. Kanały poziome i stanowiska pneumatyczne pozwalają na szybki transport próbek bez przerywania pracy reaktora.

NCBJ prowadzi usługowe NAA dla klientów zewnętrznych, w tym próbek geologicznych, biologicznych i materiałów budowlanych. Działalność ta wpisuje się w krajową infrastrukturę metrologii jądrowej. Instytut Fizyki Jądrowej PAN (IFJ PAN) w Krakowie uzupełnia tę działalność w zakresie dozymetrii termoluminescencyjnej i oceny pola promieniowania.

Akademia Górniczo-Hutnicza (Wydział Fizyki i Informatyki Stosowanej) prowadzi zajęcia dydaktyczne z NAA w ramach laboratorium jądrowego, korzystając ze źródeł Am-Be lub reaktora MARIA jako źródeł neutronów do ćwiczeń aktywacyjnych.

Poza NCBJ Polska uczestniczy w europejskiej sieci reaktorów badawczych, umożliwiającej studentom i badaczom dostęp do zagranicznych reaktorów (ILL w Grenoble, FRM-II w Monachium, HFR w Petten). Reaktor ILL osiąga maksymalny strumień termiczny około 1,5 × 10¹⁵ n/(cm²·s) i jest jednym z najintensywniejszych źródeł neutronów w Europie, stosowanym m.in. do ultra-śladowych analiz geochemicznych i standardyzacji jądrowych.

Polskie laboratoria biorą udział w międzynarodowych testach biegłości organizowanych przez IAEA (w ramach programu ALMERA — Analytical Laboratories for the Measurement of Environmental Radioactivity), weryfikując jakość wyników krajowych.

Porównanie NAA z innymi metodami wielopierwiastkowymi

Analiza aktywacyjna nie jest jedyną wielopierwiastkową metodą śladową. Jej zalety i ograniczenia najlepiej widać w porównaniu z innymi powszechnie stosowanymi technikami:

Cecha NAA (INAA) ICP-MS ICP-OES XRF
Destrukcja matrycy Nie Tak Tak Nie
Granica wykrywalności ng/g – pg/g pg/g – ng/g ng/g – µg/g µg/g – mg/g
Liczba pierwiastków jednocześnie powyżej 40 powyżej 70 powyżej 60 powyżej 20
Podatność na kontaminację Niska Wysoka Wysoka Niska
Koszt pomiaru Wysoki (czas reaktora) Średni Niski Niski
Wymagana próbka Stała lub ciekła Ciekła Ciekła Stała
Analiza na miejscu (in-situ) Nie Nie Nie Przenośne urządzenia
Możliwość analizy C, O, N, H Ograniczona (PGNAA) Nie Nie Nie

Unikalna cecha INAA — brak destrukcji próbki — sprawia, że metoda jest preferowana przy analizie rzadkich próbek historycznych, meteorytów lub materiałów biologicznych. Z kolei dla rutynowych badań środowiskowych ICP-MS jest popularniejsze ze względu na niższy koszt i szybszy czas analizy.

Interferencje i poprawki w NAA

Wynik NAA jest wiarygodny tylko wówczas, gdy potrafi się zidentyfikować i skorygować interferencje. Najważniejsze typy:

Interferencje spektrometryczne: Dwa różne radionuklidy mogą emitować linie gamma o zbliżonych energiach, nierozróżnialne nawet przez detektor HPGe. Wymagana jest identyfikacja metodą wieloliniową — jeden nukleid emituje zwykle kilka linii w znanych stosunkach intensywności.

Interferencje jądrowe (reakcje na matrycy): Aktywacja pierwiastków matrycy może wytwarzać produkty, które emitują fotony o energiach nakładających się z linią oznaczanego nuklidu. Np. w próbce bogatej w sód aktywacja 23Na → 24Na daje linię 2754 keV; przy dużej zawartości sodu może to zakłócić oznaczanie innych pierwiastków.

Reakcje sekwencyjne: W bardzo wysokich strumieniach może zajść podwójna aktywacja: pierwiastek X → produkt Y → produkt Z. Produkt Z jest źródłem interferencji, jeżeli jego linie gamma pokrywają się z liniami oznaczanego nuklidu.

Samopochłanianie neutronów: Próbka zawierająca silne absorbenty neutronów (B, Cd, Gd) pochłania część strumienia przed dotarciem do jądra tarczowego. Wyznaczony strumień jest wtedy niższy od rzeczywistego, chyba że uwzględni się samoprzesłanianie. Korekcja wymaga znajomości składu chemicznego i wymiarów próbki.

Zmienna wydajność detektora: Linie gamma produktów aktywacji mogą leżeć w szerokim zakresie energii (100–2800 keV). Krzywa wydajności detektora HPGe musi być dobrze skalibrowana, bo błąd wydajności przechodzi wprost na wynik.

Trzy przykłady obliczeniowe

Przykład 1: Wyznaczanie masy złota w próbce przez aktywację neutronową

Złota folię o masie m = 5,0 mg napromieniano przez t_a = 2 h w polu neutronów o gęstości strumienia Phi = 1 × 10¹² n/(cm²·s). Po czasie chłodzenia t_s = 4 h zmierzono aktywność izotopu 198Au.

Dane jądrowe 198Au: T₁/₂ = 2,6946 dni = 64,67 h → lambda = ln2/64,67 h = 1,072 × 10⁻² h⁻¹. Przekrój czynny sigma₀ = 98,7 × 10⁻²⁴ cm². Obfitość izotopu 197Au = 100% (złoto jednoizotopowe). Masa molowa M = 196,97 g/mol.

Liczba jąder tarczy: N = 5,0 × 10⁻³ g × 6,022 × 10²³ mol⁻¹ / 196,97 g/mol = 1,528 × 10¹⁹ jąder.

Czynnik nasycenia: S = 1 − exp(−1,072 × 10⁻² × 2) = 1 − exp(−0,02144) = 0,02121.

Aktywność bezpośrednio po aktywacji: A₀ = N × sigma × Phi × S = 1,528 × 10¹⁹ × 98,7 × 10⁻²⁴ × 1 × 10¹² × 0,02121 = 32,0 Bq.

Czynnik zaniku: D = exp(−1,072 × 10⁻² × 4) = exp(−0,04288) = 0,9580.

Aktywność po chłodzeniu: A(t_s) = 32,0 × 0,9580 = 30,7 Bq.

Spostrzeżenia: (1) krótki czas aktywacji (2 h wobec T₁/₂ = 64,67 h) sprawia, że aktywność jest daleka od nasycenia — prawie proporcjonalna do t_a; (2) 4 godziny chłodzenia mają niewielki wpływ (tylko 4,2% straty), co dowodzi, że czas chłodzenia ma znaczenie przede wszystkim przy produktach o krótkim T₁/₂.

Przykład 2: Wyznaczanie strumienia neutronów metodą aktywacji manganu

Tarczę MnO₂ o masie m = 10 mg napromieniono przez t_a = 1 h. Po chłodzeniu t_s = 30 min zmierzono aktywność 56Mn równą A(t_s) = 8 500 Bq.

Dane jądrowe 56Mn: T₁/₂ = 2,578 h → lambda = ln2/2,578 = 0,2689 h⁻¹. Dane 55Mn: sigma_th = 13,4 × 10⁻²⁴ cm², udział izotopowy f = 100%, masa molowa MnO₂ = 86,94 g/mol, udział masowy Mn = 54,94/86,94 = 63,19%.

Masa manganu: m_Mn = 10 × 10⁻³ × 0,6319 = 6,319 × 10⁻³ g.

Liczba jąder 55Mn: N = 6,319 × 10⁻³ × 6,022 × 10²³ / 54,94 = 6,926 × 10¹⁹.

Czynnik nasycenia: S = 1 − exp(−0,2689 × 1) = 0,2355.

Czynnik zaniku: D = exp(−0,2689 × 0,5) = exp(−0,1345) = 0,8742.

Mianownik N × sigma × S × D = 6,926 × 10¹⁹ × 13,4 × 10⁻²⁴ × 0,2355 × 0,8742 = 9,281 × 10⁻⁴ × 0,2059 = 1,911 × 10⁻⁴ cm².

Strumień: Phi = 8500 Bq / 1,911 × 10⁻⁴ cm² = 4,45 × 10⁷ n/(cm²·s).

Spostrzeżenie: Mn-55 ma niski stosunek I₀/sigma₀ ≈ 1,04, więc wynik jest mało wrażliwy na kształt widma nadtermicznego. Jest to monitor idealny do pomiaru składowej termicznej.

Przykład 3: Budżet niepewności wyznaczonego strumienia

Na podstawie przykładu drugiego przeliczamy budżet niepewności dla wyznaczonego Phi:

Wielkość Wartość Niepewność względna
Aktywność A(t_s) 8500 Bq ±1,5% (statystyka zliczeń)
Masa próbki m 10 mg ±0,2% (waga analityczna)
Przekrój czynny sigma 13,4 b ±1,2% (wartość tabelaryczna)
Czas aktywacji t_a 1 h ±0,1%
Czas chłodzenia t_s 0,5 h ±0,3%
Wydajność detektora ε — ±2,0% (kalibracja HPGe)

Łączona niepewność względna Phi (suma w kwadracie):
u_rel(Phi) = sqrt(1,5² + 0,2² + 1,2² + 0,1² + 0,3² + 2,0²) = sqrt(2,25 + 0,04 + 1,44 + 0,01 + 0,09 + 4,00) = sqrt(7,83) ≈ 2,80%.

Zatem Phi = (4,45 ± 0,12) × 10⁷ n/(cm²·s) na poziomie ufności 68%.

Wniosek: dominującą składową niepewności jest kalibracja wydajności detektora (2,0%) i statystyka zliczeń (1,5%). Przekrój czynny (1,2%) jest na trzecim miejscu. Poprawa kalibracji detektora i zwiększenie statystyki zliczeń przyniosłyby największe korzyści.

Niestandardowe techniki aktywacyjne

Poza klasyczną INAA istnieje kilka wariantów rozszerzających możliwości metody:

RNAA (radiochemiczna NAA): Po aktywacji próbkę roztwara się chemicznie i izoluje oznaczany pierwiastek przed pomiarem gamma. Umożliwia usunięcie interferencji matrycowych i obniżenie tła. Kosztem jest ryzyko kontaminacji i dodatkowa praca chemiczna. Stosowana dla pierwiastków o słabych liniach gamma lub w silnie interferujących matrycach.

PGNAA (NAA z natychmiastowymi fotonami gamma): Mierzy się promieniowanie gamma emitowane natychmiast po pochłonięciu neutronu (w ciągu pikosekund), bez konieczności aktywacji i chłodzenia. Wymaga pracy bezpośrednio przy strumieniu neutronów. Szczególnie przydatna dla pierwiastków dających słabe lub nietrwałe produkty aktywacji: H, B, C, N. Zastosowanie: analiza węgla w cemencie lub azotu w materiałach budowlanych.

Opóźniona emisja neutronów (DNAA): Mierzy zamiast fotonów gamma emisje opóźnionych neutronów od nuklidów o T₁/₂ od sekund do minut. Szczególnie czuła na uran i tor. Stosowana do wykrywania materiałów rozszczepialnych i ich ilościowego oznaczania.

Cykliczna NAA (CNA): Krótkie, powtarzające się cykle aktywacja-pomiar-aktywacja umożliwiają analizę nuklidów o bardzo krótkich półokresach (sekundy). Automatyczne systemy pneumatyczne transportują próbkę do reaktora i z powrotem, a detektor pracuje synchronicznie. Stosowana w analizie nuklidów 28Al (T₁/₂ = 2,24 min) dla oznaczeń glinu.

Pytania otwarte

  1. W jaki sposób można oddzielić wkład reaktywacji jąder produktu (burn-out) od prostego równania narastania aktywności? Kiedy burn-out jest istotny i jak zmienia wzór aktywności?
  2. Jak zdefiniować „efektywny strumień termiczny” w sposób niezależny od użytego monitora, jeżeli każdy monitor ma inny stosunek I₀/sigma₀?
  3. Dlaczego foliom złotą i manganową używa się razem w zestawie monitorów — jakie dodatkowe informacje o widmie neutronów można z nich łącznie wyznaczyć?
  4. Jak samopochłanianie neutronów w dużych próbkach wpływa na wynik NAA i jakie korekty stosuje się w praktyce dla próbek o masie powyżej kilku gramów?
  5. Czy NAA może być stosowana do analizy lekkich pierwiastków (C, N, O, H)? Jakie są alternatywne metody i kiedy PGNAA ma przewagę?
  6. W jaki sposób czas transportu próbki z reaktora do detektora wpływa na wynik NAA i jak modelować tę niepewność dla produktów o T₁/₂ rzędu minut?
  7. Jak weryfikuje się poprawność wyznaczenia strumienia metodą aktywacyjną przy braku niezależnej metody referencyjnej?
  8. Jakie są ograniczenia stosowania NAA w warunkach pola neutronów o silnie nieizotropowym kącie padania (np. wiązka reaktorowa), i jak próbkę orientować względem wiązki?

Podsumowanie dydaktyczne

  1. Reguła produkcji i zaniku — aktywność produktu aktywacji zależy od sześciu czynników: liczby jąder tarczy, przekroju czynnego, strumienia neutronów, czynnika nasycenia (1−e^{−λt_a}), czasu chłodzenia (e^{−λt_s}) i czasu pomiaru. Każdy z tych czynników wnosi do wyniku zarówno informację, jak i źródło błędu.

  2. Dwukierunkowe zastosowanie — ten sam wzór można rozwiązywać względem strumienia neutronów (wyznaczanie pola neutronowego monitorami aktywacyjnymi) lub względem masy pierwiastka (analiza składu próbki przy znanych warunkach aktywacji). Wybór kierunku zależy od tego, co jest lepiej znane: strumień czy skład.

  3. Widmo neutronów ma znaczenie — w obecności rezonansów i widma nadtermicznego wynik aktywacji nie opisuje tylko strumienia termicznego. Pojęcia całki rezonansowej I₀, współczynnika kadmowego CDR i czynnika g Westcotta pozwalają skwantyfikować ten wkład i wymagają uwzględnienia w rzetelnej analizie.

  4. Metoda niedestrukcyjna — INAA nie wymaga roztwarzania próbki przed aktywacją, co eliminuje ryzyko kontaminacji z odczynników. Jest to kluczowa przewaga nad technikami chemicznymi (ICP-MS, ICP-OES) przy analizie rzadkich materiałów, artefaktów historycznych lub próbek biologicznych, które nie mogą być zniszczone.

  5. Budżet niepewności jest długi — typowy wynik NAA ma co najmniej dziesięć składowych niepewności. Dominującymi często są: kalibracja wydajności detektora, przekrój czynny tabelaryczny i statystyka zliczeń. Czas aktywacji i chłodzenia są zwykle mniej istotne, o ile halflife produktu nie jest bardzo krótki względem tych czasów.

  6. Zastosowania są szerokie — od geochemii i archeologii, przez ochronę środowiska i medycynę sądową, po kontrolę materiałów jądrowych. Wiele programów certyfikacji (NIST SRM, IAEA CRM) korzysta z NAA jako metody referencyjnej przy sprawdzaniu dokładności innych technik analitycznych, co potwierdza jej status metody wzorcowej.

  7. Polska infrastruktura — reaktor MARIA w Świerku umożliwia NAA przy strumieniu do 3 × 10¹⁴ n/(cm²·s). NCBJ i IFJ PAN prowadzą usługową i naukową analizę aktywacyjną. Uczestnictwo Polski w testach biegłości IAEA ALMERA potwierdza jakość wyników krajowych laboratoriów i wpisuje NAA w szerszy system metrologii jądrowej.

  8. Połączenie z innymi metodami — NAA jest najsilniejsza jako metoda uzupełniająca. INAA wyznacza pierwiastki trudne do oznaczenia przez ICP-MS (np. Au, Hf, REE) lub niemetale przez PGNAA (B, N, H). Połączenie NAA z ICP-MS lub ICP-OES daje kompletny obraz składu wielopierwiastkowego i wzajemnie weryfikuje wyniki, co jest standardem w geochemii i naukach forensycznych.

Dodatkowe materiały multimedialne

Warto przygotować interaktywny wykres narastania i zaniku aktywności. Użytkownik zmienia półokres, czas aktywacji i czas chłodzenia, a narzędzie pokazuje czynnik nasycenia, aktywność końcową i aktywność w czasie pomiaru.

Druga miniaplikacja powinna pokazywać odwrócenie równania aktywacyjnego: użytkownik podaje zliczenia netto, tło, wydajność detektora, intensywność linii gamma i dane próbki, a program wyznacza strumień neutronów wraz z budżetem niepewności.

Gotowe narzędzia do tej sekcji: wizualizacja osi czasu aktywacji, kalkulator strumienia neutronów z aktywacji i kalkulator aktywacji.

Ćwiczenia praktyczne

Pierwsze ćwiczenie: dla syntetycznej próbki obliczyć liczbę jąder tarczy N = m f N_A / M. Przyjąć masę, udział izotopowy i masę molową, a następnie sprawdzić, jak błąd masy i udziału izotopu zmieniają wynik.

Drugie ćwiczenie: przyjąć sigma, Phi, T_1/2, t_a i t_s, a następnie obliczyć A_0 oraz A(t_s). Narysować aktywność jako funkcję czasu aktywacji i zaznaczyć punkt 1, 2, 3 i 5 półokresów.

Trzecie ćwiczenie: wygenerować syntetyczne zliczenia w fotopiku dla znanej aktywności, wydajności detektora i tła. Odjąć tło, przeliczyć zliczenia na aktywność i porównać wynik z wartością wejściową.

Czwarte ćwiczenie: odwrócić równanie i wyznaczyć Phi z aktywności po chłodzeniu. Powtórzyć obliczenie dla trzech wartości półokresu produktu i opisać, kiedy wynik jest najbardziej wrażliwy na błąd czasu chłodzenia.

Piąte ćwiczenie: przygotować dwugrupowy model neutronów: grupa termiczna i nadtermiczna. Przyjąć dwa przekroje czynne sigma_th i sigma_epi, obliczyć aktywność całkowitą, a potem sprawdzić, jak błędne założenie „wszystko jest termiczne” zniekształca wyznaczony strumień.

Szóste ćwiczenie: zastosować model korekty kadmowej wyłącznie na danych syntetycznych. Przyjąć A_total, A_Cd i F_Cd, policzyć A_th = A_total - F_Cd * A_Cd, a następnie przeprowadzić propagację niepewności dla F_Cd.

Przejdź do ćwiczenia interaktywnego

Powiązane artykuły