BeeTheory – Fizyka grawitacji i fal
Równania radialne ukrytej masy Drogi Mlecznej
Od profili gęstości do całek pierścieniowych i krzywych rotacji – matematyczne ujęcie ukrytej masy w funkcji promienia galaktyki R.
Ta strona przedstawia równania radialne używane do opisu ukrytej masy Drogi Mlecznej. Porównuje ona klasyczne profile gęstości ciemnej materii, całki pierścieniowe i powłokowe, równania masy zamkniętej, krzywe rotacji oraz interpretację brakującej masy jako możliwego efektu interferencji falowej.
~10¹² M⊙
Klasyczne oszacowanie całkowitej masy halo Drogi Mlecznej.
0,39 GeV/cm³
Typowa lokalna gęstość ciemnej materii w pobliżu Słońca.
R⊙ ≈ 8 kpc
Przybliżona odległość Słońca od Centrum Galaktyki.
~200 kpc
Przybliżona skala zewnętrzna użyta do oszacowania halo Drogi Mlecznej.
Spis treści
- Dlaczego brakuje masy?
- Profile gęstości ρ(R)
- Masa pierścienia i pierścienia dM/dR
- Zamknięta masa ciemnej materii M(<R)
- Krzywa rotacji V(R)
- Aktualne szacunki obserwacyjne
- Konkurujące hipotezy
- Perspektywa BeeTheory
1. Dlaczego brakuje masy?
W 1933 roku Fritz Zwicky zauważył, że galaktyki w gromadzie Coma poruszały się zbyt szybko, by mogły być utrzymywane razem wyłącznie dzięki swojej widocznej masie. W latach siedemdziesiątych Vera Rubin i Kent Ford zmierzyli krzywe rotacji galaktyk spiralnych i odkryli coś równie uderzającego: gwiazdy na dużych promieniach orbitują prawie tak szybko, jak te w pobliżu centrum, podczas gdy grawitacja Newtona z widocznej materii przewiduje, że powinny one zwolnić.
Dla prostej orbity Keplera wokół centralnej masy, oczekujemy:
\(V(R)\propto \frac{1}{\sqrt{R}}\)Zamiast tego obserwuje się w przybliżeniu płaską lub tylko powoli opadającą krzywą rotacji:
\(V(R)\approx V_{\infty}=\mathrm{const}\qquad \mathrm{for}\ R\gtrsim 5\,\mathrm{kpc}\)Pogodzenie tych faktów z grawitacją newtonowską wymaga dodatkowego niewidzialnego składnika masy, którego gęstość spada w przybliżeniu tak:
\(\rho(r)\propto r^{-2}\)Daje to całkowitą masę zamkniętą proporcjonalną do promienia:
\(M(<R)\propto R\)i dlatego:
\(V\propto \sqrt{\frac{M}{R}}=\mathrm{const}\)Kluczowa zagadka ilościowa
Świecąca masa barionowa Drogi Mlecznej wynosi około 5 × 10¹⁰ M⊙. Całkowita masa dynamiczna wywnioskowana z kinematyki do około 200 kpc wynosi około 10¹² M⊙. Oznacza to, że stosunek masy ciemnej do jasnej wynosi około 10 do 1.
2. Profile gęstości ρ(R)
Profil gęstości to funkcja matematyczna opisująca, w jaki sposób gęstość ciemnej materii ρ zmienia się wraz z promieniem galaktocentrycznym r lub promieniem cylindrycznym R w płaszczyźnie galaktyki.
2.1 Profil NFW
Profil NFW, wprowadzony przez Navarro, Frenka i White’a, pochodzi z kosmologicznych symulacji N-ciałowych. Tworzy on charakterystyczne podwójne prawo potęgowe z centralnym wierzchołkiem.
\(\rho_{\mathrm{NFW}}(r)=\frac{\rho_0}{\left(\frac{r}{r_s}\right)\left(1+\frac{r}{r_s}\right)^2}\)| Parametr | Symbol | Oszacowanie Drogi Mlecznej | Rola |
|---|---|---|---|
| Promień skali | rs | 15-25 kpc | Przejście między zboczem wewnętrznym i zewnętrznym |
| Gęstość charakterystyczna | ρ0 | Skalibrowany do lokalnej gęstości ciemnej materii | Ogólna normalizacja |
| Nachylenie wewnętrzne | γ | -1 | Złe zachowanie |
| Zbocze zewnętrzne | – | -3 | Gwałtowny spadek przy dużym promieniu |
2.2 Profil Einasto
Profil Einasto unika ścisłej centralnej dywergencji i wykorzystuje parametr kształtu α, który pozwala na płynną zmianę nachylenia gęstości wraz z promieniem.
\(\rho_{\mathrm{Ein}}(r)=\rho_{-2}\exp\left\{-\frac{2}{\alpha}\left[\left(\frac{r}{r_{-2}}\right)^\alpha-1\right]\right\}\)| Parametr | Symbol | Oszacowanie Drogi Mlecznej | Rola |
|---|---|---|---|
| Wskaźnik kształtu | α | Zależne od modelu | Kontroluje szybkość zmiany nachylenia |
| Promień skali | r-2 | ~18-22 kpc | Promień, przy którym nachylenie logarytmu jest równe -2 |
| Gęstość na r-2 | ρ-2 | Skalibrowany do lokalnej gęstości | Normalizacja |
Ostatnie napięcie obserwacyjne
Ostatnie badania oparte na Gaia sugerują, że krzywa rotacji Drogi Mlecznej może zmniejszać się szybciej poza promieniem Słońca niż przewidywałoby to standardowe halo NFW. To sprawia, że profile rdzeniowe lub płynnie zmieniające się profile, takie jak Einasto, są szczególnie ważne w obecnych dyskusjach.
2.3 Profil pseudo-izotermiczny
Profil pseudo-izotermiczny jest często używany jako proste analityczne przybliżenie dla rdzeniowego halo.
\(\rho_{\mathrm{iso}}(r)=\frac{\rho_0}{1+\left(\frac{r}{r_s}\right)^2}\)Przy małym promieniu gęstość zbliża się do stałej wartości. Przy dużym promieniu spada jako r-² i tworzy płaską krzywą rotacji.
\(V_{\infty}=\sqrt{4\pi G\rho_0 r_s^2}\)Problem cusp versus core
Symulacje N-ciałowe często przewidują wypukłe profile NFW, podczas gdy wiele obserwowanych galaktyk karłowatych wydaje się preferować rdzeniowe profile gęstości. Problem ten pozostaje jedną z głównych nierozwiązanych kwestii w fizyce ciemnej materii.
3. Masa pierścienia i pierścienia osadczego – dM/dR
Aby obliczyć, ile ciemnej materii znajduje się w każdym radialnym wycinku Galaktyki, całkujemy gęstość po cienkiej powłoce lub pierścieniu. Geometria zależy od tego, czy halo jest traktowane jako sferyczne czy spłaszczone.
3.1 Sferyczna cienka powłoka
Dla sferycznie symetrycznej aureoli masa w powłoce o grubości dr i promieniu r wynosi:
\(\frac{dM_{\mathrm{shell}}}{dr}=4\pi r^2\rho(r)\)3.2 Pierścień pierścieniowy z płaszczyzną dysku
Dla pierścienia leżącego w płaszczyźnie Galaktyki, o promieniu cylindrycznym R i efektywnej grubości połówkowej H(R), masa pierścienia wynosi:
\(dM_{\mathrm{ann}}=2\pi R\,\rho(R,0)\,2H(R)\,dR\)W przypadku sferycznej aureoli można to zapisać jako całkę na wysokości z:
\(\frac{dM}{dR}=2\pi R\int_{-\infty}^{+\infty}\rho\left(\sqrt{R^2+z^2}\right)dz\)W przybliżeniu sferycznym, łączy się to z powrotem do:
\(\frac{dM}{dR}\approx4\pi R^2\rho(R)\)3.3 Masa NFW na łuskę
\(\frac{dM_{\mathrm{NFW}}}{dr}=4\pi r^2\frac{\rho_0}{\left(\frac{r}{r_s}\right)\left(1+\frac{r}{r_s}\right)^2}=\frac{4\pi\rho_0 r_s r}{\left(1+\frac{r}{r_s}\right)^2}\)Funkcja ta osiąga szczyt wokół promienia skali rs, co oznacza, że znaczna część masy ciemnej materii na powłokę jest zdeponowana w pośrednim halo, a nie tylko w centrum lub na obrzeżach.
3.4 Masa Einasto na powłokę
\(\frac{dM_{\mathrm{Ein}}}{dr}=4\pi r^2\rho_{-2}\exp\left\{-\frac{2}{\alpha}\left[\left(\frac{r}{r_{-2}}\right)^\alpha-1\right]\right\}\)Masa zamknięta Einasto jest zwykle obliczana numerycznie.
Znaczenie fizyczne
Funkcja dM/dr mówi nam, który promień Galaktyki najbardziej przyczynia się do ukrytego budżetu masy. Bardziej stromy profil zewnętrzny zmniejsza wnioskowaną całkowitą masę halo, podczas gdy płytszy profil ją zwiększa.
4. Zamknięta masa ciemnej materii M(
Całkowanie elementu powłoki od 0 do R daje całkowitą masę ciemnej materii zamkniętą w promieniu R:
\(M_{\mathrm{DM}}(<R)=\int_0^R4\pi r^2\rho(r)\,dr\)
4.1 Masa zamknięta NFW
\(M_{\mathrm{NFW}}(<R)=4\pi\rho_0r_s^3\left[\ln\left(1+\frac{R}{r_s}\right)-\frac{R/r_s}{1+R/r_s}\right]\)
4.2 Masa zamknięta Einasto
\(M_{\mathrm{Ein}}(<R)=4\pi\rho_{-2}\int_0^R r^2\exp\left\{-\frac{2}{\alpha}\left[\left(\frac{r}{r_{-2}}\right)^\alpha-1\right]\right\}dr\)
4.3 Całkowity rozkład masy
Całkowitą zamkniętą masę dynamiczną można rozłożyć na składniki widoczne i ukryte:
\(M_{\mathrm{tot}}(<R)=M_{\mathrm{bulge}}(<R)+M_{\mathrm{disk}}(<R)+M_{\mathrm{DM}}(<R)\)
~3 × 10¹⁰ M⊙
Przybliżona masa ukryta wewnątrz okręgu słonecznego.
~1-2 × 10¹¹ M⊙
Przybliżona ukryta masa wewnątrz 20 kpc.
5-12 × 10¹¹ M⊙
Przybliżona masa ukryta w obszarze wirialnym, silnie zależna od modelu.
Profil masy pozostaje zależny od modelu.
Oszacowania masy halo Drogi Mlecznej silnie zależą od tego, w jaki sposób zewnętrzne halo jest ekstrapolowane poza region z silnymi ograniczeniami obserwacyjnymi.
Całkowanie elementu powłoki od 0 do R daje całkowitą masę ciemnej materii zamkniętą w promieniu R:
\(M_{\mathrm{DM}}(<R)=\int_0^R4\pi r^2\rho(r)\,dr\)4.1 Masa zamknięta NFW
\(M_{\mathrm{NFW}}(<R)=4\pi\rho_0r_s^3\left[\ln\left(1+\frac{R}{r_s}\right)-\frac{R/r_s}{1+R/r_s}\right]\)4.2 Masa zamknięta Einasto
\(M_{\mathrm{Ein}}(<R)=4\pi\rho_{-2}\int_0^R r^2\exp\left\{-\frac{2}{\alpha}\left[\left(\frac{r}{r_{-2}}\right)^\alpha-1\right]\right\}dr\)4.3 Całkowity rozkład masy
Całkowitą zamkniętą masę dynamiczną można rozłożyć na składniki widoczne i ukryte:
\(M_{\mathrm{tot}}(<R)=M_{\mathrm{bulge}}(<R)+M_{\mathrm{disk}}(<R)+M_{\mathrm{DM}}(<R)\)~3 × 10¹⁰ M⊙
Przybliżona masa ukryta wewnątrz okręgu słonecznego.
~1-2 × 10¹¹ M⊙
Przybliżona ukryta masa wewnątrz 20 kpc.
5-12 × 10¹¹ M⊙
Przybliżona masa ukryta w obszarze wirialnym, silnie zależna od modelu.
Profil masy pozostaje zależny od modelu.
Oszacowania masy halo Drogi Mlecznej silnie zależą od tego, w jaki sposób zewnętrzne halo jest ekstrapolowane poza region z silnymi ograniczeniami obserwacyjnymi.
5. Krzywa rotacji V(R)
Prędkość kołowa w promieniu R jest ustalana przez całkowitą zamkniętą masę poprzez równowagę przyciągania grawitacyjnego i przyspieszenia dośrodkowego:
\(V_c(R)=\sqrt{\frac{G\,M_{\mathrm{tot}}(<R)}{R}}\)Ponieważ niezależne składowe masy przyczyniają się do potencjału grawitacyjnego, ich składowe prędkości są często dodawane w kwadraturze:
\(V_c^2(R)=V_{\mathrm{bulge}}^2(R)+V_{\mathrm{thin\,disk}}^2(R)+V_{\mathrm{thick\,disk}}^2(R)+V_{\mathrm{gas}}^2(R)+V_{\mathrm{DM}}^2(R)\)5.1 Wkład dysku barionowego
Gwiezdny cienki dysk ma wykładniczy profil gęstości powierzchniowej:
\(\Sigma(R)=\Sigma_0\exp\left(-\frac{R}{R_d}\right)\)Odpowiednia prędkość kołowa dla dysku wykładniczego wynosi:
\(V_{\mathrm{disk}}^2(R)=\frac{2GM_d}{R_d}y^2\left[I_0(y)K_0(y)-I_1(y)K_1(y)\right],\qquad y=\frac{R}{2R_d}\)In iKn są zmodyfikowanymi funkcjami Bessela. Typowe parametry cienkiego dysku Drogi Mlecznej to Rd ≈ 2,6 kpc i Md ≈ 3,5 × 10¹⁰ M⊙.
5.2 Wkład ciemnej materii
\(V_{\mathrm{DM,NFW}}(R)=\sqrt{\frac{4\pi G\rho_0r_s^3}{R}\left[\ln\left(1+\frac{R}{r_s}\right)-\frac{R/r_s}{1+R/r_s}\right]}\)5.3 Baryoniczna relacja Tully’ego-Fishera
Barionowa zależność Tully’ego-Fishera łączy płaską prędkość rotacji galaktyki z jej całkowitą masą barionową:
\(V_{\infty}^4=G\,M_{\mathrm{bar}}\,a_0,\qquad a_0\approx1.2\times10^{-10}\,\mathrm{m/s^2}\)~230 km/s
Prędkość kołowa w pobliżu promienia Słońca.
~170-180 km/s
Możliwa malejąca wartość na dysku zewnętrznym, w zależności od danych znacznika.
~150 km/s
Przybliżona skala prędkości zewnętrznego halo na podstawie znaczników halo.
6. Aktualne szacunki obserwacyjne
Poniższa tabela podsumowuje reprezentatywne wartości gęstości i masy ciemnej materii w kluczowych promieniach Galaktyki. Dokładne wartości różnią się w zależności od zestawu danych, populacji znaczników i modelu halo.
| Promień R | Gęstość ciemnej materii | Zamknięta ciemna masa | Metoda |
|---|---|---|---|
| Centrum | Rozbieżne w NFW, skończone w modelach rdzeniowych | Zależne od modelu | Symulacje ciał N i modelowanie wewnętrznej galaktyki |
| R⊙ ≈ 8 kpc | ~0,39 GeV/cm³ | ~3 × 10¹⁰ M⊙ | Krzywa obrotu i kinematyka pionowa |
| 20 kpc | ~0,05 GeV/cm³ | ~1-2 × 10¹¹ M⊙ | Gaia i znaczniki spektroskopowe |
| 50 kpc | ~5 × 10-³ GeV/cm³ | ~3-5 × 10¹¹ M⊙ | Gromady kuliste i gwiazdy halo |
| 100-200 kpc | ≤10-³ GeV/cm³ | ~5-12 × 10¹¹ M⊙ | Galaktyki satelitarne i metody prędkości ucieczki |
Połączenie kinematyki gromad kulistych, gwiazd halo, galaktyk satelitarnych i astrometrii Gaia sugeruje, że profil zewnętrznego halo Drogi Mlecznej pozostaje niepewny. Niepewność ta ma kluczowe znaczenie dla problemu ukrytej masy.
7. Konkurujące hipotezy dotyczące brakującej masy
Kilka głównych rodzin wyjaśnień pozostaje aktywnych. Żadna z nich nie została ostatecznie potwierdzona lub wykluczona we wszystkich skalach obserwacyjnych.
7.1 Zimne cząstki ciemnej materii
Wiodącym paradygmatem pozostaje zimna ciemna materia. Cząstki kandydujące obejmują WIMP-y, neutrina sterylne i inne możliwości wykraczające poza model standardowy. Kandydaci ci tworzą rozszerzone halo, często modelowane za pomocą profili NFW lub Einasto.
\(m_{\chi}\sim10\text{–}1000\,\mathrm{GeV}\)Główne napięcie ma charakter eksperymentalny: bezpośrednia detekcja nie znalazła jeszcze potwierdzonej cząstki ciemnej materii.
7.2 Ultralekka lub rozmyta ciemna materia
Rozmyta ciemna materia wykorzystuje ultralekkie cząstki, których długość fali de Broglie’a może stać się astrofizycznie duża, tłumiąc strukturę w małej skali.
\(m_a\sim10^{-22}\,\mathrm{eV}\) \(\lambda_{\mathrm{dB}}\sim\mathrm{kpc}\)Ramy te naturalnie wytwarzają gładsze wewnętrzne rdzenie gęstości, ale są ograniczone przez dane z lasu Lyman-alfa i strukturę galaktyk karłowatych.
7.3 Zmodyfikowana dynamika newtonowska
MOND modyfikuje efektywne przyspieszenie grawitacyjne poniżej charakterystycznej skali:
\(a_0\approx1.2\times10^{-10}\,\mathrm{m/s^2}\)W reżimie głębokiego MOND efektywne przyspieszenie wynosi:
\(g_{\mathrm{eff}}=\sqrt{g_Na_0}\)MOND przewiduje barionową relację Tully-Fisher:
\(V_{\infty}^4=G\,M_{\mathrm{bar}}\,a_0\)Działa dobrze dla wielu krzywych rotacji galaktyk, ale gromady galaktyk i kosmologia pozostają trudne.
7.4 Samoistnie oddziałująca ciemna materia
Samo oddziałująca ciemna materia zakłada, że cząstki ciemnej materii oddziałują ze sobą na tyle silnie, aby zmienić profile gęstości wewnętrznego halo.
\(\frac{\sigma}{m}\sim1\text{–}100\,\mathrm{cm^2/g}\)Może to pomóc w wyjaśnieniu różnorodności rdzeni halo, ale żaden konkretny kandydat na cząstkę nie został jeszcze potwierdzony.
7.5 Pierwotne czarne dziury
Pierwotne czarne dziury powstałe we wczesnym Wszechświecie mogą stanowić część ukrytej masy. Wiele okien masy jest silnie ograniczonych przez obserwacje mikrosoczewkowania, kosmicznego mikrofalowego tła i fal grawitacyjnych.
\(10^{-16}\text{–}10^{-11}\,M_\odot\)Pozostają one spekulatywne jako pełne wyjaśnienie ukrytej masy Drogi Mlecznej.
8. Perspektywa teorii pszczół
BeeTheory proponuje, że grawitacja może być rozumiana jako efekt emergentny wynikający z zachowania falowego, a nie jako fundamentalna siła przenoszona tylko przez cząstkę lub wytwarzana tylko przez krzywiznę czasoprzestrzeni.
W tych ramach każdy masywny układ jest powiązany z funkcją falową ψ(r,t). Podstawowym kwantowym punktem wyjścia jest trójwymiarowe równanie Schrödingera:
\(i\hbar\frac{\partial\psi}{\partial t}=-\frac{\hbar^2}{2m}\nabla^2\psi+V(\mathbf{r})\psi\)Gdy dwa rozkłady masy zbliżają się do siebie, ich funkcje falowe nakładają się na siebie. Splot tych funkcji falowych można zapisać jako:
\(\Psi_{12}(\mathbf{r})=(\psi_1*\psi_2)(\mathbf{r})=\int\psi_1(\mathbf{r}’)\psi_2(\mathbf{r}-\mathbf{r}’)\,d^3r’\)BeeTheory interpretuje przyciąganie grawitacyjne jako wielkoskalową manifestację strukturalnego nakładania się fal, rezonansu i spójności pola.
8.1 Reinterpretacja ukrytej masy w teorii pszczół
W teorii pszczół to, co zwykle nazywa się ciemną materią, można interpretować jako skumulowany efekt grawitacyjny interferencji fal z wielu systemów oscylacyjnych rozmieszczonych w galaktycznym halo.
\(\rho_{\mathrm{eff}}(R)=\rho_{\mathrm{bar}}(R)+\Delta\rho_{\mathrm{wave}}(R)\)Tutaj Δρwave(R) reprezentuje dodatkową efektywną gęstość grawitacyjną wynikającą ze spójnej struktury pola falowego, a nie z bezpośrednio widocznej materii barionowej.
Termin ten musiałby odtwarzać zachowanie radialne zwykle przypisywane ciemnej materii. W szczególności musiałby on generować w przybliżeniu płaskie krzywe rotacji w odpowiednim zakresie galaktycznym.
\(\rho_{\mathrm{wave}}(R)\propto R^{-2}\)Otwarte wyzwanie ilościowe
BeeTheory musi wykazać, czy model interferencyjny oparty na falach może odtworzyć dokładny profil gęstości radialnej wymagany przez obserwowane krzywe rotacji. Musi również wyjaśnić, dlaczego efektywna masa ukryta jest często znacznie większa niż widoczna masa barionowa.
Referencje
- Ou, X., Eilers, A.-C., Necib, L., Frebel, A. – The dark matter profile of the Milky Way inferred from its circular velocity curve, MNRAS 528, 693-710, 2024.
- Navarro, J. F., Frenk, C. S., White, S. D. M. – A Universal Density Profile from Hierarchical Clustering, ApJ 490, 493, 1997.
- Einasto, J. – On the construction of a composite model for the Galaxy, Trudy 5, 87, 1965.
- Watkins, L. L., van der Marel, R. P. et al. – Evidence for an Anticorrelation between the Masses of the Milky Way and Andromeda Galaxies, ApJ 873, 111, 2019.
- Milgrom, M. – A modification of the Newtonian dynamics as a possible alternative to the hidden mass hypothesis, ApJ 270, 365, 1983.
- McGaugh, S. S. et al. – Radial Acceleration Relation in Rotationally Supported Galaxies, PRL 117, 201101, 2016.
Uwaga: jeśli strona jest wykorzystywana jako źródło cytatów naukowych, przed publikacją należy zweryfikować ostatnie lub przyszłe odniesienia.
Perspektywa BeeTheory
Problem ukrytej masy to nie tylko kwestia ilości brakującej materii. Jest to pytanie o to, jaki rodzaj struktury fizycznej wytwarza grawitację w skali galaktycznej.
Klasyczne modele ciemnej materii interpretują brakującą masę jako niewidzialną materię. BeeTheory bada uzupełniającą możliwość: część ukrytego efektu grawitacyjnego może wynikać z koherencji fal strukturalnych.
Następny krok jest matematyczny: proszę zdefiniować termin gęstości fali radialnej, wyprowadzić jego krzywą rotacji i porównać go bezpośrednio z danymi Drogi Mlecznej z okresu Gaia.