Galaktyki, krzywe rotacji i ciemna materia — ukryta masa Wszechświata

Struktura galaktyki: co możemy zobaczyć

Typowa galaktyka spiralna, jak Droga Mleczna, to grawitacyjnie związany układ zawierający 10&sup9;–10¹² gwiazd, gaz i pył międzygwiazdowy oraz (jak dziś sądzimy) dominujący komponent ciemnej materii. Morfologicznie widoczna materia jest zorganizowana w:

• Zgrubienie centralne (bulge): gęsty centralny sferoid starych, czerwonych gwiazd. Zgrubienie Drogi Mlecznej ma masę ∼2×10¹⁰ M_☉ i promień ∼3 kpc. Większość dużych spiral zawiera w centrum supermasywną czarną dziurę (Droga Mleczna: Sgr A*, 4,15×10⁶ M_☉).
• Dysk gwiazdowy: płaski, rotujący dysk gwiazd zgodny z eksponencjalnym profilem gęstości powierzchniowej Σ(R) = Σ₀ e^−R/R_d o promieniu skali R_d ∼3 kpc dla Drogi Mlecznej.
• Dysk gazowy: grubsza warstwa neutralnego wodoru (HI) i obłoków molekularnych, rozciągająca się poza dysk gwiazdowy. Emisja radiowa 21-cm mapuje rozkład HI i, co kluczowe, umożliwia pomiary krzywej rotacji daleko poza widoczny dysk optyczny.
• Halo gwiazdowe: rzadki sferyczny rozkład starych gromad kulistych i gwiazd halo, rozciągający się do ∼100 kpc.

Całkowita widoczna masa galaktyki takiej jak Droga Mleczna to ∼5×10¹⁰ M_☉ w gwiazdach plus ∼10¹⁰ M_☉ w gazie. Ale masa dynamiczna wywnioskowana z kinematyki wynosi ∼10¹² M_☉ — mniej więcej 20 razy więcej. Ta rozbieżność to problem ciemnej materii.

Krzywe rotacji — dymiący pistolet

Dla cząstki testowej na orbicie kołowej w promieniu R od rozkładu masy M(<R), grawitacja newtonowska wymaga:

W naszym Układzie Słonecznym planety idealnie podążają za tą keplerowską zależnością v ∝ R^−1/2: Merkury krąży z prędkością 47,9 km/s, Neptun z 5,4 km/s. Gdyby masa w galaktyce spiralnej była skoncentrowana w centrum jak w Układzie Słonecznym, gwiazdy w zewnętrznym dysku powinny wykazywać taki sam malejący przebieg.

Vera Rubin i Kent Ford rozpoczęli systematyczne pomiary rotacji galaktyk pod koniec lat 60. XX wieku, wykorzystując przesunięcie Dopplera linii emisyjnych Hα. Do 1980 roku ich artykuł dotyczący 21 spiral ujawnił niewątpliwy wzorzec: krzywe rotacji są płaskie przy dużych promieniach, z v(R) ≈ stałą aż do najdalszego mierzalnego punktu. Późniejsze przeglądy radiowe HI 21-cm (sięgające kilkukrotności promienia optycznego) potwierdziły płaskie lub rosnące krzywe w setkach galaktyk.

Płaska krzywa v(R) = v_flat implikuje M(<R) ∝ R — masa nadal rośnie liniowo z promieniem nawet tam, gdzie praktycznie nie ma widocznych gwiazd. Wymagana niewidoczna masa zwykle przewyższa widoczną masę czynnikiem 5–30.

Profil halo ciemnej materii NFW

Symulacje N-ciałowe formowania struktur we Wszechświecie z zimną ciemną materią (CDM) konsekwentnie wytwarzają halo ciemnej materii o charakterystycznym profilu gęstości. Julio Navarro, Carlos Frenk i Simon White (1995, 1996) odkryli, że halo obejmujące 10 rzędów wielkości masy podążają za tą samą postacią funkcyjną — profilem NFW:

Profil NFW ma dwa kontrastujące zachowania: przy małych promieniach (r « r_s) skaluje się jako ρ ∝ r²/r−1 = r⁻¹ (grot wewnętrznego rdzenia); przy dużych promieniach (r » r_s) skaluje się jako ρ ∝ r⁻³. Całkowanie daje masę zamkniętą:

Profil NFW sam w sobie nie wytwarza dokładnie płaskiej krzywej, ale w połączeniu z dyskami gwiazdowym i gazowym całkowite v(r) jest niezwykle płaskie w obserwowalnym zakresie — zgodnie z obserwacjami dobrze zbadanych spiral.

Kontrowersja grot-rdzeń

Wewnętrzny grot r⁻¹ NFW jest przewidywany przez symulacje CDM, ale obserwacje galaktyk karłowatych często sugerują płaski centralny rdzeń gęstości (ρ = const). Ten problem grot-rdzeń pozostaje aktywnym obszarem badań. Proponowane rozwiązania obejmują: sprzężenie zwrotne barionowe (eksplozje supernowych wydmuchują gaz, redukując koncentrację masy centralnej i „podgrzewając" ciemną materię), ciepłą ciemną materię (WDM) z odcięciem swobodnego strumienia, które wygładza strukturę w małej skali, lub samooddziałującą ciemną materię (SIDM), która tworzy izotermiczne rdzenie poprzez zderzenia.

Komponenty dysku Freemana i zgrubienia

Widoczny dysk gwiazdowy wnosi własny potencjał grawitacyjny. Dla cienkiego dysku eksponencjalnego (Freeman 1970) wkład prędkości kołowej w promieniu R w płaszczyźnie dysku oblicza się z całek eliptycznych, ale doskonałe przybliżenie analityczne (van der Kruit i Searle) wykorzystuje zmodyfikowane funkcje Bessela:

Zgrubienie przybliżone jako sfera Hernquista (gęstość ρ ∝ 1/[r(r+a)³]) lub sferoid de Vaucouleurs R^1/4 dodaje stromy rosnący wkład do v(r) w pobliżu centrum. Dobra dekompozycja krzywej rotacji spirali dzieli więc: zgrubienie (wewnętrzny wzrost), dysk (pośredni garb) i halo DM (zewnętrzna płaska część). Ta metoda dopasowania stosunku masy do jasności jest standardową metodą wnioskowania o frakcjach ciemnej materii w poszczególnych galaktykach.

Dowody na ciemną materię poza krzywymi rotacji

Soczewkowanie grawitacyjne

Ogólna teoria względności przewiduje, że masa zakrzywia światło. Odległe gromady galaktyk działają jak soczewki grawitacyjne, zniekształcając i powiększając galaktyki tła. Siła soczewkowania bada całkowitą masę, niezależnie od tego, czy ta masa emituje światło. Porównanie masy z soczewkowania z masą gazu rentgenowskiego (dominującym widocznym komponentem barionowym w gromadach) konsekwentnie pokazuje czynnik ∼5–7 więcej masy w ciemnej materii. Kąt ugięcia promienia światła przechodzącego z parametrem zderzenia b w pobliżu masy M wynosi:

Gromada Kuli

Najbardziej bezpośredni dowód na ciemną materię jako odrębny komponent pochodzi z Gromady Kuli (1E 0657-558), dwóch gromad galaktyk, które zderzyły się ∼100 mln lat temu. Obserwacje rentgenowskie Chandra pokazują gorący gaz (główny komponent barionowy) spowolniony i wstrząśnięty w punkcie kolizji. Ale mapy słabego soczewkowania grawitacyjnego — które śledzą całkowitą masę — pokazują masę przesuniętą przed gazem, pokrywającą się z galaktykami i gwiazdami, które przeleciały bez interakcji. Ciemna materia, podobnie jak gwiazdy, przeleciała na wprost, ponieważ ciemna materia oddziałuje tylko grawitacyjnie (i słabo, jeśli w ogóle, poprzez inne siły). To przestrzenne rozdzielenie masy ciemnej i barionowej jest zasadniczo niemożliwe do wyjaśnienia zmodyfikowaną grawitacją.

Kosmiczne mikrofalowe promieniowanie tła

Widmo mocy CMB (piki akustyczne w fluktuacjach temperatury przy z ≈ 1100) zależy wrażliwie od gęstości barionów Ω_b w stosunku do całkowitej materii Ω_m. Wyniki Planck 2018 dają:

Wszechświat bez ciemnej materii nie może odtworzyć obserwowanego drugiego i trzeciego piku akustycznego. Dane CMB niezależnie wymagają niebarionowego komponentu materii mniej więcej 5 razy liczniejszego niż zwykłe atomy.

Pierwotna nukleosynteza (Big Bang nucleosynthesis)

BBN ogranicza gęstość barionów na podstawie obserwowanych pierwotnych obfitości D, ³He, ⁴He i ⁷Li. Daje Ω_bh² ≈ 0,022 — w doskonałej zgodności z pomiarami CMB, ale znacznie poniżej całkowitej gęstości materii wymaganej przez formowanie struktur i krzywe rotacji. Ta niezależna weryfikacja z innej epoki (t ≈ 3 min) mocno ustanawia, że ciemna materia nie może być barionowa.

Czym jest ciemna materia?

Pomimo przytłaczających dowodów na grawitacyjne skutki ciemnej materii, jej natura cząstkowa pozostaje nieznana. Główne klasy kandydatów to:

• WIMP-y (Weakly Interacting Massive Particles, m ∼ 10–1000 GeV/c²): pierwszy wybór przez dekady, ponieważ supersymetryczne rozszerzenia Modelu Standardowego naturalnie przewidują stabilną najlżejszą cząstkę supersymetryczną (LSP) o w przybliżeniu właściwej gęstości reliktowej („cud WIMP-ów"). Eksperymenty detekcji bezpośredniej (LUX, PandaX, XENONnT) wykluczyły duże obszary przestrzeni parametrów bez pozytywnego sygnału.
• Aksjony (m ∼ 1–1000 μeV/c²): ultralekkie cząstki skalarne pierwotnie zaproponowane do rozwiązania problemu silnego CP w QCD. Ciemna materia aksjonowa tworzy kondensaty Bosego-Einsteina (rozmyta ciemna materia), które mogą rozwiązać problem grot-rdzeń. Eksperymenty: ADMX, CASPEr.
• Neutrina sterylne: hipotetyczne ciężkie neutrina (keV–GeV), które słabo mieszają się z neutrinami aktywnymi. Rozpadałyby się na promieniowanie rentgenowskie: wstępna linia rentgenowska 3,5 keV z gromad galaktyk przyciągnęła wiele uwagi (2014), ale pozostaje niepotwierdzona.
• PBH (pierwotne czarne dziury): czarne dziury uformowane z fluktuacji gęstości we wczesnym Wszechświecie. Przeglądy mikrosoczewkowania (MACHO, EROS, HSC) ograniczają PBH do wąskiego okna mas rzędu asteroid. Detekcje fal grawitacyjnych LIGO/Virgo odnowiły zainteresowanie ciemną materią PBH, ale obecne dane nie wspierają znaczącego wkładu.

Zmodyfikowana grawitacja: MOND

Mordehai Milgrom zaproponował w 1983 roku, że ciemna materia jest niepotrzebna, jeśli drugie prawo Newtona zostanie zmodyfikowane przy bardzo niskich przyspieszeniach. Jego Zmodyfikowana Dynamika Newtonowska (MOND):

MOND naturalnie wyjaśnia barionową zależność Tully'ego-Fishera (ścisłą empiryczną korelację v_flat⁴ ∝ M_bary), która utrzymuje się w pięciu rzędach wielkości masy galaktyk. To uderzający sukces: w paradygmacie ciemnej materii korelację trzeba wyjaśnić precyzyjnie dostrojoną zależnością między masą barionową a masą halo. MOND jednak poważnie zawodzi w gromadach galaktyk (gdzie nadal wymaga ciemnej materii — lub gorących neutrin sterylnych) i w przypadku Gromady Kuli. Pełne relatywistyczne rozszerzenie (TeVeS, Bekenstein 2004, i jego następcy) jest potrzebne dla soczewkowania, ale te teorie napotykają poważne wyzwania teoretyczne i napięcia z obserwacjami fal grawitacyjnych.

Struktura kosmologiczna: ΛCDM

Model ΛCDM (Λ zimna ciemna materia) łączy ciemną materię ze stałą kosmologiczną Λ (ciemna energia) i opisuje formowanie się wielkoskalowej struktury z małych początkowych fluktuacji kwantowych wzmocnionych przez inflację. Kluczowe parametry (Planck 2018):

W ΛCDM ciemna materia jest zimna (nierelatywistyczna w momencie rozprzęgnięcia) i bezzderzeniowa. Zapewnia grawitacyjny rusztowanie, w które opadają baryony, tworząc galaktyki. Kosmiczna sieć — filamenty, pustki i węzły gromad — to skamieniały zapis formowania struktury ciemnej materii. Symulacje (Millennium, IllustrisTNG, EAGLE) odtwarzają obserwowany rozkład galaktyk, funkcję masy gromad i funkcję korelacji dwupunktowej z niezwykłą dokładnością.

Dwa napięcia w ΛCDM

Pomimo sukcesu, ΛCDM napotyka dwa istotne napięcia w obecnych danych. Napięcie H₀: lokalne pomiary drabiny odległości (cefeidy + supernowe typu Ia) dają H₀ ≈ 73 km/s/Mpc, niezgodne z wartością wywnioskowaną z CMB wynoszącą 67 km/s/Mpc na poziomie 4–5σ. Napięcie σ₈: przeglądy słabego soczewkowania (KiDS, DES) preferują nieco mniejszą amplitudę grupowania niż przewiduje CMB. Czy te napięcia wskazują na nową fizykę poza ΛCDM, czy na błędy systematyczne, pozostaje aktywnie dyskutowane.

Gwiazdy podwójne i dynamika galaktyczna

Wewnątrz galaktyk dynamika gwiazdowa również ujawnia rozkład masy poprzez metody inne niż krzywe rotacji. Twierdzenie o wiriale (2⟨T⟩ = −⟨V⟩ dla grawitacyjnie związanego układu w równowadze) wiąże dyspersję prędkości galaktyki eliptycznej lub gromady gwiazd z jej całkowitą masą: M ≈ 5σ²R/G. Orbity gwiazd podwójnych mierzą masy gwiazdowe bezpośrednio poprzez trzecie prawo Keplera, a poprzez krzyżkorelację widm obu składowych parametry orbitalne można wyznaczyć precyzyjnie.

W największych skalach osobliwe prędkości galaktyk w gromadach, profile temperatury rentgenowskiej gazu wewnątrzgromadowego oraz efekt Siunajewa-Zeldowicza (odwrotne rozpraszanie Comptona fotonów CMB na gorących elektronach gromady) niezależnie mierzą całkowitą masę — i wszystkie dają tę samą odpowiedź zdominowaną przez ciemną materię.

Symulacje na tej platformie

Odkrywaj fizykę galaktyczną i ciemnej materii za pomocą tych interaktywnych narzędzi:

• Krzywa rotacji galaktyki — galaktyka spiralna widziana z góry z halo NFW, dyskiem Freemana, zgrubieniem; dekompozycja v(r)
• Soczewkowanie grawitacyjne — silne soczewkowanie (pierścienie Einsteina, łuki), krzywe blasku mikrosoczewkowania, potencjał Fermata
• Gwiazdy podwójne — orbity keplerowskie dwóch ciał, stosunek mas, mimośród, krzywe prędkości radialnej
• Ciemna materia — interaktywne modele halo, profile gęstości, dekompozycja krzywej rotacji
• Pierwotna nukleosynteza — pierwotne obfitości pierwiastków w funkcji gęstości barionów
• Kosmiczne mikrofalowe promieniowanie tła — widmo mocy CMB, piki akustyczne, parametry kosmologiczne