Astrofizyka – Struktura Galaktyki – 2025 r.

Masa Drogi Mlecznej: Składniki, równania i otwarte problemy

Kompletne zestawienie głównych składników masy naszej Galaktyki – od dysków gwiezdnych po centralną czarną dziurę – z radialnymi równaniami masy, symulacjami wizualnymi i otwartymi pytaniami, które pozostają nierozwiązane.

Na podstawie McMillan 2017 – Ou et al. 2024 – Bland-Hawthorn & Gerhard 2016

Spis treści

  1. Cienki dysk gwiezdny
  2. Gruby dysk gwiezdny
  3. Gaz atomowy HI
  4. Gaz cząsteczkowy H₂
  5. Wybrzuszenie i pasek
  6. Centralna czarna dziura Sagittarius A*
  7. Gwiezdna aureola
  8. Całkowita widoczna masa
  9. Brakująca masa
  10. Symulacja promieniowego profilu masy
  11. Otwarte problemy

Droga Mleczna to nasza rodzima galaktyka: spirala z poprzeczką zawierająca około stu miliardów gwiazd, duży dysk gazowy, gwiezdne halo i centralną supermasywną czarną dziurę. Pomimo tego, że jest to najlepiej zbadana galaktyka we wszechświecie, nadal pozostają fundamentalne pytania dotyczące jej całkowitej masy, zewnętrznego halo i niewidocznej masy wymaganej przez jej krzywą rotacji.

Wszystkie poniższe masy są wyrażone jako radialne masy skumulowane: całkowita masa zawarta w promieniu r od Centrum Galaktyki.

\(M(<r)\)

Jest to naturalna obserwowalna wielkość, ponieważ określa prędkość kołową zgodnie z prawem Newtona:

\(V_c(r)=\sqrt{\frac{G\,M(<r)}{r}}\) \(G=4.302\times10^{-6}\,\mathrm{kpc\,km^2\,s^{-2}\,M_\odot^{-1}}\)

1. Cienki dysk gwiezdny

Składnik 1 – Cienki dysk gwiezdny – M ≈ 3.52 × 10¹⁰ M⊙

Cienki dysk jest dominującym składnikiem gwiezdnym Drogi Mlecznej. Zawiera Słońce, ramiona spiralne, gwiazdy młode i w średnim wieku, większość gazu i pyłu międzygwiazdowego oraz główne miejsca formowania się gwiazd. Jej pionowa grubość jest niewielka w porównaniu z jej rozległością radialną.

Gęstość powierzchniowa jest modelowana jako dysk wykładniczy:

\(\Sigma_{\mathrm{thin}}(r)=\Sigma_{0,\mathrm{thin}}e^{-r/R_{d,\mathrm{thin}}}\)
ParametrSymbolWartośćŹródło
Gęstość powierzchni centralnejΣ0,thin896 M⊙ pc-²McMillan 2017
Promień skaliRd, cienki2,50 kpcMcMillan 2017
Masa całkowitaMthin3.52 × 10¹⁰ M⊙Od 2πΣ₀Rd²

Promieniowa masa skumulowana wynosi:

\(M_{\mathrm{thin}}(<r)=3.52\times10^{10}\left[1-e^{-r/2.50}\left(1+\frac{r}{2.50}\right)\right]M_\odot\)

Wzór ten pochodzi z całkowania gęstości powierzchniowej po pierścieniach kołowych. Masa cienkiego dysku gwałtownie wzrasta wewnątrz kilku kiloparseków, a następnie nasyca się w kierunku masy całkowitej.

2. Gruby dysk gwiezdny

Składnik 2 – Gruby dysk gwiezdny – M ≈ 1.05 × 10¹⁰ M⊙

Gruby dysk to starsza, bardziej rozproszona populacja gwiazd, która rozciąga się dalej powyżej i poniżej płaszczyzny Galaktyki. Gwiazdy w nim zawarte mają inną metaliczność i kinematykę niż gwiazdy w dysku cienkim i mogą być świadkami wcześniejszych fuzji lub wydarzeń związanych z ogrzewaniem w Drodze Mlecznej.

