Przejdź do narzędzia

ChainFinder buduje i wizualizuje sieć powiązań jądrowych między izotopami — odpowiada na pytanie: „jakie nuklidy można osiągnąć z danego punktu startowego przez rozpad, wychwyt neutronów lub rozszczepienie?" Wynikiem jest interaktywny graf, który pozwala szybko ocenić, czy interesujący nas produkt jest w zasięgu, jaką drogą powstaje i przez jakie pośredniki przechodzi.

Jego możliwe zastosowania wykraczają poza samą broń jądrową:

  • mapowanie łańcucha produkcji izotopu medycznego lub przemysłowego
  • identyfikacja źródła niespodziewanych zanieczyszczeń napromieniowanej tarczy
  • analiza ścieżek kumulacji aktynidów i produktów rozszczepienia w paliwie reaktorowym
  • dydaktyczne wizualizacje procesów jądrowych w reaktorze lub podczas składowania wypalonego paliwa

Uzupełnia ChainSolver — gdy wstępnie zidentyfikujemy sieć połączeń i wybierzemy interesującą ścieżkę, ChainSolver przeprowadza pełną symulację ewolucji stężeń w czasie z uwzględnieniem harmonogramu napromieniowania.1

Przeznaczenie

Typowe zastosowania: mapowanie łańcucha produkcji izotopu medycznego lub przemysłowego, identyfikacja źródła niespodziewanych zanieczyszczeń napromieniowanej tarczy, analiza ścieżek kumulacji aktynidów i produktów rozszczepienia w paliwie reaktorowym, dydaktyczne wizualizacje procesów jądrowych w reaktorze lub podczas składowania wypalonego paliwa.

Rodzaje uwzględnianych przejść

Każde przejście można włączyć lub wyłączyć niezależnie:

Typ przejścia Symbol Opis Warunek
Rozpad radioaktywny α, β⁻, EC, IT, n, p, SF Samorzutne przemiany jądrowe według tablicy nuklidów Zawsze — niezależnie od strumienia
Wychwyt neutronów (n,γ) Pochłonięcie neutronu → produkt Z, A+1; stan g/m1/m2 φ > 0
Rozszczepienie (n,f) Rozbicie na dwa ciężkie fragmenty; top-10 wydajności φ > 0 i σf > 0
Reakcje progowe (n,2n), (n,p), (n,α) Reakcje na neutronach szybkich Opcjonalnie

Pominięcie wychwytu i rozszczepień upraszcza graf do łańcucha rozpadu — przydatne przy badaniu czystego łańcucha rozpadu podczas składowania lub transportu. Pominięcie reakcji progowych jest uzasadnione w reaktorach termicznych, gdzie szybkie neutrony stanowią małą frakcję strumienia.

Audyt modelu i interpretacja ścieżek

Strona kalkulatora zawiera sekcję Audyt modelu, ponieważ sam graf nie powinien być czytany jako ranking prawdopodobieństwa. ChainFinder odpowiada przede wszystkim na pytanie: czy istnieje ścieżka połączeń jądrowych między nuklidami? Nie oznacza to jeszcze, że dana ścieżka jest dominująca w realnym reaktorze, tarczy napromienianej albo próbce po latach chłodzenia.

W interpretacji wyników trzeba rozdzielić:

  • ścieżkę formalną — przejście jest dozwolone przez dane rozpadu, wychwytu, rozszczepienia albo reakcję progową;
  • ścieżkę istotną — przejście ma duży branching ratio, przekrój czynny, wydajność rozszczepienia albo czas półrozpadu pasujący do rozważanego scenariusza;
  • ścieżkę praktyczną — warunki strumienia, energii neutronów i czasu napromieniania rzeczywiście pozwalają zgromadzić produkt w mierzalnej ilości.

Dlatego ChainFinder najlepiej traktować jako narzędzie rozpoznawcze przed obliczeniem ilościowym. Jeżeli graf wskazuje kilka możliwych dróg, ich wagę należy sprawdzić w ChainSolverze, który rozwiązuje układ równań w czasie.

Wynik ChainFindera zawiera tabelę audytu ścieżek. Dla każdego przejścia kalkulator próbuje opisać, czy jest ono mocne, słabe czy tylko formalne. Wykorzystuje do tego:

  • branching ratio dla rozpadów,
  • przekroje czynne dla wychwytu neutronów,
  • wydajności produktów rozszczepienia,
  • półokresy nuklidów pośrednich,
  • typ przejścia i jego sens fizyczny w danym kierunku wyszukiwania.

Taki ranking nie zastępuje rozwiązania równań Batemana. Ma jednak praktyczną wartość dydaktyczną: student widzi, że krótka ścieżka grafowa może być fizycznie słaba, jeśli zawiera bardzo mały branching ratio, śladowy yield albo pośrednik o półokresie niepasującym do scenariusza.

Algorytm BFS — dwa kierunki wyszukiwania

Wersja webowa implementuje BFS w dwóch trybach wybieranych przełącznikiem w formularzu:

→ Do przodu: co powstaje z X? BFS od nuklidu startowego — na każdym kroku szukamy produktów (córek) przez rozpad, wychwyt i rozszczepienie. Produkty rozszczepienia dodawane jako liście grafu, nie rozwijane dalej. Pytanie: „Co powstaje z Co-60 po 5 latach?"

← Wstecz: skąd pochodzi X? BFS od nuklidu docelowego — szukamy rodziców: nuklidów, które przez dowolne przejście dają cel jako produkt. Geometryczne odwrócenie każdego trybu, weryfikacja przez bazę danych. Pytanie: „Skąd pochodzi Cs-137 w napromieniowanym materiale?"

Odwracanie przejść (kierunek wstecz)

Dla każdego rodzaju przejścia geometryczne odwrócenie jest jednoznaczne:

Szukamy rodziców (Z, A) Kandydat geometryczny Weryfikacja
przez rozpad α (Z+2, A+4) forward α-decay → daje (Z,A)?
przez rozpad β⁻ (Z−1, A) forward β⁻ → daje (Z,A)?
przez rozpad EC (Z+1, A) forward EC → daje (Z,A)?
przez przejście IT (Z, A, s+1) lub (Z, A, s+2) forward IT → daje (Z,A,s)?
przez emisję n / p (Z, A+1) / (Z+1, A+1) forward → daje (Z,A)?
przez wychwyt (n,γ) SQL: sigmaCs WHERE z=Z, a=A−1, prod_state=s σ > 0?
przez rozszczepienie SQL: yieldTs/yieldFs WHERE product=cel Top-5 wg max(cum_yield)

Każdy geometryczny kandydat weryfikowany jest przez forward decayTransitions(): jeśli X→Y w trybie forward, to Y znajdzie X w trybie backward. Gwarantuje to spójność obu kierunków bez duplikowania logiki heurystyk.

