Спотлайт #38 – Космологія та темна матерія: РМФ, інфляція, N-тільні симуляції та великомасштабна структура

Сучасна космологія спирається на три спостережні опори: реліктовий мікрохвильовий фон (РМФ), великомасштабний розподіл галактик та поширеність легких елементів, утворених первинним нуклеосинтезом. Разом вони підтверджують модель ΛCDM (Лямбда — холодна темна матерія) — всесвіт віком 13,8 млрд років, плаский з точністю до 0,2%, що складається з 5% баріонів, 27% холодної темної матерії та 68% темної енергії. Цей допис розглядає фізику, що стоїть за кожним з обмежень.

1. Інфляційна космологія

Стандартна модель Гарячого Великого вибуху стикається з трьома проблемами тонкого налаштування: проблемою горизонту (чому температура РМФ однорідна з точністю до 1 частини на 10⁵ в областях, причинно роз'єднаних на момент рекомбінації?), проблемою площинності (чому параметр просторової кривизни |Ω_k| < 0,002?) та проблемою магнітних монополів (Великі об'єднані теорії передбачають релікт-щільність монополів, що набагато перевищує спостережувану). Алан Гут (1981) запропонував космічну інфляцію — епоху майже експоненційного розширення — яка вирішує всі три проблеми одночасно.

Рівняння Фрідмана:
  H² = (ȧ/a)² = (8πG/3)ρ − k/a² + Λ/3

Інфляційне розширення (наближення повільного скочування):
  a(t) ≈ a_i exp(H_inf t)   де H_inf ≈ константа
  Час подвоєння: t_2 = ln2 / H_inf ≈ 10−³&sup7; с  (інфляція масштабу ВОТ)

