Streszczenie

Samo osiągnięcie stanu nadkrytycznego nie mówi jeszcze, jak silny będzie wybuch jądrowy. Liczy się również to, jak szybko rośnie liczba neutronów i jak długo rdzeń pozostaje dostatecznie zwarty, zanim zacznie się hydrodynamicznie rozbiegać. Właśnie takie pytanie próbowało uporządkować wczesne równanie Bethego-Feynmana: nie „czy bomba zadziała”, lecz od czego w przybliżeniu zależy jej sprawność.1,2

Najważniejszy wniosek z tego typu modeli jest bardzo praktyczny. Wydajność rośnie silnie wraz z tempem mnożenia neutronów oraz z zapasem nadkrytyczności osiągniętym przed rozpadem układu. To dlatego w projektowaniu bomby tak ważne są: duża kompresja, mały wyciek neutronów, właściwy moment inicjacji oraz tamper, który choćby przez dziesiątki nanosekund utrzymuje rdzeń w korzystnym stanie. W skali wybuchu to wystarcza, by dodać wiele kolejnych pokoleń rozszczepień.1,3

Rozszerzenie tematu

Równanie Bethego-Feynmana należy do najwcześniejszych prób analitycznego opisu sprawności ładunku rozszczepieniowego. Przygotowali je Hans Bethe i Richard Feynman w Berkeley w 1942 roku na gruncie modelu jednorodnej ekspansji rdzenia.2 Nie było to jeszcze narzędzie do pełnego, dokładnego przewidywania uzysku, lecz sposób uporządkowania najważniejszych zależności między neutroniczną szybkością reakcji a hydrodynamicznym rozpadem układu.

Punktem wyjścia jest tempo wzrostu populacji neutronów opisane przez współczynnik mnożenia czasowego $\alpha$. Neutronową populację można zapisać jako:

$$N(t) = N_0 e^{\alpha t}$$

co oznacza, że dodatnie $\alpha$ odpowiada układowi nadkrytycznemu, a wielkość $1/\alpha$ wyznacza charakterystyczny czas jednego e-foldingu wzrostu.1 Im większe $\alpha$, tym szybciej rośnie moc i tym więcej energii układ zdąży wyzwolić, zanim jego geometria zacznie się nieodwracalnie psuć.

W przybliżeniu podanym przez Gdenarza jedna z form równania Bethego-Feynmana może być zapisana jako:

$$\eta = K R^2 \alpha^2 \delta$$

gdzie $\eta$ oznacza wydajność, $K$ jest stałą zależną od materiału rozszczepialnego, $R$ promieniem rdzenia w chwili maksymalnej nadkrytyczności, $\alpha$ tempem reakcji w tej samej chwili, a $\delta$ mierzy zapas nadkrytyczności względem promienia krytycznego.3 Sam zapis bywa w literaturze reformułowany na kilka sposobów, ale wspólny sens pozostaje ten sam: sprawność rośnie bardzo silnie wraz ze wzrostem tempa mnożenia i nadmiaru nadkrytycznego rozmiaru rdzenia.

To jest sedno, które warto zrozumieć zanim zacznie się śledzić szczegóły wyprowadzenia. Wydajność bomby nie zależy liniowo od pojedynczego parametru typu „masa materiału”. Zależy od całej sekwencji:

  1. jak szybko rdzeń przechodzi w stan natychmiastowej krytyczności,
  2. jakie maksimum $\alpha$ osiąga przy danej gęstości i geometrii,
  3. jak długo utrzymuje się ten korzystny stan zanim zacznie dominować rozbieg hydrodynamiczny,
  4. jak duży jest zapas ponad rozmiar dokładnie krytyczny.1,2,3

W tym sensie równanie Bethego-Feynmana jest pomostem między fizyką neutronów a fizyką ruchu materii. Z jednej strony potrzebny jest szybki wykładniczy wzrost liczby rozszczepień. Z drugiej strony ten wzrost sam produkuje energię, która zaczyna rozsadzać rdzeń. Wybuch jądrowy jest więc wyścigiem między mnożeniem neutronów a dezintegracją geometryczną układu. Cała sprawność zależy od tego, jak daleko reakcja zdąży zajść, zanim rdzeń przestanie być dostatecznie zwarty.2

Właśnie dlatego tak duże znaczenie ma predetonacja. Gdenarz podkreśla, że przy preinicjacji wzrost współczynnika $\alpha$ zostaje zatrzymany wcześniej, zanim układ zdąży dojść do projektowanego maksimum nadkrytyczności.3 To bezpośrednio obniża sprawność i końcową energię wybuchu. Nie chodzi więc tylko o „za wczesny start” reakcji, lecz o to, że cały rdzeń rozpada się wtedy przy zbyt małej wartości najbardziej korzystnych parametrów.

Wczesne równanie Bethego-Feynmana opierało się na dość surowym założeniu, iż rdzeń podczas deinsercji zachowuje w przybliżeniu stałą gęstość.2 To założenie upraszcza rachunek, ale ogranicza dokładność. Dlatego później pojawiały się inne równania sprawności, na przykład formy rozwijane przez Roberta Serbera, które inaczej traktowały zależność od nadkrytycznego nadmiaru i przebiegu ekspansji.2 Z dydaktycznego punktu widzenia nie zmienia to jednak najważniejszego obrazu: wszystkie te modele próbują odpowiedzieć na to samo pytanie, jak związać szybkość mnożenia neutronów z czasem rozpadu rdzenia.

