Struktura wnętrza gwiazd: Jak gwiazdy utrzymują swoją masę…

Równowaga hydrostatyczna: ciśnienie przeciwko mocy grawitacyjnej

Załóżmy, że małą sfericzną powłokę materiały gwiazdowej znajdujemy na odległości r. Sile grawitacji zacina się wewnątrz; różnica ciśnienia po obu stronach powłoki napędza ją zewnątrz. Ustalanie równości tych dwóch sił daje jedno równanie, które jest podstawą całej teorii struktury gwiazd:

dP/dr = −G·M(r)·ρ(r) / r² (równowaga hydrostatyczna) Środek Słońca: P ≈ 2,5×10¹⁶ Pa (≈250 miliarda atmosfer) T ≈ 1,5×10⁷ K Ciśnienie musi maleć zewnątrz — wnętrze gwiazdy jest naciskane do tego ciśnienia dokładnie dlatego, że musi wspierać wszystko, co znajduje się nad nim. Konsekwencją tej równowagi, twierdzenia virialnego, jest stwierdzenie, że 2E_th + E_grav = 0: gdy proto-gwiazda kontrakuje i wydziela energię, połowa wydzielonej energii grawitacyjnej ciepli wnętrze, a druga połowa ucieka jako światło. Gwiazdy stają się gorętsze, gdy stracą energię — mają negatywną pojemność ciepłową.

dP/dr = −G·M(r)·ρ(r) / r²        (hydrostatic equilibrium)

Sun's centre:  P ≈ 2.5×10¹⁶ Pa  (≈250 billion atmospheres)
               T ≈ 1.5×10⁷ K

Cztery równania, jedna gwiazda

Aby stworzyć kompleksowy sferycznie symetryczny model gwiazdy, potrzebujemy czterech powiązanych równań pierwszego rzędu: kontinuowanie masy (dM/dr = 4πr²ρ), hydrostatyczna równowaga, generacja energii (dL/dr = 4πr²ρε, wynikająca z spalenia jądrowego) i transport energii, który może być radiaacyjny lub konwective w zależności od lokalnych warunków. Aby zamknąć system, potrzebujemy równania stanu P(ρ,T) oraz przepisów dotyczących opakowania κ(ρ,T) i stopnia generacji energii jądrowej ε(ρ,T). Warunki brzegowe wynoszą: P = L = ρ = 0 na powierzchni, a M = L = 0 w centrum.

demo na żywo · powiązana symulacja● LIVE

Równanie Lane-Emdena i poltropoide

Przed komputerami astronomiowie potrzebowali skrótów analitycznych. Przyjmując równanie stanu poltropiczne P = Kρ^(1+1/n), cztery równania struktury zmniejszają się do jednego wymiarowego drugiego rzędu równania różniczkowego bezwymiarowego, zwane równaniem Lane-Emdena. To równanie rozwiązano od θ(0)=1 w stronę zewnętrznej powierzchni gwiazdy, aż do chwili, gdy bezwymiarowa gęstość θ pierwszy raz osiąga zero na powierzchni gwiazdowej.

θ″ + (2/ξ)·θ′ + θⁿ = 0 θ(0)=1, θ′(0)=0 n=0 (niepopychalny): θ = 1 − ξ²/6, ξ₁ = √6 ≈ 2.449 n=1: θ = sin(ξ)/ξ, ξ₁ = π n=5: θ = (1+ξ²/3)^(−1/2), ξ₁ → ∞ (promień nieskończony, masa skończona) n = 3/2 (γ = 5/3) opisuje gwiazdy convekcyjne pełne i nie-relatywistyczne chmury białe; n = 3 (γ = 4/3) jest model standardowy Eddingtona dla masywnych gwiazd dominujących w presji radiacyjnej, a ten sam indeks oznacza granicę mas Chandrasekhar M_Ch = 5.83 M☉/μₑ² dla relatywistycznych chmur białych. Eddington pokazał w 1926 roku, że Słońce same jest dobrze przybliżone przez n = 3, przewidując temperaturę centralną blisko 15 milionów K i gęstość centralną blisko 150 g/cm³ — te wartości potem były potwierdzone precyzyjnie seismologią słoneczną.

θ″ + (2/ξ)·θ′ + θⁿ = 0        θ(0)=1,  θ′(0)=0

n=0 (incompressible):  θ = 1 − ξ²/6,        ξ₁ = √6 ≈ 2.449
n=1:                    θ = sin(ξ)/ξ,        ξ₁ = π
n=5:                    θ = (1+ξ²/3)^(−1/2), ξ₁ → ∞ (infinite radius, finite mass)