Wyznaczanie trybów rozpadu

Tryby wyznaczane są z trzech źródeł w kolejności priorytetu:

  1. Dane bazy OOB (tabela sigma_a) — dla nuklidów z wpisem decay_type > 0
  2. SubBranchings — rozgałęzienia do stanów izomerycznych (g/m1/m2 córki)
  3. Heurystyka fizyki jądrowej — reguły empiryczne dla nuklidów bez danych:
    • Z ≥ 84, mały nadmiar neutronów → rozpad α (Z−2, A−4)
    • N > Nstab → β⁻ (Z+1, A)
    • N < Nstab → EC/β⁺ (Z−1, A)
    • Progi Z-zależne dla aktynidów: Z ≥ 92 toleruje większy nadmiar neutronów przy rozpadzie α

Wizualizacja grafu SVG

Wynik prezentowany jest jako interaktywny graf SVG renderowany przez JavaScript:

  • Węzły — prostokąty kolorowane wg dominującego trybu rozpadu (jak karta nuklidów NKE); klik otwiera szczegóły w NKE
  • Krawędzie — zakrzywione strzałki kolorowane wg typu przejścia: zielony = rozpad, niebieski = wychwyt, fioletowy = rozszczepienie
  • Grubość strzałki — proporcjonalna do branching ratio / wydajności rozszczepienia
  • Układ — BFS-poziomowy: nuklid startowy po lewej, kolejne generacje ku prawej
Graf łańcucha transmutacji wygenerowany przez ChainFinder — węzły kolorowane wg trybu rozpadu, strzałki wg typu przejścia
Graf łańcucha transmutacji wygenerowany przez ChainFinder — węzły kolorowane wg trybu rozpadu, strzałki wg typu przejścia

Przykład: łańcuchy produkcji U-234

Klasyczny przykład z dokumentacji — U-234 produkowany z U-238 lub U-235 w reaktorze uranowym wieloma równoległymi drogami:

Łańcuch trywialny (α):
U-238 →(n,2n) U-237 →β⁻ Np-237 →(n,γ) Np-238 →β⁻ Pu-238 →α U-234

Łańcuch przez β⁻ Th (krótki):
U-238 →β⁻ Th-234 →β⁻ Pa-234 →β⁻ U-234

Łańcuch przez α Pu (dłuższy):
U-235 →(n,γ) U-236 →(n,γ) U-237 →β⁻ Np-237 →(n,γ) Np-238 →β⁻ Pu-238 →α U-234

Fakt, że wszystkie te ścieżki zostaną znalezione automatycznie, potwierdza kompletność algorytmu. Uruchom wyszukiwanie dla U-234

Powiązane materiały

Streszczenie

ChainFinder to kalkulator sieci przejść jądrowych — narzędzie do odkrywania ścieżek transmutacji między izotopami. Artykuł opisuje funkcjonalność kalkulatora i wyjaśnia kontekst akademicki: teorię grafów w fizyce jądrowej, równania Batemana, kompletność sieci reakcji, algorytmy przeszukiwania i zastosowania praktyczne w reaktorach, medycynie i środowisku.

Rozszerzenie tematu

Sieć przejść jądrowych to jeden z najważniejszych pojęć w fizyce jądrowej stosowanej. Każdy reaktor, każda napromieniowana próbka, każde wypalione paliwo — jest miejscem, gdzie dziesiątki lub setki izotopów powiązanych siecią rozpadów i reakcji ewoluują jednocześnie w czasie. ChainFinder jest narzędziem do mapowania tej sieci — zanim przeprowadzi się obliczenia ilościowe (to robi ChainSolver), trzeba wiedzieć, jakie ścieżki w ogóle istnieją.

Teoria grafów a sieci nuklidów

Sieć nuklidów można naturalnie opisać jako skierowany graf (directed graph, digraph):

Węzły = nuklidy. Każdy izotop (stan jądrowy) jest węzłem. Dla izotopów z izomerami: osobny węzeł dla g, m1, m2. Łącznie baza ORIP_XXI zawiera 3 785 węzłów.

Krawędzie = przejścia. Każda możliwa przemiana (rozpad, wychwyt neutronu, rozszczepienie) jest krawędzią skierowaną: od rodzica do córki/produktu. Waga krawędzi = branching ratio, przekrój czynny, lub wydajność rozszczepienia.

Graf rozpadów. Jeśli uwzględniamy tylko samorzutne przemiany (bez neutronów), graf jest acykliczny skierowany (DAG) — nie ma cykli w rozpadach, bo każda przemiana zmniejsza energię wzbudzenia lub zmienia N/Z w kierunku stabilności.

Graf z reakcjami. Gdy dołączymy wychwyt neutronów i reakcje progowe, graf może zawierać cykle (np. A→B przez wychwyt, B→A przez rozpad) — technicznie nie jest już DAG. To komplikuje algebrę, ale ChainFinder obsługuje takie przypadki przez detekcję cykli w BFS.

Głębokość BFS a zakres sieci. Parametr "głębokości" w ChainFinder określa, ile kroków grafu maksymalnie zostanie eksplorowane od punktu startowego. Głębokość 3: do wnuków. Głębokość 6: sześć przejść. Dla grafu z 3 785 węzłami, głębokość 10 może teoretycznie objąć tysiące węzłów — dlatego parametr jest ograniczony.

Równania Batemana — dlaczego ChainFinder to nie koniec

ChainFinder pokazuje, CO jest możliwe. ChainSolver oblicza, ile KAŻDEGO nuklidu powstaje w czasie. Różnica to różnica między "mapą dróg" a "symulacją ruchu drogowego".

Równania Batemana (1910). Harry Bateman, angielski matematyk, w 1910 roku opublikował analityczne rozwiązanie dla liniowego łańcucha rozpadów (A→B→C→...→Z). Jego rozwiązanie zakłada izolowany łańcuch bez rozgałęzień i bez napromienowania.

Rozszerzenia Batemana. W rzeczywistości mamy:

  • Rozgałęzienia (izotop A rozpadająsię do B i C z różnymi branching ratios)
  • Reakcje z neutronami (wychwyt zmienia N, rozszczepienie → produkty rozszczepienia)
  • Zmienny strumień neutronów w czasie (napromieniowanie włączone/wyłączone)
  • Cykle (A→B→A przez różne ścieżki)

To wymaga ogólnego układu równań różniczkowych: dN_i/dt = ∑_j λ_ij N_j - λ_i N_i, gdzie λ_i = λ_rozpad + φ·σ_i. W formie macierzowej: dN/dt = A·N. Rozwiązanie: N(t) = exp(A·t)·N₀. Problem: macierz A jest rzadka, ale jej eksponenta macierzowa jest trudna obliczeniowo.

