Астрофізика

Опубліковано: 20 червня 2026 · 15 хв читання · Фізика зір · Політропи · Перенос енергії

Внутрішня структура зір — як зорі утримують власну вагу

Q: Що таке політропна модель зорі?

Політроп — це спрощена модель зорі, що передбачає фіксований зв'язок P = Kρ^γ між тиском і густиною по всій зорі. Індекс n = 1/(γ−1) характеризує модель: n = 3/2 відповідає адіабатичному одноатомному газу (доречно для повністю конвективних зір), n = 3 — зорям, де домінує тиск випромінювання, n = 1.5 — білим карликам у нерелятивістській границі.

Q: Що таке критерій Швацшильда для конвекції?

Критерій Швацшильда стверджує, що шар зорі є нестійким до конвекції, коли фактичний температурний градієнт (∇_rad, встановлений радіаційною дифузією) перевищує адіабатичний температурний градієнт (∇_ad). Формально: конвекція виникає при ∇_rad > ∇_ad. Конвективні зони повністю перемішані, і їх склад гомогенізується швидко порівняно з масштабами зоряної еволюції.

Q: Які чотири рівняння структури зорі?

Чотири рівняння структури зорі: (1) неперервність маси dM/dr = 4πr²ρ, (2) гідростатична рівновага dP/dr = −Gρ M/r², (3) перенос енергії dT/dr = −(3κρ/4acT³)(L/4πr²) для радіаційних зон, (4) генерація енергії dL/dr = 4πr²ρε, де ε — швидкість ядерного генерування енергії на одиницю маси.

Q: Чому Сонце має радіаційну зону й конвективну зону?

У глибокому надрі Сонця (0–0.7 сонячних радіусів) температури надзвичайно високі, непрозорість відносно низька, і випромінювання може ефективно переносити енергію — це радіаційна зона. Вище приблизно 0.7 сонячних радіусів непрозорість різко зростає через часткову іонізацію гелію, температурний градієнт перевищує адіабатичну границю, і конвекція бере верх, утворюючи конвективну оболонку, видиму як грануляція на поверхні Сонця.

Q: Що таке теорема вірілу у фізиці зір?

Теорема вірілу для зорі в гідростатичній рівновазі стверджує, що повна теплова енергія E_th дорівнює −(1/2) від гравітаційної потенціальної енергії E_grav: E_th = −E_grav/2. Важливим наслідком є те, що коли зоря стискається і вивільняє гравітаційну енергію, половина іде на нагрів надр, а інша половина випромінюється.

Q: Як астрономи визначають, що відбувається всередині зорі?

Надра зір неможливо спостерігати безпосередньо, тому астрономи використовують три підходи: (1) чисельні моделі зоряної еволюції, (2) геліосейсмологію та астросейсмологію, що інвертують спостережені частоти коливань для визначення профілів всередині, (3) нейтринні спостереження, що безпосередньо зондують ядро ядерного горіння.

Q: Що відбувається, коли зоря вичерпує водень ядра?

Коли водень ядра вичерпується, ядерне горіння в ядрі припиняється, і воно починає стискатися. Зовнішні шари опускаються всередину, нагріваючи оболонку, поки не спалахує оболонкове горіння. Вивільнена енергія змушує зорю розширюватися в червоного гіганта. Ядро стискається далі, поки не стане білим карликом, нейтронною зорею чи чорною дірою.

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

1. Гідростатична рівновага

Розглянемо тонку сферичну оболонку зоряної речовини на радіусі r, товщиною dr, густиною ρ(r) і площею 4πr². Її вага (гравітаційна сила):

dF_grav = -G M(r) ρ(r) / r² \cdot 4πr² dr

де M(r) = ∫₀ʳ 4πr'² ρ(r') dr' — маса всередині радіуса r. Сила тиску на оболонку назовні — це різниця тиску по її товщині:

dF_pressure = [P(r) - P(r+dr)] \cdot 4πr² = -(dP/dr) dr \cdot 4πr²

Прирівнюючи рівнодійну до нуля (статична рівновага), отримуємо фундаментальне рівняння гідростатичної рівноваги:

dP/dr = -G M(r) ρ(r) / r² = -ρ(r) g(r) де g(r) = G M(r) / r² — місцеве гравітаційне прискорення

Це рівняння є наріжним каменем теорії структури зір. Воно говорить, що тиск мусить спадати назовні — надра зорі перебувають під колосальним тиском, тому що мусять підтримувати вагу всього, що знаходиться вище. У центрі Сонця тиск становить приблизно 2.5 × 10¹⁶ Па — у 250 мільярдів разів більше за атмосферний, — а температура близько 1.5 × 10⁷ К.

Теорема вірілу

