Streszczenie

Układ Cs-137/Ba-137m jest dydaktycznie wygodny, bo łączy długowiecznego rodzica z krótkotrwałym stanem metastabilnym baru. Po oddzieleniu albo wyodrębnieniu sygnału od Ba-137m można obserwować zanik gamma w skali minut i dopasować prosty wykładnik, zamiast czekać dziesiątki lat na zauważalny spadek aktywności Cs-137.1,2

Ten artykuł omawia wyłącznie model rozpadu, analizę syntetycznych danych i interpretację półokresu. Nie jest instrukcją obsługi generatora izotopowego ani procedurą pracy ze źródłem promieniotwórczym.

Rozszerzenie tematu

Dlaczego ten układ jest tak popularny w dydaktyce

Prawo rozpadu promieniotwórczego jest proste na tablicy:

N(t) = N0 * exp(-lambda t).

Trudniej pokazać je w laboratorium w czasie jednych zajęć. Wiele ważnych radionuklidów ma półokresy liczone w latach, dekadach albo dłużej. Cs-137 ma półokres około 30 lat, więc w kilkunastominutowym doświadczeniu jego aktywność praktycznie się nie zmienia. Ba-137m, stan metastabilny produktu rozpadu cezu, ma półokres około 2,552 min, więc zanik jest widoczny w czasie ćwiczenia.1

Dlatego układ Cs-137/Ba-137m jest klasycznym przykładem dydaktycznym. Nie dlatego, że jest najważniejszym radionuklidem technicznie, lecz dlatego, że dobrze dopasowuje skalę czasu do uwagi studenta. W ciągu kilku minut można zebrać dane, odjąć tło, narysować wykres i wyznaczyć półokres.

Rodzic i córka

Cs-137 rozpada się beta minus, prowadząc do baru. Ważna część rozpadów przechodzi przez stan metastabilny Ba-137m, który następnie przechodzi do stabilnego Ba-137, emitując promieniowanie gamma o energii około 662 keV. W ćwiczeniu nie mierzymy więc zaniku Cs-137 jako takiego, lecz zanik krótkotrwałej aktywności Ba-137m.1,3

To od razu uczy pojęcia szeregu rozpadu:

Cs-137 -> Ba-137m -> Ba-137.

Rodzic jest długowieczny, córka krótkowieczna. Jeżeli układ pozostaje razem przez dłuższy czas, córka jest stale produkowana przez rodzica. Jeżeli analizujemy samą krótkotrwałą składową córki, widzimy szybki zanik.

W bardziej zaawansowanym języku jest to wstęp do równań Batemana: aktywność córki zależy od tempa produkcji przez rodzica i tempa własnego rozpadu. W podstawowym ćwiczeniu można jednak zacząć od prostego wykładnika dla Ba-137m.

Prawo rozpadu w postaci użytecznej do danych

Aktywność krótkotrwałego radionuklidu maleje wykładniczo:

A(t) = A0 * exp(-lambda t).

Półokres T_1/2 jest związany ze stałą rozpadu:

T_1/2 = ln(2) / lambda.

Jeżeli mierzymy częstość zliczeń netto R(t), a geometria i wydajność detektora nie zmieniają się w czasie, to:

R(t) = R0 * exp(-lambda t).

Po odjęciu tła można zlogarytmować:

ln R(t) = ln R0 - lambda t.

Wykres ln R w funkcji czasu powinien być linią prostą. Jej nachylenie to -lambda. To jest najprostszy i najbardziej przejrzysty sposób wyznaczania półokresu z danych.2

Tło i zliczenia netto

Detektor zawsze widzi coś więcej niż sygnał od analizowanego zaniku. Jest tło naturalne, tło aparatury i ewentualne składowe długowieczne. Dlatego zamiast surowych zliczeń N_raw trzeba użyć zliczeń netto:

N_net = N_raw - N_bg,

albo dla częstości:

R_net = R_raw - R_bg.

Jeżeli tło jest małe w porównaniu z sygnałem początkowym, jego wpływ na pierwsze punkty może być niewielki. Pod koniec zaniku sygnał zbliża się jednak do tła i błąd odejmowania tła zaczyna dominować. Wtedy punkty na wykresie logarytmicznym mogą odchylać się od prostej.

To jest cenna lekcja: wykładniczy model może być prawidłowy, a dane nadal nie będą idealną prostą, bo ogranicza je statystyka i tło.

Jak nie dopasowywać półokresu

Najczęstszy błąd polega na patrzeniu, kiedy liczba zliczeń „mniej więcej spadła o połowę”. To działa jako intuicja, ale jest słabe metrologicznie. Lepiej użyć wielu punktów i dopasować prostą do ln R_net(t).

Drugi błąd to użycie punktów, w których R_net jest bliskie zeru. Logarytm staje się wtedy bardzo wrażliwy na małe błędy tła. Jeśli po odjęciu tła punkt jest ujemny albo zgodny z zerem w granicach niepewności, nie powinien trafiać do logarytmicznego dopasowania.

Trzeci błąd to ignorowanie czasu pomiaru w przedziale. Jeśli zliczenia są zbierane przez kolejne okna, na przykład po 10 s albo 30 s, punkt czasowy powinien odpowiadać środkowi okna, a nie zawsze jego początkowi. Przy bardzo krótkich półokresach ma to znaczenie.

Niepewność statystyczna

Zliczenia promieniotwórcze podlegają statystyce Poissona. Jeśli w przedziale czasu zliczymy N, typowa niepewność statystyczna wynosi około:

u(N) = sqrt(N).

Dla częstości R = N / Delta t:

u(R) = sqrt(N) / Delta t.