Metody numeryczne. ChainSolver implementuje kilka metod: Crank-Nicolson (niejawna, stabilna), Padé [2/2] (dokładna dla sztywnych układów), RK4, RKF45 (adaptacyjna). Wybór metody zależy od "sztywności" układu (ang. stiffness) — stosunku najkrótszego do najdłuższego T½ w sieci.

Dlaczego ChainFinder jest potrzebny przed ChainSolverem. ChainSolver rozwiązuje układ równań dla WSZYSTKICH nuklidów w sieci jednocześnie. Dla bazy 3 785 nuklidów: macierz A jest 3785×3785. Zbudowanie tej macierzy i rozwiązanie jest możliwe, ale wolne. ChainFinder pozwala wybrać podzbiór nuklidów istotnych dla danego problemu — redukując macierz do 20–200 nuklidów.

Zastosowania sieci nuklidów w reaktorach jądrowych

Analiza sieci nuklidów w reaktorze obejmuje kilka problemów:

Wypalenie paliwa (fuel burnup). W miarę pracy reaktora: ²³⁵U ubywa (rozszczepienie), ²³⁸U pochłania neutrony → ²³⁹Pu (nowe paliwo), produkty rozszczepienia gromadzą się (trucizny: ¹³⁵Xe, ¹⁴⁹Sm). ChainFinder mapuje te ścieżki; ChainSolver oblicza skład paliwa po zadanym wypaleniu.

Aktywacja materiałów strukturalnych. Stal nierdzewna zawiera Fe, Ni, Cr, Co, Mn. Po napromieniowaniu neutronami: ⁵⁸Fe(n,γ)⁵⁹Fe, ⁵⁸Ni(n,p)⁵⁸Co, ⁵⁴Fe(n,p)⁵⁴Mn, ⁵⁹Co(n,γ)⁶⁰Co. ChainFinder szybko mapuje, które aktywacje są istotne dla danego materiału i strumienia neutronów.

Odpady radioaktywne. Co zostaje w wypalonym paliwie po 10, 100, 1 000, 10 000 lat? Bezpośrednio po wyłączeniu dominują krótkotrwałe produkty rozszczepienia. Po kilku latach — ⁹⁰Sr, ¹³⁷Cs. Po 100 latach — aktynoidy (Am-241, Np-237). ChainFinder pozwala prześledzić ścieżki transmutacji aktynoidów przez stulecia.

Ciepło powyłączeniowe (decay heat). Po wyłączeniu reaktora promieniowanie produktów rozszczepienia nadal wytwarza ciepło. To tzw. "decay heat" — ważne dla bezpieczeństwa (układ chłodzenia musi pracować nawet po wyłączeniu). ChainFinder mapuje, które izotopy dominują w decay heat w pierwszych sekundach, minutach, godzinach i dniach po wyłączeniu.

Zastosowania w medycynie nuklearnej

Produkcja izotopów medycznych jest kluczowym zastosowaniem analizy sieci nuklidów:

Mo-99 → Tc-99m. Molibden-99 jest produkowany przez: (a) wychwyt neutronów ⁹⁸Mo(n,γ)⁹⁹Mo lub (b) rozszczepienie ²³⁵U (fission molybdenum — wydajność ~6%). Technet-99m to izomer córki ⁹⁹Mo (T½ = 6 h) — generator Mo-99/Tc-99m jest wymieniany co ~1 tydzień w szpitalu. ChainFinder pozwala prześledzić ścieżkę produkcji.

Lu-177. Lutet-177 (T½ = 6,65 d) — ważny w terapii radioligandowej (np. ¹⁷⁷Lu-PSMA dla raka prostaty, ¹⁷⁷Lu-DOTATATE dla guza neuroendokrynnego). Produkowany przez: ¹⁷⁶Lu(n,γ)¹⁷⁷Lu (bezpośrednia aktywacja) lub ¹⁷⁶Yb(n,γ)¹⁷⁷Yb→¹⁷⁷Lu (pośrednia, przez iterb). Różne ścieżki dają różną czystość izotopową produktu — ważne dla dawkimetrii.

I-131. Jod-131 (T½ = 8,02 d) do terapii raka tarczycy i diagnostyki. Produkowany głównie przez rozszczepienie ²³⁵U (wysoka wydajność dla A=131). Ścieżka: ¹³¹Sn → ¹³¹Sb → ¹³¹Te → ¹³¹I. ChainFinder pokazuje ścieżkę β⁻ od cyny do jodu.

Ga-68 i PET. Gal-68 (T½ = 67,6 min) do PET. Produkowany w generatorze ⁶⁸Ge/⁶⁸Ga (T½ Ge-68 = 271 d). Ge-68 produkowany przez cyklotrony z ⁷⁰Ge(p,3n)⁶⁸Ga lub bezpośrednio z ⁶⁸Zn(p,n)⁶⁸Ga.

Sieci nuklidów w astrofizyce jądrowej

Fizyka nuklearna astrofizyczna to "grand scale" problemu sieci nuklidów:

Nukleosynteza gwiazdowa. Jak powstały ciężkie pierwiastki? W gwiazdach: reakcje jądrowe produkują pierwiastki od He do Fe przez "burning stages" (helium burning, CNO cykl, silicon burning). Sieć reakcji: setki nuklidów, tysiące przejść.

r-process. Rapid neutron capture process — szybkie pochłanianie neutronów przez jądra w środowiskach o ekstremalnym strumieniu neutronów (supernowe, zderzenia gwiazd neutronowych). Ścieżka r-process biegnie przez jądra daleko od doliny stabilności — właśnie te, które FAIR i FRIB mierzą po raz pierwszy.

s-process. Slow neutron capture process — powolne pochłanianie neutronów w gwiazdach AGB. Ścieżka biegnie przez dolina stabilności. Dane dla s-process są dobrze znane z pomiarów laboratoryjnych.

p-process. Proton-rich nuclides — niektóre izotopy nie są wytwarzane przez r- ani s-process. Mechanizm: reakcje fotojądrowe (γ,n) w gorących środowiskach.

Znaczenie dla NKE/ChainFinder. Kody astrofizyczne (np. NuGrid, CCSN codes) używają tych samych danych jądrowych (przekroje, T½) co NKE — tyle że dla tysięcy egzotycznych nuklidów i ekstremalnych temperatur. ChainFinder jest "edukacyjnym modelem" tych samych struktur grafowych.

Kody obliczeniowe: ORIGEN, FISPACT, SERPENT

Znane kody do obliczeń transmutacji i aktywacji:

ORIGEN (USA). ORIGEN (Oak Ridge Isotope GENeration) — najstarszy i najszerzej stosowany kod do obliczania aktywacji i transmutacji. ORIGEN2 (1980), ORIGEN-2.2 (2002), ORIGEN w ramach SCALE/TRITON. Używa bibliotek ORIGEN-ARP (Actinide and Reactor Products) opartych na ENDF/B. Standard w ocenach środowiskowych elektrowni jądrowych w USA.

