Energia wiązania jądra na nukleon

Streszczenie

Jeżeli narysować wykres pokazujący, jak mocno związany jest przeciętny proton lub neutron w jądrze, to bardzo szybko widać jedną rzecz: nie wszystkie jądra są równie „ciasno spięte”. Lekkie jądra, takie jak deuter, są związane słabo. Jądra o średnich liczbach masowych, w okolicy żelaza i niklu, są związane najsilniej. Bardzo ciężkie jądra, takie jak uran-235 czy pluton-239, znów schodzą niżej. To właśnie ten prosty wykres wyjaśnia zarazem sens syntezy lekkich jąder i rozszczepienia ciężkich.¹^,²

Energia wiązania jądra na nukleon nie mówi jeszcze, czy reakcja zajdzie sama. Mówi jednak, w którą stronę układ może zyskać energetycznie. Gdy produkty reakcji leżą wyżej na tej krzywej niż substraty, nadwyżka energii może zostać uwolniona jako ruch fragmentów, neutronów i promieniowania. Dlatego synteza lekkich izotopów wodoru prowadząca do helu-4 jest egzotermiczna, a rozszczepienie ciężkich aktynowców również daje dodatni bilans energetyczny.¹^,²^,³

Rozszerzenie tematu

Pojęcie energii wiązania jądra na nukleon jest najwygodniejszym skrótem do opisu stabilności jądra. Całkowita energia wiązania mówi, ile energii trzeba dostarczyć, aby rozłożyć jądro na swobodne protony i neutrony. Podzielenie tej wielkości przez liczbę masową $A$ pozwala porównywać bardzo różne jądra na jednej skali.¹ W praktyce interesuje nas więc wielkość:

$$\frac{E_w}{A}$$

gdzie $E_w$ oznacza energię wiązania, a $A$ liczbę nukleonów w jądrze. Sam mechanizm powstawania tej energii został omówiony szerzej w artykule o defekcie masy i równoważności energii, ale dopiero wykres $\frac{E_w}{A}$ pokazuje, dlaczego pewne reakcje jądrowe są opłacalne energetycznie, a inne nie.

Krzywa energii wiązania na nukleon rośnie bardzo stromo dla najlżejszych jąder. Deuter ma wartość niską, hel-4 już wyraźnie większą, a następnie wzrost zwalnia. W okolicy pośrednich liczb masowych, mniej więcej przy $A \approx 60$, pojawia się maksimum odpowiadające jądróm żelaza i niklu.¹^,² To nie jest przypadkowa własność tablicy nuklidów, tylko wynik konkurencji dwóch zjawisk. Z jednej strony działa krótkozasięgowe, lecz bardzo silne przyciąganie jądrowe między nukleonami. Z drugiej strony protony odpychają się kulombowsko, a ten składnik staje się coraz bardziej dotkliwy wraz ze wzrostem liczby protonów i rozmiaru jądra.²^,³

Dla małych jąder każdy dołożony nukleon może zyskać wielu bliskich sąsiadów i silnie obniżyć energię układu. Dlatego energia wiązania na nukleon szybko rośnie. Później zysk z dodawania kolejnych nukleonów słabnie, bo siły jądrowe mają bardzo krótki zasięg i nie „spinają” jednakowo skutecznie całego dużego jądra. W ciężkich jądrach coraz większą rolę zaczyna odgrywać odpychanie elektrostatyczne protonów, więc średnia energia wiązania na nukleon stopniowo maleje.²^,³

Na ten gładki obraz nakładają się jednak lokalne efekty struktury powłokowej. Królikowski przypomina, że niektóre jądra mają energię wiązania wyższą niż sąsiedzi właśnie dlatego, że dochodzi w nich do zamknięcia powłoki protonowej albo neutronowej, analogicznie do szczególnej trwałości gazów szlachetnych przy zamkniętych powłokach elektronowych.⁶ To nie zmienia globalnego maksimum w okolicy żelaza i niklu, ale tłumaczy, dlaczego na tle ogólnej krzywej pojawiają się lokalne „premie stabilności” związane z liczbami magicznymi i opisane szerzej w artykule o modelu powłokowym jądra atomowego.

Z tego wykresu wynikają od razu dwa najważniejsze kanały uwalniania energii jądrowej. Pierwszy to synteza lekkich jąder. Gdy z jąder słabo związanych powstaje jądro związane silniej, układ przechodzi do stanu o niższej energii wewnętrznej, a różnica musi zostać wyemitowana. W praktyce najważniejsza dla broni termojądrowej reakcja deuter-tryt prowadzi do powstania helu-4 i neutronu:

$$\mathrm{D + T \rightarrow He\text{-}4 + n + 17.588 MeV}$$

Znaczna część tego zysku bierze się stąd, że hel-4 jest jądrem wyjątkowo silnie związanym jak na tak małą liczbę nukleonów.¹^,⁴^,⁵

Drugi kanał to rozszczepienie jąder bardzo ciężkich. Uran-235, uran-233 i pluton-239 leżą już na opadającej części krzywej. Jeżeli takie jądro pęknie na dwa średnie fragmenty, produkty są przeciętnie silniej związane na każdy nukleon niż jądro wyjściowe. Stąd bierze się dodatni bilans energetyczny rozszczepienia. Gdenarz podaje, że różnica energii wiązania jest wtedy rzędu około $0{,}9 \mathrm{MeV}$ na nukleon, a Słowiński opisuje to samo zjawisko jako źródło typowego bilansu około $200 \mathrm{MeV}$ na jeden akt rozszczepienia.¹^,²

To rozróżnienie jest ważne, bo pozwala uniknąć częstego błędu interpretacyjnego. Nie jest prawdą, że ciężkie jądra „zawierają więcej energii, bo są ciężkie”. Przeciwnie: istotne jest to, że są związane mniej korzystnie niż jądra średnie. Z kolei lekkie jądra nie są atrakcyjne dlatego, że są lekkie, tylko dlatego, że część z nich może połączyć się w jądra związane silniej. Sam kierunek przemiany wynika więc z kształtu krzywej $\frac{E_w}{A}$, a nie z samej wielkości masy atomowej.¹^,²