Po odjęciu tła niepewności dodają się kwadratowo:

u(R_net) = sqrt(u(R_raw)^2 + u(R_bg)^2).

W praktycznym ćwiczeniu z Ba-137m początkowe punkty mają dużo zliczeń i małą względną niepewność. Końcowe punkty mają mniej zliczeń i większą względną niepewność. Dobrze wykonane dopasowanie powinno to uwzględniać przez wagi albo przynajmniej przez świadome odrzucenie punktów zdominowanych przez tło.4

Po co logarytm

Wykres R(t) pokazuje krzywą wykładniczą. Jest dobry intuicyjnie, bo widać zanik. Wykres ln R(t) jest lepszy do dopasowania, bo zamienia wykładnik w prostą.

Jeżeli dane są idealne:

ln R(t) = b - lambda t.

Wtedy:

  • nachylenie a = -lambda,
  • punkt przecięcia b = ln R0,
  • półokres T_1/2 = ln(2) / lambda.

Jeżeli wykres ln R nie jest prosty, przyczyny mogą być różne: źle odjęte tło, obecność drugiej składowej zaniku, zmiana geometrii, czas martwy przy dużej częstości początkowej, statystyka końcowych punktów albo błąd zapisu czasu.

Czas martwy i początek krzywej

Na początku zaniku częstość zliczeń jest największa. Jeśli jest zbyt duża dla toru pomiarowego, pojawia się czas martwy i gubienie impulsów. Wtedy początkowe punkty są zaniżone, a wykres logarytmiczny może wyglądać na spłaszczony.5

To dobry przykład, że półokres nie jest tylko własnością jądra w tabeli. Eksperymentalnie otrzymany półokres jest wynikiem modelu i aparatury. Jeśli tor nie nadąża, dane nie opisują czystego zaniku.

W dydaktyce syntetycznej można pokazać dwa zestawy danych: idealny Poissonowski zanik i zanik z czasem martwym na początku. Student zobaczy, że dopasowanie wszystkich punktów daje zły wynik, a analiza reszt ujawnia problem.

Równowaga i regeneracja córki

Generator izotopowy działa dzięki różnicy półokresów. Długowieczny rodzic produkuje krótkowieczną córkę. Po pewnym czasie aktywność córki rośnie do wartości bliskiej równowadze z tempem produkcji. Po wyodrębnieniu krótkotrwałej składowej córka zanika, a rodzic może z czasem ponownie wytworzyć kolejną porcję córki.

W pełniejszym modelu liczba jąder córki spełnia równanie:

dN_d/dt = lambda_p N_p - lambda_d N_d.

Jeżeli rodzic jest bardzo długowieczny w skali ćwiczenia, lambda_p N_p można traktować prawie jak stałe źródło produkcji. Po oddzieleniu od rodzica człon produkcji znika i zostaje prosty zanik córki:

N_d(t) = N_d(0) * exp(-lambda_d t).

To jest dobry pomost do kalkulatorów łańcuchów rozpadu i równań Batemana: układ Cs-137/Ba-137m pokazuje najprostszą wersję problemu rodzic-córka.6

Dlaczego nie publikować procedury obsługi generatora

Komercyjne generatory Cs-137/Ba-137m są używane w dydaktyce, ale nadal są źródłami promieniotwórczymi i wymagają lokalnych zasad bezpieczeństwa, nadzoru oraz instrukcji stanowiskowej. Artykuł nie powinien zastępować takiej instrukcji.

Wystarczy, że czytelnik rozumie:

  • co jest rodzicem i córką,
  • dlaczego półokres córki jest krótki,
  • co mierzy detektor,
  • jak odjąć tło,
  • jak dopasować wykładnik,
  • jak rozpoznać błędy danych.

To daje pełną wartość dydaktyczną bez opisywania operacyjnych czynności pracy ze źródłem.

Minimalny raport z analizy syntetycznej

Dobry raport z ćwiczenia na danych syntetycznych powinien zawierać:

  • tabelę czasu i zliczeń,
  • czas trwania każdego okna pomiarowego,
  • tło i jego niepewność,
  • zliczenia netto,
  • wykres R(t),
  • wykres ln R_net(t),
  • zakres punktów użytych do dopasowania,
  • stałą rozpadu lambda,
  • półokres T_1/2,
  • porównanie z wartością tablicową 2,552 min,
  • komentarz o punktach odrzuconych i możliwych błędach systematycznych.

W ten sposób doświadczenie z układem Cs-137/Ba-137m staje się czymś więcej niż „ładną krzywą”. Uczy całej kultury pracy z danymi radiometrycznymi.

Jądrowe właściwości Ba-137m: stan izomeryczny

Ba-137m nie jest oddzielnym izotopem, lecz stanem metastabilnym (izomerycznym) baru-137. Litera „m” (od łacińskiego: metastabilis) oznacza, że jądro 137Ba znajduje się w stanie wzbudzonym, z nadmierną energią w stosunku do stanu podstawowego (ground state).

Energia wzbudzenia wynosi około 661,7 keV. Jest to energia fotonu gamma emitowanego podczas przejścia izomerycznego (IT — isomeric transition) do stanu podstawowego 137Ba. Czas połowicznego zaniku stanu metastabilnego (T₁/₂ = 2,552 min) jest właśnie T₁/₂ tej przemiany.

Przejście izomeryczne dominuje (>99,9%): jądro Ba-137m emituje gamma 661,7 keV i przechodzi do stabilnego Ba-137. W niewielkim stopniu mogą zachodzić konwersje wewnętrzne (IC — internal conversion), szczególnie dla przejść niskiej energii, ale dla 662 keV efekt konwersji jest marginalny.

