Streszczenie

Siły jądrowe to oddziaływanie utrzymujące nuklony razem w jądrze atomowym. Nie są to ani grawitacja, ani siły elektromagnetyczne — rządzą nimi inne prawa i inna skala energii. Trzy właściwości wyróżniają je wśród znanych oddziaływań: krótki zasięg (~1 fm), niezależność ładunkowa (takie same dla p–p, n–n i p–n) oraz nasycenie (nukleon oddziałuje tylko z kilkoma najbliższymi sąsiadami). Z tych właściwości wynikają bezpośrednio kształt krzywej energii wiązania, warunki stabilności jąder i granica, za którą odpychanie kulombowskie ostatecznie wygrywa z przyciąganiem jądrowym.1

Zasięg i głębokość potencjału

Siły jądrowe działają na odległościach rzędu jednego fermi ($1 \mathrm{fm} = 10^{-15} \mathrm{m}$). Poza odległością $r \approx (1{,}5\text{–}2{,}2) \times 10^{-15} \mathrm{m}$ oddziaływanie zanika praktycznie do zera. Tymczasem rozmiar atomu wynosi $\sim 10^{-10} \mathrm{m}$, czyli sto tysięcy razy więcej — elektrony nie „czują" sił jądrowych, co uzasadnia całkowite rozdzielenie fizyki atomowej i jądrowej.

Głębokość studni potencjalnej tworzonej przez siły jądrowe wynosi około $U \approx 50 \mathrm{MeV}$, podczas gdy energie wiązania atomów są rzędu eV. Ta różnica pięciu rzędów wielkości wyjaśnia, dlaczego reakcje jądrowe wyzwalają milion razy więcej energii niż reakcje chemiczne.1

Na bardzo małych odległościach ($r < 0{,}5 \mathrm{fm}$) siły jądrowe stają się odpychające — jest to tzw. „twarde jądro" (ang. hard core). Odpychanie na małych odległościach zapobiega zawaleniu się jądra i odpowiada za stosunkowo stałą gęstość jądrową ($\rho \approx 10^{17} \mathrm{kg/m^3}$, niezależną od $A$). Na pośrednich odległościach $0{,}5 \mathrm{fm} < r < 2 \mathrm{fm}$ dominuje silne przyciąganie.

Niezależność ładunkowa

Siły jądrowe wykazują niezależność ładunkową: przyciąganie między protonami ($p$–$p$), między neutronami ($n$–$n$) i między protonem a neutronem ($p$–$n$) jest takie same przy tym samym stanie spinowym i tej samej odległości. Dowód pochodzi z porównania energii wiązania jąder lustrzanych — par jąder o tej samej liczbie nuklonów $A$, ale z zamienionymi rolami protonów i neutronów, np. ${}^3_2\mathrm{He}$ (2 protony, 1 neutron) i ${}^3_1\mathrm{T}$ (1 proton, 2 neutrony).

Energie wiązania tych jąder różnią się o około 0,77 MeV. Różnicę tę można wyjaśnić wyłącznie odpychaniem kulombowskim protonów w ${}^3_2\mathrm{He}$:

$$\Delta E_k = \frac{e^2}{4\pi\varepsilon_0 r_{pp}} \approx 0{,}77 \mathrm{MeV}$$

co daje odległość $r_{pp} \approx 1{,}9 \mathrm{fm}$, zgodną z rozmiarem jądra helu. Skoro energię wiązania tłumaczy samo kulombowskie poprawka, wynika z tego, że czyste siły jądrowe były dla obu jąder identyczne — niezależność ładunkowa jest potwierdzona.1

Formalnym przejawem niezależności ładunkowej jest pojęcie izospinu: proton i neutron to dwa stany tworzącej je cząstki — nukleonu — różniące się wartością izospinu ($T_z = +1/2$ dla protonu, $T_z = -1/2$ dla neutronu). Oddziaływanie jądrowe jest niezmiennicze ze względu na obroty w przestrzeni izospinowej.

Nasycenie i proporcjonalność energii wiązania do A

Siły jądrowe wykazują nasycenie: każdy nukleon oddziałuje tylko z ograniczoną liczbą najbliższych sąsiadów. Gdyby tak nie było — gdyby każdy nukleon oddziałował z każdym innym — całkowita energia wiązania rosłaby jak $A(A-1)/2 \approx A^2$, a energia wiązania na nukleon rosłaby proporcjonalnie do $A$. Obserwacje pokazują jednak, że energia wiązania na nukleon jest prawie stała ($\approx 8 \mathrm{MeV/nukleon}$) dla jąder pośrednich, co oznacza liniową zależność całkowitej energii wiązania od $A$.

Nasycenie wynika ze spinowo-izospinowej struktury sił jądrowych. Jeżeli dwa nuklony tworzą parę w stanie singletowym (antyparallelne spiny), energia wiązania jest mniejsza niż w stanie tripletowym (równoległe spiny). Wynika stąd, że układ ${}^4_2\mathrm{He}$ (dwa protony + dwa neutrony, wszystkie sparowane) jest wyjątkowo stabilny — to tzw. cząstka alfa, która pojawia się jako autonomiczna jednostka zarówno w rozpadzie α, jak i w strukturze powłokowej ciężkich jąder.

Zależność od spinu

Siły jądrowe zależą od orientacji spinów oddziałujących nuklonów. Deuteron (${}^2_1\mathrm{H}$, jądro deuteru: proton + neutron) istnieje tylko w stanie spinu całkowitego $I=1$ (spiny równoległe). Stan antyrównoległy ($I=0$) nie tworzy związanego układu. To pokazuje, że oddziaływanie p–n w stanie tripletowym jest silniejsze niż w stanie singletowym.1

Doświadczalne potwierdzenie pochodzi z różnicowego rozpraszania neutronów na orto- i parawodorze. W ortowodorze spiny protonów obu atomów są równoległe, w parawodorze — antyparallelne. Przekrój czynny na rozpraszanie neutronów przez parawodór jest 30 razy większy niż przez ortowodór, co wskazuje na zasadniczo różną siłę oddziaływania n–p w zależności od orientacji spinów.1

Siły jądrowe mają ponadto charakter tensorowy (niecentralny): potencjał oddziaływania zależy nie tylko od odległości między nuklonami, ale też od kąta między wektorem łączącym nuklony a kierunkami ich spinów. Właśnie tensorowy charakter sił powoduje, że deuteron nie ma idealnie kulistego kształtu — jego elektryczny moment kwadrupolowy jest różny od zera.

Model wymienny — mezony π i pole jądrowe

Pierwszą kwantową teorię sił jądrowych zaproponował Hideki Yukawa w 1935 roku. Analogicznie do pola elektromagnetycznego, przenoszonego przez bezmasowe fotony, siły jądrowe mają być przenoszone przez mezony π (piony) o masie $m_\pi \approx 135 \mathrm{MeV/c^2}$. Zasięg oddziaływania wiąże się z masą nośnika:

$$R \sim \frac{\hbar}{m_\pi c} \approx 1{,}4 \mathrm{fm}$$

Potencjał Yukawa dla oddziaływania wymiennego ma postać:

$$U(r) = -g_\pi \frac{e^{-r/R}}{r}$$

gdzie $g_\pi$ jest analogiem ładunku elektrycznego dla oddziaływania jądrowego. Eksperymentalne odkrycie pionów naładowanych ($\pi^\pm$) przez Lattes, Occhiallini i Powella w 1947 roku oraz pionu neutralnego ($\pi^0$) w 1950 roku potwierdziło model Yukawa. Przy wymianie jednego pionu odtwarzana jest właśnie ta część potencjału jądrowego, która odpowiada zasięgom $r \gtrsim 1 \mathrm{fm}$. Krótkozasięgowe i skomplikowane oddziaływania przy małych $r$ wymagają uwzględnienia cięższych mezonów (ρ, ω) i ostatecznie opisu w językach chromodynamiki kwantowej (QCD).3