FISPACT (Europa). FISPACT-II (UKAEA, Culham) — kod do obliczeń aktywacji materiałów, szczególnie dla fuzji (ITER). Używa bibliotek z JEFF i TENDL. Obsługuje: aktywację, wypalenie, decay heat, toxicity, clearance index.

SERPENT (Finlandia). SERPENT 2 — kod Monte Carlo do fizyki reaktorów z wbudowanym modułem wypalenia. Używa ENDF/B lub JEFF. Szczególnie popularny dla obliczeń paliwowych i projektowania rdzeni reaktorów.

OpenMC (USA). Open-source Monte Carlo (MIT). Moduł depletion (wypalenia) używa OpenMC-Depletion z ORIGEN lub własnym solverem macierzy eksponent.

ChainSolver vs ORIGEN. ChainSolver jest edukacyjnym narzędziem web, nie zastępuje ORIGEN dla obliczeń inżynierskich. Główna różnica: ORIGEN ma biblioteki dla całego reaktora (tysiące wstępnie obliczonych zestawów), ChainSolver oblicza "od zera" z danymi ORIP_XXI dla zadanego problemu.

Polska infrastruktura badań transmutacji

NCBJ i MARIA. Reaktor MARIA w Świerku umożliwia napromieniowanie próbek i produkcję izotopów. Obliczenia transmutacji materiałów (aktywacja elementów reaktora, skład paliwowy, odpady) wymagają narzędzi jak ChainFinder/ChainSolver lub pełnych kodów (ORIGEN, FISPACT).

IChTJ. Instytut Chemii i Techniki Jądrowej (Warszawa) prowadzi badania radiochemiczne, w tym produkcję izotopów medycznych i analizę skażeń środowiskowych. Obliczenia aktywacji i transmutacji są narzędziem codziennym.

Edukacja. Kalkulatory ChainFinder i ChainSolver są zaprojektowane jako narzędzia edukacyjne dla polskich doktorantów — umożliwiające zrozumienie sieciowych obliczeń transmutacji bez konieczności instalacji i obsługi pełnych kodów (ORIGEN, FISPACT).

Otwarte pytania badawcze

  1. Jak zmieniałaby się topologia sieci nuklidów (liczba węzłów, krawędzi, cykli), gdyby uwzględnić pełen zakres reakcji (n,2n), (γ,n), (p,n) zamiast tylko (n,γ)?

  2. Czy istnieje optymalna głębokość BFS dla typowych problemów produkcji izotopów medycznych — przy której sieć jest kompletna, ale obliczenia BFS nie są zbyt kosztowne?

  3. Jak porównuje się algorytm BFS ChainFindera z grafowymi bazami danych (np. Neo4j) dla zapytań o ścieżki w dużych sieciach nuklidów?

  4. Czy kombinacja ChainFinder (eksploracja topologii) + ChainSolver (obliczenia ilościowe) jest pedagogicznie skuteczniejsza od bezpośredniego korzystania z ORIGEN?

  5. Jak zmieniają się ścieżki dominujące w sieci nuklidów reaktora termicznego vs. szybkiego (SFR) — i jak to widać w grafach ChainFindera?

  6. W jaki sposób dane z FAIR i FRIB (nowe pomiary egzotycznych izotopów) zmienią grafy dla nuklidów w r-process i s-process astrofizyki jądrowej?

  7. Czy istnieje "kanoniczny" zestaw problemów transmutacji (benchmarki) dla testowania i weryfikacji kalkulatorów transmutacji?

  8. Jak obsłużyć w grafie sieciowym nuklidy z wieloma stanami izomerycznymi (g/m1/m2) — jakie wyzwania algorytmiczne to stwarza?

Słownik pojęć kluczowych

Graf skierowany (digraph) — struktura matematyczna złożona z węzłów i krawędzi ze skierowanymi strzałkami; sieć nuklidów to digraph, gdzie węzły = izotopy, krawędzie = przejścia.

BFS (Breadth-First Search) — algorytm przeszukiwania grafu "warstwa po warstwie" od węzła startowego; gwarantuje znalezienie najkrótszych ścieżek (w sensie liczby kroków).

Równania Batemana — układ liniowych ODE opisujący ewolucję stężeń nuklidów w łańcuchu (lub sieci) rozpadów i transmutacji; opublikowane przez H. Batemana w 1910.

Wypalenie paliwa (burnup) — miara zużycia paliwa jądrowego; zwykle wyrażana w MWd/tU (megawatodniach na tonę uranu); wzrost wypalenia = więcej produktów rozszczepienia, mniej ²³⁵U.

Decay heat — ciepło wytwarzane przez rozpad promieniotwórczy produktów rozszczepienia po wyłączeniu reaktora; kluczowy parametr bezpieczeństwa.

ORIGEN — Oak Ridge Isotope GENeration; najszerzej stosowany kod do obliczania aktywacji i transmutacji; standard w ocenach środowiskowych elektrowni jądrowych w USA.

Sztywność układu (stiffness) — właściwość układu ODE, gdzie najkrótszy i najdłuższy czas charakterystyczny różnią się o wiele rzędów wielkości; sieci nuklidów są klasycznie "sztywne" z T½ od nanosekund do miliardów lat.

Izotop medyczny — izotop produkowany specjalnie dla diagnostyki lub terapii medycznej; przykłady: Tc-99m (SPECT), F-18 (PET), Lu-177 (terapia radioligandowa).

r-process — rapid neutron capture process; sekwencja szybkich wychwytu neutronów w ekstremalnych środowiskach astrofizycznych (supernowe, neutron star mergers); produkuje ok. połowę ciężkich pierwiastków.

s-process — slow neutron capture process; sekwencja powolnych wychwytu w gwiazdach AGB; produkuje pierwiastki od Fe do Pb/Bi.

