Spotlight #25 — Kosmologia i wszechświat

Kosmologia to nauka o wszechświecie jako całości — jego pochodzeniu, ewolucji, wielkoskalowej strukturze i ostatecznym losie. W przeciwieństwie do większości gałęzi fizyki kosmolodzy nie mogą przeprowadzać kontrolowanych eksperymentów; obserwują jeden wszechświat w czasie kosmicznym i testują modele na podstawie jego statystycznych odcisków palców. Nasza kolekcja symulacji obejmuje teraz całą oś czasu — od pierwszych trzech minut nukleosyntezy po dzisiejszą sieć galaktyk.

1. Prawo Hubble'a — rozszerzający się wszechświat

W 1929 roku Edwin Hubble opublikował zależność, która na zawsze zmieniła kosmologię: galaktyki oddalają się od nas z prędkością proporcjonalną do odległości. Nie dlatego, że galaktyki poruszają się przez przestrzeń — lecz dlatego, że sama przestrzeń się rozciąga, co opisuje metryka FLRW.

Prawo Hubble'a:   v = H₀ · d
  v = prędkość recesji (km/s)
  d = odległość właściwa (Mpc)
  H₀ = stała Hubble'a (Planck 2018: 67,4 ± 0,5 km/s/Mpc)
       (bezpośredni pomiar SH0ES: 73,0 ± 1,0 km/s/Mpc)

Metryka FLRW:  ds² = −c²dt² + a(t)²[dr²/(1−kr²) + r²dΩ²]
  a(t) = czynnik skali (a(t₀) = 1 dzisiaj)
  k = krzywizna (0 dla płaskiego wszechświata)

Promień Hubble'a:  r_H = c/H₀ ≈ 14,4 mld lat świetlnych
  Galaktyki poza r_H oddalają się szybciej niż światło (nie jest to
  sprzeczność — żadna informacja nie jest przekazywana; rozszerza się
  przestrzeń, nie obiekty)

Równanie Friedmanna:  (ȧ/a)² = (8πG/3)ρ − kc²/a² + Λc²/3
  ρ = całkowita gęstość energii (materia + promieniowanie + ciemna energia)
  Λ = stała kosmologiczna (ciemna energia ≈ 68% ρ_całkowitej dzisiaj)

Symulacja umieszcza 80 galaktyk w rozszerzającej się powłoce i koloruje je według przesunięcia ku czerwieni — niebieskie galaktyki bliskie, czerwone odległe. Ustaw H₀ na wartość Plancka (67,4), SH0ES (73,0) lub oryginalną wartość Hubble'a z 1929 roku (500!), by zobaczyć, dlaczego „napięcie Hubble'a" wciąż zajmuje kosmologów. Sfera Hubble'a — granica, na której prędkość recesji zrównuje się z prędkością światła — pulsuje w czasie rzeczywistym.

2. Pierwotna nukleosynteza — pierwsze trzy minuty

W pierwszej sekundzie po Wielkim Wybuchu wszechświat był plazmą kwarkowo-gluonową zbyt gorącą, by mogły powstać jądra atomowe. Do trzeciej minuty ochłodził się na tyle, że protony i neutrony zaczęły łączyć się w najlżejsze jądra. Słynna książka Stevena Weinberga nazywa tę erę „Pierwsze trzy minuty" — a to właśnie ona ustaliła skład pierwiastkowy całego obserwowalnego wszechświata.

Temperatura w funkcji czasu (era zdominowana przez promieniowanie):
  T(t) ≈ 10¹⁰ K / √t   (t w sekundach)

Wymrożenie stosunku neutron-proton (T ≈ 0,8 MeV, t ≈ 1 s):
  n/p = exp(−ΔQ/kT) ≈ 1/7   (ΔQ = mₙ−mₚ = 1,293 MeV)

Kluczowe reakcje (Wagoner-Fowler-Hoyle, 1967):
  p + n  → ²H + γ           (wąskie gardło deuteru pokonane przy T ≈ 0,07 MeV)
  ²H + p → ³He + γ
  ³He + n → ⁴He + p
  ⁷Be + e⁻ → ⁷Li + νₑ

