Spotlight #29 – Fizyka Plazmy: Ekranowanie Debye'a, Zorza Polarna, Uwięzienie Fuzyjne i Magnetohydrodynamika

Plazma to gaz zjonizowany: zbiór swobodnych elektronów i dodatnich jonów, w którym dominującą rolę odgrywają siły elektromagnetyczne dalekiego zasięgu, a nie zderzenia bliskiego zasięgu. Plazma istnieje wszędzie tam, gdzie energia termiczna lub elektromagnetyczna jest wystarczająca do zjonizowania atomów — we wnętrzach gwiazd przy 10 milionach stopni, w lampach fluorescencyjnych przy 20 000 K, w kanałach błyskawic przy 30 000 K i w tokamakach fuzyjnych przy 150 milionach stopni.

Kluczową właściwością odróżniającą plazmę od losowej mieszaniny jonów i elektronów jest jej zachowanie zbiorowe: naładowane cząstki reagują na pola sąsiadów na odległościach znacznie większych niż ich rozdzielenie, tworząc fale, niestabilności oraz zdolność do uwięzienia prądów i pól magnetycznych. Zrozumienie zachowania plazmy wymaga zarówno fizyki pojedynczej cząstki (ruch cyklotronowy, dryfy), jak i teorii płynów (magnetohydrodynamika, MHD).

1. Ekranowanie Debye'a — Długość Samo-Ekranowania Plazmy

Umieść dodatni ładunek testowy w plazmie. Pobliskie elektrony są przyciągane, a jony odpychane; powstała chmura ładunku częściowo ekranuje gołe pole kulombowskie. Charakterystyczna odległość, na której ekranowanie jest skuteczne, to długość Debye'a λ_D. Poza tą odległością plazma wydaje się elektrycznie neutralna; wewnątrz niej pojedyncze cząstki doświadczają znaczących pól elektrycznych. Długość Debye'a to najbardziej fundamentalna skala długości każdej plazmy.

Potencjał elektryczny wokół ładunku testowego q w plazmie:
  φ(r) = (q / 4πε_0 r) · exp(−r/λ_D)   (potencjał Yukawy / ekranowany Coulomba)
  λ_D = długość Debye'a

Długość Debye'a (dominuje ekranowanie elektronowe):
  λ_D = √(ε_0 k_B T_e / n_e e²)
  T_e = temperatura elektronów, n_e = gęstość liczbowa elektronów

Wartości liczbowe:
  Korona słoneczna:  T_e ~ 10^6 K,  n_e ~ 10^12 m^−3  → λ_D ~ 0,07 m
  Jonosfera:    T_e ~ 1500 K,  n_e ~ 10^11 m^−3  → λ_D ~ 0,003 m
  Tokamak:       T_e ~ 10^8 K,  n_e ~ 10^20 m^−3  → λ_D ~ 7 × 10^−5 m
  Lampa fluorescencyjna: T_e ~ 2 eV, n_e ~ 10^18 m^−3  → λ_D ~ 10^−5 m

Kryteria plazmy (warunki dla zachowania zbiorowego):
  1. Quasi-neutralność: rozmiar systemu L ≫ λ_D   (masa jest neutralna; tylko warstwy graniczne odbiegają)
  2. Wiele cząstek na sferę Debye'a: N_D = (4/3)πλ_D³ n ≫ 1
     N_D ≫ 1 zapewnia ważność założenia rozkładu ciągłego.
  3. Częstość plazmowa ω_p ≫ częstość zderzeń ν_zderzeń
     ω_p = √(n_e e² / ε_0 m_e)   (częstość oscylacji plazmowych)
     Przy f < f_p = ω_p/2π fale elektromagnetyczne nie mogą się propagować (odbicie).

Warstwa graniczna plazmy (sheath):
  Na granicy quasi-neutralność załamuje się na dystansie ~kilku λ_D.
  Warstwa graniczna to cienka warstwa dodatniego ładunku przestrzennego, zapewniająca brak przepływu prądu netto.
  Fizyka warstwy granicznej rządzi interakcjami plazma-materiał zarówno w produkcji chipów, jak i w tokamakach fuzyjnych.

2. Ruch Naładowanej Cząstki — Spirale, Dryfy i Lustra

Naładowana cząstka w jednorodnym polu magnetycznym porusza się po okręgu: siła Lorentza qv×B działa prostopadle do prędkości, dając przyspieszenie dośrodkowe. W połączeniu z dowolną składową prędkości wzdłuż pola, ruch jest helisą. Pola niejednorodne generują systematyczne dryfy prostopadłe zarówno do pola, jak i jego gradientu. Te dryfy pojedynczej cząstki leżą u podstaw wielkoskalowych przepływów plazmy w magnetosferze Ziemi i urządzeniach fuzyjnych.

