Спотлайт #32 – Астрофізика та зоряна еволюція: діаграма Герцшпрунга-Рассела, будова зір, наднові та гравітаційні хвилі

Кожен атом важчий за літій було вироблено всередині зорі, а більшість водню й гелію у всесвіті було виковано впродовж перших трьох хвилин після Великого вибуху. Зоряна астрофізика перебуває на перетині ядерної фізики, гідродинаміки, перенесення випромінювання та загальної теорії відносності. Найпростіші зорі — об'єкти головної послідовності в гідростатичній рівновазі — добре описуються жменькою рівнянь. Складні ж — зорі асимптотичної гілки гігантів із багатьма оболонками горіння, масивні зорі, що обертаються та втрачають масу з вітрами, компактні подвійні системи, що зближуються до злиття — належать до найактивніших фронтирів сучасної фізики.

1. Діаграма Герцшпрунга-Рассела

Коли Ейнар Герцшпрунг та Генрі Норріс Рассел незалежно один від одного побудували графік світності зір проти температури поверхні близько 1910 року, вони виявили, що зорі розподілені не випадково: переважна більшість групується вздовж діагональної смуги (головна послідовність), із гілками гігантів і надгігантів вище та білими карликами нижче. Діаграма ГР — найкорисніша діаграма в зоряній астрофізиці, бо кожна еволюційна стадія займає окрему область.

Співвідношення Стефана-Больцмана:
  L = 4πR² σ T_eff&sup4;

Співвідношення маса-світність для головної послідовності (наближено):
  L / L_сонця ≈ (M / M_сонця)^4      для 0,5 < M/M_сонця < 50
  L / L_сонця ≈ (M / M_сонця)^2,5    для зір меншої маси

Тривалість перебування на головній послідовності:
  t_ГП ≈ t_сонця · (M/M_сонця)^{-2,5}    (сонячне t_ГП ≈ 10 млрд років)
  1 M_сонця: ~10 млрд років     5 M_сонця: ~100 млн років     25 M_сонця: ~7 млн років

Області діаграми ГР (L у L_сонця, T у К):
  Головна послідовність:  0,001–10^6 L_сонця, 2500–50 000 К
  Червоні гіганти:    10–1000 L_сонця, 3500–5500 К
  Надгіганти:   10^4–10^6 L_сонця, 3500–30 000 К
  Білі карлики:  10^-4–0,1 L_сонця, 5000–100 000 К

Спектральні класи (OBAFGKM, від найгарячіших до найхолодніших):
  O: >30 000 К   B: 10 000–30 000 К   A: 7500–10 000 К
  F: 6000–7500 К   G: 5200–6000 К   K: 3700–5200 К
  M: <3700 К (холодний червоний карлик; найчисленніший тип у Галактиці)

Позиція зорі на діаграмі ГР змінюється в міру її еволюції. Сонце покине головну послідовність приблизно через 5 млрд років, розшириться до червоного гіганта (R ~200 R_сонця), скине зовнішні шари у вигляді планетарної туманності й залишить по собі трек охолодження білого карлика. Масивні зорі класу O еволюціонують настільки швидко, що ніколи не перебувають у рівновазі з навколишнім газовим тиском: вони гинуть у вибухах наднових протягом 3–10 млн років після народження. Зоряні скупчення особливо корисні, оскільки всі їхні члени мають однаковий вік; точка сходу з головної послідовності — найяскравіша зоря головної послідовності, що ще спалює водень — дає безпосередньо вік скупчення.

2. Внутрішня будова зір та політропи

Зоря — це самогравітуюча куля газу в гідростатичній рівновазі: гравітація, що тягне всередину, і градієнт тиску, що штовхає назовні, врівноважуються на кожному шарі. Три додаткові рівняння описують генерацію енергії, перенесення енергії та зміну температури з тиском. Разом вони утворюють рівняння зоряної будови, які однозначно визначають внутрішній профіль зорі, коли задано її масу й склад (теорема Фогта-Рассела).

Гідростатична рівновага:
  dP/dr = −ρ G M(r) / r²

Неперервність маси:
  dM/dr = 4πr² ρ

Перенесення енергії (радіаційна зона):
  dT/dr = −3κρL / (64πσr²T³)
  κ = середня непрозорість Росселанда (см²/г)