Потужним наслідком гідростатичної рівноваги є теорема вірілу. Інтегрування рівняння гідростатики по всій зорі дає зв'язок між повною тепловою енергією E_th і гравітаційною потенціальною енергією E_grav:

2 E_th + E_grav = 0 Отже: E_th = -E_grav/2 і E_total = E_th + E_grav = E_grav/2 < 0

Теорема вірілу має несподіваний наслідок: коли протозоря стискається і втрачає енергію, половина вивільненої гравітаційної енергії нагріває надра, а решта випромінюється. Зорі стають гарячішими, втрачаючи енергію, — вони мають від'ємну теплоємність.

2. Чотири рівняння структури зорі

Повний опис сферично симетричної зорі в тепловій рівновазі потребує чотирьох зв'язаних ОДУ першого порядку, де незалежна змінна — радіус r:

(1) Неперервність маси: dM/dr = 4π r² ρ (2) Гідростатична рівновага: dP/dr = -G M ρ / r² (3) Генерація енергії: dL/dr = 4π r² ρ ε(ρ, T) (4) Перенос енергії: dT/dr = -(3 κ ρ L) / (16π a c T³ r²) [радіаційний] dT/dr = (1 - 1/γ) (T/P) dP/dr [конвективний]

Тут L(r) — світність, ε(ρ,T) — швидкість ядерного генерування енергії на одиницю маси, κ — непрозорість, a — константа випромінювання, c — швидкість світла, γ — адіабатичний показник. Рівняння переносу набуває однієї з двох форм залежно від того, чи сприятливі умови для радіаційного або конвективного переносу.

Щоб замкнути систему, потрібне рівняння стану P(ρ, T) і вирази для κ(ρ, T) і ε(ρ, T). Граничні умови: P = 0, L = 0, ρ = 0 при r = R і M = 0, L = 0 при r = 0. Задавши повну масу M_star, інтегруванням визначається вся структура зорі.

3. Рівняння Лейна-Емдена і політропи

До появи комп'ютерів астрономам були потрібні аналітичні розв'язки. Потужним спрощенням є політропне рівняння стану:

P = K ρ^(1 + 1/n) або рівнозначно P = K ρ^γ де γ = 1 + 1/n

де K — константа, n — політропний індекс. Вводячи безрозмірні змінні:

θ = (ρ/ρ_c)^(1/n) (безрозмірна густина), ξ = r/r_n (безрозмірний радіус) r_n = [ (n+1)K ρ_c^(1/n-1) / (4πG) ]^(1/2) (масштаб довжини)

отримуємо рівняння Лейна-Емдена:

(1/ξ²) d/dξ (ξ² dθ/dξ) = −θⁿ Граничні умови: θ(0) = 1, dθ/dξ|_{ξ=0} = 0

Аналітичні розв'язки існують для трьох особливих випадків:

n = 0 (нестислива рідина): θ = 1 − ξ²/6, ξ₁ = √6 ≈ 2.449. Рівномірна густина скрізь.
n = 1: θ = sin(ξ)/ξ, ξ₁ = π. Профіль густини — функція sinc.
n = 5: θ = (1 + ξ²/3)^(−1/2), ξ₁ → ∞. Нескінченний радіус, але скінченна маса.

Найважливіші в фізичному плані політропи:

n = 3/2 (γ = 5/3): ідеальний одноатомний газ, адіабатичний; описує повністю конвективні зорі і нерелятивістські білі карлики.
n = 3 (γ = 4/3): «стандартна модель» Еддінгтона для масивних зір; також описує релятивістські білі карлики і межу Чандрасекара.
n = 3: маса Чандрасекара M_Ch = 5.83 M_☉/μ_e² випливає з умови скінченної маси і нульового тиску на поверхні.

Стандартна модель Еддінгтона (1926): Еддінгтон показав, що Сонце добре апроксимується n = 3 політропом при домінуванні непрозорості розсіювання на електронах. Це дає центральну температуру ~15 млн К і центральну густину ~150 г/см³ — разюче близько до геліосейсмологічних вимірювань 60 років по тому.

4. Ядерне генерування енергії

Зорі живляться ядерним синтезом у своїх ядрах. Швидкість генерування енергії ε сильно залежить від температури і густини, оскільки синтез вимагає тунелювання крізь кулонівський бар'єр. Для pp-ланцюга (домінуючий процес у Сонці):

Чиста реакція: 4 ¹H \to ⁴He + 2e⁺ + 2ν_e + 26.73 МеВ Швидкість генерування енергії: ε_pp \propto ρ X² T^4 (при сонячних температурах) де X — масова частка водню

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

Цикл CNO