\(\Sigma_{\mathrm{thick}}(r)=\Sigma_{0,\mathrm{thick}}e^{-r/R_{d,\mathrm{thick}}}\)
ParametrSymbolWartość
Gęstość powierzchni centralnejΣ0,thick183 M⊙ pc-²
Promień skaliRd, grubość3.02 kpc
Masa całkowitaMthick1.05 × 10¹⁰ M⊙
\(M_{\mathrm{thick}}(<r)=1.05\times10^{10}\left[1-e^{-r/3.02}\left(1+\frac{r}{3.02}\right)\right]M_\odot\)

Łączna masa dysku gwiezdnego wynosi:

\(M_{\mathrm{disk,\star}}(<r)=M_{\mathrm{thin}}(<r)+M_{\mathrm{thick}}(<r)\) \(M_{\mathrm{disk,\star,total}}\approx4.57\times10^{10}M_\odot\)

3. Gaz atomowy – HI

Składnik 3 – Atomowy wodór – M ≈ 1,1 × 10¹⁰ M⊙

Linia radiowa 21 cm neutralnego wodoru wskazuje na duży, rozszerzony i wypaczony dysk gazowy rozciągający się daleko poza dysk gwiazdowy. W przeciwieństwie do gwiazd, HI ma centralną depresję i osiąga maksima kilka kiloparseków od Centrum Galaktyki.

\(\Sigma_{\mathrm{HI}}(r)=\Sigma_{0,\mathrm{HI}}\exp\left(-\frac{R_{m,\mathrm{HI}}}{r}-\frac{r}{R_{d,\mathrm{HI}}}\right)\)
ParametrWartośćZnaczenie
Rm,HI4.0 kpcTworzy centralny otwór
Rd,HI7.0 kpcZewnętrzna skala wykładnicza
MHI, łącznie1.1 × 10¹⁰ M⊙Całkowita masa atomowa gazu
\(M_{\mathrm{HI}}(<r)=1.1\times10^{10}\frac{\int_0^r e^{-4/R-R/7}R\,dR}{\int_0^\infty e^{-4/R-R/7}R\,dR}M_\odot\)

Szczyt rozkładu masy HI znajduje się w pobliżu r ≈ √(4 × 7) ≈ 5,3 kpc. HI jest ważny zarówno jako rezerwuar gazu, jak i wskaźnik zewnętrznego potencjału galaktycznego.

4. Gaz cząsteczkowy – H₂

Składnik 4 – Wodór cząsteczkowy – M ≈ 1.2 × 10⁹ M⊙

Wodór molekularny jest skoncentrowany w wewnętrznej Galaktyce i jest ściśle związany z gigantycznymi obłokami molekularnymi i formowaniem się gwiazd. Jest on zazwyczaj śledzony poprzez emisję CO, która wprowadza niepewność poprzez współczynnik konwersji CO na H₂.

\(\Sigma_{\mathrm{H_2}}(r)=\Sigma_{0,\mathrm{H_2}}\exp\left(-\frac{R_{m,\mathrm{H_2}}}{r}-\frac{r}{R_{d,\mathrm{H_2}}}\right)\)
ParametrWartość
Rm,H₂12.0 kpc
Rd,H₂1,5 kpc
MH₂, łącznie1.2 × 10⁹ M⊙
\(M_{\mathrm{H_2}}(<r)=1.2\times10^9\frac{\int_0^r e^{-12/R-R/1.5}R\,dR}{\int_0^\infty e^{-12/R-R/1.5}R\,dR}M_\odot\)

5. Wybrzuszenie i pręt

Składnik 5 – wybrzuszenie centralne i poprzeczka galaktyczna – M ≈ 9.23 × 10⁹ M⊙

