Model powłokowy jądra atomowego: liczby magiczne i spin-orbit

Streszczenie

Jeżeli spojrzeć tylko na model kroplowy, jądro atomowe wygląda jak zwarta kropla materii jądrowej, która może drgać, deformować się i czasem pękać. Taki obraz jest bardzo użyteczny, ale nie tłumaczy wszystkiego. Nie wyjaśnia na przykład, dlaczego pewne liczby protonów i neutronów dają jądrom wyjątkową trwałość, dlaczego poziomy wzbudzeń zmieniają się skokowo przy niektórych obsadzeniach ani czemu jądra dalekie od ścieżki stabilności potrafią mieć inny układ stanów niż ich stabilni krewniacy.¹

Model powłokowy odpowiada właśnie na te pytania. Traktuje protony i neutrony nie jako bezkształtną ciecz, lecz jako nukleony zajmujące dyskretne orbitale w średnim potencjale jądra. Dzięki temu można zrozumieć liczby magiczne, rolę sprzężenia spin-orbit oraz to, dlaczego w jądrach neutrononadmiarowych mogą pojawiać się nowe zamknięcia podpowłok, takie jak N=32, mimo że klasyczny zestaw liczb magicznych był wyprowadzony dla jąder bliższych stabilności.¹

Rozszerzenie tematu

Najprostsza intuicja jest taka: ciężkie jądro nie jest wyłącznie jedną bryłą materii, ale także układem wielu protonów i neutronów poruszających się po dozwolonych stanach kwantowych. W praktyce oznacza to, że dwa poziomy opisu działają równocześnie. Z jednej strony jądro może zachowywać się kolektywnie, jak w modelu kroplowym. Z drugiej strony jego własności zależą od tego, które orbitale są obsadzone, które pozostają puste i jak daleko energetycznie leżą od siebie sąsiednie stany.¹

To rozróżnienie ma bardzo praktyczny sens. W ciężkich aktynowcach obraz kroplowy dobrze tłumaczy sam mechanizm deformacji i rozszczepienia. Gdy jednak chcemy zrozumieć, dlaczego jedno jądro ma wyjątkowo wysoki pierwszy stan 2+, a inne nie, albo dlaczego pewne izotopy są bardziej „sztywne” niż sąsiednie, potrzebujemy języka powłok, orbitali i szczelin energetycznych. To właśnie na tej warstwie opiera się pojęcie liczb magicznych.¹

Klasyczne liczby magiczne to liczby protonów Z lub neutronów N, przy których zamyka się pewna grupa stanów jednocząstkowych: 2, 8, 20, 28, 50, 82, 126.¹ Ich znaczenie widać w wielu obserwablach naraz: w energiach wiązania, energiach pierwszych stanów wzbudzonych, przekrojach reakcji czy tendencji jąder do zachowywania kształtu bardziej kulistego i mniej podatnego na deformację. W praktyce „magiczne” oznacza tu tyle, że między ostatnim zajętym a najbliższym pustym poziomem istnieje wyraźna szczelina energetyczna, więc wzbudzenie układu kosztuje więcej energii niż w jądrze bez takiej szczeliny.¹

Tu przydaje się ostrożna analogia do atomów. W chemii szczególną stabilność mają gazy szlachetne, bo zamknięte powłoki elektronowe utrudniają dalsze wzbudzenia i reakcje. Królikowski pokazuje, że dla jąder atomowych działa podobna intuicja: niektóre jądra mają energię wiązania wyższą od sąsiadów właśnie dlatego, że dochodzi w nich do zamknięcia określonej grupy stanów protonowych albo neutronowych.³ Nie znaczy to oczywiście, że jądro jest po prostu „atomem w miniaturze”. Chodzi tylko o wspólną logikę dyskretnych stanów i szczególnej trwałości układów z zamkniętą powłoką.

Najprostsza wersja modelu powłokowego zakłada, że każdy nukleon porusza się w uśrednionym potencjale tworzonym przez wszystkie pozostałe. Rozwiązanie równania Schrödingera dla takiego potencjału daje zestaw stanów jednocząstkowych, które następnie zapełnia się zgodnie z zasadą Pauliego.¹ Sama ta idea nie wystarcza jednak do odtworzenia realnego układu liczb magicznych. Krytyczne znaczenie ma człon spin-orbit, czyli dodatkowe rozszczepienie poziomów zależne od sprzężenia spinu nukleonu z jego momentem orbitalnym.¹

Właśnie sprzężenie spin-orbit porządkuje poziomy tak silnie, że naturalnie pojawiają się dobrze znane zamknięcia powłok. Bez niego model jednocząstkowy byłby zbyt ubogi i nie odtwarzałby obserwowanej sekwencji liczb magicznych.^1,3 W sensie jakościowym działa to tak, że orbitale o tym samym l, ale różnym całkowitym momencie $j = l \pm \tfrac12$, przestają być prawie zdegenerowane. Jedne przesuwają się wyraźnie w dół energii, inne w górę, przez co w widmie poziomów otwierają się luki decydujące o „magiczności” określonych wartości N lub Z.¹

Na tym jednak nie kończy się fizyka jądra. Średni potencjał jest tylko pierwszym przybliżeniem, bo rzeczywiste nukleony nadal oddziałują ze sobą ponad to uśrednienie. Te poprawki nazywa się oddziaływaniami resztkowymi.¹ To one odpowiadają za to, że realne stany jądrowe są mieszankami konfiguracji, tworzą multiplety spinowe, mają określone przejścia gamma i nie redukują się do prostego obrazu „jedna cząstka na jednym orbitalu”. W języku praktycznym: model jednocząstkowy porządkuje mapę poziomów, a model wielocząstkowy z oddziaływaniami resztkowymi pozwala dopiero zbliżyć się do widm obserwowanych eksperymentalnie.¹