Odtwarzalność, conwecja i kryterium Schwarzschilda

Radiacja przekazuje energię w dół jako fotony, które poruszają się losowo przez plasma, odbijając się od elektronów i podlegając absorbowaniu i ponownemu emitowaniu przez iony — podróż z jądra Słońca do podstawy swojej strefy conwecji trwa około 170 000 lat. Gdy odtwarzalność κ jest za wysoka, gradient temperatury radiacyjnej potrzebny do nacisku na wymagane luminosyczności w dół staje się bardziej nachylony niż własna stopa chłodzenia adiabatycznego pakiety gazowego — kryterium Schwarzschilda, ∇_rad > ∇_ad — i conwecja zastępuje to. W Słońcu to dzieje się powyżej około 0,72 R☉, gdzie częściowe ionizowanie helu spowoduje wzrost odtwarzalności; wynikająca warstwa conwective przekształca energię resztkę drogi do powierzchni w ciągu około miesiąca, a nie 170 000 lat, tworząc granulację widoczne na powierzchni Słońca.

Heliosejsyometria: dowód z wnętrza Słońca

Słońce oscyluje w milionach trybów akustycznych, których częstotliwości, pomiarowane jako przesunięcia Dopplera na fotosferze, zawierają informacje o profilu prędkości dźwiękowej na każdej głębokości. Odwrócenie tych częstotliwości pozwala wskazać wnętrze z dokładnością lepszą niż 0,1%, potwierdzając granicę radiacyjno-konwektorską dokładnie w 0,713 R☉. Taka sama technika — heliosejsyometria — teraz stosuje się do tysięcy innych gwiazd poprzez misje Kepler i TESS. Działa nawet na rozwiązanie problemu neutrino z lat: pomiar znajdował tylko trzy czwarte przewidywanego fluxu neutrino do 2001 roku, kiedy Obserwatorium Neutrino w Sudbury udowodniło, że elektronowe neutrino przeobrażają się na inne typy podczas podróży do Ziemi — praca, która zdobyła Artowi McDonaldowi Nobla w 2015 roku.

Często zadawane pytania

Czym jest równowaga hydrostatyczna w gwiazdzie?

Równowaga hydrostatyczna to stan, w którym siła prędkościowa nacisku wyjściowego dokładnie równa się siłie przyciągania grawitacyjnej na każdym punkcie wewnątrz gwiazdy, opisanej wzorem dP/dr = −G·M(r)·ρ/r². W centrum Słońca ta równowaga wymaga nacisku około 2,5×10¹⁶ Pa, czyli około 250 miliarda razy ciśnienie atmosferyczne, oraz temperatury blisko 1,5×10⁷ K. Przerwanie tej równowagi spowoduje rozszerzenie się lub zwiędnięcie gwiazdy.

Czym jest równanie Lane-Emdena i co to jest politropo?

Politropo przyjmuje uproszczone równanie stanu P = Kρ^(1+1/n), gdzie n to indeks politropowy. Połączenie tego z równaniem równowagi hydrostatycznej oraz równaniem Poissona dla grawitacji, a następnie przejście do wymiarowej promienia ξ i gęstości θ, daje równanie Lane-Emdena θ″ + (2/ξ)·θ′ + θⁿ = 0. Ma ono dokładne rozwiązania analityczne tylko dla n=0, 1 i 5; każdy inny indeks, w tym ważny z punktu widzenia astronomii n=3/2 i n=3, wymaga numerycznej integracji.

Dlaczego politropo o indeksie n=3 nazywane jest modeliem standardowym Eddingtona?

Arthur Eddington pokazał w 1926 roku, że gwiazda dominująca przez przepuszczalność opakowania elektronowego, z luminosowością proporcjonalną do masy, jest dobrze przybliżona przez politropo o indeksie n=3. Zastosowanie tego do Słońca przewiduje temperaturę centralną blisko 15 milionów K i gęstość centralna około 150 g/cm³ — wartości, które heliosejsjometria potwierdziła z bardzo dużą dokładnością sześć dekad później. Ten sam indeks n=3 opisuje również gwiazdy o masie znacznym dominowane przez ciśnienie promieniowe oraz białe gwiazdy dżetne bliskie granicy Chandrasekhar'a.

Wypróbuj na żywo

Wszystko powyżej działa bezpośrednio w Twojej przeglądarce — otwórz Stellar Interior Structure i zmieniaj parametry podczas działania. Nic nie jest instalowane ani przesyłane na serwer, cały model działa w jednej karcie.

▶ Otwórz symulację Stellar Interior Structure

Struktura wnętrza gwiazd: Jak gwiazdy utrzymują swoją masę pod wpływem grawitacji