Podsumowanie dydaktyczne

  1. Sieć nuklidów jest naturalnym directed graphem — zrozumienie topologii tej sieci jest warunkiem wstępnym do poprawnej interpretacji obliczeń transmutacji.

  2. ChainFinder odpowiada na pytanie "jakie ścieżki istnieją" — ale nie mówi, które dominują; to wymaga ilościowych obliczeń w ChainSolverze.

  3. Równania Batemana z 1910 roku są fundamentem większości kodów transmutacji; nowoczesne metody (Padé, CRAM) to rozwinięcia tego samego formalnego aparatu.

  4. BFS jako algorytm jest optymalny dla znajdowania najkrótszych ścieżek w grafie nuklidów; głębokość eksploracji jest kluczowym parametrem kontrolującym kompletność vs. wydajność.

  5. Zastosowania praktyczne sieci nuklidów obejmują: wypalenie paliwa reaktorowego, produkcję izotopów medycznych, aktywację materiałów strukturalnych, decay heat, zarządzanie odpadami i astrofizykę jądrową.

  6. ORIGEN i FISPACT są "przemysłowymi" kodami transmutacji; ChainFinder/ChainSolver to edukacyjne narzędzia web pozwalające zrozumieć te same zasady fizyczne bez bariery instalacji.

  7. Polska ma realną infrastrukturę (reaktor MARIA, NCBJ, IChTJ) do obliczeń transmutacji — rozumienie narzędzi jak ChainFinder jest bezpośrednio aplikowalne zawodowo.

  8. Dane egzotycznych nuklidów (z FAIR, FRIB) będą sukcesywnie uzupełniać białe plamy w bazie ORIP_XXI — co rozszerzy zakres i dokładność grafów generowanych przez ChainFinder.

Historia obliczeń transmutacji — od papieru do kodu

Obliczenia transmutacji nuklidów mają długą historię, odzwierciedlającą postęp zarówno wiedzy fizycznej, jak i możliwości obliczeniowych:

Lata 1910.–1950.: obliczenia ręczne i proste łańcuchy. Harry Bateman (1910) opublikował analityczne rozwiązanie dla liniowego łańcucha A→B→C→D. Rozwiązanie jest eleganckie, ale dotyczy tylko prostych łańcuchów bez rozgałęzień i bez zewnętrznego napromienowania. Fizycy ery Projektu Manhattan (1942–1945) obliczali aktywację materiałów i produkcję plutonu metodami semi-analitycznymi — z dużym wysiłkiem ręcznym.

Lata 1950.–1960.: pierwsze programy komputerowe. Komputery IBM 704/709 umożliwiły implementację układów równań ODE. Bell Laboratory i Argonne National Laboratory pionierują kody do obliczenia aktywacji w reaktorach. Pierwsze wersje ORIGEN są rozwijane w Oak Ridge w połowie lat 1960. ORIGEN 1.0 (1973) — pierwszy szeroko stosowany kod.

Lata 1970.–1980.: standardo przemysłowe. ORIGEN 2 (1980, ORNL) staje się standardem dla obliczeń wypalenia i aktywacji w USA. W Europie rozwijany jest FISPACT (Culham). Kody te używają bibliotek stałych (pre-calculated data), nie rozwiązują układów ODE dla każdego problemu od nowa — co przyspiesza obliczenia.

Lata 1990.–2000.: Monte Carlo + depletion. SERPENT, MCNP (z modułem ORIGEN-SCALE) łączą transport Monte Carlo neutronów z obliczeniami wypalenia. To umożliwia przestrzennie niejednorodne obliczenia w geometriach reaktorowych. Dla małych elementarnych układów (jedna komórka paliwowa) obliczenia Monte Carlo + depletion są złotym standardem.

Lata 2010.–dzisiaj: open source i edukacja. OpenMC (MIT, open-source Monte Carlo) ma wbudowany solver depletion. CYCLUS (Wisconsin) — agent-based symulator cyklu paliwowego. ChainFinder/ChainSolver — web tools dla edukacji, oparte na nowoczesnych bibliotekach JavaScript i bazie ORIP_XXI.

Mechanizm wychwytów neutronów a produkcja aktynoidów

Ścieżka od ²³⁸U do transuraneów przez sekwencyjne pochłanianie neutronów jest jednym z kluczowych procesów w reaktorze:

Pu-239 z U-238.
²³⁸U + n → ²³⁹U (σ = 2,7 barns)
²³⁹U → ²³⁹Np + β⁻ (T½ = 23,5 min)
²³⁹Np → ²³⁹Pu + β⁻ (T½ = 2,36 d)

Pu-239 jest podstawowym "hodowanym" paliwem w reaktorach uranowych. Każdy reaktor LWR produkuje kilogramowe ilości Pu-239 rocznie.

Pu-240 i problem "trucizny".
²³⁹Pu + n → ²⁴⁰Pu (σ = 271 barns — wysoki!)
²⁴⁰Pu nie jest rozszczepialny przez neutrony termiczne. Wysoka zawartość ²⁴⁰Pu w plutonie reaktorowym jest problemem dla proliferatorów — ale nie dla reaktorów MOX, które radzą sobie z ²⁴⁰Pu.

Am-241 z Pu-241.
²⁴¹Pu → ²⁴¹Am + β⁻ (T½ = 14,3 lat)
Am-241 jest kluczowym długożyciowym aktynoidem w odpadach. W perspektywie 100–1000 lat dominuje aktywność alfa wielu próbek wypalonego paliwa.

Np-237 z U-235.
²³⁵U + 2n → ²³⁷U → ²³⁷Np (przez β⁻)
Lub: ²³⁶U(n,γ) → ²³⁷U → ²³⁷Np.
Np-237 (T½ = 2,14 × 10⁶ lat) jest dominującym "million-year"-scale aktynoidem w odpadach. Jego dane (przekroje, T½) mają kluczowe znaczenie dla bezpieczeństwa głębokich składowisk.

Cm-244 i dalej.
²⁴³Am(n,γ)²⁴⁴Am→²⁴⁴Cm. Cm-244 emituje duże ilości ciepła α. W wypalenym paliwie o wysokim stopniu wypalenia (>50 GWd/tU) izotopy Cm są istotnym źródłem ciepła wymagającego specjalnych środków przechowywania.

ChainFinder pokazuje te ścieżki graficznie — student widzi, że "pluton w reaktorze" to nie jeden izotop, lecz dynamicznie zmieniająca się mieszanina Pu-239, -240, -241, -242 wraz z Am i Cm.

Produkty rozszczepienia jako sieć nuklidów

Rozszczepienie ²³⁵U daje ~90 różnych bezpośrednich produktów ("primary fission fragments"), każdy bogaty w neutrony, który rozpadając się przez kilka β⁻ dochodzi do stabilnego jądra. To tworzy sieć 300–500 nuklidów:

Łańcuchy masowe. Każde A (np. A=95) ma swój "izobaryczny łańcuch" — szereg izotopów o tym samym A, rozpadających się przez β⁻ od bogatych neutronowo do stabilnych. Dla A=95: Rb-95 → Sr-95 → Y-95 → Zr-95 → Nb-95 → Mo-95 (stabilny). ChainFinder pokazuje ten łańcuch dla wybranego A.

Kryptonowe truby. Krypton-85 (T½ = 10,76 lat) jest emitowany do atmosfery w procesie PUREX (przeróbki chemicznej wypalonego paliwa). Mierzony globalnie stężenie Kr-85 w atmosferze rośnie — co koreluje z globalną skalą reprocessingu paliwa jądrowego.

