Matematyczne podstawy galaktycznej brakującej masy: Dysk, kula, gęstość, potencjał i skalowanie radialne
TL;DR: Problem brakującej masy pojawia się, gdy masa wywnioskowana z krzywych rotacji galaktyk przekracza masę bezpośrednio obserwowaną w gwiazdach, gazie i pyle. Matematycznie wymaga to połączenia gęstości powierzchniowej na dysku, gęstości objętościowej w trzech wymiarach, potencjału grawitacyjnego, przyspieszenia radialnego i masy zamkniętej.
1. Współrzędne radialne i geometria
Wyróżniamy dwie geometrie:
- Geometria dysku: widoczna materia galaktyczna jest głównie rozmieszczona w cienkim, obracającym się dysku.
- Geometria sferyczna: ciemna lub brakująca masa jest często modelowana jako z grubsza sferyczna aureola.
Element obszaru dysku:
\[ dA = R\,dR\,d\phi \]
Element objętości sferycznej:
\[ dV = 4\pi r^2\,dr \]
Ten sam symbol \(r\) jest często używany dla promienia galaktocentrycznego, ale jego znaczenie zależy od geometrii. W dysku \(R\) jest promieniem cylindrycznym. W halo \(r\) jest zazwyczaj promieniem sferycznym.
2. Widoczna masa na dysku galaktycznym
Widoczny dysk jest często przybliżany przez wykładniczą gęstość powierzchniową:
\[ \Sigma(R)=\Sigma_0 e^{-R/R_d} \]
Masa pierścienia pomiędzy \(R\) i \(R+dR\) wynosi:
\[ dM_{\rm disk}=2\pi R\Sigma(R)\,dR \]
\[ dM_{\rm disk}=2\pi R\Sigma_0e^{-R/R_d}\,dR \]
Łączna masa widocznego dysku wynosi zatem:
\[ M_{\rm disk}(R)=2\pi\int_0^R R’\Sigma_0e^{-R’/R_d}\,dR’ \]
\[ M_{\rm disk}(R)=2\pi\Sigma_0R_d^2 \left[ 1-e^{-R/R_d} \left( 1+\frac{R}{R_d} \right) \right] \]
Przy dużym promieniu:
\[ M_{\rm disk}(R)\rightarrow 2\pi\Sigma_0R_d^2 \]
Widoczna masa dysku zbliża się do skończonej wartości.
3. Masa sferyczna i gęstość objętościowa
W przypadku sferycznego rozkładu masy gęstość objętościowa \(\rho(r)\) określa masę zamkniętą:
\[ M(r)=4\pi\int_0^r \rho(r’)r’^2\,dr’ \]
Zależność odwrotna jest następująca:
\[ \rho(r)=\frac{1}{4\pi r^2}\frac{dM}{dr} \]
Zależność ta jest kluczowa dla problemu brakującej masy. Jeśli wywnioskowana masa rośnie liniowo wraz z promieniem, to odpowiadająca jej gęstość sferyczna maleje jako \(1/r^2\).
4. Masa dynamiczna z ruchu kołowego
Dla ruchu kołowego przyspieszenie grawitacyjne spełnia warunek:
\[ \frac{v(r)^2}{r}=\frac{GM(r)}{r^2} \]
W związku z tym:
\[ M_{\rm dyn}(r)=\frac{v(r)^2r}{G} \]
Dla płaskiej krzywej rotacji:
\[ v(r)\approx v_0 \]
\[ M_{\rm dyn}(r)\approx \frac{v_0^2}{G}r \]
Daje to standardowe skalowanie:
\[ M_{\rm dyn}(r)\propto r \]
\[ \rho_{\rm dyn}(r)\propto \frac{1}{r^2} \]
5. Definicja brakującej masy
Brakująca masa to różnica między masą dynamiczną a masą widzialną:
\[ M_{\rm miss}(r)=M_{\rm dyn}(r)-M_{\rm vis}(r) \]
Dla widocznego dysku wykładniczego:
\[ M_{\rm miss}(r)= \frac{v(r)^2r}{G} – 2\pi\Sigma_0R_d^2 \left[ 1-e^{-r/R_d} \left( 1+\frac{r}{R_d} \right) \right] \]
Dla \(v(r)\approx v_0\):
\[ M_{\rm miss}(r)\approx \frac{v_0^2}{G}r – 2\pi\Sigma_0R_d^2 \left[ 1-e^{-r/R_d} \left( 1+\frac{r}{R_d} \right) \right] \]
Przy dużym promieniu człon dyskowy nasyca się, podczas gdy człon dynamiczny nadal rośnie w przybliżeniu jako \(r\).
6. Potencjał grawitacyjny w 3D
Newtonowski potencjał grawitacyjny generowany przez masę punktową wynosi:
\[ \Phi(r)=-\frac{GM}{r} \]
Odpowiadające mu pole grawitacyjne jest pochodną radialną potencjału:
\[ g(r)=-\frac{d\Phi}{dr} \]
\[ g(r)=-\frac{GM}{r^2} \]
Wyjaśnia to związek między \(1/r\) i \(1/r^2\): potencjał zlokalizowanej masy spada jako \(1/r\), podczas gdy siła lub przyspieszenie spada jako \(1/r^2\).
7. Równanie Poissona
Gęstość masy i potencjał grawitacyjny są powiązane równaniem Poissona:
\[ \nabla^2\Phi=4\pi G\rho \]
W symetrii sferycznej jest to:
\[ \frac{1}{r^2}\frac{d}{dr} \left( r^2\frac{d\Phi}{dr} \right) = 4\pi G\rho(r) \]
Równanie to łączy trzy wielkości:
\[ \rho(r) \longrightarrow M(r) \longrightarrow \Phi(r) \podłużny strzałka v(r) \]
8. Potencjał rozszerzonej gęstości 3D
Dla ogólnego rozkładu gęstości 3D potencjał wynosi:
\[ \Phi(\mathbf{x}) = -G\int \frac{\rho(\mathbf{x}’)}{|\mathbf{x}-\mathbf{x}’|} \,d^3x’ \]
Jądro \(1/|\mathbf{x}-\mathbf{x}’|\) jest matematycznym źródłem potencjału \(1/r\) w trzech wymiarach.
9. Potencjał cienkiego dysku
Dla cienkiego dysku o gęstości powierzchniowej \(\Sigma(R’)\) potencjał grawitacyjny w płaszczyźnie dysku można zapisać jako:
\[ \Phi(R) = -G \int_0^\infty \int_0^{2\pi} \frac{\Sigma(R’)R’\,dR’\,d\phi} {\sqrt{R^2+R’^2-2RR’\cos\phi}} \]
Odległość między punktem pola o promieniu \(R\) a punktem źródła o promieniu \(R’\) wynosi:
\[ d(R,R’,\phi) = \sqrt{R^2+R’^2-2RR’\cos\phi} \]
Przyspieszenie radialne w dysku uzyskuje się poprzez różniczkowanie potencjału:
\[ g_R(R)=-\frac{\partial \Phi}{\partial R} \]
Prędkość obrotowa wynika z:
\[ v^2(R)=R\,|g_R(R)| \]
10. Projekcja interakcji 3D na dysk
Jeśli oddziaływanie propaguje się w trzech wymiarach, ale jest oceniane w płaszczyźnie dysku, projekcja radialna wprowadza współczynnik geometryczny. Dla dwóch punktów na dysku oddzielonych odległością \(d\) współczynnik projekcji radialnej wynosi:
\[ \cos\theta = \frac{R-R’\cos\phi}{d(R,R’,\phi)} \]
Zatem ogólne jądro radialne 3D \(K(d)\), rzutowane na dysk, wygląda następująco:
\[ K_{\rm disk}(R,R’,\phi) = K(d) \frac{R-R’\cos\phi}{d} \]
Na przykład, jądro podobne do siły Newtona ma:
\[ K(d)\propto \frac{1}{d^2} \]
\[ K_{\rm disk}\propto \frac{R-R’\cos\phi}{d^3} \]
Wykładnicze jądro 3D można zapisać jako:
\[ K(d)\propto e^{-d/\lambda} \]
\[ K_{\rm disk}\propto e^{-d/\lambda} \frac{R-R’\cos\phi}{d} \]
11. Jądra wykładnicze w trzech wymiarach
Czysto wykładniczy współczynnik radialny ma postać:
\[ e^{-r/\lambda} \]
W trójwymiarowej teorii pola, wykładniczo ekranowany potencjał często pojawia się w formie podobnej do Yukawy:
\[ \Phi_Y(r)\propto -\frac{e^{-r/\lambda}}{r} \]
Związana z tym siła promieniowa zawiera zarówno \(1/r^2\), jak i wyrazy wykładnicze:
\[ g_Y(r) = -\frac{d}{dr} \left( \frac{e^{-r/\lambda}}{r} \right) \]
\[ g_Y(r)\propto e^{-r/\lambda} \left( \frac{1}{r^2} + \frac{1}{\lambda r} \right) \]
To pokazuje, dlaczego wykładnicze zachowanie radialne w 3D nie jest niezależne od geometrii \(1/r\). Wykładniczy kontroluje tłumienie, podczas gdy \(1/r\) i \(1/r^2\) wynikają z trójwymiarowego rozprzestrzeniania się.
12. Promieniowe prawa skalowania
Problem brakującej masy jest silnie związany ze skalowaniem radialnym. Kilka ważnych praw radialnych pojawia się wielokrotnie:
| Ilość | Typowe skalowanie | Znaczenie |
|---|---|---|
| Potencjał masy punktowej | \(\Phi(r)\sim 1/r\) | 3D Zielona funkcja grawitacji |
| Siła masy punktowej | \(g(r)\sim 1/r^2\) | Pochodna \(1/r\) |
| Prędkość obrotu płaskiego | \(v(r)\sim stała\) | Obserwowana w zewnętrznych dyskach galaktycznych |
| Masa dynamiczna | \(M(r)\sim r\) | Wymagana przez płaską rotację |
| Gęstość halo | \(\rho(r)\sim 1/r^2\) | Daje \(M(r)\sim r\) |
| Dysk wykładniczy | \(\Sigma(R)\sim e^{-R/R_d}\) | Widoczny dysk szybko zanika |
| Ekranowany potencjał 3D | \(\Phi(r)\sim e^{-r/\lambda}/r\) | Tłumienie wykładnicze plus rozprzestrzenianie 3D |
13. Od gęstości do krzywej rotacji
Dla sferycznej aureoli o gęstości:
\[ \rho(r)\propto \frac{1}{r^2} \]
załączona masa wynosi:
\[ M(r)=4\pi\int_0^r \rho(r’)r’^2dr’ \]
\[ M(r)\propto r \]
Następnie:
\[ v^2(r)=\frac{GM(r)}{r} \]
\[ v^2(r)\approx stała \]
\[ v(r)\approx stała \]
Jest to matematyczny pomost między gęstością halo \(1/r^2\) a płaską krzywą rotacji galaktyki.
14. Od masy dysku do brakującej masy
Widoczna masa dysku początkowo szybko rośnie, a następnie nasyca się:
\[ M_{\rm disk}(R)\rightarrow M_d \]
Masa dynamiczna wywnioskowana z płaskiej krzywej rotacji stale rośnie:
\[ M_{\rm dyn}(r)\propto r \]
Dlatego brakująca masa zachowuje się w przybliżeniu tak:
\[ M_{\rm miss}(r)=M_{\rm dyn}(r)-M_{\rm disk}(r) \]
\[ M_{\rm miss}(r)\approx \frac{v_0^2}{G}r-M_d \]
Przy wystarczająco dużym promieniu:
\[ M_{\rm miss}(r)\sim r \]
\[ \rho_{\rm miss}(r)\sim \frac{1}{r^2} \]
15. Ostrzeżenie matematyczne: dysk i kula nie są zamienne
Równanie
\[ M(r)=\frac{v(r)^2r}{G} \]
jest dokładna dla symetrii sferycznej. W przypadku spłaszczonego dysku należy obliczyć potencjał poprzez całkowanie po dysku, a następnie wyprowadzić przyspieszenie radialne. Wyrażenie sferyczne jest często używane jako efektywne przybliżenie, szczególnie podczas omawiania masy potrzebnej do podtrzymania danej krzywej rotacji.
16. Podsumowanie kluczowych równań
Gęstość powierzchni dysku:
\[ \Sigma(R)=\Sigma_0e^{-R/R_d} \]
Element masy dysku:
\[ dM_{\rm disk}=2\pi R\Sigma(R)dR \]
Widoczna masa dysku:
\[ M_{\rm disk}(R)=2\pi\Sigma_0R_d^2 \left[ 1-e^{-R/R_d} \left( 1+\frac{R}{R_d} \right) \right] \]
Sferyczna masa objętościowa:
\[ M(r)=4\pi\int_0^r\rho(r’)r’^2dr’ \]
Gęstość od masy zamkniętej:
\[ \rho(r)=\frac{1}{4\pi r^2}\frac{dM}{dr} \]
Masa dynamiczna:
\[ M_{\rm dyn}(r)=\frac{v(r)^2r}{G} \]
Brakująca masa:
\[ M_{\rm miss}(r)=M_{\rm dyn}(r)-M_{\rm vis}(r) \]
Potencjał newtonowski:
\[ \Phi(r)=-\frac{GM}{r} \]
Równanie Poissona:
\[ \nabla^2\Phi=4\pi G\rho \]
Całka potencjału 3D:
\[ \Phi(\mathbf{x}) = -G\int \frac{\rho(\mathbf{x}’)}{|\mathbf{x}-\mathbf{x}’|} d^3x’ \]
Potencjał cienkiego dysku:
\[ \Phi(R) = -G \int_0^\infty \int_0^{2\pi} \frac{\Sigma(R’)R’dR’d\phi} {\sqrt{R^2+R’^2-2RR’\cos\phi}} \]
Prognozowane jądro radialne:
\[ \cos\theta= \frac{R-R’\cos\phi} {\sqrt{R^2+R’^2-2RR’\cos\phi}} \]
Wykładniczy potencjał ekranowany 3D:
\[ \Phi_Y(r)\propto -\frac{e^{-r/\lambda}}{r} \]
Wnioski
Problem brakującej masy jest matematycznym niedopasowaniem pomiędzy dwoma zachowaniami radialnymi. Widoczny dysk podąża za wykładniczą gęstością powierzchniową i osiąga skończoną masę skumulowaną. Masa dynamiczna wywnioskowana z w przybliżeniu płaskich krzywych rotacji rośnie w przybliżeniu liniowo wraz z promieniem. Jeśli zinterpretować ją jako sferyczne halo, odpowiada to gęstości malejącej w przybliżeniu jako \(1/r^2\). Równania całkowania dysku, powłoki sferycznej, potencjału grawitacyjnego i projekcji radialnej dostarczają języka matematycznego potrzebnego do analizy tego niedopasowania.