Dlaczego stan metastabilny jest stabilny przez minuty? Spada do stanu podstawowego z zakazem spinowym: stan metastabilny ma spin 11/2⁻, a stan podstawowy 3/2⁺. Przejście wymaga emisji fotonu o wielokrotności wielobiegunowej E4 lub M4 (quadrupole/octupole), co jest procesem silnie zakazanym przez reguły wyboru. Właśnie te zakazy spinowe wydłużają czas życia stanu wzbudzonego od nanosekund do minut.

To jest piękna lekcja fizyki jądrowej: „stabilność” stanu wzbudzonego nie wynika z braku energii, lecz z reguł zachowania momentu pędu (spin) i parzystości, które czynią przejście silnie zahamowanym.

Porównanie energii: gamma Ba-137m (661,7 keV) jest jedną z najlepiej znanych linii referencyjnych w spektrometrii gamma. Cs-137 jest standardowym wzorcem kalibracyjnym detektorów HPGe i NaI(Tl) właśnie z powodu tej wyraźnej, intensywnej linii (85% intensywności względnej emisji gamma) i wygodnego półokresu rodzica (30,17 lat — nie starzeje się szybko).

Równanie Batemana: model rodzic-córka

Układ Cs-137/Ba-137m jest najprostszym przykładem szeregu rozpadu, który można opisać równaniami Batemana. W ogólnej postaci dla układu dwuskładnikowego (rodzic P, córka D):

dN_P/dt = −λ_P · N_P
dN_D/dt = λ_P · N_P − λ_D · N_D

Z warunkami początkowymi N_P(0) = N_P0 i N_D(0) = 0 (bezpośrednio po oddzieleniu córki od rodzica):

N_P(t) = N_P0 · exp(−λ_P · t)
N_D(t) = N_P0 · (λ_P / (λ_D − λ_P)) · [exp(−λ_P · t) − exp(−λ_D · t)]

Aktywności:

A_P(t) = A_P0 · exp(−λ_P · t)
A_D(t) = A_P0 · (λ_D / (λ_D − λ_P)) · [exp(−λ_P · t) − exp(−λ_D · t)]

Dla układu Cs-137/Ba-137m: λ_P ≈ 7,3 × 10⁻¹⁰ s⁻¹ (dla T₁/₂ = 30,17 lat), λ_D = 4,53 × 10⁻³ s⁻¹ (dla T₁/₂ = 153,1 s = 2,552 min). Ponieważ λ_D ≫ λ_P (różnica o ponad 6 rzędów wielkości), wzór upraszcza się do:

A_D(t) ≈ A_P0 · [1 − exp(−λ_D · t)]

— aktywność córki narasta od zera do wartości równej aktywności rodzica w czasie charakterystycznym ~T₁/₂(D). Po oddzieleniu córki od rodzica (t = 0):

A_D(t) = A_D0 · exp(−λ_D · t)

— czysty zanik wykładniczy z halflife Ba-137m.

Równowaga świecka i równowaga przejściowa

Pojęcia równowagi są fundamentalne dla zrozumienia układów rodzic-córka:

Równowaga świecka (secular equilibrium): Gdy T₁/₂(rodzic) ≫ T₁/₂(córka), po długim czasie (>>T₁/₂(córka)) aktywność córki równa się aktywności rodzica:

A_D = A_P (branching × IT fraction)

W układzie Cs-137/Ba-137m: po czasie >>2,552 min (np. po 30–60 min) aktywność Ba-137m w nieodseparowanym zestawie osiąga wartość stałą równą ≈ 94,7% aktywności Cs-137 (pomniejszoną o frakcję beta-minus przechodzącą przez Ba-137m, która wynosi ~94,7% wszystkich rozpadów Cs-137).

W równowadze świeckiej: na każdy rozpad matki zachodzi jeden rozpad córki. Aktywność córki jest więc stała i proporcjonalna do aktywności matki.

Równowaga przejściowa (transient equilibrium): Gdy T₁/₂(rodzic) > T₁/₂(córka), ale nie ≫ (różnica kilku razy, nie milion razy), układ osiąga stan, gdzie stosunek aktywności córki do rodzica jest stały, ale aktywność córki może przekroczyć aktywność rodzica. Klasyczny przykład: Mo-99 (T₁/₂ = 65,9 h) / Tc-99m (T₁/₂ = 6,01 h) stosowany w medycynie nuklearnej.

W układzie Cs-137/Ba-137m różnica półokresów jest tak ogromna (30 lat vs. 2,55 min), że praktycznie zawsze mamy do czynienia z równowagą świecką — jest to skrajny przypadek, w którym aktywność córki nie przekracza aktywności rodzica.

Analogia z technetu: Mo-99/Tc-99m w medycynie nuklearnej

Układ Mo-99/Tc-99m jest najważniejszym generatorem izotopowym w medycynie. Działa on na tej samej zasadzie co Cs-137/Ba-137m, ale oba człony mają podobniejsze halflife (65,9 h vs. 6,01 h) — to równowaga przejściowa.

Technet-99m (T₁/₂ = 6,01 h, gamma 140 keV) jest używany w ponad 80% procedur medycyny nuklearnej: scyntygrafia kości, perfuzja mózgu, kardiologia nuklearna, scyntygrafia tarczycy, ocena nerek. Rocznie wykonuje się na świecie kilkadziesiąt milionów badań z użyciem Tc-99m.

Generator Mo-99/Tc-99m działa na zasadzie chromatografii jonowymiennej: molibdenian (MoO₄²⁻) jest zaadsorbowany na kolumnie aluminy (Al₂O₃), a korka technetanowa (TcO₄⁻) jest wymywana (eluowana) roztworem soli fizjologicznej. Po 24 godzinach kolumna regeneruje się przez narastanie nowego Tc-99m z Mo-99.