Ksenon-135 i reaktywność. Xe-135 (T½ = 9,17 h, σ = 2,6 × 10⁶ barns!) jest wytwarzany zarówno bezpośrednio w rozszczepieniu, jak i przez β⁻ z Te-135 → I-135 → Xe-135. Po wyłączeniu reaktora I-135 nadal się rozpada na Xe-135 → wzrost ilości Xe-135 → pochłanianie neutronów ("xenon pit") → przez kilkanaście godzin niemożliwy rozruch reaktora.

Sam-149 i efekt samaru. Sm-149 (T½ ≈ stabilny, σ = 40 700 barns) jest wytwarzany przez: Nd-149 → Pm-149 → Sm-149 (przez β⁻, czas ~50 h). W stanie stacjonarnym jest utrzymywana równowaga między produkcją Sm-149 i jego pochłanianiem neutronów. Po wyłączeniu reaktora: Pm-149 nadal produkuje Sm-149 (nie jest pochłaniany → gromadzi się) → wzrost reaktywności ujemnej ("samarium poisoning") po ponownym rozruchu.

Grafy i wizualizacja — narzędzia informatyczne

ChainFinder renderuje grafy w SVG przez JavaScript. Warto wiedzieć, jakie algorytmy stoją za wizualizacją:

Układ grafu. Layout grafu nuklidów (pozycjonowanie węzłów) jest problemem NP-trudnym w ogólności. ChainFinder używa układu hierarchicznego (BFS-layered layout) — węzły z tej samej "głębokości" BFS są na tej samej pionowej warstwie. To jest czytelne dla łańcuchów liniowych, mniej dla gęstych sieci z wieloma cyklami.

Algorytmy layoutu. Ogólne narzędzia do layoutu grafów: Graphviz (DOT/neato/fdp), D3.js force-directed layout, Cytoscape.js. ChainFinder używa własnego prostego BFS-layered layoutu. Bardziej zaawansowane wizualizacje (np. Cytoscape, Gephi) są dostępne dla zaawansowanych użytkowników.

Kolorowanie wg trybu rozpadu. Kolorowanie węzłów jak w Karlsruher Nuklidkarte: niebieski = stabilny, żółty = α, zielony = β⁻, różowy = EC/β⁺, fioletowy = SF/IT. To umożliwia natychmiastową wizualną identyfikację dominującego trybu rozpadu każdego nuklidu w sieci.

Szerokość strzałek. Proporcjonalna do branching ratio lub wydajności rozszczepienia. Dla wychwytu neutronów: szerokość proporcjonalna do σ(n,γ). Dla rozszczepienia: szerokość proporcjonalna do cumulative yield. Wizualnie najgrubsze strzałki = najważniejsze ścieżki.

Krytyczne spojrzenie na grafy sieci nuklidów

Grafy mają ograniczenia, o których student musi wiedzieć:

Graf nie pokazuje wartości bezwzględnych. Gruba strzałka α z Pu-238 nie mówi, ile Pu-238 jest w próbce. Mówi tylko, że ta ścieżka ma wysoki branching ratio — ale jeśli Pu-238 jest w śladowych ilościach, jego wpływ na bilans aktywności może być pomijalny.

Graf nie pokazuje czasu. Krawędź "β⁻" łączy A z B przez T½ = 2 sekundy albo T½ = 1000 lat — obie wyglądają tak samo na grafie. Dopiero ChainSolver uwzględnia czas w obliczeniach ilościowych.

Graf nie uwzględnia gęstości strumienia. Moc transmutacji przez wychwyt neutronów zależy od φ·σ — ale graf pokazuje tylko, czy ścieżka istnieje, nie jak szybko przebiega przy danym φ.

Graf może być zbyt duży. Dla głębokości 8 i nuklidu startowego ²³⁸U: dziesiątki tysięcy potencjalnych węzłów. Parametry filtrowania (minimalne branching, minimalne σ) są niezbędne dla czytelności.

Szeregi naturalne — klasyczne przykłady łańcuchów

Naturalne szeregi promieniotwórcze są historycznie pierwszymi odkrytymi łańcuchami transmutacji. Są też dobrymi przykładami do ćwiczeń z ChainFinderem:

Szereg uranu (⁴ⁿ+² lub "szereg uranowy").
²³⁸U → ²³⁴Th → ²³⁴Pa → ²³⁴U → ²³⁰Th → ²²⁶Ra → ²²²Rn → ²¹⁸Po → ²¹⁴Pb → ²¹⁴Bi → ²¹⁴Po → ²¹⁰Pb → ²¹⁰Bi → ²¹⁰Po → ²⁰⁶Pb (stabilny).
14 przejść, 14 ogniw. T½ waha się od 1,64 × 10⁻⁴ s (Po-214) do 4,47 × 10⁹ lat (U-238). To klasyczna "kaskada rozpadu".

Szereg toru (⁴ⁿ lub "szereg torowy").
²³²Th → ²²⁸Ra → ²²⁸Ac → ²²⁸Th → ²²⁴Ra → ²²⁰Rn → ²¹⁶Po → ²¹²Pb → ²¹²Bi (tutaj rozgałęzienie: 64% α→²⁰⁸Tl, 36% β⁻→²¹²Po) → ²⁰⁸Tl → ²⁰⁸Pb lub ²¹²Po → ²⁰⁸Pb.
Ważna lekcja: Bi-212 pokazuje rozgałęzienie (forked decay).

Szereg aktyniowy (⁴ⁿ+³).
²³⁵U → ²³¹Th → ²³¹Pa → ²²⁷Ac → ²²³Ra → ... → ²⁰⁷Pb (stabilny). Najdłuższy i najmniej "gęsty" — mała naturalna zawartość U-235.

Neptunowy (⁴ⁿ+¹) — zanikły z przyrody.
²³⁷Np → ²³³Pa → ²³³U → ²²⁹Th → ²²⁵Ra → ... → ²⁰⁹Bi. Ten szereg nie występuje naturalnie (T½ Np-237 = 2,14 × 10⁶ lat — zbyt krótkie w skali Ziemi) — ale jest wytwarzany w reaktorach jądrowych!

Lekcja dla ChainFindera. Każdy z tych szeregów to klasyczny "graf liniowy" (lub z małymi rozgałęzieniami). Ćwiczenie: wygeneruj szereg ²³⁸U przez ChainFinder z głębokością 15 i tylko rozpadami (bez neutronów). Porównaj z tabelą powyżej.

Systemy monitoringu CTBT a sieci nuklidów

Kompleksowy Układ o Zakazie Prób Jądrowych (CTBT) korzysta z danych nuklidowych do identyfikacji śladów testów jądrowych:

Co rejestruje IMS. International Monitoring System (IMS) CTBTO rejestruje m.in. radionuklidy atmosferyczne: ¹³³Xe, ¹³³mXe, ¹³¹mXe, ¹³⁵Xe. Ksenon jest szczególnie ważny, bo jest gazem szlachetnym — nie jest pochłaniany przez gleby i budynki.