Droga Mleczna jest galaktyką spiralną z poprzeczką. Jej centralne wybrzuszenie i poprzeczka zawierają stare gwiazdy i silnie wpływają na przepływ gazu i dynamikę gwiazd w wewnętrznej Galaktyce. Pasek jest trudny do zmierzenia z naszej pozycji wewnątrz dysku, co sprawia, że wewnętrzny rozkład masy jest niepewny.

\(\rho_{\mathrm{bulge}}(r)\propto e^{-(r/r_b)^2}\) \(r_b\approx0.5\,\mathrm{kpc}\)

Użytecznym przybliżeniem sferycznym dla masy skumulowanej jest:

\(M_{\mathrm{bulge}}(<r)\approx9.23\times10^9\left[1-e^{-r}\left(1+r+\frac{r^2}{2}\right)\right]M_\odot\)

Prawie cała masa wybrzuszenia znajduje się w odległości kilku kiloparseków. Poza obszarem baru, jego udział w masie wybrzuszenia zmienia się bardzo nieznacznie.

Problem z barem

Półdługość pręta, prędkość wzorca i orientacja pozostają niepewne. Niepewność ta przekłada się bezpośrednio na szacunki masy w odległości około 5 kpc.

6. Centralna czarna dziura – Sagittarius A*

Składnik 6 – Sagittarius A* – M = 4.0 × 10⁶ M⊙

W dynamicznym centrum Drogi Mlecznej znajduje się supermasywna czarna dziura Sagittarius A*. Jej masa jest mierzona z wysoką precyzją poprzez śledzenie orbit gwiazd w pobliżu Centrum Galaktyki.

