Streszczenie

Rozpad alfa jest jednym z najprostszych miejsc, w których mechanika kwantowa przestaje być abstrakcją. Cząstka alfa ma za małą energię, aby klasycznie przejść przez barierę Coulomba otaczającą ciężkie jądro, a jednak bywa emitowana. Wyjaśnienie wymaga tunelowania: niezerowej amplitudy przejścia przez obszar, który dla cząstki klasycznej byłby zabroniony.1

Ten artykuł prowadzi od prostego modelu bariery potencjału do rozpadu alfa, półokresu i widma alfa. Celem nie jest opis procedur pracy z emiterami alfa, lecz zrozumienie, dlaczego mała zmiana energii i wysokości bariery może oznaczać różnicę między rozpadem w mikrosekundach i rozpadem w czasie dłuższym niż wiek Ziemi.

Rozszerzenie tematu

Historia: Gamow, Gurney i Condon w 1928 roku

Kwantowa teoria rozpadu alfa jest jednym z pierwszych tryumfów mechaniki kwantowej w zastosowaniu do jąder atomowych. W 1928 roku, niezależnie od siebie, George Gamow (wówczas w Getyndze u Borna) oraz Ronald Gurney i Edward Condon (Princeton) opublikowali teorię tłumaczącą rozpad alfa przez tunelowanie.

Wcześniejszy problem: eksperymentalnie znano prawo Geigera-Nuttalla (1911) — empiryczną korelację między energią alfa a logarytmem stałej rozpadu. Nie było jednak fizycznego wyjaśnienia, dlaczego cząstka alfa może opuścić jądro z energią kilku MeV, skoro maksimum bariery Coulomba dla ciężkich jąder wynosi kilkadziesiąt MeV.

Gamow zastosował równanie Schrödingera do prostego modelu bariera Coulomba + studnia jądrowa i uzyskał wyrażenie:

T_1/2 ~ exp(π Z₁ Z₂ e² / ℏ v_α)

gdzie v_α jest prędkością cząstki alfa, a Z₁ Z₂ to produkt ładunków. To wyrażenie zwane jest czynnikiem Gamowa. Jego zależność od energii (przez v_α) odtwarza prawo Geigera-Nuttalla ze zrozumiałych zasad kwantowych, bez żadnych parametrów do dopasowania prócz masy i ładunku.

Odkrycie Gamowa stało się fundamentem fizyki jądrowej. Pokazało, że jądro można opisywać mechaniką kwantową. Gamow użył też tego samego mechanizmu do zbliżenia jąder w fuzji termojądrowej — i to było ziarno idei, które prowadziło do zrozumienia energetyki gwiazd (Bethe i Weizsäcker, 1938–1939).