Критерій конвекції (Шварцшильд):
  |dT/dr|_рад > |dT/dr|_ад   → конвективна нестійкість
  Теорія довжини перемішування: масштаб конвективної комірки ≈ α_MLT H_P  (α ~ 1,5-2)

Політропне рівняння стану:
  P = K ρ^{(n+1)/n}   (політропний показник n)
  n = 3/2 → повністю конвективна (маломасивні зорі, ядра гігантів)
  n = 3   → домінує променевий тиск (масивні зорі; Чандрасекар)

Рівняння Лейна-Емдена (безрозмірна зоряна будова):
  d/dξ[ξ² dθ/dξ] + ξ² θ^n = 0
  ξ = r / r_n (масштабований радіус),  θ = (ρ/ρ_c)^{1/n}

Надра Сонця (модель AGSS09):
  Ядро (r < 0,25 R_сонця):  T_c = 1,57×10^7 К,  ρ_c = 151 г/см³
  Радіаційна зона: 0,25–0,71 R_сонця
  Конвективна зона: 0,71–1,00 R_сонця  (глибина ~200 000 км)

Геліосейсмологія — вимірювання швидкості звукових хвиль шляхом спостереження поверхневих осциляцій — підтвердила будову надр Сонця з точністю менше 1%. Розбіжність між 3D-спектроскопічним вмістом важких елементів на Сонці (AGSS09) та геліосейсмічними обмеженнями залишається відкритою проблемою зоряної фізики — «проблемою сонячної металічності».

3. Цикли ядерного горіння

Джерело енергії, що дозволяє зорям протистояти гравітаційному колапсу, — ядерний синтез. Маломасивні зорі (M < 1,3 M_сонця) синтезують водень переважно через протон-протонний (pp) ланцюжок; масивніші зорі використовують CNO-цикл, швидкість якого дуже чутлива до температури (T^17 проти T^4 для pp). Із віком зоря нарощує гелієве ядро попелу, потім запалює гелій через потрійний альфа-процес, і — для масивних зір — продовжує через горіння вуглецю, неону, кисню й кремнію на дедалі коротших часових масштабах.

Протон-протонний ланцюжок (домінує для M ≤ 1,3 M_сонця):
  pp I:  4H → ⁴He + 2e’ + 2ν_e       (Q = 26,7 МеВ)
  pp II: гілка ⁷Be (важлива при T > 14 МК)
  pp III: гілка ⁸B (виробляє високоенергетичні нейтрино, виявлені SNO)
  Генерація енергії: ε_pp ∝ ρ X² T^4  (X = масова частка водню)

CNO-цикл (домінує для M > 1,3 M_сонця):
  ¹²C + p → ¹³N + γ → ¹³C + e’ + ν
  ¹³C + p → ¹⁴N + γ
  ¹⁴N + p → ¹⁵O + γ  ← етап, що лімітує швидкість
  ¹⁵O → ¹⁵N + e’ + ν → + p → ¹²C + ⁴He (сумарно)
  Генерація енергії: ε_CNO ∝ ρ X_CNO T^17

Потрійний альфа-процес (горіння гелію, T > 10^8 К):
  2 ⁴He → ⁸Be* (час життя резонансу 2,6×10^-16 с)
  ⁸Be + ⁴He → ¹²C* (стан Гойла, 7,6644 МеВ)  → ¹²C + 2γ
  ¹²C + ⁴He → ¹⁶O + γ  (конкурує; визначає відношення C/O)
  Резонанс Гойла передбачений Гойлом (1954), підтверджений Куком та ін.

Просунуті стадії горіння (масивні зорі, T в одиницях 10^8 К):
  Горіння He:   T ~2,   тривалість ~10^6 років   (>8 M_сонця)
  Горіння C:    T ~5,   тривалість ~10^3 років
  Горіння Ne:   T ~15,  тривалість ~1 рік
  Горіння O:    T ~20,  тривалість ~місяці
  Горіння Si:   T ~35,  тривалість ~дні
  Залізне ядро:    T ~70,  немає вивільнення енергії → колапс ядра
  Ядра групи заліза — кінцева точка: енергія зв'язку досягає піку на ⁵⁶Fe/⁵⁸Ni

4. Наднові з колапсом ядра