У зірок масивніших за ~1.3 M_☉ (з температурою ядра вище ~2 × 10⁷ К) переважає цикл CNO. Вуглець, азот і кисень є каталізаторами з ε_CNO ∝ ρ X X_CNO T¹⁸. Значно крутіша залежність від температури означає, що зорі CNO мають дуже компактні гарячі ядра і схильні до конвекції в центрі.

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

Випромінювання переносить енергію як випадкове блукання фотонів, кожен з яких поглинається та перевипромінюється незліченну кількість разів. Середня вільна пробіга фотона ℓ = 1/(κρ), де κ — непрозорість Роселанда:

1/κ_R = \int₀^\infty (1/κ_ν) (\partialB_ν/\partialT) dν / \int₀^\infty (\partialB_ν/\partialT) dν

Основні джерела непрозорості:

Розсіювання на електронах (Томсон): κ_es = 0.20(1+X) см²/г, не залежить від T. Переважає в гарячих іонізованих надрах масивних зір.
Вільно-вільне поглинання (гальмівне): κ_ff ∝ ρ T^(−7/2) (непрозорість Крамерса).
Зв'язано-вільне (фотоіонізаційне) поглинання: κ_bf ∝ ρ T^(−7/2). Важливе в зонах часткової іонізації.
Зв'язано-зв'язане (лінійне) поглинання: переважає у фотосфері й холодних зовнішніх шарах.
Непрозорість H⁻: в оболонках нижче ~8000 К негативний іон водню є домінуючим джерелом.

dT/dr = -(3 κ ρ) / (4 a c T³) \cdot L(r) / (4π r²) або в логарифмічній формі: \nabla_rad \equiv (d ln T / d ln P)_rad = (3 κ L P) / (16π a c G M T⁴)

Висока непрозорість κ означає крутий температурний градієнт. Коли він стає надто крутим, рідина стає конвективно нестійкою.

6. Конвективний перенос енергії та критерій Швацшильда

Конвекція виникає, коли фактичний температурний градієнт перевищує те, що підйомний елемент рідини може підтримати адіабатично. Критерій Швацшильда:

\nabla_rad > \nabla_ad де: \nabla_rad = (3 κ L P) / (16π a c G M T⁴) (радіаційний градієнт) \nabla_ad = 1 - 1/γ = 2/5 (для ідеального одноатомного газу з γ = 5/3)

Коли ця умова виконується, елемент рідини, що зміщується вгору, розширюється швидше за навколишнє середовище і виявляється щільнішим — він опускається назад, спричиняючи перекидну циркуляцію. Конвекція є дуже ефективним переносником: фактичний температурний градієнт у конвективній зоні залишається дуже близьким до ∇_ad.

Теорія довжини перемішування

Найбільш широко використовуваною моделлю є теорія довжини перемішування (MLT), введена Людвігом Бірманом і формалізована Ерікою Бьом-Вітензе в 1958 році. Вона розглядає конвекцію як елементи рідини, що рухаються на довжину ℓ = α H_P, де H_P = −dr/d(ln P) — масштаб висоти тиску, α ≈ 1.5–2.0 — вільний параметр:

F_conv = ρ c_P T (\nabla - \nabla_ad) (ℓ/2H_P)² (g δ / 8)^(1/2) (\nabla - \nabla_ad)^(1/2) H_P

Попри простоту, MLT відтворює треки зоряної еволюції, пульсаційні періоди і профілі швидкості звуку з вражаючою точністю.

Конвективне перевищення: реальна конвекція не зупиняється різко на межі Швацшильда. Конвективні потоки проникають трохи в стабільну зону вище, змішуючи речовину ядра в радіаційну зону. Це збільшує ефективний розмір ядра і продовжує час перебування зорі на головній послідовності.

7. Надра Сонця: практичний приклад

Сонце (маса M_☉ = 1.99 × 10³⁰ кг, радіус R_☉ = 6.96 × 10⁸ м, світність L_☉ = 3.85 × 10²⁶ Вт) складається з трьох основних внутрішніх областей:

Ядро (0–0.25 R_☉): температура ~1.5 × 10⁷ К, густина ~150 г/см³. Pp-ланцюг перетворює 600 млн тонн водню на гелій щосекунди. Перенос енергії — радіаційний.
Радіаційна зона (0.25–0.72 R_☉): температура падає від ~7 × 10⁶ К до ~2 × 10⁶ К. Випадкове блукання фотона займає ~170 000 років. Склад майже рівномірний через відсутність перемішування.
Конвективна оболонка (0.72–1.0 R_☉): вище таклоклини непрозорість різко зростає через часткову іонізацію гелію. Критерій Швацшильда порушується, і виникає конвекція, що переносить енергію за ~1 місяць. Конвективні клітини — гранули (1 Мм, 10 хв) і супергранули (30 Мм, 1 день).