Izotopy ksenonowe jako wskaźniki. Ksenon-133 (T½ = 5,24 d), Xe-133m (T½ = 2,19 d), Xe-131m (T½ = 11,9 d), Xe-135 (T½ = 9,17 h). Stosunek aktywności tych izotopów pozwala wnioskować o:

  • czasie od wybuchu (im starszy, tym więcej krótkotrwałych izotopów zanikło)
  • typie urządzenia (jądrowe vs. fuzyjne vs. reaktorowy wyciek)

Sieci nuklidów CTBT. Interpretacja danych IMS wymaga rozwiązania odwrotnego problemu: dane = pomierzone aktywności izotopów → model = kiedy i co wybuchło? Narzędzia używają sieci rozpadów z danymi nuklidowymi (jak ORIP_XXI) do symulacji ewolucji sygnatur radionuklidowych w czasie.

ChainFinder w kontekście CTBT. Ćwiczenie: Xe-135 pochodzi z łańcucha β⁻ Te-135 → I-135 → Xe-135. Znajdź te ścieżki przez ChainFinder. Jak długo po wybuchu (przy danym T½) Xe-135 będzie dominował w sygnaturze radionuklidowej?

Problemy ochrony radiologicznej a sieci nuklidów

Ochrona radiologiczna wymaga rozumienia sieci nuklidów:

Dawka od mieszaniny produktów rozszczepienia. Bezpośrednio po wybuchu: tysięce nuklidów, każdy emituje promieniowanie. Po 1 godzinie: dominują krótkotrwałe izotopy. Po 1 tygodniu: zostają te o T½ ~ kilka dni. Po 1 roku: zostają ⁹⁰Sr i ¹³⁷Cs. ChainFinder pokazuje, które ścieżki prowadzą do "przeżywających" izotopów.

Reguła "7-10". Stara heurystyka: aktywność mieszaniny produktów rozszczepienia spada 10× co 7-krotny wzrost czasu. To przybliżenie działa przez kilka lat. Bardziej precyzyjne obliczenia wymagają właśnie narzędzi typu ChainSolver z pełną siecią nuklidów.

Inhalacja vs. zewnętrzna ekspozycja. Dla dawki inhalacyjnej: ważne są izotopy wewnątrzustrojowe (jod, cez, stront, pluton). Dla zewnętrznej: gamma emittery. ChainFinder + dane dawek z ICRP-119 pozwalają wybrać izotopy kluczowe dla ochrony radiologicznej w danym scenariuszu.

Weryfikacja danych — jak sprawdzić wyniki ChainFindera

Studenci powinni umieć zweryfikować wyniki ChainFindera niezależnie:

NNDC NuDat. National Nuclear Data Center (BNL) udostępnia NuDat 3.0 — interaktywną bazę danych nuklidów online. Dane podstawowe: T½, tryby rozpadu, energie gamma. Można porównać z ORIP_XXI.

Livechart IAEA. IAEA udostępnia "Livechart of nuclides" — interaktywną kartę online z danymi z NUBASE2020 i ENSDF. Łatwe do użycia w przeglądarce, bez instalacji.

ENSDF (Evaluated Nuclear Structure Data File). Najbardziej kompletna kompilacja schematów poziomów jąder. Dostępna przez NNDC. Każdy wpis zawiera dane eksperymentalne i ocenione wartości.

Weryfikacja ścieżek. Gdy ChainFinder wskaże ścieżkę A→B→C: zweryfikuj w NuDat lub Livechart, że:

  1. A rozpadasa do B (właściwy tryb, właściwy branching)
  2. B rozpadasa do C (właściwy tryb, właściwy branching)
  3. T½ są zgodne z oczekiwaniami scenariusza

Kontekst szerszego kursu — powiązania z innymi artykułami

ChainFinder jest narzędziem "poprzecznym" — ma zastosowania w wielu działach wiedzy reprezentowanych w tym portalu:

Fizyka reaktorów. Produkty rozszczepienia, aktynidy, wypalenie paliwa → powiązane z artykułami o UF₆, cyklu paliwowym, reaktorach.

Ochrona radiologiczna. Dawka od radionuklidów → powiązane z artykułami o pomiarach dawki, monitoringu środowiskowym.

Fizykachemia odpromieniotwórczości. Łańcuchy α i β w szeregach naturalnych → powiązane z artykułami o materiałach i izotopach (uran, tor, rad, radon).

Medycyna nuklearna. Ścieżki produkcji izotopów medycznych → powiązane z artykułami o Mo-99, Am-241, I-131.

Historia. Odkrycie naturalnych szeregów promieniotwórczych (Marie Curie, Ernest Rutherford) jako historia odkrycia transmutacji → powiązane z artykułami historycznymi.

Zrozumienie sieci nuklidów przez ChainFinder daje kontekst do wszystkich tych obszarów — jest "infrastrukturą pojęciową" dla szerokiego zakresu fizyki jądrowej stosowanej.

Szczegółowy przykład: produkcja I-131 przez rozszczepienie — śledzenie przez graf

Jod-131 (T½ = 8,02 d) jest jednym z kluczowych izotopów emitowanych po awariach reaktorów jądrowych. Jego ścieżka powstania z rozszczepienia ²³⁵U jest następująca:

Bezpośrednie produkty rozszczepienia dla A=131.
Rozszczepienie ²³⁵U daje fragmenty dla A=131 o różnych Z. Typowe bezpośrednie fragmenty: ⁴⁸Cd(n,f) — bardzo bogaty w neutrony, T½ ~microsek. Potem szybkie β⁻: Sn-131 → Sb-131 → Te-131 → I-131 → Xe-131 (stabilny).

Skumulowana wydajność rozszczepienia. Yield kumulacyjny dla I-131 z ²³⁵U: ~2,9% (termiczne neutrony). To oznacza: przy każdym 100 rozszczepieniach, ok. 2,9 atomów jodu-131 powstaje w sumie (bezpośrednio + przez łańcuch β⁻ poprzedników).

Ścieżka przez Te-131. Te-131 ma T½ = 25 min (stan podstawowy) i T½ = 30 h (izomer Te-131m). Izomer Te-131m → I-131 + β⁻ (97%) lub Te-131m → Te-131 + IT (3%). To rozgałęzienie jest widoczne w ChainFinder jako dwa węzły Te-131g i Te-131m.

Implikacje dla awarii (Fukushima, Czarnobyl). Bezpośrednio po awarii reaktora: I-131 narasta przez rozkład krótkotrwałych poprzedników (Sn, Sb, Te), potem spada z T½ = 8,02 d. Po ~3 miesiącach: <0,001% pozostaje. Dlatego tarczyca wchłaniająca jod ze skażonego powietrza lub mleka jest narażona głównie w pierwszych tygodniach.