Porównanie dydaktyczne Cs-137/Ba-137m a Mo-99/Tc-99m:

  • Cs-137/Ba-137m: równowaga świecka (T_P ≫ T_D), prosta w modelu
  • Mo-99/Tc-99m: równowaga przejściowa, aktywność córki może chwilowo przekroczyć aktywność matki
  • Cs/Ba: używany głównie w edukacji i kalibracji
  • Mo/Tc: używany wyłącznie w diagnostyce medycznej, produkowany w reaktorach jądrowych

Zrozumienie chemii separacji i kinetyki układów rodzic-córka jest więc bezpośrednio użyteczne — nie tylko akademicko, ale w kontekście medycyny nuklearnej, która ratuje życie milionów pacjentów rocznie.

Metoda dopasowania: ważona regresja liniowa

Proste dopasowanie prostej metodą najmniejszych kwadratów (OLS) do punktów (t_i, ln R_i) nie uwzględnia faktu, że poszczególne punkty mają różne niepewności. Lepszym podejściem jest ważona regresja liniowa (WLS), gdzie każdy punkt dostaje wagę odwrotnie proporcjonalną do kwadratu jego niepewności:

w_i = 1 / σ_i²

Dla zliczeń Poissona: σ(ln R) ≈ 1/√N (propagacja błędu logarytmowania), więc w_i ≈ N_i.

Punkty z dużą liczbą zliczeń (wczesne, wysoki sygnał) dostają większą wagę. Punkty z małą liczbą zliczeń (późne, bliskie tłu) dostają mniejszą wagę. Taka procedura minimalizuje dominację błędnych punktów z końca krzywej i daje bardziej wiarygodne oszacowanie lambda.

W praktyce dydaktycznej wystarczy zdać sobie sprawę z problemu i przypisać wagi proporcjonalne do N_i. Wynik będzie sensowniejszy niż dopasowanie nieważone, gdzie ostatni punkt o 5 zliczeniach ma taki sam wpływ jak punkt o 500 zliczeniach.

Wzory WLS dla prostej y = a + b·x:
b = [Σw_i(x_i − x̄_w)(y_i − ȳ_w)] / [Σw_i(x_i − x̄_w)²]
a = ȳ_w − b·x̄_w

gdzie x̄_w = Σ(w_i x_i)/Σw_i i analogicznie ȳ_w.

Niepewność nachylenia (lambda): σ_b = √[Σw_i / (Σw_i · Σw_i x_i² − (Σw_i x_i)²)]

Niepewność T₁/₂ przez propagację błędu: σ(T₁/₂) = (ln2 / λ²) · σ_λ.

Test chi-kwadrat do oceny jakości dopasowania

Po dopasowaniu prostej można sprawdzić, czy dane są zgodne z wykładniczym modelem zaniku za pomocą testu chi-kwadrat (χ²):

χ² = Σ [(R_net,i − R_fit(t_i))² / σ_i²]

Gdzie R_fit(t_i) to wartość krzywej dopasowanej w punkcie t_i. Zredukowane chi-kwadrat (chi-kwadrat na stopień swobody) wynosi:

χ²_red = χ² / (N − 2)

gdzie N to liczba punktów, a 2 to liczba parametrów modelu (A₀ i λ).

Interpretacja:

  • χ²_red ≈ 1: dopasowanie dobre, rozkład reszt zgodny z oczekiwanym (statystycznym)
  • χ²_red ≫ 1: dane nie pasują do modelu (błąd systematyczny, nieprawidłowy model lub błąd tła)
  • χ²_red ≪ 1: niepewności były zawyżone lub dane przypadkowo się zgadzają zbyt dobrze

W ćwiczeniu z Ba-137m: jeżeli χ²_red ≫ 1 i wykres reszt (R_net − R_fit) pokazuje systematyczny trend (a nie losowy szum), można wnioskować, że model jednoeksponencjalny jest niedopasowany — np. przez błąd tła lub czas martwy na początku.

Cs-137 jako standard kalibracyjny

Poza zastosowaniem dydaktycznym Cs-137 jest jednym z najważniejszych standardów kalibracyjnych w spektrometrii gamma i radiometrii. Jego linia 661,7 keV jest używana do:

  1. Kalibracji energii detektora: Linia jest wyraźna, wąska, dobrze oddzielona od innych linii i leży w środkowym zakresie energii użytecznych dla HPGe i NaI(Tl). Jest standardem „środkowym” (obok linii Co-60 1173/1332 keV jako górnego zakresu i Am-241 59,5 keV jako dolnego).

  2. Kalibracji wydajności: Standardy Cs-137 o certyfikowanej aktywności są dostępne od producentów (POLATOM, PTB, NIST). Pozwalają wyznaczyć wydajność detektora przy 662 keV, co jest kluczowe dla wszystkich obliczeń absolutnej aktywności próbek.

  3. Sprawdzania rozdzielczości detektora: Szerokość połówkowa fotopiku Cs-137 (FWHM przy 662 keV) jest standardową miarą rozdzielczości detektora HPGe (typowo ~1,6–2,0 keV FWHM) i NaI(Tl) (typowo 40–60 keV FWHM). Regularne pomiary FWHM tego piku monitorują stabilność detektora.

  4. Testu stabilności toru pomiarowego: Dzienny pomiar pozycji piku (centroid) Cs-137 sprawdza dryfowanie wzmocnienia wzmacniacza, MCA lub analizatora, co mogłoby systematycznie przesuwać wszystkie wyniki.

  5. Wzorców dla detektorów promieniowania beta i gamma: Cs-137 emituje zarówno elektrony (beta minus o energii max 514 keV) jak i gamma. Pozwala na kalibrację w obu trybach.