Praktyczna konsekwencja jest bardzo konkretna. Jeżeli wydajność zależy silnie od $\alpha$, to trzeba maksymalizować wszystko, co podnosi efektywne tempo mnożenia. Oznacza to przede wszystkim:

  1. wysokie wzbogacenie lub użycie materiału o korzystnych własnościach neutronicznych,
  2. dużą kompresję, która zwiększa gęstość i zmniejsza wyciek neutronów,
  3. dobrą geometrię rdzenia,
  4. ograniczenie przedwczesnej inicjacji,
  5. użycie reflektora i tampra, które jednocześnie ograniczają wyciek neutronów i spowalniają ekspansję mechaniczną.2,3

To ostatnie jest szczególnie ważne. Tamper działa podwójnie. Neutronicznie odbija część neutronów z powrotem do rdzenia, podnosząc skuteczność reakcji. Mechanicznie dostarcza bezwładności, dzięki której fala rozprężania nie niszczy geometrii rdzenia natychmiast po rozpoczęciu wydzielania energii. W skali makroskopowej dziesiątki nanosekund brzmią jak nic. W bombie to jednak wystarczająco długo, by zajść mogło wiele dodatkowych pokoleń rozszczepień, a więc by uzysk wzrósł skokowo.2,4

Równanie Bethego-Feynmana dobrze porządkuje też intuicję, dlaczego implozja była tak cenna. Implozja nie tylko obniża masę krytyczną, ale pozwala wejść na wyższe wartości $\alpha$ i większy zapas nadkrytyczności $\delta$ zanim rozpocznie się dezintegracja układu. To właśnie dlatego silnie skompresowany rdzeń może mieć nieporównanie większą sprawność niż układ tylko nieznacznie nadkrytyczny, nawet jeśli oba formalnie „działają”.

Trzeba przy tym pamiętać, że nie jest to równanie dokładne w nowoczesnym sensie projektowym. Nie zastępuje pełnej hydrodynamiki, transportu neutronów, zmian fazowych materiału, asymetrii implozji ani strat energii. Jest jednak niezwykle użyteczne jako pierwszy model skali zjawiska. Pokazuje, dlaczego fizycy Projektu Manhattan musieli myśleć jednocześnie o dynamice rdzenia, czasie inicjacji i właściwościach materiałów, a nie tylko o samej liczbie kilogramów materiału rozszczepialnego.2,5

Najkrótsze podsumowanie jest więc takie: równanie Bethego-Feynmana jest wczesnym równaniem sprawności wybuchu jądrowego. Nie mówi wszystkiego, ale trafnie wskazuje trzy rzeczy, które trzeba maksymalizować, jeśli chce się uzyskać duży uzysk: szybkie mnożenie neutronów, duży zapas nadkrytyczności i dostatecznie powolny rozbieg hydrodynamiczny. Cała reszta konstrukcji bomby jest w dużej mierze podporządkowana właśnie temu celowi.1,2,3

Historia — Bethe i Feynman w Los Alamos

Hans Bethe i Richard Feynman trafili do Los Alamos w pierwszej połowie 1943 roku, kilka miesięcy po formalnym uruchomieniu laboratorium przez J. Roberta Oppenheimera. Bethe, urodzony w Strasburgu w 1906, był wtedy jednym z czołowych fizyków teoretycznych na świecie — znany z prac nad przejściami jadrowymi i cyklem proton-proton w gwiazdach (za co dostał Nagrodę Nobla w 1967). Feynman, młody absolwent Princeton, był dopiero na dorobku — lecz wyróżniał się wyjątkową intuicją fizyczną i umiejętnością rozwiązywania problemów, które inni uważali za za trudne.

Bethe objął kierownictwo Wydziału Teoretycznego (T-Division) Los Alamos. Feynman pracował w tym samym wydziale, początkowo jako młodszy fizyk. Równanie sprawności wybuchu rozszczepieniowego, znane dziś jako równanie Bethego-Feynmana, zostało wyprowadzone i opisane w wewnętrznym raporcie Los Alamos w 1943 roku. Raport jest do dziś sklasyfikowany w całości lub w znacznej części, choć jego główne wyniki zostały zrekon­struowane i opisane w literaturze otwartej po odtajnieniach lat 90.

