Na cenzurowanym #38 – Kosmologia i ciemna materia: CMB, inflacja, symulacje N-ciałowe i wielkoskalowa struktura

Współczesna kosmologia opiera się na trzech filarach obserwacyjnych: kosmicznym mikrofalowym promieniowaniu tła (CMB), wielkoskalowym rozkładzie galaktyk oraz obfitości lekkich pierwiastków pochodzących z pierwotnej nukleosyntezy Wielkiego Wybuchu. Razem potwierdzają one model ΛCDM (Lambda — zimna ciemna materia) — wszechświat liczący 13,8 mld lat, płaski z dokładnością do 0,2%, złożony w 5% z barionów, 27% z zimnej ciemnej materii i 68% z ciemnej energii. Ten wpis analizuje fizykę stojącą za każdym z tych ograniczeń.

1. Kosmologia inflacyjna

Standardowy model Gorącego Wielkiego Wybuchu napotyka trzy problemy dostrajania: problem horyzontu (dlaczego temperatura CMB jest jednorodna z dokładnością do 1 części na 10⁵ w obszarach przyczynowo rozłączonych w momencie rekombinacji?), problem płaskości (dlaczego parametr krzywizny przestrzennej |Ω_k| < 0,002?) oraz problem monopoli magnetycznych (Wielkie Teorie Unifikacji przewidują reliktową gęstość monopoli znacznie przewyższającą obserwacje). Alan Guth (1981) zaproponował inflację kosmiczną — epokę niemal wykładniczej ekspansji — która rozwiązuje wszystkie trzy problemy jednocześnie.

Równanie Friedmanna:
  H² = (ȧ/a)² = (8πG/3)ρ − k/a² + Λ/3

Ekspansja inflacyjna (przybliżenie powolnego staczania):
  a(t) ≈ a_i exp(H_inf t)   gdzie H_inf ≈ stała
  Czas podwojenia: t_2 = ln2 / H_inf ≈ 10−³&sup7; s  (inflacja w skali GUT)