W Polsce historia ta rezonuje ze względu na Georga Gamowa — Rosjanina, który w 1933 roku uciekł z ZSRR na Zachód na Kongresie Solvaya w Brukseli. Pracował potem w USA (George Washington University, następnie Boulder), gdzie był jednym z pierwszych popularyzatorów nauki (seria „Przygody pana Tompkinsa", polskie wydanie PWN). Część jego korespondencji naukowej przechowywana jest w Library of Congress. Gamow jest też znany z idei Wielkiego Wybuchu i teorii nukleosytezy pierwotnej (współpraca z Alphą i Hermanem, 1948).

Równolegle do Gamowa, tuż przed II wojną, Louis de Broglie i Niels Bohr rozwijali interpretacje mechaniki kwantowej. Teoria tunelowania Gamowa od razu stała się przykładem, że mechanika kwantowa nie jest tylko fizyką mikroskopową — ma mierzalne konsekwencje dla makroskopowych stałych fizycznych, takich jak wiek Ziemi (w połowie zbudowanej z ²³⁸U i ²³²Th, które rozpadają się tuż tak wolno, że trwają do dziś).

Bariera w fizyce klasycznej

W fizyce klasycznej bariera potencjału działa bezwzględnie. Jeśli cząstka o energii całkowitej E napotyka barierę o wysokości U i E < U, cząstka odbija się. Żeby znaleźć się po drugiej stronie, musiałaby w obszarze bariery mieć ujemną energię kinetyczną, co w klasycznej mechanice nie ma sensu.

Kuznetsov omawia najprostszy model: jednowymiarową barierę prostokątną o wysokości U i szerokości l.1 Klasycznie decyzja jest zero-jedynkowa: E > U oznacza przejście nad barierą, E < U oznacza odbicie. Mechanika kwantowa daje inną odpowiedź, bo opisuje cząstkę funkcją falową, a nie tylko punktem o ustalonej trajektorii.

Jeżeli E < U, funkcja falowa w obszarze bariery nie znika natychmiast. Maleje wykładniczo. Jeśli bariera jest skończonej szerokości, po drugiej stronie pozostaje mała, ale niezerowa amplituda. Kwadrat tej amplitudy daje niezerowe prawdopodobieństwo znalezienia cząstki za barierą.

Tunelowanie: nie przejście przez dziurę w ścianie

Słowo „tunelowanie” bywa mylące. Nie chodzi o to, że cząstka znajduje klasyczną szczelinę w barierze. Chodzi o to, że kwantowy opis ma niezerową amplitudę w obszarze klasycznie zabronionym. W wyniku pomiaru cząstka może pojawić się po drugiej stronie bariery, choć nie miała klasycznej energii potrzebnej do przejścia ponad nią.

Dla bariery prostokątnej prawdopodobieństwo przejścia maleje w przybliżeniu wykładniczo z szerokością bariery i z pierwiastkiem z różnicy U - E:

D ~ exp[-2 l sqrt(2m(U-E)) / hbar].

Dla bariery o dowolnym kształcie intuicja jest podobna:

D ~ exp[-(2/hbar) integral sqrt(2m(V(x)-E)) dx],

gdzie całka przebiega przez obszar klasycznie zabroniony. Ten zapis jest ważniejszy jakościowo niż rachunkowo: masa cząstki, szerokość bariery i różnica między barierą a energią wchodzą do wykładnika. Mała zmiana któregoś z tych parametrów może dać ogromną zmianę prawdopodobieństwa.

Dlaczego to pasuje do jąder atomowych

W jądrze działają dwie konkurencyjne intuicje. Siły jądrowe są silne, ale krótkiego zasięgu. Przyciągają nukleony na odległościach jądrowych. Odpychanie elektrostatyczne protonów jest słabsze na bardzo krótkiej skali, ale ma dłuższy zasięg.

Cząstka alfa, czyli jądro helu-4, może być traktowana jako zwarty klaster dwóch protonów i dwóch neutronów. W ciężkim jądrze emisja alfa może być energetycznie możliwa, jeśli masa jądra macierzystego jest większa niż suma mas jądra potomnego i cząstki alfa. W języku wartości Q oznacza to Q_alpha > 0.3

Sama dodatnia wartość Q nie wystarcza jednak do natychmiastowego rozpadu. Cząstka alfa musi wydostać się spod wpływu przyciągania jądrowego i przejść przez odpychającą barierę Coulomba. Klasycznie, jeśli jej energia jest mniejsza od maksimum tej bariery, nie mogłaby uciec. Kwantowo może tunelować.

Obraz bariery Coulomba

W uproszczonym obrazie potencjał dla cząstki alfa ma trzy obszary.

Wewnątrz jądra cząstka alfa jest związana przez oddziaływanie silne. W pobliżu powierzchni jądra pojawia się przejście między głębokim potencjałem jądrowym i odpychaniem Coulomba. Poza jądrem dominuje dodatni potencjał elektrostatyczny:

V(r) ~ (2 Z_d e^2) / (4 pi epsilon_0 r),

gdzie Z_d jest liczbą atomową jądra potomnego. Energia kinetyczna emitowanej cząstki alfa jest dodatnia, ale zwykle leży poniżej szczytu bariery. Obszar między punktem wewnętrznym i zewnętrznym, gdzie V(r) > Q_alpha, jest obszarem tunelowania.

Im wyższy ładunek jądra potomnego, tym silniejsze odpychanie Coulomba. Im większa energia alfa, tym węższy efektywny obszar zabroniony. Dlatego energie alfa i półokresy są tak silnie powiązane.

Rozpad alfa jako rozszczepienie skrajnie asymetryczne

Kuznetsov opisuje rozpad alfa jako skrajnie asymetryczne rozszczepienie: jądro macierzyste dzieli się na małą cząstkę alfa i duże jądro potomne.1 To ujęcie jest dydaktycznie wygodne, bo łączy rozpad alfa z ogólną ideą przemiany ciężkiego jądra na produkty silniej związane.

Różnica względem typowego rozszczepienia jest ogromna. W rozpadzie alfa produkt lekki jest zawsze jądrem helu-4, a proces jest zwykle pojedynczym kanałem rozpadu promieniotwórczego. W rozszczepieniu ciężkie jądro przechodzi przez deformacje, dzieli się na dwa fragmenty o rozkładzie mas i emituje neutrony. Oba zjawiska mogą jednak wymagać przejścia przez barierę potencjału, a w przypadku spontanicznych procesów bariera oznacza tunelowanie.4

To pokazuje wspólny język: energia może sprzyjać przemianie, ale bariera decyduje o szybkości.

Prawdopodobieństwo przejścia i półokres

Rozpad promieniotwórczy opisujemy stałą rozpadu lambda. Dla prostego modelu można myśleć o niej jako o iloczynie dwóch czynników:

lambda ~ f * P,

gdzie f jest częstością „prób” ucieczki klastra alfa z wnętrza jądra, a P prawdopodobieństwem przejścia przez barierę w pojedynczej próbie. To tylko model intuicyjny, ale dobrze pokazuje zależność.

Półokres wynosi:

T_1/2 = ln(2) / lambda.

Jeżeli P jest bardzo małe, półokres jest ogromny. Jeżeli bariera jest nieco niższa albo węższa, P może wzrosnąć o wiele rzędów wielkości, a półokres gwałtownie maleje. Właśnie dlatego emitery alfa mogą mieć półokresy od ułamków sekundy do miliardów lat.

Ta wykładnicza wrażliwość jest najważniejszą lekcją tunelowania. W chemii i mechanice klasycznej często spodziewamy się zmian proporcjonalnych. W tunelowaniu zmiana parametru w wykładniku potrafi całkowicie zmienić skalę czasu.

Prawo Geigera-Nuttalla jako ślad tunelowania

Historycznie zauważono, że emitery alfa o wyższej energii alfa mają zwykle krótsze półokresy. Zależność Geigera-Nuttalla wiąże logarytm półokresu z energią emitowanej cząstki alfa i ładunkiem jądra. Nie jest to przypadek empiryczny oderwany od teorii. Wynika z tunelowania przez barierę Coulomba: energia alfa zmienia szerokość obszaru klasycznie zabronionego.

Jeżeli energia alfa rośnie, zewnętrzny punkt zwrotny przesuwa się bliżej jądra. Bariera do przetunelowania staje się węższa. Prawdopodobieństwo przejścia rośnie wykładniczo, a półokres maleje. Dlatego różnice energii rzędu setek keV mogą odpowiadać ogromnym różnicom półokresu.

W dydaktyce warto to pokazywać na danych syntetycznych: nie trzeba pracować z materiałami promieniotwórczymi, aby zobaczyć, że log(T_1/2) jest znacznie bardziej naturalną osią niż samo T_1/2.

Co widzimy w spektrometrii alfa

Tunelowanie wyjaśnia, dlaczego cząstka alfa w ogóle opuszcza jądro i dlaczego czas rozpadu ma tak szeroki zakres. Spektrometria alfa mówi coś innego: jaka energia została zdeponowana w detektorze po emisji.

W artykule o widmie alfa opisano praktyczny problem: cząstki alfa bardzo łatwo tracą energię w powietrzu, warstwie preparatu i materiale detektora. Dlatego piki alfa mogą być precyzyjne, ale tylko przy dobrej geometrii, próżni i cienkim źródle.2

Połączenie obu perspektyw jest następujące:

  • tunelowanie ustala prawdopodobieństwo emisji i półokres,
  • wartość Q_alpha i poziomy jądra potomnego ustalają energie możliwych linii,
  • spektrometr i próbka decydują, jak te energie zostaną zarejestrowane.

Nie wolno mylić tych poziomów. Szeroki pik w detektorze nie oznacza, że samo tunelowanie było „rozmyte” energetycznie. Może oznaczać straty energii po emisji.

Rozpad alfa i poziomy jądra potomnego

Cząstka alfa nie zawsze zostawia jądro potomne w stanie podstawowym. Czasem przejście prowadzi do stanu wzbudzonego, który później może emitować gamma albo elektrony konwersji. Wtedy energia cząstki alfa jest mniejsza o energię wzbudzenia jądra potomnego.

Z punktu widzenia widma oznacza to, że jeden radionuklid może mieć kilka linii alfa. Z punktu widzenia tunelowania oznacza to, że różne kanały mają różne energie, a więc różne prawdopodobieństwa przejścia przez barierę. Z punktu widzenia danych jądrowych trzeba znać branching ratio, energie i niepewności.

Lokalna baza NKE ma w tym kontekście znaczenie dydaktyczne: pokazuje półokresy, tryby rozpadu, branching i energie linii, a także energie separacji takie jak S_alpha.5 Wartość separacji mówi o energetyce kanału, ale nie zastępuje obliczenia prawdopodobieństwa rozpadu. Energetycznie możliwe nie znaczy szybkie.

Tabela porównawcza emiterów alfa: Z, energia i T_1/2

Dane ilustrują prawo Geigera-Nuttalla — rosnąca energia alfa oznacza dramatycznie krótszy okres połowicznego zaniku, nawet przy podobnych wartościach Z.

Nuklid Z Q_alfa (MeV) E_alfa (MeV) T_1/2
²³²Th 90 4,082 4,013 14,05 mld lat
²³⁸U 92 4,270 4,197 4,47 mld lat
²³⁵U 92 4,679 4,397 704 mln lat
²²⁶Ra 88 4,871 4,784 1600 lat
²³⁹Pu 94 5,244 5,156 24 110 lat
²⁴¹Am 95 5,638 5,486 432 lata
²⁴⁴Cm 96 5,902 5,805 18,1 roku
²²²Rn 86 5,590 5,489 3,82 doby
²¹⁰Po 84 5,407 5,304 138,4 doby
²¹²Po 84 8,954 8,785 0,299 μs

Zakres dat: od 14 miliardów lat (²³²Th, zbliżone do wieku Wszechświata) do 0,299 mikrosekund (²¹²Po). To ponad 22 rzędy wielkości różnicy przy skromnej zmianie energii alfa i praktycznie tym samym Z. Czynnik Gamowa tłumaczy tę zmienność przez wykładniczą czułość współczynnika tunelowania — ΔE rzędu kilku MeV daje zmiany T_1/2 o dziesiątki rzędów.

Warto zwrócić uwagę na dwa izotopy polonu: ²¹⁰Po (T_1/2 = 138 d, znany jako broń użyta przeciw Aleksandrowi Litwienienko w 2006 r.) i ²¹²Po (sub-mikrosekunda). Oba emitują alfa, ale o zupełnie innych energiach i stąd tak odmiennych półokresach.

Szczegółowe widma alfa wybranych nuklidów

²¹⁰Po emituje praktycznie jednolinijkowe widmo: 5304 keV (99,999%). To izotop „modelowy" do kalibracji detektorów alfa. Brak drugorzędnych pików wynika z tego, że jądro potomne (²⁰⁶Pb) przechodzi niemal wyłącznie do stanu podstawowego. Stosowany jest w generatorach ciepła misji kosmicznych i historycznie jako inicjator neutronowy w dawnych konstrukcjach głowic.

²³⁹Pu ma bardziej złożone widmo: trzy główne linie przy 5105, 5144 i 5156 keV (razem ~99,9%). Proporcje tych linii są czułe na historię napromieniowania i wiek próbki — przez wkład ²⁴¹Am z rozpadu ²⁴¹Pu (T_1/2 = 14,3 roku). Spektrometria alfa ²³⁹Pu jest jedną z podstawowych metod analizy materiałów jądrowych w kontroli nieproliferacji.

²⁴¹Am emituje główną alfę 5486 keV (85%) i drugorzędną 5443 keV (13%). Gamma 59,5 keV towarzysząca rozpadowi umożliwia jednoczesną analizę spektrometryczną. ²⁴¹Am jest stosowany jako źródło wzorcowe detektorów alfa i jako materiał jonizacyjnych czujników dymu. Artykuł o jonizacyjnych czujnikach dymu opisuje, jak niewielka ilość ²⁴¹Am (ok. 1 μg) jest niezbędna do działania każdego takiego czujnika.

²²⁶Ra ma kilka linii: 4784 keV (94%), 4601 keV (6%). Historycznie kluczowe — Maria Curie i córka Irena pracowały z radem jako podstawowym wzorcem promieniotwórczości. Samo ²²⁶Ra rozpada się do ²²²Rn, gazu, który może opuścić preparat i pojawić się w powietrzu pomieszczenia.

²³²Th i ²³⁸U emitują alfa o niskiej energii (4–4,3 MeV), trudne do zmierzenia w próbkach środowiskowych ze względu na samopochłanianie. Ich oznaczanie wymaga zazwyczaj radiochemicznego wydzielenia. Artykuł o samopochłanianiu w spektrometrii alfa opisuje, jak głębokość preparatu rzędu kilku mikrogramów na centymetr kwadratowy istotnie zmienia kształt piku.

Czynnik Gamowa i przybliżenie WKB

Model kwantowy rozpadu alfa korzysta z przybliżenia WKB (Wentzel-Kramers-Brillouin), czyli quasi-klasycznego przybliżenia mechaniki kwantowej dla powoli zmieniającego się potencjału. Dla bariery Coulomba pomiędzy punktem wewnętrznym r₁ (powierzchnia jądra) a punktem zewnętrznym r₂ (gdzie V(r) = Q_alfa) całka WKB przyjmuje postać:

G = (2/ℏ) ∫_{r₁}^{r₂} sqrt(2μ(V(r)-Q_α)) dr

gdzie μ to masa zredukowana układu alfa + jądro potomne, a V(r) = Z_d · Z_α · e² / r (oddziaływanie Coulomba, Z_α = 2). Całka ta może być obliczona analitycznie dla potencjału Coulomba i daje:

G ≈ π Z_d e² sqrt(2μ) / (ℏ sqrt(Q_α)) · [arccos(sqrt(r₁/r₂)) - sqrt(r₁/r₂ · (1-r₁/r₂))]

W granicy r₁ ≪ r₂ (cienka bariera jądrowa w porównaniu z zewnętrznym punktem zwrotnym) wyraz w nawiasie kwadratowym zbliża się do π/2, co upraszcza wyrażenie do czynnika Gamowa:

G ≈ π Z_d e² sqrt(2μ) / (ℏ sqrt(Q_α)) = π Z_d α / v_α/c

gdzie α ≈ 1/137 to stała struktury subtelnej, a v_α/c to prędkość cząstki alfa względem prędkości światła. Dla ²¹⁰Po: Q_α = 5,4 MeV, Z_d = 82, G ≈ 24 – 27 w zależności od przyjętego r₁. Wtedy P ~ exp(-2G) ~ exp(-50) — niezwykle małe. Ale f (częstotliwość prób) wynosi ok. 10²¹ s⁻¹ (czas tranzytu przez jądro), więc λ = f·P ≈ 10²¹ · e^{-50} ≈ 10^{-2.7} s^{-1}, co daje T_1/2 ≈ 130 dni — wynik bliski obserwacji 138,4 d.

To jedno z pięknych potwierdzeń mechaniki kwantowej: model WKB bez parametrów empirycznych (poza promieniem jądra) daje właściwy rząd wielkości półokresu. Rozbieżność między ~130 d a ~138 d jest głównie wynikiem niedokładności modelu dla r₁.

Łańcuchy rozpadu i równowaga promieniotwórcza

Ciężkie nuklidy alfa-promieniotwórcze rzadko istnieją samodzielnie — wchodzą w łańcuchy szeregowe, gdzie córka jest znowu radioaktywna.

Szereg uranowy (4n+2):
²³⁸U → ²³⁴Th → ²³⁴Pa → ²³⁴U → ²³⁰Th → ²²⁶Ra → ²²²Rn → ²¹⁸Po → ²¹⁴Pb → ²¹⁴Bi → ²¹⁴Po → ²¹⁰Pb → ²¹⁰Bi → ²¹⁰Po → ²⁰⁶Pb

Łańcuch zawiera 8 rozpadów alfa i 6 beta, prowadząc do stabilnego ²⁰⁶Pb. Kluczowy element: ²²²Rn (radon) jest gazem szlachetnym — może opuścić skałę lub grunt i stać się problemem zdrowotnym w budynkach. Stężenia radonu w pomieszczeniach mierzy się właśnie przez produkty jego rozpadu (w szczególności ²¹⁴Pb i ²¹⁴Bi, emitujące gamma dobrze widoczne w spektrometrze NaI lub HPGe).

Szereg thorowy (4n):
²³²Th → ²²⁸Ra → ²²⁸Ac → ²²⁸Th → ²²⁴Ra → ²²⁰Rn → ²¹⁶Po → ²¹²Pb → ²¹²Bi → [²¹²Po + ²⁰⁸Tl] → ²⁰⁸Pb

Rozgałęzienie przy ²¹²Bi: 64% daje ²¹²Po (alfa 8,78 MeV, T_1/2 = 0,3 μs — krótki czas pozwala go identyfikować przez koincydencję z ²¹²Bi), 36% daje ²⁰⁸Tl (beta). Gamma 2614 keV z ²⁰⁸Tl jest jedną z najwyższych energii w naturalnym tle radiacyjnym i trudna do zasymulowania fałszywymi sygnałami.

Równowaga wiekowa: gdy czas od wyizolowania materiału jest długi w porównaniu z T_1/2 wszystkich krótkich ogniw, aktywność każdego członu łańcucha staje się równa aktywności ²³⁸U. Mierzysz jedno ogniwo — znasz wszystkie. To prawo Rutherforda (1902) potwierdzające, że aktywność jest addytywna wzdłuż łańcucha w równowadze.

Nieróżnoważność jest narzędziem analitycznym: w materiale geologicznym stosunek ²³⁰Th/²³⁴U odbiegający od jedności wskazuje na chemiczne frakcjonowanie (uran lepiej rozpuszcza się w wodzie niż tor). Metoda datowania U-Th (230 000–500 000 lat zasięgu) jest oparta właśnie na odchyleniach od równowagi.

Tunelowanie w chemii i biologii — pokrewne zjawiska

Efekt tunelowy nie ogranicza się do fizyki jądrowej. W chemii proton (masa ~1836 razy mniejsza niż nukleony alfa) ma dużo krótszą długość de Broglie'a niż cząstka alfa, ale przy bardzo małych barierach aktywacji zjawisko tunelowania protonów jest mierzalne eksperymentalnie.

Tunelowanie w enzymach: reakcje enzymatyczne przenoszące proton lub elektron (np. alkohole dehydrogenaza, dihydrofolian reduktaza) wykazują izotopowy efekt kinetyczny znacznie przekraczający wartości klasycznego modelu przejściowego. Stosunek stałych szybkości k_H/k_D (wodór/deuter) wynosi niekiedy ponad 10, podczas gdy klasyczna wartość Westheimer-Melander wynosi tylko ~7. To świadectwo tunelowania w warunkach biochemicznych.

STM i AFM: skaningowy mikroskop tunelowy (Scanning Tunneling Microscope, Binnig i Rohrer, Nobel 1986) działa właśnie na zasadzie tunelowania elektronu między ostrzem igły a powierzchnią próbki. Prąd tunelowy jest wykładniczo czuły na odległość (zmiana 1 Å → zmiana prądu o jeden rząd wielkości). STM pozwala obrazować pojedyncze atomy i manipulować nimi.

Dioda tunelowa Esakiego (Nobel 1973): w złączu p-n o bardzo silnym domieskowaniu pasma energetyczne nakładają się i elektrony mogą tunelować bezpośrednio przez wąską barierę złącza. Efekt daje ujemną rezystancję różniczkową (prąd maleje przy rosnącym napięciu w pewnym zakresie), co jest użyteczne w generatorach i wzmacniaczach mikrofalowych.

Efekt Josephsona: tunelowanie par Coopera (kondensatu nadprzewodnikowego) przez cienką warstwę izolatora między dwoma nadprzewodnikami. Wyróżnia się efekt DC (stały prąd bez napięcia) i AC (oscylacje prądu). Złącza Josephsona są podstawą najdokładniejszych woltomierzy kwantowych i interferometrów SQUIDowych, używanych m.in. do obrazowania aktywności mózgu.

Wszystkie te zastosowania łączy ten sam mechanizm matematyczny: wykładnicze zanikanie funkcji falowej w obszarze klasycznie zabronionym i niezerowa amplituda po drugiej stronie bariery.

Spektrometria alfa w monitoringu środowiskowym

W środowiskowych pomiarach radioaktywności spektrometria alfa odgrywa szczególną rolę dla plutonu, ameryku i nukleidów naturalnych (U, Th, Ra, Po). Gammy tych nuklidów są bardzo słabe lub ich nie ma, więc bezpośrednia spektrometria gamma jest nieskuteczna — potrzebna jest radiochemia i spektrometria alfa.

Procedura standardowa obejmuje cztery etapy:

  1. Mineralizacja próbki (kwasy, mikrofalowe trawienie) lub ługowanie mokre
  2. Wydzielenie radiochemiczne (ekstrakcja solwentową lub chromatografię jonowymienna) z dodaniem wzorca izotopowego (trac, spike) — np. ²³²U dla uranu, ²⁴³Am dla ameryku, ²⁰⁹Po dla polonu
  3. Osadzenie na dysku stalowym (elektroliza lub samodeponacja na Cu/Ag) — warstwa musi być ultracienka (< 0,1 mg/cm²)
  4. Pomiar detektorem krzemowym (pasywna implantacja lub PIPS) w próżni przez 12–72 godziny

Widmo alfa próbki środowiskowej nigdy nie jest czystą linią. Pojawia się ogon niskoenergetyczny od samopochłaniania i tło od radonu atmosferycznego. Efektywność procesu wydzielenia (wyliczana ze wzorca izotopowego) mieści się zwykle w 60–95%.

Polska sieć monitoringu: CLOR (Centralne Laboratorium Ochrony Radiologicznej) prowadzi oznaczenia alfa izotopów uranu i plutonu w glebie, wodzie i powietrzu. IFJ PAN (Kraków) specjalizuje się w oznaczeniach alfa izotopów Ra-226 i Po-210 w żywności i środowisku morskim (projekt Bałtyk). WAT (Wojskowa Akademia Techniczna) prowadzi badania metod detekcji alfa dla celów ochrony cywilnej.

Nuklearna kryminalistyka: alfa jako odcisk palca materiału jądrowego

Spektrometria alfa izotopów Pu i U jest jedną z podstawowych technik nuklearnej kryminalistyki (nuclear forensics). Skład izotopowy plutonu (stosunek ²³⁹Pu/²⁴⁰Pu/²⁴¹Pu) jest czułą funkcją historii napromieniowania — reaktora i czasu wypalenia. Płytki energetyczny i dobry recydyw izotopów (Pu-240 i Pu-239 są rozdzielane spektrometrem alfa tylko w referencyjnej geometrii próżniowej, bo ΔE ~ 11 keV) pozwalają na identyfikację. Pu-238 (T_1/2 = 87,7 lat) pochodzi głównie z absorpcji neutronów przez U-237 → Np-237 → Pu-238, a jego wysoka aktywność rozgrzewa próbkę — stąd termoelektryczne generatory (RTG) używają ²³⁸PuO₂.

W kontekście forensyki ważny jest też izotopowy wskaźnik starzenia ²⁴¹Pu/²⁴¹Am: izotop ²⁴¹Pu (beta, T_1/2 = 14,3 lat) przekształca się w ²⁴¹Am (alfa, T_1/2 = 432 lata). Stosunek ²⁴¹Am/²³⁹+²⁴⁰Pu w widmie alfa pozwala określić czas ostatniego oczyszczenia próbki. IAEA NSIL (Nuclear Security Incident Laboratory w Seibersdorfie) stosuje tę metodę standardowo do próbek przechwyconych materiałów.

Połączenie Polska: laboratoria alfa i projekty badawcze

W Polsce badania z zakresu spektrometrii alfa prowadzone są przede wszystkim w:

CLOR (Warszawa): pomiary alfa w ramach sieci wczesnego ostrzegania i badań środowiskowych; akredytowane metody oznaczania izotopów Pu, Am i Cm w próbkach środowiskowych.

IFJ PAN (Kraków): historycznie jedno z pierwszych polskich laboratoriów używających detektorów krzemowych do spektrometrii alfa (lata 70. XX w.), badania alfa w próbkach morskich i korelacja ze spowiem nuklearnym z lat 1945–1980.

NCBJ (Świerk): badania izotopowe Pu w ramach projektów safeguardowych, weryfikacja procedur NDA dla cezu i plutonu.

AGH (Kraków): metody radiochemiczne wydzielania Ra-226 i Po-210 z próbek środowiskowych, zastosowania w paleoklimatologii (korale) i oceanografii.

Wszystkie te laboratoria używają półprzewodnikowych detektorów PIPS (Partially Implanted Silicon, produkcja Mirion Technologies, wcześniej Canberra) lub ORTEC, pracy w próżni < 10⁻³ mbar i programów analizy widm (Alpha Analyst, GENIE 2000, dedykowane kody).

Spontaniczne rozszczepienie jako szersza lekcja

Tunelowanie przez barierę nie dotyczy tylko alfa. Bardzo ciężkie jądra mogą także spontanicznie rozszczepiać się przez tunelowanie przez barierę deformacji. Artykuł o mechanizmie rozszczepienia w modelu kroplowym opisuje to jako proces, w którym prawdopodobieństwo jest wykładniczo zależne od wysokości i szerokości bariery.4

W tym artykule wystarczy ogólna analogia. W rozpadzie alfa tuneluje stosunkowo dobrze zdefiniowany klaster. W rozszczepieniu spontanicznym tuneluje kolektywny stopień swobody związany z kształtem jądra. Matematyczny język jest podobny: istnieje obszar klasycznie zabroniony i niezerowe prawdopodobieństwo przejścia.

Ta analogia pomaga później zrozumieć, dlaczego same bilanse energii nie wystarczają do przewidywania trwałości jąder. Trzeba znać także bariery, poziomy kwantowe i prawdopodobieństwa przejść.

Tunelowanie a synteza jądrowa

Tunelowanie pojawia się również w reakcjach syntezy lekkich jąder. Dodatnio naładowane jądra odpychają się elektrostatycznie. Klasycznie musiałyby mieć bardzo dużą energię, aby zbliżyć się na odległość działania sił jądrowych. Kwantowo mogą z pewnym prawdopodobieństwem przenikać przez barierę Coulomba.

Kuznetsov pokazuje tę ideę na przykładzie zderzeń deuteronów: bariera elektrostatyczna jest wysoka w porównaniu z energiami termicznymi, a tunelowanie obniża wymaganą temperaturę względem prostego klasycznego oszacowania.1 W tym artykule nie rozwijamy zastosowań technicznych ani wojskowych. Dla nas ważny jest wspólny motyw: tunelowanie tłumaczy, jak reakcja może zachodzić mimo bariery, ale nie usuwa potrzeby spełnienia warunków energetycznych i statystycznych.

Mapa nukleidów i efekty powłokowe w rozpadzie alfa

Na mapie nukleidów (diagram Segré, oś N i Z) emitery alfa tworzą charakterystyczne obszary — niemal cała mapa dla Z > 82 jest alfa-promieniotwórcza lub wykazuje mieszany tryb alfa+beta. Kilka efektów struktury jądrowej modyfikuje proste przewidywania modelu tunelowego:

Magiczne liczby nukleonów: jądra o zamkniętych powłokach (N lub Z = 2, 8, 20, 28, 50, 82, 126) są silniej związane niż prognozuje gładka formuła. Jądra potomne ze strukturą blisko-magiczną mają wyższe energie wiązania, co zwiększa Q_alfa. Ale jądro macierzyste blisko-magiczne jest silniej związane niż sąsiedzi — co zmniejsza Q_alfa. Efekty te tworzą lokalne anomalie w krzywej Geigera-Nuttalla.

Nuklidy izomeryczne: wiele ciężkich nuklidów ma długo-żyjące stany wzbudzone (izomery), które mogą rozpadać się alfa niezależnie od stanu podstawowego, ale z inną energią Q i innym T_1/2. Przykładem jest ²⁴²mAm (izomer m, T_1/2 = 141 rok) obok ²⁴²Am (T_1/2 = 16 h), obydwa alfa+beta aktywne.

Efekty deformacji: jądra aktynowców są silnie zdeformowane (elipsoidalne). Cząstka alfa emitowana z jądra zdeformowanego może być emitowana z różną efektywnością wzdłuż osi symetrii i prostopadle do niej. Obserwuje się anizotropię kątową emisji alfa z spolaryzowanych jąder, badaną technikami ISOL (Isotope Separator On-Line).

Przenoszenie momentu pędu: emisja cząstki alfa może przenosić moment pędu l do jądra potomnego. Dla l > 0 bariera odśrodkowa dodaje się do bariery Coulomba (V_l(r) = ℏ²l(l+1)/(2μr²)), zmniejszając współczynnik tunelowania. Dlatego przejścia do stanów wzbudzonych jądra potomnego o wyższym spinie mają mniejszy branching niż proste prawa Geigera-Nuttalla sugerowałyby z samej energii Q.

Efekty te wyjaśniają, dlaczego dobre modele teoretyczne (np. kody Talys, FRESCO) wciąż nie odtwarzają branching ratio alfa z dokładnością lepszą niż czynnik 2–3 bez empirycznych parametrów struktury jądrowej.

Efektywna dydaktyka: łączenie poziomów abstrakcji

Nauczanie tunelowania alfa ma kilka pułapek dydaktycznych:

Pułapka 1: model prostokątny jako rzeczywistość. Bariera prostokątna jest narzędziem pedagogicznym, nie opisem fizycznym. Jądro ma potencjał Coulomba + potencjał jądrowy krótkiego zasięgu + centrifugal term + poprawki pairing. Uczący powinien wiedzieć, na którym poziomie uproszczenia operuje.

Pułapka 2: utożsamianie energii emisji z energią w detektorze. Cząstka alfa 5,4 MeV w próżni ma zasięg 3–4 cm w powietrzu, a po przejściu przez 1 μm złota traci ok. 100 keV. Widmo widziane przez detektor jest zawsze konwolucją energii emisji i strat. Osobno opisuje to artykuł o widmie alfa.

Pułapka 3: myślenie, że małe Q = niemożliwy rozpad. Wiemy już, że ²³²Th (Q = 4,08 MeV) ma T_1/2 = 14 mld lat. To nie jest „niemożliwy" rozpad — to tylko bardzo powolny. Jądra z Q < 0 dla alfa, np. ²⁰⁶Pb, naprawdę nie rozpadają się alfa; ale ²³²Th z Q > 0 rozpadnie się, zanim Słońce zgaśnie.

Pułapka 4: każde jądro ma tylko jedną linię alfa. W rzeczywistości nawet względnie proste nuklidy mają kilka linii. ²³⁸Pu ma trzy linie: 5456 (28%), 5499 (71%), 5546 (0,1%). ²⁴¹Am ma pięć linii (dominujące 5486 i 5443 keV). Im wyższe wzbudzenie dostępnych stanów jądra potomnego, tym więcej linii.

Pułapka 5: tunelowanie a spontaniczne rozszczepienie to „to samo". Mechanizm jest podobny (tunelowanie przez barierę), ale bariera rozszczepienia jest barierą względem kolektywnych stopni swobody kształtu, a nie barierą Coulomba dla ustalonego klastra. W rozszczepieniu spontanicznym nie ma preformowanej cząstki alfa — jądro deformuje się, a fragmenty tworzą się w procesie.

Dobre pytania sprawdzające: „Dlaczego ²¹²Po (T_1/2 = 0,3 μs) i ²¹⁰Po (T_1/2 = 138 d) mają tak różne T_1/2, skoro Z jest takie samo?" i „Co by się stało z T_1/2, gdyby cząstka alfa miała masę protonu?"

Dlaczego półokres jest statystyczny

Pojedyncze jądro nie ma zegarka odliczającego czas do rozpadu. Model kwantowy daje prawdopodobieństwo rozpadu w jednostce czasu. Jeśli mamy wiele identycznych jąder, liczba nierozpadniętych maleje wykładniczo:

N(t) = N_0 exp(-lambda t).

Półokres jest czasem, po którym średnio zostaje połowa jąder:

N(T_1/2) = N_0 / 2.

Tunelowanie daje mikroskopowe źródło tej losowości. Bariera nie jest pokonywana po zgromadzeniu energii jak w klasycznym modelu termicznym. Istnieje stałe prawdopodobieństwo przejścia związane z funkcją falową i stanem jądra. W dużej próbce daje to regularne prawo zaniku, choć pojedyncze zdarzenie jest losowe.

Czego nie wyjaśnia sam model prostokątny

Bariera prostokątna jest modelem szkolnym. Uczy wykładniczej zależności i niezerowej amplitudy za barierą. Jądro atomowe nie ma jednak prostokątnego potencjału. Potrzebny jest potencjał jądrowy krótkiego zasięgu, potencjał Coulomba, moment pędu, struktura powłokowa, możliwe stany potomne i prawdopodobieństwo utworzenia klastra alfa.

Dlatego dobry artykuł powinien jasno odróżniać trzy poziomy:

  • model prostokątny: pokazuje sam efekt tunelowy,
  • model bariery Coulomba: tłumaczy zależność energii alfa i półokresu,
  • dane jądrowe: podają rzeczywiste energie, branching i półokresy konkretnych nuklidów.

Jeżeli te poziomy się pomiesza, łatwo dojść do fałszywego wniosku, że jeden prosty wzór przewiduje wszystkie rozpady alfa z dużą dokładnością.

Praktyczne obliczenia z czynnikiem Gamowa: przykłady numeryczne

Poniżej dwie ilustracje numeryczne pokazujące, jak czynnik Gamowa łączy energię alfa z półokresem. Obydwa przykłady korzystają z uproszczonej formuły G ≈ π Z_d e² sqrt(2μ) / (ℏ sqrt(Q_α)), ignorując korektę na promień jądrowy.

Przykład 1: ²¹⁰Po → ²⁰⁶Pb + α

  • Q_α = 5,407 MeV = 5,407 × 1,602 × 10⁻¹³ J = 8,66 × 10⁻¹³ J
  • Z_d = 82 (ołów), μ ≈ 4 × 208 / (4 + 208) × 931,5 MeV/c² ≈ 3727 MeV/c²
  • G ≈ π × 82 × (1/137) × c / v_α, gdzie v_α/c = sqrt(2 Q_α / (m_α c²)) = sqrt(2 × 5,407/3727,4) ≈ 0,0539
  • G ≈ π × 82 / (137 × 0,0539) ≈ π × 82/7,38 ≈ 34,9
  • P ≈ exp(-2G) ≈ exp(-69,8) ≈ 2 × 10⁻³¹
  • f ≈ v_α / (2r_jądra) ≈ 0,054c / (2 × 7 × 10⁻¹⁵ m) ≈ 1,2 × 10²¹ s⁻¹
  • λ = f · P ≈ 1,2 × 10²¹ × 2 × 10⁻³¹ = 2,4 × 10⁻¹⁰ s⁻¹
  • T_1/2 = ln(2)/λ ≈ 2,9 × 10⁹ s ≈ 91 lat

Obserwacja: 138,4 doby. Różnica o dwa rzędy wynika z niedokładności przybliżenia r₁ (powierzchni jądra) i zaniedbania korekcji. Rząd wielkości jest jednak poprawny — to kluczowy wynik dydaktyczny.

Przykład 2: wpływ energii na T_1/2 przy stałym Z
Wyobraź sobie hipotetyczny emiter z Z_d = 82, ale Q_α zmienne:

  • Q_α = 5,0 MeV → G ≈ 36,3 → T_1/2 ≈ 10⁴ lat
  • Q_α = 5,5 MeV → G ≈ 34,5 → T_1/2 ≈ 30 lat
  • Q_α = 6,0 MeV → G ≈ 33,0 → T_1/2 ≈ 0,5 roku
  • Q_α = 7,0 MeV → G ≈ 30,5 → T_1/2 ≈ 2 doby
  • Q_α = 9,0 MeV → G ≈ 26,9 → T_1/2 ≈ 5 minut

Zmiana Q_α o czynnik 1,8 (od 5,0 do 9,0 MeV) daje zmianę T_1/2 o ponad 9 rzędów wielkości. To jest definicja wykładniczej wrażliwości — i jest to lekcja, której nie da się zastąpić żadnym przybliżeniem liniowym. Studentom warto narysować ten wykres i zapytać, w jakiej skali osi Y staje się on czytelny.

Bezpieczna granica dydaktyczna

Ten temat jest bezpieczny, gdy pozostaje na poziomie fizyki kwantowej, danych jądrowych i interpretacji widm. Nie wymaga opisu przygotowania źródeł alfa, separacji radiochemicznej ani pracy z materiałami toksycznymi lub kontrolowanymi.

W ćwiczeniach należy używać danych syntetycznych albo gotowych tablic rozpadu. Można liczyć wykładniczą zależność prawdopodobieństwa od szerokości bariery, porównywać półokresy na osi logarytmicznej i analizować gotowe widma alfa. Nie należy przekształcać tego w instrukcję pozyskiwania, oczyszczania lub preparatyki emiterów alfa.

Minimalna mapa pojęć

Na potrzeby dalszych artykułów warto zapamiętać:

  • klasycznie cząstka z E < U nie przechodzi przez barierę,
  • kwantowo funkcja falowa w barierze maleje, ale nie znika od razu,
  • prawdopodobieństwo tunelowania zależy wykładniczo od szerokości i wysokości bariery,
  • rozpad alfa jest ucieczką klastra alfa przez barierę Coulomba,
  • półokres wynika z prawdopodobieństwa przejścia w jednostce czasu,
  • energetycznie możliwy kanał rozpadu może być bardzo wolny.

Dodatkowe materiały multimedialne

Warto przygotować interaktywny model bariery tunelowej. Użytkownik zmienia masę cząstki, energię, wysokość i szerokość bariery, a wykres pokazuje funkcję falową, obszar klasycznie zabroniony i przybliżone log(P). Drugi tryb modelu może pokazywać barierę Coulomba dla rozpadu alfa bez przypisywania jej do pracy z konkretnym materiałem promieniotwórczym.

Powiązane kalkulatory i narzędzia

  • NKE - źródło danych o półokresach, trybach rozpadu, branching i energiach linii alfa.
  • Aktywność właściwa - pokazuje przejście od liczby atomów i masy nuklidu do aktywności.
  • ChainFinder - pomaga zobaczyć, że kanał rozpadu jest elementem sieci przemian, a nie izolowanym zdarzeniem.
  • Wizualizacja: Tunelowanie przez barierę — Animowany pakiet falowy pokazuje wykładniczą czułość prawdopodobieństwa tunelowania na masę, energię i szerokość bariery.
  • Kalkulator: Tunelowanie — Prosty współczynnik transmisji przez prostokątną barierę potencjału.

Ćwiczenia praktyczne

Pierwsze ćwiczenie: bariera prostokątna. Dla umownych wartości m, U-E i l policzyć względną zmianę D ~ exp(-a l) po podwojeniu szerokości bariery. Celem jest zobaczenie wykładniczej, a nie liniowej zależności.

Drugie ćwiczenie: półokres z prawdopodobieństwa. Przyjąć model lambda = fP z zadanym f i trzema wartościami P: 10^-10, 10^-20, 10^-30. Obliczyć T_1/2 = ln(2)/lambda i porównać skale czasu.

Trzecie ćwiczenie: energia alfa i bariera. Na danych syntetycznych narysować log(T_1/2) jako funkcję 1/sqrt(E_alpha). Student ma wyjaśnić, dlaczego zależność staje się prawie liniowa w takim układzie.

Czwarte ćwiczenie: widmo alfa a tunelowanie. Dla gotowego syntetycznego widma z dwiema liniami alfa wskazać, które informacje dotyczą procesu jądrowego, a które aparatury: energia linii, intensywność względna, FWHM, ogon niskoenergetyczny i tło.

Piąte ćwiczenie: energetyka nie wystarcza. Student dostaje trzy hipotetyczne kanały z dodatnim Q_alpha, ale różnymi barierami. Ma wskazać, dlaczego dodatnie Q nie przesądza o krótkim półokresie.

Przejdź do ćwiczenia interaktywnego

Powiązane artykuły