Feynman w swoich wspomnieniach (famously in Surely You're Joking, Mr. Feynman!) opisywał pracę w Los Alamos z mieszaniną fascynacji i trwogi: robił obliczenia dotyczące broni, nie zawsze wiedząc, jakie praktyczne konsekwencje będą mieć. Bethe pozostawił bardziej bezpośrednie opisy naukowe, wskazując, że kluczowym wyzwaniem było powiązanie tempo neutronowej lawiny z czasem utrzymania geometrycznej zwartości rdzenia.

Obaj pracowali w epoce, gdy komputery cyfrowe niemal nie istniały — w Los Alamos obliczenia wykonywano mechanicznie, na kalkulatorach, przez setki osób jednocześnie podzielonych na grupy, gdzie każda liczyła jeden krok sekwencji. Feynman był jednym z organizatorów tej „ludzkiej maszyny obliczeniowej" i wskazywał na analogie z prymitywnymi procesorami paralelnymi.6

Wyprowadzenie — podstawy analityczne

Wyjście od pierwszych zasad pozwala zrozumieć, skąd biorą się poszczególne człony równania Bethego-Feynmana.

Krok 1: Dynamika neutronów

Neutronowa lawina w rdzeniu rozszczepialnym rośnie wykładniczo ze stałą $\alpha$:

$$N(t) = N_0\,e^{\alpha t}$$

Stałą $\alpha$ można wyrazić przez parametry reaktora:

$$\alpha = \frac{(k - 1)}{\ell}$$

gdzie $k$ jest efektywnym współczynnikiem mnożenia neutronów, a $\ell$ efektywnym czasem jednego pokolenia. W broni jądrowej z neutronami prędkimi $\ell \sim 10^{-8}\,\mathrm{s}$, a $k > 1$ dzięki nadkrytycznej masie i kompresji.

Krok 2: Dezinsercja hydrodynamiczna

Rdzeń zaczyna się rozpadać hydrodynamicznie, gdy ciśnienie generowane przez energię rozszczepień przekracza siły wiążące. W bardzo uproszczonym modelu jednorodnej ekspansji czas rozpadania jest rzędu czasu, w którym fala rozprężania przebiega przez rdzeń:

$$t_{dis} \sim \frac{R}{v_s}$$

gdzie $R$ jest promieniem rdzenia, a $v_s$ prędkością dźwięku (lub fali uderzeniowej) w urozszczepialnym materiale. Dla U-235 i Pu-239 przy gęstości normalnej $v_s \sim 3{-}4\,\mathrm{km/s}$, co daje $t_{dis} \sim 1{-}10\,\mu\mathrm{s}$ dla kul o rozmiarach kilku centymetrów.

Krok 3: Energia wybuchu

Liczba rozszczepień do chwili $t_{dis}$ wynosi:

$$N_{total} = \int_0^{t_{dis}} N_0 e^{\alpha t}\,dt \approx \frac{N_0 e^{\alpha t_{dis}}}{\alpha}$$

co dla dużego $\alpha t_{dis}$ jest zdominowane przez ostatnie chwile przed dezinsercją. Energia wybuchu jest proporcjonalna do $N_{total}$, więc:

$$E \propto \frac{e^{\alpha t_{dis}}}{\alpha}$$

a sprawność (stosunek energii do maksymalnej możliwej, gdyby wszystkie atomy rozszczepiły się) zależy głównie od $\alpha$ i od tego, jak daleko rdzeń zdąży zajść przed dezinsercją.

Krok 4: Zależność od nadkrytyczności i wymiarów

W modelu Bethego-Feynmana $\alpha$ jest proporcjonalne do nadmiaru nadkrytyczności $\delta = (R - R_c)/R_c$, gdzie $R_c$ to promień krytyczny. To daje ostatecznie sprawność proporcjonalną do:

$$\eta \propto \alpha^2 \sim \delta^2$$

co wyjaśnia silną zależność od nadmiaru ponad próg krytyczny. Podwójne $\alpha^2$ wynika z tego, że zarówno intensywność eksponencjalnego wzrostu, jak i czas dostępny do wzrostu (przez wpływ $\alpha$ na dynamikę dezinsercji) zależą liniowo od $\alpha$.

Całe wyprowadzenie jest znacznie bardziej skomplikowane w szczegółach (uwzględnia zmiany $k$ podczas ekspansji, termoelastyczne właściwości materiału, wyciek neutronów z rdzenia, efekty tampera), ale główna zależność $\eta \propto \alpha^2 \delta$ pochodzi właśnie z tej podstawowej analizy.6

Porównanie z modelem Serbera

Robert Serber, fizyk Princeton pracujący w Los Alamos od początku (znany z tzw. Primer, czyli „podręcznika dla nowych pracowników"), sformułował nieco inne przybliżenia do tego samego problemu.

Model Serbera kładł większy nacisk na szczegółowe obliczenie transmisji neutronów przez rdzeń (z uwzględnieniem skończonego biegu swobodnego) oraz na czas spędzonego w warunkach nadkrytycznych podczas kompresji. W jego ujęciu sprawność była wyrażana przez tzw. sprawność fissile material (fraction of fissile material actually fissioned), zwaną dziś symbolem $f$:

$$f = \frac{E_{actual}}{E_{max}} = \frac{\mathrm{Liczba\, rozszczepionych\, jąder}}{\mathrm{Łączna\, liczba\, jąder\, rozszczepialnych}}$$

Dla Little Boy (bomba uranowa zdetonowana nad Hiroszimą) szacowana sprawność Serbera wynosiła zaledwie ok. 1–2%, co odpowiadało ~0,85 kg U-235 faktycznie rozszczepionych z ~64 kg załadowanego uranu. To katastrofalna z inżynierskiego punktu widzenia marnotrawstwo — ale wynikało z bardzo niskiej $\alpha$ w układzie działowym (gun-type), który nie osiągał maksymalnej nadkrytyczności przed dezinsercją.

Fat Man (bomba plutonowa na Nagasaki), dzięki implozji osiągał $f \approx 17{-}20\%$ — przyzwoita sprawność jak na pierwsze urządzenie implozy. Trinity test (16 lipca 1945) dał podobne wyniki.

Różnica jest bezpośrednią konsekwencją równania Bethego-Feynmana: implozja gwałtownie zwiększa $\alpha$ przez zwiększenie gęstości (mniejszy wyciek, większy $k$) i zmniejszenie czasu dezinsercji, a jednocześnie całość odbywa się w krótszym czasie, ograniczając straty przez preinicjację.6

Sprawność historycznych urządzeń

Dla kontekstu warto porównać sprawności różnych historycznych urządzeń i co z nich wynika:

Urządzenie Typ Materiał Uzysk Masa materian. Sprawność
Little Boy działowy U-235 ~15 kt 64 kg U-235 ~1,4%
Trinity/Fat Man implozyjny Pu-239 18–21 kt ~6 kg Pu ~17–20%
Mark 3 (późny) implozyjny Pu-239 21–23 kt 6 kg Pu ~20–22%
Mike (1952, test) termojądrowy Pu+fuzja D-T 10,4 Mt — wysoka

Tabela pokazuje, że sprawność nie jest stała — zależy od typu inicjacji i od projektowanych warunków nadkrytycznych. W urządzeniach termojądrowych (oba etapy: fisja + fuzja) sprawność liczona per kg Pu jest wielokrotnie wyższa, bo fuzja generuje dodatkową energię bez zużycia paliwa rozszczepialnego.6

Dezinsercja hydrodynamiczna — szczegóły

Dezinsercja hydrodynamiczna rdzenia to kluczowy element ograniczający sprawność i wymagający zrozumienia jako fizyczny mechanizm.

Gdy rozszczepienia generują energię, temperatura i ciśnienie w rdzeniu dramatycznie wzrastają. Na początku reakcji rdzeń jest w stanie stałym lub ciekłym (dla Pu-239 temperatura topnienia to ~640°C). Po emisji kilku pokoleń neutronów temperatura sięga setek tysięcy kelwinów. Generuje to ciśnienie rzędu gigapaskali — milionów atmosfer.

Równocześnie rdzeń nie jest sztywny — fala ciśnienia biegnie od środka do zewnątrz z prędkością dźwięku $v_s$. Gdy dociera do zewnętrznej powierzchni rdzenia, powoduje wyrzucenie materiału na zewnątrz (faza ablacyjna) i wytworzenie fali rozrzedzenia biegnącej do środka. Ta fala rozrzedzenia obniża ciśnienie i gęstość wewnątrz rdzenia, zmniejszając $k$ poniżej jedności — rdzeń staje się podkrytyczny.

Czas, w którym fala rozrzedzenia dociera do centrum rdzenia (pół czasu przejścia przez rdzeń), to:

$$t_{dis} \approx \frac{R}{v_s} \approx \frac{R}{3{-}4\,\mathrm{km/s}}$$

Dla $R = 5\,\mathrm{cm}$ (kula Pu-239 w warunkach normalnych) $t_{dis} \approx 1{,}5{-}2\,\mu\mathrm{s}$. W tym czasie, przy $\alpha = 10^7\,\mathrm{s^{-1}}$, liczba rozszczepień wzrasta $e^{\alpha t_{dis}} \approx e^{15} \approx 3 \times 10^6$-krotnie od wartości inicjalnej. Im większe $\alpha$, tym więcej wzrostu mieści się w tym samym oknie czasowym.

Tamper spełnia dwie role w tym kontekście. Po pierwsze, fizycznie blokuje rozprężanie zewnętrznej warstwy rdzenia — opóźnia fazę ablacyjną. Po drugie, jako gęsty materiał (wolfram, U-238) ma mniejszą prędkość fali rozrzedzenia, co wydłuża czas $t_{dis}$. Efekt ten w litaraturze opisywany jest jako „inertia confinement" — uwięzienie inercyjne. Każda dodatkowa nanosekundowa utrzymania rdzenia w konfiguracji nadkrytycznej daje wykładniczy wzrost liczby rozszczepień.6

Rola odkryć prolif­eracyjnych i eksportowych

Równanie Bethego-Feynmana i powiązane modele sprawności były przez wiele dekad objęte absolutną tajemnicą. Rząd USA sklasyfikował nie tylko szczegóły, ale nawet ogólne formy równań. Dopiero odtajnienie lat dziewięćdziesiątych (wynikające m.in. z inicjatywy ministra energii Hazel O'Leary w 1993–1994) ujawniło część tych materiałów.

Dla nieproliferacji kluczowe jest, że dysponowanie formą równania sprawności — nawet bez dokładnych danych materiałowych — pozwala na jakościową analizę, jakie parametry musi osiągnąć konstruktor broni. W połączeniu z innymi dostępnymi materiałami (dane neutronowe z ENDF/B, opisy metody implozyjnej, historia prób jądrowych) równanie staje się częścią wiedzy, która realnie opisuje, co jest technicznie trudne.

Agencja IAEA (International Atomic Energy Agency) utrzymuje specjalne grupy ds. proliferacji, które analizują, jak daleko potencjalny aktor niepaństwowy mógłby dojść bez dostępu do tajnych danych. Wyniki takich analiz (szczegóły niejawne) wskazują, że ogólne równania sprawności są możliwe do odtworzenia z literatury otwartej przez kompetentny zespół fizyków. Barierą technologiczną nie jest więc zrozumienie fizyki, lecz dostęp do materiałów rozszczepialnych i zdolności produkcyjne.

Aspekt proliferacyjny równania Bethego-Feynmana pokazuje, że fizyka nie jest neutralna — zrozumienie związku między $\alpha$, $\delta$ i sprawnością wybuchu ma konsekwencje dla kontroli zbrojeń, oceny zdolności nuklearnych i polityki eksportowej komponentów podwójnego zastosowania.6

Nowoczesne podejścia obliczeniowe

Równanie Bethego-Feynmana jest modelem analitycznym, przydatnym dydaktycznie i do szybkich szacunków, ale nowoczesne obliczenia sprawności ładunków jądrowych korzystają ze znacznie bardziej zaawansowanych narzędzi numerycznych.

Hydrodynamika wielodymensjonalna (2D/3D): kody hydrodynamiczne (jak RAGE, OVERTURE, AUTODYN — cześć z nich jawna, część tajna) obliczają ewolucję pola prędkości, gęstości, ciśnienia i temperatury we wszystkich wymiarach przestrzeni jednocześnie. Pozwala to na analizę asymetrii implozji, efektów niestabilności Rayleigha-Taylora na granicy faz, kompresji rdzenia.

Transport neutronów: Kody MCNP (Monte Carlo N-Particle), PARTISN lub DENOVO w trybie sprzężonym z hydrodynamiką obliczają transport neutronów przez zmieniający się materiał rdzenia w każdym kroku czasowym. Wydaje się to oczywiste — ale wymaga ogromu mocy obliczeniowej.

Równania stanu (EOS — Equation of State): zachowanie materii przy temperaturach milionów kelwinów i ciśnieniach gigapaskali jest opisane przez skomplikowane tabelaryczne dane, uzyskane m.in. z eksperymentów na Diamond Anvil Cell, laserowych plazmach i testach NIF (National Ignition Facility). Błędy w EOS przekładają się bezpośrednio na błędy w obliczeniach sprawności.

Sprzężenie promieniowania: Transport fotonów gamma i promieniowania X pomiędzy rdzeniem, tamperem i materiałami otaczającymi jest istotny dla termodynamiki wybuchu i stanowi osobny wielki problem numeryczny.

Całość bywa opisywana jako Multiphysics codes — sprzężone obliczenia neutronów, hydrodynamiki, promieniowania i termomechaniki. Laboratoria LANL (Los Alamos), LLNL (Lawrence Livermore) i SANL (Sandia) utrzymują kod Hydra, kody z rodziny SCF i inne niejawne systemy. Odpowiednikami w Rosji są kody z RFNC-VNIIEF (Sarow) i RFNC-VNIITF (Snieżyńsk).6

Testowanie bez testowania — synergia modeli

Jedną z konsekwencji Traktatu CTBT (Comprehensive Nuclear-Test-Ban Treaty, 1996) jest to, że mocarstwa jądrowe muszą utrzymywać wiarygodność arsenałów bez przeprowadzania testów jądrowych. USA realizuje to przez program Stockpile Stewardship (nadzorowany przez NNSA — National Nuclear Security Administration), który łączy:

  1. Eksperymenty laboratoryjne (w szczególności skrajne ciśnienia i temperatury na NIF, OMEGA laser)
  2. Eksperymenty subkrytyczne (tzw. subcritical experiments, jak seria Reston/Cygnus w Nevada Test Site podziemnie)
  3. Obliczenia numeryczne na superkomputerach (Summit, Sierra, Frontier w ORNL/LLNL — exascale)
  4. Archiwalne dane z testów jądrowych (dane radiochemiczne, hydrody­namiczne, fotometryczne z ponad 1000 prób przeprowadzonych przez USA przed 1992 r.)

Równanie Bethego-Feynmana w tym kontekście jest pierwszym warstwą modelu: szybkim narzędziem, które daje odpowiedź na pytanie „czy wynik wygląda sensownie?". Gdy kodowa symulacja hydrodynamiczna daje wynik, który znacznie odbiega od analitycznego szacunku opartego na równaniu BF, jest to sygnał do sprawdzenia, czy nie ma błędu w danych wejściowych lub w kodzie.6

One-point safety a równanie sprawności

Koncepcja one-point safety (bezpieczeństwa jednopunktowego) jest bezpośrednio związana z równaniem Bethego-Feynmana, choć drogą negatywną: chodzi o to, by sprawność była minimalna w razie przypadkowego zapalenia pojedynczego detonatora implozji.

W układzie implozyjnym detonacja wszystkich soczewek wybuchowych jednocześnie daje symetryczną falę uderzeniową i maksymalną kompresję — a więc maksymalną sprawność. Jeśli jednak odpali się tylko jeden detonator (np. wskutek awaryjnego sygnału), fala uderzeniowa jest asymetryczna. Asymetria ta powoduje, że implozja nie jest sferyczna — część paliwa wylatuje na zewnątrz bez osiągnięcia warunków kritycznych, a osiągana gęstość i $\alpha$ są dużo niższe niż projektowe.

Przepisy bezpieczeństwa USA (w szczególności PAL — Permissive Action Links i AQ&T) wymagają, by w scenariuszu one-point safety uzysk energetyczny nie przekraczał 4 funtów trotylu (ok. 1,8 kg TNT ekwiwalentu). To standardowe wymaganie wszystkich urządzeń nuklearnych w arsenale USA od lat siedemdziesiątych.

Obliczenia one-point safety korzystają wprost ze struktury równania Bethego-Feynmana: przy znacznie mniejszym $\alpha$ i braku symetrycznej implozji sprawność maleje gwałtownie — zgodnie z zależnością $\eta \propto \alpha^2$. Nawet jeśli część materiału rozszczepialnego osiągnie lokalnie stan nadkrytyczny, dezinsercja nastąpi zbyt szybko i przy za małym zapasie nadkrytyczności, by uzysk był istotny.6

Związek z koncepcją boosted fission

Jednym z kierunków zwiększenia sprawności ładunków rozszczepialnych jest tzw. boosted fission — wzmocnione rozszczepienie. Polega ono na wprowadzeniu do rdzenia plutonu lub uranu mieszaniny deuteru i trytu (D-T), które pod wpływem kompresji i temperatury wchodzą w reakcję fuzji termojądrowej, generując dodatkowe neutrony prędkie.

Reakcja D-T:

$${}^2\mathrm{H} + {}^3\mathrm{H} \rightarrow {}^4\mathrm{He} + n\quad (17{,}59\,\mathrm{MeV})$$

Każde zdarzenie fuzji generuje jeden neutron o energii 14,1 MeV — to neutron o energii znacznie wyższej niż typowe neutrony rozszczepienia (~2 MeV). Neutrony fuzji mają większy przekrój czynny na wychwyt przez U-235 lub Pu-239 w zakresie prędkim, co zwiększa $k$ i zatem $\alpha$.

Kluczowy zysk boosted fission z punktu widzenia równania Bethego-Feynmana: wzrost $\alpha$ nawet o 20–30% działa na sprawność quadratycznie ($\eta \propto \alpha^2$), więc 25% wzrost $\alpha$ daje ok. 56% wzrost sprawności. W praktyce boosted fission pozwala na zwielokrotnienie uzysku przy tej samej masie materiału rozszczepialnego lub — co ważniejsze z wojskowego punktu widzenia — zmniejszenie masy rdzenia przy tym samym uzysku.

Technologia ta jest stosowana od lat pięćdziesiątych przez USA (pierwsze boosted urządzenie: Ivy King 1952) i była podstawą miniaturyzacji głowic do rakiet balistycznych ICBM. Rosja, Wielka Brytania, Francja i Chiny posiadają podobne zdolności. Deuter i tryt (D-T) jako booster są przechowywane w osobnych rezerwuarach — T ulega rozpadowi promieniotwórczemu ($t_{1/2} = 12{,}3\,\mathrm{lat}$) i wymaga regularnej wymiany. Stąd konieczność utrzymywania zdolności produkcji trytu przez mocarstwa jądrowe.

Boosted fission nie jest termojądrową bombą H — w tej ostatniej fuzja jest zasadniczym mechanizmem energetycznym. W boosted fission fuzja służy wyłącznie jako source neutronów zwiększający $\alpha$ w rdzeniu rozszczepialnym.6

Znaczenie neutronów inicjujących — inisjator

Równanie sprawności zakłada, że reakcja łańcuchowa startuje od odpowiedniej liczby neutronów inicjujących. To nie jest trywialne — naturalny strumień neutronów (z promieniowania kosmicznego, ze spontanicznego rozszczepienia materiału rozszczepialnego) jest za mały, by zapewnić pewność inicjacji w odpowiedniej chwili, i za nieregularny, by unikać preinicjacji.

Dlatego urządzenia jądrowe zawierają inicjator (initiator, neutron initiator) — źródło neutronów, które wyzwala wymagany neutronowy „strzał inicjujący" w precyzyjnie określonym momencie względem detonacji materiałów wybuchowych.

Wczesnym inicjatorem był Urchin — urządzenie z berylem i polonem-210, wbudowane w centrum rdzenia Fat Man. Gdy fala implozji ściska centrum, Pu-239 mieszał się z Po-210 i Be, generując neutrony przez reakcję ($\alpha$,n). Mechanizm był prosty i niezawodny, ale miał poważną wadę: Po-210 ma czas połowicznego zaniku tylko 138 dni — co oznacza, że inicjatory musiały być regularnie wymieniane.

Nowocześniejsze inicjatory korzystają z elektrycznych generatorów neutronów (Deuterium-Tritium pulse tube generators) lub z zewnętrznych impulsowych źródeł neutronów. Zapewniają one powtarzalny „strzał" w kontrolowany sposób, eliminując zmienność wynikającą z rozpraszania spontanicznego.

Z punktu widzenia równania Bethego-Feynmana inicjator wyznacza $N_0$ — wyjściową liczbę neutronów w chwili $t=0$. Choć sprawność $\eta$ nie zależy logarytmicznie od $N_0$ (bo wykładniczy wzrost szybko niweluje różnice w punkcie startowym), to zbyt mała lub zbyt duża $N_0$ w złym momencie (predetonacja lub późna inicjacja) ma kolosalne konsekwencje dla uzysku.6

Predetonacja — najdroższy błąd projektowy

Predetonacja (lub fizzle) jest najważniejszym ryzkiem obniżającym sprawność bomby i jest bezpośrednio opisana przez równanie Bethego-Feynmana w aspekcie minimizacji $\alpha$ w chwili dezinsercji.

Pełna analiza predetonacji jest zawarta w odrębnym artykule Predetonacja (Fizzle), ale z perspektywy równania BF mechanizm jest jasny:

Jeśli reakcja łańcuchowa startuje za wcześnie (np. wskutek neutronu ze spontanicznego rozszczepienia Pu-240), rdzeń zaczyna generować energię, zanim implozja go wystarczająco skompresuje. Temperatura i ciśnienie narastają, ale $\alpha$ — jeszcze nie na maksymalnej wartości — jest za małe. Rdzeń dezintegruje się przy niskim $\alpha$:

$$\eta_{fizzle} \propto \alpha_{early}^2 \ll \alpha_{max}^2 = \eta_{full}$$

Wynik: uzysk energetyczny jest nieporównanie mniejszy niż projektowy — typowo kilkaset do kilku tysięcy razy mniejszy. Bomba „buch­nie", ale nie wybuchnie.

Właśnie dlatego Pu-239 o wysokiej zawartości Pu-240 (>6–7%) jest nieodpowiedni dla bomby — reaktorowy pluton zawiera 20–25% Pu-240 i jego spontaniczne rozszczepienie jest za częste. W tym kontekście bomba implozyjna była odpowiedzią na problem Pu-240: szybka implozja skraca czas między startem inicjatora a maksymalną kompresją poniżej prawdopodobnego czasu między spontanicznymi rozszczepieniam Pu-240.

Obliczenia probabilistyczne predetonacji (ile neutronów z Pu-240 jest emitowanych na sekundę, jaka jest szansa na uruchomienie reakcji łańcuchowej zanim kompresja dotrze do maksimum) były jednym z pierwszych zastosowań rachunku prawdopodobieństwa w projektowaniu broni jądrowej — Feynman i inni w Los Alamos. Wyniki prowadziły do wniosku, że bomb działowa (gun-type) z plutonem jest całkowicie nierealizowalna.6

Ograniczenia modelu i jego wartość

Równanie Bethego-Feynmana w formie $\eta = KR^2\alpha^2\delta$ jest z definicji uproszczeniem. Warto wymienić wyraźnie jego ograniczenia:

  1. Jednorodna ekspansja: Model zakłada, że rdzeń rozsadzany jest jednorodnie we wszystkich kierunkach. W rzeczywistości implozja jest nigdy idealnie symetryczna, generując niestabilności hydrodynamiczne.

  2. Stała gęstość podczas deinsercji: Gęstość rdzenia zmienia się drastycznie w trakcie wybuchu. Model ignoruje tę dynamikę, co prowadzi do błędów rzędu czynnika 2–5 w obliczeniach sprawności.

  3. Pomijanie fotonów gamma: Promieniowanie gamma z rozszczepionych jąder ogrzewa materiał rdzenia przez absorpcję, co modyfikuje EOS i prędkość dezinsercji. Model analityczny tego nie uwzględnia.

  4. Brak promieniowania X: Przy temperaturach ~10⁷ K dominującą formą transferu energii jest promieniowanie X (Planck), a nie konwencjonalne przewodzenie cieplne. To zmienia zarówno dynamikę ciśnienia, jak i efektywną prędkość dezinsercji.

  5. Wyciek neutronów niekonstantny: Model Bethe-Feynmana traktuje wyciek jako prostą funkcję geometrii, podczas gdy w trakcie dezinsercji geometria i skład izotopowy rdzenia zmieniają się dynamicznie.

Mimo tych ograniczeń model ma wielką wartość dydaktyczną: wyraża najistotniejsze zależności w zwartej formie analitycznej, pozwala na szybkie szacunki parametryczne i nadaje się do ćwiczeń studentów fizyki reaktorowej i broni jądrowej. To on nadał strukturę konceptualną całemu dalszemu rozwojowi bardziej dokładnych modeli numerycznych.1,2,6

Równania sprawności w literaturze otwartej

Kilka form równania sprawności wybuchu rozszczepialnego dostępnych jest w literaturze otwartej, co pozwala na ich porównanie i weryfikację intuicji Bethego-Feynmana:

Formuła Glasstone'a-Dolana (z The Effects of Nuclear Weapons, wyd. USGPO):

Podręcznik Glasstone'a i Dolana (wielokrotnie aktualizowany, ostatnie wydanie 1977) zawiera przystępne opisy mechanizmów wybuchu jądrowego, w tym pewne elementy analizy sprawności. Formuła tam podana w termach jakościowych jest spójna z modelem BF: sprawność rośnie ze wzrostem nadkrytyczności i maleje ze wzrostem tempa ekspansji.

Analiza Rhodes'a (z The Making of the Atomic Bomb, 1986):

Richard Rhodes w swojej nagrodzonej Pulitzerem historii Projektu Manhattan opisuje problem sprawności w sposób nieformalny, ale fizycznie trafny: kluczem jest wyścig między kolejnymi pokoleniami neutronów a hydrodynamicznym rozpadaniem rdzenia. Rhodes opierał się na rozmowach z byłymi fizykami Los Alamos i jest jednym z najlepszych dostępnych publicznie opisów.

Opracowanie Carsona Marka (John Carson Mark, były kierownik T-Division LANL, artykuły z lat 70–80. po częściowym odtajnieniu):

Carsona Mark wskazał publicznie, że ogólna postać zależności sprawności od nadkrytyczności i tempa mnożenia jest w istocie dostępna z zasad fizyki podstawowej. Modele szczegółowe (EOS, transport) wymagają tajnych danych, ale ogólna struktura równania BF jest odtwarzalna.

Te odniesienia świadczą o tym, że wiedza o równaniu Bethego-Feynmana nie jest ściśle tajna — tajna jest jej dokładna postać z konkretnymi parametrami dla konkretnych materiałów i geometrii. Ogólna forma i fizyczna interpretacja są dostępne w literaturze akademickiej. Dostęp do takich opracowań jak The Los Alamos Primer (Serber, odtajniony 1992) czy Los Alamos Science (seria specjalna) potwierdza, że fundamenty fizyczne można rekonstruować ze źródeł publicznych — choć pełne obliczenia projektowe wymagają zaawansowanych kodów numerycznych i tajnych danych o właściwościach materiałów w ekstremalnych warunkach.6

Sprawność a kaliber — reaktory, bomby, impulsory

Interesującym dopełnieniem są urządzenia impulse reactors (impulsowe reaktory), działające na zasadzie krótkich, kontrolowanych stanów nadkrytycznych — a więc w pewnym sensie na granicy między reaktorem a bombą. Przykłady to:

  • TRIGA (Training, Research, Isotopes, General Atomics): reaktory badawcze, które mogą być pulsowane do momentowej mocy tysięcy megawatów przy czasie impulsu ~20 ms, zanim ujemny współczynnik temperatury przywróci układ do podkrytyczności.
  • Godiva (seria urządzeń krytycznych w Los Alamos): gołe kule uranu-235 lub plutoniu, używane do pomiarów krytycznych i pulsów neutronowych.

W TRIGA sprawność wybuchu energetycznego (udział paliwa rozszczepiający się w impulsie) to ułamek procenta — bo ujemny współczynnik temperatury działa jak efektywna dezinsercja. W „gotowej" bombie mechanizm dezinsercji jest hydrodynamiczny, nie termiczny — rdzeń ulega dezintegracji mechanicznie, nim zdąży powoli ochłodzieć cieplnie.

To porównanie pokazuje, że dezinsercja (niezależnie od mechanizmu) jest kluczową zmienną równania Bethego-Feynmana. W reaktorze impulsowym dezinsercja jest pożądana i kontrolowana. W bombie jest nieunikniona i projektuje się cały układ tak, by nastąpiła jak najpóźniej.6

Kontekst akademicki i dydaktyczny

Równanie Bethego-Feynmana jest wykładane na kierunkach fizyki reaktorowej i nauki o broni jądrowej na uczelniach takich jak MIT (kurs 22.101 Nuclear Physics i Radiation), Naval Postgraduate School w Monterey, Texas A&M (NUEN) czy w europejskich programach ESARDA (European Safeguards Research and Development Association).

W Polsce wykłady z fizyki jądrowej na Wydziale Fizyki UW, AGH i Politechniki Warszawskiej zawierają elementy kinetyki neutronowej i sprawności wybuchu, choć z zrozumiałych względów nie na poziomie szczegółowości programów wojskowych. Podręczniki Słowińskiego (Podstawy fizyczne energetyki jądrowej) i Gdenarza (Jak się bać broni jądrowej) są polskimi zasobami na poziomie akademickim, które omawiają te tematy w granicach informacji publicznej.

Równanie BF w kontekście dydaktycznym służy przede wszystkim do budowania intuicji: dlaczego implozja jest lepsza od metody działowej, dlaczego masa krytyczna zależy od gęstości, dlaczego tamper zwiększa sprawność dwoma niezależnymi mechanizmami. Dla studentów fizyki jądrowej jest naturalnym pomostem między teorią transportu neutronów a rozumieniem podstawowej zasady działania broni jądrowej — bez omawiania tajnych szczegółów konstrukcyjnych.

Historyczną ciekawostką jest, że Feynman — po wojnie — był jednym z najbardziej otwartych krytyków tajności otaczającej fizykę jądrową. Wskazywał, że część tej tajności jest iluzoryczna: ogólna fizyka wybuchu jest dostępna z zasad pierwszych. To, co pozostaje naprawdę tajne i co jest prawdziwą barierą proliferacyjną, to nie teoria, lecz praktyczne dane o zachowaniu materiałów, wyniki pomiarów przekrojów czynnych w zakresie prędkim, szczegóły inżynierii detonatorów i opis kalibracji układów zapłonowych. Teoria jest w podręcznikach; bariery są gdzie indziej.1,2,6

Dodatkowe materiały multimedialne

Do tego artykułu nie dodano jeszcze materiałów wideo. Warto wrócić do tej sekcji dopiero wtedy, gdy uda się znaleźć materiał, który rzetelnie pokazuje związek między @@MATH_INLINE_83@@, kompresją rdzenia i dezinsercją hydrodynamiczną, zamiast sprowadzać temat do samego hasła „więcej plutonu daje większy wybuch”.

Jako rozwinięcie warto czytać ten tekst razem z masą krytyczną, predetonacją i metodą implozyjną. Wtedy dobrze widać, że równanie nie jest abstrakcyjną ciekawostką, lecz skrótem opisującym to, o co naprawdę walczy konstruktor: czas, nadkrytyczność i wykorzystanie materiału zanim układ zdąży się rozpaść.

Powiązane kalkulatory i narzędzia

  • k_eff — pokazuje, jak geometria, moderator i straty neutronów wpływają na krytyczność układu.

Ćwiczenia praktyczne

Ćwiczenie laboratoryjne powinno polegać na zbudowaniu prostego modelu obliczeniowego wydajności opartego na zależności typu $\eta = K R^2 \alpha^2 \delta$. W wariancie podstawowym należy:

  1. przyjąć realistyczne, ale umowne wartości $K$, $R$, $\alpha$ i $\delta$,
  2. policzyć względną zmianę wydajności przy zmianie każdego z parametrów osobno,
  3. sprawdzić, jak silnie wynik reaguje na niewielki spadek $\alpha$,
  4. porównać wariant z większym promieniem i gorszym mnożeniem z wariantem mniejszym, lecz mocniej skompresowanym,
  5. opisać, które parametry są najbardziej „dźwigniowe” z punktu widzenia projektu.

Celem ćwiczenia nie jest przewidywanie prawdziwej mocy konkretnej broni, lecz uchwycenie wrażliwości sprawności na parametry neutroniczne i geometryczne. W drugiej części należy pokazać, jak wynik zmienia się przy sztucznym obniżeniu $\alpha$, symulując w ten sposób przedwczesną inicjację.

Drugie ćwiczenie, teoretyczno-przemysłowe, powinno dotyczyć kompromisów konstrukcyjnych. Należy:

  1. wyjaśnić, dlaczego zwiększenie masy materiału bez poprawy kompresji nie daje proporcjonalnego wzrostu uzysku,
  2. opisać wpływ tampra na neutronikę i hydrodynamikę,
  3. porównać rolę szybkiej implozji i czystości materiału rozszczepialnego,
  4. wskazać, które parametry najłatwiej zepsuć przez asymetrię układu albo predetonację,
  5. powiązać wynik z artykułami o predetonacji, masie krytycznej i reakcji łańcuchowej.

To ćwiczenie ma pokazać, że równanie sprawności nie jest ozdobnikiem teoretycznym. Jest skrótem opisującym, dlaczego niemal każdy element konstrukcji ładunku rozszczepieniowego podporządkowano walce o kilka dodatkowych pokoleń rozszczepień przed rozpadem rdzenia.

Przejdź do ćwiczenia interaktywnego

Powiązane artykuły