Granica stabilności i dominacja kulombowska

Siły jądrowe są silne, ale krótkozasięgowe. Odpychanie kulombowskie między protonami, choć słabsze, ma zasięg nieograniczony. Dla małych jąder jądrowe przyciąganie wyraźnie dominuje, ale wraz z wzrostem $Z$ całkowita siła kulombowska rośnie jak $Z(Z-1) \propto Z^2$, podczas gdy energia wiązania rośnie tylko jak $A$. Istnieje więc $Z$ graniczne, po przekroczeniu którego żadna konfiguracja neutronów nie wystarcza do utrzymania jądra — stąd koniec tablicy nuklidów naturalnych (~$Z = 92$). Dalsze pierwiastki powstają w reaktorach i akceleratorach i mają coraz krótsze okresy połowicznego zaniku.1

Zależność ta jest wprost zakodowana w formule Weizsäckera (omówionej w artykule o mechanizmie rozszczepienia): człon kulombowski $\propto Z^2 A^{-1/3}$ jest przyczyną, dla której ciężkie jądra leżą na opadającym ramieniu krzywej energii wiązania i mogą zyskać energię przez rozszczepienie.

Rozszerzenie tematu

Fenomenologiczne potencjały jądrowe

Zamiast od podstaw wywodzić oddziaływanie nukleon-nukleon z QCD (co jest nieperturbacyjnie trudne), fizycy konstruują fenomenologiczne potencjały — wzory matematyczne dopasowane do tysięcy danych rozpraszania n-p, p-p i do własności lekkich jąder.

Ogólna struktura potencjału NN jest złożona z wielu członów:

$$V_{NN}(\vec{r}, \vec{\sigma}_1, \vec{\sigma}_2, \vec{\tau}_1, \vec{\tau}_2) = V_C(r) + V_\sigma(r)\, \vec{\sigma}_1\cdot\vec{\sigma}_2 + V_\tau(r)\, \vec{\tau}_1\cdot\vec{\tau}_2 + V_{\sigma\tau}(r)\, (\vec{\sigma}_1\cdot\vec{\sigma}_2)(\vec{\tau}_1\cdot\vec{\tau}_2) + V_T(r)\, S_{12} + \ldots$$

gdzie $S_{12} = 3(\vec{\sigma}_1\cdot\hat{r})(\vec{\sigma}_2\cdot\hat{r}) - \vec{\sigma}_1\cdot\vec{\sigma}_2$ to operator tensorowy, $\vec{\sigma}$ i $\vec{\tau}$ — operatory spinu i izospinu.

Główne potencjały w użyciu:

Potencjał Reida (1968): Jeden z pierwszych precyzyjnych potencjałów fenomenologicznych, używający sum eksponensjalnych dopasowanych oddzielnie w każdym kanale $LS$ (partial wave). Dobrze odtwarza fazy rozpraszania do $E_{lab} \approx 300 \mathrm{MeV}$ i własności deuteronu.

Potencjał Paryski (1980): Oparty na teorii wymiany mezonów (OBE — One Boson Exchange), uzupełnionej przez wymianę dwóch pionów. Lepszy opis przy krótkich odległościach.

Argonne AV18 (1995): Wysoko-precyzyjny potencjał z 18 operatorowymi strukturami ($V_{18}$). Dopasowany do 4301 danych n-n i p-p i 2025 danych p-p $poniżej 350 \mathrm{MeV}$, z $\chi^2/\mathrm{datum} \approx 1$. Zawiera łamanie izospinu ($p$-$p \neq n$-$n$). Jest standardem dla obliczeń ab initio jąder lekkich.

CD-Bonn i Nijmegen potentials: Alternatywy do AV18 z nieco innym formalizmem, zbliżoną jakością dopasowania.

Potencjał sił jądrowych (potencjał Reida) w kanale ¹S₀: odpychanie dla r<0,5 fm, głęboka studnia przyciągania przy r≈0,8 fm, zanikanie do zera przy r>2 fm. CC BY-SA 4.0 Bdushaw 2015.
Potencjał nukleon-nukleon w kanale ¹S₀ (Argonne potencjał Reida 1968): wyraźnie widoczne twarde jądro odpychające dla r < 0,5 fm, minimum potencjału ~−50 MeV przy r ≈ 0,8 fm i zanikanie poniżej 0 dla r > 2 fm. CC BY-SA 4.0.

Deuteron jako laboratorium sił jądrowych

Deuteron (${}^2_1\mathrm{H}$, proton + neutron, $I = 1^+$, $J = 1$) jest jedyną związaną parą dwu nukleonów i przez to jest idealnemu laboratorium do badania oddziaływania p-n.

Podstawowe własności deuteronu:

  • Masa: $m_d = 1875{,}613 \mathrm{MeV}/c^2$
  • Energia wiązania: $B_d = 2{,}2246 \mathrm{MeV}$ (względnie mała!)
  • Spin: $I = 1$ (stany tripletowe ${}^3S_1$, ${}^3D_1$)
  • Elektryczny moment kwadrupolowy: $Q_d = +2{,}86 \times 10^{-31} \mathrm{m^2} = 0{,}286 \mathrm{e\cdot fm^2}$
  • Magnetyczny moment: $\mu_d = 0{,}8574 \mu_N$ (prawie ale nie dokładnie $\mu_p + \mu_n = 0{,}879 \mu_N$)

Niezerowy $Q_d > 0$ (kształt prolate — jajko): świadczy o tym, że funkcja falowa deuteronu nie jest czystym ${}^3S_1$, lecz superpozycją $^3S_1$ i $^3D_1$:

$$|\psi_d\rangle = u(r) {}^3S_1 + w(r) {}^3D_1$$

Komponent D-wave (orbitalny moment pędu $l = 2$) jest wynikiem tensorowego członu potencjału jądrowego. Jego udział szacowany jest na $P_D \approx 5\text{–}6\%$.

