Matematiska grunder för galaktisk saknad massa: Disk, sfär, densitet, potential och radiell skalning
TL;DR: Problemet med den saknade massan uppstår när massan som härleds från galaktiska rotationskurvor överstiger den massa som direkt observeras i stjärnor, gas och stoft. Matematiskt kräver detta att man kopplar samman ytdensiteten på en skiva, volymdensiteten i tre dimensioner, gravitationspotentialen, den radiella accelerationen och den inneslutna massan.
1. Radialkoordinater och geometri
Vi skiljer mellan två geometrier:
- Diskgeometri: synlig galaktisk materia är huvudsakligen fördelad i en tunn roterande disk.
- Sfärisk geometri: mörk eller saknad massa modelleras ofta som en ungefärligt sfärisk halo.
Diskområde element:
\[ dA = R\,dR\,d\phi \]
Sfäriskt volymelement:
\[ dV = 4\pi r^2\,dr \]
Samma symbol \(r\) används ofta för galaktocentrisk radie, men betydelsen beror på geometrin. I en skiva är \(R\) en cylindrisk radie. I en halo är \(r\) vanligtvis en sfärisk radie.
2. Synlig massa på en galaktisk skiva
Den synliga skivan approximeras ofta med en exponentiell yttäthet:
\[ \Sigma(R)=\Sigma_0 e^{-R/R_d} \]
Massan av en ring mellan \(R\) och \(R+dR\) är:
\[ dM_{\rm disk}=2\pi R\Sigma(R)\,dR \]
\[ dM_{\rm disk}=2\pi R\Sigma_0e^{-R/R_d}\,dR \]
Den kumulativa synliga skivmassan är därför:
\[ 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 \vänster[ 1-e^{-R/R_d} \left( 1+\frac{R}{R_d} \right) \right] \]
Vid stor radie:
\[ M_{\rm disk}(R)\rightarrow 2\pi\Sigma_0R_d^2 \]
Den synliga diskmassan närmar sig ett ändligt värde.
3. Sfärisk mass- och volymdensitet
För en sfärisk massfördelning bestämmer volymdensiteten \(\rho(r)\) den inneslutna massan:
\[ M(r)=4\pi\int_0^r \rho(r’)r’^2\,dr’ \]
Det omvända förhållandet är:
\[ \rho(r)=\frac{1}{4\pi r^2}\frac{dM}{dr} \]
Detta förhållande är centralt för problemet med den saknade massan. Om den härledda massan växer linjärt med radien, minskar motsvarande sfäriska densitet som \(1/r^2\).
4. Dynamisk massa från cirkulär rörelse
För cirkulär rörelse uppfyller gravitationsaccelerationen:
\[ \frac{v(r)^2}{r}=\frac{GM(r)}{r^2} \]
Därför..:
\[ M_{\rm dyn}(r)=\frac{v(r)^2r}{G} \]
För en platt rotationskurva:
\[ v(r)\approx v_0 \]
\[ M_{\rm dyn}(r)\approx \frac{v_0^2}{G}r \]
Detta ger en standardskalning:
\[ M_{\rm dyn}(r)\propto r \]
\[ \rho_{\rm dyn}(r)\propto \frac{1}{r^2} \]
5. Definition av saknad massa
Den saknade massan är skillnaden mellan den dynamiska massan och den synliga massan:
\[ M_{\rm miss}(r)=M_{\rm dyn}(r)-M_{\rm vis}(r) \]
För en exponentiell synlig disk:
\[ M_{\rm miss}(r)= \frac{v(r)^2r}{G} – 2\pi\Sigma_0R_d^2 \vänster[ 1-e^{-r/R_d} \left( 1+\frac{r}{R_d} \right) \right] \]
För \(v(r)\approx v_0\):
\[ M_{\rm miss}(r)\approx \frac{v_0^2}{G}r – 2\pi\Sigma_0R_d^2 \vänster[ 1-e^{-r/R_d} \left( 1+\frac{r}{R_d} \right) \right] \]
Vid stor radie mättas disktermen, medan den dynamiska termen fortsätter att växa ungefär som \(r\).
6. Gravitationspotential i 3D
Den Newtonska gravitationspotentialen som genereras av en punktmassa är:
\[ \Phi(r)=-\frac{GM}{r} \]
Det motsvarande gravitationsfältet är den radiella derivatan av potentialen:
\[ g(r)=-\frac{d\Phi}{dr} \]
\[ g(r)=-\frac{GM}{r^2} \]
Detta förklarar förhållandet mellan \(1/r\) och \(1/r^2\): potentialen hos en lokaliserad massa sjunker som \(1/r\), medan kraften eller accelerationen sjunker som \(1/r^2\).
7. Poisson-ekvationen
Masstäthet och gravitationspotential hänger samman genom Poissons ekvation:
\[ \nabla^2\Phi=4\pi G\rho \]
I sfärisk symmetri blir detta:
\[ \frac{1}{r^2}\frac{d}{dr} \vänster( r^2\frac{d\Phi}{dr} \right) = 4\pi G\rho(r) \]
Denna ekvation kopplar samman tre storheter:
\[ \rho(r) \longrightarrow M(r) \longrightarrow \Phi(r) \longrightarrow v(r) \]
8. Potential för en utökad 3D-densitet
För en allmän 3D-densitetsfördelning är potentialen:
\[ \Phi(\mathbf{x}) = -G\int \frac{\rho(\mathbf{x}’)}{|\mathbf{x}-\mathbf{x}’|} \,d^3x’ \]
Kärnan \(1/|\mathbf{x}-\mathbf{x}’|\) är det matematiska ursprunget till \(1/r\)-potentialen i tre dimensioner.
9. Potential hos en tunn skiva
För en tunn skiva med ytdensitet \(\Sigma(R’)\) kan gravitationspotentialen i skivans plan skrivas som:
\[ \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}} \]
Avståndet mellan en fältpunkt med radien \(R\) och en källpunkt med radien \(R’\) är:
\[ d(R,R’,\phi) = \sqrt{R^2+R’^2-2RR’\cos\phi} \]
Den radiella accelerationen i skivan erhålls genom differentiering av potentialen:
\[ g_R(R)=-\frac{\partiell \Phi}{\partiell R} \]
Rotationshastigheten följer av detta:
\[ v^2(R)=R\,|g_R(R)| \]
10. Projektion av en 3D-interaktion på skivan
Om en interaktion fortplantar sig i tre dimensioner men utvärderas i diskplanet, medför den radiella projektionen en geometrisk faktor. För två punkter på skivan som är åtskilda med avståndet \(d\) är den radiella projektionsfaktorn:
\[ \cos\theta = \frac{R-R’\cos\phi}{d(R,R’,\phi)} \]
En generisk radiell 3D-kärna \(K(d)\), projicerad på skivan, ser alltså ut som:
\[ K_{\rm disk}(R,R’,\phi) = K(d) \frac{R-R’\cos\phi}{d} \]
Till exempel har en Newtonsk kraftliknande kärna:
\[ K(d)\propto \frac{1}{d^2} \]
\[ K_{\rm disk}\propto \frac{R-R’\cos\phi}{d^3} \]
En exponentiell 3D-kärna kan skrivas som:
\[ K(d)\propto e^{-d/\lambda} \]
\[ K_{\rm disk}\propto e^{-d/\lambda} \frac{R-R’\cos\phi}{d} \]
11. Exponentiella kärnor i tre dimensioner
En rent exponentiell radiell faktor har formen:
\[ e^{-r/\lambda} \]
I tredimensionell fältteori uppträder ofta en exponentiellt avskärmad potential i Yukawa-liknande form:
\[ \Phi_Y(r)\propto -\frac{e^{-r/\lambda}}{r} \]
Den tillhörande radialkraften innehåller både \(1/r^2\) och exponentiella termer:
\[ g_Y(r) = -\frac{d}{dr} \vänster( \frac{e^{-r/\lambda}}{r} \höger) \]
\[ g_Y(r)\propto e^{-r/\lambda} \vänster( \frac{1}{r^2} + \frac{1}{\lambda r} \right) \]
Detta visar varför exponentiellt radiellt beteende i 3D inte är oberoende av \(1/r\)-geometrin. Exponentialen styr dämpningen, medan \(1/r\) och \(1/r^2\) uppstår genom tredimensionell spridning.
12. Lagar för radiell skalning
Problemet med den saknade massan är starkt kopplat till radiell skalning. Flera viktiga radiella lagar uppträder upprepade gånger:
| Kvantitet | Typisk skalning | Betydelse |
|---|---|---|
| Potential för punktmassa | \(\Phi(r)\sim 1/r\) | 3D grön funktion av gravitation |
| Kraft från punktmassa | \(g(r)\sim 1/r^2\) | Derivatan av \(1/r\) |
| Platt rotationshastighet | \(v(r)\sim konstant\) | Observerad i yttre galaktiska skivor |
| Dynamisk massa | \(M(r)\sim r\) | Krävs av platt rotation |
| Halodensitet | \(\rho(r)\sim 1/r^2\) | Ger \(M(r)\sim r\) |
| Exponentiell skiva | \(\Sigma(R)\sim e^{-R/R_d}\) | Synlig skiva bleknar snabbt |
| Skärmad 3D-potential | \(\Phi(r)\sim e^{-r/\lambda}/r\) | Exponentiell dämpning plus 3D-spridning |
13. Från densitet till rotationskurva
För en sfärisk halo med densitet:
\[ \rho(r)\propto \frac{1}{r^2} \]
den inneslutna massan är:
\[ M(r)=4\pi\int_0^r \rho(r’)r’^2dr’ \]
\[ M(r)\propto r \]
Då så:
\[ v^2(r)=\frac{GM(r)}{r} \]
\[ v^2(r)\approx konstant \]
\[ v(r)\approx konstant \]
Detta är den matematiska bryggan mellan en halotäthet på \(1/r^2\) och en platt galaktisk rotationskurva.
14. Från diskmassa till saknad massa
Den synliga skivans massa växer snabbt till en början och mättas sedan:
\[ M_{\rm disk}(R)\rightarrow M_d \]
Den dynamiska massan som härleds från en platt rotationskurva fortsätter att växa:
\[ M_{\rm dyn}(r)\propto r \]
Därför beter sig den saknade massan ungefär som:
\[ 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 \]
Vid tillräckligt stor radie:
\[ M_{\rm miss}(r)\sim r \]
\[ \rho_{\rm miss}(r)\sim \frac{1}{r^2} \]
15. Matematisk varning: disk och sfär är inte utbytbara
Ekvationen
\[ M(r)=\frac{v(r)^2r}{G} \]
är exakt för sfärisk symmetri. För en tillplattad skiva bör man beräkna potentialen genom att integrera över skivan och sedan härleda den radiella accelerationen. Det sfäriska uttrycket används ofta som en effektiv approximation, särskilt när man diskuterar den massa som behövs för att stödja en given rotationskurva.
16. Sammanfattning av nyckelekvationer
Diskens ytdensitet:
\[ \Sigma(R)=\Sigma_0e^{-R/R_d} \]
Diskens masselement:
\[ dM_{\rm disk}=2\pi R\Sigma(R)dR \]
Synlig diskmassa:
\[ M_{\rm disk}(R)=2\pi\Sigma_0R_d^2 \Vänster[ 1-e^{-R/R_d} \left( 1+\frac{R}{R_d} \right) \right] \]
Sfärisk volymmassa:
\[ M(r)=4\pi\int_0^r\rho(r’)r’^2dr’ \]
Densitet från innesluten massa:
\[ \rho(r)=\frac{1}{4\pi r^2}\frac{dM}{dr} \]
Dynamisk massa:
\[ M_{\rm dyn}(r)=\frac{v(r)^2r}{G} \]
Saknad massa:
\[ M_{\rm miss}(r)=M_{\rm dyn}(r)-M_{\rm vis}(r) \]
Newtonsk potential:
\[ \Phi(r)=-\frac{GM}{r} \]
Poisson-ekvationen:
\[ \nabla^2\Phi=4\pi G\rho \]
3D potentialintegral:
\[ \Phi(\mathbf{x}) = -G\int \frac{\rho(\mathbf{x}’)}{|\mathbf{x}-\mathbf{x}’|} d^3x’ \]
Potential för tunn skiva:
\[ \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}} \]
Projicerad radiell kärna:
\[ \cos\theta= \frac{R-R’\cos\phi} {\sqrt{R^2+R’^2-2RR’\cos\phi}} \]
Exponentiell 3D-skärmad potential:
\[ \Phi_Y(r)\propto -\frac{e^{-r/\lambda}}{r} \]
Slutsats
Problemet med den saknade massan är en matematisk missmatchning mellan två radiella beteenden. Den synliga skivan följer en exponentiell ytdensitet och når en ändlig kumulativ massa. Den dynamiska massan som härleds från ungefär flacka rotationskurvor växer ungefär linjärt med radien. Om den tolkas som en sfärisk halo motsvarar detta en densitet som minskar ungefär som \(1/r^2\). Ekvationerna för skivintegration, sfäriska skal, gravitationspotential och radiell projektion ger det matematiska språk som behövs för att analysera denna missmatchning.