Mechanizm rozszczepienia - model kroplowy

Streszczenie

Gdy patrzy się na rozszczepienie po raz pierwszy, łatwo odnieść wrażenie, że ciężkie jądro po prostu „pęka”, bo trafił w nie neutron. To zbyt uproszczony obraz. W użyteczniejszym przybliżeniu jądro zachowuje się jak naładowana kropla cieczy: siły jądrowe próbują utrzymać ją zwartą, a odpychanie elektryczne protonów próbuje ją rozepchnąć. Dopóki kropla pozostaje prawie kulista, układ jest jeszcze stabilny. Gdy jednak wzbudzenie i deformacja staną się wystarczająco duże, pojawia się przewężenie, po którym dwa końce jądra zaczynają odpychać się tak silnie, że rozszczepienie staje się nieuniknione.¹

Model kroplowy nie daje pełnego, mikroskopowego opisu struktury jądra, ale bardzo dobrze tłumaczy samą logikę procesu. Pokazuje, skąd bierze się bariera rozszczepienia, dlaczego uran-235 łatwiej rozszczepić neutronem powolnym niż uran-238, oraz dlaczego po scission zdecydowana większość energii trafia od razu do kinetyki fragmentów rozszczepienia.^1,2,3

Rozszerzenie tematu

Model kroplowy jądra wyrósł z potrzeby uchwycenia kolektywnego zachowania bardzo ciężkich jąder, zwłaszcza wtedy, gdy nie są już niemal idealnie kuliste. Historycznie kluczowe znaczenie miała analiza opublikowana przez Nielsa Bohra i Johna Wheelera w 1939 roku, tuż po odkryciu rozszczepienia.⁵ Istota modelu polega na potraktowaniu jądra jak kropli nieściśliwej cieczy o w przybliżeniu stałej gęstości. Taki obraz nie opisuje wszystkiego, ale bardzo dobrze nadaje się do zrozumienia deformacji i rozpadu dużych jąder.

W tym przybliżeniu działają dwie konkurencyjne tendencje. Pierwsza odpowiada efektowi podobnemu do napięcia powierzchniowego: krótkozasięgowe siły jądrowe faworyzują zwarte, możliwie mało zdeformowane kształty. Druga wynika z odpychania kulombowskiego protonów, które dąży do rozpychania ładunku na większą objętość i większą odległość między jego częściami.^1,2 Dla jąder ciężkich ta konkurencja jest szczególnie ostra, bo liczba protonów jest już bardzo duża, a więc składnik kulombowski ma znaczącą wagę.

To właśnie stąd bierze się intuicja, że ciężkie jądro jest tylko warunkowo stabilne. Kiedy pozostaje blisko kształtu kulistego, siły jądrowe jeszcze utrzymują je w jednym kawałku. Jeżeli jednak coś wprowadzi do układu dostateczne wzbudzenie, jądro może wejść w tryb oscylacji kolektywnych. W języku bardziej technicznym chodzi o przejście od słabego wzbudzenia do deformacji elipsoidalnych, a następnie do kształtu z wyraźnym przewężeniem.¹

W przypadku uranu-235 schemat jest zwykle opisywany jako:

$$n + \mathrm{U\text{-}235} \rightarrow \mathrm{U\text{-}236^*}$$

czyli wychwyt neutronu i utworzenie jądra wzbudzonego.^2,3 Gwiazdka nie jest tu detalem notacyjnym, lecz kluczową informacją: jądro po wychwycie nie przechodzi od razu do stabilnego stanu, tylko dysponuje energią wzbudzenia, która może zostać rozdysponowana między różne stopnie swobody. Jeżeli wystarczy jej do pokonania bariery rozszczepienia, rozpoczyna się droga ku scission.

Bariera rozszczepienia jest jednym z najważniejszych pojęć, jakie model kroplowy porządkuje. Można ją rozumieć jako energetyczny koszt przejścia od kształtu zwartego do krytycznie zdeformowanego. Jądro nie rozpada się natychmiast po każdej niewielkiej deformacji, bo początkowo wzrost powierzchni kosztuje energię i siły jądrowe nadal utrzymują całość. Dopiero po osiągnięciu odpowiedniego stopnia wydłużenia rosnący udział odpychania elektrostatycznego zaczyna przeważać. W języku geometrycznym mówi się wtedy o przejściu przez punkt siodłowy i wejściu w obszar, w którym dalsza ewolucja ku rozpadowi staje się korzystna energetycznie.¹

To pozwala zrozumieć różnicę między uranem-235 a uranem-238. Po wychwycie neutronu przez uran-235 energia uwolniona przy przegrupowaniu jądra wystarcza do pokonania bariery bez potrzeby dodatkowego wkładu energii ruchu neutronu. Dlatego izotop ten jest wolnorozszczepialny. Dla uranu-238 pozostaje natomiast deficyt rzędu około 1 MeV, więc potrzebny jest neutron szybki, który wniesie brakującą energię kinetyczną.² Właśnie ta różnica ma potem ogromne znaczenie dla reaktorów i dla konstrukcji broni.

Po przekroczeniu bariery jądro wydłuża się dalej, tworząc kształt przypominający hantlę albo dwie krople połączone szyjką. Etap końcowy nazywa się scission, czyli przecięciem szyjki i rozdzieleniem fragmentów. To nie jest spokojny podział, lecz gwałtowne przejście, po którym dwa silnie naładowane fragmenty zaczynają odpychać się elektrostatycznie z bardzo dużą energią.^1,3 Słowiński podaje, że głównym składnikiem energii rozszczepienia jest właśnie energia kinetyczna fragmentów, rzędu około 166 MeV z całkowitego bilansu około 200 MeV.³