Przewidywane pierwotne obfitości (frakcje masowe):
  ⁴He:  Yₚ ≈ 0,247  (~24,7% masy)
  ²H:   X_D ≈ 2,5 × 10⁻⁵  (czuły na η = n_b/n_γ)
  ³He:  X ≈ 1,0 × 10⁻⁵
  ⁷Li:  X ≈ 4,7 × 10⁻¹⁰  (problem ⁷Li: teoria ~3× więcej niż obserwacje)

Stosunek barionów do fotonów η ≈ 6,1 × 10⁻¹⁰ (ustalony przez obserwacje CMB)

Symulacja uruchamia uproszczoną sieć ośmiu kluczowych reakcji, śledząc obfitości w miarę stygnięcia wszechświata od 10 MeV do 10 keV. Przesuń suwak stosunku barionów do fotonów η i obserwuj, jak końcowa wydajność helu-4 jest niemal niewrażliwa na zmiany — dowód solidności BBN — podczas gdy deuter i lit-7 zmieniają się drastycznie, dostarczając najściślejszych ograniczeń na gęstość barionową wszechświata.

3. Mikrofalowe promieniowanie tła — piki akustyczne

380 000 lat po Wielkim Wybuchu wszechświat ochłodził się na tyle, by elektrony mogły połączyć się z protonami w neutralny wodór. W tej epoce „rekombinacji" fotony oddzieliły się od materii i popłynęły swobodnie, tworząc CMB — niemal idealne ciało doskonale czarne o T = 2,725 K, wypełniające całe niebo. Ale nie jest ono idealnie jednorodne: fluktuacje temperatury rzędu jednej stutysięcznej kodują oscylacje akustyczne pierwotnej plazmy barionowo-fotonowej.

Kątowe widmo mocy:  Cₗ = (2/π) ∫ Pᵩ(k) |Δₗ(k)|² k² dk
  l = moment multipolowy ≈ 180°/θ  (l~200 odpowiada skali ~1°)
  Pᵩ(k) ∝ kⁿˢ⁻¹, nₛ ≈ 0,965  (niemal skalo-niezmiennicze, zgodnie z przewidywaniami inflacji)

Oscylacje akustyczne — horyzont dźwiękowy przy rekombinacji:
  rs = ∫₀ᵗ_rec cs dt / a ≈ 147 Mpc  (standardowy przymiar!)
  cs = c/√3(1+R),  R = 3ρ_b/(4ρ_γ)  (parametr obciążenia barionowego)

Pozycje pików (płaski ΛCDM):
  l₁ ≈ 220  (pierwszy pik kompresji — idealnie zgodny z WMAP, Planck)
  l₂ ≈ 540  (rozrzedzenie — tłumione przez obciążenie barionowe)
  l₃ ≈ 810  (drugi pik kompresji — czuły na gęstość ciemnej materii)
  Stosunek l₂/l₁ ↓ wraz ze wzrostem Ω_b (więcej barionów tłumi parzyste piki)

Kluczowe parametry kosmologiczne zakodowane w pikach:
  Ω_b h² (gęstość barionów)      → stosunek wysokości pików nieparzystych/parzystych
  Ω_c h² (zimna ciemna materia)  → obwiednia zaniku wysokości pików
  Ω_Λ    (ciemna energia)        → późny efekt ISW, plateau przy niskim l
  θₛ     (skala kątowa)          → pozycje pików (geometria)

Symulacja CMB wyświetla kątowe widmo mocy D_ℓ = ℓ(ℓ+1)C_ℓ/2π obok punktów danych Planck 2018. Dostosuj gęstość barionów Ω_b, zimną ciemną materię Ω_c i stałą Hubble'a H₀, by dopasować piki, i zobacz, dlaczego CMB określa się „kamieniem z Rosetty kosmologii" — pojedynczy obraz kodujący wszystkie sześć fundamentalnych parametrów ΛCDM.

4. Ciemna materia — krzywe rotacji galaktyk