\(\rho_{\mathrm{Sgr\,A^\ast}}(\mathbf{r})=M_{\mathrm{Sgr\,A^\ast}}\delta^{(3)}(\mathbf{r})\) \(M_{\mathrm{Sgr\,A^\ast}}(0\)

Pomimo swojej sławy, Sagittarius A* ma znikomy udział w globalnym budżecie masy. Jej znaczenie jest dynamiczne w najbardziej wewnętrznym parseku.

7. Gwiezdna aureola

Składnik 7 – halo gwiezdne – M ≈ 5 × 10⁸ do 10⁹ M⊙

Gwiezdne halo to rozproszona, z grubsza sferyczna populacja starych, ubogich w metale gwiazd otaczających dysk. Obejmuje ona gromady kuliste i strumienie gwiazd z rozbitych galaktyk karłowatych.

\(\rho_{\mathrm{halo,\star}}(r)=\rho_{0,\star}\left(\frac{r_0}{r}\right)^n,\qquad n\approx3\text{–}4\)

Dla n nie równego 3, masa skumulowana wynosi:

\(M_{\mathrm{halo,\star}}(<r)=\frac{4\pi\rho_{0,\star}r_0^n}{3-n}r^{3-n}\)

Dla n = 3:

\(M_{\mathrm{halo,\star}}(<r)=4\pi\rho_{0,\star}r_0^3\ln\left(\frac{r}{r_{\mathrm{min}}}\right)\)

Gwiezdne halo jest użyteczne jako znacznik kinematyczny, ale jego całkowita masa jest znacznie mniejsza niż niewidoczna masa wywnioskowana z krzywej rotacji.

8. Całkowita widoczna masa

Całkowita widoczna masa jest sumą dysku, gazu, wybrzuszenia, halo gwiezdnego i centralnej czarnej dziury:

\(M_{\mathrm{visible}}(<r)=M_{\mathrm{thin}}(<r)+M_{\mathrm{thick}}(<r)+M_{\mathrm{HI}}(<r)+M_{\mathrm{H_2}}(<r)+M_{\mathrm{bulge}}(<r)+M_{\mathrm{halo,\star}}(<r)+M_{\mathrm{Sgr\,A^\ast}}\)

Rozszerzona forma to:

\(M_{\mathrm{visible}}(<r)=3.52\times10^{10}\left[1-e^{-r/2.50}\left(1+\frac{r}{2.50}\right)\right]+1.05\times10^{10}\left[1-e^{-r/3.02}\left(1+\frac{r}{3.02}\right)\right]\) \(+1.1\times10^{10}f_{\mathrm{HI}}(r)+1.2\times10^9f_{\mathrm{H_2}}(r)+9.23\times10^9\left[1-e^{-r}\left(1+r+\frac{r^2}{2}\right)\right]+M_{\mathrm{halo,\star}}(<r)+4\times10^6\)
KomponentMasa całkowitaPromienie dominujące
Cienki dysk3.52 × 10¹⁰ M⊙0-15 kpc
Gruby dysk1.05 × 10¹⁰ M⊙0-15 kpc
Wybrzuszenie i pasek9.23 × 10⁹ M⊙0-4 kpc
Gaz HI1.1 × 10¹⁰ M⊙3-20 kpc
Gaz H₂1.2 × 10⁹ M⊙2-8 kpc
Gwiezdna aureola~10⁹ M⊙5-200 kpc
Strzelec A*4 × 10⁶ M⊙r = 0
Łącznie widoczne≈ 6.7 × 10¹⁰ M⊙

9. Brakująca masa – główny problem

Gdyby istniała tylko widoczna materia barionowa, prędkość rotacji spadałaby przy dużym promieniu:

\(V_{\mathrm{exp}}(r)=\sqrt{\frac{GM_{\mathrm{visible}}(<r)}{r}}\) \(r\gg R_d\quad\Longrightarrow\quad V_{\mathrm{exp}}(r)\propto\frac{1}{\sqrt{r}}\)

Zamiast tego obserwowana krzywa rotacji pozostaje w przybliżeniu płaska do dużego promienia i maleje tylko w zewnętrznych pomiarach z epoki Gaia. Masa dynamiczna wywnioskowana z kinematyki wynosi:

\(M_{\mathrm{dyn}}(<r)=\frac{rV_c^2(r)}{G}\) \(M_{\mathrm{dyn}}(<r)=2.325\times10^5\left(\frac{V_c(r)}{\mathrm{km/s}}\right)^2\left(\frac{r}{\mathrm{kpc}}\right)M_\odot\)

Niewidzialna masa to:

\(\boxed{M_{\mathrm{invisible}}(<r)=M_{\mathrm{dyn}}(<r)-M_{\mathrm{visible}}(<r)}\) \(\boxed{M_{\mathrm{invisible}}(<r)=\frac{rV_c^2(r)}{G}-M_{\mathrm{visible}}(<r)}\)

Na kole słonecznym, przy r = 8,2 kpc i Vc = 233 km/s:

\(M_{\mathrm{dyn}}(<8.2\,\mathrm{kpc})=2.325\times10^5\times233^2\times8.2\approx1.04\times10^{11}M_\odot\) \(M_{\mathrm{visible}}(<8.2\,\mathrm{kpc})\approx4.5\times10^{10}M_\odot\) \(M_{\mathrm{invisible}}(<8.2\,\mathrm{kpc})\approx5.5\times10^{10}M_\odot\)

Już przy promieniu Słońca niewidzialna masa jest porównywalna z masą widzialną. Przy większych promieniach niewidoczny składnik dominuje.

\(M_{\mathrm{Milky\ Way}}(<r)=M_{\mathrm{visible}}(<r)+M_{\mathrm{invisible}}(<r)\)

10. Promieniowe profile masy – symulacja

Poniższe wykresy obliczają przybliżone krzywe skumulowanej masy dla głównych widocznych składników, masy dynamicznej i wywnioskowanej niewidocznej masy. Porównują one również krzywą rotacji tylko barionowej ze schematyczną obserwowaną krzywą rotacji i punktami z epoki Gaia.

Skumulowana masa zamknięta M(<r) dla każdego składnika galaktycznego
Cienki dysk Gruby dysk Gaz HI Wybrzuszenie Całkowita widzialność Suma dynamiczna Masa niewidoczna
Obserwowana krzywa rotacji – widoczne składniki a pomiar dynamiczny
Tylko bariony Obserwowana Vc(r) Punkty z epoki Gaia