Krzywe Rotacji Galaktyk: Przypadek dla Ciemnej Materii

Krzywa, która nie chce spaść

Skieruj teleskop na galaktykę spiralną, zmierz przesunięcie ku czerwieni gazu wodnego lub gwiazd w rosnącej odległości od centrum i narysuj prędkość obrotową względem promienia. Newtonowska grawitacja przewiduje to stanowczo: po przekroczeniu większości widocznej masy, prędkość powinna spadać zgodnie z wzorem v(r) ∝ 1/√r, tym samym spadkiem Keplera, który rządzi planetami wokół Słońca. To jednak nie obserwujemy w galaktykach. Od około brzegu widocznego dysku na zewnątrz, v(r) utrzymuje się na płaszczyźnie – często w granicach 10% od swojej wartości szczytowej – aż po ostatnie zmierzone chmury gazu, rozciągające się na dziesiątki kiloparsecs.

Vera Rubin i Kent Ford ostatecznie to potwierdzili w latach 70. i 80. XX wieku, poszerzając wcześniejsze wskazówki z badania Fritz Zwicky'ego z 1933 roku dotyczącego gromady Coma. Rubin zmierzył krzywe obrotu dla kilkudziesięciu galaktyk spiralnych za pomocą spektroskopii optycznej i zawsze odnajdywał tę samą płaską formę, niezależnie od rozmiaru lub jasności galaktyki. Widzialne gwiazdy i gaz – to, co emituje światło – stanowią jedynie ułamek grawitacji potrzebnej do utrzymania tych zewnętrznych orbit w danej prędkości.

Obliczanie wartości

Przemieszczenie obrotowe testowej masy na orbicie okręgłej jest ustalone przez masę znajdującą się wewnątrz jej promienia: v(r) = √(G M(r) / r). Model galaktyki jako trzy komponenty – centralny bul, cienkie dyskowy płaszcz gwiezdny o profilu powierzchniowej jasności wykładniczej i gaz – oraz każdy z nich wnosi własne v(r), dodane w kwadratury. Zsumuj je i krzywa kompozytowa osiąga szczyt blisko krawędzi dysku optycznego, a następnie spada, dokładnie jak planety. Obserwowana krzywa nie opada. Do przyczyniania się do masy na dużym promieniu dochodzi coś innego, co nie emituje światła.

Halo NFW

Standardne rozwiązanie polega na tworzeniu przybliżonego sferycznego halo ciemnej materii otaczającego widoczny dysk, z znacznie większą masą niż gwiazdy i gaz razem wzięte. Kosmologiczne symulacje N-ciał zimnej ciemnej materii konsekwentnie generują halo o charakterystycznym profilu gęstości, nazwanym na cześć Navarro'ego, Frenka i White’a (1996):

ρ(r) = ρ0 / [ (r/rs) · (1 + r/rs)² ] półokość ρ0 = charakterystyczna gęstość, rs = półos wie, obszar wewnętrzny (r << rs): ρ ∝ 1/r — ostry, obszar zewnętrzny (r >> rs): ρ ∝ 1/r³ — stromy spadek. Integracja tego profilu w celu uzyskania masy zamkniętej M(r) i wprowadzenie jej do v(r) = √(GM(r)/r) generuje wkład do krzywej obrotowej, który delikatnie rośnie i następnie płaszczy się na bardzo szerokim zakresie promieni — dokładnie to brakujące element. Dodanie wkładu halo do masy dysku i bulgu daje spójną krzywą: przebiega przez centrum zdominowane przez bulg, ma powinowactwo, gdzie osiąga szczyt dysk i pozostaje płaska, ponieważ rosnąca masa zamknięta halo kompensuje rosnący promień.

ρ(r) = ρ0 / [ (r/rs) · (1 + r/rs)² ]

ρ0 = characteristic density,  rs = scale radius
inner region  (r << rs):  ρ ∝ 1/r        — cuspy
outer region  (r >> rs):  ρ ∝ 1/r³       — steep falloff

Dlaczego nie po prostu naprawić grawitacji?

Główna alternatywna teoria, Dynamika Nowo-Newtonowska (MOND), zastępuje halo zmianą w prawie siły poniżej niewielkiego progu przyspieszenia a0 ≈ 1.2×10⁻¹⁰ m/s², i imponująco pasuje do indywidualnych krzywych obrotu z jednym swobodnym parametrem na galaktykę. Problemy pojawiają się jednak wszędzie tam, gdzie MOND nie daje odpowiedzi: rozdzielone masy i światło w Grupie Kulowej (Bullet Cluster), szczyty akustyczne w mikrofalowym promieniowaniu kosmicznym oraz formowanie się dużych struktur. Obraz halo ciemnej materii przetrwał wszystkie te testy jednocześnie, dlatego pozostaje modelem roboczym, pomimo braku bezpośredniej detekcji cząstki.

Co robi symulacja

Symulacja buduje v(r) z dysku, profil wypukłości i profil masy NFW, które możesz niezależnie dostosowywać. Następnie nakłada się na to czystą predykcję Newtona, gwiazd tylko, aby w ten sposób zobaczyć przerwę. Ustaw masę halo na zero i otrzymasz upadek Keplera, którego Rubin nigdy nie widział; ustaw realistyczny udział masy halo – ciemna materia przeważa nad materią barionową około 5 do 1 w dzisiejszym budżecie kosmologicznym – a krzywa staje się płaska, odpowiadając rzeczywistym pomiarom 21 cm i krzywej obrotu optycznego.

Frequently asked questions

Wypróbuj na żywo

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