Ten rozkład energii dobrze pasuje do intuicji modelu kroplowego. W chwili scission istotna część energii potencjalnej układu jest już zmagazynowana w konfiguracji dwóch dodatnio naładowanych fragmentów znajdujących się bardzo blisko siebie. Gdy przestają być związane wspólną szyjką, odpychanie kulombowskie gwałtownie zamienia tę energię na ruch postępowy obu części. To dlatego fala uderzeniowa i ogrzewanie otoczenia biorą się przede wszystkim z zatrzymywania ciężkich fragmentów w materiale, a nie z samego aktu „pęknięcia” rozumianego jako czysto geometryczne rozerwanie.

Model kroplowy tłumaczy też, dlaczego rozszczepienie nie daje zwykle dwóch identycznych połówek. Produkty mają rozkład statystyczny skupiony wokół dwóch zakresów liczb masowych, a nie jednej symetrycznej wartości.² Oznacza to, że deformująca się kropla może przejść do wielu możliwych kanałów scission. Sam model kroplowy dobrze tłumaczy istnienie takiej rodziny kanałów, choć dokładny wybór konkretnych fragmentów wymaga już bardziej szczegółowych opisów struktury jądra.

Trzeba też jasno zaznaczyć granice tego modelu. Nie jest to pełna teoria mikroskopowa, lecz przybliżenie kolektywne. Dobrze opisuje konkurencję między energią powierzchniową i kulombowską oraz sam kształt drogi do rozszczepienia, ale słabiej radzi sobie z detalami wynikającymi ze struktury powłokowej jąder. Mimo to pozostaje bardzo użyteczny dydaktycznie i historycznie, bo właśnie on umożliwił pierwsze sensowne rozumienie tego, dlaczego pewne jądra ciężkie są tak blisko granicy rozpadu.⁴

To ma też bezpośredni związek z wcześniejszymi artykułami o energii wiązania jądra na nukleon, przekroju czynnym na rozszczepienie i modelu powłokowym jądra atomowego. Wykres energii wiązania mówi, dlaczego rozszczepienie ciężkich jąder może być energetycznie opłacalne. Model kroplowy pokazuje, jak taki rozpad może w ogóle zajść od strony deformacji i bariery energetycznej, a model powłokowy dopełnia obraz tam, gdzie trzeba zrozumieć mikroskopową strukturę poziomów i szczególne efekty stabilizacji.

Historia: odkrycie rozszczepienia i analiza Bohra-Wheelera

Zrozumienie mechanizmu rozszczepienia narastało etapami, a każdy etap był odrębnym przełomem. W grudniu 1938 roku Otto Hahn i Fritz Strassmann w Berlinie przeprowadzili eksperymenty, w których bombardowali uran neutronami. Zamiast spodziewanych transuranowców otrzymali nieoczekiwanie barium — pierwiastek o połowę lżejszy od uranu.⁷ Wynik był tak zaskakujący, że Hahn uznał go za eksperymento-chemiczny i nie podjął się fizycznego wytłumaczenia, prosząc listownie Lise Meitner — swoją wieloletnią współpracowniczkę — o interpretację.

Meitner, która wcześniej musiała uciec z nazistowskich Niemiec do Szwecji, wspólnie ze swoim siostrzeńcem Ottoem Frischem rozważyła wyniki podczas świątecznego spaceru w lesie koło Kungälv. Wykorzystując właśnie model kroplowy Bohra — intuicję, że ciężkie jądro zachowuje się jak naładowana ciecz — doszła do wniosku, że jądro uranu mogło się rzeczywiście rozciągnąć i rozpaść na dwa przybliżenie równe kawałki. Obliczyli, że energia wyzwolona w takim procesie wyniosłaby około 200 MeV, co z grubsza zgadzało się z późniejszymi pomiarami.⁷ Meitner i Frisch opublikowali w lutym 1939 roku artykuł w Nature, wprowadzając termin „fission" — zapożyczony z biologii, gdzie oznaczał podział komórki.

Niels Bohr, który właśnie przybył do Stanów Zjednoczonych, otrzymał tę informację od Frischa jeszcze przed drukiem. Natychmiast dostrzegł jej głębię i razem z Johnem Wheelerem przystąpił do systematycznej analizy teoretycznej. Ich artykuł On the Mechanism of Nuclear Fission ukazał się we wrześniu 1939 roku w Physical Review — w tym samym numerze, co wiadomości o wybuchu II wojny światowej.⁵ Był to tekst fundamentalny: po raz pierwszy w pełni sformalizował barierę rozszczepienia, wyprowadził warunki fisylności i wyjaśnił, dlaczego uran-235 reaguje na neutrony termiczne, a uran-238 tylko na szybkie.

Warto zaznaczyć, że Hahn otrzymał w 1944 roku Nagrodę Nobla w dziedzinie chemii za odkrycie rozszczepienia jądrowego, natomiast Meitner, której wkład koncepcyjny był kluczowy, nie została uhonorowana — przypadek powszechnie uznawany dziś za jedną z największych niesprawiedliwości w historii naukowych nagród. Element o liczbie atomowej 109 — meitnerium — nosi jej imię.