Ta uwaga jest ważna zwłaszcza wtedy, gdy mówi się o jądrach w pobliżu podwójnie magicznych rdzeni, takich jak 48Ca czy 208Pb. Fornal pokazuje, że właśnie tam szczególnie dobrze widać sens mówienia o nukleonach walencyjnych poza zamkniętym rdzeniem.¹ Część własności można wtedy interpretować jako wzbudzenia pojedynczych cząstek lub dziur, a część jako sprzężenie tych cząstek z drganiami i polaryzacją rdzenia. To już nie jest ani czysty model kroplowy, ani czysty model jednocząstkowy, lecz uporządkowana kombinacja obu perspektyw.

W materiałach Królikowskiego dobrze widać jeszcze jedną granicę zastosowalności. Model powłokowy jest szczególnie skuteczny dla jąder blisko zamkniętych powłok, natomiast dla jąder z większą liczbą nukleonów walencyjnych coraz silniej wchodzą do gry oddziaływania resztkowe i zjawiska kolektywne.³ Wtedy widma niskoenergetyczne zaczynają przypominać nie tyle „przeskoki jednej cząstki między orbitalami”, ile całe pasma rotacyjne albo oscylacyjne jądra zdeformowanego. Ten most rozwija osobny artykuł o wzbudzeniach jąder atomowych, przejściach gamma, pasmach rotacyjnych i oscylacyjnych.

Z punktu widzenia dydaktyki najważniejsze jest jednak coś jeszcze: układ powłok nie jest raz na zawsze dany dla wszystkich jąder. Jednym z głównych tematów nowoczesnej fizyki jądra stała się ewolucja struktury powłokowej wraz z oddalaniem się od ścieżki stabilności.¹ Innymi słowy, liczby magiczne dobrze znane dla jąder stabilnych nie muszą zachowywać swojej roli w jądrach bardzo neutrononadmiarowych lub protononadmiarowych. Zmienia się obsadzenie orbitali, a razem z nim zmieniają się średnie oddziaływania między protonami i neutronami.

W pracy Fornala wyraźnie widać tę logikę na przykładzie zamknięcia podpowłoki neutronowej przy N=32 w neutrononadmiarowych izotopach tytanu.¹ Nie chodzi tu o „obalenie” klasycznych liczb magicznych, tylko o bardziej subtelną rzecz: pojawienie się nowej szczeliny energetycznej w obszarze tablicy nuklidów, gdzie klasyczny, szkolny zestaw liczb magicznych nie dawałby jeszcze takiego wniosku. Eksperymentalna obserwacja wysoko położonego pierwszego stanu 2+ w 54Ti oraz zgodność struktury yrastowej z obliczeniami powłokowymi zostały zinterpretowane jako potwierdzenie, że przy N=32 rzeczywiście pojawia się dodatkowe zamknięcie podpowłokowe.¹

Mechanizm tego przesunięcia jest szczególnie ciekawy. Fornal wiąże je ze zmianą położenia neutronowego orbitalu $f_{5/2}$ względem $p_{3/2}$ i $p_{1/2}$ wtedy, gdy maleje obsadzenie protonowego orbitalu $f_{7/2}$.¹ Gdy odpowiednie oddziaływanie monopolowe między $\pi f_{7/2}$ i $\nu f_{5/2}$ słabnie, orbital $\nu f_{5/2}$ przesuwa się ku wyższym energiom. W połączeniu z dużym rozszczepieniem spin-orbit między $\nu p_{3/2}$ i $\nu p_{1/2}$ może to otworzyć szczelinę przy N=32 dla jąder o odpowiednio małym Z.¹ To dobry przykład tego, że model powłokowy nie jest tylko statyczną tabelką stanów, lecz narzędziem do śledzenia, jak sama mapa poziomów zależy od składu jądra.

To właśnie tutaj widać różnicę między pytaniem energetycznym a strukturalnym. Artykuł o energii wiązania jądra na nukleon odpowiada na pytanie, w którą stronę reakcje jądrowe są opłacalne energetycznie. Model powłokowy odpowiada na inne: jak zorganizowane są poziomy wewnątrz jądra i dlaczego niektóre wzbudzenia są łatwe, a inne trudne. Obie perspektywy są potrzebne równocześnie. Bez pierwszej nie zrozumiemy bilansu rozszczepienia i syntezy, a bez drugiej nie zrozumiemy, skąd biorą się szczegóły widm, magiczne liczby i zmienność struktury jąder egzotycznych.

W kontekście ciężkich jąder rozszczepialnych warto też postawić granicę zastosowalności. Model powłokowy nie zastępuje modelu kroplowego przy opisie samej deformacji prowadzącej do scission. Dobrze uzupełnia go jednak tam, gdzie trzeba wyjaśnić subtelniejsze efekty stabilizacji związane z zamknięciami powłok i strukturą stanów jednocząstkowych. W praktyce współczesna fizyka jądra korzysta z obu języków naraz: kroplowego dla ruchu kolektywnego i powłokowego dla mikroskopowej architektury poziomów.^1,3

Najkrótsze podsumowanie jest więc takie: model powłokowy porządkuje jądro od środka. Pokazuje, że protony i neutrony zajmują dyskretne orbitale, że sprzężenie spin-orbit buduje klasyczne liczby magiczne, a oddziaływania resztkowe i zmienne obsadzenie orbitali mogą daleko od stabilności tworzyć nowe podpowłoki. To nie jest konkurencja wobec modelu kroplowego, lecz brakująca druga połowa obrazu.¹

Historia: Maria Goeppert-Mayer i Hans Jensen