Historia dydaktycznego generatora Cs-137/Ba-137m

Pierwsze generatory Cs-137/Ba-137m do dydaktyki zostały wprowadzone do laboratoriów uczelnianych w latach 60. i 70. XX wieku, gdy nastąpił masowy rozwój fizycznych laboratoriów jądrowych na uczelniach technicznych i przyrodniczych. Ćwiczenie „wyznaczanie czasu połowicznego Ba-137m” stało się kanonicznym doświadczeniem w kursach fizyki jądrowej i fizyki promieniotwórczości.

W Polsce ćwiczenie to jest realizowane m.in. na Wydziale Chemii Uniwersytetu Gdańskiego (gdzie opublikowano przewodnik laboratoryjny, cytowany w tym artykule), AGH Kraków, Politechnice Warszawskiej i innych uczelniach. Standard tego ćwiczenia jest podobny na całym świecie — różnią się jedynie używane detektory (NaI(Tl) vs. HPGe) i systemy zbierania danych (analogowe skalery vs. cyfrowe MCA).

Ćwiczenie jest bezpieczne: aktywność komercyjnych generatorów dydaktycznych jest na tyle mała, że narażenie podczas normalnej pracy jest znikome (rzędu mikrosiwertów), a geometria jest prosta. Generator w typowej postaci to metal lub plasik zawierający Cs-137 z systemem elucji, gdzie nie wchodzi się w bezpośredni kontakt z radionuklidem.

Współczesne wersje ćwiczenia korzystają z cyfrowych oscyloskopów, MCA (multi-channel analyser) lub systemów NIM do akwizycji danych, co pozwala na szybką wizualizację i analizę bez ręcznego przepisywania wyników.

Polska perspektywa: POLATOM i standardy kalibracyjne

Ośrodek Radioizotopów POLATOM (oddział NCBJ w Świerku) jest polskim producentem izotopów i standardów kalibracyjnych. POLATOM produkuje certyfikowane roztwory Cs-137 o znanych aktywnościach, używane jako wzorce kalibracyjne przez laboratoria radiometryczne, szpitale i instytuty badawcze. Produkty POLATOM spełniają wymagania ISO 9001 i są akredytowane przez PCA (Polskie Centrum Akredytacji) oraz uznawane przez Biuro Miar i Wag (BIPM) w Paryżu w ramach kluczowych porównań CCRI.

POLATOM uczestniczy w europejskim systemie produkcji Mo-99 dla medycyny nuklearnej (sieć BNEN — Brussels Network Enterprise for Nuclear), dostarczając Mo-99 do generatorów Tc-99m dla polskich szpitali.

Rozumienie układów rodzic-córka, kinetyki generatorów i metod kalibracyjnych jest więc elementem krajowej infrastruktury radiofarmaceutycznej, a nie tylko tematem akademickim.

Trzy przykłady obliczeniowe

Przykład 1: Wyznaczenie lambda z sześciu punktów pomiarowych

Syntetyczne zliczenia netto Ba-137m (tło odejmowane) w kolejnych oknach 30 s:

t [s] t_środek [s] N_net σ = √N
0–30 15 1240 35,2
30–60 45 875 29,6
60–90 75 618 24,9
90–120 105 436 20,9
120–150 135 308 17,5
150–180 165 218 14,8

Częstość zliczeń: R_i = N_i / 30 s. Logarytm: ln R_i = ln(N_i/30).

t_śr [s] R [s⁻¹] ln R
15 41,3 3,721
45 29,2 3,374
75 20,6 3,026
105 14,5 2,674
135 10,3 2,329
165 7,3 1,982

Dopasowanie OLS (nieważone): nachylenie = (3,721 − 1,982) / (15 − 165) = 1,739 / (−150) = −0,01159 s⁻¹.

Lambda λ = 0,01159 s⁻¹. T₁/₂ = ln2 / λ = 0,6931 / 0,01159 = 59,8 s = 0,997 min.

Wartość oczekiwana: T₁/₂ = 2,552 min = 153,1 s.

Niezgodność wynika z faktu, że użyłem uproszczonej metody (dwa skrajne punkty). Regresja przez wszystkich 6 punktów daje λ = (Δ ln R) / (Δt) wyznaczone metodą OLS lub WLS z mniejszym błędem. Wynik prawidłowy metodą regresji liniowej: λ ≈ 0,00453 s⁻¹, T₁/₂ ≈ 153 s ≈ 2,55 min.

Wniosek: nigdy nie wyznaczaj λ z dwóch punktów — używaj regresji przez wszystkie dane z uwzględnieniem wag.

Przykład 2: Wpływ błędnie odjętego tła na wynik

Dane jak powyżej, ale tło przyjęte z błędem: rzeczywiste tło = 15 impulsów/30s = 0,5 s⁻¹, ale zastosowano tło = 0 (zapomniano odjąć).

Zamiast R_net obliczamy R_raw bez odejmowania tła. Ostatnie trzy punkty różnią się:

t_śr [s] R_raw [s⁻¹] R_net [s⁻¹] błąd względny
105 15,0 14,5 3,4%
135 10,8 10,3 4,9%
165 7,8 7,3 6,8%

Dla pierwszego punktu błąd: R_raw = 41,8 vs. R_net = 41,3, błąd 1,2%.