Pomiary krzywych rotacji galaktyk spiralnych wykonane przez Verę Rubin w latach 70. XX wieku ujawniły głęboką tajemnicę: gwiazdy na skraju galaktyki krążą tak samo szybko, jak te bliżej centrum, zamiast zwalniać, jak wymaga tego prawo Keplera. Jedynym spójnym wyjaśnieniem jest ogromne, niewidzialne halo materii, które ani nie emituje, ani nie pochłania światła — ciemna materia. Stanowi ona ≈27% gęstości energii wszechświata.

Oczekiwanie keplerowskie (tylko masa widzialna):
  v_circ(r) = √(GM(r)/r)  → v ∝ r⁻¹/² w zewnętrznym dysku

Obserwowane (płaskie krzywe rotacji):
  v_circ ≈ const ≈ 200 km/s  aż do r > 30 kpc
  ⟹  M(r) ∝ r  (gęstość ρ ∝ r⁻²)

Halo ciemnej materii NFW (Navarro-Frenk-White, 1997):
  ρ_NFW(r) = ρ₀ / [(r/rs)(1 + r/rs)²]
  rs = promień skali (zazwyczaj 15–25 kpc)
  Masa zamknięta M_NFW(r) = 4πρ₀rs³ [ln(1+r/rs) − r/(r+rs)]

Barionowa relacja Tully'ego-Fishera:
  v_flat⁴ = G · a₀ · M_bary   (a₀ ≈ 1,2 × 10⁻¹⁰ m/s² — MOND)
  ΛCDM i MOND odtwarzają relację Tully'ego-Fishera odmiennie;
  rozróżnienie między nimi wymaga przeglądów słabego soczewkowania grawitacyjnego.

Dowody na istnienie ciemnej materii (niezależne):
  • Krzywe rotacji galaktyk (Rubin i in., 1970)
  • Soczewkowanie grawitacyjne (Clowe i in., 2006 — Gromada Kuli)
  • Wysokości pików akustycznych CMB (Spergel i in., 2003)
  • Symulacje N-ciałowe wielkoskalowej struktury (Springel i in., 2005)

Symulacja ciemnej materii pozwala włączyć i wyłączyć halo NFW oraz obserwować, jak krzywa rotacji przekształca się z keplerowskiego spadku w płaską. Dostosuj parametr koncentracji halo c = r₂₀₀/rs i zobacz, jak bardziej skoncentrowane halo tworzy strome profile wewnętrzne — „problem cusp-core", który kwestionuje przewidywania czystej zimnej ciemnej materii.

5. Soczewkowanie grawitacyjne — pierścienie i łuki Einsteina

Ogólna teoria względności przewiduje, że masa zakrzywia czasoprzestrzeń, uginając światło pochodzące od źródeł w tle. Dla masy punktowej działającej jako idealna soczewka, źródło znajdujące się dokładnie za nią wygląda jak idealny pierścień — pierścień Einsteina, po raz pierwszy sfotografowany w 1988 roku. Soczewkowanie jest dziś podstawowym narzędziem do mapowania ciemnej materii w gromadach galaktyk i pomiaru stałej Hubble'a za pomocą opóźnień czasowych.

Kąt ugięcia (masa Schwarzschilda M):
  α̂ = 4GM / (c²ξ)   (ξ = odległość najbliższego zbliżenia)
  Współczynnik 2× większy niż przewidywanie newtonowskie (potwierdzone w 1919 r., Eddington)

Równanie soczewki (cienka soczewka, współrzędne kątowe):
  β = θ − θ_E²/θ    (β = rzeczywisty kąt źródła, θ = kąt obrazu)
  θ_E = √(4GM D_ls / (c² D_l D_s))  — promień Einsteina

  D_l, D_s, D_ls = kątowe odległości średnicowe
  (soczewka, źródło, soczewka-źródło)

Pierścień Einsteina:  β = 0  →  θ = θ_E  (obraz pierścieniowy)
  Dla soczewkowania galaktyka-galaktyka: θ_E ≈ 1–2 sekund kątowych

Wielokrotne obrazy:  β < θ_E  →  dwa lub więcej obrazów
  Powiększenie:  μ = |θ/β · dθ/dβ| (może wynosić ×100+ dla źródeł punktowych)