{{image-full:https://upload.wikimedia.org/wikipedia/commons/thumb/9/9f/Lise_Meitner_1924.jpg/400px-Lise_Meitner_1924.jpg:Lise Meitner (1878–1968) — fizyczka austriacko-szwedzka, która wspólnie z Ottoem Frischem jako pierwsza wyjaśniła mechanizm rozszczepienia jądrowego i obliczyła wyzwalaną energię. Element 109, meitnerium, nosi jej imię.}}

Wzór Bethego-Weizsäckera i energie składowe

Model kroplowy staje się ilościowy dopiero przez wzór Bethego-Weizsäckera (SEMF — Semi-Empirical Mass Formula), zwany też formułą ciekłokroplową. Wzór daje energię wiązania jądra jako sumę pięciu składników:

$$B(A,Z) = a_V A - a_S A^{2/3} - a_C \frac{Z(Z-1)}{A^{1/3}} - a_A \frac{(A-2Z)^2}{A} + \delta(A,Z)$$

Każdy składnik ma interpretację fizyczną:

• składnik objętościowy $a_V A$: każdy nukleon oddziaływuje ze swoimi bezpośrednimi sąsiadami, więc wkład jest proporcjonalny do liczby nukleonów; $a_V \approx 15{,}8 \mathrm{MeV}$
• składnik powierzchniowy $-a_S A^{2/3}$: nukleony na powierzchni mają mniej sąsiadów, więc ich wkład jest mniejszy; odpowiada napięciu powierzchniowemu; $a_S \approx 18{,}3 \mathrm{MeV}$
• składnik kulombowski $-a_C Z(Z-1)/A^{1/3}$: odpychanie elektrostatyczne między protonami zmniejsza energię wiązania; $a_C \approx 0{,}714 \mathrm{MeV}$
• składnik asymetrii $-a_A (A-2Z)^2/A$: odchylenie od równej liczby protonów i neutronów jest niekorzystne energetycznie; $a_A \approx 23{,}2 \mathrm{MeV}$
• składnik parzystości $\delta$: jądra parzysto-parzyste są szczególnie stabilne, a nieparzysto-nieparzyste szczególnie słabe; $\delta = +12/\sqrt{A} \mathrm{MeV}$ dla pp, $0$ dla odd-$A$, $-12/\sqrt{A} \mathrm{MeV}$ dla nn⁸

Dla fizyki rozszczepienia kluczowy jest stosunek składnika kulombowskiego do powierzchniowego. Gdy jądro deformuje się elipsoidalnie z parametrem $\varepsilon$ (małe odchylenie od kuli), energia zmienia się jak:

$$E_\mathrm{def} = E_S^{(0)}\left(1 + \frac{2}{5}\varepsilon^2\right) + E_C^{(0)}\left(1 - \frac{1}{5}\varepsilon^2\right)$$

Bariera istnieje dopóki $\Delta E = \Delta E_S + \Delta E_C > 0$, czyli dopóki wzrost energii powierzchniowej przewyższa spadek energii kulombowskiej przy deformacji. Warunek zaniku bariery to:

$$x \equiv \frac{E_C^{(0)}}{2E_S^{(0)}} = \frac{a_C Z^2/A^{1/3}}{2a_S A^{2/3}} = \frac{a_C Z^2}{2a_S A} \geq 1$$

Parametr fisylności i granica stabilności

Bezwymiarowy parametr $x = Z^2/A \cdot a_C/(2a_S)$ nosi nazwę parametru fisylności. Jego wartość określa, jak blisko granicy rozszczepienia spontanicznego znajduje się dane jądro. Ponieważ $a_C/(2a_S) \approx 0{,}0194$, granica spontanicznego rozszczepienia odpowiada $Z^2/A \approx 51{,}5$.⁸

Wartości $Z^2/A$ dla ważnych nuklidów rozszczepialnych:

Nuklidium	$Z$	$A$	$Z^2/A$	Fisylność
U-235	92	235	36,1	wolnorozszczepialne
U-238	92	238	35,6	prożno-rozszczepialne
Pu-239	94	239	37,0	wolnorozszczepialne
Pu-240	94	240	36,8	spontanicznie (niskie prawdopodobieństwo)
Cf-252	98	252	38,1	spontanicznie (znaczące)
Fm-256	100	256	39,1	spontanicznie (szybkie)
Og-294	118	294	47,6	skrajnie niestabilne

Cf-252 rozszczepienia spontanicznie z półokresem ~85 lat (95% α, 3% spontaniczne rozszczepienie) i jest ważnym laboratoryjnym źródłem neutronów — emituje 2,314 neutronu na każde rozszczepienie. Używa się go do uruchamiania reaktorów, neutronowej analizy aktywacyjnej i terapii BNCT.⁸

Dla $Z^2/A > 51{,}5$ bariera rozszczepienia zanika, a jądro rozszczepia się natychmiast po powstaniu. Tak dzieje się z najciężej zaobserwowanymi nuklidami oganesson i flerov (flerovium Z=114). Czas życia tych jąder jest mierzony w mikrosekundach lub sekundach, a ich odkrycie wymagało eksperymentów na akceleratorach (JINR Dubna).

Bariera rozszczepienia: powierzchnia energii potencjalnej

Droga od kształtu kulistego do scission nie jest prosta — prowadzi przez wielowymiarową powierzchnię energii potencjalnej (PES — Potential Energy Surface). Model kroplowy w swojej klasycznej formie opisuje tę powierzchnię za pomocą dwóch parametrów:

• β₂ — deformacja kwadrupolowa (elongacja)
• β₃ — deformacja oktupolowa (asymetria mas)

Gdy uwzględni się tylko β₂, bariera jest symetryczna i znana jako punkt siodłowy (saddle point). W rzeczywistości jednak jądra ciężkie preferują kanały asymetryczne, co widać na rozkładzie produktów rozszczepienia skupionym wokół $A \approx 95$ i $A \approx 140$, a nie wokół $A \approx 118$. Ta asymetria musi być uwzględniona przez model.

Nowoczesne obliczenia PES łączą opis makroskopowy (kroplowy) z mikroskopowym (efekty powłokowe) w formalizmu makro-mikro (Strutinsky prescription):

$$E_\mathrm{def} = E_\mathrm{LDM}(\beta) + \delta E_\mathrm{shell}(\beta)$$

Składnik $\delta E_\mathrm{shell}$ jest odpowiedzialny za to, że pewne konfiguracje deformacji są szczególnie preferowane lub szczególnie niestabilne — dokładnie tak jak w modelu powłokowym dla kulistych jąder, ale teraz dla kształtów elipsoidalnych. To właśnie ten efekt powoduje, że rozszczepienie U-235 prowadzi asymetrycznie do Ba-144 + Kr-89, a nie do dwóch jąder Sn-118.⁹

Scission: chwila podziału i jej dynamika

Punkt scission to moment, w którym szyjka łącząca dwa przedembrionalne fragmenty staje się niestabilna i pęka. Fizycznie nie jest to chwila momentalna, lecz krótki etap dynamiczny trwający rzędu $10^{-22} \mathrm{s}$. W tym czasie:

1. Szyjka przerywa się — dwa naładowane fragmenty nagle przestają być związane
2. Energia kulombowska konfiguracji dwóch bliskich ładunków natychmiast zaczyna napędzać ich wzajemne odpychanie
3. Emitowane są neutrony prędkie (prompt fission neutrons) w liczbie $\bar{\nu} \approx 2{,}4$ dla U-235 i $\approx 2{,}9$ dla Pu-239
4. Emitowane są gammy prędkie (prompt gamma rays)

Bilans energetyczny pełnego zdarzenia rozszczepienia U-235 neutronem termicznym:

Składnik energii	Wartość [MeV]
Energia kinetyczna fragmentów (TKE)	~166
Energie neutronów prędkich	~5
Energie gamm prędkich	~7
Energia β⁻ fragmentów	~7
Energie antyneutrin	~9
Energie gamm opóźnionych	~6
Razem	~200

Energia antyneutrin jest praktycznie niedostępna — neutrino przelatuje przez reaktor bez oddziaływania. W praktyce reaktorowej wykorzystywana jest suma z wyjątkiem antyneutrin: ~191 MeV per rozszczepienie. Przekłada się to na $3{,}2 \times 10^{10}$ rozszczepień na sekundę na 1 wat mocy termicznej.⁸

Asymetria rozkładu produktów rozszczepienia

Jednym z najbardziej charakterystycznych wyników rozszczepienia ciężkich jąder jest to, że produkty NIE rozkładają się symetrycznie. Rozkład masy produktów (fission yield) dla U-235 neutronami termicznymi ma kształt podwójnego garbu z maksimami przy:

• lekkie fragmenty: $A \approx 90$–$100$ (okolice Sr, Zr, Mo, Rb)
• ciężkie fragmenty: $A \approx 130$–$145$ (okolice Te, I, Cs, Ba, La, Ce)

Minimum przy $A \approx 118$ (symetryczny podział) ma wydajność $\sim 0{,}01\%$, podczas gdy maksima sięgają $\sim 6\%$. Sumaryczna wydajność obu ramion wynosi dokładnie 200% (każde rozszczepienie daje 2 fragmenty).

Fizyczną przyczyną asymetrii są efekty powłokowe. Cięższy fragment z okolicą $A \approx 140$ korzysta z magicznej lub pobliskiej magicznej liczby neutronów $N = 82$. Lekki fragment z okolicą $A \approx 95$ ma natomiast $Z \approx 38$, co jest związane z pobliską magiczną liczbą protonów $Z = 40$. Ten efekt stabilizacji powłokowej sprawia, że rozszczepienie asymetryczne jest statystycznie znacznie bardziej prawdopodobne niż symetryczne.⁹

Dla Pu-239 rozkład jest nieco różny niż dla U-235, z lekkim przesunięciem ciężkiego ramienia ku wyższym masom. Dla Cf-252 (spontaniczne rozszczepienie) rozkład jest jeszcze bardziej asymetryczny. Ta zależność od nuklidu macierzystego jest aktywnym obszarem badań — modele GEF (General description of Fission observables), CGMF i FREYA próbują odtworzyć pełny rozkład produktów dla reaktorowych obliczeń ekranowania i prognozowania produkcji izotopów medycznych.⁹

Rozszczepienie fotojądrowe i spontaniczne

Oprócz rozszczepienia wywołanego neutronem istnieją inne mechanizmy inicjowania tego procesu:

Rozszczepienie fotojądrowe (photofission): foton gamma o energii powyżej bariery rozszczepienia (typowo $>5 \mathrm{MeV}$ dla aktynowców) może być wchłonięty przez jądro i wywołać rozszczepienie. Mechanizm przebiega przez wzbudzenie Giant Dipole Resonance (GDR), omówione w artykule o wzbudzeniach jąder atomowych. Przekrój czynny dla fotofisji U-238 wzrasta gwałtownie od progu i osiąga maksimum ~200 mb przy energiach 10–20 MeV. Technika ta jest stosowana w detekcji materiałów jądrowych ukrytych w kontenerach — źródła fotonów z linijnymi akceleratorami elektronów (LINAC) mogą wzbudzać rozszczepienie nawet w ołowianych ekranach.⁹

Rozszczepienie spontaniczne: model kroplowy przewiduje, że bariera rozszczepienia nie jest nieskończona, a dla bardzo ciężkich jąder może być na tyle niska, że rozszczepienie może zajść na drodze kwantowego tunelowania. Prawdopodobieństwo tunelowania jest wykładniczo zależne od szerokości i wysokości bariery. Dla uranu spontaniczne rozszczepienie jest marginalnym kanałem rozpadu ($\lambda_{sf} \approx 10^{-17} \mathrm{s}^{-1}$ dla U-235), ale dla plutonu-240 jest już istotniejsze ($t_{1/2,sf} \approx 1{,}15 \times 10^{11} \mathrm{lat}$). Dla Cf-252 oba kanały — alfa i spontaniczne — są porównywalne.

Problem z Pu-240 ma szczególne znaczenie proliferacyjne: reaktorowy pluton zawiera do ~25% Pu-240 (Reactor Grade Plutonium, RGPu), który emituje intensywnie neutrony ze spontanicznego rozszczepienia — ok. 22 000 n/s/kg. To powoduje ryzyko predetonacji w projekcie gun-type i jest jednym z powodów, dla których do broni używa się wyłącznie Weapons Grade Plutonium (WGPu, <7% Pu-240, <1000 n/s/kg).¹⁰

Model Nilssona i powłokowo-kroplowy opis deformacji

Klasyczny model kroplowy Bohra-Wheelera był olbrzymim krokiem naprzód, ale nie mógł wyjaśnić szczegółów rozkładu masy produktów, przekrojów czynnych dla specyficznych kanałów ani dokładnych wartości barier. Kluczowym ulepszeniem było połączenie modelu kroplowego z modelem powłokowym — podejście makro-mikro.

Sven Gösta Nilsson w 1955 roku zaproponował model powłokowy dla zdeformowanych jąder, w którym poziomy energetyczne zależą od kształtu potencjału (parametru deformacji β). Diagram Nilssona pokazuje, jak poziomy jednocząstkowe przesuwają się wraz z deformacją. Przy pewnych wartościach β pojawiają się luki w spektrum poziomów — analogon magicznych liczb, ale teraz dla jąder niesferycznych. To są tzw. deformed magic numbers (np. $N = 152$ dla silnie zdeformowanych aktynowców).

Poprawka Strutinskiego (1967) sformalizowała połączenie obu podejść. Energia całkowita jądra przy danej deformacji β jest sumą gładkiej części (model kroplowy) i oscylacyjnej poprawki powłokowej:

$$E(\beta) = E_\mathrm{LDM}(\beta) + \delta U(\beta)$$

gdzie $\delta U$ jest poprawką powłokową — może być dodatnia lub ujemna, zmieniając kształt PES.⁹

To podejście wyjaśnia istnienie izomerów rozszczepialnych: stanów izomerycznych (np. Am-242m) o bardzo dużej deformacji (podwójne zagłębienie na PES) z barierą wewnętrzną, przez którą jądro musi tunelować, by spontanicznie rozszczepiać się lub wrócić do stanu podstawowego. Izomer U-236 ma drugą minimalną deformację $\beta \approx 0{,}6$, co daje czas życia $\sim 100 \mathrm{ns}$ jako izomer, a następnie spontanicznie rozszczepia się.⁹

Neutronowa statystyka rozszczepienia i multipliczność

Każde rozszczepienie emituje kilka neutronów prędkich (prompt neutrons). Ich liczba jest nieregularna i opisana rozkładem prawdopodobieństwa $P(\nu)$, gdzie $\nu = 0, 1, 2, 3, \ldots$ Parametrami tego rozkładu są:

• $\bar{\nu}$ — średnia liczba neutronów na rozszczepienie
• $\sigma^2(\nu)$ — wariancja
• rozkład ujemno-dwumianowy lub rozkład Gadgila

Wartości $\bar{\nu}$:

Nuklidium + typ neutronów	$\bar{\nu}$
U-235, neutron termiczny	2,436
U-235, neutron 1 MeV	2,520
Pu-239, neutron termiczny	2,879
Pu-239, neutron 1 MeV	2,960
U-238, neutron 2 MeV (próg ~1 MeV)	2,510
Cf-252, spontaniczne	3,757

Energia widmowa neutronów prędkich opisana jest widmem Watta:

$$\chi(E) = C \cdot e^{-E/a} \sinh\!\sqrt{bE}$$

gdzie dla U-235 parametry Watta to $a \approx 0{,}988 \mathrm{MeV}$, $b \approx 2{,}249 \mathrm{MeV}^{-1}$, a energia średnia neutronu prędkiego wynosi $\bar{E}_n \approx 2 \mathrm{MeV}$. To widmo jest bezpośrednio używane w obliczeniach reaktorów i tarcz ochronnych.⁸

Produkty rozszczepienia a bezpieczeństwo reaktorów

Model kroplowy wyjaśnia logikę podziału, ale to konkretne produkty rozszczepienia tworzą praktyczne wyzwania dla bezpieczeństwa i eksploatacji reaktorów:

Ksenon-135 ($t_{1/2} = 9{,}17 \mathrm{h}$) jest córką jodu-135 ($t_{1/2} = 6{,}57 \mathrm{h}$), który z kolei pochodzi z cezju-135 i baru-135. Ksenon-135 ma przekrój czynny na pochłanianie neutronów $\sigma = 2{,}65 \times 10^6 \mathrm{b}$ — największy spośród wszystkich nuklidów naturalnych i reaktorowych. Przy zmniejszeniu mocy reaktora lub po wyłączeniu, stężenie Xe-135 gwałtownie rośnie (ksenon pit), truciznując reaktor i uniemożliwiając jego restart przez 6–8 godzin. W reaktorze RBMK w Czarnobylu próba ominięcia tego efektu była jednym z czynników katastrofy.¹⁰

Stront-90 ($t_{1/2} = 28{,}8 \mathrm{lat}$) i cez-137 ($t_{1/2} = 30{,}2 \mathrm{lat}$) to główne długotrwałe składniki zanieczyszczeń po awariach reaktorów. Sr-90 jest analogiem wapnia i gromadzi się w kościach; Cs-137 jest analogiem potasu i jest szybko wchłaniany przez organizmy. Oba były uwolnione w znacznych ilościach w Czarnobylu i Fukushimie.¹⁰

Jod-131 ($t_{1/2} = 8{,}02 \mathrm{d}$) jest krótkotrwały, ale biologicznie aktywny w tarczycy — główne zagrożenie w pierwszych tygodniach po awarii. Tabletki jodku potasu blokują wchłanianie J-131 przez zasycenie receptorów tarczycy jodem stabilnym.

Zrozumienie tego, które nuklidy powstają i w jakiej ilości, wymaga znajomości pełnego rozkładu wydajności rozszczepienia — informacji, która jest bezpośrednio zakorzeniona w mechanizmie kroplowym i efektach powłokowych opisujących wybór kanałów scission.

Kody obliczeniowe i modelowanie rozszczepienia

Dokładne przewidywanie właściwości rozszczepienia wymaga zaawansowanych narzędzi numerycznych, bo analityczny model kroplowy jest zbyt uproszczony dla potrzeb inżynieryjnych. Współcześnie używa się kilku podejść:

GEF (General description of Fission observables, Karl-Heinz Schmidt i in., GSI Darmstadt) — kod oparty na modelu makro-mikro z empirycznymi parametrami, który odtwarza rozkłady masy, ładunku i energii kinetycznej produktów dla setek nuklidów. Używany m.in. przez ENDF/B-VIII.0 i JEFF-3.3 jako silnik generowania tabel wydajności.⁹

CGMF (Correlated Gamma and Neutron from Fission, LANL) — kod Monte Carlo symulujący pełny łańcuch emisji neutronów i gamm prędkich po scission, włącznie z korelacjami kąt-energia. Istotny dla programów detekcyjnych (DNDO).

FREYA (Fission Reaction Event Yield Algorithm, LLNL) — podobny do CGMF, z naciskiem na korelacje wielocząsteczkowe i ich zastosowania w spektroskopii.

MCNP i OpenMC — kody transportu neutronów Monte Carlo używające bibliotek wydajności rozszczepienia (ENDF/B) do symulacji pełnych środowisk reaktorowych. Poprawna biblioteka wydajności jest kluczowa dla obliczenia składu paliwa wypalonego i aktywności odpadów.

Współczesny standard to sprzężenie statystycznych PES z kodami transportu neutronów, co pozwala na pełną symulację od wyprodukowania jądra wzbudzonego przez neutron do finalnych produktów β-rozpadu, z uwzględnieniem termalizacji neutronów opóźnionych.⁸

Rozszczepienie w gwiazdach i synteza pierwiastków ciężkich

Model kroplowy ma zastosowanie nie tylko w reaktorach i broni — opisuje też procesy zachodzące w środowiskach astrofizycznych o ekstremalnej gęstości neutronów. W procesie r (rapid neutron capture process), który przebiega przy strumieniach neutronów $> 10^{20} \mathrm{n cm^{-2} s^{-1}}$, jądra mogą absorbować wiele neutronów zanim rozpadną się przez emisję beta. Jeżeli jądro staje się wystarczająco ciężkie, może ulec rozszczepieniu, a produkty stają się ziarnem dla dalszego wychwytywania — to tzw. pętla rozszczepienia w procesie r (r-process fission cycling).

Obserwacje fal grawitacyjnych i promieniowania elektromagnetycznego ze zderzenia gwiazd neutronowych GW170817 (sierpień 2017) potwierdziły po raz pierwszy bezpośrednio, że takie zderzenia produkują ogromne ilości pierwiastków transuranowych, które następnie rozszczepiają się i odpadają, zasilając kilonową emisję w zakresie podczerwieni. To odkrycie w bezprecedensowy sposób połączyło fizykę jądrową z astrofizyką obserwacyjną — bariera rozszczepienia wyznacza, które jądra mogą przeżyć do obserwowanej nukleosyntezy.⁸

Szczególnie interesujące jest to, że r-process zatrzymuje się przy $N = 126$ (magiczna liczba neutronów) — tworząc charakterystyczny pik w rozkładzie obfitości pierwiastków z $A \approx 195$. To bezpośredni odpowiednik piku powłokowego w rozkładzie produktów rozszczepienia, ale tym razem obserwowany w składzie chemicznym Galaktyki.

Rozszczepienie a proliferacja jądrowa — perspektywa MAEA

Model kroplowy jest nie tylko narzędziem dydaktycznym — ma bezpośrednie implikacje dla nieproliferacji. Parametr $Z^2/A$ mówi, które nuklidy potencjalnie mogą być użyte jako materiał rozszczepialny w broni. MAEA definiuje Significant Quantity (SQ) jako minimalną masę materiału, z której — przy zastosowaniu najlepszych znanych metod — można wyprodukować broń jądrową:

• HEU (wzbogacony >90% w U-235): SQ = 25 kg
• Pu-239: SQ = 8 kg
• U-233: SQ = 8 kg
• LEU (wzbogacony 20–90%): SQ = 75 kg (w zależności od wzbogacenia)

Wartości te są pochodną masy krytycznej i wymagają znajomości bariery rozszczepienia oraz parametrów fisylności.¹⁰ MAEA monitoruje przepływy tych materiałów w ramach systemu gwarancji (safeguards), a każde państwo-sygnatariusz NPT musi zgłaszać ich ilości. Umowy INFCIRC/153 i Additional Protocol definiują inspekcje i weryfikację deklaracji o materiałach jądrowych.

Raporty kwartalnego rozliczania (Material Balance Report, MBR) wymagają od operatorów reaktorów deklarowania każdego kilograma uranu i plutonu. Model kroplowy — przez wyznaczanie masy krytycznej — określa, ile materiału musi zostać rozliczone, by weryfikacja była sensowna technicznie.

Eksperymenty historyczne i bazy danych

Historycznym testem modelu kroplowego w praktyce były eksperymenty na krytycznych zestawach badawczych. Godiva I (1952), sferyczny zestaw HEU w Los Alamos bez reflektora, był pierwszym urządzeniem, które potwierdziło obliczenia masy krytycznej dla gołego HEU. Seria Jezebel (Pu-239), a potem Godiva II i III, pozwoliły systematycznie zmierzyć parametry kinetyki jądrowej dla różnych geometrii i konfiguracji.⁵

Szczególne znaczenie miały eksperymenty z pomiarem bariery rozszczepienia przez spektroskopię neutronów prędkich emitowanych przy rozszczepieniach. Porównanie zmierzonego widma neutronów z przewidywaniami teoretycznymi potwierdziło lub weryfikowało modele PES. Te eksperymenty były zatem jednocześnie fizyką jądrową i testem modelu kroplowego jako teorii.

Baza ICSBEP (International Criticality Safety Benchmark Evaluation Project, NEA OECD) zawiera ponad 5000 zestawów benchmarkowych od roku 1945 do współczesności. Każdy zestaw to udokumentowany eksperyment krytyczny, który można użyć do walidacji kodu obliczeniowego. Pełny katalog jest używany zarówno przez elektrownie jądrowe, jak i inspekcje MAEA.⁵

Równolegle NNDC (National Nuclear Data Center, BNL) prowadzi bazy danych wydajności rozszczepienia (ENDF/B-VIII.0, sekcja MF=8), dostępne online. Dla każdego nuklidu macierzystego zawierają tabele kumulatywnej i niezależnej wydajności poszczególnych produktów rozszczepienia, indeksowane po $A$, $Z$, izomerze i energii neutronu inicjującego. Dane te są niezbędne zarówno do analiz reaktorowych, jak i do interpretacji pomiarów MAEA.⁸

Znaczenie dydaktyczne i historyczne modelu

Model kroplowy Bohra-Wheelera pozostaje czymś wyjątkowym w fizyce jądrowej: jest to teoria, która przy stosunkowo prostym aparacie matematycznym daje prawdziwe wgląd w fizykę. Nie wymaga pełnego formalizmu kwantowego do wyjaśnienia, dlaczego pewne jądra są wolnorozszczepialne, a inne nie, dlaczego bariera istnieje i co decyduje o energii uwolnionej w procesie.

Dydaktycznie jest to jeden z najlepszych przykładów, jak analogia fizyczna (kropla cieczy) może prowadzić do ilościowej teorii niepozornie skomplikowanego zjawiska. Studenci fizykij jądrowej na całym świecie — w tym w Polsce, gdzie podręcznikami są m.in. Słowiński, Mayer-Kuckuk i Segre — zaczynają właśnie od modelu kroplowego, by mieć punkt oparcia dla wszystkich późniejszych, bardziej zaawansowanych opisów.

Historycznie model ten był też absolutnie kluczowy dla Projektu Manhattan. Leo Szilard opatentował pojęcie reakcji łańcuchowej już w 1934 roku, ale to analiza Bohra-Wheelera z 1939 roku udowodniła matematycznie, że rozszczepienie jest fizycznie możliwe i energetycznie opłacalne dla uranu-235 z neutronami termicznymi. To skierowało wysiłki Projektu Manhattan na wzbogacanie uranu i produkcję plutonu — bez tego modelu strategia wojskowa mogłaby pójść w bardzo różnych kierunkach.⁷

Dalszy kontekst: od modelu kroplowego do teorii Bosego-Einsteina i fazy jądrowej

Model kroplowy nasuwa pytanie: dlaczego materia jądrowa zachowuje się jak ciecz, a nie jak gaz lub ciało stałe? Odpowiedź leży w charakterze sił jądrowych. Krótkozasięgowe oddziaływanie silne między nukleonami ma charakter nasyconego — każdy nukleon odczuwa tylko sąsiadów. W połączeniu z zasadą Pauliego (nukleony są fermionami i nie mogą zajmować tych samych stanów kwantowych) prowadzi to do gęstości nasycenia $\rho_0 \approx 0{,}16 \mathrm{fm}^{-3}$, niezależnej od liczby masowej. To jest właśnie cecha cieczy — stała gęstość niezależnie od objętości.

Równanie stanu materii jądrowej, wyznaczone przez kombinację modelu kroplowego, modelu Nilssona i obliczeń Brücknera, definiuje tzw. symetrię energii i ściśliwość K — kluczowe dla zrozumienia zarówno jąder atomowych, jak i gwiazd neutronowych. Ściśliwość K = 200–260 MeV (wyznaczana z rezonansów monopolowych, tzw. oddychania jądra) wskazuje na to, jak sztywna jest materia jądrowa pod ciśnieniem — co bezpośrednio wpływa na równania stanu gwiazd neutronowych i wymaganą masę graniczną (masa Chandrasekara-Oppenheimera-Volkoffa).

Zaskakujące jest to, że ta sama teoria, która opisuje deformację i rozszczepienie jądra zawierającego 235 nukleonów, jest podstawą dla kosmologicznych modeli gwiazd neutronowych zawierających $\sim 10^{57}$ nukleonów. Zarówno bariera rozszczepienia jądra uranu, jak i bariera przeciwko kolapsu do czarnej dziury, wynikają z tych samych fundamentalnych właściwości materii jądrowej. Model kroplowy, choć zaproponowany w 1939 roku, pozostaje zatem żywym narzędziem — nie tylko dydaktycznym reliktem, ale aktywnym punktem odniesienia dla nowoczesnych badań jądrowych i astrofizycznych, z zastosowaniami sięgającymi od bezpieczeństwa reaktorów cywilnych aż po interpretację sygnałów grawitacyjnych ze zderzeń gwiazd neutronowych.

{{image-full:https://upload.wikimedia.org/wikipedia/commons/thumb/3/3d/Fission_chain_reaction.svg/600px-Fission_chain_reaction.svg.png:Schemat reakcji łańcuchowej rozszczepienia: neutron wchłonięty przez jądro uranu prowadzi do rozszczepienia, uwolnienia energii i kolejnych neutronów, które wchłaniane są przez sąsiednie jądra. Model kroplowy wyjaśnia każdy pojedynczy akt rozszczepienia w tym łańcuchu.}}

Warto dodać jeszcze jeden most pojęciowy. W materiałach Królikowskiego deformacje nie pojawiają się wyłącznie jako etap drogi do scission, ale także jako źródło bardziej regularnych wzbudzeń kolektywnych: pasm rotacyjnych i oscylacyjnych jąder zdeformowanych.⁶ To ważne rozróżnienie. Ten artykuł opisuje skrajną deformację prowadzącą do podziału jądra, natomiast osobny tekst o wzbudzeniach jąder atomowych, przejściach gamma i pasmach rotacyjnych opisuje łagodniejsze deformacje, w których jądro nie rozpada się, tylko przechodzi między kolejnymi stanami kolektywnymi.

Najkrótsze podsumowanie jest więc takie: model kroplowy nie zastępuje całej fizyki jądrowej, ale bardzo dobrze porządkuje najważniejszy ciąg przyczynowy. Neutron wzbudza ciężkie jądro, wzbudzone jądro deformuje się, deformacja może pokonać barierę rozszczepienia, a po scission energia kulombowska dwóch fragmentów zamienia się głównie na ich ruch. To właśnie ten łańcuch zdarzeń stoi za podstawowym mechanizmem rozszczepienia ciężkich aktynowców.^1,2,3

Dodatkowe materiały multimedialne

Materiał wizualny pokazujący intuicyjny obraz modelu kroplowego:

Jako uzupełnienie warto przejść dalej do energii wiązania jądra na nukleon, defektu masy, modelu powłokowego jądra atomowego i tekstu o reakcji łańcuchowej. Model kroplowy tłumaczy, dlaczego jądro chce się podzielić; dopiero te artykuły pokazują skąd bierze się energia, jak wygląda mikroskopowa struktura poziomów i kiedy pojedynczy akt rozszczepienia zamienia się w proces lawinowy.

Powiązane kalkulatory i narzędzia

• Energia wiązania — porównuje masy AME2020, defekt masy i przybliżenie Bethego-Weizsäckera.

Ćwiczenia praktyczne

Ćwiczenie laboratoryjne powinno polegać na zbudowaniu prostego modelu energetycznego deformacji ciężkiego jądra w duchu modelu kroplowego. W wariancie podstawowym należy:

1. przyjąć uproszczony opis energii jako sumy składnika powierzchniowego i kulombowskiego,
2. wprowadzić parametr deformacji opisujący odejście od kształtu kulistego,
3. zbadać jakościowo, jak zmienia się energia całkowita przy rosnącej deformacji,
4. zaznaczyć obszar odpowiadający barierze rozszczepienia i obszar po jej przekroczeniu,
5. porównać zachowanie modelu dla jądra „łatwiej” i „trudniej” rozszczepialnego.

Celem ćwiczenia nie jest uzyskanie idealnie dokładnych wartości liczbowych, lecz odtworzenie samej logiki bariery: początkowy wzrost energii przy deformacji, a następnie przejście do stanu, w którym odpychanie kulombowskie zaczyna dominować. W drugiej części ćwiczenia należy opisać, jak taki model pomaga zrozumieć różnicę między uranem-235 i uranem-238.

Drugie ćwiczenie, teoretyczno-przemysłowe, powinno polegać na przełożeniu tego obrazu na parametry użyteczne dla konstrukcji. Należy:

1. wyjaśnić, dlaczego samo dodatnie saldo energetyczne rozszczepienia nie wystarcza bez pokonania bariery,
2. oszacować rolę energii wzbudzenia wniesionej przez neutron,
3. porównać sytuację dla neutronów cieplnych i prędkich,
4. wskazać, jak to wpływa na wybór materiału rozszczepialnego w reaktorze i w ładunku jądrowym,
5. powiązać wynik z wymaganiami dotyczącymi przekrojów czynnych i geometrii układu krytycznego.

To ćwiczenie ma pokazać, że model kroplowy nie jest wyłącznie historyczną metaforą. Nadal porządkuje sposób myślenia o tym, kiedy ciężkie jądro tylko drga, a kiedy rzeczywiście przechodzi do nieodwracalnego rozszczepienia.

Przejdź do ćwiczenia interaktywnego

Powiązane artykuły

• Energia wiązania jądra na nukleon
• Model powłokowy jądra atomowego: liczby magiczne i spin-orbit
• Defekt masy i równoważność energii
• Pluton-239 - produkcja i właściwości
• Metoda implozyjna - mechanizm Fat Man