W praktyce technologicznej oznacza to, że wykres energii wiązania jest pierwszym filtrem sensowności reakcji. Jeżeli rozważana przemiana przesuwa układ w stronę maksimum krzywej, istnieje szansa na dodatni bilans energetyczny. To jeszcze nie wystarcza do zbudowania działającego układu technicznego, bo pozostają kwestie przekrojów czynnych, barier kulombowskich, czasu reakcji i ucieczki energii. Jednak bez tego warunku wstępnego cała reszta nie miałaby znaczenia.³^,⁴

Właśnie dlatego artykuły o reakcji łańcuchowej czy masie krytycznej warto czytać dopiero po zrozumieniu krzywej energii wiązania. Te tematy odpowiadają na pytanie, jak szybko i w jakiej geometrii uwolnić energię. Krzywa $\frac{E_w}{A}$ odpowiada wcześniej na pytanie, dlaczego w ogóle jest co uwalniać.

W rozszczepieniu nie cała energia pojawia się od razu w jednej postaci. Słowiński rozdziela ją na energię kinetyczną fragmentów, energię neutronów, promieniowanie gamma, energię późniejszych rozpadów produktów rozszczepienia oraz część wynoszoną przez antyneutrina.² Ten rozkład nie zmienia jednak podstawowego faktu: całkowity bilans dodatni bierze się z przejścia z obszaru ciężkich jąder ku obszarowi średnich jąder o wyższej energii wiązania na nukleon.

W syntezie sytuacja jest podobna, choć fizycznie trudniejsza do uruchomienia. Nie wystarczy wiedzieć, że hel-4 leży wysoko na krzywej. Trzeba jeszcze doprowadzić lekkie dodatnio naładowane jądra do tak małych odległości, aby zaczęły dominować siły jądrowe nad odpychaniem kulombowskim. Stąd ekstremalne wymagania temperaturowe i gęstościowe dla reakcji termojądrowych. Kiedy jednak taki próg zostanie przekroczony, dodatni bilans energetyczny wynika z tej samej logiki co przy rozszczepieniu: produkty są związane korzystniej niż substraty.⁴^,⁵

Warto też uważać na zbyt dosłowne czytanie maksimum krzywej. To, że jądra w pobliżu żelaza i niklu są średnio najstabilniejsze, nie znaczy, że każde przejście w ich stronę będzie samoistne i łatwe. Reakcja może wymagać dostarczenia neutronu, przekroczenia bariery kulombowskiej albo przejścia przez stan wzbudzony. Krzywa energii wiązania mówi o bilansie stanu początkowego i końcowego, nie o wszystkich przeszkodach kinetycznych po drodze.²^,³

Najkrótsze podsumowanie jest więc takie: energia wiązania na nukleon porządkuje całą fizykę energetyczną reakcji jądrowych. Tłumaczy, dlaczego aktynowce mogą dawać energię przez rozszczepienie, dlaczego lekkie izotopy wodoru mogą dawać energię przez syntezę, oraz dlaczego pośrednie jądra nie są dobrym paliwem ani dla jednego, ani dla drugiego procesu. Bez tego wykresu fizyka rozszczepienia i syntezy wyglądałaby jak zbiór osobnych faktów; z nim układa się w jeden spójny obraz.¹^,²^,⁴

Historia: defekt masy i droga do wzoru Bethego-Weizsäckera

Pomysł, że jądro atomowe może być związane — że jego masa może być mniejsza niż suma mas składowych protonów i neutronów — nie pojawił się od razu po odkryciu neutronu w 1932 roku przez Jamesa Chadwicka. Ale już w 1935 roku Carl Friedrich von Weizsäcker opublikował pierwsze ilościowe przybliżenie energii wiązania jądra jako sumy pięciu składników fizycznych. Hans Bethe i Robert Bacher rozwinęli to podejście niezależnie i spopularyzowali je w monumentalnej serii artykułów Nuclear Physics w Reviews of Modern Physics w latach 1936–1937, zwanych „Biblią Bethego".⁷

Kluczowym elementem fizyki leżącej u podstaw jest równoważność masy i energii, sformułowana przez Einsteina w 1905 roku jako $E = mc^2$. W jądrze atomowym „zaginiona" masa — różnica między sumą mas wolnych nukleonów a rzeczywistą masą jądra — jest bezpośrednio proporcjonalna do energii, jaką wyzwoliłoby rozłożenie jądra na składniki. Ta różnica nosi nazwę defektu masy $\Delta m$:

$$\Delta m = Z \cdot m_p + N \cdot m_n - M(A,Z)$$

a odpowiadająca jej energia wiązania:

$$E_w = \Delta m \cdot c^2$$

Dokładny pomiar mas atomowych jest zatem bezpośrednim pomiarem energii wiązania. Pierwsze systematyczne pomiary prowadzono techniką spektrometrii masowej już od lat 20. XX wieku (Francis Aston, Nagroda Nobla 1922), ale precyzja nie była wystarczająca do budowania modeli ilościowych. Prawdziwa rewolucja nastąpiła dopiero z rozwojem pułapek Penninga w latach 80. i 90. XX wieku.⁷

Dane liczbowe dla ważnych nuklidów

Aby ożywić abstrakcyjny wykres, warto przytoczyć konkretne wartości energii wiązania na nukleon dla nuklidów kluczowych dla tej wiki:

Nuklidium	$Z$	$A$	$E_w$ [MeV]	$E_w/A$ [MeV/nukl.]
Deuter (H-2)	1	2	2,225	1,11
Tryt (H-3)	1	3	8,482	2,83
Hel-3	2	3	7,718	2,57
Hel-4	2	4	28,296	7,07
Lit-6	3	6	31,995	5,33
Węgiel-12	6	12	92,162	7,68
Żelazo-56	26	56	492,259	8,79
Nikiel-62	28	62	545,259	8,79
Cyrkon-90	40	90	783,897	8,71
Cyna-120	50	120	1020,540	8,50
Bar-138	56	138	1158,304	8,39
Ołów-208	82	208	1636,430	7,87
Rad-226	88	226	1731,610	7,66
U-235	92	235	1783,870	7,59
U-238	92	238	1801,695	7,57
Pu-239	94	239	1806,920	7,56