Parametry powolnego staczania:
  ε = −Ḣ/H² = (M_Pl²/2)(V'/V)²   (płaskość potencjału)
  η = M_Pl² V''/V   (krzywizna potencjału)
  Inflacja: ε < 1,  kończy się gdy ε = 1

Liczba e-foldów:
  N = ∫ H dt ≈ 50–60  (wymagana do rozwiązania problemów horyzontu i płaskości)
  Obszar przyczynowy rośnie od skali subplanckowskiej do > promienia Hubble'a

Fluktuacje kwantowe → pierwotne widmo mocy:
  P_s(k) = A_s (k/k_0)^(n_s − 1)
  n_s = 1 − 6ε + 2η ≈ 0,965  (Planck 2018; n_s=1 = Harrison-Zeldovich)
  A_s ≈ 2,1 × 10−&sup9;  przy skali odniesienia k_0 = 0,05 Mpc−¹

Mody tensorowe (pierwotne fale grawitacyjne):
  Stosunek tensorowo-skalarny r = 16ε
  Planck+BICEP/Keck 2021: r < 0,036  (przy poziomie ufności 95%)
          

2. Kosmiczne mikrofalowe promieniowanie tła

CMB to promieniowanie termiczne wyemitowane, gdy wszechświat ochłodził się poniżej ~3000 K przy przesunięciu ku czerwieni z ≈ 1100, a wodór uległ rekombinacji, czyniąc wszechświat przezroczystym. Jego niemal doskonałe widmo ciała doskonale czarnego o temperaturze 2,725 K oraz niewielkie anizotropie temperatury ΔT/T ∼ 10⁻⁵ zapisują warunki początkowe wszechświata i fizykę plazmy fotonowo-barionowej.

Rekombinacja:
  z_rec ≈ 1100,  T_rec ≈ 3000 K,  t_rec ≈ 380 000 lat
  Temperatura CMB dzisiaj: T_0 = 2,7255 K

Kątowe widmo mocy C_ℓ:
  ΔT/T(θ) = Σ_(ℓm) a_(ℓm) Y_(ℓm)(θ,φ)
  C_ℓ = ⟨|a_(ℓm)|²⟩   (niezmiennik obrotowy)

Piki akustyczne (efekt Sachsa-Wolfe'a + BAO):
  1. pik: ℓ ≈ 220  (horyzont dźwiękowy przy rekombinacji, θ ≈ 1°)
    → wszechświat jest płaski: Ω_k ≈ 0
  2. pik: ℓ ≈ 540  (tłumiony obciążeniem barionowym)
    → gęstość barionów Ω_b h² ≈ 0,0224
  3. pik: ℓ ≈ 810  (ciemna materia podtrzymuje studnie potencjału)
    → dowód na Ω_DM h²

Horyzont dźwiękowy przy ostatnim rozproszeniu:
  r_s = ∫_0^(t_rec) c_s(t)/ a(t) dt ≈ 147 Mpc  (współporuszający się)
  c_s = c / √(3(1 + R))   gdzie R = 3ρ_b/(4ρ_γ)

Zintegrowany efekt Sachsa-Wolfe'a (ISW):
  Późny ISW: fotony zyskują energię w rosnących studniach potencjału
  (ΛCDM: powoduje korelację między CMB a gęstością galaktyk przy ℓ < 20)

Parametry Planck 2018 (najlepsze dopasowanie ΛCDM):
  H_0 = 67,4 ± 0,5 km/s/Mpc  (napięcie Hubble'a: wartość lokalna = 73,2 ± 1,3)
  Ω_b h² = 0,02237,  Ω_DM h² = 0,1200,  Ω_Λ = 0,685
  n_s = 0,9649,  A_s = 2,1 × 10−&sup9;,  τ = 0,054
          

3. Ciemna materia: kandydaci i detekcja

Ciemna materia stanowi ~85% całej materii we wszechświecie. Jej istnienie wnioskuje się z krzywych rotacji galaktyk, soczewkowania grawitacyjnego (najbardziej dramatyczny przykład to Gromada Kuli, gdzie gaz barionowy i ciemna materia rozdzieliły się podczas zderzenia gromad), anizotropii CMB oraz wielkoskalowej struktury. Mimo 40 lat wysiłków eksperymentalnych żaden sygnał niegrawitacyjny nie został potwierdzony.

Krzywe rotacji galaktyk (Rubin i Ford, 1970):
  Oczekiwane: v_c(r) = √(GM(r)/r) → maleje jak r−¹² poza widocznym dyskiem
  Zaobserwowane: v_c(r) ≈ const ≈ 200 km/s aż do >50 kpc
  Sugeruje M(r) ∝ r → halo ciemnej materii z ρ ∝ 1/r (izotermiczne) lub NFW

Profil halo NFW (Navarro-Frenk-White):
  ρ(r) = ρ_s / [(r/r_s)(1 + r/r_s)²]
  r_s = promień skali,  ρ_s = charakterystyczna gęstość
  Parametr koncentracji c = r_200/r_s  (c ≈ 10 dla halo o masie Drogi Mlecznej)

WIMP (słabo oddziałująca masywna cząstka):
  m ≈ 10 GeV–10 TeV,  <σv⟩ ≈ 3 × 10−²&sup6; cm³/s (relikt termiczny)
  LHC: brak sygnałów SUSY
  Detekcja bezpośrednia: LZ (LUX-ZEPLIN, 7 t ciekłego Xe): brak sygnału na 2024 rok
  Detekcja pośrednia: poszukiwania promieniowania gamma Fermi-LAT → górne granice na <σv⟩

Aksjon:
  m_a ≈ 10−&sup6;–10−³ eV (okno aksjonu QCD)
  Mechanizm Peccei-Quinn rozwiązuje silny problem CP
  Eksperyment ADMX: rezonansowa konwersja w mikrofalowej wnęce B·a → foton
  ABRACADABRA: detekcja aksjonu na elementach skupionych

Pierwotne czarne dziury (PCD):
  Powstały w erze promieniowania z fluktuacji gęstości
  Okno masowe: 10−¹³–45 M_słon. (ograniczenia mikrosoczewkowania grawitacyjnego + fal grawitacyjnych)
  PCD o masie asteroidalnej (10−¹&sup7;–10−¹¹ M_słon.): wciąż realny kandydat na ciemną materię
          

4. Formowanie struktur N-ciałowych

Pole gęstości ciemnej materii ewoluuje pod wpływem grawitacji od drobnych fluktuacji początkowych (Δρ/ρ ∼ 10⁻⁵ przy rekombinacji) do bogatej pajęczyny halo, nici, arkuszy i pustek obserwowanych dzisiaj. Współczesne kosmologiczne symulacje N-ciałowe śledzą 10¹⁰–10¹¹ cząstek, wykorzystując adaptacyjną grawitację typu tree-code oraz periodyczne warunki brzegowe.

Liniowa teoria zaburzeń (era zdominowana przez materię):
  &ddot;δ + 2Hδ̇ − 4πGρ_m δ = 0
  Mod rosnący: δ(t) ∝ D(t) (liniowy czynnik wzrostu)
  D(z) ≈ H(z) ∫_z^∞ (1+z')/H(z')³ dz'

Długość Jeansa (próg kolapsu):
  λ_J = c_s √(π/(Gρ))
  Ciemna materia: c_s ≈ 0 → kolapsują wszystkie skale (zimna CM)
  Długość swobodnego przelotu λ_FS: neutrina i ciepła CM tłumią moc na małych skalach

Funkcja masy halo Pressa-Schechtera (1974):
  dn/dM = √(2/π) (ρ_m/M²) |d ln σ/d ln M| (δ_c/σ) exp(−δ_c²/(2σ²))
  δ_c ≈ 1,686 (próg kolapsu sferycznego)
  σ(M) = fluktuacja gęstości RMS w sferze o masie M

Grawitacja TreePM:
  Krótkozasięgowa: oktadrzewo Barnesa-Huta, sumowanie sił O(N log N)
  Długozasięgowa: particle-mesh (PM) z FFT na siatce 3D
  Długość wygładzania ε ≈ 1–5% średniej odległości między cząstkami (unika rozpraszania dwuciałowego)

Symulacja Millennium (Springel, 2005):
  2160³ ≈ 10^10 cząstek ciemnej materii, sześcian 500 Mpc/h
  Illustris-TNG (2018): sześcian 300³ Mpc, pełna fizyka barionowa (gaz, gwiazdy, sprzężenie zwrotne CD)
          

5. Baryonowe oscylacje akustyczne

Przed rekombinacją bariony i fotony były ściśle sprzężone poprzez rozpraszanie Thomsona, tworząc plazmę fotonowo-barionową, która podtrzymywała fale dźwiękowe. W momencie rekombinacji fotony się oddzieliły, a fale dźwiękowe zamarły, pozostawiając zamrożony odcisk — skalę baryonowych oscylacji akustycznych r_s ≈ 147 Mpc — która działa jako „standardowa miarka” w rozkładzie galaktyk.

Skala BAO przy niskim przesunięciu ku czerwieni:
  Współporuszający się horyzont dźwiękowy r_s ≈ 147 Mpc  (taki sam jak w CMB)
  Ujawnia się jako pierścień w dwupunktowej funkcji korelacji ξ(r)
  Pik w ξ(r) w pobliżu r ≈ 150 Mpc (odległość współporuszająca się)

Kątowa odległość średnicowa D_A(z):
  D_A = (1/(1+z)) ∫_0^z c/H(z') dz'
  BAO mierzy D_A(z)/r_s oraz H(z)·r_s wzdłuż/poprzecznie do linii widzenia

Parametr Hubble'a:
  H(z)² = H_0² [Ω_m(1+z)³ + Ω_r(1+z)&sup4; + Ω_k(1+z)² + Ω_Λ f_DE(z)]
  Równanie stanu ciemnej energii: w = p/ρc² (ΛCDM: w = −1)

Wyniki BAO z pierwszego roku DESI (2024):
  Zgodne z płaskim ΛCDM z dokładnością <2%, 11,5 miliona galaktyk
  Słaba przesłanka: w_0 = −0,55, w_a = −1,32 (dynamiczna ciemna energia, 2,5σ)
  Euclid (wystrzelony w 2023): zmapuje 1,5 × 10^9 galaktyk do z < 2

Krzyżowa weryfikacja frakcji barionowej:
  f_b = Ω_b / Ω_m ≈ 0,16  (z CMB + BBN)
  Obserwowane w gromadach rentgenowskich: f_gas ≈ 0,12 (+ gwiazdy ≈ 0,04) → zgodne
          

6. Kosmiczna pajęczyna

Ciemna materia zapada się preferencyjnie wzdłuż arkuszy i nici, zanim skurczy się w zwarte halo na ich przecięciach — przybliżenie Zeldowicza przewiduje takie zachowanie typu „naleśnikowego” na podstawie wartości własnych tensora deformacji. Galaktyki formują się wewnątrz halo ciemnej materii; wynikowy wzorzec skupiania się galaktyk, obserwowany w przeglądach obejmujących miliardy lat świetlnych, ujawnia kosmiczną pajęczynę.

Przybliżenie Zeldowicza (1970):
  x(q, t) = q − D(t) ∇ψ(q)
  ψ(q) = potencjał grawitacyjny; D(t) = liniowy czynnik wzrostu
  Kolaps wzdłuż trzech osi: najpierw wzdłuż najkrótszego wektora własnego (naleśniki)

Klasyfikacja T-web (Hahn, 2007):
  Tensor pływowy T_ij = ∂²φ/∂x_i∂x_j
  # dodatnich wartości własnych: 0 → pustki, 1 → arkusze, 2 → nici, 3 → halo

Statystyka kosmicznej pajęczyny:
  Halo (węzły):      ~5% objętości, ~80% masy halo
  Nici:              ~5% objętości, ~60% rozproszonego gazu
  Arkusze:           ~15% objętości
  Pustki:            ~75% objętości, <1% materii

Baryonowy Oscylacyjny Przegląd Spektroskopowy SDSS (BOSS):
  1,3 miliona galaktyk, 0,2 < z < 0,75
  Widmo mocy P(k) zmierzone z dokładnością 1% → ograniczenia na ciemną energię

Szeroki przegląd Euclid (2025–2030):
  1/3 nieba, 10^9 źródeł, spektroskopowe z dla 35 milionów
  Cel: w ±0,01, σ_8 ±0,005
  Pierwsze wydanie danych Euclid (DR1): czerwiec 2026 (15 milionów galaktyk)
          

Wypróbuj te symulacje