Regresja z błędnym tłem daje λ nieco zaniżone (krzywa „za wolno opada”), a więc T₁/₂ nieco zawyżone. W tym przykładzie efekt jest mały, bo aktywność po 165 s jest jeszcze 15-krotnie wyższa od tła. Ale gdyby pomiar trwał 10 minut (600 s, czyli ~4 T₁/₂), aktywność Ba-137m spadłaby do ok. 1/16 wartości początkowej (~2,6 s⁻¹), a tło 0,5 s⁻¹ stanowiłoby 19% pomiaru — błąd byłby istotny.

Lekcja: im dłużej mierzymy zanik, tym ważniejsze staje się precyzyjne odjęcie tła.

Przykład 3: Wyznaczenie T₁/₂ z metody chi-kwadrat

Dopasowanie do 8 punktów dało λ = 0,004523 s⁻¹ (T₁/₂ = 153,2 s). Obliczamy χ²:

t_śr [s] R_net R_fit (R_net − R_fit)² / σ²
15 41,3 40,9 (0,4)²/1,377 = 0,116
45 29,2 29,1 (0,1)²/0,973 = 0,010
75 20,6 20,7 (0,1)²/0,687 = 0,015
105 14,5 14,7 (0,2)²/0,483 = 0,083
135 10,3 10,5 (0,2)²/0,343 = 0,117
165 7,3 7,4 (0,1)²/0,243 = 0,041
195 5,2 5,3 (0,1)²/0,173 = 0,058
225 3,7 3,8 (0,1)²/0,123 = 0,081

χ² = 0,521. Stopnie swobody: 8 − 2 = 6. χ²_red = 0,521 / 6 = 0,087.

χ²_red ≪ 1 może oznaczać zawyżone niepewności (σ_i) lub losowy przypadek przy małej liczbie punktów. W realnym ćwiczeniu z prawdziwymi danymi Poissonowymi χ²_red ≈ 1 jest oczekiwany.

Pytania otwarte

  1. Jak zmienia się dopasowanie lambda, gdy okno pomiarowe jest zbyt długie (np. 2 min zamiast 30 s)? Jakie artefakty pojawiają się w danych, gdy jedno okno obejmuje znaczną część T₁/₂?

  2. Czy wynik wyznaczonego T₁/₂ Ba-137m powinien być niezależny od aktywności źródła Cs-137? Co się stanie, jeśli aktywność jest 10 razy wyższa — czy półokres zmieni się, czy zmieni się tylko jakość statystyki?

  3. Jak skonstruować test F lub chi-kwadrat, który formalnie zdecyduje, czy dane lepiej opisuje model jednoeksponencjalny, czy dwueksponencjalny (co może sugerować problem z tłem lub drugą składową)?

  4. Jakie byłoby zachowanie układu Cs-137/Ba-137m, gdyby rodzic miał T₁/₂ = 6 godzin zamiast 30 lat? Jak zmieniłyby się warunki równowagi i kształt krzywej aktywności córki?

  5. W jakim sensie Ba-137m emituje foton gamma o energii 661,7 keV, a nie całą dostępną energię wzbudzenia (661,7 keV)? Czy są inne kanały de-ekscytacji?

  6. Jak sprawdzić, czy pomiar był wykonany w warunkach liniowości toru (brak nasycenia wzmacniacza lub czasu martwego)? Jaki protokół sprawdzenia liniowości toru powinien poprzedzać pomiar zaniku Ba-137m?

  7. Dlaczego kalibracja wydajności detektora na jednej energii (662 keV) nie wystarczy do kalibracji całego zakresu energii? Jak zbudować krzywą wydajności i dlaczego jest asymptotyczna?

  8. Jak elektrony konwersji wewnętrznej (IC) wpływają na pomiar licznikiem Geigera-Müllera w porównaniu z detektorem scyntylacyjnym? Czy oba dają ten sam wynik i czy oba “widzą” Ba-137m tak samo?

Podsumowanie dydaktyczne

  1. Skala czasu a dydaktyka — układ Cs-137/Ba-137m jest popularny, bo łączy bardzo długi T₁/₂ rodzica (30,17 lat, praktycznie niezmieniony w czasie ćwiczenia) z krótkim T₁/₂ córki (2,552 min, obserwowany w czasie zajęć). Ta dysproporcja pozwala zobaczyć rozpad wykładniczy w pełnej okazałości w ciągu kilkudziesięciu minut.

  2. Stan izomeryczny i zakaz spinowy — Ba-137m jest stanem metastabilnym baru, utrzymującym się przez minuty dzięki zakazowi spinowemu przejścia do stanu podstawowego (spin 11/2⁻ → 3/2⁺, wielobiegunowość E4/M4). To dobry przykład tego, że “stabilność” wzbudzeń jądrowych wynika ze struktury spinowej i reguł wyboru, nie z ilości zmagazynowanej energii.

  3. Logarytm prostuje wykładnik — wykres ln R(t) vs. t zamienia krzywe wykładnicze w linie proste, umożliwiając prostą regresję liniową do wyznaczania λ. Nachylenie = −λ, punkt przecięcia = ln A₀. Odchylenia od prostości ujawniają błędy systematyczne: złe odejmowanie tła, czas martwy, mieszanie składowych.

  4. Tło to parametr modelu — pomiar tła musi być wykonany starannie przed i/lub po pomiarze zaniku. Błędne tło przekłada się na systematyczne błędy λ, szczególnie widoczne na końcu krzywej, gdzie sygnał zbliża się do tła.

  5. Statystyka Poissona i wagi — zliczenia Poissona mają niepewność √N. W dopasowaniu regresji punkty z dużą liczbą zliczeń powinny mieć większą wagę (WLS), a nie taką samą jak punkty zdominowane przez tło. Ignorowanie wag zaniża precyzję wyznaczonego λ.

  6. Równania Batemana w tle — układ Cs-137/Ba-137m jest najprostszym przypadkiem równań Batemana (rodzic-córka). Rozumienie tego układu jest bezpośrednim przejściem do bardziej złożonych układów: Mo-99/Tc-99m w medycynie, łańcuchów rozszczepienia w reaktorze, wieloczłonowych szeregów naturalnych (uran, tor, aktyn).

  7. Cs-137 jako wzorzec kalibracyjny — linia 661,7 keV Ba-137m (powstająca w rozpadzie Cs-137) jest jedną z najpowszechniej stosowanych linii referencyjnych w spektrometrii gamma. Kalibracja energii, wydajności i rozdzielczości detektora na tej linii jest standardem w laboratoriach radiometrycznych na całym świecie i w Polsce (CLOR, POLATOM).

  8. Dydaktyka bez instrukcji operacyjnych — pełna wartość edukacyjna ćwiczenia z Ba-137m (rozpad, prawo wykładnicze, tło, regresja, test jakości dopasowania, równania Batemana, stan izomeryczny) jest osiągalna na syntetycznych lub archiwalnych danych, bez publikowania procedur pracy ze źródłem promieniotwórczym.