Dane pochodzą z bazy AME2020 (Atomic Mass Evaluation 2020), opublikowanej przez Huang et al. w Chinese Physics C 45, 030002–030003 (2021). Jest to najnowsza systematyczna kompilacja mas atomowych, oparta na pomiarach z ponad 100 laboratoriów.⁷

Widać wyraźnie: hel-4 ma $E_w/A = 7{,}07 \mathrm{MeV}$, co jest wyjątkowo duże jak na tak małe jądro — stąd cząstka alfa jest tak stabilna i tak często emitowana w rozpadach ciężkich jąder. Nikiel-62 i żelazo-56 osiągają maksimum ~8,79 MeV/nukl. Dla U-235 wartość spada do 7,59 MeV/nukl. Różnica na nukleon między U-235 a typowym produktem rozszczepienia (Ba-138: 8,39 MeV/nukl.) wynosi zatem $\approx 0{,}8 \mathrm{MeV}$ — co po pomnożeniu przez 235 nukleonów daje $\approx 188 \mathrm{MeV}$ na rozszczepienie. To dobrze zgadza się z mierzonym bilansem ~200 MeV.¹^,²

Wzór Bethego-Weizsäckera: rozwinięcie składników

Wzór Bethego-Weizsäckera (SEMF — Semi-Empirical Mass Formula) daje przybliżoną energię wiązania:

$$E_w(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$$

Gdzie:

$a_V \approx 15{,}8 \mathrm{MeV}$ — składnik objętościowy (siły jądrowe nasycone)
$a_S \approx 18{,}3 \mathrm{MeV}$ — składnik powierzchniowy (nukleony na powierzchni mają mniej sąsiadów)
$a_C \approx 0{,}714 \mathrm{MeV}$ — składnik kulombowski (odpychanie protonów)
$a_A \approx 23{,}2 \mathrm{MeV}$ — składnik asymetrii (kara za nierówną liczbę p i n)
$\delta$ — składnik parzystości: $+12/\sqrt{A}$ dla pp, 0 dla nieparzystego A, $-12/\sqrt{A}$ dla nn

Wzór reprodukuje ogólny kształt krzywej z dokładnością ~1 MeV/nukl., ale nie oddaje lokalnych fluktuacji przy liczbach magicznych. Dlatego jest „semi-empiryczny" — parametry dopasowuje się do danych, a nie wyprowadza z pierwszych zasad.⁷

Szczególnie ważny jest składnik parzystości δ. Jądra parzysto-parzyste (pp — parzysty Z i parzysty N) mają wyższą energię wiązania niż jądra nieparzysto-nieparzyste (nn). To dlatego większość trwałych jąder jest pp: z 274 trwałych nuklidów, 148 to pp, tylko 5 to nn. Ma to bezpośrednie konsekwencje dla fizyki reaktorów: U-235 jest nieparzyste, co zwiększa jego energię wiązania neutronu (co prowadzi do wyższej energii wzbudzenia U-236* i możliwości rozszczepienia przez neutrony termiczne).

Linia doliny stabilności i granice istnienia jąder

Nuklidy można umieścić na mapie $N$-$Z$ (mapa nuklidów, mapa Ségré). Na tej mapie istnieje wąski pas zwany doliną stabilności — jądra leżące w tym pasie mają wystarczająco wysoką energię wiązania, aby być trwałe lub co najmniej długotrwałe.

Po lewej stronie doliny (za dużo protonów) jądra rozpadają się przez emisję pozytonu lub wychwyt elektronu. Po prawej stronie (za dużo neutronów) rozpadają się przez emisję elektronu (β⁻). Poza granicami doliny, przy ekstremalnych stosunkach N/Z, leżą tzw. linie kapania (drip lines):

Proton drip line: granica, za którą kolejny proton nie jest już związany (energia wiązania ostatniego protonu = 0). Przekroczenie tej granicy prowadzi do natychmiastowej emisji protonu.
Neutron drip line: analogicznie dla neutronów. Dla lekkich jąder jest znana eksperymentalnie (np. hel-8 to jądro halo z 6 neutronami i 2 protonami, leżące blisko granicy neutronowej). Dla ciężkich jąder jest jeszcze niezbadana eksperymentalnie — szacowania teoretyczne wskazują, że neutronowa drip line dla uranu leży gdzieś przy A~300.

Jądra przy neutronowej drip line mają niezwykłe właściwości — mogą mieć tzw. halo jądrowe, gdzie jedna lub dwie pary neutronów orbitują daleko poza rdzeniem jądrowym. Klasycznym przykładem jest lit-11 (Li-11), który ma promień jądrowy zbliżony do ołowiu-208 pomimo zaledwie 11 nukleonów. Takie jądra halo mają wyjątkowo niską energię wiązania ostatnich neutronów ($\sim 0{,}1 \mathrm{MeV}$) i są laboratoriami do badania struktury jądrowej w ekstremalnych warunkach.⁸

Synteza pierwiastków we Wszechświecie i nukleosynteza

Krzywa energii wiązania ma swoje bezpośrednie odbicie w składzie chemicznym Wszechświata. Pierwiastki do żelaza ($A \leq 56$) mogą być syntetyzowane przez fuzję w gwiazdach — każdy krok ku maksimum krzywej wyzwala energię. Stąd reakcje jądrowe są źródłem energii słonecznej przez miliardy lat.

Dla pierwiastków cięższych od żelaza potrzebne są procesy, w których zewnętrzna energia jest dostarczana — np. eksplozja supernowej. W gwiazdach neutronowych i podczas zderzeń gwiazd neutronowych pierwiastki cięższe od żelaza powstają przez:

Proces s (slow neutron capture): powolny wychwyt neutronów, z czasem życia między kolejnymi wychwytami dłuższym niż czas β⁻ rozpadów. Daje pierwiastki od Fe do Ba i Pb.
Proces r (rapid neutron capture): bardzo szybki wychwyt neutronów w środowiskach o ekstremalnym strumieniu ($>10^{20} \mathrm{n/cm^2/s}$), z krótkim czasem między wychwytami. Daje pierwiastki cięższe, w tym aktynowce.
Proces p (proton capture lub foto-dysocjacja): odpowiada za ~35 nuklidów po lewej stronie linii stabilności, niedostępnych dla procesów s i r.⁸

Pięknym potwierdzeniem tej teorii jest to, że wykres obfitości pierwiastków w Układzie Słonecznym ma lokalne piki dokładnie tam, gdzie model powłokowy przewiduje magiczne liczby neutronów — bo jądra z zamkniętymi powłokami wolniej absorbują dodatkowe neutrony, więc gromadzą się jako produkty procesów s i r.

Obserwacja GW170817 (sierpień 2017) — zderzenie dwóch gwiazd neutronowych potwierdzone zarówno falami grawitacyjnymi (LIGO/Virgo), jak i promieniowaniem elektromagnetycznym — potwierdziła, że to właśnie te zderzenia produkują większość złota, platyny i uranu we Wszechświecie przez proces r.⁸

Hel-4 jako punkt szczególny na krzywej

Hel-4 zasługuje na oddzielne omówienie, bo jego wyjątkowa pozycja na krzywej ma przełomowe konsekwencje praktyczne. Przy $A=4$ energia wiązania sięga 7,07 MeV/nukl. — wyraźnie ponad lokalnym trendem dla małych jąder. Przyczyną są liczby magiczne: w He-4 zarówno protony (2), jak i neutrony (2) mają zamknięte powłoki 1s. To sprawia, że He-4 jest wyjątkowo zwartym i stabilnym konfiguracją.

Praktyczne konsekwencje:

Rozpad alfa: ciężkie jądra preferowalnie emitują właśnie cząstki alfa (jądra He-4), bo emisja właśnie tej konfiguracji jest energetycznie korzystna. Emisja protonu lub deuteronu byłaby mniej korzystna energetycznie z powodu niżej leżących wartości $E_w/A$ tych lżejszych obiektów.
Fuzja termojądrowa D+T: produktem jest He-4 (7,07 MeV/nukl.) plus neutron. Zysk energetyczny 17,588 MeV per reakcję bierze się właśnie z wyjątkowej stabilności He-4 na tle substratów D (1,11 MeV/nukl.) i T (2,83 MeV/nukl.).
Niekwenkonwencjonalne rozszczepienia: tzw. rozszczepienie ternarne (ternary fission) — emisja cząstki alfa z przedziału między dwoma głównymi fragmentami — zachodzi z prawdopodobieństwem ~0,3% dla U-235. Jest to bezpośredni efekt stabilności He-4 jako konfiguracji nukleońskiej.
Reakcja tryptonu w eksplodującej broni: w ładunkach termojądrowych boosted fission (np. W87) neutrony z reakcji D+T znacząco podwyższają liczbę neutronów wywołujących rozszczepienia — efekt jest możliwy właśnie dlatego, że D+T daje 17,588 MeV i bardzo energetyczny neutron (14,1 MeV).⁴

Pomiary mas atomowych: spektrometria masowa i pułapki Penninga

Energia wiązania jest pochodną masy atomowej — i to bardzo wrażliwą: błąd w ósmym miejscu po przecinku masy (w jednostkach u) przekłada się na błąd rzędu 1 keV w energii wiązania. Dlatego precyzja pomiaru mas jest kluczowa.

Historycznie stosowano spektrometria masową (Francis Aston, 1919 i późniejsi): wiązka jonów przechodzi przez pola elektryczne i magnetyczne, które odchylają ją w zależności od stosunku ładunek/masa. Osiągane precyzje: $\Delta m/m \sim 10^{-6}$ w klasycznych urządzeniach.

Rewolucja nastąpiła wraz z pułapkami Penninga (Penning ion trap, PIT), gdzie jon jest uwięziony między polem elektrycznym (wzdłużnym) i magnetycznym (poprzecznym), a jego masa mierzona przez wyznaczenie częstości cyklotronowej $\omega_c = qB/m$. Współczesne pułapki jak ISOLTRAP (CERN/ISOLDE), SHIPTRAP (GSI Darmstadt) czy JYFLTRAP (Jyväskylä) osiągają $\Delta m/m \sim 10^{-9}$, co odpowiada ~1 keV dla jądra o masie 100 u.⁷

Szczególną rolę odgrywa instalacja ISOLDE w CERN, gdzie pęczki radioaktywnych jąder (produkowanych bombardowaniem grubej tarczy protonami z PS Booster) są oddzielane, jonizowane i kierowane do pomiarów masowych. Właśnie tam zbadano wiele jąder na linii drip, izomerów i jąder halo, rozszerzając bazę AME2020 o kilkaset nowych punktów pomiarowych.

Energia wiązania a proliferacja i polityka nieproliferacji

Znajomość krzywej energii wiązania ma nie tylko dydaktyczne, ale i bezpośrednie implikacje polityczne. To właśnie ten wykres uzasadnia, dlaczego:

U-235 i Pu-239 są atrakcyjne jako materiały broni: leżą na opadającej gałęzi krzywej, znacznie poniżej żelaza, więc ich rozszczepienie wyzwala duże ilości energii.
LEU (nisko wzbogacony uran) jest mniej atrakcyjny: paliwo reaktorowe zawiera typowo 3–5% U-235, a reszta to U-238. Masa krytyczna dla LEU jest wielokrotnie większa niż dla HEU (wysoko wzbogaconego), a poza tym stosunek produktywności rozszczepień do tłumienia przez U-238 jest niekorzystny dla projektu broni.
Tor-232 i U-233 są alternatywnymi materiałami: U-233 z rozszczepienia torowego leży na tej samej opadającej krzywej co U-235, z podobnym $E_w/A$. Dlatego cykle torowe są potencjalnie proliferacyjnie wrażliwe i podlegają takim samym gwarancjom MAEA.

Program IAEA Nuclear Security Series (NSS) i Traktat NPT opierają swoje definicje materiałów wrażliwych właśnie na fizycznym kryterium fisylności — które jest pochodną pozycji jądra na krzywej energii wiązania.⁹

Energia wiązania w reaktorach: zmiany składu paliwa

W działającym reaktorze energia wiązania „pracuje" nieustannie, a skład paliwa zmienia się przez cały cykl. Świeże paliwo uranowe zawiera typowo ~4,5% U-235, reszta to U-238. W trakcie pracy:

U-235 ulega rozszczepieniu, wytwarzając energię
U-238 pochłania neutrony i przez kolejne przemiany ($\beta^-$) zamienia się w Pu-239 i Pu-241
Produkty rozszczepienia akumulują się i pochłaniają neutrony (trucizny jądrowe: Xe-135, Sm-149)
Izotopy ciężkie wyższe (Am, Cm, Np) narastają

Po wypaleniu paliwa do ~50 GWd/t (gigawatod na tonę) skład paliwa zawiera ~0,8% U-235, ~0,6% Pu (mix izotopów), 3–4% produktów rozszczepienia i ok. 95% U-238. Cały ten proces jest deterministycznie wyznaczony przez przekroje czynne i energię wiązania produktów kolejnych rozpadów. Obliczenia ewolucji składu paliwa (fuel depletion, burnup calculation) są przeprowadzane kodami ORIGEN (ORNL), FISPACT i SERPENT2, które numerycznie rozwiązują setki sprzężonych równań różniczkowych dla setek nuklidów.⁹

Wypalanie paliwa ma bezpośredni związek z proliferacją: wypalone paliwo jądrowe zawiera Pu-239 i Pu-241. MAEA wymaga od każdej elektrowni deklarowania ilości wypalonego paliwa i jego składu, bo z paliwa wypalonego można chemicznie wydzielić pluton (PUREX process). Krzywa energii wiązania wyznacza, ile plutonu wytworzyło się podczas pracy reaktora — i to ilości, którą MAEA musi zliczyć.

Reaktory nowej generacji i krzywa energii wiązania

Nowoczesne projekty reaktorów (Generation IV) opierają się na rozszerzeniu zakresu nuklidów używanych jako paliwo — co z punktu widzenia fizyki oznacza eksplorację nowych obszarów krzywej energii wiązania:

Reaktory prędkie (SFR, GFR): używają neutronów o energii ~1 MeV. W tym zakresie U-238 ma znaczące przekroje czynne na rozszczepienie (σ_f ~0,5 b przy 1 MeV) — w przeciwieństwie do reaktorów termicznych, gdzie jest praktycznie nierozszczepialny. Reaktory prędkie mogą więc „spalać" U-238 i długożyciowe transurany (Np, Am, Cm), przekształcając je w produkty rozszczepienia o znacznie krótszym czasie życia. To tzw. transmutacja — a jest możliwa właśnie dlatego, że te nuklidy leżą na opadającej gałęzi krzywej i ich rozszczepienie jest energetycznie korzystne.

Reaktory torowe (MSR, AHWR): U-233 produkowany z Th-232 przez wychwyt neutronu i dwa kolejne rozpady β⁻. U-233 ma podobną fisylność jak U-235 i Pu-239, a Th-232 jest ~3× bardziej obfity w przyrodzie niż uran. Kraje posiadające duże zasoby torowe (Indie, Brazylia, Norwegia) intensywnie badają ten cykl.

SMR (Small Modular Reactors): wiele projektów SMR używa HALEU (High-Assay Low-Enriched Uranium, 5–20% U-235). Wyższe wzbogacenie pozwala na mniejsze rdzenie i dłuższe cykle, co jest ekonomicznie atrakcyjne. HALEU wytwarza więcej Pu-239 na jednostkę masy niż standardowe LEU (bo U-238 jest mniej i ma mniej okazji do pochłaniania neutronów bez rozszczepienia), co wymaga dodatkowego nadzoru MAEA.⁹

Bilans energetyczny broni termojądrowej a krzywa wiązania

W broni termojądrowej krzywa energii wiązania odgrywa rolę dwukrotnie — raz dla materiału rozszczepialnego (rdzeń rozszczepieniowy), raz dla paliwa fuzyjnego. W klasycznym układzie Teller-Ulam:

Rdzeń rozszczepieniowy (Pu-239 lub HEU): energia bilansu $\approx 0{,}8 \mathrm{MeV/nukl.}$ × 239 nukleonów = ~190 MeV per rozszczepienie. Generuje promieniowanie X i neutrony inicjujące fuzję.
Paliwo fuzyjne (Li-6D, czyli lit-6 deuteryd): w obecności neutronów Li-6 ulega rozpadowi do He-4 i trytu: $n + \mathrm{Li\text{-}6} \rightarrow \mathrm{He\text{-}4} + T + 4{,}8 \mathrm{MeV}$. Następnie T i D reagują: $D + T \rightarrow \mathrm{He\text{-}4} + n + 17{,}6 \mathrm{MeV}$.

Energia fuzji na parę D+T wynosi 17,6 MeV. Dla porównania, energia rozszczepienia jednego atomu U-235 to ~200 MeV, ale masa nuklidu jest 235 razy większa niż masy reagentów w fuzji. Na kilogram materiału, fuzja D+T daje $3{,}4 \times 10^{14} \mathrm{J/kg}$ wobec $8{,}2 \times 10^{13} \mathrm{J/kg}$ dla rozszczepienia U-235 — fuzja jest ~4× bardziej energetyczna masowo. To dlatego bronie termojądrowe mogą osiągać wydajności rzędu megaton, nieosiągalne dla czystych ładunków rozszczepialnych.⁴^,⁵

Krzywa energii wiązania w edukacji jądrowej

Krzywa $E_w/A$ jako funkcja $A$ jest chyba najważniejszym pojedynczym wykresem w całej fizyce jądrowej. Można to stwierdzenie uzasadnić argumentem: zna się jeden wykres — rozumie się zasadę działania reaktorów i broni jądrowej, rozumie się genezę pierwiastków we Wszechświecie i granice istnienia jąder atomowych.

W polskim systemie edukacyjnym fizyka jądrowa pojawia się na poziomie rozszerzonym matury, ale głębokość jej omówienia bywa niewystarczająca. Podręczniki Słowińskiego i Gdernarz wprowadzają energię wiązania, ale rzadko sięgają do tablic AME i pełnego wzoru Bethego-Weizsäckera. Programy akademickie (AGH, Politechnika Warszawska, Wydział Fizyki UW) omawiają ten temat szczegółowo w kursach fizyki jądrowej i reaktorowej, często korzystając z podręczników Krane'a (Introductory Nuclear Physics, 1988) lub Mukhopadhyay'a.

Dla doktorantów specjalizujących się w fizyce jądrowej — do których adresowana jest ta wiki — kluczowe jest połączenie intuicji z wykresu z precyzją numeryczną AME2020 i zdolnością do interpretacji odchyleń od krzywej SEMF jako sygnałów struktury powłokowej. To właśnie te odchylenia — lokalne piki stabilności przy liczbach magicznych $N,Z = 2, 8, 20, 28, 50, 82, 126$ — są tematem artykułu o modelu powłokowym jądra atomowego.

{{image-full:https://upload.wikimedia.org/wikipedia/commons/thumb/5/53/Binding_energy_curve_-_common_isotopes.svg/800px-Binding_energy_curve_-_common_isotopes.svg.png:Krzywa energii wiązania na nukleon jako funkcja liczby masowej A. Widoczne maksimum przy A≈56 (żelazo-56, nikiel-62) oraz opadanie dla ciężkich aktynowców. Deuter leży wyjątkowo nisko (1,11 MeV/nukl.), hel-4 wykazuje lokalny skok ku górze (7,07 MeV/nukl.) jako efekt zamkniętych powłok. Rozszczepienie ciężkich jąder (np. U-235) i fuzja lekkich (D+T→He-4) oba przesuwają produkty w kierunku maksimum, wyzwalając energię.}}

Granica nuklearna: od izotopu do kwarkowego plazmu

Krzywa energii wiązania opisuje strukturę jąder złożonych z protonów i neutronów traktowanych jako odrębne składowe. Ta perspektywa obowiązuje w zakresie energii od keV do kilkuset MeV na nukleon. W skrajnych warunkach — przy energiach $>1 \mathrm{GeV/nukl.}$, osiągalnych na akceleratorach takich jak LHC czy RHIC — materia jądrowa przechodzi do stanu kwarkowo-gluonowej plazmy (QGP), gdzie poszczególne nukleony przestają istnieć jako odrębne obiekty.

W tym stanie energia wiązania w tradycyjnym sensie traci sens. Fizyka QCD (kwantowej chromodynamiki) zastępuje fizykę jądrową. Zderzenia jąder ołowiu przy LHC w trybie Pb+Pb (@5,02 TeV/parę nukleonów) pozwalają badać QGP eksperymentalnie — jest to temperatura $>10^{12} \mathrm{K}$, istniejąca tylko przez $\sim 10^{-23} \mathrm{s}$.

Dla kontekstu tej wiki: astrofizyczne środowiska gwiazd neutronowych zbliżają się do warunków QGP, choć temperatura jest tu niższa (~10⁹ K), a dominującą rolę odgrywa gęstość. Masa-energia wiązania jądra atomowego jest zatem tylko dolnym obszarem bardzo rozległej mapy stanu materii jądrowej.⁷

Rozszczepienie ternarne i egzotyczne kanały rozpadu

Standardowym wyjściem z rozszczepienia są dwa fragmenty (binary fission). Jednak z prawdopodobieństwem ~0,3% rozszczepienie U-235 neutronami termicznymi daje trzy produkty (ternary fission) — dwa ciężkie fragmenty i cząstkę alfa. Z jeszcze mniejszym prawdopodobieństwem mogą pojawić się deuter, tryt lub $^3$He jako trzeci produkt.

Rozszczepienie ternarne z emisją alfa jest możliwe właśnie dlatego, że alfa leży bardzo wysoko na krzywej ($E_w/A = 7{,}07 \mathrm{MeV}$) — jest więc konfiguracyjnie korzystna. Gdyby He-4 nie był tak wyjątkowo stabilny, ternarne rozszczepienie z emisją alfa byłoby znacznie mniej prawdopodobne.

Istnieją też egzotyczne kanały, obserwowane przez spektrometry masowe dołączone do instalacji ISOLDE i GANIL:

Emisja He-6, He-8 z rozszczepień ciężkich aktynowców
Emisja C-14 z Ra-226 (słynny przypadek emisji ciężkich klastrów)
Clusterowa radioaktywność — emisja $^{14}$C, $^{20}$O, $^{24}$Ne z jąder z grupy Ra-Th-U

Wszystkie te kanały można rozumieć w kategoriach krzywej energii wiązania: emitowany klaster i pozostałe jądro muszą być konfiguracyjnie korzystniejsze niż jądro macierzyste. Mapa ta pozwala przewidywać, które klastry mogą być emitowane, a które nie.⁸

Depozyty masy i energia wyzwolona w rozpadach łańcuchowych

Po rozszczepieniu fragmenty nie są od razu stabilne — są wyjątkowo bogate w neutrony (bo jądro macierzyste miało proporcję N/Z charakterystyczną dla ciężkich jąder, ~1,55 dla U-235, ale produkty o A~90-140 mają optymalną N/Z znacznie niżej, ~1,2). Stąd fragmenty rozszczepienia ulegają długim łańcuchom rozpadów β⁻, emitując elektrony i antyneutrina.

Całkowita energia rozpadów β⁻ produktów rozszczepienia to ~13 MeV na jedno zdarzenie rozszczepienia — to część bilansu energetycznego, która jest wyzwalana po rozszczepieniu, nie w jego trakcie. Ten składnik jest ważny dla termiki paliwa wypalonego i obliczeń tzw. decay heat (ciepła porozpadowego), które wyznacza minimalne chłodzenie potrzebne po wyłączeniu reaktora.

W katastrofie w Fukushimie (marzec 2011) głównym problemem było właśnie ciepło porozpadowe: reaktory zostały wyłączone natychmiast po trzęsieniu ziemi (pręty kontrolne weszły automatycznie), ale energia β⁻ rozpadów produktów rozszczepienia nadal podgrzewała rdzeń. Bez chłodzenia (pompy pompujące wodę zatrzymały się po zalaniu przez falę tsunami) temperatura rosła — i doprowadziła do stopienia rdzenia. To fundamentalne zachowanie jest wyznaczone przez krzywe rozpadu setek izotopów, ostatecznie zakorzenione w tym, gdzie leżą one na tle krzywej energii wiązania.⁹

Paliwo fuzyjne vs rozszczepienie: porównanie energetyczne

Zestawienie energii dostępnej w czystej fuzji i w czystym rozszczepieniu jest ważne dla zrozumienia różnicy możliwości obu procesów:

Materiał	Reakcja	Energia [MeV]	Energia [MJ/kg]
D + T	D + T → He-4 + n	17,59	337 000
D + D	D + D → He-3 + n (50%) / T + p (50%)	~3,6	84 000
D + ³He	D + ³He → He-4 + p	18,35	352 000
U-235	rozszczepienie	~200	83 000
Pu-239	rozszczepienie	~210	85 000
Li-6D (lit-6 deuteryd)	pełny cykl	~22 (D+T+4,8)	390 000
Węgiel (spalanie)	C + O₂ → CO₂	0,000411 eV	~33

Fuzja D+T daje ~337 000 MJ/kg, rozszczepienie U-235 ~83 000 MJ/kg. Spalanie węgla: ~33 MJ/kg. Fuzja jest ok. 4× bardziej energetyczna masowo niż rozszczepienie i ~10 000× bardziej niż spalanie.⁴^,⁵

Jednak reaktory fuzyjne nadal nie istnieją w fazie komercyjnej (NIF w grudniu 2022 osiągnął po raz pierwszy ignition — Q > 1, czyli więcej energii z fuzji niż dostarczone przez lasery — ale daleko mu do opłacalności komercyjnej). Elektrownie reaktorów fuzyjnych są prognozowane na lata 2040–2060 (ITER, ARC, SPARC). Reaktory rozszczepienia są za to dojrzałą technologią komercyjną od lat 50. — właśnie dlatego, że przekrój czynny na rozszczepienie U-235 neutronami termicznymi (584 b) sprawia, że reakcja jest inicjowalna w dostępnych warunkach, bez wymogu temperatury 100 milionów kelwinów.

Trwałość na skraju krzywej: nuklidy super- i superheavy

Poza ołowiem i bizmutem (koniec doliny stabilności dla trwałych nuklidów) istnieje szerokie spektrum jąder syntetycznych, wyprodukowanych w akceleratorach. Najcięższe znane nuklidy (oganesson Z=118, A=294) mają czas życia rzędu milisekund. Teoretycznie jednak istnieje przewidywana „wyspa stabilności" wokół $Z = 114–126$, $N = 184$ — konfiguacja, w której efekty powłokowe miałyby dawać znaczącą premię stabilności, analogiczną do tej, która sprawia, że Pb-208 ($Z=82, N=126$) jest wyjątkowo trwały.

Poszukiwania wyspy stabilności to jeden z frontów badań fizyki jądrowej — program FUSION (GSI/FAIR Darmstadt), laboratorium Dubna (JINR), RIKEN Nishina Center w Japonii i in. Gdyby wyspa istniała i gdyby jej nuklidy miały czasy życia rzędu sekund lub minut, otwierałoby to możliwość chemicznych badań pierwiastków transoganessonu — co byłoby pierwszym wglądem w chemię pierwiastków z konfigurować elektronowymi daleko poza tablicą Mendelejewa.

Z punktu widzenia krzywej energii wiązania nuklidy wyspy stabilności leżałyby dość wysoko na wykresie — jeszcze wyżej niż otaczające je nuklidy — bo premia powłokowa podnosiłaby ich $E_w/A$ powyżej płynnego trendu krzywej LDM. To byłby bezpośredni analogon magicznych jąder w środku tablicy — ale na samym jej skraju.⁷^,⁸

Podsumowanie: jeden wykres, wiele zastosowań

Krzywa energii wiązania na nukleon jest punktem startowym dla zrozumienia praktycznie każdego obszaru fizyki jądrowej:

Reaktory: bilans energetyczny rozszczepienia wynika z kształtu krzywej; składnik parzystości wyjaśnia fisylność U-235; ewolucja składu paliwa odpowiada przesuwaniu nuklidów po mapie ku minimum energii.
Broń jądrowa: zarówno ładunki rozszczepieniowe, jak i termojądrowe, wykorzystują przejście ku wyżej leżącym na krzywej jądrom produktów.
Nieproliferacja: definicje materiałów wrażliwych NPT i Significant Quantity MAEA są zakorzenione w mierzonej pozycji aktynowców na krzywej.
Astrofizyka: nukleosynteza pierwiastków od H do U jest w całości opowiedziana przez przechodzenie między stanami na krzywej.
Fizyka fundamentalna: odchylenia od gładkiej krzywej SEMF są sygnałami struktury powłokowej i kolektywnej — informacją o kwantowych korelacjach w silnie oddziałujących układach złożonych.

Ćwiczenia numeryczne dla doktorantów powinny zawsze zaczynać się od odtworzenia tej krzywej z danych AME2020 i porównania z przewidywaniami wzoru Bethego-Weizsäckera — to jest „kalibracja" oka dla każdego, kto ma rzetelnie rozumieć fizykę jądrową.

Dodatkowe materiały multimedialne

Do tego artykułu nie dodano jeszcze materiałów wideo. Warto wrócić do tej sekcji dopiero wtedy, gdy uda się znaleźć materiał dobrze pokazujący samą krzywą energii wiązania i odróżniający bilans energetyczny reakcji od ich trudności zapłonu.

Jako rozwinięcie warto od razu przejść do defektu masy i równoważności energii, rozszczepienia w modelu kroplowym, modelu powłokowego jądra atomowego oraz roli helu-4 i litu-6 w syntezie. Wtedy krzywa energii wiązania przestaje być samotnym wykresem, a staje się mapą łączącą obie wielkie klasy reakcji jądrowych oraz wewnętrzną strukturę samych jąder.

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 obliczeniowego krzywej energii wiązania na nukleon na podstawie rzeczywistych danych masowych. W wariancie podstawowym należy przygotować skrypt, który dla wybranego zestawu nuklidów:

wczyta liczbę protonów $Z$, liczbę neutronów $N$ i masę jądra,
obliczy defekt masy i całkowitą energię wiązania,
przeliczy $\frac{E_w}{A}$,
narysuje wykres w funkcji liczby masowej $A$,
zaznaczy przynajmniej deuter, hel-4, jądra w okolicy żelaza i niklu oraz uran-235.

Celem ćwiczenia nie jest samo wykonanie wykresu, lecz sprawdzenie, czy z rzeczywistych danych da się odtworzyć trzy kluczowe cechy: stromy wzrost dla najlżejszych jąder, maksimum w pobliżu $A \approx 60$ oraz spadek dla ciężkich aktynowców. Druga część tego samego ćwiczenia powinna polegać na policzeniu różnicy energii wiązania między jądrem wyjściowym a produktami dla jednego przykładowego rozszczepienia i jednej reakcji fuzyjnej prowadzącej do hel-4.

Drugie ćwiczenie, teoretyczno-przemysłowe, powinno dotyczyć skali energetycznej. Na podstawie własnego modelu i danych tablicowych należy:

oszacować energię uwalnianą przez pojedynczy akt rozszczepienia uranu-235,
przeliczyć ją na energię z 1 g i 1 kg materiału przy założeniu pełnego rozszczepienia,
porównać wynik z energią spalania paliwa chemicznego,
wykonać analogiczne oszacowanie dla reakcji $D + T$,
wskazać, jaka część różnicy między teorią a praktyką wynika nie z bilansu energetycznego, lecz z trudności realizacyjnych.

To ćwiczenie ma nauczyć rozróżniania dwóch poziomów opisu. Krzywa energii wiązania mówi, które przemiany są energetycznie korzystne. Dopiero osobne modele neutroniczne, hydrodynamiczne i cieplne mówią, czy taki bilans da się zamienić w działające urządzenie.

Przejdź do ćwiczenia interaktywnego

Powiązane artykuły

Przypisy

Jeremi Gdenarz, Jak się bać broni jądrowej. Zakres wykorzystania: wykres energii wiązania na nukleon, przykład hel-4 jako jądra silnie związanego, maksimum stabilności dla jąder pośrednich oraz oszacowanie różnicy rzędu około $0{,}9 \mathrm{MeV}$ na nukleon między ciężkimi jądrami a produktami rozszczepienia.
↩
Bronisław Słowiński, Podstawy fizyczne energetyki jądrowej. Zakres wykorzystania: interpretacja maksimum krzywej w okolicy żelaza i niklu, związek energii wiązania z egzotermicznością rozszczepienia oraz bilans około $200 \mathrm{MeV}$ energii jednego aktu rozszczepienia.
↩
Wprowadzenie do fizyki i konstrukcji broni jądrowej. Zakres wykorzystania: ogólna interpretacja stabilności jąder ciężkich, rozróżnienie jąder wolnorozszczepialnych i szybkorozszczepialnych oraz sens energetyczny rozszczepienia ciężkich aktynowców.
↩
Hybrydowe bronie rozszczepieniowo-fuzyjne. Zakres wykorzystania: rola reakcji $D + T$ jako podstawowego praktycznego kanału syntezy dla zastosowań militarnych oraz energia uwalniana w tej reakcji.
↩
Elementy projektowania broni termojądrowej. Zakres wykorzystania: zestaw głównych reakcji fuzyjnych, w tym $D + T$ i $He\text{-}3 + D$, oraz rozkład energii produktów reakcji prowadzących do hel-4.
↩
J. Królikowski, Wstęp do Fizyki Jądra Atomowego i cząstek elementarnych. VI. Wzbudzenia jąder. Model powłokowy. Deformacje. Zakres wykorzystania: uwaga o lokalnie podwyższonej energii wiązania jąder magicznych i analogii do zamkniętych powłok elektronowych.
↩
W.-J. Huang i in., The AME 2020 atomic mass evaluation (I), Chinese Physics C 45, 030002 (2021); H. A. Bethe, R. F. Bacher, Nuclear Physics A. Stationary States of Nuclei, Reviews of Modern Physics 8, 82–229 (1936). Zakres wykorzystania: historia wzoru Bethego-Weizsäckera, dane liczbowe energii wiązania, parametry SEMF, pułapki Penninga, AME2020.
↩
A. Bohr, B. R. Mottelson, Nuclear Structure vol. I (1969); B. A. Brown, The Nuclear Shell Model, Progress in Particle and Nuclear Physics 47, 517–599 (2001); R. Surman i in., r-Process Nucleosynthesis in Core-Collapse Supernovae, Annual Review of Nuclear and Particle Science 64, 481–486 (2014). Zakres wykorzystania: linia drip, jądra halo, nukleosynteza procesami r/s/p, zderzenia gwiazd neutronowych GW170817, obfitości pierwiastków.
↩
IAEA Nuclear Security Series No. 20-G, Nuclear Security Recommendations on Nuclear and Other Radioactive Material out of Regulatory Control (2011); IAEA Safeguards Glossary (wyd. 2022); D. J. Loaiza (ed.), Nuclear Criticality Safety (ANS, 2021). Zakres wykorzystania: definicje materiałów wrażliwych NPT/MAEA, HALEU, wypalanie paliwa, ORIGEN burnup, reaktory nowej generacji i proliferacja.
↩