Mikrosoczewkowanie (MACHO/WIMP):
  Charakterystyczny czas t_E = θ_E / μ_rel ≈ dni-lata
  Kształt krzywej blasku: Paczyński A = (u²+2)/[u√(u²+4)], u = β/θ_E

Umieść masę soczewki w dowolnym miejscu na interaktywnym płótnie, dostosuj jej masę oraz położenie źródła i obserwuj, jak topologia obrazu zmienia się między pojedynczym obrazem, dwoma obrazami a pierścieniem Einsteina. Siatka gwiazd tła ujawnia wzór powiększenia — kluczowe narzędzie w poszukiwaniach MACHO/ciemnej materii.

6. Gwiazdy podwójne — Kepler spotyka spektroskopię

Ponad połowa wszystkich gwiazd typu słonecznego istnieje w układach podwójnych. Układy podwójne są fundamentalne dla astrofizyki: dostarczają jedynego bezpośredniego pomiaru mas gwiazd, kalibrują skale odległości za pomocą krzywych blasku zaćmień, są prekursorami supernowych typu Ia i ostatecznie produkują źródła fal grawitacyjnych wykrywalne przez LIGO.

Problem dwóch ciał → sprowadzony do jednego ciała z μ = m₁m₂/(m₁+m₂):
  r̈ = −G(m₁+m₂) r̂ / r²

Półoś wielka i okres (trzecie prawo Keplera):
  a³/P² = G(m₁+m₂)/(4π²)   w jednostkach SI
       → a [AU]³/P [lat]² = (m₁+m₂)/M☉   w jednostkach słonecznych

Pozycje względem środka masy:
  r₁ = m₂/(m₁+m₂) · r,   r₂ = −m₁/(m₁+m₂) · r

Prędkość radialna (podwójna spektroskopowa, nachylenie i):
  K₁ = (2πa₁ sin i/P)/√(1−e²)   (półamplituda, gwiazda 1)
  Funkcja masy: f(m) = m₂³ sin³i / (m₁+m₂)² = K₁³P/(2πG)
  M sin i wyłącznie z K₁ i K₂ — bez znajomości odległości!

Czas trwania spiralowania grawitacyjno-falowego (Peters, 1964):
  T_zlania = (12c⁵/85G³) · a₀⁴(m₁+m₂)/(m₁m₂) · (1−e²)^(7/2)
  Dla dwóch gwiazd neutronowych po 1,4M☉, a₀=2R☉: T ≈ 300 mln lat

Dostosuj stosunek mas q = m₂/m₁, mimośród e i okres orbitalny oraz obserwuj, jak obie gwiazdy zakreślają swoje elipsy wokół wspólnego barycentrum. Krzywa prędkości radialnej — wielkość obserwacyjna wykorzystywana przez Verę Rubin i niezliczonych innych astronomów — aktualizuje się w czasie rzeczywistym. Ustaw wysoki mimośród, by zobaczyć, jak orbita kołowieje pod wpływem emisji fal grawitacyjnych w skalach czasu astronomicznego.

Standardowy model kosmologiczny

Wszystkie sześć symulacji łączy się z modelem ΛCDM — obecnie najlepszym opisem wszechświata. Jego sześć swobodnych parametrów (H₀, Ω_b, Ω_c, Ω_Λ, nₛ, A_s) jest wyznaczanych z dokładnością do procenta wyłącznie na podstawie CMB i niezależnie potwierdzanych przez BBN, przeglądy soczewkowania grawitacyjnego oraz grupowanie galaktyk. Pozostają jednak napięcia: napięcie Hubble'a (H₀), napięcie σ₈ (amplituda wielkoskalowej struktury) i problem litu w BBN sugerują, że model może być niepełny.

Symulacje szczególnej i ogólnej teorii względności, które stanowią teoretyczny kręgosłup kosmologii, znajdziesz w Nauka #21 — Szczególna teoria względności oraz Spotlight #16 — Szczególna teoria względności. Kwantowa strona wczesnego wszechświata — fizyka cząstek powyżej skali elektrosłabej — jest omówiona w naszym Spotlight fizyki kwantowej.