Коли масивна зоря (>8 M_сонця) накопичує залізне ядро, масивніше за межу Чандрасекара (~1,4 M_сонця), тиск виродження електронів більше не може її підтримувати. Ядро колапсує від розміру приблизно з Землю до нейтронної зорі радіусом 10 км менш ніж за секунду, вивільняючи ~3 × 10^53 ерг — гравітаційну енергію зв'язку нейтронної зорі. Це приблизно у 100 разів більше за повну електромагнітну енергію, випромінену Сонцем за весь його 10-мільярдний рік життя, вивільнену за мілісекунди.

Маса Чандрасекара (релятивістська межа виродження електронів):
  M_Ch = 5,83 Y_e² M_сонця   (Y_e = електронна частка ≈ 0,5 для ⁵⁶Fe)
  M_Ch ≈ 1,44 M_сонця    (фактичне значення ~1,2-1,4 M_сонця з поправками)

Часова шкала колапсу:
  t = 0 мс:    Залізне ядро перевищує M_Ch; починається вільне падіння
  t ~ 100 мс:  Ядро досягає ядерної густини (ρ_ядра = 2,7×10^14 г/см³)
  t ~ 110 мс:  Ядро відскакує; запускається ударна хвиля зі швидкістю ~100-200 км/с
  t ~ 200 мс:  Ударна хвиля зупиняється — фотодисоціація заліза коштує 8,8 МеВ/нуклон
  t ~ 500 мс:  Нейтринний нагрів повторно надає енергію ударній хвилі (механізм BNNS)
  t ~ 1 с:     Успішний вибух; ударна хвиля прориває поверхню через ~15-20 год
  t ~ дні:     Оптичний максимум

Нейтринний спалах (підтверджено для SN 1987A):
  Повна E_ν: ~3 × 10^53 ерг  (99% енергії колапсу)
  Тривалість: ~10 с (кельвін-гельмгольцівське охолодження нейтронної зорі)
  Kamiokande II + IMB + Baksan зареєстрували 24 нейтрино від SN 1987A (ВМХ, 160 тис. св. років)

Нуклеосинтез під час вибухової фази:
  r-процес (швидке захоплення нейтронів): виробляє ~50% ядер з A>100
  Нещодавній доказ: r-процес у кілоновій злиття НЗ GW170817 (золото, платина, європій)
  p-процес: фотодезінтеграція зернистих ядер; протон-збагачені ізотопи
  ν p-процес: нейтрино-керуване захоплення протонів у ранньому викиді

Пікова світність наднових з колапсом ядра:
  L_пік ~ 10^43 ерг/с (оптичний діапазон)
  Спад живиться ⁵⁶Ni → ⁵⁶Co → ⁵⁶Fe (t_½ = 6,08 доби, 77,2 доби)
  Яскравість ∝ синтезованій M_Ni (~0,07 M_сонця для SN 1987A)

5. Нейтронні зорі та пульсари

Реліктове ядро, що залишається після наднової з колапсом ядра — це нейтронна зоря — сфера масою ~1,4 M_сонця, стиснута до радіуса ~10 км, яку підтримує тиск виродження нейтронів та відштовхувальне жорстке ядро ядерної сили. Середня густина перевищує ядерну; чайна ложка важила б ≈ 5 × 10^8 тонн. Якщо ядро зорі-попередника перевищує ~2–3 M_сонця, колапс продовжується за межі стійкості нейтронної зорі до чорної діри.

Рівняння гідростатичної рівноваги Толмена-Оппенгеймера-Волкова (TOV) у ЗТВ:
  dP/dr = −G(M + 4πr³P/c²)(ρ + P/c²) / [r(r − 2GM/c²)]
  Зводиться до ньютонівського dP/dr = −ρGM/r² при низькій компактності