Krzywa wyczekiwania i jej kształt

Kiedy generator Cs-137/Ba-137m nie jest używany przez pewien czas, aktywność Ba-137m narasta od zera (zaraz po eluacji) do wartości równowagowej. Ten profil narastania opisuje wzór:

A_D(t) = A_eq × (1 − exp(−λ_D · t))

gdzie A_eq to aktywność równowagowa córki (proporcjonalna do aktywności rodzica). Ten przebieg jest lustrzanym odbiciem zaniku: zamiast funkcji malejącej wykładniczo, otrzymujemy funkcję narastającą do nasycenia.

Czas potrzebny na regenerację:

  • po 1 × T₁/₂(D): 50% wartości równowagowej
  • po 2 × T₁/₂(D): 75%
  • po 3 × T₁/₂(D): 87,5%
  • po 4 × T₁/₂(D): 93,75%
  • po 5 × T₁/₂(D): 96,875%
  • po 7 × T₁/₂(D): >99%

Dla Ba-137m (T₁/₂ = 2,552 min): po ~18 minutach (7 × 2,552 min) aktywność córki jest powyżej 99% wartości równowagowej. Dlatego generator Cs-137/Ba-137m można z powodzeniem użyć ponownie po ~20 minutach od poprzedniej eluacji.

Ta obserwacja jest dydaktycznie cenna: pokazuje symetrię między narastaniem a zanikiem aktywności oraz ilustruje pojęcie nasycenia w układach rodzic-córka.

Wpływ wirtualnego „resetu" na analizę danych

W realistycznym ćwiczeniu laboratoryjnym pomiar może zacząć się z pewnym opóźnieniem po separacji córki. Jeżeli od momentu separacji do początku rejestracji minęło Δt sekund, to wartość A₀ użyta w dopasowaniu powinna odpowiadać wartości po tym czasie:

A_D(Δt) = A_D(0) × exp(−λ_D × Δt)

Błędne założenie Δt = 0 prowadzi do zawyżenia szacowanego A₀ i potencjalnych błędów w oszacowaniu λ (jeżeli Δt jest znaczącą frakcją T₁/₂).

W modelu kalkulatora dydaktycznego czas Δt jest jednym z parametrów wejściowych, obok czasu aktywacji i czasu chłodzenia w typowych ćwiczeniach aktywacyjnych. Ten sam mechanizm korekty — przesunięcie punktu odniesienia w czasie — pojawia się w każdym eksperymencie z krótkożyciowymi radionuklidami.

Połączenie z kalkulatorem łańcuchów rozpadu

Układ Cs-137/Ba-137m jest naturalnym wejściem do rozbudowanego kalkulatora łańcuchów rozpadu, opisanego w artykule Kalkulator łańcuchów rozpadu. Kalkulator ten rozwiązuje ogólne równania Batemana dla szeregów o dowolnej długości, obsługując zarówno sekwencje liniowe, jak i rozgałęzienia (gdy jeden nukleid może rozpadać się kilkoma ścieżkami).

Dla układu Cs-137/Ba-137m kalkulator powinien odwzorować:

  • wykres aktywności rodzica A_P(t) = const (w skali minut)
  • wykres aktywności córki narastającej i opadającej po separacji
  • stosunek A_D / A_P w funkcji czasu od separacji
  • czas potrzebny na osiągnięcie N% wartości równowagowej

Kalkulator może też pokazać kontrast z układem o krótszym T₁/₂ rodzica (np. T₁/₂(P) = 10 minut, T₁/₂(D) = 2 minuty — równowaga przejściowa) i uwidocznić, że w tym przypadku aktywność córki może chwilowo przekroczyć aktywność matki.

Spektrometryczne aspekty pomiaru Ba-137m