Ćwiczenie z ChainFinderem. Startując od Te-131m, znajdź ścieżkę do I-131 i dalej do Xe-131 stabilnego. Ile kroków? Jakie T½ pośredników?

Obliczenia ChainFindera a programy zewnętrzne — porównanie

Studenci pytają: kiedy ChainFinder, kiedy ORIGEN?

ChainFinder. Używaj, gdy:

  • Chcesz zobaczyć, jakie ścieżki w ogóle istnieją
  • Szukasz konkretnych powiązań między nuklidami
  • Uczysz się konceptów sieci nuklidów
  • Planujesz eksperiment i chcesz sprawdzić, jakie produkty mogą powstać
  • Nie masz dostępu do ORIGEN lub FISPACT

ChainSolver. Używaj, gdy:

  • Znasz już sieć (z ChainFindera) i chcesz obliczyć ilości
  • Potrzebujesz ewolucji stężeń w czasie
  • Modelujesz aktywację materiału w reaktorze

ORIGEN/SCALE. Używaj, gdy:

  • Robisz obliczenia inżynierskie dla prawdziwego reaktora
  • Potrzebujesz certyfikowanych wyników (regulacje, bezpieczeństwo)
  • Masz dostęp do zwalidowanych bibliotek (ORIGEN-ARP, JEFF-CP)

FISPACT-II. Używaj, gdy:

  • Obliczasz aktywację materiałów dla fuzji (ITER, DEMO)
  • Potrzebujesz clearance index, radiotoksyczności
  • Masz dostęp do TEND/JEFF bibliotek

MCNP + depletion. Używaj, gdy:

  • Potrzebujesz przestrzennie niejednorodnych obliczeń
  • Modelujesz pełny rdzeń reaktora z geometrią 3D

Sieci nuklidów w bezpieczeństwie jądrowym (safeguards)

MAEA i krajowe organy regulacyjne używają sieci nuklidów w safeguards:

Analiza izotopowa paliwa. Po wyjęciu paliwa z reaktora: pomiar aktywności izotopów Cs, Nd, Eu (passywna neutronowa analiza izotopowa, PNAR). Stosunek ¹³⁷Cs/¹³³Cs, ¹⁴⁸Nd koreluje z wypaleniem. ChainFinder pozwala prześledzić, skąd te izotopy pochodzą i dlaczego korelują z wypaleniem.

Pu/U ratio. Ilość plutonu w paliwie vs. ilość uranu (depletion). Dane izotopowe Pu: ²³⁸Pu/²³⁹Pu/²⁴⁰Pu/²⁴¹Pu/²⁴²Pu. Określony "sygnałem reaktorowym" — zależy od spectrum neutronów, czasu napromieniania i historii cyklu. ChainFinder pomaga zrozumieć, dlaczego konkretne izotopy Pu mają takie proporcje.

Aktywność residualna HEU. W ramach safeguards dla materiałów jądrowych: analiza "signature" radionuklidów (m.in. ²³²U, ²³⁶U, ²³⁴U) w materiale HEU pozwala wnioskować o jego historii i ewentualnym reprocessingu. ChainFinder mapuje ścieżki tych rzadkich izotopów.

Sieci nuklidów w środowisku polonijnym — rola NCBJ Świerk

Polska infrastruktura naukowa ma bezpośredni związek z obliczeniami sieci nuklidów:

Reaktor MARIA. Reaktor badawczy MARIA (20 MW) w Otwocku-Świerku jest jednym z nielicznych reaktorów badawczych w Europie Środkowej wciąż w aktywnej eksploatacji. Napromieniowanie próbek w MARII produkuje izotopy przez: (a) wychwyt neutronów z naturalnych i wzbogaconych tarcz, (b) rozszczepienie tarcz aktynoidowych. Obliczenia aktywacji (jakie izotopy i w jakiej ilości powstaną) wymagają narzędzi opartych na sieciach nuklidów.

Produkcja izotopów medycznych. MARIA produkuje m.in.: ⁹⁹Mo → ¹³¹I → ¹⁵³Sm → ¹⁷⁷Lu (nowy, rosnące zapotrzebowanie). Każda produkcja wymaga obliczeń transmutacji dla optymalizacji czasu napromieniowania, składu tarczy i procedur radiochemicznych.

Edukacja na Świerku. NCBJ prowadzi kursy dla doktorantów z fizyki i techniki jądrowej. Kalkulatory ChainFinder i ChainSolver są narzędziami edukacyjnymi komplementarnymi do kursów prowadzonych na Świerku — pozwalają studentom eksplorować sieci nuklidów bez instalacji specjalistycznego oprogramowania.

Przykład zastosowania w Polsce. Wyobraź sobie następujący scenariusz dydaktyczny: student Wydziału Fizyki Politechniki Warszawskiej pracuje nad obliczeniami aktywacji stali nierdzewnej AISI 316L dla zastosowania w reaktorze. Korzysta z ChainFindera, aby zidentyfikować kluczowe izotopy aktywacyjne (Co-60, Fe-55, Ni-63, Mn-54, Cr-51), potem używa ChainSolvera do obliczeń ilościowych, a wyniki porównuje z danymi z ORIGEN-ARP uruchamianym na serwerze NCBJ. ChainFinder i ChainSolver są "edukacyjną bramą" do tej samej fizyki, którą oblicza ORIGEN na serwerach krajowych laboratoriów jądrowych.

Kalkulator ChainFinder reprezentuje zatem nie tylko narzędzie komputerowe, ale przede wszystkim podejście metodologiczne — od rozpoznania topologii sieci nuklidów, przez ilościowe obliczenia, do interpretacji wyników w kontekście konkretnych zastosowań inżynierskich i naukowych. Dla przyszłych polskich fizyków jądrowych, inżynierów bezpieczeństwa i specjalistów ochrony radiologicznej jest to umiejętność bazowa — fundament, na którym buduje się dalszą kompetencję w obliczeniach reaktorowych i zarządzaniu odpadami.

Dodatkowe materiały multimedialne

Narzędzia interaktywne

Ćwiczenia praktyczne

  1. Korzystając z ChainFindera, znajdź ścieżki prowadzące do produkcji Tc-99m w reaktorze neutronowym. Jakie nuklidy pośrednie muszą powstać?
  2. Porównaj grafy dla Co-60 w trybie "do przodu" i "wstecz". Które ścieżki prowadzą do Co-60 przez aktywację? Które przez rozpad?
  3. Dla rozszczepienia U-235, znajdź przez ChainFinder, jakie produkty zawierają łańcuchy β⁻ kończące się na stabilnym Cs-133. Ile kroków ma najkrótsza ścieżka?