⭐

Симулятор зоряної еволюції

Спостерігайте, як зоря старіє від народження до білого карлика або наднової на діаграмі ГР

8. Геліосейсмологія: погляд усередину Сонця

Сонце коливається в мільйонах нормальних мод — акустичні хвилі (p-моди) і поверхневі гравітаційні хвилі (f-моди), захоплені в резонансних порожнинах. Їхні частоти, що спостерігаються як доплерівські зрушення у фотосфері (прилади GOLF, MDI/HMI), несуть інформацію про швидкість звуку:

c_s(r) = \sqrt(γ k_B T / μ m_H) де μ — середня молекулярна маса, m_H — маса атома водню

Інвертуючи спостережувані частоти мод, геліосейсмологи відновлюють профіль швидкості звуку з точністю краще за 0.1%. Результати підтвердили межу радіаційно-конвективної зони на рівні 0.713 R_☉, вміст гелію (Y = 0.249) і відсутність швидко обертального ядра.

Та сама методика для інших зір — астросейсмологія — є нині великою спостережною програмою космічних телескопів Kepler і TESS NASA. Астросейсмологія визначила точні зоряні віки, маси, радіуси і швидкості обертання ядра для десятків тисяч зір.

Проблема сонячних нейтрино — розв'язана: протягом десятиліть виявляли лише ~третину передбаченого нейтринного потоку. У 2001 році Судберійська нейтринна обсерваторія довела, що нейтрино осцилюють між типами під час подорожі до Землі — відкриття, за яке Арт Макдональд отримав Нобелівську премію 2015 року. Сонячні моделі виявилися правильними від самого початку.

Часті запитання

Що таке гідростатична рівновага в зорі?

Гідростатична рівновага — стан, при якому сила градієнта тиску назовні точно врівноважує силу гравітації всередину в кожній точці всередині зорі. Вона описується рівнянням dP/dr = −ρ g(r). Порушення цього балансу змушує зорю розширюватися або стискатися.

Що таке рівняння Лейна-Емдена?

Рівняння Лейна-Емдена — це безрозмірне ОДУ другого порядку, що описує структуру політропної зорі, в якій тиск пропорційний степеню густини (P = Kρ^(1+1/n)). Воно виникає при поєднанні гідростатичної рівноваги з рівнянням Пуассона та має аналітичні розв'язки для n = 0, 1 і 5.

Що таке політропна модель зорі?

Політроп — це спрощена модель зорі з фіксованим зв'язком P = Kρ^γ між тиском і густиною. Індекс n = 1/(γ−1): n = 3/2 відповідає адіабатичному одноатомному газу, n = 3 — зорям з переважанням тиску випромінювання, n = 1.5 — нерелятивістським білим карликам.

Що таке непрозорість і чому вона важлива для структури зір?

Непрозорість κ вимірює, наскільки сильно зоряна речовина поглинає або розсіює випромінювання на одиницю маси. Висока непрозорість перешкоджає потоку радіаційної енергії, посилюючи температурний градієнт. Якщо він стає надто крутим, запускається конвекція.

Що таке критерій Швацшильда для конвекції?

Критерій Швацшильда стверджує, що шар зорі нестійкий до конвекції, коли ∇_rad > ∇_ad. Конвективні зони повністю перемішані і їх склад гомогенізується швидко порівняно з масштабами зоряної еволюції.

Які чотири рівняння структури зорі?

Чотири рівняння: (1) неперервність маси dM/dr = 4πr²ρ, (2) гідростатична рівновага dP/dr = −Gρ M/r², (3) перенос енергії dT/dr = −(3κρ/4acT³)(L/4πr²) для радіаційних зон, (4) генерація енергії dL/dr = 4πr²ρε.

Чому Сонце має радіаційну зону й конвективну зону?

У глибокому надрі Сонця (0–0.7 R_☉) радіація ефективно переносить енергію — це радіаційна зона. Вище ~0.7 R_☉ непрозорість різко зростає через часткову іонізацію гелію, і конвекція бере верх, утворюючи видиму як грануляцію конвективну оболонку.

Що таке теорема вірілу у фізиці зір?

Теорема вірілу стверджує: E_th = −E_grav/2. Коли зоря стискається, половина вивільненої гравітаційної енергії нагріває надра, а інша половина випромінюється.

Як астрономи визначають, що відбувається всередині зорі?

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

Що відбувається, коли зоря вичерпує водень ядра?

Ядерне горіння припиняється, і ядро стискається. Зовнішні шари нагрівають оболонку, поки не спалахує оболонкове горіння. Зоря розширюється в червоного гіганта. Ядро продовжує стискатися, поки не стане білим карликом, нейтронною зорею чи чорною дірою.