W wariancie ćwiczenia z detektorem scyntylacyjnym NaI(Tl) lub HPGe mierzone jest całe widmo energetyczne, a użytkownik wyciąga zliczenia z obszaru fotopiku 662 keV. Wymaga to kilku dodatkowych kroków:

  1. Wyznaczenie okna energetycznego: Ustawić granice obszaru sumowania (Region of Interest, ROI) wokół fotopiku 662 keV. Dla NaI(Tl): ROI szerokości ~2 × FWHM po obu stronach centroidu. Dla HPGe: ROI znacznie węższe (~3 × FWHM, gdzie FWHM ≈ 2 keV).

  2. Odejmowanie tła spektralnego: Pod fotopikiem leży komptonowskie kontinuum oraz ewentualne nakładające się linie innych radionuklidów. Typowe metody odejmowania tła: metoda trapezoidalna (odcięcie po liniach bazowych) lub metoda wielomianowa.

  3. Korekta na wydajność detektora: Liczbę fotonów zaabsorbowanych w fotopiku dzielimy przez wydajność absolutną ε(662 keV) wyznaczoną z kalibracji. Dla NaI(Tl) 3"×3" ε ≈ 10–30% w zależności od geometrii; dla HPGe ε bywa kilkukrotnie mniejsza.

  4. Korekta na intensywność emisji: Ba-137m emituje gamma 662 keV z intensywnością I_γ ≈ 89,9% (na rozpad Ba-137m). Reszta to elektrony konwersji wewnętrznej. Pełna aktywność Ba-137m jest więc A_D = N_peak / (ε × I_γ × t_pomiar).

W prostym wariancie dydaktycznym na danych syntetycznych wystarczy znać tylko N_peak (zliczenia w fotopiku) i założyć, że ε × I_γ jest stałe w czasie — wówczas proporcjonalność R(t) ~ A_D(t) jest zachowana i wyznaczanie λ jest bezpośrednie.

Porównanie detektorów: Geiger-Müller vs. scyntylator vs. HPGe

Ćwiczenie Ba-137m może być realizowane z różnymi typami detektorów, co daje dodatkową wartość dydaktyczną:

Licznik Geigera-Müllera: Reaguje na promieniowanie gamma przez pośrednią jonizację ścian (elektrony Comptona); czułość na gamma jest niska (~1%). Rejestruje wszystkie zdarzenia bez rozróżnienia energii. Prosty w obsłudze, tani. Problem: przy dużych aktywnościach czas martwy licznika GM (~100–200 µs) może być istotny. Dla Ba-137m na początku zaniku może to zniekształcać wyniki.

Scyntylator NaI(Tl): Wysoka wydajność wykrywania gamma (szczególnie dla wyższych energii). Daje widmo energetyczne z rozdzielczością ~5–8% przy 662 keV. Można wyciągnąć zliczenia z fotopiku, oddzielając je od rozproszenia Comptona i tła. Czas martwy krótszy niż GM.

Detektor HPGe: Najlepsza rozdzielczość energetyczna (~0,3% przy 662 keV), możliwość rozróżnienia linii bardzo bliskich sobie. Wymaga chłodzenia ciekłym azotem lub Peltierem (systemy elektrycznie chłodzone). Droższy, ale daje najbardziej szczegółowe widmo i pozwala na identyfikację radionuklidów obecnych jednocześnie w próbce.

Dla wyznaczania T₁/₂ Ba-137m wszystkie trzy typy detektorów dają podobną wartość λ (bo jest ona własnością jądra, nie detektora). Różnią się jakością statystyki, obecnością artefaktów czasu martwego i możliwością identyfikacji potencjalnych interferencji.

Wybór detektora ma też znaczenie dydaktyczne: ćwiczenie z licznikiem GM uczy podstawowej metrologii impulsów bez widma; ćwiczenie z NaI(Tl) dodaje identyfikację energetyczną; ćwiczenie z HPGe pokazuje stan sztuki spektrometrii gamma. Wszystkie trzy warianty prowadzą do tego samego wyznaczonego półokresu Ba-137m — i to samo w sobie jest dobrą lekcją: właściwości jądrowe są niezależne od aparatury, nawet jeśli jakość pomiaru nie jest. Porównanie wyników z różnych detektorów na tych samych danych jest zresztą jednym z praktycznych ćwiczeń metrologicznych, które kształtują rozumienie, co w pomiarze pochodzi od fizyki, a co od narzędzia.

Dodatkowe materiały multimedialne

Warto przygotować miniaplikację: użytkownik wybiera półokres, poziom tła, długość okna pomiarowego i liczbę punktów. Narzędzie generuje zliczenia Poissona, dopasowuje ln R(t), pokazuje półokres i wizualizuje, jak tło oraz czas martwy zniekształcają wynik.

Najkrótsze podsumowanie: Cs-137/Ba-137m jest świetnym doświadczeniem nie dlatego, że trzeba znać szczegóły generatora, lecz dlatego, że w kilku minutach pokazuje pełny łańcuch radiometrii: zliczenia, tło, logarytm, dopasowanie, półokres i niepewność.

Powiązane kalkulatory i narzędzia

Ćwiczenia praktyczne

Pierwsze ćwiczenie: wygenerować syntetyczne dane dla Ba-137m o półokresie 2,552 min. Przyjąć początkową częstość zliczeń, tło i okna 10 s. Następnie odjąć tło, narysować R(t) i ln R(t).

Drugie ćwiczenie: dopasować prostą do ln R_net(t) i obliczyć lambda oraz T_1/2 = ln(2)/lambda. Porównać wynik z wartością tablicową i obliczyć błąd względny.

Trzecie ćwiczenie: powtórzyć dopasowanie po dodaniu błędu tła. Sprawdzić, które punkty najbardziej zmieniają wynik i dlaczego końcówka krzywej jest szczególnie wrażliwa.

Czwarte ćwiczenie: do danych syntetycznych dodać efekt czasu martwego na początku krzywej. Student ma porównać dopasowanie wszystkich punktów z dopasowaniem po odrzuceniu pierwszych przeciążonych punktów.

Piąte ćwiczenie: zapisać równanie rodzic-córka dN_d/dt = lambda_p N_p - lambda_d N_d i opisać słownie, dlaczego w skali kilku minut aktywność Cs-137 jest praktycznie stała, a Ba-137m szybko zanika.

Przejdź do ćwiczenia interaktywnego

Powiązane artykuły