Odchylenie $\mu_d \neq \mu_p + \mu_n$ (o ~2,5%) jest też wynikiem D-wave i relatywistycznych poprawek (prądy wymiany mezonów, meson exchange currents). Takie poprawki są mierzone w eksperymentach elektrodezintegracji deuteronu: $d(e,e'n)p$ i $d(\gamma,pn)$.

Singletowy stan p-n (${}^1S_0$, $I = 0$) jest niemalże związany: energia własna to $E = -60 \mathrm{keV}$ poniżej progu (wirtualny stan). To "prawie"-związanie jest przyczyną wyjątkowo dużej długości rozpraszania w tym kanale: $a_s = -23{,}75 \mathrm{fm}$ (dla porównania $a_t = +5{,}42 \mathrm{fm}$ w stanie tripletowym). Ujemna długość rozpraszania oznacza właśnie wirtualny stan.

Parametryzacja przy niskich energiach — długość rozpraszania i efektywny zasięg

Przy energiach neutronów $E \ll E_R$ (poniżej pierwszego rezonansu) całkowity przekrój czynny jest opisany tylko dwoma parametrami: długością rozpraszania $a$ i efektywnym zasięgiem $r_0$.

Rozwinięcie:

$$k \cot\delta_0 = -\frac{1}{a} + \frac{1}{2}r_0 k^2 + \ldots$$

gdzie $k = \sqrt{2m_n E}/\hbar$ i $\delta_0$ — faza rozpraszania w stanie $s$.

W granicy energii zerowej: $\sigma_{tot} = 4\pi a^2$ (dla rozpraszania na sferze twardej wynosi to $4\pi R^2$; siły jądrowe dają $a \neq R$ ze względu na potencjał przyciągający).

Dla n-p (sumując po stanach spinowych):

$$\sigma_{total}(E \to 0) = \frac{3}{4}\sigma_t + \frac{1}{4}\sigma_s = 3\pi a_t^2 + \pi a_s^2 \approx 0{,}282 + 17{,}7 = 18{,}0 \mathrm{b}$$

To doskonale zgadza się z mierzonym przekrojem czynnym na rozpraszanie neutronów termicznych na protokach (wodorze) $\sigma_{np} = 20{,}4 \mathrm{b}$. (Różnica wynika z koherencji kwantowej i efektów powiązanych.)

Tunelowanie przez barierę Coulomba — teoria Gamowa

Siły jądrowe przyciągają, ale bariera Coulomba odpycha. Cząstka alfa wewnątrz jądra ciężkiego jest w studni potencjalnej: wewnątrz $V_{jad} \approx -30 \mathrm{MeV}$, na zewnątrz $V_C(r) = Z_1 Z_2 e^2/(4\pi\varepsilon_0 r) > 0$. Bariera Coulomba wynosi przy powierzchni jądra:

$$V_C(R) = \frac{Z_\alpha \cdot Z_{córka} \cdot e^2}{4\pi\varepsilon_0 R} \approx \frac{2 \cdot (Z-2) \cdot 1{,}44 \mathrm{MeV\cdot fm}}{R}$$

Dla ${}^{238}$U ($Z=92$, $R \approx 7{,}4 \mathrm{fm}$): $V_C(R) \approx 2 \times 90 \times 1{,}44 / 7{,}4 \approx 35 \mathrm{MeV}$, podczas gdy energia alfa-cząstki $E_\alpha = 4{,}27 \mathrm{MeV}$. Klasycznie: niemożliwe. Kwantowo: tunelowanie.

Gamow (1928) obliczył prawdopodobieństwo tunnelowania przez barierę Coulomba:

$$P \propto \exp\left(-2G\right), \qquad G = \frac{\sqrt{2m_\alpha}}{\hbar}\int_R^{R_{turn}} \sqrt{V_C(r) - E_\alpha}\, dr$$

gdzie $R_{turn} = Z_\alpha Z_{córka} e^2 / (4\pi\varepsilon_0 E_\alpha)$ to punkt zwrotny klasycznego odwrotu. Wynik całkowania (po aproksymacji):

$$G \approx \pi Z_\alpha Z_{córka} e^2/(4\pi\varepsilon_0 \hbar v_\alpha) \propto Z_{córka}/\sqrt{E_\alpha}$$

Stąd wynika prawo Geigera-Nuttalla (odkryte empirycznie 1911, wyjaśnione przez Gamowa 1928):

$$\log_{10} T_{1/2} = A \cdot Z / \sqrt{E_\alpha} + B$$

Praktyczna konsekwencja: niewielka zmiana $E_\alpha$ powoduje ogromną zmianę $T_{1/2}$. Przykłady:

Nuklid $E_\alpha$ (MeV) $T_{1/2}$
${}^{232}$Th 4,01 $1{,}4 \times 10^{10}$ lat
${}^{238}$U 4,27 $4{,}47 \times 10^9$ lat
${}^{226}$Ra 4,87 1600 lat
${}^{210}$Po 5,30 138 dni
${}^{214}$Po 7,69 163 μs

Zmiana $E_\alpha$ od 4 MeV do 7,7 MeV (czynnik 1,9) zmienia $T_{1/2}$ o czynnik $10^{23}$ — 23 rzędy wielkości! To niezwykły dynamiczny zakres, który mechanika kwantowa wyjaśnia naturalnie, a klasyczna fizyka nie.

Ten sam mechanizm tunelowania przez barierę Coulomba jest podstawą reakcji syntezy termojądrowej — tam nuklony muszą się do siebie zbliżyć pokonując barierę Coulomba przez tunelowanie lub energię termiczną ($kT$). Szczegóły w artykule o kontrolowanej syntezie termojądrowej.

Optyczny model oddziaływania jądra z cząstką

W zderzeniach cząstki (protonu, neutronu, jonu alfa) z jądrem, jądro zachowuje się jak "absorber" neutronów i protonów, podobny do matowej szyby dla fotonów. Model optyczny zastępuje pełny wielociałowy problem jedną cząstką poruszającą się w zespolonym potencjale optycznym:

$$U(r) = V(r) + i W(r)$$

Część rzeczywista $V(r) < 0$: potencjał przyciągający (studnia Woodsa-Saxona), reprezentujący oddziaływanie z nukleoniami jądra.

Część urojona $W(r) < 0$: pochłanianie — cząstka trafia w kanały reakcji niesprężystej (reakcje jądrowe, produkcja pre-równowagowych cząstek), co skutecznie usuwa ją ze stanu sprężystego. Pozorne "wygaszanie" amplitudy fali — stąd analogia z optycznym współczynnikiem absorpcji.

Model optyczny jest standardowym narzędziem do:

  • Obliczania elastycznych przekrojów różniczkowych $d\sigma/d\Omega$
  • Szacowania całkowitych przekrojów reakcji (przez imaginarną część)
  • Wyznaczania rozkładu gęstości jądrowej (przez kształt potencjału optycznego)
  • Wejściowych danych dla modelu jądra złożonego (czynnik transmisji $T_l$)

Parametry modelu optycznego (głębokość studni $V_0 \sim 50$ MeV, absorpcja $W \sim 5$–$15$ MeV, promień $R = r_0 A^{1/3}$, dyfuzność $a \sim 0{,}65$ fm) zbiór w atlas fenomenologicznych parametrów (RIPL — Reference Input Parameter Library, IAEA).

Pomiar rozkładu ładunku i materii jądrowej — rozpraszanie elektronów

Elektrony, jako cząstki testowe bez silnego oddziaływania, są idealnym narzędziem do "mapowania" rozkładu ładunku elektrycznego w jądrze (co odpowiada rozmieszczeniu protonów). Technika: elektrody rozpraszanie (electron scattering) — wiązka elektronów o energii kilkudziesięciu–kilkuset MeV bombarduje tarczę jądrową.

Dla wymagalnego zasięgu $\Delta r \sim 0{,}1$ fm potrzebna długość fali de Broglie'a elektronu: $\lambda = \hbar c / E_{kin} \approx 0{,}2 \mathrm{fm}$ przy $E_{kin} \approx 1 \mathrm{GeV}$.

Rozkład ładunku jądra $\rho_{ch}(r)$ jest połączony z czynnikiem formikacji $F(q)$ przez transformatę Fouriera:

$$F(\vec{q}) = \int e^{i\vec{q}\cdot\vec{r}} \rho_{ch}(\vec{r})\, d^3r$$

Przekrój czynny Motty (elektron na punkcie, ze spin): $d\sigma/d\Omega = |F(q)|^2 \cdot (d\sigma/d\Omega)_{Mott}$.

Z pomiaru $d\sigma/d\Omega$ dla wielu kątów wyznacza się $|F(q)|^2$, a przez odwrotną transformatę Fouriera — $\rho_{ch}(r)$. Parametryzacja wynikowa: 3-parametrowy rozkład Fermi-Diraca (Saxon-Woods):

$$\rho(r) = \frac{\rho_0}{1 + \exp[(r-R)/a]}$$

gdzie $R = (1{,}11\text{–}1{,}14) A^{1/3} \mathrm{fm}$ i dyfuzność $a \approx 0{,}545 \mathrm{fm}$.

Ważne wyniki elektrody rozpraszania:

  • Potwierdzenie $R \propto A^{1/3}$ i stałości gęstości centralnej
  • Precyzyjny pomiar promienia ładunkowego protonu: $r_p = 0{,}8775 \mathrm{fm}$ (kontrowersja z mu-H 2010: $r_p = 0{,}84 \mathrm{fm}$, odchylenie ~4σ, rozwiązane w 2019)
  • Rozkłady ładunku izotopów Ca-40 vs Ca-48 (detekcja neutronowej skóry)
  • Stany wzbudzone jądra przez rozpraszanie nieelastyczne ($ee'$)

Projekt MAINZ A1: precyzyjne pomiary rformfaktora protonowego; MAMI (Niemcy) jako akseleratorowe centrum tych badań.

Bariera kulombowska i synteza pierwiastków w gwiazdach

Bariera Coulomba, którą pokonuje cząstka alfa przy rozpadzie, musi też być pokonana przez zderzające się jądra w reakcjach syntezy. To ten sam mechanizm — tunelowanie — ale w odwrotnym kierunku.

W gwiazdach głównego ciągu (jak nasze Słońce) temperatura wnętrza wynosi $T_c \approx 1{,}5 \times 10^7 \mathrm{K}$, co daje $kT \approx 1{,}3 \mathrm{keV}$. Tymczasem bariera Coulomba dla dwóch protonów wynosi:

$$V_C(R = 1 \mathrm{fm}) = \frac{e^2}{4\pi\varepsilon_0 \cdot 1 \mathrm{fm}} = \frac{1{,}44 \mathrm{MeV\cdot fm}}{1 \mathrm{fm}} = 1{,}44 \mathrm{MeV} \gg 1{,}3 \mathrm{keV}$$

Termiczne ogony rozkładu Maxwella-Boltzmanna dają kilka protonów o energiach kilkudziesięciu keV — ale wciąż są 20–30 razy za wolne, by klasycznie pokonać barierę.

Reakcje zachodzą przez tunelowanie Gamowa. W rozkładzie prędkości termicznych jest "okno Gamowa" — splot gorącego ogona MB z czynnikiem przenikania bariery Coulomba przez tunelowanie. To okno leży przy:

$$E_0 = \left(\frac{\pi \alpha Z_1 Z_2 m c^2}{\sqrt{2}} kT\right)^{2/3} \propto (Z_1 Z_2)^{2/3} T^{2/3}$$

Dla p-p w Słońcu: $E_0 \approx 6 \mathrm{keV}$. Efektywna stopa reakcji jest proporcjonalna do $\exp(-3E_0/kT)$ — silnie zależna od temperatury.

Konsekwencje:

  1. Reakcja p-p w Słońcu jest ekstremalnie wolna ($\tau_p \sim 10^{10}$ lat na jeden proton!) — stąd Słońce żyje tak długo
  2. Cięższe pierwiastki (C, O, Ne, Si) potrzebują wyższej temperatury ($T \sim 10^8$–$10^9$ K), osiągalnej w masywnych gwiazdach
  3. Synteza żelaza jest granicą — powyżej tego $Q < 0$ dla fuzji i tylko supernowe/AIC mogą produkować cięższe pierwiastki przez r-process

Precyzyjne pomiary przekrojów czynnych reakcji stellarnych przy niskich energiach (poniżej okna Gamowa, ale skalowalne) wykonuje się w laboratoriach podziemnych: LUNA (Laboratorio Nazionale del Gran Sasso, Włochy, 1400 m pod ziemią — gdzie tło promieniowania kosmicznego jest 10⁶ razy mniejsze niż na powierzchni). LUNA zmierzyła np. ${}^{14}$N$(p,\gamma){}^{15}$O z precyzją 8% przy energiach 70–232 keV — kluczową reakcję dla astrofizyki słonecznej (reguluje cykl CNO).

Siły jądrowe w gwiazdach neutronowych — równanie stanu materii

Własności gwiazd neutronowych zależą bezpośrednio od zachowania sił jądrowych w ekstremalnych gęstościach. Ciśnienie i energia wewnętrzna zimnej symetrycznej materii jądrowej (równe ilości protonów i neutronów) w funkcji gęstości to równanie stanu (EoS, Equation of State) — jest to fundamentalne wyjście z teorii sił jądrowych.

W normalnej gęstości jądrowej $\rho_0 \approx 0{,}16 \mathrm{fm}^{-3}$ (2,3 × 10¹⁷ kg/m³) energia na nukleon wynosi $E/A = -16 \mathrm{MeV}$ (właśnie tyle wynosi energia wiązania na nukleon w jądrach pośrednich, zgodnie z wzorem Weizsäckera). Ciśnienie w $\rho_0$: $P = 0$ (definicja normalnej gęstości jako minimum).

Co się dzieje przy $\rho > \rho_0$ (wewnątrz gwiazd neutronowych)?

  • Siły dwuciałowe (2N): same powodują kolaps — brakuje wystarczającego ciśnienia repulsywnego.
  • Siły trójciałowe (3N): dodają silną repulsję przy $\rho > 2\rho_0$, usztywniając EoS.

Bez sił $3N$ równanie stanu jest zbyt "miękkie" — maksymalna masa gwiazdy neutronowej wychodzi $M_{max} \lesssim 1{,}8 M_\odot$. Z realistycznymi siłami $3N$ (chiralEFT): $M_{max} \sim 2{,}2$–$2{,}5 M_\odot$ — w zgodzie z obserwacjami pulsarów ($M = 2{,}01 M_\odot$ dla J0348+0432).

Energia symetrii: w asymetrycznej materii (więcej neutronów niż protonów, jak w gwiazdach neutronowych) pojawia się dodatkowy koszt energetyczny:

$$E(\rho, \delta) \approx E_0(\rho) + S(\rho) \delta^2 + \ldots, \qquad \delta = (N-Z)/A$$

gdzie $S(\rho)$ — moduł symetrii. Jego wartość i ewolucja z gęstością decydują o wielkości neutronowej skóry ${}^{208}$Pb (mierzony przez PREX-II), promieniu gwiazdy neutronowej i mechanizmie wybuchu supernowej. Precyzyjne pomiary zarówno jądrowe (przekroje ładunkowe, rozpraszanie elektrony) jak i astrofizyczne (NICER, fale grawitacyjne) starają się określić $S(\rho)$ do kilku procent.

Eksperymentalne badanie sił jądrowych — rozpraszanie nukleon-nukleon

Jak w praktyce wyznacza się potencjał jądrowy? Fundamentalnym narzędziem jest analiza fazowa (phase-shift analysis) danych rozpraszania nukleon-nukleon.

Gdy neutron bombarduje proton w tarczy, amplituda rozpraszania rozkłada się na fale cząstkowe (partial waves) charakteryzowane przez momenty pędu orbitalny $L$ i całkowity $J$. W każdej fali cząstkowej rozpraszanie przejawia się jako przesunięcie fazowe $\delta_L$ — zmiana fazy funkcji falowej spowodowana potencjałem jądrowym. Związek z przekrojem czynnym:

$$\frac{d\sigma}{d\Omega} = |f(\theta)|^2 = \left|\sum_L (2L+1) e^{i\delta_L} \sin\delta_L P_L(\cos\theta) / k\right|^2$$

Zmierzone $d\sigma/d\Omega$ w różnych kątach i energiach pozwalają metodą fitowania odtworzyć $\delta_L(E)$ dla dziesiątek kanałów. Bazę takich faz dostarczają bazy danych: Nijmegen phase shift analysis (do 350 MeV), GWU SAID (George Washington University, Seattle Arndt-Gwinn database, do 3 GeV).

Kluczowe obserwacje z danych rozpraszania:

  1. W kanale ${}^1S_0$ (singlet $L=0$): duże ujemne przesunięcie fazowe przy niskiej energii — wirtualny stan, bardzo mała energia wiązania
  2. W kanale ${}^3S_1$ (triplet $L=0$): dodatnie przesunięcie, deuteron jest stanem związanym w tym kanale
  3. Kanał ${}^3P_0$: przy ~200 MeV energia zmienia się zachowanie — wskazuje na naturę odpychania twardo-jądrowego
  4. Mieszanie kanałów ${}^3S_1$-${}^3D_1$: parametr mieszania $\epsilon_1 \neq 0$ — bezpośredni dowód na tensor forces

Precyzja pomiaru: Nowoczesne eksperymenty (np. TRIUMF, PSI, Indiana, JINR) mierzą różniczkowe przekroje czynne z dokładnością ~0,5–2%, a analizy przekrojów absolutnych z ~1%. Przy tak dużej precyzji konieczna jest korekcja elektromagnetyczna (modyfikacja Coulomba dla p-p) i efekty relatywistyczne ($\geq 100$ MeV).

Prądy wymiany mezonów — magnesy jonowe

Klasyczny model: moment magnetyczny nukleonu pochodzi wyłącznie z jego spinu i ruchu orbitalnego. Ale kwarki wewnątrz nukleonu emitują i absorbują wirtulaowe piony między sobą, tworząc prądy wymiany mezonów (MEC, Meson Exchange Currents).

Prąd wymiany pionu: kiedy jeden nukleon emituje wirtualny pion, który przepływa do sąsiedniego nukleonu, a foton jest absorbowany w trakcie tego lotu — pojawia się dodatkowy prąd elektryczny, nieobecny w klasycznym modelu. Konsekwencje:

Moment magnetyczny deuteronu: Klasyczny model przewiduje $\mu_d = \mu_p + \mu_n = 0{,}880 \mu_N$. Zmierzony: $\mu_d = 0{,}857 \mu_N$. Różnica 2,6% pochodzi z MEC i komponentu D-wave — oba wkłady są obliczalne w ramach chiralnej EFT.

Przekrój czynny na fotodezintegrację deuteronu: $\gamma + d \rightarrow p + n$ (próg 2,224 MeV, energia wiązania deuteronu). MEC wnoszą ~15–20% do przekroju czynnego, szczególnie przy wyższych energiach. Precyzyjne pomiary (Compton-GDH Collaboration, HIγS w Duke) testują te obliczenia.

Reguła sumy Gerasimova-Drell-Hearna (GDH): Integralny związek między momentem anomalnym cząstki a absorpcją spolaryzowanych fotonów na spolaryzowanym nukleonie: $\int_0^\infty (\sigma_{1/2}-\sigma_{3/2}) dE/E = -2\pi^2 \alpha \kappa^2 / m^2$. Weryfikacja na SPring-8 (Japonia) i MAMI (Niemcy): MEC są kluczowe dla zrozumienia reguły sumy.

Tensor forces — asymetria przestrzenna siły jądrowej

Standardowe siły centralne zależą tylko od $r = |\vec{r}_1 - \vec{r}_2|$. Tensor forces zależą dodatkowo od kąta między spinem a wektorem łączącym:

$$S_{12} = 3(\vec{\sigma}_1 \cdot \hat{r})(\vec{\sigma}_2 \cdot \hat{r}) - \vec{\sigma}_1 \cdot \vec{\sigma}_2$$

Własności operatora $S_{12}$:

  • $\langle S_{12} \rangle = 2$ gdy spiny są wzdłuż $\vec{r}$ (oś połączenia)
  • $\langle S_{12} \rangle = -1$ gdy spiny są prostopadle do $\vec{r}$
  • $S_{12} = 0$ dla stanów z $L = 0$ (czyste stany S)

Dlatego tensory powodują mieszanie kanałów ${}^3S_1$ i ${}^3D_1$ w deuteronie — efekt opisywany przez parametr mieszania $\epsilon_J$. Bez tensora deuteron byłby w stanie czystym $S$; z tensorem wkład D-wave wynosi ~5%.

Tensor forces mają praktyczne konsekwencje dla obliczenia energii jąder:

  • Dla ${}^4$He: tensor daje znaczną (~5 MeV) część energii wiązania
  • W materii jądrowej: tensor produkuje anizotropowy składnik ciśnienia
  • Przy wysokich gęstościach (gwiazdy neutronowe): tensor może prowadzić do zmian symetrii fazowej (np. polaryzacji spinowej materii neutronowej)

Tensor forces są obecne przy wymianie $\pi^0$ (naturalna konsekwencja pseudoskalarowego sprzężenia Yukawa) i przy wymianie mezonu ρ (z odwrotnym znakiem). Subtelna cancellacja π i ρ w kanale $T = 1$ jest jednym z powodów, dla których "czysta" wymiana π nie jest wystarczająca.

Skalowanie Levingera — bariery kulombowskie w złożonych jądrach

W klasycznych ujęciach często zakłada się, że bariera Coulomba między jądrem a cząstką alfa o ładunku $Z_\alpha = 2$ jest $V_C = 2Ze^2/(4\pi\varepsilon_0 R)$, gdzie $Z$ jest liczbą protonów jądra. Jednak cząstka alfa jest wewnątrz jądra i bariera jest tworzona przez $Z - 2$ protonów (odejmując dwa protony cząstki alfa ze składu jądra).

Dla ciężkich jąder ($A \gg 4$): $V_C(R) \approx 2Z e^2 / (4\pi\varepsilon_0 R)$ — modyfikacja jest jedynie korekcją. Ale dla lekkich jąder to ważne: dla ${}^8\mathrm{Be}$ ($Z=4$): bariera kulombowska dla cząstki alfa jest produkowana przez tylko 2 pozostałe protony ($Z - 2 = 2$), stąd $V_C$ jest znacznie mniejsza niż naiwnie sądzimy. To tłumaczy, dlaczego ${}^8$Be jest nietrwałe i natychmiast (po $10^{-16}$ s) rozpada się na dwie alfa.

Korelacja Viola-Seaborg: Analogia prawa Geigera-Nuttalla rozszerzona na wszystkie aktynowce:

$$\log_{10} T_{1/2} = (a Z / \sqrt{E_\alpha}) + (b Z + c) + h_Z$$

gdzie $h_Z$ jest korekcją spinową (spin hindrance factor) dla przejść z/do nieparzystych jąder. Korelacja ta jest używana do przewidywania periodów półtrwania niezmierzonych transaktynowców, w tym pierwiastków superciężkich blisko wyspy stabilności.

Siły jądrowe a projektowanie reaktorów i osłon

Wiedza o przekrojach czynnych neutronów (która wynika z sił jądrowych przez mechanizm potencjału jądrowego i rezonansów Breit-Wignera) jest podstawą dla każdego obliczenia reaktorowego.

Moderacja neutronów: Zderzenie elastyczne neutronów z jądrami wodoru ($A = 1$) jest najefektywniejsze energetycznie — neutron może stracić nawet 100% energii kinetycznej w centralnym zderzeniu z protonem. Skalowanie energii: $\langle \Delta E / E \rangle = 2A/(A+1)^2$. Dla wodoru: 0,5 (50% na zderzenie); dla węgla-12: 0,142; dla uranu-238: ~0,008. Stąd woda i grafit są efektywnymi moderatorami, ciężkie jądra — nie.

Rozpraszanie neutronów w moderatorach: W reaktorach termicznych neutrony wielokrotnie rozpraszają się, zanim osiągną termalizację. Liczba zderzeń do termalizacji: $n \approx [\ln(E_0/E_{thermal})]/\xi$, gdzie $\xi = 1 + [(A-1)^2/2A]\ln[(A-1)/(A+1)]$ to średni logarytmiczny ubytek energii. Dla wodoru: $\xi = 1$, $n \approx 18$; dla grafitu: $\xi = 0{,}158$, $n \approx 114$.

Minimalizacja pochłaniania przez strukturalne materiały: Reakcja $(n,\gamma)$ na nuklidach strukturalnych (żelazo, cyrkon, aluminium) pochłania neutronsy i redukuje $k_{eff}$. Stąd:

  • Rury cyrkonowe w reaktorach LWR (przekrój na pochłanianie $\sigma_a(Zr) = 0{,}18$ b)
  • Aluminium w reaktorach doświadczalnych ($\sigma_a(Al) = 0{,}23$ b)
  • Unikanie stali roztworu w rdzeniu ($\sigma_a(Fe) = 2{,}6$ b, relatywnie duże)

Szczegółowe obliczenia transportu neutronów (Monte Carlo lub deterministyczne SN) wymagają precyzyjnych przekrojów czynnych z bibliotek jądrowych — co bezpośrednio powiązuje fenomenologię sił jądrowych (przez modele optyczne i jądra złożonego) z inżynierią reaktorów.

Hideki Yukawa i historia teorii sił jądrowych

Historyczny kontekst teorii sił jądrowych jest fascynujący — to historia dwóch pomyłek i jednego wielkiego odkrycia.

1932 — odkrycie neutronów (Chadwick) i praca Heisenberga sugerująca protonowo-neutronową budowę jądra. Ale co trzyma nuklony razem? Nie grawitacja (za słaba), nie elektromagnetyzm (odpycha protony). Potrzebna jest nowa siła.

1935 — Hideki Yukawa (Osaka) proponuje teorię sił jądrowych przenoszoną przez nową cząstkę. Klucz: zasięg siły wiąże się z masą nośnika przez zasadę nieoznaczoności Heisenberga:

$$\Delta E \cdot \Delta t \geq \hbar/2 \Rightarrow m_{nosnik} c^2 \approx \frac{\hbar c}{R}$$

Dla zasięgu $R \approx 1{,}4 \mathrm{fm}$:

$$m_{nosnik} c^2 \approx \frac{197{,}3 \mathrm{MeV\cdot fm}}{1{,}4 \mathrm{fm}} \approx 141 \mathrm{MeV}$$

Yukawa przewidział mezon o masie ~140 MeV — co jest niezwykłą trafnością jak na brakujący pomiar (pion ma masę 134,98/139,57 MeV). Praca ukazała się w 1935 roku w Proceedings of the Physico-Mathematical Society of Japan i początkowo przeszła niemal niezauważona.

1937 — w promieniowaniu kosmicznym odkryto cząstkę o masie ~105 MeV (myon). Wszystkim wydawało się, że to meson Yukawa — podobna masa, ładunek elektryczny jak pion... ale myon prawie nie oddziałuje z materią jądrową! Zamieszanie trwało dekadę.

1947 — C.F. Powell z grupą z Bristolu używając emulsji fotograficznych z Pic du Midi de Bigorre (Pireneje) wyjaśnił łamigłówkę: istnieją dwa różne rodzaje mezonów. Cięższy ($\pi$, pion) jest mezonem Yukawa i szybko rozpada się do myon + neutrino. Lżejszy (myon $\mu$) to lepton, niezwiązany z siłami jądrowymi. Cząstka Yukawa wreszcie znaleziona!

Nagrody Nobla: Yukawa (1949) za teoretyczne przewidzenie mezonu; Powell (1950) za metody fotograficzne i odkrycie pionu.

Ironia historii: I.I. Rabi po odkryciu myonu zapytał retorycznie: "Who ordered that?" (Kto to zamówił?). Myon nie był przewidziany żadną teorią i jest do dziś "zbędnym" duplikatem elektronu — co było pierwszą wskazówką, że natura ma trzy generacje cząstek.

Siły jądrowe jako resztkowe oddziaływanie kolorowe

Na głębszym poziomie siły jądrowe nie są fundamentalne — są resztkowym oddziaływaniem kolorowym (residual strong force), analogicznym do sił van der Waalsa w chemii.

Dokładna analogia:

  • Atom jest elektrycznie neutralny, ale mimo to oddziałuje z sąsiednimi atomami przez fluktuacje momentów dipolowych (siły van der Waalsa, London dispersion forces). Te siły są słabsze niż wewnątrzatomowe Coulomb, ale wiążą cząsteczki w cieczy i ciała stałe.
  • Nukleon jest kolorycznie neutralny (baryon = 3 kwarki różnych kolorów = "biały"), ale kwarki z jednego nukleonu "wyłażą" trochę poza jego granicę i fluktuacje kolorowego pola gluonowego tworzą efektywne przyciąganie z sąsiednim nukleonem.

Pion (i inne mezony) są "cząstkami Goldstone'a" spontanicznie złamanej symetrii chiralnej QCD. Są najlżejszymi hadronami i dlatego mają najdłuższy zasięg (odwrotnie proporcjonalny do masy). Ich wymiana jest najważniejszą częścią efektywnego opisu sił jądrowych przy dużych odległościach.

Krótkie zasięgi ($r < 1$ fm): kolorowe oddziaływania są nieperturbacyjne i wymiana wielu gluonów między kwarkami różnych nukleonów staje się dominująca — to jest "twarde jądro". Można je modelować przez wymianę ciężkich mezonów (ρ, ω) lub przez efektywne oddziaływania kontaktowe w chiralnej EFT.

Dlaczego wymiana pionów jest wystarczającym przybliżeniem? Bo pion jest lekki (135 MeV) i ma duże sprzężenie z nukleona ($g_{\pi NN} \approx 13$, czyli $g^2/4\pi \approx 14$). Cięższie mezony mają mniejszy zasięg i mniejsze sprzężenie — ich wkład jest korektą. W precyzyjnych obliczeniach konieczna jest wymiana dwóch pionów (box diagram), który daje poprawkę ~10–20% przy zasięgach 1–2 fm.

Siły spin-orbita i ich konsekwencje

Człon spin-orbita potencjału nukleonowego:

$$V_{LS} = V_{LS}(r)\, \vec{L} \cdot \vec{S}$$

jest kluczowy dla modelu powłokowego jądra (patrz artykuł o budowie jądra). Ale skąd pochodzi to oddziaływanie?

W QFT (relatywistycznej teorii) człon spin-orbita wynika z relatywistycznej kinetyki cząstki w polu centralnym (Thomas term): cząstka poruszająca się w polu elektrycznym $\vec{E}$ "widzi" w swoim układzie spoczynkowym pole magnetyczne $\vec{B}' \propto \vec{v} \times \vec{E}$, które sprzęga się ze spinowym momentem magnetycznym.

W równaniu Diraca (relatywistycznym odpowiedniku Schrödingera dla fermionów) człon spin-orbita pojawia się naturalnie. Dla potencjału Saxon-Woods:

$$V_{SO}(r) = \frac{1}{r}\frac{dV}{dr}\, \vec{L}\cdot\vec{S}$$

Znaki i wielkości $V_{LS}$ są dobrze znane z danych rozpraszania n-p i p-p (analiza fazowa) i z własności powłokowych ciężkich jąder. W modelu powłokowym dla nuklonów (model Goeppert-Mayer) właśnie siła $V_{LS}$ jest źródłem właściwych liczb magicznych (28, 50, 82, 126) — bez niej prosty oscylator daje błędne magiczne liczby (2, 8, 20, 40, 70...).

Siły trójciałowe — klucz do zrozumienia jąder lekkich

Oddziaływanie nukleon-nukleon jest fundamentem, ale kiedy mamy trzy lub więcej nukleonów, pojawiają się nieaddytywne siły trójciałowe (3N Forces), które nie mogą być zredukowane do sumy parowych oddziaływań.

Fizyczny mechanizm sił 3N: Nukleon A emituje pion, który zostaje wchłonięty przez nukleon B, który ulega wewnętrznemu wzbudzeniu do stanu $\Delta(1232)$ (barion o spinie 3/2, masa 1232 MeV), a następnie emituje kolejny pion pochłaniany przez nukleon C. Ten proces jest kołowy i nie jest zawarty w żadnym parowym oddziaływaniu A-B, B-C czy A-C z osobna — to jest Fujita-Miyazawa 3N force (1957).

Ile siły 3N wnosi?

Dla helu-4 ($A = 4$):

  • Tylko siły 2N: $E_{B}({}^4\mathrm{He}) = -26{,}6 \mathrm{MeV}$ (zmierzone: $-28{,}3$ MeV) — za mało związane
  • Siły 2N + 3N: $E_B({}^4\mathrm{He}) = -28{,}5 \pm 0{,}5 \mathrm{MeV}$ — doskonała zgodność

Dla tlenu-16 ($A = 16$):

  • Bez 3N: energia stanu podstawowego za duża (zbyt silne wiązanie)
  • Z 3N: poprawna energia

Siły 3N są mniej poznane niż 2N: wymagają większej liczby parametrów LEC (chiralnej EFT) i mają większą niepewność. Są jednak absolutnie kluczowe dla:

  1. Poprawnego odtworzenia energii lekkich jąder
  2. Saturacji materii jądrowej (właściwego minimum energii i gęstości przy $\rho_0$)
  3. Twardości równania stanu gwiazd neutronowych
  4. Progów rozpraszania i energii wiązania

Siły czterociałowe i wyższe: prawdopodobnie istnieją, ale ich wkład jest szacowany na <0,5% dla $A \leq 12$ i jest trudny do wydzielenia eksperymentalnie.

Chiral perturbation theory — nowsze podejście do sił jądrowych

Zamiast fenomenologicznych potencjałów dopasowanych do danych, chiralna teoria perturbacyjna (ChPT, chiral EFT) wywodzi siły jądrowe systematycznie z Lagrangianu efektywnej teorii pola.

Idea jest prosta: w reżimie $E \ll \Lambda_{ChPT} \approx 1 \mathrm{GeV}$ (co jest spełnione dla fizyki jądrowej), QCD redukuje się do teorii z pionami jako dynamicznymi stopniami swobody i nukleonami jako ciężkimi źródłami. Lagrangian jest pisany w postaci rozwinięcia:

$$\mathcal{L}_{ChEFT} = \mathcal{L}_{\pi\pi} + \mathcal{L}_{\pi N} + \mathcal{L}_{NN} + \mathcal{L}_{NNN} + \ldots$$

Każdy człon jest dopuszczony przez symetrie: izospinową (wymaga izospinowej niezmienniczości), Lorentza i parzystość/odwracanie czasu. Rozwijamy w potęgach $p/\Lambda_{ChPT}$ (gdzie $p$ to pęd lub masa pionu):

  • LO (Leading Order): kontaktowe oddziaływanie ${}^1S_0$ i ${}^3S_1$-${}^3D_1$ + wymiana jednego pionu (OPEP)
  • NLO (Next-to-Leading Order): pierwsza korekcja do wymiany 2π, prądy 2N kontaktowe ze spinami
  • N2LO: pełna wymiana 2π + pierwsze siły $3N$ (Fujita-Miyazawa)
  • N3LO: wymiana 2π z pętlą, siły $3N$ (dalsze człony)
  • N4LO: wymiana 3π, siły 4N

Parametry każdego rzędu (Low Energy Constants, LECs) są wyznaczane z dopasowania do rozpraszania π-N lub N-N. W wyższych rzędach jest tych parametrów więcej, ale precyzja rośnie.

Osiągnięcia chiral EFT:

  • Reprodukcja faz rozpraszania n-p i p-p do 250 MeV z $\chi^2/\mathrm{datum} \approx 1$ przy N3LO (E. Epelbaum i in., 2002)
  • Poprawne energie i promienie jąder $A \leq 25$ z Ab Initio obliczeniami NCSM z SRG-evolved N3LO+N2LO(3N) siłami
  • Saturation of nuclear matter przy $\rho_0 = 0{,}16 \mathrm{fm}^{-3}$, $E/A = -16$ MeV bez fine-tuningu

Niepewności chiralEFT: rozwinięcie jest asymptotyczne, nie zbieżne w sensie analitycznym. Przy N3LO i wyżej, obliczenia pokazują dobrą zbieżność dla $A \leq 16$. Dla ciężkich jąder ($A > 100$) i wysokich gęstości ($\rho > 2\rho_0$) chiralEFT traci kontrol nad niepewnościami — tam potrzebne są inne podejścia (metamodel EoS lub parametryzacje).

Siły jądrowe poza pionową wymianą — rola cięższych mezonów

Wymiana jednego pionu wyjaśnia długozasięgową część potencjału ($r > 2$ fm) i tensor. Dla krótszych odległości potrzebne są cięższe mezony:

  • Meson rho ($\rho^0, \rho^\pm$, $m_\rho = 775 \mathrm{MeV}$): Wektorowy, spin 1. Wymiana $\rho$ produkuje krótszą wersję tensorowego oddziaływania (o odwrotnym znaku niż π!) i silny człon spin-orbita. Zasięg $\hbar / m_\rho c \approx 0{,}26 \mathrm{fm}$.
  • Meson omega ($\omega^0$, $m_\omega = 782 \mathrm{MeV}$): Izoskalarny wektorowy. Odpowiada za krótszą repulsję (twarde jądro) i człon $\vec{\tau}_1 \cdot \vec{\tau}_2$-niezależny. Zasięg $\sim 0{,}25 \mathrm{fm}$.
  • Meson sigma ($\sigma$ lub $f_0(500)$, $m_\sigma \approx 400$–$550 \mathrm{MeV}$): Skalarny, izoskalarny — daje średniozasięgowe przyciąganie. Jego status jako cząstki elementarnej jest dyskusyjny (może być para pionów w stanie rezonansowym).

W modelu OBE (One Boson Exchange) potencjał jądrowy jest sumą wymian π, ρ, ω, σ (i ewentualnie η, δ). Różne modele OBE (Paris, Bonn, Nijmegen) różnią się szczegółami doboru mezonów i stałych sprzężenia. Wszystkie jednak dobrze odtwarzają dane fazowe rozpraszania do energii $\sim 300$ MeV i własności deuteronu z kilkoma wolnymi parametrami.

Masy i wiązania jądra złożonego — dlaczego twardy rdzeń ma znaczenie biologiczne i reaktorowe

Krótko- i średniozasięgowe własności sił jądrowych mają bezpośrednie, praktyczne konsekwencje w kilku dziedzinach:

1. Bezpieczeństwo reaktorów jądrowych — współczynnik Dopplera

Przy wzroście temperatury paliwa uranowego (UO₂) jądra uranu drgają szybciej (efekt Dopplera). Ich prędkość relatywna w stosunku do neutronów zmienia się, co poszerza rezonanse neutronów w jądrze U-238 (rezonanse absorpcyjne w zakresie 1–10 eV). Szersze rezonanse pochłaniają więcej neutronów — reaktywność spada. To jest ujemny współczynnik Dopplera ($\alpha_T < 0$), fundamentalny mechanizm bezpieczeństwa pasywnego reaktorów lekkowodnych.

Głębsza fizyka: rezonanse w U-238 to stany złożone (Breit-Wigner) wynikające ze sił jądrowych między 239 nukleonami w złożonym jądrze U-239*. Gęstość poziomów energetycznych i szerokości Γ_n, Γ_γ wynikają wprost z siły wiązania i struktury poziomów — a te z potencjału jądrowego.

Parametr: $\alpha_T \approx -2$ do $-4 \times 10^{-5} \mathrm{K}^{-1}$ dla LWR. Obliczany kodami NJOY/SERPENT/OpenMC z bibliotek ENDF/B lub JEFF.

2. Medycyna — napromieniowanie terapeutyczne i uszkodzenia DNA

Biologiczne działanie promieniowania jonizującego — terapia protonowa, terapia jonami węgla (HIT/MedAustron) — opiera się na tym, że protony/jony węgla w tkance tracą energię według krzywej Bragga, z pikiem głęboko w guzie. Kształt krzywej Bragga zależy od:

  • Elastycznego i nieelastycznego rozpraszania protonu na jądrach tkanek (C, N, O) — rządzą nim siły jądrowe
  • Oddziaływań elektromagnetycznych (hamowanie Bethego-Blocha)
  • Reakcji jądrowych fragmentacji jonów węgla — produkty (α, p, Li) "rozmazują" dawkę poza szczyt Bragga, co jest nieodłącznym skutkiem sił jądrowych w kolizjach C-12 z jądrami tkankow

Siły jądrowe determinują też LET (Linear Energy Transfer) — miara gęstości jonizacji wzdłuż toru cząstki. Jony węgla mają LET ~70 keV/μm, ok. 30× więcej niż promieniowanie X — stąd wyższa względna biologiczna skuteczność (RBE~3) i możliwość leczenia guzów radioopornych.

3. Fizyka neutrin i oddziaływania słabe — jakie to ma związek z siłami jądrowymi?

Siły jądrowe determinują energię wiązania jądra i strukturę poziomów. To z kolei kształtuje dostępność kanałów beta-rozpadu i wychwytywania elektronów:

  • Podwójny beta-rozpad $0\nu\beta\beta$ (jeśli istnieje): wrażliwy na masę efektywną neutrina Majorany $\langle m_{\beta\beta} \rangle$ oraz na nuclear matrix element (NME) — macierz przejścia obliczaną z dokładnej znajomości funkcji falowych jądra przy pomocy QRPA lub shell model. NME zależy od sił jądrowych (szczególnie tensora i sił korekcji jądrowej — SRC)
  • Supernowa rdzeniowa (CCSN): po zapaści rdzenia ($\rho > \rho_0$) neutrina emitowane z protogwiazdy neutronowej przechodzą przez otoczkę. Oddziaływanie $\nu + A \to \nu' + A'$ (koherentne) i $\nu + n \to p + e^-$ (absorpcja) zależy od sił jądrowych
  • Detekcja neutrin geotermalnych i reaktorowych: przekroje czynne na odwrotny beta-rozpad $\bar{\nu}_e + p \to n + e^+$ są wyznaczane precyzyjnie i zależą od różnic energii wiązania

4. Technologia broni — niskie współczynniki kc w głowicach

(W kontekście czysto fizycznym, bez szczegółów technologicznych): Minimalna masa krytyczna materiału rozszczepialnego (Pu-239, U-235) wynosi kilka kg przy optymalnym odbiciu. Wartości te wynikają z równowagi między neutronami wytworzonymi ($\bar{\nu} \approx 2{,}5$ na rozszczepienie) a tracącymi się do zewnątrz. Geometria i właściwości reflektora są obliczane przez rozwiązanie równania transportu Boltzmanna — którego parametry wejściowe (przekroje czynne totalne, elastyczne, nieelastyczne) wynikają ze sił jądrowych i muszą być mierzone doświadczalnie lub obliczane z Ab Initio.

Kwantowe stany rezonansowe — stany Gamowa i ich rola w rozpadach

Klasyczne rozwiązanie równania Schrödingera daje stany związane (energia ujemna, dyskretne) i stany ciągłe (energia dodatnia, rozpraszanie). Jednak dla procesów tunelowania i rezonansów użyteczne są stany Gamowa (Gamow states, 1928) — złożone rozwiązania z kompleksową energią:

$$E = E_r - \frac{i\Gamma}{2}$$

gdzie $E_r$ to energia rezonansu, $\Gamma = \hbar / \tau$ to szerokość związana z czasem życia. Stany Gamowa spełniają warunek brzegowy tylko fal wychodzących ($e^{+ikr}$ dla dużych $r$), nie normalizują się w zwykłym sensie (rozbiegają jak $e^{\kappa r}$ dla $r \to \infty$), ale tworzą kompletną bazę Berggrena.

Fizyczne znaczenie:

  • Alfa-emiter: $E_r$ = energia kinetyczna cząstki α w nieskończoności; $\Gamma$ = szerokość poziomów, którą obserwuje się w rozpraszaniu α na jądrze córki
  • Proton-emiter: jądra bogate w protony jak Pa-151 ($T_{1/2} = 1{,}73 \mathrm{ms}$) są proton-niestabilne — stany protonowe powyżej bariery konfiguracyjnej
  • Stany neutronowe przy neutron-dripline: Li-10 to nieliniowy rezonans neutronowy o $E_r = 0{,}4 \mathrm{MeV}$, $\Gamma = 2{,}0 \mathrm{MeV}$ — szerszy niż daleki od granicy stabilności

Metody obliczeniowe: rozwiązanie równania Schrödingera z potencjałem Saxon-Woodsa używając bazy Berggrena; zastosowanie w kodach GAMOW-SHELL-MODEL, GSM. Sprzęganie z kontinuum ma szczególne znaczenie dla egzotycznych jąder blisko linii dripline, gdzie stany związane "przeciekają" do kontinuum.

Skalowanie sił jądrowych z liczbą atomową — granica jądrowa i grawitacja

Siły jądrowe są siłami nasycającymi: energia wiązania rośnie liniowo z $A$ (a nie jak $A^2$ jak van der Waals czy grawitacja). To nasycenie ma głębokie konsekwencje dla granic istnienia jąder i materii jądrowej:

Dlaczego nie ma jąder dowolnie dużych?

  1. Bariera elektrostatyczna rośnie jak $Z^2$ (odpychanie protonów), siły jądrowe rosną jak $A \cdot \delta(\text{parzyste-nieparzyste})$. Dla $Z > 92$ jądra stają się coraz mniej stabilne.
  2. Dla ekstremalnie ciężkich jąder ($Z > 120$) samo-samoistna fizjal staje się energetycznie korzystna — jądro spontanicznie się dzieli w czasach $< 10^{-21}$ s.
  3. W neutronowych gwiazdach: materia jądrowa przy $\rho \sim 2$–$3 \rho_0$ może tworzyć "nuklearne makarony" (pasta phases): cylindryczne struktury (spaghetti) i warstwy (lasagne), gdyż siły jądrowe i Coulomba minimalizują energię w tej egzotycznej konfiguracji. Obliczenia QMD (Quantum Molecular Dynamics) pokazują te fazy i ich wpływ na lepkość i przewodnictwo cieplne skorupy gwiazdy neutronowej.
  4. Granica "wyspy stabilności" przy $Z \approx 114$, $N \approx 184$: domknięte powłoki protonowe i neutronowe stabilizują jądra supercięzkie — obliczone metodami Hartree-Fock-Bogolyubov lub FRDM z mikroskopowymi siłami jądrowymi.

Hipoteza materii kwarkowej w centrum gwiazd neutronowych: jeśli $\rho > 5$–$8 \rho_0$, hadrony mogą "stopić się" do fazy kwarkowo-gluonowej. Równanie stanu (EoS) przy tych gęstościach jest niepewne, bo siły jądrowe z chiral EFT tracą ważność, a perturbacyjne QCD jeszcze nie daje dobrej kontroli. Właśnie dlatego pomiary mas i promieni gwiazd neutronowych (NICER, LIGO/Virgo) są kluczowe — dostarczają ograniczeń na EoS przy ultra-wysokich gęstościach.

Ćwiczenia praktyczne

Pierwsze ćwiczenie: przygotuj trzy krótkie definicje: krótki zasięg sił jądrowych, niezależność ładunkowa i nasycenie. Do każdej dopisz, jaki problem w budowie jądra ta cecha pomaga zrozumieć.

Drugie ćwiczenie: porównaj siłę jądrową z oddziaływaniem elektromagnetycznym w czterech punktach: zasięg, znak oddziaływania, zależność od ładunku i rola w stabilności jąder.

Trzecie ćwiczenie: wyjaśnij, dlaczego nasycenie sił jądrowych jest potrzebne do zrozumienia prawie stałej gęstości materii jądrowej. Odpowiedź powinna być jakościowa, bez udawania dokładnego modelu potencjału.

Przejdź do ćwiczenia interaktywnego