Максимальна маса нейтронної зорі:
  Межа Толмена-Оппенгеймера-Волкова ~0,7 M_сонця (чиста нейтронна речовина, м'яке РС)
  Реальна межа (жорстке ядерне РС): ~2,3–2,5 M_сонця
  Виміряно: PSR J0952-0607 = 2,35 ± 0,17 M_сонця (2022, найважча підтверджена НЗ)

Уповільнення обертання пульсара (магнітно-дипольне випромінювання):
  P_dot = −B²R^6Ω^3sin²α / (6Ic³)
  P = період обертання (найшвидший відомий: J1748-2446ad, P = 1,396 мс)
  Характерний вік: τ = P / (2 P_dot)
  Поверхневе магнітне поле:    B ~ 3,2×10^19 √(P P_dot)  Гаусс

Охолодження нейтронної зорі:
  t < 100 років:  Прямий Урка (якщо M > M_DU), швидке охолодження, T_s ~ 10^6 К
  t ~ 10^3 років: Домінує модифікований Урка, T_s ~ 5×10^5 К
  t ~ 10^5 років: Фотонне охолодження, T_s ~ 10^5 К
  Спостережені приклади: Кассіопея A (350 років) показала охолодження ~2% за десятиліття (Хо і Хайнке, 2009)

Рівняння стану (РС) щільної речовини:
  Адронне: SLy, APR4 (ядерні + тритільні сили)
  Екзотичне: гіперонне, конденсат каонів, кварк-глюонна плазма
  Обмеження радіуса (NICER): R = 12,35 ± 0,75 км для J0030+0451 (Райлі, 2019)

6. Гравітаційні хвилі від зближення компактних подвійних систем

Загальна теорія відносності передбачає, що прискорювані маси випромінюють гравітаційні хвилі — брижі кривизни простору-часу, що поширюються зі швидкістю світла. Подвійні системи компактних об'єктів (нейтронних зір або чорних дір) втрачають енергію на гравітаційне випромінювання, через що орбіта звужується, доки два тіла не зіллються. LIGO/Virgo вперше зареєстрували цей сигнал 14 вересня 2015 року (GW150914): дві чорні діри масою 36 і 29 M_сонця, що злилися на відстані ~410 Мпк, породивши сигнал «дзвону» тривалістю ~0,2 секунди.

Квадрупольна формула потужності ГХ (провідний порядок ПН):
  P_ГХ = 32G^4 m_1²m_2²(m_1+m_2) / (5c^5 a^5)
  a = орбітальна відстань

Чирп-маса ℳ (визначає еволюцію фази):
  ℳ = (m_1 m_2)^{3/5} / (m_1+m_2)^{1/5}
  Безпосередньо вимірюється зі свіп-частоти f_ГХ(t)

Орбітальний розпад (Пітерс, 1964):
  da/dt = −64G³m_1m_2(m_1+m_2) / (5c^5 a³)
  Час до злиття: t_злиття = 12a_0^4c^5 / (19G³ m_1m_2(m_1+m_2))
  PSR 1913+16 (пульсар Халса-Тейлора): t_злиття ~ 300 млн років; орбітальний розпад узгоджується із ЗТВ з точністю 0,2%

Деформація гравітаційної хвилі:
  h_+ = − (2Gℳ/c²r)(Gπfℳ/c^3)^{2/3} cos(2φ)
  h_× тієї ж амплітуди, зі зсувом фази 90°

GW150914 (перша реєстрація, вересень 2015):
  Маси компонентів: 35,6 M_сонця + 30,6 M_сонця
  Маса злиття: 63,1 M_сонця  (випромінено ΔE = 3,0 M_сонця c² ~ 5,4×10^54 ерг)
  Пікова деформація: h ~ 10^-21   (довжина плеча LIGO 4 км; ΔL ~ 4×10^-18 м, 1/1000 протона)
  Пікова світність ГХ: ~3,6×10^56 ерг/с (переважає всі видимі зорі у видимому всесвіті)

GW170817 (перше злиття НЗ-НЗ, серпень 2017):
  Маси компонентів: 1,17 + 1,36 M_сонця
  Пов'язана кілонова AT2017gfo: блакитний (Sr) + червоний (важкий r-процес) компоненти
  Гамма-спалах GRB 170817A зі зсувом 1,74 с → c_ГХ/c = 1 ± 10^-15 (перевіряє лоренц-інваріантність)
  Стала Хаббла: H_0 = 70,0^{+12}_{-8} км/с/Мпк (мультимесенджерна)

Станом на 2024 рік мережа LIGO-Virgo-KAGRA каталогізувала понад 90 кандидатів злиття компактних подвійних систем за три спостережні прогони (GWTC-3). Четвертий спостережний прогін (O4, 2023–2025) має подвоїти цю вибірку. Наступне покоління детекторів — Einstein Telescope (підземний, плечі 10 км) та Cosmic Explorer (плечі 40 км) — реєструватиме злиття подвійних чорних дір по всьому доступному спостереженню всесвіту, уможливлюючи пряме вимірювання історії космічного розширення без електромагнітних спостережень.

Пов'язані симуляції