Ruch cyklotronowy (żyracyjny):
  Równanie: m dv/dt = qv × B
  Ruch kołowy prostopadły do B z:
    Częstość cyklotronowa: ω_c = |q|B/m   (częstość kątowa)
    Promień Larmora:       r_L = mv_⊥ / |q|B = v_⊥ / ω_c
  W połączeniu z v_∥ wzdłuż B: trajektoria helikalna

Promienie Larmora (w polu 1 T):
  Elektron  (m = 9,1×10^−31 kg, T_e = 100 eV → v_⊥ = 1,87×10^8 m/s):  r_L ~ 1 mm
  Proton    (m = 1,67×10^−27 kg, T_i = 100 eV → v_⊥ = 1,38×10^5 m/s):  r_L ~ 1,4 mm
  Jon alfa (m = 6,64×10^−27 kg, T = 3,5 MeV → v_⊥ = 1,3×10^7 m/s):   r_L ~ 20 cm

Dryf E×B (w połączonych polach elektrycznym i magnetycznym):
  v_{E×B} = (E×B) / B²   (niezależny od ładunku i masy cząstki!)
  Odpowiedzialny za ogólny przepływ masy plazmy; zarówno jony, jak i elektrony dryfują razem.

Dryf gradientowy (w niejednorodnym B):
  v_{∇B} = (1/2) v_⊥ r_L (∇B × B) / B²   (zależny od ładunku)
  Kierunek: przeciwny dla jonów i elektronów → tworzy prąd netto.

Dryf krzywiznowy:
  W zakrzywionym polu magnetycznym (promień krzywizny R_c):
  v_c = (mv_∥² r_L) / (qB R_c²) (&hat;R_c × B)
  Dryfy krzywiznowe + gradientowe napędzają prąd pierścieniowy w pasach Van Allena Ziemi.

Efekt lustra magnetycznego:
  Naładowana cząstka poruszająca się z pola słabego do silnego jest odbijana, jeśli kąt padania α < α_lustra:
  sin²α_lustra = B_min / B_max   (stosunek lustrzany)
  Pasy radiacyjne Van Allena to butelki magnetyczne: cząstki odbijają się między punktami lustrzanymi.

3. Zorza Polarna — Naładowane Cząstki Malujące Niebo

Zorza polarna jest bezpośrednią wizualną konsekwencją fizyki pojedynczej cząstki plazmy. Elektrony i protony wiatru słonecznego są prowadzone wzdłuż dipolowych linii pola magnetycznego Ziemi, spiralnie zbliżając się do regionów polarnych. Tam, gdzie linie pola nurkują w górną atmosferę, energetyczne cząstki zderzają się z cząsteczkami tlenu i azotu, wzbudzając je do wyższych stanów energetycznych. Emisja światła podczas rozpadu tych stanów tworzy charakterystyczne kolory zorzy.

Parametry wiatru słonecznego (typowe):
  v_sw ~ 400–700 km/s (prędkość)
  n_sw ~ 5–15 cm^−3 (gęstość protonów przy 1 AU)
  B_IMF ~ 5–10 nT (międzyplanetarne pole magnetyczne)
  Energia protonu: ½m_p v² ~ 1–3 keV

Warstwa prądowa Chapmana-Ferraro:
  Wiatr słoneczny jest odchylany przez dipolowe pole Ziemi na magnetopauzie (~10 R_E w stronę Słońca).
  Rekoneksja magnetyczna na dziennej magnetopauzie wstrzykuje cząstki do magnetosfery.
  Rekoneksja zachodzi też w ogonie magnetosfery, przyspieszając cząstki w stronę Ziemi.

Wysokości opadania cząstek i kolory:
  O (tlen atomowy) linia zielona 557,7 nm:   90–150 km wysokości  (najczęstsza)
  O (tlen atomowy) linia czerwona  630,0 nm:   200–300 km wysokości  (miękkie elektrony, <1 keV)
  N&sub2;+ niebiesko/fioletowa 391,4 nm:             60–100 km wysokości  (twarde elektrony, >50 keV)
  Wysokość emisji koreluje z energią cząstek: twardsze cząstki penetrują głębiej.