Choć idea powłokowa dla jąder atomowych pojawiała się już w latach 30., to dopiero w 1949 roku Maria Goeppert-Mayer (Niemcy/USA) i, niezależnie, Otto Haxel, Hans Jensen i Hans Suess (Niemcy) opublikowali kluczowy artykuł wyjaśniający, że sprzężenie spin-orbit jest niezbędne do odtworzenia obserwowanych liczb magicznych. Goeppert-Mayer pracowała nad problemem od 1948 roku i — jak sama wspomina — intuicja zasugerowała jej, że sprzężenie spin-orbit musi być duże, gdy rozmawiała z Enrico Fermim o osobliwościach jąder o magicznych Z i N.⁴

Praca Goeppert-Mayer była kamieniem milowym, bo pokazała, że proste modele niemające członu spin-orbit (model harmonicznego oscylatora, model studni) dawały liczby „pseudomagiczne" 2, 8, 20, 40, 70 — nie pasujące do obserwacji dla wyższych liczb. Dopiero sil­ne sprzężenie $\vec{l} \cdot \vec{s}$ pozwoliło naturalnie wydzielić 28, 50, 82, 126 jako kolejne zamknięcia powłok.

W 1963 roku Maria Goeppert-Mayer i Hans Jensen otrzymali Nagrodę Nobla z fizyki za odkrycie struktury powłokowej jąder atomowych, podzielając ją z Eugene'em Wignerem (za inne wkłady do fizyki jądrowej). Goeppert-Mayer była drugą kobietą — po Marii Curie — uhonorowaną tą nagrodą w dziedzinie fizyki.²

Historycznie warto zaznaczyć, że model powłokowy wywodził się z pracy Waltera Elsassera (1933/1934), który jako pierwszy zwrócił uwagę na szczególne trwałości pewnych obsadzeń. Jednak bez przekonującego mechanizmu fizycznego i obliczeń, jego praca nie weszła do głównego nurtu przez ponad dekadę. Dopiero Goeppert-Mayer i Jensen dostarczyli zarówno formalizmu, jak i fizycznego uzasadnienia przez analogię do atomowej struktury elektronowej, zastępując słabe sprzężenie LS (Russell-Saunders, typowe dla elektronów) silnym sprzężeniem jj (typowym dla nukleonów, bo oddziaływanie spin-orbit jest silniejsze).⁴

Potencjał jednocząstkowy: od oscylatora harmonicznego do Woods-Saxona

Punkt startowy modelu powłokowego to wybór potencjału, w którym nukleon się porusza. Historycznie rozważano kilka opcji:

Potencjał harmonicznego oscylatora: $V(r) = \frac{1}{2} m\omega^2 r^2$. Prosta analityczna forma, dająca dokładne rozwiązania. Zamknięcia powłok przy 2, 8, 20, 40, 70, 112 — nie pasują dla wyższych liczb magicznych.

Potencjał studni sferycznej z nieskończonymi ścianami: lekko lepszy, lecz nadal nie odtwarza 50, 82, 126.

Potencjał Woods-Saxona: stosowany współcześnie. Kształt sferoidy z łagodnie zanikającą powierzchnią:

$$V(r) = -\frac{V_0}{1 + e^{(r-R)/a}}$$

gdzie $V_0 \approx 50 \mathrm{MeV}$ — głębokość potencjału, $R = r_0 A^{1/3}$ — promień jądra ($r_0 \approx 1{,}2 \mathrm{fm}$), $a \approx 0{,}65 \mathrm{fm}$ — grubość powierzchni. Ta postać dobrze odwzorowuje gęstość materii jądrowej zmierzoną przez rozpraszanie elastyczne elektronów.⁴

Człon spin-orbit: do potencjału Woods-Saxona dodaje się człon:

$$V_{so}(r) = -\lambda \frac{1}{r}\frac{dV}{dr} \vec{l} \cdot \vec{s}$$

gdzie $\lambda \approx 20–25 \mathrm{MeV \cdot fm^2}$. Człon ten rozcina każdy orbital $l > 0$ na dwie sublinie: $j = l + \frac{1}{2}$ (niżej energetycznie) i $j = l - \frac{1}{2}$ (wyżej). Rozszczepianie rośnie z $l$, dlatego dla wysokich $l$ poziomy z $j = l + \frac{1}{2}$ są przesuwane głęboko — co tworzy obserwowane luki.⁴

Człon $l^2$: dodatkowy człon $-\kappa \hbar\omega l^2$ koryguje degenerację poziomów o tym samym $n$ ale różnym $l$ w oscylatorze harmonicznym. Łącznie potencjał Nilssona (WHO+$l^2$+$ls$) jest używany do opisania zdeformowanych jąder i obliczania diagram Nilssona.

Schemat poziomów jednocząstkowych

Poniżej schematyczny schemat (uproszczony) kolejności poziomów jednocząstkowych dla neutronów w obszarze $N = 0$–$126$:

Powłoka główna	Orbitale	$j$	Pojemność	Suma N
1s	1s$_{1/2}$	1/2	2	2 ← magiczna
1p	1p$_{3/2}$, 1p$_{1/2}$	3/2, 1/2	4, 2	8 ← magiczna
2s1d	1d$_{5/2}$, 2s$_{1/2}$, 1d$_{3/2}$	5/2, 1/2, 3/2	6, 2, 4	20 ← magiczna
1f2p	1f$_{7/2}$	7/2	8	28 ← magiczna
1f2p	2p$_{3/2}$, 1f$_{5/2}$, 2p$_{1/2}$	3/2, 5/2, 1/2	4, 6, 2	40 (podpowłoka)
2d1g	1g$_{9/2}$	9/2	10	50 ← magiczna
2d1g	2d$_{5/2}$, 1g$_{7/2}$, 3s$_{1/2}$, 2d$_{3/2}$	mix	6,8,2,4	82 ← magiczna
3p1h2f	1h$_{11/2}$, dalej...	11/2 +	12 +	126 ← magiczna

