Spotlight #32 – Astrofizyka i ewolucja gwiazd: diagram HR, struktura gwiazdowa, supernowe i fale grawitacyjne

Każdy atom cięższy od litu powstał wewnątrz gwiazdy, a większość wodoru i helu we wszechświecie wykuto w pierwszych trzech minutach po Wielkim Wybuchu. Astrofizyka gwiazdowa znajduje się na przecięciu fizyki jądrowej, dynamiki płynów, transferu radiacyjnego i ogólnej teorii względności. Najprostsze gwiazdy — obiekty ciągu głównego w równowadze hydrostatycznej — są dobrze opisane przez garstkę równań. Te skomplikowane — gwiazdy asymptotycznej gałęzi olbrzymów z wieloma powłokami spalania, rotujące masywne gwiazdy tracące masę w wiatrach, zwarte układy podwójne spiralnie zbliżające się do zlania — należą do najbardziej aktywnych granic współczesnej fizyki.

1. Diagram Hertzsprunga-Russella

Gdy Ejnar Hertzsprung i Henry Norris Russell niezależnie nanieśli jasność gwiazd względem temperatury powierzchni około 1910 roku, odkryli, że gwiazdy nie są rozłożone losowo: zdecydowana większość skupia się wzdłuż ukośnego pasa (ciąg główny), z gałęziami olbrzymów i nadolbrzymów powyżej i białymi karłami poniżej. Diagram HR jest najbardziej użytecznym pojedynczym diagramem w astrofizyce gwiazdowej, ponieważ każdy etap ewolucyjny zajmuje odrębny region.

Relacja Stefana-Boltzmanna:
  L = 4πR² σ T_eff&sup4;

Relacja masa-jasność ciągu głównego (przybliżona):
  L / L_słońce ≈ (M / M_słońce)^4      dla 0,5 < M/M_słońce < 50
  L / L_słońce ≈ (M / M_słońce)^2,5    dla gwiazd o mniejszej masie

Czas życia na ciągu głównym:
  t_MS ≈ t_słońce · (M/M_słońce)^{-2,5}    (słoneczny t_MS ≈ 10 mld lat)
  1 M_słońce: ~10 mld lat   5 M_słońce: ~100 mln lat   25 M_słońce: ~7 mln lat

Regiony diagramu HR (L w L_słońce, T w K):
  Ciąg główny:  0,001–10^6 L_słońce, 2500–50 000 K
  Czerwone olbrzymy: 10–1000 L_słońce, 3500–5500 K
  Nadolbrzymy:  10^4–10^6 L_słońce, 3500–30 000 K
  Białe karły:  10^-4–0,1 L_słońce, 5000–100 000 K

Typy widmowe (OBAFGKM, od najgorętszych do najchłodniejszych):
  O: >30 000 K   B: 10 000–30 000 K   A: 7500–10 000 K
  F: 6000–7500 K   G: 5200–6000 K   K: 3700–5200 K
  M: <3700 K (chłodny czerwony karzeł; najliczniejszy typ w Galaktyce)

Pozycja gwiazdy na diagramie HR zmienia się w miarę jej ewolucji. Słońce opuści ciąg główny za ~5 mld lat, rozszerzy się w czerwonego olbrzyma (R ~200 R_słońce), utraci zewnętrzne warstwy jako mgławica planetarna i pozostawi po sobie ścieżkę chłodzenia białego karła. Masywne gwiazdy typu O ewoluują tak szybko, że nigdy nie osiągają równowagi z ciśnieniem gazu otoczenia: umierają w eksplozjach supernowych w ciągu 3–10 mln lat od narodzin. Gromady gwiazd są szczególnie użyteczne, ponieważ wszyscy członkowie mają ten sam wiek; punkt zwrotny — najjaśniejsza gwiazda ciągu głównego wciąż spalająca wodór — bezpośrednio podaje wiek gromady.

2. Struktura wnętrza gwiazdy i politropy

Gwiazda to samograwitująca kula gazu w równowadze hydrostatycznej: grawitacja ciągnąca do wewnątrz i gradient ciśnienia pchający na zewnątrz równoważą się na każdej warstwie. Trzy dodatkowe równania opisują generację energii, transport energii i sposób, w jaki temperatura zmienia się wraz z ciśnieniem. Razem tworzą równania struktury gwiazdowej, które jednoznacznie określają wewnętrzny profil gwiazdy, gdy podana jest jej masa i skład (twierdzenie Vogta-Russella).

Równowaga hydrostatyczna:
  dP/dr = −ρ G M(r) / r²

Ciągłość masy:
  dM/dr = 4πr² ρ

Transport energii (strefa radiacyjna):
  dT/dr = −3κρL / (64πσr²T³)
  κ = średnia rosselandowska nieprzezroczystość (cm²/g)