Indeks Kp i owal zorzowy:
  Kp 0–2: zorza widoczna tylko powyżej 65° szerokości geomagnetycznej
  Kp 5–6: widoczna do ~55° (Szkocja, Kanada, Niemcy)
  Kp 9 (burza): widoczna do <40° (południowe USA, Europa Środkowa, Japonia)

Cykl podburzy:
  1. Faza wzrostu (~1 h): strumień magnetyczny gromadzi się w ogonie magnetosfery
  2. Początek: eksplozywna rekoneksja w bliskim ogonie
  3. Ekspansja: jasny impuls rozprzestrzenia się na zachód wzdłuż owalu zorzowego
  4. Powrót: stopniowe rozpraszanie; pozostają spokojne łuki
  Energia wprowadzona podczas dużej burzy: ~ 10^17 J (porównywalna z tysiącami broni jądrowej)

4. Fuzja w Tokamaku — Uwięzienie Gwiazdy w Butelce

Fuzja jądrowa — źródło mocy gwiazd — wymaga temperatur plazmy rzędu 100–200 milionów stopni Celsjusza. W tych temperaturach jądra deuteru i trytu pokonują swoje odpychanie elektrostatyczne i się łączą, uwalniając hel plus energetyczny neutron. Wyzwaniem jest utrzymanie tak ekstremalnej plazmy wystarczająco długo i gorąco, by wynik energetyczny przewyższył energię wejściową. Tokamak, toroidalne urządzenie uwięzienia magnetycznego, jest obecnie wiodącym podejściem.

Reakcja fuzji D–T:
  ₁²H + ₁³H → ₂&sup4;He (3,52 MeV) + n (14,06 MeV)
  Razem: 17,59 MeV na reakcję (wobec ~200 MeV dla rozszczepienia U-235)
  Szczyt przekroju czynnego: ~5 barn przy E_cm ~ 100 keV (T ~ 10^8 K)

Kryterium Lawsona (warunek zapłonu / break-even):
  n τ_E T ≥ n τ_E T|_zapłon ~ 3 × 10^21 m^−3 s keV  (dla D–T)
  n = gęstość plazmy (m^−3)
  τ_E = czas uwięzienia energii (s)
  T = temperatura (keV, gdzie 1 keV = 11,6 miliona K)
  Moc fuzji: P_fus ~ n² ⟨σv⟩ E_fus / 4  (dla mieszaniny D–T 50–50)

Rekord świata JET (2022):
  59 MJ energii fuzji w 5 s → 11,8 MW podtrzymywanej mocy
  n τ_E T ~ 1,5 × 10^20 m^−3 s keV  (czynnik ~20 poniżej zapłonu)

ITER (w budowie):
  Promień główny R = 6,2 m, promień mały a = 2 m, objętość plazmy ~840 m³
  Cel: Q = P_fus / P_grzania ≥ 10  (500 MW na wyjściu, 50 MW na wejściu)
  Oczekiwane n τ_E T ~ 6 × 10^21 (blisko zapłonu)

Równowaga MHD (tokamak):
  Równanie Grada-Shafranova (równowaga osiowosymetryczna):
    Δ*ψ = −μ_0 R² dp/dψ − f df/dψ
    ψ = funkcja strumienia poloidalnego
    p = ciśnienie, f = RB_φ (funkcja pola toroidalnego)
  Rozwiązanie daje zagnieżdżone powierzchnie strumieniowe; plazma jest uwięziona na torusach.
  Współczynnik bezpieczeństwa: q = r B_φ / (R B_θ)  (stosunek obrotów linii pola toroidalnego do poloidalnego)
  Stabilność wymaga q > 1, by uniknąć niestabilności kink; typowe q ~ 3–5 na krawędzi plazmy.

Kluczowe niestabilności:
  Niestabilność kink (q < 1): linie pola zaginają się, prowadząc do zakłócenia
  Mod balonowy: plazma wybrzusza się na zewnątrz przy wysokim β (stosunek ciśnienia plazmy do magnetycznego)
  ELM (mody zlokalizowane na krawędzi): okresowe wybuchy na krawędzi plazmy; muszą być kontrolowane, by chronić ścianę.

5. Magnetohydrodynamika — Plazma jako Płyn Przewodzący

Gdy skale długości i czasu, które nas interesują, są znacznie większe niż długość Debye'a i okres cyklotronowy, plazmę można traktować jako pojedynczy płyn przewodzący — magnetohydrodynamikę (MHD). MHD opisuje propagację wiatru słonecznego, generację ziemskiego pola geomagnetycznego, tworzenie strumieni astrofizycznych i wielkoskalową stabilność plazmy fuzyjnej.