Luki energetyczne powyżej 2, 8, 20, 28, 50, 82, 126 nukleonów są bezpośrednimi przejawami zamknięcia powłok. Szczególnie duże luki przy 50 i 82 wynikają z przesunięcia w dół orbitali $1g_{9/2}$ i $1h_{11/2}$ przez silne sprzężenie spin-orbit przy wysokich $l$.

Jądra „podwójnie magiczne" i ich własności

Jądra, w których zarówno $Z$, jak i $N$ są magiczne, mają wyjątkowo wysoką energię wiązania i szczególne własności:

$^4$He (Z=2, N=2): najprostsze podwójnie magiczne. Cząstka alfa — wyjątkowo stabilna, emitowana przez ciężkie jądra w rozpadzie α.

$^{16}$O (Z=8, N=8): tlen-16 jest wyjątkowo trwały. Jego energia pierwszego stanu wzbudzenia ($J^\pi = 0^+_2$) leży aż przy 6,05 MeV — znacznie wyżej niż u sąsiadów. To czyni tlen-16 szczególnie ważnym benchmarkiem obliczeń ab initio.

$^{40}$Ca (Z=20, N=20): wapń-40 jest jedynym podwójnie magicznym jądrem złożonym wyłącznie z $N=Z$ (poza $^4$He i $^{16}$O). Bardzo wysoka energia pierwszego stanu $2^+$: 3,35 MeV.

$^{48}$Ca (Z=20, N=28): szczelnie zamknięte i N=28. Używany jako rdzeń dla obliczeń izotopów Ca z N=28–32. Energia pierwszego $2^+$: 3,83 MeV.

$^{132}$Sn (Z=50, N=82): podwójnie magiczny, ale daleko od stabilności (czas życia 39,7 s). Badany intensywnie przy ISOLDE (CERN) i HRIBF (Oak Ridge). Ważny dla zrozumienia struktury powłok aktynowców i procesu r.

$^{208}$Pb (Z=82, N=126): cięższy stabilny nuklid. Wyjątkowo wysoka energia pierwszego stanu $3^-$: 2,61 MeV. Używany jako tarcza w akceleratorach i jako materiał ochronny (gęsty, ale stabilny chemicznie). Ołów-208 jest benchmarkiem teorii jądrowej — każdy nowoczesny kod musi poprawnie odtwarzać jego widmo.⁴

Oddziaływania resztkowe i model wielocząstkowy (shell model diagonalization)

Model jednocząstkowy w średnim polu jest tylko punktem startowym. Oddziaływania resztkowe (residual interactions) między nukleonami walencyjnymi poza zamkniętym rdzeniem decydują o szczegółach widm energetycznych, momentach magnetycznych, momentach elektrycznych i prawdopodobieństwach przejść.

Pełne obliczenia w przestrzeni powłokowej (shell model diagonalization, SM) budują bazę Slatera ze wszystkich możliwych konfiguracji nukleonów walencyjnych na dostępnych orbitalach i diagonalizują macierz Hamiltonianu:

$$H = \sum_{i} \varepsilon_i \hat{n}_i + \sum_{i

gdzie $\varepsilon_i$ to energie jednocząstkowe (Single-Particle Energies, SPE), a $V_{ij}$ — efektywne oddziaływanie dwuciałowe (Two-Body Matrix Elements, TBME).

Główne programy obliczeniowe:

• ANTOINE (Caurier i in., IReS Strasburg): za pomocą metody Lanczosa diagonalizuje ogromne macierze, np. $10^9 \times 10^9$ dla pełnej przestrzeni $sd$-powłoki. Obsługuje przestrzenie $p$, $sd$, $pf$ i $fpg9$.
• NuShellX@MSU (Brown, Rae, MSU): interfejs użytkownika do OXBASH i powiązanych kodów.
• BIGSTICK (Johnson, USC): zoptymalizowany dla jąder z dużą liczbą cząstek.
• KSHELL (japoński kod): wydajny na wieloprocesorowych maszynach, używany w RIKEN i KEK.

Efektywne oddziaływania: USD (Universal $sd$ interaction, Wildenthal 1984), USDA/B (Brown i Richter 2006), USD+T (USD z potencjałem tensorowym), GXPF1 (Honma i in., dla $pf$-powłoki), KB3G, CCEI (coupled-cluster effective interaction) — coraz dokładniejsze przybliżenia ujednolicają opis wielu jąder w jednym regionie masy.⁴

Ewolucja struktury powłokowej — shell evolution

Jednym z najważniejszych odkryć fizyki jądra egzotycznego ostatnich dekad jest to, że liczby magiczne nie są stałe — ewoluują wraz z $N/Z$. Gdy odchodzimy od linii stabilności w kierunku jąder neutronobogatych lub protonadmiarowych, stosunek neutronów do protonów zmienia się drastycznie, co modyfikuje oddziaływania monopolowe między nukleonami walencyjnymi i przesuwa poziomy jednocząstkowe.

Przykłady znanych modyfikacji liczb magicznych:

• N=8 zanika w jądrach O-rich ($Z<8$): $^{11}$Li, $^{11}$Be, $^{12}$Be wykazują cechy niezamkniętej powłoki $N=8$. Be-12 ma stany o dużym momencie kwadrupolowym wskazujące na deformację.
• N=20 zanika w „wyspie odwrócenia" (island of inversion, $Z=10$–$12$, $N\approx20$): Ne, Na, Mg z $N=20$ wykazują duże deformacje zamiast kulistości charakterystycznej dla magicznych. Znany od końca lat 70.
• N=28 słabnie dla Ca-40 do Ca-36: badania energii stanu $2^+$ w $^{42}$Si wskazują na kolaps zamknięcia $N=28$ z dala od doliny stabilności.
• N=32 i N=34 jako nowe liczby magiczne dla neutronobogatego Ca: eksperymenty przy ISOLDE i RIKEN potwierdziły wysoko leżący stan $2^+$ w $^{52,54}$Ca — potwierdzenie nowych zamknięć podpowłokowych.

Mechanizmem tych zmian są oddziaływania monopolowe neutron-proton (tensor force i central force). Gdy liczba protonów zmienia się, zmienia się też oddziaływanie monopolowe z danym orbitalem neutronowym — przesuwając go względem sąsiadów. Ten efekt — nazywany tensor monopole shifting — jest aktywnym obszarem badań, zarówno eksperymentalnych, jak i ab initio z polem siłowym chiralnej efektywnej teorii pola.⁴

Jądra halo — ekstremalny przykład struktury powłokowej

Jądra przy neutronowej drip line mają tak słabo związane ostatnie neutrony, że ich funkcja falowa ma ogromne zasięgowe „ogony" — tworzące halo. Klasyczny przykład: lit-11 ($^{11}$Li, Z=3, N=8) ma dwuneutronowe halo: dwa ostatnie neutrony zajmują orbitale $1s_{1/2}$ i $1p_{1/2}$, których energia wiązania wynosi zaledwie ~300 keV. Promień materii $^{11}$Li jest porównywalny z $^{208}$Pb, mimo że ma tylko 11 nukleonów.

Jądra halo są laboratorium struktury powłokowej w ekstremalnych warunkach:

• Przekroje czynne na reakcje i pochłanianie są duże — wskazując na dużą sieć rozpraszania.
• Czynniki formfactor w reakcjach wytrącania (knockout) mierzą impulsowy rozkład neutronów halo.
• Momenty magnetyczne i elektryczne mierzą właściwości jąder walencyjnych bez zasłaniania przez rdzeń.

Technicznie jądra halo były dostępne dopiero po uruchomieniu instalacji radioaktywnych wiązek jonowych (Radioactive Ion Beams, RIB): ISOLDE (CERN/Genewa), TRIUMF/ISAC (Vancouver), HRIBF (Oak Ridge), RIKEN Nishina Center (Tokio), FRIB (Facility for Rare Isotope Beams, East Lansing, 2022). FRIB jest najbardziej wydajną instalacją na świecie — dostarcza izotopów ze stuprocentową ochroną cyklotronu 200 MeV/u z mocą wiązki 400 kW.⁴

Aplikacje modelu powłokowego poza fizyką akademicką

Model powłokowy ma konkretne zastosowania poza laboratorium:

Dane jądrowe dla reaktorów: przekroje czynne dla wychwytu neutronów i emisji gamma przez produkty rozszczepienia są wyznaczane przez pomiary przy intensywnych neutronowych źródłach (ILL Grenoble, HIFR ORNL) i interpretowane w kontekście modelu powłokowego. Model powłokowy pozwala interpolować i ekstrapolować przekroje dla izotopów, których nie można zmierzyć bezpośrednio.

Neutronowa aktywacja analistyczna (NAA): identyfikacja izotopów przez ich charakterystyczne widma γ wymaga tablic stanów wzbudzonych i branching ratios — opartych na obliczeniach modelu powłokowego.

Nukleosynteza: produkcja pierwiastków ciężkich (Ag, Ba, Pb, Au, Pt) przez procesy r i s jest wrażliwa na energie stanów jąder na granicy i poza doliną stabilności. Te dane — często niedostępne eksperymentalnie — muszą być dostarczone przez obliczenia powłokowe lub FRDM (Finite Range Droplet Model).⁴

Dawkimetria radiologiczna: obliczanie dawek od produktów rozszczepienia wymaga danych o widmach gamm emitowanych przy kaskadach rozpadu. Obliczenia modelu powłokowego uzupełniają tablice NNDC o dane dla jąder, których widma nie były bezpośrednio zmierzone.

Polskie wkłady do fizyki jądrowej

Polska szkoła fizyki jądrowej ma długą tradycję badań struktury jądra. Instytut Fizyki Jądrowej PAN (IFJ PAN) w Krakowie jest jednym z wiodących ośrodków eksperymentalnych i teoretycznych, z programami badawczymi przy akceleratorach europejskich (ISOLDE, GSI/FAIR, CERN/ALICE, JLab).

Bogdan Fornal (IFJ PAN), którego praca jest cytowana w artykule, pracował m.in. przy reaktorze ILL Grenoble i HIFR Oak Ridge, badając jądra neutrononadmiarowe przez reakcje β-rozpadu i głęboko nieelastycznego rozpraszania (GNIS). Jego prace o zamknięciu $N=32$ w izotopach Ti były częścią globalnego wysiłku eksperymentalnego, który zmienił rozumienie ewolucji struktury powłokowej w regionie $pf$-powłoki.

Reaktor badawczy MARIA w Świerku (NCBJ) dostarcza strumień neutronów $3{,}5 \times 10^{13} \mathrm{n \cdot cm^{-2} \cdot s^{-1}}$ i jest używany do pomiarów aktywacji, NAA, produkcji radioizotopów medycznych (Mo-99/Tc-99m) i badań struktury jądra przez rozpraszanie neutronów. Jest to jeden z ważniejszych reaktorów badawczych w Europie Środkowej.³

Korelacje między strukturą powłokową a własnościami jąder aktynowców

Aktynowce (U, Pu, Cm, Am) leżą daleko od zamkniętych powłok $Z=82$ ($^{208}$Pb) i $N=126$. Mają dużo nukleonów walencyjnych poza tymi rdzeniami, co powoduje:

• Silną deformację kwadrupolową ($\beta_2 \approx 0{,}3$–$0{,}4$)
• Bogatą strukturę rotacyjną i oscylacyjną
• Liczne izomery kształtu (deformed isomers)
• Silne kolektywne wzbudzenia GDR (Giant Dipole Resonance)

Jednak nawet dla aktynowców model powłokowy — w formie makro-mikro z poprawką Strutinskiego — dostarcza kluczowych informacji o barierze rozszczepienia (co omówiono w artykule o modelu kroplowym).

Szczególnie interesujące jest, że liczba neutronów $N=152$ odpowiada zamknięciu podpowłoki deformowanej w aktynowcach — jest to magiczna liczba dla silnie zdeformowanych jąder (deformed magic number). Jądra Cf-254 ($N=154$) i Fm-258 ($N=158$) leżą blisko i wykazują skrajnie krótkie czasy spontanicznego rozszczepienia — bezpośrednio związane z obniżoną barierą wynikającą ze struktury powłokowej przy tych N.¹

Metody eksperymentalne badania struktury powłokowej

Bezpośrednie badanie struktury powłokowej wymaga specjalizowanych technik eksperymentalnych. Najważniejsze:

Wzbudzenie Coulombowskie (Coulomb excitation, Coulex): bombardowanie ciężkiego jonu tarczą o dużym Z (np. złoto). Impuls elektromagnetyczny wzbudza niskie stany jądra projektila bez kontaktu jądrowego (przy odległościach $>1{,}2 \cdot A_P^{1/3} + 1{,}2 \cdot A_T^{1/3} \mathrm{fm}$). Pomiar prawdopodobieństwa przejść elektromagnetycznych (B(E2), B(M1)) daje bezpośredni wgląd w strukturę kolektywną i jednocząstkową. MINIBALL (ISOLDE), CHICO-2+GRETINA (NSCL), DALI2 (RIKEN) to główne detektory.

Reakcje transferu jednonukleonowego: dodawanie lub usuwanie jednego nukleonu przez reakcję typu ($d,p$), ($^3He,d$), ($p,d$), itd. pozwala mierzyć Spectroscopic Factors ($S$) — ile danej konfiguracji jednocząstkowej jest w stanie docelowym. Mierzony kąt i energia rozproszonych cząstek ($\theta$, $E$) wyznaczają moment orbitalny $l$ przenoszonego nukleonu. Technika DWBA (Distorted-Wave Born Approximation) lub CRC (Coupled Reaction Channels) przetwarza dane na parametry strukturalne.

Reakcje knockout (KO): przy energiach powyżej 100 MeV/A możliwe jest wytrącenie jednego nukleonu przez zderzenie z lekkim jądrem (C, Be). Spektrum pędów fragmentu po KO (np. $^{208}Pb(-1n)^{207}Pb$) daje rozkład pędu usuniętego nukleonu, a stąd jego orbital i spectroscopic factor. Reakcje KO mierzono przy GSI, NSCL, RIKEN.

Rozpad beta i reakcje $\gamma$: widma gamma po rozpadzie $\beta$ pozwalają na identyfikację stanów wzbudzonych przez kaskady przejść. Kluczowe instalacje: EURICA (RIKEN), GRIFFIN (TRIUMF), ISOLDE-ISOLTRAP.

Rozpraszanie elastyczne i nieelastyczne protonów: mierzy czynniki formfaktoru, potwierdza modele gęstości materii (neutronowej/protonowej skórki).⁴

Modele ab initio i chiralny potencjał jądrowy

W ostatnich 20 latach nastąpiła rewolucja w opisie teoretycznym struktury jądra: kody ab initio, obliczające własności jąder od pierwszych zasad z oddziaływania NN i 3N (trójciałowego).

Chiralny potencjał jądrowy (chiral EFT, Weinberg 1990, Entem-Machleidt 2003 i nowsze parametryzacje): oparty na chiralnej teorii perturbacji (ChPT), w której piony są mezonami mediatora oddziaływania NN. Oddziaływanie jest systematycznie rozwijane w potędze ($Q/\Lambda_\chi$)^n, gdzie $Q$ — typowy pęd nukleonu, $\Lambda_\chi \approx 1 \mathrm{GeV}$ — skala chiralnej symetrii. Kolejne rzędy (LO, NLO, NNLO, N3LO, N4LO) dodają coraz bardziej złożone diagramy i trójciałowe siły (3NF).

Kody ab initio:

• No-Core Shell Model (NCSM): rozwiązuje układ A nukleonów w bazie Slatera bez rdzenia. Dostępne dla $A \leq 20$. Kodach BIGSTICK, NTCM, MFDn.
• Coupled-Cluster (CC): rozwiązuje równania Singlesa-Doubles-Triples w bazie Slatera. Przesuwa barierę do $A \leq 100$. Kod CCSD(T) z ORCA lub NuShellX.
• In-Medium Similarity Renormalization Group (IM-SRG): transformuje Hamiltonian do quasi-diagonalnej postaci przez sekwencję unitarnych transformacji. Dostępny dla $A \leq 130$. Zaimplementowany w IMSRG++ (Vancouver).
• Self-Consistent Green's Function (SCGF): rozwiązuje dysonowe równania dla propagatora jednocząstkowego. Cenne dla gęstości stanów i funkcji Spectroscopic Factor.

Kody ab initio potwierdziły pojawienie się i zanik liczb magicznych w jądrach egzotycznych — np. pojawienie się $N=32$ i zanik $N=20$ w jądrach neonu neutrononadmiarowego — bez użycia fenomenologicznych dopasowań.⁴

Izomery kształtu i drugi minimim na PES

Ciekawym połączeniem modelu powłokowego i modelu kroplowego są izomery kształtu (shape isomers, fission isomers). Są to stany jądrowe leżące w drugim minimum energetycznym na Potential Energy Surface (PES) — przy deformacjach znacznie większych ($\beta_2 \approx 0{,}6$–$0{,}8$) niż deformacja stanu podstawowego ($\beta_2 \approx 0{,}2$–$0{,}3$).

Drugie minimum istnieje tylko w niektórych aktynowcach, gdzie efekty powłokowe przy dużej deformacji tworzą wyraźną szczelina energetyczną przy $N \approx 148$, $Z \approx 94$. Jądra w tym minimum są odizolowane barierą wewnętrzną od stanu podstawowego i barierą zewnętrzną od produktów rozszczepienia.

Przykłady:

• Am-242m: czas życia 14 ms w stanie izomerycznym; może rozszczepiać się spontanicznie.
• U-236 izomer: czas życia ~120 ns; znaleziony przez K. E. Bergkvis i in. w eksperymentach (n,γ).
• Pu-240 izomer: czas życia ~4 ns.

Izomery kształtu są mierzone przez ich charakterystyczne czasy spontanicznego rozszczepienia — drastycznie krótsze niż dla stanów podstawowych tego samego nuklidu. Na przykład Pu-240 w stanie podstawowym ma czas spontanicznego rozszczepienia $\sim 10^{11}$ lat, ale izomer Pu-240 — tylko nanosekundy. To bezpośredni efekt niższej bariery zewnętrznej dla zweiblendu minimum.^1,4

Liczby magiczne w projektowaniu bezpieczeństwa jądrowego

Choć brzmi to zaskakująco, liczby magiczne mają bezpośrednie zastosowanie w inżynierii bezpieczeństwa jądrowego. Kilka przykładów:

Trwałość pojemników do odpadów: Zr-90 ($N=50$), Sn-120 ($Z=50$) i Pb-208 ($Z=82, N=126$) są preferowanymi materiałami do pojemników na odpady jądrowe ze względu na ich wyjątkową odporność na napromieniowanie. Zamknięte powłoki oznaczają mniejsze prawdopodobieństwo pułapek dyslokacji i przebarwień radiacyjnych.

Mo-99/Tc-99m produkcja: Mo-99 jest produkowany przez rozszczepienie U-235 (jako jeden z produktów z wydajnością ~6,1%) lub przez wychwyt neutronu Mo-98. Tc-99m ($N=56$, nie magiczne, ale blisko $N=50$) ma wyjątkowe właściwości medyczne — czas połowicznego rozpadu 6,01 h i dominujące przejście 140 keV γ (bliskie optymum dla scyntygrafii). Zrozumienie widma produktów rozszczepienia i ich trwałości wymaga znajomości struktury powłokowej.

Reaktory ciężkowodne i Pb-Bi eutetyk: ołów jako chłodziwo-reflektor (reaktory LFR, BREST-OD-300 w Rosji) jest wybierany ze względu na małe pochłanianie neutronów i dużą gęstość. Jego skład izotopowy ($^{208}$Pb dominuje) i wynikające z tego małe przekroje czynne są bezpośrednią konsekwencją magiczności $N=126$.⁴

Kontekst dydaktyczny w Polsce i na świecie

Model powłokowy jest standardowym tematem kursów fizyki jądrowej na poziomie licencjatu (3. rok) i magisterium na polskich uczelniach. Podręczniki:

• K. Krane, Introductory Nuclear Physics (Wiley 1987): najpopularniejszy anglojęzyczny podręcznik na świecie, rozdział 5 poświęcony strukturze powłokowej. Dostępny elektronicznie.
• B. Burcham, M. Jobes, Nuclear and Particle Physics (Longman): bardziej zaawansowany.
• S. Hofmann, On Beyond Uranium (Taylor & Francis): o syntetycznych pierwiastkach i ich strukturze.
• Polskie podręczniki: Dostatni, Sołtan, Mayer-Kuckuk (tłumaczenie) — dostępne w bibliotekach AGH i Politechniki Warszawskiej.

Na poziomie doktoranckim kluczowe są kursy letnie: Jyväskylä Summer School (Finlandia, doroczny), TRIUMF Summer Institute (Kanada), GSI School on Condensed Matter (Niemcy). W Polsce: Szkoła Letnia Fizyki Jądrowej IFJ PAN (Kraków, co dwa lata).³

Perspektywy: wyspa stabilności i nowe zamknięcia powłok

Jednym z wielkich otwartych pytań fizyki jądrowej pozostaje: czy istnieje „wyspa stabilności" przy $Z \approx 114$–$126$, $N = 184$? Teoria modelu powłokowego i makro-mikro obliczenia przewidują dodatkowe luki energetyczne przy tych konfiguracjach — podobne do tych, które czynią Pb-208 tak trwałym. Jądra w tej wyspie miałyby czasy życia od sekund do nawet lat (gdyby luki były wystarczające duże), co pozwoliłoby na chemiczne badania pierwiastków Z=114–120.

Kluczowe pytanie brzmi: jak duże są te luki? Różne modele dają różne odpowiedzi:

• Makro-mikro z modelem Nilssona: luka przy $N=184$ wynosi ~2 MeV — wystarczająco duża dla wyraźnej wyspy.
• Relatywistyczne Mean-Field (RMF) z oddziaływaniem NL3*: luka przy $N=172$ i $N=184$ rzędu 1–2 MeV.
• Ab initio z chiralnym oddziaływaniem NN+3N: trudne do przeprowadzenia dla $A>300$, wyniki fragmentaryczne.

Eksperymentalne wskazówki: pierwiastki Fl (Z=114), Mc (Z=115), Lv (Z=116), Ts (Z=117), Og (Z=118) zostały syntetyzowane w JINR Dubna z czasami życia od milisekund do sekund. Og-294 ($N=176$) ma już czas życia ~0,89 ms — dłuższy niż oczekiwano bez premii powłokowej. To sugeru­je, że jesteśmy na brzegu wyspy, ale jej centrum ($N=184$) jest jeszcze poza zasięgiem obecnych akceleratorów.⁴

Program naukowy FAIR (Facility for Antiproton and Ion Research) w Darmstadt, którego budowa dobiega końca, ma dostarczyć wiązek ciężkich jonów wystarczających do syntezy jąder z $Z > 120$ i badania ich struktury. FRIB w East Lansing (USA) skupia się bardziej na jądrach neutronobogatych — lecz komplementarność obu instalacji otwiera możliwość zbadania zarówno wschodnich, jak i zachodnich krańców tablicy nuklidów.

Spójność modeli: od powłok do gwiazd neutronowych

Model powłokowy ma zadziwiający zasięg: te same orbitale i oddziaływania, które opisują strukturę Ca-48 i Pb-208, są ekstrapolowane do jąder neutralnych gwiazd neutronowych. Materia jądrowa przy gęstościach powyżej $\rho_0 = 0{,}16 \mathrm{fm}^{-3}$ nie jest opisana przez model powłokowy sensu stricto, ale ab initio obliczenia z chiralnym potencjałem NN+3N dostarczają równania stanu materii neutronowej (Energy of Neutron Matter, ENM), które jest warunkiem brzegowym dla modeli wnętrza gwiazd neutronowych.

Obserwacje fal grawitacyjnych ze zderzenia gwiazd neutronowych (GW170817, 2017) mierzą tzw. tidal deformability — parametr opisujący, jak łatwo gwiazda neutronowa deformuje się pod wpływem pola grawitacyjnego partnera. Ten parametr jest bezpośrednio związany z ENM i — przez to — z chiralnym potencjałem jądrowym. To znaczy, że pomiary LIGO/Virgo ograniczają parametry modelu jądrowego używanego w obliczeniach dla Ca-48 i Ni-78.

Fizyka jądrowa i astrofizyka zderzają się zatem w chiralnym potencjale jądrowym: jądra atomowe w laboratoriach ziemskich i gwiazdy neutronowe o masach ~$10^{57}$ nukleonów są opisane tą samą teorią. Model powłokowy jest jednym z węzłowych punktów tej teorii — i właśnie dlatego badania struktury jądra egzotycznego w FRIB lub RIKEN mają bezpośrednie konsekwencje dla astrofizyki obserwacyjnej.

Dodatkowe materiały multimedialne

Jako naturalne rozwinięcie warto po tym tekście wrócić do mechanizmu rozszczepienia w modelu kroplowym, energii wiązania jądra na nukleon oraz tekstów o konkretnych izotopach, takich jak uran-235 i uran-238. Wtedy widać już wyraźnie, że jedne artykuły opisują bilans energii i deformację ciężkich jąder, a ten tekst opisuje mikroskopowy układ ich stanów.

Ćwiczenia praktyczne

Ćwiczenie laboratoryjno-obliczeniowe powinno polegać na zbudowaniu prostego modelu poziomów jednocząstkowych z rozszczepieniem spin-orbit. W wariancie podstawowym należy:

1. przyjąć uporządkowaną listę orbitali dla wybranego obszaru powłokowego,
2. zaznaczyć pary stanów j = l + 1/2 i j = l - 1/2,
3. zasymulować ich rozdzielenie energetyczne po włączeniu członu spin-orbit,
4. zapełnić poziomy kolejno neutronami albo protonami zgodnie z zasadą Pauliego,
5. sprawdzić, przy których obsadzeniach pojawiają się największe szczeliny energetyczne.

Celem ćwiczenia nie jest odtworzenie pełnej fizyki jądra, lecz zobaczenie, że bez członu spin-orbit układ liczb magicznych wygląda inaczej niż w rzeczywistych danych. W drugiej części warto porównać jakościowo klasyczny zestaw zamknięć z przypadkiem, w którym jeden z orbitali przesuwa się wyżej lub niżej, aby zobaczyć, jak może powstać nowe zamknięcie podpowłokowe.

Drugie ćwiczenie, teoretyczno-przemysłowe, powinno uczyć pracy z rzeczywistymi danymi spektroskopowymi. Należy:

1. wybrać serię izotopów z jednego łańcucha, na przykład Ca albo Ti,
2. zestawić energie pierwszego stanu 2+,
3. zaznaczyć, gdzie pojawiają się skoki sugerujące większą szczelinę powłokową,
4. porównać te obserwacje z prostym schematem obsadzania orbitali,
5. opisać, które wnioski są bezpieczne, a które wymagałyby już pełnych obliczeń modelu wielocząstkowego.

To ćwiczenie ma nauczyć ostrożności interpretacyjnej. Pojedyncza energia wzbudzenia nie jest jeszcze całą teorią struktury jądra, ale bywa bardzo dobrym sygnałem, że układ poziomów zmienia się w sposób jakościowy i że prosta intuicja „stabilne jest to, co znaliśmy z jąder bliskich stabilności” przestaje wystarczać.

Przejdź do ćwiczenia interaktywnego

Powiązane artykuły

• Mechanizm rozszczepienia - model kroplowy
• Energia wiązania jądra na nukleon
• Defekt masy i równoważność energii
• Uran-235 - charakterystyka fizyczna