Kryterium konwekcji (Schwarzschild):
  |dT/dr|_rad > |dT/dr|_ad   → niestabilność konwekcyjna
  Teoria drogi mieszania: skala komórki konwekcyjnej ≈ α_MLT H_P  (α ~ 1,5-2)

Politropowe równanie stanu:
  P = K ρ^{(n+1)/n}   (wskaźnik politropowy n)
  n = 3/2 → w pełni konwekcyjna (gwiazdy małomasywne, jądra olbrzymów)
  n = 3   → zdominowana ciśnieniem promieniowania (masywne gwiazdy; Chandrasekhar)

Równanie Lane'a-Emdena (bezwymiarowa struktura gwiazdowa):
  d/dξ[ξ² dθ/dξ] + ξ² θ^n = 0
  ξ = r / r_n (skalowany promień),  θ = (ρ/ρ_c)^{1/n}

Wnętrze Słońca (model AGSS09):
  Jądro (r < 0,25 R_słońce):  T_c = 1,57×10^7 K,  ρ_c = 151 g/cm³
  Strefa radiacyjna: 0,25–0,71 R_słońce
  Strefa konwekcyjna: 0,71–1,00 R_słońce  (głębokość ~200 000 km)

Helioseismologia — pomiar prędkości fal dźwiękowych poprzez obserwację oscylacji powierzchniowych — zweryfikowała strukturę wnętrza Słońca z dokładnością poniżej 1%. Rozbieżność między trójwymiarowymi spektroskopowymi obfitościami słonecznymi (AGSS09) a ograniczeniami helioseismicznymi pozostaje otwartym problemem fizyki gwiazdowej: „problem metaliczności Słońca".

3. Cykle spalania jądrowego

Źródłem energii pozwalającym gwiazdom opierać się grawitacyjnemu zapadaniu jest fuzja jądrowa. Gwiazdy małomasywne (M < 1,3 M_słońce) fuzjują wodór głównie poprzez łańcuch proton-proton (pp); masywniejsze gwiazdy wykorzystują cykl CNO, którego tempo jest wysoce wrażliwe na temperaturę (T^17 wobec T^4 dla pp). W miarę starzenia się gwiazda buduje jądro popiołu helowego, następnie zapala hel poprzez proces potrójnego alfa, a dla masywnych gwiazd przechodzi przez spalanie węgla, neonu, tlenu i krzemu w coraz krótszych skalach czasowych.

Łańcuch proton-proton (dominujący dla M ≤ 1,3 M_słońce):
  pp I:  4H → ⁴He + 2e’ + 2ν_e       (Q = 26,7 MeV)
  pp II: gałąź ⁷Be (istotna przy T > 14 mln K)
  pp III: gałąź ⁸B (produkuje wysokoenergetyczne neutrina wykryte w SNO)
  Generacja energii: ε_pp ∝ ρ X² T^4  (X = udział masowy wodoru)

Cykl CNO (dominujący dla M > 1,3 M_słońce):
  ¹²C + p → ¹³N + γ → ¹³C + e’ + ν
  ¹³C + p → ¹⁴N + γ
  ¹⁴N + p → ¹⁵O + γ  ← etap ograniczający tempo
  ¹⁵O → ¹⁵N + e’ + ν → + p → ¹²C + ⁴He (netto)
  Generacja energii: ε_CNO ∝ ρ X_CNO T^17