Параметри повільного скочування:
  ε = −Ḣ/H² = (M_Pl²/2)(V'/V)²   (пласкість потенціалу)
  η = M_Pl² V''/V   (кривизна потенціалу)
  Інфляція: ε < 1,  закінчується коли ε = 1

Кількість е-згортань:
  N = ∫ H dt ≈ 50–60  (потрібно для вирішення проблем горизонту й площинності)
  Причинна ділянка зростає від субпланківського масштабу до > радіуса Хаббла

Квантові флуктуації → первинний спектр потужності:
  P_s(k) = A_s (k/k_0)^(n_s − 1)
  n_s = 1 − 6ε + 2η ≈ 0,965  (Planck 2018; n_s=1 = Гаррісон-Зельдович)
  A_s ≈ 2,1 × 10−&sup9;  на опорному масштабі k_0 = 0,05 Мпк−¹

Тензорні моди (первинні гравітаційні хвилі):
  Відношення тензор-до-скаляр r = 16ε
  Planck+BICEP/Keck 2021: r < 0,036  (при 95% довірчому рівні)
          

2. Реліктовий мікрохвильовий фон

РМФ — це теплове випромінювання, випущене, коли всесвіт охолов нижче ~3000 К при червоному зміщенні z ≈ 1100 і водень рекомбінував, зробивши всесвіт прозорим. Його майже ідеальний спектр абсолютно чорного тіла при 2,725 К та малі температурні анізотропії ΔT/T ∼ 10⁻⁵ кодують початкові умови всесвіту та фізику фотон-баріонної плазми.

Рекомбінація:
  z_rec ≈ 1100,  T_rec ≈ 3000 К,  t_rec ≈ 380 000 років
  Температура РМФ сьогодні: T_0 = 2,7255 К

Кутовий спектр потужності C_ℓ:
  ΔT/T(θ) = Σ_(ℓm) a_(ℓm) Y_(ℓm)(θ,φ)
  C_ℓ = ⟨|a_(ℓm)|²⟩   (обертово-інваріантний)

Акустичні піки (ефект Сакса-Вольфа + БАО):
  1-й пік: ℓ ≈ 220  (звуковий горизонт при рекомбінації, θ ≈ 1°)
    → всесвіт плаский: Ω_k ≈ 0
  2-й пік: ℓ ≈ 540  (пригнічений баріонним навантаженням)
    → щільність баріонів Ω_b h² ≈ 0,0224
  3-й пік: ℓ ≈ 810  (темна матерія підтримує потенціальні ями)
    → свідчення на користь Ω_DM h²

Звуковий горизонт при останньому розсіянні:
  r_s = ∫_0^(t_rec) c_s(t)/ a(t) dt ≈ 147 Мпк  (супутня відстань)
  c_s = c / √(3(1 + R))   де R = 3ρ_b/(4ρ_γ)

Інтегрований ефект Сакса-Вольфа (ІСВ):
  Пізній ІСВ: фотони набувають енергії у зростаючих потенціальних ямах
  (ΛCDM: спричиняє кореляцію між РМФ та густиною галактик при ℓ < 20)

Параметри Planck 2018 (найкраще узгодження ΛCDM):
  H_0 = 67,4 ± 0,5 км/с/Мпк  (напруга Хаббла: локальне значення = 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. Темна матерія: кандидати та виявлення

Темна матерія становить ~85% усієї матерії у всесвіті. Її існування виводиться з кривих обертання галактик, гравітаційного лінзування (найдраматичніше — Скупчення Кулі, де баріонний газ і темна матерія розділилися під час злиття скупчень), анізотропій РМФ та великомасштабної структури. Попри 40 років експериментальних зусиль, жодного негравітаційного сигналу підтверджено не було.

Криві обертання галактик (Рубін і Форд, 1970):
  Очікувано: v_c(r) = √(GM(r)/r) → спадає як r−¹² за межами видимого диска
  Спостережено: v_c(r) ≈ const ≈ 200 км/с аж до >50 кпк
  Свідчить про M(r) ∝ r → гало темної матерії з ρ ∝ 1/r (ізотермальне) або NFW

Профіль гало NFW (Наварро-Френк-Уайт):
  ρ(r) = ρ_s / [(r/r_s)(1 + r/r_s)²]
  r_s = масштабний радіус,  ρ_s = характерна густина
  Параметр концентрації c = r_200/r_s  (c ≈ 10 для гало маси Чумацького Шляху)

WIMP (слабко взаємодіюча масивна частинка):
  m ≈ 10 ГеВ–10 ТеВ,  <σv⟩ ≈ 3 × 10−²&sup6; см³/с (термальний релікт)
  ВАК: поки без сигналів SUSY
  Пряме виявлення: LZ (LUX-ZEPLIN, 7 т рідкого Xe): без сигналу станом на 2024
  Непряме: пошуки гамма-променів Fermi-LAT → верхні межі на <σv⟩

Аксіон:
  m_a ≈ 10−&sup6;–10−³ еВ (вікно КХД-аксіона)
  Механізм Печеї-Квінн вирішує сильну CP-проблему
  Експеримент ADMX: резонансна конверсія у мікрохвильовому резонаторі B·a → фотон
  ABRACADABRA: виявлення аксіона на зосереджених елементах

Первинні чорні діри (ПЧД):
  Сформувалися в епоху випромінювання з флуктуацій густини
  Віконце мас: 10−¹³–45 M_сонця (обмеження гравітаційного мікролінзування + ГХ)
  ПЧД астероїдної маси (10−¹&sup7;–10−¹¹ M_сонця): все ще життєздатний кандидат на темну матерію
          

4. N-тільне формування структур

Поле густини темної матерії еволюціонує під дією гравітації від крихітних початкових флуктуацій (Δρ/ρ ∼ 10⁻⁵ при рекомбінації) до багатої павутини гало, ниток, шарів і порожнин, що спостерігаються сьогодні. Сучасні космологічні N-тільні симуляції відстежують 10¹⁰–10¹¹ частинок, використовуючи адаптивну деревоподібну гравітацію та періодичні граничні умови.

Лінійна теорія збурень (епоха домінування речовини):
  &ddot;δ + 2Hδ̇ − 4πGρ_m δ = 0
  Зростаюча мода: δ(t) ∝ D(t) (лінійний фактор зростання)
  D(z) ≈ H(z) ∫_z^∞ (1+z')/H(z')³ dz'

Довжина Джинса (поріг колапсу):
  λ_J = c_s √(π/(Gρ))
  Темна матерія: c_s ≈ 0 → колапсують усі масштаби (холодна ТМ)
  Довжина вільного пробігу λ_FS: нейтрино і тепла ТМ пригнічують потужність на малих масштабах

Функція маси гало Пресса-Шехтера (1974):
  dn/dM = √(2/π) (ρ_m/M²) |d ln σ/d ln M| (δ_c/σ) exp(−δ_c²/(2σ²))
  δ_c ≈ 1,686 (поріг сферичного колапсу)
  σ(M) = середньоквадратична флуктуація густини у сфері маси M

Гравітація TreePM:
  Короткодіюча: октодерево Барнса-Хата, підсумовування сили O(N log N)
  Далекодіюча: частинка-сітка (PM) з ШПФ на 3D-сітці
  Довжина згладжування ε ≈ 1–5% від середньої відстані між частинками (уникає двотільного розсіювання)

Симуляція Millennium (Спрінгель, 2005):
  2160³ ≈ 10^10 частинок темної матерії, куб зі стороною 500 Мпк/h
  Illustris-TNG (2018): куб зі стороною 300³ Мпк, повна баріонна фізика (газ, зорі, зворотний зв'язок ЧД)
          

5. Баріонні акустичні осциляції

До рекомбінації баріони і фотони були тісно зв'язані томсонівським розсіюванням, утворюючи фотон-баріонну плазму, що підтримувала звукові хвилі. При рекомбінації фотони відокремилися, і звукові хвилі зупинилися, залишивши заморожений відбиток — масштаб баріонних акустичних осциляцій r_s ≈ 147 Мпк — який діє як «стандартна лінійка» у розподілі галактик.

Масштаб БАО при низькому червоному зміщенні:
  Супутній звуковий горизонт r_s ≈ 147 Мпк  (те саме, що й у РМФ)
  Проявляється як кільце в двоточковій кореляційній функції ξ(r)
  Пік у ξ(r) поблизу r ≈ 150 Мпк (супутня відстань)

Кутова діаметральна відстань D_A(z):
  D_A = (1/(1+z)) ∫_0^z c/H(z') dz'
  БАО вимірює D_A(z)/r_s та H(z)·r_s вздовж/поперек променя зору

Параметр Хаббла:
  H(z)² = H_0² [Ω_m(1+z)³ + Ω_r(1+z)&sup4; + Ω_k(1+z)² + Ω_Λ f_DE(z)]
  Рівняння стану темної енергії: w = p/ρc² (ΛCDM: w = −1)

Результати БАО DESI за 1-й рік (2024):
  Узгоджуються з плоскою ΛCDM з точністю <2%, 11,5 мільйона галактик
  Слабкий натяк: w_0 = −0,55, w_a = −1,32 (динамічна темна енергія, 2,5σ)
  Euclid (запущено 2023): картографує 1,5 × 10^9 галактик до z < 2

Перехресна перевірка баріонної частки:
  f_b = Ω_b / Ω_m ≈ 0,16  (з РМФ + ПНС)
  Спостережено у рентгенівських скупченнях: f_gas ≈ 0,12 (+ зорі ≈ 0,04) → узгоджується
          

6. Космічна павутина

Темна матерія колапсує переважно вздовж шарів і ниток, перш ніж стягуватися у компактні гало на їхніх перетинах — наближення Зельдовича передбачає таку поведінку «млинцювання» на основі власних значень тензора деформації. Галактики формуються всередині гало темної матерії; результуючий візерунок скупчення галактик, спостережений в оглядах, що охоплюють мільярди світлових років, розкриває космічну павутину.

Наближення Зельдовича (1970):
  x(q, t) = q − D(t) ∇ψ(q)
  ψ(q) = гравітаційний потенціал; D(t) = лінійний фактор зростання
  Колапс уздовж трьох осей: спершу вздовж найкоротшого власного вектора (млинці)

Т-класифікація павутини (Хан, 2007):
  Тензор приливних сил T_ij = ∂²φ/∂x_i∂x_j
  # додатних власних значень: 0 → порожнини, 1 → шари, 2 → нитки, 3 → гало

Статистика космічної павутини:
  Гало (вузли):     ~5% об'єму, ~80% маси гало
  Нитки:            ~5% об'єму, ~60% дифузного газу
  Шари:             ~15% об'єму
  Порожнини:        ~75% об'єму, <1% матерії

Баріонний осциляційний спектроскопічний огляд SDSS (BOSS):
  1,3 мільйона галактик, 0,2 < z < 0,75
  Спектр потужності P(k) виміряно з точністю 1% → обмеження на темну енергію

Широкий огляд Euclid (2025–2030):
  1/3 неба, 10^9 джерел, спектроскопічне z для 35 мільйонів
  Ціль: w ±0,01, σ_8 ±0,005
  Перший випуск даних Euclid (DR1): червень 2026 (15 мільйонів галактик)
          

Спробуйте ці симуляції