Idealne równania MHD (pojedynczy płyn, idealny przewodnik):
  Ciągłość:     ∂ρ/∂t + ∇·(ρv) = 0
  Pęd:       ρ Dv/Dt = −∇P + J×B
  Energia:         D/Dt (P/ρ²) = 0   (adiabatyczna)
  Indukcja:      ∂B/∂t = ∇×(v×B)   (twierdzenie o zamrożonym strumieniu: B porusza się z płynem)
  Prawo Ohma:     E + v×B = ηJ      (η = oporność, η = 0 dla idealnej MHD)
  Domknięcie:        ∇×B = μ_0 J,   ∇·B = 0

Fala Alfvéna (fala poprzeczna wzdłuż B):
  Prędkość fazowa: v_A = B / √(μ_0 ρ)   (prędkość Alfvéna)
  W koronie słonecznej: B ~ 100 G, ρ ~ 10^−12 kg/m³ → v_A ~ 3 × 10^6 m/s ~ 0,01 c
  Fale Alfvéna transportują energię wzdłuż linii pola i przypuszczalnie ogrzewają koronę.

Twierdzenie o zamrożonym strumieniu (idealna MHD):
  Linie pola magnetycznego są „zamrożone" w plazmie: gdy plazma się porusza, linie pola poruszają się razem z nią.
  Konsekwencja: topologia magnetyczna jest zachowana w idealnej MHD.

Rekoneksja magnetyczna (rezystywna MHD):
  Tam, gdzie spotykają się linie pola o przeciwnej polaryzacji, oporność η ≠ 0 pozwala liniom pola
  pęknąć i połączyć się na nowo w nowej topologii, uwalniając zmagazynowaną energię magnetyczną eksplozywnie.
  Tempo Sweeta-Parkera: v_rek ~ v_A / √(S)   (S = liczba Lundquista = μ_0 v_A L / η ~ 10^12 w koronie słonecznej)
  Obserwowane tempo rekoneksji ~ 10^−2 v_A  (znacznie szybsze: problem „szybkiej rekoneksji")
  Napędza rozbłyski słoneczne (do 10^25 J), koronalne wyrzuty masy (CME) i burze geomagnetyczne.

Ciśnienie i naprężenie magnetyczne:
  Ciśnienie magnetyczne: P_B = B² / (2μ_0)  (działa jak ciśnienie gazu, izotropowe)
  Naprężenie magnetyczne:  T_B = B² / μ_0  (działa jak napięta struna wzdłuż linii pola)
  Beta plazmy: β = P_plazma / P_B = n k_B T / (B²/2μ_0)
  Rdzeń tokamaka: β ~ 0,03–0,1;  Korona słoneczna: β ~ 0,01–0,1;  Wiatr słoneczny: β ~ 1

Rekoneksja magnetyczna to jeden z centralnych nierozwiązanych problemów fizyki plazmy. W ramach rezystywnej MHD przewidywane tempo rekoneksji jest zbyt wolne, by wyjaśnić obserwowane rozbłyski słoneczne, które mogą uwolnić swoją energię w ciągu minut. Współczesne teorie — w tym MHD Halla, fizyka skali elektronowej i rekoneksja napędzana turbulencją — mogą dopasować obserwowane szybkie tempa, ale pozostają przedmiotem aktywnych badań eksperymentalnych. Misja Parker Solar Probe, wystrzelona w 2018 roku i obecnie latająca w odległości 9 promieni słonecznych, zbiera pomiary in situ regionów rekoneksji, które bezpośrednio przetestują te teorie.

Plazma: Dominujący Stan Wszechświata

Ponad 99% całej widzialnej materii we wszechświecie znajduje się w stanie plazmy. Każda gwiazda, każde aktywne jądro galaktyczne, każda zjonizowana mgławica i każda cząstka wiatru słonecznego to plazma. A jednak na Ziemi plazma jest laboratoryjną ciekawostką, narzędziem przemysłowym (trawienie chipów, spawanie, neony) lub przejściowym zjawiskiem (błyskawice, płomienie). Symulacje w tej fali śledzą pomost od zjawisk w skali Ziemi — lampy fluorescencyjne i zasłony zorzy — do astrofizycznych procesów plazmowych, które rządzą ewolucją gwiazd, wzmocnieniem pola magnetycznego i budżetem energetycznym gromad galaktyk. Wszystko z tych samych równań: praw Maxwella i drugiej zasady dynamiki Newtona dla naładowanych cząstek.