Proces potrójnego alfa (spalanie helu, T > 10^8 K):
  2 ⁴He → ⁸Be* (czas życia rezonansu 2,6×10^-16 s)
  ⁸Be + ⁴He → ¹²C* (stan Hoyle'a, 7,6644 MeV)  → ¹²C + 2γ
  ¹²C + ⁴He → ¹⁶O + γ  (konkurencyjne; określa stosunek C/O)
  Rezonans Hoyle'a przewidziany przez Hoyle'a (1954), potwierdzony przez Cooka i in.

Zaawansowane etapy spalania (masywne gwiazdy, T w jednostkach 10^8 K):
  Spalanie He:   T ~2,   czas życia ~10^6 lat   (>8 M_słońce)
  Spalanie C:    T ~5,   czas życia ~10^3 lat
  Spalanie Ne:   T ~15,  czas życia ~1 rok
  Spalanie O:    T ~20,  czas życia ~miesiące
  Spalanie Si:   T ~35,  czas życia ~dni
  Jądro żelazne: T ~70,  brak uwalniania energii → kolaps jądra
  Jądra grupy żelaza są punktem końcowym: energia wiązania osiąga szczyt przy ⁵⁶Fe/⁵⁸Ni

4. Supernowe kolapsu jądra

Gdy masywna gwiazda (>8 M_słońce) zgromadzi jądro żelazne masywniejsze niż granica Chandrasekhara (~1,4 M_słońce), ciśnienie degeneracji elektronowej nie jest już w stanie go podtrzymać. Jądro zapada się z mniej więcej rozmiaru Ziemi do 10-kilometrowej gwiazdy neutronowej w czasie krótszym niż sekunda, uwalniając ~3 × 10^53 ergów — energię wiązania grawitacyjnego gwiazdy neutronowej. To ≈ 100 razy więcej niż całkowita energia elektromagnetyczna wypromieniowana przez Słońce w ciągu całego jego 10-miliardowego życia, uwolniona w ciągu milisekund.

Masa Chandrasekhara (relatywistyczna granica degeneracji elektronowej):
  M_Ch = 5,83 Y_e² M_słońce   (Y_e = udział elektronów ≈ 0,5 dla ⁵⁶Fe)
  M_Ch ≈ 1,44 M_słońce    (rzeczywista wartość ~1,2-1,4 M_słońce z poprawkami)

Oś czasu kolapsu:
  t = 0 ms:    Jądro żelazne przekracza M_Ch; zaczyna się swobodny spadek
  t ~ 100 ms:  Jądro osiąga gęstość jądrową (ρ_jądr = 2,7×10^14 g/cm³)
  t ~ 110 ms:  Jądro odbija się; fala uderzeniowa startuje z prędkością ~100-200 km/s
  t ~ 200 ms:  Fala uderzeniowa zatrzymuje się — fotodysocjacja żelaza kosztuje 8,8 MeV/nukleon
  t ~ 500 ms:  Ogrzewanie neutrinowe ponownie energetyzuje falę (mechanizm BNNS)
  t ~ 1 s:     Udana eksplozja; fala przebija się ~15-20 godz. później
  t ~ dni:     Maksimum optyczne

Rozbłysk neutrinowy (potwierdzony dla SN 1987A):
  Całkowita E_ν: ~3 × 10^53 erg  (99% energii kolapsu)
  Czas trwania: ~10 s (chłodzenie Kelvina-Helmholtza gwiazdy neutronowej)
  Kamiokande II + IMB + Baksan wykryły 24 neutrina z SN 1987A (LMC, 160 tys. lat świetlnych)

Nukleosynteza w fazie eksplozywnej:
  proces r (szybki wychwyt neutronów): produkuje ~50% jąder A>100
  Ostatnie dowody: proces r w kilonowej ze zlania GS GW170817 (złoto, platyna, europ)
  proces p: fotorozpad jąder zarodkowych; izotopy bogate w protony
  proces νp: napędzany neutrinami wychwyt protonów we wczesnym wyrzucie

Szczytowa jasność supernowych kolapsu jądra:
  L_szczyt ~ 10^43 erg/s (optyczna)
  Zanik napędzany przez ⁵⁶Ni → ⁵⁶Co → ⁵⁶Fe (t_½ = 6,08 d, 77,2 d)
  Jasność ∝ M_Ni zsyntetyzowana (~0,07 M_słońce dla SN 1987A)

5. Gwiazdy neutronowe i pulsary

Pozostałe po supernowej kolapsu jądra jądro to gwiazda neutronowa — sfera o masie ~1,4 M_słońce skompresowana do promienia ~10 km, podtrzymywana ciśnieniem degeneracji neutronowej i odpychającym twardym rdzeniem oddziaływania jądrowego. Średnia gęstość przekracza gęstość jądrową; łyżeczka ważyłaby ≈ 5 × 10^8 ton. Jeśli jądro progenitora przekracza ~2–3 M_słońce, kolaps kontynuuje się poza stabilność gwiazdy neutronowej aż do czarnej dziury.

Równanie Tolmana-Oppenheimera-Volkoffa (TOV) równowagi hydrostatycznej w OTW:
  dP/dr = −G(M + 4πr³P/c²)(ρ + P/c²) / [r(r − 2GM/c²)]
  Redukuje się do newtonowskiego dP/dr = −ρGM/r² przy niskiej zwartości

Maksymalna masa gwiazdy neutronowej:
  Granica Tolmana-Oppenheimera-Volkoffa ~0,7 M_słońce (czysta materia neutronowa, miękkie EOS)
  Rzeczywista granica (sztywne jądrowe EOS): ~2,3–2,5 M_słońce
  Zmierzona: PSR J0952-0607 = 2,35 ± 0,17 M_słońce (2022, najcięższa potwierdzona GN)

Spowolnienie rotacji pulsara (promieniowanie dipola magnetycznego):
  P_kropka = −B²R^6Ω^3sin²α / (6Ic³)
  P = okres rotacji (najszybszy znany: J1748-2446ad, P = 1,396 ms)
  Wiek charakterystyczny: τ = P / (2 P_kropka)
  Powierzchniowe pole B:    B ~ 3,2×10^19 √(P P_kropka)  Gauss

Chłodzenie gwiazdy neutronowej:
  t < 100 lat:  Bezpośrednia Urca (jeśli M > M_DU), szybkie chłodzenie, T_s ~ 10^6 K
  t ~ 10^3 lat: Zmodyfikowana Urca dominująca, T_s ~ 5×10^5 K
  t ~ 10^5 lat: Chłodzenie fotonowe, T_s ~ 10^5 K
  Obserwowane przykłady: Kasjopeja A (350 lat) wykazała chłodzenie ~2% na dekadę (Ho i Heinke 2009)

Równanie stanu (EOS) gęstej materii:
  Hadronowe: SLy, APR4 (siły jądrowe + trójciałowe)
  Egzotyczne: hiperonowe, kondensat kaonowy, plazma kwarkowo-gluonowa
  Ograniczenie promienia (NICER): R = 12,35 ± 0,75 km dla J0030+0451 (Riley 2019)

6. Fale grawitacyjne ze spiralnego zbliżania się układów zwartych

Ogólna teoria względności przewiduje, że przyspieszające masy emitują fale grawitacyjne — zmarszczki w krzywiźnie czasoprzestrzeni rozchodzące się z prędkością światła. Układy podwójne obiektów zwartych (gwiazdy neutronowe lub czarne dziury) tracą energię na rzecz promieniowania grawitacyjnego, powodując kurczenie się orbity, aż oba ciała się zleją. LIGO/Virgo po raz pierwszy wykryło ten sygnał 14 września 2015 r. (GW150914): dwie czarne dziury o masach 36 i 29 M_słońce zlewające się w odległości ~410 Mpc, generujące sygnał wygaszania trwający ~0,2 sekundy.

Wzór kwadrupolowy na moc GW (wiodący rząd PN):
  P_GW = 32G^4 m_1²m_2²(m_1+m_2) / (5c^5 a^5)
  a = separacja orbitalna

Masa świergotu ℳ (określa ewolucję fazy):
  ℳ = (m_1 m_2)^{3/5} / (m_1+m_2)^{1/5}
  Bezpośrednio mierzalna z przemiatania częstotliwości f_GW(t)

Zanik orbitalny (Peters 1964):
  da/dt = −64G³m_1m_2(m_1+m_2) / (5c^5 a³)
  Czas do zlania: t_zlanie = 12a_0^4c^5 / (19G³ m_1m_2(m_1+m_2))
  PSR 1913+16 (pulsar Hulse'a-Taylora): t_zlanie ~ 300 mln lat; zanik orbitalny zgadza się z OTW z dokładnością 0,2%

Odkształcenie fali grawitacyjnej:
  h_+ = − (2Gℳ/c²r)(Gπfℳ/c^3)^{2/3} cos(2φ)
  h_× o tej samej amplitudzie, przesunięcie fazowe 90°

GW150914 (pierwsza detekcja, wrzesień 2015):
  Masy składowe: 35,6 M_słońce + 30,6 M_słońce
  Masa po zlaniu: 63,1 M_słońce  (wypromieniowane ΔE = 3,0 M_słońce c² ~ 5,4×10^54 erg)
  Szczytowe odkształcenie: h ~ 10^-21   (długość ramienia LIGO 4 km; ΔL ~ 4×10^-18 m, 1/1000 promienia protonu)
  Szczytowa jasność GW: ~3,6×10^56 erg/s (przyćmiewa wszystkie widoczne gwiazdy w obserwowalnym wszechświecie)

GW170817 (pierwsze zlanie GN-GN, sierpień 2017):
  Masy składowe: 1,17 + 1,36 M_słońce
  Powiązana kilonowa AT2017gfo: składowa niebieska (Sr) + czerwona (ciężki proces r)
  Rozbłysk gamma GRB 170817A przesunięty o 1,74 s → c_GW/c = 1 ± 10^-15 (testuje niezmienniczość Lorentza)
  Stała Hubble'a: H_0 = 70,0^{+12}_{-8} km/s/Mpc (multimessenger)

Od 2024 roku sieć LIGO-Virgo-KAGRA skatalogowała ponad 90 kandydatów na zlanie się układów zwartych w trakcie trzech serii obserwacyjnych (GWTC-3). Czwarta seria obserwacyjna (O4, 2023–2025) ma podwoić tę próbkę. Kolejna generacja detektorów — Einstein Telescope (podziemny, ramiona 10 km) i Cosmic Explorer (ramiona 40 km) — wykryje zlania się podwójnych czarnych dziur w większości obserwowalnego wszechświata, umożliwiając bezpośredni pomiar historii ekspansji kosmicznej bez obserwacji elektromagnetycznych.

Powiązane symulacje